Устойчивость решений дискретных систем (5) В дискретных динамических системах могут существовать частные решения, представляющие собой стационарные, периодические, квазипериодические и хаотические последовательности, которые мы будем обозначать x n *. Рассмотрим первое частное решение дискретной системы – стационарное, или решение в виде неподвижной точки. Определение. Точка x* называется неподвижной точкой отображения F (4), если выполняется следующее условие: x* = F(x*). Это условие означает, что при графической итерации отображения движение при попадании в x* прекращается. Итерации вновь и вновь помещают очередное x n в x*. Решение в виде неподвижной точки еще называется состоянием покоя или состоянием равновесия.

Помимо самого факта существования в дискретных системах решений в виде неподвижных точек (состояний покоя или равновесия системы), важную роль играет информация об их устойчивости или неустойчивости. Слегка толкнем шарик и пронаблюдаем за движением. После совершения нескольких затухающих колебаний шарик вновь займет прежнее положение на дне ямки. Положение равновесия устойчиво: малые возмущения исходного состояния затухают во времени. При любом сколь угодно малом отклонении шарика от состояния равновесия он скатится с вершины. Положение равновесия неустойчиво: малые возмущения исходного состояния нарастают во времени. Устойчивость какого-либо состояния (движения) динамической системы определяется просто: введем небольшое отклонение (возмущение) динамической системы от исследуемого состояния и проанализируем, каким будет ее дальнейшее поведение. Если со временем система вернется в исходное состояние (возмущение затухает), то такое состояние называется устойчивым. Если начальное отклонение нарастает со временем - состояние неустойчиво.

Найдем условия устойчивости неподвижной точки x *. В связи с тем, что проблема устойчивости связана с анализом реакции системы на малое возмущение ее состояния, на первом этапе она может быть исследована в рамках линейного приближения. Пусть x = x * +, - отклонение от состояния равновесия, малая по сравнению с x * величина. Тогда отображение (4) примет вид: (6) Применяя для правой части (6) разложение в степенной ряд в окрестности неподвижной точки и отбрасывая члены порядка малости выше первого, находим: (7) С учетом (5), из соотношения (7) получаем (8) Производные берутся в точках частного решения, в данном случае в неподвижной точке x*.

Введем обозначение Здесь индексы i,k = 1,…, N – размерность отображения (4). Тогда выражение (8) можно переписать в координатной форме (или покомпонентно): (9)(9) Данное линейное точечное отображение есть линейное дискретное уравнение в вариациях. В матричной форме (9) можно записать в следующем виде: (10) Матрица А с элементами a ik является квадратичной и называется матрицей линеаризации или матрицей Якоби. Устойчивость неподвижной точки отображения определяется мультипликаторами k, которые являются собственными значениями матрицы А, или корнями характеристического уравнения: (11) E – единичная матрица.

Неподвижная точка x * исходного отображения (4) является устойчивой, если все мультипликаторы k удовлетворяют условию | k | < 1. Если среди мультипликаторов имеются такие, для которых | k | > 1, то неподвижная точка отображения (4) будет неустойчивой.

Устойчивость неподвижных точек одномерного отображения Одномерное отображение: (12) Пусть x * - неподвижная точка отображения. Введем малое возмущение : Тогда отображение (12) примет вид: Используя выражения (7), (8), приходим к новому одномерному линейному отображению, описывающему эволюцию малого возмущения неподвижной точки, которое имеет вид: (13) Из сопоставления выражений (10), (11) и (14) следует, что собственное значение матрицы линеаризации A или мультипликатор неподвижной точки одномерного отображения (12) есть = f (x*). Следовательно, условие устойчивости неподвижной точки x* сводится к выполнению неравенства | f (x*)| 1, неподвижная точка неустойчива. (14)

Возможные типы поведения неподвижной точки в зависимости от величины и знака производной f (x*). С геометрической точки зрения значения f (x*) – это тангенс угла наклона графика функции последования f( ) в точке x*. Эволюции малого возмущения n соответствует сжатие либо растяжение отрезка, соответствующее отклонению от x*. В малой окрестности неподвижной точки функцию последования можно заменить прямой А, касательной к f(x) в точке x*.

Устойчивый узел Типы неподвижных точек в одномерном отображении Состояние нейтрально. Переходные значения параметра называются критическими. Им соответствуют точки бифуркации. Неустойчивый узел Устойчивый обратный узел Неустойчивый обратный узел Если возмущение меняет знак при каждой итерации, то к названию неподвижных точек добавляют слово «обратный».

Устойчивость неподвижных точек двумерного отображения (16) (15) (16) (17) (18)

(19) (17) (20) (21) (22) (20) (18)(16) (19) (20) (21)

Чтобы найти мультипликаторы неподвижных точек отображения или собственные значения матрицы линеаризации, необходимо найти корни характеристического уравнения (11). Для двумерного отображения оно имеет вид: (22) Характеристическое уравнение (22) можно переписать в виде (23) Корни характеристического уравнения (23) (24) определяют решение системы (20), следовательно, и поведение фазовых траекторий в окрестности неподвижной точки отображения (15).

Типы неподвижных точек в двумерном отображении Как видно из выражения (24), собственные значения 1,2 могут быть действительными числами или комплексно сопряженными. Как было показано в случае одномерного отображения, характер устойчивости неподвижной точки меняется при | 1,2 | = 1. Поэтому удобно характеризовать неподвижную точку, изображая расположение собственных чисел 1 и 2 на комплексной плоскости относительно единичной окружности. Im Re Im 1 Im 2 Re 1,2 Пара комплексно сопряженных значений | 1,2 | < 1, имеющие равные действительные части, а их мнимые части различаются только знаком. Положение точек внутри окружности означает, что

Случай действительных 1 и < -1 и 2 < -1 – неустойчивый обратный узел. Траектория расходится от неподвижной точки по обоим собственным векторам, причем отклонение каждый раз меняет свой знак < 1 < 0 и 2 < -1 – обратное седло. Траектория сходится по одному из собственных векторов, но расходится по другому. Отклонение меняет знак при каждой итерации < 1 < +1 и 2 < -1 – неориентируемое седло 1. По одному из направлений траектория расходится, причем отклонение меняет знак на каждой итерации. По другому направлению отклонение монотонно убывает > +1 и 2 < -1 – неустойчивый неориентируемый узел. Траектория расходится от неподвижной точки по обоим направлениям, причем по одному из них отклонение меняет знак на каждой итерации < 1 < 0 и -1 < 2 < 0 – устойчивый обратный узел. Меняя знак отклонения на каждой итерации, возмущенная траектория сходится к неподвижной точке по любому направлению.

6. 0 < 1 < +1 и -1 < 2 < 0 – устойчивый неориентируемый узел. То же, что и 5-й случай, но по одному из собственных направлений отклонение монотонно убывает, не изменяя своего знака > +1 и -1 < 2 < 0 – неориентируемое седло 2. То же, что и 3-й случай, но отклонение меняет знак на каждой итерации для устойчивого направления и 0 < 2 < +1 – устойчивый узел. Траектория монотонно сходится к неподвижной точке по любому направлению > +1 и 0 < 2 < +1 – седло. Траектория монотонно сходится к неподвижной точке по одному направлению, но разбегается по другому > +1 и 2 > +1 – неустойчивый узел. Траектория монотонно удаляется от неподвижной точки по любому направлению. Отклонение не меняет знака.

Таким образом, при действительных 1 и 2 мы имеем 2 типа неподвижных точек: 1) узел, когда оба собственных значения по модулю либо больше, либо меньше единицы, 2) седло, когда одно из собственных значений по модулю больше единицы, а другое – меньше. В зависимости от знаков 1 и 2 к названию неподвижной точки может добавляться характеристика типа сходимости: «обратный», если оба мультипликатора отрицательны, либо «неориентируемый», если 1 и 2 имеют разные знаки.

Случай комплексно сопряженных 1 и 2. В данном случае неподвижная точка является фокусом. Если то неподвижная точка называется устойчивым фокусом. В этом случае точки итерации исходного отображения (15) лежат на скручивающейся спирали с центром в неподвижной точке, которая в этом случае устойчива. Если r > 1, то спираль раскручивается, удаляясь от неустойчивой неподвижной точки, которая в данном случае является неустойчивым фокусом. Характер сходимости возмущенной траектории к неподвижной точке меняется в зависимости от знака Re 1,2. Случай r = 1 соответствует пограничной ситуации, когда возмущенная траектория не удаляется и не приближается к неподвижной точке, вращаясь вокруг нее с некоторой угловой скоростью. Устойчивый фокус Неустойчивый фокус

Рассмотрим 2-й тип частного решения дискретной системы – периодическое решение. Определение. Последовательность точек x * 1, x * 2, …, x * m называется циклом периода m или m-циклом точечного отображения, если они удовлетворяют условиям x * n+m = F (x * n ), x * 2 = F(x * 1 ), x * 3 = F(x * 2 ), …, x * 1 = F(x * m ), причем никакие два элемента в наборе x * 1, x * 2, …, x * m не совпадают. Точки цикла x * 1, x * 2, …, x * m называют иногда m-кратными неподвижными точками и для них можно записать: x * 1 = F(x * m ) = F(F(F…F(x * 1 )…)) = F (m) (x * 1 ). Из приведенных определений следует, что неподвижная точка отображения является циклом периода 1 (когда m = 1). Циклы отображения и их устойчивость

Устойчивость m-цикла дискретного отображения можно определить, исследовав на устойчивость неподвижные точки отображения: (25) По аналогии с алгоритмом анализа на устойчивость однократной неподвижной точки (см. формулы (6)-(11)) получаем, что мультипликаторы km m-кратной неподвижной точки отображения (25) или m-цикла отображения (4) определяются как собственные значения характеристического уравнения: (26) Матрица линеаризации A m является m-периодичной и для нее справедливо следующее равенство: (27) Условие устойчивости m-цикла : (28)

В случае одномерного отображения матрица линеаризации A m для m-кратной неподвижной точки отображения или цикла периода m имеет вид: Это означает, что производная m раз примененного отображения в точке цикла – это произведение производных исходного отображения по всем точкам цикла. Таким образом, условие устойчивости цикла периода m одномерного отображения определяется неравенством: Здесь f (0) (x) = x. В случае, когда производная m раз примененной функции вычисляется для одной из точек x * i цикла периода m, получим (29) (30) (31) (32)

Сверхустойчивый цикл Из соотношения (32) видно, что устойчивость цикла в целом определяется совокупными свойствами всех его точек. При итерации на одних из них начальное отклонение может локально нарастать, тогда как на других – уменьшаться. Однако имеется особый случай, когда свойства одной точки определяют устойчивость цикла в целом, а именно, если для одной из точек x * i цикла выполняется условие то, очевидно, автоматически равно нулю и произведение производных по всем точкам. Это означает, что малое начальное отклонение от такого цикла полностью затухнет не более чем за m итераций, что и определяет название такого цикла. Поскольку с точки зрения геометрии функции последования f (x n ) данное условие означает наличие экстремума (минимум, максимум, либо точка перегиба), то можно сказать, что сверхустойчивый цикл содержит хотя бы одну критическую точку функции последования.

Пример. Рассмотрим логистическое отображение. Это одномерное квадратичное отображение, определяемое следующим образом: (33) где – управляющий параметр, а x n принадлежит интервалу [0, 1]. Данное отображение было введено еще в 1845 г. П. Ферхюльстом для описания динамики популяций в замкнутой среде. Относительная численность особей x n+1 в (n + 1)- й год пропорциональная численности особей в предыдущий год ( x n принимает значения от 0 до 1 и отражает численность популяции в n -м году), а также свободной части жизненного пространства, которая пропорциональна (1 - x n ), т.е. Положительный параметр характеризует скорость роста популяции. /4 Другой пример дает задача о банковских сбережениях при стабилизирующимся росте процента. Как было установлено, в частности, М. Фейгенбаумом, при варьировании параметра данное отображение демонстрирует довольно сложное поведение, которое становится хаотическим при больших. По мере увеличения параметра крутизна параболы плавно растет и вместе с этим будет меняться и устойчивый режим отображения. x*x*

Найдем неподвижные точки логистического отображения и исследуем их устойчивость. Неподвижная точка периода 1 определяется условием Получаем две неподвижные точки: Значение мультипликатора неподвижных точек находим из решения характеристического уравнения (для одномерного случая f (x*) - = 0): Получаем

x*x* 1) Неподвижная точка x* = 0. Ее мультипликатор =. Данная точка является устойчивой при < 1 и становится неустойчивой при 1 < < 3.

x*x* 2) Неподвижная точка x* = 1- 1/. Ее мультипликатор = 2 -. Условие устойчивости для данной неподвижной точки: | | = |2 – | < 1. Следовательно, точка устойчива при 1 3.

x1*x1* x2*x2* 3) При µ = 3 точка x* становится неустойчивой и в системе рождаются две устойчивые точки x* 1 и x* 2. Для исходного отображения (33) они образуют цикл периода 2. Эти точки удовлетворяют следующим соотношениям: что дает

Условие устойчивости цикла периода 2 определяется с помощью неравенства (32). Вычислим сначала производные функции последования в точках x * 1 и x * 2 Согласно (32) условие устойчивости для цикла периода 2 есть Подставляя в данное неравенство значения неподвижных точек, получим: Таким образом, цикл периода 2 отображения (33) будет устойчив в данном диапазоне значений параметра µ. При µ = … 2-цикл потеряет свою устойчивость и в системе родится цикл периода 4.