CN114755629B - 提高室内定位精度的方法及装置 - Google Patents
提高室内定位精度的方法及装置 Download PDFInfo
- Publication number
- CN114755629B CN114755629B CN202210410379.8A CN202210410379A CN114755629B CN 114755629 B CN114755629 B CN 114755629B CN 202210410379 A CN202210410379 A CN 202210410379A CN 114755629 B CN114755629 B CN 114755629B
- Authority
- CN
- China
- Prior art keywords
- signal data
- data
- phase
- fitting
- fitted
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 239000011159 matrix material Substances 0.000 claims abstract description 44
- 238000012545 processing Methods 0.000 claims abstract description 19
- 238000001914 filtration Methods 0.000 claims abstract description 15
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 24
- 238000010586 diagram Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radio Transmission System (AREA)
Abstract
本发明实施例的方法中提出了一种提高室内定位精度的方法及装置，在收到由接收端天线阵列连续发送的信号数据后，从中获取信号样本进行构建以生成原始信号向量；然后，使用拟合函数对原始信号向量中的I通道信号数据和Q通道信号数据进行拟合后，输出拟合相位数据；随后，基于拟合相位数据计算生成输入矩阵，并将输入矩阵代入迭代模型收敛后生成信号数据角度估计值；最后，基于输入矩阵进行滤波处理后生成信号数据相位估计值。本申请通过结合非线性递归最小二乘法和卡尔曼滤波器来提高到达角的估计精度。非线性递归最小二乘法可以消除信号数据处理过程中的多径和路径损耗效应，卡尔曼滤波器可以消除接收信号中天线切换引起的误差影响。
Description
【技术领域】
本发明涉及定位技术领域，尤其涉及一种提高室内定位精度的方法及装置。
【背景技术】
BLE(Bluetooth Low Energy，低功耗蓝牙)定位具有低成本和低功耗的优点，目前已广泛应用于诸多领域。目前，在室内定位技术领域广泛采用的算有接收强度算法、到达角度算法、到达时间算法和到达时间差算法。其中又以到达角度算法最为常用，但到达角度算法具方法对噪声以及多径和路径损耗效应极为敏感，所以关于如何提高定位精度还有许多技术难点亟待解决。
【发明内容】
有鉴于此，本发明实施例提供了一种提高室内定位精度的方法及装置。
第一方面，本发明实施例提供了一种提高室内定位精度的方法，所述方法包括：
S1、收到由接收端天线阵列连续发送的信号数据后，从中获取信号样本进行构建以生成原始信号向量；
S2、使用拟合函数对原始信号向量中的I通道信号数据和Q通道信号数据进行拟合后，输出拟合相位数据；
S3、基于拟合相位数据计算生成输入矩阵，并将输入矩阵代入迭代模型收敛后生成信号数据角度估计值；
S4、基于输入矩阵进行滤波处理后生成信号数据相位估计值。
如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述S1具体包括：
接收端天线阵列包含Ne个天线单元，其对应的信号数据的信号向量表示为：
其中，n为样本编号，/>
为样本n中第k个天线对应的信号数据；
将信号向量
表示为I通道信号数据分量和Q通道信号数据分量的形式：
其中，fc为载波频率，T为采样间隔，/>
为同相载波的幅度，/>
为正交相位载波的幅度。
如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述S2具体包括：
通过给定变量的拟合函数对I通道信号数据和Q通道信号数据进行拟合，其中，I通道信号数据通过如下函数进行拟合：
Q通道通过如下函数进行拟合：/>
/>
通过公式
计算出I通道信号数据中样本n中第k个天线对应的相位/>
通过公式
计算出样本n的加权相位，从而输出拟合相位数据，其中，/>
为样本n中第k个天线对应的方位角，/>
为标准方位角。
如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述基于拟合相位数据计算生成输入矩阵，具体包括：
根据拟合到达角数据通过公式
计算出拟合到达角θn，其中，λ为信号波长，相位差Δφn＝φn-φn-1，d为相邻天线的间距；
根据所述拟合相位数据和所述拟合到达角数据构建输入矩阵，输入矩阵表示为：
其中，
Θ＝[θ1,θ2,…,θn]T，/>
如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述输入矩阵代入迭代模型收敛后生成信号数据估计值，具体包括：
构建相互关联的模型分量
θn＝θn-1+ηnπn、/>
通过
迭代计算最小二乘解直至/>
时收敛,其中，Ψf为最小二乘解，/>
为经验阈值，f＝1,2,...N，下标f表示前f个数据，即从1-f的数据；
根据模型分量计算出收敛时的信号数据角度估计值θ。
如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述S4具体包括：
建立状态模型：
其中
为当前状态的相位估计值，/>
为上一状态的相位估计值，K(k)为凯尔曼增益，Z(k)为最新的测量值，/>
为根据先验参数计算的预测值；
将状态模型拟合至卡尔曼滤波器后进行滤波处理后得到信号数据相位估计值。
第二方面，本发明实施例提供了一种提高室内定位精度的装置，所述装置包括：
预处理模块，用于在收到由接收端天线阵列连续发送的信号数据后，从中获取信号样本进行构建以生成原始信号向量；
拟合模块，用于使用拟合函数对原始信号向量中的I通道信号数据和Q通道信号数据进行拟合后，输出拟合相位数据；
计算模块，用于基于拟合相位数据计算生成输入矩阵，并将输入矩阵代入迭代模型收敛后生成信号数据角度估计值；
滤波模块，用于基于输入矩阵进行滤波处理后生成信号数据相位估计值。
如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述预处理模块具体用于：
接收端天线阵列包含Ne个天线单元，其对应的信号数据的信号向量表示为：
其中，n为样本编号，/>
为样本n中第k个天线对应的信号数据；
将信号向量
表示为I通道信号数据分量和Q通道信号数据分量的形式：
其中，fc为载波频率，T为采样间隔，/>
为同相载波的幅度，/>
为正交相位载波的幅度。
如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述拟合模块具体用于：
通过给定变量的拟合函数对I通道信号数据和Q通道信号数据进行拟合，其中，I通道信号数据通过如下函数进行拟合：
Q通道通过如下函数进行拟合：/>
通过公式
计算出I通道信号数据中样本n中第k个天线对应的相位/>
通过公式
计算出样本n的加权相位，从而输出拟合相位数据，其中，/>
为样本n中第k个天线对应的方位角，/>
为标准方位角。
如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述计算模块具体用于：
根据拟合到达角数据通过公式
计算出拟合到达角θn，其中，λ为信号波长，相位差Δφn＝φn-φn-1，d为相邻天线的间距；
根据所述拟合相位数据和所述拟合到达角数据构建输入矩阵，输入矩阵表示为：
/>
其中，
Θ＝[θ1,θ2,…,θn]T，/>
以及，
构建相互关联的模型分量
θn＝θn-1+ηnπn、/>
通过
迭代计算最小二乘解直至/>
时收敛,其中，Ψf为最小二乘解，/>
为经验阈值，f＝1,2,...N，下标f表示前f个数据，即从1-f的数据；
根据模型分量计算出收敛时的信号数据角度估计值θ。
上述技术方案中的一个技术方案具有如下有益效果：
本发明实施例的方法中提出了一种提高室内定位精度的方法及装置，在收到由接收端天线阵列连续发送的信号数据后，从中获取信号样本进行构建以生成原始信号向量；然后，使用拟合函数对原始信号向量中的I通道信号数据和Q通道信号数据进行拟合后，输出拟合相位数据；随后，基于拟合相位数据计算生成输入矩阵，并将输入矩阵代入迭代模型收敛后生成信号数据角度估计值；最后，基于输入矩阵进行滤波处理后生成信号数据相位估计值。本申请通过结合非线性递归最小二乘法和卡尔曼滤波器来提高到达角的估计精度。非线性递归最小二乘法可以消除信号数据处理过程中的多径和路径损耗效应，卡尔曼滤波器可以消除接收信号中天线切换引起的误差影响。
【附图说明】
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。
图1是本发明实施例所提供的提高室内定位精度的方法的流程示意图；
图2为本发明实施例所提供的提高室内定位精度的装置的功能框图。
【具体实施方式】
为使本发明的目的、技术方案和优点更加清楚，下面将结合具体实施例及相应的附图对本发明的技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参考图1，其为本发明实施例所提供的提高室内定位精度的方法的流程示意图，如图所示，该方法包括以下步骤：
S1、收到由接收端天线阵列连续发送的信号数据后，从中获取信号样本进行构建以生成原始信号向量；
S2、使用拟合函数对原始信号向量中的I通道信号数据和Q通道信号数据进行拟合后，输出拟合相位数据；
S3、基于拟合相位数据计算生成输入矩阵，并将输入矩阵代入迭代模型收敛后生成信号数据角度估计值；
S4、基于输入矩阵进行滤波处理后生成信号数据相位估计值。
具体的，S1具体包括：
接收端天线阵列包含Ne个天线单元，其对应的信号数据的信号向量表示为：
其中，n为样本编号，/>
为样本n中第k个天线对应的信号数据；
将信号向量
表示为I通道信号数据分量和Q通道信号数据分量的形式：
其中，fc为载波频率，T为采样间隔，/>
为同相载波的幅度，/>
为正交相位载波的幅度。
进一步地，S2具体包括：
通过给定变量的拟合函数对I通道信号数据和Q通道信号数据进行拟合，其中，I通道信号数据通过如下函数进行拟合：
Q通道通过如下函数进行拟合：/>
通过公式
计算出I通道信号数据中样本n中第k个天线对应的相位/>
通过公式
计算出样本n的加权相位，从而输出拟合相位数据，其中，/>
为样本n中第k个天线对应的方位角，/>
为标准方位角。
进一步地，所述基于拟合相位数据计算生成输入矩阵，具体包括：
根据拟合到达角数据通过公式
计算出拟合到达角θn，其中，λ为信号波长，相位差Δφn＝φn-φn-1，d为相邻天线的间距；
根据所述拟合相位数据和所述拟合到达角数据构建输入矩阵，输入矩阵表示为：
其中，
Θ＝[θ1,θ2,…,θn]T，/>
进一步地，所述输入矩阵代入迭代模型收敛后生成信号数据估计值，具体包括：
构建相互关联的模型分量
θn＝θn-1+ηnπn、/>
/>
通过
迭代计算最小二乘解直至/>
时收敛,其中，Ψf为最小二乘解，/>
为经验阈值，f＝1,2,…N，下标f表示前f个数据，即从1-f的数据；
根据模型分量计算出收敛时的信号数据角度估计值θ。
进一步地，所述S4具体包括：
建立状态模型：
其中
为当前状态的相位估计值，/>
为上一状态的相位估计值，K(k)为凯尔曼增益，Z(k)为最新的测量值，/>
为根据先验参数计算的预测值；
将状态模型拟合至卡尔曼滤波器后进行滤波处理后得到信号数据相位估计值。
本发明实施例的方法中提出了一种提高室内定位精度的方法及装置，在收到由接收端天线阵列连续发送的信号数据后，从中获取信号样本进行构建以生成原始信号向量；然后，使用拟合函数对原始信号向量中的I通道信号数据和Q通道信号数据进行拟合后，输出拟合相位数据；随后，基于拟合相位数据计算生成输入矩阵，并将输入矩阵代入迭代模型收敛后生成信号数据角度估计值；最后，基于输入矩阵进行滤波处理后生成信号数据相位估计值。本申请通过结合非线性递归最小二乘法和卡尔曼滤波器来提高到达角的估计精度。非线性递归最小二乘法可以消除信号数据处理过程中的多径和路径损耗效应，卡尔曼滤波器可以消除接收信号中天线切换引起的误差影响。
本发明实施例进一步给出实现上述方法实施例中各步骤及方法的装置实施例。
请参考图2，其为本发明实施例所提供的提高室内定位精度的装置的功能框图，所述装置包括：
预处理模块210，用于在收到由接收端天线阵列连续发送的信号数据后，从中获取信号样本进行构建以生成原始信号向量；
拟合模块220，用于使用拟合函数对原始信号向量中的I通道信号数据和Q通道信号数据进行拟合后，输出拟合相位数据；
计算模块230，用于基于拟合相位数据计算生成输入矩阵，并将输入矩阵代入迭代模型收敛后生成信号数据角度估计值；
滤波模块240，用于基于输入矩阵进行滤波处理后生成信号数据相位估计值。
具体地，预处理模块具体用于：
接收端天线阵列包含Ne个天线单元，其对应的信号数据的信号向量表示为：
其中，n为样本编号，/>
为样本n中第k个天线对应的信号数据；
将信号向量
表示为I通道信号数据分量和Q通道信号数据分量的形式：
其中，fc为载波频率，T为采样间隔，/>
为同相载波的幅度，/>
为正交相位载波的幅度。/>
具体地，拟合模块具体用于：
通过给定变量的拟合函数对I通道信号数据和Q通道信号数据进行拟合，其中，I通道信号数据通过如下函数进行拟合：
Q通道通过如下函数进行拟合：/>
通过公式
计算出I通道信号数据中样本n中第k个天线对应的相位/>
通过公式
计算出样本n的加权相位，从而输出拟合相位数据，其中，/>
为样本n中第k个天线对应的方位角，/>
为标准方位角。
具体地，计算模块具体用于：
根据拟合到达角数据通过公式
计算出拟合到达角θn，其中，λ为信号波长，相位差Δφn＝φn-φn-1，d为相邻天线的间距；
根据所述拟合相位数据和所述拟合到达角数据构建输入矩阵，输入矩阵表示为：
其中，
Θ＝[θ1,θ2,…,θn]T，/>
以及，
构建相互关联的模型分量
θn＝θn-1+ηnπn、/>
通过
迭代计算最小二乘解直至/>
时收敛,其中，Ψf为最小二乘解，/>
为经验阈值，f＝1,2,…N，下标f表示前f个数据，即从1-f的数据；
根据模型分量计算出收敛时的信号数据角度估计值θ。
由于本实施例中的各单元模块能够执行图1所示的方法，本实施例未详细描述的部分，可参考对图1的相关说明。
在硬件层面，该装置可以包括处理器，可选地还包括内部总线、网络接口、存储器。其中，存储器可能包含内存，例如高速随机存取存储器(Random-Access Memory，RAM)，也可能还包括非易失性存储器(non-volatile memory)，例如至少1个磁盘存储器等。当然，该装置还可能包括其他业务所需要的硬件。
处理器、网络接口和存储器可以通过内部总线相互连接，该内部总线可以是ISA(Industry Standard Architecture，工业标准体系结构)总线、PCI(PeripheralComponent Interconnect，外设部件互连标准)总线或EISA(Extended Industry StandardArchitecture，扩展工业标准结构)总线等。所述总线可以分为地址总线、数据总线、控制总线等。
存储器，用于存放程序。具体地，程序可以包括程序代码，所述程序代码包括计算机操作指令。存储器可以包括内存和非易失性存储器，并向处理器提供指令和数据。
结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。
上述实施例阐明的***、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
为了描述的方便，描述以上装置时以功能分为各种单元或模块分别描述。当然，在实施本发明时可以把各单元或模块的功能在同一个或多个软件和/或硬件中实现。
本领域内的技术人员应明白，本发明的实施例可提供为方法、***、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本发明的实施例可提供为方法、***或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本发明，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本发明中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于***实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本发明的实施例而已，并不用于限制本发明。对于本领域技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。
Claims (3)
1.一种提高室内定位精度的方法，其特征在于，所述方法包括：
S1、收到由接收端天线阵列连续发送的信号数据后，从中获取信号样本进行构建以生成原始信号向量；
S2、使用拟合函数对原始信号向量中的I通道信号数据和Q通道信号数据进行拟合后，输出拟合相位数据；
S3、基于拟合相位数据计算生成输入矩阵，并将输入矩阵代入迭代模型收敛后生成信号数据角度估计值；
S4、基于输入矩阵进行滤波处理后生成信号数据相位估计值；
所述S1具体包括：
接收端天线阵列包含Ne个天线单元，其对应的信号数据的信号向量表示为：
n＝1,2,...N，其中，n为样本编号，
为样本n中第k个天线对应的信号数据；
将信号向量
表示为I通道信号数据分量和Q通道信号数据分量的形式：
其中，fc为载波频率，T为采样间隔，
为同相载波的幅度，
为正交相位载波的幅度；
所述S2具体包括：
通过给定变量的拟合函数对I通道信号数据和Q通道信号数据进行拟合，其中，I通道信号数据通过如下函数进行拟合：
Q通道通过如下函数进行拟合：
通过公式
计算出I通道信号数据中样本n中第k个天线对应的相位
通过公式
计算出样本n的加权相位，从而输出拟合相位数据，其中，
为样本n中第k个天线对应的方位角，
为标准方位角；
所述基于拟合相位数据计算生成输入矩阵，具体包括：
根据拟合到达角数据通过公式
计算出拟合到达角θn，其中，λ为信号波长，相位差
d为相邻天线的间距；
根据所述拟合相位数据和所述拟合到达角数据构建输入矩阵，输入矩阵表示为：
其中，
所述输入矩阵代入迭代模型收敛后生成信号数据估计值，具体包括：
构建相互关联的模型分量
θn＝θn-1+ηnπn、
通过
迭代计算最小二乘解直至
时收敛,其中，Ψf为最小二乘解，
为经验阈值，f＝1,2,...N，下标f表示前f个数据，即从1-f的数据；
根据模型分量计算出收敛时的信号数据角度估计值θ。
2.根据权利要求1所述的提高室内定位精度的方法，其特征在于，所述S4具体包括：
建立状态模型：
其中
为当前状态的相位估计值，
为上一状态的相位估计值，K(k)为凯尔曼增益，Z(k)为最新的测量值，
为根据先验参数计算的预测值；
将状态模型拟合至卡尔曼滤波器后进行滤波处理后得到信号数据相位估计值。
3.一种提高室内定位精度的装置，其特征在于，所述装置包括：
预处理模块，用于在收到由接收端天线阵列连续发送的信号数据后，从中获取信号样本进行构建以生成原始信号向量；
拟合模块，用于使用拟合函数对原始信号向量中的I通道信号数据和Q通道信号数据进行拟合后，输出拟合相位数据；
计算模块，用于基于拟合相位数据计算生成输入矩阵，并将输入矩阵代入迭代模型收敛后生成信号数据角度估计值；
滤波模块，用于基于输入矩阵进行滤波处理后生成信号数据相位估计值；
所述预处理模块具体用于：
接收端天线阵列包含Ne个天线单元，其对应的信号数据的信号向量表示为：
n＝1,2,...N，其中，n为样本编号，
为样本n中第k个天线对应的信号数据；
将信号向量
表示为I通道信号数据分量和Q通道信号数据分量的形式：
其中，fc为载波频率，T为采样间隔，
为同相载波的幅度，
为正交相位载波的幅度；
所述拟合模块具体用于：
通过给定变量的拟合函数对I通道信号数据和Q通道信号数据进行拟合，其中，I通道信号数据通过如下函数进行拟合：
Q通道通过如下函数进行拟合：
通过公式
计算出I通道信号数据中样本n中第k个天线对应的相位
通过公式
计算出样本n的加权相位，从而输出拟合相位数据，其中，
为样本n中第k个天线对应的方位角，
为标准方位角；
所述计算模块具体用于：
根据拟合到达角数据通过公式
计算出拟合到达角θn，其中，λ为信号波长，相位差
d为相邻天线的间距；
根据所述拟合相位数据和所述拟合到达角数据构建输入矩阵，输入矩阵表示为：
其中，
以及，
构建相互关联的模型分量
θn＝θn-1+ηnπn、
通过
迭代计算最小二乘解直至
时收敛,其中，Ψf为最小二乘解，
为经验阈值，f＝1,2,...N，下标f表示前f个数据，即从1-f的数据；
根据模型分量计算出收敛时的信号数据角度估计值θ。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210410379.8A CN114755629B (zh) | 2022-04-19 | 2022-04-19 | 提高室内定位精度的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210410379.8A CN114755629B (zh) | 2022-04-19 | 2022-04-19 | 提高室内定位精度的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114755629A CN114755629A (zh) | 2022-07-15 |
| CN114755629B true CN114755629B (zh) | 2023-03-24 |
Family
ID=82330580
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210410379.8A Active CN114755629B (zh) | 2022-04-19 | 2022-04-19 | 提高室内定位精度的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114755629B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115460693B (zh) * | 2022-09-19 | 2023-06-06 | 北京信息职业技术学院 | 基于蓝牙的定位方法及装置 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5251090B2 (ja) * | 2007-11-30 | 2013-07-31 | セイコーエプソン株式会社 | 測位方法、プログラム及び測位装置 |
| US8797213B2 (en) * | 2009-09-30 | 2014-08-05 | Broadcom Corporation | Methods and systems for estimating angle of arrival |
| CN104062628A (zh) * | 2014-07-08 | 2014-09-24 | 北京万集科技股份有限公司 | 一种车载单元obu入射方向角的计算方法及装置 |
| CN107733817B (zh) * | 2016-08-11 | 2021-10-15 | 中兴通讯股份有限公司 | 一种到达角估计的方法、装置、终端及基站 |
| CN106793060B (zh) * | 2017-03-08 | 2020-09-11 | 哈尔滨工程大学 | 一种超宽带室内定位方法 |
| CN110333478B (zh) * | 2018-03-30 | 2022-05-17 | 华为技术有限公司 | 一种到达角度、出发角度确定方法及通信装置 |
| US11079465B2 (en) * | 2018-10-30 | 2021-08-03 | Electronics And Telecommunications Research Institute | Method and apparatus for estimating location of signal source |
| CN112073895A (zh) * | 2019-06-10 | 2020-12-11 | 重庆邮电大学 | 一种基于csi高精度人员定位跟踪方法 |
| CN113050030B (zh) * | 2021-03-26 | 2021-11-05 | 江苏蓝策电子科技有限公司 | 基于到达角度测距的定位方法及装置 |
| CN114217266B (zh) * | 2021-11-30 | 2024-06-11 | 清华大学 | 基于天线数据时域拟合的降低信道误差的处理方法 |
-
2022
- 2022-04-19 CN CN202210410379.8A patent/CN114755629B/zh active Active
Non-Patent Citations (3)
| Title |
|---|
| Multi-antenna Array-based AoA Estimation Using;shuai he;《2021 the 7th International Conference on Computer and Communications》;20211231;全文 * |
| 基于最小二乘法的改进卡尔曼滤波算法;符强等;《桂林电子科技大学学报》;20191031(第05期);全文 * |
| 基于空间调制***的低复杂度天线选择算法;贾东升等;《广东通信技术》;20180131;第38卷(第2期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114755629A (zh) | 2022-07-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113050030B (zh) | 基于到达角度测距的定位方法及装置 | |
| Rife et al. | Single tone parameter estimation from discrete-time observations | |
| CN104185974B (zh) | 用于宽带rf接收器的实时i/q不平衡校正的装置和方法 | |
| CN114755629B (zh) | 提高室内定位精度的方法及装置 | |
| CN108989256B (zh) | 一种fsk/gfsk解调方法及装置 | |
| CN109729032B (zh) | 一种频偏估计值的校正方法、装置及计算机可读存储介质 | |
| CN111914506B (zh) | 天线阵测向方法、装置及电子设备 | |
| CN116232558B (zh) | 一种iq时延补偿方法、装置、电子设备和存储介质 | |
| CN111487582A (zh) | 获得蓝牙阵列天线参数标定模型、到达角的方法和装置 | |
| CN116015492A (zh) | 一种无线电台故障诊断方法、***及设备 | |
| CN113935402B (zh) | 时差定位模型的训练方法、装置及电子设备 | |
| CN113630212B (zh) | 一种数字信号的降采样方法及降采样装置 | |
| CN110429921A (zh) | 一种变步长lms自适应滤波方法及其存储介质 | |
| Wang et al. | Adaptive time delay estimation algorithm for indoor near‐field electromagnetic ranging | |
| CN110113287B (zh) | 一种信号确定方法、装置、电子设备及存储介质 | |
| CN114786127B (zh) | 一种蓝牙aoa的定位方法及装置 | |
| CN112235002B (zh) | 谐波消除方法及装置、存储介质、终端 | |
| CN115022963B (zh) | 基于蓝牙的室内定位方法及装置 | |
| US20230007488A1 (en) | Signal spoof detection at base stations | |
| CN115390113A (zh) | 一种面向电力***的bd3高精度定位方法及*** | |
| CN111723686A (zh) | 一种信号个数估计方法、装置、电子设备和存储介质 | |
| CN113411148B (zh) | 移频方法、装置、电子设备及存储介质 | |
| CN113132125B (zh) | 故障原因定位方法、装置、电子设备及存储介质 | |
| CN108510996A (zh) | 一种快速迭代自适应滤波方法 | |
| CN114514525B (zh) | 载波频率、初始相位、相位噪声的估计方法和相关设备 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |