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

Time delay estimation from the time series for optical chaos systems using deep learning

Open Access Open Access

Abstract

We propose a model-free time delay signature (TDS) extraction method for optical chaos systems. The TDS can be identified from time series without prior knowledge of the actual physical processes. In optical chaos secure communication systems, the chaos carrier is usually generated by a laser diode subject to opto-electronic/all-optical time delayed feedback. One of the most important factors to security considerations is the concealment of the TDS. So far, statistical analysis methods such as autocorrelation function (ACF) and delayed mutual information (DMI) are usually used to unveil the TDS. However, the effectiveness of these methods will be reduced when increasing the nonlinearity of chaos systems. Meanwhile, certain TDS concealment strategies have been designed against statistical analysis. In our previous work, convolutional neural network shows its effectiveness on TDS extraction of chaos systems with high loop nonlinearity. However, this method relies on the knowledge of detailed structure of the chaos systems. In this work, we formulate a blind identification method based on long short-term memory neural network (LSTM-NN) model. The method is validated against the two major types of optical chaos systems, i.e. opto-electronic oscillator (OEO) chaos system and laser chaos system based on internal nonlinearity. Moreover, some security enhanced chaotic systems are also studied. The results show that the proposed method has high tolerance to additive noise. Meanwhile, the data amount needed is less than existing methods.

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

1. Introduction

Optical or electro-optical time-delayed chaos system has unique advantages of high dimension, high dynamical complexity, and especially remarkable broadband spectrum, which offers the potential for data secure transmission at high rates. Since the early demonstrations of optical chaos secure communication (OCSC) [1], OCSC has been widely studied and become one of the most feasible hardware encryption methods [27]. In OCSC, without the matched devices or parameters settings, it is hard for an eavesdropper to recover the message embedded in chaotic carriers. At present, there are two major types of optical chaos system. One is built on the internal nonlinearity of laser [8]. Semiconductor lasers subjected to external perturbation, such as all-optical time-delayed feedback and external optical injection, can generate chaos. The other is called opto-electronic oscillator (OEO) chaos system [9,10], where chaos is generated through an external nonlinear opto-electronic time-delayed feedback. Generally, these systems can be characterized by a time delayed differential equation:

$$\dot{x}(t )= f({x(t ),\; x({t - \tau } )} ), $$
where $x(t )$ is the state variable, and $x({t - \tau } )$ is the time delayed variation. The loop delay parameter $\tau \; $(i.e., time delay signature, TDS) is usually chosen as a cryptographic key. Because the delay value can be varied over a wide range by changing the loop length of the system circuit. Meanwhile the delay values at the transmitter and receiver sides can be matched accurately using tunable optical delay lines. As a result, a large time domain key space can be obtained. Therefore, an important factor to security considerations is the concealment of TDS.

So far, much research has focused on concealing the TDS [1115]. However, it will take some time for its security analysis to mature. For a successful recovery of the TDS, one can use different strategies, such as permutation information analysis [16,17], RF spectrum analyzer [18,19], singular values fraction measure [20], optimal transformations [21], autocorrelation function (ACF), delayed mutual information (DMI) [14,15,22], local linear fitting in a low-dimensional subspace [23], extrema statistics [24], and so on. Out of the methods listed, ACF and DMI are most widely used due to their effectiveness and robustness [25]. If there is a certain degree of linear correlation between a chaotic signal $x(t )$ and its time delayed variation $x({t - s} )$ when s equals to the system loop time delay $\tau $, a sharp peak located at $s = \tau $ will exhibit in the ACF and DMI computed under variation of s. But the linear correlation will be reduced or even eliminated when increasing the nonlinearity of the loop transformation, then the peak will decrease or even disappear [26].

Recently, deep learning technique has attracted more and more attention due to its advantages in processing complex nonlinear time series analysis and nonlinear modeling [27,28]. In [29], neural network (NN) was explored for travel time prediction. Prediction of time delay chaos systems via Long Short-Term Memory neural network (LSTM-NN) was studied in [30]. Thusly deep learning is an attractive technology aiming at security analysis of OCSC systems. In [31], a TDS extraction method was proposed based on convolutional neural network (CNN). The method is effective for time-delayed systems with high loop nonlinearity, and the needed minimal data size is less than ACF and DMI. Meanwhile, the scheme shows a good robustness in terms of handling noise. Unfortunately, this kind of analysis method is only effective under the condition that the structure of the OCSC is completely public. This precondition coincides with the “Kerckhoff’s principle” in cryptanalysis. From the viewpoint of implementing complexity, repeatedly offline training of different CNN is essential for different chaos systems, and the training accuracy depends on the approximation degree between offline system and actual system. These factors reduce the applicability of the scheme.

Hence, chaotic communication cryptanalysis against the TDS is still a challenging open problem. In this work, we propose a data-driven TDS extraction method, with which any prior knowledge of the actual physical processes is no longer necessary. A LSTM-NN is employed to carry out blind inverse modeling. From experimental and simulated time series, successful TDS extraction is given to show the effectiveness of our strategy.

2. Methodology

In Eq. (1), there will be an exact functional relationship between the history and future, which can be expressed as:

$$\varphi :\;x(t )\to {\; }x({t + \tau } ),$$
and the relationship is established by the nonlinear function f and the time delay constant $\tau $. Assume that we can observe a time series $\{{{x_i}} \}$ from system (1), as shown in Fig. 1. We arbitrarily select two subsequences $\{{{x_i}} \}_n^m$ and $\{{{x_i}} \}_{n + s}^{m + s}$. Subsequence $\{{{x_i}} \}_n^m$ is used as the input of a NN, and $\{{{x_i}} \}_{n + s}^{m + s}$ is set as the expected output. Then the NN could be trained to finish a sequence-to-sequence regression, as described by Eq. (3).
$${\varphi _s}:\;\{{{x_i}} \}_n^m \to {\; }\{{{x_i}} \}_{n + s}^{m + s}. $$

 figure: Fig. 1.

Fig. 1. The principle of the proposed strategy.

Download Full Size | PDF

As a result, an approximation model of $\varphi $ could be obtained. Here ${\varphi _s}$ will approach to $\varphi $ only when $s = \tau $. Subsequently, the approximation function ${\varphi _{s = \tau }}$ will be able to predict any subsequence $\{{{x_i}} \}_{p + s}^{q + s}$ from its historical information $\{{{x_i}} \}_p^q$.

By contrast, if $s \ne \tau $, the fixed function will no longer exist. The relationship between $\{{{x_i}} \}_n^m$ and $\{{{x_i}} \}_{n + s}^{m + s}$ could be considered as random mapping, which will lead to the failure of prediction ability of model ${\varphi _s}$.

Based on these facts, a blind identification strategy for TDS is shown in Fig. 2. Without knowing the actual time delay parameter $\tau $, we assume some values of time delay candidate $s \in \{{{s_1}, \ldots ,{s_N}} \}$. The training process is shown in Fig. 2(a). For $s = {s_j}$ ($j = 1, \ldots ,N$), subsequence $\{{{x_i}} \}_n^m$ is used as the input of LSTM-NN, and $\{{{x_i}} \}_{n + {s_j}}^{m + {\; }{s_j}}$ is set as the expected output. Then the NN is trained to get a model ${\varphi _j}$. In the testing/predicting process as shown in Fig. 2(b), we put another subsequence $\{{{x_i}} \}_p^q$ into the trained NN model ${\varphi _j}$. A predicted sequence noted as $\{{{y_i}} \}_p^q$ will be obtained. The prediction reliability of model ${\varphi _j}$ is evaluated using normalized root mean square error (NRMSE) [30], which is given by Eq. (4).

$$E(s )= \frac{{\mathop \sum \nolimits_{i = p}^q {{({{x_{i + {s_j}}} - {y_i}} )}^2}}}{{\mathop \sum \nolimits_{i = p}^q {{({{x_{i + {s_j}}} - \bar{x}} )}^2}}}\quad \textrm{for}\;s = {s_j}$$
Here $\{{{x_i}} \}_{p + {s_j}}^{q + {s_j}}$ is the expected output of historical sequence $\{{{x_i}} \}_p^q$, and $\bar{x}$ is the mean value of $\{{{x_i}} \}_{p + {s_j}}^{q + {s_j}}$. For a time delay candidate set $\{{{s_1}, \ldots ,{s_N}} \}$, we can get a model set {$ {{\varphi_j}} |s = {s_j}$}. Subsequently, a NRMSE curve $E(s )$ can be obtained. As shown in Fig. 2(b), if we plot $E(s )$ versus s, there will be a minimum value at $s = \tau $, which means that the actual delay of system (1) is recovered.

 figure: Fig. 2.

Fig. 2. Data driven self-modelling TDS extraction via NN

Download Full Size | PDF

Based on the strategy, the detailed implementation method is stated as follows. Equation (1) can be approximately written as the following difference equation,

$$x({t + \Delta t} )- x(t )= \Delta tf({x(t ),x({t - \tau } )} ), $$
where $\Delta t$ is the unit time step. In order to obtain enough training set from the time series $\{{{x_i}} \}$ with limited length, the training set are formatted as
$$XTrai{n_l} = \{{{x_i}} \}_{1 + l\delta }^{L + l\delta }, $$
$$YTrai{n_l} = \{{{x_i}} \}_{1 + s + l\delta }^{L + s + l\delta } - \{{{x_i}} \}_{s + l\delta }^{L + s + l\delta - 1},$$
where $XTrai{n_l}$ and $YTrai{n_l}$ are data pairs with length L. These training data pairs moves on the time axis with a step of $\delta $. $l = 1,2, \ldots ,{\; }{M_1}$, where ${M_1}$ is the number of all the data pairs. ${\varphi _s}\; $can be obtained from these training set for each s. Similarly, the forecasting set are formatted as
$$XTes{t_l} = \{{{x_i}} \}_{L + {M_1}\delta + 1 + l\vartheta }^{L^{\prime} + L + {M_1}\delta + 1 + l\vartheta },$$
$$YTes{t_l} = \{{{x_i}} \}_{L + {M_1}\delta + 1 + l\vartheta + s}^{L^{\prime} + L + {M_1}\delta + 1 + l\vartheta + s},\; $$
where $l = 1,{\; } \ldots ,{\; }{M_2}$, and ${M_2}$ is the number of forecasting data pairs $XTes{t_l}$ and $YTes{t_l}$, $L^{\prime}$ is the data length. The forecasting data pairs moves on the time axis with a step of $\vartheta $. Put $XTest$ into the trained model ${\varphi _s}$ to get a corresponding predicted sequence $YPre$. Using $YPre$ and $YTest$, the NRMSE will be calculated according to Eq. (4).

To finish the sequence-to-sequence regression, we employ the LSTM-NN because its unique advantages in processing one dimensional time series with strong nonlinearity [30,32,33]. MATLAB2018a Neural Network Toolbox is used to carry out modelling and testing/ forecasting work in our study. A LSTM-NN is trained with the following layer specifications. The number of hidden units is set as 200, then a 50 size fully connected layer and a drop layer with a drop probability of 0.5 are created. The solver ‘Adam’ is used to train with mini-batch size of 20. Initial learning rate is set as 0.01. The gradient threshold is set as 1. The number of training iterations (MaxEpochs) is set as $\mathrm{\beth}$ given for different system.

3. TDS extraction from experimental OEO time series

Since the early laboratory demonstrations of OCSC based on OEO was shown in 2005 [1], this kind of chaotic laser has been widely studied and become one of the most feasible hardware chaotic carrier generator. Therefore, as an experimental example, OEO systems considering nonlinear intensity modulation and nonlinear phase modulation are used to generate chaotic time series to demonstrate the effectiveness of LSTM-NN based TDS extraction method.

3.1 OEO system under intensity modulation

The experimental setup of OEO system under intensity modulation is shown in Fig. 3(a). In the system, an output light from a laser diode (LD, TLG-200) is injected into a Mach-Zehnder intensity modulator (IM, EOSPACE AX-OMSS-20-PFA-PFA-LV) which works in nonlinear region. The normalized nonlinear relationship between input and output of IM can be formulated as $co{s^2}(\cdot )$. After an optical delay fiber line (DL) and a variable optical attenuator (VOA), the output light of IM is converted into an electric signal by a photodiode (PD, PD-18G-V). Then, the output of the PD is divided into two parts, where one signal is amplified by an amplifier (AMP, GT-PA-1MHz${\sim} $1 GHz) and sent back into IM, and the other signal is received and displayed on the oscilloscope (Tektronix, DSA 72504D Digital Serial Analyzer, 25 GHz, 100GS/s).

 figure: Fig. 3.

Fig. 3. (a) Schematic diagram of the OEO system, LD, laser diode; MZM, Mach-Zehnder modulator; AMP, broadband radio frequency amplifier; PD, photodiode; DL, delay line; VOA, variable optical attenuator. (b) experimental the time series $\{{{x_i}} \}$. (c) The Normalized Root Mean Square Error (NRMSE) function versus s for experimental time series $\{{{x_i}} \}$.

Download Full Size | PDF

The loop time delay is set as $\tau = 19.2ns$ (960 points), the generated experimental waveform (50 GS/s) collected from PD is shown in Fig. 3(b). A time series $\{{{x_i}} \}$ of 5160 samples (103.2 ns) is generated for TDS extraction. Varying time delay candidate s in the range $[{0.7,\; 20.2} ]ns$ with a fixed step size $h = 0.5ns$, $E(s )$ is plotted in Fig. 3(c). It is found that when s is increased to 19.2 ns, a sharp peak will emerge from the background.

As a contrast, there is no obvious peak in ACF and DMI calculated using the same time series $\{{{x_i}} \}$, as shown in Figs. 4(a) and 4(b). When the length of time series $\{{{x_i}} \}$ increases to $2 \times {10^4}$, a peak at $s = 19.2ns$ can be clearly identified in ACF and DMI, as shown in Figs. 4(c) and 4(d). The results show that the time series $\{{{x_i}} \}$ of length $N = 5160$ is not enough for ACF and DMI, but it is enough for our method.

 figure: Fig. 4.

Fig. 4. (a) ACF and (b) DMI for $s \in [{0.7,20.2} ]ns$, and a 104.2 ns time series with 5210 data points was used. (c) ACF and (d) DMI for $s \in [{0.7,20.2} ]ns$, and a 400 ns time series with $2 \times {10^4}$ data points was used.

Download Full Size | PDF

3.2 OEO system under phase modulation

The experimental setup is shown in Fig. 5(a). An output light from a laser diode (TSL-C89-P-FA-M) is modulated by a PM (KG-PM-15-10G). After an optical delay fiber line (DL), a phase-to-intensity conversion is carried out by a fiber Bragg grating (FBG, TDCMB-C000-NCBF01). The output signal of the FBG is divided into two portions after passing through an optical coupler (OC). One signal is converted into an electric signal by PD1 after passing through a VOA, and then amplified by an AMP (Centellax, OA4MVM3). The other signal is captured and displayed by the oscilloscope after PD2.

 figure: Fig. 5.

Fig. 5. (a) The time series of the laser chaos system with optical feedback and (b) $E(s )$ function versus s for TDS extraction of system (2), and a time series with 8660 data points ($86.6ns$) was used.

Download Full Size | PDF

For time delay $\tau = 73.6ns$ (7360 points), the generated waveform (100 GS/s) is shown in Fig. 5(b). A time series of 8660 samples (86.6 ns) is used for TDS extraction. The parameters in (6)-(9) are given as $L = 200,\; L^{\prime} = 100,\delta = 10,\vartheta = 10,\; {M_1} = 90,{M_2} = 10$, and MaxEpochs is set as $\mathrm{\beth} = 20$. We vary time delay estimation $s \in [{3.6,75.6} ]ns$ to get $E(s ),$ as shown in Fig. 5(c). A sharp peak is obtained at the correct value of time delay when $s = 73.6ns$. As shown in Figs. 6(a) and (b), there is no obvious peak in ACF and DMI. But when the length of time series increases to $N = 1.2 \times {10^4}$, a successful TDS extraction via ACF and DMI is shown in Figs. 6(c) and 6(d).

 figure: Fig. 6.

Fig. 6. (a) ACF and (b) DMI for $s \in [{3.6,75.6} ]ns$, and a time series with 8660 data points (86.6 ns) was used. (c) ACF and (d) DMI for $s \in [{3.6,75.6} ]ns$, and a time series with $1.2 \times {10^4}$ data points ($120ns$) was used.

Download Full Size | PDF

3.3 OEO system under strong nonlinearity

For OEO system, the nonlinearity can be measured by Lyapunov exponents, Kaplan-Yorke dimension and the Kolmogorov-Sinai entropy. These measures increase linearly with the feedback gain $\beta $ [26]. The effectiveness of mathematical or statistical based TDS extraction methods will be reduced when increasing the feedback gain $\beta $. Therefore, the identification performance of our method under strong nonlinear effect is studied in detail by numerical approach in this section.

OEO generator can be described by the following Ikeda equation [9]:

$$x(t )+ \delta \frac{{dx(t )}}{{dt}} + \frac{1}{\theta }\smallint x(s )ds = \beta co{s^2}({{x_\tau } + \gamma } )$$
where $x(t )$ is the system state, ${x_\tau } = x({t - \tau } )$ is the time delayed state, time delay $\tau = 30ns$. $\delta = 5.3ps$ and $\theta = 5.3\mu s$ are the response time of the feedback loop, $\beta $ is the feedback gain, and $\gamma ={-} \pi /4$ is the feedback phase. In our simulation, Eq. (10) is carried out by fourth order Runge-Kutta algorithm with an integration step of $0.625ps$. We focus on the influence of strengths of nonlinearity (i.e., feedback gain $\beta $) and needed minimal data size ${M_{data}}$ on the recognition performance. As shown in Table 1, ACF loses its effectiveness when $\beta $ is larger than a threshold value of 30, and when $\beta $ increases to 15, the required data size increases to $3.7 \times {10^5}$. In contrast, the effectiveness of the LSTM-NN based method is not affected by the feedback gain $\beta $. Moreover, the required data size for our method is minimal ($32.2ns$, 5830 points, slightly beyond the actual time delay value), even if the feedback gain $\beta $ increases to 30. The reason behind this phenomenon could be explained as follows.

Tables Icon

Table 1. The minimal data size ${M_{data}}$ needed to extract the TDS for different $\beta $

The ACF is defined as [19]:

$$C(\alpha)=\frac{\langle[x(t)-\langle x(t)\rangle][x(t-\alpha)-\langle x(t)\rangle]\rangle}{[x(t)-\langle x(t)\rangle]^{2}}$$
where $\Delta $ stands for time average. ACF is used to reveal the linear statistical relationship between a variable $x(t )$ and its delayed version $x({t - \tau } )$ for TDS identification. As mentioned in [26], the linear correlation will be reduced or even eliminated when increasing the nonlinearity of the loop transformation. Due to the high nonlinearly of chaotic time series, a larger number of data are needed to extract the weak linear correlation.

The DMI is defined as [19]:

$$D(\alpha )= \sum P({x(t ),x({t - \alpha } )} )ln\frac{{P({x(t ),x({t - \alpha } )} )}}{{P({x(t )} )P({x({t - \alpha } )} )}}$$
where $P({x(t )} )$ is the probability distribution function (PDF) of $x(t )$, and $P({x(t ),x({t - \alpha } )} )$ is the joint PDF. Therefore, a large amount of data is indispensable to calculate the PDF precisely.

4. TDS extraction from simulation time series generated by a laser system

In this section, another major type of optical chaos system modeled by the Lang-Kobayashi equation is studied, including a laser system with single optical feedback, a laser chaos system with unidirectional coupled optical feedback and a chaotic laser system with a phase-modulated Sagnac loop.

4.1 Laser chaos system with single optical feedback

A laser chaos system with optical feedback which can be modeled as [31]:

$$ \left\{ \begin{array}{c} \frac{{d{\boldsymbol {E}}\left( t \right)}}{{dt}} = \frac{1}{2}\left[ {\frac{{{G_N}\left( {N\left( t \right) - {N_0}} \right)}}{{1 + \epsilon {{\boldsymbol {E}}^2}\left( t \right)}} - \frac{1}{{{\tau _p}}}} \right]{\boldsymbol {E}}\left( t \right) + k{\boldsymbol {E}}\left( {t - {\tau _0}} \right)cos\theta \left( t \right) \\ \frac{{d\emptyset \left( t \right)}}{{dt}} = \frac{\alpha }{2}\left[ {\frac{{{G_N}\left( {N\left( t \right) - {N_0}} \right)}}{{1 + \epsilon {{\boldsymbol {E}}^2}\left( t \right)}} - \frac{1}{{{\tau _p}}}} \right] - k\frac{{{\boldsymbol {E}}\left( {t - {\tau _0}} \right)}}{{{\boldsymbol {E}}\left( t \right)}}sin\theta \left( t \right) \\ \frac{{dN\left( t \right)}}{{dt}} = J - \frac{{N\left( t \right)}}{{{\tau _s}}} - \frac{{{G_N}\left( {N\left( t \right) - {N_0}} \right)}}{{1 + \epsilon {{\boldsymbol {E}}^2}\left( t \right)}}{{\boldsymbol {E}}^2}\left( t \right) \\ \emptyset \left( t \right) = \omega {\tau _0} + \emptyset \left( t \right) - \emptyset \left( {t - {\tau _0}} \right) \end{array}\right.$$

The parameters are set as ${G_N} = 8.4 \times {10^{ - 13}}{m^3}{s^{ - 1}}$, ${N_0} = 1.4 \times {10^{24}}{m^{ - 3}}$, ${\tau _p} = 1.927 \times {10^{ - 12}}s$, $k = 2.75n{s^{ - 1}}$, $\alpha = 3.0$, $J = 1.098 \times {10^{33}}{m^{ - 3}}{s^{ - 1}}$, $\epsilon = 2.5 \times {10^{ - 23}}$, and ${\tau _0} = 100.8ns$. Equation (13) is carried out by fourth order Runge-Kutta algorithm with a time step of $6.25\; ps$. A time series$\; $of 21398 samples ($133.7ns$) is used for TDS extraction. The parameters in Eqs. (6)-(9) are given as $L = 60,\; L^{\prime} = 20,\delta = 100,\vartheta = 10,\; {M_1} = 49,{M_2} = 29$, and MaxEpochs is set as $\mathrm{\beth} = 25$. Figure 7(a) shows $E(s )$ when $s \in [{40.8,{\; }101.8} ]$, and a minimum is obtained at $s = 100.8ns$. Meanwhile, there is no obvious peak in ACF and DMI, as shown in Figs. 7(b) and 7(c).

4.2 Laser chaos system with unidirectional coupled optical feedback

A laser chaos system with unidirectional coupled optical feedback can be modeled as [34]:

$$\left\{\begin{array}{c} \frac{{d{E_d}\left( t \right)}}{{dt}} = \frac{1}{2}\left[ {\frac{{\textrm{g}\left( {{N_d}\left( t \right) - {N_{0,d}}} \right)}}{{1 + \varepsilon {{\left| {{E_d}\left( t \right)} \right|}^2}}} - \frac{1}{{{\tau _{p,d}}}}} \right]{E_d}\left( t \right) + {\kappa _d}{E_d}\left( {t - {\tau _d}} \right)\cos {\mathrm{\Theta }_d}\left( t \right) \\\frac{{d{{\Phi }_d}\left( t \right)}}{{dt}} = \frac{{{\alpha _d}}}{2}\left[ {\frac{{\textrm{g}\left( {{N_d}\left( t \right) - {N_{0,d}}} \right)}}{{1 + \varepsilon {{\left| {{E_d}\left( t \right)} \right|}^2}}} - \frac{1}{{{\tau _{p,d}}}}} \right] - {\kappa _d}\frac{{{E_d}\left( {t - {\tau _d}} \right)}}{{{E_d}\left( t \right)}}\sin {\mathrm{\Theta }_d}\left( t \right) \\\frac{{d{N_d}\left( t \right)}}{{dt}} = \frac{{{I_d}}}{q} - \frac{{{N_d}\left( t \right)}}{{{\tau _{s,d}}}} - \frac{{\textrm{g}\left( {{N_d}\left( t \right) - {N_{0,d}}} \right)}}{{1 + \varepsilon {{\left| {{E_d}\left( t \right)} \right|}^2}}}{\left| {{E_d}\left( t \right)} \right|^2} \\{\mathrm{\Theta }_d}\left( t \right) = {\omega _d}{\tau _d} + {{\Phi }_d}\left( t \right) - {{\Phi }_d}\left( {t - {\tau _d}} \right) \end{array}\right.$$
$$\left\{ \begin{array}{c}\frac{{d{E_r}\left( t \right)}}{{dt}} = \frac{1}{2}\left[ {\frac{{\textrm{g}\left( {{N_r}\left( t \right) - {N_{0,r}}} \right)}}{{1 + \varepsilon {{\left| {{E_r}\left( t \right)} \right|}^2}}} - \frac{1}{{{\tau _{p,r}}}}} \right]{E_r}\left( t \right) + {\kappa _{inj}}{E_d}\left( {t - {\tau _{inj}}} \right)\cos {\mathrm{\Theta }_{inj}}\left( t \right) \\ + {\kappa _r}{E_r}\left( {t - {\tau _r}} \right)\cos {\mathrm{\Theta }_r}\left( t \right) \\\frac{{d{{\Phi }_r}\left( t \right)}}{{dt}} = \frac{{{\alpha _r}}}{2}\left[ {\frac{{\textrm{g}\left( {{N_r}\left( t \right) - {N_{0,r}}} \right)}}{{1 + \varepsilon {{\left| {{E_r}\left( t \right)} \right|}^2}}} - \frac{1}{{{\tau _{p,r}}}}} \right] - {\kappa _{inj}}\frac{{{E_d}\left( {t - {\tau _{inj}}} \right)}}{{{E_r}\left( t \right)}}\sin {\mathrm{\Theta }_{inj}}\left( t \right) \\ - {\kappa _r}\frac{{{E_r}\left( {t - {\tau _r}} \right)}}{{{E_r}\left( t \right)}}\sin {\mathrm{\Theta }_r}\left( t \right) \\ \frac{{d{N_r}\left( t \right)}}{{dt}} = \frac{{{I_r}}}{q} - \frac{{{N_r}\left( t \right)}}{{{\tau _{s,r}}}} - \frac{{\textrm{g}\left( {{N_r}\left( t \right) - {N_{0,r}}} \right)}}{{1 + \varepsilon {{\left| {{E_r}\left( t \right)} \right|}^2}}}{\left| {{E_r}\left( t \right)} \right|^2} \\{\mathrm{\Theta }_{inj}}\left( t \right) = - \left( {{\omega _d} - {\omega _r}} \right)t + {\omega _d}{\tau _{inj}} + {{\Phi }_r}\left( t \right) - {{\Phi }_d}\left( {t - {\tau _{inj}}} \right) \end{array}\right.$$

 figure: Fig. 7.

Fig. 7. (a) $E(s )$ function versus s, (b) ACF and (c) DMI for $s \in [{40.8,\; 101.8} ]ns$. A time series with 21398 data points ($133.7ns$) was used.

Download Full Size | PDF

The parameters are set as $\textrm{g} = 1.5 \times {10^4}$, $\varepsilon = 5.0 \times {10^{ - 7}}$, ${N_{0,d}} = 1.5 \times {10^8}$, ${N_{0,r}} = 1.5 \times {10^8}$, ${\tau _{p,d}} = 2.0 \times {10^{ - 12}}$, ${\tau _{p,r}} = 2.0 \times {10^{ - 12}}$, ${\tau _{s,d}} = 2.0 \times {10^{ - 9}}$, ${\tau _{s,r}} = 2.0 \times {10^{ - 9}}$, ${\tau _{inj}} = 2.0 \times {10^{ - 9}}$, ${\alpha _d} = 5.0$, ${\alpha _r} = 5.0$, ${f_d} = 193550 \times {10^9}$, ${f_r} = 193563.2 \times {10^9}$, ${\kappa _d} = 15.0 \times {10^9}$, ${\kappa _{inj}} = 20.0 \times {10^9}$, ${I_d} = 22.05 \times {10^{ - 3}}$, ${I_r} = 22.05 \times {10^{ - 3}}$, $q = 1.602 \times {10^{ - 19}}$, ${\tau _d} = 30.0 \times {10^{ - 9}}$, ${\tau _r} = 10.0 \times {10^{ - 9}}$. The system is numerically solved by the fourth order Runge-Kutta method with a time step of $10ps$. Usually, the time series of laser intensity ${I^d} = E_d^2$ and ${I^r} = E_r^2$ are available for observation and recording. A time series $\{{{I^r}} \}$ of 3102 samples (31.02 ns) is used for TDS extraction. The parameters in Eqs. (6)-(9) are given as $L = 100,\; L^{\prime} = 10,\delta = 10,\vartheta = 2,\; {M_1} = 199,{M_2} = 1$, and MaxEpochs is set as $\mathrm{\beth} = 20$. Varying $s \in [{0.5,\; 11} ]ns$, Fig. 8(a) shows the $E(s )$, and a minimum is obtained at $s = {\tau _r} = 10ns$. Meanwhile, there is no obvious peak in ACF and DMI, as shown in Figs. 8(b) and 8(c).

 figure: Fig. 8.

Fig. 8. (a) $E(s )$ function versus s for TDS extraction of system (14-15), (b) ACF and (c) DMI for $s \in [{0.5,{\; }11} ]ns$. A time series with 3102 data points ($31.02ns$) was used.

Download Full Size | PDF

4.3 Chaotic laser system with a phase-modulated Sagnac loop

A chaotic laser system with a phase-modulated Sagnac loop [35] is given as:

$$\left\{ \begin{array}{l} \frac{{d{E_1}\left( t \right)}}{{dt}} = \frac{{1 + ia}}{2}\left[ {\frac{{\textrm{g}\left[ {{N_t} - {N_0}} \right]}}{{1 + \varepsilon {{\left| {{E_1}\left( t \right)} \right|}^2}}} - \frac{1}{{{\tau _p}}}} \right]{E_1}\left( t \right)\\ + {k_1}\left[ {\frac{1}{2}{E_1}\left( {t - {\tau _1}} \right)exp\left( { - i{\omega _1}{\tau _1} + i{\phi _{PM}}} \right) + \frac{1}{2}{E_1}\left( {t - {\tau _1}} \right)exp\left( { - i{\omega _1}{\tau _1}} \right)} \right]\\ + {k_2}\left[ {\frac{1}{2}{E_2}\left( {t - {\tau _{21}}} \right)\textrm{exp}\left( { - i{\omega _1}{\tau _{21}} + i{\phi _{PM}}} \right) + \frac{1}{2}{E_2}\left( {t - {\tau _{21}}} \right)\textrm{exp}\left( { - i{\omega _1}{\tau _{21}}} \right)} \right]\\ \frac{{d{E_2}\left( t \right)}}{{dt}} = \frac{{1 + ia}}{2}\left[ {\frac{{\textrm{g}\left[ {{N_t} - {N_0}} \right]}}{{1 + \varepsilon {{\left| {{E_2}\left( t \right)} \right|}^2}}} - \frac{1}{{{\tau _p}}}} \right]{E_2}\left( t \right)\\ + {k_2}\left[ {\frac{1}{2}{E_2}\left( {t - {\tau _2}} \right)\textrm{exp}\left( { - i{\omega _2}{\tau _2} + i{\phi _{PM}}} \right) + \frac{1}{2}{E_2}\left( {t - {\tau _2}} \right)\textrm{exp}\left( { - i{\omega _2}{\tau _2}} \right)} \right]\\ + {k_1}\left[ {\frac{1}{2}{E_1}\left( {t - {\tau _{12}}} \right)\textrm{exp}\left( { - i{\omega _2}{\tau _{12}} + i{\phi _{PM}}} \right) + \frac{1}{2}{E_1}\left( {t - {\tau _{12}}} \right)\textrm{exp}\left( { - i{\omega _2}{\tau _{12}}} \right)} \right]\\ \frac{{d{N_{1,2}}\left( t \right)}}{{dt}} = \frac{{{I_{1,2}}}}{q} - \frac{{{N_{1,2}}\left( t \right)}}{{{\tau _n}}} - \frac{{\textrm{g}\left[ {{N_t} - {N_0}} \right]}}{{1 + \varepsilon {{\left| {{E_{1,2}}\left( t \right)} \right|}^2}}}{\left| {{E_{1,2}}\left( t \right)} \right|^2} \end{array} \right. $$
where ${\phi _{PM}}$ is the phase-shift induced by a PM. For TDS concealment, the PM is driven by a pseudorandom binary sequence (PRBS) with $10Gb/s$ bit rate. The other parameters used in this simulation are set as $\textrm{g} = 1.5 \times {10^4}$, $\varepsilon = 5.0 \times {10^{ - 7}}$, ${N_0} = 1.5 \times {10^8}$, $\alpha = 5.0$, ${\tau _p} = 2.0 \times {10^{ - 12}}$, ${k_1} = 15.0 \times {10^9}$, ${\tau _e} = 2.0 \times {10^{ - 9}}$, ${k_2} = 15.0 \times {10^9}$, ${\lambda _1} = 1550 \times {10^{ - 9}}$, ${\lambda _2} = 1550 \times {10^{ - 9}}$, ${I_1} = 20 \times {10^{ - 3}}$, ${I_2} = 20 \times {10^{ - 3}}$, $q = 1.602 \times {10^{ - 19}}$, ${\tau _{12}} = {\tau _{21}} = ({{\tau_1} + {\tau_2}} )/2$, ${\tau _1} = 40 \times {10^{ - 9}}$, ${\tau _2} = 30 \times {10^{ - 9}}$. The above rate equations are numerically solved by the fourth order Runge-Kutta method with a time step of $10ps$. The time series of laser intensity ${I^1} = E_1^2$ and ${I^2} = E_2^2$ are usually available for observation and recording. A time series $\{{{I^1}} \}$ of 14062 samples (140.62 ns) is used for TDS extraction. The parameters in (6)-(9) are given as $L = 100,\; L^{\prime} = 10,\delta = 50,\vartheta = 2,\; {M_1} = 199,{M_2} = 1$, and MaxEpochs is set as $\mathrm{\beth} = 10$. Figure 9(a) shows $E(s )$ when $s \in [{2,\; 42} ]ns$, and a minimum is obtained at $s = 40ns$. Meanwhile, there is no obvious peak in ACF and DMI, as shown in Figs. 9(b) and 9(c).

5. Analysis on noise added to OEO time series

In this section, we focus on the influence of noise level on the recognition performance. In practical optical communication systems, optical and electrical noise should be considered simultaneously. For simplification, here we use additive white Gaussian noise to evaluate the robustness of our method.

 figure: Fig. 9.

Fig. 9. (a) $E(s )$ function versus s for TDS extraction of system (16), (b) ACF and (c) DMI for $s \in [{2,{\; }42} ]ns$. A time series with 14062 data points ($31.02ns$) was used.

Download Full Size | PDF

Simulation is carried out for the OEO generator (10). For different signal-to-noise ratio (SNR), the comparisons are shown in Table 2 when $\beta = 10$ (a relatively large value considering the limited implementation capability of actual devices). It shows that all these methods have a strong ability of handling noise. And the required data size for LSTM-NN is still minimal, and it marginally increases as the noise level increases.

Tables Icon

Table 2. The minimal data size ${M_{data}}$ needed to extract the TDS for different SNR, $\beta = 10$.

Generally, the effectiveness of both CNN and LSTM-NN based TDS extraction methods is built on the fact that there is an approximately fixed nonlinear functional relationship between the chaotic signal and its delayed version when the delay is equal to the system loop delay. In our previous work, CNN is used to identify such functional feature indirectly through the images which are generated by reconstructing the chaotic time series into a two-dimensional phase space. Therefore, prior knowledge about the chaotic system model is indispensable in the training process to generate the labelled image set. However, in our current research, LSTM-NN is used to directly model the delay function candidates ${\varphi _s}$ from the output time series, and TDS identification is established by evaluate the prediction ability of ${\varphi _s}$. As a result, model information of the chaotic system is no longer needed and blind identification can be established.

6. Conclusion

We conclude that a data driven LSTM-NN model can identify TDS of time-delayed optical chaotic systems. Compared with existing methods, including ACF, DMI, and CNN based method, LSTM-NN based method shows its advantages, which lies in the following aspects: 1) No prior knowledge about chaotic system structure or mathematical model is required; 2) Minimum data size requirements; 3) High tolerance to system nonlinearity and additive noise. In the case of chaotic OEO, the proposed method is successfully demonstrated under strong nonlinearity. The proposed strategy also has the potential to be used for detecting the vulnerability of other security enhanced OCSC systems as discussed in section 4, and this will be the focus of our research in the next stage. From the perspective of safety assessment, we suggest that the ability of resisting deep-learning-based analysis should be considered in designing physical layer secure communication schemes. Subsequently, this type of analysis method will help to promote the security level and further maturity of optical chaos secure communications.

Funding

Open Research Project of The Hubei Key Laboratory of Intelligent Geo-Information Processing (KLIGIP-2019B11); National Key Research and Development Program of China (2018YFB1801304).

Disclosures

The authors declare no conflicts of interest.

References

1. A. Argyris, D. Syvridis, L. Larger, V. Annovazzi-Lodi, P. Colet, I. Fischer, J. Garcia-Ojalvo, C. R. Mirasso, L. Pesquera, and K. A. Shore, “Chaos-based communications at high bit rates using commercial fibre-optic links,” Nature 438(7066), 343–346 (2005). [CrossRef]  

2. R. Lavrov, M. Jacquot, and L. Larger, “Nonlocal nonlinear electro-optic phase dynamics demonstrating 10gb/s chaos communications,” IEEE J. Quantum Electron. 46(10), 1430–1435 (2010). [CrossRef]  

3. J. Ke, L. Yi, Z. Yang, Y. Yang, Q. Zhuge, Y. Chen, and W. Hu, “32 Gbp/s chaotic optical communications by deep-learning-based chaos synchronization,” Opt. Lett. 44(23), 5776–5779 (2019). [CrossRef]  

4. A. K. Zhao, N. Jiang, S. Q. Liu, Y. Q. Zhang, and K. Qiu, “Generation of synchronized wideband complex signals and its application in secure optical communication,” Opt. Express 28(16), 23363–23373 (2020). [CrossRef]  

5. L. S. Wang, X. X. Mao, A. B. Wang, Y. C. Wang, Z. S. Gao, S. S. Li, and L. S. Yan, “Scheme of coherent optical chaos communication,” Opt. Lett. 45(17), 4762–4765 (2020). [CrossRef]  

6. R. H. Zhang, P. Zhou, Y. G. Yang, Q. Fang, P. H. Mu, and N. Q. Li, “Enhancing time-delay suppression in a semiconductor laser with chaotic optical injection via parameter mismatch,” Opt. Express 28(5), 7197–7206 (2020). [CrossRef]  

7. M. Li, X. Y. Zhang, Y. F. Hong, Y. Z. Zhang, Y. C. Shi, and X. F. Chen, “Confidentiality-enhanced chaotic optical communication system with variable RF amplifier gain,” Opt. Express 27(18), 25953–25963 (2019). [CrossRef]  

8. W. Y. Yang, G. Q. Xia, Z. F. Jiang, E. Jayaprasath, C. X. Hu, and Z. M. Wu, “Numerical investigations on multi-channel wideband chaotic signal generation by a multi-transverse mode vertical-cavity surface-emitting laser subject to chaotic optical injection,” Appl. Opt. 58(30), 8160–8166 (2019). [CrossRef]  

9. Y. C. Kouomou, P. Colet, L. Larger, and N. Gastaud, “Chaotic Breathers in Delayed Electro-Optical Systems,” Phys. Rev. Lett. 95(20), 203903 (2005). [CrossRef]  

10. C. Wang, Y. F. Ji, H. X. Wang, and L. Bai, “Security-Enhanced Electro-Optic Feedback Phase Chaotic System Based on Nonlinear Coupling of Two Delayed Interfering Branches,” IEEE Photonics J. 10(4), 1–15 (2018). [CrossRef]  

11. N. Q. Li, W. Pan, S. Y. Xiang, B. Luo, L. S. Yan, and X. H. Zou, “Hybrid chaos-based communication system consisting of three chaotic semiconductor ring lasers,” Appl. Opt. 52(7), 1523–1530 (2013). [CrossRef]  

12. S. Y. Xiang, A. J. Wen, W. Pan, L. Lin, H. X. Zhang, H. Zhang, X. X. Guo, and J. F. Li, “Suppression of Chaos Time Delay Signature in a Ring Network Consisting of Three Semiconductor Lasers Coupled With Heterogeneous Delays,” J. Lightwave Technol. 34(18), 4221–4227 (2016). [CrossRef]  

13. Y. N. Han, S. Y. Xiang, Y. Wang, Y. T. Ma, B. Wang, A. J. Wen, and Y. Hao, “Generation of multi-channel chaotic signals with time delay signature concealment and ultrafast photonic decision making based on a globally-coupled semiconductor laser network,” Photonics Res. 8(11), 1792–1799 (2020). [CrossRef]  

14. J. C. Bai, H. X. Wang, and Y. F. Ji, “Time-delay signature concealing electro-optic chaotic system with multiply feedback nonlinear loops,” Opt. Express 29(2), 706–718 (2021). [CrossRef]  

15. J. J. Chen, C. X. Li, X. L. Mu, L. F. Li, and Y. N. Duan, “Time-delay signature suppression of polarization-resolved wideband chaos in VCSELs with dual-path chaotic optical injections,” Appl. Opt. 59(24), 7217–7224 (2020). [CrossRef]  

16. L. Zunino, M. C. Soriano, I. Fischer, O. A. Rosso, and C. R. Mirasso, “Permutation-information-theory approach to unveil delay dynamics from time-series analysis,” Phys. Rev. E 82(4), 046212 (2010). [CrossRef]  

17. M. C. Soriano, L. Zunino, O. A. Rosso, I. Fischer, and C. R. Mirasso, “Time scales of a chaotic semiconductor laser with optical feedback under the lens of a permutation information analysis,” IEEE J. Quantum Electron. 47(2), 252–261 (2011). [CrossRef]  

18. Y. C. Wang, J. S. Liang, A. B. Wang, and M. J. Zhang, “Time-delay extraction in chaotic laser diode using RF spectrum analyser,” Electron. Lett. 46(24), 1621–1623 (2010). [CrossRef]  

19. Y. Wu, Y. C. Wang, P. Li, A. B. Wang, and M. J. Zhang, “Can fixed time delay signature be concealed in chaotic semiconductor laser with optical feedback?” IEEE J. Quantum Electron. 48(11), 1371–1379 (2012). [CrossRef]  

20. A. C. Fowler and G. Kember, “Delay recognition in chaotic time series,” Phys. Lett. A 175(6), 402–408 (1993). [CrossRef]  

21. H. Voss and J. Kurths, “Reconstruction of non-linear time delay models from data by the use of optimal transformation,” Phys. Lett. A 234(5), 336–344 (1997). [CrossRef]  

22. D. Rontani, A. Locquet, M. Sciamanna, D. S. Citrin, and S. Ortin, “Time-Delay Identification in a Chaotic Semiconductor Laser With Optical Feedback: A Dynamical Point of View,” IEEE J. Quantum Electron. 45(7), 879–1891 (2009). [CrossRef]  

23. R. Hegger, M. J. Bünner, and H. Kantz, “Identifying and Modeling Delay Feedback Systems,” Phys. Rev. Lett. 81(3), 558–561 (1998). [CrossRef]  

24. V. S. Udaltsov, L. Larger, and J. P. Goedgebuer, “Time delay identification in chaotic cryptosystems ruled by delay-differential equations,” J. Opt. Technol. 72(5), 373–377 (2005). [CrossRef]  

25. A. Elsonbaty, S. F. Hegazy, and S. S. A. Obayya, Senior, “Simultaneous Suppression of Time-Delay Signature in Intensity and Phase of Dual-Channel Chaos Communication,” IEEE J. Quantum Electron. 51(9), 1–9 (2015). [CrossRef]  

26. Silvia Ortín and L. Pesquera, “Unmasking Optical Chaotic Cryptosystems Based on Delayed Optoelectronic Feedback. Chaos Synchronization and Cryptography for Secure Communications: Applications for Encryption, edited by Santo Banerjee,” IGI Global, 386–414 (2011).

27. J. S. Zhang and X. C. Xiao, “Predicting chaotic time series using recurrent neural network,” Chin. Phys. Lett. 17(2), 88–90 (2000). [CrossRef]  

28. S. Muzaffar and A. Afshari, “Short-Term Load Forecasts Using LSTM Networks,” Energy Procedia 158, 2922–2927 (2019). [CrossRef]  

29. Y. Duan, L. V. Yisheng, and F. Wang, “Travel time prediction with LSTM neural network,” IEEE 19th International Conference on Intelligent Transportation Systems (ITSC), 1053–1058 (2016).

30. K. Yeo, “Model-free prediction of noisy chaotic time series by deep learning,” CoRR 1710.01693 (2017).

31. Y. Chen, R. Xin, M. Cheng, X. Gao, S. Li, W. Shao, L. Deng, M. Zhang, S. Fu, and D. Liu, “Unveil the time delay signature of optical chaos systems with a convolutional neural network,” Opt. Express 28(10), 15221–15231 (2020). [CrossRef]  

32. C. H. Yang and H. Y. Shen, “Analysis and Prediction of Chaotic Time Series Based on Deep Learning Neural Networks,” 2020 International Conference on System Science and Engineering (ICSSE), 2020, pp. 1–9.

33. M. Sangiorgio and F. Dercole, ““Robustness of LSTM neural networks for multi-step forecasting of chaotic time series,” Chaos, Solitons Fractals 139, 110045 (2020). [CrossRef]  

34. A. Uchida, “Optical communication with chaotic lasers: applications of nonlinear dynamics and synchronization,” (John Wiley & Sons, 2012).

35. Y. T. Ma, S. Y. Xiang, X. X. Guo, Z. W. Song, A. J. Wen, and Y. Hao, “Time-delay signature concealment of chaos and ultrafast decision making in mutually coupled semiconductor lasers with a phase-modulated Sagnac loop,” Opt. Express 28(2), 1665–1678 (2020). [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 (9)

Fig. 1.
Fig. 1. The principle of the proposed strategy.
Fig. 2.
Fig. 2. Data driven self-modelling TDS extraction via NN
Fig. 3.
Fig. 3. (a) Schematic diagram of the OEO system, LD, laser diode; MZM, Mach-Zehnder modulator; AMP, broadband radio frequency amplifier; PD, photodiode; DL, delay line; VOA, variable optical attenuator. (b) experimental the time series $\{{{x_i}} \}$. (c) The Normalized Root Mean Square Error (NRMSE) function versus s for experimental time series $\{{{x_i}} \}$.
Fig. 4.
Fig. 4. (a) ACF and (b) DMI for $s \in [{0.7,20.2} ]ns$, and a 104.2 ns time series with 5210 data points was used. (c) ACF and (d) DMI for $s \in [{0.7,20.2} ]ns$, and a 400 ns time series with $2 \times {10^4}$ data points was used.
Fig. 5.
Fig. 5. (a) The time series of the laser chaos system with optical feedback and (b) $E(s )$ function versus s for TDS extraction of system (2), and a time series with 8660 data points ($86.6ns$) was used.
Fig. 6.
Fig. 6. (a) ACF and (b) DMI for $s \in [{3.6,75.6} ]ns$, and a time series with 8660 data points (86.6 ns) was used. (c) ACF and (d) DMI for $s \in [{3.6,75.6} ]ns$, and a time series with $1.2 \times {10^4}$ data points ($120ns$) was used.
Fig. 7.
Fig. 7. (a) $E(s )$ function versus s, (b) ACF and (c) DMI for $s \in [{40.8,\; 101.8} ]ns$. A time series with 21398 data points ($133.7ns$) was used.
Fig. 8.
Fig. 8. (a) $E(s )$ function versus s for TDS extraction of system (14-15), (b) ACF and (c) DMI for $s \in [{0.5,{\; }11} ]ns$. A time series with 3102 data points ($31.02ns$) was used.
Fig. 9.
Fig. 9. (a) $E(s )$ function versus s for TDS extraction of system (16), (b) ACF and (c) DMI for $s \in [{2,{\; }42} ]ns$. A time series with 14062 data points ($31.02ns$) was used.

Tables (2)

Tables Icon

Table 1. The minimal data size M d a t a needed to extract the TDS for different β

Tables Icon

Table 2. The minimal data size M d a t a needed to extract the TDS for different SNR, β = 10 .

Equations (16)

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

x ˙ ( t ) = f ( x ( t ) , x ( t τ ) ) ,
φ : x ( t ) x ( t + τ ) ,
φ s : { x i } n m { x i } n + s m + s .
E ( s ) = i = p q ( x i + s j y i ) 2 i = p q ( x i + s j x ¯ ) 2 for s = s j
x ( t + Δ t ) x ( t ) = Δ t f ( x ( t ) , x ( t τ ) ) ,
X T r a i n l = { x i } 1 + l δ L + l δ ,
Y T r a i n l = { x i } 1 + s + l δ L + s + l δ { x i } s + l δ L + s + l δ 1 ,
X T e s t l = { x i } L + M 1 δ + 1 + l ϑ L + L + M 1 δ + 1 + l ϑ ,
Y T e s t l = { x i } L + M 1 δ + 1 + l ϑ + s L + L + M 1 δ + 1 + l ϑ + s ,
x ( t ) + δ d x ( t ) d t + 1 θ x ( s ) d s = β c o s 2 ( x τ + γ )
C ( α ) = [ x ( t ) x ( t ) ] [ x ( t α ) x ( t ) ] [ x ( t ) x ( t ) ] 2
D ( α ) = P ( x ( t ) , x ( t α ) ) l n P ( x ( t ) , x ( t α ) ) P ( x ( t ) ) P ( x ( t α ) )
{ d E ( t ) d t = 1 2 [ G N ( N ( t ) N 0 ) 1 + ϵ E 2 ( t ) 1 τ p ] E ( t ) + k E ( t τ 0 ) c o s θ ( t ) d ( t ) d t = α 2 [ G N ( N ( t ) N 0 ) 1 + ϵ E 2 ( t ) 1 τ p ] k E ( t τ 0 ) E ( t ) s i n θ ( t ) d N ( t ) d t = J N ( t ) τ s G N ( N ( t ) N 0 ) 1 + ϵ E 2 ( t ) E 2 ( t ) ( t ) = ω τ 0 + ( t ) ( t τ 0 )
{ d E d ( t ) d t = 1 2 [ g ( N d ( t ) N 0 , d ) 1 + ε | E d ( t ) | 2 1 τ p , d ] E d ( t ) + κ d E d ( t τ d ) cos Θ d ( t ) d Φ d ( t ) d t = α d 2 [ g ( N d ( t ) N 0 , d ) 1 + ε | E d ( t ) | 2 1 τ p , d ] κ d E d ( t τ d ) E d ( t ) sin Θ d ( t ) d N d ( t ) d t = I d q N d ( t ) τ s , d g ( N d ( t ) N 0 , d ) 1 + ε | E d ( t ) | 2 | E d ( t ) | 2 Θ d ( t ) = ω d τ d + Φ d ( t ) Φ d ( t τ d )
{ d E r ( t ) d t = 1 2 [ g ( N r ( t ) N 0 , r ) 1 + ε | E r ( t ) | 2 1 τ p , r ] E r ( t ) + κ i n j E d ( t τ i n j ) cos Θ i n j ( t ) + κ r E r ( t τ r ) cos Θ r ( t ) d Φ r ( t ) d t = α r 2 [ g ( N r ( t ) N 0 , r ) 1 + ε | E r ( t ) | 2 1 τ p , r ] κ i n j E d ( t τ i n j ) E r ( t ) sin Θ i n j ( t ) κ r E r ( t τ r ) E r ( t ) sin Θ r ( t ) d N r ( t ) d t = I r q N r ( t ) τ s , r g ( N r ( t ) N 0 , r ) 1 + ε | E r ( t ) | 2 | E r ( t ) | 2 Θ i n j ( t ) = ( ω d ω r ) t + ω d τ i n j + Φ r ( t ) Φ d ( t τ i n j )
{ d E 1 ( t ) d t = 1 + i a 2 [ g [ N t N 0 ] 1 + ε | E 1 ( t ) | 2 1 τ p ] E 1 ( t ) + k 1 [ 1 2 E 1 ( t τ 1 ) e x p ( i ω 1 τ 1 + i ϕ P M ) + 1 2 E 1 ( t τ 1 ) e x p ( i ω 1 τ 1 ) ] + k 2 [ 1 2 E 2 ( t τ 21 ) exp ( i ω 1 τ 21 + i ϕ P M ) + 1 2 E 2 ( t τ 21 ) exp ( i ω 1 τ 21 ) ] d E 2 ( t ) d t = 1 + i a 2 [ g [ N t N 0 ] 1 + ε | E 2 ( t ) | 2 1 τ p ] E 2 ( t ) + k 2 [ 1 2 E 2 ( t τ 2 ) exp ( i ω 2 τ 2 + i ϕ P M ) + 1 2 E 2 ( t τ 2 ) exp ( i ω 2 τ 2 ) ] + k 1 [ 1 2 E 1 ( t τ 12 ) exp ( i ω 2 τ 12 + i ϕ P M ) + 1 2 E 1 ( t τ 12 ) exp ( i ω 2 τ 12 ) ] d N 1 , 2 ( t ) d t = I 1 , 2 q N 1 , 2 ( t ) τ n g [ N t N 0 ] 1 + ε | E 1 , 2 ( t ) | 2 | E 1 , 2 ( t ) | 2
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.