kaiser1.m

很多matlab的源代码

wp = 0.2*pi;
ws = 0.3*pi;
As = 50;
tr_width = ws - wp;
M = ceil((As-7.95)/(14.36*tr_width/(2*pi))+1)+1
n = [0:1:M-1] ;
beta = 0.1102*(As-8.7)
wc = (ws+wp)/2;
hd = ideal_lp(wc,M);
w_kai = (kaiser(M,beta))';
h = hd.*w_kai;
[db,mag,pha,grd,w] = freqz_m(h,[1]);
delta_w = 2*pi/1000;
As = -round(max(db(ws/delta_w+1:1:501)))
subplot(2,2,1);
stem(n,hd);
title('Ideal Response');
axis([0 M-1 -0.1 0.3]);
xlabel('n');
ylabel('hd(n)');
subplot(2,2,2);
stem(n,w_kai);
title('Kaiser Window');
axis([0 M-1 0 1.1]);
xlabel('n');
ylabel('w(n)');
subplot(2,2,3);
stem(n,h);
title('Actual Impulse Response');
axis([0 M-1 -0.1 0.3]);
xlabel('n');
ylabel('h(n)');
subplot(2,2,4);
plot(w/pi,db);
title('Magnitude Response');
axis([0 1 -100 10]);
xlabel('frequency in pi units');
ylabel('Decibels');