【正文】 字、數(shù)字或是下劃線。數(shù)組多由一維元素構成，而矩陣中多維元素組成。則以下語句執(zhí)行后可得到有關x和y的圖形： x=1inspace（0，2 * pi, 20）； % 設定x分別為0、2p/2* 2p/3*2p/…2p y1＝sinx，y2＝cosx；％ y1，y2分別是與x對應的正弦和余弦值 plot（x, y1, x，y2）；％在同一坐標圖上分別繪制正弦和余弦曲線如果想分幾次在同一坐標圖上繪制不同的曲線，可使用hold命令： hold on；％保持坐標圖不變，后繪制的圖形疊加在原圖上 hold off；％解除對原圖的保持，將原圖清除后再繪制新圖plot命令的基本格式是：plot（x數(shù)組, y數(shù)組,‘顏色圖標’），如需要在同一圖中畫多根曲線，只需依照此基本格式往后追加其他的x和y的數(shù)組即可。3 GPS衛(wèi)星導航信號算法及其MATLAB仿真 C/A碼仿真代碼及其仿真結果 ，其算法流程歸納如下：1. 用1111111111初始化移位寄存器G1和G2；2. 根據(jù)不同的衛(wèi)星編號選擇不同的S1和S2的位置；3. 分別計算G1和G2⊕操作的和，反饋在移位寄存器的第一位，同時計算G1與G2i當前輸出的⊕結果，作為當前歷元的輸出；4. 返回到步驟2； 根據(jù)算法流程，設計MATLAB函數(shù)：function Self_CA_generate(a) %C/A碼產(chǎn)生 a是衛(wèi)星號（完整代碼見附錄C）例如：賦值a=1，產(chǎn)生1號衛(wèi)星的C/A碼前100個碼片的仿真結果： 1號衛(wèi)星的C/A碼仿真賦值a=16，產(chǎn)生16號衛(wèi)星的C/A碼前100個碼片的仿真結果： 16號衛(wèi)星的C/A碼仿真賦值a=24，產(chǎn)生24號衛(wèi)星的C/A碼前100個碼片的仿真結果：賦值a=32，產(chǎn)生32號衛(wèi)星的C/A碼前100個碼片的仿真結果： P碼的仿真代碼及其仿真結果根據(jù)第一章對P碼原理的分析，繪制如下P碼產(chǎn)生模型： P碼產(chǎn)生模型P碼由四個12位寄存器產(chǎn)生。 根據(jù)其發(fā)生器流程和寄存器間的邏輯結構，編寫MATLAB函數(shù)：function pcode=pcode_Generate1(a,NumberPCode,NumberShift)%a是衛(wèi)星號偏移，NumberPCode 是P碼碼長，NumberShift是P碼碼偏移（詳細源代碼見附錄D） 例如：賦值a=1，NumberPCode=20，NumberShift=10產(chǎn)生1號衛(wèi)星的P碼1030中的20個碼片的仿真結果：ans=1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1號衛(wèi)星的P碼 賦值a=16，NumberPCode=100，NumberShift=0產(chǎn)生16號衛(wèi)星的前100個P碼碼片的仿真結果： 16號衛(wèi)星的P碼 賦值a=24，NumberPCode=100，NumberShift=0產(chǎn)生24號衛(wèi)星的前100個P碼碼片的仿真結果： 賦值a=32，NumberPCode=100，NumberShift=0產(chǎn)生32號衛(wèi)星的前100個P碼碼片的仿真結果： GPS系統(tǒng)采用典型的CDMA體制，系統(tǒng)以碼分多址形式區(qū)分各衛(wèi)星信號。這些分析對捕獲提供了理論基礎。 本文學習了全球衛(wèi)星導航定位系統(tǒng)，重點研究GPS信號結構，設計了GPS信號M碼和P碼的生成算法，簡要介紹C/A碼；。更中要的是，在研究過程中，從隊信號處理、通信嚴厲等基礎理論的膚淺理解，到現(xiàn)在對擴頻同喜系統(tǒng)的原理和性能有了更深一步了解，熟悉了基本信號的處理方法，掌握了GPS信號實現(xiàn)的技術，培養(yǎng)了較強的自學能力和獨立工作能力，為今后的發(fā)展和從事相關行業(yè)的工作奠定了夯實的基礎。感謝四年里各位任課老師和同學在學習、工作和生活中對我的幫助，感謝對我關懷備至的父母和無私幫助我的室友同學，感謝網(wǎng)絡帶來的好處。m=ones(1,3)。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。z=rd。rb=ones(1,10)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。%ra=ra39。 %第三顆衛(wèi)星ra=ones(1,10)。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。endrd=rd39。endif a==4。rc(1)=[mod(rb(5)+rb(9),2)]。 rc(1)=[mod(rb(5)+rb(9),2)]。plot(z)save ca_code z。m=ones(1,3)。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。z=rd。rb=ones(1,10)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。%ra=ra39。 %第七顆衛(wèi)星ra=ones(1,10)。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。endrd=rd39。endif a==8。rc(1)=[mod(rb(2)+rb(9),2)]。 rc(1)=[mod(rb(2)+rb(9),2)]。plot(z)save ca_code z。m=ones(1,3)。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。z=rd。rb=ones(1,10)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。%ra=ra39。 %第十一顆衛(wèi)星ra=ones(1,10)。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。endrd=rd39。endif a==12。rc(1)=[mod(rb(5)+rb(6),2)]。 rc(1)=[mod(rb(5)+rb(6),2)]。plot(z)save ca_code z。m=ones(1,3)。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。z=rd。rb=ones(1,10)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。%ra=ra39。 %第十五顆衛(wèi)星ra=ones(1,10)。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。endrd=rd39。endif a==16。rc(1)=[mod(rb(9)+rb(10),2)]。 rc(1)=[mod(rb(9)+rb(10),2)]。plot(z)save ca_code z。m=ones(1,3)。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。z=rd。rb=ones(1,10)。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。rc(1)=[mod(rb(2)+rb(5),2)]。endif a==18。endrd=rd39。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 %第十七顆衛(wèi)星ra=ones(1,10)。%ra=ra39。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。rb=ones(1,10)。z=rd。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。m=ones(1,3)。plot(z)save ca_code z。 rc(1)=[mod(rb(7)+rb(8),2)]。rc(1)=[mod(rb(7)+rb(8),2)]。end if a==14。endrd=rd39。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 %第十三顆衛(wèi)星ra=ones(1,10)。%ra=ra39。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。rb=ones(1,10)。z=rd。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。m=ones(1,3)。plot(z)save ca_code z。 rc(1)=[mod(rb(2)+rb(3),2)]。rc(1)=[mod(rb(2)+rb(3),2)]。endif a==10。endrd=rd39。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 %第九顆衛(wèi)星ra=ones(1,10)。%ra=ra39。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。rb=ones(1,10)。z=rd。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。m=ones(1,3)。plot(z)save ca_code z。 rc(1)=[mod(rb(2)+rb(10),2)]。rc(1)=[mod(rb(2)+rb(10),2)]。endif a==6。endrd=rd39。for i=1:1023 rd(i)=mod(ra(10)+rc(1),2)。 %第五顆衛(wèi)星ra=ones(1,10)。%ra=ra39。 ra=[mod(ra(3)+ra(10),2) ra(1:9)]。rb=ones(1,10)。z=rd。 rb=[mod(rb(2)+rb(3)+rb(6)+rb(8)+rb(9)+rb(10),2) rb(1:9)]。m=ones(1,3)。plot(z)save ca_code z。 rc(1)=[mod(rb(3)+rb(7),2)]。rc(1)=[mod(rb(3)+rb(7),2)]。end if a==2。endrd=rd39。for i=1:100 rd(i)=mod(ra(10)+rc(1),2)。 %第一顆衛(wèi)星ra=ones(1,10)。在整個畢業(yè)設計過程中，我查閱了大量的相關書籍資料，對GPS全球定位系統(tǒng)有了一個較為全面的了解，提高了自己對MTATLAB軟件的使用能力。⑵通過對GPS信號的結構和特性，直接序列擴頻通信等理論的深入研究和分析，確定了GPS信號生成算法，編寫相應的MATLAB程序，實現(xiàn)GPS信號的產(chǎn)生，同時驗證其相關特性。因此一直都是過內(nèi)外GNSS研究領域的熱門話題。P碼互相關特性分析的數(shù)據(jù)說明，不同的P碼序列具有很小的相關性，其相關峰峰值與同一段P碼序列在無相位偏移時的相關峰峰值相差一個數(shù)量級。各個寄存器的抽頭系數(shù)如下式所示：x1a=[x6+x8+x11+x12]； x1b=[x1+x2+x5+x8+x9+x10+x11+x12]；x2a=[x1+x3+x4+x5+x7+x8+x9+x10+x11+x12]； ()x1b=[ x2+x3+x4+x8+x9 +x12]；這些寄存器的初始向量如下：x1aa=[0 0 0 1 0 0 1 0 0 1 0 0]；x1bb=[0 0 1 0 1 0 1 0 1 0 1 0]；x2aa=[1 0 1 0 0 1 0 0 1 0 0 1]；x2bb=[0 0 1 0 1 0 1 0 1 0 1 0]；每個寄存器循環(huán)次數(shù)分別為:x1a=3750，x1b=3749，x2a=3750，x2b=3749。grid命令用來在圖形上添加或者刪除網(wǎng)格線，它是一個切換命令，若第一次是添加，則再執(zhí)行一次就是刪除。 MATLAB特殊函數(shù)函數(shù)符號說 明zero (i,j)創(chuàng)建i行j列的全零矩陣ones (i,j)創(chuàng)建i行j列的全1矩陣eye (i,j)創(chuàng)建i行j列對角線為1的矩陣rand (i,j)創(chuàng)建i行j列的隨機矩陣 MATLAB是基本的繪圖命令有二維曲線繪圖命令plot和三維曲線繪圖命令plot3。 MATLAB所定義的特殊變量及其意義變量名意 義help在線幫助，如help quitwho列出所有定義過的變量名稱ans默認的用來表示計算結果的變量名eps極小值=pip值inf無窮大的數(shù)165。邏輯與運算= =關系運算，相等|邏輯或運算不等于邏輯非運算小于xor邏輯異或運算=小于等于……續(xù)行標志大于，分行符，結果不顯示=大于等于；分行符，結果顯示%注釋標志’矩陣轉置.’向量轉量MATLAB可以將計算結果以不同的精度輸出，列表說明如下： 數(shù)據(jù)輸出格式