%Samp2_17
dt=0.02;df=1/(6000*dt);  %采樣間隔及信號長度為2分鐘的頻率分辨率
n=0:2999;                %折疊頻率之前（取前3000個）數據進行操作
f=n*df;                  %給出頻率序列
sig=rand(1,length(n));   %運用隨機函數產生信號振幅譜
filt=[ones(1,5/df),zeros(1,(length(n)-5/df))];   %理想濾波器幅頻響應函數
[y,n1] = sigmult(filt,n,sig,n);   %信號相乘相當于對信號濾波
subplot(3,1,1),plot(n*df,sig)     %繪出原始信號頻譜
title('原始信號振幅譜') 
subplot(3,1,2),plot(n*df,filt, 'LineWidth',3)    %繪制濾波器的幅頻響應
title('濾波器振幅譜') 
subplot(3,1,3),plot(n1*df,y)      %繪制濾波后的振幅譜
title('濾波后振幅譜'),xlabel('頻率/Hz')