卫星坐标计算.doc

08测绘一班（20080754）GPS卫星坐标计算 班级：08测绘一班姓名：浦绍佼学号：20080754目录实验目的：3卫星坐标计算步骤：3具体过程：4四：运行与结果11五，心得体会：12实验目的： 根据导航文件求出卫星坐标。卫星坐标计算步骤：一:计算平均角速度： n =n0+n ;n0=GMa3；n：由导航文件给出二：规划时刻：tk=t-t0e, t0e为参考历元三：平近点角：Mk=M0+ntk;四：偏近点角：Ek=Mk+esinEk;(此处进行迭代运算)五：真近点角：Vk=tan-1sinEk1-e2cosEk-e;六：升交点角距：k=Vk+;七：摄动改正：顾及，I,n的摄动变化以及正弦改正模型的振幅项Cc，Cs;则 升交点角距： u=Cuccos2k+Cussin2k; 轨道向径：r=Crccos2k+Crssin2k; 轨道倾角：i=Ciccos2k+Cissin2k;式中：k为参考时刻的升交角距；八：改正后的升交角距：uk=k+u; 改正后的轨道向径：rk=a1-ecosEk+r; 改正后的轨道倾角：ik=i0+i+(IDOT)tk;九：卫星在升交点轨道直角坐标系的坐标： xk=rkcosuk;yk=rksinuk:十：升交点经度：t=0+(-e)t-t0e-et0e; e: 7.2921151467*1016 rad/s; :升交点赤经变化率； 0：GPS周开始时刻的升交点经度；十一：卫星在地固坐标系的空间直角坐标为： XkYkZk =R(-k)(R(-ik) xkykzkR(-k), (R(-ik)为旋转矩阵，将其代入展开后得： Xk=xkcosk-ykcosiksink; Yk=xksink+ykcosikcosk;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.Diagnostics;namespace 20080754 public partial class Form1 : Form private string stFilePath = string.Empty; public Form1() InitializeComponent(); listView1.Columns.Add(序号, 40, HorizontalAlignment.Center); listView1.Columns.Add(星历参数, 80, HorizontalAlignment.Center); listView1.Columns.Add(参数值, 130, HorizontalAlignment.Center); listView1.GridLines = true; listView1.View = View.Details; listView1.HeaderStyle = ColumnHeaderStyle.Clickable; listView1.FullRowSelect = true; listView2.Columns.Add(时间, 60, HorizontalAlignment.Center); listView2.Columns.Add(历元, 55, HorizontalAlignment.Center); listView2.Columns.Add(X, 115, HorizontalAlignment.Center); listView2.Columns.Add(Y, 115, HorizontalAlignment.Center); listView2.Columns.Add(Z, 115, HorizontalAlignment.Center); listView2.GridLines = true; listView2.View = View.Details; listView2.HeaderStyle = ColumnHeaderStyle.Clickable; listView2.FullRowSelect = true; private void listView1_SelectedIndexChanged(object sender, EventArgs e) public void drsj(object sender, EventArgs e) listView1.Items.Clear(); /读取相对路径 string str1 = AppDomain.CurrentDomain.SetupInformation.ApplicationBase; string filename = str1 + GPS卫星导航文件.txt ; /读取卫星广播星历文件 StreamReader readnavigation = new StreamReader(filename, Encoding.Default); string myinfo = readnavigation.ReadToEnd(); readnavigation.Close(); /把卫星广播星历里的D改为e，以便后续计算 string myinfor = myinfo; myinfor = myinfor.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 = 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 = 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; 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(arrsi); listView1.Items.Add(li); Mi = double.Parse(arrsi); private void 载入导航文件_Click(object sender, EventArgs e) drsj(sender, e); private void 计算_Click(object sender, EventArgs e) listView1.Items.Clear(); /读取相对路径 string str1 = AppDomain.CurrentDomain.SetupInformation.ApplicationBase; string filename = str1 + GPS卫星导航文件.txt ; /读取卫星广播星历文件 StreamReader readnavigation = new StreamReader(filename, Encoding.Default); string myinfo = readnavigation.ReadToEnd(); readnavigation.Close(); /把卫星广播星历里的D改为e，以便后续计算 string myinfor = myinfo; myinfor = myinfor.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 = 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 = 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; 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(arrsi); listView1.Items.Add(li); Mi = double.Parse(arrsi); double n0, n, tk, t, Mk, Ek, Ek0, Vk, faik, ksu, ksr, ksi, uk, Rk, Ik, xk, yk, lmt; double X, Y, Z; X = new double21; Y = new double21; Z = new double21; t = 60; for (int i = 0; i 1e-15); Vk = Math.Atan(Math.Sqrt(1 - Math.Pow(M15, 2) * Math.Sin(Ek) / (Math.Cos(Ek) - M15); double A1 = (Math.Sqrt(1 - M15 * M15) * Math.Sin(Ek); double A2 = Math.Cos(Ek) - M15; if (Vk 0) if (A1 0) Vk += 2 * Math.PI; if (A1 0 & A2 0) Vk += Math.PI; else if (A1 = 0 & A2 = 0) Vk = Vk + Math.PI; faik = Vk + M24; ksu = M14 * Math.Cos(2 * faik) + M16 * Math.Sin(2 * faik); ksr = M23 * Math.Cos(2 * faik) + M11 * Math.Sin(2 * faik); ksi = M19 * Math.Cos(2 * faik) + M21 * Math.Sin(2 * faik); uk = faik + ksu; Rk = Math.Pow(M17, 2) * (1 - M15 * Math.Cos(Ek) + ksr; Ik = M22 + ksi + M26 * tk; xk = Rk * Math.Cos(uk); yk = Rk * Math.Sin(uk); lmt = M20 + (M25 - 7.2921151467e-05) * tk - 7.2921151467e-05 * M18; Xi = xk * Math.Cos(lmt) - yk * Math.Cos(Ik) * Math.Sin(lmt); Yi = xk * Math.Sin(lmt) + yk * Math.Cos(Ik) * Math.Cos(lmt); Zi = yk * Math.Sin(Ik); ListViewItem li = new ListViewItem(); li.Text = 第 + i.ToString() + 分钟; li.SubItems.Add(M18 + i * t).ToString(); li.SubItems.Add(Xi.ToString(); li.SubItems.Add(Yi.ToString(); li.SubItems.Add(Zi.ToString(); listView2.Items.Add(li); private void button3_Click(object sender, EventArgs e) stFilePath = Application.StartupPath.Trim() + 卫星坐标计算结果 + .txt; StreamWriter swStream; if (File.Exists(stFilePath) swStream = new StreamWriter(stFilePath); else swStream = File.CreateText(stFilePath); swStream.Write( 时间 ); swStream.Write( ); swStream.Write( 历元 ); swStream.Write( ); swStream.Write( X ); swStream.Write( ); swStream.Write( Y ); s