当前位置: 首页 > 专利查询>武汉大学专利>正文

一种基于子空间技术的海杂波抑制方法技术

技术编号:19567161 阅读:51 留言:0更新日期:2018-11-25 02:36
本发明专利技术提出了一种基于子空间技术的海杂波抑制方法。本发明专利技术根据雷达系统单极子的距离多普勒谱矩阵统计单极子交叉环组阵中高信噪比谱点的位置向量;通过交叉环内对应的幅度校准值以及相位校准值对第一单极子以及第二单极子交叉环环天线的距离多普勒谱校准;通过单极子间对应的幅度校准值以及相位校准值构建校准后单极子交叉环组阵的距离多普勒谱向量;海杂波正负一阶谱区域内通过正负一阶峰边界构建正负一阶峰区域;通过校准后单极子交叉环组阵的距离多普勒谱向量构建第一以及第二正交投影矩阵,并对第一单极子交叉环以及第二单极子的距离多普勒谱检矩阵进行校准。本发明专利技术优点在于对海杂波进行有效抑制。

【技术实现步骤摘要】
一种基于子空间技术的海杂波抑制方法
本专利技术属于船载高频地波雷达杂波抑制领域,具体涉及一种基于子空间技术的海杂波抑制方法。
技术介绍
高频地波雷达(HFGWR)基于高频电磁波(3~30MHz)在海面绕射的传播特性可对指定海域进行全天候的超视距实时监测,及时发现作用范围内的海面船只与空中运动目标,并提供其距离,方位,速度等信息,受到了沿海地区的高度重视。高频地波雷达根据架设环境可分为船载和岸基,与传统的岸基高频地波雷达相比,船载不仅继承了岸基高频地波雷达的优点,还具有较强的灵活性和机动性,能够随时测量感兴趣海域和目标,不受架设场地制约,在海面目标检测和海洋气象研究等方面更具开发潜力。船载雷达相对于岸基雷达在目标检测和追踪方面有自己独特的优势,但同时也还存在特殊问题有待解决。岸基雷达海浪回波理论上多普勒频率仅仅只是两个离散的点,由于流速的影响发生了比较小的展宽,大部分船只目标均不会被海杂波淹没。对于船载雷达而言,由于平台的移动引起海洋回波的多普勒展宽,使得许多船只目标很容易被海洋一阶谱淹没,增大了目标检测的难度,给船载雷达信号处理带来了新的挑战。由于船载平台的运动让不同的方位角海浪具有不同的多普勒频移,使得海杂波的方位角与多普勒频率发生耦合,许多雷达工作者均提倡空-时自适应处理来解决杂波抑制问题。空时自适应处理(STAP)利用船载海洋回波多普勒与方位角耦合的特性,对回波进行二维联合滤波处理。但是空时自适应滤波效果取决于对杂波协方差的估计,在船载HFGWR当中应用时,受限于参考单元的匮乏,无法准确估计协方差,导致不能出现杂波凹口,滤波性能严重下降。正交加权法能够对各个通道的多普勒进行正交加权波束形成干扰方向正交加权输出为零,出现凹口。该方法大大降低了协方差维数,能够有效的抑制海杂波,但在进行目标跟踪时,方位角是由波束形成得到,其准确性主要依赖于天线口径,HFSWR雷达波长较长,这就对船载平台提出了较高要求,一般为舰载,限制了船载HFSWR雷达的应用范围。因此,一种适用于便携式HFGWR的海杂波抑制方法是尤为需要的。
技术实现思路
为了解决船载雷达目标落入一阶海杂波谱内检测困难的问题,本专利技术提出了一种基于子空间技术的海杂波抑制方法。本专利技术的技术方案为一种基于子空间技术的海杂波抑制方法,具体包括以下步骤:步骤1:构建雷达系统单极子交叉环组阵的距离多普勒谱矩阵向量,根据单极子交叉环组阵的距离多普勒谱矩阵中第一单极子或第二单极子的距离多普勒谱矩阵统计出高信噪比的谱点在相应单极子的距离多普勒谱矩阵中的位置向量;步骤2:针对步骤1得到的高信噪比的谱点的位置向量,对每个单极子交叉环的距离多普勒谱矩阵进行校准。统计每个单极子交叉环a环天线的距离多普勒谱矩阵在位置向量的元素与单极子的距离多普勒谱矩阵在位置向量的元素幅度比与相位差,计算得到每个单极子交叉环中a环天线的幅度校准值与相位校准值,通过该幅度校准值与相位校准值对每个单极子交叉环中a环天线的距离多普勒谱矩阵进行校准。统计每个单极子交叉环b环天线的距离多普勒谱矩阵在位置向量的元素与单极子的距离多普勒谱矩阵在位置向量的元素幅度比与相位差,计算得到每个单极子交叉环中b环天线的幅度校准值与相位校准值。通过该幅度校准值与相位校准值对每个单极子交叉环中b环天线的距离多普勒谱矩阵进行校准;步骤3:根据第一单极子交叉环与第二单极子的幅度比值计算单极子交叉环间的幅度校准值,根据第一单极子交叉环与第二单极子的相位差结合AIS系统计算单极子交叉环间的相位校准值,根据单极子交叉环间的幅度校准值以及相位校准值对第一单极子交叉环以及第二单极子的距离多普勒谱进行校准,根据单极子交叉环间的幅度校准值以及相位校准值对校准后第一单极子交叉环以及第二单极子的环天线的距离多普勒谱进行校准,并构建校准后单极子交叉环组阵的距离多普勒谱向量;步骤4:从第一单极子的距离多普勒谱以及第二单极子的距离多普勒谱中任意取一个单极子的距离多普勒谱的行向量,进行数据平滑以及差谱划分处理,在大致海杂波正负一阶谱区域分别得到正一阶峰边界以及负一阶峰边界,通过正一阶峰边界以及负一阶峰边界构建正负一阶峰区域;步骤5:通过校准后单极子交叉环组阵的距离多普勒谱矩阵向量分别构建第一单极子的距离多普勒谱矩阵以及第二单极子的距离多普勒谱矩阵,第一单极子的距离多普勒谱矩阵在正负一阶峰区域范围内根据相邻参考单元计算第一海杂波空域子空间,第二单极子的距离多普勒谱矩阵在正负一阶峰区域内根据相邻参考单元计算第二海杂波空域子空间,通过第一海杂波空域子空间构建第一正交投影矩阵,通过第二海杂波空域子空间构建第二正交投影矩阵,通过第一正交投影矩阵对第一单极子交叉环以及第二单极子的距离多普勒谱检矩阵进行校准,或通过第二正交投影矩阵对第一单极子交叉环以及第二单极子的距离多普勒谱检矩阵进行校准。作为优选,步骤1中所述单极子交叉环组阵的距离多普勒谱矩阵向量为:Stotal=[Smono1,Smono2,Sringa1,Sringa2,Sringb1,Sringb2]其中,Smono1表示第一单极子交叉环中单极子的距离多普勒谱矩阵,Smono2表示第二单极子交叉环中单极子的距离多普勒谱矩阵,Sringa1表示第一单极子交叉环中a环天线的距离多普勒谱矩阵,Sringa2表示第二单极子交叉环中a环天线的距离多普勒谱矩阵,Sringb1表示第一单极子交叉环中b环天线的距离多普勒谱矩阵,Sringb2表示第二单极子交叉环中b环天线的距离多普勒谱矩阵;步骤1中所述单极子交叉环组阵的距离多普勒谱向量中每个元素可以表示为:其中,R表示距离元的数量,N表示多普勒频点的数量,Sr,n(l)表示Stotal中第l个元素第r个距离元第n个多普勒频点的多普勒数据,r∈[1,R],n∈[1,N];步骤1中所述第一单极子的距离多普勒谱矩阵为Stotal(1)行向量为第二单极子的距离多普勒谱矩阵为Stotal(2)。步骤1中所述统计出高信噪比的谱点在距离多普勒谱矩阵Stotal(l)中的位置向量Pos,该位置向量Pos可从Stotal(1)或Stotal(2)中统计得到,以Stotal(1)为例,得到的位置向量Pos为:统计Stotal(1)满足条件20*log10(||Sr,n(1)||)-σ>high的元素在第一单极子的距离多普勒谱矩阵中的位置,将这些满足条件的元素的位置(r,n)组成位置向量Pos,其中σ为噪声功率,high为人为设定的信噪比门限,单位为dB;作为优选,步骤2中所述根据每个单极子交叉环a环天线的幅度比值计算每个单极子交叉环a环天线的幅度校准值以第一个单极子交叉环为例,a环天线的幅度校准值为:其中,Num为高信噪比谱点的向量位置Pos的元素数量。步骤2中所述根据第一单极子交叉环环天线的相位差计算第一单极子交叉环环天线的相位校准值:针对第一单极子交叉环高信噪比谱点的行向量位置序列Pos,根据第一单极子交叉环中a环天线的距离多普勒谱矩阵元素与第一单极子的距离多普勒谱相位差计算第一单极子交叉环a环天线的相位校准值:P1,a=mode(angle(SPos(3))-angle(SPos(1)))其中,mode(·)表示求取序列的众数,angle(·)表示求取相位针对第一单极子交叉环高本文档来自技高网
...

【技术保护点】
1.一种基于子空间技术的海杂波抑制方法,其特征在于包括以下步骤:步骤1:构建雷达系统单极子交叉环组阵的距离多普勒谱矩阵向量,根据单极子交叉环组阵的距离多普勒谱矩阵中第一单极子或第二单极子的距离多普勒谱矩阵统计出高信噪比的谱点在相应单极子的距离多普勒谱矩阵中的位置向量;步骤2:针对步骤1得到的高信噪比的谱点的位置向量,对每个单极子交叉环的距离多普勒谱矩阵进行校准,统计每个单极子交叉环a环天线的距离多普勒谱矩阵在位置向量的元素与单极子的距离多普勒谱矩阵在位置向量的元素幅度比与相位差,计算得到每个单极子交叉环中a环天线的幅度校准值与相位校准值,通过该幅度校准值与相位校准值对每个单极子交叉环中a环天线的距离多普勒谱矩阵进行校准,统计每个单极子交叉环b环天线的距离多普勒谱矩阵在位置向量的元素与单极子的距离多普勒谱矩阵在位置向量的元素幅度比与相位差,计算得到每个单极子交叉环中b环天线的幅度校准值与相位校准值,通过该幅度校准值与相位校准值对每个单极子交叉环中b环天线的距离多普勒谱矩阵进行校准;步骤3:根据第一单极子交叉环与第二单极子的幅度比值计算单极子交叉环间的幅度校准值,根据第一单极子交叉环与第二单极子的相位差结合AIS系统计算单极子交叉环间的相位校准值,根据单极子交叉环间的幅度校准值以及相位校准值对第一单极子交叉环以及第二单极子的距离多普勒谱进行校准,根据单极子交叉环间的幅度校准值以及相位校准值对校准后第一单极子交叉环以及第二单极子的环天线的距离多普勒谱进行校准,并构建校准后单极子交叉环组阵的距离多普勒谱向量;步骤4:从第一单极子的距离多普勒谱以及第二单极子的距离多普勒谱中任意取一个单极子的距离多普勒谱的行向量,进行数据平滑以及差谱划分处理,在大致海杂波正负一阶谱区域分别得到正一阶峰边界以及负一阶峰边界,通过正一阶峰边界以及负一阶峰边界构建正负一阶峰区域;步骤5:通过校准后单极子交叉环组阵的距离多普勒谱矩阵向量分别构建第一单极子的距离多普勒谱矩阵以及第二单极子的距离多普勒谱矩阵,第一单极子的距离多普勒谱矩阵在正负一阶峰区域范围内根据相邻参考单元计算第一海杂波空域子空间,第二单极子的距离多普勒谱矩阵在正负一阶峰区域内根据相邻参考单元计算第二海杂波空域子空间,通过第一海杂波空域子空间构建第一正交投影矩阵,通过第二海杂波空域子空间构建第二正交投影矩阵,通过第一正交投影矩阵对第一单极子交叉环以及第二单极子的距离多普勒谱检矩阵进行校准,或通过第二正交投影矩阵对第一单极子交叉环以及第二单极子的距离多普勒谱检矩阵进行校准。...

【技术特征摘要】
1.一种基于子空间技术的海杂波抑制方法,其特征在于包括以下步骤:步骤1:构建雷达系统单极子交叉环组阵的距离多普勒谱矩阵向量,根据单极子交叉环组阵的距离多普勒谱矩阵中第一单极子或第二单极子的距离多普勒谱矩阵统计出高信噪比的谱点在相应单极子的距离多普勒谱矩阵中的位置向量;步骤2:针对步骤1得到的高信噪比的谱点的位置向量,对每个单极子交叉环的距离多普勒谱矩阵进行校准,统计每个单极子交叉环a环天线的距离多普勒谱矩阵在位置向量的元素与单极子的距离多普勒谱矩阵在位置向量的元素幅度比与相位差,计算得到每个单极子交叉环中a环天线的幅度校准值与相位校准值,通过该幅度校准值与相位校准值对每个单极子交叉环中a环天线的距离多普勒谱矩阵进行校准,统计每个单极子交叉环b环天线的距离多普勒谱矩阵在位置向量的元素与单极子的距离多普勒谱矩阵在位置向量的元素幅度比与相位差,计算得到每个单极子交叉环中b环天线的幅度校准值与相位校准值,通过该幅度校准值与相位校准值对每个单极子交叉环中b环天线的距离多普勒谱矩阵进行校准;步骤3:根据第一单极子交叉环与第二单极子的幅度比值计算单极子交叉环间的幅度校准值,根据第一单极子交叉环与第二单极子的相位差结合AIS系统计算单极子交叉环间的相位校准值,根据单极子交叉环间的幅度校准值以及相位校准值对第一单极子交叉环以及第二单极子的距离多普勒谱进行校准,根据单极子交叉环间的幅度校准值以及相位校准值对校准后第一单极子交叉环以及第二单极子的环天线的距离多普勒谱进行校准,并构建校准后单极子交叉环组阵的距离多普勒谱向量;步骤4:从第一单极子的距离多普勒谱以及第二单极子的距离多普勒谱中任意取一个单极子的距离多普勒谱的行向量,进行数据平滑以及差谱划分处理,在大致海杂波正负一阶谱区域分别得到正一阶峰边界以及负一阶峰边界,通过正一阶峰边界以及负一阶峰边界构建正负一阶峰区域;步骤5:通过校准后单极子交叉环组阵的距离多普勒谱矩阵向量分别构建第一单极子的距离多普勒谱矩阵以及第二单极子的距离多普勒谱矩阵,第一单极子的距离多普勒谱矩阵在正负一阶峰区域范围内根据相邻参考单元计算第一海杂波空域子空间,第二单极子的距离多普勒谱矩阵在正负一阶峰区域内根据相邻参考单元计算第二海杂波空域子空间,通过第一海杂波空域子空间构建第一正交投影矩阵,通过第二海杂波空域子空间构建第二正交投影矩阵,通过第一正交投影矩阵对第一单极子交叉环以及第二单极子的距离多普勒谱检矩阵进行校准,或通过第二正交投影矩阵对第一单极子交叉环以及第二单极子的距离多普勒谱检矩阵进行校准。2.根据权利要求1所述的基于子空间技术的海杂波抑制方法,其特征在于:步骤1中所述单极子交叉环组阵的距离多普勒谱矩阵向量为:Stotal=[Smono1,Smono2,Sringa1,Sringa2,Sringb1,Sringb2]其中,Smono1表示第一单极子交叉环中单极子的距离多普勒谱矩阵,Smono2表示第二单极子交叉环中单极子的距离多普勒谱矩阵,Sringa1表示第一单极子交叉环中a环天线的距离多普勒谱矩阵,Sringa2表示第二单极子交叉环中a环天线的距离多普勒谱矩阵,Sringb1表示第一单极子交叉环中b环天线的距离多普勒谱矩阵,Sringb2表示第二单极子交叉环中b环天线的距离多普勒谱矩阵;步骤1中所述单极子交叉环组阵的距离多普勒谱向量中每个元素可以表示为:其中,R表示距离元的数量,N表示多普勒频点的数量,Sr,n(l)表示Stotal中第l个元素第r个距离元第n个多普勒频点的多普勒数据,r∈[1,R],n∈[1,N];步骤1中所述第一单极子的距离多普勒谱矩阵为Stotal(1)行向量为第二单极子的距离多普勒谱矩阵为Stotal(2);步骤1中所述统计出高信噪比的谱点在距离多普勒谱矩阵Stotal(l)中的位置向量Pos,该位置向量Pos可从Stotal(1)或Stotal(2)中统计得到,以Stotal(1)为例,得到的位置向量Pos为:统计Stotal(1)满足条件的元素在第一单极子的距离多普勒谱矩阵中的位置,将这些满足条件的元素的位置(r,n)组成位置向量Pos,其中σ为噪声功率,high为人为设定的信噪比门限,单位为dB。3.根据权利要求1所述的基于子空间技术的海杂波抑制方法,其特征在于:步骤2中所述根据每个单极子交叉环a环天线的幅度比值计算每个单极子交叉环a环天线的幅度校准值以第一个单极子交叉环为例,a环天线的幅度校准值为:其中,Num为高信噪比谱点的向量位置Pos的元素数量;步骤2中所述根据第一单极子交叉环环天线的相位差计算第一单极子交叉环环天线的相位校准值:针对第一单极子交叉环高信噪比谱点的行向量位置序列Pos,根据第一单极子交叉环中a环天线的距离多普勒谱矩阵元素与第一单极子的距离多普勒谱相位差计算第一单极子交叉环a环天线的相位校准值:P1,a=mode(angle(SPos(3))-angle(SPos(1)))其中,mode(·)表示求取序列的众数,angle(·)表示求取相位针对第一单极子交叉环高信噪比谱点的位置向量Pos,根据第一单极子交叉环中b环天线的距离多普勒谱矩阵的Pos的元素与第一单极子的距离多普勒谱矩阵的Pos的元素幅度比值计算第一单极子交叉环b环天线的幅度校准值:针对第一单极子交叉环高信噪比谱点的行向量位置序列Pos,根据第一单极子交叉环中b环天线的距离多普勒谱的行向量与第一单极子的距离多普勒谱的行向量相位差计算第一单极子交叉环b环天线的相位校准值:P1,b=mode(angle(SPos(5))-angle(SPos(1)))步骤2中所述根据第二单极子交叉环环天线的幅度校准值与相位校准值的计算方法和第二单极子交叉环相同,幅度校准值与相位校准值为:第二单极子交叉环a环天线幅度校准值为:第二单极子交叉环a环天线相位校准值为:P2,a=mode(angle(SPos(4))-angle(SPos(2)))第二单极子交叉环b环天线幅度校准值为:第二单极子交叉环b环天线相位校准值为:P2,b=mode(angle(SPos(6))-angle(SPos(2)))步骤2中所述通过第一单极子交叉环环天线的幅度校准值以及相位校准值对第一单极子交叉环环天线的距离多普勒谱校准为:第一单极子交叉环中a环天线的距离多普勒谱校准为:cSringa1=Sringa1/A1,a/exp(i·P1,a)第一单极子交叉环中b环天线的距离多普勒谱校准为:cSringb1=Sringb1/A1,b/exp(i·P1,b)步骤2中所述通过第二单极子交叉环环天线的幅度校准值以及相位校准值对第二单极子交叉环环天线的距离多普勒谱校准为:第二单极子交叉环中a环天线的距离多普勒谱校准为:cSringa2=Sringa2/A2,a/exp(i·P2,a)第二单极子交叉环中b环天线的距离多普勒谱校准为:cSringb2=Sringb2/A2,b/exp(i·P2,b)。4.根据权利要求1所述的基于子空间技术的海杂波抑制方法,其特征在于:步骤3中所述第一单极子交叉环与第二...

【专利技术属性】
技术研发人员:文必洋赵久瑞田震王思捷
申请(专利权)人:武汉大学
类型:发明
国别省市:湖北,42

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1