基于杂波协方差矩阵特征向量分析STAP降维方法

基于杂波协方差矩阵特征向量分析STAP降维方法

基于杂波协方差矩阵特征向量分析STAP降维方法

本文在对杂波协方差矩阵的特征值、阵特征向量做出分析的基础上，讨论了无信噪比损失的降维的充要条件，并给出了一种满足该条件的降维变换阵的构造方法.最后基于一组Mountaintop数据对文中提出的算法作出了验证.　　关键词:机载雷达；时空二维自适应信号处理；杂波抑制

一、引　　言　　七十年代初Reed、Mallett和Brennan提出了时空二维自适应信号处理中的样本协方差矩阵求逆算法(SMI)［1］，并对该方法的性能做了理论分析.受实际条件限制，上述算法不可能直接用于相干脉冲处理间隔(CPI)内的全维时空快拍，对实际应用更有意义的是降维STAP算法.目前较为常用的降维措施无外乎两类.第一类是以R.Klemm提出的辅助通道法［2］，保铮教授等人提出的mDT-SAP［3］，和Hong.W提出的JDL-STAP［4］为典型，在杂波谱图上直接选出时空波束(一般为二维可分离付氏基)进行降维的方法.这类方法直观且较为可靠.第二类方法是基于杂波协方差矩阵特征值的大小选取相应的降维波束.这两种降维措施在实际应用中往往是脱节的.

二、信号模型与最优自适应权　　设一均匀线阵有L阵元，其在某一被检测距离环上得到的时空快拍为X=(XT1，XT2，…，XTL)T，其中Xi=(X1i,X2i,…,XMi)T，i=1,2,…,L是第i阵元上接收到的M个脉冲构成的向量，上标‘T’表示矩阵或向量的转置.X可以分为以下两部分X=S+N.其中S是目标，N是由地(海)杂波、接收机噪声及有源干扰构成的.此时杂波协方差矩阵R=COV(N,N),考虑到接收机噪声的存在性，可以认为杂波协方差矩阵R是正定的.在输出信杂噪比最大的准则下，最优权Wopt=kR-1S,最优输出信杂噪比SINRopt=SHR-1S，k为一常数.杂波协方差矩阵还可以表示为R=UΛUH,U=(U1,U2,…,ULM),Λ=diag(λ1,λ2,…,λLM)，其中U1,U2,…,ULM是R的标准化特征向量，λ1,λ2,…,λLM是相应的特征根.此时，最优权及最优输出信杂噪比可分别表示为

其中，上标‘H’表示矩阵或向量的共轭转置.把U1，U2，…，ULM看做是彼此正交的信号通道，则λ1,λ2,…，λLM是各通道中噪声功率，而SINRopt是各通道中信杂噪比之和.

三、基于杂波协方差矩阵特征分解的降维STAP　　设B是一个LM×K的列满秩降维变换矩阵，Z=BHX，变换后杂波协方差矩阵为RD=E(ZZH)=BHRB.当降维变换后的权W=k(BHRB)-1BHS(其中k为常数)，Z中相应的最大输出信杂噪比为SHB(BHRB)-1BHS.现在有这样几个问题：1.SHB(BHRB)-1BHS与SHR-1S谁大谁小?2.它们两者可否相等?3.在什么条件下它们两者相等?　　对于确定的K,B，易证RD=BHRB是K×K正定Hermitian阵.设γi、Vi，i=1,2,…，K是RD的特征根及相应的归一化的特征向量.则RD=BHRB=VΓVH，其中Γ=diag(γ1，γ2，…,γK),V(V1,V2,…，VK).　　进一步，RD=VΓVH=BHRB=BHUΛ1/2Λ1/2UHB.从而有Γ-1/2VHBHUΛ1/2Λ1/2UHBVΓ-1/2=IK×K.令Δ=Γ-1/2VHBHUΛ1/2,则ΔΔH=IK×K.此时SHB(BHRB)-1BHS=SHBVΓ-1/2Γ-1/2VHBHS=SHUΛ-1/2ΔHΔΛ-1/2UHS.值得注意的是Λ-1/2UHS=(λ-1/21UH1S，λ-1/22UH2S，…，λ-1/2LMUHLMS)T.则SHB(BHRB)-1BHS=(λ-1/21SHU1,λ-1/22SHU2,…，λ-1/2LMSHULM)ΔHΔ(λ-1/21UH1S,λ-1/22UH2S,…,λ-1/2LMUHLMS)T　(2)由于ΔHΔΔHΔ=ΔHΔ，可见ΔHΔ是幂等的.因而，当且仅当Λ-1/2UHS属于ΔHΔ的象空间，即Λ-1/2UHS∈Im(ΔHΔ)时

SHB(BHRB)-1BHS=SHR-1S=∑LMi=1λ-1i|UHiS|2　(3)

时，V=Λ-1/2UHS∈Im(ΔHΔ).其中O1是(K-1)×1零矩阵，O2是(LM-K+1)×(K-1)零矩阵，

D1=diag(λ-11UH1S,λ-12UH2S,…，λ-1K-1UHK-1S),D2=(λ-1KUHKS,λ-1K+1UHK+1S,…，λ-1LMUHLMS)T　(5)

证明：v∈Im(ΔHΔ)当且仅当　(B(BHRB)-1BH-R-1)S=O.　(6)事实上，如果(B(BHRB)-1BH-R-1)S=O，则SHB(BHRB)-1BHS=SHR-1S.若v∈Im(ΔHΔ)，则ΔHΔv=v.于是UΛ-1/2ΔHΔΛ-1/2UHS=UΛ-1/2Λ-1/2UHS，既(B(BHRB)-1BH-R-1)S=O.由于

其中Λ1=diag(λ1，λ2，…,λK-1)，Λ2=diag(λK,λK+1,…,λLM),则

其中(P1,P2)T=UHS,且P1=(UH1S,…,UHK-1S)T,P2=(UHKS,…,UHLMS)T.由于D2=Λ-12P2，则

注释：从证明的过程中可以看出，降维变换矩阵B具有更为一般的形式.设∑m1i=1λ-1i|UHiS|2＞0,∑m2i=m1+1λ-1i|UHiS|2＞0,…,∑LMi=mK-1+1λ-1i|UHiS|2＞0,B=U(diag(D1,D2,…,DK))，其中D1=(λ-11UH1S,λ-12UH2S,…，λ-1mUHmS)T，D2=(λ-1m1+1UHm1+1S,λ-1m1+2UHm1+2S,…,λ-1m2UHm2S)T，…，DK=(λ-1mK-1+1UHmK-1+1S,λ-1mK-1+2UHmK-1+2S,…,λ-1LMUHLMS)T，　(10)则v∈Im(ΔHΔ).　　从降维变换矩阵B的构造可以看出，B是根据各Ui(i=1,2,…,LM)通道中信杂噪比的大小对Ui做线性组合后得到的，其降维变换本身就是杂波抑制的过程.　　再者，降维后的系统自由度K可取1到LM之间任意的整数.当K=1时，B=Wopt.

四、杂波协方差矩阵未知时降维自适应杂波抑制　　设X是被检测距离门上的时空快拍(基本数据)，X1,X2,…,Xm是与被检测距离门相邻的若干个距离门上的时空快拍(二次数据)，则=(1/m)∑mi=1XiXHi是杂波协方差矩阵的极大似然估计，k-1S是全维SMI最优权的估计.为达到实时处理，可采用以下两种降维STAP算法.　　(1)先利用杂波模型，根据雷达天线的发射方向图、主波束方位角及载机速度等参数计算出理想杂波协方差矩阵R′.再对R′做特征分解，并把结果带入式(4)得到降维变换矩阵B.将B作用于基本数据、二次数据，并在此基础上做自适应处理.　　(2)设Z1，Z2，…,ZQ为雷达扫描在同一方位上前数个CPI中与被检测距离门相邻的若干个距离门上的时空快拍(在此称为辅助数据).一定条件下可假设，Z1，Z2，…,ZQ与X1，X2，…,Xm相互独立，且Zi、Xj具有相同的分布.设1=(1/Q)∑Qi=1ZiZHi,则1也是杂波协方差矩阵的极大似然估计.对1做特征分解，并把结果带入式(4)得到降维变化矩阵B.将B作用于基本数据、二次数据，并在此基础上做自适应处理.　　构造上述两种降维变换矩阵B的关键是通过一定的先验知识获取对杂波协方差矩阵的估计.尽管仍需做全维的杂波协方差矩阵特征分解，但是这样做可争取到实时处理中的处理时间.分析表明，如以代替1做上述第(2)种降维自适应处理，其不仅在运算量上与全维SMI的相当，它们在性能上也是一致的［5］.

五、实测数据中的验证　　下面以Mountaintop data t38pre01v1为例验证上述算法.Mountaintop实验中，天线是由14列阵元构成的等效均匀线阵.波束宽度为26°(半功率点).每个CPI有16个脉冲.等效平台运动速度为94m/s.载波波长为0.6897m，在152公里(约在140号距离门前后)处有一目标.本文以十四个阵元的Chebyshev加权线阵做发射阵计算理想杂波协方差矩阵.降维过程中取K=30.　　图1中所示为对403个距离门做全维自适应匹配滤波得到的信噪比输出(-30dB对角加载).图2中所示为对403个距离门做降维自适应匹配滤波得到的信噪比输出.计算中还发现前一算法对对角加载很敏感，而后者对对角加载不敏感.也就是说后者的降维波束很稳定.