文档介绍：MIMO信道容量计算
实验一： MIMO信道容量计算
实验学时：3
实验类型：〔演示、验证、综合、设计、√研究〕
实验要求：〔√必修、选修〕
一、实验目的
通过本实验的学习，理解和掌握信道容量的概念和物理意义；了解4，d52=，d62=，d72=，d82=，d92=，d102=。
1、假设输入的信号总功率P=5
程序
L=input('子信道数目 L='); %输入子信道数目
d=input('各子信道噪声方差 d='); %输入各信道的噪声方差
P=input('输入信号总功率 P='); %输入信号总功率
D=sum(d); %各子信道的方差和
Pav=(P+D)*; %求平均输出频率
i=1;
p=ones(1,L); %功率初始值
boolean=1;
c=0; %关闭的信道数量
while(boolean)
while(i<=L)
if(Pav<d(i)& p(i)~=0) %找出不符合的子信道并赋值为0
p(i)=0;
D=D-d(i);
c=c+1;
end
i=i+1;
end
if(c==0) %当所有信道满足时退出循环
break;
end
if(i==L+1) %每个信道和Pav比照一遍后初始化各值，并重新计算Pav
i=1;
L=L-c;
c=0;
Pav=(P+D)*;
end
end
L=length(ai);
for i=1:L
if(p(i)~=0)
p(i)=Pav-d(i); %计算各子信道分配功率
end
end
C=0;
for i=1:L
C=C+1/2*log2(1+p(i)/d(i));
end
p
C输出
>> WF
子信道数目 L=10
各子信道噪声方差 d=[ ]
输入信号总功率 P=5
p =
Columns 1 through 8
Columns 9 through 10
C =
>> WF
子信道数目 L=10
各子信道噪声方差 d=[ ]
输入信号总功率 P=1
p =
Columns 1 through 8
0 0 0 0
Columns 9 through 10
0 0
C =
二、实验题目
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); %各子信道的方