Computation efficiency optimization in UAV-enabled mobile edge computing system with multi-carrier non-orthogonal multiple access

In this paper, we study the issue of fair resource optimization for an unmanned aerial vehicle (UAV)-enabled mobile edge computing (MEC) system with multi-carrier non-orthogonal multiple access (MC-NOMA). A computation efficiency (CE) optimization problem based on the max-min fairness principle under the partial offloading mode is formulated by optimizing the subchannel assignment, the local CPU frequency, and the transmission power jointly. The formulated problem belongs to the non-convex mixed integer nonlinear programming (MINLP), that is NP-hard to find the global optimal solution. Therefore, we design a polynomial-time algorithm based on the big-M reformulation, the penalized sequential convex programming, and the general Dinkelbach’s method, which can choose an arbitrary point as the initial point and eventually converge to a feasible suboptimal solution. The proposed algorithm framework can be also applied to computation offloading only mode. Additionally, we derive the closed-form optimal solution under the local computing only mode. Simulation results validate the convergence performance of the proposed algorithm. Moreover, the proposed partial offloading mode with the CE maximization scheme outperforms that with the computation bits (CB) maximization scheme with respect to CE, and it can achieve higher CE than the benchmark computing modes. Furthermore, the proposed MC-NOMA scheme can attain better CE performance than the conventional OFDMA scheme.

MCC is difficult to meet the computing service requirements of various real-time mobile applications [2]. Recently, Internet of Things (IoT) has been applied to many aspects, which can connect things via networks [3]. In contrast to the conventional MCC, mobile edge computing (MEC) has potential wide application ranges for IoT as a newly emerging technology [4]. The key idea of MEC is to decentralize computing servers at the wireless network edge closing to UEs so as to improve the computation capability of UEs and reduce the corresponding task delay and consumed energy greatly.
On the one hand, non-orthogonal multiple access (NOMA) is considered as a promising multiple access scheme in the future wireless communication systems [5]. The fundamental concept of NOMA is to serve more than one UE in one orthogonal resource block (ORB) by using superposition coding and multi-user detection technologies [6]. For instance, power-domain NOMA assigns different power levels to the UEs at the transmitter and exploiting successive interference cancelation (SIC) technology to decode overlapping signals at the receiver [7]. In [8], a NOMA-based resource allocation algorithm was designed for a cluster-based cognitive industrial IoT. This is very different from orthogonal multiple access (OMA) schemes, e.g., time division multiple access (TDMA), orthogonal frequency division multiple access (OFDMA), which service one UE at most in a time-domain ORB and a frequency-domain ORB, respectively. However, the system complexity and decoding delay will become high if all users are multiplexed in one ORB. To this end, multi-carrier NOMA (MC-NOMA), i.e., the combination of NOMA and OFDMA, can be employed to meet the needs of massive connectivity and improve the transmission efficiency in MEC systems [9][10][11][12], where all users are divided into several NOMA groups and one NOMA group can occupy one ORB. In [9], a MEC-NOMA technique was proposed to minimize the consumed energy of MEC users, where the users can simultaneously upload their tasks via the uplink NOMA protocol. In [10], the total computation overhead minimization problem was investigated by coalition formation game theory in a MC-NOMA based MEC system. Moreover, [11] and [12] provided their separate theories and simulation results to the energy minimization in MC-NOMA based MEC systems, which show the application of MC-NOMA in MEC can improve the system performance greatly when compared to the OMA counterpart.
On the other hand, unmanned aerial vehicle (UAV) has attracted a great quantity of research interests in the wireless communication. The benefit of UAV as the aerial base station (BS) is great when compared with the traditional ground BS, such as flexible movement, high probability of line-of-sight (LoS) channel, and low cost and easy deployment. Therefore, MEC servers can be carried to UAV for the purpose of providing flexible computing services [13][14][15][16][17][18]. In [13], the Lyapunov optimization framework was invoked for minimizing the average weighted consumed energy of a UAV-MEC system. In [14], the resource optimization for maximizing computation rate was firstly studied in a wireless powered UAV-MEC system. In [15], user association and path planning were studied in a UAV-MEC system with energy restriction. In [16], a double deep Q-network based algorithm was developed to maximize the system reward of a UAV-MEC network. Moreover, [17] studied the sum-power optimization for a MEC system with several UAVs, which was addressed by solving the decomposed subproblems iteratively. Besides, [18] minimized the consumed energy of a UAV-MEC system, which can enable fairness among all the users.
It can be seen that the delay minimization, the energy consumption minimization, and the computation bits (CB) maximization are the main focus of resource allocation in MEC systems. However, with the rapid growth of communication business, the communication industry is becoming a major energy consumer. Against this background, it is critically important to consider green communication so as to improve the efficiency of resource utilization per joule. Particularly, it is suggested that future MEC pursue the maximization of computation efficiency (CE), i.e., the ratio of the total CB to the associated energy consumption [19][20][21][22][23]. In [19], the max-min CE optimization problem was studied in a MEC system powered by wireless energy. In contrast to [19,20] maximized the weighted sum of CE in a MEC system. While TDMA is considered in both [19] and [20,21] and [22] extended the CE maximization frameworks to MEC networks under NOMA and OFDMA, respectively. Furthermore, [23] aimed to maximize the CE of a multi-UAV based MEC system. Motivated by the above discussions, we study the CE optimization problem in a UAV-MEC system with MC-NOMA, which is oriented towards computation-efficient massive connections and flexible deployment in the future MEC network.
The rest of this paper is organized as follows. The system model and problem formulation are presented in Section 2. The proposed CE optimization algorithm was detailed in Section 3. Two benchmark computing modes are discussed in Section 4. The relevant simulation results are provided in Section 5. Finally, the conclusion of this paper is made in Section 6.

Method
In this paper, we focus on a computation-efficient UAV-enabled and MC-NOMA based MEC system, where a UAV is fitted with a MEC server to supply computing services for the ground UEs via the LoS connections. Besides, MC-NOMA is employed for the uplink transmissions between all the UEs and the UAV, in which all the UEs are divided into several user groups without repetition, and the UEs in the same group can adopt uplink NOMA to upload their tasks simultaneously to the UAV-MEC server in one ORB.
On this basis, we formulate the CE optimization problem based on the max-min fairness principle under the partial offloading mode by optimizing the subchannel assignment, the UE's CPU frequency, and the transmission power jointly. Considering that the formulated problem belongs to the non-convex mixed-integer nonlinear programming (MINLP), we design a two-layer iterative algorithm based on the big-M reformulation, the penalized sequential convex programming (SCP), and the general Dinkelbach's method to enjoy a polynomial-time complexity. The proposed algorithm can start from an arbitrary point and eventually converge to a feasible suboptimal solution. Moreover, the other two typical computing mode are discussed as benchmarks, including the local computing only mode and the computation offloading only mode. Specially, the closed-form best optimal solution is derived under the local computing only mode, while a suboptimal solution is obtained using the above algorithm under the computation offloading only mode.
Finally, simulation results confirm that the proposed algorithm achieves good convergence performance and is robust to random initial points. In addition, the proposed CE maximization framework has obvious computation-efficient advantages over the CB maximization framework under the partial offloading mode. Furthermore, the CE performance of the partial offloading mode is superior to the benchmarks computing modes.
Moreover, the MC-NOMA scheme exhibits higher CE than the OFDMA scheme in the proposed UAV-MEC system.

System model and problem formulation
As illustrated in Fig. 1, we focus on a UAV-MEC system where one rotary-wing UAV with an embedded MEC server provides computing services for K ground UEs denoted by the set k ∈ K = {1, 2, . . . , K}, in which the UAV and all K UEs are each equipped with one single-antenna. We construct a three-demission (3D) Cartesian coordinate model, where the UAV and the kth UE are located at q u =[ X, Y , H] and q k =[ x k , y k , 0] , ∀k ∈ K, respectively. As a result, the distance between the UAV and the kth UE is given by

Communication model
In the proposed UAV-MEC system, the overall system bandwidth W is partitioned into N orthogonal subchannels denoted by the set n ∈ N = {1, 2, . . . , N} satisfying n∈N W n = W , where W n is the bandwidth of the nth subchannel. Besides, each subchannel is the frequency-flat block fading channel, i.e., the subchannel remains unchanged during one time block with duration T (T > 0). We assume the UAV knows the perfect channel state information (CSI) of all subchannels as well as the computation information of all the UEs. Moreover, it is assumed that all subchannels are each dominated by LoS links since the altitude of the UAV is usually much higher than that of the UEs. Hence, by adopting the commonly used free-space path loss model [13][14][15][16][17][18], the channel power gain on the nth subchannel between the kth UE and the UAV can be modeled as where μ 0 represents the channel power gain at the reference distance d 0 = 1 m. Particularly, one time block with duration T is composed of the following three stages. In the first stage, MC-NOMA is invoked for the uplink transmissions between the UEs and the UAV, i.e., each orthogonal subchannel can accommodate multiple users by employing NOMA. Specially, all the UEs are divided into different groups to upload their tasks simultaneously to the UAV-MEC server, where different groups access different subchannels and the UEs in the same group occupy one subchannel via NOMA. In the second stage, the UAV exploits SIC technique to decode the uploaded tasks from the UEs and the UAV-MEC server computes the decoded tasks. Finally, the UAV transmits the task results computed by the UAV-MEC sever to all the UEs in the third stage. Besides, each UE can perform its local computing throughout the entire time block duration T. Similar to [19][20][21][22][23], we only consider the first stage due to the strong computational capability of the UAV-MEC sever as well as the small data size of the computed results.

Local computing
For the local computing at the UEs, the total CB L loc k and total energy consumption E loc k of the kth UE during the entire time block T can be, respectively, expressed as [19][20][21][22][23] where C k denotes the CPU cycles of computing 1-bit of input task of the kth UE, f k (cycles/second) is the CPU frequency at the kth UE, and γ k is the CPU capacitance coefficient of the kth UE.

Computation offloading
To enable parallel execution, the proposed UAV-MEC system adopts the partial computation offloading mode, where the task of each UE can be divided into two parts, i.e., one part computed locally at the UE and another part uploaded to the UAV-MEC sever. Let ρ k,n denote the subchannel assignment indicator, ∀k ∈ K, n ∈ N . Specially, ρ k,n = 1 if the kth UE is assigned to the nth subchannel, otherwise ρ k,n = 0. According to the uplink NOMA protocol, the received signal at the UAV on the nth subchannel can be expressed as where s k,n represents the transmitted signal of the kth UE over the nth subchannel satisfying E{|s k,n | 2 } = 1, p k,n is the corresponding transmission power, and z n stands for the additive white Gaussian noise (AWGN) at the nth subchannel, i.e., z n ∼ CN (0, n 0 W n ), in which n 0 denotes the single-sided power spectral density of AWGN. On each subchannel, the SIC decoding order of the uplink NOMA is always from the strong UE to the weak UE so as to avoid consuming significant transmission power for compensating the path loss [24]. Then, the uplink offloading rate of the kth UE on the nth subchannel is given by where S k,n = K|g i,n < g k,n , i = k denotes the set of the UEs in the set K whose channel power gain are worse than that of the kth UE. Therefore, the number of offloaded bits of the kth UE on the nth subchannel can be calculated as and the corresponding consumed energy is where ξ is the amplifier coefficient at the UE, and p c is the fixed circuit power consumption for the computation offloading process at the UE. According to [19] and [20], the CE of the kth user can be quantified as

Problem formulation
With the purpose of guaranteeing fairness for all the UEs, the subcarrier assignment {ρ k,n }, the local CPU frequency {f k,n }, and the transmission power {p k,n }, ∀k ∈ K, n ∈ N , are jointly optimized to maximize the minimum CE among all the UEs. Therefore, the CE optimization problem based on the max-min fairness principle can be formulated as max {ρk,n,pk,n,fk,n} min where R min k in the constraint C 1 is the minimum CB per second required by the kth UE, P max k in the constraint C 2 denotes the maximum power consumption of the kth UE, f max k in the constraint C 3 denotes the maximum CPU frequency of the kth UE, the constraint C 4 indicates that the nth subchannel can support at most M n UEs, and the constraints C 5 and C 6 ensure that one UE can occupy at most one subchannel.

Proposed computation efficiency optimization algorithm
It is clearly that the problem in (8) belongs to the non-convex MINLP, that is NP-hard to find the global optimal solution in general. In this section, we design a polynomialtime iterative algorithm based on the big-M reformulation, the penalized SCP, and the generalized Dinkelbach's method to obtain a feasible suboptimal solution of (8).

Big-M reformulation
Due to the product term ρ k,n p k,n in (8), it is not easy to design an efficient algorithm for solving (8). Fortunately, we can use the big-M reformulation to decompose the product terms [25]. With the help of big-M reformulation, we first introduce auxiliary variables p k,n = ρ k,n p k,n , ∀k ∈ K, n ∈ N . Then, all the product terms ρ k,n p k,n in (8) will be replaced withp k,n , which yields Afterwards, we add the following constraints into (8): Moreover, we can rewrite the binary constraint C 6 in its equivalent continuous form as follows: Therefore, we can reformulate (8) as the following equivalent problem: where = ρ k,n , p k,n ,p k,n , f k,n , ∀k ∈ K, n ∈ N .
The big-M reformulation linearizes the term ρ k,n p k,n , such that ρ k,n is separated from p k,n , which helps to design an efficient CE optimization algorithm. Moreover, if the constraints C 6a , C 6b , and C 7 ∼ C 10 , are satisfied, thenC 1 is equivalent to C 1 , andC 2 is also equivalent to C 2 . Hereto, the original MINLP problem in (8) has been equivalently transformed into a more tractable form in (13). However, the transformed problem in (13) still belongs to the non-convex max-min fractional programming (MMFP), which can not be addressed directly. To this end, we combine fractional programming theory with sequential convex programming (SCP) [26,27].

Sequential convex programming
In particular, the core idea of SCP is to convert a non-convex optimization problem into a series of convex optimization problems and solve them iteratively, where the non-convex terms in each iteration will be substituted by appropriate inner convex terms [28].
We note thatR k,n in (13) can be expressed as the subtraction of two concave functions R k,n,1 andR k,n,2 , i.e., At the tth iteration (t ≥ 1) of SCP, we can obtain an upper boundR (t) k,n,2 ofR k,n,2 by using its first-order Taylor expansion, which is expressed as is the value ofp i,n at the (t − 1)th iteration of SCP. Similarly, we can obtain a lower upperρ (t) k,n of ρ 2 k,n by using its first-order Taylor expansion, which is expressed as where ρ (t−1) k,n is the value of ρ k,n at the (t − 1)th iteration of SCP. Furthermore, since the function min{·} is increasing, we can obtain a lower boundη (t) k ofη k at the tth iteration of SCP, i.e., As a result, a convex approximation of (13) at the tth iteration of SCP can be cast as k,n ≤ 0, ∀k ∈ K, n ∈ N ,

Generalized Dinkelbach's method
It can be observed that all the constraints in (18) are convex, i.e., the feasible region of (18) is a convex set. In addition, for eachη k , ∀k ∈ K, its numerator is non-negative and concave, meanwhile its denominator is positive and convex. Therefore, the problem in (18) belongs to the convex MMFP, which can be efficiently addressed by the generalized Dinkelbach's method [26].
Let λ * and * denote the optimal objective and the optimal solution of (18), respectively, then Moreover, * is the optimal solution of (18) if and only if the following condition is satisfied: [26] max min k∈K Based on (20), we can transform (18) into its associated subtractive form as max min where (21) is the problem to be solved at the qth iteration (q ≥ 1) of the generalized Dinkelbach's method, and λ (q−1) is the corresponding iterative parameter at the (q − 1) iteration. At the qth iteration, we compute the optimal solution (q) of (21) for the fixed λ (q−1) , then we update λ (q) by using the obtained (q) , i.e., By solving (21) iteratively and setting a very small iteration tolerance > 0, -optimal solution * = (q) of (18) is achieved if the following condition holds: [26] min k∈K Since the objective function in (18) is not smooth, we introduce a new auxiliary variable given as Accordingly, the equivalent form of (18) can be expressed as: max ,u u s.t.C 1 ,C 2 , C 3 , C 4 , C 5 , C 6a ,C 6b , C 7 ∼ C 10 , Based on the above analysis, a two-layer iterative algorithm consisting of the outer iterative procedure based on the SCP and the inner iterative procedure based on the generalized Dinkelbach's method can be employed to obtain a feasible suboptimal solution of (13). Nevertheless, the above algorithm requires a feasible point of (13) as the initial point. If the above algorithm starts from a feasible point, then all the generated points will be feasible; otherwise, it may fail at the first iteration due to the infeasibility [28]. However, finding a feasible point in the non-convex feasible region of (13) is very hard. To avoid this, we will integrate the penalized SCP [29] with the framework of the above algorithm, which can start from an arbitrary point.

Penalized sequential convex programming
To ensure that a sequence of feasible solutions can be generated, the abovementioned SCP needs a feasible initial point. Generally, an infeasible initial point misleads intermediate solutions obtained by the iterative process, which often causes an incorrect local solution. However, it is usually NP hard to find a feasible initial point for non-convex problems, e.g., the CE optimization problem in (13). To avoid the above initialization requirement, it is a good choice to adopt the penalized SCP. Specially, auxiliary variables are introduced to relax all non-convex inequality constraints, and penalty parameters are added into the objective function, so that we can randomly generate the initial point and then feasible solutions can be gradually obtained with the increasing of iterations.
It can be found that only the constraintsC 1 and C 6b are non-convex in (13). To this end, we relax the constraintC 1 via slack variables s 1,k (s 1,k ≥ 0, ∀k ∈ K) and the constraint C 6b via a slack variable s 2 (s 2 ≥ 0). Then, the constraintsC 1 andC 6b in (25) are, respectively, reformulated as In order to minimize the violations of non-convex constraints, we subtract the sum of the slack variables from the objection function of (25). Therefore, (25) can be reformulated as where (27) is the problem to be solved at the tth iteration of the penalized SCP, and τ (t−1) is the corresponding penalty factor at the (t − 1)th iteration and is updated for the next iteration by using the following equation: Here, τ max denotes the upper bound of τ (t) so as to avoid numerical problems if τ (t) is too large and provide convergence if a feasible region is not found, and μ is the predefined increasing factor [29]. Moreover, existing optimization tools, like YALMIP [30], CVX [31], CVXQUAD [32], can be used to solve the convex problem in (27). In summary, we propose a CE optimization algorithm based on the big-M reformulation, the penalized SCP, and the generalized Dinkelbach's method, which is presented in Algorithm 1.

Theorem 1 Algorithm 1 can converge to a local optimal solution of (13).
Proof See Appendix A.

Computational complexity analysis
In Algorithm 1, the computational complexity is mainly determined by solving problem (27) L 1 L 2 times, where L 1 and L 2 denote the required iterations of the inner loop and outer loop, respectively. In addition, (27) is solved by using standard optimization tools based on the interior-point method whose complexity is O (4KN + N + 2) 3.5 log(1/ε) , where O (·) denotes the big-O notation and ε represents the solution accuracy [33]. Thus, the proposed Algorithm 1 has an acceptable polynomial time computational complexity, i.e., O L 1 L 2 (4KN + N + 2) 3.5 log(1/ε) .

Computing modes as benchmarks
In this section, we discuss the other two computing modes as benchmarks, namely the local computing only mode and the computation offloading only mode. The local computing only mode requires that all the UEs only perform their local computing but never offload their tasks to the UAV-MEC server, while the computation offloading only mode Set t = t + 1.

Local computing only mode
Under the local computing only mode, the CE of the kth UE can be expressed as The corresponding CE optimization problem is formulated as In (30), we find that the optimization variables {f k , k ∈ K} are decoupled. Moreover, the objective function consists of the "min" function and K independent CE expressions, i.e., {η L,k , k ∈ K}. Without loss of optimality, we can maximize each CE expression η L,k to maximize the objective function of (30). Thus, (30) can be decomposed into K sub-problems, in which the kth sub-problem is given by By setting the first order derivative of η L,k to be zero, i.e., the zero point of ∂η L,k /∂f k can be calculated as The numerator of ∂η L,k /∂f k is a concave and monotone decreasing function with respect to f k ≥ 0, while the denominator of ∂η L,k /∂f k is positive for f k ≥ 0. Thus, the closed-form optimal solution of (31) can be derived as As a result, the best optimal solution of (30) is given by Theorem 2 The problem in (30) has many optimal solutions with the same objective value.
Proof See Appendix B.

Remark 2
Since (30) has many optimal solutions with the same objective value, "the best" means that the optimal solution in (35) can make the CEs of all the UEs as large as possible.

Computation offloading only mode
Under the computation offloading only mode, the CE of the kth UE can be expressed as The corresponding CE optimization problem is formulated as max {ρk,n,pk,n} min Since (37) is similar to (8), the framework of Algorithm 1 can be applied to solve (37). To avoid repetition, the detailed solution to (37) is omitted here.

Simulation results and discussions
In this section, we provide the relevant simulation results to first demonstrate the convergence of the proposed Algorithm 1 and then evaluate the performance of the proposed partial offloading scheme by using the local computing only scheme, the computation offloading only scheme, and the OFDMA scheme as benchmarks. In the simulations, we consider a UAV-MEC system as illustrated in Fig. 2, where K = 8 UEs are randomly located in a squared area of size 1000 m ×1000 m, and the UAV is located at Fig. 2 The locations of the UEs and the UAV in x-y axes q u =[ 0, 0, 100] [23]. To reduce the decoding complexity and time delay in one NOMA group, the maximum number of the UEs supported by each subchannel is considered as two, i.e., M n = 2, ∀n ∈ N [34]. Besides, the maximum power consumption of the kth UE is P max k = P max , ∀k ∈ K, where P max is defined as the maximum local power consumption. For convenience, the normalized CE and CB are used as the performance indicators, and their units are, "Bits/Joule/Hz" and "Bits/Hz", respectively. Except where otherwise stated, we present the remaining simulation parameters in Table 1, which are based on the typical settings in [21][22][23]. Figure 3 plots the convergence curves of the iterative objective value (t) and the sum of penalty values π (t) 1 for the proposed Algorithm 1 under the partial offloading mode, where the maximum local power is set as P max = 0.2W. As seen from Fig. 3, when the number of iterations increases, the sum of penalty values π (t) 1 decreases accordingly. Since π (t) 1 denotes the sum of violations of non-convex constraints, the generated solutions by Algorithm 1 tend to be feasible for (8). In fact, Algorithm 1 can be interpreted as an infeasible method for the original MINLP problem in (8). Specially, Algorithm 1 may generate a sequence of infeasible solutions if it starts from an infeasible point but can converge to a feasible suboptimal solution of (8) at last. Moreover, the proposed Algorithm 1 has a fast convergence rate because the iterative objective value (t) needs only 5 iterations to be stable and converge to a -optimal solution for the specified precision . Furthermore, the average iterations of Algorithm 1's inner loop and outer loop are 6.08 and 6.49, respectively. The above observations demonstrate the convergence property and the convergence efficiency of Algorithm 1. Table 2 evaluates the impact of random initial points on the CE performance of Algorithm 1 under the partial offloading mode, where P max = 0.2W. The random initial points for Algorithm 1 are given by ρ  Table 2, with these random initial points, Algorithm 1 will eventually converge to one of the three different local optimal values, denoted by, η 1 = 10.455 bits/Joule/Hz with the probability 11%, η 2 = 11.487 bits/Joule/Hz with the probability 15%, and η 3 = 12.455 bits/Joule/Hz with the probability 74%, respectively. Moreover, the probability of converging to η 3 is 74% for one run time, which means that Algorithm 1 has a 98.4% probability of converging to the best local optimal values within 3 run times.  In other words, Algorithm 1 is robust to different initial points due to the high probability of converging to the best local optimal values with very few run times. In Figs. 4 and 5, we compare the proposed partial offloading mode with the CE maximization scheme referred as "MaxCE" and the CB maximization scheme referred as "MaxCB" in terms of CB and CE with varying the maximum local power consumption P max , respectively. In particular, the CB maximization scheme is to maximize the minimum CB among all the UEs, which is implemented by using the simplified Algorithm 1 without the inner loop, i.e., λ (q) = 0, ∀q. From the above two figures, the "MaxCE" scheme has both CB and CE in common with the "MaxCB" scheme when P max is not so large, i.e., P max ≤ P th max . Due to the limitation of P max , the resource allocation strategies of the two schemes are identical when P max is small, which leads to the above result. In the case of small P max , the CE and CB of the two schemes increase with P max since the CB increases faster than the power consumption. But on the other hand, when P max > P th max , although the "MaxCB" scheme has larger CB than the "MaxCE" scheme, the "MaxCE" scheme has obvious advantages on the CE compared with the "MaxCB" scheme from Fig. 5. This is because the "MaxCB" scheme maximizing the CB will consume more inefficient power to achieve higher CB but the "MaxCE" scheme will still adopt the solution maximizing the CE instead of consuming inefficient power when P max > P th max , i.e., the solution of the "MaxCE" scheme is not affected by the increase of P max . Consequently, both CE and CB of the "MaxCE" scheme keep unchanged when P max > P th max . Figure 6 presents the CE versus the maximum local power consumption P max under different computing modes with the proposed CE maximization frameworks. From Fig. 6, it can be found that the partial offloading mode has a much higher CE performance than the other two computing mode, especially when compared to the local computing only mode. The reason is that the partial offloading scheme can offer the flexibility to adjust the resource allocation for the task offloading and the local computing. However, the local computing only scheme allows the UEs to perform their local computing only regardless of the channel conditions, and the computation offloading only mode completely uploads the tasks of the UEs to the UAV-MEC server even if the channel conditions are bad, so their CE performances are greatly compromised. Besides, the computation offloading only mode exhibits better CE performance than the local computing only mode since the former can yield better computation-efficient benefits than the latter when P max is not so small. Figure 7 shows the max-min fairness with the proposed CE maximization framework under different computing modes, where the minimum CE, the average CE, and the best CE of all the UEs are compared, and the maximum local power is set as P max = 0.2W. With respect to the above three performance indicators, the proposed partial offloading mode has the best performances, while the local computing only mode has the worst  performances. The reason for explaining this is similar to the simulation result in Fig. 6, which verifies the effectiveness of the proposed scheme in the max-min fairness. Figure 8 illustrates the CE versus the maximum local power consumption P max in the partial offloading mode with the proposed MC-NOMA scheme and the OFDMA scheme. In the OFDMA scheme, each UE can access one subchannel at most and each subchannel can support one UE at most, which is implemented by letting the maximum number of the UEs supported by one subchannel to be one, i.e., M n = 1, ∀n ∈ N . To ensure fair comparison, the number of subchannels in the OFDMA scheme is set as N = 8, but the other simulation parameters are consistent with the proposed MC-NOMA scheme. It can be observed from Fig. 8 that the proposed MC-NOMA scheme is significantly greater than the OFDMA scheme on the CE performance. This is because NOMA can service more users than OFDMA at the same frequency-domain resource block, i.e., the multiplexing gain can be obtained when NOMA is applied.

Conclusion
In this paper, we have studied a max-min fairness based CE optimization problem in a UAV-MEC system with MC-NOMA, where the joint optimization of the subchannel assignment, the local CPU frequency, and the transmission power is considered. Since it is NP-hard to find the global optimal solution of the formulated non-convex MINLP problem, we design a polynomial-time algorithm to obtain a feasible suboptimal solution. In specify, the big-M reformulation is employed firstly to obtain the equivalent continuous form of the original problem. Next, by applying the penalized SCP and the general Dinkelbach's method, a two-layer iterative algorithm is designed, which is able to start from an arbitrary initial point and finally converge to a feasible suboptimal solution. We also analyze the corresponding CE optimization problems in the benchmark computing modes. The relevant simulation results have demonstrated that the proposed Algorithm 1 enjoys a good convergence performance. Moreover, the proposed computation-efficient scheme is superior to the other benchmark schemes. Besides, when compared with the conventional OFDMA scheme, the proposed MC-NOMA scheme can get higher CE.
bounded, Algorithm 2 can gradually converge to a local optimal solution of (13). Next, it is ready to prove the convergence of Algorithm 1. If the penalty factor τ is large enough and the sum of slack variables is very close to zero, then (27) can be rewritten as (25). Therefore, Algorithm 1 is equivalent to Algorithm 2, which indicates that Algorithm 1 will converge to a local optimal of (13). Please refer to [29] for more detail.