Resource Allocation with a Rate Guarantee Constraint in Device ‐ to ‐ Device Underlaid Cellular Networks

: Device ‐ to ‐ device (D2D) communication is a crucial technique for various proximity services. In addition to high ‐ rate transmission and high spectral efficiency, a minimum data rate is increasingly required in various applications, such as gaming and real ‐ time audio/video transmission. In this paper, we consider D2D underlaid cellular networks and aim to minimize the total channel bandwidth while every user equipment (UE) needs to achieve a pre ‐ determined target data rate. The optimization problem is jointly involved with matching a cellular UE (CU) to a D2D UE (DU), and with channel assignment and power control. The optimization problem is decoupled into two suboptimization problems to solve power control and channel assignment problems separately. For arbitrary matching of CU, DU, and channel, the minimum channel bandwidth of the shared channel is derived based on signal ‐ to ‐ interference ‐ plus ‐ noise ratio (SINR) ‐ based power control. The channel assignment is a three ‐ dimensional (3 ‐ D) integer programming problem (IPP) with a triple (CU, DU, channel) . We apply Lagrangian relaxation, and then decompose the 3 ‐ D IPP into two two ‐ dimensional (2 ‐ D) linear programming problems (LPPs). From intensive numerical results, the proposed resource allocation scheme outperforms the random selection and greedy schemes in terms of average channel bandwidth. We investigate the impact of various parameters, such as maximum D2D distance and the number of channels.

shorter than the distance between a device and its corresponding base station (BS), better link quality can be acquired in D2D communications, which is advantageous in terms of high-rate transmissions and reliability. Furthermore, in D2D underlaid networks where DUs share spectrum with cellular user equipments (CUs), spectrum efficiency and network capacity can be further improved [7].
However, in D2D underlaid cellular networks, there is a possibility of severe interference due to resource sharing between CUs and DUs. Therefore, interference management must be applied to effectively handle the interference, and to guarantee a certain degree of quality of service (QoS). In the literature, various power control and resource allocation schemes have been proposed to guarantee QoS requirements in terms of signal-to-interference-plus-noise ratio (SINR) or minimum required data rate [8][9][10][11][12][13][14][15][16][17]. In [8], a power control scheme was studied to maximize the sum of the achievable data rates of a CU-DU pair. However, only the SINR of the CU is guaranteed to be greater than or equal to a pre-determined SINR threshold. In [9], multiple DUs can be matched to a single CU and an optimal energy-efficient power control scheme was proposed while guaranteeing the minimum required data rate of each user equipment (UE). However, only a single channel was taken into account and, hence, the resource allocation or channel assignment problem was not addressed. In [10], a resource allocation scheme was studied while assuming that the transmit power of UEs is fixed. Joint power control and resource allocation schemes were proposed in [11][12][13][14][15][16][17]. A single DU is assigned to only one channel, and vice versa in [11,12], whereas multiple DUs are allowed to be assigned to the same channel in [13,14]. In addition, a DU is allowed to be assigned to multiple channels in [15][16][17]. In addition to power control and resource allocation, mode selection for DUs and relay selections were considered in [18][19][20] and [21,22], respectively. In the above-mentioned studies, a fixed channel bandwidth is allocated to each channel and moreover, a channel is pre-assigned to a CU in advance, and hence, the remaining channel assignment problem is to choose the best DU for each CU.
In most previous sum rate or energy efficiency maximization schemes [8][9][10][11][12][13][14][15][16][17][18][19][20][21][22], even though SINR or minimum data rate per UE was considered as a constraint, the individual data rate may not be guaranteed with a fixed amount of channel bandwidth. When the received SINR is low due to high interference and maximum transmit power limitation, the fixed channel bandwidth may not be wide enough to guarantee the target rate. In D2D underlaid cellular networks, UEs experience diverse interference situations depending on CU-DU matching and power control, so the fixed channel bandwidth assumption may not be beneficial. In mobile communication systems, a proper modulation and coding scheme (MCS) needs to be informed during a scheduling process for successful reception [23,24]. It is also noted that resource information indicating the amount of resource and its location can be informed and, therefore, the channel bandwidth can be adaptively adjusted. Thus, it is optimal to manage the transmit power, channel bandwidth, and channel assignment jointly. However, to the best of our knowledge, there is no research on the resource allocation problem with rate guarantee constraints where the total channel bandwidth is minimized while ensuring that all the UEs can achieve their target data rates.
The contributions of this paper are as follows:  In most previous studies, each CU was assumed to have a preassigned channel, and thus, the channel assignment problem is to find the best DU for each CU in order to maximize the sum rate or energy efficiency. Then, the assigned channel information needs to be informed to the corresponding DUs. However, as indicated in [23,24], the channel information can be informed to CUs as well as DUs and, therefore, the channels can be re-assigned for CUs. Now, the channel assignment problem becomes the triple matching problem with pairs of (CU, DU, channel). This three-dimensional (3-D) channel assignment increases the possibility of choosing a better channel and can bring a performance improvement at the cost of increased computational complexity. To mitigate that extra complexity, a suboptimal algorithm is proposed by using Lagrangian relaxation and a sub-gradient-based iterative algorithm;  We consider a system where the channel bandwidth is not fixed but is adaptively determined based on the triple (CU, DU, channel). As the received SINR varies, the channel capacity changes and, therefore, the channel bandwidth required to guarantee the target data rate needs to be adjusted. In more detail, the minimum channel bandwidth to guarantee the target data rate is derived with a maximum transmit power constraint. Each CU-DU-channel match-up requires different minimum channel bandwidths. In addition, hence, the remaining problem is a 3-D assignment problem to minimize the overall channel bandwidth, while guaranteeing the target data rate;  Since the channel is shared, increasing the transmit power of CU may cause severe interference to DU and vice versa. In this paper, the exact solution to the power control problem involves iterative calculations through, for example, a root-finding algorithm. To reduce complexity, a suboptimal algorithm is proposed where power can be calculated in a closed form. In previous works, the transmit power of the UEs is determined with a fixed channel bandwidth. Even though SINR or minimum data rate per UE is considered as a constraint, some UEs may not meet the constraint, i.e., the feasible set becomes empty due to maximum transmit power limitation and fixed channel bandwidth. Herein, the channel bandwidth can be adaptively determined according to the received SINR to meet the minimum data rate. Figure 1a illustrates a D2D underlaid cellular network where the uplink channel is shared. We consider a single cell where cellular UEs (CUs) and D2D transmitters (DTs) are randomly positioned in the cell. A DT and one D2D receiver (DR) form a D2D UE (DU), and a DR is uniformly distributed within a circle whose origin is at its corresponding DT with a radius of . We assume that every CU and DT are equipped with a single transmit antenna. The BS and every DR are equipped with a single receive antenna. The network is assumed to be provided with orthogonal channels on which channel bandwidths are variable. Every UE is assigned only one channel. Moreover, every channel is either allocated to one CU alone or to a match-up of one CU and one DU for sharing. Every DU should be matched with a CU to share a channel for transmission. If channel is shared by CU ( ) and DU ( ), and transmit signals to the BS and to DR ( ), respectively, at the same time. Here, to denote the channel gain of the links from the CUs, we use the letter ℎ. Furthermore, for the index of the BS, we use . Then, the channel gain of the link from to the BS on channel is denoted as ℎ , . We use the letter to denote the channel gain for the links from the DTs, e.g., , denotes the channel gain of the link from to on channel . Note that the signal from to causes co-channel interference (CCI) with the BS, and the corresponding channel gain is denoted as , . Similarly, ℎ , denotes the channel gain of the link from to . The channel state information (CSI) of the links can be acquired by utilizing either a sounding reference signal (SRS) or reference signal (RS) associated with THE D2D link [25,26]. In this paper, we assume that the CSI of all possible combinations of the links is available at the BS, which may require indispensable CSI feedback. However, this assumption has been widely used in previous literature for performance evaluation [7][8][9][11][12][13][14][15][16][17][18][19].  Figure 1b illustrates a channel assignment scheme, where each channel with a fixed channel bandwidth is pre-assigned to one CU before the channel assignment. In each channel, there exist boxes where their size implies the channel bandwidth, and the intensity of the box color represents the sum rate of the corresponding match-up of CU and DU. For instance, channel 1 is pre-assigned to , and the color intensity varies depending on which DU is matched to share channel 1. For channel 1, the match-up of and is the best one from among possible matches, and thus, the corresponding box linked to is painted with the darkest green color. Similarly, and are selected as the best matches for and , respectively. On the rightmost side in Figure  1b, the result of the channel assignment is shown.

System Model
Meanwhile, Figure 1c illustrates the 3-D channel assignment scheme proposed in this paper. The channel bandwidth, which corresponds to the size of each box, is determined adaptively based on the channel capacity of each match-up so the CU and DU can achieve the pre-determined target data rate. Note that if a channel is shared, depending on which CU and DU are matched, the CCI and, accordingly, the channel capacity can vary on a case-by-case basis. Therefore, on each channel, there exist possible matches with different bandwidths, which are illustrated by different boxes in Figure 1c. When the channel is assigned to minimize the sum bandwidth, the best match per channel is decided based on the channel bandwidth. For instance, matching and corresponds to the smallest box on channel 1 and can be the best match on channel 1. Similarly, and can be the best match on channel 2, and and can be the best match on channel . The rightmost side in Figure 1c shows the result of the channel assignment where the channel bandwidths are different. Here, we need to note that the DU chosen as the best match in an earlier selection phase can be the best match for a later selection phase as well. For example, can also be the best match for . So, to find the best channel assignment involves searching for a total of possible matches.

Problem Formulation
Now, let us calculate the minimum required channel bandwidth, , , for an arbitrary matchup of and on channel n to guarantee the target data rate. For simplicity, the target data rates of the CU and DU are denoted as and , respectively. transmits signals with power to the BS, and transmits signals with power to . Then, given the channel bandwidth for , , , SINR at the BS can be expressed as where is the noise power density. Then, the achievable data rate, , , needs to guarantee : Similarly, with the channel bandwidth for , , , the SINR at can be written as and the following condition needs to be satisfied: Then, to guarantee both and concurrently, , should be determined as follows: We assigned the channels to minimize the total channel bandwidth, while guaranteeing and for the CUs and DUs, respectively. To minimize the total channel bandwidth, a total of possible matches need to be investigated. In addition, the transmission power of the CUs and the DTs should be optimally determined. Here, we introduce 1, … , and 1, … , to denote the sets of CUs and DUs, respectively. Moreover, we introduce a channel matching indicator: , . When and share a channel , , is set at 1; it is 0 otherwise: For simplicity, we assume that , i.e., the number of CUs is the same as that of DUs. When , CUs have no match. Although the CUs that are not matched also creates a resource allocation problem, it is trivial and straightforward. Then, the proposed optimization problem that minimizes the overall bandwidth can be formulated as follows: : min , , , , , , where Equation (7a) indicates that every CU in must be matched with one DU on one channel, Equation (7b) implies that every DU must be matched with one CU on one channel, and Equation (7c) indicates every channel must be assigned at most once. Meanwhile, Equations (7d) and (7e) are rate guarantee constraints, Equation (7f) is to satisfy the binary condition, and Equations (7g) and (7h) are power limit constraints with maximum transmission power for the CU and DU, respectively. In cases of , only the equality sign in Equation (7a) is simply replaced with the less-than inequality sign. The key notations and symbols are listed in Table A1 in Appendix A.

Rate Guarantee Resource Allocation Scheme
The optimization problem is a mixed integer non-linear programming (MINLP) problem, which is usually intractable. To reduce the computational complexity, we decompose into two suboptimization problems: 1) finding optimal transmission power control to derive the minimum channel bandwidth for arbitrary matching, and 2) solving the 3-D channel assignment problem.

Minimum Channel Bandwidth for Arbitrary Matching, , *
First, the minimum channel bandwidth for arbitrary matching (CU i, DU j, channel n) is derived for the maximum transmit power limit, . Since the channel capacity varies depending on the received SINR, the channel bandwidth needs to be adjusted to guarantee the target rate, and , for CUs and DUs, respectively. The optimal power allocation problem to minimize the required channel bandwidth can be formulated as subject to Equations (2), (4), (7g), and (7h); and are the transmission power densities, given as , and , . ( Interestingly, Equation (8) is a quasi-convex function. In addition, all of the constraints of are convex. Therefore, is a quasi-convex optimization problem. The proof is provided in Appendix B. Note, since is a quasi-convex problem, this problem can be solved iteratively via, for example, a root-finding algorithm [27].
We first investigate the achievable SINR as a function of and . The channel bandwidth is inversely proportional to the SINR, as seen in Equations (2) It should be noted that Equation (10) is quasi-concave while Equation (8) is quasi-convex. Now, we visualize the optimization problem of Equation (10) in Figure 2. We assume ℎ , 100dB, , 108 dB, ℎ , 139 dB, , 115 dB, 23 dBm, and 174 dBm Hz ⁄ .
We also assume that 1Mbps. Note that the feasible region is determined by substituting Equation (9) into (2), (4), (7g) and (7h), which corresponds to the transmission power limit to guarantee the target data rate. In particular, if either or (or both) uses , , is located on the boundary, which is illustrated as a black bold line.
Let us consider three representative points in Figure 2. Table 1  . Then, the simultaneous improvement of , and , is no longer available. Therefore, point is the optimum point for power allocation, which minimizes , . In our example, it is that uses the full transmission power, i.e., the transmission power of , , 23 dBm, as seen in Table 1.
It is easy to see that the condition in Equation (11) holds, even for . Now the remaining problem is to increase the transmit power maximally until either the CU or the DT reaches its maximum transmit power while satisfying Equation (11). Since the channel bandwidth depends on the received SINR, which also depends on the channel bandwidth as well as the transmit power of the CU and DT, the optimum power can be obtained in an iterative fashion, e.g., via the root-finding algorithm.
In order to reduce the computational complexity, we provide a suboptimal solution by assuming that , ≫ 1 and , ≫ 1 at the optimum point. Then, Equation 11 can be approximated as Figure 2. Geometrical interpretation of optimal power densities.
It is easy to see that the CU transmits at when 1, i.e., ; then, can be obtained with (12), (11), and (4). When 1, equals instead. Note, if Newton's method is used to solve , a good approximation of the complexity is , where is the complexity of calculating inverse matrix of Jacobian matrix with -digit precision [28]. Note, if Gaussian elimination is applied, is approximated to . Meanwhile, because (12) is a closed-form solution, the corresponding complexity to obtain the suboptimal solution is 1 .

Three-Dimensional Channel Assignment
Once we obtain , * for all ∈ Ω , ∈ Ω , and channels, the channel assignment problem can be interpreted as two joint problems of 1) which CU each DU is matched with, and 2) which CU-DU match each channel is allocated to so that the overall bandwidth can be minimized. This can be formulated as an integer programming problem (IPP) as follows: , * ∈ ∈ (13) subject to Equations (7a), (7b), (7c), and (7f).
Note that is a 3-D channel assignment problem with respect to triple , , of , for all ∈ , ∈ , and 1, … , . The optimal solution can be found by the branch-and-bound algorithm. However, is a non-deterministic polynomial-time (NP)-hard problem with huge computational complexity [29].
We propose an algorithm for searching out the suboptimal solution to with much lower complexity. Here, to express the problem in canonical form, we introduce two column vectors , … , and , … , . The 1 1 element equals the , , th element of the corresponding 3-D matrix. Then, Equation (13) can be rewritten as * min (14) and , Equation (7c) can be rewritten as We relax Equation (14c) by introducing Lagrange multipliers, , … , , where the elements are all positive, i.e., 0, ∀ . Then, we can formulate the relaxed problem of , which is given by subject to Equations (14a), (14b), and (7f).
It is interesting to note that the matrix in Equation (14a) and the matrix in Equation (14b) are totally unimodular (TU). Then, the binary variables can be transformed into continuous variables that turn into a linear programming problem (LPP) for which the solution can be found with polynomial time complexity by a projective interior method [30]. The proof that is an LPP is provided in the Appendix C. Given , let be the solution to . Due to Lagrangian relaxation, the constraint that every channel cannot be used more than once is not strictly satisfied. Thus, from , we only extract the matching information, , , which is 1 if and are matched. Otherwise, , 0: , , Note that , carries no information about which channel is assigned for that match. Hence, the complete channel assignment can be decomposed as where is the channel assignment indicator, and 1 if channel is assigned to . Otherwise, 0. Given , , according to Equation (17), can be rewritten with respect to . To do that, let denote the channel bandwidth assigned to on channel n. Then, it can be obtained by , , Let and be column vectors of size , derived from and , respectively. Now, can be reformulated as the LPP to find the optimal channel for each match, as follows: : min (19) subject to where Equation (19a) is the constraint that only one channel is allocated to a single CU-DU match, and Equation (19b) is the constraint that a channel cannot be used more than once. Note that because the matrix in Equations (19a) and (19b) is totally unimodular (TU), the binary constraint corresponding to is relaxed, as in Equation (19c). Given , is a lower bound of * owing to positive Lagrangian multipliers, i.e., * . In contrast, is no smaller than * . In this paper, we use the subgradient algorithm to update . A detailed explanation of the subgradient algorithm can be found in [31]. Since the subgradient of at is , can be updated as where ∑ is the step size. The detailed procedures for the proposed channel assignment scheme are given in Algorithm 1. Finally, we obtain the complete channel assignment information, , * , according to Equation (17). 20. return * .

Performance Evaluation
The simulation parameters are shown in Table 2. Unless otherwise stated, we assume 1Mbps and the maximum D2D link distance, 50m. The number of CUs is the same as the number of DUs, i.e., . For simplicity, the receiver noise figure is 0dB by assuming an ideal receiver. We used Matlab 2019a to obtain the solutions to the IPP in and the LPPs in and . We investigated the cumulative distribution function (CDF) of the channel bandwidth in Figure  3. Before channel assignment, the minimum channel bandwidth is collected for all the arbitrary matchings by solving . Through the channel assignment process, channels are assigned to better CU-DU matches in order to minimize the total channel bandwidth, and hence, CU-DU matches requiring a wider channel bandwidth are sorted out from the candidates. This results in steep CDF after channel assignment. Before channel assignment, the average channel bandwidth is approximately 155 kHz, regardless of . With 4, the average channel bandwidth is reduced to approximately 101 kHz after channel assignment, and is further reduced by 16% when 10. In terms of the median SINR, compared to before channel assignment, it increased by 9dB when 4 and by 15 dB when 10, which results in the reduction in channel bandwidth after channel assignment, as seen in Figure 3. We need to note that this increase in the SINR after channel assignment validates the assumptions that , ≫ 1 and , ≫ 1 used in Equation (12) Figure 4 shows a snapshot of the iterative operation of Algorithm 1 to obtain the minimum total channel bandwidth. Here, we assume 6. For comparative purposes, the optimal value, * , is given by solving the original IPP problem, . As Algorithm 1 iterates, the Lagrangian multiplier is updated so the output of increases closer to * . Since Lagrangian relaxation is applied, multiple CU-DU pairs are allowed to be matched to a channel when solving , which is resolved by solving , where a channel is assigned solely by a single CU-DU match. Since channels are only re-assigned when the CU-DU matching information obtained from solving is kept intact, the output of may not decrease, and will remain at the same value. This is because although the Lagrangian multiplier is updated to bring an increase in , the CU-DU matching information may be kept unchanged. At the seventh iteration, the algorithm terminates, since the gap between two outputs becomes less than a pre-determined threshold.
(c)  Figure 5 shows the average channel bandwidths of different channel assignment schemes when varying the number of CUs. For Optimal, the transmission power is determined optimally, and the channels are assigned by solving . Suboptimal adopts the suboptimal transmission power scheme and Algorithm 1 for the channel assignment. Note that the performance of Suboptimal matches Optimal quite well. For comparison, we implemented simple channel assignment schemes: random selection and greedy. In random selection, CU-DU-channel matchings are determined randomly. If random selection is applied, the average channel bandwidth slightly increases as increases. When the greedy scheme is used, a CU selects the best DU-channel pair in consecutive order. Once the first CU selects a DU-channel pair from among choices, the next CU repeats the same selection from among 1 choices, and so on. Therefore, each CU does not calculate the possibility that the current best selection may not be best in terms of total channel bandwidth, i.e., it is beneficial to select the second best one for now, since the best one is more beneficial for the next selection in terms of total channel bandwidth. The greedy scheme is computationally beneficial at the cost of performance degradation. When 10, the average channel bandwidth of the proposed scheme can be reduced by approximately 25% over the greedy scheme.  Figure 6 compares two cases, depending on whether each CU is pre-assigned a channel before channel assignment, as seen in Figures 1b and c. As the number of CUs, , increases, a CU-DU pair can have a better channel with a higher probability. For instance, when the maximum D2D distance is 50 m, the performance gain increases from approximately 9.5% to 17.5% as increases from two to 10 without channel pre-assignment. When each CU is pre-assigned a channel before channel assignment (denoted as pre-assignment), the remaining channel assignment problem is to find the best matches for CUs and DUs to minimize the total channel bandwidth while guaranteeing the target data rate. Then, the problem reduces to two-dimensional channel assignment. Due to more candidate matches being available with three-dimensional assignment, the proposed 3-D assignment algorithm outperforms the pre-assigned case, and this performance gain is higher with a larger and more CUs. For example, with = 50 m and 10 CUs, a 21% wider channel bandwidth is required for the pre-assigned case than the proposed scheme. We also simulated the proposed scheme under Rician fading for 20 m. In the simulations, we assumed the K factor was 9dB, and used the same path loss model applied in [33]. Under Rician fading, channel pre-assignment is shown to be still disadvantageous in terms of the average channel bandwidth, causing approximately 15% more bandwidth consumption. Furthermore, because DUs experience comparatively better channel quality under Rician fading, DTs tend to use less transmission power to achieve their target data rates. It reduces the interference to the BS, yielding better channel quality for the links from CUs to the BS. In our simulations, the average channel bandwidth decreases approximately by 10~13% if Rician fading is adopted instead of Rayleigh fading.
(c) Figure 6. Average channel bandwidth versus the number of CUs with and without channel preassignment. Figure 7 investigates the CDFs of the achieved data rate of the proposed rate guarantee scheme and those of the two previous schemes with different objectives; maximization of spectral efficiency (Max. SE) [11,15] and maximization of energy efficiency (Max. EE) [17,19]. For each previous scheme, only one channel is assigned to a DU in Single case, whereas multiple channels can be assigned to a DU and multiple DUs can be assigned to a channel in a Multiple case. In the previous schemes, the bandwidth of every channel is set and fixed at the same value as the average channel bandwidth of the proposed rate guarantee scheme, and the SINR threshold is appropriately determined to meet the minimum data rate with the given channel bandwidth. Figure 7a shows the CDFs of the achieved data rate of CUs. By adaptively adjusting the channel bandwidth, the proposed rate guarantee scheme can guarantee the target data rate. In contrast, approximately 10% CUs cannot meet their target data rate in the previous schemes. Figure 7b shows the CDFs of the achieved data rate of DUs. As expected, DUs are also able to achieve their pre-determined target data rates in the proposed scheme. In the previous schemes, when the target data rate of both UEs cannot be simultaneously guaranteed, DUs are forced to decrease their transmission power. This brings less interference to the CUs sharing the same channel. In addition, the CUs can achieve the higher achievable data rate at the sacrifice of the data rate of DUs. For instance, when DUs are assigned only single channel, approximately 70% of DUs cannot meet the target data rate. Obviously, if DUs are allowed to be assigned more than one channel, this probability is reduced. Nevertheless, owing to the fixed bandwidth limitation, about 30% of DUs are still not able to achieve the pre-determined target data rate.
The computational complexity of the proposed channel assignment algorithm is compared with that of exhaustive search [34] and the iterative Hungarian method [35] in Table 3. For the evaluation of the complexity, the total number of operations is evaluated by Big notation. The problem is a 3-dimensional (3-D) assignment problem, and the complexity is ! if the exhaustive search is used. The iterative Hungarian method decomposes the 3-D problem into two integer programming problems (IPPs). Then, the Hungarian method is iteratively applied to obtain the solutions of the IPPs until they converge. In this scheme, the complexity per iteration is . In the proposed scheme, by exploiting the properties of the totally unimodular (TU) matrix, the 3-D problem is decomposed into two linear programming problems (LPPs) which can be solved with much less complexity than the IPPs. If the primal-dual interior point method is used, the upper bound (U.B.) and the lower bound (L.B.) of the complexity are log and log , respectively [36]. Here, is the accuracy parameter used to relax the complementary slackness condition, and we use 0.01 [37]. As increases, the complexity of the proposed scheme decreases significantly compared to the previous schemes. For instance, when 4, the complexity can be reduced by about 50% compared to the iterative Hungarian method. However, the complexity reduction reaches up to approximately 93% when 10.  We proposed a resource allocation scheme that minimizes the total channel bandwidth while guaranteeing the target data rate of the individual UE. The channel bandwidth and transmit power were jointly optimized together with the 3-D channel assignment problem using triples ( CUs, DUs, and channels). The 3-D channel assignment problem involves finding the best match among possible matches. This large number of candidates enables us to provide improved performance over two-dimensional (2-D) channel assignment where each channel is pre-assigned. However, the proposed resource allocation problem becomes an MINLP. An optimal and iterative power control scheme and a suboptimal method were suggested to derive the minimum channel bandwidth for an arbitrary match. By exploiting the property of a TU matrix and with Lagrangian relaxation, the 3-D IPP channel assignment problem could be solved by iteratively solving two 2-D LPP problems based on a subgradient algorithm. The proposed scheme outperforms the existing schemes, including the random selection and greedy schemes, in terms of the average channel bandwidth. Moreover, it guarantees every UE achieves the target data rate, whereas the existing schemes cannot. In addition, we investigated the impact of the maximum D2D distance and the number of CUs.  Required bandwidth for if and are matched on channel is also convex. It is easily proved that the intersection of any two convex sets is also a convex set. Therefore, the feasible region of is a convex set.