Energy and Spectrally Efficient Modulation Scheme for IoT Applications

Due to the Internet of Things (IoT) requirements for a high-density network with low-cost and low-power physical (PHY) layer design, the low-power budget transceiver systems have drawn momentous attention lately owing to their superior performance enhancement in both energy efficiency and hardware complexity reduction. As the power budget of the classical transceivers is envisioned by using inefficient linear power amplifiers (PAs) at the transmitter (TX) side and by applying high-resolution analog to digital converters (ADCs) at the receiver (RX) side, the transceiver architectures with low-cost PHY layer design (i.e., nonlinear PA at the TX and one-bit ADC at the RX) are mandated to cope with the vast IoT applications. Therefore, in this paper, we propose the orthogonal shaping pulses minimum shift keying (OSP-MSK) as a multiple-input multiple-output (MIMO) modulation/demodulation scheme in order to design the low-cost transceiver architectures associated with the IoT devices. The OSP-MSK fulfills a low-power budget by using constant envelope modulation (CEM) techniques at the TX side, and by applying a low-resolution one-bit ADC at the RX side. Furthermore, the OSP-MSK provides a higher spectral efficiency compared to the recently introduced MIMO-CEM with the one-bit ADC. In this context, the orthogonality between the in-phase and quadrature-phase components of the OSP are exploited to increase the number of transmitted bits per symbol (bps) without the need for extra bandwidth. The performance of the proposed scheme is investigated analytically and via Monte Carlo simulations. For the mathematical analysis, we derive closed-form expressions for assessing the average bit error rate (ABER) performance of the OSP-MSK modulation in conjunction with Rayleigh and Nakagami-m fading channels. Moreover, a closed-form expression for evaluating the power spectral density (PSD) of the proposed scheme is obtained as well. The simulation results corroborate the potency of the conducted analysis by revealing a high consistency with the obtained analytical formulas.


Introduction
Recently, the dawn of the Internet of Things (IoT) has intensified the research efforts to fulfill the escalating demand for low-power budget wireless transceiver systems. Overnight, the communication industry experienced an explosion in the number of sensors and terminals that connect mobile users to the internet [1,2]. However, low-cost, low-power, and low-bitrate sensors and terminals are mandated

•
The OSP-MSK Modulation/Demodulation Scheme: The OSP-MSK is proposed as a low-power and low-cost modulation scheme. The OSP-MSK relaxes both the power budget and the hardware complexity by using nonlinear PAs at the TX side and by deploying a low-resolution one-bit ADC at the RX side. Furthermore, the OSP-MSK achieves a higher spectral efficiency than the recent MIMO-CEM system by increasing the number of bits per symbol without additional bandwidth. • OSP-MSK Performance Analysis: Firstly, a generic closed form expression for assessing the average bit error rate (ABER) performance of the proposed OSP-MSK modulation over an additive white Gaussian noise (AWGN) is presented thoroughly. Afterwards, this expression is used to derive closed-form expressions for the ABER performance of the OSP-MSK over the generalized Rayleigh and Nakagami-m fading channels. • OSP-MSK Bandwidth Utilization: The bandwidth utilization of the OSP-MSK is evaluated by deriving a closed form expression of its power spectral density (PSD). This PSD is used to assess the ability of the OSP-MSK to boost the spectral efficiency compared to the MIMO-CEM system.
The rest of this paper is structured as follows. Section 2 addresses the system model of the proposed OSP-MSK modulation. The performance analysis of the OSP-MSK is driven in Section 3. The simulation results are introduced in Section 4. Finally, the paper is concluded in Section 5.
Notation: Throughout this paper, vectors are represented using bold lowercase letters and matrices are represented using bold capital letters. Complex Gaussian distributed random variables (RVs) are expressed by CN µ, σ o 2 , where µ and σ o 2 denote the mean and the variance, respectively. We use p Ψ (Ψ) to refer to the probability density function (PDF) of a random variable Ψ. Moreover, the notations log (.), |.|, and (.) are used to indicate the natural logarithm, the magnitude, and the , and Γ(.) are used to indicate the Q-function, the complementary error function, and the gamma function, respectively. Finally, we use 2F 1 (a, b; c; z) to indicate the Gaussian hypergeometric function, where a, b, and c are Gaussian hypergeometric function arbitrary parameters, and z is the function argument. Figure 1 depicts the MIMO orthogonal shaping pulses minimum shift keying (OSP-MSK) modulation system model. As shown in the figure, the input binary data are firstly encoded by using a convolutional encoder accompanied with an interleaver to make the transmitted signal of the OSP-MSK more robust against the prospective disruptions (i.e., noise, interference, multipath propagation, etc.). The convolutionally encoded data are then partitioned into two distinct sets, and each dataset is divided into in-phase components (i.e., even bits) and quadrature-phase components (i.e., odd bits). These components are then modulated using constant envelope MSK modulation, which results in a signal with a unitary peak to average power ratio (PAPR). Therefore, an energy efficient nonlinear PA is utilized to amplify each modulated dataset. Hence, a superior enhancement in the power consumption is achieved by using the OSP-MSK modulation [26]. each dataset is divided into in-phase components (i.e., even bits) and quadrature-phase components (i.e., odd bits). These components are then modulated using constant envelope MSK modulation, which results in a signal with a unitary peak to average power ratio (PAPR). Therefore, an energy efficient nonlinear PA is utilized to amplify each modulated dataset. Hence, a superior enhancement in the power consumption is achieved by using the OSP-MSK modulation [26].

The OSP-MSK Modulation System Model
As such, the constant envelope signal of each dataset is given by: which can be rewritten as follows: where is the gain of the first amplifier, is the energy per bit, is the bit period, and is the carrier frequency.
( ) and (t) with = 1, 2 are the in-phase and the quadrature-phase components of each MSK signal, respectively. In the same way, ( ) is given by: where denotes the gain of the second amplifier. Therefore, the transmitted signal of the proposed OSP-MSK is expressed as follows: (5) As such, the constant envelope signal of each dataset is given by: which can be rewritten as follows: where A 1 is the gain of the first amplifier, E b is the energy per bit, T is the bit period, and f c is the carrier frequency. S iI (t) and S iQ (t) with i = 1, 2 are the in-phase and the quadrature-phase components of each MSK signal, respectively. In the same way, S 2 (t) is given by: where A 2 denotes the gain of the second amplifier. Therefore, the transmitted signal of the proposed OSP-MSK is expressed as follows: Confusion may be arisen here that both S 1 (t) and S 2 (t) are two constant envelope signals (i.e., unitary PAPR signals), but their sum (i.e., S(t)) is not a constant envelope signal. The OSP-MSK transmitted signal S(t) is actually not a constant envelope signal, but it is a low PAPR signal (see Appendix C), because either S 1 (t) or S 2 (t) is amplified separately before it is combined with the other signal. This is in contrast to the MIMO-OFDM system, where all of the subcarriers are combined (summed) before amplifying [27]. Hence, the low PAPR property of the OSP-MSK transmitted signal enables the one-bit ADC to be beneficially used with the OSP-MSK at the RX side.
On its path, the OSP-MSK transmitted signal undergoes a fading channel, and it is corrupted by AWGN. At the RX side, a coherent demodulation is applied to the received signal, as shown in Figure 1. To extract the MSK in-phase components S 1I and S 2I , the OSP-MSK received signal S(t) given by Equation (5) is firstly multiplied by cos(2π f c t). Then, a low pass filter (LPF) is applied to the resulting signal in order to detect the summation of the two in-phase components S I as follows: Finally, the MSK in-phase components are separated by exploiting the orthogonality between the shaping pulses i.e., cos πt 2T and cos . By multiplying Equation (6) by 2 T cos πt 2T and applying the integration of the result over 2T, the second term in Equation (6) will vanish. Thus, the first in-phase component S 1I can be easily extracted by using a low resolution one-bit ADC. The same procedure is applied to S 2I , but Equation (6) is multiplied by 2 T cos instead of 2 T cos πt 2T . In the same way, the MSK quadrature-phase components S 1Q and S 2Q can be extracted by multiplying the OSP-MSK received signal S(t) by sin(2π f c t). Then, a filtering process is applied to detect the summation of the two quadrature-phase components S Q as follows: The quadrature-phase component of each dataset is then obtained by the same orthogonality principle used above. The extracted MSK in-phase and quadrature-phase components are then parallelly to serially converted. Then, a convolutional decoder and deinterleaver are applied in pursuit of estimating the input binary data.
It should be noted here that the use of low-resolution one-bit ADCs to detect the in-phase and the quadrature-phase components at the RX side of the OSP-MSK removes most of the analog stages (e.g., automatic gain control circuit (AGC), analog filters, etc.). Moreover, the use of the coherent demodulation along with the one-bit ADCs at the RX side of the OSP-MSK eliminates the need for the complex maximum-likelihood sequence estimation (MLSE) equalizer mandated at the RX side of the MIMO-CEM system [21][22][23][24][25]. Hence, with the OSP-MSK modulation, a worthwhile reduction in both the power consumption and the hardware complexity is achieved over the MIMO-CEM system.

OSP-MSK Performance Analysis
In this section, the analytical average bit error rate (ABER) performance of the proposed OSP-MSK modulation is derived. Furthermore, to evaluate the spectral efficiency of the OSP-MSK, its analytical power spectral density (PSD) is obtained in detail.

OSP-MSK ABER Perfrmance
The ABER performance of the OSP-MSK modulation relies mainly on the correlation receiver (i.e., four-arms correlation receiver), as shown in Figure 1. For further illustration, one arm correlation receiver of the OSP-MSK demodulator is shown in Figure 2. It is assumed here that the received signals S I (t) and S Q (t) at the inputs of the four-arms correlator receiver are corrupted by AWGN distributed as n(t) ∼ CN 0, N o 2 .

OSP-MSK ABER Perfrmance
The ABER performance of the OSP-MSK modulation relies mainly on the correlation receiver (i.e., four-arms correlation receiver), as shown in Figure 1. For further illustration, one arm correlation receiver of the OSP-MSK demodulator is shown in Figure 2. It is assumed here that the received signals ( ) and ( ) at the inputs of the four-arms correlator receiver are corrupted by AWGN distributed as ( )⁓ (0, ). Accordingly, the output of the sampler can be expressed as follows: (2 ) = ( ) ( ) + ( ) ( ) .
However, in Equation (11), is a deterministic signal, but is an AWGN distributed as ⁓ (0, ). Therefore, the output of the sampler (2 ) provided by Equation (11) can be rewritten as follows: (2 ) = + ℎ "1" and the decision threshold can be expressed as follows: Hence, to correctly recover "0" when is transmitted, the term − (2 ) should be positive. That is: Therefore, the noise should satisfy the following inequality: and the error will occur when the noise satisfies the following inequality: Accordingly, the probability of error (i.e., the error in detecting "0" when is transmitted) is given by: Accordingly, the output of the sampler can be expressed as follows: However, in Equation (11), S o is a deterministic signal, but n o is an AWGN distributed as n o ∼ CN 0, N o 2 . Therefore, the output of the sampler r I (2T) provided by Equation (11) can be rewritten as follows: and the decision threshold ξ can be expressed as follows: Hence, to correctly recover "0" when S o2 is transmitted, the term ξ − r I (2T) should be positive. That is: Therefore, the noise n o should satisfy the following inequality: and the error will occur when the noise n o satisfies the following inequality: Accordingly, the probability of error (i.e., the error in detecting "0" when S o2 is transmitted) is given by: where σ o 2 is the noise variance, which is equal to N o 2 . By letting r = q , the probability of error can be expressed as follows: However, in the proposed OSP-MSK, S o in Equation (11) can be obtained by letting φ(t) as follows: substituting Equation (21) in Equation (11), S o can be expressed as follows: Hence, S o will take one of two values as follows: Substituting Equation (23) in Equation (20), the probability of error can be expressed as follows: The probability of error in Equation (24) is the probability of error of the first arm correlation receiver of the OSP-MSK demodulator. Likewise, the probability of error of the second arm correlation receiver can be acquired, and it is equal to the probability of error of the first arm. However, it should be mentioned here that the second arm correlation receiver constitutes a quadrature multiplexed version of the first arm. Therefore, the probability of correctness due to the two arms can be expressed as follows: Hence, the probability of error due to the first two arms can be expressed as follows: Likewise, the probability of error due to the other two arms p e 3,4 can be acquired, and it equals to the probability of error due to the first two arms, p e 1,2 .
Lastly, the OSP-MSK demodulator can be interpreted as two independent correlation receivers, each with a probability of error of p e 1,2 and p e 3,4 , respectively. Therefore, the ABER of the proposed OSP-MSK over AWGN can be expressed as follows: To obtain the ABER of the proposed OSP-MSK modulation over different fading channels, the previous analysis can be expanded by averaging the ABER of AWGN provided by Equation (27) over the PDF of fading channels [28]. More specifically, when the fading is present, the received signal amplitude is attenuated by the fading amplitude ε, and the instantaneous signal power will be attenuated by ε 2 . It is should be noted here that ε is a random variable (RV) with a mean square value β = ε 2 , and its PDF depends upon the type of fading channel. As such, the signal to noise ratio (SNR) per bit in the presence of fading can be expressed as follows: and the average SNR per bit can be expressed as follows: Therefore, the ABER performance of the OSP-MSK over any fading channel can be determined by integrating the channel fading effect in Equation (28) in the formula of the ABER over AWGN in Equation (27) to acquire the conditional ABER (i.e., conditioned on the channel). Then, the ABER due to the presence of fading can be expressed as follows: where ABER E|H denotes the conditional ABER, and p Ψ (Ψ) is the PDF of the fading channel.

Rayleigh Fading Channel
Following the aforementioned procedure, the conditional ABER E|H over the Rayleigh fading channel of the proposed OSP-MSK can be obtained by replacing the term E b /N o in Equation (27) by Ψ in Equation (28) as follows: Substituting Equation (31) into Equation (30), the ABER of the OSP-MSK over Rayleigh fading channel (ABER Ray ) can be expressed as follows: where p Ψ (Ψ) is the PDF of the Rayleigh fading channel. Using the well-known PDF of the Rayleigh fading, the ABER Ray can be expressed as follows: Integral I 1 in Equation (33) is determined in detail in Appendix A, and it can be expressed as follows: However, the second integral I 2 in Equation (33) can be written in the following form: which can be easily solved as in [29] as follows: which can be easily solved as in [29] as follows: where ℽ is given as follows: Hence, the second integral 2 in Equation (33) is determined by substituting 1 = 1 and 2 = 1/ᴪ ̅ in Equations (36) and (37) as follows: Substituting Equations (34) and (38) into Equation (33), the ABER over the Rayleigh fadin channel can be expressed as follows:

Nakagami-M Fading Channel
Using the same procedure, the ABER in the case of the Nakagami-m fading channel can b expressed as follows: Integral 3 in Equation (40) is determined in detail in Appendix B, and it can be expressed as follows: where ( ᴪ ̅ ) can be expressed as follows: Integral 4 in Equation (40) is in the form of the integral in Equation (35), but by substitutin 1 = and 2 = ᴪ ̅ . Therefore, it can be determined as in Equation (36) as follows: which can be easily solved as in [29] as follows: where ℽ is given as follows: Hence, the second integral 2 in Equation (33) is determined by substituting 1 = 1 and 1/ᴪ ̅ in Equations (36) and (37) as follows: Substituting Equations (34) and (38) into Equation (33), the ABER over the Rayleigh channel can be expressed as follows:

Nakagami-M Fading Channel
Using the same procedure, the ABER in the case of the Nakagami-m fading channel expressed as follows: (40) is determined in detail in Appendix B, and it can be exp as follows: ) can be expressed as follows: (40) is in the form of the integral in Equation (35), but by subst . Therefore, it can be determined as in Equation (36) as follows: where Sensors 2018, 18, x FOR PEER REVIEW 9 of 20 which can be easily solved as in [29] as follows: where ℽ is given as follows: Hence, the second integral 2 in Equation (33) is determined by substituting 1 = 1 and 2 = 1/ᴪ ̅ in Equations (36) and (37) as follows: Substituting Equations (34) and (38) into Equation (33), the ABER over the Rayleigh fading channel can be expressed as follows:

Nakagami-M Fading Channel
Using the same procedure, the ABER in the case of the Nakagami-m fading channel can be expressed as follows: Integral 3 in Equation (40) is determined in detail in Appendix B, and it can be expressed as follows: where ( ᴪ ̅ ) can be expressed as follows: Integral 4 in Equation (40) is in the form of the integral in Equation (35), but by substituting . Therefore, it can be determined as in Equation (36) as follows: is given as follows: Sensors 2018, 18, x FOR PEER REVIEW which can be easily solved as in [29] as follows: where ℽ is given as follows: Hence, the second integral 2 in Equation (33) is determined by substituting 1/ᴪ ̅ in Equations (36) and (37) as follows: Substituting Equations (34) and (38) into Equation (33), the ABER over the channel can be expressed as follows:

Nakagami-M Fading Channel
Using the same procedure, the ABER in the case of the Nakagami-m fadin expressed as follows:

Integral
3 in Equation (40) is determined in detail in Appendix B, and it as follows: ) can be expressed as follows: Hence, the second integral I 2 in Equation (33) is determined by substituting k 1 = 1 and k 2 = 1/Ψ in Equations (36) and (37) as follows: Substituting Equations (34) and (38) into Equation (33), the ABER over the Rayleigh fading channel can be expressed as follows:

Nakagami-M Fading Channel
Using the same procedure, the ABER in the case of the Nakagami-m fading channel can be expressed as follows: Integral I 3 in Equation (40) is determined in detail in Appendix B, and it can be expressed as follows: where U Ψ m can be expressed as follows: Integral I 4 in Equation (40) is in the form of the integral A in Equation (35), but by substituting k 1 = m and k 2 = m Ψ . Therefore, it can be determined as in Equation (36) as follows: where Z is given as follows: Substituting Equations (41) and (43) into Equation (40), the ABER over the Nakagami-m fading channel can be expressed as follows: when m is non intger (45)

OSP-MSK Power Spectral Density
In the following, a closed form expression for the PSD of the OSP-MSK is analytically derived to evaluate the effectiveness of the proposed OSP-MSK in boosting the spectral efficiency compared to the MIMO-CEM. As depicted in Figure 1, the only main difference between the OSP-MSK and the conventional MSK is the second branch (i.e., S 2 (t)) in Equation (4). Thus, the PSD of the OSP-MSK mainly relies on Equation (4). As, cos = cos πt 2T , then Equation (4) can be rewritten as follows: Therefore, the spectrum of S 1 (t) or S 2 (t) mainly relies on the bandwidth occupancy of the orthogonal shaping pulses cos πt 2T and sin πt 2T . Accordingly, the PSD of the OSP-MSK can be obtained by deriving the PSD of the complex envelope representation of either S 1 (t) provided by Equation (2) or S 2 provided by Equation (4) [30]. As such, to obtain the PSD of S 1 (t), we can rewrite Equation (2) as follows: where S even (t) and S odd (t) are the in-phase and quadrature-phase components, respectively. Therefore, the complex representation of S 1 (t) can be expressed as follows: However, S even (t) and S odd (t) can be expressed as follows: where S 1I k (t), and S 1Q k (t) are the in-phase and quadrature-phase bits, respectively, and g(t) is the half cycle sine shaping pulse, which can be expressed as follows: Since the in-phase and quadrature-phase components of the OSP-MSK signal are independent and the PSDs of S even (t) and S odd (t) are equal, the baseband PSD of S 1 (t) can be expressed as follows: where G( f ) is the Fourier transform of the half-cycle shaping pulse g(t). However, the magnitude squared of G( f ) can be expressed as follows: Substituting Equation (53) into Equation (52), the baseband PSD of S 1 (t) can be expressed as follows: Using some mathematical manipulations, the baseband PSD of S 1 (t) can be expressed as follows: (55)

Simulation Results
Firstly, the ABER performances of the proposed OSP-MSK over AWGN, Rayleigh fading channel, and Nakagami-m fading channel are evaluated analytically and via comprehensive Monte Carlo simulations. To compare meaningfully, the analytical formulas and the simulation results are assessed in both uncoded and coded scenarios. It should be mentioned here that the relation between the code symbol energy per N o and the bit energy per N o is expressed as follows [31]: in which n k denotes the channel encoder rate. Secondly, the ABER performance of the proposed MIMO OSP-MSK modulation is evaluated and compared to the ABER performance of the MIMO-CEM system using the traditional GMSK modulation, and by using the quadrature multiplexed Gaussian minimum shift keying (QMGMSK) modulation [32,33]. The QMGMSK modulation is a special form of the MSK modulation that boosts up the spectral efficiency of the MIMO-CEM by reducing the main lobe spectrum (i.e., the bandwidth) of the MSK modulation. The ABER comparison is performed over the Rayleigh fading channel and under different scenarios based upon the simulation parameters listed in Table 1.
Finally, the PSD of the proposed OSP-MSK modulation is evaluated analytically and via numerical simulation and compared with the PSDs of the above-mentioned schemes in order to evaluate the spectral capabilities of the proposed OSP-MSK modulation. (27) compared to the numerical simulation results for the uncoded and coded scenarios. As shown in Figure 3, the simulation results manifest a close match to the analytical formula over a wide range of SNR values.  Figure 4 depicts the analytical ABER of the OSP-MSK modulation over the Rayleigh fading channel in Equation (39) and over the Nakagami-m fading channel in Equation (45) (e.g., m = two and m = four) compared to the numerical simulation results for the uncoded scenario. Furthermore, Figure 5 considers the coded scenario, where the convolutional encoder rate is set to 1 2 . As depicted in Figures 4 and 5, again, the simulation results are very akin to the analytical results for a wide range of SNR values.

Figure 3 depicts the analytical ABER of the proposed OSP-MSK over AWGN in Equation
The high consistency between the analytical formulas and the simulation results dates primarily to the absence of any mathematical assumptions in deriving the ABER of the OSP-MSK modulation. This, without a doubt, evinces the potency of the mathematical analysis conducted throughout this paper.

Channel decoder Hard decision Viterbi decoder Number of transmit antennas
Two antennas Number of receive antennas Two antennas Figure 4 depicts the analytical ABER of the OSP-MSK modulation over the Rayleigh fading channel in Equation (39) and over the Nakagami-m fading channel in Equation (45) (e.g., m = two and m = four) compared to the numerical simulation results for the uncoded scenario. Furthermore, Figure  5 considers the coded scenario, where the convolutional encoder rate is set to ½. As depicted in figures 4 and 5, again, the simulation results are very akin to the analytical results for a wide range of SNR values.
The high consistency between the analytical formulas and the simulation results dates primarily to the absence of any mathematical assumptions in deriving the ABER of the OSP-MSK modulation. This, without a doubt, evinces the potency of the mathematical analysis conducted throughout this paper.     Figure 6 depicts the ABER performance of the single-input single-output orthogonal shaping pulses minimum shift keying (SISO-OSP-MSK) modulation compared to the ABER performance of the SISO-CEM system by using GMSK and QMGMSK modulation schemes [32,33]. As depicted in Figure 6, the proposed SISO-OSP-MSK captures the same ABER performance of the conventional SISO-CEM with GMSK modulation, while the ABER performance of the SISO-CEM system is degraded by about four dB of SNR value when the QMGMSK modulation is utilized. Figure 7 depicts the ABER performance of the MIMO-OSP-MSK modulation compared to the ABER performance of the MIMO-CEM system with the GMSK and QMGMSK modulation schemes, where a 2×2 transmit and receive antenna configuration is utilized to achieve a multiplexing gain for both of the considered schemes. As depicted in Figure 7, a comparable result is obtained in the MIMO scenario. Again, the proposed MIMO-OSP-MSK captures the same ABER performance of the conventional MIMO-CEM with GMSK modulation, while the ABER performance of the MIMO-CEM system is degraded by about five dB of SNR value when the QMGMSK modulation is utilized.
In Figure 8, the ABER performance of the MIMO-OSP-MSK modulation is evaluated and compared to the ABER performance of the conventional MIMO-CEM system in the presence of the quantization errors (i.e., one-bit ADC is utilized) and in the absence of the quantization errors (i.e.,    Figure 6 depicts the ABER performance of the single-input single-output orthogonal shaping pulses minimum shift keying (SISO-OSP-MSK) modulation compared to the ABER performance of the SISO-CEM system by using GMSK and QMGMSK modulation schemes [32,33]. As depicted in Figure 6, the proposed SISO-OSP-MSK captures the same ABER performance of the conventional SISO-CEM with GMSK modulation, while the ABER performance of the SISO-CEM system is degraded by about four dB of SNR value when the QMGMSK modulation is utilized. Figure 7 depicts the ABER performance of the MIMO-OSP-MSK modulation compared to the ABER performance of the MIMO-CEM system with the GMSK and QMGMSK modulation schemes, where a 2×2 transmit and receive antenna configuration is utilized to achieve a multiplexing gain for both of the considered schemes. As depicted in Figure 7, a comparable result is obtained in the MIMO scenario. Again, the proposed MIMO-OSP-MSK captures the same ABER performance of the conventional MIMO-CEM with GMSK modulation, while the ABER performance of the MIMO-CEM system is degraded by about five dB of SNR value when the QMGMSK modulation is utilized.
In Figure 8, the ABER performance of the MIMO-OSP-MSK modulation is evaluated and compared to the ABER performance of the conventional MIMO-CEM system in the presence of the quantization errors (i.e., one-bit ADC is utilized) and in the absence of the quantization errors (i.e.,  Figure 6 depicts the ABER performance of the single-input single-output orthogonal shaping pulses minimum shift keying (SISO-OSP-MSK) modulation compared to the ABER performance of the SISO-CEM system by using GMSK and QMGMSK modulation schemes [32,33]. As depicted in Figure 6, the proposed SISO-OSP-MSK captures the same ABER performance of the conventional SISO-CEM with GMSK modulation, while the ABER performance of the SISO-CEM system is degraded by about four dB of SNR value when the QMGMSK modulation is utilized. Figure 7 depicts the ABER performance of the MIMO-OSP-MSK modulation compared to the ABER performance of the MIMO-CEM system with the GMSK and QMGMSK modulation schemes, where a 2 × 2 transmit and receive antenna configuration is utilized to achieve a multiplexing gain for both of the considered schemes. As depicted in Figure 7, a comparable result is obtained in the MIMO scenario. Again, the proposed MIMO-OSP-MSK captures the same ABER performance of the conventional MIMO-CEM with GMSK modulation, while the ABER performance of the MIMO-CEM system is degraded by about five dB of SNR value when the QMGMSK modulation is utilized.
In Figure 8, the ABER performance of the MIMO-OSP-MSK modulation is evaluated and compared to the ABER performance of the conventional MIMO-CEM system in the presence of the quantization errors (i.e., one-bit ADC is utilized) and in the absence of the quantization errors (i.e., infinite-bit ADC is utilized) for the uncoded and coded scenarios (e.g., soft Viterbi decoder is utilized at the RX side). It should be noted that when the infinite-bit ADC is utilized with the OSP-MSK modulation, the use of a MLSE equalizer is mandated at the RX side in order to efficiently recover the input binary data.
As shown in Figure 8, for the uncoded scenario, about an eight-dB difference in the SNR value is observed between the one-bit and the infinite-bit cases. While, for the coded scenario, the difference is about three dB. The prior results without a doubt manifest the effect of the harsh quantization error caused by the one-bit ADC, which severely degrades the performance of the proposed scheme as the conventional massive MIMO systems [11,12]. Furthermore, the results emphasize the importance of using channel coding techniques in both the MIMO-CEM and the proposed OSP-MSK modulation. infinite-bit ADC is utilized) for the uncoded and coded scenarios (e.g., soft Viterbi decoder is utilized at the RX side). It should be noted that when the infinite-bit ADC is utilized with the OSP-MSK modulation, the use of a MLSE equalizer is mandated at the RX side in order to efficiently recover the input binary data. As shown in Figure 8, for the uncoded scenario, about an eight-dB difference in the SNR value is observed between the one-bit and the infinite-bit cases. While, for the coded scenario, the difference is about three dB. The prior results without a doubt manifest the effect of the harsh quantization error caused by the one-bit ADC, which severely degrades the performance of the proposed scheme as the conventional massive MIMO systems [11,12]. Furthermore, the results emphasize the importance of using channel coding techniques in both the MIMO-CEM and the proposed OSP-MSK modulation.  infinite-bit ADC is utilized) for the uncoded and coded scenarios (e.g., soft Viterbi decoder is utilized at the RX side). It should be noted that when the infinite-bit ADC is utilized with the OSP-MSK modulation, the use of a MLSE equalizer is mandated at the RX side in order to efficiently recover the input binary data. As shown in Figure 8, for the uncoded scenario, about an eight-dB difference in the SNR value is observed between the one-bit and the infinite-bit cases. While, for the coded scenario, the difference is about three dB. The prior results without a doubt manifest the effect of the harsh quantization error caused by the one-bit ADC, which severely degrades the performance of the proposed scheme as the conventional massive MIMO systems [11,12]. Furthermore, the results emphasize the importance of using channel coding techniques in both the MIMO-CEM and the proposed OSP-MSK modulation.    Figure 9 depicts the analytical PSD of the OSP-MSK modulation compared to the PSD of the OSP-MSK conducted by the numerical simulation. Furthermore, the PSD of the OSP-MSK is compared with the PSDs of the different constant envelope modulation schemes (i.e., GMSK and QMGMSK) in order to corroborate the bandwidth effectiveness of the OSP-MSK modulation. As shown in Figure 9, the spectrum of the QMGMSK [32,33] occupies about half of the main lobe spectrum of the traditional GMSK modulation. Thus, the QMGMSK approach can be used to double the spectral efficiency of the MIMO-CEM. However, unfortunately, this occurs with deterioration in the ABER performance, as shown Figure 6 and Figure 7. This dates back to the inapplicability of the QMGMSK with its non-coherent differential demodulation with the one-bit ADC in the MIMO-CEM system [32,33]. However, the spectrum of the OSP-MSK modulation occupies the same main lobe spectrum as the GMSK modulation. Therefore, the proposed OSP-MSK increases the number of transmitted bits per symbol of the conventional MSK/GMSK from two bits/symbol to four bits/symbol without the need for additional bandwidth. This undoubtedly proves the effectiveness of the proposed OSP-MSK modulation in boosting the spectral efficiency over the MIMO-CEM system.  Figure 9 depicts the analytical PSD of the OSP-MSK modulation compared to the PSD of the OSP-MSK conducted by the numerical simulation. Furthermore, the PSD of the OSP-MSK is compared with the PSDs of the different constant envelope modulation schemes (i.e., GMSK and QMGMSK) in order to corroborate the bandwidth effectiveness of the OSP-MSK modulation.  Figure 9 depicts the analytical PSD of the OSP-MSK modulation compared to the PSD of the OSP-MSK conducted by the numerical simulation. Furthermore, the PSD of the OSP-MSK is compared with the PSDs of the different constant envelope modulation schemes (i.e., GMSK and QMGMSK) in order to corroborate the bandwidth effectiveness of the OSP-MSK modulation. As shown in Figure 9, the spectrum of the QMGMSK [32,33] occupies about half of the main lobe spectrum of the traditional GMSK modulation. Thus, the QMGMSK approach can be used to double the spectral efficiency of the MIMO-CEM. However, unfortunately, this occurs with deterioration in the ABER performance, as shown Figure 6 and Figure 7. This dates back to the inapplicability of the QMGMSK with its non-coherent differential demodulation with the one-bit ADC in the MIMO-CEM system [32,33]. However, the spectrum of the OSP-MSK modulation occupies the same main lobe spectrum as the GMSK modulation. Therefore, the proposed OSP-MSK increases the number of transmitted bits per symbol of the conventional MSK/GMSK from two bits/symbol to four bits/symbol without the need for additional bandwidth. This undoubtedly proves the effectiveness of the proposed OSP-MSK modulation in boosting the spectral efficiency over the MIMO-CEM system. As shown in Figure 9, the spectrum of the QMGMSK [32,33] occupies about half of the main lobe spectrum of the traditional GMSK modulation. Thus, the QMGMSK approach can be used to double the spectral efficiency of the MIMO-CEM. However, unfortunately, this occurs with deterioration in the ABER performance, as shown Figures 6 and 7. This dates back to the inapplicability of the QMGMSK with its non-coherent differential demodulation with the one-bit ADC in the MIMO-CEM system [32,33]. However, the spectrum of the OSP-MSK modulation occupies the same main lobe spectrum as the GMSK modulation. Therefore, the proposed OSP-MSK increases the number of transmitted bits per symbol of the conventional MSK/GMSK from two bits/symbol to four bits/symbol without the need for additional bandwidth. This undoubtedly proves the effectiveness of the proposed OSP-MSK modulation in boosting the spectral efficiency over the MIMO-CEM system.

Conclusions
In this paper, the multiple-input multiple-output orthogonal shaping pulses minimum shift keying (MIMO-OSP-MSK) has been presented as a low-cost modulation/demodulation scheme. The OSP-MSK fulfills a low-power budget by using a nonlinear PA at the TX side and by applying a low resolution one-bit ADC at the RX side. Moreover, the OSP-MSK harnesses in a new fashion the orthogonality between the shaping pulses to boost up the spectral efficiency. By means of this orthogonality, the OSP-MSK increases the number of transmitted bits per symbol of the MSK/GMSK without additional cost for extra bandwidth. Mathematical closed-form expressions for assessing the ABER performance of the OSP-MSK modulation over the AWGN channel, Rayleigh fading channel, and Nakagami-m fading channel were derived and analyzed thoroughly. In addition, a mathematical closed-form expression for the power spectral density of the proposed OSP-MSK was also derived and studied in detail. The simulation results corroborate the effectiveness of the conducted analysis by revealing a close match to the obtained analytical formulas. Furthermore, the simulation results evinced that the proposed OSP-MSK outweighs the MIMO-CEM with the one-bit ADC, in sense of the spectral efficiency and without any deterioration in the ABER performance.
letting u = tan(θ) = sin(θ) cos(θ) and using some mathematical manipulations, the integral I 1 can be expressed as follows: Using partial fraction, the integral I 1 can be rewritten as follows: Using formulas in [29], integral I 1 can be given as the formula in Equation (34).

Appendix B
Integral I 3 in Equation (40) can be rewritten by using the Q-function alternative as follows: But by letting C as follows: By using formulas in [29], the integral C can be expressed as follows: Substituting Equation (A12) into Equation (A10), integral I 3 can be expressed as follows: The integral I 3 in Equation (A13) is definite integral whose closed-form solution can be found in [28,34] as Equation (41).

Appendix C
The peak to average power ratio (PAPR) of the OSP-MSK transmitted signal S(t) can be expressed as follows: PAPR S(t) = 10log 10 P peak P avg = 10log 10 max{|S(t) where P peak and P avg are the peak and the average power of the OSP-MSK transmitted signal, respectively. However, S 2 (t) can be given as follows: S 2 (t) = a 2 + a 2 sin πt T sin(4π f c t) + a 2 2 sin πt T + a 2 2 sin(4π f c t), where a = 2E b T . Hence, to find the maximum of {|S(t)| 2 }, we let d dt S 2 (t) = 0. That is: d dt S 2 (t) = πa 2 T cos πt T sin(4π f c t) + 4πa 2 f c sin πt T cos(4π f c t) + πa 2 2T cos πt T + 2πa 2 f c sin πt T cos(4π f c t) = 0. (A16) For simplicity, let t = T 2 and f c = 1 4T ; therefore, the right and the left-hand sides of Equation (A16) are equalized. As such, the peak power of S(t) can be obtained as follows: However, the average power of S(t) is given by: which can be rewritten as follows: 2T 0 a 2 dt + (A20) Therefore, the average power P avg in Equation (A20) is given as follows: Substituting equations (A17) and (A21) into Equation (A14), the PAPR of the OSP-MSK transmitted signal is given by: PAPR S(t) = 10log 10 P peak P avg = 10log 10 (3) = 4.77dB.
The value of the PAPR of the transmitted signal S(t) provided by Equation (A22) is low compared to the high PAPR values of the OFDM system, which can be greater than 10 dB [27,35]. However, it should be mentioned here that the low PAPR property of the OSP-MSK signal doesn't contradict with the use of nonlinear PAs at the TX side of the OSP-MSK modulation, since both S 1 (t) and S 2 (t) are constant envelope signals (i.e., unitary PAPR signals).