S S symmetry
Article Indoor Wavelet OFDM VLC-MIMO System: Performance Evaluation
Waleed K. Badawi , Marwa G. El-Hossary and Moustafa H. Aly *
Department of Electronics and Communication Engineering, Collage of Engineering and Technology, Arab Academy for Science, Technology and Maritime Transport, Alexandria 1029, Egypt; [email protected] (W.K.B.); [email protected] (M.G.E.-H.) * Correspondence: [email protected]
Abstract: Both light emitting diode (LED) characteristics for illumination and communication si- multaneously have made visible light communication-orthogonal frequency division multiplexing (VLC-OFDM) a strong competitive to radio frequency (RF). In this juncture, to improve signal to noise ratio (SNR) and coverage contour, the wavelet-OFDM is suggested for indoor VLC sys- tems. In this paper, a wavelet VLC-OFDM is proposed for imaging multiple-input multiple-output (MIMO) systems. The proposed wavelet-OFDM is exploited for a hybrid space-frequency domain pre-equalization technique instead of the traditional fast Fourier transform (FFT)-OFDM technique. The Meyer filter is selected and employed in the proposed technique. A comparable achievement is elaborated for several numbers of channels to achieve the enhanced performance in terms of bit rate and coverage contour. In addition, a useful comparison is executed between our wavelet VLC-OFDM and the traditional FFT-OFDM for a hybrid space-frequency domain pre-equalization technique. The simulation results emphasize the superiority point of wavelet VLC-OFDM MIMO system by −3 improving the coverage contour by ~20% over the traditional OFDM at a 10 bit error rate (BER) target. Hence, the proposed technique can be potentially executed in indoor VLC-MIMO systems.
Citation: Badawi, W.K.; El-Hossary, M.G.; Aly, M.H. Indoor Wavelet Keywords: visible light communication (VLC); wavelet-OFDM; orthogonal frequency division OFDM VLC-MIMO System: multiplexing (OFDM); imaging receiver (ImR); imaging angle diversity receiver (ImADR); pre- Performance Evaluation. Symmetry equalization; Meyer wavelet; Dmey wavelet 2021, 13, 270. https://doi.org/ 10.3390/sym13020270
Academic Editor: Palle E. 1. Introduction T. Jorgensen 1.1. State of the Art Regarding FFT-OFDM and DWT-OFDM Received: 17 January 2021 White light-emitting diodes (LEDs) are considered a successful technology for the next Accepted: 2 February 2021 Published: 5 February 2021 decade. Due to low energy consumption, high performance, and high durability, white LEDs can be used in a large variety of communication systems to increase performance
Publisher’s Note: MDPI stays neutral and minimize costs. Visible Light Communication (VLC) has the ability to carry data by with regard to jurisdictional claims in modulating light. This limits the extra cost of VLC techniques which use LEDs. Addi- published maps and institutional affil- tionally, these benefits enable people to access the Internet through the same visible light. iations. Consequently, the demand for VLC has increased rapidly with the increase of LED power. Major research was undertaken to establish high data indoor VLC systems [1]. VLC has an uncontrolled spectrum, precisely from 400 to 700 µm, providing a huge bandwidth of communication that leads to high data rates. The VLC spectrum in adjacent communication cells can be reused. Intensity modulation/direct detection (IM/DD) is Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. added to VLC systems, where the transmit messages can be modulated by light without This article is an open access article taking phase information into account. Besides, it can be directly sensed by the detector. distributed under the terms and The essential properties of light made VLC more effective, more power-efficient because the conditions of the Creative Commons LEDs are used in illumination and simultaneously are used in communication. Moreover, Attribution (CC BY) license (https:// the VLC is more secure, does not need license and has the ability to deliver high data rates creativecommons.org/licenses/by/ compared to RF communication. Hence, the proposed system has the above-mentioned 4.0/). properties and consequently has superiority over radio frequency systems.
Symmetry 2021, 13, 270. https://doi.org/10.3390/sym13020270 https://www.mdpi.com/journal/symmetry Symmetry 2021, 13, 270 2 of 19
Several VLC studies are involved in efficient modulation schemes like Orthogonal Frequency Division Multiplexing (OFDM) and Multiple Input Multiple Output (MIMO) to raise the data rate by maximizing the bandwidth. MIMO is used in many new technologies to increase data rate capacity and spectral efficiency [2–4]. MIMO systems contain multiple transmitter antennas and multiple receiver antennas. There are many types of MIMO techniques that are used for transmitting data over an indoor optical wireless channel to provide high link capacity and spectral efficiency, including spatial modulation (SM) and spatial multiplexing (SMP) [5]. The bandwidth limitation of LEDs is considered as one of the big challenges of achieving high efficiency in VLC systems [6]. Optical MIMO could be willing to offer a large data transmission rate by transferring data in the transponder and spatial spectrum [7]. Different pre-frequency domain equalization (Pre-FDE) mechanisms could be reported to improve LEDs bandwidth, such as analog Pre-FDE [6], digital Pre-FDE [8] in addition to the adaptive digital Pre-FDE [9]. VLC OFDM MIMO technique requires a detector for separating data from specific LEDs. Generally, there are two types of receivers with different concentrator designs that can be used in VLC-MIMO systems to increase VLC systems capacity [10]. These are categorized into two major optical categories; one of them is referred to imaging receiver (ImR) and the other is called a non-imaging receiver (NImR). The space division multiplexing VLC (SDM-VLC) technique is a promising technique for high speed indoor wireless communication. C. Chen et al. [11] suggested two protocols to boost an indoor SDM-VLC contour. There are two VLC link configurations: line-of-sight (LoS) which is mostly used in commercial systems and non-line-of-sight (Non-LoS) [12]. Thanks to OFDM for its high spectral efficiency, high capacity, and its ability to resolve multi-path channels. However, OFDM suffers from sensitivity to carrier frequency offset and synchronization errors. There is a lot of work that has been done to develop this, such as the short-time Fourier transform (STFT), the wavelet transform (WT), and wavelet packets (WP). The traditional OFDM techniques utilize inverse fast Fourier transform/fast Fourier transform (IFFT/FFT) systems at the transceiver for data multiplexing to simultaneously relay them through subcarrier numbers. The cyclic prefix (CP) is added before data transmission to minimize the inter-symbol-interference/inter-channel-interference (ISI/ICI) and to increase the bandwidth efficiency. Furthermore, CP reduces the channels of spectral inclusion. Wavelet (WT) fairly reflects a new term. It is used instead of FFT in OFDM. Wavelet- OFDM is named by orthogonal wavelet division multiplexing (OWDM), which is related to the inverse discrete wavelet transform (IDWT) rather than inverse discrete Fourier transform (IDFT) [13,14]. The great reason for utilizing DWT based on the OFDM system is the length of the basic functions which combat the narrowband interference better and inherently further resistant to ICI. R. Mishra et al. [15] studied the overview and performance estimation results of OFDM techniques including the traditional-OFDM and wavelet-OFDM. The authors in [16] proposed a new wavelet-OFDM analytical system which confirms Wi-Fi standard hiring the transmitter’s windowing function (rectangular waveform). Authors in [17] investigated a wavelet-OFDM for 4G of wireless communication. R. Asif et al. [18] elaborated the FFT-OFDM performance technique in opposite to WT based multicarrier system utilizing a simple zero-forcing (ZF) equalization in the time domain. Moreover, the wavelet has been applied in other different fields such as civil structures [19,20]. The Meyer WT is both orthogonal and symmetric [14]. This provides a numerical ap- proach to the DWT/IDWT through the use of a low pass filter (LPF)/high pass filter (HPF). To implement a wavelet transformation, these filters must satisfy the orthonormal basis, which means that they must be normalized and orthogonal to each other. Wael. H et al. utilized (16-QAM, 32-QAM, 64-QAM, and 128-QAM) and different types of wavelet filters Daubechies (Db-3, Db-5, Db-8, Db-10) and Haar filter to compare the performance between the FFT-OFDM and WT-OFDM [21]. Symmetry 2021, 13, 270 3 of 19
1.2. Frame Work This section shows, in brief, the novelty of our work and difference between it and the related previous work. Furthermore, it explains clearly our framework. For once, an imaging VLC-MIMO technique based on traditional FFT-OFDM has been proposed by C. Chen et al. [22] with a raw 1.2 Gbps bit rate for various channels. The technique used the dual property of the space/frequency domain pre-equalization scheme. Towards the aim of Signal to Noise Ratio (SNR) and Bit Error Rate (BER) improvement, a pre-space domain equalization (Pre-SDE) is performed after the Pre-FDE at a target BER of 10−3. In comparison with the system which uses only Pre-FDE, a communication coverage improvement by 52.6% is achieved in our system. In this paper, through VLC, we propose a wavelet-OFDM, instead of FFT-OFDM [22], for both the space and frequency domain pre-equalization technique for the first time. In the proposed system, due to the nature of WT overlapping symbols, the CP is eliminated. The removal of CP, which can be used to avoid multipath besides the Inter-Symbol In- terference (ISI), enables the WT to make full use of its spectral efficiency. This leads to improving system performance in terms of the coverage contour and data rate. In this eval- uation, we executed a Monte Carlo protocol to efficiently estimate an indoor VLC-MIMO performance characterized by N-channels. The obtained results reveal that the suggested system significantly enhances the SNR performance. Moreover, a performance comparison is performed between the suggested system based on the discrete Meyer (Dmey) wavelet (DWT-OFDM) and the FFT-OFDM for imaging VLC-MIMO systems for both the space and frequency domain pre-equalization technique. This comparison reflects the superiority of the DWT-OFDM proposed system over the FFT-OFDM one, in terms of the coverage contour for transmission and data rate. Thus, the suggested system can be a good applicant for light fidelity (Li-Fi) systems. The remainder of this paper is organized as follows. Section2 introduces the system model and numerical analysis of the proposed system involving all equations that are used to estimate technique performance. The simulation results are displayed and discussed in Section3. Section4 is devoted to the main conclusions.
2. System Model The distribution of the wavelet-OFDM MIMO system is illustrated within this section. This depends on the transmitter, the channel, in addition to the receiver which is utilizing the multiplexed method. The number of transmitting and receiving antennas is N and M, respectively.
2.1. The Transmitter We consider a wavelet-OFDM indoor imaging VLC-MIMO transmitter technique which utilizes a hybrid space-frequency domain pre-equalization system with a 4-LED modules transmitter (1, ... , N). Using LED is not sufficient to illuminate a room with these dimensions. However, in the proposed system, LED modules are used, and each module could consist of at least 60 distributed LEDs. The block diagram for the transmitter is illustrated in Figure1. Symmetry 2021, 13, x FOR PEER REVIEW 4 of 19
Symmetry 2021, 13, 270 4 of 19 Symmetry 2021, 13, x FOR PEER REVIEW 4 of 19
Figure 1. A wavelet-Orthogonal Frequency Division Multiplexing (OFDM) indoor imaging Visible Light Communication- Multiple Input Multiple Output (VLC-MIMO) transmitter utilizing the hybrid Pre Space Domain Equalization (Pre-SDE) system.
1. Firstly, the series input information is separated into N streams (according to N-trans- mitters), the positions of the LED modules are listed in Table 1. The data stream is converted from the data line to the data array, via a serial-to-parallel (S/P) converter.
2. The data is mapped to symbols by using quadrature amplitude modulation (QAM). FigureFigure 1. A wavelet-Orthogonal 1. A wavelet-Orthogonal Frequency FrequencyThe Division transmitter Division Multiplexing Multiplexing utilizes (OFD a 64 (OFDM)M) QAM indoor indoordigital imaging imagingmo Visibledulation Visible Light to Light Communication-map Communication- the serial bits Multiple into the Input Multiple Output (VLC-MIMO) transmitter utilizing the hybrid Pre Space Domain Equalization (Pre-SDE) system. Multiple Input Multiple Output (VLC-MIMO)OFDM symbols transmitter Χ (i) utilizing as N is the the LED hybrid modules Pre Space number Domain and Equalization 0≤i≤s−1, (Pre- while s is SDE) system. the number of the parallel data stream. 1. Firstly, the series input information is separated into N streams (according to N-trans- 3. Then, Hermitian symmetry (HS) is enforced to produce real values data. mitters), the positions of the LED modules are listed in Table 1. The data stream is 1. 4. Firstly,The Pre-FDE the series is inputperformed information by employing is separated the electrical-optical-electrical into N streams (according (EOE) to N- fre- converted from the data line to the data array, via a serial-to-parallel (S/P) converter. transmitters),quency response. the positions The LED of themodules LED moduleslocations are are listed listed in in Table Table1. 1 The [22] data and stream the meas- 2. The data is mapped to symbols by using quadrature amplitude modulation (QAM). is convertedured EOE fromfrequency the data response line to theof a data commercially array, via a available serial-to-parallel LED with (S/P) 50-MHz converter. modu- The transmitter utilizes a 64 QAM digital modulation to map the serial bits into the 2. Thelation data BW is mapped is illustrated to symbols in Ref. by [23]. using quadrature amplitude modulation (QAM). Χ (i) 0≤i≤s−1, OFDMThe transmitter symbols utilizes as aN 64 is QAM the LED digital modules modulation number to and map the serial bits while into thes is the number of the parallel data stream. TableOFDM 1. Light symbols EmittingX (Diodei) asN (LED) is the modules LED modules locations number (m) [22].and 0 ≤ i ≤ s − 1, while s is N 3. Then,the number Hermitian of the symmetry parallel data(HS) stream.is enforced to produce real values data. N = 1 (2.5, 2.5, 3) 3.4. TheThen, Pre-FDE Hermitian is performed symmetry by (HS) employing is enforced the to electrical-optical-electrical produce real values data. (EOE) fre- 4. NquencyThe = 2 Pre-FDE response. (1.5, is2.5, performed The 3) (3.5, LED 2.5, modules by 3) employing locations the are electrical-optical-electrical listed in Table 1 [22] and (EOE)the meas- fre- Nuredquency = 3 EOE response. (1.5,frequency 1.5, 3) The response(3.5, LED 1.5, 3) modulesof (2.5, a commercially 3.5, locations3) available are listed LED in Tablewith 150-MHz[ 22] and modu- the Nlationmeasured = 4 BW (1.5,is EOE illustrated 1.5, frequency 3) (3.5, in 1.5,Ref. response 3) [23]. (1.5, 3.5, of a 3) commercially (3.5, 3.5, 3) available LED with 50-MHz modulation BW is illustrated in Ref. [23]. 5.Table5. Wavelet 1.Wavelet Light Emitting transform transform Diode for for(LED) modulation modulation modules is locations performed.is perfor (m)med. [22] As. As shown shown in Figurein Figure1, the 1, blockthe block N = 1named named (2.5, IDWT IDWT 2.5, 3) is ais functiona function in thein the MATLAB MATLAB toolbox toolbox which which utilizes utilizes Low Low Pass Pass Filter Filter N = 2(LPF) (LPF) (1.5, and and High2.5, High 3) Pass(3.5, Pass Filter2.5, Filter 3) (HPF). (HPF). WT WT utilizes utilizes filters filters like like a vector a vector of approximation of approximation coefficientcoefficient (CA) (CA) and and a a vector vector of of detail detail coefficientscoefficients (CD). The signal signal is is divided divided into into sub- N = 3 (1.5, 1.5, 3) (3.5, 1.5, 3) (2.5, 3.5, 3) bands which are divided into low and high frequencies. OFDM symbol 〖Χ〗_N (i) N = 4sub-bands (1.5, 1.5, which 3) (3.5, are divided1.5, 3) (1.5, into 3.5, low 3) and(3.5, high 3.5, 3) frequencies. OFDM symbol _N (i)is is converted converted from from parallel parallel to to serial. serial. It It has has a a vector vector YY YY that that may may be be transposed transposed into into CA CAas as represented represented in in Figure Figure 22 which which is is called called approximated approximated coefficients. coefficients. 5. Wavelet transform for modulation is performed. As shown in Figure 1, the block named IDWT is a function in the MATLAB toolbox which utilizes Low Pass Filter Table 1. Light Emitting Diode (LED) modules locations (m) [22]. (LPF) and High Pass Filter (HPF). WT utilizes filters like a vector of approximation N =coefficient 1 (CA) and a vector (2.5, of 2.5, detail 3) coefficients (CD). The signal is divided into sub- 〖 〗 N =bands 2 which are divided (1.5,into 2.5, low 3) and (3.5, 2.5,high 3) frequencies. OFDM symbol Χ _N (i) is converted from parallel to serial. It has a vector YY that may be transposed into CA N = 3 (1.5, 1.5, 3) (3.5, 1.5, 3) (2.5, 3.5, 3) as represented in Figure 2 which is called approximated coefficients. N = 4 (1.5, 1.5, 3) (3.5, 1.5, 3) (1.5, 3.5, 3) (3.5, 3.5, 3)
SymmetrySymmetry2021 2021, 13, 13, 270, x FOR PEER REVIEW 55 of of 19 19
FigureFigure 2. 2.Proposed Proposed transmitter transmitter model model of of Discrete Discrete Wavelet Wavelet Transform Transform (DWT)-Orthogonal (DWT)-Orthogonal Frequency Fre- Divisionquency Division Multiplexing Multiplexing (OFDM). (OFDM).
TheThe signals signals are are up-sampled up-sampled by by LPF LPF coefficients coefficients while while the the HPF HPF filter filter can can be be performed performed byby the the convolution convolution with with Zero Zero Padding Padding (ZP) (ZP) data data Detail Detail Coefficient Coefficient (CD). (CD). This This is done is done to get to high-frequencyget high-frequency data. data. The LPFThe LPF filter filter contains contains approximation approximation coefficients. coefficients. Data Data is simulated is simu- bylated a MATLAB by a MATLAB toolbox. toolbox. The transmitted The transmitted signal si mustgnal be must in the be discretein the discrete domain domain x[k]. x[k]. VariousVarious wavelet wavelet families families have have several several filters filters containing containing different different values values of of CA CA as as well well asas CD.CD. AllAll ofof themthem seekseek toto satisfy satisfy the the WT WT bases. bases. TheThe CPCP isis notnot utilizedutilized forfor WT-OFDMWT-OFDM whichwhich isis due due to to the the wavelet wavelet orthogonality orthogonality overlapping overlapping attribute attribute and and its its better better stopband stopband attenuation.attenuation. UnlikeUnlike FFT-OFDM,FFT-OFDM, CP CP is is important important to to overcome overcome the the ISI. ISI. Each Each one one of of them them is is builtbuilt toto form form the the sum sum of Mof waveformsΜ waveformsφ(κ )ϕ,(whereκ), whereφ(κ )ϕis(κ the) is complex the complex exponential exponential basic functionsbasic functions and Φm andis the Φ scaling is the function. scaling function. The transmitted The tran signalsmitted may signal be shown may inbe ashown discrete in form.a discrete The transmittedform. The symboltransmitted is built symbol by performing is built by inverse performing DWT (IDWT)inverse [DWT14,24, 25(IDWT)]: [14,24,25]: [ ] = M−1 [ − ] x k ∑τ∑m =0 ajΦm k sM . (1) 𝑥[k] = ∑∑ a Φ [k − sM]. (1) ForFor performing performing the the wavelet-OFDM wavelet-OFDM technique technique mentioned mentioned in (1), in we (1), execute we execute the Meyer the wavelet.Meyer wavelet. The Meyer The wavelet Meyer is wavelet a limited is frequency a limited orthogonalfrequency wavelet.orthogonal Meyer wavelet. wavelet Meyer can bewavelet implemented can be implemented in fast wavelet in transform fast wavelet (FWT) transform and DWT (FWT) because and ofDWT its symmetry because of and its orthogonality, in addition to its fast ability of localization and decaying from the central symmetry and orthogonality, in addition to its fast ability of localization and decaying peak than any inverse polynomial [14]. from the central peak than any inverse polynomial [14]. OWDM is a program which is very flexible and simple. It has lower complexity with OWDM is a program which is very flexible and simple. It has lower complexity with a low filter than FFT processors. Furthermore, the filter type selected must be dependent a low filter than FFT processors. Furthermore, the filter type selected must be dependent on the channel situation or the information source. A Dmey is a discrete form of the Meyer on the channel situation or the information source. A Dmey is a discrete form of the Meyer wavelet. The Meyer wavelet and scaling function is defined in the frequency domain [14]. wavelet. The Meyer wavelet and scaling function is defined in the frequency domain [14]. 1.1. AA parallel-to-serialparallel-to-serial (P/S)(P/S) operation operation is performed to obtain the time domain signal. 2.2. Pre-SDEPre-SDE isis implementedimplemented byby normalizingnormalizing allall datadata streamsstreams totoensure ensure that that the the detected detected N-data have the same value of SNR. The Pre-SDE concept of an optical DCO-OFDM N-data have the same value of SNR. The Pre-SDE concept of an optical DCO-OFDM system is based on the VLC-MIMO N-channel imaging technique. The output optical system is based on the VLC-MIMO N-channel imaging technique. The output optical power and the modulation index for every LED are Popt and ξ, respectively. Pre-SDE power and the modulation index for every LED are Popt and ξ, respectively. Pre-SDE is performed by using power allocation matrix A = diag (a1; a2; ······ ; aN) which is performed by using power allocation matrix A = diag (a1; a2; · · ·; aN) which adjusts adjusts all electrical powers modulating data. all electrical powers modulating data.
Symmetry 2021, 13, 270 6 of 19
Furthermore, the transmitted optical vector S could be calculated by Equation (2), where x is the modulated data vector [22]:
S = Popt (1 + ξ A x). (2)
The modulation index (ξ) is supposed to be the same for all LEDs. Because of A, the electrical powers of different modulating data may be diverse. To preserve a fixed average electrical power as soon as exciting Pre-SDE, power restrictions are required to the diagonal elements of matrix [22]. The signal is normalized and a digital-to-analog (D/A) conversion has been performed. At every stream, a DC bias is applied to get a unipolar signal. Hence, the complex-valued OFDM signal is given by: x(t) = xRe + j xIm, (3) T where x = (x1 x2, ······ xN) is a vector of modulating signals and xi (i = 1, 2, ······ , N) is 2 the normalized bipolar OFDM signal with variance σ xi = 1. The electrical powers of different modulating signals may be different. To maintain a constant average electrical power when performing Pre-SDE, a power constraint is imposed on the diagonal elements of matrix A. These elements, ai, satisfy the following equation [22]: 1 N a2 = 1. (4) N ∑i=1 i
2.2. Channel Estimation 2.2.1. White Gaussian Noise Estimation
The equalizer output includes a Gaussian noise with a total variance GN, which is the sum of shot noise and thermal noise besides the ICI contributions using the optical path difference. We disregard the noise contributions from gate leakage current and 1/f noise:
2 2 2 2 GN = σshot + σthermal + γ σrISI. (5)
The received power PrISI by ICI is [26]:
Z ∞ LEDs PrISI = ∑ = hi(t) ⊗ X(t) dt. (6) T i 1 The shot noise variance is given by [26]:
2 σshot = 2qγ PrSignal + PrISI B + 2qIbgI2B, (7)
where q is the electronic charge, B is the equivalent noise bandwidth, Ibg is background current. The noise bandwidth factors I2 = 0.562. The thermal noise variance is calculated by [27]:
2 2 8πkTk 2 16π kTkΓ 2 2 3 σthermal = ηAI2B + η A I3B . (8) G gm Furthermore, both sides clarified noise from feedback-resistor, and noise from FET channels, respectively. Similarly, k is Boltzmann constant, TK is the absolute temperature, G is the open-loop voltage gain, η is the constant photodetector capacitance (PD) per unit area, Γ is the FET channel noise factor, gm is the FET trans-conductance, I3 = 0.0868, 2 Tk = 295 K, Γ is the FET channel noise, G = 10 dB, gm = 30 mS, η = 1.5, Γ = 112 pF/cm , even B = 100 Mbps. In our work, the surrounding current comes from the direct sum light [27].
2.2.2. Non-Imaging LOS Channel Gain A generalized Lambertian radiation pattern can model the LOS irradiance of LED as shown in Figure3. The vector υRS can be written as υRS = [a, b, c] = [XR,YR,ZR]–[XS, Symmetry 2021, 13, x FOR PEER REVIEW 7 of 19
ZS], where [XS, YS, ZS] and [XR, YR, ZR] are the positions of transmitter and receiver, respec- tively [28]. The optical LOS channel gain between the tth LED lamp and the rth PD is Symmetry 2021, 13, 270 7 of 19 proved by:
( ) μη cos (∅rt) cos(θrt) , 0 ≤ θ ≤ϕ h (0) = . (9) YS,ZS], where [XS,YS,ZS0,] and [X R ,Y R ,Z R ] are the positions of transmitter θ >ϕ and receiver, respectively [28]. The optical LOS channel gain between the tth LED lamp and the rth PD To conclude, m = −ln2/ln (cos Ψ1/2) is the emitted lighting order with Ψ1/2 is the semi- is proved by: angle of the transmitter at half power, APD is the PD active area, drt is the distance between
the tth LED lamp and( the(m +rth1)A PD,PD μ andm η are the optical filter gain likewise the optical 2 µη cos (∅rt) cos(θrt), 0 ≤ θ ≤ φ lens, respectively,hrt(0) =ϕrt is the2 πemissiondrt angle, and θrt is the incident angle. Moreover,. if θrt is (9) θ > φ outside the Field of View 0,(FOV) of the receiver, the optical gain hrt is zero and when an imaging non-diversity receiver (ImADR) is used, θrt becomes zero.
FigureFigure 3. 3.LineLine of ofSight Sight (LOS) (LOS) channel channel configuration configuration [28] [28. ].
2.2.3. Non-ImagingTo conclude, LOS m = Channel−ln2/ln Matrix (cos Ψ 1/2) is the emitted lighting order with Ψ1/2 is the semi-angleThe general of the 4 × transmitter 4 MIMO channel at half matrix power, H A isPD drivenis the by: PD active area, drt is the distance between the tth LED lamp and the rth PD, µ and η are the optical filter gain likewise the h ⋯h optical lens, respectively, φrt is the emission angle, and θrt is the incident angle. Moreover, H = ⋮⋱⋮ . (10) if θrt is outside the Field of View (FOV) of the receiver, the optical gain hrt is zero and when h ⋯h an imaging non-diversity receiver (ImADR) is used, θrt becomes zero. Furthermore, the channel gain h is progressed, where r and t refer to the numbers of2.2.3. transmitters Non-Imaging and receivers, LOS Channel respectively. Matrix Similarly,The general the 4ICI× is4 MIMOdeserted channel by ImR matrix or ImADR H is driven use, the by: matrix of the channel H of N-channel imaging VLC-MIMO technique is evidenced in [26]: h11 ··· h14 = h ,h ,…….,h H diag ( . . . ). (11) H = . .. . . (10) h41 ··· h44 2.3. The Receiver WeFurthermore, consider a VLC-MIMO the channel indoor gain hrt imagingis progressed, receiver where system r and based t refer on towavelet-OFDM the numbers of usingtransmitters one receiver and receivers,(M) of the respectively. hybrid space-frequency domain pre-equalization system. The receiverSimilarly, location the ICI (scenario-room is deserted by corner) ImR or is ImADR given by use, (0, the0, 0.85 matrix m). ofThe the system channel block H of diagramN-channel is shown imaging in Figure VLC-MIMO 4. technique is evidenced in [26]: • The produced N serials of unipolar data have been used to drive N LED modules, H = diag (h11, h22, ...... , hNN). (11) respectively, as shown in Figure 1. The LED radiation light is sensed by an ImR or an 2.3.ImADR The Receiver at the detector side. The recipient side of the wavelet-OFDM VLC-MIMO techniqueWe consider is illustrated a VLC-MIMO in Figure indoor 4. imaging receiver system based on wavelet-OFDM • usingAn oneA/D receiverconversion (M) is of utilized the hybrid to convert space-frequency the analog signals domain to pre-equalizationdigital. The 3-dB system.mod- Theulation receiver bandwidth location (scenario-roomof LED is adjusted corner) as 50 is MHz. given by (0, 0, 0.85 m). The system block diagram is shown in Figure4.
Symmetry 2021, 13, x FOR PEER REVIEW 8 of 19
Symmetry 2021, 13, 270 The specific processes to perform digital Pre-FDE were fully illustrated in [22]. In8 ofthe 19 OFDM transmitter (Tx), a Pre-FDE has performed adaptively according to the estimated SNR information.
FigureFigure 4. 4. AA Wavelet-OFDM Wavelet-OFDM indoor indoor imaging VLC-MIMO rece receiveriver by utilizing the hybrid Pre-SDE system.
2.3.1. Received Wavelet-OFDM Data • The produced N serials of unipolar data have been used to drive N LED modules, DWTrespectively, receiver as is shownthe reverse in Figure operation1. The of LED IDWT. radiation It performs light a is DWT sensed according by an ImR to the or MATLABan ImADR toolbox at. theThe detector transmitted side. signal The recipient ∪[k] is sidethe front-end of the wavelet-OFDM receiver successive VLC-MIMO digi- tally techniquemodulated is symbols illustrated to inbe Figure transformed,4. the data is decomposed across filters, LPF and• HPFAn A/Drelated conversion to CA in addition is utilized to the to CD. convert CA is the the analog output signals signal of to the digital. approximation The 3-dB coefficientsmodulation or LPF bandwidth and CD are of the LED output is adjusted data of as the 50 MHz.detailed coefficients or HPF to exe- cute thatThe specificoperation. processes The data to performis transposed digital before Pre-FDE converting were fully from illustrated S/P and in the [22 wavelet]. In the familyOFDM used transmitter is Demy (Tx), wavelet. a Pre-FDE has performed adaptively according to the estimated [ ] SNR∪ information.k contains some zeros elements which are decomposed from CD information produced in the transmitter comparing to CA data. So, the CD data will be considered as imaginary2.3.1. Received data Wavelet-OFDMand the CA will Databe considered as real data to avoid any loss of data. That is all DWTgo viareceiver the DWT-OFDM is the reverse receiver operation which of can IDWT. be seen It performs in Figure a 5. DWT according to the MATLAB toolbox. The transmitted signal ∪ [k] is the front-end receiver successive digitally modulated symbols to be transformed, the data is decomposed across filters, LPF and HPF related to CA in addition to the CD. CA is the output signal of the approximation coefficients or LPF and CD are the output data of the detailed coefficients or HPF to execute that operation. The data is transposed before converting from S/P and the wavelet family used is Demy wavelet. ∪[k] contains some zeros elements which are decomposed from CD information produced in the transmitter comparing to CA data. So, the CD data will be considered as imaginary data and the CA will be considered as real data to avoid any loss of data. That is all go via the DWT-OFDM receiver which can be seen in Figure5.
Figure 5. Proposed receiver model of DWT-OFDM.
Symmetry 2021, 13, x FOR PEER REVIEW 8 of 19
The specific processes to perform digital Pre-FDE were fully illustrated in [22]. In the OFDM transmitter (Tx), a Pre-FDE has performed adaptively according to the estimated SNR information.
Figure 4. A Wavelet-OFDM indoor imaging VLC-MIMO receiver by utilizing the hybrid Pre-SDE system.
2.3.1. Received Wavelet-OFDM Data DWT receiver is the reverse operation of IDWT. It performs a DWT according to the MATLAB toolbox. The transmitted signal ∪[k] is the front-end receiver successive digi- tally modulated symbols to be transformed, the data is decomposed across filters, LPF and HPF related to CA in addition to the CD. CA is the output signal of the approximation coefficients or LPF and CD are the output data of the detailed coefficients or HPF to exe- cute that operation. The data is transposed before converting from S/P and the wavelet family used is Demy wavelet. ∪ [k] contains some zeros elements which are decomposed from CD information Symmetry 2021, 13, 270 produced in the transmitter comparing to CA data. So, the CD data will be considered9 of19 as imaginary data and the CA will be considered as real data to avoid any loss of data. That is all go via the DWT-OFDM receiver which can be seen in Figure 5.
FigureFigure 5. 5.Proposed Proposed receiver receiver model model of of DWT-OFDM. DWT-OFDM.
The forward DWT is used to reconstruct the data symbol. The DWT-OFDM now could be obtained by a weighted sum of wavelet similar to scale carriers, which can be expressed as: ( ) = ( ) ( ) + ( ) W t ∑j≤J∑kwj,k t .Ψj,k t ∑k aj,kΦJ,k t , (12) th where wj,k(t) is the wavelet coefficients existing at the k position from the scale j, Ψj,k(t) is the IDWT modulator as a sequence of wavelet, ΦJ,k(t) is the scaling function, and aj,k is the approximation coefficient. The output of DWT consists of two vectors: the approximation coefficients (AC) and the detail coefficients (DC). The wavelet modulation (IDWT) in addition to the wavelet demodulation is obtained according to the Mallet algorithm [14]. The IDWT depends on up-sampling and DWT depends on down-sampling by a factor of two to each other, and filtering the approximated coefficients as shown in Equation (10), relative to LPF/HPF, respectively.
2.3.2. Pre-SDE Equations After propagation, the LED light is detected by an ImR or an ImADR that is supposed to calculate the photodetector responsivity, R. Firstly, the DC components are separated from the detected signals. Likewise, the electrical received signal vector y is acquired as [28]: y = RPopt ξ HA x + n, (13) where n acts as the noise with variance and zero mean mentioned in [26,29]. H is the matrix of the channel, as well as n, is the additive vector of the noise. The detector will provide both LOS [28,30] and diffuse components in exemplary indoor environments. Note, in the traditional indoor areas, the detector can obtain both LOS and diffuse components. Once the receiver is positioned at room corners, the diffuse components are highest, while it is still at least 7 dB lower than the weakest LOS components in electric power. This is the reason for considering LOS only as our concerning study.
2.3.3. SNR Estimation The calculated electrical signals of four SNRs are mentioned in [11]:
2 Rpopt ξaihii SNR0 = = a2SNR , i = 1, 2, ··· , N, (14) i σ2 1 i ni Symmetry 2021, 13, 270 10 of 19
2 2 2 where SNRi = (R Popt ξ hii) σ xi /σ ni, i = 1, 2, ... .., N. R is the photodiode responsivity, Popt is the LED output power, ξ is the modulation index, and hii are the elements of the H matrix. Therefore, Pre-SDE could be effectively employed [22]. To guarantee that the four received signals have the same SNR, as per Equation (4), the following requirement must be satisfied: 2 2 2 a 1SNR1 = a 2 SNR2 = ······ = a NSNRN. (15) By solving Equations (2) and (5), the N diagonal elements of matrix A can be obtained by [14,24]: s N ai = . (16) SNR ∑N 1 i i=1 SNRi The diagonal elements of A have been provided based on the SNR value of every channel. So, the evaluated SNR values must be supposed to be transferred back to the ceiling via the uplink channels as feedback information. For the most part, the optical channel remains static, while users of VLC are usually in fixed positions or moving with a slightly slower speed. Hence, the uplink delay will provide a relative effect on Pre-SDE performance. However, it should be observed that the electrical power of every channel ought to be flexible to adapt to a confirmed extent to efficiently excite the Pre-SDE, in the pre-coded VLC-MIMO techniques [31]. Thus, for VLC-MIMO systems, the Pre-SDE can be supposed like a precoding system. The SNR calculations are carried out in the wavelet-OFDM receiver (Rx), and the acquired SNR information is back to the wavelet-OFDM transmitter (Tx). By using the returns feedback of the SNR data, a sufficient modulation bandwidth could be offered to the OFDM data [22].
2.3.4. Received Data Estimation In the VLC-MIMO techniques, the performance cannot be optimized as a stable modulation form. In addition, just LoS channel gains through the transmitter to the receiver are utilized as channel matrix coefficients. However, the non-LoS is proposed to show the difference between the BER performance of OFDM and wavelet-OFDM. The nonlinearity and limited bandwidth properties of LEDs have not been considered to achieve the significant enhancement of channel diversity. To decode and recover data from the obtained signals, a channel estimation matrix technique is used. MIMO de-multiplexing is performed because of its low complication, zero-forcing (ZF) is implemented using basic channel observation [10]. Therefore, in order to get the transferred data estimation Xest and decompose the detected streams of Y, a matrix inversion and multiplication should be required: −1 Xest = H Y. (17) After ZF based MIMO de-multiplexing, we have [11]:
−1 −1 ex = H y = RPopt ξ HA x + H n, (18)
where H is the channel matrix for the 4 × 4 SDM-VLC system given by Equation (11) and n is the additive noise vector. Likewise, ni (i = 1, 2, ······ , N) is represented as a 2 zero-mean real-valued AWGN through variance σ ni. Both shot and thermal noises, given − by Equations (7) and (8), respectively, are jointly added to total noise. The notation (.) 1 indicates the inverse operator. diag(.) symbolized to the diagonal items, which are the results, are given within the rounded brackets. All blocks except the DWT-OFDM ones are explained in Ref. [22]. Here, we focus on DWT-OFDM blocks in detail. Symmetry 2021, 13, x FOR PEER REVIEW 10 of 19
pre-coded VLC-MIMO techniques [31]. Thus, for VLC-MIMO systems, the Pre-SDE can be supposed like a precoding system. The SNR calculations are carried out in the wavelet-OFDM receiver (Rx), and the acquired SNR information is back to the wavelet-OFDM transmitter (Tx). By using the returns feedback of the SNR data, a sufficient modulation bandwidth could be offered to the OFDM data [22].
2.3.4. Received Data Estimation In the VLC-MIMO techniques, the performance cannot be optimized as a stable mod- ulation form. In addition, just LoS channel gains through the transmitter to the receiver are utilized as channel matrix coefficients. However, the non-LoS is proposed to show the difference between the BER performance of OFDM and wavelet-OFDM. The nonlinearity and limited bandwidth properties of LEDs have not been considered to achieve the sig- nificant enhancement of channel diversity. To decode and recover data from the obtained signals, a channel estimation matrix technique is used. MIMO de-multiplexing is per- formed because of its low complication, zero-forcing (ZF) is implemented using basic channel observation [10]. Therefore, in order to get the transferred data estimation Xest and decompose the detected streams of Y, a matrix inversion and multiplication should be required: