 Research
 Open Access
 Published:
Optimal linear precoding for opportunistic spectrum sharing under arbitrary input distributions assumption
EURASIP Journal on Advances in Signal Processing volume 2013, Article number: 59 (2013)
Abstract
Cognitive radio network with multipleinput multipleoutput is an effective method to improve not only spectrum efficiency, but also energy efficiency. In this article, a linear precoding matrix optimization algorithm, named gradientaided mutual information optimization (GAMIO), is designed to maximize the secondary users’ spectrum efficiency. Unlike the previous algorithms which were developed under a specific input assumption, the GAMIO algorithm can work without imposing any input assumption. Furthermore, a framework is also proposed to develop the energyefficient algorithm which can work with arbitrary spectrumefficient algorithm. In this way, an energyefficient algorithm, which can work under arbitrary input assumption, be developed based on the GAMIO algorithm (EEGAMIO). Numerical results indicate that either the GAMIO algorithm or the EEGAMIO algorithm shows the best performance at the present time.
1 Introduction
Fixed spectrum allocation is a traditional spectrum allocation methodology for wireless communication systems. This method causes less spectrum efficiency due to the fact that most of the available spectrum is forbidden for the users except for the licensed user (primary user). Besides the spectrum efficiency, we also need an energyefficient system in order to protect our environment and avoid the green house effect because, in recent years, the using and production of information and communication technology (ICT) contribute an increasing share to the global green house gas emissions. ICT, especially mobile telecommunication network, shows exponentially increasing energy consumption and will no doubt become a major part of energy consumption in the future. Fortunately, with the development of cognitive radio (CR) [1] network and multipleinput multipleoutput (MIMO), people find new opportunities to simultaneously improve spectrum efficiency and energy efficiency. Besides the traditional radio resource allocation research of CR network which focused on time or frequency domains, MIMO CR network can allocate the resources of space domain and offers the secondary user more freedom degrees to trade off between maximizing the spectrum/energy efficiency and minimizing the interference at the primary user.
Among all of the new technologies for MIMO CR network, linear precoding [2–6] is the most popular research topic. There are several contributions in this field: Haykin [7] gives an overview of the green wireless communication via cognitive dimension but it mainly focuses on the concept and pays less attention to the concrete application. Huang et al. [8] and Xing et al. [9] utilize game theory to maximize the sum of secondary users’ utility functions under the interferencepower constraint. In [10, 11], resource allocation of the secondary users is studied by applying the graphtheoretic models. Zhang and Liang [12] research the subspace algorithm performance in CR network by singularvalue decomposition (SVD) theory. However, all of the aforementioned works rely on the impractical Gaussian input assumption, which often leads to substantial performance degradation in real application. In order to overcome this shortage, a parameterized iterative algorithm [13], which can work under equipprobable discrete input assumption, was proposed. Though the parameterized iterative algorithm cannot directly be used in CR network because of the interference constraint, an algorithm, named branch and boundaided mutual information optimization (BAMIO), was proposed in [14, 15] to optimize the precoding matrix in CR network. However, these algorithms can work only under the equipprobable finite alphabet inputs assumption. Thanks to the contribution of [16–18], a fundamental relation between the mutual information and the minimum mean square error (MMSE) in Gaussian channel was unveiled. Furthermore, an algorithm named mercury water filling (M/WF) [19] was proposed to maximize the mutual information over parallel channel under arbitrary inputs assumption. But M/WF can work only when the channel matrix is diagonal and cannot be used in the CR network. This constraint limited its application. All of the abovementioned drawbacks motivate the research of this article. Our goal is to address a linear precoding algorithm of spectrum/energy efficiency which can be used in CR network without imposing any input assumption.
The main contributions of this article are summarized as follows: First, this article formulates the spectrum and energyefficient problem for CR network. Second, in the case of spectrum efficiency, a new algorithm based on the gradientaided mutual information optimization (GAMIO) is proposed. Different from the previous works, the GAMIO algorithm does not limit any input assumption. It is more realistic than the previous algorithm which developed under the Gaussian or any other specific input assumption. Third, in the case of energy efficiency, this article proposes a framework to develop the energyefficient algorithm which can work with any spectrumefficient algorithm. In this way, it is easy to develop an energyefficient algorithm (EEGAMIO) which can inherit the advantage of GAMIO algorithm and work under arbitrary input assumption.
The remainder of this article is organized as follows: In Section 2, the MIMO CR network system model is proposed. Based on this model, we formulate the spectrum and energyefficient optimization problems. Two typical channel capacity expressions are also proposed to examine the algorithm performance in simulation. Linear precoding algorithm is given in Section 3 to solve the proposed optimization problem. The performance evaluation is shown in Section 4. Finally, conclusions are drawn in Section 5.
Notation: Bold face uppercase letters denote matrices, bold face lower case letters denote vectors. The superscript (.)^{†} and (.)^{∗} stand for the conjugate transpose operator and the conjugate operator, respectively. The operator d i a g ( a ) denotes a diagonal matrix with elements given by a. Trace(a) denotes the trace operation of matrix a. The operator E ( . ) denotes the statistical expectation and E _{ a } ( . ) denotes the statistical expectation with respect to a. The operator ∥.∥ and . denote the l _{2} norm and the matrices determinant, respectively.
2 System model and problem formulation
As shown in Figure 1, a MIMO CR network with K primary receivers and a single pair of secondary user is considered. All of the primary users and secondary users share the same spectrum for transmission. We consider the scenario where each user is equipped with multiple antennas. The channel state information (CSI) from the secondary transmitter to every receiver is perfectly known at the secondary transmitter. Under such assumptions, the secondary transmitter is able to adapt linear precoding matrix to optimize spectrum/energy efficiency.
In practice, the CSI between primary transmitter and secondary receiver can be obtained at the secondary transmitter by sensing the primary transmitteremitted signal. On the other hand, the CSI between the secondary transmitter and secondary receiver can be obtained by periodical training when the timedivisionduplexing model is employed. Of course, it is difficult to get the perfect CSI for the secondary transmitter. Hence, the achieved spectrum/energy efficiency is an upper bounder.
Based on the model which we mentioned above, we further assume that a Gaussian channel composed N _{ t } transmitted antennas and N _{ r } received antennas is considered. Then, the MIMO signal can be formulated as
In (1), y’ and x denote the received and transmitted signal vectors, respectively, and ${\mathbf{H}}^{\prime}\in {C}^{{N}_{r}\times {N}_{t}}$ denotes the secondary users channel (in this article, we always assume that the channel matrix is full rank). z’ is the Gaussian additive noise vector. It is assumed that z ^{′}∼C N(0,σ ^{2}). B is the linear precoding matrix which we want to optimize. We denote the transmit covariance matrix of secondary users signal as ${\mathbf{R}}_{\mathbf{x}{\mathbf{x}}^{\u2020}}$, ${\mathbf{R}}_{\mathbf{x}{\mathbf{x}}^{\u2020}}=\mathbf{E}\left[\mathbf{x}{\mathbf{x}}^{\u2020}\right]$, where the expectation E[.] can be referenced from the codebook. In order to simplify the representation, we normalize y’ and H’ with the noise vector z’. So the equivalent signal can be represented as
The variables y, H, and z are $\mathbf{y}={\mathbf{R}}_{{\mathbf{z}}^{\prime}{\mathbf{z}}^{\prime}}^{0.5}{\mathbf{y}}^{\prime},\mathbf{H}={\mathbf{R}}_{{\mathbf{z}}^{\prime}{\mathbf{z}}^{\prime}}^{0.5}{\mathbf{H}}^{\prime}$, and $\mathbf{z}={\mathbf{R}}_{{z}^{\prime}{z}^{\prime}}^{0.5}{\mathbf{z}}^{\prime}$, respectively. ${\mathbf{R}}_{{\mathbf{z}}^{\prime}{\mathbf{z}}^{\prime}}$ is the covariance matrix of z’. It is clear that z∼C N(0,1). And the covariance matrix of z is an identity matrix I. The total transmitted power for secondary transmitter is denoted as P _{ t }. So the power constraint must be held as
Just like we mentioned above, we assume that there are K primary receivers in the CR network and each primary receiver is equipped with N _{ r } antennas. For every primary receiver, there must be an interference power constraint over all received antennas. The interference power constraint can be represented as
G _{ k } represents the channel matrix between secondary transmitter and K th primary receiver. From [7], it is clear that the obtained upper bound for the capacity loss of primary user is a function of the interference power Γ _{ k }. So the capacity loss of any primary user due to the secondary transmission can arbitrarily be small by controlling the interference power Γ _{ k }.
Now, it is time to formulate the efficient problem. The first goal is to design the linear precoding matrix of secondary user to maximize its spectrum efficiency under its own transmit power constraint together with a set of interference power constraints at primary users. The problem can be represented as (Pr1)
We notice that many researchers usually develop their algorithm based on the Gaussian input assumption even though the input signal is not always a Gaussian signal, because it has a very simple and elegant expression (6) and can be solved effectively by convex optimization algorithm [20, 21].
Actually, the channel capacity has different expressions in different situations. For example, when we consider the equipprobable discrete signaling constellations, such as Mary PSK, PAM, and QAM, the channel capacity is more likely to be represented as [12, 20, 21]
where M is the number of points in the signal constellation. Of course, there are many other channel capacity expressions in practical scenarios (e.g., TCM). Due to the limited space, we do not further enumerate them. It will obviously result in performance loss if we ignore the real input signal distribution. Thus, a spectrumefficient optimization algorithm which can adapt to various input signals is needed. In Section 4, we will examine GAMIO algorithm under the Gaussian input assumption and 2symbol equipprobable discrete input assumption. The simulation results show a consistent conclusion.
The second goal is to design the linear precoding matrix of secondary user to maximize its energy efficiency under not only its own transmit power constraint, but also a set of interference power constraints at primary users. The problem can be represented as (Pr2)
where C _{ t } represents the minimal channel capacity for the quality of service. The energyefficient problem has a close relationship with the spectrumefficient problem. The algorithm for Pr1 can be modified to deal with Pr2 directly.
3 Efficient optimization algorithm under arbitrary input assumption
In this section, we present the linear precoding optimization algorithm. First, we give the gradient of mutual information with respect to precoding matrix B.
3.1 Gradient of mutual information with respect to the precoding matrix B
To start with, we give the gradient of mutual information with respect to the linear precoding matrix B, which was initially proposed in [17]. Since the channel is AWGN, the conditional output of y can be represented as
where n represents the dimension of x. The mutual information is
Then the derivation of (10) is
From (9), we know that the derivation of p _{ y  x }(y) is
Therefore, (11) can be rewritten as
Notice that
Then (13) can be represented as
We also notice that p _{ x  y }<<p _{ x } for almost every y. Based on this condition and $\frac{\partial {p}_{\mathbf{x}y}}{\partial {p}_{\mathbf{x}}}(x,y)=\frac{{p}_{\mathbf{y}x}(yx)}{{p}_{\mathbf{y}}(y)}$
(RadonNikodym derivative), we can rewrite (16) as
Now, using
we can get
From the definition of conditional expectation $\mathbf{E}\left[xy\right]=\int \mathit{\text{xp}}(xy)\mathit{\text{dx}}$, we have
Because y = H B x + z, (20) can be represented as
Similarly, we also have
From (19), (21), and (22), we can get
We also notice that
We denote (24) as T which represents the MMSE matrix. From (23), (24), and the chain rule of derivation, we know that
Of course, MMSE matrix have different expressions under different input assumptions. Under the Gaussian input assumption, we have [14]
If we assume that the signal follows the finite equipprobable discrete distribution, the MMSE matrix is represented as
3.2 Optimization algorithm of Pr1
In this section, we propose an iterative algorithm to solve Pr1. The linear precoding matrix B is updated by the cuttingplane method. The idea is proposed in [22, 23]. First, the algorithm localizes a set of candidate B in a closed set. Second, half of the closed set will be eliminated from the candidate set by choosing appropriate direction. The iteration will not stop until the candidate set is small enough. Usually, a common choice of the update direction is the gradient of objective function which we want to optimized. So, from (25), We know it is
Unfortunately, we cannot promise that the whole candidate set satisfies the constraint condition. So, we have to examine them and choose the right direction in order to eliminate the bad candidates which does not satisfy the constraint. The feasible update direction is (29) which is the derivation of (3) and (4) [24]. If the constraint (3) or (4) cannot be satisfied, the GAMIO algorithm should take (29) instead of (28) for good candidates which satisfy the constraint condition.
Beside the update direction, how to choose the initial candidate set is also an important element of GAMIO algorithm. A reasonable choice is the minimal size ellipsoid contained the optimal precoding matrix B ^{∗}. So, we should take the initial set a and center B _{ 0 } as (30) based on the transmitted power constraint (3).
According to the property of gradient, the GAMIO algorithm will converge to the global optimal value if the Pr1 is a convex/concave problem. Of course, not all of input signal distribution promise Pr1 is a convex/concave problem. The GAMIO algorithm performance will greatly depend on the initial value if the Pr1 is not a convex/concave problem. Usually, it will converge to a local optimal value. By the way, it should be noted that the Gaussian input assumption leads Pr1 to a concave problem. So, the GAMIO algorithm can achieve the global optimal value. It has the same performance as the interiorpoint algorithm which is considered as the best algorithm under the Gaussian input assumption.
From [25], we know that the GAMIO, being a kind of cuttingplane algorithm, is a polynomialtime algorithm. The iteration complexity is decided by the threshold ε _{1} and initial set a. Without loss of generality, we assume that the channel capacity is Lipschitz. Then the iteration complexity is $O(2{n}^{2}ln(\frac{\mathit{\text{LR}}}{{\epsilon}_{1}}))$. L represent the Lipschitz condition constant. R is decided by the initial set a and represents the volume of ellipsoid. The most complex part of GAMIO algorithm is to estimate the MMSE matrix based on the current precoding matrix. In order to reduce the timeconsuming operations, we take the trick which employed in [15] and update the MMSE matrix recursively. Let s be the number of operations in the MMSE matrix evaluation. Then the computational complexity of GAMIO algorithm is $O(2s{n}^{2}ln(\frac{\mathit{\text{LR}}}{\epsilon}))$. Of course, s is different under different input assumption and can be adjusted according to the algorithm precision. The GAMIO algorithm is concluded as follows.
GAMIO algorithm under arbitrary input assumption

1.
Initial the candidate set:
$$\begin{array}{l}\mathbf{A}=\mathbf{diag}([{P}_{t},{P}_{t},\dots ,{P}_{t}]),\phantom{\rule{2em}{0ex}}\\ {B}_{0}=\mathbf{diag}([{P}_{t}^{\frac{1}{2}},{P}_{t}^{\frac{1}{2}},\dots ,{P}_{t}^{\frac{1}{2}}]),\phantom{\rule{2em}{0ex}}\\ \phantom{\rule{3em}{0ex}}\left\{\begin{array}{l}{d}_{1}={H}^{\u2020}\mathit{\text{HBT}}\\ {d}_{2}=\mathbf{B}{R}_{x{x}^{\u2020}}\\ {d}_{3,k}={\mathbf{G}}_{k}^{\u2020}{\mathbf{G}}_{k}\mathbf{B}{R}_{x{x}^{\u2020}},k=1,\dots ,K;\end{array}\right.\phantom{\rule{2em}{0ex}}\end{array}$$ 
2.
Repeat:

(a)
If Trace $({\mathit{\text{BR}}}_{{\mathbf{xx}}^{\u2020}}{\mathbf{B}}^{\u2020})>{P}_{t},$
$$\phantom{\rule{3em}{0ex}}d={d}_{2}.$$
Else ifTrace (
Else

(b)
Update as
$$\begin{array}{l}\phantom{\rule{3em}{0ex}}t={(\text{Trace}({}_{d}^{\u2020}\mathbf{A}d))}^{\frac{1}{2}}\\ \phantom{\rule{3em}{0ex}}{B}_{0}={B}_{0}\frac{\mathbf{A}}{N+1}\frac{d}{t}\\ \phantom{\rule{3em}{0ex}}\mathbf{A}=\frac{{N}^{2}}{{N}^{2}1}(\mathbf{A}\frac{2}{N+1}\mathbf{A}\frac{d}{t}{(\frac{d}{t})}^{T}{\mathbf{A}}^{T});\end{array}$$ 
3.
Until t<ϵ _{1}, end.
3.3 Optimization algorithm of Pr2
Just as we mentioned before, Pr2 closely coupled with Pr1. In general, Pr2 is not a concave problem. So, we cannot solve it by traditional optimization method directly. Thus, we try to find a suboptimal solution instead. We notice that the Pr2 decrease as its denominator value increases, while it will increase with the raise of its numerator value. So, we can minimize the denominator value and maximize the numerator value separately to optimize Pr2. It can be done by combining arbitrary spectrum efficient optimization algorithm with bisection method. So, the Pr2 can be separated into two parts. The first part can be represented as (Pr3)
The only difference between Pr3 and Pr1 is the constraint power ${P}_{i}^{\ast}$ which represents current power constraint during i th iteration. Obviously, ${P}_{t}\ge {P}_{i}^{\ast}$. If we have $\mathcal{I}({x}^{\ast},{y}^{\ast})\ge {C}_{t}$, then
Otherwise, if I(x ^{∗} , y ^{∗})<C _{ t }, then
where $\mathcal{I}({x}^{\ast},{y}^{\ast})$ represents the optimal value of Pr3. The initial values of P _{Min} and P _{Max} are 0 and P _{ t }, respectively. The iteration will not stop until P _{Min} and P _{Max} are close enough. The computational complexity of this new iterative approach is tightly related with the spectrumefficient algorithm for Pr3. Without any doubt, although arbitrary spectrumefficient algorithm can work in this frame to optimize the energy efficiency, the energyefficient algorithm which is developed based on the GAMIO algorithm (EEGAMIO) is a better choice than any other algorithm to deal with various input signal. From [23], we know that the iteration complexity of bisection method is
where ϵ _{2} represents the threshold of bisection method. Based on the conclusion of Section 3.2, we can get the iteration and computational complexity of EEGAMIO energyefficient algorithm are $O(2{n}^{2}\phantom{\rule{1pt}{0ex}}ln\phantom{\rule{1pt}{0ex}}(\frac{\mathit{\text{LR}}}{{\epsilon}_{1}})\phantom{\rule{1pt}{0ex}}{log}_{2}\phantom{\rule{1pt}{0ex}}\frac{{P}_{t}}{{\epsilon}_{2}})$ and $O(2s{n}^{2}\phantom{\rule{1pt}{0ex}}ln\phantom{\rule{1pt}{0ex}}(\frac{\mathit{\text{LR}}}{{\epsilon}_{1}})\phantom{\rule{1pt}{0ex}}{log}_{2}\phantom{\rule{1pt}{0ex}}\frac{{P}_{t}}{{\epsilon}_{2}})$ respectively. Table 1 summarize the complexity of the proposed algorithms and the framework of the energyefficient algorithm is concluded as follows.
The framework of energyefficient algorithm

1.
Initial the searching set:
$$\begin{array}{l}{P}_{min}=0,\begin{array}{ll}{P}_{max}={P}_{t},& {P}_{0}^{\ast}=\frac{{P}_{min}+{P}_{max}}{2};\end{array}\end{array}$$ 
2.
Repeat:

(a)
Solve the problem
$$\begin{array}{cc}max:& \mathcal{I}(\mathbf{x},y)\\ \mathrm{s.t}:& \begin{array}{c}\text{Trace}(\mathbf{B}{R}_{x{x}^{\u2020}}{B}^{\u2020})\le {P}_{i}^{\ast}\\ \begin{array}{cc}\text{Trace}({\mathbf{G}}_{k}\mathbf{B}{R}_{x{x}^{\u2020}}{B}^{\u2020}{\mathbf{G}}_{k}^{\u2020})\le {\Gamma}_{k},& k=1\mathrm{...K};\end{array}\end{array}\end{array}$$ 
(b)
If $\mathcal{I}({x}^{\ast},{y}^{\ast})\ge {C}_{t}$,
$$\begin{array}{ll}{P}_{max}={P}_{i}^{\ast},& {P}_{i+1}^{\ast}=\frac{{P}_{min}+{P}_{max}}{2};\end{array}$$
Else if $\mathcal{I}({x}^{\ast},{y}^{\ast})<{C}_{t}$,

3.
Until P _{max}−P _{min}≤ϵ _{2}, end.
4 Simulation results
In this section, we demonstrate some simulation results. The simulation is separated into two parts. For simplicity, we assume that there is only one primary user in the CR network.
4.1 Spectrum efficiency
In first part, we test the spectrum efficiency of GAMIO algorithm. We take (6) and (7) as the channel capacity for the Gaussian input assumption and the finite equipprobable discrete input assumption, respectively. The channel matrix from secondary transmitter to primary receiver is $G=\left[\begin{array}{ll}2& 1\\ 3& 4\end{array}\right]$and the channel matrix between secondary user is $H=\left[\begin{array}{ll}2& 1\\ 1& 1\end{array}\right]$. The average signalnoiseradio (SNR) of CR MIMO system is given by $\text{SNR}=\frac{{P}_{t}}{{\sigma}^{2}}$. The total transmit power P _{ t } is 2 W. The interferencepower constraint Γ is 1 W.
In Figures 2 and 3, we test the achievable channel capacity of secondary user under different situations. In Figure 2, we assume that the input signal follows the Gaussian distribution. So, the channel capacity can be expressed as (6). The line denoted with circle signs represents the performance of interiorpoint algorithm. It has the best performance. The line with ’+’ represents the performance of GAMIO algorithm. It has the same performance as interiorpoint algorithm because the channel capacity (6) is a concave problem under the Gaussian input assumption. The line with ’ ∗’ represents the performance of subspace algorithm which is a wellknown suboptimal algorithm [12]. The line with ’ ×’ represents the performance of BAMIO algorithm and much lower than the others because the BAMIO algorithm is developed under equipprobable discrete input assumption.
In Figure 3, we assume that the input signal follows the 2symbol equipprobable discrete distribution (e.g., BPSK). So the channel capacity can be expression as (7) with M=2. Just as what Figure 2 shows, the first line with ’+’ represents the performance of GAMIO algorithm and the second line with ’ ×’ represents the performance of BAMIO algorithm. They almost have the same performance. Similarly, the third line denoted with circle signs represents the performance of interiorpoint algorithm and the fourth line with ’ ∗’ represents the performance of subspace algorithm. Figure 3 shows that the performance of GAMIO and BAMIO algorithm is much better than interiorpoint algorithm and subspace algorithm. It almost reaches the best capacity, 2 bit/s, when the SNR is 20 dB. On the other hand, the capacity of interiorpoint and subspace algorithm lower than 1 bit/s because they intend to allocate too much energy to the best channel and degrade the channel capacity.
Figures 2 and 3 show that, unlike the traditional spectrumefficient algorithm which developed under specific input assumption, the GAMIO algorithm considers the specific input signal distribution and achieves the best performance under various input assumption. This result is consistent with our analysis and previous research results [8, 10].
Figure 4 gives the convergence performance results of GAMIO algorithm with different SNR. As shown in Figure 4, the algorithm converges to the optimal value before the iteration ends. On the other hand, we can also find that the convergence rate improves with SNR increasing.
In order to unveil the meaning of different channel capacity expression, we present a BPSK MIMO system as Figure 5 shows. Figure 6 gives the BER of this system with different SNR, Figure 6 indicates that the BER of GAMIO and BAMIO algorithms is obviously lower than the interiorpoint algorithm and SVD algorithm. This result is consistent with Figure 3. And Figure 2 cannot reflect this property. Figures 2, 3, and 6 show that the input signal distribution affect the real performance greatly in practical scenarios. So, the GAMIO algorithm is more useful than any other spectrumefficient algorithm, which can work only under specific input assumption, for various practical system design.
4.2 Energy efficiency
In the second part, we simulate the energy efficiency problem which is defined as Pr2. Just like the first part, we also take into account two kinds of channel capacity expression. The first one is based on the Gaussian input assumption and the other one based on 2symbol finite equipprobable discrete input assumption respectively. All of the assumptions are the same as the first part except adding a minimum channel capacity constraint, C _{ t }= 1 bit/s.
In Figure 7, we test the achievable energy efficiency under different situations. Under the Gaussian input assumption, we develop energyefficient algorithm based on the interiorpoint algorithm (EEinteriorpoint), subspace algorithm (EESVD), and GAMIO algorithm (EEGAMIO), respectively. The results are noted with plus, star, and circle, respectively. As Figure 7 shown, the performances of EEinteriorpoint algorithm and EEGAMIO algorithm have no significant difference under Gaussian input assumption. EESVD algorithm cannot satisfy the capacity constraint when the SNR is less than 12 dB. But it is almost the same as the other algorithms when SNR is beyond 12 dB. Similarly, we also develop energyefficient algorithm based on the BAMIO algorithm and GAMIO algorithm under the 2symbol equipprobable discrete input assumption. The performance line is annotated by diamond and invertedtriangle, respectively. There are no performance results for EEinteriorpoint algorithm and EESVD algorithm because they cannot satisfy the capacity constraint under the 2symbol equipprobable discrete input assumption. Figure 7 shows that, under either Gaussian or 2symbol equipprobable discrete input assumption, the EEGAMIO algorithm can reach the best performance.
The last line denoted with square sign represents the energy efficiency of traditional interiorpoint algorithm under the Gaussian input assumption. It is much lower than all of the others. This result shows that the proposed framework of energyefficiency is effective.
Figure 8 shows the energy efficiency of MIMO system which Figure 5 figured. The lines which denoted with plus, star, circle, diamond, and invertedtriangle sign are same as Figure 7. As Figure 8 shown, although the Gaussian input assumption can keep high energy efficiency, it brings more loss in the BER performance. So the Gaussian input assumption is more difficult to satisfy the channel capacity constraint than the 2symbol equipprobable discrete input assumption in the MIMO system as Figure 5. It also reveal that the input signal distribution affect the energy efficiency greatly.
5 Conclusion
In this article, linear precoding matrix optimization for the secondary user in a CR network is considered. We set both spectrum and energy efficiency as our goals. Unlike previous work, we do not limit the input signal following specific distribution. Instead, we propose a spectrumefficient algorithm (GAMIO), which can work well under not only the Gaussian input assumption but also any other kind of input assumption. With this advantage, the GAMIO algorithm is more practical than existing algorithms, such as interiorpoint algorithm, SVD algorithm, or BAMIO algorithm. In the case of energy efficiency, we propose a framework to optimize the energy efficiency. With this framework, we can get different energyefficient algorithms based on different spectrumefficient algorithms. In order to deal with various signal, the EEGAMIO algorithm is a smart choice. The simulation results show that these algorithms can significantly increase the energy efficiency.
References
 1.
Mitola J, Maguire GQ: Cognitive radios: making software radios more personal. Personal Commun 1999, 6(4):1318. 10.1109/98.788210
 2.
Scaglione A, Staica P, Barbarossa S: Optimal designs for spacetime linear precoders and decoders. IEEE Trans. Signal Process 2002, 50(5):10511064. 10.1109/78.995062
 3.
Sampath H, Paulraj A: Linear precoding for spacetime coded systems with known fading correlations. IEEE Commun. Lett 2004, 6(6):234241.
 4.
Zhou S, Giannakis GB: Optimal transmitter eigen beamforming and spacetime block coding based on channel correlations. IEEE Trans. Commun 2003, 50(5):16731690.
 5.
Kermoal JP: A stachastic MIMO radio channel model with expermental validation. IEEE J. Sel. Areas Commun 2002, 20(6):1211—1226.
 6.
Costa MHM: Writing on dirty paper. IEEE Trans. Inf. Theory 1983, 29(3):439441. 10.1109/TIT.1983.1056659
 7.
Haykin S: Cognitive radio: brainempowered wireless communications. J. Sel. Areas Commun 2005, 23(2):201220.
 8.
Huang J, Berry R, Honig ML: Auctionbased spectrum sharing. ACM/Springer Mob. Netw. Appl. J. (MONET) 2006, 11(3):405418. 10.1007/s110360065192y
 9.
Xing Y, Mathur CN, Haleem MA, Chandramouli R, Subbalakshmi KP: Dynamic spectrum access with QoS and interference temperature constraints. Trans. Mob. Comput 2007, 6(4):423433.
 10.
Zheng H, Peng C: Collaboration and fairness in opportunistic spectrum access. In Proceedings of the IEEE International Conference on Communication. USA: IEEE Press; 2005:31323136.
 11.
Wang W, Liu X: Listcoloring based channel allocation for openspectrum wireless networks. In Proceedings of the IEEE Vehicular Technology Conference (VTC). USA: IEEE Press; 2005:690694.
 12.
Zhang R, Liang YC: Exploiting multiantennas for opportunistic spectrum sharing in cognitive radio networks. IEEE J. Sel. Topic Signal Process 2008, 2: 88102.
 13.
Xiao C, Zheng YR, Ding Z: Globally optimal linear precoders for finite alphabet signals over complex vector Gaussian channels. IEEE Trans. Signal Process 2011, 59(7):33013314.
 14.
Zeng W, Xiao C, Lu J: Globally optimal precoder design with finitealphabet input for cognitive radio networks. IEEE J. Sel. Topics Commun 2012, 30(10):18611874.
 15.
Zeng W, Xiao C, Wang M, Lu J: Linear precoding for finitealphabet inputs over MIMO fading channels with statistical CSI. IEEE Trans. Signal Process 2012, 60(6):31343148.
 16.
Guo D, Shitz SS, Verdu S: information, Mutual minimum meansquare error in Gaussian channels. IEEE Trans. Inf. Theory 2005, 51(4):12611282. 10.1109/TIT.2005.844072
 17.
Palomar DP, Verdu S: Gradient of mutual information in linear vector Gaussian channels. IEEE Trans. Inf. Theory 2006, 52: 141154.
 18.
Kay SM: Fundamentals of Statistical Signal Processing: Estimation Theory (PrenticeHall). NJ: Englewood Cliffs; 1993.
 19.
Lozano A, Tulino AM, Verdu S: Optimum power allocation for parallel Gaussian channels with arbitrary input distributions. IEEE Trans. Inf. Theory 2006, 52(7):30333051.
 20.
Duncan TE: On the calculation of mutual information. SIAM. J. Appl. Math 1970, 19: 215220.
 21.
Jengren G, Skoglund M, Ottersen B: Combining beamforming and orthogonal spacetime block coding. IEEE Trans. Inf. Theory 2002, 48(3):611617. 10.1109/18.985950
 22.
Cover TM, Thomas JA: Elements of Information Theory. New York: Wiley; 1991.
 23.
Boyd S, Vandenberghe L: Convex Optimization. UK: Cambridge University Press Cambridge; 2004.
 24.
Magnus JR, Neudecker H: Matrix Differential Calculus With Applications in Statistics and Econometrics. New York: Wiley; 1999.
 25.
Boyd S: EE392o Course Notes Stanford University. 2004.http://www.stanford.edu/class/ee392o/elp.pdf
Acknowledgements
The authors would like to thank the associate editor and reviewers for helpful comments and suggestions that greatly improved the quality of the paper. This study was supported by the Commission of National “863” project (NO. 2012AA011402), National “973” project (NO. 2012CB316000), National Foundation Science of China (60832008), National Major Project (NO. 2012ZX03004004002), Independent Research Project (NO. 2010TH20302), and Research project (Research of Dynamic Spectrum Allocation Technology in Cognitive Radio).
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Zhu, R., Zhao, Y., Li, Y. et al. Optimal linear precoding for opportunistic spectrum sharing under arbitrary input distributions assumption. EURASIP J. Adv. Signal Process. 2013, 59 (2013). https://doi.org/10.1186/16876180201359
Received:
Accepted:
Published:
Keywords
 Energy Efficiency
 Primary User
 Channel State Information
 Channel Capacity
 Secondary User