数字信号处理课程设计报告基于inversesincFIR的数字低通滤波器.docx
- 文档编号:23142939
- 上传时间:2023-05-08
- 格式:DOCX
- 页数:19
- 大小:1.87MB
数字信号处理课程设计报告基于inversesincFIR的数字低通滤波器.docx
《数字信号处理课程设计报告基于inversesincFIR的数字低通滤波器.docx》由会员分享,可在线阅读,更多相关《数字信号处理课程设计报告基于inversesincFIR的数字低通滤波器.docx(19页珍藏版)》请在冰豆网上搜索。
数字信号处理课程设计报告基于inversesincFIR的数字低通滤波器
南华大学
数字信号处理课程设计
学院:
电气工程学院
学生姓名:
潘艺
专业班级:
电子091班
题目编号:
0204
设计题目:
inversesinc技术低通滤波器
起止时间:
2012年12月26日至2013年1月11日
指导老师:
陈忠泽
2013年1月15日
一、设计要求与目的
了解和掌握设计InverseSincFIR的数字低通滤波器的原理和低通滤波器的原理。
基于InverseSincLowpass思想设计一个数字低通滤波器。
数字低通滤波器的设计指标:
(1)通带截止频率错误!
未找到引用源。
=0.49192πrad,
(2)过渡带宽度错误!
未找到引用源。
=0.09πrad,
(3)滚降错误!
未找到引用源。
=36dB,id=45。
二、InverseSincFIR的数字低通滤波器设计思想
1、基本原理
(1)FIR滤波器简介
FIR滤波器(有限脉冲响应滤波器)是一种既能够保证幅度特性满足技术要求,又能够做到有严格的线性相位特性。
其系统函数为:
H(Z)=
H(Z)是
的N-1的多次项,他在Z平面内有N-1个零点,在原点有一个N-1重极点。
因此H(Z)永远稳定。
稳定和线性相位是FIR滤波器最突出的特点。
数字滤波器一般需要用递归模型来实现,因而有时也称之为递归滤波器。
FIR滤波器的冲激响应只能延续一定时间,在工程实际中可以采用递归的方式实现,也可以采用非递归的方式实现。
3、设计内容
1确定“数字低通滤波器”设计指标;
2设计InverseSincLowpass;
而数字低通滤波器的性能指标:
(1)通带截止频率错误!
未找到引用源。
=0.49192πrad,
(2)过渡带宽度错误!
未找到引用源。
=0.09πrad,
(3)滚降错误!
未找到引用源。
=36dB,id=45。
实际带宽减去理想带宽/理想带宽=滚降系数,阻带滚降越厉害,阻带内的最小衰减就越大。
总结数字低通的设计步骤:
设fs=8000Hz
通带截止频率ωpc=0.49192πrad,
=1.96768kHz
过渡带宽度错误!
未找到引用源。
,
阻带起始频率
=0.09πrad+0.49192πrad=0.58192πrad,
b).因为H
(e
)={
首先由所需低通滤波器的过渡带求截止频率
=
(
+
)=2147.68Hz
其对应的数字频率为
=
=2
=0.53692
通带最大衰减为
,阻带最小衰减为
通带最大衰减ap=1dB
阻带最小衰减as=36dB
%Fp —frequencyatthestartofthepassband.Specifiedinnormalizedfrequencyunits.AlsocalledFpass.
%Fst —frequencyattheendofthestopband.Specifiedinnormalizedfrequencyunits.AlsocalledFstop.
d=fdesign.isinclp(0.49192,0.58192,1,36);%前两个频率需归一化和后两个单位dB
hd=design(d,'equiripple');
fvtool(hd);
求得h(n)长度为33,系数如下:
h=[0.006378391-0.015779155-0.0242883170.0011287610.017887846-0.010458761
-0.022624740.019467660.026849823-0.033645542-0.0310519360.057605389
0.035526688-0.108429032-0.0454730890.3219756940.5529945870.321975694
-0.045473089-0.1084290320.0355266880.057605389-0.031051936-0.0336455420.0268498230.01946766-0.02262474-0.0104587610.0178878460.001128761
-0.024288317-0.0157791550.006378391]
其中,
h(0)=0.0063783905466236416=h(32)h
(1)=-0.015779155415281004=h(31)
h
(2)=.024*********=h(30)h(3)=0.0011287610293565845=h(29)
……
h(15)=0.32197569395966075=h(17)h(16)=0.55299458675125945=h(16)
则设计的InverseSincFIR的数字低通滤波器的系统函数为:
根据上述所求h(n)即可求得
脉冲响应
阶跃响应
极零点图
三、滤波器的不同结构对性能指标的影响
在理想状态下,对于同一个传递函数几乎对应着无数种等效结构,然而这些结构却并不一定都能实现。
在无限参数字长的情况下,所有能实现传递函数的结构之间,其表现完全相同。
然而,在实际中,由于参数字长有限的限制,各实现结构的表现并不相同。
下面我们就将对比直接型线性相位结构和直接型多相滤波结构在本例中对性能指标的影响。
在MATLAB中可以利用FDATOOL工具箱构建不同类型的数字滤波器。
在此为了使对比效果明显,我们不妨先将将上述初步设计的InverseSincFIR的数字低通滤波器的设计参数的字长(即转移函数中的系数)进行保留小数点后3位的进一步的缩减。
缩减后的参数如下:
Numerator=[0.0280.0530.0710.0530.028]
Denominator=[1.000-2.0262.148-1.1590.279]
将上述参数输入FDATOOL中的filtercoefficients工具中如下图所示。
filtercoefficients工具工作界面
1、利用直接型结构构建数字滤波器
Ghp(z)=
Y(n)=[0.028x(n)+0.053x(n-1)+0.071x(n-2)+0.053x(n-3)+0.028x(n-4)]+[x(n)-2.026x(n-1)+2.148x(n-2)-1.159x(n-3)+0.279x(n-4)]
x(n)0.0281y(n)
0.0532.026
0.071-2.148
0.0531.159
0.028-0.279
直接型I型结构流图
直接型线性相位结构滤波器的典型频率响应表达式为:
选择filterstructure选项框中的Direct-FormITransposed选项,点击窗口下方的ImportFilter按钮,构建直接型线性相位结构的InverseSincFIR的数字低通滤波器,结果如图所示:
Direct-FormI结构的滤波器幅频响应图
读图可以得直接型多相滤波结构的滤波器技术指标(fpc,
ftz,单位为kHz;
单位为dB)如表1所示:
表1Direct-FormI结构滤波器对性能指标的影响
性能指标
初始设计指标
直接型线性相位结构
△
fpc
1.96768
1.967285
-0.000395
ftz
0.36
0.346192
-0.013808
36
37.14677
1.14677
从表1中可以看出,采用直接型多相滤波结构构成的滤波器:
滤波器幅频曲线在通带和阻带内波动较均匀。
出现了衰减
上升了1.14677dB,
fpc,
ftz分别较初始设计分别减少了0.000395kHz和0.013808kHz。
2、利用级联结构构建数字滤波器
将Ghp(z)=
因式分解,
symsz
factor((0.028*z^0+0.053*z^(-1)+0.071*z^(-2)+0.053*z^(-3)+0.028*z^(-4))/(1.000*z^0-2.026*z^(-1)+2.148*z^(-2)-1.159*z^(-3)+0.279*z^(-4)))
Ghp(z)=0.028(1+1.5464*z^(-1)+z^(-2))(1+0.3466*z^(-1)+z^(-2))/((1-1.0646*z^(-1)+0.7549z^(-2))(1-1.0646*z^(-1)+0.3696z^(-2))
Ghp(z)=
得到:
x(n)10.02811y(n)
1.06460.04329921.06460.3466
-0.75490.0280.36961
级联型的结构流图
选择Edit下拉菜单中点击ConverttoSecond-orderSections选项,将构建好的Direct-FormI结构的InverseSincFIR的数字低通滤波器转换为级联滤波器,结果如下图所示。
级联型结构滤波器幅频响应图
读图可以得直接型线性相位结构的滤波器技术指标(fpc,
ftz,单位为kHz;
单位为dB)如表2所示:
表2级联型结构滤波器对性能指标的影响
性能指标)
初始设计指标
直接型线性相位结构
△
fpc
1.96768
1.967773
0.000093
ftz
0.36
0.355469
-0.004531
36
36.12895
0.12895
由图3和表1可以看出:
滤波器幅频曲线在通带和阻带内波动较均匀。
出现了衰减
上升了0.12895dB,
fpc,
ftz分别较初始设计分别增加了0.000093kHz和减少了0.004531kHz。
两种结构滤波器对指标影响比较与原因分析
比较表1和表2发现:
在参数字长仅保留了小数点后3位的情况下,两种结构的滤波器较初始设计在性能指标方面均有误差。
但是级联型结构滤波器结构误差受有限参数字长的影响较Direct-FormI更小一些,主要表现在:
直接型线性相位结构和
与设计要求相应的性能间的差的绝对值。
此外,直接型和级联型的幅频响应曲线的通带的波动均不稳定,这是其一大缺点。
造成这一现象的原因是:
直接型线性相位结构滤波器的系数不是直接决定单个零极点,因而不能很好的进行滤波器性能的控制;此外直接型线性相位结构滤波器的极点对参数的变化过于敏感,从而使得系统的频率响应对参数的变化也特别敏感,也就是对参数的有限字长运算过于灵敏,容易出现不稳定或产生较大误差。
小结
在受实际条件限制,参数字长有限的情况下,级联型结构滤波器结构对参数变化的反应要比Direct-FormI结构的更低,性能指标误差更小,滤波器失真更小,滤波效果更好,更能符合设计指标的要求。
四、参数字长对性能指标的影响
在实际的数字滤波器的设计中,由于计算机或DSP芯片等的字长和存储空间有限,所以也只能对设计参数取有限的字长进行设计。
然而,如果字长太短,则设计的滤波器误差就会太大,造成滤波效果不佳。
下面就将以上述性能指标为依据,初始设计的InverseSincFIR的数字低通滤波器为例,研究不同参数字长对性能指标的影响。
并为合适参数字长的确定探索规律。
将计算获得的低通数字滤波器的系数输入FDATOOL的filtercoefficients工具中,并点击ImportFilter按钮,生成数字滤波器。
在Direct-FormITransposed结构滤波器下,运用FDATOOL左下侧一列按键中的
按钮,在filterarithmetic下拉菜单下选择FixedPoint选项,进入下图所示的界面。
Setquantizationparameters工作界面
通过改变numeratorwordlength的值便可以改变参与构建滤波器的参数字长。
1、参数字长取2位对性能指标的影响
将numeratorwordlength、numeratorfrac.length的值均改为2,点击下方的Apply按钮,此时设计的滤波器幅频响应曲线和性能指标如下图所示:
参数字长取2位时的滤波器幅频响应曲线图
图7中的虚线为供参考的理想字长下生成的滤波器的幅频响应曲线,图中实线为参数字长取为2位时的滤波器幅频响应曲线。
从图中可以看出:
字长为2位时,滤波器幅频响应曲线是条直线,远远不能满足设计指标的要求。
2、参数字长取4位对性能指标的影响
将numeratorwordlength、numeratorfrac.length的值均改为4,点击下方的Apply按钮,此时设计的滤波器幅频响应曲线和性能指标如下图所示:
参数字长取4位时的滤波器幅频响应曲线图
虚线为供参考的理想字长下生成的滤波器的幅频响应曲线,图中实线为参数字长取为4位时的滤波器幅频响应曲线。
从图中可以看出:
字长为4位时,滤波器的各项性能指标离设计指标偏差很大,滤波器失真明显,滤波效果很差,远远不能满足设计指标的要求。
3、参数字长取6位对性能指标的影响
将numeratorwordlength、numeratorfrac.length的值均改为6,点击下方的Apply按钮,此时设计的滤波器幅频响应曲线和性能指标如下图所示:
参数字长取6位时的滤波器幅频响应曲线图
由图可以看出,当参数字长取为6位时,幅频曲线失真进一步减小,已经初具低通滤波器的形制。
但性能指标远未达到设计要求,误差还是很大。
4、参数字长取8位对性能指标的影响
将numeratorwordlength、numeratorfrac.length的值均改为8,点击下方的Apply按钮,此时设计的滤波器幅频响应曲线和性能指标如下图所示:
参数字长取8位时的滤波器幅频响应曲线图
由图可以看出,当参数字长取为8位时,幅频曲线失真进一步减小,但仍可以看出。
性能指标与设计要求的差值继续减小,也已不明显;通带最大衰减频率和阻带最小衰减频率与设计指间标误差还是有一定的差距,距离设计要求仍可以进一步接近。
5、参数字长取10位对性能指标的影响
将numeratorwordlength、numeratorfrac.length的值均改为10,点击下方的Apply按钮,此时设计的滤波器幅频响应曲线和性能指标如下图所示:
参数字长取10位时的滤波器幅频响应曲线图
由图看出,当参数字长取为10位时,幅频曲线失真进一步减小,但程度已经很小,失真几乎可以忽略。
截止频率改善的程度也变得很小;通带最大衰减和阻带最小衰减仍然和设计要求仍有一定的差距,不可以忽略。
6、参数字长取12位对性能指标的影响
将numeratorwordlength、numeratorfrac.length的值均改为12,点击下方的Apply按钮,此时设计的滤波器幅频响应曲线和性能指标如下图所示:
参数字长取12位时的滤波器幅频响应曲线图
由图看出,当参数字长取为12位时,幅频曲线失真进一步减小,但程度已经很小,失真几乎可以忽略。
截止频率改善的程度也变得很小;通带最大衰减和阻带最小衰减仍然和设计要求仍有一定的差距,不可以忽略。
7、参数字长取14位对性能指标的影响
将numeratorwordlength、numeratorfrac.length的值均改为14,点击下方的Apply按钮,此时设计的滤波器幅频响应曲线和性能指标如下图所示:
参数字长取14位时的滤波器幅频响应曲线图
由图看出,当参数字长取为14位时,幅频曲线失真进一步减小,但程度已经很小,失真几乎可以忽略。
截止频率改善的程度也变得很小;通带最大衰减和阻带最小衰减仍然和设计要求仍有一定的差距,不可以忽略。
8、参数字长取16位对性能指标的影响
将numeratorwordlength、numeratorfrac.length的值均改为16,点击下方的Apply按钮,此时设计的滤波器幅频响应曲线和性能指标如下图所示:
参数字长取16位时的滤波器幅频响应曲线图
由图可以看出,当参数字长取为16位及以上时,幅频曲线失真几乎为零,设计的曲线与要求的曲线几乎完全重合。
截止频率,同组带衰减也与设计要求几乎完全相同。
设计的滤波器各项性能指标达到设计要求。
小结
总结以上8个步骤发现:
参数字长越长,设计出的滤波器就越符合设计指标要求,误差越小,稳定性越好。
且当参数字长达到16位及以上位时,设计的滤波器便可达到设计性能指标。
总结
通过本课程设计,我掌握了InverseSincLowpass滤波器的设计原理以及步骤,学会了应用MATLAB编写程序,学会了其强大的fadtool工具和InverseSincFilterDesignDialogBox,熟悉了它的使用方法与注意事项。
在设计过程中,发现有很多知识点存在盲点,在查阅了相关资料后方有所感悟,数字信号处理博大精深,不是一朝一夕就能完全掌握,还有待于老师的孜孜教诲,学习需要与实践相结合,自己推公式是个很有必要的过程,通过本次设计,我受益匪浅,好好学习,知识体系结构还有待完善。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 课程设计 报告 基于 inversesincFIR 数字 滤波器