channelcorrelationmatrix.m

Capacity_and_Channel for mimo channel

% MyChannelCorrelationMatrix
clear;
global DOA sigma

AngleSpread=30/180*pi;
sigma=AngleSpread/sqrt(12);
sigma=AngleSpread;
DOA=20/180*pi;

CorrelationMatrix=zeros(4,4);

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

temp1=DOA-AngleSpread/2;
temp2=DOA+AngleSpread/2;
temp1=-4*pi;
temp2=4*pi;
Coefficient_1_2=quadv(@CorrCoe_laplace_1_2,temp1,temp2);
Coefficient_1_3=quadv(@CorrCoe_laplace_1_3,temp1,temp2);
Coefficient_1_4=quadv(@CorrCoe_laplace_1_4,temp1,temp2);
Coefficient_4_1=quadv(@CorrCoe_laplace_4_1,temp1,temp2);
Coefficient_4_2=quadv(@CorrCoe_laplace_4_2,temp1,temp2);
Coefficient_4_3=quadv(@CorrCoe_laplace_4_3,temp1,temp2);

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

CorrelationMatrix=abs(CorrelationMatrix);
lhw=CorrelationMatrix