帮我写一个matlab代码，研究弹簧摆的运动，要求做出弹簧摆的模拟运动图像并画出摆锤的运动轨迹

以下是一个基本的弹簧摆的Matlab代码，可以绘制出弹簧摆的模拟运动图像并画出摆锤的运动轨迹：

% 声明初始值
g = 9.8; % 重力加速度
L = 1; % 摆长
k = 1; % 弹簧系数
m = 1; % 摆锤质量
theta0 = pi/4; % 初始摆角
v0 = 0; % 初始角速度
tspan = [0 10]; % 时间范围

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

% 解微分方程
[t, y] = ode45(f, tspan, [theta0; v0]);

% 计算摆锤位置
x = L*sin(y(:,1));
y = -L*cos(y(:,1));

% 绘制模拟运动图像
figure;
for i = 1:length(t)
    plot([0 x(i)], [0 y(i)], 'r-', 'LineWidth', 2);
    hold on;
    plot(x(i), y(i), 'ro', 'MarkerSize', 10, 'MarkerFaceColor', 'r');
    xlim([-L-0.5 L+0.5]);
    ylim([-L-0.5 0.5]);
    xlabel('x');
    ylabel('y');
    title('弹簧摆模拟运动');
    pause(0.1);
    if i < length(t)
        clf;
    end
end

% 绘制摆锤运动轨迹
figure;
plot(x, y, 'r-', 'LineWidth', 2);
hold on;
plot(x(1), y(1), 'ro', 'MarkerSize', 10, 'MarkerFaceColor', 'r');
plot(x(end), y(end), 'r*', 'MarkerSize', 10, 'MarkerFaceColor', 'r');
xlim([-L-0.5 L+0.5]);
ylim([-L-0.5 0.5]);
xlabel('x');
ylabel('y');
title('弹簧摆运动轨迹');

运行该代码可以得到弹簧摆的模拟运动图像和摆锤的运动轨迹。可以根据需要调整初始值和时间范围等参数进行模拟。