1.电磁偶极子对阵列的四元数ESPRIT参数估计方法，其特征在于，包括以下步骤：K个互不相关完全极化横电磁波信号同时入射到电磁矢量传感器阵列上，所述阵列阵元是由一个电偶极子和一个磁偶极子组成的电磁偶极子对，步骤一、对电磁矢量传感器阵列的输出信号进行M次采样得到第一组采样数据X，延时ΔT后同步采样M次，得到第二组采样数据Y；

其中，xne(m)表示第n个阵元的电偶极子输出信号的第m次采样数据，xnh(m)表示第n个阵元的磁偶极子输出信号的第m次采样数据，yne(m)表示延时ΔT后第n个阵元的电偶极子输出信号的第m次采样数据，ynh(m)表示延时ΔT后第n个阵元的磁偶极子输出信号的第m次采样数据；

步骤二、将第一组采样数据X和第二组采样数据Y分别按照同阵元的电偶极子和磁偶极子的同次快拍数据叠加构成第一组接收四元数数据矩阵Z1和第二组接收四元数数据矩阵Z2，构造全阵列接收数据矩阵Z；

其中， 表示由第n个阵元的电偶极子和磁偶极子输出信号的第m次采样数据叠加构成的四元数数据， 表示延时ΔT由第n个阵元的电偶极子和磁偶极子输出信号的第m次采样数据叠加构成的四元数数据；

第一组采样数据X构成的四元数数据矩阵Z1＝A1S+N1，其中，A1＝[a1(θ1,φ1,γ1,η1),…,a1(θk,φk,γk,ηk),…,a1(θK,φK,γK,ηK)]是第一组采样数据对应的阵列导向矢量，a1(θk,φk,γk,ηk)＝ckq(θk,φk)， 是第k个入射信号的第一组采样数据中电磁场的四元数数据表示， 和 分别为第k个入射信号在坐标原点处沿z轴方向的电场分量和磁场分量，q(θk,φk)为全阵列相位中心的空域导向矢量，θk是第k个入射信号的俯仰角，φk是第k个入射信号的方位角，γk是第k个入射信号的辅助极化角，ηk第k个入射信号的极化相位差，N1是第一组采样数据的高斯白噪声矢量，S为入射信号构成的幅度矩阵；

第二组采样数据Y构成的四元数数据矩阵Z2＝A2S+N2＝A1ΦS+N2，其中，A2＝A1Φ是第二组采样数据对应的阵列导向矢量，Φ为时延矩阵，N2是第二组采样数据的高斯白噪声矢量；

构造全阵列接收数据矩阵 其中， 是全数据对

应的阵列导向矢量矩阵， 是全数据噪声矩阵；

步骤三、计算全阵列接收数据矩阵的自相关矩阵Rz，对自相关矩阵进行四元数特征分解，得到第一组采样数据对应的阵列导向矢量估计值 第二组采样数据对应的阵列导向矢量估计值 和全数据对应的阵列导向矢量矩阵估计值其中， 为入射信号的自相关函数，σ2为噪声的方差，I为单位矩阵，(·)H表示转置复共轭操作；

对自相关矩阵Rz进行四元数矩阵特征分解得到信号子空间Es，根据子空间原理，存在K×K的非奇异矩阵T，且Es＝AT，分别取Es的前N行及后N行分别组成的矩阵E1和E2，由信号子空间的定义，A1、A2与E1、E2之间满足E1＝A1T，E2＝A2T＝A1ΦT，则有对矩阵 进行四元数特征分解，K个大特征值构成延时矩阵估计值 特征值对应的特征矢量构成非奇异矩阵估计值 根据 得到步骤四、计算信号到达角的估计值；

根据 计算相邻两阵元间的相位差矢量 其中，A1

(2:N,k)表示A1的第k列的第2到第N个元素，A1(1:N-1,k)表示A1的第k列的第1到第N-1个元素，./表示对应的元素相除；

计算相位矩阵 arg(·)表示取相位；

根据 计算第k个入射信号的x轴方向的方向余弦估计值 和y轴方向的方向余弦估计值 [W]#是位置矩阵W的伪逆矩阵；

根据方向余弦的估计值 得到信号到达角的估计值：

步骤五、由 的实部和三个虚部重构Z轴方向的电偶极子子阵的阵列导向矢量估计值和磁偶极子子阵的阵列导向矢量估计值 根据子阵导向矢量间的旋转不变关系得到信号极化参数的估计值；

是 的实部， 是 的三个虚部，根据步骤二中

四元数矩阵的构成， 重构Z轴方向的电偶极子子阵的导向矢量估计值和磁偶极子子阵导向矢量估计值 两子阵导向矢量间的旋转

不变关系为 为两子阵间的旋转不变矩阵，根据 计算信号极化参数的估计值：其中， 表示旋转不变矩阵 的第k行第k列元素；

前述步骤中的n＝1,…,N，m＝1,…,M，N为阵列的阵元数，M为采样次数， 为四元数的3个虚数单位。

2.根据权利要求1所述的电磁偶极子对阵列的四元数ESPRIT参数估计方法，其特征在于：所述阵列为圆环形阵列，电偶极子的轴线和磁偶极子的轴线平行于z轴，N个阵元均匀分布在圆环上，坐标原点位于圆环的圆心。