sy32.m

数字信号处理工具箱

wp1=0.35*pi;ws1=0.2*pi;
wp2=0.65*pi;ws2=0.8*pi;
M=45;
beta=5.658
n=[0:1:M-1];
wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;
hd=ideal_lp(wc2,M)-ideal_lp(wc1,M);
w_kai=(kaiser(M,beta))';
h=hd.*w_kai;
[db,mag,pha,grd,w]=freqz_m(h,[1]);
delta_w=2*pi/1000;
Rp=-min(db(wp1/delta_w+1:1:wp2/delta_w));
As=-round(max(db(ws2/delta_w+1:1:501)));
subplot(2,2,1);plot(w/pi,db);xlabel('Frequency(单位：pi）');ylabel('|H|:（db)');grid;
subplot(2,2,2);plot(w/pi,pha);xlabel('Frequency(单位：pi）');ylabel('phase');
subplot(2,2,3);stem(n,h);ylabel('h(n)');axis([0 M-1 -0.4 0.5])
subplot(2,2,4);stem(n,w_kai);ylabel('win-kaiser');axis([0 M-1 0 1.1]);