作者narwhalgary (独角鲸葛瑞)
看板MATLAB
标题[问题] 矩阵+微分方程+ode45的问题
时间Sat Dec 17 01:27:53 2011
我想要解一个微分方程如下
G和L都是complex的矩阵 L是20*20
dG(t)/dt = L(t)*G(t)
一开始是用G(t)=exp(L*t)的方法去叠代解出G
但是学长说用ode45应该可以更快解出来
问题如下
1. 解出来後发现G是4*4为什麽?
2. 如何直接用ode45去解这个题目??
学长说ode45只能处理vector
所以他做了这样的事情
function dg = myfun(...)
G = reshape(g,20,20);
G = L*G;
dg = G(:);
end
===========================
G_ini=ones(1,400)
ode45(@myfun,[0,10],G_ini)
因为不管怎麽做都怪怪的
请问有错或是有更好的算法吗?
关於矩阵的微分方程??
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.247.201
1F：→ Rasin:你需要的是Mathematica 12/17 13:06