
книги из ГПНТБ / Свешников А.А. Вероятностные методы в прикладной теории гироскопов
.pdf§ 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 |
о |
Iо |
> |
|
Если 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Щ + п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) |
является случайной |
|||
функцией |
времени. |
|
|
числа |
ГУ сводится |
Итак, |
исследование точности большого |
к анализу вероятностных свойств решения неоднородного линей ного уравнения второго порядка с постоянными коэффициентами,