MATLAB 二阶微分方程求解：ode45 函数详解及示例

MATLAB 可以使用 ode45 函数来描述并求解二阶微分方程。以下是一个完整的示例：

假设有一个二阶微分方程：

y'' + 2y' + 2y = sin(t)

其中 y(0) = 0，y'(0) = 0。

可以将其转化为一组一阶微分方程：

z = y'

z' + 2z + 2y = sin(t)

y' = z

然后，可以使用 ode45 函数进行求解：

function dydt = eqns(t,y)
    dydt = zeros(2,1);
    dydt(1) = y(2);
    dydt(2) = -2*y(2) - 2*y(1) + sin(t);
end

[t,y] = ode45(@eqns,[0,10],[0,0]);

plot(t,y(:,1))

运行以上代码，可以得到该二阶微分方程的数值解。

注意，ode45 函数需要提供一个 ODE 方程的句柄（即 @eqns）和初值（即 [0,0]）。同时，该函数的返回值是时间 t 和相应的解 y。在此示例中，我们只需要绘制 y 的值，因此使用 y(:,1) 来表示 y 的值。