Решение задачи

Система уравнений, описывающая движение фазовыми координатами

, (3)

где

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

- вектор фазовых координат.

Записываем функцию Гамильтона для нашей задачи:

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

Для вспомогательных переменных построим гамильтонову систему:

,

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

Найдем общее решение системы для .

Мы получили однородное дифференциальное уравнение 2-ого порядка. Составим характеристическое уравнение:

Так как корни характеристического уравнения действительные и не равны имеем:

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

Отсюда следует, что:

при максимум функции Гамильтона достигается при u=1.5,

при максимум функции Гамильтона достигается при u=-1, т.е.

Из максимума функции Н видно, что каждое оптимальное управление u(t), где является кусочно-постоянной функцией принимающие значение 1.5 или -1,и имеющей, не более двух интервалов знако-постоянства, так как поведение функции определяет экспонента, которая не более одного раза пересекает t.

В конечном итоге получим:

и поведение управления (рис.1).

Рис.1

1.Для отрезка времени на котором мы имеем систему ДУ:

Найдем из 1-ого уравнения нашей системы:

найдем из 2-ого уравнения нашей системы:

,

получаем: ,

отсюда находим:

,

проинтегрировав, получим общее решение:

обозначим: , получим:

, где S принимает любое значение.

1. Пусть S>0, тогда имеем:

- это семейство гипербол с центром в точке (3.33;0).

a= – действительная полуось.

b= – мнимая полуось.

Уравнение обеих асимптот имеет вид:

, отсюда

Точка пересечения асимптот O1(0,3.33).

Это семейство гипербол показано на рис.2. Для определения движения по гиперболе точки с координатами х и х2 берем уравнение:

А) Если наша правая часть , имеем: , если , то или возрастает, т.е. движение вверх.

Б) Если наша производная меньше 0, имеем:

, если , то или убывает, т.е. движение вниз.

2. Пусть S<0. Тогда уравнение преобразуется в уравнение:

- это семейство гипербол с центром в точке (3.33;0).

a= – действительная полуось.

b= – мнимая полуось.

Уравнение асимптот такое же, что и в случае S>0, и имеет вид:

; Т.к. гиперболы сопряженные (сопряженные гиперболы – гиперболы, уравнения которых имеют вид: , а нашем случае: и ), а сопряженные гиперболы имеют одинаковые уравнения асимптот.

Точка пересечения асимптот O1(0,3.33).

Это семейство гипербол показано на рис.2. Для определения движения по гиперболе точки с координатами х и х2 берем уравнение:

А) Если наша правая часть , имеем: , если , то или возрастает, т.е. движение вверх.

Б) Если наша производная меньше 0, имеем:

, если , то или убывает, т.е. движение вниз.

2. Рассмотрим интервал, на котором .Тогда из (3) имеем:

,

найдем из 1-ого уравнения нашей системы:

найдем из 2-ого уравнения нашей системы:

,

получаем:

,

отсюда находим:

,

проинтегрировав, получим:

обозначим ,

получим:

1. Пусть S>0, тогда

- это семейство гипербол с центром в точке (-2.22;0).

a= – действительная полуось,

b= – мнимая полуось.

Уравнение обеих асимптот имеет вид:

, отсюда

Точка пересечения асимптот O2(0,-2.22).

Это семейство гипербол показано на рис.3. Для определения движения по гиперболе точки с координатами х и х2 берем уравнение:

А) Если наша правая часть больше 0, имеем:

, если , то или возрастает, т.е. движение вверх.

Б) Если наша левая часть меньше 0, имеем: , если , то или убывает, т.е. движение вниз.

2. Пусть S<0. Тогда уравнение преобразуется в уравнение:

- это семейство гипербол с центром в точке (-2.22;0).

a= – действительная полуось,

b= – мнимая полуось.

Уравнение асимптот такое же, что и в случае S>0, и имеет вид:

Точка пересечения асимптот O2(0,-2.22).

Это семейство гипербол показано на рис.3. Для определения движения по гиперболе точки с координатами х и х2 берем уравнение:

А) Если наша правая часть больше 0, имеем:

, если , то или возрастает, т.е. движение вверх.

Б) Если наша левая часть меньше 0, имеем: , если , то или убывает, т.е. движение вниз.

В начало координат приходит две траектории без переключения. Каждое оптимальное управление u(t) является кусочно-постоянной функцией принимающей -1 и 1.5. Если u(t) сначала в течении некоторого времени было равно +1.5, а затем -1, то фазовая траектория состоит из двух кусков гиперболы примыкающих друг к другу, причем второй кусок гиперболы лежит на гиперболе, которая проходит через (0;0).

Уравнение траектории, приходящее в начало координат c u=-1 без переключения имеет вид:

, x1=0, x2=0, отсюда имеем:

, S=1.11, следовательно:

ОА: - уравнение траектории приходящее в начало координат при u=-1 без переключения.

Уравнение траектории, приходящее в начало координат с u=1.5 без переключения имеет вид:

, x1=0, x2=0, отсюда имеем:

, S=2.4975, следовательно:

ОВ: - уравнение траектории приходящее в начало координат при u=1.5 без переключения.

В начало координат мы можем придти по траекториям без переключения и с одним переключением. Без переключения мы приходим по ОА и ОВ, с одним переключением мы должны по траектории попасть на ОА или ОВ(рис.4)

Рис.4

Рис.2

Рис.3


Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  



double arrow
Сейчас читают про: