Abstract
Simultaneously transmitting and reflecting (STAR) reconfigurable intelligent surface (RIS) 比年来成为提升传统仅反射RIS性能的一项有远景的技能。思量到比较不同STAR-RIS工作模式的无线系统的复杂性，以及由离散选择引起的性能下降，本文提出了一个统一的优化框架，用于处置惩罚来自不同STAR-RIS工作模式和离散相位系数的约束。通过奥妙引入惩罚项，该框架将原题目转化为两个迭代子题目：一个处置惩罚选择型约束（selection-type constraints），另一个处置惩罚其他无线资源。在温和的条件下，整个算法的收敛点至少是一个驻点。作为示例，本文将该框架应用于下行链路传输中的和速率最大化题目。仿真结果表明，基于该框架的算法不仅优于其他为不同STAR-RIS场景设计的现有算法，还能实现不同工作模式下的公平统一比较。此外，研究还发现，使用4个乃至2个离散相位的STAR-RIS，其和速率性能几乎与一连相位设置雷同，这初次表明离散相位并不肯定会导致显著的性能下降。
通过前述的讨论（preceding deliberations）可以看出，STAR-RIS中的不同约束条件是通过多种优化技能来处置惩罚的，这使得不同范例的STAR-RIS之间的比较变得困难。此外，当多个约束条件同时出现时，如今尚不清楚哪些现有方法可以推广到这种情形。为填补这一空白，本文初次提出了一种创新的统一框架，用于处置惩罚操作模式和离散相位约束。通过引入辅助变量来表示STAR-RIS系数，我们奥妙地将原题目转化为两个不同的子题目。一个子题目专门处置惩罚涉及离散选择的约束，另一个子题目则处置惩罚总功率约束以及额外的无线资源约束。这种战略性的解耦使得能够推导出与选择约束相关的子题目的闭式全局最优解，从而有助于证明所提出框架的解质量。具体而言，收敛的解在宽松条件下包管至少是一个驻点。据我们所知，这是初次在各种STAR-RIS配置下，纵然在离散相位条件下，也提供了解质量包管的工作。
STAR-RIS有三种约束范例，即操作模式约束、无损功率约束和相位约束。以下是这些约束的具体分析。
Operating mode constraint. STAR-RIS 通常有三种操作模式 [15]：
- ES模式（Energy Splitting mode）
在这种模式下，入射信号被分为反射信号和透射信号，对应的约束条件为：
∣ v m t ∣ , ∣ v m r ∣ ∈ [ 0 , 1 ] |v_m^t|, |v_m^r| \in [0, 1] ∣vmt∣,∣vmr∣∈[0,1]
- MS模式（Mode Switching mode）
在这种模式下，每个 STAR-RIS 单元要么专用于反射，要么专用于透射，对应的约束条件为：
∣ v m t ∣ , ∣ v m r ∣ ∈ { 0 , 1 } |v_m^t|, |v_m^r| \in \{0, 1\} ∣vmt∣,∣vmr∣∈{0,1}
- TS模式（Time Switching mode）
在这种模式下，传输时间被分别为两个子区间，其中一个子区间专用于反射，另一个专用于透射。令 λ t ≥ 0 \lambda^t \geq 0 λt≥0 和 λ r ≥ 0 \lambda^r \geq 0 λr≥0 分别表示分配给透射时段和反射时段的时间比例，则有：
λ t + λ r = 1 \lambda^t + \lambda^r = 1 λt+λr=1对于 ES 和 MS 模式，由于透射和反射均占据整个时间区间，可以设置： λ t = λ r = 1 \lambda^t = \lambda^r = 1 λt=λr=1
Lossless power constraint. 通常假设超表面是无损的。因此，在 ES模式 和 MS模式 下，反射能量与透射能量之和必须即是入射信号能量，对应的约束条件为：
∣ v m t ∣ 2 + ∣ v m r ∣ 2 = 1 |v_m^t|^2 + |v_m^r|^2 = 1 ∣vmt∣2+∣vmr∣2=1
在 TS模式 下，由于在某一时间隔断内仅举行透射或反射操作，无损约束意味着：
∣ v m t ∣ = ∣ v m r ∣ = 1 |v_m^t| = |v_m^r| = 1 ∣vmt∣=∣vmr∣=1
图2展示了 STAR-RIS 在不同操作模式下的工作情况以及对应的约束条件。
Phase constraint ：
在现实应用中，STAR-RIS 的相位通常无法具有无限分辨率。在这种情况下：
∠ v m t , ∠ v m r ∈ { 0 , 2 π / L , ⋯ , 2 π ( L − 1 ) / L } \angle v_m^t, \angle v_m^r \in \{0, 2\pi/L, \cdots, 2\pi (L-1)/L\} ∠vmt,∠vmr∈{0,2π/L,⋯,2π(L−1)/L}
其中， L L L 是可用相位的数量。近来，有研究提出了一种新的 耦合相位模子（coupled-phase model） [17]-[19]，该模子指出物理可实现的 STAR-RIS 必须满足以下条件：
∣ v m t ∣ ∣ v m r ∣ cos ( ∠ v m t − ∠ v m r ) = 0 |v_m^t| |v_m^r| \cos(\angle v_m^t - \angle v_m^r) = 0 ∣vmt∣∣vmr∣cos(∠vmt−∠vmr)=0
在 ES模式 下，这等价于：
∠ v m t − ∠ v m r ∈ { π / 2 ( mod 2 π ) , − π / 2 ( mod 2 π ) } \angle v_m^t - \angle v_m^r \in \{\pi/2 \ (\text{mod } 2\pi), -\pi/2 \ (\text{mod } 2\pi)\} ∠vmt−∠vmr∈{π/2 (mod 2π),−π/2 (mod 2π)}
不同范例的 STAR-RIS 是上述约束条件的混合与匹配的结果，其总结见表 I。
需要留意的是，带有耦合相位的 MS 模式现实上只是基本的 MS STAR-RIS，由于 MS 模式要求：
∣ v m t ∣ , ∣ v m r ∣ ∈ { 0 , 1 } |v_m^t|, |v_m^r| \in \{0, 1\} ∣vmt∣,∣vmr∣∈{0,1}
结合无损功率约束： ∣ v m t ∣ 2 + ∣ v m r ∣ 2 = 1 |v_m^t|^2 + |v_m^r|^2 = 1 ∣vmt∣2+∣vmr∣2=1。可以得出，必须满足以下条件之一： ∣ v m t ∣ = 0 |v_m^t| = 0 ∣vmt∣=0 或 ∣ v m r ∣ = 0 |v_m^r| = 0 ∣vmr∣=0。这使得耦合相位约束： ∣ v m t ∣ ∣ v m r ∣ cos ( ∠ v m t − ∠ v m r ) = 0 |v_m^t||v_m^r| \cos(\angle v_m^t - \angle v_m^r) = 0 ∣vmt∣∣vmr∣cos(∠vmt−∠vmr)=0 自动建立。
对于 TS 模式，由于在某一时间段内，每个 STAR-RIS 元件只使用反射相位或透射相位，因此耦合相位约束不可能存在。因此，带有耦合相位（无论是否带有离散相位约束）的 MS 和 TS STAR-RIS 是不可行的，因此表 I 中未列出这些配置。
令 z {\boldsymbol z} z 表示需要优化的其他通讯资源 v t = [ v 1 t , ⋯ , v M t ] T {\boldsymbol v}^{\mathcal{t}}=\left[{v_{1}^{\mathcal{t}},\cdots,v_{M}^{\mathcal{t}}} \right]^{T} vt=[v1t,⋯,vMt]T， v r = [ v 1 r , ⋯ , v M r ] T {\boldsymbol v}^{\mathcal{r}}=\left[{v_{1}^{\mathcal{r}},\cdots,v_{M}^{\mathcal{r}}} \right]^{T} vr=[v1r,⋯,vMr]T 分别作为聚集 { v m t } m = 1 M \{v_m^t\}_{m=1}^M {vmt}m=1M 和 { v m r } m = 1 M \left\{{v_{m}^{\mathcal{r}}}\right\}_{m=1}^{M} {vmr}m=1M 的简写表示，则一个涉及 STAR-RIS 的通用优化题目可以表述为：
min { z , v t , v r , λ t , λ r } F ( z , v t , v r , λ t , λ r ) s . t . λ t ∣ v m t ∣ 2 + λ r ∣ v m r ∣ 2 = 1 , { ∣ v m t ∣ = ∣ v m r ∣ = 1 , { λ t , λ r } ≥ 0 , i f T S ∣ v m t ∣ , ∣ v m r ∣ ∈ { 0 , 1 } , λ t = λ r = 1 , i f M S ∠ v m t − ∠ v m r ∈ { π / π 2 2 ( m o d 2 π ) , − π / π 2 2 ( m o d 2 π ) } , λ t = λ r = 1 , i f E S 000000 ∠ v m t , ∠ v m r ∈ { 0 , 2 π / 2 π L L , ⋯ , 2 π ( L − 1 ) / 2 π ( L − 1 ) L L } , ( z , v t , v r , λ t , λ r ) ∈ Ω . \begin{align*} & \tag{1a} \mathop{\min}\limits_{\left\{{{{\boldsymbol z}},{{\boldsymbol v}}^{{\mathcal{t}}}, {{\boldsymbol v}}^{{\mathcal{r}}},{{\lambda}^{{\mathcal{t}}}}, {{\lambda}^{{\mathcal{r}}}}}\right\}}\;\;{\rm{{\cal F}}}\left({{{\boldsymbol z}}, {{\boldsymbol v}}^{{\mathcal{t}}},{{\boldsymbol v}}^{{\mathcal{r}}},{{\lambda}^{{\mathcal{t} }}},{{\lambda}^{{\mathcal{r}}}}}\right) \\ & \tag{1b}\;\;\;\;\;\;\;\;s.t.\;\;\;{\lambda^{\mathcal{t}}}{\left|{v_{m}^{ \mathcal{t}}}\right|^{2}}+{\lambda^{\mathcal{r}}}{\left|{v_{m}^{\mathcal{r}}} \right|^{2}}=1, \\ &\qquad \left\{ \begin{array}{l} \left|{v_{m}^{\mathcal{t}}}\right|=\left|{v_{m}^{\mathcal{r}}} \right|=1,\left\{{{{\lambda}^{{\mathcal{t}}}},{{\lambda}^{{\mathcal{r} }}}}\right\}\geq 0,\qquad {\rm if\;TS} \\[3pt] \left|{v_{m}^{\mathcal{t}}}\right|,\left|{v_{m}^{\mathcal{r}}} \right|\in\left\{{0,1}\right\},{{\lambda}^{{\mathcal{t}}}}={{ \lambda}^{{\mathcal{r}}}}=1,\quad {\rm if\;MS} \\[3pt] \begin{array}{l}\angle v_{m}^{\mathcal{t}}-\angle v_{m}^{\mathcal{r}}\in \left\{{{\pi\mathord{\left/{\vphantom{\pi 2}}\right.\kern-1.2pt}2} \left({\bmod 2\pi}\right),}\right. \\ \left.-{\pi\mathord{\left/{\vphantom{\pi 2}}\right.\kern-1.2pt}2}\left({\bmod 2 \pi}\right)\right\},{\lambda^{\mathcal{t}}}={\lambda^{\mathcal{r}}}=1, \end{array} \quad {\rm if\;ES} \end{array}\right. \phantom{000000} \tag{1c}\\ &\qquad \angle v_{m}^{\mathcal{t}},\angle v _{m}^{\mathcal{r}}\in\left\{{0,{{2\pi}\mathord{\left/{\vphantom{{2\pi} L}}\right.\kern-1.2pt}L},\cdots,{{2\pi\left({L-1}\right)}\mathord{\left/{ \vphantom{{2\pi\left({L-1}\right)}L}}\right.\kern-1.2pt}L}}\right\}, \tag{1d}\\ &\qquad \left({{\boldsymbol z},{{{\boldsymbol v}}^{\mathcal{t}}},{{{\boldsymbol v}}^{\mathcal{r}}},{{\lambda}^{\mathcal{t}}},{{\lambda}^{\mathcal{r}}} }\right)\in{\boldsymbol{\Omega}}. \tag{1e} \end{align*} {z,vt,vr,λt,λr}minF(z,vt,vr,λt,λr)s.t.λt vmt 2+λr∣vmr∣2=1,⎩ ⎨ ⎧∣vmt∣=∣vmr∣=1,{λt,λr}≥0,ifTS∣vmt∣,∣vmr∣∈{0,1},λt=λr=1,ifMS∠vmt−∠vmr∈{π/π22(mod2π),−π/π22(mod2π)},λt=λr=1,ifES000000∠vmt,∠vmr∈{0,2π/2πLL,⋯,2π(L−1)/2π(L−1)LL},(z,vt,vr,λt,λr)∈Ω.(1a)(1b)(1c)(1d)(1e)
在该优化题目中， F \mathcal F F 是目的函数，并假设其有下界（这是一个显而易见的假设，由于题目是一个最小化题目）。公式 (1) 中的目的函数可以表示系统性能的不同形式，比方：功耗 [15], [19]、均方误差函数 [26] 。
另一方面，以下性能指标也常被用作优化目的，但由于这些指标需要最大化而非最小化，使用它们时需要在公式 (1) 中添加负号：总速率（Sum-rate）[14], [20] 、频谱服从（Spectral efficiency）[29]、保密容量（Secrecy capacity）[30], [34]
其中， Ω \boldsymbol{{\Omega}} Ω 是 z z z、 v t {\boldsymbol v}^{\mathcal{t}} vt、 v r {\boldsymbol v}^{\mathcal{r}} vr、 λ t \lambda^t λt、 λ r \lambda^r λr 的耦合约束聚集。约束条件 (1b) 是涵盖所有三种 STAR-RIS 模式的通用表达形式，具体如下：
B. A Penalty-based Reformulation of (1)
留意到题目 (1) 中，由于约束 (1c) 包罗二元选择，因此求解 STARRIS 系数 v t \boldsymbol v_t vt 和 v r \boldsymbol v_r vr 具有肯定的寻衅性。更具体地说，(1c) 中的第二个约束和第三个约束分别对振幅和两个相位之差举行了二元选择。此外，(1d) 中可能出现的离散相位也使得题目 (1) 成为一个混合整数优化题目（mixed integer optimization problem）。
为了解决这个题目，我们提出采用辅助向量 φ t , φ r ∈ C M × 1 \boldsymbol \varphi_t, \boldsymbol \varphi_r \in \mathbb{C}^{M \times 1} φt,φr∈CM×1，并结合一个惩罚项来处置惩罚约束 (1c) 和 (1d)。重新表述后的题目可以写成如下形式：
min { z , v t , v r , φ t , φ r , λ t , λ r } F ( z , v t , v r , λ t , λ r ) + γ 2 ∑ p = t , r ∣ v p − φ p ∣ 2 2 s . t . { ∣ φ m t ∣ = ∣ φ m r ∣ = 1 , { λ t , λ r } ≥ 0 , i f T S ∣ φ m t ∣ , ∣ φ m r ∣ ∈ { 0 , 1 } , λ t = λ r = 1 , i f M S ∠ φ m t − ∠ φ m r ∈ { π / π 2 2 ( m o d 2 π ) , − π / π 2 2 ( m o d 2 π ) } , λ t = λ r = 1 , i f E S ∠ φ m t , ∠ φ m r ∈ { 0 , 2 π / 2 π L L , ⋯ , 2 π ( L − 1 ) / 2 π ( L − 1 ) L L } , ( 1b ) , ( 1e ) \begin{align*} &\mathop{\min}\limits_{\left\{{\scriptstyle{ \boldsymbol{z}},{{\boldsymbol{v}}^{{\mathcal{t}}}},{{\boldsymbol{v}}^{{ \mathcal{r}}}},{\boldsymbol{\varphi}}^{{\mathcal{t}}},\atop\scriptstyle{ {\boldsymbol{\varphi}}^{{\mathcal{r}}}},{\lambda}^{{\mathcal{t}}},{\lambda^{{ \mathcal{r}}}}} \right\}}{{\rm{{\cal F}}}}\left({{\boldsymbol{z}},{{ \boldsymbol{v}}^{\mathcal{t}}},{{\boldsymbol{v}}^{\mathcal{r}}},{{\lambda}^{ \mathcal{t}}},{{\lambda}^{\mathcal{r}}}}\right)+\frac{\gamma}{2}\sum\limits_{{ \mathcal{p}}={\mathcal{t}},{\mathcal{r}}}{\left|{{{\boldsymbol{v}}^{\mathcal{p }}} {-} {{\boldsymbol{\varphi}}^{\mathcal{p}}}}\right|_{2}^{2}}\tag{2a} \\ &\;\;\;\;\;\;\;\;s.t.\;\;\;\left\{\begin{array}{l}\left|{\varphi_{m}^{\mathcal{t}}}\right|=\left|{\varphi _{m}^{\mathcal{r}}}\right|=1,\left\{{{{\lambda}^{{\mathcal{t}}}},{{\lambda}^{{ \mathcal{r}}}}}\right\}\geq 0,\quad\;\;\;{\rm if\;TS}\\ \left|{\varphi_{m}^{\mathcal{t}}}\right|,\left|{\varphi_{m}^{ \mathcal{r}}}\right|\in\left\{{0,1}\right\},{{\lambda}^{{\mathcal{t}}}}={{ \lambda}^{{\mathcal{r}}}}=1,\;\;\;\;{\rm if\;MS}\\ \begin{array}{l}\angle\varphi_{m}^{\mathcal{t}}-\angle\varphi_{m }^{\mathcal{r}}\in\left\{{{\pi\mathord{\left/{\vphantom{\pi 2}}\right.}2}\left({\bmod 2\pi}\right),}\right.\\ \left.{-{\pi\mathord{\left/{\vphantom{\pi 2}}\right.}2}\left({\bmod 2 \pi}\right)}\right\},{\lambda^{\mathcal{t}}}={\lambda^{\mathcal{r}}}=1,\end{array}\;\;\;\;\;\;{\rm if\;ES}\end{array}\right.\tag{2b} \\ &\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\angle\varphi _{m}^{\mathcal{t}},\angle\varphi_{m}^{\mathcal{r}}\in\left\{{0,{{2\pi}\mathord {\left/{\vphantom{{2\pi}L}}\right.}L},\cdots,{{2\pi\left({L-1} \right)}\mathord{\left/{\vphantom{{2\pi\left({L-1}\right)}L}}\right.}L}}\right\}, (\text{1b}), (\text{1e})\tag{2c} \end{align*} {φr,λt,λrz,vt,vr,φt,}minF(z,vt,vr,λt,λr)+2γp=t,r∑∣vp−φp∣22s.t.⎩ ⎨ ⎧∣φmt∣=∣φmr∣=1,{λt,λr}≥0,ifTS∣φmt∣,∣φmr∣∈{0,1},λt=λr=1,ifMS∠φmt−∠φmr∈{π/π22(mod2π),−π/π22(mod2π)},λt=λr=1,ifES∠φmt,∠φmr∈{0,2π/2πLL,⋯,2π(L−1)/2π(L−1)LL},(1b),(1e)(2a)(2b)(2c)
其中， γ \gamma γ 是惩罚系数。当惩罚系数增大时，RIS 系数向量 v t \boldsymbol v_t vt 和 v r \boldsymbol v_r vr 将被迫取与辅助向量 φ t , φ r \boldsymbol \varphi_t, \boldsymbol \varphi_r φt,φr 雷同的值，从而使 v t \boldsymbol v_t vt 和 v r \boldsymbol v_r vr 满足约束 (1c) 和 (1d)。
留意到 (2) 中 { φ t , φ r } \{\boldsymbol{\varphi}_t, \boldsymbol{\varphi}_r\} {φt,φr} 和 { z , v t , v r , λ t , λ r } \{\boldsymbol{z}, \boldsymbol{v}_t, \boldsymbol{v}_r, {\lambda}_t, {\lambda}_r\} {z,vt,vr,λt,λr} 的约束并未耦合，因此可以采用 BCD 框架来处置惩罚该题目，该框架涉及交替求解以下两个子题目：
P 1 : min { φ t , φ r } ∣ v t − φ t ∣ 2 2 + ∣ v r − φ r ∣ 2 2 s . t . { ∣ φ m t ∣ = ∣ φ m r ∣ = 1 , i f T S ∣ φ m t ∣ , ∣ φ m r ∣ ∈ { 0 , 1 } , i f M S ∠ φ m t − ∠ φ m r ∈ { π / π 2 2 ( m o d 2 π ) , − π / π 2 2 ( m o d 2 π ) } , i f E S ∠ φ m t , ∠ φ m r ∈ { 0 , 2 π / 2 π L L , ⋯ , 2 π ( L − 1 ) / 2 π ( L − 1 ) L L } . \begin{align*} \begin{split}\displaystyle{\rm P1}:\mathop{\min}\limits_{\left\{{{{\boldsymbol\varphi} ^{\mathcal{t}}},{{\boldsymbol\varphi}^{\mathcal{r}}}}\right\}}&\displaystyle{\left|{{{ \boldsymbol v}^{\mathcal{t}}} {-} {{\boldsymbol\varphi}^{\mathcal{t}}}}\right|_{2}^{2}{+} \left|{{{\boldsymbol v}^{\mathcal{r}}} {-} {{\boldsymbol\varphi}^{\mathcal{r}}}}\right|_{2}^ {2}}\\ \displaystyle s.t.&\displaystyle\left\{\begin{array}{l}\left|{\varphi_{m}^{\mathcal{t}}}\right|=\left|{\varphi_{m}^{\mathcal{r}} }\right|=1,\qquad\qquad\qquad\;\;\;\;\;{\rm if\;TS}\\ \left|{\varphi_{m}^{\mathcal{t}}}\right|,\left|{\varphi_{m}^{ \mathcal{r}}}\right|\in\left\{{0,1}\right\},\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; \;\;\;\;\;\;{\rm if\;MS}\\ \begin{array}{l}\angle\varphi_{m}^{\mathcal{t}}-\angle\varphi_{m}^{\mathcal{ r}}\in\left\{{{\pi\mathord{\left/{\vphantom{\pi 2}}\right.\kern-1.2pt}2}\left ({\bmod 2\pi}\right),}\right.\\ \left.{-{\pi\mathord{\left/{\vphantom{\pi 2}}\right.\kern-1.2pt}2}\left({\bmod 2 \pi}\right)}\right\},\end{array}\;\;\;\;\;\;{\rm if\;ES}\end{array}\right.\\ &\displaystyle\;\angle\varphi_{m}^{\mathcal{t}},\angle\varphi_{m}^{\mathcal{r} }\in\left\{{0,{{2\pi}\mathord{\left/{\vphantom{{2\pi}L}}\right.\kern-1.2pt}L}, \cdots,{{2\pi\left({L-1}\right)}\mathord{\left/{\vphantom{{2\pi\left({L-1} \right)}L}}\right.\kern-1.2pt}L}}\right\}.\end{split} \end{align*} P1:{φt,φr}mins.t. vt−φt 22+∣vr−φr∣22⎩ ⎨ ⎧∣φmt∣=∣φmr∣=1,ifTS∣φmt∣,∣φmr∣∈{0,1},ifMS∠φmt−∠φmr∈{π/π22(mod2π),−π/π22(mod2π)},ifES∠φmt,∠φmr∈{0,2π/2πLL,⋯,2π(L−1)/2π(L−1)LL}.
III. OPTIMIZING AUXILIARY VARIABLES φ \boldsymbol φ φ IN P1
如表 I 所示，公式 (3) 的约束条件的不同组合导致了不同范例的 STAR-RIS。以下将讨论分为两种情况。起首，我们思量没有耦合相位约束的情况，这对应于表 I 中的第 1 至第 6 种情况。然后，我们讨论表 I 中的第 7 和第 8 种具有耦合相位的情况。
无耦合相位的 STAR-RIS： 由此产生的子题目为：
min φ m t , φ m r ∣ v m t − φ m t ∣ 2 + ∣ v m r − φ m r ∣ 2 s . t . { ∣ φ m t ∣ = ∣ φ m r ∣ = 1 , i f T S ∣ φ m t ∣ , ∣ φ m r ∣ = { 0 , 1 } o r { 1 , 0 } , i f M S φ m t , φ m r ∈ C , i f E S ∠ φ m t , ∠ φ m r ∈ { 0 , 2 π / 2 π L L , ⋯ , 2 π ( L − 1 ) / 2 π ( L − 1 ) L L } . \begin{align*} \mathop{\min}\limits_{\varphi_{m}^{\mathcal{t}},\varphi_{m}^{ \mathcal{r}}} & \left|{v_{m}^{\mathcal{t}}-\varphi_{m}^{\mathcal{t}}}\right|^{2}+ \left|{v_{m}^{\mathcal{r}}-\varphi_{m}^{\mathcal{r}}}\right|^{2} \tag{4a}\\ s.t. & \left\{\begin{array}{l}\left|{\varphi_{m}^{ \mathcal{t}}}\right|=\left|{\varphi_{m}^{\mathcal{r}}}\right|=1,\;\;\;\;\;\; \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\rm if\;TS}\\[3pt] \left|{\varphi_{m}^{\mathcal{t}}}\right|,\left|{\varphi_{m}^{ \mathcal{r}}}\right|=\left\{{0,1}\right\}{\rm\;or\;}\left\{{1,0}\right\},\;\;\;\;\;{\rm if\;MS}\\[3pt] \varphi_{m}^{\mathcal{t}},\varphi_{m}^{\mathcal{r}}\in{\mathbb{C}},\;\;\;\;\; \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\rm if\;ES} \end{array}\right. \tag{4b}\\ & \; \angle\varphi_{m}^{\mathcal{t}},\angle\varphi_{m}^{\mathcal{r}} \in\left\{{0,{{2\pi}\mathord{\left/{\vphantom{{2\pi}L}}\right.\kern-1.2pt}L},\cdots,{{2\pi\left({L-1}\right)}\mathord{\left/{\vphantom{{2\pi\left({L -1}\right)}L}}\right.\kern-1.2pt}L}}\right\}. \tag{4c} \end{align*} φmt,φmrmins.t. vmt−φmt 2+∣vmr−φmr∣2⎩ ⎨ ⎧∣φmt∣=∣φmr∣=1,ifTS∣φmt∣,∣φmr∣={0,1}or{1,0},ifMSφmt,φmr∈C,ifES∠φmt,∠φmr∈{0,2π/2πLL,⋯,2π(L−1)/2π(L−1)LL}.(4a)(4b)(4c)
一般来说，这是一个包罗整数变量的题目，由于相位是离散的，因此难以求解。然而，由于公式 (4) 中的幅度约束和相位约束是分离的，我们可以分别对其求解，并通过以下引理获得公式 (4) 的闭式解。
引理 1. 定义 α m t = Proj Θ ( ∠ v m t ) \alpha_m^t = \text{Proj}_\Theta (\angle v_m^t) αmt=ProjΘ(∠vmt)， α m r = Proj Θ ( ∠ v m r ) \alpha_m^r = \text{Proj}_\Theta (\angle v_m^r) αmr=ProjΘ(∠vmr)， β m t = ∣ v m t ∣ cos ( α m t − ∠ v m t ) \beta_m^t = |v_m^t| \cos (\alpha_m^t - \angle v_m^t) βmt=∣vmt∣cos(αmt−∠vmt)， β m r = ∣ v m r ∣ cos ( α m r − ∠ v m r ) \beta_m^r = |v_m^r| \cos (\alpha_m^r - \angle v_m^r) βmr=∣vmr∣cos(αmr−∠vmr)。公式 (4) 的最优解为：
φ m t = e j α m t , φ m r = e j α m r , i f T S φ m t = 1 + s g n ( β m t − β m r ) 2 e j α m t , φ m r = 1 + s g n ( β m r − β m t ) 2 e j α m r , i f M S φ m t = β m t e j α m t , φ m r = β m r e j α m r , i f E S \begin{align*} \begin{array}{l}\varphi_{m}^{\mathcal{t}}={e^{j\alpha_{m}^{ \mathcal{t}}}},\varphi_{m}^{\mathcal{r}}={e^{j\alpha_{m}^{\mathcal{r}}}},\;\; \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; {\rm if\;TS}\\[3pt] \varphi_{m}^{\mathcal{t}}=\frac{{1+{\mathop{\rm sgn}}\left({\beta_ {m}^{\mathcal{t}} -\beta_{m}^{\mathcal{r}}}\right)}}{2}{e^{j\alpha_{m}^{ \mathcal{t}}}},\varphi_{m}^{\mathcal{r}}=\frac{{1+{\mathop{\rm sgn}}\left({ \beta_{m}^{\mathcal{r}} -\beta_{m}^{\mathcal{t}}}\right)}}{2}{e^{j\alpha_{m}^{ \mathcal{r}}}}, {\rm if\;MS}\\[3pt] \varphi_{m}^{\mathcal{t}}=\beta_{m}^{\mathcal{t}}{e^{j\alpha_{m}^{\mathcal{t}} }},\varphi_{m}^{\mathcal{r}}=\beta_{m}^{\mathcal{r}}{e^{j\alpha_{m}^{\mathcal{ r}}}},\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; {\rm if \;ES} \end{array}\tag{5} \end{align*} φmt=ejαmt,φmr=ejαmr,ifTSφmt=21+sgn(βmt−βmr)ejαmt,φmr=21+sgn(βmr−βmt)ejαmr,ifMSφmt=βmtejαmt,φmr=βmrejαmr,ifES(5)
引理 2. If L > 2 and is an even number, the optimal solution of (6) is
φ m t = ∣ v m t ∣ cos ( θ m t − ∠ v m t ) e j θ m t , φ m r = ∣ v m r ∣ ∣ sin ( θ m t − ∠ v m r ) ∣ e j ( θ m t − π 2 s g n ( sin ( θ m t − ∠ v m r ) ) ) , \begin{align*} \begin{array}{l}\varphi_{m}^{\mathcal{t}}=\left|{v_{m}^{\mathcal{t}}}\right| \cos\left({\theta_{m}^{\mathcal{t}}-\angle v_{m}^{\mathcal{t}}}\right){e^{j \theta_{m}^{\mathcal{t}}}},\\ \varphi_{m}^{\mathcal{r}}=\left|{v_{m}^{\mathcal{r}}}\right|\left|{\sin\left({ \theta_{m}^{\mathcal{t}}-\angle v_{m}^{\mathcal{r}}}\right)}\right|{e^{j\left ({\theta_{m}^{\mathcal{t}}-\frac{\pi}{2}{\mathop{\rm sgn}}\left({\sin\left({ \theta_{m}^{\mathcal{t}}-\angle v_{m}^{\mathcal{r}}}\right)}\right)}\right)}}, \end{array}\tag{7} \end{align*} φmt=∣vmt∣cos(θmt−∠vmt)ejθmt,φmr=∣vmr∣∣sin(θmt−∠vmr)∣ej(θmt−2πsgn(sin(θmt−∠vmr))),(7)
其中 θ m t = P r o j Θ ( ∠ v m t − b m / − b m 2 2 + π / π 2 2 ) \theta_{m}^{{\mathcal{t}}}={\rm Proj}_{\boldsymbol{\Theta}}\left({\angle v_{m} ^{\mathcal{t}}{-{b_{m}}\mathord{\left/{\vphantom{{-{b_{m}}}2}}\right.\kern-1.2 pt}2}+{\pi\mathord{\left/{\vphantom{\pi 2}}\right.\kern-1.2pt}2}}\right) θmt=ProjΘ(∠vmt−bm/−bm22+π/π22)， b m = − j ln ( [ ∣ v m r ∣ 2 cos ( 2 ∠ v m t − 2 ∠ v m r ) + ∣ v m t ∣ 2 ] + j ∣ v m r ∣ 2 sin ( 2 ∠ v m t − 2 ∠ v m r ) [ ∣ v m r ∣ 2 cos ( 2 ∠ v m t − 2 ∠ v m r ) + ∣ v m t ∣ 2 ] 2 + [ ∣ v m r ∣ 2 sin ( 2 ∠ v m t − 2 ∠ v m r ) ] 2 ) {b_{m}}=-j\ln\left({{\frac{{\left[{{\left|{v_{m}^{\mathcal{r}}}\right|}^{2}}\cos\left({2 \angle v_{m}^{\mathcal{t}}-2\angle v_{m}^{\mathcal{r}}}\right)+\left|{v_{m}^{ \mathcal{t}}}\right|^{2}\right]+j{{\left|{v_{m}^{\mathcal{r}}}\right|}^{2}} \sin\left({2\angle v_{m}^{\mathcal{t}}-2\angle v_{m}^{\mathcal{r}}}\right)}}{{ \sqrt{{{\left[{{{\left|{v_{m}^{\mathcal{r}}}\right|}^{2}}\cos\left({2\angle v_ {m}^{\mathcal{t}}-2\angle v_{m}^{\mathcal{r}}}\right)+{{\left|{v_{m}^{\mathcal {t}}}\right|}^{2}}}\right]}^{2}}+{{\left[{{{\left|{v_{m}^{\mathcal{r}}}\right| }^{2}}\sin\left({2\angle v_{m}^{\mathcal{t}}-2\angle v_{m}^{\mathcal{r}}} \right)}\right]}^{2}}}}}}}\right) bm=−jln([∣vmr∣2cos(2∠vmt−2∠vmr)+∣vmt∣2]2+[∣vmr∣2sin(2∠vmt−2∠vmr)]2 [∣vmr∣2cos(2∠vmt−2∠vmr)+∣vmt∣2]+j∣vmr∣2sin(2∠vmt−2∠vmr))。
Proof. See Appendix C in supplementary material.
这里我必须要吐槽一下作者，这里看似 b m {b_{m}} bm 的结果没有题目，现实上根本是无法仿真的， ln \ln ln 函数中有虚数在MATLAB中是无法运行的。同时我下载了所谓 supplementary material，我们发现Appendix C中对于 χ m \chi_{m} χm 的推导是错误的，精确的结果应该是 χ m = [ ∣ v m t ∣ 2 cos ( 2 ∠ v m t − 2 ∠ v m t ) − ∣ v m t ∣ 2 ] 2 + [ ∣ v m r ∣ 2 sin ( 2 ∠ v m t − 2 ∠ v m r ) ] 2 \chi_m=\left[\left|v_m^t\right|^2\cos\left(2\angle v_m^t-2\angle v_m^t\right)-\left|v_m^t\right|^2\right]^2+\left[\left|v_m^r\right|^2\sin\left(2\angle v_m^t-2\angle v_m^r\right)\right]^2 χm=[∣vmt∣2cos(2∠vmt−2∠vmt)−∣vmt∣2]2+[∣vmr∣2sin(2∠vmt−2∠vmr)]2。此外，此处相角的导出也非常反人类的直觉，他是根据 e j b m = cos b m + j sin b m e^{jb_m}=\cos b_m+j\sin b_m ejbm=cosbm+jsinbm，得到 b m = − j ln ( cos b m + j sin b m ) b_{m}=-j\ln(\cos b_m+j\sin b_m) bm=−jln(cosbm+jsinbm)，只能说bro 肯定觉得自己很幽默。
引理 1 和引理 2 涵盖了离散和一连相位的 STAR-RIS。对于后者情况，相当于取 L → ∞ L \to \infty L→∞，并且引理 1 和引理 2 中的投影函数可以省略。需要留意的是，对于最简单的 ES STAR-RIS（表 I 中的第五种情况），由于公式 (3) 中没有涉及相位和幅度的约束，辅助变量 φ t \varphi^t φt 和 φ r \varphi^r φr 的最优解将始终即是 v t v^t vt 和 v r v^r vr。因此，惩罚循环只需实行一次，这就简化为很多现有文献中的传统无惩罚设计 [28], [44]。
由于 P1 的解覆盖了所有现有范例的 STAR-RIS，一旦解决了 P2，就可以轻松比较所有 STAR-RIS 范例的系统性能。这是第一次以统一的方式启用这样的比较。这一点很重要，由于如果我们分别在不同的 STAR-RIS 模式下优化资源分配题目 (1)，性能的差异可能不仅来自操作模式的不同，还来自算法的不同。因此，所提出的框架为通讯系统研究人员提供了一种统一的方法，以决定哪种范例的 STAR-RIS 最恰当特定场景或应用。此类比较的示例将在第 V 节中提供。
V. SIMULATION RESULTS AND DISCUSSIO
在本部分中，我们通过仿真评估下行链路的总速率传输性能。所有题目实例均在Windows x64桌面平台上使用Matlab-R2023a举行仿真，硬件配置为2.8 GHz CPU和16 GB RAM。仿真结果通过对100次仿真试验的均匀值获得。仿真中使用了并行处置惩罚指数 P = 4 P=4 P=4。
在仿真设置中，基站（BS）和STAR-RIS分别位于坐标 ( 0 , 20 m ) (0, 20\text{m}) (0,20m)和 ( 40 m , 0 ) (40\text{m}, 0) (40m,0)。STAR-RIS沿 y y y轴放置并垂直于地面。反射和透射用户在间隔STAR-RIS两侧8米范围内均匀分布。参数 { G , h l , d l } \{G, h_l, d_l\} {G,hl,dl}被建模为Rician衰落信道，包罗视距（LoS）和非视距（NLoS）分量[50]。以 h l h_l hl为例，Rician衰落信道模子为：
h l = ν h l κ h + 1 ( κ h h l LoS + h l NLoS ) h_l = \sqrt{\frac{\nu h_l}{\kappa_h + 1}} \left( \sqrt{\kappa_h} h_l^{\text{LoS}} + h_l^{\text{NLoS}} \right) hl=κh+1νhl (κh hlLoS+hlNLoS)
其中每个参数的具体定义如下：
- 路径消耗参数
ν h l = L 0 ( d h l d 0 ) − α h \nu h_l = L_0 \left(\frac{d_{hl}}{d_0}\right)^{-\alpha_h} νhl=L0(d0dhl)−αh表示从RIS到第 l l l个用户的间隔相关路径消耗，其中 L 0 = − 30 dB L_0 = -30\text{dB} L0=−30dB表示参考间隔 d 0 = 1 m d_0=1\text{m} d0=1m处的路径消耗， d h l d_{hl} dhl是STAR-RIS与第 l l l个用户之间的间隔， α h = 2.2 \alpha_h = 2.2 αh=2.2为RIS-用户链路的路径消耗指数。相应地，BS-RIS链路和BS-用户链路的路径消耗指数分别为2.2和3.6。
- Rician因子
κ h \kappa_h κh为RIS-用户链路的Rician因子，较高的Rician因子表示较强的LoS分量。当Rician因子为0时，表明不存在LoS信号，信道简化为Rayleigh衰落。在仿真中，RIS-用户链路的Rician因子 κ h \kappa_h κh设定为5，相应地，BS-RIS链路和BS-用户链路的Rician因子分别为5和0。
- LoS分量
h l LoS h_l^{\text{LoS}} hlLoS建模为阵列响应的导向矢量，其第 m m m个元素为：
h l LoS [ m ] = e j 2 π ( m − 1 ) d A sin ( ω ) / λ h_l^{\text{LoS}}[m] = e^{j2\pi(m-1)d_A \sin(\omega)/\lambda} hlLoS[m]=ej2π(m−1)dAsin(ω)/λ
其中， ω \omega ω表示到达角（AoA）或脱离角（AoD）。在仿真中，设 d A / λ = 1 / 2 d_A/\lambda=1/2 dA/λ=1/2，且 ω \omega ω在 [ 0 , 2 π ) [0, 2\pi) [0,2π)范围内均匀分布。而 h l NLoS h_l^{\text{NLoS}} hlNLoS表示NLoS分量信号，其每个元素服从归一化复高斯分布。
为了制止重复形貌图例中的参数， ( M , N , K r , K t , L , P B S , σ l 2 ) (M, N, K_r, K_t, L, P_{BS}, \sigma^2_l) (M,N,Kr,Kt,L,PBS,σl2)的设置详见各图的图例分析。
在实现过程中，采用指数增长的 γ \gamma γ更新方式： γ ← γ × 1.2 \gamma \gets \gamma \times 1.2 γ←γ×1.2，这是现有涉及惩罚方法的研究中广泛采用的[51]-[54]。惩罚参数不停增长，直至交替优化中的 v t , v r v_t, v_r vt,vr与 ϕ t , ϕ r \phi_t, \phi_r ϕt,ϕr的差距低于预定义的容差范围。已知只要惩罚参数在每次迭代中的增长幅度保持在适度范围内，性能将保持不变。 仿真结果如图4(a)和4(b)所示，分别对应于ES STAR-RIS和MS STAR-RIS。
图5比较了八种STAR-RIS的总速率性能，采用了提出的通用优化框架。起首可以看到，额外的离散相位约束仅略微降低了网络吞吐量（ES的两离散相位约束下降2.94%，MS和TS分别下降6.16%和8.31%），这与文献[55]中的结论不同。文献[55]认为稀疏相位（少于8相位）会显著影响性能，原因在于[55]的结论基于从一连相位解到近来离散相位的量化。在允许的离散相位数少少（比方仿真中仅允许2个相位）的情况下，一连相位与近来离散值差异很大。由于量化是独立应用于每个RIS单元的，与一连优化解相比，累积性能退化将非常显著。相比之下，本文提出的算法中辅助变量 ϕ t , ϕ r \phi_t, \phi_r ϕt,ϕr仅受到相位约束的影响，这可以被视为对 v t , v r v_t, v_r vt,vr的相位校正，而不仅仅是相位约束的直接作用。通过本文提出的统一惩罚框架，初次揭示了这一洞见，并提供了便于比较的方式。
与直接将 v t , v r v_t, v_r vt,vr量化为 ϕ t , ϕ r \phi_t, \phi_r ϕt,ϕr不同，本文引入了惩罚项来实现间接量化，这为 v v v和 ϕ \phi ϕ在搜索更优解时提供了更大的自由度。这种方式制止了显著的性能损失，表明量化是导致性能退化的更重要因素，而不是离散相位本身。
其次，需要留意的是，MS和TS模式中的离散相位约束导致的性能损失比ES模式更显着。这可能是由于MS和TS STAR-RIS系数上额外施加了约束条件(1c)。在没有幅度约束的情况下，ES STAR-RIS可以通过优化幅度来弥补离散相位约束带来的损失。
第三（Thirdly），从图5的放大部分可以观察到，与非耦合相位情况相比，STAR-RIS中一连耦合相位约束（continuous coupled-phase constraint）对系统吞吐量几乎没有影响（仅0.89%的损失），这一结果与文献[20]的结论一致。相反，在MS模式下，由0-1幅度约束引入的吞吐量下降更加显着。这一洞见初次被揭示，这是由于提出的统一惩罚框架使得比较更加容易。
图6(a)和6(b)聚焦于具有耦合相位的ES STAR-RIS（即表1中的案例7和8），比较了本文提出的算法与以下算法的收敛性能：逐元素交替优化算法（elementwise-AO）[19]、耦合相位STAR-RIS框架（CP框架）[20]以及基于惩罚的保密波束成形算法（PSB算法）[34]。固然原始PSB算法是为保密总速率设计的，但我们可以对其举行修改以处置惩罚总速率最大化题目。
对于图6(a)中的一连耦合相位情况，本文提出的算法在最短的实行时间内收敛到最高的总速率。相比之下，PSB算法和CP框架的性能略逊一筹，而逐元素交替优化（elementwise-AO）的表现最差。如命题1（第二节B部分）之前所述，CP框架和PSB算法所需的惩罚权重比本文提出的算法更大，从而减缓了它们的收敛速率。另一方面，逐元素交替优化需要在STAR-RIS逐元素级别上求解非凸子题目，这些题目通过SCA和CVX解决，因此计算量非常大。
为了进一步分析本文提出的算法和对比算法的复杂性，表3总结了各种算法的复杂度，其中 I A O I_{AO} IAO表示交替优化（AO）迭代次数。留意，本文框架中的 ς \varsigma ς为二分搜索的精度。从表3可以看出，忽略所有方法共有的项后，本文框架的复杂度在每次迭代中与 M M M呈线性关系，即 O ( M ) O(M) O(M)。相比之下，CP框架和PSB算法在每次迭代中的复杂度至少为 O ( M 2 ) O(M^2) O(M2)。
由于STAR-RIS的单元数量 M M M通常大于基站天线数量 N N N以及用户数量 K r K^r Kr和 K t K^t Kt，CP框架和PSB算法的复杂度会远高于本文提出的框架。固然逐元素交替优化（elementwise AO）算法没有外部惩罚循环，其复杂度包罗 I A O ⋅ O ( M 3.5 ) I_{AO} \cdot O(M^{3.5}) IAO⋅O(M3.5)项，这显著降低了其整体实行速率。
由于本文方法在每次迭代中的复杂度较低，如图6(a)所示，与CP框架和PSB算法相比，本文算法将计算时间减少了近10倍；与逐元素交替优化算法相比，计算时间减少了近100倍。需要留意的是，由于CP框架、PSB算法和本文提出的算法引入了惩罚项，总速率可能不会随着迭代次数单调增长，由于单调性仅适用于包罗惩罚项的目的函数。
免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！更多信息从访问主页：qidao123.com:ToB企服之家，中国第一个企服评测及商务社交产业平台。 |
