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

hua99dushu

鐵蟲 (小有名氣)

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

[求助] 求版主幫忙跑下面的1stopt的程序，參數(shù)估計(jì)。已有1人參與

[/code]Title "test 1";
Parameters a1[0,],a2[0,],a3[0,],b1[0,],b2[0,],b3[0,],c1[0,],c2[0,],c3[0,],M[1.5,],C[y,],E[z,],N[1.5,];
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：
x=[1.0000,1.0777,1.1443,1.3497,1.5881,1.9814,2.1905,2.5799,2.9763,3.2445,3.5765,3.9315];
y=[1.0000,1.0597,1.0982,1.2035,1.3566,1.4610,1.4771,1.4862,1.6052,1.6094,1.8892,2.0334];
z=[1.0000,1.0920,1.2033,1.4031,1.6920,1.9382,2.2159,2.6186,2.9489,3.1536,3.5979,4.0233];

感謝小木蟲論壇，感謝斑竹們的耐心幫忙。

回復(fù)此樓

» 猜你喜歡

重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有159人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有28人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有21人回復(fù)

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

跪求1stOpt2.5以上版本，或者請高人幫忙擬合包含復(fù)數(shù)的函數(shù)的參數(shù) 已經(jīng)有8人回復(fù)
求一個(gè)真真注冊版的1stOpt 已經(jīng)有5人回復(fù)

1樓 2014-02-20 10:38:45

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

hua99dushu

鐵蟲 (小有名氣)

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

引用回帖:

2樓: Originally posted by dingd at 2014-02-20 20:40:16
參考：

Parameters a1=,a2=,a3=,b1=,b2=,b3=,c1=,c2=,c3=,M=,C=,E=,N=;
Variable t,x,y,z;
ODEFunction 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;
Dat ...

之前程序里參數(shù)的設(shè)置我沒有講清楚。
是這樣的：M[1.5,],C[1,],E[1,],N[1.5,].
幫忙再來一遍。。急著對比結(jié)果啊

贊一下

回復(fù)此樓

4樓2014-02-24 10:09:35

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

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

dingd

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

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
hua99dushu: 金幣+30, ★★★很有幫助, 麻煩斑竹啦，版主好人。。 2014-02-24 10:06:49

參考：

CODE:

Parameters a1=[0,],a2=[0,],a3=[0,],b1=[0,],b2=[0,],b3=[0,],c1=[0,],c2=[0,],c3=[0,],M=[0,],C=[0,],E=[0,],N=[0,];

Variable t,x,y,z;

ODEFunction  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;

Data;

t=[0:1:12];

x=[1.0000,1.0777,1.1443,1.3497,1.5881,1.9814,2.1905,2.5799,2.9763,3.2445,3.5765,3.9315];

y=[1.0000,1.0597,1.0982,1.2035,1.3566,1.4610,1.4771,1.4862,1.6052,1.6094,1.8892,2.0334];

z=[1.0000,1.0920,1.2033,1.4031,1.6920,1.9382,2.2159,2.6186,2.9489,3.1536,3.5979,4.0233];

均方差(RMSE):0.054186898245058
殘差平方和(SSE):0.096895258066869
相關(guān)系數(shù)(R): 0.990468602959854
相關(guān)系數(shù)之平方(R^2): 0.981028053449245
確定系數(shù)(DC): 0.981119298577706
F統(tǒng)計(jì)(F-Statistic): -68.497450670183

參數(shù) 最佳估算
-------------------- -------------
a1 0.934922929556056
a2 43.1779982925144
a3 0.634247035099968
b1 0.885424066820235
b2 1.72323467449167E-14
b3 1.64537787581792
c1 7.17385064510645E-11
c2 47.0215125361918
c3 0.0444107007235367
m 0.4861243914792
c 0.00874309356262975
e 201.823147971463
n 5.99589073558176

====== 結(jié)果輸出 ======

文件: 數(shù)據(jù)文件-1
No t 目標(biāo) x 計(jì)算 x 目標(biāo) y 計(jì)算 y 目標(biāo) z 計(jì)算 z
1 1 1.0777 1.03073467525285 1.0597 1.06691160587989 1.092 1.06234022501214
2 2 1.1443 1.18178677578727 1.0982 1.1387423298674 1.2033 1.22168768708318
3 3 1.3497 1.39065370562373 1.2035 1.21472123201098 1.4031 1.43111557008028
4 4 1.5881 1.63651472709221 1.3566 1.29454185100693 1.692 1.67502754880883
5 5 1.9814 1.91230103994829 1.461 1.37802854578025 1.9382 1.94813144944932
6 6 2.1905 2.21483623814851 1.4771 1.46501620348821 2.2159 2.24783814736362
7 7 2.5799 2.54118886530888 1.4862 1.5553142562366 2.6186 2.57123967133884
8 8 2.9763 2.88710686907347 1.6052 1.64870518739119 2.9489 2.91357690174086
9 9 3.2445 3.2461270702916 1.6094 1.74495815962043 3.1536 3.26712262655401
10 10 3.5765 3.60891777176678 1.8892 1.84385184177429 3.5979 3.62014395661707
11 11 3.9315 3.96272367218441 2.0334 1.9451949237521 4.0233 3.95587454727496

贊一下

回復(fù)此樓

2樓2014-02-20 20:40:16

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

hua99dushu

鐵蟲 (小有名氣)

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

引用回帖:

有幾個(gè)參數(shù)是這樣設(shè)置的：
M[1.5,],C[1,],E[1,],N[1.5,]。
求幫著再跑一下啊。。

贊一下

回復(fù)此樓

3樓2014-02-21 14:54:25

已閱回復(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í)
蟲號: 291104
注冊: 2006-10-28

參考下：

均方差(RMSE):0.0743428646356062
殘差平方和(SSE):0.182386430233526
相關(guān)系數(shù)(R): 0.990490604242565
相關(guān)系數(shù)之平方(R^2): 0.981071637092801
確定系數(shù)(DC): 0.97907392333002
F統(tǒng)計(jì)(F-Statistic): -17.029556181244

參數(shù) 最佳估算
-------------------- -------------
a1 1.30190551784766
a2 4.06003441093189
a3 0.869490127930172
b1 9.42946312355037
b2 1.57560369388824
b3 5.9857438663248
c1 0.0249019284686679
c2 13.4132204832955
c3 9.74396616982886E-6
m 332.144120934994
c 22.8549320352576
e 45695.8913468096
n 1.78908285577578

====== 結(jié)果輸出 ======

文件: 數(shù)據(jù)文件-1
No t 目標(biāo) x 計(jì)算 x 目標(biāo) y 計(jì)算 y 目標(biāo) z 計(jì)算 z
1 1 1.0777 1.13953501723234 1.0597 1.05984607578657 1.092 1.13757057269921
2 2 1.1443 1.30021884916371 1.0982 1.12347532732008 1.2033 1.29904743776794
3 3 1.3497 1.48448399130672 1.2035 1.19141949444282 1.4031 1.48688285600025
4 4 1.5881 1.69476430026936 1.3566 1.26432646601749 1.692 1.70323733352619
5 5 1.9814 1.93339002993677 1.461 1.34297964122961 1.9382 1.94974242279121
6 6 2.1905 2.20245049650347 1.4771 1.42831743928452 2.2159 2.22719556655783
7 7 2.5799 2.50362263043615 1.4862 1.52145197784229 2.6186 2.5351804050362
8 8 2.9763 2.8379668528597 1.6052 1.62368427128187 2.9489 2.87161074307044
9 9 3.2445 3.20569653533687 1.6094 1.73651368452624 3.1536 3.23220275868481
10 10 3.5765 3.60593368625635 1.8892 1.8616395305816 3.5979 3.60988868394394
11 11 3.9315 4.03647149994331 2.0334 2.00094968015341 4.0233 3.99419780445337

贊一下(1人)

回復(fù)此樓

5樓2014-02-24 11:52:03

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

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

[考研] 一志愿西安交大材料學(xué)碩（英一數(shù)二）347，求調(diào)劑到高分子/材料相關(guān)專業(yè) +3	zju51 2026-03-31	5/250	2026-03-31 18:26 by 253863592
[考研] 290分調(diào)劑求助 +10	吉祥止止陳 2026-03-25	10/500	2026-03-31 17:54 by 544594351
[考研] 08開頭看過來�。�！ +3	wwwwffffff 2026-03-31	5/250	2026-03-31 17:45 by 星光/
[考研] 286分調(diào)劑 +11	Faune 2026-03-30	13/650	2026-03-31 17:28 by michael2011
[考研] 一志愿西交大080500材料學(xué)碩349 +5	jqx1258 2026-03-31	6/300	2026-03-31 16:16 by 記事本2026
[考研] 353求調(diào)劑 +3	江上楓_26 2026-03-28	3/150	2026-03-31 15:53 by jp9609
[考研] 08工科，295，接受跨專業(yè)調(diào)劑 +4	lmnlzy 2026-03-31	4/200	2026-03-31 15:25 by 西京學(xué)院招辦
[考研] 266分，求材料冶金能源化工等調(diào)劑 +8	哇呼哼呼哼 2026-03-27	10/500	2026-03-31 13:35 by Huaxue_Wang
[考研] 303求調(diào)劑 +7	DLkz1314. 2026-03-30	7/350	2026-03-30 21:07 by peike
[考研] 310求調(diào)劑 +10	爭取九點(diǎn)睡 2026-03-30	10/500	2026-03-30 16:45 by ztnimte
[考研] 071010 323 分求調(diào)劑 +3	Baekzhy 2026-03-27	3/150	2026-03-30 14:24 by andresqi
[考研] 329求調(diào)劑 +8	星野? 2026-03-26	8/400	2026-03-30 13:41 by chemdavid
[考研] 343求調(diào)劑085601 +3	要努力學(xué)習(xí)x 2026-03-29	3/150	2026-03-29 18:35 by wxiongid
[考研] 329求調(diào)劑 +10	鈕恩雪 2026-03-25	10/500	2026-03-29 13:32 by peike
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6	芝士雪baoo 2026-03-28	8/400	2026-03-29 08:16 by 松花缸1201
[考研] 283求調(diào)劑 +7	A child 2026-03-28	7/350	2026-03-28 12:05 by zllcz
[考研] 340求調(diào)劑 +5	jhx777 2026-03-27	5/250	2026-03-28 04:18 by fmesaito
[考研] 295求調(diào)劑 +5	1428151015 2026-03-27	6/300	2026-03-28 04:04 by fmesaito
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3	askeladz 2026-03-26	3/150	2026-03-26 09:08 by 給你你注意休息