matlab用ode45函数解常系数非齐次线性微分方程组写出范例代码

以下是一个用 Matlab 的 ode45 函数解常系数非齐次线性微分方程组的范例代码：

function y = solve_ode45()
    % 定义初始条件
    y0 = [0; 0];
    % 定义时间步长和终止时间
    tspan = [0 10];
    % 定义常系数非齐次线性微分方程组
    f = @(t, y) [2 -1; 1 -2] * y + [1; 0.5*sin(t)];
    % 用 ode45 函数求解微分方程组
    [~, y] = ode45(f, tspan, y0);
    % 返回结果
    y = y';
end

这个函数解的是如下的常系数非齐次线性微分方程组：

y1' = 2*y1 - y2 + 1
y2' = y1 - 2*y2 + 0.5*sin(t)

其中 y1(0) = 0，y2(0) = 0。