| 查看: 1181 | 回复: 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) 拟合时设定了初始值，但无论初始值设为多少，拟合结果就是多少，请问问题出在那里了呢？ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
-
为什么中国大学教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
-
“人文社科而论，许多学术研究还没有达到民国时期的水平”
已经有5人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有11人回复
-
什么是人一生最重要的？
已经有4人回复
-
版面费该交吗
已经有3人回复
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