capacity_ccdf.m

关于mimo的程序

%% Channel Capacity of MIMO 

temp=MyAngleSpread/2/180*pi;
Coefficient_1_2=quad(@CorrCoe_1d,-temp,temp)/2/temp;
Coefficient_1_3=quad(@CorrCoe_2d,-temp,temp)/2/temp;
Coefficient_1_4=quad(@CorrCoe_3d,-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_1_2;
Re_CorrelationMatrix(2,3)=Coefficient_1_2;
Re_CorrelationMatrix(2,4)=Coefficient_1_3;
Re_CorrelationMatrix(3,1)=Coefficient_1_3;
Re_CorrelationMatrix(3,2)=Coefficient_1_2;
Re_CorrelationMatrix(3,4)=Coefficient_1_2;
Re_CorrelationMatrix(4,1)=Coefficient_1_4;
Re_CorrelationMatrix(4,2)=Coefficient_1_3;
Re_CorrelationMatrix(4,3)=Coefficient_1_2;

Re_CorrelationMatrix=real(Re_CorrelationMatrix);

Capacity=zeros(1,1);

X_Capacity=0:0.1:7;
ccdf=zeros(1,max(size(X_Capacity)));

for nn=1:1500
     WhiteMatrix=sqrt(1/2)*randn(4)+j*sqrt(1/2)*randn(4);
     ChannelMatrix=Re_CorrelationMatrix*WhiteMatrix;
     Capacity=log2(det(eye(4)+25*(Re_CorrelationMatrix)'...
            *WhiteMatrix*eye(4)*WhiteMatrix'));
     Capacity=real(Capacity)/4;
    for nnn=1:max(size(X_Capacity))
        if Capacity>=X_Capacity(nnn);
           ccdf(nnn)=ccdf(nnn)+1;
        else 
           ccdf(nnn);
        end
    end
end

ccdf=ccdf/1500;