Impact of the Chaotic Synchronization’s Stability on the Performance of QCPSK Communication System

The current work presents a study of the implementation of a quadrature chaos phase-shift keying communication system (QCPSK) based on the employment of different chaos oscillators. The research takes two directions, with one being the study of the chaos synchronization’s noise immunity for several chaos oscillators that are the potential core blocks of the QCPSK system. The correlation coefficient over time is used to estimate the synchronization noise immunity. The second direction is the estimation of the QCPSK system’s baseband model performance in the AWGN propagation channel using the bit error ratio (BER) as the estimation method for several chaos oscillators employed as the core of the QCPSK system’s model.


Introduction
The reason for the increasing interest in the employment of chaos phenomenon as a building block in the design of communication systems is supported by the characteristics of the chaos-based communication systems: non-periodic, noise-like signal waveform, wide frequency bandwidth, and sensitive dependence on initial conditions. These properties allow the additional physical layer of security for the information transfer. In the ongoing digital transformation of the world in such directions as the internet of things (IoT) and wireless sensor networks (WSNs), there is a demand for highly efficient, wideband, secure communication, that could be ensured by the use of chaos-based data transmission systems [1,2]. Having grasped the benefits of the employment of chaos in data communication, the various chaos-based modulation schemes were introduced and studied. These modulation schemes are generally divided into two broad groups: noncoherent modulation schemes and coherent modulation schemes.
In the case of coherent modulation schemes, the in-phase synchronization between the receiver and the transmitter must be sustained. The modulation schemes that refer to this group are chaos-shift keying (CSK) [3,4], quadrature chaos shift-keying (QCSK) [5], chaos-based code division multiple access (CDMA) [6]. In the CSK system, two or more chaos generators with different parameters are located in the transmitter. The information is transmitted utilizing switching between signals from the corresponding chaos generators. The QCSK system adds the employment of the orthogonal signals, increasing the rate of band usage.
In the case of noncoherent modulation schemes, in-phase synchronization is not required. The modulation schemes of this group are chaos on-off keying (COOK) [7], differential chaos-shift keying (DCSK) [8], frequency-modulated differential chaos-shift keying (FM-DCSK) [9]. In COOK, the chaos generator is turned on when the "1" bit is transmitted and off when the "0" bit is transmitted. The receiver then estimates the energy of the received signal. For the received "1" bit, the energy of the signal is notably higher.

Communication System Based on Quadrature Chaos Phase-Shift Keying
The current section presents the model of a communication system first proposed with the description of the applied elements and their functions in [16]. The block diagram for the communication system under investigation is shown in Figure 1.
Electronics 2021, 10, x FOR PEER REVIEW 2 of 15 the received signal. For the received "1" bit, the energy of the signal is notably higher. The DCSK modulation scheme uses two consecutive parts of the chaotic signal in one symbol transmission, the first part is a reference signal, and the second part bears the information. The FM-DCSK utilizes frequency modulation performed by chaotic signals. Complementing many possible chaos modulations, several approaches in implementing chaos signal generators exist. The first type of chaotic generators are known as chaotic maps, also referred to as chaotic functions [10,11]. Another approach for generating chaotic signals is to employ chaotic systems defined with a set of differential equations. Chaotic systems of such type are implemented using either analog circuit arithmetic [12] or analog circuit state variables [13][14][15].
In [16], we proposed a model of quadrature chaos phase-shift keying-based communication system (QCPSK) employing chaos oscillators (defined by state variable differential equations). The current paper further investigates the proposed system using some approaches for the study developed in Reference [17].
The paper is organized in the following way: the model of the communication system based on chaos phase-shift keying and the quadrature transceiver is presented in Section 2. Section 3 describes the chaos oscillators employed in the research of the presented communication system. In Section 4, the first step of the research is to study the synchronization properties of the proposed chaos oscillators using simulation, hardware experiments, and analysis of the results acquired in this step. Section 5 describes the second step of the research, which is focused on investigating the performance of the proposed communication system with the employment of the chaos oscillators studied in the first step. The performance is estimated using the baseband model of the communication system in the AWGN propagation channel. Section 5 also includes the analysis of the results acquired in the second step of the research. The final section presents conclusions.  [16]. X and Y in Figure 1 represent two systems exhibiting chaotic behavior. These systems are described by the state variables, which are 1 X , 2 X , 3 X for the X system, Y for the Y system. The state variables of chaotic systems are also referred to as chaotic signals. The "master" device utilizes 1 X and 2 X state variables for data transmission. The state variable 1 X is used to encode the binary message ( ) m t , ensured by a corresponding data-driven switch. When the "1" bit is transmitted, the output of the switch is equal to 1 X , and when the "0" bit is transmitted, the output of the switch is equal to 1 X − . In the "master" device synchronization and the encoded message signals are then fed into the quadrature modulator. The 2 X , which is used as a synchronization X and Y in Figure 1 represent two systems exhibiting chaotic behavior. These systems are described by the state variables, which are X 1 , X 2 , X 3 for the X system, and Y 1 , Y 2 , Y 3 for the Y system. The state variables of chaotic systems are also referred to as chaotic signals. The "master" device utilizes X 1 and X 2 state variables for data transmission. The state variable X 1 is used to encode the binary message m(t), ensured by a corresponding datadriven switch. When the "1" bit is transmitted, the output of the switch is equal to X 1 , and when the "0" bit is transmitted, the output of the switch is equal to −X 1 . In the "master" device synchronization and the encoded message signals are then fed into the quadrature modulator. The X 2 , which is used as a synchronization signal between chaotic systems X and Y, is the in-phase component (commonly referred to as "I") of the baseband signal. It is up-converted to the passband by multiplying it with a cosine waveform having a frequency w 0 , which is regarded as the carrier frequency. In turn, the message signal is the quadrature component (commonly referred to as "Q") of the baseband signal and it is up-converted to the passband by multiplying it with a sine waveform having the same carrier frequency. The master device transmits the modulated harmonic signal, which is the sum of up-converted I and Q components. The environment between the transmitter and the receiver is modeled as a direct additive white Gaussian noise (AWGN) channel [16].
In the "slave" device, the quadrature demodulator and the low-pass filters (LFP) are used to down-convert the received passband signal from the carrier frequency and separate synchronization signal carried by the "I" component, and encoded message signal carried by the "Q" component of the baseband signal. The synchronization signal X 2 substitutes Y 2 , thereby causing Y 1 repeat behavior of X 1 and Y 3 , repeat behavior of X 3 . This type of chaotic synchronization is called "substitution-based chaos synchronization". Studies [10,[17][18][19] demonstrated that replacing a single state variable of a chaotic system Y with a similar state variable from a similar chaotic system X results in Y system, following the behavior of X. In this configuration, the X system acts as "master" and Y system acts as "slave". Moreover, Section 4 introduces differential equitation systems of investigated chaotic oscillators showing how state variables impact each other. Differential equitation systems of "master" and "slave" systems are absolutely the same, except substitution.
Y To decode the received message signal, a time-windowed correlation coefficient β is calculated for a pair of Y 1 and the received encoded message signal (which is either X 1 or −X 1 ). The correlation coefficient β is greater than 0 when X 1 is received, and less than 0 when −X 1 is received. At perfect conditions β = 1 for X 1 , and β = −1 for −X 1 , β = 1 means Y 1 and X 1 are precisely identical, while β = −1 shows that Y 1 and X 1 are precisely opposite because of −X 1 . β = 0 indicates the absence of correlation between the two signals. The final block interprets the results of the correlation coefficient as binary "1" if β > 0 and "0" if β < 0, thereby acquiring the decoded message signal m(t) [16].

Chaos Oscillators
As shown in the previous section, the most important components of the proposed communication system defining its performance are chaotic oscillators. The chaos oscillators selected for this research were: the "Simple Two-Transistor Single-Supply Resistor-Capacitor Chaotic Oscillator" [15], "Simple RC chaotic oscillator" [14] and "A simple chaotic oscillator for educational purposes" [13] (also referred to as "Vilnius chaotic oscillator" [17]). For the sake of simplicity, the following denotation is used for the mentioned chaos oscillators throughout this paper: "RC1 chaos oscillator", "RC2 chaos oscillator" and "Vilnius chaos oscillator." This section further describes each of the mentioned chaos oscillators starting with the RC1. The circuit diagram of every mentioned chaos oscillator is shown in Figure 2. From Figure 2a, the mathematical model of the RC1 oscillator is represented with a set of differential equations, Equation (1) [15,16], where the state variables are the voltages across the capacitors 1 C -4 C .
( ) From Figure 2a, the mathematical model of the RC1 oscillator is represented with a set of differential equations, Equation (1) [15,16], where the state variables are the voltages across the capacitors C 1 -C 4 .
where R A = R 1 = R 7 = 5.6 kΩ; R B = R 2 = R 3 = 10 kΩ; R 4 = 15 kΩ; R 5 = 33 kΩ; R 6 = 47 kΩ; C = C 1 = C 2 = C 3 = 1 nF; C 4 = 330 pF; v C1 -voltage across the capacitor C 1 ; v C2 -voltage across the capacitor C 2 ; i C1 -current through the capacitor C 1 ; i C4 -current through the capacitor C 4 ; v BE1 -base-emitter voltage of the 2N3904 transistor Q 1 (is equal to the voltage across the capacitor C 3 ); v BE2 -base-emitter voltage of the 2N3904 transistor Q 2 (is equal to the voltage across the capacitor C 4 ); V 1 -5.8 V power supply voltage [15,16]. The circuit diagram of the RC2 chaos oscillator, shown in Figure 2b, could be described by the differential equations presented in Equation (2) [14]. From Equations (2)-(4), it is seen that the state variables are the voltages across the capacitors C 1 , C 2 , C 3 , as well as the output voltages of the operational amplifiers OA 1 and OA 2 .
The coefficients of the differential equations for this chaos oscillator are defined from the Figure 2b circuit parameters, summarized in [14]: For the circuit in Figure 2b, two additional equations describing the outputs of the operational amplifiers OA 1 and OA 2 are provided as follows: In [14].
The circuit diagram of the Vilnius chaos oscillator is shown in Figure 2c with the mathematical model described by a system of differential equations in Equation (5) [13,17], defining the state variables for the Vilnius chaos oscillator as voltages across the capacitors C 1 and C 2 , and a current through the inductor L 1 . where C 1 = 1 nF; C 2 = 150 pF; L 1 = 1 mH; R = R 1 = 1 kΩ; R 2 = R 3 = 10 kΩ; R 4 = 20 kΩ; R 5 = 2.9 Ω; v C -voltage across the capacitor C 1 ; v C 2 -voltage across the capacitor C 2 ; i L -current through the inductor L 1 ; i D -current through the 1N4148 diode D 1 ; i 0 -current through the resistor R 4 ; k-is a gain of the TL082 operational amplifier OA [17]. In Figure 2c,

Study of the Properties of Chaos Oscillators
The current section describes the first step in the investigation of the proposed communication system's performance with the employment of selected chaos oscillators. The first step is the study of noise immunity of chaos oscillator's synchronization. The importance of this step is approved by the fact that the model of the communication system in Figure 1 is based on the synchronization of two chaotic systems, and the noise immunity of the chaotic synchronization will directly affect the performance of the entire system. This section considers the model of the experiment, complemented with simulation and hardware implementations. The block diagram of the experimental study is depicted in Figure 3.
-current through the resistor 4 R ; k -is a gain of the TL082 operational amplifier OA [17]. In Figure 2c,

Study of the Properties of Chaos Oscillators
The current section describes the first step in the investigation of the proposed communication system's performance with the employment of selected chaos oscillators. The first step is the study of noise immunity of chaos oscillator's synchronization. The importance of this step is approved by the fact that the model of the communication system in Figure 1 is based on the synchronization of two chaotic systems, and the noise immunity of the chaotic synchronization will directly affect the performance of the entire system. This section considers the model of the experiment, complemented with simulation and hardware implementations. The block diagram of the experimental study is depicted in Figure 3.  The experimental model in Figure 3 is resembling the model of the communication system from Figure 1, excluding the message encoder, quadrature modulator-demodulator, and the decision-making device. The model consists of two chaotic systems X (master) and Y (slave). State variables X 2 and Y 2 are used for chaos synchronization, much like in Figure 1. For the study of synchronization's noise immunity, the environment between the two oscillators is an additive white Gaussian noise channel with a limited frequency bandwidth. The use of band-limited noise signals reflects the influence of low-pass filters (LPF) in Figure 1. This means that the noise signal with the same frequency band as the received chaotic signals will affect the Y system. The frequency bandwidth of the noise signals n 1 and n 2 is equal to the frequency band of the corresponding signals X 1 and X 2 . The noise signal n 2 is added to the synchronization signal X 2 and results in a signal X 2 + n 2 . This signal is applied to the Y system for the Y 2 substitution. The other signal that is being transmitted to the Y system through the noisy environment is X 1 that results in the X 1 + n 1 signal. The simultaneous presence of the noise in the channel for both transmitted signals is assumed as the worst-case scenario. Synchronization of the two systems is estimated using the 100 µs window correlation coefficient β for the signals X 1 + n 1 and Y 1 . The correlation coefficient's values are in the range from 0 to 1, as only the direct X 1 signal is taken to calculate correlation with Y 1 . The higher value β indicates better synchronization. The window of the correlation coefficient has a fixed time frame like in  Figure 1. The equation used for calculating the correlation coefficient between the given signals is as follows: where β-correlation coefficient; m, k-array indexes. The n 1 , n 2 noise signal levels are measured in SNR (signal to noise ratio) from 20 dB to 0 dB with a step of 5 dB. The experiment is applied to all of the three previously shown chaos oscillators.

Simulation-Based Study
This section presents the simulation setup for the experiment shown in Figure 3. The simulation is done using LTspice software for circuit simulation and MATLAB for the required noise signals generation and calculation of the correlation coefficient. The simulation setup for this experiment is shown in Figure 4 using the Vilnius chaos oscillator as an example [17].  Elements in Figure 4 have the following values: OA are TL082 operational amplifiers. All the operational amplifiers in Figure 4 are supplied with ±3 V. In Figure 4, voltage signal applied to the capacitor 2S C , like the signal 2 2 X n + in Figure 3. The signal 1 1 X n + is obtained summing up the noise signal 1 n with the capacitor's 1M C voltage signal [17]. Equation (7) shows the state variable differential equation of the "slave" circuit with . Equation (7)  follows C v 1M -the "slave" oscillator follows the "master" oscillator. Synchronization for chaos oscillators performs exactly as previously described for X and Y chaotic systems in Figure 1. In Equation (7), the influence of noise 2 V on the "slave" chaos oscillator during synchronization can be seen-increasing noise level will deviate the "slave" circuit from following the "master." Elements in Figure 4 have the following values: R 6 = R 8 = R 9 = 1 kΩ; R 7 = 2 kΩ; R 10 = R 11 = 200 Ω; OA 1 -OA 4 are TL082 operational amplifiers. All the operational amplifiers in Figure 4 are supplied with ±3 V. In Figure 4, X 1 is the capacitor's C 1M voltage, X 2 is the capacitor's C 2M voltage, Y 1 is the capacitor's C 1S voltage, Y 2 is the capacitor's C 2S voltage. A signal V 2 is the noise n 2 voltage that is added to the capacitor's C 2M voltage signal applied to the capacitor C 2S , like the signal X 2 + n 2 in Figure 3. The signal X 1 + n 1 is obtained summing up the noise signal n 1 with the capacitor's C 1M voltage signal [17]. Equation (7) shows the state variable differential equation of the "slave" circuit with C 2S voltage v C 2S substituted by C 1M voltage with noise v C 2M + V 2 . Equation (7) demonstrates substitution-based chaos synchronization in action: v C 2S state variable of the "slave" chaos oscillator is replaced with v C 2M ; assuming V 2 = 0 V. As v L 1S , v C 1S , v C 2S are related, v L 1S follows the behavior of v L 1M , and v C 1S follows v C 1M -the "slave" oscillator follows the "master" oscillator. Synchronization for chaos oscillators performs exactly as previously described for X and Y chaotic systems in Figure 1. In Equation (7), the influence of noise Electronics 2021, 10, 640 7 of 14 V 2 on the "slave" chaos oscillator during synchronization can be seen-increasing noise level will deviate the "slave" circuit from following the "master." In the case of RC1 chaos oscillator, X 1 is the voltage across the C 1M capacitor and the Y 1 is the voltage across C 1S (C 1 in Figure 2a), X 2 is the voltage of the C 2M capacitor, and the Y 2 is the voltage of the C 2S (C 2 in Figure 2a). In the case of the RC2 chaos oscillator, X 1 is the voltage signal at the output of the operational amplifier OA 2M and Y 1 is the voltage at the output of the OA 2S amplifier (OA 2 in Figure 2b). X 2 is the voltage of the C 2M and the Y 2 is the voltage of the C 2S (C 2 in Figure 2b).

Hardware-Based Study
This section presents the hardware implementation (see Figure 5) of the experimental model shown in Figure 3. For the hardware experiment, the chaos oscillators are implemented in specially-designed PCBs. A pair of chaos oscillators is powered by two separate independent power supplies. The chaotic signal X 2 of the master oscillator is applied to the slave oscillator's signal Y 2 through the synchronization circuit, which adds the n 2 noise signal to the applied X 2 signal. The design of the synchronization circuit is identical to the one used in the simulation setup in Figure 4. The Analog Discovery 2 (later referred to as AD2) devices were used in this experiment, labeled as AD2-X and AD2-Y [17].

Hardware-Based Study
This section presents the hardware implementation (see Figure 5) of the experimental model shown in Figure 3. For the hardware experiment, the chaos oscillators are implemented in specially-designed PCBs. A pair of chaos oscillators is powered by two separate independent power supplies. The chaotic signal 2 X of the master oscillator is applied to the slave oscillator's signal 2 Y through the synchronization circuit, which adds the 2 n noise signal to the applied 2 X signal. The design of the synchronization circuit is identical to the one used in the simulation setup in Figure 4. The Analog Discovery 2 (later referred to as AD2) devices were used in this experiment, labeled as AD2-X and AD2-Y [17]. AD2-X is employed to generate the 1 n noise signal that is added to the 1 X (W1), to apply power to the synchronization circuit (V+ and V-), and to record the 1 X and 1 Y voltage signals (Ch1 and Ch2). The AD2-Y purpose is to generate the 2 n noise signal that is added to the 2 Y signal in the synchronization circuit (W1) and to record the 1 n noise signal generated by the AD2-X (Ch1). To ensure the synchronous performance of the two AD2 devices, a dedicated application was created and run on the host computer, providing the execution of the mentioned tasks by AD2-X and AD-Y [17]. MATLAB software was used to sum the recorded 1 n and 1 X signals, resulting in the 1 1 X n + signal, as well as to calculate the 100 μs window correlation coefficient between the 1 1 X n + and the 1 Y signals. The noise pair 1 n , 2 n signal levels were measured from 20 dB to 0 dB SNR with a step of 5 dB. It is also important to notice that the 1 X , 2 X , 1 Y , 2 Y , for each chaos oscillator were the same as in the case of the simulation study.

Result Analysis
The current section summarizes and analyzes the results acquired in both simulationbased ( Figure 4) and hardware-based ( Figure 5) studies applied to RC1, RC2, and Vilnius chaos oscillators. The results for corresponding oscillators are presented in Figures 6-8.  AD2-X is employed to generate the n 1 noise signal that is added to the X 1 (W1), to apply power to the synchronization circuit (V+ and V−), and to record the X 1 and Y 1 voltage signals (Ch1 and Ch2). The AD2-Y purpose is to generate the n 2 noise signal that is added to the Y 2 signal in the synchronization circuit (W1) and to record the n 1 noise signal generated by the AD2-X (Ch1). To ensure the synchronous performance of the two AD2 devices, a dedicated application was created and run on the host computer, providing the execution of the mentioned tasks by AD2-X and AD-Y [17]. MATLAB software was used to sum the recorded n 1 and X 1 signals, resulting in the X 1 + n 1 signal, as well as to calculate the 100 µs window correlation coefficient between the X 1 + n 1 and the Y 1 signals. The noise pair n 1 , n 2 signal levels were measured from 20 dB to 0 dB SNR with a step of 5 dB. It is also important to notice that the X 1 , X 2 , Y 1 , Y 2 , for each chaos oscillator were the same as in the case of the simulation study.

Result Analysis
The current section summarizes and analyzes the results acquired in both simulationbased ( Figure 4) and hardware-based ( Figure 5) studies applied to RC1, RC2, and Vilnius chaos oscillators. The results for corresponding oscillators are presented in Figures 6-8.       In Figures 6-8, the 100 μs window correlation coefficient  over time sinc ginning of synchronization is shown for different signal-to-noise ratios of the noi pair 1 n and 2 n . The first value of  is at 100 μs, as  has a 100 μs long window in all of the mentioned figures, in general, the achievable values of  tended to In Figures 6-8, the 100 µs window correlation coefficient β over time since the beginning of synchronization is shown for different signal-to-noise ratios of the noise signal pair n 1 and n 2 . The first value of β is at 100 µs, as β has a 100 µs long window. As seen in all of the mentioned figures, in general, the achievable values of β tended to increase with the increasing of SNR, approaching perfect synchronization indicated by β = 1. This showed that lower noise levels allow better synchronization between the two chaos oscillators. If we compare the results of the three chaos oscillators only by the maximum achieved value of β at specific SNR levels, then the RC2 chaos oscillator was superior from this point-proven by peak values of β at SNR levels 0 dB and 5 dB. However, in Figures 6-8, β curves fluctuate in simulation and hardware-based studies. These fluctuations are caused by the sufficiently high level of noise and the stability of oscillators themselves. The amplitude of fluctuations shows that some oscillators are more susceptible to noise. Vilnius chaos oscillator is superior from this point, as Figure 8 shows the least fluctuations. It is important to consider both: the maximum achievable value of β at a certain SNR level, as well as the fluctuations of the β curve at this SNR level.
By comparing the simulation results with those acquired in the hardware-based study for each chaos oscillator, the correlation among the results was obvious. Some minor differences in β curves for the same SNR levels in simulation and hardware were caused by the tolerances of the elements used in the hardware experiments and the incomplete models of the elements used in the simulation.
Observing the simulation and hardware-based results individually for each chaos oscillator, in the case of the RC1 chaos oscillator, the hardware results for SNR of 0 and 5 dB in Figure 6 are better than the simulation results for the corresponding SNR levels. The explanation is that for the low SNR levels, the fluctuation of the β curves is greater, which means the initial state of the system before the synchronization plays a crucial part in synchronization. In the case of the RC2 chaos oscillator in Figure 7, the fluctuation of the β curves for the SNR levels of 0 dB and 5 dB is even more severe. Even for the SNR levels 10-20 dB, the fluctuation is noticeable in simulation results and more explicitly in the hardware results. In the case of the Vilnius chaos oscillator in Figure 8, the fluctuations of the β curves are the least in comparison to all the chaos oscillators under research, which also means better synchronization noise immunity than for the other two chaos oscillators.
The conclusions from the results acquired in this study state that the level of synchronization of the chaos oscillators is decent even for the low signal to noise ratio and that the use of 100 µs as a time window for the correlation coefficient is sufficient for the determination of the correlation between two chaotic signals. None of the selected chaos

Simulation-Based Study
The experimental model, presented in Figure 9, utilizes the baseband model of the communication system shown in Figure 1, excluding the quadrature modulator and demodulator. The noise being permanently present for 1PM X and 2 X , assuming the worstcase scenario. 1PM X represents state variable 1 X manipulated by ( ) m t message signal.
When the "1" bit is transmitted, , and when the "0" bit is transmitted, bit to noise power spectral density ratio (Eb/N0) levels. The frequency bandwidth of the noise signals 1 n and 2 n correspond to those of the signals 1PM X and 2 X (see Figure   3). The state variables 1 X , 2 X , 1 Y , and 2 Y were the same for each chaos oscillator like in the previous experimental study. The propagation of both signals was simulated in the LTspice circuit simulator, whereas MATLAB was utilized for the calculations and analysis of the results.

Hardware-Based Study
This section presents the hardware implementation (see Figure 10.) of the experimental model shown in Figure 9. The model is similar to the one shown in Figure 5 with some additions. A hardware binary phase shifter was introduced to manipulate the phase of the 1 X signal. AD2-X had an additional task of generating and applying information signals ( ) m t to the chaos binary phase shifter (W2). AD2-X records the 2 X phase manipulated and 1 Y voltage signals (Ch1 and Ch2). AD2-Y records the 1 n noise signal that For this experiment, the length of the message signal m(t) was 10 5 bits. The length of the correlation coefficient's β window was 100 µs. The correlation coefficient's β values are between −1 and 1; β = 1 when X 1 is received, and β = −1 when −X 1 is received. The correlation coefficient is used to decode the message signal, the decoded received message signal m(t) = "1" if β > 0, and m(t) = "0" if β < 0. The performance of the system was estimated by calculating the bit error ratio (BER) for the transmitted message signal m(t) and the decoded received message signal m(t) for several energy per bit to noise power spectral density ratio (Eb/N0) levels. The frequency bandwidth of the noise signals n 1 and n 2 correspond to those of the signals X 1PM and X 2 (see Figure 3). The state variables X 1 , X 2 , Y 1 , and Y 2 were the same for each chaos oscillator like in the previous experimental study. The propagation of both signals was simulated in the LTspice circuit simulator, whereas MATLAB was utilized for the calculations and analysis of the results.

Hardware-Based Study
This section presents the hardware implementation (see Figure 10.) of the experimental model shown in Figure 9. The model is similar to the one shown in Figure 5 with some additions. A hardware binary phase shifter was introduced to manipulate the phase of the X 1 signal. AD2-X had an additional task of generating and applying information signals m(t) to the chaos binary phase shifter (W2). AD2-X records the X 2 phase manipulated and Y 1 voltage signals (Ch1 and Ch2). AD2-Y records the n 1 noise signal that is generated by the AD2-X (Ch1), similarly as in the previous hardware experiment. The binary phase shifter's circuit used for the hardware experiments is shown ure 11. The phase manipulation was based on the use of the ADG1219 analog sw cuit. Elements in Figure 11 have the following nominals:  This experiment was based on the hardware-in-the-loop approach since MATLAB was employed to sum the recorded n 1 and X 1PM signals, giving the X 1PM + n 1 signal. Moreover, MATLAB was employed to decode the message signal m(t) by applying a 100 µs window correlation coefficient between the X 1PM + n 1 and Y 1 signals. BER is acquired by comparing decoded m(t) and original m(t) message signals. BERs were measured for all three chaos oscillators up to 11 dB Eb/N0, due to the hardware limitations.
The binary phase shifter's circuit used for the hardware experiments is shown in Figure 11. The phase manipulation was based on the use of the ADG1219 analog switch circuit. Elements in Figure 11 have the following nominals: R 1 = R 2 = R 3 = 10 kΩ; C 1 = 10 nF; OA 1 -OA 3 are TL082 operational amplifiers. The operational amplifiers were utilized to acquire an inverted and direct X 1 . C 1 and R 1 form a high-pass filter that removes the DC component from X 1 signal. OA 1 is employed as a buffer to isolate the high-pass filter from the chaos oscillator circuit. ADG1219 serves as a switch between inverted and direct X 1 according to the m(t); thus, the X 1PM signal forms on the output of the ADG1219. OA is employed as a buffer to isolate the high-pass filter from the chaos oscillator circuit. ADG1219 serves as a switch between inverted and direct 1 X according to the ( ) m t ; thus, the 1PM X signal forms on the output of the ADG1219.

Result Analysis
This section is dedicated to the analysis of the results acquired in simulation and hardware experimental studies. The results are presented in Figure 12 for the proposed Figure 11. Chaos binary phase shifter.

Result Analysis
This section is dedicated to the analysis of the results acquired in simulation and hardware experimental studies. The results are presented in Figure 12 for the proposed communication system's baseband experimental study of the model in the AWGN channel.  In Figure 12, the BERs over energy per bit to noise power spectral density ratio Eb/N0 are shown for all three selected chaos oscillators employed as chaotic X and Y systems. Figure 12 reflects superior simulation results of the system using the RC2 chaos oscillators compared to the other two chaos oscillators. For Eb/N0 12-17 dB, the system based on the RC1 chaos oscillators demonstrated better results than in the case of the other two chaos oscillators. At the point of Eb/N0 around 22 dB, the RC1 and RC2 based systems provided error-free communication. Overall, the RC2 chaos oscillator demonstrated the best performance in the simulation.
The results in hardware naturally demonstrated a greater number of error bits than the circuit simulation results, which is expected, as hardware implementation of chaos oscillators and other functional blocks introduced many parasitic parameters and imprecisions, compared to the simulation. It is important to notice that hardware results showed greater BER where simulation demonstrated error-free communication. However, the hardware results follow that the RC2 chaos oscillator demonstrated the worst performance among the three chaos oscillators, whereas the RC1 chaos oscillator had superior performance compared to other ones. Results in Figures 6-8 show that the lower SNR allows better synchronization, which in Figure 12 corresponds to a lower number of error bits and better BER (SNR and Eb/N0 are closely related). As noise level increases, the quality of synchronization decreases-depicted by lower  values and increased fluctuations in Figures 6-8. This corresponds to a greater number of error bits and worse BER in Figure  12. The results of the two studies correlate and give a better understanding of how the In Figure 12, the BERs over energy per bit to noise power spectral density ratio Eb/N0 are shown for all three selected chaos oscillators employed as chaotic X and Y systems. Figure 12 reflects superior simulation results of the system using the RC2 chaos oscillators compared to the other two chaos oscillators. For Eb/N0 12-17 dB, the system based on the RC1 chaos oscillators demonstrated better results than in the case of the other two chaos oscillators. At the point of Eb/N0 around 22 dB, the RC1 and RC2 based systems provided error-free communication. Overall, the RC2 chaos oscillator demonstrated the best performance in the simulation.
The results in hardware naturally demonstrated a greater number of error bits than the circuit simulation results, which is expected, as hardware implementation of chaos oscillators and other functional blocks introduced many parasitic parameters and imprecisions, compared to the simulation. It is important to notice that hardware results showed greater BER where simulation demonstrated error-free communication. However, the hardware results follow that the RC2 chaos oscillator demonstrated the worst performance among the three chaos oscillators, whereas the RC1 chaos oscillator had superior performance compared to other ones. Results in Figures 6-8 show that the lower SNR allows better synchronization, which in Figure 12 corresponds to a lower number of error bits and better BER (SNR and Eb/N0 are closely related). As noise level increases, the quality of synchronization decreases-depicted by lower β values and increased fluctuations in Figures 6-8. This corresponds to a greater number of error bits and worse BER in Figure 12. The results of the two studies correlate and give a better understanding of how the noise impacts the proposed CQPSK system.

Conclusions
Research done in the current paper was aimed at investigating the impact of the chaotic synchronization on the performance of the quadrature transceiver-based communication system that utilizes chaotic signals for data transmission. For this study, we chose RC1, RC2, and Vilnius chaos oscillators and studied their synchronization noise immunity as a first step of the research. The results have shown that the synchronization's noise immunity properties of the chaos oscillators are sufficient in potential employment for data transmission in low-complexity secure communication systems. We performed a comparison between the communication system hardware implementation with simulation to estimate the influence of the chaos oscillator's properties on the performance of the quadrature transceiver-based communication system in the AWGN channel. The results in this step revealed notable differences between simulation and hardware BERs. The employment of the 3 different chaos oscillators yielded comparable results, with the RC2 chaos oscillator showing superior performance in simulation, whereas the RC1 chaos oscillator showed superior performance in hardware.