结构动力学使用中心差分法计算单自由度体系动力反应的matlab程序

以下是使用中心差分法计算单自由度体系动力反应的Matlab程序：

% 定义初始条件
m = 1; %质量
k = 10; %弹性系数
c = 0.1; %阻尼系数
F = 1; %外力
dt = 0.01; %时间步长
t = 0:dt:10; %时间范围

% 定义初始值
x = zeros(size(t)); %位移
v = zeros(size(t)); %速度
a = zeros(size(t)); %加速度

% 中心差分法计算动力反应
for i = 2:length(t)
    x(i) = x(i-1) + v(i-1)*dt + 0.5*a(i-1)*dt^2;
    v(i) = v(i-1) + 0.5*(a(i-1)+a(i))*dt;
    a(i) = (-k*x(i)-c*v(i)+F)/m;
end

% 绘制位移-时间曲线
plot(t,x);
xlabel('时间 (s)');
ylabel('位移 (m)');
title('单自由度体系动力反应');

该程序使用中心差分法计算了单自由度体系的动力反应，并将其绘制为位移-时间曲线。程序中定义了初始条件和初始值，并使用for循环计算了每个时间步长的位移、速度和加速度。最后，程序使用plot函数将位移-时间曲线绘制出来。