Gyre Precoding and T-Transformation-Based GFDM System for UAV-Aided mMTC Network

: In this paper, an unmanned aerial vehicle (UAV)-aided multi-antenna conﬁgured downlink mmWave cooperative generalized frequency division multiplexing (GFDM) system is proposed. To provide physical layer security (PLS), a 3D controlled Lorenz mapping system is introduced. Furthermore, the combination of T-transformation spreading codes, walsh Hadamard transform, and discrete Fourier transform (DFT) techniques are integrated with a novel linear multi-user multiple-input multiple-output (MU-MIMO) gyre precoding (GP) for multi-user interference reduction. Furthermore, concatenated channel-coding with multi-user beamforming weighting-aided maximum-likelihood and zero forcing (ZF) signal detection schemes for an improved bit error rate (BER) are also used. The system is then simulated with a single base station (BS), eight massive machine-type communications (mMTC) users, and two UAV relay stations (RSs). Numerical results reveal the robustness of the proposed system in terms of PLS and an achievable ergodic rate with signal-to-interference-plus-noise ratio (SINR) under the implementation of T-transformation scheme. By incorporating the 3D mobility model, brownian perturbations of the UAVs are also analyzed. An out-of-band (OOB) reduction of 320 dB with an improved BER of 1 × 10 − 4 in 16-QAM for a signal-to-noise ratio, E b / N 0 , of 20 dB is achieved.


Introduction
An uncrewed aircraft handled by remote control or embedded computer programs is commonly known as an unmanned aerial vehicle (UAV) or drone. Recently, there has been tremendous amount of interest growing to develop UAV-ground communications using low-cost massive UAVs under existing 5G as well as future-generation (B5G/6G) cellular networks. To maintain secure and reliable flight operation, UAVs can exchange safetycritical information with remote pilots, closest aerial vehicles, and air traffic controllers with an assistance of control and non-payload communication (CNPC). UAVs are also delivering goods and improving the throughput of 5G networks [1][2][3].
Additionally, mission-oriented data, such as high-resolution video, data packets, and aerial images, are possible to seamlessly transfer with UAVs in payload communication. Even in an emergency situation when terrestrial mobile stations go through unexpected discontinuity due to disasters, UAVs can be utilized as a base station (BS) or as a relay the complexity. In [14], a orthogonal frequency division multiplexing (OFDM) combined index modulation signaling scheme is incorporated with UAV communication systems and results in higher spectral efficiency, higher energy efficiency, and lower BER. Generally, the OFDM system is associated with drawbacks related to issues such as high peak-to-average power ratios (PAPR), bandwidth loss associated with the cyclic prefix (CP), and high outof-band (OOB) emissions. To meet the challenges of fifth generation (5G) and B5G wireless communication networks to handle peak data rates of at least 1 Tb/s, over-the-air latency of 10-100 µs, and user-experienced data rates of 1 Gb/s, generalized frequency division multiplexing (GFDM), non-orthogonal multicarrier modulation (MCM) scheme may be preferred in lieu of OFDM. The future generation wireless networks should target ensuring low power consumption for machine-to-machine communication and low latency for tactile internet and internet-of-things (IoT) [15]. By reducing inter-carrier interference (ICI) and inter-symbol interference (ISI) through a properly designed better than raised cosine (BTRC) pulse shaping filter, improved sinc power (ISP), parametric linear pulse (PLP), and parametric exponential (PEXP), bit error rate (BER) performance is improved in the GFDM system [16]. In [17], a complementary GFDM repeated transmission scheme is proposed to enhance BER performance for mMTC and IoT. A system model is presented by the authors of [18] to develop an Einstein product of a tensor-based precoding scheme suitable for MIMO GFDM systems. In [19], a joint data symbol detection and phase noise compensation algorithm is proposed to reduce the effects of phase noise in GFDM systems. In this present study, a cooperative UAV-aided millimeter wave (mmWave) CP-free GFDM system is proposed considering data transmission for ground mMTC. In broadband wireless communication systems, one of the most popular modulation schemes, CP-OFDM, utilizes the CP to make the frequency-selective multipath fading channel appear as flat-fading [20]. As the CP causes overhead in the spectrum efficiency and power efficiency, the use of a CP is being avoided, and guard intervals (GIs) are being introduced between two consecutive GFDM symbols. In this UAV relay-assisted terrestrial networking system, more focus is attributed to enhancing secure transmission using Lorenz's hyperchaos mapping system addressed in [10].
This paper mainly deals with a multi-antenna-configured UAV-aided cooperative CPless GFDM system for secured B5G wireless networks. The physical size of each antenna of the multi-antenna-mounted UAV is very small since mmWave transmission frequency (28 GHz) is considered for this scheme. Generally, MIMO antenna-configured systems are preferably used to enhance the diversity gain and data rate of the system concerned. As the total transmission power of the UAV is limited, a single power-amplifying section for each transmitting layer is utilized here instead of using large RF chains, and thus, any negative effect on operational life cycle of the UAV is being minimized.
NOMA, GFDM, universal filtered multicarrier (UFMC), and filter bank multicarrier (FBMC) has always been regarded as an excellent contender waveform for 5G networks. Among the mentioned, GFDM is a symbol-based multicarrier multiplexing technique which keeps the signal well confined in time and frequency domains. It utilizes circular filter for effective reduction in the low out-of-band (OOB) emission and ensures simultaneous transmission of multiple symbols at different time slots. By taking all these properties into account, a GFDM-aided system is presented in place of the OFDM-based signaling technique. As far as the novelty of this work is concerned, designing a system for the mMTC network utilizing a CP-less UAV-aided GFDM system has not been well investigated before. Significant contributions of the paper are outlined as follows: • T-transformation spreading codes in combination with the walsh Hadamard transform (WHT) and discrete Fourier transform (DFT) is proposed along with recently the introduced gyre precoding (GP) [21] technique. This reduces multi-user interference (MUI) and computational complexity. • The BTRC pulse shaping filter and concatenated channel-coding techniques are utilized to reduce OOB power emission and improve BER performance.
• Shen et al.'s proposed a 3D controlled Lorenz mapping system [10] is implemented as a cryptographic technique in this proposed system for PLS encryption. • The 3D mobility model is also introduced to observe 3D positional uncertainty of the UAVs due to Brownian motion and analyzed its impact on the achievable ergodic rate. • Simulation results verify the significant enhancement of the achievable ergodic rate, bit error rate (BER) performance, OOB reduction at a reasonably acceptable PAPR.
The structure of this paper is as follows. Section 2 describes the system model, including the block diagram and insights of the signal model considered. In Section 3, simulation parameters, numerical results, and outcomes of the proposed system are presented. Finally, a brief discussion on final remarks and future directions are provided in Section 4.

Block Diagram
The conceptual block diagram of the secured cooperative UAV-based CP-less GFDM system is shown in Figure 2. In its ground segment section, initial significant signal processing operations are executed on the ground BS. Ground BS is set to be comprised of N T (=16) transmitting antennas, which should serve eight ground mMTC users via one main UAV and two cooperative UAVs, and all the UAVs are N T × N T multi-antenna configured. The synthetically generated binary data assigned for individual mMTC users are processed with concatenated channel-coding schemes based on (3,2) single parity check (SPC), low density parity check (LDPC), and repeat and accumulate (RA) [22,23] and subsequently encrypted with improved 3D controlled Lorenz mapping-based encryption technique [10]. Initially, complex symbols are produced via the digital modulation of the encrypted binary data [24]. After that, the generated complex symbols undergo T-transformation with a combined effect of the WHT and DFT [25] and zero padded prior to the execution of the GFDM modulation. In the GFDM modulation, the zero padded complex symbols are circularly convolved with the BTRC filter, followed by frequency up conversion [16,26]. The GFDM modulated symbols are power scaled and subsequently precoded with channeldependent multi-user transmit GP schemes [21]. The transmission of the spatially multiplexed precoded signals from the ground BS to the main UAV with ending signal processing operations in the ground segment section of this proposed system.
In the UAV segment, the main UAV decodes the transmitted signal with the ZF linear signal detection technique [27], and the power of the detected signal is boosted up to make it compatible with the power of the ground BS. In the first time phase, the main UAV transmits the detected signal simultaneously to the cooperative UAVs working in half-duplex mode and the eight ground mMTC users, each of which is equipped with N R (=2) receiving antennas. In the second time phase, the cooperative UAVs simply adopt the amplifyand-forward (AF) relaying strategy [28] and forwards the signal to the ground mMTC users with identical transmission power of the main UAV. On the basis of the estimated effective channel between the UAVs and the individual mMTC user in two consecutive time phases and their respective beamforming weight, mMTC user's own transmitted signal is detected at the receiving end of each mMTC user [29]. The detected signals are power scaled again to be restored at the original power. The GFDM demodulation process is then initiated by the removal of the padded zeros, and the multiplication with inverse T-transformation process yields despreaded signals. Despreaded signals are then passed through the phases of digital demodulation, channel decoding, and improved 3D controlled Lorenz mapping-based decryption to eventually recover the transmitted data.

Signal Model
In the previous section, various signal processing techniques have been addressed by presenting their applicability in this proposed system. In this section, a significant emphasis has been put on the signal model to describe various useful processing techniques. This can be divided in two segments first, which are: (1) the Ground Segment and (2) the UAV Segment.

Ground Segment
To improve both transmission and the performance in terms of security, the 3D controlled Lorenz mapping system addressed in [10] can be represented as: where the initial values of and x, y, z variables and key parameter values (x, y, z, ρ, ζ, r, µ, β) have been assigned to (1.99, 2.2, 2.0, 9, 2, 35, 25, 8/3). The solution component {z(t)} 768 t=1 can be used to produce 768 × 1 matrix-sized primary key K 0 as: where · is indicative of rounding operation to the nearest integer. The primary key K 0 is disseminated into eight encrypted keys with each key containing 96 elements. The eight keys for the eight mMTC users can be considered as: K 1 = 0.5K 0 (1 : 96, 1) K 2 = 1.0K 0 (97 : 192, 1) K 3 = 1.5K 0 (193 : 288, 1) K 4 = 2.0K 0 (289 : 384, 1) K 5 = 2.5K 0 (385 : 480, 1) K 6 = 3.0K 0 (481 : 576, 1) K 7 = 3.5K 0 (577 : 672, 1) where the first to 96th elements of K 0 is designated by K 0 (1 : 96, 1), 97th to 192th elements of K 0 is designated by K 0 (97 : 192, 1), and so on. The elemental values of all the keys K 1 through K 8 fall within the range, and each elemental value provides eight binary bits. In binary form, each of the generated keys K 1 through K 8 for eight mMTC users is of 1 × 768 matrix in size. By using the repmat function in MATLAB, the length of the binary data of each of the generated keys K 1 through K 8 is made identical to the data length of the concatenated channel encoded data. The encrypted concatenated channel encoded binary data vector b¯k of data lengthN for mMTC userk can be represented as: where the symbol ⊕ is indicative of the XOR operation, andK¯k is the ultimately generated key for the mMTC userk. With the data lengthÑ, the conversion of binary data vectorb¯k as a complex symbol vector yields c¯k through digital modulation. The T-transformed 8 × 8 sized matrix based on combining the WHT and DFT can be represented as: where ω = e j2π/8 . The signal model presented in Equation (5) is applicable row wise for different mMTC usersk. The first row is applicable to mMTC user 1, the second row is applicable to mMTC user 2, and so on. The application of the T-transformation spreading technique on the modulated complex symbol vector ck for mMTC userk yields the complex symbol vectorck of data length 8Ñ.
In GFDM structured block generation, the data vectorck is rearranged into a K × M matrixsized GFDM block, where M represents the number of subsymbols and K denotes the number of subcarriers. For every GFDM block with N(=KM) samples, its first and last subsymbol contain null subcarriers to avoid collision in between two consecutive GFDM blocks, eliminate the effect of the fading channel, and additionally provide significant OOB emissions reduction for neighboring channel interference reduction. Considering data symbol d k,m,l in each GFDM block for mMTC user k, the lth GFDM block for mMTC userk can be seen as a result of circular convolution with BTRC, g[n], followed by frequency up conversion and can be written as [26]: where denotes circular convolution,k = 1, 2, . . . , 8 is the mMTC user identification number, l = 1, 2, . . . , L is the GFDM block identification number, δ[n − mK] is the Dirac function with n = 0, 1, 2, . . . , N − 1, and m = 0, 1, 2, . . . , M − 1.
The BTRC pulse shaping filter addressed in [21] and used in signal model of Equation (5) can be written in the time domain as: where β = ln2 αB , T = M∆T = N F s is the GFDM symbol period, α is the filter roll-off factor of the BTRC pulse shaping filter, ∆T = K F s is the period of each subsymbol, F s is the sampling frequency, ∆ f = 1 ∆T represents the subcarrier spacing, and B = K∆ f = K ∆T represents the GFDM signal bandwidth. The lth discrete time domain GFDM signal (Equation (6)) for mMTC userk is rescaled such that E|xk ,l | 2 = 1, and the power rescaled signal vector of unity power is represented byxk ,l [n]. Concatenation of all the signal vectors ofxk ,l [n] yields the desired signals for the mMTC userk and are denoted by a data matrix Dk of size N × L. All the elements of Dk are stacked within a single column vector and further arranged into a matrixDk of N R × N L (= N × L)/N R in size, and theDk can be represented as:Dk From Equation (8), all the data for eight mMTC users can be confined into a single data matrix D of size N T × N L .
Prior to the application of proper transmit beamforming weights to the matrixed data, D, it is desirable to design transmit precoders utilizing the 3D geometrical MIMO flat-fading channels, which are based on the probabilistic path loss model, and in such estimated fading channels, the 3D mobility model has been incorporated to consider 3D positional uncertainty of the main UAV and cooperative UAV's positions relative to their fixed positions due to Brownian motion. The 3D mobility model can be used to present positioning errors using homogeneous stochastic differential equations in the 3D spatial (x, y, z) coordinate system as: ; dw 1t , dw 2t , dw 3t are the differential form of the Brownian motion in three mutually perpendicular x, y, and z directions and are equally applicable for all UAVs. In Brownian motion, it can be assumed The (x t,s , y t,s , z t,s ) are the 3D positional coordinates of the s-th UAV (s = 0 for main UAV, s = 1 for cooperative UAV #1 and s = 2 for cooperative UAV #2). The (x t,s ,ỹ t,szt,s ) are the 3D errored positional coordinates of the s-th UAV due to Brownian motion [30,31]. With the consideration of Brownian motion, the x, y, and z coordinates of the main UAV and the cooperative first and second UAVs are (x t,0 ,ỹ t,0zt,0 ), (x t,1 ,ỹ t,1zt,1 ), and (x t,2 ,ỹ t,2zt,2 ), respectively. The 3D distance d k,j between main/cooperative UAVj and the ground mMTC userk with its 3D coordinates (x¯k, yk, zk) can be represented as: the elevation angle θ k,j can be represented as: . The LOS probability from UAVs to ground mMTC communications can be represented as [32]: where s = 0, 1, 2, and ψ (=11.95) and β 0 (=0.14) are the functions of the carrier frequency and environment. The path loss between the main/cooperative UAVj and the ground mMTC userk can be represented as: where α is the path loss exponent, f c (=28 GHz) is the carrier frequency, η 1 (=10 3/10 ) represents the excessive path loss coefficient for LOS link, η 2 (=10 23/10 ) denotes the excessive path loss coefficient for NLOS link, and c (=3 × 10 8 ) m/s represents the speed of light. The simple form of NLOS probability LOS . Between the main/cooperative UAVsj and the ground mMTC userk, the average path loss equation can be represented as: With its 3D coordinates [x 0 , y o , z o ] to the main UAV, the path loss L 0 of the MIMO flat-fading channel H 0 for the case of uplink transmission from ground BS can be represented as: where .
Considering very high LOS probability and a low excessive path loss for transmission between the main UAV and the cooperative UAVs, the path losses L 01 and L 02 for the MIMO flat-fading channels H 01 and H 01 existed between the main UAV and first UAV, and the main UAV and second UAV can be represented as: where d 1,m and d 2,m are the corresponding 3D distances between the main UAV and cooperative UAVs. Considering the applicability of the signal models represented in Equations (11) where (·) T indicates the non-conjugate transpose operator. Considering the GP, a novel linear multiuser multiple-input multiple-output (MU-MIMO) precoding scheme cited at [21], the ZF precoding of matrix size N T × N T can be represented as: where (·) H indicates the conjugate transpose operator.

UAV Segment
Due to precoding, signal power varies from rowed to rowed data of X 0 , and to make identical signal power in every rows of X 0 , it is to be multiplied by a normalized signal power P 1 such that the transmitted power from the ground BS would be P gr and the arrived signal at the main UAV can be represented as: where n m ∼ Cℵ(0 N T , σ 2 m I N T ) denotes the additive white Gaussian noise (AWGN). By introducing the ZF signal detection scheme cited at [27], the decoded signal can be represented as:X The signal model presented in Equation (22) is rescaled to make it compatible to a desirable UAV transmission power P uav of matrix size N T × N T by multiplying with a normalized signal power P 2 , and the power scaled signal can be represented as: In the first time phase, the main UAV transmits its signalX 0 to all cooperative UAVs and ground mMTC users. The signal received at the ground mMTC userk can be represented as: where n 1 k ∼ Cℵ(0 N R , σ 2 k I N R ) denotes the additive white Gaussian noise (AWGN). The signal received at the two cooperative UAVs are given by: where n 1 01 ∼ Cℵ(0 N T , σ 2 01 I N T ) and n 1 02 ∼ Cℵ(0 N T , σ 2 02 I N T ) denote the AWGN, and the superscript 1 on different variables indicates the first time phase. The amplifying gain in this system at the two cooperative UAVs can be represented as: In the second time phase, both the cooperative UAVs transmit its amplified signals to the ground mMTC users. The signal received at the ground mMTC userk from cooperative UAVs can be represented as: where n 2 k1g ∼ Cℵ(0 N R , σ 2 k1g I N R ) and n 2 k2g ∼ Cℵ(0 N R , σ 2 k2g I N R ) are the effective AWGN in transmission from the first and second cooperative UAV, respectively.
In the second phase, the total signal received by the ground mMTC userk can be represented as: where the superscript 2 on different variables indicates the second time phase. By combining the signals received at the ground mMTC userk in both phases [33], the following can be written: where Hk = Hk (H 1k G1H 01 + H 2k G2H 02 ) is the effective channel between the UAVs and ground mMTC userk and is of 2N R × N T matrix in size, while the combined signal Yk d and AWGN noise nk ∼ Cℵ(0 2N R , σ 2 k I 2N R ) are of 2N R × N L matrix in size. On further simplification, Equation (29) can be re-written as: where Hk = Hk √ P 1 P 2 is the effective channel for the signal model presented in Equation (30) and it is also of 2N R × N T matrix in size. Generally, such non-symmetric complex channel matrix is found to be sparse in nature and due to sparsity of matrix, its pseudo inverse operation for achieving exact solution from signal model presented in Equation (30) cannot be obtained. Using the regularized ZF equalization technique [34] to invert the effect of the channel Hk, the detected signal in receiverSk d of N T × N L matrix in size can be represented as: where H = 1.0 × 10 −19 is the regularization parameter, andnk is the AWGN noise. From Equation (31), it is quite obvious that the received signal in ground mMTC userk contains its own and interference signals from other mMTC users. Multiplying Equation (31) by Hk, a modified form of the received signalSk d of N R × N L matrix in size is obtained, and it can be represented as: In Equation (32), all the components are of N R × N L matrix in size, and the noise component represents the AWGN noise. Again, ||H¯kW¯k|| 2 represents the instantaneous total signal power in two consecutive time phases for the mMTC userk, the total instantaneous interference power of the signal received by mMTC user k for thej-th mMTC user is represented by ||H¯kW¯k|| 2 , and N R σ 2 k denotes the AWGN noise power. The received signal-to-interference-plus-noise ratio (SINR) for the mMTC userk(SI NRk) after gyre decoding can be represented as: The logarithmic function of SINR represents the achievable ergodic rate Rk for mMTC userk, and it can be written as after gyre decoding [35]: However, from the signal model presented in Equation (32), it is quite obvious that the reduction in the interfering signal power is very much significant as the GP removes MUI by designing GP beamforming weight Wk assigned for mMTC userk to fall in the null space of the total MIMO flat-fading channel matrix H and producesHkWj = 0. Considering HkWk =H¯k and multiplying Equation (32) (H H kHk ) −1HH k , a new form of the processed signal for mMTC userk is obtained as: where the first, second, and third terms of Equation (35) are indicative of signal, interference, and noise components n = (H H kHk ) −1HH k nk and is ∼ Cℵ(0 N R , σk I N R ), respectively. With further processing for necessary steps with the execution of the inverse T-transformation, achievable ergodic rate, SINR, and SNR values can be estimated.

Numerical Results and Relative Discussion
This section presents the numerical outcomes for the secured mmWave cooperative UAV-aided CP-less GFDM system. The considered simulation parameters are presented in Table 1. Based on the noise power of −80 dBm, which corresponds to 1 × 10 −11 W, the Additive white Gaussian noise (AWGN) component is introduced at different receive antennas of the main UAV [36]. In this scenario, there is a ground BS, one main and two cooperative UAVs, which are connected in an integrated UAV-terrestrial network to serve eight massive machine-type communications (mMTC) users, which are 1. Considering the heights of the buildings and ground base station, the cooperative UAV's trajectory path is preferably being considered as circular and not random. Thus, the chance of any conflict with existing infrastructure is absolutely minimal. Three-dimensional Brownian perturbation model-based [30] trajectory planning combines the randomness with directed motion, and thus, the role of UAVs in successful signal transmission is being observed. The energy consumption of each UAV associates both communication-related energy consumption and propulsion energy consumption. Here, a single UAV controlling ground node (GN) based energy-efficient rotary-wing UAV communications is considered, and a 3D-based Brownian perturbation mobility model is applied, which controls UAV speed/deviation from the UAV target position, which, in turn, reduces propulsion energy consumption [37]. Basically, the proposed scheme can importantly balance the additional consumption of power of UAVs. The transmitting UAVs and ground mMTC users have different channel state information (CSI). In such case, MIMO flat-fading channels are estimated utilizing the probabilistic path loss model cited in [32].
In this study, the implementation of the PLS technique is based on the utilization of Lorenz's hyperchaos mapping system presented in Figure 3. Chaotic systems have been extensively used in physical layer encryption because of their high initial sensitivity and good randomness. In chaotic encryption schemes, sharing of the initial value between the transmitter and receiver is an important factor of system security, which makes it difficult for illegal users to steal information through malicious attacks. To achieve higher security performance, Lorenz's hyperchaos mapping system is considered by introducing new state variables in the encryption. Figure 3 shows that the behavior of the system is fully dependent on the considered parameter values. In Figure 4, a 3D plot of various keys assigned to individual mMTC user is shown, which confirms that assigned keys are different for different mMTC users.
In Figure 5, the cumulative distribution function (CDF) of the 3D positioning errors for both the main and two cooperative UAVs have been presented in a static scene. It can be seen from the CDF that the error at 95% confidence is 0.08 m in vertical Z direction and average 0.06 m in both horizontal X and Y directions for the main UAV. In the case of the cooperative UAV #1 with identical confidence level, the positional error is 0.06 m in vertical Z direction and 0.12 m in horizontal X direction and 0.05 m in horizontal Y direction. In the case of cooperative UAV #2, the positional error is 0.06 m in vertical Z direction and 0.02 m in horizontal X direction and 0.04 m in horizontal Y direction. In all cases, the 3D positional errors of the UAVs are reasonably acceptable in the context of state-of-the-art.   The impact of both gyre decoding and the T-transformation scheme on the received signal-tonoise ratio (SNR), received SINR, and achievable ergodic rate with the variation in UAV transmission power for eight mMTC users have been illustrated through Figures 6-8. In Figure 6a, the estimated average received SNR values for eight mMTC users without the application of gyre decoding and inverse T-transformation schemes are presented. From Figure 6a, it is seen by analyzing all the estimated SNR values for mMTC user 1 through mMTC user 8, the maximum average SNR value is approximately 30.533 dB, and the minimum average SNR value is approximately 20.529 dB. In Figure 6b, the estimated average SNR values for eight mMTC users with the implementation of both gyre decoding and inverse T-transformation schemes at the receiver are presented. There is no significant difference observed while comparing with the SNR values in the absence of gyre decoding and T-transformation schemes, as shown in Figure 6a. It is seen from Figure 6b by estimating all SNR values for mMTC user 1 through mMTC user 8 that the maximum average SNR value is approximately 29.663 dB, and the minimum average SNR value is approximately 19.659 dB.  Figure 7b shows that due to the simultaneous implementation of both gyre decoding and inverse T-transformation schemes based orthogonal spreading codes, the interfering signal power is totally nullified, causing improved SINR performance as compared to the results presented in Figure 7a. The estimated SINR values for mMTC user 1 through mMTC user 8 are identical to SNR values presented in Figure 6b. In Figure 8a, it is observable that even when the UAV transmission power is varied, the estimated achievable ergodic rate for different mMTC users is comparatively low and almost constant. The estimated achievable ergodic rate for mMTC user 1 through mMTC user 8 are 1.2345 bps/Hz, 0.5120 bps/Hz, 0.7976 bps/Hz, 0.8790 bps/Hz, 0.3894 bps/Hz, 0.4243 bps/Hz, 0.2263 bps/Hz, and 0.2731 bps/Hz, respectively. In Figure 8b, it is observable that the estimated achievable ergodic rates for different mMTC users are very much improved with the application of gyre decoding and T-transformation scheme than in the absence of them, as shown in Figure 8a. It is also noticeable from the figure that the achievable ergodic rate values decrease when the UAV transmission power is varied. By analyzing all the estimated achievable ergodic rates for mMTC user 1 through mMTC user 8, the maximum and minimum average achievable ergodic rate is found to have values 9.856 bps/Hz and 6.553 bps/Hz, respectively. Considering Brownian motion, it is observable that for the x, y, z coordinate variation of 0.2 m relative to UAV's fixed positions, the mean error in estimated achievable ergodic rates at different E b /N 0 values is found to be of 0.26 bps/Hz. On the other hand, in the case of x, y, z coordinate variation of −0.2 m relative to UAV's fixed positions, the mean error in estimated achievable ergodic rates is of −0.83 bps/Hz. In Figure 9, the estimated achievable ergodic rate values for ground mMTC users confined within a circular region of a 50 m radius are presented. In such case, ground BS transmit power (46 dBm = 39.81 W) and UAV base station transmit power (30 dBm = 1 W) from Table 10.3 of [5] has been considered. From the 3D displayed achievable ergodic rate of Figure 9 with identical heights of all the UAVs at 100 m, 200 m, and 300 m, it is seen that with increase in UAV height, the estimated achievable ergodic rate values decrease. The maximum estimated achievable ergodic rate in the case of UAV heights of 100 m, 200 m, and 300 m are found to have values of 14.26 bps/Hz, 13.19 bps/Hz, and 12.21 bps/Hz.
In Figure 10, estimated maximum achievable ergodic rate values at different normalized signalto-noise ratios (E b /N 0 ) are compared with the similar works of different authors from [38][39][40]. It is quite clear from the figure that the proposed system outperforms all the other similar works utilizing different systems. Performance in terms of achievable ergodic rate is even better with gyre decoding and T-transformation scheme utilized in tandem than merely with gyre decoding.
An acceptable out-of-band ( Figure 11 is a typical case of OOB power reduction performance of mMTC user 1.   To illustrate the efficacy of the proposed system further, BER performance is also analyzed by implementing concatenated channel coding with multi-user beamforming weighting-aided maximumlikelihood and zero forcing (ZF) detection technique adopting 16-QAM digital modulation. It is quite easily remarkable from the simulation outcomes presented in Figures 12-14 that the BER performance of the proposed system implementing (3,2) single parity check (SPC) and repeat and accumulate (RA) concatenated channel-coding scheme is better than other two concatenated channel-coding schemes tested in this system.   It is seen from Figure 12 that the BER varies from 33.30% to 29.40% for a typically assumed E b /N 0 value of 10 dB with the system utilizing low density parity check (LDPC) and (3,2) SPC concatenated channel-coding technique. Overall, BER performance is not quite satisfactory for any of the mMTC users. In SNR region of 11 dB to 21 dB, some mMTC users are getting their signals with almost identical BER.
In Figure 13, comparatively improved system performance is noticeable with the system tested under the LDPC and RA concatenated channel-coding technique. It is seen by observing simulated data of all eight users in Figure 13 that the BER varies from approximately 24.39% to 12.00% for a typically granted E b /N 0 value of 10 dB.
In Figure 14, BER performance of the system under concatenated channel-coding scheme with a combination of (3,2) SPC and RA is quite satisfactory. Over the wide range of SNR from 0 dB to 30 dB, the estimated maximum and minimum BER are found to have values 0.14190 (14.19%) and 0 (0%), respectively. It is seen from Figure 14 that the BER varies from 2.75% to 0.51% for a typically assumed E b /N 0 value of 10 dB. It is also observed that in the case of mMTC user 1, BER is reduced from 5.90% to 0.68% with an increase in the SNR value from 6 dB to 16 dB.
In Figure 15, BER performance of the proposed system utilizing (3,2) SPC and RA concatenated channel-coding technique has been compared with other GFDM systems introduced in previous works. From a brief description perspective of the systems presented in [39,[41][42][43], it can be said that the authors in [41] did a thorough analysis on the applicability of wireless energy harvested GFDM-based cooperative network and presented numerical results in terms of BER. In [42], the effectiveness of localized discrete gabor transform (LDGT) algorithm implemented the GFDM system with varying windowing lengths was studied accounting the BER performance. Authors in [39] conducted a performance evaluative study for a GFDM-assisted NOMA system with respect to BER and achievable rate analysis. In [43], the authors highlighted the effectiveness of a designed pulse shaping filter based on quadratic programming in reducing the OOB radiation and BER performance of the GFDM system. Now, it is quite clearly visible from Figure 15 that the proposed system is showing comparatively better performance in terms of BER against E b /N 0 than cooperative GFDM, LDGT-aided GFDM, GFDM-based NOMA, and designed pulse shaping filter implemented GFDM systems. Figure 15. Comparative BER performance analysis between the proposed system and other B5G compatible multicarrier systems.
In Figure 16, it is clear from the high BER values that the improved 3D controlled Lorenz mapping system-aided PLS encryption technique is very much effective in offering secure multi-user data transmission. In addition, it is difficult to identify the transmitted data correctly at the receiver without any knowledge of assigned parameter values used in such an encryption technique.
Due to the utilization of a channel-dependent precoding technique in this proposed system, the power of different signals at different transmitting channels is varied. As a consequence, PAPR values also vary. It is observable from Figure 17 that at a complementary CDF (CCDF) of 10 −3 , the maximum and minimum PAPR values are 9.5 dB (Transmitting antenna 2) and 8 dB (Transmitting antenna 4), respectively. As it is seen that the variation among all the transmitting channels is well within the PAPR value of 1.5 dB, curves for different channels coincides at some stage.
In Figure 18, it is observable that the estimated PAPR values at different ground transmitting channels in the proposed system are reasonably acceptable if they are considered with respect to the state-of-the-art.

Conclusions
A UAV-aided mmWave downlink cooperative CP-less GFDM system is proposed in this paper. UAVs are combined with a terrestrial cellular network, and emphasis is placed on physical layer security for massive machine-type communication (mMTC) users in such a network. In this study, the implementation of an efficient low-complexity T-transformation spreading code as well as MU-MIMO Gyre precoding resulted in the reduction in both MUI and the enhancement of the achievable ergodic rate. Utilizing null subcarriers at both ends of a GFDM block for sub-time symbols and better than raised cosine (BTRC) filter resulted in significant reduction in out-of-band (OOB) spectrum power by eliminating the effects of the fading channel at mmWave. The 3D mobility model resulted in reducing the effects of the Brownian motion of UAVs. In this proposed system, the performance in terms of BER is not quite satisfactory. However, the BER performance is reasonably acceptable with the combined implementation of (3,2) SPC and RA channel-coding schemes and outperforms other proposals found in the literature. In the future, introducing massive MIMO to the proposed system can improved it and make it more robust.