| 查看: 2954 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
需要用matlab拟合数据的外行人急需帮助
|
|||
|
我有一函数y=f(x),里面有个未知数k，还有个x,y数值的txt文件，用这些条件求k 首先建立了一个M文件， function F=myfun(x,xdata) 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))); 然后就运行下面的命令： >> load data.txt >> xdata=data(:,1); >> ydata=data(:,2); >> x=lsqcurvefit(@myfun,0.05,xdata,ydata) 然后就有错误了 ??? 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] = ...  这是哪里的错误呢， lsqcurvefit 的运用中有木有语法错误，哪位大侠可以帮我解答一下，感激涕零，如果能帮我跑一下程序就最好不过了…… |
回复此楼» 猜你喜欢
时间戳又变了
已经有8人回复
-
言之凿凿，8月21号(明天)放榜
已经有16人回复
-
时间戳变了，能看出什么问题？
已经有21人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
只有每年这种时候来逛逛小木虫
已经有21人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 非线性拟合估计方程参数（matlab，spss，origin都可以）
已经有4人回复
- 请教matlab拟合
已经有9人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
3楼2011-12-26 14:36:57
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 程序强帖: 2
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
2楼2011-12-26 10:48:56
4楼2011-12-26 17:33:50
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2011-12-26 17:35:31
