离散信道的线性分组编码和译码实现

clear
%编码
G=input('请输入生成矩阵G\n G=');
[k,n]=size(G);
r=n-k;
m=input('请输入需传送信息m\n m=');
l=length(m);
if (mod(l, k))
    disp('输入的信息有误');
else
    g=l/k;
%将输入序列转化成矩阵m
t1=[];
for i=1:g
    t1(i,:)=m(k*(i-1)+1:i*k) ;
end
m=t1;
%求校验矩阵H
c=mod(m*G, 2);
A=G(:,k+1:n);
H=[A',eye(r)];
disp('校验矩阵');H
disp('译码矩阵' );c
end
disp('敲回车键继续'); pause
%解码
y= input('输入接收序列y\n y=');
t2=[]; 
for i=1:g
    t2(i,:)=y(1,n*(i-1)+1:i*n);
end
y=t2
s=mod(y*H',2) ;
e=s*pinv(H');
for i=1:g
    for j=1:n
        if(e(i,j)>0.5-eps)
            e(i,j)=1;
        else
            e(i,j)=0;
        end
    end
end
c1=mod(y+e,2);
disp('差错图样'); e
disp('译码序列'); c1