Файл: Брандин, В. Н. Основы экспериментальной космической баллистики.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 17.10.2024

Просмотров: 188

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Естественно, при этом увеличивается ошибка определения матрицанта. Эту ошибку можно сделать допустимой в каждой конкрет­ ной задаче выбором соответствующего числа т интервалов раз­ биения основного интервала (Y0, /]. Понятие мультипликативного интеграла может быть применено также для определения при­ ближенного решения неоднородных линеаризованных уравнений вида (2.6.17).

§ 2.9. МЕТОДЫ РАСЧЕТА МАТРИЦАНТА

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

точные аналитические;

приближенные аналитические;

численные.

2.9.1. Точные аналитические методы

Точные аналитические методы основываются на возможности получения модели движения космического объекта в замкнутой форме. Такая возможность связана с принятием предположения о кеплеровом движении центра масс объекта и движении Эйле­ ра-— Пуансо относительно центра масс. Наиболее полно точные аналитические методы применительно к кеплеровому движению объекта рассмотрены в работе [67]. Этим же вопросам посвяще­ ны работы [18, 64, 72—74]. В указанных работах уделяется вни­ мание выбору наиболее приемлемой системы координат с целью получения более простого и удобного представления матрицанта. В качестве таких систем используются орбитальные или эквато­ риальные вращающиеся и невращающиеся прямоугольные сис­ темы, а также система цилиндрических координат. Существенное влияние на вид конечных зависимостей, используемых для рас­ чета матрицанта, оказывает выбор параметров, определяющих взаимное положение исходной и конечной точек на опорной тра­ ектории объекта. В качестве таковых применяются: время, ис­ тинная и эксцентрическая аномалии, разность времен или ука­ занных аномалий для конечной и исходной точек.

Выбор формы окончательного представления матрицанта необходимо связывать с решением конкретной задачи, и здесь вряд ли уместны какие-либо общие рекомендации. Для построе­ ния матрицанта с использованием точных аналитических методов можно использовать два пути. Первый из них основывается на соответствующих преобразованиях обычных моделей кеплерова движения или движения Эйлера — Пуансо, второй — использует линеаризованные модели указанных движений. Один из сущест­ венных недостатков точных аналитических методов связан преж­ де всего с погрешностями матрицанта, возникающими из-за от­ личия возмущенного движения космического объекта от опорного

56


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

2.9.2. Приближенные аналитические методы

Приближенные аналитические методы базируются на разло­ жении в ряд Тейлора интегралов обычной [63] или линеаризован­ ной [52] модели движения космического объекта. Рассмотрим их суть на примере определения матрицанта для поступательного движения центра масс. Пусть на интервале [U, /] движение цент­ ра масс объекта описывается моделью (2.1.12)

«

^ = —

f6

x x = Vy\

 

 

 

 

 

 

v 2= — !Л0

*2

X2 = V2\

(2.9.1)

 

v* —

' П) '

£з_

X , = V :

 

 

 

/-3

з-

 

Пусть искомый вектор х параметров движения удовлетворя­ ет модели (2.9.1) и соответствует вектору начальных условий ж0, заданных для момента времени t = t0:

х т^Wx -Xl ^ v ^ V s I;

(2.9.2)

•* 0 — |] * 10* 20* 30^10^20 ^ 3 0 II

(2 .9 .3 )

Рассматривая правые части модели (2.9.1) как аналитические функции, разложим каждую компоненту Xk {k=\, 2, ..., 6) векто­ ра х интегралов системы (2.9.1) в ряд Тейлора:

dxk

d2xt

 

+ ... (k = \ , 2 ,

6),

** = -%)+ dt {к)

dt2

( к )

(2.9.4)

где точками обозначены члены, содержащие время t в более вы­

57


сокой степени. Удерживая в разложении (2.9.4) члены, содержа­ щие время t в первой степени, запишем

-----' H j g ----- fig

■*10

t\

x xx w-\-vmt\

 

 

=

^20 - 1 * 0

~ ~

t\

x 2=

л 20 + v 2(J ;

(2.9.5)

®з= ^зо

 

 

 

 

V30

 

 

 

Дифференцируя (2.9.5)

по компонентам xj0 (j= 1, 2, ...,

6),

по­

лучим матрицант, представленный в виде матрицы Якоби

 

 

 

 

 

дхх

 

дхх

 

дхх

 

 

 

 

 

дх10

 

dx2Q

 

dv3Q

 

 

Ф(*. *о)=

 

dv3

 

dv3

 

dv3

(2.9.6)

 

 

 

 

 

 

 

 

 

 

дх10

 

дх20

 

dv30

 

 

п

дх\

дхх

 

dv3

легко находятся

на

где производные

 

дх20

 

dv30

дх10

 

 

 

 

--------- -

 

 

 

изложенный подход можно при­

менить для получения в разложении (2.9.4) члена, содержащего t , и последующего уточнения элементов матрицанта (2.9.6). По­ скольку разложение (2.9.4) при ограниченном числе членов яв­ ляется приближенным, то рассмотренный метод применяют для небольших интервалов [t0, t}. Если возникает необходимость в определении матрицанта для точек, значительно удаленных от начальной, то интервал (70, t] может быть разбит на ряд корот­ ких интервалов при помощи узловых точек г|?0. Тогда на каждом из коротких интервалов можно применить разложение (2.9.4) относительно узловой точки и получить для короткого интервала матрицант Ф(^, /ф0), где t — время, принадлежащее ф-му интер­ валу. В соответствии со вторым свойством матрицанта имеем

*о)= ®(*, ^ „ )ф (^0, 4,).

(2.9.7)

где Ф (^в, t) может быть рассчитан точным аналитическим ме­ тодом или численным, если нельзя пренебречь влиянием возму­ щении. Ьолее общим (в смысле полноты учитываемых возмуще­ нии) приближенным аналитическим методом по сравнению с рассмотренным выше является метод, основанный на разложе­ нии в ряд Тейлора вектор-функции x(t), удовлетворяющей одно­ родной линеаризованной модели (2.7.1), соответствующей неод­ нородной модели вида (2.6.17). Предполагается, что векторфункция x{t) является аналитической для всех точек t, принад-

58


лежащих интервалу |Y0, f\. Тогда в окрестности произвольной точки решение уравнения (2.7.1) может быть приближенно представлено в виде т + 1 членов ряда Тейлора

тs

(2.9.8)

5= 0

где (/ф)— вектор производных порядка s в момент Ц. Последовательно дифференцируя (2.7.1) и полагая

ДО)

М)

 

(г!ф)= лгф; xls) (^)= дг^

Л (/ф )= Л ф; A {s){t^)=A\s\

получим

 

•K'l) - - А ф.дгф,

 

(2.9.9)

Д 2>_ дП)*.

 

-Л.фид£;(1)ф )-

 

(2.9.10)

ф

л\ф X ф ‘

 

 

х <а)--=А%)х* + 2 А Р Л ' )4-А*фЛфхт

и т. п.

(2.9.11)

Подставляя теперь последовательно (2.9.9) в

(2.9.10) и

(2.9.10) в (2.9.11), имеем

 

 

 

 

 

Д2) ( А ^

Aty) Хуш,

 

(2.9.12)

4 3)= ( Лф+

\1[)А ф+

ЛфЛ^1+

АI) дгф.

(2.9.13)

Подстановка (2.9.9),

(2.9.12) и

(2.9.13)

в (2.9.8)

для s = 3

дает

 

 

 

 

 

* ( 0 = [ я + ( * - Л ) л <

 

- (Лф1^ а \)-\-

 

— ( Лф )-(-2Аф1)Лф+ ЛфЛф1^ Лф)1 АГф

(2.9.14)

6

 

 

 

 

 

или з общем случае

 

 

 

 

 

S— 1

 

 

 

 

 

x{t) = ^ i Cks_lAw { t)x ^ - k-^{t)

( s = l,

2 ,..., m),

(2.9.15)

ft-0

 

 

 

 

 

где С*§_i — биномиальные коэффициенты;

Л (ft)(*) — матрицы, по­

лученные ^-кратным дифференцированием A{t).

Имея в виду, что

 

лг(/) = Ф {t, ^)jc+,

(2.9.16)