循環(huán)碼編譯碼matlab程序

1、精品文檔循環(huán)碼編譯碼matlab 程序循環(huán)碼編碼程序function C = cyclic_encoder( Si )%8循環(huán)編碼的輸出編碼結(jié)果%寸 xA8+1 進行模 2 因式分解得到：xA8+1=(xA3+xA2+x+1)*(xA5+xA4+x+1)y=size(Si,2);%y 表示 Si 的列數(shù)，即輸入碼元的個數(shù)M=ceil(y/5);%將信息碼元分成 M幀，一幀5個信息碼元n=8;%循環(huán)編碼的一幀碼長k=5；%信息位的個數(shù)r=n-k;% 監(jiān)督位的個數(shù)gx=1,1,1,1;%(8,5) 循環(huán)碼的生成多項式g(x)=xA3+xA2+x+1Ai=zeros(1,8*M);%Ai 用來存放所

2、輸入的碼元經(jīng)過循環(huán)編碼后的碼字Axi=zeros(1,8);%Axi 用來表示循環(huán)編碼后的一幀的編碼輸出碼字mi=zeros(1,5);%mi 用來存放每一幀的信息碼元for i=1:Mfor j=1:5mi(j)=Si(j+(i-1)*5);endAxi(4:8)=mi(1:5);Axi=circshift(Axi',-r)'%實現(xiàn) (xA(n-k)*m(x), 其中 m(x) 的系數(shù)由mi 決定qx,rx=deconv(Axi,gx);% 實現(xiàn) (xA(n-k)*m(x)/g(x), 得到商 q(x) 和余數(shù) r(x)Axi=Axi+rx;% 實現(xiàn) Axi(x)=Axi(x

3、)+r(x), 得到的 Axi 就是循環(huán)編碼的編碼輸出碼字Ai(8*i-4:8*i)=Axi(1:5);Ai(8*i-7:8*i-5)=Axi(6:8);end%for 循環(huán)是為了實現(xiàn)模2 相加，使循環(huán)編碼的輸出碼字Ai 中只有 0,1for i=1:8*Mif rem(abs(Ai(i),2)=0Ai(i)=0;elseAi(i)=1;endendC=Ai;%循環(huán)編碼的輸出碼字C=Aiend循環(huán)碼譯碼程序function So = cyclic_decoder( R )喻入R為經(jīng)AWGN1道傳輸后的二進制信息，So為循環(huán)譯碼器的譯碼結(jié)果R=1*(R>0.5);% 對接收到的信號進行抽樣

4、判決y=size(R,2);%y表示R的列數(shù)，即輸入碼元的個數(shù)M=ceil(y/8);% 將接收到的碼元 R分成M幀，一幀8個碼元So=zeros(1,5*M);% 用來存放糾檢錯之后的譯碼結(jié)果Axo=zeros(1,8);% 用來存放糾檢錯之后每一幀的譯碼結(jié)果%對接收信號進行糾檢錯譯碼n=8;%循環(huán)編碼的碼長k=5；%信息位的個數(shù)r=n-k;% 監(jiān)督位的個數(shù)s=zeros(1,3);% 用來存放校正子sh=1,1,0,0,1,1;% 監(jiān)督多項式 h(x)=xA5+xA4+x+1hn=;%hn(x) 為監(jiān)督多項式h(x) 的逆多項式，hn 用來放 hn(x) 的各項系數(shù)for i=2:5hn(

5、i)=h(k+2-i);endhn(1)=h(1);hn(6)=h(6);%計算監(jiān)督矩陣HH=zeros(r,n);% 監(jiān)督矩陣H 為 r*n 階矩陣H0=zeros(1,n);%用來存放hn(x) 的系數(shù)H1=zeros(1,n);%用來存放x*hn(x) 的系數(shù)H2=zeros(1,n);%用來存放xA2*hn(x)的系數(shù)H0(3:8)=hn(1:6);H1(2:7)=hn(1:6);H2(1:6)=hn(1:6);H(1,:)=H2(1:8);H(2,:)=H1(1:8);H(3,:)=H0(1:8);flag=0;% 出錯的標(biāo)志，為1 表示檢測出錯誤for i=1:Mfor j=1:8Axo(j)=R(j+(i-1)*8);ends=Axo*H'% 計算校正子sfor k=1:3if rem(abs(s(k),2)=0 s(k)=0;elses(k)=1;endendif s=0 0 0if flag=0flag=0;endelseflag=1;endfor k=1:8if rem(abs(Axo(k),2)=0Axo(k)=0;elseAxo(k)=1;endendSo(1+(i-1)*5:5+(i-1)*5)=A