1.一种均匀圆阵超分辨率空间谱估计方法，其特征在于：包括以下步骤：

S1：由均匀圆阵的阵元位置坐标构造分别位于x轴和y轴的两个虚拟均匀线阵作为参考阵列；

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

S3：将方位角全角域均匀划分为L个区间，即设定方位角搜索步长为τ＝2π/L；

S4：确定搜索方位角并将该角度指定为相位补偿角，计算出基于该相位补偿角的两个相位补偿矩阵，用两个相位补偿矩阵分别对步骤S2中估计得到的信号子空间进行相位补偿，得到两个经过相位补偿后的信号子空间，接着对经过相位补偿后的两个信号子空间施行经典ESPRIT算法，得到两个子空间拟合矩阵，然后对两个子空间拟合矩阵进行特征值分解，获得相应的两个特征向量矩阵和对应的特征值；

S5：用其中一个特征向量矩阵的共轭转置矩阵与另一个特征向量矩阵相乘得到一个乘积矩阵，确定该乘积矩阵模值最大的数据元素，并通过该数据确定对应的两个拟合矩阵的特征值；由这两个特征值联立解出一个角度估值，再联合最大模数据元素，两个特征值计算出该相位补偿角上的扩展ESPRIT空间谱值；

S6：变更搜索方位角，重复步骤S4-S5，获得方位角全角域内的空间谱，确定该空间谱最大的N个谱峰位置，对应的即为N个信号的DoA估计。

2.根据权利要求1所述的均匀圆阵超分辨率空间谱估计方法，其特征在于：步骤S1包括以下步骤：S11：设定实际阵列阵元数为M，所构造的两个参考阵列分别由放置在x轴的虚拟均匀线阵A和y轴上的虚拟均匀线阵B组成，每个虚拟均匀线阵的阵元数都和实际阵列阵元数相同；

x轴上的虚拟均匀线阵的两个端阵元的x坐标分别等于实际阵元中x坐标的最大值和最小值，其余虚拟阵元在x轴上均匀分布；y轴上的虚拟均匀线阵的两个端阵元的y坐标分别等于实际阵元中y坐标的最大值和最小值，其余虚拟阵元在y轴上均匀分布；

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

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

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

3.根据权利要求1所述的均匀圆阵超分辨率空间谱估计方法，其特征在于：所述步骤S3包括：设定在方位角[0 2π)区间的搜索点数目为L，则搜索步长为τ＝2π/L，因此所有的方位角搜索角度为0，2π/L，4π/L，...，2π(L-1)/L。

4.根据权利要求1所述的均匀圆阵超分辨率空间谱估计方法，其特征在于：在步骤S4中，所述的确定搜索方位角并将该角度指定为相位补偿角，计算出基于该相位补偿角的两个相位补偿矩阵，用两个相位补偿矩阵分别对步骤S2中估计得到的信号子空间进行相位补偿，得到两个经过相位补偿后的信号子空间，包括以下步骤：S41：确定第l(l＝1,2,..,L)个方位角搜索角度，并把它设为相位补偿角

再根据实际阵列与x轴上虚拟参考均匀线阵的阵元对应关系确定一个M×M

选择矩阵JA，该选择矩阵每一行每一列都仅有一个数据元素为1，其余元素为0，元素“1”的行号表示以该行号编号的虚拟参考阵列阵元，其列号表示对应的实际阵元位置编号，即JA＝[e1,eM,e2,eM-1,...]T其中，ej表示第j个元素为1的M×1维标准单位列向量，然后用选择矩阵与实际阵列位置矩阵相乘，再减去x轴上虚拟均匀线阵位置矩阵，得到实际阵列与x轴参考均匀线阵的位置差矩阵ΔA：ΔA＝[ΔAx,ΔAy]：＝JA[x,y]-[xA,yA]

其中ΔAx为实际阵列与位于x轴的参考虚拟均匀线阵x坐标的位置差向量，ΔAy为实际阵列与位于x轴的参考虚拟均匀线阵y坐标的位置差向量，再根据该位置差矩阵和相位补偿角 定义一个对角矩阵作为x轴参考阵的相位补偿矩阵：其中，λ为信号工作波长；

S42：根据实际阵列与y轴上参考虚拟线阵的阵元对应关系确定另一个M×M的选择矩阵JB，该选择矩阵每一行每一列也仅有一个数据元素为1，其余元素为0，元素“1”的行号表示以该行号编号的虚拟参考阵列阵元，其列号表示对应的实际阵元位置编号，即JB＝[eP,eP-1,eP+1,eP-2,eP+2...]T其中P表示实际阵列中最大y坐标阵元的编号，然后用选择矩阵与实际阵列位置矩阵相乘，再减去y轴上虚拟均匀线阵位置矩阵，得到实际阵列与y轴参考阵列的位置差矩阵ΔB：ΔB＝[ΔBx,ΔBy]：＝JB[x,y]-[xB,yB]

其中，ΔBx为实际阵列与位于y轴的参考虚拟均匀线阵x坐标位置差向量，ΔBy为实际阵列与位于y轴的参考虚拟均匀线阵y坐标位置差向量，再根据该位置差矩阵和相位补偿角定义一个对角矩阵作为y轴参考阵的相位补偿矩阵：S43：将对应x轴参考阵的相位补偿矩阵ΓA先乘上对应x轴上参考均匀线阵的选择矩阵JA，再与S2中的信号子空间相乘完成基于x轴参考阵的相位补偿操作；将对应y轴参考阵的相位补偿矩阵ΓB先乘上对应y轴上参考线阵的选择矩阵JB，再与步骤S2中的信号子空间相乘完成基于y轴参考阵的相位补偿操作。

5.根据权利要求1所述的均匀圆阵超分辨率空间谱估计方法，其特征在于：在步骤S4中，所述对经过相位补偿后的两个信号子空间施行经典ESPRIT算法，得到两个子空间拟合矩阵，然后对两个子空间拟合矩阵进行特征值分解，获得相应的两个特征向量矩阵和对应的特征值，包括以下步骤：S44：选择S43中已经完成的对应x轴上参考阵列相位补偿矩阵相位更新后的信号子空间的前M-1个行向量作为x轴参考阵列子阵列1的信号子空间，选择该相位更新信号子空间的后M-1个行向量作为x轴参考阵列子阵列2的信号子空间，然后利用最小二乘方法计算这两个子阵列信号空间的拟合矩阵，该拟合矩阵即为基于x轴参考阵的拟合矩阵；

S45：选择S43中已经完成的对应y轴上参考阵列相位补偿矩阵相位更新后的信号子空间的前M-1个行向量作为y轴参考阵列子阵列1的信号子空间，选择该相位更新信号子空间的后M-1个行向量作为y轴参考阵列子阵列2的信号子空间，然后利用最小二乘方法计算这两个子阵列信号空间的拟合矩阵，该拟合矩阵即为基于y轴参考阵的拟合矩阵；

S46：对两个子空间拟合矩阵进行特征值分解，每一个拟合矩阵分解得到的所有特征向量组成的矩阵即为特征向量矩阵并且得到这些特征向量对应的特征值。

6.根据权利要求1所述的均匀圆阵超分辨率空间谱估计方法，其特征在于：所述步骤S5包括以下步骤：S51：空间信号数量为N，S4中获得的两个特征向量矩阵均为N×N维矩阵，用一个特征矩阵的共轭转置矩阵与另一个特征向量矩阵相乘，获得的乘积矩阵有N2个元素，确定这N2个数据中模值最大的那个元素Ωuv同时确定对应的两个特征值 和S52：联立特征值 和 解得一个角度估计

其中， 表示对应特征值的相角；

S53：利用最大模值数据Ωuv，特征值 和 以及角度估计 计算出方位角 处的空间谱值：