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))));