matlab单自由度振动响应源代码

以下是一个简单的MATLAB单自由度振动响应的源代码示例：

% 定义系统参数
m = 1; % 质量（kg）
k = 10; % 刚度（N/m）
c = 0.1; % 阻尼系数（Ns/m）

% 定义初始条件
x0 = 1; % 初始位移（m）
v0 = 0; % 初始速度（m/s）

% 定义时间范围和步长
t = 0:0.01:10; % 时间范围（s）
dt = t(2) - t(1); % 步长（s）

% 初始化数组
x = zeros(size(t)); % 位移数组
v = zeros(size(t)); % 速度数组
a = zeros(size(t)); % 加速度数组

% 初始条件
x(1) = x0; % 设置初始位移
v(1) = v0; % 设置初始速度

% 基于欧拉法求解
for i = 1:length(t)-1
    a(i) = (-k*x(i) - c*v(i))/m; % 计算加速度
    v(i+1) = v(i) + a(i)*dt; % 更新速度
    x(i+1) = x(i) + v(i)*dt; % 更新位移
end

% 绘制位移-时间图
plot(t, x);
xlabel('时间（s）');
ylabel('位移（m）');
title('单自由度振动响应');

在这个例子中，我们定义了一个质量为1kg的系统，具有10N/m的刚度和0.1Ns/m的阻尼系数。我们通过欧拉法求解差分方程来计算系统在给定条件下的振动响应。在求解过程中，我们迭代计算加速度、速度和位移，并将结果存储在相应的数组中。最后，我们使用`plot`函数绘制了位移-时间图形。