版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

caisinong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.3
帖子: 16
在線: 7.9小時
蟲號: 2514684
注冊: 2013-06-20
專業(yè): 自然語言理解與機器翻譯

[求助] 微分方程組擬合用1stopt，不是不響應(yīng)就是不開始。。大神們看下是哪里出錯了？已有1人參與

Title "Type your title here";
//Parameters ;
//Variable ;
//Function ;
Variable s,i,d;
Parameters a,b,c;
Function s'=-a*s*i;
d'=c*i;
i'=a*s*i-b*i-c*i;
RawData;
11744951,29,20;
11744897,66,37;
11744873,83,44;
11744843,101,56;
11744803,122,75;
11744792,136,72;
11744779,146,75;
11744769,155,76;
11744767,157,76;
11744742,174,84;
11744719,186,95;
11744656,215,129;
11744628,236,136;
11744602,264,134;
11744610,270,120;
11744587,303,110;
11744591,309,100;
11744585,314,101;
11744540,339,121;
11744515,358,127;
11744494,373,133;
11744457,394,149;
11744393,406,201;
11744229,494,277;
11744139,557,304;
11744058,601,341;
11743978,635,387;
11743843,710,447;
11743702,768,530;
11743528,843,629;
11743460,904,636;
11743402,981,617;
11743269,1041,690;
11743122,1142,736;
11743029,1192,779;
11742866,1260,874;
11742836,1327,837;
11742708,1428,864;
11742584,1525,891;
11742403,1607,990;
11742293,1708,999;
11742225,1781,994;
11742194,1814,992;
11742129,1876,995;
11742083,1910,1007;

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有296人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有72人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

用1stopt軟件做微分方程組的參數(shù)估計已經(jīng)有13人回復(fù)
這個編號后為什么1stopt沒有反應(yīng)，不能執(zhí)行已經(jīng)有5人回復(fù)
求助，用1stopt解方程組，兩次得到的解不唯一？什么原因，大神幫忙已經(jīng)有6人回復(fù)
1stopt4階非線性常微分方程，幫小弟看一下，跪求。已經(jīng)有7人回復(fù)
1stOpt求解常微分方程邊值問題已經(jīng)有5人回復(fù)

1樓 2015-02-09 04:23:58

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

回帖支持 ( 顯示支持度最高的前 50 名 )

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

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

★ ★ ★ ★ ★
caisinong: 金幣+5, ★★★★★最佳答案, 謝謝大神無私幫助！ 2015-02-09 12:39:41

原代碼中關(guān)鍵詞 RawData 也是錯誤的，正確的為 RowData。此外，方便起見，我做了變量命名轉(zhuǎn)換，因為原定義中因變量 i 可能與虛數(shù)產(chǎn)生歧義。
代碼如下：

CODE:

Parameters a,b,c;

Variable t,y1,y2,y3;

ODEFunction y1'=-a*y1*y2;

y3'=c*y2;

y2'=a*y1*y2-b*y2-c*y2;

Data;

1    11744951          29          20

           2    11744897          66          37

           3    11744873          83          44

           4    11744843         101          56

           5    11744803         122          75

           6    11744792         136          72

           7    11744779         146          75

           8    11744769         155          76

           9    11744767         157          76

          10    11744742         174          84

          11    11744719         186          95

          12    11744656         215         129

          13    11744628         236         136

          14    11744602         264         134

          15    11744610         270         120

          16    11744587         303         110

          17    11744591         309         100

          18    11744585         314         101

          19    11744540         339         121

          20    11744515         358         127

          21    11744494         373         133

          22    11744457         394         149

          23    11744393         406         201

          24    11744229         494         277

          25    11744139         557         304

          26    11744058         601         341

          27    11743978         635         387

          28    11743843         710         447

          29    11743702         768         530

          30    11743528         843         629

          31    11743460         904         636

          32    11743402         981         617

          33    11743269        1041         690

          34    11743122        1142         736

          35    11743029        1192         779

          36    11742866        1260         874

          37    11742836        1327         837

          38    11742708        1428         864

          39    11742584        1525         891

          40    11742403        1607         990

          41    11742293        1708         999

          42    11742225        1781         994

          43    11742194        1814         992

          44    11742129        1876         995

          45    11742083        1910        1007

計算結(jié)果：

CODE:

Parameters       Best Estimate

-------------------- -------------

a       1.48929558416951E-8

b       0.0120499576251067

c       0.0641192905655217

====== Output Results =====

File: Data file - 1

No Obs. y1 Cal. y1 Obs. y2 Cal. y2 Obs. y3 Cal. y3

11744897 11744945.668502 66 32.0098471641041 37 21.954366

2249076

11744873 11744939.7836624 83 35.3320768680856 44 2

1115717780767

11744843 11744933.2880524 101 38.9991096933077 56 2

4926687112677

11744803 11744926.1182829 122 43.0467308193201 75 2

1208939388352

11744792 11744918.2043863 136 47.5144391972183 72 3

021895998748

11744779 11744909.4691336 146 52.4458329259643 75 3

223985328777

11744769 11744899.8272811 155 57.8890346108615 76 3

7584105134752

11744767 11744889.1847386 157 63.8971608504808 76 4

6596631963224

11744742 11744877.4376513 174 70.5288404279365 84 4

9658146308117

11744719 11744864.4713866 186 77.8487862564284 95 5

7188871526054

11744656 11744850.1594159 215 85.9284266519954 129 5

9652641951778

11744628 11744834.3620802 236 94.8466020834877 136 6

7561428468951

11744602 11744816.9252277 264 104.690334186439 134 6

1480333618917

11744610 11744797.6787099 270 115.555674529899 120 7

2033104944093

11744587 11744776.4347221 303 127.548641400163 110 8

9908220308955

11744591 11744752.9859716 309 140.78625372005 100 92.586560

45099

11744585 11744727.103656 314 155.397672165231 101 102.07440

426294

11744540 11744698.5352315 339 171.52545857897 121 112.54693

6236917

11744515 11744667.0019504 358 189.326965933505 127 1

10634650769

11744494 11744632.1961429 373 208.975872351017 133 1

86542934334

11744457 11744593.7782171 394 230.663874091683 149 1

948683287116

11744393 11744551.3733482 406 254.602553953918 201 1

493525383908

11744229 11744504.567825 494 281.025443226915 277 183.65163

1310685

11744139 11744452.9050167 557 310.190297204049 304 2

590414452341

11744058 11744395.8809228 601 342.381606324876 341 2

494658312249

11743978 11744332.9392622 635 377.913367282602 387 2

568318133379

11743843 11744263.4660526 710 417.132140933892 447 2

036509245148

11743702 11744186.7836299 768 460.420426601673 530 3

147701458794

11743528 11744102.1440477 843 508.200385394373 629 3

176140827076

11743460 11744008.7217945 904 560.937948504155 636 3

424522281206

11743402 11743905.6057565 981 619.147350122214 617 4

226939079813

11743269 11743791.7903502 1041 683.396128653765 690 4

952137673838

11743122 11743666.1657379 1142 754.310644364555 736 4

007109532617

11743029 11743527.5070319 1192 832.582166483306 779 5

841054716886

11742866 11743374.4623843 1260 918.973588162441 874 5

949755554095

11742836 11743205.5398465 1327 1014.32683360825 837 6

880402702559

11742708 11743019.0928735 1428 1119.57102818053 864 7

236920219392

11742584 11742813.304333 1525 1235.7315093835 891 803.685841011576

11742403 11742586.1688676 1607 1363.93976448183 990 8

962789292442

11742293 11742335.4734398 1708 1505.44438903799 999 9

879632433323

11742225 11742058.7758766 1781 1661.62317004587 994 1

33237094253

11742194 11741753.3812077 1814 1833.99640759864 992 1

30984227545

11742129 11741416.3155756 1876 2024.24160024928 995 1

90332183955

11742083 11741044.2974683 1910 2234.20963147612 1007 1

3171129704

因為計算較慢，所以只迭代計算了幾次，看來效果是不錯的，見附圖。
至于計算出來的結(jié)果，比如參數(shù)取值范圍之類的，樓主可自行調(diào)節(jié)。

因為樓主的參數(shù)只有3個，所以可以用 1stOpt 高版本的試用版：
http://www.gaoyang168.com/bbs/viewthread.php?tid=8050245

微分方程組擬合用1stopt，不是不響應(yīng)就是不開始。。大神們看下是哪里出錯了？

附圖1.png

微分方程組擬合用1stopt，不是不響應(yīng)就是不開始。。大神們看下是哪里出錯了？-1

附圖2.png

微分方程組擬合用1stopt，不是不響應(yīng)就是不開始。。大神們看下是哪里出錯了？-2

附圖3.png

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

8樓2015-02-09 11:33:18

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

普通回帖

caisinong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.3
帖子: 16
在線: 7.9小時
蟲號: 2514684
注冊: 2013-06-20
專業(yè): 自然語言理解與機器翻譯

查了好多東西，發(fā)現(xiàn)應(yīng)該是因為自己的版本太低。。不知道高版本怎么獲��？網(wǎng)上都查不到。
或者如果有前輩可以幫忙跑下也好，不勝感激！

贊一下

回復(fù)此樓

2樓2015-02-09 04:38:21

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

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

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

1. 低版本不支持常微分方程擬合。
2. 代碼也存在錯誤。樓主定義的變量有s,i,d 三個，那么自變量是哪個呢？樓主定義的常微分方程有三個 s'=-a*s*i; d'=c*i; i'=a*s*i-b*i-c*i; 這樣看來，s i d都是因變量，那么缺少自變量的定義。在data中，只給出了三個變量的數(shù)據(jù)，缺少一組數(shù)據(jù)（應(yīng)該是自變量的數(shù)據(jù)）。
3. 常微分方程擬合，關(guān)鍵詞不再是 Function，而是ODEFunction。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

3樓2015-02-09 09:03:46

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

caisinong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.3
帖子: 16
在線: 7.9小時
蟲號: 2514684
注冊: 2013-06-20
專業(yè): 自然語言理解與機器翻譯

引用回帖:

3樓: Originally posted by 月只藍 at 2015-02-09 09:03:46
1. 低版本不支持常微分方程擬合。
2. 代碼也存在錯誤。樓主定義的變量有s,i,d 三個，那么自變量是哪個呢？樓主定義的常微分方程有三個 s'=-a*s*i; d'=c*i; i'=a*s*i-b*i-c*i; 這樣看來，s i d都是因變量，那么缺少 ...

謝謝回答。。我也是新手，本來用matlab做神經(jīng)網(wǎng)絡(luò)和模擬退火，效果都很差，您可以看下我另一篇求助貼里有你和效果圖。自變量就是t，昨天趕得太急沒有細看，我再看看

贊一下

回復(fù)此樓

4樓2015-02-09 10:12:46

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

引用回帖:

4樓: Originally posted by caisinong at 2015-02-09 10:12:46
謝謝回答。。我也是新手，本來用matlab做神經(jīng)網(wǎng)絡(luò)和模擬退火，效果都很差，您可以看下我另一篇求助貼里有你和效果圖。自變量就是t，昨天趕得太急沒有細看，我再看看...

自變量是t，也需要給出t的數(shù)據(jù)。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

5樓2015-02-09 10:21:56

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

caisinong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.3
帖子: 16
在線: 7.9小時
蟲號: 2514684
注冊: 2013-06-20
專業(yè): 自然語言理解與機器翻譯

引用回帖:

5樓: Originally posted by 月只藍 at 2015-02-09 10:21:56
自變量是t，也需要給出t的數(shù)據(jù)。...

t就是時間序列，添加了這一列的數(shù)據(jù)為：
1,11744951,29,20;
2,11744897,66,37;
3,11744873,83,44;
4,11744843,101,56;
5,11744803,122,75;
6,11744792,136,72;
7,11744779,146,75;
8,11744769,155,76;
9,11744767,157,76;
10,11744742,174,84;
11,11744719,186,95;
12,11744656,215,129;
13,11744628,236,136;
14,11744602,264,134;
15,11744610,270,120;
16,11744587,303,110;
17,11744591,309,100;
18,11744585,314,101;
19,11744540,339,121;
20,11744515,358,127;
21,11744494,373,133;
22,11744457,394,149;
23,11744393,406,201;
24,11744229,494,277;
25,11744139,557,304;
26,11744058,601,341;
27,11743978,635,387;
28,11743843,710,447;
29,11743702,768,530;
30,11743528,843,629;
31,11743460,904,636;
32,11743402,981,617;
33,11743269,1041,690;
34,11743122,1142,736;
35,11743029,1192,779;
36,11742866,1260,874;
37,11742836,1327,837;
38,11742708,1428,864;
39,11742584,1525,891;
40,11742403,1607,990;
41,11742293,1708,999;
42,11742225,1781,994;
43,11742194,1814,992;
44,11742129,1876,995;
45,11742083,1910,1007;

我這里低版本看來是不行了，能麻煩您跑下嗎？如果驗證這個可行，我再找辦法跑別的數(shù)據(jù)。。

贊一下

回復(fù)此樓

6樓2015-02-09 10:35:37

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

caisinong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.3
帖子: 16
在線: 7.9小時
蟲號: 2514684
注冊: 2013-06-20
專業(yè): 自然語言理解與機器翻譯

引用回帖:

5樓: Originally posted by 月只藍 at 2015-02-09 10:21:56
自變量是t，也需要給出t的數(shù)據(jù)。...

s',i',d'分別為ds/dt,di/dt,dd/dt

贊一下

回復(fù)此樓

7樓2015-02-09 10:36:42

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

caisinong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.3
帖子: 16
在線: 7.9小時
蟲號: 2514684
注冊: 2013-06-20
專業(yè): 自然語言理解與機器翻譯

引用回帖:

8樓: Originally posted by 月只藍 at 2015-02-09 11:33:18
原代碼中關(guān)鍵詞 RawData 也是錯誤的，正確的為 RowData。此外，方便起見，我做了變量命名轉(zhuǎn)換，因為原定義中因變量 i 可能與虛數(shù)產(chǎn)生歧義。
代碼如下：

Parameters a,b,c;
Variable t,y1,y2,y3;
ODEFunction ...

先感謝您了，擬合效果很好，就是還沒太懂參數(shù)為什么跟預(yù)期差那么大。我自己再試試！

贊一下

回復(fù)此樓

9樓2015-02-09 12:41:03

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

caisinong

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 4.3
帖子: 16
在線: 7.9小時
蟲號: 2514684
注冊: 2013-06-20
專業(yè): 自然語言理解與機器翻譯

引用回帖:

有點郁悶。。用那個試用版，直接復(fù)制您給的代碼也不跑，不知道為什么

贊一下

回復(fù)此樓

10樓2015-02-09 13:04:39

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 caisinong 的主題更新

返回列表

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

[考研] 277材料科學(xué)與工程080500求調(diào)劑 +7	自由煎餅果子 2026-03-16	7/350	2026-03-22 22:40 by ACS Nano——
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +3	jiajunX 2026-03-22	3/150	2026-03-22 19:32 by brblmd
[考研] 269專碩求調(diào)劑 +6	金恩貝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[考研] 318求調(diào)劑 +4	plum李子 2026-03-21	7/350	2026-03-22 14:17 by ColorlessPI
[考研] 286求調(diào)劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[基金申請] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 材料學(xué)求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷

24小時熱門版塊排行榜

[求助] 微分方程組擬合用1stopt，不是不響應(yīng)就是不開始。。大神們看下是哪里出錯了？ 已有1人參與

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

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

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

[求助] 微分方程組擬合用1stopt，不是不響應(yīng)就是不開始。。大神們看下是哪里出錯了？已有1人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: