先谢谢你,本人大一,初次接触matlab,做的mimo遍历容量仿真,老师给的程序如下,能否尽量每行注解或者关键处注解,目前我注解了部分,如有不对,希望你帮忙修改一下,谢谢了。
有些缩写我注明一下:PDF:probability density function,概率密度函数
nT、nR:应该是收发天线数。
CDF: cumulative distribution function,累计分布函数
% mimo_bianli_CDF.m
clear all, close all
SNR_dB=10; %定义信噪比10dB
SNR_linear=10.^(SNR_dB/10.);%数字化定义信噪比
N_iter=50000;
sq2=sqrt(0.5);
grps = ['b:'; 'b-'];
for Icase=1:2
if Icase==1, nT=2; nR=2; % 2x2
else nT=4; nR=4; % 4x4
end
n=min(nT,nR); I = eye(n);
for iter=1:N_iter
H = sq2*(randn(nR,nT)+j*randn(nR,nT));
C(iter) = log2(real(det(I+SNR_linear/nT*H'*H)));
end
[PDF,Rate] = hist(C,50);
PDF = PDF/N_iter;
for i=1:50
CDF(Icase,i) = sum(PDF([1:i]));
end
plot(Rate,CDF(Icase,:),grps(Icase,:)); hold on
end
xlabel('Rate[bps/Hz]')
ylabel('CDF')
axis([1 18 0 1]); grid on; set(gca,'fontsize',10);
legend('{\it N_T}={\it N_R}=2','{\it N_T}={\it N_R}=4');