CN112639515B - 一种信号处理方法及装置 - Google Patents

一种信号处理方法及装置 Download PDF

Info

Publication number
CN112639515B
CN112639515B CN202080004821.6A CN202080004821A CN112639515B CN 112639515 B CN112639515 B CN 112639515B CN 202080004821 A CN202080004821 A CN 202080004821A CN 112639515 B CN112639515 B CN 112639515B
Authority
CN
China
Prior art keywords
probability distribution
waveform
determining
peak point
probability
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
Application number
CN202080004821.6A
Other languages
English (en)
Other versions
CN112639515A (zh
Inventor
黄龙
石现领
秦博雅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huawei Technologies Co Ltd
Original Assignee
Huawei Technologies Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Huawei Technologies Co Ltd filed Critical Huawei Technologies Co Ltd
Publication of CN112639515A publication Critical patent/CN112639515A/zh
Application granted granted Critical
Publication of CN112639515B publication Critical patent/CN112639515B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/483—Details of pulse systems
    • G01S7/486—Receivers
    • G01S7/4865—Time delay measurement, e.g. time-of-flight measurement, time of arrival measurement or determining the exact position of a peak
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/483—Details of pulse systems
    • G01S7/486—Receivers
    • G01S7/487—Extracting wanted echo signals, e.g. pulse detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/06—Systems determining position data of a target
    • G01S17/08—Systems determining position data of a target for measuring distance only
    • G01S17/10—Systems determining position data of a target for measuring distance only using transmission of interrupted, pulse-modulated waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/483—Details of pulse systems
    • G01S7/486—Receivers
    • G01S7/4865—Time delay measurement, e.g. time-of-flight measurement, time of arrival measurement or determining the exact position of a peak
    • G01S7/4866—Time delay measurement, e.g. time-of-flight measurement, time of arrival measurement or determining the exact position of a peak by fitting a model or function to the received signal

Abstract

本申请实施例提供一种信号处理方法及装置，涉及激光探测技术领域，可应用于自动驾驶或辅助驾驶。本申请实施例能够提高接收脉冲信号的峰值点信号强度和峰值点接收时间的准确性。该方法包括：激光雷达确定M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布；确定预设波形集合中的各波形为接收脉冲信号的第五概率分布；根据第一概率分布和第五概率分布，确定接收脉冲信号为预设波形集合中的各波形的第六概率分布；根据第六概率分布，从预设波形集合中确定第一预设波形；根据第一预设波形确定接收脉冲信号的峰值点接收时间和峰值点信号强度。

Description

一种信号处理方法及装置
技术领域
本申请实施例涉及激光探测技术领域，尤其涉及一种信号处理方法及装置。
背景技术
目前，激光雷达可以通过向目标对象发射激光脉冲信号，并记录激光雷达接收目标对象反射生成的接收脉冲信号的峰值点接收时间、峰值点信号强度等信息，以根据这些信息对该目标对象进行探测。例如，可以根据激光雷达发射激光脉冲信号的发射时间T1和接收脉冲信号的峰值点接收时间T2，计算激光的飞行时间(time of flight，ToF)为τ＝T2-T1。进而根据飞行时间τ以及光速，可以计算出目标对象到激光雷达的距离。再例如，还可以根据激光雷达所发出的激光脉冲信号的信号强度和接收脉冲信号的峰值点信号强度的关系，判断目标对象的表面反射率，进而判断目标对象的材质等信息。
为了确定接收脉冲信号的峰值点接收时间和峰值点信号强度，通常需要通过对接收脉冲信号进行采样，得到各采样点对应的接收时间和信号强度。然后根据多个采样点，确定接收脉冲信号的波形，进而确定接收脉冲信号的峰值点接收时间和峰值点信号强度。
例如，图1中的(a)所示为一种拟合出的接收脉冲信号的波形示意图。通过对接收脉冲信号进行采样，得到各采样点对应的接收时间和信号强度。其中，图1中每个圆圈表示一次采样。然后，根据各采样点对应的接收时间和信号强度，便可以拟合出接收脉冲信号的波形，如，图1中的(a)所示。根据拟合出的接收脉冲信号，可以确定出接收脉冲信号的峰值点接收时间和峰值点信号强度，其中峰值点接收时间为t，峰值点信号强度为P。
但在实际应用中，激光雷达在采样接收脉冲信号的过程中，可能会由于接收脉冲信号的信号强度超出激光雷达的动态范围(dynamic range)，或激光雷达中电子元件的非线性记忆效应等原因导致采集到的采样点的采样参数失真，影响探测结果的准确性。
例如，当接收脉冲信号的信号强度超出激光雷达的动态范围时，根据所采集到的采样点拟合出的接收脉冲信号的波形就会出现削顶现象，如图1中的(b)所示。结合图1中的(b)，假设在时刻t1和时刻t2之间的区域内接收脉冲信号的信号强度超出激光雷达的动态范围，则在这一区域内的采样点对应的信号强度为激光雷达动态范围的最大信号强度，而非接收脉冲信号的真实信号强度。这样，根据这些采样点拟合出的接收脉冲信号的波形就会出现削顶现象。
因此，如何才能避免由采样参数失真导致探测结果不准确的这种情况，这是保证激光雷达正常工作的关键之一。
发明内容
本申请实施例提供一种信号处理方法及装置，用于解决由于接收脉冲信号的采样点的采样参数失真而导致无法准确获取接收脉冲信号的峰值点接收时间和峰值点信号强度，进而导致探测结果不准确的问题。
第一方面，提供一种信号处理方法，包括：对接收脉冲信号进行N次采样，得到N个采样点。其中，N个采样点中的每个采样点的采样参数包括采样点对应的接收时间和信号强度，N个采样点包括M个上升沿采样点以及P个削顶采样点，N、M和P分别为大于或等于1的整数；确定M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布；其中，预设波形集合中包括至少一个预设波形；根据P个削顶采样点确定接收脉冲信号的峰值点接收时间的第二概率分布和峰值点信号强度的第三概率分布；确定接收脉冲信号的脉冲宽度的第四概率分布；根据第二概率分布、第三概率分布以及第四概率分布确定预设波形集合中的各波形为接收脉冲信号的第五概率分布；根据第一概率分布和第五概率分布，确定接收脉冲信号为预设波形集合中的各波形的第六概率分布；根据第六概率分布，从预设波形集合中确定第一预设波形；根据第一预设波形确定接收脉冲信号的峰值点接收时间和峰值点信号强度。
本实施例所提供上述技术方案，在确定接收脉冲信号的波形的过程中，首先根据采样到的上升沿采样点确定出关于接收脉冲信号的波形的似然函数(即第一概率分布)。然后，根据接收脉冲信号的先验信息，确定出接收脉冲信号的波形的先验分布(即第五概率分布)，其中先验信息包括影响接收脉冲信号波形的各种参数的概率分布，具体包括：接收脉冲信号的峰值点接收时间的概率分布(即第二概率分布)、峰值点信号强度的概率分布(即第三概率分布)以及接收脉冲信号脉冲宽度的概率分布(即第四概率分布)。进而，可以计算出接收脉冲信号的波形的后验分布(即第六概率分布)。然后，根据第六概率分布，可以从预设波形集合中确定出可能是接收脉冲信号的波形(即第一预设波形)。从而根据第一预设波形，即可以确定出接收脉冲信号的峰值点接收时间和峰值点信号强度。通过上述方法，可以更加准确的确定出峰值点信号强度和峰值点接收时间。
在一种可能的设计中，根据P个削顶采样点确定接收脉冲信号的峰值点接收时间的第二概率分布和峰值点信号强度的第三概率分布，包括：根据P个削顶采样点所对应的接收时间中的最小值，确定峰值点接收时间的取值范围的最小值；根据P个削顶采样点所对应的接收时间中的最大值，确定峰值点接收时间的取值范围的最大值；根据峰值点接收时间的取值范围，确定第二概率分布；根据P个削顶采样点对应的信号强度，确定峰值点信号强度的取值范围的最小值；根据峰值点信号强度的取值范围，确定第三概率分布。基于上述设计，可以根据削顶采样点确定出峰值点接收时间的取值范围以及峰值点信号强度的取值范围，进而基于确定出的取值范围可以确定出更加准确的第二概率分布和第三概率分布。
在一种可能的设计中，根据峰值点接收时间的取值范围，确定第二概率分布，包括：确定峰值点接收时间在峰值点接收时间的取值范围内满足均匀分布的情况下，峰值点接收时间的概率分布，得到第二概率分布。基于上述设计，本申请可以简化计算过程，快速得到第二概率分布。
在一种可能的设计中，根据峰值点信号强度的取值范围，确定第三概率分布，包括：确定峰值点信号强度在峰值点信号强度的取值范围内满足均匀分布的情况下，峰值点信号强度的概率分布，得到第三概率分布。基于上述设计，本申请可以简化计算过程，快速得到第三概率分布。
在一种可能的设计中，根据峰值点信号强度的取值范围，确定第三概率分布，包括：确定在反射接收脉冲信号的对象表面的反射率的概率分布为均匀分布的情况下，峰值点信号强度在峰值点信号强度的取值范围内的概率分布，得到第三概率分布。基于上述设计，可以通过根据反射接收脉冲信号的对象表面的反射率的概率分布得到第三概率分布，以得到更加准确的第三概率分布。
在一种可能的设计中，接收脉冲信号的脉冲宽度的第四概率分布，包括：确定脉冲宽度在脉冲宽度的取值范围内满足均匀分布的情况下，脉冲宽度的概率分布，得到第四概率分布。基于上述设计，本申请可以简化计算过程，快速得到第四概率分布。
在一种可能的设计中，确定M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布，包括：确定M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布；根据M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布，确定第一概率分布。基于上述设计，可以利用M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布，确定出第一概率分布。
在一种可能的设计中，确定M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布，包括：根据在采样目标采样点时噪声的强度的概率分布，确定目标采样点在预设波形集合中的各波形上分布的概率分布；其中，目标采样点包括M个上升沿采样点中任一个。基于上述设计，通过利用噪声对采样结果的影响，来确定M个采样点中每个采样点在预设波形集合中的各波形上分布的概率分布，进而可以确定出更加准确的概率分布结果。
在一种可能的设计中，目标采样点在预设波形集合中的各波形上分布的概率分布，满足以下公式一：
Figure BDA0002955266620000031
其中，pk表示目标采样点k在预设波形集合中目标波形上分布的概率，yk表示目标采样点k对应的信号强度，sk表示目标波形在采集目标采样点k时的信号强度。在上述设计中，通过利用高斯分布来反映采样时的噪声强度的概率分布，进而提出利用上述公式一来计算采样点在预设波形集合中的各波形上分布的概率分布。通过上述公式一，可以快速计算出各采样点在预设波形集合中的各波形上分布的概率分布。
在一种可能的设计中，根据第一概率分布和第五概率分布，确定接收脉冲信号为预设波形集合中的各波形的第六概率分布，包括：根据第一概率分布以及第五概率分布，构建以预设波形集合为状态空间的马尔可夫链，得到马尔可夫链的样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}；其中，样本序列中各样本分别用于表征一种波形；样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}中的后n项样本{zm+1,...,zm+n-1,zm+n}用于表征第六概率分布；根据第六概率分布，从预设波形集合中确定第一预设波形，包括：计算后n项样本{zm+1,...,zm+n-1,zm+n}的期望，得到第一预设波形。在上述设计中，通过构建马尔可夫链，可以方便快捷的确定出第六概率分布，进而得到第一预设波形。
在一种可能的设计中，根据第一概率分布以及第五概率分布，构建以预设波形集合为状态空间的马尔可夫链，包括：利用Metropolis-Hastings算法，构建接受分布α(z,z')满足以下公式二的马尔可夫链：
Figure BDA0002955266620000041
其中，z为马尔可夫链的样本序列中的一项样本，z'为马尔可夫链的样本序列中z后一项样本，p(y|z')为第一概率分布中z'对应的概率，p(y|z)为第一概率分布中z对应的概率，p(z')为第五概率分布中z'对应的概率，p(z)为第五概率分布中z对应的概率，q(z',z)和q(z,z')为马尔可夫链的建议分布。
第二方面，本申请提供一种信号处理装置，该信号处理装置可以实现上述第一方面或第一方面中可能的设计中的功能。这些功能可以通过硬件实现，也可以通过硬件执行相应的软件实现。硬件或软件包括一个或多个上述功能相应的模块。如：该信号处理装置可以包括：采样单元，用于对接收脉冲信号进行N次采样，得到N个采样点，N个采样点中的每个采样点的采样参数包括采样点对应的接收时间和信号强度，N个采样点包括M个上升沿采样点以及P个削顶采样点，N、M和P分别为大于或等于1的整数；第一概率确定单元，用于确定M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布；其中，预设波形集合中包括至少一个预设波形；第二概率确定单元，用于根据P个削顶采样点确定接收脉冲信号的峰值点接收时间的第二概率分布和峰值点信号强度的第三概率分布；第二概率确定单元，还用于确定接收脉冲信号的脉冲宽度的第四概率分布；第三概率确定单元，用于根据第二概率分布、第三概率分布以及第四概率分布确定预设波形集合中的各波形为接收脉冲信号的第五概率分布；第四概率确定单元，用于根据第一概率分布和第五概率分布，确定接收脉冲信号为预设波形集合中的各波形的第六概率分布；波形确定单元，用于根据第六概率分布，从预设波形集合中确定第一预设波形；参数确定单元，用于根据第一预设波形确定接收脉冲信号的峰值点接收时间和峰值点信号强度。当然，该信号处理装置还可以包括更多或更少的单元，用于实现其他的功能。
第三方面，提供一种信号处理装置，该信号处理装置包括一个或多个处理器，该一个或多个处理器和一个或多个存储器耦合；一个或多个存储器存储有计算机指令；当一个或多个处理器执行计算机指令时，使得信号处理装置执行如上述第一方面所提供的信号处理方法。
第四方面，提供一种芯片，芯片包括处理电路和接口；处理电路用于从存储介质中调用并运行存储介质中存储的计算机程序，以执行如上述第一方面所提供的信号处理方法。
第五方面，提供一种激光雷达，包括上述第二或第三方面所提供的信号处理装置，或者包括上述第四方面提供的芯片。
第六方面，提供一种计算机可读存储介质，该计算机可读存储介质中存储有指令，当该指令运行时，执行上述第一方面或者第一方面中可能的设计中的信号处理方法。
第七方面，提供一种包含指令的计算机程序产品，当其在计算机上运行时，使得计算机可以执行上述第一方面或者第一方面中可能的设计中的信号处理方法。
其中，第二方面至第七方面中任一种设计方式所带来的技术效果可以参见上述第一方面中的不同设计方式所带来的技术效果，在此不再赘述。
附图说明
图1为本申请实施例提供的一种拟合出的接收脉冲信号的波形示意图之一；
图2为本申请实施例提供的一种激光雷达的结构示意图；
图3为本申请实施例提供的一种拟合出的接收脉冲信号的波形示意图之二；
图4为本申请实施例提供的一种拟合出的接收脉冲信号的波形示意图之三；
图5为本申请实施例提供的一种信号处理方法的流程示意图之一；
图6为本申请实施例提供的一种信号处理方法的流程示意图之二；
图7为本申请实施例提供的一种对接收脉冲信号进行采样得到的采样点的示意图；
图8为利用现有技术拟合接收脉冲信号的波形示意图；
图9为利用本申请实施例所提供的信号处理方法拟合接收脉冲信号的波形示意图；
图10为本申请实施例所提供的一种信号处理装置的结构示意图之一；
图11为本申请实施例所提供的一种信号处理装置的结构示意图之二。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行描述。为了便于清楚描述本申请实施例的技术方案，在本申请的实施例中，采用了“第一”、“第二”等字样对功能和作用基本相同的相同项或相似项进行区分。本领域技术人员可以理解“第一”、“第二”等字样并不对数量和执行次序进行限定，并且“第一”、“第二”等字样也并不限定一定不同。同时，在本申请实施例中，“示例性的”或者“例如”等词用于表示作例子、例证或说明。本申请实施例中被描述为“示例性的”或者“例如”的任何实施例或设计方案不应被解释为比其它实施例或设计方案更优选或更具优势。确切而言，使用“示例性的”或者“例如”等词旨在以具体方式呈现相关概念，便于理解。
首先，对本申请实施例所涉及的相关技术进行介绍：
激光雷达是以发射激光束探测目标对象的位置、速度等特征量的雷达系统。其工作原理是向目标对象发射激光脉冲信号，然后记录激光雷达接收到从目标对象反射回来的信号(下文称为接收脉冲信号)的峰值点接收时间、峰值点信号强度。然后根据接收脉冲信号的峰值点接收时间、峰值点信号强度以及发送的激光脉冲信号的信息，便可以获取目标对象的相关信息。例如，目标对象到激光雷达的距离、目标对象表面反射率等参数。
如图2所示，为一种激光雷达的结构示意图。其中，在激光雷达10中，波形发生器(waveform generator)101用于控制发射机(transmitter)102发射预设波形的激光脉冲信号。另外，波形发生器101和发射机102还将发射的激光脉冲信号的信号参数(如发射时间、信号强度等)发送至信号处理模块108。偏转镜(fine steering)103用于控制激光脉冲信号的方向，以使激光脉冲信号透过发射端透镜104，射向待检测的目标对象。然后，激光雷达10通过接收端透镜105接收目标对象反射的接收脉冲信号，并利用偏转镜106调整接收脉冲信号的方向以使探测模块107接收到该接收脉冲信号。然后探测模块107通过对接收脉冲信号进行检测、采样得到多个采样点，其中各采样点的采样参数可以包括采样点对应的接收时间和信号强度等等。之后探测模块107将各采样点的采样参数发送至信号处理模块108。信号处理模块108根据来自探测模块107的采样点和来自波形发生器101、发射机102的信号参数，确定目标对象的相关信息。例如，信号处理模块108根据上述采样点，拟合出接收脉冲信号的波形，以确定接收脉冲信号的峰值点接收时间和峰值点信号强度。然后根据接收脉冲信号的峰值点接收时间和峰值点信号强度，以及来自波形发生器101、发射机102的信号参数，计算目标对象到激光雷达的距离等相关信息。
需要说明的是，在本实施例中，采样点对应的接收时间可以是指在该采样点采样到的接收脉冲信号的接收时间，采样点对应的信号强度可以是指在该采样点采样到的接收脉冲信号的信号强度。
上述图2所示激光雷达10仅为了便于读者理解本申请技术方案，示例性的给出了激光雷达的原理结构。容易理解的是，其中激光雷达10中各功能模块可以由硬件、软件或软硬件结合的方式来实现。例如，其中发射机102可以为激光二极管(laser diode)，探测模块107可以包括雪崩光电二极管(Avalanche Photon Diode，APD)、模拟数字转换器(analogto digital converter，ADC)以及检测器(detector)。
目前，激光雷达在采样接收脉冲信号的过程中，若存在激光雷达发射的激光脉冲信号强度较大、目标对象的表面反射率较高，或者目标对象与激光雷达之间距离较近这些情况时，就可能使接收脉冲信号的强度较高。当接收脉冲信号的信号强度超出激光雷达的动态范围时，例如在图2所示激光雷达10中当接收脉冲信号的强度超出探测单元107的正常工作范围时，探测单元107将无法采样到大于正常工作范围的信号强度，进而根据采样点拟合出的接收脉冲信号的波形会出现削顶现象。
其中，需要说明的是，本申请实施例中所称激光雷达的动态范围，可以理解为激光雷达中用于检测接收脉冲信号的功能单元(例如图2中探测单元107)的线性工作范围。当接收脉冲信号的信号强度没有超出线性工作范围时，采样点对应的信号强度可以反映当前接收脉冲信号的信号强度，即随着接收脉冲信号的信号强度的变大，所采集的采样点对应的信号强度也相应变大。当接收脉冲信号的信号强度超出线性工作范围的最大信号强度时，激光雷达在采样点采集到的信号强度就会发生饱和，即此时采集到的采样点对应的信号强度不再随着接收脉冲信号强度的增大而增大，而是会保持线性工作范围的最大信号强度。
另外，由于激光雷达中电子元件的非线性记忆效应，如图1中的(c)所示，可以看出在接收脉冲信号的下降沿部分的采样点处，采集到的信号强度要大于实际接收脉冲信号的信号强度，进而拟合出的接收脉冲信号的波形出现拖尾现象；另外在下降沿部分的采样点处，采集到的信号强度还会出现大幅度波动，进而拟合出的接收脉冲信号的波形出现振铃现象。
进一步的，当激光雷达在各采样点采集采样参数时，若同时存在削顶现象、拖尾现象以及振铃现象时，采集到的采样点的信号强度如图1中的(d)所示。可以看出，此时除上升沿的部分采样点的信号强度能够反映接收脉冲的信号强度之外，其他采样点的信号强度都发生了失真。
另外，在连续两次发射激光脉冲信号的情况下，当两次激光脉冲信号的发射时间间隔较短时，在接收脉冲信号的下降沿进行采样时，还可能受到下一个接收脉冲信号的干扰。
可以看出，当激光雷达在各采样点采集采样参数时，若出现削顶现象、拖尾现象、振铃现象以及受到其他接收脉冲信号的干扰等情况，都可能导致在接收脉冲信号的后一部分进行采样时采样点的信号强度失真。而采样点采集的信号强度失真，就会进一步导致无法准确拟合出接收脉冲信号的波形，也就无法确定接收脉冲信号的峰值点信号强度以及峰值点接收时间。最终影响激光雷达的探测结果的准确性。因此，如何才能避免由采样参数失真导致探测结果不准确的这种情况，这是保证激光雷达正常工作的关键之一。
为了解决上述技术问题，可以通过在出现削顶现象之前的n个上升沿采样点采集采样参数，例如图3中的(a)所示在a、b、c、d和e五个采样点采集采样参数。然后根据这n个采样点进行波形拟合，拟合出接收脉冲信号的波形，例如可以利用最小二乘法拟合出接收脉冲信号的波形。然后根据拟合出的接收脉冲信号，确定接收脉冲信号的峰值点信号强度和峰值点接收时间，进而得到探测结果。但该方法中，由于只有上升沿能提供有效拟合点，因此拟合出的波形误差较大。例如，根据上述a、b、c、d和e五个采样点，利用最小二乘法拟合出的接收脉冲信号的波形可以为图3中的(b)中虚线波形所示。可以看出，其中拟合出的波形与接收脉冲信号实际波形之间存在较大误差。因此该方法无法得出准确的探测结果。
另外，还可以先对接收脉冲信号进行采样，然后设置一个滑窗，用于获取滑窗内采样点所采集到的的信号强度。当滑窗内的信号强度持续大于预设阈值时，则确定滑窗内是接收脉冲信号。然后，再在找到的接收脉冲信号中找到波形形心，波形形心的接收时间即作为接收脉冲信号的峰值点接收时间。如图4所示，在滑窗(图中虚线框所示)内采集到的信号强度持续保持在最大信号强度P，进而确定滑窗内是接收脉冲信号，然后将滑窗的中点作为接收脉冲信号的峰值点。可以看出，这种方法对滑窗的宽度有要求，如图4中当滑窗的宽度较小时，根据滑窗中点确定出的接收脉冲信号的峰值点接收时间是Xc，而实际的峰值点接收时间则是Xc’，存在误差。另外，这种方式还存在只能估计峰值点接收时间，无法确定该峰值点信号强度的问题。
针对上述问题，本申请实施例提供一种信号处理方法，在确定接收脉冲信号的波形的过程中，通过引入接收脉冲信号的先验信息，其中先验信息包括：接收脉冲信号的峰值点接收时间的概率分布、峰值点信号强度的概率分布以及接收脉冲信号脉冲宽度的概率分布，然后利用这些先验信息对接收脉冲信号的波形的限制，进而确定出更加准确的接收脉冲信号的波形，从而确定出更加准确的峰值点信号强度和峰值点接收时间。
本申请实施例提供的信号处理方法可以应用于激光雷达或其他需要确定接收脉冲信号的峰值点接收时间和峰值点信号强度的设备中，以使得上述设备确定出更加准确的接收脉冲信号。以下以激光雷达为例，对该方法进行介绍，如图5所示，该方法包括以下S201-S207：
S201、激光雷达对接收脉冲信号进行N次采样，得到N个采样点。
例如图2所示激光雷达10中，在向目标对象发射激光脉冲信号后，激光脉冲信号遇到目标对象反射生成接收脉冲信号。然后接收脉冲信号通过接收端透镜105、偏转镜106至探测模块107。然后探测模107便可以在接收到接收脉冲信号后，对接收脉冲信号进行采样得到采样点。
其中，N个采样点中每个采样点的采样参数包括该采样点对应的接收时间和信号强度。例如N个采样点包括a、b和c三个采样点，其中采样点a的采样参数包括采样点a对应的接收时间和信号强度，采样点b的采样参数包括采样点b对应的接收时间和信号强度，以及采样点c的采样参数包括采样点c对应的接收时间和信号强度。
另外，N个采样点包括M个上升沿采样点以及P个削顶采样点。其中，N、M和P分别为大于或等于1的正整数。
其中，在实际实施过程中，采样点对应的信号强度，可以为在采样点处接收脉冲信号的光功率等反映光信号强度的参数，也可以为在采样点处响应于接收脉冲信号的光信号强度产生的电压值或电流值。对于信号强度所使用的实际物理量，本申请可以不做限制。
需要说明的是，本申请实施例中，上升沿采样点指在接收脉冲信号的上升沿的采样点。例如，图1中的(d)中，包括a、b、c和d四个上升沿采样点。削顶采样点指在出现削顶现象的部分的采样点。例如，图1中的(d)中，包括e、f、g和h四个削顶采样点。
S202、激光雷达确定M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布。
其中，预设波形集合中包括至少一个预设波形。
通过上文对激光雷达的介绍可知，激光雷达在接收脉冲信号的上升沿进行采样时失真的可能性较小，因此可以利用上升沿采样点的采样参数，确定出预设波形集合中的各波是接收脉冲信号的波形的概率分布，即M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布。
示例性的，考虑到接收脉冲信号通常为高斯脉冲，也就是说接收脉冲信号的当前信号强度s与当前时间x满足以下公式(1)：
Figure BDA0002955266620000081
其中，当波形参数A、t、σ取不同值时对应的波形不同。下文中为了便于描述，用θ表示A、t、σ的集合，即θ＝(A,t,σ)。
因此，可以认为满足上述公式(1)的波形，即为可能是接收脉冲信号的波形。因此，可以利用满足上述公式(1)的波形构建预设波形集合，即预设波形集合中包括满足上述公式(1)的至少一个预设波形。
当然，可以理解的是，在本申请实施例所提供技术方案在实施过程中，也可以按照其他方式确定预设波形集合中的波形。也就是说，本领域技术人员可以根据实际需要选择合适的预设波形集合。对于预设波形集合中包括哪些波形，本申请可以不做限制。
在一种实现方式中，考虑到可以根据M个事件分别发生的概率来确定M个事件同时发生的概率，例如将M个事件分别发生的概率求积即可得到M个事件同时发生的概率。因此也可以根据M个上升沿采样点中各采样点在预设波形集合中各波形上分布的概率分布，来确定M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布。因此，如图6所示，本申请S202具体可以包括以下S2021-S2022：
S2021、激光雷达确定M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布。
例如，M个上升沿采样点中共包括五个采样点，分别为：采样点1、采样点2、采样点3、采样点4、采样点5。五个采样点的采样参数分别包括该采样点对应的接收时间和信号强度。另外，假设预设波形集合中包括满足上述公式(1)的波形。然后，针对每个采样点，根据该采样点的接收时间和信号强度，便可以计算出该采样点在预设波形集合中的各波形上分布的概率。
在一种可能的设计中，考虑到在接收脉冲信号的上升沿部分对采样点处采样参数产生干扰的主要是噪声。也就是说在采样点处采集到的信号强度，实际上反映了接收脉冲信号的当前信号强度与当前噪声的强度之和。因此，在计算目标采样点在预设波形集合中的各波形上分布的概率时，上述S2021可以包括：根据在采样目标采样点时噪声的强度的概率分布，确定目标采样点在预设波形集合中的各波形上分布的概率分布。其中目标采样点可以为M个上升沿采样点中任一采样点。
例如，假设目标采样点所采集到的信号强度为y，在采样目标采样点时噪声的强度为w，预设波形集合中波形(如称为目标波形)在目标采样点的信号强度为s。那么，若接收脉冲信号为目标波形，则需要满足w＝y-s。另外，因为y已知，所以根据在采样目标采样点时噪声的强度w的概率分布，便可以计算出w＝y-s的概率，即就是目标采样点在目标波形上分布的概率。例如，假设目标采样点所采集到的信号强度为5(信号强度的物理单位可根据实际情况确定，可以是光功率、电压或电流等)，目标波形在目标采样点的信号强度为4，在采样目标采样点时噪声的强度为1的概率为20％，则可以确定目标采样点在目标波形上分布的概率为20％。
进一步的，在一种可能的设计中，考虑到噪声的强度的概率分布通常符合均值为0，方差为N0的高斯分布，其中N0为噪声功率，即噪声的强度的概率分布p(w)满足以下公式(3)：
Figure BDA0002955266620000091
因此，本申请实施例所提供技术方案中，在确定M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布时，可以根据以下公式(4)确定目标采样点在预设波形集合中的各波形上分布的概率分布，也就是说目标采样点在预设波形集合中的各波形上分布的概率分布满足以下公式(4)：
Figure BDA0002955266620000092
其中，pk表示目标采样点k在预设波形集合中目标波形上分布的概率，yk表示目标采样点k对应的信号强度，sk表示目标波形在采集目标采样点k时的信号强度。N0的取值可根据实际需要进行设定。
S2022、激光雷达根据M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布，确定第一概率分布。
例如，在确定M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布之后，可以将M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布进行求积运算，进而确定出第一概率分布。
例如，当目标采样点在预设波形集合中的各波形上分布的概率分布满足上述公式(4)，预设波形集合中包括满足上述公式(1)的至少一个预设波形的情况下，则可以确定第一概率分布满足以下公式(5)：
Figure BDA0002955266620000101
其中，p(y|θ)表示M个上升沿采样点在波形参数A、t、σ为θ的目标波形上分布的概率，yk表示目标采样点k对应的信号强度，sk(θ)表示波形参数A、t、σ为θ的目标波形在采集目标采样点k时的信号强度。
S203、激光雷达根据P个削顶采样点，确定接收脉冲信号的峰值点接收时间的第二概率分布和峰值点信号强度的第三概率分布。
其中，考虑到若存在削顶现象，则说明接收脉冲信号的峰值点接收时间在削顶部分对应的时间范围内，并且说明接收脉冲信号的峰值点信号强度比削顶部分采样到的信号强度大。例如在图1中的(b)中，可以确定该接收脉冲信号的峰值点接收时间在t1-t2之间，并且接收脉冲信号的峰值点信号强度大于信号强度P。那么利用P个削顶采样点所采集到的采样参数，便可以确定出接收脉冲信号的峰值点接收时间和峰值点信号强度的取值范围，进而确定出峰值点接收时间和峰值点信号强度的概率分布(即第二概率分布和第三概率分布)。
因此，在一种实现方式中，如图6所示，上述根据P个削顶采样点确定第二概率分布，可以包括：
S2031、激光雷达根据P个削顶采样点所对应的接收时间中的最小值和最大值，确定峰值点接收时间的取值范围。
其中，可以根据P个削顶采样点所对应的接收时间中的最小值，确定峰值点接收时间的取值范围的最小值。另外，可以根据P个削顶采样点所对应的接收时间中的最大值，确定峰值点接收时间的取值范围的最大值。
例如，图1中的(b)中，可以看出采样点1和采样点2对应的接收时间，分别为削顶采样点所对应的接收时间中的最小值和最大值。那么根据采样点1和采样点2对应的接收时间，便可以确定出峰值点接收时间的取值范围的最小值和最大值。
S2032、激光雷达根据峰值点接收时间的取值范围，确定第二概率分布。
具体的，在确定峰值点接收时间的取值范围后，便可以根据激光雷达本身的特性等先验信息，确定出在取值范围内峰值点接收时间的概率分布(即第二概率分布)。
在一种可能的设计中，为了便于计算，可以认为峰值点接收时间在取值范围内的概率分布为均匀分布。进而确定峰值点接收时间在峰值点接收时间的取值范围内满足均匀分布的情况下，峰值点接收时间的概率分布，得到第二概率分布。
例如，以接收脉冲信号的波形满足公式(1)为例，则第二概率分布可以反映为波形参数t的概率分布p(t)，表示为：
Figure BDA0002955266620000111
其中，t1表示峰值点接收时间的取值范围的最小值，t2表示峰值点接收时间的取值范围的最大值。
需要说明的是，还可以利用其他方式确定峰值点接收时间的第二概率分布，例如可以通过对激光雷达进行多次测试，将每次测试得到的接收脉冲信号的峰值点接收时间进行统计，然后根据统计结果确定峰值点接收时间的第二概率分布，对此本申请可以不做限制。
另外，在一种实现方式中，如图6所示，上述根据P个削顶采样点确定第三概率分布，可以包括：
S2033、激光雷达根据P个削顶采样点对应的信号强度，确定峰值点信号强度的取值范围。
例如，图1中的(b)中，可以看出由于削顶现在的存在，因此削顶采样点的信号强度普遍在信号强度P附近，另外由于可以确定峰值点信号强度的最小值要大于削顶采样点的信号强度，因此可以根据P个削顶采样点对应的信号强度确定峰值点信号强度的取值范围的最小值。例如，可以将P个削顶采样点对应的信号强度中最大的信号强度作为峰值点信号强度的取值范围的最小值。
然后，在确定了峰值点信号强度的取值范围的最小值后，便可以进一步确定出峰值点信号强度的取值范围。例如，可以将激光雷达发射的激光脉冲信号的峰值点信号强度作为峰值点信号强度的最大值。然后，已知峰值点信号强度的最大值和最小值，便可以确定出峰值点信号强度的取值范围。
S2034、激光雷达根据峰值点信号强度的取值范围，确定第三概率分布。
与上述S2032类似，在确定出峰值点信号强度的取值范围后，便可以根据激光雷达本身的特性等先验信息，确定出在取值范围内峰值点信号强度的第三概率分布。
在一种可能的设计中，为了便于计算，可以认为峰值点信号强度在取值范围内的概率分布为均匀分布。进而可以通过确定峰值点信号强度在峰值点信号强度的取值范围内满足均匀分布的情况下，峰值点信号强度的概率分布，得到第三概率分布。
例如，以接收脉冲信号的波形满足公式(1)为例，则第三概率分布可以反映为波形参数A的概率分布p(A)，表示为：
Figure BDA0002955266620000112
其中，A2表示峰值点信号强度的取值范围的最小值，A1表示峰值点信号强度的取值范围的最大值。
在另一种可能的设计中，考虑到接收脉冲信号的信号强度与反射接收脉冲信号的对象表面的反射率存在关联关系，例如对象表面的反射率越高对应接收脉冲信号的信号强度通常越强。另外，由于激光雷达所探测的对象多种多样，因此可以将反射接收脉冲信号的对象表面的反射率视为均匀分布。进而可以通过确定在反射接收脉冲信号的对象表面的反射率的概率分布为均匀分布的情况下，峰值点信号强度在峰值点信号强度的取值范围内的概率分布，得到第三概率分布。
另外，还可以利用其他方式确定第三概率分布，例如可以通过对激光雷达进行多次测试，将每次测试得到的接收脉冲信号的峰值点信号强度进行统计，然后根据统计结果确定峰值点信号强度的第三概率分布，对此本申请可以不做限制。
S204、激光雷达确定接收脉冲信号的脉冲宽度的第四概率分布。
例如，可以通过测试激光雷达发出的激光脉冲信号在不同使用场景下得到的接收脉冲信号的脉冲宽度的取值范围以及脉冲宽度在取值范围内的第四概率分布。
在一种可能的设计中，为了便于计算，可以认为脉冲宽度在取值范围内的概率分布为均匀分布。进而可以通过确定在反射接收脉冲信号的对象表面的反射率的概率分布为均匀分布的情况下，峰值点信号强度在峰值点信号强度的取值范围内的概率分布，得到第四概率分布。
例如，以接收脉冲信号的波形满足公式(1)为例，则第四概率分布可以反映为波形参数σ的概率分布p(σ)，表示为：
Figure BDA0002955266620000121
其中，σ1表示脉冲宽度的取值范围的最小值，σ2表示脉冲宽度的取值范围的最大值。
S205、激光雷达根据第二概率分布、第三概率分布以及第四概率分布，确定预设波形集合中各波形为接收脉冲信号的第五概率分布。
例如，预设波形集合中一种波形(如目标波形)为接收脉冲信号的概率，可以表示为同时满足以下三点的概率：
1)该目标波形的峰值点接收时间为接收脉冲信号的峰值点接收时间；
2)该目标波形的峰值点信号强度为接收脉冲信号的峰值点信号强度；
3)该目标波形的脉冲宽度为接收脉冲信号的脉冲宽度。
也就是说，预设波形集合中一种波形为接收脉冲信号的概率，可以表示为上述三点各自分别发生的概率之积，即p(θ)＝p(t)p(A)p(σ)，其中，p(θ)表示预设波形集合中波形参数A、t、σ为θ的波形为接收脉冲信号的概率。即就是说第五概率分布可以表示为第二概率分布、第三概率分布以及第四概率分布的乘积。
S206、激光雷达根据第一概率分布和第五概率分布，确定接收脉冲信号为预设波形集合中的各波形的第六概率分布。
例如，根据贝叶斯公式可知：第一概率分布可以认为是给定M个上升沿采样点(即M个在接收脉冲信号的波形上分布的采样点)时，接收脉冲信号的波形的概率分布p(y|θ)，即第一概率分布是关于接收脉冲信号的波形的似然函数L(θ|y)；另外第五概率分布可以认为是接收脉冲信号的波形的先验分布，即p(θ)。那么，在已知L(θ|y)和p(θ)的情况下，便可以利用贝叶斯公式
Figure BDA0002955266620000131
计算出接收脉冲信号为预设波形集合中的各波形的第六概率分布P(θ|y)，即接收脉冲信号的波形的后验分布。其中，P(y)可以根据实际需要设定或在计算中可以利用相关算法约去P(y)。
在一种实现方式中，可以利用马尔可夫链蒙特卡洛(Markov Chain Monte Carlo，MCMC)算法来确定第六概率分布，具体的，上述S206可以包括：
根据第一概率分布以及第五概率分布，构建以预设波形集合为状态空间的马尔可夫链，得到马尔可夫链的样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}。
其中，马尔可夫链的样本序列中各样本分别用于表征一种波形。另外，根据MCMC算法可知，该马尔可夫链的样本序列中样本取值的概率会逐渐向第六概率分布收敛，因此，样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}中的后n项样本{zm+1,...,zm+n-1,zm+n}用于表征第六概率分布。
在一种可能的设计中，可以利用MCMC算法中的Metropolis-Hastings算法，来构建马尔可夫链。其中，该Metropolis-Hastings算法中的接受分布α(z,z')满足以下公式(6)：
Figure BDA0002955266620000132
其中，z为马尔可夫链的样本序列中的一项样本，z'为马尔可夫链的样本序列中z后一项样本，马尔可夫链的样本序列中各样本分别用于表征一种波形，p(y|z')为第一概率分布中波形z'对应的概率，p(y|z)为第一概率分布中波形z对应的概率，p(z')为第五概率分布中波形z'对应的概率，p(z)为第五概率分布中波形z对应的概率，q(z',z)和q(z,z')为马尔可夫链的建议分布。
示例性的，以接收脉冲信号为高斯脉冲，接收脉冲信号的当前信号强度s与当前时间x满足公式(1)：
Figure BDA0002955266620000133
为例，对利用Metropolis-Hastings算法构建马尔可夫链的过程进行说明，该过程可以包括以下步骤1-步骤4：
步骤1、随机设置一个初始波形，即马尔可夫链的第一项样本Z1。
具体的，由于接收脉冲信号满足公式(1)，即可以用一个波形参数θ表示一种波形。因此，初始波形Z1可以表示为θ1＝(A1,t1,σ1)。
步骤2、设置Metropolis-Hastings算法的建议分布。
在实际实施过程中，可以参照现有技术中的相关内容设置建议分布。其中，不同建议分布有着不同的收敛速度。对于采用哪种建议分布，本申请可以不做限制。
步骤3、根据建议分布，抽取一个θ2＝(A2,t2,σ2)。根据接受分布α(z,z')，计算α(θ1,θ2)。
步骤4、从区间(0,1)中随机抽取一个数u。若u＜α(θ1,θ2)，则确定马尔可夫链的第二项样本Z2为θ2；否则，重复步骤3，直至确定出马尔可夫链的第二项样本Z2。
通过迭代上述步骤3-步骤4，依次确定马尔可夫链的各项样本，进而得到马尔可夫链的样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}。
S207、根据第六概率分布，从预设波形集合中确定第一预设波形。根据第一预设波形确定接收脉冲信号的峰值点接收时间和峰值点信号强度。
在确定接收脉冲信号为预设波形集合中的各波形的第六概率分布之后，则根据第六概率分布，可以从预设波形集合中确定出可能是接收脉冲信号的波形(即第一预设波形)，例如可以计算第六概率分布的期望得到第一预设波形，或者选择第六概率分布中概率最大的波形作为第一预设波形。然后根据第一预设波形，即可以确定出接收脉冲信号的峰值点接收时间和峰值点信号强度。
在一种可能的设计中，当利用MCMC算法确定第六概率分布的情况下，则第六概率分布，从预设波形集合中确定第一预设波形，可以包括：计算马尔可夫链的样本序列中后n项样本{zm+1,...,zm+n-1,zm+n}的期望，得到第一预设波形。
本实施例所提供信号处理方法，在确定接收脉冲信号的波形的过程中，首先根据采样到的上升沿采样点确定出关于接收脉冲信号的波形的似然函数(即第一概率分布)。然后，根据接收脉冲信号的先验信息，确定出接收脉冲信号的波形的先验分布(即第五概率分布)，其中先验信息包括影响接收脉冲信号波形的各种参数的概率分布，具体包括：接收脉冲信号的峰值点接收时间的概率分布(即第二概率分布)、峰值点信号强度的概率分布(即第三概率分布)以及接收脉冲信号脉冲宽度的概率分布(即第四概率分布)。进而，可以计算出接收脉冲信号的波形的后验分布(即第六概率分布)。然后，根据第六概率分布，可以从预设波形集合中确定出可能是接收脉冲信号的波形(即第一预设波形)。从而根据第一预设波形，即可以确定出接收脉冲信号的峰值点接收时间和峰值点信号强度。通过上述方法，可以更加准确的确定出峰值点信号强度和峰值点接收时间。
以下结合实际仿真结果，对本实施例所提供信号处理方法的有益效果进行介绍：如图7所示，为激光雷达对接收脉冲信号进行采样得到的采样点(图中圆圈所示)以及利用采样点拟合出的接收脉冲信号的波形(图中采样点依次连接得到的波形)。可以看出，图中拟合的接收脉冲信号在虚线圈出的部分存在削顶现象。因此，根据图7所拟合的接收脉冲信号无法准确确定出接收脉冲信号的峰值点信号强度和峰值点接收时间。
然后，分别使用两种方法来拟合接收脉冲信号的波形。第一种，根据多个上升沿采样点，利用最小二乘法拟合出接收脉冲信号的波形，并根据拟合出的波形确定接收脉冲信号的峰值点接收时间和峰值点信号强度。第二种，利用本实施例所提供的信号处理方法，确定接收脉冲信号的波形，并根据拟合出的波形确定接收脉冲信号的峰值点接收时间和峰值点信号强度。其中两种方法所使用的接收脉冲信号的上升沿部分的采样点相同，均为图7中a、b、c和d四个采样点。
在分别利用上述两种方法进行100次蒙特卡洛仿真的过程中，在利用第一种方法拟合接收脉冲信号的波形时，时常出现图8所示的拟合出的接收脉冲信号的波形。可以看出，其中拟合出的接收脉冲信号的波形与接收脉冲信号之间存在较大误差。而图9为利用本实施例所提供的信号处理方法(即第二种方法)，拟合出的接收脉冲信号的波形。可以看出，其中利用本实施例所提供的信号处理方法的拟合表现良好。
进一步的，下表1为分别利用上述两种方法进行100次蒙特卡洛仿真后，利用仿真结果计算出的ToF的统计结果，表2为分别利用上述两种方法进行100次蒙特卡洛仿真后，利用仿真结果计算出的接收脉冲信号的幅度(即峰值点信号强度)的统计结果。
表1
ToF的平均偏差 ToF的标准差
第一种方法 27.81 250.59
第二种方法 0.14 0.31
表2
幅度的平均偏差 幅度的标准差
第一种方法 2.9886e+23 2.9886e+24
第二种方法 0.0046 0.0344
从上表1和表2可以看出，利用本申请所提供方法(即表中第二种方法)的激光雷达可以大幅提高探测结果的准确性。
在另一种实施例中，如图10所示，为本申请实施例提供的一种信号处理装置的组成示意图。该信号处理装置30可以为上述方法实施例中激光雷达的一部分，当该信号处理装置30运行时，可以使激光雷达执行上述实施例中所提供的信号处理方法。例如，在具体实现时，信号处理装置30具体可以是图2所示激光雷达中的信号处理模块108。信号处理装置30还可以为激光雷达中的芯片或片上系统。当然，信号处理装置30也可以应用于其他需要确定接收脉冲信号的峰值点接收时间和峰值点信号强度的设备中，以使得上述设备确定出更加准确的接收脉冲信号。
作为一种实现方式，如图10所示，该信号处理装置30可以包括：采样单元301，用于对接收脉冲信号进行N次采样，得到N个采样点，该N个采样点中的每个采样点的采样参数包括采样点对应的接收时间和信号强度，N个采样点包括M个上升沿采样点以及P个削顶采样点，N、M和P分别为大于或等于1的整数。
第一概率确定单元302，用于确定M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布；其中，预设波形集合中包括至少一个预设波形。
第二概率确定单元303，用于根据P个削顶采样点确定接收脉冲信号的峰值点接收时间的第二概率分布和峰值点信号强度的第三概率分布。
第二概率确定单元303，还用于确定接收脉冲信号的脉冲宽度的第四概率分布。
第三概率确定单元304，用于根据第二概率分布、第三概率分布以及第四概率分布确定预设波形集合中的各波形为接收脉冲信号的第五概率分布。
第四概率确定单元305，用于根据第一概率分布和第五概率分布，确定接收脉冲信号为预设波形集合中的各波形的第六概率分布。
波形确定单元306，用于根据第六概率分布，从预设波形集合中确定第一预设波形。
参数确定单元307，用于根据第一预设波形确定接收脉冲信号的峰值点接收时间和峰值点信号强度。
在一种可能的设计中，第二概率确定单元303，具体用于：根据P个削顶采样点所对应的接收时间中的最小值，确定峰值点接收时间的取值范围的最小值；根据P个削顶采样点所对应的接收时间中的最大值，确定峰值点接收时间的取值范围的最大值。根据峰值点接收时间的取值范围，确定第二概率分布。根据P个削顶采样点对应的信号强度，确定峰值点信号强度的取值范围的最小值。根据峰值点信号强度的取值范围，确定第三概率分布。
在一种可能的设计中，根据峰值点接收时间的取值范围，确定第二概率分布，包括：确定峰值点接收时间在峰值点接收时间的取值范围内满足均匀分布的情况下，峰值点接收时间的概率分布，得到第二概率分布。
在一种可能的设计中，根据峰值点信号强度的取值范围，确定第三概率分布，包括：确定峰值点信号强度在峰值点信号强度的取值范围内满足均匀分布的情况下，峰值点信号强度的概率分布，得到第三概率分布。
在一种可能的设计中，根据峰值点信号强度的取值范围，确定第三概率分布，包括：确定在反射接收脉冲信号的对象表面的反射率的概率分布为均匀分布的情况下，峰值点信号强度在峰值点信号强度的取值范围内的概率分布，得到第三概率分布。
在一种可能的设计中，第二概率确定单元303，具体用于确定脉冲宽度在脉冲宽度的取值范围内满足均匀分布的情况下，脉冲宽度的概率分布，得到第四概率分布。
在一种可能的设计中，第一概率确定单元302，具体用于：确定M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布。根据M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布，确定第一概率分布。
在一种可能的设计中，确定M个上升沿采样点中每个采样点在预设波形集合中的各波形上分布的概率分布，包括：根据在采样目标采样点时噪声的强度的概率分布，确定目标采样点在预设波形集合中的各波形上分布的概率分布；其中，目标采样点包括M个上升沿采样点中任一个。
在一种可能的设计中，目标采样点在预设波形集合中的各波形上分布的概率分布，满足以下公式一：
Figure BDA0002955266620000161
其中，pk表示目标采样点k在预设波形集合中目标波形上分布的概率，yk表示目标采样点k对应的信号强度，sk表示目标波形在采集目标采样点k时的信号强度。
在一种可能的设计中，第四概率确定单元305，具体用于：根据第一概率分布以及第五概率分布，构建以预设波形集合为状态空间的马尔可夫链，得到马尔可夫链的样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}；其中，样本序列中各样本分别用于表征一种波形；样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}中的后n项样本{zm+1,...,zm+n-1,zm+n}用于表征第六概率分布。波形确定单元306，具体用于计算后n项样本{zm+1,...,zm+n-1,zm+n}的期望，得到第一预设波形。
在一种可能的设计中，根据第一概率分布以及第五概率分布，构建以预设波形集合为状态空间的马尔可夫链，包括：
利用Metropolis-Hastings算法，构建接受分布α(z,z')满足以下公式二的马尔可夫链：
Figure BDA0002955266620000171
其中，z为马尔可夫链的样本序列中的一项样本，z'为马尔可夫链的样本序列中z后一项样本，p(y|z')为第一概率分布中z'对应的概率，p(y|z)为第一概率分布中z对应的概率，p(z')为第五概率分布中z'对应的概率，p(z)为第五概率分布中z对应的概率，q(z',z)和q(z,z')为马尔可夫链的建议分布。
在另一种实施例中，如图11所示为另一种数据传输装置的组成示意图。数据传输装置40包括：一个或多个处理器401以及一个或多个存储器402。一个或多个处理器401与一个或多个存储器402耦合，存储器402用于存储计算机执行指令。示例性地，在一些实施例中，当处理器401执行存储器402存储的指令时，使得该数据传输装置40执行如图5所示的S201-S207以及激光雷达需要执行的其他操作。
本申请实施例还提供一种计算机可读存储介质，该计算机可读存储介质中存储有指令，当该指令运行时，执行执行本申请实施例所提供的方法。
本申请实施例还提供一种包含指令的计算机程序产品。当其在计算机上运行时，使得计算机可以执行本申请实施例所提供的方法。
另外，本申请实施例还提供一种芯片。该芯片包括处理器。当处理器执行计算机程序指令时，使得芯片可以执行本申请实施例提供的方法。该指令可以来自芯片内部的存储器，也可以来自芯片外部的存储器。可选的，该芯片还包括作为通信接口的输入输出电路。
在上述实施例中的功能或动作或操作或步骤等，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件程序实现时，可以全部或部分地以计算机程序产品的形式来实现。该计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行计算机程序指令时，全部或部分地产生按照本申请实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一个计算机可读存储介质传输，例如，所述计算机指令可以从一个网站站点、计算机、服务器或者数据中心通过有线(例如同轴电缆、光纤、数字用户线(digital subscriber line，DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包括一个或多个可以用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质(例如，软盘、硬盘、磁带)，光介质(例如，DVD)、或者半导体介质(例如固态硬盘(solid state disk，SSD))等。
尽管结合具体特征及其实施例对本申请进行了描述，显而易见的，在不脱离本申请的精神和范围的情况下，可对其进行各种修改和组合。相应地，本说明书和附图仅仅是所附权利要求所界定的本申请的示例性说明，且视为已覆盖本申请范围内的任意和所有修改、变化、组合或等同物。显然，本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样，倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内，则本申请也意图包括这些改动和变型在内。

Claims (26)

1.一种信号处理方法，其特征在于，包括：
对接收脉冲信号进行N次采样，得到N个采样点，所述N个采样点中的每个采样点的采样参数包括所述采样点对应的接收时间和信号强度，所述N个采样点包括M个上升沿采样点以及P个削顶采样点，N、M和P分别为大于或等于1的整数；
确定所述M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布；其中，所述预设波形集合中包括至少一个预设波形；
根据所述P个削顶采样点确定所述接收脉冲信号的峰值点接收时间的第二概率分布和峰值点信号强度的第三概率分布；
确定所述接收脉冲信号的脉冲宽度在脉冲宽度的取值范围内的第四概率分布；
根据所述第二概率分布、所述第三概率分布以及所述第四概率分布确定所述预设波形集合中的各波形为所述接收脉冲信号的第五概率分布；
根据所述第一概率分布和所述第五概率分布，确定所述接收脉冲信号为所述预设波形集合中的各波形的第六概率分布；
根据所述第六概率分布，从所述预设波形集合中确定第一预设波形；
根据所述第一预设波形确定所述接收脉冲信号的峰值点接收时间和峰值点信号强度。
2.根据权利要求1所述的方法，其特征在于，所述根据所述P个削顶采样点确定所述接收脉冲信号的峰值点接收时间的第二概率分布和峰值点信号强度的第三概率分布，包括：
根据所述P个削顶采样点所对应的接收时间中的最小值，确定所述峰值点接收时间的取值范围的最小值；根据所述P个削顶采样点所对应的接收时间中的最大值，确定所述峰值点接收时间的取值范围的最大值；
根据所述峰值点接收时间的取值范围，确定所述第二概率分布；
根据所述P个削顶采样点对应的信号强度，确定所述峰值点信号强度的取值范围的最小值；
根据所述峰值点信号强度的取值范围，确定所述第三概率分布。
3.根据权利要求2所述的方法，其特征在于，所述根据所述峰值点接收时间的取值范围，确定所述第二概率分布，包括：
确定所述峰值点接收时间在所述峰值点接收时间的取值范围内满足均匀分布的情况下，所述峰值点接收时间的概率分布，得到所述第二概率分布。
4.根据权利要求2或3所述的方法，其特征在于，所述根据所述峰值点信号强度的取值范围，确定所述第三概率分布，包括：
确定所述峰值点信号强度在所述峰值点信号强度的取值范围内满足均匀分布的情况下，所述峰值点信号强度的概率分布，得到所述第三概率分布。
5.根据权利要求2或3所述的方法，其特征在于，所述根据所述峰值点信号强度的取值范围，确定所述第三概率分布，包括：
确定在反射所述接收脉冲信号的对象表面的反射率的概率分布为均匀分布的情况下，所述峰值点信号强度在所述峰值点信号强度的取值范围内的概率分布，得到所述第三概率分布。
6.根据权利要求1-3任一项所述的方法，其特征在于，所述确定所述接收脉冲信号的脉冲宽度在脉冲宽度的取值范围内的第四概率分布，包括：
确定所述脉冲宽度在所述脉冲宽度的取值范围内满足均匀分布的情况下，所述脉冲宽度的概率分布，得到所述第四概率分布。
7.根据权利要求1-3任一项所述的方法，其特征在于，所述确定所述M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布，包括：
确定所述M个上升沿采样点中每个采样点在所述预设波形集合中的各波形上分布的概率分布；
根据所述M个上升沿采样点中每个采样点在所述预设波形集合中的各波形上分布的概率分布，确定所述第一概率分布。
8.根据权利要求7所述的方法，其特征在于，所述确定所述M个上升沿采样点中每个采样点在所述预设波形集合中的各波形上分布的概率分布，包括：
根据在采样目标采样点时噪声的强度的概率分布，确定所述目标采样点在所述预设波形集合中的各波形上分布的概率分布；其中，所述目标采样点包括所述M个上升沿采样点中任一个。
9.根据权利要求8所述的方法，其特征在于，所述目标采样点在所述预设波形集合中的各波形上分布的概率分布，满足以下公式一：
Figure FDA0003166567520000021
其中，pk表示所述目标采样点k在所述预设波形集合中目标波形上分布的概率，yk表示所述目标采样点k对应的信号强度，sk表示所述目标波形在采集所述目标采样点k时的信号强度，N0表示噪声功率。
10.根据权利要求1-3任一项所述的方法，其特征在于，所述根据所述第一概率分布和所述第五概率分布，确定所述接收脉冲信号为所述预设波形集合中的各波形的第六概率分布，包括：
根据所述第一概率分布以及所述第五概率分布，构建以所述预设波形集合为状态空间的马尔可夫链，得到所述马尔可夫链的样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}；其中，所述样本序列中各样本分别用于表征一种波形；所述样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}中的后n项样本{zm+1,...,zm+n-1,zm+n}用于表征所述第六概率分布；
所述根据所述第六概率分布，从所述预设波形集合中确定第一预设波形，包括：
计算所述后n项样本{zm+1,...,zm+n-1,zm+n}的期望，得到所述第一预设波形。
11.根据权利要求10所述的方法，其特征在于，所述根据所述第一概率分布以及所述第五概率分布，构建以所述预设波形集合为状态空间的马尔可夫链，包括：
利用Metropolis-Hastings算法，构建接受分布α(z,z')满足以下公式二的马尔可夫链：
Figure FDA0003166567520000031
其中，z为所述马尔可夫链的样本序列中的一项样本，z'为所述马尔可夫链的样本序列中z后一项样本，p(y|z')为所述第一概率分布中z'对应的概率，p(y|z)为所述第一概率分布中z对应的概率，p(z')为所述第五概率分布中z'对应的概率，p(z)为所述第五概率分布中z对应的概率，q(z',z)和q(z,z')为所述马尔可夫链的建议分布。
12.一种信号处理装置，其特征在于，包括：
采样单元，用于对接收脉冲信号进行N次采样，得到N个采样点，所述N个采样点中的每个采样点的采样参数包括所述采样点对应的接收时间和信号强度，所述N个采样点包括M个上升沿采样点以及P个削顶采样点，N、M和P分别为大于或等于1的整数；
第一概率确定单元，用于确定所述M个上升沿采样点在预设波形集合中的各波形上分布的第一概率分布；其中，所述预设波形集合中包括至少一个预设波形；
第二概率确定单元，用于根据所述P个削顶采样点确定所述接收脉冲信号的峰值点接收时间的第二概率分布和峰值点信号强度的第三概率分布；
所述第二概率确定单元，还用于确定所述接收脉冲信号的脉冲宽度在脉冲宽度的取值范围内的第四概率分布；
第三概率确定单元，用于根据所述第二概率分布、所述第三概率分布以及所述第四概率分布确定所述预设波形集合中的各波形为所述接收脉冲信号的第五概率分布；
第四概率确定单元，用于根据所述第一概率分布和所述第五概率分布，确定所述接收脉冲信号为所述预设波形集合中的各波形的第六概率分布；
波形确定单元，用于根据所述第六概率分布，从所述预设波形集合中确定第一预设波形；
参数确定单元，用于根据所述第一预设波形确定所述接收脉冲信号的峰值点接收时间和峰值点信号强度。
13.根据权利要求12所述的装置，其特征在于，所述第二概率确定单元，具体用于：
根据所述P个削顶采样点所对应的接收时间中的最小值，确定所述峰值点接收时间的取值范围的最小值；根据所述P个削顶采样点所对应的接收时间中的最大值，确定所述峰值点接收时间的取值范围的最大值；
根据所述峰值点接收时间的取值范围，确定所述第二概率分布；
根据所述P个削顶采样点对应的信号强度，确定所述峰值点信号强度的取值范围的最小值；
根据所述峰值点信号强度的取值范围，确定所述第三概率分布。
14.根据权利要求13所述的装置，其特征在于，所述根据所述峰值点接收时间的取值范围，确定所述第二概率分布，包括：确定所述峰值点接收时间在所述峰值点接收时间的取值范围内满足均匀分布的情况下，所述峰值点接收时间的概率分布，得到所述第二概率分布。
15.根据权利要求13或14所述的装置，其特征在于，所述根据所述峰值点信号强度的取值范围，确定所述第三概率分布，包括：确定所述峰值点信号强度在所述峰值点信号强度的取值范围内满足均匀分布的情况下，所述峰值点信号强度的概率分布，得到所述第三概率分布。
16.根据权利要求13或14所述的装置，其特征在于，所述根据所述峰值点信号强度的取值范围，确定所述第三概率分布，包括：确定在反射所述接收脉冲信号的对象表面的反射率的概率分布为均匀分布的情况下，所述峰值点信号强度在所述峰值点信号强度的取值范围内的概率分布，得到所述第三概率分布。
17.根据权利要求12-14任一项所述的装置，其特征在于，所述第二概率确定单元，具体用于确定所述脉冲宽度在所述脉冲宽度的取值范围内满足均匀分布的情况下，所述脉冲宽度的概率分布，得到所述第四概率分布。
18.根据权利要求12-14任一项所述的装置，其特征在于，所述第一概率确定单元，具体用于：
确定所述M个上升沿采样点中每个采样点在所述预设波形集合中的各波形上分布的概率分布；
根据所述M个上升沿采样点中每个采样点在所述预设波形集合中的各波形上分布的概率分布，确定所述第一概率分布。
19.根据权利要求18所述的装置，其特征在于，所述确定所述M个上升沿采样点中每个采样点在所述预设波形集合中的各波形上分布的概率分布，包括：
根据在采样目标采样点时噪声的强度的概率分布，确定所述目标采样点在所述预设波形集合中的各波形上分布的概率分布；其中，所述目标采样点包括所述M个上升沿采样点中任一个。
20.根据权利要求19所述的装置，其特征在于，所述目标采样点在所述预设波形集合中的各波形上分布的概率分布，满足以下公式一：
Figure FDA0003166567520000041
其中，pk表示所述目标采样点k在所述预设波形集合中目标波形上分布的概率，yk表示所述目标采样点k对应的信号强度，sk表示所述目标波形在采集所述目标采样点k时的信号强度，N0表示噪声功率。
21.根据权利要求12-14任一项所述的装置，其特征在于，所述第四概率确定单元，具体用于根据所述第一概率分布以及所述第五概率分布，构建以所述预设波形集合为状态空间的马尔可夫链，得到所述马尔可夫链的样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}；其中，所述样本序列中各样本分别用于表征一种波形；所述样本序列{z1,z2,...,zm,zm+1,...,zm+n-1,zm+n}中的后n项样本{zm+1,...,zm+n-1,zm+n}用于表征所述第六概率分布；
所述波形确定单元，具体用于计算所述后n项样本{zm+1,...,zm+n-1,zm+n}的期望，得到所述第一预设波形。
22.根据权利要求21所述的装置，其特征在于，所述根据所述第一概率分布以及所述第五概率分布，构建以所述预设波形集合为状态空间的马尔可夫链，包括：
利用Metropolis-Hastings算法，构建接受分布α(z,z')满足以下公式二的马尔可夫链：
Figure FDA0003166567520000051
其中，z为所述马尔可夫链的样本序列中的一项样本，z'为所述马尔可夫链的样本序列中z后一项样本，p(y|z')为所述第一概率分布中z'对应的概率，p(y|z)为所述第一概率分布中z对应的概率，p(z')为所述第五概率分布中z'对应的概率，p(z)为所述第五概率分布中z对应的概率，q(z',z)和q(z,z')为所述马尔可夫链的建议分布。
23.一种信号处理装置，其特征在于，所述信号处理装置包括一个或多个处理器，所述一个或多个处理器和一个或多个存储器耦合；所述一个或多个存储器存储有计算机指令；
当所述一个或多个处理器执行所述计算机指令时，使得所述信号处理装置执行如权利要求1-11中任一项所提供的信号处理方法。
24.一种芯片，其特征在于，所述芯片包括处理电路和接口；所述处理电路用于从存储介质中调用并运行所述存储介质中存储的计算机程序，以执行如权利要求1-11中任一项所提供的信号处理方法。
25.一种激光雷达，其特征在于，所述激光雷达包括权利要求12-23任一项所提供的信号处理装置，或者，所述激光雷达包括权利要求24所提供的芯片。
26.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有指令；当所述指令运行时，执行上述权利要求1-11任一项所提供的信号处理方法。
CN202080004821.6A 2020-07-20 2020-07-20 一种信号处理方法及装置 Active CN112639515B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2020/103124 WO2022016341A1 (zh) 2020-07-20 2020-07-20 一种信号处理方法及装置

Publications (2)

Publication Number Publication Date
CN112639515A CN112639515A (zh) 2021-04-09
CN112639515B true CN112639515B (zh) 2021-10-26

Family

ID=75291192

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202080004821.6A Active CN112639515B (zh) 2020-07-20 2020-07-20 一种信号处理方法及装置

Country Status (4)

Country Link
EP (1) EP4184212A4 (zh)
JP (1) JP2023534534A (zh)
CN (1) CN112639515B (zh)
WO (1) WO2022016341A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023183632A1 (en) * 2022-03-25 2023-09-28 Innovusion, Inc. A method for accurate time-of-flight calculation on saturated and non-saturated lidar receiving pulse data

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104181508A (zh) * 2014-08-31 2014-12-03 西安电子科技大学 基于压缩感知的威胁雷达信号检测方法
CN106840179A (zh) * 2017-03-07 2017-06-13 中国科学院合肥物质科学研究院 一种基于多传感器信息融合的智能车定位方法
US9954714B2 (en) * 2016-04-29 2018-04-24 The Boeing Company System and method for blindly acquiring frequency hopped spread spectrum signals
CN108732554A (zh) * 2017-04-18 2018-11-02 百度在线网络技术（北京）有限公司 激光雷达标定方法和装置
CN109061568A (zh) * 2018-04-28 2018-12-21 南京航空航天大学 一种基于目标特性的窄带雷达波形设计方法
CN109270427A (zh) * 2018-08-15 2019-01-25 南京国电南自电网自动化有限公司 一种基于相位比较的采样坏点识别和校正方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2001280758A1 (en) * 2000-07-24 2002-02-05 University Of Kansas Center For Research, Inc. Laser radar system and method
CN1988383A (zh) * 2006-12-22 2007-06-27 中兴通讯股份有限公司 一种信号处理中波形识别的方法
KR101078466B1 (ko) * 2007-01-19 2011-10-31 아사히 가세이 가부시키가이샤 각성 상태 판정 모델 생성 장치, 각성 상태 판정 장치, 경고 장치, 차량, 각성 상태 판정 모델 생성 방법, 각성 상태 판정 방법, 각성 상태 판정 모델 생성 프로그램을 기록한 컴퓨터 판독 가능 기록매체, 및 각성 상태 판정 프로그램을 기록한 컴퓨터 판독 가능 기록매체
US20090015447A1 (en) * 2007-03-16 2009-01-15 Daniel Kilbank Method for processing data using quantum system
CN104280725A (zh) * 2013-09-25 2015-01-14 中国科学院光电研究院 全波形激光雷达数据波形分解方法
CN105044702B (zh) * 2015-09-18 2017-06-13 宁波华仪宁创智能科技有限公司 脉冲波形的拟合方法
CN108732553B (zh) * 2018-06-01 2022-02-01 北京航空航天大学 一种激光雷达波形时刻鉴别方法与在线测距系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104181508A (zh) * 2014-08-31 2014-12-03 西安电子科技大学 基于压缩感知的威胁雷达信号检测方法
US9954714B2 (en) * 2016-04-29 2018-04-24 The Boeing Company System and method for blindly acquiring frequency hopped spread spectrum signals
CN106840179A (zh) * 2017-03-07 2017-06-13 中国科学院合肥物质科学研究院 一种基于多传感器信息融合的智能车定位方法
CN108732554A (zh) * 2017-04-18 2018-11-02 百度在线网络技术（北京）有限公司 激光雷达标定方法和装置
CN109061568A (zh) * 2018-04-28 2018-12-21 南京航空航天大学 一种基于目标特性的窄带雷达波形设计方法
CN109270427A (zh) * 2018-08-15 2019-01-25 南京国电南自电网自动化有限公司 一种基于相位比较的采样坏点识别和校正方法

Also Published As

Publication number Publication date
EP4184212A1 (en) 2023-05-24
JP2023534534A (ja) 2023-08-09
WO2022016341A1 (zh) 2022-01-27
EP4184212A4 (en) 2023-09-20
CN112639515A (zh) 2021-04-09

Similar Documents

Publication Publication Date Title
US10488495B2 (en) Single laser LIDAR system
EP3936895A1 (en) Distance measurement method, apparatus and device
US9110158B2 (en) Biphoton ranging with HOM interference
CN104502899A (zh) 一种自适应的恒虚警率目标检测方法
CN112639515B (zh) 一种信号处理方法及装置
US20230341529A1 (en) Target detection method, lidar and storage medium
CN113050071A (zh) 激光雷达数据处理方法、装置、设备及存储介质
US20110040525A1 (en) Method and system for calculating elementary symmetric functions of subsets of a set
CN116540250B (zh) 一种激光测距方法及系统
Ghani et al. Pulse repetition interval analysis using decimated Walsh-Hadamard transform
CN112904292A (zh) 基于重要性采样的宽带雷达检测门限确定及目标检测方法
US20200293860A1 (en) Classifying information using spiking neural network
CN115980689A (zh) 基于点云检测的辐射源信号分选方法、装置、设备及介质
WO2023111675A1 (en) Interference-resilient lidar waveform and estimation method thereof
CN113433563B (zh) 测距方法及测距装置
Ghani et al. Statistical modeling for missing and spurious pulses in pulse repetition interval (PRI) Analysis
RU2390724C2 (ru) Способ светолокационного определения дальности
Saeed et al. Design and implementation of a new CFAR based on weighted and statistical algorithms
KR102625993B1 (ko) 외부 간섭신호 처리 기반의 거리 측정 방법 및 그를 위한 장치
RU2791275C1 (ru) Способ обнаружения воздушных объектов при зондировании сверхкороткими радиоимпульсами
CN115186426B (zh) 射频探测系统主要目标识别全链路建模方法、装置、设备
CN116047533B (zh) 环境光干扰抑制方法、装置及其应用
CN117310636B (zh) 一种固定脉冲重复间隔测量方法、设备和介质
RU2292114C2 (ru) Способ управления порогом в оптимальном обнаружителе
CN117686995A (zh) 激光雷达回波时刻的鉴别方法、装置及通信设备

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