matlab生成重物质量m不同时的响应曲线（L=1m），满足一部分等幅振荡一个周期趋于稳定，另一部分斜坡后趋于稳定的代码

以下是生成重物质量m不同时的响应曲线的MATLAB代码：

clear all;
close all;

L = 1; % 长度
tmax = 20; % 时间长度
dt = 0.01; % 时间步长
t = 0:dt:tmax; % 时间
n = length(t); % 时间步数

% 重物质量m不同时的响应曲线
figure; hold on;
for m = [1, 5, 10] % 重物质量
    x = zeros(n, 1); % 位移
    v = zeros(n, 1); % 速度
    a = zeros(n, 1); % 加速度
    f = zeros(n, 1); % 力
    w0 = sqrt(9.81/L); % 自然角频率
    f0 = 1; % 力的振幅
    f(1) = f0; % 初值
    for i = 1:n-1
        a(i+1) = (f(i) - m*w0^2*x(i))/m; % 计算加速度
        v(i+1) = v(i) + a(i+1)*dt; % 计算速度
        x(i+1) = x(i) + v(i+1)*dt; % 计算位移
        if x(i+1) > 0.5 % 斜坡后趋于稳定
            f(i+1) = 0; % 施加外力为0
        else % 一部分等幅振荡一个周期趋于稳定
            f(i+1) = f0*sin(2*pi*f0*t(i+1));
        end
    end
    plot(t, x, 'LineWidth', 1.5);
end
xlabel('时间');
ylabel('位移');
title('重物质量m不同时的响应曲线');
legend('m=1', 'm=5', 'm=10');

其中，当位移超过0.5时，施加外力为0，即斜坡后趋于稳定；当位移小于等于0.5时，施加外力为等幅振荡的正弦函数，即一部分等幅振荡一个周期趋于稳定。可以通过修改重物质量m、力的振幅f0等参数来生成不同的响应曲线。