燕大数字课设低通fir滤波器设计与应用.docx
- 文档编号:30760440
- 上传时间:2023-08-20
- 格式:DOCX
- 页数:17
- 大小:228.22KB
燕大数字课设低通fir滤波器设计与应用.docx
《燕大数字课设低通fir滤波器设计与应用.docx》由会员分享,可在线阅读,更多相关《燕大数字课设低通fir滤波器设计与应用.docx(17页珍藏版)》请在冰豆网上搜索。
燕大数字课设低通fir滤波器设计与应用
燕山大学
课程设计说明书
题目:
低通FIR滤波器设计与应用
学院(系):
电气工程学院
年级专业:
10级检测二班
学号:
学生姓名:
***********************
教师职称:
讲师
电气工程学院《课程设计》任务书
课程名称:
数字信号处理课程设计
基层教学单位:
仪器科学与工程系指导教师:
王娜
学号
1
学生姓名
(专业)班级
检测二班
设计题目
5、低通FIR滤波器设计与应用
设
计
技
术
参
数
令x(n)为一正弦信号[f=(1/16Hz)]加白噪声,长度为500,y=fftfilt(h,x)实现卷积运算
设
计
要
求
设计低通FIR滤波器完成滤波运算。
参
考
资
料
数字信号处理方面资料
MATLAB方面资料
周次
前半周
后半周
应
完
成
内
容
收集消化资料、学习MATLAB软件,进行相关参数计算
编写仿真程序、调试
指导教
师签字
基层教学单位主任签字
说明:
1、此表一式四份,系、指导教师、学生各一份,报送院教务科一份。
2、学生那份任务书要求装订到课程设计报告前面。
电气工程学院教务科
摘要··········································································2
第一章滤波器简介···························································3
1.1FIR滤波器简介·····················································3
1.2窗函数简介·························································3
1.3Matlab简介·························································5
1.4窗函数设计FIR滤波器的步骤·········································6
第二章滤波器参数设计·······················································7
第三章程序设计与仿真·······················································7
3.1本设计中使用的Matlab指令···········································7
3.2程序设计····························································8
3.3仿真结果····························································9
3.4仿真结果分析·······················································12
第四章总结与体会··························································12
参考文献····································································13
摘要
数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。
根据其单位冲激响应函数的时域特性可分为两类:
无限脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器。
与IIR滤波器相比,FIR的实现是非递归的,总是稳定的;更重要的是,FIR滤波器在满足幅频响应要求的同时,可以获得严格的线性相位特性。
因此,它在高保真的信号处理,如数字音频、图像处理、数据传输、生物医学等领域得到广泛应用。
MATLAB是“矩阵实验室”(MATrixLABoratoy)的缩写,它是由美国Mathworks公司于1984年正式推出的,是一种面向科学和工程计算的语言,它集数值分析、矩阵运算、信号处理和图形显示于一体,具有编程效率高、调试手段丰富、扩充能力强等特点。
MATLAB的信号处理工具箱具有强大的函数功能,它不仅可以用来设计数字滤波器,还可以使设计达到最忧化,是数字滤波器设计的强有力工具。
关键字:
数字滤波器、信号处理、窗函数、FIR、MATLAB
第一章滤波器的简介
1.1FIR滤波器简介
有限长脉冲响应序列FIR熟悉滤波器有严格的线性相位,任意的幅度特性,工作稳定的优点。
另外,任何非因果有限长序列经过一定的延时都能变成因果的有限长序列,因而FIR滤波器总能用因果系统来实现,而且FIR滤波器由于单位脉冲响应是有限长的,可以用FFT算法来计算,从而大大提高了运算效率。
设h(n)(n=0,1,2,…,N-1)为滤波器的冲激响应,输入信号为x(n),则FIR滤波器就是要实现下列差分方程:
上式就是FIR滤波器的差分方程。
FIR滤波器的最主要的特点是没有反馈回路,因此它是无条件稳定系统。
它的单位脉冲响应h(n)是一个有限长序列。
由上面的方程可见,FIR滤波算法实际上是一种乘法累加运算,它不断地输入样本x(n),经延时做乘法累加,再输出滤波结果。
对式
(1)进行Z变换,整理后可得FIR滤波器的传递函数为:
FIR数字滤波器的设计方法主要有窗函数法和频率抽样设计法,其中窗函数法是基本而有效的设计方法。
FIR数字滤波器的显著特点是线性相位,下面来介绍线性相位的条件。
满足第一类线性相位条件为:
h(n)是实序列且关于n=(N-1)/2偶对称,即
h(n)=h(N-n-1)
满足第二类线性相位条件为:
h(n)是实序列且关于n=(N-1)/2奇对称,即
h(n)=-h(N-n-1)
1.2窗函数简介
窗函数设计法是一种通过截短和计权的方法使无限长非因果序列成为有限长脉冲响应序列的设计方法。
就是根据给定的滤波器技术指标,选择滤波器的阶数N和合适的窗函数w(n)。
即用一个有限长度的窗口函数序列w(n)来截取一个无限长的序列hd(n)获得一个有限长序列h(n)。
通常在设计滤波器之前,应该先根据具体的工程应用确定滤波器的技术指标。
在大多数实际应用中,数字滤波器常常被用来实现选频操作,所以指标的形式一般为在频域中以分贝值给出的相对幅度响应和相位响应。
常用的窗函数有如下几种:
(1)、矩形窗(RectangleWindow):
窗函数为:
其频率响应为:
的主瓣宽度为4
/N,第一旁瓣比主瓣低13dB,最小阻带衰减为21dB。
(2)、汉宁(Hanning)窗,又称升余弦窗:
窗函数为:
当N>>1时,N—1
N,此时,汉宁窗的幅度函数为:
可见,
是由三部分相加得出的,能量更集中在主瓣中。
主瓣宽度为8
/N,第一旁瓣比主瓣低31dB,阻带最小衰减为44dB。
(3)、哈明(Hamming)窗,又称改进的升余弦窗:
窗函数为:
当N>>1时,其幅度响应为:
这种改进的升余弦窗,可将99.96%的能量集中在窗谱的主瓣内,与汉宁窗相比,主瓣宽度相同,但第一旁瓣比主瓣低41dB,最小主瓣衰减为53dB。
(4)、布莱克曼窗(Blackman):
窗函数为:
幅度函数为:
其幅度函数由五部分组成,他们都是移位不同,且幅度也不同的函数,使旁瓣进一步抵消。
阻带衰减进一步增加主瓣宽度是矩形窗的3倍,即12
/N。
(5)凯塞-贝塞尔(Kaiser-Basel)窗:
窗函数为:
,
式中,
,是一个可选参数,用来选择主瓣宽度和旁瓣衰减之间的交换关系,一般说来,β越大,过渡带越宽,阻带越小衰减也越大。
选择窗函数要满足以下两个条件:
(1)窗谱主瓣尽可能地窄,以获得较陡的过渡带;
(2)尽量减少窗谱的最大旁瓣的相对幅度,也就是能量尽量集中于主瓣,使峰肩和纹波减小,就可增多阻带的衰减。
这就给窗函数序列的形状和长度选择。
表1几种常用窗函数的基本参数
1.3Matlab简介
Matlab是一套用于科学计算的可视化高性能语言与软件环境。
它集数值分析、矩阵运算、信号处理和图形显示于一体构成了一个界面友好的用户环境。
它的信号处理工具箱包含了各种经典的和现代的数字信号处理技术是一个优秀的算法研究与辅助设计的工具。
Matlab包含的内容非常丰富,功能强大,可以概括为以下几个方面:
a.可以在多种操作系统下运行,如DOS、Windows95/98/2000/2000/NT、CompaqAlpha、LinuxSunSolaris等。
b.有超过500种的数学、统计、科学及工程方面的函数,使用简单快捷,并且有很强的用户自定义函数的能力。
c.有强大的图形绘制和可视化功能,可以进行视觉数据处理和分析,进行图形、图像的显示及编辑,能够绘制二维、三维图形,使用户可以制作高质量的图形,从而写出图文并茂的文章。
d.有和用其他高级语言(如C,C++,FORTRAN,JAVA)编写的外部程序相接口的能力,也可把MATLAB程序转换成上述高级语言的子程序。
e.有从外部文件及外部硬件设备读入数据的能力。
f.有丰富的网络资源,从相关的Web网站可以直接获得全套的MATLAB联机帮助文件和说明书的电子文档,还可以获得各类技术支持与帮助。
g.有丰富的工具箱〔toolbox〕。
各个领域的专家学者将众多学科领域中常用的算法编写为一个个子程序,即m文件,这些m文件包含在一个个工具箱中。
其工具箱可以分为两大类,即功能性工具箱和科学性工具箱。
功能性工具箱主要用来扩充MATLAB的符号计算、图形可视化、建模仿真、文字处理等功能以及与硬件实时交互的功能。
学科性工具箱是按学科领域来分类的,如信号处理、控制、通信、神经网络图像处理、系统辨识、鲁棒控制、模糊逻辑、小波等工具箱。
1.4窗函数设计FIR滤波器的步骤
下面介绍用窗函数设计FIR滤波器的步骤
1)根据技术要求确定待求滤波器的单位取样响应
,
2)根据对过渡带及阻带衰减的要求,选择窗函数的形式w(n),并估计窗口长度N。
3)计算滤波器的单位取样响应h(n)=w(n)hd(n)。
4)检验指标是否符合要求。
第二章滤波器的设计思路及参数计算
根据任务书中要求,正弦信号的频率为1/16Hz,白噪声长度为500,y=fftfilt(h,x)实现卷积运算。
采用窗函数设计滤波器。
首先,确定截止频率ωc,给出的输入信号为连续信号x(t)=sin2πft,要对其进行采样,x(n)=sin(2πfn/fs),取采样频率fs=21HZ时,得到的有限序列x(n)=sin(2πfn/fs)正好有1.5个周期显示在窗口内,取截止频率ωc=0.007。
长度为N=500。
要求对输入信号加入白噪声,利用Matlab,加入白噪声,带噪信号就是输入滤波器前的信号。
调用五种窗函数,比较其幅度特性。
滤波器的单位取样响应可以直接调用Matlab中的fir指令求得。
进行y(n)=h(n)*x(n)卷积运算。
输出y(n),即为滤波后的信号波形。
第三章程序设计与仿真
本次设计分别应用了矩形窗、三角窗、汉宁窗、布莱克曼窗和凯瑟窗五种窗函数设计FIR低通滤波器,利用Matlab进行滤波器仿真,对带噪信号滤波后的情况进行比对。
3.1本设计中使用的Matlab指令
Pi 圆周率
Plot 线性图形
Figure 建立图形,系统自动从1,2,3,4...来建立图形,数字代表第几幅图形
Randn 正态分布的随机数矩阵
fir1窗函数法设计FIR滤波器
freqz计算数字滤波器的频率响应
holdon 将不同图形在同一坐标系里显示
fftfilt实现卷积运算
legend对显示图形进行字符标注
beta用于控制凯瑟窗的旁瓣宽度,窗函数长度n一定时,beta越大,其频谱的旁瓣就越小,但主瓣宽度相应的增加;当beta一定时,n发生变化,其旁瓣高度不会发生变化。
b0=boxcar(n1)、b1=triang(n1)、b2=hanning(n1)、b3=blackman(n1)、b4=kaiser(n1,beta)分别表示调用矩形窗、三角窗、汉宁窗、布莱克曼窗和凯瑟窗。
3.2程序设计
N=500;%设置白噪声的长度
n=0:
N-1;%设置采样点数
f=1/16;
fs=21;%设置信号的采样频率为21HZ
signal=sin(2*pi*f*n/fs);%设置初始采样正弦信号
noise=randn(1,500);%设置噪声信号长度为500
x=signal+noise;%设置带噪信号,为滤波前的输入信号
figure
(1);
plot(n,signal,'co');%输出初始的正弦信号波形
figure
(2);
plot(n,x);%输出加噪后的正弦信号波形
n1=11;%设置窗函数阶数
beta=10;
b0=boxcar(n1);%分别表示调用矩形窗、三角形窗、汉宁窗、布
b1=triang(n1);莱克曼窗和凯瑟窗
b2=hanning(n1);
b3=blackman(n1);
b4=kaiser(n1,beta);
h0=fir1(n1-1,0.007,b0);%调用fir1函数,窗函数设计截止频率为0.007的h1=fir1(n1-1,0.007,b1);低通滤波器
h2=fir1(n1-1,0.007,b2);
h3=fir1(n1-1,0.007,b3);
h4=fir1(n1-1,0.007,b4);
[h0,f0]=freqz(b0,1,512,2);%计算窗函数的频率响应
[h1,f1]=freqz(b1,1,512,2);
[h2,f2]=freqz(b2,1,512,2);
[h3,f3]=freqz(b3,1,512,2);
[h4,f4]=freqz(b4,1,512,2);
figure(3);%绘制窗函数的幅度谱,在不同坐标系内表示
subplot(2,3,1);
plot(f0,20*log10(abs(h0)));
title('矩形窗');
grid;
subplot(2,3,2);
plot(f1,20*log10(abs(h1)));
title('三角窗');
grid;
subplot(2,3,3);
plot(f2,20*log10(abs(h2)));
title('汉宁窗');
grid;
subplot(2,3,4);
plot(f3,20*log10(abs(h3)));
title('布莱克曼窗');
grid;
subplot(2,3,5);
plot(f4,20*log10(abs(h4)));
title('凯瑟窗');
grid;
figure(4);%计算卷积,绘制家噪信号在滤波后的图像,在不y0=fftfilt(h0,x);同坐标系内表示
y1=fftfilt(h1,x);
y2=fftfilt(h2,x);
y3=fftfilt(h3,x);
y4=fftfilt(h4,x);
plot(n,y0/max(y0));
title('矩形窗');
grid;
subplot(2,3,2);
plot(n,y1/max(y1));
title('三角窗');
grid;
subplot(2,3,3);
plot(n,y2/max(y2));
title('汉宁窗');
grid;
subplot(2,3,4);
plot(n,y3/max(y3));
title('布莱克曼窗');
grid;
subplot(2,3,5);
plot(n,y4/max(y4));
title('凯瑟窗');
grid;
3.3仿真结果
一、输入的正弦信号
图1长度为500的初始正弦信号
图2输入滤波器前的带噪信号
二、五种窗函数的幅度特性图
图4五种窗函数的幅频特性图
四、滤波后的信号
图5滤波后的信号
图6五个输出信号在同一坐标下的对比图
3.4仿真结果分析
由仿真结果可以看出,对于相同的窗口长度N,虽然带噪信号经过低通FIR滤波器的滤波,消除了一部分的噪声,但是输出信号还有一定的失真。
并且,通过不同窗函数设计的FIR低通滤波器,对带噪信号的滤波效果也不相同。
矩形窗FIR低通滤波器的滤波效果最差,布莱克曼和凯瑟窗低通FIR滤波器的滤波器的滤波效果较好。
第四章总结与体会
通过本次课设,我们研究了低通FIR数字滤波的理论知识,为系统整体设计奠定了理论基础,并且研究了MATLAB软件在数字信号处理中的使用方法。
通过利用MATLAB的强大函数运算力和仿真,设计低通FIR数字滤波器,并且通过仿真结果,切实观察到了滤波器对输入信号的滤波效果。
并且,通过相应的程序设计,对MATLAB的设计语言和一部分功能有了初步的了解。
MATLAB信号处理工具箱为滤波器设计及分析提供了非常优秀的辅助设计工具,在设计数字滤波器时,善于应用MATLAB进行辅助设计,能够大大提高设计效率。
数字滤波只是数字信号处理中应用比较广泛的一个领域,在以后的工作和学习生活中,我们需要综合运用模拟电子、数字电子和DSP基本原理等课程中所学到的理论知识去独立完成一个设计项目的设计。
在对本课题的学习和研究的过程,我从中取得了一些成绩,理论水平得到了一定程度的提高,为以后的相关工作积累了经验。
当遇到一个陌生的课题的时候,首先要去查阅大量的相关资料,对课题整体有个初步认知,然后针对具体的问题,进行细致的分析。
另外,要有足够的知识储备,理论一定要扎实,这需要平时的积累。
当然,必要的时候,要虚心向别人请教。
以上是我做课程设计的一些体会,这将对我今后的学习和工作有很大的帮助。
同时感谢指导老师的耐心指导。
由于时间和能力的限制,还存在一些问题和不足之处,有待今后的工作中去研究和思考。
并且希望以后能够获得更多这样实践的机会。
参考文献:
1谢平,王娜,林洪彬著信号处理原理及应用[M].机械工业出版社2008.10
2杨可心,王宝顺,贺志国著MATLB7.X数字信号处理[M]人民邮电出版社2006.12
3维纳·K·英格尔,约翰·G·普洛克斯著刘树棠译数字信号处理(MATLAB版)[M]西安交通大学出版社2008.1
4张德丰著详解MATLAB数字信号处理[M]电子工业出版社2010.6
5丁玉美,高西全著数字信号处理[M]西安电子科技大学出版社2008.6
燕山大学课程设计评审意见表
指导教师评语:
①该生学习态度(认真较认真不认真)
②该生迟到、早退现象(有无)
③该生依赖他人进行设计情况(有无)
平时成绩:
指导教师签字:
2012年6月30日
图面及其它成绩:
答辩小组评语:
①设计巧妙,实现设计要求,并有所创新。
②设计合理,实现设计要求。
③实现了大部分设计要求。
④没有完成设计要求,或者只实现了一小部分的设计要求。
答辩成绩:
组长签字:
2012年6月30日
课程设计综合成绩:
答辩小组成员签字:
2012年6月30日
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字 课设低通 fir 滤波器 设计 应用