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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 请大侠帮忙看一下matlab利用lsqcurvefit曲线拟合问题
61/1返回列表
查看: 1506  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕306英一数二 +4 z1z2z3879 2026-03-16 4/200 2026-03-16 13:53 by laoshidan
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 梁成伟老师课题组欢迎你的加入 +6 一鸭鸭哟 2026-03-14 7/350 2026-03-15 22:12 by Winj1e
[考研] 材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 一志愿郑大070303，338分，求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 土木第一志愿276求调剂，科研和技能十分丰富，求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
信息提示
关闭
请填处理意见
关闭