GPS衛(wèi)星坐標(biāo)計算

1、衛(wèi)星定位技術(shù)與方法-根據(jù)播送星歷參數(shù)計算衛(wèi)星坐標(biāo)作業(yè)報告指導(dǎo)教師：熊 永良班 級： 測繪二班學(xué)生姓名：段 海東學(xué)生學(xué)號:2 0 0 8 0 7 8 3作業(yè)日期：2021年 12月 08日1 .數(shù)據(jù)2.2 .計算步驟2.1 .平均角速度 (mean angular speed)22 .規(guī)化時亥U (normal time) 3.3 .平近點角(mean anomaly)34 .偏近點角(eccentric anomaly) 35 .真近點角(true anomaly) 36 .升交距角(argument of ascending node) 37 . 軌道向徑(Orbital radius )

2、38 .擾動改正(Perturbed correction) 410 .衛(wèi)星在開交點軌道直角坐標(biāo)系中的坐標(biāo) 411 . 升交點經(jīng)度(Longitude of ascending node ) 53 .源程序5.4 .程序運行結(jié)果1.4七.作業(yè)體會15根據(jù)播送星歷參數(shù)計算衛(wèi)星坐標(biāo)數(shù)據(jù)：根據(jù)以下的播送星歷參數(shù)計算 UTC2004年1月30日8點0分00秒一20分00秒,每隔一分鐘的PRN7的衛(wèi)星坐標(biāo).Compute the coordinate of PRN7 with interval of 1 minute.Navigation data:2NAVITIQII MTARINEX VERSION

3、 / TEPEEFHTOt解£3-L 版本3Q-JAN-04 碌;41 PCM / RUH BY / DATEHMD OF HEADEB4 i i sc, s o a a- a, 5/ ：a21332157D-04-0. eaz 17 55 it? 23> 11 a. oa maoa oaaoo choo D. 13SDOOD0OOQ(nm)3 QB6先町50叩0加02 0.14755430 13TD-O0 (J. LDO13O<7339H)1 D,號g然BD05 Q 朋0026源B360Mz 飆 5口密5309995巾45 0,5L535T77KM口,4B 口 BOOO

4、OOOOmm)6 0. SieeOSFSaGgtD-OT 0. 305054077a39t>HJL 0.仇 966365662 D07 Q,照l旄7例計00 0.27946875000ai>+03-0.10S701071O0M0-0,8L86L98K33D-08 J. i；L_J.L.7_L.liJT, Il QJQ 口 QOWW® 用 I 0 一,一 一二.,二刀匚 CO工加/ D, QDO(XX)(X)(XKOIM)0 Q OQOOTOOOOWaCHQO-O- 60時蹣糕0 AOS 0F 13 90000 OOOOOOX D, QEX)瞅MQtXDM)口7 4 1 a

5、0 9 0 Oi 0 Q62932S3&7313M3-0 54K96鴕1064D-11 QQOQQWQIXKKKJtHOO D,31(JWJQXKMKDM)3 Q 5肝312印叩0研02 0.4350059387771)-09-0,31<9904TfiaT3IHOl 0. 2754B5ZL7571>-0S 0. 125S56591351D-01 0. 90375542640715-05 (J. 5 L5SF3Q 38319(04 口, 4W瞅(XKXXKMHOg UZ067碼611冊處能 0.193198315097W1-0.7254313062001)-07 D, 93了

6、曲094"5網(wǎng)口 0. ISierSODOOdOMa-O. L85137234909DKIL-0. 7SlflS»109M2IH)8 -0,】3醛4"42alM)9 CH口QOOOOOOOQg+Ql 0. L25EOOOOOCK)OOKH Q. 000000000000 1nm D, ODODOODOOWDQHIO 0, OOOOOOOOOOOOIHK-O. Lafi2M514923D-O8 0.2100000000001X)3 D. QQOOCOO3QOOCWQ衛(wèi)星導(dǎo)航文件格式:PRN購/H工skit (iicotid?登二財弓K)IX出比小r radian 曲

7、MXtsifirs)四儂艱t駛際回：不以：這如U sec cfGPStek口神必埼孰礴由國年 Iraier s)而誣皿j|港徊好；OMGA DC 7 (2也 5 穌 tDOT值蒯噌膽；C 湎 ctt L2 diannelGP5AV«k-LiPlitafla?SV anzrzcy .底:gSVhtlrETGDfirafc'33C cfDi：iTstfm：<s：cti at M mt*sa«>fcondispare5pafispine計算步驟 ：The steps for satellite coordinates1 .平均角速度 (mean angular

8、 speed):n n0 nGMn由播送星歷獲得,GM=3.986005e+142 .規(guī)化時刻(normal time):tkt tot0(由播送星歷獲得),t為GP姍秒3 .平近點角(mean anomaly):Mk Mo n tkM0(由播送星歷獲得)4 .偏近,電角(eccentric anomaly) :Ek Mk e sin Ek迭代求解：初始值取E=M以弧度為單位5 .真近點角(true anomaly):1e2 sin EkVkarctancosEk e6 .升交品巨角(argument of ascending node)0 V k 近地點角距(argument of peri

9、gee)7 .軌道向徑(Orbital radius ):rka (1 e cosEk)8 .擾動改正(Perturbed correction ):升交角距(Argument of ascending nodeu= Gc cos2CuS sin 2 0軌道向徑(Orbital radius)r C © CoS 2CrS sin 2 0?軌道頃角(Orbital inclination)i C iC cos 2CiS Sin 2 00 是升交角品巨 (the argument of ascending node)9 .改正后升交角距、軌道向徑、軌道傾角改正后開交角品E ( Corre

10、cted argument of ascending nodu k 0 u改正后的軌道向徑(Corrected orbital radius)rk a (1 e cosEk) r改正后的軌道傾角(Corrected orbital inclination :ik i. i (IDOT)tk10 .衛(wèi)星在開交點軌道直角坐標(biāo)系中的坐標(biāo)：如下列圖所示xk rkcosukYk rksinuk11. 升交,電經(jīng)度(Longitude of ascending node )：如下列圖所示t 0 ( e)tket 0e_ -5e 7.2921151467 10 rad/s12. 在地固坐標(biāo)系中的衛(wèi)星位置 (

11、 Expressed in spheric coordinate systemk Xk cos k Yk cosik sin k Yk Xk sin k Yk cosik cos k Zk Yksinik三.源程序：using System;using System.Collections.Generic;using System.ComponentModel;using System.Data;using System.Drawing;using System.Text;using System.Windows.Forms;using System.IO;using System.Text.

12、RegularExpressions;namespace Test(public partial class Forml : Form(public Form1()(InitializeComponent();listView1.Columns.Add("listView1.Columns.Add("listView1.Columns.Add("listView1.GridLines = true; /listView1.View = View.Details;/listView1.HeaderStyle頭進(jìn)行設(shè)置listView2.Columns.Add(&qu

13、ot;listView2.Columns.Add("xlistView2.Columns.Add("ylistView2.Columns.Add("zlistView2.GridLines = true; /listView2.View = View.Details;/listView2.HeaderStyle頭進(jìn)行設(shè)置序號", 40);星歷內(nèi)容",130);導(dǎo)航數(shù)據(jù)", 130);顯示表格線顯小表格細(xì)節(jié)=ColumnHeaderStyle.Clickable;/ 對表時間,60);坐標(biāo)",150);坐標(biāo)",150

14、);坐標(biāo)",150);顯示表格線顯小表格細(xì)節(jié)=ColumnHeaderStyle.Clickable;/ 對表private void button1_Click(object sender, EventArgs e)(/ 讀取相對路徑string str1 =AppDomain.CurrentDomain.SetupInformation.ApplicationBase;string filename = strl + "navigation data.txt "/讀取衛(wèi)星播送星歷文件StreamReader myreader = new StreamReade

15、r(filename, Encoding.Default);string myinfo = myreader.ReadToEnd();myreader.Close();/ 把衛(wèi)星播送星歷里的D改為e,以便后續(xù)計算string mystring = myinfo;string myinfor = mystring;myinfor = mystring.Replace("D", "e");/把數(shù)據(jù)分開讀入到一個數(shù)組中存儲string split = new string " " ;string arrs = myinfor.Split(s

16、plit, StringSplitOptions.RemoveEmptyEntries);double M = new doublearrs.Length;/ 星歷代碼string N = new stringarrs.Length;N0= "PRN" N1 = "Yer" N2= "Mon" N3 = "day" N4="H" N5 = "M" N6 = "sec" N7 = "a0" N8 = "a1" N9 =

17、 "a2"N10 = "IODE" N11 = "Crs" N12 = "delta-n" N13="M0"N14= "Cuc" N15 = "e" N16 = "Cus" N17 = "sqrt(a)"N18= "t0e"N19 = "Cic" N20 = "omega0" N21 = "Cis"N22 = "i0"

18、 N23 = "Crc" N24 = "omega" N25= "omega-spot"N26 = "IDOT" N27 = "Codes on L2 channel" N28 = "GPSWeek" N29 = "L2 P data flag"N30 = "SV accuracy" N31 = "SV health" N32 = "TGD"N33 = "IODC Issue of Da

19、ta"N34 = "Transmission time of message"/把衛(wèi)星播送星歷讀入到數(shù)組中for (int i = 0; i < arrs.Length; i+)ListViewItem li = new ListViewItem();li.Text = (i + 1).ToString();li.SubItems.Add(Ni);li.SubItems.Add(arrsi);listView1.Items.Add(li);Mi = double.Parse(arrsi);double t0e=460800.00;double t;doubl

20、e XK = new double21;double YK = new double21;double ZK = new double21;for (int l = 0; l < 21;l+ )/計算平均角速度double GM = 3986004.418e008;double n, n0;n0 = Math.Sqrt(GM / (Math.Pow(M17, 6);n = M12 + n0;/規(guī)劃時刻t=t0e+l*60;double tk = t-t0e;double Mk = M13 + n * tk;/迭代計算平近點角的計算double Ek, Ek1;Ek = Mk;Ek1 =

21、Mk + M15 * Math.Sin(Ek);do(Ek = Ek1;Ek1 = Mk + M15 * Math.Sin(Ek);)while (Math.Abs(Ek1 - Ek) > 1e-15);/計算真近點角double Vk = Math.Atan(Math.Sqrt(1 - M15 * M15) * Math.Sin(Ek) / (Math.Cos(Ek) - M15);/反正切值的象限處理double A1 = Math.Atan(Math.Sqrt(1 - M15 * M15) * Math.Sin(Ek);double A2 = Math.Cos(Ek) - M15;

22、if (Vk < 0)(if (A1 < 0 && A2 > 0)Vk += 2 * Math.PI;)if (A1 > 0 && A2 < 0)(Vk += Math.PI;)else(if (A1 <= 0 && A2 <= 0)Vk = Vk + Math.PI;)/計算升交角距double fk = Vk + M24;/攝動改正值的計算doubledu = M14 * Math.Cos(2 * fk) + M16 * Math.Sin(2* fk);doubledr = M23 * Math.Co

23、s(2 * fk) + M11 * Math.Sin(2* fk);doubledi = M19 * Math.Cos(2 * fk) + M21 * Math.Sin(2* fk);/攝動改正double uk = fk + du;doublerk = Math.Pow(M17, 2) * (1 - M15 * Math.Cos(Ek)+ dr;double ik = M22 + di + M26 * tk;/計算衛(wèi)星在開交點軌道直角坐標(biāo)系的坐標(biāo)double xk = rk * Math.Cos(uk);double yk = rk * Math.Sin(uk);/計算開交點經(jīng)度double

24、 we = 7.2921151467e-05;double jdt = M20 + (M25 - we) * tk - we * M18;/衛(wèi)星在地固坐標(biāo)系中的空間直角坐標(biāo)XKl = xk * Math.Cos(jdt) - yk * Math.Cos(ik) * Math.Sin(jdt);YKl = xk * Math.Sin(jdt) + yk * Math.Cos(ik) * Math.Cos(jdt);ZKl = yk * Math.Sin(ik);)/輸出衛(wèi)星在地固坐標(biāo)系中的空間直角坐標(biāo)for (int i = 0; i < 21; i+)(ListViewItem list

25、 = new ListViewItem();list.Text ="第"+ i.ToString() + " 分鐘"list.SubItems.Add(XKi.ToString();list.SubItems.Add(YKi.ToString();list.SubItems.Add(ZKi.ToString();listView2.Items.Add(list);)private void button2_Click(object sender, EventArgs e)(/ 讀取相對路徑string str1 =AppDomain.CurrentDom

26、ain.SetupInformation.ApplicationBase;string filename = str1 + "navigation data.txt "/讀取衛(wèi)星播送星歷文件StreamReader myreader = new StreamReader(filename, Encoding.Default);string myinfo = myreader.ReadToEnd();myreader.Close();/把衛(wèi)星播送星歷里的D改為e,以便后續(xù)計算string mystring = myinfo;string myinfor = mystring;

27、myinfor = mystring.Replace("D", "e");/把數(shù)據(jù)分開讀入到一個數(shù)組中存儲string口 split = new string " " ;string arrs = myinfor.Split(split,StringSplitOptions.RemoveEmptyEntries);/ 星歷代碼string N = new stringarrs.Length;N0= "PRN" N1 = "Yer" N2= "Mon" N3 = "d

28、ay" N4="H" N5 = "M" N6 = "sec" N7 = "a0" N8 = "a1" N9 = "a2"N10 = "IODE" N11 = "Crs" N12 = "delta-n" N13= "M0"N14= "Cuc" N15 = "e" N16 = "Cus" N17 = "sqrt(a)&qu

29、ot;N18= "t0e" N19 = "Cic" N20 = "omega0" N21 = "Cis"N22 = "i0" N23 = "Crc" N24 = "omega" N25= "omega-spot"N26 = "IDOT" N27 = "Codes on L2 channel" N28 = "GPS Week" N29 = "L2 P data flag

30、"N30 = "SV accuracy" N31 = "SV health" N32 = "TGD"N33 = "IODC Issue of Data"N34 = "Transmission time of message"/把衛(wèi)星播送星歷文件內(nèi)容輸出到listView進(jìn)行查看for (int i = 0; i < arrs.Length; i+)ListViewItem li = new ListViewItem();li.Text = (i + 1).ToString();l

31、i.SubItems.Add(Ni);li.SubItems.Add(arrsi);listView1.Items.Add(li);四.程序運行結(jié)果星歷內(nèi)容的讀取、顯示:顯示衛(wèi)星星所文件內(nèi)容序號星歷內(nèi)容導(dǎo)航數(shù)據(jù)1IEUT72Yer43Hor1qday3DL5KSe070.0S0 E293253T7313-03日-0. 5<1569&32L0&4e-11100-0000000X00 O-e+OO11工ODE0 210000000000312Crs0 53西12500口叩電何£13del ta.-si0.4B5305023777e-0&14-0, Z1499

32、2476B72e+0115Cue0 2T54a5E17571e-0515g0 12553691351-01ITCm0 903755426407-0513sqrt (a)0 5153T3口羽田麻田419tOe0 4&口80口00000以卻&20Cie0 20&7536J1565e-0&El阿0 193190315097122Cis-0. T264316O32®e-0723100 93769009471 &eMJ0廿衛(wèi)星坐標(biāo)計算結(jié)果:訐SOM里粽("in r"ii iiirfin!""! mr'iin

33、 llll,F,i-時間支坐標(biāo)建標(biāo)第咐鑰1-13113645.4&3S0051609ZCE7.2&E309ZieSZ6OZ4.39269T&第1分鑰1-13W56BT,T535£57110556,南1303165(ng 31Kx4第渤鑰1T3557口16,270341316329,1953Z2216385603.3302314第3分鑰1-13力7在乙5T 6556916H6S96.633TOOT16 煙 543.7533164第嗎鑰1-13396496.3567713IB1&5239.29E3028161302&3,475666&第5分

34、鑰1-1014335.424OCK16163EL7,46401116012771.1636663第兇朝1-14132225,721007516202471.5293033l5S8607B.75&04第俏智1-1 <1243651.321409416221241.635Z33T15756138.3209013第吩釣1-14X4134.430894716240116.579861915 陷 116.4 10132第涉智1-1447363T,333356G16259091.9QEB50Z154 兆 B69.35 比 364第 1 口分鐘 -1459216£. 8&7口

35、子曰16278151,6&923S515357457.1625691第 1L 分鐘 -14704702.3T5622216297236 47124T3152348S8.9691003第1工分鐘 -14S葭印.彩934491631&49J.6500163 117% &114355第 13分鐘 -14926799.201Q52216335751.137534914966323.8101145第 14分鐘 -15036341. EE2E51216353056.S561E9T14B303銀41總咨第 15 分鐘 -15144869,931133163796.63057614日9

36、3252.羽 6792第 16 分鐘 -15E523T6.1rB35 國 9163937S5.6&9765514555051.472E第17分鐘 T5358朋1,592.1516413146.2&85T021441s7弭 02 筆 543第歸分鐘 -15<1&4310,4211616432535.359609914ZT53T0,1200T4B第1日分鐘-15%仃級電721,16451-9L6.61T5113KL33910.02=36661第匆分鐘 -15&T20S5,4436鞏816411231.4580113991334. J7S2333堇亞俵一運行界面

37、:品三p星早歷卻嘻用弓里歷的考刷藏苗Jk汴同理奇百標(biāo)理標(biāo)1FFJ7?Yer43Mcil14d 班305JiE6MQ7see0.d8W訊眈Q區(qū)小亂31Alfl 54300210-11L 0衣Q. OODDCOXDaUDafaLLICBE0.2：OtO：CCOOCiDeKiO12Crs0.5JT!H?5nOOOn*T2L3Me(M853J5Q£8TTEL4MO-O.EL4,39e4TW7Ze+01l 與CtocL&gOJ2S3DOCG1351.-OlITCusQ 03F55®ftinTC&L8%r 七1>0.5L5OT3O398lB«HXL9tCeO.460»LCC00ODel£fiZDCicD 20&?3BI5«5*<S2LmN0.193l'K315097cKil2之Ci-U,準(zhǔn)弘31刖道:口 gUT2310D. Mli' b J JJ“T ：第.分腫以*1*1工*&町口1,&C