一、基本理论

二、MATLAB实现

        令式（1.3）中A0 = 2，omega0 = 30，omega = 40，matlab程序如下：

clear;
clc;
close all;

A0 = 2;
omega0 = 30;
omega = 40;
t = 0:0.02:5;
y = A0 * sin( (omega0 - omega) * t /2) .* sin( (omega0 + omega) * t /2);
plot(t, y,'b','LineWidth',1.3);
xlabel('Time/s');
ylabel('Amp/m');
set(gca, 'FontSize', 16);
set(0,'defaultfigurecolor','w');

         当omega = 40时质量块随时间振动的波形如图1所示。

        令式（1.4）中A0 = 1（代指式中幅值部分），omega = 45，matlab程序如下：

clear;
clc;
close all

A0 = 1;
omega = 45;
t = 0:0.02:5;
y = A0 * t.* sin( omega * t) ;
plot(t, y,'r','LineWidth',1.3);
xlabel('Time/s');
ylabel('Amp/m');
set(gca, 'FontSize', 16);
set(0,'defaultfigurecolor','w');

        图2所示为外力频率和固有频率相等时，质量块的振动波形。