Исследование области притяжения нелинейной системы в условиях интервальной неопределенности | Статья в сборнике международной научной конференции

Библиографическое описание:

Галикеева Ю. Р. Исследование области притяжения нелинейной системы в условиях интервальной неопределенности [Текст] // Технические науки в России и за рубежом: материалы Междунар. науч. конф. (г. Москва, май 2011 г.). — М.: Ваш полиграфический партнер, 2011. — С. 48-51. — URL https://moluch.ru/conf/tech/archive/3/723/ (дата обращения: 20.08.2018).

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

Рассмотрим систему дифференциальных уравнений, описывающих продольное движение ЛА [1].

(1)

где V – скорость самолёта истребителя, m – его масса,, – проекции суммы внешних силы, действующих на ЛА, и тяги на соответственно оси Ox, Oy связанной системы координат, – проекция момента внешних сил на ось Oz, – угловая скорость относительно оси Oz, – угол тангажа, – момент инерции относительно оси Oz.

Для разложения в ряд Тейлора уравнений (1) необходимо знать аналитическое выражение величин , , . Данные зависимости сложны и определяются из теоретических выводов и опытных испытаний [1-3]. Результаты опытных испытаний имеют графическое представление [2], для определения их аналитических зависимостей предлагается применить аппроксимацию полиномами Чебышева [4]:

,

, (2)

,

где – искомая аналитическая зависимость аппроксимированной функции определённой на интервале ;

– полиномы Чебышева,

– дискретные значения функции, полученные из графических зависимостей

– сглаженные дискретные значения функции,

x – область определения искомой функциональной зависимости.

Для численных расчётов, были использованы характеристики на самолёт МИГ 29.

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

Рис. 1. Аппроксимация зависимости дополнительного коэффициента лобового сопротивления от числа Маха ортогональными полиномами Чебышева

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

. (3)

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

, (4)

где .

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

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

В общем случае под областью притяжения &#; системы (3) понимают наибольшую область начальных отклонений x(0), при которых решения этой системы являются асимптотически устойчивыми, или, иначе, устойчивые решения, удовлетворяющие условиям

(5)

принадлежат области притяжения нулевого положения равновесия системы (3).

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

. (6)

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

. (7)

Дальнейшее исследование предполагает построение интервальной матрицы вида

, , (8)

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

где , ; (9)

В противном случае

где . (10)

Задача нахождения минимума или максимума вещественной функции решается двумя методами: метод Хука и Дживса и метод генетического алгоритма.

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

  1. Генерация начальной популяции. Начальная популяция задаётся экспертным путём.

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

  3. Двухточечный кроссинговер. Случайным образом выбираются две точки разрыва i и j из диапазона [1, n-1]. Пусть при этом i<j. Тогда потомки определяются путем обмена соответствующих частей родительских хромосом:
    H1 = {c11, c12,.., c2i, c2i+1,.., c2j, c1j+1,.., c1n},
    H2 = {c21, c22,.., c1i, c1i+1,.., c1j, c2j+1,.., c2n}.

  4. Мутация. Вероятность 10 %

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

Область притяжения для ЛА МИГ 29, описывается неравенствами

рад/с, м/с, (11)

рад, рад

Применив метод ГА и метод Хука-Дживса, получены следующие результаты:

Матрица А и В, найденная методом Хука-Дживса:


Матрица А и В, найденная методом ГА:

Результирующие матрицы:


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

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


Литература:
  1. Бочкарёв А.Ф., Адреевский В.В. Аэромеханика самолёта–М: Изд-во Машиностроение, 1985.–357с.

  2. Учебное пособие, Практическая аэродинамика самолёта МИГ-29 –М., 1995.– 316с.

  3. Д. Гроп. Методы идентификации систем – М., 1978.

  4. Остославский И. В. Аэродинамика самолёта – М., 1957. – 560с.: ил.

  5. Исследование сложных технических систем на основе анализа областей притяжения / Р. Р. Акчурин, В. Н. Ефанов // 10-я Международная конференция по компьютерным наукам и информационным технологиям CSIT’2008: сб. науч. ст., Анталия, 2008, Т. 3. с. 131-136

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

Похожие статьи

Некоторые соображения о корректности и точности линейной...

Исследование области притяжения нелинейной системы... Рис. 1. Аппроксимация зависимости дополнительного коэффициента лобового сопротивления от числа Маха ортогональными полиномами Чебышева.

Генетический алгоритм для нахождения коэффициентов...

Библиографическое описание: Пожарский Д. А., Золотов Н. Б., Семенов И. Е. Генетический алгоритм для нахождения коэффициентов аппроксимации функции в контактных задачах

Таким образом, необходимо найти значения коэффициентов , при которых будет наименьшей.

Определение наилучших вероятностей мутации генетического...

Поиск особей с наилучшими значениями функции полезности производится по следующему алгоритму: Шаг 1. Запуск алгоритма. Создание начальной популяции.

Алгоритм интервального оценивания параметров нелинейных...

m — число параметров функции-модели; n — число точек регрессии; step0, step — начальный и рабочий шаг исходного многогранника (симплекса)

alarm — Булевская переменная, принимающая значение True в случае, когда количество вычислений ЦФ в одном цикле (I1)...

Математическая модель популяции, подверженной промыслу

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

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

Исследование функций преобразования емкостных уровнемеров...

Особенностью применения ЕУ в системах управления и контроля является малая начальная емкость датчиков, которая в большинстве случаев лежит в пределах от 10 до

Значение этого сопротивления зависит от качества изоляционных материалов

В общем случае, когда...

Применение метода линейного программирования для решения...

В этой точке функция принимает максимальное значение, равное 127.

Строим вектор целевой функции , координаты которого равны коэффициентам целевой

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

Применение генетического алгоритма для решения задачи...

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

Многомерная интерполяция сеточной вектор-функции

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

В области Dn определена сеточная вектор-функция Fc(xc) с размерностью m, т.е. заданы значения m компонент функции Fc в узлах xc

Обсуждение

Социальные комментарии Cackle

Похожие статьи

Некоторые соображения о корректности и точности линейной...

Исследование области притяжения нелинейной системы... Рис. 1. Аппроксимация зависимости дополнительного коэффициента лобового сопротивления от числа Маха ортогональными полиномами Чебышева.

Генетический алгоритм для нахождения коэффициентов...

Библиографическое описание: Пожарский Д. А., Золотов Н. Б., Семенов И. Е. Генетический алгоритм для нахождения коэффициентов аппроксимации функции в контактных задачах

Таким образом, необходимо найти значения коэффициентов , при которых будет наименьшей.

Определение наилучших вероятностей мутации генетического...

Поиск особей с наилучшими значениями функции полезности производится по следующему алгоритму: Шаг 1. Запуск алгоритма. Создание начальной популяции.

Алгоритм интервального оценивания параметров нелинейных...

m — число параметров функции-модели; n — число точек регрессии; step0, step — начальный и рабочий шаг исходного многогранника (симплекса)

alarm — Булевская переменная, принимающая значение True в случае, когда количество вычислений ЦФ в одном цикле (I1)...

Математическая модель популяции, подверженной промыслу

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

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

Исследование функций преобразования емкостных уровнемеров...

Особенностью применения ЕУ в системах управления и контроля является малая начальная емкость датчиков, которая в большинстве случаев лежит в пределах от 10 до

Значение этого сопротивления зависит от качества изоляционных материалов

В общем случае, когда...

Применение метода линейного программирования для решения...

В этой точке функция принимает максимальное значение, равное 127.

Строим вектор целевой функции , координаты которого равны коэффициентам целевой

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

Применение генетического алгоритма для решения задачи...

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

Многомерная интерполяция сеточной вектор-функции

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

В области Dn определена сеточная вектор-функция Fc(xc) с размерностью m, т.е. заданы значения m компонент функции Fc в узлах xc

Задать вопрос