通信原理实验报告(MatLab下循环码).docVIP

通信原理实验报告 ----------基于MatLab下的循环码仿真 关键词：生成多项式 信息码组 生成码组 分析： 循环码（n，k）是一种高效而重要的线性分组码，通过生成多项式g（x）可以得到一个完整的循环码组，在此我们使用matlab软件对其进行程序编写，程序由主程序，编码子函数，解码子函数组成，其中测试用的码为（7,4）循环码组，生成多项式为g（x）=x^3+x+1, 实验环境和运行平台 WindowsXP ＋matlab命令窗口 代码 主函数 （xunhuan.m） g=[1 0 1 1]; u1=[0 0 0 1]; u2=[0 0 1 0]; u=[1 0 0 1]; n=7;k=4; code=circue(u,g,n); code1=circue(u1,g,n); code2=circue(u2,g,n); y1=decode(code,g,k); 编码函数 调用div.m编码 （circue.m） function code=circue(cin,g,n) [tmp k]=size(cin); jiaoyan=div(g,cin,n,k); jiaoyan=jiaoyan(k+1:end); code=[cin jiaoyan]; 子函数 （div.m）生成效验码 function jiaoyan=div(g,u,n,k) [x len_g]=size(g); %g=[g zeros(1,n-k)]; u=[u zeros(1,n-k)]; [y len_u]=size(u); if(len_glen_u) u=[zeros(1,len_g-len_u) u]; else g=[zeros(1,len_u-len_g) g]; end uu=find(u); gg=find(g); um=uu(1); gm=gg(1); beichu=u; while(um=gm) tmp=gm-um; %gtmp=[g zeros(tmp)]; gtmp=zuoyi(g,tmp); check=xor(beichu,gtmp); uutmp=find(check); umtmp=uutmp(1); beichu=check; um=umtmp; end jiaoyan=beichu; 解码函数 调用shang.m解码 （decode.m） function y=decode(in,g,k) n=length(in); h=shang(n,g,k); hh=zeros(n-k,n); for i=1:(n-k) htmp=zuoyi(h,i-1); hh(n-k-i+1,:)=htmp; end s=in*hh; s=mod(s,2); e=s/hh; ss=sum(e); aa=find(ss==max(ss)); aa=aa(1); e=e(aa,:); y=xor(in,e); 子函数1 产生错误图像 （shang.m） function h=shang(n,g,k); x=zeros(1,n+1); x(1)=1;x(n+1)=1; [t xm]=size(g); g=[zeros(1,n+1-xm) g]; beichu=x; h1=zeros(1,n+1); h=zeros(1,n); a=find(beichu); am=a(1); b=find(g); bm=b(1); hci=bm-am; while(am=bm) tmp=bm-am; h1(n+1-tmp)=1; gtmp=zuoyi(g,tmp); check=xor(beichu,gtmp); beichu=check; a=find(beichu); if(~isempty(a)) am=a(1); else break; end end h1=h1(2:end); for i=n-k:n h(n-i+n-k)=h1(i); end 子函数2 左移函数 （zuoyi.m） function aa=zuoyi(a,b) aa=[a(b+1:end) zeros(1,b)];