【文章內(nèi)容簡介】 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%將生成的偽隨機(jī)值存放到對應(yīng)的矩陣中%%%%%%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(i1)) y1(k)=x1(i)。 t1(k)=t(i)。 　　k=k+1。 else y1(k)=x1(i1)。 y1(k+1)=x1(i)。 t1(k)=t(i)。 t1(k+1)=t(i)+。 %如果該時刻階躍變化， k=k+2。 end if (x2(i)==x2(i1)) y2(m)=x2(i)。 t2(m)=t(i)。 m=m+1。 else y2(m)=x2(i1)。 y2(m+1)=x2(i)。 t2(m)=t(i)。 　　 t2(m+1)=t(i)+。 　m=m+2。 end if (x3(i)==x3(i1)) y3(n)=x3(i)。 t3(n)=t(i)。 n=n+1。 else y3(n)=x3(i1)。 y3(n+1)=x3(i)。 t3(n)=t(i)。 t3(n+1)=t(i)+。 n=n+2。 endend%%%%%%畫圖說明產(chǎn)生的三組偽隨機(jī)碼%%%%%subplot(3,1,1)。plot(t1,y1,′y′)。grid on。axis(［1,1001,］)。subplot(3,1,2)。plot(t2,y2,′y′)。grid on。axis(［1,1001,］)。subplot(3,1,3)。plot(t3,y3,′y′)。grid on。axis(［1,1001,］)。 whitebg(′black′)。%將三組偽隨機(jī)值存放到新矩陣并傳送到Workspace中， 方便以后調(diào)用數(shù)據(jù)%simin1=［t1。y1］′。simin2=［t2。y2］′。simin3=［t3。y3］′。 輸入上面的文本后，運(yùn)行得到如圖529所示的偽隨機(jī)碼。 其結(jié)果還可以通過Workspace查看，之所以要傳送到Workspace中， 是因為我們在后面的仿真中還要用到這些偽隨機(jī)碼。 圖 模擬三組偽隨機(jī)碼的仿真示意圖 獲取導(dǎo)航信息的仿真在這個仿真中，我們要得到衛(wèi)星的導(dǎo)航信息，其中忽略了頻率偏移，且衛(wèi)星信號的數(shù)據(jù)是基于前面章節(jié)中得到的衛(wèi)星數(shù)據(jù)。在運(yùn)行程序之前， 。 clear。 %清除所有的變量clc。 %關(guān)閉所有指令窗close all。 %關(guān)閉所有未隱藏的窗口load % %n=length(GPSsignal)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%SvNum=12。%， 獲得衛(wèi)星編號為12的C/A碼Temp=fGenerateCAcode3(SvNum)。index1=find(Temp==0)。 %找出C/A碼中的低電平， 形成列向量存放在index1中Temp(index1)=ones(1,length(index1))。 %返回一個與index1相同長度的一維數(shù)組 SinWave=sin(［0:2*pi/8:2*pi*7/8］)。 %產(chǎn)生一個步進(jìn)2*pi/8的正弦波SinWave=single(SinWave)。 %將SinWave矩陣轉(zhuǎn)換成單精度矩陣GpsMatch=zeros(1,1)。 %生成一個1*1的零矩陣SinWave=［SinWave SinWave SinWave SinWave SinWave］。  　%生成一個5 ms的連續(xù)信號 for i=1:length(Temp) GpsMatch=［GpsMatch Temp(1,i)*SinWave］。 end GpsMatch=GpsMatch(2:length(GpsMatch))。 %生成本地匹配信號并去掉第一個元素n=length(GpsMatch)。m=50000。for i=1:m Res(i)=GpsMatch*GPSsignal(1,i:i+n1)′。  　　　%將本地生成匹配信號GpsMatch與GPSsignal循環(huán)卷積 endplot(1:m,Res)。 %圖示卷積的幅值%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Res2=abs(Res)。for i=1:100 %該循環(huán)目的是找出Res中的最大值［C I］=max(Res2)。Res2(1,I)=0。Index(1,i)=I。end。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%下面是要找到GPS signal中的導(dǎo)航電文% %方法是檢驗GPS signal(已被C/A code調(diào)制過)的每個最大值處是否具有相位偏移%w=length(GpsMatch)。m=1。for i=Index(1,1):w:(length(GPSsignal)w+1) NavigationBit(m)=(GPSsignal(i:i+w1)*GpsMatch′)/w。 m=m+1。endNavigationCode=zeros(1,1)。NavigationCode(1,1)=NavigationBit(1,1)。m=2。Count=0。 for i=2:length(NavigationBit) if(NavigationBit(1,i)～=NavigationBit(1,i1)) 　　NavigationCode(1,m)=NavigationBit(1,i)。 m=m+1。 Count=0。 else Count=Count+1。 if(Count=5) 　NavigationCode(1,m)=NavigationBit(1,i)。 m=m+1。 Count=0。 end endend %%%%%我們已得到導(dǎo)航碼NavigationCode,下面將它轉(zhuǎn)化成數(shù)字%%%%NavigationCode=NavigationCode./abs(NavigationCode)。index1=find(NavigationCode0)。NavigationCode(index1)=zeros(1,length(index1))。Table1=［0 0 0 0。 　 0 0 0 1。 　0 0 1 0。 　0 0 1 1。 　 0 1 0 0。 　0 1 0 1。 　 0 1 1 0。 　　　　 0 1 1 1。 1 0 0 0。 1 0 0 1。 1 0 1 0。 1 0 1 1。 1 1 0 0。 1 1 0 1。 1 1 1 0。 1 1 1 1。］。 Result=″。 for i=1:4:length(NavigationCode)  TT=NavigationCode(i:i+3)。 MatchTable=Table1(:,1:4)［TT。 TT。 TT。 TT。 TT。 TT。TT。 TT。 TT。TT。 TT。 TT。TT。 TT。 TT。TT］。 MatchTable=sum(abs(MatchTable)′)。 TempChar=′0′。 BestMatch=find(MatchTable==min(MatchTable))。 %用查表的方式將導(dǎo)航碼NavigationCode轉(zhuǎn)換成十六進(jìn)制數(shù) switch (BestMatch) case 1  TempChar=′0′。  case 2  TempChar=′1′。 case 3  TempChar=′2′。 case 4  TempChar=′3′。 case 5  TempChar=′4′。 case 6 TempChar=′5′。 case 7  TempChar=′6′。 case 8 TempChar=′7′。 case 9 TempChar=′8′。 case 10 TempChar=′9′。 case 11 TempChar=′a′。 case 12 TempChar=′b′。 case 13 TempChar=′c′。 case 14 TempChar=′d′。 case 15 TempChar=′e′。 case 16 TempChar=′f′。 otherwise end  Result=strcat(Result, TempChar)。 %將Result與TempChar兩行的字符連接到一起 end Num1=hex2num(Result(1:16))。  %將十六進(jìn)制數(shù)轉(zhuǎn)換成二進(jìn)制數(shù)，存放在Num1中，