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

登錄注冊

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

24小時(shí)熱門版塊排行榜

返回列表

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é)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有155人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有78人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(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的回帖

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é)方法

引用回帖:

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

有幾個(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的回帖

hua99dushu

鐵蟲 (小有名氣)

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

引用回帖:

之前程序里參數(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的回帖

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

hua99dushu

鐵蟲 (小有名氣)

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

引用回帖:

5樓: Originally posted by dingd at 2014-02-24 11:52:03
參考下：

均方差(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-Statist ...

真的很感謝你如此迅速的回復(fù)，多次認(rèn)真的解答我的問題。
根據(jù)實(shí)際的條件，已經(jīng)確定了一些參數(shù)。
[/CODE]
Parameters a1=[0,],a2=[0,],a3=[0,],b1=[0,],b2=[0,],b3=[0,],c1=[0,],c2=[0,],c3=[0,];
Variable t,x,y,z;
ODEFunction x'=a1*x*(y/33.214-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/22.854)+b3*z*(1-z/45.695);
z'=c1*x*(x/1.789-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];
求幫忙再跑一下上面的代碼，還有就是是采用的什么算法�。�

贊一下

回復(fù)此樓

6樓2014-02-24 14:46:18

已閱回復(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é)方法

有一行錯(cuò)了，正確的是：ODEFunction x'=a1*x*(y/13.214-1)-a2*y+a3*z;
y/M，m=13.214.

贊一下

回復(fù)此樓

7樓2014-02-24 15:15:05

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

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

返回列表

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

[考研] 293求調(diào)劑 +4	加一一九 2026-03-24	4/200	2026-03-24 11:38 by 星空星月
[考研] 一志愿211 初試270分求調(diào)劑 +4	谷雨上岸 2026-03-23	5/250	2026-03-23 21:18 by 不惑可樂
[考研] 303求調(diào)劑 +4	元夕元 2026-03-20	4/200	2026-03-23 19:00 by macy2011
[考研] 350求調(diào)劑 +6	weudhdk 2026-03-19	6/300	2026-03-23 15:47 by tangyuan0840221
[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-20	3/150	2026-03-22 16:00 by ColorlessPI
[考研] 生物學(xué)調(diào)劑 +5	Surekei 2026-03-21	5/250	2026-03-22 14:39 by tcx007
[考研] 269專碩求調(diào)劑 +6	金恩貝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[考研] 302求調(diào)劑 +12	呼呼呼。。。。 2026-03-17	12/600	2026-03-21 17:29 by ColorlessPI
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-18	6/300	2026-03-21 00:32 by JourneyLucky
[考研] 一志愿西南交大，求調(diào)劑 +5	材化逐夢人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne

24小時(shí)熱門版塊排行榜

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

» 猜你喜歡

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

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

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

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