一种利用GNSS海面反射信号相关功率探测目标的方法技术

一种利用GNSS海面反射信号相关功率探测目标的方法，它有两大步骤：步骤一：计算海面反射信号的时延-多普勒二维相关功率DDM；包括GNSS直射信号的捕获和跟踪；反射信号时延和多普勒频移估计；探测区域的划分和计算反射信号二维相关功率；步骤二：利用压缩感知理论对退化后的DDM进行恢复运算；包括建立基于压缩感知理论的退化后DDM的复原模型；DDM的稀疏性分解；选取高斯型随机矩阵作为测量矩阵和DDM的复原运算。本发明专利技术在遥感技术应用领域里有广阔地应用前景。

【技术实现步骤摘要】
一种利用GNSS海面反射信号相关功率探测目标的方法
本专利技术涉及一种基于GNSS海面反射信号的目标探测方法，更特别地说，是一种利用GNSS海面反射信号相关功率的时延-多普勒映射(DelayDopplerMapping，DDM)探测目标的方法，属于遥感技术应用领域。
技术介绍
目前的目标探测主要通过雷达来实现，利用定向天线将调制信号的电磁波能量发射到特定区域进行目标搜索，通过雷达接收机处理目标反射的回波，得到目标的特征信息。随着电子技术的不断发展，传统雷达受到电子干扰、反辐射导弹、各种隐身武器的威胁。现实环境中，存在着大量FM广播信号、通信基站信号、电台信号和GNSS卫星信号等，可为建立非合作式双(多)基地无源雷达提供丰富的外辐射源，以应对传统雷达所面临的威胁。GNSS为用户提供精确导航、定位和授时服务的同时，全天时、全天候地提供高度稳定、长期可用、安全性较高且覆盖全球的微波信号源。随着GNSS的不断发展，利用海面反射后的GNSS信号实现海面特征参数的提取是近年来GNSS-R(GNSSReflectometer)技术的研究热点。GNSS信号经海面反射后，信号的波形、极化、功率及频率等参量都受反射面的特征所影响，与反射面直接相关；因此可以利用GNSS卫星作为辐射源，建立非合作式双(多)基地无源雷达，通过对海面反射信号的处理来实现目标的探测和定位。在利用GNSS海面反射信号探测目标的双(多)基无源雷达系统中，当发射机和接收机速度、高度、天线增益等参数固定时，探测区域的反射信号相关功率主要由海面的雷达散射截面(RadarCrossSection，RCS)决定，而海面有无目标时的RCS差异较大，因此可以通过海面反射信号相关功率的时延-多普勒映射(DDM)来实现目标的探测与定位。
技术实现思路
本专利技术的目的是提供一种利用GNSS海面反射信号相关功率探测目标的方法，通过对海面反射区域的GNSS反射信号进行相关处理，得到反射信号的时延-多普勒二维相关功率分布，再利用图像处理的思想，将得到的DDM看作是原图像受到某种干扰后的退化图像，采用压缩感知理论对DDM进行恢复运算，得到退化之前的原图像，进而从原图像中提取出目标是否存在和目标位置的信息。一种利用GNSS海面反射信号相关功率探测目标的方法，包括以下几个步骤：步骤一：计算海面反射信号的时延-多普勒二维相关功率(DDM)；具体内容包括：(1)GNSS直射信号的捕获和跟踪；通过对直射信号的捕获跟踪，得到信号的C/A码码相位和载波多普勒频移，并完成导航定位解算，获得卫星和接收机的位置信息。(2)反射信号时延和多普勒频移估计；根据步骤(1)中卫星和接收机的位置信息，计算镜面反射点的位置，并估计出反射信号相对于直射信号的时延和多普勒频移。(3)探测区域的划分；以步骤(2)中镜面反射点的位置为中心，设定探测区域的范围，并对探测区域进行单元格划分。(4)计算反射信号二维相关功率；将步骤(3)中划分后的探测区域从空间域转换到时延-多普勒域，根据步骤(2)中反射信号相对于直射信号的时延和多普勒频移信息，利用步骤(1)中直射信号的码相位和载波多普勒辅助产生本地反射信号，并与接收到的反射信号进行相关运算，得到反射信号的DDM。步骤二：利用压缩感知理论对退化后的DDM进行恢复运算；具体内容包括：(1)建立基于压缩感知理论的退化后DDM的复原模型；在压缩感知理论框架下建立DDM退化模型，将模糊函数作为退化系统的冲击响应函数，把DDM转换为一个列向量；再建立DDM的复原模型，对原DDM进行估计。(2)DDM的稀疏性分解；对退化后的DDM进行小波变换，设定小波变换的层数，选择二维小波基矩阵作为稀疏基，对DDM进行稀疏表示。(3)选取高斯型随机矩阵作为测量矩阵；(4)DDM的复原运算；利用贪婪迭代算法对原始DDM稀疏系数进行估计，将其与二维小波基矩阵相乘，并对所得结果进行小波逆变换，得到原始DDM，再从时延-多普勒域转换到空间域，从反射信号相关功率分布的变化中可以探测到该区域是否有无目标，如果有目标，则可以得到目标所在单元格的位置信息。本专利技术的优点为：1、作为收发分置的双(多)基雷达，具有抗反辐射导弹攻击的能力，利用前向散射信号，可对隐形目标实现有效探测；2、将图像处理思想运用到时延-多普勒映射中，用于目标探测，较好的契合了反射信号DDM的特点；3、利用压缩感知理论对DDM进行复原运算，与传统的直接反卷积运算相比，减少了运算量和复杂度。附图说明图1是本专利技术的方法流程图；图2是本专利技术中探测区域划分实例图；图3是本专利技术中压缩感知理论框架下的DDM复原模型；图中符号说明如下：图1中的τE和fE分别表示反射信号相对于直射信号的时延和多普勒频移估计值；DDMde表示反射信号的二维相关功率分布；Φ为压缩感知理论框架下的测量矩阵；θ为对原始DDM稀疏分解后的稀疏系数矩阵，是一个一维的列向量矩阵；W表示对原始DDM稀疏分解的稀疏基矩阵，是一个二维小波基矩阵。图3中m,n分别表示时延和多普勒划分个数，为退化后的DDM，即DDMde；G(m,n)为通过测量矩阵Φ对DDMde在加性噪声的干扰下的观测结果；n表示加性噪声；H-1[·]表示估计过程；为对原始DDM的估计结果。具体实施方式下面结合附图对本专利技术的具体实施方式进行详细说明。本专利技术的一种利用GNSS海面反射信号相关功率探测目标的方法，流程如图1所示，包括以下几个步骤：步骤一：计算海面反射信号的时延-多普勒二维相关功率(DDM)；由于海面反射信号比较微弱，要直接捕获跟踪反射信号比较困难，本专利技术采用直射信号辅助反射信号的模式，首先对GNSS直射信号进行捕获跟踪，通过定位解算，得到卫星和接收处理平台的位置信息，计算出镜面反射点的位置，估计出反射信号相对于直射信号的时延和多普勒，结合直射信号的C/A码和载波信息，产生本地复制的反射信号；然后以镜面反射点为中心，对探测区域进行空间域的划分，并转换到时延和多普勒域，再与接收的反射信号进行相关处理，得到海面反射信号的时延-多普勒二维相关功率。本专利技术的反射信号相关功率计算方法具体包括：(1)GNSS直射信号的捕获跟踪；根据GNSS空间星座图和目标探测区域的几何关系选择要捕获跟踪的卫星，利用传统的方法，通过在时延和多普勒域并行搜索，完成GNSS直射信号的捕获，然后对其进行精确跟踪。接收到的直射信号S(t)表示为：S(t)＝A(t)D[t-τ(t)]C[t-τ(t)]cos[φ(t)](1)其中A(t)表示接收信号幅度电平，D[·]为导航电文数据，C[·]表示卫星的PRN码，φ(t)表示接收信号的载波相位，τ(t)表示码相位的空间传播延迟。在对GNSS直射信号精确跟踪后，进行导航定位解算，获得卫星和接收平台的位置和姿态信息。(2)反射信号时延和多普勒频移估计；根据步骤(1)中GNSS卫星和接收平台的位置和速度信息，可以得到接收平台相对于基准水平面的高度hR和卫星仰角利用反射事件的几何关系计算出镜面反射点的位置，并且估计出反射信号相对于直射信号的时延τE为其中c为光速；反射信号相对于直射信号的载波多普勒频移fE为fE＝[vt·m-vr·n-(vt-vr)utr]/λ(3)式中，vt和vr分别表示GNSS卫星和接收平台的运行速度；m本文档来自技高网...

【技术保护点】
一种利用GNSS海面反射信号相关功率探测目标的方法，其特征在于：包括以下几个步骤：步骤一：计算海面反射信号的时延‑多普勒二维相关功率DDM；(1)GNSS直射信号的捕获和跟踪；通过对直射信号的捕获跟踪，得到信号的C/A码码相位和载波多普勒频移，并完成导航定位解算，获得卫星和接收机的位置信息；(2)反射信号时延和多普勒频移估计；根据步骤(1)中卫星和接收机的位置信息，计算镜面反射点的位置，并估计出反射信号相对于直射信号的时延和多普勒频移；(3)探测区域的划分；以步骤(2)中镜面反射点的位置为中心，设定探测区域的范围，并对探测区域进行单元格划分；(4)计算反射信号二维相关功率；将步骤(3)中划分后的探测区域从空间域转换到时延‑多普勒域，根据步骤(2)中反射信号相对于直射信号的时延和多普勒频移信息，利用步骤(1)中直射信号的码相位和载波多普勒辅助产生本地反射信号，并与接收到的反射信号进行相关运算，得到反射信号的DDM；步骤二：利用压缩感知理论对退化后的DDM进行恢复运算；(1)建立基于压缩感知理论的退化后DDM的复原模型；在压缩感知理论框架下建立DDM退化模型，将模糊函数作为退化系统的冲击响应函数，把DDM转换为一个列向量；再建立DDM的复原模型，对原DDM进行估计；(2)DDM的稀疏性分解；对退化后的DDM进行小波变换，设定小波变换的层数，选择二维小波基矩阵作为稀疏基，对DDM进行稀疏表示；(3)选取高斯型随机矩阵作为测量矩阵；(4)DDM的复原运算；利用贪婪迭代算法对原始DDM稀疏系数进行估计，将其与二维小波基矩阵相乘，并对所得结果进行小波逆变换，得到原始DDM，再从时延‑多普勒域转换到空间域，从反射信号相关功率分布的变化中探测到该区域是否有无目标，如果有目标，则得到目标所在单元格的位置信息。...

【技术特征摘要】
1.一种利用GNSS海面反射信号相关功率探测目标的方法，其特征在于：包括以下几个步骤：步骤一：计算海面反射信号的时延-多普勒二维相关功率；(1)GNSS直射信号的捕获和跟踪；通过对直射信号的捕获跟踪，得到信号的C/A码码相位和载波多普勒频移，并完成导航定位解算，获得卫星和接收机的位置信息；(2)反射信号时延和多普勒频移估计；根据步骤(1)中卫星和接收机的位置信息，计算镜面反射点的位置，并估计出反射信号相对于直射信号的时延和多普勒频移；(3)探测区域的划分；以步骤(2)中镜面反射点的位置为中心，设定探测区域的范围，并对探测区域进行单元格划分；(4)计算反射信号二维相关功率；将步骤(3)中划分后的探测区域从空间域转换到时延-多普勒域，根据步骤(2)中反射信号相对于直射信号的时延和多普勒频移信息，利用步骤(1)中直射信号的码相位和载波多普勒频移辅助产生本地反射信号，并与接收到的反射信号进行相关运算，得到反射信号的时延-多普勒二维相关功率；步骤二：利用压缩感知理论对退化后...

【专利技术属性】
技术研发人员：杨东凯，江治东，李伟强，王强，王峰，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11