CN103837884A - 基于时域分析的数字核脉冲信号梯形成形算法 - Google Patents

基于时域分析的数字核脉冲信号梯形成形算法 Download PDF

Info

Publication number
CN103837884A
CN103837884A CN201410065940.9A CN201410065940A CN103837884A CN 103837884 A CN103837884 A CN 103837884A CN 201410065940 A CN201410065940 A CN 201410065940A CN 103837884 A CN103837884 A CN 103837884A
Authority
CN
China
Prior art keywords
trapezoidal
impulse response
digital
formation system
pulse signal
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
CN201410065940.9A
Other languages
English (en)
Other versions
CN103837884B (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 CN201410065940.9A priority Critical patent/CN103837884B/zh
Publication of CN103837884A publication Critical patent/CN103837884A/zh
Application granted granted Critical
Publication of CN103837884B publication Critical patent/CN103837884B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Radiation (AREA)

Abstract

本发明公开了一种基于时域分析的数字核脉冲信号梯形成形算法，包括：根据模拟梯形成形系统的电路，在时域中推导出模拟梯形成形系统的冲激响应；根据采样定理将模拟梯形成形系统由模拟域转到数字域，得到数字梯形成形系统的冲激响应；在时域中，将数字核脉冲信号与数字梯形成形系统的冲激响应进行卷积和运算，实现数字核脉冲的梯形成形。

Description

基于时域分析的数字核脉冲信号梯形成形算法
技术领域
本发明涉及放射性测量中数字核脉冲信号的梯形成形，尤其涉及一种基于时域分析的数字核脉冲信号梯形成形算法。
背景技术
在核能谱测量系统中，为了减小噪声、弹道亏损等对能谱测量的影响, ,需对核脉冲信号进行滤波成形。梯形成形可减小或消除弹道亏损，当梯形平顶宽度大于探测器最大电荷收集时间时，可避免弹道亏损。因此，梯形成形是一种用于核脉冲信号滤波成形的重要方法。与模拟梯形成形相比,数字梯形成形在调整梯形的宽度时，不需对硬件进行调整，具有更高的灵活性和稳定性。文中基于模拟梯形成形电路，在时域中推导出了数字梯形成形系统的冲激响应，并在时域中用卷积和运算实现对核脉冲信号的滤波成形，为数字核脉冲信号梯形成形的实现提供了一种实现方式。
发明内容
为解决上述技术问题，本发明的目的是提供一种基于冲激响应不变法的数字核脉冲梯形成形方法。该方法克服了核脉冲模拟梯形成形的不足，实现数字核脉冲信号的梯形成形。
本发明的目的通过以下的技术方案来实现：
基于时域分析的数字核脉冲信号梯形成形算法，包括：
根据模拟梯形成形系统的电路，在时域中推导出模拟梯形系统的冲激响应；
对模拟梯形成形系统的冲激响应进行傅里叶变换，得到模拟梯形成形系统的频率响应；
根据采样定理及模拟梯形成形系统的频宽确定对模拟梯形成形系统的冲激响应进行采样的频率后，对模拟梯形成形系统的冲激响应进行采样，得到数字梯形成形系统的冲激响应；
对数字梯形成形系统的冲激响应进行截取，冲激响应的长度越长，成形梯形越宽；
将数字核脉冲信号与选定长度的数字梯形成形系统的冲激响应进行卷积和运算，实现数字核脉冲梯形成形。
 与现有技术相比，本发明的一个或多个实施例可以具有如下优点：
有效克服模拟梯形成形系统的不足，将数字核脉冲信号成形为准梯形波形，在调节参数值时，不需对硬件进行调整。使用中可根据对成形波形的需要灵活调节参数值, 以满足不同的测量需求。 
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例共同用于解释本发明，并不构成对本发明的限制。在附图中：
图1是基于时域分析的数字核脉冲信号梯形成形算法的流程图；
图2是模拟梯形成形系统的电路原理图；
图3（a）是模拟梯形成形系统的频率响应的幅度谱；
（b）是数字梯形成形系统的冲激响应h(n)的波形；
图4是对实测核脉冲信号的梯形成形结果。
具体实施方式
容易理解，根据本发明的技术方案，在不变更本发明的实质精神下，本领域的一般技术人员可以提出本发明的多个结构方式和制作方法。因此以下具体实施方式以及附图仅是本发明的技术方案的具体说明，而不应当视为本发明的全部或者视为本发明技术方案的限定或限制。
下面结合实施例及附图对本发明作进一步详细的描述。
如图1所示，为基于时域分析的数字核脉冲信号梯形成形算法，该方法包括以下步骤：
步骤10 根据模拟梯形成形系统的电路，在时域中推导出模拟梯形系统的冲激响应h(t)，包括如下步骤A-D ：
A 上述模拟梯形成形系统的电路原理图如图2所示，整个系统由三个子系统级联组成，系统1的输入信号f(t)与输出信号y1(t)之间的关系为：
Figure 696670DEST_PATH_IMAGE001
    （1）
在时域中对（1）式的微分方程进行求解，得系统1的冲激响应h1(t)为： 
Figure 810119DEST_PATH_IMAGE002
(2)
B 同理，由系统2的电路得：
          (3)
系统2的冲激响应为： 
Figure 897341DEST_PATH_IMAGE004
      (4)
C 由系统3的电路得： 
Figure 958838DEST_PATH_IMAGE005
       (5)
Figure 180872DEST_PATH_IMAGE006
(6)
系统3的冲激响应为： 
Figure 593398DEST_PATH_IMAGE007
 (7)
D 模拟高斯成形系统的冲激响应为三个子系统冲激响应的卷积积分，即：
Figure 621397DEST_PATH_IMAGE008
 (8)
将（2）式、（4）式、7（式）代入（8）式得：
Figure 802980DEST_PATH_IMAGE009
（9）
其中，
Figure 830159DEST_PATH_IMAGE011
 
Figure 396269DEST_PATH_IMAGE012
Figure 697938DEST_PATH_IMAGE013
 
步骤20 对模拟梯形成形系统的冲激响应进行傅里叶变换，得到模拟梯形成形系统的频率响应。
对（9）式进行傅里叶变换得到模拟高斯成形系统的频率响应H(ω)为：
Figure 497660DEST_PATH_IMAGE014
（10）
步骤30 根据采样定理及模拟梯形成形系统的频宽，确定对模拟梯形成形系统的冲激响应进行采样的频率后，对模拟梯形成形系统的冲激响应进行采样，得到数字梯形成形系统的冲激响应。
根据(10)式在matlab中画出模拟成形系统频率响应的幅度谱，以确定H(ω)的频宽，其频宽与参数值有关，以k1=10，k2=1，k3=200，k4=1，k5=2为例，模拟成形系统频率响应的幅度谱见图3（a）所示，由图可知，频宽约10000rad/s，因此，采样的角频率ω应大于20000rad/s，则采样频率f应大于
Figure 884779DEST_PATH_IMAGE015
。取采样频率为20MHz，即采样周期为Ts为0.05μs, 对模拟梯形成形系统的冲激响应进行采样，得到数字梯形成形系统的冲激响应，即使数字梯形成形系统的冲激响应h(n)等于模拟梯形成形系统的冲激响应h(t)的采样值，即：
  
Figure 254580DEST_PATH_IMAGE016
 (11)
根据(9)式和（11）式得：
Figure 145176DEST_PATH_IMAGE017
  （12）
在Matlab中画出h(n)的波形见图3（b）所示。
步骤40 对数字梯形成形系统的冲激响应进行截取，冲激响应的长度越长，成形梯形越宽。由图3（b）可知，h(n)的长度较长，在250000点左右衰减为0，实际使用时，对其进行截取。令 
Figure 879914DEST_PATH_IMAGE018
（13）
步骤50 将数字核脉冲信号与选定长度的数字梯形成形系统的冲激响应进行卷积和运算，实现数字核脉冲梯形成形。对于数字梯形成形系统，输入信号f(n)经过该系统的响应y(n)为： 
 
Figure 488750DEST_PATH_IMAGE019
(14)
图4为实测60Co的核脉冲及其滤波成形后的波形，其中，N为h(n)的长度。各参数取值分别为：k1=10、k2=1、k3=200、k4=1,。由图中可以看出，成形后的信号在滤除噪声的同时，将核脉冲信号成形为梯形波形。N值越大，成形后的波形越宽；在改变N值时，需同时改变k5的值，否则梯形顶部不平坦。应用中可根据实际需求选择合适的参数值，以满足不同的测量需求。
 虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (4)

1.基于时域分析的数字核脉冲信号梯形成形算法，其特征在于，所述方法包括：
根据模拟梯形成形系统的电路，在时域中推导出模拟梯形系统的冲激响应；
对模拟梯形成形系统的冲激响应进行傅里叶变换，得到模拟梯形成形系统的频率响应，根据采样定理及模拟梯形成形系统的频宽确定对模拟梯形成形系统的冲激响应进行采样的频率后，对模拟梯形成形系统的冲激响应进行采样，得到数字梯形成形系统的冲激响应；
对数字梯形成形系统的冲激响应进行截取，冲激响应的长度越长，成形梯形越宽；
将数字核脉冲信号与选定长度的数字梯形成形系统的冲激响应进行卷积和运算，实现数字核脉冲信号的梯形成形。
2.根据权利要求1所述的基于时域分析的数字核脉冲信号梯形成形算法，其特征在于，所述模拟梯形系统冲激响应的推导在时域中进行。
3.根据权利要求1所述的基于时域分析的数字核脉冲信号梯形成形算法，其特征在于，数字梯形成形系统冲激响应的长度越长，成形梯形越宽。
4.根据权利要求1所述的基于时域分析的数字核脉冲信号梯形成形算法，其特征在于，改变数字梯形成形系统冲激响应的长度，需同时调整相关参数的取值，使成形梯形的顶部平坦。
CN201410065940.9A 2014-02-26 2014-02-26 基于时域分析的数字核脉冲信号梯形成形算法 Expired - Fee Related CN103837884B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410065940.9A CN103837884B (zh) 2014-02-26 2014-02-26 基于时域分析的数字核脉冲信号梯形成形算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410065940.9A CN103837884B (zh) 2014-02-26 2014-02-26 基于时域分析的数字核脉冲信号梯形成形算法

Publications (2)

Publication Number Publication Date
CN103837884A true CN103837884A (zh) 2014-06-04
CN103837884B CN103837884B (zh) 2016-02-10

Family

ID=50801583

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410065940.9A Expired - Fee Related CN103837884B (zh) 2014-02-26 2014-02-26 基于时域分析的数字核脉冲信号梯形成形算法

Country Status (1)

Country Link
CN (1) CN103837884B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104123473A (zh) * 2014-07-31 2014-10-29 成都理工大学 时域内核脉冲信号放大电路的数值分析方法
CN106772545A (zh) * 2017-02-09 2017-05-31 成都理工大学 一种采用尖顶成形算法的数字脉冲幅度分析器
CN107783173A (zh) * 2017-10-24 2018-03-09 东华理工大学 一种数字化核能谱测量系统中的脉冲矩形成形方法
CN107817514A (zh) * 2017-10-24 2018-03-20 东华理工大学 一种数字化核能谱测量系统中的脉冲阶梯成形方法
CN110276441A (zh) * 2019-05-20 2019-09-24 成都理工大学 一种基于深度学习的梯形重叠核脉冲估计方法
CN110347965A (zh) * 2019-07-17 2019-10-18 四川新先达测控技术有限公司 核脉冲信号处理方法及装置
CN112327347A (zh) * 2020-10-29 2021-02-05 中广核久源（成都）科技有限公司 曲率可调的数字化核脉冲成形系统
CN114252899A (zh) * 2022-03-02 2022-03-29 四川新先达测控技术有限公司 一种核信号的级联冲激卷积成形方法和装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08189974A (ja) * 1995-01-11 1996-07-23 Toshiba Corp 放射線エネルギースペクトル測定装置
CN101710183A (zh) * 2009-12-31 2010-05-19 中国原子能科学研究院 用于核谱学及核电子学的数字符合多道系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08189974A (ja) * 1995-01-11 1996-07-23 Toshiba Corp 放射線エネルギースペクトル測定装置
US5781142A (en) * 1995-01-11 1998-07-14 Kabushiki Kaisha Toshiba Method and apparatus for converting an analog measurement signal to a digital signal having reduced conversion error
CN101710183A (zh) * 2009-12-31 2010-05-19 中国原子能科学研究院 用于核谱学及核电子学的数字符合多道系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ZHOU JIANBIN ET AL.: "Study of time-domain digital pulse shaping algorithms for nuclear signals", 《NUCLEAR SCIENCE AND TECHNIQUES》 *
吴军龙等: "基于 NI MUlTISIM 的核信号滤波成形电路研究", 《核电子学与探测技术》 *
肖无云等: "数字化多道脉冲幅度分析技术研究", 《核技术》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104123473A (zh) * 2014-07-31 2014-10-29 成都理工大学 时域内核脉冲信号放大电路的数值分析方法
CN106772545A (zh) * 2017-02-09 2017-05-31 成都理工大学 一种采用尖顶成形算法的数字脉冲幅度分析器
CN107817514B (zh) * 2017-10-24 2020-01-07 东华理工大学 一种数字化核能谱测量系统中的脉冲阶梯成形方法
CN107783173A (zh) * 2017-10-24 2018-03-09 东华理工大学 一种数字化核能谱测量系统中的脉冲矩形成形方法
CN107817514A (zh) * 2017-10-24 2018-03-20 东华理工大学 一种数字化核能谱测量系统中的脉冲阶梯成形方法
CN107783173B (zh) * 2017-10-24 2019-06-11 东华理工大学 一种数字化核能谱测量系统中的脉冲矩形成形方法
CN110276441B (zh) * 2019-05-20 2021-05-04 成都理工大学 一种基于深度学习的梯形重叠核脉冲估计方法
CN110276441A (zh) * 2019-05-20 2019-09-24 成都理工大学 一种基于深度学习的梯形重叠核脉冲估计方法
CN110347965A (zh) * 2019-07-17 2019-10-18 四川新先达测控技术有限公司 核脉冲信号处理方法及装置
CN110347965B (zh) * 2019-07-17 2023-02-03 四川新先达测控技术有限公司 核脉冲信号处理方法及装置
CN112327347A (zh) * 2020-10-29 2021-02-05 中广核久源（成都）科技有限公司 曲率可调的数字化核脉冲成形系统
CN114252899A (zh) * 2022-03-02 2022-03-29 四川新先达测控技术有限公司 一种核信号的级联冲激卷积成形方法和装置
CN114252899B (zh) * 2022-03-02 2022-05-20 四川新先达测控技术有限公司 一种核信号的级联冲激卷积成形方法和装置

Also Published As

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

Similar Documents

Publication Publication Date Title
CN103837884A (zh) 基于时域分析的数字核脉冲信号梯形成形算法
MX2017012567A (es) Método para procesar datos sísmicos proporcionando recopilaciones de imágenes comunes de apertura de superficie.
GB2556237A (en) Adaptive phase-distortionless magnitude response equalization for beamforming applications
CN103777228B (zh) 基于iir滤波器的数字核脉冲信号高斯成形方法
CN102176667B (zh) 一种基于频率跟踪的工频数字陷波装置与方法
CN102998528B (zh) 动态正弦畸变信号中频率分量的提取方法及其实现装置
CN104242877B (zh) 基于模拟cr‑rc电路的核脉冲信号数字高斯成形方法
CN106679659B (zh) 一种基于参数可调非线性跟踪微分器的信号去噪方法
CN104062498B (zh) 对称频谱信号的中心频率的质心估计方法
CN103871422A (zh) 基于三因子迭代联合块对角化的时域混叠盲信号分离方法
CN103941280B (zh) 基于冲激响应不变法的数字核脉冲高斯成形方法
Tsiklauri et al. Causality and delay and physics in real systems
Domoshnitsky et al. Nonoscillation and positivity of solutions to first order state-dependent differential equations with impulses in variable moments
CN103064119B (zh) 一种时域瞬变电磁探测仪器的固定高采样率设计方法
CN109506762A (zh) 基于滤波器的水听器接收信号修正方法
CN202364188U (zh) 一种用于发电机转子电压测量的滤波电路
US9077316B2 (en) Transmitter finite impulse response characterization
CN103777221B (zh) 基于窗函数法的数字核脉冲信号高斯成形方法
CN106814394B (zh) 台站式检波器与常规检波器联合低频恢复方法
CN104545913A (zh) 一种应用于环境磁场抑制的自适应补偿方法及系统
CN104237637A (zh) 阻波器阻塞电阻值的测量方法
Zhao Modeling and simulation of digital filter
Guo Nonlinear impulsive fractional differential equations in Banach spaces
CN108241171B (zh) 一种采用复值Gauss积分滤波器对地震数据进行滤波和提取三瞬属性的方法
CN103630935B (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