[问题] 联立微分方程式求数值解

楼主: 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