循环码编码

M=textread（‘M1.txt‘）   %从文件中读取信息元序列M
H=textread（‘H1.txt‘）   %从文件中读取一致校验矩阵H
m=length（H（:1））       %求解校验位m
k=length（M）      %求解信息位数k
n=k+m            %求解码长n

%生成发送码字C
G=[eye（k）H（:[1:k]）‘]   %计算生成矩阵G
C=M*G;           %产生码字，该码字为发送码字C
C=mod（C2）          %求模2运算，即变成二进制数

%生成接收码字R
E=randint（1n）  %产生只错一位的随机错误图样E
sum=0;
for i=1:n
      sum=sum+E（i）
end
R=xor（CE）;      %由发送码字C和错误图样E进行异或运算 
R=mod（R2）          %求模2运算，即变成二进制数

%计算伴随式S
S=R*H‘;    %计算伴随式S，H‘为H的转置
S