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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 请大侠帮忙看一下matlab利用lsqcurvefit曲线拟合问题
61/1返回列表
查看: 1505  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 一志愿安徽大学材料工程专硕313分，求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 求调剂，一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 一志愿中科院，化学方向，295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 一志愿西南交大，材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] （081700）化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 283求调剂，材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
信息提示
关闭
请填处理意见
关闭