[問題] 矩陣+微分方程+ode45的問題
我想要解一個微分方程如下
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
→
12/17 13:06, , 1F
12/17 13:06, 1F