## Abstract

This paper discusses the time-of-arrival (TOA) based indoor visible light communication (VLC) positioning system in a non-line-of-sight environment. The propagation delay is assumed to be gamma distributed. The generalized Cramer–Rao lower bound for multipath propagation is derived as the theoretical accuracy limitation. The performance of the positioning system is affected by the shape parameter and the scale parameter of gamma distribution. The influences on positioning accuracy of multipath effects are analyzed through discussing the physical meaning of the gamma distribution parameters. It is concluded that the lower bound of positioning accuracy is attained when variance of the non-line-of-sight propagation-induced path lengths is zero. The simulation result proves that the theoretical positioning accuracy is in the order of centimeters with the given scenario.

© 2015 Chinese Laser Press

## 1. INTRODUCTION

Recently, visible light communication (VLC) has been considered as the most attractive solution for indoor positioning because it has many features such as no electromagnetic interference and better positioning accuracy than radio signals. VLC technology could realize high-precision positioning of the intelligent mobile robots, which would enable them to engage in more sophisticated and complex work [1]. VLC positioning technology is also suitable for the blind navigation, which is aimed at supporting visually impaired people who travel indoors [2]. These scenarios require the high-precision positioning of submeters, which traditional radio signal positioning cannot afford. With all of these advantages, the VLC positioning technology is becoming a strong candidate for future indoor positioning applications.

The analysis of theoretical accuracy often provides a guide to practical applications, which could help people understand the performance limit and discover a proper method to approach it. But only few articles have analyzed the theoretical accuracy of the indoor VLC positioning system. The authors in [3] gave an indoor VLC positioning system based on the received signal strength indicator (RSSI) method and analyzed the impact of the LED parameters on positioning accuracy. It is concluded that, if the diffuse channel gain is measured previously in a certain environment and the modulation speed is far less than the channel cut-off frequency, the theoretical accuracy limit of VLC RSSI algorithm is not affected by multipath link. Reference [4] calculated the theoretical accuracy of time of arrival (TOA) algorithm and analyzed the influence of the signal parameters on distance estimates.

In the above articles, there are no studies of the impact of multipath effects on theoretical accuracy of a TOA-based indoor VLC positioning system in a non-line-of-sight (NLOS) environment, whereas the effect of NLOS propagation delays cannot be omitted in practical scenarios. Furthermore, in the condition that the statistical distribution of propagation delays in a NLOS environment is known, the range of theoretical accuracy limit has not been obtained. Thus, this paper has the important guiding significance in practical applications.

In this paper, the VLC optical channel models are set up at first. The TOA positioning system based on white lighting LEDs is described. The expression of a generalized Cramer–Rao lower bound (G-CRLB) based on a VLC TOA algorithm in a NLOS environment is described in detail. Then, we obtain the upper and lower bounds of the theoretical accuracy. The physical meaning of the gamma distribution parameters is analyzed. The simulation results show that the G-CRLB in the level of centimeters could be achieved with typical parameter values.

## 2. G-CRLB of TOA

In the VLC positioning system, we call the LEDs used for positioning base stations (BSs) and the object needed to locate is the mobile station (MS). We assume the optical filter with ideal bandpass, which causes no loss in the passband and completely blocks any optical frequencies in the stop band.

For ranging estimates based on intensity-modulated signals transmitted by LEDs, the received optical signal of a single path despite whether it is line-of sight (LOS) or NLOS is [4]

where $R$ is the responsivity of the photodiode, $\tau $ is the time taken for the light to travel from the transmitter to the receiver, $x(t)$ is the transmitted ranging optical signal, $n(t)$ is the noise signal, which is mainly composed by shot noise and can be modeled as additive white Gaussian noise (AWGN) with single-sided spectral density ${N}_{0}$. $\alpha $ is the attenuation of the optical channel. We assume that the optical axis of LEDs and receivers parallel to each other [3]: where $S$ is the detector effective area, $m$ is the Lambertian order, $d$ is the distance between the LED and the MS, and $h$ is the height of the room. We assume that the illumination distribution of LEDs conforms to approximate Lambertian model.In order to meet the condition of non-negativity and “smooth” for G-CRLB, $x(t)$ should be windowed as [4]

Let $\mathbb{Z}=\{1,2,\dots ,M,\dots ,Z\}$ be the set of indices of all BSs whose locations are known. Denote the set of indices of $M$ BSs that receive NLOS signals by $\mathbb{N}\ell =\{1,2,\dots ,M\}$. The complement of $\mathbb{N}\ell $, denoted as $\ell (=\mathbb{Z}\setminus \mathbb{N}\ell )$, is the set of indices of LOS stations, with its cardinality being $L=Z-M$. The estimated NLOS propagation-induced path lengths are ${({l}_{1},{l}_{2},\dots ,{l}_{M})}^{T}$. We define an $M$ dimensional vector $\theta $ as

The time-delay estimation can be approximated asLet ${\tau}_{z}$ be the time delay of the received signal at the $z$-th BS:

#### A. Case Where No Prior NLOS Statistics Are Available

The G-CRLB is independent of NLOS signals when no prior NLOS statistics are available, where we use CRLB to obtain the theoretical positioning accuracy [5]. The CRLB is defined as

where ${J}_{1}$ is the Fisher information matrix (FIM).And

Note that the window function $(1+\mathrm{cos}(2\pi t/{T}_{\mathrm{TOA}}))$ is a narrowband signal whose bandwidth is much lower than the frequency $f$. Thus, $\beta $ and $f$ have the relationship of [4]

#### B. Case Where Prior NLOS Statistics Are Available

In some cases, we may have the probability distribution of the prior NLOS statistics based on some statistical scattering models. Reference [6] points out that the gamma distribution is the best found for diffuse delay spreads, and varying the reflectivity of the walls did not affect the distribution fit tests.

Here, the amplifier gain in a receiver and its effect to signal power and noise power are not considered in this paper [4]. For simplicity, we assume that the NLOS propagation-induced path lengths are independent [5]. Then, we use gamma distribution to model the NLOS propagation-induced path lengths [6]:

In this case, we use the G-CRLB to obtain the theoretical positioning accuracy [5]. The G-CRLB is defined as

where ${J}_{2}$ is the FIM consisting of two components: Subscripts $D$ and $P$ stand for “data” and “prior” information.Noting that ${J}_{1}$ is actually independent of $l$, we reach

For $b=1$,

For the other cases of $0<b<1$ and $1<b<2$, we may not be able to find a closed-form expression for Eq. (18).

For $l>0$, $a>0$, $b>2$, we can obtain

Thus, another boundary of the theoretical accuracy for the MS is

Because ${J}_{P}\ge 0$, the bound of Eq. (13) is greater than that of Eq. (24). Thus, we can draw the conclusion that Eqs. (13) and (24) give the upper and lower bounds of the G-CRLB for the MS position estimate, respectively:

## 3. NUMERICAL RESULTS

In this section, the empirical parameters of VLC rms delay spread are first obtained. Next, we evaluate the influences of gamma distribution parameters and noise on G-CRLB by simulation.

In previous work, models for the rms delay spread of diffuse channels were developed [6]. We use the data of reference [6] for simulation. Figure 1 gives the result. As summarized, the gamma distribution is the best empirical models for diffuse delay spreads.

Figure 2 plots the G-CRLB versus the emitted optical power in rooms with different sizes. Here, only the lower bound is discussed. We select the parameters of ${T}_{\mathrm{TOA}}=0.01\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{s}$, $m=1$, $S=1\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{cm}}^{2}$, $R=0.4\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mA}/\mathrm{mW}$, and ${N}_{0}={\sigma}^{2}=9.2\times {10}^{-26}$ [4,7,8,9]. According to the experiments, the distribution fits the gamma model well when the data rate of optical wireless communication is below $100\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{Mb}/\mathrm{s}$. We choose direct modulation and $f=50\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{MHz}$. Here, the pre-equalization technique is considered to extend the LED modulation bandwidth.

Figure 2 shows that the G-CRLB increases with increasing $a$ and $b$. This is because the variance of ${l}_{z}$ decreases when $a$ increases, and the influence of $a$ is higher than that of $b$. We can obtain that the bound is better than 1 cm for an average transmitted optical power of 1 W in a $12\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}\times 12\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}\times 6\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}$ room.

The changes of the upper and lower bounds of the G-CRLB with noise are analyzed. The single-sided spectral density of the noise is [4]

where $q$ denotes the charge on an electron, $p$ is the background spectral irradiance, $\mathrm{\Delta}\lambda =360\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$ is the bandwidth of the optical filter in front of the photodiode, and $R$ and $S$ are the same as above. The optical power of LED is 1 W. $p$ varies from $5\times {10}^{-6}$ to $8\times {10}^{-6}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{\mu W}/{\mathrm{cm}}^{2}\xb7\mathrm{nm}$ [10].In Fig. 3, we can see that the influences of noise on the upper bound are significantly higher than on the lower bound. It is probably because we know more information about multipath information and the variance of NLOS delay has more influence than noise for the lower bound.

## 4. CONCLUSION

In this paper, we derive the G-CRLB model of the VLC positioning system based on a TOA algorithm in a NLOS environment. The G-CRLB range of the theoretical position estimation is obtained. For the lower bound, the variance of the NLOS delay is zero, which implies that the NLOS-induced path lengths are exactly obtained. In the upper bound, infinitely large variance means no information on NLOS delays at all. Models for the rms delay spread of diffuse channels are developed. We obtain the value of shape and scale parameters of gamma distribution in different rooms from experiment data and evaluate their influences on G-CRLB. It is concluded that the bound is in the order of centimeters for an average transmitted optical power of 1 W. Changes of the upper and lower bounds of G-CRLB with noise are also analyzed. We can see that the influences of noise on the upper bound are significantly higher than the lower bound.

In future works, further detailed research on the influence of incidence angles of light on the G-CRLB will be a priority. In addition, positioning algorithms for mobile tracking in a NLOS environment are planned for study.

## Acknowledgment

This work was supported by the National Key Basic Research Program of China (973 program) under grant 2013CB329204 and the Natural Science Foundation of China under grant 61375083.

## REFERENCES

**1. **R. Murai, T. Sakai, H. Kawano, Y. Matsukawa, Y. Kitano, Y. Honda, and K. C. Campbell, “A novel visible light communication system for enhanced control of autonomous delivery robots in a hospital,” in *IEEE/SICE International Symposium on System Integration (SII)* (2012), Vols. 510–516, pp. 16–18.

**2. **M. Nakajima and S. Haruyama, “Indoor navigation system for visually impaired people using visible light communication and compensated geomagnetic sensing,” in *1st IEEE International Conference on Communications in China (ICCC)* (2012), Vols. 524–529, pp. 15–17.

**3. **X. Zhang, J. Duan, Y. Fu, and A. Shi, “Theoretical accuracy analysis of indoor visible light communication positioning system based on received signal strength indicator,” J. Lightwave Technol. **32**, 3578–3584 (2014).

**4. **T. Q. Wang, Y. A. Sekercioglu, and A. Neild, “Position accuracy of time-of-arrival based ranging using visible light with application in indoor localization systems,” J. Lightwave Technol. **31**, 3302–3308 (2013). [CrossRef]

**5. **Y. Qi, “Wireless geolocation in a non-line-of-sight environment,” Ph.D. dissertation (Princeton University, 2003).

**6. **J. B. Carruthers and S. M. Carroll, “Statistical impulse response models for indoor optical wireless channels,” Int. J. Commun. Syst. **18**, 267–284 (2005). [CrossRef]

**7. **L. Zeng, D. C. O’Brien, H. L. Minh, G. E. Faulkner, K. Lee, D. Jung, Y. Oh, and E. T. Won, “High data rate multiple input multiple output (MIMO) optical wireless communications using white LED lighting,” IEEE J. Sel. Areas Commun. **27**, 1654–1662 (2009). [CrossRef]

**8. **J. Grubor, S. Randel, K. D. Langer, and J. W. Waleski, “Broadband information broadcasting using LED-based interior lighting,” J. Lightwave Technol. **26**, 3883–3892 (2008). [CrossRef]

**9. **T. Komine and M. Nakagawa, “Fundamental analysis for visible-light communication system using LED lights,” IEEE Trans. Consum. Electron. **50**, 100–107 (2004). [CrossRef]

**10. **Q. Jeffery, F. Greg, B. Patrick, and E. Eric, “Quantifying the indoor light environment-testing for light stability in retail & residential environments,” in *International Conference on Digital Printing Technologies NIP & Digital Fabrication Conference* (2004), pp. 689–698.