2012-11-20 83 views
1

好,我想解決複雜的ODE的形式的系統中:解決複雜的微分方程系統在八度

$ I \ HBAR \ dfrac {\局部\ RHO} {\局部噸} = \左[ H,\ RHO \右] $

http://mathbin.heroku.com/rq65idQ

其中$ \ $ RHO和$ H $是$ n×n的$矩陣。

我試圖定義使用一般功能:

function xdot = f(x,t) 

i*xdot(1)= 
i*xdot(2)= 

endfunction 
x0=[0;0]; 
t=linspace(0,20,200); 
y=lsode("f",x0,t). 

,但我得到了有關一個錯誤消息。那麼我的問題是:如何在八度中求解一個複數微分方程?

回答

0

我認爲這個問題是因爲我是等號的LHS而不是RHS。

你試過了嗎?

xdot(1) = <whatever your RHS expression is>/i; 
xdot(2) = <whatever your RHS expression is>/i;