[問題] 聯立微分方程式求數值解
各位好，
因為不曾寫過聯立微分方程式，所以帶不出結果來 Q_Q
還請多指教。
我先是定義一些常數與四則運算式
Kc:= 0.05
Ft:= Fa+Fb+Fc
k:= 0.7
Ct0:= 0.2
ra:=-k*Ct0*(Fa/Ft - Ct0/K* Fb/Ft* Fc/Ft)
kc:= 0.2
接著有三條變數都是V的微分方程式，分別是：
d(Fa)/dV = ra
d(Fb)/dV = -ra - (kc*Ct0*Fb)/Ft
d(Fa)/dV = -ra
個別起始條件是 Fa(V=0) = 10, Fb(V=0) = Fc(V=0) = 0
而我欲想的寫法是：
DSovle[{Fa'[V]==ra, Fb'[V]==-ra-(kc*Ct0*Fb)/Ft, Fc'[V]==-ra,
Fa[0]==10, Fb[0]==0, Fc[0]==0}, {Fa, Fb, Fc}, V]
不過很明顯是跑不出結果 (V=??)
請問我的寫法錯誤在哪呢？
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 112.105.125.179
推
10/04 13:06, , 1F
10/04 13:06, 1F
推
10/04 13:07, , 2F
10/04 13:07, 2F
推
10/04 13:10, , 3F
10/04 13:10, 3F
推
10/04 13:11, , 4F
10/04 13:11, 4F
→
10/12 21:49, , 5F
10/12 21:49, 5F
→
10/15 00:27, , 6F
10/15 00:27, 6F
→
10/15 00:28, , 7F
10/15 00:28, 7F
→
10/15 00:31, , 8F
10/15 00:31, 8F
→
10/15 00:32, , 9F
10/15 00:32, 9F
Mathematica 近期熱門文章
PTT數位生活區 即時熱門文章
