一种幅度定比弹道亏损校正方法与流程

1.本发明属于核脉冲信号处理技术领域，具体涉及一种幅度定比弹道亏损校正方法。


背景技术：

2.电荷灵敏前置放大器在高分辨测量领域应用十分广泛，图1为阻容反馈式电荷灵敏前置放大器的原理图，其中rf的存在决定了这一类电荷灵敏前置放大器的一个固有属性，即信号幅度弹道亏损，图2中(u
m-u1)即为幅度亏损的部分。
3.总体上，弹道亏损分为两类：一致性亏损和非一致性亏损。一致性亏损是指：如果所有信号上升时间一致，则所有信号的亏损比例也是一致的，此时对于能谱来说就不需要校正了，如闪烁体探测器的信号，其上升时间主要是光传输到光阴极的时间与电子倍增时间之和，这两者变化很小，所以信号幅度不需要校正。非一致性亏损是指：对于半导体探测器来说，粒子在晶体内作用位置不一样，电荷运动到电极被收集的时间就会有明显差异，具体与探测器的尺寸和载流子运动速度有关，这样就会带来不一致的弹道亏损。因此，该类探测器如果要获得好的能量分辨率，必须充分考虑弹道亏损的不一致性，只有在校正了弹道亏损后才能实现较高的能量分辨率，如高纯锗和czt探测器等。
4.由于亏损大小与电荷被收集的时间相关，因此早期的弹道亏损校正技术是由f.goulding等(f.goulding等,large coaxial germanium detectors-correction for ballistic deficit and trapping losses,ieee transactions on nuclear science ns-37，1990)提出的在硬件电路上对信号上升时间进行补偿。随着数据信号处理技术的兴起，发展了包括梯形滤波在内的多种数字滤波成形技术(v.radeka,trapezoidal filtering of signals from large germanium detectors at high rates,ieee transactions on nuclear science ns-19，1972；v.t.jordanov等,digital synthesis of pulse shapes in real time for high resolution radiation spectroscopy,nuclear instruments and methods a 345，1994；a.kalinin等，pulse shaping for ge-spectrometers optimized for ballistic deficit and electronic noise,nuclear instruments and methods a 538，2005；shefali saxena等，digital pulse deconvolution with adaptive shaping for real-time high-resolution high-throughput gamma spectroscopy，nuclear inst.and methods in physics research,a，2018)。2016年，成都理工大学chuan-yun xiong等人(chuan-yun xiong等，ballistic deficit compensation method for a large-volume hpge detector at high count rates，nucl sci tech(2016)27:67)提出以数字信号处理的方式将信号上升时间补偿法应用到高纯锗探测器的弹道亏损校正上。
5.上述技术均采用绝对补偿的方式进行弹道亏损校正，即亏损多少就校正多少。亏损校正时采用逐点微分式校正方式，存在单次校正量小、计算量大以及对低能射线校正能力有限等问题，对于直接从前置放大器获取数字化信号进行处理的模式尤为显著，并不能
适用于高、低能射线的弹道亏损校正。


技术实现要素：

6.本发明的目的是解决现有弹道亏损校正技术均采用绝对补偿的方式，亏损校正时采用逐点微分式校正存在单次校正量小、计算量大以及对低能射线校正能力有限等问题，而提供了一种幅度定比弹道亏损校正方法。
7.本发明的思路是：先通过单延迟线成形去除核脉冲的长拖尾，然后对成形后的脉冲进行积分获得与理论最大幅度成定比例关系的归一化积分面积，从而可准确表征理论脉冲幅度，达到弹道亏损校正目的。该幅度定比弹道亏损校正方法显著减小了计算量，对于低能射线的能量校正适用性更强。
8.通过数学变换，使所有的信号幅度均可归结为kum的形式，k为恒定的比例系数，这样，所有的信号均与自身的理论幅度值形成了固定比例关系，这与一致性亏损原理是相同的。实际上，绝对校正技术是幅度定比校正技术中k＝1的特例。
9.为达到上述目的，本发明采用的技术方案为：
10.一种幅度定比弹道亏损校正方法，其特殊之处在于，包括以下步骤：
11.步骤1、利用单延迟线成形技术，以特定的脉冲延迟时间，对单个核脉冲进行成形，得到成形后的脉冲波形；
12.步骤2、计算成形后的脉冲波形总的积分面积，得到单个核脉冲信号幅值的比例系数k；
13.步骤3、选定增益系数z，z＞0，进行信号幅度归一化处理，求解得到校正系数k0，完成核脉冲的幅度定比弹道亏损校正。
14.进一步地，步骤1中，所述单延迟线成形计算为：
[0015][0016]
其中，u
t
为成形后的脉冲波形t时刻的数据；v
t
为原始脉冲波形t时刻数据；
[0017]
td为脉冲延迟时间，td≥tr；tr为信号上升时间；
[0018]
为原始脉冲波形(t-td)时刻的数据；
[0019]
λ为原始脉冲波形的衰减常数。
[0020]
进一步地，步骤2具体为：
[0021]
单延迟线成形后的脉冲波形，包括信号上升段，信号衰减段和信号下降段；
[0022]
2.1、信号上升段1的积分面积s1为：
[0023][0024]
其中，um为理论最大信号幅度；
[0025]
2.2、信号衰减段2的积分面积s2为：
[0026][0027]
2.3、信号下降段3的积分面积s3为：
[0028][0029]
2.4、成形后的脉冲波形总的积分面积s
sum
为：
[0030][0031]
2.5、单个核脉冲信号幅值的比例系数k：
[0032]
令求解得到：
[0033][0034]
进一步地，步骤3中，校正后的信号幅度u
out
：
[0035][0036]
求解得到校正系数k0：
[0037][0038]
与现有技术相比，本发明具有的有益技术效果如下：
[0039]
本发明提出的幅度定比弹道亏损校正方法，通过单延迟线成形去除核脉冲的长拖尾，然后对成形后的脉冲进行积分获得与理论最大幅度成定比例关系的归一化积分面积，从而可准确表征理论脉冲幅度，达到弹道亏损校正目的，使所有的信号幅度均可归结为与理论最大信号幅度之间成恒定的比例系数关系，且比例系数与脉冲上升时间无关；利用本方法在弹道亏损校正时计算量显著减小，并且本方法对高、低能射线的弹道亏损的准确校正均适用。
附图说明
[0040]
图1为电荷灵敏前置放大器的基本结构示意图；
[0041]
图2为恒定电流收集模式理论脉冲及弹道亏损后的脉冲示意图；
[0042]
图3为本发明幅度定比弹道亏损校正方法流程图；
[0043]
图4为本发明实施例单延迟线成形示意图。
具体实施方式
[0044]
为使本发明的目的、优点和特征更加清楚，以下结合附图和具体实施例对本发明提出的一种幅度定比弹道亏损校正方法作进一步详细说明。本领域技术人员应当理解的是，这些实施方式仅仅用来解释本发明的技术原理，目的并不是用来限制本发明的保护范围。
[0045]
以电荷灵敏前置放大器的输出核脉冲为研究对象，具体阐述采用本发明方法实现
核脉冲幅度值的弹道亏损准确校正的过程。
[0046]
如图1和图2所示，电荷灵敏前置放大器的基本结构示意图及其输出的脉冲信号图。
[0047]
理论上电荷灵敏前置放大器输出的原始脉冲信号为恒定电流收集模式，如图2中曲线5，理论最大信号幅度为um，电荷灵敏前置放大器的时间常数为τ＝r
fcf
，其中rf和cf分别为反馈电阻和反馈电容，则原始脉冲波形的衰减常数为λ＝1/τ，信号上升时间为tr。实际上，输出的原始脉冲信号存在弹道亏损，实际脉冲最大信号幅度为u1，而且该弹道亏损随着脉冲信号上升时间tr而改变，如图2中曲线4。
[0048]
以上为初始条件，现通过幅度定比弹道亏损校正方法，将各种亏损不一致的实际脉冲幅度u1变换成统一的幅度u
out
，且u
out
＝kum，其中比例系数k与信号上升时间tr无关，也就是说，无论上升时间tr如何变化，校正后的信号幅度u
out
始终与理论最大信号幅度um成定比例系数k。
[0049]
如图3所示，本实施例应用幅度定比弹道亏损校正方法对电荷灵敏前置放大器的输出核脉冲信号处理时，包括以下步骤：
[0050]
步骤1、利用单延迟线成形技术，以特定的脉冲延迟时间td，对电荷灵敏前置放大器的输出的单个核脉冲进行成形，得到成形后的脉冲波形；
[0051]
单延迟线成形计算为：
[0052][0053]
其中，u
t
为成形后的脉冲波形t时刻的数据；v
t
为原始脉冲波形t时刻数据；
[0054]
td为脉冲延迟时间，td≥tr；tr为信号上升时间；
[0055]
为原始脉冲波形(t-td)时刻的数据；
[0056]
λ为原始脉冲波形的衰减常数。步骤2、单延迟线成形后的脉冲波形，如图4中的实线波形所示，包括信号上升段1，信号衰减段2和信号下降段3，计算成形后的脉冲波形总的积分面积，得到单个核脉冲信号幅值的比例系数k；
[0057]
2.1、信号上升段1的积分面积s1为：
[0058][0059]
2.2、信号衰减段2的积分面积s2为：
[0060][0061]
2.3、信号下降段3的积分面积s3为：
[0062][0063]
2.4、成形后的脉冲波形总的积分面积s
sum
为：
[0064]
[0065]
2.5、单个核脉冲信号幅值的比例系数k：
[0066]
令求解得到：
[0067][0068]
λ对一个信号电路来说是定值，td为选定的固定值，可见，k为定值，所以，对于单个核脉冲信号，无论其最大幅值u1亏损有多大，其成形后的积分面积与其理论最大的幅值um始终成定比例k的关系，且与脉冲上升时间tr的无关。
[0069]
步骤3、选定合适的增益系数z，z＞0，本实施例中选定z＝td，进行信号幅度进行归一化处理，得到校正后的信号幅度u
out
：
[0070][0071]
求解得到校正系数k0：
[0072][0073]
可以看出，k0也是一个定值，与脉冲上升时间tr无关。
[0074]
于是校正后的信号幅度u
out
与理论最大信号幅度um之间就是一个定比例关系k0，不受信号核脉冲上升时间tr的变化而变化，即无论原始信号弹道亏损有多大，校正后的幅度相对于自己理论最大脉冲幅度的比例均是固定的，从而达到了校正弹道亏损的目的。
[0075]
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明技术方案的范围。