Next Article in Journal
Privacy-Preserving Distributed Analytics in Fog-Enabled IoT Systems
Next Article in Special Issue
Modeling and Performance Analysis of Opportunistic Link Selection for UAV Communication
Previous Article in Journal
The Danger of Walking with Socks: Evidence from Kinematic Analysis in People with Progressive Multiple Sclerosis
Previous Article in Special Issue
Interference Spreading through Random Subcarrier Allocation Technique and Its Error Rate Performance in Cognitive Radio Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Energy-Efficient Unselfish Spectrum Leasing Scheme for Cognitive Radio Networks

by
Denis Bilibashi
,
Enrico M. Vitucci
,
Vittorio Degli-Esposti
and
Andrea Giorgetti
*
Department of Electrical, Electronic, and Information Engineering “Guglielmo Marconi” (DEI), CNIT, University of Bologna, 40126 Bologna, Italy
*
Author to whom correspondence should be addressed.
Submission received: 1 October 2020 / Revised: 23 October 2020 / Accepted: 26 October 2020 / Published: 29 October 2020
(This article belongs to the Special Issue Advances in Cognitive Radio Networks)

Abstract

:
Cooperative Communications in Cognitive Radio (CR) have been introduced as an essential and efficient technique to improve the transmission performance of primary users and offer transmission opportunities for secondary users. In a typical multiuser Cooperative Communication in CR, each primary user can choose one secondary user as a relay node. To encourage the cooperative behavior of the secondary users, primary users lease a fraction of their allocated spectrum to the relay secondary users to transmit their data packets. In this work, a novel unselfish spectrum leasing scheme in CR networks is proposed that offers an energy-efficient solution minimizing the environmental impact of our network. A network management architecture is introduced, and resource allocation is proposed as a constrained sum energy efficiency maximization problem. The optimization problem is formulated and solved using non-linear programming methods and based on a modified Kuhn-Munkres bipartite matching algorithm. System simulations demonstrate an increment in the energy efficiency of the primary users’ network compared with previously proposed algorithms.

1. Introduction

The exponential growth of wireless services in the last decades has called for additional spectrum to accommodate the ever-increasing traffic flow [1]. Besides, the current static spectrum policy, used by most countries in the world, has almost resulted in licensed spectrum exhaustion of some frequency bands [2]. On the contrary, recent spectrum utilization measurements have shown that the available spectrum opportunities are severely underutilized [3]. Cognitive radio (CR), with its ability to sense available frequency bands and adaptively adjust transmission frequency, has attracted considerable attention as one of the most promising solutions to spectrum underutilization [4,5,6,7]. Traditionally in CR networks (CRNs), there two categories of users, namely, primary users (PUs) and secondary users (SUs). PUs are licensed users, and they have the exclusive right to use their respective channels, while SUs are unlicensed, and they wish to use the underutilized channels [8]. To pay off the use of a licensed band, cooperative CRNs (CCRNs) have been proposed in which SUs act as relays to provide a better transmission performance to the PUs [9,10,11].
Extensive research has been conducted regarding spectrum leasing in CCRN. To address the challenges several approaches have been adopted. The game-theoretical approaches, such as Stackelberg game [12], are used to achieve the equilibrium state (e.g., Nash Equilibrium [13,14]) and involves PUs and SUs as players of the game. In Reference [15], a scenario, in which PUs compete for cooperation with available SUs under shared constraint set by SU’s requirements, is modeled as a generalized Nash Equilibrium. The authors of Reference [16], propose a noncooperative game to find the optimal power of SUs consumed in a scenario where SUs are used as relays in cooperation to aid PUs’ transmission. In Reference [17], a spectrum leasing technique to improve the spectral utility of the PUs has been proposed. Nash bargaining is used to determine the optimal division of temporal resources between relaying and reimbursement. In Reference [18], a novel approach is presented in which the spectrum leasing problem is formulated as a Stackelberg game, while the relay selection is solved by the Hungarian method. The authors of Reference [19], derive a cooperative game to find the optimal spectrum leasing time and allocated power of SUs.
Another popular approach is convex optimization. An adaptive spectrum leasing scheme is proposed in Reference [20] to improve the overall network throughput. The authors of Reference [21], present a spectrum leasing scheme for an overlay CRN in which PUs’ transmission is assisted by a cognitive multi-hop multi-branch network. The branch selection and spectrum allocation are jointly optimized such that the sum-rate of the SUs is maximized. A multi-objective optimization framework to incorporate average sum power minimization and the leased time minimization is proposed in Reference [22]. In Reference [23], a cooperative sharing model for cognitive wireless powered networks is proposed, to accomplish PU’s transmission early and vacate opportunity for SU’s transmission. An energy efficiency and spectral efficiency trade-off in an energy harvesting cooperative CRN is addressed in Reference [24]. The objective is to maximize energy and spectral efficiency via joint optimization of spectrum sensing time and self as well as cooperative transmission gains.
Different other approaches have been chosen to address the challenges of spectrum leasing in CCRN. In Reference [25], an opportunistic relaying based spectrum leasing is presented in which PUs lease their frequency band to the SUs and, at the same time, the primary network benefits in terms of outage probability. A cooperative relaying scheme is investigated in Reference [26], where one of the inactive SUs is used as a relay for the PUs. By adopting an aspiration level to model SUs behavior, the SUs’ strategy dynamics is formulated as a discrete-time Markov chain [27].
However, the majority of the research works regarding CRNs, are considered mostly PU-centered, with PUs that decide if and when to exploit SUs as relays based on certain conditions that can improve PUs’ transmission. This approach cannot guarantee full cooperation between PUs and SUs because will be dominated by the selfish behavior of PUs. In addition, networks are required nowadays to be more efficient not only in terms of spectrum use but also in terms of optimal power while guaranteeing a satisfying level of quality of service (QoS). In References [15,25,26], a scenario where a single pair of PUs or SUs is considered. Single user transmissions in CCRN compared to multiuser ones, can result incomplete because a full description of the efficiency of a cell cannot be given. It is proven that the rising number of users in a CCRN can significantly affect the overall performance of the network. Furthermore, spectrum leasing schemes in which PUs or SUs are assumed to have fixed transmission power like in Reference [18], cannot fully exploit the possibility of a more energy efficient network. Besides, energy-efficient systems are required not only from an engineering point of view but also to decrease the environmental impact through networks with a low carbon footprint [28].
Motivated by the aforementioned observations, we aim at jointly minimize the overall power consumption for energy-efficient communication in multiuser CRN. In addition, to provide improved cooperation between PUs and SUs and to secure an unselfish behavior of each user, we propose a resource management network architecture for a considered cell. This architecture, composed of various entities, will be responsible for making decisions regarding the usage of relay and for applying the optimal strategy to achieve the maximum energy efficiency of the cell. Our proposed CRN consists of multiple PUs and multiple SUs in which SUs are used as relays to enhance the PUs’ transmission performance and, in return, PUs lease a fraction of their bandwidth to the corresponding SUs. We formulate an energy efficiency based resource allocation, power allocation, and relay selection as an optimization problem. Our analysis shows that this problem can be divided into two subproblems: (1) power allocation, (2) transmission mode and relay selection. To solve the power allocation subproblem, we choose the maximization of the sum of energy efficiencies as the optimization target. This method was initially introduced in Reference [29] but due to its sum-of-ratio form can become difficult to solve. The method is further developed in our work by transforming power allocation into a subtractive form optimization problem with methods similar to those used to solve classical fractional programming problems [30,31,32]. Therefore, the optimal powers closed-form expressions can be derived through a two-layer optimization. Furthermore, the transmission mode and relay selection subproblem is solved based on a modified Kuhn-Munkres bipartite matching algorithm. System simulations performed in a referenced multiuser case validate the effectiveness of our approach proven by the proposed algorithm’s convergence and the increment in the efficiency of the PUs’ network performance when channel conditions get worse.
To summarize, the main contributions of the paper are:
  • A novel unselfish resource allocation scheme for CRNs. Unlike most of works where PUs decide whether or not to cooperate with SUs, the decisions are made by a centralized approach with the goal of improving the overall energy efficiency of the cell.
  • A resource management network architecture responsible for maximizing the energy efficiency of the cell.
  • We propose a two-stage three-dimensional matching algorithm to maximize the energy efficiency of the cell, in which the joint optimization problem is decoupled into two subproblems and solved separately in two stages.
The rest of this paper is organized as follows. In Section 2, the architecture and the system model for CCRN are introduced. In Section 3, we formulate the energy efficiency based resource allocation, power allocation, and relay selection as an optimization problem. The non-convex energy efficiency for CCRN optimization problem is formulated and solved via an iterative algorithm in Section 4. Section 5 presents numerical performance results, and finally, conclusions are provided in Section 5.

2. Proposed Architecture and System Model for CCRN

In CCRN, one of the goals is to use the available spectrum resources in an efficient and coordinated way to guarantee a satisfying level of QoS for all users and achieve performance enhancement of the whole network. To this end, in this section, we propose a resource management network architecture consisting of a number of user resource management entities (URMEs), local resource management entities (LRMEs), and one cell resource management entity (CRME), as shown in Figure 1. The main functions of UMREs, LRMEs, and CRME are described as follows.
URME: It is a functional module embedded in each PU and SU, used to store channel state information, device characteristics and service requirements, and so forth. Through contacting associated LRMEs, URMEs send the collected information to the LRMEs and receive the resource allocation strategy accordingly also from CRMEs through LRMEs.
LRME: It is deployed in each primary base station (PBS) or secondary base station (SBS), being responsible for managing local resource status through interacting with the associated URMEs and the CRME. More specifically, receiving the network and service information from URMEs and then forwarding to the CRME, and receiving the resource management strategy from the CRME and forwarding to the URMEs.
CRME: By interacting with the LRME, the CMRE receives the network status, channel state information, and user service requirement information of all the users within the cell; then performs the proposed resource allocation and relay selection algorithm to obtain the optimal strategy for the cell and communicated to the associated LRMEs.

System Model

In this paper, we consider a CCRN consisting of multiple PUs, multiple SUs, one PBS, and one SBS. Assume PUs are allowed to transmit to the PBS simultaneously, by orthogonal frequency division multiple access (OFDMA). Further, assume that PUs may transmit to the PBS in direct transmission mode or one-hop relay transmission mode, while the relay SUs may also transmit their own data packets to the SBS. Figure 1 illustrates the scenario considered in this paper.
To encourage SUs to relay data packets for the PUs, we propose a spectrum leasing scheme where the PUs lease part of their allocated spectrum to relay SUs so that they can transmit their own data to the SBS exploiting the licensed spectrum [33,34,35]. In this paper, we assume that the decode-and-forward (DF) scheme is applied at each relay node [36].
Let M and K denote the number of PUs and the number of SUs, respectively, B m = ρ m B denote the allocated bandwidth of the mth PU, m = 1 , , M , where ρ m [ 0 , 1 ] denotes the spectrum fraction of the mth PU for transmitting the data packets of the PU in the relay transmission mode, and as a consequence B k = ( 1 ρ m ) B , k , is the remaining spectrum available for the relay SU to transmit its own data to the SBS. To support relay communication, the transmission time slot T is divided into two periods. For the first part, the PUs transmit their data packets to the corresponding relay SUs, then the SUs forward the received data packets to the PBS during the remaining time. Meanwhile, the SUs transmit their own data packets to the SBS for the whole time period. Figure 2 shows the time and spectrum division mode for relay transmission of the mth PU.
We assume, all channels undergo flat Rayleigh fading and log-normal shadowing and we consider that LRMEs know all channel state information. The distance between mth PU and kth SU relay is denoted as d m , k ( p , s ) . Then the channel power gain is represented as h m , k ( p , s ) = k 0 ( d m , k ( p , s ) ) γ g m , k ( p , s ) s m , k ( p , s ) , where k 0 denotes the average channel gain at a reference distance d 0 = 1 m and this gain depends on carrier frequency, transmitter and receiver antenna gains and propagation characteristics, γ > 0 is the path-loss exponent, g m , k ( p , s ) is a Rayleigh distributed r.v., and s m , k ( p , s ) is a log-normal r.v. with shadowing parameter σ s . The same channel model is considered for the other links h m ( p , d ) , h m , k ( p , r ) , h m , k ( s ) , representing respectively the channel gains of the links between the mth PU and PBS in direct transmission mode, the kth SU and PBS in relay mode, and the kth SU when transmitting its own data to the SBS on the fraction of bandwidth leased by the mth PU.
In the case of multiple relay SUs being available, the optimal relay selection scheme should be designed. For each PU-SU pair, the transmit power of the PU, P m , k ( p , s ) and the transmit power of the relay SU, consumed during relay P m , k ( p , r ) and during its own data transmission P m , k ( s ) , should be designed for performance optimization of the network.

3. Proposed Joint Resource Allocation and Relay Selection Scheme

In the following subsections, we formulate an optimization problem, and we present a new solution to the energy efficiency based resource allocation, power allocation, and relay selection.

3.1. Energy Efficiency of PUs

The energy efficiency of all the PUs can be expressed as
η ( p ) = m = 1 M η m ( p ) ,
where η m ( p ) denotes the energy efficiency of the mth PU. As PUs may choose direct transmission or cooperative transmission through a relay SU for information transmission to the PBS, the energy efficiency of the mth PU can be calculated as
η m ( p ) = β m ( d ) η m ( p , d ) + k = 1 K β m , k ( c ) η m , k ( p , c ) ,
where η m ( p , d ) is the energy efficiency of the mth PU when using direct mode, while β m ( d ) { 0 , 1 } denotes the transmission indicator of the mth PU for the direct mode. More specifically, if β m ( d ) = 1 , the mth PU uses direct transmission mode, while for β m ( d ) = 0 , the PU exploits SU relay assisted transmission. The term η m , k ( p , c ) denotes the energy efficiency of the mth PU when using the kth SU as relay node for cooperative transmission mode, while β m , k ( c ) { 0 , 1 } is the corresponding transmission and relay selection variable, that is, β m , k ( c ) = 1 indicates that the mth PU uses the kth SU as relay node for cooperative transmission, and β m , k ( c ) = 0 , otherwise. Note that β m ( d ) and β m , k ( c ) are mutually exclusive. In the following subsections, we derive the expression of η m ( p , d ) and η m , k ( p , c ) , respectively.

3.1.1. Direct Mode

The energy efficiency of the mth PU in direct mode, that is, η m ( p , d ) can be calculated as
η m ( p , d ) = R m ( p , d ) P m ( p , d ) + P c ( p ) ,
where P m ( p , d ) is the power consumed by the mth PU when transmitting in direct mode, P c ( p ) denotes the circuit power consumption of the mth PU, which is assumed to be a constant for all the PUs in this paper, and R m ( p , d ) denotes the data rate achieved by the mth PU in direct mode, which can be expressed as
R m ( p , d ) = B log 2 1 + P m ( p , d ) h m ( p , d ) σ 2 ,
where σ 2 is the noise power of the link between the mth PU and the PBS. Without loss of generality, the noise power of all the transmission links are assumed to be the same in this paper.

3.1.2. Relay Mode

The energy efficiency of the mth PU when using the kth SU as relay node for cooperative transmission, denoted by η m , k ( p , c ) , can be expressed as
η m , k ( p , c ) = R m , k ( p , c ) P m , k ( p , c ) .
The denominator P m , k ( p , c ) is the power consumption of the mth PU when using the kth SU as a relay node and can be calculated as
P m , k ( p , c ) = P m , k ( p , s ) + P m , k ( p , r ) + P c ( p ) + P c ( s ) ,
where P c ( s ) denotes the circuit power consumption of the SU which is assumed to be a constant for all the SUs. The numerator R m , k ( p , c ) represents the data rate of the mth PU when using the kth SU as relay node for cooperative transmission. R m , k ( p , c ) can be expressed as
R m , k ( p , c ) = min ( R m , k ( p , s ) , R m , k ( p , r ) ) ,
where R m , k ( p , s ) and R m , k ( p , r ) are, respectively, the data rate of the link from the mth PU to the kth SU
R m , k ( p , s ) = t 1 ρ m B log 2 1 + P m , k ( p , s ) h m , k ( p , s ) σ 2
and that from the kth SU to the PBS when the kth SU is chosen as the relay node of the mth PU
R m , k ( p , r ) = ( 1 t 1 ) ρ m B log 2 1 + P m , k ( p , r ) h m , k ( p , r ) σ 2 .

3.2. Energy Efficiency of SUs

The energy efficiency of all the SUs in the CCRN, denoted by η ( s ) , can be calculated as
η ( s ) = m = 1 M k = 1 K β m , k ( c ) η m , k ( s ) ,
where η m , k ( s ) is the energy efficiency of the kth SU when transmitting its own data on the subchannel leased by the mth PU, which can be calculated as
η m , k ( s ) = R m , k ( s ) P m , k ( s ) + P c ( s ) ,
where R m , k ( s ) denote the data rate of the kth SU when transmitting its own data on the subchannel leased by the mth PU
R m , k ( s ) = ( 1 ρ m ) B log 2 1 + P m , k ( s ) h m , k ( s ) σ 2 .

4. Energy Efficiency for CCRN Problem Formulation and Solution

The total energy efficiency of the PUs and the SUs can be formulated as:
η = η ( p ) + η ( s ) .
In the following, we propose a joint resource allocation, power allocation, and relay selection strategy to maximize the total energy efficiency (13).

4.1. Optimization Problem Formulation

The energy efficiency based power allocation, transmission mode and relay selection scheme can be obtained by solving the following optimization problem:
max [ β m ( d ) , β m , k ( c ) , P m ( p , d ) , P m , k ( p , s ) , P m , k ( p , r ) , P m , k ( s ) ] η s . t . C 1 : β m ( d ) + k = 1 K β m , k ( c ) 1 , m = 1 , 2 , . . . , M C 2 : m = 1 M β m , k ( c ) 1 , C 3 : P m ( p , d ) P m ( p , max ) , if β m ( d ) = 1 , C 4 : P m , k ( p , s ) P m ( p , max ) , if β m , k ( c ) = 1 , C 5 : P m , k ( p , r ) + P m , k ( s ) P k ( s , max ) , if β m , k ( c ) = 1 , C 6 : R m ( p , d ) R m ( p , min ) , if β m ( d ) = 1 , C 7 : R m , k ( p , s ) R m ( p , min ) , if β m , k ( c ) = 1 , C 8 : R m , k ( p , r ) R m ( p , min ) , if β m , k ( c ) = 1 , C 9 : R m , k ( s ) R m , k ( s , min ) , if β m , k ( c ) = 1 .
C1 represents the transmission mode indicator constraint since it is assumed that every PU can only choose direct transmission mode or relay transmission mode. Similarly, C2 specifies the relay mode indicator constraint as it is assumed that every PU can choose only one SU as its relay node and each SU can only forward packets for one PU. C3, C4 and C5 denote the maximum power constraint where P m ( p , max ) and P k ( s , max ) denote respectively the maximum transmit power of the mth PU and the kth SU. C6-C9 represent the minimum data rate constraints where R m ( p , min ) and R k ( s , min ) denote respectively the minimum required data rate of the mth PU and the kth SU.
The optimization problem in (14) can be classified as a nonlinear binary fractional program and therefore difficult to solve using traditional optimization methods [31]. From the optimization constraints C3–C9 given in (14), it can be shown that such an optimization problem can be transformed equivalently into two subproblems: (1) power allocation, (2) transmission mode and relay selection. The solution for each subproblem is presented in the following subsections.

4.2. Power Allocation Subproblem

Assuming that the mth PU selects direct transmission mode, that is, β m ( d ) = 1 , the energy efficient optimal power allocation problem can be formulated as:
max [ P m ( p , d ) ] R m ( p , d ) P m ( p , d ) + P c ( p ) s . t . C 1 : P m ( p , d ) P m ( p , max ) , C 2 : R m ( p , d ) R m ( p , min ) .
For the relay transmission mode, when the mth PU selects the kth SU as its relay node, that is, β m , k ( c ) = 1 , the energy efficiency of the PU-SU pair can be expressed as
η m , k ( p , c ) + η m , k ( s ) .
Therefore, the optimal power allocation problem for the mth PU and the kth SU can be formulated as:
max [ P m , k ( p , s ) , P m , k ( p , r ) , P m , k ( s ) ] R m , k ( p , c ) P m , k ( p , c ) + R m , k ( s ) P m , k ( s ) + P c ( s ) s . t . C 1 : P m , k ( p , s ) P m ( p , max ) , C 2 : P m , k ( p , r ) + P m , k ( s ) P k ( s , max ) , C 3 : R m , k ( p , s ) R m ( p , min ) , C 4 : R m , k ( p , r ) R m ( p , min ) , C 5 : R m , k ( s ) R m , k ( s , min ) .

4.2.1. Equivalent Problem Transformation

In order to solve the optimization problem in (17), we exploit the fractional programming approach [37]. As we mentioned earlier, we are considering DF at the relay node. In case R m , k ( p , s ) < R m , k ( p , r ) , our optimization problem becomes:
max [ P m , k ( p , s ) , P m , k ( p , r ) , P m , k ( s ) ] R m , k ( p , s ) P m , k ( p , c ) + R m , k ( s ) P m , k ( s ) + P c ( s ) s . t . C 1 : P m , k ( p , s ) P m ( p , max ) , C 2 : P m , k ( p , r ) + P m , k ( s ) P k ( s , max ) , C 3 : R m , k ( p , s ) < R m , k ( p , r ) , C 4 : R m , k ( p , s ) R m ( p , min ) , C 5 : R m , k ( p , r ) R m ( p , min ) , C 6 : R m , k ( s ) R m , k ( s , min ) .
To proceed, we rewrite the problem in (18) into an equivalent form [38]:
max [ P m , k ( p , s ) , P m , k ( p , r ) , P m , k ( s ) , α 1 , α 2 ] α 1 + α 2 s . t . C 1 C 6 in ( 18 ) , C 7 : R m , k ( p , s ) P m , k ( p , c ) α 1 , C 8 : R m , k ( s ) P m , k ( s ) + P c ( s ) α 2 ,
where C7 and C8 denote the energy efficiency constraints.
Theorem 1.
If the set ( P m , k ( p , s ) * , P m , k ( p , r ) * , P m , k ( s ) * , α 1 * , α 2 * ) is the solution for (19), then there exist ν and κ such that ( P m , k ( p , s ) * , P m , k ( p , r ) * , P m , k ( s ) * ) satisfies that Karush-Kuhn-Tucker (KKT) conditions for ν = ν ¯ , κ = κ ¯ , α 1 = α 1 * , α 2 = α 2 * of the following problem [39]
max [ P m , k ( p , s ) , P m , k ( p , r ) , P m , k ( s ) ] ϕ s . t . C 1 C 6 in ( 18 ) ,
where ϕ = ν ( R m , k ( p , s ) α 1 P m , k ( p , c ) ) + κ ( R m , k ( s ) α 2 ( P m , k ( s ) + P c ( s ) ) ) . Also the set of solutions satisfies a system of equations for ν = ν ¯ , κ = κ ¯ , α 1 = α 1 * , and α 2 = α 2 * , with
ν ¯ = 1 P m , k ( p , c ) κ ¯ = 1 P m , k ( s ) + P c ( s ) α 1 * = R m , k ( p , s ) P m , k ( p , c ) α 2 * = R m , k ( s ) P m , k ( s ) + P c ( s ) .
Proof for Theorem 1 is provided in Appendix A.
In Theorem 1 we proved that by satisfying the system of equations in (21) among the solutions of (20), the optimization problem in (18) can be solved. It shows as well that for an optimization problem with an objective function in fractional form, there exists an equivalent in subtractive form, that is, ϕ in our case. As a result, we use the equivalent objective function for the rest of the paper.

4.2.2. Energy Efficiency Maximization

In this section we present an iterative algorithm in order to find solutions for the problem in (20). After transforming our problem (18) into (20), it can be easily proven that the tranformed problem is convex and can be solved using the Lagrange dual method [40,41]. Thus, by relaxing the constraints the Lagrange function can be represented as
L ( P m , k ( p , r ) , P m , k ( p , s ) , P m , k ( s ) , λ , δ , ϵ , ξ , θ , μ , ν , κ ) = + ν ( R m , k ( p , s ) α 1 ( P m , k ( p , c ) ) ) + κ ( R m , k ( s ) α 2 ( P m , k ( s ) + P c ( s ) ) ) + λ ( P m ( ( p , max ) P m , k ( p , s ) ) + δ ( P k ( s , max ) P m , k ( p , r ) P m , k ( s ) ) + ϵ ( R m , k ( p , r ) R m , k ( p , s ) ) + ξ ( R m , k ( p , s ) R m ( p , min ) ) + δ ( R m , k ( p , r ) R m ( p , min ) ) + μ ( R m , k ( s ) R k ( s , min ) ) ,
where λ , δ > 0 and ϵ , ξ , θ , μ > 0 represent the Lagrange multipliers associated with power constraints and minimum data rate requirements, respectively. The corresponding Lagrange dual method can be expressed as follows:
min [ λ , δ , ϵ , ξ , θ , μ ] max [ P m , k ( p , r ) , P m , k ( p , s ) , P m , k ( s ) ] L s . t . λ 0 , δ 0 , ϵ 0 , ξ 0 , θ 0 , μ 0 .
The above dual problem can be solved by optimizing the transmit power for a fixed set of Lagrange multipliers, and updating the Lagrange multipliers iteratively as well as our parameters ν , κ , α 1 , α 2 . For a given set of Lagrange multipliers λ , δ , ϵ , ξ , θ , μ , the locally optimal transmit powers can be obtained through calculating the derivative of L ( P m , k ( p , r ) , P m , k ( p , s ) , P m , k ( s ) , λ , δ , ϵ , ξ , θ , μ ) over the transmit powers and setting it to zero thus obtaining
P m , k ( p , s ) * = t 1 ρ m B ( ν + ξ ϵ ) ( λ + ν α 1 ) l n 2 σ 2 h m , k ( p , s ) +
P m , k ( p , r ) * = ( 1 t 1 ) ρ m B ( ϵ + θ ) ( δ + ν α 1 ) l n 2 σ 2 h m , k ( p , r ) +
P m , k ( s ) * = ( 1 ρ m ) B ( μ + κ ) ( δ + κ α 2 ) l n 2 σ 2 h m , k ( s ) + ,
where [ z ] + = max ( 0 , z ) .
In case R m , k ( p , s ) R m , k ( p , r ) , a similar approach can be followed, leading to the optimization problem:
max [ P m , k ( p , s ) , P m , k ( p , r ) , P m , k ( s ) ] R m , k ( p , s ) P m , k ( p , c ) + R m , k ( s ) P m , k ( s ) + P c ( s ) s . t . C 1 : P m , k ( p , s ) P m ( p , max ) , C 2 : P m , k ( p , r ) + P m , k ( s ) P k ( s , max ) , C 3 : R m , k ( p , s ) R m , k ( p , r ) , C 4 : R m , k ( p , s ) R m ( p , min ) , C 5 : R m , k ( p , r ) R m ( p , min ) , C 6 : R m , k ( s ) R m , k ( s , min ) .
By applying to (28) the same transformation explained in Theorem 1 and applying Lagrange dual method to the transformed problem, we are able to calculate the transmit powers:
P m , k ( p , s ) * = t 1 ρ m B ( ξ + ϵ ) ( λ + ν α 1 ) l n 2 σ 2 h m , k ( p , s ) +
P m , k ( p , r ) * = ( 1 t 1 ) ρ m B ( ν ϵ + θ ) ( δ + ν α 1 ) l n 2 σ 2 h m , k ( p , r ) +
P m , k ( s ) * = ( 1 ρ m ) B ( μ + κ ) ( δ + κ α 2 ) l n 2 σ 2 h m , k ( s ) + .
Based on these solutions, we propose a two-layer optimization algorithm. In the inner layer the Lagrange multipliers in (24)–(26) and (28)–(30) can be updated through using the subgradient method [42], that is,
λ ( n + 1 ) = [ λ ( n ) ψ ( P m ( p , max ) P m , k ( p , s ) ) ] +
δ ( n + 1 ) = [ δ ( n ) ψ ( P k ( s , max ) P m , k ( p , r ) P m , k ( s ) ) ] +
ϵ ( n + 1 ) = [ ϵ ( n ) ψ | R m , k ( p , r ) R m , k ( p , s ) | ] +
ξ ( n + 1 ) = [ ξ ( n ) ψ ( R m , k ( p , s ) R m ( p , min ) ) ] +
θ ( n + 1 ) = [ θ ( n ) ψ ( R m , k ( p , r ) R m ( p , min ) ) ] +
μ ( n + 1 ) = [ μ ( n ) ψ ( R m , k ( s ) R k ( s , min ) ) ] + ,
where ψ denotes the learning rate.
In the outer layer the parameters ν , κ , α 1 , α 2 can be updated using Newton-Raphson method [43], that is,
ν ( n + 1 ) = ν ( n ) χ ν ( n ) P m , k ( p , c ) 1 P m , k ( p , c )
κ ( n + 1 ) = κ ( n ) χ κ ( P m , k ( s ) + P c ( s ) ) 1 P m , k ( s ) + P c ( s )
α 1 ( n + 1 ) = α 1 ( n ) χ α 1 ( n ) P m , k ( p , c ) τ P m , k ( p , c )
α 2 ( n + 1 ) = α 2 ( n ) χ α 2 ( n ) ( P m , k ( s ) + P c ( s ) ) R m , k ( s ) P m , k ( s ) + P c ( s ) ,
where χ denotes the learning rate and τ = min ( R m , k ( p , s ) , R m , k ( p , r ) ) . Detailed steps of our algorithm are shown in Algorithm 1.
Algorithm 1. Power Allocation Optimization
1:Set tollerance ζ and maximum number of iterations n max
Let n = 0 , choose ψ ( 0 , 1 ) , χ ( 0 , 1 )
Choose arbitrarily P m , k ( p , s ) ( 0 ) , P m , k ( p , r ) ( 0 ) , P m , k ( s ) ( 0 ) such that satisfies
power constraints
Calculate R m , k ( p , s ) ( 0 ) , R m , k ( p , r ) ( 0 ) and R m , k ( s ) ( 0 )
if R m , k ( p , s ) < R m , k ( p , r )
then τ = R m , k ( ( p , s )
else τ = R m , k ( p , r )
end
Let
ν ( n ) = 1 P m , k ( p , c ) ( n ) κ ( n ) = 1 P m , k ( s ) ( n ) + P c ( s ) α 1 ( n ) = τ P m , k ( p , c ) ( n ) α 2 ( n ) = R m , k ( s ) ( n ) P m , k ( s ) ( n ) + P c ( s )
Initialize the Lagrange multipliers λ , δ , ϵ , ξ , θ , μ
2:if τ = R m , k ( p , s )
update P m , k ( p , s ) and P m , k ( p , r ) with (24) and (25) and
obtain P m , k ( p , s ) * and P m , k ( p , r ) *
else
update P m , k ( p , s ) and P m , k ( p , r ) with (28) and (29) and
obtain P m , k ( p , s ) * and P m , k ( p , r ) *
end
Update P m , k ( s ) with (26) and obtain P m , k ( s ) *
3:Update the Lagrange multipliers using (31)–(36)
4:Compute ϕ * with P m , k ( p , s ) * , P m , k ( p , r ) * and P m , k ( s ) *
if | ϕ * ϕ ( n ) | ζ
P m , k ( p , s ) ( n + 1 ) = P m , k ( p , s ) * , P m , k ( p , r ) ( n + 1 ) = P m , k ( p , r ) * , P m , k ( s ) ( n + 1 ) = P m , k ( s ) *
exit program (convergence reached)
else
P m , k ( p , s ) ( n ) = P m , k ( p , s ) * , P m , k ( p , r ) ( n ) = P m , k ( p , r ) * , P m , k ( s ) ( n ) = P m , k ( s ) *
ϕ ( n ) = ϕ *
update ν , κ , α 1 , α 2 using (37)–(40)
Update R m , k ( p , s ) ( n + 1 ) , R m , k ( p , r ) ( n + 1 ) , R m , k ( s ) ( n + 1 )
n n + 1
end
go to step 2 unless n reached n max
For the problem in (15), a similar approach can be followed. By transforming it using Theorem 1 and applying the Lagrange dual method, the optimal transmit power for the direct mode can be obtained.

4.3. Transmission Mode and Relay Selection Subproblem

Given the optimal transmit power of the PUs and the SUs, the total energy efficiency of the network can be calculated as:
η * = m = 1 M k = 1 K ( β m ( d ) η m * + β m , k ( c ) η m , k * ) .
The transmission mode and relay selection subproblem can be formulated as:
max [ β m ( d ) , β m , k ( c ) ] η * s . t . C 1 : β m ( d ) + k = 1 K β m , k ( c ) 1 , C 2 : m = 1 M β m , k ( c ) 1 ,
which is a linear binary optimization problem, that can be solved using graph-based optimization method [44]. To implement such optimization problem, we construct two tables: Table 1a, whose columns contain the maximum energy efficiency when the mth PU chooses direct transmission mode; Table 1b, which contain the optimal energy efficiency when the mth PU chooses the kth SU as a relay node for cooperation mode.
It can be seen that the optimal transmission mode and relay selection solution of (42) is equivalent to finding the maximum sum of the energy efficiency elements which are chosen from various rows and columns, except for the direct mode. From (42), we can see that when the energy efficiency of the cell where only direct transmission is available is larger than all the energy efficiency of a cell where relay transmission is adopted, the PU should choose the direct mode. Hence, to solve the optimization problem formulated in (42), for the mth row, 1 m M , we first compare the element in direct mode with the other elements in the same row number in Table 1b; if the energy efficiency of the PU obtained in direct mode is the largest, then we set β m ( d ) = 1 , β m , k ( c ) = 0 , and delete the corresponding row in Table 1b. For the remaining, we can solve the optimal relay selection subproblem for cooperation transmission mode.
Given the constraints on both PUs and SUs, the optimal relay selection subproblem can be described as a bipartite graph, and the problem of finding the optimal relay can be regarded as an optimal matching algorithm in the bipartite graph that can be solved by a modified Kuhn-Munkres (K-M) algorithm [45,46]. A brief description of this algorithm and of its complexity is presented in Appendix B.
Based on the optimal relay selection results, we obtain the corresponding energy efficiency for cooperative mode, which is compared with that obtained from the direct mode for each PU. If the latter is larger, we set in direct mode the corresponding PU and rerun the K-M algorithm until the optimal energy efficiency obtained for cooperative mode is larger than that from the direct mode for all the PUs. In Algorithm 2, we present a brief description of the algorithm.
Algorithm 2. Transmission Mode and Relay Selection Optimization
1:Solving optimal power allocation subproblem
to obtain η m * and η m , k * , construct Table 1
2:for the mth PU
   if η m * > η m , k * , 1 k K
   then
      set β m ( p , d ) = 1
      delete the mth row in Table 1b
end
end
Table 1b is updated
3:repeat
4:Apply K-M algorithm on Table 1b to find
optimal relay node obtaining β m , k ( c )
5:if β m , k ( c ) η m , k * < η m * , 1 m M
      set β m ( p , d ) = 1
      set β m , k ( c ) = 0
Update Table 1b removing the mth row
end
6:until β m , k ( c ) η m , k * > η m * ,∀ PU m in Table 1b.

5. Results

In this section, we evaluate the system performance and effectiveness of the proposed solution throughout Monte Carlo simulations. A single PUs’ cell of 250 m radius is considered (see Figure 1), and a summary of the scenario parameters are shown in Table 2. Results are obtained averaging over 5000 independent configuration snapshots, each one with different PUs’ and SUs’ positions within the cell (PUs and SUs are uniformly distributed within the area), and independent realizations of shadowing and fading. Regarding the time slot subdivision between PU and SU, we initially assume t 1 = T / 2 , which is the choice maximizing energy efficiency, as shown in Section 5.6.

5.1. Convergence of the Iterative Algorithm

In Figure 3, we show the energy efficiency in deciBel, that is, η ( dB ) = 10 log 10 η , for different values of PUs’ and SUs’ circuit power, P c , as a function of the number of iterations, to illustrate the convergence of the proposed algorithm. It can be seen that the iterative algorithm generates a nondecreasing sequence and converges to a stable point within 15 iterations. In addition, we can observe the effect of P c to energy efficiency. As expected, the energy efficiency decreases with the increase of P c .

5.2. Energy Efficiency Versus Maximum Transmit Power

In Figure 4, the average energy efficiency varying the maximum transmit power, P max , is presented. The number of PUs and SUs is chosen randomly at every adaptation with a maximum number of 10 PUs and 10 SUs. From the figure, it can be observed that the energy efficiency increases with an increase of P max up to approximately 0.1 W, beyond which the energy efficiency saturates. Besides, we compare the performance of our solution with three different algorithms:
  • A random choice algorithm in which PUs choose cooperative relays randomly.
  • A non-cooperative approach in which a similar system model is considered but PUs choose relays to increase their own energy efficiency, instead of maximizing the total energy efficiency of the system.
  • A system in which there are only PUs, and therefore only direct transmission is possible.
As can be observed, the proposed algorithm outperforms the other three solutions. In addition, compared to a system where only PUs are present, we can show that with our approach having the same available spectrum, higher energy efficiency can be achieved, and also more users are able to use the resources available.

5.3. Energy Efficiency Versus Number of Pus

Figure 5 illustrates the effect that an increasing number of PUs have on energy efficiency when the number of SUs is set to 10. When the number of PUs is smaller than the number of SUs, we can observe a higher energy efficiency achieved. However, increasing the number of PUs, the energy efficiency decreases because of the increased contention of PUs to choose the most convenient cooperative relay. When the number of PUs becomes larger than the number of SUs, the decline of the curve becomes sharper. This comes from the fact that not every PU can choose a cooperative relay and is bounded to perform a direct transmission. As we have shown before, the energy efficiency based only on direct transmission is lower than the one based on the cooperative relay.

5.4. Energy Efficiency Versus Path Loss Exponent

In Figure 6, the average energy efficiency varying the path-loss exponent, γ , is presented. With the increase of γ , the energy efficiency decreases; this is due to the increased attenuation, which reflects in a lower data rate. Also, it can be seen that our algorithm outperforms the energy efficiency obtained from the scheme proposed in Reference [18], in which only PUs’ power is optimized, the non-cooperative approach, as well as a system where only PUs are present. In harsh propagation environments with γ > 3.5 , the energy efficiency of the proposed algorithm outperforms the direct transmission scenario by approximately 110% and the proposed scheme in Reference [18] by 10%.

5.5. Energy Efficiency Versus Shadowing Standard Deviation

In Figure 7, we illustrate the average energy efficiency varying the shadowing standard deviation. As can be seen, also considering the previous figure, the proposed algorithm is more robust to shadowing and attenuation hence offers better performance in harsh propagation environments. The proposed scheme outperforms the direct transmission scenario by 72% and the proposed scheme in Reference [18] by nearly 20%.

5.6. Energy Efficiency Versus Time Slot Division

In Figure 8, the average energy efficiency varying the time slot division chosen in the two relay transmission links is illustrated. As we notice in the figure, when both links have the same duration, we can achieve the highest value of energy efficiency. In addition, it can be observed that the variation between the maximum, for t 1 / T = 0.5 , and minimum value, for t 1 / T = 0.9 , of energy efficiency is not very high. This fact confirms the robustness of the proposed algorithm to the choice of the time slot division in the relay transmission.

6. Conclusions

In this paper, a novel energy-efficient architecture for CRN is proposed where each PU can choose one SU as a relay node. To encourage the cooperative behavior of the SUs, PUs lease a fraction of their allocated spectrum to the relay SUs to transmit their data. We proposed a centralized resource management network architecture to achieve a performance enhancement of the network. The resource allocation problem is formulated as a constrained sum energy efficiency maximization problem. Our analysis shows that this problem can be divided into two subproblems: (1) power allocation, (2) transmission mode and relay selection. To tackle the sum-of-ratio form of the power allocation subproblem, we transform it into a subtractive form optimization problem, and then closed-form optimal powers are calculated through a two-layer optimization. A modified K-M algorithm bipartite matching algorithm is used to solve the transmission mode and relay selection subproblem. System simulations performed in a typical multiuser case show that the proposed algorithm converges to the solution within a small number of iterations. In addition, simulations show an increment in the energy efficiency of the PUs network performance compared to previously proposed algorithms when channel conditions get worse. In addition our proposed algorithm shows a robustness to the choice of the time division in the relay transmission. Furthermore, in the system simulations we illustrated the effect that an increasing number of PUs have on the energy efficiency of the cell. Future work will deal with the use of 3D image-based algorithms to acquire more realistic channel characteristics in more realistic scenarios. In addition, we will consider comparing our algorithm with some existing non-CCRNs-based schemes.

Author Contributions

Conceptualization, D.B.; Investigation, D.B.; Software, D.B.; Supervision, A.G., E.M.V. and V.D.-E.; Validation, D.B., A.G., E.M.V. and V.D.-E.; Visualization, A.G., E.M.V., V.D.-E.; Writing—original draft, D.B.; Writing—review and editing, A.G., E.M.V. and V.D.-E. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Italian Ministry of Education, University and Research (MIUR) under the program “Dipartimenti di Eccellenza (2018–2022)—Precise CPS”, and by the EU COST action CA15104 “Inclusive Radio Communications (IRACON)”.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of Theorem 1

To solve the problem in (19), Lagrange dual method is used [41]. For this purpose, we first need the Lagrange function of this problem. The Lagrangian can be written as:
L ( P m , k ( p , r ) , P m , k ( p , s ) , P m , k ( s ) , α 1 , α 2 , λ , δ , ϵ , ξ , θ , μ , ν , κ ) = α 1 + α 2 + λ ( P m ( p , max ) P m , k ( p , s ) ) + δ ( P k ( s , max ) P m , k ( p , r ) P m , k ( s ) ) + ϵ ( R m , k ( p , r ) R m , k ( p , s ) ) + ξ ( R m , k ( p , s ) R m ( p , min ) ) + δ ( R m , k ( p , r ) R m ( p , min ) ) + μ ( R m , k ( s ) R k ( s , min ) ) + ν ( R m , k ( p , s ) α 1 ( P m , k ( p , c ) ) ) + κ ( R m , k ( s ) α 2 ( P m , k ( s ) + P c ( s ) ) ) .
If the set of solutions ( P m , k ( p , s ) * , P m , k ( p , r ) * , P m , k ( s ) * , α 1 * , α 2 * ) is the solution of the problem in (19), then exist the Lagrange multipliers λ ¯ , δ ¯ , ϵ ¯ , ξ ¯ , θ ¯ , μ ¯ , ν ¯ , κ ¯ in order to satisfy the KKT conditions [39] for the problem (19). The conditions are as follows:
L P m , k ( p , s ) = ν ¯ ( P m , k ( p , s ) α 1 * ) λ ¯ ϵ ¯ R m , k ( p , s ) = 0 L P m , k ( p , r ) = ( ϵ ¯ + θ ¯ ) R m , k ( p , r ) ν ¯ α 1 * δ ¯ = 0 L P m , k ( s ) = ( μ ¯ + κ ¯ ) R m , k ( s ) κ ¯ α 2 * δ ¯ = 0
L α 1 = 1 ν ¯ P m , k ( p , c ) = 0
L α 2 = 1 κ ¯ ( P m , k ( s ) + P c ( s ) ) = 0
L λ = P m p , m a x P m , k ( p , s ) = 0
L δ = P k s , m a x P m , k ( p , r ) P m , k ( s ) = 0
L ϵ = R m , k ( p , r ) R m , k ( p , s ) = 0
L ξ = R m , k ( p , s ) R m ( p , min ) = 0
L θ = R m , k ( p , r ) R m ( p , min ) = 0
L ξ = R m , k ( s ) R k ( s , min ) = 0
L ν = R m , k ( p , s ) α 1 * P m , k ( p , c ) = 0
L κ = R m , k ( s ) α 2 * ( P m , k ( p , c ) + P c ( s ) ) = 0
Considering that P m , k ( p , c ) > 0 and P m , k ( p , c ) + P c ( s ) > 0 then from (A3) and (A4) we have
ν ¯ = 1 P m , k ( p , c ) κ ¯ = 1 P m , k ( s ) + P c ( s )
and (A11) and (A12) are equivalent to
α 1 * = R m , k ( p , s ) P m , k ( p , c ) α 2 * = R m , k ( s ) P m , k ( s ) + P c ( s ) .
The system Equations (A2), (A5)–(A10) are the KKT conditions of the optimization problem in (20) for ν = ν ¯ , κ = κ ¯ , α 1 = α 1 * , α 2 = α 2 * .

Appendix B. Summary of the Kuhn-Munkres Algorithm

For the sake of completeness, we briefly introduce some definitions and a theorem regarding the K-M algorithm.
Bipartite: A bipartite graph is a graph whose vertices can be divided into two disjoint sets U and V, such that every edge connects a vertex in U to one in V. It can be represented as G = ( U , V , E ) with E denoting the edges of the graph.
Weighted bipartite: A weighted bipartite is a bipartite in which each edge ( u , v ) has a weight factor w ( u , v ) .
Matching: A matching in a graph is a subset H G . If H and G share the same vertex set, then H is called a complete matching. The size of a matching is denoted as | H | which equals to the number of edges in H.
Feasible vertex labeling: A feasible vertex labeling in G is a real-valued function l on U V such that for u U and v V ,
l ( u ) + l ( v ) w ( u , v ) .
Equality subgraph: If l is a feasible labeling, we denote a subgraph of G as G l which contains a number of edges and the endpoints of these edges. If the edges of G l meet the condition l ( u ) + l ( v ) w ( u , v ) , then G l is called the equality subgraph for l.
Theorem A1.
If l is a feasible vertex labeling for G, and H is a complete matching of U to V with H G l , then H is an optimal assignment of U to V.

Appendix B.1. Solving Optimal Relay Selection Problem Based on K-M Algorithm

Applying K-M algorithm to solve the optimal relay selection problem of the PUs, a weighted bipartite graph G with a bipartite division G 0 = ( U , V , E ) is constructed, where the set of vertices U represents the set of PUs, that is, U = [ PU 1 , PU 2 , , PU M ] and the set of vertices V represents the set of SUs, V = [ SU 1 , SU 2 , , SU K ] . The weight of the edge ( PU m , SU k ) in the weighted bipartite graph can be defined as the joint energy efficiency of the mth PU and the kth SU, that is, η m , k ( p , c ) + η m , k ( s ) , m = 1 , , M and k = 1 , , K .
The steps of solving the optimal relay selection problem based on K-M algorithm can be described as follows:
  • Find initial feasible vertex labeling and determine G l 0 and choose an arbitrary matching H in G l 0 .
  • If H is a maximum matching for G, then the optimization problem is solved. Otherwise, the label having not being allotted by the distribution H is selected in G l 0 . Set S = U , and T = Φ , which denotes the empty set.
  • N G l 0 ( S ) denotes the collection of points which connect with S in G l 0 . If N G l 0 ( S ) T , go to step (2). Otherwise, N G l 0 ( S ) = T . Find
    Δ = min ( l ( u ) + l ( v ) w ( u , v ) | u S , v V T )
    and replace existing labeling l with l by
    l ( u ) = l ( u ) Δ , u S l ( u ) + Δ , u T l ( u ) , otherwise .
The process continues until an equal subgraph consisting a complete match is obtained.

Appendix B.2. Complexity of the K-M Algorithm

Theoretically, the K-M algorithm is guaranteed to reach the global optimal solution [47]. It is the first proved polynomial algorithm and it was observed to be a strongly polynomial by Munkres in Reference [48]. Further, Edmonds and Karp in Reference [49] noticed that it can be modified to O ( n 3 ) running time.
In our case in each of the steps (1) or (2) the algorithm adds one edge to the matching and this happens O ( | V | ) . In addition, it needs O ( | V | ) to find the right vertex (if there is one) and improving the labelling takes O ( | V | ) time to find delta and to update the labelling accordingly. We might have to improve the labelling O ( | V | ) times until we matching occurs. This makes for a total of O ( | V 2 | ) . In all, there are O ( | V | ) iterations each taking O ( | V | ) work, leading to a total running time of O ( | V 3 | ) .

References

  1. Bartoli, G.; Marabissi, D.; Pucci, R.; Ronga, L.S. AI Based Network and Radio Resource Management in 5G HetNets. J. Signal Process. Syst. 2017, 89, 133–143. [Google Scholar] [CrossRef]
  2. Federal Communications Commission. ET Docket No. 03–322. Notice of Proposed Rule Making and Order. Available online: https://docs.fcc.gov/public/attachments/FCC-03-322A1.pdf (accessed on 28 October 2020).
  3. Bartoli, G.; Fantacci, R.; Gei, F.; Marabissi, D.; Micciullo, L. Dynamic PSS and commercial networks spectrum sharing in the time domain. Trans. Emerg. Telecommun. Technol. 2016, 28, e3083. [Google Scholar] [CrossRef]
  4. Mitola, J.; Maguire, G.Q. Cognitive radio: Making software radios more personal. IEEE Pers. Commun. 1999, 6, 13–18. [Google Scholar] [CrossRef] [Green Version]
  5. Sithamparanathan, K.; Giorgetti, A. Cognitive Radio Techniques: Spectrum Sensing, Interference Mitigation, and Localization; Artech House: Boston, MA, USA, 2012. [Google Scholar]
  6. Giorgetti, A.; Varrella, M.; Chiani, M. Analysis and performance comparison of different cognitive radio algorithms. In Proceedings of the 2009 Second International Workshop on Cognitive Radio and Advanced Spectrum Management, Aalborg, Denmark, 18–20 May 2009; pp. 127–131, Invited Paper. [Google Scholar]
  7. Fantacci, R.; Marabissi, D. Cognitive Spectrum Sharing: An Enabling Wireless Communication Technology for a Wide Use of Smart Systems. Future Internet 2016, 8, 23. [Google Scholar] [CrossRef] [Green Version]
  8. Aqeel, S.; Kok-Lim, Y. Spectrum Leasing in Cognitive Radio Networks: A Survey. Int. J. Distrib. Sens. Netw. 2014, 10, 329235. [Google Scholar]
  9. Amin, O.; Ikki, S.S.; Uysal, M. On the performance analysis of multi-relay cooperative diversity systems with channel estimation errors. IEEE Trans. Veh. Technol. 2011, 60, 2050–2059. [Google Scholar] [CrossRef]
  10. Amin, O.; Lampe, L. Opportunistic Energy Efficient Cooperative Communication. IEEE Wirel. Commun. Lett. 2012, 1, 412–415. [Google Scholar] [CrossRef]
  11. Amin, O.; Ikki, S.S.; Uysal, M. Adaptive bit loading for multi-relay cooperative orthogonal frequency division multiple with imperfect channel estimation. IET Commun. 2012, 6, 1821–1828. [Google Scholar] [CrossRef]
  12. Nie, P.Y.; Zhang, P. A note on Stackelberg games. In Proceedings of the Chinese Control and Decision Conference, Yantai, China, 2–4 July 2008; pp. 1201–1203. [Google Scholar]
  13. Scutari, G.; Pang, J. Joint Sensing and Power Allocation in Nonconvex Cognitive Radio Games: Nash Equilibria and Distributed Algorithms. IEEE Trans. Inf. Theory 2013, 59, 4626–4661. [Google Scholar] [CrossRef]
  14. Xu, H.; Gao, H.; Zhou, C.; Duan, R.; Zhou, X. Resource Allocation in Cognitive Radio Wireless Sensor Networks with Energy Harvesting. Sensors 2019, 19, 5115. [Google Scholar] [CrossRef] [Green Version]
  15. Elkourdi, T.; Simeone, O. Spectrum Leasing via Cooperation With Multiple Primary Users. IEEE Trans. Veh. Technol. 2012, 61, 820–825. [Google Scholar] [CrossRef]
  16. Toroujeni, S.M.M.; Sadough, S.M.; Ghorashi, S.A. Spectrum Leasing for OFDM-Based Cognitive Radio Networks. IEEE Trans. Veh. Technol. 2013, 62, 2131–2139. [Google Scholar] [CrossRef]
  17. Hafeez, M.; Elmirghani, J.M.H. Dynamic Spectrum Leasing for Bi-Directional Communication: Impact of Selfishness. IEEE Trans. Commun. 2016, 64, 2427–2437. [Google Scholar] [CrossRef] [Green Version]
  18. Tang, Y.; Grace, D.; Wang, L. Spectrum Leasing and Relay Selection in Multiuser Cooperative Cognitive Radio Networks. J. Commun. 2014, 9, 198–207. [Google Scholar] [CrossRef]
  19. Liu, Z.; Zhao, M.; Xie, Y.; Li, X.; Ma, K. Relay selection and power allocation for spectrum leasing in cooperative cognitive radio networks. In Proceedings of the Chinese Control Conference (CCC), Dalian, China, 26–28 July 2017; pp. 8958–8963. [Google Scholar]
  20. Zhai, C.; Zhang, W. Adaptive Spectrum Leasing with Secondary User Scheduling in Cognitive Radio Networks. IEEE Trans. Wirel. Commun. 2013, 12, 3388–3398. [Google Scholar] [CrossRef]
  21. Bastami, A.H.; Kazemi, P. Cognitive Multi-Hop Multi-Branch Relaying: Spectrum Leasing and Optimal Power Allocation. IEEE Trans. Wirel. Commun. 2019, 18, 4075–4088. [Google Scholar] [CrossRef]
  22. Ma, Y.; Yin, X.; Yang, X.; Liu, Q. Resource allocation for spectrum-leasing based CRN with delay-sensitive traffic. EURASIP J. Wirel. Commun. Netw. 2017, 108. [Google Scholar] [CrossRef] [Green Version]
  23. Yin, S.; Qu, Z.; Wang, Z.; Li, L. Energy-Efficient Cooperation in Cognitive Wireless Powered Networks. IEEE Commun. Lett. 2017, 21, 128–131. [Google Scholar] [CrossRef]
  24. Chatterjee, S.; Maity, S.P.; Acharya, T. Energy-Spectrum Efficiency Trade-Off in Energy Harvesting Cooperative Cognitive Radio Networks. IEEE Trans. Cogn. Commun. Netw. 2019, 5, 295–303. [Google Scholar] [CrossRef]
  25. Kong, H.Y.; Koo, I. Opportunistic relaying based spectrum leasing for cognitive radio networks. J. Commun. Netw. 2011, 13, 50–55. [Google Scholar] [CrossRef]
  26. Kader, M.F.; Chowdhury, M. Cooperative secondary user selection as a relay for the primary system in underlay cognitive radio networks. In Proceedings of the 15th International Conference on Computer and Information Technology (ICCIT), Chittagong, Bangladesh, 22–24 December 2012; pp. 275–278. [Google Scholar]
  27. Shao, C.; Roh, H.; Lee, W. Aspiration Level-Based Strategy Dynamics on the Coexistence of Spectrum Cooperation and Leasing. IEEE Commun. Lett. 2014, 18, 70–73. [Google Scholar] [CrossRef]
  28. Marabissi, D.; Bartoli, G.; Fantacci, R.; Micciullo, L. Energy efficient cooperative multicast beamforming in ultra dense networks. IET Commun. 2017, 12, 573–578. [Google Scholar] [CrossRef]
  29. Bilibashi, D.; Chai, R. Energy Efficiency Based Resource Allocation Scheme for Cooperative Cognitive Radio Networks. Int. J. Electron. Electr. Eng. 2017, 5, 7–14. [Google Scholar] [CrossRef]
  30. Jagannathan, R. On some properties of programming problems in parametric form pertaining to fractional programming. Manag. Sci. 1966, 12, 609–615. [Google Scholar] [CrossRef]
  31. Dinkelbach, W. On nonlinear fractional programming. Manag. Sci. 1967, 13, 492–498. [Google Scholar] [CrossRef]
  32. Dür, M.; Horst, R.; Thoai, N. Solving sum-of-ratios fractional programs using efficient points. Optimization 2001, 49, 447–466. [Google Scholar] [CrossRef]
  33. Darsena, D.; Gelli, G.; Melito, F.; Verde, F. Performance Analysis of Amplify-and-Forward Multiple-Relay MIMO Systems with ZF Reception. IEEE Trans. Veh. Technol. 2015, 64, 3274–3280. [Google Scholar] [CrossRef]
  34. Verde, F.; Scaglione, A.; Darsena, D.; Gelli, G. An Amplify-and-Forward Scheme for Spectrum Sharing in Cognitive Radio Channels. IEEE Trans. Wirel. Commun. 2015, 14, 5629–5642. [Google Scholar] [CrossRef]
  35. Amin, O.; Uysal, M. Adaptive power loading for multi-relay OFDM Regenerative networks with relay selection. IEEE Trans. Commun. 2012, 60, 614–619. [Google Scholar] [CrossRef]
  36. Bagheri, S.; Verde, F.; Darsena, D.; Scaglione, A. Randomized Decode-and-Forward Strategies for Two-Way Relay Networks. IEEE Trans. Wirel. Commun. 2011, 12, 4214–4225. [Google Scholar] [CrossRef]
  37. Bertsekas, D.P. Nonlinear Programming, 2nd ed.; Athena Scientific: Belmont, MA, USA, 1999. [Google Scholar]
  38. He, S.; Huang, Y.; Yang, L.; Ottersten, B. Coordinated Multicell Multiuser Precoding for Maximizing Weighted Sum Energy Efficiency. IEEE Trans. Signal Process. 2014, 62, 741–751. [Google Scholar] [CrossRef]
  39. Borwein, J.; Lewis, A. Convex Analysis and Nonlinear Optimization: Theory and Examples; Springer: New York, NY, USA, 2006. [Google Scholar]
  40. Ng, D.W.K.; Lo, E.S.; Schober, R. Energy-Efficient Resource Allocation in OFDMA Systems with Large Numbers of Base Station Antennas. IEEE Trans. Wirel. Commun. 2012, 11, 3292–3304. [Google Scholar] [CrossRef]
  41. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  42. Boyd, S.; Mutapcic, A. Subgradient Methods, Notes for EE364; Standford University: Stanford, CA, USA, 2006. [Google Scholar]
  43. Lee, E.T.; Wang, J.W. Statistical Methods for Survival Data Analyses, 3rd ed.; John Wiley and Sons, Inc.: Hoboken, NJ, USA, 2003; pp. 428–432. [Google Scholar]
  44. Samanta, S.; Sarkar, B. Isomorphism on generalized fuzzy graphs and image visualizations. Soft Comput. 2020, 24, 14401–14409. [Google Scholar] [CrossRef]
  45. Kuhn, H. The hungarian method for the assignment problem. Naval Res. Logist. Q. 1955, 2, 83–97. [Google Scholar] [CrossRef] [Green Version]
  46. Mensah, K.K.; Chai, R.; Bilibashi, D.; Gao, F. Energy efficiency based joint cell selection and power allocation scheme for HetNets. Digit. Commun. Netw. 2016, 2, 184–190. [Google Scholar] [CrossRef]
  47. Hong, C.; Jingjing, Z.; Chunfeng, C.; Qinyu, C. Solving large-scale assignment problems by Kuhn-Munkres algorithm. In Proceedings of the International Conference on Advances in Mechanical Engineering and Industrial Informatics, Hangzhou, China, 9–10 April 2016. [Google Scholar]
  48. Bourgeois, F.; Lassalle, J.C. An extension of the munkres algorithm for the assignment problem to rectangular matrices. Commun. ACM 1971, 14, 802–804. [Google Scholar] [CrossRef]
  49. Edmonds, J.; Karp, R.M. Theoretical improvements in algorithmic efficiency for network flow problems. Commun. ACM 1972, 19, 248–264. [Google Scholar] [CrossRef]
Figure 1. Proposed architecture and scenario for the cell.
Figure 1. Proposed architecture and scenario for the cell.
Sensors 20 06161 g001
Figure 2. Bandwidth allocation and time slot division for primary user (PU) and secondary user (SU).
Figure 2. Bandwidth allocation and time slot division for primary user (PU) and secondary user (SU).
Sensors 20 06161 g002
Figure 3. Energy efficiency versus number of iterations with M = 10 and K = 10 , for different circuit power consumption P c .
Figure 3. Energy efficiency versus number of iterations with M = 10 and K = 10 , for different circuit power consumption P c .
Sensors 20 06161 g003
Figure 4. Energy efficiency varying the maximum transmit power. Three different algorithms are compared. M and K in this case change randomly from 1 to 10.
Figure 4. Energy efficiency varying the maximum transmit power. Three different algorithms are compared. M and K in this case change randomly from 1 to 10.
Sensors 20 06161 g004
Figure 5. Energy efficiency varying number of PUs. M varies from 2 to 18 and K = 10 .
Figure 5. Energy efficiency varying number of PUs. M varies from 2 to 18 and K = 10 .
Sensors 20 06161 g005
Figure 6. Energy efficiency increment with respect to direct mode varying the path loss exponent. M and K change randomly from 1 to 15 while σ s = 8dB.
Figure 6. Energy efficiency increment with respect to direct mode varying the path loss exponent. M and K change randomly from 1 to 15 while σ s = 8dB.
Sensors 20 06161 g006
Figure 7. Energy efficiency increment with respect to direct mode versus shadowing standard deviation. M and K in this case change randomly within 1 and 15 while γ = 3.
Figure 7. Energy efficiency increment with respect to direct mode versus shadowing standard deviation. M and K in this case change randomly within 1 and 15 while γ = 3.
Sensors 20 06161 g007
Figure 8. Energy efficiency versus time slot division in the relay transmission link. The parameters M and K are random variables uniformly distributed from 1 to 10.
Figure 8. Energy efficiency versus time slot division in the relay transmission link. The parameters M and K are random variables uniformly distributed from 1 to 10.
Sensors 20 06161 g008
Table 1. Tables of direct mode and cooperative mode.
Table 1. Tables of direct mode and cooperative mode.
(a) Direct Mode(b) Cooperative Mode
Direct ModeCoop. with SU 1 Coop. with SU 2 Coop. with SU K
PU 1 η 1 * η 1 , 1 * η 1 , 2 * η 1 , K *
PU 2 η 2 * η 2 , 1 * η 2 , 2 * η 2 , K *
PU M η M * η M , 1 * η M , 2 * η M , K *
Table 2. System parameters.
Table 2. System parameters.
Cell radius 250 m
Small scale fading distributionRayleigh fading
Carrier frequency 3.5 GHz
Bandwidth 50 MHz
Noise power per link, σ 2 90 dBm
Shadowing, s m , k ( p , s ) Log-normal with standard deviation of 8 dB
Average channel gain at 1 m, k 0 39 dB
Path-loss exponent, γ 3
Circuit power per UE (PU or SU), P c 20 dBm
Maximum transmit power per UE, P max 24 dBm
Minimum data rate constraint, R m ( p , min ) 100 Mbit/s
Fraction of bandwidth, ρ m 0.66
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bilibashi, D.; Vitucci, E.M.; Degli-Esposti, V.; Giorgetti, A. An Energy-Efficient Unselfish Spectrum Leasing Scheme for Cognitive Radio Networks. Sensors 2020, 20, 6161. https://0-doi-org.brum.beds.ac.uk/10.3390/s20216161

AMA Style

Bilibashi D, Vitucci EM, Degli-Esposti V, Giorgetti A. An Energy-Efficient Unselfish Spectrum Leasing Scheme for Cognitive Radio Networks. Sensors. 2020; 20(21):6161. https://0-doi-org.brum.beds.ac.uk/10.3390/s20216161

Chicago/Turabian Style

Bilibashi, Denis, Enrico M. Vitucci, Vittorio Degli-Esposti, and Andrea Giorgetti. 2020. "An Energy-Efficient Unselfish Spectrum Leasing Scheme for Cognitive Radio Networks" Sensors 20, no. 21: 6161. https://0-doi-org.brum.beds.ac.uk/10.3390/s20216161

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop