代码拉取完成,页面将自动刷新
%% Clean up
clear all;
clc;
%% Init
if ~isdeployed
addpath('./codes');
end
% Constellation size
M = 4;
% LDPC config
blkSize = 256;
codeRate = '5/6';
% Get LDPC struct
LDPC = ldpcGet(blkSize, codeRate);
% Simulation parameters
ebno = 8;
numIter = 1e2;
numErr = 0;
% Convert E_b/N_0 to some SNR
snr = ebno + 10*log10(log2(M)) + 10*log10(str2num(codeRate));
%% Simulate
for i = 1:numIter
% Generate random data
data = randi([0 1], 1, LDPC.numInfBits);
% Encode
dataEnc = ldpcEncode(data, LDPC);
% QAM mapping
dataMod = qammod(dataEnc(:), M, 'InputType', 'bit', 'UnitAveragePower', true);
% AWGN
dataRx = awgn(dataMod, snr);
% LLR demapping
dataLlr = qamdemod(dataRx, M, 'OutputType', 'llr', 'UnitAveragePower', true);
% Decode
dataHat = ldpcDecode(dataLlr', LDPC);
% Count number of bit errors
numErr = numErr + sum(abs(dataHat - data));
end
%% BER
ber = numErr / (numIter * LDPC.numInfBits)
此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。