MIMO系统的信道容量分析 及Matlab仿真 下载本文

本科毕业论文(设计)

第 21 页

correlation matrix approach [J] . IEEE Communications Lett ers, 2002, 16 (1) : 34- 35. [10] D Gesbert , H Bolcskei, D A Gore, et al. Performance evaluat ion for scattering MIMO channel models[A] . 34th Asilomar Conference on Signal,System and Comput ers[C] .USA. Nov. 2000, 1(1) : 748- 752.

[11] 程健, 陈明, 程时昕.无线通信领域中空时编码技术[J].电路与系统学报, 2002, 7( 1) : 66- 71.

[12] 孙丹,张晓光.MIMO系统信道容量研究[D].徐州:无线通信实验室,2009. [13] P.Driessen and G.Foschini,On the capacity of multiple input-multiple output wireless channels:a geometric interpretation,IEEETrans.Commun.,vol.47,no.2, Feb.1999.

附 录

1.信道容量计算程序 程序代码如下:

SNR=0; %定义信噪比为0dB A=10^(SNR/10); %信噪比的单位转换关系式

Im=eye(1); %这里的\为发射天线和接收天线两者数

目少的根数

sum=0; %初始化sum为0 for i=1:1000 %程序循环1000次 H=randn(1,4); %产生高斯信号源 Q=H*H'; %求威沙特矩阵

sum=sum+log2(det(Im+A*Q/4)); 00次循环得到的信道容量总和 end

C=sum/1000 00次循环的平均值

本科毕业论文(设计)

第 22 页

2.发射天线数nT为4,接收天线nR逐渐增大,4种不同信噪比下的信道容量仿真程序

程序代码如下:

x=[1,2,3,4,5,6,7,8,9,10,20,30];

y1=[0.9267,1.7682,2.5436,3.2482,3.8977,4.5062,5.0686,5.5634,6.0534,6.5123,9.4206, 11.0734];

plot(x,y1 ,'r') hold on

y2=[1.8720,3.5384,5.0226,6.3022,7.4272,8.4336,9.2894,10.0768,10.7693,11.4027, 15.1084,17.0578];

plot(x,y2 ,'g') hold on

y3=[3.1524,5.9879,8.3914,10.4579,12.1479,13.5424,14.6978,15.6897,16.5322,17.0927, 20.9439,22.4942];

plot(x,y3 ,'b') hold on

y4=[4.6711,8.9050,12.5479,15.4908,17.6953,19.4381,20.7939,21.9145,22.8506,23.6258,27.0859,28.8846];

plot(x,y4 ,'k'),grid on xlabel('接收天线数nR'); ylabel('容量/(bit/s/Hz)');

3.接收天线nR为4,发射天线数nT逐渐增大,4种不同信噪比下的信道容量仿真程序 程序代码如下:

x=[1,2,3,4,5,6,7,8,9,10,20,30];

y1=[2.1177,2.7428,3.0877,3.2612,3.3472,3.4298,3.5241,3.5681,3.5964,3.6174,3.7560, 3.7899];

plot(x,y1 ,'r') hold on

y2=[3.4491,4.9309,5.9981,6.4075,6.6011,6.7730,6.9885,7.1010,7.2006,7.3069,7.6129,

本科毕业论文(设计)

7.6546];

plot(x,y2 ,'g') hold on

第 23 页

y3=[5.0124,7.7014,9.5183,10.5682,11.1156,11.5987,11.9483,12.1547,12.3633,12.5237, 13.0893,13.3259];

plot(x,y3 ,'b') hold on

y4=[6.6068,10.7773,13.7736,15.5987,16.5676,17.2934,17.7721,18.1088,18.3576, 18.5360,19.2798,19.4907];

plot(x,y4 ,'k'),grid on xlabel('发射天线数nT'); ylabel('容量/(bit/s/Hz)');