脉冲波形的拟合方法技术

本发明专利技术提供了脉冲波形的拟合方法，所述脉冲波形的拟合方法包括以下步骤：(A1)分别获得脉冲波形PObs及模板波形PModel，其中，模板波形的采样率SDaq与脉冲波形的采样率SModel之比大于10；(A2)用隐函数y＝f(x)表示所述模板波形PModel，通过函数的平移及缩放操作，得到目标函数：y＝A·f[B(x-C)]，A、B及C为待定因子；(A3)将所述脉冲波形PObs的采样点的数值代入所述目标函数，并利用非线性最小二乘法获得所述待定因子A、B及C。本发明专利技术在脉冲激光测距及脉冲激光三维测量方面提供了新方法，适用于任意波形形状的全波形激光雷达系统，同时适用于其它需要进行波形拟合的电子系统。

【技术实现步骤摘要】
脉冲波形的拟合方法
本专利技术涉及波形拟合，特别涉及脉冲波形的拟合方法。
技术介绍
全波形激光雷达(Waveform-DigitizingLiDAR)将发射脉冲信号和回波脉冲信号均以很小的采样间隔进行采样并记录，用户根据实际应用需求，对记录的波形数据进行处理和分析，相比传统激光雷达，可以得到更丰富的激光回波次数和目标特征信息。波形数据处理和分析的关键内容是如何进行准确稳定的波形分解。模型拟合法是最常用的波形分解方法，其假设激光雷达的脉冲波形符合某个数学模型，然后使用非线性最小二乘法计算出数学模型的具体参数。其中，最常用的数学模型是高斯函数或广义高斯函数，然而受脉冲激光器和光电探测器输出特性影响，激光脉冲信号符合类高斯函数，因此使用高斯函数或广义高斯函数作为数学模型，并不准确。
技术实现思路
为解决上述现有技术方案中的不足，本专利技术提供了一种精度高、结构简单、低成本的复杂环境中弱吸收气体的检测装置。本专利技术的目的是通过以下技术方案实现的：脉冲波形的拟合方法，所述脉冲波形的拟合方法包括以下步骤：(A1)分别获得脉冲波形PObs及模板波形PModel，其中，模板波形的采样率SDaq与脉冲波形的采样率SModel之比大于10；(A2)用隐函数y＝f(x)表示所述模板波形PModel，通过函数的平移及缩放操作，得到目标函数：y＝A·f[B(x-C)]，A、B及C为待定因子；(A3)将所述脉冲波形PObs的采样点的数值代入所述目标函数，并利用非线性最小二乘法获得所述待定因子A、B及C。根据上述的脉冲波形的拟合方法，优选地，在步骤(A1)中，所述脉冲波形PObs的获得方式为：电脉冲模拟信号送入采样率为SDaq的数据采集装置转换得到数字信号PDaq＝(pDaq(1,v1),...,pDaq(i,vi),...,pDaq(m,vm))，其中m表示采样点数，vi表示第i个采样点的电压值；pDaq(i,vi)为数字信号PDaq的峰值点，取该峰值点及其前m1个和后m2个共m1+1+m2个点组成脉冲波形m1和m2的选取规则为：floor()表示向零方向取整函数，trise和tfall分别表示电脉冲信号上升沿和下降沿持续时间；取峰值点pDaq(i,vi)前l1个和后l2个点之间的采样时间总和为脉冲波形的半高宽τHw_Daq：l1和l2的选取规则为：根据上述的脉冲波形的拟合方法，优选地，在步骤(A1)中，所述模板波形PModel的获得方式为：电脉冲信号送入采样率为SModel的数据采集装置转换得到数字信号PRaw＝(pRaw(1,u1),...,pRaw(j,uj),...,pRaw(n,un))，n为采样点数，插值倍数为nInter，uj表示第j个采样点的电压值；pRaw(j,uj)为波形的峰值点，取该峰值点及其前n1个点和后n2个点共(n1+1+n2)个点组成模板波形，n1和n2的选取规则为：取峰值点pRaw(j,uj)前k1个和后k2个之间的采样时间总和为模板波形的半高宽τHw_Model：k1和k2的选取规则为：将模板信号沿时间轴向原点方向整体移动j个采样间隔，得到模板波形根据上述的脉冲波形的拟合方法，优选地，在步骤(A3)，所述待定因子A、B及C的获得方式为：设定待定因子A、B及C的初始值分别为：根据所述脉冲波形的采样点的数值计算隐函数导数f′est和目标函数值yest；Aest、Best和Cest分别表示待定因子A、B及C的中间迭代值；待定因子偏导数和的计算方式为：拟合得到的最优待定因子为Alm、Blm和Clm，则拟合得到的脉冲波形峰值电压vpk、脉宽τhw和峰值时刻tpk分别为：与现有技术相比，本专利技术具有的有益效果为：本专利技术创造性地使用隐函数模型表示模板波形，完全消除了显函数(如高斯函数或广义高斯函数)模型的欠准确缺点，通过隐函数平移与缩放操作，实现了脉冲波形的拟合。由于模板波形可以为任意形状的波形，因此该方法在波形拟合方面提供了新方法，适用于任意波形形状拟合的电子系统。附图说明参照附图，本专利技术的公开内容将变得更易理解。本领域技术人员容易理解的是：这些附图仅仅用于举例说明本专利技术的技术方案，而并非意在对本专利技术的保护范围构成限制。图中：图1为本专利技术的脉冲波形与模板波形拟合示意图；图2为本专利技术的非线性最小二乘法的隐函数导数和目标函数值计算流程图。具体实施方式图1~2和以下说明描述了本专利技术的可选实施方式以教导本领域技术人员如何实施和再现本专利技术。为了教导本专利技术技术方案，已简化或省略了一些常规方面。本领域技术人员应该理解源自这些实施方式的变型或替换将在本专利技术的范围内。本领域技术人员应该理解下述特征能够以各种方式组合以形成本专利技术的多个变型。由此，本专利技术并不局限于下述可选实施方式，而仅由权利要求和它们的等同物限定。实施例：本专利技术实施例的脉冲波形的拟合方法，所述脉冲波形的拟合方法包括以下步骤：(A1)获得脉冲波形，具体为以下方式：(a)电脉冲模拟信号送入采样率为SDaq的数据采集装置转换得到数字信号PDaq＝(pDaq(1,v1),...,pDaq(i,vi),...,pDaq(m,vm))，其中m表示采样点数，vi表示第i个采样点的电压值；数据采集装置的采样率需符合采样定律；(b)设pDaq(i,vi)为数字信号PDaq的峰值点，取该峰值点及其前m1个和后m2个共m1+1+m2个点组成脉冲波形m1和m2的选取规则为：其中，floor()表示向零方向取整函数，trise和tfall分别表示电脉冲信号上升沿和下降沿持续时间；(c)取峰值点pDaq(i,vi)前l1个和后l2个点之间的采样时间总和为脉冲波形的半高宽τHw_Daq：其中l1和l2的选取规则为：获得模板波形，具体为以下方式：(a)使用采样率为SModel的数据采集装置采集电脉冲信号，具备并使用插值功能，插值倍数为nInter；为消除脉冲能量不稳定性和探测器电噪音对电脉冲信号的影响，数据采集装置具备并使用平均功能；为采集到更精确的模板波形，采样率SModel至少比采样率SDaq高10倍；(b)电脉冲信号经数据采集装置采集得到数字信号PRaw＝(pRaw(1,u1),...,pRaw(j,uj),...,pRaw(n,un))，其中n为采样点数，uj表示第j个采样点的电压值。设pRaw(j,uj)为波形的峰值点，取该峰值点及其前n1个点和后n2个点共(n1+1+n2)个点组成模板波形，n1和n2的选取规则为：(c)取峰值点pRaw(j,uj)前k1个和后k2个之间的采样时间总和为模板波形的半高宽τHw_Model：其中k1和k2的选取规则为：(d)将模板信号沿时间轴向原点方向整体移动j个采样间隔，得到模板波形根据函数平移与缩放原理，波形拟合具体实现过程如下：(A2)使用隐函数y＝f(x)表示模板波形PModel，通过函数平移与缩放操作，得到目标函数：y＝Af[B(x-C)](7)其中，A表示沿y轴的缩放因子，B表示沿x轴的缩放因子，C表示沿x轴的平移因子，如图1所示；(A3)使用非线性最小二乘法确定待定因子A、B和C。非线性最小二乘法的具体原理和算法不做详细介绍，这里以脉冲波形PObs中采样点pDaq(i,vi)为例，详细说明与本专利技术相关的内容：计算待定因子初始值A0、B0和C本文档来自技高网...

【技术保护点】
脉冲波形的拟合方法，所述脉冲波形的拟合方法包括以下步骤：(A1)分别获得脉冲波形PObs及模板波形PModel，其中，模板波形的采样率SDaq与脉冲波形的采样率SModel之比大于10；(A2)用隐函数y＝f(x)表示所述模板波形PModel，通过函数的平移及缩放操作，得到目标函数：y＝A·f[B(x‑C)]，A、B及C为待定因子；(A3)将所述脉冲波形PObs的采样点的数值代入所述目标函数，并利用非线性最小二乘法获得所述待定因子A、B及C。

【技术特征摘要】
1.脉冲波形的拟合方法，所述脉冲波形的拟合方法包括以下步骤：(A1)分别获得脉冲波形PObs及模板波形PModel，其中，模板波形的采样率SDaq与脉冲波形的采样率SModel之比大于10；所述脉冲波形PObs的获得方式为：电脉冲模拟信号送入采样率为SDaq的数据采集装置转换得到数字信号PDaq＝(pDaq(1,v1),...,pDaq(i,vi),...,pDaq(m,vm))，其中m表示采样点数，vi表示第i个采样点的电压值；pDaq(i,vi)为数字信号PDaq的峰值点，取该峰值点及其前m1个和后m2个共m1+1+m2个点组成脉冲波形m1和m2的选取规则为：floor()表示向零方向取整函数，trise和tfall分别表示电脉冲信号上升沿和下降沿持续时间；取峰值点pDaq(i,vi)前l1个和后l2个点之间的采样时间总和为脉冲波形的半高宽τHw_Daq：l1和l2的选取规则为：(A2)用隐函数y＝f(x)表示所述模板波形PModel，通过函数的平移及缩放操作，得到目标函数：y＝A·f[B(x-C)]，A、B及C为待定因子；(A3)将所述脉冲波形PObs的采样点的数值代入所述目标函数，并利用非线性最小二乘法获得所述待定因子A、B及C；所述模板波形PModel的获得方式为：电脉冲信号送入采样率为SModel的数据采集装置转换得到数字信号PRaw＝(pRaw(1,u1),...,pRaw(j,uj),...,pRaw(n,un))，n为采样点数，插值倍数为nInter，uj表示第j个...

【专利技术属性】
技术研发人员：陈瑞强，闻路红，
申请(专利权)人：宁波华仪宁创智能科技有限公司，
类型：发明
国别省市：浙江;33