ldpc_demo.m

实用的LDPC编解码MATLAB程序

%clear all;
%clc;
tic

rows=256;
cols=512;


s=round(rand(1, cols-rows));

%产生H矩阵
H=genH(rows,cols);

%使用H矩阵进行LDPC编码
[u,P,rearranged_cols]=ldpc_encode(s,H);

SNR=2;
amp=1;
tx_waveform=bpsk(u,amp);
rx_waveform=awgn(tx_waveform,SNR);
scale(1:length(u))=1;  %No fading.

%LDPC译码
[uhat vhat]=ldpc_decode(rx_waveform,SNR,amp,scale,H,rearranged_cols);

errmax=find(s~=uhat);
nerr=length(errmax);
ber=nerr/cols
toc