Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике ДЛЯ СТУДЕНТОВ И ШКОЛЬНИКОВ ЕСТЬ
ZADANIA.TO
Рассмотрим неподвижные точки квадратичного отображения (7.2.4) и их линейную устойчивость. Неподвижные точки находятся из условия: которое соответствует пересечению кривой $f(x)$, определяемой (7.2.4б), и прямой $f=x$. Решение (7.2.12) имеет два корня: Устойчивость неподвижных точек определяется из линеаризованного отображения ( $\$ 3.3$ ). Подстазляя в (7.2.4) и удерживая только линейные члены, получаем где собственное значение $\lambda_{1}=f^{\prime}\left(x_{1}\right)$. Следовательно, $x_{1}$ — устойчивая (притягивающая) неподвижная точка, если $\left|\lambda_{1}\right|<1$, и неустойчивая (отталкивающая), если $\left|\lambda_{1}\right|>1$. Отсюда следует, что точка $x_{11}$ устойчива при или для $1 / 2<C<3 / 2$. Аналогичғо неподвижная точка $x_{10}$ устойчива при или для $-1 / 2<C<1 / 2$. Чтобы понять, что происходит при $C<-1 / 2$, нужно рассмотреть периодические точки периода 2 , которые находятся из условия: Это соответствует пересечению кривой с прямой $f=x$, как показано на рис. $7.11, a$ и б для $C>-1 / 2$ и $C<-1 / 2$ соответственно. Если $C>-1 / 2$, то оба пересечения соответствуют неподвижным точкам $x_{11}=1 / 2-$ и $x_{10}=0$, которые, очевидно, удовлетворяют и (7.2.16). В этом случае наклон $f_{2}(x)$ в нуле меньше единицы, так как а $\left|\lambda_{1}\right|<1$ для $C>-1 / 2$. Если же параметр $C$ становится меньше — $1 / 2$, то наклон увеличивается и рождается пара неподвижных точек отображения $f_{2}$ (или периодических точек отображения $f$, Рис. 7.11. Возникновение двух неподвижных точек $x_{2}$ с п периодом 2 для квадратичного отображения. рис. 7.11, б). Они удовлетворяют условиям Позже мы получим явные выражения для $x_{2+}$ и $x_{2-}$ как функций параметра $C$. Устойчивость $x_{2 \pm}$ можно исстедовать, как обычно, с помощью замены $x_{2, n}=x_{2 \pm}+\Delta x_{2, n}$ и линеаризации отображения. Для точки $x_{2-}$ находим где Аналогично для $x_{2+}$ : Видно, что наклоны (собственные значения) одинаковы для обеих точек. Как отмечено в п. 3.За, для канонических отображений Из рис. 7.11, б следует, что если значение $C$ чуть меньше $C_{1}=$ $=-1 / 2$, то $\left|\lambda_{2}\left(x_{2-}\right)\right|=\left|\lambda_{2}\left(x_{2+}\right)\right|<1$. Таким образом, как только неподвижная точка периода 1 становится неустойчивой, появляется пара устойчивых неподвижных точек с удвоенным периодом. Это явление иллюстрируется на рис. 7.12, где показана Рис. 7.12. Положение неподвижных точек $x_{10}, x_{11}, x_{2}+$ и $x_{2-}$ в зависимости от $C$. зависимость координат неподвижных точек от параметра С. Сплошная линия обозначает устойчивость точки, а штриховая — неустойчивость. Рождение пары устойчивых точек удвоенного периода при $C<-1 / 2$ является примером бифуркации удвоения (см. §7.1). На рис. 7.2, в показана также зеркальная бифуркация удвоения, соответствуюцая зеркальному отображению (7.2.7). При дальнейшем уменьшении параметра $C$ обе неподвижные точки периода 2 становятся неустойчивыми. Поскольку в их окрестности отображение $f_{2}(x)$ также квадратично, следует ожидать аналогичной бифуркации, при которой возникают четыре неподвижные точки периода 4 и т. д. Означает ли это, что при любом значении параметра имеется либо устойчивое решение, либо (устойчивая) бифуркация? Оказывается, что нет. Последовательность бифуркаций удвоения периода кончается, достигая бесконечного периода при конечном значении параметра $C=C_{\infty}$. За этим значением лежат области хаоса. Метод ренормализации. Для более глубокого понимания явления последовательных бифуркаций рассмотрим ренормализацию отображения с одним максимумом при переходе от одной бифуркации к следующей. Строгая теория ренормализации $[83,122]$, применимая к любым отображениям с одним максимумом, основана на решении некоторых функциональных уравнений. Соответствующие математические методы сложны и выходят за рамки этой монографии. Вместо этого, следуя Хеллеману [180-182], мы используем более простой, хотя и менее общий, метод локальной аппроксимации отображения с одним квадратичным максимумом. Это позволяет построить приближенную теорию ренормализации на основе алгебраических, а не функциональных уравнений. Для квадратичного отображения (7.2.4), как мы видели выше, устойчивая неподвижная точка $x_{10}$ периода 1 возникает с уменьшением $C$ при $C=$ $=C_{0}=1 / 2$, а устойчивые точки периода 2 — при $C_{1}=-1 / 2$. Последние получаются из (7.2.18): где Подставляя в (7.2.4), получаем где и начальные условия выбираются таким образом, чтобы при четных $n$ величина $\Delta x$ была близка к $x_{2+}$, а при нечетных — к $x_{2-}$. Исключая $\Delta x_{n+1}$ и удерживая члены по $\Delta x$ до квадратичных включительно, получаем для четных $n$ Путем изменения масштаба приводим (7.2.24) к виду где а $b$ определяется согласно (7.2.21Е). Так как (7.2.26) имеет тот же самый вид, что и исходное отображение (7.2.4), его неподвижные точки испытывают бифуркацию удвоения периода при $C^{\prime}=-1 / 2$. С учетом (7.2.27) это значение соответствует при котором возникают неподвижные точки периода 4 . Из (7.2.21) следует, что $x_{2+} \approx 0,466$ и $x_{2-} \approx-0,241$ при $C=C_{2}$. Отметим, что точка $x_{2+}$ лежит вблизи экс.ремума $f$ при $x^{*}=-C 2$; этот факт мы используем ниже при нахождении закона подобия (7.2.36). Последовательность бифуркаций сходится к. предельному значению $C^{\prime}=C=C_{\infty}$, которое находится из (7.2.27): отсюда Численное моделирование дает $C_{\infty} \approx-0,78497$, что находится в хорошем согласии с результатом ренормализации. Скорость сходимости $C$ к $C_{\infty}$ можно приближенно найти, предполагая, что асимптотически сходимость происходит по закону ${ }^{1}$ ) Подставляя это выражение в (7.2.27) и замечая, что $k$-я бифуркация по $C^{\prime}$ соответствует ( $k+1$ )-й по $C$, т. е. получаем Численное решение точного уравнения ренормализации, найденное впервые Фейгенбаумом [122], дает для отображения с одним квадратичным максимумом $\delta \approx 4,6692$. Закон подобия (7.2.30) можно представить в виде Постоянная $\delta$ не зависит от выбора параметра. Действительно, введем новый параметр $P=g(C)$. Считая, что $g$ — обратимая функция, разложим $P(C)$ вблизи $C_{\infty}$ : Решая это уравнение относительно $C_{k}$ и подставіяя решение в (7.2.32), находим универсальное соотношение Фактически из точной теории ренормализации следует, что $\delta-$ универсальная константа для всех отображений с одним квадратичным максимумом. Принимая, чтс асимптотический закон (7.2.30) справедлив и для $k=0$ и $k=1$, получаем полезную оценку для предельного значения параметра (точки сгущения): где использовано приближенное значение (7.2.31) для $\delta$. Точная теория ренормализации дает $\alpha \approx-2,5029$. Параметр $\alpha$ определяет (асимптотически) изменение масштаба $x$ при последовательных бифуркациях. Иначе говоря, при увеличении в $\alpha$ раз вблизи $x^{*}=-C / 2$ очередная бифуркация будет выглядеть точно так же, как и предыдущая. В частности, расстояние между неподвижными точками при последовательных бифуркациях вблизи $x^{*}$ подчиняется закону подобия Это соотношение справедливо для ветви, идущей от $x_{k+}$. Потученное соотношение выполняется для значений $x_{k}$ на ветви, идущей от $x_{k-}$. Таким образом, для половины неподвижных точек периода $2^{k}$ масштаб изменяется как $\alpha$, а для остальных как $\alpha^{2}$ [см. (7.2.38)]. Рис. 7.13 поясняет это поведение. На рис. 7.13 , a отрезок $\Delta x_{1}$ вблизи минимума отображения переходит в отрезок $\left|\Delta x_{2}\right| \propto\left(\Delta x_{1}\right)^{2}$. Однако затем $\Delta x_{2}$ переходит в $\Delta x_{3}$ снова вблизи минимума, причем $\left|\Delta x_{3}\right| \propto\left|\Delta x_{2}\right|$. При каждом отображении знак $\Delta x$ изменяется. На рис. 7.13 , б показана полная структура неподвижных точек Рис. 7.13. Иллюстрация масштабного преобразования по $x$ для квадратичного отображения. для трех первых бифуркаций. Цифры в круглых скобках показывают последовательность движения для траекторий с периодом 4 и 8 , начиная с верхней точки. Отметим также, что расстояния между соответствующими ветвями, идущими от $x_{k+}$ и от $x_{k-}$, отєличаются в $\alpha$ раз. При любом обратимом преобразовании от $x$ к новой переменной y соотношения (7.2.36) и (7.2.37) сохраняются. Следовательно, Рис. 7.14. Последовательность бифуркаций удвоения периода для квадратичного отображения в двойном лотарифмическом масштабе (по данным работы [417]). константа $\alpha$ также является универсальной, в том же смысле, что и $\delta$. Помимо последовательности бифуркаций, при уменьшении $C$ имеется зеркальная последовательность бифуркаций согласно $(7.2 .6)-(7.2 .8)$. Для нее $C$ увеличивается, начиная с $1 / 2$, и стремится к своей точке сгущения, которая определяется из (7.2.8) и $(7.2 .29)$ : Для отображения (7.2.5) зеркальная последовательность бифуркаций приходится как раз на обычно рассматриваемый интервал $0<\mu<4$. Используя (7.2.8) и соотношение $\bar{\mu}=2 \bar{C}$, находим, что первая бифуркация наступает при $\bar{\mu}_{1}=3$, а точка сгущения равна что хорошо согласуется с численным значением ${ }^{\circ} \bar{\mu}_{\infty}=3,5700$. Как известно, спектр периодического движения дискретный. Когда происходит бифуркация удвоения периода, в спектре появляются субгармоники основной частоты отображения. Для дальнейшего анализа введем непрерывное время $t$ и обозначим через $x^{(n)}(t)$ решение после $n$-й бифуркации, а через $T_{n}$ — его период. Чтобы найти спектр, нам нужно общее соотношение для положения аттракторов. Из рис. 7.13, б и пояснений в тексте следует, что ветви бифуркаций делятся на две группы, расстояния между которыми удовлетворяют рекуррєнтному соотношению: где $t=0$ соответствует верхней неподвижной точке. С помощью сдвига на половину периода интегрирование можно провести по периоду $T_{n}$ : Для четных $l=2 k$ при $n \gg 1$ Это означает, что фурье-амплитуда для данной частоты остается неизменной при всех последующих бифуркациях. Для нечетных $l=2 k+1$ подстановка (7.2.38) в (7.2.39) дает Сдвигая пределы во втором интеграле, имеем Подставляя сюда разложение Фурье получаем рекуррентное соотношение де В пределе больших $n$, полагая $2 k+1=\xi$ и $2 k^{\prime}+1=\xi^{\prime}$ и заменяя сумму в (7.2.42) интегралом по $\xi^{\prime}$, находим Поэтому для амплитудного спектра вблизи точки сгущения выполняется универсальный закон подобия ${ }^{1}$ ): Для приближенного значения (7.2.35) для $\alpha$ имеем $\gamma=5,79$; точное значение $\alpha=-2,5029$ дает $\gamma=6,57$. Согласно (7.2.44), чтобы получить амплитуду новой гармоники (вдали от точки бифуркации), появившейся в результате $(n+1)$-й бифуркации, нужно взять уменьшенное в 6,57 раза значение амплитуды гармоники, появившейся в результате $n$-й бифуркации. Это теоретическое предсказание подтверждается экспериментально (см. работу [123] и рис. 7.33). Мы обсудим эти-экспериментальные данные в $\S 7.4$. Другие периодические траектории. При уменьшении параметра $C$ от $C_{0}$ до $C_{\infty}$ происходит бесконечное число бифуркаций удвоения. При $C<C_{\infty}$ также имеются области периодического движения. Периодические траектории рождаются здесь парами (устойчивая и неустойчивая) в результате тангенциальной бифуркации. В качестве примера ${ }^{2}$ ) рождения траектории с периодом 3 на рис. 7.15 отложена зависимость для $f(x)$ согласно (7.2.4б). Неподвижные точки периода 3 удовлетворяют уравнению Для $C$ незначительно больше $C_{0}^{(3)}$ это уравнение имеет только два (действительных) корня, которые дают неподвижные точки периода 1 (рис. $7.15, a$ ). Если же $C<C_{0}^{(3)}$, то уравнение имеет шесть корней (рис. 7.15, б), которые ссответствуют двум разным траекториям с периодом 3. Легко показать, что одна траектория устойчива (рис. 7.15, б, темные кружки), а другая-неустойчива (светлые кружки). Критическое значение параметра $C_{0}^{(3)}$, при котором рождаются эти траектории, равно Рис. 7.15. Рождение пары траекторий периода 3 при тангенциальной бифуркации. Поскольку вблизи неподвижных точек отображение $f_{3}$ лока.тьно квадратично, то можно ожидать, что при дальнейшем уменьшении $C$ будут возникать бифуркации удвоения с периодами $3,6,12,24 \ldots$ Их точку сгущения $C_{\infty}^{(3)}$ можно найти с помощью описанной выше ренормализации или же путем численного моделирования. В последнем случае $C_{\infty}^{(3)} \approx-0,92475$. Существуют и зеркальные бифуркации удвоения при $\bar{C}=1-C$. Точно так же с помощью отображений $f_{4}, f_{5}, f_{6}, \ldots$ можно найти устойчивые траектории с основным периодом $n=4,5,6$, … I. каждой из этих траекторий существует своя последовательность бифуркаций удвоения с начальной точкой $C_{0}^{(n)}$ и точкой сгущения $C_{\infty}^{(n)}$. Порядок, в котором появляются основные периоды при уменьшении параметра, является универсальным для всех отображений с одним квадратичным максимумом. Например, первые шесть периодов упорядочены при уменьшении $C$ следующим образом: 1,6 , $5,3,5,6,4,6,5,6$. Можно также найти и их общее число [296]. Например, имеются 202 основнее траектории с периодом до 11 включительно; их упорядочение по параметру исследовано Метрополисом и др. [300]. Недавно Гейзал и Нирветберг [151] показали, что для всех этих траекторий ренормализация имеет универсальную структуру.
|
1 |
Оглавление
|