當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

hua99dushu

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 74.1
帖子: 161
在線: 25.5小時(shí)
蟲號(hào): 2078718
注冊(cè): 2012-10-22
性別: MM
專業(yè): 控制論中的數(shù)學(xué)方法

[求助] 用1stopt軟件做微分方程組的參數(shù)估計(jì) 已有1人參與

近一個(gè)月里才知道1stopt這個(gè)軟件，可以做微分方程組的參數(shù)估計(jì)問(wèn)題。
之前一直想用matlab的神經(jīng)網(wǎng)絡(luò)來(lái)做，由于水平有限，也做不出來(lái)。
網(wǎng)上下載的1stopt軟件大多是無(wú)法寫代碼或者未注冊(cè)，根據(jù)一些實(shí)例寫了我的代碼，可能有錯(cuò)誤。
請(qǐng)幫忙改正后，跑一下軟件。

Title "test 1";
Parameters a1,a2,a3,b1,b2,b3,c1,c2,c3,M,C,E,N[0,];
Sharedmodel;
Variable t, x,y,z;
Function x'=a1*x*(y/M-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E);
z'=c1*x*(x/N-1)-c2*y-c3*z;
RowData;
0,1,2,3,4,5,6,7,8,9,10;
265,291,325,378,453,557,647,770,911,1006,1201;
7217,7344,7800,8796,10607,13266,14893,16219,17303,17677,19786;
8554,9457,10607,12443,15004,18599,21742,26018,30981,34457,41425;

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有131人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有27人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有47人回復(fù)

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

尋求方程組的解法，急急急��！已經(jīng)有7人回復(fù)
origin曲線擬合求參數(shù) 出錯(cuò) 已經(jīng)有7人回復(fù)
請(qǐng)教一個(gè)二元回歸問(wèn)題，求高手解決，謝謝。已經(jīng)有4人回復(fù)
求助，matlab/ 1stOpt 自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來(lái)求方程中的參數(shù) 已經(jīng)有4人回復(fù)
兩個(gè)自變量怎么擬合求參數(shù)值呢？已經(jīng)有7人回復(fù)
解帶參數(shù)方程，參數(shù)在變已經(jīng)有6人回復(fù)
微分方程組參數(shù)擬合的問(wèn)題（多參數(shù)動(dòng)力系統(tǒng)參數(shù)估計(jì)）已經(jīng)有11人回復(fù)
求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有9人回復(fù)
關(guān)于matlab微分、及非線性擬合問(wèn)題已經(jīng)有5人回復(fù)
MATLAB微分方程參數(shù)擬合問(wèn)題，求大神已經(jīng)有7人回復(fù)
1stopt4階非線性常微分方程，幫小弟看一下，跪求。已經(jīng)有7人回復(fù)
matlab求解方程中的參數(shù) 已經(jīng)有21人回復(fù)
形如1-(1-x)^a=by的方程怎么做擬合? 已經(jīng)有11人回復(fù)
請(qǐng)問(wèn)matlab擬合二元方程的問(wèn)題已經(jīng)有3人回復(fù)
一階微分方程求解，并優(yōu)化參數(shù)值已經(jīng)有14人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較已經(jīng)有4人回復(fù)
高金求助matlab解微分方程組已經(jīng)有12人回復(fù)
求回歸方程參數(shù)估計(jì)的最優(yōu)化算法matlab代碼已經(jīng)有9人回復(fù)
【緊急求助】帶有微分平方項(xiàng)的ODE求解已經(jīng)有11人回復(fù)
重金加人民幣求擬合wilson方程和NRTL方程參數(shù)的程序已經(jīng)有17人回復(fù)
【求助】求助一個(gè)一階常微分方程的初值問(wèn)題已經(jīng)有16人回復(fù)

1樓 2014-02-10 13:32:35

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

hua99dushu

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 74.1
帖子: 161
在線: 25.5小時(shí)
蟲號(hào): 2078718
注冊(cè): 2012-10-22
性別: MM
專業(yè): 控制論中的數(shù)學(xué)方法

引用回帖:

3樓: Originally posted by dingd at 2014-02-10 22:03:32
感覺(jué)你的公式有些問(wèn)題，似乎可以簡(jiǎn)化，下面結(jié)果參考下：

均方差(RMSE):666.824975358025
殘差平方和(SSE):13339666.4328369
相關(guān)系數(shù)(R): 0.99281958012936
相關(guān)系數(shù)之平方(R^2): 0.985690718688239
確定系數(shù) ...

此結(jié)果是唯一的嗎？看擬合結(jié)果是有點(diǎn)問(wèn)題，可以設(shè)定參數(shù)的其他約束條件嗎？

贊一下

回復(fù)此樓

4樓2014-02-11 10:18:42

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

查看全部 14 個(gè)回答

hua99dushu

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 74.1
帖子: 161
在線: 25.5小時(shí)
蟲號(hào): 2078718
注冊(cè): 2012-10-22
性別: MM
專業(yè): 控制論中的數(shù)學(xué)方法

所有的參數(shù)都大于0，不知道程序那么定義參數(shù)對(duì)不對(duì)。。

贊一下

回復(fù)此樓

2樓2014-02-10 13:54:34

已閱回復(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
hua99dushu: 金幣+30, ★★★很有幫助, 辛苦了，能否幫忙再做一下？付報(bào)酬。 2014-02-13 12:37:55
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-02-17 12:03:19

感覺(jué)你的公式有些問(wèn)題，似乎可以簡(jiǎn)化，下面結(jié)果參考下：

均方差(RMSE):666.824975358025
殘差平方和(SSE):13339666.4328369
相關(guān)系數(shù)(R): 0.99281958012936
相關(guān)系數(shù)之平方(R^2): 0.985690718688239
確定系數(shù)(DC): 0.752302595710559
F統(tǒng)計(jì)(F-Statistic): -0.355256071121647

參數(shù) 最佳估算
-------------------- -------------
a1 0.063244391232729
a2 26.6537700533949
a3 425.865501005945
b1 3256273550.99063
b2 1385.92982355728
b3 190853579.985187
c1 1967.65266459696
c2 7557.60987652871
c3 174107.388897955
m 24291.4733546288
c 0.10182112015089
e 1.27923766610631E22
n 0.0972369125973638

文件: 數(shù)據(jù)文件-1
No t 目標(biāo) x 計(jì)算 x 目標(biāo) y 計(jì)算 y 目標(biāo) z 計(jì)算 z
1 1 291 295.919384351006 7344 7409.69953342258 9457 8964.92474439791
2 2 325 329.82579568465 7800 8490.90903548783 10607 10769.8722203566
3 3 378 365.520171155061 8796 9774.32316004032 12443 13050.8936628585
4 4 453 403.287476232365 10607 11128.0651594216 15004 15713.5200315771
5 5 557 443.400579346621 13266 12544.2720271355 18599 18789.0601075391
6 6 647 486.105229396296 14893 14028.6589510091 21742 22331.0481756248
7 7 770 531.633932176283 16219 15588.082085985 26018 26401.9444748202
8 8 911 580.212516031404 17303 17228.8956295004 30981 31071.3625784165
9 9 1006 632.066439410749 17677 18957.0942374736 34457 36416.4169412114
10 10 1201 687.422431267362 19786 20778.230707146 41425 42522.166618056

贊一下(1人)

回復(fù)此樓

3樓2014-02-10 22:03:32

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

hua99dushu

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 74.1
帖子: 161
在線: 25.5小時(shí)
蟲號(hào): 2078718
注冊(cè): 2012-10-22
性別: MM
專業(yè): 控制論中的數(shù)學(xué)方法

★
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-02-17 12:03:30

引用回帖:

4樓: Originally posted by hua99dushu at 2014-02-11 10:18:42
此結(jié)果是唯一的嗎？看擬合結(jié)果是有點(diǎn)問(wèn)題，可以設(shè)定參數(shù)的其他約束條件嗎？...

Title "test 1";
Parameters a1[0,],a2[0,],a3[0,],b1[0,],b2[0,],b3[0,],c1[0,],c2[0,],c3[0,],M[8000,],C[y,],E[z,],N[300,];
Sharedmodel;
Variable t, x,y,z;
Function x'=a1*x*(y/M-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E);
z'=c1*x*(x/N-1)-c2*y-c3*z;
RowData;
0,1,2,3,4,5,6,7,8,9,10;
265,291,325,378,453,557,647,770,911,1006,1201;
7217,7344,7800,8796,10607,13266,14893,16219,17303,17677,19786;
8554,9457,10607,12443,15004,18599,21742,26018,30981,34457,41425;

贊一下

回復(fù)此樓

5樓2014-02-11 10:31:18

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

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

[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[教師之家] 焦慮 +8	水冰月月野兔 2026-03-13	11/550	2026-03-18 13:10 by 水冰月月野兔
[考研] 266求調(diào)劑 +4	陽(yáng)陽(yáng)哇塞 2026-03-14	8/400	2026-03-18 12:43 by Linda Hu
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛(ài)生物生物愛(ài)?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 268求調(diào)劑 +6	簡(jiǎn)單點(diǎn)0 2026-03-17	6/300	2026-03-18 09:04 by 無(wú)際的草原
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 085600材料與化工 +4	安全上岸！ 2026-03-16	4/200	2026-03-17 14:02 by 勇敢太監(jiān)王公公
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 270求調(diào)劑 085600材料與化工專碩 +3	YXCT 2026-03-11	3/150	2026-03-13 10:13 by houyaoxu