CN114252899B - 一种核信号的级联冲激卷积成形方法和装置 - Google Patents

一种核信号的级联冲激卷积成形方法和装置 Download PDF

Info

Publication number
CN114252899B
CN114252899B CN202210196100.0A CN202210196100A CN114252899B CN 114252899 B CN114252899 B CN 114252899B CN 202210196100 A CN202210196100 A CN 202210196100A CN 114252899 B CN114252899 B CN 114252899B
Authority
CN
China
Prior art keywords
signal
impulse
formula
cascade
convolution
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.)
Active
Application number
CN202210196100.0A
Other languages
English (en)
Other versions
CN114252899A (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.)
Sichuan Xstar Measurement Control Technology Co ltd
Original Assignee
Sichuan Xstar Measurement Control Technology Co ltd
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 Sichuan Xstar Measurement Control Technology Co ltd filed Critical Sichuan Xstar Measurement Control Technology Co ltd
Priority to CN202210196100.0A priority Critical patent/CN114252899B/zh
Publication of CN114252899A publication Critical patent/CN114252899A/zh
Application granted granted Critical
Publication of CN114252899B publication Critical patent/CN114252899B/zh
Priority to US18/087,853 priority patent/US20230281266A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/17—Circuit arrangements not adapted to a particular type of detector
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Complex Calculations (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本发明公开了一种核信号的级联冲激卷积成形方法和装置，针对现有高斯滤波算法比较复杂，构造核信号的实时数字高斯滤波器有难度的问题，本发明的技术方案包括：通过探测器获取到探测器信号；探测器信号多级级联成形系统先与高斯信号进行卷积，然后再进行双指数冲激成形，经过这种多级级联成形系统，生成用于分析的脉宽很窄的高斯形态脉冲信号。根据复杂系统的多级级联卷积可以交换卷积顺序的特性，探测器信号可以先经过级联逆系统形成冲激信号，然后将冲激信号与高斯信号进行卷积，也能生成级联冲激卷积信号。本方法可以扩展到三指数或四指数信号进行高斯、梯形、CUSP、余弦平方分布成形和柯西分布成形等。

Description

一种核信号的级联冲激卷积成形方法和装置
技术领域
本发明属于核信号处理技术领域，具体涉及一种核信号的级联冲激卷积成形方法和装置。
背景技术
核信号携带了各种信息，比如辐射粒子的能量，辐射粒子的种类、辐射事件的发生时刻等。从核信号中提取的核信息可以用于核性质、核结构和核衰变等基础科学研究。在核科学技术中要得到精准的核信息，往往需要用电子学方法探测核信号，并从核信号中提取核信息。随着高速数字处理芯片及高速ADC的发展，核信号数字化及数字处理技术逐渐成熟。
现有技术中，对核信号数字处理方法主要以研究数字梯形、CUSP、高斯滤波器为主流，锯齿滤波器作为PSD（脉冲形状甄别）研究，冲激滤波器等用于高计数率研究的补充方法。
现有技术存在以下技术问题：
1.现有高斯滤波器的抗噪能力好，但是算法比较复杂，较耗费硬件资源，构造实时的数字高斯滤波器有难度；
2. 数字芯片上部署数字算法技术门槛高，开发周期长，浮点运算能力差，无法部署复杂算法。同时，数字芯片资源有限，算法的实现受到严重的限制。
发明内容
针对上述现有技术中存在的问题，本发明提出了一种核信号的级联冲激卷积成形方法和装置，其目标为：通过将信号采用精细的双指数冲激成形，再级联卷积标准数字高斯信号实现了信号的高斯成形，并将数字高斯信号直接卷积双指数冲激成形滤波器信号实现了核信号数字高斯成形滤波器，并且可以扩展到三指数或四指数高斯成形、COS平方分布成形和柯西分布成形等。
为实现上述目标本发明所采用的技术方案是：提供一种核信号的级联冲激卷积成形方法，包括：
S1：通过探测器获取到探测器信号；
S2：将探测器信号作为输入信号，将输入信号经过多级级联成形系统，输入信号先通过级联卷积系统与目标信号进行卷积，然后再通过级联逆系统进行冲激成形，生成用于分析的级联冲激卷积信号，并得到级联冲激卷积信号的函数表达式；
S3:然后基于多级级联成形系统的级联卷积支持交换卷积顺序的特性，将输入信号经过级联逆系统进行冲激成形，形成冲激信号，并得到输入信号冲激成形的系统函数表达式；
S4：将冲激信号通过级联卷积系统与目标信号进行卷积，生成级联冲激卷积成形信号，并得到多级级联成形系统的函数表达式。
较优的，本发明S2中：
所述目标信号包括：标准高斯信号、余弦平方信号、柯西分布信号和梯形信号，将冲激信号与目标信号进行卷积，然后再通过级联逆系统进行冲激成形，生成探测器信号的级联冲激卷积信号，其得到函数表达式如公式（23）所示：
Figure 992503DEST_PATH_IMAGE001
*h[n]（23）；
其中，
Figure 311489DEST_PATH_IMAGE002
表示级联冲激卷积信号的函数表达式，
Figure 858926DEST_PATH_IMAGE003
表示输入信号的函数表达式，
Figure 605165DEST_PATH_IMAGE004
表示标准高斯信号表达式，h[n]表示双指数冲激脉冲成形系统函数表达式，n为采集的输入信号点序列。
较优的，本发明S3具体为：
S3.1：定义输入信号
Figure 100868DEST_PATH_IMAGE005
为双指数信号，将输入信号
Figure 223545DEST_PATH_IMAGE005
输入第一级INV_RC系统，输出单指数衰减信号
Figure 121094DEST_PATH_IMAGE006
，所述输入信号为
Figure 772655DEST_PATH_IMAGE005
由式（1）所示，所述单指数衰减信号
Figure 83551DEST_PATH_IMAGE006
由式（2）所示：
Figure 682022DEST_PATH_IMAGE007
(1)
Figure 496395DEST_PATH_IMAGE008
(2)；
其中，m和M是双指数信号的系统参数，n为采集的输入信号点序列；INV是inverse反向，RC指的是电路中的电阻R和电容C，连起来指的是逆RC,就是通过逆向运算，把电路中的RC影响去掉。
S3.2：将所述单指数衰减信号
Figure 787699DEST_PATH_IMAGE006
输入第二级INV_RC系统，输出得到冲激响应信号
Figure 320311DEST_PATH_IMAGE009
，冲激响应信号
Figure 519211DEST_PATH_IMAGE009
由式（3）所示：
Figure 125773DEST_PATH_IMAGE010
(3)
S3.3：根据INV_RC算子的数字解推导式(2)和式(3)，得到公式(4)和式(5)，如下所示：
Figure 853558DEST_PATH_IMAGE011
(4)
Figure 607887DEST_PATH_IMAGE012
(5)
这里INV_RC中的INV是inverse反向，RC指的是电路中的电阻R和电容C，连起来指的是逆RC,就是通过逆向运算，把电路中的RC影响去掉。
将公式（4）代入公式（5），得到由输入信号
Figure 610478DEST_PATH_IMAGE005
通过级联逆系统变为冲激响应信号
Figure 337126DEST_PATH_IMAGE009
的数字变换表达式（6），如下所示：
Figure 32549DEST_PATH_IMAGE013
(6)
对公式（6）两边求差分整理得式（7）、式（8）和式（9），如下所示：
Figure 710393DEST_PATH_IMAGE014
(7)
Figure 516675DEST_PATH_IMAGE015
(8)
Figure 363408DEST_PATH_IMAGE017
(9)
S3.4：对公式（8）进行整理，可得公式（10）的形式，如下所示：
Figure 964154DEST_PATH_IMAGE019
Figure 100002_DEST_PATH_IMAGE020
(10)
通过公式（10）得到双指数信号冲激成形的系统函数表达式如式（11）所示：
h1[n]=
Figure 896338DEST_PATH_IMAGE021
(11)。
较优的，本发明当定义输入信号
Figure 100002_DEST_PATH_IMAGE022
为单指数信号时，m=0，通过公式（11）得到单指数信号冲激成形的系统函数表达式如式（12）所示：
Figure 709573DEST_PATH_IMAGE023
(12)。
较优的，本发明S3具体为：
S3.1：当定义输入信号
Figure 207550DEST_PATH_IMAGE005
为带反冲的双指数信号时，将输入信号
Figure 916880DEST_PATH_IMAGE005
输入第一级INV_RC系统，第一级INV_RC系统的输出信号为
Figure 398677DEST_PATH_IMAGE006
，所述输入信号为
Figure 484445DEST_PATH_IMAGE005
由式（13）所示，所述输出信号
Figure 102508DEST_PATH_IMAGE006
由式（14）所示：
Figure 248319DEST_PATH_IMAGE024
(13)
Figure 951833DEST_PATH_IMAGE025
(14)
S3.2：可知式（14）为第一级INV_RC系统的函数表达式，通过式（11）的得到过程推导出带反冲的双指数信号冲激成形的系统函数表达式如式（15）所示：
Figure DEST_PATH_IMAGE027
Figure 575712DEST_PATH_IMAGE028
(15)。
较优的，本发明S4中：
所述目标信号包括：标准高斯信号、余弦平方信号、柯西分布信号和梯形信号，将冲激信号通过级联卷积系统与目标信号进行卷积，生成级联冲激卷积成形信号，其函数表达式如式（16）所示：
Figure DEST_PATH_IMAGE029
（16）；
其中，
Figure 15658DEST_PATH_IMAGE030
表示级联冲激卷积成形信号的函数表达式，
Figure DEST_PATH_IMAGE031
表示双指数冲激脉冲成形系统函数表达式，
Figure 332370DEST_PATH_IMAGE004
表示标准高斯信号。
较优的，本发明当目标信号为标准高斯信号时，将冲激信号与标准高斯信号进行卷积，生成输入信号的级联冲激卷积成形信号，分别将式（11）、式（12）和式（15）代入式（16），然后基于多级级联成形系统的级联卷积支持交换卷积顺序的特性，得到多级级联成形系统的函数表达式如式（17）、式（18）和式（19）所示：
Figure 523180DEST_PATH_IMAGE032
=
Figure 950750DEST_PATH_IMAGE033
（17）
Figure 543406DEST_PATH_IMAGE034
=
Figure 31019DEST_PATH_IMAGE035
（18）
Figure 100002_DEST_PATH_IMAGE036
=
Figure 709125DEST_PATH_IMAGE037
（19）；
其中，
Figure 100002_DEST_PATH_IMAGE038
表示目标信号的函数表达式，此处目标信号采用标准高斯信号。
较优的，本发明将标准高斯信号替换成余弦平方信号或柯西分布信号，并将其与冲激信号进行卷积，生成探测器信号的冲激余弦平方成形信号或探测器信号的冲激柯西分布成形信号，其中，数字高斯卷积的系数由公式（20）确定，余弦平方卷积的系数由公式（21）确定，数字柯西卷积的系数由公式(22)来确定，如下所示：
Figure 940386DEST_PATH_IMAGE039
(20)
Figure 100002_DEST_PATH_IMAGE040
(21)
Figure 590810DEST_PATH_IMAGE041
(22);
其中，n为采集的输入信号点序列，H为公式对应信号的半宽度。这里的C1[n]、C2[n]和C3[n]相当于我们常用的f（x），代表函数式，n是变量。
本发明还提出一种核信号的级联冲激卷积成形装置，包括：
数据采集单元：实时采集探测器信号并通过芯片传输到ARM处理器；
冲激成形单元：将探测器信号经过级联逆系统进行冲激成形，形成冲激信号；
卷积成形单元：将目标信号与冲激成形单元形成的冲激信号进行卷积，生成级联冲激卷积成形信号；
TCP/IP网络：根据不同的探测器和信号调整电路进行参数设置。
较优的，本发明所述目标信号包括：标准高斯信号、余弦平方信号、柯西分布信号和梯形信号。
相比现有技术，本发明的技术方案具有如下优点/有益效果：
1.本发明基于探测器信号冲激成形后再高斯卷积的数字高斯滤波方法，或者高斯卷积后再冲激成形的数字高斯滤波方法，针对不同的探测器信号，高斯成形的参数易于调节。
2.本发明成形后去掉了脉冲的拖尾，使其对称性更好，脉冲宽度变窄更适合于高计数率时的能谱测量。
3.本发明缩减了乘法器的规模，并且推广到了双指数以及更加复杂信号的高斯成形。
4.本发明设计的实时数字高斯成形方法相对要简单很多，与AD采样系统可以高速并行，对于闪烁探测器可以在有几十个以上乘法器的中端FPGA器件上部署，对于具有更高分辨的半导体探测器需要采用有几百个以上乘法器或者DSP的中高端FPGA器件来部署。
5.本发明实现了核信号的数字化级联冲激卷积成形滤波器，并且可以扩展到三指数信号或四指数信号进行高斯成形、梯形成形、CUSP成形、余弦平方分布成形和柯西分布成形等。
附图说明
为了更清楚地说明本发明实施方式的技术方案，下面将对实施方式中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1是本发明连续阶跃信号与高斯信号一阶导模拟(c=60/√2)(输入信号)示意图。
图2是本发明连续阶跃信号与高斯信号一阶导卷积模拟（输出信号）示意图。
图3是本发明基于级联逆系统的冲激脉冲数字成形示意图。
图4是本发明探测器信号双指数冲激成形模拟示意图。
图5是本发明探测器冲激信号高斯成形模拟示意图。
图6是本发明高斯短上升沿模拟示意图。
图7是本发明高斯长上升沿模拟示意图。
图8是本发明COS平方短上升沿模拟示意图。
图9是本发明COS平方长上升沿模拟示意图。
图10是本发明探测器冲激信号直接卷积高斯成形模拟示意图。
图11是本发明阶跃信号与余弦平方分布信号一阶导卷积模拟(H（半高宽）=128个采样点)示意图。
图12是本发明连续阶跃信号与柯西分布信号一阶导卷积模拟示意图。
图13是本发明单指数信号数字梯形成形卷积信号模拟示意图。
图14是本发明单指数信号数字梯形成形卷积信号模拟示意图。
图15是本发明双指数信号数字梯形成形卷积信号模拟示意图。
图16是本发明高斯成形探测器信号测试(NaI探测器, 1.65μs数字脉冲宽度H=16,65点Gauss)示意图。
图17是本发明高斯成形能谱测试(NaI探测器 Cs-137 FWHM:6.81% ,1.65μs数字脉冲宽度)示意图。
图18是本发明高斯成形能谱测试（Cs-137+K-40+Th-232）示意图。
图19为本发明实施例1的基于级联逆系统的卷积冲激成形示意图。
具体实施方式
为使本发明目的、技术方案和优点更加清楚，下面对本发明实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式是本发明的一部分实施方式，而不是全部的实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。因此，以下提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施方式。
因为采用高斯信号与探测器信号直接卷积，处理后的信号拖尾很长，图1是采用高斯信号一阶导与阶跃信号的卷积模拟，图2是模拟结果。由此可见，可以采用把输入信号变换成连续阶跃信号来处理。
实施例1：
如图所示，本发明提出一种核信号的级联冲激卷积成形方法，包括：
S1：通过探测器获取到探测器信号；
S2：如图19所示，S2：将探测器信号作为输入信号，将输入信号经过多级级联成形系统，输入信号先通过级联卷积系统与目标信号进行卷积，然后再通过级联逆系统进行冲激成形，生成用于分析的级联冲激卷积信号，并得到级联冲激卷积信号的函数表达式；S2中：
所述目标信号包括：标准高斯信号、余弦平方信号、柯西分布信号和梯形信号，将冲激信号与目标信号进行卷积，然后再通过级联逆系统进行冲激成形，生成探测器信号的级联冲激卷积信号，其得到函数表达式如公式（23）所示：
Figure 249325DEST_PATH_IMAGE001
*h[n]（23）；
其中，
Figure 414727DEST_PATH_IMAGE002
表示级联冲激卷积信号的函数表达式，
Figure 184100DEST_PATH_IMAGE003
表示输入信号的函数表达式，
Figure 485768DEST_PATH_IMAGE004
表示标准高斯信号表达式，h[n]表示双指数冲激脉冲成形系统函数表达式，n为采集的输入信号点序列。
如图1和图2所示，将高斯信号与连续阶跃信号进行卷积成形生成便于分析的脉宽很窄的高斯形态脉冲信号，由于连续阶跃信号的微分是冲激信号，因此先将探测器信号变换为冲激信号。图3为本发明的级联逆系统，本发明通过级联逆系统将双指数信号转换成冲激信号为现有技术。
S3:然后基于多级级联成形系统的级联卷积支持交换卷积顺序的特性，将输入信号经过级联逆系统进行冲激成形，形成冲激信号，并得到输入信号冲激成形的系统函数表达式；本发明S3具体为：
S3.1：定义输入信号
Figure 813719DEST_PATH_IMAGE005
为双指数信号，将输入信号输入第一级INV_RC系统，输出单指数衰减信号
Figure 466417DEST_PATH_IMAGE006
，所述输入信号为
Figure 773902DEST_PATH_IMAGE005
由式（1）所示，所述单指数衰减信号
Figure 195656DEST_PATH_IMAGE006
由式（2）所示：
Figure 195973DEST_PATH_IMAGE007
(1)
Figure 70388DEST_PATH_IMAGE008
(2)；
其中，m和M是双指数信号的系统参数，n为采集的输入信号点序列；
INV是inverse反向，RC指的是电路中的电阻R和电容C，连起来指的是逆RC,就是通过逆向运算，把电路中的RC影响去掉。
S3.2：将所述单指数衰减信号
Figure 447143DEST_PATH_IMAGE006
输入第二级INV_RC系统，输出得到冲激响应信号
Figure 723403DEST_PATH_IMAGE009
，冲激响应信号
Figure 691359DEST_PATH_IMAGE009
由式（3）所示：
Figure 990753DEST_PATH_IMAGE010
(3)
S3.3：根据INV_RC算子的数字解推导式(2)和式(3)，得到公式(4)和式(5)，如下所示：
Figure 702358DEST_PATH_IMAGE011
(4)
Figure 36387DEST_PATH_IMAGE012
(5)
将公式（4）代入公式（5），得到由输入信号
Figure 440823DEST_PATH_IMAGE005
通过级联逆系统变为冲激响应信号
Figure 227514DEST_PATH_IMAGE009
的数字变换表达式（6），如下所示：
Figure 742809DEST_PATH_IMAGE013
(6)
本发明成为单指数衰减信号后，幅度缩小到原来的1/M, 单指数信号经过第二级INV_RC系统后成为冲激信号，幅度扩大到原来的M倍,因此幅度必须缩小到1/M。
对公式（6）两边求差分整理得式（7）、式（8）和式（9），如下所示：
Figure 728082DEST_PATH_IMAGE014
(7)
Figure 241103DEST_PATH_IMAGE015
(8)
Figure 100002_DEST_PATH_IMAGE042
(9)
S3.4：对公式（8）进行整理，可得公式（10）的形式，如下所示：
Figure 296782DEST_PATH_IMAGE019
Figure 615768DEST_PATH_IMAGE020
(10)
通过公式（10）得到双指数信号冲激成形的系统函数表达式如式（11）所示：
h1[n]=
Figure 658811DEST_PATH_IMAGE021
(11)。
当定义输入信号
Figure 405050DEST_PATH_IMAGE022
为单指数信号时，m=0，通过公式（11）得到单指数信号冲激成形的系统函数表达式如式（12）所示：
Figure 900753DEST_PATH_IMAGE043
(12)。
当定义输入信号为带反冲的双指数信号时，将输入信号
Figure 23430DEST_PATH_IMAGE005
输入第一级INV_RC系统，第一级INV_RC系统的输出信号为
Figure 717716DEST_PATH_IMAGE006
，所述输入信号为
Figure 572540DEST_PATH_IMAGE005
由式（13）所示，所述输出信号
Figure 883436DEST_PATH_IMAGE006
由式（14）所示：
Figure 481907DEST_PATH_IMAGE024
(13)
Figure 30700DEST_PATH_IMAGE025
(14)
可知式（14）为第一级INV_RC系统的函数表达式，通过式（11）的得到过程推导出带反冲的双指数信号冲激成形的系统函数表达式如式（15）所示：
Figure 100002_DEST_PATH_IMAGE044
Figure 587583DEST_PATH_IMAGE028
(15)。
S4：将冲激信号通过级联卷积系统与目标信号进行卷积，生成级联冲激卷积成形信号，并得到多级级联成形系统的函数表达式。图4是探测器信号的双指数冲激成形信号的效果图，图5是利用图4中的冲激成形数据再与高斯信号卷积成形的模拟图。
本发明S4中：
所述目标信号包括：标准高斯信号、余弦平方信号、柯西分布信号和梯形信号，将冲激信号通过级联卷积系统与目标信号进行卷积，生成级联冲激卷积成形信号，其函数表达式如式（16）所示：
Figure 120196DEST_PATH_IMAGE029
（16）；
其中，
Figure 256779DEST_PATH_IMAGE030
表示级联冲激卷积成形信号的函数表达式，
Figure 191237DEST_PATH_IMAGE031
表示双指数冲激脉冲成形系统函数表达式，
Figure 653443DEST_PATH_IMAGE004
表示标准高斯信号。
当目标信号为标准高斯信号时，将冲激信号与标准高斯信号进行卷积，生成输入信号的级联冲激卷积成形信号，分别将式（11）、式（12）和式（15）代入式（16），然后基于多级级联成形系统的级联卷积支持交换卷积顺序的特性，得到多级级联成形系统的函数表达式如式（17）、式（18）和式（19）所示：
Figure 407772DEST_PATH_IMAGE032
=
Figure 410363DEST_PATH_IMAGE033
（17）
Figure 635546DEST_PATH_IMAGE034
=
Figure 65390DEST_PATH_IMAGE035
（18）
Figure 510278DEST_PATH_IMAGE036
=
Figure 316560DEST_PATH_IMAGE037
（19）；
其中，
Figure 897714DEST_PATH_IMAGE038
表示目标信号的函数表达式，此处目标信号采用标准高斯信号。
较优的，本发明将标准高斯信号替换成余弦平方信号或柯西分布信号，并将其与冲激信号进行卷积，生成探测器信号的冲激余弦平方成形信号或探测器信号的冲激柯西分布成形信号，其中，数字高斯卷积的系数由公式（20）确定，余弦平方卷积的系数由公式（21）确定，数字柯西卷积的系数由公式(22)来确定，如下所示：
Figure 764039DEST_PATH_IMAGE039
(20)
Figure 430643DEST_PATH_IMAGE040
(21)
Figure 775037DEST_PATH_IMAGE041
(22);
其中，H为公式对应信号的半宽度。
图6和图7为式（23）卷积得到的不同上升时间双指数信号-高斯卷积成形信号图。图8和图9为式（23）卷积生成的不同上升时间双指数信号-COS平方分布成形卷积信号。图10为双指数信号的冲激高斯卷积成形信号的模拟。
图11为阶跃信号与余弦平方分布信号一阶导卷积模拟(H（半高宽）=128个采样点)。
图12为连续阶跃信号与柯西分布信号一阶导卷积模拟。
本发明按照同样的原理将梯形信号与单指数冲激系统信号卷积构建出单指数信号数字梯形成形卷积信号，图14模拟的是输入信号的时间常数与梯形上升的点数相等时的效果，图13是单指数信号变换到冲激的卷积信号，与梯形信号做卷积的结果，即采用单指数信号作为卷积信号也能实现正常的单指数信号梯形成形，按照同样的原理将梯形信号与双指数冲激系统信号卷积构建出双指数信号数字梯形成形卷积信号如图15。
本发明还提出一种核信号的级联冲激卷积成形装置，包括：
数据采集单元：实时采集探测器信号并通过芯片传输到ARM处理器；
冲激成形单元：将探测器信号经过级联逆系统进行冲激成形，形成冲激信号；
卷积成形单元：将目标信号与冲激成形单元形成的冲激信号进行卷积，生成级联冲激卷积成形信号；所述目标信号包括：标准高斯信号、余弦平方信号、柯西分布信号和梯形信号。
TCP/IP网络：根据不同的探测器和信号调整电路进行参数设置。
对构建的核信号级联冲激卷积数字滤波器（即核信号的级联冲激卷积成形装置）进行测试，图16为取高斯信号半宽度为16（1.65us）共65点NaI探测器信号数字高斯成形，由图可见信号对称性很好，与高斯信号的近似度很高，噪声很小；图17为高斯成形Cs-137能谱测试(Φ75╳100 NaI探测器，1.65μs数字脉冲宽度。）FWHM=6.81%；且图17中的低能部分出来了半个峰（之前都作为噪声去掉了，看不到。说明该方法对信号和噪声的分辨能力很强。），一般的都看不见。分辨率能够提高0.1-0.2左右。图18为Cs-137+K-40+Th-232能谱测试，能谱线性很好。
以上仅是本发明的优选实施方式，应当指出的是，上述优选实施方式不应视为对本发明的限制，本发明的保护范围应当以权利要求所限定的范围为准。对于本技术领域的普通技术人员来说，在不脱离本发明的精神和范围内，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (9)

1.一种核信号的级联冲激卷积成形方法，其特征在于，包括：
S1：通过探测器获取到探测器信号；
S2：将探测器信号作为输入信号，将输入信号经过多级级联成形系统，输入信号先通过级联卷积系统与目标信号进行卷积，然后再通过级联逆系统进行冲激成形，生成用于分析的级联冲激卷积信号，并得到级联冲激卷积信号的函数表达式；
S3:然后基于多级级联成形系统的级联卷积支持交换卷积顺序的特性，将输入信号经过级联逆系统进行冲激成形，形成冲激信号，并得到输入信号冲激成形的系统函数表达式；
S3具体为：
S3.1：当输入信号
Figure DEST_PATH_IMAGE001
为双指数信号时，将输入信号
Figure 825039DEST_PATH_IMAGE001
输入第一级INV_RC系统，输出单指数衰减信号
Figure DEST_PATH_IMAGE002
，所述输入信号为
Figure 273338DEST_PATH_IMAGE001
由式（1）所示，所述单指数衰减信号
Figure 584234DEST_PATH_IMAGE002
由式（2）所示：
Figure DEST_PATH_IMAGE003
(1)
Figure DEST_PATH_IMAGE004
(2)；
其中，m和M是双指数信号的系统参数，n为采集的输入信号点序列，INV_RC表示逆RC,INV是inverse反向，RC指的是电路中的电阻R和电容C，即通过逆向运算，把电路中的RC影响去掉；
S3.2：将所述单指数衰减信号
Figure 245022DEST_PATH_IMAGE002
输入第二级INV_RC系统，输出得到冲激响应信号
Figure DEST_PATH_IMAGE005
，冲激响应信号
Figure 590553DEST_PATH_IMAGE005
由式（3）所示：
Figure DEST_PATH_IMAGE006
(3)
S3.3：根据INV_RC算子的数字解推导式(2)和式(3)，得到公式(4)和式(5)，如下所示：
Figure DEST_PATH_IMAGE007
(4)
Figure DEST_PATH_IMAGE008
(5)
将公式（4）代入公式（5），得到由输入信号
Figure 475332DEST_PATH_IMAGE001
通过级联逆系统变为冲激响应信号
Figure 40568DEST_PATH_IMAGE005
的数字变换表达式（6），如下所示：
Figure DEST_PATH_IMAGE009
(6)
对公式（6）两边求差分整理得式（7）、式（8）和式（9），如下所示：
Figure DEST_PATH_IMAGE010
(7)
Figure DEST_PATH_IMAGE011
(8)
Figure DEST_PATH_IMAGE013
(9)
S3.4：对公式（8）进行整理，可得公式（10）的形式，如下所示：
Figure DEST_PATH_IMAGE015
Figure DEST_PATH_IMAGE017
(10)
通过公式（10）得到双指数信号冲激成形的系统函数表达式如式（11）所示：
h1[n]=
Figure DEST_PATH_IMAGE019
(11)；
S4：将冲激信号通过级联卷积系统与目标信号进行卷积，生成级联冲激卷积成形信号，并得到多级级联成形系统的函数表达式。
2.根据权利要求1所述的一种核信号的级联冲激卷积成形方法，其特征在于，S2中：
所述目标信号包括：标准高斯信号、余弦平方信号、柯西分布信号和梯形信号，将输入信号与目标信号进行卷积，然后再通过级联逆系统进行冲激成形，生成探测器信号的级联冲激卷积信号，其得到函数表达式如公式（23）所示：
Figure DEST_PATH_IMAGE020
*h[n]（23）；
其中，
Figure DEST_PATH_IMAGE021
表示级联冲激卷积信号的函数表达式，
Figure DEST_PATH_IMAGE022
表示输入信号的函数表达式，
Figure DEST_PATH_IMAGE023
表示目标信号表达式，h[n]表示冲激成形系统函数表达式，n为采集的输入信号点序列。
3.根据权利要求1所述的一种核信号的级联冲激卷积成形方法，其特征在于，当输入信号
Figure DEST_PATH_IMAGE024
为单指数信号时，m=0，通过公式（11）得到单指数信号冲激成形的系统函数表达式如式（12）所示：
Figure DEST_PATH_IMAGE025
(12)。
4.根据权利要求3所述的一种核信号的级联冲激卷积成形方法，其特征在于，S3具体为：
S3.1：当输入信号
Figure 990201DEST_PATH_IMAGE001
为带反冲的双指数信号时，将输入信号
Figure 393500DEST_PATH_IMAGE001
输入第一级INV_RC系统，第一级INV_RC系统的输出信号为
Figure 183601DEST_PATH_IMAGE002
，所述输入信号为
Figure 672352DEST_PATH_IMAGE001
由式（13）所示，所述输出信号
Figure 206101DEST_PATH_IMAGE002
由式（14）所示：
Figure DEST_PATH_IMAGE026
(13)
Figure DEST_PATH_IMAGE028
(14)
S3.2：可知式（14）为第一级INV_RC系统的函数表达式，通过式（11）的得到过程推导出带反冲的双指数信号冲激成形的系统函数表达式如式（15）所示：
Figure DEST_PATH_IMAGE030
Figure DEST_PATH_IMAGE032
(15)。
5.根据权利要求4所述的一种核信号的级联冲激卷积成形方法，其特征在于，S4中：
所述目标信号包括：标准高斯信号、余弦平方信号、柯西分布信号和梯形信号，将冲激信号通过级联卷积系统与目标信号进行卷积，生成级联冲激卷积成形信号，其函数表达式如式（16）所示：
Figure DEST_PATH_IMAGE034
（16）；
其中，
Figure DEST_PATH_IMAGE035
表示级联冲激卷积成形信号的函数表达式，
Figure DEST_PATH_IMAGE036
表示冲激成形系统函数表达式，
Figure 588541DEST_PATH_IMAGE023
表示目标信号表达式。
6.根据权利要求5所述的一种核信号的级联冲激卷积成形方法，其特征在于，
当目标信号为标准高斯信号时，将冲激信号与标准高斯信号进行卷积，生成输入信号的级联冲激卷积成形信号，分别将式（11）、式（12）和式（15）代入式（16），然后基于多级级联成形系统的级联卷积支持交换卷积顺序的特性，得到多级级联成形系统的函数表达式如式（17）、式（18）和式（19）所示：
Figure DEST_PATH_IMAGE038
=
Figure DEST_PATH_IMAGE040
（17）
Figure DEST_PATH_IMAGE042
=
Figure DEST_PATH_IMAGE044
（18）
Figure DEST_PATH_IMAGE046
=
Figure DEST_PATH_IMAGE048
（19）；
其中，
Figure DEST_PATH_IMAGE049
表示目标信号的函数表达式，此处目标信号采用标准高斯信号。
7.根据权利要求6所述的一种核信号的级联冲激卷积成形方法，其特征在于，将标准高斯信号替换成余弦平方信号或柯西分布信号，并将其与冲激信号进行卷积，生成探测器信号的冲激余弦平方成形信号或探测器信号的冲激柯西分布成形信号，其中，数字高斯卷积的系数由公式（20）确定，余弦平方卷积的系数由公式（21）确定，数字柯西卷积的系数由公式(22)来确定，如下所示：
Figure DEST_PATH_IMAGE051
(20)
Figure DEST_PATH_IMAGE053
(21)
Figure DEST_PATH_IMAGE055
(22);
其中，n为采集的输入信号点序列，H为公式对应信号的半宽度。
8.一种核信号的级联冲激卷积成形装置，其特征在于，包括：
数据采集单元：实时采集探测器信号并通过芯片传输到ARM处理器；
冲激成形单元：将探测器信号经过级联逆系统进行冲激成形，形成冲激信号；
具体为：
步骤1：当输入信号
Figure 695082DEST_PATH_IMAGE001
为双指数信号，将输入信号
Figure 733445DEST_PATH_IMAGE001
输入第一级INV_RC系统，输出单指数衰减信号
Figure 742990DEST_PATH_IMAGE002
，所述输入信号为
Figure 183198DEST_PATH_IMAGE001
由式（1）所示，所述单指数衰减信号
Figure 252785DEST_PATH_IMAGE002
由式（2）所示：
Figure 778445DEST_PATH_IMAGE003
(1)
Figure 591680DEST_PATH_IMAGE004
(2)；
其中，m和M是双指数信号的系统参数，n为采集的输入信号点序列，INV_RC表示逆RC,INV是inverse反向，RC指的是电路中的电阻R和电容C，即通过逆向运算，把电路中的RC影响去掉；
步骤2：将所述单指数衰减信号
Figure 387860DEST_PATH_IMAGE002
输入第二级INV_RC系统，输出得到冲激响应信号
Figure 893927DEST_PATH_IMAGE005
，冲激响应信号
Figure 641304DEST_PATH_IMAGE005
由式（3）所示：
Figure 258230DEST_PATH_IMAGE006
(3)
步骤3：根据INV_RC算子的数字解推导式(2)和式(3)，得到公式(4)和式(5)，如下所示：
Figure 407451DEST_PATH_IMAGE007
(4)
Figure 84420DEST_PATH_IMAGE008
(5)
将公式（4）代入公式（5），得到由输入信号
Figure 319093DEST_PATH_IMAGE001
通过级联逆系统变为冲激响应信号
Figure 474130DEST_PATH_IMAGE005
的数字变换表达式（6），如下所示：
Figure 743438DEST_PATH_IMAGE009
(6)
对公式（6）两边求差分整理得式（7）、式（8）和式（9），如下所示：
Figure 591308DEST_PATH_IMAGE010
(7)
Figure 313276DEST_PATH_IMAGE011
(8)
Figure 272005DEST_PATH_IMAGE013
(9)
步骤4：对公式（8）进行整理，可得公式（10）的形式，如下所示：
Figure DEST_PATH_IMAGE015A
Figure 661398DEST_PATH_IMAGE017
(10)
通过公式（10）得到双指数信号冲激成形的系统函数表达式如式（11）所示：
h1[n]=
Figure 680170DEST_PATH_IMAGE019
(11)；
卷积成形单元：将目标信号与冲激成形单元形成的冲激信号进行卷积，生成级联冲激卷积成形信号；
当目标信号为标准高斯信号时，基于多级级联成形系统的级联卷积支持交换卷积顺序的特性，得到多级级联成形系统的函数表达式如式（17）所示：
Figure 125320DEST_PATH_IMAGE038
=
Figure 950056DEST_PATH_IMAGE040
（17）；
其中，
Figure 866060DEST_PATH_IMAGE049
表示目标信号的函数表达式，此处目标信号采用标准高斯信号；
TCP/IP网络：根据不同的探测器和信号调整电路进行参数设置。
9.根据权利要求8所述的一种核信号的级联冲激卷积成形装置，其特征在于，所述目标信号包括：标准高斯信号、余弦平方信号、柯西分布信号和梯形信号。
CN202210196100.0A 2022-03-02 2022-03-02 一种核信号的级联冲激卷积成形方法和装置 Active CN114252899B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210196100.0A CN114252899B (zh) 2022-03-02 2022-03-02 一种核信号的级联冲激卷积成形方法和装置
US18/087,853 US20230281266A1 (en) 2022-03-02 2022-12-23 Cascaded impulse convolution shaping method and apparatus for nuclear signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210196100.0A CN114252899B (zh) 2022-03-02 2022-03-02 一种核信号的级联冲激卷积成形方法和装置

Publications (2)

Publication Number Publication Date
CN114252899A CN114252899A (zh) 2022-03-29
CN114252899B true CN114252899B (zh) 2022-05-20

Family

ID=80797318

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210196100.0A Active CN114252899B (zh) 2022-03-02 2022-03-02 一种核信号的级联冲激卷积成形方法和装置

Country Status (2)

Country Link
US (1) US20230281266A1 (zh)
CN (1) CN114252899B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116774273B (zh) * 2023-08-16 2023-10-31 成都理工大学 一种核脉冲信号展开为冲激信号的方法及系统

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2870603A1 (fr) * 2004-05-19 2005-11-25 Commissariat Energie Atomique Mesure et traitement d'un signal comprenant des empilements d'impulsions elementaires
WO2009121130A1 (en) * 2008-03-31 2009-10-08 Southern Innovation International Pty Ltd Screening method and apparatus
CN103837884A (zh) * 2014-02-26 2014-06-04 成都理工大学 基于时域分析的数字核脉冲信号梯形成形算法
CN103941280A (zh) * 2013-11-29 2014-07-23 成都理工大学 基于冲激响应不变法的数字核脉冲高斯成形方法
CN104242877A (zh) * 2014-07-16 2014-12-24 成都理工大学 基于模拟cr-rc电路的核脉冲信号数字高斯成形方法
CN106019357A (zh) * 2016-05-23 2016-10-12 成都理工大学 基于rc逆变换的核脉冲信号处理方法
CN106772545A (zh) * 2017-02-09 2017-05-31 成都理工大学 一种采用尖顶成形算法的数字脉冲幅度分析器
CN106772541A (zh) * 2016-12-16 2017-05-31 成都理工大学 探测器输出信号反卷积处理方法研究
CN107300713A (zh) * 2017-05-23 2017-10-27 西安交通大学 基于反卷积的自给能中子探测器延迟效应消除方法
CN109507709A (zh) * 2018-11-21 2019-03-22 成都理工大学 核脉冲信号处理方法及装置
CN109873622A (zh) * 2019-03-27 2019-06-11 成都理工大学 核信号数字处理方法、装置及数字化核仪器
CN109885903A (zh) * 2019-01-29 2019-06-14 吉林大学 一种基于模型的地面核磁共振信号尖峰噪声去除方法
CN111404542A (zh) * 2020-03-31 2020-07-10 四川新先达测控技术有限公司 一种双指数核信号计数方法及装置
CN111934690A (zh) * 2020-08-04 2020-11-13 东华理工大学 基于自适应压缩感知的核信号重构新方法
CN112462675A (zh) * 2021-01-27 2021-03-09 泛华检测技术有限公司 一种级联型双指数核脉冲信号发生装置及其控制方法
CN112462676A (zh) * 2021-01-27 2021-03-09 泛华检测技术有限公司 一种可模拟重叠核脉冲信号发生装置及其控制方法
CN113011287A (zh) * 2021-03-02 2021-06-22 四川新先达测控技术有限公司 一种提高高斯信号分析精确度的方法及系统
CN113189634A (zh) * 2021-03-02 2021-07-30 四川新先达测控技术有限公司 一种类高斯成形方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018091434A1 (en) * 2016-11-15 2018-05-24 Thermo Fisher Scientific Messtechnik Gmbh System and method of neutron radiation detection
CN208384130U (zh) * 2018-07-21 2019-01-15 东华理工大学 一种基于fpga/stm32的数字多道谱仪
CN111273336B (zh) * 2020-02-13 2022-02-11 东华理工大学 一种数字核脉冲信号高斯成形方法
CN113568032B (zh) * 2021-07-02 2024-02-23 四川新先达测控技术有限公司 一种基于z变换的负指数核脉冲信号处理方法及系统

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2870603A1 (fr) * 2004-05-19 2005-11-25 Commissariat Energie Atomique Mesure et traitement d'un signal comprenant des empilements d'impulsions elementaires
WO2009121130A1 (en) * 2008-03-31 2009-10-08 Southern Innovation International Pty Ltd Screening method and apparatus
CN103941280A (zh) * 2013-11-29 2014-07-23 成都理工大学 基于冲激响应不变法的数字核脉冲高斯成形方法
CN103837884A (zh) * 2014-02-26 2014-06-04 成都理工大学 基于时域分析的数字核脉冲信号梯形成形算法
CN104242877A (zh) * 2014-07-16 2014-12-24 成都理工大学 基于模拟cr-rc电路的核脉冲信号数字高斯成形方法
CN106019357A (zh) * 2016-05-23 2016-10-12 成都理工大学 基于rc逆变换的核脉冲信号处理方法
CN106772541A (zh) * 2016-12-16 2017-05-31 成都理工大学 探测器输出信号反卷积处理方法研究
CN106772545A (zh) * 2017-02-09 2017-05-31 成都理工大学 一种采用尖顶成形算法的数字脉冲幅度分析器
CN107300713A (zh) * 2017-05-23 2017-10-27 西安交通大学 基于反卷积的自给能中子探测器延迟效应消除方法
CN109507709A (zh) * 2018-11-21 2019-03-22 成都理工大学 核脉冲信号处理方法及装置
CN109885903A (zh) * 2019-01-29 2019-06-14 吉林大学 一种基于模型的地面核磁共振信号尖峰噪声去除方法
CN109873622A (zh) * 2019-03-27 2019-06-11 成都理工大学 核信号数字处理方法、装置及数字化核仪器
CN111404542A (zh) * 2020-03-31 2020-07-10 四川新先达测控技术有限公司 一种双指数核信号计数方法及装置
CN111934690A (zh) * 2020-08-04 2020-11-13 东华理工大学 基于自适应压缩感知的核信号重构新方法
CN112462675A (zh) * 2021-01-27 2021-03-09 泛华检测技术有限公司 一种级联型双指数核脉冲信号发生装置及其控制方法
CN112462676A (zh) * 2021-01-27 2021-03-09 泛华检测技术有限公司 一种可模拟重叠核脉冲信号发生装置及其控制方法
CN113011287A (zh) * 2021-03-02 2021-06-22 四川新先达测控技术有限公司 一种提高高斯信号分析精确度的方法及系统
CN113189634A (zh) * 2021-03-02 2021-07-30 四川新先达测控技术有限公司 一种类高斯成形方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Counting-loss correction for X-ray spectroscopy using unit impulse pulse shaping;Xu Hong et al.;《Journal of Synchrotron Radiation》;20180301;第505-513页 *
Peak tailing cancellation techniques for digital CR-(RC)n filter;Xu Hong et al.;《Applied Radiation and Isotopes》;20201011;第109471-1-6页 *
Study of Digital Gaussian Shaping for a Semiconductor X-ray Detector in Real Time;Zhou Jianbin et al;《IEEE 2011 10th International Conference on Electronic Measurement & Instruments》;20111010;第163-165页 *
基于Z变换的C-R逆算子分析及其应用研究;金幼贤 等;《核技术》;20210730;第44卷(第7期);第070403-1-7页 *
数字化多道中双指数脉冲梯形成形算法研究;孙超 等;《核电子学与探测技术》;20170730;第37卷(第7期);第752-756页 *
核脉冲实时读出系统及数字成形算法研究;王明;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技II辑》;20210415(第4期);第C040-41页 *
高计数率下的实时梯形成形算法;黄宇雁 等;《清华大学学报（自然科学版）》;20170515;第57卷(第5期);第521-524页 *

Also Published As

Publication number Publication date
CN114252899A (zh) 2022-03-29
US20230281266A1 (en) 2023-09-07

Similar Documents

Publication Publication Date Title
Stein et al. X-ray detectors with digitized preamplifiers
CN110609050B (zh) 一种x射线荧光光谱峰拖尾消除方法及系统
Saxena et al. Investigation of FPGA-based real-time adaptive digital pulse shaping for high-count-rate applications
CN114252899B (zh) 一种核信号的级联冲激卷积成形方法和装置
Regadío et al. Implementation of a real-time adaptive digital shaping for nuclear spectroscopy
CN111600584B (zh) 一种核脉冲信号处理方法及系统
Liu et al. Implementation of real-time digital CR–RCm shaping filter on FPGA for gamma-ray spectroscopy
CN113189634B (zh) 一种类高斯成形方法
CN113568032B (zh) 一种基于z变换的负指数核脉冲信号处理方法及系统
CN111697952B (zh) 一种基于数字pzc系统调节脉冲宽度的方法及系统
Arsenault et al. Comparison of techniques for extracting signals from a strong background
CN106772540A (zh) 一种j‑text托卡马克装置的硬x射线通量检测系统
CN111008356B (zh) 一种基于WTSVD算法扣除背景的γ能谱集分析方法
CN112255665A (zh) 一种闪烁体发光衰减时间数字化测量系统及测量方法
CN113934359B (zh) 信号处理器、信号处理的方法及装置、可读存储介质
CN114254495B (zh) 一种基于蒙特卡洛方法的模拟能谱生成方法
Paul et al. Implementation of FPGA based real-time digital DAQ for high resolution, and high count rate nuclear spectroscopy application
KalaiPriya et al. Vlsi implementation of nonlinear variable cutoff high pass filter algorithm
CN113507279B (zh) 一种高精度星间距离观测数据的降速率滤波方法
Van den Bogaert et al. Quantification of chromatographic data using a matched filter: robustness towards noise model errors
Xie et al. A digital energy spectroscopy based on FIR filter
Raad et al. Efficient algorithms for pulse parameter estimation, pulse peak localization and pileup reduction in gamma ray spectroscopy
Jiang et al. Trapezoidal shaping algorithm based on FPGA
Li et al. Review of research on nuclear signal pulse shaping
RU211396U1 (ru) Устройство для цифровой коррекции нелинейности

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant