CN103777221A - 基于窗函数法的数字核脉冲信号高斯成形方法 - Google Patents

基于窗函数法的数字核脉冲信号高斯成形方法 Download PDF

Info

Publication number
CN103777221A
CN103777221A CN201410066013.9A CN201410066013A CN103777221A CN 103777221 A CN103777221 A CN 103777221A CN 201410066013 A CN201410066013 A CN 201410066013A CN 103777221 A CN103777221 A CN 103777221A
Authority
CN
China
Prior art keywords
window
window function
pulse signal
filter
frequency
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410066013.9A
Other languages
English (en)
Other versions
CN103777221B (zh
Inventor
葛青
葛良全
罗耀耀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN201410066013.9A priority Critical patent/CN103777221B/zh
Publication of CN103777221A publication Critical patent/CN103777221A/zh
Application granted granted Critical
Publication of CN103777221B publication Critical patent/CN103777221B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Manipulation Of Pulses (AREA)

Abstract

本发明公开了一种基于窗函数法的数字核脉冲信号高斯成形方法，包括以下步骤：根据模拟高斯成形滤波器的频率响应的幅度谱，确定FIR数字滤波器的性能指标，即FIR数字滤波器通带截止频率、阻带截止频率以及阻带最小衰减；根据阻带最小衰减来选择窗函数，常用的窗函数有：矩形窗、三角窗、汗宁窗、海明窗、布莱克曼窗、凯塞窗等；根据选择的窗函数、过渡带宽以及理想FIR数字滤波器的截止频率，在MATLAB中计算出对应的FIR数字滤波器的系数；数字核脉冲信号通过得到的FIR数字滤波器后，被滤波成形为准高斯波形。该方法克服了核脉冲信号模拟高斯成形的不足，成形后的波形具有较好的准高斯特性，并且可以实现线性相位。

Description

基于窗函数法的数字核脉冲信号高斯成形方法
技术领域
本发明涉及放射性测量中数字核脉冲信号的高斯成形，尤其涉及一种基于窗函数法的数字核脉冲信号高斯成形方法。
背景技术
在核仪器中，通常要采用滤波成形电路,将核探测器输出的核脉冲信号进行滤波并成形为需要的波形,以满足后续分析的需求。因为高斯型波形具有弹道亏损较小、脉冲顶部较平坦等优点,所以核脉冲的滤波成形常以成形为高斯波形或准高斯波形为目标。Sallen-Key低通滤波器是一种核仪器中常用的模拟滤波成形电路，可以将核脉冲信号滤波成形为准高斯波形。由于数字滤波器具有更高的灵活性、稳定性等优点，用数字滤波器实现对核脉冲信号的滤波成形愈来愈受到关注。有限冲激响应（FIR）数字滤波器的设计方法有窗函数法、频率采样法等，窗函数法由于有公式可循，使用方便，成为工程实际中广泛应用的方法。基于窗函数法对数字核脉冲信号进行滤波成形，在满足幅频响应要求的同时，可以获得线性相位。
发明内容
本发明的目的在于公开一种基于窗函数法的数字核脉冲信号高斯成形方法，该方法克服了核脉冲信号模拟高斯成形的不足，在滤除噪声的同时，将数字核脉冲信号成形为准高斯波形，为实现核脉冲信号的数字高斯成形提供了一种实现方法。
本发明是通过以下技术方案实现的，具体包括以下步骤：
根据模拟高斯成形系统的电路，得到系统电路输入信号与输出信号的微分方程，将微分方程在频域中进行求解，得到模拟高斯成形系统的频率响应；
根据模拟高斯成形系统的频率响应的幅度谱，确定FIR数字滤波器的指标，即FIR数字滤波器的通带截止频率、阻带截止频率以及阻带最小衰减；
根据阻带最小衰减来选择窗函数，常用的窗函数有：矩形窗、三角窗、汗宁窗、海明窗、布莱克曼窗、凯塞窗等；
根据选择的窗函数、过渡带宽以及理想数字滤波器的截止频率，在MATLAB中计算出对应的FIR数字滤波器的系数； 
数字核脉冲信号通过得到的FIR数字滤波器后，被滤波成形为准高斯波形。
 与现有技术相比，本发明的一个或多个实施例可以具有如下优点：
有效克服模拟高斯成形系统的不足，在滤除噪声的同时，实现对数字核脉冲信号的成形，成形后的波形具有较好的准高斯特性，并且可以获得线性相位。
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例共同用于解释本发明，并不构成对本发明的限制。在附图中：
图1是基于FIR滤波器的数字核脉冲信号高斯成形方法流程图；
图2是模拟高斯成形系统的电路原理图；
图3是模拟高斯成形系统频率响应的幅度谱；
图4是实测核脉冲信号通过矩形窗设计的FIR数字滤波器后的响应y(n)的波形。
具体实施方式
容易理解，根据本发明的技术方案，在不变更本发明的实质精神下，本领域的一般技术人员可以提出本发明的多个结构方式和制作方法。因此以下具体实施方式以及附图仅是本发明的技术方案的具体说明，而不应当视为本发明的全部或者视为本发明技术方案的限定或限制。
下面结合实施例及附图对本发明作进一步详细的描述。
如图1所示，为基于窗函数法的数字核脉冲信号高斯成形方法，该方法包括以下步骤：
步骤10 根据模拟高斯成形系统的电路，得到系统电路输入信号与输出信号的微分方程，将微分方程在频域中进行分析，推导出模拟高斯成形系统的频率响应，即模拟Sallen-Key低通滤波器的频率响应，包括如下步骤A1-B1：
A1 根据模拟Sallen-Key低通滤波器的电路原理图（如图2所示），列出图中输入信号f(t)与输出信号y(t)之间的数学关系为:
Figure 181676DEST_PATH_IMAGE001
                 (1)
B1 令f(t)的傅里叶变换为F(Ω)，y(t)的傅里叶变换为Y(Ω)，对方程两边同时进行傅里叶变换，得到      
Figure 679654DEST_PATH_IMAGE002
              (2) 
则，模拟Sallen-Key低通滤波器的频率响应函数H(Ω)为：
                                 
Figure 451301DEST_PATH_IMAGE003
    (3)
    步骤20 根据模拟Sallen-Key低通滤波器的频率响应的幅度谱，确定FIR数字滤波器的技术指标，即确定FIR数字滤波器的通带截止频率、阻带截止频率以及阻带最小衰减；FIR数字滤波器指标的确定包括如下步骤A2-C2：
A2 根据（3）式，在MATLAB中画出模拟Sallen-Key低通滤波器的频率响应的幅度谱，该幅度谱的频宽随RC的增大而减小，当RC分别为0.000005、0.000009时，模拟高斯成形系统的频率响应的幅度谱见图3所示。                                     
B2 根据模拟Sallen-Key低通滤波器的频率响应的幅度谱，确定FIR模拟滤波器的指标，即通带截止频率Ωp、阻带截止频率Ωs以及阻带最小衰减αss，以图3中RC=0.000005时的模拟高斯成形系统频率响应的幅度谱为例，指标如下：
Ωp=250000 rad/s
Ωs=630000 rad/s
Figure 464256DEST_PATH_IMAGE004
C2 令采样频率为fs，当采样频率为fs=20MHz时，根据模拟Sallen-Key低通滤波器的指标，得到FIR数字滤波器的通带截止频率ωp、阻带截止频率ωs以及阻带最小衰减αs分别为：
ωp=Ωp /fs=0.0125
ωs=Ωs /fs=0.0315
Figure 346761DEST_PATH_IMAGE005
步骤30 根据阻带最小衰减来确定选用的窗函数，常用的窗函数有：矩形窗、三角窗、汗宁窗、海明窗、布莱克曼窗、凯塞窗等，在本发明的具体实施例中，由于阻带最小衰减为20dB，则上述窗函数均能满足阻带最小衰减的要求，因此均可选择；
步骤40 根据选择的窗函数、过渡带宽以及理想数字滤波器的截止频率，在MATLAB中计算出对应的FIR数字滤波器的系数，包括如下步骤A3-C3：
A3过渡带宽Δω及理想数字滤波器的截止频率ωc的计算：
Δω=ωs-ωp              (4)
ωc=(ωs-ωp)/2       (5)
    B3 根据Δω及选用的窗函数，计算滤波器的阶数N，根据选用的窗函数，查窗函数的基本参数表得到A值，则： 
       
Figure 699245DEST_PATH_IMAGE006
              (6)
C3 根据选用的窗函数、N以及ωc，在MATLAB中计算出滤波器的系数b，其中b是一个长度为N+1的向量，则FIR数字滤波器的系统函数为：
Figure 641793DEST_PATH_IMAGE007
                          (7)
步骤50 数字核脉冲信号通过得到的FIR数字滤波器后，被滤波成形为准高斯波形。数字核脉冲信号f(n)通过FIR数字滤波器的响应y(n)为： 
  
Figure 142045DEST_PATH_IMAGE008
      (8)
     图4为实测60Co的数据通过选用矩形窗设计的FIR数字滤波器后的响应y(n)的波形，从图中可以看出数字核脉冲信号被滤波成形为准高斯信号，成形后的波形具有较好的准高斯特性。
虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (2)

1.基于窗函数法的数字核脉冲信号高斯成形方法，其特征在于，所述方法包括：
（1）根据模拟高斯成形滤波器的电路，得到系统电路输入信号与输出信号的微分方程，在频域中求解出模拟高斯成形滤波器的频率响应；
（2）根据模拟高斯成形滤波器的频率响应的幅度谱，确定FIR数字滤波器的技术指标，即确定FIR数字滤波器的通带截止频率、阻带截止频率以及阻带最小衰减；
（3）根据阻带最小衰减来选择窗函数，常用的窗函数有：矩形窗、三角窗、汗宁窗、海明窗、布莱克曼窗、凯塞窗等；
（4）根据选择的窗函数、过渡带宽以及理想FIR数字滤波器的截止频率，在MATLAB中计算出对应的FIR数字滤波器的系数； 
（5）数字核脉冲信号通过得到的FIR数字滤波器后，被滤波成形为准高斯波形。
2.根据权利要求1所述的基于窗函数法的数字核脉冲信号高斯成形方法，其特征在于，所述（3）中根据阻带最小衰减来选择窗函数，可以选择的窗函数有：矩形窗、三角窗、汗宁窗、海明窗、布莱克曼窗、凯塞窗等，选择不同的窗函数计算出的滤波器的系数不同，因此，推导出的FIR数字滤波器的系统函数不同，但都可实现对数字核脉冲信号的高斯成形。
CN201410066013.9A 2014-02-26 2014-02-26 基于窗函数法的数字核脉冲信号高斯成形方法 Expired - Fee Related CN103777221B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410066013.9A CN103777221B (zh) 2014-02-26 2014-02-26 基于窗函数法的数字核脉冲信号高斯成形方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410066013.9A CN103777221B (zh) 2014-02-26 2014-02-26 基于窗函数法的数字核脉冲信号高斯成形方法

Publications (2)

Publication Number Publication Date
CN103777221A true CN103777221A (zh) 2014-05-07
CN103777221B CN103777221B (zh) 2016-02-10

Family

ID=50569683

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410066013.9A Expired - Fee Related CN103777221B (zh) 2014-02-26 2014-02-26 基于窗函数法的数字核脉冲信号高斯成形方法

Country Status (1)

Country Link
CN (1) CN103777221B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103941279A (zh) * 2014-05-08 2014-07-23 东华理工大学 基于脉宽和幅度双参数的核脉冲信号数字高斯成形方法
CN112327347A (zh) * 2020-10-29 2021-02-05 中广核久源（成都）科技有限公司 曲率可调的数字化核脉冲成形系统
CN113189634A (zh) * 2021-03-02 2021-07-30 四川新先达测控技术有限公司 一种类高斯成形方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007025993A2 (fr) * 2005-08-30 2007-03-08 Thales Procede et dispositif de controle de puissance crete et de largeur d'impulsion d'un emetteur de forte puissance rf impulsionnel gaussien large bande
CN101236243A (zh) * 2007-02-02 2008-08-06 株式会社东芝 脉冲信号发送装置、其波形调整方法及dme地面站装置
CN101309092A (zh) * 2008-07-14 2008-11-19 哈尔滨工业大学 分数傅立叶变换域超宽带脉冲信号调制与解调方法
CN102983838A (zh) * 2012-12-05 2013-03-20 天津光电通信技术有限公司 一种基于fpga实现高斯滤波器数字逻辑电路的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007025993A2 (fr) * 2005-08-30 2007-03-08 Thales Procede et dispositif de controle de puissance crete et de largeur d'impulsion d'un emetteur de forte puissance rf impulsionnel gaussien large bande
CN101236243A (zh) * 2007-02-02 2008-08-06 株式会社东芝 脉冲信号发送装置、其波形调整方法及dme地面站装置
CN101309092A (zh) * 2008-07-14 2008-11-19 哈尔滨工业大学 分数傅立叶变换域超宽带脉冲信号调制与解调方法
CN102983838A (zh) * 2012-12-05 2013-03-20 天津光电通信技术有限公司 一种基于fpga实现高斯滤波器数字逻辑电路的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
何小河: "基于窗函数的FIR滤波器的设计", 《四川理工学院学报（自然科学版）》 *
周伟 等: "基于Sallen-Key滤波器的数字高斯成形方法的仿真", 《系统仿真学报》 *
祁中 等: "基于低通S-K滤波器的核脉冲成形电路", 《兰州大学学报（自然科学版）》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103941279A (zh) * 2014-05-08 2014-07-23 东华理工大学 基于脉宽和幅度双参数的核脉冲信号数字高斯成形方法
CN112327347A (zh) * 2020-10-29 2021-02-05 中广核久源（成都）科技有限公司 曲率可调的数字化核脉冲成形系统
CN113189634A (zh) * 2021-03-02 2021-07-30 四川新先达测控技术有限公司 一种类高斯成形方法
CN113189634B (zh) * 2021-03-02 2022-10-25 四川新先达测控技术有限公司 一种类高斯成形方法

Also Published As

Publication number Publication date
CN103777221B (zh) 2016-02-10

Similar Documents

Publication Publication Date Title
CN103777228B (zh) 基于iir滤波器的数字核脉冲信号高斯成形方法
CN103185837B (zh) 一种电力系统频率测量的方法
CN103837884B (zh) 基于时域分析的数字核脉冲信号梯形成形算法
CN108764073B (zh) 一种结合频谱能量形态拟合的加速度滤噪和积分方法
CN103777221B (zh) 基于窗函数法的数字核脉冲信号高斯成形方法
CN109142863B (zh) 一种电力系统测频方法及系统
Singh et al. Analysis of digital IIR filter with Labview
CN105785124A (zh) 一种采用谱估计和互相关的电力系统谐波和间谐波测量方法
Patel et al. Design technique of bandpass FIR filter using various window function
CN103823177A (zh) 基于窗函数设计的滤波器性能检测方法和系统
CN104991119B (zh) 一种消除伪峰、谱泄漏效应的互素谱分析方法及其装置
CN104242877B (zh) 基于模拟cr‑rc电路的核脉冲信号数字高斯成形方法
CN104502702A (zh) 检测电力信号的频率的方法和系统
CN111948454A (zh) 一种抗衰减直流分量干扰的同步相量测量方法及系统
CN109460614B (zh) 基于瞬时带宽的信号时间-频率分解方法
CN111505375A (zh) 一种频率测量方法及装置
CN104504272B (zh) 基于理想重构滤波器冲激响应拟合的数字积分器实现方法
CN110285881A (zh) 一种基于全相位滤波的密集谱频率估计法
CN102338890B (zh) 一种地球物理勘探中的圆形窗带通保幅滤波数据处理方法
CN103941280A (zh) 基于冲激响应不变法的数字核脉冲高斯成形方法
CN107609305A (zh) 基于vc的滤波器辅助设计系统
CN109327203B (zh) 一种基于二次对称滤波的数字下变频方法
Song Low-pass filter design and sampling theorem verification
CN106767940B (zh) 一种抑制供电频率干扰的方法
CN109728860B (zh) 一种基于采集终端检测装置的通讯干扰抑制方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160210

Termination date: 20170226

CF01 Termination of patent right due to non-payment of annual fee