一种基于压缩感知算法的多径参数估计方法技术

本发明专利技术公开了一种基于压缩感知算法的多径参数估计方法。使用本发明专利技术能够实现较高的参数估计精度，且算法收敛快、性能稳定。本发明专利技术利用多径信号在时间轴上的稀疏特性，基于压缩感知算法，对多径参数进行估计，有效克服传统最小二乘估计方法的无偏性局限，明显改善对病态数据的拟合效果，具备较高的数字稳定性，相比传统最大似然估计方法的噪声敏感性缺陷，采用的压缩感知方法能够有效抑制噪声影响，稳定数字估计，故能实现较高的参数估计精度，且算法收敛快、性能稳定，克服了信号域多径信号参数估计技术对多路模型的多径信号分离/消除效能不高的问题。

【技术实现步骤摘要】
一种基于压缩感知算法的多径参数估计方法
本专利技术属于卫星导航
，涉及卫星/伪卫星导航信号测距技术，具体涉及一种基于压缩感知算法的多径信号参数估计方法。
技术介绍
多径误差作为目前卫星导航定位系统及其增强系统的主要误差来源，严重影响接收机伪距测量准确度，是高精度定位亟待解决的技术难题。多径抑制手段，在终端设计层面，体现于天线设计、信号域基带算法设计、测量域数据处理等三大阶段。其中，天线设计只能有效处理地面所反射的多径信号，而对来自天线上方的多径信号意义不大；测量域数据处理需基于长时观测，利用多径信号在时间和空间上的弱相关特性予以处理，不适用实时处理；信号域基带算法设计基于信号相关函数特性检测多径信号，既保证锁定各形态多径信号又兼顾实时性需求。信号域基带算法设计主要分为码相关参考波形(CCRW)和多径估计技术。码相关参考波形，诸如：窄相关、Double-Delta、Strobe、HRC技术等，均基于无限带宽假设前提设计抗多径算法，其性能在有限带宽的现实应用条件下，难以充分发挥性能优势；多径估计技术，最具代表的多径延迟锁定环(MEDLL)和多径消除技术(MMT)，基于多径信号模型假设，以最大似然估计为准则检测估计多径信号参数，以分离单/多路强多径信号，效能优势明显，但算法复杂度和硬件资源占用，会随多径信号路径数目的增加而急剧上升，故一般应用于单路径强多径信号的检测和抑制。
技术实现思路
有鉴于此，本专利技术提供了一种基于压缩感知算法的多径参数估计方法，能够实现较高的参数估计精度，且算法收敛快、性能稳定。本专利技术的基于压缩感知算法的多径参数估计方法，包括如下步骤：步骤1，基于伪随机码多相关器对接收到的导航信号自相关函数进行采样；根据多相关器采样点的自相关函数幅值，构建多径分量稀疏基矩阵方程；其中，多相关器的数量由信道传播情况确定；步骤2，对步骤1构建的稀疏基矩阵方程进行求解，获得多径信号幅值估计值；步骤3，对步骤2获得的多径信号幅值估计值进行判断，确定多径信号的有无及位置、幅值；其中，A)若多径信号幅值则判定在对应的样值估计点βm处无多径信号；m＝1,2，……，M，M为多相关器的个数；B)若且及则判定在样值估计点βm处存在多径信号，该多径信号幅值为C)若且及都接近于0(即小于或等于10-2量级)，则判定在样值估计点βm与βm+1之间，存在一条多径信号；该多径信号的位置及幅值由βm与βm+1样值点的位置、幅值线性拟合确定；D)若连续多个均大于0，且其余的多径信号幅值均小于0，则降低当前多相关器样值点采样周期，重复步骤1～4，直至获得多径信号的位置及幅值。进一步的，所述步骤2包括如下子步骤：步骤2.1，利用凸优化算法根据基追踪/基追踪去噪算法模型，将步骤1构建的多径分量稀疏基矩阵方程的多径信号幅值求解问题转化为最小化L1范数的优化求解问题，并获得多径信号幅值的最优解形式；步骤2.2，根据步骤2.1的最优解形式，将最优解的求解问题转化为一个边界约束最优化的二次规划问题，并进行求解，获得多径信号幅值估计值。有益效果：(1)多径参数估计精度高本专利技术利用多径信号在时间轴上的稀疏特性，基于压缩感知算法，对多径参数进行估计，有效克服传统最小二乘估计方法的无偏性局限，明显改善对病态数据的拟合效果，具备较高的数字稳定性，故能实现较高的参数估计精度，克服了信号域多径信号参数估计技术对多路模型的多径信号分离/消除效能不高的问题。(2)抗噪性能优势明显本专利技术所设计压缩感知重构模型，能有效对抗稀疏基观测矩阵多重共线性(矩阵病态)，相比传统最大似然估计方法的噪声敏感性缺陷，压缩感知方法能够有效抑制噪声影响，稳定数字估计。(3)算法高效本专利技术算法收敛快且性能稳定，可基于少量观测数据实现多径信号重构，并保证较高估计精度。(4)本专利技术利用BP/BPDN模型并结合二次规划对构建的稀疏基矩阵方程的求解过程进行优化，求解速度快，且精度高。附图说明图1为本专利技术基于压缩感知算法的多径参数估计方法结构图。图2为本专利技术基于压缩感知算法的多径参数估计方法流程图。图3为多径信号自相关函数波形示意图。图4为多相关器采样示意图。具体实施方式下面结合附图并举实施例，对本专利技术进行详细描述。本专利技术提供了一种基于压缩感知算法的多径参数估计方法，基于多径信号信道传播模型，利用压缩感知算法，刻画多径信号时间轴稀疏特性，估计多径信号参数，依据模型假设分离多径分量，实现多径抑制/消除。本专利技术多径参数估计方法的框架图和原理图分别如图1和图2所示，具体步骤如下：步骤1，接收机终端天线接收导航信号，基于伪随机码多相关器对信号自相关函数进行采样，根据采样点信息构建多径信号分量估计模型及稀疏基矩阵方程。具体的，导航信号受传播信号物理特征影响，诸如：地形地貌、周围建筑物等，发生散射和反射现象，从而引发多径效应。为简化分析，不考虑导航电文数据影响，接收机接收到的信号表达式为：其中，A为直达信号幅值，C为信号所调制伪随机码序列，f0为信号载波标称频点，fd为信号多普勒频移，τ为直达信号路径传播时延，为直达信号载波相位，n为多径信号索引号，N为多径信号路径数目，an为第n路多径信号幅值，τn为第n路多径信号延迟，为第n路多径信号载波相位，noise为信道高斯白噪声，t为时间。接收机对接收信号进行去载波操作，并在本地复现多路伪随机码，基于多相关器原理对信号自相关函数(如图3所示)进行时延采样(见图4所示)，得到信号正交解调后相干积分结果表达式为：其中，fres为残余载波频率，m为采样点索引，M为多相关器数目(相关函数采样点数)，βm为采样点对应时延，R为相关函数，sinc为辛格函数。由于多径信号分量呈时延簇分布，各路径时延间隔远大于时延搜索间隔(多相关器间隔)，使得时延在时域基满足稀疏特性，由此，可利用压缩感知算法对多径分量参数进行模型估计。根据多相关器采样点相关函数幅值，构建多径分量稀疏基矩阵方程，其表达式为：R＝Ha+V(3)其中，a为可能的多径信号幅值，表达式为：a＝[a1,a2,…,aM](4)H为可能的多径信号观测矩阵(稀疏基矩阵)，表达式为：V为信道估计噪声，服从均值为0，方差为σ2的加性高斯白噪声：V＝[v1,v2,…,vM](6)可知，多径信号恢复重构问题可转化为求解线性方程组(3)中向量a的问题。上述量测模型(3)，可以根据信道传播情况、待估计多径参数精度需求，配置多径信号路径数目及相关器个数，加以调整，其中，多径信号路径数越多，需要的相关器个数就越多，这样估计出来的路径参数才最精确。步骤2，对步骤1构建的稀疏基矩阵方程进行求解，获得多径信号幅值估计值；其中，可以利用BP/BPDN模型联合二次规划对构建的稀疏基矩阵方程的求解过程进行优化，求解速度快，且精度高，具体包括如下子步骤：步骤2.1，根据约束条件调整，优化稀疏基矩阵方程求解形式。由于上述稀疏基矩阵方程(式(3))为欠定方程组，故存在无穷多解。基于稀疏重构理论，矩阵H满足约束等距性(RestrictedIsometryProperty,RIP)等稀疏重构条件，可对公式(3)做进一步改进。利用凸优化算法根据基追踪/基追踪去噪算法(BasisPursuit/BasisPursuitDe-Noising,BP/BPDN)模型，本文档来自技高网...

【技术保护点】
1.一种基于压缩感知算法的多径参数估计方法，其特征在于，包括如下步骤：步骤1，基于伪随机码多相关器对接收到的导航信号自相关函数进行采样；根据多相关器采样点的自相关函数幅值，构建多径分量稀疏基矩阵方程；其中，多相关器的数量由信道传播情况确定；步骤2，对步骤1构建的稀疏基矩阵方程进行求解，获得多径信号幅值估计值；步骤3，对步骤2获得的多径信号幅值估计值进行判断，确定多径信号的有无及位置、幅值；其中，A)若多径信号幅值

【技术特征摘要】
1.一种基于压缩感知算法的多径参数估计方法，其特征在于，包括如下步骤：步骤1，基于伪随机码多相关器对接收到的导航信号自相关函数进行采样；根据多相关器采样点的自相关函数幅值，构建多径分量稀疏基矩阵方程；其中，多相关器的数量由信道传播情况确定；步骤2，对步骤1构建的稀疏基矩阵方程进行求解，获得多径信号幅值估计值；步骤3，对步骤2获得的多径信号幅值估计值进行判断，确定多径信号的有无及位置、幅值；其中，A)若多径信号幅值则判定在对应的样值估计点βm处无多径信号；m＝1,2，……，M，M为多相关器的个数；B)若且及则判定在样值估计点βm处存在多径信号，该多径信号幅值为C)若且及都接近于0，则判定在样值估计点βm与βm+1之间，存在一...

【专利技术属性】
技术研发人员：罗瑞丹，徐颖，袁洪，袁杰，史雨薇，唐阳阳，
申请(专利权)人：中国科学院光电研究院，
类型：发明
国别省市：北京,11