[发明专利]一种定位信号滤波方法、装置及设备

权利要求书

1.一种定位信号滤波方法，其特征在于，包括：

接收定位信号；

对所述定位信号进行信号转换处理得到连续N个相干积分结果，其中，N为正整数；

基于所得到的相干积分结果和预设的最大似然鉴别算法，并利用参数估计迭代算法，计算所述定位信号的残余载波频率和相位；

根据预设的卡尔曼滤波模型，对所述残余载波频率和相位进行平滑处理；

根据平滑处理后得到的残余载波频率和相位，对所述定位信号进行滤波处理；

所述对所述定位信号进行信号转换处理得到连续N个相干积分结果的步骤，包括：

对所述定位信号进行下变频处理，得到数字中频信号；

对所述数字中频信号进行捕获、跟踪处理，得到连续N个相干积分结果；

所述预设的最大似然鉴别算法，通过以下表达式表示：

其中，M(θ)表示最大似然代价函数，V_corr,n为第n次相干积分结果，T_coh为相干积分时间，Δf，分别为定位信号的残余载波频率和相位，real()和imag()分别为取实部函数和取虚部函数；

所述利用参数估计迭代方法计算所述定位信号的残余载波频率和相位的步骤，包括：

获得和λ的迭代初始值，为2×1的状态向量i为迭代次数i＝0,1,2……，λ为2×2的对角阵；

根据所获得的迭代初始值，获得矩阵H_i，其中，

判断H_i+λ是否为正定矩阵，如果否，则增大λ的值，直至H_i+λ为正定矩阵；

根据λ的值和预设的迭代最优解公式，确定的值，并判断是否满足如果否，则增大λ的值，直至满足其中，所述预设的迭代最优解公式为：G_i为2×1的梯度向量；

按照以下表达式，确定梯度向量G_i+1的值，作为目标值，其中，

当所述目标值小于第一预设阈值时，以预设数值为增量增大λ的值，并根据λ的值计算判断所计算的与的差值是否小于第二预设阈值，如果是，根据确定所述定位信号的残余载波频率和相位，如果否，令i＝i+1，并返回所述根据所获得的迭代初始值，获得矩阵H_i的步骤；

当所述目标值大于第三预设阈值时，按照所述预设数值减小λ的值，并根据λ的值计算判断所计算的与的差值是否小于所述第二预设阈值，如果是，根据确定所述定位信号的残余载波频率和相位，如果否，令i＝i+1，返回所述根据所获得的迭代初始值，获得H_i的步骤，其中，所述第三预设阈值大于所述第一预设阈值；

当所述梯度向量G_i+1的值小于第四预设阈值时，根据确定所述定位信号的残余载波频率和相位，或者当迭代次数i超过预设最大次数时，返回获得和λ的迭代初始值的步骤；

所述预先设定的卡尔曼滤波模型为以状态方程和观测方程描述的模型；其中，

所述状态方程为：X_k＝φX_k-1+W_k-1，

其中，X表示状态向量，X为：K-1表示时刻K的前一时刻，为载波相位估计值，ω_d＝2πΔf为角频率估计值，ω'_d为角频率变化率，W_k-1为输入的均值为零、方差为Q的高斯白噪声，φ为状态转移矩阵形式如下：

T＝NT_coh，T表示环路更新周期，T_coh为相干积分时间，

方差Q表示为:

所述观测方程为：Y_k＝HX_k+V_k，其中H为观测矩阵，形式如下：

V_k为均值为零、方差阵为R的观测噪声，

方差阵R表示为：式中I为单位矩阵。