example3.m

来自「MATLAB数字信号处理的入门学习例程及实验结果」· M 代码 · 共 23 行

wp=0.2*pi;ws=0.3*pi;
tr_width=ws-wp;
M=ceil(6.6*pi/tr_width)+1;
N=[0:1:M-1];
wc=(ws+wp)/2;
hd=ideal_lp(wc,M);
w_ham=(hamming(M))';
% hd=fir1(M-1,wc,hamming(M));
% h=hd;
h=hd.*w_ham;
[db,mag,pha,grd,w]=freqz_m(h,[1]);
delta_w=2*pi/1000;
Rp=-(min(db(1:1:wp/delta_w+1)));
As=-round(max(db(ws/delta_w+1:1:501)));
close all;
subplot(2,2,1);stem(hd);title('理想冲击响应')
axis([0 M-1 -0.1 0.3]);ylabel('hd[n]');
subplot(2,2,2);stem(w_ham);title('汉明窗');
axis([0 M-1 0 1.1]);ylabel('w[n]');
subplot(2,2,3);stem(h);title('实际冲击响应');
axis([0 M-1 -0.1 0.3]);ylabel('h[n]');
subplot(2,2,4);plot(w/pi,db); title('衰减幅度');
axis([0 1 -100 10]);ylabel('Decibles');