channelcapacity.m

Capacity_and_Channel for mimo channel

%% Channel Capacity of MIMO 

temp=MyAngleSpread/2/180*pi;
Coefficient_1_2=quadv(@CorrCoe_1_2,-temp,temp)/2/temp;
Coefficient_1_3=quadv(@CorrCoe_1_3,-temp,temp)/2/temp;
Coefficient_1_4=quadv(@CorrCoe_1_4,-temp,temp)/2/temp;
Coefficient_4_1=quadv(@CorrCoe_4_1,-temp,temp)/2/temp;
Coefficient_4_2=quadv(@CorrCoe_4_2,-temp,temp)/2/temp;
Coefficient_4_3=quadv(@CorrCoe_4_3,-temp,temp)/2/temp;

Re_CorrelationMatrix(1,2,:)=Coefficient_1_2;
Re_CorrelationMatrix(1,3,:)=Coefficient_1_3;
Re_CorrelationMatrix(1,4,:)=Coefficient_1_4;
Re_CorrelationMatrix(2,1,:)=Coefficient_4_3;
Re_CorrelationMatrix(2,3,:)=Coefficient_1_2;
Re_CorrelationMatrix(2,4,:)=Coefficient_1_3;
Re_CorrelationMatrix(3,1,:)=Coefficient_4_2;
Re_CorrelationMatrix(3,2,:)=Coefficient_4_3;
Re_CorrelationMatrix(3,4,:)=Coefficient_1_2;
Re_CorrelationMatrix(4,1,:)=Coefficient_4_1;
Re_CorrelationMatrix(4,2,:)=Coefficient_4_2;
Re_CorrelationMatrix(4,3,:)=Coefficient_4_3;

Re_CorrelationMatrix=real(Re_CorrelationMatrix);

Capacity=zeros(1,51);

for mm=1:51
    for nn=1:1500
        WhiteMatrix=sqrt(1/2)*randn(4)+j*sqrt(1/2)*randn(4);
        Capacity(mm)=Capacity(mm)+log2(det(eye(4)+25*(Re_CorrelationMatrix(:,:,mm))'...
            *WhiteMatrix*eye(4)*WhiteMatrix'));
    end
    Capacity(mm)=Capacity(mm)/1500;
end

Capacity=real(Capacity)';