CN111610498A - 一种高自由度去耦合的空时自适应信号处理方法 - Google Patents
一种高自由度去耦合的空时自适应信号处理方法 Download PDFInfo
- Publication number
- CN111610498A CN111610498A CN202010574966.1A CN202010574966A CN111610498A CN 111610498 A CN111610498 A CN 111610498A CN 202010574966 A CN202010574966 A CN 202010574966A CN 111610498 A CN111610498 A CN 111610498A
- Authority
- CN
- China
- Prior art keywords
- virtual
- clutter
- array
- space
- frequency
- 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
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 17
- 238000003672 processing method Methods 0.000 title claims abstract description 13
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 59
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 56
- 239000013598 vector Substances 0.000 claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims description 65
- 238000003491 array Methods 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 8
- 238000006880 cross-coupling reaction Methods 0.000 claims description 4
- 230000002123 temporal effect Effects 0.000 claims description 4
- 101100134058 Caenorhabditis elegans nth-1 gene Proteins 0.000 claims description 2
- 230000008878 coupling Effects 0.000 abstract description 19
- 238000010168 coupling process Methods 0.000 abstract description 19
- 238000005859 coupling reaction Methods 0.000 abstract description 19
- 238000001514 detection method Methods 0.000 abstract description 7
- 238000004422 calculation algorithm Methods 0.000 description 17
- 230000001808 coupling effect Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 6
- 238000005070 sampling Methods 0.000 description 5
- 239000002131 composite material Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 102000001694 Signal Transducing Adaptor Proteins Human genes 0.000 description 1
- 108010029228 Signal Transducing Adaptor Proteins Proteins 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Radar, Positioning & Navigation (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Remote Sensing (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种高自由度去耦合的空时自适应信号处理方法,针对互质结构单频工作模式下,舍弃差分合成结构中非均匀虚拟部分引起的空时自适应信号滤波器权矢量估计性能损失的问题,提出了一种利用双频工作模式提升互质结构空时自适应信号权矢量估计精度的方法,通过选择合适的单一额外工作频率获取额外不同位置的差分合成虚拟阵元和差分合成虚拟脉冲,进而填充参考工作频率下互质结构差分合成阵和差分合成脉冲中两个缺失的虚拟阵元和脉冲。相比于单频工作模式,双频工作模式获取了更多均匀分布的虚拟阵元和脉冲,滤波器的自由度得到了提升,改善了信号检测性能和角度多普勒分辨率,产生的互质结构可以缓解相互耦合,增强滤波器的虚拟自由度。
Description
技术领域
本发明属于雷达技术领域,尤其涉及一种高自由度去耦合的空时自适应信号处理方法。
背景技术
空时自适应处理(STAP)可以提高抑制杂波和检测目标的能力,在机载雷达中起着至关重要的作用。传统的算法,如降阶、降维和参数化模型,通常基于精确的空时转向矢量为研究基础,在阵元之间相互耦合的情况下,它们的性能显著下降,甚至无法工作。在实际雷达系统中,需要考虑相互耦合效应,因此,考虑相互耦合的STAP是阵列信号处理中的热点问题之一。
一般来说,研究人员主要集中在均匀线性阵列(ULA)上,在以往的研究中,互耦校准和补偿通常是通过硬件实现的,如增加校准元件和采用低耦合元件,然而,这种方法由于成本相对较高,精度较低,在许多应用中都不容易实现,所以互耦的校准和补偿逐渐转化为阵列参数估计。过去,国内外学者在相互耦合下提出了许多到达估计方法的方向,然而,许多算法需要多维搜索和多参数优化,故不能保证全局收敛性,众所周知,阵元间的感应电流在ULA中产生电磁耦合,从而改变阵元附近的磁场,单元间间距越小,相互耦合效应越强,基于稀疏模型的互质阵列广泛应用于无源位置,其阵元间间距大于半波长,阵元间距的增加使互耦效应显著减小,同时,阵列孔径也得到了有效的扩展,自由度(DOF)也得到了相应的提高,正因为如此,近几年来,互质阵列已成为阵列信号处理领域的研究热点。
互质阵列是一种特殊的非均匀稀疏阵列,相比于最小冗余阵列,互质阵列的布局设计更为简单,而相比于嵌套阵列,互质阵列的阵元间距更大,阵元间互耦影响更小。互质阵列对应的差分合成阵列中虚拟阵元并不完全均匀连续分布,这意味着存在缺失的虚拟阵元即空洞。根据布局特征,两个对称的空洞将互质阵列的差分合成阵列划分为三个子差分合成阵,一个是位于中间的均匀子差分合成阵,另外两个是位于两端的非均匀子差分合成阵,其中,均匀子差分合成阵是差分合成阵的主体部分。然而,由于其空洞的存在,其应应用中舍弃了位于分合成阵两端的两个非均匀子差分合成阵,这意味着自由度的损失,也就是说,互质阵列提供的全部自由度并未被完全利用。
发明内容
针对现有技术中的上述不足,本发明提供的一种高自由度去耦合的空时自适应信号处理方法,利用两个工作频率获得了比单个工作频率更多的均匀差分合成虚拟阵列和脉冲,带来了滤波器自由度的提升,增加了系统的检测性能。
为了达到以上目的,本发明采用的技术方案为:
本方案提供一种高自由度去耦合的空时自适应信号处理方法,其特征在于,包括以下步骤:
S1、在参考工作频率下,利用雷达系统的互质结构获取低互耦的原始杂波协方差矩阵;
S2、根据所述原始杂波协方差矩阵,通过添加额外工作频率填充参考工作频率下差分合成结构中的空洞,增加差分合成结构中虚拟阵元和虚拟脉冲的个数;
S3、根据增加的虚拟阵元和虚拟脉冲的个数,利用矩阵增秩的方法填充原始杂波协方差矩阵得到虚拟杂波协方差矩阵;
S4、根据所述虚拟杂波协方差矩阵计算得到虚拟滤波器权矢量,完成对高自由度去耦合的空时自适应信号的处理。
本发明的有益效果是:本发明通过选择合适的单一额外工作频率获取额外不同位置的差分合成虚拟阵元和差分合成虚拟脉冲,进而填充参考工作频率下互质结构差分合成阵和差分合成脉冲中两个缺失的虚拟阵元和脉冲。相比于单频工作模式,双频工作模式获取了更多均匀分布的虚拟阵元和脉冲,滤波器的自由度得到了提升,改善了信号检测性能和角度和多普勒分辨率,由此产生的互质结构可以缓解相互耦合,增强滤波器的虚拟自由度。
进一步地,所述步骤S1中原始协方差矩阵的表达式如下:
P(w)=diag([p1,p2,...,pNc]T),pi=E(|ac,i|2)
fc,i(w)=(2vT cos(θi))/λ
其中,Rc(w)表示原始协方差矩阵,V(w)表示杂波空时导向矩阵,P(w)表示杂波功率矩阵,VH(w)表示V(w)的共轭转置,RNM-1,NM-1(w)表示杂波协方差矩阵中第NM-1行第NM-1列的元素,表示杂波协方差矩阵中第l1M+r1-1行第l2M+r2-1列的元素,且l1,l2=0,…,N-1,r1,r2=1,…,M,Nc表示方位域内独立杂波块的数目,pk表示第K块杂波功率,j表示虚数,k表示波数,和nl2分别表示第和nl2个阵列位置,d表示间距,θi表示第i个杂波块的角度,和分别表示第和个脉冲位置,T表示最小脉冲重复间隔,v表示雷达平台速度,表示第Nc个杂波块的空时导向矢量,表示第Nc块杂波功率,pi表示第i个杂波块,E(·)表示求期望运算,|ac,k|2表示第K块杂波块幅度均方,和v(fc,i(w))分别表示工作频率为w时第i个杂波块的空间和时间导向矢量,sNM-1,i(w)表示第i个杂波块的空时导向矢量中第NM-1个元素,nl表示第nl个阵列位置,和fc,i分别表示第i个杂波块的空间和时间归一化频率,λ表示雷达波长。
上述进一步方案的有益效果是:本发明通过计算得到原始协议差矩阵,为添加额外工作频率提供良好的条件,以使能填充参考工作频率下的差分结构中的空洞。
再进一步地,所述步骤S2中虚拟阵元的额外工作频率的表达式如下:
其中,wh表示虚拟阵元的额外工作频率,N1和N2分别表示互质阵列中第1个和第2个子阵列的阵元个数,w0表示参考工作频率。
上述进一步方案的有益效果是:本发明中附加频率的使用在差分合成阵列中指定位置增加了额外的虚拟阵元,因此,合理选择附加频率能够产生期望的虚拟阵元出现在差分合成阵列的空洞位置。
再进一步地,所述步骤S2中虚拟脉冲的额外工作频率的表达式如下:
其中,wh2表示脉冲的额外工作频率,M1和M2分别表示互质脉冲串中第1个和第2个子脉冲串的脉冲个数,w0表示参考工作频率。
上述进一步方案的有益效果是:本发明中附加频率的使用在脉冲中指定位置增加了额外的脉冲,因此,合理选择附加频率能够产生期望的脉冲出现在差分合成脉冲的空洞位置。
再进一步地,所述步骤S3中虚拟杂波协方差矩阵的表达式如下:
D=(2N1N2+4N1-1)(2M1M2+4M1-1)
其中,Rc(w0)表示参考频率下的虚拟杂波协方差矩阵,RD-1,D-1(w)表示虚拟杂波协方差矩阵第D-1行第D-1列的元素,D表示矩阵维度,M1和M2分别表示互质脉冲串中第1个和第2个子脉冲串的脉冲个数,N1和N2分别表示互质阵列中第1个和第2个子阵列的阵元个数。
上述进一步方案的有益效果是:本发明利用矩阵增秩的方法来获得用于STAP分析的虚拟杂波协方差矩阵,能有效地避免对进行空间平滑操作带来的复杂计算量。
再进一步地,所述步骤S4中虚拟滤波器权矢量的表达式如下:
Ru(w0)=Rc(w0)+σn 2I(w0)
其中,w表示虚拟滤波器权矢量,vt表示目标空时导向矢量,vt H表示vt的共轭转置,Ru(w0)表示杂波加噪声协方差矩阵,σn 2表示噪声功率,Rc(w0)表示参考频率下的虚拟杂波协方差矩阵,I(w0)表示单位阵。
上述进一步方案的有益效果是:本发明利用求得的杂波加噪声协方差矩阵能够求出滤波器权矢量。
附图说明
图1为本发明的方法流程图。
图2为本实施例中的互质结构图。
图3为本实施例中差分合成阵列和差分合成脉冲示意图。
图4为本实施例中的自由度DOF比率图。
图5为本实施例中的无互偶时空时的波束对比图。
图6为本实施例中的有互偶时空时的波束对比图。
图7为本实施例中多普勒域和角度域的对比图。
图8为本实施例中的输出SINR示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
实施例
减少互质阵列自由度损失的关键是填充差合成阵列的空洞,进而联合非均匀差分子合成阵的虚拟阵元形成阵元数目更多、孔径更大的均匀差分子合成阵。考虑到互质阵列差分合成阵中空洞的分布特点,在参考工作频率之外,本申请采用一个额外的工作频率来扩展合成阵的布局结构,用以填充工作在参考频率时差分合成阵中均匀差分子合成阵与两个非均匀差分子合成阵之间两个对称的空洞,这样包含被填充空洞在内的2N1个虚拟阵元被吸收到均匀差分子合成阵中。因此,互质阵列通过双频工作模式形成了包含2N1N2+4N1-1个虚拟阵元的均匀差分子合成阵,获得的自由度为2N1N2+4N1-1。而对于脉冲串采用同样的工作方式,同样可以填充空洞,增加均匀虚拟脉冲,提高系统时间自由度,当阵列和脉冲参用相同的布局结构时,二者的差分结构空洞的位置相同,故其所需额外的工作频率一致。相比于单频工作模式,双频模式采用一个额外的工作频率获得的时间和空间自由度提升分别为2N1和2M1,系统总自由度提升为4N1N2M1+4N1M1M2+12N1M1-2N1-2M1。此外,为了避免复杂的空间平滑处理,本申请采用了矩阵增秩的方法来构造用于STAP滤波器权矢量估计的虚拟杂波协方差矩阵,数值仿真结果验证了本申请提出算法的有效性。如图1所示,本发明提供了一种高自由度去耦合的空时自适应信号处理方法,其实现方法如下:
S1、在参考工作频率下,利用雷达系统的互质结构获取低互耦的原始杂波协方差矩阵;
S2、根据所述原始杂波协方差矩阵,通过添加额外工作频率填充参考工作频率下差分合成结构中的空洞,增加差分合成结构中虚拟阵元和虚拟脉冲的个数;
S3、根据增加的虚拟阵元和虚拟脉冲的个数,利用矩阵增秩的方法填充原始杂波协方差矩阵得到虚拟杂波协方差矩阵;
S4、根据所述虚拟杂波协方差矩阵计算得到虚拟滤波器权矢量,完成对高自由度去耦合的空时自适应信号的处理。
本实施例中,首先,我们将给出关于互质结构的空时自适应信号处理的一些基本准备工作。假设侧视机载相控阵雷达具有N个接受阵元和M个发射脉冲,雷达平台速度为v,波长为λ,接收阵列由两个子阵列组成,一个是具有N2个位于{N1p2d,0≤p2≤N2-1}的阵元组成,另一个由具有2N1-1个位于{N2p1d,1≤p1≤2N1-1}的阵元组成位置的传感器组成,N1和N2是互质的整数且N1<N2,d=λ/2是最小阵元间距,如图2(a)所示,其为互质阵列。类似的,有具有时间间隔为{M1q2T,0≤q2≤M2-1}和{M2q1T,1≤q1≤2M1-1}的两个子发射脉冲组成的脉冲串,脉冲总个数为M,如图2(b)所示,其为互质脉冲串,这里M1和M2是互质的整数且M1<M2,T最小脉冲重复间隔,工作频率为w时,接收到的杂波加噪声数据为:
这里:
nα∈{N1p2∪N2p1},α=1,…,N-1 (5)
mβ∈{M1q2∪M2q1},β=1,…,M-1 (6)
对应的第i个杂波块的空间和时间归一化频率分别为:
fc,i(w)=(2vTcos(θi))/λ (8)
式中,θi第i个杂波块的角度。根据(7)式和(8)式,(3)式和(4)式可以重写为:
进一步:
式中:
式中,l=0,…,N-1,r=1,…,M,i=1,…,Nc。假设不同的杂波块是独立的,杂波加噪声协方差矩阵(CNCM)可以表示为:
Ru(w)=E[xu(w)xu H(w)]=Rc(w)+σn 2INM(w) (15)
式中,Rc(w)=V(w)P(w)VH(w)代表杂波协方差矩阵。其中,杂波空时导向矩阵为:
式中:
式中,l1,l2=0,…,N-1,r1,r2=1,…,M。
本实施例中,在实际雷达系统中,由于阵元之间的电磁耦合,接收信号发生了变化。根据电磁耦合原理,阵元间间距越大,耦合效应越小,如果阵元间间距大于几个波长,则可以忽略耦合效应,因此,阵列的权重函数w(m)可以描述互偶效应的大小,定义为:
本实施例中,互质阵列前三个w(m)都是2,然而均匀线性阵列的前三个权重函数分别是w(1)=N-1,w(2)=N-2,w(3)=N-3,N≥3,其中,N是阵元个数。当N>5,均匀线性阵列的权重函数大于互质阵列的,这就意味着均匀线性阵列互偶很重。
本实施例中,本申请通过双频来加强互质阵列和脉冲的配置,以扩展连续的差分合成阵列和脉冲序列,提高时间和空间自由度。
其中,连续的部分只从-(N1N2+N1-1)到(N1N2+N1-1)。如图3(b)表明,互质阵列结构的两个不连续差分子阵列完全对称,不连续差分协方差可以细分为一个单孔子阵列和几个孔子阵列。
如图3(b)表明,差分脉冲串的结构和性质与差分阵列相同。如果差分阵和脉冲中的孔可以被填充,可以显著增加可用的虚拟阵元和脉冲的数量,以在空间和时间域获得更大的DOF。图中,实心圆代表物理传感器,空心圆代表空空间。实心矩形表示物理发送脉冲,空心矩形表示空空间。
本实施例中,根据上述内容可知:互质采样结构在单频工作时,只利用了差分结构中的均匀虚拟子阵列和子脉冲,在空间域和时间域中分别提供了2N1N2+2N1-1和2M1M2+2M1-1的自由度,系统滤波器具有(2N1N2+2N1-1)(2M1M2+2M1-1)的虚拟自由度。也就是说,互质结构STAP算法的检测能力依赖于差分合成结构中均匀的虚拟阵元和脉冲的个数,舍弃非均匀差分合成子阵列和子脉冲而导致系统分辨和检测能力下降的关键是差分合成结构中空洞的存在。因此,本申请提出使用双频来填充互质采样结构中差分合成阵列和脉冲串中的空洞,扩大均匀差分合成子阵列和子脉冲的范围,从而获得更多的滤波器虚拟自由度,实现系统性能的提高。假设信号源的频率范围足以覆盖互质采样结构工作的两个频率,即参考频率w0和用于填充差分合成结构空洞的额外频率wh。考虑图2所示的互质采样结构,假设参考频率w0对应的阵元间隔为d0=λ/2和T=T0,k=k0。当互质采样结构工作在额外频率wh=αhw0时,接收信号表示为:
其中,l=0,…,N-1,r=1,…,M,i=1,…,Nc.kh=wh/c是wh时的波数。由于kh=αhk0,式(26)可改写为:
对比(14)和(27)式可知,wh对应的空时导向矩阵等效于尺度变换后w0对应的空时导向矩阵,尺度变换后第i个阵元位置为αhnid,第j个脉冲位置为αhmjT。因此,互质结构工作在wh时的差分合成阵列和脉冲也是工作在w0时差分合成阵列和脉冲对应的尺度变换,对于式(1)所示的w0时差分合成阵列和脉冲中对应的虚拟阵元和脉冲的位置,当工作于wh时尺度变换分别为:
当wh>w0时,及αh>1时,差分合成阵列和脉冲会扩展,反之则缩小。换句话说,附加频率的使用在差分合成阵列和脉冲中指定位置增加了额外的虚拟阵元和脉冲。因此,合理选择附加频率能够产生期望的虚拟阵元和脉冲出现在差分合成结构的空洞位置。对于图3(a)所示的差分合成阵列中均匀子差分合成阵两侧对称的处于±(N1N2+N1)的两个空洞,对其填充所使用的额外工作频率可选择为:
也就是说,与空洞外侧相邻的两个虚拟阵元向内侧收缩到空洞位置,从而实现了空洞填充。事实上,填充±(N1N2+N1)的两个空洞可用任意的虚拟阵元进行尺度变换来实现,考虑到额外频率使用对系统复杂度和信号源特性起伏的影响,额外频率的选择要求接近并且小于参考频率,因此,式(11)中的额外频率是填充±(N1N2+N1)的两个空洞的最佳选择,处于±(N1N2+N1)的两个空洞被填充后,获得的均匀差分合成虚拟阵元扩张到±(N1N2+2N1-1),获得的自由度达到2N1N2+4N1-1,相比于单频工作时2N1N2+2N1-1,一个额外频率的使用带来了虚拟自由度的提升。同理,对于图3(b)所示的差分合成脉冲中均匀子合成脉冲串两侧对称的处于±(M1M2+M1)的两个空洞同样可以使用上述方法填充,对其填充所使用的额外工作频率可选择为:
本实施例中,当阵列和脉冲采用相同配置方式的互质结构时,其空洞所处位置相同,换言之,添加一个额外频率,可以同时填充差分合成阵列和脉冲中的空洞,相对而言系统复杂度较低,本申请研究的就是此结构。
本实施例中,针对双频工作时获得的均匀分布于-(N1N2+2N1-1)至(N1N2+2N1-1)的均匀差分合成虚拟阵元和-(M1M2+2M1-1)至(M1M2+2M1-1)的均匀差分合成虚拟脉冲。为了避免对进行空间平滑操作带来的复杂计算量,本申请采用矩阵增秩的方法来获得用于STAP分析的虚拟杂波协方差矩阵,即表示为:
式中,D=(2N1N2+4N1-1)(2M1M2+4M1-1)。此虚拟杂波协方差矩阵等效于以阵元间距为d0阵元数量为2N1N2+4N1-1的均匀线性阵列和固定重复周期为T0的2M1M2+4M1-1个脉冲的脉冲串所对应的杂波协方差矩阵,相比于单频工作模式,互质结构使用双频工作模式将自由度提升了4N1N2M1+4N1M1M2+12N1M1-2N1-2M1。依此类推,使用多个额外的工作频率能够填充互质结构中差分合成阵列和脉冲中非均匀子差分阵和脉冲多个空洞部分的多个空洞,但此时一个额外频率带来的自由度提升量小于4N1N2M1+4N1M1M2+12N1M1-2N1-2M1。相比于本申请使用一个额外频率填充非均匀子差分阵和脉冲单个空洞部分的一个空洞,在同等代价下对自由度的提升有限,因此,本申请重点针对单个空洞部分的一个空洞设计互质结构双频工作模式。
本实施例中,所述虚拟滤波器权矢量的表达式如下:
Ru(w0)=Rc(w0)+σn 2I(w0) (34)
其中vt表示目标空时导向矢量,Ru(w0)表示杂波加噪声协方差矩阵,σn 2表示噪声功率,I(w0)表示单位阵(假设噪声功率不随工作频率发生变化)。
本实施例中,本申请选择了三种STAP法:传统的STAP(T-STAP)、传统的互质STAP(C-STAP)以及本申请所提出的算法(DFC-STAP),然后比较了它们的性能。当存在互偶时,T-STAP、C-STAP和DFC-STAP分别简写为MC-T-STAP、MC-C-STAP和MC-DFC-STAP。雷达系统中有N=6物理脉冲和M=6物理传感器,其中,λ=0.05m,T0=0.25ms,Nc=361,v=50m/s和对于互质STAP我们设置参数为N1=M1=2和N2=M2=3,假设系统的工作频率为参考频率w0,相应的阵元间间距为d0=λ/2和脉冲间距为T0,目标归一化角和多普勒频率分别为0.1和-0.2,杂噪声比为30dB,信噪比为0dB,所有模拟结果为100次结果取均值。为了直观地比较各种算法获得DOF的能力,我们将DOF比率定义为:
γ(K)=L(K)/LDFC(K) (35)
式中,LDFC(K)表示当物理传感器和脉冲的总数同时为K时,所提出的算法的最大DOF,γ(K)=L(K)/LDFC(K)与DOF的容量在一定程度上是成正比的。图4说明了C-STAP和DFC-STAP比T-STAP具有更高的DOF,然而,T-STAP的DOF明显低于DFC-STAP,特别是,在较小K的时候比其他方法有明显的增强。
本实施例中,本申请描述了不同方法的波束图,图5和图6给出了有无相互耦合时的空时波束图。令c=[1,0.5ejπ/4,0.25ej0.7π,0.5ej0.7π/3]T和B=3,由于相互耦合,传统算法不能再抑制主杂波区域的杂波,性能严重恶化。虽然传统的互质算法可以抑制主杂波并检测信号,但其空时分辨率低于本申请所提的算法,无论是否存在相互耦合,FDC-STAP都可以在目标位置形成最佳的空时波束,在主杂波区域形成一个缺口,然而,T-STAP具有最严重的相互耦合效应,由于C-STAP减少了相互耦合,与T-STAP相比,即使C-STAP的DOF小于FDC-STAP,其空间和多普勒分辨率低于FDC-STAP,如图7所示。这是由于传统互质算法的互质结构中差分合成阵和脉冲中±8的位置存在虚拟阵元和脉冲缺失,均匀的15个虚拟阵元和15个虚拟脉冲处于-7至7的区间,因此,从参考频率w0对应的空时导向矩阵中能够提取-7~7连续相关阵元和脉冲的相关元素,作为15个均匀分布虚拟阵元和脉冲的等效数据,此时可获得225系统自由度,超过了互质结构中物理阵元和脉冲所获得的原始自由度(原始算法)。本申请所提算法根据式(31)选取额外的工作频率wh=8/9w0以此填充处于±8位置的空洞,此时参考频率w0时虚拟阵元和脉冲位置±9收缩为±8,用于实现空洞的填充,因此,利用两个频率可分别获得从-9至9的19个均匀连续的虚拟阵元和脉冲,提供的自由度提升到361。本申请相对传统算法,利用一个额外的工作频率,提升了136的自由度。最后,图8给出了有无互偶下三种算法的输出SINR性能与归一化多普勒频率对比图。本申请所提由于采用双频工作模式填充了差分结构中的空洞,提高了系统的自由度,所以其输出SINR性能最优,依次为C-STAP、T-STAP,同时由于互质阵列相对均匀线性阵列阵元间距较大,互偶效应较弱,互偶条件下其SINR性能几乎不变,而T-STAP则受到了互偶的影响,恶化了输出SINR的性能。
本实施例中,本申请提出的双频互质结构STAP算法,利用两个工作频率获得了比单个工作频率更多的均匀差分合成虚拟阵列和脉冲,带来了滤波器自由度的提升,增加了系统的检测性能。对于具有N2+2N1-1个物理阵元和M2+2M1-1)个物理脉冲的互质结构,双频工作比单频工作即增加使用一个额外的工作频率能够带来4N1N2M1+4N1M1M2+12N1M1-2N1-2M1个自由度的提升。由于一个额外工作频率的使用,双频工作增加了系统的工作带宽并引入了信号源特性起伏,从而在一定程度上增加了系统的复杂程度与成本以及STAP滤波器权矢量估计误差,本申请选择低于且最接近于参考工作频率的额外工作频率,能够最大程度上减小双频工作带来的不利影响。
Claims (6)
1.一种高自由度去耦合的空时自适应信号处理方法,其特征在于,包括以下步骤:
S1、在参考工作频率下,利用雷达系统的互质结构获取低互耦的原始杂波协方差矩阵;
S2、根据所述原始杂波协方差矩阵,通过添加额外工作频率填充参考工作频率下差分合成结构中的空洞,增加差分合成结构中虚拟阵元和虚拟脉冲的个数;
S3、根据增加的虚拟阵元和虚拟脉冲的个数,利用矩阵增秩的方法填充原始杂波协方差矩阵得到虚拟杂波协方差矩阵;
S4、根据所述虚拟杂波协方差矩阵计算得到虚拟滤波器权矢量,完成对高自由度去耦合的空时自适应信号的处理。
2.根据权利要求1所述的高自由度去耦合的空时自适应信号处理方法,其特征在于,所述步骤S1中原始杂波协方差矩阵的表达式如下:
fc,i(w)=(2vTcos(θi))/λ
其中,Rc(w)表示原始杂波协方差矩阵,V(w)表示杂波空时导向矩阵,P(w)表示杂波功率矩阵,VH(w)表示V(w)的共轭转置,RNM-1,NM-1(w)表示杂波协方差矩阵中第NM-1行第NM-1列的元素,表示杂波协方差矩阵中第l1M+r1-1行第l2M+r2-1列的元素,且l1,l2=0,…,N-1,r1,r2=1,…,M,Nc表示方位域内独立杂波块的数目,pk表示第K块杂波功率,j表示虚数,k表示波数,和nl2分别表示第和nl2个阵列位置,d表示间距,θi表示第i个杂波块的角度,和分别表示第和个脉冲位置,T表示最小脉冲重复间隔,v表示雷达平台速度,表示第Nc个杂波块的空时导向矢量,表示第Nc块杂波功率,pi表示第i个杂波块,E(·)表示求期望运算,|ac,k|2表示第K块杂波块幅度均方,sNM-1,i(w)表示第i个杂波块的空时导向矢量中第NM-1个元素,nl表示第nl个阵列位置,和fc,i(w)分别表示第i个杂波块的空间和时间归一化频率,λ表示雷达波长。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010574966.1A CN111610498B (zh) | 2020-06-22 | 2020-06-22 | 一种高自由度去耦合的空时自适应信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010574966.1A CN111610498B (zh) | 2020-06-22 | 2020-06-22 | 一种高自由度去耦合的空时自适应信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111610498A true CN111610498A (zh) | 2020-09-01 |
CN111610498B CN111610498B (zh) | 2023-04-07 |
Family
ID=72195417
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010574966.1A Active CN111610498B (zh) | 2020-06-22 | 2020-06-22 | 一种高自由度去耦合的空时自适应信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111610498B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116299401A (zh) * | 2023-05-19 | 2023-06-23 | 成都航空职业技术学院 | 基于目标散射点位置的恒虚警方法、设备及其存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6400306B1 (en) * | 1999-12-17 | 2002-06-04 | Sicom Systems, Ltd | Multi-channel moving target radar detection and imaging apparatus and method |
US20040090906A1 (en) * | 2002-11-08 | 2004-05-13 | Shimon Moshavi | Reduced complexity MMSE multiuser detection for a multirate CDMA link |
CN104977571A (zh) * | 2015-06-25 | 2015-10-14 | 西安电子科技大学 | 基于俯仰频率分集stap的距离模糊杂波抑制方法 |
CN106338723A (zh) * | 2016-09-12 | 2017-01-18 | 深圳大学 | 一种基于互质脉冲重复间隔的空时自适应处理方法及装置 |
CN107329108A (zh) * | 2017-05-03 | 2017-11-07 | 浙江大学 | 基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法 |
CN108802705A (zh) * | 2018-04-24 | 2018-11-13 | 深圳大学 | 一种基于稀疏的空时自适应处理方法及系统 |
-
2020
- 2020-06-22 CN CN202010574966.1A patent/CN111610498B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6400306B1 (en) * | 1999-12-17 | 2002-06-04 | Sicom Systems, Ltd | Multi-channel moving target radar detection and imaging apparatus and method |
US20040090906A1 (en) * | 2002-11-08 | 2004-05-13 | Shimon Moshavi | Reduced complexity MMSE multiuser detection for a multirate CDMA link |
CN104977571A (zh) * | 2015-06-25 | 2015-10-14 | 西安电子科技大学 | 基于俯仰频率分集stap的距离模糊杂波抑制方法 |
CN106338723A (zh) * | 2016-09-12 | 2017-01-18 | 深圳大学 | 一种基于互质脉冲重复间隔的空时自适应处理方法及装置 |
CN107329108A (zh) * | 2017-05-03 | 2017-11-07 | 浙江大学 | 基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法 |
CN108802705A (zh) * | 2018-04-24 | 2018-11-13 | 深圳大学 | 一种基于稀疏的空时自适应处理方法及系统 |
Non-Patent Citations (5)
Title |
---|
MINGXIN LIU ET AL: "Robust STAP With Reduced Mutual Coupling and Enhanced DOF Based on Super Nested Sampling Structure", 《IEEE ACCESS》 * |
王乐宁 等: "双频机载预警雷达系统动态杂波自由度研究", 《现代雷达》 * |
阳召成 等: "基于空时功率谱稀疏性的空时自适应处理技术研究进展", 《电子学报》 * |
陈川 等: "基于双频互质阵列的DOA估计方法", 《计算机应用研究》 * |
马源斌: "互质空时压缩采样机载雷达杂波抑制技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116299401A (zh) * | 2023-05-19 | 2023-06-23 | 成都航空职业技术学院 | 基于目标散射点位置的恒虚警方法、设备及其存储介质 |
CN116299401B (zh) * | 2023-05-19 | 2023-10-17 | 成都航空职业技术学院 | 基于目标散射点位置的恒虚警方法、设备及其存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN111610498B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107329108B (zh) | 基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法 | |
CN107315160B (zh) | 基于内插虚拟阵列信号原子范数最小化的互质阵列波达方向估计方法 | |
Guerci | Theory and application of covariance matrix tapers for robust adaptive beamforming | |
CN107102291B (zh) | 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 | |
CN109471086B (zh) | 基于多采样快拍和集阵列信号离散傅里叶变换的互质mimo雷达波达方向估计方法 | |
CN109143190B (zh) | 一种零陷展宽的宽带稳健自适应波束形成方法 | |
CN104020469B (zh) | 一种mimo雷达距离-角度二维超分辨率成像算法 | |
CN101887117B (zh) | 基于三迭代的机载mimo雷达空时降维自适应处理方法 | |
CN109765521B (zh) | 一种基于子阵划分的波束域成像方法 | |
CN107104720B (zh) | 基于协方差矩阵虚拟域离散化重建的互质阵列自适应波束成形方法 | |
CN107092007A (zh) | 一种虚拟二阶阵列扩展的波达方向估计方法 | |
CN104656073B (zh) | 三维成像声纳波束形成方法及在多核处理器上的实现方法 | |
CN107219511A (zh) | 波束‑多普勒方向图稀疏约束的stap方法及装置 | |
CN111610498B (zh) | 一种高自由度去耦合的空时自适应信号处理方法 | |
CN110082789B (zh) | 基于圆阵的空时域自适应宽线性降秩波束形成方法 | |
CN108828586B (zh) | 一种基于波束域的双基地mimo雷达测角优化方法 | |
CN114185007A (zh) | 基于多核dsp的空时联合降维自适应杂波抑制方法 | |
CN109471087B (zh) | 基于互质mimo雷达差集和集信号快速傅里叶变换的波达方向估计方法 | |
CN108614234A (zh) | 基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法 | |
CN109490846B (zh) | 基于空时联合优化的多输入多输出雷达波形设计方法 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
Li et al. | Joint design of transmit beamforming and receive filter for transmit subaperturing MIMO STAP radar | |
CN110850421A (zh) | 基于混响对称谱的空时自适应处理的水下目标检测方法 | |
CN107167782B (zh) | 基于信杂噪比最大的雷达三维异构阵稀疏重构方法 | |
CN113406586B (zh) | 基于约束张量分解的mimo雷达二维波达方向估计方法 |
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 |