fire码matlab实现

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<