Training sequence based frequency-domain channel estimation for indoor diffuse wireless optical communications

Channel estimation is a key technology for wireless optical communication (WOC) systems. Based on the training sequence (TS), this article develops three frequency domain (FD) channel estimation approaches for diffuse wireless optical channels. Considering the propagation property of light and the characteristics of optical modulation, this article establishes a link model for the indoor WOC systems. Using the established link model, three FD channel estimation methods, such as the LS method, the minimum mean square error (MMSE) method and the ORL-MMSE method are proposed. The computational complexity analysis gives time complexities of the three channel estimation methods in FD. This article also evaluates the performance of the proposed three methods by computer simulation, measuring in terms of MSE and BER, respectively. Further more, based on the Cramer–Rao Bound theorem, optimal TSs are found and tabulated for different channel responses and TS lengths. By using the optimal TSs, channel estimation errors measured by BER performance are also investigated.


Introduction
Wireless optical communication (WOC) is an interesting technology combining the radio frequency communications with the traditional infrared communications [1].Using a light emitting diode (LED) for data transmission, WOC systems can offer several potential advantages over their radio frequency and infrared counterparts, such as security at the physical level, high-data-rate, cost-effective, wide-bandwidth [2].Therefore, there has been a growing interest in WOC systems in recent years.Despite these advantages, WOC systems for indoor applications are still having many critical issues to be resolved.Generally speaking, influenced by the complicated optical propagation environment, the propagation path between transmitter and receiver is not always unobstructed.Transmit lights were diffusely reflected by reflecting surface, such 710071, China Full list of author information is available at the end of the article as the ceiling, the wall or some other shiny surfaces before they reached to the receiver [3,4].Therefore, the indoor WOC system can be termed as diffusive propagation system.The diffusive WOC system image is drawn in Figure 1.Such diffuse propagation may cause the intersymbol-interference (ISI), which is a primary impediment to communication at high speeds [5,6].In order to eliminate the effects of ISI, the instantaneous channel impulse response must be obtained before the recovery of the transmitted information.Obviously, channel estimation plays a critical role in diffuse WOC systems.
To the best of the authors' investigation, a preamble training sequence (TS) is often used to estimate the channel impulse response in digital communication systems (as seen in Figure 2).For TS-based channel estimation in high speed wireless communications, there mainly are two different ways to estimate the channel impulses, i.e., time-domain (TD) method and frequency-domain (FD) method [7,8].As a typical channel estimation method, TD channel estimation has been discussed extensively for conventional radio communications [9,10].However, http://jwcn.eurasipjournals.com/content/2012/1/326LED LED only a few works have been done for the channel estimation for diffused WOC systems [11][12][13].In [11], TS-based channel estimation for WOC systems has been investigated, and the optimal TSs for a better channel estimation performance have been found by computer search in [12].Based on the "two-sided" channel estimation, Wang et al. [13] focuses on the complementary sets of sequences for estimating a channel impulse response in WOC systems.Since channel estimation methods in TD are time consuming and the search strategy for optimal TSs based on TD method incurs high computational complexity, a recently proposed channel estimation method based on FD techniques has received enormous attention in radio frequency communications [14,15].Different with the channel estimation in TD, FD channel estimation denotes the channel impulse response was estimated in the FD, followed by inverse discrete Fourier transform (DFT).To the best of the authors' knowledge, no similar results for channel estimation in WOC systems are available in the open literature.
Motivated by the research [7][8][9][10][11][12][13][14][15], this article focuses on the TS-based channel estimation methods in FD for WOC systems.Considering the propagation property of light and the characteristics of optical modulation, this article establishes the mathematical channel model for indoor WOC systems.Further, based on the established link model, three FD channel estimation methods, such as the LS method, the minimum mean square error (MMSE) method and the modified MMSE method which can be called optimal low-rank MMSE (OLR-MMSE) are proposed.In addition, this article compares the complexity of the three methods.Furthermore, the performance of the proposed three methods are studied and compared, measuring in terms of Mean square error (MSE) and bit error ratio (BER), respectively.Based on the Cramer-Rao Bound (CRB) theorem, optimal TSs are found and tabulated for different channel responses and TS lengths.By using the optimal TSs, channel estimation errors measured by BER performance are also investigated.
The rest of this article is organized as follows: the mathematical model for WOC systems is described in Section 2. Three channel estimation methods in FD for WOC channels using a know sequence are proposed in Section 3. Section 4 analyzes the computational complexity of the proposed FD channel estimation algorithms.The simulation results of performance of the proposed three channel estimation methods are provided in Section 5. Based on the CRB theorem, Section 6 investigates a valid performance metric for determining optimal TSs.In Section

Mathematical model in frequency domain (FD) for diffuse WOCs
Consider the IM/DD modulation for WOC systems, in which the desired waveform is modulated onto the instantaneous power of the carrier and a photo-detector produces a current proportional to the received instantaneous power.The binary information is modulated into "0" and "1".As illustrated in Figure 1, the received signal of the WOC systems is given by [1] where * denotes convolution, a(t) ∈ R + ∪ {0} denotes the instantaneous power in the transmitter, h(t) ∈ R + ∪ {0} is the impulse response of the channel, and n (t) ∈ R is the optical channel noise.Digital sampling the continuous signals in the Equation ( 1), the signal of the ith sampling instant with the sampling period T s can be expressed as Assuming the channel is represented as a finite impulse response filter with J taps.The channel impulse is given by where δ(t) is the Dirac delta function, h m is mth channel impulse response.Let the initial transmit sequence is {a m ∈ {0, 1} |m = 0, 1, . . ., N − 1 }.Therefore, the received samples described by (2) can be further expressed as where M = N + J − 1 and {n i } denotes the channel noise.Since, a DFT approach is used to estimate {h m }, {h m } and {a m } are augmented such that where j = √ −1.Similarly, the DFT of {h m }, y m and {n m } are {H k }, {Y k } and {N k }, respectively.Assuming the M-point DFT of the received signal samples can be denoted as Y = [Y 0 , Y 1 , . . ., Y M−1 ] T .Defining the input matrix Ā = diag({A k }), the system model in (4) can be described in FD as where is the Gaussian noise with variance σ 2 .

TS based FD channel estimation method
In data aided channel estimation, the known information to the receiver is inserted in information symbols so that the current channel state can be estimated.By analyzing the relationship between the known TS and the received symbols, the instantaneous channel impulse response can be estimated [16].

Least-squares (LS) estimation
LS channel estimation is commonly used in radio frequency communications when a TS is available.The optimality criterion of this method is to minimize the least square errors to find an optimal estimator to the unknown parameters [9].It is widely used in practice due to its salient feature that no probabilistic assumptions are made about the parameters.For the channel statistical model in FD as described in (6) for WOC systems, the sum of square error of channel estimation in FD can be defined as where • denotes the norm of a vector, (x) H is the Hermitian transposition of x.Therefore, the optimality criterion can be defined as By differentiating with respect to each channel coefficient H and setting the result to zero, i.e.
The closed-form expression for the LS channel estimate in FD can be obtained as Substituting Equation ( 6) into Equation ( 10), the LS channel estimate can be expressed as where H is the objective of the channel estimation, Ā−1 N denotes the estimated error.It can be observed that, the estimated channel impulse response is free from the interference between taps.In other words, LS channel estimation method is independent of the correlation properties of TS.However, the estimation error variance is largely depending on the channel noise.http://jwcn.eurasipjournals.com/content/2012/1/326

MMSE estimation
Unlike the LS approach, the optimality criterion of MMSE method is to minimize the MSEs to find an optimal estimator to the unknown parameters.Define the linear estimator that achieves the MMSE as ĤMMSE = VY (12) where V is the linear estimation operator [17].From the Equation ( 12), the MSE of the channel estimation for WOC systems can be described as where E [•] denotes the expectations operator.To obtain the expression of the FD channel estimation, it is expected that the estimation operator V must be determined at first.Obviously, the optimal operator V should be satisfying the function as follows By differentiating with respect to V and setting the result to zero, i.e.
∂ ∂V After solving the Equation ( 15), the close form of expression of operator V can be written as where R H = E HH H denotes the autocovariance matrix of channel coefficient in FD, R N = E NN H is the covariance matrix of the noise, (•) −1 denotes the inverse operation of a matrix.Therefore, the MMSE channel estimate in FD can be expressed as follows From the Equation ( 17), the MMSE approach can alleviate the effect of channel noise in some degree when compared with the LS approach.However, the MMSE channel estimation requires prior knowledge of the channel coefficient covariance matrix R H and noise covariance matrix R H , which means that the MMSE method is more complicated than the LS approach.

OLR-MMSE estimation
The MMSE estimator, as described in (17), is of considerable complexity since a matrix inversion is needed every time the data in Ā changes.We can reduce the complexity of this estimator by averaging over the transmitted data [15].
The first simplification of MMSE estimator is to replace the term Ā ĀH −1  in (17) with its expectation E Ā ĀH −1 .Assuming the transmitted symbols draw equal-probably.Then, we have where I is the identity matrix.Defining the average signal to noise ratio (SNR) as The estimator can be simplified by where To further reduce the complexity of the estimator, we proceed with the OLR approximations below.Rank reduction is achieved by using the singular value decomposition (SVD).The SVD of the channel autocovariance matrix is where U is a unitary matrix containing the singular vectors and is a diagonal matrix containing the singular values λ 1 ≥ λ 2 ≥ • • • ≥ λ M on its diagonal.The SVD also dramatically reduces the calculation complexity of matrices.
The system first determines the number of ranks required by the estimator, denoted by p, which should be no smaller than J + 1 [18].Then, given the noise variance and the channel autocovariance matrix R H , the receiver pre-calculates β, SNR, the unitary matrix U, and the singular values λ k s.It thus obtains the M × M diagonal matrix p with entries During the transmission, using the transmitted pilot Ā and Y, the ĤLS is calculated according to (10), and the OLR-MMSE estimator with rank p is given by The OLR-MMSE estimator can be interpreted as first projecting the LS estimates onto a subspace and then performing the estimation.Because the subspace has a small dimension (as small as J + 1) and still describes the channel well, the complexity of OLR-MMSE estimator is much lower than MMSE estimator with a good performance.http://jwcn.eurasipjournals.com/content/2012/1/326

Computational complexity analysis
Even when an algorithm is computationally feasible in principle, it may not be realizable in practice if the algorithm requires an inordinate amount of resources (such as running time, storage space, and communication bandwidth) to obtain the solution.Thus, the computational complexity of an algorithm is an important factor for the realization.Moreover, running time is of primary concern in most cases.Generally, the running time needed by an algorithm expressed as a function of the size of the input parameters is called the time complexity of the algorithm [19].In this section, the computational complexities of the proposed three FD channel estimation algorithms are analyzed in terms of time complexity.
To avoid the effect of running speed, the running time of an algorithm is usually measured by the number of operations needed by the algorithm to obtain the solution.However, in most cases, it is very difficult to determine the exact number of operations of an algorithm.Therefore, using the big-O notation, the time complexity of an algorithm is usually expressed in an asymptotic upper bound on the number of operations [20].
For simplicity, only major computational blocks of the proposed FD estimation algorithms are considered to evaluate the time complexity.Moreover, except for the operations related to the data processing, operations required for matrix formulation, loop counters, assignments of values to variables, etc., are not taken into account.In other words, the time complexity will be measured based on the numbers of multiplication and addition in the major computational blocks.Assuming that the length of TS is N, the time complexity is analyzed in Table 1.For a large-enough size of input parameters, the lower order terms of time complexity can be ignored.Therefore, from Table 1, we observe that the LS method is the most time efficient method among the FD algorithms.MMSE algorithm with a highest complexity and the OLR-MMSE is efficient than MMSE by reducing the rank of the channel correlation matrix.

Quality measure of FD channel estimate
To evaluate the performance of the proposed FD channel estimation schemes, the proposed three methods are studied and compared by computer simulation.Here, two comparisons are considered.Firstly, the proposed three methods are compared in terms of MSE, which denotes

FD-LS O(M)
O( 1) the accuracy of the channel estimation.Secondly, the three methods are compared, measuring by BER, which denotes the system performance influenced by channel estimation.
In the first comparison, one way of quantifying channel estimation performance is by MSE matrix [9], also referred to as the error covariance.The MSE matrix can be defined as where h and ĥ are the IDFT of H and Ĥ, respectively.Without lose of generality, we assume that the length of TS used in channel estimation N = 26 and the number of channel taps J = 3.
Figure 3 depicts the MSE performance versus the SNR with different FD channel estimation methods.It can be seen from Figure 3 that the MSE performance is decreases with the increase of the SNR.In addition, the best and worst MSE performances are achieved by the proposed MMSE channel estimation and the LS estimation algorithm, respectively.Also, for OLR-MMSE estimator, a degraded performance is shown due to the loss of channel information by reducing the rank of the channel correlation matrix.Moreover, since the LS method is sensitive to the signal noise power, the required SNR is 5 dB higher than the OLR-MMSE method with the MSE constraint and is set to be 10 −2 .
Figure 4 shows the MSE performance versus SNR with different estimation errors.With the LS channel estimation method, the MSE performances are free from the effects caused by the coefficient covariance matrix R H .However, the MSE increases significantly with the existence of the estimation error of R H when used MMSE channel estimation method.With the increase of the error variance of R H , the performance getting worse.Even though the SNR is high, the MSE performance for MMSE estimation is worse than LS estimation when there exits the estimation error for R H .
Figure 5 shows the second comparison.In Figure 5, the BER performance versus the SNR is depicted for the proposed three FD channel estimation schemes over the diffuse channel of WOC systems.It can be seen from Figure 5 that the BER performance is decreases with the increase of the SNR.The MMSE estimation method yields the best performance, and LS yields the worst.Also, for the OLR-MMSE estimator, a BER floor is shown due to loss of channel information by reducing the rank of the channel correlation matrix.

TS optimization
The CRB provides a lower bound on the statistical variance of any unbiased estimator, which means that the lower bound on CRB is the minimum variance for an estimator can obtain.The CRB for random channels has been wildly used as a performance measure for the design of optimal TSs [17,18].In this section, we first derive a closed-form expression for the CRB for the estimation of H.
Based on the CRB theorem, the Fisher information matrix (FIM) I (H) is given as follows [21] where x H is the Hermitian transposition of x, p (Y, H) is the joint probability density function of received signal sequences and channel coefficients, the derivatives are evaluated at the true value of H and the expectation is taken with respect to p (Y, H).
As described in the CRB theorem, it is sometimes possible to determine the minimum variance unbiased (MVU) estimator if the equality constraints of the CRB theorem are satisfied.Ĥ = g (Y) will be the MVU estimator if for some function g.Furthermore, the covariance matrix of Ĥ will be I −1 (H).For the linear channel model of ( 6), we have then the left part of the Equation (27) can be derived as where σ 2 is the variance of the Gaussian noise.Using the identities

Table 2 Complexity comparison
Optimal TS search method Multiplication Addition The simpler expression of Equation ( 28) can be expressed as Since, ĀH Ā is invertible, (30) can be rewritten as which is exactly in the form of (26) with Hence, the MVU estimator of H is given by (32a), and its covariance matrix is Further, the CRB of channel estimation in FD is given by  where tr( ) means the sum of the main diagonal elements of matrix .Obviously, (34) can be used to evalute the performance of different TSs.Therefore, the optimal TS can be obtained as Based on (35) optimal binary sequences are found by exhaustive computer search.The resulting sequences offer the best performance of error covariance at the output of the channel estimator.
Table 2 compares the time complexity of the search methods in FD and TD, respectively.For simplicity, only major computational block of the proposed FD search strategy and TD search strategy are considered to evaluate the time complexity.Comparing the search strategies described in [12] for TD channel estimation method with the FD method described in (35), we can see from the table that the optimal TSs search strategy in FD can significantly reduces the complexity of the search.In other words, the procedure of searching for optimum sequences in FD is time consuming than that of TD method.

Enumerative computer search for optimal TSs
As the placement of N clusters of pilot symbols is shown in Figure 2, a computer search has been performed to find the FD optimal TSs satisfying the Equation (35) for a given number of channel order J and TS length N.
Table 3 lists the optimal sequences found up to length 18.The sequences are binary and are represented in hexadecimal form for compactness.The most significant bit corresponds to the first symbol of the TS to be transmitted.The length of the preamble is uniquely determined by its location in the table.If the sequence length is not a multiple of four, as many zeros as is necessary are added or neglected to obtain the correct preamble length.To evaluate the performance, it is desired to compare our proposed optimal TSs in FD with the optimal TSs in TD for FD-LS channel estimation in terms of the BER via computer simulations, respectively.Although no TSs in TD have been proposed in published articles for WOC system, in our former research, we have investigated the optimal TSs for TD channel estimation for WOC systems [14].
Figure 6 compares the BER versus the SNR for the FD-LS channel estimation methods with proposed optimal TS and optimal TS in TD, respectively.The number of channel taps are J= 4 and J= 5.The TS length is N= 12.It can be seen from Figure 6 that the BER decreases with the increase of the SNR.With the same length of channel taps, the proposed FD optimal TS with the FD channel estimation can obtained a better BER performance compared with the optimal TS in TD method.In other words, the optimal TS in FD can get a better performance with the same channel state.It should be noted that the maximum value is obtained by the channel estimation method with a higher channel taps when the TS length hold fixed.
Figure 7 depicts the BER versus the SNR by using optimal TS in FD with the same number of channel taps but different number of TS symbols.The number of channel taps was J= 3, but the TS lengths are N = 6, 8, 10, respectively.It can be seen from Figure 7 that the BER decreases with the increase of the SNR.It should be noted that, the optimal TS with the less number of symbols can obtain the better BER performance under circumstances of the number of channel taps held fixed.

Conclusions
This article investigates the TS-based channel estimation in FD for indoor WOC systems.Firstly, the mathematical model is established for indoor WOC systems.Based on the established link model, three channel estimation methods, such as the LS method, the MMSE method and the ORL-MMSE method are proposed.Moreover, this article evaluates the performance of the proposed three methods by computer simulation, measuring in terms of MSE and BER, respectively.The mathematical analysis and the simulation results show that the LS estimator performs the worst of all methods and it has a lowest computational complexity, however.The MMSE estimator yields the best performance with its relatively higher complexity.The OLR-MMSE estimator gives the best trade off between performance and complexity.Based on the CRB theorem, a performance measure has been proposed to assess the quality of pilot symbols for FD channel estimation.Then, optimal TSs are found by a computer search.It has been demonstrated that the proposed method for searching the optimal TSs based on FD techniques invites a lower complexity search strategy compared with TD techniques.It should be noted that the optimal TSs in FD have a better BER performance compared to the TSs in TD.The minimum value is obtained by the channel estimation method with a lower channel taps when the TS length hold fixed.Moreover, with the number of channel taps held fixed, the longer the TS length is, the better the channel estimation performance can be obtained.

Figure 2
Figure 2The structure of data frame with TS.

Figure 3
Figure 3 MSE performance versus the SNR with different FD channel estimation methods.

Figure 4 Figure 5
Figure 4 MSE performance versus SNR with different estimation errors.

Figure 6 Figure 7
Figure6BER performance versus the SNR by using proposed optimal TS and optimal TS in TD, respectively.N = 12.