24小时热门版块排行榜    

查看: 1406  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

beina

新虫 (正式写手)

[求助] 关于matlab的一个基础问题

请懂matlab的高手指点,不胜感激!
我新建了一个M函数
function F=myfun(x,xdata)
F=0.6518 + .7713./atan(xdata) - .0911.*atan(xdata) + ((.9121.* xdata + 1.572)/(.9121.* xdata + 1.59 + 2.3026./ xdata + 0.0637*log(1 + 15.7./ xdata)) - 1)/((1 + 5.04*x.* xdata)*(1 + (x^0.674).* (xdata.^0.173)));
end
这个函数书写上哪里有错误,是不是在点乘点除上的错误?
我运行了下面的命令后,会有如下错误
>> load data.txt
>> xdata=data(:,1);
>> ydata=data(:,2);
>> x=lsqcurvefit(@myfun,0.05,xdata,ydata)
??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==> myfun at 4
F=0.6518 + .7713./atan(xdata) - .0911.*atan(xdata) + ((.9121.* xdata + 1.572)/(.9121.* xdata + 1.59 + 2.3026./ xdata + 0.0637*log(1 + 15.7./ xdata)) - 1)/((1 + 5.04*x.* xdata)*(1
+ (x^0.674).* (xdata
Error in ==> lsqncommon at 152
            fuser = funfcn{3}(xargin{:},varargin{:});

Error in ==> lsqcurvefit at 186
[x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...

Caused by:
    Failure in initial user-supplied objective function evaluation. LSQCURVEFIT cannot continue.
难道lsqcurvefit的命令中书写也有错误?请大侠指点,感激涕零!
回复此楼
流泪撒种的,必欢呼收割
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beina

新虫 (正式写手)

送鲜花一朵
引用回帖:
: Originally posted by renchong at 2011-12-24 17:10:54:
想了一下 LZ的xdata和ydata在程序中没见到有具体的值啊,还有LZ的点乘有问题,请参考一下这个吧(我是在添加了xdata和ydata值后做的)能算出结果 F=0.6518 + 0.7713./atan(xdata) - 0.0911*atan(xdata) + ((0.912 ...

按照你的写法果然对了,但是又有错误了,呜呜
??? Error using ==> qr
Complex sparse QR is not yet available.

Error in ==> aprecon at 57
   RPCMTX = qr(TM(:,p));

Error in ==> trdog at 47
         [R,permR] = feval(pcmtx,H,pcoptions,DM,DG,varargin{:});

Error in ==> snls at 346
      [sx,snod,qp,posdef,pcgit,Z] = trdog(x,g,A,D,delta,dv,...

Error in ==> lsqncommon at 285
    [x,FVAL,LAMBDA,JACOB,EXITFLAG,OUTPUT,msg]=...

Error in ==> lsqcurvefit at 186
[x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...
xdata ydata有个txt文件的值,没有写上去
流泪撒种的,必欢呼收割
5楼2011-12-25 09:51:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

renchong

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
想了一下 LZ的xdata和ydata在程序中没见到有具体的值啊,还有LZ的点乘有问题,请参考一下这个吧(我是在添加了xdata和ydata值后做的)能算出结果 F=0.6518 + 0.7713./atan(xdata) - 0.0911*atan(xdata) + ((0.9121* xdata + 1.572)./(0.9121* xdata + 1.59 + 2.3026./ xdata + 0.0637*log(1 + 15.7./ xdata)) - 1)./((1 + 5.04*x* xdata).*(1 + (x^0.674)* (xdata.^0.173)));

» 本帖已获得的红花(最新10朵)

3楼2011-12-24 17:10:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

persistent06

金虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
就是上面的除法错了,应该用点除
4楼2011-12-24 17:43:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工学调剂 +11 用户573181 2026-03-20 16/800 2026-03-24 20:36 by :@sbs123
[考研] 0854AI CV方向招收调剂 +3 章小鱼567 2026-03-23 3/150 2026-03-24 20:25 by 汪!?!
[考研] 招08考数学 +7 laoshidan 2026-03-20 16/800 2026-03-24 17:52 by 乌拉儿山脉
[考研] 081700 调剂 267分 +9 迷人的哈哈 2026-03-23 9/450 2026-03-24 11:58 by 544594351
[考研] 一志愿华东理工大学081700,初试分数271 +5 kotoko_ik 2026-03-23 6/300 2026-03-24 10:29 by 学术搬砖er
[考研] 材料专业求调剂 +11 hanamiko 2026-03-18 11/550 2026-03-23 23:12 by peike
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 303求调剂 +4 元夕元 2026-03-20 4/200 2026-03-23 19:00 by macy2011
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生(材料、化学、电化学,环境) +3 我爱学电池 2026-03-23 3/150 2026-03-23 17:16 by AZMK
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
信息提示
请填处理意见