channelcorrelationcoefficient_gauss.m

Capacity_and_Channel for mimo channel

%% MIMO_ChannelCorrelationCoefficien
clear;
clc;
global sigma
results=zeros(5,51);
x=0;
k=1/sqrt(12);
% AngleSpread=30/180*pi;
% temp=15/180*pi;
% sigma=k*AngleSpread;
% results(1,:)=quadv(@CorrCoe_Gauss,-temp,temp);

AngleSpread=20/180*pi;
temp=10/180*pi;
sigma=k*AngleSpread;
results(2,:)=quadv(@CorrCoe_Gauss,-temp,temp);

% AngleSpread=10/180*pi;
% temp=5/180*pi;
% sigma=k*AngleSpread;
% results(3,:)=quadv(@CorrCoe_Gauss,-temp,temp);
% 
% AngleSpread=5/180*pi;
% temp=2.5/180*pi;
% sigma=k*AngleSpread;
% results(4,:)=quadv(@CorrCoe_Gauss,-temp,temp);
% 
% AngleSpread=1/180*pi;
% temp=0.5/180*pi;
% sigma=k*AngleSpread;
% results(5,:)=quadv(@CorrCoe_Gauss,-temp,temp);

d2lamda=0:50;
newresults=real(results).';
plot (d2lamda,newresults, 'DisplayName','newresults', 'YDataSource', 'newresults'); figure(gcf)
axis([0 50 -0.2 1.2])
grid