Файл: Свешников А.А. Вероятностные методы в прикладной теории гироскопов.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).