应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 请大侠帮忙看一下matlab利用lsqcurvefit曲线拟合问题
61/1返回列表
查看: 1674  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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-')

不知道以上为啥不对？各位大侠有没有别的合适 的曲线拟合的方法？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
加油加油~
1楼 2012-11-30 11:25:35
已阅   回复此楼   关注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 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 14:27 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 14:24 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:36 by LH5NkK5Ud4bh
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 05:51 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 02:23 by DgNGHc3h5tPl
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:59 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
信息提示
关闭
请填处理意见
关闭