Полный текст
1. Введение. Волновые уравнения, возникающие в задачах управления распределенными колебательными процессами, имеют неизменный теоретический интерес и существенное практическое значение (см. [1–5, 7–9, 11–18, 21, 22]). В задачах математического моделирования часто возникает необходимость генерации желаемой формы колебания или стабилизации колебания. Решение данной проблемы реализуется исследователями, как правило, с помощью задач граничного управления (см. [2, 4, 5, 7–9, 11–18, 21, 22]). Благодаря многочисленным приложениям многоточечные краевые задачи управления и оптимального управления динамикой являются активно развиваемым направлением в современной теории управления. В этих задачах, наряду с классическими краевыми условиями (начальными и конечными), дополнительно заданы многоточечные условия в промежуточные фиксированные моменты времени.
Задачам управления (в том числе и оптимального) динамикой разнородных составных систем посвящены, в частности, работы [1, 2, 4, 5, 7, 8, 11–18, 26]. Применительно к распределенной колебательной системе, включающей два кусочно однородных участка, эта задача была впервые сформулирована А. Г. Бутковским и исследована в [12]. Серия работ академика В. А. Ильина (см., например, [7, 8]) и работы [4, 5, 11, 13–16, 26] посвящены проблемам граничного управления (оптимального управления) процессами, которые моделируются одномерным волновым уравнением, состоящим из двух участков с разными физическими свойствами. Длины таких участков выбирались исходя из предположения, что время прохождения колебаний по каждому из них является одинаковым. Авторами указанных работ были изучены и выведены формулы типа Даламбера, при этом задачи исследовались методом бегущих волн.
В данной статье рассматривается серия задач граничного управления и оптимального граничного управления динамикой распределенной неоднородной колебательной системы, причем в промежуточные моменты времени известны различные состояния колебательного процесса, который состоит из двух кусочно однородных участков. Считаем, что физические характеристики этих участков удовлетворяют сделанным выше предположениям. Будем осуществлять управление и оптимальное управление за счет смещения одного конца (при закрепленном противоположном конце), а также за счет одновременного смещения обоих концов с заданными условиями: в начальный и конечный моменты времени, а также в разные определенные промежуточные моменты времени. Критерий качества в задачах оптимального граничного управления задан на всем интервале времени.
Сформулированные в данной работе задачи отличаются от существующих постановок тем, что помимо стандартных краевых условий заданы дополнительно многоточечные условия в промежуточные моменты времени, а именно: на функции колебания (прогиба), на их производную (функции скоростей точек), а также одновременно на функции колебания и производную функции колебания. При исследовании этих задач используется метод разделения переменных (метод Фурье).
Цель данной работы состоит в создании аналитического подхода и разработке алгоритма построения функции граничного управления и оптимального управления одномерными колебательными системами, обладающими неоднородными свойствами, динамика которых (под действием сформированного закона управления) за конечный отрезок времени переходит из определенного начального состояния через многоточечные промежуточные состояния в известное (желаемое) конечное состояние.
2. Постановка задачи. Пусть кусочно однородная среда состоит из двух участков с соответствующими длинами и (т.е. , ), скорость прохождения волны по -му участку, где линейная плотность, модуль Юнга, . При этом имеет место равенство
(2.1)
так что время прохождения волны по участкам разной длины совпадает. Пусть состояние неоднородной распределенной системы описывается функцией , , , а отклонения от состояния равновесия можно представить в виде волнового уравнения следующего вида:
(2.2)
с граничными условиями двух видов:
1. (2.3)
<p >2. (2.4)
Функции и управляющие воздействия (граничные управления) с условиями сопряжения в точке соединения участков
(2.5)
Распределенный кусочно однородный процесс (2.2) можно охарактеризовать как динамическую систему переменной структуры (см. [19]). Уравнение (2.2) характеризует математическая модель продольных (либо поперечных) колебаний стержня (струны) соответственно, где плотность, модуль упругости (натяжение струны).
Пусть классические условия (начальные и конечные) имеют вид
(2.6)
(2.7)
Пусть также в некоторые определенные моменты времени :
известны промежуточные значения функции колебания (прогиба) и значения ее производной (скоростей точек системы) в следующем виде:
(2.8)
(2.9)
(2.10)
В промежуточных условиях (2.10) предполагается, что четное число.
Замечание 1. Промежуточные значения функции колебания и значения производной функции колебания в условиях (2.10) можно задавать в любой очередности.
Для задач управления с граничными условиями (2.3) будем рассматривать функционал вида
(2.11)
а для задач управления с граничными условиями (2.4) будем рассматривать функционал вида
(2.12)
Предполагается, что функция , где , а функции удовлетворяют условиям и . Кроме того, полагаем, что выполнены условия согласования
(2.13)
где , , .
Для уравнения (2.2) с условиями (2.6) и (2.7) на отрезке сформулированы шесть задач граничного управления с граничными условиями (2.3) и (2.4), с заданными различными условиями (2.8)(2.10) на функцию колебания и ее производную в фиксированные промежуточные значения.
Номером 1 обозначим задачи, в которых управление реализуется за счет перемещения только одного конца (для определенности, левого) при закрепленном другом конце.
Задача 1A. Граничное управление колебательным процессом с заданными значениями функции колебания в выделенные промежуточные моменты времени (2.8).
Задача 1B. Граничное управление колебательным процессом с заданными значениями производной функции колебания в выделенные промежуточные моменты времени (2.9).
Задача 1С. Граничное управление колебательным процессом с заданными значениями функции колебания и ее производной в выделенные разные промежуточные моменты времени (2.10).
Сформулируем перечисленные задачи управления 1A, 1B, 1C с указанными граничными условиями (2.3).
Требуется найти такое граничное управление , , (см. (2.3)), под влиянием которого колебания системы (2.2) с условиями сопряжения (2.5) из известного состояния (2.6) в начале отрезка переходят в состояние (2.7) в конце отрезка , обеспечивая выполнение следующих значений:
A. функции колебания в фиксированные промежуточные значения времени (см. (2.8));
B. производной функции колебания в фиксированные промежуточные значения времени (см. (2.9));
C. функции колебания и ее производной в фиксированные промежуточные значения и (см. (2.10)).
Номером 2 обозначим далее задачи, в которых управление реализуется за счет перемещения обоих концов системы.
Задача 2A. Граничное управление колебательным процессом с заданными значениями функции колебания в выделенные промежуточные моменты времени (2.8).
Задача 2B. Граничное управление колебательным процессом с заданными значениями производной функции колебания в выделенные промежуточные моменты времени (2.9).
Задача 2C. Граничное управление колебательным процессом с заданными значениями функции колебания и ее производной в выделенные разные промежуточные моменты времени (2.10).
Сформулируем перечисленные задачи управления 2A, 2B, 2C с указанными граничными условиями (2.4).
Требуется найти такие граничные управления и , (см. (2.4)), под влиянием которых колебания системы (2.2) с условиями сопряжения (2.5) из известного состояния (2.6) в начале отрезка переходят в состояние (2.7) в конце отрезка , обеспечивая выполнение следующих значений:
A. функции колебания в фиксированные промежуточные значения времени (см. (2.8));
B. производной функции колебания в фиксированные промежуточные значения времени (см. (2.9));
C. функции колебания и ее производной в фиксированные промежуточные значения и (см. (2.10)).
Для уравнения (2.2) с начальными (2.6) и конечными (2.7) условиями на отрезке времени и функционалами (2.11), (2.12) сформулированы шесть задач оптимального граничного управления с граничными условиями (2.3) и (2.4), с заданными различными условиями (2.8)(2.10) на функцию колебания и ее производную в определенные промежуточные значения из временного интервала. Сохраняя принятую выше нумерацию задач, отметим задачи оптимального граничного управления дополнительно верхним индексом <<0>>.
Задача 1 A. Оптимальное граничное управление колебательным процессом с заданными значениями функции колебания в выделенные промежуточные моменты времени (2.8) и минимизирующие функционал (2.11).
Задача 1 B. Оптимальное граничное управление колебательным процессом с заданными значениями производной функции колебания в выделенные промежуточные моменты времени (2.9) и минимизирующие функционал (2.11).
Задача 1 C. Оптимальное граничное управление колебательным процессом с заданными значениями функции колебания и ее производной в выделенные разные промежуточные моменты времени (2.10) и минимизирующие функционал (2.11).
Сформулируем перечисленные задачи 1 A, 1 B, 1 C оптимального граничного управления с условиями (2.3).
Требуется найти оптимальное граничного управление , (см. (2.3)) под воздействием которого колебания системы (2.2) с условиями сопряжения (2.5) из известного состояния (2.6) в начале отрезка переходят в состояние (2.7) в конце отрезка , обеспечивая минимум функционала (2.11) и выполнение следующих заданных значений:
A. функции колебания в фиксированные промежуточные значения времени (см. (2.8));
B. производной функции колебания в фиксированные промежуточные значения времени (см. (2.9));
C. функции колебания и ее производной в фиксированные промежуточные значения и (см. (2.10)).
Задача 2 A. Оптимальное граничное управление колебательным процессом с заданными значениями функции колебания в выделенные промежуточные моменты времени (2.8) и минимизирующие функционал (2.12).
Задача 2 B. Оптимальное граничное управление колебательным процессом с заданными значениями производной функции колебания в выделенные промежуточные моменты времени (2.9) и минимизирующие функционал (2.12).
Задача 2 C. Оптимальное граничное управление колебательным процессом с заданными значениями функции колебания и ее производной в выделенные разные промежуточные моменты времени (2.10) и минимизирующие функционал (2.12).
Сформулируем перечисленные задачи 2 A, 2 B, 2 C оптимального граничного управления с условиями (2.4).
Требуется найти такие оптимальные граничные управления и , (см. (2.4)), под воздействием которых колебания системы (2.2) с условиями сопряжения (2.5) и граничными условиями из известного состояния (2.6) в начале отрезка переходит в известное состояние (2.7) в конце отрезка , обеспечивая минимум функционала (2.12) и выполнение следующих заданных значений:
A. функции колебания в фиксированные промежуточные значения времени (см. (2.8));
B. производной функции колебания в фиксированные промежуточные значения времени (см. (2.9));
C. функции колебания и ее производной в фиксированные промежуточные значения и (см. (2.10)).
Замечание 2 .Так как во всех задачах управления и оптимального управления в отдельные промежуточные моменты времени ( ) заданы или только значения функции колебания, или только значения производной функции колебания (значения скоростей точек), то использовать подход поэтапного исследования задач нецелесообразно.
В данной работе для всех перечисленных задач по единой схеме предлагается конструктивный подход решения, в котором учитывается специфика промежуточных условий.
Схема построения решений сформулированных задач включает следующие шаги:
Шаг 1. Задачи сводятся к задачам управления с распределенными воздействиями с нулевыми граничными условиями.
Шаг 2. При помощи метода разделения переменных полученные задачи сводятся к задачам управления и оптимального управления для обыкновенных дифференциальных уравнений с заданными начальными, конечными и многоточечными промежуточными условиями.
Шаг 3. При помощи методов теории управления и оптимального управления конечномерными системами с многоточечными промежуточными условиями для произвольного числа первых гармоник строятся граничные управления и оптимальные граничные управления, которые представляются в явном аналитическом виде.
3. Сведение исходных задач к задачам с нулевыми граничными условиями. Для выполнения шага 1 из схемы построения решения перейдем к новой переменной:
(3.1)
что позволит реализовать растяжение или сжатие отрезка относительно точки . При этом с учетом (2.1) вместо отрезка будем иметь отрезок .
Для функции получим на отрезках одинаковой длины одинаковое уравнение
или
(3.2)
с граничными условиями
(3.3)
(3.4)
начальными условиями
(3.5)
промежуточными условиями
(3.6)
(3.7)
(3.8)
конечными условиями
(3.9)
и с условиями сопряжения в точке соединения участков:
(3.10)
Отметим, что для простоты и удобства после замены переменной (3.1) все обозначения функций сохранены.
3.1. Сведение неоднородных граничных условий к нулевым граничным условиям. Поскольку граничные условия (3.3), (3.4) неоднородны, будем строить решение уравнения (3.2) в следующем виде:
(3.11)
где требующая определения функция с однородными граничными условиями
(3.12)
а решение уравнения (3.2) с неоднородными граничными условиями
(3.13)
(3.14)
Функция для условий (3.3) и (3.4) представляется в виде
(3.15)
(3.16)
Подстановка (3.11) в (3.2) и учет (3.15), (3.16) дают следующее уравнение для :
(3.17)
где
(3.18)
(3.19)
Функция удовлетворяет условию (3.10) в точке . Следует отметить, что согласно (3.1) имеем
(3.20)
3.2. Сведение начальных, промежуточных и конечных условий к соответствующим условиям для неоднородного уравнения. Учитывая выражения (3.15), (3.16) для функции и условия согласования (3.20), из известных начальных (3.5), промежуточных (3.6)(3.8) и конечных условий (3.9) получим соответствующие условия для функции . Для задач граничного управления колебаниями смещением левого конца при закрепленном правом конце, т.е. для функции , получим следующие условия: начальные
(3.21)
промежуточные
(3.22)
(3.23)
(3.24)
конечные
(3.25)
Для задач граничного управления колебаниями смещением двух концов для функции получим следующие условия: начальные
(3.26)
промежуточные
(3.27)
(3.28)
(3.29)
конечные
(3.30)
Итак, приходим к задачам управления колебаниями, моделируемыми уравнением (3.17) с однородными граничными условиями (3.12), которые формулируются следующим образом:
Задачи управления с нулевыми граничными условиями.
. Смещение левого конца при закрепленном правом конце. Требуется найти такое граничное управление , , которое переводит колебание системы, описываемое уравнением (3.17), (3.18) с граничными условиями (3.12), из известного начального состояния (3.21) в конечное состояние (3.25), обеспечивая выполнение следующих промежуточных условий: A (3.22); B (3.23); C (3.24).
. Смещение двух концов. Требуется найти такие граничные управления и , , которые переводят колебание системы, описываемое уравнением (3.17), (3.19) с граничными условиями (3.12), из известного начального состояния (3.26) в конечное состояние (3.30), обеспечивая выполнение следующих промежуточных условий: A (3.27); B (3.28); C (3.29).
Задачи оптимального управления с нулевыми граничными условиями.
. Смещение левого конца при закрепленном правом конце. Требуется найти такое оптимальное граничное управление , , которое переводит колебание системы, описываемое уравнением (3.17), (3.18) с граничными условиями (3.12), из известного начального состояния (3.21) в конечное состояние (3.25), обеспечивая выполнение следующих промежуточных условий: A (3.22); B (3.23); C (3.24), и которое минимизирует функционал (2.11).
. Смещение двух концов. Требуется найти такие оптимальные граничные управления и , , которые переводят колебание системы, описываемое уравнением (3.17), (3.19) с граничными условиями (3.12), из известного начального состояния (3.26) в конечное состояние (3.30), обеспечивая выполнение следующих промежуточных условий: A (3.27); B (3.28); C (3.29), и которые минимизируют функционал (2.12).
4. Решение задачи. Применение метода разделения переменных. Перейдем к выполнению шагов 2 и 3. Будем искать решение уравнения (3.17) в виде
(4.1)
Функции , и представим в виде рядов Фурье в базисе , . Подставим далее их значения вместе с функцией в уравнения (3.17)(3.19) и в условия (3.21)(3.30). В результате получим
(4.2)
(4.3)
(4.4)
Для задач под номером 1 (смещение одного конца при закрепленном другом конце) начальные, промежуточные и конечные условия запишутся в виде
(4.5)
(4.6)
(4.7)
(4.8)
(4.9)
Здесь , и коэффициенты Фурье, соответствующие функциям , и .
Для задач под номером 2 (смещение обоих концов) начальные, промежуточные и конечные условия запишутся в виде
(4.10)
(4.11)
(4.12)
(4.13)
(4.14)
Общее решение уравнения (4.2) имеет вид
(4.15)
Учитывая начальные (4.5) (или (4.10)), промежуточные (4.6)(4.8) (или (4.11)(4.13)) и конечные (4.9) (или (4.14)) условия, из (4.15) получим, что функции для каждого должны удовлетворять интегральным соотношениям в виде
(4.16)
(4.17)
(4.18)
(4.19)
где приняты обозначения
(4.20)
(4.21)
(4.22)
(4.23)
Отметим, что задачам управления и оптимального управления с условиями A соответствуют интегральные соотношения (4.16), (4.17), задачам с условиями B соотношения (4.16), (4.18), а задачам с условиями C соотношения (4.16), (4.19). Приведем дальнейшее построение решения (шаг 3 схемы) для задач граничного управления колебаниями, выделяя построение смещением левого конца при закрепленном правом конце и смещением двух концов.
4.1. Построение решения задач граничного управления колебаниями смещением левого конца при закрепленном правом конце. Подставляя выражение функции из (4.3) в соотношения (4.16)(4.19) и интегрируя по частям с учетом условий согласования (2.13), получим из (4.16) следующие соотношения:
(4.24)
а из (4.17), (4.18) и (4.19) получим следующие соотношения:
(4.25)
(4.26)
(4.27)
где
а также
(4.28)
(4.29)
Введем следующие обозначения:
(4.30)
(4.31)
(4.32)
Тогда, с учетом введенных обозначений (4.30)(4.32), соотношения (4.24)(4.27) запишутся следующим образом:
(4.33)
Здесь в верхнем индексе обозначения соответствуют задачам со смещением левого конца при закрепленном правом конце с условиями A, B и C.
На практике, как правило, выбираются несколько первых гармоник колебаний и решается задача синтеза управлений, используя методы теории управления конечномерными системами. Поэтому
(4.34)
с размерностями , при всех .
Для первых гармоник соотношение (4.33), с учетом (4.34), запишется в виде
(4.35)
Из (4.35) вытекает справедливость следующего утверждения.
Теорема 1. Первые гармоник системы (4.2), (4.3) с условиями (4.5)(4.9) вполне управляемы тогда и только тогда, когда для любого вектора (4.34) можно найти управление , , удовлетворяющее условию (4.35).
Для произвольного числа первых гармоник управляющее воздействие , удовлетворяющее интегральному соотношению (4.35), имеет вид (см. [6, 20])
(4.36)
где транспонированная матрица, такая вектор-функция, что
Здесь известная матрица размерностью , при .
Из формулы (4.36) следует, что существует множество управляющих функций, решающих задачи граничных управлений.
Учитывая обозначения (4.29), функции управления представляются в виде
(4.37)
Подставляя из (4.36) (или из (4.37)) управление в (4.3), а найденное для выражение в (4.15), получим функцию , . Далее, из формулы (4.1) будем иметь
где
(4.38)
а функция колебания , , для первых гармоник запишется в виде
(4.39)
Учитывая обозначения (3.1), функция при представляется в виде:
(4.40)
4.2. Построение решения задач граничного управления колебаниями смещением двух концов.Подставим значение функции в виде (4.4) в соотношения (4.16)(4.19). Интегрируя их по частям с учетом условий согласования (2.11)(2.13), из (4.16) получим, что функции и для каждого должны удовлетворять интегральным соотношениям в виде
(4.41)
а из (4.17), (4.18) и (4.19) получим следующие интегральные соотношения:
(4.42)
(4.43)
(4.44)
где
(4.45)
Отметим, что выражения для , , , совпадают с приведенными в формулах (4.20)(4.23), а выражения для , , , приведены в (4.28).
Введем следующие обозначения:
(4.46)
(4.47)
(4.48)
Тогда с учетом введенных обозначений (4.46)(4.48) соотношения (4.41)(4.44) запишутся следующим образом:
(4.49)
Здесь через , , обозначены вектор-функции управления и оптимального управления для задач смещением двух концов, т.е. для задач 2A, 2B, 2C и 2 A, 2 B, 2 C соответственно.
Таким образом, для поиска функции , , для всех перечисленных задач получили бесконечные интегральные соотношения, которые представлены в единой записи (4.49). Введем для первых гармоник следующие обозначения блочных матриц:
(4.50)
размерностей и соответственно. Для первых гармоник с учетом (4.50) соотношение (4.49) запишется в виде
(4.51)
Из (4.51) следует утверждение, аналогичное теореме 1: первые гармоник системы (4.2), (4.4) с условиями (4.10)(4.14) вполне управляемы тогда и только тогда, когда для любого вектора из (4.50) можно найти управление , , удовлетворяющее условию (4.51).
Для произвольного числа первых гармоник управляющее воздействие , удовлетворяющее интегральному соотношению (4.51), имеет вид (см. [6, 20])
(4.52)
где транспонированная матрица, такая вектор-функция, что
(4.53)
Здесь внешнее произведение, известная матрица размерности , для которой предполагается, что , при .
Здесь также из формулы (4.52) следует, что для задач 2A, 2B, 2C существует множество управляющих функций, решающих задачи граничных управлений.
Подставляя из (4.52) величины и в (4.4), а найденное для выражение в (4.15), получим функцию , . Далее, из формулы (4.1) будем иметь
где
(4.54)
а функция колебания , для первых гармоник запишется в виде
(4.55)
Учитывая обозначения (3.1), представим функцию колебания при в следующем виде:
(4.56)
4.3. О дальнейшем построении решения задач оптимального граничного управления колебаниями. В ходе построения решения задач оптимального граничного управления колебаниями, для первых гармоник в случае управления смещением левого конца при закрепленном правом конце получено интегральные соотношения в виде (4.35), а в случае управления смещением двух концов интегральные соотношения (4.51). Ясно, что левая часть соотношения (4.35) или (4.51) линейная операция, порожденная функцией управления на промежутке времени , а функционалы (2.11) или (2.12) являются нормой соответствующего нормированного пространства .
Таким образом, задачу оптимального управления с интегральными условиями (4.35) при функционале (2.11) или с интегральными условиями (4.51) при функционале (2.12) можно рассматривать как проблему моментов, а решение этих задач следует строить с помощью алгоритма решения проблемы моментов (см. [10]).
5. Заключение. Используя методы разделения переменных, теории управления и оптимального управления конечномерными системами с многоточечными промежуточными условиями, предложен конструктивный подход построения граничного управления и оптимального управления неоднородной колебательной системой с заданными значениями функции колебания и производной функции колебания в разные промежуточные моменты времени. Предложенный для одномерного неоднородного волнового уравнения подход можно распространить на другие одномерные и неодномерные колебательные системы.