[发明专利]一种均匀圆阵二维波达方向快速估计方法

权利要求书

1.一种均匀圆阵二维波达方向快速估计方法，其特征在于：该方法包括以下步骤：

S1：由均匀圆阵的阵元位置坐标构造一个十字正交虚拟参考阵列；

S2：将二维角度参数空间均匀划分成多个网格，确定每个网格的中心方向；

S3：由接收到的阵列快拍数据获得信号子空间的估计；

S4：将步骤S2中划分得到的网格进行顺序编号，选定最小编号网格；

S5：指定选定网格的中心方向为相位补偿角，并记为对应该网格的初始相位补偿角；

S6：计算出基于相位补偿角的两个相位补偿矩阵，并用两个相位补偿矩阵分别对步骤S3中估计得到的信号子空间进行相位补偿，得到两个经过相位补偿后的信号子空间，接着对经过相位补偿后的两个信号子空间施行经典2D ESPRIT算法，得到两个子空间拟合矩阵，然后对两个子空间拟合矩阵进行特征值分解，获得相应的两个特征向量矩阵和对应的特征值；

S7：用其中一个特征向量矩阵的共轭转置矩阵与另一个特征向量矩阵相乘得到一个乘积矩阵，判断该乘积矩阵的元素的实部最大值是否大于设定的阈值a，若大于阈值a则结合步骤S6中得到的两个对应特征值联立解出一个二维方向角对并转向S8，否则选定下一个网格并转向步骤S5；

S8：比较步骤S7中估得的二维方向角与步骤S6中计算相位补偿矩阵的二维补偿角绝对差，与设定的另一个阈值b的大小，若该绝对差大于阈值b则将估得的该二维方向角设定为新的相位补偿角并转向步骤S6，否则转向步骤S9；

S9：判断步骤S7中估得的二维方向角与步骤S5中设置的初始相位补偿角的方位角绝对差是否大于方位角方向网格间隔的一半，并且判断俯仰角绝对差是否大于俯仰角方向上网格间隔的一半，若两者之一不成立则选定下一个网格并转向步骤S5，若两者同时成立则输出该二维方向角为信号的一个二维波达方向估计，同时选定下一个网格并转向步骤S5。

2.根据权利要求1所述的一种均匀圆阵二维波达方向快速估计方法，其特征在于：所述步骤S1具体为：

S101：设定实际阵列阵元数为M，十字正交虚拟参考阵列由放置在x轴的虚拟均匀线阵A和y轴上的虚拟均匀线阵B组成，每个虚拟均匀线阵的阵元数都和实际阵列阵元数相同；x轴上的虚拟均匀线阵的两个端阵元的x坐标分别等于实际阵元中x坐标的最大值和最小值，其余虚拟阵元在x轴上均匀分布；y轴上的虚拟均匀线阵的两个端阵元的y坐标分别等于实际阵元中y坐标的最大值和最小值，其余虚拟阵元在y轴上均匀分布；

S102：x轴上的虚拟均匀线阵的阵元与实际阵元的对应关系按x坐标的大小顺序关系一一对应，若实际阵元中有多个阵元x坐标相同，则先对应这些x坐标相同阵元中最小y坐标的阵元，再对应这些x坐标相同阵元中最大y坐标的阵元，然后再对应这些x坐标相同阵元中倒数第二小y坐标的阵元，即相同x坐标则按y坐标大小交叉安排对应；

S103：y轴上的虚拟均匀线阵的阵元与实际阵元的对应关系按y坐标的大小关系一一对应，若实际阵元中有多个阵元y坐标相同，则先对应这些y坐标相同阵元中最小x坐标的阵元，再对应这些y坐标相同阵元中最大x坐标的阵元，然后再对应这些y坐标相同阵元中倒数第二小x坐标的阵元，即相同y坐标则按x坐标大小交叉安排对应；

S104：计算得出x轴上虚拟参考线阵相邻阵元之间的距离d_x，即最大x坐标减去最小x坐标，差值再除以间隔个数，即M-1，再计算得出y轴上虚拟参考线阵相邻阵元之间的距离d_y，即最大y坐标减去最小y坐标，差值再除以间隔个数，即M-1。

3.根据权利要求1所述的一种均匀圆阵二维波达方向快速估计方法，其特征在于：所述步骤S2具体为：根据预先设定的方位角方向间隔参数和俯仰角方向间隔参数，将方位角区间[0,360°)等间距划分，同时将俯仰角区间[0,90°)等间距划分，即空间角被均匀划分为多个网格，每个网格中心点对应的方位角和俯仰角组成(方位角，俯仰角)对，即为该网格中心方向。