samp7_5.m

MATLAB7.x数字信号处理 光盘内容

%Samp7_5
wp=0.5*pi;ws=0.66*pi;       %滤波器的边界频率
wdelta=ws-wp;            %过渡带宽度
N=ceil(8* pi/wdelta);%求解滤波器的最小阶数，根据Hanning 窗主瓣宽
Wn=(0.5+0.66)*pi/2;%截止频率取通带和阻带边界频率的中点
b=fir1(N,Wn/pi,hanning(N+1));%设计FIR滤波器,注意fir1要求输入归一化频率
[H,f]=freqz(b,1,512,50); %采用50Hz的采样频率求出频率响应
subplot(2,1,1),plot(f,20*log10(abs(H)))
xlabel('频率/Hz');ylabel('振幅/dB');grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel('频率/Hz');ylabel('相位/^o');grid on;