# Energy-Efficient Layered Video Multicast over OFDM-Based Cognitive Radio Systems

Authors: Wenjun Xu; Shengyu Li; Yue Xu; Zhiyong Feng; Jiaru Lin
Journal: International Journal of Distributed Sensor Networks
 (2015)

---

## Abstract

An energy-efficient layered video multicast (LVM) scheme for “bandwidth-hungry” video services is studied in OFDM-based cognitive radio (CR) systems, where the video data is encoded into a base layer and several enhancement layers with the former intended for all subscribers to guarantee the basic quality of reconstructed video and the latter aiming at the quality improvement for the promising users with good channel conditions. Moreover, in order to balance user experience maximization and power consumption minimization, a novel performance metric energy utility (EU) is proposed to measure the sum achieved quality of reconstructed video at all subscribers when unit transmit power is consumed. Our objective is to maximize the system EU by jointly optimizing the intersession/interlayer subcarrier assignment and subsequent power allocation. For this purpose, we first perform subcarrier assignment for base layer and enhancement layers using greedy algorithm and then present an optimal power allocation algorithm to maximize the achievable EU using fractional programming. Simulation results show that the proposed algorithms can adaptively capture the state variations of licensed spectrum and dynamically adjust the video transmission to exploit the scarce spectrum and energy resources adequately. Meanwhile, the system EU obtained in our algorithms is greatly improved over traditional spectrum efficiency (SE) and energy efficiency (EE) optimization models.


---

## Body

## 1. Introduction

Recently, the fifth generation (5G) mobile wireless system has been under heated discussion [1, 2]. It is reported that the wireless traffic volume will increase by 1000-fold over the next decade [3], and hence there is an urgent need to design novel spectrum-efficient transmission paradigms. Cognitive radio (CR) [4] is one of the best technologies to improve the spectrum efficiency (SE) and has attracted many researchers’ attention. The basic idea of CR is to bear data transmission among secondary users (SUs) by reusing licensed spectrum without harming the benefits of authorized users (also known as primary users, PUs). Orthogonal frequency division multiplexing (OFDM) supports a flexible spectrum management by dividing the available spectrum into fine-granularity subcarriers and hence is recognized as a promising technology for spectrum reusing [5]. As a result, it is meaningful to combine CR and OFDM together and investigate new transmission paradigm in OFDM-based CR systems.Meanwhile, energy efficiency (EE) is also a key metric for 5G, in which energy consumption needs to be reduced on the order of several magnitudes [6]. Note that extra power consumption incurred by spectrum sensing makes the power saving issue more critical in CR systems [7, 8]. Until now, large amounts of researches have been conducted to study energy-efficient transmission in OFDM-based CR systems. For example, in our early work [9], the EE metric measured by the achieved transmission bits per Joule is adopted, and the optimal power allocation for EE maximization is derived usingfractional programming. Subsequently, the model is improved with the minimum rate guarantee and subcarrier assignment taken into account in [10, 11], respectively. Then, the authors in [12] further consider channel uncertainty and study the EE maximization problem with a probabilistic interference control policy.However, all these researches mentioned above focus on unicast transmission. Along with the proliferation of smart phones, mobile multimedia services, especially mobile video services, have been in the explosive growth [13]. CR can effectively alleviate the more and more serious spectrum scarcity issue and is one of the key candidate technologies in 5G. Hence, it is almost an inevitable trend to deliver the increasing popular video services in the future communication system without licensed spectrum. Even for networks which have already been allocated with some spectrum bands, integrating CR function into the networks, for example, Licensed-Assisted Access (LAA) [14], can provide more competent video transmission, improving the user experience greatly. Moreover, as secondary networks can only access the authorized spectrum opportunistically, its transmission capacity is limited by the prioritized access mode. Hence, it is more challenging to transmit “bandwidth-hungry” video services with QoS guarantee in CR networks. The achieved results can provide a good guidance for how to bear other types of services in CR networks.Video multicast has become an indispensable part for wireless networks, and hence it is of great significance to study how to scalably multicast video in CR systems. Several kinds of multicast schemes have been proposed for video transmission in the literature, including conventional multicast (CM) [13, 15], multiple description coding multicast (MDCM) [16, 17], and layered video multicast (LVM) [18, 19]. In CM [15], all subscribers in a multicast group receive the intended content with the identical quality, and the transmission rate is limited by the least channel gain of all subscribers. To cope with this issue, MDCM and LVM introduce source coding to support distinguished video transmission for different subscribers. In MDCM [16, 17], the video data is encoded into multiple descriptions and transmitted at different rates. For subscribers of various channel conditions, different sets of descriptions are received to jointly recover video with different resolutions. Despite being attractive in terms of system throughput, MDCM cannot guarantee the successful reception of key information and hence applies poorly in practice. By comparison, in LVM [18, 19], the video data is encoded into a base layer (BL) and several enhancement layers (ELs). The BL is intended to all subscribers at a low rate and hence can guarantee a basic recovered video quality, while the ELs are transmitted at incremental rates and opportunistically received by subscribers with promising channel conditions to persistently improve the video quality.To the best our knowledge, existing researches on video multicast in OFDM-based CR systems mainly focus on CM [13, 15] or MDCM [16, 17], and what is more, only SE maximization model is studied. In this paper, joint intersession/interlayer subcarrier assignment and power allocation problem for energy-efficient LVM are studied in OFDM-based CR systems. For existing energy-efficient transmission models [9–12],full-buffer traffic model is assumed; that is, there is always infinite data waiting for transmission, and EE maximization is studied with no consideration of service characteristic. For LVM, the final recovered video quality is not linear with the receiving rate. Therefore, a utility function is introduced in our model to depict the relationship between the user experience, that is, recovered video quality, and its receiving rate, and a more accurate performance metric,energy utility (EU), is designed to measure the sum user experience achieved per Joule.In detail, the main contributions of this paper can be summarized as follows.(1)
EU-Based Optimization Model. A novel optimization model is established for energy-efficient LVM in OFDM-based CR systems which aims at maximizing the system EU to balance the total recovered video quality and power consumption while guaranteeing multiple interference constraints for PUs.(2)
Spectrum Assignment Method for LVM. Both BL and EL subcarrier assignment algorithms are proposed to execute the intersession/interlayer subcarrier assignment in LVM. The BL subcarrier assignment aims at guaranteeing the basic video qualities for all multicast sessions with as fewer subcarriers as possible, while the EL subcarrier assignment tries to optimize the system EU by assigning subcarriers to the proper subscribers.(3)
Optimal Power Allocation Method for EU Maximization. For the multiconstrained EU-maximization problem, an optimal power allocation algorithm is presented by jointly utilizing fractional programming and subgradient method, which can be considered as the framework of optimizing energy-aware video multicast in OFDM-based CR systems.The rest of the paper is organized as follows. In Section2, we build the EU-based optimization model for LVM over OFDM-based CR systems. The subcarrier assignment for BL as well as EL and the EU-based power allocation are proposed and discussed by Section 3. Finally, simulation results are shown in Section 4, and conclusions are drawn in Section 5.


## 2. System Model and Problem Formulation

In this section, the spectrum division manner and the mutual-interference model are depicted firstly, and then the energy-efficient LVM transmission model is formulated with the objective of EU maximization.


### 2.1. OFDM-Based CR System

The considered CR system is composed of a primary network and a CR network, which are both deployed in a cellular fashion. The primary network consists of a primary base station (BS) andL PUs, while the CR network is made up of a secondary BS and K SUs. The whole licensed spectrum spans B Hz, and each PU occupies a disjoint fraction of the spectrum, denoted as Bl. For the finer-grained spectrum utilization, the OFDM technology is adopted to divide the whole sensed spectrum B into N subcarriers, with each subcarrier spanning B0=B/N Hz [5, 9]. Figure 1 illustrates the details of CR system and spectrum distribution.Figure 1
System scenario and spectrum distribution.
(a)
CR system: a primary network and a CR network coexist
(b)
Spectrum division by OFDM and spectrum occupation of PUsSince the sensed spectrum by CR network is licensed to the primary network, the privilege of PUs to use the spectrum must be guaranteed, which typically necessitates the interference control. Therefore, a common interference evaluation model is introduced from [20].In [19], the power spectrum density (PSD) of subcarrier n is written as(1)φnf=pnTssin⁡πfTsπfTs2,where pn is the transmitted power of signal on subcarrier n and Ts is the duration of OFDM symbol. Then, the interference introduced by the signal on subcarrier n into PU l can be expressed as(2)I~l,nSP=GlSP∫dl,n-Bl/2dl,n+Bl/2φnfdf=Il,npn,where dl,n denotes the spectral distance between subcarrier n and the center frequency of PU l, GlSP denotes the channel gain from the secondary BS to PU l, and Il,n represents the interference caused by the normalized power on subcarrier n to PU l.Meanwhile, the interference caused by the primary BS to SUk on subcarrier n is calculated as [20](3)I~k,nPS=∑l=1LGk,nPS∫dl,n-Δf/2dl,n+Δf/2EIlωdω,where Gk,nPS denotes the channel gain from primary BS to SU k on subcarrier n, and E[·] is the expectation operator. E[Il(ω)] is the signal PSD of PU l after N-Fast-Fourier-transform (FFT) processing, and it is represented as(4)EIlω=12πN∫-ππΦlejϕsin⁡ω-ϕN/2sin⁡ω-ϕ/22dϕ,where Φl(ejϕ) is the PSD of the transmitted signal by the primary BS to PU l, which is generally modeled as an elliptically filtered white noise process [20].


### 2.2. LVM Transmission Model

As shown in Figure1(a), all K users are partitioned into G multicast groups according to the video contents they are interested in. The set and number of SUs in group g are denoted as Kg and Kg, respectively. Let Hk∣g,nSS denote the channel gain between the secondary BS and SU k in group g, and the corresponding channel gain to interference-plus-noise ratio (CINR) can be expressed as(5)γk∣g,n=Hk∣g,nSSσn2+I~k∣g,nPS,where σn2 is the noise power on subcarrier n and I~k∣g,nPS is calculated according to (3) with an index mapping from k∣g to k. Based on the Shannon formula, the achievable data rate of SU k in group g on subcarrier n is calculated by(6)rk∣g,n=B0log2⁡1+pnγk∣g,n.For CM, the lowest rate of SUs in a group is conservatively adopted to ensure the correct data reception of all SUs [15]. Thus, the data transmission rate of group g on subcarrier n is expressed as(7)rg,nc=mink∈Kg⁡rk∣g,n=B0log2⁡1+pnmink∈Kg⁡γk∣g,n.In wireless transmission, the heterogeneity of the receiving channels for different SUs will seriously limit the performance for CM. In order to overcome the shortcomings, LVM [18, 19] is introduced and modeled in this paper. In LVM, video data are transmitted resiliently on different subcarriers to adapt to the diverse channel conditions. This is accomplished by coding the source data into a BL and several ELs, and as long as the BL is received, SU can decode the video stream with the basic quality. If more ELs are received, the decoded video quality is increasingly improved [21]. The essential difference between CM and LVM lies in that the former requires that all SUs in a group receive the video with the identical quality, whilst the latter allows the differential reception, depending on the individual channel quality. Therefore, LVM provides a new degree of freedom, that is, the transmission rates on subcarriers, to exploit different channel conditions.Specifically, it is assumed that each videog, which is received by group g, is encoded into one BL with rate Rgb and one EL with rate Rge (see Figure 1(a)). Rgmax=Rgb+Rge denotes the maximal rate for video g. As in [19], we also assume that the fine grained scalability (FGS) coding technique is adopted, so that the EL can be truncated at any bit location with all the remaining bits still being useful at the decoder [22].In LVM, the BL data is of great importance to reconstruct the source video, and it is imperative that the BL can be received by all the SUs in a group. Hence, if subcarriern is used by group g for the BL transmission, the transmission rate is equal to that of CM; that is,(8)rg,nb=rg,nc.If subcarrier n is used by group g for the EL transmission at a rate of rg,ne, the receiving rate of SU k,k∈Kg on subcarrier n is written as(9)σrk∣g,n≥rg,nerg,ne,where σ(T) is 1 if T is true, otherwise 0. Define two indicators ρg,n and sn to perform the intersession and interlayer subcarrier assignment as follows. ρg,n∈{1,0} denotes whether subcarrier n is assigned to group g or not, and sn denotes whether subcarrier n is used for the BL transmission (sn=1) or the EL transmission (sn=0). The achieved rate of SU k,k∈Kg is calculated as(10)Rk=∑n=1Nρg,n1-snσrk∣g,n≥rg,nerg,ne+snrg,nb.


### 2.3. EU-Based Problem Formulation

In LVM, the achieved rate cannot reflect the quality directly, and some metrics, such as peak signal to noise ratio (PSNR) and mean square error (MSE), may be more accurate to evaluate the received video quality. For generalizing the expression, the utility functionUk(Rk) is defined to denote the change relationship of reconstructed video quality with the achieved rate. Moreover, Uk(Rk) is assumed to be nondecreasing and concave according to the measured or theoretical results [18, 19]. Then, the total weighted-utility is summed as(11)Utotal=∑k=1KwkUkRk,where wk(wk>0) is the weight of SU k and can be used to reflect the relative priorities among SUs.Apart from the video quality, the energy cost of video transmission over OFDM-based CR systems should be also considered. The total energy consumption in a timeslot includes three parts: the sense energyPs, that is, the energy consumed to sense the available spectrum, the transmission energy Pt, that is, the energy consumed to transmit data on the sensed spectrum, and the circuit energy Pc. Note that the concepts of “energy” and “power” are not strictly distinguished in this paper, since only the constant duration is multiplied or not. The total energy consumption is accumulated as(12)Etotal=Pt+Ps+Pc=∑n=1Npn+ξ∑n=1Npn+Pc,where ξ denotes the ratio of the sense energy Ps to the transmission energy Pt [15]. The formula accounts for the required energy consumption for downlink (DL) video multicast in OFDM-based CRNs. In practice, the total system power consumption includes the power consumption at both BS and end device sides. The power consumption of end devices during DL transmission mainly results from multiple active circuit modules, such as the channel estimation and feedback module and the baseband signal processing module. However, the aggregate value is relatively small compared with the power consumption of BS, whose radio frequency (RF) transmit power is often very large due to the large-scale pathloss fading. As a result, the power consumption of end devices is not considered in formula (12), in which the total system energy consumption is approximated as the power consumption at the secondary BS [7, 9–12].In the existing work, EE is maximized to optimize the achieved rate by the unit energy consumption [9–12]. For the video transmission, it makes more sense to shift EE to EU. Therefore, the EU-based video transmission problem is formulated as(13)OP0:maxρg,n,pn,sn,rg,ne,rg,nb⁡ηEU=UtotalEtotalC1:∑n=1Npn≤PtotalC2:∑n=1NIl,npn≤Ithl,1≤l≤LC3:pn≥0,1≤n≤NC4:∑n=1Nρg,nsnrg,nb≥Rgb,1≤g≤GC5:∑g=1Gρg,n≤1,1≤n≤NC6:ρg,n∈0,1,1≤n≤N,1≤g≤GC7:sn∈0,1,1≤n≤N,where Ptotal is the total transmitted power by the secondary BS and Ithl is the interference threshold for PU l. In OP0, C2   denotes the total interference constraint from each PU, C4 requires that the BL data can be received by all the SUs in each group to guarantee the basic quality of reconstructed video, and C5 and C6 restrict that one subcarrier is exclusively used by only one group, since different groups may be interested in distinct video contents from the secondary BS. It is assumed that the secondary BS has the perfect knowledge of multicast grouping, the bandwidth of B, channel gains, and interference coefficients.


## 2.1. OFDM-Based CR System

The considered CR system is composed of a primary network and a CR network, which are both deployed in a cellular fashion. The primary network consists of a primary base station (BS) andL PUs, while the CR network is made up of a secondary BS and K SUs. The whole licensed spectrum spans B Hz, and each PU occupies a disjoint fraction of the spectrum, denoted as Bl. For the finer-grained spectrum utilization, the OFDM technology is adopted to divide the whole sensed spectrum B into N subcarriers, with each subcarrier spanning B0=B/N Hz [5, 9]. Figure 1 illustrates the details of CR system and spectrum distribution.Figure 1
System scenario and spectrum distribution.
(a)
CR system: a primary network and a CR network coexist
(b)
Spectrum division by OFDM and spectrum occupation of PUsSince the sensed spectrum by CR network is licensed to the primary network, the privilege of PUs to use the spectrum must be guaranteed, which typically necessitates the interference control. Therefore, a common interference evaluation model is introduced from [20].In [19], the power spectrum density (PSD) of subcarrier n is written as(1)φnf=pnTssin⁡πfTsπfTs2,where pn is the transmitted power of signal on subcarrier n and Ts is the duration of OFDM symbol. Then, the interference introduced by the signal on subcarrier n into PU l can be expressed as(2)I~l,nSP=GlSP∫dl,n-Bl/2dl,n+Bl/2φnfdf=Il,npn,where dl,n denotes the spectral distance between subcarrier n and the center frequency of PU l, GlSP denotes the channel gain from the secondary BS to PU l, and Il,n represents the interference caused by the normalized power on subcarrier n to PU l.Meanwhile, the interference caused by the primary BS to SUk on subcarrier n is calculated as [20](3)I~k,nPS=∑l=1LGk,nPS∫dl,n-Δf/2dl,n+Δf/2EIlωdω,where Gk,nPS denotes the channel gain from primary BS to SU k on subcarrier n, and E[·] is the expectation operator. E[Il(ω)] is the signal PSD of PU l after N-Fast-Fourier-transform (FFT) processing, and it is represented as(4)EIlω=12πN∫-ππΦlejϕsin⁡ω-ϕN/2sin⁡ω-ϕ/22dϕ,where Φl(ejϕ) is the PSD of the transmitted signal by the primary BS to PU l, which is generally modeled as an elliptically filtered white noise process [20].


## 2.2. LVM Transmission Model

As shown in Figure1(a), all K users are partitioned into G multicast groups according to the video contents they are interested in. The set and number of SUs in group g are denoted as Kg and Kg, respectively. Let Hk∣g,nSS denote the channel gain between the secondary BS and SU k in group g, and the corresponding channel gain to interference-plus-noise ratio (CINR) can be expressed as(5)γk∣g,n=Hk∣g,nSSσn2+I~k∣g,nPS,where σn2 is the noise power on subcarrier n and I~k∣g,nPS is calculated according to (3) with an index mapping from k∣g to k. Based on the Shannon formula, the achievable data rate of SU k in group g on subcarrier n is calculated by(6)rk∣g,n=B0log2⁡1+pnγk∣g,n.For CM, the lowest rate of SUs in a group is conservatively adopted to ensure the correct data reception of all SUs [15]. Thus, the data transmission rate of group g on subcarrier n is expressed as(7)rg,nc=mink∈Kg⁡rk∣g,n=B0log2⁡1+pnmink∈Kg⁡γk∣g,n.In wireless transmission, the heterogeneity of the receiving channels for different SUs will seriously limit the performance for CM. In order to overcome the shortcomings, LVM [18, 19] is introduced and modeled in this paper. In LVM, video data are transmitted resiliently on different subcarriers to adapt to the diverse channel conditions. This is accomplished by coding the source data into a BL and several ELs, and as long as the BL is received, SU can decode the video stream with the basic quality. If more ELs are received, the decoded video quality is increasingly improved [21]. The essential difference between CM and LVM lies in that the former requires that all SUs in a group receive the video with the identical quality, whilst the latter allows the differential reception, depending on the individual channel quality. Therefore, LVM provides a new degree of freedom, that is, the transmission rates on subcarriers, to exploit different channel conditions.Specifically, it is assumed that each videog, which is received by group g, is encoded into one BL with rate Rgb and one EL with rate Rge (see Figure 1(a)). Rgmax=Rgb+Rge denotes the maximal rate for video g. As in [19], we also assume that the fine grained scalability (FGS) coding technique is adopted, so that the EL can be truncated at any bit location with all the remaining bits still being useful at the decoder [22].In LVM, the BL data is of great importance to reconstruct the source video, and it is imperative that the BL can be received by all the SUs in a group. Hence, if subcarriern is used by group g for the BL transmission, the transmission rate is equal to that of CM; that is,(8)rg,nb=rg,nc.If subcarrier n is used by group g for the EL transmission at a rate of rg,ne, the receiving rate of SU k,k∈Kg on subcarrier n is written as(9)σrk∣g,n≥rg,nerg,ne,where σ(T) is 1 if T is true, otherwise 0. Define two indicators ρg,n and sn to perform the intersession and interlayer subcarrier assignment as follows. ρg,n∈{1,0} denotes whether subcarrier n is assigned to group g or not, and sn denotes whether subcarrier n is used for the BL transmission (sn=1) or the EL transmission (sn=0). The achieved rate of SU k,k∈Kg is calculated as(10)Rk=∑n=1Nρg,n1-snσrk∣g,n≥rg,nerg,ne+snrg,nb.


## 2.3. EU-Based Problem Formulation

In LVM, the achieved rate cannot reflect the quality directly, and some metrics, such as peak signal to noise ratio (PSNR) and mean square error (MSE), may be more accurate to evaluate the received video quality. For generalizing the expression, the utility functionUk(Rk) is defined to denote the change relationship of reconstructed video quality with the achieved rate. Moreover, Uk(Rk) is assumed to be nondecreasing and concave according to the measured or theoretical results [18, 19]. Then, the total weighted-utility is summed as(11)Utotal=∑k=1KwkUkRk,where wk(wk>0) is the weight of SU k and can be used to reflect the relative priorities among SUs.Apart from the video quality, the energy cost of video transmission over OFDM-based CR systems should be also considered. The total energy consumption in a timeslot includes three parts: the sense energyPs, that is, the energy consumed to sense the available spectrum, the transmission energy Pt, that is, the energy consumed to transmit data on the sensed spectrum, and the circuit energy Pc. Note that the concepts of “energy” and “power” are not strictly distinguished in this paper, since only the constant duration is multiplied or not. The total energy consumption is accumulated as(12)Etotal=Pt+Ps+Pc=∑n=1Npn+ξ∑n=1Npn+Pc,where ξ denotes the ratio of the sense energy Ps to the transmission energy Pt [15]. The formula accounts for the required energy consumption for downlink (DL) video multicast in OFDM-based CRNs. In practice, the total system power consumption includes the power consumption at both BS and end device sides. The power consumption of end devices during DL transmission mainly results from multiple active circuit modules, such as the channel estimation and feedback module and the baseband signal processing module. However, the aggregate value is relatively small compared with the power consumption of BS, whose radio frequency (RF) transmit power is often very large due to the large-scale pathloss fading. As a result, the power consumption of end devices is not considered in formula (12), in which the total system energy consumption is approximated as the power consumption at the secondary BS [7, 9–12].In the existing work, EE is maximized to optimize the achieved rate by the unit energy consumption [9–12]. For the video transmission, it makes more sense to shift EE to EU. Therefore, the EU-based video transmission problem is formulated as(13)OP0:maxρg,n,pn,sn,rg,ne,rg,nb⁡ηEU=UtotalEtotalC1:∑n=1Npn≤PtotalC2:∑n=1NIl,npn≤Ithl,1≤l≤LC3:pn≥0,1≤n≤NC4:∑n=1Nρg,nsnrg,nb≥Rgb,1≤g≤GC5:∑g=1Gρg,n≤1,1≤n≤NC6:ρg,n∈0,1,1≤n≤N,1≤g≤GC7:sn∈0,1,1≤n≤N,where Ptotal is the total transmitted power by the secondary BS and Ithl is the interference threshold for PU l. In OP0, C2   denotes the total interference constraint from each PU, C4 requires that the BL data can be received by all the SUs in each group to guarantee the basic quality of reconstructed video, and C5 and C6 restrict that one subcarrier is exclusively used by only one group, since different groups may be interested in distinct video contents from the secondary BS. It is assumed that the secondary BS has the perfect knowledge of multicast grouping, the bandwidth of B, channel gains, and interference coefficients.


## 3. Optimized Video Multicast Transmission

In order to achieve the optimal performance, spectrum assignmentρg,n, sn, power allocation pn, and the data transmission rate selection  rg,ne, rg,nb should be jointly determined. Nevertheless, the tight coupling among these variables will incur a prohibitively high computation complexity. Therefore, allowing for the features of LVM and CRN, we decompose OP0 into three steps, that is, BL subcarrier assignment, EL subcarrier assignment, and power allocation. Additionally, how to deal with the spectrum scarcity is also discussed.


### 3.1. BL Subcarrier Assignment

For unifying the expression, letI0,n=1, 1≤n≤N, Ith0=Ptotal. In the traditional OFDM systems, where no primary network exists, equal power Ptotal/N is generally assumed for subcarrier assignment [23]. In OFDM-based CRNs, the interference to PUs has to be considered, and the interference coefficient Il,n varies from one subcarrier to another. The coarse assumption of equal power deviates from the actual situations and may degrade the performance. Instead, the ladder-profile power is assumed in this paper:(14)pn=min⁡Ith0NI0,n,Ith1NI1,n,…,IthLNIL,n.This idea is borrowed from [24], which supposes that each subcarrier produces the identical amount of interference to PUs. Obviously, the power pn in (14) rests with the interference coefficients and may differ among subcarrier individuals. It can be verified that all the power and interference constraints C1 and C2 are satisfied by (14); that is, ∑n=1NIl,npn≤∑n=1NIl,nIthl/NIl,n≤Ithl for 0≤l≤L. Furthermore, a factor β≤1 is calculated to scale the power pn as(15)β=max0≤l≤L⁡∑n=1NIl,npnIthl,pn=pnβ,such that at least one constraint in C1, C2 is active for maintaining the reasonable power level.With (14) and (15), the achievable rate rg,nb of BL is readily computed according to (7) and (8). The details of BL subcarrier assignment are provided in Algorithm 1. The objective of Algorithm 1 is to meet the BL rate requirements with as small number of subcarriers as possible. At each iteration, the most unsatisfied group g with the largest gap from the BL rate requirement is selected, and then the best subcarrier n with the largest achievable rate is assigned to this group. Subsequently, the state of whether the rate requirement is fulfilled or not is inspected for group g, and if satisfied, it is immediately removed from the set of unsatisfied groups G to save more subcarriers for EL data transmission.Algorithm 1:Algorithm for BL subcarrier assignment.
Input:
(1)Ptotal, γm∣g,n for m∈Kg, 1≤g≤G, pn for 1≤n≤N, Il,n for 1≤l≤L, 1≤n≤N, Ithl for 1≤l≤L, and Rgb for 1≤g≤G;
Initialization:
(2)  Initialize the set of BL-rate-unsatisfied groupsG={1,2,…,G}, the set of unassigned subcarriers A={1,2,…,N},
the achieved rate  R~g=0 for group g, the set of assigned subcarriers Ωk=⌀ for SU k, the assignment indicator
ρg,n=0, 1≤n≤N, 1≤g≤G;
Iteration:
(3)whileA≠⌀andG≠⌀do
(4)  Find groupg∈G satisfying R~g-Rgb≤R~i-Rib, ∀i∈G;
(5)  Select subcarriern to maximize the achieved rate as n=arg⁡maxj⁡rg,jb, j∈A;
(6)  Assign subcarriern to group g for BL data transmission as ρg,n=1, A=A-{n}, sn=1, and update R~g, Ωk as
R~g=R~g+rg,nb, Ωk=Ωk∪{n}, ∀k∈Kg;
(7)ifR~g≥Rgbthen
(8)   Delete groupg as G=G-{g};
(9)end if
(10)end while
Output:
(11)G,A,R~g,Ωk,ρg,n,sn.EitherA=⌀ or G=⌀ will terminate Algorithm 1. However, only when A≠⌀, that is, there is a surplus of subcarriers, EL subcarrier assignment is needed, which progressively improves the quality of video.


### 3.2. EL Subcarrier Assignment

For EL subcarrier assignment, besides the task of assigning subcarriers to the appropriate groups, the data transmission raterg,ne should also be determined jointly. In OP0, rg,ne is limited to be a nonnegative real number. It seems that the derivative method has to be adopted to obtain the optimal rg,ne as the value of rg,ne is continuous. However, by a close observation, rg,ne can only take finite K+1 values. Specifically, the following theorem is presented to reduce the complexity of subcarrier assignment for EL.Theorem 1.
With the nondecreasing utility functions, for any subcarriern assigned to transmit the EL data, the optimal rate rg,ne is equal to 0 or rg,ne∈Bn=rm∣g,n∣1≤g≤G,m∈Kg; that is, rg,ne must be zero or one of the K achievable rates in Bn.Proof.
See AppendixA.Theorem1 indicates that if subcarrier n is assigned to group g, then only the Kg achievable rates may be the optimal solution. In total, all K=∑g=1GKg achievable rates constitute the candidates for the optimal solution if subcarrier n is not assigned yet. This contracts the domain of rg,ne from the nonnegative real number to K+1 discrete values and thereby eases the determination process of rg,ne. Define SU m as the barrier user (BU) on subcarrier n if its achievable rate is selected for data transmission, that is, rg,ne=rm∣g,n, since its rate marks the watershed between the success and failure in data reception. According to Theorem 1, the rate of SU k in (10) is rewritten as(16)Rk=∑n=1Nρg,n1-snσrk∣g,n≥rm∣g,nrm∣g,n+snrg,nb,where BU m may be different among subcarriers and needs to be decided by Algorithm 2.Algorithm 2:Algorithm for EL subcarrier assignment.
Input:
(1)wk,ξ,Pc,Ptotal, γm∣g,n for m∈Kg, 1≤g≤G, n∈A, pn for n∈A, Il,n for 1≤l≤L, n∈A and Ithl for 1≤l≤L;
Initialization:
(2) Initialize the set of unassigned subcarriersA, the set of assigned subcarriers Ωk for SU k, and the assignment indicator
ρg,n,n∉A,1≤g≤G according to the output of Algorithm 1. Set ρg,n=0, n∈A, 1≤g≤G, and the achieved rate
Rk=R~g for SU k, k∈Kg;
Iteration:
(3)whileA≠⌀do
(4)  Select subcarriern,n∈A in order, and find group g and BU m to maximize the achieved EU as
g,m=arg⁡maxρi,n,rj∣i,n⁡ηEU, 1≤i≤G, j∈Kg;
(5)   Assign subcarriern to group g for EL data transmission as ρg,n=1, A=A-{n}, sn=0, and update Rk,Ωk as
Rk=Rk+σ(rk∣g,n≥rm∣g,n)rm∣g,n for k∈Kg,Ωk=Ωk∪{n} for rk∣g,n≥rm∣g,n, k∈Kg;
(6)end while
Output:
(7)Ωk,ρg,n,sn,rm∣g,n.The objective of Algorithm2 is to maximize the EU via subcarrier assignment for EL. At each iteration, subcarrier n is selected in order, and then, based on Theorem 1, the traversal of K+1 values for rg,ne is carried out in order to achieve the maximal EU ηEU. In Algorithm 2, for 1≤i≤G, j∈Kg, to maximize ηEU is equivalent to maximizing the total utility as(17)arg⁡maxρi,n,rj∣i,n⁡ηEU=arg⁡maxρi,n,rj∣i,n⁡∑k=1KwkUkRk,since the power consumption values are equal for the denominator of ηEU when determining which group subcarrier n is assigned to and which SU is selected as the BU.


### 3.3. EU-Based Power Allocation

With the determined subcarrier assignment, the rate of SUk,k∈Kg in (16) can be further expressed as(18)Rk=∑n∈Ωk1-snrm∣g,n+snrg,nb=∑n∈Ωkrn=∑n∈ΩkB0log2⁡1+pnγn,where rn and γn are the achievable rate and CINR, respectively, for the BU related to subcarrier n in group g. Note that if subcarrier n is used for EL, the BU is specified in Algorithm 2. Otherwise, for BL, the BU is the SU in group g with the worst channel condition on subcarrier n.When determining subcarrier assignment, the ladder-profile power is simply assumed in (14) and (15) to evaluate the rate or utility in Algorithms 1 and 2. Thus, the power needs to be reallocated to maximize the EU for multicast video transmission, which is formulated as(19)OP1:maxpn⁡ηEUp=UtotalpEtotalpC1:∑n=1Npn≤PtotalC2:∑n=1NIl,npn≤Ithl,1≤l≤LC3:pn≥0,1≤n≤NC4:∑n∈Θgbrn≥Rgb,1≤g≤G,where Θgb is the set of subcarriers assigned to group g for the BL transmission and p denotes the vector of pn, 1≤n≤N. Even though the optimization variables associated with subcarrier assignment in OP0 have been settled, OP1 is still intractable due to the fractional objective function as well as a large number of constraints on the power and BL rate.As in our previous work [9],fractional programming [25] is employed to deal with the objective function issue. With the positive parameter α, a new function g(p,α) is defined as(20)gp,α=Etotalp-αUtotalp.Then, another problem is formulated as(21)OP2:minpn⁡gp,αC1,C2,C3,C4.LetS denote the feasible region of OP1 and OP2, and define the optimal value and solution of OP2 as(22)Fα=minp⁡gp,α∣p∈S,fα=arg⁡minp⁡gp,α∣p∈S.The following lemma can relate OP1 and OP2, and the detailed proof of Lemma 2 can be found in [25].Lemma 2.
p
∗ and 1/α∗ correspond to the optimal solution and value of problem OP1 in (19); that is, 1/α∗=Utotal(p∗)/Etotal(p∗)=maxpηEUp∣p∈S, if and only if(23)Fα∗=0,fα∗=p∗.Lemma2 indicates that the optimal solution to OP2 is also the optimal solution to OP1, provided that (23) is satisfied. Hence, solving OP1 can be realized by finding the optimal power allocation of OP2 for a given α, and then update α until (23) is established. The problem at hand is how to achieve the optimal power allocation for OP2 with a given α. To this end, the following theorem is presented.Theorem 3.
If the utility function is nondecreasing and concave for each SU, problemOP2 with a fixed α belongs to convex optimization ones.Proof.
See AppendixB.The convexity ofOP2 with a fixed α enables the solution in the dual domain without the dual gap [26]. Namely, the optimal solution to the dual problem is exactly the optimal solution to the primal problem OP2. The Lagrangian of OP2 is defined as(24)L=Etotalp-αUtotalp+∑l=0Lλl∑n=1NIl,npn-Ithl+∑g=1GυgRgb-∑n∈Θgbrn,where λ=[λ0,λ1,…,λL]⪰0 and υ=[υ1,…,υG]⪰0 are vectors of dual variables. Therefore, its dual function is defined as(25)hλ,υ,α=minpn⁡L.According to the Karush-Kuhn-Tucker (KKT) conditions [26], the optimal power allocation should satisfy(26)∂L∂pn=-α∑k∈ΦnwkdUkRkdRk-snυgB0ln⁡2γn1+pnγn+∑l=0LλlIl,n+ξ+1=0,where Φn denotes the set of SUs who can receive data from subcarrier n. By some manipulations, the optimal pn is derived as(27)pn=B0ln⁡2α∑k∈ΦnwkdUkRk/dRk+snυg∑l=0LλlIl,n+ξ+1-1γn+,where x+ denotes max⁡{0,x}. In general, the analytical expression of pn is unavailable, except that Uk(Rk) is linear with Rk for any SU k. In such cases, dUk(Rk)/dRk is independent of pn. Otherwise, n nonlinear equations with n unknowns in (26) have to be solved numerically to find the optimal pn.Once the optimalpn is achieved, the dual function h(λ,υ,α) in (25) can be computed for the fixed tuple (λ,υ). Now, the dual problem needs to be solved as(28)maxλ⪰0,υ⪰0⁡hλ,υ,α.The optimal solution to (28) can be steadily obtained in a subgradient method [27], which iteratively updates the dual variables (λ,υ) in the subgradient direction until they converge. Specifically, the updates for (28) can be performed as(29)λlt+1=λlt-εtIthl-∑n=1NIl,npn+,0≤l≤L,υgt+1=υgt-ωt∑n∈Θgbrn-Rgb+,1≤g≤G,where εt,ωt>0 form two sequences of step sizes. In general, the step size εt,ωt can be selected as κ/t or κ/t based on the diminishing rules [27], where κ>0 is a constant, and t is the number of iterations. The subgradient method can guarantee that (λ,υ) converges to the optimal point as long as εt and ωt are sufficiently small [27]. Due to the convexity depicted by Theorem 3, the optimal solution to OP2 can be achieved according to (26) or (27) once the dual optimal solution (λ∗,υ∗) is reached.With the optimal power allocationp for the given α, update α=Etotal(p)/Utotal(p) and solve problem OP2 again. The process is repeated until the condition of Lemma 2 is satisfied. The iteration convergence of α update is proved in [24]. The details of power allocation are provided in Algorithm 3, where x denotes the 2-norm of vector x.Algorithm 3:EU-based power allocation.
Input:
(1)wk,ξ,Pc,Ptotal, γm∣g,n for m∈Kg, 1≤g≤G, 1≤n≤N, Il,n for 1≤l≤L, 1≤n≤N,
Ithl for 1≤l≤L, Ωk for 1≤k≤K, Φn, sn for 1≤n≤N and Θgb,Rgb for 1≤g≤G;
Initialization:
(2)  Initializeα, and the tolerable error ϵ. Set F(α)=∞;
Iteration:
(3)whileF(α)>ϵdo
(4)  Initialize the tolerable errorϵ~, the dual variables λl, 0≤l≤L, υg, 1≤g≤G, and e=∞;
(5)whilee>ϵ~do
(6)   Calculatepn according to (26) or (27);
(7)   Updateλl,υg according to (29);
(8) Compute e as e=Δλ/λ(t)+Δυ/υ(t) with Δλ=λ(t+1)-λ(t), Δυ=υ(t+1)-υ(t);
(9)end while
(10)Compute F(α), and update α=Etotal(p)/Utotal(p);
(11)end while
Output:
(12)pn and the final EU ηEU=1/α.The objective of Algorithm3 is to maximize the EU by allocating power among subcarriers, whilst guaranteeing the BL rate requirements. In the inner loop, power allocation p and dual variables (λ,υ) are alternately updated to attain the optimal solution for a fixed α. In the outer loop, α is repeatedly updated until it converges. It can be asserted that Algorithm 3 can reap the optimal power for OP1, owing to the factor that both the inner-loop and outer-loop iterations can converge to the respective optimal solutions.In a practical system, the video is encoded into data layers at first, and then the BL data is mapped onto subcarriers by Algorithm1. Sequentially, the EL data is arranged on the unused subcarriers by Algorithm 2, and the power, corresponding to the transmission rate, for each subcarrier is ultimately settled by Algorithm 3. Through all three algorithms, the EU is consecutively optimized to deliver the video to diverse receivers with the detected spectrum and limited energy.


### 3.4. Spectrum Scarcity Discussion

Typically, the video service is “bandwidth-hungry,” and wireless network states, including channel quality, service requests, and the numbers of SUs and PUs, remain ever-changing. Sometimes the sensed spectrum is not enough for supporting the BL rates of all the groups. For example, the available bandwidth is overly narrow due to the full occupation of the primary network, or the interference threshold is excessively tight because the CR network is geographically close to the primary network. For such cases, problem OP0 has no feasible solutions, which is usually reflected by the phenomenon that Algorithm 1 stops with A=⌀.In order to address the issue, theaggressive solution is to drop some SUs according to the service emergency or the requested rate of BL data, and in the next timeslot, more transmission opportunities are offered to these sacrificed SUs. Theconservative solution is to relax the BL rate constraint C4 in OP0 and adaptively determine the transmission rates of each SU by its contribution to the EU ηEU; that is, the SU who can increase the total utility most by the unit energy consumption will achieve the transmission opportunity. This solution can restrain parts of SUs from data transmission in an “EU-competitive” fashion.


## 3.1. BL Subcarrier Assignment

For unifying the expression, letI0,n=1, 1≤n≤N, Ith0=Ptotal. In the traditional OFDM systems, where no primary network exists, equal power Ptotal/N is generally assumed for subcarrier assignment [23]. In OFDM-based CRNs, the interference to PUs has to be considered, and the interference coefficient Il,n varies from one subcarrier to another. The coarse assumption of equal power deviates from the actual situations and may degrade the performance. Instead, the ladder-profile power is assumed in this paper:(14)pn=min⁡Ith0NI0,n,Ith1NI1,n,…,IthLNIL,n.This idea is borrowed from [24], which supposes that each subcarrier produces the identical amount of interference to PUs. Obviously, the power pn in (14) rests with the interference coefficients and may differ among subcarrier individuals. It can be verified that all the power and interference constraints C1 and C2 are satisfied by (14); that is, ∑n=1NIl,npn≤∑n=1NIl,nIthl/NIl,n≤Ithl for 0≤l≤L. Furthermore, a factor β≤1 is calculated to scale the power pn as(15)β=max0≤l≤L⁡∑n=1NIl,npnIthl,pn=pnβ,such that at least one constraint in C1, C2 is active for maintaining the reasonable power level.With (14) and (15), the achievable rate rg,nb of BL is readily computed according to (7) and (8). The details of BL subcarrier assignment are provided in Algorithm 1. The objective of Algorithm 1 is to meet the BL rate requirements with as small number of subcarriers as possible. At each iteration, the most unsatisfied group g with the largest gap from the BL rate requirement is selected, and then the best subcarrier n with the largest achievable rate is assigned to this group. Subsequently, the state of whether the rate requirement is fulfilled or not is inspected for group g, and if satisfied, it is immediately removed from the set of unsatisfied groups G to save more subcarriers for EL data transmission.Algorithm 1:Algorithm for BL subcarrier assignment.
Input:
(1)Ptotal, γm∣g,n for m∈Kg, 1≤g≤G, pn for 1≤n≤N, Il,n for 1≤l≤L, 1≤n≤N, Ithl for 1≤l≤L, and Rgb for 1≤g≤G;
Initialization:
(2)  Initialize the set of BL-rate-unsatisfied groupsG={1,2,…,G}, the set of unassigned subcarriers A={1,2,…,N},
the achieved rate  R~g=0 for group g, the set of assigned subcarriers Ωk=⌀ for SU k, the assignment indicator
ρg,n=0, 1≤n≤N, 1≤g≤G;
Iteration:
(3)whileA≠⌀andG≠⌀do
(4)  Find groupg∈G satisfying R~g-Rgb≤R~i-Rib, ∀i∈G;
(5)  Select subcarriern to maximize the achieved rate as n=arg⁡maxj⁡rg,jb, j∈A;
(6)  Assign subcarriern to group g for BL data transmission as ρg,n=1, A=A-{n}, sn=1, and update R~g, Ωk as
R~g=R~g+rg,nb, Ωk=Ωk∪{n}, ∀k∈Kg;
(7)ifR~g≥Rgbthen
(8)   Delete groupg as G=G-{g};
(9)end if
(10)end while
Output:
(11)G,A,R~g,Ωk,ρg,n,sn.EitherA=⌀ or G=⌀ will terminate Algorithm 1. However, only when A≠⌀, that is, there is a surplus of subcarriers, EL subcarrier assignment is needed, which progressively improves the quality of video.


## 3.2. EL Subcarrier Assignment

For EL subcarrier assignment, besides the task of assigning subcarriers to the appropriate groups, the data transmission raterg,ne should also be determined jointly. In OP0, rg,ne is limited to be a nonnegative real number. It seems that the derivative method has to be adopted to obtain the optimal rg,ne as the value of rg,ne is continuous. However, by a close observation, rg,ne can only take finite K+1 values. Specifically, the following theorem is presented to reduce the complexity of subcarrier assignment for EL.Theorem 1.
With the nondecreasing utility functions, for any subcarriern assigned to transmit the EL data, the optimal rate rg,ne is equal to 0 or rg,ne∈Bn=rm∣g,n∣1≤g≤G,m∈Kg; that is, rg,ne must be zero or one of the K achievable rates in Bn.Proof.
See AppendixA.Theorem1 indicates that if subcarrier n is assigned to group g, then only the Kg achievable rates may be the optimal solution. In total, all K=∑g=1GKg achievable rates constitute the candidates for the optimal solution if subcarrier n is not assigned yet. This contracts the domain of rg,ne from the nonnegative real number to K+1 discrete values and thereby eases the determination process of rg,ne. Define SU m as the barrier user (BU) on subcarrier n if its achievable rate is selected for data transmission, that is, rg,ne=rm∣g,n, since its rate marks the watershed between the success and failure in data reception. According to Theorem 1, the rate of SU k in (10) is rewritten as(16)Rk=∑n=1Nρg,n1-snσrk∣g,n≥rm∣g,nrm∣g,n+snrg,nb,where BU m may be different among subcarriers and needs to be decided by Algorithm 2.Algorithm 2:Algorithm for EL subcarrier assignment.
Input:
(1)wk,ξ,Pc,Ptotal, γm∣g,n for m∈Kg, 1≤g≤G, n∈A, pn for n∈A, Il,n for 1≤l≤L, n∈A and Ithl for 1≤l≤L;
Initialization:
(2) Initialize the set of unassigned subcarriersA, the set of assigned subcarriers Ωk for SU k, and the assignment indicator
ρg,n,n∉A,1≤g≤G according to the output of Algorithm 1. Set ρg,n=0, n∈A, 1≤g≤G, and the achieved rate
Rk=R~g for SU k, k∈Kg;
Iteration:
(3)whileA≠⌀do
(4)  Select subcarriern,n∈A in order, and find group g and BU m to maximize the achieved EU as
g,m=arg⁡maxρi,n,rj∣i,n⁡ηEU, 1≤i≤G, j∈Kg;
(5)   Assign subcarriern to group g for EL data transmission as ρg,n=1, A=A-{n}, sn=0, and update Rk,Ωk as
Rk=Rk+σ(rk∣g,n≥rm∣g,n)rm∣g,n for k∈Kg,Ωk=Ωk∪{n} for rk∣g,n≥rm∣g,n, k∈Kg;
(6)end while
Output:
(7)Ωk,ρg,n,sn,rm∣g,n.The objective of Algorithm2 is to maximize the EU via subcarrier assignment for EL. At each iteration, subcarrier n is selected in order, and then, based on Theorem 1, the traversal of K+1 values for rg,ne is carried out in order to achieve the maximal EU ηEU. In Algorithm 2, for 1≤i≤G, j∈Kg, to maximize ηEU is equivalent to maximizing the total utility as(17)arg⁡maxρi,n,rj∣i,n⁡ηEU=arg⁡maxρi,n,rj∣i,n⁡∑k=1KwkUkRk,since the power consumption values are equal for the denominator of ηEU when determining which group subcarrier n is assigned to and which SU is selected as the BU.


## 3.3. EU-Based Power Allocation

With the determined subcarrier assignment, the rate of SUk,k∈Kg in (16) can be further expressed as(18)Rk=∑n∈Ωk1-snrm∣g,n+snrg,nb=∑n∈Ωkrn=∑n∈ΩkB0log2⁡1+pnγn,where rn and γn are the achievable rate and CINR, respectively, for the BU related to subcarrier n in group g. Note that if subcarrier n is used for EL, the BU is specified in Algorithm 2. Otherwise, for BL, the BU is the SU in group g with the worst channel condition on subcarrier n.When determining subcarrier assignment, the ladder-profile power is simply assumed in (14) and (15) to evaluate the rate or utility in Algorithms 1 and 2. Thus, the power needs to be reallocated to maximize the EU for multicast video transmission, which is formulated as(19)OP1:maxpn⁡ηEUp=UtotalpEtotalpC1:∑n=1Npn≤PtotalC2:∑n=1NIl,npn≤Ithl,1≤l≤LC3:pn≥0,1≤n≤NC4:∑n∈Θgbrn≥Rgb,1≤g≤G,where Θgb is the set of subcarriers assigned to group g for the BL transmission and p denotes the vector of pn, 1≤n≤N. Even though the optimization variables associated with subcarrier assignment in OP0 have been settled, OP1 is still intractable due to the fractional objective function as well as a large number of constraints on the power and BL rate.As in our previous work [9],fractional programming [25] is employed to deal with the objective function issue. With the positive parameter α, a new function g(p,α) is defined as(20)gp,α=Etotalp-αUtotalp.Then, another problem is formulated as(21)OP2:minpn⁡gp,αC1,C2,C3,C4.LetS denote the feasible region of OP1 and OP2, and define the optimal value and solution of OP2 as(22)Fα=minp⁡gp,α∣p∈S,fα=arg⁡minp⁡gp,α∣p∈S.The following lemma can relate OP1 and OP2, and the detailed proof of Lemma 2 can be found in [25].Lemma 2.
p
∗ and 1/α∗ correspond to the optimal solution and value of problem OP1 in (19); that is, 1/α∗=Utotal(p∗)/Etotal(p∗)=maxpηEUp∣p∈S, if and only if(23)Fα∗=0,fα∗=p∗.Lemma2 indicates that the optimal solution to OP2 is also the optimal solution to OP1, provided that (23) is satisfied. Hence, solving OP1 can be realized by finding the optimal power allocation of OP2 for a given α, and then update α until (23) is established. The problem at hand is how to achieve the optimal power allocation for OP2 with a given α. To this end, the following theorem is presented.Theorem 3.
If the utility function is nondecreasing and concave for each SU, problemOP2 with a fixed α belongs to convex optimization ones.Proof.
See AppendixB.The convexity ofOP2 with a fixed α enables the solution in the dual domain without the dual gap [26]. Namely, the optimal solution to the dual problem is exactly the optimal solution to the primal problem OP2. The Lagrangian of OP2 is defined as(24)L=Etotalp-αUtotalp+∑l=0Lλl∑n=1NIl,npn-Ithl+∑g=1GυgRgb-∑n∈Θgbrn,where λ=[λ0,λ1,…,λL]⪰0 and υ=[υ1,…,υG]⪰0 are vectors of dual variables. Therefore, its dual function is defined as(25)hλ,υ,α=minpn⁡L.According to the Karush-Kuhn-Tucker (KKT) conditions [26], the optimal power allocation should satisfy(26)∂L∂pn=-α∑k∈ΦnwkdUkRkdRk-snυgB0ln⁡2γn1+pnγn+∑l=0LλlIl,n+ξ+1=0,where Φn denotes the set of SUs who can receive data from subcarrier n. By some manipulations, the optimal pn is derived as(27)pn=B0ln⁡2α∑k∈ΦnwkdUkRk/dRk+snυg∑l=0LλlIl,n+ξ+1-1γn+,where x+ denotes max⁡{0,x}. In general, the analytical expression of pn is unavailable, except that Uk(Rk) is linear with Rk for any SU k. In such cases, dUk(Rk)/dRk is independent of pn. Otherwise, n nonlinear equations with n unknowns in (26) have to be solved numerically to find the optimal pn.Once the optimalpn is achieved, the dual function h(λ,υ,α) in (25) can be computed for the fixed tuple (λ,υ). Now, the dual problem needs to be solved as(28)maxλ⪰0,υ⪰0⁡hλ,υ,α.The optimal solution to (28) can be steadily obtained in a subgradient method [27], which iteratively updates the dual variables (λ,υ) in the subgradient direction until they converge. Specifically, the updates for (28) can be performed as(29)λlt+1=λlt-εtIthl-∑n=1NIl,npn+,0≤l≤L,υgt+1=υgt-ωt∑n∈Θgbrn-Rgb+,1≤g≤G,where εt,ωt>0 form two sequences of step sizes. In general, the step size εt,ωt can be selected as κ/t or κ/t based on the diminishing rules [27], where κ>0 is a constant, and t is the number of iterations. The subgradient method can guarantee that (λ,υ) converges to the optimal point as long as εt and ωt are sufficiently small [27]. Due to the convexity depicted by Theorem 3, the optimal solution to OP2 can be achieved according to (26) or (27) once the dual optimal solution (λ∗,υ∗) is reached.With the optimal power allocationp for the given α, update α=Etotal(p)/Utotal(p) and solve problem OP2 again. The process is repeated until the condition of Lemma 2 is satisfied. The iteration convergence of α update is proved in [24]. The details of power allocation are provided in Algorithm 3, where x denotes the 2-norm of vector x.Algorithm 3:EU-based power allocation.
Input:
(1)wk,ξ,Pc,Ptotal, γm∣g,n for m∈Kg, 1≤g≤G, 1≤n≤N, Il,n for 1≤l≤L, 1≤n≤N,
Ithl for 1≤l≤L, Ωk for 1≤k≤K, Φn, sn for 1≤n≤N and Θgb,Rgb for 1≤g≤G;
Initialization:
(2)  Initializeα, and the tolerable error ϵ. Set F(α)=∞;
Iteration:
(3)whileF(α)>ϵdo
(4)  Initialize the tolerable errorϵ~, the dual variables λl, 0≤l≤L, υg, 1≤g≤G, and e=∞;
(5)whilee>ϵ~do
(6)   Calculatepn according to (26) or (27);
(7)   Updateλl,υg according to (29);
(8) Compute e as e=Δλ/λ(t)+Δυ/υ(t) with Δλ=λ(t+1)-λ(t), Δυ=υ(t+1)-υ(t);
(9)end while
(10)Compute F(α), and update α=Etotal(p)/Utotal(p);
(11)end while
Output:
(12)pn and the final EU ηEU=1/α.The objective of Algorithm3 is to maximize the EU by allocating power among subcarriers, whilst guaranteeing the BL rate requirements. In the inner loop, power allocation p and dual variables (λ,υ) are alternately updated to attain the optimal solution for a fixed α. In the outer loop, α is repeatedly updated until it converges. It can be asserted that Algorithm 3 can reap the optimal power for OP1, owing to the factor that both the inner-loop and outer-loop iterations can converge to the respective optimal solutions.In a practical system, the video is encoded into data layers at first, and then the BL data is mapped onto subcarriers by Algorithm1. Sequentially, the EL data is arranged on the unused subcarriers by Algorithm 2, and the power, corresponding to the transmission rate, for each subcarrier is ultimately settled by Algorithm 3. Through all three algorithms, the EU is consecutively optimized to deliver the video to diverse receivers with the detected spectrum and limited energy.


## 3.4. Spectrum Scarcity Discussion

Typically, the video service is “bandwidth-hungry,” and wireless network states, including channel quality, service requests, and the numbers of SUs and PUs, remain ever-changing. Sometimes the sensed spectrum is not enough for supporting the BL rates of all the groups. For example, the available bandwidth is overly narrow due to the full occupation of the primary network, or the interference threshold is excessively tight because the CR network is geographically close to the primary network. For such cases, problem OP0 has no feasible solutions, which is usually reflected by the phenomenon that Algorithm 1 stops with A=⌀.In order to address the issue, theaggressive solution is to drop some SUs according to the service emergency or the requested rate of BL data, and in the next timeslot, more transmission opportunities are offered to these sacrificed SUs. Theconservative solution is to relax the BL rate constraint C4 in OP0 and adaptively determine the transmission rates of each SU by its contribution to the EU ηEU; that is, the SU who can increase the total utility most by the unit energy consumption will achieve the transmission opportunity. This solution can restrain parts of SUs from data transmission in an “EU-competitive” fashion.


## 4. Simulation Results

In this section, channel fading gains from the secondary BS to SUs and PUs are modeled as consisting of six independent Rayleigh multipaths with an average channel power gain of 0 dB. The sensed spectrum is divided intoN=64 subcarriers with the subcarrier spacing B0=15 KHz in frequency band. L=4 PUs are assumed to be active in the primary network, and each PU occupies a fraction of spectrum with an equal bandwidth of 2B0. The locations of PU bands are assumed to evenly span the whole frequency band.The noise power plus the interference power is identically set asσ02 for each subcarrier, and the total transmission power is Ptotal=1 W. The average signal to interference-plus-noise ratio (SINR), which is defined as Ptotal/(Nσ02), ranges from 8dB to 28dB. The interference threshold of each PU Ithl=0.5mW for 1≤l≤L. The weight wk=1, 1≤k≤K, is designated to maximize the system sum-utility. The circuit energy consumption Pc=0.1 W, and the ratio of the sense energy to the transmission energy is set as ξ=0.1.In the CR network, two video sequences, that is,Foreman andHarbour in [18], are dedicated to G=2 groups. The piecewise linear utility function [18] is selected to evaluate the PSNR performance as(30)Ukr=Ug0+χgr-Rgbr≤RgbUg0+τgr-RgbRgb<r≤Rgmax⁡Ug0+τgRgmax⁡-Rgbr>Rgmax⁡,where Ug0, k∈Kg, denotes the PSNR when only the BL data is received and χg and τg are the slopes of the two pieces. The parameters Ug0, Rgb, Rgmax, χg, and τg are set according to the video sequencesForeman andHarbour in [18]. With the specific utility function (30), the optimal pn in (27) can be rewritten as(31)pn=B0ln⁡2α∑k∈ΦnLwkχg+∑k∈ΦnMwkτg+snυg∑l=0LλlIl,n+ξ+1-1γn+,where ΦnL and ΦnM are the subsets of Φn and denote the sets of SUs in Φn, whose total rate Rk, respectively, satisfies Rk≤Rgb and Rgb<Rk≤Rgmax. For comparison, another utility function Uk(r)=r is also simulated to maximize the weighted-rate. All simulation results are averaged over 1000 channel realizations.Figure2 shows the convergence process of dual variables λl, l=0,1,2,3,4, and υg, g=1,2, in the EU-based power allocation algorithm, where the number of SUs is 4 in each group, and the average SINR is equal to 12dB. It can be observed that the dual variables rapidly converge to the stable points within 10~80 iterations. Table 1 lists the average iteration numbers of auxiliary variable α. No matter which objective function is optimized, weighted-PSNR or weighted-rate, approximately 2~4 iterations are sufficient to achieve the converged α, which indicates that the high convergence rate of α is little impacted by the objective function. Combining Figure 2 with Table 1, it is concluded that the EU-based power allocation algorithm possesses the favourable iteration complexity.Table 1
Average number of iterations for auxiliary variableα.
Optimization objective function
Weighted-PSNR
Weighted-rate
Average number of iterations
2.1
3.6Figure 2
The convergence process of dual variables in the energy-utility-based power allocation algorithm.In Figure3, the EEs of the CM [14] and LVM are demonstrated with the increase of the SU number in each group. That is, the utility function is chosen as Uk(r)=r. As the number of SUs increases, both multicast methods are capable of achieving the higher EE, due to the fact that more SUs can benefit from a single transmission. In the meantime, the EE gap between them continuously enlarges, because the LVM enables the differential data reception and breaks through the worst channel restraint in CM. For example, with SINR = 12 dB, when the number of SUs is 10 in each group, the EE gain is nearly up to 40%.Figure 3
The energy efficiency comparison between the conventional multicast and layered video multicast.In Figure4, we compare three optimization objective functions, that is, the EU (the weighted-PSNR divided by the total energy consumption), the EE (the weighted-rate divided by the total energy consumption) [9–12], and the weighted-PSNR [18, 19], for the video transmission in terms of the achieved PSNR by the unit energy consumption. From Figure 4, it can be seen that because of integrating the energy cost into consideration the EE optimization achieves a significant performance gain over the weighted-PSNR optimization. Moreover, the EU optimization is considerably superior to the EE optimization thanks to the direct PSNR metric for the video quality. For example, with SINR = 15 dB, when the number of SUs is 4 in each group, the EU optimization can achieve about 120% and 20% performance gain compared to the weighted-PSNR and EE optimizations, respectively. It means that the selection of EU as the optimization objective is much qualified for green video transmission over OFDM-based CR systems.Figure 4
The energy utility comparison among different optimization objective functions.The proportion distribution is shown in Figure5 for the SUs, who are classified into three categories according to the total received rate; that is, Rk<Rgb, Rgb≤Rk<Rgb+Rge, and Rk=Rgb+Rge, k∈Kg. As the average SINR rises, the proportion of SUs, whose rates are lower than the BL rate, decreases evidently, while the proportion of SUs, who achieve the saturated rate Rgb+Rge, increases remarkably. It indicates that the “bandwidth-hungry” video services can fully utilize the sensed spectrum by our proposed algorithms. In addition, the proportion (the middle green bars) of the SUs, whose rates satisfy Rgb≤Rk<Rgb+Rge, always dominates the proportion distribution except when the average SINR equals 10 dB, which reveals that the relative fairness among SUs can be ensured as a result of the BL rate requirements and the selected PSNR utility function.Figure 5
The proportion distribution of SUs within three rate ranges.


## 5. Conclusion

This paper studies LVM transmission over OFDM-based CR systems, where multiple interference constraints are necessitated to carry out the performance protection for PUs. An EU-based optimization model, which is well tailored for green video delivery, is formulated, and then the BL and EL spectrum assignments are separately proposed for the complexity reduction. Furthermore, the EU-based power allocation algorithm is also proposed by jointly employing fractional programming and subgradient method. Simulation results show that the proposed power allocation algorithm can converge to the optimal solution rapidly, and LVM greatly outperforms CM, which is attributed to allowing the elastic data reception. Via the proposed algorithms, the sensed spectrum can be fully exploited by “bandwidth-hungry” video services, and the EU-based optimization notably surpasses the EE-based optimization in considering the video quality and energy consumption simultaneously. As the future work, multi-EL modeling and the sequential receiving issue for the video transmission will be considered to further explore the video multicast potential over OFDM-based CR systems.



---

*Source: 138328-2015-10-12.xml*