数学建模社区-数学中国
标题:
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; _; e$ a( z
tdata=1970:1:2008
- L9 n# ~# s' P. Z/ A( h
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]]
9 G7 `5 Q* g0 N0 O4 X0 u7 p7 R) {! f
x=lsqcurvefit(fun,[5000,0.2],tdata,cdata)
( z8 `1 Q* ], N1 _% R1 g9 G- |% n
f=fun(x,tdata)
& d7 S- {3 l! Y1 b3 D
norm(feval(fun,cdata,tdata)-cdata)^2
# T$ l7 F( S' N: c) ^! W5 t
plot(tdata,cdata)
1 m: p+ v8 g4 A
1 h+ Y5 ~. J- |/ 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