lab7_1

clear all; close all;
clc
Nr=[1 2 4 8]; % число элементов приемной АР
Nt=[1 2 4 8]; % число элементов передающей АР
for j=1:length(Nr) % цикл по числу элементов АР
 I=eye(Nr(j)); % инициализация единичной матрицы
 g=100; % ОСШ
 count = 1e5; % число случайных реализаций
 for n=1:count
 % инициализация канальной матрицы H
 H=sqrt(1/2)*randn(Nr(j),Nt(j))+j*sqrt(1/2)*randn(Nr(j),Nt(j));
 % Оценка ПС по формуле (7.32)
 C(n)=real(log2(det(I+(g/Nt(j))*(H*H'))));
 end
 % оценка ПВ ПС
 Cbin_size=100;
 [Cfreq,Cbin]=hist(C,Cbin_size);
 Cstep=(max(Cbin)-min(Cbin))./(Cbin_size-1);
 Cpdf=Cfreq./(sum(Cfreq)*Cstep);
 % оценка ИФР ПС
 for i=1:length(Cbin)
 Ccdf_range=length(find(C<Cbin(i)));
 ICcdf(i)= 1-Ccdf_range/length(C);
 end
 % построение графиков
 subplot(length(Nr),2,2*j-1); bar(Cbin,Cpdf);
 ylabel('ПВ'); grid on;
 title(['N_{tx} = ',num2str(Nt(j)),', N_{rx} = ',num2str(Nr(j))]);
 legend(['M[C] = ', num2str(round(mean(C))),' (бит/с/Гц)']);
 if j==length(Nr)
 xlabel(['Случайная реализация ПС (бит/с/Гц) при ОСШ = ',...
 num2str(10*log10(g)), ' [дБ]']);
 end
 subplot(length(Nr),2,2*j); plot(Cbin,ICcdf,'r-', 'linewidth',2);
 ylabel('ИФР'); grid on;
 title(['N_{tx} = ',num2str(Nt(j)),', N_{rx} = ',num2str(Nr(j))]);
 legend(['C_{90%} = ',...
 num2str(round(mean(Cbin(find(abs(ICcdf-0.9)<Cstep/2))))),'(бит/с/Гц)']);
 if j==length(Nr)
 xlabel(['Случайная реализация ПС (бит/с/Гц) при ОСШ = ',...
 num2str(10*log10(g)),' [дБ]']);
 end
end