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

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

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

§ 4.2]

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

187

жим, что угол отклонения физического маятника і (t) имеет кор­ реляционную функцию вида

Кг (х) = a^éf1,7.|т| (cos Х^х + р sin \ г I XIj .

(4.98)

Подстановка выражения (98) в (97) позволяет проследить влия­ ние Xна погрешности ГВ, вызванные погрешностями физического маятника. Сделав указанную подстановку и выполнив интегриро­ вание, найдем для первого слагаемого в формуле (97)

СО

(4.99)

О

Полученная формула показывает, что коэффициент пропорцио-

 

 

 

 

V 2

 

 

 

 

 

 

нальности к —

 

1 + 7*7

 

у

дисперсии

погрешности

 

 

(_1+^+§) +4 (і?)

безразмерных параметров:

физического маятника зависит от двух

х/рх =

ж и 4х/рх =

у.

При значениях

 

 

 

 

 

 

 

 

 

<

| [ - 7 2 + 1 +

2Ѵ1 +

У2+ г/4]

 

 

этот

коэффициент растет, при (—

 

у2 +

1 -)- 2

\ / 1 +

г/ 2 + г/4]

достигает

максимума, а затем монотонно убывает, стремясь к нулю

при х /^ -*

со. На рис. 4.1, а и б представлена зависимость коэф­

фициента к от x/jj-x =

z при различных значениях ^х/пх — У- Анализ

кривых рис. 4.1

позволяет выбрать

значение

удельной

скорости

коррекции X (или постоянной времени

Т =

1/х), если

считать, что

основное влияние на дисперсию погрешности

гировертикали ока­

зывает первое слагаемое в формуле

(97).

 

 

 

 

Для вычисления второго слагаемого в этой формуле его целесо­

образно

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

по

частям преобразовать к виду

^ 2 \

arcsin &в(т)гіт

Qi

I

m

 

■fcj(x)dx. (4.100)

У-2 / / 2

'

KtfiҢЪ 5

О

 

 

 

 

О

\/l + k\ (t)

Последний интеграл удобнее вычисляется численным интегри­ рованием, чем интеграл (97). Его значение зависит от вида кор­ реляционной функции Ы (т). Если корреляционная функция угла Ѳ(t) имеет вид (2.13), то

(4.101)

кд(х) = е ’‘'в**1(c o s Xöt — р sin Хѳ | х .

188

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

[ Г Л . 4

6)

Рис. 4.1. Графики зависимости коэффициента к от величины х = х /ц 7 при раз­ личных значениях у = Ру

Рис. 4.2. Графики зависимости 1 -\-F (х, у) от величины х х/цѳ при раз­ личных значениях у=\д/рд.

§ 4.2]

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

189

Введя в этом случае в (100) новую переменную интегрирова­

ния I =

(іѳт, замечаем, что интеграл зависит

от

двух безразмер­

ных параметров: х/р0

и Хѳ/[х8, т.

е.

 

 

 

 

 

2Ql

 

 

=

J f s f 1 +

^ ( ^

, ^ ) ] .

(4-102)

 

ТІХ. / / 2 ] в-**arcsin к„ (т) *

где функция

 

 

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

f ( ~ ,

Ре/

71 J

 

#ѳ (т) dz

 

(4.103)

 

\Рѳ

 

 

 

 

 

может быть найдена численно.

 

 

 

-f-.F (х/р„, Хѳ/р9)

Для примера на рис. 4.2 приведены графики 1

для нескольких значений l fJp.b=y

как функции х/у^—х. Графики

показывают, что выражение

(1 0 2 )

довольно сильно зависит от х.

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

выборе удельной скорости коррекции нужно учитывать

зависи­

мость от X обеих слагаемых в формуле (97), дающей значение дис­

персии погрешности ГВ.

 

3. Гиротахометр. Для гиротахометра в соответствии с (76)

X(t) = *(t) + ±-M,

(4.104)

а а1= х = і / Т , где

ш (t)

— измеряемая

угловая скорость.

имеющей

Будем считать

со (t)

случайной

функцией времени,

математическое ожидание ш (t)

и

корреляционную

функцию

Кш(т). (Случайный характер угловой скорости может быть выз­ ван, например, случайными перемещениями объекта, по которому производится визирование.)

Задачей гиротахометра является определение текущего зна­ чения угловой скорости. Следовательно, ошибкой прибора в дан­

ном случае будет не угол (3(f), а разность

 

-«>(*)■

(4.105)

Подставляя в последнее равенство вместо

ß (t) выражение (80),

с учетом (81) и (104) получим формулу для ошибки гиротахометра

t t

 

t

t-т

 

T <0(х) <^ - —

+

5 е

T M(z)dz,

(4.106)

о

 

о

 

 

где возмущающий момент М (t) будем считать стационарной слу­ чайной функцией; также стационарной функцией будем считать разность [ м (t)—<со (£)].

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

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

 

 

 

*

_ч_

 

 

 

 

 

Z(t) = Y ^ e

T&(t — x1 )dx1

m(t).

(4.107)

 

 

 

о

 

 

 

 

 

Предположил!, как ото обычно и бывает в действительности,

что математическое ожидание <г>(t) достаточно плавно

меняется

со временем, так что в интервале т,,

соизмеримом с постоянной

времени Т, функцию о>(t—тх) можно

разложить в ряд Тейлора

около точки t, т. е. написать

 

 

 

 

(О(t -- t:) =

(t) --- di»(t)

J_

( 0

T 2

(4.108)

 

 

 

 

~~dt~

2

dt*

1

 

Подставляя

(108) в (107) и интегрируя по т1;

получим

 

£ (t) ~ — е

— Т

1

0+f>

dto (£)

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ Т1

 

 

 

 

d2o) (г)

(4.109)

 

 

< * + т - & У '

 

dt2 +

Если t T (переходный процесс закончился), то последняя формула упрощается и принимает вид

£(г) = _ г ^ і і ) +

Т2 ^ М -------

(4.110)

О С

CL L

 

Формулы (109) и (110) показывают, что при переменной угловой скорости ГТ дает систематическую ошибку, зависящую от угло­ вого ускорения dw/dt и производных от угловой скорости ш более высоких порядков, а также от постоянной времени Т гиротахо­ метра, уменьшаясь с уменьшением Т.

Предположим, что момент М (t) обусловлен наличием сухого трения в осях подвеса. Тогда, находя дисперсию обеих частей равенства (106), учитывая (38) и рассуждая так же, как при выводе формулы (97), для установившегося процесса получим

 

 

 

со

Т

 

 

+

-jrpy-

^ е

r arcsin /Cg (х) dx.

(4.111)

СО

Т

 

о

 

 

D (()J = Y j Л

[А ш (0) -

Кю(т) 1

dx +

 

 

0

 

 

 

 

 

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

§ 4.2]

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

191

скорости и связанную с тем, что вследствие инерционности, при­ сущей ГТ, это изменение отрабатывается с искажением. Второе слагаемое дает ошибку, вызываемую моментом сил сухого трения. Формула (111) показывает, что дисперсия ошибки измерения ско­ рости зависит от постоянной времени ГТ сложным образом, по­ скольку Т входит кроме множителей перед интегралами в подын­ тегральные выражения. Поэтому вывод о целесообразности умень­ шения Т для снижения систематической ошибки, полученный на основании анализа формулы (109), должен быть уточнен путем исследования зависимости D [е (£)] от Т.

4. Поплавковый интегрирующий гироскоп. В качестве по­ следнего примера дифференциального уравнения первого порядка, содержащего зависимую переменную, рассмотрим уравнение (78) для поплавкового интегрирующего гироскопа. В этом случае

t

+ f j M ^ d t , ,

(4.112)

о

т. е. в правой части равенства стоит линейная комбинация стацио нарной случайной функции if (t) — угла рыскания основания гиро-

t

скопа и интеграла ^ М (fj) dt1 от возмущающего'момента, который

о

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

Находя математическое ожидание обеих частей (112) и учиты­ вая (82), получим

t т

t

т(t~т

1

 

р (t) = Y j е 1у (t — х) dx -f -L j e T M m (tj) dxJdx.

(4.113)

>o

о

>

 

Если if (t) и M (t) стационарны, то интегрирование в последней формуле может быть выполнено и мы получим

t + pL(^— т) + Те

(4.114)

или, для установившегося процесса,

р (t) = fof -f- р (t Т) tri.

(4.115)

Таким образом, наличие интеграла в правой части выражения (1 1 2 ) приводит к возникновению накапливающейся систематиче­ ской ошибки, если математическое ожидание возмущающего мо­ мента отлично от нуля.

ff)2 ГУ, ОПИСЫВАЕМЫЕ ЛИНЕЙНЫ М И УРАВНЕИИИМИ (ГЛ. 4

Для дисперсии (3 (t)

в соответствии с (84), (87) и (112) имеем

D|ß(f)] = * j

'2 £г—XTК(t

— х) dt -f-

 

 

 

I

I ІІ+Т.Д

 

 

 

 

 

+

r K»(t ~~xv

* — Ts) <M X*>

(4.116)

 

О

О

 

 

 

 

где через Kn(tv

t2) обозначена корреляционная функция интеграла

 

 

 

t

 

 

 

 

 

N {t)= \i M(tx) dtv

(4.117)

 

 

О

 

 

 

В соответствии с формулой (1.88) имеем

 

t%

 

 

tx

Кт(Х) d x ~

 

Кп(*П *2) = S (*s Х) Кт (Х) dx +

S

(*1 х)

 

О

 

 

О

 

 

 

 

 

и-іі

 

 

(4.118)

 

 

S

(*s — *і х )Кт(х) ^х-

 

 

 

о

 

 

 

Подстановка последнего выражения в (116) после простых пре­ образований дает

D[ß(f)]=

т 1

1^ 1

2

/гг

I&1

*;(*-T)dT+p2{[

 

О

 

 

 

 

 

і

/

і \ X

27’(l — е

Т)_

Те Т — Те

( - 2

+ е- г ) е Ц к т(г)

 

 

 

 

 

 

 

К,п(Х) dx-

 

 

 

 

 

 

 

(4.119)

Формула (119)

показывает,

что из-за наличия в правой части

 

 

 

 

 

 

 

t

уравнения

(78)

нестационарного слагаемого

р ^ М (tj) dtx диспер-

 

ß (t) в

 

 

 

 

 

о

сия угла

этом случае

будет расти со временем примерно

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

D[ß(i)l = e* + b,

(4.120)

§ 4.2]

 

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

193

где

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

 

а =

2 р2 J

(х) dx,

 

 

О

 

(4.121)

 

 

 

■ 2p*J (jc + г + 1 е- * ] Я я (х )< гх + -£ Je“ * /^(xjdx.

Определение закона распределения случайной функции я (t), являющейся решением уравнения (79), осуществляется просто только, если правая часть этого уравнения X (t) является нормаль­ ной функцией. В этом случае я (t) будет также нормальной, и для определения законов распределения этой функции любого числа измерений достаточно знать только я (t) и Ка(£l5 12), формулы для которых получены выше. В том случае, когда X (t) не является нормальной, так же как и для уравнения (1 ), закон распределения ординат функции я (t) может быть выражен практически с любой точностью через моменты этих ординат более высокого порядка. Однако нахождение этих моментов усложняется тем, что для их вычисления нужно располагать моментами соответствующих по­ рядков случайной функции X (t), которые для закона распреде­ ления, не являющегося нормальным, уже не могут быть выражены через Кх (tb t2). Таким образом, задача определения закона рас­ пределения я (t) уже выходит за рамки корреляционной теории случайных функций.

Формулы, связывающие моменты ординат я (t) более высокого порядка с моментами ординат случайной функции X (t), могут быть получены тем же методом, что и формулы для ä(t) и Ка (tv t2). Так, например, для третьего и четвертого начальных моментов будем иметь

t

i

t

 

mW = ^

J

J e-"* (n+4 +-t3) X

 

ooo

 

 

 

 

X M [X (t — Tj) X (t t2) X (t t3)] dxjdxgdxg,

(4.122)

t

t

t t

 

mW ~ J

J J ^ e-Oi (ч+4+ѵи,) X

 

o o o o

 

X M [X (t — Xj) X (t — x2) X (t — x3) X (t — x4)] dxjdxjdxj^.

(4.123)

Для моментов более высокого порядка эти формулы еще более усложняются.

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

194

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

[ГЛ. 4

§4.3. Линейное уравнение второго порядка

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

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

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

ä -f- ajâ -f- a2a = X (t),

(4.124)

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

рых постоянные ах и а2, а также и случайная

функция X (t)

имеют различные значения.

служить уравнения

Примерами уравнения типа (124) могут

(3. 89), определяющие

угловые отклонения

а и

ß гиромаятника

с демпфированием

 

 

 

ä -f- 2xa -[- (х2 -f- к2) a =

-i- (х2 -f- к2) (тф — zÖ) +

(т)с — zÖ), (4.125)

Р + 2xß + (х2 + * 2) Р =

у (Ч</ - z'Ö),

 

(4.126)

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

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

X + 2^пХ + п2Х = K w (t),

(4.127)

где компонента ускорения w (t) в том случае, когда плоскость, качания маятника параллельна плоскости симметрии объекта, на котором установлен маятник (например, диаметральной пло­ скости корабля), определяется формулой [см. (3.56) ]

W (t) = Іо у? + ztj>,

(4.128)

а в том случае, когда плоскость качания перпендикулярна пло­ скости симметрии (например, лежит в плоскости шпангоута ко­ рабля), определяется формулой [см. (3.55)]

w (£) = % + щ — zÖ.

(4.129)

Уравнение (3.109), определяющее отклонение (3инерциальной вертикали при наличии демпфирования, принадлежит к тому же типу, так как имеет вид

(4.130)

Р + - ^ + *8Р = ! ^ + й-лЗг(<).

§ 4.3]

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

195

Аналогичным уравнением является и уравнение (3.112) для инерциальной вертикали при наличии внешней информации, поступающей с ошибками:

 

+

+

(4Л31>

где

bS (t) — ошибка счисления

пути движущегося

объекта,

а

Ъѵ (t) — ошибка определения

абсолютной скорости

объекта.

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

объекта определяется

с помощью допплеровской системы, то

Ъѵ (t) можно считать

стационарной функцией времени. Если

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

К этому же типу относятся и уравнения (3.122)

и (3.132)

для гиротахометра, имеющие вид

 

 

ß + + п2ß =

+ Рі (Мг + М)

(4.132)

и

ГѲ — pn^,

(4.133)

Гр + 2С Гр + ß == &ос +

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

выми частями

равенств.

 

 

 

 

К уравнению второго порядка с постоянными коэффициентами

сводится и уравнение гиротахоакселерометра (3.143)

Я2а +

all (ci^2 +

c2^l) ®+ с-^ЛЩа. = Нсг1

+ # 2ф (4.134)

и уравнение (3.149)

для

вибрационного гироскопа

 

 

Ä-)- 2Сиа м2а =

kjW', (t) sin Qt.

(4.135)

Последнее уравнение, не отличаясь по типу от уравнения (124),

имеет ту особенность,

что случайная функция

о> {{) входит в пра­

вую часть равенства

с коэффициентом sin

изменяющимся во

времени с заданной частотой £2 .

 

 

 

Наконец, к рассматриваемому типу относится и уравнение,

определяющее

интеркардинальную

погрешность

однороторного

гирокомпаса, которое для случая бортовой качки имеет вид [см.

систему (3.293)] (при

7 = ^ = 0 )

 

 

* 4- 2№х +

**« = -

Ѳ2 (t) sin 2К.

(4.136)

Частным случаем линейного уравнения второго порядка яв­ ляется уравнение, не содержащее первой производной от зависимой переменной, т. е. уравнение вида

й + я2а = х (*),

(4.137)

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

решение которого обладает рядом свойств, отличных от свойств

решения уравнения

(124).

характеризует поведение

ряда ГУ.

Уравнение типа

(137)

К этому типу принадлежит

уравнение (3.81) для гиромаятника

без демпфирования

й ~Г к'2а qkMi +

(4.138)

 

где Му и М 2 — возмущающие моменты по наружной и соответ­ ственно внутренней осям карданова подвеса.

Аналогичными являются уравнения (3.95) для отклонений инерциальной вертикали с периодом М. Шулера (без демпфиро­ вания)

й +

v2a — qM„,

(4.139)

ß +

^ = qMv

(4.140)

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

Наконец, к рассматриваемому типу относятся уравнения, определяющие ошибку инерциальной вертикали вследствие раз­ личных инструментальных ошибок системы. Для ошибок, вызы­ ваемых отклонениями кинетического момента Н гироскопа и коэф­ фициента у от расчетных их значений Н0 и у0 (условие невозмущаемости [i0/H0=l/R), уравнение (3.105) дает

 

Р+

 

й)*

 

(4-141)

где Н = Н 0-\-ЬН и p = f*o+V-

 

 

 

акселерометра

Для ошибок, вызываемых смещением нуля

ЬаѴі и дрейфом первого интегратора е,

в соответствии с уравнением

(3.107) имеем

 

 

 

 

 

P + v2ß = ^ ( 4 ,

+

s)-

 

(4-142)

Дрейф гироскопа в этом случае дает ошибку, определяемую

уравнением (3.108), т. е.

 

 

 

 

 

Р + ѵ2р = йд(г),

 

 

(4.143)

где угловая скорость дрейфа гироскопа

<і>д (t)

является случайной

функцией

времени.

 

 

числа

ГУ сводится

Итак,

исследование точности большого

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

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