基于labview的数字滤波器仿真研究.docxVIP

基于labview的数字滤波器仿真研究 数字核能谱仪精度高，干扰好，重建能力强，是新一代核材料的发展趋势。在高分辨率的核能谱测量中,对核信号进行滤波成形是一项重要的环节。该环节的任务一是提高信噪比,二是确定信号形状以满足后续多道测量分析的要求。数字梯形滤波作为一种常用的滤波成形方法,能够减少脉冲堆积、适应成形要求,从而有效地替代复杂的模拟滤波成形电路,改善系统的稳定性、灵活性及适应性。 1 数字滤波的正交试验 根据核能谱测量系统中最佳滤波器的理论,最佳滤波器的输出波形应是一种无限宽尖顶脉冲。这种波形具有最佳信噪比,因此测量核脉冲信号( 含噪声的指数衰减信号) 幅度的最好方法是将其滤波成形为无限宽尖顶脉冲,再进行幅度测量及分析。但无限宽尖顶脉冲是一种理想的波形而无法物理实现,因为它具有无限长尾。实际的滤波成形器都采用了其他成形波形。其中三角波具有较好的信噪比和窄脉宽。在数字滤波成形中,将指数衰减信号进行简单的递归处理便可得到三角波,而三角波是梯形波中特殊的一种。三角波的峰值采样要求较高,而平顶宽度可调的梯形波能够很好地解决这个问题,并能够充分适配电荷收集时间变化的半导体探测器。 梯形成形滤波是将输入的指数衰减信号滤波成形为斜坡宽度和平顶宽度可调的梯形波。如图1 所示。呈单指数衰减的输入信号及梯形滤波成形输出信号的离散序列表达式分别由式( 1) 和式( 2) 所示。 梯形滤波器的系统函数可表示为: 2 atla理论 数字梯形滤波器的仿真主要采用了LabVIEW中的控制设计与仿真工具。该工具可以有效地对连续线性、非线性系统及离散系统进行仿真,操作较为简单。用户通过前面板输入仿真信号的参数,包括仿真时间、时间常数、采样周期。输入信号的噪声标准差默认为0. 05,幅值为3V。同时,输入梯形成形参数,包括梯形斜坡宽度和平顶宽度。经公式节点计算出各项系数确定系统函数对输入信号进行数字滤波,然后输出滤波后的信号。为了对滤波器进行一定的分析,该程序使用了Lab VIEW与mat lab的混合化编程,在Lab VIEW中编写mat lab脚本语言可以直接调用其相关的服务。主要通过mat lab计算滤波器的零、极点及幅频响应。仿真前面板如图2 所示。 从梯形滤波器的系统函数表达式可以看出,其系统函数随梯形成形脉冲斜坡宽度、平顶宽度以及时间常数变化而变化。图3 给出了采样周期为0. 05 μs、时间常数为1 μs、平顶宽度为1 μs的滤波器输出随斜坡宽度D的变化。图4 给出了斜坡宽度为1 μs的滤波器输出随平顶宽度D的变化。随着平顶宽度或斜坡宽度的增大,噪声对系统的影响越小。但是取值过大会出现输出脉冲的堆积或重叠。一般情况下梯形平顶宽度与探测器收集电荷的时间有关,因此在满足后续测量要求的情况下适当减小输出平顶宽度可以减小脉冲堆积的概率。 图5 为采样周期0. 05 μs、斜坡宽度1 μs、平顶宽度0. 5 μs的滤波器随时间常数 τ 的变化。其取值偏离输入信号的时间常数时,输出幅度会发生变化,并伴有过零趋势。取值越大,幅度亏损和下冲现象越明显。 频率响应在频域表征系统特性,可反映滤波类型及滤波效果。图6 给出了采样周期为0. 05 μs,平顶宽度为0. 5 μs,时间常数为1 μ时滤波器的幅频响应随ta值的变化。梯形斜坡时间较短时,梯形滤波器的幅频响应类似于积分梳状滤波器的幅频特性,具有不规律的振荡旁瓣,振荡周期较长、过渡带较大。梯形斜坡时间越长,其通带特性越明显,旁瓣振荡周期小、过渡带越窄,表明了滤波器对高频噪声作用范围越广、抑制能力越强,但在一定程度上加大了脉冲堆积的概率。实际应用中可根据噪声特性来确定斜坡宽度。 3 数字视频滤波系统的原理 采用Lab VIEW平台实现了数字梯形滤波器的仿真研究,在一定程度上对滤波器性能进行了有效分析,验证了滤波器参数对系统的影响。在实际应用中,可利用该软件平台观察与分析系统的输入输出信号,对滤波器参数进行配置和优化,进而使用FPGA或DSP等数字信号处理设备来对核信号进行数字滤波。 其中:A为信号幅值;τ为衰减时间常数;T为采样周期;na为经时间离散化的梯形波斜坡宽度(na=ta/T,即梯形斜坡时间区间内的采样点数)。设梯形波平顶宽度为D,则nb=na+D,nc为梯形波底边宽度,nc=na+nb。 利用z域微分定理,对式(1)、(2)进行z变换,可得: 由实数位移定理Z{s(t-k T)}=z-kS(z),对(6)式进行z反变换,可得滤波器的差分方程为: 梯形滤波器的输出到输入存在有反馈,可归结为一种IIR滤波器。其系统函数零点数目为(2ta+tb)/T+1,即(2na+nb+1),且在单位圆上存在有2个极点,可判定梯形滤波系统处于临界稳定,为非因果稳定系统。系统频率响应能反映滤波器的抗噪特性,根据系统函数可