
книги из ГПНТБ / Свешников А.А. Вероятностные методы в прикладной теории гироскопов
.pdfI 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) примет вид
•^г.эР — Hâ 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) и