GPS C/A Code Signal Structure - Fundamentals of Global Positioning System Receivers: A Software Approach - page 102

%p51.m generates MLS and G2 outputs and check their delay
time
% input
% k1, k2: positions of two taps
% k3: delay time
%
k=input('enter [k1 k2 k3] = ');
inp1=-ones(1,10); % initial condition of register
for j=1:1023;
mlsout(j)=inp1(10); % MLS output
modulo=inp1(2)*inp1(3)*inp1(6)*inp1(8)*inp1(9)*inp1(10);
inp1(2:10)=inp1(1:9);
inp1(1)=modulo;
g2(j)=inp1(k(1))*inp1(k(2)); % G2 output
end
if mlsout==g2([k(3):1023 1:k(3)-1])
disp('OK')
else
disp('not match')
end