Матика 2 курс / Дифференциальные уравнения 2б
.pdf111
ектория - замкнутая линия вокруг центра. Когда энергия достигает критического значения Е = 0 частица движется по особой траекто-
рии - сепаратрисе (смысл этого названия будет прояснен ниже).
Наконец, когда энергия Е становится больше нуля, частица способ-
на преодолеть барьер, разделяющий ямы потенциала, и траектория становится замкнутой линией, обходящей все три положения рав-
новесия системы (рис. 4.6). Теперь становится понятным смысл названия критической траектории Е = 0 «сепаратриса». Эта траек-
тория разделяет траектории на два качественно различных вида -
локализованные (около конкретного центра) и делокализованные
(обходящие оба центра и седло).
Особенностью рассматриваемой системы является возмож-
ность точно проинтегрировать уравнение сепаратрисы:
|
|
|
|
dq |
|
|
, |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
(4.19) |
||||
q2 |
q4 |
/ 2 |
|||||||
|
|
|
|||||||
|
|
|
|
и получить его решение в виде
q( ) |
2 |
|
ch . |
(4.20) |
Сепаратриса начинается вблизи седла при τ = - ∞ и вновь при-
ближается к седлу при τ = + ∞, проходя точку максимального уда-
ления от седла при τ = 0 . Особенность сепаратрисы заключается в следующем. Если для обхода любой из замкнутых траекторий (рис.
4.6) требуется конечное время, то для обхода по сепаратрисе - бес-
конечное время (в силу теоремы единственности, фазовые траекто-
рии не пересекаются).
112
q
Рис. 4.6. Фазовые траектории в отсутствии затухания и вынуж-
дающей силы
ЗАДАЧА 4.8. Проинтегрируйте (4.19), используя замену u = 1/q.
Рассмотрим изменение фазовых траекторий при включении слабого затухания 0 < γ << 1 . В этом случае система уравнений имеет вид
dqdt p ,
dp |
|
(4.21) |
|
|
p q q3. |
||
dt |
|||
|
|
Его положения равновесия расположены там же, где и у системы без затухания, однако их характер изменяется: р = q = 0 по-
прежнему седло; р = 0, q = ± 1 - устойчивые фокусы.
ЗАДАЧА 4.9. Исследуйте тип положений равновесия (4.21).
В присутствии затухания энергия убывает, поскольку в силу
113
уравнения (4.21)
dHdt p dpdt (q3 q) dqdt p2 0.
Отсутствие сохранения энергии приводит к разрушению за-
мкнутости сепаратрисы и разбиению фазового пространства вдали от фокусов на систему двух чередующихся качественно различных слоев: система «притягивается» к определенному фокусу только в том случае, если она находится во вполне определенном слое фазо-
вого пространства (рис. 4.7).
q
Рис. 4.7. Разбиение фазового пространства на слои в присутствии затухания незамкнутыми сепаратрисами
Исследуем, как изменяется характер фазовых траекторий в от-
сутствии затухания при действии на систему слабой 0 < f << 1 вы-
нуждающей силы частоты ω. Уравнения движения имеют вид
114
dqdt p ,
dp |
q3 q f cos . |
|
|
|
(4.22) |
||
dt |
|||
|
Будем искать поправки, обусловленные внешней силой, к не-
возмущенной траектории (4.18) (q0 (τ), p0 (τ)) . Представляя исход-
ное решение в виде
q( ) q0 ( ) q; |
p( ) p0 ( ) p |
и, выполняя линеаризацию в (4.22), получаем
|
|
|
|
d ( q) |
|
d (p) |
, |
|
|
|
|
|
dt |
|
dt |
||
|
|
|
|
|
|
|
||
|
d (p) |
|
( 3q 2 |
( ) 1) q f cos t |
||||
|
|
|||||||
|
dt |
0 |
|
|
|
|
||
|
|
|
|
|
|
|
||
или |
|
|
|
|
|
|
d 2 ( q) |
|
2 ( ) q f cos , |
|||
dt 2 |
|||||
|
|
(4.23) |
|||
2 |
|
2 |
|||
( ) 3q |
( ) 1. |
||||
|
|
|
0 |
|
Это уравнение описывает вынужденные колебания линейной системы с переменной частотой. В отсутствии резонанса исходная фазовая траектория получает, согласно (4.23), приращение, про-
порциональное амплитуде вынуждающей силы. Если исходная тра-
ектория замкнутая, за один период невозмущенных колебаний воз-
мущенная траектория пересекает возмущенную конечное число раз.
Если невозмущенная траектория - сепаратриса, то соответ-
ствующая ей возмущенная траектория пересечет сепаратрису бес-
конечное число раз. Это связано с тем фактом, что обход по сепара-
115
трисе требует бесконечного времени. Следовательно, вынуждаю-
щая сила превращает сепаратрису исходной системы из гладкой кривой в «дрожащую», с амплитудой пропорциональности f (рис.
4.8).
Рассмотрим качественную картину возникновения хаотиче-
ского режима при совместном действии вынуждающей силы и за-
тухания. Под воздействием вынуждающей силы сепаратрисы си-
стемы с затуханием (рис. 4.7), разделяющие фазовое про-
странство на слои, ведущие к различным фокусам, начинают
«дрожать».
~ f
q
Рис. 4.8. Сепаратриса возмущенная вынуждающей силой f cos
Увеличение затухания γ «раздвигает» сепаратрисы, а возрас-
тание амплитуды силы f усиливает их «дрожание». Если удается установить, что сепаратрисы имеют точку пересечения, это будет означать, что они пересекаются бесконечное число раз. Следова-
тельно, необходимо иметь некоторый критерий пересечения сепа-
116
ратрис.
Отметим, что в настоящее время нет законченных критериев возникновения хаоса, дающих стопроцентный результат. Все суще-
ствующие критерии помогают лишь установить некоторые грани-
цы, приближенно указывающие, когда возникновение хаоса воз-
можно.
Мы опишем ниже критерий Мельникова, связывающий воз-
можность образования хаотического режима с изменением инте-
грала энергии вдоль невозмущенной сепаратрисы.
Согласно уравнению (4.15) скорость изменения энергии (4.17)
на сепаратрисе имеет вид
dH |
P |
2 ( ) f P ( ) cos( ( |
|
)) |
, |
(4.24) |
|
0 |
|||||
dt |
S |
S |
|
|||
|
|
|
|
|
|
где введена начальная фаза вынуждающей силы 0 , фиксирую-
щая начало ее воздействия на невозмущенное движение. Вычисляя
ps dqs / dt по уравнению сепаратрисы (4.20) |
и интегрируя по , |
|||||||||||||
получаем изменение энергии вдоль сепаратрисы: |
|
|||||||||||||
|
sh2 |
|
|
|
|
|
|
sh cos( ( 0 |
)) |
|
||||
H S ( 0 ) 2 |
d f 2 |
|
d |
|||||||||||
|
|
|
|
|
|
|
|
|||||||
ch2 |
|
|
ch2 |
|
||||||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
4 |
|
|
|
|
sin( 0 ) |
|
(4.25) |
||||
|
|
|
|
|
|
|
|
|
||||||
|
|
f |
2 |
|
|
|||||||||
|
|
|
|
. |
|
|
||||||||
|
3 |
ch( / 2) |
|
|
ЗАДАЧА 4.10. Получите формулу (4.25), воспользовавшись подходящими таблицами интегралов.
Смысл формулы (4.25) понятен. Первое слагаемое – 4 γ / 3
117
описывает убывание энергии за счет затухания, второе - ее измене-
ние за счет действия вынуждающей силы. Знак второго слагаемого чувствителен к начальной фазе ωτ0. Знак H ( 0 ) отрицательный,
если
|
|
|
|
|
|
4 |
|
f 2 |
|||
|
|
|
|
||
3 |
сh ( / 2) |
||||
|
и является неопределенным в противоположном случае. Именно условие неопределенности знака H ( 0 ) и служит критерием хао-
тического режима:
|
|
|
|
|
|
ch / 2 |
. |
|
f |
|
|
2 |
|
||||
|
|
|
|
|
(4.26) |
|||
|
|
3 / 2 |
||||||
|
|
Как показывает сравнение критерия (4.26) с результатами фи-
зических и численных экспериментов, он дает нижнюю границу области хаотических режимов. Таким образом, критерий хаоса
(4.26) следует понимать в следующей трактовке.
Если параметры системы γ, f, ω не удовлетворяют критерию
(4.26), движение системы строго предсказуемо.
Если γ, f, ω удовлетворяют (4.26), возможно возникновение хаотического режима, при котором движение системы становится непредсказуемым.
118
ПРИЛОЖЕНИЕ 1.
НЕКОТОРЫЕ СВЕДЕНИЯ ИЗ ЛИНЕЙНОЙ АЛГЕБРЫ
Наиболее важными примерами линейных пространств могут
служить множество векторов размерности n : Rn ; множество мат-
риц порядка n : an ; множество функций, непрерывных на от-
резке. Для любых элементов этих множеств определены линейные операции - умножение на число и сложение.
Базисом в линейном пространстве Rn называют упорядочен-
ную систему линейно независимых векторов (элементов простран-
|
|
|
ства): E1 e1, e2 |
, , en . Число векторов базиса совпадает с размер- |
ностью пространства. Любой элемент пространства (вектор) может быть выражен через элементы базиса единственным образом:
|
|
|
|
|
|
X |
x1e1 |
x2e2 |
... xnen xiei , |
||
где числа x1, x2 ,..., xn |
называют координатами вектора в данном |
базисе. При этом вектор задается своими координатами, например,
в виде матрицы-столбца:
x1x2
X... .
....
xn
Впространстве Rn определено бесконечное количество бази-
119
|
|
|
|
|
также базис в пространстве Rn и |
|||||
сов. Пусть E2 e |
, e2 |
, , en |
||||||||
|
1 |
|
|
|
|
|
|
|
|
|
вектор X имеет в этом базисе координатное представление |
||||||||||
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
|
... |
. |
|
|
||
|
|
|
|
|
|
.... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
Тогда координаты вектора в базисах |
E1 |
и E2 связаны мат- |
|||||||
ричными уравнениями |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
TX ; |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
T 1X , |
|
|
||||
где |
T T 1 T 1T E . |
|
|
|
|
|
|
|
||
|
Здесь T матрица перехода от базиса E1 |
к базису E2 . |
||||||||
|
|
|
t11 |
|
t12 |
t1n |
|
|||
|
|
|
|
|
|
|
|
t2n |
|
|
|
|
|
t21 |
t22 |
|
|||||
|
|
|
T |
|
|
|
|
. |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
tn2 |
|
|
|
|
|
|
|
tn1 |
tnn |
||||||
|
Столбцы этой матрицы задают координаты базисных векторов |
|||||||||
E2 |
в базисе E1 . |
|
|
|
|
|
|
|
|
|
|
При известной матрице Т представления одной и той же мат- |
|||||||||
рицы А в различных базисах связаны преобразованием подобия |
||||||||||
|
|
|
A |
T 1 |
A T; |
|
|
|||
|
|
|
E2 |
|
|
|
E1 |
|
|
120
A T A T |
1 |
. |
(П 1.1) |
|
E1 |
E2 |
|
||
|
|
|
Матрицы, связанные преобразованием подобия, называют эк-
вивалентными друг другу:
A A .
E1 E2
Можно выбрать базис, в котором матрица А имеет наиболее про-
стой вид. Таким базисом является базис из собственных векторов.
|
Rn называют собственным вектором |
||
Ненулевой вектор X в |
|||
матрицы А, если найдется число такое, что |
|
||
|
|
|
|
|
AX |
X . |
(П.1.2) |
Число называют собственным значением матрицы А. Раскрывая выражение (П.1.2) получаем систему однородных уравнений
a11 x1 a12 x2 a1n xn 0 a21x1 a22 x2 a2n xn 0
................................................... (П.1.3) an1x1 an2 x2 ann xn 0,
которая имеет ненулевое решение тогда и только тогда, когда опре-
делитель ее равен нулю:
|
|
|
|
|
a11 |
a12 |
... |
a1n |
|
|
|
|
|
|
|
|
|
||||
Det |
|
A E |
|
|
a21 |
a22 |
... |
a2n |
0 . |
(П.1.4) |
|
|
|||||||||
|
|
......... ......... ......... |
..... |
|||||||
|
|
|
|
|||||||
|
|
|
|
|
an1 |
an2 |
... |
ann |
|
|
Уравнение (П.1.4) называют характеристическим уравнением матрицы А. Собственные значения матрицы А являются корнями