автореферат диссертации по авиационной и ракетно-космической технике, 05.07.09, диссертация на тему:Методика идентификации характеристик лёгкого самолёта для мониторинга природных и техногенных катастроф

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

У

-ЬГ'

Буй Куанг Ли

□□34В1645

МЕТОДИКА ИДЕНТИФИКАЦИИ ХАРАКТЕРИСТИК ЛЁГКОГО САМОЛЁТА ДЛЯ МОНИТОРИНГА ПРИРОДНЫХ И ТЕХНОГЕННЫХ

КАТАСТРОФ

Специальность 05.07.09 Динамика, баллистика, управление движением летательных аппаратов

АВТОРЕФЕРАТ

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

1 2 ФЕБ

Москва-2008

003461645

Работа выполнена на кафедре «Системный анализ и управление» Московского авиационного института (государственного технического университета).

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

Официальные оппоненты:

Доктор технических наук, профессор, Баранов Вячеслав Николаевич.

Доктор технических наук, профессор, Дивеев Асхат Ибрагимович. Кандидат технических наук, доцент, Чернышев Андрей Валерьевич.

Ведущая организация:

ФГУП «НПО им. С.А.Лавочкина».

_часов на заседании

Защита диссертации состоится " дЬ/'Жми^ 2009 г. в . диссертационного совета Д.212.125.12 Московского авиационного института (государственного технического университета) по адресу: 125993, Москва, А-80, ГСП, Волоколамское шоссе, д.4.

С диссертацией можно ознакомиться в библиотеке Московского авиационного института (государственного технического университета).

Автореферат разослан «Z3 » мЩугт г.

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

В.В.Дарнопых

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТА

Актуальность темы диссертации

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

Объект диссертационного исследования - легкий самолет, предназначенный для мониторинга природных и техногенных катастроф.

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

Цель работы - выявление основных характеристик легкого самолета.

Задачи работы, решаемые для достижения поставленной цели:

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

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

3. Создание методики обработки информации на основе моделей факторного эксперимента.

4. Разработка модифицированных алгоритмов статистической обработки информации в виде обобщенных фильтров Калмана.

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

Методы исследования

В диссертационной работе основу исследований составляют теория вероятностей и случайных процессов, включая методы формирующих фильтров, методы оптимальной статистической обработки информации, современная теория оптимального управления летательными аппаратами, динамика полета и вычислительная математика пакета МАТЬАВ.

Научная новизна результатов диссертационной работы состоит в следующем:

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

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

3. Разработан безитерационный алгоритм формирования заданных режимов полета легкого самолета при летных испытаниях.

Практическая значимость результатов работы

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

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

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

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

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

Основные положения диссертационной работы, выносимые на защиту

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

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

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

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

Апробация работы и публикации

Результаты диссертационной работы докладывались на 12-й и 13-й Международной конференции «Системный анализ и управление». Крым, Евпатория, 2007 и 2008 , а так же на научном семинаре кафедры 604 Московского авиационного института (государственного технического университета).

Результаты работы опубликованы в 5-ых печатных изданиях, в том числе в 2-х изданиях из перечня, рекомендованного ВАК-ом Минобрнауки России.

Объем и структура работы

Диссертационная работа состоит из введения, 3 глав, заключения по работе, списка литературы и приложения. Работа содержит 141 страниц, 52 рисунков.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

dV 1 nV2

= — (Р cos a cos/?-С, z—S-G sin©); dt т "2

dpi nV1 nV1

m VK = (P sin a+C"yjx ^S) cos + (P cos a sin /?-СЦр S) sin y„ - G cos 0 ;

ич> nV2 * nV2

-mVK cos0-= (Psma+C"a——S)sinya-(Pcosasinfl-C? /3——S)cosy0\

dt 2 2

dav . .

1 + t-I ,)=My;

(1.1)

dip 1 / -V

— =--K cos у - со, sin y)\

dt cosS y

d3_

dt

- = 0)y sin у + Ú), cos y\

^ = a, - tg&{e>y cos у - юг sin у);

— = Frcos0cos4'; dt к

— = K.sin0; dtK

— = -Vr cosOsin^; dt K

sin © = cos a eos /í sin 5 - (sin а sin /5 eos /+sin /? sin /) eos 3;

sin 4* eos 0 = cos a eos físiny/ eos 9+sin a eos ¡3 (eos 4/ sin y+eos y sin 4/ sin S) -

— sin /? (eos v cosy - sin y sin 5 sin^);

sin eos 0 = eos a sin /?sin 5 - (sin «sin /? eos / - eos /3 sin y) eos <9,

где

m, G- масса и вес самолета соответственно; VK , V - земная и воздушная скорость; L , z - продольная и боковая дальность полета; у - высота полета; Р - тяга двигателя;

0 - угол наклона траектории к горизонту;

Ч1, у/ - угол пути и рыскания;

а, р - углы атаки и скольжения;

3 - угол тангажа;

у - угол крена;

уа - скоростной угол крена;

со,, ау, а, - угловые скорости;

1Х, 1у, 1г - моменты инерции;

Мх, Му, Мг - моменты сил;

С^ - коэффициент аэродинамической продольной силы;

' cí ' частные производные коэффициентов аэродинамической нормальной и поперечной силы по углу атаки и скольжения соответственно;

В качестве одного из основных атмосферных возмущений рассматривается ветер. Вектор скорости ветра характеризуется тремя проекциями на оси земной системы координат W(Wx,Wy,W2). Наличие ветра приводит к необходимости рассмотрения воздушной скорости и соответственно угла скольжения и дополнительного угла атаки - да„,.

Воздушная скорость определяется следующим образом:

K^v^v^vj, (1.2)

где

VK=V.cosO + Wx; VWr=V.sin Q + Wy\ F^ =F.cos0.sin V + W,. Угол скольжения и дополнительный угол атаки могут быть найдены по следующим приближенным соотношениям:

K2 + K2-Wr2cos24>~w2 Дат,,. = arceos-—;

2ГЛ

. V2 + V2cos2Q-W2-W^ .. „

0 = arceos—---1-(1.3)

2VrVKc os©

i

где Ув = {v,l + V^ + V¡ eos2 0 sin2 Ч*)2;

Исходные уравнения используются при моделировании возмущенного движения самолета. Уравнения движения, полученные на основе линеаризации,

используются при формировании алгоритмов идентификации аэродинамических характеристик легкого самолета.

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

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

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

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

1. Математическая модель движения самолета для решения задач идентификации.

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

(1(0: М6г.8 + М".а№ + АС .со,

Л = /.

а»

Необходимо отметить, что исследовалось «плоское» движение самолета для возможности наиболее четкой и простой организации летного эксперимента.

2. Предварительное определение коэффициентов аэродинамических моментов статической устойчивости и демпфирования т°, т"' на основе анализа переходных процессов.

Л к

Л к ¿о А,

(2.1)

Предварительное определение коэффициентов тт"• осуществляется на основе исследования переходных процессов с использованием линеаризованной модели продольного движения самолета.

+ ат Д Ук + ат Да + ам Д0 = -ат АЗ+а05 Хш;

ш

... с1гА9 с1А9 ■ ¿Да А . с1А8 ^15 (2.2)

°40д Кл+а42Да + а„Д© - = + а45К;

-Д5+Да + Д0 = О, где ал - динамические коэффициенты.

На основе модели (2.2) может быть получено уравнение, описывающее переходный процесс по углу атаки:

Да+а|1.Да+а|2.Да = а13.Д5, (2.3)

Где я„, ап динамические коэффициенты, определяемые формулами:

а»- 1г -—Г, у I ' {2А)

При этом величина периода колебаний Т и логарифмического декремента

затуханий £ связаны с динамическими коэффициентами следующими соотношениями:

Г = = —(2.5)

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

(2.6)

п a(tn )

где п - число колебаний, a(t0), a(tn) - значения угла атаки в моменты времени t0,

Таким образом, определение m", m"- осуществляется в следующей последовательности: на основе анализа переходного процесса определяются величина Т и по формуле (2.6) величина £. Далее с использованием выражений (2.5) находятся динамические коэффициенты оп, а12 и окончательно из соотношений (2.4) определяются m% .

В качестве иллюстрации были проведены расчеты применительно к гипотетическому малоразмерному самолету со следующими аэродинамическими характеристиками:

С* = 0,08; = 0,03; m" =-0,48; =-0,4; m6/ =-0,72. На рис.2.1 показано изменение a(t) по времени при угле отклонения руля высоты S = -10°.

— ■ 1— а |

Рис. 2.1. График изменения a(t) при угле отклонения руля высоты S = -10°

Анализ переходного процесса (рис.2.1) дает следующие результаты:

Л, =2,8;

© = 2л-/= 17,58 (рад 1с)\ ап= со1 =309,18;

-<^1 = 309,18;

< = '-^=-0.472; q-S.bA

а0 =24,050;

а4 = 17,660;

^Нй=о'о77:

я,, =2.^7^7 = 2,71;

<'=-^ = -0,413. q.S.bA

Необходимо отметить, что угол атаки непосредственно не измеряется, а может быть определен расчетным путем на основе обработки данных инерциальных измерителей. Однако в этом случае удается найти угол между вектором VK и продольной осью самолета, отличающийся от угла атаки на некоторую величину, обусловленную наличием ветра. В связи с этим дополнительно выполнялись расчеты коэффициентов т" и т"- с учетом возмущающего воздействия ветра. При этом использовалась процедура метода статистических испытаний (дисперсия скорости ветра принималась равной l,s[m2/c2]). Результаты расчетов выборочных

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

1 30

»ь =-0,473

—» i ,0

ш/=—.У <-(i) = - 0,412 30 tr

i 30 _

=0,00031

Таким образом, предлагаемое предварительное определение коэффициентов и т'*' может быть выполнено с достаточной степенью точности.

3. Предварительное определение аэродинамических коэффициентов Сх, С° .

Предварительные оценки Сх, С" получаются на основе обработки показаний акселерометров с помощью метода наименьших квадратов. Вектор измерений 1 имеет вид:

п.

2 =

(2.7)

Р í ^ г^а • \ РУ S

где л, =—+I-C, cosa + C,asina)---;

* G v ' y' ' 2 G

n =ÍCr sina + C°acosa)^-—. Г \ y. ' 2 G

Для возможности использования аналитической реализации метода

наименьших квадратов производится линеаризация уравнения измерений:

AZ = Ah.AX + и, (2.8)

или в развернутом виде:

ГДл

ДZ =

Дл„

д с:

. дпх рУ2 5 . , дпх . рУ2 Б п

дСх 2 С оСл 2 О

. рУ2 5 , рК2 5 Л,, =—— = зта---; и,, =—^ = асоза---,

дс2 0 5с;а 2 в

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

Оценки ДСх и ДС° получаются в соответствии с процедурой метода наименьших квадратов:

ДЛ- = (длг.дй)"' .Дйг.Д2, где (дйг.да) -дисперсионная матрица.

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

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

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

(2.10)

где /1=(IV/Ж; /2=dв/dt; Я = ¿¿М; /4=4-М; /5 = dP|dt; /6=Ау. /<й; /7=<М/<й.

Расширенный вектор состояния оцениваемой динамической системы имеет следующий вид:

х = Х-Х1КШ =[АУ а© Ду Ар Дй>2 А,9 Дот" Д/я°г ДС_ АС*]Г. (2.11) При этом для А/и/, дсг<) ДС£ используются уравнения формальных формирующих фильтров:

с/Дот," /Л = 0; ¿Д<- Дй = 0; </ДСг> /Л = 0; ¿ДС^ /Л = 0. В итоге получаем линеаризованную динамическую модель системы:

— = F(t)-x + Ви, dt

где

р/2 £

В - матрица (11x11) с отличным от нуля элементом = —--;

2 т

и - вектор управления с отличным от нуля элементом иь=8ь. Матрица коэффициентов линеаризации Б имеет следующий вид:

(2.12) (2.13)

Пх) =

/2 0 0 / 5 0 0 0 0 А. ю 0

fll /22 0 0 /25 0 0 0 0 0 /.11

/з, /,2 0 0 0 0 0 0 0 0 0

Л. /42 0 0 0 0 0 0 0 0 0

/я /52 0 0 /55 0 0 0 0 0 0

/«1 0 0 0 /б5 /« 0 /б8 /« 0 0

0 0 0 0 0 /76 0 0 0 0 0

0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 0 0 0

(2.14)

где

f^dfi/dv-f.^dfjd®-,.......f^dfljdm:- ;/U0=df,/dCXc -J^dfijdCl; (2.15)

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

(Z^L)

Z =

2г = У

(2.16)

На основе линеаризации уравнений (2.16) получаем линейную модель измерений:

(2.17)

(2.18)

Z(0 = h(t)-x(t) + u(t),

Матрица коэффициентов линеаризации имеет следующий вид:

'0 0 /г,, ООО 0 0 0 0 о л

0 0 0 А« 0 0 0 0 0 0 0

А(0 = К 0 0 0 Л,5 0 0 0 0 0

К 0 0 0 hi5 0 0 0 0 0 K\s

1° 0 ООО h* 0 0 0 0 0 J

где

К

dL l,>hi ду Ь

и ЛЬ-- и ^ 8V'^5 dp

az«

" dV ' 45 dp

Kk

h,„ =

az3

5Z4

эс:

az,

5or

2- = l.

Оптимальные статистические оценки для вектора состояния динамической системы (2.13) при использовании измерений (2.17) получаются на основе решения уравнений фильтра Калмана:

йх(1Ч)

dt dP{t) dt

- = F(l)x(t/t) + k(t)

Z(t)-h{l)x(t!t)

(2.19)

=F(0P(0+p(t)rmr-mmTR'[(№t)P(t),

где

m=p«)hT №-'(.»;

P(t) = £ )*(</'/]• - ковариационная матрица ошибок оценок.

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

В том случае рассматривалась наблюдаемая динамическая система вида:

Z„ = hx„+vn,

где

fl> = / + F-Ai.

Уравнение фильтрации рассматривались в следующем виде:

Кы = Ф*„-,/„-, + [Z„ - ^v,/,,-, ];

Р„,„ = [P;L + hTR-yh]~' = />„_, - Pn,„.,hT [hP„,„_,hT + Я]"' hP„M.

(2.20)

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

В качестве исследуемого режима полета легкого самолета был выбран режим набора высоты при ступенчатом отклонении руля высоты.

Интенсивность белых шумов измерений характеризуется диагональной матрицей Л, в которой:

Я = [о,25 1 КГ* 10"4 10"5].

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

Шя Р0=[1 О,ЗЛО"3 4 9 0.01 0,МО3 0,3.10^ 0,25.10'2 0,25.10~2 0,25.10"2 0,25.10""2]. Остальные элементы матрицы Р0 полагаются равными нулю.

Результаты моделирования в виде ошибок оценок ех и их дисперсий йх показаны на рис.2.2.

Рис. 2.2. Результаты моделирования процесса оптимальной обработки информации с помощью фильтра Калмана

Анализ представленных результатов моделирования показывает, что фильтр Калмана позволяет получить итоговые оценки характеристик легкого самолета т°, т"-, С^, С* с высокой степенью точности.

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

1. Решение задачи режгшов полета на основе принципа минимума Понтрягина

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

Математическая модель движения самолета представлена в форме (1.1).

Критерий оптимальности представляется в виде следующего оптимизирующего функционала:

J = R = [(¿. - if + {у. -у)г+(:. (3.1)

где L, у., z. - координаты самолета, соответствующие заданной программной траектории полета L, у, г- координаты самолета на возмущенной траектории реального полета.

Для решения поставленной задачи воспользуемся условиями оптимальности управления в форме принципа минимума. С этой целью составим гамильтониан системы:

„ ,r, , dV . d® , d4 . di . dy . <t ,,

/ = (3.2)

где ÁT = - сопряженный вектор.

Учитывая, что величина угла атаки на практике не превышает величину порядка 20°, можно положить Р sin а = Ра, Pcosa = Р. Кроме того, при формировании алгоритмов управления не будем рассматривать угол скольжения, который в дальнейшем будет учитываться при проверке полученных алгоритмов на основе моделирования процесса управления самолетом с учетом атмосферных возмущений, в частности, ветра.

Выделим в выражении (3.3) слагаемые, содержащие управление

,__L

■ mVA 2 I ° mVK eos©

Pa + C°a^S\svcira (3.3)

Структура оптимального управления определяется из условия минимума гамильтониана:

« = argimntf(A\A,fi), (3.4)

где X - фазовый вектор размерности (6x1), й - произвольное управление. Нетрудно видеть, что угол атаки входит в выражение (3.3) линейно, причем

оУг Л

P+C^^-S 1*0 т.е. особое управление отсутствует. Таким образом, оптимальное

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

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

о

Ш{уа) = Х1со&г.—уа. (3.5)

COS©

Введем обозначения

(3.6)

cos©

Из условия минимума выражения (3.5) по уа определяется структура оптимального управления . В этом случае фактически осуществляется минимизация скалярного произведения вектора Л (А,, Я,) и вектора-орта Kcos^.sin^).

Д =0;

+ 4>0; (3.7)

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

X = —, Х(10)~ задано дЛ

Л'=-—, Mt,) = —(3.8) дХ ' дХ '

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

2. Прямое решение задачи оптимизации управления (безитерационный алгоритм управления)

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

В связи с этим рассматривается следующая задача, эквивалентная исходной задаче оптимизации управления:

'/ >/ < min R = R0+ Jmin R{t)dt = R^ + R^ (tf -10 ) + J Jmin R(z, t)dzdt ;

r° '» 'o 'o r'

argmin(^fel). (3.9)

Использование второй производной позволяет получить соотношения, содержащие непосредственно угол крена и в конечном итоге определить . С этой целью составим выражения для R (/) и Л", (/), в которых в целях сокращения записи зависимость от времени не указывается:

^ ALAL + АуАу + АгАг И 10")

R '

где

AL = (L- L), Ау = (у. - у), Лг = (г.-г);

ДL = (L~L),Ay = (у.-у), А:' = (г.-г);

. Л(Д/.'2 + ALAL + Ау2 + АуАу + Л;2 + A='Az) - R(ALAL + АуАу + ДгДг)

R =-:-—--, (3.11)

где

A¿"=(¿- L), Ay = (у.; -/), А--' = (г; - г"). В выражении (3.11) от угла крена будет зависеть только часть этого выражения, содержащая вторые производные Д£", Ду", Дг". Эту часть в целях упрощения дальнейших выкладок будем рассматривать отдельно, обозначив следующим образом

R" (уа) = AL All + АуАу + Az Az (3.12)

В этом случае условие оптимальности (3.9) можно записать в виде

Ус,» = argmin(Ä^) (3.13)

г

Используя кинематические уравнения для L, у , : из (1.1), получим:

L = VK cos 0 cos «P - (^.0' sin 0.cos Ч' + ^Ч'зтЧ'cos0)

y' = VK sin0 + ©VA.cos0 (3.14)

= -VK cos 0sin4J + (Vl:Q sin 0sin cos 4* cos 0)

Используя динамические уравнения системы (1.1) для V', 0 , Ч* и подставляя их в (3.14), можно получить выражения для L, у", г", содержащие угол крена в явном виде:

(У.) - 003 У а sin 0 cos Ч* - Ду cos 0 - Дг sin 0 sin 40 (3.15)

- sin (At sin У + Дт cos 4*) Введем, как и ранее, следующие обозначения

А^ = (A£sin0cos4'-4ycos0-Azsin0sm4') В2 = -(AL sin У + Az cos Ч>) Тогда выражение (3.16) может быть представлено в виде, аналогичном (3.5), (3.6), т.е. можно рассматривать в виде скалярного произведения векторов Л(А2,В2) и ¡(eos у,sin/). Поэтому структура оптимального управления будет аналогичной (3.7).

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

У„._ =

-xsign(B2)+arcíg ß

arctg(-¡-),A2< 0.

4

— \,Л2 > 0; (3.17)

где

Аг =(£.-1)зт0созЧ/-(.у.-у)соз0-(2,-г^тОвтЧ*

(3.1 о)

В2 = - ¿^¡п У - (г. -г)соеЧ» Выражения (3.17), (3.18), как нетрудно видеть, фактически соответствуют решению задачи синтеза оптимального управления, формируемому на основе знания лишь фазовых координат. При этом, как показывают результаты расчетов, связанных с реализацией алгоритмов (3.7) и (3.17), численные значения выражений для (3.6) и (3.18) оказываются одинаковыми.

С динамической точки зрения алгоритм (3.17) обеспечивает максимальное значение проекции вектора подъемной силы на вектор Я(г), что фактически приводит к решению (3.9).

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

3. Результаты моделирования безитерационного алгоритма управления

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

- атмосферные возмущения в виде вариаций плотности атмосферы и ветра;

- случайные отклонения от расчетных значений начальных условий полета самолета;

- отклонения массово-инерционных и аэродинамических характеристик самолета от расчетных значений.

Вариации плотности атмосферы учитывались в следующем виде:

Р = РСТ( 1 + ДЙ, (3.19)

где - рст значение плотности, соответствующее стандартной атмосфере;

— Ар

- Ар = —i— относительное отклонение плотности от стандартного значения.

Рст

При моделировании возмущенного движения самолета использовались исходные уравнения движения (1.1). Результаты моделирования, характеризующие отклонения возмущенной траектории от номинальной при использовании алгоритма управления (3.17), показаны на рис.3.1.

V(m/c) О(рад)

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

условиях возмущений

- - полет по номинальной траектории.

......................- полет по возмущенной траектории с использованием алгоритма

управления (3.17).

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

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

4. Решение краевой задачи для канонической системы уравнений характеристик

Учитывая, что гамильтониан не зависит от дальности полета и бокового отклонения, а также пренебрегая зависимостью плотности атмосферы от высоты, можно считать что

где С4, С5, С6 — константы, подлежащие определению.

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

т.е.

(¡1 ¿1 ' Л ф ' Л <ь Я4 = С4,Л5 = С5Л = с6

(3.20)

(3.21)

(3.22)

Л

Таким образом, система уравнений характеристик имеет следующий вид:

—=собОСОЗЧ' ;

А* Л

Л

(3.23)

Л * '

Г—1-г:

а 'дч>\си) 6'9¥\ЛУ

где

dOKdt ) VK

8 (dV) ~ 0.sin0 (0 . na pV1 „, .

— — =4"« ---;— P.sma + C" .a.-—.S .siny„;

d®\dt J m.VK.cos2&{ >' 2 '

—f ■—I = Le = -VK .0. sin 0. cos 4; 8Q\dt J

—f—I = ¿,4- = -K..cos0.4''.sm4'; —i —| = r,f =-K,..cos0.lF'.cos4/.

5. Сравнительный анализ алгоритмов управления

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

1) Прямое решение задачи оптимизации управления в соответствии с алгоритмом (3.17).

На рис.3.4 в виде сплошной линии показано изменения величины расстояния между текущим положением самолета и заданной точкой на номинальной траектории при реализации алгоритма управления (3.17).

При расчетах использовались следующие граничные условия: ^(г0) = 50(л,/с) 0ft) = 1 (град.) Ч%) = 2 (град.) L(Í0) = SO(m) у('о)= 100(л<) Г(0 = 32(Л<) L(tj)=m(M) у(1,) = 2Щм) г(//) = 132(ж) t0 - начальный момент времени. tf - конечный момент времени.

2) Нахождение оптимального управления на основе решения двухточечной краевой задачи для системы уравнений характеристик (3.23), причем «начальные» значения сопряженных переменных Aj, Л3 принимались равными:

= (AL. sin 0. cos Ч" - Ay. cos © - A=. sin 0. sin 4>);

A, = -cos0(A£.sin4/ + Ar.cos(P).

Математическая модель движения _Х=ДХ,и)_

Оптимизирующий функционал J = R{t,u)

Структура оптимального управления и = arg min Н(X, Л, й)

"1 = К„ = . у„„ ШОЛО))

Прямое решение

(безитерационный алгоритм

управления)

Оптимизирующий функционал

Решение

=0;

-xsign(B2) + arag^j, Аг > 0

ß arctg(-f),A2< 0. Ä2

4 А

Двухточечная краевая задача с- ЗЯ

X = —,X(t„) - задано дЛ

дН_

' 8Х'

дХ

Функция штрафа

<7, =£,</,= = г; К-номер итерации

Декомпозиция сопряженного вектора

- Функции

МО

Л» - ~ - с6

времени

- константы,

определяемые методом деформируемого многогранника

Сравнение результатов моделирования двух алгоритмов управления

Рис. 3.2. Блок-схема решения оптимизационных задач управления

Сопряженные переменные в виде постоянных величин С4, С5, С6 находились

методом деформируемого многогранника. Для приведенных исходных данных

величины С4, С5, С6 оказались равными: С4 =0.23092, С5 = 0.066427, С6 =-0.060789.

На рис.3.3 показано изменение ключевых сопряженных переменных Хг и [м]

50

40 30 20 10

О 10

О 0.6 1 1.5 2 2.5 3

Рис.3.3. Изменение сопряженных переменных при решении краевой задачи

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

Как видно из рис.3.4 прямое решение достаточно хорошо совпадает с решением, получаемым на основе принципа минимума.

[м]

Рис.3.4. Изменение расстояния между самолетом и фиксированной точкой на номинальной траектории (сплошная линия соответствует прямому решению задачи оптимизации управления, пунктирная - решению на основе принципа

минимума)

ВЫВОДЫ ПО РАБОТЕ

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

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

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

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

СПИСОК ОСНОВНЫХ ПУБЛИКАЦИЙ

1. Буй Куанг Ли, Баранов В.Н. Комплексная методика определения движения и характеристик легкого самолета. Труды 12-й Международной конференции «Системный анализ и управление». Крым, Евпатория, 2007. - с. 56-57.

2. Буй Куанг Ли. Методика определения аэродинамических характеристик легкого самолета на основе результатов летных испытаний. Труды 12-й Международной конференции «Системный анализ и управление». Крым, Евпатория, 2007. - с. 57-58.

3. Буй Куанг Ли, Баранов В.Н. Алгоритм оптимального управления легким самолетом при летных испытаниях. Труды 13-й Международной конференции «Системный анализ и управление». Крым, Евпатория, 2008. - с. 112.

4. Буй Куанг Ли, Баранов В.Н. Оптимальное управление легким самолетом при летных испытаниях. «Вестник МАИ», № 4,2008. - с. 62-66.

5. Буй Куанг Ли, Баранов В.Н. Прямое решение задачи оптимизации управления пространственным движением атмосферных летательных аппаратов. Труды Института системного анализа Российской Академии Наук. Динамка неоднородных систем». -М: КомКнига, 2007, Том 31(2). - с. 7-13.

Введение.

Глава 1. Разработка и анализ обобщенных математических моделей движения самолета в возмущенной атмосфере.

1.1. Уравнения движения легкого самолета.

1.1.1. Уравнения движения центра масс.

1.1.2. Уравнения движения легкого самолета вокруг центра масс.

1.2. Линеаризация уравнений движения.

1.2.1. Общая система уравнений движения центра масс самолета.

1.2.2. Линеаризация уравнения движения.

1.3. Математические модели атмосферных возмущений.

1.3.1. Канонические разложения вариаций плотности и ветра.

1.3.2. Построение формирующих фильтров атмосферных возмущений.

1.3.3. Учет ветровых возмущений при испытаниях легкого самолета.

1.4. Выводы по главе 1.

Глава 2. Определение характеристик легкого самолета на основе результатов летных испытаний.

2.1. Основные методы обработки информации.

2.1.1. Фильтр Калмана и проблема расходимости процесса фильтрации.

2.1.2. Тактика обработки информации.

2.2. Предварительное определение характеристик легкого самолета на основе анализа переходных процессов

2.2.1. Математическая модель движения самолета для решения задач идентификации.

2.2.2. Предварительное определение коэффициентов аэродинамических моментов статической устойчивости и демпфирования т", т°г на основе анализа переходных процессов.

2.2.3. Предварительное определение коэффициентов А Сх0, А С".

2.2.4. Итоговое уточнение оценок на основе фильтра Калмана.

2.3. Выводы по главе II.

Глава 3. Формирование заданных режимов полета легкого самолета при летных испытаниях

3.1. Проблема реализации режимов и маршрута полета самолета при летных 83 испытаниях.

3.2. Решение задачи оптимизации управления на основе принципа минимума 84 Понтрягина.

3.3. Прямое решение задачи оптимизации управления (безитерационный алгоритм управления).

3.4. Результаты моделирования безитерационного алгоритма управления.

3.5. Решение краевой задачи для канонической системы уравнений характеристик

3.6. Сравнительный анализ алгоритмов управления.

3.7. Техническая реализация алгоритма оптимального управления.

3.8. Выводы по главе III.

Актуальность темы диссертации

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

- информация о лесных пожарах.

- поиск заблудившихся в лесу (джунглях), на море и при наводнениях.

- обнаружение и исследование объектов с вредными выбросами.

- наблюдение и поддержка при спасательных работах.

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

Вопросам отработки самолетов и их систем управления в процессе летных испытаний посвящено большое количество работ, частности [8], [17], [36], [43]. Однако, в большинстве из них в качестве объекта экспериментальных исследований рассматриваются пилотируемые самолеты различного назначения. Проблемы летной отработки малых беспилотных самолетов в литературе освещены менее значительно [38].

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

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

Объект диссертационного исследования - легкий самолет, предназначенный для мониторинга природных и техногенных катастроф.

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

Цель работы - выявление основных характеристик легкого самолета.

Задачи работы, решаемые для достижения поставленной цели:

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

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

3. Создание методики обработки информации на основе моделей факторного эксперимента.

4. Разработка модифицированных алгоритмов статистической обработки информации в виде обобщенных фильтров Калмана.

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

Методы исследования

В диссертационной работе основу исследований составляют теория вероятностей и случайных процессов, включая методы формирующих фильтров, методы оптимальной статистической обработки информации, современная теория оптимального управления летательными аппаратами, динамика полета и вычислительная математика пакета МАТЬАВ.

Научная новизна результатов диссертационной работы состоит в следующем:

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

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

3. Разработан безитерационный алгоритм формирования заданных режимов полета легкого самолета при летных испытаниях.

Практическая значимость результатов работы

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

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

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

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

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

Основные положения диссертационной работы, выносимые на защиту

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

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

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

4. Результаты численного тестирования разработанного в среде МАТЬАВ программно-математического обеспечения, реализующего предложенные алгоритмы оптимальной статистической обработки информации и управления легким самолетом.

Апробация работы

Результаты диссертационной работы докладывались на 12-й и 13-й Международной конференции «Системный анализ и управление». Крым, Евпатория, 2007 и 2008 , а так же на научном семинаре кафедры 604 Московского авиационного института (государственного технического университета).

Объем и структура работы

Диссертационная работа состоит из введения, 3 глав, заключения по работе и списка литературы. Работа содержит 141 страниц, 52 рисунков.

ВЫВОДЫ К ГЛАВЕ III

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

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

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

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

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

98

ЗАКЛЮЧЕНИЕ

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

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

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

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

100

1. Адлер Ю.П., Маркова Е.В., Грановский Ю.В. Планирование эксперимента при поиске оптимальных условий. - М.: Наука, 1971, 284 с.

2. Бажинов И.К., Почукаев В.Н. Оптимальное планирование навигационных измерений в космическом полете. М.: Машиностроение, 1976, 288 с.

3. Баранов В.Н., Михалев A.A. Методика применения датчиков угловых скоростей Analog Devices в системе управления малогабаритных летательных аппаратов. Труды 11-й Международной конференции «Системный анализ и управление». Крым, Евпатория, 2006

4. Баранов В.Н., Горбатенко Д.С. Динамическое проектирование системы управления сверхлегкого дистанционно пилотируемого вертолета. Труды 10-й Международной конференции «Системный анализ и управление». Крым, Евпатория, 2005, стр.87, 88

5. Баранов В.Н., Гавриков В.Г. К задаче обработки и информации при управлении спуском в атмосфере. Космические исследования АН СССР, Т. ХУ1, выпуск 4, 1978, 0.8 п.л.

6. Баранов В.Н., Гавриков В.Г. Приближенное решение задачи фильтрации при управлении полетом КА в атмосфере. Космические исследования АН СССР, Т. ХУШ, выпуск 4, 1980,1 п.л.

7. Берестов Л.М., Поплавский Б.К., Мирошниченко Л.Я. Частотные методы идентификации летательных аппаратов. М., Машиностроение, 1985. -184 с. ил.

8. Берестов Л.М., Горин В.В. Моделирование динамики управляемого полета на летающих лабораториях. М.: Машиностроение, 1988. 110 с.

9. Болтянский В.Г. Математические методы оптимального управления. М.-Наука. 1967.

10. Брандин В.Н., Васильев A.A., Куницкий A.A. Экспериментальная баллистика космических аппаратов. М.: Машиностроение, 1984, 264 с.

11. Буй Куанг Ли, Баранов В.Н. Комплексная методика определения вижения и характеристик легкого самолета. Труды 12-й Международной конференции «Системный анализ и управление». Крым, Евпатория, 2007, с. 56-57.

12. Буй Куанг Ли. Методика определения аэродинамических характеристик легкого самолета на основе результатов летных испытаний. Труды 12-й Международной конференции «Системный анализ и управление». Крым, Евпатория, 2007, с. 57-58.

13. Буй Куанг Ли, Баранов В.Н. Алгоритм оптимального управления легким .самолетом при летных испытаниях. Труды 13-й Международной конференции «Системный анализ и управление». Крым, Евпатория, 2008, с. 112.

14. Буй Куанг Ли, Баранов В.Н. Прямое решение задачи оптимизации управления пространственным движением атмосферных летательных аппаратов. «Труды ИСА РАН, Динамка неоднородных систем», Том 31(2), 2007, с. 7-13.

15. Буй Куанг Ли, Баранов В.Н. Оптимальное управление беспилотным легким самолетом при летных испытаниях. «Вестник МАИ», №5, 2008, с. 6266.

16. Буй Куанг Ли, Малышев В.В., Баранов В.Н., Ли Вэй. Методика определения аэродинамических характеристик малоразмерного самолета на основе результатов летных испытаний. «Вестник МАИ», 2009

17. Васильченко К.К., Леонов В.А., Пашковский И.М., Поплавский Б.К. Летные испытания самолетов. Учебник для студентов высших учебных заведений. М.: Машиностроение, 1996 - 720 с.

18. Войтишек A.B. Основы метода Монте-Карло в алгоритмах и задачах. -Новосибирск, НГУ, 1997.

19. Голубев И.С., Янкевич Ю.И. Основы устройства, проектирования, конструирования и производства летательных аппаратов. М.: Изд-во МАИ, 2006. - 528 с. ил.

20. Гришин A.M. Физика лесных пожаров. Томск: Изд-во Томского университета, 1994.

21. Джонсон. Н., Лион. Ф. Статистика и планирование эксперимента в технике и науке. Методы обработки данных. М.: Мир, 1980, 610 с.

22. Долженков H.H., Янкевич Ю.И. Современные комплексы БЛА и перспективы их развития. Вестник воздушного флота - 2004, № 10.

23. Ермаков С.М., Жиглевский A.A. Математическая теория оптимального эксперимента. М.: Наука, 1987 - 320 с.

24. Ермаков С.М., Михайлов Г.А. Статистическое моделирование. М.: Наука, 1982.

25. Ершов A.A. Стабильные методы оценки параметров (обзор). Автоматика и телемеханика. 1978, № 8. с 66-100.

26. Иванов Н.М., Лысенко Л.М., Мартынов А.И. Методы теории систем в задачах управления космическим аппаратом. М.: Машиностроение, 1981, 254 с.

27. Калман, Бьюси. Новые результаты в линейной фильтрации и теории предсказания. Техническая механика. Т.83, Серия Д, 1961, № 21.

28. Кринецкий Е.И (ред.). Летные испытания ракет и космических аппаратов. М.: Машиностроение, 1979 - 464 с.

29. Лебедев A.A., Красильщиков М.Н., Малышев В.В. Оптимальное управление движением космических летательных аппаратов. М.: Машиностроение, 1974, 199 с.

30. Лебедев A.A., Чернобровкин Л.С. Динамика полета беспилотных летательных аппаратов. М.: Машиностроение, 1973, 615 с.

31. Лебедев A.A., Бобронников В.Т., Красильщиков М.Н, Малышев В.В. Статистическая динамика управляемого полета. М.: Машиностроение, 1978, 240 с.

32. Леондес К.Т. Фильтрация и стохастическое управление в динамических системах. М.: Мир, 1980.

33. Летов A.M. Динамика полета и управление. М.: Наука, 1969.

34. Р. Ли. Оптимальные оценки, определения характеристик и управление. -М.: Наука, 1966.

35. Махонькин Ю.Е., Павлова З.А., Фальков А.И., Корачков В.И. Автоматизированная обработка результатов измерений при летных испытаниях. М.: Машиностроение, 1983,112 с ил.

36. Р. Мехра. Идентификация и адаптивная фильтрация Калмана. Механика (сборник статей), 1971, №3, с.34-51.

37. Миронов А.Д., и др. Методы исследований на летающих моделях. М.: Машиностроение. 1988, 144 с.

38. Остославский И.В., Стражева И.В. Динамика полета (траектории летательных аппаратов). М.: Машиностроение, 1969, 499 с.

39. Пугачев В.Н. Теория случайных функций и ее приложение к задачам автоматического управления. М.: Физматгиз , 1960, 883 с.

40. Рачкулик В.И., Ситникова М.В. Отражательные свойства и состояние растительного покрова. Л.: Гидрометеоиздат, 1981.

41. Соболь И.М. Численные методы Монте-Карло. М.: Наука, 1973.

42. Сухорученков Б.И., Меньшиков В.А. Методы анализа характеристик летательных аппаратов. М.: Машиностроение, 1995, 368 с.

43. Д. Химмельблау. Прикладное нелинейное программирование. М.: Мир, 1975.-534 с.

44. Шкадов Л.М., Буханова P.C., Илларионов В.Ф. Механика оптимального пространственного движения летательных аппаратов в атмосфере. М.: Машиностроение. 1972.

45. Школьный Е.П., Майборода Л.А. Атмосфера и управление движением летательный аппаратов. Л.: Гидрометеоиздат, 1973. - 307 с.

46. Шли, Стендиш, Тода. Расходимость фильтрации по методу Калмана. -Ракетная техника, № 6, 1967.

47. ГОСТ 20058-80. Динамика летательных аппаратов в атмосфере. Термины, определения и обозначения.