| 查看: 1155 | 回复: 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) 拟合时设定了初始值，但无论初始值设为多少，拟合结果就是多少，请问问题出在那里了呢？ |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
求推荐博导
已经有4人回复
-
为什么nbs上溴 没有产物点出现呢
已经有5人回复
-
存款400万可以在学校里躺平吗
已经有34人回复
-
求助一下有机合成大神
已经有4人回复
-
最失望的一年
已经有16人回复
-
求推荐英文EI期刊
已经有5人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
疑惑？
已经有5人回复
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
4楼2016-11-18 10:11:38