绘制单摆角度与时间曲线的matlab代码

以下是绘制单摆角度与时间曲线的 MATLAB 代码：

% 定义常量
g = 9.81;   % 重力加速度
L = 1;      % 单摆长度
dt = 0.01;  % 时间步长

% 定义初值
theta0 = 0.2;      % 初始角度
omega0 = 0;        % 初始角速度

% 定义时间和状态变量
t = 0:dt:10;       % 时间范围
theta = zeros(size(t));   % 角度
omega = zeros(size(t));   % 角速度

% 初始值
theta(1) = theta0;
omega(1) = omega0;

% 迭代计算
for i = 2:length(t)
    omega(i) = omega(i-1) - (g/L)*theta(i-1)*dt;
    theta(i) = theta(i-1) + omega(i)*dt;
end

% 画图
plot(t, theta);
xlabel('时间');
ylabel('角度');
title('单摆角度与时间曲线');

运行代码后，会得到一个单摆角度与时间曲线的图像。注意，这里只是一个简单的示例，实际上，单摆的运动是一个非线性系统，如果要更准确地描述单摆的运动，需要使用更高级的数学方法和更复杂的模型。