[发明专利]基于电磁矢量传感器的MIMO-Y型天线阵列形成方法

权利要求书

1.基于电磁矢量传感器的MIMO-Y型天线阵列形成方法，其特征在于，包括以下步骤：

步骤一、计算接收端多天线为Y型阵列时的入射信号空间导向矢量，

所述Y型阵列由Y₁、Y₂、Y₃三个子线形阵列组成，相邻子线形阵列夹角为120°，相邻阵元间距为d，假设第k个天线阵元的坐标为(x_k，y_k)(k＝1，2，...，3M)，则

xk=(k-1)dk≤M-0.5(k-M-1)dM+1≤k≤2M-0.5(k-2M-1)d2M+1≤k≤3M]]>

yk=0k≤M0.53(k-M-1)M+1≤k≤2M-0.53(k-2M-1)2M+1≤k≤3M]]>

以三维坐标中的原点为参考点，所述Y型阵列的导向矢量为：

α(θ，ψ)_UYA＝[α₁(θ，ψ)^T，α₂(θ，ψ)^T，α₃(θ，ψ)^T]^T

其中，

α1(θ,ψ)=[1,ejkωτ2,...,ejkwτM]Tα2(θ,ψ)=[1,ejkwτM+2,...,ejkwτ2M]Tα3(θ,ψ)=[1,ejkwτ2M+2,...,ejkwτ3M]T]]>

其中k_ω＝2π/λ，λ表示接收信号的波长，τ_k＝x_kcosψsinθ+y_ksinψsinθ，得到：

α(θ,ψ)UYA=[1,...,ejkwd(k-1)cosψsinθ,...,1,...,ejkwd(k-M-1)(-12cosψ+32sinψ)sinθ,...,1,...,ejkwd(k-2M-1)(-12cosψ+32sinψ)sinθ,...]T;]]>

步骤二、计算基于电磁矢量传感器的所述Y型阵列导向矢量，

电磁矢量传感器EVS的导向矢量为：

αEVS(Θ)=[eTbT]T=Ψ(θ,ψ)Ω(γ,η)=cosθcosψ-sinψcosθψsinψcosψ-sinθ0-sinψ-cosθψcosψcosψ-cosθψsinψ0sinθsinγejηcosγ]]>

引入电磁矢量传感器的所述Y型阵列导向矢量表示为：

αEYA(Θ)=αEVS(Θ)⊗αUYA(θ,ψ);]]>

步骤三、计算所述Y型阵列空间衰落相关性函数，

在MIMO多天线阵列中，阵元m和阵元n之间的空间衰落相关定义为：

其中E[·]为数学期望，(·)^*表示复数共轭，为阵元m接收信号能量均值，和分别为阵元m和n的导向矢量，为波达信号的三维空间概率分布函数，在所述Y型阵列中ρ(m，n)的实部和虚部分别为

Re[ρ(m,n)]=-1sinc(Δθ)sin(θ0)×{Σk=0∞Σm=0k(-1)2k+m+1Z2k(22kk!)22k+1m×sinc[(2k+1-2m)Δθ]sin[(2k+1-2m)θ0]+2Σk=1∞Σm=0∞Σl=0k+m(-1)k+2m+l+1m!Γ(2k+m+1)×sinc(2kΔψ)cos(2k(ψ0+ξ))×sinc[[2(k+m-l)+1]Δθ]2(k+m)+1l×sin[[2(k+m-l)+1]θ0](Z4)2(k+m)}]]>

Im[ρ(m,n)]=2sinc(Δθ)sin(θ0)Σk=0∞Σm=0∞(-1)mm!Γ(2k+m+2)×sinc[(2k+1)Δψ]×sin[(2k+1)(ψ0+ξ)](Z2)2(k+m)+1×{122(k+m+1)2(k+m+1)k+m+1+Σl=0k+m(-1)k+m+l+122(k+m)+1×sinc[2(k+m-l+1)Δθ]cos[2(k+m-l+1)θ0]2(k+m+1)l}]]>

其中，sinc(x)＝sin(x)/x，Z=Zx2+Zy2,ξ=tan-1ZxZy,]]>且

Z_x＝k_w(x_m-x_n)，Z_y＝k_w(y_m-y_n)；

步骤四、计算所述Y型阵列的阵元为电磁矢量传感器单元的空间衰落相关性，

所述Y型阵列中第m个电磁矢量传感器的第p个空间极化分量与第n个电磁矢量传感器的第q个空间极化分量的空间衰落相关函数为

ρ(m,n,p,q)=γ(m,n,p,q)ΦpΦq=∫ψ0-Δψψ0+Δψ∫θ0-Δθθ0+Δθ(Ψ1,pΨ1,q+Ψ2,pΨ2,q)ejkwsinθ[(xm-xn)cosψ+(ym-yn)sinψ]sinθdθdψΦpΦq,]]>

其中ΦpΦqT=∫ψ0-Δψψ0+Δψ∫θ0-Δθθ0+Δθ(Σi=12Ψi,p2(θ,ψ))sinθdθdψ∫ψ0-Δψψ0+Δψ∫θ0-Δθθ0+Δθ(Σi=12Ψi,q2(θ,ψ))sinθdθdψ,]]>

化简可得

Re1[γ(m,n,p,q)]=Σk=0∞(-1)k(k!)2(z2)2k(Ra1a200Sb1b2100+Rc1c200Sd1d2100)]]>

Re2[γ(m,n,p,q)]=2Σk=1∞Σm=0∞(-1)mm!(2k+m+2)(z2)2(k+m)×(Ra1a210Sb1b2010+Rc1c210Sc1c2010)]]>

Re[γ(m，n，p，q)]＝Re₁[γ(m，n，p，q)]+R₂[γ(m，n，p，q)]

Im[γ(m,n,p,q)]=2Σk=0∞Σm=0∞(-1)mm!Γ(2k+m+2)(z2)2(k+m)+1×(Ra1a201Sb1b2001+Rc1c201Sd1d2001)]]>

其中，Rtuvw=∫ψ0+ξ-Δψψ0+ξ+Δψsost(β)sinu(β)cosv(2kβ)cosw[(2k+1)β]dβ]]>

Stuvwx=∫θ0-Δθθ0+Δθcost(θ)sinu(θ)sinv(2k+1)(θ)sinw[2(k+l)+1](θ)×sinx[2(k+l+1)](θ)dθ,]]>

而cos(2nx)=1-4n22!sin2x+4n2(4n2-22)4!sin4x+4n2(4n2-2)(4n2-42)6!sin6x+···,]]>

sin(2nx)=2ncosx{sinx-4n2-223!sin3x+(4n2-22)(4n2-42)5!sin5x-···}.]]>

2.根据权利要求1所述的基于电磁矢量传感器的MIMO-Y型天线阵列形成方法，其特征在于：

所述步骤四中得到ρ(m，η，p，q)的过程为，将所述步骤三中的ρ(m，η)采用极坐标形式表达

ρ(p,q)=∫ψ∫θ∫γ∫ηαq(Θ)αp*(Θ)sin(θ)p(Θ)dηdγdθdψ∫ψ∫θ∫γ∫η|αp(Θ)|2sin(θ)p(Θ)dηdγdθdψ×1∫ψ∫θ∫γ∫η|αq(Θ)|2sin(θ)p(Θ)dηdγdθdψ,]]>

则αEYA(m,p)(Θ)=αEVSp(Θ)⊗αUYAm(θ,ψ)]]>

=[Ψ1,psinγejn+Ψ2,pconγ]αUYAm(θ,ψ)]]>

其中，α_UYAm(θ，ψ)是所述Y型阵列第m个阵元的导向矢量，Ψ_1，p和Ψ_2，p分别为中Ψ(θ，ψ)的第一列和第二列的第p个元素；假设p(Θ)对于的四个参数θ、ψ、η、γ是相互独立的，即p(Θ)＝p(θ)p(ψ)p(η)p(γ)且p(η)在[-π，π]上服从均匀分布，p(γ)在[0，π/2]上服从均匀分布，即得到所述步骤四中的ρ(m，η，p，q)。