| 查看: 1222 | 回复: 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一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
论文撤稿了
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
河北省自然科学基金
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有8人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有5人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有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
10