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

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

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

返回列表

julescesar

鐵桿木蟲 (知名作家)

傻子中的牛虻

應(yīng)助: 85 (初中生)
金幣: 7405.9
散金: 89
紅花: 5
帖子: 5225
在線: 552.7小時(shí)
蟲號(hào): 548890
注冊(cè): 2008-04-19
性別: GG
專業(yè): 環(huán)境地球化學(xué)

[求助] 求高版本1stopt跑一個(gè)程序已有1人參與

程序如下，如果不合語(yǔ)法，請(qǐng)高手幫忙修正后再跑，謝謝啦
Title "lump6";
Parameter k1,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11,k12,k13,k14;
ParameterDomain=[0,];
Variable t,y(6),z(6),p(6),q(6),r(6);
ODEFunction
   y1'=-(k1+k2+k3+k4+k5)*y1/(y1/500+y2/300+y3/150+y4/50+y5/20+y6/400)/(1+50*k6)/(1+0.5*k7);
   y2'=(k1*y1-(k8+k9+k10+k11)*y2)/(y1/500+y2/300+y3/150+y4/50+y5/20+y6/400)/(1+50*k6)/(1+0.5*k7);
   y3'=(k2*y1+k8*y2-(k12+k13+k14)*y3)/(y1/500+y2/300+y3/150+y4/50+y5/20+y6/400)/(1+50*k6)/(1+0.5*k7);
   y4'=(k3*y1+k9*y2+k12*y3)/(y1/500+y2/300+y3/150+y4/50+y5/20+y6/400)/(1+50*k6)/(1+0.5*k7);
   y5'=(k4*y1+k10*y2+k13*y3)/(y1/500+y2/300+y3/150+y4/50+y5/20+y6/400)/(1+50*k6)/(1+0.5*k7);
   y6=100-y1-y2-y3-y4-y5;
   z1'=-(k1+k2+k3+k4+k5)*z1/(z1/500+z2/300+z3/150+z4/50+z5/20+z6/400)/(1+50*k6)/(1+0.4*k7);
   z2'=(k1*z1-(k8+k9+k10+k11)*z2)/(z1/500+z2/300+z3/150+z4/50+z5/20+z6/400)/(1+50*k6)/(1+0.4*k7);
   z3'=(k2*z1+k8*z2-(k12+k13+k14)*z3)/(z1/500+z2/300+z3/150+z4/50+z5/20+z6/400)/(1+50*k6)/(1+0.4*k7);
   z4'=(k3*z1+k9*z2+k12*z3)/(z1/500+z2/300+z3/150+z4/50+z5/20+z6/400)/(1+50*k6)/(1+0.4*k7);
   z5'=(k4*z1+k10*z2+k13*z3)/(z1/500+z2/300+z3/150+z4/50+z5/20+z6/400)/(1+50*k6)/(1+0.4*k7);
   z6=100-z1-z2-z3-z4-z5;
   p1'=-(k1+k2+k3+k4+k5)*p1/(p1/500+p2/300+p3/150+p4/50+p5/20+p6/400)/(1+50*k6)/(1+0.3*k7);
   p2'=(k1*p1-(k8+k9+k10+k11)*p2)/(p1/500+p2/300+p3/150+p4/50+p5/20+p6/400)/(1+50*k6)/(1+0.3*k7);
   p3'=(k2*p1+k8*p2-(k12+k13+k14)*p3)/(p1/500+p2/300+p3/150+p4/50+p5/20+p6/400)/(1+50*k6)/(1+0.3*k7);
   p4'=(k3*p1+k9*p2+k12*p3)/(p1/500+p2/300+p3/150+p4/50+p5/20+p6/400)/(1+50*k6)/(1+0.3*k7);
   p5'=(k4*p1+k10*p2+k13*p3)/(p1/500+p2/300+p3/150+p4/50+p5/20+p6/400)/(1+50*k6)/(1+0.3*k7);
   p6=100-p1-p2-p3-p4-p5;
   q1'=-(k1+k2+k3+k4+k5)*q1/(q1/500+q2/300+q3/150+q4/50+q5/20+q6/400)/(1+50*k6)/(1+0.2*k7);
   q2'=(k1*q1-(k8+k9+k10+k11)*q2)/(q1/500+q2/300+q3/150+q4/50+q5/20+q6/400)/(1+50*k6)/(1+0.2*k7);
   q3'=(k2*q1+k8*q2-(k12+k13+k14)*q3)/(q1/500+q2/300+q3/150+q4/50+q5/20+q6/400)/(1+50*k6)/(1+0.2*k7);
   q4'=(k3*q1+k9*q2+k12*q3)/(q1/500+q2/300+q3/150+q4/50+q5/20+q6/400)/(1+50*k6)/(1+0.2*k7);
   q5'=(k4*q1+k10*q2+k13*q3)/(q1/500+q2/300+q3/150+q4/50+q5/20+q6/400)/(1+50*k6)/(1+0.2*k7);
   q6=100-q1-q2-q3-q4-q5;
   r1'=-(k1+k2+k3+k4+k5)*r1/(r1/500+r2/300+r3/150+r4/50+r5/20+r6/400)/(1+50*k6)/(1+0.1*k7);
   r2'=(k1*r1-(k8+k9+k10+k11)*r2)/(r1/500+r2/300+r3/150+r4/50+r5/20+r6/400)/(1+50*k6)/(1+0.1*k7);
   r3'=(k2*r1+k8*r2-(k12+k13+k14)*r3)/(r1/500+r2/300+r3/150+r4/50+r5/20+r6/400)/(1+50*k6)/(1+0.1*k7);
   r4'=(k3*r1+k9*r2+k12*r3)/(r1/500+r2/300+r3/150+r4/50+r5/20+r6/400)/(1+50*k6)/(1+0.1*k7);
   r5'=(k4*r1+k10*r2+k13*r3)/(r1/500+r2/300+r3/150+r4/50+r5/20+r6/400)/(1+50*k6)/(1+0.1*k7);
   r6=100-r1-r2-r3-r4-r5;
Data;
t y1 y2  y3  y4  y5 y6  z1 z2  z3 z4 z5 z6  p1  p2  p3  p4  p5 p6  q1 q2 q3 q4 q5 q6  r1 r2 r3  r4 r5  r6
0  93.7 6.3  0 0 0 0  93.7 6.3 0 0  0 0  93.7  6.3  0 0  0 0 93.7 6.3  0  0    0 0 93.7 6.3  0 0  0  0
1 5.34 16  47  20 2  9.66  4.95 15 48 19  3  10.05 4.42 13 48  20  3  11.58 3.79  13 47  21 3.3 11.91  3.58 12 47 22  3.3 12.12

回復(fù)此樓

» 猜你喜歡

請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有139人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有76人回復(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日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)

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

1stopt解方程組已經(jīng)有12人回復(fù)
用1stopt軟件做微分方程組的參數(shù)估計(jì) 已經(jīng)有13人回復(fù)
求用高版本的1stOpt擬合方程已經(jīng)有8人回復(fù)
有1stopt的蟲友給跑下程序，求！已經(jīng)有5人回復(fù)
我的MS低版本打不開高版本，怎么辦啊，急求助！�。� 已經(jīng)有9人回復(fù)
誰(shuí)用1stopt這個(gè)軟件啊，離散數(shù)據(jù)擬合方程，方程未知已經(jīng)有5人回復(fù)
求大神幫忙用1stOpt 擬合復(fù)數(shù)數(shù)據(jù) 已經(jīng)有13人回復(fù)
rmse，sse計(jì)算求助已經(jīng)有3人回復(fù)
求高版本1stopt擬合，已經(jīng)有11人回復(fù)
1stopt中使用sharedmodel擬合曲線出現(xiàn)問題，求指示�。� 已經(jīng)有6人回復(fù)
微分方程組參數(shù)擬合的問題（多參數(shù)動(dòng)力系統(tǒng)參數(shù)估計(jì)）已經(jīng)有11人回復(fù)
求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有9人回復(fù)
求大俠幫助跑一下1stopt 已經(jīng)有12人回復(fù)
誰(shuí)有1stopt，幫忙跑一下啊已經(jīng)有6人回復(fù)
哪位能用1stopt的高版本幫我走個(gè)模型啊！謝謝了！我的1.5的不能跑兩模型共參數(shù)啊已經(jīng)有3人回復(fù)
哪位大哥能幫我用1stopt3.0以上版本跑一個(gè)擬合程序已經(jīng)有11人回復(fù)

一個(gè)罪犯也有權(quán)檢舉另一個(gè)罪犯，以事實(shí)為采信依據(jù)；反之同理。神州法治move-forward！

1樓 2014-04-08 09:03:48

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

julescesar

鐵桿木蟲 (知名作家)

傻子中的牛虻

應(yīng)助: 85 (初中生)
金幣: 7405.9
散金: 89
紅花: 5
帖子: 5225
在線: 552.7小時(shí)
蟲號(hào): 548890
注冊(cè): 2008-04-19
性別: GG
專業(yè): 環(huán)境地球化學(xué)

請(qǐng)有資源的朋友幫個(gè)忙，謝謝

贊一下

回復(fù)此樓

一個(gè)罪犯也有權(quán)檢舉另一個(gè)罪犯，以事實(shí)為采信依據(jù)；反之同理。神州法治move-forward！

2樓2014-04-08 11:21:39

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

julescesar

鐵桿木蟲 (知名作家)

傻子中的牛虻

應(yīng)助: 85 (初中生)
金幣: 7405.9
散金: 89
紅花: 5
帖子: 5225
在線: 552.7小時(shí)
蟲號(hào): 548890
注冊(cè): 2008-04-19
性別: GG
專業(yè): 環(huán)境地球化學(xué)

求腳印。。。。

回復(fù)此樓

一個(gè)罪犯也有權(quán)檢舉另一個(gè)罪犯，以事實(shí)為采信依據(jù)；反之同理。神州法治move-forward！

3樓2014-04-08 12:24:12

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

julescesar

鐵桿木蟲 (知名作家)

4樓2014-04-08 15:09:13

已閱回復(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

將就參考下：

參數(shù) 最佳估算
-------------------- -------------
k1 2.0920685752888
k2 2.03145877517177E-9
k3 0.248632668721487
k4 1.35390776847532E-17
k5 6.03374135447111E-9
k6 3.61924495624649E-9
k7 0.689252026420965
k8 2.16983628672973
k9 2.98161186429235E-10
k10 4.39628051135967E-8
k11 1.57735517991293E-9
k12 0.278468868858522
k13 0.0748917320017696
k14 0.283316627567864

贊一下

回復(fù)此樓

5樓2014-04-08 15:30:28

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

julescesar

鐵桿木蟲 (知名作家)

傻子中的牛虻

應(yīng)助: 85 (初中生)
金幣: 7405.9
散金: 89
紅花: 5
帖子: 5225
在線: 552.7小時(shí)
蟲號(hào): 548890
注冊(cè): 2008-04-19
性別: GG
專業(yè): 環(huán)境地球化學(xué)

引用回帖:

5樓: Originally posted by dingd at 2014-04-08 15:30:28
將就參考下：

參數(shù) 最佳估算
-------------------- -------------
k1 2.0920685752888
k2 2.03145877517177E-9
k3 0.248632668721487
k4 1.35390776847532E-17
k5 6.03374135447111E-9
k6 3.61 ...

謝謝版主。。。
能不能把你修正過的程序代碼和計(jì)算后的實(shí)測(cè)值與模擬值比較部分的數(shù)據(jù)貼出來(lái)？謝謝啦

贊一下

回復(fù)此樓

一個(gè)罪犯也有權(quán)檢舉另一個(gè)罪犯，以事實(shí)為采信依據(jù)；反之同理。神州法治move-forward！

6樓2014-04-08 15:52:30

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

julescesar

鐵桿木蟲 (知名作家)

傻子中的牛虻

應(yīng)助: 85 (初中生)
金幣: 7405.9
散金: 89
紅花: 5
帖子: 5225
在線: 552.7小時(shí)
蟲號(hào): 548890
注冊(cè): 2008-04-19
性別: GG
專業(yè): 環(huán)境地球化學(xué)

引用回帖:

k2的數(shù)值太小了點(diǎn)兒，感覺不大對(duì)

贊一下

回復(fù)此樓

一個(gè)罪犯也有權(quán)檢舉另一個(gè)罪犯，以事實(shí)為采信依據(jù)；反之同理。神州法治move-forward！

7樓2014-04-08 16:15:22

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

julescesar

鐵桿木蟲 (知名作家)

傻子中的牛虻

應(yīng)助: 85 (初中生)
金幣: 7405.9
散金: 89
紅花: 5
帖子: 5225
在線: 552.7小時(shí)
蟲號(hào): 548890
注冊(cè): 2008-04-19
性別: GG
專業(yè): 環(huán)境地球化學(xué)

引用回帖:

請(qǐng)版主指點(diǎn)一二，小蟲外行

贊一下

回復(fù)此樓

一個(gè)罪犯也有權(quán)檢舉另一個(gè)罪犯，以事實(shí)為采信依據(jù)；反之同理。神州法治move-forward！

8樓2014-04-08 21:55:34

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

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

返回列表

普通表情龍兔虎貓

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

[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +10	手機(jī)用戶 2026-03-17	11/550	2026-03-23 16:30 by lingjue
[考研] 306求0703調(diào)劑一志愿華中師范 +6	紙魚ly 2026-03-21	6/300	2026-03-23 12:08 by Iveryant
[考研] 333求調(diào)劑 +6	87639 2026-03-21	10/500	2026-03-23 10:41 by Iveryant
[考研] 276求調(diào)劑 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在風(fēng)里
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚.. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11	楓橋ZL 2026-03-18	13/650	2026-03-22 20:26 by edmund7
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3	這么名字咋樣 2026-03-22	4/200	2026-03-22 17:56 by 云民大李老師
[考研] 308求調(diào)劑 +3	墨墨漠 2026-03-21	3/150	2026-03-22 16:54 by i_cooler
[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 0805材料320求調(diào)劑 +3	深海物語(yǔ) 2026-03-20	3/150	2026-03-21 15:46 by 無(wú)際的草原
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 279分求調(diào)劑一志愿211 +14	chaojifeixia 2026-03-19	15/750	2026-03-21 13:24 by zhukairuo
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 301求調(diào)劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky

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

[求助] 求高版本1stopt跑一個(gè)程序 已有1人參與

» 猜你喜歡

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

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

[求助] 求高版本1stopt跑一個(gè)程序已有1人參與

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