| 查看: 1400 | 回复: 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的命令中书写也有错误？请大侠指点，感激涕零！ |
回复此楼» 猜你喜欢
中国科学院-材料计算小组寻求2026级研究生（包括调剂生/联合培养等）
已经有132人回复
-
378分一志愿大连理工大学车辆工程调剂
已经有0人回复
-
金属材料论文润色/翻译怎么收费?
已经有263人回复
-
08工科求调剂
已经有4人回复
-
天津高校航空材料先进制造技术与装备课题组拟接收调剂硕士研究生5-6名。
已经有0人回复
-
294 英二数二物化 求调剂
已经有6人回复
-
机械352分求调剂，本科太原理工大学
已经有1人回复
-
英国格拉斯哥大学招生机器人方向博士生
已经有0人回复
-
找金属材料类期刊
已经有2人回复
-
283求调剂
已经有0人回复
-
西华大学材料学院表面科学与工程技术科研团队2026年招收研究生
已经有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
