Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Joint power allocation and orientation for uniform illuminance in indoor visible light communication

Open Access Open Access

Abstract

In this paper, we propose a technique called joint power allocation and orientation of light-emitting diodes (LEDs) to improve the illuminance quality of an indoor visible light communication (VLC) system. The optimization problem for LEDs, considering power allocation coefficients and normal vectors, is solved through an iterative algorithm. The simulation results with various LED geometries show that the proposed scheme can achieve a lower received power variance and a higher quality factor compared to previous power allocation schemes.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Recently, visible light communication (VLC) with the widespread deployment of light-emitting diodes (LEDs) has attracted a lot of attention around the world. The optical intensity is utilized to modulate the information bits, which can help alleviate the shortage of available radio frequencies. Moreover, VLC has various applications that cover a wide area and have a high demand value. Furthermore, multiple-input multiple-output (MIMO) techniques can improve the inherently limited modulation bandwidth of VLC via its spectral efficiency [1,2]. Consequently, multiple LEDs and photodiodes (PDs) have been widely employed in MIMO communication systems, thus being an attractive research topic in VLC [3,4].

On the other hand, a uniform illuminance condition is desirable in every VLC system [1]. Firstly, the inherent purpose of any indoor lighting system is to provide illuminance quality as evenly as possible while providing a dimming function with a wide range of brightness. Secondly, the uniformity of light signals around the room is also an important feature in indoor optical wireless communication scenarios since it helps ensure full coverage and improves system performance [5]. However, since it is difficult to get a line-of-sight (LOS) signal in many locations, the illuminance distribution and the uniform illumination condition is generally not taken into account in VLC research [6].

Recently, to achieve uniform irradiance for VLC systems, various techniques have been proposed [711]. One approach to realize a uniform illuminance condition is to specifically design the locations of the LEDs in the array. In [11], uniform illumination is achieved by using an ab initio design of an LED array. Moreover, in [8], a numerical optimization method was proposed based on a local search algorithm to design an LED array for a highly uniform illumination distribution around a room. However, specific configurations of the algorithm may lead to local extremum wherein the algorithm cannot obtain a global solution for the design problem, thus reducing the effectiveness of the proposed method. Another approach is the allocation of power for a specific LED geometry. In [7], the uniform irradiance on the receiver plane is achieved by using a simple heuristic power allocation (PA) scheme for a random LED array. Specifically, by modeling the LEDs locations using stochastic geometry with random parameters, the proposed algorithm helps to properly allocate the power for all LEDs to generate average uniform illumination. Recently, in [9], an optimal power allocation scheme was proposed. By considering the variance of the received power on the receiver plane as a metric and reforming the objective function as a convex optimization problem, uniform illuminance can be accomplished. Moreover, various LED array geometries have been considered such as square, circle-square, binomial point process (BPP) [7], and Matern type II hard-core point process (HCPP) [9] to achieve better lighting quality.

In this paper, we aim to further improve the uniform illuminance feature of a VLC system by introducing an additional degree of optimization. In particular, we form and effectively solve a convex optimization problem considering normal vectors and power allocation coefficients of LEDs. By altering the orientation of LEDs in conjunction with optimal power allocation, the variance of the received signals can be further decreased. Moreover, the proposed method can increase the quality factor of the system compared with previous scheme [9]. On the other hand, to facilitate the dimming control capacity in the system, the simulation result with various thresholds on the maximum total power of the system is also considered.

2. System model and the LEDs normal vectors

In Fig.  1, an indoor VLC system with $N$ LEDs is considered. A set of $K$ light received PDs is considered at the receiver plane. Moreover, all the parameters and properties of the system are well defined in [7,9,12]. Here, we only consider a VLC link with LOS characteristics. The channel gain between the $j$-th LED at the $i$-th location, can be simplified as [6]

$$h_{ij}=\frac{\left( \gamma +1 \right) A_pg}{2\pi d_{ij}^{2}}\cos ^{\gamma}\phi _{ij}\cos \theta _{ij},$$
where $\gamma$, $\phi _{ij}$, $\theta _{ij}$, and $d_{ij}$ denote the Lambertian order, the irradiance angle at the $j$-th LED with respect to the $i$-th PD, the incident angle at the PD with respect to the $j$-th LED, and the distance between the $j$-th LED and the $i$-th PD, respectively. Moreover in this paper, as an initial exploration, with the assumption of the general Lambertian coefficient [79,13] $\gamma =1$, the channel coefficient is simplified as
$$h_{ij}=\frac{A_pg}{\pi d_{ij}^{2}}\cos\phi _{ij}\cos \theta _{ij}.$$

In the other hand, we have

$$\cos \phi _{ij}=\frac{\overrightarrow{T}_j\cdot \overrightarrow{V}_{ij}}{\lVert \overrightarrow{T}_j \rVert \lVert \overrightarrow{V}_{ij} \rVert},$$
and
$$\cos \theta _{ij}=\frac{-\overrightarrow{V}_{ij}\cdot \overrightarrow{U_i}}{\lVert \overrightarrow{V}_{ij} \rVert \lVert \overrightarrow{U_i} \rVert},$$
where $\overrightarrow {T}_j$ represents the normal vector of the $j$-th LED in irradiance direction, $\overrightarrow {V}_{ij}$ denotes the vector from the $j$-th LED to the $i$-th PD, and $\overrightarrow {U}_i$ is the normal vector of the $i$-th PD in the incident direction, where generally $\overrightarrow {U}_i=\left [ \textrm {0 0 }1 \right ] ^T$. Moreover, Eqs. (3) and (4) can be represented as
$$\cos \phi _{ij}=\frac{\mathbf{t}_{j}^{T}\mathbf{v}_{ij}}{\lVert \mathbf{t}_j \rVert \lVert \mathbf{v}_{ij} \rVert}, $$
$$ \cos \theta _{ij}=\frac{-\mathbf{v}_{ij}^{T}\mathbf{u}_i}{\lVert \mathbf{v}_{ij} \rVert \lVert \mathbf{u}_i \rVert}. $$
where $\mathbf {t}_j$, $\mathbf {v}_{ij}$, and $\mathbf {u}_i$ are the position vectors of $\overrightarrow {T}_j$, $\overrightarrow {V}_{ij}$, and $\overrightarrow {U}_i$, respectively. It is assumed that the locations of the LEDs and PDs are fixed. The orientation of the $j$-th LED can be specified by a normal vector with parameter $\overrightarrow {T_j}\left ( \alpha _j,\beta _j \right )$ with $0\leqslant \alpha _j\leqslant \pi$ and $0\leqslant \beta _j\leqslant 2\pi$, as illustrated in Fig.  2, where $\beta _j$ and $\alpha _j$ denote the elevation angle and the azimuth angle, respectively. In practice, parameters $\beta _j$ and $\alpha _j$ can be used to adjust the orientation of the LEDs to improve the system illumination quality. Therefore, the normal vector of the $j$-th LED can be expressed as [14]
$$\mathbf{t}_j\equiv \overrightarrow{T_j}\left( \alpha _j,\beta _j \right) \triangleq \begin{cases} x_j=\lVert \overrightarrow{T_j} \rVert \sin \alpha _j\cos \beta _j\\ y_j=\lVert \overrightarrow{T_j} \rVert \sin \alpha _j\sin \beta _j\\ z_j=\lVert \overrightarrow{T_j} \rVert \cos \alpha _j\\ \end{cases},$$
In our paper, by optimizing the normal vectors $\mathbf {t}_j$ of the $j$-th LED, we can convert $\mathbf {t}_j$ back to spherical coordinates $\overrightarrow {T_j}\left ( \alpha _j,\beta _j \right )$.

3. Problem formulation

The channel coefficient with the help of (5) can be expressed as

$$h_{ij}=\frac{A_pg}{\pi d_{ij}^{2}}\frac{-\mathbf{v}_{ij}^{T}\mathbf{u}_i}{\lVert \mathbf{v}_{ij} \rVert \lVert \mathbf{u}_i \rVert}\frac{\mathbf{v}_{ij}^{T}\mathbf{t}_j}{\lVert \mathbf{v}_{ij} \rVert \lVert \mathbf{t}_j \rVert}.$$
Consequently, $h_{ij}$ can be expressed as a function of the LED normal vectors $\mathbf {t}_j$
$$h_{ij}=q_{ij}\left( \frac{\mathbf{v}_{ij}^{T}\mathbf{t}_j}{\lVert \mathbf{v}_{ij} \rVert \lVert \mathbf{t}_j \rVert} \right),$$
with the constants $q_{ij}=\frac {A_pg}{\pi d_{ij}^{2}}\frac {-\mathbf {v}_{ij}^{T}\mathbf {u}_i}{\lVert \mathbf {v}_{ij} \rVert \lVert \mathbf {u}_i \rVert }$. Moreover, by defining $\mathbf {\hat {v}}_{j}^{T}=\frac {\mathbf {v}_{j}^{T}}{\lVert \mathbf {v}_j \rVert }$ and $\mathbf {\hat {t}}_j=\left [ \frac {x_j}{\lVert \mathbf {t}_j \rVert }\,\,\frac {y_j}{\lVert \mathbf {t}_j \rVert }\,\,\frac {z_j}{\lVert \mathbf {t}_j \rVert } \right ] ^T$, (8) can be rewritten as
$$h_{ij}=q_{ij}\mathbf{\hat{v}}_{ij}^{T}\mathbf{\hat{t}}_j.$$
Therefore, from (9), the channel coefficient vector from all LEDs to location $i$-th where $\mathbf {h}_i=\left [ h_{i1}\cdots h_{iN} \right ] ^T$ is represented as
$$\mathbf{h}_i=\mathbf{q}_i\circ \mathbf{p}_i,$$
where $\circ$ is the element-wise product notation. In (10), $\mathbf {q}_i$ is the $N\times 1$ vector where each element $q_{ij}$ is the relation between the $i$-th PD and the $j$-th LED, and $\mathbf {p}_i$ is the $N\times 1$ vector where each element $p_{ij}$ is defined as $p_{ij}=\mathbf {\hat {v}}_{ij}^{T}\mathbf {\hat {t}}_j$. Furthermore, $\mathbf {p}_i$ can be calculated as
$$\mathbf{p}_i=\mathbf{\hat{V}}_i\mathbf{\hat{t}},$$
where $\mathbf {\hat {V}}_i=\textrm {blkdiag}\left ( \mathbf {\hat {v}}_{i1}^{T},\ldots ,\mathbf {\hat {v}}_{iN}^{T} \right ) =\left [ \begin {matrix}{} \mathbf {\hat {v}}_{i1}^{T}& \cdots & \mathbf {0}\\ \vdots & \ddots & \vdots \\ \mathbf {0}& \cdots & \mathbf {\hat {v}}_{iN}^{T}\\ \end {matrix} \right ]$ is an $N\times 3N$ matrix and $\mathbf {\hat {t}}$ is a $3N\times 1$ column vector that includes all the normal vectors for the LEDs $\mathbf {\hat {t}}=\left [ \mathbf {\hat {t}}_1;\mathbf {\hat {t}}_2;\cdots ;\mathbf {\hat {t}}_{N} \right ]$. In here, ’$\textrm {blkdiag}$’ is the block diagonalization matrix function. After inserting (11) into (10), we can derive
$$\mathbf{h}_i=\mathbf{q}_i\circ \mathbf{p}_i=\mathbf{q}_i\circ \left( \mathbf{\hat{V}}_i\mathbf{\hat{t}} \right) \\ =\mathbf{b}_i\mathbf{\hat{t}},$$
where the constant vector $\mathbf {b}_i=\left [ \begin {matrix} q_{i1}\mathbf {\hat {v}}_{i1}^{T}& \cdots & \mathbf {0}\\ \vdots & \ddots & \vdots \\ \mathbf {0}& \cdots & q_{iN}\mathbf {\hat {v}}_{iN}^{T}\\ \end {matrix} \right ]$.

 figure: Fig. 1.

Fig. 1. The system model.

Download Full Size | PDF

 figure: Fig. 2.

Fig. 2. The coordinates of an LED normal vector.

Download Full Size | PDF

 figure: Fig. 3.

Fig. 3. SNR profiles for different geometries.

Download Full Size | PDF

In the other hand, considering power allocation for the LEDs with the power allocation vector $\mathbf {w}=\left [ w_1\cdots w_N \right ] ^T$ of size $N\times 1$, the received power at the $i$-th PD can be expressed as

$$\mathbf{r}_i=\mathbf{w}\circ\mathbf{h}_i=\mathbf{w}\circ \left( \mathbf{b}_i\mathbf{\hat{t}} \right).$$
Therefore, total received power at the $i$-th PD as the sum of all elements of $\mathbf {r}_i$ can be rewritten as
$$R_i=\mathbf{w}^{\boldsymbol{T}}\mathbf{b}_i\mathbf{\hat{t}}.$$
Moreover, as in [9], the quality factor to evaluate the illuminance performance of the LED array can be expressed as
$$F_{\varLambda}=\frac{\bar{\varLambda}}{2\sqrt{\textrm{var}\left( \varLambda \right)}},$$
where $\varLambda _i=\frac {R_{i}^{2}}{\sigma _{i}^{2}}$ denotes the received electrical signal-to-noise ratio (SNR) of the $i$-th PD while $\bar {\Lambda }$ and $\textrm {var}\left ( \Lambda \right )$ denote the mean and variance of $\Lambda _i, i=\textrm {1},\ldots ,K$, respectively. Intuitively, to achieve uniform illuminance in VLC a system, the power received at any photo detector on the receiver plane should be as even as possible. Similar to [9], the variance of the received power should be minimized. Therefore, uniform illumination can be obtained by solving the following optimization problem
$$\textrm{P1: }\underset{\mathbf{w,\;t}}{\min\textrm{ }}\textrm{var}\left( R_i \right) $$
$$\textrm{s.t. }\textrm{sum}\left( \mathbf{w} \right) =P_{\textrm{total}}, $$
$$\mathbf{t}_{j}^{T}\mathbf{t}_j=1, $$
$$0\leqslant \mathbf{\hat{V}}_j\mathbf{t}\leqslant1, $$
$$ \mathbf{w}\geqslant 0, $$
where constraint (16b) requires the total power of all LEDs to equal $P_{\textrm {total}}$. Constraint (16c) with $\mathbf {t}=\left [ \mathbf {t}_1;\cdots ;\mathbf {t}_j;\cdots ;\mathbf {t}_N \right ]$ states that the lengths of all LED normal vectors are one, hence we can replace $\mathbf {\hat {t}}$ with $\mathbf {t}$ in the auxiliary constrain (16c). The constraint (16d) assures that $0\leqslant \cos \phi _{ij}\leqslant \pi /2$, or to put it another way, that there exists a LOS link between any PD and any LED. The last constraint indicates non negative power allocation coefficients.

4. Proposed solution

It can be seen that the constraints of problem (P1) are separable. Therefore, an iterative algorithm can be used to address problem (P1) with the block variables $\mathbf {t}$ and $\mathbf {w}$ [15]. Specifically, each time we optimize one variable while the other remains fixed [16]. This leads to the following two sub-problems in each iteration.

1 The subproblem w.r.t. $\mathbf {t}$

With fixed power allocation vector $\mathbf {w}$, the variable $\mathbf {t}$ is updated by solving the following optimization problem

$$\textrm{P2: }\underset{\mathbf{t}}{\min}\mathbf{t}^T\mathbf{\Pi t} $$
$$\textrm{s.t. }\mathbf{t}_{j}^{T}\mathbf{t}_j=1, $$
$$ 0\leqslant \mathbf{\hat{V}}_i\mathbf{t}\leqslant 1, $$
where the formation of (P2) can be referred to Appendix A. Moreover, by relaxing and linearizing constraint (17b) to $2\left ( \mathbf {t}_{j}^{\left ( k-1 \right )} \right ) ^T\mathbf {t}_j-\left ( \mathbf {t}_{j}^{\left ( k-1 \right )} \right ) ^T\mathbf {t}_{j}^{\left ( k-1 \right )}\geqslant 1$, the optimization problem (P2) can be solved by
$$\textrm{P3: }\underset{\mathbf{t}}{\min}\mathbf{t}^T\mathbf{\Pi t} $$
$$\textrm{s.t. }2\left( \mathbf{t}_{j}^{\left( k-1 \right)} \right) ^T\mathbf{t}_j-\left( \mathbf{t}_{j}^{\left( k-1 \right)} \right) ^T\mathbf{t}_{j}^{\left( k-1 \right)}\geqslant 1 $$
$$ 0\leqslant \mathbf{\hat{V}}_i\mathbf{t}\leqslant 1. $$

2 The subproblem w.r.t. $\mathbf {w}$

With fixed LED normal vectors $\mathbf {t}$, we can reform the object function in problem (P1) according to the variable $\mathbf {w}$. The variable $\mathbf {w}$ is updated by solving the following optimization problem

$$\textrm{P4: }\underset{\mathbf{w}}{\min}\mathbf{w}^{\boldsymbol{T}}\mathbf{\Delta w} $$
$$\textrm{s.t. }\textrm{sum}\left( \mathbf{w} \right) =P_{\textrm{total}}, $$
$$ \mathbf{w}\geqslant 0, $$
where the formation of (P4) can be referred to Appendix B. For clarification, the optimization algorithm is summarized in Algorithm 1, which is a two-stage iterative algorithm. The first stage is employed to solve problem (P3) to find the optimal LED normal vectors, while the second stage is employed to update the optimal LED power allocation coefficients. Each iteration of the first stage mainly consists of three steps. In the first step, where the input is the previous value of $\mathbf {w}$, the coefficient matrix $\mathbf {\Pi }$ is constructed and the optimization problem is formed and solved where the output is the vector $\mathbf {t}$. Similarly, in the second stage, the same steps are employed in solving the optimization problem (P4), where the input is the previous vector $\mathbf {t}$. The iteration is terminated when the convergence threshold is met.

oe-27-20-28575-i001

As mentioned in [9], since the variances of received power $\mathbf {t}^T\mathbf {\Pi t}$ and $\mathbf {w}^{\boldsymbol {T}}\mathbf {\Delta w}$ are always positive, $\mathbf {\Pi }$ and $\mathbf {\Delta }$ are positive-definite matrices. Moreover, since all the constraints are linear functions and hence, are convex. Thus the optimization problem (P3) and (P4) with quadratic objective functions are convex optimization problem [15]. Consequently, the optimization problem (P3) and (P4) can be efficiently solved using the CVX solver [17]. Using the solution obtained from previous iterations as the new center point, a new approximated problem can be solved. This procedure is convergent because the sequential iteration is guaranteed to reduce the objective value after each iteration. Moreover, the length of normal vectors in the relaxed constraint (18b) usually around one. Therefore, after each iteration, the normalization step would ensure $\lVert \mathbf {t}_{j}^{\left ( k \right )} \rVert =1$ without significant degradation in the quality of the solution. However, due to the relaxing and linearizing from constraint (17b) to (18b), the joint algorithm in general can only find the local solution of the problem. Therefore, the initialization of input values in the Algorithm 1 needs to be carefully considered. Moreover, through our simulation, we observe that a wise choice of initial vector $\mathbf {t}^{\left ( k \right )}$ can give a favorable solution and usually, with $\mathbf {t}_{j}^{\left ( k \right )}=\left [ \textrm {0} \quad 0 \quad -1 \right ] ^T$, the proposed method gives better performance in comparison with the PA in [9]. In the case of $\gamma \ne 1$, the simplified equation from (5) to (7) should be modified by employing the first order Taylor expansion of the channel coefficient $h_{ij}$ around $\mathbf {t}^{(k)}_j$ and expressing the received power variance in a similar manner. The modified optimization problem can be easily solved for specific $\gamma$ values. Therefore, in this paper, we consider the general case of $\gamma = 1$ which has been employed in various results [79,13].

5. Numerical results

In this section, the simulation results are provided for various LED geometries. The system is located inside a room with dimensions of $5\times 5\times 3$m. The receiver plane is 0.85m above the floor. The number of LEDs in each case is 16 with a semiangle of $60^\circ$. The total power of all LEDs is set to $P_{\textrm {total}}=1.1$W. Moreover, for fair comparison, a number of $K=40401$ PDs are uniformly distributed around the room as a square array of $201\times 201$ with each side is composed of 201 PDs. The distance between PDs is $0.025$m, the field of view of PD is $60^o$, and the physical area of the receiver PD is $1\textrm {cm}^2$. For comparison, the LED arrays of square, circle-square, BPP, and HCPP are considered with uniform power allocation and the power allocation scheme in [9]. In Fig.  3, the SNR profiles of different geometries for equal power allocation, the power allocation method in [9], and the proposed method are shown. It can be observed that the joint power allocation and LED orientation can achieve more uniform light since it provides an additional degree of optimization.

The advantage of the joint optimization solution can be seen in Table  1 and Table  2, where the variance of the received power and the quality factor of the schemes with various geometries are presented. The proposed scheme achieves excellent performance in terms of the variance of received power. Meanwhile, the quality factors of the proposed scheme are higher compared to the previous schemes. Moreover, as mentioned in [9], the circle-square geometry and HCPP yield the best quality factors, while the BPP geometry yields the lowest quality factor in both methods. On the other hand, the dimming capability is a very important feature of any VLC system. In the different conditions of lighting intensity required, the total power of the system $P_{\textrm {total}}$ is variously adjusted, and is generally many times higher than the upper threshold of individual LED power. In this section, we take into consideration the maximum value of individual LED power. By adding the constraint $\mathbf {w}\leqslant P_{\textrm {LED}}$ into (P4), we have the modified optimization problem

$$\textrm{P4*: }\underset{\mathbf{w}}{\min}\mathbf{w}^{\boldsymbol{T}}\mathbf{\Delta w} $$
$$\textrm{s.t. }\textrm{sum}\left( \mathbf{w} \right) =P_{\textrm{total}}, $$
$$ 0\leqslant \mathbf{w}\leqslant P_{\textrm{LED}}. $$

Tables Icon

Table 1. Variance comparison of the received power of the proposed scheme

Tables Icon

Table 2. Quality factor comparison of the proposed scheme

Consequently, the Algorithm 1 iteratively solves (P3) and (P4*) to obtain the solution for $\mathbf {w}$ and $\mathbf {t}$. By setting $P_{\textrm {LED}}=\textrm {0.125W}$, where the total power of the LEDs is $P_{\textrm {total}}=1.1$W, the SNR profile results are shown in Fig.  4. On the other hand, the variance of received power and the quality factor are shown in Table  3 and Table  4, respectively. It can be seen that the maximum power threshold has impact on both the power allocation scheme in [9] and the proposed scheme. This results in a significant increase in the variance and a decrease in the quality factor but the proposed scheme still has more advantage in general. Moreover, the joint optimization scheme is less sensitive to the maximum power threshold and consequently, results in a smaller change in the results. The restriction on the individual LED power also has greatly impact on the performance of HCPP geometry. Moreover, to clearly show the advantage of the proposed scheme when the dimming control is realized, in Fig.  5, we present the simulation results with various maximum total powers for the circle-square geometry. The maximum power of individual LEDs is set to $0.125$W, where the total power of the system $P_{\textrm {total}}$ is considered to be $0.03$W to $2$W since the maximum total power of all 16 LEDs is equal to $16\times 0.125$W or $2$W. Overall, the proposed scheme still provides lower variance and a higher quality factor with various values of $P_{\textrm {total}}$. Moreover, with the maximum value at less than $1.15$W, significant gains are observed in the quality factor for both the power allocation and the proposed method. On the other hand, a higher value of the maximum power threshold $P_{\textrm {total}}$ results in the convergence of the quality factor of the three schemes. Therefore, in many scenarios, the maximum total power of the system and the maximum power of each LED should be carefully considered to achieve optimal uniform illumination with high quality factor. Finally, we investigate the advantage of proposed scheme in term of achievable communication performance; more specifically, the instantaneous channel capacity. Unlike the conventional bit error rate evaluation, the capacity analysis is not limited to any specific channel coding or modulation schemes, and can be used to analyze the performance of VLC systems. With the power allocation vector $\mathbf {w}$, the instantaneous unconstrained channel capacity $C_i$ of the $i$-th PD can be expressed as [18]

$$C_i=\log _2\left[ 1+\frac{1}{\sigma ^2B}\left( \mathbf{w}\circ \mathbf{h}_i \right) ^{\boldsymbol{T}}\left( \mathbf{w}\circ \mathbf{h}_i \right) \right],$$
where $\sigma ^2$ is the noise variance and $B$ is the bandwidth of LED. In Fig.  6, the channel capacity performances of BPP and HCPP geometries are given for the case of equal power, power allocation [9], and proposed scheme. It can be observed that, with equal power level between LEDs, the channel capacity value is only acceptable in center of the room; while at locations near the corners, the value is only around 10 to 11 bps/Hz. On the other hand, with the power allocation and especially our proposed scheme, the channel capacity performances of all locations are over 13 bps/Hz, which are uniformly distributed across the room. From Fig.  6, the higher the received power from LEDs, the higher the channel capacity and the communication performance. Therefore, our proposed joint LED orientation and power allocation scheme can be a promising technique to provide highly uniform illumination with acceptable communication performance for indoor VLC systems.

 figure: Fig. 4.

Fig. 4. SNR profiles for different geometries with a threshold value of 0.125W.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. Performance of the circle square geometry with different maximum power values.

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. Channel capacity profiles for different geometries.

Download Full Size | PDF

Tables Icon

Table 3. Variance of received power of the proposed scheme with the threshold

Tables Icon

Table 4. Quality Factor of the proposed scheme with the threshold

6. Conclusion

In this paper, the joint optimization of LED power allocation and orientation is proposed for uniform illumination in VLC systems. By optimizing the power allocation coefficients with normal vectors for LEDs and solving the problem through an iterative algorithm, the optimal uniform illumination condition can be achieved. Moreover, the optimization with a maximum power threshold for individual LEDs is also considered. The simulation results show that the proposed scheme can achieve lower received power variance and higher quality factor values compared to previous schemes.

Appendix A. Formulation of optimization problem (P2)

We have

$$\mathbb{E}\left[ \left( P_i \right) ^2 \right] =\frac{1}{K}\mathbf{t}^T\left( \sum_{i=1}^K{\left( \mathbf{b}_i \right) ^T\mathbf{ww}^{\boldsymbol{T}}\mathbf{b}_i} \right) \mathbf{t},$$
$$\begin{aligned}\left( \mathbb{E}\left[ P_i \right] \right) ^2 &=\left( \mathbb{E}\left[ \mathbf{w}^{\boldsymbol{T}}\mathbf{b}_i\mathbf{t} \right] \right) ^2=\left[ \frac{1}{K}\mathbf{w}^{\boldsymbol{T}}\left( \sum_{i=1}^K{\mathbf{b}_i} \right) \mathbf{t} \right] ^2 \nonumber\\ &=\frac{1}{K^2}\mathbf{t}^T\left( \sum_{i=1}^K{\mathbf{b}_i} \right) ^T\mathbf{ww}^{\boldsymbol{T}}\left( \sum_{i=1}^K{\mathbf{b}_i} \right) \mathbf{t}.\end{aligned}$$
Therefore
$$\begin{aligned} \textrm{var}\left( P_i \right) & =\mathbb{E}\left[ \left( P_i \right) ^2 \right] -\left( \mathbb{E}\left[ P_i \right] \right) ^2 \\ & =\frac{1}{K}\mathbf{t}^T\left( \sum_{i=1}^K{\left( \mathbf{b}_i \right) ^T\mathbf{ww}^{\boldsymbol{T}}\mathbf{b}_i} \right) \mathbf{t}-\frac{1}{K^2}\mathbf{t}^T\left( \sum_{i=1}^K{\mathbf{b}_i} \right) ^T\mathbf{ww}^{\boldsymbol{T}}\left( \sum_{i=1}^K{\mathbf{b}_i} \right) \mathbf{t} \\ & =\mathbf{t}^T\left[ \frac{1}{K}\sum_{i=1}^K{\left( \mathbf{b}_i \right) ^T\mathbf{ww}^{\boldsymbol{T}}\mathbf{b}_i}-\frac{1}{K^2}\left( \sum_{i=1}^K{\mathbf{b}_i} \right) ^T\mathbf{ww}^{\boldsymbol{T}}\left( \sum_{i=1}^K{\mathbf{b}_i} \right) \right] \mathbf{t} \\ & =\mathbf{t}^T\mathbf{\Pi t}. \end{aligned}$$
with $\mathbf {\Pi }=\left [ \frac {1}{K}\sum _{i=1}^K{\left ( \mathbf {b}_i \right ) ^T\mathbf {ww}^{\boldsymbol {T}}\mathbf {b}_i}-\frac {1}{K^2}\left ( \sum _{i=1}^K{\mathbf {b}_i} \right ) ^T\mathbf {ww}^{\boldsymbol {T}}\left ( \sum _{i=1}^K{\mathbf {b}_i} \right ) \right ]$.

Appendix B. Formulation of optimization problem (P4)

We have

$$\mathbb{E}\left[ \left( P_i \right) ^2 \right] =\mathbb{E}\left[ \mathbf{w}^{\boldsymbol{T}}\mathbf{b}_i\mathbf{tt}^T\left( \mathbf{b}_i \right) ^T\mathbf{w} \right] =\frac{1}{K}\mathbf{w}^{\boldsymbol{T}}\left( \sum_{i=1}^K{\mathbf{b}_i\mathbf{tt}^T\left( \mathbf{b}_i \right) ^T} \right) \mathbf{w}.$$
Moreover, we have
$$\left( \mathbb{E}\left[ P_i \right] \right) ^2=\left( \mathbb{E}\left[ \mathbf{w}^{\boldsymbol{T}}\mathbf{b}_i\mathbf{t} \right] \right) ^2=\left[ \frac{1}{K}\mathbf{w}^{\boldsymbol{T}}\left( \sum_{i=1}^K{\mathbf{b}_i} \right) \mathbf{t} \right] ^2 \\ =\frac{1}{K^2}\mathbf{w}^{\boldsymbol{T}}\left( \sum_{i=1}^K{\mathbf{b}_i} \right) \mathbf{tt}^T\left( \sum_{i=1}^K{\mathbf{b}_i} \right) ^T\mathbf{w}.$$
Therefore
$$\begin{aligned} \textrm{var}\left( P_i \right) & =\mathbb{E}\left[ \left( P_i \right) ^2 \right] -\left( \mathbb{E}\left[ P_i \right] \right) ^2 \\ & =\frac{1}{K}\mathbf{w}^{\boldsymbol{T}}\left( \sum_{i=1}^K{\mathbf{b}_i\mathbf{tt}^T\left( \mathbf{b}_i \right) ^T} \right) \mathbf{w}-\frac{1}{K^2}\mathbf{w}^{\boldsymbol{T}}\left( \sum_{i=1}^K{\mathbf{b}_i} \right) \mathbf{tt}^T\left( \sum_{i=1}^K{\mathbf{b}_i} \right) ^T\mathbf{w} \\ & =\mathbf{w}^{\boldsymbol{T}}\left[ \frac{1}{K}\sum_{i=1}^K{\mathbf{b}_i\mathbf{tt}^T\left( \mathbf{b}_i \right) ^T-\frac{1}{K^2}\left( \sum_{i=1}^K{\mathbf{b}_i} \right) \mathbf{tt}^T\left( \sum_{i=1}^K{\mathbf{b}_i} \right) ^T} \right] \mathbf{w} \\ & =\mathbf{w}^{\boldsymbol{T}}\mathbf{\Delta w}. \end{aligned}$$
with $\mathbf {\Delta }=\frac {1}{K}\sum _{i=1}^K{\mathbf {b}_i\mathbf {tt}^T\left ( \mathbf {b}_i \right ) ^T-\frac {1}{K^2}\left ( \sum _{i=1}^K{\mathbf {b}_i} \right ) \mathbf {tt}^T\left ( \sum _{i=1}^K{\mathbf {b}_i} \right ) ^T}$.

Funding

National Research Foundation of Korea (NRF-2016R1D1A1B03934653, NRF-2019R1A2C1005920).

References

1. A. Jovicic, J. Li, and T. Richardson, “Visible light communication: opportunities, challenges and the path to market,” IEEE Commun. Mag. 51(12), 26–32 (2013). [CrossRef]  

2. T. Fath and H. Haas, “Performance comparison of MIMO techniques for optical wireless communications in indoor environments,” IEEE Trans. Commun. 61(2), 733–742 (2013). [CrossRef]  

3. P. H. Pathak, X. Feng, P. Hu, and P. Mohapatra, “Visible light communication, networking, and sensing: A survey, potential and challenges,” IEEE Commun. Surv. Tutorials 17(4), 2047–2077 (2015). [CrossRef]  

4. P. W. Berenguer, D. Schulz, J. Hilt, P. Hellwig, G. Kleinpeter, J. K. Fischer, and V. Jungnickel, “Optical Wireless MIMO Experiments in an Industrial Environment,” IEEE J. Select. Areas Commun. 36(1), 185–193 (2018). [CrossRef]  

5. N. Chi, LED-Based Visible Light Communications, Signals and Communication Technology (Springer-Verlag, Berlin Heidelberg, 2018).

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

7. G. V. S. S. P. Varma, R. Sushma, V. Sharma, A. Kumar, and G. V. V. Sharma, “Power allocation for uniform illumination with stochastic LED arrays; 00008.,” Opt. Express 25(8), 8659–8669 (2017). [CrossRef]  

8. A. Sewaiwar, S. V. Tiwari, and Y. H. Chung, “Smart LED allocation scheme for efficient multiuser visible light communication networks,” Opt. Express 23(10), 13015–13024 (2015). [CrossRef]  

9. G. V. S. S. P. Varma, “Optimum power allocation for uniform illuminance in indoor visible light communication,” Opt. Express 26(7), 8679–8689 (2018). [CrossRef]  

10. D. Ramane and A. Shaligram, “Optimization of multi-element LED source for uniform illumination of plane surface,” Opt. Express 19(S4), A639–A648 (2011). [CrossRef]  

11. S. Pal, “Optimization of LED array for uniform illumination over a target plane by evolutionary programming,” Appl. Opt. 54(27), 8221–8227 (2015). [CrossRef]  

12. Z. Wang, C. Yu, W.-D. Zhong, J. Chen, and W. Chen, “Performance of a novel LED lamp arrangement to reduce SNR fluctuation for multi-user visible light communication systems,” Opt. Express 20(4), 4564–4573 (2012). [CrossRef]  

13. Z. Wang, Q. Wang, S. Chen, and L. Hanzo, “An adaptive scaling and biasing scheme for OFDM-based visible light communication systems,” Opt. Express 22(10), 12707 (2014). [CrossRef]  

14. A. Nuwanpriya, S.-W. Ho, and C. S. Chen, “Indoor MIMO visible light communications: Novel angle diversity receivers for mobile users,” IEEE J. Select. Areas Commun. 33(9), 1780–1792 (2015). [CrossRef]  

15. S. P. Boyd and L. Vandenberghe, Convex optimization (Cambridge University, 2004).

16. Y. Huang, J. Zhang, and M. Xiao, “Constant Envelope Hybrid Precoding for Directional Millimeter-Wave Communications,” IEEE J. Select. Areas Commun. 36(4), 845–859 (2018). [CrossRef]  

17. “CVX: Matlab software for disciplined convex programming,” (2013).

18. S. Sugiura and H. Iizuka, “Element-by-element full-rank optical wireless MIMO systems using narrow-window angular filter designed based on one-dimensional photonic crystal,” J. Lightwave Technol. 34(24), 5601–5609 (2016). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1.
Fig. 1. The system model.
Fig. 2.
Fig. 2. The coordinates of an LED normal vector.
Fig. 3.
Fig. 3. SNR profiles for different geometries.
Fig. 4.
Fig. 4. SNR profiles for different geometries with a threshold value of 0.125W.
Fig. 5.
Fig. 5. Performance of the circle square geometry with different maximum power values.
Fig. 6.
Fig. 6. Channel capacity profiles for different geometries.

Tables (4)

Tables Icon

Table 1. Variance comparison of the received power of the proposed scheme

Tables Icon

Table 2. Quality factor comparison of the proposed scheme

Tables Icon

Table 3. Variance of received power of the proposed scheme with the threshold

Tables Icon

Table 4. Quality Factor of the proposed scheme with the threshold

Equations (40)

Equations on this page are rendered with MathJax. Learn more.

h i j = ( γ + 1 ) A p g 2 π d i j 2 cos γ ϕ i j cos θ i j ,
h i j = A p g π d i j 2 cos ϕ i j cos θ i j .
cos ϕ i j = T j V i j T j V i j ,
cos θ i j = V i j U i V i j U i ,
cos ϕ i j = t j T v i j t j v i j ,
cos θ i j = v i j T u i v i j u i .
t j T j ( α j , β j ) { x j = T j sin α j cos β j y j = T j sin α j sin β j z j = T j cos α j ,
h i j = A p g π d i j 2 v i j T u i v i j u i v i j T t j v i j t j .
h i j = q i j ( v i j T t j v i j t j ) ,
h i j = q i j v ^ i j T t ^ j .
h i = q i p i ,
p i = V ^ i t ^ ,
h i = q i p i = q i ( V ^ i t ^ ) = b i t ^ ,
r i = w h i = w ( b i t ^ ) .
R i = w T b i t ^ .
F Λ = Λ ¯ 2 var ( Λ ) ,
P1:  min   w , t var ( R i )
s.t.  sum ( w ) = P total ,
t j T t j = 1 ,
0 V ^ j t 1 ,
w 0 ,
P2:  min t t T Π t
s.t.  t j T t j = 1 ,
0 V ^ i t 1 ,
P3:  min t t T Π t
s.t.  2 ( t j ( k 1 ) ) T t j ( t j ( k 1 ) ) T t j ( k 1 ) 1
0 V ^ i t 1.
P4:  min w w T Δ w
s.t.  sum ( w ) = P total ,
w 0 ,
P4*:  min w w T Δ w
s.t.  sum ( w ) = P total ,
0 w P LED .
C i = log 2 [ 1 + 1 σ 2 B ( w h i ) T ( w h i ) ] ,
E [ ( P i ) 2 ] = 1 K t T ( i = 1 K ( b i ) T w w T b i ) t ,
( E [ P i ] ) 2 = ( E [ w T b i t ] ) 2 = [ 1 K w T ( i = 1 K b i ) t ] 2 = 1 K 2 t T ( i = 1 K b i ) T w w T ( i = 1 K b i ) t .
var ( P i ) = E [ ( P i ) 2 ] ( E [ P i ] ) 2 = 1 K t T ( i = 1 K ( b i ) T w w T b i ) t 1 K 2 t T ( i = 1 K b i ) T w w T ( i = 1 K b i ) t = t T [ 1 K i = 1 K ( b i ) T w w T b i 1 K 2 ( i = 1 K b i ) T w w T ( i = 1 K b i ) ] t = t T Π t .
E [ ( P i ) 2 ] = E [ w T b i t t T ( b i ) T w ] = 1 K w T ( i = 1 K b i t t T ( b i ) T ) w .
( E [ P i ] ) 2 = ( E [ w T b i t ] ) 2 = [ 1 K w T ( i = 1 K b i ) t ] 2 = 1 K 2 w T ( i = 1 K b i ) t t T ( i = 1 K b i ) T w .
var ( P i ) = E [ ( P i ) 2 ] ( E [ P i ] ) 2 = 1 K w T ( i = 1 K b i t t T ( b i ) T ) w 1 K 2 w T ( i = 1 K b i ) t t T ( i = 1 K b i ) T w = w T [ 1 K i = 1 K b i t t T ( b i ) T 1 K 2 ( i = 1 K b i ) t t T ( i = 1 K b i ) T ] w = w T Δ w .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.