Улучшенный метод Эйлера

                   

 

Классические методы Рунге-Кутта

Методы Рунге-Кутта являются одноточечными методами, и для вычисления очередного значения функции используют только одно значение функции в предыдущей точке интегральной кривой и на каждом шаге правая часть дифференциального уравнения вычисляется несколько раз. Эти методы позволяют изменять величину шага интегрирования в любой точке и имеют достаточно большую область устойчивости.

Решение задачи Коши методами Рунге-Кутта описываются соотношениями:

                  

где r – порядок метода, h – шаг интегрирования,

αn, βrm, prn – постоянные коэффициенты.

 

 

Для примера классический метод Рунге-Кутта 4-го порядка имеет вид:

                            

                         

 

В данном разделе представлена математическая модель самолёта в виде системы дифференциальных уравнений поступательного движения центра масс и вращательного движения относительно центра масс как жесткого материального тела.

Уравнения Эйлера углового движения ЛА:

Уравнения поступательного движения ЛА в ЗСК:

                    

Уравнения вращательного движения ЛА в ССК:

 

Кинематические уравнения Пуассона

Кинематические уравнения Пуассона и условия ортогональности осей:

       

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

 

Из уравнений Пуассона находим углы Эйлера:

;                                         

;                                           

;                                           

В целях сохранения нормы матрицы и ортогональности осей на направляющие косинусы накладываются условия:

                                        

Путевая скорость:      

Угол наклона траектории (УНТ) имеет вид: Q = arcsin(Wy/W) или Q = arctg(Wy/Wx).

Географические координаты долгота и широта вычисляются из дифференциальных уравнений:

                     

где R – радиус Земли, H – высота полёта, wз – угловая скорость вращения Земли.

 


Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  



double arrow
Сейчас читают про: