CN107883955A - 结合椭圆aco‑ofdm的led定位方法及系统 - Google Patents
结合椭圆aco‑ofdm的led定位方法及系统 Download PDFInfo
- Publication number
- CN107883955A CN107883955A CN201710979861.2A CN201710979861A CN107883955A CN 107883955 A CN107883955 A CN 107883955A CN 201710979861 A CN201710979861 A CN 201710979861A CN 107883955 A CN107883955 A CN 107883955A
- Authority
- CN
- China
- Prior art keywords
- mrow
- led
- msub
- aco
- signal
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
- G01C21/206—Instruments for performing navigational calculations specially adapted for indoor navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optical Communication System (AREA)
Abstract
本发明公开了一种结合椭圆ACO‑OFDM的LED定位方法及系统，用于减轻室内环境中多径失真的影响，LED用作具有唯一ID的发射机以区分其位置，并且移动接收机基于RSS定位算法以获得自身在水平位置上的定位，并进一步采用信赖域技术实现三维定位，定位精度高。
Description
技术领域
本发明涉及位置定位领域，尤其涉及LED定位方面，更具体地说，涉及一种结合椭圆ACO-OFDM的LED定位方法及系统。
背景技术
位置服务通过移动终端和无线网络的配合，确定出移动用户的实际地理位置，从而提供用户需要的与位置相关的信息服务。全球定位系统(Global SatellitePositioning,GPS)是目前应用最为广泛的定位系统，但卫星信号不能穿透建筑物，且室内环境存在严重的多径和非视距干扰，在室内难以实现定位.为了提高室内定位的精度，人们采用了红外线定位、超声波定位、射频识别定位等技术，但传统的室内定位技术在系统成本、稳定性、定位精度等方面不能完全满足用户的需求。随着发光二极管(light-emitting-diode,LED)照明技术的发展，可见光通信技术(Visible Light Communication,VLC)也得到发展，基于可见光通信的室内定位导航技术成为了研究的热点。与传统的射频定位和其它无线定位技术相比，基于LED可见光通信的室内定位具有发射功率高、定位精度高、使用场合广、无电磁干扰和节约能源等优点。
然而在实际应用中，LED定位系统都会受到室内环境中多径干扰的影响，期刊《Journal of Lightwave Technology》于2016年第34卷第10期刊登的文章《Impact ofMultipath Reflections on the Performance of Indoor Visible Light PositioningSystems》研究了多径干扰对定位精度的影响，结果表明多径干扰会大大降低定位精度，所以多径干扰需考虑到定位算法中.正交频分复用技术(orthogonal frequency divisionmultiplexing,OFDM)一直被普遍应用于抵抗多径干扰引起的码间串扰问题，并提高传输速率。与传统射频无线通信不同的是，VLC使用强度调制直接检测(intensity-modulationdirect-detection,IM/DD)，要求以光强为载体的发送信号必须是非负实数。直流偏置光正交频分复用(direct-current-biased optical OFDM,DCO-OFDM)通过叠加直流偏置获得非负传输信号，较小的直流偏置会引起较大的限幅噪声，而较大的直流偏置则会降低功率效率。非对称限幅光正交频分复用(asymmetrically clipped optical OFDM,ACO-OFDM)仅调制奇载波，脉冲幅度调制-离散多音频调制(pulse amplitude modulated discretemultitone,PAM-DMT)仅调制载波虚部，两者均通过零值限幅获得非负传输信号，限幅噪声只影响未调制载波部分，但零值限幅使发送端一半能量耗费在未调制载波部分，致使ACO-OFDM和PAM-DMT相较于传统双极性OFDM有3dB的性能损失。期刊《Journal of LightwaveTechnology》于2013年第31卷第7期刊登的文章《Armstrong J.Comparison of ACO-OFDM,DCO-OFDM and ADO-OFDM in IM/DD Systems》提出了一种将ACO-OFDM和DCO-OFDM混合传输的方法，提高了传输速率。为了进一步提高频谱利用率并降低传输复杂性，期刊《Journalof Optical Communications and Networking》于2014年第6卷第4期刊登的文章《HybridAsymmetrically Clipped OFDM-Based IM/DD Optical Wireless System》提出了将ACO-OFDM和PAM-DMT混合传输的方法。
目前大多数LED定位系统都使用LED作为发射机，固定在目标上的光探测器作为接收机，并利用收发机之间的距离来计算接收机的坐标，该距离的获取方法有接收时间(time-of-arrival,TOA)(可参考期刊《Journal of Lightwave Technology》于2013年第31卷第20期刊登的文章《Accuracy of Time-of-Arrival Based Ranging Using VisibleLight With Application in Indoor Localization Systems》)、接收时间差(time-difference-arrival,TDOA)(可参考期刊《IEEE Transactions on ConsumerElectronics》于2011年第57卷第4期刊登的文章《TDOA-based optical wireless indoorlocalization using LED ceiling lamps》)和接收信号的强度(received-signal-strength,RSS)(可参考期刊《IEEE Transactions on Systems,Man,and Cybernetics,Part C(Applications and Reviews)》于2007年第37卷第6期刊登的文章《Survey ofWireless Indoor Positioning Techniques and Systems》以及期刊《Optics Express》于2015年第23卷第16期刊登的文章《Exploring the effect of diffuse reflection onindoor localization systems based on RSSI-VLC》)三种。RSS通过发射信号的能量衰减来计算距离，由于其原理简单而被广泛使用。期刊《Optical Engineering》于2016年第55卷第5期刊登的文章《Indoor location estimation with optical-based orthogonalfrequency division multiplexing communications》提出了将ACO-OFDM与定位系统结合的算法，但其定位精度还有待进一步提高，且其具有ACO-OFDM系统普遍存在的发送端一半能量耗费在未调制载波部分的问题，会产生3dB性能损失。
发明内容
本发明要解决的技术问题在于，针对上述的技术缺陷，提供了一种结合椭圆ACO-OFDM的LED定位方法及系统。
根据本发明的其中一方面，本发明为解决其技术问题，提供了一种结合椭圆ACO-OFDM的LED定位方法，包含：
发射端处理步骤：
A1、各LED发射机分别将包含自身ID编码及位置信息的输入信号分别进行M元正交调幅调制在奇载波上，生成长度为N/4的复数数据，偶载波被全部置零，其中，M为星座图的大小，N为子载波数，并将第一和第N/2个子载波也置零，得到频域信号()*为表示向量的复数共轭；
A2、各LED发射机分别将各自的频域信号C作为输入进行N点IFFT变换，得到时域信号对应的时域信号c＝{caco[0],caco[1],caco[2],…,caco[N-1]}；
A3、各LED发射机分别将各自的时域信号c＝{caco[0],caco[1],caco[2],…,caco[N-1]}发送出去，其中任意caco[n]在被发送时，只传送其对应的虚部其中n＝0、1、2…、N-1，为caco[n]的相位，α为一个大小为的换算因子以将约束在区间[-π/β,π/β]内，σc为caco[n]的标准偏差，β为一个大于2的正数，a和b分别为caco[n]在当焦点落在实轴时椭圆的长轴与短轴；
接收端处理步骤：
B1、接接收各LED发射机发射的信号，获取各LED发射机的ID编码及位置信息并分别计算出各LED发射机的计算信号；
B2、根据与ID编码对应的计算信号、位置信息以及发射功率的大小，采用RSS方法，计算出接收机的水平坐标。
在本发明的LED定位方法的步骤B2中，接收机的水平坐标通过下述公式计算得出：
Q＝(MTM)-1MTN，
其中，且Q为所求的接收机坐标(x,y)，Xi、Yi为第i个发射机的水平面的横坐标和纵坐标，ri为接收机与第i个发射机之间的水平距离，i＝1、2、…、n，n为发射机的个数，ri为通过下述公式所描述：
式中，Pt为第i个LED的重构后的PAM-DMT信号的功率的取平均值，Pri为第i个LED的PAM-DMT信号的发射功率的平均值，h为接收机与接收到的发射光所对应LED灯的垂直距离，朗伯阶数为m＝-ln2/ln(cosΨ1/2)，Ψ1/2是用于发射信号的LED的半功率角，A是用于接收信号的光电探测器的物理面积，Di为接收机与第i个发射机之间的距离，θi为第i盏灯发射角，为第i个接收机相对于其法线方向的入射角，为第i个LED的光滤波器的传播系数，是第i个LED的光集中器的增益，大小由式下述公式确定：
其中，Φc为光集中器视角的半角度，n为光集中器的折射率。
在本发明的LED定位方法中，接收端处理步骤还包括步骤B3，步骤B3包含：
B31、初始化接收机坐标的初始值h0为接收机与接收到的发射光所对应LED灯的垂直距离h的初始值，H为接收到的发射光所对应LED灯所在房间的高度，其中，接收机的水平面中坐标(x,y)通过下述步骤B2得出：
B32、根据计算得到初始定位误差根据初始信赖域半径R0，取以为圆心、以R0为半径的信赖域内部的点更新通过更新
其中定位误差计算公式为：
(Xi,Yi,Zi)为接收到的发射光所对应的第i个LED灯的三维坐标，为计算时的接收机的三维坐标；
B33、进行定位误差判决，当定位误差大于设定的定位误差时，重新选取信赖域内部的点更新当定位误差小于等于设定的定位误差时，输出最新的
B34、进行迭代终止判决，将最新的中与上一次迭代得到的中的进行比较，当二者差值大于预设值时，以最新的中的确定下一次和更新迭代的初始值h0，返回步骤B31开始下一次迭代，否则终止迭代，输出最终的作为接收机的三维位置坐标。
在本发明的LED定位方法的步骤B1中，包含对接收到的每个LED发射机的信号分别进行如下处理的步骤：
对接收到的信号进行数模转换，然后移除循环前缀，对移除循环前缀后的信号进行椭圆ACO-OFDM逆变换，对变换后的结果进行FFT变换，然后估算出信号。
在本发明的LED定位方法中，发射机的信号发射方向和接收机的信号接收方向均于竖直方向。
根据本发明的另一方面，本发明为解决其技术问题，还提供了一种结合椭圆ACO-OFDM的LED定位系统，包含：
发射端处理模块，包含：
输入调制子模块，用于将各LED发射机分别将包含自身ID编码及位置信息的输入信号分别进行M元正交调幅调制在奇载波上，生成长度为N/4的复数数据，偶载波被全部置零，其中，M为星座图的大小，N为子载波数，并将第一和第N/2个子载波也置零，得到频域信号()*为表示向量的复数共轭；
时域变换子模块，用于将各LED发射机分别将各自的将频域信号C作为输入进行N点IFFT变换，得到时域信号对应的时域信号c＝{caco[0],caco[1],caco[2],…,caco[N-1]}；
虚部处理子模块，用于将各LED发射机分别将各自的时域信号c＝{caco[0],caco[1],caco[2],…,caco[N-1]}发送出去，其中任意caco[n]在被发送时，只传送其对应的虚部其中n＝0、1、2…、N-1，为caco[n]的相位，α为一个大小为的换算因子以将约束在区间[-π/β,π/β]内，σc为caco[n]的标准偏差，β为一个大于2的正数，a和b分别为caco[n]在当焦点落在实轴时椭圆的长轴与短轴；
接收端处理模块，包含：
接收估计子模块，用于接收各LED发射机发射的信号，获取各LED发射机的ID编码及位置信息并分别计算出各LED发射机的计算信号；
2D坐标计算子模块，用于根据与ID编码对应的计算信号、位置信息以及发射功率的大小，采用RSS方法，计算出接收机的水平坐标。
在本发明的LED定位系统的水平坐标计算子模块中，接收机的水平坐标通过下述公式计算得出：
Q＝(MTM)-1MTN，
其中，且Q为所求的接收机坐标(x,y)，Xi、Yi为第i个发射机的水平面的横坐标和纵坐标，ri为接收机与第i个发射机之间的水平距离，i＝1、2、…、n，n为发射机的个数，ri为通过下述公式所描述：
式中，Pt为第i个LED的重构后的PAM-DMT信号的功率的取平均值，Pri为第i个LED的PAM-DMT信号的发射功率的平均值，h为接收机与接收到的发射光所对应LED灯的垂直距离，朗伯阶数为m＝-ln2/ln(cosΨ1/2)，Ψ1/2是用于发射信号的LED的半功率角，A是用于接收信号的光电探测器的物理面积，Di为接收机与第i个发射机之间的距离，θi为第i盏灯发射角，为第i个接收机相对于其法线方向的入射角，为第i个LED的光滤波器的传播系数，是第i个LED的光集中器的增益，大小由式下述公式确定：
其中，Φc为光集中器视角的半角度，n为光集中器的折射率。
在本发明的LED定位系统中，接收端处理模块还包括3D坐标计算子模块，3D坐标计算子模块包含：
初始化单元，用于初始化接收机坐标的初始值h0为接收机与接收到的发射光所对应LED灯的垂直距离h的初始值，H为接收到的发射光所对应LED灯所在房间的高度，其中，接收机的水平面中坐标(x,y)通过下述步骤B2得出：
更新单元，用于根据计算得到初始定位误差根据初始信赖域半径R0，取以为圆心、以R0为半径的信赖域内部的点更新通过更新
其中定位误差计算公式为：
(Xi,Yi,Zi)为接收到的发射光所对应的第i个LED灯的三维坐标，为计算的接收机的三维坐标；
定位误差判决单元，用于当定位误差大于设定的定位误差时，重新选取信赖域内部的点更新当定位误差小于等于设定的定位误差时，输出最新的
迭代终止判单元，用于将最新的中与上一次迭代得到的中的进行比较，当二者差值大于预设值时，以最新的中的确定下一次和更新迭代的初始值h0，返回步骤B31开始下一次迭代，否则终止迭代，输出最终的作为接收机的三维位置坐标。
在本发明的LED定位系统的接收估计子模块中，包含对接收到的每个LED发射机的信号分别进行如下处理的单元：
对接收到的信号进行数模转换，然后移除循环前缀，对移除循环前缀后的信号进行椭圆ACO-OFDM逆变换，对变换后的结果进行FFT变换，然后估算出信号。
在本发明的LED定位系统的接收估计子模块中，发射机的信号发射方向和接收机的信号接收方向均于竖直方向。
本发明公开了一种结合椭圆ACO-OFDM的LED定位方法及系统，用于减轻室内环境中多径失真的影响，LED用作具有唯一ID的发射机以区分其自己的位置，并且移动接收机基于RSS定位算法以获得定位自身在水平位置的定位，并进一步采用信赖域技术实现三维定位，定位精度高。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1是本发明的结合椭圆ACO-OFDM的LED定位方法的一实施例的流程图；
图2是本发明的结合椭圆ACO-OFDM的LED定位系统模型图；
图3是本发明的椭圆映射示意图；
图4是本发明的结合椭圆ACO-OFDM的LED定位方法的另一实施例流程图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
参考图1及图2，本实施例中的结合椭圆ACO-OFDM(EA-OFDM)的LED定位系统可适用于室内定位，下述将结合图2的室内使用情况对图1中的LED定位方法进行说明。
图2中，长方体为一个长6米、宽6米、高4.2米的房间，房间的天花板上有四个LED灯泡，每个灯作为一个独立的光发射机，发射唯一的ID(身份识别)编码及位置坐标，本实施例的位置坐标采用坐标(x,y，z)表示，x表示x轴的坐标，y表示y轴的坐标，z表示z轴的坐标，除包含在水平面的坐标(x,y)外，还包含高度坐标z。
在发射端的每个LED发射机中，包含LED的ID编码及位置信息的输入信号被M元正交调幅(M-ary quadrature amplitude modulation,M-QAM)调制在奇载波上，生成长度为N/4的复数数据，偶载波被全部置零，其中，M为星座图的大小，N为子载波数。此外，为了满足埃尔米特对称法则以确保输出信号为真，第一和第N/2个子载波也置零。此时，生成的信号为该信号为N点快速傅里叶反变换(inverse fastFourier transform,IFFT)的输入信号，其中()*为表示向量的复数共轭。
经过N点IFFT变换后，可以得到时域信号为c＝{caco[0],caco[1],caco[2],…,caco[N-1]}，该信号用于映射到椭圆上。如图3所示，椭圆的相位信息可由式(1)计算得到：
其中，α为一个大小为的换算因子，将约束在区间[-π/β,π/β]内，σc为caco[n]的标准偏差，β为一个大于2的正数。(re[n],im[n])表示在椭圆上的映射坐标，可由式(2)计算得到：
其中a和b分别为当焦点落在实轴时椭圆的长轴与短轴，由此，(re[n],im[n])可以转化为：
在本实施例的EA-OFDM系统中，只传输虚部通过零值限幅获得非负传输信号后，加入循环前缀(CP)用于削弱码间串扰的影响。
接收机对接收到的每个LED发射机的信号分别进行如下处理：对接收到的信号进行数模转换，然后移除循环前缀，对移除循环前缀后的信号进行椭圆ACO-OFDM逆变换，对变换后的结果进行FFT变换，然后估算出信号，获取各LED发射机的ID编码及位置信息并分别计算出各LED发射机的计算信号，根据与ID编码对应的计算信号、位置信息以及发射功率的大小，采用RSS方法，计算出接收机的水平坐标(x，y)。如何计算出接收机的水平坐标(x，y)可采用下述下述2D算法实现：
获取重构后的第i个LED的PAM-DMT信号的功率的取平均值Pt、第i个LED的PAM-DMT信号的发射功率的平均值Pri、朗伯阶数m＝-ln2/ln(cosΨ1/2)、用于发射信号LED的半功率角Ψ1/2、用于接收信号的光电探测器的物理面积A，Di为接收机与第i个发射机之间的距离，θi为第i盏LED发射角，为第i个接收机相对于其法线方向的入射角，为第i个LED的光滤波器的传播系数，是第i个LED的光集中器的增益，大小由式下述公式确定：
其中，Φc为光集中器视角的半角度，n为光集中器的折射率。
然后根据下述公式计算出ri：
式中，h为接收机与接收到的发射光所对应LED灯的垂直距离h。
然后计算出接收机在水平面中的坐标：
Q＝(MTM)-1MTN，
其中，且Q为所求的接收机坐标，Xi、Yi为第i个发射机的水平面的横坐标和纵坐标，ri为接收机与第i个发射机之间的水平距离，i＝1、2、…、n，n为发射机的个数。
在上述Di 2的计算公式中，Ψ1/2、m、Φc、Ts、n、A为系统参数，其一般通过人为预先设定得到，而cos(θi)、会根据发射机的发射器的发射方向和接收机的接收器的接收方向的不同而不同。
在本实施例中，发射机的发射器的发射方向和接收机的接收器的接收方向都与天花板垂直，此时有：
由于高度h固定，可得到收发机之间的距离的简化公式为：
由上式(5)可得：
其中，C为将下表中系统参数在公式(4)中计算相应部分所得的常数，其中C具体的计算式如下：
在本实施例中，部分参数的大小如下表所示。
上述介绍了在h固定情况下利用2D算法求解接收机的水平坐标(x,y)，下述在上述2D算法求解(x,y)上，进一步求解z坐标，具体的可参照图4，包含：
A1、初始化接收机坐标，的初始值其中h0为接收机与接收到的发射光所对应LED灯的垂直距离h的初始值，H为接收到的发射光所对应LED灯所在房间的高度；
A2、根据计算得到初始定位误差根据初始信赖域半径R0，取以为圆心、以R0为半径的信赖域内部的点更新通过更新
其中定位误差计算公式为：
(Xi,Yi,Zi)为接收到的发射光所对应的第i个LED灯的三维坐标，为计算的接收机的三维坐标；
A3、进行定位误差判决，当定位误差大于设定的定位误差时，重新选取信赖域内部的点更新当定位误差小于等于设定的定位误差时，输出最新的
A4、进行迭代终止判决，将最新的中与上一次迭代得到的中的进行比较，当二者差值大于预设值时，以最新的中的确定下一次和更新迭代的初始值h0，返回步骤A1开始下一次迭代，否则终止迭代，输出最终的作为接收机的三维位置坐标。再次返回A1时，需要利用新的h0再次进行2D算法求解(x，y)，并利用新的(x，y)进行后续的步骤。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。
Claims (10)
1.一种结合椭圆ACO-OFDM的LED定位方法，其特征在于，包含：
发射端处理步骤：
A1、各LED发射机分别将包含自身ID编码及位置信息的输入信号分别进行M元正交调幅调制在奇载波上，生成长度为N/4的复数数据，偶载波被全部置零，其中，M为星座图的大小，N为子载波数，并将第一和第N/2个子载波也置零，得到频域信号()*为表示向量的复数共轭；
A2、各LED发射机分别将各自的频域信号C作为输入进行N点IFFT变换，得到时域信号对应的时域信号c＝{caco[0],caco[1],caco[2],…,caco[N-1]}；
A3、各LED发射机分别将各自的时域信号c＝{caco[0],caco[1],caco[2],…,caco[N-1]}发送出去，其中任意caco[n]在被发送时，只传送其对应的虚部其中n＝0、1、2…、N-1，θ[n]＝αcaco[n]为caco[n]的相位，α为一个大小为的换算因子以将θ[n]约束在区间[-π/β,π/β]内，σc为caco[n]的标准偏差，β为一个大于2的正数，a和b分别为caco[n]在当焦点落在实轴时椭圆的长轴与短轴；
接收端处理步骤：
B1、接接收各LED发射机发射的信号，获取各LED发射机的ID编码及位置信息并分别计算出各LED发射机的发射信号；
B2、根据各LED发射机的发射信号中与ID编码对应的位置信息以及发射功率的大小，采用RSS方法，计算出接收机的水平坐标。
2.根据权利要求1所述的LED定位方法，其特征在于，所述步骤B2中，接收机的水平坐标通过下述公式计算得出：
Q＝(MTM)-1MTN，
其中，且Q为所求的接收机坐标(x,y)，Xi、Yi为第i个发射机的水平面的横坐标和纵坐标，ri为接收机与第i个发射机之间的水平距离，i＝1、2、…、n，n为发射机的个数，ri为通过下述公式所描述：
<mrow>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mroot>
<mfrac>
<mrow>
<mi>C</mi>
<mo>&CenterDot;</mo>
<mi>h</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mn>4</mn>
</mroot>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>D</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>C</mi>
<mo>&CenterDot;</mo>
<mi>h</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>,</mo>
</mrow>
式中，Pt为第i个LED的重构后的PAM-DMT信号的功率的取平均值，Pri为第i个LED的PAM-DMT信号的发射功率的平均值，h为接收机与接收到的发射光所对应LED灯的垂直距离，朗伯阶数为m＝-ln2/ln(cosΨ1/2)，Ψ1/2是用于发射信号的LED的半功率角，A是用于接收信号的光电探测器的物理面积，Di为接收机与第i个发射机之间的距离，θi为第i盏灯发射角，为第i个接收机相对于其法线方向的入射角，为第i个LED的光滤波器的传播系数，是第i个LED的光集中器的增益，大小由式下述公式确定：
其中，Φc为光集中器视角的半角度，n为光集中器的折射率。
3.根据权利要求2所述的LED定位方法，其特征在于，接收端处理步骤还包括步骤B3，步骤B3包含：
B31、初始化接收机坐标的初始值h0为接收机与接收到的发射光所对应LED灯的垂直距离h的初始值，H为接收到的发射光所对应LED灯所在房间的高度，其中，接收机的水平面中坐标(x,y)通过下述步骤B2得出：
B32、根据计算得到初始定位误差根据初始信赖域半径R0，取以为圆心、以R0为半径的信赖域内部的点更新通过更新
其中定位误差计算公式为：
<mrow>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<munder>
<mi>argmin</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>,</mo>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>4</mn>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
(Xi,Yi,Zi)为接收到的发射光所对应的第i个LED灯的三维坐标，为计算时的接收机的三维坐标；
B33、进行定位误差判决，当定位误差大于设定的定位误差时，重新选取信赖域内部的点更新当定位误差小于等于设定的定位误差时，输出最新的
B34、进行迭代终止判决，将最新的中与上一次迭代得到的中的进行比较，当二者差值大于预设值时，以最新的中的确定下一次和更新迭代的初始值h0，返回步骤B31开始下一次迭代，否则终止迭代，输出最终的作为接收机的三维位置坐标。
4.根据权利要求1所述的LED定位方法，其特征在于，所述步骤B1中，包含对接收到的每个LED发射机的信号分别进行如下处理的步骤：
对接收到的信号进行数模转换，然后移除循环前缀，对移除循环前缀后的信号进行椭圆ACO-OFDM逆变换，对变换后的结果进行FFT变换，然后估算出信号。
5.根据权利要求1所述的LED室内3D定位方法，其特征在于，发射机的信号发射方向和接收机的信号接收方向均于竖直方向。
6.一种结合椭圆ACO-OFDM的LED定位系统，其特征在于，包含：
发射端处理模块，包含：
输入调制子模块，用于将各LED发射机分别将包含自身ID编码及位置信息的输入信号分别进行M元正交调幅调制在奇载波上，生成长度为N/4的复数数据，偶载波被全部置零，其中，M为星座图的大小，N为子载波数，并将第一和第N/2个子载波也置零，得到频域信号()*为表示向量的复数共轭；
时域变换子模块，用于将各LED发射机分别将各自的频域信号C作为输入进行N点IFFT变换，得到时域信号对应的时域信号c＝{caco[0],caco[1],caco[2],…,caco[N-1]}；
虚部处理子模块，用于将各LED发射机分别将各自的时域信号c＝{caco[0],caco[1],caco[2],…,caco[N-1]}发送出去，其中任意caco[n]在被发送时，只传送其对应的虚部其中n＝0、1、2…、N-1，θ[n]＝αcaco[n]为caco[n]的相位，α为一个大小为的换算因子以将θ[n]约束在区间[-π/β,π/β]内，σc为caco[n]的标准偏差，β为一个大于2的正数，a和b分别为caco[n]在当焦点落在实轴时椭圆的长轴与短轴；
接收端处理模块，包含：
接收估计子模块，用于接收各LED发射机发射的信号，获取各LED发射机的ID编码及位置信息并分别计算出各LED发射机的计算信号；
2D坐标计算子模块，用于根据各LED发射机的发射信号中与ID编码对应的位置信息以及发射功率的大小，采用RSS方法，计算出接收机的水平坐标。
7.根据权利要求1所述的LED定位系统，其特征在于，所述水平坐标计算子模块中，接收机的水平坐标通过下述公式计算得出：
Q＝(MTM)-1MTN，
其中，且Q为所求的接收机坐标(x,y)，Xi、Yi为第i个发射机的水平面的横坐标和纵坐标，ri为接收机与第i个发射机之间的水平距离，i＝1、2、…、n，n为发射机的个数，ri为通过下述公式所描述：
<mrow>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mroot>
<mfrac>
<mrow>
<mi>C</mi>
<mo>&CenterDot;</mo>
<mi>h</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mn>4</mn>
</mroot>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>D</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>C</mi>
<mo>&CenterDot;</mo>
<mi>h</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>i</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>,</mo>
</mrow>
式中，Pt为第i个LED的重构后的PAM-DMT信号的功率的取平均值，Pri为第i个LED的PAM-DMT信号的发射功率的平均值，h为接收机与接收到的发射光所对应LED灯的垂直距离，朗伯阶数为m＝-ln2/ln(cosΨ1/2)，Ψ1/2是用于发射信号的LED的半功率角，A是用于接收信号的光电探测器的物理面积，Di为接收机与第i个发射机之间的距离，θi为第i盏灯发射角，为第i个接收机相对于其法线方向的入射角，为第i个LED的光滤波器的传播系数，是第i个LED的光集中器的增益，大小由式下述公式确定：
其中，Φc为光集中器视角的半角度，n为光集中器的折射率。
8.根据权利要求7所述的LED定位系统，其特征在于，接收端处理模块还包括3D坐标计算子模块，3D坐标计算子模块包含：
初始化单元，用于初始化接收机坐标的初始值h0为接收机与接收到的发射光所对应LED灯的垂直距离h的初始值，H为接收到的发射光所对应LED灯所在房间的高度，其中，接收机的水平面中坐标(x,y)通过下述步骤B2得出：
更新单元，用于根据计算得到初始定位误差根据初始信赖域半径R0，取以为圆心、以R0为半径的信赖域内部的点更新通过更新
其中定位误差计算公式为：
<mrow>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<munder>
<mi>argmin</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>,</mo>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>4</mn>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
(Xi,Yi,Zi)为接收到的发射光所对应的第i个LED灯的三维坐标，为计算的接收机的三维坐标；
定位误差判决单元，用于当定位误差大于设定的定位误差时，重新选取信赖域内部的点更新当定位误差小于等于设定的定位误差时，输出最新的
迭代终止判单元，用于将最新的中与上一次迭代得到的中的进行比较，当二者差值大于预设值时，以最新的中的确定下一次和更新迭代的初始值h0，返回步骤B31开始下一次迭代，否则终止迭代，输出最终的作为接收机的三维位置坐标。
9.根据权利要求1所述的LED定位系统，其特征在于，所述接收估计子模块中，包含对接收到的每个LED发射机的信号分别进行如下处理的单元：
对接收到的信号进行数模转换，然后移除循环前缀，对移除循环前缀后的信号进行椭圆ACO-OFDM逆变换，对变换后的结果进行FFT变换，然后估算出信号。
10.根据权利要求1所述的LED室内3D定位系统，其特征在于，发射机的信号发射方向和接收机的信号接收方向均于竖直方向。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710979861.2A CN107883955A (zh) | 2017-10-19 | 2017-10-19 | 结合椭圆aco‑ofdm的led定位方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710979861.2A CN107883955A (zh) | 2017-10-19 | 2017-10-19 | 结合椭圆aco‑ofdm的led定位方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107883955A true CN107883955A (zh) | 2018-04-06 |
Family
ID=61782023
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710979861.2A Pending CN107883955A (zh) | 2017-10-19 | 2017-10-19 | 结合椭圆aco‑ofdm的led定位方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107883955A (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20130085024A (ko) * | 2012-01-18 | 2013-07-26 | 경희대학교 산학협력단 | 실내 가시광 통신을 이용한 3차원 위치 측정 방법 및 그 시스템과, 위치 기반 서비스 제공 방법 |
| CN106443585A (zh) * | 2016-09-09 | 2017-02-22 | 中国地质大学（武汉） | 一种结合加速度计的led室内3d定位方法 |
-
2017
- 2017-10-19 CN CN201710979861.2A patent/CN107883955A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20130085024A (ko) * | 2012-01-18 | 2013-07-26 | 경희대학교 산학협력단 | 실내 가시광 통신을 이용한 3차원 위치 측정 방법 및 그 시스템과, 위치 기반 서비스 제공 방법 |
| CN106443585A (zh) * | 2016-09-09 | 2017-02-22 | 中国地质大学（武汉） | 一种结合加速度计的led室内3d定位方法 |
Non-Patent Citations (3)
| Title |
|---|
| BILAL RANJHA ,ET AL: "Hybrid Asymmetrically Clipped OFDM-Based IM/DD Optical Wireless System", 《JOURNAL OF OPTICAL COMMUNICATIONS AND NETWORKING》 * |
| MOHAMMADREZA AMINIKASHANI,ET AL: "Indoor location estimation with optical-based orthogonal frequency division multiplexing communications", 《OPTICAL ENGINEERING》 * |
| TIANQI MAO，ET AL: "Ellipse-based DCO-OFDM for visible light communications", 《OPTICS COMMUNICATIONS》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107769851A (zh) | 一种基于混合ofdm的led定位方法与系统 | |
| CN110830097B (zh) | 一种基于反射面的主被动互惠共生传输通信系统 | |
| CN102421189B (zh) | 一种基于物理层的无线网络室内定位的实现方法及其装置 | |
| Aminikashani et al. | Indoor positioning with OFDM visible light communications | |
| CN104813189B (zh) | 减小用于位置感知的接收信号强度指示符（rssi）测量的变化的系统和方法 | |
| CN110224965A (zh) | 一种基于ofdm反向散射通信系统半盲接收机设计方法 | |
| CN104618300B (zh) | 双重限制下的dco-ofdm直流偏置和功率快速优化方法 | |
| CN104181503A (zh) | 一种基于rss的led定位系统及方法 | |
| CN102299894A (zh) | 一种基于叠加周期序列的aco-ofdm无线光通信系统信道估计方法 | |
| CN112751798B (zh) | 用于太赫兹通信感知一体化的nmw-ofdm系统 | |
| CN104753836A (zh) | 一种载波频偏估计方法 | |
| CN107167767A (zh) | 一种基于sc‑fdma的可见光室内高精度定位方法 | |
| CN107942288A (zh) | 一种结合ofdm的led定位方法及系统 | |
| CN109061557A (zh) | 室内可见光通信定位集成方法 | |
| CN101834705B (zh) | 一种无线通信系统的信号映射与解映射方法 | |
| Wang et al. | VLC-based indoor positioning algorithm combined with OFDM and particle filter | |
| CN107018105A (zh) | 水下可见光通信方法及系统 | |
| Aminikashani et al. | Indoor location estimation with optical-based orthogonal frequency division multiplexing communications | |
| CN107883955A (zh) | 结合椭圆aco‑ofdm的led定位方法及系统 | |
| CN115087095B (zh) | 一种基于csi的加权knn的可见光室内定位方法 | |
| CN104467915B (zh) | 一种非相干超宽带通信系统的相位噪声分析方法 | |
| CN101032076B (zh) | 估算接收信号大小、噪声功率及信号噪声比的方法 | |
| CN114374435B (zh) | 一种基于ofdm的可见光通信与定位一体化方法及系统 | |
| Choi et al. | Influence of optical path difference on visible light communication systems | |
| CN107040313B (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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180406 |
|
| WD01 | Invention patent application deemed withdrawn after publication |