24小时热门版块排行榜    

查看: 2177  |  回复: 7

liulian1122

新虫 (初入文坛)

[求助] origin非线性拟合出错,求助 已有2人参与

大家好,我用origin进行非线性拟合(2个自变量,1个因变量,5个待求参数),设置了参数控制(均>0),但是出来的拟合结果中,有一个参数值是负值,不知道为什么,ps:本次只是用17组数据对origin拟合进行测试
不知道是不是因为测试数据太少,求大神指点一番。
以下数据分别是y,v1,v2
0.205          1        1.11707
0.258   2        1.23767
0.222          3        1.37491
0.19          4        2.40285
0.247          5        3.55603
0.313   6        3.60568
0.721          7        5.72943
0.348   8        12.68831
0.376          9        12.51504
0.404          10        13.68702
0.25          11        15.00665
0.279          12        15.00868
0.277          13        15.02606
0.327          14        16.51203
0.223          15        17.62735
0.326          16        18.26764
0.267          17        19.27053

拟合函数是y=(b+(a-b)*exp(-v1/p1))*(1-exp(-v2/p2))+c
参数控制:c>0;c<b<a;a<1;p1>0;p2>0;

拟合出的参数值如下:
a     0.7339
b     0.13606
c     -0.0476
p1   13.07059
p2   2.53515
回复此楼

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

ScientificPlots

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
1stOpt求解:
CODE:
Parameter [a,b]=[0,1],c>0,p1>0,p2>0;
Variable y,v1,v2;
Function  y=(b+(a-b)*exp(-v1/p1))*(1-exp(-v2/p2))+c;
          c<b<a;
Data;
//y,v1,v2
0.205          1        1.11707
0.258   2        1.23767
0.222          3        1.37491
0.19          4        2.40285
0.247          5        3.55603
0.313   6        3.60568
0.721          7        5.72943
0.348   8        12.68831
0.376          9        12.51504
0.404          10        13.68702
0.25          11        15.00665
0.279          12        15.00868
0.277          13        15.02606
0.327          14        16.51203
0.223          15        17.62735
0.326          16        18.26764
0.267          17        19.27053

均方差(RMSE): 0.0994761927092236
残差平方和(SSE): 0.168223719570684
相关系数(R): 0.544686868702746
相关系数之平方(R^2): 0.296683784937203
修正R平方(Adj. R^2): 0.134380042999634
确定系数(DC): 0.295959704287254
卡方系数(Chi-Square): 0.198174425821472
F统计(F-Statistic): 1.14478859461284
约束条件: b-c = 4.533449283E-014
          b-a = -0.6370260749

参数        最佳估算
----------        -------------
a        0.637026074853859
b        5.27397200790168E-14
c        7.40522724589439E-15
p1        18.1979059607059
p2        2.67995326078984
2楼2016-08-11 09:09:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liulian1122

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2016-08-11 09:09:01
1stOpt求解:

Parameter =,c>0,p1>0,p2>0;
Variable y,v1,v2;
Function  y=(b+(a-b)*exp(-v1/p1))*(1-exp(-v2/p2))+c;
          c<b<a;
Data;
//y,v1,v2
0.205          1        1.11707 ...

谢谢大神回复,正打算试试1stop
我用origin出现的问题,您知道是怎么回事吗?有点想不通了。
3楼2016-08-11 13:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
3楼: Originally posted by liulian1122 at 2016-08-11 13:19:24
谢谢大神回复,正打算试试1stop
我用origin出现的问题,您知道是怎么回事吗?有点想不通了。...

Origin拟合功能比1stOpt差很远,但也不至于出现范围越界的情况,是否哪里操作有误。
4楼2016-08-11 13:56:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
3楼: Originally posted by liulian1122 at 2016-08-11 13:19:24
谢谢大神回复,正打算试试1stop
我用origin出现的问题,您知道是怎么回事吗?有点想不通了。...

Origin拟合功能比1stOpt差很远,但也不至于出现范围越界的情况,是否哪里操作有误。
5楼2016-08-11 16:35:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liulian1122

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by dingd at 2016-08-11 16:35:46
Origin拟合功能比1stOpt差很远,但也不至于出现范围越界的情况,是否哪里操作有误。...

新手实在是找不到哪里出了问题
6楼2016-08-11 16:41:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

温柔一刀5900

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
初始值设置的不好,你把这位仁兄给你设置的初始值带进去算算看,应该不会错,另外数据拟合的质量好不好可以看看最后的偏差。

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

7楼2016-08-11 21:57:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liulian1122

新虫 (初入文坛)

送红花一朵
引用回帖:
7楼: Originally posted by 温柔一刀5900 at 2016-08-11 21:57:54
初始值设置的不好,你把这位仁兄给你设置的初始值带进去算算看,应该不会错,另外数据拟合的质量好不好可以看看最后的偏差。

是测试数据的问题,导致得不到好的拟合方案,谢谢各位大神的解答
8楼2016-08-31 11:14:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liulian1122 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见