Исследование и оптимизация интегрированной системы измерения параметров полета летательного аппарата вблизи поверхности - davaiknam.ru o_O
Главная
Поиск по ключевым словам:
страница 1
Похожие работы
Название работы Кол-во страниц Размер
Параметрическая оптимизация информационно-измерительной системы определения... 1 279.88kb.
Нейросетевой супервизор автопилота летательного аппарата в докладе... 1 77.85kb.
Рейтинг за работу Преподаватель к т. н., доцент / И. С. Зегжда /... 1 17.83kb.
Пилот сверхлегкого летательного аппарата: кто это 1 104.7kb.
Исследование гидродинамических параметров и молекулярного состава... 1 195.78kb.
Искусственные спутники Земли 1 15.3kb.
Направления создания интегрированной информационной системы управления... 1 35.38kb.
Пожалуй нет более важных понятий, характеризующих летные свойства... 1 105.01kb.
Закон «О техническом регулировании» 1 18.92kb.
Методика обработки наблюдений полета болида для определения его траектории 1 46.86kb.
Технический паспорт космического летательного аппарата 1 11.16kb.
Метеорологические элементы и явления погоды, определяющие условия... 1 127.81kb.
Направления изучения представлений о справедливости 1 202.17kb.

Исследование и оптимизация интегрированной системы измерения параметров полета летательного - страница №1/1



На правах рукописи

Румянцева Елизавета Анатольевна



ИССЛЕДОВАНИЕ И ОПТИМИЗАЦИЯ ИНТЕГРИРОВАННОЙ

СИСТЕМЫ ИЗМЕРЕНИЯ ПАРАМЕТРОВ ПОЛЕТА

ЛЕТАТЕЛЬНОГО АППАРАТА ВБЛИЗИ ПОВЕРХНОСТИ


Специальность 05.13.01

«Системный анализ, управление и обработка информации»

(в технике и технологиях)


Автореферат

Диссертации на соискание ученой степени

кандидата технических наук

Санкт-Петербург

2008

Работа выполнена в государственном образовательном учреждении высшего профессионального образования «Санкт-Петербургский государственный университет аэрокосмического приборостроения»



Научный руководитель:

Доктор технических наук,

профессор, Заслуженный

деятель науки РФ Небылов Александр Владимирович




Официальные оппоненты:
Доктор технических наук,

Профессор Коновалов Александр Сергеевич

Кандидат технических наук Розенгауз Михаил Борисович


Ведущая организация: Федеральное государственное унитарное предприятие Санкт-Петербургское опытно-конструкторское бюро «Электроавтоматика»

Защита диссертации состоится “___ ” ______________ 2008 г. в ____ часов на заседании диссертационного совета Д 212.233.02 при Государственном образовательном учреждении высшего профессионального образования «Санкт-Петербургский государственный университет аэрокосмического приборостроения» по адресу: 190000, Санкт-Петербург, ул. Большая Морская, д. 67, ГУАП.

С диссертацией можно ознакомиться в библиотеке ГУАП.
Автореферат разослан “___ ” ______________ 2008 г.

Ученый секретарь диссертационного совета,

Доктор технических наук, профессор Л.А. Осипов
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. Задачи сбора, обработки и анализа информации о параметрах полета летательного аппарата (ЛА) являются одними из важнейших для обеспечения безаварийного полета, повышения точности и надежности измерительной системы. Многообразие современных летательных аппаратов определяется конструктивными особенностями и областями их применения. Повышенные требования предъявляются к полету на низкой высоте. Актуальность задачи повышения точности навигации для обеспечения безаварийного и экономичного движения аппарата в диапазоне малых высот связана с совершенствованием неводоизмещающих морских аппаратов, к которым относятся суда на управляемых подводных крыльях, воздушной подушке и экранопланы, а также применением системы управления движением на гидросамолетах и экранопланах, осуществляющих полет на предельно малых высотах при решении определенных функциональных задач. Диапазон рассматриваемых малых высот может характеризоваться единицами метров. Конкретная высота движения выбирается и стабилизируется как минимально допустимая из условий достаточно малой вероятности столкновения с препятствиями с учетом характера неровности подстилающей поверхности, например, балльности и характеристик трехмерности морского волнения. Этапы взлета и посадки также характеризуются повышенными требованиями к точности измерения малой высоты и параметров неровностей.

Большинство публикаций по проблеме измерения параметров морских волн посвящено изучению вопросов, связанных с исследованием методов и средств определения интегральных характеристик волнения. Такие исследования наиболее эффективны при расположении аппаратуры на авиационных или космических носителях либо при визировании морской поверхности судовыми радиолокаторами под малыми углами. Все они предполагают фактически не отслеживание профиля волны, а оценку интегральных характеристик волнения по косвенным признакам. Кроме того, в ряде работ рассмотрены волнографы контактного типа, а также принципы построения океанографической аппаратуры для специальных научно-исследовательских судов, где датчики приспосабливают для работы в условиях дрейфующего судна. Очевидно, что все эти исследования не могут служить основой для создания аппаратуры измерения характеристик морских волн, пригодной для полета на предельно низких высотах. В настоящее время наиболее эффективным представляется вариант совместной обработки показаний неконтактных датчиков геометрической высоты и инерциальных датчиков вертикального ускорения. Важнейшими из достоинств этого метода являются его сравнительная простота и доступность средств измерения.


Цель диссертационной работы состоит в разработке и совершенствовании теоретических основ построения интегрированной системы датчиков высоты; исследовании точности измерения малых высот при полете над земной и морской поверхностью при совместной обработке показаний локационных и инерциальных датчиков.
Основные задачи. Для достижения поставленной цели в диссертационной работе решались следующие основные задачи:

  1. Исследование и общая классификация существующих методов измерения параметров полета ЛА.

  2. Выбор и обоснование состава датчиков параметров низковысотного полета с учетом современных тенденций их совершенствования.

  3. Исследование основных вероятностных свойств и спекрально-корреляционных характеристик морского волнения.

  4. Разработка алгоритмов расчета спектра погрешности локационного датчика (радиовысотомера) при полете на малой высоте над взволнованным морем.

  5. Оптимизация интегрированной системы измерения параметров низковысотного полета по критериям точности и отказоустойчивости.

  6. Сравнительный анализ точности измерения высоты полета при построении оптимальной и робастной систем.

  7. Экспериментальные исследования точности измерения высоты полета ЛА над морской поверхностью.


Методы исследования. При выполнении диссертационных исследований использовались: общие методы системного анализа, теория комплексирования разнотипных датчиков, общая теория случайных процессов и теория систем автоматического регулирования, методы цифровой обработки экспериментальных данных и методы математического моделирования.
Научная новизна. В процессе выполнения исследований были получены следующие новые научные результаты:

  1. Результаты расчетно-экспериментальных исследований временных и пространственных характеристик морского волнения.

  2. Предложена обобщенная математическая модель погрешности датчиков и модель схемы фильтрации выходного сигнала.

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

  4. Выведены аналитические выражения для определения значения спектра морского волнения на нулевой частоте, который представляет при исследовании точности наибольший интерес.

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



Практическая значимость. В диссертационной работе предложена модель интегрированной измерительной системы, позволяющая выполнять исследования в области повышения точности измеряемой высоты полета над морской и земной поверхностями, а также повышения уровня отказоустойчивости системы. Полученные при этом результаты дают возможность изучения детальных характеристик каждого датчика и системы в целом, изучения характеристик морского волнения, условий полета и качества функционирования датчиков, установленных на борту ЛА.

Значимость работы определяется возможностью использования полученных результатов при синтезе комплексированных измерителей малых высот движения для информационно-управляющих комплексов гидросамолетов и экранопланов при решении задач повышения точности измерения высоты, а также надежности системы измерения высоты полета и безопасности полетов.


Основные положения, выносимые на защиту:

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

  2. Алгоритм расчета и результаты анализа спектра погрешности локационного высотомера при полете на малой высоте над взволнованным морем.

  3. Методы синтеза интегрированной системы измерения параметров низковысотного полета с учетом возможности отказов первичных датчиков.

  4. Оптимизация интегрированной системы измерения параметров низковысотного полета по критериям точности и отказоустойчивости.


Внедрение результатов. Результаты диссертационной работы внедрены на предприятии ЗАО «НПФ «ТИРС» (Закрытое акционерное общество «Научно-производственная фирма «Телевидение и радиосвязь») при решении комплекса задач по повышению точности оценки измеряемых величин для разрабатываемых на предприятии радиосистем. Также результаты применялись при выполнении теоретических исследований по грантам РФФИ № 06-08-00550-а «Методы гарантирования приемлемой точности систем навигации и управления подвижными объектами», № 07-08-00293-а «Методы исследования нежестких структур и предотвращения критических режимов движения» (НИР № 5361-673-1, № 5361-678-1). Результаты работы использованы в Международном институте передовых аэрокосмических технологий ГУАП при создании экспериментальной системы записи параметров полета и параметров морского волнения для 12-метрового экраноплана. Кроме того, результаты диссертационной работы внедрены в учебный процесс Санкт-Петербургского государственного университета аэрокосмического приборостроения (по кафедре компьютерного проектирования аэрокосмических измерительно-вычислительных комплексов).

Внедрение результатов диссертационной работы подтверждено соответствующими актами, копии которых даны в приложениях к диссертации.


Апробация работы. Основные положения работы докладывались на следующих научно-технических конференциях и семинарах:

-XVI IFAC World Congress – 16-ый Всемирный Конгресс Международной федерации аэрокосмического управления, Прага, 2005;

-IV международная Конференция «Авиация и космонавтика», Москва, 2005;

-17th Symposium IFAC in Automatic Control in Aerospace – 17-ый Симпозиум по Автоматическому управлению в аэрокосмических системах Международной федерации аэрокосмического управления, Тулуза, Франция, 2007;

-XXXIV Всероссийская конференция «Управление движением морских судов и специальных аппаратов», Институт проблем управления им. В.А. Трапезникова РАН, Москва, 2007 г.;

-VII, VIII, IX Конференции молодых ученых ФГУП ЦНИИ “Электроприбор”, 2005, 2006, 2007 гг.;

-VIII, IX, Х Научные сессии ГУАП, 2005, 2006, 2007 гг.
Публикации. По результатам диссертационных исследований опубликовано 12 печатных работ, в числе которых:

-9 научных статей, одна из которых напечатана в журнале, рекомендованном для публикации ВАК России;

-3 публикации тезисов докладов.
Структура и объем диссертации. Диссертация состоит из введения, пяти глав, заключения, списка литературы (101 наименование) и двух приложений. Объем основной части – 150 страниц машинописного текста. Работа содержит 41 рисунок.


КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ ПО РАЗДЕЛАМ


Во введении обоснована актуальность и практическая значимость проводимых исследований, определена цель работы, основные задачи исследования, сформулированы научные и практические результаты, выносимые на защиту.
Первый раздел носит обзорный характер. В нем проведен обзор и сравнительный анализ локационных и инерциальных датчиков параметров низковысотного полета, даны рекомендации по выбору датчиков, поставлены требования к интегрированной системе. Обоснована актуальность применения интегрированной системы на современных экранопланах и гидросамолетах с учетом новейших направлений применения авиационной техники, представлены особенности расположения датчиков на ЛА таким образом, чтобы обеспечить точность, надежность и отказоустойчивость системы.

В настоящее время известны различные способы измерения высоты полета: радиотехнический, лазерный, акустический, радиоизотопный, барометрический, инерциальный. Если ЛА при полете не имеет непосредственного контакта с водной или какой-либо другой подстилающей поверхностью, то измерить высоту его движения можно лишь косвенно, путем оценки какой-либо другой величины, связанной с высотой определенной зависимостью. Более или менее упорядоченная зависимость естественных свойств воздушной среды от высоты в области высот, составляющих единицы метров, практически отсутствует. Следовательно, единственно возможным способом измерения геометрической высоты полета следует признать метод активной локации. Этот метод наиболее эффективен при совместной обработке показаний датчиков геометрической высоты, где наиболее пригодным является фазовый радиовысотомер, и инерциальных датчиков вертикального ускорения.

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

Минимальное количество пар датчиков высотомер-акселерометр равно трем. Их разнесение по корпусу ЛА дает возможность измерять углы крена и тангажа, и, как следствие, повысить точность определения координат (рис. 1).



Рис.1. Схема расположения датчиков на корпусе экраноплана


Исследуемая интегрированная система позволяет решить задачи измерения малых высот полета при ее использовании на современных ЛА, повысив точность показаний и обеспечив отказоустойчивость системы. Система может быть установлена на таких ЛА, как пожарный самолет-амфибия Бе-200ЧС, выполняющий аварийно-спасательные и поисковые операции на малых высотах, или на современных грузопассажирских экранопланах – скоростных низколетящих судах, использующих при движении эффект экрана. Одной из новейших разработок ЛА является экраноплан Buchon M-1, разработанный в UNEFA University в Венесуэле при техническом содействии Международного института передовых аэрокосмических технологий ГУАП.

Во втором разделе рассмотрены методы повышения точности измерений, проанализирована погрешность локационного высотомера и инерциального датчика акселерометра. Решена задача исследования влияния морского волнения на точность измерения. Временные и пространственные характеристики морского волнения необходимо учитывать в связи с тем, что при низковысотном полете погрешность измерения высоты ЛА относительно среднего уровня волновой поверхности практически полностью обусловлена отслеживанием высотомером профиля морских волн, а влияние собственной погрешности высотомера, связанной с особенностями его построения и с флюктуациями информационного параметра отраженного радиосигнала, оказывается пренебрежимо малым. Проведен анализ полученных выражений для спектральной плотности волнения, который позволяет найти наиболее неблагоприятные направления движения аппарата. Изучение уровня спектральной плотности на нулевой частоте приводит к задаче рассмотрения системы в установившемся режиме. В случае использования винеровской фильтрации необходимо знать весь спектр, а не только его уровень на нулевой частоте, но в рассматриваемом случае спектр практически равномерен, поэтому исследуется вопрос о том, какую реальную точность дает винеровская фильтрация с заданными спектральными плотностями воздействий.

Спектральные составляющие погрешности локационного датчика можно считать пренебрежимо малыми по отношению к спектру морского волнения. Ошибка в канале радиовысотомера будет практически полностью обусловлена динамическими искажениями за счет фильтрации отраженного сигнала, содержащего волновую составляющую.

Спектры морского волнения при использовании его трехмерной модели будут различаться не только положением на оси частот и шириной, но и формой. Сложность возрастает и с уменьшением скорости полета, так как при этом спектр погрешности становится более низкочастотным, а также со снижением высоты, когда пропадает эффект пространственного сглаживания волн в пределах угловой апертуры высотомера.

Для скоростных морских аппаратов максимум спектральной плотности на нулевой частоте при сильном морском волнении достигается при движении в направлении, близком к перпендикулярному относительно генерального направления распространения волн. Следовательно, такое направление в рамках поставленной задачи является наиболее неблагоприятным.

Для построения математической модели морского волнения необходимо произвести пересчет его спектральной плотности из неподвижной системы координат в движущуюся. Сложность такой задачи для трехмерной нерегулярной модели волнения обусловлена тем, что элементарные волновые составляющие с различными частотами распространяются не только с различными скоростями, но и в различных направлениях. Это обстоятельство требует определенного вклада в искомый спектр каждой элементарной составляющей волнения с последующим интегрированием этих вкладов по всем возможным частотам и направлениям распространения. Переход в движущуюся систему координат связан с частотным преобразованием вида
,
где V – горизонтальная скорость движения аппарата, g – ускорение силы тяжести,  – курсовой угол на волну,  – направление распространения элементарной волновой составляющей. Таким образом, в случае принятия гипотезы о взаимной независимости законов углового и частотного распределения мощности волнения справедлива формула:
,
где S1() – спектр морского волнения в движущейся системе координат; – спектр морского волнения в неподвижной системе координат, принятый в качестве стандартного; – закон углового распределения мощности волнения;
h3% – высота морских волн трехпроцентной обеспеченности.

При переходе в движущуюся систему координат за основу принимается спектр Пирсона-Мошковица. В идеальном случае подстилающая поверхность может являться горизонтальной плоскостью с зеркальным отражением, в таком случае в сигнале радиовысотомера присутствовала бы только когерентная составляющая, что свидетельствует о наличии в таком сигнале исключительно двух погрешностей: методической и инструментальной. В реальной погрешности датчика присутствует также некогерентная составляющая сигнала при отражении от подстилающих поверхностей. Получаемые параметры измерения и информационной параметр имеют флюктуирующий характер, что соответственно снижает точность измерения. Аппроксимировав полученное выражение, можно получить следующую приближенную формулу для спектральной плотности S1():



с параметрами



,



,

, ,
где ; ; ;

– уровень спектральной плотности на нулевой частоте, который может быть с приемлемой точностью рассчитан по методу прямоугольников при разбиении интервала интегрирования на 100 подынтервалов.

Даже при использовании оптимальной фильтрации, получение приемлемой точности обеспечивается только для некоторых режимов работы измерительной системы. Весьма важной информацией является уровень спектральной плотности погрешности радиовысотомера. Эту погрешность в пределах полосы пропускания можно рассматривать как белый шум. Чтобы найти предельное состояние фильтра для исследования точности системы, необходимо рассматривать уравнение фильтра в установившемся режиме. Такой фильтр можно трактовать как фильтр Винера.

Спектральные плотности полезного сигнала (измеряемой высоты) , ошибки высотомера и ошибки акселерометра (приведенной к полезному сигналу) составляют:
, , ,
где – параметры, входящие в корреляционные функции составляющих измеряемых параметров; , – интенсивность белого шума в каналах радиовысотомера и акселерометра.

В зависимости от соотношений между параметрами сигналов решения приобретают те или иные качественные особенности. Наиболее характерные варианты взаимного расположения кривых условных спектральных плотностей (при ), которые выполняются для ЛА, приведены на рис. 2.




1. , .



2. , .



3. , .





4. , .



5. , .






Рис. 2. Взаимное расположение кривых спектральных плотностей

Таким образом, необходимо учитывать только флюктуационные ошибки источников информации типа белых шумов. В результате решения системы N уравнений вида:


,
где - условные взаимные спектральные плотности выходных сигналов измерителей, соответствующие взаимно корреляционным функциям ; - условные взаимные спектральные плотности желаемого сигнала и сигналов измерителей; - неизвестные функции, не имеющие особенностей в верхней полуплоскости комплексной плоскости .

Проведенный анализ показывает, что одним из требований к системе должно быть выполнение условия инвариантности относительно полезного сигнала. Обеспечение конечности установившейся дисперсии ошибки является актуальным при решении задачи преобразования сигналов двух разнотипных датчиков. В исследованном установившемся режиме интервалы работы системы превышают по длительности наибольшие постоянные времени.


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

Любая структурная избыточность в измерительной системе, каждый дополнительный датчик расширяют возможности повышения качества измерений. Чтобы сделать сопоставимыми характеристики погрешностей датчиков различных типов, в том числе позиционных датчиков, целесообразно условно привести погрешность инерциального датчика к размерности перемещения. Собственную погрешность датчика можно представить как сумму погрешности типа медленного ухода нуля шкалы датчика, погрешности за счет отклонения коэффициента передачи датчика от номинального значения и погрешности вследствие неточности стабилизации оси чувствительности датчика . Будем считать, что сигнал ускорения на выходе датчика дважды интегрируется, в результате чего получается сумма полезного сигнала перемещения и погрешности с размерностью перемещения. Вторая производная такой погрешности есть погрешность на выходе датчика ускорения. Погрешность трех локационных высотомеров, разнесенных по корпусу ЛА, на низких частотах (конкретно в диапазоне примерно ) имеет равномерную спектральную плотность, максимальный возможный уровень которой составляет . При исследовании интегрированной измерительной системы наиболее важными являются спектральные плотности «геометрической» погрешности измерения высоты и погрешности собственно высотомера на частоте . Необходимо определить какой точностью должен обладать локационный высотомер для выполнения условия << . Считая погрешность экспоненциально коррелированным шумом со спектральной плотностью =, для дисперсии этой погрешности запишем формулу:


,
где – уровень спектральной плотности погрешности на нулевой частоте; – спектральная плотность «геометрической» погрешности измерения высоты; – практически максимальная погрешность; = 0,1с – постоянная времени радиовысотомера.

При пятибалльном морском волнении и скорости полета , уровень спектральной плотности волновой (геометрической) составляющей погрешности . Отсюда получим при условии, что уровень спектральной плотности погрешности собственно высотомера составлял не более 2% от . Полученный результат доказывает правомерность пренебрежения погрешностью собственно высотомера при исследовании результирующей погрешнрости.

Исследуя вторую и третью составляющие результирующей погрешности датчика вертикального ускорения, поперечные ускорения можно принять пренебрежимо малыми. Большое количество случайных факторов, определяющих характер изменения функций второй и третьей производных погрешности типа медленного ухода нуля шкалы датчика и , дает основание предложить, что их значения распределены по законам, близким к нормальному. Поэтому для дисперсий функций и можно принять соответствующие оценки: , .

При анализе дисперсий производных составляющих погрешности измерения параметров будем считать известными дисперсии функций: вертикального ускорения объекта , вертикальной скорости объекта ; угловой ошибки стабилизации оси чувствительности датчика (в вертикальной плоскости) ; скорости изменения этой угловой ошибки . При оценке количественных характеристик погрешностей будем использовать следующие численные значения: , , , .

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

Применительно к дисперсии первой производной этой погрешности получено: . В канале вертикальных перемещений погрешность необходимо дополнить составляющей, вызванной ускорением силы тяжести и имеющей квазипостоянную вторую производную с дисперсией . Третью производную можно считать практически нулевой, т.е. , а первую производную – неограниченно большой.

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

Поэтому ее дисперсия, учитывая взаимную независимость случайных процессов и , равна .

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

Таким образом, функцию можно представить в виде суммы двух взаимно независимых составляющих, у одной из которых ограничена только вторая производная, а у другой – вторая и первая производные. Представим дисперсию второй производной составляющей формулой:
.
Составляющая имеет дисперсию второй производной и первой производной .

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

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

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



Рис. 3. Погрешность интегрированного измерителя с механическим (2)

и микромеханическим (1) акселерометром
Результирующая погрешность при использовании погрешности микромеханических акселерометров превышают допустимую величину ошибок, следовательно, их использование при высокой балльности волнения (свыше 5 баллов) является нецелесообразным в интегрированных измерительных системах, поскольку не позволяет достичь требуемой точности измерения высоты.
В четвертом разделе рассмотрены принципы построения информационно-управляющих комплексов, рассмотрен вопрос об обеспечении нормального и эффективного функционирования интегрированной системы обработки сигналов измерителей с учетом имеющихся сведений о свойствах сигналов, предусмотрена автоматическая перестройка фильтров при наличии структурной избыточности. Цифровая реализация алгоритмов обработки измерений позволяет построить структурную схему системы, состоящую из цифровых и непрерывных динамических звеньев, определенным образом соединенных. Осуществляется выбор цены младшего разряда в каналах датчиков. Задача технического диагностирования должна не только устанавливать, что на борту есть неисправный датчик, но определять какой именно датчик отказал и перестраивать измерения таким образом, чтобы не использовать показания неисправного датчика.

При исследовании принципиальных возможностей того или иного варианта интегрированной системы по обнаружению отказов первичных датчиков и при определении максимальных ошибок измерения, накапливающихся к моменту обнаружения отказа, важен характер изменения погрешности отказавшего датчика. Можно контролировать три уровня точности измерений высоты полета в следующих условиях: 1. Все датчики исправны, обеспечивается максимальная точность измерения; 2. Один датчик отказал, но отказ не обнаружен – наиболее опасный случай; 3. Отказ датчика обнаружен, система реконфигурирована.

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

Рис. 4. Трехканальная система независимого вычисления

координат по сигналам шести измерителей

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



Если две из разностей, например и , превосходят пороговое значение и при этом третья разность мала, то отсюда следует, что ошибки второго измерителя недопустимо велики и что его показания, по-видимому, недостоверны. Устройство обработки данных в этом случае должно отключать измеритель и производить соответствующую коммутацию в основном блоке .

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



В пятом разделе приведены результаты математического моделирования интегрированной системы измерения параметров полета. Построены математические модели и исследованы переходные процессы в каналах датчиков, исследованы зависимость среднеквадратической ошибки измерения от направления движения и скорости ЛА. Проведено сравнение эффективности применения оптимальной и робастной фильтрации, приведены результаты математического моделирования при различных режимах работы датчиков, что позволяет сделать вывод о целесообразности применения исследованной интегрированной системы измерения параметров низковысотного полета.

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


,
где - коэффициенты, подлежащие определению.

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

Каждый из фильтров с дробно-рациональной передаточной функцией может быть реализован независимо от других фильтров, что позволит получить схему измерителя в виде параллельного соединения шести обособленных фильтров. Ширина разрядной сетки должна иметь величину 24 разряда. Учитывая единообразие знаменателей передаточных функций , целесообразным является реализовать вычислитель как единую схему фильтра с шестью входами, не разделенной на обособленные каналы и содержащей меньшее число элементов, чем схема фильтрации. Передаточная функция схемы в виде замкнутого контура имеет вид: . Фильтр должен выделять погрешность одного датчика из ее аддитивной смеси с погрешностью другого датчика .

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

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

Запишем подлежащие параметрической оптимизации передаточные функции каналов локационного датчика , канала инерциального датчика , а также передаточную функцию фильтров в виде


,

,

,
где .
Оптимизация параметров дала следующие результаты:



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

По результатам математического моделирования проведен сравнительный анализ ошибок измерений для исследуемой системы, результаты изображены на рис. 5. и на рис. 6.






Рис.5. Результаты моделирования интегрированной системы измерения высоты


Рис. 6. Зависимость среднеквадратической ошибки измерений от высоты волн трехпроцентной обеспеченности:

Где  - направление движения, град.;

е – среднеквадратическое отклонение, отн.ед.;



V – скорость полета.

1 – точность, которую дает винеровская фильтрация;

2 – точность, которую дает робастная фильтрация при , V=60м/с.



Суммарная ошибка измерений представляет собой низкочастотный случайный процесс, что обеспечивает сравнительно благоприятные условия функционирования аппарата, использующего такой сигнал в системе управления. Рассматривалось движение со скоростью от 10 м/с до 100 м/с в разных направлениях. Точностные характеристики измерительной системы определяются в основном уровнем спектральной плотности волнения на нулевой частоте. В связи с этим при низких скоростях движения наиболее неблагоприятным является направление, совпадающее с генеральным направлением распространения волн, а при высоких скоростях – направление, близкое к перпендикулярному относительно волны. Точность измерений в значительной степени зависит от взаимного расположения векторов скорости объекта и генерального направления распространения волн. Видно, что в диапазоне высот волн трехпроцентной обеспеченности ( м) система с постоянными коэффициентами обеспечивает более низкую точность по сравнению с требуемой (>0,3отн.ед.). Такой случай недопустим по критериям безопасности, и решением проблемы может быть наложение ограничений на направление движения аппарата или запрет полета в условиях повышенного морского волнения, поскольку это может привести к аварийной ситуации.

Исследование показывает, что робастный подход дает более точную оценку величин в широком диапазоне условий функционирования, чем в случае построения фильтра Винера при одинаковых условиях полета. Оптимальный фильтр, синтезированный для номинальной модели, теряет свойства при отклонениях от реальных характеристик, что подтверждает целесообразность использования робастной фильтрации при построении интегрированной системы измерения низковысотного полета.

Результаты математического моделирования в трех режимах полета (1 –все датчики исправны, обеспечивается максимальная точность измерений, 2 – отказ датчика обнаружен, система реконфигурирована, 3 – один датчик отказал, но отказ не обнаружен – наиболее опасный случай) изображены на рис. 7 при и моделировании движения ЛА в направлении .



Рис. 7. Результаты математического моделирования



в различных режимах работы датчиков
Наибольшая точность обеспечивается при наличии трех исправных датчиков на борту ЛА. Когда один датчик отказал, но отказ не обнаружен, полет является наиболее опасным и может привести к аварийной ситуации. В этом случае погрешность превышает максимально возможную, наблюдается довольно большой разброс параметров. В третьем случае, когда датчик отказал, отказ обнаружен и система реконфигурирована, показания стабильны относительно истинных, но дают меньшую точность измерения, что свидетельствует о целесообразности использования трех пар датчиков в исследуемой системе.
В заключении сформулированы основные результаты диссертационной работы.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ
Основные результаты диссертационной работы кратко можно перечислить в следующем виде:


  1. Выполнена общая классификация задач, относящихся к области измерения малых высот полета, представлен обзор существующих методов измерения параметров полета, алгоритмов автоматического обнаружения отказов датчиков, выделены основные проблемы и ограничения, характерные для задач обработки информации. Рассмотрены отличительные особенности построения интегрированных систем.

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

  3. Проведен анализ спектра погрешности локационного высотомера при полете на малой высоте над взволнованным морем, исследованы основные вероятностные свойства и спектрально-корреляционные характеристики морского волнения, решена задача синтеза и оптимизации системы по критериям точности и отказоустойчивости.

  4. Предусмотрена цифровая реализация интегрированной измерительной системы.

  5. Разработан подход к синтезу интегрированной системы, который обеспечивает гарантирование приемлемого качества оценивания и управления в условиях неполной априорной информации о погрешностях датчиков и возмущениях с учетом всего многообразия режимов движения летательного аппарата. Использование новых возможностей совершенствования алгоритмов навигации и управления движением вблизи поверхности, создаваемых современными средствами информационного обеспечения и ресурсами бортовых вычислительных сетей, может существенно способствовать развитию и повышению функциональной эффективности экранопланов, гидросамолетов и других летательных аппаратов.

  6. Решена задача оптимизации динамических характеристик фильтров, установлены критерии оценки точности, верхняя и нижняя границы невыхода параметров за допустимые пределы оптимальными фильтрами, что позволило значительно улучшить характеристики.

  7. Предложенная структурная избыточность обеспечивает повышение точности и надежности системы, а система обработки информации обеспечивает более быстрое обнаружение неисправностей и перестройку параметров.

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



СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ


  1. Небылов А.В., Румянцева Е.А. Исследование интегрированного измерителя высоты полета над взволнованной морской поверхностью //Материалы VIII Научной сессии, -СПб.: СПбГУАП, 2005. С. 121-124.

  2. Nebylov A.V., Rumyantseva E.A., Nebylov V.A., Danilov S.N. Wing-In-Ground Flight Automatic Control System //Proceedings of XVI IFAC Congress, Prague, 2005. Р. 247-252.

  3. Nebylov A.V., Rumyantseva E.A. Accuracy and reliability insuring for low altitude flight integrated measurеment system //Тез.докл. IV Международной Конференции «Авиация и космонавтика», -М.:МАИ, 2005. С. 113.

  4. Лопарев А.В., Румянцева Е.А. Синтез интегрированного измерителя малых высот для морских объектов //Материалы VII Конференции молодых ученых «Навигация и управление движением» под редакцией академика РАН В.Г. Пешехонова, -СПб.: ГНЦ РФ ЦНИИ «Электроприбор», 2005. С. 87-93.

  5. Румянцева Е.А., Греков А.И. Исследование интегрированной системы измерения низковысотного полета при наиболее неблагоприятных свойствах воздействий //Материалы IX Научной сессии, -СПб.: СПбГУАП, 2006. С. 72-74.

  6. Румянцева Е.А., Греков А.И. Исследование интегрированной системы измерения параметров низковысотного полета по критериям точности и отказоустойчивости //Тез.докл. VIII Конференции молодых ученых «Навигация и управление движением» под редакцией академика РАН В.Г. Пешехонова, -СПб.: ГНЦ РФ ЦНИИ «Электроприбор», 2006. С. 281.

  7. Румянцева Е.А., Греков А.И. Методы анализа отказоустойчивости интегрированной системы измерения параметров полета //Материалы IX Научной сессии. -СПб.: СПбГУАП, 2006. С. 24-26.

  8. Небылов А.В., Румянцева Е.А., Греков А.И. Робастные алгоритмы прогнозирования в задачах управления движением //Российская Академия Наук, Институт проблем управления им. В.А. Трапезникова//Материалы XXXIY Всероссийской конференции «Управление движением морских судов и специальными аппаратами», -М.: РАН, 2007. С. 26-30.

  9. Румянцева Е.А., Небылов А.В. Разработка методики синтеза интегрированной измерительной системы с учетом требований отказоустойчивости //Тез.докл. IX-ой Конференции молодых ученых «Навигация и управление движением» под редакцией академика РАН В.Г. Пешехонова, -СПб.: ГНЦ РФ ЦНИИ «Электроприбор», 2007. С. 25.

  10. Румянцева Е.А.. Разработка и исследование вариантов построения интегрированной системы измерения параметров полета на малой высоте //Материалы X Научной сессии ГУАП, -СПб.: СПбГУАП, 2007. С. 47-52.

  11. Nebylov A.V., Rumyantseva E.A. Comparative Analysis of Design Variants for Low Altitude Flight Parameters Measuring System //Proceedings of 17th Symposium IFAC in Automatic Control in Aerospace, -Tolouse, 2007. Р. 445-449.

  12. Румянцева Е. А. Оценка гарантированной точности //Известия высших учебных заведений. Технические науки. -Поволжский регион, 2008, №2. С. 115-122.








Раньше учтивый мужчина уступал женщине место в автобусе. Сегодня высшее проявление вежливости — уступить женщине свое рабочее место. Патрик О’Рурк
ещё >>