liulian1122

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 81.3
帖子: 11
在線: 8.3小時(shí)
蟲號(hào): 3288386
注冊(cè): 2014-06-23
專業(yè): 大氣遙感和大氣探測(cè)

[求助] origin非線性擬合出錯(cuò)，求助已有2人參與

大家好，我用origin進(jìn)行非線性擬合（2個(gè)自變量，1個(gè)因變量，5個(gè)待求參數(shù)），設(shè)置了參數(shù)控制（均>0），但是出來的擬合結(jié)果中，有一個(gè)參數(shù)值是負(fù)值，不知道為什么，ps：本次只是用17組數(shù)據(jù)對(duì)origin擬合進(jìn)行測(cè)試
不知道是不是因?yàn)闇y(cè)試數(shù)據(jù)太少，求大神指點(diǎn)一番。
以下數(shù)據(jù)分別是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

擬合函數(shù)是y=(b+(a-b)*exp(-v1/p1))*(1-exp(-v2/p2))+c
參數(shù)控制：c>0;c<b<a;a<1;p1>0;p2>0;

擬合出的參數(shù)值如下：
a    0.7339
b    0.13606
c    -0.0476
p1 13.07059
p2 2.53515

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +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
相關(guān)系數(shù)(R): 0.544686868702746
相關(guān)系數(shù)之平方(R^2): 0.296683784937203
修正R平方(Adj. R^2): 0.134380042999634
確定系數(shù)(DC): 0.295959704287254
卡方系數(shù)(Chi-Square): 0.198174425821472
F統(tǒng)計(jì)(F-Statistic): 1.14478859461284
約束條件: b-c = 4.533449283E-014
b-a = -0.6370260749

參數(shù) 最佳估算
---------- -------------
a 0.637026074853859
b 5.27397200790168E-14
c 7.40522724589439E-15
p1 18.1979059607059
p2 2.67995326078984

贊一下

回復(fù)此樓

2樓2016-08-11 09:09:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liulian1122

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 81.3
帖子: 11
在線: 8.3小時(shí)
蟲號(hào): 3288386
注冊(cè): 2014-06-23
專業(yè): 大氣遙感和大氣探測(cè)

引用回帖:

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 ...

謝謝大神回復(fù)，正打算試試1stop
我用origin出現(xiàn)的問題，您知道是怎么回事嗎？有點(diǎn)想不通了。

贊一下

回復(fù)此樓

3樓2016-08-11 13:19:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

引用回帖:

3樓: Originally posted by liulian1122 at 2016-08-11 13:19:24
謝謝大神回復(fù)，正打算試試1stop
我用origin出現(xiàn)的問題，您知道是怎么回事嗎？有點(diǎn)想不通了。...

Origin擬合功能比1stOpt差很遠(yuǎn)，但也不至于出現(xiàn)范圍越界的情況，是否哪里操作有誤。

贊一下

回復(fù)此樓

4樓2016-08-11 13:56:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

引用回帖:

Origin擬合功能比1stOpt差很遠(yuǎn)，但也不至于出現(xiàn)范圍越界的情況，是否哪里操作有誤。

贊一下

回復(fù)此樓

5樓2016-08-11 16:35:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liulian1122

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 81.3
帖子: 11
在線: 8.3小時(shí)
蟲號(hào): 3288386
注冊(cè): 2014-06-23
專業(yè): 大氣遙感和大氣探測(cè)

引用回帖:

5樓: Originally posted by dingd at 2016-08-11 16:35:46
Origin擬合功能比1stOpt差很遠(yuǎn)，但也不至于出現(xiàn)范圍越界的情況，是否哪里操作有誤。...

新手實(shí)在是找不到哪里出了問題

贊一下

回復(fù)此樓

6樓2016-08-11 16:41:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

溫柔一刀5900

鐵桿木蟲 (正式寫手)

應(yīng)助: 15 (小學(xué)生)
金幣: 9105.3
散金: 116
紅花: 11
帖子: 871
在線: 651.6小時(shí)
蟲號(hào): 628256
注冊(cè): 2008-10-16
性別: GG
專業(yè): 生物、醫(yī)學(xué)光子學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

初始值設(shè)置的不好，你把這位仁兄給你設(shè)置的初始值帶進(jìn)去算算看，應(yīng)該不會(huì)錯(cuò)，另外數(shù)據(jù)擬合的質(zhì)量好不好可以看看最后的偏差。

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

liulian1122

7樓2016-08-11 21:57:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liulian1122

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 81.3
帖子: 11
在線: 8.3小時(shí)
蟲號(hào): 3288386
注冊(cè): 2014-06-23
專業(yè): 大氣遙感和大氣探測(cè)

送紅花一朵

引用回帖:

7樓: Originally posted by 溫柔一刀5900 at 2016-08-11 21:57:54
初始值設(shè)置的不好，你把這位仁兄給你設(shè)置的初始值帶進(jìn)去算算看，應(yīng)該不會(huì)錯(cuò)，另外數(shù)據(jù)擬合的質(zhì)量好不好可以看看最后的偏差。

是測(cè)試數(shù)據(jù)的問題，導(dǎo)致得不到好的擬合方案，謝謝各位大神的解答

贊一下

回復(fù)此樓

8樓2016-08-31 11:14:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 材料292調(diào)劑 +12	橘頌思美人 2026-03-23	12/600	2026-03-27 15:44 by caszguilin
[考研] 283求調(diào)劑（080500） +4	A child 2026-03-27	4/200	2026-03-27 15:34 by XPU李慶
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級(jí)，六級(jí)521，求調(diào)劑 +4	美味可樂雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 265求調(diào)劑 +4	小木蟲085600 2026-03-27	4/200	2026-03-27 14:26 by laoshidan
[考研] 化學(xué)308分求調(diào)劑 +8	你好明天你好 2026-03-23	9/450	2026-03-27 14:01 by 楊光于青云
[考研] 085601 材料工程 313分求調(diào)劑 +5	Ong3 2026-03-27	5/250	2026-03-27 12:24 by goldfish51
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +10	77chaselx 2026-03-20	10/500	2026-03-27 11:28 by wangjy2002
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +6	小小甜面團(tuán) 2026-03-25	6/300	2026-03-26 23:01 by 不吃魚的貓
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4	柳溪邊 2026-03-26	4/200	2026-03-26 21:14 by XPU李慶
[考研] 321求調(diào)劑 +6	wasdssaa 2026-03-26	6/300	2026-03-26 20:57 by sanrepian
[考研] 生物學(xué) 296 求調(diào)劑 +4	朵朵- 2026-03-26	6/300	2026-03-26 19:01 by 不吃魚的貓
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	7/350	2026-03-26 16:07 by 不吃魚的貓
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4	為學(xué)666 2026-03-24	4/200	2026-03-25 11:27 by BruceLiu320
[考研] 化學(xué)調(diào)劑 +6	yzysaa 2026-03-21	6/300	2026-03-25 09:27 by aa331100
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +4	這么名字咋樣 2026-03-22	5/250	2026-03-25 08:20 by mx.yue
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7	小羊睡著了? 2026-03-23	10/500	2026-03-25 02:24 by 清風(fēng)拂揚(yáng)。 m
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 306求調(diào)劑 +5	來好運(yùn)來來來 2026-03-22	5/250	2026-03-22 16:17 by BruceLiu320
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO