Lecture Notes-Multiple Antennas For MIMO Communications - Channel Correlation
Lecture Notes-Multiple Antennas For MIMO Communications - Channel Correlation
EE6832
NUS/ECE
EE6832
NUS/ECE
EE6832
ij , pq
E hij h*pq
i, p 1,2,, N
j, q 1,2,, M
(1)
NUS/ECE
EE6832
= AOD
= AOA
Receiving
array
Transmitting
array
Hon Tat Hui
Scatterers
Scatterers
NUS/ECE
EE6832
Example 1
Find the spatial correlation, 11,21, of the channels h11 and h21
of a MIMO system with N = 2 and M = 1. All the antennas
are dipole antennas. The channels are random with a
Gaussian distribution (zero mean and unit variance). Assume
that the AOA at the receiver is 360 on the plane (H-plane)
perpendicular to the dipole antennas and the radiation patterns
of the dipole antennas are omni-directional. Furthermore,
assume that the incident fields at the receiver are polarization
dr
matched.
h11
Z
V
g
o1
Vin ~
Hon Tat Hui
h21
Vo2
NUS/ECE
EE6832
Solutions
As there is only one transmitting antenna, the AOD is not
relevant for the calculation of the spatial correlation.
We define a channel as the open-circuit voltage Vo developed
at a receiving antenna to the excitation voltage Vin at a
transmitting antenna. Therefore,
Vo1
Vo 2
h11 , h21
Vin
Vin
Note that Vo1 and Vo2 are random complex numbers because
the channels h11 and h21 are random. However, Vin is
deterministic. Thus the correlation coefficient 11,21 can be
written as:
Hon Tat Hui
NUS/ECE
11,21
E h h* E h h*
E V V * E V V *
*
E h11h21
11 11
EE6832
E Vo1Vo*2
21 21
o1 o1
o2 o2
NUS/ECE
EE6832
E1, E2
Receiving dipoles
(top view)
NUS/ECE
EE6832
1
*
E Vo1Vo 2 E I z E1 dzd
I m 0 0
*
1
I I z E2 dzd
m 0 0
where I(z) is the current distribution on a dipole antennas
when it is in the transmission mode, E1() and E2() are the
incident fields on the receiving dipole antennas. Note that
E1() and E2() are random complex Gaussian numbers due
to the random nature of the channels.
2
NUS/ECE
EE6832
We can write
E V V*
o1 o 2
E1 E0
E2 E0 e jkdr cos
Costant
2
2
1
*
*
2 I z I z dzdz E E1 d E2 d
Im 0 0
0
0
1
*
2 I z I z dzdz E E1 E2* d
Im 0 0
0
1
2
*
2 I z I z dzdz E E0
Im 0 0
CJ 0 kd r
jkd r cos
E0 = path gain
from transmitter
to receiver (a
Gaussian random
number with
each scatterer)
1
J 0 kd
2
jkd r cos
e
d
10
NUS/ECE
EE6832
1
*
*
C
I
z
I
z
dzdz
E
E
2
Im 0 0
E Vo1Vo*1 E Vo 2Vo*2 C
Hence the correlation coefficient is then:
*
E Vo1Vo 2
CJ 0 kd r
11,21
J 0 kd r
CC
*
*
E V V E V V
o1 o1
o2 o2
11
NUS/ECE
EE6832
Example 2
Similar to Example 1 but now find the spatial correlation,
11,12, of the channels h11 and h12 of a MIMO system with N =
1 and M = 2, i.e., one receiving antenna and two transmitting
antennas. Assume that the AOD at the transmitter is 360.
Zg
Vin ~
dt
Zg
Vin ~
h11
Vo
h12
12
NUS/ECE
EE6832
Solutions
Now there is only one receiving antenna, the AOA is not
relevant for the calculation of the spatial correlation. The
channels are now:
Vo
Vo
h11 , h12
Vin
Vin
Thus the correlation coefficient 11,12 is:
11,12
E h h* E h h*
E V V * E V V *
*
E h11h12
11 11
E VoVo*
12 12
o o
o o
13
NUS/ECE
EE6832
Scatterers in the
far-field region of
the transmitter
(assuming in a
circular form)
plane waves
travelling to
the receiver
e1 , e2
transmitting dipoles
(top view)
1
*
E VoVo E I z ge1 dzd
I m 0 0
* e1,e2 = far fields
2
1
generated by the
I I z ge 2 dzd transmitting antennas
m 0 0
Hon Tat Hui
Multiple Antennas for MIMO Communications - Channel Correlation
14
NUS/ECE
EE6832
e1 e0 e2 e0 e jkdt cos
E VoVo
2
2
1
*
*
*
2 I z I z dzdz E ge1 d g e2 d
Im 0 0
0
0
1
2 I z I * z dzdz E ge1 g*e*2 d
Im 0 0
0
1
2
*
2 I z I z dzdz E g
Im 0 0
C J 0 kd t
jkd t cos
A constant
Hon Tat Hui
15
NUS/ECE
EE6832
Similarly,
E VoVo* E VoVo* C
Hence the correlation coefficient is then:
11,12
J kd
*
*
E V V E V V
*
E VoVo
o o
11,21
o o
16
NUS/ECE
EE6832
Example 3
Similar to Examples 1 and 2 but now find the spatial
correlation, 11,22, of the channels h11 and h22 of a MIMO
system with N = 2 and M = 2, i.e., two receiving antennas and
two transmitting antennas. Assume that the AOD at the
transmitter and AOA at the receiver are both 360.
Zg
Vin ~
dt
Vo1
dr
Zg
Vin ~
h11
h22
Vo2
17
NUS/ECE
EE6832
Solutions
Now the output voltages at the two receiving antennas Vo1
and Vo2 can be expressed in terms of the channels as:
11,22
E h h* E h h*
E V V * E V V *
*
E h11h22
11 11
*
E Vo11Vo 22
22 22
o11 o11
o 22 o 22
18
NUS/ECE
EE6832
where Vo11 and Vo22 are the partial output voltages at antenna
1 and antenna 2 that are due to signals passed through,
respectively, channels h11 and h22.
Combining the
expressions in Examples 1 and 2, we have:
E Vo11Vo*22
1 2
E I z ge1 d E1 dzd
0
Im 0 0
*
2
2
1
I I z ge 2 d E2 dzd
0
m 00
Hon Tat Hui
19
NUS/ECE
EE6832
E Vo11Vo*22
1
* z dzdz
I
z
I
2
Im 0 0
2
2
2
2
*
*
*
E ge1 d g e2 d E1 d E2 d
0
0
0
1
*
I z I z dzdz
2
Im 0 0
2
2
jkd t cos
E
g
e
d
E
E
0
C J 0 kd t J 0 kd r
jkd r cos
e
d
20
NUS/ECE
EE6832
Similarly,
*
*
E Vo11Vo11 E Vo 22Vo 22 C
Hence the correlation coefficient is then:
11,22
E Vo11Vo*22
E Vo11Vo*11 E Vo 22Vo*22
J 0 kd t J 0 kd r 11,12 11,21
21
NUS/ECE
EE6832
Notes:
In a MIMO system with arbitrary numbers of transmitting
(M) and receiving (N) dipole antennas and the antenna
separations are dt in the transmitter and dr in the receiver, the
correlation coefficients can be calculated two-by-two at a
time. The general formula is:
ij ,k J 0 kd t j J 0 kd r k i
22
NUS/ECE
EE6832
h
N 1 hN 2
Hon Tat Hui
h1M
h2 M
hNM
(2)
23
NUS/ECE
EE6832
hN 1
h
12
vec(H )
(the dimension of vec(H ) is NM ? )
hN 2
(3)
h1M
hNM
Hon Tat Hui
24
NUS/ECE
EE6832
(4)
(5)
25
NUS/ECE
EE6832
26
NUS/ECE
EE6832
Example 4
Write a Matlab program to obtain the channel matrix of a 33
MIMO system equipped with dipole antennas aligned as
uniform linear arrays (ULAs). The antenna separations at the
transmitter and receiver are 0.2 and 0.15, respectively. The
AOD at the transmitter and the AOA at the receiver of the
multipath signals are all 360. Assume that the channels are
Gaussian random channels with a unit variance and a zero
mean, and the antenna mutual coupling can be ignored.
Calculate the channel capacity when the SNR = 20dB.
Solutions
dt = 0.2, dr = 0.15
Hon Tat Hui
27
NUS/ECE
EE6832
h11 h12
H h21 h22
h31 h32
h13
h23 , hij CN 0,1
h33
28
NUS/ECE
*
E
h
h
*
E
h
h
E h h*
11 11
r E h21h11*
E h h*
31 11
J 0 0.3
J 0 0.6
*
E h11h21
21 21
31 21
J 0 0.3
1
J 0 0.3
EE6832
*
E h11h31
*
E h21h31
*
E h31h31
J 0 0.6
J 0 0.3
29
NUS/ECE
EE6832
*
E
h
h
E h h*
E h h*
11 11
*
t E h12 h11
E h h*
13 11
*
E h11h12
12 12
13 12
*
E h11h13
1
J 0 0.4 J 0 0.8
* J 0.4
1
E h12 h13
J 0 0.4
0
J 0.8 J 0.4
1
0
0
E h13h13
30
NUS/ECE
EE6832
31
NUS/ECE
EE6832
for i=1:M;
for j=1:M;
pt(i,j)=bessel(0,k*dt*abs(j-i));
end;
end;
RH=kron(pt,pr);
[V,D] = eig(RH);
G=V*sqrt(D);
%-----------channel matrix generation
snrdB=20;
snr=10^(snrdB/10);
for n=1:5000;
r=sqrt(0.5)*(randn(N,M)+1j*randn(N,M));
32
NUS/ECE
EE6832
for j=1:M;
for i=1:N;
vec_r(i+(j-1)*N)=r(i,j);
end;
end;
vec_H=G*vec_r';
for j=1:M;
for i=1:N;
H(i,j)=vec_H(i+(j-1)*N);
end;
end;
%-----------capacity calculation
C(n)=log2(real(det(eye(N)+snr/M*(H'*H))));
end;
cdfplot(C)
Average_C=mean(C)
Hon Tat Hui
33
NUS/ECE
EE6832
c df(C)
0.6
0.5
0.4
0.3
0.2
0.1
0
4
10
12
14
C (bits /s /Hz)
16
18
20
34
NUS/ECE
EE6832
35
NUS/ECE
EE6832
Z g 2 Z 22
Z gM Z MM
Z 21
Z2M
1
Z Z
(7)
Z
Z
Zt g1
11
gM
MM
ZM 2
ZM1
Z g1 Z11 Z g 2 Z 22
36
NUS/ECE
EE6832
v c Z r 1v u
(8)
where Zr is the mutual impedance matrix containing the
receiving mutual impedances (see Lecture Notes on
Mutual Coupling in Antenna Arrays):
Z t21
Zr Z L
N1
Z
t
Z L
Hon Tat Hui
Z t12
ZL
1
Z tN 2
ZL
Z t1N
ZL
Z t2 N
ZL
(9)
37
NUS/ECE
EE6832
ZL
vu
v oc
Zin Z L
(10)
ZL
vc
Z r 1v oc
Zin Z L
Hon Tat Hui
(11)
38
NUS/ECE
EE6832
Zin Z L
voc
vc
ZL
(12)
voc Z r1HZt1v s v n
Hon Tat Hui
(13)
39
NUS/ECE
EE6832
40
NUS/ECE
EE6832
Example 5
Re-do Example 4 but now take the antenna mutual coupling
into account. It is given that the mutual impedance between
two transmitting antennas are:
dt = 0.2, Z12 = 25.91-j15.34 , Z21 = 25.28-j15.78
dt = 0.4, Z12 = -0.90-j20.30 , Z21 = -1.42-j20.11
The mutual impedance between two receiving antennas are:
dr = 0.15, Zt12 = -17.73+j2.75 , Zt21 = -17.48+j2.94
dr = 0.30, Zt12 = -8.29+j10.44 , Zt21 = -7.96+j10.51
The internal impedance of the dipole antennas is:
Zin = 39.00+j7.17 = Z11 = Z22 = = ZMM
The terminal load impedance of the dipole antennas is:
ZL = 50
Hon Tat Hui
41
NUS/ECE
EE6832
Solutions
N=M=3
dt = 0.2, dr = 0.15
1
0.28-j0.19 -0.03-j0.23
1
0.28-j0.19
Z t 0.27-j0.20
1
-0.03-j0.22 0.27-j0.20
Z r 0.35-j0.06
0.16-j0.21
Hon Tat Hui
0.35-j0.05 0.17-j0.21
1
0.35-j0.06
0.35-j0.05
1
42
NUS/ECE
EE6832
are
shown
below
(filename:
clear all;
M=3; % number of transmit antennas
N=3; % number of receive antennas
k=2*pi;
dr=0.15 %lambda
dt=0.2 %lambda
%-----------spatial channel correlations generation
for i=1:N;
for j=1:N;
pr(i,j)=bessel(0,k*dr*abs(j-i));
end;
end;
Hon Tat Hui
43
NUS/ECE
EE6832
for i=1:M;
for j=1:M;
pt(i,j)=bessel(0,k*dt*abs(j-i));
end;
end;
RH=kron(pt,pr);
[V,D] = eig(RH);
G=V*sqrt(D);
44
NUS/ECE
EE6832
z12=25.9059+1j*(-15.3365);
z21=25.2796+1j*(-15.7831);
z13=-0.8920+1j*(-20.3036);
z31=-1.4192+1j*(-20.1113);
zt12=-17.73449488+1j*(2.74569212);
zt21=-17.47727875+1j*(2.94131405);
zt13=-8.28960286+1j*(10.43902986);
zt31=-7.96114038+1j*(10.50848904);
zt=[
1
z21/(zl+zin)
z31/(zl+zin)
]
z12/(zl+zin)
1
z12/(zl+zin)
z13/(zl+zin)
z21/(zl+zin)
1
45
NUS/ECE
EE6832
zr=[
1
-zt21/zl
-zt31/zl
]
-zt12/zl
1
-zt12/zl
-zt13/zl
-zt21/zl
1
46
NUS/ECE
EE6832
vec_H=G*vec_r';
for j=1:M;
for i=1:N;
H(i,j)=vec_H(i+(j-1)*N);
end;
end;
H1=inv(zr)*H*inv(zt);
%-----------capacity calculation
C(n)=log2(real(det(eye(N)+snr/M*(H1*H1'))));
end;
cdfplot(C)
Average_C=mean(C)
47
NUS/ECE
EE6832
cdf (bits/s/Hz)
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
4
10
12
14
C (bits/s/Hz)
16
18
20
48
NUS/ECE
EE6832
Final remarks:
Computer modeling and simulation of channels plays an
important role in MIMO system design. Except the channel
modulation method as discussed in Examples 4 and 5 (in
which the scatterer effect and the antenna mutual coupling
effect are modeled separately), we can also adopt a holistic
modeling approach in which the scatterer effect and the
antenna mutual coupling effect are built into the channel
matrix H using a single-step modeling method, without the
need for separate calculations of the spatial correlation and
antenna mutual coupling. This holistic approach combines
antenna EM simulation with random channel modeling.
Interested students can read reference [4].
Hon Tat Hui
49
NUS/ECE
EE6832
References:
[1] J. W. Wallace and M. A. Jensen, Modeling the indoor MIMO wireless
channel, IEEE Transactions on Antennas and Propagation, vol. 50, no. 5, pp.
591-599, 2002.
[2] R. Janaswamy, Effect of element mutual coupling on the capacity of fixed
length linear arrays, IEEE Antennas and Wireless Propagation Letters, vol. 1,
pp. 157-160, 2002.
[3] H. T. Hui, "Influence of antenna characteristics on MIMO systems with
compact monopole arrays," IEEE Antennas and Wireless Propagation Letters,
vol. 8, pp. 133-136, 2009.
[4] H. T. Hui and Xuan Wang, Building antenna characteristics into MIMO
channel simulation, International Journal of Electronics, vol. 97, no. 6, pp.
703-714, 2010.
50