fcaiyang2.m

FIR数字滤波器频率采样法设计（基于matlab的软件设计）

N=input('输入FIR滤波器的长度N=');
wp=input('输入FIR滤波器的通带截止频率(单位：弧度)wp=');
ws=input('输入FIR滤波器的阻带截止频率(单位：弧度)ws=');
wc=(wp+ws)/2;                      %确定理想低通滤波器的截止频率
k=fix(wc*N/(2*pi)+1);
alpha=(N-1)/2;
Hrs=[ones(1,k),zeros(1,N-2*k+1),ones(1,k-1)];
k1=0:floor((N-1)/2);k2=floor((N-1)/2)+1:N-1;
angH=[-alpha*(2*pi)/N*k1,alpha*(2*pi)/N*(N-k2)];
H=Hrs.*exp(j*angH);
h=real(ifft(H,N));                    %求DFT反变换得到单位冲激响应h(n)
[db,mag,pha,w]=freqz_m(h,[1]);
subplot(1,1,1);
subplot(2,1,1);
n=0:N-1;
stem(n,h);                              %绘制单位冲激响应波形
title('用频率采样法直接设计的低通滤波器脉冲响应');
axis([-1 N-1 -0.1 0.3]);
xlabel('n');ylabel('h(n)');
subplot(2,1,2);
plot(w/pi,db);                        %绘制FIR滤波器的幅频响应
title('滤波器的幅度响应');
axis([0 1 -60 10]);grid;
xlabel('以pi为单位的频率');
ylabel('分贝数／dB');