lowpass_fir.m

来自「利用matlab实现fir滤波器的源程序」· M 代码 · 共 31 行

clear;close all;
N=4;
wc=pi/5;                          %理想低通滤波器参数
n=0:N-1;r=(N-1)/2;
hdn=sin(wc*(n-r))/pi./(n-r);      %计算理想低通单位脉冲响应 
if rem(N,2)~=0 hdn(r+1)=wc/pi;end  %N为奇数时，处理n=r点的0/0型
wn1=boxcar(N);                     %矩形窗
hn1=hdn.*wn1';	    %加窗
wn2=hamming(N);	%hamming窗
hn2=hdn.*wn2';	%加窗
w=2*[0:511]/512;
hw1=fft(hn1,512);
hw2=fft(hn2,512);
subplot(2,2,1);                    %作图
stem(n,hn1,'.');grid;
xlabel('n'); ylabel('h(n)');
title('矩形窗设计的h(n) ');        	%图形划分及标注语句
subplot(2,2,2);
plot(w,20*log10(abs(hw1)));grid;
xlabel('w/pi'); ylabel('Magnitude(dB)');
title('幅频特性（dB）');
subplot(2,2,3);
 stem(n,hn2,'.');grid;
xlabel('n'); ylabel('h(n)');
title('hamming窗设计的h(n) ');
subplot(2,2,4);
plot(w,20*log10(abs(hw2)));grid;
xlabel('w/pi'); ylabel('Magnitude(dB)');
title('幅频特性（dB）');