CC BY
35
УДК 517.977.1
А. Н. Рогалев
ВОПРОСЫ РЕАЛИЗАЦИИ ГАРАНТИРОВАННЫХ МЕТОДОВ ВКЛЮЧЕНИЯ ВЫЖИВАЮЩИХ ТРАЕКТОРИЙ УПРАВЛЯЕМЫХ СИСТЕМ
Описано применение гарантированных методов, позволяющих получать границы траекторий управляемых динамических систем, сохраняющихся (выживающих) до некоторого момента времени в пределах заданного множества. Такие количественные оценки являются границами множеств достижимости, удовлетворяющими условиям выживаемости. Приведены примеры расчетов границ выживающих траекторий.
Ключевые слова: гарантированные границы множеств достижимости, выживающие траектории.
Рассмотрим систему обыкновенных дифференциальных уравнений (ОДУ), в которую входят члены, описывающие параметрические возмущения, или неточности задания коэффициентов системы. Известны лишь границы, в которых они могут лежать. Среди математических описаний подобных задач можно выделить задачи проверки гарантированных условий безопасности и задачи построения множеств достижимости, т. е. задачи выживания (термин взят из [1]).
Задачи выживания для управляемых динамических систем включают в себя большое число вполне конкретных приложений, интерес к которым не ослабевает с конца 50-х гг. прошлого столетия. К числу таких прикладных задач относятся задачи об обходе препятствия, о построении управления, удерживающего траектории системы в заранее заданном множестве, в частности на заданном многообразии, некоторые задачи математической экономики и многое другое.
Поведение управляемого объекта выражается системой дифференциальных уравнений
^ = /(/,у,и), и еР, t е[^,0], to <0<да, (1)
где у - да-мерный вектор системы; и - управление; Р -компакт в евклидовом пространстве Рк. При этом предполагается, что выполнены следующие условия:
- функция /(/,у,и) непрерывна по совокупности
переменных в I х Рт х Р, а для любой ограниченной
замкнутой области В с Рт существует постоянная Ь = Ь(В) е (0, да), такая что
||/(/,у*,и) - /(Г,у*,и)|| < Ь||у. - у*||
V(^,у*, и), (/,у*,и) е I х Б х Р;
- существует постоянная це (0, да), такая что
||/(I, у, и)|| < ц(1 +|у||) V(t, у, и) е I х Рт х Р.
Допустимое управление и(Г) - измеримая по Лебегу функция, удовлетворяющая включению
и^) е Р, t е I.
Решением уравнения (1), порожденным допустимым управлением, является абсолютно непрерывная
функция у^), t е I, такая что — = /^, у(t), и^))
Л
почти всюду на I.
Множество достижимости У (^ определяется как множество всех точек из фазового пространства, в которые можно перейти на отрезке времени ^0,Т] из всех точек начального множества фазовых состояний М0 по решениям системы (1) с начальным условием уЦ0) е М0 и с допустимым управлением и^). Точное или приближенное описание множества достижимости управляемой системы позволяет оценить предельные возможности системы управления и выбрать оптимальное управление [2-6].
Включением (гарантированной оценкой) множества достижимости ОДУ У (/) называется многозначная функция У (Г), такая что для любого уЦ0) е М0 и t е ^0, Т] решение у(^ у0) е У^). Выбор возможных реализаций управляющих воздействий и стеснен ограничениями и е Q, t е [0 ,Т], отражающими особенности рассматриваемой задачи. Для многих задач ограничения на воздействия могут носить только геометрический характер. Это значит, что в каждый момент времени t е ^0, Т] значение и ^) может быть любым из некоторого выпуклого компактного множества Q. Задача гарантированного оценивания состоит в определении множества (или его границ)
У^) э и у^, t0, у°). (2)
у0еУ 0
Чтобы выделить выживающую траекторию у ^) е N при ^ <t < П,необходимо контролировать выполнение включения У(п) е М для любого движения у(-), исходящего из точек области допустимых
начальных позиций У0 и заданных множествах N, М при переборе всех управляющих воздействий, удовлетворяющих ограничению и е Q [6]. Таким образом
предполагается, что неопределенные факторы в математической модели управляемой системы не имеют вероятностного описания, а известны лишь с точностью до множеств, их содержащих. Одной из основ-
ных проблем теории управления, т. е. оценивания в условиях неопределенности, является проблема нахождения информационных множеств неопределенных динамических систем (2), а также подмножеств указанного семейства, состоящих из всех траекторий, удовлетворяющих дополнительному условию выживаемости. Таким образом, задача гарантированного оценивания состоит в описании множества решений системы (2) (траекторного пучка) и соответствующих сечений в момент времени t (сечений этого множества).
Выполнение гарантированных методов, основанных на аппроксимации оператора сдвига вдоль траектории [7-14], разделено на два этапа: предиктор и корректор.
На первом этапе (предикторе) происходит построение (запись) символьных формул приближенных решений как векторных функций Б" (У0) о Б"-1(У0) о... о £ '(У0), где вектор У0 - вектор начальных значений, рассматриваемых как символьные величины. Затем вычисляется область значений Бу этой формулы.
На втором этапе (корректоре) определяется гарантированная оценка глобальной ошибки приближенного решения Р(У0, г).
Выполнение преобразований над символьными формулами является началом построения гарантированных границ решений задачи, после которого нужно провести числовые расчеты. Некоторые системы компьютерной алгебры позволяют исходя из построенных формул генерировать программы на подходящем для таких расчетов языке программирования. Однако получаемые формулы, как правило, являются громоздкими и непосредственные вычисления с их помощью неэкономны, особенно если речь идет о вычислениях в цикле. Кроме того, предоставляемая системами компьютерной алгебры возможность проводить вычисления с числами большой разрядности делает арифметические операции дорогими, так что проблема экономии вычислений стоит в этом случае еще более остро.
Используя последовательные подстановки и приведение подобных членов, символьную формулу можно преобразовать в выражение, зависящее только от У0. При этом затраты памяти для хранения получающихся символьных формул будут крайне высокими.
Поэтому в общем случае для гарантированных методов предлагается алгоритм вычислений (преобразований и вычислений) символьных формул, основанный на поэтапном статичном хранении информации и ее преобразовании на завершающем этапе алгоритма. Таким образом, формула будет представлять рекурсивную структуру, размер которой изменяется. Для записи такой формулы в компьютере используются линейные динамические структуры. В силу этого модель вычислений (преобразований и вычислений) символьных формул осуществляется без явного вы-
писывания суперпозиций компонентов формулы, определяемых на каждом шаге. Связь между этими компонентами определяется посредством задания механизма адресации. Ссылки на адреса различных уровней хранятся в стековой памяти в виде дерева. Генерация кода вычислений по символьной формуле осуществляется в процессе обхода этого дерева, начиная с вершин.
В алгоритме получения символьных формул
У" = р" ^0 ^ У0 У1 У" ) =
= Б" (У0) о... о £ 2(У‘-1) о £ '(У1')
используется приведенная ниже методика обработки последовательности символьных формул. Пусть
ф(у0 ) - это однозначное отображение единичного интервала из Р1 на гиперкуб из Р", которое каждой точке t е Р сопоставляет точку у = ф^). При помощи такого отображения можно построить алгоритм исполнения, который для каждой точки t е Р позволяет
определить формулу отображения р(У10, У20,..., У°) и процесс ее сборки по адресам. Для этого в качестве отображения ф(у0) предлагается использовать непрерывное, однозначное отображение единичного интервала на "-мерный куб, известное также под названием кривой Пеано, заполняющей пространство. Фактически кривая Пеано представляет собой непрерывную, нигде не дифференцируемую кривую, которая проходит через все точки единичного гиперкуба в
пространстве Р". Изобразить кривую Пеано нельзя, можно лишь дать последовательность кривых, которая в пределе сходится к ней [15]. Каждая такая кривая называется приближением кривой Пеано и имеет номер, определяющий ее номер в последовательности кривых.
Таким образом, 2да-приближение можно рассматривать как некоторую аппроксимацию да-функции в рекурсивной формуле. Это соответствие задано отображением элементов конечного множества отрезков из единичного интервала и элементами конечного
множества гиперкубов, входящих в Р".
Такое соответствие строится следующим образом. Гиперкуб разбивается на 2" гиперкубов, называемых квантами первого разбиения. Длина ребра каждого такого кванта равна 1/2, кванты нумеруются индексом / от 2" -1 так, чтобы номера квантов, имеющих общую грань, отличались на 1. Кванты первого разбиения будем обозначать через г (/1), / = 2" -1.
Каждый квант первого разбиения таким же образом, как и гиперкуб в Р", разбивается на кванты второго разбиения с длиной ребра 1/4, которые нумеруются по тому же принципу, что и кванты первого разбиения. При этом нулевой квант второго разбиения, входящий в г, должен иметь общую грань с (2" - 1)-м
квантом второго разбиения, входящим в квант г (/ -1). Кванты второго разбиения обозначим через
г (/1, /2), где /2 - номер кванта второго разбиения, входящего в квант г (/1).
Аналогично получаем кванты любого разбиения с
номером г (/ , ^ , ..., 'да ), 0 < О < 2". Способ соединения
между собой квантов, номера которых отличаются на 1 в лексикографическом порядке, лежит в основе алгоритма, связывающего адреса, по которым хранятся компоненты символьной формулы.
Совокупность гиперкубов в пространстве Р" , рекурсивно квантованных да раз, носит название дискретного пространства да разбиения и обозначается Ода. Если в каждом разбиении число частей, на которые делится квант предыдущего разбиения, выбрать равным 2", то при любом да количество квантов в Ода1
и Ода" совпадет. Кванты да-разбиения Ода1 , полученные из одного кванта (да - 1)-го разбиения, нумеруются слева направо индексом 1да ,0 < 1да < 2" и обозначаются д(1г, /2,...,/да). Сопоставив кванты О1да и Опт с одинаковыми наборами индексов 11, /2,..., 1да, получим взаимно однозначное отображение фда : Р1да ^ Рпда, для отдельных квантов фда (д(1и /2,..., 1да)). Отображение ф да можно также интерпретировать как соответствие
между равноотстоящими точками из Р1 (центрами квантов да-разбиения) и узлами равномерной прямоугольной решетки в Р", где каждый узел также является центром кванта да-разбиения. Итогом работы описанного алгоритма будет возможность в любой точке ^ построить символьную формулу решения и получить на основе этой формулы численные значения решений.
Приведем примеры расчетов границ выживающих траекторий.
Первый пример - определение области при проти-воаварийном управлении.
При больших возмущениях в электроэнергетических системах (ЭЭС) возникают электромеханические переходные процессы, характеризующиеся взаимными движениями роторов синхронных машин в этих системах, существенными изменениями напряжений в узлах и токов в элементах ЭЭС. С электромеханическими переходными процессами связана проблема динамической устойчивости. Если динамическая устойчивость не нарушается, то в процессе электромеханических колебаний взаимные движения роторов синхронных машин происходят в некоторой ограниченной области. При нарушении динамической устойчивости разности углов, определяющих пространственное положение роторов синхронных машин, неограниченно возрастают. При развитии аварийных ситуаций в сложных ЭЭС имеют место длительные переходные процессы, связанные с отклонениями частот в отдельных частях этих систем. В случае не-
благоприятного развития таких ситуаций эти отклонения могут оказаться недопустимыми по условиям работы оборудования системы и потребителей. Подобные процессы часто происходят при системных авариях. С ними связано понятие живучести ЭЭС, определяемой как свойство системы противостоять возмущениям, не допуская их каскадного развития с массовым нарушением питания потребителей. При исследовании электромеханических переходных процессов система дифференциальных уравнений, описывающих динамику состояний ЭЭС, включает в себя уравнения следующих динамических элементов системы: движения роторов синхронных и асинхронных машин; электромагнитных переходных процессов в обмотках роторов: динамики регуляторов возбуждения и скорости синхронных машин; динамики регуляторов возбуждения и скорости синхронных машин; динамики регуляторов выпрямительных нагрузок и электропередач постоянного тока.
Границы множества решений в простейшей модели энергосистемы «станции - шины бесконечной мощности» [16] вычислялись при постоянно действующих возмущениях
= у2 ,
^2
= -16,52(у2 +1,134) 81п( у1 +1,256) + + 2,16Ып2( у1 +1,256) +16,05,
= -0,257у3 + 0,067 4 соб(у1 +1,256) + + 0,19и + 0,023 553 8.
Результаты расчета приведены ниже (рис. 1).
Рис. 1
Второй пример - расчет двойного маятника с общей управляемой точкой подвеса:
= у2 ,
dy2
12 = - у + у +и,
ёу1
= У4,
ёу±
Этот маятник можно рассматривать как систему из двух тел единичной массы, соединенных пружиной с единичным коэффициентом упругости и скользящих без трения вдоль неподвижного горизонтального стержня. Уравнения модели совпадают с линеаризованными уравнениями, описывающими поведение двойного маятника.
Управляющее воздействие и є Я прикладывается к левому телу для компенсации влияния внешних возмущений ^ є Я, воздействующих на правое тело. Возмущения предполагаются произвольными, но ограниченными в любой момент времени: |^(ґ)| < 1.
Обозначим через у1, у2 координату и скорость левого тела, а через у3, у4 - координату и скорость пра-
< Уі ^
У2 Уз
чУ 4 У
зового состояния данной динамической системы, полностью описывающей движение двойного маят-
вого тела. Тогда у =
R4 является вектором фа-
ника. В качестве выхода возьмем вектор
харак-
теризующийся величиной управления и координатой второго тела, на которой действуют внешние возмущения. Полученные в результате расчета включения множеств точных решений компоненты у2, у3 представим в графическом виде (рис. 2).
Таким образом, результаты применения методов, основанных на аппроксимации оператора сдвига вдоль траектории, дают хорошие оценки множеств решений, отклонение которых от точных решений стремится к нулю при уменьшении шага сетки.
Библиографические ссылки
1. Aubin J.-P. Viability theory. Boston ; Basel ; Berlin : Birkhauser, 1991.
2. Куржанский А. Б., Филиппова Т. Ф. Дифференциальные включения с фазовыми ограничениями. Метод возмущений // Оптимальное управление и дифференциальные уравнения : тр. Мат. ин-та Рос. акад. наук. М., 1995. Т. 211. С. 304-315.
3. Незнахин А. А., Ушаков В. Н. Сеточный метод построения ядра выживаемости для дифференциального включения // Журн. вычисл. математики и мат. физики. 2001. Т. 41, № 6. С. 895-898.
4. Незнахин А. А. О построении ядра выживаемости для обобщенной динамической системы // Прикл. математика и механика. 2006. Т. 70, № 5. С. 781-790.
5. Незнахин А. А., Ушаков В. Н. О построении ядра выживаемости для нелинейной управляемой системы при наличии целевого множества // Прикл. математика и механика. 2005. Т. 69, № 6. С. 976-985.
6. Куржанский А. Б., Филиппова Т. Ф. Об описании пучка выживающих траекторий управляемой системы // Дифференц. уравнения. 1987. Т. 23, № 8. С. 1303-1315.
7. Новиков В. А., Рогалев А. Н. Построение сходящихся верхних и нижних оценок решений систем обыкновенных дифференциальных уравнений // Журн. вычисл. математики и мат. физики. 1993. Т. 33, № 2. С. 219-231.
8. Рогалев А. Н. Исследование практической устойчивости при постоянно действующих возмущениях // Вычисл. технологии. 2002. Т. 7, № 5. С. 148-150.
9. Рогалев А. Н. Гарантированные методы решения систем обыкновенных дифференциальных уравнений на основе преобразования символьных формул // Вычисл. технологии. 2003. Т. 8, № 5. С. 102-116.
10. Рогалев А. Н. Гарантированные оценки безопасного функционирования технических и электроэнергетических систем // Современные методы математического моделирования природных и антропогенных катастроф : тр. Всерос. конф. с междунар. участием / Ин-т вычисл. моделирования Сиб. отд-ния Рос. акад. наук. Красноярск, 2003. С. 42-48.
11. Рогалев А. Н. Границы множеств решений систем обыкновенных дифференциальных уравнений с интервальными начальными данными // Вычисл. технологии. 2004. Т. 9, № 1. С. 86-93.
12. Рогалев А. Н. Символьные вычисления в гарантированных методах, выполненные на нескольких процессорах // Вестник НГУ. Серия : Информационные технологии. 2006. Т. 4, № 1. С. 56-62.
13. Рогалев А. Н. Вопросы устойчивости ансамблей дифференциальных уравнений // Вычисл. технологии. 2008. Т. 13, № 3. С. 111-117.
Рис. 2
14. Rogalyov A. N. Computation of reachable sets guaranteed bounds // Control, Diagnostics, and Automation (ACIT-CDA-2010) : Proc. of the IASTED Intern. Conf. on Automation, Control, and Inform. Technology. Calgary, Canada, 2010. P. 132-139.
15. Вирт Н. Алгоритмы + структуры данных = программы. М. : Мир, 1985.
16. Руденко Ю. И., Ушаков И. А. К вопросу оценки живучести сложных систем энергетики // Изв. АН СССР. Серия : Энергетика и транспорт. 1979. № 1. С. 14-19.
A. N. Rogalyov
IMPLEMENTATION OF CONTRIBUTED METHODS OF CONTAINMENT OF VIABLE TRAJECTORIES OF CONTROLLED SYSTEMS
The article describes implementation of contributed methods, which give the bounds of the trajectory of controlled systems, remaining (viable) to a certain time within a given set. Such quantitative estimations are bounds of reachable sets, satisfying the conditions of survival. Examples of computations of surviving trajectories bounds are presented.
Keywords: guaranteed bounds of the reachable sets, viable trajectories.
© Рогалев А. Н., 2011
УДК 551.510.534
Т. В. Рублева, В. Б. Кашкин, В. М. Мясников, С. А. Селиверстов
ИЗМЕНЧИВОСТЬ ПОЛЯРНОГО СТРАТОСФЕРНОГО ОЗОНА ЮЖНОГО ПОЛУШАРИЯ ПО СПУТНИКОВЫМ ДАННЫМ (АППАРАТУРА TOMS, OMI)
Рассмотрены пространственно-временные особенности весенней озоновой аномалии в Антарктиде (озоновой дыры). На основе спутниковых данных и с помощью корреляционно-экстремального алгоритма рассчитаны зональные скорости масс озона в южной полярной стратосфере. Приведены результаты расчетов площади озоновой дыры для области, где общее содержание озона (ОСО) меньше 220 единиц Добсона, и дефицита массы озона в стратосфере над Антарктидой. Обнаружена тесная связь дефицита массы озона в озоновой дыре и зональной скорости области накопления общего содержания озона в циркумполярном вихре за 1996-2009 гг.
Ключевые слова: истощение озонового слоя, динамика нижней стратосферы.
Исследование вариаций параметров озоносферы и поиски физических причин этих изменений до настоящего времени остаются проблемой, актуальность которой определяется существенной ролью озона в сложной системе «Солнце - атмосфера - поверхность Земли - космос». Поглощая ИК-излучение Земли в полосе 9,57 мкм, озон вносит вклад в изменение температуры атмосферных слоев, а поглощая УФ-излучение Солнца в диапазоне 200... 320 нм, он влияет на температурную стратификацию стратосферы, регулируя в этом слое планетарные тепловые и циркуляционные процессы.
С экологической точки зрения наиболее важен стратосферный озон, максимальная концентрация которого приходится на область высот 16.25 км (нижнюю стратосферу). Если собрать весь озон в слой вблизи поверхности Земли при давлении 760 мм рт. ст. и температуре 0 °С, то толщина этого слоя, т. е. общее содержание озона (ОСО), составит всего около 3 мм, или 300 единиц Добсона (ЕД) (1 ЕД = 10-5 м).
Несмотря на большое количество экспериментальных и теоретических исследований причины уменьшения общего содержания озона до сих пор не известны. По оценкам [1], до введения запрета на
производство озоноразрушающих веществ озоновый слой в северном полушарии истощался на 3,4 % за 10 лет, в Южном - на 4,8 %. После введения этого запрета в северном полушарии наблюдается уменьшение озона на 2,9 %, а в южном - на 1,0 %. Погрешность оценки скорости деградации ОСО не превышает 0,02 % за десятилетие. Эти оценки свидетельствуют о том, что уменьшение озонового слоя над северным полушарием происходит быстрее, чем над южным.
Проблема истощения озонового слоя обратила внимание на существующие межполушарные различия в сезонном ходе полярного стратосферного озона, особенно в весенний период [2]. Так, весной в северном полушарии ОСО снижается с 420 до 340 ЕД, а в южном - с 360 до 200 ЕД (рис. 1). Различия в значениях ОСО обусловлены особенностями циркуляционных процессов нижней стратосферы, а также климатическими изменениями в этих регионах.
Для Арктики характерна наибольшая межгодовая и внутрисезонная изменчивость суммарного озона. Сильное и устойчивое уменьшение ОСО в этом регионе началось в конце 1980-х гг., продолжалось до 1998 г. и сменилось периодом колебаний озона. Наи-
