一种基于拉格朗日插值的脉冲幅度获取方法与流程

本发明涉及核探测技术领域，尤其涉及的是一种基于拉格朗日插值的脉冲幅度获取方法。

背景技术：

脉冲幅度获取是核探测技术里最常见的测量需求，用于实现该功能的仪器称为多道脉冲分析器。

目前，脉冲幅度获取手段主要有以下两种方案：

1）模拟化方案：首先采用峰值保持电路对脉冲信号进行峰位的保持和展宽，然后利用低速率的模数转换器(adc)对峰位电平进行采样，采样完成后再对峰值保持电路进行复位，以便响应下一个脉冲；

2）数字化方案：直接采用高速率的adc直接采样脉冲信号，然后对采样数据进行数字化的处理获得幅度信息，例如寻找最大值，复杂的方案如梯形滤波等。

模拟化方案的优点是由于采用了峰值保持电路，峰位电平已被锁住，因此不需要高速率的adc；但缺点是需要采用额外的硬件电路，例如峰值保持电路。

数字化方案的优点是可直接对信号进行采样，不需要模拟方案中的峰值保持电路；但缺点是需要高速的adc，成本高、功耗大；而且，高速adc还需搭配可编程逻辑器件(fpga)或者数字信号处理器(dsp)来进行数据处理，导致电路结构变得复杂。

因此，现有技术尚有待改进和发展。

技术实现要素：

为解决上述技术问题，本发明提供一种基于拉格朗日插值的脉冲幅度获取方法，电路结构更加简单，不需要采用额外的硬件电路，且成本低、功耗小。

本发明的技术方案如下：一种基于拉格朗日插值的脉冲幅度获取方法，包括以下步骤：

a、采样脉冲信号波形峰值附近至少3个连续数据点的时刻和幅度，且这些数据点的幅度根据其时间顺序同时满足从小到大再到小的关系；

b、将所述数据点的时刻和幅度代入拉格朗日插值函数式中，得到新的时刻t对应的幅度y(t)；

c、从幅度最大的第n个数据点开始，以yn作为起始峰值ypeak，向前指定时间间隔δt1，代入所述拉格朗日插值函数式，得到(tpeak-δt1)时刻数据点的幅度yt-；同时，向后指定时间间隔δt2，代入所述拉格朗日插值函数式，得到(tpeak+δt2)时刻数据点的幅度yt+；

d、比较yt-、ypeak和yt+之间的大小：若ypeak最大，则获得脉冲信号波形的峰值并结束；若ypeak＜yt-，则进入步骤e之后结束；若yt-＜ypeak＜yt+，则进入步骤f之后结束；

e、以yt-作为新的峰值ypeak，继续向前指定时间间隔δt1，代入所述拉格朗日插值函数式，得到(tpeak-δt1)时刻数据点的幅度yt-，并比较ypeak和yt-的大小，直到当ypeak≥yt-时，获得脉冲信号波形的峰值幅度；

f、以yt+作为新的峰值ypeak，继续向后指定时间间隔δt2，代入所述拉格朗日插值函数式，得到(tpeak+δt2)时刻数据点的幅度yt+，并比较ypeak和yt+的大小，直到当ypeak≥yt+时，获得脉冲信号波形的峰值。

所述的基于拉格朗日插值的脉冲幅度获取方法，其中：在所述步骤a中，采样脉冲信号波形峰值附近5个连续数据点的时间和幅度，且这5个数据点的幅度根据其时间顺序同时满足yn-2＜yn-1…yn-1≤yn…yn≥yn+1…yn+1＞yn+2的关系式。

所述的基于拉格朗日插值的脉冲幅度获取方法，其中：所述指定时间间隔δt1和δt2相等，且均与采样率相关。

所述的基于拉格朗日插值的脉冲幅度获取方法，其中：对于2mhz采样率的adc，重建采样率20mhz的数据序列，对脉冲信号波形峰值段进行拟合，指定时间间隔δt1和δt2均为50ns。

所述的基于拉格朗日插值的脉冲幅度获取方法，其中：所述脉冲信号波形为高斯滤波成形信号。

所述的基于拉格朗日插值的脉冲幅度获取方法，其中：所述高斯滤波成形信号采用cr-rc²成形电路输出，成形时间1.5μs。

本发明所提供的一种基于拉格朗日插值的脉冲幅度获取方法，由于采用了拉格朗日插值算法重建脉冲信号幅度，只需采用低速率的模数转换器，使得电路的结构更加简单，既不需要采用峰值保持电路等额外的硬件电路，也不需要搭配可编程逻辑器件或者数字信号处理器来进行数据处理，成本低、功耗小。

附图说明

图1是本发明基于拉格朗日插值的脉冲幅度获取方法实施例的流程图；

图2是本发明基于拉格朗日插值的脉冲幅度获取方法实施例的重建波形图；

图3是本发明基于拉格朗日插值的脉冲幅度获取方法实施例所用测试波形图；

图4是本发明没有经过任何插值处理的脉冲幅度谱图；

图5是本发明经过拉格朗日插值处理的脉冲幅度谱图。

具体实施方式

以下将结合附图，对本发明的具体实施方式和实施例加以详细说明，所描述的具体实施例仅用以解释本发明，并非用于限定本发明的具体实施方式。

如图1所示，图1是本发明基于拉格朗日插值的脉冲幅度获取方法实施例的流程图，本发明基于拉格朗日插值的脉冲幅度获取方法包括以下步骤：

步骤s110、找出脉冲信号波形峰值附近多个连续的数据点用来进行插值；以5个数据点为例，如果采样数据中有5个连续数据点的幅度根据其时间顺序同时满足yn-2＜yn-1…yn-1≤yn…yn≥yn+1…yn+1＞yn+2的关系式，即这些数据点的幅度根据其时间顺序同时满足从小到大再到小的关系，就可以采用这5个数据点的时刻和幅度数据进行插值，其中，yn代表采样数据中第n个数据点的幅度；

理论上讲，用于进行插值的数据点越多，插值重建的精度就越高，但是所需的运算时间也就越长，运算的时间与k²(k是用来插值的数据点个数)成正比关系，而且，随着数据点个数的增加，精度的提高程度也越来越不明显；因此，对于核探测里最常用的高斯滤波成形信号，选用5个数据点已能够满足插值重建的精度要求，而选用3个数据点或者4个数据点进行插值重建的结果则要稍差一点；

另外，在上面的关系式中，由于脉冲信号的波形都是从最大值(即峰位)处向两边递减的，因此不管采样时刻如何，总能在脉冲波形峰值左右找到5个满足上述关系式的数据点。

步骤s120、将这5个数据点的时刻和幅度数据代入拉格朗日插值函数式中，得到新的时刻t对应的幅度y(t)；其中，interpolate()为拉格朗日插值函数，t为时间变量，k是用来插值的数据点个数，(tj，yj)表示第j个数据点的时刻和幅度，tm表示第m个数据点的时刻，0≤m≤k，m≠j；

从幅度最大的第n个数据点开始，以该数据点的幅度yn作为起始峰值ypeak，向前指定时间间隔δt1例如50ns，代入所述拉格朗日插值函数式，得到(tpeak-δt1)时刻数据点的幅度yt-；

同时，向后指定时间间隔δt2例如50ns，代入所述拉格朗日插值函数式，得到(tpeak+δt2)时刻数据点的幅度yt+。

步骤s130、比较yt-、ypeak和yt+之间的大小，这里分成三种情况：

1)若ypeak最大，即ypeak≥yt+并且ypeak≥yt-，则获得脉冲信号波形的峰值并结束(即步骤s160)；

2)若ypeak＜yt-，说明脉冲波形峰值在左侧，则进入步骤s140；

3)若yt-＜ypeak＜yt+，说明脉冲波形峰值在右侧，则进入步骤s150。

步骤s140、以yt-作为新的峰值ypeak，继续向前指定时间间隔δt1例如50ns，代入所述拉格朗日插值函数式，得到(tpeak-δt1)时刻数据点的幅度yt-；比较ypeak和yt-的大小，直到当ypeak≥yt-时(即步骤s145)，获得脉冲信号波形的峰值并结束(即步骤s160)。

步骤s150、以yt+作为新的峰值ypeak，继续向后指定时间间隔δt2例如50ns，代入所述拉格朗日插值函数式，得到(tpeak+δt2)时刻数据点的幅度yt+，并比较ypeak和yt+的大小，直到当ypeak≥yt+时(即步骤s155)，获得脉冲信号波形的峰值并结束(即步骤s160)。

在本发明基于拉格朗日插值的脉冲幅度获取方法的优选实施方式中，以2mhz采样率的adc为例，目前市面上的许多单片机的内置adc都可达到该采样水平，采样时间间隔500ns，以核探测里最常用的高斯滤波成形信号作为测试信号，得到由2mhz采样信号点p1(即空心点)组成的原始信号波形s1(即实心线)图如图2所示，图2是本发明基于拉格朗日插值的脉冲幅度获取方法实施例的重建波形图，横坐标为时间time，单位是μs，纵坐标为电压信号，单位mv，由图2可知，该测试信号波形的达峰时间约2μs，波形底部宽度约10μs，可以满足几十khz计数率以下的测量需求；本实施例中将拉格朗日插值算法的数据间隔δt1和δt2均设定为50ns，即重建一个20mhz采样率的数据序列，并最终得到由拉格朗日插值的数据点(即实心点)组成的峰值位置处的波形s2(即点化线)，显然所述指定时间间隔δt1和δt2均与采样率相关；不难看出，尽管采用了低速率的adc，使得原始采样数据难以保证准确采样到峰值电平，但是经过拉格朗日插值算法之后，在不增加峰位保持电路等额外硬件电路的情况下，仍然可以准确地重建脉冲信号幅度，且成本低、功耗小。

最后，以核电子学里常用的cr-rc2成形电路为例，成形时间1.5μs，其输出波形如图3所示，图3是本发明基于拉格朗日插值的脉冲幅度获取方法实施例所用测试波形图，横坐标为时间，单位是μs，纵坐标为幅度，无单位；本实施例测试了10000个脉冲信号，对每个脉冲信号从随机位置开始采样，adc采样频率2mhz，如果没有后续的拉格朗日插值处理，直接将采样最大值作为幅度，获得的幅度谱图如图4所示，图4是本发明没有经过任何插值处理的脉冲幅度谱图，横坐标为幅度，纵坐标为计数；而经过5点拉格朗日插值处理之后获得的幅度谱图如图5所示，图5是本发明经过拉格朗日插值处理的脉冲幅度谱图，横坐标为幅度，纵坐标为计数；可以看出，没有后续的任何插值处理或者不增加峰位保持电路，确实难以准确重建脉冲信号幅度，而经过后续的拉格朗日插值处理，即使不增加峰位保持电路，仍然可以准确地重建脉冲信号幅度，非常适合用于计数率要求不高的测试场合：采用低速率的adc对信号直接进行采样，然后利用拉格朗日插值算法获取幅度信息。低速率的adc采样，并不能保证可以准确采集到脉冲的峰值电平，但是经过拉格朗日插值算法之后，却可以很好地还原信号的幅度。

应当理解的是，以上所述仅为本发明的较佳实施例而已，并不足以限制本发明的技术方案，对本领域普通技术人员来说，在本发明的精神和原则之内，可以根据上述说明加以增减、替换、变换或改进，而所有这些增减、替换、变换或改进后的技术方案，都应属于本发明所附权利要求的保护范围。