winkaeser.m

用matlab编写的数字信号处理程序

clc,clear all,close all;
ws1=0.15*pi;           %给定各边界频率
wp1=0.3*pi;
wp2=0.5*pi;
ws2=0.65*pi;
As=40;                 %As为阻带最小衰减
tr_width=min((wp1-ws1),(ws2-wp2));           %过渡带，两个过渡带应该相等
M=(As-7.95)/(2.285*tr_width)+1;              %滤波器长度M
beta=0.5842*(As-21)^0.4+0.07886*(As-21);     %求出beta值
w_kai=(kaiser(M,beta))';                     %给出凯泽窗
n=[0:1:M];
wc1=(ws1+wp1)/2;                             %两个截止频率
wc2=(ws2+wp2)/2;
hd=ideal_lp(wc2,M+1)-ideal_lp(wc1,M+1);      %ideal_lp函数用来生成低通滤波器单位冲击响应，两个低通想减得到带通
h=hd.*w_kai;                                 %单位冲击与凯泽窗相乘
[db,mag,pha,grd,w]=freqz_m(h,[1]);           %产生绝对值的幅度响应和DB的响应
%依次作出理想脉冲响应，凯泽窗，实际冲击响应和最终的DB响应
subplot(2,2,1);
stem(n,hd);
title('ideal impulse response');
xlabel('n');ylabel('hd(n)');
subplot(2,2,2);
stem(n,w_kai);
title('kaiser window');
xlabel('n');ylabel('kaiser(n)');
subplot(2,2,3);
stem(n,h);
title('actual impulse response');
xlabel('n');ylabel('h(n)');
subplot(2,2,4);
plot(w/pi,db);
title('magnitude response in db');
grid;
xlabel('frequency in pi units');
ylabel('Decibels');