| 查看: 1402 | 回复: 4 | ||
[求助]
关于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的命令中书写也有错误？请大侠指点，感激涕零！ |
回复此楼» 猜你喜欢
机械学硕调剂互助
已经有6人回复
-
求调剂
已经有0人回复
-
机械工程论文润色/翻译怎么收费?
已经有293人回复
-
面上参与人员重要吗？找不到参与人
已经有5人回复
-
一志愿北京科技大学 085501 348分 四级已过 求调剂
已经有0人回复
-
一志愿北京科技大学 085501 348分 四级已过 求调剂
已经有0人回复
-
求调剂
已经有3人回复
-
085502车辆工程281分接受夸专业，求b区08专业收留甘愿当牛做马
已经有0人回复
-
机械专硕376求调剂，本科湖南大学
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于Matlab作图的问题
已经有14人回复
- 关于matlab的一个诡异问题
已经有7人回复
- 一个关于matlab的遗传算法的问题
已经有3人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 关于matlab 2010b 的一个小问题
已经有7人回复
- 关于MATLAB和OpenCV对相同算法的运算速度对比问题
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】问一个关于matlab中legend的问题
已经有11人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
 |
2楼2011-12-24 10:54:23
renchong
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 743.6
- 散金: 530
- 帖子: 165
- 在线: 48.1小时
- 虫号: 1522738
- 注册: 2011-12-04
- 性别: GG
- 专业: 传动机械学
【答案】应助回帖
感谢参与，应助指数 +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
persistent06
金虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 2658.7
- 散金: 136
- 红花: 3
- 帖子: 1144
- 在线: 342.5小时
- 虫号: 1212755
- 注册: 2011-02-25
- 性别: GG
- 专业: 应用有机化学
4楼2011-12-24 17:43:29
送鲜花一朵 |
按照你的写法果然对了，但是又有错误了，呜呜 ??? 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
5