cas_ergodic.m

改程序实现了mimo的遍历容量在仿真区域中的三维分布情况

clear
clc
M=4;
N=1;L=16;
xx=401;yy=401;
Pt=1;
snr_dB=20;
snr=10^(snr_dB/10);
I=eye(M,M);
dota=26;
sigma=8;
Pout=10/100;
for x=1:dota:801
    for y=1:dota:801
        d=sqrt((x-xx).^2+(y-yy).^2);
        for count=1:1000
            H=[];
            for i=1:N
                u(i)=10*log10(Pt/(N*L))-4*10*log10(d/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_ergodic(floor(x/dota)+1,floor(y/dota)+1)=(1/count)*sum(C);
    end
end
C_ergodic=abs(C_ergodic);
[m,n]=size(C_ergodic);
x=1:dota:m*dota;
y=1:dota:n*dota;
mesh(x-401,y-401,C_ergodic)
CAS_max=max(max(C_ergodic));
CAS_min=min(min(C_ergodic));
axis([-400 400 -400 400 0 CAS_max])
C_ergodic_mean=sum(sum(C_ergodic))/(m*n)

xlabel('x')
ylabel('y')
zlabel('ergodic capaciy of CAS')