Файл: Свешников А.А. Вероятностные методы в прикладной теории гироскопов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 212

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

5 9 . 6 ] С У Щ Н О С Т Ь М Е Т О Д А О П Т И М А Л Ь Н О Й Ф И Л Ь Т Р А Ц И И К Л Л М Л Н А 495

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

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

Другой способ состоит в использовании заранее вычисленных коэффициентов усиления.

Достоинствами рассматриваемого метода оптимальной фильт­ рации являются:

1)Он дает, как и другие методы оптимизации, наилучшие в смысле минимума погрешности возможные оценки, исходя из из­ вестных вероятностных характеристик как входных переменных, так и ошибок измерений.

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

3)Этот метод позволяет получить при решении оптимальной задачи непосредственно структуру оптимального фильтра, кото­ рая практически может быть осуществлена.

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

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

для конечного и бесконечного времени наблюдения.

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

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


496

В Ы Б О Р О П Т И М А Л Ь Н Ы Х С Х Е М И П А Р А М Е Т Р О В ГУ

[ Г Л . 9

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

Здесь они применяются при решении различных задач, в част­ ности, при первоначальной выставке гиростабилизированных площадок, при осуществлении демпфирования ИНС по скорости и коррекции координат места, для улучшения динамических свойств различных устройств ИНС при отсутствии внешней инфор­ мации, для определения источников ошибок различных инерциаль­ ных приборов и осуществления их коррекции, при решении задач гирокомпасирования (определение курса объекта без использо­ вания гирокомпаса), при разделении действительных ускорений объекта от «ускорений», вызванных инструментальными ошиб­ ками инерциальных приборов с целью уменьшения их влияния на точность навигации и т. д.

§9.7. Примеры определения оптимальных схем

ипараметров ГУ

Рассмотрим несколько примеров определения оптимальных систем, поясняющих решение задач трех типов, указанных в § 9.1:

а) определение оптимальной динамической системы без ка­ ких-либо ограничений, наложенных на ее структуру;

б) нахождение параметров ГУ при заданной его структурной схеме;

в) определение корректирующего устройства при заданной неизменяемой части гироскопической системы (гироскоп и другие элементы).

Пример 9.1. Определить оптимальную передаточную функцию измерителя угловой скорости бортовой качки корабля (гиротахо­ метра), если память системы не ограничена, угол крена 6(t) — стационарная случайная функция, а на вход системы поступает сумма полезного сигнала £/(£)= Ѳ(2) и помехи V(t), возникающей вследствие различных возмущающих факторов. Дано, что й—г>=0,

а спектральная плотность

S u ( со)

определяется

формулой

(2.15)

SuИ =

2 о2р

Ъ1

 

 

( £ 2 =

Р2 + X2, а*=

\х\— \ 3).

(9.77)

0)4 -f- 2 а-со2

 

 

В целях упрощения задачи примем сперва, что V(t) можно

считать белым шумом, т.

е.

положим

 

 

 

 

 

5

»

=

- é r c ’

 

<9-78)

где С — интенсивность

белого

шума.

S v ( ш) является

Сделанное

предположение

о

постоянстве

приближением,

которое

тем

 

лучше соответствует действитель-



§ 9.7]

ПРИМ ЕРЫ ОП РЕДЕЛЕН И Я СХЕМ И ПАРАМЕТРОВ ГУ

497

ности,

чем быстрее затухает К ѵ(т) с ростом t сравнительно с за­

туханием і£а(т). Это предположение в рассматриваемой задаче

обычно может быть принято.

 

и

V(t) взаимно

Предположим, что случайные функции U(t)

не коррелированьи

 

 

 

сигнала X(t)

Р е ш е н и е .

По условию задачи для входного

имеем

X(t) = U(t)+V(t).

 

[(9.79)

 

 

Требуемый сигнал Z(t) в данном случае есть

производная от

U(t), т. е.

 

 

 

 

 

 

Z{t) = ± U { t ) .

 

(9.80)

Сигнал Y(t)~на выходе

системы связан с X(t) соотношением

 

Y<t) =

LX(t),

 

1(9.81)

где L — искомый оптимальный оператор.^

 

рассматривае­

Для определения передаточной функции Ь(ш)

мой динамической системы, соответствующей оператору L, имеем

окончательное

выражение

(19),

справедливое

при отсутствии

запаздывания.

Согласно (79) спектральная плотность5 Sx ( ш) вследствие не­ связанности случайных функций U(t) и V(t) определяется со­ отношением

Sx (^) — Su(o>) +

St (o)).

 

(9.82)

При этом взаимная спектральная

плотность

Sxz (со) имеет вид

s x, H = iiüSu(«>)•

 

(9.83)

Обозначив

 

 

(9.84)

G = 4з>,

 

 

перепишем (77) следующим образом:

 

 

 

С ( \ ______________ ______________

(9.85)

0 9 — :2к )2 + 2 ш + 62)

_ 2Хо) + 62)

 

 

*

Учитывая (85) и (78), имеем

(9.86)

что после приведения к общему знаменателю дает

„ , , _ 1

С \<оі +

2 (62 2 1 2 ) М2 4 - щ +

G62

(9.87)

2%

(со* +

2Хш + 62) (са2 — 2U +

б«) *

 

Обозначая

 

 

 

(9.88)

 

 

ix.2 — X2

 


498 ВЫ БОР ОПТИМАЛЬНЫХ СХЕМ И ПАРАМЕТРОВ ГУ [ГЛ. 9

и учитывая

(84), преобразуем (87) к виду

 

 

 

 

1

Сші +

2 6 , 0 ) 2 +

62 (G + С62)

(9.89)

 

S » = - 271

)2 + 2Хсо + 62) (и)2 _

2Хо) + 62)

Найдем нули функции Sx (ш); имеем

 

 

 

 

 

С

+

2 Ѵ 2 + Ь 2 2 +

| -

) 1

= 0.

(9.90)

Обозначим

 

 

w — со2;

 

 

 

(9.91)

 

 

 

 

 

 

тогда уравнение (90) перепишем в виде

 

 

 

 

 

w2

2 bxw -f- Ъ2 2 -}-

 

=

О,

(9.92)

откуда

 

 

 

 

 

 

 

 

 

 

w'і.* = - Ь і ± \

 

% - & ( & + %)•

(9.93)

Принимая во внимание (84)

и (8 8 ),

получим

 

}/ГЪ\ -

Ь2 2 +

| - ) =

/([X 2 -

^ 2) 2 -

-2

+

X2) ([X2 +

X2 + £ ) ,

откуда следует, что при р <( X подкоренное выражение будет отри­

цательным.

 

 

 

 

 

 

 

 

 

Обозначая

 

 

 

 

 

 

 

 

 

 

d = y

- Ъ \ + Ѵ-(Ъ*+^),

 

(9.94)

имеем

 

 

w12 = b1+ id.

 

 

 

 

 

 

 

 

 

Принимая во внимание (91), находим

 

 

 

 

 

 

U),joj3>4==: +

V—b1 +

Id

(9.95)

Пользуясь формулой Муавра, получим

 

 

 

 

 

со, =

(6 2 +

d2)'/<(cos

----і sin у ) ,

 

 

ü>2 =

— (b\ -f d2)'l<('cos у — i sin у ) ,

(9.96)

 

(ö3 =

(b2 +

d2),/i (cos

-f- г sin у ) ,

 

 

o)4 — — + d2)4>(cos у -f і sin y ) ,

где

(9.97)

cp = arctg