文章目录
为了表述功率约束,定义发射波形的协方差为
R
=
E
(
x
[
n
]
x
H
[
n
]
)
.
\begin{equation*} \boldsymbol {R} = {\textbf E} \left(\boldsymbol {x} [n] \boldsymbol {x} ^ H [n] \right) . \tag{5} \end{equation*}
R=E(x[n]xH[n]).(5)
将(1)-(4)代入(5),得到协方差 R \boldsymbol {R} R 为
R \boldsymbol {R} R as R = W r W r H + W c W c H . \begin{equation*} \boldsymbol {R} = \boldsymbol {W}_r \boldsymbol {W}_r ^ H + \boldsymbol {W}_c \boldsymbol {W}_c ^ H. \tag{6} \end{equation*} R=WrWrH+WcWcH.(6)
The per-antenna power constraint implies that for each m = 1 , … , M \ m = 1,\ldots,M m=1,…,M it holds that [ R ] m , m = [ W r W r H + W c W c H ] m , m = P t / M , \begin{equation*} \left[\boldsymbol {R}\right]_{m,m} = \left[\boldsymbol {W}_r \boldsymbol {W}_r ^ H + \boldsymbol {W}_c \boldsymbol {W}_c ^ H\right]_{m,m} = P_t / M, \tag{7} \end{equation*} [R]m,m=[WrWrH+WcWcH]m,m=Pt/M,(7)
where P t P_t Pt is the total transmit power.
III. PERFORMANCE METRICS OF RADAR AND COMMUNICATION
基于联合发射波束形成的信号模型,我们的预编码器设计遵循以下指导原则:对于MIMO雷达,预编码器设计用于合成针对感兴趣雷达目标的发射波束;对于多用户MIMO通信,预编码器的设计是为了保证通信用户接收信噪比。MIMO雷达和多用户MIMO通信的这些性能指标分别在第III-A节和第III-B节中进行了适当的阐述。
A. MIMO Radar Performance
MIMO雷达波束形成的主要目的是将发射波束引导到几个给定的方向,以便获得这些波束照射目标的更多信息。这些方向通常是发射机所知道的:当雷达工作在跟踪模式时,波束方向是根据先前观测获得的目标方向推断出来的;当雷达工作在搜索模式时,波束方向由角扇形感兴趣中心(angular sector-of-interest)给出。因此,为了制定与MIMO雷达波束形成相关的性能指标,我们首先表示每个方向上的发射信号,然后开发评估发射波束方向的损失函数。将损失函数和天线功率约束相结合,实现了一个关于雷达性能的优化问题。
在 DFRC 系统中,通信信号可用于感知,因为雷达接收机具有发射通信波形的完整知识。这样,通信信号在雷达接收机上就不被认为是干扰。假设发射波形为窄带,传播路径为视距(LoS), θ θ θ 方向的基带信号可以表示为
y [ n ; θ ] = a H ( θ ) x [ n ] , \begin{equation*} y[n; \theta ] = \boldsymbol {a} ^ H (\theta) \boldsymbol {x} [n], \tag{8} \end{equation*} y[n;θ]=aH(θ)x[n],(8)
其中
a
(
θ
)
\boldsymbol {a}(\theta)
a(θ) 是
θ
θ
θ 方向的阵列导向向量。当发射波形被位于角方向
θ
θ
θ 的点目标反射时,接收信号可表示为
y
[
n
]
=
β
a
∗
(
θ
)
a
H
(
θ
)
x
[
n
−
n
′
]
+
z
[
n
]
,
\begin{equation*} \boldsymbol {y} [n] = \beta \boldsymbol {a} ^ *(\theta) \boldsymbol {a} ^ H (\theta) \boldsymbol {x} [n - n^{\prime }] + \boldsymbol {z} [n], \tag{9} \end{equation*}
y[n]=βa∗(θ)aH(θ)x[n−n′]+z[n],(9)
其中 β \beta β 与目标的雷达-散射截面(RCS)成比例的复振幅, n ′ n^{\prime } n′ 代表离散时间延迟,and z [ n ] \boldsymbol {z}[n] z[n] 为协方差 R z \boldsymbol {R} _ z Rz 的加性零均值时间白噪声。
根据文献[50]中关于MIMO雷达探测信号设计的指导原则,MIMO雷达发射波束形成的预期目标包括:
- 优化给定方向的发射功率,或通常匹配期望的波束方向图;
- 降低多个给定方向信号之间的互相关模式,这对自适应MIMO雷达技术的性能至关重要。
这里,角方向 θ θ θ 处的发射功率(波束方向图)为
P ( θ ; R ) = E ( ∣ y [ n ; θ ] ∣ 2 ) = E ( a H ( θ ) x [ n ] x H [ n ] a ( θ ) ) = a H ( θ ) R a ( θ ) , \begin{align*} & P(\theta ; \boldsymbol {R}) = {\textbf E} \left(\left| y[n;\theta ] \right|^2 \right) \\ &= {\textbf E} \left(\boldsymbol {a} ^ H (\theta) \boldsymbol {x} [n] \boldsymbol {x} ^ H [n] \boldsymbol {a} (\theta) \right) = \boldsymbol {a} ^H (\theta) \boldsymbol {R} \boldsymbol {a}(\theta), \tag{10} \end{align*} P(θ;R)=E(∣y[n;θ]∣2)=E(aH(θ)x[n]xH[n]a(θ))=aH(θ)Ra(θ),(10)
θ 1 \theta _1 θ1 and θ 2 \theta _2 θ2 之间的互相关模式(cross correlation pattern)定义为
P c ( θ 1 , θ 2 ; R ) = E ( y ∗ [ n ; θ 1 ] y [ n ; θ 2 ] ) = E ( a H ( θ 2 ) x [ n ] x H [ n ] a ( θ 1 ) ) = a H ( θ 2 ) R a ( θ 1 ) . \begin{align*} & P_c (\theta _1, \theta _2; \boldsymbol {R}) = {\textbf E} \left(y ^ *[ n; \theta _1 ] y[n; \theta _2] \right) \\ &= {\textbf E} \left(\boldsymbol {a} ^ H (\theta _2) \boldsymbol {x} [n] \boldsymbol {x} ^ H [n] \boldsymbol {a} (\theta _1) \right) = \boldsymbol {a} ^H (\theta _2) \boldsymbol {R} \boldsymbol {a}(\theta _1). \tag{11} \end{align*} Pc(θ1,θ2;R)=E(y∗[n;θ1]y[n;θ2])=E(aH(θ2)x[n]xH[n]a(θ1))=aH(θ2)Ra(θ1).(11)
由式(10)和式(11)可知,发射波束方向图和相互关联方向图均由协方差 R \boldsymbol {R} R 决定,因此,通过设计协方差矩阵 R \boldsymbol {R} R [50],[54],可以实现MIMO雷达波形的适当波束形成。
为此,我们使用在[50],[54]中提出的损失函数来评估雷达性能,它是波束方向图误差和互相关两部分的加权和。具体来说,第一部分是得到的波束方向图与某些期望波束方向图之间的均方误差(MSE),由
L r , 1 ( R , α ) = 1 L ∑ l = 1 L ∣ α d ( θ l ) − P ( θ l ; R ) ∣ 2 , \begin{equation*} L_{r,1}(\boldsymbol {R}, \alpha) = \frac{1}{L} \sum _ {l =1} ^ L \left| \alpha d(\theta _l) - P(\theta _l; \boldsymbol {R}) \right| ^ 2, \tag{12} \end{equation*} Lr,1(R,α)=L1l=1∑L∣αd(θl)−P(θl;R)∣2,(12)
其中 α \alpha α 是缩放因子, d ( θ ) d(\theta) d(θ) 是给定的期望波束方向图, { θ l } l = 1 L \lbrace \theta _l \rbrace _{l=1} ^ {L} {θl}l=1L 是采样的角度网格。第二部分为均方互相关模式,表示为
L r , 2 ( R ) = 2 P 2 − P ∑ p = 1 P − 1 ∑ q = p + 1 P ∣ P c ( θ ‾ p , θ ‾ q ; R ) ∣ 2 , \begin{equation*} L_{r,2}(\boldsymbol {R}) = \frac{2 }{P^2 - P} \sum _{p = 1}^{P-1} \sum _{q = p+1}^{P} \left| P_c(\overline{ \theta } _p, \overline{ \theta } _q ; \boldsymbol {R}) \right| ^ 2, \tag{13} \end{equation*} Lr,2(R)=P2−P2p=1∑P−1q=p+1∑P Pc(θp,θq;R) 2,(13)
其 { θ ‾ p } p = 1 P \lbrace \overline{ \theta } _p \rbrace _{p=1} ^ {P} {θp}p=1P 为目标的给定方向。(13)中的求和被 2 P 2 − P \frac{2 }{P^2 - P} P2−P2 归一化 ,因为集合 { θ ‾ p } \lbrace \overline{ \theta } _p \rbrace {θp}中存在 P 2 − P 2 \frac{P^2 - P}{2} 2P2−P对不同的方向。雷达的损失函数为
L r ( R , α ) = L r , 1 ( R , α ) + w c L r , 2 ( R ) , \begin{equation*} L_{r}(\boldsymbol {R}, \alpha) = L_{r,1}(\boldsymbol {R}, \alpha) + w_c L_{r,2}(\boldsymbol {R}), \tag{14} \end{equation*} Lr(R,α)=Lr,1(R,α)+wcLr,2(R),(14)
其中 w c w_c wc 为加权因子。如[50]、[54]所述,损失函数 L r ( R , α ) L_{r}(\boldsymbol {R}, \alpha) Lr(R,α) 可以写成 R \boldsymbol {R} R 和 α \alpha α 的正半定二次函数。
结合(12)和(13)中的损失函数,在没有通信约束的情况下,即在只有雷达的设定下,可以根据每个天线功率约束下的总体雷达目标[50],设计发射信号的协方差。
min R , α L r ( R , α ) subject to R ∈ S M + , [ R ] m , m = P t / M , m = 1 , … , M . \begin{align*} \mathop{\text{min}} _{\boldsymbol {R}, \alpha } \ \ \ \ & \ L_r(\boldsymbol {R}, \alpha) \tag{15a} \\ \text{subject to} \ & \ \boldsymbol {R} \in \mathcal {S}_{M}^+, \tag{15b} \\ & \ \left[\boldsymbol {R}\right]_{m,m} = P_t / M, \ m = 1,\ldots,M. \tag{15c} \end{align*} minR,α subject to Lr(R,α) R∈SM+, [R]m,m=Pt/M, m=1,…,M.(15a)(15b)(15c)
我们用 R 0 \boldsymbol {R}_0 R0 表示该问题的最优协方差。通常,如果发射波形的协方差为 R 0 \boldsymbol {R}_0 R0,则无法满足第III-B节所述的多用户MIMO通信的性能要求。换句话说,与仅雷达情况相比,由于与通信共享频谱,存在固有的雷达性能损失。为了解决我们的DFRC系统的通信性能,我们将在下一小节中讨论通信度量指标。
B. Multiuser MIMO Communication Performance
多用户广播通信的常见性能指标是SINR,它直接关系到在降低复杂度解码下(reduced complexity decoding)的可实现速率[55,第8章]。因此,我们设计了MIMO传输的预编码器,以优化用户的SINR。为此,我们首先提出了通信信号模型,并推导了SINR的表达式。
考虑一个下行链路多用户MIMO传输场景,其中 K < M K < M K<M 单天线用户观察频率平坦的高斯噪声信道的输出。在时间实例 n n n 处的 K K K 个用户的信道输出,通过 K × 1 K \times 1 K×1 向量 r [ n ] {\boldsymbol r}[n] r[n] 表示,由
r [ n ] = H W c c [ n ] + H W r s [ n ] + v [ n ] , \begin{equation*} \boldsymbol {r} [n] = \boldsymbol {H} \boldsymbol {W}_c \boldsymbol {c} [n] + \boldsymbol {H} \boldsymbol {W}_r \boldsymbol {s} [n] + \boldsymbol {v} [n], \tag{16} \end{equation*} r[n]=HWcc[n]+HWrs[n]+v[n],(16)
式中, H \boldsymbol {H} H 为 K × M K\times M K×M 窄带信道矩阵(narrow-band channel matrix), v [ n ] \boldsymbol {v}[n] v[n] 为协方差为 σ 2 I K \sigma ^2 \boldsymbol {I}_K σ2IK 的加性高斯噪声(AWGN)。
在多用户发射波束形成中,预编码器的设计应保证用户处有一定的SINR。这里,假设发射阵列知道瞬时下行信道 H \boldsymbol {H} H,例如,可以通过在时分双工模式下工作时利用无线信道互易性获得,即下行信道是通过上行信道估计获得的。或者,在频分双工模式下,下行信道可以通过用户的信道反馈获得,参见,例如[56]。将等效雷达到用户信道和等效用户间信道矩阵定义为
F r = H W r , F c = H W c , \begin{align*} \boldsymbol {F}_r = \boldsymbol {H} \boldsymbol {W}_r,\tag{17a} \\ \boldsymbol {F}_c = \boldsymbol {H} \boldsymbol {W}_c, \tag{17b} \end{align*} Fr=HWr,Fc=HWc,(17a)(17b)
分别。由于用户通常不能相互合作,因此 F c \boldsymbol {F}_c Fc 的非对角线元素会导致用户间干扰,应该通过预编码来缓解。同时,由于用户一般没有任何关于雷达波形的先验信息, F r \boldsymbol {F}_r Fr 会导致来自雷达的干扰。在第 k k k 个用户处,信号功率为
E ( ∣ [ F c ] k , k c k ( t ) ∣ 2 ) = ∣ [ F c ] k , k ∣ 2 , \begin{equation*} {\textbf E} \left(| [ \boldsymbol {F}_c ]_{k,k} c_k(t) | ^ 2 \right) = \left| [ \boldsymbol {F}_c ]_{k,k} \right| ^ 2, \tag{18} \end{equation*} E(∣[Fc]k,kck(t)∣2)=∣[Fc]k,k∣2,(18)
用户间干扰功率为
E ( ∣ [ F c ] k , k c k ( t ) ∣ 2 ) = ∣ [ F c ] k , k ∣ 2 , \begin{equation*} {\textbf E} \left(| [ \boldsymbol {F}_c ]_{k,k} c_k(t) | ^ 2 \right) = \left| [ \boldsymbol {F}_c ]_{k,k} \right| ^ 2, \tag{18} \end{equation*} E(∣[Fc]k,kck(t)∣2)=∣[Fc]k,k∣2,(18)
雷达的干扰强度是
E ( ∑ i = 1 M ∣ [ F r ] k , i s i ( t ) ∣ 2 ) = ∑ i = 1 M ∣ [ F r ] k , i ∣ 2 . \begin{equation*} {\textbf E} \left(\sum _{i = 1}^{M} | [ \boldsymbol {F}_r ]_{k,i} s_i(t) | ^ 2 \right) = \sum _{i = 1}^{M} \left| [ \boldsymbol {F}_r ]_{k,i} \right|^2. \tag{20} \end{equation*} E(i=1∑M∣[Fr]k,isi(t)∣2)=i=1∑M∣[Fr]k,i∣2.(20)
因此,第 k k k 个用户处的SINR表示为
γ k = ∣ [ F c ] k , k ∣ 2 ∑ i ≠ k ∣ [ F c ] k , i ∣ 2 + ∑ i = 1 M ∣ [ F r ] k , i ∣ 2 + σ 2 . \begin{equation*} \gamma _k = \frac{| [ \boldsymbol {F}_c ]_{k,k} | ^ 2 }{\sum _{i \ne k} | [ \boldsymbol {F}_c ]_{k,i} |^2 + \sum _{i = 1}^{M} | [ \boldsymbol {F}_r ]_{k,i} |^2 + \sigma ^ 2 }. \tag{21} \end{equation*} γk=∑i=k∣[Fc]k,i∣2+∑i=1M∣[Fr]k,i∣2+σ2∣[Fc]k,k∣2.(21)
多用户波束成形的两个典型设计准则是[57]、[58]:
- 吞吐量:最大化和速率
C ( γ ) = ∑ k = 1 K log 2 ( 1 + γ k ) . \begin{equation*} C(\boldsymbol {\gamma }) = \sum _{k=1}^{K} \log _2 \left(1 + \gamma _k \right). \tag{22} \end{equation*} C(γ)=k=1∑Klog2(1+γk).(22) - 公平性:最大化最小SINR,指的是
F ( γ ) = min { γ 1 , … , γ K } . \begin{equation*} F(\boldsymbol {\gamma }) = \mathop{\text{min}} \lbrace \gamma _1, \ldots, \gamma _K \rbrace . \tag{23} \end{equation*} F(γ)=min{γ1,…,γK}.(23)
这里,
γ
=
[
γ
1
,
…
,
γ
K
]
T
\boldsymbol {\gamma } = [\gamma _1, \ldots, \gamma _K] ^ T
γ=[γ1,…,γK]T。在这项工作中,我们使用公平性SINR
F
(
γ
)
F(\boldsymbol {\gamma })
F(γ) 作为多用户通信的性能指标,并要求它高于给定的阈值
Γ
\Gamma
Γ,以保证每个用户的通信服务质量的最低水平,即
γ
k
≥
Γ
,
k
=
1
,
…
,
K
.
\begin{equation*} \gamma _ k \geq \Gamma, \ k = 1,\ldots,K. \tag{24} \end{equation*}
γk≥Γ, k=1,…,K.(24)
公平性波束成形在计算复杂度上更简单,可以在多项式时间内求解,而最优吞吐量波束成形是NP难的[57]。我们注意到,在Section IV公式化的联合波束形成问题中,公平性SINR要求可以扩展为具有K个单独的SINR约束[59],即 γ k ≥ Γ k \gamma _ k \geq \Gamma _k γk≥Γk 对于 k = 1 , … , K k = 1,\ldots,K k=1,…,K,其中 Γ k \Gamma_k Γk 为第K个用户处的SINR阈值。
由于目标回波(echoes)具有双向传播损耗,雷达信号的发射功率通常远高于典型的通信发射机。在DRFC方案中,高功率雷达干扰会导致通信接收机恢复发射符号的能力显著下降。因此,提出了SINR约束(24)来缓解雷达对用户的干扰。当通信用户和目标的角度方向足够明显时,雷达波形可以对目标形成高度定向的波束,同时通过独立通信和雷达波形的发射波束形成方法消除通信用户处的雷达干扰。当只传输如图2(b)所示的通信波形时,这种波束形成能力可能无法实现,这促使使用遵循图2(a)所示模型的单独预编码通信和雷达波形。
然而,当用户和雷达目标位于相似的角度方向时,取消雷达干扰可能导致无法正确照射雷达目标。然而,DFRC系统完全了解通信波形,因此其反射信号也可以用于目标检测,这表明可以通过形成具有通信波形的发射波束同时覆盖目标和用户进行联合雷达传感和数据传输来克服这一挑战。这可以通过仅传输通信波形来实现感知和通信,如图2(b)所示,这是图2(a)中联合波束形成模型的一个特殊情况。因此,基于单独预编码的雷达和通信波形的联合波束形成模型,预计将适用于广泛的场景,包括图2(b)中的模型适用的场景,以及其他具有挑战性的设置,如第v节中数值演示的那样。为了实现这种改进的、鲁棒的双重功能能力,在下一节中,我们根据图2(b)的信号模型推导出联合波束形成方案。
IV.Joint Transmit Beamforming
根据提出的MIMO雷达和通信性能指标,我们现在转向设计DFRC联合波束形成方案。我们首先将联合发射波束成形描述为关于第IV-A节中的预编码矩阵的优化问题。为了解决这个问题,我们在第IV-B节中提出了一种基于半定松弛(SDR)的优化方案,以及在第IV-C节中提出了一种消除用户间干扰和雷达干扰的迫零(ZF)方法。
A. Problem Formulation
我们联合DFRC波束形成的目标是在发射功率和通信服务质量约束下优化雷达波束方向。特别是,在单天线功率约束(7)和每个下行链路用户公平SINR约束(24)下,我们最小化了(14)中定义的雷达波束方向图上的损失函数。
设 W = [ W c , W r ] \boldsymbol {W} = [\boldsymbol {W}_c, \boldsymbol {W}_r] W=[Wc,Wr] 为整体预编码矩阵。预编码矩阵可以通过求解下面的优化问题得到
min W , α L r ( R , α ) subject to R = W W H ∈ S M + , [ R ] m , m = P t / M , m = 1 , … , M , γ k ≥ Γ , k = 1 , … , K , \begin{align*} \mathop{\text{min}} _{\boldsymbol {W}, \alpha } \ & \ L_r(\boldsymbol {R}, \alpha) \tag{25a} \\ \text{subject to} \ & \ \boldsymbol {R} = \boldsymbol {W}\boldsymbol {W}^ H \in \mathcal {S}_{M}^+, \tag{25b} \\ & \ \left[\boldsymbol {R}\right]_{m,m} = P_t / M, \ m = 1,\ldots,M, \tag{25c} \\ & \ \gamma _ k \geq \Gamma, \ k = 1,\ldots, K, \tag{25d} \end{align*} minW,α subject to Lr(R,α) R=WWH∈SM+, [R]m,m=Pt/M, m=1,…,M, γk≥Γ, k=1,…,K,(25a)(25b)(25c)(25d)
其中(25a) - (25c)来自(15)处理雷达性能,(25d)来自考虑公平SINR要求(24)。
阈值的选择 Γ \Gamma Γ 影响通信质量和雷达性能之间的权衡。什么时候 Γ = 0 \Gamma = 0 Γ=0, (25d)成立,雷达-通信联合波束形成问题(25)简化为只针对雷达的优化问题(15)。什么时候 Γ > 0 \Gamma > 0 Γ>0,与(15)中仅雷达发射波束形成问题相比,(17)中指定等效信道的预编码器 W \boldsymbol {W} W 受到(25d)中SINR约束的限制。因此,与仅使用雷达的情况相比,可能存在由于需要满足通信性能保证而导致的固有雷达性能损失。如果更高 Γ \Gamma Γ 设置,则预计在用户端观测到更高的信号功率和更少的干扰,进一步限制预编码矩阵。因此,MIMO雷达的性能损失变得更加显著。
由于(25b)中的二次等式约束,优化问题(25)是非凸的,因此难以求解。尽管如此,我们在第IV-B节中表明,它可以使用SDR进行重制,使得可解松弛问题的解也是原始非凸问题的全局优化(25),即松弛是紧的。为了进一步降低计算复杂度,我们在 Section IV-C 中提出了一种次优的迫零波束成形策略,该策略被证明能够接近第五节中提出的数值研究中(25)的全局解的性能。
B. Joint Transmit Beamforming via SDR
在本小节中,我们使用SDR策略[60],[61]来解决非凸问题(25)。为了实现这个目标,我们首先显式地将关系(25b)写成关于 W \boldsymbol {W} W 的每一列的二次约束,让 w i \boldsymbol {w} _ i wi 表示 W \boldsymbol {W} W 的第 i i i 列,对于 i = 1 , … , M + K i = 1, \ldots, M + K i=1,…,M+K,那么(25b)变成
R = ∑ i = 1 M + K w i w i H . \begin{equation*} \boldsymbol {R} = \sum _{i=1}^{M+K} \boldsymbol {w} _ i \boldsymbol {w} _ i ^ H. \tag{26} \end{equation*} R=i=1∑M+KwiwiH.(26)
定义 R i = w i w i H \boldsymbol {R}_i = \boldsymbol {w}_i \boldsymbol {w}_i ^ H Ri=wiwiH,我们有
R = ∑ i = 1 M + K R i , \begin{equation*} \boldsymbol {R} = \sum _{i=1}^{M+K} \boldsymbol {R} _ i, \tag{27} \end{equation*} R=i=1∑M+KRi,(27)
我们省略了秩为一约束。(25d)中的SINR约束可以转换为关于秩为1的矩阵 { R i } \lbrace \boldsymbol {R}_{i} \rbrace {Ri} 的线性约束。让 h k H \boldsymbol {h} _ {k} ^ H hkH 表示 H \boldsymbol {H} H 的第 k k k 行,对于 k = 1 , … , k k = 1,\ldots, k k=1,…,k。则等效信道矩阵的项可写成 [ F r ] k , i = h k H w i + K [ \boldsymbol {F} _ r]_{k,i} = \boldsymbol {h} _ k ^ H \boldsymbol {w} _ {i+K} [Fr]k,i=hkHwi+K 和 [ F c ] k , i = h k H w i [ \boldsymbol {F} _ c]_{k,i} = \boldsymbol {h} _ k ^ H \boldsymbol {w} _ i [Fc]k,i=hkHwi。因此,SINR约束变成
γ k = h k H w k w k H h k ∑ 1 ≤ i ≤ M + K , i ≠ k h k H w i w i H h k + σ 2 = h k H R k h k ∑ 1 ≤ i ≤ M + K , i ≠ k h k H R i h k + σ 2 = h k H R k h k h k H R h k − h k H R k h k + σ 2 ≥ Γ . \begin{align*} \gamma _k &= \frac{ \boldsymbol h _ k ^ H \boldsymbol w _ k \boldsymbol w _ k ^ H \boldsymbol h _ k }{ \sum _{1\leq i \leq M+K, i \ne k} \boldsymbol h _ k ^ H \boldsymbol w _ i \boldsymbol w _ i ^ H \boldsymbol h _ k + \sigma ^2} \\ &= \frac{ \boldsymbol h _ k ^ H \boldsymbol R _ k \boldsymbol h _ k }{ \sum _{1\leq i \leq M+K, i \ne k} \boldsymbol h _ k ^ H \boldsymbol R _ i \boldsymbol h _ k + \sigma ^2} \\ &= \frac{ \boldsymbol h _ k ^ H \boldsymbol R _ k \boldsymbol h _ k }{ \boldsymbol h _ k ^ H \boldsymbol R \boldsymbol h _ k - \boldsymbol h _ k ^ H \boldsymbol R _ k \boldsymbol h _ k + \sigma ^2} \geq \Gamma . \tag{28} \end{align*} γk=∑1≤i≤M+K,i=khkHwiwiHhk+σ2hkHwkwkHhk=∑1≤i≤M+K,i=khkHRihk+σ2hkHRkhk=hkHRhk−hkHRkhk+σ2hkHRkhk≥Γ.(28)
我们现在将(25)转换为具有秩一约束的等价的二次半定规划(QSDP)
min R , { R i } , α L r ( R , α ) subject to R = ∑ i = 1 M + K R i ∈ S M + , [ R ] m , m = P t / M , m = 1 , … , M , R i ∈ S M + , r a n k ( R i ) = 1 , i = 1 , … , K + M , ( 1 + Γ − 1 ) h k H R k h k ≥ h k H R h k + σ 2 , k = 1 , … , K , \begin{align*} \mathop{\text{min}} _{ \boldsymbol {R},\lbrace \boldsymbol {R} _i \rbrace, \alpha } \ & \ L_r(\boldsymbol {R}, \alpha) \tag{29a} \\ \text{subject to} \ & \ \boldsymbol {R} = \sum _{i=1}^{M+K} \boldsymbol R _ i \in \mathcal {S}_{M}^+, \tag{29b} \\ & \ \left[\boldsymbol {R}\right]_{m,m} = P_t / M, \ m = 1,\ldots,M, \tag{29c} \\ & \ \boldsymbol {R}_i \in \mathcal {S}_{M}^+, \ \mathrm{rank} (\boldsymbol R _ i) = 1, \ i = 1,\ldots,K+M, \tag{29d} \\ \ \left(1 + \Gamma ^{-1} \right) &\boldsymbol h _ k ^ H \boldsymbol R _ k \boldsymbol h _ k \geq \boldsymbol h _ k ^ H \boldsymbol R \boldsymbol h _ k + \sigma ^2, \ k = 1,\ldots,K, \tag{29e} \end{align*} minR,{Ri},α subject to (1+Γ−1) Lr(R,α) R=i=1∑M+KRi∈SM+, [R]m,m=Pt/M, m=1,…,M, Ri∈SM+, rank(Ri)=1, i=1,…,K+M,hkHRkhk≥hkHRhk+σ2, k=1,…,K,(29a)(29b)(29c)(29d)(29e)
其中(29e)来源于(28)。我们观察到,在问题(29)中,独立的矩阵 { R i } i ≥ K + 1 \lbrace \boldsymbol R _ i \rbrace _{i \geq K+1} {Ri}i≥K+1 对SINR约束没有影响,只是封装在整体协方差矩阵R中。因此,我们可以去掉变量 { R i } i ≥ K + 1 \lbrace \boldsymbol R _ i \rbrace _{i \geq K+1} {Ri}i≥K+1 将(29)转换为
min R , R 1 , … , R K , α L r ( R , α ) subject to R ∈ S M + , R − ∑ k = 1 K R k ∈ S M + , [ R ] m , m = P t / M , m = 1 , … , M , R k ∈ S M + , r a n k ( R k ) = 1 , k = 1 , … , K , ( 1 + Γ − 1 ) h k H R k h k ≥ h k H R h k + σ 2 , k = 1 , … , K . \begin{align*} \mathop{\text{min}} _{ \boldsymbol {R}, \boldsymbol {R}_1, \ldots, \boldsymbol {R}_K, \alpha } \ & \ L_r(\boldsymbol {R}, \alpha) \tag{30a} \\ \text{subject to} \ & \ \boldsymbol {R} \in \mathcal {S}_{M}^+, \ \boldsymbol {R} - \sum _{k=1}^{K} \boldsymbol R _ k \in \mathcal {S}_{M}^+, \tag{30b} \\ & \ \left[\boldsymbol {R}\right]_{m,m} = P_t / M, \ m = 1,\ldots,M, \tag{30c} \\ & \ \boldsymbol {R}_k \in \mathcal {S}_{M}^+, \ \mathrm{rank} (\boldsymbol R _ k) = 1, \ k = 1,\ldots,K, \tag{30d} \\ \ \left(1 + \Gamma ^{-1} \right) &\boldsymbol h _ k ^ H \boldsymbol R _ k \boldsymbol h _ k \geq \boldsymbol h _ k ^ H \boldsymbol R \boldsymbol h _ k + \sigma ^2, \ k = 1,\ldots,K. \tag{30e} \end{align*} minR,R1,…,RK,α subject to (1+Γ−1) Lr(R,α) R∈SM+, R−k=1∑KRk∈SM+, [R]m,m=Pt/M, m=1,…,M, Rk∈SM+, rank(Rk)=1, k=1,…,K,hkHRkhk≥hkHRhk+σ2, k=1,…,K.(30a)(30b)(30c)(30d)(30e)
由于秩为1的约束,优化问题(30)仍然是非凸的。省略这些约束会导致以下松弛:
min R , R 1 , … , R K , α L r ( R , α ) subject to R ∈ S M + , R − ∑ k = 1 K R k ∈ S M + , [ R ] m , m = P t / M , m = 1 , … , M , R k ∈ S M + , k = 1 , … , K , ( 1 + Γ − 1 ) h k H R k h k ≥ h k H R h k + σ 2 , k = 1 , … , K . \begin{align*} \mathop{\text{min}} _{ \boldsymbol {R}, \boldsymbol {R}_1, \ldots, \boldsymbol {R}_K, \alpha } \ & \ L_r(\boldsymbol {R}, \alpha) \tag{31a} \\ \text{subject to} \ & \ \boldsymbol {R} \in \mathcal {S}_{M}^+, \ \boldsymbol {R} - \sum _{k=1}^{K} \boldsymbol R _ k \in \mathcal {S}_{M}^+, \tag{31b} \\ & \ \left[\boldsymbol {R}\right]_{m,m} = P_t / M, \ m = 1,\ldots,M, \tag{31c} \\ & \ \boldsymbol {R}_k \in \mathcal {S}_{M}^+, \ k = 1,\ldots,K, \tag{31d} \\ \ \left(1 + \Gamma ^{-1} \right) & \boldsymbol h _ k ^ H \boldsymbol R _ k \boldsymbol h _ k \geq \boldsymbol h _ k ^ H \boldsymbol R \boldsymbol h _ k + \sigma ^2, \ k = 1,\ldots,K. \tag{31e} \end{align*} minR,R1,…,RK,α subject to (1+Γ−1) Lr(R,α) R∈SM+, R−k=1∑KRk∈SM+, [R]m,m=Pt/M, m=1,…,M, Rk∈SM+, k=1,…,K,hkHRkhk≥hkHRhk+σ2, k=1,…,K.(31a)(31b)(31c)(31d)(31e)
松弛的优化模型(31)是一个凸QSQP,因为目标函数是一个半正定二次型,所有的约束要么是线性的,要么是半正定的。(31)的全局最优可以在多项式时间内通过凸优化工具箱[62]-[65]得到。
定理1:对于(31)存在全局最优,用 R ~ , R ~ 1 , … , R ~ K \tilde{\boldsymbol{R}},\tilde{ \boldsymbol {R} }_{1}, \ldots, \tilde{\boldsymbol{R}}_{K} R~,R~1,…,R~K 表示,满足
r
a
n
k
(
R
~
k
)
=
1
,
k
=
1
,
…
,
K
.
\begin{equation*} \mathrm{rank} (\tilde{ \boldsymbol {R} }_k) = 1, \ k = 1,\ldots,K. \end{equation*}
rank(R~k)=1, k=1,…,K.
我们注意到定理1只陈述了秩为1的全局最优存在。一般来说,(31)的全局最优可能不是唯一的,凸优化软件可能不会给出一个秩一解。一旦最优解 R ^ , R ^ 1 , … , R ^ K \hat{\boldsymbol {R}}, \hat{\boldsymbol {R}}_1,\ldots,\hat{\boldsymbol {R}} _ K R^,R^1,…,R^K 得到后,我们用它来寻找一个秩一最优解 R ~ 1 , … , R ~ K \tilde{\boldsymbol {R}}_1,\ldots,\tilde{\boldsymbol {R}} _ K R~1,…,R~K 以及相应的最优预编码器 w ~ 1 , … , w ~ K \tilde{\boldsymbol {w}}_1,\ldots,\tilde{\boldsymbol {w}} _ K w~1,…,w~K,如附录a所述。
首先,我们计算
R
~
1
,
…
,
R
~
K
\tilde{\boldsymbol {R}}_1,\ldots,\tilde{\boldsymbol {R}} _ K
R~1,…,R~K 和
R
~
,
w
~
1
,
…
,
w
~
K
\tilde{\boldsymbol R}, \tilde{\boldsymbol {w}}_1,\ldots,\tilde{\boldsymbol {w}} _ K
R~,w~1,…,w~K 通过
R
~
=
R
^
,
w
~
k
=
(
h
k
H
R
^
k
h
k
)
−
1
/
2
R
^
k
h
k
,
R
~
k
=
w
~
k
w
~
k
H
,
\begin{equation*} \tilde{\boldsymbol R}= \hat{\boldsymbol{R}}, \ \tilde{ \boldsymbol w }_k = \left(\boldsymbol h_k ^ H \hat{\boldsymbol{R}} _ k \boldsymbol h _ k \right) ^ {-1/2} \hat{\boldsymbol{R}} _ k \boldsymbol h _ k, \ \tilde{\boldsymbol{R}}_k = \tilde{ \boldsymbol w }_k \tilde{ \boldsymbol w }_k ^ H, \tag{32} \end{equation*}
R~=R^, w~k=(hkHR^khk)−1/2R^khk, R~k=w~kw~kH,(32)
k = 1 , … , K k = 1,\ldots,K k=1,…,K 。根据定理1的证明, R ~ , R ~ 1 , … , R ~ K \tilde{\boldsymbol {R}}, \tilde{\boldsymbol {R}}_1,\ldots,\tilde{\boldsymbol {R}} _ K R~,R~1,…,R~K 是最优的。为了证明 R ~ , R ~ 1 , … , R ~ K \tilde{\boldsymbol {R}}, \tilde{\boldsymbol {R}}_1,\ldots,\tilde{\boldsymbol {R}} _ K R~,R~1,…,R~K 对(29)也是最优的,我们构造了rank-one 矩阵 { R ~ i } i ≥ K + 1 \lbrace \tilde{\boldsymbol{R}} _ i \rbrace _{i \geq K+1} {R~i}i≥K+1as R ~ i = w ~ i w ~ i H \tilde{\boldsymbol {R}}_i= \tilde{\boldsymbol {w}} _ i\tilde{\boldsymbol {w}} _ i ^ H R~i=w~iw~iH,其中向量 w ~ i \tilde{\boldsymbol {w}}_i w~i for i > K i > K i>K 是由Cholesky分解计算的[66]
W r W r H = R ~ − ∑ k = 1 K w ~ k w ~ k H , \begin{equation*} \boldsymbol W _ r \boldsymbol W _ r ^ H = \tilde{\boldsymbol{R}} - \sum _{k=1}^{K} \tilde{ \boldsymbol w} _ k \tilde{\boldsymbol {w}} _ k ^ H, \tag{33} \end{equation*} WrWrH=R~−k=1∑Kw~kw~kH,(33)
其中 W r = [ w ~ K + 1 , … , w ~ K + M ] \boldsymbol W _ r = [\tilde{ \boldsymbol w} _ {K+1}, \ldots, \tilde{\boldsymbol w}_{K+M} ] Wr=[w~K+1,…,w~K+M] 是一个下三角矩阵。从(33)可以验证,约束(29b)适用于 R ~ , R ~ 1 , … , R ~ K + M \tilde{\boldsymbol R}, \tilde{\boldsymbol R}_1, \ldots, \tilde{\boldsymbol R}_{K+M} R~,R~1,…,R~K+M 。因此, R ~ , R ~ 1 , … , R ~ K + M \tilde{\boldsymbol R}, \tilde{\boldsymbol R}_1, \ldots, \tilde{\boldsymbol R}_{K+M} R~,R~1,…,R~K+M 是对(29)的可行解,因此也是对(29)的最优解。此外,预编码矩阵 W ~ = [ w ~ 1 , … , w ~ K + M ] \tilde{\boldsymbol W} = [ \tilde{ \boldsymbol w }_1,\ldots, \tilde{ \boldsymbol w }_{K+M} ] W~=[w~1,…,w~K+M] 是对(25)的最优解。
总结了算法1中预编码矩阵 W \boldsymbol {W} W 的计算过程。算法1的主要计算负担来自求解QSDP(31)。具体来说,给定一个解精度 ϵ \epsilon ϵ,用[67],[68]中的原始对偶内点算法求解QSDP(31)的最坏情况复杂性是 O ( K 6.5 M 6.5 log ( 1 / ϵ ) ) \mathcal {O} (K ^ {6.5} M ^ {6.5} \log (1 / \epsilon)) O(K6.5M6.5log(1/ϵ))。
C. Joint Transmit Beamforming via ZF
通过算法1获得预编码器的计算负担促使寻求降低复杂度的次优波束形成策略。在本小节中,我们重点讨论ZF波束形成。ZF方法有助于获得封闭形式、易于处理和可解释的预编码器[58],[69]。除了相对简单之外,从通信的角度来看,ZF波束形成已知会渐近地接近广播信道中的和容量[70],这表明它有可能在涉及多用户通信的设置中接近最佳性能。
我们设计了预编码器来消除通过限制Fc toa获得的用户间干扰和雷达干扰
对角矩阵和Fr为零矩阵,即,
这里,pik 为第k个用户处的信号功率,取值为1≤k≤k。强制消除干扰有助于在用户处实现高SINR值。在第V节的数值研究中,我们证明了在附加ZF约束下可实现的性能接近于(25)的全局最优解,当SINR阈值很高时,通过算法1增加计算负担获得的全局最优解。
在ZF波束形成中,SINR 约束(25d)被重新表示为 p k ≥ Γ σ 2 p_k \geq {\Gamma } \sigma ^ 2 pk≥Γσ2,对应的优化问题(25)变为
min W , α L r ( R , α ) subject to R = W W H ∈ S M + , [ R ] m , m = P t / M , m = 1 , … , M , H W = [ d i a g ( p ) , 0 K × M ] , p k ≥ Γ σ 2 , k = 1 , … , K , \begin{align*} \mathop{\text{min}} _{\boldsymbol {W}, \alpha } \ & \ L_r(\boldsymbol {R}, \alpha) \tag{35a} \\ \text{subject to} \ & \ \boldsymbol {R} = \boldsymbol {W}\boldsymbol {W}^ H \in \mathcal {S}_{M}^+, \tag{35b} \\ & \ \left[\boldsymbol {R}\right]_{m,m} = P_t / M, \ m = 1,\ldots,M, \tag{35c} \\ & \ \boldsymbol {H W} = \left[ \mathrm{diag}(\sqrt{\boldsymbol {p}}),\boldsymbol {0}_{K \times M} \right], \tag{35d} \\ & \ {p}_k \geq {\Gamma } \sigma ^ 2, \ k = 1,\ldots,K, \tag{35e} \end{align*} minW,α subject to Lr(R,α) R=WWH∈SM+, [R]m,m=Pt/M, m=1,…,M, HW=[diag(p),0K×M], pk≥Γσ2, k=1,…,K,(35a)(35b)(35c)(35d)(35e)
利用(38),全局最优的 R \boldsymbol R R 到
min R , α L r ( R , α ) subject to R ∈ S + M , H R H H = d i a g ( p ) , [ R ] m , m = P t / M , m = 1 , … , M , p k ≥ Γ σ 2 , k = 1 , … , K . \begin{align*} \mathop{\text{min}} _{\boldsymbol {R}, \alpha } \ & \ L_r(\boldsymbol {R}, \alpha) \tag{39a} \\ \text{subject to} \ & \ \boldsymbol {R} \in \mathcal {S}_+ ^ M, \ \boldsymbol {H R H}^ H = \mathrm{diag} \left(\boldsymbol {p}\right), \tag{39b} \\ & \ \left[\boldsymbol {R}\right]_{m,m} = P_t / M, \ m = 1,\ldots,M, \tag{39c} \\ & \ {p}_k \geq {\Gamma } \sigma ^ 2, \ k = 1,\ldots,K. \tag{39d} \end{align*} minR,α subject to Lr(R,α) R∈S+M, HRHH=diag(p), [R]m,m=Pt/M, m=1,…,M, pk≥Γσ2, k=1,…,K.(39a)(39b)(39c)(39d)
与(31)类似,优化(39)是一个凸QSDP,因为目标函数是正半定二次形式,所有约束要么是线性的,要么是半定的。(39)的最优值可以在多项式时间内得到。正如我们在后续文章中所示,ZF波束形成的总体复杂性大大低于通过算法1得到全局最优的复杂性。
(39)的解,即矩阵 R ~ \tilde{\boldsymbol R} R~ 和向量 p ~ \tilde{\boldsymbol p} p~,用于构造最优预编码矩阵 W ~ \tilde{\boldsymbol W} W~,如定理2的证明中所详细说明的那样。在这里,我们简要给出最终的表达式。首先,我们恢复一个 M × M M\times M M×M 矩阵,它满足 R ~ = L r L r H \tilde{\boldsymbol R}= \boldsymbol L_r \boldsymbol L_r^H R~=LrLrH。这可以使用,例如,Cholesky分解得到,尽管 L r \boldsymbol {L}_r Lr 不一定是三角形的,任何满足 R ~ = L r L r H \tilde{\boldsymbol R}= \boldsymbol L _r \boldsymbol L_r^H R~=LrLrH 的矩阵都可以用来计算 W ~ \tilde{\boldsymbol W} W~。然后,得到的预编码器 W ~ \tilde{\boldsymbol W} W~ 是
W ~ = L r Q h H [ Q f T ] 1 : M T , \begin{equation*} \tilde{\boldsymbol W} = \boldsymbol L_r \boldsymbol{Q}_h^H \left[\boldsymbol {Q}_f^T\right]^T_{{1:M}}, \tag{40} \end{equation*} W~=LrQhH[QfT]1:MT,(40)