Рис. 22.14
и автономной систем со спект ральными плотностями 5р(ю) и Sa(tO).
Предположим, что спектраль ные плотности Sp(co) и Sa(to) из вестны (рис. 22.15), причем энер гия помехи «р(0 распределена в широкой полосе частот, а энер
гия помехи на(г) сосредоточена в низкочастотной области. Для измерителей, построенных на различных физических принципах, это предположение на практике обычно выпол няется. Как будет показано, оно является необходимым условием высокой эффективности комплексирования.
Пользуясь структурной схемой, приведенной на рис. 22.14, получим сигнал на выходе фильтра с передаточной функцией Н(р):
na(t) = na(t)H(p) -np(t)H(p), |
(22.59) |
где p = d/dz — оператор, дифференцирования.
Очевидно, что для повышения точности комплексной системы необходимо передаточную функцию фильтра Н(р) выбрать таким образом, чтобы фильтр наилучшим образом воспроизводил на своем выходе низкочастотную помеху автономной системы ла(г), рассматриваемую в качестве полезного сигнала. При этом шумом измерения является широкополосная помеха радиоканала np(f). Для полной априорной определенности структуру рассматриваемого фильтра можно получить используя методы теории калмановской фильтрации.
Таким образом, сигнал на выходе комплексной систе мы можно записать следующим образом:
z(r) = X(O+«a(O[l-H(p)]+np(t)H(p). |
(22.60) |
Анализ этого выражения показывает, что если бы |
фильтр с передаточной функцией |
Н(р) абсолютно точно |
воспроизводил погрешность автономного измерителя иа(г)
и полностью подавлял помехи радиоканала np(t), то выходной сигнал комплексной системы идеально точно воспроизводил бы полезный навигационный параметр Х(?). Фактически же из-за того, что спектры помех перекры ваются и идеальную фильтрацию осуществить невозможно, в выходном сигнале помимо полезного навигационного параметра Х(?) содержится ошибка е(г), которая, согласно
(22.60), |
имеет вид |
|
40 = z(0 - Х(0 = Ha(0 [1 - Н(р) ] + np(l)H(p). |
(22.61) |
Дисперсия этой ошибки с учетом обычно выполняюще |
гося условия независимости помех иа(0. и ир(0 |
|
= |
J [|l-H(7Co)|2Sa(co)+|/7(7-co)|25p(co)]dco. |
(22.62) |
2л - оо
Из соотношения (22.62) следует, что если спектральные плотности 5а(со) и Sp(co) сосредоточены в различных частотных областях, то дисперсия а2 существенно меньше дисперсий погрешностей каждого измерителя.
Отметим также, что при использовании способа ком пенсации в комплексной системе не возникает динамических ошибок, так как ошибка системы £(/), как это следует из выражения (22.61), не зависит от Х(/). Поэтому на выбор параметров системы не влияет модель процесса Х(?), что является большим достоинством этого способа при отсутствии априорной информации относительно Х(г). Конечно, отказ от априорных данных ухудшает точностные характеристики комплексной системы.
Системы, в которых ошибки не зависят от характе ристик полезного сообщения Х(?), называются инвариант ными по отношению к X(z). Таким образом, комплексная навигационная система, построенная по способу компенса ции, является инвариантной относительно измеряемого навигационного параметра Х(г). Необходимое' условие осу ществления инвариантной системы — наличие в ней по крайней мере двух измерителей.
Способ фильтрации. Другим возможным способом комплексирования навигационных измерителей является объединение их по схеме фильтрации (рис. 22.16). В этом
случае выходной сигнал |
комплексной |
системы |
N |
N |
N |
|
40 = Z УМН((Р)=1(0 X Н,(р) + £ |
(22.63) |
1=1 |
> = 1 |
1=1 |
|
Для того чтобы этот сигнал не содержал динамических |
ошибок, необходимо |
выполнить условие |
N |
|
|
Хад=1, |
(22.64) |
i = 1 |
условие инвариантно |
т. e. |
сти |
комплексной |
системы |
относительно навигационно го параметра Х(?). При его выполнении выходной сиг нал системы
=+ £Х(')Я,.(р).(22.65)
i= 1
Следовательно, результирующая ошибка |
измерения |
имеет только флуктуационную составляющую |
|
N |
|
E(O = z(r)-X(r)= £ |
(22.66) |
При Я=2 и выполнении условия инвариантности (22.64) выражения (22.66) и (22.61) для ошибок измерения на выходе комплексных систем, организованных по спо собам фильтрации и компенсации, совпадают.
Способ комплексирования с введением дополнительной информации внутрь контура слежения радиотехнического измерителя. Введение информации о параметрах движения цели в контур слежения (рис. 22.17, а) позволяет улучшить характеристики точности и помехоустойчивости радиотех нических измерителей в режиме слежения; уменьшить вероятности ложных захватов и срыва слежения; обеспечить возможность сопровождения при пропадании радиосигна лов; сократить время поиска.
Найдем выходной сигнал z(t) комплексной системы, приведенной на рис. 22.17, а. Так как рассматриваемая система линейна, то, воспользовавшись принципом супер позиции, получим
z(0 = Я, (р)ур(?) + Я2(р)уа(0, |
(22.67) |
где |
|
|
Н^(Р)Н^(Р) . |
К(р)^(р) |
(22,68) |
Я1(р) |
1{р) 1+я(1’(р)я<2’(р)' |
Нетрудно видеть, |
что выражение (22.67) |
совпадает |
с (22.63) и, следовательно, данная схема может быть представлена в виде двухканальной схемы фильтрации, приведенной на рис. 22.17, б. При этом условие инвари антности (22.64) выполняется в том случае, когда
К(р)Н^(р)=\. |
|
(22.69) |
С учетом условия (22.69) и формул .(22.58) выражение |
(22.67) принимает вид |
|
|
Н<1Чр)Н<2Чр) |
1 |
(22.70) |
^0-Х(04-,+Д)Д-Л(0 |
+ , + д,,,м„,а(;,)-..(0. |
где второе слагаемое в правой части представляет собой ошибки оценивания навигационного параметра X(z) за счет помех в радиотехническом измерителе, а третье—за счет помех в автономном. Таким образом, при выполнении условия инвариантности (22.69) эффективную ширину по лосы пропускания радиотехнического измерителя следует выбирать только исходя из условия минимума суммы дисперсий этих ошибок.
Синтез комплексных систем. Из изложенного следует,
что к задаче синтеза комплексных систем возможны два основных подхода. Сформулируем их, ограничиваясь рам ками теории линейной фильтрации.
Первый подход предполагает заданными модели из меряемого навигационного параметра k(r) и погрешностей автономных датчиков иа(Г) в виде уравнений состояний и измерений. Это позволяет использовать для синтеза оптимальной комплексной системы результаты теории калмановской фильтрации. Соответствующий многомерный фильтр полностью определяется уравнениями (4.88) — (4.90). Достоинством этого подхода является возможность получения прн указанных условиях оптимальной структуры и параметров комплексной системы. В то же время необходимо отметить, что предположение о точном знании модели навигационного параметра k(z) редко выполняется на практике. К тому же размерность оптимального фильтра совместной обработки данных оказывается до вольно высокой, что приводит к значительному росту объема вычислений при^реализации фильтра на ЭВМ.
Поэтому часто для синтеза комплексных систем при меняют второй подход, при котором данные о модели к(?) исключаются из рассмотрения и комплексирующий фильтр оказывается инвариантным относительно парамет ров движения объекта. Использование этого подхода позволяет упростить алгоритм фильтрации и снизить требования к необходимому объему априорных данных.
Синтез комплексных систем на основе принципа инва риантности можно произвести различными способами. Для этой цели можно применить рассмотренный ранее метод компенсации (см. рис. 22.14) и синтезировать фильтр Калмана, осуществляющий оптимальное оценивание погреш ностей автономного датчика иа(?), наблюдаемых на фоне помех радиотехнического измерителя. Однако при програм мной реализации алгоритмов комплексирования на ЭВМ принцип инвариантности удобнее осуществить путем сос тавления специальной системы разностных уравнений, в которой исключается ряд неизвестных, описывающих движение объекта. При этом сигналы тех измерительных систем, математические модели которых неизвестны или малодостоверны, рассматриваются как компоненты вектора управления. В качестве примера синтезируем инвариантную комплексную систему, содержащую импульсный радиодаль номер и датчик воздушной скорости, установленные на борту движущегося объекта.
Прежде всего необходимо задать модель погрешностей датчика воздушной скорости. Анализ работы таких датчи ков показывает, что модель погрешностей измерителя скорости может быть описана экспоненциально коррели рованным процессом ла(Г), определяемым стохастическим дифференциальным уравнением вида
(22.71)
<и Га Га
где Та—постоянная времени, зависящая от состояния атмосферы и скорости полета; ^а(?) — белый гауссовский шум со спектральной плотностью 2afTa; ста—дисперсия погрешности измерения скорости в этом датчике.
Составим уравнение состояния системы. Так как модель движения объекта предполагается неизвестной, то в качестве компонентов вектора состояния целесообразно выбрать текущую дальность объекта и погрешность авто номного датчика скорости, т. е. задать вектор состояния Х.(г) в виде Х(г) = ||D(r)wa(r)||T. Йстинное значение скорости объекта Г(?) удовлетворяет следующему соотношению:
H0=-^-=J>(0-”a(O- |
(22.72) |
|
Объединив уравнения (22.71) и (22.72) в систему при |
матричной форме записи, получим |
|
(22.73)
В этом уравнении состояния выходной сигнал датчика автономной скорости уа(г) выступает как известное управ ляющее воздействие (вектор управления). Поэтому уравне ние измерений описывает только процесс формирования
отсчетов |
дальности: |
|
|
уР(0 = 11 1 |
oil |
о(г) |
+«р(0- |
(22.74) |
|
|
«а (О |
|
Структурная схема рассматриваемой модели процесса и измерений приведена на рнс. 22.18.
Приступим непосредственно к синтезу дискретного алгоритма комплексирования, реализуемого на ЭВМ. Для этого, используя соотношения (22.45)— (22.47), осуществим переход от дифференциального уравнения (22.73) к уравне нию в конечных разностях. Вычислим переходную матрицу дискретной системы. По определению [13],
Рис. 22.18
В; = В(г;, fi_1) = exp{F(r1 —ri_1)} = £“1{||pl—F||-1}, |
(22.75) |
где £-1{-} — оператор обратного преобразования Лапласа.
Подставив |
в (22.75) |
выражение для |
матрицы |
F из |
(22.73), |
найдем |
|
|
|
|
|
|
Bf = |
О1 |
т |
О |
|
|
|
|
|
(22.76) |
где |
w = exp{—Тд/Га}; |
Тл—период дискретизации. |
|
|
Из |
выражения |
(22.45) |
с |
учетом соотношений |
(22.73) |
и (22.76) следует, |
что |
|
|
|
|
к, — В,к, _ J + |
|
|
1 |
d* |
Уа(.-1) + |
|
|
|
|
|
|
|
О |
|
|
|
|
|
|
|
О1 |
|
^a(i- 1)* |
|
|
(22.77) |
|
Вычислив входящие в (22.77) интегралы, окончательно |
получим |
|
|
|
|
|
|
|
Dt |
|
О |
т |
|
Л-1 |
|
|
|
”ai |
|
О |
|
^a(i- 1) |
|
|
2>a((-1) |
|
|
|
|
|
|
(22.78) |
|
|
|
|
|
|
|
|
|
где |
|
|
I, — эквивалентные дискретные белые |
шумы. |
|
Уравнение измерений отличается от (22.74) только |
дискретным |
временем |
|
|
|
|
|
ypi=ll 1 |
он |
Di |
|
|
|
|
|
(22.79) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Выражения (22.78) и (22.79) позволяют построить на |
основе |
(4.88) — (4.90) |
соответствующий |
фильтр Калмана, |
особенностью которого является использование в нем для вычисления оценок экстраполяции вектор управления II ^аО-^О Г.
Таким образом, искомый алгоритм фильтрации для рассматриваемой инвариантной комплексной системы мо жет быть записан следующим образом:
Di |
+ |
ТцУa(i - 1) + |
|
|
Wna(i-l) |
о |
Рис. 22.19
Структурная схема фильтра, реализующего этот алго ритм, приведена на рис. 22.19. На вход фильтра поступают отсчеты дальности ypj и автономной скорости уа;. Погреш ности измерения последней корректируются путем вычита ния из них оценок погрешностей йа1 с соответствующими масштабными коэффициентами. Эта схема наряду с опти мальной фильтрацией обеспечивает надежное формирова ние оценок при кратковременных пропаданиях радиосигна лов, поскольку при нулевых значениях коэффициентов Glf и G2i происходит интегрирование автономных данных о скорости yai с помощью дискретного интегратора, вырабатывающего оценку дальности Dt. На практике алгоритм (22.81) реализуется программным путем.
§ 22.4. НАВИГАЦИОННЫЕ КОМПЛЕКСЫ САМОЛЕТОВ И СУДОВ
Под навигационным комплексом понимают совокупность бортовых измерительных средств и вычислителей, позво ляющих определять местоположение и скорость самолета (судна) относительно Земли. Ни один из существующих навигационных измерителей не может полностью решить эти задачи, так как каждый из них в отдельности не обеспечивает необходимой точности, помехозащищенности или надежности.
Задачи, решаемые навигационным комплексом, много образны. Среди них одной из важнейших является счисле ние пути, обеспечивающее непрерывное измерение коорди наты объекта. Основным недостатком систем счисления является ухудшение точности определения координат с уве личением времени работы. Поэтому для получения требу-
Рис. 22.20
емой точности счислимые координаты необходимо непре рывно или периодически корректировать на основании информации, поступающей от радиотехнических измерите лей, т. е. осуществлять комплексную обработку данных.
Структурная схема типового навигационного комплек са самолета приведена на рис. 22.20 [13]. Основу этого комплекса составляет инерциальная навигационная система (ИНС) на гиростабилизированной платформе. Она измеря ет как угловое положение самолета (углы крена, тангажа, рыскания и их производные), так и составляющие ускоре ния и скорости. Скорость самолета измеряется также с помощью ДИСС и датчика воздушной скорости, вхо дящего в состав системы воздушных сигналов (СВС). В качестве вспомогательного измерителя курса использу ется система курсовертикали (СКВ). Высота и скорость ее изменения измеряются с помощью радиовысотомеров (РВ). Сигналы этих устройств обрабатываются в вычис лительном устройстве, являющемся частью распределенной бортовой вычислительной системы. В качестве систем коррекции координат местоположения самолета использу ются данные радиотехнических систем ближней (РСБН) и дальней (РСДН) навигации (таких, как «Омега», «Лоран- С» или системы с использованием ИСЗ), бортовых' РЛС,
корреляционно-экстремальных систем, а также данные, получаемые с выхода других измерителей, например аст рономических ориентиров, оптических или электронно оптических визиров.
В навигационных комплексах с более высокой степенью интеграции оборудования используются обратные связи (показаны на рис. 22.20 пунктирными линиями). За счет этих связей обеспечиваются коррекция положения гироплат формы ИНС, предварительная настройка ДИСС по данным датчика воздушной скорости или ИНС, установка визиров в предполагаемое местоположение ориентиров и т. п. Так как системы, входящие в навигационный комплекс, опре деляют навигационные параметры в собственной системе координат, в алгоритмах навигационного вычислительного устройства предусмотрена процедура пересчета данных этих систем в основную систему координат, в которой осуществляется счисление пути.
Навигационный комплекс является составной частью пилотажно-навигационного комплекса (ПНК), который вк лючает в себя также систему автоматического управления самолетом и систему индикации и отображения пилотаж но-навигационной информации. ПНК предназначен для навигации и пилотирования самолета на всех этапах полета. В круг задач, решаемых ПНК, помимо непрерыв ного определения координат местоположения самолета, счисления пути и его коррекции входят программирование маршрута полета, вычисление и передача в САУ управ ляющих сигналов, выдача информации системам отобра жения и индикации, автоматический контроль исправности бортовых устройств и систем ПНК, а также автоматическая стабилизация и управление самолетом во всех режимах полета.
Навигационные комплексы морских судов имеют схо жую структуру. На рис. 22.21 приведена структурная схема интегрированного навигационного комплекса «Data Bridge» норвежской фирмы «Norcontrol» *, предназначенного для автоматизации судовождения и предотвращения столкно вений. Счисление пути в этом комплексе осуществляется по данным лага и гирокомпаса. В качестве систем коррекции координат местоположения используются на
вигационные |
системы |
Декка |
(непрерывная |
коррекция |
в условиях прибрежного плавания), «Омега» , |
«Лоран-С», |
* См.: |
Москвич Г. |
За |
электронными |
картами — |
будущсе.//Морской флот. |
1986. №11. С. 21—24. |
|