generatecode.m

gps c/a码程序以及模拟卫星定位的源代码、地球模型

close all;
clear;
clc;
t=0:5:1000;


x1=rand(1,length(t));
x1=round(x1);
index1=find(x1==0);
x1(index1)=-ones(1,length(index1));


x2=rand(1,length(t));
x2=round(x2);
index1=find(x2==0);
x2(index1)=-ones(1,length(index1));


x3=rand(1,length(t));
x3=round(x3);
index1=find(x3==0);
x3(index1)=-ones(1,length(index1));


y1(1)=x1(1);y2(1)=x2(1);y3(1)=x3(1);
t1(1)=t(1);t2(1)=t(1);t3(1)=t(1);
k=2;m=2;n=2;

for i=2:length(t)
    if (x1(i)==x1(i-1))
        y1(k)=x1(i);
        t1(k)=t(i);
        k=k+1;
        
    else
        y1(k)=x1(i-1);
        y1(k+1)=x1(i);
        t1(k)=t(i);
        t1(k+1)=t(i)+0.01;
        k=k+2;
    end
    if (x2(i)==x2(i-1))
        y2(m)=x2(i);
        t2(m)=t(i);
        m=m+1;
        
    else
        y2(m)=x2(i-1);
        y2(m+1)=x2(i);
        t2(m)=t(i);
        t2(m+1)=t(i)+0.01;
        %t2(k)=t(i);
        %t2(k+1)=t(i)+0.01;
        m=m+2;
    end
    if (x3(i)==x3(i-1))
        y3(n)=x3(i);
        t3(n)=t(i);
        k=k+1;
        
    else
        y3(n)=x3(i-1);
        y3(n+1)=x3(i);
        t3(n)=t(i);
        t3(n+1)=t(i)+0.01;
        n=n+2;
    end
end

subplot(3,1,1);
plot(t1,y1,'y');
grid on;
axis([-1,1001,-1.5,1.5]);


subplot(3,1,2);
plot(t2,y2,'y');
grid on;
axis([-1,1001,-1.5,1.5]);

subplot(3,1,3);
plot(t3,y3,'y');
grid on;
axis([-1,1001,-1.5,1.5]);

whitebg('black');

simin1=[t1;y1]';
simin2=[t2:y2]';
simin3=[t3:y3]';