實(shí)驗(yàn)6-BCH循環(huán)碼的編碼與譯碼的matlab實(shí)現(xiàn)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上實(shí)驗(yàn)6 BCH循環(huán)碼的編碼與譯碼一、實(shí)驗(yàn)內(nèi)容用VC或Matlab軟件編寫循環(huán)BCH碼的編碼與譯碼程序。利用程序?qū)炭茣睦}做一個(gè)測(cè)試。二、實(shí)驗(yàn)環(huán)境1. 計(jì)算機(jī)2. Windows 2000 或以上3. Microsoft Visual C+ 6.0 或以上4. Matlab 6.0或以上三、實(shí)驗(yàn)?zāi)康?. 通過BCH循環(huán)碼的編碼與譯碼程序的編寫，徹底了解并掌握循環(huán)BCH的編碼與譯碼原理2. 通過循環(huán)BCH碼的編碼與譯碼程序的編寫，提高編程能力。四、實(shí)驗(yàn)要求1. 提前預(yù)習(xí)實(shí)驗(yàn)，認(rèn)真閱讀實(shí)驗(yàn)原理以及相應(yīng)的參考書。2. 對(duì)不同信道的進(jìn)行誤碼率分析。特別是對(duì)稱信道，畫出誤碼性

2、能圖。即信道誤碼率與循環(huán)漢明碼之間的關(guān)系。3. 認(rèn)真填寫實(shí)驗(yàn)報(bào)告。五、實(shí)驗(yàn)原理1. 循環(huán)BCH的編碼與譯碼原理(略)2. 循環(huán)BCH的程序?qū)崿F(xiàn)。六、實(shí)驗(yàn)步驟bch_en_decode.m文件function bch_en_decode() code=bch155 code=code+randerr(5,15,1:3); code=rem(code,2); code=gf(code) %隨機(jī)產(chǎn)生1-3位錯(cuò)誤 decode=debch155(code)end function decode=debch155(code) code=gf(code); M=4; code = gf(code.x,M)

3、; 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)

4、' 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:lengt

5、h(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) %定義函數(shù)k=5; %信息碼位，BCH（15，5）if nargin<1 x2=randint(5,k); n=5; msg=x2 %判斷輸入信息 ,若未輸入,系統(tǒng)自動(dòng)產(chǎn)生5組信息碼，并顯示出信息位elseif rem(length(x

6、),k)=0;n=length(x)/k;x2=; %判斷msg是否為K的整數(shù)倍，并把輸入碼員分組 for i=0:n-1 x2=x2;x(i*k+1) x(i*k+2) x(i*k+3) x(i*k+4) x(i*k+5); end if rem(length(x),k)>0 %把輸入碼員補(bǔ)零并分組 x=x,zeros(size(1:k-rem(length(x),k); n=length(x)/k;x2=; for i=0:n-1 x2=x2;x(i*k+1) x(i*k+2) x(i*k+3) x(i*k+4) x(i*k+5); end endendik=eye(5) zeros

7、(5,10); %輸入信息碼擴(kuò)展x3=x2*ik;yout=;for i=1:n g=1 0 1 0 0 1 1 0 1 1 1; w,yo=deconv(x3(i,:),g); %產(chǎn)生余式 yo=abs(rem(yo,2); yout=yout;yo;endyout=yout+x3; %產(chǎn)生信息碼end運(yùn)行結(jié)果：msg = 1 1 0 1 0 0 1 1 0 1 0 1 0 0 0 1 1 1 0 1 0 0 1 0 0code = 1 1 0 1 0 1 1 0 0 1 0 0 0 1 1 0 1 1 0 1 1 1 0 0 0 0 1 0 1 0 0 1 0 0 0 1 1 1 1 0

8、1 0 1 1 0 1 1 1 0 1 0 1 1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 1 1 1 0 1 0 1 1 code = GF(2) array. Array elements = 1 1 0 1 0 1 0 1 0 1 0 0 0 1 1 0 1 1 1 1 1 1 0 0 0 0 1 0 1 0 0 0 0 0 0 1 1 1 1 0 0 0 1 1 0 1 1 0 0 1 0 1 1 1 0 1 0 0 0 1 0 0 1 1 0 0 1 1 1 1 0 1 0 1 1 decode = GF(24) array. Primitive polynomial = D4+D+1 (19 decimal) Array elements = 1 1 0 1 0 1 1 0 0 1 0 0 0 1