gpssignal.m

fGenerateNavigationData.m C_A_code_generator.m GPSsignal.m 依次产生导航电文,C/A码,GPS最终数据

    x1=10.23;y1=6.24;z1=17.8;    %假设的卫星位置参数
    Code1=fGenerateNavigationData(x1,y1,z1);%产生导航电文,通过调用子函数fGenerateNavigationData实现
    %将产生的Code1中的0找出转换为-1
    index1=find(Code1==0);
    Code1(index1)=-ones(1,length(index1));
    SvNum=12;            %设定卫星编号为12 
    Code2=zeros(1,1);    %定义Code2的初值为0
    Temp=fGenerateCAcode3(SvNum); %将编号为SvNum的卫星通过调用子函数fGenerateCAcode3生成C/A码
    %将Temp中的0找出并转换为-1
    index1=find(Temp==0);
    Temp(index1)=-ones(1,length(index1)); 
    Temp=[Temp Temp];
    %生成Code2
    for  i=1:length(Code1)
        Code2=[Code2  Code1(1,i)*Temp];
    end
    Code2=Code2(2:length(Code2));
    %每位数据通过正弦波来调制
    SinWave=sin([0:2*pi/8:2*pi*3/8]);
    SinWave=single(SinWave);
    GPSsignals=zeros(1,1);
    SinWave=[SinWave  SinWave];
      for i=1:length(Code2)
        GPSsignals=[GPSsignals  Code2(1,i)*SinWave];  
      end
      whos;
     GPSsignals=GPSsignals(2:length(GPSsignals));
     figure(2);
     plot(GPSsignals);axis([100,4000,-1.5 1.5]);
     figure(3);
     plot(10*log10(abs(fft(GPSsignals))));