多普勒计程仪的测频方法技术

本发明专利技术公开了一种多普勒计程仪的测频方法，包括：从数据集S中选取s个点组成一个样本，并根据所述样本获得样本模型，其中，所述数据集S为所述多普勒计程仪接收到的回波信号的复相关或所述复相关中的一段；利用所述样本模型确定距离阈值dth内的数据点集Si，其中，Si是S的一个子集，所述Si中的数据点为内点；根据Si中的内点数量与阈值Nth的大小关系对所述样本模型进行重估，并在重估后的所述样板模型中估计直线以获得多普勒频移。通过仿真验证了本发明专利技术提供的多普勒计程仪的测频方法相比现有技术具有更高的稳定性，从而能够减小多普勒计程仪的测频误差。

Doppler log frequency measuring method

The invention discloses a method of frequency measurement, a Doppler log includes selecting s points form a sample from the S data set, and according to the sample to obtain the sample model, among them, the data set for the S echo signal received by Doppler log complex correlation or the complex period in determining the distance threshold; the DTH data set Si, using the sample model in which Si is a subset of the S, the Si of the data points for the interior point; the reassessment of the sample model according to the relation between the number of Si points and Nth threshold, and estimation straight to the Doppler shift in the model after the revaluation. Through the simulation of Doppler log frequency measuring method provided by the invention has higher stability compared with the existing technology, which can reduce the frequency measurement errors of Doppler log.

【技术实现步骤摘要】
多普勒计程仪的测频方法
本专利技术涉及多普勒计程仪
，特别是涉及一种多普勒计程仪的测频方法。
技术介绍
多普勒计程仪(DopplerVelocityLog，DVL)是一种测量船舶或水下平台相对于水底航速的仪器，其利用向水底发射声波和接收水底反射回波之间的多普勒频移，计算平台在横向和纵向的速度。多普勒计程仪进行信号处理的关键是多普勒频率的测量，其直接影响了测速精度。早期测频方法主要是过零点检测方法，包括双计数器过零点检测、数字式过零点检测方法等。当前测频主要采用复相关法，其主要依据是回波信号复相关的相位是多普勒频率的一次函数。根据采用的信号形式，多普勒计程仪可分为三类，一是脉冲非相干技术，即采用CW脉冲信号；二是脉冲相干技术，采用CW脉冲对；三是宽带编码技术，采用M序列或其它伪随机序列调制信。不论采用何种信号形式，对回波的复相关处理后，计算其相位时会得到多普勒频移的一次函数，对其斜率的估计误差将直接影响多普勒测频精度，从而影响测速精度。影响相位斜率估计精度的因素有三点，一是噪声的引入；二是波束开角造成区域多目标回波的叠加；三是复相关处理时对回波的截取方式，这涉及到回波的底跟踪处理。复相关相位斜率估计的传统方法是使用最小二乘法，其本质是使误差的均方根达到最小。最小二乘法适用于误差较小时，但当误差较大时其估计结果与真值差异较大，因此，现有技术的多普勒计程仪的测频方法存在测试误差较大的问题。
技术实现思路
本专利技术旨在至少解决现有技术中存在的技术问题之一。为此，本专利技术的目的在于提出一种测试误差小的多普勒计程仪的测频方法。本专利技术提供一种多普勒计程仪的测频方法，包括：从数据集S中选取s个点组成一个样本，并根据所述样本获得样本模型，其中，所述数据集S为所述多普勒计程仪接收到的回波信号的复相关或所述复相关中的一段；利用所述样本模型确定距离阈值dth内的数据点集Si，其中，Si是S的一个子集，所述Si中的数据点为内点；根据Si中的内点数量与阈值Nth的大小关系对所述样本模型进行重估，并在重估后的所述样板模型中估计直线以获得多普勒频移。根据本专利技术提供的多普勒计程仪的测频方法，采用RANSAC(RandomSampleConsensus，随机采样一致性)算法的技术，即从一个样本集中估计数学模型，可以适用于样本集中含有大量异常数据的情况，用于多普勒侧频中，对于复相关相位中存在较大误差、野值时有很好的稳定性，最终通过仿真验证了本专利技术提供的多普勒计程仪的测频方法相比现有技术具有更高的稳定性，从而能够减小多普勒计程仪的测频误差。另外，根据本专利技术上述的多普勒计程仪的测频方法，还可以具有如下附加的技术特征：进一步地，所述根据Si中的内点数量与阈值Nth的大小关系对所述样本模型进行重估，并在重估后的所述样板模型中估计直线以获得多普勒频移的步骤包括：若Si中的内点数量大于阈值Nth，则使用Si的所有内点的内点集重新估计模型，以获得重估后的所述样板模型；若Si中的内点数量小于阈值Nth，则选择一个新的样本子集，并判断新的样本子集中的内点数量是否大于阈值Nth，若否，则重复选择一个新的样本子集的步骤，直至新的样本子集中的内点数量大于阈值Nth。进一步地，所述方法还包括：从所述数据集S中选取选择两点组成一个样本，根据选择的两点确定一条直线；将到该直线距离小于所述距离阈值dth的点确定为当前内点，其余点设为当前外点，若所述当前内点的数量小于所述阈值Nth，则重新选择样本；若所述当前内点的数量大于所述阈值Nth，则取样结束，使用所述当前内点组成的内点集，并利用最小二乘法重新估计所述直线，估计得到的直线斜率即为所述多普勒频移。进一步地，所述复相关包括CW单脉冲的复相关。进一步地，所述CW单脉冲的复相关的获取方法包括以下步骤：当发射信号为CW单脉冲时，发射信号为：其中T为发射信号脉冲宽度，f0为信号频率，为初始相位；接收回波经过时延τ和多普勒频移fd，其信号为：其中，同相支路是回波信号与cos(2πf0t)相乘，具体为：正交支路是回波信号与-sin(2πf0t)相乘，具体为正交解调后信号中均包括一个高频项和一个低频项，低通滤波后，信号分别为：两路信号合并为复数信号，即：最终复相关为：进一步地，所述复相关包括CW脉冲对的复相关。进一步地，所述CW脉冲对的复相关的获取方法包括以下步骤：当发射声波是两个接近的CW脉冲组成的脉冲对时，发射信号为：其中T为发射信号脉冲宽度，f0为信号频率，为初始相位，Tτ为两个脉冲的起始时刻之间的时间间隔。接收回波经过时延τ和多普勒频移fd，并假设在短时内回波中两个脉冲的幅度仍然保持一致，且相位也没有差异，其信号形式可记为：经正交解调、低通滤波后形成的复数信号形式为：最终复相关为：附图说明本专利技术的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：图1是根据本专利技术一实施例的多普勒计程仪的测频方法的流程示意图；图2是发射波束几何示意图；图3是发射信号波形图；图4是接收回波信号波形图；图5是复相关相位的分布图；图6是严格提取互相关段落的估计速度的分布图；图7是模糊提取互相关段落的估计速度的分布图。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。请参阅图1，本专利技术一实施例提出的多普勒计程仪的测频方法，至少包括以下步骤：S101，从数据集S中选取s个点组成一个样本，并根据所述样本获得样本模型，其中，所述数据集S为所述多普勒计程仪接收到的回波信号的复相关或所述复相关中的一段；S102，利用所述样本模型确定距离阈值dth内的数据点集Si，其中，Si是S的一个子集，所述Si中的数据点为内点；具体的，Si即为该样本模型的一个支撑。S103，根据Si中的内点数量与阈值Nth的大小关系对所述样本模型进行重估，并在重估后的所述样板模型中估计直线以获得多普勒频移。上述方法其主要是采用RANSAC(RandomSampleConsensus，随机采样一致性)算法。作为一个具体示例，步骤S103中，所述根据Si中的内点数量与阈值Nth的大小关系对所述样本模型进行重估，并在重估后的所述样板模型中估计直线以获得多普勒频移的步骤包括：若Si中的内点数量大于阈值Nth，则使用Si的所有内点的内点集重新估计模型，以获得重估后的所述样板模型；若Si中的内点数量小于阈值Nth，则选择一个新的样本子集，并判断新的样本子集中的内点数量是否大于阈值Nth，若否，则重复选择一个新的样本子集的步骤，直至新的样本子集中的内点数量大于阈值Nth。具体实施时，可以在经过N次试验选择最大的支撑集Si，使用Si的所有点重估模型。作为一个具体示例，本专利技术一实施例还提出的一种多普勒计程仪的测频方法，包括：首先从所述数据集S中选取选择两点组成一个样本，根据选择的两点确定一条直线；将到该直线距离小于所述距离阈值dth的点确定为当前内点，其余点设为当前外点，若所述当前内点的数量小于所述阈值Nth本文档来自技高网...

【技术保护点】
一种多普勒计程仪的测频方法，其特征在于，包括：从数据集S中选取s个点组成一个样本，并根据所述样本获得样本模型，其中，所述数据集S为所述多普勒计程仪接收到的回波信号的复相关或所述复相关中的一段；利用所述样本模型确定距离阈值dth内的数据点集Si，其中，Si是S的一个子集，所述Si中的数据点为内点；根据Si中的内点数量与阈值Nth的大小关系对所述样本模型进行重估，并在重估后的所述样板模型中估计直线以获得多普勒频移。

【技术特征摘要】
1.一种多普勒计程仪的测频方法，其特征在于，包括：从数据集S中选取s个点组成一个样本，并根据所述样本获得样本模型，其中，所述数据集S为所述多普勒计程仪接收到的回波信号的复相关或所述复相关中的一段；利用所述样本模型确定距离阈值dth内的数据点集Si，其中，Si是S的一个子集，所述Si中的数据点为内点；根据Si中的内点数量与阈值Nth的大小关系对所述样本模型进行重估，并在重估后的所述样板模型中估计直线以获得多普勒频移。2.根据权利要求1所述的多普勒计程仪的测频方法，其特征在于，所述根据Si中的内点数量与阈值Nth的大小关系对所述样本模型进行重估，并在重估后的所述样板模型中估计直线以获得多普勒频移的步骤包括：若Si中的内点数量大于阈值Nth，则使用Si的所有内点的内点集重新估计模型，以获得重估后的所述样板模型；若Si中的内点数量小于阈值Nth，则选择一个新的样本子集，并判断新的样本子集中的内点数量是否大于阈值Nth，若否，则重复选择一个新的样本子集的步骤，直至新的样本子集中的内点数量大于阈值Nth。3.根据权利要求1所述的多普勒计程仪的测频方法，其特征在于，所述方法还包括：从所述数据集S中选取选择两点组成一个样本，根据选择的两点确定一条直线；将到该直线距离小于所述距离阈值dth的点确定为当前内点，其余点设为当前外点，若所述当前内点的数量小于所述阈值Nth，则重新选择样本；若所述当前内点的数量大于所述阈值Nth，则取样结束，使用所述当前内点组成的内点集，并利用最小二乘法重新估计所述直线，估计得到的直线斜率即为所述多普勒频移。4.根据权利要求1至3任意一项所述的多普勒计程仪的测频方法，其特征在于，所述复相关包括CW单脉冲的复相关。5.根据权利要求4所述的多普勒计程仪的测频方法，其特征在于，所述CW单脉冲的复相关的获取方法包括以下步骤：当发射信号为CW单脉冲时，发射信号为：其中T为发射信号脉冲宽度，f0为信号频率，为初始相位；接收回波经过时延τ和多普勒频移fd，其信号为：其中，同相支路是回波信号与cos(2πf0t)相乘，具体为：正交支路是回波信号与-sin(2πf0t)相乘，具体为正交解调后信号中均包括一个高频项和一个...

【专利技术属性】
技术研发人员：杨海东，邱震钰，白书华，郭鹏飞，黄斌，
申请(专利权)人：南昌理工学院，
类型：发明
国别省市：江西,36