[发明专利]非共点电磁矢量传感器阵列的参数估计四元数方法

摘要

非共点电磁矢量传感器阵列的参数估计四元数方法，步骤如下对接收阵列的接收数据进行采样得到采样数据，将采样数据分别按照同阵元的电偶极子和磁偶极子的采样数据叠加，构成四元数数据；计算全阵列接收四元数数据的自相关矩阵，并进行四元数特征分解，得到阵列导向矢量的估计值；由上层子阵的阵列导向矢量估计值得到上层子阵的相位中心的相对于坐标原点的相位差构成的空域导向矢量估计值，根据阵元位置矩阵与相位矩阵的关系，计算信号到达角的估计值；由上层子阵的阵列导向矢量估计值重构电偶极子子阵的阵列导向矢量估计值和磁偶极子子阵的阵列导向矢量估计值，根据子阵的阵列导向矢量间的旋转不变关系得到辅助极化角估计值和极化相位差估计值。

主权项

非共点电磁矢量传感器阵列的参数估计四元数方法，其特征在于，包括以下步骤：K个完全极化单位功率电磁波信号入射到接收阵列上，阵列的阵元由一个电偶极子和一个磁偶极子构成，步骤一、对接收阵列的接收数据进行M次采样得到采样数据，采样数据包括上层电偶极子子阵的采样数据Xue、下层电偶极子子阵的采样数据Xle、上层磁偶极子子阵的采样数据Xuh以及下层磁偶极子子阵的采样数据Xlh；步骤二：将四组采样数据分别按照同阵元的电偶极子和磁偶极子的采样数据叠加，构成四元数数据，由上层子阵和下层子阵的两组四元数数据构成全阵列接收四元数数据；第k个完全极化单位功率电磁波信号入射到第n个阵元时，该阵元的电偶极子的接收数据x'ek(n)＝ekqn(θk,φk)qe(θk)+ne(n)，该阵元的磁偶极子的接收数据x'hk(n)＝hkqn(θk,φk)qh(θk)+nh(n)，将第n个阵元的电偶极子和磁偶极子的接收数据x'ek(n)和x'hk(n)叠加，得到该阵元的四元数电磁场数据是第n个阵元接收的噪声的四元数数据；其中，ek为坐标原点处的电场，qn(θk,φk)为第n个阵元的阵元中心相对于坐标原点的相位差，qe(θk)为电偶极子相对于其阵元中心的相位，hk为坐标原点处的磁场，qh(θk)为磁偶极子相对于其阵元中心的相位，ne(n)为第n个阵元的电偶极子接收的高斯白噪声，nh(n)是第n个阵元的磁偶极子接收的高斯白噪声，e'k(n)为第n个阵元的电偶极子接收的z轴方向的电场，h'k(n)第n个阵元的磁偶极子接收的z轴方向的磁场，ck为坐标原点处z轴方向的电场和磁场叠加形成的电磁场的四元数表示；将2N个阵元的接收数据叠加成四元数数据，下层电偶极子子阵和下层磁偶极子子阵的采样数据叠加构成下层子阵的四元数数据矩阵Xl＝AlS+Nl，其中Al＝[al(θ1,φ1,γ1,η1),…,al(θk,φk,γk,ηk),…,al(θK,φK,γK,ηK)]为下层子阵的阵列导向矢量，al(θk,φk,γk,ηk)＝ckql(θk,φk)，ql(θk,φk)为下层子阵的相位中心的相对于坐标原点的相位差构成的空域导向矢量，Nl为下层子阵接收的四元数噪声矢量，S为入射信号构成的幅度矩阵；上层子阵的四元数数据矩阵Xu＝AuS+Nu，其中，Au＝[au(θ1,φ1,γ1,η1),…,au(θk,φk,γk,ηk),…,au(θK,φK,γK,ηK)]为上层子阵的阵列导向矢量，au(θk,φk,γk,ηk)＝ckqu(θk,φk)，qu(θk,φk)为上层子阵的相位中心的相对于坐标原点的相位差构成的空域导向矢量，Nu为上层子阵接收的四元数噪声矢量，Au＝AlΦ，Φ为空间相位差因子矩阵；全阵列接收四元数数据为：其中，为全阵列导向矢量矩阵，为全阵列接收的四元数形式的噪声矩阵；步骤三：计算全阵列接收四元数数据Z的自相关矩阵Rz，并对自相关矩阵进行四元数特征分解，得到阵列导向矢量的估计值；Rz=1MZZH=ARsAH+σ2I;]]>其中，(·)H表示转置复共轭操作，为入射信号的自相关函数，I为单位矩阵，σ2为阵元的电偶极子、磁偶极子接收到的高斯白噪声的方差；对Rz进行四元数矩阵特征分解得到信号子空间Es，根据子空间原理，存在K×K的非奇异矩阵T满足：Es＝AT，取Es的前2N行组成矩阵Eu，取后2N行组成矩阵El，由信号子空间的定义，Au、Al与Eu、El之间满足El＝AlT,Eu＝AuT＝AlΦT，通过矩阵运算得到表示El的伪逆矩阵；对矩阵进行四元数特征分解，K个大特征值构成相位差因子矩阵估计值其对应的特征矢量构成非奇异矩阵估计值从而得到上层子阵的阵列导向矢量估计值和下层子阵的阵列导向矢量估计值及全阵列导向矢量估计值步骤四：由上层子阵的阵列导向矢量估计值得到上层子阵相邻阵元间的相位差构成的空域导向矢量估计值根据阵元位置矩阵与相位矩阵的关系，利用最小二乘方法计算信号到达角的估计值；q^(θk,φk)=A^u(2:N,k)./A^u(1:N-1,k);]]>其中，表示的第k列的第2到第N个元素，表示的第k列的第1到第N‑1个元素，./表示对应的元素相除；对取相位得到相位矩阵其中，为x轴方向的方向余弦，为y轴方向的方向余弦，W为阵元的位置矩阵；根据相位矩阵Ω和位置矩阵W的关系得到方向余弦估计值：α^kβ^k=[W]#Ω,]]>式中的[W]#是W的伪逆矩阵；根据方向余弦估计值得到第k个信号的俯仰角估计值和方位角估计值θ^k=arcsin(α^k2+β^k2)φ^k=arctan(α^kβ^k),β^k≥0φ^k=π+arctan(α^kβ^k),β^k<0;]]>步骤五、由上层子阵的阵列导向矢量估计值的实部和三个虚部重构电偶极子子阵的阵列导向矢量估计值和磁偶极子子阵的阵列导向矢量估计值根据子阵的阵列导向矢量间的旋转不变关系得到辅助极化角估计值和极化相位差估计值上层子阵的阵列导向矢量估计值是的实部，是的三个虚部，是四元数的3个虚数单位，根据步骤二中四元数矩阵的构成，重构电偶极子子阵的阵列导向矢量估计值和磁偶极子子阵的阵列导向矢量估计值根据及将步骤四得到信号俯仰角的估计值代入Φz，得到和其中，是阵元中心处z轴方向的电偶极子子阵的阵列导向矢量估计值，是阵元中心处z轴方向的磁偶极子子阵的阵列导向矢量估计值；根据阵元中心的电偶极子子阵和磁偶极子子阵导向矢量间的关系其中为两子阵间的旋转不变矩阵，得到信号极化参数的估计值：γ^k=tan-1(|Ψ^(k,k)|),]]>η^k=arg(-Ψ^(k,k)),]]>其中表示旋转不变矩阵第k行第k列的元素，arg(·)表示取相位；前述步骤中的θk为第k个入射信号的俯仰角，φk为第k个入射信号的方位角，γk为第k个入射信号的辅助极化角，ηk为第k个入射信号的极化相位差，k＝1,…K，n＝1,…,2N。