单自由度震荡matlab

以下是单自由度粘性阻尼系统的自由振动的MATLAB程序示例：

% 定义系统参数
m = 1; % 质量
k = 1; % 刚度
c = 0.1; % 阻尼系数

% 定义初始条件
x0 = 1; % 初始位移
v0 = 0; % 初始速度

% 定义时间范围和时间步长
tspan = [0 50];
dt = 0.01;

% 定义ODE函数
odefun = @(t, y) [y(2); -c/m*y(2)-k/m*y(1)];

% 求解ODE
[t, y] = ode45(odefun, tspan, [x0; v0]);

% 绘制位移-时间图像
plot(t, y(:, 1));
xlabel('Time (s)');
ylabel('Displacement (m)');
title('Single Degree of Freedom Viscously Damped System Free Vibration');

该程序定义了单自由度粘性阻尼系统的质量、刚度和阻尼系数，并定义了初始条件和时间范围。然后，它定义了ODE函数并使用ode45求解ODE。最后，它绘制了位移-时间图像。