Question

matlab用ode45求解多元常微分方程组

Dx = 2*x-3*y+3*z
Dy= 4*x - 5*y + 3*z
Dz= 4*x + 4*y + 2*z
应该怎样编写matlab的ode方程并求解呢？

Answer 1

题主给出的多元常微分方程组缺初始条件，为了解题，我们自行给出（到是题主修改一下就可以了）。

x0=[0.01,0.01,0.01]; %初值

func=@(t,x)[2*x(1)-3*x(2)+3*x(3);4*x(1)-5*x(2)+3*x(3);4*x(1)+4*x(2)+2*x(3)];

[t,x]=ode45(func,[0 2],x0);

disp('          t         x(t)         y(t)         z(t)')

disp([t x])