Интерполяция направления неустойчивости

Зависимость направления неустойчивости от координат X, Y КА образует поверхность, проекции которой представлены на рис. 36-38.

Рис. 36. Точки, для которых рассчитано направление неустойчивости.

Рис. 37. Зависимость направления неустойчивости от координаты X.

Рисунок 38. Зависимость направления неустойчивости от координаты Y.

Для использования этих данных при моделировании движения КА необходимо произвести их интерполяцию. Для этого была разработана функция Matlab. Входными параметрами данной функции являются текущие координаты X, Y КА и угол между радиус-вектором КА и осью X в плоскости эклиптики (Alpha). Выходной параметр - это направление неустойчивости в точке с координатами X, Y.

Точки, для которых рассчитано направление неустойчивости, образуют при проекции на плоскость XY 360 лучей. По параметру Alpha в функции определяются 2 луча, между которыми лежит исследуемая точка. Затем на этих лучах ищутся 3 ближайшие к ней точки, через которые проводится плоскость, после чего в уравнение данной плоскости подставляются координаты исследуемой точки. Полученное значение является интерполированным значением направления неустойчивости.