Может ли кто-нибудь объяснить мне Matlab код подробно??

[johnnycage65]

Это Matlab кодекса, которая дает UA Символ ошибки вероятность кривой на 64-QAM модуляция.
Теперь я хочу кого-то объяснить мне этот код подробно с самого начала.Как это работает??Буду очень признателен за это.

Вот код:
-------------------------------------------------- ------------------------
N = 7 * 10 ^ 5;%% число символов
M = 64;%% количества точек созвездия

K = SQRT (1 / ((2 / 3) * (М-1)));%% нормализации фактор

M = [1: SQRT (M) / 2];%% алфавиты
alphaMqam = [- (2 * M-1) 2 * M-1];

Es_N0_dB = [0:30];%% Es/N0 несколько значений

ipHat = нуля (1, N);%% инициализации

за II = 1: длина (Es_N0_dB)

IP = randsrc (1, N, alphaMqam) J * randsrc (1, N, alphaMqam);
S = K * IP;%% нормализации энергии на 1
N = 1/sqrt (2) * [randn (1, N) J * randn (1, N)];% guassian белый шум, 0dB разница

Y = S 10 ^ (-Es_N0_dB (II) / 20) * N;%% аддитивного гауссовского белого шума

%% демодуляции
y_re = Real (Y) / K;% реальной части
y_im IMAG = (Y) / K;% мнимой части

%% Округления до ближайшего алфавиту
%% 0 2 -> 1
% 2 4 -> 3
% 4 6 -> 5 и
т.д.ipHat_re = 2 * этаж (y_re / 2) 1;
ipHat_re (Find (ipHat_re> Max (alphaMqam))) = Max (alphaMqam);
ipHat_re (Find (ipHat_re <мин (alphaMqam))) = мин (alphaMqam);

%% Округления до ближайшего алфавиту
%% 0 2 -> 1
%% 2 4 -> 3
%% 4 6 -> 5 и
т.д.ipHat_im = 2 * этаж (y_im / 2) 1;
ipHat_im (Find (ipHat_im> Max (alphaMqam))) = Max (alphaMqam);
ipHat_im (Find (ipHat_im <мин (alphaMqam))) = мин (alphaMqam);

ipHat = ipHat_re J * ipHat_im;
nErr (II) = размер (найти ([IP-ipHat]), 2);%% подсчете количества ошибок

конец

simSer = nErr / N;
theorySer = 2 * (1-1/sqrt (М)) * erfc (K * SQRT ((10. ^ (Es_N0_dB/10)))) ...
- (1-2/sqrt (M) 1 / M) * (erfc (K * SQRT ((10. ^ (Es_N0_dB/10 ))))).^ 2;
закрыть все
цифра
semilogy (Es_N0_dB, theorySer, 'BS-', 'LineWidth', 2);
держитесь
semilogy (Es_N0_dB, simSer, 'м *-',' Linewidth', 1);
Ось ([0 30 10 ^ -5 1])
по сетке
Легенда ( 'теория', 'моделирование');
xlabel ( 'Es / No, дБ')
ylabel ( 'Символ ошибки Оценить')
заголовок ( 'Символ ошибки вероятность кривой на 64-QAM модуляцией')
-------------------------------------------------- ------------------------------