CN113873430A - 一种基于时频分析的终端簇室内定位方法及定位系统 - Google Patents

Abstract

本发明公开了一种基于时频分析的终端簇室内定位方法及定位系统，该方法针对的基于终端簇的室内定位系统包含：室内接入节点T_α以及若干终端Q_p，p为终端编号，α为节点编号，且至少含有四个不在同一平面上的节点T₁、T₂、T₃和T₄，节点为待定位环境周围的基站，以终端Q₁表示区域内若干终端中任意一待定位终端，终端Q₁具有读写功能；该定位方法包含：(S100)节点以混合预编码方式发送信息；(S200)终端Q₁对接收到的来自节点T_α的信号进行检测；(S300)对终端Q₁进行三维空间位置估计。本发明的方法基于时频分析利用时间频率域来表征并且进行处理和分析信号的方法，摆脱了传统的傅里叶变换单一频域描述信号特征的局限性。

Description

一种基于时频分析的终端簇室内定位方法及定位系统

技术领域

本发明涉及一种室内定位方法，具体涉及一种基于时频分析的终端簇室内定位方法及定位系统。

背景技术

随着数据业务和多媒体业务的快速增加，人们对定位的需求日益增大，尤其在复杂的室内环境，如机场大厅、超市、图书馆、地下停车场等环境中，常常需要确定移动终端或其持有者、设施与物品在室内的位置信息。当前的大多数定位算法只是针对无线二维网络展开研究的，然而在实际应用中无线网络节点也会经常处于三维环境中，这些场合下需要提供移动终端的三维位置信息，目前研究人员提出了许多基于射频识别的室内定位解决方案。

基于典型的非测距DV-Hop定位算法具有成本低和对部署环境要求低的优势，但是定位过程中很大程度上地依赖网络中信标节点的分布，如果网络中锚节点分布不均匀，在估算未知节点坐标时会存在一定的定位误差。

发明内容

本发明的目的是提供一种基于时频分析的终端簇室内定位方法及定位系统，基于时频分析利用时间频率域来表征并且进行处理和分析信号的方法，摆脱了传统的傅里叶变换单一频域描述信号特征的局限性。

为了达到上述目的，本发明提供了一种基于时频分析的终端簇室内定位方法，该方法针对的基于终端簇的室内定位系统包含：室内接入节点T_α以及若干终端Q_p，p为终端编号，α为节点编号，且至少含有四个不在同一平面上的节点T₁、T₂、T₃和T₄，节点为待定位环境周围的基站，以终端Q₁表示区域内若干终端中任意一待定位终端，终端Q₁为具有读写功能的装置。

该定位方法包含：

(S100)节点以混合预编码方式发送信息：将任意一个节点T_α覆盖区域中的终端总数设为K，将该节点T_α覆盖区域中的终端基于信道时频分布划分为N个终端簇，设第q簇的终端数为K_q，q∈[1,2,…,N]，且

通过节点T₁,T₂,T₃,T₄端分别进行基于混合预编码的发射，节点T₁,T₂,T₃,T₄发送的信息包括节点的id和位置信息等；

(S200)终端Q₁对接收到的来自节点T_α的信号进行检测；

(S300)对终端Q₁进行三维空间位置估计：通过加权的距离向量-跳段方法求得终端Q₁与节点T₁、T₂、T₃和T₄之间的距离l₁、l₂、l₃和l₄，分别以节点T₁、T₂、T₃和T₄为圆心，距离l₁、l₂、l₃和l₄为半径得到四个球，实际的终端Q₁在这四个球所围成的空间区域中；设定四个节点T₁、T₂、T₃和T₄的坐标分别为(x₁,y₁,z₁)、(x₂,y₂,z₂)、(x₃,y₃,z₃)和(x₄,y₄,z₄)，建立三维空间球面方程组，采用最小二乘估计法解方程组，则终端Q₁的位置坐标为：

其中，上角标’表示求转置，上角标-1表示求逆。

在步骤(S100)中，节点以混合预编码方式发送信息的方法，包含：

(S110)对节点到终端的信道做乔-威廉姆斯变换得到CWD时频分布

为：

式中，μ为时间积分参数一，τ为时延积分参数二；β为衰减系数；t为工作时间；f为多普勒频率；

为节点T_α到第r簇第K_r个终端的信道；上角标#表示求共轭；j表示虚数单位；

(S120)基于CWD时频分布划分终端簇：设定第q簇终端簇的工作时间为(t_阈值下限,t_阈值上限)且对应的多普勒频率为(f_阈值下限,f_阈值上限)，根据终端p的CWD时频分布的三个坐标时间t_p、多普勒频率f_p和时频功率谱Power_p对终端进行分簇，为：

当时频功率谱Power_p>0，且

则终端p划为第q簇终端簇，记为终端簇

(S130)根据终端簇估计混合预编码：针对终端簇

采用几何平均分解获得的酉矩阵作为模拟预编码c_q,RF，采用规则化迫零构造数字预编码c_q,BB，根据模拟预编码和数字预编码得到混合预编码为：c_q＝c_q,RF·c_q,BB；

(S140)基于终端分簇的混合预编码发射：节点T_α的原始信号S通过混合预编码进行加权，并映射到相应的天线端口，发射信号为：c_q·s。

在步骤(S300)中，估计所述终端Q₁到各节点的距离的方法，包含：

(S321)获得终端Q₁与节点T₁、T₂、T₃和T₄的最小跳数，为：

其中，

是广义节点T_j+的信号强度因子，取为该广义节点T_j+与上一跳广义节点T_j间的信道矩阵的谱范数的平方：

j为广义节点号，j+为表示节点T_j的相邻上一跳节点的节点号，广义节点包括节点和伪节点，伪节点为除终端Q₁外的其它待定位终端；INT_α为终端Q₁到节点T_α的最小整数跳数；

(S322)估计距离终端Q₁最近节点T_i的加权平均跳距HopSizeavg_Q，为：

其中，HopSize_Q为采用均方根误差法求得的距离终端Q₁最近节点T_i的平均跳距；ζ_α为加权系数；

(S323)计算终端Q₁到四个节点T_α的距离l_α，为：

l_α＝Hop_αQ·HopSizeavg_Q。

优选地，在步骤(S130)中，根据终端簇估计混合预编码的方法，包含：

(S131)估计模拟预编码c_q,RF：针对第q簇终端簇

设定其聚合信道为

对聚合信道

进行几何平均分解，为：

将酉矩阵Γ作为模拟预编码，c_q,RF＝Γ；其中，U、Γ为酉矩阵，满足U^*U＝Γ^*Γ＝I_η×η，上角标*表示求共轭转置，η是聚合信道

的秩，

I_η×η是η行η列的单位阵；Λ为具有相等对角元素的实上三角矩阵，且其对角元素为

λ₁≥λ₂≥…≥λ_η≥0是相关矩阵

的非零特征值，所有对角元素值都等于

特征值的几何均值，分解能够获得具有相等增益的子信道；

(S132)估计数字预编码c_q,BB：采用规则化迫零构造数字预编码，为：

式(8)中，I为单位阵；P_α为节点T_α的发射功率；上角标*表示求共轭转置；K为节点T_α覆盖区域中的终端总数；δ_q为功率归一化因子，且

其中，V_q为中间变量，

上角标-1表示求逆；

(S133)根据模拟预编码和数字预编码得到混合预编码，为：c_q＝c_q,RF·c_q,BB。

优选地，在步骤(S322)中，采用均方根误差法求得的距离终端Q₁最近节点T_i的平均跳距为：

其中，Hop_αi为节点T_α到节点T_i的最小跳数，α≠i。

优选地，在步骤(S322)中，加权系数ζ_α为：

优选地，在步骤(S110)中，0.1<β<10。

本发明的另一目的是提供一种基于时频分析的终端簇室内定位系统，该室内定位系统包含：室内接入节点T_α以及若干终端Q_p，p为终端编号，α为节点编号，且至少含有四个不在同一平面上的节点T₁、T₂、T₃和T₄，节点为待定位环境周围的基站，以终端Q₁表示区域内若干终端中任意一待定位终端，终端Q₁为具有读写功能的装置；该室内定位系统中任意一待定位终端Q₁通过所述的方法进行定位。

本发明的基于时频分析的终端簇室内定位方法及定位系统，具有以下优点：

本发明的方法基于时频分析利用时间频率域来表征并且进行处理和分析信号的方法，摆脱了传统的傅里叶变换单一频域描述信号特征的局限性。时频域的分析比单一时域的信号冲击响应或单一频域的信号功率谱密度更为有效。通常把时域和频域结合起来描述和观察信号的时频特征，构成信号的时频谱图，通过时频分析算法分析信号的频谱含量怎样随时间变化。

此外，采用基于几何平均分解的混合预编码，可以将信道矩阵分解为具有相同对角元素值的信道，所有对角元素值均等于信道矩阵特征值的几何均值，使得子信道获得等增益，避免复杂的比特分配。

基于加权的距离向量-跳段测距法，通过引入信号强度因子对最小跳数及平均跳距加权，减小了常规距离向量-跳段中最小跳数和平均跳距带来的误差，从而定位更为精确。

附图说明

图1为本发明的基于时频分析的终端簇室内定位方法的流程图。

图2为本发明节点以混合预编码方式发送信息的方法的流程图。

图3为本发明方法的定位框图。

图4为本发明基于信道时频分布划分终端簇的示意图。

图5为本发明时频分布等高线图。

图6为本发明时频分布图。

图7为本发明基于终端分簇的混合预编码发射的示意图。

图8为本发明加权的DV-Hop测距方法示意图。

具体实施方式

下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例1

一种基于时频分析的终端簇室内定位方法，如图3示，为本发明方法的定位框图，该方法针对的基于终端簇的室内定位系统包含：室内接入节点T_α以及若干终端Q_p，p为终端编号，α为节点编号，且至少含有四个不在同一平面上的节点T₁、T₂、T₃和T₄，节点为待定位环境周围的基站，以终端Q₁表示区域内若干终端中任意一待定位终端，终端Q₁为具有读写功能的装置。如图1所示，为本发明的基于时频分析的终端簇室内定位方法的流程图，该定位方法包含：

(S100)节点以混合预编码方式发送信息：将任意一个节点T_α覆盖区域中的终端总数设为K，将该节点T_α覆盖区域中的终端基于信道时频分布划分为N个终端簇，设第q簇的终端数为K_q，q∈[1,2,…,N]，且

通过节点T₁,T₂,T₃,T₄端分别进行基于混合预编码的发射，节点T₁,T₂,T₃,T₄发送的信息包括节点的id和位置信息等；

(S200)终端Q₁对接收到的来自节点T_α的信号进行检测：在实际检测时，假设终端Q₁接收到T₁、T₂、T₃和T₄四个节点的信号，从中分离出各个信号，然后进行处理；

(S300)对终端Q₁进行三维空间位置估计：通过加权的距离向量-跳段方法求得终端Q₁与各节点T_α(α＝1,2,3,4)之间的距离l_α(α＝1,2,3,4)，分别以各节点T_α(α＝1,2,3,4)为圆心，距离l_α(α＝1,2,3,4)为半径得到四个球，实际的终端Q₁在这四个球所围成的空间区域中；再采用最小二乘估计法解方程组，则终端Q₁的位置坐标为：

其中，上角标’表示求转置，上角标-1表示求逆。

在步骤(S100)中，节点以混合预编码方式发送信息的方法，参见图2，包含：

(S110)对节点到终端的信道做乔-威廉姆斯变换得到CWD时频分布。

如图4-6所示，图4为本发明基于信道时频分布划分终端簇的示意图，根据终端与节点的无线传输信道的时频分布对终端进行分簇。

在室内立体停车库中，有步行的人流手持的终端Q₁、Q₂、Q₃(工作时间为(t_Ⅱ,t_Ⅲ)，对应的多普勒频率为f_Ⅰ)、层间升降的车辆乘载的终端Q₄、Q₅(工作时间为(t_Ⅰ,t_Ⅱ)，对应的多普勒频率为f_Ⅱ)和地面行驶的车辆乘载的终端Q₆、Q₇、Q₈和Q₉(工作时间为(0,t_Ⅲ)，对应的多普勒频率为f_Ⅲ)等，上述三种情形因工作时间、移动速度和方向的不同，使得终端具有不同的时频功率谱。可以将步行的人流、地面行驶的车辆和层间升降的车辆承载的终端划分为I、II和III等三个簇。图4中终端总数K＝9，终端簇数Q＝3，第I簇、第II簇和第III簇终端数分别为3个、2个和4个。

本发明终端簇中终端的选择方法，假设终端p的工作时间为t_p，多普勒频率为f_p，时频功率谱为Power_p，这里p＝1,2,3,…9。

参见图5、图6，分别为乔-威廉姆斯(Choi-Williams)时频分布等高线图和时频分布图。假设H_α,r1为节点T_α到第r簇(r≠q)第1个终端的信道，

为节点T_α到第r簇第2个终端的信道，…，

为节点T_α到第r簇第K_r个终端的信道；即对信道H_α,r1做乔-威廉姆斯变换，得到CWD时频分布，CWD时频分布揭示了信道的时频域功率谱，能够清晰的分辨信道频谱在具体时间点的变化过程。CWD时频分布为：

式中，μ为时间积分参数一，τ为时延积分参数二；β为衰减系数，β越小，对交叉项的抑制能力越大，但β的取值要在分辨率和交叉项抑制之间取折中，一般0.1<β<10；t为工作时间；f为多普勒频率；

为节点T_α到第r簇第K_r个终端的信道；上角标#表示求共轭；j表示虚数单位，

(S120)基于CWD时频分布划分终端簇：根据CWD时频分布的三个坐标时间t、多普勒频率f和时频功率谱Power_p对终端进行分簇，具体为：

(1)当时频功率谱Power_p>0，且

p＝1,2,3，则划为第I类终端簇，记为

(2)当Power_p>0，且

p＝4,5，则划为第II类终端簇，记为

(3)当Power_p>0，且

p＝6,7,8,9，则划为第III类终端簇，记为

其中，频域阈值f_Ⅰ、f_Ⅱ和f_Ⅲ可根据具体应用场景进行选取。例如，室内停车库场景，终端传输信号频率f_z＝2GHz＝2×10⁹Hz，光速c＝3×10⁸m/s，人员手持终端步行速度假设为3.6km/h，即v＝3.6×1000/3600m/s，节点与终端的连线与速度方向的夹角θ＝0，根据多普勒频率公式，则

车辆垂直升降速度为10km/h，节点与终端的连线与速度方向的夹角θ＝π/4，则f_Ⅱ＝13.1Hz；车辆水平行进速度为15km/h，节点与终端的连线与速度方向的夹角θ＝0，则f_Ⅲ＝27.9Hz。

时域阈值t_Ⅰ、t_Ⅱ和t_Ⅲ可根据终端实际工作时间进行选取。

(S130)根据终端簇估计预编码

针对终端簇，利用GMD分解进行模拟预编码和数字预编码的混合预编码设计。通过信道矩阵的GMD分解获得的酉矩阵做为模拟预编码，用于抑制终端簇之间的干扰；通过迫零处理设计的数字预编码，用于降低终端簇内干扰。具体包含：

(S131)估计模拟预编码：

针对第q簇终端，假设

为节点T_α到第q簇中所有终端之间的聚合信道；H_α,q1为节点T_α到第q簇第1个终端的信道，H_α,q2为节点T_α到第q簇第2个终端的信道，…，H_α,qKq为节点T_α到第q簇第K_q个终端的信道。

对聚合信道

进行几何平均分解，即

其中U、Γ为酉矩阵，即满足U^*U＝Γ^*Γ＝I_η×η，上角标*表示求共轭转置，η是聚合信道

的秩，即

I_η×η是η行η列的单位阵；Λ为具有相等对角元素的实上三角矩阵，且其对角元素为

λ₁≥λ₂≥…≥λ_η≥0是相关矩阵

的非零特征值。所有对角元素值都等于

特征值的几何均值，从而分解能够获得具有相等增益的子信道。

则，选择酉矩阵Γ作为模拟预编码，即：c_q,RF＝Γ。

(S132)估计数字预编码

设聚合信道

与第q簇终端模拟预编码构成的等效信道为：

采用规则化迫零构造内层预编码，则：

式(8)中，I为单位阵；P_α为节点T_α的发射功率；上角标*表示求共轭转置；K为节点T_α覆盖区域中的终端总数；δ_q为功率归一化因子，且

其中，V_q为中间变量，

上角标-1表示求逆。

(S133)根据模拟预编码和数字预编码得到混合预编码，为：

c_q＝c_q,RF·c_q,BB (9)

同理，得到对应于终端簇

和

的预编码分别为：c₂和c₃。

(S140)基于终端分簇的混合预编码发射

如图7所示，为基于终端分簇的混合预编码发射的示意图，对于终端簇

分别采用混合预编码方式发送信息，即采用对应的预编码c₁、c₂和c₃，具体为节点T₁的原始信号S通过混合预编码进行加权，并映射到相应的天线端口，即发射信号为：c₁·s，c₂·s，……，c₃·s，图中为避免图形混淆仅给出了终端簇

的发射图例。

(S300)终端进行三维空间位置估计

如图8所示，为加权的DV-Hop测距方法示意图，设终端Q₁的坐标为(x,y,z)，四个节点T₁、T₂、T₃和T₄的坐标已知分别为(x₁,y₁,z₁)、(x₂,y₂,z₂)、(x₃,y₃,z₃)和(x₄,y₄,z₄)，且终端周围有伪节点T_w1,T_w2,…,T_w10(如其它待定位终端)，伪节点能够进行信息传输但位置坐标未知。则终端进行三维空间位置估计的方法，包含：

(S310)验证四个节点T_α(α＝1,2,3,4)的坐标不在同一平面上

四个坐标确定唯一的三棱锥，进而确定出唯一的移动终端位置坐标；实际上，四个节点T₁、T₂、T₃和T₄是预先布置的且不在同一个平面上。

(S320)利用加权的DV-Hop方法估计终端Q₁到各节点的距离，具体包含：

(S321)获得终端与所有节点间的最小跳数

利用信号强度因子对相邻节广义点间(包括节点、伪节点)的跳数进行加权修正，即将与广义节点直接通信的相邻广义节点的信号强度因子作为基准，将第一跳记为1，并在广义节点传输的信息包中加入基准信号强度因子。其余广义节点接收到该带有基准信号强度因子的信息包后用自身广义节点的信号强度因子与基准信号强度因子作比值处理，并将上一跳的跳数与加权跳数的和作为本广义节点的跳数，设终端Q₁到四个节点T₁、T₂、T₃和T₄的最小跳数分别为Hop_1Q、Hop_2Q、Hop_3Q和Hop_4Q，则Hop_1Q、Hop_2Q、Hop_3Q和Hop_4Q具体如下：

其中，ρ_TjTj+(j为广义节点号)是广义节点T_j+的信号强度因子，取为该广义节点T_j+与上一跳广义节点T_j间的信道矩阵的谱范数的平方：

例如，伪节点T_w1与相邻的上一跳节点T₁间的信道矩阵

则伪节点T_w1的信号强度因子为：

(S322)估计终端的加权平均跳距

已知各节点的坐标信息，则节点T₃与节点T₁、T₂和T₄间的距离分别为：d₁₃、d₂₃和d₄₃，终端Q₁仅记录距离其最近节点的平均跳距，此时节点T₃距离终端Q₁为最近的节点，因此将节点T₃的平均跳距作为Q₁的平均跳距，采用均方根误差法求得平均跳距：

其中，Hop_α3为节点T_α(α≠3)到节点T₃的最小跳数。

加权系数，为：

则，Q₁的加权平均跳距：

(S323)计算终端与节点间的距离

则终端Q₁到四个节点T₁、T₂、T₃和T₄的距离分别为：

l₁＝Hop_1Q·HopSizeavg_Q

l₂＝Hop_2Q·HopSizeavg_Q

l₃＝Hop_3Q·HopSizeavg_Q

l₄＝Hop_4Q·HopSizeavg_Q

(S330)建立三维空间球面方程组并求解

三维空间以节点T₁、T₂、T₃和T₄的坐标(x₁,y₁,z₁)、(x₂,y₂,z₂)、(x₃,y₃,z₃)和(x₄,y₄,z₄)为球心，分别以l₁、l₂、l₃和l₄为半径的球面方程：

根据极大似然估计法求解上述(10)-(13)方程组，得到终端Q₁的三维坐标的估计值为：

其中，上角标’表示求转置，上角标-1表示求逆，

尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。

Claims (6)

1.一种基于时频分析的终端簇室内定位方法，其特征在于，该方法针对的基于终端簇的室内定位系统包含：室内接入节点T_α以及若干终端Q_p，p为终端编号，α为节点编号，且至少含有四个不在同一平面上的节点T₁、T₂、T₃和T₄，节点为待定位环境周围的基站，以终端Q₁表示区域内若干终端中任意一待定位终端，终端Q₁为具有读写功能的装置；

该定位方法包含：

(S100)节点以混合预编码方式发送信息：将任意一个节点T_α覆盖区域中的终端总数设为K，将该节点T_α覆盖区域中的终端基于信道时频分布划分为N个终端簇，设第q簇的终端数为K_q，q∈[1,2,…,N]，且

通过节点T₁,T₂,T₃,T₄端分别进行基于混合预编码的发射，节点T₁,T₂,T₃,T₄发送的信息包括节点的id和位置信息等；

(S200)终端Q₁对接收到的来自节点T_α的信号进行检测；

(S300)对终端Q₁进行三维空间位置估计：通过加权的距离向量-跳段方法求得终端Q₁与节点T₁、T₂、T₃和T₄之间的距离l₁、l₂、l₃和l₄，分别以节点T₁、T₂、T₃和T₄为圆心，距离l₁、l₂、l₃和l₄为半径得到四个球，实际的终端Q₁在这四个球所围成的空间区域中；设定四个节点T₁、T₂、T₃和T₄的坐标分别为(x₁,y₁,z₁)、(x₂,y₂,z₂)、(x₃,y₃,z₃)和(x₄,y₄,z₄)，建立三维空间球面方程组，采用最小二乘估计法解方程组，则终端Q₁的位置坐标为：

其中，上角标’表示求转置，上角标-1表示求逆。

在步骤(S100)中，节点以混合预编码方式发送信息的方法，包含：

(S110)对节点到终端的信道做乔-威廉姆斯变换得到CWD时频分布

为：

式中，μ为时间积分参数一，τ为时延积分参数二；β为衰减系数；t为工作时间；f为多普勒频率；

为节点T_α到第r簇第K_r个终端的信道；上角标#表示求共轭；j表示虚数单位；

(S120)基于CWD时频分布划分终端簇：设定第q簇终端簇的工作时间为(t_阈值下限,t_阈值上限)且对应的多普勒频率为(f_阈值下限,f_阈值上限)，根据终端p的CWD时频分布的三个坐标时间t_p、多普勒频率f_p和时频功率谱Power_p对终端进行分簇，为：

当时频功率谱Power_p>0，且

则终端p划为第q簇终端簇，记为终端簇

(S130)根据终端簇估计混合预编码：针对终端簇

采用几何平均分解获得的酉矩阵作为模拟预编码c_q,RF，采用规则化迫零构造数字预编码c_q,BB，根据模拟预编码和数字预编码得到混合预编码为：c_q＝c_q,RF·c_q,BB；

(S140)基于终端分簇的混合预编码发射：节点T_α的原始信号S通过混合预编码进行加权，并映射到相应的天线端口，发射信号为：c_q·s；

在步骤(S300)中，估计所述终端Q₁到各节点的距离的方法，包含：

(S321)获得终端Q₁与节点T₁、T₂、T₃和T₄的最小跳数，为：

其中，

是广义节点T_j+的信号强度因子，取为该广义节点T_j+与上一跳广义节点T_j间的信道矩阵的谱范数的平方：

j为广义节点号，j+为表示节点T_j的相邻上一跳节点的节点号，广义节点包括节点和伪节点，伪节点为除终端Q₁外的其它待定位终端；INT_α为终端Q₁到节点T_α的最小整数跳数；

(S322)估计距离终端Q₁最近节点T_i的加权平均跳距HopSizeavg_Q，为：

其中，HopSize_Q为采用均方根误差法求得的距离终端Q₁最近节点T_i的平均跳距；ζ_α为加权系数；

(S323)计算终端Q₁到四个节点T_α的距离l_α，为：

l_α＝Hop_αQ·HopSizeavg_Q。

2.根据权利要求1所述的基于时频分析的终端簇室内定位方法，其特征在于，在步骤(S130)中，根据终端簇估计混合预编码的方法，包含：

(S131)估计模拟预编码c_q,RF：针对第q簇终端簇

设定其聚合信道为

对聚合信道

进行几何平均分解，为：

将酉矩阵Γ作为模拟预编码，c_q,RF＝Γ；其中，U、Γ为酉矩阵，满足U^*U＝Γ^*Γ＝I_η×η，上角标*表示求共轭转置，η是聚合信道

的秩，

I_η×η是η行η列的单位阵；Λ为具有相等对角元素的实上三角矩阵，且其对角元素为

λ₁≥λ₂≥…≥λ_η≥0是相关矩阵

的非零特征值，所有对角元素值都等于

特征值的几何均值，分解能够获得具有相等增益的子信道；

(S132)估计数字预编码c_q,BB：采用规则化迫零构造数字预编码，为：

式(8)中，I为单位阵；P_α为节点T_α的发射功率；上角标*表示求共轭转置；K为节点T_α覆盖区域中的终端总数；δ_q为功率归一化因子，且

其中，V_q为中间变量，

上角标-1表示求逆；

(S133)根据模拟预编码和数字预编码得到混合预编码，为：c_q＝c_q,RF·c_q,BB。

3.根据权利要求1所述的基于时频分析的终端簇室内定位方法，其特征在于，在步骤(S322)中，采用均方根误差法求得的距离终端Q₁最近节点T_i的平均跳距为：

其中，Hop_αi为节点T_α到节点T_i的最小跳数，α≠i。

4.根据权利要求1所述的基于时频分析的终端簇室内定位方法，其特征在于，在步骤(S322)中，加权系数ζ_α为：

5.根据权利要求1所述的基于时频分析的终端簇室内定位方法，其特征在于，在步骤(S110)中，0.1<β<10。

6.一种基于时频分析的终端簇室内定位系统，其特征在于，该室内定位系统包含：室内接入节点T_α以及若干终端Q_p，p为终端编号，α为节点编号，且至少含有四个不在同一平面上的节点T₁、T₂、T₃和T₄，节点为待定位环境周围的基站，以终端Q₁表示区域内若干终端中任意一待定位终端，终端Q₁为具有读写功能的装置；该室内定位系统中任意一待定位终端Q₁通过如权利要求1-5中任意一项所述的方法进行定位。

Images