基于频率抽样法的FIR数字低通滤波器的设计
- 格式:doc
- 大小:564.94 KB
- 文档页数:12
实验四FIR数字滤波器的设计
FIR(有限冲击响应)数字滤波器是一种常见的数字信号处理器件,
可以用于滤波、降噪等应用。
下面是一种FIR数字滤波器的设计流程:
1.确定滤波器的需求:首先确定需要滤除的频率范围和滤波的类型,
例如低通、高通、带通、带阻等等。
2.设计滤波器的频率响应:根据滤波器的需求,设计其理想的频率响应。
可以使用窗函数、最小二乘法等方法获得一个理想的滤波器响应。
3.确定滤波器的阶数:根据设计的频率响应,确定滤波器的阶数。
阶
数越高,滤波器的响应越陡峭,但计算复杂度也会增加。
4.确定滤波器的系数:根据滤波器的阶数和频率响应,计算滤波器的
系数。
可以使用频域窗函数或时域设计方法。
5.实现滤波器:根据计算得到的滤波器系数,实现滤波器的计算算法。
可以使用直接形式、级联形式、传输函数形式等。
6.评估滤波器的性能:使用所设计的FIR滤波器对输入信号进行滤波,评估其滤波效果。
可以使用频率响应曲线、幅频响应、群延时等指标进行
评估。
7.调整滤波器设计:根据实际的滤波效果,如果不满足需求,可以调
整滤波器的频率响应和阶数,重新计算滤波器系数,重新实现滤波器。
以上是FIR数字滤波器的基本设计流程,设计过程中需要考虑滤波器
的性能、计算复杂度、实际应用需求等因素。
频率采样法设计fir滤波器
频率采样法设计FIR滤波器是一种在实际应用中非常有用的方法,它可以有效地实现滤波器的设计,并且能够得到良好的性能。
这种方法通过采样系统的输入信号来确定最佳滤波器设计,这些采样点是通过测量输入信号的功率谱密度函数(PSD)来确定的。
在频率采样法设计FIR滤波器的过程中,首先需要测量输入信号的PSD,这一步就是确定采样点的关键,因为这些采样点将作为滤波器设计的基石。
然后,需要使用Fourier变换来根据所采样的PSD来计算滤波器的频率响应,这一步也是决定滤波器特性的重要环节。
最后,需要使用反向FT算法来计算所需的滤波器系数,以实现滤波器的设计。
在频率采样法设计FIR滤波器的过程中,通常使用大量的采样点,以便能够更准确地表示信号的PSD,从而让滤波器的性能更好。
当采样点越多时,滤波器的响应就会变得更加精确,而且可以得到更低的相位延迟,从而使其具有更好的性能。
在实际应用中,频率采样法设计FIR滤波器通常能够得到很好的效果,但也存在一些不足之处。
首先,它所需要的采样点数量可能会比较多,这可能会增加设计的复杂
度,从而降低滤波器的性能。
其次,由于实际信号的PSD 可能受到噪声的影响,因此采样点的准确性也可能会受到影响,从而影响滤波器的性能。
总之,频率采样法设计FIR滤波器是一个实用的方法,它可以有效地实现滤波器的设计,但也存在一些不足之处,因此在实际应用中,必须根据实际情况来进行适当的取舍。
FIR数字滤波器的设计与实现介绍在数字信号处理中,滤波器是一种常用的工具,用于改变信号的频率响应。
FIR (Finite Impulse Response)数字滤波器是一种非递归的滤波器,具有线性相位响应和有限脉冲响应。
本文将探讨FIR数字滤波器的设计与实现,包括滤波器的原理、设计方法和实际应用。
原理FIR数字滤波器通过对输入信号的加权平均来实现滤波效果。
其原理可以简单描述为以下步骤: 1. 输入信号经过一个延迟线组成的信号延迟器。
2. 延迟后的信号与一组权重系数进行相乘。
3. 将相乘的结果进行加和得到输出信号。
FIR滤波器的特点是通过改变权重系数来改变滤波器的频率响应。
不同的权重系数可以实现低通滤波、高通滤波、带通滤波等不同的滤波效果。
设计方法FIR滤波器的设计主要有以下几种方法:窗函数法窗函数法是一种常用简单而直观的设计方法。
该方法通过选择一个窗函数,并将其与理想滤波器的频率响应进行卷积,得到FIR滤波器的频率响应。
常用的窗函数包括矩形窗、汉宁窗、哈密顿窗等。
不同的窗函数具有不同的特性,在设计滤波器时需要根据要求来选择合适的窗函数。
频率抽样法频率抽样法是一种基于频率抽样定理的设计方法。
该方法首先将所需的频率响应通过插值得到一个连续的函数,然后对该函数进行逆傅里叶变换,得到离散的权重系数。
频率抽样法的优点是可以设计出具有较小幅频纹波的滤波器,但需要进行频率上和频率下的补偿处理。
最优化方法最优化方法是一种基于优化理论的设计方法。
该方法通过优化某个性能指标来得到最优的滤波器权重系数。
常用的最优化方法包括Least Mean Square(LMS)法、Least Square(LS)法、Parks-McClellan法等。
这些方法可以根据设计要求,如通带波纹、阻带衰减等来得到最优的滤波器设计。
实现与应用FIR数字滤波器的实现可以通过硬件和软件两种方式。
硬件实现在硬件实现中,可以利用专门的FPGA(Field-Programmable Gate Array)等数字集成电路来实现FIR滤波器。
FIR数字滤波器设计实验_完整版
在FIR数字滤波器设计实验中,我们需要完成以下步骤:
1.确定滤波器的规格:包括滤波器的类型(低通、高通、带通或带阻)、截止频率、通带波纹、阻带衰减等。
2.选择适当的滤波器设计方法:常见的设计方法包括窗函数法、频率抽样法等。
3.根据选择的设计方法,计算滤波器的系数。
4.实现滤波器:根据计算得到的系数,编写程序在计算机或嵌入式系统中实现滤波器。
5.对输入信号进行滤波处理:将需要滤波的信号输入到滤波器中,获得滤波后的输出信号。
6.评估滤波效果:通过对比输入和输出信号,评估滤波器的性能,包括频率响应、相位响应、时域响应等。
完成FIR数字滤波器设计实验需要具备一定的信号处理和数字滤波器设计的知识,以及一些编程和实验能力。
实验中通常会使用MATLAB、Python等工具进行滤波器设计和信号处理的仿真和实现。
这样的实验对于学习信号处理和数字滤波器设计非常有帮助,可以加深对理论知识的理解,并锻炼实际应用的能力。
实验八频率采样法设计FIR数字滤波器一、实验目的掌握频率取样法设计FIR数字滤波器的原理及具体方法。
二、实验设备与环境计算机、MATLAB软件环境三、实验基础理论1.基本原理频率取样法从频域出发,把理想的滤波器等间隔取样得到,将作为实际设计滤波器的,N-1得到以后可以由来唯一确定滤波器的单位脉冲响应,()D_Dd___________ðϨϨ________________求得其中为内插函数由求得的频率响应来逼近。
如果我们设计的是线性相位FIR滤波器,则的幅度和相位一定满足线性相位滤波器的约束条件。
我们将表示成如下形式当为实数,则由此得到即以k=N/2为中心呈偶对称。
再利用线性条件可知,对于1型和2型线性相位滤波器对于3型和4型线性相位滤波器其中,表示取小于该数的最大的整数。
2.设计步骤(1)由给定的理想滤波器给出和。
(2)由式求得。
(3)根据求得和。
四、实验内容1.采用频率采样设计法设计FIR数字低通滤波器,满足以下指标(1)取N=20,过渡带没有样本。
(2)取N=40,过渡带有一个样本,T=0.39。
(3)取N=60,过渡带有两个样本,T1=0.5925,T2=0.1009。
(4)分别讨论采用上述方法设计的数字低通滤波器是否能满足给定的技术指标。
实验代码与实验结果(1)N=20 过渡带没有样本N=20;alpha=(N-1)/2;l=0:N-1;wl=(2*pi/N)*l;Hrs=[1,1,1,zeros(1,15),1,1]; *对理想幅度函数取样得到取样样本Hdr=[1,1,0,0];wdl=[0,0.25,0.25,1]; *用于绘制理想函数幅度函数的曲线k1=0:floor((N-1)/2);k2=floor((N-1)/2)+1:N-1;angH=[-alpha*(2*pi)/N*k1,alpha*(2*pi)/N*(N-k2)];H=Hrs.*exp(j*angH); *计算H(k)h=ifft(H,N); *计算h(n)w=[0:500]*pi/500;H=freqz(h,1,w); *计算幅度响应[Hr,wr]=zerophase(h); *计算幅度函数subplot(221);plot(wdl,Hdr,wl(1:11)/pi,Hrs(1:11),'o');axis([0,1,-0.1,1.1]);xlabel('\omega(\pi)');ylabel('Hr(k)');subplot(222);stem(l,h,'filled');axis([0,N-1,-0.1,0.3]);xlabel('n');ylabel('h(n)');subplot(223);plot(wr/pi,Hr,wl(1:11)/pi,Hrs(1:11),'o');axis([0,1,-0.2,1.2]);xlabel('\omega(\pi)');ylabel('Hr(w)');subplot(224);plot(w/pi,20*log10((abs(H)/max(abs(H)))));axis([0,1,-50,5]);grid;xlabel('\omega(\pi)');ylabel('dB');0.51ω(π)H r (k )051015nh (n )0.51ω(π)H r (w )0.51-40-20ω(π)d B(2)N=40 过渡带有一个样本,T=0.39 N=40;alpha=(N-1)/2;l=0:N-1;wl=(2*pi/N)*l;Hrs=[1,1,1,1,1,0.39,zeros(1,29),0.39,1,1,1,1]; *设置过渡带样本 Hdr=[1,1,0.39,0,0];wdl=[0,0.2,0.25,0.3,1]; k1=0:floor((N-1)/2); k2=floor((N-1)/2)+1:N-1;angH=[-alpha*(2*pi)/N*k1,alpha*(2*pi)/N*(N-k2)]; H=Hrs.*exp(j*angH); h=ifft(H,N);w=[0:500]*pi/500; H=freqz(h,1,w);[Hr,wr]=zerophase(h); subplot(221);plot(wdl,Hdr,wl(1:21)/pi,Hrs(1:21),'o'); axis([0,1,-0.1,1.1]); xlabel('\omega(\pi)'); ylabel('Hr(k)'); subplot(222);stem(l,h,'filled'); axis([0,N-1,-0.1,0.3]); xlabel('n');ylabel('h(n)'); subplot(223);plot(wr/pi,Hr,wl(1:21)/pi,Hrs(1:21),'o'); axis([0,1,-0.2,1.2]); xlabel('\omega(\pi)'); ylabel('Hr(w)'); subplot(224);plot(w/pi,20*log10((abs(H)/max(abs(H))))) axis([0,1,-80,5]);grid;xlabel('\omega(\pi)'); ylabel('dB');0.51ω(π)H r (k )0102030nh (n )0.51ω(π)H r (w )0.51-80-60-40-200ω(π)d B(3)N=60 过渡带有两个样本 T1=0.5925,T2=0.1009 N=60;alpha=(N-1)/2;l=0:N-1;wl=(2*pi/N)*l;Hrs=[1,1,1,1,1,1,1,0.5925,0.1099,zeros(1,43), 0.1099, 0.5925, 1,1,1,1,1,1]; *设置过渡带样本Hdr=[1,1,0.5925,0.1099,0,0];wdl=[0,0.2,0.2+1/30,0.3-1/30,0.3,1]; k1=0:floor((N-1)/2); k2=floor((N-1)/2)+1:N-1;angH=[-alpha*(2*pi)/N*k1,alpha*(2*pi)/N*(N-k2)]; H=Hrs.*exp(j*angH); h=ifft(H,N);w=[0:500]*pi/500;H=freqz(h,1,w);[Hr,wr]=zerophase(h); subplot(221);plot(wdl,Hdr,wl(1:31)/pi,Hrs(1:31),'o'); axis([0,1,-0.1,1.1]);xlabel('\omega(\pi)');ylabel('Hr(k)'); subplot(222);stem(l,h,'filled'); axis([0,N-1,-0.1,0.3]); xlabel('n');ylabel('h(n)'); subplot(223);plot(wr/pi,Hr,wl(1:31)/pi,Hrs(1:31),'o'); axis([0,1,-0.2,1.2]); xlabel('\omega(\pi)'); ylabel('Hr(w)'); subplot(224);plot(w/pi,20*log10((abs(H)/max(abs(H))))) axis([0,1,-120,5]);grid;xlabel('\omega(\pi)');ylabel('dB');0.5100.51ω(π)H r (k )02040-0.10.10.20.3nh (n )0.5100.51ω(π)H r (w )0.51-100-50ω(π)d B(4)答:由实验结果第四个图可知,当时,阻带增益都没有达到-50dB,阻带增益有所减低,所以设计结果不能满足最初的设计要求。
FIR数字滤波器的设计
FIR(有限冲激响应)数字滤波器的设计主要包括以下几个步骤:
1.确定滤波器的要求:根据应用需求确定滤波器的类型(如低通、高通、带通、带阻等)和滤波器的频率特性要求(如截止频率、通带波动、阻带衰减等)。
2.确定滤波器的长度:根据频率特性要求和滤波器类型,确定滤波器的长度(即冲激响应的系数个数)。
长度通常根据滤波器的截止频率和阻带宽度来决定。
3.设计滤波器的冲激响应:使用一种滤波器设计方法(如窗函数法、频率抽样法、最小二乘法等),根据滤波器的长度和频率特性要求,设计出滤波器的冲激响应。
4.计算滤波器的频率响应:将设计得到的滤波器的冲激响应进行傅里叶变换,得到滤波器的频率响应。
可以使用FFT算法来进行计算。
5.优化滤波器的性能:根据频率响应的实际情况,对滤波器的冲激响应进行优化,可以通过调整滤波器的系数或使用优化算法来实现。
6.实现滤波器:将设计得到的滤波器的冲激响应转化为差分方程或直接形式,并使用数字信号处理器(DSP)或其他硬件进行实现。
7.验证滤波器的性能:使用测试信号输入滤波器,检查输出信号是否满足设计要求,并对滤波器的性能进行验证和调整。
以上是FIR数字滤波器的一般设计步骤,具体的设计方法和步骤可能因应用需求和设计工具的不同而有所差异。
在实际设计中,还需要考虑滤波器的实时性、计算复杂度和存储资源等方面的限制。
信息工程学院课程设计报告课程:数字信号处理学号:2009012347 2009012387姓名:班级:通信0903 教师:日期:2012-5-22摘要:随着计算机和信息科学的飞速发展,数字信号处理逐渐发展成为一门独立的学科,成为信息科学的重要组成部分。
在语音处理、雷达、图形处理、通信、生物医学工程等众多的领域中得到广泛的应用。
下面主要介绍数字信号处理的使用软件 MATLAB 以及使进行语音信号的分析处理、应用。
MATLAB 在信号与系统中的应用主要包括符号运算和数值计算仿真分析。
由于信号与系统课程的许多内容都是基于公式演算,而借助符号数学工具箱提供的符号运算功能基本满足信号与系统课程的需求。
MATLAB主要包括五大通用的功能:数值计算功能(Nemeric)符号运算功能,数据可视化功能(Graphic);数据图形文字统一处理功能(Notebook);建模仿真可视化功能(Simulink)本次数字信号处理课程设计采用Matlab和建模仿真可视化功能(Simulink)本次数字信号处理课程设计采用Matlab软件对数。
例如,解微分方程、傅里叶正反变换、拉普拉斯正反变换、z正反变换等。
MATLAB 在信号与系统中的另一主要应用是数值计算与仿真分析,主要包括函数波形绘制、函数运算、冲激响应与阶跃响应仿真分析、信号的时域分析、信号的频谱分域分析、零极点图绘制等内容。
数值计算仿真分析可以帮助我们能更深入分析系统S域分析、零极点图绘制等内容。
并为将来使用 MATLAB 进行信号处理领域的各种分析和实际应用打下基础设计名称:基于频率采样法FIR数字滤波器的设计设计内容要求:采取最优化的设计方法通过频率采样实现FIR滤波器,这种方法是为了得到更大的衰减,增大M值,并使过渡带中的样本成为自由样本,这个问题可以采用线性规划技术来解决。
在下面的例子中,仍设低通滤波器指标为:0.2,0.250.3,50p p s p R dBR dB ωπωπ====设计原理与说明:FIR 滤波器设计方法:频率采样法的基本思想是使所设计的FIR 数字滤波器的频率特性在某些离散频率点上的值准确地等于所需滤波器在这些频率点处的值,在其它频率处的特性则有较好的逼近。
基于频率抽样法的FIR 数字低通滤波器的设计1 设计目的熟悉频率采样法的理论及其应用;掌握频率采样法设计FIR 数字滤波器的方法。
了解FIR 数字滤波器的频率特性和相位特性,观察过渡带取样点对滤波器幅频特性的影响。
掌握用频率采样法设计线性相位FIR 低通数字滤波器的方法,并掌握该方法的matlab 编程和仿真。
2 FIR 数字滤波器设计的原理2.1频率抽样设计法FIR 低通滤波器的设计一般方法有两种,即频率抽样法和窗函数法,频率抽样法设计不同于窗函数法,窗函数是从时域出发,把理想的()d h n 用一定形状得窗函数截取成有限长的()h n ,以此()h n 来近似理想的()d h n ,这样得到的频率响应()jw H e 逼近于所要求的理想的频率响应()jw d H e 。
频率抽样法则是从频域出发,把给定的理想频率响应()jw d H e 加以等间隔抽样,即2()|()jw d d w k N H e H k π==然后以此()d H k 作为实际FIR 数字滤波器的频率特性的抽样值()H k ,即令2()()()|0,1,,1jw d d w k N H k H k H e k N π====-,知道()H k 后,由DFT 定义,可以用频域的这N 个抽样值()H k 来唯一确定有限长序列()h n ,而由()X z 的内插公式知道,利用这N 个频域抽样值()H k 同样可求得FIR 滤波器的系统函数()H z 及频率响应()jw H e 。
这个()H z 或()jw H e 将逼近()d H z 或()jw d H e ,()H z 和()jw H e 的内插公式为1101()()1NN k k N z H k H z NW z ----=-=-∑ (2.1) 102()()()N jwk H e H k w k Nπ-==Φ-∑ (2.2) 其中()w Φ是内插函数1()2sin()12()sin()2N jw wN w e w N --Φ= (2.3) 将式(2.3)代入(2.2)式,化简后可得11()20sin()12()()sin()2k N N j j w jw N k wN H e e H k e w k N Nππ----==-∑ (2.4) 即 11(1)()20sin[()]12()()sin()2k N N j N j w jw N k w k N N H e e H k e w k N N πππ----=-=••-∑ (2.5) 从内插公式(2.2)看到,在各频率抽样点上,滤波器的实际频率响应是严格地和理想频率响应数值相等,即22()()()()j k j k N N d d H e H k H k H e ππ===。
FIR低通滤波器的设计低通滤波器是一种常见的信号处理工具,它可以将高频信号从输入信号中滤除,只留下低频信号。
在很多应用中,低通滤波器被用于去除噪声、平滑信号、降低带宽等。
设计一个低通滤波器需要考虑多个因素,包括滤波器类型、阶数、截止频率、群延迟等。
以下是设计低通滤波器的步骤:1. 确定滤波器类型:首先需要选择滤波器的类型,常见的低通滤波器有巴特沃斯(Butterworth)、切比雪夫(Chebyshev)、椭圆(Elliptic)等。
每种类型的滤波器有不同的特性和设计参数,选择合适的类型取决于具体的应用需求。
2.确定滤波器阶数:滤波器的阶数与其滤波特性的平滑程度有关,阶数越高,滤波曲线越陡峭。
一般来说,阶数越高,滤波器设计越复杂,实现难度也越大。
选择适当的阶数需要在设计要求和性能之间进行平衡。
3.确定截止频率:截止频率是指滤波器在此频率以下开始滤除高频信号的频率。
确定截止频率需要考虑到信号中的有用频率范围以及滤波器对信号的影响。
需要注意的是,低通滤波器的截止频率应该小于采样频率的一半,否则会导致混叠效应。
4.根据以上参数进行滤波器设计:根据选择的滤波器类型、阶数和截止频率,可以利用不同的设计方法进行滤波器设计。
常用的设计方法有频率变换法、零极点法、传递函数设计法等。
这些方法可以通过数学计算或者使用相关软件进行设计。
5.实现滤波器:设计好滤波器后,需要将其实现到具体的系统中。
这通常涉及到电子电路、数字信号处理器(DSP)或者软件实现。
具体的实现方式取决于应用要求和所使用的平台。
在设计低通滤波器时,还需要考虑一些其他因素,例如群延迟、通带波动、阻带抑制等。
群延迟是指滤波器对不同频率的信号引起的延迟差异,通常希望群延迟尽可能平均,以避免引起相移问题。
通带波动是指滤波器在通带内的幅频响应变化情况,阻带抑制是指滤波器在阻带内对高频信号的抑制能力。
总结来说,低通滤波器设计是一个综合考虑信号需求、滤波器特性和实现条件的过程。
F I R数字低通滤波器设计(总14页)--本页仅作为文档封面,使用时请直接删除即可----内页可以根据需求调整合适字体及大小--第1章 绪论设计的作用、目的课程设计是理论学习的延伸,是掌握所学知识的一种重要手段,对于贯彻理论联系实际、提高学习质量、塑造自身能力等于有特殊作用。
本次课程设计一方面通过MATLAB 仿真设计内容,使我们加深对理论知识的理解,同时增强其逻辑思维能力,另一方面对课堂所学理论知识作一个总结和补充。
设计任务及要求通过课程设计各环节的实践,应使学生达到如下要求:1.掌握双线性变换法及脉冲响应不变法设计IIR 数字滤波器以及窗函数法 设计FIR 数字滤波器的原理、具体方法及计算机编程。
2.观察双线性变换法、脉冲响应不变法及窗函数法设计的滤波器的频域特性,了解各种方法的特点。
3.用MATLAB 画出三种方法设计数字滤波器的幅频特性曲线,记带宽和衰减量,检查结果是否满足要求。
设计内容设计题目:FIR 数字滤波器的设计 设计内容:(1)设计一线性相位FIR 数字低通滤波器,截止频率π2.0=Ωf ,过渡带宽度 π4.0≤∆Ω,阻带衰减dB A s 30>。
(2)设计一线性相位FIR 数字低通滤波器,截止频率π2.0=Ωf ,过渡带宽度π4.0≤∆Ω,阻带衰减dB A s 50>。
第2章FIR 数字低通滤波器的原理数字低通滤波器的设计原理FIR 数字滤波器传统的设计方法有窗函数法、频率抽样法和等波纹逼近法。
用窗函数设计FIR 数字滤波器就是用有限长的脉冲相应逼近序列,其基本设计思想为:首先选定一个理想的选频滤波器,然后截取它的脉冲响应得到线性相位。
滤波器(filter ),是一种用来消除干扰杂讯的器件,将输入或输出经过过滤而得到纯净的直流电。
对特定频率的频点或该频点以外的频率进行有效滤除的电路,就是滤波器,其功能就是得到一个特定频率或消除一个特定频率。
就是允许某一部分频率的信号顺利的通过,而另外一部分频率的信号则受到较大的抑制,它实质上是一个选频电路。
基于频率抽样法的低通FIR滤波器的设计频率抽样法是一种常见的低通FIR滤波器设计方法。
FIR滤波器是一种无反馈的数字滤波器,通过对输入信号的每个采样点进行加权平均来输出滤波后的信号。
在频率抽样法中,滤波器的频率响应是通过在频率域进行采样和插值得到的。
首先,我们需要确定滤波器的截止频率。
截止频率是指滤波器在此频率以上的信号会被抑制,而在此频率以下的信号会被保留。
截止频率的选择要根据实际应用需求进行。
通常情况下,截止频率的选择是基于信号的带宽及噪声的影响。
接下来,我们要确定滤波器的阶数。
阶数是指FIR滤波器中加权延迟线的数量。
阶数的选择会影响滤波器的截止带宽、过渡带宽和频率响应的幅频特性。
然后,我们需要确定频率抽样点的数量。
频率抽样点的选择会影响滤波器的频率响应的分辨率。
通常情况下,频率抽样点的数量越多,分辨率越高,但是计算复杂度也会增加。
在确定了上述参数后,我们可以按照以下步骤进行基于频率抽样法的低通FIR滤波器的设计:1.在频率域中定义一个理想低通滤波器的频率响应,该频率响应在截止频率以下的频率范围内具有幅度为1的增益,而在截止频率以上的频率范围内具有幅度为0的增益。
2.计算频率抽样点的间隔,即频率抽样点之间的频率间隔。
频率抽样点的数量可以根据需要进行选择,一般选择足够大的数值来保证滤波器的精度和性能。
3.在频率域中选择频率抽样点,这些频率点应该尽量均匀地分布在滤波器的所需频率范围内。
4.对于每个频率抽样点,计算出其对应的滤波器的幅度响应。
幅度响应可以通过将理想低通滤波器的频率响应与频率抽样点进行卷积求得。
5.按照频率抽样点的顺序,将得到的幅度响应值转换为时域中的滤波器的冲激响应。
可以使用逆傅里叶变换来实现这一转换。
在计算逆傅里叶变换时,需要使用快速傅里叶变换算法来加快计算速度。
6.最后,根据所得到的滤波器的冲激响应,可以通过对输入信号的每个采样点进行加权平均来实现滤波器的操作。
频率取样法设计fir滤波器频率取样法是一种设计FIR滤波器的方法,它的基本思想是将滤波器的频率响应看作是离散时间傅里叶变换(DTFT)的样本,然后通过对这些样本进行插值来得到滤波器的系数。
在这个过程中,采样频率和滤波器阶数是两个重要的参数。
首先,确定采样频率。
采样频率应该高于信号中最高频率成分的两倍,以避免混叠现象。
通常情况下,选择一个稍微高于两倍最高频率成分的采样频率即可。
其次,确定滤波器阶数。
滤波器阶数决定了滤波器能够抑制多少不需要的频率成分。
一般情况下,阶数越高,抑制能力越强,但也会导致计算量增加和延迟增加。
因此需要在抑制不需要的频率成分和保留需要的频率成分之间做出权衡。
接着,在选定采样频率和滤波器阶数后,可以开始进行滤波器设计。
具体步骤如下:1. 确定截止频率:根据实际需求确定所需截止频率,并将其归一化到采样频率的一半。
2. 确定通带和阻带:根据截止频率确定通带和阻带的边界。
3. 确定通带和阻带增益:根据实际需求确定通带和阻带的增益,以及过渡区域的斜率。
4. 计算DTFT样本:根据上述参数计算DTFT样本,即在通带内取样一些点,然后通过插值得到整个DTFT。
常用的插值方法有线性插值、三次样条插值等。
5. 计算滤波器系数:通过对DTFT进行反离散时间傅里叶变换(IDTFT),可以得到滤波器系数。
需要注意的是,IDTFT是一个复杂的计算过程,可以利用FFT快速算法来加速计算。
6. 实现滤波器:将得到的滤波器系数应用于差分方程中,就可以实现FIR滤波器了。
最后需要注意的是,在实际应用中,由于离散时间傅里叶变换具有周期性,因此在设计FIR滤波器时需要使用周期延拓技术来避免边界效应。
常用的周期延拓方式有零填充、循环延拓等。
总之,频率取样法是一种简单而有效的FIR滤波器设计方法,它可以根据实际需求灵活地选择采样频率和滤波器阶数,并通过插值和反离散时间傅里叶变换得到滤波器系数。
在实际应用中,需要注意周期延拓技术来避免边界效应。
基于频率抽样法的FIR 数字低通滤波器的设计1 设计目的熟悉频率采样法的理论及其应用;掌握频率采样法设计FIR 数字滤波器的方法。
了解FIR 数字滤波器的频率特性和相位特性,观察过渡带取样点对滤波器幅频特性的影响。
掌握用频率采样法设计线性相位FIR 低通数字滤波器的方法,并掌握该方法的matlab 编程和仿真。
2 FIR 数字滤波器设计的原理2.1频率抽样设计法FIR 低通滤波器的设计一般方法有两种,即频率抽样法和窗函数法,频率抽样法设计不同于窗函数法,窗函数是从时域出发,把理想的()d h n 用一定形状得窗函数截取成有限长的()h n ,以此()h n 来近似理想的()d h n ,这样得到的频率响应()jw H e 逼近于所要求的理想的频率响应()jw d H e 。
频率抽样法则是从频域出发,把给定的理想频率响应()jw d H e 加以等间隔抽样,即2()|()jw d d w k NH e H k π==然后以此()d H k 作为实际FIR 数字滤波器的频率特性的抽样值()H k ,即令2()()()|0,1,,1jw d d w k NH k H k H e k N π====-,知道()H k 后,由DFT 定义,可以用频域的这N 个抽样值()H k 来唯一确定有限长序列()h n ,而由()X z 的内插公式知道,利用这N 个频域抽样值()H k 同样可求得FIR 滤波器的系统函数()H z 及频率响应()jw H e 。
这个()H z 或()jw H e 将逼近()d H z 或()jw d H e ,()H z 和()jw H e 的内插公式为1101()()1NN k k Nz H k H z NW z ----=-=-∑ (2.1) 102()()()N jwk H e H k w k Nπ-==Φ-∑ (2.2) 其中()w Φ是内插函数1()2sin()12()sin()2N jw wNw e w N --Φ= (2.3)将式(2.3)代入(2.2)式,化简后可得11()2sin()12()()sin()2kN N j j wjw N k wNH e e H k e w k N N ππ----==-∑ (2.4) 即 11(1)()20sin[()]12()()sin()2k N N j N j w jw N k w k N N H e e H k e w k N Nπππ----=-=••-∑ (2.5) 从内插公式(2.2)看到,在各频率抽样点上,滤波器的实际频率响应是严格地和理想频率响应数值相等,即22()()()()j k jk NNd d H eH k H k H eππ===。
但是在抽样点之间的频率响应则是由各抽样点的加权内插函数的延伸叠加而形成的,因而有一定的逼近误差,误差大小取决于理想频率响应曲线形状,理想频率响应特性变化越平缓,则内插值越接近理想值,逼近误差越小,如图2.1梯形理想频率特性所示。
反之,如果抽样点之间的理想频率特性变化越陡,则内插值与理想值之误差就越大,因而在理想频率特性的不连续点附近,就会产生肩峰和波纹,如图2.2矩形理想频率特性所示。
图2.1 梯形理想频率特性图2.2 矩形理想频率特性2.2线性相位的约束对于第二类线性相位FIR 滤波器,由于()h n 偶对称、N 为偶数时, 公式是1()2()()N j w jwH e H w e--=,其中幅度函数()H w 应为奇对称的,()(2)H w H w π=--,如果抽样值2()()jk NH k H eπ=也用幅值k H (纯标量)与相角k Q 表示,则为:22()()()K K j k jQ jQ Nk H k H e H k e H e Nππ=== (2.6)其中k Q 必须为: 121()(1)2k N Q k k N Nππ-=-=-- (2.7)k H 必须满足奇对称,即k N k H H -=-。
2.3线性相位第一种频率抽样21()()N jnk Nn H k h n eπ--==∑ (2.8)当()h n 为实数时,满足**()(())()()N N H k H N k R k H N k =-=-,由此得出()()H k H N k =-,()()Q k Q N k =--,也就是说,()H k 的模()H k 以2k N =为对称中心呈偶对称,()H k 的相角()Q k 以2k N =为对称中心呈奇对称。
再利用线性相位的条件1()2jw N Q e w -=-,即可得到(N 为偶数): 210,,122()0221()1,,122N N k k N NQ k k N NN k k N Nππ⎧-⎛⎫⎛⎫-=- ⎪⎪⎪⎝⎭⎝⎭⎪⎪==⎨⎪⎪-⎛⎫⎛⎫-=+- ⎪ ⎪⎪⎝⎭⎝⎭⎩(2.9)21()221()()2()0,,12()02()1,,12N j k N N j N k N N H k e k NH k k N H N k e k N ππ----⎧⎛⎫=-⎪ ⎪⎝⎭⎪⎪==⎨⎪⎪⎛⎫-=+-⎪ ⎪⎝⎭⎩(2.10)()12(0)sin 2sin 2N j j N H H e e N ωωωω--⎧⎛⎫ ⎪⎪⎪⎝⎭=⎨⎛⎫⎪ ⎪⎪⎝⎭⎩()121sin sin 22sin sin 22Nk k k N N H k N N k k N N N ωπωπωπωπ-=⎫⎡⎤⎡⎤⎡⎤⎛⎫⎛⎫-+⎪ ⎪ ⎪⎢⎥⎢⎥⎢⎥⎪⎝⎭⎝⎭⎣⎦⎣⎦⎢⎥++⎬⎢⎥⎛⎫⎛⎫⎪-+ ⎪ ⎪⎢⎥⎪⎝⎭⎝⎭⎣⎦⎭∑(2.11)2.4过渡带抽样的优化设计为了提高逼近质量,使逼近误差更小,也就是减小在通带边缘由于抽样点的陡然变化而引起的起伏振荡,这种起伏振荡使阻带内最小衰减变小。
这里是使理想频率响应的不连续点的边缘加上一些过渡的抽样点,从而增加过渡带,减小频率边缘的突变,也就减小了起伏振荡,增大了阻带最小衰减。
这些周扬点上的取值不同,效果也就不同,如果精心设计过渡带的抽样值,就有可能使它的有用频带(通带、阻带)的纹波得以减小。
从而设计出较好的滤波器。
一般过渡带取一、二、三点抽样值即可得到满意的结果。
在低通设计中,不加过渡抽样点时,阻带最小衰减为-20dB,一点过渡抽样的最优设计,阻带最小衰减可提高到-40dB到-54dB左右,二点过渡抽样的最优设计可达-60dB 到-75dB左右,而加三点过渡抽样的最优设计则可达-80dB到-95dB左右。
加过渡抽样点的示意图如图2.3所示。
图2.3 加过渡抽样点3 FIR 数字滤波器设计方法3.1设计任务本次课程设计的任务是利用频率抽样法设计一个低通FIR 数字低通滤波器,其理想频率特性是矩形的,即(3.1)给定抽样频率为42 1.510(/sec)s rad πΩ=⨯⨯,通带截止频率的32 1.610(/sec)p rad πΩ=⨯⨯,阻带起始频率为32 3.110(/sec)st rad πΩ=⨯⨯,通带波动11dB σ≤,阻带衰减250dB σ≥。
幅度特性如图3.1所示。
图3.1 ()H j Ω幅度特性曲线3.2参数计算通带的截止频率为342 1.610220.2132 1.510p pp s s w f πππππΩΩ⨯⨯===⨯=Ω⨯⨯ 阻带的起始频率为342 3.110220.4132 1.510st st st s s w f πππππΩΩ⨯⨯===⨯=Ω⨯⨯理想低通截止频率31()2 2.3510(/sec)2c p st rad πΩ=Ω+Ω=⨯⨯其对应的数字频率342 2.3510220.3132 1.510c c w s πππππΩ⨯⨯==⨯=Ω⨯⨯ 过渡带带宽为0.2st p w w w π∆=-=,设抽样点数为N ,由于23w Nπ∆=⨯,故算得抽样点数N 等于30。
3.3设计方法根据指标,可画出频率抽样后的()H k 序列,由于()H k 是对称于w π=的,我们又只对0w π≤≤即015k ≤≤的区间感兴趣,故可将2w ππ≤≤即1630k ≤≤的图略去不画。
截止频率0.313c w π=,截止频率抽样点的位置应为:0.313302 4.6954ππ⨯=≈,按第一种频率抽样方式来设计,N=30,则10Int 42()0Int 151522c c N k H k N N k ωπωπ⎧⎡⎤≤≤=⎪⎢⎥⎪⎣⎦=⎨⎡⎤⎪+=≤≤=⎢⎥⎪⎣⎦⎩(3.2) []Int •表示取整数部分,将这些值代入式(2-11)可得:()515130sin 30sin 30sin 230230230sin 30sin 30sin 2230230j j k k k H e e k k ωωωπωπωωωπωπ-=⎧⎫⎡⎤⎡⎤⎡⎤⎛⎫⎛⎫⎛⎫-+⎪⎪ ⎪ ⎪⎢⎥ ⎪⎢⎥⎢⎥⎪⎪⎝⎭⎝⎭⎝⎭⎣⎦⎣⎦⎢⎥=++⎨⎬⎢⎥⎛⎫⎛⎫⎛⎫⎪⎪-+ ⎪ ⎪ ⎪⎢⎥⎪⎪⎝⎭⎝⎭⎝⎭⎣⎦⎩⎭∑ (3.3)按此式计算1020log ()jw H e 的结果如图3.2所示。
由图看出,过渡带宽为230π,而最小阻带衰减约为-20dB 。
这一衰减在大多数情况下是不令人满意的,也不符合这次课程设计的要求.图3.2 幅度频率特性曲线为了改善频率特性,以满足指标要求,可在通带和阻带交界处安排一个或者几个不等于1的抽样值。
在本次课设中用优化算法使K =5处值为0.5886,即(5)0.5886H =,和K =6处值为0.1065,即(6)0.1065H =。
则得到图3.3所示的结果,过渡带带宽为230.230ππ⨯=,最小阻带衰减约为-60dB 左右,如图3.4所示。
H k图3.3 增加过渡点后的()图3.4 过渡带抽样优化后的幅频特性综上所述,采用第一种频率抽样方式来设计,抽样点数取N=30,优化抽样H=和过渡带采用在通带和阻带交界处安排两个不等于1的点,即(5)0.5886H=,最后设计出来的基于频率抽样法FIR数字低通滤波器满足设计(6)0.1065要求。
4 程序设计4.1 两种频率抽样方法幅频特性对比程序close all;clear;N=30;H=[ones(1,4),zeros(1,22),ones(1,4)];k=0:(N/2-1);k1=(N/2+1):(N-1);k2=0;A=[exp(-j*pi*k*(N-1)/N),exp(-j*pi*k2*(N-1)/N),exp(j*pi*(N-k1)*(N-1)/N)];HK=H.*A;hn=ifft(HK);fs=15000;[c,f3]=freqz(hn,1);f3=f3/pi*fs/2;figure(1)subplot(2,1,1)plot(f3,20*log10(abs(c)));title('频谱特性');xlabel('频率/HZ');ylabel('衰减/dB');grid;N=30;H=[ones(1,4),zeros(1,22),ones(1,4)];H(1,5)=0.5886;H(1,26)=0.5886;H(1,6)=0.1065;H(1,25)=0.1065;k=0:(N/2-1);k1=(N/2+1):(N-1);k2=0;A=[exp(-j*pi*k*(N-1)/N),exp(-j*pi*k2*(N-1)/N),exp(j*pi*(N-k1)*(N-1)/N)];HK=H.*A;hn=ifft(HK);fs=15000;[c,f3]=freqz(hn,1);f3=f3/pi*fs/2;subplot(2,1,2)plot(f3,20*log10(abs(c)));title('频谱特性');xlabel('频率/HZ');ylabel('衰减/dB');grid;4.2 FIR数字低通滤波器程序close all;clear;N=30;H=[ones(1,4),zeros(1,22),ones(1,4)];H(1,5)=0.5886;H(1,26)=0.5886;H(1,6)=0.1065;H(1,25)=0.1065;k=0:(N/2-1);k1=(N/2+1):(N-1);k2=0;A=[exp(-j*pi*k*(N-1)/N),exp(-j*pi*k2*(N-1)/N),exp(j*pi*(N-k1)*(N-1)/N)]; HK=H.*A;hn=ifft(HK);fs=15000;[c,f3]=freqz(hn,1);f3=f3/pi*fs/2;figure(1);plot(f3,20*log10(abs(c)));title('频谱特性');xlabel('频率/HZ');ylabel('衰减/dB');grid;figure(2);stem(real(hn),'.');line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');t=(0:100)/fs;W=sin(2*pi*t*750)+sin(2*pi*t*3000)+sin(2*pi*t*6500);q=filter(hn,1,W);[a,f1]=freqz(W);f1=f1/pi*fs/2;[b,f2]=freqz(q);f2=f2/pi*fs/2;figure(4);plot(f1,abs(a));title('输入波形频谱图');xlabel('频率');ylabel('幅度')grid;figure(5);plot(f2,abs(b));title('输出波形频谱图');xlabel('频率');ylabel('幅度')grid;5 仿真结果及分析分析:程序运行后结果如下,图5.1显示的是两种抽样设计方法幅频响应特性对比图,由图可以明显看到,过渡带抽样优化后的频率特性明显优于第一种普通抽样方法。