举个matlab求解差分方程的例子

在MATLAB中，你可以使用`ode45`函数来求解常微分方程（ODE），它是一个四阶龙格-库塔方法。这里给出一个简单的例子，假设我们要解一阶线性常微分方程 dy/dt = f(t, y)，其中f是一个描述系统动态的函数。

% 定义微分方程的右侧函数
function dydt = my_diffeq(t, y)
    % 这里y是状态向量，例如 y(1)代表y值，y(2)代表dy/dt等
    dydt = [y(2); -0.1 * y(1)]; % 假设是个衰减振荡的例子，d^2y/dt^2 = -0.1*y

% 初始条件
y0 = [1; 1]; % 假设初始位置和速度为1

% 时间范围
tspan = [0 10]; % 求解从0到10秒的时间区间

% 解方程
[t, y] = ode45(@my_diffeq, tspan, y0);

% 绘制结果
plot(t, y(:, 1)); % 只画出y值随时间的变化
xlabel('Time');
ylabel('State variable');
title('Solution of the differential equation');

%