BCH码编译码.zip_BCH 误码率_BCH码性能比较_BCH码比较_BCH码误码率_bch码性能仿真

function decodmsg = bolicanpu(t,k,shoumazi) %% 构造伽罗华域 m=6; n=2^m-1; p=(de2bi(primpoly(m),m+1)); a_gf=gf(0:(2^(length(p)-1)-1),length(p)-1,bi2de(p)); p_test=seqreverse(p_tongyu(p)); p_test_de(:,2)=bi2de(seqreverse(p_test)); p_test_de(:,1)=-1:n-1; %GF(2^6) %% 译码 for q=1:size(shoumazi,1) %% 计算伴随式 s=zeros(2*t,m);s_order=zeros(2*t,1); for j=1:2*t for i=1:n if shoumazi(q,i)~=0 s(j,:)=gfadd(s(j,:),p_test(mod(j*(n-i),n)+2,:),2); end end s_order(j,1)=find(bi2de(seqreverse(s(j,:)))==p_test_de(:,2))-2; %减2是为了保证s_order中存储的是x的阶数 end %% 求解delta(x) s=bi2de(seqreverse(s))'; %伴随式转成十进制 f=a_gf(1:length(s)); for i=1:length(s) f(i)=a_gf(s(i)+1); end s=f;e=length(s)/2; S_buf=zeros(e,e+1); S_buf=gf(S_buf,m); for i=1:e S_buf(i,:)=s(e+i:-1:i); end S_buf=S_buf'; S=[S_buf(2:e+1,1:e);S_buf(1,:)]'; [S_type,e] = zg_jie(S,e); [~,b]=size(S_type); delta=S_type(:,b)'; delta_x_final=zeros(1,length(delta)); buf=delta.x; for i=1:length(delta) a=find(buf(1,i)==p_test_de(:,2)); delta_x_final(1,i)=a-2; end %% 钱氏搜索法计算delta(x)=0的根以及错误位置 delta_x_root=zeros(length(delta_x_final),1); hang=zeros(length(delta_x_final),1); dxs=zeros(length(delta_x_final),m); cishu=zeros(1,length(delta_x_final)); j=1; for i=1:n zhi=zeros(1,m); for z=1:length(delta_x_final) cishu(1,z)=mod(delta_x_final(1,z)+z*i,n); end for z=1:length(delta_x_final) hang(z,:)=find(cishu(1,z)==p_test_de(:,1)); end for z=1:length(delta_x_final) dxs(z,:)=seqreverse(de2bi(p_test_de(hang(z,:),2),m)); end for z=1:length(delta_x_final) zhi=mod(gfadd(zhi,dxs(z,:),2),2); end zhi=bi2de(seqreverse(zhi)); if (zhi==1) delta_x_root(j,:)=i; j=j+1; end end error_location=n-delta_x_root; %确定最后错误位置数 %% 译码结果 e_x=zeros(1,n); for v=1:length(delta_x_root) if((delta_x_root(v))~=0) e_x=mod(seqreverse((de2bi(2^(error_location(v)),n)))+e_x,2); end end shoumazi(q,:)=mod(shoumazi(q,:)+e_x,2); end %% 提取信息位 for q=1:size(shoumazi,1) decodmsg(1,(k*(q-1)+1):k*q)=shoumazi(q,1:k); end end