无线定位参数估计方法、装置、计算机设备及存储介质

摘要

本申请涉及一种无线定位参数估计方法、装置、计算机设备和存储介质。该方法包括：对待定位设备发送的定位信号进行信道估计得到信道响应参数，对信道响应参数进行时延域超分辨谱估计得到时延域超分辨谱信息，并对时延域超分辨谱信息进行空间谱估计得到时延域超分辨谱‑空间谱信息，进而从该时延域超分辨谱‑空间谱信息中得到各传输路径对应的定位参数估计值。采用本方法能够先进行时延域超分辨谱估计再空间谱估计，通过级联方式得到二维谱信息并获取定位参数估计值，该方法避免了采用二维超分辨算法时所涉及到的高维度矩阵运算，从而降低了定位参数估计算法的复杂度，减少了算法的运算量，缩短了定位参数的估计周期，提高了定位参数的估计效率。

权利要求

1.一种无线定位参数估计方法，其特征在于，所述方法包括：对待定位设备发送的定位信号进行信道估计，得到信道响应参数；

对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息；

对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，所述时延域超分辨谱‑空间谱信息包括各传输路径对应的定位参数估计值；

其中，所述对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，包括：采用数字波束形成算法，对所述时延域超分辨谱信息进行空间谱估计，得到所述时延域超分辨谱‑空间谱信息。

2.根据权利要求1所述的方法，其特征在于，所述方法还包括：从所述时延域超分辨谱‑空间谱信息中确定直达径对应的目标定位参数估计值，所述目标定位参数估计值包括目标到达角度和目标传播时延。

3.根据权利要求2所述的方法，其特征在于，所述从所述时延域超分辨谱‑空间谱信息中确定直达径对应的所述目标定位参数估计值，包括：提取所述时延域超分辨谱‑空间谱信息中的谱峰，得到所述谱峰对应的衰落系数、所述定位信号的到达角度以及传播时延；

根据所述谱峰对应的衰落系数以及所述传播时延，对所述直达径进行识别，并输出所述直达径对应的所述目标到达角度以及所述目标传播时延。

4.根据权利要求1所述的方法，其特征在于，所述对待定位设备发送的定位信号进行信道估计，得到信道响应参数，包括：获取所述待定位设备发送的所述定位信号；

对所述定位信号进行时频转换，得到所述待定位设备的频域定位信号；

通过信道估计模型对所述频域定位信号进行信道估计，得到所述信道响应参数。

5.根据权利要求4所述的方法，其特征在于，所述信道估计模型包括所述信道响应参数与所述频域定位信号之间的映射关系。

6.根据权利要求1所述的方法，其特征在于，所述对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息，包括：采用加权最小二乘法对各信号接收通道的信道响应参数进行时延域超分辨谱估计，得到所述时延域超分辨谱信息。

7.一种无线定位参数估计装置，其特征在于，所述装置包括：信道估计模块，用于对待定位设备发送的定位信号进行信道估计，得到信道响应参数；

第一谱估计模块，用于对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息；

第二谱估计模块，用于对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，所述时延域超分辨谱‑空间谱信息包括各传输路径对应的定位参数估计值；

其中，所述第二谱估计模块具体用于采用数字波束形成算法，对所述时延域超分辨谱信息进行空间谱估计，得到所述时延域超分辨谱‑空间谱信息。

8.根据权利要求7所述的装置，其特征在于，所述装置还包括：目标估计值确定模块，其中：

所述目标估计值确定模块，用于从所述时延域超分辨谱‑空间谱信息中确定直达径对应的目标定位参数估计值。

9.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至6中任一项所述方法的步骤。

10.一种存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至6中任一项所述方法的步骤。

说明书

无线定位参数估计方法、装置、计算机设备及存储介质

技术领域

[0001] 本申请涉及无线定位技术领域，特别是涉及一种无线定位参数估计方法、装置、计算机设备及存储介质。

背景技术

[0002] 随着工业互联网、物联网和车联网的快速发展，高精度定位成为智能机器人、无人车等移动终端不可或缺的关键支撑服务。为了能够在受遮蔽的环境和室内环境提供更好的
导航定位服务，常用的定位技术包括：蜂窝网定位、无线局域网定位、蓝牙定位和超宽带定
位等；其中，这些定位技术大多采用了宽带发射信号以及阵列天线。而在室内定位系统中，
受限于移动终端尺寸，定位基站的天线阵孔径通常不会太大，空域分辨能力有限。因此，如
何提高定位系统在复杂环境中的分辨能力，成为复杂环境中定位系统的一核心问题。

[0003] 传统技术中，SpotFi无线局域网定位系统采用二维超分辨算法，对多径的到达角度（Angle of Arrival，AoA）和传播时延等定位参数进行估计，使得能够分辨多于阵元数的
多径信号，其中传播时延一般也称为到达时间（Time of Arrival，ToA）。但是，传统的定位
参数估计方式需要进行二维参数空间搜索，从而导致定位参数的估计效率较低。

发明内容

[0004] 基于此，有必要针对上述技术问题，提供一种能够提高定位参数估计效率的无线定位参数估计方法、装置、计算机设备及存储介质。

[0005] 一种无线定位参数估计方法，所述方法包括：

[0006] 对待定位设备发送的定位信号进行信道估计，得到信道响应参数；

[0007] 对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息；

[0008] 对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息；

[0009] 其中，所述时延域超分辨谱‑空间谱信息包括各传输路径对应的定位参数估计值。

[0010] 在其中一个实施例中，所述从所述时延域超分辨谱‑空间谱信息中确定目标定位参数估计值，包括：从所述时延域超分辨谱‑空间谱信息中确定直达径对应的所述目标定位
参数估计值，所述目标定位参数估计值包括目标到达角度和目标传播时延。

[0011] 在其中一个实施例中，所述从所述时延域超分辨谱‑空间谱信息中确定直达径对应的所述目标定位参数估计值，包括：

[0012] 提取所述时延域超分辨谱‑空间谱信息中的谱峰，得到所述谱峰对应的衰落系数、所述定位信号的到达角度以及传播时延；

[0013] 根据所述谱峰对应的衰落系数以及所述传播时延，对所述直达径进行识别，并输出所述直达径对应的所述目标到达角度以及所述目标传播时延。

[0014] 在其中一个实施例中，所述对待定位设备发送的定位信号进行信道估计，得到信道响应参数，包括：

[0015] 获取所述待定位设备发送的所述定位信号；

[0016] 对所述定位信号进行时频转换，得到所述待定位设备的频域定位信号；

[0017] 通过信道估计模型对所述频域定位信号进行信道估计，得到所述信道响应参数。

[0018] 在其中一个实施例中，所述信道估计模型包括所述信道响应参数与所述频域定位信号之间的映射关系。

[0019] 在其中一个实施例中，所述对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息，包括：采用加权最小二乘法对各信号接收通道的信道响应参数进行
时延域超分辨谱估计，得到所述时延域超分辨谱信息。

[0020] 在其中一个实施例中，所述对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，包括：采用数字波束形成算法，对所述时延域超分辨谱信息进
行空间谱估计，得到所述时延域超分辨谱‑空间谱信息。

[0021] 一种无线定位参数估计装置，所述装置包括：

[0022] 信道估计模块，用于对待定位设备发送的定位信号进行信道估计，得到信道响应参数；

[0023] 第一谱估计模块，用于对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息；

[0024] 第二谱估计模块，用于对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息；

[0025] 其中，所述时延域超分辨谱‑空间谱信息包括各传输路径对应的定位参数估计值。

[0026] 一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：

[0027] 对待定位设备发送的定位信号进行信道估计，得到信道响应参数；

[0028] 对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息；

[0029] 对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息；

[0030] 其中，所述时延域超分辨谱‑空间谱信息包括各传输路径对应的定位参数估计值。

[0031] 一种存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：

[0032] 对待定位设备发送的定位信号进行信道估计，得到信道响应参数；

[0033] 对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息；

[0034] 对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，所述时延域超分辨谱‑空间谱信息包括各传输路径对应的定位参数估计值。

[0035] 上述无线定位参数估计方法、装置、计算机设备和存储介质，接收站可以对待定位设备发送的定位信号进行信道估计，得到信道响应参数，对信道响应参数进行时延域超分
辨谱估计，得到时延域超分辨谱信息，并对时延域超分辨谱信息进行空间谱估计，得到时延
域超分辨谱‑空间谱信息，进而从该时延域超分辨谱‑空间谱信息中得到各传输路径对应的
定位参数估计值；该方法可以先进行时延域超分辨谱估计，再空间谱估计，通过这样的级联
方式得到二维谱信息，并获取定位参数估计值，该方法避免了采用二维超分辨算法执行并
行估计的处理方式时所涉及到的高维度矩阵运算，降低了定位参数估计算法的复杂度，减
少了算法的运算量，缩短了定位参数的估计周期，提高了定位参数的估计效率。

附图说明

[0036] 图1为一个实施例中无线定位参数估计方法的应用环境图；

[0037] 图2为一个实施例中无线定位参数估计方法的流程示意图；

[0038] 图3为一个实施例中获取信道响应参数的具体方法流程示意图；

[0039] 图4为另一个实施例中确定直达径对应的目标定位参数估计值的具体方法流程示意图；

[0040] 图5为另一个实施例中目标到达角度估计精度对比图；

[0041] 图6为另一个实施例中传播时延估计精度对比图；

[0042] 图7为另一个实施例中算法运行时间对比图；

[0043] 图8为一个实施例中无线定位参数估计装置的结构框图；

[0044] 图9为一个实施例中计算机设备的内部结构图。

具体实施方式

[0045] 为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不
用于限定本申请。

[0046] 本申请提供的无线定位参数估计方法，可以适用于图1所示的无线定位参数估计系统。如图1所示，该系统包括接收站和待定位设备。其中，接收站和待定位设备之间可以通
过无线电或者电磁波传播信号。可选的，上述接收站不限于4G/5G、无线局域网、超宽带定位
系统中已标定位置的接收设备，也就是4G/5G基站、无线局域网接入点以及超宽带锚点；上
述待定位设备可以为PC、便携式设备、服务器等具有数据处理功能的电子设备。本实施例可
以适用于单快拍场景，该单快拍场景可以理解为接收站在接收到单次定位信号时就能够实
现无线定位参数估计，且无需在估计过程中进行传统技术中的平滑操作就能够直接处理由
于多径传播造成的相干入射信号。需要说明的是，本实施例对待定位设备的具体形式不做
限定。

[0047] 在一个实施例中，如图2所示，提供了一种无线定位参数估计方法，以该方法应用于图1中的接收站为例进行说明，包括以下步骤：

[0048] S100、对待定位设备发送的定位信号进行信道估计，得到信道响应参数。

[0049] 具体的，接收站可以对待定位设备发送的定位信号进行信道估计。该定位信号可以为时域定位信号，也可以为频域定位信号。若上述定位信号为频域定位信号，则接收站可
以直接对该定位信号进行信道估计；若上述定位信号为时域定位信号，则接收站可以先对
定位信号进行预处理，然后对预处理后的定位信号进行信道估计。可选的，预处理可以为时
频转换处理，也可以为先对定位信号进行区间数据截取处理之后再进行时频转换处理，当
然，还可以为先对定位信号进行时频转换处理再进行区间数据截取处理等等。

[0050] 可以理解的是，接收站可以基于导频符号、基于判决反馈对定位信号进行信道估计，还可以基于训练序列和导频序列对定位信号进行信道估计，当然还可以通过其它方法
对定位信号进行信道估计。可选的，该信道估计实际上可以理解为从定位信号中将假定的
某个信道模型的模型参数估计出来的过程。其中，信道响应参数可以为信道响应矩阵。

[0051] 其中，为了减少定位信号的资源占用，待定位设备通常可以按照一定的图样将定位信号向天线阵对应的子带映射，此时，接收站可以根据映射关系直接提取相应子带的频
域接收信号，也就是，这种情况下接收站直接接收的是频域定位信号的情况。

[0052] S200、对信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息。

[0053] 具体的，接收站可以对信道响应参数进行时延域超分辨谱估计，得到多条传输路径对应的时延域超分辨谱信息。可选的，传输路径可以为定位信号从待定位设备传输至接
收站的路径，这些路径可以包括直达径和非直达径；该直达径可以理解为定位信号从待定
位设备传输至接收站的一条直线路径，非直达径可以有多条，每条非直达径可以理解为定
位信号从待定位设备传输至接收站多条不在同一直线上的路径组合在一起的一条曲折路
径。

[0054] S300、对时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，时延域超分辨谱‑空间谱信息包括各传输路径对应的定位参数估计值。

[0055] 具体的，接收站可以对获取到的时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，即二维谱信息。可选的，定位参数估计值可以包括定位信号的到
达角度、传播时延、到达时间差以及传输速率等等，但在本实施例中，上述定位参数估计值
包括定位信号的到达角度和传播时延。可选的，二维谱信息中包括直达径对应的定位参数
估计值，还包括非直达径对应的定位参数估计值。

[0056] 上述无线定位参数估计方法中，接收站可以对待定位设备发送的定位信号进行信道估计，得到信道响应参数，对信道响应参数进行时延域超分辨谱估计，得到时延域超分辨
谱信息，并对时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，进
而从该时延域超分辨谱‑空间谱信息中得到各传输路径对应的定位参数估计值；该方法可
以先进行时延域超分辨谱估计，再空间谱估计，通过这样的级联方式得到二维谱信息，并获
取定位参数估计值，该方法避免了采用二维超分辨算法执行并行估计的处理方式时所涉及
到的高维度矩阵运算，从而降低了定位参数估计算法的复杂度，减少了算法的运算量，从而
缩短了定位参数的估计周期，提高了定位参数的估计效率，进一步提高定位参数估计的实
时性。

[0057] 在一些场景中，时延域超分辨谱‑空间谱信息包含多组定位参数估计值，为了提高定位参数估计值的准确性，可以将其中一组定位参数估计值作为目标定位参数估计值，在
其中一个本实施例中，上述无线定位参数估计方法还可以包括：从时延域超分辨谱‑空间谱
信息中确定目标定位参数估计值，目标定位参数估计值包括：目标到达角度以及目标传播
时延。

[0058] 在本实施例，上述目标定位参数估计值可以包括目标到达角度以及目标传播时延。可选的，时延域超分辨谱‑空间谱信息可以通过信号谱图表示。接收站可以从时延域超
分辨谱‑空间谱信息中确定目标定位参数估计值，也就是可以将时延域超分辨谱‑空间谱信
息中直达径对应的定位参数估计值作为目标定位参数估计值，也可以将时延域超分辨谱‑
空间谱信息中非直达径中的任意一条非直达径对应的定位参数估计值作为目标定位参数
估计值。还可以理解为，接收站可以从时延域超分辨谱‑空间谱信息对应的信号谱图中选取
任意一点，并将该点对应的定位参数估计值作为目标定位参数估计值。

[0059] 其中，为了能够精确获取定位参数估计值，上述从时延域超分辨谱‑空间谱信息中确定目标定位参数估计值具体可以包括：从时延域超分辨谱‑空间谱信息中确定直达径对
应的目标定位参数估计值。

[0060] 在本实施例中，上述时延域超分辨谱‑空间谱信息可以通过一个二维谱图像展示。接收站可以遍历二维谱图中的每个点，判断当前点相对于其它所有邻接点的谱强度，若当
前点的谱强度大于其它所有邻接点的谱强度，则判定当前点为一个谱峰值点，按照这种方
式找出二维谱图的所有谱峰值点，之后按照谱峰强度进行排序，提取出其中最大的K1个谱
峰，这K1个谱峰对应为1个直达径和K1‑1个反射径对应的信号分量。接收站可以基于直达径
相比于其它非直达径（即反射径）的传播时间更短或者直达径比其它非直达径能量更强等
基本准则，从时延域超分辨谱空间谱信息中确定直达径对应的目标定位参数估计值，作为
最优的目标定位参数估计值。

[0061] 上述无线定位参数估计方法可以从时延域超分辨谱空间谱信息对应的信号谱图中确定直达径对应的目标定位参数估计值，从而提高了定位参数估计值的准确性。

[0062] 作为其中一个实施例，如图3所示，上述S100中对待定位设备发送的定位信号进行信道估计，得到信道响应参数的步骤，可以通过以下步骤实现：

[0063] S110、获取待定位设备发送的定位信号。

[0064] 具体的，接收站可以接收待定位设备发送的定位信号。在本实施例中，该定位信号可以为时域信号。

[0065] 在本实施例中，接收站可以通过天线阵的信号接收通道接收待定位设备发送的定位信号。可选的，该定位信号可以包括待定位设备的经纬度信息、方位角等等。可选的，上述
天线阵可以为点阵、面阵，还可以为线阵等，根据类型还可以为圆形天线阵、方形天线阵、菱
形天线阵等，本实施例对天线阵的类型本实施例不做任何限定。可选的，若天线阵共有N个
阵元，则每个阵元可以对应一个信号接收通道。

[0066] 若定位信号占用的子带数量为M，则天线阵从信号接收通道n接收的定位信号可以表示为向量 ， ，其中，Xm,n表示第n个信号接收通道第m个
子带接收的定位信号， 表示复数空间， 表示 维复数空间，也就是M维复向量空
间。可选的，本实施例中的所述向量均可以指列向量。其中，接收站的所有信号接收通道接
收到的定位信号矩阵可以表示为 ， ，若在M个子带上发送的定
位信号序列为 ， ，发送定位信号的中心载频为 ，定位信号的对应波长
为 （ ）， 表示真空中的光速。

[0067] S120、对定位信号进行时频转换，得到待定位设备的频域定位信号。

[0068] 在本实施例中，接收站可以对定位信号进行时频转换，得到待定位设备的频域定位信号。可选的，时频转换的方法可以为傅里叶变换，还可以为快速傅里叶变换，还可以为
傅里叶变换和快速傅里叶变换的组合变换等等。

[0069] 在本实施例中，上述时频转换方法可以为快速傅里叶变换，以减少算法的运算量，缩短定位参数的估计周期。为了不失一般性，假设M个子带呈均匀分布，分布间隔为 ，接
收的天线阵为等距线阵，阵元间距为d，且定位信号经由K条路径传播至等距线阵，第k条路
径的传播时延、方位角以及衰落系数分别表示为 、 和 ， 表示定位信号的入射方向
与等距线阵法线方向的夹角；在本实施例中，定位信号的传播时延可以表示定位信号的传
播距离，该传播时延与传播距离之间可以通过光速 相互转化，因此，多个信号接收通道接
收到的定位信号矩阵X可以表示为：

[0070]                  （1）；

[0071] 式（1）为经过快速傅里叶变换后得到的定位信号矩阵X的表达式，且式中S表示待定位设备发送的定位信号的对角矩阵，即 ， 可以表示以
向量中的每个元素作为主对角线元素获得对角矩阵； （ ）表示时延
域匹配矢量函数，该时延域匹配矢量函数的输入可以为传播时延 ，输出可以为M维传播时
延域匹配向量，表示时延域匹配矢量函数的作用域，且 为所有可能传输路径对应的传播
时延 的集合，即 ， 表示实数空间，传播时延域匹配向量的第m个元素可以表
示定位信号传播时延在第m个子带造成的相位偏移，从而，有
，表示虚数单位，定义为 ， （
）表示接收阵列导向矢量函数，接收阵列导向矢量函数的输入可以为定
位信号的到达角度 ，接收阵列导向矢量函数的输出可以为对应到达角度的阵列导向矢
量， 表示接收阵列导向矢量函数的作用域为 ，对应的值域为N维复数向量，为
所有可能的信号到达角度所组成的空间，即 ；若接收站的接收阵为等距线阵
时，函数值 的第n 个元素为 ， ，
为噪声矩阵，W中第m行第n列元素表示第n个信号接收通道第m个子带上的噪声
分量。

[0072] S130、通过信道估计模型对频域定位信号进行信道估计，得到信道响应参数。

[0073] 具体的，上述信道估计可以为盲估计法、半盲估计法，还可以为结合盲估计与基于训练序列估计法等等。

[0074] 其中，信道估计模型包括信道响应参数与频域定位信号之间的映射关系。

[0075] 在本实施例中，接收站可以通过预先设置的信道估计的模型，对频域定位信号进行信道估计，以得到信道的频域响应，进而为后续处理提供所需的信道状态信息。可选的，
信道估计的方法可以为最小二乘算法、最陡下降法或者最小均方误差法等等。在本实施例
中，上述信道估计的模型可以为基于参考信号的信道估计模型，具体信道估计的模型可以
自定义设置，包括信道响应参数与频域定位信号之间的映射关系。可选的，信道响应参数与
频域定位信号之间可以为正相关关系。

[0076] 但在本实施例中，上述信道估计模型包括：

[0077]                        （2）；

[0078] 式（2）中， 的第n列表示接收站的第n个信号接收通道的信道响应参数， 表示定位信号的对角矩阵， 表示频域定位信号矩阵。

[0079] 可以理解的是，假设已知接收端接收到的定位信号在频域上的具体形式，则上述信道估计可以采用最小二乘算法，得到的信道响应参数 具体表示为：

[0080]              （3）；

[0081] 式（3）中，（ ）的第n列表示第n个信号接收通道的信道响应参数； （）表示信道响应参数中的噪声分量， 的第m行第n列元素表示第n个
通道第m个子带信道估计值中的噪声分量。

[0082] 进一步地，上述对信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息的步骤，具体可以包括：采用加权最小二乘法对各信号接收通道的信道响应参数进行
时延域超分辨谱估计。

[0083] 具体的，接收站可以通过信道响应参数，对每个信号接收通道进行时延域超分辨谱估计。

[0084] 假设 （ ）表示信道响应参数H的第n列元素，也就是第n个信号接收通道的信道响应参数，则 可以表示为：

[0085]               （4）；

[0086] 式（4）中， 表示第n个接收阵元对第k条路径传输的定位信号的响应， 为向量 的第n个元素， （ ）表示第n个信号接收通道的噪声向量且为矩
阵 的第n列。

[0087] 对于定位信号的频域等间隔采样序列，不模糊的传播时延范围为 ，假设将该传播时延范围划分为P+1份，为了降低算法复杂度，通常设定 ，这P个扫描栅格点
对应的传播时延分别为 （ ）；记 为每个扫描栅格点上的衰落系数，当
时， ，在其它P‑K个扫描栅格点上， 。记 为
扫描栅格点集合上的衰落系数向量，且 （ ）表示扫描栅
格点集合上的传播时延匹配矩阵，则有第n个信号接收通道的信道响应参数可以表示为：

[0088]                     （5）。

[0089] 可以理解的是，接收站可以采用加权最小二乘法、基于矩阵特征空间分解法或者改进的矩阵特征空间分解法，对各信号接收通道的信道响应参数进行时延域超分辨谱估
计。在本实施例中，采用加权最小二乘法对各信号接收通道的信道响应参数进行时延域超
分辨谱估计可以表示为：

[0090]                    （6）；

[0091] 式（6）中 （ ）表示向量 的加权 范数， 表示第n个信号接收通道第 个扫描栅格点处的干扰协方差矩阵，该干扰由当前扫描栅格点 以外的信号分量
构成，即 可以表示为：

[0092]                   （7）；

[0093] 其中， 表示第n个信号接收通道信道响应参数的协方差矩阵，可以表示为， 表示取期望值。

[0094] 通常，为了提高分辨能力，降低谱旁瓣，可以通过对干扰协方差矩阵逆加权的最小二乘问题求解，以在估计当前扫描栅格点的时延域超分辨谱信息时，对其它扫描栅格点上
潜在的分量进行了抑制，该加权的最小二乘问题的解为：

[0095]                        （8）；

[0096] 式（8）中，由于 为未知量，因此无法直接计算；通常，可以采用迭代法迭代地对和 进行估计，每次计算时在公式（8）中代入上一次的迭代结果。其中，可以采用谐波
模型对 进行估计，具体方式如：

[0097]                           （9）；

[0098] 式（9）中的 表示功率估计值组成的矩阵，可以表示为：

[0099]                 （10）。

[0100] 现有的加权最小二乘问题的在迭代式求解过程中，涉及大量矩阵的乘法运算，导致运算量较大，本实施例利用传播时延匹配矩阵 的结构特点，可以采用快速傅里叶变换
对算法进行加速。详细步骤如下，针对每个信号接收通道n依次可以执行以下步骤，获取 ，
具体步骤如下：

[0101] 1）确定传播时延的扫描栅格点数P以及对应的传播时延栅格 ，为了保证精度，传播时延扫描间隔 一般小于定位系统固有传播时延分辨率的1/10，即
。为了使用快速傅里叶变换进行算法加速，P可以选为2的整数次
幂，即P可以表示为：

[0102]                         （11）；

[0103] 其中， 表示向上取整。

[0104] 2）令 ，可以使用快速傅里叶逆变换计算周期图谱，作为迭代初值：

[0105]                      （12）；

[0106] 其中， 表示对向量作P点快速傅里叶逆变换，记 表示第i次迭代的功率向量，其各元素表示向量 对应元素的模平方。

[0107] 3）可以使用快速傅里叶逆变换实现协方差矩阵第一行元素的快速计算，具体方法如 。

[0108] 4）利用矩阵 的范德蒙特性和 为对角矩阵的特性，可得协方差矩阵 为对角矩阵，因此，当前迭代步骤i的协方差矩阵 可以通过其第一行元素构建，得到
，其中， 表示取向量 的前 个元素， 表示
Toeplitz化算子，表示由向量产生相应的Toeplitz矩阵。

[0109] 5）计算协方差矩阵逆 ， 。

[0110] 6）与步骤7）共同完成谱值计算公式中的分子部分，首先计算 。

[0111] 7）利用矩阵 的范德蒙特性以及其元素为快速傅里叶变换因子的特性，计算。

[0112] 8）和步骤9）和步骤10）共同完成谱值计算公式中的分母部分，同样利用 的范德蒙性以及其各项为复指数的特性实现分母的快速计算，首先计算：

[0113]        （13）；

[0114] 9）将 排列为矩阵 ，则有

[0115] 。

[0116] 10）计算 。

[0117] 11）更新时延域超分辨谱估计值 ，其第p个元素为 。

[0118] 12）令 ，重复步骤3）至步骤11）直至 无明显改善，此时，令 。

[0119] 另外，上述对时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息具体可以包括：采用数字波束形成算法，对时延域超分辨谱信息进行空间谱估计，得到
时延域超分辨谱‑空间谱信息。

[0120] 具体的，接收站可以采用数字波束形成算法，对各个传播时延扫描栅格点上N个信号接收通道的时延域超分辨谱信息进行空间谱估计。接收站先进行时延域超分辨谱估计，
再进行空间谱估计，通过这样的级联方式得到时延域超分辨谱‑空间谱信息，即二维谱信
息，以避免采用二维超分辨算法执行并行估计的处理方式时所涉及到的高维度矩阵运算。

[0121] 其中，若第p个时延单元，N个信号接收通道的时延域超分辨谱信息对应的向量表示为 ，则有 ，进而对时延域超分辨谱信息对应的向量 依次进行
空间谱估计，具体可以通过以下步骤实现：

[0122] （1）接收站可以确定定位信号到达角度的扫描栅格点集合。假设使用等间隔 将天线阵覆盖的到达角度范围 划分为Q份，对应的到达角度集合为 。其中，天
线阵覆盖范围由天线阵的阵元方向图所确定，到达角度间隔 一般选择为阵列波束宽度的
1/20至1/10。记到达角度的扫描栅格点集合上的阵列流型矩阵为 。

[0123] （2）接收站可以对时延域超分辨谱估计值向量 依次进行空间谱估计，记时延域超分辨谱‑空间谱信息为 ， ，则有 ；

[0124] 记 ，其第q行第p列元素 表示对时延 到达角度 处的信道衰落系数的估计值。

[0125] 上述无线定位参数估计方法避免了采用二维超分辨算法执行并行估计的处理方式时所涉及到的高维度矩阵运算，同时，该方法在进行时延域超分辨谱估计的过程中采用
快速傅里叶变换，对谱求解的过程进行了加速，从而减少了算法的运算量，缩短了定位参数
的估计周期，进一步提高定位参数估计的实时性。

[0126] 在一些场景中，为了能够精确获取定位参数估计值，如图4所示，上述从时延域超分辨谱‑空间谱信息中确定直达径对应的目标定位参数估计值的步骤具体可以通过以下步
骤实现：

[0127] S400、提取时延域超分辨谱‑空间谱信息中的谱峰，得到谱峰对应的衰落系数、定位信号的到达角度以及传播时延。

[0128] 具体地，接收站可以将时延域超分辨谱‑空间谱信息通过一个二维谱图像展示，遍历二维谱图中的每个点，判断当前点相对于其它所有邻接点的谱强度，若当前点谱强度大
于其它所有邻接点的谱强度，则判定当前点为一个谱峰值点，按照这种方式找出二维谱图
的其它所有谱峰值点，之后，按照谱峰强度进行排序，提取出其中最大的K1个谱峰，这K1个谱
峰对应为1个直达径和K1‑1个反射径对应的信号分量，从而可以得到相应路径的衰落系数
、定位信号的到达角度 和传播时延 ，其中， 。

[0129] S500、根据谱峰对应的衰落系数以及传播时延，对直达径进行识别，并输出直达径对应的目标到达角度以及目标传播时延。

[0130] 具体的，接收站可以基于直达径相比于其它非直达径（即反射径）的传播时间更短或者直达径比其它非直达径能量更强等基本准则，根据K1个谱峰对应的衰落系数和传播时
延，确定K1个谱峰对应的直达径的时延域超分辨谱‑空间谱信息，从而通过获取到的直达径
的时延域超分辨谱‑空间谱信息得到直达径对应的目标到达角度和目标传播时延。

[0131] 示例性的，基于sub‑6G频段的5G系统的仿真平台对本实施例的方法进行验证，并和SpotFi算法结果对比，5G定位系统的仿真参数如下：

[0132]

[0133] SpotFi算法和本实施例中到达角度和传播时延估计误差的累计分布函数图分别如图5和图6所示，图5为到达角度估计精度对比图，图6为传播时延估计精度对比图，其中传
播时延已换算成距离，以68%分为点误差作为评估标准，SpotFi算法和本实施例中的到达角
度误差分别为：0.13°和0.32°；距离估计误差分别为：0.040m和0.077m。其中，如图7所示为
SpotFi算法和本实施例在每一次实验中的运行时间对比图，综合以上精度对比结果，可以
看到，本实施例以少量的精度损失为代价，运行时间相比于SpotFi节省了25倍以上，因此，
本实施例还非常适用于实时定位系统。

[0134] 上述无线定位参数估计方法可以从时延域超分辨谱‑空间谱信息对应的信号谱图中确定直达径对应的目标定位参数估计值，从而提高了定位参数估计值的准确性。

[0135] 应该理解的是，虽然图2‑4的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤
的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2‑4中的至少一
部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行
完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而
是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。

[0136] 在一个实施例中，如图8所示，提供了一种无线定位参数估计装置，包括：信道估计模块11、第一谱估计模块12和第二谱估计模块13，其中：

[0137] 信道估计模块11，用于对待定位设备发送的定位信号进行信道估计，得到信道响应参数；

[0138] 第一谱估计模块12，用于对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息；

[0139] 第二谱估计模块13，用于对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，所述时延域超分辨谱‑空间谱信息包括各传输路径对应的定位参
数估计值。

[0140] 本实施例提供的无线定位参数估计装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。

[0141] 在其中一个实施例中，无线定位参数估计装置还包括：目标估计值确定模块，其中：

[0142] 目标估计值确定模块，用于从所述时延域超分辨谱‑空间谱信息中确定直达径对应的所述目标定位参数估计值。

[0143] 本实施例提供的无线定位参数估计装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。

[0144] 在其中一个实施例中，目标估计值确定单元包括：信息提取子单元和识别子单元，其中：

[0145] 信息提取子单元，用于提取所述时延域超分辨谱‑空间谱信息中的谱峰，得到所述谱峰对应的衰落系数、所述定位信号的到达角度以及传播时延；

[0146] 识别子单元，用于根据所述谱峰对应的衰落系数以及所述传播时延，对所述直达径进行识别，并输出所述直达径对应的所述目标到达角度以及所述目标传播时延。

[0147] 本实施例提供的无线定位参数估计装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。

[0148] 在其中一个实施例中，信道估计模块11包括：定位信号获取单元、时频转换单元和信道估计单元，其中：

[0149] 定位信号获取单元，用于获取所述待定位设备发送的所述定位信号；

[0150] 时频转换单元，用于对所述定位信号进行时频转换，得到所述待定位设备的频域定位信号；

[0151] 信道估计单元，用于通过信道估计模型对所述频域定位信号进行信道估计，得到所述信道响应参数。

[0152] 其中，信道估计模型包括所述信道响应参数与所述频域定位信号之间的映射关系。

[0153] 本实施例提供的无线定位参数估计装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。

[0154] 在其中一个实施例中，第一谱估计模块12具体用于采用加权最小二乘法对各信号接收通道的信道响应参数进行时延域超分辨谱估计，得到所述时延域超分辨谱信息。

[0155] 本实施例提供的无线定位参数估计装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。

[0156] 在其中一个实施例中，第二谱估计模块13具体用于采用数字波束形成算法，对所述时延域超分辨谱信息进行空间谱估计，得到所述时延域超分辨谱‑空间谱信息。

[0157] 本实施例提供的无线定位参数估计装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。

[0158] 关于无线定位参数估计装置的具体限定可以参见上文中对于无线定位参数估计方法的限定，在此不再赘述。上述无线定位参数估计装置中的各个模块可全部或部分通过
软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处
理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上
各个模块对应的操作。

[0159] 在一个实施例中，提供了一种计算机设备，该计算机设备可以是服务器，其内部结构图可以如图9所示。该计算机设备包括通过系统总线连接的处理器、存储器和网络接口。
其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易
失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该
内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备
的数据库用于存储定位信号。该计算机设备的网络接口用于与外部的端通过网络连接通
信。该计算机程序被处理器执行时以实现一种无线定位参数估计方法。

[0160] 本领域技术人员可以理解，图9中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备
可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。

[0161] 在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：

[0162] 对待定位设备发送的定位信号进行信道估计，得到信道响应参数；

[0163] 对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息；

[0164] 对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，所述时延域超分辨谱‑空间谱信息包括各传输路径对应的定位参数估计值。

[0165] 在一个实施例中，提供了一种存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：

[0166] 对待定位设备发送的定位信号进行信道估计，得到信道响应参数；

[0167] 对所述信道响应参数进行时延域超分辨谱估计，得到时延域超分辨谱信息；

[0168] 对所述时延域超分辨谱信息进行空间谱估计，得到时延域超分辨谱‑空间谱信息，所述时延域超分辨谱‑空间谱信息包括各传输路径对应的定位参数估计值。

[0169] 本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机
可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，
本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可
包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器（Read‑
Only Memory，ROM）、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器
（Random Access Memory，RAM）或外部高速缓冲存储器。作为说明而非局限，RAM可以是多种
形式，比如静态随机存取存储器（Static Random Access Memory，SRAM）或动态随机存取存
储器（Dynamic Random Access Memory，DRAM）等。

[0170] 以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛
盾，都应当认为是本说明书记载的范围。

[0171] 以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来
说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护
范围。因此，本申请专利的保护范围应以所附权利要求为准。