Джеймс Глейк - Хаос. Создание новой науки

К 1980 г. он понял: что-то идет не так, разрушая его модель. Выяснилось, что ключевая роль принадлежала виду, который не был учтен при построении модели, — муравьям. Некоторые коллеги ученого подозревали, что все дело в необычной зимней погоде, другие — в летней или прочих неучтенных факторах. Шаффер принялся обдумывать, как бы учесть дополнительные параметры. Все же ученый казался удрученным. Его аспиранты утверждали, что тем летом он выглядел мрачным и напряженно работал.

Потом все изменилось. Случайно обнаружив препринт статьи о химическом хаосе в сложном лабораторном эксперименте, Шаффер почувствовал, что ее авторы столкнулись с тем же явлением, что и он сам. Выявить и описать десятки продуктов реакции в пробирке оказалось так же невозможно, как и учесть все многообразие видов в горах Аризоны. Все же химикам удалось достичь успеха там, где эколог потерпел крах. Шаффер принялся читать о реконструкции фазового пространства, познакомился наконец с работами Лоренца, Йорка и других исследователей. Университет Аризоны выделил средства на серию лекций «Порядок внутри хаоса». Читать их пригласили Гарри Суинни, а он мог многое рассказать о практических опытах. Когда Суинни объяснил, что вызывает хаос в химии, продемонстрировал странный аттрактор и заявил, что «это реальные данные», Шаффер покрылся холодным потом.

«Внезапно я понял, что это судьба», — вспоминал позже Шаффер, которому предстоял год академического отпуска. Он отозвал свою заявку из Национального научного фонда, куда обращался с просьбой о финансировании, и начал все снова. Высоко в горах Аризоны популяция муравьев росла и уменьшалась, пчелы с жужжанием кружились в воздухе, облака медленно плыли по небу, а Шаффер постигал новую науку. Он больше не мог работать как прежде.

* Компьютерная программа, воспроизводящая систему Мандельбро, нуждается в разъяснении нескольких существенных деталей. Главный ее механизм состоит в том, что выбирается начальное комплексное число и к нему применяется арифметическое правило. Для рассматриваемой ниже системы правило таково: z → z ² + с , где z начинается с нуля, а с представляет собой комплексное число, соответствующее тестируемой точке. Итак, возьмем нуль, умножим его на самого себя, прибавим начальное число; взяв результат (начальное число), умножим его на самое себя и прибавим начальное число; возьмем новый результат, опять умножим его на самого себя и прибавим начальное число. Арифметика с комплексными числами ведет нас прямо вперед. Комплексное число состоит из двух частей, например: 2 + Зi (местоположение точки: 2 к востоку и 3 к северу на комплексной плоскости). Чтобы сложить два комплексных числа, надо лишь сложить действительные части для получения новой действительной части и мнимые — для получения новой мнимой части:

Чтобы перемножить два комплексных числа, нужно умножить каждую часть одного из них на каждую часть другого (но правилам перемножения двучленов) и сложить получившиеся четыре результата. Поскольку і, умноженное на самое себя, дает -1, то в силу первоначального определения мнимых чисел один член результата переходит в другой:

Чтобы прекратить движение по петле, программа должна отслеживать текущий итог. Если результат стремится к бесконечности, все более и более удаляясь от центра плоскости, выбранная точка не принадлежит к системе. В том случае, когда итог превышает 2 или становится меньше -2 либо в действительной, либо в мнимой части, результат, бесспорно, стремится к бесконечности и работа программы может продолжаться. Коль скоро она выполняет одни и те же вычисления много раз, не превышая 2, точка является частью системы. (Число раз зависит от степени увеличения. Для масштаба, доступного персональному компьютеру, ста или двухсот раз часто бывает достаточно, а тысяча повторений дает полную гарантию.) Программа должна повторить данный процесс для каждой из тысяч точек решетки. Масштаб можно увеличить. Затем программа должна показать полученный результат. Точки, входящие в систему, могут быть обозначены черным цветом, а не принадлежащие к ней — белым. Для получения более живого изображения белый цвет можно заменить оттенками других цветов. В частности, если итерация прекращается после десяти повторений, программа должна выдать красную точку, после двадцати — оранжевую, после сорока — желтую и т. д. Выбор цветов и момент остановки расчета точек программист может выбрать сам. Цвета надлежащим образом обозначают контуры, оставшиеся за пределами системы.

Перевод С. Я. Маршака.

В Северной Америке ураганами принято называть тропические циклоны. (Примеч. ред.)

Удобства ради в данной весьма абстрактной модели численность особей выражена через дробь, которая больше нуля, но меньше единицы, причем нуль обозначает вымирание, а единица — наиболее высокую численность животных (в данном случае рыб), достижимую в пределах популяции. Итак, начнем: произвольно выберем значение параметра r , скажем 2,7, и начальную численность популяции, к примеру 0,02. Отнимем от единицы 0,02 и получим 0,98, умножим 0,98 на 0,02 и получим в итоге 0,0196. Теперь умножим полученный результат на 2,7 и получим 0,0529. Крошечная начальная численность популяции выросла более чем в два раза. Повторим процедуру, используя только что полученную численность особей в качестве исходных данных, и получим 0,1353. С небольшим калькулятором, в который можно ввести определенную программу, для получения такого результата нужно лишь нажимать одну и ту же клавишу снова и снова. Популяция увеличивается до 0,3159, затем до 0,5835; 0,6562 — рост численности замедляется. Далее, по мере того как смертность «догоняет» воспроизводство, численность достигает 0,6092; 0,6428, 0,6199, 0,6362, 0,6249. Значения в числовом ряду скачут: то возрастают, то уменьшаются. Впрочем, заканчивается он строго определенным значением: 0,6328, 0,6273, 0,6312, 0,6285, 0,6304, 0,6291, 0,6300, 0,6294, 0,6299, 0,6295, 0,6297, 0,6296, 0,6297, 0,6296, 0,6296, 0,6296, 0,6296, 0,6296, 0,6296, 0,6296. Это явный успех. Когда все расчеты выполнялись вручную, и даже во времена механических счетных машинок с ручным вводом, дальше подобных вычислений дело не шло.

Скажем, при r = 3,5 и начальной численности популяции 0,4 Мэй увидел следующий числовой ряд: 0,4000; 0,8400; 0,4704; 0,8719; 0,3908; 0,8332; 0,4862; 0,8743; 0,3846; 0,8284; 0,4976; 0,8750; 0,3829; 0,8270; 0,4976; 0,8750; 0,3829; 0,8270; 0,5008; 0,8750; 0,3828; 0,8269; 0,5009; 0,8750; 0,3828; 0,8269; 0,5009; 0,8750 и т. д.