%clear all;
%clc;
tic

rows=128;
cols=256;


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

%產生H矩陣
H=genH(rows,cols);

%使用H矩陣進行LDPC編碼
[u,P,rearranged_cols]=ldpc_encode(s,H);

SNR=10;
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)

toc