【文章內(nèi)容簡介】 round(x1)。 %將 x1 湖北大學(xué)本科畢業(yè)論文（設(shè)計(jì)） index1=find(x1==0)。 %找出 x1數(shù)組中為 0的元素 x1(index1)=ones(1,length(index1))。 %返回一個 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=zeros(1,1)。 y2=zeros(1,1)。 y3=zeros(1,1)。 %生成三個 1*1的零矩陣，存放偽 t1=zeros(1,1)。 t2=zeros(1,1)。 t3=zeros(1,1)。 %生成三個 1*1 的零矩陣，存放時 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%將生成的偽隨機(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。 湖北大學(xué)本科畢業(yè)論文（設(shè)計(jì)） 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中， 是因?yàn)槲覀冊诤竺娴姆抡嬷羞€要用到這些偽隨機(jī)碼。 湖北大學(xué)本科畢業(yè)論文（設(shè)計(jì)） 圖 模擬三組偽隨機(jī)碼的仿真示意圖 獲取導(dǎo)航信息的仿真 在這個仿真中，我們要得到衛(wèi)星 的導(dǎo)航信息，其中忽略了頻率偏移，且衛(wèi)星信號的數(shù)據(jù)是基于前面章節(jié)中得到的衛(wèi)星數(shù)據(jù)，存放在 。在運(yùn)行程序之前， 我們必須將 Workspace中。 clear。 % clc。 % close all。 % load %從 Workspace中 %n=length(GPSsignal)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% SvNum=12。 %調(diào)用 ， 獲得衛(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 湖北大學(xué)本科畢業(yè)論文（設(shè)計(jì)） SinWave=single(SinWave)。 %將 SinWave GpsMatch=zeros(1,1)。 %生成一個 1*1 SinWave=［ SinWave SinWave SinWave SinWave SinWave］ 。 %生成一個 5 ms 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)航電文 % %方法是檢驗(yàn) 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。 湖北大學(xué)本科畢業(yè)論文（設(shè)計(jì)） 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 switch (BestMatch) case 1 TempChar=′ 0′ 。 case 2 TempChar=′ 1′ 。 case 3 TempChar=′ 2′ 。 case 4 湖北大學(xué)本科畢業(yè)論文（設(shè)計(jì)） 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