24小时热门版块排行榜

返回列表

liulian1122

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 81.3
帖子: 11
在线: 8.3小时
虫号: 3288386
注册: 2014-06-23
专业: 大气遥感和大气探测

[求助] 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

回复此楼

1楼2016-08-10 22:21:17

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

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

感谢参与，应助指数 +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

205          1        1.11707

258   2        1.23767

222          3        1.37491

19          4        2.40285

247          5        3.55603

313   6        3.60568

721          7        5.72943

348   8        12.68831

376          9        12.51504

404          10        13.68702

25          11        15.00665

279          12        15.00868

277          13        15.02606

327          14        16.51203

223          15        17.62735

326          16        18.26764

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

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 81.3
帖子: 11
在线: 8.3小时
虫号: 3288386
注册: 2014-06-23
专业: 大气遥感和大气探测

引用回帖:

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的回帖