| 2 | 1/1 | 返回列表 |
| 查看: 276 | 回復(fù): 1 | ||
martinwgc金蟲 (小有名氣)
|
[求助]
求助:這個(gè)matlab的代碼有問題嗎? 已有1人參與
|
|
我這里有一個(gè)偏微分方程組。希望得到它在t=0,3,6,12,24,48,72,120,168的時(shí)候的各個(gè)變量的值。為啥我這個(gè)運(yùn)行出來屏幕上啥也沒有?請(qǐng)幫我看看代碼,謝謝了! 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ā)表 | |
|---|---|---|---|---|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|