ldpc_demo1.m

实现ldpc的编译码功能

clear all;
clc;
rows=1024;
cols=2048;


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

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

%使用H矩阵进行LDPC编码
[u,P,rearranged_cols]=ldpc_encode(s,H);
SNR=5;
amp=1;
tx_waveform=bpsk(u,amp);
rx_waveform=awgn(tx_waveform,SNR);

scale(1:length(u))=1;  %No fading.

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

errmax=find(s~=uhat);
nerr=length(errmax);
BER=nerr/(cols-rows);

%LDPC性能图
plot(SNR,BER,'-.r*')
set(gca,'Yscale','log');
hold on

grid on;
ylabel('BER');
xlabel('SNR');
title('LDPC在AWGN下的性能曲线');