channelcapacity_batch.m

本人详细分析后的mimo信道容量的仿真代码

%% Channel Capacity of MIMO 
clear;
clc;
results=zeros(51,5);

Re_CorrelationMatrix=zeros(4,4,51);

for mm=1:4
    Re_CorrelationMatrix(mm,mm,:)=1;
end

Coefficient_1_2=zeros(1,51);
Coefficient_1_3=zeros(1,51);
Coefficient_1_4=zeros(1,51);
Coefficient_4_1=zeros(1,51);
Coefficient_4_2=zeros(1,51);
Coefficient_4_3=zeros(1,51);
%以下求在不同扩展角范围条件下的信道容量，扩展角依次是40°，20°，10°，5°，1°
MyAngleSpread=40; %扩展角为40°，这样积分范围不一样，得到的相关矩阵也就不一样
ChannelCapacity;  %求扩展角为40°的条件下，信道容量随天线间距/波长的变化情况，天线间距/波长从0以0.1为步长到5
results(:,1)=Capacity;

MyAngleSpread=20;%扩展角为20°
ChannelCapacity; %求扩展角为20°的条件下，信道容量随天线间距/波长的变化情况，天线间距/波长从0以0.1为步长到5
results(:,2)=Capacity;

MyAngleSpread=10;%扩展角为10°
ChannelCapacity;%求扩展角为10°的条件下，信道容量随天线间距/波长的变化情况，天线间距/波长从0以0.1为步长到5
results(:,3)=Capacity;

MyAngleSpread=5;%扩展角为5°
ChannelCapacity;%求扩展角为5°的条件下，信道容量随天线间距/波长的变化情况，天线间距/波长从0以0.1为步长到5
results(:,4)=Capacity;

MyAngleSpread=1;%扩展角为1°
ChannelCapacity;%求扩展角为1°的条件下，信道容量随天线间距/波长的变化情况，天线间距/波长从0以0.1为步长到5
results(:,5)=Capacity;
%做出的图形表示的是不同扩展角条件下信道容量随天线间距的变化情况。
dd=0:0.1:5;
plot (dd,results, 'DisplayName','results', 'YDataSource', 'results'); figure(gcf)
grid