实验6 BCH循环码的编码与译码的matlab实现.doc

实验6 BCH循环码的编码与译码 一、实验内容 用VC或Matlab软件编写循环BCH码的编码与译码程序。利用程序对教科书的例题做一个测试。 二、实验环境 计算机 Windows 2000 或以上 Microsoft Visual C++ 6.0 或以上 Matlab 6.0或以上 三、实验目的 通过BCH循环码的编码与译码程序的编写，彻底了解并掌握循环BCH的编码与译码原理 通过循环BCH码的编码与译码程序的编写，提高编程能力。 四、实验要求 提前预习实验，认真阅读实验原理以及相应的参考书。 对不同信道的进行误码率分析。特别是对称信道，画出误码性能图。即信道误码率与循环汉明码之间的关系。 认真填写实验报告。 五、实验原理 循环BCH的编码与译码原理(略) 循环BCH的程序实现。 六、实验步骤 bch_en_decode.m文件 function bch_en_decode() code=bch155 code=code+randerr(5,15,1:3); code=rem(code,2); code=gf(code) %随机产生1-3位错误 decode=debch155(code) end function decode=debch155(code) code=gf(code); M=4; code = gf(code.x,M); [m , n]=size(code);decode=[]; code1=[]; for i=1:m ;code1=code(i,:); M=code1.m;T2=6;N=15; S = code1* ((gf(2,M,code1.prim_poly)).^([N-1:-1:0]*([1:T2]))); LambdaX = gf([1 zeros(1,T2)],M,code1.prim_poly); Tx = [0 1 zeros(1,T2-1)]; L=0; for k = 1:T2; LambdaXTemp = LambdaX; Delta = S(k) - LambdaXTemp(1+[1:L])*(S(k-[1:L])); if Delta.x; LambdaX = LambdaXTemp - Delta*Tx; if 2*L k; L = k-L; Tx = LambdaXTemp/Delta; end; end; Tx = [0 Tx(1:T2)]; end; LambdaXValue = LambdaX.x; LambdaX = gf(LambdaXValue(1:max(find(LambdaXValue))), M, code1.prim_poly); errLoc_int = roots(LambdaX);errLoc = log(errLoc_int); for i = 1:length(errLoc); errorMag = 1; code1(N-errLoc(i)) = code1(N-errLoc(i)) - errorMag; end; decode=[decode;code1]; end; ccode = gf(decode.x); decoded = ccode(:,1:5); end function [yout]=bch155(x) %定义函数 k=5; %信息码位，BCH（15，5） if nargin1 x2=randint(5,k); n=5; msg=x2 %判断输入信息 ,若未输入,系统自动产生5组信息码，并显示出信息位 elseif rem(length(x),k)==0;n=length(x)/k;x2=[]; %判断msg是否为K的整数倍，并把输入码员分组 for i=0:n-1