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

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

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

I 4.3]

УРАВНЕНИЕ 2-ГО ПОРЯДКА

211

т. е. после интегрирования

 

 

ß (0 =

^

г- sin2 т + £ (* — т sin ѵ0 •

(4.246)

Определяя математическое ожидание и дисперсию

выражения

(246), пользуясь общими формулами (1.38) и (1.39), получим

ß (*) =

^

Sin2 £ + i§(« -

I sin v f),

(4.247)

D [ß (*)J =

^

sin4 у ofe^ + 1 ^

i sin vi) 2 < .

(4.248)

Выражения (247) и (248) существеннощтличаются от аналогич­ ных выражений (231) и (234), полученных для того случая, когда в правой части уравнения типа (137) стоит случайная функция времени.

Особенно существенное отличие имеет место для дисперсии: при наличии в правой части уравнения стационарной случайной функции дисперсия, помимо осциллирующих членов, содержит слагаемые, зависящие от времени линейно, в то время как наличие случайной величины ЬаѴ[ приводит к возникновению только коле­ бательных слагаемых, а случайная величина ех, входящая в правую часть уравнения (142) в виде множителя у і, (т. е. s—ехі) добавляет

квыражению дисперсии слагаемое, квадратичное относительно t.

Втом случае, когда ошибку е (і) нельзя представить в виде произведения случайной величины на время, а нужно рассматри­ вать как нестационарную функцию времени, математическое ожи­ дание и корреляционная функция ß (і) могут быть найдены по общим формулам (153) и (156), в которых в данном случае следует положить р =0, ѵ0 =ѵ. Упрощение вычислений в этом случае в об­ щем виде уже невозможно, и для получения окончательного от­ вета необходимо располагать видом корреляционной функции и

математического ожидания правой части уравнения.

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

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

218 ГУ, ОПИСЫВАЕМЫЕ ЛИНЕЙНЫМИ УРАВНЕНИЯМИ [ГЛ. 4

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

Если правая часть уравнения (124) не является нормальной, то закон распределения решения уравнения также не будет нор­ мальным и его определение возможно только приближенное, например, с помощью использования ряда Шарлье (1.44), коэф­ фициенты которого содержат моменты ординат решения уравне­ ния (124) более высокие, чем вторые. Выражения для этих моментов могут быть получены из формулы (152) (или общей формулы (150)) тем же методом, каким были получены формулы для математиче­ ского ожидания и корреляционной функции (дисперсии) решения уравнения, однако кратность интегралов, через которые выра­ жаются эти моменты, повышается с ростом порядка момента, а в подынтегральные функции входят моменты правой части уравнения до того же порядка, что и вычисляемый момент.

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

§ 4.4. Система двух уравнений первого порядка, сводимая к одному уравнению первого порядка

скомплексными коэффициентами

1.Примеры ГУ, описываемых уравнениями данного |типа.

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

Уравнения этого типа в общем виде можно записать следую­ щим образом:

Ri

.

F 2 (*), J

(4.249)

Р +

«і« + #2а =

 

где коэффициенты ах и а2 и правые части уравнений Fx (t) и F2 (t) являются случайными или детерминированными функциями вре­ мени.

§ 4.4] СИСТЕМА ДВУХ УРАВНЕНИЙ 1-ГО ПОРЯДКА 219

Умножив второе уравнение системы (249) на мнимую единицу і

и сложив уравнения, получим

 

 

 

(1 + а2і) (А -|- ф) +

іаг (а -f- iß) =

F1(t) +

iF2 (t).

(4.250)

Если ввести обозначения

а 1 ( a 2 + 0

 

 

В=

а + iß,

 

 

« 1

+ 1

 

 

 

 

 

 

V

(F\ + а<^\) + 1 ( ^ 2

— «2 ^ 1 )

 

(4.251)

А W — ----------- ------------------

 

 

то вместо уравнения (250) получим

 

 

 

 

&(*) +

а (08(0 = Х(*),

 

(4.252)

где коэффициент а (t), функция X (t) и функция

8 (()

имеют ве­

щественные и мнимые части.

 

 

порядка

Таким образом,

система

(249) уравнений первого

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

Ксистеме уравнений (249), а следовательно, и к уравнению

(252)сводятся уравнения движения различных ГУ.

Так, например, система уравнений (3.84) для гиромаятника в рамках прецессионной теории с учетом возмущающих момен­

тов М х и М2

имеет вид

* “ { * [ ! + 7 ^ (0 - < p (0 }P = - 7 ^ 4W + am s,

 

(4.253)

P + H 1 + 7 ^ W ? w} <*= y W ( ( t ) + hMv

где (t),

(t) и Wc (t) — составляющие линейного ускорения

точки подвеса маятника.

Система уравнений (253) аналогична системе (249) и, следова­ тельно, может быть сведена к уравнению типа (252).

Если в ГМ для осуществления демпфирования используется коррекция от физического маятника с линейной характеристикой, определяющего углы отклонения оси гироскопа от вертикали с ошибками ул (t) и у2 (t), то система уравнений движения гиромаят­ ника принимает вид [см. (3.86)]

* + x a - [ Ä ( l + l W c) - « p (0]ß =

- | W 4 + xXl(i)>

 

(4.254)

^ + х Р + [ л ( і + і ^ с) - ф ( і ) ] « =

+

Эта система по типу не отличается от

системы (253),

220

ГУ, ОПИСЫВАЕМЫЕ ЛИНЕЙНЫМ И УРАВНЕНИЯМИ

[ГЛ. 4

 

Вторым примером подобной системы может служить

система

уравнений *), описывающая нутационное движение гиромаятника

 

 

T$ — â = — qM2,

 

 

(4.255)

 

 

 

+ Р=

,

 

 

 

 

 

 

 

 

которая путем замены

переменных

 

 

 

 

 

ccj=

&,

 

ßj= ß

 

 

(4.256)

сводится

к системе типа (249) и, следовательно,

к уравнению (252).

Аналогичным образом

преобразуется

и

система

уравнений

(3. 193)

для гироскопического

интегратора линейных

ускорений

(если в ней положить —mlwx— М2 = —Af.T]; Мл— Му)\

 

 

 

h . $ — Hâ. = ~ M Xi, I

 

(4.257)

 

 

 

+ Щ = M9l,

[

 

так как,

положив

 

 

 

 

 

 

 

 

 

 

получим

а1 = Ѵ^ГТ)

ßl — V

г. э Р>

(4.258)

 

 

 

 

 

 

 

 

 

 

ѴJ7~1T-

 

ѵт~ M,Ji’

 

(4.259)

 

 

f\J

Г . 3

 

J Г7]

 

 

 

Pi'

__H_

 

 

■MXl,

 

 

 

 

 

'JTH

 

 

 

 

 

r7j,f r. э

 

 

 

 

 

T.e. систему типа (249).

Ктакой же системе уравнений сводится система уравнений (3. 6 ) для трехстепенного астатического гироскопа

/г. J — Но. cos ß0 = — MXl,

(4.260)

Щ cos ß0 = M,.

Если имеет место сухое трение и для моментов сил трения спра­ ведливы формулы (2.103)

Mrx = M0TX+ QxSignQ(t), Мч = М І у — (>,sign*(f), (4.261)

то система (260), после замены переменных (258), не отличается от системы (257) и, следовательно, может быть приведена к стандарт­ ному виду (252).

Если имеет место жидкостное трение, то в соответствии с (2.101)

= п2+ Ö).

М ч = — Щ (â + tp)

(4.262)

*) Она может быть получена из

(3.77), если в ней принять 7 = 0, не

учитывать моментов инерции колец подвеса (7Гіэ = / г7, = / э)

и обозначить

/„

1

~Н==Т0

Я = ? '

§ 4.4] СИСТЕМА ДВУ Х УРАВНЕНИЙ 1-ГО ПОРЯДКА 221

и система (260) примет вид

•^г.эР — cos ß0 + n2ß = —n j,

(4.263)

/ Г|Д 4 - #ß cos ß0 + rijâ = —«jep.

Эта система также может быть сведена к уравнению вида (252), если коэффициенты nY и п2 пропорциональны моментам инерции / г э и / г?, т. е. если

Пі __ п2

(4.264)

7 гС ] г. э

 

Вэтом случае, использовав замену переменных (258), систему

(263)можно представить в виде

 

Н cos ßo

 

ß, +

Ш»! :

 

Ѵ^гС-^г. э

 

 

Р,

Н cos

ai +

raßi:

с'г.э

 

 

 

 

-<P.

(4.265)

:0 .

Умножая последнее из полученных равенств на і и складывая, получим

ft __ Н cosßo йі + иді =

В .

тВ

(4.266)

гИ г. г

ѵ'/.с/

 

 

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

 

 

(4.267)

8і = ®і + ißi-

 

 

Уравнение (266) не отличается по типу от уравнения (252), причем в этом случае

а = ti­ ll cos ßo

X (t) = — n 9

+

.

(4.268)

^ T . в

\ Вѵ ^ г . э

Iv'B^ гС

 

Так же к уравнению (252) сводится и система уравнений (3.19), описывающая прецессионное движение ГН на подвижном основании

а — нВ ==—в.

Г]Г с

Р + и <х.= —в5■

М

н

(4.269)

л

где u^, ип и в? — проекции переносной угловой скорости системы отсчета, а МХі и МѴі — возмущающие моменты, имеющие обычно случайный характер. Если объект перемещается с постоянной линейной скоростью ѵ параллельно земной поверхности и оси системы отсчета связаны с его траекторией, то

и^ = -—U сок <рsin К — , нт = U cos cp cos К, u^— U sin <р,

(4.270)

222

ГУ, ОПИСЫВАЕМЫЕ ЛИНЕЙНЫ М И УРАВНЕНИЯМ И

[ГЛ. 4

где

U — угловая скорость вращения Земли, а закон

изменения

широты места tp и курса К зависит от характера движения объекта, который для (270) соответствует ортодромии.

К уравнению (252) сводится и система уравнений ГВ с маят­ никовой коррекцией, которая, например, при наличии жидкост­ ного трения и одинаковых коэффициентов трения для обеих осей

(n1 =ra2 =ra) имеет вид [см.

(3.70)]

â + xoc— ^

р = *хі (О+ 'тг'ИО.

(4.271)

+= %(*) + 7г0(0, .

что не отличается по типу от системы уравнений (265) и, следо­ вательно, сводится к уравнению

8

+ а

8

= Х(г)

8

*ß),

(4.272)

 

 

( =<x +

 

где в данном случае

(4.273)

Для ГВ с учетом сухого трения в осях подвеса и при наличии рыскания корабля с угловой скоростью <р(t) справедлива система (3.72), имеющая вид

* +

*« + ?(*)£ =

*Хі У) +

7 Г

+

& siS n Ф)>

(4.274)

 

 

 

7

 

.

ß +

— Ф(0 * =

Ч і (*)+

+

Qy sign 0)> .

 

которая по типу не отличается от системы (271) и, следовательно, также сводится к уравнению (252).

В качестве последнего примера системы уравнений рассматри­ ваемого типа укажем на уравнения гировертикали (§ 3.2, п. 4), основанной на использовании трехстепенного астатического гиро­ скопа без коррекции, при наличии вращения подвеса

(4.275)

4.4] СИСТЕМА ДВУ Х УРАВНЕНИЙ 1-ГО ПОРЯДКА 223

и для гировертикали с линейной маятниковой коррекцией

\

â + ха — coß = -JJ- М ѵ

(4.276)

ß + xß + <оа — ~ң М ѵ

где М х и М 2 — моменты, вызванные различными возмущениями. И система (275) и система (276) путем введения комплексной

переменной 8 = a - fiß сводятся к уравнению типа (252).

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

t

 

 

 

j а(Ц)dtx

 

(4.277)

8 (* )= 8 (0 О

+

X (tx) dt„

 

где 8 (0) — начальное значение переменной 8 (t) на основании (251) связано с начальными значениями углов а и ß соотношением

8(0) = a(0) + Jß(0).

(4.278)

Уста овим связь между первыми двумя моментами ординат случайной функции 8 (t) и моментами ординат случайных функ­ ций а (г) и ß (t).

Находя математическое ожидание обеих частей равенства (251),

имеем

 

8 (0 = *(0 + Ф(0.

(4.279)

т. е. A(t) есть вещественная, а ß(t)—-мнимая части комплексного выражения 8 (t):

ä(t) = Reb(t), ß (t) = Im 8 (t).

(4.280)

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

ведения [8*(^) — 8* (іг)] [8 (£2) — 8 (£2)], как это делается в соответ­ ствии с формулой (1.62) при определении корреляционной функ­ ции K s (tj, t 2) комплексной случайной функции, или математическое

ожидание произведения [8 (£х) — 8 (^)] [8 (t2) 8 (£,)], которое можно рассматривать как взаимную корреляционную функцию 8* (іг) и 8 (t.2),

т. е. как Rm{tlt £2). Выражая в обоих этих произведениях 8 (t)

224 ГУ, ОПИСЫВАЕМЫЕ ЛИНЕЙНЫ М И УРАВНЕНИЯМ И ІГЛ. 4

и §(£)

через

а(І), ß (£), &(t)

и ß (t)

по формулам (251)

и (279),

получим

 

 

 

 

 

 

tt)=M ({[«(*!)—«(*!)]—

 

 

 

-

I [р (4) -

ß С^)]> {[«(У -

а (*,)] +

aß (а) -

ß (а)]}).

(4-281)

# 8 * 3 ( 4 , І 2) = М ( { [ « ( І 2) — « ( 4 ) J +

 

 

 

+

i [ß (а) -

ß (а)]} {[* ( g -

* (ія)] + 1tß ( g -

ß (g j» .

(4.282)

Раскрывая в последних выражениях фигурные скобки под знаком математического ожидания и находя математическое ожидание суммы с учетом формул (1.62) и (1.63), получим

# s(4 , t2) = K a(tv g

+

^ p (g g +

*#«P(4. g

— # g ( 4 >g> (4.283)

д»**(4>tt) = K a(tlt g

Kß(tlt g

+

i#«ß(g g

+

«#РЛ4 >g -

(4.284)

Отделяя в последних равенствах

вещественную и

мнимую части,

будем иметь

 

 

 

 

^ ѳ # з ( 4 > g> )

 

 

( 4 > 4 ) 4 ~ # ß ( 4 >

4 ) =

 

 

# a(4> а)

(4> ^2) —

#5*s (4>

g > }

 

 

# 0ß(4i

а)

Rp*(ßu

4) — i m #s(4>

g>

i

/4 2яні

#«ß(4. g +

#p«(g

g

=

im #*.8(«1,

g .

j

i

i

Произведя сложение и вычитание каждой пары уравнений системы (285) и сложение уравнений (286), получим окончатель­ ные выражения для корреляционных функций Ка (tv t2), К» (tlt t2) и взаимной корреляционной функции # я„ (£х, г2):

# а (4>

а )— 2

(4> а)4-#5*з(а> а)]’

#ß (4>

g = l R e [ t f s (tu

а)

# 8*8 (4>

д]>

#aß (4>

g =

-f Im [ ^ 8

(а,

g

+ # 8*8 (^i!

а)]-

Введем обозначения

U (£),

Im а (£) == V (t),

 

Re а (£) =

 

n e X ( t ) = X 1{t),

Im X(t) =

X a(t),

 

t

 

t

 

 

t

 

 

cp (£) = ^ а (ij) dfj = J U (fx) dt1+

i ^

F (£x) d£x =

o

o

 

 

 

о

 

 

 

 

 

 

= <Pi(*) + *<Ps (*)•

Тогда формулу (277) можно переписать в виде

 

 

 

t

 

 

 

 

 

8 (f) =

8 (0 )<r*(O-j- j

 

 

 

 

 

(4.287)

(4.288)

(4.289)

(4.290)

(4.291)

0

§ 4.4І

СИСТЕМА ДВУХ УРАВНЕНИЙ 1-ГО ПОРЯДКА

225

Дальнейшее рассмотрение целесообразно провести отдельно для варианта, когда коэффициент а (t) уравнения (252), а следо­ вательно, и функция <р(t), определяемая равенством (290), яв­ ляются детерминированными (неслучайными) функциями времени

идля другого варианта, когда а (t) — случайная функция.

Впервом случае, находя математическое ожидание обеих ча­

стей равенства (291), получим

t

8 (t) = 8 (0) е—(й -(- ^ е-?Н)+?(Ц)х (^) dtv (4.292)

о

Следовательно, отделяя вещественную и мнимую части в по­ следнем равенстве, учитывая (280) и обозначения (290), имеем

а =

е-?> [а (0 ) cos <р2 (*) -f ß (0 ) sin ср2 (*)] +

 

 

 

І

 

 

 

+

j е"*«)+?.«,)

(^ )cos[tp2 (*x) <f»s (f)] —

 

 

 

0

 

 

(ij) sin [cp2 (*x) — «Pa («)]} * 1.

(4.293)

В== е-и (O [—ä (0) sin cp2 (f) +

ß (0) cos cp2 (i)] +

 

 

+

j e m W +ъМ {Xl (k) sin p2(k) «p2 (0] +

 

 

 

0

 

 

 

 

+

*s (*l) cos [?2 (*l) — ?2 (0 1 } d t v

 

Находя корреляционную функцию для функции Ъ(і) [см. (281)] (будем считать начальные значения a (t) и ß (t) некоррелирован­ ными друг с другом и с ординатами функции X (t)), получим

к , {h, к) = (D [а (0)] + D [ß (0)]} e-T*(K)+9(W +

/ 2 Л

+ j j е-^(5)+т*(ч)-р(Щ +т(ч)^(Хі> г2)йт^т2. (4.294)

о о

Аналогичным образом вычисление Rs*s(t1, t2) дает

Rm (tv t2) = (D [Д (0 )] — D[ß (0)]}

+

 

h и

 

 

+ 5

I е-'Р(<')+<Р(и)-т(«+т(ч)дЛ

(Ті) T

(4.295)

0

0

 

 

где через Rx*x(t\* t2) обозначено

 

 

Rx*x (k,

t2) = M {[X (k) 3- ( t , ) ] [X (к) x (k)]}.

(4.296)

15 А, А. Свешников, С. С. Ривкин

226

ГУ, ОПИСЫВАЕМЫЕ ЛИНЕЙНЫ М И УРАВНЕНИЯМ И

[ГЛ. 4

Формулы (294) и (296), совместно с формулами (287), (288), (289), определяют корреляционные функции и взаимную корреля­ ционную функцию случайных функций а {£) и ß (t), характеризую­ щих поведение ГУ в данном случае. Входящие в эти формулы кор­ реляционная функция Кх (tv t2) и взаимная корреляционная функ­ ция R x*x (tv t2) в соответствии с обозначениями (290) могут быть выражены через корреляционные функции KXl (tv t2) и КХг (fx, i2) вещественной и мнимой части правой части уравнения (252) и взаимную корреляционную функцию RX]Xi (tu t2). Выполнив про­ стые преобразования, имеем

^ж(^І>

^2) == К х, (t j , i2) -f-/ѵЖі(і], £2) -f-

 

 

 

+ і[Д *л(гі> t2) — Rx^ { t v

<2)],

(4.297)

Rx*x{t\i

t‘>) = K.Ti (£ц £.,)— Kx2 (i], £,)+

 

 

 

+ i\Rxtx2(£и ^,) + ^x2x,(^|,

£•,)!•

(4.298)

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

Положение меняется, когда коэффициент а (t) является случай­ ной функцией времени.

Формулы (287), (288), (289) справедливы и в этом случае; также справедливой остается и формула (291), дающая общее решение уравнения (252), однако связь между моментами 8 (t) и моментами ординат функций X (£) и а (t) в этом случае значительно услож­ няется. Действительно, так как функция <р(t), определяемая фор­ мулами (290), теперь является случайной, то из формулы (291) уже не следуют формулы (292), (294) и (295), вычисление корреля­ ционных функций K s (tv t2) и i?8*8 (tv t2) в общем случае значи­ тельно усложняется и простые результаты удается получить только, когда функции U (t) я V (t) являются нормальными.

Рассмотрим этот случай подробнее.

Так как в соответствии с (290) ^ (<) и <р2 (t) являются интегра­ лами от нормальных случайных функций U (t) и V (t), то они также являются нормальными, а закон распределения их ординат полностью определяется их математическими ожиданиями, корре­ ляционными функциями и взаимной корреляционной функцией, т. е. величинами, связанными в соответствии с (290), (1.69) и

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