Полный текст

Введение В этой статье мы продолжаем исследование задач управления для непрерывно-дискретных (гибридных) динамических моделей с последействием и полиэдральными ограничениями на управление, существенно используя результаты предыдущих работ [1-5]. Для модели, рассмотренной в работе [5], на основе полученного в ней описания множества достижимых значений целевых функционалов предлагаются конструкции и алгоритмы построения управляющих воздействий, позволяющих достигать заданных значений целевых показателей, если они принадлежат множеству достижимости. О роли множеств достижимости в теории управления можно судить по многочисленным публикациям, из которых мы здесь упомянем работы [6-15]. Следует отметить, что в подавляющем большинстве работ достижимость понимается по отношению к значениям координат фазового вектора в конечный момент времени. Для задач управления моделями экономических систем характерно более широкое понимание целей управления, при котором эти цели задаются с помощью конечной системы целевых показателей, в число которых входят не только локальные, сосредоточенные в отдельных точках, но и интегральные характеристики поведения системы на всем промежутке управления в целом. Динамика системы управления описывается совокупностью дифференциальных уравнений с запаздыванием и разностных уравнений с дискретным аргументом. Используемое в настоящей работе описание процессов (см. работу [5]) мотивируется наличием взаимодействия переменных, имеющих различный характер изменения - непрерывный (непрерывное производство) и дискретный (финансирование). Для систем с непрерывным временем и кусочно-постоянным запаздывающим аргументом, определяющим дискретную составляющую памяти входящих в уравнения операторов, подробное обсуждение макроэкономических моделей можно найти в работе [15] и цитированных в ней источниках. 1. Постановка задачи Напомним кратко данное в работе [5] описание рассматриваемой экономико-математической модели взаимодействия производственной подсистемы, описываемой уравнением (1) и финансовой подсистемы, описываемой уравнением (2) Здесь фазовая переменная зависит от непрерывного времени , а фазовая переменная зависит от дискретного времени Конкретный характер взаимодействия определяется заданными матричными коэффициен-тами Управления и реализуются с помощью интегральных операторов с ядрами и соответственно. Функции и интерпретируются как внешние воздействия на систему, например непредвиденные потери или возможные погрешности моделирования. Начальное состояние системы (1)-(2) считается заданным: (3) Цель управления задается с помощью определенного на компонентах и вектор-функционала (4) и представляет собой совокупность N равенств (5) где - вектор заданных целевых значений. Постоянные -матрица -матрицы и -матрица с измеримыми и ограниченными элементами считаются заданными. Общая форма (4) вектор-функционала со значениями впозволяет охватить разнообразные конкретные случаи целевых условий, возникающих в прикладных задачах. Соответствующие примеры приведены в работе [5]. Задача достижения целей (5) решается в предположении, что управления и стеснены ограничениями (6) где - заданная матрица размерностью предполагается, что множество решений системы линейных неравенств непусто и ограничено. Цель этой статьи - предложить для задачи (1)-(6) конструкции и алгоритмы построения управлений и , приводящих к выполнению целевых условий для достижимых значений показателей . 2. Основные конструкции Определим основные пространства для траекторий и управлений. Пусть - пространство суммируемых функций с нормой (- норма в ); - гильбертово пространство функций с нормой, порожденной скалярным произведением (- символ транспонирования); - пространство абсолютно непрерывных функций с нормой Зафиксируем множество и обозначим через пространство функций с нормой . Дальнейшие построения этой работы используют представление общего решения системы (1)-(2) в операторной форме [3]: (7) Здесь - фундаментальная матрица, - оператор Коши, компоненты оператора Коши являются ограниченными операторами. Предполагается, что операторы также являются ограниченными операторами. Обозначим далее . (8) Объединим компоненты и в общий вектор управлений Как показано в работе [5], целевые условия (5) приводятся к виду (9) где - моментная -матрица с элементами из . В терминах этой матрицы дается описание множества достижимых целевых значений в задаче (1)-(6). Приведем для этой матрицы определяющие соотношения, воспользовавшись представлением общего решения (7). Используя представление (4) целевого вектор-функционала представление (7) и обозначения (8), получаем где для -мерного вектор-столбца через обозначен i-й -мерный вектор-столбец. Отсюда следует, что в блочном представлении моментной матрицы блоки и определяются равенствами (10) и (11) после выполнения всех операций в правых частях и приведения подобных под знаками интегралов. Проведем в качестве примера соответствующие преобразования для равенства (10). Для этого в случае первого слагаемого воспользуемся представлением оператора и его свойствами [3]. Имеем Обозначая через выражение в фигурных скобках, получаем окончательно для первого слагаемого в равенстве (10): . (12) Проводя аналогичные преобразования, для второго слагаемого в равенстве (10) получаем (13) где - элементы матричного представления оператора - характеристическая функция отрезка . Таким образом, из равенств (12), (13) следует представление Напомним теперь некоторые обозначения из работы [5], используемые для описания множества достижимых значений целевых показателей. Для каждого вводится функция (14) Далее для каждого и каждого фиксируется совокупность точек множества доставляющих значение функционалу и вводится функция (15) В таком случае множество достижимых значений вектора состоит из тех и только тех точек для которых неравенство (16) выполняется для всех Дадим описание алгоритма построения программных управлений и с использованием построенной моментной матрицы 1. Задание конечного набора векторов, определяющих при каждом вместе с матрицей градиент целевой функции в задаче (14). 2. Задание разбиения основного промежутка [0,T] точками 3. Для фиксированного последовательно в точках решается задача 4. На каждом промежутке значение условно оптимального (оптимального по направлению кусочно-постоянного управления определяется равенством В случае неоднозначности операции argmax (экстремальное значение функционала достигается на грани многогранника) используется конструкция (15). 5. На основном промежутке [0, T] k-е условно оптимальное управление определяется равенством . 6. Для каждого вычисляется По построению, каждое является достижимым значением для оно достигается на построенном в п. 5 кусочно-постоянном управлении 7. Совокупность точек определяет выпуклый многогранник. Любая точка этого многогранника может быть представлена в виде выпуклой линейной комбинации точек 8. Управление под действием которого достигается целевое значение определяется равенством Первые компоненты этого -мерного управления определяют компоненты управления , остальные - компоненты управления 3. Иллюстрирующий пример Построим программные управления для системы (15)-(16) работы [5]: (17) (18) с нулевым начальным состоянием (19) и целевыми условиями (20) Управления и стеснены следующими ограничениями: (21) Приведем представления компонент оператора Коши для системы (17)-(18): где - характеристическая функция отрезка В рассматриваемом случае для нахождения элементов моментной матрицы следует выполнить замену и подставить выражения для компонент решения в компоненты заданного в этом примере целевого вектор-функционала (20). Действуя так, приходим к следующим выражениям: Таким образом, для элементов моментной матрицы имеем: Далее с помощью описанного выше алгоритма для и находятся достижимые значения компонент целевого вектор-функционала, они показаны на рис. 1. При этом черным цветом выделены три точки, образующие треугольник, содержащий заданный вектор целевых значений (точка с координатами (10; 6)), и указаны номера соответствующих направлений градиента при поиске условно оптимальных управлений. Выделенным точкам соответствуют следующие управления: (здесь и ниже индекс компонент управления соответствует номеру направления оптимизации); ; . Рис. 1. Достижимые и заданные значения целевых показателей Для нахождения искомого управления остается найти коэффициенты выпуклой линейной комбинации вершин упомянутого треугольника для внутренней точки (10; 6). Приведем координаты вершин треугольника с точностью до 0,001: (18,607; 6,760), (2,653; 4,634), (5,847; 7,168). Упомянутые коэффициенты выпуклой линейной комбинации определяются равенствами Таким образом, поставленная задача решается управлением где Рис. 2. Графики компонент построенного программного управления Графики компонент этого управления представлены на рис. 2, где более тонкая линия соответствует компоненте

Об авторах

В. П Максимов

Пермский государственный национальный исследовательский университет

Список литературы

Статистика

Просмотры

Аннотация - 55

PDF (Russian) - 47

Ссылки

• Ссылки не определены.