用逐步积分法计算单自由度结构在一条地震波作用下的位移响应matlab代码

以下是一个简单的MATLAB代码，使用逐步积分法计算单自由度结构在一条地震波作用下的位移响应：

% 定义常数
m = 1;      % 质量
k = 10;     % 刚度
c = 0.2;    % 阻尼系数
dt = 0.01;  % 时间步长
T = 10;     % 总共计算的时间

% 定义初始条件
u = 0;      % 位移
v = 0;      % 速度
a = 0;      % 加速度

% 定义地震波
g = 9.8;    % 重力加速度
A = 0.1;    % 地震波振幅
w = 2*pi;   % 地震波角频率

% 计算时间步数
t = 0:dt:T;
n = length(t);

% 初始化位移、速度和加速度数组
u_array = zeros(1,n);
v_array = zeros(1,n);
a_array = zeros(1,n);

% 使用逐步积分法计算位移、速度和加速度
for i = 1:n
    % 计算当前时刻的加速度
    a = -c/m*v - k/m*u + A*g*sin(w*t(i))/m;
    
    % 计算当前时刻的速度
    v = v + a*dt;
    
    % 计算当前时刻的位移
    u = u + v*dt;
    
    % 存储当前时刻的位移、速度和加速度
    u_array(i) = u;
    v_array(i) = v;
    a_array(i) = a;
end

% 绘制位移响应曲线
plot(t,u_array);
xlabel('时间 (s)');
ylabel('位移 (m)');
title('单自由度结构在一条地震波作用下的位移响应');

请注意，这只是一个简单的示例代码，实际情况可能会更加复杂。在实际应用中，您可能需要考虑更多的因素，例如地震波的形状、结构的初始条件等。