Генерирование М-последовательности

clear;
% Массив с элементами М-последовательности
a = [];
% Веса характеристического полинома
c = [1 1 0 0 0 0 1];
% Откидываем отвод обратной связи из формулы
c = c(2: length(c));
% Степень полинома
n = length(c);
% Длина периода М-последовательности
M = 2 ^ n - 1;
% Заполнение нулями первой части массива
a = zeros(1,n-1);
a(n) = 1;
% Цикл расчета элементов одного периода
for j = n + 1 : M + n
   temp = 0;
   for k = 1 : n
       temp = xor(temp, c(k) * a(j-k));
   end;
   a(j) = temp;
end;
% Отсечение первых n элементов
a = a(n+1:length(a));
% Сохранение результатов в файл pnseq.txt
F = fopen('pnseq.txt', 'w');
fprintf(F, ' %d', a);
fclose(F);

% Строим сигнал из элементов {-1; +1}
for k = 1 : M
  if a(k) == 0 signal(k) = 1;
  else signal(k) = -1;
  end;
end;
% Cогласованный с сигналом фильтр
filter = conj(signal);
% Cтроим матрицу А для расчета АКФ
for k = 1 : M
   for l = 1 : M
       A(k, l) = signal(k) * filter(l);
   end;
end;
% Цикл расчета вектора АКФ
for z = -M : M
   index_time = z + M + 1;
   times(index_time) = z;
   K(index_time) = sum(diag(A, z));
end;
% Вывод графика АФК на экран
time_shift = [-M : M];
plot(time_shift, K);
title('Autocorrelation function of PN-sequence');
xlabel('tau'); ylabel('R(tau)')
grid;