凯塞窗设计线性相位FIR滤波器(实验4).pdfVIP

Aaaaa实验报告

学院：大数据与信息工程学院

专业：通信工程班级：通信162

姓名学号实验组

实验时间指导教师成绩

2018.12.

26

实验项目名称用凯塞窗设计线性相位FIR滤波器

实

验

通过本实验掌握窗函数法设计线性相位FIR低通滤波器的方法。

目

的

实1、实验预习：前期实验和教材相关内容。

验2、列出本实验编写的所有文件及各项实验结果曲线，加注必要的说明。

要3、对给定指标要求的低通滤波器，理论计算用凯塞窗设计所需的滤波器长

求度和形状参数。

4、总结实验体会及实验中存在的问题。

1、窗函数设计法原理参见《数字信号处理》教材。

2、凯塞窗是对于给定的阻带衰减，是主瓣具有最大能量意义下的最佳窗

函数，因此具有最陡的过渡带。所以通常选用凯塞窗来进行FIR滤波器

的设计。

凯塞窗函数：

其中B为形状参数，IOx是第二类修正零阶贝塞尔函数，其幂函数展开

式为：

计算机一台；

Matlab软件

实

验

仪

器

1、利用凯塞窗来设计线性相位的FIR滤波器。

2、绘出常用窗函数的时域和频域图，可用wvtool实现

实

验

实

内

验

容

步

骤

实一、实验代码

验1.%firlvbo函数，调用dlp.m和kwind.m文件Wp=0.3*pi;

数Wst=0.5*pi;as=40;

据

[wkn1,N]=kwind(as,Wst,Wp);hdn=dlp(Wp,Wst,N);

firjuanji=hdn.*wkn1;figure(1)%时域与频域wvtool(firjuanji)figure(2)%单位冲

击响应impz(firjuanji)figure(3)%幅度谱与相位谱freqz(firjuanji)

2.%计算理想低通滤波器单位取样响应的函数functionhdn=

dlp(Wp,Wst,N)%UNTITLED此处显示有关此函数的摘要%此处显示详细说

明Wc=(Wp+Wst)/2;

a=(N-1)/2;forn=0:N-1if(n==a)

hdn(n+1)=Wc/pi;

else

hdn(n+1)=sin(Wc*(n-a))/(pi*(n-a));

end

end

end

3.

%计算凯塞窗的函数function[wknN]=kwind(as,Wst,Wp)%UNTITLED2此

处显示有关此函数的摘要%此处显示详细说明DW=Wst-Wp;

N=(as-7.95)/(2.286*DW);

%N需要取奇数，才如此选择if(mod(N,2)==0)

N=N+1;

endif(as=50)

B=0.1102*(as-8.7);

elseif(as21)B=0.5842*(as-21)^0.4+0.07886*(as-21);

elseB=0;endendforn=0:N-1

t1=B*sqrt(1-(1-(2*n/(N-1)))^2);IB1=besl(t1);

IB2=besl(B);wkn(n+1)=IB1/IB2;

end

end

4.%计算贝塞尔的函数functions=besl(x)eps=10^(-9);

n=1;

s=1;

D=1;whileD(eps*s)T=x/(2*n);n=n+1;

D