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.mfbz.cn/a/111322.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

UE4 使用材质后期 制作玻璃有雨效果

效果展示,其实这是一个动画效果 以上为所有逻辑 拿到TexCoord给到Panner,Time和Speed都是通过下面计算而来,后面讲,再拿到时间和速度值过后,加上扰动值,最后取G值,因为雨事从上而下的动&#xf…

【AIFEM案例操作】电器盒谐响应分析

AIFEM是由天洑自主研发的一款通用的智能结构仿真软件,助力用户解决固体结构相关的静力学、动力学、振动、热力学等实际工程问题,软件提供高效的前后处理工具和高精度的有限元求解器,帮助用户快速、深入地评估结构的力学性能,加速产…

微信小程序项目案例之导游证考试刷题小程序

前言 很多计算机专业的同学在做毕设选题时不知道该如何选题,有的同学是已经选择了要开发一款小程序,但是又不知道开发哪类小程序。本篇将为大家介绍一个小程序的开发方向,考试刷题类小程序是目前比较火的小程序项目之一,在小程序…

JavaScript基础知识19——循环结构:while循环

哈喽,你好,我是雷工。 本节学习JavaScript基础语法的循环结构:while循环,以下为学习笔记。 while循环 循环概念:重复执行一些操作; 循环特征:不断地重复; while:在…期间…

面向Three.js开发者的3D自动纹理化开发包

DreamTexture.js 是面向 three.js 开发者的 3D 模型纹理自动生成与设置开发包,可以为 webGL 应用增加 3D 模型的快速自动纹理化能力。 图一为原始模型, 图二图三为贴图后的模型。提示词: city, Realistic , cinematic , Front view ,Game scene graph 1、…

Kitex踩坑 [Error] KITEX: processing request error,i/o timeout

报错问题 2023/010/28 17:20:10.250768 default_server_handler.go:234: [Error] KITEX: processing request error, remoteService, remoteAddr127.0.0.1:65425, errordefault codec read failed: read tcp 127.0.0.1:8888->127.0.0.1:65425: i/o timeout 分析原因 Hert…

基于dockerfile搭建lnmp

目录 任务需求: 一、规划: 二、准备: 三、部署nginx容器(172.18.0.10): 1.编写Dockerfile构建镜像: 2.准备nginx配置文件: 3.构建镜像,启动nginx容器: 四…

【QT】事件分发器

event事件分发器,用于分发事件,在这里也可以做拦截,返回值boo。如果返回的是true代表拦截处理,不再向下分发。 展示事件拦截 上一段代码:【QT】鼠标常用事件-CSDN博客 代码 // 事件分发器 // 拦截鼠标按下 // QEven…

四十二、【进阶】

目录 1、覆盖索引 2、案例分析 (1)select * 查询 (2)使用字段查询 (3)性能差异原因 3、分析 (1)主键id查询 (2)覆盖索引 1、覆盖索引 简单点说&#x…

Nginx 部署多个安全域名,多个服务【工作记录】

以下是本人通过Docker 部署的Nginx挂载出来的文件目录 先看下 nginx.conf 配置文件内容:如下 ps:当前文件就是安装后的初始内容,无修改。主要关注最后一行 include /etc/nginx/conf.d/*.conf;表示引入其他目录下的.conf配置文件;…

【Arduino环境下驱动合宙esp32c3单片机基本外设】

【esp32c3基本外设驱动】 1. GPIO调试1.1 源码分享2.2 实验效果 2. ADC调试2.1 源码分享2.2 实验效果 3. WS2812驱动3.1 源码分享3.2 实验效果 4. 旋转编码器4.1 源码分享4.2 测试效果 5. SSD1306屏幕驱动5.1 源码分享5.2 测试效果 6. 双cpu同时工作测试6.1 源码分享6.2 测试效…

Arrays,Arrays重载的sort方法

Arrays -1的原因.因为返回正数不就是表示存在只能是负数 Arrays重载的sort方法 //这个方法只能给引用数据类型排序 //如果是基本数据类型需要转化为对应的包装类 public class arrays {public static void main(String[] args) {Integer arr[]{2,1,4,6,3,5,8,7,9};Arrays.s…

学习笔记3——JVM基础知识

学习笔记系列开头惯例发布一些寻亲消息 链接:https://baobeihuijia.com/bbhj/contents/3/196593.html JVM(Write Once,Run Anywhere) 以下是一些学习时有用到的资料,只学习了JVM的基础知识,对JVM整体进…

python爬虫—使用xpath方法进行数据解析

1. 背景信息 爬取安居客二手房源信息 URL地址:https://wuhan.anjuke.com/sale/?fromnavigation 2. 代码实现 import requests from lxml import etreeif __name__ __main__:# 1.指定URLurl "https://wuhan.anjuke.com/sale/?fromnavigation"# 2.U…

Apache Doris (四十八): Doris表结构变更-替换表

🏡 个人主页:IT贫道_大数据OLAP体系技术栈,Apache Doris,Clickhouse 技术-CSDN博客 🚩 私聊博主:加入大数据技术讨论群聊,获取更多大数据资料。 🔔 博主个人B栈地址:豹哥教你大数据的个人空间-豹哥教你大数据个人主页-哔哩哔哩视频 目录

在Jetpack Compose中使用Paging 3实现无限滚动

在Jetpack Compose中使用Paging 3实现无限滚动 本文将介绍在Jetpack Compose中进行分页加载。分页加载意味着一次只加载应用程序中的小数据块。 假设您在服务器上有大量数据,并且您希望在UI上显示这些数据。显然,您不希望一次性加载所有数据。您希望每次…

reactor (百万并发服务器) -- 1

为了从点滴开始,文章会先从一些基础socket去补充一些经常发生但是没有很深入去思考的细节。然后我们再开始去设计reactor的设计,可以选择跳过起过前面部分。 为了能从0开始去设计,测试,优化...整个过程会分为2-3篇文章输出&#x…

Redis队列Stream

1 缘起 项目中处理文件的场景: 将文件处理请求放入队列, 一方面,缓解服务器文件处理压力; 另一方面,可以根据文件大小拆分到不同的队列,提高文件处理效率。 这是Java开发组Leader佳汇提出的文件处理方案&a…

网络安全—小白自学

1.网络安全是什么 网络安全可以基于攻击和防御视角来分类,我们经常听到的 “红队”、“渗透测试” 等就是研究攻击技术,而“蓝队”、“安全运营”、“安全运维”则研究防御技术。 2.网络安全市场 一、是市场需求量高; 二、则是发展相对成熟…

QT实现用本地资源管理器来打开文件夹

QString path"文件夹路径";QDesktopServices::openUrl(QUrl("file:"path, QUrl::TolerantMode)); 在windows中QT编程,使用资源管理器来打开指定本地文件夹的方法: 第一种:使用Qprocess命令(相当于在cmd命令管…