 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 752 | 回復(fù): 1 | |||
[交流]
【求助】解剛性微分方程組 已有1人參與
|
|
我要解一個(gè)44維剛性微分方程組，使用ode23s，但計(jì)算結(jié)果明顯不對，在此求助大家！微分方程組是關(guān)于化學(xué)反應(yīng)動力學(xué)方面的 代碼如下： A=zeros(5,1); E=zeros(5,1); tspan=[0 2820]; y0=zeros(1,44); y0(1,44)=0.003425; [t,y]=ode23s('rateequation', tspan, y0,[]); function ydot=rateequation(t,y) A1=3.98E+12; E1=281.4E+3; A2=1.58E+10; E2=49.14E+3; A3=1.00E+14; E3=130.2E+3; A4=3.16E+7; E4=57.54E+3; A5=1.00E+10; E5=0; T=653.15; R=8.314; %J/mol/K k1=A1*exp(E1/R/T); %求解動力學(xué)速率常數(shù) k2=A2*exp(E2/R/T); k3=A3*exp(E3/R/T); k4=A4*exp(E4/R/T); k5=A5*exp(E5/R/T); ydot=zeros(44,1); ydot(1)=k3*y(42); %所求微分方程組，44維，y（）為某一組分的濃度 ydot(2)=k3*y(43); ydot(3)=k3*y(43); ydot(4)=k3*y(42); ydot(5)=k3*y(41); ydot(6)=k3*y(40); ydot(7)=k3*y(39); ydot(8)=28*k4*y(44)*y(24)+2*k5*(y(15)*y(19)+2*y(16)*y(18)+2*y(17)*y(17)); ydot(9)=14*k4*y(44)*y(26)+2*k5*(y(15)*y(21)+2*y(16)*y(19)+2*y(17)*y(18)); ydot(10)=28*k4*y(44)*y(28)+2*k5*(y(15)*y(23)+2*y(16)*y(21)+2*y(17)*y(19)+2*y(18)*y(18)); ydot(11)=28*k4*y(44)*y(30)+2*k5*(y(15)*y(25)+2*y(16)*y(23)+2*y(17)*y(21)+2*y(18)*y(19)); ydot(12)=28*k4*y(44)*y(32)+2*k5*(y(15)*y(27)+2*y(16)*y(25)+2*y(17)*y(23)+2*y(18)*y(21)+2*y(19)*y(19)); ydot(13)=28*k4*y(44)*y(34)+2*k5*(y(15)*y(29)+2*y(16)*y(27)+2*y(17)*y(25)+2*y(18)*y(23)+2*y(19)*y(21)); ydot(14)=28*k4*y(44)*y(36)+2*k5*(y(15)*y(31)+2*y(16)*y(29)+2*y(17)*y(27)+2*y(18)*y(25)+2*y(19)*y(23)+2*y(21)*y(21)); ydot(15)=k1*y(44)-2*k5*y(15)*(y(19)+y(21)+y(23)+y(25)+y(27)+y(29)+y(31))+k3*y(39); ydot(16)=k1*y(44)-4*k5*y(16)*(y(18)+y(19)+y(21)+y(23)+y(25)+y(27)+y(29))+k3*y(37); ydot(17)=k1*y(44)-4*k5*y(17)*(y(17)+y(18)+y(19)+y(21)+y(23)+y(25)+y(27))+k3*y(41); ydot(18)=k1*y(44)-4*k5*y(18)*(y(16)+y(17)+y(18)+y(19)+y(21)+y(23)+y(25))+k3*y(42); ydot(19)=k1*y(44)-2*k5*y(19)*(y(15)+2*y(16)+2*y(17)+2*y(18)+2*y(19)+2*y(21)+2*y(23))+k3*y(43)-4*k2*y(19); ydot(20)=4*k2*y(19)-28*k4*y(20)*y(44); ydot(21)=k1*y(44)-2*k5*y(21)*(y(15)+2*y(16)+2*y(17)+2*y(18)+2*y(19)+2*y(21))+k3*y(43)-4*k2*y(21); ydot(22)=4*k2*y(21)-28*k4*y(22)*y(44); ydot(23)=k1*y(44)-2*k5*y(23)*(y(15)+2*y(16)+2*y(17)+2*y(18)+2*y(19))+k3*y(42)-4*k2*y(23); ydot(24)=4*k2*y(23)-28*k4*y(24)*y(44); ydot(25)=k1*y(44)-2*k5*y(25)*(y(15)+2*y(16)+2*y(17)+2*y(18))+k3*y(41)-4*k2*y(25); ydot(26)=4*k2*y(25)-28*k4*y(26)*y(44); ydot(27)=k1*y(44)-2*k5*y(27)*(y(15)+2*y(16)+2*y(17))+k3*y(40)-4*k2*y(27); ydot(28)=4*k2*y(27)-28*k4*y(28)*y(44); ydot(29)=k1*y(44)-2*k5*y(29)*(y(15)+2*y(16))+k3*y(39)-4*k2*y(29); ydot(30)=4*k2*y(29)-28*k4*y(30)*y(44); ydot(31)=k1*y(44)-2*k5*y(15)*y(31)+k3*y(38)-4*k2*y(31); ydot(32)=4*k2*y(31)-28*k4*y(32)*y(44); ydot(33)=k1*y(44)+k3*y(37)-4*k2*y(33); ydot(34)=4*k2*y(33)-28*k4*y(34)*y(44); ydot(35)=k1*y(44)-4*k2*y(35); ydot(36)=4*k2*y(35)-28*k4*y(36)*y(44); ydot(37)=2*y(44)*k4*(2*y(24)+y(26)+2*y(28)+2*y(30)+2*y(32)+2*y(34)+2*y(36))-k3*y(37); ydot(38)=2*y(44)*k4*(2*y(24)+y(26)+2*y(28)+2*y(30)+2*y(32)+2*y(34)+2*y(36))-k3*y(38); ydot(39)=2*y(44)*k4*(2*y(24)+y(26)+2*y(28)+2*y(30)+2*y(32)+2*y(34)+2*y(36))-k3*y(39); ydot(40)=2*y(44)*k4*(2*y(24)+y(26)+2*y(28)+2*y(30)+2*y(32)+2*y(34)+2*y(36))-k3*y(40); ydot(41)=2*y(44)*k4*(2*y(24)+y(26)+2*y(28)+2*y(30)+2*y(32)+2*y(34)+2*y(36))-k3*y(41); ydot(42)=2*y(44)*k4*(2*y(24)+y(26)+2*y(28)+2*y(30)+2*y(32)+2*y(34)+2*y(36))-k3*y(42); ydot(43)=2*y(44)*k4*(2*y(24)+y(26)+2*y(28)+2*y(30)+2*y(32)+2*y(34)+2*y(36))-k3*y(43); ydot(44)=-7*k1*y(44)-14*y(44)*k4*(2*y(24)+y(26)+2*y(28)+2*y(30)+2*y(32)+2*y(34)+2*y(36)); |
回復(fù)此樓» 猜你喜歡
070300化學(xué)求調(diào)劑
已經(jīng)有3人回復(fù)
-
311求調(diào)劑
已經(jīng)有3人回復(fù)
-
生物學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
269專碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有8人回復(fù)
-
尋找調(diào)劑
已經(jīng)有3人回復(fù)
-
291 求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有6人回復(fù)
-
318求調(diào)劑
已經(jīng)有7人回復(fù)
信彼南山
木蟲 (著名寫手)
- 應(yīng)助: 33 (小學(xué)生)
- 金幣: 4142.9
- 散金: 1221
- 紅花: 16
- 帖子: 1178
- 在線: 233.5小時(shí)
- 蟲號: 1133529
- 注冊: 2010-10-27
- 專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)
2樓2011-04-07 00:32:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 nadou 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)求調(diào)劑 +3 | 苑豆豆 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
|
[考研] 尋找調(diào)劑 +3 | 倔強(qiáng)芒? 2026-03-21 | 3/150 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+3
|
酥酥魚.. 2026-03-21 | 3/150 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
|
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3 | 1911236844 2026-03-17 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
|
[論文投稿]
申請回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|