реализация Р-кода GPS

clear
clc
nul=zeros(1,4083); 
InitX1A=[0 0 1 0 0 1 0 0 1 0 0 0];%кодовый вектор инициализации регистра Х1А
InitX1B=[0 1 0 1 0 1 0 1 0 1 0 0];%кодовый вектор инициализации регистра Х1В
X1A=[nul InitX1A];% массивы для последовательностей Х1А и Х1В 
X1B=[nul InitX1B];
% формирую первые 4095 бит Х1А и Х1В 
for t = 4083:-1:1 
 x1a=X1A(t+6)+X1A(t+8)+X1A(t+11)+X1A(t+12); 
 x1a=mod(x1a,2);
 X1A(t)=x1a;
 x1b=X1B(t+1)+X1B(t+2)+X1B(t+5)+X1B(t+8)+X1B(t+9)+X1B(t+10)+X1B(t+11)+X1B(t+12); 
 x1b=mod(x1b,2);
 X1B(t)=x1b;
end
% обрезаю массивы,формируя первые кадры Х1А и Х1В 
X1A=X1A(4:4095);
X1B=X1B(3:4095);
% создаю последовательность из 15345000 бит, т.е. 3750 кадров Х1А 
for i=1:1:12
    X1A=[X1A X1A];
end
X1A=X1A(1:15345000);
% создаю последовательность из 15344657 бит, т.е. 3749 кадров Х1В 
for i=1:1:12
    X1B=[X1B X1B];
end
X1B=X1B(344:15345000);
% смотрю какое значение на выходе регистра Х1В в конце 3749 кадра, это нолик
z=X1B(1);
z=zeros(1,343);
% добавляю ноликов до 15345000 в последовательность X1B
X1B=[z X1B];
% суммирую две последовательности, получаю эпоху X1
X1=(X1A+X1B);
X1=mod(X1,2);
%аналогично формирую эпоху X2i, с той лишь разницей что последовательность
%Х2 циклически сдвигаю на единичку
nul=zeros(1,4083); 
InitX2A=[1 0 0 1 0 0 1 0 0 1 0 1];
InitX2B=[0 1 0 1 0 1 0 1 0 1 0 0];
X2A=[nul InitX2A];
X2B=[nul InitX2B];
for t = 4083:-1:1 
 x2a=X2A(t+1)+X2A(t+3)+X2A(t+4)+X2A(t+5)+X2A(t+7)+X2A(t+8)+X2A(t+9)+X2A(t+10)+X2A(t+11)+X2A(t+12); 
 x2a=mod(x2a,2);
 X2A(t)=x2a;
 x2b=X2B(t+2)+X2B(t+3)+X2B(t+4)+X2B(t+8)+X2B(t+9)+X2B(t+12); 
 x2b=mod(x2b,2);
 X2B(t)=x2b;
end
X2A=X2A(4:4095);
X2B=X2B(3:4095);
for i=1:1:12
    X2A=[X2A X2A];
end
X2A=X2A(1:15345000);
for i=1:1:12
    X2B=[X2B X2B];
end
X2B=X2B(344:15345000);
z=X2B(1);
z=zeros(1,343);
X2B=[z X2B];
X2=(X2A+X2B);
X2=mod(X2,2);
% сдвигаю Х2 на единицу
X2i=circshift(X2,[0 +1]);
%складываю две эпохи, 
P=X1+X2i;
P=mod(P,2);