二元序列的Berlekamp-Massey算法的matlab程序

clearclc
a = [10011011000111010100];%二元序列
l = zeros（1length（a）+1）;
d = zeros（1length（a）+1）;
n0 = find（a（:）~=01）;
fx = cell（1length（a）+1）;
for i = 1:n0
    fx{i}=1;
end
fx{n0+1} = zeros（1n0+1）;
fx{n0+1}（1）=1;
fx{n0+1}（n0+1）=1;
l（2:n0）=0;
d（n0）=a（n0）;
l（n0+1）=n0;
for i = n0+2:length（a）+1
    fx{i}=（mod（fx{i-1}2））;
    d（i-1） = mod（sum（a（i-1:-1:i-length（fx{i}））.*fx{i}）2）;
    if d（i-1） == 0