example9_6.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 21 行

%例9-6程序，example9_6
%设计指标
Ar=60;beta=0.1102*(Ar-8.7);
M=45;n=[0:1:M-1];
wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;        %理想截止频率
%
w_kai=(kaiser(M,beta))';         %选择凯泽窗
wc1=pi/3;wc2=2*pi/3;
hd=ideal_lp(wc1,M)+ideal_lp(pi,M)-ideal_lp(wc2,M);        %设计滤波器
h=hd.*w_kai;               %单位脉冲响应
[db,mag,pha,grd,w]=freqz_m(h,[1]);%幅度响应
%画图
subplot(1,1,1)
subplot(2,2,1);stem(n,hd);title('理想脉冲响应')
axis([-1 M -0.2 0.8]);xlabel('n');ylabel('hd(n)');
subplot(2,2,2);stem(n,w_kai);title('凯泽窗')
axis([-1 M 0 1.1]);xlabel('n');ylabel('w(n)');
subplot(2,2,3);stem(n,h);title('实际脉冲响应')
axis([-1 M -0.2 0.8]);xlabel('n');ylabel('h(n)');
subplot(2,2,4);plot(w/pi,db);title('幅度响应（单位：dB）')
axis([0 1 -80 10]);xlabel('pi');ylabel('分贝');