the fixed-point equation theory. To proceed, we introduce some
auxiliary variables {v
j
}
K
j=1
which represent the virtual noise
variance. Then the problem (6) with intractable per-BS power
constraints is transformed as follows
min
v
j
max
{w
m
,p
m
}
min
m
g
m
s.t.
I
m=1
v
m/I
⌈⌉
p
m
≤
K
j=1
v
j
P
j
, p
m
≥ 0, ||w
m
|| = 1
⎧
⎪
⎪
⎨
⎪
⎪
⎩
(7)
The existing results in [9–12] have revealed that the problem (7) is
dual to the following problem
min
v
j
max
{w
m
,
l
m
}
min
m
g
m
s.t.
I
m=1
l
m
≤
K
j=1
v
j
P
j
,
l
m
≥ 0, ||w
m
|| = 1.
⎧
⎪
⎪
⎨
⎪
⎪
⎩
(8)
where
l
m
denotes the transmit power of the virtual uplink user-m,
and
g
m
is the virtual uplink SINR for user-m and is given by
g
m
=
l
m
h
H
m,m
w
m
2
I
n=1,n=m
l
n
h
H
m,n
w
m
2
+ v
m/I
⌈⌉
. (9)
It is easy to see that the optimisation problem (8) is convex in {v
j
},
the problem (7) with respect to variables {v
j
} can be solved via a
subgradient projection method, and the subgradient of {v
j
}is
given by g = P
1
−
I
k=1
p
1,k
, ..., P
K
−
I
k=1
p
K,k
. Fixing {
l
n
}
and {v
j
}, the optimal receive beamformer is the MMSE receiver
which has an analytic structure, given by
w
opt
m
=
n=m
l
n
h
m,n
h
H
m,n
+ v
m/I
⌈⌉
I
M
m/I
⌈⌉
−1
h
m,m
n=m
l
n
h
m,n
h
H
m,n
+ v
m/I
⌈⌉
I
M
m/I
⌈⌉
−1
h
m,m
. (10)
According to the results obtained in [12], the uplink transmit power
with given
g
m
can be iteratively updated as follows
l
m
=
g
m
I
n=1,n=m
l
n
h
H
m,n
w
m
2
+ v
m/I
⌈⌉
h
H
m,m
w
m
2
, ∀m. (11)
Once the dual uplink problem (8) has been solved, according to the
uplink–downlink duality theory, the optimal uplink beamforming
vectors are the beamforming solution of the original downlink
problem, while the power solution of the downlink problem can
also be obtained from the uplink power allocation. Defining an
extended power vector
˜
p =
p
1
, and an extended coupling matrix
Q =
DG D1
I
1
P
max
˜
v
T
DG
1
P
max
˜
v
T
D1
I
,
⎡
⎣
⎤
⎦
(12)
where P
max
=
K
j=1
v
j
P
j
, the matrices G and D are, respectively,
given by
G
m,n
=
0 m = n
h
H
n,m
w
n
2
, m = n,
(13)
D
m,n
=
g
m
h
H
m,m
w
m
2
, m = n,0,m = n,
.
(14)
and
˜
v = v
1
···v
1
!" #
I
v
2
···v
2
!" #
I
···v
K
···v
K
!" #
I
⎛
⎝
⎞
⎠
T
. The conclusions in [9]
have revealed that the optimal power vector p is obtained as the
first
I components of the dominant eigenvector of Q, which can be
scaled so that its last component equals one.
In summary, the distributed algorithm computing the optimal
solution of the problem (7) is given in Algorithm 1 (Fig. 1) where
f
max
is the maximum eigenvector of the extended coupling matrix
of Q,
6
is the update step-size, g is the subgradient of the virtual
noise vector v, ɛ denotes the threshold of inner iteration, and δ
denotes the threshold of outer iteration. It is easy to see that the
main computational complexity in Algorithm 1 (Fig. 1) is the
matrix inversion and the singular value decomposition (SVD) in
steps 4 and 6. The complexity of the matrix inversion and the
SVD of M × M matrix is O M
2.736
and O M
3
, respectively [22].
Therefore, in Algorithm 1 (Fig. 1), the complexity of the matrix
inversion is
f
K
j=1
O M
2.736
j
and the complexity of the SVD is
w
O KI + 1()
3
, where
f
and j denote, respectively, the number
of the iterations of steps 4 and 6.
4 Algorithm for large-scale system
The beamforming vectors can be efficiently calculated via analytical
expression, however, the iterative update of both the beamformers
and the power allocation is based on the instantaneous CSI in
Algorithm 1 (Fig. 1) designed for finite-scale system. The
Fig. 1 Multicell MU-MISO max–min SINR optimisation
IET Commun., 2016, Vol. 10, Iss. 17, pp. 2380–2390
2382
&
The Institution of Engineering and Technology 2016