| 查看: 1292 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
为什么拟合结果总为所设初始值？
|
||
|
用这样一个函数：yy=1./(1+5.*b./(2.*k)+2.*b.^2./k.^2+b.^3./(2.*k.^3)).*b.*exp(b.*xx).*(xx<0)+(b./k+2.*b.^2./k.^2+b.^3./k.^3)./(4+10.*b./k+8.*b.^2./k.^2+2.*b.^3./k.^3).*k.*(k.*xx+2./(b./k+1)).^2.*exp(-k.*xx).*(xx>=0)去拟合（xx,yy） myfittype=fittype('1./(1+5.*b./(2.*k)+2.*b.^2./k.^2+b.^3./(2.*k.^3)).*b.*exp(b.*xx).*(xx<0)+(b./k+2.*b.^2./k.^2+b.^3./k.^3)./(4+10.*b./k+8.*b.^2./k.^2+2.*b.^3./k.^3).*k.*(k.*xx+2./(b./k+1)).^2.*exp(-k.*xx).*(xx>=0)','dependent',{'yy'},'independent',{'xx'},'coefficients',{'b','k'}); myfit=fit(xx',yy',myfittype,'start',[4.2e9,1.4e9]) plot(myfit,xx,yy) myfit = General model: myfit(xx) = 1./(1+5.*b./(2.*k)+2.*b.^2./k.^2+b.^3./(2.*k.^3)).*b.*exp(b.*xx) .*(xx<0)+(b./k+2.*b.^2./k.^2+b.^3./k.^3)./(4+10.*b./k+ 8.*b.^2./k.^2+2.*b.^3./k.^3).*k.*(k.*xx+2./(b./k+1) ).^2.*exp(-k.*xx).*(xx>=0) Coefficients (with 95% confidence bounds): b = 4.2e+09 (3.72e+09, 4.68e+09) k = 1.4e+09 (1.369e+09, 1.431e+09) 拟合时设定了初始值，但无论初始值设为多少，拟合结果就是多少，请问问题出在那里了呢？ |
回复此楼» 猜你喜欢
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
两块石头
已经有8人回复
-
初秋的晨风
已经有6人回复
-
学科评审组评审是指会评吗？
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有7人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有17人回复
-
Ei源刊怎么投
已经有4人回复
-
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生
已经有3人回复
4楼2016-11-18 10:11:38
2楼2016-10-28 21:44:33
dlli2014
新虫 (文坛精英)
- 应助: 0 (幼儿园)
- 金币: 33610.4
- 红花: 4
- 帖子: 24263
- 在线: 286.6小时
- 虫号: 3542502
- 注册: 2014-11-17
- 专业: 凝聚态物性 II ：电子结构
3楼2016-10-28 23:38:53