CN111273336B - 一种数字核脉冲信号高斯成形方法 - Google Patents
一种数字核脉冲信号高斯成形方法 Download PDFInfo
- Publication number
- CN111273336B CN111273336B CN202010091332.0A CN202010091332A CN111273336B CN 111273336 B CN111273336 B CN 111273336B CN 202010091332 A CN202010091332 A CN 202010091332A CN 111273336 B CN111273336 B CN 111273336B
- Authority
- CN
- China
- Prior art keywords
- gaussian
- digital
- pulse signal
- circuit
- forming
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 230000005540 biological transmission Effects 0.000 claims abstract description 5
- 238000007493 shaping process Methods 0.000 claims description 28
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 230000007547 defect Effects 0.000 abstract description 12
- 230000000694 effects Effects 0.000 abstract description 6
- 238000005316 response function Methods 0.000 abstract 2
- 238000010586 diagram Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 229910052770 Uranium Inorganic materials 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Measurement Of Radiation (AREA)
- Manipulation Of Pulses (AREA)
Abstract
本发明公开了一种数字核脉冲信号高斯成形方法,基于三阶Sallen‑Key电路的高斯成形的基本思想,获取其在拉普拉斯域内的传输函数,根据传输函数推导了其数值递推模型,并给出了该模型电路在时域中高斯成形的响应函数,将数字核脉冲信号输入到数值递推函数模型,可实现其高斯成形输出,通过调整成形参数k,可改变高斯成形输出的脉冲幅度与脉冲宽度,根据数值递推模型推导其在Z域中的幅频响应函数,获取其幅频响应曲线,分析高斯成形的效果;本发明克服了基于CR‑RCm电路的数字高斯成形过程中出现过零下冲、成形输出幅度小的缺点;克服基于Sallen‑Key电路数字高斯成形输出对称性不足、滤波性能不佳的缺点,为数字核脉冲信号快速、有效地高斯成形提供了一种新的实现方式。
Description
技术领域
本发明涉及核脉冲信号数字化处理技术领域,特别涉及一种数字核脉冲信号高斯成形方法。
背景技术
数字核能谱测量系统由于稳定性、灵活性、抗干扰性等方面具有明显的优越性,是当前核能谱测量与获取的主要研究热点,在核电运行与安全、铀资源勘查与采冶、环境监测与辐射防护、国土安全与核反恐、工业检测与辐射成像等领域得到了广泛应用。
核脉冲信号数字化处理的成形方法直接影响着脉冲信号的噪声抑制、堆积识别以及幅度提取等功能,进而影响所测能谱的能量分辨率等主要性能指标,高斯成形作为当前数字核脉冲信号的主要成形方法,由于其成形简单,噪声抑制与幅度提取等综合性能良好被广泛用于数字核脉冲的成形处理中。
现今常用的数字核脉冲信号高斯成形方法包括:基于CR-RCm电路的数字高斯成形与基于Sallen-Key电路的数字高斯成形两种,其中数字CR-RCm采用一级的CR微分成形,对于输入为负指数下降沿的核脉冲信号,经过一级的CR微分之后,出现过零下冲信号,影响了脉冲幅度信号的准确提取,同时m级的RC成形(为了保证信噪比,m一般取3或是4),使得成形输出的脉冲幅度过小,也会影响脉冲幅度的准确提取;一级数字Sallen-Key高斯成形输出对称性较差,滤波效果不好,多级的数字Sallen-Key高斯成形输出幅度过小,均会影响脉冲幅度的有效提取。
针对上述问题,本申请提供一种数字核脉冲信号高斯成形方法,克服了基于CR-RCm电路的数字高斯成形过程中出现过零下冲、成形输出幅度小的缺点;克服基于Sallen-Key电路数字高斯成形输出对称性不足、滤波性能不佳的缺点,为数字核脉冲信号快速、有效地高斯成形提供了一种新的实现方式。
发明内容
本发明的目的在于提供一种数字核脉冲信号高斯成形方法,克服了基于CR-RCm电路的数字高斯成形过程中出现过零下冲、成形输出幅度小的缺点;克服基于Sallen-Key电路数字高斯成形输出对称性不足、滤波性能不佳的缺点,为数字核脉冲信号快速、有效地高斯成形提供了一种新的实现方式。
本发明提供了一种数字核脉冲信号高斯成形方法,所述方法包括:
S1:根据三阶Sallen-Key滤波成形电路,得到该电路在拉普拉斯域内的传输函数,其表达式如下:
其中,R为电路中的电阻,C为电路中的电容,Vi(s)为输入信号Vi(t)的拉普拉斯变换,Vo(s)为输出信号Vo(t)的拉普拉斯变换;
S2:根据传输函数,获取在时域内的输入信号与成形输出信号之间的微分方程,输入信号Vi(t)与输出信号Vo(t)之间的微分方程表达式如下:
R1R2R3C1C2C3·Vo(t)”'+(R2R3C2C3+R1R3C2C3+R1R2C1C3+R1R3C1C3)Vo(t)”+(R3C3+R2C3+R1C3+R1C1)·Vo(t)'+Vo(t)=Vi(t)
S3:根据微分方程,建立在时域内的数值递推算式,将Vi(t)转换成数字序列x[n],Vo(t)转换成数字序列y[n],令R1=R2=R3=R,C1=C2=C3=C,则数值递推算式为:
其中,k=RC/Δt,RC为输入信号Vi(t)的时间常数,Δt为离散化间隔时间;
S4:将数字化的核脉冲信号x[n]带入数值递推算式,实现数字核脉冲信号高斯成形输出。
进一步地,根据电路的数值递推算式,确定其在Z域中的传输函数,并获取不同参数下的幅频响应曲线;
根据不同参数下的幅频响应曲线,分析基于该电路的数值递推算式的高斯成形的合理性。
进一步地,分析高斯成形的合理性,具体包括:
k值越大,成形后的脉冲越宽,峰位越低,成形结果越趋近于高斯形。
进一步地,所述方法实现电路的幅频响应为低通滤波器。
与现有技术相比,本发明具有如下显著优点:
本发明提供了一种数字核脉冲信号高斯成形方法,本发明通过改变相关参数的取值调节成形输出的脉冲宽度,脉冲幅度等信息,克服了基于CR-RCm电路的数字高斯成形输出出现过零下冲、幅度过小的缺点,克服基于Sallen-Key电路的数字高斯成形输出对称性不好,滤波效果差的缺点,为数字核脉冲信号快速、有效地高斯成形提供了一种新的实现方式。
附图说明
图1为本发明实施例提供的一种数字核脉冲信号高斯成形方法的数字核脉冲信号高斯成形流程图;
图2为本发明实施例提供的三阶Sallen-Key电路原理图;
图3为本发明实施例提供的m阶CR-RCm电路原理图;
图4为本发明实施例提供的Sallen-Key电路原理图;
图5为本发明实施例提供的不同参数下的幅频响应曲线图;
图6为本发明实施例提供的仿真数字核脉冲信号高斯成形输出图;
图7为本发明实施例提供的实际采样核脉冲信号高斯成形输出图;
图8为本发明实施例提供的高斯成形方法与现有技术的效果对比图。
具体实施方式
下面结合本发明中的附图,对本发明实施例的技术方案进行清楚、完整的描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都应属于本发明保护的范围。
为了便于理解和说明,参照图1-8,本发明提供了一种数字核脉冲信号高斯成形方法,所述方法包括:
S1:根据三阶Sallen-Key滤波成形电路,得到该电路在拉普拉斯域内的传输函数,其表达式如下:
其中,R为电路中的电阻,C为电路中的电容,Vi(s)为该输入信号Vi(t)的拉普拉斯变换,Vo(s)为该输出信号Vo(t)的拉普拉斯变换;
S2:根据传输函数,获取在时域内的输入信号与成形输出信号之间的微分方程,输入信号Vi(t)与输出信号Vo(t)之间的微分方程表达式如下:
R1R2R3C1C2C3·Vo(t)”'+(R2R3C2C3+R1R3C2C3+R1R2C1C3+R1R3C1C3)Vo(t)”+(R3C3+R2C3+R1C3+R1C1)·Vo(t)'+Vo(t)=Vi(t)
S3:根据微分方程,建立在时域内的数值递推算式,将Vi(t)转换成数字序列x[n],Vo(t)转换成数字序列y[n],令R1=R2=R3=R,C1=C2=C3=C,则数值递推算式为:
其中,k=RC/Δt,RC为输入信号Vi(t)的时间常数,Δt为离散化间隔时间(对应于ADC采样间隔时间);
S4:将输入信号带入数值递推算式,实现数字核脉冲信号高斯成形输出。
进一步地,根据电路的数值递推算式,确定其在Z域中的传输函数,并获取不同参数下的幅频响应曲线;
根据不同参数下的幅频响应曲线,分析基于该电路的数值递推算式的高斯成形的合理性。
进一步地,分析高斯成形的合理性,具体包括:
k值越大,成形后的脉冲越宽,峰位越低,成形结果越趋近于高斯形。之后并验证高斯成形效果。由数值递推算式可得出:
将z=ejw带入上式求模可得:
A=10k6+40k5+120k4+98k3+40k3+8k+1
B=4k(2k5+8k4+38k3+33k2+12k+2)
C=4k2(k4+4k3-2k2-4k-1)
D=4k3(2k3+8k2+6k+1)
E=2k3(k3+4k2+4k+1)
结合图5可以看出,该方法的数值递推函数具有低通滤波器特性,可用于核脉冲信号的数字成形处理。其中k越大,成形输出后的脉冲越宽,峰位越低,越趋近于高斯形。
进一步地,所述方法实现电路的幅频响应为低通滤波器。
如图6所示,输入信号为带有噪声的仿真核脉冲信号,在不同的成形参数下的成形输出可知,随着k的增加,成形输出幅度越小,上升沿越平缓,越趋近于高斯形,且具有较好的噪声抑制性能。
如图7所示,输入信号为实际采样获取的核脉冲信号,在不同的参数下的成形输出可知,随着k的增加,成形输出幅度越小,越趋近于高斯形,滤波效果越好,且能有效将输入信号成形为高斯形。
如图8所示,本发明中的高斯成形方法同其他2种现有的高斯成形方法对比。相同参数下的不同方法的高斯成形输出,基于CR-RC3电路的数字高斯成形由于CR的微分特性,使得成形输出出现过零下冲,而且脉冲幅度太小;基于Sallen-Key电路的数字高斯成形输出脉冲对称性不好,滤波性能较差;本发明实施的数字高斯成形输出在成形脉冲幅度、脉冲形状,滤波等方面综合性能最佳。
综上所述,本发明通过改变相关参数的取值调节成形后脉冲幅度、波形等信息,克服了基于CR-RCm电路的数字高斯成形出现过零下冲、成形输出幅度过小的缺点,克服基于Sallen-Key电路的数字高斯成形输出脉冲对称性不好、滤波性能不佳的缺点,为数字核脉冲信号快速、有效地高斯成形提供了一种新的实现方式。
以上公开的仅为本发明的几个具体实施例,但是,本发明实施例并非局限于此,任何本领域的技术人员能思之的变化都应落入本发明的保护范围。
Claims (4)
1.一种数字核脉冲信号高斯成形方法,其特征在于,所述方法包括:
S1:根据三阶Sallen-Key滤波成形电路,得到该电路在拉普拉斯域内的传输函数,其表达式如下:
其中,R为电路中的电阻,C为电路中的电容,Vi(s)为输入信号Vi(t)的拉普拉斯变换,Vo(s)为输出信号Vo(t)的拉普拉斯变换;
S2:根据传输函数,获取在时域内的输入信号与成形输出信号之间的微分方程,输入信号Vi(t)与输出信号Vo(t)之间的微分方程表达式如下:
R1R2R3C1C2C3·Vo(t)”'+(R2R3C2C3+R1R3C2C3+R1R2C1C3+R1R3C1C3)Vo(t)”+(R3C3+R2C3+R1C3+R1C1)·Vo(t)'+Vo(t)=Vi(t)
S3:根据微分方程,建立在时域内的数值递推算式,将Vi(t)转换成数字序列x[n],Vo(t)转换成数字序列y[n],令R1=R2=R3=R,C1=C2=C3=C,则数值递推
算式为:
其中,k=RC/Δt,RC为输入信号Vi(t)的时间常数,Δt为离散化间隔时间;
S4:将数字化的核脉冲信号x[n]带入数值递推算式,实现数字核脉冲信号高斯成形输出。
2.如权利要求1所述的一种数字核脉冲信号高斯成形方法,其特征在于,根据电路的数值递推算式,确定其在Z域中的传输函数,并获取不同参数下的幅频响应曲线;
根据不同参数下的幅频响应曲线,分析基于该电路的数值递推算式的高斯成形的合理性。
3.如权利要求2所述的一种数字核脉冲信号高斯成形方法,其特征在于,分析高斯成形的合理性,具体包括:
k值越大,成形后的脉冲越宽,峰位越低,成形结果越趋近于高斯形。
4.如权利要求1所述的一种数字核脉冲信号高斯成形方法,其特征在于,所述方法实现电路的幅频响应为低通滤波器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010091332.0A CN111273336B (zh) | 2020-02-13 | 2020-02-13 | 一种数字核脉冲信号高斯成形方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010091332.0A CN111273336B (zh) | 2020-02-13 | 2020-02-13 | 一种数字核脉冲信号高斯成形方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111273336A CN111273336A (zh) | 2020-06-12 |
CN111273336B true CN111273336B (zh) | 2022-02-11 |
Family
ID=70999506
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010091332.0A Expired - Fee Related CN111273336B (zh) | 2020-02-13 | 2020-02-13 | 一种数字核脉冲信号高斯成形方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111273336B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112134545B (zh) * | 2020-09-25 | 2022-04-22 | 中国核动力研究设计院 | 一种基于最佳滤波器的梯形成形方法、系统、终端及介质 |
CN112327347B (zh) * | 2020-10-29 | 2022-11-18 | 中广核久源(成都)科技有限公司 | 曲率可调的数字化核脉冲成形系统 |
CN113189634B (zh) * | 2021-03-02 | 2022-10-25 | 四川新先达测控技术有限公司 | 一种类高斯成形方法 |
CN114252899B (zh) * | 2022-03-02 | 2022-05-20 | 四川新先达测控技术有限公司 | 一种核信号的级联冲激卷积成形方法和装置 |
CN116466384B (zh) * | 2023-06-15 | 2023-11-10 | 苏州瑞派宁科技有限公司 | 闪烁脉冲的处理方法、装置、电子设备以及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110460313A (zh) * | 2019-08-22 | 2019-11-15 | 湘潭大学 | 一种用于辐射探测器的微弱信号读出电路 |
-
2020
- 2020-02-13 CN CN202010091332.0A patent/CN111273336B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110460313A (zh) * | 2019-08-22 | 2019-11-15 | 湘潭大学 | 一种用于辐射探测器的微弱信号读出电路 |
Non-Patent Citations (4)
Title |
---|
Front-end and back-end of design of Gaussian pulse shaper;vatsala Sharma et al.;《International Journalof electronics,electrical and computation system》;20170531;第6卷(第5期);第18-23页 * |
study of sallen-key digital filters in nuclear pulse signal processing;张怀强 等;《NUCL SCI TECH》;20190927;第144-151页 * |
基于FPGA的伽马射线脉冲数字高斯成形算法研究;王海涛 等;《核电子学与探测技术》;20131130;第33卷(第11期);第1334-1337页 * |
基于Sallen-Key滤波器的数字高斯成形方法的仿真;周伟等;《系统仿真学报》;20130108(第01期);第195-196页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111273336A (zh) | 2020-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111273336B (zh) | 一种数字核脉冲信号高斯成形方法 | |
US7860688B2 (en) | Signal baseline processing device and processing method thereof | |
CN105223482B (zh) | 局部放电特高频信号波形的小波分解二值去噪方法 | |
CN105785324B (zh) | 基于mgcstft的线性调频信号参数估计方法 | |
CN107817514B (zh) | 一种数字化核能谱测量系统中的脉冲阶梯成形方法 | |
CN107395157A (zh) | 基于小波变换和加权移动平均的接地网电位差滤波方法 | |
CN105842728A (zh) | 数字化核能谱测量系统中的脉冲基线估计方法 | |
CN111181634B (zh) | 一种分布式光纤振动信号快速定位方法 | |
CN117434337B (zh) | 一种宽频电压信号测量方法及系统 | |
CN108089162A (zh) | 一种脉冲干扰信号的检测与抑制方法 | |
CN112332807A (zh) | 一种弱包络信号检测方法和系统 | |
CN103823177A (zh) | 基于窗函数设计的滤波器性能检测方法和系统 | |
CN113568032B (zh) | 一种基于z变换的负指数核脉冲信号处理方法及系统 | |
CN101106487A (zh) | 一种检测网络流量异常的方法及装置 | |
CN108252753B (zh) | 一种滤除干扰信号的方法、装置和设备 | |
CN113566852A (zh) | 基线更新、相对状态检测方法和系统、电子设备 | |
CN110287853B (zh) | 一种基于小波分解的暂态信号去噪方法 | |
CN110333506B (zh) | 一种提取索力测量雷达的拉索位置参数的方法 | |
CN113218520B (zh) | 一种激光脉宽的优化神经网络提取方法 | |
CN102546116B (zh) | 基于压缩感知的跳频信号盲检测方法 | |
CN113988131A (zh) | 一种基于模拟量的数据滤波方法、系统、设备及存储介质 | |
D’Orazio et al. | Measuring reverberation time using preprocessed energy detection | |
CN110598617B (zh) | 煤矿安全监控系统瓦斯异常数据辨识与重构方法 | |
CN110582697B (zh) | 用于测量和跟踪流体中组分的数量或浓度的设备和方法 | |
Biroth et al. | A low-noise and fast pre-amplifier and readout system for SiPMs |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220211 |