作者powerDNA (Mayday)
看板MATLAB
标题[运算] 用ode45解联立常微分方程
时间Tue May 28 19:30:30 2013
我的研究碰到此下列微分方程式：
d Ka(t)/dt=Ka(t)[ie(t)-ra0]+ga(0)v(t)+Sum[Caa'm(t),{a'=L,R},{m,1,100}]-Sum[
Daa'(t),{a'=L,R}]; (1)
d Qam(t)/dt=Qam[ie(t)-ram]-nam+Sum[dagger[Ca'am(t)],{a'=L,R}]; (2)
d Caa'm(t)/dt=-ga(0)dagger[Qa'm(])]-[ra0+ram']Caa'm(t); (3)
d Daa'(t)=ga(0)dagger[Ka'(t)]+Ka(t)ga'(0)-[ra0+ra'0]Daa'(t); (4)
这里的m=1~100,a,a'=Left or Right. 对某一组特定的m, 先解(2),(4), 解出来的
Caa'm(t), 然後这100组的Caa'm(t)全部求出来. 再把它带入(1)作求和Sum的计算,
再把(1),(4)一起解. 我是用ode45, 请问要怎麽把计算出来numerical的Caa'm(t)求和算
出来output结果变成(1)式的input function求Ka(t)和Daa'(t)？因为ode45好像只能接受
除了被解的函数以外其他在微分方程里面的项必须写进去必须是exact function, 或者
其他人有更好地解决方式也可以. 谢谢！
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.25.105