Математическая модель движения спутника по орбите

Федеральное государственное бюджетное  образовательное учреждение высшего  профессионального образования

 

Ульяновский Государственный Университет

 

 

 

           Факультет Математики и Информационных технологий

          

           Кафедра     Информационной безопасности и теории управления

                  

 

 

 

КУРСОВАЯ  РАБОТА

 

              Математическая модель движения спутника по орбите.

            -------------------------------------------------------------------------------------

                                                      (название темы)

 

 

 

Прикладная  математика и информатика  01.05.01

Математика 01.01.01

            --------------------------------------------------------------------------------------

                                   (наименование и номер специальности)

 

 

 

ПМ-о-09/1                                          Морозов И.В.

   Проект выполнил студент      ____________       _____________        ___________________

                                                   группа                        подпись, дата                                    Ф.И.О.

профессор, д.ф.-м.н.                                Андреев А.С.

   Научный руководитель    _____________________     _____________     _________________

                                                 должность                               подпись, дата                           Ф.И.О.

 

_________________

оценка

 

 

 

 

 

У Л Ь Я  Н О В С К

 

2013 г.

 

Содержание.

Введение

 

Глава 1. Кеплерова орбита спутника.

 

Глава 2. Разработка программы.

 

Заключение

 

Библиографический список

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введение.

 

С древних времен и до наших дней вид звездного неба будоражил воображение людей своей красотой и величием. Исторически одна из первых задач астрономии заключалась в разумном объяснении движении планет, случайным образом изменяющих свое положение относительно звезд. Над поиском ответа на эту задачу ученые безуспешно бились на протяжении почти двух тысячелетий вплоть до Кеплера. Законы Кеплера математически выражают поразительную простоту планетных движений, наблюдаемых гелиоцентрической системе отсчета. Это был поистине фантастический прорыв в понимании Природы. Но и поныне движения небесных тел – малых и больших планет Солнечной системы, их спутников, комет, астероидов, а в наше время – также рукотворных космических кораблей и искусственных спутников – дают наиболее впечатляющие опытные подтверждения законов. Современное общество не представляет свое существование без искусственных спутников, потому что в настоящее время это рабочие элементы глобальных связных телекоммуникационных систем, систем прогнозирования погоды, навигационных систем и организуется новая система обнаружения космических объектов, представляющих угрозу земле от столкновения с ними. Также ИСЗ используют для военных целей - разведки, слежения за состоянием военных сил противника. Со спутников ведется разведка полезных ископаемых. Ну а интернет не смог бы иметь такие масштабы развития без ИСЗ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 Глава 1

 

  Кеплерова орбита  спутника.

Если спутник выведен  на заданное расстояние r от центра (масс) Земли и его двигатели отключены, то в данной точке пространства он приобрел известную нас скорость v, а также направление движения – угол x между направлением движения спутника и радиус-вектором r. Этих данных достаточно, чтобы определить кеплеров эллипс, по которому спутник будет двигаться вокруг Земли. Здесь мы исходим из положения, что вся масса Земли сосредоточена в ее центре и на спутник не действуют прочие возмущающие массы (Луна, Солнце, планеты).

Используя исходные данные r, v, и x, определим постоянные величины эллипса рис.1.

Большая полуось орбиты

                                              (1)

( G - гравитационная постоянная, M - масса Земли, GM – геоцентрическая гравитационная постоянная),

эксцентриситет

      (2)

 

Рис. Кеплеров эллипс

О – центр; , - фокусы; a, b – большая и малая полуоси; - фокальный параметр; - расстояние от центра до фокуса; П, А,- точки перицентра и апоцентра; А-П – линия апсид; GM – центрическая постоянная; r – радиус – вектор тела m; - расстояние от m до второго фокуса; v – скорость тела m; - средняя скорость; х – угол между радиусом- вектором и касательной; , , - эксцентрическая, средняя и истинная аномалии.

 

 

     2

Малая полуось орбиты :

                           (3)

Расстояние от центра Земли до центра эллипса:

                                 (4)  

Фокальный параметр эллипса:

                          (5)

Средняя скорость спутника (если бы спутник двигался по кругу  с центром в центре эллипса):

                              (6)

Расстояние от центра Земли до ближайшей точки эллипса  – перигея:

                              (7)

Расстояние от центра Земли до наиболее  удаленной  точки эллипса – апогея:

                               (8)

Скорости спутника в  перигее:

                (9)

Скорость спутника в  апогее:

                (10)

Период обращения спутника по эллипсу:

          (11)

Площадь эллипса:

                       (12)

Геоцентрическая гравитационная постоянная:

                  (13)

( - круговая скорость тел на расстоянии r от центра масс Земли).

 

Таким образом на основе минимальных данных v, r, x определены постоянные параметры эллипса, являющегося траекторией невозмущенного движения спутника.

Далее находим переменные величины, характеризующие спутник  на орбите в момент отключения двигателей: истинная аномалия (угол между направлением на перигей и радиус-вектором r) рис.

      (14)

Эксцентрическая аномалия :

               (15)

 

 

3

Средняя аномалия :

                            (16)

Равенство (16) – это  трансцендентное уравнение Кеплера. Средняя аномалия представляет собой угол поворота большой полуоси за некоторый промежуток времени в предложении, что орбита является круговой. А в связи с тем, что вычисления по (16) должны производиться в радианной мере, то следовательно, это и дуга круга, которую

 

описывает орбитальное тело за данное время. Именно благодаря средней  аномалии , всегда пропорциональной времени, находится промежуток времени, необходимый телу при движении от линии апсид до заданной точки эллипса

                       (17)

Однако надо сказать, что при вычислении орбит планет вначале, исходя из заданного времени, определяется средняя аномалия по равенству (17), а затем уже методом последовательных приближений находится эксцентрическая аномалия по уравнению Кеплера (16). Знание вспомогательной величины дает возможность находить другие необходимые переменные величины на заданный момент времени и таким образом уравнение (16) является важнейшим в решении задач движения тел по орбитам.

По второму закону Кеплера определим  площадь, описываемую радиусом-вектором за известный промежуток времени 

                                  (18)

Определим и расстояние от спутника до второго фокуса эллипса

                                      (19)

Проверим точность найденных  углов  и по одной из формул Кеплера.

                                 (20)

или по следующим равенствам

        (21)

После того, как установлены  основные параметры спутника на момент выключения двигателей, необходимо определить его положение на орбите через какой-либо другой промежуток времени, прошедший от исходного момента. Для этого к известному промежутку (за который спутника проходит по орбите от перигея до точки отключения двигателей) прибавим интересующее нас время :

                                 (22)

По суммарному времени  находиться средняя аномалия по (17), а по уравнению Кеплера (16) методом последовательной итерации определяется эксцентрическая аномалия . Именно благодаря трансцендентному уравнению Кеплера (16) можно получать характеристики тела в любой из моментов его движения по эллипсу.

Определим параметры спутника на заданный момент времени:

Угол между радиусом-вектором и  направлением движения спутника

                     (23)

 

 

 

4

угол между радиусом-вектором и направлением на перигей – истинная аномалия

                              (24)

расстояние от центра Земли до спутника

                          (25)

Вычислим скорость спутника в заданный момент времени по одной из формул Кеплера (позже названной интегралом энергии)

                                     (26)

Наконец по второму закону Кеплера (18) вычислим площадь, описываемую  радиус-вектором за суммарный промежуток времени (22).

Однако, если мы намерены вначале задать не промежуток времени , а угол между радиусом-вектором и направлением на перигей и определить параметры спутника, когда он находиться под заданным углом, то вначале по углу определим расстояние от центра масс Земли до спутника

                                        (27)

Также по углу находиться и угол между направлением движения спутника в данной точке пространства и радиусом-вектором

                            (28)

затем по интегралу энергии (26) находиться скорость спутника в  данной точке.

Иногда можно исключить использование  вспомогательных углов  и из наших формул, хотя в действительности новые формулы выводятся благодаря уравнению Кеплера (9.16).

 

При очень малой скорости по эллипсу, пересекающему поверхность Земли, спутник упадет на Землю. Если двигатели  спутника остановлены при скорости, превышающей , то орбита космического аппарата станет параболической и он преодолеет тяготение Земли и перейдет в сферу гравитационного действия Солнца.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

Глава 2.

 

Разработка программы.

Код, находящийся в  главной форме Form1.cs

using System;

using System.Collections.Generic;

using System.ComponentModel;

using System.Data;

using System.Drawing;

using System.Drawing.Drawing2D;

using System.Linq;

using System.Text;

using System.Threading.Tasks;

using System.Windows.Forms;

 

namespace Kepler

{

    public partial class Form1 : Form

    {

        public Form1()

        {

            InitializeComponent();

            graphic = this.CreateGraphics();

           // _DrawGraphics();

        }

        private Graphics graphic;

        private Geomtries.Ellipse Ellipse = new Geomtries.Ellipse(0,0,1,1);

        private Geomtries.Ellipse Earth = new Geomtries.Ellipse(10,10,1,1);

        Geomtries.Ellipse Satelite = new Geomtries.Ellipse(10, 10, -10, 10);

        /// <summary>

        /// Угол наклона радиус вектора

        /// </summary>

        private int arcAngleOffset = 1;       

        /// <summary>

        /// Глобальный таймер - на каждом отсчёте фиксируется точка и линии от фокусов до неё

        /// </summary>

        private Timer GlobalTimer = new Timer();

        private double m_speed;

        private double speed

        {

            get

            {

              //  double exs = Math.Sqrt(1 - (((Ellipse.Height / 2) * (Ellipse.Height / 2)) / ((Ellipse.Width / 2) * (Ellipse.Width / 2))));

                double r1 = Distance(Satelite.X,Earth.X,Satelite.Y,Earth.Y);

               // double r2 = (Ellipse.Width / 2) / (1 - exs*Math.Cos(arcAngleOffset));

                m_speed = Math.Sqrt(6.67384 * 5.9736 * Math.Pow(10, 13)*

                    ((2/r1) - (1/(Ellipse.Width/2))));

                return m_speed;

            }

            set { m_speed = value; }

        }

        private double Distance(double x1, double x2, double y1, double y2)

        {

           return Math.Sqrt(Math.Abs((x2 - x1) * (x2 - x1) + (y2 - y1) * (y2 - y1)));

        }

        /// <summary>

        /// Нарисовать при прорисовки формы

        /// </summary>

        /// <param name="sender"></param>

        /// <param name="e"></param>

        public void Form1_Paint(object sender, System.Windows.Forms.PaintEventArgs e)

                                         6

{

            _DrawGraphics();                     

        }

 

        /// <summary>

        /// Нарисовать при изменении размера формы

        /// </summary>

        /// <param name="sender"></param>

        /// <param name="e"></param>

        private void Form1_SizeChanged(object sender, EventArgs e)

        {

            _DrawGraphics();

        }

 

        /// <summary>

        /// Основная функция прорисовка осей, земли и орбиты

        /// </summary>

        private void _DrawGraphics()

        {

            graphic = this.CreateGraphics();

            // Задаём новые значение ширины и высоты

            Ellipse.Height = this.Size.Height-40; Ellipse.Width = this.Size.Width-20;

            Earth.Height = 50; Earth.Width = 50;

            Earth.X = (float)Ellipse.Focus - 20;

            Earth.Y = (this.Size.Height / 2) - 20 - (int)Earth.Height / 2;

            #region Работа с label'ами

            lEllipseWidth.Text = "Ellipse width=" + Ellipse.Width;

            lEllipseHeight.Text = "Ellipse height=" + Ellipse.Height;

            lFocus.Text = "Focus=" + Ellipse.Focus;

            lFormHeight.Text = "Form Height=" + this.Size.Height;

            lFormWidth.Text = "Form Width=" + this.Size.Width;

            #endregion

           

           // graphic.Clear(System.Drawing.Color.WhiteSmoke);

            // Ось OX

            graphic.DrawLine(new Pen(Color.White, 1), 1, (this.Size.Height / 2) - 20, this.Size.Width, (this.Size.Height / 2) - 20);

            // Эллипс

            graphic.DrawEllipse(new Pen(Color.White, 1), Ellipse.X, Ellipse.Y, (float)Ellipse.Width, (float)Ellipse.Height);

            // Земля

            graphic.DrawEllipse(new Pen(Color.Aqua, 2), Earth.X, Earth.Y, (float)Earth.Width, (float)Earth.Height);

            graphic.Save();

        }

 

        private void bStart_Click(object sender, EventArgs e)

        {           

            Satelite.Y = (this.Size.Height / 2) - 20 - (float)Satelite.Height / 2;