数学建模社区-数学中国
标题:
matlab拟合问题~不知道错在哪里~
[打印本页]
作者:
guizhanbu
时间:
2009-3-22 23:14
标题:
matlab拟合问题~不知道错在哪里~
fun=inline('x(1)./(1+(x(1)./3315.84-1).*exp(-x(2).*tdata))','x','tdata');
7 J$ }4 J8 t: q- N0 Y5 X1 b
tdata=1970:1:2008
% \- S7 p' Q( K. G
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]]
6 U0 j* R2 ] Z* V" _: J0 y6 X
x=lsqcurvefit(fun,[5000,0.2],tdata,cdata)
0 K8 S; u; r' X2 P/ b
f=fun(x,tdata)
$ ^. k# {. v p; f7 V# I: \
norm(feval(fun,cdata,tdata)-cdata)^2
- k4 b! \+ ?1 y
plot(tdata,cdata)
2 H8 ~4 m5 C" @5 F0 V# e+ ^' ?
8 {! a3 B( ?# N/ Y3 s! A2 {
为什么怎么改变,[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