CC BY
59
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
УДК 629.
№ 284 Декабрь 2004
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
В.В. Кашковский, И.И. Тихий, Ю. Н. Шишкин
ИДЕНТИФИКАЦИЯ ПАРАМЕТРОВ ДИНАМИЧЕСКОЙ МОДЕЛИ ЛЕТЧИКА ПО ДАННЫМ БОРТОВОГО УСТРОЙСТВА РЕГИСТРАЦИИ
В настоящее время летная подготовка осуществляется по жесткой программе. Для повышения эффективности летного обучения необходим переход к индивидуальным планам летной подготовки. Осуществление индивидуального плана возможно только при условии объективной оценки качества пилотирования.
Индивидуальный план летной подготовки предполагает объективную оценку качества пилотирования. Существующий оперативный контроль не позволяет в полной мере решить эту задачу. Главный недостаток оперативного контроля - невозможность прогноза. То есть результатом оперативного контроля является сообщение о нарушении летного ограничения после того как оно произошло.
Вопросы оценки качества пилотирования активно исследовались в институте космической медицины [1]. В 1976 и 1978 гг. академик Н.И. Фролов доказал, что одной из характеристик качества работы летчика является структура движения ручкой управления и время активной работы стабилизатора на посадке. Однако данная работа выполнена в условиях научно-исследовательского института с применение специального исследовательского оборудования. Применить результаты этой работы для практической работы в летных школах и училищах невозможно. Поэтому предлагается разработка методики и программного обеспечения параметрической идентификации
модели летчика для оценки качества пилотирования летчиков в летных школах и училищах по данным бортовых устройств регистрации (БУР) для совершенствования индивидуальной летной подготовки.
Использование объективной информации БУР для оценки качества пилотирования является обязательным слагаемым летной подготовки, но возможности БУР на этом далеко не исчерпаны. Поэтому разработка новых алгоритмов оценки качества пилотирования остается актуальной задачей. Одному из таких подходов и посвящена данная статья.
Рассмотрим пилотирование самолета по глиссаде в ди-ректорном режиме. Для самолета фронтовой авиации продолжительность этого участка полета составляет примерно 100 с. При пилотировании по глиссаде летчик крайне ограничен в свободе своих действий, так как допустимое отклонение от нее не превышает нескольких метров по высоте и курсу. Это позволяет предельно упростить его модель в контуре пилотирования самолетом. Возможная обобщенная модель [2, 3] представлена на рис. 1.
Рис. 1
В [1, 2] доказана возможность оценки качества пилотирования при заходе на посадку по отклонению ручки летчика в продольном канале управления само-
летом. Исходя из этого и модели на рис. 1 рассмотрим структурную схему, представленную на рис 2.
тренд
8В
8с
ДЗ
Д©
3 £ 5РУСтг кЮ: ю2а (т© +1) —► 1 1
г Р 2 + 2^аЮа Р +юа р г т©+1
БУР
Алгоритм идентификации
оценка
Рис. 2
Тренд на рис. 2 - это реализация мысленного сигнала траекторного управления самолетом по Д© или концептуальная модель действий летчика. Выходным сигналом ^л(р) является отклонение ручки летчика по тангажу
канале управления Жл(р) является форсирующим звеном.
Отклонение 8Р
измеряется датчиком и регистрируется
8р
Исследования [3] показывают, что в продольном
БУР для обработки на земле. С помощью цифрового фильтра с частотой настройки юа сигнал 8РУС^ можно разделить на тренд и высокочастотную составляющую
ю
2
тг
тг
тг
8рустг . Не участвующая в траекторном управлении высокочастотная составляющая отражает процессы перерегулирования при реализации концептуальной модели. Высокочастотную составляющую сигнала 8РУС^ будем называть
исполнительской моделью действий летчика. Очевидно, что перерегулирование при исполнении концептуальной модели является одним из показателей качества пилотирования, которое статистически можно связать с уровнем подготовки данного летчика. Количественно перерегулирование может быть описано через параметры Wл(p). Возникает задача оценки параметров №л(р) по высокочастотной составляющей 8РУСтг.
Сложность данной задачи заключается в том, что мы можем наблюдать выходной сигнал звена !¥л(р) на
ограниченном участке времени —н и в ограниченном частотном диапазоне. Нижняя граница частотного диапазона равна юа верхняя определяется частотой регистрации 8РУСг и для большинства самолетов составляет 2 Гц. Входной сигнал №л(р) нам неизвестен, так как вычисленный тренд является таковым только условно. Анализ полетных данных показывает, что спектр случайной составляющей выходного сигнала !¥л(р) имеет множество локальных максимумов, на фоне которых трудно выделить максимум, равный частоте собственных колебаний Жл(р) (рис. 3), кроме того, выходной
О
сигнал 8рус„ имеет функцию плотности распределения, не являющуюся Гауссовой (рис. 4).
Рис. 3
Для количественного выражения степени отклонения плотности распределения дискретных измерений выходного
О
сигнала 8рус1Т от нормального закона распределения предлагается использовать алгоритм, основанный на свойствах срединной и средней арифметической ошибок нормального закона распределения. Для этого исходную последовательность измерений X(г) = {х(1),х(2),...,х(п -1), х(п)} сигнала
8 рустг преобразуем в последовательность У (г) = {у(1), у(2),.... ...,у(п -1),у(п)}, где у(г) = |х(г)-дах|. Затем последовательность У (г) ранжируем по возрастанию, и вычисляем коэффициенты а и Ь:
у(п /2 + 1)п ь = У(п /2 + 1)
Рис. 4.
/V О ~г'^ О
К 8 руСтг (т)=2|5 8 русТг (/)со8(2п/т)С/ , (1)
где 5 8 руСтг (/)= Иш — т„ ^да —
спектр мощ-
ности 8русТг (т), вычисляемый методом дискретного преобразования Фурье (ДПФ).
Из-за малого времени наблюдения оценка (1) имеет мультипликативную погрешность и функционально зависит от времени наблюдения —н. Исходя из этого,
оценку К 8русТг (т) предлагается вычислять через нор-
а =-
0,8453£ У (г)
г = 1
0,67451
Е У (г)2
г=1
п-1
мированную корреляционную функцию г 8 русТг (т)
Ко (т)= Б О Го (т) =
8 РУС тг 8 РУС тг 8 РУС тг
+ да
I Ф О (/ )со§(2 п/т ) С/ (2)
В соответствии с любой реализацией 8рус1Т можно поставить числовое значение коэффициентов а и Ь, количественно выражающих отклонение плотности распределения дискретных измерений этой реализации от нормального закона распределения Спектр и гистограмма выходного
О
сигнала 8русТг являются информативными, но при этом трудно установить их взаимосвязь с параметрами Жл(р).
Были выполнены предварительные исследования свойств статистической оценки корреляционной функции, полученной по теореме Винера - Хинчина:
= Б о
8 РУС
где Ф о
э 8руСтг
IФ о (/) й/
0 э 8Р___
—н О
| 8русТг (¿)в-j 1%/сИ
энергетический спектр
О 1 8русТг , вычисляемый методом ДПФ; БО =-х
8русТг п -1
хЕ
г=1
8 РУСтг (г—д )
оценка дисперсии вычисленная по
+ да
2
реализации 8русТг . Анализ полетных данных показал, что график корреляционной функции (2) имеет устойчивые отличительные признаки, позволяющие количественно характеризовать динамические свойства летчика. На рис. 5 - 6 показаны нормированные корреляционные функции 2-х реализаций посадочного маневрирования, выполненных летчиком.
1
0,75
0,5
0,25
0
-0,25
-0,5
л д --
V V V 'и 1 ’ / V 'V
ч
-10 -5
0
Рис. 5.
1
0,75
0,5
0,25
0
-0,25
-0,5
-10
-5 0
Рис. 6.
5
--—
5
Однако, как и в случае со спектром выходного сигнала остается неизвестная зависимость корреляцион-
О
ной функции сигнала 8 рус„ с параметрами модели летчика Wл (р).
Таким образом, возникла практическая задача идентификации параметров ^мфл, ^мл, |мл модели летчика:
К 2п^л (Р + 2п^мфл
Wмп (Р)=
Рмфл (Р2 + 4п^млРмл Р + 4п2
22
мл
по корреляционной функции К 8 русТг (т). Задача усложняется тем, что форма К 8русТг (т) зависит не только от параметров модели летчика, но и от частоты регистрации 8РУСг и частоты настройки фильтра юа.
Решение этой задачи предлагается методом имитационного моделирования [4]. Для выполнения исследований предлагается три имитационных модели летчика (рис 7 а, б, в), далее называемых цифровыми, чтобы отличать их от модели летчика Wл (р).
В начале были выполнены предварительные исследования модели б (рис. 7), в которой летчик представлен колебательным звеном. Цель исследований - выбор такого имитационного входного сигнала цифровых моделей, при котором спектр и коэффициенты а и Ь входили в диапазон отклонений аналогичных парамет-
О
ров, полученных по множеству реализаций 8 мрусТг .
Рис. 7
Исследования модели а (рис. 7), где летчик представлен усилительным звеном, позволили выявить детерминированную взаимосвязь экстремумов производной корреляционной функции К 8 русТг (т) с частотой
О
дискретизации 8 рус„ и частотой настройки фильтра
юа. Графически эту зависимость можно выразить в виде «сетки» опорными точками которой являются нули
и максимумы производной К 8русТг (т) (рис. 8).
Наилучшие результаты были достигнуты при использовании входного сигнала, плотность распределения дискретных измерений которого показана на рис. 9.
Исследования показывают, что «сетки» экстремумов производной корреляционной функции К 8 русТг (т) модели летчика, полученной по реальным полетным данным и цифровых моделей летчика б и в не совпадают с сеткой экстремумов модели а. Это позволяет предположить, что данное расхождение зависит только от параметров модели летчика Wл(p).
а
Рис. 8
0,095 0,00 0,085 0,08 0,075 0,07 0,065 0,06 0,055 0,05 0,045 0,04 0,035 0,03 0,025 0,02 0,015 0,01 0,005 0
-60 -40 -20 0 20 40 60 60
Рис. 9. Плотность распределения входного сигнала цифровых моделей и нормальный закон распределения
Целью дальнейших исследований цифровых моде- «сетки» цифровой модели а и «сетки», полученной по
лей б и в (рис. 7) является количественное выражение результатам обработки полетных данных.
параметров ^мфл, ^мл, |мл через разность опорных точек
ЛИТЕРАТУРА
1. Фролов Н.И. Пути изучения работоспособности летчика в полете // Косм. биол. 1978. № 1. С. 3-10.
2. Авиационная медицина: Руководство / Под ред. Н.М.Рудного, П.В.Васильева, С.А. Гозулова. М.: Медицина, 1986. 580 с., ил.
3. Ефремов А.В., Оглоблин А.В., Предтеченский А.Н., Родченко В.В. Летчик как динамическая система. М.: Машиностроение, 1992. 336 с.
4. Информационное и математическое обеспечение задач управления организационно-техническими системами: Учебн. пособие // Под ред.
В.Я. Головина. М.: ВАТУ, 2002. 239 с., ил.
Статья представлена кафедрой эксплуатации авиаоборудования Иркутского военного авиационного инженерного института (г. Иркутск), поступила в научную редакцию «Кибернетика» 15 августа 2004 г.
