任意噪声和约束下的最佳数字滤波器设计

2008-05-07 14:33:20编辑:孙树宾 关键字:FIR滤波器  噪声  线性组合  电路参数  时间常数  采样系统  采样周期  指数衰减

  摘 要:为提高数字谱仪的能量分辨率,利用约束最优化方法提出了一个在任意平稳噪声背景下直接设计最佳数字有限冲击响应(FIR)滤波器的方法。该方法易于附加任意时域和频域约束(平顶、基线恢复等)。详细给出了最佳滤波器的设计推导过程,并在设计实例基础上,分析了滤波器的时域、频域响应和不同约束条件下的最大信噪比曲线。结果证明,该方法能够有效去除基线漂移和弹道亏损,得到了约束条件下的最大信噪比,适合在未知系统噪声特性情况下的高计数率、高精度谱仪中使用。

  关键调:最佳滤波器;数字谱仪;基线恢复

  在数字谱仪中,探测器输出信号直接进行模数转换,经过数字滤波可以改善模数转换器(ADC)的线性和精度,因此数字滤波器的设计成为提高整个系统性能指标的关键。谱仪滤波器的设计,除尽可能提高信噪比外,通常还要考虑成形脉冲形状的要求,例如窄成形脉冲以减少堆积,渡形顶部平坦以减少弹道亏损等。

  文提出了一种可以附加任意约束、在任意噪声背景下设计最佳模拟滤波器的方法,由于设计出来的是权函数,无法直接用于数字滤波,而且通过抽样最佳模拟滤波器得到的数字滤波器不是最优滤波器。文提出了移动窗口解卷积的设计方法,此方法无法考虑噪声影响,仅适合指数衰减形式的脉冲。文提出了惩罚最小方差方法设计最佳滤波器,其约束必须写成二次型的形式,添加约束比较困难,并且在多个约束的条件下,无法综合出严格的平顶。

  有限冲击响应(FIR)滤波器没有反馈通路,某个时刻的输出直接由此时刻之前固定时间长度的输入决定,其结构特点易于通用数字信号处理器(DSP)芯片或在线可编程逻辑阵列(FPGA)实现。本文提出了一种基于信噪比准则的直接设计最佳数字FIR滤波器的方法。

  l 问题的提出

  假定噪声是和输入信号不相关的平稳随机过程,在已知输入信号形状噪声自相关函数的条件下,寻找满足约束条件的FIR滤波器,使输出信号的信噪比最大。先求得以滤波器系数为变量的信噪比表达式,然后添加适当约束,把求解最佳滤波器的任务转化成求解一个约束最优化问题。

  1.1 SNR的计算

  设噪声的自相关函数为{r(i)},待求的N阶FIR滤波器的冲击响应为h={h(0),h(1),…,h(N一1)}T,采样得到的信号为{s(i)},由已知波形的输入脉冲信号{x(i)}和噪声信号{n(i)}相加而成。从{x(i)}选取N个连续数据,使N个数的平方和最大,重新排列标号,使得选取的数据为{x(1),…,x(N)}。

  

  

  1.2 时域约束的添加

  滤波后的脉冲为输入和FIR滤波器的卷积,即成形后脉冲各点是输入以h为系数的线性组合,这样把在时域的约束转化成了对h的约束。以添加平顶约束为例,设要求输出脉冲的平顶宽度为m,即要求y(N)直到y(N一m)相等,则有:

  

  

  1.3 频域约束的添加

  FIR的频率响应为

  式中DTFT表示离散Fourier变换。

  频率响应函数为FIR滤波器系数的线性组合,因此把在频域的约束也转化成对h的约束。为了去除基线漂移,FIR的频率响应在0 Hz处为0,即

  平顶和去除基线这两个约束构成了一个线性方程组。这些约束相互之间是独立的,否则可以从这些方程中选取出最大无关向量组来构成A。综合式(1)和(2),构成一个线性约束最优化问题。

  2 问题的求解

  在有约束情况下,约束线性方程组把N维线性空间直和分解分成两个空间:滤波器允许空间(任意向量都满足约束条件)及此空间的正交补。把无约束下的解向滤波器允许空间投影就得到了约束条件下的最佳解。

  2.1 无约束的情况

  噪声的相关矩阵为实对称正定矩阵,可以共轭对角化为

  

  

  2.2 有约束的情况

  在有约束的情况下,h1受约束条件的限制,约束条件转化成一个线性方程组,即

  

  假设共有m个独立约束,则约束方程组的解构成一个N一m维线性空间,称为滤波器允许空间,记做V,则V垂直于矩阵的行空间。把作QR分解得到

  由于P为单位正交阵,P的前m列构成的行空间一个正交基底,后N一m列构成V的一个正交基底。把分解成V中的向量和与V垂直的向量,有

  

  

 

  当约束方程中的b不为0时,则h1被限制在一个仿射空间内,把h1平移到此空间中的任意一点,就变化成和b为0一样的解法。在约束为非线性时,可以采用数值最优化方法求解。

  2.3 求解结果的分析

  无约束的情况:在白噪声情况下,自相关矩阵为噪声功率乘以单位阵,此时的解简化为h=kx。也就是说,滤波器与已知信号的形状成正比,对应于模拟情况下的匹配滤波器。在有色噪声情况下,相当于x先对输入噪声进行白化处理,得到白噪声,因此,有色噪声情况下的最佳滤波器是白化滤波器与其后的白噪声匹配滤波器的级联。

  在有约束的情况下,给定的约束条件限制了滤波器的范围,最佳滤波器为无约束条件下的最佳滤波器向滤波器允许空间的投影。

  2.4 输波形和噪声自相关函数的计算

  在已知前放和抗混叠滤波器参数的情况下,可以直接计算得到输入信号。由于RC参数的不精确,计算得到的波形和实际的波形有差别,可以通过估计极零点的方法加以修正。或者平均多次ADC的采样来减小噪声,从而得到输入信号的波形。

  在已知噪声功率谱的情况下,可以通过Fourier反变换得到连续时间下的自相关函数。在很多情况下,可能不知道准确的系统噪声特性,可以通过采样无信号输入时ADC的输出来计算系统的自相关函数。如果采样得到的噪声序列是{x(i)},那么怙计自相关函数为

  3 设计结果

  以下面的参数为例进行分析:输入脉冲有3个相等的实极点,时间常数为200 ns,幅度归一化为1,采样频率为50 MHz,噪声为白噪声。

  3.1 滤波器长度的选取

  滤波器的长度越长,系统能够得到的信噪比越大。但成形后的脉冲宽度也越宽,因而容易产生堆积。在设计实际的滤波器时,要根据计数率和噪声等因素综合考虑,选择合适的长度。一般在选取输入波形数据段时,如果段数据平方包含的面积和整个波形平方后的总面积相差不多,则设计出的最佳滤波器的信噪比就很接近理论上的最大值。在下面讨论中,把实际滤波器的信噪比以滤波器无限长时的最大信噪比进行了归一化。在上述输入条件下,FIR滤波器长度取70时的信噪比已达到0.998。

  3.2 设计结果

  取滤波器的长度为70,在无约束、有平顶宽度为20、去除基线和既有平顶叉去除基线4种约束情况下分别求解,FIR滤波器的设计结果如图l所示。有平顶约束的FIR滤波器其尾部出现振荡,振荡的长度基本上就是平顶约束的长度,在附加去除基线漂移的约束下,滤波器向下移动,总面积为0。

 

  

 

  

  4种不同约束条件下滤波成形后的脉冲波形如图2所示,采用此设计方法成形后的脉冲有很平坦的顶部,可以减小弹道亏损。采用FIR滤波器,成形后的脉冲宽度是滤波器的至少2倍宽,有一定的前后沿,但由于每个输出采样点的数据只是前面一定长度的数据的卷积,输入脉冲的时间间隔如果大于滤波器的宽度,就不会发生峰堆积。

 

  

 

  

  图3为3种约束情况下滤波器的频率响应。附加去除基线约束后,可以看到,在极低频部分有很大的衰减,对直流的抑制达到400 dB,但输入信号的能量很大一部分也在低频部分,成形后的脉冲幅度也降低很多,导致信噪比变小。

 

 

  

 

 

  

  3.3 约束对信噪比的影响

  约束条件下得到的最大信噪比要小于无约束的情况,图4是不同滤波器长度N在各种约束条件下的信噪比曲线。从频率响应图可以看出,在附加平顶约束后,高频部分的振荡幅度加大,对高频噪声的抑制降低,但在中低频它和无约束频率响应曲线基本一致,信噪比损失比基线约束小。根据是否有基线约束,图4上的SNR曲线分成了明显的两组,附加基线约束后的信噪比要降低很多。无基线约束的条件下,在滤波器长度N=70时,SNR对N的增加已不敏感,N>70的3条曲线几乎重合。

  

  

  实际应用中,选取的平顶宽度大于电荷收集时间一个采样周期可以避免弹道亏损。是否去除基线可根据系统的需要来选择,在低计数率的条件下,采用单独的基线估计器可能会有更好的信噪比。

  4 结 论

   本文提出了在约束条件下直接设计最佳数字滤波器的方法。整个求解过程稳定,可在很短的时间内完成,适合实时应用。通过添加不同约束,可以方便地设定成形脉冲平坦顶部宽度;设计出的滤波器可以直接去除基线漂移,适合在系统温度、电路参数变化大的环境中使用。还可通过采样系统噪声估算噪声的自相关函数,综合出最佳滤波器,实现自动配置,减少人工调节,适合在高计数率的高精度谱仪中应用。

 

关键字:FIR滤波器  噪声  线性组合  电路参数  时间常数  采样系统  采样周期  指数衰减

来源: 老古开发网 引用地址:http://www.eeworld.com.cn/afdz/2008/0507/article_553.html
本网站转载的所有的文章、图片、音频视频文件等资料的版权归版权所有人所有,本站采用的非本站原创文章及图片等内容无法一一联系确认版权者。如果本网所选内容的文章作者及编辑认为其作品不宜公开自由传播,或不应无偿使用,请及时通过电子邮件或电话通知我们,以迅速采取适当措施,避免给双方造成不必要的经济损失。

上一篇:一种隔离型双向软开关DC/DC变换器
下一篇:基于共享技术的电荷泵电路

关注eeworld公众号 快捷获取更多信息
关注eeworld公众号
快捷获取更多信息
关注eeworld服务号 享受更多官方福利
关注eeworld服务号
享受更多官方福利

推荐阅读

详解FIR滤波器和IIR滤波器的区别

数字滤波器广泛应用于硬件电路设计,在离散系统中尤为常见,一般可以分为FIR滤波器和IIR滤波器,那么他们有什么区别和联系呢。FIR滤波器定义:FIR滤波器是有限长单位冲激响应滤波器,又称为非递归型滤波器,是数字信号处理系统中最基本的元件,它可以在保证任意幅频特性的同时具有严格的线性相频特性,同时其单位抽样响应是有限长的,因而滤波器是稳定的系统。特点:FIR滤波器的最主要的特点是没有反馈回路,稳定性强,故不存在不稳定的问题;FIR具有严格的线性相位,幅度特性随意设置的同时,保证精确的线性相位;FIR设计方式是线性的,硬件容易实现;FIR相对IIR滤波器而言,相同性能指标时,阶次较高,对CPU的性能要去较高。图1 FIR滤波原理图
发表于 2017-05-08 17:32:19
详解FIR滤波器和IIR滤波器的区别

LabVIEW系列——FIR滤波器与IIR滤波器比较

FIR滤波器 FIR——Finite impulse response,有限冲激响应。 FIR滤波器,也可称为非递归滤波器,卷积滤波器,滤波器当前输出与当前输入值和以前输入值有关。公式如下:   x是输入序列,y是滤波后的输出序列,h是FIR滤波器系数。   IIR滤波器 IIR——Infinite impulse response,无限冲激响应。 IIR滤波器,也可称为递归滤波器,滤波器当前输出不仅与当前输入值,以前输入值,还与以前输出值有关。因此IIR滤波器的冲激响应趋向与0,当永远不为0,这是和FIR滤波器的区别之一。 公式如下
发表于 2015-10-16 08:27:17

DSP in FPGA:FIR滤波器(二)

  这一节主要讲解一下转置型FIR滤波器实现。  FIR滤波器的单位冲激响应h(n)可以表示为如下式:    对应转置型结构的FIR滤波器,如图1所示,抽头系数与上一节中讲解直接型FIR滤波器的实例相同,滤波器阶数为10。    图1  可以发现转置型结构不对输入数据寄存,而是对乘累加后的结果寄存,这样关键路径上只有1个乘法和1个加法操作,相比于直接型结构,延时缩短了不少。  综合得到结果如下:  Number of Slice Registers: 1  Number of Slice LUTs: 18  Number of DSP48E1s: 11  Minimum period: 4.854ns{1} (Maximum
发表于 2014-12-27 22:57:13
DSP in FPGA:FIR滤波器(二)

基于FPGA+MATLAB的串行多阶FIR滤波器设计

数字滤波器是用于过滤时间离散信号的数字系统,通过对抽样数据进行数学处理达到频域滤波的目的。根据单位冲激响应函数的时域特性可分为两类:无限冲激响应(Infinite Impulse Response,IIR)滤波器和有限冲激响应(Finite Inpulse Response,FIR)滤波器。与IIR滤波器相比,FIR滤波器的实现是非递归的,较稳定;且FIR滤波器可获得严格的线性相位特性。因此,应用领域较广。FIR滤波器具有成熟的结构,使用Matlab强大的功能使得本来繁重的计算工作变得轻松。在具体应用时,要根据工程当中信号的特点:采样速率、信号带宽等具体使用Matlab得到FIR滤波器系数。滤波器的结构实现可以使用PC软件、单片机
发表于 2014-05-24 10:59:46
基于FPGA+MATLAB的串行多阶FIR滤波器设计

基于DSP的FIR滤波器的设计(二)

初始状态、输入1个和2个样本后的存储情况。   图1 8级循环缓冲区结构。   3 仿真结果   下面是一组信号采样序列样本x (n), 其中存在有高频干扰。现以x (n) 作为输入序列, 然后滤除其中的干扰成分。   {x (n)} = {-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,*,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0}   本设计的线性相位低通FIR数字滤波器的截止频率ωc为0.2
发表于 2013-09-24 10:45:50
基于DSP的FIR滤波器的设计(二)

基于DSP的FIR滤波器的设计(一)

引言   数字信号处理现已在通信与信息系统、信号与信息系统、自动控制、需达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理应用中, 滤波占有十分重要的地位, 如对信号的过滤、检测、预测等, 都要广泛地用到滤波器。IIR数字滤波器的设计保留了一些典型模拟滤波器优良的幅度特性, 但所涉及的滤波器相位特性一般是非线性的, 而FIR滤波器则可在保证幅度特性并满足技术要求的同时, 也很容易做到严格的线性相位特性。   1 基于窗函数法的FIR滤波器设计1.1 单位冲激响应   首先应根据技术要求确定待求滤波器的单位冲激响应hd (n)。如果给出待求滤波器的频率为Hd (ej), 那么单位取样响应则可用下式求出
发表于 2013-09-24 10:45:06
基于DSP的FIR滤波器的设计(一)

小广播

About Us 关于我们 客户服务 联系方式 器件索引 网站地图 最新更新 手机版

站点相关: 视频监控 智能卡 防盗报警 智能管理 处理器 传感器 其他技术 综合资讯 安防论坛

北京市海淀区知春路23号集成电路设计园量子银座1305 电话:(010)82350740 邮编:100191

电子工程世界版权所有 京ICP证060456号 京ICP备10001474号 电信业务审批[2006]字第258号函 京公海网安备110108001534 Copyright © 2005-2018 EEWORLD.com.cn, Inc. All rights reserved