作者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