Отображения на примере школьной задачи

Кузнецов С. П., Кузнецов А. П.
Саратовская группа теоретической нелинейной динамики
Источник: http://sgtnd.narod.ru/wts/rus/kruzh/maps.htm

Практически в любом сборнике олимпиадных задач по физике можно найти задачу о бесконечной цепочке резисторов. Она формулируется так: чему равно сопротивление изображенной на рис. 1 цепочки резисторов, которая состоит из одинаковых звеньев (рис.1)?

Рис. 1

Сопротивления всех резисторов одинаковы, и мы будем полагать их равными единице. Традиционное решение этой задачи основано на соображении о том, что добавление еще одного звена к бесконечной системе не изменяет ее, т.е. данная схема эквивалентна показанной на рис. 2.

Рис. 2

В результате приходим к схеме, показанной на рис. 3.

Рис. 3

Здесь x₀ — искомое сопротивление. Нетрудно получить, что

, или .

Тогда

.

Откуда следует ответ к задаче:

Описанный метод в журнале «Квант» в статье, посвященной методам вычисления сопротивлений различных схем, назвали методом Ийона Тихого — по имени известного героя Станислава Лема, сумевшего разместить постояльца в полностью заполненной гостинице с бесконечным числом номеров.

Такое решение, хотя и весьма изящно, но оставляет некоторое чувство неудовлетворения, поскольку остаются открытыми некоторые вопросы. В первую очередь: какое отношение к реальности имеет эта задача и это решение? Ведь реальная цепочка будет содержать конечное число звеньев. Можно сформулировать вопрос так: если мы можем измерять сопротивление с заданной точностью, то сколько звеньев должна содержать цепочка, чтобы считаться бесконечной? Кроме того, если вдуматься, то заранее не очевидно, что добавление новых звеньев в реальной конечной цепочке будет приближать результат измерения сопротивления к значению x₀. Вдруг малые погрешности, которые вносят реальные резисторы, будут вносить нарастающий вклад, ведь цепочка бесконечная.

Оказывается, на все эти вопросы можно ответить, привлекая изящный математический аппарат дискретных отображений, являющийся частью современной теории динамических систем.

Итак, обратимся к нашей схеме, отсчитаем конечное число звеньев n и n+1 от правого конца (рис. 4).

Рис. 4

Из рисунка хорошо видно, что схема эквивалентна такой:

Рис. 5

Отсюда легко получаем

Это и есть простейший пример дискретного отображения. В общем виде одномерное дискретное отображение задается соотношением

(В нашем случае ) Мы назвали отображение одномерным, поскольку в него входит одна переменная — x.

Дискретное отображение x_n+1=f(x_n) является, по-видимому, простейшим примером динамической системы. Смысл этого термина раскрывается просто: отображение x_n+1=f(x_n) по заданному начальному значению x₁ позволяет определить все последующие значения переменной x₂, x₃ и т.д. Действительно:

Свойства дискретных отображений удобно иллюстрировать на так называемой итерационной диаграмме. Для ее построения надо прежде всего изобразить график функции f(x). В нашем случае это перевернутая гипербола, смещенная влево на две единицы и на одну единицу вверх (рис. 6). (Заметим, что поскольку сопротивления положительны, то смысл имеет лишь кусок графика в первой четверти.)

Рис. 6 — Итерационная диаграмма

На графике изображается еще и биссектриса. Задавшись теперь начальным значением x₁, можно найти x₂=f(x₁) по графику. Затем это значение переносится на биссектрису и процедура повторяется. Возникает своеобразная лесенка, иллюстрирующая ход итераций.

Как видно из графика, наше отображение имеет, как говорят, неподвижную точку, т.е. такую, для которой x₀=f(x₀). Она как раз соответствует решению задачи по методу Ийона Тихого.

По рисунку получается, что последовательные итерации сходятся к неподвижной точке. Как это доказать более строго? Каков характер сходимости (быстрый, медленный)? На эти вопросы можно ответить в общей формулировке.

Для этого исследуем поведение системы в случае, когда значений переменной близки к предельному значению x₀. Положим поэтому и , где знаком «тильда» сверху обозначены малые добавки к x₀. Тогда из имеем

Таким образом, если имеется некоторая маленькая добавка к значению x₀, то после первой итерации она умножается на постоянное число C= f '(x₀), после второй — на C², после третьей — на C³ и т.д. Это означает, что переменная x приближается к неподвижной точке по закону геометрической прогрессии с показателем C. Отметим, что наше рассмотрение на итерационной диаграмме соответствует тому, что мы аппроксимируем f(x) касательной (вспомните геометрический смысл производной!) в окрестности x₀ Соответствующая итерационная диаграмма и дает геометрическую прогрессию (рис. 7).

Рис. 7. Слева — аппроксимация нашего графика прямой с помощью касательной в неподвижной точке. Справа — итерационная диаграмма. Приводящая к геометрической прогрессии.

Теперь по свойству геометрической прогрессии автоматически получаем, что если

• , то итерации сходятся;
• , то итерации расходятся.

Это, как говорят, позволяет судить об устойчивости неподвижной точки.

В первом случае неподвижную точку называют устойчивой, а во втором — неустойчивой.

Сходящийся процесс мы уже изобразили на итерационной диаграмме, а расходящийся может выглядеть так (рис. 8):

Рис. 8.

Вернемся от отображения общего вида к нашему случаю. Тогда

.

Итак, как мы видим C=0.145900. Это означает, что итерации сходятся, причем поскольку f '(x₀) мало, значит, сходятся очень быстро. В этом убеждаемся, итерируя соотношение

Результат представлен в следующей таблице.

Число звеньев в цепочке, N	xN	%
1	x1=1
2	x2=2/3≈0,666667	7,86
3	x3=5/8≈0,625000	1,12
4	x4=13/21≈0,619048	0,16
5	x5=34/55≈0,618182	0,02
…	…	…
∞	x0=	0

Мы убеждаемся, что увеличение числа звеньев действительно приводит в неподвижную точку. Из таблицы видно, что оказывается цепочку всего из трех звеньев можно считать бесконечной с точностью до 1%, а из пяти — уже с точностью всего 0,02%!

Устойчивость неподвижной точки этого отображения снимает и еще один физический вопрос: о том, не «испортят» ли возможные дефекты в цепочке результат нашего решения. Продумайте сами, почему это так и изобразите соответствующие итерационные диаграммы.

Нам осталось сказать, что в силу большой важности величины f '(x₀) она носит специальное название — мультипликатор, и обозначается обычно .