帮我写一个matlab单摆的简谐振动代码,要求结果有图像

下面是一个简单的 Matlab 单摆的简谐振动代码，其中使用了欧拉法求解微分方程。在代码中，可以通过修改初始条件、时间步长和模拟时间等参数，观察单摆的运动情况，最终得到单摆振动的图像。

g = 9.81; % 重力加速度
L = 1; % 单摆长度
theta0 = 0.2; % 初始角度
omega0 = 0; % 初始角速度

% 定义微分方程
f = @(t, x) [x(2); -g/L*sin(x(1))];

% 模拟参数
dt = 0.01; % 时间步长
tmax = 10; % 模拟时间

% 初始化
t = 0:dt:tmax;
n = length(t);
y = zeros(n, 2);
y(1, :) = [theta0, omega0];

% 使用欧拉法求解微分方程
for i = 2:n
    y(i, :) = y(i-1, :) + dt*f(t(i-1), y(i-1, :));
end

% 绘制图像
figure;
plot(t, y(:, 1));
xlabel('时间');
ylabel('角度');
title('单摆的简谐振动');

运行这段代码，即可得到一个简单的单摆简谐振动的图像，如下所示：


在这个图像中，我们可以看到单摆的角度随时间的变化而周期性地振荡，符合简谐振动的特点。如果需要观察更多的运动情况，可以尝试修改初始条件、时间步长和模拟时间等参数，以得到不同的结果。