‘拍’现象的MATLAB模拟

clear 
clc
fs=10000
t=0:1/fs:20;
f1=151;%输入振动1的频率f1
w1=2*pi*f1;   
f2=152;%输入振动2的频率f2
w2=2*pi*f2;  
y1=2*cos(w1*t);
y2=2*cos(w2*t);
y=y1+y2;

figure
subplot(3,1,1)      %显示振动1的图像
plot(t,y1)
axis([0,0.5,-2,2])
ylabel('y1')
sound(y1,fs);       %播放振动1的声音
pause(10);

subplot(3,1,2)     %显示振动2的那图像
plot(t,y2)
axis([0,0.5,-2,2])
ylabel('y2')
sound(y2,fs);       %播放振动2的声音
pause(10);

subplot(3,1,3)     %显示合振动的图像
plot(t,y)
axis([0,3,-4,4])
ylabel('y')
xlabel('t')

sound(y,fs)         %播放合振动的声音

拍现象：频率较大而频率之差很小的两个同方向简谐运动的合成，其合振动的振幅时而加强时而减弱的现象.

拍的数学表达：