24小时热门版块排行榜    

查看: 1482  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思:在爱情中有过遗憾吗? +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错,点赞! +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题,请指教! 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
请填处理意见