线性分组码的信道编码和译码.doc

《线性分组码的信道编码和译码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、clear;clc;%编码G=input('请输入生成矩阵G,例如:G=[10111;01101]G=');[k,n]=size(G);r=n-k;m=input('请输入需传送信息m,如m=[00011011]m=');l=length(m);if(mod(l,k))disp('输入的信息有误');elsege=l/k;%将输入序列转化成矩阵mtemp1=[];fori=1:getemp1(i,:)=m(k*(i-1)+1:i*k);endm=temp1;%求校验矩阵Hc=mod(m*G,2);A=G(:,

2、k+1:n);H=[A',eye(r)];disp('校验矩阵');Hdisp('译码矩阵');cenddisp('敲回车键继续');pause%解码y=input('输入接收序列y,如:y=[00000011011011110010]y=');temp2=[];fori=1:getemp2(i,:)=y(1,n*(i-1)+1:i*n);endy=temp2s=mod(y*H',2);e=s*pinv(H');fori=1:geforj=1:nif(e(i,j)>0.5-eps)e(i,j)=1;elsee(i

3、,j)=0;endendendcc=mod(y+e,2);%cc=xor(y,e)sc=cc(:,1:2);disp('差错图样');edisp('估计值');ccdisp('译码序列');sc