cdf_vs_snr_cas.m

改程序实现了mimo无线通信的容量的CDF、PDF等分布曲线

clear
clc
hold on
M=4;
N=1;L=16;
xx=[401 401 401 401];
yy=[401 401 401 401];
x=1;y=1;
Pt=1;
snr_dB=21;
snr=10^(snr_dB/10);
I=eye(M,M);
Pout=10/100;
dota=1;
sigma=8;
    for i=1:N
        d(i)=sqrt((x-xx(i)).^2+(y-yy(i)).^2);
    end
    for count=1:10000
        H=[];
        for i=1:N
            u(i)=10*log10(Pt/(N*L))-4*10*log10(d(i)/1000);
            s_dB(i)=u(i)+sigma*randn(1,1);
            s(i)=10^(s_dB(i)/10);
            Hi=sqrt(s(i)/2)*randn(M,L)+sqrt(-1)*sqrt(s(i)/2)*randn(M,L);
            H=[H Hi];
        end
        C(count)=log2(det(I+(snr/(N*L))*H*H'));
    end
    C=abs(C);

M=100;
[x,pdfx,cdfx,ccdfx]=pdf_cdf(C, M);
plot(x,cdfx,'--')
%xlabel('Capacity: bit/s/Hz')
%ylabel('CDF')
%grid