tic
m=252;n=504;Eb_N0=2;
BER=0;
R=(n-m)/n;
sigma_2=1/(2*10^(Eb_N0/10)*R);
for num=1:500
s=round(rand(1,n-m));                 %隨機產生長為(n-m)的信息序列
load G   
c=mod(s*G,2);                         %LDPC編碼

waveform=bpsk(c);                     %BPSK調制           


y=waveform+sqrt(sigma_2)*randn(1,n);  %加性高斯白噪聲信道

maxiter=100;                        %設置最大譯碼迭代次數maxiter
[v]=BP1(y,H,sigma_2,maxiter);      %LDPC譯碼(概率域(SPA1)和對數域上(SPA2)的和積算法,最小和算法(MSA))

v0=v(m+1:n);
err_max=find(s~=v0);               %尋找錯誤信息位
num_eer=length(err_max);            %求出錯誤信息位位數
BER0=num_eer/(n-m)                 %計算比特誤碼率BER
BER=BER+BER0;
end
BER=BER/500
toc
t=toc