24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1994  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 化学专业申博 +3 赵子羊 2026-05-23 4/200 2026-05-24 18:10 by 工大学长
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 1rx34o113h 2026-05-23 3/150 2026-05-24 17:41 by 0i3mu4vkjz
[基金申请] 评审有感 +16 popular289 2026-05-18 27/1350 2026-05-24 17:34 by hhs666
[教师之家] 论文撤稿了 +4 bjvtcliu 2026-05-24 7/350 2026-05-24 17:29 by bjvtcliu
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:10 by 75ui6h7z2t
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 17:01 by 75ui6h7z2t
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 5/250 2026-05-24 16:21 by hhx1yx9evi
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 4/200 2026-05-24 16:16 by hhx1yx9evi
[基金申请] 河北省自然科学基金 +6 Peterchao 2026-05-18 9/450 2026-05-24 16:02 by 130067131
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 pmo95bazuy 2026-05-23 8/400 2026-05-24 15:56 by 1uy1ht2y9r
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[论文投稿] 投稿求助,期刊 +4 希冀,有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀,有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
请填处理意见