如何用MATLAB进行高通FIR滤波器仿真

1、实际应用背景：一个信号由5Hz余弦信号与15Hz余弦信号相加构成，如何滤除15Hz分量，得到5Hz分量。

2、产生原始信号：y=cos(2π*5t)+cos(2π*15t)；Matlab程序： f1=5;%第一个点频信号分量频率 f2=15;%第二个点频信号分量频率 fs=100;%采样率 T=2;%时宽 B=10;%FIR截止频率 n=round(T*fs);%采样点个数 t=linspace(0,T,n); y=cos(2*pi*f1*t)+cos(2*pi*f2*t);

3、观察原始信号时域与频谱：频谱中可以观察到信号的5Hz分量和15Hz分量怡觎现喾Matlab程序： figure; plot(t,y); title烫喇霰嘴('原始信号时域'); xlabel('t/s'); ylabel('幅度'); figure; fft_y=fftshift(fft(y)); f=linspace(-fs/2,fs/2,n); plot(f,abs(fft_y)); title('原始信号频谱'); xlabel('f/Hz'); ylabel('幅度'); axis([ 0 50 0 100]);

4、生成高通FIR滤波器冲击响应：MATLAB程序： b=fir1(80, B/(fs/2),'high'); %滤波产生指定带宽的噪声信号 figure; freqz(b);

5、原始信号经过高通滤波器： MATLAB程序： y_after_fir=filter(b,1,y);

6、滤波后信号时域及频谱：滤波后得到15Hz的疙熳阊涓信号，5Hz信号被滤除。MATLAB程序： figure; plot(t,y_after_fir); title('滤波后信号时域'); xlabel('t/s'); ylabel('幅度'); fft_y1=fftshift(fft(y_after_fir)); f=linspace(-fs/2,fs/2,n); figure; plot(f,abs(fft_y1)); title('滤波后信号频谱'); xlabel('f/Hz'); ylabel('幅度'); axis([ 0 50 0 100]);