Файл: Свешников А.А. Вероятностные методы в прикладной теории гироскопов.pdf

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

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

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

Добавлен: 10.04.2024

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

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

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

§ 5.i 1 УРАВНЕНИЯ, СОДЕРЖАЩИЕ СЛУЧАЙНЫЕ ВЕЛИЧИНЫ 309

V j + k j a t j , VJ + 1 ,

. . . , v j . Разложив

это решение по степеням к ^ Ѵ}.

и ограничиваясь при этом только

первыми степенями этого

пара­

метра, получим

 

 

 

 

 

да (/)

 

 

 

 

 

 

 

(5.53)

 

 

 

 

 

 

дѵ 4 kJaT

 

 

 

 

 

 

 

Следовательно, искомая производная определяется равенством

 

da(t)

_

_

aj { t ) — â{t)

(5.54)

 

дѵj

a>>j

 

 

k j

 

 

 

 

 

 

где в соответствии с (50) a(t;

vlt . . .,

vm) заменено на a (t).

 

Подстановка

(54), например,

в (52)

дает

 

 

Dfa ( 0 1

=

 

 

а -( t) — а (г)

(5.55)

 

2

 

 

к,

Таким образом, мы получаем следующую схему определения дисперсии решения уравнения: исходное уравнение должно быть решено (пг+1 ) раз — первый раз при замене всех случайных ве­ личин их математическими ожиданиями, остальные гп раз — при замене /-й случайной величины на v .-\-kj^Vj, а остальных случай­

ных величин — их математическими ожиданиями. Результат пер­ вого решения дает <х(t), остальные решения позволяют вычислить D [ а (t) ] по формуле (55) (и по аналогичной формуле для Ка (tv t2)).

В формулу (54) входит неопределенная постоянная^., от значения

которой, по смыслу решаемой задачи, не должна зависеть левая часть этой формулы. Поэтому при выборе &. нужно руководство­

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

тем меньше будет зависеть (54) от значений этих постоянных, так как полученные выше формулы имеют в своей основе при­ ближенную формулу (53), которая тем более точна, чем меньше к . ш

С другой стороны, очень малые значения &. являются нежелатель­

ными, поскольку в окончательную формулу (55) входят разности (£)—ä (t), которые при малых к . будут также малы, и следо­

вательно, для получения достаточной точности необходимо реше­ ния уравнения находить с большой точностью. Поэтому значения /с.

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


310

У Р А В Н Е Н И Я ГУ СО С Л У Ч А Й Н Ы М И К О Э Ф Ф И Ц И ЕН ТА М И

[ГЛ . 5

Рассмотренная выше схема расчета позволяет учесть в разло­ жении (49) только линейные члены. Путем усложнения схемы можно учесть и следующие слагаемые ряда Тейлора [25].

Кроме рассмотренного выше способа учета случайного харак­ тера коэффициентов уравнения возможно получение разложения решения а (£; Ѵг, Ѵ2, . . ., Ѵт) по степеням ( V ѵ}) путем приме­

нения метода малого параметра, сущность которого заключается в том, что решение ищется в виде ряда по степеням «малых» откло­ нений (Vj ѵ), коэффициентами которого являются функции

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

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

 

[L (р) +

{V V) Q (/?)] а (t)=X (t),

 

(5.56)

где L (р) и Q (р) — полиномы с постоянными коэффициентами

оператора дифференцирования p = dldt, не содержащие

случай­

ных величин;

X (t) — случайная функция с известными

харак­

теристиками,

не зависящая от случайной величины V, а начальные

условия для

простоты будем считать нулевыми. В

соответствии

с общей идеей будем искать решение в виде

 

 

а (t) = а0 (t) +

(t) (F — ѵ) + а2 (f) (F — ^ ) 2 +

• ■•

(5.57)

Подставляя (57) в (56), приравнивая правой части уравнения сла­ гаемые, не содержащие разности (V—г?), и приравнивая нулю сла­ гаемые, содержащие эту разность в первой степени, во второй степени и т. п., получим систему линейных дифференциальных урав­ нений с постоянными коэффициентами

L (Р) а о (t) = х (0.

L {р) et] (t) = —Q (р) а0 (t),

(5 .5 8 )

L (р) а.. (t) = Q (р) ay_t (t),


§ 5.1]

УРАВНЕНИЯ, СОДЕРЖАЩИЕ СЛУЧАЙНЫЕ ВЕЛИЧИНЫ

ЗИ

Решение первого из этих уравнений может быть явно выражено через X (t) по формуле

t

 

a0( t ) = \ l ( T ) X ( t - z ) d r .

(5.59)

О

 

Подставляя (59) во второе уравнение системы (58), получим

возможность выразить через X (t) функцию

(t) и т. д.

Таким образом, все функции а (t) линейно выражаются через

правую часть исходного уравнения X (t), и следовательно, мо­ менты ординат этих функций могут быть выражены через моменты ординат функции X (t) того же порядка.

Возвращаясь к (57), замечаем, что любые моменты а (t) просто выражаются через моменты случайной величины V и моменты орди­ нат случайных функций а . (t). Например, для математического

ожидания и дисперсии имеем (считаем V не зависящей от X (t))

ÄW = äo( 0 + «2 (*)a2 + - ••>

 

I

D [a (i)l = D К (t)] + {D [aj («)] + к

(t)f +

(5.60)

+

{t, t)} D [H] +

• • • j

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

5. Гироскоп направления. Рассмотрим применение выведен­ ных выше общих формул к исследованию отдельных ГУ.

Начнем с рассмотрения азимутального ухода ГН, описы­ ваемого уравнением (2). Это уравнение типа (12), где положим

(3(0) = 0 ,

в = ^ ,

(5.61)

X ( 0 = l + y W 5 (i).

(5.62)

Будем считать кинетический момент Н и смещение центра тяжести Іг гироскопа независимыми случайными величинами. Допуская


312

У Р А В Н Е Н И Я ГУ СО С Л У Ч А Й Н Ы М И К О Э Ф Ф И Ц И Е Н Т А М И

[ГЛ. 5

линеаризацию (61) относительно отклонения (Н К) и учитывая,

что

М [Іг1=0, получим

 

 

 

4 = 0 ,

=

(5-Г.З)

Принимая для вертикального ускорения точки установки ГН выражение (2. 34), получим

X(t) = 1 + j [Сс (t) +

Z. (t) + уЬ (t) - хф (t)].

(5.64)

Следовательно, учитывая обозначение (14), имеем

 

Y ( t ) = t - b y fCс ( 0 + zB(t) + ф (t) — хф (г)] —

 

 

 

[Сс (0 ) +

і в(0 ) + ф (0 ) -

хф (0 )|, '

(5-65)

y{t) = t,

 

 

 

 

 

 

D [Y (г)] = - j r

{[^сс (0) +

К,я(0 ) + у % (0 ) +

(0 )] -

 

 

\к,с (t) +

(t) + y2Kf) (t) +

(01}-

(5.66)

Применение к

данному

случаю

формулы

(17) для

дисперсии

ухода ГН а (t)

дает

 

 

 

 

 

D K *)] = | f

(l

 

 

+

(5.67)

T. e. дисперсия ухода ГН, вызванная случайной статической не­ уравновешенностью гироскопа, растет со временем по квадратич­ ному закону, будучи пропорциональной дисперсии небаланса гироскопа Iг. Наличие случайных ускорений места установки ГН приводит только к постоянному слагаемому в величине дисперсии а (t), т. е. не вызывает нарастания ошибки со временем.

Аналогичным образом исследуются и уравнения (1) и другие уравнения ГУ, характеризуемые одночленными линейными урав­ нениями первого порядка.

6 . Гировертикаль. Уравнения (3) для ГВ с маятниковой кор­ рекцией являются двучленными линейными уравнениями пер­ вого порядка, т. е. принадлежат к уравнениям типа (2 0 ), где в дан­ ном случае (рассмотрим первое из уравнений (3))

A i = W

^ = 1 -- ^(<)= Хі (0 + 4 - № + <?»^пф(0]. (5.68)

Для

нахождения â (t) и D

[a (t) ] применим разложение (27)

и следующие из него формулы

(28) и (29), однако при этом будем

учитывать,

что в данном случае В = А Ѵ и поэтому величины А х

и В нельзя считать независимыми.


§ 5.1]

У Р А В Н Е Н И Й , С О Д Е РЖ А Щ И Е С Л У Ч А Й Н Ы Е В Е Л И Ч И Н Ы

313

Линеаризуя выражение для А х относительно отклонения

(Н h),

получим

 

Находя математическое ожидание и дисперсию выражения (27), с учетом только членов второго порядка относительно ойі получим (считаем а ( 0 ) = 0 )

 

* (0 = ЬУ (!) aWi (t) + Y

(t),

 

(5.70)

D [а (t)] =

62Ä ff (t,

t) + ol[Ky (t, f) + i/ 2 (£) +

(C

*) +

 

+ b2yl(t) t ö R yWi(t,

t) 25y (t) уг (t) +

b \ R y y t (t,

* )],

(5.71)

где Y (t),

Fj (t) и

F2 (t)

определяются

(30)

и, следовательно, у (t),

9i(t), УА*)’

К уЬ, t), R „„,(*>

*) и КѵА*і г) могут

быть

выражены через математическое ожидание и корреляционную функцию X (t), определяемую формулой (6 8 ).

Уравнение инерциальной вертикали (7) является линейным уравнением второго порядка, не содержащим первой производной. Следовательно, при его исследовании можно использовать фор­

мулы, полученные выше для уравнения (35), если положить

 

ja=

= 0, v =

j / ^ - g 0,

X ( t ) = v N(t){^ — ± ) .

(5.72)

Подставляя (72) в (36),

для данного^ случая

получим

 

“ W =

( і г - т ) S * *

- х) * ™

Ѵ Ч г х d x '

<5-73)

где производную vN (t) будем считать неслучайной функцией. Разлагая последнее выражение по степеням отклонений (;лх—p.j) и К), сохраняя при этом члены не выше второго порядка мало­ сти и учитывая, что при выполнении условия невозмущаемости

 

 

 

йо

 

(5.74)

 

 

 

Я

 

 

 

 

 

 

получим

 

 

 

 

 

1

йо

До_В Я___1

((J.Q йо)

 

■(*) = ^Rgo

Йо

Я

2

Йо

 

 

. 1

(Н—К)2

*

 

xj d-z +

 

 

Я2

j Vn (t — т) Sin (^J/

 

 

)

 

 

 

 

 

 

 

 

 

 

 

t)t cos

(5.75)

+ т ( Ѵ

- £ т £ ) , - т ! » » ( ‘