Question

怎么用matlAB解二阶常微分方程组

现在有四个方程，二阶的。怎么用龙格库塔法解，谢谢

Answer 1

必须先转换成一阶微分方程组，才可以解。

追问

怎么同时解4个方程呢

追答

关键不在几个，关键在必须是一阶微分方程。
举个例子给你：
function test_solvede[T, Y] = ode45(@dmyfun, [0 50], [0.02 0.98]);plot(T, Y);legend f gfunction dy = dmyfun(t, y)dy = zeros(2, 1);dy(1) = y(1)*(y(2)-0.3);dy(2) = -y(1)*y(2);
这是两个方程构成的一阶微分方程组。

关键不在几个，关键在必须是一阶微分方程。
举个例子给你：
function test_solvede
[T, Y] = ode45(@dmyfun, [0 50], [0.02 0.98]);plot(T, Y);
legend f g

function dy = dmyfun(t, y)
dy = zeros(2, 1);
dy(1) = y(1)*(y(2)-0.3);
dy(2) = -y(1)*y(2);

这是两个方程构成的一阶微分方程组。