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

Матика 2 курс / Дифференциальные уравнения 2б

.pdf
Скачиваний:
22
Добавлен:
10.05.2015
Размер:
2.24 Mб
Скачать

111

ектория - замкнутая линия вокруг центра. Когда энергия достигает критического значения Е = 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) называют характеристическим уравнением матрицы А. Собственные значения матрицы А являются корнями