1 引言
在窄带多径环境下,由于多径干扰会对消期望信号,传统波束形成算法性能将急剧下降甚至失效[1]。为解决此问题,多径干扰波束形成算法得到广泛地研究,目前,根据对多径干扰的处理方式不同,可分为多径干扰抑制波束形成算法和多径信号接收波束形成算法。
多径干扰抑制波束形成算法主要是通过解相干处理或线性约束实现期望信号的有效接收。一类典型的解相干处理算法是空间平滑算法[2, 3, 4],该算法可以有效地实现解相干处理,但是牺牲了阵列的有效孔径,且稳健性较差。此外,学者还提出了如:多零陷约束算法[5]、时域加权主成分分析算法[6]以及结合最差性能约束、软约束和Duvall结构的稳健算法等[7, 8]。
多径信号接收波束形成算法[9, 10, 11, 12]对多径信号联合接收,目的是为充分利用多径信号信息。文献[11]首先利用变换矩阵去除阵列接收数据中非相关干扰信息,并估计多径信号的合成导向矢量,然后基于合成导向矢量和阵列协方差矩阵的特征结构进行波束形成,该算法有效实现了多径信号的接收,且具有较快的收敛速度,但需估计非相关干扰的信号来向。文献[12]通过在期望信号和多径干扰来向附近施加最差性能约束实现了多径信号的联合接收,且具有较好的稳健性,但该算法需估计多径干扰来向,且阵列输出性能易受多径信号初始相位差的影响。
本文首先分析多径环境下期望信号对消的原因:对应期望信号和多径干扰的阵列输出间的相位差始终在 π附近变化。然后基于反对角单位阵构造新的数据协方差矩阵和导向矢量,利用最小方差无畸变准则进行波束输出,该算法可以有效地补偿期望信号和多径干扰对应的阵列输出间的相位差,并可通过选取合适的阵元个数获取更优的阵列输出性能。最后,仿真实验分析验证了算法的有效性。
2 信号模型和Capon波束形成器
考虑P个远场窄带平面波si(t)(i=0,1,⋯,P−1)以角度θi入射到M阵元均匀直线阵。假定s0(t)是期望信号,其余P−1个信号与s0(t)非相关,且s0(t)存在Q个多径干扰,其来波方向为ˆθq(q=1,2,⋯,Q),每个多径干扰由期望信号经幅度衰减和相位延迟得到,即:
其中
ˆsq(t)表示第
q个多径信号,
ρq表示第
q个多径信号相对期望信号的幅度衰减因子,
ˆϕq为多径传播引起的相位差。
假定阵元间距为d,波达方向为θ的信号的导向矢量a(θ)可表示为:
a(θ)=[1e−j2πdλsinθ⋯e−j(M−1)2πdλsinθ]T
|
(2)
|
其中
λ为信号波长。
t时刻阵列接收数据矢量为:
x(t)=[a(θ0)+Q∑q=1a(ˆθq)ρqe−jˆϕq]s0(t)+P−1∑i=1a(θi)si(t)+n(t)
|
(3)
|
其中
n(t)为阵列接收噪声矢量。
阵列输出可表示为:
其中
w为复加权矢量。
传统Capon波束形成器[13]在满足阵列输出功率最小化的同时约束期望方向阵列响应为1,可表示为如式(5)优化问题:
minwwHRws.t.wHa(θ0)=1}
|
(5)
|
其中
R是
M×M维协方差矩阵。实际应用中,
R用采样协方差矩阵代替
其中
N为采样点数。利用拉格朗日乘子法,上述优化问题的解为:
wopt=ˆR−1a(θ0)aH(θ0)ˆR−1a(θ0)=μˆR−1a(θ0)
|
(7)
|
其中
在理想条件下,传统Capon波束形成器能够有效地抑制非相关干扰。但当多径干扰存在时,多径干扰会引起期望信号对消,波束形成器性能将急剧下降,甚至失效。下面给出期望信号对消的具体原因。
3 多径环境下期望信号对消原因
假定M=2,P=1,Q=1,期望信号为复指数信号,其离散形式为:
s0(n)=Aej2πftn=Aejϕ0(n)
|
(8)
|
其中A为期望信号幅度,tn为第n个采样点对应的采样时刻。则多径干扰可表示为:
ˆs1(n)=ρ1e−jˆϕ1s0(n)=ρ1Aej(ϕ0(n)−ˆϕ1)=ρ1Aejϕ1(n)
|
(9)
|
期望信号和多径干扰信号的导向矢量a(θ0)和a(ˆθ1)分别为:
a(θ0)=[1e−j2πdλsinθ0]T=[1ejϕθ0]T
|
(10)
|
a(ˆθ1)=[1e−j2πdλsinˆθ1]T=[1ejϕˆθ1]T
|
(11)
|
忽略噪声项n(t)的影响,阵列接收数据矢量可表示为:
x(n)=A[ejϕ0(n)+ρ1ejϕ1(n)ej(ϕ0(n)+ϕθ0)+ρ1ej(ϕ1(n)+ϕˆθ1)]T
|
(12)
|
根据式(6)和式(12)可得阵列接收数据协方差矩阵
R=A2[1+ρ21+ρ1ejˆϕ1+ρ1e−jˆϕ1e−jϕθ0+ρ21e−jϕˆθ1+ρ1ej(ˆϕ1−ϕˆθ1)+ρ1e−j(ˆϕ1+ϕθ0)ejϕθ0+ρ21ejϕˆθ1+ρ1e−j(ˆϕ1−ϕˆθ1)+ρ1ej(ˆϕ1+ϕθ0)1+ρ21+ρ1ej(ˆϕ1+ϕθ0−ϕˆθ1)+ρ1ej(ϕˆθ1−ˆϕ1−ϕθ0)]
|
(13)
|
将式(10)和式(13)代入式(7)得最优加权矢量
wopt=μA2|R|[ρ21−ρ21ej(ϕθ0−ϕˆθ1)−ρ1e−jˆϕ1+ρ1ej(ϕˆθ1−ˆϕ1−ϕθ0)ρ21ejϕθ0−ρ21ejϕˆθ1+ρ1ej(ϕθ0−ˆϕ1)+ρ1ej(ϕˆθ1−ˆϕ1)]
|
(14)
|
其中
|R|为
R的行列式值。根据式(4),对应期望信号的阵列输出为:
y0(n)=wHopta(θ0)s0(n)=μA3ρ21ejϕ0(n)|R|[2−ej(ϕθ0−ϕˆθ1)−ej(ϕˆθ1−ϕθ0)]
|
(15)
|
对应多径干扰的阵列输出为:
y1(n)=wHopta(ˆθ1)ˆs1(n)=−μA3ρ21ejϕ0(n)|R|[2−ej(ϕθ0−ϕˆθ1)−ej(ϕˆθ1−ϕθ0)]=y0(n)e±jπ
|
(16)
|
对比式(15)和式(16)可知,y0(n)和y1(n)的幅值相等,相位差为π。阵列的输出为:
此时,期望信号被完全对消,因此Capon波束形成器的性能将急剧下降。上述分析是基于2个阵元,仅存在1个多径干扰条件下进行,随着阵元个数干扰数目的增加,我们发现阵列输出对应期望信号和多径干扰间的相位差仍在
π附近变化,由于推导过程较为复杂,本文在第5节给出了仿真验证。对窄带随机信号的理论推导见附录。
4 基于反对角单位阵的多径信号接收算法
基于上述分析,若能求取最优加权矢量使对应期望信号和多径干扰的阵列输出间的相位差不等于π,则可避免期望信号的完全对消,且随着相位差的逐渐减小,接收信号的幅值逐渐增大。
构造反对角单位阵
J=[0⋯01⋮⋱⋱00⋱⋱⋮10⋯0]M×N
|
(18)
|
使
易知
J=JT=JH=J−1,因此新的阵列接收协方差矩阵为:
ˆRJ=E{ˆx(n)ˆxH(n)}=E{Jx(n)xH(n)J}=JRJ
|
(20)
|
式(2)左侧乘以J可得新的阵列导向矢量为:
aJ(θ)=Ja(θ)=e−j(M−1)2πdλsinθ[1,ej2πdλsinθ,⋯,ej(M−1)2πdλsinθ]T=e−j(M−1)2πdλsinθa∗(θ)
|
(21)
|
由式(21)可知,aJ(θ)和a∗(θ)间仅存在一个相位差,因此可认为aJ(θ)与a∗(θ)等价。将式(20)和式(21)代入式(7)可得最优加权矢量为:
wopt_J=ˆR−1JaJ(θ0)aHJ(θ0)ˆR−1JaJ(θ0)=JR−1a(θ0)aH(θ0)R−1a(θ0)=Jwopt
|
(22)
|
为使输出波束的波峰和零陷与原输入信号的方向对应,可以选择如下两种方式进行波束形成:
|Pout(θ)|=|wHopt_JaJ(θ)|
|
(23)
|
|Pout(θ)|=|wTopt_Ja∗J(θ)|=|wTopt_Ja(θ)|
|
(24)
|
其中
|Pout(θ)|为
θ方向输出波束的模值。若按式(23)进行波束形成,阵列输出为:
yout(n)=wHopt_Jˆx(n)=wHoptJHJx(n)=wHoptx(n)
|
(25)
|
此时,阵列输出与传统Capon波束形成器输出相同,期望信号将被完全对消。根据式(24),可得
yout(n)=wTopt_Jx(n)=wToptJx(n)
|
(26)
|
将x0(n)=a(θ0)s0(n)和x1(n)=a(ˆθ1)ˆs1(n)代入式(26)可得对应期望信号的多径干扰的阵列输出分别为:
ˆy0(n)=wToptJx0(n)=μA3ρ21ej(ϕ0(n)−ϕˆθ1±π)|R|(ejϕθ0−ejϕˆθ1)2
|
(27)
|
ˆy1(n)=wToptJx1(n)=μA3ρ21ej(2ϕ1(n)−ϕ0(n)−ϕθ0)|R|(ejϕθ0−ejϕˆθ1)2
|
(28)
|
对比式(27)和式(28)可得ˆy0(n)和ˆy1(n)的相位差为:
与传统Capon算法对比,上述算法在对应期望信号和多径干扰的阵列输出间引入相位差Δˆϕ=2ˆϕ1+ϕθ0−ϕˆθ1,Δˆϕ与2kπ,k∈Z的差值越大,阵列输出性能越好。上述分析基于两个阵元和一个多径干扰,其中ˆϕ1为多径干扰和期望信号间的初始相位差,当多径干扰数目增加时,ˆϕ1将随着变化,而ϕθ0−ϕˆθ1由多径信号入射到阵列的角度及阵列元个数决定,因此,Δˆϕ随着阵元个数和多径干扰数目的变化而变化,将Δˆϕ重新表述为:
其中
Δˆϕp表示路径传播引起的相位差,当入射信号的传播路径固定时,
Δˆϕp为一固定值,
Δˆϕθ表示阵列变化引起的相位差。传统多径信号接收算法的阵列输出主要由
Δˆϕp决定,当
Δˆϕp接近
2kπ,k∈Z时,阵列输出性能将急剧下降。但是,本文算法中
Δˆϕθ随着阵列信息的变化而变化,实际应用中,可以调整选取的阵元数目以获取相对传统多径信号接收算法较优的阵列输出性能。基于上述分析,基于反对角单位阵的多径信号接收算法可表示为如下优化问题:
minwwHˆRJws.t.wHa∗(θ0)=1}
|
(31)
|
并基于式(26)求取阵列输出。
5 仿真实验
实验1 阵列输出性能随阵元个数变化性能分析
考虑一均匀直线阵,期望信号来向为15º,信噪声比(SNR)为10 dB。4个非相关干扰分别以角度-20º,30º,-55º和60º入射到阵列,干扰与噪声比(INR)为30 dB。两个多径干扰相对期望信号的幅度衰减因子ρ1=ρ2=0.5,相位延迟在[0,2π]内随机产生,来向分别为-35º和45º。以第1个阵元为参考阵元,阵元间距为期望信号最高频率对应波长的一半,快拍数为512。每次试验均做200次蒙特卡罗试验求平均值。
图 1给出了阵元个数由10变化到64时,传统Capon算法和本文算法输出信号与干扰噪声比(SINR)的变化曲线。由图 1可知,本文算法的输出SINR远大于传统Capon算法,且随着阵元个数的增加在20 dB上下波动,与文中分析的阵元个数对阵列输出性能的影响一致。其中SINR由式(32)给出
SINRoutput=wTopt_JRs+ciw∗opt_JwTopt_JRui+nw∗opt_J
|
(32)
|
其中
Rs+ci表示期望信号和多径干扰对应的阵列接收数据协方差,
Rui+n为非相关干扰和噪声对应的阵列接收数据协方差。
图 2(a)和图 2(b)分别给出了阵元个数分别为28和29时,传统Capon算法和本文算法输出信号实部与期望信号的对比图。对比图 1和图 2,当阵元个数为29时,阵列输出性能和输出信号的实部值均低于阵元个数为28时,但是输出信号实部值均大于期望信号,而传统Capon算法无法输出有效的信号。由图 1和图 2可知,阵列输出性能变化趋势随着阵元个数的变化存在不规律的周期性,因此在实际应用中,可在满阵元个数的基础上,减少进行波束形成的阵元个数,以对比获取较优的阵列输出性能。
实验 2 阵列输出SINR随采样点数变化性能分析
设定阵元个数为28,改变采样点数由10到1000,其余仿真条件同实验1。图 3给出了不同采样点数下本文算法、文献[12]算法、空间平滑算法和传统Capon算法的输出SINR变化曲线。由图 3可知,本文算法的输出SINR随采样点数的增加较快达到收敛值,收敛速度大于文献[12]算法和传统Capon算法,但是略慢于空间平滑算法。
实验 3 输出SINR随输入SNR的变化性能分析
设定阵元个数为28,改变输入SNR由-10 dB到20 dB,其余仿真条件同实验1。图 4给出了不同输入SNR条件下本文算法、文献[12]算法、空间平滑算法和传统Capon算法的输出SINR变化曲线。由图 4可知,本文算法输出SINR优于文献[12]算法、空间平滑算法和传统Capon算法,且差值随着输入SNR的增加而增大。
实验 4 对应期望信号和多径干扰的阵列输出间的相位差变化分析
表 1给出了不同阵元个数和入射角度下本文算法、文献[12]算法和传统Capon算法对应期望信号和多径干扰的阵列输出间的相位差,其中初始相位差为多径传播路径延迟引起。由表 1中数据可知,传统Capon算法的输出相位差始终在π附近变化,进一步验证了多径环境下传统Capon算法中期望信号对消的原因;文献[12]算法的输出相位差与初始相位差一致,阵列输出性能将多径传播引起的相位差决定,当初始相位差接近π时,阵列输出性能将急剧下降;本文算法的输出相位差不仅受初始相位差影响,而且受阵元个数和多径信号入射角度的影响,因此,当入射到阵列的多径信号确定时,可以通过调整选取的阵元的个数修正阵列输出间的相位差,从而提高阵列输出性能。
表 1(Tab.1)
表 1 对应期望信号和多径干扰的阵列输出间的相位差Tab.1 The phase difference of the array output corresponding to the SOI and multipath interferences
DOAs | 阵元个数 | 24 | 25 | 5º | 25º | -10º | 15º | 5º | 25º | -10º | 15º | 初始相位差 | 2.6928 | -0.5984 | 2.6928 | -0.5984 | 2.6928 | -0.5984 | 2.6928 | -0.5984 | 传统Capon算法 | -3.1408 | 3.1410 | 3.1406 | -3.1408 | -3.1411 | 3.1410 | 3.1414 | -3.1411 | 文献[12]算法 | 2.6928 | -0.5984 | 2.6928 | -0.5984 | 2.6928 | -0.5984 | 2.6928 | -0.5984 | 本文算法 | 0.1956 | -0.1022 | 0.5897 | 0.2886 | 1.1996 | 0.9012 | 1.8828 | 1.5829 |
|
表 1 对应期望信号和多径干扰的阵列输出间的相位差
Tab.1 The phase difference of the array output corresponding to the SOI and multipath interferences |
6 结束语
本文针对多径环境下传统Capon波束形成算法中期望信号对消的问题,分析了期望信号对消的原因,并提出一种基于反对角单位阵的多径信号接收算法。理论分析和仿真实验表明,该算法能够有效地避免期望信号的对消,且能够获得优于传统多径信号接收算法的输出性能。
附录
窄带随机信号的复指数形式可表示为:
其中
M(t)为
s(t)的复包络,
ejw0t为复载频。其多径形式可表示为:
s(t−τ)=ρM(t−τ)ejw0(t−τ)
|
(A-2)
|
结合式(10)和式(11)可得阵列接收数据为:
x(t)=[M(t)ejw0(t)+ρM(t−τ)ejw0(t−τ)M(t)ejw0(t)ejϕθ0+ρM(t−τ)ejw0(t−τ)ejϕˆθ1]T
|
(A-3)
|
则
t时刻阵列接收数据协方差矩阵为:
R(t)=[M2(t)+ρ2M2(t−τ)+ρM(t)M∗(t−τ)ejw0τ+ρM(t−τ)M∗(t)e−jw0τM2(t)e−jϕθ0+ρ2M2(t−τ)e−jϕˆθ1+ρM(t)M∗(t−τ)ej(w0τ−ϕˆθ1)+ρM(t−τ)M∗(t)e−j(w0τ+ϕθ0)M2(t)+ρ2M2(t−τ)+ρM(t)M∗(t−τ)ej(w0τ+ϕθ0−ϕˆθ1)+ρM(t−τ)M∗(t)e−j(w0τ+ϕθ0−ϕˆθ1)M2(t)ejϕθ0+ρ2M2(t−τ)ejϕˆθ1+ρM(t−τ)M∗(t)ej(w0τ−ϕˆθ1)+ρM(t)M∗(t−τ)e−j(w0τ+ϕθ0)]
|
(A-4)
|
记
E[M2(t)]=E[M2(t−τ)]=aE[M(t)M∗(t−τ)]=E[M(t−τ)M∗(t)]=b}
|
(A-5)
|
故
ˆR=E[R(t)]=[a+ρ2a+ρbejw0τ+ρbe−jw0τa+ρ2a+ρbej(w0τ+ϕθ0−ϕˆθ1)+ρbe−j(w0τ+ϕθ0−ϕˆθ1)ae−jϕθ0+ρ2ae−jϕˆθ1+ρbej(w0τ−ϕˆθ1)+ρbe−j(w0τ+ϕθ0)aejϕθ0+ρ2aejϕˆθ1+ρbej(w0τ−ϕˆθ1)+ρbe−j(w0τ+ϕθ0)]
|
(A-6)
|
对比式(13)和式(A-6)可知,式(A-6)各项的常数和复指数项与式(13)具有相同的形式,又因为最优权矢量仅与R和期望信号导向矢量有关,故阵列输出对应期望信号和多径干扰间的相位差为±π。