24小时热门版块排行榜    

查看: 1519  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +3 CrisMessi 2026-08-18 3/150 2026-08-18 23:38 by 必须成功
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 27/1350 2026-08-18 21:44 by Noways
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见