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

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

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

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

Добавлен: 17.10.2024

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

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

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

или

РАФ

г>

 

— Г -

smcReK — \ т \ и е

 

п М ) = -

О (г)

(3.9.5)

 

 

где G(r) — сила притяжения Земли.

Если ньютонометр является интегрирующим или дважды ин­ тегрирующим, то для таких приборов измеряемые функции име­

ют вид

 

t

 

W (*)= j1n{t)dt\

(3.9.6)

о

 

t t

 

S{t)= Cj n[t)dtdt.

(3.9.7)

В заключение заметим, что уравнение (3.4.10) или эквива­ лентное ему уравнение (3.4.3) называют основным уравнением инерциальной навигации. Они используются как исходные при определении координат и скорости объекта по показаниям нью­ тонометров и гироскопических приборов.

§3.10. ЛИНЕАРИЗОВАННЫЕ ИЗМЕРЯЕМЫЕ ФУНКЦИИ

3.10.1.Общие зависимости, используемые для линеаризации измеряемых функций

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

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

клинеаризации измеряемых функций. Не нарушая общности, приведем здесь выводы, относящиеся к линеаризации измеряе­ мых функций первого рода. Пусть измеряемая функция

yi(t) = u,(x, t)

(3.10.1)

соответствует опорному движению, а функция

yi{t)= Ui{x, t)

(3.10.2)

возмущенному движению космического объекта. Причем

х = х - \ - а х ,

(3.10.3)

где х, х — векторы параметров возмущенного и опорного дви­ жений, заданные в опорной системе координат.

103


Отклонение Ayi(t) I-й измеряемой функции относительно зна­ чения этой функции для опорного движения

 

Д (/) = «,(*,

t) — ut(x, t)

 

(3.10.4)

или

 

 

 

 

 

 

kyi{t) = ui(x J r кх, Ь)— щ { х , /).

'

(3.10.5)

Считая измеряемую функцию непрерывной и дифференцируе­

мой,

разложим щ (х + Ддс, t) в

ряд

Тейлора

около

значения

ui(x,t):

 

 

 

 

 

 

П '

_

 

 

 

Uiix+ ь х , t) = ut ( x , 0 +

 

-

+

(3.10.6)

где

1 ’■■— частные производные

/-Й измеряемой

функции

 

dxk .

 

 

 

 

по текущим параметрам Xk (&=1, 2, ..., п), вычисленные для опор­

ного движения; ^

— остаточный член разложения.

Принимая отклонения возмущенного движения от опорного

малыми и отбрасывая в разложении

(3.10.6)

остаточный член

Rup получим для малых отклонений измеряемых функций

 

 

 

П

_

дУ ^ и ^ х +

А х , t ) - Ul(х ,

t) =

^

ьх* (ЗЛ0-7)

или

 

 

 

 

 

Д0/(0 = 2

£,*(')***.

(3.10.8)

 

k=\

 

 

 

где

 

 

 

 

 

^ ( 0 = - - ^

(3.10.9)

 

 

дхъ

 

 

Выражение (3.10.8) можно записать в векторном виде

 

=

 

 

(3.10.10)

где!;! (t ) — транспонированный векторпервых частных производ­

ных'

 

 

;И*)= 1М*)

«1Х„;

(3.10.11)

104


Ах — вектор отклонений параметров движения

 

Ах{

Дл: =

А х 2

(3.10.12)

 

Ахп

Для каждого момента времени t вектор Ц (t) можно рас­ сматривать как градиент измеряемой функции в пространстве текущего вектора x(t). Этот вектор назовем текущим градиен­ том измеряемой функции в отличие от начального градиента

4»1(0> определяемого в пространстве начальных значений x(t0) вектора x(t). Распространим теперь полученный резуль­

тат на

малые отклонения измеряемой

вектор-функции y(t) —

= и (х, t). Очевидно,

 

 

 

 

 

 

 

§1(0 АХ

 

55(0

‘1Д-*||лх1 —

§5(0 а х

 

ДУ(0 =

(3.10.13),

 

5» (0 тх п

 

 

§Ш(0 АХ

или

 

 

 

 

 

дз>(/) = Е(*)длг,

(3.10.14)

 

§п (0

§12 (0 • • • §1,(0

где

Н (о = §21 (0

§22 (0 • •

§2,(0

 

§„l(0

§»2(0..

§»,(0

прямоугольная матрица (т Х п ) частных производных измеряе­ мых функций по текущим параметрам движения. Матрицу 5 (t) назовем текущей градиентной матрицей. Полезно отметить одно из свойств матрицы Е (0 , состоящее в том, что для прямых (не­ посредственных) измерений

 

1

0

0 .

.0

гп

0

1

0 .

.0

II а X 3 II

 

 

 

 

0

0

0 .

.1

(3.10.15)

При решении задач определения движения космических объек­ тов, сводящихся к отысканию по экспериментальным данным на­ чальных условий движения x(tQ) = x 0, разложение измеряемой функции yi(t) —щ(х, t) в ряд Тейлора проводят около значения

щ(хо, t), соответствующего начальным условиям *о опорного

105


движения. Для этого случая выражение для малых

отклонений

измеряемой функции запишется в виде

 

Д#/(0 = Фг(0 Д-*о.

(3.10.16)

где ф/ (/)—'Транспонированный вектор первых частных произ­ водных l-й измеряемой функции по начальным условиям движе­ ния или начальный градиент измеряемой функции:

4»1(0= IIфц (*)Ф/,(0.- • - .Ф/д WllixH

, tl)\

t> Ь О

II н о

I оМ

_

(3.10.17)

(3.10.18)

(3.10.19)

Хо,Хо — векторы начальных условий возмущенного и опорного движении.

Малые отклонения измеряемой вектор-функции имеют вид

ДУ(0=^(<)Д*о. (3.10.20)

где

Фи(0

 

 

Ф21W

< Ы Л -Л 2Л 0

(3.10.21)

ФтЛ*)

'Ы Л - Л т Л О

 

прямоугольная матрица (т Х п ) частных производных измеряе­ мых функций по начальным условиям движения. Матрицу ф(1) назовем начальной градиентной матрицей. Начальная градиент­ ная матрица может быть представлена в виде произведения те­ кущей градиентной матрицы Н(^) и матрицанта Ф(^, to), вычис­ ленных для опорного движения, заданного начальными условия­ ми x{to) —Xq. В самом деле, поскольку

Длг= Ф(^, 70) длг0,

из выражения (3.10.14) следует

ду = Е (^) Ф (if, t0) b x 0.

Поэтому

= tQ).

Для t = t0

О

T [II 0

(3.10.22)

(3.10.23)

(3.10.24)

(3.10.25)

106


3.10.2. Пример определения текущей и начальной градиентных матриц измеряемой вектор-функции положения ya{t) = un(x, t)

Пусть из базисной точки / (см. рис. 3.6.3) измеряются на­ клонная дальность р(^) и три направляющих косинуса cosQj (/=1, 2, 3). Найдем для указанного состава измерений элементы текущей градиентной матрицы. Воспользуемся для этого резуль­ татом записи выражений для малых отклонений измеряемых функций. Из равенства (3.6.2) варьированием получаем

 

др(*) =

дг(0.

 

 

(3.10.26)

На основании (3.6.3)

имеем

 

 

 

 

 

Д р№ = Д р(0ер +

Р(0Л£Р-

 

(3.10.27)

Умножим выражение (3.10.27) на р (*)= Р (0 е

и

учтем, что*

ердер= 0 . Тогда запишем

 

 

 

 

 

 

Р(<)ДР(*)=Р(*)'ДР(0-

 

(3.10.28)

Подставим в равенство

(3.10.28)

выражение

(3.10.26):

Р (*) др (/)=

Р (t ) Д р (*)= Р {t) Дг (t).

 

(3.10.29)

Отсюда

 

 

 

 

 

 

 

др(0 = ^ - L r [ t ) = e ^ r { t )

 

(3.10.30)

или

 

р(0

 

 

 

 

 

 

 

 

 

 

 

 

ДР(*)=Г(^

■*/-—Д r{t).

 

(3.10.31)

 

 

р(0

 

 

 

 

Варьируя (3.6.13), получим

 

 

 

 

 

Д C O S f}j =

Ap(Qg/

 

p ( t )

Ар (t ) e j

 

 

 

 

РЧО

 

 

 

 

Р 0 )

 

 

 

 

Выражая здесь Л р (0

через Д г(/)

и используя

результат

(3.10.31), получим

 

 

 

r ( Q - r y (Q

 

 

д cos Q} ( t ) = -----

eJj~ cos в . (t )

 

(3.10.32)

 

L r ( t ) .

p(0 L.

 

 

 

P(0

 

 

Представим результаты (3.10.33) и (3.10.32) в матричной

форме. Обозначим для этого

 

 

 

 

 

д cos b j(t) =

A y a j {t)

 

( / =

И 2, 3);

 

(3.10.33)

 

 

 

 

 

 

 

* Поскольку ( e p )2 = 1, то 2cp Дср= о.

107