s41.m

数字信号处理的MATLAB仿真加实验指导

%a24.m
wp=0.2*pi;wr=0.4*pi;tr_width=wr-wp
M=ceil(8*pi/tr_width)+1
n=[0:1:M-1];wc=(wr+wp)/2;
hd=ideal_lp(wc,M);
w_ham=(hamming(M))';
h=hd.*w_ham;
[db,mag,pha,w]=freqz_m(h,[1]);
delta_w=2*pi/1000;
Ap=-(min(db(1:1:wp/delta_w+1)))
Ar=-round(max(db(wr/delta_w+1:1:501)))
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_ham);title('海明窗')
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('分贝数');