RIS辅助MIMO广播信道容量

RIS辅助MIMO广播信道容量

  • 摘要
  • RIS辅助的BC容量
  • 矩阵形式的泰勒展开学习
  • 舒尔补

RIS-Aided Multiple-Input Multiple-Output
Broadcast Channel Capacity论文阅读记录

基于泰勒展开求解了上行容量和最差用户的可达速率,学习其中的展开方法。


摘要

Scalable algorithms are conceived for obtaining the
sum-rate capacity of the reconfigurable intelligent surface (RIS)-
aided multiuser (MU) multiple-input multiple-output (MIMO)
broadcast channel (BC), where a multi-antenna base station (BS)
transmits signals to multi-antenna users with the help of an
RIS equipped with a massive number of finite-resolution pro-
grammable reflecting elements (PREs). As a byproduct, scalable
path-following algorithms emerge for determining the sum-rate
capacity of the conventional MIMO BCs, closing a long-standing
open problem of information theory. The paper also develops
scalable algorithms for maximizing the minimum rate (max-min
rate optimization) of the users achieved by the joint design of
RIS’s PRE and transmit beamforming for such an RIS-aided BC.
The simulations provided confirm the high performance achieved
by the algorithms developed, despite their low computational
complexity.

可扩展算法旨在获得可重构智能表面 (RIS) 辅助多用户 (MU) 多输入多输出 (MIMO) 广播信道 (BC) 的总速率容量,其中多天线基站 (BS) 进行传输借助配备大量有限分辨率可编程反射元件 (PRE) 的 RIS,向多天线用户发送信号。作为副产品,可扩展的路径跟踪算法出现了,用于确定传统 MIMO BC 的总速率容量,解决了信息论中长期存在的开放问题。该论文还开发了可扩展算法,通过 RIS PRE 的联合设计和针对 RIS 辅助 BC 的传输波束成形来最大化用户的最小速率(最大-最小速率优化)。所提供的模拟证实了所开发的算法所实现的高性能,尽管其计算复杂度较低。

RIS辅助的BC容量

K User
RIS N elements
Nr antenna UEs
Nt antenna BS
BS-RIS channel : LoS
RIS-UEs channel : NLoS 因为 RIS 通常位于高层建筑的显眼位置

从 BS 和 RIS 到 UE k 的准静态平坦衰落信道:
在这里插入图片描述
C N r × N t ∋ H k ( θ ) ≜ H ~ R , k R R , k 1 / 2 diag ⁡ ( e ȷ θ ) H ~ B , R + H ~ B , k = H ~ B R , k diag ⁡ ( e ȷ θ ) H B , R + H ~ B , k \begin{aligned} \mathbb{C}^{N_{r} \times N_{t}} \ni \mathcal{H}_{k}(\boldsymbol{\theta}) & \triangleq \tilde{H}_{R, k} R_{R, k}^{1 / 2} \operatorname{diag}\left(e^{\jmath \boldsymbol{\theta}}\right) \tilde{H}_{B, R}+\tilde{H}_{B, k} \\ & =\tilde{H}_{B R, k} \operatorname{diag}\left(e^{\jmath \boldsymbol{\theta}}\right) H_{B, R}+\tilde{H}_{B, k} \end{aligned} CNr×NtHk(θ)H~R,kRR,k1/2diag(eθ)H~B,R+H~B,k=H~BR,kdiag(eθ)HB,R+H~B,k
在这里插入图片描述

θ n ∈ B ≜ { ν 2 π 2 b , ν = 0 , 1 , … , 2 b − 1 } \boldsymbol{\theta}_{n} \in \mathcal{B} \triangleq\left\{\nu \frac{2 \pi}{2^{b}}, \nu=0,1, \ldots, 2^{b}-1\right\} θnB{ν2b2π,ν=0,1,,2b1}

for  n ∈ N ≜ { 1 , … , N } , i.e.  \text { for } n \in \mathcal{N} \triangleq\{1, \ldots, N\} \text {, i.e. }  for nN{1,,N}, i.e. 

θ ∈ B N \boldsymbol{\theta} \in \mathcal{B}^{N} θBN

injective,单射的
脏纸编码:

关于脏纸编码,理论基础是Costa 在1987 的一篇论文 Writing on Dirty paper。 Paper 中介绍了一种思路。在日常中,一张纸上假如已经写了很多内容或者已经有很多污点了,即脏纸,但是再在上面写上一些信息的话,这份信息在未知墨水颜色与位置下还是可读的。因此Costa 想把这种思路拓展且应用到通信领域中。为此他引用了一种新的信道模型,即:

需要传输的信息或数据为w,然后通过一定的编码机制即编码器,此编码器是已知一部分信道干扰s的(s符合一般正态分布),也就是说,编码器与这部分干扰通过一定机制对信息w编码,在编码机制后便会得到需要被传输的信息x,信息x将通过AWGN信道,其干扰为n且符合二项分布。最终在接收端得到信息y。通过一定的解码机制便能将w的判断值算出来。在一堆的数学分析与运算后,Costa便发现,在一定编码译码机制下,受到抑制干扰s与未知干扰n的信道容量可以达到只受到n一样的信道容量。当然,信道是受香农定理约束的。这就提供了一种可能,假如,把信息制作成s,且开发出相应的机制,就能将信息很完美且不损耗信道容量地隐藏,或者复用,亦或者移动通信的 MIMO技术。但,当下对于脏纸编码还是指处于相对的理论研究阶段。对于具体的实施方案也有了相应的论文,但是都不是特别具体的机制(也有可能只是笔者没找到),因此这项科技虽然被讨论很多,却没有还没有实质进展。理论的尝试倒是已经很多了。笔者因为毕业论文也尝试了一些方法尝试。利用CDMA方法进行拓频,利用拓频编码的特性进行累加且mod加法,但是,得到的结果不是特别理想,因此。项目就暂置了。因此当下只写个大概的内容,之后有进展便会继续拓展本博文。也希望假如有前辈从事类似方向的能够一起讨论。
————————————————
原文链接:https://blog.csdn.net/rogerwutao/article/details/25669613





矩阵形式的泰勒展开学习

log ⁡ det ⁡ ( I n d + H d K x H d H ) − log ⁡ det ⁡ ( K v ′ + H e K x H e H ) ≥ log ⁡ det ⁡ ( I n d + H d K x H d H ) − log ⁡ det ⁡ ( K v ′ ( 0 ) + H e K x ( 0 ) H e H ) − Tr ⁡ ( H e H ( K v ′ ( 0 ) + H e K x ( 0 ) H e H ) − 1 H e K x ) + Tr ⁡ ( H e H ( K v ′ ( 0 ) + H e K x ( 0 ) H e H ) − 1 H e K x ( 0 ) ) . \begin{array}{l} \log \operatorname{det}\left(\mathbf{I}_{n_{d}}+\mathbf{H}_{d} \mathbf{K}_{x} \mathbf{H}_{d}^{H}\right)-\log \operatorname{det}\left(\mathbf{K}_{v^{\prime}}+\mathbf{H}_{e} \mathbf{K}_{x} \mathbf{H}_{e}^{H}\right)\\ \ge \log \operatorname{det}\left(\mathbf{I}_{n_{d}}+\mathbf{H}_{d} \mathbf{K}_{x} \mathbf{H}_{d}^{H}\right)-\log \operatorname{det}\left(\mathbf{K}_{v^{\prime}}^{(0)}+\mathbf{H}_{e} \mathbf{K}_{x}^{(0)} \mathbf{H}_{e}^{H}\right) \\ -\operatorname{Tr}\left(\mathbf{H}_{e}^{H}\left(\mathbf{K}_{v^{\prime}}^{(0)}+\mathbf{H}_{e} \mathbf{K}_{x}^{(0)} \mathbf{H}_{e}^{H}\right)^{-1} \mathbf{H}_{e} \mathbf{K}_{x}\right) \\ +\operatorname{Tr}\left(\mathbf{H}_{e}^{H}\left(\mathbf{K}_{v^{\prime}}^{(0)}+\mathbf{H}_{e} \mathbf{K}_{x}^{(0)} \mathbf{H}_{e}^{H}\right)^{-1} \mathbf{H}_{e} \mathbf{K}_{x}^{(0)}\right) . \end{array} logdet(Ind+HdKxHdH)logdet(Kv+HeKxHeH)logdet(Ind+HdKxHdH)logdet(Kv(0)+HeKx(0)HeH)Tr(HeH(Kv(0)+HeKx(0)HeH)1HeKx)+Tr(HeH(Kv(0)+HeKx(0)HeH)1HeKx(0)).

log ⁡ det ⁡ ( K v ′ + H e K x H e H ) ≤ log ⁡ det ⁡ ( K v ′ ( 0 ) + H e K x ( 0 ) H e H ) + Tr ⁡ ( H e H ( K v ′ ( 0 ) + H e K x ( 0 ) H e H ) − 1 H e K x ) − Tr ⁡ ( H e H ( K v ′ ( 0 ) + H e K x ( 0 ) H e H ) − 1 H e K x ( 0 ) ) \log \operatorname{det}\left(\mathbf{K}_{v^{\prime}}+\mathbf{H}_{e} \mathbf{K}_{x} \mathbf{H}_{e}^{H}\right) \le \log \operatorname{det}\left(\mathbf{K}_{v^{\prime}}^{(0)}+\mathbf{H}_{e} \mathbf{K}_{x}^{(0)} \mathbf{H}_{e}^{H}\right) \\ +\operatorname{Tr}\left(\mathbf{H}_{e}^{H}\left(\mathbf{K}_{v^{\prime}}^{(0)}+\mathbf{H}_{e} \mathbf{K}_{x}^{(0)} \mathbf{H}_{e}^{H}\right)^{-1} \mathbf{H}_{e} \mathbf{K}_{x}\right) \\ -\operatorname{Tr}\left(\mathbf{H}_{e}^{H}\left(\mathbf{K}_{v^{\prime}}^{(0)}+\mathbf{H}_{e} \mathbf{K}_{x}^{(0)} \mathbf{H}_{e}^{H}\right)^{-1} \mathbf{H}_{e} \mathbf{K}_{x}^{(0)}\right) logdet(Kv+HeKxHeH)logdet(Kv(0)+HeKx(0)HeH)+Tr(HeH(Kv(0)+HeKx(0)HeH)1HeKx)Tr(HeH(Kv(0)+HeKx(0)HeH)1HeKx(0))

function:
f = log ⁡ ( d e t ( H ⋅ X ⋅ Z + Y ) ) f = \log(\mathrm{det}(H\cdot X\cdot Z+Y)) f=log(det(HXZ+Y))

gradient:
∂ f ∂ X = ( Z ⋅ i n v ( Y + H ⋅ X ⋅ Z ) ⋅ H ) ⊤ \frac{\partial f}{\partial X} = (Z\cdot \mathrm{inv}(Y+H\cdot X\cdot Z)\cdot H)^\top Xf=(Zinv(Y+HXZ)H)

张贤达《矩阵分析与应用》P160
在这里插入图片描述


多元函数的泰勒展开式
摘抄自https://zhuanlan.zhihu.com/p/33316479

实际优化问题的目标函数往往比较复杂。为了使问题简化,通常将目标函数在某点附近展开为泰勒(Taylor)多项式来逼近原函数。

一元函数在点 x k x_k xk处的泰勒展开式为

f ( x ) = f ( x k ) + ( x − x k ) f ′ ( x k ) + 1 2 ! ( x − x k ) 2 f ′ ′ ( x k ) + o n f(x)=f\left(x_{k}\right)+\left(x-x_{k}\right) f^{\prime}\left(x_{k}\right)+\frac{1}{2 !}\left(x-x_{k}\right)^{2} f^{\prime \prime}\left(x_{k}\right)+o^{n} f(x)=f(xk)+(xxk)f(xk)+2!1(xxk)2f′′(xk)+on

二元函数在点 ( x k , y k ) (x_k,y_k) (xk,yk)处的泰勒展开式为:
f ( x , y ) = f ( x k , y k ) + ( x − x k ) f x ′ ( x k , y k ) + ( y − y k ) f y ′ ( x k , y k ) + 1 2 ! ( x − x k ) 2 f x x ′ ′ ( x k , y k ) + 1 2 ! ( x − x k ) ( y − y k ) f x y ′ ′ ( x k , y k ) + 1 2 ! ( x − x k ) ( y − y k ) f y x ′ ′ ( x k , y k ) + 1 2 ! ( y − y k ) 2 f y y ′ ′ ( x k , y k ) + o n \begin{array}{c} f(x, y)=f\left(x_{k}, y_{k}\right)+\left(x-x_{k}\right) f_{x}^{\prime}\left(x_{k}, y_{k}\right)+\left(y-y_{k}\right) f_{y}^{\prime}\left(x_{k}, y_{k}\right) \\ +\frac{1}{2 !}\left(x-x_{k}\right)^{2} f_{x x}^{\prime \prime}\left(x_{k}, y_{k}\right)+\frac{1}{2 !}\left(x-x_{k}\right)\left(y-y_{k}\right) f_{x y}^{\prime \prime}\left(x_{k}, y_{k}\right) \\ +\frac{1}{2 !}\left(x-x_{k}\right)\left(y-y_{k}\right) f_{y x}^{\prime \prime}\left(x_{k}, y_{k}\right)+\frac{1}{2 !}\left(y-y_{k}\right)^{2} f_{y y}^{\prime \prime}\left(x_{k}, y_{k}\right) \\ +o^{n} \end{array} f(x,y)=f(xk,yk)+(xxk)fx(xk,yk)+(yyk)fy(xk,yk)+2!1(xxk)2fxx′′(xk,yk)+2!1(xxk)(yyk)fxy′′(xk,yk)+2!1(xxk)(yyk)fyx′′(xk,yk)+2!1(yyk)2fyy′′(xk,yk)+on

在这里插入图片描述

MM算法中的一阶泰勒展开代函数
在这里插入图片描述

舒尔补

Schur Complement

在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/125661.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

什么是神经网络,它的原理是啥?(1)

参考:https://www.youtube.com/watch?vmlk0rddP3L4&listPLuhqtP7jdD8CftMk831qdE8BlIteSaNzD 视频1: 简单介绍神经网络的基本概念,以及一个训练好的神经网络是怎么使用的 分类算法中,神经网络在训练过程中会学习输入的 pat…

Pmdarima实现单变量时序预测与交叉验证

目录 1. pmdarima实现单变量时间序列预测 2. 时间序列交叉验证 2.1 滚动交叉验证(RollingForecastCV) 2.2 滑窗交叉验证(SildingWindowForecastCV) 1. pmdarima实现单变量时间序列预测 Pmdarima是以statsmodel和autoarima为基础、封装研发出的Python时序分析库、也是现在市…

故障诊断模型 | Maltab实现GRU门控循环单元故障诊断

文章目录 效果一览文章概述模型描述源码设计参考资料效果一览 文章概述 故障诊断模型 | Maltab实现GRU门控循环单元故障诊断 模型描述 利用各种检查和测试方法,发现系统和设备是否存在故障的过程是故障检测;而进一步确定故障所在大致部位的过程是故障定位。故障检测和故障定位…

3ds Max2022安装教程(最新最详细)

目录 一.简介 二.安装步骤 网盘资源见文末 一.简介 3DS Max是由Autodesk公司开发的一款专业三维建模、动画和渲染软件,广泛应用于影视、游戏、建筑和工业设计等领域。 3DS Max的主要特点和功能包括: 三维建模:3DS Max提供了各种强大的建…

如何用思维导图开会

在办公室和会议室使用思维导图会有无数好处。今天我们就聊聊思维导图在开会中的作用? 为什么要在会议中使用思维导图? 思维导图可以帮助我们整理思路。会议通常涉及到复杂的议题和讨论,使用思维导图可以帮助整合和梳理参与者的思路和观点。通…

Vue项目搭建及使用vue-cli创建项目、创建登录页面、与后台进行交互,以及安装和使用axios、qs和vue-axios

目录 1. 搭建项目 1.1 使用vue-cli创建项目 1.2 通过npm安装element-ui 1.3 导入组件 2 创建登录页面 2.1 创建登录组件 2.2 引入css(css.txt) 2.3 配置路由 2.5 运行效果 3. 后台交互 3.1 引入axios 3.2 axios/qs/vue-axios安装与使用 3.2…

priority_queue 的模拟实现

priority_queue 的底层结构 我们已经学习过栈和队列了,他们都是用一种容器适配出来的。今天我们要学习的 prority_queue 也是一个容器适配器。在 priority_queue 的使用部分我们已经知道想要适配出 priority_queue,这个底层的容器必须有以下接口&#x…

安装Python环境

Python 安装包下载地址:https://www.python.org/downloads/ 打开该链接,可以看到有两个版本的 Python,分别是 Python 3.x 和 Python 2.x,如下图所示: Python下载页面截图 图 1 Python 下载页面截图(包含…

基于单片机设计的防煤气泄漏装置

一、前言 煤气泄漏是一个严重的安全隐患,可能导致火灾、爆炸以及对人体健康的威胁。为了提高家庭和工业环境中煤气泄漏的检测和预防能力,设计了一种基于单片机的防煤气泄漏装置。 单片机选择STC89C52作为主控芯片。为了检测煤气泄漏,采用了…

cocosCreator 调用wxAPI 及后台授权设置、获取用户昵称和头像

版本: 3.8.0 语言: TypeScript 环境: Mac 官方文档: 微信官方文档 - 开放能力 微信 API 小游戏环境 在cocosCreator的3.x版本项目开发中,TypeScript最终会被转换为JavaScript语言。 JavaScript的运行时调用的API…

联发科MT6893(天玑1200)_MTK5G芯片规格参数性能_安卓手机主板方案

联发科天玑1200集成MediaTek 5G调制解调器,通过包含6大维度、72个场景测试的德国莱茵TV Rheinland认证,支持高性能5G连接,带给用户全场景的高品质5G连网体验。 进入5G时代,AI多媒体成为主流应用,天玑1200以强劲的平台…

【高效开发工具系列】你真的会使用Mac吗?

💝💝💝欢迎来到我的博客,很高兴能够在这里和您见面!希望您在这里可以感受到一份轻松愉快的氛围,不仅可以获得有趣的内容和知识,也可以畅所欲言、分享您的想法和见解。 推荐:kwan 的首页,持续学…

uniapp实现瀑布流

首先我们要先了解什么是瀑布流: 瀑布流(Waterfall Flow)是一种常见的网页布局方式,也被称为瀑布式布局或砌砖式布局。它通常用于展示图片、博客文章、商品等多个不同大小和高度的元素。 瀑布流布局的特点是每个元素按照从上到下…

现代的简洁,诠释轻奢的精致!福州中宅装饰,福州装修

轻奢风是一种生活新时尚 优雅、低调、舒适、简单,不断地推陈出新 站在时尚的前沿,引领潮流 中宅装饰集团轻奢风格产品 追求高品质生活细节 以设计精致的空间构造营造出 一种优雅、时尚生活氛围 将低调奢华之美注入现代家居设计中 客厅|The Sitt…

世界电信日 | 人大金仓助力中国移动租赁核算系统升级上线

世界电信日 5月17日恰逢第五十四个世界电信日,运营商作为新型基础设施建设以及维护网信安全的主力军,掌握关键核心技术,实现科技自立自强刻不容缓。 作为数据库领域国家队,人大金仓坚持原始创新,低难度、低成本、低风…

早安心语微语早读,保持一颗平常心,坐看云起落花开谢得之淡然

1、保持一颗平常心,坐看云起落花开谢得之淡然,失之坦然,让生命中每一天都充满着阳光和希望! 2、每个人都一样,都有一段独行的日子,或长或短,这都是无可回避的。不必总觉得生命空空荡荡&#xf…

ArcGIS计算土地现状容积率

本文讲解在ArcGIS中,基于建筑数据和地籍边界数据,计算土地容积率。 一、容积率介绍 容积率(Plot Ratio/Floor Area Ratio/Volume Fraction)是指一个小区的地上建筑总面积与净用地面积的比率。又称建筑面积毛密度。 二、数据分析 (1)建筑数据(dwg) (2)地籍边界数据…

利用AI Chat 将电子书自动截屏并保存成pdf文件

电子书如果要下载下来,无非就两种类型的方法,一种是从内部破解,通常是某些极客将软件破解成免费版,但是风险也大。另一种是从外部破解,就是截屏保存,然后将所有图片拼成pdf文件。 如果要将整本电子书截屏保…

LeetCode刷题---简单组(六)

文章目录 🍒题目一 69. x 的平方根🍒解法一🍒解法二🍒题目二 70. 爬楼梯🍒解法一 🍒题目一 69. x 的平方根 🍒解法一 class Solution(object):def mySqrt(self, x):""":type x:…

制作一个简单的C语言词法分析程序

1.分析组成 C语言的程序中,有很单词多符号和保留字。一些单词符号还有对应的左线性文法。所以我们需要先做出一个单词字符表,给出对应的识别码,然后跟据对应的表格来写出程序 2.程序设计 程序主要有循环判断构成。不需推理即可产生的符号我…