MATLAB-振动问题：单自由度阻尼振动系统受迫振动

一、基本理论

二、MATLAB实现 

        单自由度阻尼振动系统受迫振动，MATLAB代码如下：

clear;
clc;
close all

A = 1;
psi = 0;
F0 = 10;
D = 20;
Rm = 0.5;
M = 1;
omega = 2;
delta = Rm / (2*M);
omega0 = sqrt(D / M);
Omega = sqrt(omega0^2 - delta^2);
Zm = Rm + i * (M * omega - D/M);
Zm0 = abs(Zm);
phi = angle(Zm);
t = 0:0.01:20;
y = A * exp(-delta * t) .* cos(Omega * t + psi) + (F0/(omega * Zm0)) * sin(omega * t - phi);

plot(t, y,'r','LineWidth',1.5);
xlabel('Time/s');
ylabel('Amp/m');
set(gca, 'FontSize', 16);
set(0,'defaultfigurecolor','w');

        当阻力系数Rm = 0.5时，质量块随时间振动的波形如图1所示。

图1