| 2 | 1/1 | 返回列表 |
| 查看: 277 | 回復(fù): 1 | |||
martinwgc金蟲 (小有名氣)
|
[求助]
求助:這個(gè)matlab的代碼有問題嗎? 已有1人參與
|
|
我這里有一個(gè)偏微分方程組。希望得到它在t=0,3,6,12,24,48,72,120,168的時(shí)候的各個(gè)變量的值。為啥我這個(gè)運(yùn)行出來屏幕上啥也沒有?請幫我看看代碼,謝謝了! function Model1 clear;clc x0=[0,0,0,100,0.4368,1.1432,0.5389,0.0511,2156.08,946.60,0,0]; [t,x]=ode45(@sac1,[0,3,6,12,24,48,72,120,168],x0); x; function dx=sac1(t,x); global E1max; global K1ad; global E1T; global E2T; global S0; global K2ad; global E2max; global K1r; global K1IG2; global K1IG; global K2r; global K2IG2; global K2IG; global K3r; global K3M; global K3IG; global KG3plus; global KG3minus; global Gcrtri; E1max=0.03257; K1ad=1.238; E1T=1.58; E2T=0.59; S0=100; K2ad=1.865; E2max=0.00102; K1r=18.86; K1IG2=0.0042; K1IG=0.0947; K2r=7.926; K2IG2=147.87; K2IG=0.01509; K3r=228.264; K3M=3.1740; K3IG=0.9617; KG3plus=0.0075; KG3minus=0.1197; Gcrtri=40; dx=zeros(12,1); x(6)=E1max*K1ad*x(5)*x(4)/(1+K1ad*x(5)); x(8)=E2max*K2ad*x(7)*x(4)/(1+K2ad*x(7)); x(9)=K1r*x(6)*x(4)^2/(S0*(1+x(2)/K1IG2+x(1)/K1IG)); x(10)=K2r*(x(6)+x(8))*x(4)^2/(S0*(1+x(2)/K2IG2+x(1)/K2IG)); x(11)=K3r*x(7)*x(2)/(K3M*(1+x(1)/K3IG)+x(2)); dx(4)=-x(9)-x(10); dx(2)=1.056*x(9)-x(11); dx(1)=1.111*x(10)+1.053*x(11)-1.071*x(12); dx(3)=x(12); x(5)=E1T-x(6); x(7)=E2T-x(8); x(12)=(KG3plus*x(1)-KG3minus*x(3))*(1/(1+exp(Gcrtri-x(1)))); |
銅蟲 (初入文壇)
|
把第一個(gè)函數(shù)里邊x后面的分號(hào)去掉。以下是顯示的結(jié)果 x = 1.0e+03 * 0 0 0 0.1000 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0178 0.0013 0 0.0828 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0255 0.0013 0 0.0759 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0354 0.0011 0 0.0671 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0470 0.0008 0 0.0570 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0591 0.0005 0 0.0464 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0658 0.0004 0 0.0405 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0736 0.0002 0 0.0336 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 0.0782 0.0001 0 0.0295 0.0004 0.0011 0.0005 0.0001 2.1561 0.9466 0 0 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +6 | 十三加油 2026-03-21 | 6/300 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +3 | vv迷 2026-03-22 | 3/150 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 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 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|