卫星轨道计算

【实例截图】

【核心代码】

namespace 卫星坐标计算
{
    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
        }

        private Control FindControl(string controlName)
        {
            return this.Controls[controlName];
        }


        double[] M = new double[36];
        private void button1_Click(object sender, EventArgs e)
        {
            openFileDialog1.Title = "Selete a File";
            openFileDialog1.FileName = "";
            openFileDialog1.Filter = "Text File(*.txt)|*.txt";
            if (openFileDialog1.ShowDialog() == DialogResult.OK)
            {
                StreamReader read = new StreamReader(openFileDialog1.FileName);
                string infor1 = read.ReadToEnd();
                read.Close();

                string mystring =infor1;
                string infor = mystring;
                infor = mystring.Replace("D", "e");

                //把数据分开读入到一个数组中存储
                string[] split = new string[] { " " };
                string[] arrs = infor.Split(split, StringSplitOptions.RemoveEmptyEntries);

                for (int i = 0; i < arrs.Length; i )
                {
                    M[i] = double.Parse(arrs[i]);
                    richTextBox1.Text = arrs[i] " ";
                }
                textBox1.Text = arrs[1].ToString(); textBox2.Text = arrs[2].ToString(); textBox3.Text = arrs[3].ToString();
                textBox4.Text = arrs[4].ToString(); textBox5.Text = arrs[5].ToString(); textBox6.Text = arrs[6].ToString();
                textBox7.Text = arrs[7].ToString(); textBox8.Text = arrs[8].ToString(); textBox9.Text = arrs[9].ToString();
                textBox10.Text = arrs[10].ToString(); textBox11.Text = arrs[11].ToString(); textBox12.Text = arrs[12].ToString();
                textBox13.Text = arrs[13].ToString(); textBox14.Text = arrs[14].ToString(); textBox15.Text = arrs[15].ToString();
                textBox16.Text = arrs[16].ToString(); textBox17.Text = arrs[17].ToString(); textBox18.Text = arrs[18].ToString();
                textBox19.Text = arrs[19].ToString(); textBox20.Text = arrs[20].ToString(); textBox21.Text = arrs[21].ToString();
                textBox22.Text = arrs[22].ToString(); textBox23.Text = arrs[23].ToString(); textBox24.Text = arrs[24].ToString();
                textBox25.Text = arrs[25].ToString(); textBox26.Text = arrs[26].ToString(); textBox27.Text = arrs[27].ToString();
                textBox28.Text = arrs[28].ToString(); textBox29.Text = arrs[29].ToString(); textBox30.Text = arrs[30].ToString();
                textBox31.Text = arrs[31].ToString(); textBox32.Text = arrs[32].ToString();
                textBox33.Text = arrs[33].ToString(); textBox34.Text = arrs[34].ToString();


            }


        }


      
        private void button2_Click(object sender, EventArgs e)
        {
            richTextBox1.Text=" ";


      //清空放在groupBox里的所有文本框里的数据  
            foreach (Control Ctrol in groupBox1.Controls)
            {
                if (Ctrol is TextBox)
                {
                    Ctrol.Text = "";
                }
            }
            
           

        }

        private void button3_Click(object sender, EventArgs e)
        {
            double[] Xk = new double[100];
            double[] Yk = new double[100];
            double[] Zk = new double[100];

            double f = double.Parse(textBox35.Text);
            if (f > 100)
                MessageBox.Show("坐标数超限，请重新输入");

            else
            {
               
                double n, n0;
                double GM = 3986004.418e008, a;
                a = Math.Pow(M[17], 2);
                n0 = Math.Sqrt(GM / (a * a * a));
                n = n0 M[12];

                for (int j = 0; j < f; j )
                {
                    //规划时刻
                    double t2, t, t1,tk;
                    t2 = M[18] j * 30;
                    t1 = M[7] M[8] * (t2 - M[18]) M[9] * Math.Pow((t2 - M[18]), 2);
                    t = t2 - t1;
                    tk = t - M[18];

                    //平近点角
                    double Mk;
                    Mk = M[13] n * tk;

                    //偏近点角

                    double Ek, Ek1;
                    Ek = Mk;

                    Ek1 = Ek M[15] * Math.Sin(Ek);
                    do
                    {

                        Ek = Ek1;
                        Ek1 = Mk M[15] * Math.Sin(Ek);
                    }
                    while (Math.Abs(Ek1 - Ek) > 1.0e-12);

                    //真近点角
                    double Vk;
                    Vk = Math.Atan(((Math.Sqrt(1 - M[15] * M[15])) * Math.Sin(Ek)) / (Math.Cos(Ek) - M[15]));


                    //升交角距
                    double ok;
                   ok=Vk M[24];

                    //摄动改正
                   double ou, or, oi;
                   ou = M[14] * Math.Cos(2 * ok) M[16] * Math.Sin(2 * ok);
                   or = M[23] * Math.Cos(2 * ok) M[11] * Math.Sin(2 * ok);
                   oi = M[19] * Math.Cos(2 * ok) M[21] * Math.Sin(2 * ok);
                  

                    //改正后的升交角距
                   double uk, rk, ik;
                   uk = ok ou;
                   rk = a * (1 - M[15] * Math.Cos(Ek)) or;
                   ik = M[22] oi M[26] * tk;


                   //卫星在升交点轨道直角坐标系的坐标//
                   double xk, yk;
                   xk = rk * Math.Cos(uk);
                   yk = rk * Math.Sin(uk);
                   //升交点经度//
                   double xt, we = 7.2921151467e-005;
                   xt = M[20] (M[25] - we) * (t - M[18]) - we * M[18];

                   //卫星在地固坐标系中的空间直角坐标

                   Xk[j] = xk * Math.Cos(xt) - yk * Math.Cos(ik) * Math.Sin(xt);
                   Yk[j] = xk * Math.Sin(xt) yk * Math.Cos(ik) * Math.Cos(xt);
                   Zk[j] = yk * Math.Sin(ik);

                   ListViewItem li = new ListViewItem();

                    li.Text = t2.ToString();
                    li.SubItems.Add(Xk[j].ToString("f3"));
                    li.SubItems.Add(Yk[j].ToString("f3"));
                    li.SubItems.Add(Zk[j].ToString("f3"));

                    listView1.Items.Add(li);







                }


            }

        }

        private void button4_Click(object sender, EventArgs e)
        {
            string filePath = @"卫星数据.txt";
            
            StreamWriter SW = File.AppendText(filePath);
            SW.Write(" 历元/s" "                 X/m" "                    Y/m" "                     Z/m" "\r\n");

            foreach (ListViewItem LVI in this.listView1.Items)
            {

                SW.Write(LVI.SubItems[0].Text "             ");
                SW.Write(LVI.SubItems[1].Text "             ");
                SW.Write(LVI.SubItems[2].Text "             ");
                SW.WriteLine(LVI.SubItems[3].Text);

            }
            MessageBox.Show("输出成功，请在|卫星坐标计算|bin |Debug路径下查找txt文件");
            SW.Close();


        }

        private void groupBox1_Enter(object sender, EventArgs e)
        {

        }
    }
}