matlab中的梯形成形,基于MATLAB的数字化梯形滤波器的研究

第 28 卷 　第 1 期 核电子学与探测技术 Vol. 28 　No. 1 　 　 2008 年 　 1 月 Nuclear Electronics & Detection Technology Jan. 　 2008 　　 基于 MATLAB的数字化梯形滤波器的研究 周清华 ,张软玉 ,李泰华 (四川大学辐射物理及技术教育部重点实验室 ,四川大学物理科学与技术学院 ,成都 　 610064) 　　摘要 :为方便而快捷地开展数字化成形系统的研究 ,本文建立了一种运用 MATLAB 实现梯形滤波 器的参数优化的方法 ,并用此方法讨论了各参数对滤波效果的影响。 关键词 :MATLAB ;梯形成形 ;数字滤波器 ;数字核信号处理 中图分类号 : 　TN713 　　文献标识码 : 　A 　　文章编号 : 　 025820934(2008) 010085204收稿日期 :200604208作者简介 :周清华 ,四川大学硕士研究生 ,研究方向为测控技术及应用。 　　与模拟核仪器相比 ,数字化核仪器系统具有更好的稳定性、抗干扰性和保真性 ,是新一代核仪器发展的趋势。核信号处理方法是核仪器研究的重要内容 ,梯形滤波器是核信号处理的重要方法之一。本文介绍了数字化核信号处理梯形滤波器分析软件 ,使用该软件可以根据测量需要在对话框中配置各参数并即时显示处理结果 ,实现对滤波器的优化。 1 　数字梯形滤波器的算法设计 数字化梯形成形不仅可以提高核谱仪系统的稳定性 ,还可以改善系统的灵活性和自适应性 ,它能很好地适应核谱仪的脉冲成形要求 ,易于满足不同应用的需要。 梯形滤波算法是将输入指数衰减信号滤波成形为斜坡时间和平顶宽度均可调整的等腰梯形 ,如图 1 所示 ,它在时域的表达式为 : ν 0 = ∑ 4 i =1 yi ( t) y1 = A tu( t) / ta , y2 = - A ( t - ta) u( t - ta) / ta , y3 = - A (t - tb) u( t - tb) / ta , y4 = A ( t - tc) u( t -tc) / ta ,其中 A 为信号脉冲幅度 ,ta 为梯形斜坡 图 1 　梯形脉冲 时间 ,D 为梯形平顶宽度 ,tb = ta + D , tc 为波形总宽度 tc = ta + tb , u( t) 为阶跃函数。 设 Ts 为采样周期 ,令 na = ta/ Ts , nb = tb/ Ts , nc = tc/ Ts ,对ν o 进行 Z 变换得 V o ( z) = Az (1 - z- na - z- nb + z- nc ) / na ( z - 1) 2 　　单指数信号ν i = Ae- t/τ,τ为时间常数 ,进 行 Z变换得 V i ( z) = A z/ ( z - e- Ts/τ) , 令 a = e- Ts/τ ,则单指数信号梯形成形滤波器的传递函数为 H( z) = [ z- 1 - z- ( na+1) - z- ( nb+1) + z- ( nc+1) - a( z- 2 - z- ( na+2) - z- ( nb+2) + z- ( nc +2) ) ]/ na (1 - 2z- 1 + z- 2 ) 2 　软件的设计 根据以上的算法 ,在 MA TLAB 下编写了单指数信号梯形滤波器程序 TRASfilter. m ,为 58 了设计出能根据需要直接调节各参数的梯形滤波器软件 ,本文按如图 2 所示的流程图设计了梯形滤波器分析程序。 图 2 　程序流程图 该软件可以实现以下