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

waveform=bpsk(c);                     %BPSK調(diào)制           


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

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

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