GPS卫星坐标计算

1、精选优质文档-倾情为你奉上卫星定位技术与方法-根据广播星历参数计算卫星坐标作业报告指导教师： 熊 永 良 班 级： 测 绘 二 班 学生姓名： 段 海 东 学生学号： 2 0 0 8 0 7 8 3 作业日期： 2010 年 12月 08 日 目 录 根据广播星历参数计算卫星坐标一. 已知数据： 根据以下的广播星历参数计算UTC2004年1月30日8点0分00秒20分00秒，每隔一分钟的PRN7的卫星坐标。Compute the coordinate of PRN7 with interval of 1 minute.Navigation data:卫星导航文件格式：二. 计算步骤：The s

2、teps for satellite coordinates1.平均角速度 (mean angular speed)：Dn由广播星历获得, GM=3.e+142.规化时刻(normal time)：t0已知（由广播星历获得），t为GPS周秒3.平近点角(mean anomaly)： M0已知（由广播星历获得）4.偏近点角(eccentric anomaly)：迭代求解：初始值取E=M，以弧度为单位5.真近点角(true anomaly)：6.升交距角(argument of ascending node)：w近地点角距(argument of perigee)7. 轨道向径（Orbital r

3、adius）：8. 扰动改正（Perturbed correction）： 升交角距（Argument of ascending node） 轨道向径（Orbital radius） 轨道顷角（Orbital inclination） 是升交角距 (the argument of ascending node)9. 改正后升交角距、轨道向径、轨道倾角改正后升交角距（Corrected argument of ascending nod）改正后的轨道向径（Corrected orbital radius)改正后的轨道倾角（Corrected orbital inclination） 10.卫星在

4、升交点轨道直角坐标系中的坐标：如下图所示 11. 升交点经度（Longitude of ascending node）：如下图所示12. 在地固坐标系中的卫星位置（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.I

5、O;using System.Text.RegularExpressions;namespace Test public partial class Form1 : Form public Form1() InitializeComponent(); listView1.Columns.Add("序号", 40); listView1.Columns.Add("星历内容", 130); listView1.Columns.Add("导航数据", 130); listView1.GridLines = true; /显示表格线 list

6、View1.View = View.Details;/显示表格细节 listView1.HeaderStyle = ColumnHeaderStyle.Clickable;/对表头进行设置 listView2.Columns.Add("时间", 60); listView2.Columns.Add("x坐标", 150); listView2.Columns.Add("y坐标", 150); listView2.Columns.Add("z坐标", 150); listView2.GridLines = true;

7、 /显示表格线 listView2.View = View.Details;/显示表格细节 listView2.HeaderStyle = ColumnHeaderStyle.Clickable;/对表头进行设置 private void button1_Click(object sender, EventArgs e) /读取相对路径 string str1 = AppDomain.CurrentDomain.SetupInformation.ApplicationBase; string filename = str1 + "navigation data.txt "

8、/读取卫星广播星历文件 StreamReader myreader = new StreamReader(filename, Encoding.Default); string myinfo = myreader.ReadToEnd(); myreader.Close(); /把卫星广播星历里的D改为e，以便后续计算 string mystring = myinfo; string myinfor = mystring; myinfor = mystring.Replace("D", "e"); /把数据分开读入到一个数组中存储 string split

9、 = new string " " ; string arrs = myinfor.Split(split, StringSplitOptions.RemoveEmptyEntries); double M = new doublearrs.Length; /星历代码 string N = new stringarrs.Length; N0 = "PRN" N1 = "Yer" N2 = "Mon" N3 = "day" N4 = "H" N5 = "M"

10、 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 = "Ci

11、c" 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" N30 = "SV accuracy"

12、; N31 = "SV health" N32 = "TGD" N33 = "IODC Issue of Data" N34 = "Transmission time of message" /把卫星广播星历读入到数组中 for (int i = 0; i < arrs.Length; i+) ListViewItem li = new ListViewItem(); li.Text = (i + 1).ToString(); li.SubItems.Add(Ni); li.SubItems.Add(arrs

13、i); listView1.Items.Add(li); Mi = double.Parse(arrsi); double t0e=.00; double t; double XK = new double21; double YK = new double21; double ZK = new double21; for (int l = 0; l < 21;l+ ) /计算平均角速度 double GM = .418e008; double n, n0; n0 = Math.Sqrt(GM / (Math.Pow(M17, 6); n = M12 + n0; /规划时刻 t=t0e+

14、l*60; double tk = t-t0e; double Mk = M13 + n * tk; / 迭代计算平近点角的计算 double Ek, Ek1; Ek = Mk; Ek1 = 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); /反

15、正切值的象限处理 double A1 = Math.Atan(Math.Sqrt(1 - M15 * M15) * Math.Sin(Ek); double A2 = Math.Cos(Ek) - M15; 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; /计算升交角距 doubl

16、e fk = Vk + M24; /摄动改正值的计算 double du = M14 * Math.Cos(2 * fk) + M16 * Math.Sin(2 * fk); double dr = M23 * Math.Cos(2 * fk) + M11 * Math.Sin(2 * fk); double di = M19 * Math.Cos(2 * fk) + M21 * Math.Sin(2 * fk); /摄动改正 double uk = fk + du; double rk = Math.Pow(M17, 2) * (1 - M15 * Math.Cos(Ek) + dr; do

17、uble ik = M22 + di + M26 * tk; /计算卫星在升交点轨道直角坐标系的坐标 double xk = rk * Math.Cos(uk); double yk = rk * Math.Sin(uk); /计算升交点经度 double we = 7.e-05; double jdt = M20 + (M25 - we) * tk - we * M18; /卫星在地固坐标系中的空间直角坐标 XKl = xk * Math.Cos(jdt) - yk * Math.Cos(ik) * Math.Sin(jdt); YKl = xk * Math.Sin(jdt) + yk *

18、 Math.Cos(ik) * Math.Cos(jdt); ZKl = yk * Math.Sin(ik); /输出卫星在地固坐标系中的空间直角坐标 for (int i = 0; i < 21; i+) ListViewItem list = new ListViewItem(); list.Text = "第" + i.ToString() + "分钟" list.SubItems.Add(XKi.ToString(); list.SubItems.Add(YKi.ToString(); list.SubItems.Add(ZKi.ToStr

19、ing(); listView2.Items.Add(list); private void button2_Click(object sender, EventArgs e) /读取相对路径 string str1 = AppDomain.CurrentDomain.SetupInformation.ApplicationBase; string filename = str1 + "navigation data.txt " /读取卫星广播星历文件 StreamReader myreader = new StreamReader(filename, Encoding.D

20、efault); string myinfo = myreader.ReadToEnd(); myreader.Close(); /把卫星广播星历里的D改为e，以便后续计算 string mystring = myinfo; string myinfor = mystring; myinfor = mystring.Replace("D", "e"); /把数据分开读入到一个数组中存储 string split = new string " " ; string arrs = myinfor.Split(split, StringSp

21、litOptions.RemoveEmptyEntries); /星历代码 string N = new stringarrs.Length; N0 = "PRN" N1 = "Yer" N2 = "Mon" N3 = "day" N4 = "H" N5 = "M" N6 = "sec" N7 = "a0" N8 = "a1" N9 = "a2" N10 = "IODE" N1

22、1 = "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" N23 = "Crc" N24 = "omega" N25 = "omega-spot" N26 = "IDOT" N27 = "Codes on L2 channel