CN105866738B - 一种用于提高室内可见光定位***定位精度的迭代算法 - Google Patents
一种用于提高室内可见光定位***定位精度的迭代算法 Download PDFInfo
- Publication number
- CN105866738B CN105866738B CN201610245106.7A CN201610245106A CN105866738B CN 105866738 B CN105866738 B CN 105866738B CN 201610245106 A CN201610245106 A CN 201610245106A CN 105866738 B CN105866738 B CN 105866738B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- mtr
- 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
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/16—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using electromagnetic waves other than radio waves
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/11—Arrangements specific to free-space transmission, i.e. transmission through air or vacuum
- H04B10/114—Indoor or close-range type systems
- H04B10/116—Visible light communication
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Navigation (AREA)
Abstract
一种用于提高室内可见光定位***定位精度的迭代算法,包括基于接收信号强度的粗定位以及结合接收信号强度和估计反射信号强度的迭代定位两部分:首先,根据接收光信号的强度进行粗定位;随后根据粗定位的结果和预先计算的在接收平面网格节点上的反射信号强度对粗定位位置处的反射信号强度进行估算;再结合接收信号强度和估算的反射信号强度,通过对定位结果和定位方程进行多次迭代,就可以得到最终的精确定位结果。由于该算法能够对来自墙面及其他物体表面的漫反射光信号进行估计并抑制了这些反射信号的影响,因此,应用于实际场景中能实现较为准确的室内可见光定位。
Description
技术领域
本发明涉及照明集成的室内可见光定位技术,尤其涉及一种用于提高室内可见光定位***定位精度的迭代算法,属于可见光通信技术领域。
背景技术
随着定位导航技术的广泛应用和基于位置服务的不断发展,室内定位技术逐渐进入人们的视野,相关需求也不断增长。在定位与导航市场中,GPS和北斗等基于卫星信号的定位***占据了绝大多数。但是与此同时,由于在城市中卫星信号会被高楼遮挡;而且在室内,卫星信号经过穿墙之后会被严重衰减,因此GPS和北斗这类基于卫星信号的定位技术在室内环境下并不实用。近年来,随着绿色照明的发展,LED以其高发光效率、低能耗、低成本、高显色性和环境友好等优势逐渐得到普及。而基于LED光源的可见光通信和可见光定位技术也不断发展,出现了多种原理的与可见光定位相关的论文和专利。
与基于Wi-Fi和蓝牙的定位技术相比,可见光定位技术不会产生电磁辐射,且不受电磁信号的干扰。因此其定位精度高于基于Wi-Fi和蓝牙定位,且更加稳定。而与基于超声、红外和超宽带(UWB)的室内定位技术相比,可见光定位技术可直接利用室内的照明LED灯,无需额外假设信号发射装置。因此其成本更低,更易于推广和实用。
对于当前的室内可见光定位技术,出于成本的考虑,很多都是在接收端采用光电探测器,同时利用接收到的来自各个光源的光信号强度进行定位,例如专利“一种基于频移键控的室内可见光定位及信息推送方法”(申请号201510379394.0)和专利“基于可见光信号强度的三维空间室内定位方法”(申请号201510063426.6)。但是对于这一类的定位技术而言,都只考虑到了直接来自光源的直射信号,而忽略了来自于墙面或家具的反射光信号,因此在实际应用中存在较大误差。
发明内容
本发明的目的在于解决目前室内可见光定位技术受到反射光信号干扰的问题。提出一种用于提高室内可见光定位***定位精度的迭代算法。
本发明的技术核心为:其利用基于接收信号强度的粗定位以及结合接收信号强度和估计反射信号强度的迭代定位来进行基于可见光信号的室内定位,并通过不断迭代来提高定位精度;本发明一种用于提高室内可见光定位***定位精度的迭代算法,具体步骤如下:
步骤一、预处理,具体包括:
1.1在室内天花板上放置N个光源,每个光源都采用强度调制的方式进行驱动,同时每个光源通过光信号发送自己的空间坐标,在完成照明功能的同时作为可见光定位***的信号源;
其中,所述的N>3;
1.2根据室内环境,将可见光定位***的接收端所处的接收平面分为R个网格节点,并求出或测量出在各个网格节点处分别来自各个光源的经墙面和家具所反射光信号的强度;
其中,所述的可见光定位***的接收端位于光源下方的某个平面上,假设接收端与光源的垂直距离为h,且竖直向上放置;
所述的R的选择依据主要为***定位精度和定位时间要求,且优选的R值密度范围为:100个/平方米到400个/平方米之间;
1.3初始化迭代次数K为1,预设迭代次数M;
步骤二、可见光定位***的接收端根据接收到的来自各个光源的光信号强度进行粗定位,具体过程为:
假设第i个光源所发射的原始光信号强度为Pti,接收端接收到的来自第i个光源的光信号强度为Pri,第i个光源的坐标为(xi,yi);通过求解以下超定方程组的最小二乘解,可以获得粗定位结果的坐标(xc,yc),并将其作为迭代定位算法的初值;
步骤三、寻找距离步骤二输出的粗定位结果坐标(xc,yc)最近的S个网格节点:结合预处理中所得的各个网格节点处分别来自各个光源的经墙面和家具所反射光信号的强度,可以利用插值的方法,估计在粗定位结果的位置上来自第i个光源的经墙面和家具所反射光信号的强度,记为Prdi;
其中,S的大小由网格节点的分布形式决定,通常可为3、4或6,分别对应
三角形排布矩形排布以及六边形排布;
步骤四、根据步骤三估计的粗定位坐标附近的反射光信号强度,求解以下方程组,重新计算接收端的位置坐标(x,y);
步骤五、判断迭代次数K是否到达预设迭代次数M,并根据判断结果进行如下操作:
5.1当K<M,即没有达到预定迭代次数,将K值加1,并将步骤四重新计算的接收端坐标(x,y)作为新的粗定位坐标(xc,yc),跳至步骤三;
5.2当K=M,即已经达到预定迭代次数,步骤四得出的接收端坐标(x,y)即为最终的定位结果;
其中,所述的迭代次数K由***的定位准确度和复杂度要求决定。
有益效果
本发明同现有技术中“一种基于频移键控的室内可见光定位及信息推送方法”和“基于可见光信号强度的三维空间室内定位方法”相比,具有以下有益效果:
1.考虑并从一定程度上消除了反射光对室内可见光定位***的影响,在同等条件下的定位精度更高,对于实施例1中的情况,可以将室内的平均定位误差减小10倍以上;
2.利用迭代算法在计算中不断提高定位精度,通过设置不同的迭代次数,可在定位精度和算法复杂度之间任意取舍。
附图说明
下面结合附图与具体实施方式对本发明作更进一步的说明:
图1为本发明一种用于提高室内可见光定位***定位精度的迭代算法的算法流程图;
图2为本发明一种用于提高室内可见光定位***定位精度的迭代算法实施例1中步骤S1网格节点示意图;
图3为本发明一种用于提高室内可见光定位***定位精度的迭代算法实施例1中步骤S5反射信号强度估计算法中网格节点示意图。
具体实施方式
下面,结合附图详细描述本发明的优选实施方式。为了便于描述和突出显示本发明,附图中省略了现有技术中已有的相关部件,并将省略对这些公知部件的描述。
实施例1
具体结合图1、图2和图3对本发明一种用于提高室内可见光定位***定位精度的迭代算法进行细致阐述。图1为本发明一种用于提高室内可见光定位***定位精度的迭代算法及实施例中的算法流程图;图2为本发明一种用于提高室内可见光定位***定位精度的迭代算法及实施例1中步骤S1网格节点示意图;图3为本发明一种用于提高室内可见光定位***定位精度的迭代算法实施例1中步骤S5反射信号强度估计算法中网格节点示意图。
在本实施例中,假设在室内房间中进行定位,房间的大小为5米×5米×3米,天花板上有4个LED灯在提供照明的同时作为可见光定位***的信号源。这4个LED灯采用不同频率的开关键控(OOK)信号调制,并循环发送自己的位置坐标,分别为(1,1,3),(1,3,3),(3,1,3)和(3,3,3)。与此同时,接收端位于光源下方且与光源垂直距离为3米的平面上,并竖直向上放置。图1中所描述的是本发明中定位算法的总体流程图,其主要步骤如下:
步骤S1:根据室内环境,将接收端所处的接收平面均匀分为2601个网格节点,每个节点与相邻节点的间距为10厘米(如图2所示),此时节点密度为100个/平方米;随后利用可见光信道传输模型求出在各个网格节点处分别来自各个光源的经墙面和家具所反射光信号的强度;初始化迭代次数K为1,预设迭代次数M为3;
步骤S2:在接收端接收来自4个不同LED光源的光信号。根据不同的频率对其进行识别并分别测量光信号强度,分别记为Pr1,Pr2,Pr3和Pr4;与此同时,对分别来自4个LED光源的光信号进行解调,并获取他们的位置坐标,即(1,1,3),(1,3,3),(3,1,3)和(3,3,3),并将光信号强度和光源位置坐标一一对应;
步骤S3:在定位时,首先根据接收到的来自4个光源的光信号强度进行粗定位,假设第i个光源所发射的原始光信号强度为Pti,列出超定方程组为:
步骤S4:求解步骤S3中所给出的超定方程组的最小二乘解,可以获得粗定位结果坐标,并将其作为迭代定位算法的初值;
步骤S5:在获取了粗定位结果的坐标后,寻找距离其最近的4个网格节点(如图3所示),假设对于第m个节点,来自光源i的经墙面和家具所反射光信号的强度为Prdim;利用双线性插值的方法,利用下式可以估计在粗定位结果的位置上,来自光源i的经墙面和家具所反射光信号的强度为Prdi:
其中,Aci为以粗定位位置与节点i之间连线作为对角线的矩形的面积;
步骤S6:根据估计的粗定位位置的反射光信号强度,求解以下方程组,重新计算接收端的位置坐标;
步骤S7:统计迭代次数K;
步骤S8:对于定位精度(定位速度)优先的***,判定当前已经进行过的迭代次数K是否等于预设的M值,并进行相应操作:
S8.1若等于M=3次,将定位结果(即接收端的位置)输出,定位结束;
S8.2若少于M=3次,则将步骤S6中求得的接收端位置坐标代替粗定位结果,作为新的迭代初值,将迭代次数K加1,并跳转至步骤S5,继续进行迭代定位的计算。
实施例2
将实施例1中的4个LED灯改成N(N>4)个灯,该定位***依然有效。
实施例3
将实施例1中的LED灯改成其他固体光源如OLED,该定位***依然有效。为了更好地理解本发明,对有关内容补充说明如下:
关于反射光信号的预处理:在本算法中,需要预先获取接收平面网格节点上来自各个光源的反射信号光强。这些光强可由多种仿真算法或软件对室内环境仿真计算得出,也可通过在实际场景中的实地测量获得。上述实施方式中未述及的有关技术内容采取或借鉴已有技术即可实现,例如实施例1中发射端的的可见光信号调制驱动电路和接收端的数据采集,均可采用专利“一种基于频移键控的室内可见光定位及信息推送方法”(申请号201510379394.0)中所述的方法进行。
需要说明的是,本说明书所述的只是本发明的3种较佳具体实施例,以上实施例仅用于说明本发明的技术方案而非对本发明的限制。凡本领域技术人员依本发明的构思通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在本发明的范围之内。
Claims (5)
1.一种用于提高室内可见光定位***定位精度的迭代算法,其特征在于:具体步骤如下:
步骤一、预处理,具体包括:
1.1在室内天花板上放置N个光源,每个光源都采用强度调制的方式进行驱动,同时每个光源通过光信号发送自己的空间坐标,在完成照明功能的同时作为可见光定位***的信号源;
1.2根据室内环境,将可见光定位***的接收端所处的接收平面分为R个网格节点,并求出或测量出在各个网格节点处分别来自各个光源的经墙面和家具所反射光信号的强度;
其中,所述的可见光定位***的接收端位于光源下方的某个平面上,假设接收端与光源的垂直距离为h,且竖直向上放置;
1.3初始化迭代次数K为1;
步骤二、可见光定位***的接收端根据接收到的来自各个光源的光信号强度进行粗定位,步骤二具体过程为:
假设第i个光源所发射的原始光信号强度为Pti,接收端接收到的来自第i个光源的光信号强度为Pri,第i个光源的坐标为(xi,yi);通过求解以下超定方程组的最小二乘解,可以获得粗定位结果的坐标(xc,yc),并将其作为迭代定位算法的初值;
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mn>3</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mn>3</mn>
</mrow>
</msub>
</mrow>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>N</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
步骤三、寻找距离步骤二输出的粗定位结果坐标(xc,yc)最近的S个网格节点;结合预处理中所得的各个网格节点处分别来自各个光源的经墙面和家具所反射光信号的强度,具体为:利用插值的方法,估计在粗定位结果的位置上来自第i个光源的经墙面和家具所反射光信号的强度,记为Prdi;
其中,S的大小由网格节点的分布形式决定,通常可为3、4或6,分别对应三角形排布矩形排布以及六边形排布;
步骤四、根据步骤三估计的反射光信号的强度,重新计算接收端的位置坐标;
步骤五、判断迭代次数K是否到达预设迭代次数M,并根据判断结果进行如下操作:
5.1当K<M,即没有达到预定迭代次数,将K值加1,并将步骤四重新计算的接收端坐标(x,y)作为新的粗定位坐标(xc,yc),跳至步骤三;
5.2当K=M,即已经达到预定迭代次数,步骤四得出的接收端坐标(x,y)即为最终的定位结果;
至此,步骤一至步骤五完成了一种用于提高室内可见光定位***定位精度的迭代算法。
2.如权利要求1所述的一种用于提高室内可见光定位***定位精度的迭代算法,其特征还在于:
步骤一中所述的N>3。
3.如权利要求1所述的一种用于提高室内可见光定位***定位精度的迭代算法,其特征还在于:
步骤一中所述的R的选择依据主要为***定位精度和定位时间要求,且优选的R值密度范围为:100个/平方米到400个/平方米之间。
4.如权利要求1所述的一种用于提高室内可见光定位***定位精度的迭代算法,其特征还在于:
步骤四中重新计算接收端的位置坐标(x,y)采用解如下方程组得出;
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>d</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mn>3</mn>
</mrow>
</msub>
</mrow>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<msqrt>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>N</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>d</mi>
<mi>N</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>d</mi>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
5.如权利要求1所述的一种用于提高室内可见光定位***定位精度的迭代算法,其特征还在于:
步骤五中所述的迭代次数K由***的定位准确度和复杂度要求决定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610245106.7A CN105866738B (zh) | 2016-04-19 | 2016-04-19 | 一种用于提高室内可见光定位***定位精度的迭代算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610245106.7A CN105866738B (zh) | 2016-04-19 | 2016-04-19 | 一种用于提高室内可见光定位***定位精度的迭代算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105866738A CN105866738A (zh) | 2016-08-17 |
CN105866738B true CN105866738B (zh) | 2017-08-25 |
Family
ID=56633181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610245106.7A Active CN105866738B (zh) | 2016-04-19 | 2016-04-19 | 一种用于提高室内可见光定位***定位精度的迭代算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105866738B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106441304B (zh) * | 2016-10-08 | 2019-05-21 | 东南大学 | 一种利用可见光通信的三维室内定位方法 |
CN106595662B (zh) * | 2016-11-25 | 2019-05-14 | 大连理工大学 | 一种可见光定位中基于光源阵列旋转特性的抗干扰方法 |
CN108732537B (zh) * | 2018-05-08 | 2022-04-05 | 北京理工大学 | 一种基于神经网络和接收信号强度的室内可见光定位方法 |
CN112394322B (zh) * | 2020-05-27 | 2023-06-20 | 南京邮电大学 | 一种基于两步法实现室内可见光通信***精确定位的方法 |
CN112462330B (zh) * | 2020-07-10 | 2021-10-01 | 杭州石峦科技有限公司 | 特定区域内物距定位***与方法 |
CN113405551B (zh) * | 2021-06-09 | 2022-10-04 | 北京科技大学 | 一种基于区域划分的自适应可见光室内定位方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101241793B1 (ko) * | 2012-11-01 | 2013-03-15 | (주)위니텍 | 실내 측위용 비콘 모듈을 구비한 형광등 소켓 커넥터 및 이를 이용한 실내 측위 시스템 |
CN103199923B (zh) * | 2013-04-22 | 2016-04-20 | 中国矿业大学 | 一种基于可见光通信的井下移动目标光指纹定位跟踪方法 |
CN103389488B (zh) * | 2013-07-30 | 2015-07-29 | 中国科学院深圳先进技术研究院 | 一种基于光强度的多光源室内定位装置与方法 |
CN104991228B (zh) * | 2015-02-06 | 2016-11-02 | 北京理工大学 | 基于可见光信号强度的三维空间室内定位方法 |
CN105068047A (zh) * | 2015-07-01 | 2015-11-18 | 北京理工大学 | 一种基于频移键控的室内可见光定位及信息推送方法 |
-
2016
- 2016-04-19 CN CN201610245106.7A patent/CN105866738B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105866738A (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105866738B (zh) | 一种用于提高室内可见光定位***定位精度的迭代算法 | |
CN106100734B (zh) | 一种基于人工神经网络的高精度室内可见光定位方法 | |
Lv et al. | High accuracy VLC indoor positioning system with differential detection | |
CN104965510B (zh) | 一种基于可见光定位与导航的物流车控制*** | |
Kouhini et al. | LiFi positioning for industry 4.0 | |
CN106646368B (zh) | 一种基于指纹匹配的可见光通信场景中三维定位方法 | |
CN107063266A (zh) | 基于室内led灯和智能手机摄像头实现室内定位的方法 | |
CN102360076B (zh) | 一种多维定位方法和装置 | |
CN107796400B (zh) | 基于可见光通信的室内三维定位方法 | |
CN105738868A (zh) | 一种基于双图像传感器和可见光通信技术的高精度室内定位方法 | |
US20140085642A1 (en) | Lighting apparatus for measuring the position of a mobile terminal, and position measuring system using same | |
CN105629080A (zh) | 一种电磁分布态势图或路径电磁分布图的绘制方法 | |
CN104360313B (zh) | 一种基于多led光源定位的方法 | |
CN110018508A (zh) | 一种定位方法及装置 | |
CN104391273A (zh) | 一种基于圆形投影的可见光定位方法及*** | |
CN107786266A (zh) | 基于蚁群算法的可见光通信室内三维定位的方法及*** | |
CN115032585B (zh) | 一种基于rss测距和射线追踪技术的非视距场景定位方法 | |
CN203416427U (zh) | 基于ZigBee技术的自组网定位*** | |
KR101031129B1 (ko) | 휴대용 주차위치 표시기 | |
CN105652236A (zh) | 一种基于ZigBee技术的商场室内无线定位方法和*** | |
CN206892313U (zh) | 一种抗nlos干扰的可见光定位硬件*** | |
CN107919910A (zh) | 一种抗nlos干扰的可见光定位方法 | |
KR20190001386A (ko) | 저전력 블루투스 기반의 실내 측위 방법 및 장치 | |
CN106772505A (zh) | 基于伪卫星信号载噪比的导航定位方法 | |
Lakafosis et al. | From single-to multihop: The status of wireless localization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |