%% 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;
agd.zip_Windows编程_bandfy7_unknown9sp_信道容量
版权申诉
29 浏览量
2022-09-24
10:24:33
上传
评论
收藏 9KB ZIP 举报
agd.zip （27个子文件） 
ChannelCapacity_batch.m 769B ChannelCorrelationMatrix.m 1KB ChannelCorrelationCoefficient_Gauss.m 881B CorrCoe_Gauss.m 202B 9Capacity_ccdf.m 1KB ChannelCorrelationCoefficient.m 673B CorrCoe_laplace.m 129B CorrCoe_laplace_1_4.m 155B CorrCoe_3d.m 78B CorrCoe_laplace_4_1.m 154B CorrCoe_1_4.m 120B CorrCoe_laplace_4_3.m 154B CorrCoe_4_1.m 119B CorrCoe_laplace_1_2.m 155B ChannelCapacity.m 1KB CorrCoe_2d.m 78B CorrCoe_1d.m 78B CorrCoe_1_2.m 118B CorrCoe_4_3.m 119B CorrCoe_laplace_4_2.m 152B CorrCoe.m 116B CorrCoe_4_2.m 119B mimo_channel_batch.m 864B FCapacity_ccdf_batch.m 269B rccdf_test.m 281B CorrCoe_laplace_1_3.m 153B CorrCoe_1_3.m 120B
评论0