Эквализация минимальной среднеквадратической ошибки (MMSE)
[ применяет эквализацию MMSE к извлеченным элементам ресурса физического канала eqSym,csi]
= nrEqualizeMMSE(rxSym,hest,nVar)rxSym и возвращает компенсируемые символы в eqSym. Процесс эквализации использует предполагаемую информацию о канале hest и оценка полученного шумового отклонения nVar. Функция также возвращает мягкую информацию о состоянии канала csi.
Выполните эквализацию MMSE на извлеченных элементах ресурса физического канала телевизионного вещания (PBCH).
Создайте символы и индексы для передачи PBCH.
ncellid = 146; v = 0; E = 864; cw = randi([0 1],E,1); pbchTxSym = nrPBCH(cw,ncellid,v); pbchInd = nrPBCHIndices(ncellid);
Сгенерируйте пустой массив ресурса для одной антенны передачи. Заполните массив с символами PBCH при помощи сгенерированных индексов PBCH.
P = 1; txGrid = zeros([240 4 P]); txGrid(pbchInd) = pbchTxSym;
Выполните модуляцию OFDM.
txWaveform = ofdmmod(txGrid,256,[22 18 18 18],[1:8 249:256].');
Создайте матрицу канала и примените канал к переданной форме волны.
R = 4; H = dftmtx(max([P R])); H = H(1:P,1:R); H = H / norm(H); rxWaveform = txWaveform * H;
Создайте оценку канала.
hEstGrid = repmat(permute(H.',[3 4 1 2]),[240 4]); nEst = 0.1;
Выполните демодуляцию OFDM.
rxGrid = ofdmdemod(rxWaveform,256,[22 18 18 18],0,[1:8 249:256].');
Чтобы подготовиться к декодированию PBCH, используйте nrExtractResources чтобы извлечь символы из полученного и канала оценивают сетки. Постройте полученное созвездие PBCH.
[pbchRxSym,pbchHestSym] = nrExtractResources(pbchInd,rxGrid,hEstGrid); figure; plot(pbchRxSym,'o:'); title('Received PBCH Constellation');
Декодируйте PBCH с извлеченными элементами ресурса. Постройте компенсируемое созвездие PBCH.
[pbchEqSym,csi] = nrEqualizeMMSE(pbchRxSym,pbchHestSym,nEst); pbchBits = nrPBCHDecode(pbchEqSym,ncellid,v); figure; plot(pbchEqSym,'o:'); title('Equalized PBCH Constellation');
