武汉理工大学 基于矩形窗、三角窗、海明窗、汉宁窗、布拉克曼窗的FIR数字滤波器设计.pdfVIP

武汉理工大学基于矩形窗、三角窗、海明窗、汉宁窗、布拉克曼窗的FIR数字滤波器设计--第1页

1.基于矩形窗的FIR数字低通滤波器设计

参数如下：

主程序：%基于矩形窗的FIR数字低通滤波器设计

wp=2*pi*1.75/15;

ws=2*pi*3.25/15

w=(wp+ws)/2;

N=ceil(4*pi/(ws-wp));

n=0:1:(N-1);

alpha=(N-1)/2;

m=n-alpha+eps;

hd=sin(w*m)./(pi*m);%得到理想低通滤波器

B=boxcar(N);

string=[Boxcar,N=,num2str(N)];

h=hd.*(B);%得到FIR数字滤波器

[H,w]=freqz(h,[1],1024);%求其频率响应

db=20*log10(abs(H)+eps);%得到幅值

pha=angle(H);%得到相位

figure(1);

subplot(2,2,1);

stem(n,h,r);

axis([0,N-1,-0.1,0.35]);

xlabel(n);

ylabel(h(n));

title(实际低通滤波器的h(n));

text((0.3*N),0.275,string);

subplot(2,2,2);

plot(w/pi,db,r);

axis([0,1,-100,5]);

xlabel(w/pi);

ylabel(dB);

title(衰减特性(dB));

grid;

subplot(2,2,3);

plot(w,pha);

holdon;

plot(0:4,zeros(5),k);

title(相频特性);

xlabel(频率(rad));

ylabel(相位(rad));

axis([0,3.2,-4,4]);

subplot(2,2,4);

plot(w,abs(H));

title(频率特性);

xlabel(频率W(rad));

ylabel(幅值);

axis([0,3.15,0,1.5]);

text(0.9,1.3,string);

2.基于矩形窗的带通滤波器的设计

参数如下：

武汉理工大学