我的研究碰到此下列微分方程式：
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