Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: in multi - user mimo systems , one main challenge is the increased complexity and energy consumption of the signal processing to mitigate the interferences between multiple co - channel users . to achieve energy efficient transmission , the multi - user mimo system with very large antenna arrays at each base station ( known as `` massive mimo '' system ) has been advocated recently @xcite . the key result is that , with very large antenna arrays at each base station , both the intracell and intercell interferences can be substantially reduced with simple linear beamforming ( bf ) processing @xcite .    on the other hand , full - duplex relaying ( fdr ) is a promising approach to improve the spectral efficiency ( se ) of relaying network while retains the merits of half - duplex relaying ( hdr ) ( e.g. , path loss reduction ) . in fdr , the relay transmits and receives simultaneously at the same frequency and time , but at the cost of a strong echo interference ( ei ) due to signal leakage between the relay output and input @xcite . to mitigate ei , three approaches have been investigated , i.e. , 1 ) passive cancellation , 2 ) time - domain cancellation @xcite and 3 ) spatial suppression @xcite . the passive cancellation relies on a combination of path loss , cross - polarization and antenna directionality @xcite . the time - domain cancellation is based on the fact that ei signal is known at full - duplex node . thus , it can be regenerated and removed in time - domain @xcite . in spatial suppression , ei is mitigated with the multiple transmit / receive antennas by approaches such as null - space projection @xcite . inspired by these works , a number of works have dedicated to the study of fdr protocol on both theory and testbed ( see @xcite and the references therein ) . to achieve spectral and energy efficient transmissions of multiple source - destination pairs , recent works sought to incorporate both hdr @xcite and fdr @xcite-@xcite with massive mimo .    however , the aforementioned works on massive mimo are actually based on the assumption that the base stations or relays are equipped with a large number of high - quality transmit / receive radio frequency ( rf ) chains ( which are expensive and power - hungry ) . in contrast to conventional mimo system ( e.g. , at most 8 antennas in lte system ) , massive mimo must be built with low - cost components @xcite since the deploy cost and energy consumption of circuits will increase dramatically as the number of antennas grows very large . such low - cost components are prone to hardware imperfections ( e.g. , phase noise , nonlinear power amplifier , i / q imbalance , nonlinear low - noise amplifier and adc impairments ) , which must be considered in the design of practical massive mimo system . this paper focuses on the transceiver design for massive mimo full - duplex relaying ( mm - fdr ) with hardware impairments . the effect of hardware impairments is modeled using transmit / receive distortion noises @xcite . there are several challenges in the design of practical transceiver scheme in the considered system . the first is : _ how to deal with the ei cancellation without instantaneous ei channel ? _ ei cancellation is a critical problem in mm - fdr transceiver design which is not only important to reduce ei power , but also useful to reduce distortion noises caused by hardware impairments at the relay ( as will be shown in section iv ) . different from fdr with small - scale relay antenna arrays @xcite-@xcite , the instantaneous ei channel is usually not easy to obtain in mm - fdr . this is because the learning of ei channel requires training sequence with length not less than the number of relay antennas @xcite , which is prohibitive in mm - fdr as the channel coherent time is limited . the lack of ei channel makes traditional ei cancellation techniques ( e.g. , time - domain cancellation and spatial suppression ) difficult to apply . although passive cancellation does not relies on ei channel estimation , it usually can not provide satisfactory performance when used alone @xcite . another problem is : _ how to suppress distortion noises caused by hardware impairments at sources and destinations . _ different from multi - user interference ( mui ) , the distortion noise caused by transmit imperfection of source can be viewed as an interference signal with the same channel as transmit data . thus , it can not be suppressed by relay antenna arrays during coherent combining ( similar problem appears in reception at destinations ) . this causes performance ceiling on achievable rate as the number of relay antennas grows large , which degrades the gain of massive mimo significantly .    in this paper , we propose practical transceiver scheme for mm - fdr with hardware impairments considering the above problems . different from the traditional mm - fdr protocol @xcite ( where sources and destinations are equipped with single antenna ) , we consider a general model where sources and destinations are allowed to equip with multiple antennas . the contributions are summarized as follows :    * we first examine the limitation of traditional mm - fdr protocol under hardware impairments . in particular , we derive the upper bound on end - to - end achievable rate for traditional mm - fdr protocol with linear processing at the relay . the bound reveals that the achievable rate is limited by the hardware impairments at the sources and destinations , and performance ceiling appears as the number of relay antenna grows large . the result also implies that it is impossible to cancel the `` ceiling effect '' with linear processing , if sources and destinations are only equipped with single antenna . * based on the upper bound analysis , we propose a hardware impairments aware transceiver scheme ( hia scheme ) to mitigate the distortion noises by exploiting the statistical channel knowledge and antennas arrays of sources / destinations . the scheme needs no instantaneous knowledge of ei channel . the asymptotic end - to - end achievable rate of mm - fdr with hia scheme ( hia - mm - fdr ) is derived and the scaling behaviors of mui , ei and distortion noises are determined . * a joint degree of freedom and power optimization ( jdpo ) algorithm is presented to further improve the se of hia - mm - fdr . the paper is organized as follows . in section ii , we review the related work . the system model is described in section iii . the upper bound of achievable rate is analyzed in section iv . the hia scheme is proposed in section v and the jdpo algorithm is presented in section vi . simulation results are presented in section vii and some conclusions will be drawn at last . : @xmath0 and @xmath1 denote the expectation and variance . @xmath2 is @xmath3 identity matrix . @xmath4 , @xmath5 and @xmath6 denote conjugate , transpose and conjugate - transpose , respectively . @xmath7 , @xmath8 , @xmath9_{ij}$ ] , @xmath10 and @xmath11 denote the spectral radius , trace , @xmath12th element of matrix , @xmath13th largest eigenvalue and eigenvector with respect to @xmath13th largest eigenvalue of matrix @xmath14 , respectively . @xmath15 scales with @xmath16 means @xmath17 . the design and performance of using unlimited number of antennas at the base station in cellular system were first considered for independent antennas @xcite , and soon extended to the scenarios with spatial correlated antennas @xcite-@xcite . the asymptotic sinr for single cell cellular system was analyzed in @xcite . it has been shown that , with very large antenna arrays at the base station , a deterministic sinr ( also called the `` deterministic equivalent '' of sinr ) which depends only on the large - scale fading of channels can be achieved , if the transmit power is scaled by @xmath18 with perfect csi and @xmath19 with imperfect csi ( @xmath18 denotes the number of base station antennas ) . the authors in @xcite-@xcite have done considerable work to derive the `` deterministic equivalent '' of sinr for massive mimo system with spatial correlated antennas .    in the field of massive mimo relaying , the se and ee of hdr with very large relay antenna arrays were investigated in @xcite . the mm - fdr with decode - and - forward ( df ) relay was first introduced in @xcite , and analyzed in end - to - end achievable rate as linear processing is employed . the asymptotic performance of amplify - and - forward based mm - fdr was considered in @xcite and the scaling behavior of the infinitely repeating echo interference was determined .    only a few works considered the effect of hardware impairments on massive mimo system . a constant envelope signal design has been proposed in @xcite to facilitate the use of power - efficient rf power amplifiers . the authors in @xcite presented a low peak - to - average - power ratio ( par ) precoding solution to enable efficient implementation using non - linear rf components in massive mimo system . in fdr with small - scale antennas , the optimal precoding under limited adc dynamic range was studied in @xcite . however , these works can not provide much insight for the effect of hardware impairments on mm - fdr . consider the network with @xmath20 source - destination ( s - d ) pairs and a single full - duplex relay @xmath21 , where source @xmath22 wishes to communicate with destination @xmath23 ( @xmath24 ) with the help of @xmath21 . the relay adopts the df policy . it is assumed that the sources and destinations are equipped with @xmath25 and @xmath26 antennas respectively , while the relay is equipped with @xmath27 antennas ( @xmath28 for reception and @xmath29 for transmission ) . we are interested in the large-@xmath30 regime , i.e. , @xmath31 . @xmath25 , @xmath26 and @xmath20 can be either fixed or scale with @xmath32 . let @xmath33 be the channel matrix from @xmath22 to receive antenna array of relay and let @xmath34 be the channel matrix from @xmath23 to transmit array of relay . let @xmath35 denote the ei channel matrix between transmit and receive arrays of relay . the spatial correlation of each mimo channel is characterized by the kronecker model @xcite . thus , @xmath36 can be expressed as @xmath37 , where @xmath38 represents the large - scale fading . @xmath39 and @xmath40 characterize the spatial correlation of received signals across receive array of relay and that of transmitted signals across transmit array of @xmath22 . @xmath41 consists of the random components of channels whose elements are i.i.d with distribution @xmath42 . based on the same model , the channel matrices @xmath43 and @xmath44 can be expressed as @xmath45 and @xmath46 . to facilitate the analysis , we assume the following conditions on correlation matrices , i.e. , @xmath47 , where @xmath48 ,    * * a1 * : the spectral radius of @xmath49 is bounded by a constant , i.e. , @xmath50 . * * a2 * : @xmath49 is a hermitian and teoplitz matrix and has unit diagonal elements .    the former is a common assumption in the studies of massive mimo which follows from energy conservation @xcite and * a2 * corresponds to the case of uniform linear array ( ula ) @xcite .    to characterize the effect of hardware imperfections , we adopt the new signal model from @xcite . we model the effect of imperfect transmit rf chain by adding , per transmit antenna , an independent zero - mean gaussian `` transmit distortion noise '' , whose power is proportional to the signal power transmitted at that antenna . the experimental results in @xcite have shown that the independent gaussian distortion noise model closely captures the joint effect of imperfect components in transmit rf chain . let @xmath51 \in \mathbb c^{n_s\times 1}$ ] and @xmath52\in \mathbb c^{n_t\times 1}$ ] denote the transmit vectors of source @xmath22 and relay at time instant @xmath53 . based on the above model , the distorted transmit signals can be expressed as @xmath54 = { { \mathbf{x}}_{s , k}}\left [ u \right ] + { { \mathbf{t}}_{s , k}}\left [ u \right ] \hspace{.5em } { \rm and } \hspace{.5em } { { \mathbf{\tilde x}}_r}\left [ u \right ] = { { \mathbf{x}}_r}\left [ u \right ] + { { \mathbf{t}}_r}\left [ u \right]\ ] ] where the distortion noises @xmath55\sim { \cal cn}\left ( { 0,{\nu _ { s , k}}{\text{diag}}\left ( { \mathbb{e}\left [ { { { \mathbf{x}}_{s , k}}\left [ u \right]{\mathbf{x}}_{s , k}^h\left [ u \right ] } \right ] } \right ) } \right)$ ] and @xmath56\sim{\cal c}{\cal n}\left ( 0,\nu _ r \times \right.$ ] @xmath57{{\bf{x}}_r^h}\left [ u \right ] } \right ] } \right ) \right)$ ] . note that @xmath58 ( @xmath59 ) characterizes the level of transmit imperfection . for example , @xmath60 ( @xmath61 ) corresponds to the conventional assumption of perfect transmit rf chains . the quality of transmit rf chains degrades as @xmath62 ( @xmath63 ) increases . we model the effect of imperfect receive rf chain by adding , per receive antenna , an independent zero - mean gaussian `` receive distortion noise '' , whose variance is proportional to the signal power received at that antenna . more precisely , let @xmath64\in \mathbb c^{n_r\times 1}$ ] and @xmath65 \in \mathbb c^{n_d\times 1}$ ] be the undistorted received signals of the relay and destination @xmath23 at time instant @xmath53 , the distorted received signals can be expressed as @xmath66 = { { \bf{y}}_r}\left [ u \right ] + { { \bf{r}}_r}\left [ u \right ] \hspace{.5em }   { \rm and } \hspace{.5em } { { { \bf{\tilde y}}}_{d , k}}\left [ u \right ] = { { \bf{y}}_{d , k}}\left [ u \right ] + { { \bf{r}}_{d , k}}\left [ u \right ]      \end{aligned}\ ] ] where the distortion noises @xmath67 \sim { \cal c}{\cal n } ( { 0,{\mu _ r}{\rm{diag } } ( { { \mathbb{e } } [ { { { \bf{y}}_r } [ u ] { { \bf{y}}_r^h } [ u ] } ] } ) } ) $ ] and @xmath68 \sim{\cal cn } ( 0,{\mu _ { d , k } } { \text{diag } } ( { \mathbb{e } [ { { { \mathbf{y}}_{d , k } } [ u ] { \mathbf{y}}_{d , k}^h [ u ] } ] } ) ) $ ] . @xmath69 ( @xmath70 ) characterizes the level of receive imperfection . the experimental studies in @xcite have shown that the independent gaussian noise model is a good approximation to the joint effect of imperfect components in receive rf chain .    at time instant @xmath53 , all sources transmit signals @xmath71 $ ] ( @xmath72 ) to the relay simultaneously . meanwhile , the relay broadcasts the decoded signals to the destinations . based on the models in ( [ 2 - 3 ] ) and ( [ 2 - 4 ] ) , the received signals at the relay and @xmath23 can be expressed as @xmath73 = \sum\limits_{k = 1}^k { { \mathbf{h}}_{sr , k}}\left ( { { { \mathbf{x}}_{s , k}}\left [ u \right ] + { { \mathbf{t}}_{s , k}}\left [ u \right ] } \right ) \\          & + { { \mathbf{h}}_{ei}}\left ( { { { \mathbf{x}}_r}\left [ u \right ] + { { \mathbf{t}}_r}\left [ u \right ] } \right ) + { { \mathbf{r}}_r}\left [ u \right ] + { { \mathbf{n}}_r}\left [ u \right ] \\          & { { \mathbf{\tilde y}}_{d , k}}\left [ u \right ] = { \mathbf{h}}_{rd , k}^h\left ( { { { \mathbf{x}}_r}\left [ u \right ] + { { \mathbf{t}}_r}\left [ u \right ] } \right ) + { { \mathbf{r}}_{d , k}}\left [ u \right ] + { { \mathbf{n}}_{d , k}}\left [ u \right ]      \end{aligned}\ ] ] to keep the complexity low , it is assumed that a singe data steam is transmitted at each source and each node employs only linear processing . in particular , @xmath22 transmits the unit - power symbol @xmath74 $ ] using the unitary bf vector @xmath75 . therefore , the transmit vector of @xmath22 can be expressed as @xmath71 = \sqrt{e_{s , k}}{\bf p}_{s , k } s_k[u]$ ] , where @xmath76 denotes the transmit power . relay combines the received signal by multiplying the receive bf matrix @xmath77 $ ] , i.e. , @xmath78 = { \bf{w}}_r^h{{{\bf{\tilde y}}}_r}\left [ u \right]$ ] . the @xmath79th element of @xmath78 $ ] @xmath80   = \sqrt{e_{s , k}}{\mathbf{w}}_{r , k}^h{{\mathbf{h}}_{sr , k}}{{\mathbf{p}}_{s , k}}{s_k}\left [ u \right ] \\          & +   \sqrt{e_{s , j } } \sum\limits_{j = 1,j \ne k}^k { { \mathbf{w}}_{r , k}^h{{\mathbf{h}}_{sr , j}}{{\mathbf{p}}_{s , j}}{s_j}\left [ u \right ] }   + { \mathbf{w}}_{r , k}^h\sum\limits_{j = 1}^k { { { \mathbf{h}}_{sr , j}}{{\mathbf{t}}_{s , j}}\left [ u \right ] } \\          & + { \bf{w}}_{r , k}^h{{\bf{h}}_{ei}}\left ( { { { \bf{x}}_r}\left [ u \right ] + { { \bf{t}}_r}\left [ u \right ] } \right ) + { \bf{w}}_{r , k}^h{{\bf{r}}_r}\left [ u \right ] + { \bf{w}}_{r , k}^h{{\bf{n}}_r}\left [ u \right ]      \end{aligned}\ ] ] is used to decode the symbol of @xmath22 . the relay forwards the decoded symbol using transmit bf matrix @xmath81 $ ] . the transmit vector of the relay can be expressed as @xmath82\ ! = \ ! { { \mathbf{w}}_t}{\mathbf{\lambda } _ r^{1/2}}{\mathbf{s}}\left [ { u - d } \right]$ ] , where @xmath83\!=\![{{s_1}}\left [ { u - d } \right ] , \cdots,{{s_k}}\left [ { u - d } \right]]^t$ ] and @xmath84 is the power allocation matrix at the relay . @xmath85 denotes the processing delay of the relay . to meet the relay s power constraint , @xmath82 $ ] must satisfy @xmath86{\bf{x}}_r^h\left [ u \right ] } \right ] } \right ) \!\ ! = \!\ ! \sum_{l = 1}^k { { e_{r , l } } } \!\!\le\ ! \!e_r^{\rm max}$ ] . @xmath23 uses the unitary bf vector @xmath87 to combine the received signal @xmath88 $ ] . the combined signal is expressed as @xmath89 =   \sqrt{e_{r , k}}{\mathbf{p}}_{d , k}^h{\mathbf{h}}_{rd , k}^h{{\mathbf{w}}_{t , k}}{s_k}\left [ { u - d } \right ] \\           & + { \mathbf{p}}_{d , k}^h\sum\limits_{j = 1,j \ne k}^k \sqrt{e_{r , j } } { { \mathbf{h}}_{rd , k}^h{{\mathbf{w}}_{t , j}}{s_j}\left [ { u - d } \right ] }   \\          & + { \mathbf{p}}_{d , k}^h{\mathbf{h}}_{rd , k}^h{{\mathbf{t}}_r}\left [ u \right ] + { \mathbf{p}}_{d , k}^h{{\mathbf{r}}_{d , k}}\left [ u \right ] + { \mathbf{p}}_{d , k}^h{{\mathbf{n}}_{d , k}}\left [ u \right ]      \end{aligned}\ ] ] this section analyzes the upper bound on achievable rate of mm - fdr with single antenna sources and destinations . the goal is to reveal the fundamental limitation of previous mm - fdr protocol with single antenna sources and destinations in combating distortion noises . moreover , the bound provides us important insight on the design of practical transceiver scheme to mitigate the distortion noises , when multiple antennas are available at sources and destinations ( see section v ) .      as @xmath90 , we set @xmath91 , replace @xmath92 ( @xmath93 ) with @xmath94 , and let @xmath95 $ ] and @xmath96 = { \left [ { { \mathbf{t}}_{s,1}^t\left [ u \right ] , \cdots , { \mathbf{t}}_{s , k}^t\left [ u \right ] } \right]^t}$ ] . we consider a block - fading channel with coherence time @xmath97 ( symbol times ) , where @xmath98 are used for uplink training of each source / destination , and the remaining @xmath99 are used for data transmission . the upper bound on achievable rate is obtained by assuming perfect knowledge of @xmath100 , @xmath101 and @xmath102 can be provided with pilot signals , and meanwhile , the muis ( i.e. , the second terms of right - hand side of ( [ 2 - 6 ] ) and ( [ 2 - 7 ] ) ) and ei term @xmath103 $ ] can be somehow cancelled . this gives us the following upper bound on end - to - end achievable rate @xmath104 where @xmath105 and @xmath106 denote the achievable rates of @xmath107 and @xmath108 channels , respectively . according to ( [ 2 - 6 ] ) , ( [ 2 - 7 ] ) and the assumptions in the above , we have @xmath109 } \right\|}^2 } } \right ] } + 1 } } } \right ) \end{aligned}\ ] ] where @xmath110 is given as @xmath111 with @xmath112{\bf t}^h_s[u]\right]$ ] , @xmath113{\bf t}^h_r[u]\right]$ ] and @xmath114{\bf r}^h_s[u]\right]$ ] denoting the covariance matrices of distortion noises @xmath115 $ ] , @xmath116 $ ] and @xmath117 $ ] , respectively . _ 1 ) achievable rate of @xmath118 channel : _ since @xmath119 is only related with @xmath105 , it can be optimized separately to maximize @xmath105 . from ( [ 3 - 0 - 2 ] ) , the optimization problem is equivalent to the generalized rayleigh quotient problem , which can be solved as @xmath120 . the resultant upper bound on @xmath105 ( denoted by @xmath121 ) can be expressed as @xmath122 where @xmath123 is defined as @xmath124 and is independent of @xmath100 . the second step follows from the matrix inversion lemma @xcite . [ t1 ] assuming that the assumptions * a1 * and * a2 * hold . as @xmath125 , the achievable rate of @xmath118 channel in the large-@xmath30 regime is bounded as @xmath126 where @xmath127 , @xmath128 is determined by the following fixed - point algorithm with @xmath129 @xmath130    see appendix - b . the convergence of the fixed - point algorithm in theorem [ t1 ] has been proved in @xcite . _ 2 ) achievable rate of @xmath108 channel : _ from ( [ 3 - 0 - 2 ] ) , the achievable rates of @xmath118 and @xmath108 channels are coupled through @xmath113{\bf t}^h_r[u]\right]= { \rm diag}\left(\mathbb e\left[{\bf w}_t { \bf \lambda}_r{\bf w}_t^h\right]\right)$ ] , which makes the design of @xmath131 very challenging . however , with the following theorem we show that this coupling disappears in the large-@xmath30 regime , which allows a simple upper bound for achievable rate of @xmath108 channel . [ t2 ] assuming that the assumptions * a1 * and * a2 * hold . as @xmath125 , the achievable rate of @xmath118 channel is independent of @xmath131 . the achievable rate of @xmath108 channel in the large-@xmath30 regime is bounded as ( [ 3 - 3 - 1 ] ) , which is achieved by the eigen bf @xmath132 . see appendix - c . @xmath133    ' '' ''    with theorems [ t1 ] and [ t2 ] , the upper bound of end - to - end achievable in the large-@xmath30 regime can be expressed as @xmath134    [ rm1 ] the bound in ( [ 3 - 4 ] ) is derived by assuming linear processing and fixed transmit powers at sources and relay , which is in general not capacity achieving . one exception is when @xmath135 and @xmath31 . from theorem [ t2 ] , the design of receive and transmit bf matrices is decoupled as @xmath31 . moreover , the distortion noises are assumed to be circularly symmetric complex gaussian distributed and independent of the desired signal . thus , based on the results in @xcite , linear processing along with power control is sufficient to achieve the capacity . when @xmath135 , the terms @xmath136 and @xmath137 in ( [ 3 - 3 ] ) vanish . by neglecting the low order term ( in @xmath32 ) , the bound in ( [ 3 - 4 ] ) reduces to @xmath138 , which is independent of @xmath76 and @xmath139 . thus , ( [ 3 - 4 ] ) is the upper bound for capacity of @xmath140 channel as @xmath135 and @xmath31 . different from remark [ rm1 ] , we consider the general case in which @xmath20 can be either fixed or scales with @xmath32 . with theorem [ t1 ] and theorem [ t2 ] , the effect of hardware impairments is still hard to analysis since @xmath141 is not in closed - form . to make the analysis tractable , we assume @xmath142 @xmath143 . using eigenvalue decomposition @xmath144 on @xmath145 in ( [ 3 - 3 ] ) , we have @xmath146 by examining ( [ 3 - 5 ] ) and the expression of @xmath147 in ( [ 3 - 3 ] ) , @xmath148 is an @xmath149 term . note that this is valid for arbitrary @xmath150 . thus , we can conclude that @xmath151 . replacing @xmath152 in @xmath153 with the symbol @xmath149 and substituting the result into ( [ 3 - 2 ] ) , @xmath105 can be written as ( [ 3 - 6 ] ) . @xmath154    ' '' ''    for convenience , we further assume @xmath155 . note that this assumption does not effect the basis conclusions of the analysis . by neglecting the low order terms ( in @xmath32 ) in ( [ 3 - 6 ] ) and ( [ 3 - 3 - 1 ] ) , the upper bound on end - to - end achievable rate reduces to @xmath156 one can make several observations from ( [ 3 - 7 ] ) :    \1 ) as @xmath157 , the effect of hardware impairments at the relay and ei disappears in the large-@xmath30 regime and the upper bound converges to that in remark [ rm1 ] . the end - to - end achievable rate is limited by distortion noises caused by hardware impairments at sources and destinations . the explanation is that the transmit distortion can be viewed as an interference signal the same channel as data of @xmath22 from ( [ 2 - 5 ] ) . thus , it can not be suppressed after combining by @xmath119 . the power of distortion due to receive imperfection of @xmath23 is proportional to that of desired receive signal ( after eigen bf ) and also can not be reduced by the transmit bf scheme of relay . thus , there is a finite ceiling on end - to - end achievable rate as @xmath31 . in fact , this poses a major limitation on mm - fdr with single antenna sources and destinations , which degrades the gain of massive array of relay greatly . \2 ) as @xmath20 scales with @xmath32 , the effect of hardware impairments at the relay and ei is not negligible . since the ei channel is typically stronger than desired channels , the @xmath107 channel becomes the bottleneck of end - to - end achievable rate . this indicates that it is of great importance to suppress ei to avoid the bottleneck effect . in this section , we consider a general model where sources and destinations are equipped with an arbitrary number of antennas . we first present a low complexity hardware impairments aware transceiver scheme ( referred to as hia scheme ) . then the achievable rate of mm - fdr with proposed hia scheme ( hia - mm - fdr ) is analyzed in the large-@xmath30 regime and the effect of hardware impairments is discussed .      the optimal transceiver scheme for mm - fdr to optimize achievable rate is hard to find since the transmissions of @xmath158 and @xmath159 channels are deeply coupled due to the presence of ei and distortion noises . even without ei and distortion noises , the problem has been proved np - hard @xcite . as @xmath102 is an @xmath160 matrix , the learning of @xmath102 requires training sequence with length no less than @xmath32 . this task is prohibitive as the duration of channel coherent time is limited . thus , we assume that passive ei cancellation @xcite has been used and instantaneous ei channel is not available . a low complexity hia scheme is proposed to mitigate ei ( after passive cancellation ) and distortion noises . as will be shown in simulations , the hia scheme along with passive ei cancellation provides satisfactory gain compared to hdr .    _ 1 ) transceiver design at the relay : _ as the variance of ei channel is typically stronger than desired channels , it is of great importance to control the ei power . from section iv - b , ei suppression is also important in reducing distortion noise due to hardware impairments of relay . we consider a two - stage bf scheme at the relay . the receive and transmit bf matrices are expressed as @xmath161 and @xmath162 .    _ outer bf matrices : _ the outer bf matrices @xmath163 ( @xmath164 ) and @xmath165 ( @xmath166 ) are designed to suppress ei . @xmath167 and @xmath168 are design parameters . due to the lack of @xmath102 , @xmath169 and @xmath170 are allowed to depend only on the statistical knowledge of ei channels . with ( [ 2 - 6 ] ) , we design an average ei power minimization problem , i.e. , @xmath171   \end{aligned}\ ] ] in ( [ 4 - 1 ] ) , the design of @xmath169 and @xmath170 is coupled with the inner bf matrix , which makes closed - form solution inaccessible . thus , it is of great interest to find new target function to decouple the problem . for such , using cauchy - schwarz inequality twice , an upper bound on average ei power can be obtained as @xmath172 then the problem ( [ 4 - 1 ] ) can be rewritten approximately as @xmath173 \\    & \bigg . + {   { \nu _ r}\sum\limits_{l = 1}^k { { e_{r , k } } } \mathbb e\left [ { { \text{tr}}\left ( { { \mathbf{p}}_r^h{{\mathbf{h}}_{ei}}{\mathbf{h}}_{ei}^h{{\mathbf{p}}_r } } \right ) } \right ] } \bigg\ }   \\    & = { \text{arg}}\mathop { { { \min}}}\limits_{{{\mathbf{p}}_r},{{\mathbf{p}}_t } } \left\ { { \text{tr}}\left ( { { { \mathbf{\underline{w}}}_t}{\lambda _ r}{\mathbf{\underline{w}}}_t^h } \right){\text{tr}}\left ( { { \mathbf{p}}_t^h{{{\mathbf{\tilde c}}}_{ei}}{{\mathbf{p}}_t } } \right)\right . \\    & \bigg . \times{\text{tr}}\left ( { { \mathbf{p}}_r^h{{\mathbf{c}}_{ei}}{{\mathbf{p}}_r } } \right)+ { \nu _ r}{n_t}\sum\limits_{l = 1}^k { { e_{r , k } } } { \text{tr}}\left ( { { \mathbf{p}}_r^h{{\mathbf{c}}_{ei}}{{\mathbf{p}}_r } } \right ) \bigg\ } \end{aligned}\ ] ] the multiplicative term @xmath174 has been removed since it has no effect on the optimal solution . the second step follows from a similar derivation of ( [ a3 ] ) . from ( [ 4 - 3 ] ) , the columns of @xmath175 ( @xmath176 ) are composed of the eigenvectors corresponding to the @xmath167 ( @xmath168 ) smallest eigenvalues of @xmath177 ( @xmath178 ) , i.e. , @xmath179 \\ & { \bf \dot{p}}_t = \left [ { { { \mathbf{u}}_{{n_t } - { a_t } + 1}}\left ( { { { { \mathbf{\tilde c}}}_{ei } } } \right),{{\mathbf{u}}_{{n_t } - { a_t } + 2}}\left ( { { { { \mathbf{\tilde c}}}_{ei } } } \right ) , \cdots , { { \mathbf{u}}_{{n_t}}}\left ( { { { { \mathbf{\tilde c}}}_{ei } } } \right ) } \right ] \end{aligned}\ ] ] in ( [ 4 - 3 - 1 ] ) , @xmath167 and @xmath168 can be viewed as parameters to balance the allowable ei power and available degree of freedom ( dof ) for data transmission , which should be optimized with respect to specific metric . this will be considered in section vi .    _ inner bf matrix : _ the inner bf matrices @xmath180 and @xmath181 are designed to realize the multi - user communication . @xmath182 and @xmath183 can be designed with different criteria , e.g. , maximizing the desired signal power which corresponding to the eigen bf , or minimizing the mui which corresponding to the zero - forcing ( zf ) scheme @xcite . we adopt the latter one since the zf scheme is known to approach the asymptotic limit ( in @xmath32 ) of achievable rate faster as the number of relay antennas increases @xcite . for given bf vectors at sources and destinations , define effective channel vectors of @xmath107 and @xmath108 channels as @xmath184 and @xmath185 , and let @xmath186 $ ] and @xmath187 $ ] , the inner bf matrices can be written as @xmath188 where @xmath189 is a diagonal normalized matrix with @xmath190_{l , l}={{\mathbf{e}}_l^h{{\left ( { { \mathbf{\underline{h}}}_{rd}^h{{\mathbf{\underline{h}}}_{rd } } } \right)}^ { - 1}}{{\mathbf{e}}_l}}$ ] . in practice , @xmath191 and @xmath192 should be estimated in order to compute ( [ 4 - 4 ] ) . this will be considered in section v - b . note that ( [ 4 - 4 ] ) requires @xmath193 so that @xmath167 and @xmath168 can be selected to ensure the invertibility of @xmath194 and @xmath195 . _ 2 ) transceiver design at the sources and destinations : _ similar to @xcite , we assume no instantaneous knowledge of channels at sources and destinations . this is reasonable since the amount of feedback could be very huge and unaffordable in mm - fdr . however , it is assumed that the local statistical knowledge of channels ( i.e. , @xmath196 for @xmath22 and @xmath197 for @xmath23 ) can be obtained . as observed from ( [ 3 - 7 ] ) , when @xmath125 , the achievable rates of @xmath118 channel and @xmath108 channel are limited by the transmit distortion noise at @xmath22 and receive distortion noise at @xmath23 . this motivates us to design @xmath75 and @xmath87 to suppress these negative factors with the antenna arrays of @xmath22 and @xmath23 . _ design of @xmath75 : _ intuitively , according to ( [ 2 - 5 ] ) we can design the following problem @xmath198}}{{\mathbb e\left [ { { \rm{tr}}\left ( { { \bf{h}}_{sr , k}{\bf{\theta } } _ { s , k}^\textsf{t}{{\bf{h}}_{sr , k}^h } } \right ) } \right]}}\\ & = \arg \mathop { \max}\limits_{\|{{\bf{p}}_{s , k}}\|=1 } \frac{{{\bf{p}}_{s , k}^h{{{\bf{\tilde c}}}_{sr , k}}{{\bf{p}}_{s , k}}}}{{{\nu _ { s , k}}{\rm{tr}}\left ( { { { { \bf{\tilde c}}}_{sr , k}}{\rm{diag}}\left ( { { \bf{p}}_{s , k}^h{{\bf{p}}_{s , k } } } \right ) } \right ) } } \\ & = \arg \mathop { \max}\limits_{\|{{\bf{p}}_{s , k}}\|=1 } \frac{1}{{{\nu _ { s , k}}}}{\bf{p}}_{s , k}^h{{{\bf{\tilde c}}}_{sr , k}}{{\bf{p}}_{s , k } } \end{aligned}\ ] ] where @xmath199 is the covariance matrix of @xmath200 , i.e. , @xmath201{\bf t}^h_{s , k}[u]]= \nu_{s , k}{\rm diag}\left({\bf p}_{s , k } { \bf p}^h_{s , k}\right)$ ] . the second step follows from a similar derivation with that in ( [ a3 ] ) and the last step is based on assumption * a2*. from ( [ 2 - 5 ] ) , the right - hand side of ( [ 4 - 5 - 0 ] ) can be interpreted as the average signal to transmit distortion noise ratio at @xmath22 . the solution of problem ( [ 4 - 5 - 0 ] ) is @xmath202    _ design of @xmath87 : _ similarly , based on ( [ 2 - 5 ] ) , the received bf vector at @xmath23 can be derived by solving the problem @xmath203}}{{{\mu _ { d , k}}{\mathbf{p}}_{d , k}^h{\text{diag}}\left ( { \mathbb{e}\left [ { { \mathbf{h}}_{rd , k}^h{{\mathbf{h}}_{rd , k } } } \right ] } \right){{\mathbf{p}}_{d , k}}}}$ ] , which results in @xmath204      during a training phase , each source / destination transmits pilot sequence sequentially which allows the relay to compute the estimates of channels . with hia scheme , it is sufficient to estimate the effective channel vectors @xmath205 and @xmath206 . this allows a reduced dimension estimation scheme . for brevity , we consider the estimation of @xmath205 . let @xmath207 and @xmath208 denote the pilot sequence and power of each pilot symbol . @xmath209 is multiplied by @xmath210 and transmitted by @xmath22 . the received pilot matrix at the relay is expressed as @xmath211 where @xmath212 denotes the transmit distortion noise of @xmath22 , whose @xmath13th column has distribution ( based on model ( [ 2 - 3 ] ) ) @xmath213 . @xmath214 denotes receive distortion noise at the relay . with model ( [ 2 - 4 ] ) , the @xmath13th column of @xmath215 has distribution @xmath216 ) } ) $ ] . @xmath217 is the awgn matrix , whose elements are i.i.d and distributed as @xmath42 . by multiplying both side of ( [ 4 - 5 ] ) with @xmath218 , we have @xmath219 from ( [ 4 - 6 ] ) , as a merit of hia scheme , the total length of pilot sequences to obtain the estimates of all effective channels can be as less as @xmath220 ( when @xmath98 is set to 1 ) . [ t3 ] the lmmes estimator of effective channel @xmath205 can be expressed as @xmath221 where @xmath222 is the covariance matrix of @xmath205 , which can be expressed as @xmath223 . @xmath224 is given by @xmath225 and @xmath226    the real effective channel @xmath227 can be decomposed as @xmath228 with @xmath229 denoting the estimation error . @xmath230 and @xmath231 are uncorrelated whose covariance matrices can be expressed as @xmath232 and @xmath233 , respectively . a sketch of the proof is presented in appendix - d . different from that with ideal hardware @xcite , the lmmse estimator is not equivalent to mmse estimator , since the received pilot signal is corrupted by the term @xmath234 ( due to transmit imperfection of @xmath22 ) , which is not independent with the channel to be estimated . there might exist non - linear estimator that results in smaller mse . however , the difference should be small since the distortion noises are relatively weak .    the lmmse estimator of @xmath206 can be obtained using the same approach . for further analysis , we define @xmath235 and @xmath236 as the covariance matrix and lmmse estimates of @xmath206 , respectively . moreover , we define @xmath237 as the covariance matrix of @xmath236 , where @xmath238 is given by @xmath239      this subsection analyzes the achievable rate of hia - mm - fdr . the achievable rate expressions are derived based on the bounding technique in @xcite .    _ 1 ) achievable rate of @xmath240 channel : _ by treating @xmath241}s_k[u]$ ] as the desired signal at @xmath242 channel , and approximating the remaining terms in ( [ 2 - 6 ] ) , i.e. , @xmath243 \!-\!\sqrt{e_{s , k}}{\mathbb e [ { { \bf{w}}_{r , k}^h{{\bf{h}}_{sr , k}}{{\bf{{\dot{p}}}}_{s , k } } } ] } s_k[u]$ ] , using the worst - case uncorrelated additive gaussian noise with the same variance , the rate in ( [ 4 - 11 ] ) is achievable on @xmath244 channel , where @xmath245 , @xmath246 , @xmath247 and @xmath248 denote the ei , mui , effective distortion due to transmit imperfection of sources and effective distortion due to receive imperfection of relay respectively , which are given by    @xmath249 } \right|}^2}}}{{{e_{s , k}}{\rm{var}}\left ( { { \bf{w}}_{r , k}^h{{\bf{h}}_{sr , k}}{{\bf{{\dot{p}}}}_{s , k } } } \right){+ } \mathbb e\left [ { { { \textsf{ei}}}{{{}}_k } } \right ] + \mathbb e\left [ { { \rm{\textsf{mui}}}{{{}}_{r , k } } } \right ] + \mathbb e\left [ { { { \textsf{d}}}_{r , k}^{{\textsf{t } } } } \right ] + \mathbb e\left [ { { { \textsf{d}}}_{r , k}^{{\textsf{r } } } } \right ] + \mathbb e\left [ { { { \left\| { { \bf{w}}_{r , k } } \right\|}^2 } } \right ] } } } \right)\ ] ]    ' '' ''    @xmath250{\mathbf{t}}_r^h\left [ u \right ] } \right ] } \right){\bf{h}}_{ei}^h{{\bf{w}}_{r , k } } \\ { { \textsf{mui}}_{r , k } } = \sum\limits_{j = 1,j \ne k}^k { { e_{s , j}}{{\left| { { \bf{w}}_{r , k}^h{{\bf{h}}_{sr , j}}{{\bf{{\dot{p}}}}_{s , j } } } \right|}^2 } } \\ \textsf{d}_{r , k}^\textsf{t } = \sum\limits_{j = 1}^k { { \bf{w}}_{r , k}^h{{\bf{h}}_{sr , j } } \mathbb   e\left [ { { { \bf{t}}_{s , j}}[u]{\bf{t}}_{s , j}^h[u ] } \right]{\bf{h}}_{sr , j}^h{{\bf{w}}_{r , k } } } \\ \textsf{d}_{r , k}^\textsf{r } = { \bf{w}}_{r , k}^h \mathbb e\left [ { { { \mathbf{r}}_r}\left [ u \right]{\mathbf{r}}_r^h\left [ u \right ] } \right]{{\bf{w}}_{r , k } } \end{array}\ ] ] wherein , the bf matrices at the relay are given by @xmath251 where @xmath252 $ ] ( @xmath93 ) and @xmath253 is a diagonal normalized matrix with @xmath254_{l , l}={{\mathbf{e}}_l^h{{\left ( { { \mathbf{\underline{\hat{h}}}}_{rd}^h{{\mathbf{\underline{\hat{h}}}}_{rd } } } \right)}^ { - 1}}{{\mathbf{e}}_l}}$ ] . note that we have replaced the effective channels in inner bf matrices with their estimates obtained in theorem [ t3 ] . [ t4 ] assume that the assumptions * a1 * and * a2 * hold , and @xmath167 is selected so that @xmath255 . with hia scheme , the achievable rate of @xmath118 channel in the large-@xmath30 regime is given by @xmath256 , where the powers of desired signal and awgn are @xmath257 } \right|}^2 } } } = e_{s , k}$ ] and @xmath258}= { ( { { \rm{tr } } ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k } } } ) } ) ^ { - 1 } } $ ] , and @xmath259 = \frac{{\beta _ { ei}}{\rm{tr}}\left ( { { { { \bf{\tilde c}}}_{ei}}{{\bf \omega } _ r } } \right){\rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k}}{\bf{{\dot{p}}}}_r^h{{\bf{c}}_{ei}}{{\bf{{\dot{p}}}}_r } } \right)}{\left ( { { \rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } \right)^2 } \\ & \mathbb e[{{\textsf{mui}}_{r , k } } ] = \sum_{j = 1,j \ne k}^k \frac{{e_{s , j}}{\rm{tr}}\left ( { \left ( { { { \bf{\underline{c}}}_{sr , j } } - { { { \bf{\hat { \underline{c}}}}}_{sr , j } } } \right){{{\bf{\hat { \underline{c}}}}}_{sr , k } } } \right)}{\left ( { { \rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } \right)^ { 2}}\\ & \mathbb e[\textsf{d}_{r , k}^\textsf{t } ] = \frac{{\nu _ { s , k}}{e_{s , k } } g_{kk } + \sum\limits_{j = 1,j \ne k}^k { { \nu _ { s , j}}{e_{s , j } } g_{kj } { { } } } } { \left ( { { \rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } \right)^ {   2}}\\ & \mathbb e \left[\textsf{d}_{r , k}^\textsf{r } \right ] = \frac{\mu _ r } { { { \rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } } \\ & \times { { \sum\limits_{j = 1}^k { { e_{s , j}}{\beta _ { sr , j}}{\rm{tr}}\left ( { { { { \bf{\tilde c}}}_{sr , j}}{{\bf{\omega } } _ { s , j } } } \right ) }   + { \beta _ { ei}}{\rm{tr}}\left ( { { { { \bf{\tilde c}}}_{ei}}{{\bf{\omega } } _ r } } \right ) + 1 } } \end{aligned}\ ] ] where @xmath260 is defined as @xmath261 $ ] and is given by lemma [ lm1 ] in appendix - a . @xmath262 and @xmath263 are given by ( [ d10 ] ) and ( [ d12 ] ) in the appendix - e . @xmath264 and @xmath265 are expressed as @xmath266 the scaling behaviors of above factors are shown in table [ tab1 ] . see appendix - e .        [ cols= " < , < , < " , ]     several important results on the effect of ei , mui and hardware impairments on achievable rate of @xmath118 channel can be obtained from theorem [ t4 ] and table [ tab1 ] . _ effect of ei : _ with the proposed hia scheme , the scaling behavior of ei is high correlated with @xmath267 and @xmath268 . with the assumption * a1 * , we have @xmath269 and @xmath270 . in fact , under specific spatial correlation model ( e.g. , the physical channel model with a fixed number of angular bins @xcite ) , it is possible to make @xmath271 and @xmath272 by selecting proper @xmath167 and @xmath168 . in this case , the power of ei can decrease faster than @xmath273 in the large-@xmath30 regime . note that this is in contrast with the result in @xcite , which shown that the power of ei decreases exactly with @xmath273 . _ effect of mui : _ the scaling behavior of mui confirms well with the classic result in mm - fdr with perfect hardware @xcite , i.e. , the mui diminishes as @xmath274 if @xmath275 . _ effect of transmit imperfection of sources : _ ( i ) when @xmath25 is fixed , an interesting observation is that the term @xmath276 , which is widely convinced @xmath277 terms with perfect hardware , scales as @xmath278 . as a result , the achievable rate of @xmath118 channel is limited by the joint effect of @xmath276 and the effective distortion due to the transmit imperfection of sources , i.e. , @xmath247 . in fact , the variation of scaling behavior of @xmath276 can also be viewed as a `` negative effect '' of transmit imperfection of sources . to see this , inserting the expression of @xmath260 ( given by lemma [ lm1 ] in appendix - a ) into ( [ 4 - 13 ] ) and letting @xmath274 , we can obtain @xmath279 , which demonstrates that @xmath280 scales as @xmath281 as long as @xmath282 . ( ii ) when @xmath25 scales with @xmath28 and @xmath283 , it is seen that the proposed hia scheme reduces the effect of transmit imperfection at sources by a factor of @xmath284 . the reason is that the transmit distortion is pre - suppressed by the transmit bf vector @xmath285 . _ effect of receive imperfection of the relay : _ by comparing the result in table [ tab1 ] and @xmath286 , it is seen that the hia scheme does not change the scaling behavior of distortion noise caused by receive imperfection of the relay . however , the receive distortion noise power is still reduced , since ei ( which is a main contributor of receive distortion noise ) is suppressed . _ 2 ) achievable rate of @xmath108 channel : _ according to ( [ 2 - 7 ] ) and the bounding technique in @xcite , the achievable rate of @xmath287 channel can be expressed as @xmath288 } \right|}^2}}}{{{e_{r , k}}{\mathop{\rm var } } \left ( { { \bf{{\dot{p}}}}_{d , k}^h{\bf{h}}_{rd , k}^h{{\bf{w}}_{t , k } } } \right ) + \mathbb e\left [ { { \textsf{mui}}_{d , k } } \right ] + \mathbb e\left [ { \textsf{d}_{d , k}^\textsf{t } } \right ] + \textsf{d}_{d , k}^\textsf{r}}+1 } } \right)\end{array }   \end{aligned}\ ] ] where @xmath289 , @xmath290 and @xmath291 denote the mui , effective distortion due to transmit imperfection of relay and effective distortion due to receive imperfection of destinations , respectively , which are given by @xmath292    [ t5 ] assume that assumptions * a1 * and * a2 * hold , and @xmath168 is selected so that @xmath293 . with hia scheme , the achievable rate of @xmath108 channel in the large-@xmath30 regime is given by @xmath294 , where the power of desired signal is @xmath295 } \right|^2 } = { e_{r , k } } { { { \rm{tr } } ( { { \bf{\hat { \underline{c}}}}_{rd , k}^ { } } ) } } $ ] , and @xmath296 = \sum\limits_{j = 1,j \ne k}^k \frac{{e_{r , j}}{\rm{tr}}\left ( { \left ( { { { \bf{\underline{c}}}_{rd , k } } - { { { \bf{\hat { \underline{c}}}}}_{rd , k } } } \right){{{\bf{\hat { \underline{c}}}}}_{rd , j } } } \right ) } { { { { \rm{tr}}\left ( { { \bf{\hat { \underline{c}}}}_{rd , j}^ { } } \right ) } } } \\ & \mathbb e\left [ \textsf{d}_{d , k}^\textsf{t } \right ] = { \nu _ r}\sum\limits_{l = 1}^k \frac{{e_{r , l}}{\rm{tr}}\left ( { { { \bf{c}}_{rd , k}}{\rm{diag}}\left ( { { { \bf{{\dot{p}}}}_t}{{{\bf{\hat { \underline{c}}}}}_{rd , l}}{\bf{{\dot{p}}}}_t^h } \right ) } \right)}{\left({\lambda _ 1}\left ( { { { { \bf{\tilde c}}}_{rd , k } } } \right)\right)^{-1 } { { { \rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{rd , l } } } \right ) } } } \\ & \textsf{d}_{d , k}^\textsf{r }   = \frac{{{\mu _ { d , k}}{e_{r , k}}}}{{{\text{tr}}\left ( { { { { \mathbf{\hat { \underline{c}}}}}_{rd , k } } } \right)}}\left ( { \frac{{{\mathbf{{\dot{p}}}}_{d , k}^h{\text{diag}}\left ( { { { \mathbf{{\dot{p}}}}_{d , k}}{\mathbf{{\dot{p}}}}_{d , k}^h } \right){{\mathbf{{\dot{p}}}}_{d , k}}}}{{{{\left ( { { \text{tr}}\left ( { { { { \mathbf{\hat { \underline{c}}}}}_{rd , k } } } \right ) } \right)}^ { - 2 } } } } } \right . \\ &   +   \frac {   \frac{{{\mu _ r}}}{\tau } \left ( { \frac{1}{{{e_t } } } { + \frac{1}{{{\mu _ r}{e_t}}}}+ { \beta _ { rd , k}}\left ( { { \lambda _ 1}\left ( { { { { \mathbf{\tilde c}}}_{rd , k } } } \right ) + { \nu _ { d , k } } } \right ) } \right )     } { \left({{\text{tr}}\left ( { { \mathbf{\underline{c}}}_{rd , k}^3{\mathbf{\gamma } } _ { rd , k}^2 } \right)}\right)^{-1}{{\lambda _ 1}\left ( { { { { \mathbf{\tilde c}}}_{rd , k } } } \right ) } } \\ & \left . { + \sum\limits_{j = 1,j \ne k}^k { \frac{{{e_{r , j}}{\beta _ { rd , j}}}}{{{e_{r , k}}}}{\text{tr}}\left ( { { { \mathbf{c}}_{rd , k}}{{\mathbf{{\dot{p}}}}_t}{{{\mathbf{\hat { \underline{c}}}}}_{rd , j}}{\mathbf{{\dot{p}}}}_t^h } \right ) }   + \frac{{{\text{tr}}\left ( { { { { \mathbf{\hat { \underline{c}}}}}_{rd , k } } } \right)}}{{{e_{r , k } } } } } \right ) \end{aligned}\ ] ] @xmath297 is defined as @xmath298 $ ] and is given by lemma [ lm1 ] in appendix - a . the scaling behaviors of above factors are shown in table [ tab2 ] . the proof is similar with that for theorem [ t4 ] , and thus it is neglected .    from theorem [ t5 ] and table [ tab2 ] , one can observe a strong similarity between the effects of hardware impairments on @xmath287 channel and that on @xmath118 channel . in particular , the effective distortion due to the transmit imperfection of relay @xmath290 is suppressed by the large transmit array of the relay , and it diminishes as @xmath299 and @xmath300 . moreover , when @xmath26 is fixed , the achievable rate of @xmath287 channel is limited by the term @xmath301 and distortion noise caused by receive imperfection of @xmath23 . however , as @xmath26 scales with @xmath29 and @xmath300 , these factors are suppressed by @xmath302 . this is because that applying @xmath303 can be viewed as a post - suppression on receive distortion noise .    with hia scheme , the following end - to - end rate for @xmath304 channel is achievable in the large-@xmath30 regime @xmath305 different from that for upper bound in the last section , @xmath306 can not be expressed in a simple form like ( [ 3 - 7 ] ) . however , based on table [ tab1 ] and table [ tab2 ] , one can simply deduce that @xmath307 will converge to a finite ceiling if @xmath25 and @xmath26 are fixed . however , with hia scheme , the end - to - end achievable rate grows without bound as @xmath308 if @xmath25 and @xmath26 scale with @xmath309 and @xmath310 . this is the same to the situation with ideal hardware @xcite . based on the achievable rate expressions , we discuss the hardware design of sources , destinations and relay in this subsection . although the achievable rate expressions hold for arbitrary @xmath311 , we will restrict our analysis to @xmath157 . this is a fundamental condition to ensure the benefit of massive mimo , i.e. , the mui is suppressed by a large surplus of degrees of freedom @xcite . in mm - fdr with hardware impairments , it is also essential to control the number of s - d pairs ( by , e.g. , some user scheduling scheme ) to suppress the ei and distortions .      as @xmath25 and @xmath26 are fixed and small ( due to e.g. , size limitation ) , from theorem [ t4 ] and [ t5 ] , it is of great importance to reduce @xmath312 , @xmath313 and @xmath314 in order to alleviate the `` ceiling effect '' on achievable rate . this indicates that it is beneficial to use high - quality hardware at sources and destinations in this case . note that , different from the relay side , the use of high - quality hardware at sources and destinations is affordable , since the increased cost scales with @xmath25 and @xmath26 , instead of @xmath32 . when @xmath25 and @xmath26 scale with @xmath32 , from theorem [ t4 ] and [ t5 ] , the distortion noise powers due to hardware impairments at sources and destinations are lower order terms when compared with desired signal powers . this means that the sources and destinations can decrease the hardware quality to some extent without degrading the performance . later an example will be given to show how to achieve this . similar to that for sources and destinations when @xmath25 and @xmath26 scale with @xmath32 , the distortion noises caused by hardware impairments at the relay are lower order terms when compared with desired signal powers , which makes it possible to decrease the hardware quality without hurting the performance greatly . _ example : _ to get a clear insight , we consider the achievable rate of @xmath315 channel given by theorem [ t4 ] . we assume that there is only one s - d pair and let @xmath316 . when @xmath25 scales with @xmath28 , the achievable rate of @xmath107 channel reduces to @xmath317 where @xmath318 and @xmath319 are positive constants independent of @xmath320 , @xmath321 and @xmath322 . in ( [ 4 - 19 ] ) , we have replaced @xmath25 with the product of @xmath28 and some constant . assume we wish to achieve rate @xmath323 . according to ( [ 4 - 19 ] ) , we have ( as @xmath324 ) @xmath325 @xmath326 . the result is encouraging since it implies that we can increase @xmath327 and @xmath328 linearly as @xmath28 increase ( with @xmath25 scaling with @xmath28 ) without degrading the achievable rate . as a result , inexpensive mm - fdr is possible . in this section , we propose a low - complexity jdpo algorithm to maximize the se ( defined as the sum of all destinations achievable rates ) of hia - mm - fdr , subject to the maximum power constrains . the achievable rate expressions obtained in section v - c are utilized in the proposed algorithm . the algorithm needs only statistical knowledge of channels . therefore , it can be computed offline at a central node ( e.g. , relay ) and then broadcasts to the other nodes . let @xmath329 and @xmath330 be the maximum transmit power constraints at @xmath22 and relay , respectively , the se optimization problem can be formulated as follows @xmath331 where @xmath332 and @xmath333 . according to theorem 4 and theorem 5 , we rewrite @xmath334 and @xmath335 as @xmath336 and @xmath337 , where @xmath338 and @xmath339 denote the effective received sinrs of @xmath107 and @xmath108 channels respectively , which can be expressed as @xmath340 where @xmath341 @xmath342    with the above results and the formula @xmath343 @xmath344 ( @xmath345 is defined as @xmath346 ) , ( [ 5 - 2 ] ) can be rewritten as @xmath347 by using some algebra manipulations on the inequality constraints , @xmath348 , @xmath349 can be rewritten as @xmath350    problem ( [ 5 - 4 ] ) is a combinatorial optimization problem which is in general np hard . to solve ( [ 5 - 4 ] ) , we propose a jdpo algorithm to find suboptimal solution . our strategy is as follows . first , using the similar approach as that in @xcite , we show that the power control problem with fixed @xmath167 and @xmath168 can be approximated as a geometric programming ( gp ) and solved efficiently using the convex optimization tools @xcite . then we present a heuristic approach based on sequential optimization to solve @xmath167 and @xmath168 .    _ power control with fixed @xmath167 and @xmath168 : _ from ( [ 5 - 4 ] ) and ( [ 5 - 5 ] ) , the inequality constraints are posynomial functions @xcite . as shown in @xcite , for any @xmath351 , @xmath352 can be approximated using a monomial function @xmath353 near a point @xmath354 , where @xmath355 and @xmath356 . by using this , the target function of ( [ 5 - 4 ] ) can be approximated as a monomial function , i.e. , @xmath357 . as a result , the solution for power control problem ( with fixed @xmath167 and @xmath168 ) can be obtained by solving several gps . the algorithm is summarized in algorithm [ alg1 ] .    . * compute @xmath358 and @xmath359 . * solve the following gp problem + l _ e_s,1 , , e_s , k;e_r,1 , , e_r , k ; _ k = 1^k _ k^ - 1_k^ - _ k s.t . c_1 ~c_5 if @xmath360 or @xmath361 , stop . otherwise , set @xmath362 and go back to step 2 . [ alg1 ]    _ optimization of @xmath167 and @xmath168 : _ an obvious approach to obtain @xmath167 and @xmath168 is via 2-d search . however , the complexity is very high since the power control problem described in the above should be computed for each @xmath363 . in this work , we present a heuristic approach to find suboptimal solution . in this approach , @xmath167 ( @xmath168 ) is searched over a subset @xmath364 ( @xmath365 ) , whose elements are sampled from @xmath366 ( @xmath367 ) . for example , @xmath364 can be selected as @xmath368 . the motivation of this idea is that the target function does not change dramatically as @xmath167 and @xmath168 increase / decrease by a small step as observed from the numerical results . moreover , instead of 2-d search , @xmath366 and @xmath367 are optimized sequentially over the sampled subsets @xmath364 and @xmath365 via 1-d search . the process repeats several times ( much less than @xmath369 ) to improve the solution . a summary of the jdpo algorithm is presented in algorithm [ alg2 ] . the algorithm converges to a local optimum since the target function is improved in each iteration . the complexity of algorithm [ alg2 ] is upper bounded as @xmath370 , where @xmath371 is the complexity to solve gp when @xmath372 and @xmath373 . usually , gp is solved using inner point method with polynomial time . the exact expression of @xmath371 is quite difficult and related with the structure of the problem . some insights on @xmath371 can be found in [ 37 , sec . 11.5 ] . * initialization : * set an initial @xmath168 , i.e. , @xmath374 , and a maximum repeat times @xmath375 . select the subsets @xmath364 and @xmath365 . . update @xmath376 . [ alg2 ]    [ rm2 ] _ ( jdpo for ee optimization ) _ instead of the se optimization , we can also formulate an ee ( defined as the se divided by total transmit power @xcite ) optimization problem subject to a target se @xmath377 @xmath378 using ( [ 5 - 3 ] ) , the above problem is equivalent to @xmath379 @xmath380 and @xmath381 . with the technique in @xcite , the equality constraint can be converted to a monomial , and the power control problem ( with fixed @xmath167 and @xmath168 ) becomes a gp . thus , ( [ 5 - 6 ] ) can be solving by using a similar jdpo algorithm as that in algorithm [ alg2 ] . this section presents the simulation results to verify the analyses in the previous sections . throughout this section , we set @xmath382 , @xmath383 , @xmath384 for convenience . the correlation matrices of desired channels are generated with the exponential correlation model @xcite @xmath385_{l , j } } = \left\ { { \begin{array}{*{20}{l } } { r_{i , k}^{j - l},l \le j}\\ { { { \left ( { r_{i , k}^{l - j } } \right)}^*},l > j } \end{array } } \right.,i \in \left\ { { sr , rd } \right\}\ ] ] the model approximates the property of ula , where the correlation between adjacent antennas is @xmath386 $ ] and the phase of @xmath387 describes the angle of arrival / departure as seen from the array . @xcite shows how to map some of the parameters of ula to this model . the correlation matrices of ei channel @xmath388 and @xmath389 are generated similarly with a parameter @xmath390 . for convenience , we let @xmath391 , @xmath392 . the phases of @xmath387 and @xmath390 are uniformly selected from @xmath393 $ ] . we assume that 25@xmath39435db ei cancellation can be provided by passive ei suppression techniques ( more than 40db cancellation has been reported by using such techniques for infrastructure node @xcite ) . the variances of ei channel ( after passive cancellation ) and desired channels are selected as @xmath395$]db . with the path loss model in @xcite , the above range corresponds to the setup with the distances from sources / destinations to relay varying from 250 m to 500 m and 10 m segregation between relay transmit and receive arrays . this subsection considers the effect of hardware impairments on se of mm - fdr . the channel coherent time is set to @xmath396 and the length of pilot sequence is @xmath397 .    , @xmath398 , @xmath399db , @xmath400db , @xmath401 , @xmath402 , @xmath403 , @xmath404 , @xmath405.,width=321 ]    .,width=321 ]    fig . [ fig1a ] shows the se of mm - fdr with single antenna at sources / destinations with different levels of hardware impairments . the ses based on transceiver scheme in section iv ( which achieves the upper bound when @xmath125 ) and hia scheme ( with @xmath406 ) are simulated . from fig . [ fig1a ] , the se is more sensitive to hardware impairments at sources and destinations . when @xmath407 , the se approaches to a finite ceiling quickly as the number of relay antennas increases . similar results can be observed when sources and destinations are equipped with multiple but fixed number of antennas in fig . [ fig1b ] . however , the result changes when @xmath25 and @xmath26 scale with @xmath32 and hia scheme is used . from fig . [ fig1b ] , as @xmath408 and @xmath409 , similar ses are achieved as sources / destinations or relay employ low - quality hardware , and no performance ceiling appears . this demonstrates the validness of hia scheme . at last , it is seen that the asymptotic results in theorem [ t1 ] & [ t2 ] and theorem [ t4 ] & [ t5 ] match well with exact results .    , @xmath398 , @xmath410db , @xmath400db , @xmath411 , @xmath412db , @xmath402 , @xmath403 , @xmath404 , @xmath405 , @xmath408 , @xmath409.,width=321 ]    as the performance is affected by impairments of both transmit and receive rf chains at the relay , we compare the effect of transmit and receive imperfections on se in fig . we assume that the sources and destinations use high - quality hardware ( @xmath413 ) . from the figure , the effect of receive imperfection is more detrimental . the reason is that , with ei at the relay , the power of distortion noise caused by receive imperfection is much stronger than that due to transmit imperfection . however , the performance difference decreases as @xmath28 increases from @xmath414 to @xmath415 , since the power of effective receive distortion scales as @xmath273 ( see table [ tab1 ] ) . this implies that use relatively higher - quality hardware or more antennas at the receive side of relay is beneficial . using a similar setup , one can obtain a parallel conclusion for destination , i.e. , the receive imperfection is more harmful than transmit imperfection . this is because that the transmit imperfection of destination only induces larger channel estimator errors , which is a part of the received signal at the destination . thus , based on model ( [ 2 - 4 ] ) , the distortion noise due to receive imperfection will be more detrimental .      , where @xmath416 , @xmath417 , @xmath410db , @xmath400db , @xmath411 , @xmath412db , @xmath418 , @xmath419 , @xmath420 , @xmath421 , @xmath404 , @xmath405.,width=321 ]    fig . [ fig4 ] shows the se as a function of the number of s - d pairs @xmath20 , where the length of pilot sequence of each source or destination is set to 1 . when @xmath396 , the figure reveals that the se will not increase without bound as @xmath20 increases since se is ultimately limited by the channel coherent time . in fact , as the number of s - d pair approaches to @xmath422 , the se converges to zero since all time is allocated to pilot phase . meanwhile , the se per s - d pair is a strictly decreasing function of @xmath20 . similar result can be observed as the channel coherent time is very long ( @xmath423 ) . this is because that the mui , ei and distortion noises become limiting factors when @xmath20 is large as seen from the scaling behaviors in table [ tab1 ] and table [ tab2 ] . this implies that the number of s - d pairs should be limited in order to ensure a se guarantee for each s - d pair , and this number decreases as the quality of hardware degrades .      in this subsection , the se of hia - mm - fdr is compared with the massive mimo hdr ( mm - hdr ) @xcite and mm - fdr with zf - based transceiver ( zf - mm - fdr ) @xcite . the channel coherent time is set to @xmath396 and the length of pilot sequence is @xmath397 . the power of pilot symbol is @xmath400db . for hia - mm - fdr , the maximum power constraints at sources and relay are set to @xmath424db and @xmath425 . without jdpo , we set @xmath426 , @xmath427db . moreover , @xmath167 and @xmath168 are set to @xmath404 and @xmath428 . when jdpo is applied , @xmath76 , @xmath139 , @xmath167 and @xmath168 are determined by algorithm [ alg2 ] . the maximum repeat time @xmath375 of jdpo is set to 3 and the elements of subset @xmath429 are picked uniformly from @xmath430 with step @xmath431 . @xmath432 is obtained with the similar approach . the above parameters are chosen so that the performance loss due to the suboptimal search approach to obtain @xmath167 and @xmath168 is negligible .    , @xmath411 , @xmath412db , @xmath408 , @xmath409 , @xmath418 , @xmath419.,width=321 ]    fig . [ fig4 ] simulates the number of relay antennas required to achieve 3bit / s / hz se per s - d pair ( which ideally can support the 64-qam transmission with 1/2 channel code ) with different levels of hardware quality , where we set @xmath416 and @xmath417 . the figure reveals a tradeoff between the number of antennas and hardware quality , that is ( as discussed in section v - d ) , by increasing @xmath28 and @xmath29 , we can reduce @xmath433 and @xmath434 linearly without degrading the se . meanwhile , it is seen that the proposed hia - mm - fdr reduces the required number of relay antennas significantly when compared with zf - mm - fdr . moreover , the hia - mm - fdr outperforms the mm - hdr for large @xmath433 and @xmath434 . this is because that the distortion noises become the main limiting factor in this case when compared with ei . , @xmath408 , @xmath409 , @xmath411 , @xmath412db , @xmath418 , @xmath419 , @xmath435 , @xmath436.,width=321 ]    the effect of asymmetric numbers of transmit and receive antennas at the relay is shown in fig . [ fig5 ] , where we set @xmath437 . it is shown that allocating more antennas to the receive side of relay is beneficial for hia - mm - fdr ( without jdpo ) and zf - mm - fdr . the reason is that , in the considered setup , the limiting factors of system performance are ei and distortion noise due to receive imperfection at the relay . from theorem [ t4 ] , these factors can be suppressed by receive antenna array of relay . however , when jdpo is applied , it is optimal to set @xmath438 . this is because that the jdpo algorithm in fact tries to balance the achievable rates of each hop by power control and adjusting @xmath167 and @xmath168 . this makes it unnecessary to allocate more antennas to receive side .    , where @xmath398 , @xmath408 , @xmath409 , @xmath411 , @xmath418 , @xmath419 , @xmath435 , @xmath436.,width=321 ]    fig . [ fig6 ] shows the se of hia - mm - fdr as a function of variance for ei channel @xmath439 . as expected , there exists a switching point between hia - mm - fdr and mm - hdr as @xmath439 increases . by increasing @xmath28 and @xmath29 , the constraint on @xmath439 for hia - mm - fdr to achieve a performance gain relaxes , which indicates that hia - mm - fdr becomes more attractive when the number of relay antennas is large . moreover , fig . [ fig6 ] demonstrates that the proposed jdpo algorithm can reduce the constraint on @xmath439 significantly . in particular , as @xmath440 , a se gain of 7.5bit / s / hz can be achieved by hia - mm - fdr when compared to mm - hdr as @xmath439 is @xmath441db .    , @xmath442db , @xmath408 , @xmath409 , @xmath402 , @xmath403 , @xmath435 , @xmath436.,width=321 ]    fig . [ fig7 ] considers the se - ee tradeoff of different schemes . the large - scale fading coefficients of channels are set to @xmath443 which is a realization generated with the model in @xcite . the ees of hia - mm - fdr and zf - mm - fdr are optimized by solving the problem in remark [ rm2 ] with jdpo algorithm and [ 12 , algorithm 1 ] , respectively . it is observed that hia - mm - fdr achieves better se - ee tradeoff when compared with zf - mm - fdr . the gain is mainly due to the optimization of @xmath167 and @xmath168 . this reason is that , if we fix @xmath372 and @xmath373 , the jdpo algorithm for ee optimization is similar to [ 12 , algorithm 1 ] . the only difference is that the power for each steam at the relay , i.e. , @xmath139 , is optimized in jdpo algorithm and [ 12 , algorithm 1 ] optimizes only the total power of relay . this paper considers the transceiver design of mm - fdr with hardware impairments . a low complexity hia scheme is proposed to mitigate the distortion noises by exploiting the statistical knowledge of channels and antenna arrays at sources and destinations . a joint degree of freedom and power optimization algorithm is presented to further optimize the se of hia - mm - fdr . the analytic results demonstrate that the proposed scheme can mitigate the  ceiling effect " appears in traditional mm - fdr protocol , if the numbers of antennas at sources and destinations can scale with that at the relay . moreover , simulation results show that the hia - mm - fdr outperforms mm - fdr with traditional transceiver scheme . [ lm1 ] let @xmath444 and @xmath445 be the lmmse estimates of @xmath446 and @xmath447 , respectively . define @xmath448 $ ] and @xmath449 $ ] . in the large-(@xmath450 ) regime , we have @xmath451 the expression of @xmath452 can be obtained by replacing @xmath453 , @xmath454 , @xmath455 , @xmath312 and @xmath285 in ( [ x1 ] ) with @xmath456 , @xmath457 , @xmath458 , @xmath313 and @xmath303 , respectively . we show the proof for @xmath260 and the derivation for @xmath297 is similar . according to the uncorrlation between @xmath459 and @xmath460 , we have @xmath461 = \mathbb e [ { { { | { { \mathbf{{\underline{h}}}}_{sr , k}^h{{{\mathbf{\hat { \underline{h}}}}}_{sr , k } } } |}^2 } } ] - { ( { { \text{tr } } ( { { { { \mathbf{\hat { \underline{c}}}}}_{sr , k } } } ) } ) ^2}$ ] . with ( [ 4 - 7 ] ) and ( [ 4 - 8 ] ) , @xmath462 $ ] can be written as @xmath463 = \mathbb e\left [ { { { \left| { { \mathbf{{\underline{h}}}}_{sr , k}^h{{\mathbf{{\underline{c}}}}_{sr , k}}{{\mathbf{\gamma } } _ { sr , k}}{{\mathbf{{\underline{h}}}}_{sr , k } } } \right|}^2 } } \right ] \\    & + { \left ( { \frac{1}{{\tau { e_t } } } } \right)^2 } \mathbb e\left [ { { { \left| { { \mathbf{{\underline{h}}}}_{sr , k}^h{{\mathbf{{\underline{c}}}}_{sr , k}}{{\mathbf{\gamma } } _ { sr , k}}{\mathbf{{\dot{p}}}}_r^h{{\mathbf{h}}_{sr , k}}{\mathbf{t}}_{sr , k}^{}\phi _ { } ^ * } \right|}^2 } } \right ] \\    &   + { \left ( { \frac{1}{{\tau { e_t } } } } \right)^2 } \mathbb e\left [ { { { \left| { { \mathbf{{\underline{h}}}}_{sr , k}^h{{\mathbf{{\underline{c}}}}_{sr , k}}{{\mathbf{\gamma } } _ { sr , k}}{\mathbf{{\dot{p}}}}_r^h{{\mathbf{r}}_{sr , k}}\phi _ { } ^ * } \right|}^2 } } \right ] \\    & + { \left ( { \frac{1}{{\tau { e_t } } } } \right)^2}\mathbb e\left [ { { { \left| { { \mathbf{{\underline{h}}}}_{sr , k}^h{{\mathbf{\underline{c}}}_{sr , k}}{{\mathbf{\gamma } } _ { sr , k}}{\mathbf{{\dot{p}}}}_r^h{{\mathbf{n}}_{sr , k}}\phi _ { } ^ * } \right|}^2 } } \right ] \\ \end{aligned}\ ] ] in the large-(@xmath450 ) regime , the result in [ 17 , lemma 4 ( iv ) ] shows that @xmath464 = { \left ( { { \text{tr}}\left ( { { { { \mathbf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } \right)^2}$ ] . replacing @xmath465 with @xmath466 and using the following large @xmath28 approximation ( obtained based on [ 17 , lemma 4 ( ii ) ] ) @xmath467 it can be shown that @xmath468 \\    & = { \nu _ { s , k}}\tau e_t^2\frac{{{\bf{{\dot{p}}}}_{s , k}^h{\bf{\tilde c}}_{sr , k}^{}{\rm{diag}}\left ( { { \bf{{\dot{p}}}}_{s , k}^h{{\bf{{\dot{p}}}}_{s , k } } } \right){\bf{\tilde c}}_{sr , k}^{}{{\bf{{\dot{p}}}}_{s , k}}}}{{\left ( { { \rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } \right)^{-2}}{{{\left ( { { \lambda _ 1}\left ( { { { { \bf{\tilde c}}}_{sr , k } } } \right ) } \right)}^2 } } } \end{aligned}\ ] ] using the similar approach on the remaining terms in ( [ x3 ] ) , lemma [ lm1 ] is obtained . [ lm2 ] let @xmath469 and @xmath470 ( @xmath471 ) be the covariance matrices of effective channel @xmath472 and its estimates @xmath473 given by theorem [ t3 ] . if @xmath474 and @xmath475 , we have @xmath476 and @xmath477 . moreover , @xmath478 and @xmath479 ( i.e. , @xmath167 and @xmath168 scale with @xmath480 and @xmath481 , respectively ) . we present the proof for @xmath482 and @xmath167 . using definition of @xmath483 in theorem [ t3 ] and eigenvalue decomposition @xmath484 , we have @xmath485 where @xmath486 and @xmath487 are positive constants independent of @xmath167 . from ( [ x6 ] ) , @xmath488 is bounded as @xmath489 in ( [ x7 ] ) , the upper bound is an @xmath480 term . with the definition of @xmath222 in the theorem [ t3 ] , we have @xmath490 the last step is based on assumption * a1*. since @xmath491 , we can conclude that @xmath492 ( otherwise @xmath493 ) . thus , the lower bound in ( [ x7 ] ) also scales with @xmath480 . moreover , based on the definition of @xmath222 in the theorem [ t3 ] , we have @xmath494 where @xmath495 is the @xmath13th column of @xmath496 . from ( [ x9 ] ) , @xmath167 must scale @xmath480 . otherwise , @xmath497 based on assumption * a1*.      to facilitate analysis , we first derive the covariance matrices of distortion noises @xmath115 $ ] and @xmath117 $ ] . based on the model ( [ 2 - 3 ] ) , we have @xmath498 . moreover , define @xmath499 = { \left [ { { \mathbf{x}}_{s,1}^t\left [ u \right ] , \cdots , { \mathbf{x}}_{s , k}^t\left [ u \right ] } \right]^t}$ ] , @xmath500 can be expressed based on model ( [ 2 - 4 ] ) as @xmath501 } \right ) \\    & + { \mu _ r}{\text{diag}}\left ( { \mathbb{e}\left [ { { { \mathbf{h}}_{ei}}{{\mathbf{\omega } } _ r}{\mathbf{h}}_{ei}^h } \right ] } \right ) + { \mu _ r}{{\mathbf{i}}_{{n_r } } } \end{aligned}\ ] ] where @xmath502 + { { \bf{t}}_s}\left [ u \right ] } \right){{\left ( { { \bf{s}}\left [ u \right ] + { { \bf{t}}_s}\left [ u \right ] } \right)}^h } } \right ] \\ \hspace{1.7em } = { \rm { diag}}\left ( { \left ( { 1 + { \nu _ { s,1 } } } \right){e_{s,1 } } , \cdots , \left ( { 1 + { \nu _ { s , k } } } \right){e_{s , k } } } \right)\\ { { \bf{\omega } } _ r } = \mathbb{e}\left [ { \left ( { { { \bf{x}}_r}\left [ u \right ] + { { \bf{t}}_r}\left [ u \right ] } \right){{\left ( { { { \bf{x}}_r}\left [ u \right ] + { { \bf{t}}_r}\left [ u \right ] } \right)}^h } } \right ] \\ \hspace{1.7em } = \mathbb{e}\left [ { { { \mathbf{x}}_r}\left [ u \right]{\mathbf{x}}_r^h\left [ u \right ] } \right ] + { \mathbf{\theta } } _ r^\textsf{t } \end{array}\ ] ] for further analysis , we approximate @xmath503{\mathbf{x}}_r^h\left [ u \right ] } \right]$ ] in ( [ a1 - 1 ] ) with @xmath504{\mathbf{x}}_r^h\left [ u \right ] } \right ] ) = \frac{1}{\nu_r } { \mathbf{\theta } } _ r^\textsf{t}$ ] . note that the approximation results in a new upper bound on @xmath121 ( the new bound will also be referred to as @xmath121 for convenience ) . straight - forward computations yield to @xmath505 & = \sum_{l = 1}^k { { \nu _ { s , l}}{e_{s , l}}\mathbb e\left [ { { { \mathbf{h}}_{sr , l}}{\mathbf{h}}_{sr , l}^h } \right ] }   \\ & = \sum_{l = 1}^k   { { ( 1+\nu _ { s , l})}{e_{s , l}}\beta_{sr , l}{{\mathbf{c}}_{sr , l } } } \end{aligned}\ ] ] following the spatial correlation model of @xmath506 in section iii and the approximation under ( [ a1 - 1 ] ) , the expression of @xmath507}$ ] can be obtained as in ( [ a3 ] ) , where the third step follows from the formula @xmath508 and the fourth step is due to the fact that the elements of @xmath509 are i.i.d with distribution @xmath42 . inserting ( [ a2 ] ) , ( [ a3 ] ) into ( [ a1 ] ) , the expression of @xmath500 in ( [ 3 - 3 ] ) is obtained . @xmath510 & \approx \left(1+\frac{1}{\nu_r}\right){\beta _ { ei } } \mathbb e\left [ { { \mathbf{c}}_{ei}^{1/2}{{\mathbf{x}}_{ei}}{\mathbf{\tilde c}}_{ei}^{1/2}{{\mathbf{\theta } } _ r^\textsf{t}}{\mathbf{\tilde c}}_{ei}^{1/2}{\mathbf{x}}_{ei}^h{\mathbf{c}}_{ei}^{1/2 } } \right ] \hfill \\    & = \left(1+\frac{1}{\nu_r}\right){\beta _ { ei}}\sum\limits_{l = 1}^{{n_t } } { \mathbb{e}\left [ { { \mathbf{c}}_{ei}^{1/2}{{\mathbf{x}}_{ei}}{\mathbf{\tilde c}}_{ei}^{1/2}{{\left ( { { \mathbf{\theta } } _ r^\textsf{t } } \right)}^{1/2}}{{\mathbf{e}}_l}{\mathbf{e}}_l^t{{\left ( { { \mathbf{\theta } } _ r^\textsf{t } } \right)}^{1/2}}{\mathbf{\tilde c}}_{ei}^{1/2}{\mathbf{x}}_{ei}^h{\mathbf{c}}_{ei}^{1/2 } } \right ] }   \hfill \\    & = \left(1+\frac{1}{\nu_r}\right){\beta _ { ei}}\sum\limits_{l = 1}^{{n_t } } { \mathbb{e}\left [ { { { \left\| { { { \left ( { { \mathbf{\tilde c}}_{ei}^{1/2}{{\left ( { { \mathbf{\theta } } _ r^\textsf{t } } \right)}^{1/2}}{{\mathbf{e}}_l } } \right)}^t } \otimes { \mathbf{c}}_{ei}^{1/2}{\text{vec}}\left ( { { { \mathbf{x}}_{ei } } } \right ) } \right\|}^2 } } \right ] }   \hfill \\    & = \left(1+\frac{1}{\nu_r}\right){\beta _ { ei}}\sum\limits_{l = 1}^{{n_t } } { { { \left\| { { { \left ( { { \mathbf{\tilde c}}_{ei}^{1/2}{{\left ( { { \mathbf{\theta } } _ r^\textsf{t } } \right)}^{1/2}}{{\mathbf{e}}_l } } \right)}^t } \otimes { \mathbf{c}}_{ei}^{1/2 } } \right\|}^2}}\hfill = \left(1+\frac{1}{\nu_r}\right){\beta _ { ei}}{\text{tr}}\left ( { { { { \mathbf{\tilde c}}}_{ei}}{\mathbf{\theta } } _ r^\textsf{t } } \right){\mathbf{c}}_{ei }   \hfill \\ \end{aligned}\ ] ]      then we derive the asymptotic bound in theorem [ t1 ] . by replacing @xmath511 in @xmath512 with @xmath513 , we approximate @xmath512 as @xmath514 the approximation results in a new upper bound on @xmath121 , since it reduces the power of residual ei after combining by @xmath119 ( the new bound will also be referred to as @xmath121 for convenience ) . based on [ 17 , lemma 4 ( ii ) ] , we can replace @xmath515 with its deterministic equivalence in the large-@xmath29 regime @xmath516 . finally , the expression of @xmath121 in theorem [ t1 ] can be obtained by first applying [ 17 , lemma 4 ( ii ) ] on @xmath517 , and then using [ 17 , theorem 1 ] . since the scenario of interest is the large-@xmath30 regime , we need only consider the effect of @xmath131 on @xmath518 ( given by theorem [ t1 ] ) , more precisely , the term @xmath519 . with the assumption * a2 * , we have @xmath520 . thus , using the model ( [ 2 - 3 ] ) , we have @xmath521}{\mathbf{x}}_r^h[u ] } \right ] } \right ) } \right ) = { \nu _ r}\sum_{l=1}^k e_{r , l}$ ] , which is independent of @xmath131 .    then we derive the upper bound in theorem [ t2 ] . using the model ( [ 2 - 4 ] ) , the power of received distortion at the relay can be derived as @xmath522 } \right\|}^2 } } \right ] = { \mu _ { d , k}}\left ( { { n_t}{\beta _ { rd , k}}{e_{r , k } } + { \nu _ r}{n_t}{\beta _ { rd , k}}\sum_{l = 1}^k { { e_{r , l } } }   + 1 } \right)$ ] . moreover , in the large @xmath29 regime , the term @xmath523 in ( [ 3 - 0 - 2 ] ) approaches to @xmath524 [ 17 , lemma 4 ( ii ) ] . therefore , it is sufficient to design @xmath131 to maximize the numerator of ( [ 3 - 0 - 2 ] ) , which is exactly the eigen bf scheme . the resultant upper bound on achievable rate can be obtained by applying [ 17 , lemma 4 ( ii ) ] and assumption * a2 * on the numerator of ( [ 3 - 0 - 2 ] ) .      with the expression of lmmse estimator @xcite , we have @xmath525 { ( { \mathbb{e } [ { { { { \bf{\tilde { z}}}}_{sr , k}}{\bf { \tilde{z}}}_{sr , k}^h } ] } ) ^ { - 1}}{{{\bf { \tilde{z}}}}_{sr , k}}$ ] . using the independence between @xmath526 and distortion noises , we have @xmath527 & = \mathbb e\left [ { { \bf{{\dot{p}}}}_r^h{{\bf{h}}_{sr , k}}{{\bf{{\dot{p}}}}_{s , k}}{\bf{{\dot{p}}}}_{s , k}^h{\bf{h}}_{sr , k}^h{{\bf{{\dot{p}}}}_r } } \right]\\ & = { \beta _ { sr , k}}{\lambda _ 1}\left ( { { { { \bf{\tilde c}}}_{sr , k } } } \right){\bf{{\dot{p}}}}_r^h{{\bf{c}}_{sr , k}}{{\bf{{\dot{p}}}}_r } \\   & \buildrel \delta \over = { { \bf{\underline{c}}}_{sr , k } } \end{aligned}\ ] ] where the second step follows from a similar derivation with that in ( [ a3 ] ) . moreover , @xmath528}$ ] can be expressed as @xmath529= { { \bf{\underline{c}}}_{sr , k}}\\ & + { \left ( { \frac{1}{{\tau { e_t } } } } \right)^2}\mathbb{e}\left [ { { \mathbf{{\dot{p}}}}_r^h{{\mathbf{h}}_{sr , k}}{{\mathbf{t}}_{sr , k}}{\phi ^*}{\phi ^t}{\mathbf{t}}_{sr , k}^h{\mathbf{h}}_{sr , k}^h{{\mathbf{{\dot{p}}}}_r } } \right]\\ &   + { \left ( { \frac{1}{{\tau { e_t } } } } \right)^2}\mathbb{e}\left [ { { \mathbf{{\dot{p}}}}_r^h{{\mathbf{r}}_{sr , k}}{\phi ^*}{\phi ^t}{\mathbf{r}}_{sr , k}^h{{\mathbf{{\dot{p}}}}_r } } \right ] \\ & + { \left ( { \frac{1}{{\tau { e_t } } } } \right)^2}\mathbb{e}\left [ { { \mathbf{{\dot{p}}}}_r^h{{\mathbf{n}}_{sr , k}}{\phi ^*}{\phi ^t}{\mathbf{n}}_{sr , k}^h{{\mathbf{{\dot{p}}}}_r } } \right ] \end{aligned}\ ] ] according to the independence between @xmath191 and @xmath530 , the second term of right - hand side of ( [ c2 ] ) can be rewritten as @xmath531 \\ & = { \nu _ { s , k}}\tau e_t^2\mathbb{e}\left [ { { \mathbf{{\dot{p}}}}_r^h{{\mathbf{h}}_{sr , k}}{\text{diag}}\left ( { { { \mathbf{{\dot{p}}}}_{s , k}}{\mathbf{{\dot{p}}}}_{s , k}^h } \right){\mathbf{h}}_{sr , k}^h{{\mathbf{{\dot{p}}}}_r } } \right]\\ &   = { \nu _ { s , k}}\tau e_t^2 \beta_{sr , k } { \rm{tr}}\left ( { { { { \bf{\tilde c}}}_{sr , k}}{\rm{diag}}\left ( { { { \bf{{\dot{p}}}}_{s , k}}{\bf{{\dot{p}}}}_{s , k}^h } \right ) } \right){\bf{{\dot{p}}}}_r^h{{\bf{c}}_{sr , k}}{{\bf{{\dot{p}}}}_r } \end{aligned}\ ] ] using the similar approach on the remaining terms of ( [ c2 ] ) , and substituting the resultant expression of ( [ c2 ] ) and ( [ c1 ] ) into the expression of lmmse estimator , the result in theorem [ t3 ] is obtained . using ( [ 4 - 7])-([4 - 8 ] ) , the second order statistics of @xmath532 and @xmath229 in theorem [ t3 ] can be easily verified . note that @xmath535_{l , j } } = { \bf{\hat { \underline{h}}}}_{i , l}^h{{{\bf{\hat { \underline{h}}}}}_{i , j } } = { \mathbf{x}}_{i , l}^h{\mathbf{\hat { \underline{c}}}}_{i , l}^{1/2}{\mathbf{\hat { \underline{c}}}}_{i , j}^{1/2}{{\mathbf{x}}_{i , j}}$ ] ( @xmath471 ) . based on lemma [ lm2 ] , the dimensions of @xmath536 and @xmath537 ( i.e. , @xmath167 and @xmath168 ) approach to infinity as @xmath538 . thus , according to [ 17 , lemma 4 ( ii ) and ( iii ) ] , in the large-(@xmath450 ) regime , we can deduce @xmath539      _ ( a ) asymptotic expression : _ based on the model ( [ 2 - 3 ] ) , we have @xmath541{\mathbf{t}}_r^h\left [ u \right ] } \right ] } = { \nu _ r}{\rm{diag}}\left ( { \mathbb e\left [ { { { \bf{w}}_t}{\lambda _ r}{\bf{w}}_t^h } \right ] } \right)$ ] . since @xmath44 and @xmath542 are independent , @xmath245 can be rewritten as @xmath543 = \mathbb e [ { \bf{w}}_{r , k}^h{{\bf{h}}_{ei}}{{\bf{\omega } } _ r}{\bf{h}}_{ei}^h{{\bf{w}}_{r , k } } ]   = \mathbb e [ { { { \| { { \bf{e}}_k^t { { ( { { \bf{\hat { \underline{h}}}}_{sr}^h{{{\bf{\hat { \underline{h}}}}}_{sr } } } ) } ^ { - 1}}{{{\bf{\hat { \underline{h}}}}}_{sr}}{{\bf{{\dot{p}}}}_r^h}{{\bf{h}}_{ei}}{\bf{\omega } } _ r^{1/2 } } \|^2 } } } ] $ ] , where @xmath544 is defined as @xmath545 + { { \bf{t}}_r}\left [ u \right ] } \right){{\left ( { { { \bf{x}}_r}\left [ u \right ] + { { \bf{t}}_r}\left [ u \right ] } \right)}^h } } ] = \mathbb e\left [ { { { \bf{w}}_t}{\lambda _ r}{\bf{w}}_t^h } \right ] + { \nu _ r}{\rm{diag}}\left ( { \mathbb e\left [ { { { \bf{w}}_t}{\lambda _ r}{\bf{w}}_t^h } \right ] } \right)$ ] . inserting ( [ 4 - 12 - 1 ] ) into the expression of @xmath544 and using ( [ d1 ] ) , we can obtain @xmath546 in the large-@xmath30 regime , using ( [ d1 ] ) , we have @xmath547 &   = \frac{\mathbb e\left [ { { { { \bf{\hat { \underline{h}}}}}_{sr , k}}{{\bf{{\dot{p}}}}_r^h}{{\bf{h}}_{ei}}{{\bf{\omega } } _ r}{\bf{h}}_{ei}^h{{\bf{{\dot{p}}}}_r}{\bf{\hat { \underline{h}}}}_{sr , k}^h } \right]}{\left ( { { \rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } \right)^ { 2 } }   \\ & = \frac{\mathbb e\left [ { { \rm{tr}}\left ( { { \bf{\hat { \underline{c}}}}_{sr , k}^{1/2}{{\bf{{\dot{p}}}}_r^h}{{\bf{h}}_{ei}}{{\bf{\omega } } _ r}{\bf{h}}_{ei}^h{{\bf{{\dot{p}}}}_r}{\bf{\hat { \underline{c}}}}_{sr , k}^{1/2 } } \right ) } \right]}{\left ( { { \rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } \right)^ {   2 } } \end{aligned}\ ] ] the second step is based on the property @xmath548 and independence between @xmath549 and @xmath44 . finally , using a similar derivation with ( [ a3 ] ) , @xmath543 $ ] in theorem [ t4 ] is obtained .    _ ( b ) scaling behavior : _ substituting ( [ d6 ] ) into the expression of @xmath543 $ ] in theorem [ t4 ] and neglecting the terms that have no effect on the scaling behavior , we have @xmath550 = { \beta _ { ei}}{\left ( { { \text{tr}}\left ( { { { { \mathbf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } \right)^ { - 2}}\sum\limits_{l = 1}^k { { e_{r , l}}{{\left ( { { \text{tr}}\left ( { { { { \mathbf{\hat { \underline{c}}}}}_{rd , l } } } \right ) } \right)}^ { - 1 } } }    \\ \times { \text{tr}}\left ( { { { { \mathbf{\hat { \underline{c}}}}}_{rd , l}}{\mathbf{{\dot{p}}}}_t^h{{{\mathbf{\tilde c}}}_{ei}}{{\mathbf{{\dot{p}}}}_t } } \right){\text{tr}}\left ( { { { { \mathbf{\hat { \underline{c}}}}}_{sr , k}}{\mathbf{{\dot{p}}}}_r^h{{\mathbf{c}}_{ei}}{{\mathbf{{\dot{p}}}}_r } } \right ) \hfill \\     \le { \beta _ { ei}}\sum\limits_{l = 1}^k \frac{{{e_{r , l}}{\lambda _ { { n_r } - { a_r } + 1}}\left ( { { { \mathbf{c}}_{ei } } } \right ) } { \lambda _ { { n_t } - { a_t } + 1}}\left ( { { { { \mathbf{\tilde c}}}_{ei } } } \right ) } { { { \text{tr}}\left ( { { { { \mathbf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } } \end{aligned}\ ] ] where the second step follows from the definition of @xmath496 and @xmath551 . based on lemma [ lm2 ] , @xmath552 . thus we can conclude that @xmath543 \le { \cal o}({\lambda _ { { n_r } - { a_r } + 1}}\left ( { { { \mathbf{c}}_{ei } } } \right ) { \lambda _ { { n_t } - { a_t } + 1 } } ( { { { { \mathbf{\tilde c}}}_{ei } } } ) kn_r^ { - 1})$ ] .      _ ( a ) asymptotic expression : _ based on the model ( [ 2 - 3 ] ) , the covariance matrix of @xmath553 can be derived as @xmath554{\mathbf{t}}_{s , j}^h\left [ u \right ] } \right ] } = { \nu _ { s , j}}{e_{s , j}}{\rm{diag}}\left ( { { { \bf{{\dot{p}}}}_{s , j}}{\bf{{\dot{p}}}}_{s , j}^h } \right)$ ] . thus , in the large-@xmath30 regime , using ( [ d1 ] ) @xmath534 $ ] can be expressed as @xmath555 = { \left ( { { \rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } \right)^ { - 2}}{\nu _ { s , j}}{e_{s , k}}\\    & \times\underbrace { \mathbb e\left [ { { \bf{\hat { \underline{h}}}}_{sr , k}^h{\bf{{\dot{p}}}}_r^h{{\bf{h}}_{sr , k}}{\rm{diag}}\left ( { { { \bf{{\dot{p}}}}_{s , k}}{\bf{{\dot{p}}}}_{s , k}^h } \right){\bf{h}}_{sr , k}^h{{\bf{{\dot{p}}}}_r}{{{\bf{\hat { \underline{h}}}}}_{sr , k } } } \right]}_{{g_{kk } } } \\     & + { \left ( { { \rm{tr}}\left ( { { { { \bf{\hat { \underline{c}}}}}_{sr , k } } } \right ) } \right)^ { - 2}}\sum\limits_{j = 1,j \ne k}^k { \nu _ { s , j}}{e_{s , j } } \\    & \times\underbrace { \mathbb e\left [ { { \bf{\hat { \underline{h}}}}_{sr , k}^h{\bf{{\dot{p}}}}_r^h{{\bf{h}}_{sr , j}}{\rm{diag}}\left ( { { { \bf{{\dot{p}}}}_{s , j}}{\bf{{\dot{p}}}}_{s , j}^h } \right){\bf{h}}_{sr , j}^h{{\bf{{\dot{p}}}}_r}{{{\bf{\hat { \underline{h}}}}}_{sr , k } } } \right]}_{{g_{kj } } } \end{aligned}\ ] ] according to the independence between @xmath556 and @xmath557 when @xmath558 , @xmath559 can be derived as @xmath560 moreover , substituting ( [ 4 - 7 ] ) and ( [ 4 - 8 ] ) in theorem [ t3 ] into ( [ d9 ] ) and using the large-@xmath28 approximation in ( [ x5 ] ) , it is straightforward but tedious to show that @xmath262 can be derived as @xmath561 where @xmath562 . inserting ( [ d10 ] ) and ( [ d12 ] ) into ( [ d9 ] ) , the expression of @xmath534 $ ] is obtained . _ ( b ) scaling behavior : _ using the similar approach as that for @xmath540 $ ] , it can be shown that , as @xmath25 is fixed , the first and second terms on the right - hand side of ( [ d9 ] ) scale as @xmath149 and @xmath273 , respectively . thus , @xmath534 = { \cal o}(1)$ ] since @xmath563 . as @xmath25 scales with @xmath28 , it is shown in @xcite that the eigenvectors of @xmath564 form a unitary dft matrix in the large @xmath25 regime . in this case , we have @xmath565 . applying this property and lemma [ lm2 ] on ( [ d12 ] ) , it can be verified that the first term on the right - hand side of ( [ d9 ] ) scales as @xmath566 as @xmath25 scales with @xmath28 . this completes the proof . h. a. suraweera , i. krikidis , g. zheng , c. yuen and p. j. smith , `` low - complexity end - to - end performance optimization in mimo full - duplex relay systems , '' _ ieee trans . wireless commun . _ , vol . 13 , pp . 913 - 927 , feb . 2014 .    c. y. a. shang , p. j. smith , g. k. woodward and h. a. suraweera , `` linear transceivers for full duplex mimo relays , '' in proc australian communications theory workshop ( ausctw 2014 ) , sydney , australia , feb . 2014 , pp . d. nguyen , l .- n . tran , p. pirinen , and m. latva - aho , `` precoding for full duplex multiuser mimo systems : spectral and energy efficiency maximization , '' _ ieee trans . signal process . 16 , pp . 4038 - 4050 , aug . 2013 . a. sabharwal , p. schniter , d. guo , d. w. bliss , s. rangarajan , and r. wichman , `` in - band full - duplex wireless : challenges and opportunities , '' _ ieee j. sel . areas commun . 1637 - 1652 , sep . 2014 . h. q. ngo , h. a. suraweera , m. matthaiou , and e. g. larsson , `` multipair full - duplex relaying with massive arrays and linear processing , '' _ ieee j. sel . areas commun . _ , 9 , pp . 1721 - 1737 , sep . 2014 . j. zhang , c. k. wen , s. jin , x. gao , k. wong , `` on capacity of large - scale mimo multiple access channels with distributed sets of correlated antennas '' _ ieee j. sel . areas commun . 133 - 148 , feb . 2013 . h. yin , d. gesbert , l. cottatellucci , `` dealing with interference in distributed large - scale mimo systems : a statistical approach , '' _ ieee j. sel . topics signal process . _ , vol . 8 , no . 942 - 953 , aug . 2014 . g. santella and f. mazzenga , `` a hybrid analytical - simulation procedure for performance evaluation in m - qam - ofdm schemes in presence of nonlinear distortions , '' _ ieee trans . 142 - 151 , feb . 1998 . h. suzuki , t. tran , i. b. collings , and et al , `` transmitter noise effect on the performance of a mimo - ofdm hardware implementation achieving improved coverage , '' _ ieee j. sel . areas commun . _ , vol . 26 , pp . 867 - 876 , aug . 2008 . j. w. silverstein and z. d. bai , `` on the empirical distribution of eigenvalues of a class of large dimensional random matrices , '' _ journal of multivariate analysis _ , 175 - 192 , 1995 . s. wagner , r. couillet , m. debbah , and d. t. m. slock , `` large system analysis of linear precoding in correlated miso broadcast channels under limited feedback , '' _ ieee trans . inf . theory _ 58 , no . 7 , pp . 4509 - 4537 , jul . 2012 . e. bj@xmath567rnson , p. zetterberg , m. bengtsson , and b. ottersten , `` capacity limits and multiplexing gains of mimo channels with transceiver impairments , '' _ ieee commun . 91 - 94 , 2013 . p. c. weeraddana , m. codreanu , m. latva - aho , and a. ephremides , `` resource allocation for cross - layer utility maximization in wireless networks , '' _ ieee trans . 2790 - 2809 , 2011 . e. bj@xmath567rnson , d. hammarwall , and b. ottersten , `` exploiting quantized channel norm feedback through conditional statistics in arbitrarily correlated mimo systems , '' _ ieee trans . signal process . 4027 - 4041 , 2009 . xiaochen xia received his b.e . degree in electronic science and technology from tianjin university ( tju ) in 2010 , and m.s . degree in communication and information system from pla university of science and technology ( plaust ) , in 2013 . he is currently working toward the ph.d . degree in institution of communications engineering , plaust . his research interests include relaying network , full - duplex communication , network coding , mimo techniques . he received the 2013 excellent master degree dissertation award of jiangsu province , china . dongmei zhang received the b.s . degree in institute of communications engineering , nanjing , china in 1994 , and the m.s . degree in communications and information system from the institute of communications engineering , nanjing , china in 1997 . she is currently an associate professor with the pla university of science and technology , nanjing , china . her research interests include broadband wireless communications , cognitive radio , wireless resource allocation . kui xu received the b.s . degree and ph.d . degree from the pla university of science and technology ( plaust ) , nanjing , china in 2004 and 2009 . he is currently a lecturer in the plaust . his research interests include broadband wireless communications , signal processing for communications and wireless communication networks . he is the author of about 50 papers in refereed journals and conference proceedings and holds 5 patents in china . xu currently served on the technical program committee of the ieee wcsp 2014 . he received the ursi young scientists award in 2014 and the 2010 ten excellent doctor degree dissertation award of plaust . wenfeng ma received the b.s . degree in microwave communication from the institute of communications engineering , nanjing , china in 1995 , and the m.s . degree in communications and information system from the institute of communications engineering , nanjing , china in 1998 , and the ph.d . degree in communications and information system from the pla university of science and technology , nanjing , china , in 2002 . he is currently an associate professor with the pla university of science and technology , nanjing , china . his research interests include wireless communication networks , broadband wireless communications . youyun xu was graduated from shanghai jiao tong university with ph.d . degree in information and communication engineering in 1999 . he is currently a professor with nanjing institute of communication engineering , china . he is also a part - time professor with the institute of wireless communication technology of shanghai jiao tong university ( sjtu ) , china . professor xu has more than 20 years professional experience of teaching and researching in communication theory and engineering . now , his research interests are focusing on new generation wireless mobile communication system ( imt - advanced and related ) , advanced channel coding and modulation techniques , multi - user information theory and radio resource management , wireless sensor networks , cognitive radio networks , etc .
Abstract:

Answer: this paper studies the massive mimo full - duplex relaying ( mm - fdr ) , where multiple source - destination pairs communicate simultaneously with the help of a common full - duplex relay equipped with very large antenna arrays . different from the traditional mm - fdr protocol , a general model where sources / destinations are allowed to equip with multiple antennas is considered . in contrast to the conventional mimo system , massive mimo must be built with low - cost components which are prone to hardware impairments . in this paper , the effect of hardware impairments is taken into consideration , and is modeled using transmit / receive distortion noises . we propose a low complexity hardware impairments aware transceiver scheme ( named as hia scheme ) to mitigate the distortion noises by exploiting the statistical knowledge of channels and antenna arrays at sources and destinations . a joint degree of freedom and power optimization algorithm is presented to further optimize the spectral efficiency of hia based mm - fdr . the results show that the hia scheme can mitigate the  ceiling effect " appears in traditional mm - fdr protocol , if the numbers of antennas at sources and destinations can scale with that at the relay . massive mimo full - duplex relaying , hardware impairments , transceiver design , joint degree of freedom and power optimization , achievable rate .