CN107064904A - 一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法 - Google Patents
一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法 Download PDFInfo
- Publication number
- CN107064904A CN107064904A CN201710428463.1A CN201710428463A CN107064904A CN 107064904 A CN107064904 A CN 107064904A CN 201710428463 A CN201710428463 A CN 201710428463A CN 107064904 A CN107064904 A CN 107064904A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtr
- mtd
- mtable
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/411—Identification of targets based on measurements of radar reflectivity
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法,本发明涉及舰载高频地波雷达目标方位高分辨方法。本发明解决了传统合成孔径方法不能应用在海态环境、传统舰载高频地波雷达信号处理方法不能在方位向区分邻近目标的问题。本发明先使用正交投影算法对回波信号在多普勒域进行处理,检测出目标以及目标所在的多普勒频率;再将阵列划分为三个子阵列,利用子阵列在对目标所在的多普勒频率的信号使用正交投影算法进行海杂波抑制,并通过求和将子阵列信号等效成单个阵元信号;最后使用重叠相关算法对等效的三个等效阵元的接收信号做虚拟孔径处理和波束形成,得到目标方位的高分辨结果。本发明应用于雷达目标方位估计和分辨领域。
Description
技术领域
本发明属于雷达目标探测领域,具体涉及舰载高频地波雷达目标方位高分辨方法。
背景技术
舰载高频地波雷达(High Frequency Surface Wave Radar,HFSWR)是以载舰为平台,使用高频波段的电磁波对海面目标或掠海飞行的低空目标进行探测的一种新型雷达。它工作在3-30MHz的频带范围,由垂直极化的天线来发射高频波段的电磁波,利用垂直极化电磁波沿海面传播衰减小来实现大范围探测;利用高频波段电磁波的曲面绕射效应,能突破地球曲率的限制,达到超视距的效果。不过受到舰船自身长度的约束,舰载高频地波雷达的天线阵列的孔径较小,使得雷达在探测时形成的波束宽度很宽,方位分辨能力极低;同时,受舰船平台移动特性的影响,舰载高频地波雷达的一阶Bragg峰出现展宽,使得运动目标在多普勒域被淹没在展宽后的一阶海杂波中,影响目标的检测与估计。
如何基于有限长度的平台,在强大的一阶海杂波背景下实现目标的方位高分辨,是舰载高频地波雷达研究和研制中一个需要解决的问题,目前,有一些应用于其他体制雷达的方位分辨方法,如超分辨算法、合成孔径方法等。
(一)超分辨方法。这类方法是依据优化约束准则或者子空间分解,来区分信号、噪声以及干扰,达到方位超分辨的目的。超分辨算法的提出最早可以追溯到六十年代,由BurgJ.P.提出的最大熵谱估计算法(MEM),以及由Capon J.提出的最小方差谱估计算法(MVM)。而在八十年代由Schmit R.O.等人提出的多重信号分类(MUSIC)算法以及随后由Paulraj等人提出的ESPRIT算法最为突出。这几种算法在理想情况下的估计性能能逼近Cramer-Rao界,前提是必须准确地知道阵列流型以及背景噪声模型。
(二)合成孔径方法。这类方法是利用雷达平台和被测目标间的相对运动,在一定的积累时间内,将雷达在不同空间位置上接收的回波信号进行相干处理,合成比雷达真实天线孔径大得多的合成阵列,实现方位向的高分辨。合成孔径雷达(Synthetic ApertureRadar,SAR)的概念在二十世纪五十年代就已经被提出了,1951年,Carl.Wiley最早提出用频率分析的方法可以提高目标方位分辨率;1971年,美国学者Raney首次提出动目标检测技术(Ground Moving Target Indication,GMTI)与SAR***二者结合,实现GMTI和合成孔径一体化,可以在实现对静止目标高分辨的基础上完成运动目标检测和参数估计。随着技术发展,多通道SAR处理方法逐渐出现,1991年Lightstone L等人在Multiple phase centreDPCA for airborne radar[C].Proceedings of IEEE National Radar Conference,1991:36-40.中提出对空雷达的相位中心偏置(Displaced Phase Center Antenna,DPCA)技术,Kim D等人在Remote sensing of ocean waves and currents using NASA(JPL)AIRSAR along-track interferometry(ATI)[C].IGARSS,2002,2:931-933.中提出沿航迹干涉(Along-Track Interferometry,ATI)处理方法,Barbarossa S等人在Space-time-frequency processing of synthetic aperture radar signals[J].IEEE Transactionson Aerospace and Electronic Systems,1994,30(2):341-358.中提出空时二维自适应处理处理(sky-time adaptive processing,STAP)技术。
仿照雷达领域的机载与星载合成孔径,在声呐探测领域,提出了被动合成孔径的方法。1976年,Williams R E在Williams R E.Creating an acoustic syntheticaperture in the ocean[J].Journal of the Acoustical Society of America,1976,60(1):60-73.中提出了声呐领域的被动合成孔径方法;1989年由Stergiopoulos S和Sullivan E J等人在Stergiopoulos S,Sullivan E J.Extended towed arrayprocessing by an overlap correlator[J].Journal of the Acoustical Society ofAmerica,1989,86(86):158-171.中提出扩展拖曳阵列尺寸(Extended Towed ArrayMeasurements,ETAM)算法;在近年,为了改善ETAM算法的性能,出现了不少针对ETAM的改进算法,如应用线性预测技术的被动合成孔径处理方法、基于非重叠阵元的NOEA算法等。
通常定义雷达实阵列的方位分辨率为式中λ为雷达信号的波长,D为雷达实阵列的孔径大小,高分辨方法是指方位分辨率高于实阵列的分辨方法。
发明内容
本发明的目的是为了解决现有技术提高舰载高频地波雷达在进行海面目标探测时,受实际阵列孔径的限制,在进行多个海面目标的方位分辨时方位分辨力低,不能很好地区分出邻近目标的缺点,以及由于一阶海杂波影响,使得传统高分辨方法不能直接应用于舰载条件下的缺点,而提出一种载高频地波雷达海面目标方位高分辨方法。
一种载高频地波雷达海面目标方位高分辨方法包括以下步骤:
步骤一:设舰载高频地波雷达的天线阵列为线阵,阵元数为M,M为大于3的正整数,雷达的天线阵列接收到的信号为X(q),对矩阵X(q)中的每一个脉冲周期内的q1、q2、q3、……、qc采样时刻的信号分别作FFT,得到每个采样时刻的M个阵元的回波的多普勒谱Fq1(s),Fq2(s),…,Fqc(s);所述FFT为快速傅里叶变换,q为采样时刻序列,q=[q1,q2,…,qc],其中q1、q2…qc为采样时刻,c为每一个脉冲周期内时域采样总数;
步骤二:对步骤一得到的q1采样时刻(可以取q1、q2、q3、……、qc采样时刻中的任一时刻进行步骤二至步骤三的过程)M个阵元的回波的多普勒谱Fq1(s),使用正交投影算法在频域进行一阶海杂波抑制,得到抑制后的回波信号多普勒矩阵FDPRS(s);
步骤三:对FDPRS(s)矩阵按列求和,得到求和后的序列检测序列中的峰值得到l个峰值在序列中的序号Ipeak(1)、Ipeak(2)、…、Ipeak(l)以及峰值对应的频率speak(1)、speak(2)、…、speak(l);
步骤四:选取正整数L作为天线阵列的子阵阵元个数,将步骤一中得到的q1采样时刻M个阵元的回波的多普勒矩阵(多普勒谱)Fq1(s)按照第1行到第L行、第2r+1行到第(2r+L)行、第4r+1行到第(4r+L)行、第6r+1行到第(6r+L)行……的规律重新写成[(M-L)/2r+1]([(M-L)/2r+1]也可以表示为)个部分,由[(M-L)/2r+1]个部分组成天线阵列所有子阵的回波的多普勒矩阵Fsub(s);令多普勒矩阵Fsub(s)中多普勒频率为speak(1)的列上的数据保持不变,其余位置的值为0,对Fsub(s)中的[(M-L)/2r+1]个部分别使用正交投影算法进行一阶海杂波抑制,得到抑制后的天线阵列所有子阵的回波的多普勒矩阵其中L≤2r,r<M/2,r为天线阵列在相邻两个脉冲周期之间走过的阵元个数;[(M-L)/2r+1]中的[]表示向下取整;
步骤五:对步骤四中得到的抑制后的天线阵列所有子阵的回波的多普勒矩阵的[(M-L)/2r+1]个部分的每一部分分别按列求和,得到一个[(M-L)/2r+1]×R大小的中间变量矩阵对的每一行做逆傅里叶变换(IFFT)得到q1采样时刻的等效阵元的接收信号xΣq1(n),其中n表示第n个脉冲周期,1≤n≤N,N为积累脉冲数;
步骤六:依次取每个脉冲周期内q2、q3、……、qc采样时刻的M个阵元的回波的多普勒矩阵Fq2(s),Fq3(s),…,Fqc(s),对speak(1),重复执行步骤四和步骤五,得到q2、q3、……、qc采样时刻的等效阵元的接收信号xΣq2(n)、xΣq3(n)、……、xΣqc(n),根据xΣq1(n)、xΣq2(n)、……、xΣqc(n)得到海杂波抑制后的等效阵列的接收信号XΣ(n);
步骤七:使用重叠相关算法对信号矩阵XΣ(n)进行处理,得到多普勒频率为speak(1)的目标的虚拟孔径波束形成结果y1(θ),其中θ表示角度,角标1表示第一个目标;
步骤八:对步骤三中得到的频率为speak(2)、……、speak(l),重复执行步骤四至步骤七,直至得到第2、3、……、l个目标的虚拟孔径波束形成结果y2(θ)、y3(θ)、……、yl(θ),完成对舰载高频地波雷达的l个目标的方位高分辨。
本发明的有益效果为:
本发明一种舰载高频地波雷达方位高分辨方法,为舰载高频地波雷达对目标进行高分辨以及定位提供了一种有效方法。相较于被动合成孔径算法,本发明可以在存在一阶海杂波的情况下实现目标的方位高分辨;相较于传统的舰载高频地波雷达多普勒域处理方法,本发明可以在分辨目标的同时对目标实现方位估计。
本发明的目的是为了解决现有舰载高频地波雷达信号处理方法不能实现邻近目标在方位向上的高分辨,以及传统合成孔径方法不能在一阶海杂波背景下适用的缺点。
本发明通过使用正交投影加权法抑制了一阶海杂波,再通过将四阵元子阵杂波抑制后求和等效成一个等效阵元,解决了传统合成孔径方法不能在一阶海杂波背景下适用的缺点。正交投影加权法是依据目标回波矢量和一阶海杂波矢量在向量空间中方向不重叠,利用正交加权将回波投影到与一阶海杂波正交的方向上去,目标回波投影后不为0,一阶海杂波投影后变成了0,达到了抑制一阶海杂波的效果;相较于传统合成孔径方位高分辨方法,本发明处理后的信号的信杂比提升了10-20dB。
通过一阶海杂波抑制后的多普勒域分选以及重叠相关算法,解决了舰载条件下邻近目标在方位向的定位和高分辨。重叠相关算法实质是利用阵列平台在空间中的移动,通过重叠相关因子补偿将空时分布的阵列对齐到同一时刻,等效成一个在空间分布的大孔径虚拟阵列,用这个虚拟阵列进行波束形成就可以实现邻近目标在方位向的定位和高分辨;相较于传统的实阵列波束形成方法,本发明的方法在分辨能力上比传统方法分辨能力高出10-25倍。传统的实阵列波束形成方法对40度和50度上两个目标的方位分辨结果如图8所示,本发明方法对40度和50度上两个目标的方位分辨结果如图9所示。
附图说明
图1是回波信号的多谱勒谱f1(s)的示意图;
图2是使用正交投影算法进行海杂波抑制后的回波信号的角度-多谱勒矩阵FDPRS(s)示意图;
图3是步骤四中求和后的序列的曲线示意图;
图4是步骤五中多谱勒分选并进行海杂办抑制后的的曲线示意图;
图5是对的每一行做逆傅里叶变换(IFFT),得到xΣq1(n)的曲线示意图;
图6是使用本发明方法,对多普勒频率为-0.205Hz的这一个目标处理后的方位分辨结果y1(θ)的示意图;
图7是使用本发明方法,对多普勒频率为-0.155Hz的这一个目标处理后的方位分辨结果y2(θ)的示意图;
图8是传统的实阵列波束形成方法对40度和50度上两个目标的方位分辨结果图;
图9是本发明方法对40度和50度上两个目标的方位分辨结果图;
具体实施方式
具体实施方式一:一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法包括以下步骤:
步骤一:设舰载高频地波雷达的天线阵列为线阵,阵元数为M,M为大于3的正整数,雷达的天线阵列接收到的信号为X(q),对矩阵X(q)中的每一个脉冲周期内的q1、q2、q3、……、qc采样时刻的信号分别作FFT,得到每个采样时刻的M个阵元的回波的多普勒谱Fq1(s),Fq2(s),…,Fqc(s);所述FFT为快速傅里叶变换,q为采样时刻序列,q=[q1,q2,...,qc],其中q1、q2…qc为采样时刻,c为每一个脉冲周期内时域采样总数;阵元之间的间距为d,雷达信号的波长为λ,脉冲宽度为τ,d需要满足d≤λ/2,阵列平台速度为v,且存在正整数r、u(r<M/2),使得abs(vuT1-rd)的值最小,其中abs表示取绝对值,T1为脉冲重复周期,令Tp=uT1,积累脉冲数为N,采样频率为fs,q=[q1,q2,...,qc]∈[0,fsτ];
步骤二:对步骤一得到的q1采样时刻(可以取q1、q2、q3、……、qc采样时刻中的任一时刻进行步骤二至步骤三的过程)M个阵元的回波的多普勒谱Fq1(s),使用正交投影算法在频域进行一阶海杂波抑制,得到抑制后的回波信号多普勒矩阵FDPRS(s);
步骤三:对FDPRS(s)矩阵按列求和,得到求和后的序列检测序列中的峰值得到l个峰值在序列中的序号Ipeak(1)、Ipeak(2)、…、Ipeak(l)以及峰值对应的频率speak(1)、speak(2)、…、speak(l);
步骤四:选取正整数L作为天线阵列的子阵阵元个数,将步骤一中得到的q1采样时刻M个阵元的回波的多普勒矩阵Fq1(s)按照第1行到第L行、第2r+1行到第(2r+L)行、第4r+1行到第(4r+L)行、第6r+1行到第(6r+L)行……的规律重新写成[(M-L)/2r+1]([(M-L)/2r+1]也可以表示为)个部分,“[]”表示向下取整。由[(M-L)/2r+1]个部分组成天线阵列所有子阵的回波的多普勒矩阵Fsub(s);令多普勒矩Fsub(s)中多普勒频率为speak(1)的列上的数据保持不变,其余位置的值为0,对Fsub(s)中的[(M-L)/2r+1]个部分别使用正交投影算法进行一阶海杂波抑制,得到抑制后的天线阵列所有子阵的回波的多普勒谱其中L≤2r,r<M/2,r为天线阵列在相邻两个脉冲周期之间走过的阵元个数;[(M-L)/2r+1]中的[]表示向下取整;
步骤五:对步骤四中得到的抑制后的天线阵列所有子阵的回波的多普勒矩阵的[(M-L)/2r+1]个部分的每一部分分别按列求和,得到一个[(M-L)/2r+1]×R大小的中间变量矩阵对的每一行做逆傅里叶变换(IFFT)得到q1采样时刻的等效阵元的接收信号xΣq1(n),其中n表示第n个脉冲周期,1≤n≤N,N为积累脉冲数;
步骤六:依次取每个脉冲周期内q2、q3、……、qc采样时刻的M个阵元的回波的多普勒矩阵Fq2(s),Fq3(s),…,Fqc(s),对speak(1),重复执行步骤四和步骤五,得到q2、q3、……、qc采样时刻的等效阵元的接收信号xΣq2(n)、xΣq3(n)、……、xΣqc(n),根据xΣq1(n)、xΣq2(n)、……、xΣqc(n)得到海杂波抑制后的等效阵列的接收信号XΣ(n);
步骤七:使用重叠相关算法对信号矩阵XΣ(n)进行处理,得到多普勒频率为speak(1)的目标的虚拟孔径波束形成结果y1(θ),其中θ表示角度,角标1表示第一个目标;
步骤八:对步骤三中得到的频率为speak(2)、……、speak(l),重复执行步骤四至步骤七,直至得到第2、3、……、l个目标的虚拟孔径波束形成结果y2(θ)、y3(θ)、……、yl(θ),完成对舰载高频地波雷达的l个目标的方位高分辨。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中得到每个采样时刻的M个阵元的回波的多普勒矩阵Fq1(s),Fq2(s),…,Fqc(s)的具体过程为:
设天线阵列的阵元之间的间距为d,雷达信号的波长为λ,脉冲宽度为τ,d需要满足d≤λ/2,阵列平台速度为v,雷达的天线阵列接收到的信号为:
xn,m(q)为第n个脉冲重复周期内,雷达天线阵列的第m个阵元的接收信号
取矩阵X(q)中的第一列,即取第1个阵元的接收信号:
取每一个脉冲周期内q1、q2……、qc采样时刻的信号分别作快速傅里叶变换(FFT),得到第一个阵元接收的回波信号所有采样时刻的多普勒谱:
f1,1(s)=FFT([x1,1(q1),x2,1(q1),...,xN,1(q1)])
……
f1,c(s)=FFT([x1,1(qc),x2,1(qc),...,xN,1(qc)])
式中,f1,1(s)、f1,2(s)、……、f1,c(s)分别为第1个阵元的q1、q2……、qc采样时刻对应的多普勒谱,s为频率采样序列,s=[s1,s2,...,sR],其中s1、s2…sR为频域采样,R为总的多普勒单元数,由于FFT前后序列个数不变,因此多普勒单元总数等于时域脉冲积累总数,因此R=N;s=[s1,s2,...,sN]∈[0,fs];
依次取矩阵X(q)中的第2列、第3列、……、第M列,得到每个采样时刻的M个阵元的回波的多普勒矩阵Fq1(s),Fq2(s),…,Fqc(s):
……
式中,Fq1(s)、Fqc(s)分别为雷达阵列在q1、qc采样时刻上的多普勒谱矩阵,fm,1(s)、fm,c(s)分别为第m个阵元在q1、qc采样时刻的回波信号的多普勒谱。
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤二中对步骤一得到的q1采样时刻M个阵元的回波的多普勒矩阵Fq1(s),使用正交投影算法在频域进行一阶海杂波抑制,得到抑制后的回波信号多普勒矩阵FDPRS(s)的具体过程为:
式中,FDPRS(s)为使用正交投影算法抑制后的回波信号多普勒矩阵,w(1)P×M、w(2)P×M、……、w(N)P×M分别为第1、2、……、R个多普勒单元所对应的正交加权因子,p为角度序列,p=[p1,p2,…,pk]∈[0,180],k为角度序列的元素总个数,[fp1(s1),fp2(s1),...,fpk(s1)]T=w(1)p×M×[f1,1(s1),f2,1(s1),...,fm,1(s1),...,fM,1(s1)]T,fp1(s1),fp2(s1),...,fpk(s1)分别为多普勒频率为s1的信号在p1、p2……、pk角度上的波束形成响应。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤三中具体为:
为fp1(s1),fp2(s1),...,fpk(s1)求和后的值。
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述步骤四中由[(M-L)/2r+1]个部分组成天线阵列的所有子阵的回波的多普勒矩阵Fsub(s)具体为:
其它步骤及参数与具体实施方式一至四之一相同。
具体实施方式六:本实施方式与具体实施方式一至五之一不同的是:所述步骤四中对Fsub(s)中的[(M-L)/2r+1]个部分别使用正交投影算法进行一阶海杂波抑制,得到抑制后的天线阵列所有子阵的回波的多普勒矩阵的具体为:
针对speak(1)这一多普勒频率,令矩阵Fsub(s)中多普勒频率为speak(1)的列上的数据保持不变,其余位置的值为0,再令j=Ipeak(1),对Fsub(s)中的[(M-L)/2r+1]个部分别使用正交投影算法进行一阶海杂波抑制,得到抑制后的天线阵列所有子阵的回波的多普勒谱
式中,j=Ipeak(1),w(j)p×M为第j个多普勒单元所对应的正交加权因子,sj为第j个多普勒单元对应的频率,频率为sj的信号在p1角度上第一个子阵的波束形成响应,频率为sj的信号在p1角度上第二个子阵的波束形成响应。
其它步骤及参数与具体实施方式一至五之一相同。
具体实施方式七:本实施方式与具体实施方式一至六之一不同的是:所述步骤五中对步骤四中得到的抑制后的天线阵列所有子阵的回波的多普勒矩阵的[(M-L)/2r+1]个部分的每一部分分别按列求和,得到一个[(M-L)/2r+1]×R大小的中间变量矩阵具体为:
对的每一行做逆傅里叶变换(IFFT),得到q1采样时刻的等效阵元的接收信号xΣq1(n);
其中IFFT为快速逆傅里叶变换,x1,q1(1),x1,q1(2),…,x1,q1(N)为第一个阵元在第1、2、……、N个脉冲周期内q1采样时刻的接收信号。
其它步骤及参数与具体实施方式一至六之一相同。
具体实施方式八:本实施方式与具体实施方式一至七之一不同的是:所述步骤六中根据得到的q2、q3、……、qc采样时刻的等效阵元的接收信号xΣq2(n)、xΣq3(n)、……、xΣqc(n)得到海杂波抑制后的等效阵列的接收信号XΣ(n)(该等效阵列阵元数为[(M-L)/2r+1],脉冲积累总数为N,每一个脉冲周期内的时域采样总为c)的具体过程为:
XΣ(n)=[xΣq1(n),xΣq2(n),…,xΣqi(n),…,xΣqc(n)]c×N×[(M-L)/2r+1]
公式xΣqi(n)中,角标1、2、3、……、[(M-L)/2r+1]表示等效的阵元序号,括号内数字1、2、……、N表示脉冲周期,角标qi表示每一个脉冲周期内的第i个时域采样。
其它步骤及参数与具体实施方式一至七之一相同。
具体实施方式九:本实施方式与具体实施方式一至八之一不同的是:所述步骤七中重叠相关算法的具体过程为:
步骤七一:将XΣ(n)按照每个阵元接收信号写成[(M-L)/2r+1]个二维矩阵的形式:
……
式中,X1,Σ(n)、X2,Σ(n)、……、X[(M-L)/2r+1],Σ(n)分别表示第1、2、……、[(M-L)/2r+1]个等效阵元在各个脉冲周期内的信号,矩阵的每一行表示一个脉冲周期的信号,依次用第j+1号等效阵元的第n个脉冲周期的接收信号与第j号等效阵元的第n+1个脉冲周期的接收信号进行相关运算并对所有阵元取均值,求得相邻周期间两个信号之间的相位差φ;其中1≤j≤[(M-L)/2r+1],n取1、2、……、N-1;
步骤七二:用步骤七一中求得的相关因子分别对X2,Σ(n)、X3,Σ(n)、……、X[(M-L)/2r+1],Σ(n)中的第2个脉冲周期至第N个脉冲周期的信号进行相位补偿,得到
……
取X1,Σ(n)、的第1行,的第2行到[(M-L)/2r+1]行,构成一个新的信号矩阵X(m,q);
式中,m表示行,1≤m≤[(M-L)/2r+N],q表示列,1≤q≤c;
步骤三:取X(m,q)中的第一列数据,按照阵元数为[(M-L)/2r+N],阵元间距为4d,波长为λ进行波束形成,得到波束形成结果y1(θ);其中θ表示角度,0≤θ≤180;
其它步骤及参数与具体实施方式一至八之一相同。
采用以下实施例验证本发明的有益效果:
实施例一:
步骤一:仿真参数如下,信号频率f0=5MHz,阵元数为3M=12,阵元之间的间距为d=8,雷达信号的波长为λ=60m,脉冲宽度τ=0.4μs,阵列平台速度为v=4m/s,脉冲重复周期Tp=4s,积累脉冲数为N=50。采样频率fs=500MHz。目标1角度为40度,径向速度为3.2m/s,目标2角度为50度,径向速度为2.1m/s,信噪比为3dB。雷达的天线阵接收到的信号为:
式中,X(q)为雷达天线阵的接收信号,q为采样时刻序列,q=[q1,q2,...,q201]∈[0,fsτ],其中q1、q2…q201为采样时刻,时域采样总数为201。取矩阵X(q)中的第一列:
取每一个脉冲周期内q1、q2……、qc时刻的信号分别作快速傅里叶变换(FFT),得到第一个阵元接收的回波信号所有时刻的多普勒谱:
f1,1(s)=FFT([x1,1(q1),x2,1(q1),...,x50,1(q1)])
……
f1,c(s)=FFT([x1,1(qc),x2,1(qc),...,x50,1(q201)])
式中,f1,1(s)、f1,2(s)、……、f1,c(s)分别为第1个阵元的q1、q2……、qc时刻对应的多普勒谱,s为频率采样序列,s=[s1,s2,...,s50],其中s1、s2…s50为频域采样,50为总的多普勒单元数。f1(s)的曲线图如图1所示。
依次取矩阵X(q)中的第2列、第3列、……、第12列,按照步骤一中的方法,得到所有阵元的回波的多普勒谱:
……
式中,Fq1(s)、Fq201(s)分别为雷达阵列在q1、q201时刻上的多普勒谱矩阵,fm,1(s)、fm,201(s)分别为第m个阵元在q1、q201时刻的回波信号的多普勒谱
式中,Fq1(s)为雷达阵列接收信号的多普勒谱矩阵,fm,1(s)为第m个阵元在q1时刻的的回波信号的多普勒谱。
步骤二:从0到180度,以2度为步长,使用正交投影算法在频域进行一阶海杂波抑制,得到抑制后的回波信号多普勒矩阵:
式中,FDPRS(s)为使用正交投影算法抑制后的回波信号多普勒矩阵,w(1)91×12、w(2)91×12、……、w(N)91×12分别为第1、2、……、N个多普勒单元所对应的正交加权因子,p为角度序列,p=[p1,p2,…,p91]=[0,2,4,...,180],91为角度序列的元素总个数,[fp1(s1),fp2(s1),...,fp91(s1)]T=w(1)91×12*[f1,1(s1),f2,1(s1),...,fm,1(s1),...,f12,1(s1)]T。FDPRS(s)的结果图如图2所示
步骤三:对FDPRS(s)矩阵按列求和,得到一个新的序列:
式中为求和后的新序列,的曲线图如图3所示。检测序列中的峰值,得到2个峰值在序列中的序号14、18以及所对应的频率-0.205Hz、-0.155Hz。
步骤四:将步骤二中所得的矩阵F(s)分为三部分
针对-0.205Hz这一多普勒频率,令矩阵F(s)中多普勒频率为-0.205Hz的列上的数据保持不变,其余位置的值为0,对上式中的三个部分别使用正交投影算法进行一阶海杂波抑制:
式中,w(7)91×4为第7个多普勒单元所对应的正交加权因子,-0.205Hz为第7个多普勒单元对应的频率,的曲线图如图4所示。
步骤五:对步骤五中得到的三个部分的每一部分分别按列求和,得到一个3×50大小的矩阵其中50为总的多普勒单元数。
对的每一行做逆傅里叶变换(IFFT),得到q1时刻的等效阵元的接收信号xΣq1(n),其中n表示第n个脉冲周期,1≤n≤50。
xΣq1(n)的实部曲线如图5所示。
步骤六:依次取每个脉冲周期内q2、q3、……、qc时刻的12个阵元的回波的多普勒谱Fq2(s),Fq3(s),…,Fq201(s),对-0.205这个多普勒频率,重复执行步骤四和步骤五,得到q2、q3、……、q201时刻的等效阵元的接收信号;xΣq2(n)、xΣq3(n)、……、xΣqi(n)、……、xΣq201(n),其中201为时域采样总数。
XΣ(n)=[xΣq1(n),xΣq2(n),...,xΣqi(n),...,xΣq201(n)]201×25×3
其中
式中,角标1、2、3表示等效的阵元序号,括号内数字1、2、……、50表示脉冲周期,角标qi表示每一个脉冲周期内的第i个时域采样。
步骤七:使用重叠相关算法对信号矩阵XΣ(n)进行处理,得到多普勒频率为-0.205Hz的这一个目标的虚拟孔径波束形成结果y1(θ),其中θ表示角度,角标1表示第一个目标,y1(θ)的曲线示意图如图6所示。
步骤八:对步骤三中得到的频率为-0.1600Hz这个多普勒频率,重复执行步骤四至步骤七,即可得到第2个目标的虚拟孔径波束形成结果y2(θ),y1(θ)的曲线示意图如图7所示。至此完成了对舰载高频地波雷达的多个目标的方位高分辨。
重叠相关算法对信号矩阵XΣ(n)的进行处理的步骤如下所示:
步骤一:将XΣ(n)按照每个阵元接收信号写成三个二维矩阵的形式:
式中,X1,Σ(n)、X2,Σ(n)、X3,Σ(n)分别表示第1、2、3个等效阵元在各个脉冲周期内的信号,矩阵的一行表示一个脉冲周期的信号。依次用第1个等效阵元的第n个(n取1、2、……、49)脉冲周期的接收信号与第3个等效阵元的第n+1个脉冲周期的接收信号进行相关运算,求得两个信号之间的相位差φ。
步骤二:用步骤一中求得的相关因子分别对X1,Σ(n)、X2,Σ(n)、X3,Σ(n)中的第二个脉冲周期至第N个脉冲周期的信号进行相位补偿,得到
取X1,Σ(n)、的第1行、的第2行到第50行,构成一个新的信号矩阵X(m,q)。
式中,m表示行,1≤m≤52;q表示列,1≤q≤201,201为时域采样总数。
步骤七三:取X(m,q)中的第一列数据,按照阵元数为52,阵元间距为32,波长为60这些参数进行波束形成,得到波束形成结果y1(θ),其中θ表示角度,0≤θ≤180。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (9)
1.一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法,其特征在于:所述基于虚拟孔径的舰载高频地波雷达方位高分辨方法包括以下步骤:
步骤一:设舰载高频地波雷达的天线阵列为线阵,阵元数为M,M为大于3的正整数,雷达的天线阵列接收到的信号为X(q),对矩阵X(q)中的每一个脉冲周期内的q1、q2、q3、……、qc采样时刻的信号分别作FFT,得到每个采样时刻的M个阵元的回波的多普勒谱Fq1(s),Fq2(s),…,Fqc(s);所述FFT为快速傅里叶变换,q为采样时刻序列,q=[q1,q2,...,qc],其中q1、q2…qc为采样时刻,c为每一个脉冲周期内时域采样总数;
步骤二:对步骤一得到的q1采样时刻M个阵元的回波的多普勒谱Fq1(s),使用正交投影算法在频域进行一阶海杂波抑制,得到抑制后的回波信号多普勒矩阵FDPRS(s);
步骤三:对FDPRS(s)矩阵按列求和,得到求和后的序列检测序列中的峰值得到l个峰值在序列中的序号Ipeak(1)、Ipeak(2)、…、Ipeak(l)以及峰值对应的频率speak(1)、speak(2)、…、speak(l);
步骤四:选取正整数L作为天线阵列的子阵阵元个数,将步骤一中得到的q1采样时刻M个阵元的回波的多普勒矩阵Fq1(s)按照第1行到第L行、第2r+1行到第(2r+L)行、第4r+1行到第(4r+L)行、第6r+1行到第(6r+L)行……的规律重新写成[(M-L)/2r+1]个部分,由[(M-L)/2r+1]个部分组成天线阵列所有子阵的回波的多普勒矩阵Fsub(s);令多普勒矩阵Fsub(s)中多普勒频率为speak(1)的列上的数据保持不变,其余位置的值为0,对Fsub(s)中的[(M-L)/2r+1]个部分别使用正交投影算法进行一阶海杂波抑制,得到抑制后的天线阵列所有子阵的回波的多普勒矩阵其中L≤2r,r<M/2,r为天线阵列在相邻两个脉冲周期之间走过的阵元个数;[(M-L)/2r+1]中的[]表示向下取整;
步骤五:对步骤四中得到的抑制后的天线阵列所有子阵的回波的多普勒矩阵的[(M-L)/2r+1]个部分的每一部分分别按列求和,得到一个[(M-L)/2r+1]×R大小的中间变量矩阵对的每一行做逆傅里叶变换得到q1采样时刻的等效阵元的接收信号xΣq1(n),其中n表示第n个脉冲周期,1≤n≤N,N为积累脉冲数;
步骤六:依次取每个脉冲周期内q2、q3、……、qc时刻的M个阵元的回波的多普勒矩阵Fq2(s),Fq3(s),…,Fqc(s),对speak(1),重复执行步骤四和步骤五,得到q2、q3、……、qc采样时刻的等效阵元的接收信号xΣq2(n)、xΣq3(n)、……、xΣqc(n),根据xΣq1(n)、xΣq2(n)、……、xΣqc(n)得到海杂波抑制后的等效阵列的接收信号XΣ(n);
步骤七:使用重叠相关算法对信号矩阵XΣ(n)进行处理,得到多普勒频率为speak(1)的目标的虚拟孔径波束形成结果y1(θ),其中θ表示角度,角标1表示第一个目标;
步骤八:对步骤三中得到的频率为speak(2)、……、speak(l),重复执行步骤四至步骤七,直至得到第2、3、……、l个目标的虚拟孔径波束形成结果y2(θ)、y3(θ)、……、yl(θ),完成对舰载高频地波雷达的l个目标的方位分辨。
2.根据权利要求1所述的一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法,其特征在于:所述步骤一中得到每个采样时刻的M个阵元的回波的多普勒矩阵Fq1(s),Fq2(s),…,Fqc(s)的具体过程为:
设天线阵列的阵元之间的间距为d,雷达信号的波长为λ,脉冲宽度为τ,d需要满足d≤λ/2,阵列平台速度为v,雷达的天线阵列接收到的信号为:
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mi>M</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
xn,m(q)为第n个脉冲重复周期内,雷达天线阵列的第m个阵元的接收信号取矩阵X(q)中的第一列:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
取每一个脉冲周期内q1、q2……、qc采样时刻的信号分别作快速傅里叶变换(FFT),得到第一个阵元接收的回波信号所有采样时刻的多普勒谱:
f1,1(s)=FFT([x1,1(q1),x2,1(q1),...,xN,1(q1)])
……
f1,c(s)=FFT([x1,1(qc),x2,1(qc),...,xN,1(qc)])
式中,f1,1(s)、f1,2(s)、……、f1,c(s)分别为第1个阵元的q1、q2……、qc采样时刻对应的多普勒谱,s为频率采样序列,s=[s1,s2,...,sR],其中s1、s2…sR为频域采样,R为总的多普勒单元数,R=N;
依次取矩阵X(q)中的第2列、第3列、……、第M列,得到每个采样时刻的M个阵元的回波的多普勒矩阵Fq1(s),Fq2(s),…,Fqc(s):
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>......</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
2
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中,Fq1(s)、Fqc(s)分别为雷达阵列在q1、qc采样时刻上的多普勒谱矩阵,fm,1(s)、fm,c(s)分别为第m个阵元在q1、qc采样时刻的回波信号的多普勒谱。
3.根据权利要求2所述的一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法,其特征在于:所述步骤二中对步骤一得到的q1采样时刻M个阵元的回波的多普勒矩阵Fq1(s),使用正交投影算法在频域进行一阶海杂波抑制,得到抑制后的回波信号多普勒矩阵FDPRS(s)的具体过程为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>D</mi>
<mi>P</mi>
<mi>R</mi>
<mi>S</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<mrow>
<mi>w</mi>
<msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>&times;</mo>
<mi>M</mi>
</mrow>
</msub>
<mo>*</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>w</mi>
<msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>&times;</mo>
<mi>M</mi>
</mrow>
</msub>
<mo>*</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>w</mi>
<msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>&times;</mo>
<mi>M</mi>
</mrow>
</msub>
<mo>*</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>k</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>k</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>k</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,FDPRS(s)为使用正交投影算法抑制后的回波信号多普勒矩阵,w(1)P×M、w(2)P×M、……、w(N)P×M分别为第1、2、……、R个多普勒单元所对应的正交加权因子,p为角度序列,p=[p1,p2,...,pk],k为角度序列的元素总个数,[fp1(s1),fp2(s1),...,fpk(s1)]T=w(1)p×M×[f1,1(s1),f2,1(s1),...,fm,1(s1),...,fM,1(s1)]T,fp1(s1),fp2(s1),...,fpk(s1)分别为多普勒频率为s1的信号在p1、p2……、pk角度上的波束形成响应。
4.根据权利要求3所述的一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法,其特征在于:所述步骤三中具体为:
<mrow>
<msub>
<mover>
<mi>F</mi>
<mo>&Sigma;</mo>
</mover>
<mrow>
<mi>D</mi>
<mi>P</mi>
<mi>R</mi>
<mi>S</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
求和后的值。
5.根据权利要求4所述的一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法,其特征在于:所述步骤四中由[(M-L)/2r+1]个部分组成天线阵列的所有子阵的回波的多普勒矩阵Fsub(s)具体为:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>s</mi>
<mi>u</mi>
<mi>b</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mi>L</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
<mtr>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
<mtr>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>L</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>L</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>L</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
<mtr>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
<mtr>
<mtd>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mi>L</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mi>L</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>f</mi>
<mrow>
<mn>4</mn>
<mi>r</mi>
<mo>+</mo>
<mi>L</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>R</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
6.根据权利要求5所述的一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法,其特征在于:所述步骤四中对Fsub(s)中的[(M-L)/2r+1]个部分别使用正交投影算法进行一阶海杂波抑制,得到抑制后的天线阵列所有子阵的回波的多普勒矩阵的具体为:
式中,j=Ipeak(1),w(j)p×M为第j个多普勒单元所对应的正交加权因子,sj为第j个多普勒单元对应的频率,频率为sj的信号在p1角度上第一个子阵的波束形成响应,频率为sj的信号在p1角度上第二个子阵的波束形成响应。
7.根据权利要求6所述的一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法,其特征在于:所述步骤五中对步骤四中得到的抑制后的天线阵列所有子阵的回波的多普勒矩阵的[(M-L)/2r+1]个部分的每一部分分别按列求和,得到一个[(M-L)/2r+1]×R大小的中间变量矩阵具体为:
<mrow>
<msub>
<mover>
<mi>F</mi>
<mo>*</mo>
</mover>
<mrow>
<mo>&Sigma;</mo>
<mo>-</mo>
<mi>D</mi>
<mi>P</mi>
<mi>R</mi>
<mi>S</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mn>0</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mover>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>*</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mn>0</mn>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mn>0</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mover>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>*</mo>
<mo>*</mo>
</mrow>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mn>0</mn>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mn>0</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mover>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>*</mo>
<mo>*</mo>
<mo>*</mo>
</mrow>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mn>0</mn>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>......</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
对的每一行做逆傅里叶变换,得到q1采样时刻的等效阵元的接收信号xΣq1(n);
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mo>&Sigma;</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mn>0</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mover>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>*</mo>
</mover>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mn>0</mn>
<mo>&rsqb;</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mn>0</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mover>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>*</mo>
<mo>*</mo>
</mrow>
</mover>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mn>0</mn>
<mo>&rsqb;</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mn>0</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mover>
<msub>
<mi>f</mi>
<mrow>
<mi>p</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>*</mo>
<mo>*</mo>
<mo>*</mo>
</mrow>
</mover>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mn>0</mn>
<mo>&rsqb;</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中IFFT为快速逆傅里叶变换,x1,q1(1),x1,q1(2),...,x1,q1(N)为第一个阵元在第1、2、……、N个脉冲周期内q1采样时刻的接收信号。
8.根据权利要求7所述的一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法,其特征在于:所述步骤六中根据得到的q2、q3、……、qc采样时刻的等效阵元的接收信号xΣq2(n)、xΣq3(n)、……、xΣqc(n)得到海杂波抑制后的等效阵列的接收信号XΣ(n)的具体过程为:
XΣ(n)=[xΣq1(n),xΣq2(n),...,xΣqi(n),...,xΣqc(n)]c×N×[(M-L)/2r+1]
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>&Sigma;</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
公式xΣqi(n)中,角标1、2、3、……、[(M-L)/2r+1]表示等效的阵元序号,括号内数字1、2、……、N表示脉冲周期,角标qi表示每一个脉冲周期内的第i个时域采样。
9.根据权利要求8所述的一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法,其特征在于:所述步骤七中重叠相关算法的具体过程为:
步骤七一:将XΣ(n)按照每个阵元接收信号写成[(M-L)/2r+1]个二维矩阵的形式:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mo>&Sigma;</mo>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mo>&Sigma;</mo>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>......</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mo>&Sigma;</mo>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>2</mn>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,X1,Σ(n)、X2,Σ(n)、……、X[(M-L)/2r+1],Σ(n)分别表示第1、2、……、[(M-L)/2r+1]个等效阵元在各个脉冲周期内的信号,矩阵的每一行表示一个脉冲周期的信号,依次用第j+1号等效阵元的第n个脉冲周期的接收信号与第j号等效阵元的第n+1个脉冲周期的接收信号进行相关运算并对所有阵元取均值,求得相邻周期间两个信号之间的相位差φ;其中1≤j≤[(M-L)/2r+1],n取1、2、……、N-1;
步骤七二:用步骤七一中求得的相关因子分别对X2,Σ(n)、X3,Σ(n)、……、X[(M-L)/2r+1],Σ(n)中的第2个脉冲周期至第N个脉冲周期的信号进行相位补偿,得到
取X1,Σ(n)、的第1行,的第2行到[(M-L)/2r+1]行,构成一个新的信号矩阵X(m,q);
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>3</mn>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>j&phi;</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>j&phi;</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>q</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>j&phi;</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>,</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>,</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>,</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>,</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>,</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>,</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>,</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>,</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>,</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>,</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>,</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>,</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>&lsqb;</mo>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mi>N</mi>
<mo>&rsqb;</mo>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>&lsqb;</mo>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mi>N</mi>
<mo>&rsqb;</mo>
<mo>,</mo>
<mn>2</mn>
<mo>)</mo>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>&lsqb;</mo>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mi>N</mi>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,m表示行,1≤m≤[(M-L)/2r+N],q表示列,1≤q≤c;
步骤七三:取X(m,q)中的第一列数据,按照阵元数为[(M-L)/2r+N],阵元间距为4d,波长为λ进行波束形成,得到波束形成结果y1(θ);
<mrow>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>2</mn>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mfrac>
<mrow>
<mn>4</mn>
<mi>n</mi>
<mi>d</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>&lambda;</mi>
</mfrac>
</mrow>
</msup>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
9
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710428463.1A CN107064904B (zh) | 2017-06-08 | 2017-06-08 | 一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710428463.1A CN107064904B (zh) | 2017-06-08 | 2017-06-08 | 一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107064904A true CN107064904A (zh) | 2017-08-18 |
CN107064904B CN107064904B (zh) | 2019-09-13 |
Family
ID=59617565
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710428463.1A Active CN107064904B (zh) | 2017-06-08 | 2017-06-08 | 一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107064904B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107807352A (zh) * | 2017-09-30 | 2018-03-16 | 武汉大学 | 一种海上平台高频雷达阵列不变波束形成方法 |
CN108594190A (zh) * | 2018-04-20 | 2018-09-28 | 西安电子科技大学 | 一种高分辨海杂波的仿真方法 |
RU2713503C1 (ru) * | 2019-08-19 | 2020-02-05 | Федеральное государственное унитарное предприятие "Ростовский-на-Дону научно-исследовательский институт радиосвязи" (ФГУП "РНИИРС") | Способ углового сверхразрешения в приемных цифровых антенных решетках |
WO2022252311A1 (zh) * | 2021-06-04 | 2022-12-08 | 惠州市德赛西威汽车电子股份有限公司 | 一种估算目标模糊速度的方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011033344A (ja) * | 2009-07-29 | 2011-02-17 | Mitsubishi Electric Corp | レーダ装置 |
CN103616690A (zh) * | 2013-12-11 | 2014-03-05 | 哈尔滨工业大学 | 基于船载高频地波超视距雷达的海面风向提取方法 |
US20140197984A1 (en) * | 2013-01-15 | 2014-07-17 | Raytheon Canada Limited | Virtual antenna extension for sampled aperture arrays |
CN103983952A (zh) * | 2014-05-16 | 2014-08-13 | 哈尔滨工程大学 | 一种非圆信号双基地mimo雷达低复杂度收发角度联合估计方法 |
CN104155632A (zh) * | 2014-07-18 | 2014-11-19 | 南京航空航天大学 | 一种基于局部相关性的改进子空间海杂波抑制方法 |
CN104391279A (zh) * | 2014-11-24 | 2015-03-04 | 哈尔滨工业大学 | 基于电离层传播特性的相径扰动抑制方法 |
CN105445718A (zh) * | 2015-11-19 | 2016-03-30 | 哈尔滨工业大学 | 一种基于阵列重构的分布式多载舰超视距雷达的doa估计方法 |
CN106443633A (zh) * | 2016-12-01 | 2017-02-22 | 哈尔滨工业大学 | 船载高频地波雷达海杂波时域抑制方法 |
-
2017
- 2017-06-08 CN CN201710428463.1A patent/CN107064904B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011033344A (ja) * | 2009-07-29 | 2011-02-17 | Mitsubishi Electric Corp | レーダ装置 |
US20140197984A1 (en) * | 2013-01-15 | 2014-07-17 | Raytheon Canada Limited | Virtual antenna extension for sampled aperture arrays |
CN103616690A (zh) * | 2013-12-11 | 2014-03-05 | 哈尔滨工业大学 | 基于船载高频地波超视距雷达的海面风向提取方法 |
CN103983952A (zh) * | 2014-05-16 | 2014-08-13 | 哈尔滨工程大学 | 一种非圆信号双基地mimo雷达低复杂度收发角度联合估计方法 |
CN104155632A (zh) * | 2014-07-18 | 2014-11-19 | 南京航空航天大学 | 一种基于局部相关性的改进子空间海杂波抑制方法 |
CN104391279A (zh) * | 2014-11-24 | 2015-03-04 | 哈尔滨工业大学 | 基于电离层传播特性的相径扰动抑制方法 |
CN105445718A (zh) * | 2015-11-19 | 2016-03-30 | 哈尔滨工业大学 | 一种基于阵列重构的分布式多载舰超视距雷达的doa估计方法 |
CN106443633A (zh) * | 2016-12-01 | 2017-02-22 | 哈尔滨工业大学 | 船载高频地波雷达海杂波时域抑制方法 |
Non-Patent Citations (2)
Title |
---|
李高鹏 等: "基于海杂波保持的高频地波雷达干扰抑制方法研究", 《电波科学报》 * |
李高鹏 等: "高频地波雷达干扰抑制方法研究", 《电子学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107807352A (zh) * | 2017-09-30 | 2018-03-16 | 武汉大学 | 一种海上平台高频雷达阵列不变波束形成方法 |
CN107807352B (zh) * | 2017-09-30 | 2020-07-10 | 武汉大学 | 一种海上平台高频雷达阵列不变波束形成方法 |
CN108594190A (zh) * | 2018-04-20 | 2018-09-28 | 西安电子科技大学 | 一种高分辨海杂波的仿真方法 |
CN108594190B (zh) * | 2018-04-20 | 2021-04-27 | 西安电子科技大学 | 一种高分辨海杂波的仿真方法 |
RU2713503C1 (ru) * | 2019-08-19 | 2020-02-05 | Федеральное государственное унитарное предприятие "Ростовский-на-Дону научно-исследовательский институт радиосвязи" (ФГУП "РНИИРС") | Способ углового сверхразрешения в приемных цифровых антенных решетках |
WO2022252311A1 (zh) * | 2021-06-04 | 2022-12-08 | 惠州市德赛西威汽车电子股份有限公司 | 一种估算目标模糊速度的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107064904B (zh) | 2019-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Klemm | Principles of space-time adaptive processing | |
CN101858976B (zh) | 一种基于单通道sar多子孔径的运动目标检测方法 | |
Huang et al. | Multichannel sea clutter modeling for spaceborne early warning radar and clutter suppression performance analysis | |
CN107064904B (zh) | 一种基于虚拟孔径的舰载高频地波雷达方位高分辨方法 | |
CN101556328A (zh) | 基于杂波协方差矩阵的机载雷达空时二维滤波器构建方法 | |
CN104977571B (zh) | 基于俯仰频率分集stap的距离模糊杂波抑制方法 | |
CN101482610A (zh) | 相控阵雷达自适应抑制特殊干扰技术 | |
CN104360325A (zh) | 机载前视阵雷达的空时自适应处理方法 | |
CN107966677B (zh) | 一种基于空间稀疏约束的圆阵模态域方位估计方法 | |
CN105223560A (zh) | 基于杂波俯仰方位谱稀疏恢复的机载雷达目标检测方法 | |
CN1752771A (zh) | 单通道相控阵接收信号重构及空间信号处理方法 | |
CN109375179A (zh) | 一种基于稀疏表示的机载雷达近程强杂波抑制方法 | |
CN109765529B (zh) | 一种基于数字波束形成的毫米波雷达抗干扰方法及*** | |
CN1088199C (zh) | 空时二维多波束自适应信号处理方法 | |
CN103760529A (zh) | 一种无源探测高效级联空时自适应处理方法 | |
CN110082744A (zh) | 多普勒分级复用的mimo机载双基雷达杂波抑制方法 | |
CN103760540B (zh) | 基于重构信号和1-范数的动目标检测与参数估计方法 | |
Zhou et al. | Signal model and moving target detection based on MIMO synthetic aperture radar | |
CN109597034B (zh) | 一种基于欧几里得距离的空时自适应处理方法 | |
CN104035078A (zh) | 一种基于阵元阶数递推的降维空时自适应权值计算方法 | |
Zenghui et al. | Local degrees of freedom of clutter for airborne space–time adaptive processing radar with subarrays | |
Riddolls | Comparison of linear and planar arrays for auroral clutter control in an over-the-horizon radar | |
Zhu et al. | Application of frequency diverse array to resolve range ambiguity for SAR imaging | |
Mei et al. | Nonuniform three-dimensional configuration distributed SAR signal reconstruction clutter suppression | |
Klemm | Adaptive airborne MTI with two-dimensional motion compensation |
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 |