文章目录
II. SYSTEM MODEL AND PROBLEM FORMULATION
我们考虑如图1所示的IRS辅助ISAC系统,该系统由一个具有
M
>
1
M > 1
M>1 个发射天线的BS,一个具有1个接收天线的CU,一个具有
N
>
1
N > 1
N>1 个单元的 ULA-IRS 以及BS NLoS区域的多个潜在目标组成。为了进行初步研究,假设BS完全知道全局信道状态信息(CSI)和目标的大致方向。
为了方便ISAC,我们认为BS使用发射波束形成来发送信息和专用雷达传感信号。设 s s s 为CU所需的信息信号, w ∈ C M × 1 {\bf w} \in {\mathbb{C}^{M \times 1}} w∈CM×1 为相应的发射波束形成矢量,其中设 s s s 为均值为零,单位方差为零的随机变量。设 x 0 ∈ C M × 1 {{\mathbf{x}}_0} \in {\mathbb{C}^{M \times 1}} x0∈CM×1 表示专用的传感信号,是一个随机生成的序列与s无关,具有零均值和协方差矩阵 R 0 ≜ E ( x 0 x 0 H ) ≽ 0 {{\mathbf{R}}_0} \triangleq \mathbb{E}\left( {{x_0}x_0^{\text{H}}} \right) \succcurlyeq {\mathbf{0}} R0≜E(x0x0H)≽0 。注意我们考虑了一般的多波束感知信号传输以0≤rank(R0)≤M,以便为感知提供更多的自由度。其中rank(R0)对应BS发送的感应光束数,每个感应光束可通过R0的特征值分解(EVD)得到。将信息信号与传感信号结合,BS发射的信号x为
x = w s + x 0 . \begin{equation*}{\mathbf{x}} = {\mathbf{ws}} + {{\mathbf{x}}_0}.\tag{1}\end{equation*} x=ws+x0.(1)
设P0表示BS点处的最大发射功率。我们有发射功率约束为
E ( ∥ x ∥ 2 ) = ∥ w ∥ 2 + tr ( R 0 ) ≤ P 0 . \begin{equation*}\mathbb{E}\left(\|\boldsymbol{x}\|^{2}\right)=\|\boldsymbol{w}\|^{2}+\operatorname{tr}\left(\boldsymbol{R}_{0}\right) \leq P_{0}.\tag{2}\end{equation*} E(∥x∥2)=∥w∥2+tr(R0)≤P0.(2)
首先,我们考虑从BS到CU的信息传输。设 G ∈ C N × M , h r ∈ C N × 1 {\boldsymbol{G}} \in {\mathbb{C}^{N \times M}},{{\boldsymbol{h}}_{\text{r}}} \in {\mathbb{C}^{N \times 1}} G∈CN×M,hr∈CN×1, h d ∈ C M × 1 {{\boldsymbol{h}}_{\text{d}}} \in {\mathbb{C}^{M \times 1}} hd∈CM×1 分别表示BS-IRS链路的信道矩阵,IRS-CU和BS-CU链路的信道向量。注意,通道矩阵/矢量 G \boldsymbol{G} G、 h r \boldsymbol{h}_r hr 和 h d \boldsymbol{h}_{\mathrm{d}} hd 可能包含相应链路中的LoS和NLoS信号路径。更进一步,设 ϕ n ∈ ( 0 , 2 π ] \phi_{n} \in(0,2 \pi] ϕn∈(0,2π] 表示反射元件 n ∈ { 1 , 2 , ⋯ , N } n \in\{1,2, \cdots, N\} n∈{1,2,⋯,N} 在IRS处的相移, Φ = diag ( e j ϕ 1 , … , e j ϕ N ) \boldsymbol{\Phi}=\operatorname{diag}\left(e^{j \phi_{1}}, \ldots, e^{j \phi_{N}}\right) Φ=diag(ejϕ1,…,ejϕN) 表示对应的反射波束形成矩阵。将直接BS-CU链路传输的信号与反射BS-IRS-CU链路相结合,在CU处接收到的信号为
y = ( h r H Φ G + h d H ) ( w s + x 0 ) + n (3) y=\left(\boldsymbol{h}_{\mathrm{r}}^{\mathrm{H}} \boldsymbol{\Phi} \boldsymbol{G}+\boldsymbol{h}_{\mathrm{d}}^{\mathrm{H}}\right)\left(\boldsymbol{w} s+\boldsymbol{x}_{0}\right)+n\tag{3} y=(hrHΦG+hdH)(ws+x0)+n(3)
其中 n ∼ C N ( 0 , σ 2 ) n\sim {\mathcal{C}}\mathcal{N}\left( {0,{\sigma ^2}} \right) n∼CN(0,σ2) 表示CU接收器处的加性高斯白噪声(AWGN)。需要注意的是,在CU接收端,接收信息信号 s s s 可能会受到感测信号 x 0 \boldsymbol{x}_{0} x0 的干扰。然而,由于 x 0 \boldsymbol{x}_{0} x0 可以离线生成,因此CU在传输之前就知道,因此CU可以预先取消来自 x 0 \boldsymbol{x}_{0} x0 的干扰,以便于 s s s 的解码[7]。经过这样的处理,CU接收到的SNR为
γ = ∣ ( h r H Φ G + h d H ) w ∣ 2 σ 2 (4) \gamma=\frac{\left|\left(\boldsymbol{h}_{\mathrm{r}}^{\mathrm{H}} \boldsymbol{\Phi} \boldsymbol{G}+\boldsymbol{h}_{\mathrm{d}}^{\mathrm{H}}\right) \boldsymbol{w}\right|^{2}}{\sigma^{2}}\tag{4} γ=σ2 (hrHΦG+hdH)w 2(4)
接下来,我们考虑雷达对 BS NLoS 区域潜在目标的感知。在这种情况下,我们使用由IRS反射产生的 virtual LoS links 来感知它们。因此,我们采用IRS在期望的传感角度上的波束方向增益作为传感性能指标。设 d IRS d_{\text {IRS }} dIRS 表示IRS处连续反射元件之间的间距, λ λ λ 表示波长。在IRS处,离开角(AoD) θ θ θ 的导向矢量
a ( θ ) = [ 1 , e j 2 π d I R S λ sin θ , … , e j 2 π ( N − 1 ) d I R S λ sin θ ] T \boldsymbol{a}(\theta)=\left[1, e^{j \frac{2 \pi d_{\mathrm{IRS}}}{\lambda} \sin \theta}, \ldots, e^{j \frac{2 \pi(N-1) d_{\mathrm{IRS}}}{\lambda} \sin \theta}\right]^{\mathrm{T}} a(θ)=[1,ejλ2πdIRSsinθ,…,ejλ2π(N−1)dIRSsinθ]T
( P 1 ) : max w , R 0 , Φ min l ∈ L a H ( θ l ) Φ G ( w w H + R 0 ) G H Φ H a ( θ l ) s.t. ∣ ( h r H Φ G + h d H ) w ∣ 2 σ 2 ≥ Γ , ∥ w ∥ 2 + tr ( R 0 ) ≤ P 0 , R 0 ⪰ 0 , Φ = diag ( e j ϕ 1 , … , e j ϕ N ) , (7) \begin{aligned} (\mathrm{P} 1): \max _{\boldsymbol{w}, \boldsymbol{R}_{0}, \boldsymbol{\Phi}} & \min _{l \in \mathcal{L}} \boldsymbol{a}^{\mathrm{H}}\left(\theta_{l}\right) \boldsymbol{\Phi} \boldsymbol{G}\left(\boldsymbol{w} \boldsymbol{w}^{\mathrm{H}}+\boldsymbol{R}_{0}\right) \boldsymbol{G}^{\mathrm{H}} \boldsymbol{\Phi}^{\mathrm{H}} \boldsymbol{a}\left(\theta_{l}\right) \\ \text { s.t. } & \frac{\left|\left(\boldsymbol{h}_{\mathrm{r}}^{\mathrm{H}} \boldsymbol{\Phi} \boldsymbol{G}+\boldsymbol{h}_{\mathrm{d}}^{\mathrm{H}}\right) \boldsymbol{w}\right|^{2}}{\sigma^{2}} \geq \Gamma ,\\ & \|\boldsymbol{w}\|^{2}+\operatorname{tr}\left(\boldsymbol{R}_{0}\right) \leq P_{0} ,\\ & \boldsymbol{R}_{0} \succeq 0 ,\\ & \boldsymbol{\Phi}=\operatorname{diag}\left(e^{j \phi_{1}}, \ldots, e^{j \phi_{N}}\right),\tag{7} \end{aligned} (P1):w,R0,Φmax s.t. l∈LminaH(θl)ΦG(wwH+R0)GHΦHa(θl)σ2 (hrHΦG+hdH)w 2≥Γ,∥w∥2+tr(R0)≤P0,R0⪰0,Φ=diag(ejϕ1,…,ejϕN),(7)
其中 Γ Γ Γ 表示CU的最小SNR阈值。由于发射波束形成器( w \boldsymbol{w} w 和 R 0 \boldsymbol{R}_{0} R0)和反射波束形成器( Φ \boldsymbol{\Phi} Φ)之间的耦合,问题(P1)高度非凸,因此非常难以被优化解决。
备注1:值得注意的是,对于没有IRS的传统ISAC系统,波束模式匹配(beampatten matching)是另一个广泛采用的传感设计标准(例如[7],[8],[18]),其目标是通过允许BS以全功率发射,使实现的发射波束方向图与预定波束方向图之间的均方匹配误差最小化。然而,这种设计可能不适用于我们感兴趣的 IRS-assisted ISAC 系统。这是因为在波束模式匹配设计下,BS可以选择引导能量正交于IRS,并使用直接LoS链路为CU服务,这可能导致最小(甚至为零)匹配误差,但在IRS处的波束模式增益非常小(甚至为零),这不是传感所希望的。
III. PROPOSED JOINT BEAMFORMING SOLUTION TO PROBLEM (P1)
本节提出了一种基于交替优化原理的求解问题(P1)的有效算法,其中BS处的发射波束形成器( w \boldsymbol{w} w 和 R 0 \boldsymbol{R}_{0} R0)和IRS处的反射波束形成器( Φ \boldsymbol{\Phi} Φ)通过将对方视为给定,以交替方式进行优化。对所提出的联合波束形成算法的细节进行了总结在算法1中。下面,我们将分别关注发射和反射波束形成的优化。
A. Transmit Beamforming Optimization at BS
首先,我们在任何给定的反射波束形成器 Φ \boldsymbol{\Phi} Φ 下优化问题(P1)中的发射波束形成器 w \boldsymbol w w 和 R 0 \boldsymbol{R}_0 R0。这个问题被表述为
( P 2 ) : max w , R 0 min l ∈ L a H ( θ l ) Φ G ( w w H + R 0 ) G H Φ H a ( θ l ) s.t. (7a), (7b), and (7c). \begin{align*} (\mathrm{P} 2): &\max _{\boldsymbol{w}, \boldsymbol{R}_{0}} \min _{l \in \mathcal{L}} \boldsymbol{a}^{\mathrm{H}}\left(\theta_{l}\right) \boldsymbol{\Phi} \boldsymbol{G}\left(\boldsymbol{w} \boldsymbol{w}^{\mathrm{H}}+\boldsymbol{R}_{0}\right) \boldsymbol{G}^{\mathrm{H}} \boldsymbol{\Phi}^{\mathrm{H}} \boldsymbol{a}\left(\theta_{l}\right) \\ &\text { s.t. } \quad \text { (7a), (7b), and (7c). } \end{align*} (P2):w,R0maxl∈LminaH(θl)ΦG(wwH+R0)GHΦHa(θl) s.t. (7a), (7b), and (7c).
为此,我们定义 W = w w H \boldsymbol{W}=\boldsymbol{w} \boldsymbol{w}^{\mathrm{H}} W=wwH,其中 W ⪰ 0 \boldsymbol{W} \succeq \mathbf{0} W⪰0,rank( W \boldsymbol{W} W)≤1。令 h = G H Φ H h r + h d \boldsymbol{h}=\boldsymbol{G}^{\mathrm{H}} \boldsymbol{\Phi}^{\mathrm{H}} \boldsymbol{h}_{\mathrm{r}}+\boldsymbol{h}_{\mathrm{d}} h=GHΦHhr+hd 表示从BS到CU的组合信道矢量,代之以 W = w w H \boldsymbol{W}=\boldsymbol{w} \boldsymbol{w}^{\mathrm{H}} W=wwH,将问题(P2)中的发射波束形成优化重新表述为
( P 2.1 ) : max W , R 0 min l ∈ L a H ( θ l ) Φ G ( W + R 0 ) G H Φ H a ( θ l ) s.t. tr ( h h H W ) ≥ Γ σ 2 tr ( W + R 0 ) ≤ P 0 R 0 ⪰ 0 , W ⪰ 0 rank ( W ) ≤ 1. \begin{align*} (\mathrm{P} 2.1): &\max _{\boldsymbol{W}, \boldsymbol{R}_{0}} & &\min _{l \in \mathcal{L}} \boldsymbol{a}^{\mathrm{H}}\left(\theta_{l}\right) \boldsymbol{\Phi} \boldsymbol{G}\left(\boldsymbol{W}+\boldsymbol{R}_{0}\right) \boldsymbol{G}^{\mathrm{H}} \boldsymbol{\Phi}^{\mathrm{H}} \boldsymbol{a}\left(\theta_{l}\right) \\ &\text { s.t. } & & \operatorname{tr}\left(\boldsymbol{h} \boldsymbol{h}^{\mathrm{H}} \boldsymbol{W}\right) \geq \Gamma \sigma^{2}\tag{8a} \\ & & & \operatorname{tr}\left(\boldsymbol{W}+\boldsymbol{R}_{0}\right) \leq P_{0}\tag{8b} \\ & & & \boldsymbol{R}_{0} \succeq 0, \boldsymbol{W} \succeq 0\tag{8c} \\ & & & \operatorname{rank}(\boldsymbol{W}) \leq 1 .\tag{8d} \end{align*} (P2.1):W,R0max s.t. l∈LminaH(θl)ΦG(W+R0)GHΦHa(θl)tr(hhHW)≥Γσ2tr(W+R0)≤P0R0⪰0,W⪰0rank(W)≤1.(8a)(8b)(8c)(8d)
然而,问题(P2.1)是非凸的,因为 W \boldsymbol{W} W 在(8d)上的秩一约束。为了解决这个问题,我们放宽了rank-one constraint ,得到问题的SDR版本(P2.1)为
( S D R 2.1 ) : max W , R 0 min l ∈ L a H ( θ l ) Φ G ( W + R 0 ) G H Φ H a ( θ l ) s.t. ( 8 a ) , ( 8 b ) , and (8c). \begin{aligned} (\mathrm{SDR} 2.1): &\max _{\boldsymbol{W}, \boldsymbol{R}_{0}} && \min _{l \in \mathcal{L}} \boldsymbol{a}^{\mathrm{H}}\left(\theta_{l}\right) \boldsymbol{\Phi} \boldsymbol{G}\left(\boldsymbol{W}+\boldsymbol{R}_{0}\right) \boldsymbol{G}^{\mathrm{H}} \boldsymbol{\Phi}^{\mathrm{H}} \boldsymbol{a}\left(\theta_{l}\right) \\ &\text { s.t. } && (8 \mathrm{a}),(8 \mathrm{~b}), \text { and (8c). } \end{aligned} (SDR2.1):W,R0max s.t. l∈LminaH(θl)ΦG(W+R0)GHΦHa(θl)(8a),(8 b), and (8c).
可以观察到,问题(SDR2.1)是一个半确定规划(SDP),可以通过CVX等凸求解器进行最优求解[19]。设 W ∗ \boldsymbol{W}^{*} W∗ 和 R 0 ∗ \boldsymbol{R}_0^{*} R0∗ denote为问题(SDR2.1)的最优解,其中 W ∗ \boldsymbol{W}^{*} W∗ 通常是高阶的。基于 W ∗ \boldsymbol{W}^{*} W∗,我们可以重构一个等价的秩一解,并据此找到问题(P2.1)的最优解(从而找到问题(P2)),如下面的命题所示。
Proposition 1. The optimal solution to problem (P2.1) is
W
^
=
w
^
w
^
H
,
R
^
0
=
R
0
∗
+
W
∗
−
W
^
,
\begin{align*} & \hat {\boldsymbol{W}} = \hat {\boldsymbol{w}}{ \hat{\boldsymbol{w}}^{\text{H}}},\tag{9} \\ & { \hat{\boldsymbol{R}}_0} = {\boldsymbol{R}}_0^{\ast} + {{\boldsymbol{W}}^{\ast}} - {\boldsymbol{\hat W}},\tag{10}\end{align*}
W^=w^w^H,R^0=R0∗+W∗−W^,(9)(10)
式中, w ^ = ( h H W ∗ h ) − 1 / 2 W ∗ h \hat{\boldsymbol{w}}=\left(\boldsymbol{h}^{\mathrm{H}} \boldsymbol{W}^{*} \boldsymbol{h}\right)^{-1 / 2} \boldsymbol{W}^{*} \boldsymbol{h} w^=(hHW∗h)−1/2W∗h 表示在BS处对应的发射波束形成矢量。相应地, w ^ \hat{\boldsymbol{w}} w^ 和 R ^ 0 \hat{\boldsymbol{R}}_0 R^0 成为问题(P2)的最优解。
注2:值得注意的是,对于(P2.1)的最优解,我们通常有 rank ( R ^ 0 ) ≥ 1 ( R ^ 0 ≠ 0 ) \operatorname{rank}\left(\hat{\boldsymbol{R}}_{0}\right) \geq 1\left(\hat{\boldsymbol{R}}_{0} \neq \mathbf{0}\right) rank(R^0)≥1(R^0=0),特别是当 G \boldsymbol G G 是随机生成的(例如,遵从瑞利衰落之后),这将在第四节的数值结果中显示。这表明使用专门的传感波束来增强IRS ISAC性能的必要性。这与[7]中没有IRS的情况不同,在那里可能不需要专用感知波束。
备注3:当BS与ULA一起部署,并且从BS到IRS的信道矩阵为LoS时,也值得讨论一种特殊情况,即
G = a ( θ I R S ) b H ( θ B S ) \boldsymbol{G}=\boldsymbol{a}\left(\theta_{\mathrm{IRS}}\right) \boldsymbol{b}^{\mathrm{H}}\left(\theta_{\mathrm{BS}}\right) G=a(θIRS)bH(θBS)
其中 θ I R S \theta_{\mathrm{IRS}} θIRS 和 θ B S \theta_{\mathrm{BS}} θBS 分别表示BS-IRS链路在IRS和BS处的到达角(AoA)和离开角(AoD), a ( θ I R S ) \boldsymbol{a}\left(\theta_{\mathrm{IRS}}\right) a(θIRS) 表示式(5)中IRS处的导向矢量, b ( θ B S ) \boldsymbol{b}\left(\theta_{\mathrm{BS}}\right) b(θBS) 表示BS处的导向矢量,即
b ( θ B S ) = [ 1 , e j 2 π d B S λ sin θ B S , … , e j 2 π ( M − 1 ) d B S λ sin θ B S ] T , \boldsymbol{b}\left(\theta_{\mathrm{BS}}\right)=\left[1, e^{j \frac{2 \pi d_{\mathrm{BS}}}{\lambda} \sin \theta_{\mathrm{BS}}}, \ldots, e^{j \frac{2 \pi(M-1) d_{\mathrm{BS}}}{\lambda} \sin \theta_{\mathrm{BS}}}\right]^{\mathrm{T}}, b(θBS)=[1,ejλ2πdBSsinθBS,…,ejλ2π(M−1)dBSsinθBS]T,
在这种情况下,可以类似于[7]表明, R 0 ∗ = 0 \boldsymbol{R}_{0}^{*}=\mathbf{0} R0∗=0 是问题(P2.1)的一般最优,这意味着在这种特殊情况下不需要 dedicated sensing beams。
B. Reflective Beamforming Optimization at IRS
接下来,我们在任意给定的主动发射波束形成器 w \boldsymbol{w} w 和 R 0 \boldsymbol{R}_{0} R0 下优化反射波束形成器 Φ \boldsymbol{\Phi} Φ 问题(P1)。设 v = [ e j ϕ 1 , … , e j ϕ N ] H \boldsymbol{v}=\left[e^{j \phi_{1}}, \ldots, e^{j \phi_{N}}\right]^{\mathrm{H}} v=[ejϕ1,…,ejϕN]H 表示IRS处的反射相移矢量。然后,从IRS到角度 θ θ θ 的传感波束图增益重写为
P ( θ ) = v H R 1 ( θ ) v (14) \mathcal{P}(\theta)=\boldsymbol{v}^{\mathrm{H}} \boldsymbol{R}_{1}(\theta) \boldsymbol{v}\tag{14} P(θ)=vHR1(θ)v(14)
与
R
1
(
θ
)
=
diag
(
a
H
(
θ
)
)
G
(
R
0
+
W
)
G
H
diag
(
a
(
θ
)
)
\boldsymbol{R}_{1}(\theta)=\operatorname{diag}\left(\boldsymbol{a}^{\mathrm{H}}(\theta)\right) \boldsymbol{G}\left(\boldsymbol{R}_{0}+\boldsymbol{W}\right) \boldsymbol{G}^{\mathrm{H}} \operatorname{diag}(\boldsymbol{a}(\theta))
R1(θ)=diag(aH(θ))G(R0+W)GHdiag(a(θ))
更进一步,定义
R 2 ( θ ) = [ R 1 ( θ ) 0 N × 1 0 1 × N 0 ] , v ‾ = [ v 1 ] \boldsymbol{R}_{2}(\theta)=\left[\begin{array}{cc} \boldsymbol{R}_{1}(\theta) & \mathbf{0}_{N \times 1} \\ \mathbf{0}_{1 \times N} & 0 \end{array}\right], \overline{\boldsymbol{v}}=\left[\begin{array}{l} \boldsymbol{v} \\ 1 \end{array}\right] R2(θ)=[R1(θ)01×N0N×10],v=[v1]
将(16)代入(14),我们得到
P
(
θ
)
=
v
‾
H
R
2
(
θ
)
v
‾
\mathcal{P}(\theta)=\overline{\boldsymbol{v}}^{\mathrm{H}} \boldsymbol{R}_{2}(\theta) \overline{\boldsymbol{v}}
P(θ)=vHR2(θ)v。进一步,让
H
=
diag
(
h
r
H
)
G
∈
C
N
×
M
\boldsymbol{H}=\operatorname{diag}\left(\boldsymbol{h}_{\mathrm{r}}^{\mathrm{H}}\right) \boldsymbol{G} \in \mathbb{C}^{N \times M}
H=diag(hrH)G∈CN×M,(4)中接收到的信号功率重写为
∣
(
h
r
H
Φ
G
+
h
d
H
)
w
∣
2
=
∣
(
v
H
H
+
h
d
H
)
w
∣
2
\left|\left(\boldsymbol{h}_{\mathrm{r}}^{\mathrm{H}} \boldsymbol{\Phi} \boldsymbol{G}+\boldsymbol{h}_{\mathrm{d}}^{\mathrm{H}}\right) \boldsymbol{w}\right|^{2}=\left|\left(\boldsymbol{v}^{\mathrm{H}} \boldsymbol{H}+\boldsymbol{h}_{\mathrm{d}}^{\mathrm{H}}\right) \boldsymbol{w}\right|^{2}
(hrHΦG+hdH)w
2=
(vHH+hdH)w
2。(7a)中的信噪比约束表示为
(
v
H
H
+
h
d
H
)
W
(
H
H
v
+
h
d
)
≥
Γ
σ
2
\left(\boldsymbol{v}^{\mathrm{H}} \boldsymbol{H}+\boldsymbol{h}_{\mathrm{d}}^{\mathrm{H}}\right) \boldsymbol{W}\left(\boldsymbol{H}^{\mathrm{H}} \boldsymbol{v}+\boldsymbol{h}_{\mathrm{d}}\right) \geq \Gamma \sigma^{2}
(vHH+hdH)W(HHv+hd)≥Γσ2
等价于
v
‾
H
R
3
v
‾
+
h
d
H
W
h
d
≥
Γ
σ
2
\overline{\boldsymbol{v}}^{\mathrm{H}} \boldsymbol{R}_{3} \overline{\boldsymbol{v}}+\boldsymbol{h}_{\mathrm{d}}^{\mathrm{H}} \boldsymbol{W} \boldsymbol{h}_{\mathrm{d}} \geq \Gamma \sigma^{2}
vHR3v+hdHWhd≥Γσ2
其中
R
3
=
[
H
W
H
H
H
W
h
d
h
d
H
W
H
H
0
]
.
\boldsymbol{R}_{3}=\left[\begin{array}{cc} \boldsymbol{H} \boldsymbol{W} \boldsymbol{H}^{\mathrm{H}} & \boldsymbol{H} \boldsymbol{W} \boldsymbol{h}_{\mathrm{d}} \\ \boldsymbol{h}_{\mathrm{d}}^{\mathrm{H}} \boldsymbol{W} \boldsymbol{H}^{\mathrm{H}} & 0 \end{array}\right] .
R3=[HWHHhdHWHHHWhd0].
因此,
Φ
\boldsymbol{\Phi}
Φ 问题(P1)的优化就变成了
v
‾
\overline{\boldsymbol{v}}
v 的优化,这就是下面的问题:
(P3):
max
v
‾
min
l
∈
L
v
‾
H
R
2
(
θ
l
)
v
‾
s.t.
v
‾
H
R
3
v
‾
+
h
d
H
W
h
d
≥
Γ
σ
2
,
∣
v
‾
n
∣
=
1
,
∀
n
∈
{
1
,
…
,
N
+
1
}
.
\begin{aligned} \text { (P3): } \max _{\overline{\boldsymbol{v}}} & \min _{l \in \mathcal{L}} \overline{\boldsymbol{v}}^{\mathrm{H}} \boldsymbol{R}_{2}\left(\theta_{l}\right) \overline{\boldsymbol{v}} \\ \text { s.t. } & \overline{\boldsymbol{v}}^{\mathrm{H}} \boldsymbol{R}_{3} \overline{\boldsymbol{v}}+\boldsymbol{h}_{\mathrm{d}}^{\mathrm{H}} \boldsymbol{W} \boldsymbol{h}_{\mathrm{d}} \geq \Gamma \sigma^{2}, \\ & \left|\overline{\boldsymbol{v}}_{n}\right|=1, \forall n \in\{1, \ldots, N+1\} . \end{aligned}
(P3): vmax s.t. l∈LminvHR2(θl)vvHR3v+hdHWhd≥Γσ2,∣vn∣=1,∀n∈{1,…,N+1}.
接下来,我们用 V ‾ = v ‾ v ‾ H \overline{\boldsymbol{V}}=\overline{\boldsymbol{v}} \overline{\boldsymbol{v}}^{\mathrm{H}} V=vvH 和 rank ( V ‾ ) = 1 \operatorname{rank}(\overline{\boldsymbol{V}})=1 rank(V)=1。请注意 v ‾ H R 2 ( θ ) v ‾ = tr ( R 2 ( θ ) V ‾ ) \overline{\boldsymbol{v}}^{\mathrm{H}} \boldsymbol{R}_{2}(\theta) \overline{\boldsymbol{v}}=\operatorname{tr}\left(\boldsymbol{R}_{2}(\theta) \overline{\boldsymbol{V}}\right) vHR2(θ)v=tr(R2(θ)V) 和 v ‾ H R 3 v ‾ = tr ( R 3 V ‾ ) \overline{\boldsymbol{v}}^{\mathrm{H}} \boldsymbol{R}_{3} \overline{\boldsymbol{v}}=\operatorname{tr}\left(\boldsymbol{R}_{3} \overline{\boldsymbol{V}}\right) vHR3v=tr(R3V)。通过将 V ‾ = v ‾ v ‾ H \overline{\boldsymbol{V}}=\overline{\boldsymbol{v}} \overline{\boldsymbol{v}}^{\mathrm{H}} V=vvH 代入,问题(P3)中的反射波束形成优化问题被重新表述为
(P3.1): max V ‾ min l ∈ L tr ( R 2 ( θ l ) V ‾ ) s.t. tr ( R 3 V ‾ ) + h d H W h d ≥ Γ σ 2 , V ‾ n , n = 1 , ∀ n ∈ { 1 , … , N + 1 } , V ‾ ⪰ 0 , rank ( V ‾ ) = 1. \begin{aligned} \text { (P3.1): } \max _{\overline{\boldsymbol{V}}} & \min _{l \in \mathcal{L}} \operatorname{tr}\left(\boldsymbol{R}_{2}\left(\theta_{l}\right) \overline{\boldsymbol{V}}\right) \\ \text { s.t. } & \operatorname{tr}\left(\boldsymbol{R}_{3} \overline{\boldsymbol{V}}\right)+\boldsymbol{h}_{\mathrm{d}}^{\mathrm{H}} \boldsymbol{W} \boldsymbol{h}_{\mathrm{d}} \geq \Gamma \sigma^{2}, \\ & \overline{\boldsymbol{V}}_{n, n}=1, \forall n \in\{1, \ldots, N+1\}, \\ & \overline{\boldsymbol{V}} \succeq \mathbf{0}, \\ & \operatorname{rank}(\overline{\boldsymbol{V}})=1 . \end{aligned} (P3.1): Vmax s.t. l∈Lmintr(R2(θl)V)tr(R3V)+hdHWhd≥Γσ2,Vn,n=1,∀n∈{1,…,N+1},V⪰0,rank(V)=1.
然而,问题(P3.1)是非凸的,因为 V ‾ \overline{\boldsymbol{V}} V 在(21d)上的秩一约束。为了解决这个问题,我们放宽了秩一约束,得到问题的SDR版本(P3.1)为
(SDR3.1) : max V ‾ min l ∈ L tr ( R 2 ( θ l ) V ˉ ) s.t. ( 21 a ) , ( 21 b ) , and ( 21 c ) . \begin{aligned} \text { (SDR3.1) : } \max _{\overline{\boldsymbol{V}}} & \min _{l \in \mathcal{L}} \operatorname{tr}\left(\boldsymbol{R}_{2}\left(\theta_{l}\right) \bar{\boldsymbol{V}}\right) \\ \text { s.t. } & (21 \mathrm{a}),(21 \mathrm{~b}), \text { and }(21 \mathrm{c}) . \end{aligned} (SDR3.1) : Vmax s.t. l∈Lmintr(R2(θl)Vˉ)(21a),(21 b), and (21c).
问题(SDR3.1)是一个SDP,可以通过诸如CVX[19]等凸求解器进行最优求解。设 V ˉ ∗ \bar{\boldsymbol{V}}^{*} Vˉ∗ 表示得到的问题(SDR3.1)的最优解,它通常是高阶的。然后用高斯随机化(Gaussian randomization)构造一个近似的秩一解。具体来说,我们首先生成若干随机化 r ∼ C N ( 0 , V ˉ ∗ ) \boldsymbol{r} \sim \mathcal{C N}\left(\mathbf{0}, \bar{\boldsymbol{V}}^{*}\right) r∼CN(0,Vˉ∗),并相应地构造问题(P3)的候选可行解为 v = e j arg ( [ r r N + 1 ] ( 1 : N ) ) \boldsymbol{v}=e^{j \arg \left(\left[\frac{\boldsymbol{r}}{\boldsymbol{r}_{N+1}}\right]_{(1: N)}\right)} v=ejarg([rN+1r](1:N))。由独立生成高斯随机向量 r \boldsymbol{r} r 进行多次运算,将目标值近似为所有随机实现中的最大值。注意,高斯随机化应该实现多次,以确保每次迭代的目标值都增加。
将得到的(P2.1)和(P3.1)问题的解结合起来,就完成了求解(P1)问题的交替优化算法。注意,在交替优化的每次迭代中,问题(P2.1)都会得到最优解,因此会得到一个不递减的目标值。还需要注意的是,如果随机化的数量足够多,那么解决问题(P3.1)后的目标值也将是单调不递减的。因此,所提出的基于交替优化的求解问题算法(P1)的收敛性得到了保证。
IV. NUMERICAL RESULTS
本节提供数值结果来评估我们提出的 IRS-assisted ISAC设计的性能。我们考虑BS-IRS和IRS-CU链路的 Rician 衰落(Rician因子为0.5)和BS-CU链路的 Rayleigh 衰落。距离相关的路径损耗建模为 K 0 ( d d 0 ) − α K_{0}(\frac{d}{d_{0}})^{-\alpha} K0(d0d)−α,其中 K 0 = − 30 d B K_0=−30 \mathrm{~dB} K0=−30 dB为参考距离 d 0 = 1 m d_0= 1 \mathrm m d0=1m 处的路径损耗,BS-IRS、IRS-CU和BS-CU链路的路径损耗指数 α α α 分别设为2.5、2.5和3.5。
h = C 0 ( d d 0 ) α ⋅ ( β β + 1 h L o S + 1 β + 1 h N L o S ) h=\sqrt{C_0\left(\frac{d}{d_0}\right)^\alpha}\cdot\left(\sqrt{\frac{\beta}{\beta+1}}h^{LoS}+\sqrt{\frac{1}{\beta+1}}h^{NLoS}\right) h=C0(d0d)α⋅(β+1βhLoS+β+11hNLoS)这里我们根据要求可以设置 β = 0.5 \beta=0.5 β=0.5
omega = 1;
% Normalization
% To generate a Rayleigh fading channel gain
h_rayleigh = sqrt(omega/2) *(randn(1,1) + 1i*randn(1,1));
% Rician factor
k = 10;
% To generate a Rician fading channel gain
h_rician = sqrt(k/(k+1)) + sqrt(1/(k+1)) * h_rayleigh;
由于潜在的障碍,我们考虑了BS-CU链路的额外阴影衰落,标准偏差为 10 d B 10 \mathrm {~dB} 10 dB。BS、CU和IRS分别位于坐标(0,0)、(50 m, 0)和(18 m, 2 m)。所需的感测角度在 [ − 6 1 ◦ , − 5 9 ◦ ] [−61^◦,−59^◦] [−61◦,−59◦], [ − 3 1 ◦ , − 2 9 ◦ ] [−31^◦,−29^◦] [−31◦,−29◦], [ − 1 ◦ , 1 ◦ ] [−1^◦,1^◦] [−1◦,1◦], [ 2 9 ◦ , 3 1 ◦ ] [29^◦,31^◦] [29◦,31◦],和 [ 5 9 ◦ , 6 1 ◦ ] \left[59^◦,61^◦\right] [59◦,61◦]进行采样,其中采样间隔为 0.2 5 ◦ 0.25^◦ 0.25◦。我们设BS点天线数 M = 8 M = 8 M=8,IRS反射单元数 N = 64 N = 64 N=64。设 P 0 = 20 d B m P_0 = 20 \mathrm{~dBm} P0=20 dBm, σ 2 = − 80 d B m σ^2=−80 \mathrm{~dBm} σ2=−80 dBm。
图2显示了我们提出的求解问题(P1)的交替优化算法的收敛行为,其中 Γ = 10 d B Γ= 10\mathrm{~dB} Γ=10 dB。结果表明,所提出的基于交替优化的算法在大约10次迭代内收敛,从而验证了其有效性。
接下来,我们比较了我们提出的 IRS-assisted ISAC设计与以下基准方案的性能。
-
Information beamforming only:BS使用信息波束 w \boldsymbol w w 进行通信和传感。这种情况下的联合波束成形设计对应于R0= 0的问题(P1)。在该方案中,得到的问题(SDR2.1)的解通常具有高秩,并使用高斯随机化来重建近似的秩一解。
-
分离波束形成设计:该方案分别优化了BS处的发射波束形成和IRS处的反射波束形成。首先,对反射波束形成器Φ进行优化,使IRS的最小通道规范朝向所需的传感角度最大化,即:
( P4 ) : max Φ min l ∈ L E ( ∥ a H ( θ l ) Φ G ∥ 2 ) s . t . Φ = d i a g ( e j ϕ 1 , . . . , e j ϕ N ) . \begin{aligned}(\text{P4}):&\max_{\Phi}& &\min_{l\in\mathcal{L}}\mathbb{E}(\|a^{\mathrm{H}}(\theta_l)\Phi G\|^2)\\&\mathrm{s.t.}& &\Phi=\mathrm{diag}(e^{j\phi_1},...,e^{j\phi_N}).\end{aligned} (P4):Φmaxs.t.l∈LminE(∥aH(θl)ΦG∥2)Φ=diag(ejϕ1,...,ejϕN).
然后,利用问题(P4)中的反射波束形成器Φobtained,通过求解问题(P2)来设计联合传输信息和感知波束形成。
图3显示了不同方案可实现的波束方向图增益,其中Γ= 10 dB。观察到,所提出的设计产生多个信号波束从IRS到所需的传感角度。所提出的基于交替优化的设计被观察到在最差角度下获得更高的波束图增益,比其他两种基准方案。
图4显示了在传感角度下的,最小波束方向图增益与接收到的CU信噪比阈值 Γ Γ Γ 的关系。可以观察到,所提出的基于交替优化的算法在与两个基准方案对比,取得了显著的性能增益。当 Γ Γ Γ 小(例如 Γ < 6 d B \Gamma<6\mathrm{~dB} Γ<6 dB)时,
观察到单独的波束成形设计的表现接近所提出的基于交替优化的设计,并优于仅信息波束成形方案。这是因为在这种情况下,CU处的SNR可以很容易地得到满足,因此在单独的波束形成设计中,面向传感的反射波束形成(即问题(P4))变得可取。此外,可以观察到,与仅使用信息波束成形的情况相比,使用感知信号会导致显著的感知性能增强。由于考虑了随机信道g,这与备注2一致。相反,当Γis高(例如Γ> 12 dB)时,可以观察到,仅信息设计与基于交替优化的设计接近,单独波束形成设计变得不可行。这是因为在这种情况下,为了满足CU处的SNR要求,需要将大部分能量分配给信息信号,从而有利于信息波束形成的设计。