Файл: Свешников А.А. Вероятностные методы в прикладной теории гироскопов.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 10.04.2024
Просмотров: 265
Скачиваний: 1
§ 4.1 І У Р А В Н Е Н И Е 1-ГО П О Р Я Д К А Б Е З ЗАВИСИМ ОЙ П Е Р Е М Е Н Н О Й |
179 |
являются еще и нормальными, поскольку это предположение обычно может быть принято. В этом случае система случайных функций wVj(t) и wz{t) полностью определяется корреляционными функциями KWy{x), KWz(x) и взаимной корреляционной функцией
Wz ( т)-
Находя математическое ожидание (64), получим
|
<7 Ся, с |
(с, СУі) |
wz (0 ). |
|
(4.65) |
|
У1 |
' |
|
|
|
Определение корреляционной функции с учетом (1.62) и (1.41) |
|||||
дает |
|
|
|
|
|
к Ah, h) = ^ ( c , |
Су)2ѴК«У^)КшЛХ) + |
|
Уі' |
(4.66) |
В том случае, когда имеют место все три рассмотренные выше причины, вызывающие уход гироскопа, в формулы (8 ) и (9) вместо f f i {£) и Кт(т) нужно подставить суммы [т(t)+y (t) ] и [ІІГт (т)+ + ЙГ (х)-{-Кг ( т) ] соответственно.
6 . |
Выводы из рассмотрения уравнения (1). Подведем некото |
|
рые итоги рассмотрения уравнения (1). Первые два момента |
||
случайной функции |
а(7), характеризующие ошибку ГУ, выража |
|
ются формулами (5), (6 ), (7), содержащими первые два момента |
||
случайных функций |
Н (t) и М (t). Моменты случайной функции |
|
Н (t) определяются |
конкретными причинами, вызывающими не |
|
стабильность кинетического момента, и не могут быть вычислены |
||
в общем |
виде. |
|
Учет |
неуравновешенности гироскопа может быть выполнен |
в рамках корреляционной теории. Для учета упругих деформаций потребовалось использование двумерных законов распределения (например, предположение о нормальном законе распределения при выводе формулы (6 6 )).
Моменты закона распределения сил жидкостного трения М ТХ определяются моментами закона распределения (того же порядка) угловой скорости вращения основания ГУ, а для определения моментов в случае сухого трения необходимо использовать пер вые два закона распределения угловой скорости. Однако и в этом случае вычисление т (і) и Кт (tu t2) не связано с принципиальными трудностями.
К уравнению первого порядка, не содержащему зависимую переменную, сводится и ряд других задач прикладной теории гироскопов: исследование влияния неаксиальности установки ротора [см. (2.88)] на ошибки ГН, исследование ошибки гиро интегратора [см. (3. 186) ] и ряд других задач.
Во всех этих случаях уравнение, описывающее движение гиро скопа, имеет такой же вид, что и уравнение (1 ), но только в его
12*
180 |
Г у , |
о п и с ы в а е м ы я л и н е й н ы м и |
у р а в н е н и я м и |
[Гл. 4 |
правой |
части |
стоит не момент трения, |
а случайные |
величины |
и функции, имеющие другую физическую природу. Исследование этих уравнений может быть выполнено так же, как и уравнения (1 ), и не требует специальных пояснений.
7. Определение закона распределения а(Д). В заключение« данного параграфа рассмотрим вопрос об определении закона распределения случайной функции a . ( t ) . Будем исходить для простоты из уравнения (1 ), считая кинетический момент И не случайным.
Если имеется жидкостное трение, момент которого определя ется формулой (1 1 ), то уравнение (1 ) интегрируется и мы получим
(при а(0 ) = 0 ) |
|
«(0 = - ^ 0 (О + | 0(0). |
(4.67) |
Следовательно, если законы распределения ординат случай ной функции 0 (<) известны до второго порядка включительно, то задача сводится к определению закона распределения суммы зависимых слагаемых по формуле (1.42). Если t достаточно ве лико, чтобы Ѳ(0) и 0 (t) можно было считать независимыми вели чинами, то мы будем иметь обычную задачу о композиции законов распределения, для решения которой в данном случае достаточно знать только одномерный закон распределения ординаты Ѳ(і). Таким образом, в случае жидкостного трения задача определения закона распределения ординат случайной функции a ( t ) решается без каких-либо принципиальных затруднений для любого закона распределения ординат случайной функции Ѳ(£).
Эта задача еще более упрощается, если функция Ѳ(£) нормаль
ная. |
В этом случае |
вследствие линейной зависимости между |
|||
а (t) |
и Ѳ( t ) |
закон распределения ординат а ( t ) |
будет также нормаль |
||
ным, и следовательно, найденные выше â(t) |
и D[a(i)l |
однозначно |
|||
определяют этот закон распределения. |
является |
жидкостным |
|||
Задача |
усложняется, если трение не |
||||
и момент |
сил трения |
определяется, например, формулой (19). |
В этом случае соотношение между а (() и Ѳ(t) уже перестает быть
линейным и, следовательно, закон |
распределения а ( t ) уже не |
будет нормальным, даже когда Ѳ(t) |
нормально. |
Для сухого трения можно воспользоваться приближенным разложением плотности вероятности /(а ) в ряд Шарлье (1. 44), для чего требуется определение моментов а (t) более высоких порядков, чем второй. Оно может быть выполнено тем же методом, каким были определены первые два момента, однако с повышением порядка момента вычислительная сложность расчетов увеличи вается.
В качестве второго пути можно указать применение марков ских процессов. Этот способ пригоден, когда случайная функция
5 4.11 |
У Р А В Н Е Н И Е |
1-ГО П О Р Я Д К А Г.ЕЗ ЗАВИСИМ ОЙ П Е Р Е М Е Н Н О Й |
181 |
d(t) |
нормальна и |
имеет дробно-рациональную спектральную |
плотность.
Для иллюстрации этого метода предположим сначала, что
А'іі (т) -- ф |
|
(4.68) |
|
и, следовательно |ем. (1 . 1 2 2 )|, |
|
|
|
S6 (ü>) ~ |
" ( » 2 + Р2)' |
|
(4.69) |
|
|
||
Как следует из формул (1. |
137) и (1. |
138), |
Ѳ(it) в этом случае |
является марковским процессом и может рассматриваться как решение дифференциального уравнения
(*)+ |
рб = <*(*), |
(4.70) |
где I (t) — белый шум, a c2 = 2 |
a?fA. |
|
Рассматривая уравнение (1), с учетом (19) и (70), как совмест ную систему уравнений, замечаем, что функции U1{t)= a(t) и tfs(f)=Ö( 0 являются компонентами двумерного марковского процесса. Определяя для этого случая коэффициенты уравнений Колмогорова по формулам, аналогичным (1. 132) и (1. 133), получим
Qx ■ |
а2= — рг/2, |
|
ai = jf |
sign у2, |
|
^ 1 1 -- ^ 1 2 |
-- |
^22 -- С^. |
Следовательно, второе уравнение Колмогорова для плот ности вероятности / системы случайных величин Y 1= U l ( т), Y2—U2( т) будет иметь вид
df I |
д |
^psign (if,)/]— |
д , |
г _ |
£2 |
(4.71) |
|
|
|
2 ду\ |
|||||
д і ' д у х |
|
|
|
|
|
||
Решение |
последнего |
уравнения |
при |
начальных условиях |
|||
/|т_0 = 8 (у1 —хг) §(г/2 —х2) |
даст двумерную |
плотность вероятности |
/(Уі,У2)- Искомая плотность вероятности/ ( а) может быть полу чена путем интегрирования / (уъ у2) по у2:
/(« )= = /(г/з) = J І(Уі, У2)(ІУг- |
(4.72) |
—СО |
|
Таким образом, если Ѳ(і!) — нормальная случайная функция, имеющая спектральную плотность (69), то задача сводится к инте грированию уравнения (71). Хотя при решении этого уравнения применим метод расщепления Фурье, решение получается доста точно сложным. Еще более сложным будет решение, если в знаме нателе £б(й>) стоит полином ш более высокой степени, чем второй.
182 |
ГУ, О П И С Ы В А Е М Ы Е Л И Н Е И Н Ы М И У Р А В Н Е Н И Я М И |
[ГЛ. 4 |
Формула (6 8 ) обычно слишком плохо аппроксимирует кор реляционную функцию угловой скорости основания. Более хоро ший результат можно получить, если положить
А ’е (х) == о|е_11'т| ^ros Хх — у- sin X |х |). |
(4.73) |
Последняя формула, например, часто применяется для харак теристики угловой скорости качки корабля, колебаний самолета и в ряде других случаев. Корреляционной функции (73) соответ ствует спектральная плотность [см. (1.111) и (1.124)]
2алр.ш-
|
S b Н = |
к |
+ fj.2 + Х2)2 |
_ 4 Х 2Ы2] • |
( ' i - 7 ^ ) |
|
Следовательно, |
если |
Ѳ(t) — нормальная случайная |
функция, |
|||
то Ѳ(t) |
будет компонентой |
двумерного марковского |
процесса, |
|||
а а(t) |
является |
уже |
компонентой |
трехмерного марковского |
процесса. Для этого процесса уравнение Колмогорова также может быть легко составлено [см. (1.147)], однако его решение еще более усложняется, поскольку появляется новая независимая переменная.
В большинстве практических задач необходимости в точном определении / ( а) обычно не возникает, так как ошибки ГУ вхо дят составной частью в ошибки более сложных систем, закон рас пределения на выходе которых уже с высокой точностью можно считать нормальным, так как они складываются из большого числа независимых слагаемых и имеются условия практической при менимости центральной предельной теоремы. В этом случае достаточно знать первые два момента ошибок ГУ, определение которых было рассмотрено выше.
§ 4.2. Линейное уравнение первого порядка, содержащее зависимую переменную
1. Вывод общих формул. Уравнение первого порядка, со держащее зависимую переменную, описывает поведение различ ных ГУ. Примером может служить уравнение движения гиро вертикали с маятниковой коррекцией, имеющей линейную ха рактеристику [см. (3.58)]
“ + х2а ~ хгХі W + у/ ^ 2 > |
(4-75) |
где XiW —угол отклонения от вертикали физического маятника, х2 — удельная скорость коррекции, а М 2 — суммарный возмуща ющий момент, возникающий вследствие наличия момента трения, момента из-за неуравновешенности гироскопа и т. п.
§ 4.2] У Р А В Н Е Н И Е 1-ГО П О Р Я Д К А С ЗАВИСИМ ОЙ П Е Р Е М Е Н Н О Й 183
Аналогичный вид имеет и упрощенное уравнение гиротахо
метра с тремя степенями свободы [см. (3.141)]: |
|
||||
fS+ l ß |
= co + |
i м, |
(4.76) |
||
гиромагнитного компаса [см. |
(3. 321)| |
|
|
||
d -|- ха = |
— и{ + |
X (е -}- 8), |
(4.77) |
||
поплавкового интегрирующего |
гироскопа [см. (3.173) |
и (3.166) |
|||
при со^ = ф] |
|
|
t |
|
|
ß - f i ß = |
± c p + f |
(4.78) |
|||
о |
|||||
|
|
|
|
и ряда других гироскопических устройств.
Все приведенные выше уравнения принадлежат к одному
математическому типу и могут быть записаны в виде |
|
|
о . а ха. = X (t) |
(aj)>0), |
(4-79) |
где постоянная аг и случайная функция X(t) в различных слу чаях могут иметь различные значения и различный физический
смысл. |
Так, |
например, |
для |
гировертикали |
в |
соответствии |
||
с (75) |
X (t) |
является |
суммой |
двух |
(обычно |
стационарных) |
||
случайных функций x ^ |
(t) |
и -ң М2 (t), а |
в случае |
интегрирую |
щего гироскопа, согласно (78), содержит интеграл от стационар ной функции и, следовательно, не является стационарной. Для гиротахометра, как ясно из уравнения (76), случайная функция X (t) является суммой полезного сигнала <a(f) и помехи; в других случаях полезная составляющая отсутствует.
Рассмотрим получение моментов ординат случайной функции а (t). Для этой цели запишем решение уравнения (79) в явном виде
t
a(t) — a (0) е~а'1 ^ I (t — т) X (т) dz, |
(4.80) |
о |
|
где импульсная переходная функция l ( z ) в данном случае имеет вид экспоненты
I (т) = |
е_а‘т, |
(4.81) |
а начальное значение а (0 ) угла |
а в дальнейшем мы не |
будем |
учитывать, рассматривая в качестве ошибки ГУ отклонение,
возникающее |
за |
время |
t. |
а (t) является результатом |
Формула |
(80) |
показывает, что |
||
применения |
линейного |
оператора |
к случайной функции X (t). |