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

Enhanced performance of MIMO multi-branch hybrid neural network in single receiver MIMO visible light communication system

Open Access Open Access

Abstract

For the single receiver multiple-input-multiple-output (SR-MIMO) visible light communication (VLC) system, the superposing of two transmitters will introduce severe distortion in the time-domain and frequency-domain. In this paper, we first proposed a MIMO multi-branch hybrid neural network (MIMO-MBNN) as the post-equalizer in the SR-MIMO pulse amplitude magnitude eight levels (PAM8) VLC system. Compared with the traditional single-input-single-output least mean square equalizer with Volterra series (SISO-LMS) and SISO deep neural network (SISO-DNN), MIMO-MBNN can achieve at most 3.35 dB Q factor improvement. Furthermore, the operation range of MIMO-MBNN is at least 2.33 times of SISO-DNN and SISO-LMS among the measured signal peak to peak voltage. At last, 2.1 Gbps data rate is achieved by MIMO-MBNN below the 7% hard-decision forward error correction (HD-FEC) threshold. As far as we know, this is the highest data rate in the SR-MIMO VLC system.

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

1. Introduction

With the advent of the 5G era, indoor mobile traffic will show explosive growth. VLC is a promising wireless communication technology to solve this problem [1]. VLC uses a light-emitting diode (LED) as the transmitter and PIN or avalanche photon diode (APD) as the receiver to implement data transmission [2,3]. Its many advantages, such as no spectrum limitation, energy and cost-efficient, non-line of sight (NLOS) encryption invalid, and the combination with lighting make it a hotspot in the field of optical wireless research [4,5].

However, the 3dB bandwidth of LED is just 10∼20 MHz, which is a big obstacle to realize a high-speed VLC system [6,7]. In recent years, MIMO VLC has become a hotspot because it can effectively increase the capacity of the system at a low cost. Furthermore, the concept of MIMO accommodates well with the concept of massive access, uplinks, and downlinks in 5G. [8,9]. MIMO-VLC can be divided into two fields: imaging and non-imaging [10]. The imaging MIMO can provide a full-rank channel matrix. However, the requirement for the alignment between LED and detector is strict, which restricts the practical utilization of imaging MIMO. Therefore, most of the current MIMO-VLC research focuses on non-imaging MIMO. However, the correlation of the VLC signal is strong, and the decorrelation may produce an ill-conditioned channel matrix [11]. Traditional MIMO algorithms, such as singular value decomposition (SVD), orthogonal circulant matrix transform (OCT) require a perfect channel matrix [12,13]. Therefore, not all data streams can be decoded in the MIMO-VLC system. Interestingly, the superposition coding modulation (SCM) strategy provides us with another possibility, and it has been proved to be feasible in the MIMO-VLC system. Furthermore, in some scenarios such as the internet of vehicles, the car lights operate at high power. The nonlinear effect of the LED makes the transmitter unable to support high order modulation format. By utilizing two LEDs and SCM, we can reduce the order of modulation format on each transmitter without sacrificing the overall spectrum efficiency. In [14], the author superposes two non-equal probability PAM4 signal to generate an equiprobable PAM7 signal, whose data rate is much higher than the space-time block coding (STBC) PAM4 signal. However, the PAM4 signal is superposed in one dimension, and if the source entropy continues to increase, the system performance will deteriorate dramatically. In [15], the author used QPSK and 16QAM to realize a 64QAM SCM constellation and achieved 1.5 Gbps data transmission with the aid of the look-up table and successive interference cancellation. In [16], the author used a geometrical shaping 8QAM square constellation point combined with 4QAM to implement 32QAM superposition coded modulation, which showed its superiority to the other 32QAM SCM constellation points. However, QAM-based superposition modulation requires both receivers to adopt a differentiated design. For high-order 64QAM, it can only work at a power difference of more than four times. At the same time, big power diversity will introduce power competition between two signals, and the signal-to-noise ratio (SNR) of the signal with lower power will further deteriorate. If the SCM is conducted by two PAM8 streams, we can get a 64QAM constellation at the receiver side as well, and the optimal power ratio is 1 to 1, which will not introduce additional power competition. However, the research of PAM-based SCM signals in MIMO-VLC systems is still lacking.

On the other hand, the nonlinear effect in the VLC system will distort the signal severely [17]. For the SCM signal, it is more susceptible to the nonlinear distortion of the VLC channel. However, the transfer function of the VLC channel is complicated, and it is difficult to be expressed by analytical equations [18]. At this time, with the help of the powerful nonlinear mapping ability of the neural network, the performance of the received signal will be improved to a big degree. In [19], 1.5 Gbps PAM8 transmission over 1.2 m underwater VLC channel is implemented by utilizing a SISO deep neural network (SISO-DNN) as a post equalizer. In [20], the author designed a two tributary heterogeneous neural network (TTHNET) to emulate the VLC channel. The parameters of TTHNET are only 0.8% of the regular multilayer perceptron (MLP). Therefore, it is foreseeable that after the distortion of the MIMO VLC channel, the received SCM signal could be compensated by the neural network. Unfortunately, there is a lack of research on the neural network as a post-equalizer in the MIMO VLC system.

In this paper, we proposed a novel MIMO multi-branch hybrid neural network (MIMO-MBNN) post equalizer and implemented it in a 64QAM SCM and carrier-less amplitude-phase (CAP) modulation single-receiver MIMO (SR-MIMO) VLC system. The superposed 64QAM constellation at the receiver consists of two PAM8 signals from two transmitters. The BER performance is continuous in all the operation range, which is different from that in [15]. Furthermore, the proposed MIMO-MBNN post equalizer shows great superiority to traditional SISO least mean square with Volterra series (SISO-LMS) and SISO-DNN post equalizers when the power of two transmitters are unbalanced. The innovation points are listed as follow:

  • (1) This article marks the first time that MIMO-MBNN is proposed and utilized in 64QAM superposed coded modulation SR-MIMO VLC system.
  • (2) The Q performance of MIMO-MBNN is superior to SISO-LMS and SISO-DNN in all the operation range, and it achieved at most 3.35 dB Q factor improvement. When Vpp1=250 mV, 350 mV, and 600 mV, the valid operation range above the 2×10−2 HD-FEC BER threshold of MIMO-MBNN is at least 2.33 times that of SISO-LMS and SISO-DNN.
  • (3) The highest data rate of 2.1 Gbps of the SR-MIMO VLC system is achieved in this paper. Compared with the data rate of 1.5 Gbps in [15], the highest data rate in the SR-MIMO VLC system is improved by 40 percent.

2. Principles and experimental setup

2.1 Principles of superposed 64QAM constellation

The principles of SCM 64QAM constellations are shown in Fig. 1. In this paper, we define the signal peak to peak voltage of LED1 and LED2 as Vpp1 and Vpp2, respectively. The power ratio ${\mathrm{\beta}} \textrm{ = }{V_{pp2}}/{V_{pp1}}$. The PAM8 SCM scheme is shown in Fig. 1(a). The in-phase part and quadrature part of SCM 64QAM constellation at the receiver side is transmitted by LED1 and LED2, respectively. The original bitstreams of LED1 and LED2 can be modulated into PAM8 symbols, which can be expressed as:

$${\mathbf T} = [{T_1}(t),{T_2}(t)]\;,\;{T_k}(t) \in \{{\pm} 1, \pm 3, \pm 5, \pm 7\}. $$

Here ${T_1}(t)$ and ${T_2}(t)$ are the PAM8 symbol column vectors of LED1, and LED2 respectively. $t = 1,2,\ldots ,N$ and N is the total number of transmitted symbols. As the VLC system is intensity modulation and direct detection (IM/DD) system and the signal are modulated by carrier-less amplitude magnitude (CAP) modulation format, the PAM8 signal in LED1 and LED2 are multiplied by cosine subcarrier and sine subcarrier respectively to make them orthogonal. Therefore, the superposed signal at the receiver is a QAM constellation. If we consider no distortion and just additive white gaussian noise in the channel, the received constellation for each transmitted symbol after post equalization in Fig. 1(a) can be expressed as:

$$Y(t) = {Y_1}(t) + {\mathrm{\beta}} {Y_2}(t) \cdot i. $$

 figure: Fig. 1.

Fig. 1. Principles of SCM 64QAM schemes. (a) PAM8. (b) QPSK and 16QAM in [15].

Download Full Size | PDF

As a comparison, the QPSK and 16QAM based 64QAM SCM scheme in [15] is shown in Fig. 1(b). From Fig. 1, we can conclude the differences between these two schemes as follow: (i) The optimal power ratio between LED1 and LED2 of the PAM8 SCM 64QAM scheme in Fig. 1(a) is 1, while that in Fig. 1(b) is 4. This indicates that the system in Fig. 1(a) is more conducive to practical application, as the two transmitters do not require a differentiated design. (ii) The valid theoretical value of ${\mathrm{\beta}} $ in Fig. 1(b) is below 0.5 or above 3. When the power ratio is not in the valid operation range, the signal cannot be demodulated. While that in Fig. 1(a) is not limited. Furthermore, the power competition between QPSK and 16QAM will further deteriorate the overall system performance. Therefore, we believe that the PAM8 SCM 64QAM scheme is a better choice for the SR-MIMO VLC system.

Furthermore, we choose the CAP instead of the OFDM as the modulation format because the CAP has a relatively lower peak to average power ratio (PAPR) than the OFDM. As two signals will superpose at the receiver, higher PAPR will clip the signal dramatically and deteriorate the system performance. Therefore, the CAP is superior to the OFDM in the SCM VLC system. Additionally, by utilizing a square-root raised cosine low-pass filter, the CAP signal is more adapted to the high-frequency fading VLC channel than the conventional OFDM signal. Therefore, we choose the CAP as the modulation format in this paper.

2.2 Data flow and experimental setup of SR-MIMO VLC system

The data flow and experimental setup of the SR-MIMO VLC system are shown in Fig. 2. To simplify the expression, we abbreviate Transmitter 1 as Tx1, Transmitter 2 as Tx2, and receiver as Rx. The transmitted signals before up-conversion can be expressed as:

$$\left[ {\begin{array}{l} {{S_1}(t)}\\ {{S_2}(t)} \end{array}} \right] = \left[ {\begin{array}{c} {{g_1}(t) \otimes T_1^{up}(t)}\\ {{g_2}(t) \otimes T_2^{up}(t)} \end{array}} \right]. $$
where ${S_1}(t)$ and ${S_2}(t)$ are the 1st and 2nd transmitted signal before up-conversion respectively. The ${g_1}(t)$ and ${g_1}(t)$ are the square-root raised cosine low-pass filters in [4], and the roll-off factor is 0.205, which is selected according to the results in [21,22]. The ${\otimes} $ represent for convolution operation. $T_1^{UP}(t)$, $T_2^{UP}(t)$ are the PAM8 symbols after upsampling. In this paper, the upsample times is 4. Then after up-conversion, the transmitted signals from LED1 and LED2 can be written as:
$$\left[ {\begin{array}{c} {{X_1}(t)}\\ {{X_2}(t)} \end{array}} \right] = \left[ {\begin{array}{cc} {\cos (2\pi {f_0}t)}&0\\ 0&{ - \sin (2\pi {f_0}t)} \end{array}} \right]\left[ {\begin{array}{c} {{S_1}(t)}\\ {{S_2}(t)} \end{array}} \right]. $$
where ${X_1}(t)$ and ${X_2}(t)$ are the transmitted data streams loaded by the arbitrary waveform generator (AWG). ${f_0}$ is the up-conversion center frequency. Then the two data streams are generated by the AWG M9502. A self-designed T-bridge passive hardware equalizer is connected to the output of AWG to pre-emphasize the high-frequency components of the electrical signals. In this way the system bandwidth can be increased from several ten MHz to 350 MHz. The pre-equalized signal is amplified by an electrical amplifier (EA) before it injected into the LED1 and LED2.

 figure: Fig. 2.

Fig. 2. Data flow and experimental setup of the SR-MIMO CAP VLC system. AWG: arbitrary waveform generator; Eq: hardware equalizer; EA: electrical amplifier; OSC: oscilloscope.

Download Full Size | PDF

The received signal after the VLC channel can be expressed as:

$$\left[ {\begin{array}{c} {{R_1}(t)}\\ {{R_2}(t)} \end{array}} \right] = \left[ {\begin{array}{cc} {{H_{11}}}&{{H_{12}}}\\ {{H_{21}}}&{{H_{22}}} \end{array}} \right]\left[ {\begin{array}{c} {{f_1}({X_1}(t))}\\ {{f_2}({X_2}(t))} \end{array}} \right]\textrm{ + }\left[ {\begin{array}{c} {{N_1}(t)}\\ {{N_2}(t)} \end{array}} \right]. $$
where ${R_1}(t)$ and ${R_2}(t)$ are the received signal from LED1 and LED2 respectively. ${H_{ij}}$ represent for the power distortion from the ith transmitter to the jth transmitter (i=1,2; j=1,2). ${N_1}(t)$ and ${N_2}(t)$ are the noise of Tx1 and Tx2 at the receiver. ${f_1}(t)$ and ${f_2}(t)$ represent the amplitude magnitude (AM) response from LED1 and LED2 to Rx, respectively. To simplify the analyze, we treat ${f_1}(x) = {f_2}(x) = f(x)$, which can be expressed as:
$$f(x) = \left\{ {\begin{array}{l} x\\ {({a_1}x + {a_2}{x^2} + {a_3}{x^3})/b} \end{array}\begin{array}{l} {,\;low\;{V_{pp}}}\\ {,\;high\;{V_{pp}}} \end{array}} \right.. $$

The typical value for ${a_1}$, ${a_2}$, ${a_3}$, and b are 1.259, 0.01373,-0.4576 and 0.8165 respectively, which are measured at the 1.4V signal peak to peak voltage (Vpp). After power normalization, ${H_{11}} = 1$ and ${H_{22}}\textrm{ = }{\mathrm{\beta}} $. ${N_1}(t)$ and ${N_2}(t)$ denote the system noise. To reduce the complexity of the analysis, we think that the mutual influence of LED1 and LED2 at the receiver is 0. Therefore, Eq. (5) can be rewritten as:

$$\left[ {\begin{array}{c} {{R_1}(t)}\\ {{R_2}(t)} \end{array}} \right] = \left[ {\begin{array}{cc} 1&0\\ 0&{\mathrm{\beta}} \end{array}} \right]\left[ {\begin{array}{c} {{X_1}(t) \otimes {h_1}(t)}\\ {{X_2}(t) \otimes {h_2}(t)} \end{array}} \right]\textrm{ + }\left[ {\begin{array}{c} {{N_1}(t)}\\ {{N_2}(t)} \end{array}} \right]. $$
where ${h_1}(t)$ and ${h_2}(t)$ are the time-domain response of the VLC channel from LED1 and LED2 to PIN respectively. To simplify the analysis, we treat ${h_1}(t) = {h_2}(t) = h(t)$, which means the two LEDs have the same time-domain response. In the following part, we treat the system operating at low Vpp. Therefore, the superposed signal $Y(t)$ at the receiver can be written as:
$$\begin{array}{l} Y(t) = f({R_1}(t) + {R_2}(t)) + N(t)\\ \;\;\;\;\;\;{\kern 1pt} \textrm{ = }{\kern 1pt} {\kern 1pt} f({X_1}(t) \otimes h(t) + {\mathrm{\beta}} \cdot {X_2}(t) \otimes h(t)) + N(t) \end{array}. $$
where $N(t)$ is the overall noise power. If we consider the AM response from LED to PIN is linear and ${f_y}(y) = y$, then the signal ${Y_1}(t)$ and ${Y_2}(t)$ after the multiplication of cosine subcarrier and negative sine subcarrier respectively can be expressed as:
$$\begin{array}{l} \left[ {\begin{array}{l} {{Y_1}(t)}\\ {{Y_2}(t)} \end{array}} \right] = \left[ {\begin{array}{l} {\cos (2\pi {f_0}t)}\\ { - \sin (2\pi {f_0}t)} \end{array}} \right] \cdot Y(t)\\ \;\;\;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} {\kern 1pt} = \left[ {\begin{array}{l} {(({S_1}(t) \cdot \cos ({{\mathrm{\omega}}_0}t)) \otimes h(t) + {\mathrm{\beta}} \cdot ({S_2}(t) \cdot \sin ({{\mathrm{\omega}}_0}t)) \otimes h(t)) \cdot \cos ({{\mathrm{\omega}}_0}t)}\\ { - (({S_1}(t) \cdot \cos ({{\mathrm{\omega}}_0}t)) \otimes h(t) + {\mathrm{\beta}} \cdot ({S_2}(t) \cdot \sin ({{\mathrm{\omega}}_0}t)) \otimes h(t)) \cdot \sin ({{\mathrm{\omega}}_0}t)} \end{array}} \right] + \left[ {\begin{array}{l} {{N_1}(t)}\\ {{N_2}(t)} \end{array}} \right] \end{array}. $$
where ${{\mathrm{{\mathrm{\omega}}}} _\textrm{0}}\textrm{ = 2}\pi {f_0}$. After deduction and simplification which are shown in the Appendix, Eq. (9) can be rewritten as:
$$\begin{array}{l} {Y_1}(t) = \frac{1}{2}{S_1}(t) \otimes (h(t)\cos ({{\mathrm{\omega}} _0}t)) + \frac{{\mathrm{\beta}} }{2}{S_2}(t) \otimes (h(t)\sin ({{\mathrm{\omega}} _0}t))\\ \;\;\;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} - \frac{1}{2}\sqrt {{S_1}{{(t)}^2} + {{({\mathrm{\beta}} {S_2}(t))}^2}} \sin (2{{\mathrm{\omega}} _0}t + {{\mathrm{\alpha}} _1}) \otimes (h(t)\cos ({{\mathrm{\omega}} _0}t))\\ \;\;\;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} - \frac{1}{2}\sqrt {{S_1}{{(t)}^2} + {{({\mathrm{\beta}} {S_2}(t))}^2}} \sin (2{{\mathrm{\omega}} _0}t + {{\mathrm{\alpha}} _2}) \otimes (h(t)\sin ({{\mathrm{\omega}} _0}t)) + {N_1}(t) \end{array}. $$
$$\begin{array}{l} {Y_2}(t) ={-} \frac{1}{2}{S_1}(t) \otimes (h(t)\sin ({{\mathrm{\omega}} _0}t)) + \frac{{\mathrm{\beta}} }{2}{S_2}(t) \otimes (h(t)\cos ({{\mathrm{\omega}} _0}t))\\ \;\;\;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} - \frac{1}{2}\sqrt {{S_1}{{(t)}^2} + {{({\mathrm{\beta}} {S_2}(t))}^2}} \sin (2{{\mathrm{\omega}} _0}t + {{\mathrm{\alpha}} _2}) \otimes (h(t)\cos ({{\mathrm{\omega}} _0}t))\\ \;\;\;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} + \frac{1}{2}\sqrt {{S_1}{{(t)}^2} + {{({\mathrm{\beta}} {S_2}(t))}^2}} \sin (2{{\mathrm{\omega}} _0}t + {{\mathrm{\alpha}} _1}) \otimes (h(t)\sin ({{\mathrm{\omega}} _0}t)) + {N_2}(t) \end{array}. $$

Here $\tan {{\mathrm{\alpha}} _1} ={-} \frac{{{S_1}(t)}}{{{\mathrm{\beta}} {S_2}(t)}},\tan {{\mathrm{\alpha}} _2} = \frac{{{\mathrm{\beta}} {S_2}(t)}}{{{S_1}(t)}}$. The first two terms of Eq. (10) and Eq. (11) are linear in-band frequency response, which can be treated as the convolution of origin signal ${S_1}(t),{S_2}(t)$ and channel response. As we usually consider the VLC system as a non-time-varying system, the channel response can be treated as a finite impulse response (FIR) filter. Therefore, the first term can be equalized by a linear equalizer. Additionally, for ${Y_1}(t)$, the second term contains ${S_2}(t)$ can be considered as noise and vice versa. Therefore, the value of ${\mathrm{\beta}} $ will significantly affect the SNR of ${S_1}(t)$ and ${S_2}(t)$. The bigger the ${\mathrm{\beta}} $ is, the higher the SNR of ${S_2}(t)$ is. At the same time, the SNR of ${S_1}(t)$ deteriorates. For ${S_1}(t)$, the condition is the contrary. The last two terms are high-frequency components, which can be removed by the square-root raised cosine low-pass filter. According to Eqs. (6), (10) and (11), the expression of ${Y_1}(t)$ and ${Y_2}(t)$ can be easily derived by replacing ${S_1}(t)$ with ${a_2}{S_1}^2(t)/b$ or ${a_3}{S_1}^3(t)/b$ when the system operates at high Vpp. Here we define Eq. (10) and Eq. (11) as:

$$\begin{array}{l} {Y_1}(t) = {F_1}({S_1}(t),{S_2}(t))\\ {Y_2}(t) = {F_2}({S_2}(t),{S_1}(t)) \end{array}$$
Then the expression of ${Y_1}(t)$ and ${Y_2}(t)$ can be written as:
$$\begin{array}{l} {Y_1}(t) = \frac{{{a_1}}}{b}{F_1}({S_1}(t),{S_2}(t)) + \frac{{{a_2}}}{b}{F_1}(S_1^2(t),S_2^2(t)) + \frac{{{a_3}}}{b}{F_1}(S_1^3(t),S_2^3(t))\\ {Y_2}(t) = \frac{{{a_1}}}{b}{F_2}({S_2}(t),{S_1}(t)) + \frac{{{a_2}}}{b}{F_2}(S_2^2(t),S_1^2(t)) + \frac{{{a_3}}}{b}{F_2}(S_2^3(t),S_1^3(t)) \end{array}$$
Therefore, according to Eqs. (6), (10) and (11), the ${Y_1}(t)$ without the high-frequency components can be depicted as:
$$\begin{array}{l} {Y_1}(t) = \underbrace{{\frac{{{a_1}}}{b}(\frac{1}{2}{S_1}(t) \otimes (h(t)\cos ({{\mathrm{\omega}} _0}t))}}_{{linear\;distortion}} + \underbrace{{\frac{{\mathrm{\beta}} }{2}{S_2}(t) \otimes (h(t)\sin ({{\mathrm{\omega}} _0}t))}}_{{linear\;crosstalk}} + \\ \;\;\;\;\;\;\;\;\;\underbrace{{\frac{{{a_2}}}{b}(\frac{1}{2}S_1^2(t) \otimes (h(t)\cos ({{\mathrm{\omega}} _0}t)) + \frac{{\mathrm{\beta}} }{2}S_2^2(t) \otimes (h(t)\sin ({{\mathrm{\omega}} _0}t))}}_{{nonlinear\;distortion\;and\;crosstalk}} + \\ \;\;\;\;\;\;\;\;\;\underbrace{{\frac{{{a_3}}}{b}(\frac{1}{2}S_1^3(t) \otimes (h(t)\cos ({{\mathrm{\omega}} _0}t)) + \frac{{\mathrm{\beta}} }{2}S_2^3(t) \otimes (h(t)\sin ({{\mathrm{\omega}} _0}t))}}_{{nonlinear\;distortion\;and\;crosstalk}} \end{array}$$

The expression of ${Y_2}(t)$ is similar to that of Eq. (14). From the above equation, we can see that the received signal after multiplying by $\cos ({{\mathrm{\omega}} _0}t)$ and $\textrm{ - }\sin ({{\mathrm{\omega}} _0}t)$ can be divided into three parts: linear distortion, linear crosstalk, nonlinear distortion and crosstalk. The linear distortion can be successfully compensated by a linear equalizer. For the mitigating of the nonlinear distortion and crosstalk, an additional effective nonlinear equalizer must be considered. Our experimental setup consists of two independent input signals. Each signal should be compensated by two neural networks, one is applied for the linear distortion, the other is for nonlinear distortion. In this way, a total 4 tributaries NN should be used for the signal compensation. Since two input signals are combined together at the single receiver after 1.2m free-space transmission, we assume using only one neural network to compensate the nonlinear distortion for the two signals simultaneously, thus a total triple branches NN can be implemented for post-equalization. Our following experiment results evidently prove the viability of our assumption.

After the processing of three equalizers, the two data streams are loaded by the match filter pairs ${m_1}(t)$ and ${m_2}(t)$ in the CAP modulation. Then the two signals are equalized by two LMS equalizers. Different from the SISO-LMS in case 1, these additional equalizers don’t have the Volterra series, which cannot compensate for the nonlinear distortion. They are utilized to mitigate residual noise. The input taps of the additional two LMS equalizers are 33, and the step is set to 0.007. This way of compensating the received signal with two-stage equalizers has been proven effective in the VLC system [22,23]. Finally, the signals are demodulated and recovered to bitstreams again. The details about the components in Fig. 2 are listed in Table 1.

Tables Icon

Table 1. Details of components in Fig. 2.

2.3. Structure of MIMO-MBNN

Figure 3 is the forward propagation of MIMO-MBNN in the SR-MIMO VLC system. The symbols ${Y_i}(t)$, ${S_i}(t)$ and ${P_i}(t)$ are consistent with that in Fig. 2 (i=1,2). For the robustness of the neural network, we should normalize the signals at first. In this paper, we utilize the ‘max absolute value’ normalization function, which can be written as:

$${F_{\max (abs)}}(\overrightarrow x ) = \overrightarrow x /\max (abs(\overrightarrow x )). $$
Here $\overrightarrow x $ is the signal vector. Unless otherwise specified, all vectors default to column vectors in the following part.

 figure: Fig. 3.

Fig. 3. Structure of MIMO-MBNN

Download Full Size | PDF

To simplify the principle chart of MIMO-MBNN in Fig. 3, we set the input length of each data stream to 3. In the experiment, the optimal input data length is 53. Compared to the input length of 191 in [20], the input length in this paper is acceptable. Here we define the input vector as $\overrightarrow {Y(t)} = {[{Y_1}(t - 1),Y{(t)_1},{Y_1}(t + 1),{Y_2}(t - 1),{Y_2}(t),{Y_2}(t + 1)]^T}$ and $t = 2,3,\ldots N - 1$, where N is the length of the training data set. Afterward, the data streams are concatenated and loaded by the 1st, 2nd linear branch, and the 3rd nonlinear branch. Then the linear and nonlinear output of 1st, 2nd, and 3rd branches can be expressed as:

$$\left\{ {\begin{array}{l} {{Y_{1,L}}(t) = {\mathbf W}_1^{(1)\textrm{T}} \mathop{{Y(t)}}\limits^\rightharpoonup + \overrightarrow {b_1^{(1)}} }\\ {\left[ {\begin{array}{l} {{Y_{1,NL}}(t)}\\ {{Y_{2,NL}}(t)} \end{array}} \right] = {\mathbf W}_3^{(3)T}R ({\mathbf W}_3^{(2)T}R ({\mathbf W}_3^{(1)\textrm{T}}\overrightarrow {Y(t)} + \overrightarrow {b_3^{(1)}} ) + \overrightarrow {b_3^{(2)}} ) + \overrightarrow {b_3^{(3)}} }\\ {{Y_{2,L}}(t) = {\mathbf W}_2^{(1)\textrm{T}} \mathop{{Y(t)}}\limits^\rightharpoonup + \overrightarrow {b_2^{(1)}} } \end{array}} \right.$$
where ${\mathbf W}_k^{(n)\textrm{T}}$ and $\overrightarrow {b_k^{(n)}}$ are the nth layer’s weight matrix and bias vector of the kth branch respectively. T represents the matrix transpose. $R(x)$ is the famous ‘relu’ function, and it is expressed as:
$$R(x) = \left\{ {\begin{array}{l} {x,{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} x > 0}\\ {0,{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} x \le 0} \end{array}} \right.. $$
Another activation function mentioned in this paper is ‘tanh’, which is written as:
$$\tanh (x) = \frac{{{e^x} - {e^{ - x}}}}{{{e^x} + {e^{ - x}}}}. $$
The 1st and 2nd output of the MIMO-MBNN can be written as:
$${P_1}(t) = {Y_{1,L}}(t) + {Y_{1,NL}}(t),{P_2}(t) = {Y_{2,L}}(t) + {Y_{2,NL}}(t). $$
The loss function of MIMO-MBNN is the mean square error (MSE) function. It will calculate the differences between the received signal and transmitted signal before up-conversion in Fig. 3, which can be expressed as:
$${\mathbf W}_k^{(n)},\overrightarrow {b_k^{(n)}} = \mathop {\arg \min }\limits_{{\mathbf W}_k^{(n)},\overrightarrow {b_k^{(n)}} } \frac{1}{N}\sum\limits_{t = 2}^{N - 1} {{{||{{P_1}(t) - {S_1}(t)} ||}^2}} + {||{{P_2}(t) - {S_2}(t)} ||^2}. $$
The optimizer utilized in this paper is adaptive moment estimation (Adam). We utilize the ‘Checkpoint’ tool in Keras to record the training information and utilizing the model with the smallest validation loss to predict the signal. If the neural network began to overfit, the training is actually ‘finished’ and the model used for prediction is the model before overfitting. Therefore, ‘Checkpoint’ can also prevent overfitting to some extent.

3. Results and discussion

Firstly, to show the nonlinear response of the system and the interaction between the two signals, we measured the AM response between the transmitted signal and the received signal. For each point in Fig. 4, the X-axis value is the amplitude of normalized ${X_1}(t)$ in Fig. 2, while the Y-axis value is that of normalized $Y(t)$ in Fig. 2. The normalization way is shown in Eq. (15). Figure 4(a) is the AM response without the interference from Tx2. From Fig. 4(a), we can see that at the lower Vpp1 of 300 mV, the AM response is almost linear. However, with the increase of Vpp1, the AM response at high signal voltage shows apparent nonlinear response between ${X_1}(t)$ and $Y(t)$. Figure 4(b) shows the AM response between Tx1 and Rx at different values of Tx1: Tx2 when Vpp1 is 900 mV. We can easily recognize that with the increase of power ratio of Tx1 to Tx2, the noise power brought by Tx2 is continuously decreasing. This indicates that the signal from Tx1 will inevitably be distorted by that from Tx2 and vice versa.

 figure: Fig. 4.

Fig. 4. amplitude magnitude response. (a) Tx1 and Rx without Tx2 at different Vpp. (b) Tx1 and Rx with varying values of Tx1: Tx2 (Vpp1=900 mV)

Download Full Size | PDF

Afterward, the bit error rate (BER) performance versus hidden node numbers and layers in MIMO-MBNN are shown in Fig. 5. We measured the BER performance of MIMO-MBNN with different activation functions and different hidden layers. The MIMO-MBNN model at each point is not overfitting, and ‘Checkpoint’ is utilized to ensure optimal performance of MIMO-MBNN. Additionally, the loss function we utilized is the mean square error (MSE). Figure 5(a) is the BER performance with just 1 hidden layer. The best BER performance is achieved by ‘relu’, and the 1st node number is 136. Therefore, we set the node number and activation function of 1st layer to 136 and ‘relu’, respectively. Figure 5(b) shows the BER performance with two hidden layers under different activation function combinations. The 1st layer node number is fixed to 136. We can conclude from Fig. 5(b) that the optimal choice of the activation function is ‘relu’ + ‘relu’, and the optimal node number of the 2nd hidden layer is 104. We have verified that further increasing the number of nodes and layers will not significantly improve the system performance. Therefore, the layer number, 1st hidden layer node number, 2nd hidden layer node number, and activation function are fixed to 2, 136, 104, and ‘relu’. The parameter details are shown in Table 1. The epoch, batch size, training set ratio, and total sample number are 40, 512, 50%, and 73728, respectively. Therefore, the training set length and testing set length for SISO-LMS, SISO-DNN, and MIMO-MBNN are all 36864. After training, the parameters of MIMO-MBNN, SISO-LMS, and SISO-DNN are fixed. Then we utilize a new set of data to measure the BER or Q performance of the above equalizers and the length is 73728. The computation complexity in terms of the total parameters can be found in Table 2, where N is the input taps of SISO-LMS and SISO-DNN, ${N_v}$ is the input taps of the Volterra series, ${H_1}$, ${H_2}$, and ${H_3}$ are the node number of hidden layers. Clearly, the value of the total parameters depends on the input taps and the node number of hidden layers. In the SISO-LMS and SISO-DNN cases, the optimum input taps is 53, while the optimum input taps for the Volterra series is obtained at 9. Hence the total parameters of SISO-LMS, SISO-DNN, and MIMO-MBNN are 689, 46257, and 29224 respectively. It should be noted that although LMS boasted the lower complexity, the achieved BER performance is above the FEC threshold, which will result in failure detection. The complexity of MBNN is merely 63% of SISO-DNN, clearly validating the feasibility of the proposed scheme.

 figure: Fig. 5.

Fig. 5. BER performance versus different node numbers and activation function of MIMO-MBNN. (a) BER performance of 1 hidden layer. (b) BER performance of 2 hidden layers

Download Full Size | PDF

Tables Icon

Table 2. Detailed parameters of SISO-LMS, SISO-DNN, and MIMO-MBNN.

Then we measured the frequency response and spectrum mismatch of $Y(t)$, SISO-LMS, SISO-DNN, and MIMO-MBNN, which are shown in Figs. 6(a) and 6(b) respectively. The spectrum mismatch $M({\mathrm{\omega}} )$ is depicted as:

$$M({\mathrm{\omega}} ) = {|{X({\mathrm{\omega}} ) - P({\mathrm{\omega}} )} |^2}$$
where $P({\mathrm{\omega}} )$ is the frequency-domain response of $P(t)$ in Fig. 2. $X({\mathrm{\omega}} )$ is the frequency-domain response of the sum of ${X_1}(t)$ and ${X_2}(t)$. The spectrum mismatch in Fig. 6(b) shows that MIMO-MBNN outperforms SISO-LMS and SISO-DNN.

 figure: Fig. 6.

Fig. 6. Frequency and spectrum mismatch of the received signal, SISO-LMS, SISO-DNN, and MIMO-MBNN. (a) Frequency response. (b) spectrum mismatch

Download Full Size | PDF

The Q performance versus data rate is shown in Fig. 7. The Q in dB can be expressed as:

$$Q(dB) = 20 \cdot {\log _{10}}(\sqrt 2 \cdot erf{c^{ - 1}}(2BER))$$
The $erf{c^{\textrm{ - }1}}(x)$ is the inverse function of the complementary error function. The Vpp1 and Vpp2 are equal to 350 mV, which is the optimal voltage operating point in the SR-MIMO VLC system. The Q performance of MIMO-MBNN outperforms SISO-DNN and SISO-LMS, and the operation range of MIMO-MBNN above the 3.8×10−3 hard-decision forward error correction BER threshold (HD-FEC) is 165 Mbps, which is bigger than that of SISO-DNN and SISO-LMS. Additionally, the highest data rate of 2.1 Gbps is achieved by MIMO-MBNN. To the best of our knowledge, this is the highest data rate of the SR-MIMO VLC system [15]. To make the results more convincing, We also utilized the LMS with Volterra series after demodulation using $\cos ({{\mathrm{\omega}} _0}t)$ and $\textrm{ - }\sin ({{\mathrm{\omega}} _0}t)$, which named as SISO-LMS-Complex in Fig. 7. For SISO-LMS-Complex, the two signal streams after multiplying by $\cos ({{\mathrm{\omega}} _0}t)$ and $\textrm{ - }\sin ({{\mathrm{\omega}} _0}t)$ are combined into a complex signal stream and equalized by the LMS equalizer with Volterra series. From the results in Fig. 7, we can see that the Q performance of SISO-LMS is comparable to that of SISO-LMS-Complex. In this paper, we choose SISO-LMS to measure the system performance. As the nonlinear distortion in the SR-MIMO VLC system is severe, the SD-FEC is not very suitable for the SR-MIMO VLC system. Therefore, we choose HD-FEC as the figure of merit to measure the system performance.

 figure: Fig. 7.

Fig. 7. Q factor versus data rate

Download Full Size | PDF

Afterward, we measured the system $- {\log _{10}}(BER)$ performance at different Vpp1 and Vpp2. The results are shown in the contours in Fig. 8. The data rate is 2.1 Gbps. We can easily recognize that the system performance after equalization of MIMO-MBNN is better than that of SISO-DNN and SISO-LMS. When Vpp1 is between 250 mV and 350 mV and Vpp2 is between 350 mV and 450 mV, the system performance equalized by MIMO-MBNN owns an operation range above the 3.8×10−3 HD-FEC, while SISO-LMS and SISO-DNN own no operation range. Besides, the operation range above the 2×10−2 HD-FEC threshold of MIMO-MBNN is much bigger than that of SISO-LMS and SISO-DNN as well.

 figure: Fig. 8.

Fig. 8. Contour of system performance SISO-LMS, SISO-DNN, and MIMO-MBNN (from 1 to 3 in each row). (a1)-(a3). Average BER performance. (b1)-(b3). Tx1 BER performance. (c1)-(c3). Tx2 BER performance.

Download Full Size | PDF

To further investigate the system performance at different operation points, we measured the system Q performance of LED1, LED2, and the average system Q performance of LED1 and LED2 versus ${\mathrm{\beta}} $ when Vpp1=250 mV, 350 mV, and 600 mV. The results are shown in Fig. 9. When Vpp1=250 mV as Fig. 9(i)–9(iii) show, the operation range above the 2×10−2 HD-FEC threshold of MIMO-MBNN is from 0.88 to 1.68, which is 8 times and 2.76 times of SISO-LMS and SISO-DNN respectively. In Figs. 9(iv) to 9(vi), when Vpp1 raises to 350 mV, the operation range of MIMO-MBNN is from 0.70 to 1.63, which is 2.74 times and 2.33 times of SISO-LMS and SISO-DNN respectively. Besides, MIMO-MBNN owns the operation range from 0.98 to 1.08 above the 3.8×10−3 HD-FEC threshold, while SISO-LMS and SISO-DNN don’t have. When Vpp1 raises to 600 mV as Figs. 9(vii)–9(ix) depict, the system is distorted by the nonlinear AM response, and the operation range of MIMO-MBNN decreases from 0.78 to 0.92, while SISO-LMS and SISO-DNN own no operation range. Additionally, the average Q performance of SISO-LMS, SISO-DNN, and MIMO-MBNN at (c1), (c2), and (c3) in Figs. 9(iv)–9(vi) are 0.94, 2.77, and 6.12 dB respectively. Therefore, MIMO-MBNN shows a Q improvement of 3.35 dB compared with SISO-LMS and SISO-DNN.

 figure: Fig. 9.

Fig. 9. Q performance of SISO-LMS (left column), SISO-DNN (middle column), and MIMO-MBNN (right column) at different ${\mathrm{\beta}} $. (i)-(iii) Vpp1=250 mV. (iv)-(vi) Vpp1=350 mV. (vi)-(ix) Vpp1=600 mV.

Download Full Size | PDF

At last, we demonstrate the SCM 64QAM constellations and eye-diagrams of Tx1 and Tx2 at the receiver after post equalization of SISO-LMS, SISO-DNN, and MIMO-MBNN. The marks in Fig. 10 correspond to that in Fig. 9. In Figs. 10(a1)–10(a3), the ${\mathrm{\beta}} = 0.39$, which means the power of LED2 at Rx is about 0.4 times of LED1. Therefore, the SNR of Tx1 is bigger than that of Tx2, which could be easily recognized by the eye diagrams and constellations. When the power of Tx1 and Tx2 are almost equal from Figs. 10(b1) to 10(b3), the SCM 64QAM constellations at the receiver shows the shape of the standard 64QAM, and the eye diagrams of LED1 and LED2 behave similarly to each other. As Figs. 10(c1) to 10(c3) shows, increasing ${\mathrm{\beta}} $ to 1.76 will make Tx2 outperform Tx1, which is contrary to Figs. 10(a1) to 10(a3) and accommodates the analysis in section 2. Besides, the eye diagrams and constellations demonstrate that MIMO-MBNN is superior to SISO-LMS and SISO-DNN.

 figure: Fig. 10.

Fig. 10. Constellation and eye diagrams of (a1)-(c3) in Fig. 9.

Download Full Size | PDF

4. Conclusion

In this paper, we proposed for the first time a new structure of neural network post equalizer—MIMO-MBNN in the SR-MIMO VLC system, which modulated by SCM 64QAM constellation and CAP. Experimental results show that the operation range of MIMO-MBNN is much larger than that of SISO-DNN and SISO-LMS above the 2×10−2 HD-FEC threshold. Additionally, MIMO-MBNN can improve the Q performance by at most 3.35 dB compared with SISO-DNN and SISO-LMS. The data rate of 2.1 Gbps above the 3.8×10−3 HD-FEC threshold is achieved by MIMO-MBNN. To the best of our knowledge, this is the highest data rate in the SR-MIMO VLC system, and it is 40% higher than the original record. Further research will focus on the simplification of MIMO-MBNN.

Appendix

Deduction of Eq. (10) is shown as follows. the noise terms ${N_1}(t)$, and ${N_2}(t)$ are ignored:

In this paper, we treat the Fourier transfer pairs as $X(t) \leftrightarrow X({\mathrm{\omega}} )$. Therefore,

$$\begin{array}{l} {X_1}(t) = {S_1}(t) \cdot \cos ({{\mathrm{\omega}} _0}t) \leftrightarrow {X_1}({\mathrm{\omega}} ) = \frac{1}{2}[{S_1}({\mathrm{\omega}} - {{\mathrm{\omega}} _0}) + {S_1}({\mathrm{\omega}} + {{\mathrm{\omega}} _0})]\\ {X_2}(t) ={-} {S_2}(t) \cdot \sin ({{\mathrm{\omega}} _0}t) \leftrightarrow {X_2}({\mathrm{\omega}} ) = \frac{1}{{2i}}[{S_2}({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) - {S_2}({\mathrm{\omega}} - {{\mathrm{\omega}} _0})] \end{array}$$
The ${Y_1}(t)$ at the receiver can be expressed as: ${Y_1}(t) = Y(t) \cdot \cos ({{\mathrm{\omega}} _0}t)$. The frequency response of ${Y_1}(t)$ can be written as:
$$\begin{array}{l} {Y_1}({\mathrm{\omega}} ) = \frac{1}{{2\pi }}[{X_1}({\mathrm{\omega}} ) \cdot H({\mathrm{\omega}} ) + {\mathrm{\beta}} {X_2}({\mathrm{\omega}} ) \cdot H({\mathrm{\omega}} )] \otimes (\pi [\delta ({\mathrm{\omega}} - {{\mathrm{\omega}} _0}) + \delta ({\mathrm{\omega}} + {{\mathrm{\omega}} _0})])\\ {\kern 1pt} {\kern 1pt} \;\;\;\;\;\;\; = \frac{1}{2}[{X_1}({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) \cdot H({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) + {X_1}({\mathrm{\omega}} - {{\mathrm{\omega}} _0}) \cdot H({\mathrm{\omega}} - {{\mathrm{\omega}} _0})] + \\ \;\;\;\;\;\;\;{\kern 1pt} \;\;\;\frac{{\mathrm{\beta}} }{2}[{X_2}({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) \cdot H({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) + {X_2}({\mathrm{\omega}} - {{\mathrm{\omega}} _0}) \cdot H({\mathrm{\omega}} - {{\mathrm{\omega}} _0})] \end{array}$$
While
$$\begin{array}{l} {X_1}({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) = \frac{1}{2}[{S_1}({\mathrm{\omega}} ) + {S_1}({\mathrm{\omega}} + 2{{\mathrm{\omega}} _0})],{X_1}({\mathrm{\omega}} - {{\mathrm{\omega}} _0}) = \frac{1}{2}[{S_1}({\mathrm{\omega}} - 2{{\mathrm{\omega}} _0}) + {S_1}({\mathrm{\omega}} )]\\ {X_2}({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) = \frac{1}{{2i}}[{S_2}({\mathrm{\omega}} + 2{{\mathrm{\omega}} _0}) - {S_2}({\mathrm{\omega}} )],{X_2}({\mathrm{\omega}} - {{\mathrm{\omega}} _0}) = \frac{1}{{2i}}[{S_2}({\mathrm{\omega}} ) - {S_2}({\mathrm{\omega}} - 2{{\mathrm{\omega}} _0})] \end{array}$$
Therefore, the ${Y_1}({\mathrm{\omega}} )$ can be written as:
$$\begin{array}{l} {Y_1}({\mathrm{\omega}} ) = \frac{1}{4}[{S_1}({\mathrm{\omega}} )H({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) + {S_1}({\mathrm{\omega}} + 2{{\mathrm{\omega}} _0})H({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) + {S_1}({\mathrm{\omega}} - 2{{\mathrm{\omega}} _0})H({\mathrm{\omega}} - {{\mathrm{\omega}} _0})\\ \;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} \;\;\; + {S_1}({\mathrm{\omega}} )H({\mathrm{\omega}} - {{\mathrm{\omega}} _0})] + \frac{{\mathrm{\beta}} }{{4i}}[{S_2}({\mathrm{\omega}} + 2{{\mathrm{\omega}} _0})H({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) - {S_2}({\mathrm{\omega}} )H({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) + \\ \;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} \;\;\;{S_2}({\mathrm{\omega}} )H({\mathrm{\omega}} - {{\mathrm{\omega}} _0}) - {S_2}({\mathrm{\omega}} - 2{{\mathrm{\omega}} _0})H({\mathrm{\omega}} - {{\mathrm{\omega}} _0})]\\ \;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} = \frac{1}{4}[{S_1}({\mathrm{\omega}} )(H({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) + H({\mathrm{\omega}} - {{\mathrm{\omega}} _0}))] + \frac{1}{4}{S_1}({\mathrm{\omega}} + 2{{\mathrm{\omega}} _0})H({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) + \\ \;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} \;\;\;\frac{1}{4}{S_1}({\mathrm{\omega}} - 2{{\mathrm{\omega}} _0})H({\mathrm{\omega}} - {{\mathrm{\omega}} _0}) + \frac{{\mathrm{\beta}} }{{4i}}[{S_2}({\mathrm{\omega}} )(H({\mathrm{\omega}} - {{\mathrm{\omega}} _0}) - H({\mathrm{\omega}} + {{\mathrm{\omega}} _0}))] + \\ \;\;\;\;\;\;\;{\kern 1pt} {\kern 1pt} \;\;\;\frac{{\mathrm{\beta}} }{{4i}}{S_2}({\mathrm{\omega}} + 2{{\mathrm{\omega}} _0})H({\mathrm{\omega}} + {{\mathrm{\omega}} _0}) - \frac{{\mathrm{\beta}} }{{4i}}{S_2}({\mathrm{\omega}} - 2{{\mathrm{\omega}} _0})H({\mathrm{\omega}} - {{\mathrm{\omega}} _0}) \end{array}$$
Then the inverse Fourier transform can be expressed as:
$$\begin{array}{l} {Y_1}(t) = \frac{1}{4}[{S_1}(t) \otimes (h(t) \cdot 2\cos ({{\mathrm{\omega}} _0}t))] + \frac{1}{4}({S_1}(t){e^{ - i2{{\mathrm{\omega}} _0}t}}) \otimes (h(t){e^{ - i{{\mathrm{\omega}} _0}t}})\\ \;\;\;\;\;\;\;\;\; + \frac{1}{4}({S_1}(t){e^{i2{{\mathrm{\omega}} _0}t}}) \otimes (h(t){e^{i{{\mathrm{\omega}} _0}t}}) + \frac{{\mathrm{\beta}} }{{4i}}[{S_2}(t) \otimes (h(t) \cdot 2i \cdot \sin ({{\mathrm{\omega}} _0}t))]\\ \;\;\;\;\;\;\;\;\; + \frac{{\mathrm{\beta}} }{{4i}}({S_2}(t){e^{ - i2{{\mathrm{\omega}} _0}t}}) \otimes (h(t){e^{ - i{{\mathrm{\omega}} _0}t}}) - \frac{{\mathrm{\beta}} }{{4i}}({S_2}(t){e^{i2{{\mathrm{\omega}} _0}t}}) \otimes (h(t){e^{i{{\mathrm{\omega}} _0}t}})\\ \;\;\;\;\;\;\; = \frac{1}{2}[{S_1}(t) \otimes (h(t) \cdot \cos ({{\mathrm{\omega}} _0}t))] + \frac{{\mathrm{\beta}} }{2}[{S_2}(t) \otimes (h(t) \cdot \sin ({{\mathrm{\omega}} _0}t))]\\ \;\;\;\;\;\;\;\;\; + \frac{1}{4}({S_1}(t)(\cos (2{{\mathrm{\omega}} _0}t) - i \cdot \sin (2{{\mathrm{\omega}} _0}t))) \otimes (h(t)(\cos ({{\mathrm{\omega}} _0}t) - i \cdot \sin ({{\mathrm{\omega}} _0}t)))\\ \;\;\;\;\;\;\;\;\; + \frac{1}{4}({S_1}(t)(\cos (2{{\mathrm{\omega}} _0}t) + i \cdot \sin (2{{\mathrm{\omega}} _0}t))) \otimes (h(t)(\cos ({{\mathrm{\omega}} _0}t) + i \cdot \sin ({{\mathrm{\omega}} _0}t)))\\ \;\;\;\;\;\;\;\;\; + \frac{{\mathrm{\beta}} }{{4i}}({S_2}(t)(\cos (2{{\mathrm{\omega}} _0}t) - i \cdot \sin (2{{\mathrm{\omega}} _0}t))) \otimes (h(t)(\cos ({{\mathrm{\omega}} _0}t) - i \cdot \sin ({{\mathrm{\omega}} _0}t)))\\ \;\;\;\;\;\;\;\;\; - \frac{{\mathrm{\beta}} }{{4i}}({S_2}(t)(\cos (2{{\mathrm{\omega}} _0}t) + i \cdot \sin (2{{\mathrm{\omega}} _0}t))) \otimes (h(t)(\cos ({{\mathrm{\omega}} _0}t) + i \cdot \sin ({{\mathrm{\omega}} _0}t)))\\ \;\;\;\;\;\;\; = \frac{1}{2}[{S_1}(t) \otimes (h(t) \cdot \cos ({{\mathrm{\omega}} _0}t))] + \frac{{\mathrm{\beta}} }{2}[{S_2}(t) \otimes (h(t) \cdot \sin ({{\mathrm{\omega}} _0}t))]\\ \;\;\;\;\;\;\;\;\; + \frac{1}{2}{S_1}(t)\cos (2{{\mathrm{\omega}} _0}t) \otimes h(t)\cos ({{\mathrm{\omega}} _0}t) - \frac{1}{2}{S_1}(t)\sin (2{{\mathrm{\omega}} _0}t) \otimes h(t)\sin ({{\mathrm{\omega}} _0}t)\\ \;\;\;\;\;\;\;\;\; - \frac{{\mathrm{\beta}} }{2}{S_2}(t)\cos (2{{\mathrm{\omega}} _0}t) \otimes h(t)\sin ({{\mathrm{\omega}} _0}t) - \frac{{\mathrm{\beta}} }{2}{S_2}(t)\sin (2{{\mathrm{\omega}} _0}t) \otimes h(t)\cos ({{\mathrm{\omega}} _0}t)\\ \;\;\;\;\;\;\; = \frac{1}{2}[{S_1}(t) \otimes (h(t) \cdot \cos ({{\mathrm{\omega}} _0}t))] + \frac{{\mathrm{\beta}} }{2}[{S_2}(t) \otimes (h(t) \cdot \sin ({{\mathrm{\omega}} _0}t))]\\ \;\;\;\;\;\;\;\;\; - \frac{1}{2}\sqrt {{S_1}{{(t)}^2} + {{({\mathrm{\beta}} {S_2}(t))}^2}} \sin (2{{\mathrm{\omega}} _0}t + {{\mathrm{\alpha}} _1}) \otimes (h(t)\cos ({{\mathrm{\omega}} _0}t))\\ \;\;\;\;\;\;\;\;\; - \frac{1}{2}\sqrt {{S_1}{{(t)}^2} + {{({\mathrm{\beta}} {S_2}(t))}^2}} \sin (2{{\mathrm{\omega}} _0}t + {{\mathrm{\alpha}} _2}) \otimes (h(t)\sin ({{\mathrm{\omega}} _0}t)) \end{array}$$
Here $\tan {{\mathrm{\alpha}} _1} ={-} \frac{{{S_1}(t)}}{{{\mathrm{\beta}} {S_2}(t)}},\tan {{\mathrm{\alpha}} _2} = \frac{{{\mathrm{\beta}} {S_2}(t)}}{{{S_1}(t)}}$. The deduction of Eq. (11) is similar to that of Eq. (10).

Funding

National Key Research and Development Program of China (2017YFB0403603); National Natural Science Foundation of China (No.61925104).

Acknowledgments

This work was partially supported by the National Key Research and Development Program of China (2017YFB0403603), the NSFC project (No.61925104).

Disclosures

The authors declare no conflicts of interest.

References

1. S. Wu, H. Wang, and C. Youn, “Visible light communications for 5G wireless networking systems: from fixed to mobile communications,” IEEE Network 28(6), 41–45 (2014). [CrossRef]  

2. N. Chi, H. Haas, M. Kavehrad, T. Little, and X. Huang, “Visible Light Communications:Demand Factors, Benefits and Opportunities,” IEEE Wireless Commun. 22(2), 5–7 (2015). [CrossRef]  

3. 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]  

4. N. Chi, Y. Zhou, S. Liang, F. Wang, J. Li, and Y. Wang, “Enabling Technologies for High-Speed Visible Light Communication Employing CAP Modulation,” J. Lightwave Technol. 36(2), 510–518 (2018). [CrossRef]  

5. C. Kottke, S. Nerreter, K. Langer, and J. W. Walewski, “513 Mbit / s Visible Light Communications Link Based on DMT-Modulation of a White LED,” J. Lightwave Technol. 28(24), 3512–3518 (2010). [CrossRef]  

6. X. Huang, Z. Wang, J. Shi, Y. Wang, and N. Chi, “1 . 6 Gbit / s phosphorescent white LED based VLC transmission using a cascaded pre-equalization circuit and a differential outputs PIN receiver,” Opt. Express 23(17), 22034–22042 (2015). [CrossRef]  

7. X. Huang, J. Shi, J. Li, Y. Wang, Y. Wang, and N. Chi, “750Mbit/s Visible Light Communications employing 64QAM-OFDM Based on Amplitude Equalization Circuit,” in 2015 Optical Fiber Communications Conference and Exhibition (OFC) (2015), pp. 15–17.

8. K. Ying, H. Qian, R. J. Baxley, and S. Yao, “Joint Optimization of Precoder and Equalizer in MIMO VLC Systems,” IEEE J. Select. Areas Commun. 33(9), 1949–1958 (2015). [CrossRef]  

9. H. Burchardt, N. Serafimovski, D. Tsonev, S. Videv, and H. Haas, “VLC: Beyond Point-to-Point Communication,” IEEE Commun. Mag. 52(7), 98–105 (2014). [CrossRef]  

10. Y. Wang and N. Chi, “Demonstration of High-Speed 2 × 2 Non-Imaging MIMO Nyquist Single Carrier Visible Light Communication With Frequency Domain Equalization,” J. Lightwave Technol. 32(11), 2087–2093 (2014). [CrossRef]  

11. L. Zeng, D. C. O’Brien, H. Le Minh, G. E. Faulkner, K. Lee, D. Jung, Y. Oh, and E. T. Won, “High data rate Multiple Input Multiple Output (MIMO) optical wireless communications using white LED lighting,” IEEE J. Select. Areas Commun. 27(9), 1654–1662 (2009). [CrossRef]  

12. Y. Hong, J. Xu, and L.-K. Chen, “Experimental investigation of multi-band OCT precoding for OFDM-based visible light communications,” Opt. Express 25(11), 12908 (2017). [CrossRef]  

13. Y. Hong, L. K. Chen, and J. Zhao, “Experimental demonstration of performance-enhanced MIMO-OFDM visible light communications,” 2017 Opt. Fiber Commun. Conf. Exhib. OFC 2017 - Proc. (1), 2–4 (2017).

14. M. Shi, M. Zhang, F. Wang, M. Zhao, and N. Chi, “Equiprobable pre-coding PAM7 Modulation for nonlinearity mitigation in underwater 2 × 1 MISO visible light communications,” J. Lightwave Technol. 36(22), 5188–5195 (2018). [CrossRef]  

15. L. Qiao, X. Lu, S. Liang, J. Zhang, and N. Chi, “Performance analysis of space multiplexing by superposed signal in multi-dimensional VLC system,” Opt. Express 26(16), 19762 (2018). [CrossRef]  

16. X. Guo and N. Chi, “Superposed 32QAM Constellation Design for 2 × 2 Spatial Multiplexing MIMO VLC Systems,” J. Lightwave Technol. 38(7), 1702–1711 (2020). [CrossRef]  

17. K. Ying, Z. Yu, R. J. Baxley, H. Qian, G.-K. Chang, and G. T. Zhou, “Nonlinear distortion mitigation in visible light communications,” IEEE Wireless Commun. 22(2), 36–45 (2015). [CrossRef]  

18. F. Miramirkhani and M. Uysal, “Visible Light Communication Channel Modeling for Underwater Environments with Blocking and Shadowing,” IEEE Access 6, 1082–1090 (2018). [CrossRef]  

19. N. Chi, Y. Zhao, M. Shi, P. Zou, and X. Lu, “Gaussian kernel-aided deep neural network equalizer utilized in underwater PAM8 visible light communication system,” Opt. Express 26(20), 26700 (2018). [CrossRef]  

20. Y. Zhao, P. Zou, W. Yu, and N. Chi, “Two tributaries heterogeneous neural network based channel emulator for underwater visible light communication systems,” Opt. Express 27(16), 22532 (2019). [CrossRef]  

21. N. Chi, M. Zhang, J. Shi, and Y. Zhao, “Spectrally efficient multi-band visible light communication system based on Nyquist PAM-8 modulation,” Photonics Res. 5(6), 588 (2017). [CrossRef]  

22. N. Chi and M. Shi, “Advanced modulation formats for underwater visible light communications [Invited],” Chin. Opt. Lett. 16(12), 120603 (2018). [CrossRef]  

23. Y. Wang and N. Chi, “8-Gb / s RGBY LED-Based WDM VLC System Employing High-Order CAP Modulation and Hybrid Post Equalizer System Employing High-Order CAP Modulation and Hybrid Post Equalizer,” IEEE Photonics J. 7(6), 1–7 (2015). [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 (10)

Fig. 1.
Fig. 1. Principles of SCM 64QAM schemes. (a) PAM8. (b) QPSK and 16QAM in [15].
Fig. 2.
Fig. 2. Data flow and experimental setup of the SR-MIMO CAP VLC system. AWG: arbitrary waveform generator; Eq: hardware equalizer; EA: electrical amplifier; OSC: oscilloscope.
Fig. 3.
Fig. 3. Structure of MIMO-MBNN
Fig. 4.
Fig. 4. amplitude magnitude response. (a) Tx1 and Rx without Tx2 at different Vpp. (b) Tx1 and Rx with varying values of Tx1: Tx2 (Vpp1=900 mV)
Fig. 5.
Fig. 5. BER performance versus different node numbers and activation function of MIMO-MBNN. (a) BER performance of 1 hidden layer. (b) BER performance of 2 hidden layers
Fig. 6.
Fig. 6. Frequency and spectrum mismatch of the received signal, SISO-LMS, SISO-DNN, and MIMO-MBNN. (a) Frequency response. (b) spectrum mismatch
Fig. 7.
Fig. 7. Q factor versus data rate
Fig. 8.
Fig. 8. Contour of system performance SISO-LMS, SISO-DNN, and MIMO-MBNN (from 1 to 3 in each row). (a1)-(a3). Average BER performance. (b1)-(b3). Tx1 BER performance. (c1)-(c3). Tx2 BER performance.
Fig. 9.
Fig. 9. Q performance of SISO-LMS (left column), SISO-DNN (middle column), and MIMO-MBNN (right column) at different ${\mathrm{\beta}} $. (i)-(iii) Vpp1=250 mV. (iv)-(vi) Vpp1=350 mV. (vi)-(ix) Vpp1=600 mV.
Fig. 10.
Fig. 10. Constellation and eye diagrams of (a1)-(c3) in Fig. 9.

Tables (2)

Tables Icon

Table 1. Details of components in Fig. 2.

Tables Icon

Table 2. Detailed parameters of SISO-LMS, SISO-DNN, and MIMO-MBNN.

Equations (27)

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

T = [ T 1 ( t ) , T 2 ( t ) ] , T k ( t ) { ± 1 , ± 3 , ± 5 , ± 7 } .
Y ( t ) = Y 1 ( t ) + β Y 2 ( t ) i .
[ S 1 ( t ) S 2 ( t ) ] = [ g 1 ( t ) T 1 u p ( t ) g 2 ( t ) T 2 u p ( t ) ] .
[ X 1 ( t ) X 2 ( t ) ] = [ cos ( 2 π f 0 t ) 0 0 sin ( 2 π f 0 t ) ] [ S 1 ( t ) S 2 ( t ) ] .
[ R 1 ( t ) R 2 ( t ) ] = [ H 11 H 12 H 21 H 22 ] [ f 1 ( X 1 ( t ) ) f 2 ( X 2 ( t ) ) ]  +  [ N 1 ( t ) N 2 ( t ) ] .
f ( x ) = { x ( a 1 x + a 2 x 2 + a 3 x 3 ) / b , l o w V p p , h i g h V p p .
[ R 1 ( t ) R 2 ( t ) ] = [ 1 0 0 β ] [ X 1 ( t ) h 1 ( t ) X 2 ( t ) h 2 ( t ) ]  +  [ N 1 ( t ) N 2 ( t ) ] .
Y ( t ) = f ( R 1 ( t ) + R 2 ( t ) ) + N ( t )  =  f ( X 1 ( t ) h ( t ) + β X 2 ( t ) h ( t ) ) + N ( t ) .
[ Y 1 ( t ) Y 2 ( t ) ] = [ cos ( 2 π f 0 t ) sin ( 2 π f 0 t ) ] Y ( t ) = [ ( ( S 1 ( t ) cos ( ω 0 t ) ) h ( t ) + β ( S 2 ( t ) sin ( ω 0 t ) ) h ( t ) ) cos ( ω 0 t ) ( ( S 1 ( t ) cos ( ω 0 t ) ) h ( t ) + β ( S 2 ( t ) sin ( ω 0 t ) ) h ( t ) ) sin ( ω 0 t ) ] + [ N 1 ( t ) N 2 ( t ) ] .
Y 1 ( t ) = 1 2 S 1 ( t ) ( h ( t ) cos ( ω 0 t ) ) + β 2 S 2 ( t ) ( h ( t ) sin ( ω 0 t ) ) 1 2 S 1 ( t ) 2 + ( β S 2 ( t ) ) 2 sin ( 2 ω 0 t + α 1 ) ( h ( t ) cos ( ω 0 t ) ) 1 2 S 1 ( t ) 2 + ( β S 2 ( t ) ) 2 sin ( 2 ω 0 t + α 2 ) ( h ( t ) sin ( ω 0 t ) ) + N 1 ( t ) .
Y 2 ( t ) = 1 2 S 1 ( t ) ( h ( t ) sin ( ω 0 t ) ) + β 2 S 2 ( t ) ( h ( t ) cos ( ω 0 t ) ) 1 2 S 1 ( t ) 2 + ( β S 2 ( t ) ) 2 sin ( 2 ω 0 t + α 2 ) ( h ( t ) cos ( ω 0 t ) ) + 1 2 S 1 ( t ) 2 + ( β S 2 ( t ) ) 2 sin ( 2 ω 0 t + α 1 ) ( h ( t ) sin ( ω 0 t ) ) + N 2 ( t ) .
Y 1 ( t ) = F 1 ( S 1 ( t ) , S 2 ( t ) ) Y 2 ( t ) = F 2 ( S 2 ( t ) , S 1 ( t ) )
Y 1 ( t ) = a 1 b F 1 ( S 1 ( t ) , S 2 ( t ) ) + a 2 b F 1 ( S 1 2 ( t ) , S 2 2 ( t ) ) + a 3 b F 1 ( S 1 3 ( t ) , S 2 3 ( t ) ) Y 2 ( t ) = a 1 b F 2 ( S 2 ( t ) , S 1 ( t ) ) + a 2 b F 2 ( S 2 2 ( t ) , S 1 2 ( t ) ) + a 3 b F 2 ( S 2 3 ( t ) , S 1 3 ( t ) )
Y 1 ( t ) = a 1 b ( 1 2 S 1 ( t ) ( h ( t ) cos ( ω 0 t ) ) l i n e a r d i s t o r t i o n + β 2 S 2 ( t ) ( h ( t ) sin ( ω 0 t ) ) l i n e a r c r o s s t a l k + a 2 b ( 1 2 S 1 2 ( t ) ( h ( t ) cos ( ω 0 t ) ) + β 2 S 2 2 ( t ) ( h ( t ) sin ( ω 0 t ) ) n o n l i n e a r d i s t o r t i o n a n d c r o s s t a l k + a 3 b ( 1 2 S 1 3 ( t ) ( h ( t ) cos ( ω 0 t ) ) + β 2 S 2 3 ( t ) ( h ( t ) sin ( ω 0 t ) ) n o n l i n e a r d i s t o r t i o n a n d c r o s s t a l k
F max ( a b s ) ( x ) = x / max ( a b s ( x ) ) .
{ Y 1 , L ( t ) = W 1 ( 1 ) T Y ( t ) + b 1 ( 1 ) [ Y 1 , N L ( t ) Y 2 , N L ( t ) ] = W 3 ( 3 ) T R ( W 3 ( 2 ) T R ( W 3 ( 1 ) T Y ( t ) + b 3 ( 1 ) ) + b 3 ( 2 ) ) + b 3 ( 3 ) Y 2 , L ( t ) = W 2 ( 1 ) T Y ( t ) + b 2 ( 1 )
R ( x ) = { x , x > 0 0 , x 0 .
tanh ( x ) = e x e x e x + e x .
P 1 ( t ) = Y 1 , L ( t ) + Y 1 , N L ( t ) , P 2 ( t ) = Y 2 , L ( t ) + Y 2 , N L ( t ) .
W k ( n ) , b k ( n ) = arg min W k ( n ) , b k ( n ) 1 N t = 2 N 1 | | P 1 ( t ) S 1 ( t ) | | 2 + | | P 2 ( t ) S 2 ( t ) | | 2 .
M ( ω ) = | X ( ω ) P ( ω ) | 2
Q ( d B ) = 20 log 10 ( 2 e r f c 1 ( 2 B E R ) )
X 1 ( t ) = S 1 ( t ) cos ( ω 0 t ) X 1 ( ω ) = 1 2 [ S 1 ( ω ω 0 ) + S 1 ( ω + ω 0 ) ] X 2 ( t ) = S 2 ( t ) sin ( ω 0 t ) X 2 ( ω ) = 1 2 i [ S 2 ( ω + ω 0 ) S 2 ( ω ω 0 ) ]
Y 1 ( ω ) = 1 2 π [ X 1 ( ω ) H ( ω ) + β X 2 ( ω ) H ( ω ) ] ( π [ δ ( ω ω 0 ) + δ ( ω + ω 0 ) ] ) = 1 2 [ X 1 ( ω + ω 0 ) H ( ω + ω 0 ) + X 1 ( ω ω 0 ) H ( ω ω 0 ) ] + β 2 [ X 2 ( ω + ω 0 ) H ( ω + ω 0 ) + X 2 ( ω ω 0 ) H ( ω ω 0 ) ]
X 1 ( ω + ω 0 ) = 1 2 [ S 1 ( ω ) + S 1 ( ω + 2 ω 0 ) ] , X 1 ( ω ω 0 ) = 1 2 [ S 1 ( ω 2 ω 0 ) + S 1 ( ω ) ] X 2 ( ω + ω 0 ) = 1 2 i [ S 2 ( ω + 2 ω 0 ) S 2 ( ω ) ] , X 2 ( ω ω 0 ) = 1 2 i [ S 2 ( ω ) S 2 ( ω 2 ω 0 ) ]
Y 1 ( ω ) = 1 4 [ S 1 ( ω ) H ( ω + ω 0 ) + S 1 ( ω + 2 ω 0 ) H ( ω + ω 0 ) + S 1 ( ω 2 ω 0 ) H ( ω ω 0 ) + S 1 ( ω ) H ( ω ω 0 ) ] + β 4 i [ S 2 ( ω + 2 ω 0 ) H ( ω + ω 0 ) S 2 ( ω ) H ( ω + ω 0 ) + S 2 ( ω ) H ( ω ω 0 ) S 2 ( ω 2 ω 0 ) H ( ω ω 0 ) ] = 1 4 [ S 1 ( ω ) ( H ( ω + ω 0 ) + H ( ω ω 0 ) ) ] + 1 4 S 1 ( ω + 2 ω 0 ) H ( ω + ω 0 ) + 1 4 S 1 ( ω 2 ω 0 ) H ( ω ω 0 ) + β 4 i [ S 2 ( ω ) ( H ( ω ω 0 ) H ( ω + ω 0 ) ) ] + β 4 i S 2 ( ω + 2 ω 0 ) H ( ω + ω 0 ) β 4 i S 2 ( ω 2 ω 0 ) H ( ω ω 0 )
Y 1 ( t ) = 1 4 [ S 1 ( t ) ( h ( t ) 2 cos ( ω 0 t ) ) ] + 1 4 ( S 1 ( t ) e i 2 ω 0 t ) ( h ( t ) e i ω 0 t ) + 1 4 ( S 1 ( t ) e i 2 ω 0 t ) ( h ( t ) e i ω 0 t ) + β 4 i [ S 2 ( t ) ( h ( t ) 2 i sin ( ω 0 t ) ) ] + β 4 i ( S 2 ( t ) e i 2 ω 0 t ) ( h ( t ) e i ω 0 t ) β 4 i ( S 2 ( t ) e i 2 ω 0 t ) ( h ( t ) e i ω 0 t ) = 1 2 [ S 1 ( t ) ( h ( t ) cos ( ω 0 t ) ) ] + β 2 [ S 2 ( t ) ( h ( t ) sin ( ω 0 t ) ) ] + 1 4 ( S 1 ( t ) ( cos ( 2 ω 0 t ) i sin ( 2 ω 0 t ) ) ) ( h ( t ) ( cos ( ω 0 t ) i sin ( ω 0 t ) ) ) + 1 4 ( S 1 ( t ) ( cos ( 2 ω 0 t ) + i sin ( 2 ω 0 t ) ) ) ( h ( t ) ( cos ( ω 0 t ) + i sin ( ω 0 t ) ) ) + β 4 i ( S 2 ( t ) ( cos ( 2 ω 0 t ) i sin ( 2 ω 0 t ) ) ) ( h ( t ) ( cos ( ω 0 t ) i sin ( ω 0 t ) ) ) β 4 i ( S 2 ( t ) ( cos ( 2 ω 0 t ) + i sin ( 2 ω 0 t ) ) ) ( h ( t ) ( cos ( ω 0 t ) + i sin ( ω 0 t ) ) ) = 1 2 [ S 1 ( t ) ( h ( t ) cos ( ω 0 t ) ) ] + β 2 [ S 2 ( t ) ( h ( t ) sin ( ω 0 t ) ) ] + 1 2 S 1 ( t ) cos ( 2 ω 0 t ) h ( t ) cos ( ω 0 t ) 1 2 S 1 ( t ) sin ( 2 ω 0 t ) h ( t ) sin ( ω 0 t ) β 2 S 2 ( t ) cos ( 2 ω 0 t ) h ( t ) sin ( ω 0 t ) β 2 S 2 ( t ) sin ( 2 ω 0 t ) h ( t ) cos ( ω 0 t ) = 1 2 [ S 1 ( t ) ( h ( t ) cos ( ω 0 t ) ) ] + β 2 [ S 2 ( t ) ( h ( t ) sin ( ω 0 t ) ) ] 1 2 S 1 ( t ) 2 + ( β S 2 ( t ) ) 2 sin ( 2 ω 0 t + α 1 ) ( h ( t ) cos ( ω 0 t ) ) 1 2 S 1 ( t ) 2 + ( β S 2 ( t ) ) 2 sin ( 2 ω 0 t + α 2 ) ( h ( t ) sin ( ω 0 t ) )
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.