| 查看: 1431 | 回复: 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的命令中书写也有错误？请大侠指点，感激涕零！ |
回复此楼» 猜你喜欢
重庆长江师范学院材料工程接收材料类和机械类考生调剂硕士研究生
已经有26人回复
-
重庆长江师范学院材料工程接收材料类考生调剂硕士研究生
已经有0人回复
-
机械工程论文润色/翻译怎么收费?
已经有79人回复
-
11408 369分调剂求助
已经有2人回复
-
天津大学化学工程专业专业课
已经有0人回复
-
二氧化铂的XRD和XPS 检测结果完全对不上，这是为什么？？？？
已经有0人回复
-
厦大公卫招收人工智能等相关专业调剂
已经有0人回复
-
第一次投Proc IMechE Part C被拒，推荐一下其他期刊
已经有0人回复
-
332，085601求调剂
已经有14人回复
-
英语一来了就有
已经有0人回复
-
福建理工大学招收英语一 数学一或者二
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于Matlab作图的问题
已经有14人回复
- 关于matlab的一个诡异问题
已经有7人回复
- 一个关于matlab的遗传算法的问题
已经有3人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 关于matlab 2010b 的一个小问题
已经有7人回复
- 关于MATLAB和OpenCV对相同算法的运算速度对比问题
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】问一个关于matlab中legend的问题
已经有11人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
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))); |
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
