如何在Matlab中绘制弹性弹簧摆的摆动

% 定义模型参数
initial_length = 1; % 弹簧的初始长度
spring_constant = 10; % 弹簧的劲度系数
mass = 1; % 质点的质量
gravity = 9.8; % 重力加速度

% 定义微分方程
dydt = @(t, y) [y(2); -gravity - (spring_constant/mass) * (y(1) - initial_length)];

% 设置时间步长和模拟时间
dt = 0.01; % 时间步长
tspan = 0:dt:10; % 模拟时间范围

% 求解微分方程
[t, y] = ode45(dydt, tspan, [initial_length; 0]);

% 绘制摆动曲线
plot(t, y(:, 1));
xlabel('时间');
ylabel('位置');
title('弹性弹簧摆的摆动');