CN108872948A - 一种高频地波雷达电离层杂波抑制方法 - Google Patents
一种高频地波雷达电离层杂波抑制方法 Download PDFInfo
- Publication number
- CN108872948A CN108872948A CN201810865310.8A CN201810865310A CN108872948A CN 108872948 A CN108872948 A CN 108872948A CN 201810865310 A CN201810865310 A CN 201810865310A CN 108872948 A CN108872948 A CN 108872948A
- Authority
- CN
- China
- Prior art keywords
- clutter
- accumulation regions
- matrix
- ionospheric
- radar
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种高频地波雷达电离层杂波抑制方法,本发明涉及高频地波雷达电离层杂波抑制方法。本发明为了解决现有方法较多地用来抑制高频雷达中的窄带射频干扰以及瞬态冲击干扰,不能很好的抑制电离层的干扰,且不宜应用于宽波束雷达,以及会造成一阶Bragg峰信息的缺失以及海浪回波能量的减弱的问题。本发明相对传统杂波抑制算法,利用电离层多聚集区杂波之间的互相关性来完成对消处理,在抑制电离层杂波的同时,保留海浪回波及一阶Bragg峰的信息,减少正常海洋回波能量的下降。经过本发明方法处理后,雷达回波谱中的强电离层杂波被明显抑制,杂波抑制幅度提高,回波信杂噪比获得显著改善。本发明用于高频地波雷达探测领域。
Description
技术领域
本发明涉及高频地波雷达探测领域,具体涉及高频地波雷达电离层杂波抑制方法。
背景技术
高频地波雷达具有独特的超视距和全天候探测能力,广泛应用于舰船和飞机等运动目标的探测和海洋观测。在实际应用中,由于地波雷达发射天线做不到完全理想地使电磁波沿海面传播,因此虽然大部分能量能够沿水平面发射,但仍不可避免地存在小部分能量向上发射,并经由电离层传播反射回来,被雷达所接收形成了电离层杂波。
电离层杂波能量很强,由于电离层电子浓度和高度的时变特性,杂波通常在距离维和多普勒维均有扩展,表现出明显的非平稳信号特征,这些杂波经常使雷达在相应距离上性能急剧下降甚至失效,严重影响了雷达的探测性能。在海态反演应用中,电离层杂波能量有时甚至超过海浪回波的能量,对于海洋信息的提取十分不利,其不确定性也为杂波消除带来了极大的困难。因此,对高频地波雷达电离层杂波抑制一直是研究的难点和亟待解决的问题。
目前现有的电离层杂波抑制方法包括自适应波束形成技术、极化滤波技术等。自适应波束形成技术是通过让阵列天线方向图的主瓣指向所需方向,使波束的副瓣与外部干扰的分布相适应来抑制干扰,但它对天线主瓣波束宽度有一定的限制,因而不宜应用于宽波束雷达。极化滤波技术是利用目标和干扰在极化域的特征来进行滤波处理,由于该技术需要较多的水平天线使***变得复杂,故在高频地波雷达的实际工程应用中很少采用。常用于高频地波雷达杂波抑制的特征值分解法以及基于时频分析、基于(多普勒)频域滤波的干扰检测和抑制方法,则较多地用来抑制高频雷达中的窄带射频干扰以及瞬态冲击干扰(主要是流星余迹干扰和雷电干扰)。采用特征值分解法抑制高频地波雷达电离层杂波的同时,也会造成一阶Bragg峰信息的缺失以及海浪回波能量的减弱。
电离层一般包括D层、E层和F层。高频地波雷达所接收的电离层杂波可能存在的传播模式主要包括Es(或E)层直接回波、F层直接回波以及Es(或E)层海面模式回波等,根据其杂波分布情况可分为聚集区和扩散区。经由天波、海洋调制的混合路径产生的电离层杂波,存在于电离层高度之上且存在较大的距离扩展,由于闪烁等原因,其时频分布显得模糊散乱,有时强于海洋回波,这样的杂波区域称为扩展区。对应于天顶方向的Es层直接回波距离分布较集中且能量通常很强,有着聚集的时频分布,形成的集中分布条带状区域称为聚集区,通常出现在雷达距离-多普勒谱(R-D谱)杂波区的最底部。因Es层的强反射作用,在R-D谱的上方存在Es层杂波的一跳和二跳回波等,这些杂波区域能量一般较强且很集中,也属于聚集区。由于电离层Es层的厚度非常薄,通常情况下为百米量级,而高频地波雷达的距离分辨率为数公里,因此Es层杂波及其多跳回波在回波谱图上出现的距离相对集中,能量也主要集中在邻近的几个距离单元上(多为3至5个距离元)。聚集区杂波对同距离的有用信号形成压制性影响,必须对其进行抑制。
发明内容
本发明的目的是为了解决现有方法较多地用来抑制高频雷达中的窄带射频干扰以及瞬态冲击干扰,不能很好的抑制电离层的干扰,且不宜应用于宽波束雷达,以及会造成一阶Bragg峰信息的缺失以及海浪回波能量的减弱的问题,而提出一种高频地波雷达电离层杂波抑制方法。
一种高频地波雷达电离层杂波抑制方法包括以下步骤:
步骤一:由雷达距离-多普勒谱(R-D谱),根据电离层杂波分布特性,按距离由近至远确定对应每个电离层杂波聚集区的杂波中心距离元ki,i=0,1,…,p-1,对应Es层杂波及其一跳回波、二跳回波等,p为电离层杂波聚集区的个数;
步骤二:根据步骤一确定的杂波中心距离元ki,计算电离层杂波在杂波中心距离元ki和与其相邻的距离元ki+1上的距离相关系数
步骤三:设置电离层杂波聚集区的相关系数参考阈值Th1,比较由步骤二得到的距离相关系数的值与参考阈值Th1的大小关系,如果则判定该杂波中心距离元对应杂波区域为聚集区,数据有效;否则,判定该杂波中心距离元对应杂波区域不是聚集区,数据无效,舍弃;筛选出所有满足条件的杂波中心距离元ki,按从小到大的顺序依次排列并赋值给dc,c=0,1,…,m-1,表示聚集区的杂波中心距离元,m为有效数据个数,执行步骤四;
步骤四:利用电离层Es层杂波与其多跳回波之间的互相关性,计算电离层Es层聚集区的杂波中心距离元d0与第j跳回波聚集区的杂波中心距离元dj的互相关系数j=1,2,…,m-1;
步骤五:设置电离层杂波聚集区的互相关系数参考阈值Th2,将步骤四求得的每个互相关系数与参考阈值Th2作比较,如果则认为第j跳回波聚集区数据有效;否则,舍弃第j跳回波聚集区数据;筛选出所有满足条件的杂波中心距离元dj,对应区域为有效聚集区,执行步骤六;
步骤六:由步骤五筛选得到的每个满足条件的杂波中心距离元dj,对应的雷达接收回波谱信号记为n=1,…,N,其中N为采样点数,dj为有效聚集区的杂波中心距离元;使用的上下相邻r个距离元信号构造聚集区信号矩阵使用杂波中心距离元d0对应的信号的上下相邻r个距离元信号构造聚集区信号矩阵r的取值由聚集区内距离元的个数确定(r可以为二分之一的距离元个数);
步骤七:对应每个满足条件的dj,由步骤六得到相应的聚集区信号矩阵取和每个构造多聚集区联合参考矩阵X(和分别取第1至2r+1列);
步骤八:根据步骤七得到的多聚集区联合参考矩阵X,求得X的相关矩阵的最大似然估计R;
步骤九:对步骤八求得的最大似然估计R进行特征值分解,得到降序排列的特征值λs和与λs相对应的特征向量vs,s为特征值的个数;
步骤十:在R中,电离层杂波的能量所占比重较大,在步骤十得到的特征向量vs中,认为前L个大的特征值所对应的特征向量为电离层杂波的特征向量,取前L(L为经验值)个大的特征值所对应的特征向量为电离层杂波的特征向量,构成矩阵,形成电离层杂波子空间特征向量矩阵V;
步骤十一:利用步骤十得到的特征向量矩阵V,构造正交投影矩阵P;
步骤十二:通过将步骤七得到的多聚集区联合参考矩阵X,投影到由步骤十一确定的正交投影矩阵P的正交子空间上,来消除待处理距离元上与多聚集区联合参考矩阵相关的电离层杂波,求出待处理聚集区经过正交投影处理后的杂波抑制信号矩阵Y。
本发明的有益效果为:
本发明的目的是提供一种基于多聚集区距离相关性的高频地波雷达电离层杂波抑制方法,以解决现有背景技术中存在的问题,更好的实现电离层杂波抑制,提高高频地波雷达的探测性能。
本发明相对传统杂波抑制算法,利用电离层多聚集区杂波之间的互相关性来完成对消处理,在抑制电离层杂波的同时,保留海浪回波及一阶Bragg峰的信息,减少正常海洋回波能量的下降。
经过本发明方法处理后,雷达回波谱中的强电离层杂波被明显抑制,杂波抑制幅度提高,回波信杂噪比获得显著改善。
实施例一中定量分析本发明的抑制效果,以两个一阶Bragg峰为信号,信号以外的其余部分为杂波和噪声基底,分别计算算法处理前后的信杂噪比,计算两者差值可以得到杂波抑制幅度,它反映了算法的真实性能。取121.5km处的回波数据计算,得到处理前回波谱的信噪比为14.74dB,处理后回波谱的信噪比为30.65dB,信噪比改善幅度提高15.91dB,杂波抑制的效果明显。
实施例二中定量分析本发明的抑制效果,分别计算算法处理前后的信杂噪比。取119km处的回波数据计算,得到处理前回波谱的信噪比为13.01dB,处理后回波谱的信噪比为25.03dB,信噪比改善幅度提高12.02dB,杂波抑制的效果明显。
附图说明
图1为实施例一电离层杂波严重的原始回波R-D谱;
图2为实施例一杂波聚集区的原始回波谱;
图3为实施例一杂波对消处理后的回波R-D谱;
图4为实施例一聚集区杂波对消处理前后的回波谱对比图;
图5为实施例二电离层杂波严重的原始回波R-D谱;
图6为实施例二杂波聚集区的原始回波谱;
图7为实施例二杂波对消处理后的回波R-D谱;
图8为实施例二聚集区杂波对消处理前后的回波谱对比图;
图9为本发明流程图。
具体实施方式
具体实施方式一:一种高频地波雷达电离层杂波抑制方法包括以下步骤:
步骤一:由雷达距离-多普勒谱(R-D谱),根据电离层杂波分布特性,按距离由近至远确定对应每个电离层杂波聚集区的杂波中心距离元ki,i=0,1,…,p,对应Es层杂波及其一跳回波、二跳回波等,p为杂波中心距离元的个数;
步骤二:根据步骤一确定的杂波中心距离元ki,计算电离层杂波在杂波中心距离元ki和与杂波中心距离元相邻的距离元ki+1上的距离相关系数
步骤三:研究表明,正常海洋回波信号的距离相关系数为0.6至0.8,而聚集区内杂波的距离相关系数则超过0.88。因此可以利用聚集区内邻近距离元回波信号的强相关性,来区别正常海洋回波信号与聚集区杂波。
根据高频地波雷达回波信号的距离相关系数分布情况,设置电离层杂波聚集区的相关系数参考阈值Th1,比较由步骤二得到的距离相关系数的值与参考阈值Th1的大小关系,如果则判定该距离元对应区域为聚集区,数据有效;否则,判定该距离元对应区域不是聚集区,数据无效,舍弃;筛选出所有满足条件的杂波中心距离元ki,按从小到大的顺序依次排列并赋值给dc,c=0,1,…,m-1,表示聚集区的杂波中心距离元,m为有效数据个数,执行步骤四;
步骤四:电离层杂波不仅在每个聚集区内有很强的相关性,在Es层杂波与其多跳回波形成的多个聚集区之间也存在着很强的互相关性。由实测数据的相关性计算结果可知,Es层杂波与其多跳回波之间的互相关系数为0.65至0.8,而正常海洋回波信号与Es层杂波及多跳回波的互相关系数小于0.45,即前者的互相关性大于后者的互相关性。因此可以利用Es层杂波与其多跳回波形成的多个聚集区之间的互相关性,来实现对多个聚集区中包含强电离层干扰雷达回波信号的杂波抑制。
利用电离层Es层杂波与其多跳回波之间的互相关性,计算电离层Es层聚集区的杂波中心距离元d0与第j跳回波聚集区的杂波中心距离元dj的互相关系数j=1,2,…,m-1;
步骤五:根据Es层杂波与其多跳回波形成的多个聚集区之间的互相关系数分布情况,设置电离层杂波聚集区的互相关系数参考阈值Th2,将步骤四求得的每个互相关系数与参考阈值Th2作比较,如果则认为第j跳回波聚集区数据有效;否则,舍弃第j跳回波聚集区数据;筛选出所有满足条件的杂波中心距离元dj,对应区域为有效聚集区,执行步骤六;
步骤六:由步骤五筛选得到的每个满足条件的杂波中心距离元dj,对应的雷达接收回波谱信号记为n=1,…,N,其中N为采样点数,dj为有效聚集区的杂波中心距离元;使用的上下相邻r个距离元信号构造聚集区信号矩阵使用杂波中心距离元d0对应的信号的上下相邻r个距离元信号构造聚集区信号矩阵r的取值由聚集区内距离元的个数确定(r可以为二分之一的距离元个数);
步骤七:对应每个满足条件的dj,由步骤六得到相应的聚集区信号矩阵取和每个构造多聚集区联合参考矩阵X(和分别取第1至2r+1列);
步骤八:根据步骤七得到的多聚集区联合参考矩阵X,求得X的相关矩阵的最大似然估计R;
步骤九:对步骤八求得的最大似然估计R进行特征值分解,得到降序排列的特征值λs和与λs相对应的特征向量vs,s为特征值的个数;
步骤十:在R中,电离层杂波的能量所占比重较大,在步骤十得到的特征向量vs中,认为前L个大的特征值所对应的特征向量为电离层杂波的特征向量,取前L(L为经验值)个大的特征值所对应的特征向量为电离层杂波的特征向量,构成矩阵,形成电离层杂波子空间特征向量矩阵V;
步骤十一:利用步骤十得到的特征向量矩阵V,构造正交投影矩阵P;
步骤十二:通过将步骤七得到的多聚集区联合参考矩阵X,投影到由步骤十一确定的正交投影矩阵P的正交子空间上,来消除待处理距离元上与多聚集区联合参考矩阵相关的电离层杂波,求出待处理聚集区经过正交投影处理后的杂波抑制信号矩阵Y。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤二中计算电离层杂波在杂波中心距离元ki和与杂波中心距离元相邻的距离元ki+1上的距离相关系数的具体过程为:
其中X(ki)为雷达在第ki个距离元上的回波数据,Cov[X(ki),X(ki+1)]为第ki个距离元与第ki+1个距离元上回波数据的协方差,Var[X(ki)]为X(ki)的方差; 越大,电离层杂波在两个相邻距离元ki与ki+1上的相关性越强。
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤四中计算电离层聚集区的杂波中心距离元d0与第j跳回波聚集区的杂波中心距离元dj的互相关系数的具体过程为:
其中X(d0)为雷达在Es层聚集区的杂波中心距离元d0上的回波数据,X(dj)为雷达在第j跳回波聚集区的杂波中心距离元dj上的回波数据,Cov[X(d0),X(dj)]为电离层Es层杂波与其第j跳回波数据的协方差,Var[X(dj)]为X(dj)的方差; 越大,电离层Es层杂波与其第j跳回波之间的互相关性越强。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤六中使用每个有效距离元信号的上下相邻r个距离元信号构造聚集区信号矩阵具体为:
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述步骤七中取和每个构造多聚集区联合参考矩阵X具体为:
其它步骤及参数与具体实施方式一至四之一相同。
具体实施方式六:本实施方式与具体实施方式一至五之一不同的是:所述步骤八中根据步骤七得到的多聚集区联合参考矩阵X,求得X的相关矩阵的最大似然估计R具体为:
其中XH为X的共轭转置。
其它步骤及参数与具体实施方式一至五之一相同。
具体实施方式七:本实施方式与具体实施方式一至六之一不同的是:所述步骤九中对步骤八求得的最大似然估计R进行特征值分解,得到降序排列的特征值λs和与λs相对应的特征向量vs的具体过程为:
其中为vs的共轭转置。
其它步骤及参数与具体实施方式一至六之一相同。
具体实施方式八:本实施方式与具体实施方式一至七之一不同的是:所述步骤十中取前L个大的特征值所对应的特征向量为电离层杂波的特征向量,构成矩阵,形成电离层杂波子空间特征向量矩阵V具体为:
V=[v1,…,vL]
其它步骤及参数与具体实施方式一至七之一相同。
具体实施方式九:本实施方式与具体实施方式一至八之一不同的是:所述步骤十一中利用步骤十得到的特征向量矩阵V,构造正交投影矩阵P的具体为:
P=I-VVH
其中VH为V的共轭转置。
其它步骤及参数与具体实施方式一至八之一相同。
具体实施方式十:本实施方式与具体实施方式一至九之一不同的是:所述步骤十二中求出待处理聚集区经过正交投影处理后的杂波抑制信号矩阵Y的具体过程为:
Y=XPH
其中PH为P的共轭转置;
得到的矩阵Y是杂波抑制后各聚集区回波的联合矩阵,即:
其中为杂波中心距离元d0对应聚集区经过杂波抑制后的回波信号矩阵,为第j跳回波聚集区经过杂波抑制后的回波信号矩阵。
其它步骤及参数与具体实施方式一至九之一相同。
实施例一:
实验中采用的数据来源于哈尔滨工业大学电子所某雷达实验站于2016年4月10日下午18时采集的实测数据,表1给出雷达的相关参数。
表1雷达***参数
根据相关要求,部分参数不予给出,以*代替。
原始雷达回波数据经两次FFT运算、一次DBF运算后得到回波数据的方位-距离-多普勒谱(ARD谱),由ARD谱可得到每个波束的距离-多普勒谱(R-D谱)和每个距离元的多普勒谱。当高频地波雷达受到很强的电离层杂波干扰时,其第4波束上接收回波的R-D谱和聚集区的回波谱如图1和图2所示。图1中120km,240km和350km附近呈横向条带状分布的杂波分别为Es层杂波、一跳回波和二跳回波,相应区域为杂波集中分布的聚集区,图2选取杂波聚集区内一个距离门上的数据,在多普勒频率-0.2Hz和0.2Hz附近相对零频对称的两个尖峰是一阶Bragg峰,可以看出聚集区内的电离层杂波很强,对需要提取的一阶Bragg峰信息造成了很大影响。
采用基于多聚集区距离相关性的高频地波雷达电离层杂波对消算法,取r=4,L=8,
图3给出了利用多聚集区距离相关性对消处理后得到的回波R-D谱,从图中可以看出,经过算法处理后,聚集区的电离层杂波被明显抑制,同时保留了一阶Bragg峰的信息,正常海洋回波能量也没有明显减弱。
图4给出了利用多聚集区距离相关性对消处理前后得到的回波谱对比图,从图中可以看出,经过算法处理后,两个一阶Bragg峰之间及Bragg峰附近的强电离层杂波被明显抑制。
定量分析本发明的抑制效果,以两个一阶Bragg峰为信号,信号以外的其余部分为杂波和噪声基底,分别计算算法处理前后的信杂噪比,计算两者差值可以得到杂波抑制幅度,它反映了算法的真实性能。取121.5km处的回波数据计算,得到处理前回波谱的信噪比为14.74dB,处理后回波谱的信噪比为30.65dB,信噪比改善幅度提高15.91dB,杂波抑制的效果明显。
实施例二:
实验中采用的数据来源于哈尔滨工业大学电子所某雷达实验站于2016年4月10日下午19时采集的实测数据,表2给出雷达的相关参数。
表2雷达***参数
根据相关要求,部分参数不予给出,以*代替。
原始雷达回波数据在第4个波束上接收回波的R-D谱和聚集区的回波谱如图5和图6所示。图5中120km,240km和350km附近为杂波集中分布的聚集区。图6选取杂波聚集区内一个距离门上的数据,在多普勒频率-0.2Hz和0.2Hz附近相对零频对称的两个尖峰是一阶Bragg峰。
图7给出了利用多聚集区距离相关性对消处理后得到的回波R-D谱,从图中可以看出,经过算法处理后,聚集区的电离层杂波被明显抑制,同时保留了一阶Bragg峰的信息,正常海洋回波能量也没有明显减弱。
图8给出了利用多聚集区距离相关性对消处理前后得到的回波谱对比图,从图中可以看出,杂波对消处理前,正一阶Bragg峰被电离层杂波淹没,算法处理后,两个一阶Bragg峰之间及Bragg峰附近的强电离层杂波被明显抑制,被淹没的正一阶Bragg峰显现出来。
定量分析本发明的抑制效果,分别计算算法处理前后的信杂噪比。取119km处的回波数据计算,得到处理前回波谱的信噪比为13.01dB,处理后回波谱的信噪比为25.03dB,信噪比改善幅度提高12.02dB,杂波抑制的效果明显。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种高频地波雷达电离层杂波抑制方法,其特征在于:所述高频地波雷达电离层杂波抑制方法包括以下步骤:
步骤一:由雷达距离-多普勒谱,确定对应每个电离层杂波聚集区的杂波中心距离元ki,i=0,1,…,p-1,p为电离层杂波聚集区的个数;
步骤二:根据步骤一确定的杂波中心距离元ki,计算电离层杂波在杂波中心距离元ki和与其相邻的距离元ki+1上的距离相关系数
步骤三:设置电离层杂波聚集区的相关系数参考阈值Th1,比较由步骤二得到的距离相关系数的值与参考阈值Th1的大小关系,如果则判定该杂波中心距离元对应杂波区域为聚集区,数据有效;否则,判定该杂波中心距离元对应杂波区域不是聚集区,数据无效,舍弃;筛选出所有满足条件的杂波中心距离元ki,按从小到大的顺序依次排列并赋值给dc,c=0,1,…,m-1,表示聚集区的杂波中心距离元,m为有效数据个数,执行步骤四;
步骤四:计算杂波中心距离元d0与第j跳回波聚集区的杂波中心距离元dj的互相关系数
步骤五:设置电离层杂波聚集区的互相关系数参考阈值Th2,将步骤四求得的每个互相关系数与参考阈值Th2作比较,如果则认为第j跳回波聚集区数据有效;否则,舍弃第j跳回波聚集区数据;筛选出所有满足条件的杂波中心距离元dj,对应区域为有效聚集区,执行步骤六;
步骤六:由步骤五筛选得到的每个满足条件的杂波中心距离元dj,对应的雷达接收回波谱信号记为其中N为采样点数,dj为有效聚集区的杂波中心距离元;使用的上下相邻r个距离元信号构造聚集区信号矩阵使用杂波中心距离元d0对应的雷达接收回波谱信号的上下相邻r个距离元信号构造聚集区信号矩阵
步骤七:取和构造多聚集区联合参考矩阵X;
步骤八:根据步骤七得到的多聚集区联合参考矩阵X,求得X的相关矩阵的最大似然估计R;
步骤九:对步骤八求得的最大似然估计R进行特征值分解,得到降序排列的特征值λs和与λs相对应的特征向量vs,s为特征值的个数;
步骤十:取前L个大的特征值所对应的特征向量为电离层杂波的特征向量,构成矩阵,形成电离层杂波子空间特征向量矩阵V;
步骤十一:利用步骤十得到的特征向量矩阵V,构造正交投影矩阵P;
步骤十二:通过将步骤七得到的多聚集区联合参考矩阵X,投影到由步骤十一确定的正交投影矩阵P的正交子空间上,来消除待处理距离元上与多聚集区联合参考矩阵相关的电离层杂波,求出待处理聚集区经过正交投影处理后的杂波抑制信号矩阵Y。
2.根据权利要求1所述一种高频地波雷达电离层杂波抑制方法,其特征在于:所述步骤二中计算电离层杂波在杂波中心距离元ki和与其相邻的距离元ki+1上的距离相关系数的具体过程为:
其中X(ki)为雷达在第ki个距离元上的回波数据,Cov[X(ki),X(ki+1)]为第ki个距离元与第ki+1个距离元上回波数据的协方差,Var[X(ki)]为X(ki)的方差; 越大,电离层杂波在两个相邻距离元ki与ki+1上的相关性越强。
3.根据权利要求2所述一种高频地波雷达电离层杂波抑制方法,其特征在于:所述步骤四中计算杂波中心距离元d0与第j跳回波聚集区的杂波中心距离元dj的互相关系数的具体过程为:
其中X(d0)为雷达在Es层聚集区的杂波中心距离元d0上的回波数据,X(dj)为雷达在第j跳回波聚集区的杂波中心距离元dj上的回波数据,Cov[X(d0),X(dj)]为电离层杂波与其第j跳回波数据的协方差,Var[X(dj)]为X(dj)的方差; 越大,电离层杂波与其第j跳回波之间的互相关性越强。
4.根据权利要求3所述一种高频地波雷达电离层杂波抑制方法,其特征在于:所述步骤六中使用每个有效距离元信号的上下相邻r个距离元信号构造聚集区信号矩阵具体为:
5.根据权利要求4所述一种高频地波雷达电离层杂波抑制方法,其特征在于:所述步骤七中取和构造多聚集区联合参考矩阵X具体为:
6.根据权利要求5所述一种高频地波雷达电离层杂波抑制方法,其特征在于:所述步骤八中根据步骤七得到的多聚集区联合参考矩阵X,求得X的相关矩阵的最大似然估计R具体为:
其中XH为X的共轭转置。
7.根据权利要求6所述一种高频地波雷达电离层杂波抑制方法,其特征在于:所述步骤九中对步骤八求得的最大似然估计R进行特征值分解,得到降序排列的特征值λs和与λs相对应的特征向量vs的具体过程为:
其中为vs的共轭转置。
8.根据权利要求7所述一种高频地波雷达电离层杂波抑制方法,其特征在于:所述步骤十中取前L个大的特征值所对应的特征向量为电离层杂波的特征向量,构成矩阵,形成电离层杂波子空间特征向量矩阵V具体为:
V=[v1,…,vL]。
9.根据权利要求8所述一种高频地波雷达电离层杂波抑制方法,其特征在于:所述步骤十一中利用步骤十得到的特征向量矩阵V,构造正交投影矩阵P的具体为:
P=I-VVH
其中VH为V的共轭转置。
10.根据权利要求9所述一种高频地波雷达电离层杂波抑制方法,其特征在于:所述步骤十二中求出待处理聚集区经过正交投影处理后的杂波抑制信号矩阵Y的具体过程为:
Y=XPH
其中PH为P的共轭转置;
得到的矩阵Y是杂波抑制后各聚集区回波的联合矩阵,即:
其中为杂波中心距离元d0对应聚集区经过杂波抑制后的回波信号矩阵,为第j跳回波聚集区经过杂波抑制后的回波信号矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810865310.8A CN108872948B (zh) | 2018-08-01 | 2018-08-01 | 一种高频地波雷达电离层杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810865310.8A CN108872948B (zh) | 2018-08-01 | 2018-08-01 | 一种高频地波雷达电离层杂波抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108872948A true CN108872948A (zh) | 2018-11-23 |
CN108872948B CN108872948B (zh) | 2022-04-12 |
Family
ID=64307376
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810865310.8A Active CN108872948B (zh) | 2018-08-01 | 2018-08-01 | 一种高频地波雷达电离层杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108872948B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109444840A (zh) * | 2018-12-04 | 2019-03-08 | 南京航空航天大学 | 一种基于机器学习的雷达杂波抑制方法 |
CN109597034A (zh) * | 2018-12-12 | 2019-04-09 | 哈尔滨工业大学 | 一种基于欧几里得距离的空时自适应处理方法 |
CN109655793A (zh) * | 2018-12-07 | 2019-04-19 | 南京邮电大学 | 高频雷达从瞬态干扰中提取回波的方法 |
CN109765529A (zh) * | 2018-12-30 | 2019-05-17 | 成都汇蓉国科微***技术有限公司 | 一种基于数字波束形成的毫米波雷达抗干扰方法及*** |
CN110045337A (zh) * | 2019-05-10 | 2019-07-23 | 武汉大学 | 基于张量子空间投影的高频地波雷达射频干扰抑制方法 |
CN111175709A (zh) * | 2019-12-27 | 2020-05-19 | 中国电子科技集团公司第十四研究所 | 一种基于误差抑制的面向大范围气象雷达的拼图方法 |
CN112363122A (zh) * | 2020-11-23 | 2021-02-12 | 哈尔滨工业大学(威海) | 高频地波雷达电离层噪声中弱谐波信号的提取方法及应用 |
CN112882018A (zh) * | 2021-01-13 | 2021-06-01 | 哈尔滨工业大学(威海) | 一种海洋和电离层一体化探测高频雷达***及其控制方法 |
CN114660537A (zh) * | 2022-04-15 | 2022-06-24 | 武汉大学 | 基于Es层的电离层F层天波俯仰角估计方法及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226236A (zh) * | 2007-12-18 | 2008-07-23 | 西安电子科技大学 | 基于旁瓣约束的天波超视距雷达自适应干扰抑制方法 |
CN104133198A (zh) * | 2014-08-13 | 2014-11-05 | 武汉大学 | 一种高频地波雷达中电离层干扰抑制方法 |
CN104155632A (zh) * | 2014-07-18 | 2014-11-19 | 南京航空航天大学 | 一种基于局部相关性的改进子空间海杂波抑制方法 |
-
2018
- 2018-08-01 CN CN201810865310.8A patent/CN108872948B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226236A (zh) * | 2007-12-18 | 2008-07-23 | 西安电子科技大学 | 基于旁瓣约束的天波超视距雷达自适应干扰抑制方法 |
CN104155632A (zh) * | 2014-07-18 | 2014-11-19 | 南京航空航天大学 | 一种基于局部相关性的改进子空间海杂波抑制方法 |
CN104133198A (zh) * | 2014-08-13 | 2014-11-05 | 武汉大学 | 一种高频地波雷达中电离层干扰抑制方法 |
Non-Patent Citations (4)
Title |
---|
JIANG WEI 等: "Characteristic study of ionospheric clutter in high-frequency over the horizon surface wave radar", 《2009 IEEE YOUTH CONFERENCE ON INFORMATION, COMPUTING AND TELECOMMUNICATION》 * |
XU XING-AN 等: "Ionospheric clutter analysis and suppression in a distributed MIMO-HFSWR system", 《OCEANS 2014 - TAIPEI》 * |
吴敏: "高频地波雷达干扰与杂波抑制方法的研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
周浩 等: "高频地波雷达中电离层Es层杂波分析及其抑制", 《华中科技大学学报(自然科学版)》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109444840B (zh) * | 2018-12-04 | 2020-12-25 | 南京航空航天大学 | 一种基于机器学习的雷达杂波抑制方法 |
CN109444840A (zh) * | 2018-12-04 | 2019-03-08 | 南京航空航天大学 | 一种基于机器学习的雷达杂波抑制方法 |
CN109655793A (zh) * | 2018-12-07 | 2019-04-19 | 南京邮电大学 | 高频雷达从瞬态干扰中提取回波的方法 |
CN109597034A (zh) * | 2018-12-12 | 2019-04-09 | 哈尔滨工业大学 | 一种基于欧几里得距离的空时自适应处理方法 |
CN109597034B (zh) * | 2018-12-12 | 2021-08-31 | 哈尔滨工业大学 | 一种基于欧几里得距离的空时自适应处理方法 |
CN109765529A (zh) * | 2018-12-30 | 2019-05-17 | 成都汇蓉国科微***技术有限公司 | 一种基于数字波束形成的毫米波雷达抗干扰方法及*** |
CN109765529B (zh) * | 2018-12-30 | 2020-11-10 | 成都汇蓉国科微***技术有限公司 | 一种基于数字波束形成的毫米波雷达抗干扰方法及*** |
CN110045337A (zh) * | 2019-05-10 | 2019-07-23 | 武汉大学 | 基于张量子空间投影的高频地波雷达射频干扰抑制方法 |
CN110045337B (zh) * | 2019-05-10 | 2023-07-21 | 武汉大学 | 基于张量子空间投影的高频地波雷达射频干扰抑制方法 |
CN111175709B (zh) * | 2019-12-27 | 2023-02-24 | 中国电子科技集团公司第十四研究所 | 一种基于误差抑制的面向大范围气象雷达的拼图方法 |
CN111175709A (zh) * | 2019-12-27 | 2020-05-19 | 中国电子科技集团公司第十四研究所 | 一种基于误差抑制的面向大范围气象雷达的拼图方法 |
CN112363122A (zh) * | 2020-11-23 | 2021-02-12 | 哈尔滨工业大学(威海) | 高频地波雷达电离层噪声中弱谐波信号的提取方法及应用 |
CN112882018B (zh) * | 2021-01-13 | 2022-07-26 | 哈尔滨工业大学(威海) | 一种海洋和电离层一体化探测高频雷达***及其控制方法 |
CN112882018A (zh) * | 2021-01-13 | 2021-06-01 | 哈尔滨工业大学(威海) | 一种海洋和电离层一体化探测高频雷达***及其控制方法 |
CN114660537B (zh) * | 2022-04-15 | 2022-11-15 | 武汉大学 | 基于Es层的电离层F层天波俯仰角估计方法及设备 |
CN114660537A (zh) * | 2022-04-15 | 2022-06-24 | 武汉大学 | 基于Es层的电离层F层天波俯仰角估计方法及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN108872948B (zh) | 2022-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108872948A (zh) | 一种高频地波雷达电离层杂波抑制方法 | |
US7145503B2 (en) | Surface wave radar | |
CN101226236B (zh) | 基于旁瓣约束的天波超视距雷达自适应干扰抑制方法 | |
DE60304692T2 (de) | System und verfahren zum erzeugen eines spektrums in radargeräten | |
CN104155648B (zh) | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 | |
CN106707247B (zh) | 一种基于紧凑天线阵的高频海洋雷达目标检测方法 | |
CN107561512B (zh) | 一种脉冲多普勒雷达抗压制式拖曳干扰的极化对消方法 | |
CN101881822B (zh) | 一种针对共享谱雷达同频干扰的抑制方法 | |
CN104155632A (zh) | 一种基于局部相关性的改进子空间海杂波抑制方法 | |
CN108931766A (zh) | 一种基于稀疏重构的非均匀stap干扰目标滤除方法 | |
CN112379333B (zh) | 基于空时维正交投影滤波的高频雷达海杂波抑制方法 | |
Droszcz et al. | Beamforming of LOFAR radio telescope antennas used as sensors in passive radiolocation system | |
Bosse et al. | Model-based multifrequency array signal processing for low-angle tracking | |
Zhang et al. | Improved main-lobe cancellation method for space spread clutter suppression in HFSSWR | |
Wen-long et al. | Ionospheric clutter mitigation for high-frequency surface-wave radar using two-dimensional array and beam space processing | |
Davis et al. | A maximum-likelihood beamspace processor for improved search and track | |
McDonald et al. | Elevation angle-of-arrival determination for a standard and a modified SuperDARN HF radar layout | |
Bosse et al. | Improved radar tracking using a multipath model: maximum likelihood compared with eigenvector analysis | |
CN114814738A (zh) | 高频雷达射频干扰抑制的自适应参考距离单元选择方法 | |
Zhao et al. | Using sky-wave echoes information to extend HFSWR's maximum detection range | |
CN109061597A (zh) | 基于盲源分离与时频脊波域滤波的电离层杂波抑制方法 | |
Wan et al. | Experimental investigation of directional characteristics for ionospheric clutter in HF surface wave radar | |
CN108474845A (zh) | 用于对含盐液体介质中的杂波进行谱估计的方法 | |
Zhang et al. | Ionospheric clutter extraction in HFSWR based on range-doppler spectral image processing | |
Parent | A frequency averaging method to improve sea-state measurements with a HF skywave radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |