bch译码算法

Max_SNR=10;
N_trials=10000;
N=255;
k=239;
Eb=1;
ber_m1=0;
ber_m2=0;
ber_m3=0;
for trial=1:1:N_trials
    msg=round(rand(1,k));
    s1=1-msg.*2;
    n1=randn(1,k)+j.*randn(1,k);
    
    s0=bchenc(gf(msg),N,k);
    s=s0.x;
    s2=1-double(s)*2;
    n2=randn(1,N)+j.*randn(1,N);
    
    bit_s=0;%奇偶校验位
    for i=1:1:k
        bit_s=bitxor(bit_s,s(i));
    end
    s3=1-double([s,bit_s])*2;
    n3=randn(1,N+1)+j.*randn(1,N+1);
    
    ber_v1=[];
    ber_v2=[];
    ber_v3=[];
    for snr_dB=2:1:Max_SNR
        %无bch编码
        snr1=10.^(snr_dB./10);
        N01=Eb./snr1;
        sgma1=sqrt(N01./2);
        y=sqrt(Eb).*s1+sgma1.*n1;
        y1=sign