24小时热门版块排行榜    

查看: 1335  |  回复: 5

taiyanghef

铁虫 (著名写手)


[求助] origin如何非线性拟合这个公式呢? 已有1人参与

本人用origin9.1非线性拟合上述公式,式中有7个待求参数A1-A7,一个自变量x,一个因变量Y,控制参数是:A1和A5>0,但是我自己求解出来的拟合过程中总量出错,不知道为什么,求大神指点一番。

以下数据分别是x和Y的10组实验数据:
1        32.52
2        43.03
3        48.04
4        51.26
5        53.65
6        56.35
7        58.6
8        59.75
9        61.7
10        62.95

origin如何非线性拟合这个公式呢?
11.png
回复此楼

» 收录本帖的淘帖专辑推荐

ScientificPlots

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ferryshaobo

新虫 (初入文坛)

你这个参数太多,拟合是根据最小化定理拟合的。你的点只有十个而待求参数多达七个,最好的办法是将系数的范伟进一步缩小,先固定五个参数求另外两个否则很难求。如果方便把系数变动范围收小,我可以帮你算一下。

发自小木虫Android客户端
2楼2016-08-17 12:47:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
taiyanghef(月只蓝代发): 金币+5, 感谢指导! 2016-08-17 21:05:46
taiyanghef: 金币+5, ★★★★★最佳答案, 谢谢你的帮助,我下午也用你这软件拟合了,跟你的拟合数据有很大差别,还有就是你的拟合前提是所有参数都大于0了,我只要求参数1和5是大于0的,不过还是很感谢你,另外,你这最后给出的a1-a3最佳估算是怎么确定的?一会我重新发了个帖子,麻烦帮我再看看,谢谢了 2016-08-17 21:47:18
1stOpt求解:
CODE:
ParameterDomain = [0,];
Function y=(A1*150^A2*x^A3*exp(-A4/190))/(A3+1)+A5*150^A6*x*exp(-A7/190);
Data;
1        32.52
2        43.03
3        48.04
4        51.26
5        53.65
6        56.35
7        58.6
8        59.75
9        61.7
10        62.95

公式有过拟合现象,目标函数(最小二乘)一样但参数组值不唯一:
1:
均方差(RMSE): 1.09722266644578
残差平方和(SSE): 12.0389757976239
相关系数(R): 0.992592323979301
相关系数之平方(R^2): 0.985239521622631
修正R平方(Adj. R^2): 0.981022242086239
确定系数(DC): 0.985123443847963
卡方系数(Chi-Square): 0.151169189754915
F统计(F-Statistic): 32.4087280049164

参数        最佳估算
----------        -------------
a1        19.8887680296124
a2        0.184915005560769
a3        0.260451660623446
a4        24.2481065803997
a5        1.70095357493163E-15
a6        0.0455380528899979
a7        20.1615099619736

2:
均方差(RMSE): 1.09722266644578
残差平方和(SSE): 12.0389757976239
相关系数(R): 0.992592323979793
相关系数之平方(R^2): 0.985239521623606
修正R平方(Adj. R^2): 0.981022242087494
确定系数(DC): 0.985123443847963
卡方系数(Chi-Square): 0.151169189697482
F统计(F-Statistic): 32.4087279807105

参数        最佳估算
----------        -------------
a1        3.71448551969752
a2        0.507884339036319
a3        0.260451660610369
a4        12.9177550813263
a5        1.65575495639403E-15
a6        0.0233122548072396
a7        3.34121107925197

其实公式可简化为:y=A1*x^A2+A3*x;
效果完全一样,参数组值唯一:
均方差(RMSE): 1.09722266644578
残差平方和(SSE): 12.0389757976239
相关系数(R): 0.9925923240161
相关系数之平方(R^2): 0.985239521695682
修正R平方(Adj. R^2): 0.981022242180163
确定系数(DC): 0.985123443847963
卡方系数(Chi-Square): 0.151169191013433
F统计(F-Statistic): 226.861094248425

参数        最佳估算
----------        -------------
a1        35.0790846880362
a2        0.260451659644489
a3        1.94918389396105E-15
3楼2016-08-17 14:49:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

taiyanghef

铁虫 (著名写手)


引用回帖:
2楼: Originally posted by ferryshaobo at 2016-08-17 12:47:14
你这个参数太多,拟合是根据最小化定理拟合的。你的点只有十个而待求参数多达七个,最好的办法是将系数的范伟进一步缩小,先固定五个参数求另外两个否则很难求。如果方便把系数变动范围收小,我可以帮你算一下。
...

真心谢谢你的乐助,我用另个软件1stopt拟合出来了,但是有多组拟合结果,都不唯一,不知道如何确定最佳参数,我现在重新发了个帖子,希望能帮忙看看,谢谢了
4楼2016-08-17 21:44:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

taiyanghef

铁虫 (著名写手)


引用回帖:
3楼: Originally posted by dingd at 2016-08-17 14:49:05
1stOpt求解:

ParameterDomain = ;
Function y=(A1*150^A2*x^A3*exp(-A4/190))/(A3+1)+A5*150^A6*x*exp(-A7/190);
Data;
1        32.52
2        43.03
3        48.04
4        51.26
5        53.65
...

我明白你的a1-a3的意思了,是简化公式的结果,我也想过简化,但是我需要的是上述公式的七个参数,如果只有这简化中的三个公式,又要反求另外4个,肯定会带来一定误差,所以直接拟合
5楼2016-08-17 22:19:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

taiyanghef

铁虫 (著名写手)


引用回帖:
5楼: Originally posted by taiyanghef at 2016-08-17 22:19:52
我明白你的a1-a3的意思了,是简化公式的结果,我也想过简化,但是我需要的是上述公式的七个参数,如果只有这简化中的三个公式,又要反求另外4个,肯定会带来一定误差,所以直接拟合...

对了,你的拟合公式里面,是x^(a3+1),不是x^a3
6楼2016-08-17 22:29:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 taiyanghef 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856材料专硕353求调剂 +4 NIFFFfff 2026-03-20 4/200 2026-03-22 09:49 by 2026paper
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
请填处理意见