24小时热门版块排行榜    

查看: 1375  |  回复: 5

gaojinggigi

新虫 (正式写手)

[求助] 请大侠帮忙看一下matlab利用lsqcurvefit曲线拟合问题

function f = curvefun1(a,lambda)
f=1+a(1)./(lambda.^2-a(2))+a(3)./(lambda.^2-a(4))+a(5)./(lambda.^2-a(6))+a(7)./(lambda.^2-a(8));
end

lambda=0.45:0.05:1.55;
ndata=[2.96005 2.81205 2.69584 2.62716 2.58355 2.55385 2.53057 2.51262 2.4984 2.48737 2.47786 2.46993 2.46352 2.4578 2.45288 2.44862 2.44506 2.44179 2.43901 2.43643 2.43413 2.43216 2.4303];
a0=ones(8,1);
a=lsqcurvefit('curvefun1',a0,lambda,ndata);
f=curvefun1(a,lambda);
plot(lambda,ndata,'*',lambda,f,'r-')

不知道以上为啥不对?各位大侠有没有别的合适 的曲线拟合的方法?
回复此楼
加油加油~
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
fegg7502: 金币+2, 鼓励交流 2012-11-30 15:52:18
gaojinggigi: 金币+1, 有帮助 2012-12-04 18:03:58
a=lsqcurvefit(@curvefun1,a0,lambda,ndata);
另外,你的lambda=1时,分母为零,会报错的。
showmethemoney
2楼2012-11-30 14:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
fegg7502: 金币+2, 鼓励交流 2012-11-30 15:52:29
gaojinggigi: 金币+1, ★★★很有帮助 2012-12-04 18:04:08
下面的结果供参考,用1stOpt计算的:

均方差(RMSE): 0.000173297428049338
残差平方和(SSE): 6.90735967075853E-7
相关系数(R): 0.999999129682542
相关系数之平方(R^2): 0.999998259365841
决定系数(DC): 0.999998259365841
卡方系数(Chi-Square): 1.35117472958674E-7
F统计(F-Statistic): 1231078.41564268

参数        最佳估算
----------        -------------
a1        0.0178095668462072
a2        0.271198708521368
a3        -0.0178898339024012
a4        0.270640309966547
a5        3315.37310259053
a6        -2356.99656143336
a7        0.057829693695241
a8        0.0974042234170173
3楼2012-11-30 14:47:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaojinggigi

新虫 (正式写手)

引用回帖:
2楼: Originally posted by csgt0 at 2012-11-30 14:32:17
a=lsqcurvefit(@curvefun1,a0,lambda,ndata);
另外,你的lambda=1时,分母为零,会报错的。

那初始值应该怎么选呢?比如我把a0全改成-1,拟合出来差的也好大啊,怎么办呢?
加油加油~
4楼2012-12-01 13:33:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaojinggigi

新虫 (正式写手)

引用回帖:
3楼: Originally posted by dingd at 2012-11-30 14:47:05
下面的结果供参考,用1stOpt计算的:

均方差(RMSE): 0.000173297428049338
残差平方和(SSE): 6.90735967075853E-7
相关系数(R): 0.999999129682542
相关系数之平方(R^2): 0.999998259365841
决定系数(DC): 0 ...

书上的结果是1.612 0.0284 1.922 0.0571 0.884 0.110  0.127 0.209
虽然我的实验数据跟他的可能有点差距,但是我觉得我的趋势不对,我换了几种材料的实验结果,拟合结果都不对呢,用origin做的拟合fit
did not converge,唉
加油加油~
5楼2012-12-01 13:39:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

你将参数结果带回去验证一下就知道孰对孰错了。
6楼2012-12-01 15:22:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gaojinggigi 的主题更新
信息提示
请填处理意见