矩阵单调性优化框架
总论
本文是北理工邢成文老师于2015年的一份工作,主要讨论通过函数的单调性来达到降低优化问题维度,拿到最优解形式的目的。
参考资料:
C. Xing, S. Ma and Y. Zhou, “Matrix-Monotonic Optimization for MIMO Systems,” in IEEE Transactions on Signal Processing, vol. 63, no. 2, pp. 334-348, Jan.15, 2015, doi:10.1109/TSP.2014.2373332.
文章结构
首先在绪论中,作者介绍了针对矩阵优化问题的一般解决思路,并引出了达到这个思路的两种方法:基于KKT条件的优化方法以及基于最优化理论的优化方法。在分析了这两种方法的本质和优缺点后,作者提出了一个思考:**是否可以将函数的单调性像凸性一样变成解优化问题的一个切入点呢?**本文既是在这个思考下做出的一些尝试。
文章结构如下:
在II节用一系列MIMO例子引出MVOP(matrix-variate optimization problems)问题的普适性;
在III节介绍MMOP(matrix-monotonic optimization problems)优化框架;
在IV节介绍MVOP如何转化为MMOP问题(基于酉阵的转化方法);
在V节扩展到多变量问题。
下面按照文章内容进行介绍。
I.绪论
绪论主要有3个问题:
- 矩阵优化问题的挑战及解决思路,并介绍了两种实现方法并分析了本质及缺陷;
- 本文的思考出发点
- 本文相比其他工作的不同
- 贡献
矩阵优化问题
解决思路
矩阵优化问题的挑战主要是优化变量维度过高,造成很大的复杂度。所以希望能降低优化问题的维度,而实现的方式是利用最优解的结构信息来降低问题维度。而最优解结构信息的获得主流有两种方式:基于KKT条件的优化方法以及基于有优化理论的优化方法。
基于KKT条件的优化方法
KKT条件提供了最优解的必要条件,当优化问题满足强对偶条件时则进化为充要条件,所以满足KKT条件的最优解的共同结构则是我们所需要的最优解结构。
KKT条件方法的主要问题是当问题过于复杂,我们可能无法通过该方法获得最优结构。
基于最优化理论的优化方法
该方法的思路是将原始的性能指标换为一个MSE矩阵对角线单元的函数。
而该思路是system-dependent的,并且要求MSE函数必须是关于优化变量的函数。
本文思考出发点
单调性和凸性一样是函数的基本性质,所以希望能从单调性出发给出一个优化矩阵函数的优化方法。
本文相比其他工作的不同
2点不同:
- 目标函数是一个H阵,而不是H阵的标量函数(从这里看貌似适用范围变小了?)
- 最终效果是给出了一个最优闭式解,而不是给出了最优解的数值计算方法(这是一个进步)
这两点不同的比较对象见原论文,具体在第二页右上角。
贡献
- 给出了MMOP优化框架;
- 在上述框架基础上讨论了Pareto最优解的实现(即在该最有结构基础上,通过注水法可得到资源的最佳分配,Pareto最优解是资源分配问题,原问题指在系统中让无人受损基础上有人收益更大,稳定之后即资源最优分配);
- MVOP如何转化为MMOP(基于酉阵的转化方法);
- 扩展到多变量。
II.问题统一
本小节作者提出了一系列MIMO中的基本问题,并转化为MVOP问题,旨在说明MVOP问题应用的广泛性。
MIMO中问题
互信息最大
Case 1:
min
X
−
log
∣
X
H
H
H
K
X
−
1
H
X
+
N
∣
s.t.
K
X
=
Tr
(
X
X
H
Ψ
)
Σ
+
σ
n
2
I
Tr
(
X
X
H
)
≤
P
(1)
\begin{aligned} \text { Case 1: } \min _{\mathbf{X}} &-\log \left|\mathbf{X}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{X}}^{-1} \mathbf{H X}+\mathbf{N}\right| \\ \text { s.t. } & \mathbf{K}_{\mathbf{X}}=\operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}}\right) \leq P \end{aligned} \tag{1}
Case 1: Xmin s.t. −log∣∣XHHHKX−1HX+N∣∣KX=Tr(XXHΨ)Σ+σn2ITr(XXH)≤P(1)
其中
H
,
N
,
Ψ
,
Σ
\mathbf{H},\mathbf{N},\boldsymbol{\Psi},\boldsymbol{\Sigma}
H,N,Ψ,Σ为恒定矩阵,
N
,
Ψ
,
Σ
\mathbf{N},\boldsymbol{\Psi},\boldsymbol{\Sigma}
N,Ψ,Σ为半正定矩阵。同样问题的不同表述为:
Case 2:
min
X
−
log
∣
A
H
X
H
H
H
K
X
−
1
H
X
A
+
I
∣
s.t.
K
X
=
Tr
(
X
X
H
Ψ
)
Σ
+
σ
n
2
I
Tr
(
X
X
H
)
≤
P
,
(2)
\begin{aligned} \text { Case 2: } \min _{\mathbf{X}} &-\log \left|\mathbf{A}^{\mathrm{H}} \mathbf{X}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{X}}^{-1} \mathbf{H} \mathbf{X} \mathbf{A}+\mathbf{I}\right| \\ \text { s.t. } & \mathbf{K}_{\mathbf{X}}=\operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}}\right) \leq P, \end{aligned} \tag{2}
Case 2: Xmin s.t. −log∣∣AHXHHHKX−1HXA+I∣∣KX=Tr(XXHΨ)Σ+σn2ITr(XXH)≤P,(2)
A
,
N
\mathbf{A},\mathbf{N}
A,N都满秩。
MSE最小化问题
Case 3:
min
X
Tr
[
(
X
H
H
H
K
X
−
1
H
X
+
N
)
−
1
]
s.t.
K
X
=
Tr
(
X
X
H
Ψ
)
Σ
+
σ
n
2
I
Tr
(
X
X
H
)
≤
P
(3)
\begin{aligned} \text { Case 3: } \min _{\mathbf{X}} & \operatorname{Tr}\left[\left(\mathbf{X}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{X}}^{-1} \mathbf{H} \mathbf{X}+\mathbf{N}\right)^{-1}\right] \\ \text { s.t. } & \mathbf{K}_{\mathbf{X}}=\operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}}\right) \leq P \end{aligned} \tag{3}
Case 3: Xmin s.t. Tr[(XHHHKX−1HX+N)−1]KX=Tr(XXHΨ)Σ+σn2ITr(XXH)≤P(3)
同样问题的不同表述为:
Case 4:
min
X
Tr
[
(
(
X
H
H
H
K
X
−
1
H
X
)
⊗
M
+
N
⊗
M
)
−
1
]
s.t.
K
X
=
Tr
(
X
X
H
Ψ
)
Σ
+
σ
n
2
I
Tr
(
X
X
H
)
≤
P
(4)
\begin{aligned} \text { Case 4: } \min _{\mathbf{X}} & \operatorname{Tr}\left[\left(\left(\mathbf{X}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{X}}^{-1} \mathbf{H} \mathbf{X}\right) \otimes \mathbf{M}+\mathbf{N} \otimes \mathbf{M}\right)^{-1}\right] \\ \text { s.t. } & \mathbf{K}_{\mathbf{X}}=\operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}}\right) \leq P \end{aligned} \tag{4}
Case 4: Xmin s.t. Tr[((XHHHKX−1HX)⊗M+N⊗M)−1]KX=Tr(XXHΨ)Σ+σn2ITr(XXH)≤P(4)
Dual-Hop AF MIMO Relaying Systems
Case 5:
min
X
log
∣
A
H
(
X
H
H
H
K
X
−
1
H
X
+
I
)
−
1
A
+
N
∣
s.t.
K
X
=
Tr
(
X
X
H
Ψ
)
Σ
+
σ
n
2
I
Tr
(
X
X
H
)
≤
P
(5)
\begin{aligned} \text { Case 5: } \min _{\mathbf{X}} & \log \left|\mathbf{A}^{\mathrm{H}}\left(\mathbf{X}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{X}}^{-1} \mathbf{H} \mathbf{X}+\mathbf{I}\right)^{-1} \mathbf{A}+\mathbf{N}\right| \\ \text { s.t. } & \mathbf{K}_{\mathbf{X}}=\operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}}\right) \leq P \end{aligned} \tag{5}
Case 5: Xmin s.t. log∣∣∣AH(XHHHKX−1HX+I)−1A+N∣∣∣KX=Tr(XXHΨ)Σ+σn2ITr(XXH)≤P(5)
Case 6:
min
X
Tr
[
A
H
(
X
H
H
H
K
X
−
1
H
X
+
I
)
−
1
A
]
s.t.
K
X
=
Tr
(
X
X
H
Ψ
)
Σ
+
σ
n
2
I
Tr
(
X
X
H
)
≤
P
(6)
\begin{aligned} \text { Case 6: } \min _{\mathbf{X}} & \operatorname{Tr}\left[\mathbf{A}^{\mathrm{H}}\left(\mathbf{X}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{X}}^{-1} \mathbf{H} \mathbf{X}+\mathbf{I}\right)^{-1} \mathbf{A}\right] \\ \text { s.t. } & \mathbf{K}_{\mathbf{X}}=\operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}}\right) \leq P \end{aligned} \tag{6}
Case 6: Xmin s.t. Tr[AH(XHHHKX−1HX+I)−1A]KX=Tr(XXHΨ)Σ+σn2ITr(XXH)≤P(6)
Additive Majorization Theory Based Optimization Problems
Case 7:
min
X
f
A-Schur
(
d
[
(
X
H
H
H
K
X
−
1
H
X
+
I
)
−
1
]
)
s.t.
K
X
=
Tr
(
X
X
H
Ψ
)
Σ
+
σ
n
2
I
Tr
(
X
X
H
)
≤
P
(7)
\begin{aligned} \text { Case 7: } \min _{\mathbf{X}} & f_{\text {A-Schur }}\left(\mathbf{d}\left[\left(\mathbf{X}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{X}}^{-1} \mathbf{H} \mathbf{X}+\mathbf{I}\right)^{-1}\right]\right) \\ & \text { s.t. } \mathbf{K}_{\mathbf{X}}=\operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}}\right) \leq P \end{aligned} \tag{7}
Case 7: XminfA-Schur (d[(XHHHKX−1HX+I)−1]) s.t. KX=Tr(XXHΨ)Σ+σn2ITr(XXH)≤P(7)
其中
f
A-Schur
(
∙
)
f_{\text {A-Schur }}(\bullet)
fA-Schur (∙)为增函数且为加性舒尔凹(或凸)。
Multiplicative Majorization Theory Based Optimization Problems
Case 8: min f M-Schur ( d 2 [ L ] ) s.t. ( X H H H K X − 1 H X + I ) − 1 = L L H K X = Tr ( X X H Ψ ) Σ + σ n 2 I Tr ( X X H ) ≤ P (8) \begin{aligned} \text { Case 8: } \min & f_{\text {M-Schur }}\left(\mathbf{d}^{2}[\mathbf{L}]\right) \\ \text { s.t. } &\left(\mathbf{X}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{X}}^{-1} \mathbf{H} \mathbf{X}+\mathbf{I}\right)^{-1}=\mathbf{L} \mathbf{L}^{\mathrm{H}} \\ & \mathbf{K}_{\mathbf{X}}=\operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{X X}^{\mathrm{H}}\right) \leq P \end{aligned} \tag{8} Case 8: min s.t. fM-Schur (d2[L])(XHHHKX−1HX+I)−1=LLHKX=Tr(XXHΨ)Σ+σn2ITr(XXH)≤P(8)
MVOP问题
以上MIMO优化问题都可以写为如下形式:
MVOP 1:
min
X
f
Matrix
(
K
X
−
1
/
2
H
X
)
s.t.
K
X
=
Tr
(
X
X
H
Ψ
)
Σ
+
σ
n
2
I
Tr
(
X
X
H
)
≤
P
(9)
\begin{aligned} \text { MVOP 1: } \min _{\mathbf{X}} & f_{\text {Matrix }}\left(\mathbf{K}_{\mathbf{X}}^{-1 / 2} \mathbf{H X}\right) \\ \text { s.t. } & \mathbf{K}_{\mathbf{X}}=\operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{X} \mathbf{X}^{\mathrm{H}}\right) \leq P \end{aligned} \tag{9}
MVOP 1: Xmin s.t. fMatrix (KX−1/2HX)KX=Tr(XXHΨ)Σ+σn2ITr(XXH)≤P(9)
由辅助酉阵
Q
\mathbf{Q}
Q,用
F
\mathbf{F}
F换
X
\mathbf{X}
X,原问题转化为:
MVOP 2:
min
F
f
Vector
[
λ
(
F
H
H
H
K
F
−
1
H
F
)
]
s.t.
K
F
=
Tr
(
F
F
H
Ψ
)
Σ
+
σ
n
2
I
Tr
(
F
F
H
)
≤
P
,
(10)
\begin{aligned} \text { MVOP 2: } \min _{\mathbf{F}} & f_{\text {Vector }}\left[\boldsymbol{\lambda}\left(\mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H F}\right)\right] \\ \text { s.t. } & \mathbf{K}_{\mathbf{F}}=\operatorname{Tr}\left(\mathbf{F} \mathbf{F}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{F} \mathbf{F}^{\mathrm{H}}\right) \leq P, \end{aligned} \tag{10}
MVOP 2: Fmin s.t. fVector [λ(FHHHKF−1HF)]KF=Tr(FFHΨ)Σ+σn2ITr(FFH)≤P,(10)
其中
λ
(
Z
)
=
[
λ
1
(
Z
)
,
…
,
λ
N
(
Z
)
]
T
\boldsymbol{\lambda}(\mathbf{Z})=\left[\lambda_{1}(\mathbf{Z}), \ldots, \lambda_{N}(\mathbf{Z})\right]^{\mathrm{T}}
λ(Z)=[λ1(Z),…,λN(Z)]T为
Z
\mathbf{Z}
Z特征值向量的降序排列;
f
A-Schur
(
∙
)
f_{\text {A-Schur }}(\bullet)
fA-Schur (∙)是递减的向量函数。这个过程则是MVOP问题基于酉阵的转化过程。且最优的酉阵有如下功效:
f
Matrix
(
K
F
−
1
/
2
H
F
Q
)
≥
f
Matrix
(
K
F
−
1
/
2
H
F
Q
opt
)
=
f
Vector
[
λ
(
F
H
H
H
K
F
−
1
H
F
)
]
(11)
\begin{aligned} f_{\text {Matrix }}\left(\mathbf{K}_{\mathbf{F}}^{-1 / 2} \mathbf{H F Q}\right) & \geq f_{\text {Matrix }}\left(\mathbf{K}_{\mathbf{F}}^{-1 / 2} \mathbf{H F Q}_{\text {opt }}\right) \\ &=f_{\text {Vector }}\left[\boldsymbol{\lambda}\left(\mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H F}\right)\right] \end{aligned} \tag{11}
fMatrix (KF−1/2HFQ)≥fMatrix (KF−1/2HFQopt )=fVector [λ(FHHHKF−1HF)](11)后续也有寻找酉阵的过程,而最优的酉阵可以完成从问题
9
9
9到问题
10
10
10的转化。
III.MMOP优化框架
MMOP (Matrix-Monotonic Optimization Problem)
首先什么是MMOP问题,所谓矩阵单调,是定义在广义不等式上的矩阵函数运算。矩阵单调性如下: N ⪰ M ⪰ 0 → f ( N ) ≤ f ( M ) (12) \boldsymbol{N} \succeq \boldsymbol{M} \succeq \mathbf{0} \rightarrow f(\boldsymbol{N}) \leq f(\boldsymbol{M}) \tag{12} N⪰M⪰0→f(N)≤f(M)(12)由矩阵单调性的定义和半正定矩阵的形式可知,上述的 f Vector [ λ ( ∙ ) ] f_{\text {Vector }}[\boldsymbol{\lambda}(\boldsymbol{\bullet})] fVector [λ(∙)]有如下的两个性质,这两个性质可以保证我们基于单调性得到最优解:
- f Vector [ λ ( M ) ] f_{\text {Vector }}[\boldsymbol{\lambda}(\boldsymbol{\boldsymbol{M}})] fVector [λ(M)]关于 M \boldsymbol{M} M单调递减
- f Vector [ λ ( M ) ] f_{\text {Vector }}[\boldsymbol{\lambda}(\boldsymbol{\boldsymbol{M}})] fVector [λ(M)]与 M \boldsymbol{M} M特征值对应的特征向量无关。
基于此,上述
(
10
)
(10)
(10)的MVOP问题可以被转化为MMOP问题:
M
M
O
P
1
:
max
F
F
H
H
H
K
F
−
1
H
F
s.t.
K
F
=
Tr
(
F
F
H
Ψ
)
Σ
+
σ
n
2
I
Tr
(
F
F
H
)
≤
P
(13)
\begin{aligned} \mathbf{M M O P} 1: \max _{\mathbf{F}} & \mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H} \mathbf{F} \\ \text { s.t. } & \mathbf{K}_{\mathbf{F}}=\operatorname{Tr}\left(\mathbf{F} \mathbf{F}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left(\mathbf{F} \mathbf{F}^{\mathrm{H}}\right) \leq P \end{aligned} \tag{13}
MMOP1:Fmax s.t. FHHHKF−1HFKF=Tr(FFHΨ)Σ+σn2ITr(FFH)≤P(13)从
(
10
)
(10)
(10)到
(
13
)
(13)
(13),我们将外层函数
f
Vector
[
λ
(
∙
)
]
f_{\text {Vector }}[\boldsymbol{\lambda}(\boldsymbol{\bullet})]
fVector [λ(∙)]去掉了,目标函数从一个标量函数变成了一个矩阵函数,这是这一步变换最重要的特征。而
(
13
)
(13)
(13)在优化中算多目标优化问题,在Boyd教材中有所讨论。
这一步等价变换有一个很重要的问题是,从
(
11
)
(11)
(11)到
(
13
)
(13)
(13)该问题被泛化了,换而言之,对于一个特定的外层函数
f
Vector
[
λ
(
∙
)
]
f_{\text {Vector }}[\boldsymbol{\lambda}(\boldsymbol{\bullet})]
fVector [λ(∙)],其可行域是问题
(
13
)
(13)
(13)的一个子集。所有满足上述条件的外层函数,若其内核与限制条件不变,解出的矩阵最优解应该是一样的。
下面对问题 ( 13 ) (13) (13)进行讨论。问题 ( 13 ) (13) (13)有如下两个性质:
- 问题 ( 13 ) (13) (13)的目标函数,最大化 F H H H K F − 1 H F \mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H} \mathbf{F} FHHHKF−1HF等价于最大化其特征值向量,即 λ ( F H H H K F − 1 H F ) \boldsymbol{\lambda}\left(\mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H} \mathbf{F}\right) λ(FHHHKF−1HF)。
- 问题 ( 13 ) (13) (13)具有酉不变性,即如果有最优解 F opt \mathbf{F}_{\text {opt }} Fopt ,对其做任意酉变换 F o p t U \mathbf{F}_{\mathrm{opt}} \mathbf{U} FoptU仍然是最优解。
这两个性质使得问题 ( 13 ) (13) (13)相对较友好,也是下面讨论的基础。基于上述性质,以及下面的两个引理:
- 如果 f Monot ( ∙ ) f_{\text {Monot }}(\bullet) fMonot (∙)是关于半正定阵的单调函数,最优解应在其边界上,所以最优解的必要条件是: Tr ( F F H ) = P (14) \operatorname{Tr}\left(\mathbf{F} \mathbf{F}^{\mathrm{H}}\right)=P \tag{14} Tr(FFH)=P(14)
- 定义一个辅助变量: η f ≜ Tr ( F F H Ψ ) α + σ n 2 with α = λ min ( Σ ) (15) \eta_{f} \triangleq \operatorname{Tr}\left(\mathbf{F} \mathbf{F}^{\mathrm{H}} \boldsymbol{\Psi}\right) \alpha+\sigma_{n}^{2} \text { with } \alpha=\lambda_{\min }(\boldsymbol{\Sigma}) \tag{15} ηf≜Tr(FFHΨ)α+σn2 with α=λmin(Σ)(15)则 ( 14 ) (14) (14)可以等价于 Tr [ F F H ( α P Ψ + σ n 2 I ) ] / η f = P \operatorname{Tr}\left[\mathbf{F} \mathbf{F}^{\mathrm{H}}\left(\alpha P \boldsymbol{\Psi}+\sigma_{n}^{2} \mathbf{I}\right)\right] / \eta_{f}=P Tr[FFH(αPΨ+σn2I)]/ηf=P
基于上述引理1与引理2,优化问题
(
13
)
(13)
(13)可以被转化为
MMOP 2:
max
F
F
H
H
H
K
F
−
1
H
F
s.t.
K
F
=
Tr
(
F
F
H
Ψ
)
Σ
+
σ
n
2
I
Tr
[
F
F
H
(
α
P
Ψ
+
σ
n
2
I
)
]
/
η
f
=
P
(16)
\begin{aligned} \text { MMOP 2: } \max _{\mathbf{F}} & \mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H F} \\ \text { s.t. } & \mathbf{K}_{\mathbf{F}}=\operatorname{Tr}\left(\mathbf{F} \mathbf{F}^{\mathrm{H}} \boldsymbol{\Psi}\right) \boldsymbol{\Sigma}+\sigma_{n}^{2} \mathbf{I} \\ & \operatorname{Tr}\left[\mathbf{F} \mathbf{F}^{\mathrm{H}}\left(\alpha P \boldsymbol{\Psi}+\sigma_{n}^{2} \mathbf{I}\right)\right] / \eta_{f}=P \end{aligned} \tag{16}
MMOP 2: Fmax s.t. FHHHKF−1HFKF=Tr(FFHΨ)Σ+σn2ITr[FFH(αPΨ+σn2I)]/ηf=P(16)为了进一步将约束条件2简化,定义如下辅助变量:
F
~
≜
1
/
η
f
(
α
P
Ψ
+
σ
n
2
I
)
1
/
2
F
(17)
\tilde{\mathbf{F}} \triangleq 1 / \sqrt{\eta_{f}}\left(\alpha P \boldsymbol{\Psi}+\sigma_{n}^{2} \mathbf{I}\right)^{1 / 2} \mathbf{F} \tag{17}
F~≜1/ηf(αPΨ+σn2I)1/2F(17)则问题
(
16
)
(16)
(16)进一步转化为
max
F
~
F
~
H
(
α
P
Ψ
+
σ
n
2
I
)
−
1
/
2
H
H
(
K
F
/
η
f
)
−
1
H
×
(
α
P
Ψ
+
σ
n
2
I
)
−
1
/
2
F
~
s.t.
Tr
(
F
~
F
~
H
)
=
P
(18)
\begin{array}{l} \max _{\tilde{\mathbf{F}}} \tilde{\mathbf{F}}^{\mathrm{H}}\left(\alpha P \boldsymbol{\Psi}+\sigma_{n}^{2} \mathbf{I}\right)^{-1 / 2} \mathbf{H}^{\mathrm{H}}\left(\mathbf{K}_{\mathbf{F}} / \eta_{f}\right)^{-1} \mathbf{H} \times\left(\alpha P \boldsymbol{\Psi}+\sigma_{n}^{2} \mathbf{I}\right)^{-1 / 2} \tilde{\mathbf{F}} \\ \text { s.t. } \operatorname{Tr}\left(\tilde{\mathbf{F}} \tilde{\mathbf{F}}^{\mathbf{H}}\right)=P \end{array} \tag{18}
maxF~F~H(αPΨ+σn2I)−1/2HH(KF/ηf)−1H×(αPΨ+σn2I)−1/2F~ s.t. Tr(F~F~H)=P(18)到这里作者进一步做了一个简化,当
Ψ
\boldsymbol{\Psi}
Ψ与
Σ
\boldsymbol{\Sigma}
Σ不正比于单位阵
I
\boldsymbol{I}
I时,问题会变得很复杂,所以做了
Ψ
∝
I
or
Σ
∝
I
\boldsymbol{\Psi} \propto \mathbf{I} \text { or } \boldsymbol{\Sigma} \propto \mathbf{I}
Ψ∝I or Σ∝I这样的假设,并且举例说明了这样的假设在实际工程应用中是有实际意义的。
基于上面正比于单位阵的假设,进一步定义下面两个恒定矩阵:
K
Ψ
≜
P
λ
max
(
Ψ
)
P
λ
max
(
Ψ
)
α
+
σ
n
2
Σ
+
σ
n
2
P
λ
max
(
Ψ
)
α
+
σ
n
2
I
Π
≜
K
Ψ
−
1
/
2
H
(
α
P
Ψ
+
σ
n
2
I
)
−
1
/
2
(19)
\begin{aligned} \mathbf{K}_{\Psi} & \triangleq \frac{P \lambda_{\max }(\boldsymbol{\Psi})}{P \lambda_{\max }(\boldsymbol{\Psi}) \alpha+\sigma_{n}^{2}} \boldsymbol{\Sigma}+\frac{\sigma_{n}^{2}}{P \lambda_{\max }(\boldsymbol{\Psi}) \alpha+\sigma_{n}^{2}} \mathbf{I} \\ \boldsymbol{\Pi} & \triangleq \mathbf{K}_{\boldsymbol{\Psi}}^{-1 / 2} \mathbf{H}\left(\alpha P \boldsymbol{\Psi}+\sigma_{n}^{2} \mathbf{I}\right)^{-1 / 2} \end{aligned} \tag{19}
KΨΠ≜Pλmax(Ψ)α+σn2Pλmax(Ψ)Σ+Pλmax(Ψ)α+σn2σn2I≜KΨ−1/2H(αPΨ+σn2I)−1/2(19)问题
(
18
)
(18)
(18)进一步可以简化:
MMOP 3:
max
F
~
F
~
H
Π
H
Π
F
~
s.t.
Tr
(
F
~
F
~
H
)
≤
P
.
(20)
\begin{aligned} \text { MMOP 3: } & \max _{\tilde{\mathbf{F}}}{\tilde{\mathbf{F}}}^{\mathrm{H}} \boldsymbol{\Pi}^{\mathrm{H}} \boldsymbol{\Pi} \tilde{\mathbf{F}} \\ & \text { s.t. } \operatorname{Tr}\left(\tilde{\mathbf{F}} \tilde{\mathbf{F}}^{\mathrm{H}}\right) \leq P . \end{aligned} \tag{20}
MMOP 3: F~maxF~HΠHΠF~ s.t. Tr(F~F~H)≤P.(20)该问题目标函数关于半正定的优化变量
F
~
\tilde{\mathbf{F}}
F~是凸的,所以松弛之后有如下问题:
max
ρ
,
F
~
ρ
s.t.
Tr
(
F
~
F
~
H
)
≤
P
F
~
H
Π
H
Π
F
~
=
ρ
F
~
i
n
H
Π
H
Π
F
~
i
n
(21)
\begin{array}{ll} \max _{\rho, \tilde{\mathbf{F}}} & \rho \\ \text { s.t. } & \operatorname{Tr}\left(\tilde{\mathbf{F}} \tilde{\mathbf{F}}^{\mathrm{H}}\right) \leq P \\ & \tilde{\mathbf{F}}^{\mathrm{H}} \mathbf{\Pi}^{\mathrm{H}} \mathbf{\Pi} \tilde{\mathbf{F}}=\rho \tilde{\mathbf{F}}_{\mathrm{in}}^{\mathrm{H}} \mathbf{\Pi}^{\mathrm{H}} \mathbf{\Pi} \tilde{\mathbf{F}}_{\mathrm{in}} \end{array} \tag{21}
maxρ,F~ s.t. ρTr(F~F~H)≤PF~HΠHΠF~=ρF~inHΠHΠF~in(21)其中
F
~
i
n
\tilde{\mathbf{F}}_{\mathrm{in}}
F~in是范数球的一个内点,即
Tr
(
F
~
in
F
~
in
H
)
<
P
\operatorname{Tr}\left(\tilde{\mathbf{F}}_{\text {in }} \tilde{\mathbf{F}}_{\text {in }}^{\mathrm{H}}\right)<P
Tr(F~in F~in H)<P。由一系列的结论和引理(具体过程见原论文),该问题的解的结构由如下定理给出:
若有 Ψ ∝ I or Σ ∝ I \boldsymbol{\Psi} \propto \mathbf{I} \text { or } \boldsymbol{\Sigma} \propto \mathbf{I} Ψ∝I or Σ∝I, ( 13 ) (13) (13)中的最优解有如下形式: F o p t = η f ( α P Ψ + σ n 2 I ) − 1 / 2 V Π Λ F U A r b H (22) \mathbf{F}_{\mathrm{opt}}=\sqrt{\eta_{f}}\left(\alpha P \boldsymbol{\Psi}+\sigma_{n}^{2} \mathbf{I}\right)^{-1 / 2} \mathbf{V}_{\mathbf{\Pi}} \boldsymbol{\Lambda}_{\mathbf{F}} \mathbf{U}_{\mathrm{Arb}}^{\mathrm{H}} \tag{22} Fopt=ηf(αPΨ+σn2I)−1/2VΠΛFUArbH(22)其中 η f \eta_{f} ηf有 η f = P / Tr [ ( α P Ψ + σ n 2 I ) − 1 V Π Λ F Λ F H V Π H ] \eta_{f}=P / \operatorname{Tr}\left[\left(\alpha P \boldsymbol{\Psi}+\sigma_{n}^{2} \mathbf{I}\right)^{-1} \mathbf{V}_{\mathbf{\Pi}} \boldsymbol{\Lambda}_{\mathbf{F}} \boldsymbol{\Lambda}_{\mathbf{F}}^{\mathrm{H}} \mathbf{V}_{\mathbf{\Pi}}^{\mathrm{H}}\right] ηf=P/Tr[(αPΨ+σn2I)−1VΠΛFΛFHVΠH], V Π \mathbf{V}_{\mathbf{\Pi}} VΠ有 Π = U Π Λ Π V Π H with Λ Π ↘ \mathbf{\Pi}=\mathbf{U}_{\mathbf{\Pi}} \boldsymbol{\Lambda}_{\mathbf{\Pi}} \mathbf{V}_{\mathbf{\Pi}}^{\mathrm{H}} \text { with } \boldsymbol{\Lambda}_{\Pi} \searrow Π=UΠΛΠVΠH with ΛΠ↘, U A r b \mathbf{U}_{\mathrm{Arb}} UArb为任意酉阵。
所以下面的问题就是两个,第一 Λ F \boldsymbol{\Lambda}_{\mathbf{F}} ΛF如何求解;第二,由问题 ( 10 ) (10) (10)转为问题 ( 13 ) (13) (13)的酉变换怎么找。
Λ F \boldsymbol{\Lambda}_{\mathbf{F}} ΛF的求解
对于
Λ
F
\boldsymbol{\Lambda}_{\mathbf{F}}
ΛF的求解大部分基于KKT条件或者注水法。
首先基于上面的引理结论:
Π
=
U
Π
Λ
Π
V
Π
H
with
Λ
Π
↘
\boldsymbol{\Pi}=\mathbf{U}_{\mathbf{\Pi}} \boldsymbol{\Lambda}_{\mathbf{\Pi}} \mathbf{V}_{\mathbf{\Pi}}^{\mathrm{H}} \text { with } \boldsymbol{\Lambda}_{\mathbf{\Pi}} \searrow
Π=UΠΛΠVΠH with ΛΠ↘
F
~
o
p
t
=
V
Π
Λ
F
U
A
r
b
H
\tilde{\mathbf{F}}_{\mathrm{opt}}=\mathbf{V}_{\mathbf{\Pi}} \boldsymbol{\Lambda}_{\mathbf{F}} \mathbf{U}_{\mathrm{Arb}}^{\mathrm{H}}
F~opt=VΠΛFUArbH且
Λ
F
\boldsymbol{\Lambda}_{\mathbf{F}}
ΛF满足
Λ
F
H
Λ
Π
H
Λ
Π
Λ
F
↘
\boldsymbol{\Lambda}_{\mathbf{F}}^{\mathrm{H}} \boldsymbol{\Lambda}_{\mathbf{\Pi}}^{\mathrm{H}} \boldsymbol{\Lambda}_{\mathbf{\Pi}} \boldsymbol{\Lambda}_{\mathbf{F}} \searrow
ΛFHΛΠHΛΠΛF↘,问题
(
20
)
(20)
(20)可以被写为:
MOLP
:
max
{
f
n
2
}
[
λ
Π
,
1
2
f
1
2
,
…
,
λ
Π
,
N
2
f
N
2
]
T
s.t.
λ
Π
,
1
2
f
1
2
≥
λ
Π
,
2
2
f
2
2
⋯
≥
λ
Π
,
N
2
f
N
2
∑
n
=
1
N
f
n
2
≤
P
.
(23)
\begin{aligned} \operatorname{MOLP}: \max _{\left\{f_{n}^{2}\right\}} &\left[\lambda_{\mathbf{\Pi}, 1}^{2} f_{1}^{2}, \ldots, \lambda_{\mathbf{\Pi}, N}^{2} f_{N}^{2}\right]^{\mathrm{T}} \\ \text { s.t. } & \lambda_{\mathbf{\Pi}, 1}^{2} f_{1}^{2} \geq \lambda_{\mathbf{\Pi}, 2}^{2} f_{2}^{2} \cdots \geq \lambda_{\mathbf{\Pi}, N}^{2} f_{N}^{2} \\ & \sum_{n=1}^{N} f_{n}^{2} \leq P . \end{aligned} \tag{23}
MOLP:{fn2}max s.t. [λΠ,12f12,…,λΠ,N2fN2]TλΠ,12f12≥λΠ,22f22⋯≥λΠ,N2fN2n=1∑Nfn2≤P.(23)为了拿到
Λ
F
\boldsymbol{\Lambda}_{\mathbf{F}}
ΛF闭式解,进一步转化为如下形式:
min
{
f
n
2
}
[
1
λ
Π
,
1
2
f
1
2
+
1
,
…
,
1
λ
Π
,
N
2
f
N
2
+
1
]
T
s.t.
∑
n
=
1
N
f
n
2
≤
P
(24)
\begin{array}{l} \min _{\left\{f_{n}^{2}\right\}}\left[\frac{1}{\lambda_{\mathbf{\Pi}, 1}^{2} f_{1}^{2}+1}, \ldots, \frac{1}{\lambda_{\mathbf{\Pi}, N}^{2} f_{N}^{2}+1}\right]^{\mathrm{T}} \\ \text { s.t. } \sum_{n=1}^{N} f_{n}^{2} \leq P \end{array} \tag{24}
min{fn2}[λΠ,12f12+11,…,λΠ,N2fN2+11]T s.t. ∑n=1Nfn2≤P(24)改写为该式的好处是其中线性分式函数是关于优化变量的凸函数,对于凸的向量优化问题解Pareto最优解,基本方法是标量化:
min
{
f
n
2
}
∑
n
w
n
λ
Π
,
n
2
f
n
2
+
1
s.t.
∑
k
=
1
K
f
k
2
≤
P
,
(25)
\begin{array}{l} \min _{\left\{f_{n}^{2}\right\}} \sum_{n} \frac{w_{n}}{\lambda_{\Pi, n}^{2} f_{n}^{2}+1} \\ \text { s.t. } \sum_{k=1}^{K} f_{k}^{2} \leq P, \end{array} \tag{25}
min{fn2}∑nλΠ,n2fn2+1wn s.t. ∑k=1Kfk2≤P,(25)其中
w
n
w_n
wn为权重系数,
f
n
f_{n}
fn的大小关系被带入
w
n
w_n
wn的选取中考虑,该问题可由注水法得闭式解:
f
n
2
=
(
w
n
λ
Π
,
n
2
μ
−
1
λ
Π
,
n
2
)
+
(26)
f_{n}^{2}=\left(\sqrt{\frac{w_{n}}{\lambda_{\mathbf{\Pi}, n}^{2} \mu}}-\frac{1}{\lambda_{\mathbf{\Pi}, n}^{2}}\right)^{+} \tag{26}
fn2=(λΠ,n2μwn−λΠ,n21)+(26)权重系数选取考虑下式:
w
1
λ
Π
,
1
2
≥
w
2
λ
Π
,
2
2
⋯
≥
w
N
λ
Π
,
N
2
(27)
\sqrt{w_{1} \lambda_{\mathbf{\Pi}, 1}^{2}} \geq \sqrt{w_{2} \lambda_{\mathbf{\Pi}, 2}^{2}} \cdots \geq \sqrt{w_{N} \lambda_{\mathbf{\Pi}, N}^{2}} \tag{27}
w1λΠ,12≥w2λΠ,22⋯≥wNλΠ,N2(27)
酉变换的确定
这里酉变换有两种方法,基于矩阵不等式的方法以及基于Majorization theory的方法。具体推到证明过程见"A. W. Marshall and I. Olkin, Inequalities: Theory of Majorization andIts Applications. New York, NY, USA: Academic, 1979."及参考文献,此处仅记录一些常用问题的结论:
首先有:
A
=
U
A
Λ
A
V
A
H
with
Λ
A
↘
A
N
−
1
A
H
=
U
ANA
Λ
ANA
U
ANA
H
with
Λ
ANA
N
=
U
‾
N
Λ
ˉ
N
U
‾
N
H
with
Λ
ˉ
N
↗
F
H
H
H
K
F
−
1
H
F
=
U
F
H
F
Λ
F
H
F
U
F
H
F
H
with
Λ
F
H
F
↘
(28)
\begin{aligned} \mathbf{A} &=\mathbf{U}_{\mathbf{A}} \Lambda_{\mathrm{A}} \mathbf{V}_{\mathrm{A}}^{\mathrm{H}} \text { with } \boldsymbol{\Lambda}_{\mathbf{A}} \searrow \\ \mathbf{A} \mathbf{N}^{-1} \mathbf{A}^{\mathrm{H}} &=\mathbf{U}_{\text {ANA }} \Lambda_{\text {ANA }} \mathbf{U}_{\text {ANA }}^{\mathrm{H}} \text { with } \Lambda_{\text {ANA }} \\ \mathbf{N} &=\overline{\mathbf{U}}_{\mathbf{N}} \bar{\Lambda}_{\mathrm{N}} \overline{\mathbf{U}}_{\mathbb{N}}^{\mathrm{H}} \text { with } \bar{\Lambda}_{\mathrm{N}} \nearrow \\ \mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H F} &=\mathbf{U}_{\mathbf{F H F}} \Lambda_{\mathrm{FHF}} \mathbf{U}_{\mathbf{F H F}}^{\mathrm{H}} \text { with } \Lambda_{\mathrm{FHF}} \searrow \end{aligned} \tag{28}
AAN−1AHNFHHHKF−1HF=UAΛAVAH with ΛA↘=UANA ΛANA UANA H with ΛANA =UNΛˉNUNH with ΛˉN↗=UFHFΛFHFUFHFH with ΛFHF↘(28)结论:
- Case 1: f Matrix ( ∙ ) = − log ∣ Q H F H H H K F − 1 H F Q + N ∣ Solution: Q o p t = U F H F U ‾ N H (29) \begin{array}{l} \text { Case 1: } f_{\text {Matrix }}(\bullet)=-\log \left|\mathbf{Q}^{\mathrm{H}} \mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H F Q}+\mathbf{N}\right| \\ \text { Solution: } \quad \mathbf{Q}_{\mathrm{opt}}=\mathbf{U}_{\mathbf{F H F}} \overline{\mathbf{U}}_{\mathbf{N}}^{\mathrm{H}} \end{array} \tag{29} Case 1: fMatrix (∙)=−log∣∣QHFHHHKF−1HFQ+N∣∣ Solution: Qopt=UFHFUNH(29)
- Case 2: f Matrix ( ∙ ) = − log ∣ A H Q H F H H H K F − 1 H F Q A + I ∣ Solution: Q o p t = U F H F U A H . (30) \begin{aligned} \text { Case 2: } f_{\text {Matrix }}(\bullet) &=-\log \left|\mathbf{A}^{\mathrm{H}} \mathbf{Q}^{\mathrm{H}} \mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H F Q} \mathbf{A}+\mathbf{I}\right| \\ \text { Solution: } \mathbf{Q}_{\mathrm{opt}} &=\mathbf{U}_{\mathbf{F H F}} \mathbf{U}_{\mathbf{A}}^{\mathrm{H}} . \end{aligned} \tag{30} Case 2: fMatrix (∙) Solution: Qopt=−log∣∣AHQHFHHHKF−1HFQA+I∣∣=UFHFUAH.(30)
- Case 3: f Matrix ( ∙ ) = Tr [ ( Q H F H H H K F − 1 H F Q + N ) − 1 ] Solution: Q o p t = U F H F U ‾ N H (31) \begin{array}{l} \text { Case 3: } f_{\text {Matrix }}(\bullet)=\operatorname{Tr}\left[\left(\mathbf{Q}^{\mathrm{H}} \mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H F Q}+\mathbf{N}\right)^{-1}\right] \\ \text { Solution: } \quad \mathbf{Q}_{\mathrm{opt}}=\mathbf{U}_{\mathbf{F H F}} \overline{\mathbf{U}}_{\mathbf{N}}^{\mathrm{H}} \end{array} \tag{31} Case 3: fMatrix (∙)=Tr[(QHFHHHKF−1HFQ+N)−1] Solution: Qopt=UFHFUNH(31)
- Case 4: f Matrix ( ∙ ) = Tr [ ( Q H F H H H K F − 1 H F Q ⊗ M + N ⊗ M ) − 1 ] Solution: Q o p t = U F H F U N H (32) \begin{aligned} \text { Case 4: } & f_{\text {Matrix }}(\bullet) \\ =& \operatorname{Tr}\left[\left(\mathbf{Q}^{\mathrm{H}} \mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H F Q} \otimes \mathbf{M}+\mathbf{N} \otimes \mathbf{M}\right)^{-1}\right] \\ \text { Solution: } & \mathbf{Q}_{\mathrm{opt}}=\mathbf{U}_{\mathbf{F} H F}{\mathbf{U}}_{\mathbf{N}}^{\mathrm{H}} \end{aligned} \tag{32} Case 4: = Solution: fMatrix (∙)Tr[(QHFHHHKF−1HFQ⊗M+N⊗M)−1]Qopt=UFHFUNH(32)
- Case 5: f Matrix ( ∙ ) = log ∣ A H ( Q H F H H H K F − 1 H F Q + I ) − 1 A + N ∣ Solution: Q o p t = U F H F U A N A H (33) \begin{aligned} \text { Case 5: } f_{\text {Matrix }}(\bullet) \\ =& \log \left|\mathbf{A}^{\mathrm{H}}\left(\mathbf{Q}^{\mathrm{H}} \mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H F Q}+\mathbf{I}\right)^{-1} \mathbf{A}+\mathbf{N}\right| \\ \text { Solution: } & \mathbf{Q}_{\mathrm{opt}}=\mathbf{U}_{\mathbf{F H F}} \mathbf{U}_{\mathbf{A N A}}^{\mathrm{H}} \end{aligned} \tag{33} Case 5: fMatrix (∙)= Solution: log∣∣∣AH(QHFHHHKF−1HFQ+I)−1A+N∣∣∣Qopt=UFHFUANAH(33)
- Case 6: f Matrix ( ∙ ) = Tr [ A H ( Q H F H H H K F − 1 H F Q + I ) − 1 A ] Solution: Q o p t = U F H F U A H (34) \begin{aligned} \text { Case 6: } f_{\text {Matrix }}(\bullet) \\ =& \operatorname{Tr}\left[\mathbf{A}^{\mathrm{H}}\left(\mathbf{Q}^{\mathrm{H}} \mathbf{F}^{\mathrm{H}} \mathbf{H}^{\mathrm{H}} \mathbf{K}_{\mathbf{F}}^{-1} \mathbf{H F Q}+\mathbf{I}\right)^{-1} \mathbf{A}\right] \\ \text { Solution: } & \mathbf{Q}_{\mathrm{opt}}=\mathbf{U}_{\mathbf{F H F}} \mathbf{U}_{\mathbf{A}}^{\mathrm{H}} \end{aligned} \tag{34} Case 6: fMatrix (∙)= Solution: Tr[AH(QHFHHHKF−1HFQ+I)−1A]Qopt=UFHFUAH(34)