数学建模社区-数学中国
标题:
matlab拟合问题~不知道错在哪里~
[打印本页]
作者:
guizhanbu
时间:
2009-3-22 23:14
标题:
matlab拟合问题~不知道错在哪里~
fun=inline('x(1)./(1+(x(1)./3315.84-1).*exp(-x(2).*tdata))','x','tdata');
u% c% o; U' |' z
tdata=1970:1:2008
0 P2 A1 V1 P2 |0 C; {
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]]
. E. h7 I0 C) r1 L- E1 c
x=lsqcurvefit(fun,[5000,0.2],tdata,cdata)
$ R7 U5 V* k- A8 Q( p
f=fun(x,tdata)
# q* F _3 i5 R& Y1 r- `# ^* l, P/ }
norm(feval(fun,cdata,tdata)-cdata)^2
/ z. X8 ^. C/ r0 r! p
plot(tdata,cdata)
) {5 n S1 e1 o6 Z# O0 y$ T
) A0 g6 O* x: @
为什么怎么改变,[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