一种室内无线光定位方法、装置及系统

1.本发明涉及室内定位技术领域，尤其是指一种室内无线光定位方法、装置及系统。


背景技术：

2.随着智能手机数量的不断增长，物联网、虚拟现实等技术的不断发展，需要追踪室内环境中密切接触人的位置，高精度室内定位的变得越来越重要。用于室外定位的卫星定位技术，由于其信号在穿越墙壁时会经历严重的衰减，无法在室内环境中提供精准定位。近年来，许多新型室内定位技术不断发展，包括wi-fi定位、超声波传感器定位、蓝牙定位、rfid定位、基于无线光通信的可见光定位等。
3.基于指纹库识别法的可见光定位方案是一种常见的方案，该方案预先采集各个待定位点的特征信息，建立成库。在定位过程中，通过比对采集到的特征信息和库中存储的特征信息，进行定位。基于指纹库识别法的定位方法分为两个过程，第一个过程是在定位空间中选择一些特定位置作为参考点，每个参考点的特征不同，用这些特征构建指纹库。第二个过程是在实际定位时，用待定位点的特征和指纹库中的特征对比，从而估算出带定位点的位置。这种方法，前期构建指纹库的工作十分繁琐，耗时长，且应用环境是特定的，可移植性差。
4.另外一种广泛应用的定位方案是几何测量法。这种定位方法又分为基于角度的定位方法和基距离的定位方法。基于到达角度(aoa)的定位方法需要测量发射端和接收端之间的夹角。基于距离的定位方法的关键是找到发射端和接收端的距离关系。距离关系可以通过到达时间(toa)，接收信号强度 (rss)等方法确定。几何测量法中基于距离的定位方法，在采用toa获取距离时，需要收发机之间或光源之间保持严格的时间同步，这增加了系统的实现复杂度。几何测量法中基于aoa的定位方法，需要对接收机阵列进行一定改造以获取角度，硬件复杂度较高，且较小的角度测量误差会带来较大的定位误差。
5.基于图像传感器的定位方案也是一种常用的方案。该方法的接收端是由图像传感器和一个透镜组成，透镜的中心就是待测点。由于不同led经过透镜后在图像传感器上显示的位置不同，通过对图像传感器上的led的成像和发射端led的几何关系分析，完成定位。算法复杂度较高，响应速率也相对较慢。
6.无线光室内定位技术以其成本低，无电磁干扰，安全性好，定位精度高
7.(可达厘米级别)等优点，在这些室内定位技术中脱颖而出。然而当前的无线光室内定位技术大多基于时分复用或者频分复用的方法，以此区分不同灯的信号，完成定位。时分复用至少需要三个时隙发送信号，频分复用至少需要三个不同的频率发送信号，并且未将通信和定位结合在一起考虑，只是单独考虑定位，功能单一，忽略了通信这一重要需求，浪费了通信资源。


技术实现要素：

8.为此，本发明所要解决的技术问题在于克服现有技术中未充分利用通信资源、可
移植性差、硬件复杂度高、成本大的问题。
9.为解决上述技术问题，本发明提供了一种室内无线光定位方法，包括：
10.s1:选取不同的两组均匀分布在室内天花板上的光源发射机，并至少保证每组中的光源发射机数量不小于3，且不在同一条直线上；
11.s2:根据第一组光源发射机在第一时隙中发送至接收机的通信信号，进行信道估计，并将信道估计的频域信息转变为时域冲激响应的形式，得到第一信道估计结果；
12.s3:根据所述第一信道估计结果，在时域中选取峰值最大的三个视距链路脉冲作为第一组脉冲，并根据所述第一组脉冲计算其对应的三个光源发射机所发出的信号，从发送到被所述接收机接收的时间之间的第一组时间差值；
13.s4:根据所述第一组时间差值，利用tdoa算法对所述接收机进行定位，求解得到第一组接收机坐标；
14.s5:对在第二时隙发送通信信号的第二组光源发射机重复s2-s4的步骤，得到第二组接收机坐标；
15.s6:找到所述第一组接收机坐标和所述第二组接收机坐标之间距离最近的一对坐标，并选择这对坐标中任意一个坐标作为接收机最终定位。
16.优选地，所述根据第一组光源发射机在第一时隙中发送至接收机的通信信号，进行信道估计，并将信道估计的频域信息转变为时域冲激响应的形式，得到第一信道估计结果包括：
17.当所述接收机接收到所述光源发射机发送的通信信号后，将所述通信信号转换为电信号并进行模数转换和滤除直流偏置，得到接收信号；
18.对加窗后的接收信号y(n)进行离散傅里叶变换，得到频域中的接收信号 y(k)＝dft[y(n)]＝x(k)*h(k) w(k)k＝1,2....n，其中，加窗长度为n，n表示dft 的点数，和子载波个数相同，y(k)表示第k个子载波接收的导频符号，x(k)表示第k个子载波发送的导频符号，w(k)是噪声；
[0019]
根据所述频域中的接收信号计算频域响应的信道估计根据所述频域中的接收信号计算频域响应的信道估计其中，是y(k)加上噪声后的估计值，j为虚数单位，q为加窗位置偏移量；
[0020]
取频域响应的前一半序列进行共轭对称：
[0021][0022]
并进行离散傅里叶逆变换得到时域冲激响应的信道估计：
[0023]
其中，
[0024]
表示的是时刻i对应的信道冲激响应，n为采样点数。
[0025]
优选地，所述根据第一组光源发射机在第一时隙中发送至接收机的通信信号，进行信道估计，并将信道估计的频域信息转变为时域冲激响应的形式，得到第一信道估计结果后还包括：
[0026]
对所述第一信道估计结果进行均值除噪。
[0027]
优选地，所述根据所述第一信道估计结果，在时域中选取峰值最大的三个视距链路脉冲作为第一组脉冲，并根据所述第一组脉冲计算其对应的三个光源发射机所发出的信号，从发送到被所述接收机接收的时间之间的第一组时间差值包括：
[0028]
根据所述第一信道估计结果，在时域中选取峰值最大的三个视距链路脉冲，将所述三个视距链路脉冲对应的三个光源发射机记为第一光源发射机l1，第二光源发射机l2，第三光源发射机l3，将所述三个视距链路脉冲对应的采样点下标记为第一采样点下标n1，第二采样点下标n2，第三采样点下标n3；
[0029]
假设n1》n2》n3，计算从l1发出通信信号到接收机接收到该通信信号的时间和l2发出通信信号到接收机接收到该通信信号的时间之间的时间差值δt1＝ (n
1-n2)*ts，其中，ts为采样周期；
[0030]
计算从l1发出通信信号到接收机接收到该通信信号的时间和l3发出通信信号到接收机接收到该通信信号的时间之间的时间差值δt2＝(n
1-n3)*ts。
[0031]
优选地，所述根据所述第一组时间差值，利用tdoa算法对所述接收机进行定位，求解得到第一组接收机坐标包括：
[0032]
根据所述第一组时间差值列写方程组：
[0033][0034]
其中，c表示光速，x
i i∈1，2，3表示l1，l2，l3的横坐标，y
i i∈1，2，3表示l1，l2，l3的纵坐标，z
i i∈1，2，3表示l1，l2，l3的高度坐标，xr表示接收机的横坐标，yr表示接收机的纵坐标，zr表示接收机的高度坐标；
[0035]
从第一组光源发射机中随机选出三个光源发射机，并进行排列组合，得到多组l1，l2，l3的三维坐标值；
[0036]
根据所述多组l1，l2，l3的三维坐标值，在只考虑二维定位的情况下，利用chan’s方法或fang’s方法求解所述方程组得到个解，即第一组接收机坐标，其中，m1为第一组光源发射机的数量。
[0037]
优选地，所述利用chan’s方法求解所述方程组包括：
[0038]
将所述方程组转换为线性方程组：
[0039][0040]
并进行变形得到xr，yr，r1之间的关系：
[0041][0042]
其中，r
2,1
＝c*δt1，表示l1到接收机的距离和l2到接收机的距离之差，r
3,1
＝ c*δt2，表示l1到接收机的距离和l3到接收机的距离之差，r1表示l1到接收机的距离；
[0043]
将所述xr，yr，r1之间的关系带入公式：
[0044]
[0045]
得到一个关于r1的一元二次方程进行求解得到r1的值，并带入所述线性方程组，求解得到接收机的坐标，其中，a，b，c1均为常数。
[0046]
优选地，所述利用fang’s方法求解所述方程组包括：
[0047]
建立一个新的坐标系，使得l1位于坐标原点(0，0)，使得l2位于(x2，0),l3 位于(x3，y3)，则l1到接收机的距离r1可表示为：
[0048][0049]
l2，l3到接收机的距离可表示为：
[0050][0051]
其中，r
i,1
表示li到接收机的距离和l1到接收机的距离之差；
[0052]
将式(1)代入式(2)中可得：
[0053][0054][0055]
根据式(3)和式(4)找到接收机横坐标和纵坐标之间的关系：
[0056]
yr＝g
·
xr p
ꢀꢀ
(5)
[0057]
其中，系数g＝(r
3,1
x2/r
2,1-x3)/y3，系数系数
[0058]
将式(5)代入式(3)，得到最终待求解的方程：对其进行求解得到接收机坐标；
[0059]
其中，系数d＝-{1-(x2/r
2,1
)2 g2}，系数e＝x2{1-(x2/r
2,1
)2}-2
gp
，系数系数
[0060]
本发明还提供了一种室内无线光定位装置，包括：
[0061]
光源发射机组选取模块，用于选取不同的两组均匀分布在室内天花板上的光源发射机，并至少保证每组中的光源发射机数量不小于3，且不在同一条直线上；
[0062]
信道估计模块，用于根据第一组光源发射机在第一时隙中发送至接收机的通信信号，进行信道估计，并将信道估计的频域信息转变为时域冲激响应的形式，得到第一信道估计结果；根据第二组光源发射机在第二时隙中发送至接收机的通信信号，进行信道估计，并将信道估计的频域信息转变为时域冲激响应的形式，得到第二信道估计结果；
[0063]
时间差提取模块，用于根据所述第一信道估计结果，在时域中选取峰值最大的三个视距链路脉冲作为第一组脉冲，并根据所述第一组脉冲计算其对应的三个光源发射机所发出的信号，从发送到被所述接收机接收的时间之间的第一组时间差值；根据所述第二信道估计结果，在时域中选取峰值最大的三个视距链路脉冲作为第二组脉冲，并根据所述第二组脉冲计算其对应的三个光源发射机所发出的信号，从发送到被所述接收机接收的时间之间的第二组时间差值；
[0064]
坐标计算模块，用于利用tdoa算法对所述接收机进行定位，根据所述第一组时间差值求解得到第一组接收机坐标，根据所述第二组时间差值求解得到第二组接收机坐标；
[0065]
定位模块，用于找到所述第一组接收机坐标和所述第二组接收机坐标之间距离最近的一对坐标，并选择这对坐标中任意一个坐标作为接收机最终定位。
[0066]
本发明还提供了一种室内无线光定位系统，包括：
[0067]
均匀安装在室内天花板上的多个光源发射机；
[0068]
如上述的室内无线光定位装置，用于利用所述多个光源发射机定位设置有接收机的物品或设备。
[0069]
优选地，所述室内无线光定位系统应用于地下停车场的汽车定位，所述光源发射机为led灯，所述接收机为安装于汽车上的光电探测器。
[0070]
本发明的上述技术方案相比现有技术具有以下优点：
[0071]
本发明所述的室内无线光定位方法，利用信道状态信息实现无线光室内定位，即利用两组灯在两个不同时隙通信时，对通信过程中的信道信息进行信道估计，之后通过到达时间差的方法(tdoa)进行定位，实现精确的无线光室内定位。本发明将无线光通信和定位相融合，利用信道估计得到的信息进行定位，充分利用了通信资源，两组光源使用时分复用的方式，轮流和接收机进行通信，只需要两个时隙，完成定位，使用时无需占用额外的通信资源来传输定位信息，节省了通信资源；可应用于不同的室内场景，只需通过信道估计得到时域冲激响应，无论室内环境有多复杂，都能实现精确的无线光室内定位；利用到达时间差，使用tdoa的方法进行定位，降低了对发射端和接收端的时间同步的需求，系统可移植性好，大大降低了硬件的复杂度，降低了成本。
附图说明
[0072]
为了使本发明的内容更容易被清楚的理解，下面根据本发明的具体实施例并结合附图，对本发明作进一步详细的说明，其中：
[0073]
图1是本发明所述的室内无线光定位的实现流程图；
[0074]
图2是一种实施例提供的室内无线光定位系统图；
[0075]
图3是一种实施例提供的信道的频域响应示意图；
[0076]
图4是一种实施例提供的时域冲激响应示意图；
[0077]
图5是二维定位结果示意图；
[0078]
图6为本发明实施例提供的一种室内无线光定位装置的结构框图。
具体实施方式
[0079]
本发明的核心是提供一种室内无线光定位方法、装置及系统，充分利用了通信资源，可移植性好，大大降低了硬件的复杂度，降低了成本。
[0080]
为了使本技术领域的人员更好地理解本发明方案，下面结合附图和具体实施方式对本发明作进一步的详细说明。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0081]
在通信的过程中，每组光源发射机在不同的时隙中和接收机通信，向接收机发送通信信息。接收机利用接收到的通信信息，进行信道估计。我们对信道估计的结果进行离散傅里叶逆变换(idft)，将信道估计的频域信息变为时域冲激响应的形式。在时域中，选取峰值最大的三个视距链路脉冲，利用这三个脉冲，计算出不同光源发射机到接收机的时间差。之后通过到达时间差的方法(tdoa)进行定位，算出第一组解。接着利用第二个时隙进行同
样的过程，计算出第二组解。找到这两组解之间最接近的一个坐标，该坐标就是室内定位的结果。
[0082]
请参考图1，图1为本发明所提供的室内无线光定位方法的实现流程图；具体操作步骤如下：
[0083]
s101:选取不同的两组均匀分布在室内天花板上的光源发射机，并至少保证每组中的光源发射机数量不小于3，且不在同一条直线上；
[0084]
发射端光源均匀分布在天花板上，随着房间的增大，适当增加光源，以同时满足通信、定位和照明的需求；
[0085]
将m个光源分成两组，第一组有光源m1个，第二组有光源m2个。这两组光源使用时分复用的方式，轮流和接收机进行通信，只需要两个时隙，完成定位；
[0086]
在使用tdoa进行定位的过程中，在时隙1时，从m1个光源中选取3个光源，同样的，时隙2中，从m2个光源中选取3个光源，因此，每组中的光源发射机数量不小于3，且不在同一条直线上，以保证tdoa定位的方法是有解的，此外两次时隙选取的光源不能完全一样且需要能保证排除对称解。
[0087]
s102:根据第一组光源发射机在第一时隙中发送至接收机的通信信号，进行信道估计，并将信道估计的频域信息转变为时域冲激响应的形式，得到第一信道估计结果；
[0088]
随机生成一组发射端(即光源发射机)和接收端(即接收机)都已知的长度为n/2-1的二进制序列，其中，n为ofdm子载波个数；
[0089]
对所述二进制序列进行调制，得到导频序列序列中的元素表示调制后的一个符号，如果调制方式是qam调制，就是1个qam符号；
[0090]
对所述导频序列进行共轭对称，生成共轭对成序列：
[0091][0092]
对所述共轭对成序列进行离散傅里叶逆变换得到导频序列在时域中对应的离散序列，并对所述离散序列加入循环前缀，得到发射端ofdm符号x
cp
(n)，对所述发射端ofdm符号进行数模(d/a)转换，得到发射端ofdm符号的波形x(t)，对所述波形加上直流偏置，生成用于调制光源发射机光强度的第一电信号x
dc
(t)，将所述第一电信号经过光源发射机生成光信号，所述光信号即为所述通信信号；
[0093]
所述通信信号在空间中的信道传输，在接收端由光电转换器将光信号转换为第二电信号其中，w(t)是噪声的强度，r是光电转换系数，h(t)是无线光信道的传递函数，可表示为：其中，n
l
表示通信过程中使用的光源发射机个数，由于之后的求解是分两个时隙，在两个时隙下，n
l
的值分别为m1和m2，hi(0)表示第i个灯到接收机的直流增益， c表示光速，di表示第i个灯到接收机的距离，t表示时间，发射端开始发送信息的时刻时间为0，在时，h(t)不为0，是因为有个发射端的信息被接收端接收到了。在没有噪声的情况下，其余时间t，h(t)均为0；
[0094]
在接收端对y
dc
(t)进行模数(a/d)转换，以ts为采样周期进行采样(t)进行模数(a/d)转换，以ts为采样周期进行采样其中，为卷积符号，n为采样点的个数。
[0095]
为了方便表述，将光电转换系数置为1。在接收端滤除直流偏置，得到接收信号y
cp
(n)；
[0096]
对y
cp
(n)加窗(假设接收端加窗位置和发射端是完美的)，对加窗之后的信号y(n)进行离散傅里叶变换dft，将时域离散的信号变到频域中处理，得到频域中的接收信号:y(k)＝dft[y(n)]＝x(k)*h(k) w(k) k＝1,2....n，其中，加窗长度为n(实际情况中加窗位置不可能是完美的，一定存在偏移。由于循环卷积， y(n)可以看成以长度n为周期的序列，所以只要y(n)加窗的长度为n，就能保证截取一个周期的信息。只要加窗长度为n，即使存在偏移，也不影响频域响应的频谱的幅值，只影响频域响应的相位。对频域响应相位的影响，最终在时域中表现为信道的时域冲激响应的整体偏移，这是由dft的时域循环位移定理确定的)，n表示dft的点数，和子载波个数相同，y(k)表示第k个子载波接收的导频符号，x(k)表示第k个子载波发送的导频符号，w(k)是噪声；
[0097]
根据所述频域中的接收信号计算频域响应的信道估计根据所述频域中的接收信号计算频域响应的信道估计其中，是y(k)加上噪声后的估计值，j为虚数单位，q为加窗位置偏移量；
[0098]
对于式其中第一项和n/2项，由于发送的序列为0，除以0会是无穷大，这些位置在信道估计时置为0，取频域响应的前一半序列进行共轭对称：
[0099][0100]
并进行离散傅里叶逆变换得到时域冲激响应的信道估计：
[0101]
其中，
[0102]
表示的是时刻i对应的信道冲激响应，n为采样点数。
[0103]
s103:根据所述第一信道估计结果，在时域中选取峰值最大的三个视距链路脉冲作为第一组脉冲，并根据所述第一组脉冲计算其对应的三个光源发射机所发出的信号，从发送到被所述接收机接收的时间之间的第一组时间差值；
[0104]
假设本系统有m个光源发射机，第一通信时隙中的m1个光源发射机到接收机的时域冲激响应的估计值，是包含噪声的。m1个光源发射机到接收机的时域冲激响应表现为幅值较高的脉冲信号，噪声表现为幅值较低的脉冲信号。因此，对第一信道估计结果进行均值除噪，此时每个灯到接收机的时域冲激响应和噪声的幅值比例相差更大。均值除噪后选取幅值最大的三个视距链路脉冲，这三个脉冲是最接近接收机的三个光源发射机发出的。
[0105]
根据所述第一信道估计结果，在时域中选取峰值最大的三个视距链路脉冲，将所述三个视距链路脉冲对应的三个光源发射机记为第一光源发射机l1，第二光源发射机l2，第三光源发射机l3(目前还不知道这三个灯对应的是m1个光源发射机中的哪三个光源发射机)，将所述三个视距链路脉冲对应的采样点下标记为第一采样点下标n1，第二采样点下标n2，第三采样点下标n3；
[0106]
假设n1》n2》n3(实际情况中，可能不是n1》n2》n3，因此需要通过排列组合，得到多种
情况下的解)，计算从l1发出通信信号到接收机接收到该通信信号的时间和l2发出通信信号到接收机接收到该通信信号的时间之间的时间差值δt1＝(n
1-n2)*ts，其中，ts为采样周期；
[0107]
计算从l1发出通信信号到接收机接收到该通信信号的时间和l3发出通信信号到接收机接收到该通信信号的时间之间的时间差值δt2＝(n
1-n3)*ts。
[0108]
s104:根据所述第一组时间差值，利用tdoa算法对所述接收机进行定位，求解得到第一组接收机坐标；
[0109]
对在第二时隙发送通信信号的第二组光源发射机重复s102-s104的步骤，得到第二组接收机坐标；
[0110]
s105:找到所述第一组接收机坐标和所述第二组接收机坐标之间距离最近的一对坐标，并选择这对坐标中任意一个坐标作为接收机最终定位。
[0111]
本发明所述的室内无线光定位方法，利用信道状态信息实现无线光室内定位，即利用两组灯在两个不同时隙通信时，对通信过程中的信道信息进行信道估计，之后通过到达时间差的方法(tdoa)进行定位，实现精确的无线光室内定位。本发明将无线光通信和定位相融合，利用信道估计得到的信息进行定位，充分利用了通信资源，两组光源使用时分复用的方式，轮流和接收机进行通信，只需要两个时隙，完成定位，使用时无需占用额外的通信资源来传输定位信息，节省了通信资源；可应用于不同的室内场景，只需通过信道估计得到时域冲激响应，无论室内环境有多复杂，都能实现精确的无线光室内定位；利用到达时间差，使用tdoa的方法进行定位，降低了对发射端和接收端的时间同步的需求，系统可移植性好，大大降低了硬件的复杂度，降低了成本；可扩展性好，通过增加灯的数量，即可稳定扩大定位的范围。
[0112]
基于以上实施例，本实施例对步骤s104-s105进行进一步详细说明：
[0113]
根据所述第一组时间差值列写方程组：
[0114][0115]
其中，c表示光速，x
i i∈1，2，3表示l1，l2，l3的横坐标，y
i i∈1，2，3表示l1，l2，l3的纵坐标，z
i i∈1，2，3表示l1，l2，l3的高度坐标，xr表示接收机的横坐标，yr表示接收机的纵坐标，zr表示接收机的高度坐标，本系统暂时只考虑二维定位，所以方程组中z1，z2，z3，zr是已知的，且由于灯在同一天花板上，所以z1，z2，z3相等；
[0116]
从第一组光源发射机中随机选出三个光源发射机，并进行排列组合，得到多组l1，l2，l3的三维坐标值；
[0117]
到达时间差的方法本质上是求双曲线的交点，根据所述多组l1，l2，l3 的三维坐标值，在只考虑二维定位的情况下，利用chan’s方法或fang’s方法求解所述方程组得到个解，即第一组接收机坐标，其中，m1为第一组光源发射机的数量。
[0118]
利用chan’s方法求解所述方程组包括：
[0119]
将所述方程组转换为线性方程组：
[0120]
[0121]
并进行变形得到xr，yr，r1之间的关系：
[0122][0123]
其中，r
2,1
＝c*δt1，表示l1到接收机的距离和l2到接收机的距离之差，r
3,1
＝ c*δt2，表示l1到接收机的距离和l3到接收机的距离之差，经过计算r
2,1
和r
3,1
都是已知量，r1表示l1到接收机的距离；
[0124]
将所述xr，yr，r1之间的关系带入公式：
[0125][0126]
得到一个关于r1的一元二次方程进行求解得到r1的值，并带入所述线性方程组，求解得到接收机的坐标(每一组光源发射机的三维坐标可求得两个解(x
r1
,y
r1
)和(x
r2
,y
r2
))，其中，a，b，c1均是由光源发射机坐标和r
2,1
，r
3,1
求出的常数。
[0127]
利用fang’s方法求解所述方程组包括：
[0128]
建立一个新的坐标系，使得l1位于坐标原点(0，0)，使得l2位于(x2，0),l3 位于(x3，y3)，则l1到接收机的距离r1可表示为：
[0129][0130]
l2，l3到接收机的距离可表示为：
[0131][0132]
其中，r
i,1
表示li到接收机的距离和l1到接收机的距离之差；
[0133]
将式(1)代入式(2)中，经过计算合并同类项，可得：
[0134][0135][0136]
根据式(3)和式(4)找到接收机横坐标和纵坐标之间的关系：
[0137]
yr＝g
·
xr p
ꢀꢀ
(5)
[0138]
其中，系数g＝(r
3,1
x2/r
2,1-x3)/y3，系数系数均为已知量；
[0139]
将式(5)代入式(3)，得到最终待求解的方程：对其进行求解得到接收机坐标(每一组光源发射机的三维坐标可求得两个解(x
r1
,y
r1
) 和(x
r2
,y
r2
))；
[0140]
其中，系数d＝-{1-(x2/r
2,1
)2 g2}，系数e＝x2{1-(x2/r
2,1
)2}-2
gp
，系数系数均为已知量。
[0141]
chan’s和fang’s这两种方法算出的接收机坐标都有两个值，记为(x
r1
,y
r1
) 和(x
r2
,y
r2
)。由于在实际情况下并不确定距离接收机最近的三个光源发射机是天花板上的哪三个光源发射机，即l1，l2，l3的具体坐标可能是天花板上任意一个光源发射机的坐标。在m1个光源发射机中任意的选出3个光源发射机来对应l1，l2，l3，灯排列组合有种，每种灯的排列有两个解，所以在通信的第一时隙总共有个解：
[0142]
同样的，在在通信的第二时隙总共有个解：个解：
[0143]
两组解的最近两个坐标点，即为真实解:
[0144][0145]
基于以上实施例，本实施例为了评估所提出的无线光定位系统及方法的性能，在一个具体的房间内进行定位：
[0146]
房间尺寸为5m
×
5m
×
3m(长
×
宽
×
高)。光源均匀分布在天花板上，为led 灯，接收机为光电探测器，如图2中所示。在这次示例中，令m1＝3，选取led1、 led7、led9为第一组灯，在时隙t1中和接收机pd进行通信。同时令m2＝3，led1、 led4、led6为第二组灯，在时隙t2中和接收机pd进行通信。在通信时，需要进行信道估计，信道估计的关键仿真参数如下：
[0147]
房间大小5m
×
5m
×
3m(长
×
宽
×
高)时间分辨率0.1nsfft点数1024带宽125mhz导频序列个数100
[0148]
图3显示了，时隙t1下，接收机位置在(2m，1.5m，1m)的位置时，led1、 led7、led9和接收机pd进行通信时，信道的频域响应。图4显示了该频域响应进行idft后，对应的时域冲激响应。根据时域冲激响应，提取出用于定位的时间差值。利用时间差值，结合tdoa的方法，计算出12个可能解。再通过第二组灯在时隙t2和pd进行通信，重复t1时隙的步骤，也得到12个可能解。最后进行对称解的排除，最终实现对接收机的定位。
[0149]
我们在不同的位置进行测试，待测点所在的接收平面高度为1m，这些位置的间隔均为0.5m，如图5中的十字线所示。圆圈表示的是用提出的无线光定位系统和方法估计出的结果，在这次示例中，平均定位误差为2.13厘米。
[0150]
请参考图6，图6为本发明实施例提供的一种室内无线光定位装置的结构框图；具体装置可以包括：
[0151]
光源发射机组选取模块100，用于选取不同的两组均匀分布在室内天花板上的光源发射机，并至少保证每组中的光源发射机数量不小于3，且不在同一条直线上；
[0152]
信道估计模块200，用于根据第一组光源发射机在第一时隙中发送至接收机的通信信号，进行信道估计，并将信道估计的频域信息转变为时域冲激响应的形式，得到第一信道估计结果；根据第二组光源发射机在第二时隙中发送至接收机的通信信号，进行信道估计，并将信道估计的频域信息转变为时域冲激响应的形式，得到第二信道估计结果；
[0153]
时间差提取模块300，用于根据所述第一信道估计结果，在时域中选取峰值最大的三个视距链路脉冲作为第一组脉冲，并根据所述第一组脉冲计算其对应的三个光源发射机所发出的信号，从发送到被所述接收机接收的时间之间的第一组时间差值；根据所述第二信道估计结果，在时域中选取峰值最大的三个视距链路脉冲作为第二组脉冲，并根据所述第二组脉冲计算其对应的三个光源发射机所发出的信号，从发送到被所述接收机接收的时间之间的第二组时间差值；
[0154]
坐标计算模块400，用于利用tdoa算法对所述接收机进行定位，根据所述第一组时间差值求解得到第一组接收机坐标，根据所述第二组时间差值求解得到第二组接收机坐标；
[0155]
定位模块500，用于找到所述第一组接收机坐标和所述第二组接收机坐标之间距离最近的一对坐标，并选择这对坐标中任意一个坐标作为接收机最终定位。
[0156]
本实施例的室内无线光定位装置用于实现前述的室内无线光定位方法，因此室内无线光定位装置中的具体实施方式可见前文室内无线光定位方法的实施例部分，例如，光源发射机组选取模块100，信道估计模块200，时间差提取模块300，坐标计算模块400，定位模块500，分别用于实现上述无线光定位方法中步骤s101，s102，s103，s104和s105，所以，其具体实施方式可以参照相应的各个部分实施例的描述，在此不再赘述。
[0157]
本发明还提供了一种室内无线光定位系统，包括：
[0158]
均匀安装在室内天花板上的多个光源发射机；
[0159]
如上述的室内无线光定位装置，用于利用所述多个光源发射机定位设置有接收机的物品或设备。
[0160]
该室内无线光定位系统，可应用于地下停车场的汽车定位，所述光源发射机为led灯，所述接收机为安装于汽车上的光电探测器。
[0161]
将接收机集成在手机中，该室内无线光定位系统还可应用于室内人群追踪。
[0162]
显然，上述实施例仅仅是为清楚地说明所作的举例，并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。