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