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