Exemplary embodiments of the invention will hereafter be described with reference to the accompanying drawings, wherein like numerals will denote like elements.
a is a contour plot of a virtual channel power matrix for a first channel configuration in accordance with an exemplary embodiment.
b illustrates the first channel configuration in accordance with an exemplary embodiment.
a is a contour plot of a virtual channel power matrix for a second channel configuration in accordance with an exemplary embodiment.
b illustrates the second channel configuration in accordance with an exemplary embodiment.
a is a contour plot of a virtual channel power matrix for a third channel configuration in accordance with an exemplary embodiment.
b illustrates the third channel configuration in accordance with an exemplary embodiment.
With reference to
Multiple antenna arrays may be used for transmitting data in wireless communication systems. For example, multiple antennas may be used at both the transmitter and at the receiver as shown with reference to the exemplary embodiment of
A virtual channel representation that provides an accurate and analytically tractable model for physical wireless channels is utilized where H denotes an N×N virtual channel matrix representing N antennas at the transmitter and the receiver. The virtual representation is analogous to representing the channel in beamspace or the wavenumber domain. Specifically, the virtual representation describes the channel with respect to spatial basis functions defined by virtual fixed angles that are determined by the spatial resolution of the arrays. With reference to
The dominant non-vanishing entries of the virtual channel matrix reveal the statistically independent degrees of freedom (DoF), D, in the channel, which also represent the number of resolvable paths in the scattering environment. For sparse channels, D<N2. With reference to
A family of channels is described by two parameters (p,q), D=pq that represent different configurations of the D<N2 DoF. For all feasible (p,q), the MIMO capacity of the corresponding channel configuration is accurately approximated by
C(N,ρ,D,p)≈p log(1+ρD/p2) (1)
where ρ denotes the transmit SNR (can be interpreted as the nominal received SNR if an attenuation factor is included to reflect path loss relating the total power at the receiver to the total transmitted power), p represents the multiplexing gain (MG) or the number of parallel channels (number of independent data streams transmitted at the transmitting communication device), q represents the DoF per parallel channel, and ρD/p2=ρrx denotes the received SNR per parallel channel. From equation (1), increasing p comes at the cost of ρrx and vice versa. Based on an analysis of equation (1), on one extreme, beamforming channels (BF) in which the channel power is distributed to maximize ρrx at the expense of p result, and, on the other extreme, multiplexing channels (MUX) which favor p over ρrx result. The ideal channel (IDEAL) lies in between and corresponds to an optimal distribution of channel power to balance p and ρrx. p reflects the number of independent data streams, and hence the rate of transmission, whereas ρrx reflects the received SNR, and hence the reliability of decoding a particular data stream at the receiver. Maximizing capacity (maximum number of data streams that can be reliably communicated) involves optimally balancing p as a function of the operating SNR. The BF, MUX, and IDEAL configurations reflect the capacity-maximizing configurations at low, high, and medium SNRs, respectively. Precise values of low, high, and medium SNRs can be determined through measureed channel parameters, such as the number of dominant non-vanishing virtual channel entries and the total average power contributed by the spatial multi-path channel (the sum of the average powers of the dominant non-vanishing virtual channel entries).
With reference to
Three canonical antenna array configurations are sufficient for near-optimum performance over the entire SNR range as illustrated with reference to
C
α
=N
α log(1+ρNγ−2α) (2)
corresponding to D(N)=Nγ, p(N)=Nα, and q(N)=Nγ−α. Cα is plotted for 10 equally spaced values of αε[0,1] for γ=1 and N=25. With reference to
In a single-user MIMO system with a uniform linear array of Nt transmit and Nr receive antennas. The transmitted signal s and the received signal x are related by x=Hs+n where H is the MIMO channel matrix and n is the additive white Gaussian noise (AWGN) at the receiver. A physical multi-path channel can be accurately modeled as
where the transmitter and receiver arrays are coupled through L propagation paths with complex path gains {βl}, angles of departure (AoD) {θt,l} and angles of arrival (AoA) {θr,l}. In equation (3), αr(θr) and αt(θt) denote the receiver response and transmitter steering vectors for receiving/transmitting in the normalized direction θr/θt, where θ is related to the physical angle (in the plane of the arrays) φε[−π/2,π/2] as θ=d sin(φ)/λ, d is the antenna spacing and λ is the wavelength of propagation. Both αr(θr) and αt(θt) are periodic in θ with period one.
The virtual MIMO channel representation characterizes a physical channel via coupling between spatial beams in fixed virtual transmit and receive directions
where
are fixed virtual receive and transmit angles that uniformly sample the unit θ period and result in unitary discrete fourier transform matrices At and Ar. Thus, H and Hv are unitarily equivalent: Hv=ArHHAt. The virtual representation is linear and is characterized by the matrix Hv.
Virtual path partitioning relates the virtual coefficients to the physical paths gains
where Sr,m and St,n are the spatial resolution bins of size 1/Nr and 1/Nt corresponding to the m-th receive and n-th transmit virtual angle. Thus, Hv(m,n) is approximately the sum of the gains of all paths whose transmit and receive angles lie within the (m, n)-th resolution bin. If there are no paths in a particular resolution bin, the corresponding Hv(m,n)≈0. Each Hv(m,n) is associated with a disjoint set of physical paths and is approximately equal to the sum of the gains of the corresponding paths. It follows that the virtual channel coefficients are approximately independent. The virtual channel coefficients can be assumed to be statistically independent zero-mean Gaussian random variables in a Rayleigh fading environment. For a Rician environment (with a line-of-sight path or non-random reflecting paths), the virtual channel coefficients corresponding to line-of-sight (reflecting) paths can be modeled with an appropriate non-zero mean.
In Rayleigh fading, the statistics of H are characterized by the virtual channel power matrix Ψ:Ψ(m,n)=E└|Hv(m,n)|2┘. The matrices Ar and At constitute the matrices of eigenvectors for the transmit and receive covariance matrices, respectively: E[HHH]=AtΛtAtH and E[HHH]=ArΛrArH, where Λt=E[HvHHv] and Λr=E[HvHvH] are the diagonal matrices of transmit and receive eigenvalues (correlation matrices in the virtual domain). Ψ is the joint distribution of channel power as a function of the transmit and receive virtual angles. Λt and Λr are the corresponding marginal distributions:
An N×N Hv is sparse if it contains D<N2 non-vanishing coefficients. Each non-vanishing coefficient reflects the power contributed by the unresolvable paths associated with it. D reflects the statistically independent DoF in the channel and the channel power
In general, the sparser the Hv in the virtual domain, the higher the correlation in the antenna domain H. A sparse Hv can be modeled as
H
v
=M•H
iid (6)
where • denotes an element-wise product, Hiid is an iid matrix with CN(0,1) entries, and M is a mask matrix with D unit entries and zeros elsewhere. Under these assumptions, Ψ=M and the entries of Λr and Λt represent the number of non-zero elements in the rows and columns of M, respectively.
The ergodic capacity of a MIMO channel, assuming knowledge of H at the receiver, is given by
where ρ is the transmit SNR, and Q=E[ssH] is the transmit covariance matrix. The capacity-maximizing Qopt is diagonal. Furthermore, for general correlated channels, Qopt is full-rank at high SNR's, whereas it is rank-1 at low SNR's. As ρ is increased from low to high SNR's, the rank of Qopt increases from 1 to N.
The capacity of a sparse virtual channel matrix Hv depends on three fundamental quantities: 1) the transmit SNR ρ, 2) the number of DoF, D<N2, and 3) the distribution of the D DoF in the available N2 dimensions. For any ρ, there is an optimal configuration of the DoF characterized by an optimal mask matrix Mopt that yields the highest capacity at that ρ. The corresponding MIMO channel can be termed the IDEAL MIMO channel, and the resulting capacity can be termed the ideal MIMO capacity at that ρ.
Consider a fixed N and D<N2 and let M(D) denote the set of all N×N mask matrices with D non-zero (unit) entries. For any ρ, the ideal MIMO capacity is defined as
and an Mopt that achieves Cid(N,D,ρ) defines the IDEAL MIMO Channel at that ρ.
Mopt is not unique in general. The family of mask matrices is defined by two parameters (p,q) such that D=pq. For D=Nγ,γε[0,2], the matrices can be further parameterized via p=Nα,αε[αmin,αmax] where αmin=max(γ−1,0) and αmax=min(γ,1), and q=D/p.
For a given D=Nγ,γε[0,2], and any p=Nα,αε[αmin,αmax] the mask matrix M(D,p) is an N×N matrix, but its non-zero entries are contained in a non-zero sub-matrix of size r×p,r=max(q,p), consisting of P non-zero columns, and q non-zero (unit) entries in each column. The corresponding r×p virtual sub-matrices {tilde over (H)}v defined by equation (6) satisfy ρc=D and their transmit and receive correlation matrices are given by
Since each {tilde over (H)}v defines a regular channel, the capacity maximizing input allocates uniform power over the non-vanishing transmit dimensions,
and no power in the remaining dimensions. The channel capacity for any M(D,p) is characterized by equation (1) which was derived for large N, but yields accurate estimates even for relatively small N. For sufficiently large N, the capacity of the MIMO channel defined by the mask M(D,p) is accurately approximated as a function of ρ by
For a given ρ, the IDEAL MIMO Channel is characterized by M(D,popt)popt where
Different values of p reveal a multiplexing gain (MG) versus received SNR tradeoff. In equation (11),
is the received SNR per parallel channel. Thus, increasing the MG comes at the cost of a reduction in ρrx and vice versa. For ρ<ρlow, the optimal BF configuration (
The ratio ρhigh/ρlow=(pmax/pmin)2 attains its largest value, N2, for γ=1(D=N), whereas it achieves its minimum value of unity for γ=0(D=1) or γ=2(D=N2). Thus, the MG-ρrx tradeoff does not exist for the extreme cases of highly correlated (γ=0) and iid (γ=2) channels. On the other hand, the impact of the MG-ρrx tradeoff on capacity is maximum for γ=1 corresponding to D=N.
An antenna spacing at the transmitter is denoted dt and at the receiver is denoted dr. Consider D=Nγ, γε[1,2) (since for γε(0,1), it is advantageous to use fewer antennas to effectively increase γ to 1). For a given array dimension N, a class H(D) of channels is said to be randomly sparse with D DoF if it contains L=D<N2 resolvable paths that are randomly distributed over the maximum angular spreads or some sufficiently large antenna spacings dt,max and dr,max; that is, (θr,l,θt,l)ε[−½,½]×[−½,½] in equation (3).
The maximum antenna spacings correspond to the choice p=pmax=N (MUX configuration); that is, (dt,max,dr,max)pmax. For any p,pmin≦p≦pmax define the antennas spacings
where r=max(q,p) and q=D/p. As a result, for each p, the non-vanishing entries of the resulting Hv are contained within an r×p sub-matrix {tilde over (H)}v with power matrix
By way of a proof, for a given scattering environment, the channel power does not change with antenna spacing. By assumption we have ρc=tr(E[HvHvH])=D. Also by assumption, the D randomly distributed paths cover maximum angular spreads (AS's) at the maximum spacings. Since θ=d sin(φ)/λ, where φ is the physical angle associated with a path (which remains unchanged), the dr and dt in (11) result in smaller AS's: {−p/2N,p/2N] at the transmitter and [−r/2N,r/2N] at the receiver. Since the spacing between virtual angles is Δθ=1/N, it follows that only p=p/N/Δθ virtual angles lie within the reduced AS at the transmitter and only r virtual angles lie within the reduced angular spread at the receiver. Thus, the non-zero entries in Hv are contained in a sub-matrix {tilde over (H)}v of size r×p. The channel power ρc=D is uniformly distributed over its entries so that
where the expectation is over the statistics of the D non-vanishing coefficients as well as their random locations. The power matrix of the reconfigured channel corresponding to the spacings in (13) satisfies: Ψ=M(D,p) for p≦√{square root over (D)}(q≧p), but Ψ≠M(D,p) for p>√{square root over (D)}(q<p).
In randomly sparse physical channels, the virtual channel matrix generated by reconfiguring antenna spacings has identical statistics (marginal and joint) to those generated by the mask matrix M(D,p) for P≦q, but only the marginal statistics are matched for p>q. It follows that the reconfigured channel achieves the capacity corresponding to M(D,p) for P≦q, but the capacity may deviate a little for p>q especially at high SNR's since the reconfigured channel always has a kronecker (separable) structure whereas M(D,p) is non-separable for p>q. With this qualification, in randomly sparse physical channels, the (capacity maximizing) IDEAL MIMO channel at any transmit SNR can be created by choosing dr,opt and dt,opt in (13) corresponding to popt defined in (12).
Three channel configurations are highlighted in pbf=pmin=1, the transmitter array is in a low-resolution configuration (
pid=√{square root over (D)}=√{square root over (N)}, both the transmitter and receiver arrays are in a medium-resolution configuration (
pmux=pmax=N, both the transmitter and receiver arrays are in a high-resolution configuration (
With reference to
The effect of decreasing dt with ρ is to concentrate channel power in fewer non-vanishing transmit dimensions. As a result, the number of non-vanishing transmit eigenvalues is reduced, but their size is increased. This reflects a form of source-channel matching: the rank of the optimal input is better-matched to the rank of Hv. As a result, less channel power (none for regular channels) is wasted as the multiplexing gain is optimally reduced through dtp. Physically, as dt is decreased, fewer data streams (p) are transmitted over a corresponding number of spatial beams, whereas the width of the beams gets wider (see
With reference to
T/R signal processor 102 forms the transmitted signals s(t) transmitted from each antenna of the plurality of antennas 22. 22 in the transmitting device. In a receiving device, the processor 102 determines the way in which the signals received on the plurality of antennas 22 are processed to decode the transmitted signals from the transmitting device, for example, based on the modulation and encoding used at the transmitting device. Actuator 104 adjusts an antenna spacing of the plurality of antennas 22 based on the determined optimum antenna spacing. For example, actuator 104 adjusts the antenna spacing based on equation (13). The actuators may be based on any available or emerging technology for reconfiguring the array configuration (antenna spacing for uniform linear arrays), such as MEMS technology. For arrays other than uniform linear arrays, the array reconfiguration may be based on other mechanisms related to the possible radiation patterns (e.g., appropriate excitations in a fractal array).
Memory 106 stores antenna spacing application 110, in addition to other information. Device 100 may have one or more memories 106 that uses the same or a different memory technology. Memory technologies include, but are not limited to, random access memory, read only memory, flash memory, etc. In an alternative embodiment, memory 106 may be implemented at a different device.
Processor 108 executes instructions that may be written using one or more programming language, scripting language, assembly language, etc. The instructions may be carried out by a special purpose computer, logic circuits, or hardware circuits. Thus, processor 108 may be implemented in hardware, firmware, software, or any combination of these methods. The term “execution” is the process of running an application or the carrying out of the operation called for by an instruction. Processor 108 executes antenna spacing application 110 and/or other instructions. Device 100 may have one or more processors 108 that use the same or a different processing technology. In an alternative embodiment, processor 108 may be implemented at a different device.
Antenna spacing application 110 is an organized set of instructions that, when executed, cause device 100 to determine an antenna spacing. Antenna spacing application 110 may be written using one or more programming language, assembly language, scripting language, etc. In an alternative embodiment, antenna spacing application 110 may be executed and/or stored at a different device.
Determining the capacity-optimal channel configuration may include use of channel sounding. Two channel parameters can be determined through channel sounding: 1) the total received signal power as a function of the total transmitted signal power to determine the operating SNR (this accounts for the path loss encountered during propagation and the total power contributed by the multiple paths in the scattering environment), and 2) the number of dominant non-vanishing entries in the virtual channel matrix. Knowledge of 2) can lead to the determination of 1). With reference to 2), a variety of channel sounding/estimation methods may be used. For example, in the method proposed in Kotecha and Sayeed, “Transmit Signal Design for Optimal Estimation of Correlated MIMO Channlels,” IEEE Transactions on Signal Processing, February 2002, training signals are transmitted sequentially on different virtual transmit beams at the first (transmitting device) and the entries in the corresponding column of the virtual channel matrix Hv are estimated by processing the signals in the different virtual beam directions at the second (receiving) device. In this fashion, channel coefficients in different columns of the virtual channel matrix are sequentially estimated at the receiving device from the sequential transmissions in different virtual directions from the transmitting device.
By performing channel sounding (estimation of the virtual channel matrix entries) a sufficient number of times, the average power in different virtual channel coefficients can be estimated to form an estimate of the virtual channel power matrix Ψ. Once the virtual channel power matrix Ψ is estimated, the effective operating SNR can be directly estimated from the total channel power (sums of all the entries in the power matrix) and includes the impact of path loss by comparing the total transmitted power to the total received power. From the virtual channel power matrix Ψ, the dominant number of entries in the power matrix can be estimated by comparing to an appropriately chosen threshold (to discount virtual channel coefficients with insignificant power) yielding the number of degrees of freedom D in the channel.
Based on knowledge of D, the optimal array configurations can be determined at any desired operating SNR via equations (12) and (13). The maximum antenna spacings in equation (13) defining the reference MUX configuration are determined from the physical angular spread of the scattering environment (the spacings are adjusted so that the physical channel exhibits maximum angular spread in the virtual (beamspace) domain). The estimation of the channel power matrix is performed at the receiving device, and the value of D is transmitted back to the transmitting device so that the optimum transmit array configuration can be chosen for a given operating SNR. The receiving device also configures its array configuration according to D and the operating SNR. The procedure discussed above is implicitly based on a scattering environment in which the paths are randomly and uniformly distributed over the angular spreads. Appropriate modifications may be made for non-uniform distribution of scattering paths by those knowledgeable in the art for further enhancements in performance.
The foregoing description of exemplary embodiments of the invention have been presented for purposes of illustration and of description. It is not intended to be exhaustive or to limit the invention to the precise form disclosed, and modifications and variations are possible in light of the above teachings or may be acquired from practice of the invention. The embodiments were chosen and described in order to explain the principles of the invention and as practical applications of the invention to enable one skilled in the art to utilize the invention in various embodiments and with various modifications as suited to the particular use contemplated. It is intended that the scope of the invention be defined by the claims appended hereto and their equivalents.
This invention was made with United States government support awarded by the following agencies: NSF 0431088. The United States government has certain rights in this invention.