#### Abstract

With the further research in communication systems, especially in wireless communication systems, a statistical model called Nakagami-*m* distribution appears to have better performance than other distributions, including Rice and Rayleigh, in explaining received faded envelopes. Therefore, the Nakagami-*m* quantile function plays an important role in numerical calculations and theoretical analyses for wireless communication systems. However, it is quite difficult to operate numerical calculations and theoretical analyses because Nakagami-*m* quantile function has no exact closed-form expression. In order to obtain the closed-form expression that is able to fit the curve of Nakagami-*m* quantile function as well as possible, we adopt the method of curve fitting in this paper. An efficient expression for approximating the Nakagami-*m* quantile function is proposed first and then a novel heuristic optimization algorithm—generalized opposition-based quantum salp swarm algorithm (GO-QSSA)—which contains quantum computation, intelligence inspired by salp swarm and generalized opposition-based learning strategy in quantum space, to compute the coefficients of the proposed expression. Meanwhile, we compare GO-QSSA with three swarm intelligence algorithms: artificial bee colony algorithm (ABC), particle swarm optimization algorithm (PSO), and salp swarm algorithm (SSA). The comparing simulation results reveal that GO-QSSA owns faster convergence speed than PSO, ABC, and SSA. Moreover, GO-QSSA is capable of computing more accurately than traditional algorithms. In addition, the simulation results show that compared with existing curve-fitting-based methods, the proposed expression decreases the fitting error by roughly one order of magnitude in most cases and even higher in some cases. Our approximation is proved to be simple and efficient.

#### 1. Introduction

It is well known that radiowave propagation is complicated under the environment of wireless channels and various effects exert different impacts on the characteristics of radiowave propagation. Considerable research focusing on building the statistical models and characterizing these effects has been undertaken. Multipath fading is of great significance among the effects and is caused by the combination of various signal components containing randomly diffracted, reflected, delayed, and scattered components. There are different models used under different radio propagation environments, such as Rayleigh, Rice, Hoyt, and Weibull distributions [1], and these four models are capable of describing the corresponding statistical characterization of the envelope caused by multipath fading. Propagation paths can be modelled by Rayleigh distribution if line-of-sight (LOS) component does not exist [2]. Propagation paths containing many random weaker components and one strong direct LOS component are usually modelled by the Rice distribution, such as suburban land-mobile and microcellular urban components with the first resolvable LOS paths [3]. Among these models, the best statistical model that is capable of defining different multipath fading conditions is the Nakagami-*m* distribution. Nakagami-*m* distribution is diffusely adopted in wireless communication systems due to its two main advantages. One is its flexibility feature; in fact, Rayleigh distribution is essentially the special case when the fading parameter of Nakagami-*m* distribution is set as 1 [1]. Furthermore, the channel modelled with Nakagami-*m* distribution can converge to a nonfading AWGN channel, and different parameters even make it a close approximation to the Hoyt and Rice distribution. Another significant feature is that the Nakagami-*m* distribution is a perfect statistical model to fit scintillating ionospheric radio links, indoor-mobile multipath propagation, and land-mobile multipath propagation [4–8].

Numerical calculations and theoretical analyses in wireless communications under the Nakagami-*m* fading environment can be operated with inverse and normal CDFs. However, it is difficult to utilize Nakagami-*m* quantile function firsthand to operate these calculations because the exact closed-form expression of Nakagami-*m* quantile function does not exist. In recent years, a new method to approximate the quantile function of Nakagami-*m* distribution by dint of the curve-fitting method has been illustrated, and scholars have utilized several metaheuristic optimization algorithms to compute the coefficients of their expressions correspondingly, such as the genetic algorithm [4], backtracking search optimization algorithm [5], and artificial bee colony optimization [6].

In this paper, an expression with lower computational complexity and better performance for approximating the Nakagami-*m* quantile function is proposed firstly, and then we propose a novel algorithm called generalized opposition-based quantum salp swarm algorithm (GO-QSSA) to compute coefficients of the proposed expression. Here, we propose the generalized opposition-based learning strategy in quantum space, and we apply it to QSSA. Finally, our simulation results under different parameters of Nakagami-*m* quantile function prove the accuracy of the proposed expression and the efficiency of GO-QSSA.

The contributions of our work are summarized as follows:(1)We propose a briefer and more efficient expression to approximate the quantile function of the Nakagami-*m* model. Compared with the approaches in [9], our curve-fitting-based method is superior both in the computational complexity and the number of expression coefficients. Among all the existing curve-fitting-based methods in [4], the fitting equation is characterized by a division comprised of two polynomial equations and five parameters in total. However, the expression proposed in this paper is just a polynomial containing three simple elementary functions and four parameters which we need to compute. Compared with References [4–6], our simulation results prove the strengths of our approximation.(2)We propose the generalized opposition-based learning strategy in quantum space, which is first proposed in the literature. Meanwhile, we apply this strategy to QSSA to improve its convergence speed.(3)We propose the quantum intelligent algorithm which is named GO-QSSA. It is the first time that quantum swarm intelligence is combined with salp swarm algorithm and quantum evolutionary strategy, and the results prove the accuracy of the proposed expression and the efficiency of GO-QSSA.

The remaining part is structured as follows: Section 2 is the related work. In Section 3, the Nakagami-*m* model and approximation problem model are introduced. In Section 4, the methodology is illustrated, the expression proposed to approximate Nakagami-*m* quantile function is presented, we propose generalized opposition-based quantum salp swarm algorithm, and then we use GO-QSSA to solve this problem. The enormous simulation results and numerical analyses are shown in Section 5. In Section 6, finally, the main work of this paper is summarized.

#### 2. Related Work

Many researchers have focused on numerical calculations and theoretical analyses under the Nakagami-*m* fading environment with its inverse and normal CDFs. Das et al. in [8] considered the mobility of nodes and derived the CDF and PDF of the received power at the mobile nodes in the Nakagami-*m* fading environment. Eventually, they deduced the expressions for average bit error rate (BER) and coverage probability. In [10], a finite network of drones (also called unmanned aerial vehicles, UAVs) serving a given region was considered, and the downlink coverage probability was derived assuming all wireless links were compliant with Nakagami-*m* distribution. Hou et al. [11] studied the outage performance considering Nakagami-*m* fading with different fading parameters and the fixed power allocation of nonorthogonal multiple access (NOMA). In addition, a communication network assisted by a wireless powerbeacon (PB) was investigated in [12]. In this network, one receiver, multiple transmitters, and multiple PBs were considered. These transmitters harvest energy from the PBs and transmit their data, respectively, to the receiver. Fading coefficients were assumed to be Nakagami-*m* random variables during the transmission. The authors derived the expression of achievable throughput and an approximated closed-form expression of throughput when the signal-to-noise ratio was high. Moreover, in [13], under the Nakagami-*m* and Hoyt fading channels, analytical expressions of the symbol error rate (SER) system were deduced, and the simulation results confirmed the GFDM/OQAM model proposed by the authors. These illustrations above enlighten us to the fact that the Nakagami-*m* quantile function is vital to calculate some performance indicators of wireless communication systems such as BER, throughput, coverage probability, and outage probability. However, it is difficult to utilize Nakagami-*m* quantile function firsthand to operate these calculations because the exact closed-form expression for Nakagami-*m* quantile function does not exist. Some methods have been proposed to acquire the approximated Nakagami-*m* quantile function such as Gaussian distribution and Hastings’ approach [9], whereas these approaches have some defects including their high computational complexity and the complex form of their expressions.

In recent years, a new method to approximate the quantile function of Nakagami-*m* distribution by dint of the curve-fitting method has been illustrated, and scholars have utilized several metaheuristic optimization algorithms to compute the coefficients of their expressions correspondingly, such as the genetic algorithm [4], backtracking search optimization algorithm [5], and artificial bee colony optimization [6]. Bilim and Develi[4] proposed an expression consisting of a cubic polynomial and a linear polynomial and then used genetic algorithm to compute its total five coefficients. Kabalci [5, 6] investigated a more concise expression with four coefficients. Meanwhile, he adopted backtracking search optimization algorithm and artificial bee colony optimization to solve this curve-fitting problem, respectively. In this paper, a novel and more efficient expression and a novel heuristic algorithm called GO-QSSA are proposed to approximate the Nakagami-*m* quantile function. Related theoretical background and methodology will be illustrated next.

#### 3. Theoretical Background and Problem Modelling

##### 3.1. Nakagami-*m* Model

Concepts and properties have been interpreted in [9] in detail. In this section, for clarity, we restate the concepts and main properties of Nakagami-*m* distribution. Meanwhile, the method for calculating Nakagami-*m* quantile function is presented in detail. Nakagami-*m* distribution has two parameters: a fading parameter, , denoting the fading effect of the wireless communication channel, and a scaling parameter, , controlling spread. Assuming is a random variable complying with the Nakagami-*m* distribution, its PDF is defined as follows:where is the symbol of the Gamma function. The fading parameter and scaling parameter can also be estimated aswhere and , respectively, denote variance and expectation operators.

The Nakagami-*m* distribution shows the largest range of AF (amount of fading) among all the distributions describing multipath fading characterization considered in [1]. It contains two special cases, one is the one-sided Gaussian distribution when the fading parameter is set as , and the other is Rayleigh distribution when the fading parameter is set as . Furthermore, when , it is close to the Rice distribution; when , it is capable of approximating the Hoyt distribution. Finally, the Nakagami-*m* distribution is regarded as the best statistic model to fit scintillating ionospheric radio links, indoor-mobile multipath propagation, and land-mobile multipath propagation.

According to the relationship between PDF and CDF in probability theory, the CDF is defined as

Hence, the quantile function of Nakagami-*m* can be acquired by dint of the CDF of Nakagami-*m* distribution:where is defined aswhere is the second moment of and is a transformation function. Finally, we acquire the Nakagami-*m* quantile function with the variable transformation , However, it is obvious that the consequent expression of Nakagami-*m* quantile function involves a complex integral. Moreover, it has been proven that the exact closed-form expression for the Nakagami-*m* quantile function does not exist except when [1, 9]. Therefore, we propose a novel expression to approximate the Nakagami-*m* quantile function in Section 3.2.

##### 3.2. Modelling for Approximating Nakagami-*m* Quantile Function

We have obtained the exact value of the Nakagami-*m* quantile function for a certain in the former section, and in fact, represents the Nakagami-*m* CDF value. We use to denote our proposed expression for approximating the Nakagami-*m* quantile function, which will be shown in Section 4. Assuming that the expression entirely contains coefficients denoted as , respectively, we aim at computing the optimal coefficients satisfying the best curve-fitting performance. A generally used criterion to evaluate the performance of a fitting equation is the root-mean-square error (RMSE). Here, we use the RMSE to measure the approximation error between the exact value of Nakagami-*m* quantile function and the value predicted by the proposed expression . In this paper, we use to denote the RMSE and our fitness function; then, the RMSE is calculated aswhere denotes the *i*th sampling point, is the number of sampling points, is the approximation value at , and is the exact value of Nakagami-*m* quantile function at . Here, we can model the curve fitting problem aswhere indicates the lower bound of the *j*th coefficient and indicates the upper bound of the *j*th coefficient, respectively, and represent optimal coefficients which enable the RMSE to achieve the minimum.

So far, we have modelled the approximation problem. Next, we will expound the methodology on both expression and algorithm aspects.

#### 4. Methodology

There are two main aspects to be addressed when solving this problem. One is the expression designed to approximate Nakagami-*m* quantile function, which should be not only brief but also efficient, that is to say, the proposed expression is supposed to minimize the RMSEs. Meanwhile, it should not increase excessive computational complexity; the other is the heuristic algorithm adopted to calculate coefficients of the expression. The algorithm needs to satisfy the requirements of faster convergence speed and higher accuracy. Here, we propose GO-QSSA to achieve this goal.

##### 4.1. Expression for Approximating Nakagami-*m* Quantile Function

It is well known that the foremost aspect is the mathematical expression for a curve-fitting problem, and as mentioned above, we propose a new approximation model with higher accuracy as follows:where is the inverse hyperbolic tangent function, is the exponential function, and , and are coefficients.

In order to compare clearly with available approximations previously reported in the literature, we introduce them briefly. Beaulieu and Cheng [9] followed Hastings’ approach and proposed an approximation as follows:where is an ancillary variable, is the fading parameter, , and are coefficients, and represents an approximated value of the Nakagami-*m* quantile function. Apparently, this approximation, containing two equations and five coefficients, is quite complex. In order to simplify the expression, Bilim and Develi [4] applied the curve fitting method to this problem, and they proposed the following expression:

This is much simpler than (9), although it also has five coefficients, , and . Later, Kabalci [5, 6] investigated a more concise expression as follows:

There are only four coefficients in this expression, with an inverse hyperbolic tangent function and two exponential functions. However, in practice, equation (11) does not always have better performance than (10) under various situations, which can be observed in those figures of simulation results in their papers. The defects of (11) motivate us to further investigate this problem. To some extent, the expression we propose, (8), is inspired from (11). The inverse hyperbolic tangent function is the core for both expressions (8) and (11). The deeper cause for the choice of the inverse hyperbolic tangent function is that this function is similar in shape to the Nakagami-*m* quantile function, which is illustrated in Figures 1(a) and 1(b). It can be easily observed that the Nakagami-*m* quantile function has the same trend with the inverse hyperbolic tangent function in the domain , when the fading parameter is small, and when grows, the Nakagami-*m* quantile function has the same trend with the inverse hyperbolic tangent function in the whole domain . Hence, it is convenient to adopt the inverse hyperbolic tangent function as the basis of the expression for approximating the Nakagami-*m* quantile function; meanwhile, we compensate for the deviations between the values of inverse hyperbolic tangent function and the exact values of Nakagami-*m* quantile function by dint of the product of a power function and an exponential function. Finally, we acquire the proposed expression (8). Comparing (8) with (11), we abandon the first exponential function and replace the constant in (11) with a power function. The first exponential function rarely contributes to the deviation compensation in fact, and it can be merged with the second exponential function into one, which is simpler. Another drawback is the existence of the constant . We can see from Figure 1(b) that the values of the Nakagami-*m* quantile function gradually decrease to zero when the values of Nakagami-*m* CDF decrease from 1 to 0. However, for (11), which means that the constant blocks the tendency, and this flaw can be observed in the simulation results of [5, 6]. The expression we propose meets this requirement , and the excellent performance will be presented in Section 5.

**(a)**

**(b)**

##### 4.2. Approximation for Nakagami-*m* Quantile Function Based on GO-QSSA

###### 4.2.1. Generalized Opposition-Based Quantum Salp Swarm Algorithm

SSA is inspired by salps’ swarming behavior [14] and has been used in several aspects including wireless sensor networks [15], feature selection [16–18], parameter estimation [19], and clustering [20]. Salps usually exist in the form of a swarm called salp chain in deep oceans. There are two groups called the leader and followers in the population according to their positions in the salp chain. The first one of the chain is regarded as the leader; yet, the remaining salps of the chain act as followers. The leader guides the swarm to move in the direction to the food source, whereas each follower follows the one in the front of itself. The principle of SSA is quite simple and easy to understand. But the original SSA easily traps into local optimal solutions, which leads to inexact solutions. Another flaw of SSA is that the designed equations rely on the maximum number of iterations severely, which means that SSA needs the maximum number of iterations to be large enough so that it can explore and exploit well in the search space.

According to the theory of SSA and quantum computing [21–25], the GO-QSSA is proposed in this paper. The GO-QSSA not only takes advantage of salps’ swarm intelligence but also designs brand new mathematical equations to overcome the drawbacks of SSA.

There is a salp swarm containing salps which possess their own positions. Thus, we define the *r*th salp’s position at *t*th iteration as mapped by the *r*th salp’s quantum position , where and is dimension of the salp’s position, and it is also the number of coefficients of the proposed expression. The mapping function from a salp’s quantum position to its position is defined as follows:where , , and , respectively, denote the lower and upper bounds of the *d*th dimension.

The position of each salp denotes a *D*-dimensional potential solution, and we use the fitness function to evaluate the quality of a potential solution at the *t*th generation. Meanwhile, the position with best fitness value is regarded as the position of the food source. Here, we use and to denote the food source’s position and quantum position, respectively, and we use to denote the best quantum position that the *r*th salp acquired by the end of *t*th iteration.

According to the theory of quantum computing, the quantum rotation angles are updated through quantum positions of the last generation, and then, we use the updated quantum rotation angles to update and get new quantum positions. Meanwhile, considering the structure of the salp chain, the leader should guide the swarm to move in the direction to the food source and the followers follow the one in front of themselves, we propose three strategies to update their quantum rotation angles. The first and second strategies are designed for the leader and can be respectively expressed as

In equations (13) and (14), and are parameters controlling quantum rotation angles, denotes a random number which is uniformly distributed in , denotes a standard normal random number distributed in , and is an random integer generated from . The leader randomly selects (13) or (14) to update its quantum rotation angles . These two strategies are selected by the leader with a probability of 50%. As for followers, we use the third strategy to update their quantum rotation angles. Equations (15) and (16) give the third strategy, we first compute an auxiliary quantum rotation angle with (15), and then the quantum rotation angle is the average value of the auxiliary quantum rotation angle and the quantum rotation angle of the salp in front of current salp .

Here, we have obtained the quantum rotation angles of all salps, and the next step is to update their quantum positions with equation (17), in which is the absolute value symbol:

Opposition-based learning was first proposed to improve algorithms in machine intelligence by Tizhoosh [26], and after that, Wang et al. [27] designed a new scheme called generalized opposition-based learning (GOBL). Researchers have applied these two schemes to several heuristic optimization algorithms since their proposing and have succeeded in improving performance of these algorithms [28–30]. Inspired by GOBL, we propose the corresponding concept in quantum space, GOBL in quantum space, and use it to improve QSSA later. Related definitions are given below.

*Definition 1 (opposite point in real-number space). *Let be the *i*th point in *D*-dimensional space and , then the opposite point of is defined by , in which is defined as follows:

*Definition 2 (GOBL in real number space). *Let be the position of *i*th individual in the population at the *t*th generation, is *d*th dimension of , then the opposite position of is defined by , in which is defined as follows:where is a random number uniformly distributed in , and , respectively, denote the minimum and maximum of *d*th dimension that the the *i*th candidate solution ever reached by the end of *t*th generation. Evaluating and , if the fitness of is better than the fitness of , the position of *i*th individual remains unchanged. Otherwise, the position of *i*th individual is changed to .

*Definition 3 GOBL in quantum space. *Let be the corresponding quantum position of *i*th individual’s position , , thus the opposite quantum position of is defined by , in which is defined as follows:where, is a random number uniformly distributed in . Here, and , respectively, denote the minimum and maximum of *d*th dimension that the *i*th individual’s quantum position ever reached by the end of *t*th generation and the corresponding opposite position of is denoted as . Then, evaluating and , if the fitness of is better than the fitness of , the quantum position of *i*th individual remains unchanged. Otherwise, the quantum position of *i*th individual is changed to . In GO-QSSA, each individual operates GOBL with a certain probability of .

###### 4.2.2. GO-QSSA for Approximating Nakagami-*m* Quantile Function

In this section, we give the method to approximate Nakagami-*m* quantile function with equation (8) and GO-QSSA. In Section 3.2, we have modelled the approximation problem as optimizing the coefficients of equation (8), and our aim is to get proper coefficients so as to minimize the RMSE, which is a continuous optimization problem. In the GO-QSSA, quantum positions of all salps are updated through their corresponding strategies until the terminating condition is satisfied. The GO-QSSA is summarized as follows. *Step 1*. Set the parameters of GO-QSSA: the maximum number of iterations , the population size , , , and GOBL probability . *Step 2.* Randomly generate the quantum position of each salp , and record the initial quantum position as the local best quantum position of the *r*th salp acquired so far. At beginning, set . *Step 3.* Evaluate each salp’s position with the fitness function . Select the position with the best fitness value as the position of food source and record the quantum position of food source . *Step 4.* For the leader, update its quantum rotation angles and quantum position according to the first or second strategy. For the followers, update their quantum rotation angles and quantum positions according to the third strategy.* **Step 5*. For each salp, generate a number which is randomly distributed in . If , then compute its opposite quantum position of the updated quantum position with GOBL scheme in quantum space and record its corresponding opposite position and position , if , GOBL is not operated for the *r*th salp. *Step 6.* Evaluate each salp’s position and evaluate *r*th salp’s opposite position if the *r*th salp operated GOBL. If the fitness of is better than the fitness of , choose as *r*th salp’s new quantum position. Otherwise, choose as *r*th salp’s new quantum position. Update the quantum position and position of food source if there is a better salp’s position. Update . *Step 7*. Stop and output the position of food source if . Otherwise, return to Step 4; meanwhile, set .

#### 5. Simulation Results and Discussion

Experiments under several combinations of and are conducted, and the simulation results show the excellent performance of our approximation against that of the existing approximation methods.

Here, we give the results optimized by total four algorithms—PSO, ABC, SSA, and GO-QSSA. There are numerous heuristic algorithms in the literature, including cuckoo search algorithm [31, 32], firefly algorithm [33], teaching-learning-based optimization [34], biogeography-based optimization [35], and so on. We choose PSO because it is the most well-known and classic swarm intelligent optimization algorithm, and it has been successfully applied to near all engineering domains. We choose ABC because it is the most efficient algorithm in the publishing literatures when solving this approximation problem. We choose SSA because it is a novel optimization algorithm and has great potentiality to improve the approximation problem of the Nakagami-*m* model. Compared to PSO, ABC, and SSA, it is easy to understand the superiority of the proposed algorithm in approximating Nakagami-*m* quantile function. Additionally, PSO is a well-known and classic algorithm based on particle swarm intelligence, ABC is a typical algorithm inspired by the swarm intelligence of bees, and SSA is the source of inspiration of GO-QSSA. The common parameters of the four algorithms are and , which are set as 40 and 5,000, respectively. Here, we set to be large enough to ensure a stable solution is achieved for each algorithm. We set the inertia constant as 1.0, and the two acceleration constants are both equal to 2.0 in PSO, and in GO-QSSA, we set , . The parameters of ABC and SSA are the same as in [6, 14], respectively.

##### 5.1. Performance of GO-QSSA

In this paper, we conduct independent simulations 100 times for each algorithm under each combination of fading and scale parameters. Convergence curves are shown in Figures 2(a)–2(d) when scaling parameter is set as 1 and fading parameter is set as 2.5, 3.0, 3.5, and 4.0, respectively, and Figures 3(a)–3(d) give convergence curves when scaling parameter is set as 2 and fading parameter is set as 7.0, 8.0, 9.0, and 10.0, respectively. Tables 1 and 2, respectively, give the corresponding statistical data including the minimum, maximum, average, and standard deviation of RMSEs.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

As we can see in these figures, PSO and ABC quickly trap into local optimal solutions. There is some probability that SSA can escape from some local optimal solutions, but SSA fails to compute the global optimal solutions in the end. The statistical data in Tables 1 and 2 tell us that the minimum, maximum, and average of RMSEs computed by GO-QSSA are almost on the same level under each combination of fading and scaling parameters, and standard deviation values of RMSEs are much smaller than those of other three algorithms, which means that GO-QSSA can steadily find the best solution. It is obvious that GO-QSSA prevails over ABC, PSO, and SSA in both convergence speed and convergence accuracy, that is to say, GO-QSSA is capable of balancing its exploring and exploiting abilities better than ABC, PSO, and SSA so that it has the ability to jump out of local optima and exploit more in the search space.

##### 5.2. Efficiency of the Proposed Expression

Another main aspect in the approximation problem is the expression for approximating Nakagami-*m* quantile function. Here, the comparisons of the performance of the proposed expression with that of the existing expressions are shown in Tables 3 and 4, from which we can see that, compared with the results in [4–6], the proposed expression decreases the RMSEs by roughly one order of magnitude in most cases () and even higher in some cases (). Lower RMSEs of the proposed expression mean that the proposed expression is capable of fitting the curve of Nakagami-*m* quantile function much better than those expressions in [4–6]. Meanwhile, Tables 3 and 4 reveal that for a certain scaling parameter such as or , the RMSE increases when the fading parameter grows, which means that the difficulty of fitting the curve of Nakagami-*m* quantile function increases as the fading parameter grows.

Considering the simulation results and analyses above, we finally compute the coefficients of the proposed expression with the assistance of GO-QSSA. The coefficients of the proposed expression computed by GO-QSSA under various scaling parameters ( and ) and fading parameters varying from 0.5 to 10 are shown in Tables 5 and 6. In order to show the results more intuitively, comparisons between the approximated values of the proposed expression with coefficients shown in Tables 5 and 6 and the exact values under different fading and scaling parameters are depicted in Figures 4(a)–4(d), from which we can see that our approximation has much lower biases compared with the results in [4–6] under the same parameters. That is to say, our approach to approximate Nakagmi-*m* quantile function is more efficient. Additionally, we give the plot of mean approximating error against Nakagami-*m* CDF value under different parameters in Figure 5. Here, the mean approximating error means the average of the difference between the exact and approximating values of Nakagami-*m* quantile function. We can see from Figure 5 that the approximating error is higher when and when is close to 0 or 1, which means that it is more challenging to approximate Nakagami-*m* quantile function when the scaling parameter increases and when Nakagami-*m* CDF value is close to 0 or 1.

**(a)**

**(b)**

**(c)**

**(d)**

#### 6. Conclusions

In this paper, a simple and efficient expression for approximating the Nakagami-*m* quantile function with higher accuracy containing four coefficients is proposed by dint of three elementary functions including an inverse hyperbolic tangent function, an exponential function, and a power function. Then, we combine SSA with quantum swarm intelligence and propose QSSA. Additionally, we expand generalized opposition-based learning strategy from real-number space to quantum space, and we use GO-QSSA to compute coefficients of the proposed expression. Finally, enormous simulation results under different combinations of fading parameters and scaling parameters validate the accuracy of the proposed expression and effectiveness of GO-QSSA.

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 61571149), China Postdoctoral Science Foundation (Grant no. 2015T80325), and Fundamental Research Funds for the Central Universities (HEUCFP201808).