samp2_17.m

来自「MATLAB7.x数字信号处理 光盘内容」· M 代码 · 共 14 行

%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')