实验一： MIMO信道容量计算
实验学时：3
实验类型：（演示、验证、综合、
设计
领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计
、√研究）
实验要求：（√必修、选修）
一、实验目的
通过本实验的学习，理解和掌握信道容量的概念和物理意义；了解多天线系统信道容量的计算方法；采用计算机编程实现经典的注水算法。
二、实验内容
MIMO信道容量；
注水算法原理；
采用计算机编程实现注水算法。
三、实验组织运行要求
以学生自主训练为主的开放模式组织教学
四、实验条件
（1）微机
（2）MATLAB编程工具
五、实验原理、方法和手段
MIMO（MIMO，Multiple Input Multiple Output）技术利用多根天线实现多发多收，充分利用了空间资源，在有限的频谱资源上可以实现高速率和大容量，已成为4G通信系统以及未来无线通信系统的关键技术之一。
图1平坦衰弱MIMO信道模型
1．MIMO信道模型
MIMO指多输入多输出系统，当发送信号所占用的带宽足够小的时候，信道可以被认为是平坦的，即不考虑频率选择性衰落。平坦衰弱的MIMO信道可以用一个的复数矩阵描述：
（1）
其中为发送端天线数，为接收端天线数，的元素表示从第根发射天线到第根接收天线之间的空间信道衰落系数。
窄带MIMO信道模型（如图1所示）可以描述为：
（2）
其中，为发送信号；为接收信号；为加性高斯白噪声。
2．MIMO信道容量
假设服从均值为0，协方差为单位阵的复高斯分布。根据信道容量的定义，可以
证明
住所证明下载场所使用证明下载诊断证明下载住所证明下载爱问住所证明下载爱问
当服从高斯分布时，达到MIMO信道容量。令的协方差矩阵为，则MIMO信道容量可表示为：
（3）
其中上标‘H’表示复共轭，为单位阵，det表示取行列式。表示单位带宽下的MIMO信道传输速率，单位为Nat/sec。
发射机的传输功率可以表示为：
其中，表示向量的模，Tr表示取矩阵的迹，E表示求期望。
假设发射机的最大传输功率为，则功率约束下的MIMO信道容量计算问题可以描述为：
（4）
即在功率约束下找到的最佳分布使得信道容量最大。此处表示为半正定矩阵。
3．奇异值分解及问题转化
将信道矩阵进行奇异值分解，即可分解为，其中和为酉矩阵满足
，
为的矩形对角矩阵。
利用恒等式和酉矩阵性质，得到
（5）
利用恒等式和酉矩阵性质，得到
（6）
定义，利用（5）和（6），问题（4）可写为：
（7）
经过变量替换，令，问题（7）（i.e., (4)）等价为
（8）
根据Hadamard不等式（参考《信息论基础》Thomas M. Cover，定理16.8.2），问题（8）的最优解必为对角阵。因此问题（8）可简化为
（9）
其中，r为（或者）的秩，和分别为和的对角元素。
4．注水（water-filling）算法 求解问题（9）
对问题（9）使用拉格朗日乘子法，可知最优解必须满足：
其中为拉格朗日乘子。为求得，需先确定，它满足
(10)
注意左边为的递增函数，因此满足（10）的唯一。
找到满足（10）的的方法可形象地称为注水（warter-filling）。这是因为，我们可以将看做是第片区域的水平线，然后对整个区域注水，使其具有深度，如图1所示。所需总水量为，不断注水，直至总水量为,。第个区域的水位深度即为最优的。
图2 注水算法图示。每片区域的高度为。注水总量为，对整个区域
注水使其高度达到。每片区域上水的高度（阴影部分所示）即为最优的。
六、实验数据及具体步骤
1）随机产生,每个信道系数服从均值为0，方差为1的高斯分布（对Rayliegh分布）
2) 奇异值分解（利用matlab SVD函数）
3）计算
4）注水法求解问题（9）
5）计算，根据（3）计算容量（或者根据问题（9）的最优值确定最大容量）。
七、实验程序和结果
一、有一并联高斯加性信道，各子信道的方差为d12=0.1，d22=0.2，d32=0.3，d42=0.4，d52=0.5，d62=0.6，d72=0.7，d82=0.8，d92=0.9，d102=1.0。
1、若输入的信号总功率P=5
程序
L=input('子信道数目 L='); %输入子信道数目
d=input('各子信道噪声方差 d='); %输入各信道的噪声方差
P=input('输入信号总功率 P='); %输入信号总功率
D=sum(d); %各子信道的方差和
Pav=(P+D)*1.0/L; %求平均输出频率
i=1;
p=ones(1,L); %功率初始值
boolean=1;
c=0; %关闭的信道数量
while(boolean)
while(i<=L)
if(Pav
> WF
子信道数目 L=10
各子信道噪声方差 d=[0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0]
输入信号总功率 P=5
p =
Columns 1 through 8
0.9500 0.8500 0.7500 0.6500 0.5500 0.4500 0.3500 0.2500
Columns 9 through 10
0.1500 0.0500
C =
6.0661
>> WF
子信道数目 L=10
各子信道噪声方差 d=[0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0]
输入信号总功率 P=1
p =
Columns 1 through 8
0.4000 0.3000 0.2000 0.1000 0 0 0 0
Columns 9 through 10
0 0
C =
2.3514
二、实验题目
clear all;
nt=[2 4 8];
nr=[2 4 8];
v=0; %给定1/v初始值（也就是这里的v）
for count=1:3
y=0;
for PTdb=0:5:30
y=y+1;
PT=10^(PTdb/10); %把PT转化为10进制
for k=1:1000
H=sqrt(1/2)*complex(randn(nr(count),nt(count)),randn(nr(count),nt(count))); %随机产生H ,每个信道系数服从均值为0，方差为1的高斯分布（对Rayliegh分布）(sqrt(1/2)为了使模也为1)
[U,D,VH]=svd(H); %奇异值分解 H=U*D*VH
A=D'*D; %求DH*D
ai=sort((1./diag(A))'); %求ai并按从小到大排序
L=length(ai);
AI=sum(ai); %各子信道的方差和
Pav=(PT+AI)/L; %求平均输出频率
i=1;
p=ones(1,L); %功率初始值
boolean=1;
ct=0; %关闭的信道数量
while(boolean)
while(i<=L)
if(Pav
本文档为【MIMO信道容量计算】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483