clc;
clear;
n=35;
k=27;
sourcecode=[1?0?0?1?0?1?0?1?0?1?1?0?1?0?1?0?0?1?1?0?1?0?1?0?1?0?1];
g=[1?0?1?1?0?1?0?1?1];
disp（sourcecode）
a=zeros（1n+1）;%x.^n-1
a（1）=1;
a（n+1）=-1;
G=zeros（kn）;
for?i=1:k
????G（ii:n-k+i）=fliplr（g）;
end
f_code=sourcecode*G;
f_code=rem（f_code2）;

[hv]=deconv（ag）;
v=mod（v2）;
v_len=length（v）;
t=zeros（1v_len）;
if?isequal（tv）
????h=mod（h2）;
else
????disp（‘generate?polynomial?is?wrong‘）;
end

H=zeros（n-kn）;
for?i=1:n-k
????H（ii:k+i）=h;
end
lens=3*n-3;
s1=eye（n）;
s2=zeros（n-1n）;
for?i=1:n-1
????s2（ii:i+1）=[1?1];
end
s3=zeros（n-2n）;
for?i=1:n-2
????s3（ii:i+2）=[111];
end
b=[s1;s2;s3];
err=mod（b*H.‘2）;

%?%?error_add?module%?%
rat=0.1;
x=f_code;
for?i=1:n
????if（rand<