Effect of Antenna Correlation On Massive MIMO (2) : Rohit Budhiraja MIMO Wireless Communications (EE677)
Effect of Antenna Correlation On Massive MIMO (2) : Rohit Budhiraja MIMO Wireless Communications (EE677)
Effect of Antenna Correlation On Massive MIMO (2) : Rohit Budhiraja MIMO Wireless Communications (EE677)
Rohit Budhiraja
Nov 5, 2019
For uncorrelated fading R = IM . Correlation matrices are normalized such that Tr(R) = M
For uncorrelated fading R = IM . Correlation matrices are normalized such that Tr(R) = M
Spatial channel correlation makes 30 of the 100 eigenvalues larger than in uncorrelated case
For uncorrelated fading R = IM . Correlation matrices are normalized such that Tr(R) = M
Spatial channel correlation makes 30 of the 100 eigenvalues larger than in uncorrelated case
Spatial channel correlation reduces the remaining eigenvalues to substantially smaller values
For uncorrelated fading R = IM . Correlation matrices are normalized such that Tr(R) = M
Spatial channel correlation makes 30 of the 100 eigenvalues larger than in uncorrelated case
Spatial channel correlation reduces the remaining eigenvalues to substantially smaller values
We can generate any random channel vector h ∼ NC (0M , R)
For uncorrelated fading R = IM . Correlation matrices are normalized such that Tr(R) = M
Spatial channel correlation makes 30 of the 100 eigenvalues larger than in uncorrelated case
Spatial channel correlation reduces the remaining eigenvalues to substantially smaller values
We can generate any random channel vector h ∼ NC (0M , R)
1 1
perfom eigen value decomposition of R = UDUH ; take its square root R 2 = UD 2 UH
For uncorrelated fading R = IM . Correlation matrices are normalized such that Tr(R) = M
Spatial channel correlation makes 30 of the 100 eigenvalues larger than in uncorrelated case
Spatial channel correlation reduces the remaining eigenvalues to substantially smaller values
We can generate any random channel vector h ∼ NC (0M , R)
1 1
perfom eigen value decomposition of R = UDUH ; take its square root R 2 = UD 2 UH
1
generate ě ∼ NC (0M , IM ) and perform h = R 2 ě
For uncorrelated fading R = IM . Correlation matrices are normalized such that Tr(R) = M
Spatial channel correlation makes 30 of the 100 eigenvalues larger than in uncorrelated case
Spatial channel correlation reduces the remaining eigenvalues to substantially smaller values
We can generate any random channel vector h ∼ NC (0M , R)
1 1
perfom eigen value decomposition of R = UDUH ; take its square root R 2 = UD 2 UH
1
generate ě ∼ NC (0M , IM ) and perform h = R 2 ě
H
note that E{hh } = R
MIMO Wireless Communications (Rohit Budhiraja, IITK) Correlation 2
Basic impact of spatial channel correlation (1)
Consider uplink of a single-cell with K users. Assume hk ∼ (0M , Rk ), k = 1, . . . , K
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
factor K normalizes the average channel gain such that βk = M1 Tr(Rk ) = 1
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
factor K normalizes the average channel gain such that βk = M1 Tr(Rk ) = 1
note that Eke k k2 = M/K (will later use this fact)
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
factor K normalizes the average channel gain such that βk = M1 Tr(Rk ) = 1
note that Eke k k2 = M/K (will later use this fact)
Channel model implies that each UE has a strongly spatially correlated channel
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
factor K normalizes the average channel gain such that βk = M1 Tr(Rk ) = 1
note that Eke k k2 = M/K (will later use this fact)
Channel model implies that each UE has a strongly spatially correlated channel
with only M/K rather than M non-zero eigenvalues of R
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
factor K normalizes the average channel gain such that βk = M1 Tr(Rk ) = 1
note that Eke k k2 = M/K (will later use this fact)
Channel model implies that each UE has a strongly spatially correlated channel
with only M/K rather than M non-zero eigenvalues of R
Eigenspaces of individual correlation matrices are all orthogonal.
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
factor K normalizes the average channel gain such that βk = M1 Tr(Rk ) = 1
note that Eke k k2 = M/K (will later use this fact)
Channel model implies that each UE has a strongly spatially correlated channel
with only M/K rather than M non-zero eigenvalues of R
Eigenspaces of individual correlation matrices are all orthogonal.
implies although the UEs’ channels are random, they “live” in mutually orthogonal subspaces
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
factor K normalizes the average channel gain such that βk = M1 Tr(Rk ) = 1
note that Eke k k2 = M/K (will later use this fact)
Channel model implies that each UE has a strongly spatially correlated channel
with only M/K rather than M non-zero eigenvalues of R
Eigenspaces of individual correlation matrices are all orthogonal.
implies although the UEs’ channels are random, they “live” in mutually orthogonal subspaces
K
X
y = hi si + n
i =1
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
factor K normalizes the average channel gain such that βk = M1 Tr(Rk ) = 1
note that Eke k k2 = M/K (will later use this fact)
Channel model implies that each UE has a strongly spatially correlated channel
with only M/K rather than M non-zero eigenvalues of R
Eigenspaces of individual correlation matrices are all orthogonal.
implies although the UEs’ channels are random, they “live” in mutually orthogonal subspaces
K
X
y = hi si + n
i =1
K
!
X
UHk y = UHk hi si + n
i =1
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
factor K normalizes the average channel gain such that βk = M1 Tr(Rk ) = 1
note that Eke k k2 = M/K (will later use this fact)
Channel model implies that each UE has a strongly spatially correlated channel
with only M/K rather than M non-zero eigenvalues of R
Eigenspaces of individual correlation matrices are all orthogonal.
implies although the UEs’ channels are random, they “live” in mutually orthogonal subspaces
K
X
y = hi si + n
i =1
K
! K
X X √
UHk y = UHk hi si + n = K UHk Ui e i si + UHk n
i =1 i =1
here Uk ∈ CM×M/K are tall unitary matrices (UHk Uk = IM/K ) and UHk Uj = 0M , ∀k 6= j
factor K normalizes the average channel gain such that βk = M1 Tr(Rk ) = 1
note that Eke k k2 = M/K (will later use this fact)
Channel model implies that each UE has a strongly spatially correlated channel
with only M/K rather than M non-zero eigenvalues of R
Eigenspaces of individual correlation matrices are all orthogonal.
implies although the UEs’ channels are random, they “live” in mutually orthogonal subspaces
K
X
y = hi si + n
i =1
K
! K
X X √ √
UHk y = UHk hi si + n = K UHk Ui e i si + UHk n = K e k sk + ň k
i =1 i =1
Kpk ke k k2
E{SINRk } = E 2
σUL
Signal power of each UE gets a boost of M – antenna array captures the same amount of energy
Signal power of each UE gets a boost of M – antenna array captures the same amount of energy
but this energy is concentrated on a subset of the eigen/spatial directions
Signal power of each UE gets a boost of M – antenna array captures the same amount of energy
but this energy is concentrated on a subset of the eigen/spatial directions
We compare the above case with the one where all UEs have same R
Signal power of each UE gets a boost of M – antenna array captures the same amount of energy
but this energy is concentrated on a subset of the eigen/spatial directions
We compare the above case with the one where all UEs have same R
R = K UUH
Signal power of each UE gets a boost of M – antenna array captures the same amount of energy
but this energy is concentrated on a subset of the eigen/spatial directions
We compare the above case with the one where all UEs have same R
√
R = K UUH ⇒ hk = K Ue k ,
Signal power of each UE gets a boost of M – antenna array captures the same amount of energy
but this energy is concentrated on a subset of the eigen/spatial directions
We compare the above case with the one where all UEs have same R
√
R = K UUH ⇒ hk = K Ue k , where e k ∼ NC 0M/K , IM/K
Signal power of each UE gets a boost of M – antenna array captures the same amount of energy
but this energy is concentrated on a subset of the eigen/spatial directions
We compare the above case with the one where all UEs have same R
√
R = K UUH ⇒ hk = K Ue k , where e k ∼ NC 0M/K , IM/K
We have
K
X
y = h i si + n
i=1
We have
K
X
y = h i si + n
i=1
K
!
X
H H
U y = U h i si + n
i=1
We have
K
X
y = h i si + n
i=1
K
! K √
X X
H H
U y = U h i si + n = K UH Ue i si + UH
kn
i=1 i=1
We have
K
X
y = h i si + n
i=1
K
! K √ K
H H
X X √ X
U y = U h i si + n = K UH Ue i si + UH
k n = K e i si + ň
i=1 i=1 i=1
We have
K
X
y = h i si + n
i=1
K
! K √ K
H H
X X √ X
U y = U h i si + n = K UH Ue i si + UH
k n = K e i si + ň
i=1 i=1 i=1
We have
K
X
y = h i si + n
i=1
K
! K √ K
H H
X X √ X
U y = U h i si + n = K UH Ue i si + UH
k n = K e i si + ň
i=1 i=1 i=1
We have
K
X
y = h i si + n
i=1
K
! K √ K
H H
X X √ X
U y = U h i si + n = K UH Ue i si + UH
k n = K e i si + ň
i=1 i=1 i=1
We have
K
X
y = h i si + n
i=1
K
! K √ K
H H
X X √ X
U y = U h i si + n = K UH Ue i si + UH
k n = K e i si + ň
i=1 i=1 i=1
We have
K
X
y = h i si + n
i=1
K
! K √ K
H H
X X √ X
U y = U h i si + n = K UH Ue i si + UH
k n = K e i si + ň
i=1 i=1 i=1
Received signal at BS is sum of Npath multipath components, where Npath is a large number.
Received signal at BS is sum of Npath multipath components, where Npath is a large number.
Suppose scattering is localized around UE. BS is elevated and has no scatterers in its near-field
Received signal at BS is sum of Npath multipath components, where Npath is a large number.
Suppose scattering is localized around UE. BS is elevated and has no scatterers in its near-field
ϕ̄ = ϕ + δ, where ϕ is a deterministic nominal angle and
Received signal at BS is sum of Npath multipath components, where Npath is a large number.
Suppose scattering is localized around UE. BS is elevated and has no scatterers in its near-field
ϕ̄ = ϕ + δ, where ϕ is a deterministic nominal angle and
δ is a random deviation from the nominal angle with standard deviation σϕ
Received signal at BS is sum of Npath multipath components, where Npath is a large number.
Suppose scattering is localized around UE. BS is elevated and has no scatterers in its near-field
ϕ̄ = ϕ + δ, where ϕ is a deterministic nominal angle and
δ is a random deviation from the nominal angle with standard deviation σϕ
2
Gaussian distributed δ ∼ N (0, σϕ ),
Received signal at BS is sum of Npath multipath components, where Npath is a large number.
Suppose scattering is localized around UE. BS is elevated and has no scatterers in its near-field
ϕ̄ = ϕ + δ, where ϕ is a deterministic nominal angle and
δ is a random deviation from the nominal angle with standard deviation√σϕ
2
Gaussian distributed δ ∼ N (0, σϕ ),Laplacian distributed δ ∼ Lap(0, σϕ / 2),
Received signal at BS is sum of Npath multipath components, where Npath is a large number.
Suppose scattering is localized around UE. BS is elevated and has no scatterers in its near-field
ϕ̄ = ϕ + δ, where ϕ is a deterministic nominal angle and
δ is a random deviation from the nominal angle with standard deviation√σϕ
2
Gaussian distributed δ ∼ N (0, σϕ ),Laplacian distributed δ ∼ Lap(0, σϕ / 2), Uniformily distributed
√ √
deviations δ ∼ U[− 3σϕ , 3σϕ ]
Received signal at BS is sum of Npath multipath components, where Npath is a large number.
Suppose scattering is localized around UE. BS is elevated and has no scatterers in its near-field
ϕ̄ = ϕ + δ, where ϕ is a deterministic nominal angle and
δ is a random deviation from the nominal angle with standard deviation√σϕ
2
Gaussian distributed δ ∼ N (0, σϕ ),Laplacian distributed δ ∼ Lap(0, σϕ / 2), Uniformily distributed
√ √
deviations δ ∼ U[− 3σϕ , 3σϕ ]
Standard deviation σϕ is measured in radians and is called the angular standard deviation (ASD)
Received signal at BS is sum of Npath multipath components, where Npath is a large number.
Suppose scattering is localized around UE. BS is elevated and has no scatterers in its near-field
ϕ̄ = ϕ + δ, where ϕ is a deterministic nominal angle and
δ is a random deviation from the nominal angle with standard deviation√σϕ
2
Gaussian distributed δ ∼ N (0, σϕ ),Laplacian distributed δ ∼ Lap(0, σϕ / 2), Uniformily distributed
√ √
deviations δ ∼ U[− 3σϕ , 3σϕ ]
Standard deviation σϕ is measured in radians and is called the angular standard deviation (ASD)
since it determines how large the deviations from the nominal angle are
Received signal at BS is sum of Npath multipath components, where Npath is a large number.
Suppose scattering is localized around UE. BS is elevated and has no scatterers in its near-field
ϕ̄ = ϕ + δ, where ϕ is a deterministic nominal angle and
δ is a random deviation from the nominal angle with standard deviation√σϕ
2
Gaussian distributed δ ∼ N (0, σϕ ),Laplacian distributed δ ∼ Lap(0, σϕ / 2), Uniformily distributed
√ √
deviations δ ∼ U[− 3σϕ , 3σϕ ]
Standard deviation σϕ is measured in radians and is called the angular standard deviation (ASD)
since it determines how large the deviations from the nominal angle are
A reasonable value of σϕ in urban cellular networks is 10o
MIMO Wireless Communications (Rohit Budhiraja, IITK) Correlation 6
Local Scattering Spatial Correlation Model (2)
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
h iT
Array response of each path an = gn 1 e 2πj dH sin(ϕ̄n ) · · · e 2πj dH (M−1) sin(ϕ̄n )
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
h iT
Array response of each path an = gn 1 e 2πj dH sin(ϕ̄n ) · · · e 2πj dH (M−1) sin(ϕ̄n )
gn accounts for gain of nth path. dH is antenna spacing
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
h iT
Array response of each path an = gn 1 e 2πj dH sin(ϕ̄n ) · · · e 2πj dH (M−1) sin(ϕ̄n )
gn accounts for gain of nth path. dH is antenna spacing
gn are i.i.d. random variables with zero-mean and variance E{|gn |2 }.
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
h iT
Array response of each path an = gn 1 e 2πj dH sin(ϕ̄n ) · · · e 2πj dH (M−1) sin(ϕ̄n )
gn accounts for gain of nth path. dH is antenna spacing
gn are i.i.d. random variables with zero-mean and variance E{|gn |2 }.
variance represents the average gain of the nth path
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
h iT
Array response of each path an = gn 1 e 2πj dH sin(ϕ̄n ) · · · e 2πj dH (M−1) sin(ϕ̄n )
gn accounts for gain of nth path. dH is antenna spacing
gn are i.i.d. random variables with zero-mean and variance E{|gn |2 }.
variance represents the average gain of the nth path
PNpath
total average gain of the multipath components is β = n=1
E{|gn |2 }
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
h iT
Array response of each path an = gn 1 e 2πj dH sin(ϕ̄n ) · · · e 2πj dH (M−1) sin(ϕ̄n )
gn accounts for gain of nth path. dH is antenna spacing
gn are i.i.d. random variables with zero-mean and variance E{|gn |2 }.
variance represents the average gain of the nth path
PNpath
total average gain of the multipath components is β = n=1
E{|gn |2 }
assume angles ϕ̄n are i.i.d. random variables with angular pdf f (ϕ̄)
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
h iT
Array response of each path an = gn 1 e 2πj dH sin(ϕ̄n ) · · · e 2πj dH (M−1) sin(ϕ̄n )
gn accounts for gain of nth path. dH is antenna spacing
gn are i.i.d. random variables with zero-mean and variance E{|gn |2 }.
variance represents the average gain of the nth path
PNpath
total average gain of the multipath components is β = n=1
E{|gn |2 }
assume angles ϕ̄n are i.i.d. random variables with angular pdf f (ϕ̄)
Channel response Npath
X
h= an with
n=1
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
h iT
Array response of each path an = gn 1 e 2πj dH sin(ϕ̄n ) · · · e 2πj dH (M−1) sin(ϕ̄n )
gn accounts for gain of nth path. dH is antenna spacing
gn are i.i.d. random variables with zero-mean and variance E{|gn |2 }.
variance represents the average gain of the nth path
PNpath
total average gain of the multipath components is β = n=1
E{|gn |2 }
assume angles ϕ̄n are i.i.d. random variables with angular pdf f (ϕ̄)
Channel response Npath
X
an with R = E hhH
h=
n=1
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
h iT
Array response of each path an = gn 1 e 2πj dH sin(ϕ̄n ) · · · e 2πj dH (M−1) sin(ϕ̄n )
gn accounts for gain of nth path. dH is antenna spacing
gn are i.i.d. random variables with zero-mean and variance E{|gn |2 }.
variance represents the average gain of the nth path
PNpath
total average gain of the multipath components is β = n=1
E{|gn |2 }
assume angles ϕ̄n are i.i.d. random variables with angular pdf f (ϕ̄)
Channel response
Npath NX
path
X
H
an aH
h= an with R = E hh =E n
n=1 n=1
Each multipath components results in a plane wave that reaches array from a particular angle ϕ̄n
h iT
Array response of each path an = gn 1 e 2πj dH sin(ϕ̄n ) · · · e 2πj dH (M−1) sin(ϕ̄n )
gn accounts for gain of nth path. dH is antenna spacing
gn are i.i.d. random variables with zero-mean and variance E{|gn |2 }.
variance represents the average gain of the nth path
PNpath
total average gain of the multipath components is β = n=1
E{|gn |2 }
assume angles ϕ̄n are i.i.d. random variables with angular pdf f (ϕ̄)
Channel response
Npath NX
path
X
H
an aH
h= an with R = E hh =E n
n=1 n=1