Хаотическая динамика и классическая теория устойчивости движения
.pdfЧАСТОТНЫЕ ОЦЕНКИ ПЕРИОДА КОЛЕБАНИЙ СИСТЕМ |
151 |
Введем передаточную функцию системы (11.1)
W (p) = c (A − pI)−1b,
где p — комплекная переменная.
Теорема 11.1. Пусть для натурального числа N выполнено неравенство
1 |
+ ReW |
exp |
2π |
ij 0 |
(11.3) |
|
|||||
k |
|
|
N |
|
|
при всех j = 0, . . . , N −1. Тогда не существует N-периодической последовательности x(t), удовлетворяющей уравнению (11.1) и включению
c x(t) |
|
Ω, |
|
t |
|
Z. |
(11.4) |
|
|
|
|
|
ДОКАЗАТЕЛЬСТВО ТЕОРЕМЫ.
Рассмотрим последовательность u(t) t = 0, . . . , N − 1 и ее дискретное преобразование Фурье
|
|
|
|
|
N −1 |
|
|
|
2π |
|
|
|
|
v(j) = u(t) exp |
|
ijt . |
(11.5) |
||||||
|
|
|
|
||||||||
|
|
|
|
|
t=0 |
|
|
|
N |
|
|
Напомним, что |
|
|
|
|
|
|
|
|
|
|
|
N − |
1 |
|
2π |
|
|
N |
при k = 0 |
|
|||
exp |
ikt |
= |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|||||
t=0 |
|
N |
0 при k = 1, . . . , N − 1 . |
|
|||||||
Поэтому для |
преобразований |
Фурье |
последовательностей |
u1(t) |
|||||||
и u2(t) справедливо равенство |
|
|
|
|
|
|
|||||
|
|
N −1 |
|
|
|
|
N −1 |
|
|||
|
|
v1(j) v2(j) = N u1(t) u2(t). |
(11.6) |
||||||||
|
|
j=0 |
|
|
|
|
t=0 |
|
|||
152 ГЛАВА 11
В том случае, если u(t + N) = u(t), t Z, имеем формулу
N −1 |
|
2π |
|
|
2π |
|
|
u(t + 1) exp |
|
ijt |
= exp − |
ij v(j). |
(11.7) |
||
|
|
|
|||||
t=0 |
|
N |
|
N |
|
||
Обозначим через y(j) и w(j) преобразования Фурье-последова- тельностей x(t) и ϕ(t, c x(t)), где t = 0, 1, . . . , N − 1. Предполагая, что x(t + N) = x(t), из (11.7) получим равенство
c y(j) = |
|
|
|
− 2π ij w(j). |
(11.8) |
|
− |
W exp |
|
||||
|
|
|
N |
|
|
|
|
|
|
|
|
||
Используя (11.6) и (11.8), получим следующую цепочку равенств:
N −1
t=0
ϕ(t, c x(t))(c x(t) |
|
ϕ(t, c x(t) |
) = |
|
|
|
|
||||||||||
− |
|
k |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
N −1 |
|
|
|
|
|
2π |
1 |
|
|
||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|||||
= |
|
j=0 w(j) |
−W |
|
exp − |
|
ij |
|
− k |
w(j) = |
|||||||
N |
|
N |
|||||||||||||||
|
1 |
N −1 |
|
|
2π |
1 |
|
||||||||||
|
|
|
|
|
|
|
|
2 |
|||||||||
|
|
= − |
|
|
|
ReW |
exp |
|
ij |
+ k |
|w(j)| . |
||||||
|
|
N |
j=0 |
N |
|||||||||||||
Из условия (11.2) следует, что первая сумма в этой цепочке равенств положительна, а из условия (11.3) вытекает, что последнее выражение в этой цепочке неположительно. Полученное противоречие доказы-
вает утверждение теоремы. |
|
ПРИМЕР 11.1. Рассмотрим уравнение |
|
θ(t + 1) = θ(t) − α sin θ(t), |
(11.9) |
описывающее динамику дискретной системы фазовой автоподстройки частоты [40, 41, 98]. Здесь α — положительный параметр.
Представим уравнение (11.9) в виде (11.1) с c = 1, b = −1,
A = 1 − ε,
ϕ(t, σ) = α sin σ − εσ.
ЧАСТОТНЫЕ ОЦЕНКИ ПЕРИОДА КОЛЕБАНИЙ СИСТЕМ |
153 |
В этом случае
1
W (p) = p − 1 + ε .
Хорошо известно [40], что при α < α1, где α1 — корень уравнения
*α2 − 1 = π + arccos α1 ,
отрезок [−π, π] отображается в себя: если θ(t) [−π, π], то θ(t + 1)[−π, π]. Кроме того, в этом случае для θ(0) (−π, π) имеют место соотношения
|
|
|
|
|
|
|
" |
|
|
|
|
|
1 |
" |
|
|
|
|
θ(t) |
| |
|
|
= |
"* |
α2 |
|
1 |
|
arccos |
" |
|
(11.10) |
|
lim |
θ |
0 |
− |
− |
, |
|||||||||||
t + |
∞ |
| |
|
|
" |
|
|
|
α |
" |
|
|
||||
→ |
|
|
|
|
|
" |
|
|
|
|
|
|
" |
|
|
|
при α 2 и
lim |θ(t)| = 0
t→+∞
при α < 2.
В дальнейшем будем рассматривать случай α > 2. Определим множество Ω следующим образом:
5
Ω = [−θ0, 0) (0, θ0].
Из оценки (11.10) следует, что для ненулевого N-периодического решения уравнения (11.9) с начальными данными θ0 (−π, π) справедливо включение θ(t) Ω, t Z, t 0.
При k = α − ε справедливы неравенства (11.2), если
|
|
|
|
ε < |
α sin θ0 |
. |
|
(11.11) |
|
|
|
|
|
|
|
||||
|
|
|
|
|
θ0 |
|
|
||
Рассмотрим здесь случай N = 3. Неравенство (11.3) примет вид |
|||||||||
|
1 |
|
+ |
cos ω − 1 + ε |
0 |
(11.12) |
|||
α − ε |
|||||||||
|
|
(cos ω − 1 + ε)2 + (sin ω)2 |
|
|
|||||
при ω = 0, ω = 2π/3, ω = 4π/3.
154 |
|
|
|
ГЛАВА 11 |
|
Ясно, что это неравенство выполнено, если |
|||||
1 |
|
|
ε − 23 |
||
|
|
|
+ |
|
0. |
|
α − ε |
(ε − 23 )2 + 43 |
|||
Последнее неравенство можно переписать следующим образом:
|
23 α − 3 |
(11.13) |
|
ε > |
|
. |
|
|
|||
|
α − 23 |
|
|
При некотором ε условия (11.11) и (11.13) выполнены, если имеет место неравенство
α sin θ0 |
|
23 α − 3 |
(11.14) |
|
|
> |
|
. |
|
|
|
|||
θ0 |
|
α − 23 |
|
|
Неравенство (11.14) выполнено при α (2, 3.7).
Таким образом, при α (2, 3.7) уравнение (11.9) не может иметь цикла периода 3. Заметим, что на этом интервале для уравнения (11.9) наблюдается переход к хаосу через бифуркации удвоения
периода [40]. |
|
ПРИМЕР 11.2. Рассмотрим логистическое уравнение [70, 112]. |
|
θ(t + 1) = μθ(t)(1 − θ(t)), θ(0) [0, 1], |
(11.15) |
где μ — положительный параметр.
Будем рассматривать здесь случай μ (1, 4). (При μ < 1 все траектории уравнения (11.15) θ(t) стремятся при t → +∞ к устойчивому состоянию равновесия θ(t) ≡ 0). В этом случае уравнение (11.15) имеет два состояния равновесия θ(t) ≡ 0, θ(t) ≡ θ0 = 1 − μ−1.
Очевидно, что при θ(0) (0, 1) выполнено включение
θ(1) (0, μ/4].
ЧАСТОТНЫЕ ОЦЕНКИ ПЕРИОДА КОЛЕБАНИЙ СИСТЕМ |
155 |
||||||||||||||||
Отсюда и из неравенства |
|
|
|
|
|
|
|
|
|
||||||||
|
|
μθ(1 − θ) > θ, θ (0, θ0 ) |
|
||||||||||||||
вытекает оценка |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
μ2 |
|
||||||
|
|
|
|
|
lim |
θ(t) > |
|
|
|
(4 − μ). |
(11.16) |
||||||
|
|
|
|
|
16 |
||||||||||||
|
|
|
|
|
t→+∞ |
|
|
|
|
|
|
|
|
|
|||
В дальнейшем будем предполагать, что выполнено неравенство |
|||||||||||||||||
|
|
|
|
|
|
μ2 |
1 |
|
|
||||||||
|
|
|
|
|
|
|
(4 − μ) > μ . |
(11.17) |
|||||||||
|
|
|
|
|
|
16 |
|||||||||||
В этом случае из (11.16) и (11.17) имеем оценку |
|
||||||||||||||||
|
|
|
μθ(t)(1 − θ(t)) − θ0 |
< −a1 |
(11.18) |
||||||||||||
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
θ(t) − θ0 |
|
|||||||||||
для любой нетривиальной N — периодической последовательно- |
|||||||||||||||||
сти θ(t). Здесь |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
2 |
(4 − μ)) − (1 − μ1 ) |
|
|
||||||
|
|
μ |
(4 − μ)(1 − |
μ |
. |
|
|||||||||||
a1 = |
|
16 |
16 |
|
|||||||||||||
|
|
|
|
2 |
|
|
|||||||||||
|
|
|
|
|
|
1 − μ1 |
|
|
|
|
|||||||
|
|
|
|
|
|
− |
μ |
(4 − μ) |
|
||||||||
|
|
|
|
|
|
16 |
|
||||||||||
Поскольку для такой последовательности выполнено неравенство θ(t) μ/4, получим оценку
|
μθ(t)(1 − θ(t)) − θ0 |
> −a2, |
(11.19) |
|||||||
|
|
|
|
|
|
|||||
|
θ(t) − θ0 |
|
|
|
|
|||||
где |
|
|
|
|
|
|
|
|
||
|
|
1 − μ1 |
2 |
|
|
|
|
|||
|
|
− |
μ |
(4 − μ) |
|
|||||
|
a2 = |
16 |
|
|||||||
|
|
|
|
|
|
|
. |
|
||
|
μ1 |
|
|
μ |
|
|
|
|||
|
|
+ |
− 1 |
|
||||||
|
|
4 |
|
|||||||
156 ГЛАВА 11
Применяя теорему с A = −a1, b = −1, c = 1, σ = θ − θ0, ϕ(t, σ) = = −(μθ(1 − θ) − θ0) + a1(θ − θ0) и
1
W (p) = p + a1 ,
получим, что для любого нетривиального N — периодического решения выполнено условие (11.2) с k = a2 − a1. Условие (11.3) примет вид
1 |
|
a1 + cos |
2π j |
|
|
|
|
|
|
||||
+ |
|
2π |
|
N |
|
0, |
|
j = 0, . . . , |
N − |
1. |
|||
k |
|
(a1 + cos |
j)2 + (sin |
2π |
j)2 |
|
|
|
|||||
|
|
N |
N |
|
|
|
|
||||||
Таким образом, для N = 3 получим следующее условие отсутствия циклов периода 3:
1 |
|
a1 − 21 |
|
|
+ |
|
0. |
a2 − a1 |
(a1 − 21 )2 + 43 |
||
Последнее неравенство можно записать следующим образом:
2(1 + a1a2) a1 + a2. |
|
Это неравенство выполнено при |
|
μ 3,635. |
(11.20) |
Таким образом, если выполнено условие (11.20), то уравнение (11.15) не имеет циклов периода 3.
Хорошо известно [70,112], что при бифуркационном параметре μ = 3 состояние равновесия θ(t) ≡ θ0 теряет устойчивость и рождается цикл периода 2. При μ = 3,44948 . . . цикл периода 2 теряет устойчивость и рождается цикл периода 4. Далее на интервале (3,52, 3,57) происходят аналогичные бифуркации удвоения пери-
ода. В результате численного анализа уравнения (11.15) при μ = 1 + |
|
+ √8 = 3,8284 . . . обнаружен цикл периода 3 [70, 112]. |
|
Литература
1. Айзерман М. А., Гантмахер Ф. Р. Абсолютная устойчивость регулируемых систем. — М.: Изд-во АН СССР, 1963. — 140 c.
2. Анищенко В. С., Астахов В. В., Вадивасова Т. Е., Нейман А. Б., Стрелкова Г. И., Шиманский-Гайер Л. Нелинейные эффекты в хаотических и стохастических системах. Институт компьютерных исследований. Серия «Компьютинг в математике, физике, биологии». — МоскваИжевск, 2003. — 529 с.
3. Анищенко В. С. Сложные колебания в простых системах. — М.: Наука, 1990. — 311 c.
4. Бабин Ф. В., Вишик М. И. Аттракторы эволюционных уравнений с частными производными и оценка их размерности // Успехи математических наук. — 1983. — т. 38, вып. 4. — с. 133–187.
5. Бабин А. В., Вишик М. И. Аттракторы эволюционных уравнений. — М.: Наука, 1989. — 293 c.
6. Барабанов Н. Е, Гелиг А. Х., Леонов Г. А., Лихтарников А. Л., Матвеев А. С., Смирнова В. Б., Фрадков А. Л. Частотная теорема (лемма Якубовича – Калмана) в теории управления // Автоматика и телемеханика. — 1996. — № 10. — c. 3–40.
7. Барбашин Е. А., Табуева В. А. Динамические системы с цилиндрическим фазовым пространством. — М.: Наука, 1969. — 299 с.
8. Бойченко В. А., Леонов Г. А. Об оценках размерности аттракторов отображения Хенона // Вестник Санкт-Петербургского университета. — 2000. — сер. 1, вып. 3. — с. 8–13.
9. Былов Б. Ф., Виноград Р. Э., Гробман Д. М., Немыцкий В. В. Теория показателей Ляпунова и ее приложение к вопросам устойчивости. — М.: Наука, 1966. — 576 с.
10. Вышнеградский И. А. О регуляторах прямого действия // Изв. СанктПетербургского технологического института. — 1877. — c. 21–62.
158 |
ЛИТЕРАТУРА |
11. Гелиг А. Х., Леонов Г. А., Якубович В. А. Устойчивость нелинейных систем с неединственным состоянием равновесия. — M.: Наука, 1978. — 400 c.
12. Демидович Б. П. Лекции по математической теории устойчивости. — М.: Наука, 1967. — 472 с.
13. Демидович Б. П. Орбитальная устойчивость решений автономных систем. I, II // Дифференциальные уравнения. — 1968, № 4. — с. 575–588, № 8. — с. 1359–1373.
14. Дмитриев А. С., Кислов В. Я. Стохастические колебания в радиофизике и электронике. — М.: Наука, 1989. — 277 с.
15. Жуковский Н. Е. О прочности движения. Собрание сочинений. Т. 1. — М.: Гостехиздат, 1948. — c. 67–160.
16. Жуковский Н. Е. О прочности движения // Ученые Записки Московского Университета. — 1882. — № 4.
17. Заславский Г. М. Стохастичность динамических систем. — М.: Мир, 1984. — 271 c.
18. Зорич В. А. Математический анализ. Ч. I, II. — М.: Наука, 1981. — 543.; 1984. — 640 с.
19. Ильяшенко Ю. С., Вейгу Ли. Нелокальные бифуркации. — М.: МЦНМО-Ро, 1999. — 415 c.
20. Кузнецов С. П. Динамический хаос. — М.: Физматлит, 2001. — 296 с.
21. Ладыженская О. А. О нахождении минимальных глобальных аттракторов для уравнений Навье – Стокса и других уравнений с частными производными // Успехи математических наук. — 1987. — т. 42, вып. 6. — с. 25–60.
22. Леонов Г. А. Об одной модификации контрпримера Перрона // Дифференциальные уравнения. — 2003. — т. 39, № 11. — с. 1566–1567.
23. Леонов Г. А. Критерии неустойчивости по первому приближению для нестационарных линеаризаций // Прикладная математика и механика. — 2004. — № 6.
24. Леонов Г. А. Проблема обоснования первого приближения в теории устойчивости движения // Успехи механики. — 2003. — № 3.
25. Леонов Г. А. Формулы ляпуновской размерности аттракторов Хенона и Лоренца // Алгебра и анализ. — 2001. — т. 13, вып. 3. — с. 155–170.
|
ЛИТЕРАТУРА |
159 |
26. |
Леонов Г. А. Странные аттракторы и классическая теория устойчиво- |
|
|
сти движения // Успехи механики. — 2002. — т. 1, № 3. — с. 3–42. |
|
27. |
Леонов Г. А. Математические проблемы теории управления. Мотива- |
|
|
ция к анализу. — СПб.: Изд-во С.-Петербургского ун-та, 1999. — 160 с. |
|
28. |
Леонов Г. А., Пономаренко Д. В. Критерии орбитальной устойчиво- |
|
|
сти траекторий динамических систем // Известия высших учебных |
|
|
заведений. Математика. — 1993. — № 4. — c. 88–94. |
|
29. Леонов Г. А., Полтинникова М. С. О ляпуновской размерности аттрактора диссипативного отображения Чирикова // Труды Санкт-Петер- бургского математического общества. — 2002. — т. 10. — с. 186–198.
30. Леонов Г. А. Об устойчивости по первому приближению // Прикладная математика и механика. — 1998. — т. 69, вып. 4. — с. 548–555.
31. |
Леонов Г. А. Неустойчивость по первому приближению для неста- |
|
ционарных линеаризаций // Прикладная математика и механика. — |
|
2002. — т. 66, вып. 2. — с. 330–333. |
32. |
Леонов Г. А. Об оценках хаусдорфовой размерности аттракторов. |
|
Вестник Ленинградского университета. — 1991. — cер. 1, вып. 3. — |
|
с. 41–44. |
33. |
Леонов Г. А. Об одном способе исследования глобальной устойчиво- |
|
сти нелинейных систем // Вестник Ленинградского университета. — |
|
1991. — сер. 1, вып. 4. — с. 11–14. |
34. Леонов Г. А. Построение специальной внешней меры Каратеодори для оценки хаусдорфовой размерности аттракторов // Вестник Санкт-Пе-
|
тербургского университета. — 1995. — cер. 1, вып. 4. — с. 24–30. |
35. |
Леонов Г. А. Верхние оценки хаусдорфовой размерности аттракто- |
|
ров // Вестник Санкт-Петербургского университета. — 1998. — сер. 1, |
|
вып. 1. — с. 19–22. |
36. |
Леонов Г. А. О нижних оценках ляпуновских экспонент и верх- |
|
них оценках хаусдорфовой размерности аттракторов // Вестник |
|
Санкт-Петербургского университета. — 1996. — сер. 1, вып. 4. — с. 18– |
|
24. |
37. |
Леонов Г. А. Оценки ляпуновских экспонент дискретных систем |
|
// Вестник Санкт-Петербургского университета. — 1997. — сер. 1, |
|
вып. 3. — с. 49–56. |
38. |
Леонов Г. А. О нижних оценках размерности аттракторов дискретных |
|
систем // Вестник Санкт-Петербургского университета. — 1998. — |
|
сер. 1, вып. 4. — с. 45–48. |
160 |
ЛИТЕРАТУРА |
39. Леонов Г. А. Частотный критерий существования странных аттракторов дискретных систем // Автоматика и телемеханика. — 1999. — № 5. — c. 113–122.
40. Леонов Г. А., Селеджи С. М. Системы фазовой синхронизации в аналоговой и цифровой схемотехнике. — СПб.: Невский диалект, 2002. — 110 c.
41. Леонов Г. А., Смирнова В. Б. Математические проблемы теории фазовой синхронизации. — СПб.: Наука, 2000. — 399 с.
42. Ляпунов А. М. Общая задача об устойчивости движения.. — Харьков, 1892. — 250 c.
43. Максвелл Д. К., Вышнеградский И. А. и Стодола A. Теория автоматического регулирования. Линеаризованные задачи. — М.: Изд-во АН
СССР, 1949. — 431 с.
44. Малкин И. Г. Теория устойчивости движения. — М.: Наука, 1966. — 525 с.
45. Массера Х. Л. К теории устойчивости. Сб. переводов “Математика”. — ИЛ, 1957. — т. 1, № 4, c. 81–101.
46. Монин А. С. О природе турбулентности. // Успехи физических наук, 1978. — т. 125, вып. 1, с. 97–122.
47. Натансон И. П. Теория функций вещественной переменной. — М.: Наука, 1974. — 480 с.
48. Неймарк Ю. И., Ланда П. С. Стохастические и хаотические колебания. — М.: Наука, 1987. — 423 с.
49. Немыцкий В. В., Степанов В. В. Качественная теория дифференциальных уравнений. — М.-Л.: ГИТТЛ, 1947. — 448 c.
50. Оселедец В. И. Мультипликативная эргодическая теорема. Характеристические показатели Ляпунова динамических систем // Труды Московского математического общества. — 1968. — т. 19. — с. 179–210.
51. Первозванский А. А. Курс теории автоматического управления. — М.: Наука, 1986. — 625 с.
52. Персидский К. П. О характеристических числах дифференциальных уравнений // Изв. АН Казахской ССР. Сер. матем. и мех. — 1947. — вып. 1.
53. Песин Я. Б. Характеристики размерностного типа для инвариантных множеств динамических сисем // Успехи математических наук. — 1988. — т. 43, вып. 4. — с. 95–128.
