《GPS卫星坐标计算》由会员分享，可在线阅读，更多相关《GPS卫星坐标计算（17页珍藏版）》请在金锄头文库上搜索。
1、0卫星定位技术与方法-根据广播星历参数计算卫星坐标作业报告指导教师： 熊 永 良 班 级： 测 绘 二 班 学生姓名： 段 海 东 学生学号： 2 0 0 8 0 7 8 3 作业日期： 2010 年 12 月 08 日 1目 录一. 已知数据 .2二. 计算步骤 .21.平均角速度 (mean angular speed) .22.规化时刻(normal time) .33.平近点角(mean anomaly) .34.偏近点角(eccentric anomaly) .35.真近点角(true anomaly) .36.升交距角(argument of ascending node) .37
2、. 轨道向径（Orbital radius） .38. 扰动改正（Perturbed correction） .410.卫星在升交点轨道直角坐标系中的坐标 .411. 升交点经度（Longitude of ascending node） .5三. 源程序 .5四 程序运行结果 .14七.作业体会 .152根据广播星历参数计算卫星坐标一. 已知数据： 根据以下的广播星历参数计算 UTC2004 年 1 月 30 日 8点 0 分 00 秒20 分 00 秒，每隔一分钟的 PRN7 的卫星坐标。Compute the coordinate of PRN7 with interval of 1 mi
3、nute.Navigation data:卫星导航文件格式：二. 计算步骤：The steps for satellite coordinates1.平均角速度 (mean angular speed)： nn0 30aGMn 由广播星历获得 , GM=3.986005e+1432.规化时刻(normal time)： 0ttkt0 已知（由广播星历获得） ，t 为 GPS 周秒3.平近点角(mean anomaly)： kk tnM0M0 已知（由广播星历获得）4.偏近点角(eccentric anomaly)： kkk EeMEsin迭代求解：初始值取 E=M，以弧度为单位5.真近点角(t
4、rue anomaly)： eEVkkk cosin1artn26.升交距角(argument of ascending node)： 近地点角距(argument of perigee)7. 轨道向径（Orbital radius）： )cos1( kk Eeark048. 扰动改正（Perturbed correction）： 升交角距（Argument of ascending node） 00u2sin2cosuSuCC 轨道向径（Orbital radius） 00r 2sin2cosrSrC 轨道顷角（Orbital inclination） 00i 2sin2cos iSiCC0
5、是升交角距 (the argument of ascending node)9. 改正后升交角距、轨道向径、轨道倾角改正后升交角距（Corrected argument of ascending nod）改正后的轨道向径（Corrected orbital radius)改正后的轨道倾角（Corrected orbital inclination）10.卫星在升交点轨道直角坐标系中的坐标：如下图所示rkk)cosEe1(aru0kki0kt)IDOT(ikksinurycox511. 升交点经度（Longitude of ascending node）：如下图所示12. 在地固坐标系中的卫星位
6、置（Expressed in spheric coordinate system）三. 源程序：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.RegularExpressions;s/rad10295467.t)(5-e eke0t kkkiyZixYsncos
7、ncok6namespace Testpublic partial class Form1 : Formpublic Form1()InitializeComponent();listView1.Columns.Add(序号, 40);listView1.Columns.Add(星历内容, 130);listView1.Columns.Add(导航数据, 130);listView1.GridLines = true; /显示表格线listView1.View = View.Details;/显示表格细节listView1.HeaderStyle = ColumnHeaderStyle.Cli
8、ckable;/对表头进行设置listView2.Columns.Add(时间, 60);listView2.Columns.Add(x坐标, 150);listView2.Columns.Add(y坐标, 150);listView2.Columns.Add(z坐标, 150);listView2.GridLines = true; /显示表格线listView2.View = View.Details;/显示表格细节listView2.HeaderStyle = ColumnHeaderStyle.Clickable;/对表头进行设置private void button1_Click(o
9、bject sender, EventArgs e)/读取相对路径string str1 = AppDomain.CurrentDomain.SetupInformation.ApplicationBase;7string filename = str1 + navigation data.txt ;/读取卫星广播星历文件StreamReader myreader = new StreamReader(filename, Encoding.Default);string myinfo = myreader.ReadToEnd();myreader.Close();/把卫星广播星历里的D改为e，
10、以便后续计算string mystring = myinfo;string myinfor = mystring;myinfor = mystring.Replace(D, e);/把数据分开读入到一个数组中存储string split = new string ;string arrs = myinfor.Split(split, StringSplitOptions.RemoveEmptyEntries);double M = new doublearrs.Length;/星历代码string N = new stringarrs.Length;N0 = PRN; N1 = Yer; N2
11、 = Mon; N3 = day; 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);N18 = t0e; N19 = Cic; N20 = omega0; N21 = 8Cis;N22 = i0; N23 = Crc; N24 = omega; N25 = omega-spot;N26 = IDOT; N27 = Codes on L2 channel; N28 = GPS Week; N29 = L2 P data flag;N30 = SV accuracy; N31 = SV health; N32 = TGD; N33 = IODC Issue of Data;N34
