CN111707986A - 一种基于稀疏面阵的三维参数估计方法 - Google Patents

Abstract

本发明提出了一种基于稀疏面阵的三维参数估计方法。首先，将面阵所在平面按照信号入射角的范围分为8个区域。在每一个区域内，通过分析稀疏面阵，将稀疏面阵按照信号的入射方向映射为非均匀虚拟线阵，构造稀疏面阵的方向矢量并计算稀疏面阵和非均匀虚拟线阵之间的相位差。其次，将面阵方向矢量分别乘上面阵与虚拟线阵之间的相位差，得到虚拟线阵的方向矢量。在此基础上，构造入射信号并利用三维参数估计算法进行AoA(Arrival of Angle)、ToF(Time of Flight)和DFS(Doppler Frequency Shift)的联合参数估计，并利用谱函数搜索得到一系列的峰值。最后，通过分析稀疏面阵映射为非均匀虚拟线阵的几何关系，利用角度搜索得到正确的入射角对应的峰值。本发明克服了业务天线稀疏面阵排列规则不满足空间采样定理而导致的无法进行参数估计的问题，为实际应用中基于业务天线的室内跟踪定位等应用奠定了理论基础。

Description

一种基于稀疏面阵的三维参数估计方法

技术领域

本发明属于参数估计方法，具体涉及在Wi-Fi系统下，一种利用稀疏面阵对室内目标和收发机组成的系统的参数估计方法。

背景技术

近年来，目标不携带任何设备的无源被动跟踪定位技术备受关注，旨在对室内的人进行定位跟踪，具体的应用包括老年人、病人安全监测、智能家居和许多其他基于物联网(IoT)的应用。现有的各种室内目标跟踪方法大多要求目标携带专用设备或可穿戴设备，但是在某些情况下携带设备是很不方便的。

目前，基于Wi-Fi的跟踪定位系统研究一直是人们关注的焦点，这种系统不需要额外的基础设施，只需要Wi-Fi接入点(AP)和一个或几个支持Wi-Fi协议(如802.11n/ac)的接收设备并分别布置于不同的环境。检测环境中存在的人体会对Wi-Fi信号的传输环境造成一定程度的影响，而且CSI(Channel State Information)可以细粒化的记录Wi-Fi信号的变化情况，通过目标上的信号反射来提取基本运动和位置信息。但是无源跟踪比有源无线发射器的定位更具挑战性，因为由人体反射的反射信号通常比直接路径信号弱几个数量级，并且通常与强直接路径信号以及从墙壁、家具和其他附近的杂物反射的信号叠加。因此很难从反射信号中提取有用的、准确的定位信息。而且现有的利用WiFi进行参数估计的方法都是基于满足空间采样定理的均匀阵列。但是在实际的应用场景中，大部分的天线阵列都是天线分布不规则不满足空间采样定理的业务天线阵列，当天线阵列不满足空间采样定理时，其估计出的参数会存在较大的伪峰，造成参数估计正确的峰值的误判，从而影响参数估计的结果。针对上述问题，本发明在WiFi系统下，设计了一种将稀疏面阵映射为非均匀线阵来进行参数估计的方法。在现有的均匀阵列的参数估计的基础上用稀疏面阵进行参数估计，更具有实用性和普适性。同时利用面阵映射为线阵的角度去除稀疏阵列导致的伪峰，从而提高了参数估计的精度。

发明内容

本发明的目的是在Wi-Fi系统下，提供一种利用稀疏面阵对室内目标和收发机组成的系统的参数估计方法，它能在不满足空间采样定理的业务天线阵列上进行精确的参数估计，对室内目标进行跟踪定位。

本发明所述的一种基于稀疏面阵的三维联合参数估计方法，具体包括以下步骤：

步骤一：假设在Wi-Fi系统中有D个信号源，N个子载波，由于室内Wi-Fi单个AP业务天线数通常为4根，因此，接收端为由4根天线组成2×2的方形面阵，接收来自信号源的P个数据包中的CSI信息，天线间距为λ，不满足空间采样定理。将4根天线从上至下从左到右依次编号为天线1、2、3、4。θ_i表示信号源i∈{1,2,…,D}的入射角，不同信号源入射角度向量可以表示为Θ＝{θ₁,θ₂,…,θ_D}。

步骤二：将面阵所在平面按照信号入射方向进行分块，将信号入射角θ的大小分为0°～45°、45°～90°、90°～135°、135°～180°、180°～225°、225°～270°、270°～315°、315°～360°8个区域，分别对每个区域进行映射。

步骤三：沿着信号入射方向将稀疏面阵映射为非均匀得线阵。

首先当信号从区域一入射到阵列时，将天线2和天线3沿着信号入射方向投影到天线1和天线4的连线上，投影之后的虚拟天线与天线1和天线4构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_subcarrier(t)＝[1 e^-j2πft … e^-j2π(N-1)ft]^T

Φ_array1(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差，Φ_subcarrier(t)表示信号由不同的子载波传输导致的相位差，Φ_packet(v)表示不同的数据包之间的相位差。

以天线1作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

由于入射角未知θ_i，而面阵和线阵之间的相位差ΔΦ₁(θ_i)与入射角有关，因此，我们让入射角θ_i从0°～45°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₁(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

当信号从区域二入射到阵列时，将天线2和天线3沿着信号入射方向投影到天线1和天线4的连线上，投影之后的虚拟天线与天线1和天线4构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_array2(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差。

以天线1作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

让入射角θ_i从45°～90°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₂(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

当信号从区域三入射到阵列时，将天线1和天线4沿着信号入射方向投影到天线2和天线3的连线上，投影之后的虚拟天线与天线2和天线3构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_array3(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差。

以天线3作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

让入射角θ_i从90°～135°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₃(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

当信号从区域四入射到阵列时，将天线1和天线4沿着信号入射方向投影到天线2和天线3的连线上，投影之后的虚拟天线与天线2和天线3构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_array4(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差。

以天线3作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

让入射角θ_i从135°～180°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₄(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

信号从区域5、6、7、8入射到整列上和信号从区域1、2、3、4入射到阵列上是相反的方向，映射方式完全相同。

步骤四：当信号从区域i,(1≤i≤8)入射到阵列时，得到虚拟线阵的方向矢量

将得到的虚拟线阵的方向矢量构造入射信号

计算入射信号的协方差矩阵：

其中H表示共轭转置，S为入射信号，

为M×D维的阵列流型矩阵，M和D分别为天线和入射信号的个数。W为高斯白噪声，其概率密度服从正态分布，

均值μ＝0，方差为σ²。

对矩阵S进行特征分解，得到NMP个特征值，其中有NMP-D个较小的特征值等于噪声的方差σ²，NMP-D个较小的特征值只与噪声有关，假设λ_min为最小特征值，S₀是矩阵S的特征向量，则有λ_minS₀＝σ²I。NMP-D个较小的特征值对应的特征向量就构成了噪声子空间。

当入射信号的个数D少于子载波数量和天线阵列数和数据包数的乘积NMP时，

的秩小于NMP，因此：

在S＝λS₀两边同时乘e_i，e_i为最小特征值对应的特征向量，得到：

Se_i＝λ_iS₀e_i,i＝D+1,...,NMP

因为：

则：

即：

由于

S_F≠0，所以：

由上式可以得到矩阵S的所有最小特征向量e_i与矩阵

的各列向量正交。最后得到相互正交的噪声子空间和信号子空间。将噪声子空间中的特征向量构造一个NMP×(NMP-D)维的噪声特征向量矩阵E_N，a(θ,τ,v)为与映射得到的线阵的阵列流型相同的方向矢量，θ、τ、v为待估计的参数，分别为入射信号的到达角、飞行时间和速度。计算E_N与矢量a(θ,τ,v)构成的谱函数：

其中，H为共轭转置符号。搜索谱函数f(θ,τ,v)得到一系列的峰值。

步骤五：因为稀疏面阵的天线空间排列不满足空间采样定理，因此通过步骤四搜索得到的峰值中存在很多的伪峰值。只有当搜索的角度θ_i和面阵映射为线阵的角度相等时，得到的相位差ΔΦ(θ_i)才能正确地表达稀疏面阵和虚拟非均匀线阵之间的关系。此时，通过参数求解得到的结果才是正确的且等于θ_i。因此，在得到的峰值中，将每一个峰值通过映射的几何关系进行反解，我们就可以得出入射角θ_i。首先通过峰值搜索得到了一系列的峰值SP(i,j,k)，其中i，j，k分别代表峰值中多普勒速度、到达角和飞行时间对应的值。根据稀疏面阵和映射得到的虚拟非均匀线阵之间的几何关系，可以计算出信号入射到稀疏面阵上的角度和信号入射到虚拟非均匀线阵上的角度之间的关系。当信号从区域一、二入射时，信号入射到稀疏面阵上的角度和信号入射到虚拟非均匀线阵上的角度之间的关系为：

当信号从区域三、四入射时，信号入射到稀疏面阵上的角度和信号入射到虚拟非均匀线阵上的角度之间的关系为：

其中

为信号入射到映射后的非均匀虚拟线阵上的角度，θ表示信号入射到稀疏面阵上的角度。信号入射到非均匀虚拟线阵上的角度

包含在峰值搜索得到的一系列峰值SP(i,j,k)中，将每一个峰值对应的角度与角度

进行对比就能得到正确估计的AoA(Arrival of Angle)、ToF(Time of Flight)和DFS(Doppler Frequency Shift)三个参数。

有益效果

首先，为了克服稀疏面阵的天线排列规则不满足空间采样定理导致的参数估计的准确率不高的问题，分析了稀疏面阵的空间排列规则以及将面阵映射为线阵的几何关系以及现有的参数估计方法，证明了该方法对提高用稀疏面阵进行参数估计的准确性有效。本方法首先通过分析稀疏面阵，将稀疏面阵按照信号的入射方向映射为非均匀虚拟线阵，构造稀疏面阵的方向矢量并计算稀疏面阵和非均匀虚拟线阵之间的相位差。其次，将面阵方向矢量分别乘上面阵与虚拟线阵之间的相位差，得到虚拟线阵的方向矢量。在此基础上，构造入射信号并利用三维参数估计算法进行AoA、ToF和DFS的联合参数估计，并利用谱函数搜索得到一系列的峰值。最后，通过分析稀疏面阵映射为非均匀虚拟线阵的几何关系，利用角度搜索得到正确的入射角对应的峰值。本发明克服了业务天线稀疏面阵排列规则不满足空间采样定理而导致的无法进行参数估计的问题，为实际应用中基于业务天线的室内跟踪定位等应用奠定了理论基础。

附图说明

图1为基于稀疏面阵的三维联合估计的具体实施流程。

图2为信号入射到稀疏面阵示意图。

图3为参数估计结果图。

具体实施方案

步骤一：假设在Wi-Fi系统中有D个信号源，N个子载波，由于室内Wi-Fi单个AP业务天线数通常为4根，因此，接收端为由4根天线组成2×2的方形面阵，接收来自信号源的P个数据包中的CSI信息，天线间距为λ，不满足空间采样定理。将4根天线从上至下从左到右依次编号为天线1、2、3、4。θ_i表示信号源i∈{1,2,…,D}的入射角，不同信号源入射角度向量可以表示为Θ＝{θ₁,θ₂,…,θ_D}。

步骤二：将面阵所在平面按照信号入射方向进行分块，将信号入射角θ的大小分为0°～45°、45°～90°、90°～135°、135°～180°、180°～225°、225°～270°、270°～315°、315°～360°8个区域，分别对每个区域进行映射。

步骤三：沿着信号入射方向将稀疏面阵映射为非均匀得线阵。

首先当信号从区域一入射到阵列时，将天线2和天线3沿着信号入射方向投影到天线1和天线4的连线上，投影之后的虚拟天线与天线1和天线4构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_subcarrier(t)＝[1 e^-j2πft … e^-j2π(N-1)ft]^T

Φ_array1(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差，Φ_subcarrier(t)表示信号由不同的子载波传输导致的相位差，Φ_packet(v)表示不同的数据包之间的相位差。

以天线1作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

由于入射角未知θ_i，而面阵和线阵之间的相位差ΔΦ₁(θ_i)与入射角有关，因此，我们让入射角θ_i从0°～45°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₁(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

当信号从区域二入射到阵列时，将天线2和天线3沿着信号入射方向投影到天线1和天线4的连线上，投影之后的虚拟天线与天线1和天线4构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_array2(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差。

以天线1作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

让入射角θ_i从45°～90°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₂(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

当信号从区域三入射到阵列时，将天线1和天线4沿着信号入射方向投影到天线2和天线3的连线上，投影之后的虚拟天线与天线2和天线3构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_array3(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差。

以天线3作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

让入射角θ_i从90°～135°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₃(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

当信号从区域四入射到阵列时，将天线1和天线4沿着信号入射方向投影到天线2和天线3的连线上，投影之后的虚拟天线与天线2和天线3构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_array4(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差。

以天线3作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

让入射角θ_i从135°～180°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₄(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

信号从区域5、6、7、8入射到整列上和信号从区域1、2、3、4入射到阵列上是相反的方向，映射方式完全相同。

步骤四：当信号从区域i,(1≤i≤8)入射到阵列时，得到虚拟线阵的方向矢量

将得到的虚拟线阵的方向矢量构造入射信号

计算入射信号的协方差矩阵：

其中H表示共轭转置，S为入射信号，

为M×D维的阵列流型矩阵，M和D分别为天线和入射信号的个数。W为高斯白噪声，其概率密度服从正态分布，

均值μ＝0，方差为σ²。

对矩阵S进行特征分解，得到NMP个特征值，其中有NMP-D个较小的特征值等于噪声的方差σ²，NMP-D个较小的特征值只与噪声有关，假设λ_min为最小特征值，S₀是矩阵S的特征向量，则有λ_minS₀＝σ²I。NMP-D个较小的特征值对应的特征向量就构成了噪声子空间。

当入射信号的个数D少于子载波数量和天线阵列数和数据包数的乘积NMP时，

的秩小于NMP，因此：

在S＝λS₀两边同时乘e_i，e_i为最小特征值对应的特征向量，得到：

Se_i＝λ_iS₀e_i,i＝D+1,...,NMP

因为：

则：

即：

由于

S_F≠0，所以：

由上式可以得到矩阵S的所有最小特征向量e_i与矩阵

的各列向量正交。最后得到相互正交的噪声子空间和信号子空间。将噪声子空间中的特征向量构造一个NMP×(NMP-D)维的噪声特征向量矩阵E_N，a(θ,τ,v)为与映射得到的线阵的阵列流型相同的方向矢量，θ、τ、v为待估计的参数，分别为入射信号的到达角、飞行时间和速度。计算E_N与矢量a(θ,τ,v)构成的谱函数：

其中，H为共轭转置符号。搜索谱函数f(θ,τ,v)得到一系列的峰值。

步骤五：因为稀疏面阵的天线空间排列不满足空间采样定理，因此通过步骤四搜索得到的峰值中存在很多的伪峰值。只有当搜索的角度θ_i和面阵映射为线阵的角度相等时，得到的相位差ΔΦ(θ_i)才能正确地表达稀疏面阵和虚拟非均匀线阵之间的关系。此时，通过参数求解得到的结果才是正确的且等于θ_i。因此，在得到的峰值中，将每一个峰值通过映射的几何关系进行反解，我们就可以得出入射角θ_i。首先通过峰值搜索得到了一系列的峰值SP(i,j,k)，其中i，j，k分别代表峰值中多普勒速度、到达角和飞行时间对应的值。根据稀疏面阵和映射得到的虚拟非均匀线阵之间的几何关系，可以计算出信号入射到稀疏面阵上的角度和信号入射到虚拟非均匀线阵上的角度之间的关系。当信号从区域一、二入射时，信号入射到稀疏面阵上的角度和信号入射到虚拟非均匀线阵上的角度之间的关系为：

当信号从区域三、四入射时，信号入射到稀疏面阵上的角度和信号入射到虚拟非均匀线阵上的角度之间的关系为：

其中

为信号入射到映射后的非均匀虚拟线阵上的角度，θ表示信号入射到稀疏面阵上的角度。信号入射到非均匀虚拟线阵上的角度

包含在峰值搜索得到的一系列峰值SP(i,j,k)中，将每一个峰值对应的角度与角度

进行对比就能得到正确估计的AoA(Arrival of Angle)、ToF(Time of Flight)和DFS(Doppler Frequency Shift)三个参数。

Claims (3)

1.一种基于稀疏面阵的三维参数估计方法，包括以下步骤：

步骤一：假设在Wi-Fi系统中有D个信号源，N个子载波，由于室内Wi-Fi单个AP业务天线数通常为4根，因此，接收端为由4根天线组成2×2的方形面阵，接收来自信号源的P个数据包中的CSI信息，天线间距为λ，不满足空间采样定理。将4根天线从上至下从左到右依次编号为天线1、2、3、4。θ_i表示信号源i∈{1,2,L,D}的入射角，不同信号源入射角度向量可以表示为Θ＝{θ₁,θ₂,L,θ_D}。

步骤二：将面阵所在平面按照信号入射方向进行分块，将信号入射角θ的大小分为0°～45°、45°～90°、90°～135°、135°～180°、180°～225°、225°～270°、270°～315°、315°～360°8个区域，分别对每个区域进行映射。

步骤三：对每一个区域分别沿着信号入射方向进行映射，并计算投影得到的非均匀虚拟阵天线之间的间距d。然后计算非均匀线阵的方向矢量

θ_i从-90°～90°以步长为1°变化，每变化一次，就得到一个稀疏面阵与非均匀虚拟线阵对应的相位差ΔΦ(θ_i),

将面阵方向矢量分别乘上面阵与线阵之间的相位差，得到虚拟线阵的方向矢量

步骤四：将得到的虚拟线阵的方向矢量构造入射信号

将构造的入射信号用三维参数估计联合估计算法进行三维参数估计。

步骤五：通过稀疏面阵和非均匀线阵的几何关系去除三维参数估计得到的伪峰值，得到正确的信号参数AoA、ToF、DFS。

2.根据权利要求1所述的一种基于稀疏面阵的三维参数估计方法，所述步骤三，计算虚拟线阵的方向矢量，包括以下步骤：

当信号从区域一入射到阵列时，将天线2和天线3沿着信号入射方向投影到天线1和天线4的连线上，投影之后的虚拟天线与天线1和天线4构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_subcarrier(t)＝[1 e^-j2πft L e^-j2π(N-1)ft]^T

Φ_array1(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差，Φ_subcarrier(t)表示信号由不同的子载波传输导致的相位差，Φ_packet(v)表示不同的数据包之间的相位差。

以天线1作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

由于入射角未知θ_i，而面阵和线阵之间的相位差ΔΦ₁(θ_i)与入射角有关，因此，我们让入射角θ_i从0°～45°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₁(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

当信号从区域二入射到阵列时，将天线2和天线3沿着信号入射方向投影到天线1和天线4的连线上，投影之后的虚拟天线与天线1和天线4构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_array2(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差。

以天线1作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

让入射角θ_i从45°～90°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₂(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

当信号从区域三入射到阵列时，将天线1和天线4沿着信号入射方向投影到天线2和天线3的连线上，投影之后的虚拟天线与天线2和天线3构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_array3(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差。

以天线3作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

让入射角θ_i从90°～135°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₃(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

当信号从区域四入射到阵列时，将天线1和天线4沿着信号入射方向投影到天线2和天线3的连线上，投影之后的虚拟天线与天线2和天线3构成一个非均匀的虚拟直线阵列。

构造的面阵方向矢量为：

其中

Φ_array4(θ_i)表示信号以角度θ_i入射到稀疏面阵上由天线之间的间距导致的相位差。

以天线3作为参考天线，计算投影得到的非均匀线虚拟阵天线之间的间距：

稀疏面阵映射得到的非均匀线阵的阵列流型为：

其中

表示信号以角度θ_i入射到虚拟非均匀线阵上由天线之间的间距导致的相位差，入射角为θ_i时，稀疏面阵和映射得到的非均匀线阵之间的相位差为：

其中

表示入射角为θ_i的信号入射到稀疏面阵和非均匀虚拟线阵之间的由于空间位置导致的相位差。

让入射角θ_i从135°～180°以步长为1°变化，θ_i每变化一次，就得到一个对应的相位差ΔΦ₄(θ_i)。将稀疏面阵的方向矢量与稀疏面阵和非均匀虚拟线阵之间的相位差相乘，得到虚拟非均匀线阵的方向矢量：

信号从区域5、6、7、8入射到整列上和信号从区域1、2、3、4入射到阵列上是相反的方向，映射方式完全相同。

3.根据权利要求1所述的一种基于稀疏面阵的三维参数估计方法，所述步骤五，通过稀疏面阵和非均匀线阵的几何关系去除三维参数估计得到的伪峰值，得到正确的信号参数AoA、ToF、DFS，包括以下步骤：

因为稀疏面阵的天线空间排列不满足空间采样定理，因此通过步骤四搜索得到的峰值中存在很多的伪峰值。只有当搜索的角度θ_i和面阵映射为线阵的角度相等时，得到的相位差ΔΦ(θ_i)才能正确地表达稀疏面阵和虚拟非均匀线阵之间的关系。此时，通过参数求解得到的结果才是正确的且等于θ_i。因此，在得到的峰值中，将每一个峰值通过映射的几何关系进行反解，我们就可以得出入射角θ_i。首先通过峰值搜索得到了一系列的峰值SP(i,j,k)，其中i，j，k分别代表峰值中多普勒速度、到达角和飞行时间对应的值。根据稀疏面阵和映射得到的虚拟非均匀线阵之间的几何关系，可以计算出信号入射到稀疏面阵上的角度和信号入射到虚拟非均匀线阵上的角度之间的关系。当信号从区域一、二入射时，信号入射到稀疏面阵上的角度和信号入射到虚拟非均匀线阵上的角度之间的关系为：

当信号从区域三、四入射时，信号入射到稀疏面阵上的角度和信号入射到虚拟非均匀线阵上的角度之间的关系为：

其中

为信号入射到映射后的非均匀虚拟线阵上的角度，θ表示信号入射到稀疏面阵上的角度。信号入射到非均匀虚拟线阵上的角度

包含在峰值搜索得到的一系列峰值SP(i,j,k)中，将每一个峰值对应的角度与角度

进行对比就能得到正确估计的AoA、ToF和DFS。

Images