
книги из ГПНТБ / Свешников А.А. Вероятностные методы в прикладной теории гироскопов
.pdf
388 Н Е Л И Н Е Й Н Ы Е УРАВНЕНИЯ ПРИКЛАДНОЙ ГИРОСКОПИИ [ГЛ. 6
Учитывая вышеприведенные соотношения и формулу (109), уравнение (106) можно представить в виде
дЕ . . |
|
> дЕ . , |
|
. <?£■ - |
|
|
|
||
dz |
"Ь ( ! Ѵ і — |
Р Л ) ~äz----- Ь ( V i + |
vi z 2 — z3> ~ лГ + |
|
|
|
|||
|
|
|
+ (c22 2 + |
ftez4) |
+ (pTz5 — k2zx) |
+ |
|
|
|
|
|
+ (teiȀz2 + |
°і!Ѵ4 + |
Ф / 5) E + |
^0 = ^ + 1 + |
|
|
||
■ c 2 |
dE |
|
J |
E(zv z2, z3 + z, z4, z6,ze)-^- = |
0, |
(6 .1 1 0 ) |
|||
|
2 |
dz. |
|
||||||
|
|
|
|
|
|
|
|
||
где |
dz. |
обозначает |
частную производную |
от |
в |
которой |
|||
аргумент z6 заменен на z6 + |
1 соответственно. |
|
|
|
|||||
|
Дифференцируя |
это |
равенство по различным переменным z . |
и полагая потом все переменные z^.=0 , получим ряд рекуррентных соотношений, связывающих последовательные моменты случайной
величины и г (т)=а (т) (считаем, |
что все |
появляющиеся |
таким об |
|
разом моменты существуют). |
|
|
одному разу по |
|
Так, продифференцировав|последовательно по |
||||
z2 и z3 и положив потом %. = 0 , получим два уравнения, |
содержа |
|||
щие â (т) и В (т), решая которые, |
находим |
|
|
|
t |
|
___ |
|
|
а (t) = сх I е_|1‘ <і!_т)п;к(т) dz, |
ß(£) = |
0 . |
(6 .1 1 1 ) |
|
о |
|
|
|
|
Дифференцируя (110) по другим аргументам характеристической функции Zj и комбинируя получающиеся после обращения в нуль всех переменных Zj уравнения, можно получить ряд соотношений между моментами случайных величин U . (т) более высокого порядка.
При исследовании ГИ обычно представляет интерес не только ошибка ускорения, пропорционального â (г), но и ошибка опреде
ления самой скорости, пропорциональная отклонению угла а (t) t
от значения — ^ wK(^) dtv которое должен показывать интегратор
о
в том случае, когда ошибки отсутствуют. Определение ординат функции а (t) можно выполнить аналогичным образом, используя уравнения (1 0 0 ), однако расчеты при этом усложняются.
Исследование ГИ, приведенное выше в качестве иллюстрации применения теории марковских процессов к анализу нелинейных гироскопических устройств, может быть проведено и каким-либо другим приближенным методом, например, методом последователь ных приближений.
390 Н ЕЛ И Н Е Й Н Ы Е УРАВНЕНИЯ ПРИКЛАДНОЙ ГИРОСКОПИИ [ГЛ. 6
Учитывая, что значение t предполагается большим, в послед нем выражении сохраним только слагаемое, пропорциональное t, а верхний предел интегрирования положим равным бесконечности. После этого получим
СО |
|
|
D[a(t)]*»f S cos ^ [ Я ,,( т ) + |
£ ,,(*)]* . |
(6.115) |
о |
|
|
Для определения корреляционных функций |
|
|
Zx (t) = т)2 sign 6 (t) и Z2(t) = |
тд sign 9 (f) |
|
в соответствии с (4. 38) имеем
АД (т) = |
2 |
\ |
— т]| arc sin kb(т), |
(6.116) |
|
|
|
|
АД (Д = |
a r o s i l 1 h i '1)’ |
, |
где нормированные корреляционные функции кь (") и к? (т) в соот ветствии с (1 1 2 ) определяются равенствами
кь (т) = е~^е1т 1(c o s Хѳт ------ |
s j n Х8 1т |^ , |
(6.117)
/сф(т) = е '^ І т| (1 — [а |т|).
Дальнейшие вычисления приходится вести численно. Подстав ляя (116) в (115), учитывая (4.36) и производя численное инте грирование, получим
D[a(ü)] = 5 • IO“ 10 t рад2!сек. |
(6.118) |
Например, при t= l0 мин D [а (t) ]= 3 -ІО- 7 рад2.
Для сравнения рассмотрим решение этой задачи методом ста тистической линеаризации. Согласно этому методу [25] нелиней
ные выражения sign Ѳи sign сЬнужно заменить |
линейными функ |
циями |
(6.119) |
sign Q= kfl, sign cp = /Cjtp, |
где коэффициенты k2 и А;'подобраны так, чтобы дисперсии обеих ча стей равенства (119) были одинаковы, т. е. в нашем случае, учиты вая (117), имеем
(6.120)
_!/" |
і _ |
1 |
• |
У |
АГФ(0) |
|
Подставляя (119) в (4), для определения a(t) и ß (t) получим си стему линейных уравнений с постоянными коэффициентами, в пра
§ 6.4] |
ПРИМЕРЫ НА НЕПРИВОДИМЫЕ Н ЕЛ И Н Е Й Н Ы Е УРАВНЕНИ Я |
393 |
|||||||||||
ajj = |
0 , 6 - 1 0 рад2; |
= |
0,042 |
1 /сек; |
\ = |
0,42 |
1 /сек; |
а° = |
30°; |
||||
х = S/H = 0,05 |
1 /сек; |
х' = |
xz/g = |
0,051 |
сек [см. (3.38)]. Начальные |
||||||||
значения ошибок а (f) |
и ß (£) |
равны нулю. |
|
|
|
|
|||||||
Р е ш е н и е . |
|
Система |
уравнений |
для |
рассматриваемого |
ГН |
|||||||
имеет вид (3. 37), |
т. е. |
|
|
|
|
|
|
|
|
|
|||
|
|
|
d — ßO cos (a0 -f- а) ==0 , |
|
|
(6.132) |
|||||||
|
ß + |
xß = |
— (Ѳ—f- хѲ —j- х'Ѳ) sin (а0 -j- а). |
||||||||||
|
|
|
|||||||||||
В нулевом приближении в соответствии с (62) имеем |
|
|
|||||||||||
|
d0 — ß0Öcos а0 = |
0 ,' |
|
|
|
|
|
(6.133) |
|||||
|
|
|
ß0 -f- xß0 = |
— (0 -f- хѲ -f- х'Ѳ) sin а0. |
|||||||||
|
|
|
|
|
|||||||||
Из второго уравнения находим |
|
|
|
|
|
|
|||||||
|
ß0(t) = |
|
|
і |
|
|
|
|
xö(x)-j-x'Ö(x)]dz, |
|
|
||
|
— sin a° ^e~x |
|
[0(x)+ |
(6.134) |
|||||||||
T. e. |
|
|
|
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Po ( 0 |
= |
o . |
|
|
|
|
|
|
|
M [ß0 i t ) 6 (0 ] = Sin a° |
t |
|
|
(x) - |
х £ ѳ (x) - |
x 'Ä 8 (x)] d z . |
|
|||||
Äßoé = |
j e - " [ K , |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
(6.135) |
Интегрируя по частям и ограничиваясь рассмотрением установив
шегося процесса, |
получим |
|
|
|
||
|
*М = |
~ |
№ + Х?) ^ |
+ 2^ e l |
(6-136) |
|
Находя математическое ожидание первого равенства системы |
||||||
(133), имеем |
|
|
|
|
|
|
а0 = kßj cos а° = |
— 2 T(- ^ j 3 + Щ W + *І) W + |
|
|
|||
и |
|
|
+ Xg -f- 2 х(гѳ) sin2 a° |
(6.137) |
||
|
|
|
|
|
|
|
ао — аог — |
2 f(x + рѳр + X[5] (^ 0 + ^ |
(^ 0 + x 0 + |
|
|
||
|
|
|
+ 2xfi,) sin 2а° = |
—0,23 • 10~4t |
1/сек. |
(6.138) |
Поправка |
(t) |
к |
математическому |
ожиданию |
В(t), найденному |
в первом приближении согласно (71), равна нулю. Следующие же приближения дают в данном случае несущественные поправки. Аналогичным образом можно убедиться, что â0 (t) в данном примере дает математическое ожидание d (t) с достаточной для практики точностью-
394 |
Н Е Л И Н Е Й Н Ы Е УРАВНЕНИ Я ПРИКЛАДНОЙ ГИРОСКОПИИ |
[ГЛ. 6 |
§ 6.5. Гироскопические устройства, при исследовании которых возникают нелинейные задачи, не связанные
свидом уравнений ГУ
Вприкладной теории гироскопов возникает ряд нелинейных задач, не связанных с видом уравнений ГУ. К подобным задачам можно отнести исследование кардановых погрешностей гироскопа направления [23], [82], определение взаимного азимутального по ворота стабилизированных площадок на корабле при различной ориентации осей их кардановых подвесов [23], определение сред него времени пребывания ошибки ГУ в данном интервале и др.
Рассмотрим несколько примеров на решение подобных задач. Пример 6.4. Найти математическое ожидание и дисперсию кардановой погрешности ГН, установленного на самолете и слу жащего для определения угла рыскания ф самолета. Ось гироскопа (при отсутствии наклонов основания) составляет с продольной
осью самолета угол х.
Дано: углы крена у(£), тангажа &(t) и рыскания ф(і) само лета — независимые нормальные стационарные_ случайные функ ции времени, их математические ожидания у=& = ф = 0 , а средние квадратические отклонения 0^= oft= о = 1 °.
Р е ш е н и е . Карданова погрешность Дф ГН определяется соотношением Р5]
|
Дф = х + ф — arctg[tg(x + ф) Щ Ь + tgO siny]. |
(6.139) |
||||||
При |
«продольной» |
установке ГН, |
когда х = |
0, имеем |
|
|
||
|
Дфх=0 = Ф— arctg (tg ф |
+ tg » sin у ). |
(6.140) |
|||||
При |
«поперечной» |
установке ГН, |
когда х = |
90°, |
из (139) |
сле |
||
дует |
|
|
|
|
|
|
|
|
|
л I |
t |
. |
t g Ф COS Ь |
.— . |
/ „ |
. , . . |
|
|
Дфх=90о = |
ф — arctg ---------т-Ц—. о |
(6.141) |
|||||
|
т |
т |
6 cos у — tg ф Sin 0 SH1 у |
|
' |
' |
||
Считая углы ф, 9 |
и j малыми (поскольку их средние |
квад |
||||||
ратические отклонения равны |
1°), вместо (140) |
и |
(141) |
прибли |
||||
женно имеем |
Дфх=0 = |
—&у, |
|
|
|
(6.142) |
||
|
|
|
|
|
||||
|
Дфх_90о = |
_ ± ф ( т2_& 2). |
|
|
(6.143) |
|||
Учитывая независимость случайных величин & и у, |
из |
формулы |
||||||
(142) следует |
|
|
|
|
|
(6.144) |
||
|
|
Дфх=о = —&у = 0 , |
|
|
||||
|
|
°itx=0 — 0»°у |
|
|
|
(6.145) |
8.53 |
ПРИМЕРЫ ДРУГИХ |
Н Е Л И Н Е Й Н Ы Х ЗАДАЙ |
|
|
|
395 |
||||||
Подставляя числовые данные примера, получим |
|
|
|
|
|
|||||||
|
сДфх=0 = |
0,3045 • Ю~ 3 рад = 147. |
|
|
|
|
|
|||||
Аналогичным образом, учитывая отсутствие зависимости |
ф от |
|||||||||||
у и & и принимая |
во внимание, |
что ф = 0 , согласно (143) |
получим |
|||||||||
|
Дф.-0П° = |
- у ?(D [tJ - |
D ІА]} = |
0 , |
|
|
(6.146) |
|||||
°ІФх-воо = |
Т °Ф[ 3 К |
+ |
— 2зК і = |
|
|
|
|
|
|
|
||
|
|
|
|
= т - [ 2 Н + 1 |
) + |
(а? - 1 |
)2]- (6-147) |
|||||
По формуле (147) вычисляем а Дфх=90О: |
|
|
|
|
|
|
||||||
|
адфх=90о = 0,531 |
- ІО- 5 |
рад = |
1",1. |
|
|
|
|
||||
Сопоставляя полученные |
выше |
значения, замечаем, |
что аДфх=доо |
|||||||||
• ^ адфх=0! |
это обстоятельство указывает на целесообразность «по |
|||||||||||
перечной» установки ГН на самолете. |
|
угла |
|
Представим |
||||||||
Найдем а Дф для |
произвольного значения |
х . |
||||||||||
для этой цели формулу (139) приближенно |
в виде |
|
|
|
|
|||||||
|
|
|
Дф = —&у cos2 (х -f- ф), |
|
|
|
|
(6.148) |
||||
откуда при малом ф находим |
|
|
|
|
|
|
|
|
||||
|
Дф = |
—&у (cos2 X— ф sin 2х — ф2 cos 2х). |
|
(6.149) |
||||||||
Так как § = у = |
0, то из (149) получим |
|
|
|
|
|
|
|||||
|
|
|
|
Дф = 0. |
|
|
|
|
|
(6.150) |
||
Для дисперсии |
D [Дф] имеем |
|
|
|
|
|
|
|
||||
|
D [Дф] = М {[Дф]2} — [Дф] 2 = М {[Дф]2)- |
|
|
(6.151) |
||||||||
Подставляя сюда (149), |
находим |
|
|
|
|
|
|
|
||||
D [Дф] = |
D [&] D [у] {cos4 X4- D [ф] sin2 2х -f- |
|
|
|
|
|
|
|||||
или |
|
+ |
3 (D [ф] ) 2 cos2 2х — 2D [ф] cos2 х cos 2х} |
(6.152) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
аДф = а|а2 [cos4 х + |
а | (sin2 2х — 2 cos2 х cos[2x) -[- Заф cos2 2х]. |
(6.153) |
||||||||||
Для аф = а&= |
|
= 1 ° и различных величин установочного угла |
||||||||||
X по формуле (153) вычислены |
значения адф и зависимость |
аДф = |
||||||||||
= / (х) построена на |
рис. 6.3, из которого следует, |
что |
аДф |
дос |
||||||||
тигает наименьшего |
значения при х = |
90°. |
|
|
|
|
|
|
||||
Пренебрегая в (153) тремя малыми последними слагаемыми, |
||||||||||||
приближенно имеем |
|
|
o2cos4 X, |
|
|
|
|
(6.154) |
||||
|
|
|
|
а ІФ = |
|
|
|
|
396 Н Е Л И Н Е Й Н Ы Е У РАВНЕНИ Я ПРИКЛАДНОЙ ГИРОСКОПИИ [ГЛ. 6
откуда следует, что при х = 90° в первом приближении
Одф = 0.
Рассмотрим, насколько изменится полученный результат, если учесть отброшенные в (153) слагаемые. Вводя обозначения
|
COS 2 У . ~ Х , |
а д ф = < р ( ж ) , а§ат : |
(6.155) |
|
перепишем (153) |
в виде |
|
|
|
ср (х) — к |ж2 |
+ |
(1 — х) х — -тгх(2х-— 1 ) -|- Ззф (2х — 1 )2J, |
||
откуда |
|
|
|
|
|
|
= 2 к [X (1 - 8 4 + 1 2 4 ) + 3 4 — 6 а |] . |
(6 .1 5 6 ) |
|
Так как 4 |
<С о2<^1, |
то при любых значениях х |
в интер |
вале (0,1) dtp (x)/dx )> 0. Это означает, что функция tp (х) не имеет
в |
этом интервале экстремума, |
а |
наименьшее значение функции |
tp (X) будет при наименьшем зна чении аргумента, т. е. при х=0
или при х=90° [см. (155)]. Сле довательно, оптимальным значе нием установочного угла является х=90°.
Наиболее полное устранение кардановой погрешности ГН может быть достигнуто путем установки его на стабилизированной относи тельно плоскости горизонта пло щадке, которая может представ
лять собой площадку силового гирогоризонта или площадку, стабилизируемую от центральной гировертикали с помощью сле дящей системы.
Пример 6.5. Определить вероятность того, что нормаль к ста билизированной площадке корабельного гирогоризонта силового типа в случайный момент времени t не отклонится от истинной вертикали места больше чем на угол у0, обеспечивающий условия нормальной работы установленного на площадке прибора (секстан, гравиметр).
Дано: срединные значения *) ошибок стабилизации площадки по осям подвеса (рис. 6.4) Ел= 3'; Е^=3'; математические ожида ния этих ошибок <х= р = 0 ; ось вращения наружного кольца кар-
*) Для нормального закона распределения срединное значение Е ошибки связано со средним квадратическим значением ее а соотношением Е —
= р Ѵ/2'а(р ==0,4769...).