MIMO信道容量推导(奇异值分解法)

2024-02-05 10:38

本文主要是介绍MIMO信道容量推导(奇异值分解法),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

[转]MIMO信道容量推导

(奇异值分解法)
Massive MIMO,目前是5G的一项关键技术。那么求解它的信道容量,对于我们研究它的属性尤为重要。今天,我们就求算一下信道容量。
在通信系统中我们学到信道容量的一个计算式:
C = B l o g 2 ( 1 + S N R ) C=Blog_2(1+SNR) C=Blog2(1+SNR)
其中 B B B是带宽, S N R SNR SNR是信噪比, C C C是信道容量。
那么其实我们就是需要求解 S N R SNR SNR
假设M个接收天线,N个发送天线,那么MIMO模型的表达式如下:
y = Hx + n \textbf{y}=\textbf{H}\textbf{x}+\textbf{n} y=Hx+n
其中, y \textbf{y} y M × 1 M\times 1 M×1输出的信息, H \textbf{H} H是一个只能被接收机感知的 M × N M\times N M×N矩阵, x \textbf{x} x是输入信息的 N × 1 N\times 1 N×1列向量, n \textbf{n} n是AWGN噪声组成的 M × 1 M\times 1 M×1列向量。

先了解一下奇异值分解(SVD):
SVD也是对矩阵进行分解,但是和特征分解不同,SVD并不要求要分解的矩阵为方阵。那么我们定义矩阵H的SVD为:
H = UD V H \textbf{H}=\textbf{U}\textbf{D}\textbf{V}^H H=UDVH
其中 U \textbf{U} U是一个 M × M M\times M M×M的矩阵, D \textbf{D} D是一个 M × N M\times N M×N的矩阵,除了主对角线上的元素以外全为0,主对角线上的每个元素都称为奇异值, V \textbf{V} V是一个 N × N N\times N N×N的矩阵。 U \textbf{U} U V \textbf{V} V都是酉矩阵,即满足:
I = X T X \textbf{I}=\textbf{X}^T\textbf{X} I=XTX

那么,考虑使用奇异值分解法处理 H \textbf{H} H矩阵,得到:
H = UD V H \textbf{H}=\textbf{U}\textbf{D}\textbf{V}^H H=UDVH
y = UD V H x + n \textbf{y}=\textbf{U}\textbf{D}\textbf{V}^H\textbf{x}+\textbf{n} y=UDVHx+n
U H ( y − n ) = D V H x \textbf{U}^H(\textbf{y}-\textbf{n})=\textbf{D}\textbf{V}^H\textbf{x} UH(yn)=DVHx
令: y ^ = U H y \hat\textbf{y}=\textbf{U}^H\textbf{y} y^=UHy x ^ = V H x \hat\textbf{x}=\textbf{V}^H\textbf{x} x^=VHx n ^ = U H n \hat\textbf{n}=\textbf{U}^H\textbf{n} n^=UHn
则有,
y ^ = D x ^ + n ^ \hat\textbf{y}=\textbf{D}\hat\textbf{x}+\hat\textbf{n} y^=Dx^+n^
令奇异值为 λ i \sqrt{\lambda_i} λi ,最多 r a n k ( H ) = min ⁡ ( M , N ) rank(\textbf{H})=\min(M,N) rank(H)=min(M,N)个非零奇异值,那么接收信号为:
y ^ i = λ i x ^ i + n ^ i , i = 1 , 2 , . . . , r a n k \hat{y}_i=\sqrt{\lambda_i}\hat{x}_i+\hat{n}_i,i=1,2,...,rank y^i=λi x^i+n^i,i=1,2,...,rank y ^ i = n ^ i , i = r a n k + 1 , . . . , M \hat{y}_i=\hat{n}_i,i=rank+1,...,M y^i=n^i,i=rank+1,...,M
由上述可知,MIMO其实是rank个平行的等效子信道组成的。幅度增益就是奇异值,那么单个信道的信噪比就是 p i σ 2 \frac{p_i}{\sigma^2} σ2pi,其中信息功率为平均下来每个天线的接收功率,即 p i = P T M p_i=\frac{P_T}{M} pi=MPT,而后对单个信道进行叠加,得到:
C = B log ⁡ 2 ∏ i = 1 r a n k ( 1 + λ i P T M σ 2 ) C=B\log_2\prod \limits_{i=1}^{rank}\left(1+\frac{{\lambda_i}P_T}{M\sigma^2}\right) C=Blog2i=1rank(1+Mσ2λiPT)
即最终MIMO的信道容量。

这篇关于MIMO信道容量推导(奇异值分解法)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/680649

相关文章

uva 10014 Simple calculations(数学推导)

直接按照题意来推导最后的结果就行了。 开始的时候只做到了第一个推导,第二次没有继续下去。 代码: #include<stdio.h>int main(){int T, n, i;double a, aa, sum, temp, ans;scanf("%d", &T);while(T--){scanf("%d", &n);scanf("%lf", &first);scanf

特征值分解(EVD)和奇异值分解(SVD)—应用于图片压缩

特征值分解(EVD)和奇异值分解(SVD)—应用于图片压缩 目录 前言 一、特征值分解 二、应用特征值分解对图片进行压缩 三、矩阵的奇异值分解 四、应用奇异值分解对图片进行压缩 五、MATLAB仿真代码 前言         学习了特征值分解和奇异值分解相关知识,发现其可以用于图片压缩,但网上没有找到相应代码,本文在学习了之后编写出了图片压缩的代码,发现奇异值分

【无线通信发展史⑧】测量地球质量?重力加速度g的测量?如何推导单摆周期公式?地球半径R是怎么测量出来的?

前言:用这几个问答形式来解读下我这个系列的来龙去脉。如果大家觉得本篇文章不水的话希望帮忙点赞收藏加关注,你们的鼓舞是我继续更新的动力。 我为什么会写这个系列呢? 首先肯定是因为我本身就是一名从业通信者,想着更加了解自己专业的知识,所以更想着从头开始了解通信的来源以及在每一个时代的发展进程。 为什么会从头开始写通信? 我最早是学习了中华上下五千年,应该说朝代史,这个算个人兴趣,从夏

CF #278 (Div. 2) B.(暴力枚举+推导公式+数学构造)

B. Candy Boxes time limit per test 1 second memory limit per test 256 megabytes input standard input output standard output 题目链接: http://codeforces.com/contest/488/problem/B There

HDU2524(规律推导)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2524 解题思路: 暴力推出矩阵,以n = 2 , m = 4为例: 1 3  6  10 3 9 18 30 可以发现第一行和第一列都是有规律的,彼此相差2、3、4·····,其他元素为相应行第一个元素乘以第一列元素的积。预处理之后,我们O(1)就可以输出g[n][m]的值。 另外,

Ural1209(数学推导)

题目链接:点击打开链接 解题思路: 此题甚好。推导公式,首先观察序列110100100010000·····,我们把为1的下标单独拿出来看。依次为1、2、4 、7、 11·····,可以分解为1+(0) 、1+(0+1)、1+(0+1+2)、1+(0+1+2+3)、1+(0+1+2+3+4),可以推导出规律1 + x * (x - 1) / 2。 那么对于每个n,我们只要判断是否存在x

Python中如何实现列表推导式(List Comprehension)

Python中的列表推导式(List Comprehension)是一种简洁且高效的方式来创建列表。它不仅让代码更加简洁,而且通常比使用循环和条件语句生成列表更快。列表推导式的基本形式允许你从现有的列表或其他可迭代对象中创建新的列表,同时应用过滤和转换操作。下面我将详细解释列表推导式的概念、基本语法、高级用法以及其在实际应用中的优势。 一、列表推导式的基本概念 列表推导式是Python中的一种

对极约束及其性质 —— 公式详细推导

Title: 对极约束及其性质 —— 公式详细推导 文章目录 前言1. 对极约束 (Epipolar Constraint)2. 坐标转换 (Coordinate Transformations)3. 像素坐标 (Pixel Coordinates)4. 像素坐标转换 (Transformations of Pixel Coordinates)5. 本质矩阵 (Essential Matr

连分数因子分解法——C语言实现

参考网址:连分数分解法寻找整数的因子(Python)-CSDN博客 大数运算:C语言实现 大数运算 加减乘除模运算 超详细_64编程 加减乘除取模 复杂运算-CSDN博客 ‌连分数因子分解法‌是一种用于大整数因子分解的算法,它是计算数论中的一个重要方法。连分数因子分解法通过寻找x2≡y2 (mod p)x2≡y2 (mod p)的形式来分解N。具体来说,这种方法涉及到计算N的简单连分数展开,并

时序预测|变分模态分解-双向时域卷积-双向门控单元-注意力机制多变量时间序列预测VMD-BiTCN-BiGRU-Attention

时序预测|变分模态分解-双向时域卷积-双向门控单元-注意力机制多变量时间序列预测VMD-BiTCN-BiGRU-Attention 文章目录 一、基本原理1. 变分模态分解(VMD)2. 双向时域卷积(BiTCN)3. 双向门控单元(BiGRU)4. 注意力机制(Attention)总结流程 二、实验结果三、核心代码四、代码获取五、总结 时序预测|变分模态分解-双向时域卷积