【正文】 ,cube,cube,1)。 tempx=0。 tempy=6400。 tempz=0。 cube=100。 boxplot3(tempx,tempy,tempz,cube,cube,cube,2)。 tempx=0。 tempy=0。 tempz=6400。 cube=100。 boxplot3(tempx,tempy,tempz,cube,cube,cube,4)。 endfunction drawsatellite(moveX,moveY,moveZ,color)。Length=500。Width=500。Height=500。boxplot3(moveX,moveY,moveZ,Length,Width,Height,color)。boxplot3(moveX,moveY,moveZ,Length,Width*5,Height/10,color)。 它的作用就是畫出一個在上述地球坐標圖當中的衛(wèi)星。30仿真程序一：繪制衛(wèi)星的軌道平面程序見附錄：圖 52 衛(wèi)星軌道平面仿真圖仿真程序二：單顆衛(wèi)星不同時刻的動態(tài)仿真程序見附錄：地球 質(zhì)心31圖 53 單顆衛(wèi)星動態(tài)仿真圖由上圖 10 可知，衛(wèi)星的天空瞬時位置是隨著時間的變動而發(fā)生變化的。對于我們在地球上面的一個用戶來說，一天當中的不同時刻看到的衛(wèi)星是不相同的。這個仿真程序的功能實際上就是仿真了在一個軌道上面的衛(wèi)星在不同的時候（這里以一個小時為一個觀測時元，進行動態(tài)的在屏幕上顯示器位置）通過改變時間可以顯示出不同時間內(nèi)的衛(wèi)星的瞬時位置） ，而且，通過 matlab 當中的三維旋轉(zhuǎn)圖標，我們可以很清晰地從不同角度看到衛(wèi)星和坐標原點（這里用黑點表示的的相對位置的變化） 。哪么，在其他軌道面上的衛(wèi)星的運動也可以類似的模擬出來，這里不在重復！ 仿真程序三：衛(wèi)星在某個時刻的全軌道平面的分布和可見衛(wèi)星程序見附錄：32圖 54 全軌道平面的圖形其在軌坐標分別如下：SatellitePosition = +004 * 33 這個程序仿真了在某個時刻（在程序里面是在時刻 timenow=0）的全部 24顆衛(wèi)星的軌道圖形，以及對于用戶來說在這個時刻可以看到的衛(wèi)星。在程序當中，我們假定了用戶的位置坐標是（6400，3352，5410） 。通過改變用戶的不同位置，可以在同一個時刻看到不同（4 到 11 顆）的可見衛(wèi)星。對于在這個程序當中的可見衛(wèi)星的及時在軌坐標如下：34圖 55 用戶可見的衛(wèi)星分布SatellitePosition = +004 * 0 0 0 0 0用戶位置地心位置35 0 0 0 0 0 0 0 0 0 仿真程序四：用可見衛(wèi)星計算用戶位置程序見附錄多于四顆的衛(wèi)星用最小二乘法逼近計算。對于前面的假設用戶位置（6400，3352，5410）進行計算，結(jié)果如下：36圖 56 用戶位置的計算由圖 13 可見，根據(jù)最小二乘法的原理，計算出來了用戶在一個時刻的位置坐標，并把它表示在了這個天球坐標系當中。calculaterecord = +003 * 0 0 0 先后經(jīng)過了六次的迭代算法，吧用戶的計算位置一步一步的逼近了用戶的實際位置，根據(jù)部同的精度要求，我們運算的量的大小也有不同。這點，可以根據(jù)程序當中的參數(shù) Error 的設定而有所出入。圖 57 用戶位置的計算通過對用戶仿真計算的進一步放大，可以很清楚地看到：在每一次迭代的過程當中，我們都實時的把每一步迭代的值也表示在了同一個坐標軸里面：目的就是很直觀的反應出我們基于最小二乘法的偽距算法的主要思路。如下圖所37示：圖 58 用戶位置的計算 上圖 15 當中的黑色箭頭從上到下一次表示的是第一次到第六次計算出來的用戶的位置，第七箭頭表示的用戶的實際位置（用白色的柱體表示） 。38結(jié) 論通過這次的設計，學習到了很多，因為 GPS 衛(wèi)星運動及定位以前沒有接觸過，所以本次設計難度很大，但在最后都還是通過查找資料或者老師同學幫助也都一一克服。GPS 衛(wèi)星是一個很復雜的系統(tǒng)，其內(nèi)容很多，比如它的坐標系、時間系統(tǒng)等，其坐標系分為天球坐標和地球坐標。天球坐標描述衛(wèi)星的運動位置和狀態(tài)，地球系統(tǒng)則描述地面測站的位置。GPS 衛(wèi)星作為一個高空動態(tài)已知點，其位置是隨時間不斷變化的，因此，在給出衛(wèi)星運行位置的同時，必須給出相應的瞬時時刻。所以時間系統(tǒng)是 GPS 衛(wèi)星計算里面很重要的一環(huán)。要計算衛(wèi)星的運動軌跡及位置，所需要的參數(shù)很多，衛(wèi)星在空間繞地球運動時，除了受到地球重力場的引力作用外，還受到太陽、月亮和其他天體引力以及太陽光壓、大氣阻力和地球潮汐等因素影響。通常把作用于衛(wèi)星上的各種力，按其影響的大小分為兩類：中心力和攝動力。因為此次設計要求的精度不高，所以在設計中攝動力忽略不計。衛(wèi)星在地球引力場中無攝運動，也稱為開普勒運動。所以本次設計中主要用開普勒定律計算其位置軌跡，在運算中還加入了最小二乘法計算參數(shù)數(shù)據(jù)。此次設計中衛(wèi)星的動態(tài)仿真，計算中主要遇到的問題是對于特定軌道的衛(wèi)星在不同時刻的近地點幅角的計算。所以，根據(jù)衛(wèi)星軌道的星座參數(shù)，找出在某時刻的近地點幅角（一般是零時刻） ，然后在此基礎上以時間為基礎進行疊加，確定出每個時間點的瞬時近地點幅角值?；诿嫦?qū)ο蠛统绦蛟O計思想，利用 matlab 編程了 GPS 定位仿真程序，該程序從自定義的星座數(shù)據(jù)文件中建立起衛(wèi)星星座。設計一個子函數(shù)，命名為。其功能為其他程序得出一個全局變量 satelliteposotion，把它和其他程序放入同一個文件夾中，以便其他程序使用。此次設計程序主要包含四個模塊：1 繪制衛(wèi)星的軌道平面，以便于直觀了解其軌道；2 單顆衛(wèi)星的動態(tài)仿真，了解其不同時刻的位置變化，通過 matlab 當中的三維旋轉(zhuǎn)圖標，可以很清晰的從不同角度看到衛(wèi)星和坐標原點；3 衛(wèi)星在某時刻的全軌道平面分布和可見衛(wèi)星，該程序仿真了在某時刻的全部 24 顆衛(wèi)星的軌道圖形，以及對于用戶來說在該時刻可見衛(wèi)星位置；394 利用可見衛(wèi)星計算用戶位置，利用最小二乘法原理，計算用戶某一時刻的位置坐標。通過 matlab 仿真衛(wèi)星運動及軌跡，對其有了更加深入的認識和了解，對GPS 也有了更直觀的認識，對于以后遇到 GPS 的應用時有了一個簡單的基礎。40致謝在論文完成之際，我懷著激動和喜悅的心情向所有關(guān)心和幫助過我的老師和同學表示萬分的感謝！首先要感謝的是幫助我完成本次設計的葛青老師。在本次設計中，我遇到了許多知識和技術(shù)方面的疑問，葛老師都一一幫我解答，并提出了許多具有創(chuàng)新意義的建議。她淵博的知識，嚴謹?shù)墓ぷ鲬B(tài)度，對工作盡心盡職的精神和平易近人的品格都深深感染了我，將會讓我受益終身。在此，我向葛老師表示深深的感謝，并致以崇高的敬意！同時還要感謝我的班主任孟玲華老師、周偉老師。大學四年期間，在大一大二時孟老師作為我們班主任，她給以我無私的關(guān)懷，是她幫助我成長，讓我從一個無知的高中生脫變到具有獨立思考能力的大學生，并提供給我許多鍛煉自身能力的機會。在大三大四的時候周偉老師出任我們班主任，他作為一個專業(yè)課老師，給了我在學習上很大的幫助。在后此外還要感謝所有教育和關(guān)懷過我的老師們，還要感謝成都理工大學核技術(shù)與自動化工程學院的領(lǐng)導們四年來給予的多方面的關(guān)心和幫助! 感謝所有的朋友和同學，讓我度過了一個快樂而又充實的大學時光！最后要感謝我的父母和家人，他們總是盡最大能力提供給我最好的條件，讓我無憂無慮的去學習。是他們的鼓勵和支持，是他們深深的關(guān)愛，讓我得以順利的完成了學業(yè)。成都理工大學畢業(yè)設計( 論文)41參考文獻[1] 周忠謨, 衛(wèi)星測量原理與應用 [M] .北京:測繪出版社,1997.[2] 李征航, 測量與數(shù)據(jù)處理 [M] .武漢:武漢大學出版社,2022.[3] Elliott . Understanding GPS: principles and applications [M] .邱致和,王萬義,京:電子工業(yè)出版社,2022.[4] 衛(wèi)星全球定位系統(tǒng)簡介[J]. 科技經(jīng)濟市場 2022（8）[5] 郭秋英, 衛(wèi)星坐標的計算 [J] .全球定位系統(tǒng) ,2022 (3) .[6] 夏振純,金志華, MATLAB 的 GPS 選星分析及實驗研究 [J] .計算機仿真,2022,22(1) :5154.[7] 楊俊， 基本原理及其 Matlab 仿真[M].西安：西安電子科技大學出版社[8]張勤， 測量原理及應用[M].北京：科學出版社[9] 衛(wèi)星定位仿真分析與系統(tǒng)程序設計[D]. 同濟大學工學碩士學位論文[10]劉瓊。岳昊， 在 GPS 衛(wèi)星幾何分布計算中的應用[C]. 全國測繪科技信息網(wǎng)中南分網(wǎng)第二十四次學術(shù)信息交流會論文集[11] 衛(wèi)星定位系統(tǒng)簡介[Z].陜西國防工業(yè)職業(yè)技術(shù)學院學報第二十卷第一期。[12]郭秋英， 衛(wèi)星坐標的計算[J].2022[13] GPS 及在現(xiàn)代交通中的應用[J].[14]樓益棟，施闖，葛茂榮， 衛(wèi)星實時精密定軌及初步結(jié)果分析[J].[15] GPS 衛(wèi)星定位系統(tǒng)，提高車輛利用率[J].[16]趙齊樂. GPS 導航星座及低軌衛(wèi)星的精密定軌理論和軟件研究[D] .武漢: 武漢大學 ,2022[17]Chen K. Real2Time Precise Point Positioning andIt s Potential Applications [ C] . ION GNSS 17th In2ternational Technical Meeting of the Satellite Divi2sion , Long Beach , CA ,2022[18]Liu Jingnan , Ge Maorong. PANDA Sof tware andIt s Preliminary Result of Positioning and Orbit Det2ermination[J ] . Wuhan University Journal of NaturalSciences , 2022 ,8 (B2) : 603609[19] Kouba J , Heroux P. GPS Precise Point PositioningUsing IGS Orbit Projuct s[J ] . GPS Solution , 2022 ,5 (2) : 1228[20]成都理工大學畢業(yè)設計( 論文)42附錄程序一（繪制衛(wèi)星的軌道平面）代碼function drawsatelliteorbita=26560。% 衛(wèi)星軌道的長半軸.e=。%e是橢圓的偏心率.E=[0::2*pi]。x=a*(cos(E)e)。y=a*sqrt((1e^2))*sin(E)。z=0*E。DtoR=2*pi/360。A1=[ ]。%衛(wèi)星星座數(shù)據(jù) .for k=1:6 A=A1(k)*DtoR。%升交點的經(jīng)度 B=55*DtoR。%軌道的傾角 C=pi/100。%近地點幅角 %總共有6個衛(wèi)星軌道平面 R1=[cos(A) sin(A) 0。 sin(A) cos(A) 0。 0 0