数学建模社区-数学中国
标题:
matlab拟合问题~不知道错在哪里~
[打印本页]
作者:
guizhanbu
时间:
2009-3-22 23:14
标题:
matlab拟合问题~不知道错在哪里~
fun=inline('x(1)./(1+(x(1)./3315.84-1).*exp(-x(2).*tdata))','x','tdata');
$ Z2 k2 t( j$ E5 `1 \0 `6 x
tdata=1970:1:2008
% Q1 P: v/ p ]8 j8 A& |4 L
cdata=[[3315.84,3389.89,3450.79,3512.92,3560.87,3614.47,3662.82,3707.10,3750.96,3792.33,3826.58,3871.51,3924.32,3963.10,3993.09,4029.56,4070.07,4121.19,4169.85,4208.88,4234.91,4269.5,4304.4,4334.8,4363.7,4389,4413,4434.8,4456.2,4475.4,4679.91,4697.27,4730.76,4763.46,4803.48,4898,4980,5060,5120]]
+ [1 J: h) H8 F- l- y
x=lsqcurvefit(fun,[5000,0.2],tdata,cdata)
, G- x1 F" V5 Q3 v
f=fun(x,tdata)
$ f7 g9 p4 j8 V7 P1 c
norm(feval(fun,cdata,tdata)-cdata)^2
) O: j0 G: e' B" I3 ? }. }
plot(tdata,cdata)
2 J: h) @+ _6 c
3 b- i) O" |" P( i$ `
为什么怎么改变,[5000,0.2],这个初值,但结果还是不变呢,有没有高手解小生心中的疑惑啊~
作者:
mcm-dlu-edu
时间:
2009-3-23 13:13
用cftool可做拟合!!!很方便!!!!!!
作者:
hecuan1123
时间:
2009-3-23 23:55
顶一个再顶一个
作者:
guizhanbu
时间:
2009-3-24 12:51
cftool,有没有这样的课件啊~我不知道哪里有~或者发我的邮箱,谢谢这位大哥了啊~
guizhanbu@126.com
作者:
guizhanbu
时间:
2009-3-24 22:52
已经做出来了,用cftool做是很简单,谢谢这个大哥了~
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5