new frontiers in feedback for interference...

45
The Wireless Networking and Communications Group New Frontiers in Feedback for Interference Alignment Robert W. Heath Jr. & Omar El Ayach The University of Texas at Austin http://www.profheath.org November 2011 Also associated with Kuma Signals LLC and MIMO Wireless Inc. Performed by ONR grant N000141010337 Friday, November 18, 2011

Upload: others

Post on 07-Aug-2020

2 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

The Wireless Networking and Communications Group

New Frontiers in Feedback for Interference Alignment

Robert W. Heath Jr. & Omar El Ayach

The University of Texas at Austinhttp://www.profheath.org

November 2011

Also associated with Kuma Signals LLC and MIMO Wireless Inc.

Performed by ONR grant N000141010337

Friday, November 18, 2011

Page 2: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

WNCG Enhances UT Visibility Consistently win large multi-PI grants, e.g.

$6.5m DARPA IT-MANET (PI: Andrews) $1m Intel / Cisco grand challenge on video networks (PI: Heath) $1.3m on cross-layer delay-tolerant nets (PI: Shakkottai)

Seven major best paper awards in last four yearsWNCG Impacts Industry Four definitive textbooks on wireless

Widely-cited magazine articles on hot topics

Developed key features of wireless standards

Software packages and toolkits

Host popular annual conference. www.twsummit.com

WNCG:  The  USA’s  Premier  Wireless  Research  Center

1

WNCG Affiliates

WN

CG

Fac

ulty

16 faculty over 3 departments, all actively involved in center activities

10 NSF Career Awards Comprehensive wireless expertise $4m/year external funding

WN

CG

Stu

dent

s

120 PhD students in pooled space Many co-advised students 64%/yr intern at affiliates in last 4 yrs

• Students receive perks, special awards and travel funds for help with affiliates• Staff, space & other resources shared efficiently amongst all faculty/students

Affiliates champion large federal proposals, provide technical input/feedback, unrestricted gift funds

WNCG provides pre-prints, pre-competitive research ideas, vast expertise, first access to students

Wireless  Networking  andCommunica?ons  Group

2

Friday, November 18, 2011

Page 3: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

OutlineIntroduction to interference alignment

Analog feedback for interference alignment

Limited feedback for interference alignment

Conclusions

3

Friday, November 18, 2011

Page 4: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

Wireless Systems

4

Interference limits performance

Cellular MANET LAN

Friday, November 18, 2011

Page 5: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

K-User Interference Channel

5

transmitter receiverdirect channel

interferencereceived

interferencecreated

Friday, November 18, 2011

Page 6: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Conventional Approach

6

transmitter receiver

Multiple access protocol enables sharing

Friday, November 18, 2011

Page 7: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Miracle of Interference Alignment

7

transmitter receiver

Space occupied by interferers is reduced

interferencealigned

Friday, November 18, 2011

Page 8: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Interference Alignment: Subspaces

8

H!!

H""

H##

H$$

Alignment in space (MIMO) or frequency (OFDM)Friday, November 18, 2011

Page 9: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

Hot Topics

9

Theory Practice

Random access protocols

Proof-of-concept in real channelsAlgorithms

Prototyping

Impact of estimation error

Practical feedback strategies

Degrees of freedom

Blind alignment

Capacity scaling

Relays

Feasibility

NeutralizationImpact of spatial correlation

Realistic assessment of overhead

Realistic assessment of overheadSuccessive alignment

Interferencepricing

IA tutorial http://newport.eecs.uci.edu/~syed/papers/fntfinaltutorial.pdfIA bibliography http://www.profheath.org/research/interference-alignment/

System level performance

Partial connectivity

Alignment in cellular systems

Friday, November 18, 2011

Page 10: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Practical Challenge I:Real Propagation Channels

10

Antenna Correlation

Line of sightCorrelation across users

Good BadSome randomness

IA works!Indoor & Outdoor Testbed

3 MIMO users

IA designed for

random channels

IA will be used in

“real” channels

but

with a small performance

gap0 5 10 15 20 25 30 35 40

0

5

10

15

20

25

30

35

40

SNR (dB)

Su

m R

ate

b/s

/Hz

Interference Alignment

Time Division

Older Strategies

Real Channels

O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility of Interference Alignment of Measured MIMO-OFDM Channels'', IEEE Trans. on Veh. Tech.. vol. 59. No. 9. pp. 4309-4321, Nov. 2010.

Friday, November 18, 2011

Page 11: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Practical Challenge II:Channel Uncertainty

11

IA designed for

ideal channels

Only have

estimated channels

butDelay

Pilot contamination

Estimation error

Training Sequences

Good Bad

0 10 20 30 400

5

10

15

20

25

30

35

40

45

50

SNR (dB)

Sum

Rat

e

6 Pairs

3 Pairs

IncreasingEstimation

Error

75%

Red

uctio

n

70%

Redu

ctio

n

! = 0.02

! = 0.2

Real Systems

estimate =channel + error

channel estimate creates self-interference

0 20 40 60 80 100 120 140 160 180 200

0

2

4

6

8

time

magnitude

Errorfloorsevere for

more users

B. Nosrat-Makouei, J. G. Andrews, and  R. W. Heath, Jr., ``MIMO Interference Alignment Over Correlated Channels with Imperfect CSI," IEEE Trans. on Signal Processing, vol. 59, no. 6, pp. 2783-2794, June 2011

Friday, November 18, 2011

Page 12: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Practical Challenge III:Overhead

12

Approach

IA designed for

no overhead overhead cost

but

TX

TX

RX

RX

RX

TX

Partition the network

FRAME k FRAME k+1

Training + Feedback

Training + Feedback

Data (all nodes)

L(M,K) symbols

T - L(M,K) symbols

Account for overheads Channel training Channel feedback

Insights

Rates 0 with many users

1 2 3 4 5 60

5

10

15

20

25

30

35

P (# of groups)

Sum

Rat

e b/

s/H

z = 0.9

= 0.7

= 0.5

= 0.3

= 0.1

Smaller partitions usually optimum

will have to pay

# partitions

sum

rat

e

S. W. Peters and R. W. Heath, Jr., ``User Partitioning for Less Overhead in MIMO Interference Channels'', submitted to IEEE Trans. on Wireless, see arxiv.

Friday, November 18, 2011

Page 13: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Practical Challenge IV:Channel State Information Feedback

13

Channel estim

ation

Cal

cula

te P

reco

ders

CSI feedback

Send payload data

1

2

3

4

Modulation

Pilot Symbols

Channel Estimation

SimpleProcessing

Pilot Symbols

Channel Estimation

Channel Quantization Channel

Index

Modulation0110100

Limited feedback

Limited feedbackBetter compression through structure

Analog feedbackGood for fast & simple CSI transfer

Analog feedback

Performed by ONR grant N000141010337

CSI is the hidden IA “killer”

Main Challenges

Feedback overhead scales with K2

Feedback quality must increase at high SNR

Friday, November 18, 2011

Page 14: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Related Work on CSI FeedbackExploit reciprocity [GomCadJaf08, ShiBerHon09]

Limited feedback (CSI quantization)

Grassmannian limited feedback for SISO IA over frequency [ThuBol09]

Grassmannian limited feedback for MIMO IA over frequency [KriVar10]

Exploiting network topology

User grouping algorithms [PetHea10, MaLiLiuChe11]

Partially connected channels could affect feedback [ChoChu08, GuiGes11]

Blind IA transmission strategies [Jaf10, GouJaf11, WanPapRamCai11]

14

Contribution: New approaches for analog and limited feedback

Friday, November 18, 2011

Page 15: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

OutlineIntroduction to interference alignment

Analog feedback for interference alignment

Limited feedback for interference alignment

Conclusions

15

Friday, November 18, 2011

Page 16: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

MIMO System Model

Interference alignment with MIMO signal model

For submission to IEEE Trans. Wireless Commun. 5

Under our assumptions, the received signal at sink node i can be written as

yi =

�P

diHi,iFisi +

� �=i

�P

d�Hi,�F�s� + vi,

where yi is the Nr × 1 received signal vector, Hi,� is the Nr ×Nt channel matrix from source �

to sink i, Fi is the Nt× di unitary precoding matrix used at source i, si is the di× 1 transmitted

symbol vector at node i, with unit norm elements, i.e. E [�si�2] = di, and vi is a complex vector

of i.i.d circularly symmetric white Gaussian noise with covariance matrix σ2INr .

We place no assumption on the reciprocity of the forward and reverse channel. This is similar

to a frequency division duplexed (FDD) system model, for example, in which the forward and

reverse channels are uncorrelated. For the reverse or feedback channel we write the received

signal at source node i as

←−y i =

�Pf

NtGi,i

←−x i +

��=i

�Pf

NtGi,�

←−x � + νi, (1)

where Pf is the transmit power used to transmit pilot and feedback symbols, Gi,� is the Nt×Nr

reverse channel between sink node i and source node � with i.i.d CN (0, 1) elements1,←−x i is

the symbol vector with unit norm elements sent by sink i, and νi is a complex vector of i.i.d

circularly symmetric white Gaussian noise with covariance matrix σ2INt .

B. Interference Alignment

Interference alignment for the MIMO interference channel is a linear precoding technique

that can achieve the maximum multiplexing gain, or degrees of freedom defined as limP→∞

Rsumlog2 P

,

which in this case isKNr2 , when coding over infinitely many channel extensions [10]. While

the maximum multiplexing gain may not be achieved without time extensions, interference

alignment for the constant MIMO channel can still provide an increase in sum rate [21]. To

do so, given global channel knowledge, interference alignment computes the transmit precoders

Fi to align interference at all receivers in a strict subspace of the received signal space, thus

leaving interference free dimensions for the desired signal. While interference alignment is only

one of the many precoding strategies for the interference channel [3], [22], [23], some of which

1Note the effective reversal in the indexing of the channel, while the indexing still has the form “sink,source” the transmit

receive roles have been switched.

Desired Signal Interference

IA precoders push interference to a low-dimensional subspace

16

Resulting sum rate with zero-forcing receiver

For submission to IEEE Trans. Wireless Commun. 7

signals and treating interference as noise, is

Rsum =K�

i=1

di�

m=1

Rmi =

K�

i=1

di�

m=1

log2

�1 +

Pdi|(wm

i )∗Hi,if

mi |2

I1i,m + I2

i,m + σ2

�, (5)

where I1i,m and I2

i,m are the inter-stream and inter-user interference, respectively. These sum

interference terms are given by

I1i,m =

��=m

P

di

��(wmi )

∗Hi,if

�i

��2 ,

I2i,m =

k �=i

dk�

�=1

P

dk

��(wmi )

∗Hi,kf

�k

��2 .

In the presence of perfect channel knowledge, and for an achievable degree of freedom vector

d = [d1, d2, . . . , dK ], equations (2), (3), and (4) can be satisfied with probability one and thus

I1i = I2

i = 0. This gives

limP→∞

Rsum

log2 P= lim

P→∞

�i,m

log2

�1 +

Pdi|(wm

i )∗Hi,ifmi |2

σ2

log2 P

=K�

i=1

di ≤KNr

2.

It is not immediately clear, however, if the same sum rate scaling behavior can be expected

from a network with only imperfect knowledge of the channel derived from noisy feedback.

Results on single user MIMO prove an acceptable constant loss in sum rate due to imperfect

channel state information [25]. In multi-user scenarios, however, the cost of imperfect channel

knowledge may be much higher, potentially resulting in the loss of the channel’s multiplexing

gain [26] which saturates achieved sum rate at high SNR [27]. In Section III, we show that such

performance can be expected from a realistic system via interference alignment provided that the

quality of channel knowledge scales sufficiently with transmit power, or effectively the forward

channel’s SNR. This is similar to the results presented in [13] and [14] for SISO channels where

the feedback scaling is in terms of codebook size. We discuss several advantages of using the

feedback scheme proposed in this paper over limited feedback quantization.

=Leakage interference depends on CSI quality

Friday, November 18, 2011

Page 17: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Problem Statement

17

Design low overhead feedback strategy for MIMO IA

Prior work on analog feedback:Analog feedback for MISO broadcast [MarHoc06, SamMan06]

Comparison with quantized feedback [CaiJinKobRav07, KobJibCai10]

Solution: Analog channel state feedback

Objectives

Fast MIMO CSI transfer

Achieves multiplexing gains

Fixed complexity (independent of SNR)

Existing problems

CSI is quality and delay sensitive

Exploding codebooks in

limited feedback

Prior work cannot handle MIMO

Receivers estimate MIMO channel

Transmitters calculate the transmit directions

Receivers feedback the quantized channels

Transmitters send payload data

Friday, November 18, 2011

Page 18: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Limited feedback

Analog vs. Limited Feedback

18

Pilot Symbols

Channel Estimation

Channel Quantization Channel

Index

Codebook

Modulation0110100Hk,l

ModulationPilot

SymbolsChannel

EstimationSimple

ProcessingHk,l Hk,l

Analog feedback

Friday, November 18, 2011

Page 19: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Assumptions

Any number of antennas

Any number of streams

FDD (no reciprocity)

Block fading model

SNR is known

Analog feedback IA

Transmit power (P)

Feedback power (Pf)

Variable training and feedback

Analog Feedback

2. Train reverse channel (Gik)

3. Feedback Hik as QAM symbols

4. Use CSI estimates to do IA

For submission to IEEE Trans. Wireless Commun. 8

III. INTERFERENCE ALIGNMENT WITH ANALOG FEEDBACK

In this section we propose a feedback and transmission strategy based on analog feedback and

interference alignment, which uses the estimated channels as if they were the true propagation

channels.

A. Analog Feedback

To feedback the forward channel matrices Hi,� reliably across the feedback channel given in

(1), we propose dividing the feedback stage into two phases. First, each source learns its reverse

channels. Second, the forward channels are fed back and estimated. We neglect the initial training

phase in which sinks learn the forward channels and, thus, assume they have been estimated

perfectly. This is since imperfect CSI at the receiver will adversely affect performance of all

feedback schemes and is not exclusive to analog feedback. In fact, any error in estimating the

forward channels will only add an error term to the forward channels in (7), which also decays

with transmit power, and thus similar results can be shown.

1) Reverse Link Training: To learn all reverse links, the K sink nodes must transmit pilot

symbols over a period τp ≥ KNr. Similar to the analysis done in [15], we let the sinks collectively

transmit a KNr×τp unitary matrix of pilots Φ shown to be optimal in [28]. Such reverse channel

training requires no assumptions other than training sequences be known to all sink nodes, to

guarantee orthogonality. Of course, if we enforce orthogonality in time, this assumption reduces

to synchronization, and exact sequences need not be known to all sinks.

Let←−Y i =

�←−y i[1] . . . ←−

y i[τp]�

be the Nt × τp received training matrix at source node i, and

let←−Yp =

�←−Y

∗1,←−Y

∗2, . . . ,

←−Y

∗K

�∗be the composite received training matrix. We write the received

training as←−Yp =

�τpPf

NrGΦ+V,

where G is the KNt × KNr composite reverse channel matrix between all sinks and sources

and V is a KNt× τp matrix with i.i.d CN (0, σ2INr) elements. Using the received training, each

source derives an MMSE estimate of its channels all of which can be together written as

�G =

�τpPf

Nr

σ2 + τpPf

Nr

←−YpΦ

∗. (6)

For submission to IEEE Trans. Wireless Commun. 8

III. INTERFERENCE ALIGNMENT WITH ANALOG FEEDBACK

In this section we propose a feedback and transmission strategy based on analog feedback and

interference alignment, which uses the estimated channels as if they were the true propagation

channels.

A. Analog Feedback

To feedback the forward channel matrices Hi,� reliably across the feedback channel given in

(1), we propose dividing the feedback stage into two phases. First, each source learns its reverse

channels. Second, the forward channels are fed back and estimated. We neglect the initial training

phase in which sinks learn the forward channels and, thus, assume they have been estimated

perfectly. This is since imperfect CSI at the receiver will adversely affect performance of all

feedback schemes and is not exclusive to analog feedback. In fact, any error in estimating the

forward channels will only add an error term to the forward channels in (7), which also decays

with transmit power, and thus similar results can be shown.

1) Reverse Link Training: To learn all reverse links, the K sink nodes must transmit pilot

symbols over a period τp ≥ KNr. Similar to the analysis done in [15], we let the sinks collectively

transmit a KNr×τp unitary matrix of pilots Φ shown to be optimal in [28]. Such reverse channel

training requires no assumptions other than training sequences be known to all sink nodes, to

guarantee orthogonality. Of course, if we enforce orthogonality in time, this assumption reduces

to synchronization, and exact sequences need not be known to all sinks.

Let←−Y i =

�←−y i[1] . . . ←−

y i[τp]�

be the Nt × τp received training matrix at source node i, and

let←−Yp =

�←−Y

∗1,←−Y

∗2, . . . ,

←−Y

∗K

�∗be the composite received training matrix. We write the received

training as←−Yp =

�τpPf

NrGΦ+V,

where G is the KNt × KNr composite reverse channel matrix between all sinks and sources

and V is a KNt× τp matrix with i.i.d CN (0, σ2INr) elements. Using the received training, each

source derives an MMSE estimate of its channels all of which can be together written as

�G =

�τpPf

Nr

σ2 + τpPf

Nr

←−YpΦ

∗. (6)

For submission to IEEE Trans. Wireless Commun. 6

marginally outperform it at low SNR [3], it is analytically tractable. Its complete interference

suppression properties make it especially amenable to the study of performance with feedback

and imperfect CSI.

To express the conditions for alignment, consider the K-user interference channel with pre-

coding presented in Section II-A and any corresponding achievable degree of freedom allocation

vector d = [d1 d2 . . . dK ]. Source node i sends its di spatial streams along the columns f�i of

the precoder Fi, resulting in a transmitted symbol

xi =1√diFisi =

1√di

di�

�=1

f�i s

�i i = 1, . . . , K

where we note that �f �i �2 = 1 and��s�i

��2 = 1 to satisfy the total power constraint with equality. We

assume equal power allocation since the gain observed from water-filling is at most a constant

and thus is negligible at high SNR [24].

While in general interference alignment can be used with any receiver design, the discussion

and proofs in this paper assume a linear zero-forcing receiver in which sink node i projects its

received signal on to the columns, w�i , of the Nr × di combiner Wi. Simulations in Section V

indicate that the same performance can be expected from an optimal receiver.

Writing the per stream input-output relation after projection gives

(wmi )

∗yi = (wm

i )∗√PHi,if

mi smi +

� �=m

(wmi )

∗√PHi,if

�i s

�i+

k �=i

dk�

�=1

(wmi )

∗√PHi,kf

�ks

�k+(wm

i )∗vi,

for m ∈ {1, . . . , di} and i ∈ {1, . . . , K}, where �wmi �2 = 1. At the output of these linear

receivers w�i , the conditions for perfect interference alignment can be restated as

(wmi )

∗Hi,if

�i = 0, ∀i, � �= m (2)

(wmi )

∗Hi,kf

�k = 0, ∀k �= i, and ∀m, � (3)

|(wmi )

∗Hi,if

mi | ≥ c > 0, ∀i,m (4)

where interference alignment is guaranteed by the first two conditions, and the third ensures the

decodability of the di desired streams.

The suboptimal sum rate achieved by the linear zero-forcing receiver, assuming Gaussian input

1. Train forward channels (Hik)

TX k-1

TX k

TX k+1

RX k-1

RX k

RX k+1

19

Friday, November 18, 2011

Page 20: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Analog Feedback

Done anyway for coherent detection Not specific to analog feedback

Let’s assume it’s done perfectly

Hk-1,k-1

Hk,k-1

TX k-1

TX k

TX k+1

RX k-1

RX k

RX k+1

20

1. Train forward channels (Hik)

Though we really don’t need to

Friday, November 18, 2011

Page 21: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Analog FeedbackGk-1,k-1

Gk,k-1

Feedback Power

Training Time

Training Sequence

For submission to IEEE Trans. Wireless Commun. 8

III. INTERFERENCE ALIGNMENT WITH ANALOG FEEDBACK

In this section we propose a feedback and transmission strategy based on analog feedback and

interference alignment, which uses the estimated channels as if they were the true propagation

channels.

A. Analog Feedback

To feedback the forward channel matrices Hi,� reliably across the feedback channel given in

(1), we propose dividing the feedback stage into two phases. First, each source learns its reverse

channels. Second, the forward channels are fed back and estimated. We neglect the initial training

phase in which sinks learn the forward channels and, thus, assume they have been estimated

perfectly. This is since imperfect CSI at the receiver will adversely affect performance of all

feedback schemes and is not exclusive to analog feedback. In fact, any error in estimating the

forward channels will only add an error term to the forward channels in (7), which also decays

with transmit power, and thus similar results can be shown.

1) Reverse Link Training: To learn all reverse links, the K sink nodes must transmit pilot

symbols over a period τp ≥ KNr. Similar to the analysis done in [15], we let the sinks collectively

transmit a KNr×τp unitary matrix of pilots Φ shown to be optimal in [28]. Such reverse channel

training requires no assumptions other than training sequences be known to all sink nodes, to

guarantee orthogonality. Of course, if we enforce orthogonality in time, this assumption reduces

to synchronization, and exact sequences need not be known to all sinks.

Let←−Y i =

�←−y i[1] . . . ←−

y i[τp]�

be the Nt × τp received training matrix at source node i, and

let←−Yp =

�←−Y

∗1,←−Y

∗2, . . . ,

←−Y

∗K

�∗be the composite received training matrix. We write the received

training as←−Yp =

�τpPf

NrGΦ+V,

where G is the KNt × KNr composite reverse channel matrix between all sinks and sources

and V is a KNt× τp matrix with i.i.d CN (0, σ2INr) elements. Using the received training, each

source derives an MMSE estimate of its channels all of which can be together written as

�G =

�τpPf

Nr

σ2 + τpPf

Nr

←−YpΦ

∗. (6)

For submission to IEEE Trans. Wireless Commun. 8

III. INTERFERENCE ALIGNMENT WITH ANALOG FEEDBACK

In this section we propose a feedback and transmission strategy based on analog feedback and

interference alignment, which uses the estimated channels as if they were the true propagation

channels.

A. Analog Feedback

To feedback the forward channel matrices Hi,� reliably across the feedback channel given in

(1), we propose dividing the feedback stage into two phases. First, each source learns its reverse

channels. Second, the forward channels are fed back and estimated. We neglect the initial training

phase in which sinks learn the forward channels and, thus, assume they have been estimated

perfectly. This is since imperfect CSI at the receiver will adversely affect performance of all

feedback schemes and is not exclusive to analog feedback. In fact, any error in estimating the

forward channels will only add an error term to the forward channels in (7), which also decays

with transmit power, and thus similar results can be shown.

1) Reverse Link Training: To learn all reverse links, the K sink nodes must transmit pilot

symbols over a period τp ≥ KNr. Similar to the analysis done in [15], we let the sinks collectively

transmit a KNr×τp unitary matrix of pilots Φ shown to be optimal in [28]. Such reverse channel

training requires no assumptions other than training sequences be known to all sink nodes, to

guarantee orthogonality. Of course, if we enforce orthogonality in time, this assumption reduces

to synchronization, and exact sequences need not be known to all sinks.

Let←−Y i =

�←−y i[1] . . . ←−

y i[τp]�

be the Nt × τp received training matrix at source node i, and

let←−Yp =

�←−Y

∗1,←−Y

∗2, . . . ,

←−Y

∗K

�∗be the composite received training matrix. We write the received

training as←−Yp =

�τpPf

NrGΦ+V,

where G is the KNt × KNr composite reverse channel matrix between all sinks and sources

and V is a KNt× τp matrix with i.i.d CN (0, σ2INr) elements. Using the received training, each

source derives an MMSE estimate of its channels all of which can be together written as

�G =

�τpPf

Nr

σ2 + τpPf

Nr

←−YpΦ

∗. (6)

For submission to IEEE Trans. Wireless Commun. 8

III. INTERFERENCE ALIGNMENT WITH ANALOG FEEDBACK

In this section we propose a feedback and transmission strategy based on analog feedback and

interference alignment, which uses the estimated channels as if they were the true propagation

channels.

A. Analog Feedback

To feedback the forward channel matrices Hi,� reliably across the feedback channel given in

(1), we propose dividing the feedback stage into two phases. First, each source learns its reverse

channels. Second, the forward channels are fed back and estimated. We neglect the initial training

phase in which sinks learn the forward channels and, thus, assume they have been estimated

perfectly. This is since imperfect CSI at the receiver will adversely affect performance of all

feedback schemes and is not exclusive to analog feedback. In fact, any error in estimating the

forward channels will only add an error term to the forward channels in (7), which also decays

with transmit power, and thus similar results can be shown.

1) Reverse Link Training: To learn all reverse links, the K sink nodes must transmit pilot

symbols over a period τp ≥ KNr. Similar to the analysis done in [15], we let the sinks collectively

transmit a KNr×τp unitary matrix of pilots Φ shown to be optimal in [28]. Such reverse channel

training requires no assumptions other than training sequences be known to all sink nodes, to

guarantee orthogonality. Of course, if we enforce orthogonality in time, this assumption reduces

to synchronization, and exact sequences need not be known to all sinks.

Let←−Y i =

�←−y i[1] . . . ←−

y i[τp]�

be the Nt × τp received training matrix at source node i, and

let←−Yp =

�←−Y

∗1,←−Y

∗2, . . . ,

←−Y

∗K

�∗be the composite received training matrix. We write the received

training as←−Yp =

�τpPf

NrGΦ+V,

where G is the KNt × KNr composite reverse channel matrix between all sinks and sources

and V is a KNt× τp matrix with i.i.d CN (0, σ2INr) elements. Using the received training, each

source derives an MMSE estimate of its channels all of which can be together written as

�G =

�τpPf

Nr

σ2 + τpPf

Nr

←−YpΦ

∗. (6)

Step 2: Train the reverse channel Send a pilot matrix

For submission to IEEE Trans. Wireless Commun. 8

III. INTERFERENCE ALIGNMENT WITH ANALOG FEEDBACK

In this section we propose a feedback and transmission strategy based on analog feedback and

interference alignment, which uses the estimated channels as if they were the true propagation

channels.

A. Analog Feedback

To feedback the forward channel matrices Hi,� reliably across the feedback channel given in

(1), we propose dividing the feedback stage into two phases. First, each source learns its reverse

channels. Second, the forward channels are fed back and estimated. We neglect the initial training

phase in which sinks learn the forward channels and, thus, assume they have been estimated

perfectly. This is since imperfect CSI at the receiver will adversely affect performance of all

feedback schemes and is not exclusive to analog feedback. In fact, any error in estimating the

forward channels will only add an error term to the forward channels in (7), which also decays

with transmit power, and thus similar results can be shown.

1) Reverse Link Training: To learn all reverse links, the K sink nodes must transmit pilot

symbols over a period τp ≥ KNr. Similar to the analysis done in [15], we let the sinks collectively

transmit a KNr×τp unitary matrix of pilots Φ shown to be optimal in [28]. Such reverse channel

training requires no assumptions other than training sequences be known to all sink nodes, to

guarantee orthogonality. Of course, if we enforce orthogonality in time, this assumption reduces

to synchronization, and exact sequences need not be known to all sinks.

Let←−Y i =

�←−y i[1] . . . ←−

y i[τp]�

be the Nt × τp received training matrix at source node i, and

let←−Yp =

�←−Y

∗1,←−Y

∗2, . . . ,

←−Y

∗K

�∗be the composite received training matrix. We write the received

training as←−Yp =

�τpPf

NrGΦ+V,

where G is the KNt × KNr composite reverse channel matrix between all sinks and sources

and V is a KNt× τp matrix with i.i.d CN (0, σ2INr) elements. Using the received training, each

source derives an MMSE estimate of its channels all of which can be together written as

�G =

�τpPf

Nr

σ2 + τpPf

Nr

←−YpΦ

∗. (6)

For submission to IEEE Trans. Wireless Commun. 8

III. INTERFERENCE ALIGNMENT WITH ANALOG FEEDBACK

In this section we propose a feedback and transmission strategy based on analog feedback and

interference alignment, which uses the estimated channels as if they were the true propagation

channels.

A. Analog Feedback

To feedback the forward channel matrices Hi,� reliably across the feedback channel given in

(1), we propose dividing the feedback stage into two phases. First, each source learns its reverse

channels. Second, the forward channels are fed back and estimated. We neglect the initial training

phase in which sinks learn the forward channels and, thus, assume they have been estimated

perfectly. This is since imperfect CSI at the receiver will adversely affect performance of all

feedback schemes and is not exclusive to analog feedback. In fact, any error in estimating the

forward channels will only add an error term to the forward channels in (7), which also decays

with transmit power, and thus similar results can be shown.

1) Reverse Link Training: To learn all reverse links, the K sink nodes must transmit pilot

symbols over a period τp ≥ KNr. Similar to the analysis done in [15], we let the sinks collectively

transmit a KNr×τp unitary matrix of pilots Φ shown to be optimal in [28]. Such reverse channel

training requires no assumptions other than training sequences be known to all sink nodes, to

guarantee orthogonality. Of course, if we enforce orthogonality in time, this assumption reduces

to synchronization, and exact sequences need not be known to all sinks.

Let←−Y i =

�←−y i[1] . . . ←−

y i[τp]�

be the Nt × τp received training matrix at source node i, and

let←−Yp =

�←−Y

∗1,←−Y

∗2, . . . ,

←−Y

∗K

�∗be the composite received training matrix. We write the received

training as←−Yp =

�τpPf

NrGΦ+V,

where G is the KNt × KNr composite reverse channel matrix between all sinks and sources

and V is a KNt× τp matrix with i.i.d CN (0, σ2INr) elements. Using the received training, each

source derives an MMSE estimate of its channels all of which can be together written as

�G =

�τpPf

Nr

σ2 + τpPf

Nr

←−YpΦ

∗. (6)

For submission to IEEE Trans. Wireless Commun. 8

III. INTERFERENCE ALIGNMENT WITH ANALOG FEEDBACK

In this section we propose a feedback and transmission strategy based on analog feedback and

interference alignment, which uses the estimated channels as if they were the true propagation

channels.

A. Analog Feedback

To feedback the forward channel matrices Hi,� reliably across the feedback channel given in

(1), we propose dividing the feedback stage into two phases. First, each source learns its reverse

channels. Second, the forward channels are fed back and estimated. We neglect the initial training

phase in which sinks learn the forward channels and, thus, assume they have been estimated

perfectly. This is since imperfect CSI at the receiver will adversely affect performance of all

feedback schemes and is not exclusive to analog feedback. In fact, any error in estimating the

forward channels will only add an error term to the forward channels in (7), which also decays

with transmit power, and thus similar results can be shown.

1) Reverse Link Training: To learn all reverse links, the K sink nodes must transmit pilot

symbols over a period τp ≥ KNr. Similar to the analysis done in [15], we let the sinks collectively

transmit a KNr×τp unitary matrix of pilots Φ shown to be optimal in [28]. Such reverse channel

training requires no assumptions other than training sequences be known to all sink nodes, to

guarantee orthogonality. Of course, if we enforce orthogonality in time, this assumption reduces

to synchronization, and exact sequences need not be known to all sinks.

Let←−Y i =

�←−y i[1] . . . ←−

y i[τp]�

be the Nt × τp received training matrix at source node i, and

let←−Yp =

�←−Y

∗1,←−Y

∗2, . . . ,

←−Y

∗K

�∗be the composite received training matrix. We write the received

training as←−Yp =

�τpPf

NrGΦ+V,

where G is the KNt × KNr composite reverse channel matrix between all sinks and sources

and V is a KNt× τp matrix with i.i.d CN (0, σ2INr) elements. Using the received training, each

source derives an MMSE estimate of its channels all of which can be together written as

�G =

�τpPf

Nr

σ2 + τpPf

Nr

←−YpΦ

∗. (6)Estimate G

Normalization: Lets not worry about it!!

TX k-1

TX k

TX k+1

RX k-1

RX k

RX k+1

21

2. Train reverse channel (Gik)

Friday, November 18, 2011

Page 22: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Analog Feedback

Step 3: Feedback CSI and estimate H

For submission to IEEE Trans. Wireless Commun. 9

Since �G is an MMSE estimate of Gaussian random variables corrupted by Gaussian noise, this

results in �G ∼ CN�0,

τpPfNr

σ2+τpPfNr

�and �G = G − �G with i.i.d CN

�0, σ2

σ2+τpPfNr

�elements. No

collaboration between sources is needed in (6); each node estimates its channels based on its

own received signals which we have concatenated in (6) and written in terms of the composite

received training matrix,←−Yp, to conserve space.

2) Analog CSI Feedback: After training the reverse channels, each sink node i directly sends

its unquantized uncoded estimates of Hi,� ∀ � over a period τc. To have the sink nodes feedback

their CSI simultaneously while maintaining orthogonality, each node post multiplies its Nr×KNt

feedback matrix [Hi,1 . . .Hi,K ] with a unitary KNt × τc matrix Ψ∗i such that Ψ∗

iΨ� = IKNtδi,�,

a general orthogonal feedback structure that can capture the case of orthogonality in time [15].

This requires τc ≥ KNtNr. The transmitted Nr × τc feedback matrix←−X i from node i can be

written as←−X i =

�τcPf

KNtNr[Hi,1 . . .Hi,K ]Ψ

∗i .

The concatenated, or composite received feedback KNt × τc matrix is then given by

←−Yc =

�τcPf

KNtNr

K�

i=1

Gi,1

...

Gi,K

[Hi,1 . . .Hi,K ]Ψ

∗i +V

where V is now a KNt × τc noise matrix.

To estimate the forward channels Hi,�, the source nodes first isolate the training from sink

node i by post multiplying their received training by Ψi which gives

←−YcΨi =

�τcPf

KNtNr

Gi,1

...

Gi,K

� �� �Gi

[Hi,1 . . .Hi,K ]� �� �Hi

+VΨi. (7)

To simplify the analysis in Section III-B, we assume sources share the complete matrix←−Yc, and

Feedback

Spread feedback on orthogonal sequences

For submission to IEEE Trans. Wireless Commun. 9

Since �G is an MMSE estimate of Gaussian random variables corrupted by Gaussian noise, this

results in �G ∼ CN�0,

τpPfNr

σ2+τpPfNr

�and �G = G − �G with i.i.d CN

�0, σ2

σ2+τpPfNr

�elements. No

collaboration between sources is needed in (6); each node estimates its channels based on its

own received signals which we have concatenated in (6) and written in terms of the composite

received training matrix,←−Yp, to conserve space.

2) Analog CSI Feedback: After training the reverse channels, each sink node i directly sends

its unquantized uncoded estimates of Hi,� ∀ � over a period τc. To have the sink nodes feedback

their CSI simultaneously while maintaining orthogonality, each node post multiplies its Nr×KNt

feedback matrix [Hi,1 . . .Hi,K ] with a unitary KNt × τc matrix Ψ∗i such that Ψ∗

iΨ� = IKNtδi,�,

a general orthogonal feedback structure that can capture the case of orthogonality in time [15].

This requires τc ≥ KNtNr. The transmitted Nr × τc feedback matrix←−X i from node i can be

written as←−X i =

�τcPf

KNtNr[Hi,1 . . .Hi,K ]Ψ

∗i .

The concatenated, or composite received feedback KNt × τc matrix is then given by

←−Yc =

�τcPf

KNtNr

K�

i=1

Gi,1

...

Gi,K

[Hi,1 . . .Hi,K ]Ψ

∗i +V

where V is now a KNt × τc noise matrix.

To estimate the forward channels Hi,�, the source nodes first isolate the training from sink

node i by post multiplying their received training by Ψi which gives

←−YcΨi =

�τcPf

KNtNr

Gi,1

...

Gi,K

� �� �Gi

[Hi,1 . . .Hi,K ]� �� �Hi

+VΨi. (7)

To simplify the analysis in Section III-B, we assume sources share the complete matrix←−Yc, and

Estimate H

For submission to IEEE Trans. Wireless Commun. 10

effectively compute a common least squares estimate �Hi of Hi given by

�Hi =

�KNtNr

τcPf

��G∗

i�Gi

�−1 �G∗i

←−YcΨi

= Hi����Real Channel

+ �Hi����Error

,

where �Gi is the estimate of Gi; a KNt ×Nr block taken from �G obtained in (6). It is realized

that such node cooperation is not practical; we return to this assumption at the end of this section

and provide alternative approaches that we show in Section V perform close to this special case.

The error in the estimates of Hi,� can then be written as

�Hi =��G∗

i�Gi

�−1 �G∗i

��KNtNr

τcPfV − �GiHi

�.

which makes it clear that the error in the estimate consists of two error terms: the first due to

noisy feedback and the second due to a noisy estimate of the feedback channel. To quantify the

effect of the error on the achieved sum rate, we derive the variance of the error term introduced

by analog feedback. Recall that the elements of Hi,� are CN (0, 1), those of V are CN (0, σ2),

and those of �Gi are CN (0, σ2

σ2+τpPfNr

). As a result, the error term �GiHi due to the reverse channel

estimation has independent elements with a variance of Nrσ2

σ2+τpPfNr

. Similarly to [15] we see that

the covariance of each columns of �Hi denoted �H(�)i , conditioned on �Gi is

Cov(�H(�)i |�Gi) =

�KNtNrσ2

τcPf+

Nrσ2

σ2 + τpPf

Nr

���G∗

i�Gi

�−1.

Since the elements of the MMSE estimate �Gi are Gaussian and uncorrelated, the diagonal

elements of��G∗

i�Gi

�−1are reciprocals of scaled chi-squared random variables with 2(KNt −

NR+1) degrees of freedom [15]. As a result, the mean square error, σ2f , in the elements of �Hi,�

is given by

σ2f =

σ2

(KNt −Nr)Pf

�N2

r

τp+

KNtNr

τc

�1 +

Nrσ2

τpPf

��. (8)

At high SNR this gives

σ2f ≈

σ2�

N2r

τp+ KNtNr

τc

(KNt −Nr)Pf. (9)Variance of error

Feedback Power

Feedback Time

SpreadingSequence

For submission to IEEE Trans. Wireless Commun. 8

III. INTERFERENCE ALIGNMENT WITH ANALOG FEEDBACK

In this section we propose a feedback and transmission strategy based on analog feedback and

interference alignment, which uses the estimated channels as if they were the true propagation

channels.

A. Analog Feedback

To feedback the forward channel matrices Hi,� reliably across the feedback channel given in

(1), we propose dividing the feedback stage into two phases. First, each source learns its reverse

channels. Second, the forward channels are fed back and estimated. We neglect the initial training

phase in which sinks learn the forward channels and, thus, assume they have been estimated

perfectly. This is since imperfect CSI at the receiver will adversely affect performance of all

feedback schemes and is not exclusive to analog feedback. In fact, any error in estimating the

forward channels will only add an error term to the forward channels in (7), which also decays

with transmit power, and thus similar results can be shown.

1) Reverse Link Training: To learn all reverse links, the K sink nodes must transmit pilot

symbols over a period τp ≥ KNr. Similar to the analysis done in [15], we let the sinks collectively

transmit a KNr×τp unitary matrix of pilots Φ shown to be optimal in [28]. Such reverse channel

training requires no assumptions other than training sequences be known to all sink nodes, to

guarantee orthogonality. Of course, if we enforce orthogonality in time, this assumption reduces

to synchronization, and exact sequences need not be known to all sinks.

Let←−Y i =

�←−y i[1] . . . ←−

y i[τp]�

be the Nt × τp received training matrix at source node i, and

let←−Yp =

�←−Y

∗1,←−Y

∗2, . . . ,

←−Y

∗K

�∗be the composite received training matrix. We write the received

training as←−Yp =

�τpPf

NrGΦ+V,

where G is the KNt × KNr composite reverse channel matrix between all sinks and sources

and V is a KNt× τp matrix with i.i.d CN (0, σ2INr) elements. Using the received training, each

source derives an MMSE estimate of its channels all of which can be together written as

�G =

�τpPf

Nr

σ2 + τpPf

Nr

←−YpΦ

∗. (6)

For submission to IEEE Trans. Wireless Commun. 9

Since �G is an MMSE estimate of Gaussian random variables corrupted by Gaussian noise, this

results in �G ∼ CN�0,

τpPfNr

σ2+τpPfNr

�and �G = G − �G with i.i.d CN

�0, σ2

σ2+τpPfNr

�elements. No

collaboration between sources is needed in (6); each node estimates its channels based on its

own received signals which we have concatenated in (6) and written in terms of the composite

received training matrix,←−Yp, to conserve space.

2) Analog CSI Feedback: After training the reverse channels, each sink node i directly sends

its unquantized uncoded estimates of Hi,� ∀ � over a period τc. To have the sink nodes feedback

their CSI simultaneously while maintaining orthogonality, each node post multiplies its Nr×KNt

feedback matrix [Hi,1 . . .Hi,K ] with a unitary KNt × τc matrix Ψ∗i such that Ψ∗

iΨ� = IKNtδi,�,

a general orthogonal feedback structure that can capture the case of orthogonality in time [15].

This requires τc ≥ KNtNr. The transmitted Nr × τc feedback matrix←−X i from node i can be

written as←−X i =

�τcPf

KNtNr[Hi,1 . . .Hi,K ]Ψ

∗i .

The concatenated, or composite received feedback KNt × τc matrix is then given by

←−Yc =

�τcPf

KNtNr

K�

i=1

Gi,1

...

Gi,K

[Hi,1 . . .Hi,K ]Ψ

∗i +V

where V is now a KNt × τc noise matrix.

To estimate the forward channels Hi,�, the source nodes first isolate the training from sink

node i by post multiplying their received training by Ψi which gives

←−YcΨi =

�τcPf

KNtNr

Gi,1

...

Gi,K

� �� �Gi

[Hi,1 . . .Hi,K ]� �� �Hi

+VΨi. (7)

To simplify the analysis in Section III-B, we assume sources share the complete matrix←−Yc, and

For submission to IEEE Trans. Wireless Commun. 9

Since �G is an MMSE estimate of Gaussian random variables corrupted by Gaussian noise, this

results in �G ∼ CN�0,

τpPfNr

σ2+τpPfNr

�and �G = G − �G with i.i.d CN

�0, σ2

σ2+τpPfNr

�elements. No

collaboration between sources is needed in (6); each node estimates its channels based on its

own received signals which we have concatenated in (6) and written in terms of the composite

received training matrix,←−Yp, to conserve space.

2) Analog CSI Feedback: After training the reverse channels, each sink node i directly sends

its unquantized uncoded estimates of Hi,� ∀ � over a period τc. To have the sink nodes feedback

their CSI simultaneously while maintaining orthogonality, each node post multiplies its Nr×KNt

feedback matrix [Hi,1 . . .Hi,K ] with a unitary KNt × τc matrix Ψ∗i such that Ψ∗

iΨ� = IKNtδi,�,

a general orthogonal feedback structure that can capture the case of orthogonality in time [15].

This requires τc ≥ KNtNr. The transmitted Nr × τc feedback matrix←−X i from node i can be

written as←−X i =

�τcPf

KNtNr[Hi,1 . . .Hi,K ]Ψ

∗i .

The concatenated, or composite received feedback KNt × τc matrix is then given by

←−Yc =

�τcPf

KNtNr

K�

i=1

Gi,1

...

Gi,K

[Hi,1 . . .Hi,K ]Ψ

∗i +V

where V is now a KNt × τc noise matrix.

To estimate the forward channels Hi,�, the source nodes first isolate the training from sink

node i by post multiplying their received training by Ψi which gives

←−YcΨi =

�τcPf

KNtNr

Gi,1

...

Gi,K

� �� �Gi

[Hi,1 . . .Hi,K ]� �� �Hi

+VΨi. (7)

To simplify the analysis in Section III-B, we assume sources share the complete matrix←−Yc, and

TX k-1

TX k

TX k+1

RX k-1

RX k

RX k+1

22

3. Feedback Hik as QAM symbols

Friday, November 18, 2011

Page 23: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Performance Analysis

For submission to IEEE Trans. Wireless Commun. 11

Having computed feedback error, we return to the assumption on node cooperation. Source

cooperation simplifies analysis by effectively making a common �H known to all users, however,

this cannot be assumed in practical systems. The first alternative to cooperation is to have one user

calculate precoders and combiners, and feed them forward to the other users as in [29]. The other

alternative is for each node to calculate its precoder and combiner based on the feedback only

it receives. We refer to this as the “distributed processing” approach. Due to independent noise

across users in the feedback stage, distributed processing implies that users will calculate vectors

based on mismatched information, the effect of which can also be bounded. We elaborate on

this after presenting our results in Sections III-B and V and show that the cost of no cooperation

is limited.

B. Multiplexing Gain with Analog Feedback

To characterize the performance of interference alignment with analog feedback, we examine

the mean loss in sum-rate [26] incurred by naive interference alignment where the channel

estimates obtained at the sources are used, as if they were the true channels, to calculate the

columns of the precoders, fmi ∀i,m and combiners w

mi ∀i,m. Transmit and receive vectors are

calculated to satisfy (2), (3), and (4) using the estimated channels.

The mean loss in sum-rate is defined as

∆Rsum � EHRsum − EHRsum, (10)

where EHRsum is the average sum rate from interference alignment with perfect CSI, with

instantaneous rate given in (5), and EHRsum is the average rate with CSI obtained via feedback.

Theorem 1: Interference alignment on the K-user Nr×Nt interference channel with imperfect

channel state information obtained via the analog feedback strategy described in Section III-A

achieves the same average sum-rate scaling observed with perfect interference alignment as long

as the feedback power Pf scales with the transmit power P . Thus, the original multiplexing gain

is preserved. Moreover, the mean loss in sum rate ∆Rsum is O(1).

Proof: Let the K-user Nr × Nt interference channel use the analog feedback scheme

presented to achieve a vector of multiplexing gains d. Let the transmit precoding and receive

Rate with perfect CSI Rate with CSI from feedback

For submission to IEEE Trans. Wireless Commun. 12

combining vectors be calculated to satisfy

(wmi )

∗ �Hi,if�i = 0, ∀i, � �= m (11)

(wmi )

∗ �Hi,k f�k = 0, ∀k �= i, and ∀m, � (12)

���(wmi )

∗ �Hi,ifmi

��� ≥ c > 0, ∀i,m. (13)

Using these precoding and combining vectors, the input-output relationship at the output of a

linear zero-forcing receiver is

(wmi )

∗yi = (wm

i )∗√PHi,if

mi smi +

� �=m

(wmi )

∗√PHi,if

�i s

�i+

k �=i

dk�

�=1

(wmi )

∗√PHi,k f

�ks

�k+(wm

i )∗vi.

(14)

Using the received signal in (2), the instantaneous rate expression in (5) , and the sum rate loss

defined in (10), this gives the following upper bound on mean loss in sum rate:

∆Rsum = EH

i,m

log2

�1 +

Pdi|(wm

i )∗Hi,if

mi |2

σ2

�− E

H,H

i,m

log2

1 +

Pdi

���(wmi )

∗Hi,if

mi

���2

I1i,m + I2

i,m� �� �Ii,m

+σ2

= EH

i,m

log2

�1 +

Pdi|(wm

i )∗Hi,if

mi |2

σ2

− EH,H

i,m

log2

1 +Ii,m + P

di

���(wmi )

∗Hi,if

mi

���2

σ2

+ EH,H

i,m

log2

�1 +

Ii,m

σ2

(a)

≤ EHH

i,m

log2

�1 +

Ii,m

σ2

where (a) is due to the fact that wmi , wm

i , fmi , and fmi are independent of Hi,i and therefore

Pdi

���(wmi )

∗Hi,if

mi

���2

and Pdi|(wm

i )∗Hi,if

mi |2 are identically distributed. As a result P

di

���(wmi )

∗Hi,if

mi

���2+

Ii,m stochastically dominates Pdi|(wm

i )∗Hi,if

mi |2 [26]. We now apply Jensen’s inequality to the

Derived with perfect CSIDerived with imperfect CSI

Feedback error =

Leakage Interference

Recall: No interferenceIA conditions

23

What happens to sum rate?

What is the mean sum rate loss relative to perfect CSI?

Friday, November 18, 2011

Page 24: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

For submission to IEEE Trans. Wireless Commun. 13

upper bound in (a) to get

∆Rsum ≤�

i,m

log2

�1 +

EH,HIi,m

σ2

�. (15)

Since (11), (12), (13) are satisfied, however, the total interference term Ii,m can be simplified to

include only residual interference due to the channel estimation errors �Hi,�. Equation (15) can

be further upper bounded by noticing that���(wm

i )∗(Hi,k + �Hi,k)f

�k

���2=

���(wmi )

∗(�Hi,k)f�k

���2

∀k, ∀� �= m

≤ ��Hi,k�2F ,(16)

which gives

∆Rsum ≤�

i,m

log2

1 +

K��=1

Pd�

(d� − δi,�)EH,H�Hi,��2F

σ2

. (17)

From (8), however, we have EH,H�Hi,��2F = NtNrσ2

f = cPf

, where c is a constant, independent

of Pf at high enough SNR, given by

c = NtNr

σ2�

N2r

τp+ KNtNr(1+�)

τc

(KNt −Nr). (18)

Combining (17) and (18) gives the final upper bound on throughput loss

∆Rsum ≤�

i,m

log2

1 +P

σ2di

(di − 1)c

Pf+

� �=i

d�c

Pf

.

Therefore, if we have Pf = α−1P we get

∆Rsum ≤�

i

di log2

�1 +

(�d�1 − 1)αc

diσ2

�. (19)

The bound has been presented at high SNR for simplicity of exposition only; (18) can be adapted

for any SNR > 0 by using (8) instead of (9).

In summary, Theorem 1 states that if feedback power is equal to any constant fraction of

transmit power, the cost of imperfect CSI is a constant number of bits, independent of SNR. Since

transmit and feedback power are likely to be comparable in practice, this result is promising.

Using analog feedback allows the system to overcome the problem of exploding complexity and

For submission to IEEE Trans. Wireless Commun. 12

combining vectors be calculated to satisfy

(wmi )

∗ �Hi,if�i = 0, ∀i, � �= m (11)

(wmi )

∗ �Hi,k f�k = 0, ∀k �= i, and ∀m, � (12)

���(wmi )

∗ �Hi,ifmi

��� ≥ c > 0, ∀i,m. (13)

Using these precoding and combining vectors, the input-output relationship at the output of a

linear zero-forcing receiver is

(wmi )

∗yi = (wm

i )∗√PHi,if

mi smi +

� �=m

(wmi )

∗√PHi,if

�i s

�i+

k �=i

dk�

�=1

(wmi )

∗√PHi,k f

�ks

�k+(wm

i )∗vi.

(14)

Using the received signal in (2), the instantaneous rate expression in (5) , and the sum rate loss

defined in (10), this gives the following upper bound on mean loss in sum rate:

∆Rsum = EH

i,m

log2

�1 +

Pdi|(wm

i )∗Hi,if

mi |2

σ2

�− E

H,H

i,m

log2

1 +

Pdi

���(wmi )

∗Hi,if

mi

���2

I1i,m + I2

i,m� �� �Ii,m

+σ2

= EH

i,m

log2

�1 +

Pdi|(wm

i )∗Hi,if

mi |2

σ2

− EH,H

i,m

log2

1 +Ii,m + P

di

���(wmi )

∗Hi,if

mi

���2

σ2

+ EH,H

i,m

log2

�1 +

Ii,m

σ2

(a)

≤ EHH

i,m

log2

�1 +

Ii,m

σ2

where (a) is due to the fact that wmi , wm

i , fmi , and fmi are independent of Hi,i and therefore

Pdi

���(wmi )

∗Hi,if

mi

���2

and Pdi|(wm

i )∗Hi,if

mi |2 are identically distributed. As a result P

di

���(wmi )

∗Hi,if

mi

���2+

Ii,m stochastically dominates Pdi|(wm

i )∗Hi,if

mi |2 [26]. We now apply Jensen’s inequality to the

For submission to IEEE Trans. Wireless Commun. 13

upper bound in (a) to get

∆Rsum ≤�

i,m

log2

�1 +

EH,HIi,m

σ2

�. (15)

Since (11), (12), (13) are satisfied, however, the total interference term Ii,m can be simplified to

include only residual interference due to the channel estimation errors �Hi,�. Equation (15) can

be further upper bounded by noticing that���(wm

i )∗(Hi,k + �Hi,k)f

�k

���2=

���(wmi )

∗(�Hi,k)f�k

���2

∀k, ∀� �= m

≤ ��Hi,k�2F ,(16)

which gives

∆Rsum ≤�

i,m

log2

1 +

K��=1

Pd�

(d� − δi,�)EH,H�Hi,��2F

σ2

. (17)

From (8), however, we have EH,H�Hi,��2F = NtNrσ2

f = cPf

, where c is a constant, independent

of Pf at high enough SNR, given by

c = NtNr

σ2�

N2r

τp+ KNtNr(1+�)

τc

(KNt −Nr). (18)

Combining (17) and (18) gives the final upper bound on throughput loss

∆Rsum ≤�

i,m

log2

1 +P

σ2di

(di − 1)c

Pf+

� �=i

d�c

Pf

.

Therefore, if we have Pf = α−1P we get

∆Rsum ≤�

i

di log2

�1 +

(�d�1 − 1)αc

diσ2

�. (19)

The bound has been presented at high SNR for simplicity of exposition only; (18) can be adapted

for any SNR > 0 by using (8) instead of (9).

In summary, Theorem 1 states that if feedback power is equal to any constant fraction of

transmit power, the cost of imperfect CSI is a constant number of bits, independent of SNR. Since

transmit and feedback power are likely to be comparable in practice, this result is promising.

Using analog feedback allows the system to overcome the problem of exploding complexity and

What happens to sum rate?

What is the mean sum rate loss relative to perfect CSI?

24

Performance Analysis

Sum rate loss is a function of interference power

Feedback power counteracts interference

power

For Pfb =O(P)Mux. gain preserved!!

Constant loss in sum rate!!

Friday, November 18, 2011

Page 25: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Simulation Results

Perfe

ct CSI

Pfb =P/2

Pfb =constant

Multiplexing gainpreserved

Multiplexing gain = 0

SNR independent gap & Tight lower bound

• 3 User• 2x2 MIMO links• Rayleigh fading

25

0 5 10 15 20 25 30 35 400

10

20

30

40

50

60

70

SNR (dB)

Su

m R

ate

(b

its/s

/Hz)

IA with Perfect CSI

IA with Pf=P/2

IA with Dist. Processing & Pf=P/2

Approx. Lower Bound (with c2(!

p!c) )

Lower Bound (with c(!p!c) )

IA with Pf = P

0.5

IA with Fixed Feedback Quality

Partial multiplexing gains achieved

Friday, November 18, 2011

Page 26: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Part I: SummaryAnalog feedback enables interference alignment

Flexible approach for providing high resolution channel estimates

Gives a scalable multiplexing optimal feedback strategy

Can be improved through optimization of various parameters

Further work

Feedback in partial interference channels

Analysis of distributed performance with uncertainty

26

Friday, November 18, 2011

Page 27: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

OutlineIntroduction to interference alignment

Analog feedback for interference alignment

Limited feedback for interference alignment

Conclusions

27

Friday, November 18, 2011

Page 28: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

OFDM System ModelK-user wide band interference channel, single antenna

Transmitters use OFDM to obtain virtual MIMO channels

Time domain channel taps are temporally correlated

28

TX

k-1

Hk-1,k-1

Hk,k

Hk+1,k+1

Hk,k-1

Hk+1,k

Hk-1,k

Hk,k+1Hk-1

,k+1

Hk+1,k-1

TX

k

TX

k+1

RX

k-1

RX

k

RX

k+1

Figure 1: K-User SISO interference channel model

matrix input-output relationship is

yk[t] = Hk,k[t]xk[t] +!

! !=k

Hk,![t]x![t] + vk[t]. (1)

where xk[t] is the OFDM symbol sent by userk at time t with the average power constraintE"!xk[t]!2

#= NscP , the Nsc " Nsc matrix

Hk,![t] = diag$FNsc

%h"k,![t], 0Nsc#L

&"'represents

the channel frequency response between transmitter ! andreceiver k at time t, and vk[t] is the i.i.d. CN (0,"2

nINsc)thermal noise observed by user k. The system modelassumes perfect time and frequency synchronization, and acyclic prefix longer than all channel responses.

The channels seen by the t-th OFDM symbol are as-sumed to remain constant. The channels over consecutiveOFDM symbols, however, are temporally correlated suchthat E [|hk,![t+ 1]hk,![t]|] = # # 1.

3. INTERFERENCE ALIGNMENT IN FREQUENCY

In this section we review the concept of IA over frequencyextensions and summarize the effect of imperfect CSI on theperformance of IA.

3.1 SISO IA via Frequency ExtensionsIA for the SISO interference channel can achieve the max-imum degrees of freedom when coding over infinite chan-nel extensions [1]. Using IA over Nsc frequency extensions,each transmitter k sends dk < Nsc symbols along the pre-coding vectors fmk [t]. As a result the transmitted symbol is

xk[t] =dk!

m=1

fmk [t]xmk [t], (2)

where xmk [t] are the symbols transmitted by user k at time

t. To satisfy the power constraint, we set !fmk [t]!2 = 1,and E

"xmk [t]2

#= NscP/dk, such that the total power in

each Nsc subcarriers is NscP . The transmit directions fmk [t]are calculated such that the interference from K $ 1 users isaligned at all receivers, leaving interference free dimensionsfor the desired signal.

In this section, we restrict our attention to IA with a zero-forcing receiver. At the output of the linear receiver, the re-ceived signal is

wmk [t]"yk[t] =wm

k [t]"Hk,k[t]fmk [t]xm

k [t]

+!

(i,!) !=(k,m)

wmk [t]"Hk,i[t]f

!i [t]x

!i [t]

+wmk [t]"vk[t],

(3)

for m % {1, . . . , dk} and k % {1, . . . ,K}, where!wm

k [t]!2 = 1. With a linear receiver, the conditions forperfect IA can be restated as

wmk [t]"Hk,i[t]f

!k[t] = 0, (i, !) &= (k,m) (4)

|wmk [t]"Hk,k[t]f

mk [t]| ' c > 0, (k,m (5)

where alignment is achieved by (4), and (5) ensures the de-codability of the dk desired streams.

The achievability proof in [1] showed that if fading is in-dependent on all subcarriers, then the vectors f !k[t] can befound to satisfy (4) and (5), if dk’s are chosen as in [1].Fortunately, [5] has claimed that fading on each subcarrierneed not be independent provided that the channel impulseresponse is long enough.

3.2 The Effect of Imperfect CSI FeedbackWith imperfect or limited CSI feedback, condition (4) is notsatisfied, resulting in residual interference. As the trans-mit power increases, so does the leakage interference power,which saturates the sum rate at high SNR.

In [6], it is shown that if imperfect CSI is used to calculatethe IA precoders and combiners, (fmk [t] and (wm

k [t], the meanloss in sum rate is upper bounded by

!Rsum #!

k,m

1

Nsclog2

)

*1 +EH

%I1k,m + I2

k,m

&

"2

+

, , (6)

where I1i,m[t] =

-! !=m

NscPdk

... (wmk [t]"Hk,k[t](f !k[t]

...2, and

I2i,m[t] =

-i !=k

di-!=1

NscPdi

... (wmk [t]"Hk,i[t](f !i [t]

...2

are the inter-

stream and inter-user interference respectively. The objec-tive of the feedback algorithm then becomes minimizing thetotal leakage interference by improving effective CSI accu-racy. Using the result from [5], the individual interferenceterms can be upper bounded as

!!! "wmk [t]!Hk,i[t]"f !i [t]

!!!2!

""wmk [t] # "f !i [t]"2"hk,i[t]"2

#1$

!!!!!hk,i[t]

!"hk,i[t]

"hk,i[t]"""hk,i[t]"

!!!!!

2$.

(7)

From (7) we see that leakage interference is directly relatedto the angle between the normalized channel impulse re-sponse, hk,i[t]

$hk,i[t]$ , and its quantized version!hk,i[t]

$!hk,i[t]$. As a

result, to limit performance degradation an efficient feedbackstrategy must attempt to minimize the angle between the ac-tual and quantized channels.

Desired Signal Interference

TX

k-1

Hk-1,k-1

Hk,k

Hk+1,k+1

Hk,k-1

Hk+1,k

Hk-1,k

Hk,k+1Hk-1

,k+1

Hk+1,k-1

TX

k

TX

k+1

RX

k-1

RX

k

RX

k+1

Figure 1: K-User SISO interference channel model

matrix input-output relationship is

yk[t] = Hk,k[t]xk[t] +!

! !=k

Hk,![t]x![t] + vk[t]. (1)

where xk[t] is the OFDM symbol sent by userk at time t with the average power constraintE"!xk[t]!2

#= NscP , the Nsc " Nsc matrix

Hk,![t] = diag$FNsc

%h"k,![t], 0Nsc#L

&"'represents

the channel frequency response between transmitter ! andreceiver k at time t, and vk[t] is the i.i.d. CN (0,"2

nINsc)thermal noise observed by user k. The system modelassumes perfect time and frequency synchronization, and acyclic prefix longer than all channel responses.

The channels seen by the t-th OFDM symbol are as-sumed to remain constant. The channels over consecutiveOFDM symbols, however, are temporally correlated suchthat E [|hk,![t+ 1]hk,![t]|] = # # 1.

3. INTERFERENCE ALIGNMENT IN FREQUENCY

In this section we review the concept of IA over frequencyextensions and summarize the effect of imperfect CSI on theperformance of IA.

3.1 SISO IA via Frequency ExtensionsIA for the SISO interference channel can achieve the max-imum degrees of freedom when coding over infinite chan-nel extensions [1]. Using IA over Nsc frequency extensions,each transmitter k sends dk < Nsc symbols along the pre-coding vectors fmk [t]. As a result the transmitted symbol is

xk[t] =dk!

m=1

fmk [t]xmk [t], (2)

where xmk [t] are the symbols transmitted by user k at time

t. To satisfy the power constraint, we set !fmk [t]!2 = 1,and E

"xmk [t]2

#= NscP/dk, such that the total power in

each Nsc subcarriers is NscP . The transmit directions fmk [t]are calculated such that the interference from K $ 1 users isaligned at all receivers, leaving interference free dimensionsfor the desired signal.

In this section, we restrict our attention to IA with a zero-forcing receiver. At the output of the linear receiver, the re-ceived signal is

wmk [t]"yk[t] =wm

k [t]"Hk,k[t]fmk [t]xm

k [t]

+!

(i,!) !=(k,m)

wmk [t]"Hk,i[t]f

!i [t]x

!i [t]

+wmk [t]"vk[t],

(3)

for m % {1, . . . , dk} and k % {1, . . . ,K}, where!wm

k [t]!2 = 1. With a linear receiver, the conditions forperfect IA can be restated as

wmk [t]"Hk,i[t]f

!k[t] = 0, (i, !) &= (k,m) (4)

|wmk [t]"Hk,k[t]f

mk [t]| ' c > 0, (k,m (5)

where alignment is achieved by (4), and (5) ensures the de-codability of the dk desired streams.

The achievability proof in [1] showed that if fading is in-dependent on all subcarriers, then the vectors f !k[t] can befound to satisfy (4) and (5), if dk’s are chosen as in [1].Fortunately, [5] has claimed that fading on each subcarrierneed not be independent provided that the channel impulseresponse is long enough.

3.2 The Effect of Imperfect CSI FeedbackWith imperfect or limited CSI feedback, condition (4) is notsatisfied, resulting in residual interference. As the trans-mit power increases, so does the leakage interference power,which saturates the sum rate at high SNR.

In [6], it is shown that if imperfect CSI is used to calculatethe IA precoders and combiners, (fmk [t] and (wm

k [t], the meanloss in sum rate is upper bounded by

!Rsum #!

k,m

1

Nsclog2

)

*1 +EH

%I1k,m + I2

k,m

&

"2

+

, , (6)

where I1i,m[t] =

-! !=m

NscPdk

... (wmk [t]"Hk,k[t](f !k[t]

...2, and

I2i,m[t] =

-i !=k

di-!=1

NscPdi

... (wmk [t]"Hk,i[t](f !i [t]

...2

are the inter-

stream and inter-user interference respectively. The objec-tive of the feedback algorithm then becomes minimizing thetotal leakage interference by improving effective CSI accu-racy. Using the result from [5], the individual interferenceterms can be upper bounded as

!!! "wmk [t]!Hk,i[t]"f !i [t]

!!!2!

""wmk [t] # "f !i [t]"2"hk,i[t]"2

#1$

!!!!!hk,i[t]

!"hk,i[t]

"hk,i[t]"""hk,i[t]"

!!!!!

2$.

(7)

From (7) we see that leakage interference is directly relatedto the angle between the normalized channel impulse re-sponse, hk,i[t]

$hk,i[t]$ , and its quantized version!hk,i[t]

$!hk,i[t]$. As a

result, to limit performance degradation an efficient feedbackstrategy must attempt to minimize the angle between the ac-tual and quantized channels.

Revision 3 on August 22, 2011 6

Throughout this paper, we consider temporally correlated Gaussian channels according to the

P -order autoregressive channel model defined in [23]. For such processes, each time series is

generated as

hk,![t] =P!

m=1

!mhk,![t!m] + "zk,![t], (3)

where zk,l[t] = CN (0,Rhk,!). The coefficients !m and " are calculated by fitting an autoregressive

model to the original Doppler spectrum proposed by Clarke [24] given by E[|h!k,!(t!m)hk,![t]|] =

J0(2#fDTsm) where fDTs is the channel’s normalized Doppler spread1 and J0 is the 0-th order

Bessel function of the first kind [25]. In some cases, for simplicity of exposition, we restrict our

attention to the special case of the first order autoregressive model where channels follow the

relation hk,![t] = "fhk,![t!m] +"1! "2fzk,![t] with "f = J0(2#fDTs).

III. INTERFERENCE ALIGNMENT IN FREQUENCY

In this section we review the concept of IA over frequency extensions when perfect channel

state information is available at the transmitter (CSIT), and summarize the effect of imperfect

CSIT on the performance of IA.

A. IA with Perfect CSI at the Transmitter

IA for the SISO interference channel can achieve the maximum degrees of freedom defined aslim

P"#1

Nsc

Rsum

log2 P= K

2 when coding over infinite channel extensions [2]. Using IA over N frequency

extensions, each transmitter k at time t sends dk < N symbols, xdk[t], along the N "1 precoding

vectors fk,d[t]. As a result, the input-output relation is

xk[t] =dk!

d=1

fk,d[t]xdk[t], (4)

where #fk,d[t]#2 = 1, and E#|xd

k[t]|2$= NP/dk, such that the total power in each N subcarriers

is NP . The transmit directions fk,d[t] are calculated such that the interference from K ! 1 users

is aligned at all receivers, leaving interference free dimensions for the desired signal.

Note that users need not code over all Nsc subcarriers, due to increasing algorithm complexity

and marginal gains with increasing dimensions. Users could potentially treat each N subcarriers

1fD is the channel’s Doppler spread, and Ts is the feedback interval, Ts = 1 OFDM symbol time, for example

Revision 3 on August 22, 2011 6

Throughout this paper, we consider temporally correlated Gaussian channels according to the

P -order autoregressive channel model defined in [23]. For such processes, each time series is

generated as

hk,![t] =P!

m=1

!mhk,![t!m] + "zk,![t], (3)

where zk,l[t] = CN (0,Rhk,!). The coefficients !m and " are calculated by fitting an autoregressive

model to the original Doppler spectrum proposed by Clarke [24] given by E[|h!k,!(t!m)hk,![t]|] =

J0(2#fDTsm) where fDTs is the channel’s normalized Doppler spread1 and J0 is the 0-th order

Bessel function of the first kind [25]. In some cases, for simplicity of exposition, we restrict our

attention to the special case of the first order autoregressive model where channels follow the

relation hk,![t] = "fhk,![t!m] +"1! "2fzk,![t] with "f = J0(2#fDTs).

III. INTERFERENCE ALIGNMENT IN FREQUENCY

In this section we review the concept of IA over frequency extensions when perfect channel

state information is available at the transmitter (CSIT), and summarize the effect of imperfect

CSIT on the performance of IA.

A. IA with Perfect CSI at the Transmitter

IA for the SISO interference channel can achieve the maximum degrees of freedom defined aslim

P"#1

Nsc

Rsum

log2 P= K

2 when coding over infinite channel extensions [2]. Using IA over N frequency

extensions, each transmitter k at time t sends dk < N symbols, xdk[t], along the N "1 precoding

vectors fk,d[t]. As a result, the input-output relation is

xk[t] =dk!

d=1

fk,d[t]xdk[t], (4)

where #fk,d[t]#2 = 1, and E#|xd

k[t]|2$= NP/dk, such that the total power in each N subcarriers

is NP . The transmit directions fk,d[t] are calculated such that the interference from K ! 1 users

is aligned at all receivers, leaving interference free dimensions for the desired signal.

Note that users need not code over all Nsc subcarriers, due to increasing algorithm complexity

and marginal gains with increasing dimensions. Users could potentially treat each N subcarriers

1fD is the channel’s Doppler spread, and Ts is the feedback interval, Ts = 1 OFDM symbol time, for example

OFDM symbol

Note: Feedback strategy will be independent of correlation model

Can now use closed form solutions for diagonal channels or run our favorite algorithm

Friday, November 18, 2011

Page 29: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Problem Statement

29

Exploit time variation to improve CSI quality

Receivers quantize time varying channel

Transmitters calculate the transmit directions

Receivers feedback the quantized channels

Transmitters send payload data

Prior work on differential feedback:Differential feedback of input covariances [SacKal09]

Predictive Grassmannian feedback for MISO broadcast [InoHea11]

Solution: Differential quantization and feedback

Friday, November 18, 2011

Page 30: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

What CSI is actually needed?

Loss in sum rate due to imperfect CSI

Manifold Structure

30

Revision 3 on August 22, 2011 9

design of efficient, low overhead feedback strategies that limit the effect of imperfect CSI and

allow good IA performance.

In [12], [13], it was shown that quantizing wideband channel taps, in a manner similar to single

user MIMO quantization [28], and scaling the number of feedback bits with SNR, preserves the

system’s multiplexing gain. Scaling feedback bits, and effectively codebook size, however, leads

to increasing quantization complexity at high SNR. In [14] the authors showed that analog

feedback can be used instead of quantization based algorithms to maintain the IA sum rate

perfomance while keeping feedback complexity and overhead constant at all SNRs. This is

made possible by the fact that, in analog feedback, CSI quality improves automatically with

SNR and keeps leakage interference power bounded in the desired signal space.

While both these feedback strategies propose valid ways to maintain system performance, both

neglect the temporal correlation of the channel. This results in the inefficient use of feedback

resources. The analysis done in [12]–[14], [29], however, extends to arbitrary feedback strategies

and we therefore summarize their findings to motivate the problem and refer the reader to [12],

[14] for a more detailed derivation.

In [12], [14], it is shown that if the transmitters use the estimates ofHk,!, !k, ! " {1, 2, . . . , K},

denoted by !Hk,!, to calculate IA precoders, !fmk [t], and combiners, !wmk [t], then the mean loss in

sum rate due to imperfect CSI can be upper bounded by [14],

!Rsum #"

k,m

1

Nsclog2

#

1 +EH

$I1k,m + I2

k,m

%

"2

&

. (13)

The objective then becomes minimizing the sum leakage interference term I1k,m + I2

k,m. Using

the result from [12], the individual interference terms can be upper bounded by

NscP

dk

'''!wmk [t]

!Hk,i[t]!f !i [t]'''2=NscP

dk

'''(!wm

k [t] $!f!i [t]

)!hfk,i[t]

'''2

#NscP

dk%!wm

k [t] $ !f !i [t]%2%hfk,i[t]%

2

*

+1&

'''''hfk,i[t]

!!hfk,i[t]

%hfk,i[t]%%!h

fk,i[t]%

'''''

2,

-

=NscP

dk%!wm

k [t] $ !f !i [t]%2%hk,i[t]%2

*

+1&

'''''hk,i[t]!!hk,i[t]

%hk,i[t]%%!hk,i[t]%

'''''

2,

- ,

(14)

Revision 3 on August 22, 2011 9

design of efficient, low overhead feedback strategies that limit the effect of imperfect CSI and

allow good IA performance.

In [12], [13], it was shown that quantizing wideband channel taps, in a manner similar to single

user MIMO quantization [28], and scaling the number of feedback bits with SNR, preserves the

system’s multiplexing gain. Scaling feedback bits, and effectively codebook size, however, leads

to increasing quantization complexity at high SNR. In [14] the authors showed that analog

feedback can be used instead of quantization based algorithms to maintain the IA sum rate

perfomance while keeping feedback complexity and overhead constant at all SNRs. This is

made possible by the fact that, in analog feedback, CSI quality improves automatically with

SNR and keeps leakage interference power bounded in the desired signal space.

While both these feedback strategies propose valid ways to maintain system performance, both

neglect the temporal correlation of the channel. This results in the inefficient use of feedback

resources. The analysis done in [12]–[14], [29], however, extends to arbitrary feedback strategies

and we therefore summarize their findings to motivate the problem and refer the reader to [12],

[14] for a more detailed derivation.

In [12], [14], it is shown that if the transmitters use the estimates ofHk,!, !k, ! " {1, 2, . . . , K},

denoted by !Hk,!, to calculate IA precoders, !fmk [t], and combiners, !wmk [t], then the mean loss in

sum rate due to imperfect CSI can be upper bounded by [14],

!Rsum #"

k,m

1

Nsclog2

#

1 +EH

$I1k,m + I2

k,m

%

"2

&

. (13)

The objective then becomes minimizing the sum leakage interference term I1k,m + I2

k,m. Using

the result from [12], the individual interference terms can be upper bounded by

NscP

dk

'''!wmk [t]

!Hk,i[t]!f !i [t]'''2=NscP

dk

'''(!wm

k [t] $!f!i [t]

)!hfk,i[t]

'''2

#NscP

dk%!wm

k [t] $ !f !i [t]%2%hfk,i[t]%

2

*

+1&

'''''hfk,i[t]

!!hfk,i[t]

%hfk,i[t]%%!h

fk,i[t]%

'''''

2,

-

=NscP

dk%!wm

k [t] $ !f !i [t]%2%hk,i[t]%2

*

+1&

'''''hk,i[t]!!hk,i[t]

%hk,i[t]%%!hk,i[t]%

'''''

2,

- ,

(14)

Power of “leakage” interference

Magnitude of quantized channel can be “normalized”

Phase of channel estimate doesn’t matter

Leakage interference has terms that look like

Required CSI lives on the Grassmann manifold

Structure can be used compress and improve feedback

Friday, November 18, 2011

Page 31: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

h!"#[t]

h!"#[t-1]

e!"#[t]

h!"#[t]

h!"$[t]

Grassmannian Differential Feedback

31

Old quantized channelNew channelTangent vector between channels

Geodesic path between channelsNew quantized channel

Tangent codebook

Channels can be related by tangent vectors and

geodesic paths

^

^

^CalculateTangent

QuantizeTangent

UpdateSystem State

Observechannel

Tangentindices

Move AlongGeodesic

Move AlongGeodesic

ReceiveQuantizedTangent

QuantizedChannel

UpdateSystem State

Differential feedback decoderDifferential feedback encoder

Quantize tangent magnitude and direction separately

Friday, November 18, 2011

Page 32: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Scalar quantization problem

Propose uniform quantization

Adaptive range to avoid error floors

Performs well in simulations

32

Initialization: Synchronous operation of the Grassman-nian differential feedback algorithm is ensured by the factthat at each iteration, both transmitter and receiver calculate aquantized channel vector based on the same commonly avail-able knowledge. For this to hold, however, both transmitterand receiver need a common initial vector, !g(0), as input tothe algorithm, otherwise the time series observed by trans-mitter and receiver will not be coupled. This vector can bebased on a memoryless quantization of the channel [4] or ini-tialized with a common random vector.Tangent Magnitude Quantization: The tangent vector

calculated in (9) is decomposed naturally into a tangent mag-nitude and a unit norm tangent direction. In this paper, themagnitude and direction are quantized separately as it canbe shown in simulation that there is little to gain from jointquantization.

The problem of quantizing the tangent magnitude is thatof quantizing a positive scalar and is done as follows

!emag = argminei!Cmag

|!e[t]! " ei| , (11)

where Cmag is the magnitude quantization codebook. The in-dex of the minimizer is then sent to the transmitter via a delayand error free link which requires log2(|Cmag|) bits. Find-ing the exact probability density function of the magnitudesis intractable, and thus we do not seek to find an optimalquantization codebook. One solution is to uniformly quan-tize a range of magnitudes, !e[t]! # [0, 1], which is subopti-mal. For example, [9] has observed from simulations that inhighly correlated channels, where such feedback strategiesare most useful, quantization error in the magnitude domi-nates the error in tangent direction.

Motivated by the correlation between magnitudes in con-secutive iterations, we propose to adapt the quantizationrange to the dynamics of the system. Given the magni-tude of a tangent at time t, !!e[t]!, the codebook at timet+1 becomes a uniform quantization codebook in the range[!!!e[t]!,min {"!!e[t]!,#/2}], where 0 < ! < 1 < " arefixed parameters of the codebook. This allows the feedbackalgorithm to accurately track the statistics of the magnitudeand quantize the current range of magnitudeswith higher res-olution. In static channels, this allows our approach to con-verge to perfect CSI.Tangent Direction Quantization: The problem of quan-

tizing the tangent direction vector is that of quantizing a unitnorm vector which lies in the tangent space orthogonal tothe base vector !g[t " 1]. General vector quantization code-books, such as a random vector codebook, can not be usedto quantize the tangent directly for several reasons. First,traditional codebooks quantize the full L dimensional spacewhereas the tangent vector is of lower dimension. Further,traditional codebooks do not enforce the structural constraintthat requires the tangent direction codewords to be orthogo-nal to the base vector !g[t]. With such a non-orthogonal tan-gent vector, the geodesic path is undefined and the output ofG(!g[t"1],!e[t], $) does not lie on the manifold. Finally, notethat the tangent space changes for each base vector, whichnecessitates an adaptive codebook.

To respect the varying tangent space geometry and or-thogonality constraints, we propose to use a canonical gen-erating codebook to be adapted at each iteration. The code-book design provided allows perfectly projecting a canonical

codebook onto the tangent plane at each iteration. This en-sures that the output of the Grassmannian differential feed-back algorithm remains on the manifold. We define a canon-ical tangent codebook as Cgen which has |Cgen| = N unitnorm vector entries x1, . . . ,xN . This vector codebook canbe any L dimensional vector codebook whose entries spanthe full L-dimensional space, such as the random codebook.At each iteration we form a codebook, Cdir, with entries or-thogonal to the base vector by using a projection operation.

Definition 1 The normalized projection matrix function

P(x,xb) =IL " xbx

"b"

1" (x"xb)2x,

computes the closest unit vector to x that is also orthogonalto the base vector xb.

Definition 2 The tangent direction codebook, Cdir(!g[t"1]),for the base point !g[t" 1] is

Cdir(!g[t" 1]) = {P(x1, !g[t" 1]), . . . , P(xN , !g[t" 1])} .

To construct good tangent codebooks, note that if the changein the channel is assumed to be isotropic, then it can be shownthat the tangent direction vector is also isotropically dis-tributed in the tangent space. This motivates finding canon-ical codebooks that lead to an isotropic distribution in thetangent space. Further improving the direction codebook de-sign, or constructing an optimal one, is left for future work.

To formalize the tangent direction quantization, recallthat the quantized channel in the next time instant will becalculated as G(!g[t " 1],!e[t], 1). Given that the loss in sumrate is related to the chordal distance between the actual andquantized channel, the quantized tangent direction will begiven as

!edir = argminxi!Cdir(!g[t#1])

d(G(!g[t" 1], !emagxi, 1),g[t]) (12)

where the tangent magnitude, !emag, is given by the output ofthe magnitude quantization step.

5. SIMULATION RESULTSIn this section we present simulation results to demonstratethe performance of IA when channel knowledge at the trans-mitter is obtained via the Grassmannian differential feedbackstrategy detailed in Section 4. To remove the limitation of aper-stream receiver, we calculate the sum rate of a decoderwhich considers all desired symbols jointly and treats leak-age interference as colored Gaussian noise. Since the fre-quency extended system can be viewed as a virtualNsc$NscMIMO system, the sum rate achieved is given by,

Rsum =K!

k=1

1Nsc

log2

"""I+#!2I+Rk

$!1

(Hk,kFkF"

kH"

kk)""" ,

where Rk =#

i$=k Hk,iFiF"iH

"k,i is the interference co-

variance matrix and the precoders, Fk =$f1k , f

2k , . . . , f

dk

k

%,

are calculated given ideal or quantized CSI. For the resultsin this section, we use the IA algorithm in [2]. Althougha closed form solution for the IA precoders exists for theSISO frequency extended interference channel in [1], it can

Initialization: Synchronous operation of the Grassman-nian differential feedback algorithm is ensured by the factthat at each iteration, both transmitter and receiver calculate aquantized channel vector based on the same commonly avail-able knowledge. For this to hold, however, both transmitterand receiver need a common initial vector, !g(0), as input tothe algorithm, otherwise the time series observed by trans-mitter and receiver will not be coupled. This vector can bebased on a memoryless quantization of the channel [4] or ini-tialized with a common random vector.Tangent Magnitude Quantization: The tangent vector

calculated in (9) is decomposed naturally into a tangent mag-nitude and a unit norm tangent direction. In this paper, themagnitude and direction are quantized separately as it canbe shown in simulation that there is little to gain from jointquantization.

The problem of quantizing the tangent magnitude is thatof quantizing a positive scalar and is done as follows

!emag = argminei!Cmag

|!e[t]! " ei| , (11)

where Cmag is the magnitude quantization codebook. The in-dex of the minimizer is then sent to the transmitter via a delayand error free link which requires log2(|Cmag|) bits. Find-ing the exact probability density function of the magnitudesis intractable, and thus we do not seek to find an optimalquantization codebook. One solution is to uniformly quan-tize a range of magnitudes, !e[t]! # [0, 1], which is subopti-mal. For example, [9] has observed from simulations that inhighly correlated channels, where such feedback strategiesare most useful, quantization error in the magnitude domi-nates the error in tangent direction.

Motivated by the correlation between magnitudes in con-secutive iterations, we propose to adapt the quantizationrange to the dynamics of the system. Given the magni-tude of a tangent at time t, !!e[t]!, the codebook at timet+1 becomes a uniform quantization codebook in the range[!!!e[t]!,min {"!!e[t]!,#/2}], where 0 < ! < 1 < " arefixed parameters of the codebook. This allows the feedbackalgorithm to accurately track the statistics of the magnitudeand quantize the current range of magnitudeswith higher res-olution. In static channels, this allows our approach to con-verge to perfect CSI.Tangent Direction Quantization: The problem of quan-

tizing the tangent direction vector is that of quantizing a unitnorm vector which lies in the tangent space orthogonal tothe base vector !g[t " 1]. General vector quantization code-books, such as a random vector codebook, can not be usedto quantize the tangent directly for several reasons. First,traditional codebooks quantize the full L dimensional spacewhereas the tangent vector is of lower dimension. Further,traditional codebooks do not enforce the structural constraintthat requires the tangent direction codewords to be orthogo-nal to the base vector !g[t]. With such a non-orthogonal tan-gent vector, the geodesic path is undefined and the output ofG(!g[t"1],!e[t], $) does not lie on the manifold. Finally, notethat the tangent space changes for each base vector, whichnecessitates an adaptive codebook.

To respect the varying tangent space geometry and or-thogonality constraints, we propose to use a canonical gen-erating codebook to be adapted at each iteration. The code-book design provided allows perfectly projecting a canonical

codebook onto the tangent plane at each iteration. This en-sures that the output of the Grassmannian differential feed-back algorithm remains on the manifold. We define a canon-ical tangent codebook as Cgen which has |Cgen| = N unitnorm vector entries x1, . . . ,xN . This vector codebook canbe any L dimensional vector codebook whose entries spanthe full L-dimensional space, such as the random codebook.At each iteration we form a codebook, Cdir, with entries or-thogonal to the base vector by using a projection operation.

Definition 1 The normalized projection matrix function

P(x,xb) =IL " xbx

"b"

1" (x"xb)2x,

computes the closest unit vector to x that is also orthogonalto the base vector xb.

Definition 2 The tangent direction codebook, Cdir(!g[t"1]),for the base point !g[t" 1] is

Cdir(!g[t" 1]) = {P(x1, !g[t" 1]), . . . , P(xN , !g[t" 1])} .

To construct good tangent codebooks, note that if the changein the channel is assumed to be isotropic, then it can be shownthat the tangent direction vector is also isotropically dis-tributed in the tangent space. This motivates finding canon-ical codebooks that lead to an isotropic distribution in thetangent space. Further improving the direction codebook de-sign, or constructing an optimal one, is left for future work.

To formalize the tangent direction quantization, recallthat the quantized channel in the next time instant will becalculated as G(!g[t " 1],!e[t], 1). Given that the loss in sumrate is related to the chordal distance between the actual andquantized channel, the quantized tangent direction will begiven as

!edir = argminxi!Cdir(!g[t#1])

d(G(!g[t" 1], !emagxi, 1),g[t]) (12)

where the tangent magnitude, !emag, is given by the output ofthe magnitude quantization step.

5. SIMULATION RESULTSIn this section we present simulation results to demonstratethe performance of IA when channel knowledge at the trans-mitter is obtained via the Grassmannian differential feedbackstrategy detailed in Section 4. To remove the limitation of aper-stream receiver, we calculate the sum rate of a decoderwhich considers all desired symbols jointly and treats leak-age interference as colored Gaussian noise. Since the fre-quency extended system can be viewed as a virtualNsc$NscMIMO system, the sum rate achieved is given by,

Rsum =K!

k=1

1Nsc

log2

"""I+#!2I+Rk

$!1

(Hk,kFkF"

kH"

kk)""" ,

where Rk =#

i$=k Hk,iFiF"iH

"k,i is the interference co-

variance matrix and the precoders, Fk =$f1k , f

2k , . . . , f

dk

k

%,

are calculated given ideal or quantized CSI. For the resultsin this section, we use the IA algorithm in [2]. Althougha closed form solution for the IA precoders exists for theSISO frequency extended interference channel in [1], it can

Initialization: Synchronous operation of the Grassman-nian differential feedback algorithm is ensured by the factthat at each iteration, both transmitter and receiver calculate aquantized channel vector based on the same commonly avail-able knowledge. For this to hold, however, both transmitterand receiver need a common initial vector, !g(0), as input tothe algorithm, otherwise the time series observed by trans-mitter and receiver will not be coupled. This vector can bebased on a memoryless quantization of the channel [4] or ini-tialized with a common random vector.Tangent Magnitude Quantization: The tangent vector

calculated in (9) is decomposed naturally into a tangent mag-nitude and a unit norm tangent direction. In this paper, themagnitude and direction are quantized separately as it canbe shown in simulation that there is little to gain from jointquantization.

The problem of quantizing the tangent magnitude is thatof quantizing a positive scalar and is done as follows

!emag = argminei!Cmag

|!e[t]! " ei| , (11)

where Cmag is the magnitude quantization codebook. The in-dex of the minimizer is then sent to the transmitter via a delayand error free link which requires log2(|Cmag|) bits. Find-ing the exact probability density function of the magnitudesis intractable, and thus we do not seek to find an optimalquantization codebook. One solution is to uniformly quan-tize a range of magnitudes, !e[t]! # [0, 1], which is subopti-mal. For example, [9] has observed from simulations that inhighly correlated channels, where such feedback strategiesare most useful, quantization error in the magnitude domi-nates the error in tangent direction.

Motivated by the correlation between magnitudes in con-secutive iterations, we propose to adapt the quantizationrange to the dynamics of the system. Given the magni-tude of a tangent at time t, !!e[t]!, the codebook at timet+1 becomes a uniform quantization codebook in the range[!!!e[t]!,min {"!!e[t]!,#/2}], where 0 < ! < 1 < " arefixed parameters of the codebook. This allows the feedbackalgorithm to accurately track the statistics of the magnitudeand quantize the current range of magnitudeswith higher res-olution. In static channels, this allows our approach to con-verge to perfect CSI.Tangent Direction Quantization: The problem of quan-

tizing the tangent direction vector is that of quantizing a unitnorm vector which lies in the tangent space orthogonal tothe base vector !g[t " 1]. General vector quantization code-books, such as a random vector codebook, can not be usedto quantize the tangent directly for several reasons. First,traditional codebooks quantize the full L dimensional spacewhereas the tangent vector is of lower dimension. Further,traditional codebooks do not enforce the structural constraintthat requires the tangent direction codewords to be orthogo-nal to the base vector !g[t]. With such a non-orthogonal tan-gent vector, the geodesic path is undefined and the output ofG(!g[t"1],!e[t], $) does not lie on the manifold. Finally, notethat the tangent space changes for each base vector, whichnecessitates an adaptive codebook.

To respect the varying tangent space geometry and or-thogonality constraints, we propose to use a canonical gen-erating codebook to be adapted at each iteration. The code-book design provided allows perfectly projecting a canonical

codebook onto the tangent plane at each iteration. This en-sures that the output of the Grassmannian differential feed-back algorithm remains on the manifold. We define a canon-ical tangent codebook as Cgen which has |Cgen| = N unitnorm vector entries x1, . . . ,xN . This vector codebook canbe any L dimensional vector codebook whose entries spanthe full L-dimensional space, such as the random codebook.At each iteration we form a codebook, Cdir, with entries or-thogonal to the base vector by using a projection operation.

Definition 1 The normalized projection matrix function

P(x,xb) =IL " xbx

"b"

1" (x"xb)2x,

computes the closest unit vector to x that is also orthogonalto the base vector xb.

Definition 2 The tangent direction codebook, Cdir(!g[t"1]),for the base point !g[t" 1] is

Cdir(!g[t" 1]) = {P(x1, !g[t" 1]), . . . , P(xN , !g[t" 1])} .

To construct good tangent codebooks, note that if the changein the channel is assumed to be isotropic, then it can be shownthat the tangent direction vector is also isotropically dis-tributed in the tangent space. This motivates finding canon-ical codebooks that lead to an isotropic distribution in thetangent space. Further improving the direction codebook de-sign, or constructing an optimal one, is left for future work.

To formalize the tangent direction quantization, recallthat the quantized channel in the next time instant will becalculated as G(!g[t " 1],!e[t], 1). Given that the loss in sumrate is related to the chordal distance between the actual andquantized channel, the quantized tangent direction will begiven as

!edir = argminxi!Cdir(!g[t#1])

d(G(!g[t" 1], !emagxi, 1),g[t]) (12)

where the tangent magnitude, !emag, is given by the output ofthe magnitude quantization step.

5. SIMULATION RESULTSIn this section we present simulation results to demonstratethe performance of IA when channel knowledge at the trans-mitter is obtained via the Grassmannian differential feedbackstrategy detailed in Section 4. To remove the limitation of aper-stream receiver, we calculate the sum rate of a decoderwhich considers all desired symbols jointly and treats leak-age interference as colored Gaussian noise. Since the fre-quency extended system can be viewed as a virtualNsc$NscMIMO system, the sum rate achieved is given by,

Rsum =K!

k=1

1Nsc

log2

"""I+#!2I+Rk

$!1

(Hk,kFkF"

kH"

kk)""" ,

where Rk =#

i$=k Hk,iFiF"iH

"k,i is the interference co-

variance matrix and the precoders, Fk =$f1k , f

2k , . . . , f

dk

k

%,

are calculated given ideal or quantized CSI. For the resultsin this section, we use the IA algorithm in [2]. Althougha closed form solution for the IA precoders exists for theSISO frequency extended interference channel in [1], it can

Initialization: Synchronous operation of the Grassman-nian differential feedback algorithm is ensured by the factthat at each iteration, both transmitter and receiver calculate aquantized channel vector based on the same commonly avail-able knowledge. For this to hold, however, both transmitterand receiver need a common initial vector, !g(0), as input tothe algorithm, otherwise the time series observed by trans-mitter and receiver will not be coupled. This vector can bebased on a memoryless quantization of the channel [4] or ini-tialized with a common random vector.Tangent Magnitude Quantization: The tangent vector

calculated in (9) is decomposed naturally into a tangent mag-nitude and a unit norm tangent direction. In this paper, themagnitude and direction are quantized separately as it canbe shown in simulation that there is little to gain from jointquantization.

The problem of quantizing the tangent magnitude is thatof quantizing a positive scalar and is done as follows

!emag = argminei!Cmag

|!e[t]! " ei| , (11)

where Cmag is the magnitude quantization codebook. The in-dex of the minimizer is then sent to the transmitter via a delayand error free link which requires log2(|Cmag|) bits. Find-ing the exact probability density function of the magnitudesis intractable, and thus we do not seek to find an optimalquantization codebook. One solution is to uniformly quan-tize a range of magnitudes, !e[t]! # [0, 1], which is subopti-mal. For example, [9] has observed from simulations that inhighly correlated channels, where such feedback strategiesare most useful, quantization error in the magnitude domi-nates the error in tangent direction.

Motivated by the correlation between magnitudes in con-secutive iterations, we propose to adapt the quantizationrange to the dynamics of the system. Given the magni-tude of a tangent at time t, !!e[t]!, the codebook at timet+1 becomes a uniform quantization codebook in the range[!!!e[t]!,min {"!!e[t]!,#/2}], where 0 < ! < 1 < " arefixed parameters of the codebook. This allows the feedbackalgorithm to accurately track the statistics of the magnitudeand quantize the current range of magnitudeswith higher res-olution. In static channels, this allows our approach to con-verge to perfect CSI.Tangent Direction Quantization: The problem of quan-

tizing the tangent direction vector is that of quantizing a unitnorm vector which lies in the tangent space orthogonal tothe base vector !g[t " 1]. General vector quantization code-books, such as a random vector codebook, can not be usedto quantize the tangent directly for several reasons. First,traditional codebooks quantize the full L dimensional spacewhereas the tangent vector is of lower dimension. Further,traditional codebooks do not enforce the structural constraintthat requires the tangent direction codewords to be orthogo-nal to the base vector !g[t]. With such a non-orthogonal tan-gent vector, the geodesic path is undefined and the output ofG(!g[t"1],!e[t], $) does not lie on the manifold. Finally, notethat the tangent space changes for each base vector, whichnecessitates an adaptive codebook.

To respect the varying tangent space geometry and or-thogonality constraints, we propose to use a canonical gen-erating codebook to be adapted at each iteration. The code-book design provided allows perfectly projecting a canonical

codebook onto the tangent plane at each iteration. This en-sures that the output of the Grassmannian differential feed-back algorithm remains on the manifold. We define a canon-ical tangent codebook as Cgen which has |Cgen| = N unitnorm vector entries x1, . . . ,xN . This vector codebook canbe any L dimensional vector codebook whose entries spanthe full L-dimensional space, such as the random codebook.At each iteration we form a codebook, Cdir, with entries or-thogonal to the base vector by using a projection operation.

Definition 1 The normalized projection matrix function

P(x,xb) =IL " xbx

"b"

1" (x"xb)2x,

computes the closest unit vector to x that is also orthogonalto the base vector xb.

Definition 2 The tangent direction codebook, Cdir(!g[t"1]),for the base point !g[t" 1] is

Cdir(!g[t" 1]) = {P(x1, !g[t" 1]), . . . , P(xN , !g[t" 1])} .

To construct good tangent codebooks, note that if the changein the channel is assumed to be isotropic, then it can be shownthat the tangent direction vector is also isotropically dis-tributed in the tangent space. This motivates finding canon-ical codebooks that lead to an isotropic distribution in thetangent space. Further improving the direction codebook de-sign, or constructing an optimal one, is left for future work.

To formalize the tangent direction quantization, recallthat the quantized channel in the next time instant will becalculated as G(!g[t " 1],!e[t], 1). Given that the loss in sumrate is related to the chordal distance between the actual andquantized channel, the quantized tangent direction will begiven as

!edir = argminxi!Cdir(!g[t#1])

d(G(!g[t" 1], !emagxi, 1),g[t]) (12)

where the tangent magnitude, !emag, is given by the output ofthe magnitude quantization step.

5. SIMULATION RESULTSIn this section we present simulation results to demonstratethe performance of IA when channel knowledge at the trans-mitter is obtained via the Grassmannian differential feedbackstrategy detailed in Section 4. To remove the limitation of aper-stream receiver, we calculate the sum rate of a decoderwhich considers all desired symbols jointly and treats leak-age interference as colored Gaussian noise. Since the fre-quency extended system can be viewed as a virtualNsc$NscMIMO system, the sum rate achieved is given by,

Rsum =K!

k=1

1Nsc

log2

"""I+#!2I+Rk

$!1

(Hk,kFkF"

kH"

kk)""" ,

where Rk =#

i$=k Hk,iFiF"iH

"k,i is the interference co-

variance matrix and the precoders, Fk =$f1k , f

2k , . . . , f

dk

k

%,

are calculated given ideal or quantized CSI. For the resultsin this section, we use the IA algorithm in [2]. Althougha closed form solution for the IA precoders exists for theSISO frequency extended interference channel in [1], it can

Magnitude codebook at time t+1

Constants chosen by user

Quantizing the Tangent Magnitude

Friday, November 18, 2011

Page 33: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Tangent direction codebook considerationsUnit norm vector quantization

Tangent space changes with channel

Orthogonality can not be violated

33

Result: Can not use constant

“Grassmannian codebooks”

Need to adapt

Need to be on tangent space

Quantizing the Tangent Direction

Friday, November 18, 2011

Page 34: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Canonical direction codebook

Projection Function

Resulting codebook:

34

Initialization: Synchronous operation of the Grassman-nian differential feedback algorithm is ensured by the factthat at each iteration, both transmitter and receiver calculate aquantized channel vector based on the same commonly avail-able knowledge. For this to hold, however, both transmitterand receiver need a common initial vector, !g(0), as input tothe algorithm, otherwise the time series observed by trans-mitter and receiver will not be coupled. This vector can bebased on a memoryless quantization of the channel [4] or ini-tialized with a common random vector.Tangent Magnitude Quantization: The tangent vector

calculated in (9) is decomposed naturally into a tangent mag-nitude and a unit norm tangent direction. In this paper, themagnitude and direction are quantized separately as it canbe shown in simulation that there is little to gain from jointquantization.

The problem of quantizing the tangent magnitude is thatof quantizing a positive scalar and is done as follows

!emag = argminei!Cmag

|!e[t]! " ei| , (11)

where Cmag is the magnitude quantization codebook. The in-dex of the minimizer is then sent to the transmitter via a delayand error free link which requires log2(|Cmag|) bits. Find-ing the exact probability density function of the magnitudesis intractable, and thus we do not seek to find an optimalquantization codebook. One solution is to uniformly quan-tize a range of magnitudes, !e[t]! # [0, 1], which is subopti-mal. For example, [9] has observed from simulations that inhighly correlated channels, where such feedback strategiesare most useful, quantization error in the magnitude domi-nates the error in tangent direction.

Motivated by the correlation between magnitudes in con-secutive iterations, we propose to adapt the quantizationrange to the dynamics of the system. Given the magni-tude of a tangent at time t, !!e[t]!, the codebook at timet+1 becomes a uniform quantization codebook in the range[!!!e[t]!,min {"!!e[t]!,#/2}], where 0 < ! < 1 < " arefixed parameters of the codebook. This allows the feedbackalgorithm to accurately track the statistics of the magnitudeand quantize the current range of magnitudeswith higher res-olution. In static channels, this allows our approach to con-verge to perfect CSI.Tangent Direction Quantization: The problem of quan-

tizing the tangent direction vector is that of quantizing a unitnorm vector which lies in the tangent space orthogonal tothe base vector !g[t " 1]. General vector quantization code-books, such as a random vector codebook, can not be usedto quantize the tangent directly for several reasons. First,traditional codebooks quantize the full L dimensional spacewhereas the tangent vector is of lower dimension. Further,traditional codebooks do not enforce the structural constraintthat requires the tangent direction codewords to be orthogo-nal to the base vector !g[t]. With such a non-orthogonal tan-gent vector, the geodesic path is undefined and the output ofG(!g[t"1],!e[t], $) does not lie on the manifold. Finally, notethat the tangent space changes for each base vector, whichnecessitates an adaptive codebook.

To respect the varying tangent space geometry and or-thogonality constraints, we propose to use a canonical gen-erating codebook to be adapted at each iteration. The code-book design provided allows perfectly projecting a canonical

codebook onto the tangent plane at each iteration. This en-sures that the output of the Grassmannian differential feed-back algorithm remains on the manifold. We define a canon-ical tangent codebook as Cgen which has |Cgen| = N unitnorm vector entries x1, . . . ,xN . This vector codebook canbe any L dimensional vector codebook whose entries spanthe full L-dimensional space, such as the random codebook.At each iteration we form a codebook, Cdir, with entries or-thogonal to the base vector by using a projection operation.

Definition 1 The normalized projection matrix function

P(x,xb) =IL " xbx

"b"

1" (x"xb)2x,

computes the closest unit vector to x that is also orthogonalto the base vector xb.

Definition 2 The tangent direction codebook, Cdir(!g[t"1]),for the base point !g[t" 1] is

Cdir(!g[t" 1]) = {P(x1, !g[t" 1]), . . . , P(xN , !g[t" 1])} .

To construct good tangent codebooks, note that if the changein the channel is assumed to be isotropic, then it can be shownthat the tangent direction vector is also isotropically dis-tributed in the tangent space. This motivates finding canon-ical codebooks that lead to an isotropic distribution in thetangent space. Further improving the direction codebook de-sign, or constructing an optimal one, is left for future work.

To formalize the tangent direction quantization, recallthat the quantized channel in the next time instant will becalculated as G(!g[t " 1],!e[t], 1). Given that the loss in sumrate is related to the chordal distance between the actual andquantized channel, the quantized tangent direction will begiven as

!edir = argminxi!Cdir(!g[t#1])

d(G(!g[t" 1], !emagxi, 1),g[t]) (12)

where the tangent magnitude, !emag, is given by the output ofthe magnitude quantization step.

5. SIMULATION RESULTSIn this section we present simulation results to demonstratethe performance of IA when channel knowledge at the trans-mitter is obtained via the Grassmannian differential feedbackstrategy detailed in Section 4. To remove the limitation of aper-stream receiver, we calculate the sum rate of a decoderwhich considers all desired symbols jointly and treats leak-age interference as colored Gaussian noise. Since the fre-quency extended system can be viewed as a virtualNsc$NscMIMO system, the sum rate achieved is given by,

Rsum =K!

k=1

1Nsc

log2

"""I+#!2I+Rk

$!1

(Hk,kFkF"

kH"

kk)""" ,

where Rk =#

i$=k Hk,iFiF"iH

"k,i is the interference co-

variance matrix and the precoders, Fk =$f1k , f

2k , . . . , f

dk

k

%,

are calculated given ideal or quantized CSI. For the resultsin this section, we use the IA algorithm in [2]. Althougha closed form solution for the IA precoders exists for theSISO frequency extended interference channel in [1], it can

Initialization: Synchronous operation of the Grassman-nian differential feedback algorithm is ensured by the factthat at each iteration, both transmitter and receiver calculate aquantized channel vector based on the same commonly avail-able knowledge. For this to hold, however, both transmitterand receiver need a common initial vector, !g(0), as input tothe algorithm, otherwise the time series observed by trans-mitter and receiver will not be coupled. This vector can bebased on a memoryless quantization of the channel [4] or ini-tialized with a common random vector.Tangent Magnitude Quantization: The tangent vector

calculated in (9) is decomposed naturally into a tangent mag-nitude and a unit norm tangent direction. In this paper, themagnitude and direction are quantized separately as it canbe shown in simulation that there is little to gain from jointquantization.

The problem of quantizing the tangent magnitude is thatof quantizing a positive scalar and is done as follows

!emag = argminei!Cmag

|!e[t]! " ei| , (11)

where Cmag is the magnitude quantization codebook. The in-dex of the minimizer is then sent to the transmitter via a delayand error free link which requires log2(|Cmag|) bits. Find-ing the exact probability density function of the magnitudesis intractable, and thus we do not seek to find an optimalquantization codebook. One solution is to uniformly quan-tize a range of magnitudes, !e[t]! # [0, 1], which is subopti-mal. For example, [9] has observed from simulations that inhighly correlated channels, where such feedback strategiesare most useful, quantization error in the magnitude domi-nates the error in tangent direction.

Motivated by the correlation between magnitudes in con-secutive iterations, we propose to adapt the quantizationrange to the dynamics of the system. Given the magni-tude of a tangent at time t, !!e[t]!, the codebook at timet+1 becomes a uniform quantization codebook in the range[!!!e[t]!,min {"!!e[t]!,#/2}], where 0 < ! < 1 < " arefixed parameters of the codebook. This allows the feedbackalgorithm to accurately track the statistics of the magnitudeand quantize the current range of magnitudeswith higher res-olution. In static channels, this allows our approach to con-verge to perfect CSI.Tangent Direction Quantization: The problem of quan-

tizing the tangent direction vector is that of quantizing a unitnorm vector which lies in the tangent space orthogonal tothe base vector !g[t " 1]. General vector quantization code-books, such as a random vector codebook, can not be usedto quantize the tangent directly for several reasons. First,traditional codebooks quantize the full L dimensional spacewhereas the tangent vector is of lower dimension. Further,traditional codebooks do not enforce the structural constraintthat requires the tangent direction codewords to be orthogo-nal to the base vector !g[t]. With such a non-orthogonal tan-gent vector, the geodesic path is undefined and the output ofG(!g[t"1],!e[t], $) does not lie on the manifold. Finally, notethat the tangent space changes for each base vector, whichnecessitates an adaptive codebook.

To respect the varying tangent space geometry and or-thogonality constraints, we propose to use a canonical gen-erating codebook to be adapted at each iteration. The code-book design provided allows perfectly projecting a canonical

codebook onto the tangent plane at each iteration. This en-sures that the output of the Grassmannian differential feed-back algorithm remains on the manifold. We define a canon-ical tangent codebook as Cgen which has |Cgen| = N unitnorm vector entries x1, . . . ,xN . This vector codebook canbe any L dimensional vector codebook whose entries spanthe full L-dimensional space, such as the random codebook.At each iteration we form a codebook, Cdir, with entries or-thogonal to the base vector by using a projection operation.

Definition 1 The normalized projection matrix function

P(x,xb) =IL " xbx

"b"

1" (x"xb)2x,

computes the closest unit vector to x that is also orthogonalto the base vector xb.

Definition 2 The tangent direction codebook, Cdir(!g[t"1]),for the base point !g[t" 1] is

Cdir(!g[t" 1]) = {P(x1, !g[t" 1]), . . . , P(xN , !g[t" 1])} .

To construct good tangent codebooks, note that if the changein the channel is assumed to be isotropic, then it can be shownthat the tangent direction vector is also isotropically dis-tributed in the tangent space. This motivates finding canon-ical codebooks that lead to an isotropic distribution in thetangent space. Further improving the direction codebook de-sign, or constructing an optimal one, is left for future work.

To formalize the tangent direction quantization, recallthat the quantized channel in the next time instant will becalculated as G(!g[t " 1],!e[t], 1). Given that the loss in sumrate is related to the chordal distance between the actual andquantized channel, the quantized tangent direction will begiven as

!edir = argminxi!Cdir(!g[t#1])

d(G(!g[t" 1], !emagxi, 1),g[t]) (12)

where the tangent magnitude, !emag, is given by the output ofthe magnitude quantization step.

5. SIMULATION RESULTSIn this section we present simulation results to demonstratethe performance of IA when channel knowledge at the trans-mitter is obtained via the Grassmannian differential feedbackstrategy detailed in Section 4. To remove the limitation of aper-stream receiver, we calculate the sum rate of a decoderwhich considers all desired symbols jointly and treats leak-age interference as colored Gaussian noise. Since the fre-quency extended system can be viewed as a virtualNsc$NscMIMO system, the sum rate achieved is given by,

Rsum =K!

k=1

1Nsc

log2

"""I+#!2I+Rk

$!1

(Hk,kFkF"

kH"

kk)""" ,

where Rk =#

i$=k Hk,iFiF"iH

"k,i is the interference co-

variance matrix and the precoders, Fk =$f1k , f

2k , . . . , f

dk

k

%,

are calculated given ideal or quantized CSI. For the resultsin this section, we use the IA algorithm in [2]. Althougha closed form solution for the IA precoders exists for theSISO frequency extended interference channel in [1], it can

Initialization: Synchronous operation of the Grassman-nian differential feedback algorithm is ensured by the factthat at each iteration, both transmitter and receiver calculate aquantized channel vector based on the same commonly avail-able knowledge. For this to hold, however, both transmitterand receiver need a common initial vector, !g(0), as input tothe algorithm, otherwise the time series observed by trans-mitter and receiver will not be coupled. This vector can bebased on a memoryless quantization of the channel [4] or ini-tialized with a common random vector.Tangent Magnitude Quantization: The tangent vector

calculated in (9) is decomposed naturally into a tangent mag-nitude and a unit norm tangent direction. In this paper, themagnitude and direction are quantized separately as it canbe shown in simulation that there is little to gain from jointquantization.

The problem of quantizing the tangent magnitude is thatof quantizing a positive scalar and is done as follows

!emag = argminei!Cmag

|!e[t]! " ei| , (11)

where Cmag is the magnitude quantization codebook. The in-dex of the minimizer is then sent to the transmitter via a delayand error free link which requires log2(|Cmag|) bits. Find-ing the exact probability density function of the magnitudesis intractable, and thus we do not seek to find an optimalquantization codebook. One solution is to uniformly quan-tize a range of magnitudes, !e[t]! # [0, 1], which is subopti-mal. For example, [9] has observed from simulations that inhighly correlated channels, where such feedback strategiesare most useful, quantization error in the magnitude domi-nates the error in tangent direction.

Motivated by the correlation between magnitudes in con-secutive iterations, we propose to adapt the quantizationrange to the dynamics of the system. Given the magni-tude of a tangent at time t, !!e[t]!, the codebook at timet+1 becomes a uniform quantization codebook in the range[!!!e[t]!,min {"!!e[t]!,#/2}], where 0 < ! < 1 < " arefixed parameters of the codebook. This allows the feedbackalgorithm to accurately track the statistics of the magnitudeand quantize the current range of magnitudeswith higher res-olution. In static channels, this allows our approach to con-verge to perfect CSI.Tangent Direction Quantization: The problem of quan-

tizing the tangent direction vector is that of quantizing a unitnorm vector which lies in the tangent space orthogonal tothe base vector !g[t " 1]. General vector quantization code-books, such as a random vector codebook, can not be usedto quantize the tangent directly for several reasons. First,traditional codebooks quantize the full L dimensional spacewhereas the tangent vector is of lower dimension. Further,traditional codebooks do not enforce the structural constraintthat requires the tangent direction codewords to be orthogo-nal to the base vector !g[t]. With such a non-orthogonal tan-gent vector, the geodesic path is undefined and the output ofG(!g[t"1],!e[t], $) does not lie on the manifold. Finally, notethat the tangent space changes for each base vector, whichnecessitates an adaptive codebook.

To respect the varying tangent space geometry and or-thogonality constraints, we propose to use a canonical gen-erating codebook to be adapted at each iteration. The code-book design provided allows perfectly projecting a canonical

codebook onto the tangent plane at each iteration. This en-sures that the output of the Grassmannian differential feed-back algorithm remains on the manifold. We define a canon-ical tangent codebook as Cgen which has |Cgen| = N unitnorm vector entries x1, . . . ,xN . This vector codebook canbe any L dimensional vector codebook whose entries spanthe full L-dimensional space, such as the random codebook.At each iteration we form a codebook, Cdir, with entries or-thogonal to the base vector by using a projection operation.

Definition 1 The normalized projection matrix function

P(x,xb) =IL " xbx

"b"

1" (x"xb)2x,

computes the closest unit vector to x that is also orthogonalto the base vector xb.

Definition 2 The tangent direction codebook, Cdir(!g[t"1]),for the base point !g[t" 1] is

Cdir(!g[t" 1]) = {P(x1, !g[t" 1]), . . . , P(xN , !g[t" 1])} .

To construct good tangent codebooks, note that if the changein the channel is assumed to be isotropic, then it can be shownthat the tangent direction vector is also isotropically dis-tributed in the tangent space. This motivates finding canon-ical codebooks that lead to an isotropic distribution in thetangent space. Further improving the direction codebook de-sign, or constructing an optimal one, is left for future work.

To formalize the tangent direction quantization, recallthat the quantized channel in the next time instant will becalculated as G(!g[t " 1],!e[t], 1). Given that the loss in sumrate is related to the chordal distance between the actual andquantized channel, the quantized tangent direction will begiven as

!edir = argminxi!Cdir(!g[t#1])

d(G(!g[t" 1], !emagxi, 1),g[t]) (12)

where the tangent magnitude, !emag, is given by the output ofthe magnitude quantization step.

5. SIMULATION RESULTSIn this section we present simulation results to demonstratethe performance of IA when channel knowledge at the trans-mitter is obtained via the Grassmannian differential feedbackstrategy detailed in Section 4. To remove the limitation of aper-stream receiver, we calculate the sum rate of a decoderwhich considers all desired symbols jointly and treats leak-age interference as colored Gaussian noise. Since the fre-quency extended system can be viewed as a virtualNsc$NscMIMO system, the sum rate achieved is given by,

Rsum =K!

k=1

1Nsc

log2

"""I+#!2I+Rk

$!1

(Hk,kFkF"

kH"

kk)""" ,

where Rk =#

i$=k Hk,iFiF"iH

"k,i is the interference co-

variance matrix and the precoders, Fk =$f1k , f

2k , . . . , f

dk

k

%,

are calculated given ideal or quantized CSI. For the resultsin this section, we use the IA algorithm in [2]. Althougha closed form solution for the IA precoders exists for theSISO frequency extended interference channel in [1], it can

The tangent space’s normal vector

Projection adapts the codebook to the channel

Can be a random codebook

Tangent Codebook Construction

Friday, November 18, 2011

Page 35: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

24

0 50 100 150 200 250 3000

0.05

0.1

0.15

0.2

0.25

0.3

0.35

Channel Realization

Chor

dal D

istan

ce

GPC (6+2) bitsRVQ with 8 bits

Proposed Algorithm (6+2) bits

Fig. 2. Chordal distance, d(!g(t), g(t)), plotted over time for a slowly channel with normalized Doppler fDTs = 0.003.This shows the high quality CSI acheived by the proposed algorithm with 7 and 3 bits for channel direction and magnitude

respectively, compared to [15] and memoryless quantization with a random codebook of 10 bits.

10!4 10!3 10!2 10!110!3

10!2

10!1

100

fd

Ts (Normalized Doppler)

Chor

dal D

istan

ce (Q

uant

izat

ion

Erro

r)

Chordal Distance vs. Bits for Magnitude

Proposed algorithm

GPC

Fig. 3. This shows the average chordal distance, or quantization error, versus the number of bits allocated for the tangent

magnitude (1 to 10 bits), while keeping direction bits fixed. This shows the good performance of the proposed algorithm when

compared with [15].

35

Previous work:Performance limited by magnitude quantization

Proposed work:No use spending more than 1 bit on magnitude

Very close to perfect magnitude knowledge

[AyaHea11]

• Simulation parameters‣ 3 tap channel‣ uniform power profile

Note: previous work in [InoHea11]

Performance Evaluation

Friday, November 18, 2011

Page 36: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

25

10!4 10!3 10!2 10!110!3

10!2

10!1

100

fD

Ts (Normalized Doppler)

Chor

dal D

istan

ce (Q

uant

izat

ion

Erro

r)

Chordal Distance vs. Direction Bits (3 to 10)

Proposed algorithm

GPC

Fig. 4. This shows the average chordal distance, or quantization error, versus the number of bits allocated for the tangent

direction, while keeping magnitude bits fixed. This shows the good performance of the proposed algorithm when compared with

[15].

0 5 10 15 20 25 30 35 400

2

4

6

8

10

12

14

16

18

SNR (dB)

Sum

Rat

e (b

its/s

/Hz)

IA w/ Perfect CSIIA with GDC f

dT

s=0.01

IA with GDC fd

Ts=0.001

IA with GDC fd

Ts=0.03

IA with GDC fd

Ts=0.05

IA with RVQRandom Beamforming

Fig. 5. This figure shows the performance of IA with imperfect CSI obtained through the proposed algorithm, as well as

random vector quantization. This shows that for slowly varying channels, the proposed algorithm allows interference alignment

networks to come very close to the perfect CSI upper bound.

36

Previous work:Error floor obvious

at 0.01 Doppler

Proposed work:No error floor

for > 3 bits

Slightly diminishing returns with random

codebooks

Benefits from more bits

• Simulation parameters‣ 3 tap channel‣ uniform power profile We can also analytically characterize performance

(see journal paper)

Performance Evaluation

Friday, November 18, 2011

Page 37: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011 37

Approaches perfect CSI

Outperforms competingpredictive algorithms

Outperforms memoryless quantization

at significant Doppler

• 3 tap channel• uniform power profile• (7+3) feedback bits• 16 channel extensions• IA via alternating min.

Proposed strategy works well even at 40dB

Codebook sizes are manageable

Performance Evaluation

Friday, November 18, 2011

Page 38: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011 38

• 3 tap channel• uniform power profile• 16 channel extensions• IA via alternating min.

IA w/ Perfect CSI

IA w/ analog feedback (example)IA with RVQ

IA w/ GDC fDTs =0.01IA w/ GDC fDTs=0.03IA w/ GDC fDTs=0.05

IA w/ GDC fDTs=0.001

0 5 10 15 20 25 30 35 402

4

6

8

10

12

14

16

18

20

SNR (dB)

Sum

Rat

e (b

/s/H

z)

Incomple

te Compar

ison

Comparison depends onFeedback powerFeedback channel capacityCoding schemeFeedback bit errors

Here we use(7+3) feedback bits

orPfb =P/2

Could do analog differential feedback...

Performance Comparison

Friday, November 18, 2011

Page 39: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Proposed Grassmannian differential feedback

Exploits temporal correlation

Reduces feedback overhead, making limited feedback more practical

Approaches perfect CSI in slow fading channels

Outperforms other feedback schemes that exploit correlation

Main limitation: training (overhead) is still incurred

Further work

Improved magnitude and direction codebooks

Analytical performance results and characterization

39

Part #2: Summary

Friday, November 18, 2011

Page 40: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

Introduction to interference alignment

Analog feedback for interference alignment

Limited feedback for interference alignment

Conclusions

40

Outline

Friday, November 18, 2011

Page 41: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

ConclusionsHighlighted the CSI requirements in IA networks

Presented methods to reduce overhead (not just feedback)

Proposed analog feedback for MIMO IA

Exploits optimality of analog transmission

Provides a low complexity feedback method

Proposed Grassmannian differential feedback for OFDM IA

Exploits CSI structure and correlation

Works well with very small codebooks

41

http://www.profheath.org/research/interference-alignment/

What are the remaining IA killers...

Friday, November 18, 2011

Page 42: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

References[CadJaf08] V. Cadambe and S. Jafar, “Interference alignment and degrees of freedom of the K-user interference channel,” IEEE Transactions on Information Theory, vol. 54, no. 8, pp. 3425–3441, August 2008.

[MadMotKha06] M. Maddah-Ali, A. Motahari, and A. Khandani, “Signaling over MIMO multi-base systems: combination of multi-access and broadcast schemes,” Proc. of IEEE International Symposium on Information Theory, pp. 2104–2108, July 2006.

[GomCadJaf08] K. Gomadam, V. Cadambe, and S. Jafar, “Approaching the capacity of wireless networks through distributed interference alignment,” Proc. of IEEE Global Telecommunications Conference, pp. 1–6, December 2008.

[PetHea09] S. W. Peters and R. W. Heath, Jr., “Interference alignment via alternating minimization,” Proc. of IEEE International Conference on Acoustics, Speech, and Signal Processing, April 2009.

[ThuBol09] I. Thukral and H. Bolcskei, “Interference alignment with limited feedback,” Proc. of IEEE International Symposium on Information Theory, pp. 1759–1763, July 2009.

42

Friday, November 18, 2011

Page 43: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

[YetJafKay09] C. M. Yetis, S. A. Jafar, and A. H. Kayran, “Feasibility conditions for interference alignment,” CoRR, vol. abs/0904.4526, 2009.

[PetHea10] S. W. Peters and R. W. Heath, Jr., “Cooperative algorithms for the MIMO interference channel,” to appear in IEEE Trans. on Sig. Proc.

[KriVar10] R. Krishnamachari, M. Varanasi, “Interference alignment under limited feedback for MIMO interference channels”, on Arxiv, 2009.

[AyaHea10] O. El Ayach, S.W. Peters, R. W. Heath, Jr., “The feasibility of Interference alignment over measured MIMO-OFDM channels”, to appear in IEEE Trans. on Veh. Tech, 2010, also available on Arxiv.

[AyaHea10-2] O. El Ayach, R.W. Heath, Jr., “Interference alignment with analog CSI feedback” IEEE Conference on Military Communication (MILCOM), Oct. 2010.

[GolPerKat09] S. Gollakota, S. Perli, D. Katabi “Interference alignment and cancellation. In Proceedings of the ACM SIGCOMM 2009 Conference on Data Communication (Barcelona, Spain, August 16 - 21, 2009). SIGCOMM '09. ACM, New York, NY,

43

References

Friday, November 18, 2011

Page 44: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

References[HoGes10] Z. Ho, D. Gesbert, “Balancing egoism and altruism on interference channel: the MIMO case” CoRR, vol. abs/0910.1688, 2010.

[Liu09] A. Liu, A. Sabharwal, Y. Liu, H. Xiang, W. Luo, “Distributed MIMO network optimization based on duality and local message passing,” In Proceeding of 47th Annual Allerton Conference, Oct 2009.

[Sri08] S. Sridharan, A. Jafarian, S. Vishwanath, S. Jafar, and S. Shamai, “A layered lattice coding scheme for a class of three user Gaussian interference channels,” in Proc. Allerton Conf. Commun. Ctrl. Cmpt., Sept. 2008, pp. 531–538.

[Sant10] I. Santamaria, O. Gonzalez, R. W. Heath, Jr., and S. W. Peters, ``Maximum Sum-Rate Interference Alignment Algorithms for MIMO Channels,'' to appear in the Proc. of IEEE Global Telecommunications Conference (GLOBECOM), Miami, FL, December 6-10, 2010.

[Luo10] Zhi-Quan Luo and Meisam Razaviyayn, “Linear Transceiver Design for Interference Alignment,” preprint.

[Pap10] D. Papailiopoulos and A.G. Dimakis, “Interference Alignment as a Rank Constrained Rank Minimization,” Proceedings of IEEE GLOBECOM, 2010.

44

Friday, November 18, 2011

Page 45: New Frontiers in Feedback for Interference Alignmentusers.ece.utexas.edu/~rheath/presentations/2011/Heath_IA_2011.pdf · O. El Ayach, S. W. Peters, and R. W. Heath, Jr., ``Feasibility

(c) Robert W. Heath Jr. 2011

References[WanJaf10] C. Wang, T. Gou and S. A. Jafar Aiming perfectly in the dark - blind interference alignment through staggered antenna switching, preprint available online at arXiv:1002.2720v1

[Jaf09] S. A. Jafar Exploiting channel correlations - simple interference alignment schemes with no CSIT, preprint available online at arXiv:0910.0555v1

[Yetis2009] Yetis, C.M.; Tiangao Gou; Jafar, S.A.; Kayran, A.H.; , "Feasibility Conditions for Interference Alignment," Global Telecommunications Conference, 2009. GLOBECOM 2009. IEEE , vol., no., pp.1-6, Nov. 30 2009-Dec. 4 2009

45

Friday, November 18, 2011