[問題] 聯立微分方程式求數值解

作者: theye (Momiji)   2013-10-03 18:38:45
各位好,
因為不曾寫過聯立微分方程式,所以帶不出結果來 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=??)
請問我的寫法錯誤在哪呢?
作者: Frobenius (▽.(▽×▽φ)=0)   0000-00-00 00:00:00
Ft是Ft:=Fa+Fb+Fc=10+0+0=10還是Ft:=Fa[V]+Fb[V]+Fc[V]?k跟K一樣嗎?其他式中的Fa,Fb,Fc是Fa[V],Fb[V],Fc[V]還是常數?DSolve[...,{Fa, Fb, Fc},V]應改為Fa[V],Fb[V],Fc[V]
作者: theye (Momiji)   0000-00-00 00:00:00
謝謝費伯尼斯!藉費伯尼斯大大的建議,我的方程式已是可運算的狀態但這幾天還是試不出:給定最終條件,其他函數的值如給定Vf:=500,找 Fa[Vf]=??看課本都寫得很簡單的樣子 哈哈

Links booklink

Contact Us: admin [ a t ] ucptt.com