Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Свешников А.А. Вероятностные методы в прикладной теории гироскопов

.pdf
Скачиваний:
29
Добавлен:
27.10.2023
Размер:
18.8 Mб
Скачать

§ 4 . 6 ]

ПРИМЕРЫ

НА ИССЛЕДОВАНИЕ ГУ

297

Воспользовавшись формулами

 

 

 

 

 

sh а sh ß =

[ ch + ß) — ch (а — ß)]/2 ,

 

ch а ch ß =

[ch (а — ß) +

ch (а -f- ß)]/2 ,

 

sin а sin ß =

[cos (а — ß) — cos (а +

ß)]/2 ,

 

cos а cos ß =

[cos (a — ß) -|- cos (а -f- ß)J/2 ,

последнюю формулу можно переписать в виде

 

D [а (i)l =

t

t

 

 

 

 

 

 

{[ch v2 (Xl + x2) — ch v2 (x2 — xx)] X

 

0

0

 

 

 

 

 

X [cos V,

(X j + X ,) + cos Vj (x2 X j) ] К ъ (x2

X j)-f [ch v2 (xL -j- x2) +

+ ch v2 (x2 — Xj)][cos vj (x2 — Xj) — cos Vj (x2 +

x^]

(x2 — xx)} dxxdx2.

Переходя от переменных интегрирования х2 и х2 к перемен­

ным интегрирования х = х2 — х1;

£ = х2 -[-х1, получим

 

t

(2t—т

COS

 

 

 

D [а (0.1 =

Ш 5

— ТМ f ^ Ch

+

 

 

 

О

I

X

 

 

 

 

-)- ch v2£ cos VjX— cos vj ch v2x — cos v2x ch v2x) Kb (x) -)- (ch v2£ cos v j—

— ch v2£ cos VjX-(- cos Vjl ch v2x — cos ѵхх ch v2x) (x)] d l J dx.

Интегрирование по £ может быть выполнено нѳзависимо?от

вида корреляционных функций К ь (х)

и

К-, (х).

Если 'при' этом

учесть малость |ѵ2| сравнительно с

ѵ1(

уа и ут,

то послеіинте-

грирования

D [а (0]

примет вид

 

 

 

D[a(f)] = -

Im 0! HFH cos Vj (2t x) —

 

 

jlL

 

 

 

 

2x

2

2

л

 

— Щ cos v,x

^ sin VX(2t x) +

-g- Sin vxxj K b (x) +

+(sin Vj (2t — x) — sin Vjx) — 4 (t — x) cos VjxJ ä ^ (x)J dx.

При заданном в примере

значении времени t, которое зна­

чительно превосходит время

корреляции

случайных процессов

&(f) и f (t), верхний предел в последнем

интеграле может быть

положен равным сю. Поэтому, принимая во внимание, что для

принятого вида

корреляционных

функций Кь (х) и

(х) после

Кь(х) =

—а$ (у! + Х|) e-w>M (cos Хэх — g

sin Xe | x |) ,

K.<(x) =

—o| (u.2 -f- X?)

(cos XTx

r

smXT

 

 

 

 

 

298

ГУ , О П И С Ы В А Е М Ы Е

Л И Н Е И Н Ы М И

У Р А В Н Е Н И Я М И

[ГЛ. 4

выполнения интегрирования по т получим

 

 

 

 

 

D [а (0 ] =

{J[(4fsin2V

+ ^ cos 2Ѵ +

т)

 

 

 

— 2 sin 2vjf • 7(ö>-f- ^4г cos 2vji

 

sin 2 v ^

7(ö) —

 

 

 

 

2 ( 1 + cos 2 vxt) /(») +

■(cos 2Vji — 1) 7(f) -f-

 

 

где

 

+

i

(sin 2vjt -

2v/) 7(f) +

47(f)]},

(4.595)

 

 

 

 

 

 

 

 

 

 

 

 

7f)

 

sin VjTdx =

 

 

 

 

 

 

 

 

 

 

_64 + 4)

 

f4 +

 

+ 4 ,

~~ 4 ~ 4

 

 

 

 

 

2 Xa

 

f4 + (vi + ^a) 2

4 + (7i — ^a) 2 .

7(9j __ — ^ К ъ (т) cos

 

=

 

 

 

 

 

 

 

 

 

о

_ °\ 64 + 4)

 

 

 

 

щ ~

ь

 

 

 

 

- ѵі

 

 

 

 

 

Щ

 

4 + (vi + ^s) 2

4 + Ha —vi)2

7(ö) =

■— ^ TCa (t)t sin v1 t^t=

 

 

 

 

 

 

 

 

 

__ 2aI (p-I H-

}

 

p.9____________ fT9_______

 

 

 

 

2l*

I 4 + (Vi + )ч>)?

4 + Hi - h)2

 

 

 

 

 

 

_

^a 64 + У , + 4)

p9 (^x9Ц

 

 

 

 

 

 

[p& + (vi + ^a)2] 2

Ha + Hi — ^a) 2] 2

J ’

 

00

 

 

 

 

 

 

 

 

 

 

 

7(*) =

— ^ Äs (t) t cos vxxdx =

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

_

4 6*5 + 4) I

— У!_____ |_

 

2 \д — Уд

.

 

 

 

 

 

1 Ра + Hi + ^s) 2

4 + На —vi) 2

 

 

 

 

 

 

 

 

2 vjfia

 

 

2

(2 Xa — vj) 4

I

 

 

 

 

 

[4 + Hi ■+ ^a) 2] 2

[4 + На — vi) 2] 2

J

a 7(f),

7|f),

7(f) и 7(f)

отличаются от

7(9),

J ^ \

J ^ \ 7(®) соот­

ветственно только тем, что всюду индекс «&» нужно заменить на индекс «у».

Подставляя величины 7(s), 7(&), 73(9), 7(e), 7jf), 7|f), 7(f), 7(f) в (595), получим окончательное выражение для D [<* (£)].

Г Л А В А Й

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

§ 5Л. Линейные уравнения, коэффициенты которых являются случайными величинами

1. ГУ, характеризуемые линейными уравнениями со случай­ ными коэффициентами. В прикладной теории гироскопов часто приходится иметь дело с системами линейных дифференциальных уравнений, коэффициенты которых являются случайными вели­ чинами.

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

Простейшим примером уравнений ГУ, содержащих случайные величины, являются уравнения (3. 1 0 ) трехстепенного астатиче­ ского гироскопа в рамках прецессионной теории, имеющие вид

(5.1)

где Н и MXii M'. можно считать случайными величинами.

К такому же типу относится и уравнение (3. 29), характеризую­ щее азимутальное движение ГН при статической неуравновешен­ ности гироскопа

(5.2)

где вес гироскопа Р, кинетический момент гироскопа Н и смеще­ ние Іг центра тяжести гироскопа относительно точки подвеса являются случайными величинами.

300

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

Несколько более сложны уравнения (3. 62) ГВ с маятниковой коррекцией, которые после учета моментов сил трения в осях под­ веса примут вид

Т Ь + а = Х і

т§+ Р = * 2

где

(0 + № + Q, signФ(<)1.

(5.3)

(0 + 7/-[Щу + Qysignб (01.

(5.4)

а Я

и S — случайные величины.

(3. 164)

Уравнение поплавкового интегрирующего гироскопа

(при

отсутствии возмущающих

воздействий)

 

 

Jи. гР +

ЩНш^,

(5.5)

где J„ т, b и Н — случайные величины, с математической точки зрения отличается от (3) только тем, что является уравнением первого порядка не относительно ß, а относительно ß.

Из ГУ, характеризуемых дифференциальным уравнением вто­ рого порядка (или системой двух уравнений первого порядка), можно указать гиротахометр, уравнение которого (3. 116) имеет вид (считаем <і>=0, М т= 0, М = 0)

К. J +

bß + Cß = //cV

(5.6)

инерциальную вертикаль [см.

(3.105)]

(положим р =

р0)

 

 

 

(5.7)

вибрационный гироскоп (см. (3. 146))

 

 

/.ft -f- Ъд. -)- са = Нш

(t) sin Qt,

(5.8)

гиромаятник в рамках прецессионной

теории [см. (3. 80)]

 

 

 

(5.9)

и ряд других уравнений, коэффициенты которых являются случай­ ными величинами.

Из уравнений более высокого порядка, чем второй, можно ука­ зать на инерциальную вертикаль при наличии статической не­ уравновешенности гироскопа [см. (3. 1 0 1 ) 1

§ 5.1]

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

301

где

 

 

 

к = і г ’ ѵ8 = §->

(5Л1)

а I, Р и Н — случайные величины.

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

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

2.Уравнение первого порядка со случайными коэффициен­

тами. Перейдем к рассмотрению различных типов уравнений со случайными коэффициентами.

Уравнение первого порядка, не содержащее зависимую пере­

менную, в общем случае может быть представлено в виде

 

я = B X (t),

(5.12)

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

Будем считать В независимой от X (t). Интегрирование урав­ нения (1 2 ) дает

a(t) = a0 + B Y (t),

(5.13)

где обозначено

 

t

 

Y ( t ) = \ X ( t l)dtl.

(5.14)

о

 

Применяя к (13) операцию нахождения математического ожи­ дания и учитывая взаимную независимость а0, В и Y (t), получим

<*(*)= «о

(5.15)

Аналогичным образом, находя корреляционную функцию (13), получим

Ku{t,, і2) = 0 Ы + {І? + 0 \ В ] } К ^ і ѵ Q + y itJ y iQ D iB ] . (5.16)

Полагая в последнем равенстве tt=t^=t, для дисперсии угла а (t) будем иметь

D [а (*)]= D ы + + D \В\] D [Г (t)\ + [y{t)f D [В]. (5.17)

302

У Р А В Н Е Н И Я

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

[ГЛ.

5

 

Выражение для

D [ а (і) | в том случае,

когда коэффициент

В

в уравнении

(1 2 ) не является случайной

величиной и

равен

Ъ,

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

равенством

 

 

 

 

 

D [ a ( i) l - D |a 0] + ö2 D |F (0].

(5.18)

Сравнивая (17) с (18), находим, что увеличение дисперсии угла а (t), возникающее вследствие случайности коэффициента В , определяется формулой

ÖD (t)]= D [В\ 2 (t) + D [У (01),

(5.19)

где выражения для D [Y (£)] были получены в § 4.1.

Уравнение первого порядка, содержащее искомую переменную,

имеет вид

 

А -j- A-p = ВХ (t),

(5.20)

где Аг и В — случайные величины, математические ожидания кото­

рых

йл и Ъ и элементы корреляционной матрицы ки = D [)4, |,

к22=

D [В ] и кп = М [(i4j — (Д) — 5)] будем считать известными.

Решение уравнения (20) в соответствии с (4.80)

и (4.81)

имеет вид

 

 

t

 

 

а (t) — а (0) -f- В ^ е~А'хХ (t х) dx,

(5.21)

 

о

 

где, в отличие от (4. 81), А 1ш В являются случайными величинами. Для определения моментов ординат случайной функции а (t)

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

t

 

А (t) = А (0) 4- ЪJ E ai (ix) x(t т) dx,

(5.22)

о

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

Еа, (и) — exp I —у olp2-f iäpij

(5.23)

и формула (2 2 ) принимает вид

t

A (t) = А (0) -)- В j е 2 °а' a^ x ( t x) dx.

(5-24)

о

§ 5 . І ] УРАВНЕНИЯ, СОДЕРЖАЩИЕ СЛУЧАЙНЫЕ ВЕЛИ ЧИ НЫ 303

Аналогичным образом вычисление дисперсии (21) дает

D [а (0] = D (.« (0)] +

t t

+ (D [В] + В2} j j Ещ(ix, + ix2) Kx (t — x„ t x2) cM x2 +

0 о

tt

+D [B] ^ Je1(ix, + ix2) x (t — X , ) Ж(i — X 2 ) dx,dx2+

0 0

t t

-f 62 jJ[2?(ix,+ ix2)— ■E (ix,)jf?(ix2)Iж (t — X,)X(t x2)dx,dx2 (5.25)

0 0

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

t

t

 

 

 

D [ос (t)] = D [а (0)] + (D [B\ + £*}$$ exp

<

(x, + x2) 2 -

 

0

0

 

 

 

— ä] (x,+

x2)j Kx(x2

X,)dx,dx2.

(5.26)

Формулы (24) и (25) требуют пояснения, так как интегралы, входящие в эти формулы, растут с увеличением t, а для наиболее часто встречающихся случаев постоянного математического ожи­ дания X и экспоненциального затухания корреляционной функции Кх (т) растут неограниченно при t -> со. Этот результат про­ тиворечит физическому смыслу задачи и связан с тем, что предпо­ ложение о нормальном законе распределения случайной вели­ чины Л, в данном случае неприемлемо и должно быть отброшено. Действительно, предполагая Л , нормальной, мы тем самым с ко­ нечной вероятностью допускаем, что эта величина может прини­ мать отрицательные значения. Однако при Л, <С 0 уравнение (20) соответствует неустойчивой динамической системе, для которой обобщенная координата а (і) будет расти со временем. Поэтому, производя при вычислении і (і) и D [ а. (і) ] усреднение по различ­ ным значениям А „ мы и получим расходящийся результат. В дейст­ вительности в реальных гироскопических системах коэффициент Л, всегда положителен, т. е. рассматриваемая система всегда устой­ чива. Таким образом, предположение о нормальном законе рас­ пределения, которое обычно является допустимым при решении подобных задач, в данном случае неприемлемо и должно быть отброшено. Полученный результат, представляющий интерес с общетеоретической точки зрения, в задачах прикладной теории гироскопов обычно не играет существенного значения, так как разброс коэффициентов уравнения (2 0 ) вызывается неточностью

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

изготовления прибора и является малым в том смысле, что реше­ ние уравнения (2 0 ) можно разложить по степеням разности (Нх—йх) и воспользоваться более простыми приближенными формулами, применимость которых основана на допущении, что ни при каких реальных значениях А г система не может потерять устойчивость.

Итак,

считая

разность (Лх— ал) малой и разложив е~А'т=

= <г“‘т •

в ряд по степеням (Лх — «x), вместо (2 1 ) будем иметь

 

t

t

а (t) = а (0) -f- В I е~а'тХ (t — х) dx х — а:) В ^ e~a'zxX (t — х) dx -f-

о

u

 

t

 

 

+ y (A, ä,)2 В \ e-ä^ X

(t — x ) d x ~ . . .

(5.27)

u

Полученный ряд позволяет вычислить моменты а (t) с любой точностью. Так, например, ограничиваясь членами второго по­ рядка малости относительно {А±—вх) (т. е. сохраняя в окончатель­ ных выражениях центральные моменты случайной величины A t не выше второго), получим

 

 

 

t

 

 

 

t

 

 

â (t) =

£ (0) -f- b ^ e~ä^x (t — x)dx + Y

^ e~ä‘Tx2x (t x) dx,

(5.28)

D [<* m

=

D [a (0)3 +

D [в ] (D [F (*)] +

y2{t) +

 

 

 

+

D Ш

{ D [ F x(0] +

F ( t ) +

(t , t) +

g (t) y 2 m

+

+

52 (D lF(i)l +

D [ni]{D [F1 (i)] +

y ? (0 + /? w,(^ *)}).

(5-29)

где введены обозначения

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

F (f) =

j e-ä^X {t — x) dx,

Y1 {t) =

 

 

 

 

0

t

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

= J e~“‘TxX (t — x) dx,

Y2(t) =

^ e~ä'*x2X (t — x) dx

(5.30)

 

 

о

 

 

 

 

0

 

 

и, следовательно,

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у (t) == ^ er^x (t x) dx,

 

(5.31)

 

 

 

 

 

t

о

 

 

 

 

 

 

 

 

t

 

 

 

D [F (*)] =

*,(*,

0 =

i

\ e - ^ ^ K

x( t - x v

t - x 2)dx1dx2 (5.32)

о0

иT. Д. Полученные формулы существенно упрощаются в частных

случаях.

Например,

если а(0) = 0,

аь = 0 , имеем

 

 

t

 

 

x (t — x) dx,

 

â(t) = b j е_а‘т 1

+

(5.33)

 

о

 

 

 

 

D |а (0) =

V {Ку (t,

t) +

о*, [Ку, (t,

t) + y\ (t) + R m (t, *)]}•

(5.34)

§ 5 . 1 ]

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

305

3.Уравнение второго порядка со случайными коэффициен­

тами. Перейдем к рассмотрению уравнения

второго порядка,

т. е. уравнения вида

 

ос -(- /Ijd -j- А2а = X (t),

(5.35)

где коэффициенты A t я А 2 являются случайными величинами, не­ зависимыми от ординат случайной функции X (t). Будем считать для простоты начальные условия нулевыми. Тогда в соответствии с (4. 147) и (4. 150) решение уравнения (35) мо?кет быть предста­ влено в виде

t

 

 

ос (t) = - ^ I е~ѵ-тХ (t — т) sin ѵт dr,

(5.36)

о

 

 

где

________

 

р = 4 л г,

Ѵ= ] / Л 2 - І Л ? .

(5.37)

Если закон распределения системы случайных величин А ѵ

А 2 известен, то можно считать известной и плотность вероятности

/(р, ѵ) системы случайных величин р, ѵ. Поэтому все моменты ординат случайной функции а (t) могут быть вычислены. Например,

для математического ожидания а (t) имеем

t

СО

 

 

а (t) =

^ ^ -і- е~ѵ-тх (t — т) sin ѵт • / (р,

ѵ) dp. dv dz.

(5.38)

0

— со

 

 

Однако формулы такого типа практически оказываются не­

пригодными для

вычисления, так как нахождение

плотности

/ (р, ѵ) по заданному закону распределения

системы случайных

величин (Аг, Л2) связано с определенными трудностями (даже для нормальной системы величин получается слитком громоздкое выражение), а вычисление интегралов типа (38) обычно возможно только численное.

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

ням «малых» величин

(А х—ах) и (А2—а2),

получим

 

І* = у

йі + у ( 4 і ~ ®і)>

 

(5.39)

 

< -• У ^

Ш

л і

‘ (5.40)

 

2 Jj

! ѵ9^

1

 

j= 1

20 А . А . С в е ш н и к о в . С. С. Р и в к и ң

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

где введены обозначения

 

 

 

■V

1

Ло = Ап а„ ■

 

(5-41)

T ß

Поскольку моменты случайных величин А ги А 2 предполагаются заданными, можно вычислить и все моменты системы случайных величин (р, ѵ). Например,

 

 

 

 

СО

 

«к

=

ѵ" =

ѵо - 2 М К - ^ | т г М ( Л І 1 ,

(5.42)

где

 

 

 

j=1

 

 

 

 

 

 

М И 3] =

-

і < ,

M f^] =

a^ + iâ f a li + i â 1 ^ . ) + 1 L !A(«1)

(5.43)

и т. д.,

где

— центральный момент порядка / случайной вели­

чины Л 1( а последняя формула написана для случая, когда

и

А 2 — независимые

случайные величины. Подстановка (39) и

(40) в (36)

позволяет получить разложение а (t) по степеням слу-

чайных величин у (Лх—■äj) и Л3, которые в данной задаче можно считать малыми. Производя это разложение и сохраняя только

первые степени «малых» случайных величин, получим

 

 

t

)

 

* (t) «

j с

2 а,Т {sin ѵ0т + 2 ^ ( ѵ fos ѵ0т — sin ѵ0х) л з —

 

 

о

— у х (Лх — äj) sin v0xl X (f — x) dt.

(5.44)

 

 

Находя математическое ожидание и дисперсию последнего выражения, учитывая при этом (43) и сохраняя в окончательном

выражении

центральные моменты случайных величин А г и А %

не выше вторых, получим

 

 

 

 

â (i)= — J e

8

(sin v0x+ ^

(—v„tcos v0t+

sin v0x)| X * (t - x) dx,

 

 

 

 

 

 

 

(5.45)

D L« (*)] =

4

f f e

2 a,(T'+T2)G (xlf x2) Kx (t — x„

t — x2) dxjdx2 — [â (i

 

 

 

 

 

 

 

*(*)]2,

 

 

 

 

 

 

 

(5.46)

где G (xj,

x2) — тригонометрический полином, содержащий

sin ѵ0х,,

cos v0x1; sin""v0x2, cos v0x2, в коэффициенты которого входят

ові, a^,

äj и â2.

D [X (£)]

можно считать малой величиной,

как и диспер­

Еслй

сии ajj и

a^,

то,

отбрасывая

члены, содержащие

произведения

Соседние файлы в папке книги из ГПНТБ