CN106019357B - 基于rc逆变换的核脉冲信号处理方法 - Google Patents
基于rc逆变换的核脉冲信号处理方法 Download PDFInfo
- Publication number
- CN106019357B CN106019357B CN201610347460.0A CN201610347460A CN106019357B CN 106019357 B CN106019357 B CN 106019357B CN 201610347460 A CN201610347460 A CN 201610347460A CN 106019357 B CN106019357 B CN 106019357B
- Authority
- CN
- China
- Prior art keywords
- pulse signal
- core pulse
- negative exponent
- signal
- fingers number
- 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
Links
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)
Abstract
本发明提供一种基于RC逆变换的核脉冲信号处理方法,包括以下步骤:S1、将前置放大器电路中对负指数核脉冲信号有影响的分布电容和分布电阻等效为RC电路;S2、通过RC电路微分数值分析方法得到负指数核脉冲信号与双指数核脉冲信号的函数表达式;S3、通过ADC采样得到双指数核脉冲信号,利用双指数核脉冲信号通过步骤S2中函数表达式计算得到恢复上升沿的负指数核脉冲信号;S4、取恢复上升沿的负指数核脉冲信号的上升沿和双指数核脉冲信号的下降沿构成新的核脉冲信号。本发明消除分布电容和分布电阻对探测器输出信号上升沿的影响;对恢复上升沿后的核脉冲信号做梯形成形数字处理,改善了核脉冲信号滤波成形效果,扩大了系统的电压测量范围。
Description
技术领域
本发明涉及信号处理技术领域,具体涉及一种基于RC逆变换的核脉冲信号处理方法。
背景技术
放射性测量仪器中探测器和前置电路输出的脉冲信号是一个基于时间常数的指数衰减信号,经过FPGA电路处理后,将该信号滤波成形为特定形状,再进行多道脉冲幅度分析,最后得出能量谱,输出至上位机显示测量结果。
在进行核脉冲信号滤波成形之前,理论上无损耗的采集信号是一个标准的负指数核脉冲信号,但是由于前置放大器电路中分布电容和分布电阻的影响,信号经过前置放大器后变为一个具有上升沿的双指数核脉冲信号,如果直接采用该双指数核脉冲信号进行后续处理,会影响核脉冲信号的滤波成形效果。
发明内容
本发明的目的在于针对现有技术的不足,提供一种基于RC逆变换的核脉冲信号处理方法,该方法可以得到一个新的恢复上升沿的核脉冲信号,从而改善核脉冲信号滤波成形效果。
为达到上述要求,本发明采取的技术方案是:提供一种基于RC逆变换的核脉冲信号处理方法,该方法包括以下步骤:
S1、将前置放大器电路中对负指数核脉冲信号有影响的分布电容和分布电阻等效为RC电路;
S2、通过RC电路微分数值分析方法得到负指数核脉冲信号与双指数核脉冲信号的函数表达式;
S3、通过ADC采样获取前置放大器输出的双指数核脉冲信号,并根据双指数核脉冲信号和函数表达式计算得到负指数核脉冲信号;
S4、采用计算得到的所述负指数核脉冲信号的上升沿和双指数核脉冲信号的下降沿构建新的核脉冲信号。
该基于RC逆变换的核脉冲信号处理方法具有的优点如下:将ADC采样得到的双指数核脉冲信号进行RC电路逆运算,得到恢复上升沿的负指数核脉冲信号,消除分布电容和分布电阻对探测器输出信号上升沿的影响;对恢复上升沿后的核脉冲信号做梯形成形数字处理,改善了核脉冲信号滤波成形效果,扩大了系统的电压测量范围。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,在这些附图中使用相同的参考标号来表示相同或相似的部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1示出了根据本申请一个实施例的流程示意图;
图2示出了根据本申请一个实施例的等效RC电路的结构示意图;
图3示出了根据本申请一个实施例的ADC采集的双指数核脉冲信号;
图4示出了根据本申请一个实施例的恢复上升沿的负指数核脉冲信号;
图5示出了根据本申请一个实施例构成的新的核脉冲信号;
图6为恢复上升沿的核脉冲信号滤波成形效果和没有做逆运算滤波成形效果的对比图;
图7为恢复上升沿的负指数核脉冲信号滤波后形成的能谱图;
图8为没有做逆运算的双指数核脉冲信号滤波后形成的能谱图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,以下结合附图及具体实施例,对本申请作进一步地详细说明。
在以下描述中,对“一个实施例”、“实施例”、“一个示例”、“示例”等等的引用表明如此描述的实施例或示例可以包括特定特征、结构、特性、性质、元素或限度,但并非每个实施例或示例都必然包括特定特征、结构、特性、性质、元素或限度。另外,重复使用短语“根据本申请的一个实施例”虽然有可能是指代相同实施例,但并非必然指代相同的实施例。
为简单起见,以下描述中省略了本领域技术人员公知的某些技术特征。
根据本申请的一个实施例,提供一种基于RC逆变换的核脉冲信号处理方法,如图1所示,包括以下步骤:
S1、将前置放大器电路中对负指数核脉冲信号有影响的分布电容和分布电阻等效为RC电路,RC电路的电路图参考图2;
S2、通过RC电路微分数值分析方法得到负指数核脉冲信号与双指数核脉冲信号的函数表达式,RC电路Vi与Vo的关系可依据基尔霍夫电流定律用下式表示:
经过整理得到
其中,Vi为负指数核脉冲信号,Vo为双指数核脉冲信号,R为核探测仪器中ADC之前的电路阻值,C为核探测仪器中ADC之前的容值,t为时间。
S3、通过ADC采样获取前置放大器输出的双指数核脉冲信号,并根据所述双指数核脉冲信号和所述函数表达式计算得到负指数核脉冲信号;
本步骤具体为:通过ADC采样,将前置放大器输出的双指数核脉冲信号转化为具有一定时间间隔的离散信号,并将所述函数表达式根据微分定义转化为微分函数表达式,即令V0=y[n],Vi=x[n]根据微分的定义dV0=y[n]-y[n-1],dt=Δt则该微分函数表达式为:
再令可得出表达式:
x[n]-y[n]=k(y[n]-y[n-1])
即,
x[n]=(k+1)y[n]-ky[n-1];
其中,如图3所示,y[n]为双指数核脉冲信号经ADC采样获取的具有一定时间间隔的离散信号,x[n]为负指数核脉冲信号对应的具有一定时间间隔的离散信号,R为核探测仪器中ADC之前的电路阻值,C为核探测仪器中ADC之前的容值,t为时间。
然后根据采集到的离散信号y[n]和所述微分函数表达式计算得到负指数核脉冲信号x[n],该负指数核脉冲信号为恢复上升沿的单指数信号,如图4所示。
S4、取恢复上升沿的负指数核脉冲信号的上升沿和双指数核脉冲信号的下降沿构成新的核脉冲信号,如图5所示,恢复后的信号上升沿变陡,从双指数信号变成了单指数信号。
图6中峰值较高的曲线为没有经过逆运算的核脉冲信号的滤波成形效果,其平顶是倾斜的;峰值较低的曲线为采用本方法得到的新的核脉冲信号的滤波成形效果,改善了平顶倾斜的问题。
从图7、图8中可以看出,恢复上升沿的负指数核脉冲信号滤波后形成的能谱图的谱峰要比有经过逆运算的核脉冲信号形成的谱峰靠后一些,说明通过RC逆运算,扩大了系统的电压测量范围。
以上所述实施例仅表示本发明的几种实施方式,其描述较为具体和详细,但并不能理解为对本发明范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明保护范围。因此本发明的保护范围应该以所述权利要求为准。
Claims (1)
1.基于RC逆变换的核脉冲信号处理方法,其特征在于,包括以下步骤:
S1、将影响前置放大器电路中负指数核脉冲信号的分布电容和分布电阻等效为RC电路;
S2、通过RC电路微分数值分析方法得到负指数核脉冲信号与双指数核脉冲信号的函数表达式,所述函数表达式为:
其中Vi为负指数核脉冲信号,Vo为双指数核脉冲信号,R为核探测仪器中ADC之前的电路阻值,C为核探测仪器中ADC之前的容值,t为时间;
S3、通过ADC采样获取前置放大器输出的双指数核脉冲信号,并根据所述双指数核脉冲信号和所述函数表达式计算得到负指数核脉冲信号,计算得到所述负指数核脉冲信号为单指数信号,通过ADC采样,将前置放大器输出的双指数核脉冲信号转化为具有一定时间间隔的离散信号,并将所述函数表达式根据微分定义转化为微分函数表达式,根据所述离散信号和所述微分函数表达式计算得到负指数核脉冲信号;所述微分函数表达式为:
其中,y[n]为双指数核脉冲信号经ADC采样获取的具有一定时间间隔的离散信号,x[n]为负指数核脉冲信号对应的具有一定时间间隔的离散信号,R为核探测仪器中ADC之前的电路阻值,C为核探测仪器中ADC之前的容值,t为时间;
S4、采用计算得到的所述负指数核脉冲信号的上升沿和双指数核脉冲信号的下降沿构建新的核脉冲信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610347460.0A CN106019357B (zh) | 2016-05-23 | 2016-05-23 | 基于rc逆变换的核脉冲信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610347460.0A CN106019357B (zh) | 2016-05-23 | 2016-05-23 | 基于rc逆变换的核脉冲信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106019357A CN106019357A (zh) | 2016-10-12 |
CN106019357B true CN106019357B (zh) | 2019-08-13 |
Family
ID=57092982
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610347460.0A Active CN106019357B (zh) | 2016-05-23 | 2016-05-23 | 基于rc逆变换的核脉冲信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106019357B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109471151B (zh) * | 2018-11-21 | 2020-07-21 | 成都理工大学 | 信号处理方法、装置及用户终端 |
CN109507709B (zh) * | 2018-11-21 | 2020-09-25 | 成都理工大学 | 核脉冲信号处理方法及装置 |
CN110071706B (zh) * | 2019-05-06 | 2021-03-19 | 成都理工大学 | 核脉冲信号上升沿恢复方法、装置及数字化核仪器 |
CN110276441B (zh) * | 2019-05-20 | 2021-05-04 | 成都理工大学 | 一种基于深度学习的梯形重叠核脉冲估计方法 |
CN111404542B (zh) * | 2020-03-31 | 2023-06-23 | 四川新先达测控技术有限公司 | 一种双指数核信号计数方法及装置 |
CN111600584B (zh) * | 2020-06-01 | 2023-05-26 | 四川新先达测控技术有限公司 | 一种核脉冲信号处理方法及系统 |
CN111969982A (zh) * | 2020-08-19 | 2020-11-20 | 成都理工大学 | 一种脉冲波形变换方法 |
CN112596097B (zh) * | 2020-12-11 | 2022-10-28 | 成都理工大学 | 一种基于权冲激函数的核信号前端处理系统 |
CN113466924B (zh) * | 2021-07-01 | 2023-05-05 | 成都理工大学 | 一种对称类弹头脉冲成形装置及方法 |
CN113568032B (zh) * | 2021-07-02 | 2024-02-23 | 四川新先达测控技术有限公司 | 一种基于z变换的负指数核脉冲信号处理方法及系统 |
CN114252899B (zh) * | 2022-03-02 | 2022-05-20 | 四川新先达测控技术有限公司 | 一种核信号的级联冲激卷积成形方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104242877A (zh) * | 2014-07-16 | 2014-12-24 | 成都理工大学 | 基于模拟cr-rc电路的核脉冲信号数字高斯成形方法 |
CN104483557A (zh) * | 2014-12-24 | 2015-04-01 | 中国科学技术大学 | 一种减少计数损失的脉冲幅度测量电路及方法 |
CN104635254A (zh) * | 2015-01-30 | 2015-05-20 | 陕西迪泰克新材料有限公司 | 一种便携式γ辐射谱仪 |
-
2016
- 2016-05-23 CN CN201610347460.0A patent/CN106019357B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104242877A (zh) * | 2014-07-16 | 2014-12-24 | 成都理工大学 | 基于模拟cr-rc电路的核脉冲信号数字高斯成形方法 |
CN104483557A (zh) * | 2014-12-24 | 2015-04-01 | 中国科学技术大学 | 一种减少计数损失的脉冲幅度测量电路及方法 |
CN104635254A (zh) * | 2015-01-30 | 2015-05-20 | 陕西迪泰克新材料有限公司 | 一种便携式γ辐射谱仪 |
Non-Patent Citations (3)
Title |
---|
基于冲激响应不变法的核信号数字滤波成形算法研究;葛青;《核电子学与探测技术》;20140731;第34卷(第7期);913-916 |
核电子学滤波电路数字模拟分析研究;喻杰;《核电子学与探测技术》;20150831;第35卷(第8期);817-820 |
运用VBA技术仿真信号处理电路的研究;周伟;《实验科学与技术》;20160229;第14卷(第1期);23-27 |
Also Published As
Publication number | Publication date |
---|---|
CN106019357A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106019357B (zh) | 基于rc逆变换的核脉冲信号处理方法 | |
CN107024711B (zh) | 一种闪烁脉冲数字化信号的拟合方法 | |
US20160324498A1 (en) | On-line energy coincidence method and system for all-digital pet system | |
CN104748844B (zh) | 四象限光电探测器信号处理系统 | |
CN104614756B (zh) | 一种核素能谱数字化采集系统及其采集方法 | |
CN106772545B (zh) | 一种采用尖顶成形算法的数字脉冲幅度分析器 | |
CN107817514A (zh) | 一种数字化核能谱测量系统中的脉冲阶梯成形方法 | |
CN208384130U (zh) | 一种基于fpga/stm32的数字多道谱仪 | |
CN107783173A (zh) | 一种数字化核能谱测量系统中的脉冲矩形成形方法 | |
CN103364614A (zh) | 自适应宽量程电流电压转换装置 | |
CN103236825A (zh) | 一种用于高精度数据采集系统的数据校正方法 | |
CN107526096A (zh) | 一种修正温度对SiPM增益影响的方法 | |
CN104181601A (zh) | 基于arm7的核辐射探测系统 | |
CN103713310A (zh) | 一种核脉冲信号高速数字分离方法、系统及数字化谱仪 | |
CN108964637B (zh) | 一种基于fpga的超高通过率电流型数字化脉冲处理方法 | |
CN104198826A (zh) | 一种脉冲信号检测系统、方法及粒子计数器 | |
CN111600584A (zh) | 一种核脉冲信号处理方法及系统 | |
CN105759305A (zh) | 一种自适应级联消褶积梯形合成的数字脉冲幅度分析器 | |
CN113189634B (zh) | 一种类高斯成形方法 | |
CN113568032A (zh) | 一种基于z变换的负指数核脉冲信号处理方法及系统 | |
CN104168022B (zh) | 一种基于离散时间增量型σδadc的x射线ccd读出系统 | |
CN106291652A (zh) | 一种数字类高斯脉冲成形递推算法 | |
CN103903235A (zh) | 细胞图像超分辨率重构方法 | |
CN104462679B (zh) | 辐射探测前端读出系统数字滤波成形电路设计方法 | |
US20230281266A1 (en) | Cascaded impulse convolution shaping method and apparatus for nuclear signal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |