Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Математические основы теории систем.-2

.pdf
Скачиваний:
19
Добавлен:
05.02.2023
Размер:
2.17 Mб
Скачать

Tp = p1 x1 + p2 x2 +...+ xi xi xi

По теореме Эйлера [4] имеем

 

T

T

2T

= x

x

+ x

 

x

2

x

1

 

 

x1

 

 

x2

 

 

 

 

 

T

 

 

 

 

 

xn .

 

x

+... + p

n

 

 

 

x

 

x

 

i

 

 

i

 

 

T

 

 

+ ...+ xn

x

,

 

 

 

xn

 

что, учитывая соотношение (6.6.16), сведется к уравнению

2T= p1 x1 + p2 x2 +... + pn xn

x

Дифференцируя последнее выражение по xi

2 Tp = p1 x1 + p2 x2 +... + xi xi xi

=2Tp.

,имеем

xpn xin .

Вычитая выражение (6.6.21) из выражения (6.6.23), получим

 

 

 

 

 

 

 

 

 

Tp

 

= −

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xi

 

 

 

 

 

 

Теперь продифференцируем выражение (6.6.19) по pi :

 

Tp

=

Tx

 

x1

+

Tx

 

x2

+...+

 

Tx

 

xn

,

 

pi

x1

pi

x2

pi

 

xi

pi

 

 

 

 

 

 

 

 

 

 

 

 

 

или, с учетом соотношения (6.6.16),

 

 

 

 

 

 

 

 

 

 

 

 

 

Tp

 

= p

x1 + p

x2 +... + p

 

xn .

 

 

 

 

pi

 

 

 

 

 

1 pi

 

 

 

2 pi

 

 

 

 

n pi

 

 

 

(6.6.21)

(6.6.22)

(6.6.23)

(6.6.24)

(6.6.25)

305

Возьмем частную производную по pi

от выражения (6.6.22):

 

2

Tp

= p x1 + p

x2

+... + x

+...+ p

xn .

(6.6.26)

pi

2 pi

 

1 pi

i

 

n pi

 

Вычитая выражение (6.6.25) из равенства (6.6.26), получим

 

 

 

 

Tp

 

 

 

 

 

 

 

= xi .

 

 

(6.6.27)

 

 

 

pi

 

 

Используя формулы (6.6.27) и (6.6.16), уравнение Лагранжа (6.6.14) приведем к виду

Tp

 

 

pi +

 

= ri .

(6.6.28)

xi

Две системы уравнений (6.6.27) и (6.6.28)) называются уравнениями Гамильтона (для кинетической энергии).

В случае консервативной системы входное воздействие определяется выражением

 

 

ri = −

V ,

(6.6.29)

 

 

 

xi

 

 

 

 

 

где V =V (x) – потенциальная энергия, не зависящая от x .

 

Так как лагранжиан определяется формулой (6.6.8)

 

 

 

 

L xi , xi

= Txi , xi V (xi ),

 

 

 

x

 

 

то его частные производные равны

L

 

T

 

L

 

T

V .

=

x

и

=

x

 

 

 

 

xi

xi

xi

xi

 

xi

 

 

 

 

 

 

 

 

306

Сучетом этого, уравнение Эйлера – Лагранжа (6.6.9) можно записать

ввиде

 

 

 

d

 

T

 

T

 

V

 

 

 

 

 

 

 

x

 

=

x

,

 

 

 

 

 

 

 

 

x

 

 

 

 

dt

 

 

x

 

 

 

 

 

 

 

 

xi

 

i

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

d

(grad

T

)= gradx

T

gradx V.

(6.6.30)

 

dt

 

 

 

x

x

 

 

 

 

x

 

 

 

Функция, выражающая полную энергию системы через координаты х и импульсы р, называют функцией Гамильтона Н. То есть

H = Tp +V = H (p,x) .

Дифференцируя выражение (2.6.31), получаем

H = Tp + V xi xi xi

и

H = Tp . pi pi

(6.6.31)

(6.6.32)

(6.6.33)

Подставляя равенства (6.6.29), (6.6.32) и (6.6.33) в уравнения (6.6.28) и (6.6.27), окончательно имеем

dpi

= −

H

,

(6.6.34)

dt

 

 

 

x

 

 

 

 

 

 

i

 

 

dxi

=

H

 

(6.6.35)

dt

 

 

 

p

 

 

 

 

 

 

i

 

 

Две системы уравнений (6.6.34) и (6.6.35) носят название канониче-

ских уравнений Гамильтона.

307

Поскольку при движении консервативной системы ее полная энергия остается неизменной, Н – функция Гамильтона не зависит от времени и

dHdt = 0 . Это действительно так, поскольку

dH

n

 

H dpi +

H

 

,

=

dxi

dt

i=1

 

pi

dt

xi

dt

 

а выражение в скобках, согласно уравнениям (6.6.34) и (6.6.35), равно нулю, если Н явно не зависит от времени.

6.6.3.Уравнение Гамильтона – Якоби

Во многих случаях решения уравнений (6.6.34), (6.6.35) найти не удается. Один из путей решения этих уравнений состоит в переходе от координат (p,x) к другой системе координат (α,β) , относительно которых

преобразованные уравнения имеют более простой вид. Такие преобразования, результатом которых являются новые уравнения все в той же канонической форме, называются каноническими преобразованиями. Отсюда, если

pi = pi (α,β) и xi = xi (α,β)

(6.6.36)

есть канонические преобразования, то уравнения движения в новой системе координат (α,β) будут иметь все тот же канонический вид

dαi = − H , dt ∂βi

(6.6.37)

dβi = H , dt ∂αi

где H – гамильтониан, выраженный в новой системе (α,β) .

Смысл преобразований (6.6.36) и перехода к каноническим уравнениям Гамильтона (6.6.37) состоит в том, чтобы гамильтониан H являлся бы только функцией переменных α и не зависел бы от β . Такую цель позволяет достичь преобразование

308

p =

S

,

β

 

=

S

,

 

i

 

i

xi

 

 

 

∂αi

 

 

 

 

 

 

 

где функция S(α,x) называется производящей функцией. Уравнения движения в этом случае принимают вид

dαi = − H = 0 , dt ∂βi

dβi = H = const . dt ∂αi

(6.6.38)

(6.6.39)

(6.6.40)

Из первой системы уравнений (6.6.39) вытекает, что все αi – констан-

ты. Вторая система уравнений (6.6.40) следует из того факта, что H зависит только от αi , а все αi – константы. Системы уравнений (6.6.39) и

(6.6.40) много проще, чем уравнения (6.6.34) и (6.6.35). Дело за малым: нужно определить производящую функцию S(α,x) , удовлетворяющую

дифференциальному уравнению в частных производных

H

 

S

, x

 

 

 

(α ) .

 

= H

(6.6.41)

 

 

i

 

xi

 

 

i

 

 

 

 

 

 

 

 

Чтобы немного упростить уравнение (6.6.41), проведем следующие

рассуждения.

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как β

 

=

S

, то

 

 

 

 

 

 

 

 

 

 

i

∂αi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dβi

 

 

S

 

2

S

 

 

 

 

 

 

=

 

=

 

.

(6.6.42)

 

 

 

 

dt

 

t

∂αi

t∂αi

 

Найдем производную по t от уравнения (6.6.41), предполагается пока, что H зависит также и от βi

H = H dαi + H dβi . t ∂αi dt ∂βi dt

309

Используя равенства (6.6.37) и (6.6.42), из последнего уравнения имеем

H

 

 

 

 

 

dαi

 

dαi

 

 

 

 

 

 

 

dαi

 

 

 

 

 

=

H

 

 

 

2S

=

H

 

d

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(6.6.43)

t

 

∂αi

 

dt

 

dt

 

t∂αi

 

∂αi

 

dt

 

dt

t

 

Чтобы гамильтониан H не зависел от βi , первое слагаемое в уравне-

нии (6.6.43), согласно соотношению (6.6.39), должно быть равно нулю. Это будет выполняться, если справедливо уравнение

H

= −

d

S

t

 

 

.

 

 

dt

t

Из последнего соотношения получаем уравнение для S

H = −

S

,

(6.6.44)

t

 

 

 

которое называется уравнением Гамильтона – Якоби.

К каноническим уравнениям Гамильтона и Гамильтона – Якоби приходят при синтезе оптимальных систем методом максимума Понтрягина или методом динамического программирования Беллмана [14].

Контрольные вопросы

1.Какой вид имеет стандартная форма (каноническая форма фазовой переменной) записи уравнений состояния?

2.Что такое нормальные координаты системы?

3.Как можно перейти к нормальной (канонической) форме уравнений состояния?

4.Какой вид имеют матрицы в канонической форме уравнений состо-

яния?

5.Что такое переходная (фундаментальная) матрица?

6.Какие методы существуют для вычисления переходной матрицы?

7.Назовите основные свойства переходной матрицы.

8.Что такое сопряженная система и сопряженный оператор?

9.В каком случае переходная матрица нестационарной системы представляет матричную экспоненту?

10.Что такое матрицант и как он вычисляется?

310

11.Каков физический смысл лагранжиана системы?

12.Сформулируйте принцип Даламбера.

13.Запишите простейшее уравнение Эйлера – Лагранжа.

14.В каких случаях возможно аналитическое решение уравнения Эйлера – Лагранжа?

15.Каков физический смысл функции Гамильтона?

16.В чем смысл перехода к каноническим уравнениям Гамильтона – Якоби?

311

ЛИТЕРАТУРА

1.Карпов А.Г. Математические основы теории систем. Часть 1: Учебное пособие. Томск: Изд-во НТЛ, 2007. 184 с. ISBN 978-5-89503-357-9.

2.Карпов А.Г. Математические основы теории систем. Часть 2: Учеб-

ное пособие. Томск: Томский межвузовский центр дистанционного образования, 2002. 138 с.

3.Карпов А.Г. Теория автоматического управления. Часть 2: Учебное пособие. Томск: Изд-во ТМЛ-Пресс, 2012. 264 с. ISBN 978-5-9130- 2136-6.

4.Карпов А.Г. Цифровые системы автоматического управления (Ос-

новы теории): Учебное пособие. Томск: Изд-во НТЛ, 2007. 288 с. ISBN 978-5-89503-358-6.

5. Перегудов Ф.И., Тарасенко Ф.П., Основы системного анализа. – Томск: Изд-во НТЛ, 2003. – 396 с. ISBN 5-89503-004-1.

6. Вунш Г. Теория систем. – М.: Советское радио, 1978. – 288 с.

7. Кузнецов О.П., Адельсон-Вельский Г.М. Дискретная математика для инженера. – М.: Энергоатомиздат, 1988. – 480 с. ISBN 5-283-01563-7.

8. Закревский А.Д. Алгоритмы синтеза дискретных автоматов. – М.: Наука, 1971. – 512 с.

9.Мелихов А.Н. Ориентированные графы и конечные автоматы. – М.: Наука, 1971. – 416 с.

10.Деруссо П., Рой Р., Клоуз Ч. Пространство состояний в теории управления. – М.: “Наука”, 1970. – 620с.

11.Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления. – М.: Госиздат, 1951.

12.Диткин В.А., Прудников А.П. Интегральные преобразования и операционное исчисление. – М.: “Наука”, 1974.

13.Гантмахер Ф.Р. Теория матриц. – М.: “Наука”, 1966.

14.Ту Ю. Современная теория управления. – M.: “Машиностроение”, 1971. – 472c.

312

Учебное издание

Александр Георгиевич Карпов

Математические основы теории систем

Учебное пособие

Издание подготовлено

вавторской редакции

6.7Коорректор Г.И. Иванченко

Верстка макета и дизайн обложки

Редактор

Верстка

Изд. лиц. Подписано к печати.

Формат 60×841/16. Бумага офсетная. Печать офсетная. Гарнитура «Times».

Усл. п. л.. Уч.-изд. л.. Тираж 100 экз.

ООО «Издательство ТМЛ-Пресс» 634050, г. Томск, ул. Советская, 33, оф. 10, тел. (382-2) 52-87-15

Отпечатано Томск. гос. ун-т систем управления и радиоэлектроники, г. Томск, пр. Ленина, 40

125