1.一种幂迭代辅助Root‑MUSIC的快速MIMO测向方法，其特征在于，具体包括以下步骤：
步骤1、初始化接收机结构的基本参数，建立大规模MIMO接收机的系统模型；
步骤2、选取全阵列的一部分天线对辐射源方向进行初步估计，通过ESPRIT算法估计入射波信号的初始角度，基于估计值拓展构建全阵列初始阵列流形；
步骤3、在全阵列初始阵列流形的基础上，对全部阵元采用幂迭代方法求解准确的阵列流形；
步骤4、根据准确的阵列流形，采用Root‑MUSIC算法估计最终的到达方向DOA；
步骤1中初始化接收机结构的基本参数，包括：
初始化接收远场的窄带入射波信号为 其中s(t)为基带信号，wc＝2πfc，fc为信号的载波频率，初始输入的到达角为θ；
步骤1中建立大规模MIMO接收机的系统模型，具体如下：
接收机采用全数字天线阵列结构，接收机阵列含M根天线，经过射频链路和模数转换，单个发射源情形下，接收信号形式表达为：
T
z(t)＝[z1(t),z2(t),…,zM(t)]＝a(θ)s(t)+n(t)
式中，z1(t),z2(t),…,zM(t)依次为第1,2,…,M根天线的接收信号；n(t)是均值为0、方
2
差为σ的高斯白噪声；a(θ)为阵列流形，表示为：
式中，λ是载波频率的波长，d等于载波频率的半波长即
步骤2所述选取全阵列的一部分天线对辐射源方向进行初步估计，通过ESPRIT算法估计入射波信号的初始角度，基于估计值拓展构建全阵列初始阵列流形，具体如下：
取全接收阵列的一部分 根天线来进行初步估计DOA，M0根天线对应的接收信号向量 表示为：
2
式中， 为M0根天线对应的阵列流形，n(t)为M0根天线对应的均值为0、方差为σ 的高斯白噪声；
的样本协方差矩阵和相应的特征分解由下式给出：
式中，i表示第i个样本，L为样本数量，U能够分解成信号子空间US和噪声子空间UN，Λ为特征值构成的对角阵；
通过ESPRIT算法估计入射波信号的初始角度θ0：矩阵U和RM0的主特征值对应的部分组成US，抽取US的前M0‑1行组成矩阵U1，后M0‑1行组成矩阵U2，构造一个矩阵Ψ并进行特征值分解：
通过矩阵Ψ的主特征值 得到入射波信号的初始角度θ0：
式中，为特征值的相位；
基于估计的入射波信号的初始角度，将阵列流形扩展到M维，构造全阵列初始阵列流形a(θ0)如下
2.根据权利要求1所述的幂迭代辅助Root‑MUSIC的快速MIMO测向方法，其特征在于，步骤3所述在全阵列初始阵列流形的基础上，对全部阵元采用幂迭代方法求解准确的阵列流形，具体如下：
将a(θ0)作为幂迭代方法的初始迭代向量x0，得：
n n
xn＝Rxn‑1＝Rx0＝Ra(θ0)
式中，xn为第n次迭代结果，R为协方差矩阵；
假设上述幂迭代在迭代n次时收敛，则V＝xn是准确的阵列流形。
3.根据权利要求2所述的幂迭代辅助Root‑MUSIC的快速MIMO测向方法，其特征在于，步骤4所述根据准确的阵列流形，采用Root‑MUSIC算法估计最终的到达方向DOA，具体如下：
采用Root‑MUSIC算法估计最终的到达方向DOA，求根多项式表示为：
M‑1 T
式中p(z)＝[1,z,…,z ]，矩阵UN为噪声子空间；
通过求取p(z)的全部根值，在单位圆内选取距离最近的根 得到最终的到达方向DOA为
4.一种幂迭代辅助Root‑MUSIC的快速MIMO测向装置，其特征在于，包括：
模型构建模块，用于初始化接收机结构的基本参数，建立大规模MIMO接收机的系统模型；
ESPRIT算法模块，用于选取全阵列的一部分天线对辐射源方向进行初步估计，通过ESPRIT算法估计入射波信号的初始角度，基于估计值拓展构建全阵列初始阵列流形；
幂迭代求解模块，用于在全阵列初始阵列流形的基础上，对全部阵元采用幂迭代方法求解准确的阵列流形；
Root‑MUSIC算法模块，用于根据准确的阵列流形，采用Root‑MUSIC算法估计最终的到达方向DOA；
初始化接收机结构的基本参数，包括：
初始化接收远场的窄带入射波信号为 其中s(t)为基带信号，wc＝2πfc，fc为信号的载波频率，初始输入的到达角为θ；
建立大规模MIMO接收机的系统模型，具体如下：
接收机采用全数字天线阵列结构，接收机阵列含M根天线，经过射频链路和模数转换，单个发射源情形下，接收信号形式表达为：
T
z(t)＝[z1(t),z2(t),…,zM(t)]＝a(θ)s(t)+n(t)
式中，z1(t),z2(t),…,zM(t)依次为第1,2,…,M根天线的接收信号；n(t)是均值为0、方
2
差为σ的高斯白噪声；a(θ)为阵列流形，表示为：
式中，λ是载波频率的波长，d等于载波频率的半波长即
选取全阵列的一部分天线对辐射源方向进行初步估计，通过ESPRIT算法估计入射波信号的初始角度，基于估计值拓展构建全阵列初始阵列流形，具体如下：
取全接收阵列的一部分 根天线来进行初步估计DOA，M0根天线对应的接收信号向量 表示为：
2
式中， 为M0根天线对应的阵列流形，n(t)为M0根天线对应的均值为0、方差为σ 的高斯白噪声；
的样本协方差矩阵和相应的特征分解由下式给出：
式中，i表示第i个样本，L为样本数量，U能够分解成信号子空间US和噪声子空间UN，Λ为特征值构成的对角阵；
通过ESPRIT算法估计入射波信号的初始角度θ0：矩阵U和 的主特征值对应的部分组成US，抽取US的前M0‑1行组成矩阵U1，后M0‑1行组成矩阵U2，构造一个矩阵Ψ并进行特征值分解：
通过矩阵Ψ的主特征值 得到入射波信号的初始角度θ0：
式中，为特征值的相位；
基于估计的入射波信号的初始角度，将阵列流形扩展到M维，构造全阵列初始阵列流形a(θ0)如下
5.一种移动终端，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至3中任一项的幂迭代辅助Root‑MUSIC的快速MIMO测向方法。
6.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现如权利要求1至3中任一项所述幂迭代辅助Root‑MUSIC的快速MIMO测向方法中的步骤。