【OFDM】Joint Communications and Sensing Design for Multi-Carrier MIMO Systems

1-模型

1.1-预编码矩阵

预编码矩阵的维度

• $k$ 子载波索引
• $\mathbf{F}$ 预编码矩阵，而非向量

• 类比波束赋形，给第$k$个用户发送的beamformer向量是 $\mathbf{w}[k]\in {\mathbb{C}}^{N_{\text{t}}\times 1}$，每个天线是个复数值
• 也可以把$K$个用户的矢量拼成矩阵 $[\mathbf{w}[1],\cdot \cdot \cdot ,\mathbf{w}[K]]\in {\mathbb{C}}^{N_{\text{t}}\times K}$

• $N_{\text{s}}$ 数据流（data streams）的数量. 单纯从beamformer发射端的角度来说，在OFDM中有多个子载波，所以增加子载波维度

• 第 $k$ 子载波的预编码矩阵
• 第 $k$ 子载波有$N_{\text{s}}$ 数据流，需将其映射到$N_{\text{t}}$根天线上
• 似乎是先DFT，再预编码

$$\mathbf{F}[k]\in {\mathbb{C}}^{N_{\text{t}}\times N_{\text{s}}}$$

• 对每个子载波进行功率约束
  $${\Vert \mathbf{F}[k]\Vert }_{\mathcal{F}}^2=P_{\mathrm{BS}}$$

1.1-预编码矩阵

接收向量

• 接收向量的维度应和发送向量$\mathbf{s}[k]$一致，即$N_{\text{t}}\times 1$

$$\mathbf{y}[k]=\mathbf{W}[k{]}^{\mathsf{H}}\mathbf{H}[k]\mathbf{F}[k]\mathbf{s}[k]+\mathbf{W}[k{]}^{\mathsf{H}}\mathbf{n}[k]$$

1.3-可达码率

$${\mathcal{R}}_k={\log }_2\left|{\mathbf{I}}_{N_{\mathrm{r}}}+\frac{P_{\mathrm{BS}}}{{\sigma }_{\mathrm{n}}^2{N}_{\mathrm{s}}}\mathbf{W}[k{]}^{†}\mathbf{H}[k]\mathbf{F}[k]\mathbf{F}[k{]}^{\mathrm{H}}\mathbf{H}[k{]}^{\mathrm{H}}\mathbf{W}[k]\right|$$

2-模型

探测信号的协方差矩阵 $\mathbf{R}[k]$

minimize { R [ k ] } ∑ k ∈ J ∑ t = 1 T ∣ P d ( θ t , f k ) − a ( θ t , f k ) H R [ k ] a ( θ t , f k ) ∣  subject to  [ R [ k ] ] n , n = P B S / N t , n ∈ { 1 , … , N t } , R [ k ] ⪰ 0 , R [ k ] = R [ k ] H , \begin{aligned} & \underset{\substack{ \{ \mathbf{R}[k] \} }}{\textrm{minimize}} \sum_{k \in \mathcal{J}} \sum_{t=1}^T\left|\mathcal{P}_{\mathrm{d}}\left(\theta_t, f_k\right)-\mathbf{a}\left(\theta_t, f_k\right)^{\mathrm{H}} \mathbf{R}[k] \mathbf{a}\left(\theta_t, f_k\right)\right| \\ & \text { subject to }[\mathbf{R}[k]]_{n, n}=P_{\mathrm{BS}} / N_{\mathrm{t}}, \quad n \in\left\{1, \ldots, N_{\mathrm{t}}\right\}, \\ & \mathbf{R}[k] \succeq \mathbf{0}, \mathbf{R}[k]=\mathbf{R}[k]^{\mathrm{H}},\end{aligned} ​{R[k]}​minimize​k∈J∑​t=1∑T​ ​Pd​(θt​,fk​)−a(θt​,fk​)HR[k]a(θt​,fk​) ​ subject to [R[k]]n,n​=PBS​/Nt​,n∈{1,…,Nt​},R[k]⪰0,R[k]=R[k]H,​

复现

直接运行会报错，需要按照下面要求安装manopt

函数或变量 'obliquecomplexfactory' 无法识别

Matlab添加 Manopt toolbox