作者shoduit (比一千個太陽還亮)
看板MATLAB
標題[問題] ode45解方程式
時間Thu Jan 12 19:59:01 2012
最近遇到一個問題 總是覺得很怪
所以想請教大家~
問題是這樣的:
我用ode45解一條方程式
大致上是 u_z = - 1i*D*u - 1i*2*|U0|^2 u
D 是一個矩陣, U0 = U0(z,t)
問題就在於這個U0.
當我給U0一個解析的函數(寫個function的副程式,讓電腦去call U0)
如 U0(z,t) = sech(t).*exp(1i*z/2);
算出來的答案是正確的.
但是當我把U0化成矩陣代到ode45去解的時候
得出的答案就不一樣(錯的答案)
請問如果U0只能得到矩陣, 無法寫成解析的函數時
該如何去修正程式 才能得到正確的解?
麻煩大家幫幫忙....多謝 ~
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.114.213.62