本科毕业论文(设计)

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仿真相关推荐

  1. matlab仿真直流输电接地极系统,轻型高压直流输电系统的MATLAB仿真

    第 28 卷 第 6 期 2005 年 6 月 合肥 工 业 大 学 学 报 ( 自 然 科 学 版 ) JOU RNAL O F HEFE IUN IV ERSITY O F TECHNOLO GY ...

  2. matlab求解系统的差分方程图形,用MATLAB仿真散系统差分方程.doc

    用MATLAB仿真散系统差分方程 HEFEI UNIVERSITY 信号与系统项目设计报告 系 别 电子信息与电子工程系 题 目 项目第十题 专 业 电子信息工程 班 级 11电子信息工程(2)班 小 ...

  3. matlab仿真采样时间,系统中的采样时间 - MATLAB Simulink - MathWorks 中国

    纯离散系统 纯离散系统完全由离散模块组成,可以使用固定步长或可变步长求解器进行建模.要对离散系统进行仿真,需要仿真器在每个采样时间点执行一个仿真步.对于多速率离散系统 - Simulink® 以不同速 ...

  4. matlab用于系统框图建模的函数,MATLAB产品家族中文

    MATLAB MATLAB 是一种用于算法开发.数据可视化.数据分析以及数值计算的高级技术计算语言和交互式环境.使用 MATLAB,您可以较使用传统的编程语言(如 C.C++ 和 Fortran)更快 ...

  5. matlab励磁系统接不了,基于MATLAB的同步发电机励磁系统仿真.pdf

    基于MATLAB的同步发电机励磁系统仿真,发电机matlab仿真,同步发电机励磁系统,同步发电机励磁控制,无刷励磁同步发电机,同步发电机励磁,同步发电机励磁方式,matlab同步电机仿真,同步发电机励 ...

  6. 利用MATLAB进行系统时域分析,实验二 利用matlab进行系统的时域分析

    实验二 利用matlab进行系统的时域分析 实验二 利用MATLAB进行系统的时域分析 1.实验目的 在理论学习的基础上,通过本实验熟悉LTI连续时间系统的时域分析方法, 熟悉系统的零输入响应.零状态 ...

  7. Hypermesh+dyna整车碰撞模型( 仿真+试验对标)-加计算分析报告 碰撞仿真与碰撞试验的对标模型及报告

    Hypermesh+dyna整车碰撞模型( 仿真+试验对标)-加计算分析报告 国外专业机构做的碰撞仿真与碰撞试验的对标模型及报告, 基本款与轻量化款对比,数据详实,不可多得的好资料,懂的工程师联系我. ...

  8. matlab 减震,系统减震论文,关于基于MATLAB汽车减震系统相关参考文献资料-免费论文范文...

    导读:这是一篇与系统减震论文范文相关的免费优秀学术论文范文资料,为你的论文写作提供参考. (东北大学秦皇岛分校 河北秦皇岛 066000) 摘 要:汽车的减震系统性能直接决定汽车的舒适度.本研究采用弹 ...

  9. 计算机仿真实训系统的企业排名,闵行区虚拟仿真实训平台企业-哪家不错

    闵行区虚拟仿真实训平台企业-哪家不错,尤其是对于现实世界被VR所改变后的景象,描绘之精彩令人叫绝.VR(VirtualReality)技术,中文即虚拟现实,就是通过一系列技术手段来制造虚拟场景并模拟各 ...

  10. matlab一直系统函数画脉冲响应,【 MATLAB 】使用 impz 函数计算并画出脉冲响应

    这篇博文在于讨论 impz 函数的使用. 我在 MATLAB 中查看 impz 的帮助文档时,始终看的不是太明白这个函数的使用,于是我根据一个例子,对这个函数做了一点分析,解决了一些我的疑惑,记录于此 ...

最新文章

  1. 独家 | 如何让情绪智力驱动创新?
  2. oracle监听 客户 实例,oracle 数据库实例 监听
  3. 吕帅 计算机学院,刘磊-吉林大学计算机科学与技术学院
  4. java实现人字拼,人字拼地板拼法大全
  5. dubbo k8s 服务发现_将Dubbo微服务迁移到k8s集群环境中前的思考与落地
  6. 如何操作反射中构造方法、属性和普通方法?
  7. php 图片裁剪保存,PHP图片裁剪与缩放示例(无损裁剪图片)
  8. node.js 安装详细步骤教程
  9. linux命令控制器在哪打开,在Linux上使用PS4控制器的两种方法:安装Steam和DS4DRV
  10. Eclipse中添加外部Web浏览器
  11. 火车头采集html5游戏,火车头采集工具Wed发布配置方法
  12. 组成原理 - 网卡与HBA卡的区别
  13. html头像裁剪仿微信,仿微信群头像css布局
  14. 雷达传统的干扰方式(二次更新)
  15. 20-Netty TCP 粘包和拆包及解决方案
  16. Unity3D教程:手游开发常用排序算法 -下
  17. python录制视频和声音_moviepy 视频编辑添加MP3没声音的问题
  18. macOS-允许打开任何来源的程序
  19. r语言 相关性作图_R语言:多个基因的相关性分析与展示
  20. Python100行-贪吃蛇小游戏

热门文章

  1. BM3D算法相关笔记
  2. MDCC 2016参会总结
  3. 使用Tiled编辑铁锈战争自定义地图
  4. 越过山丘才发现无人等候 倒计时84
  5. 【Cubieboard2】配置编译内核支持SPI全双工通信驱动
  6. rstudio中读取数据_用RStudio导入数据
  7. paip 突破 网站 手机 验证码 的 破解 总结
  8. 数据分析新人如何面对繁杂且突然的数据需求
  9. AIOT:基于智能家居谈AIOT
  10. 第五篇 《小强升职记》