MIMO信道容量仿真MATLAB,MIMO系统的信道容量分析 及Matlab仿真
本科毕业论文(设计)
2.3.2 MIMO系统信道容量推导
1、 MIMO系统的瞬时信道容量的推导
第 9 页
这部分我们将给出MIMO信道容量的一般性表达。根据上面的信道模型,我们可以得到收发信号关系为:
y?H x ? n (2.10)
首先假设信道的加性噪声n是服从协方差为Rnn = E{nnH}的零均值循环对称复高斯分布(ZMCSCG),即:n~CN(0N,N0IN)并且n与x之间是不相关的。发送信号x是服从零均值、协方差为RXX ={xxH }、概率密度函数为fs(x)的分布,总的发送功率限制为Tr{Rxx}=P,其中Ex为在一个符号周期内总的发送能量。
在下面的推导过程中我们假设信道矩阵H在接收端已经完全已知,但是它是随机的,因此我们可以得到瞬时信道容量为:
C(H)?maxI(x;y) (2.11)
fs(x)其中I(x;y)为在已知信道H的情况下输入x与输出y之间的互信息量,有:
I(x;y)?H(y)?H(y|x) (2.12)
H(y)是y的差分熵,H(y|x)是给定x条件下y的差分熵,由于发送信号与噪声之间是独立的,因此有H(y|x)=H(n),所以式(2.12)可以重新写为:
I(x;y)?H(y)?H(n) (2.13)
因为接收信号的协方差矩阵为:
Ryy?E{yyH}?HRxxHH?N0IN
(2.14)
对于输出信号y的差分熵,根据Neeser的分析,在给定协方差矩阵Ryy的条件下,H(y)只有在y是也服从ZMCSCG分布情况下才可以达到最大值,所以发送信号x也应该服从ZMCSCG分布,此时的y与n的差分熵分别为:
(2.15) H(y)?log2{det(?eRyy)} bits/s/ H (2.16) H(n)?log2{det(?eN0IN)} bits/s/ H 本科毕业论文(设计)
所以我们可以得到信道瞬时交互信息I(x;y)为:
I(x;y)?log2{det[IN?
第 10 页
1 (2.17) HRxxHH]} bits/s/ HN0我们知道信道容量是最大的输入输出交互信息,所以(2.11)可以重新写为:
C(H)?maxlog2{det[IN?
Tr{Rxx}?p1HRxxHH]} bits/s/Hz (2.18) N0上面得到的瞬时容量值是随不同的信道条件而不断变化的,由于信道矩阵H的随机性,所以我们可以知道瞬时信道容量也是随机的,是一个随机变量。
2、利用矩阵理论推导MIMO系统容量
对于MIMO无线信道,信道是极其复杂的。因此原始的信道矩阵也就显得复杂,不便于分析,而且一般矩阵不经过处理计算行列式很困难。这就自然想到在信源端对发射信号做某种预处理,使得经过预处理的信号经过的信道变得简单易分析,而且具体实现也变得简单。对于信道矩阵来说,对角矩阵是最简单的,所以自然就想到把信道矩阵分解,利用矩阵理论中的奇异值分解可以达到这种目的。下面就矩阵的奇异值分解来计算MIMO的信道容量。首先,假设信道矩阵在发射端为未知,在接收端为已知。
由奇异值分解(SVD)理论,任何一个nR×nT矩阵H可以写成
H?UDVH (2.19)
式中,D是nR×nT非负对角矩阵;U和V分别是 nR×nR和 nT×nT 的酉矩阵。则有UUH=InR和VVH=InT,其中InR和InT分别是nR×nR和nT×nT单位阵。D的对角元素是矩阵HHH的特征值的非负平方根。HHH的特征值(用λ表示)定义为
HHHy??y,y?0
(2.20)
式中,y是与λ相对应的nR×1维矢量,称为特征矢量。
特征值的非负平方根也称为H的奇异值,而且U的列矢量是HHH的特征矢量,V的列矢量是HHH的特征矢量。把把(2.19)代入(2.10),可以得到接收矢量r
r?UDVHx?n
(2.21)
引入下列变换:
本科毕业论文(设计)
第 11 页
r??UHr
x??VHx (2.2
2)
n??UHn
U和V是可逆的。显然,式(2.22)中定义的矩阵r 、x和n与相应矩阵的乘积仅有一个缩放比例的效果。矢量n′是一个零均值高斯随机变量,其实部和虚部独立同分布。这样,前面讨论的信道与下式所描述的信道是等价的。
r??Dx??n (2.23)
矩阵HHH的非零特征值的数量等于矩阵H的秩,用r表示。对nR×nT矩阵H,秩的最大值为m=min(nR,nT),也就是说,至多有m个奇异值是非零的。用?i表示H的奇异值。将?i代入式(2.23),得到接收信号元素为
???ri??ixi?ni?i?1,2,?,r?
??ri?ni?i?r?1,r?2,?,nR? (2.24)
?式(2.24)显示,接收元素ri(i?r?1,r?2,?,nR)并不依赖于发射信号,即信道增
?益是零。另一方面,接收元素ri(i?1,2?,r)仅仅取决于发射元素Xi’。因此,可以认为,通过(2.23)得到的等效MIMO信道是由r个去耦平行子信道组成的。为每个子信道分配的矩阵H的奇异值,相当于信道幅度增益。因此,信道功率增益等于矩阵HHH的特征值。例如,如果nT > nR,由于H的秩不可能比nR高,那么式(2.24)显示了在等效的MIMO信道中,最多有nR个非零增益子信道。
另一方面,如果nR > nT,在等效的MIMO信道中,最多有nT个非零增益子信道。特征值谱是对MIMO信道的一种描述方式,适用于对最佳发射路径进行估计。
由式(2.22),可以导出信号r′、x′和n′的协方差矩阵和它们的迹
Rr?r??UHRrrU
Rx?x??VHRxxV
(2.25)
本科毕业论文(设计)
第 12 页
Rn?n??UHRnnU
r?)?tr(rrR) tr(r?R?x?)?tr(xxR) (2.26) tr(xR?n?)?tr(nnR) tr(nR以上关系显示,r'、x' 和n' 的协方差矩阵有相等的对角元素和,从而有相等的功率;而对于原始信号和r、x和n,它们是各不相等的。
考虑到(2.22)所描述的等价MIMO信道模型中,子信道是去耦的,因此其容量可以直接相加。假设在等效MIMO信道中,每根天线的发射功率为P/nT,运用香农公式,可以估算出总的信道容量(用C表示)为:
?P? C?W?log2?1?ri? (2.27)2???i?1
r式中,W是每个子信道的带宽;Pri是在第i个子信道中接收的信号功率,由下式给出:
Pri?
?iP (2.28) nT式中,?i是信道矩阵H的奇异值。因此信道容量可以写成:
r???iP??iP??? (2.29) C?W?log2?1??Wlog1?2??2?2?nT??i?li?1??nT??
r下面说明信道容量是如何与信道矩阵H相关的。假定m=min(nR,nT),(2.20)定义了特征值-特征矢量的关系,可重新写为:
??Im?Q?y?0
y?0 (2.30)
式中,Q是威沙特(Wishart)矩阵,定义为:
n?nT?HHHQ??H R (2.31)
n?nRT?HH
即当且仅当?Im?Q是奇异矩阵时,λ是Q的一个特征值。因此?Im?Q的行列式必定为零,即
MIMO信道容量仿真MATLAB,MIMO系统的信道容量分析 及Matlab仿真相关推荐
- matlab仿真直流输电接地极系统,轻型高压直流输电系统的MATLAB仿真
第 28 卷 第 6 期 2005 年 6 月 合肥 工 业 大 学 学 报 ( 自 然 科 学 版 ) JOU RNAL O F HEFE IUN IV ERSITY O F TECHNOLO GY ...
- matlab求解系统的差分方程图形,用MATLAB仿真散系统差分方程.doc
用MATLAB仿真散系统差分方程 HEFEI UNIVERSITY 信号与系统项目设计报告 系 别 电子信息与电子工程系 题 目 项目第十题 专 业 电子信息工程 班 级 11电子信息工程(2)班 小 ...
- matlab仿真采样时间,系统中的采样时间
- MATLAB Simulink
- MathWorks 中国
纯离散系统 纯离散系统完全由离散模块组成,可以使用固定步长或可变步长求解器进行建模.要对离散系统进行仿真,需要仿真器在每个采样时间点执行一个仿真步.对于多速率离散系统 - Simulink® 以不同速 ...
- matlab用于系统框图建模的函数,MATLAB产品家族中文
MATLAB MATLAB 是一种用于算法开发.数据可视化.数据分析以及数值计算的高级技术计算语言和交互式环境.使用 MATLAB,您可以较使用传统的编程语言(如 C.C++ 和 Fortran)更快 ...
- matlab励磁系统接不了,基于MATLAB的同步发电机励磁系统仿真.pdf
基于MATLAB的同步发电机励磁系统仿真,发电机matlab仿真,同步发电机励磁系统,同步发电机励磁控制,无刷励磁同步发电机,同步发电机励磁,同步发电机励磁方式,matlab同步电机仿真,同步发电机励 ...
- 利用MATLAB进行系统时域分析,实验二 利用matlab进行系统的时域分析
实验二 利用matlab进行系统的时域分析 实验二 利用MATLAB进行系统的时域分析 1.实验目的 在理论学习的基础上,通过本实验熟悉LTI连续时间系统的时域分析方法, 熟悉系统的零输入响应.零状态 ...
- Hypermesh+dyna整车碰撞模型( 仿真+试验对标)-加计算分析报告 碰撞仿真与碰撞试验的对标模型及报告
Hypermesh+dyna整车碰撞模型( 仿真+试验对标)-加计算分析报告 国外专业机构做的碰撞仿真与碰撞试验的对标模型及报告, 基本款与轻量化款对比,数据详实,不可多得的好资料,懂的工程师联系我. ...
- matlab 减震,系统减震论文,关于基于MATLAB汽车减震系统相关参考文献资料-免费论文范文...
导读:这是一篇与系统减震论文范文相关的免费优秀学术论文范文资料,为你的论文写作提供参考. (东北大学秦皇岛分校 河北秦皇岛 066000) 摘 要:汽车的减震系统性能直接决定汽车的舒适度.本研究采用弹 ...
- 计算机仿真实训系统的企业排名,闵行区虚拟仿真实训平台企业-哪家不错
闵行区虚拟仿真实训平台企业-哪家不错,尤其是对于现实世界被VR所改变后的景象,描绘之精彩令人叫绝.VR(VirtualReality)技术,中文即虚拟现实,就是通过一系列技术手段来制造虚拟场景并模拟各 ...
- matlab一直系统函数画脉冲响应,【 MATLAB 】使用 impz 函数计算并画出脉冲响应
这篇博文在于讨论 impz 函数的使用. 我在 MATLAB 中查看 impz 的帮助文档时,始终看的不是太明白这个函数的使用,于是我根据一个例子,对这个函数做了一点分析,解决了一些我的疑惑,记录于此 ...
最新文章
- 独家 | 如何让情绪智力驱动创新?
- oracle监听 客户 实例,oracle 数据库实例 监听
- 吕帅 计算机学院,刘磊-吉林大学计算机科学与技术学院
- java实现人字拼,人字拼地板拼法大全
- dubbo k8s 服务发现_将Dubbo微服务迁移到k8s集群环境中前的思考与落地
- 如何操作反射中构造方法、属性和普通方法?
- php 图片裁剪保存,PHP图片裁剪与缩放示例(无损裁剪图片)
- node.js 安装详细步骤教程
- linux命令控制器在哪打开,在Linux上使用PS4控制器的两种方法:安装Steam和DS4DRV
- Eclipse中添加外部Web浏览器
- 火车头采集html5游戏,火车头采集工具Wed发布配置方法
- 组成原理 - 网卡与HBA卡的区别
- html头像裁剪仿微信,仿微信群头像css布局
- 雷达传统的干扰方式(二次更新)
- 20-Netty TCP 粘包和拆包及解决方案
- Unity3D教程:手游开发常用排序算法 -下
- python录制视频和声音_moviepy 视频编辑添加MP3没声音的问题
- macOS-允许打开任何来源的程序
- r语言 相关性作图_R语言:多个基因的相关性分析与展示
- Python100行-贪吃蛇小游戏