example9_4.m

原名：数字信号处理基础及MATLAB实现,解压缩密码：[e-mail]

%例9-4程序，example9_4
%设计指标
wp=0.2*pi;ws=0.3*pi;
tr_width=ws-wp;           %过渡带宽度
M=100;n=[0:1:M-1];
wc=(ws+wp)/2;             %理想截止频率
%
hd=ideal_lp(wc,M);        %设计理想低通滤波器
w_bl=(Blackman(M))';         %选择布拉克曼窗
h=hd.*w_bl;               %单位脉冲响应
[db,mag,pha,grd,w]=freqz_m(h,[1]);%幅度响应
delta_w=2*pi/1000;
%画图
subplot(1,1,1)
subplot(2,2,1);stem(n,hd);title('理想脉冲响应')
axis([0 M-1 -0.1 0.3]);xlabel('n');ylabel('hd(n)');
subplot(2,2,2);stem(n,w_bl);title('布拉克曼窗(k=3)')
axis([0 M-1 0 1.1]);xlabel('n');ylabel('w(n)');
subplot(2,2,3);stem(n,h);title('实际脉冲响应')
axis([0 M-1 -0.1 0.3]);xlabel('n');ylabel('h(n)');
subplot(2,2,4);plot(w/pi,db);title('幅度响应（单位：dB）')
axis([0 1 -100 10]);xlabel('pi');ylabel('分贝');
%检验
Rp=-(min(db(1:1:wp/delta_w+1)))         %检验实际的通带波动
Ar=-round(max(db(ws/delta_w+1:1:501)))  %检验最小阻带衰减