вместо X. Естественно задать вопрос: какова вероятность того, что в данный момент времени концентрация вещества Х имеет то или иное значение? Ясно, что эта вероятность флуктуирует, поскольку флуктуирует число столкновений между молекулами различных веществ, участвующих в реакции. Нетрудно выписать уравнение, описывающее, как изменяется распределение вероятности Р (X, t) в результате процессов рождения и уничтожения молекул X. Для равновесных или стационарных систем это распределение вероятности можно вычислить. Начнем с результатов, которые удается получить для равновесных систем.
В равновесных условиях мы, по существу, открываем заново одно из классических распределений вероятности, известное под названием распределения Пуассона. Оно описано в любом учебнике теории вероятностей, поскольку выполняется в огромном числе самых различных случаев: например, по Пуассону, распределены количество вызовов, поступающих на телефонную станцию, время ожидания в ресторане, флуктуации концентрации частиц в жидкости или газе. Математическая формула, задающая распределение Пуассона, для нас сейчас не имеет значения. Мы хотели бы лишь подчеркнуть два аспекта этого важного распределения. Во-первых, оно приводит к закону больших чисел именно в том виде, в каком он сформулирован в предыдущем разделе; следовательно, в большой системе флуктуации допустимо считать пренебрежимо малыми. Во-вторых, закон больших чисел позволяет нам вычислять корреляции между числом молекул Х в двух точках пространства, находящихся на заданном расстоянии друг от друга. Как показывают вычисления, в равновесных условиях такая корреляция не существует. Вероятность одновременно найти молекулу Х в точке r и молекулу X' в точке r' (отличной от точки r) равна произведению вероятности найти молекулу X в точке r и вероятности найти молекулу X' в точке r' (мы рассматриваем случай, когда расстояние между точками r и r' велико по сравнению с радиусом межмолекулярного взаимодействия).
Один из наиболее неожиданных результатов недавних исследований состоял в том, что в неравновесной области ситуация резко изменяется. Во-первых, при подходе вплотную к точкам бифуркации флуктуации стано-
239
вятся аномально сильными и закон больших чисел нарушается. Этого следовало ожидать, так как в сильно неравновесной области система при прохождении точек бифуркации «выбирает» один из различных возможных режимов. Амплитуды флуктуаций имеют такой же порядок величины, как и средние макроскопические значения. Следовательно, различие между флуктуациями и средними значениями стирается. Кроме того, в случае нелинейных химических реакций того типа, который мы рассматривали в гл. 5, появляются дальнодействующие корреляции. Частицы, находящиеся на макроскопических расстояниях друг от друга, перестают быть независимыми. «Отзвуки» локальных событий разносятся по всей системе. Интересно отметить3, что такие дальнодействующне корреляции появляются в самой точке перехода от равновесного состояния к неравновесному. В этом смысле потеря устойчивости равновесным состоянием напоминает фазовый переход, с той лишь особенностью, что амплитуды дальнодействующих корреляций сначала малы, а затем по мере удаления от равновесного состояния нарастают и в точках бифуркаций могут обращаться в бесконечность.
Мы считаем, что такой тип поведения представляет особый интерес, поскольку позволяет подвести «молекулярную основу» под обсуждавшуюся ранее при рассмотрении химических часов проблему связи между частицами. Дальнодействующие корреляции организуют систему еще до того, как происходит макроскопическая бифуркация. Мы снова возвращаемся к одной из главных идей нашей книги: к неравновесности как источнику порядка. В данном случае ситуация особенно ясна. В равновесном состоянии молекулы ведут себя независимо: каждая из них игнорирует остальные. Такие независимые частицы можно было бы назвать гипнонами («сомнамбулами»). Каждая из них может быть сколь угодно сложной, но при этом «не замечать» присутствия остальных молекул. Переход в неравновесное состояние пробуждает гипноны и устанавливает когерентность, совершенно чуждую их поведению в равновесных условиях. Аналогичную картину рисует и микроскопическая теория неравновесных процессов, с которой мы познакомимся в гл.9.
Активность материи связана с неравновесными условиями, порождаемыми самой материей. Так же как и
240
в макроскопическом поведении, законы флуктуации и корреляций в равновесных условиях (когда мы обнаруживаем распределение Пуассона) носят универсальный характер. При переходе границы, отделяющей равновесную область от неравновесной, они утрачивают универсальность и обретают сильнейшую зависимость от типа нелинейности системы.
3. Усиление флуктуаций
Рассмотрим сначала два примера, на которых во всех подробностях можно проследить за ростом флуктуаций, предшествующим образованию новой структуры. Первый пример — образование колонии коллективных амеб, стягивающихся при угрозе голода в единую многоклеточную массу. В гл. 5 мы уже упоминали об этом ярком примере самоорганизации. Другой иллюстрацией роли флуктуации может служить первая стадия постройки гнезда термитами. Она была впервые описана Грассе, а Денюбург исследовал ее с интересующей нас точки зрения4.
Процесс самоорганизации в популяции насекомых
Личинки жука Dendroctonus micans [Scol.] первоначально случайным образом распределены между двумя горизонтальными стеклянными пластинками с зазором 2 мм. С боковых сторон пространство между пластинками открыто. Площадь поверхности 400 см2.
Скопление личинок происходит под влиянием конкуренции двух факторов: случайных движений личинок и их реакции на особое химическое вещество феромон, синтезируемое личинками из терпенов, содержащихся в дереве, которым они питаются. Личинки испускают феромоновые сигналы с частотой, зависящей от степени насыщения. Феромон диффундирует в пространстве, и личинки перемещаются в направлении, задаваемом градиентом его концентрации. Такая реакция является автокаталитическим механизмом, поскольку скопление личинок увеличивает притягательность соответствующей области. Чем выше локальная плотность личинок в данной области, тем выше градиент концентрации феромона и тем сильнее тенденция других личинок к сползанию в точку скопления.
Как показывает эксперимент, плотность популяцни личинок определяет не только скорость, но и эффективность процесса самоорганизации, т.е. число личинок в скоплении на его конечном этапе. При большой плотности (рис. А) скопление возникает и быстро растет в центре экспериментальной установки. При очень малых плотностях устойчивое скопление не образуется (рис. В).
241
Рис. А. Самоорганизация при большой плотности. Распределение личинок через 0 и 21 мин после начала эксперимента.
242
Рис В. Самоорганизация при малой плотности. Распределение личинок через 0 и 22 мин после начала эксперимента.
243
Рис. С. Доля личинок в центральном скоплении (в процентах) от общего числа личинок как функция времени при трех различных плотностях.
Рис. D. Распад начальных ядер из 10 личинок. Общая численность популяции в каждом эксперименте 80 личинок, N — число личинок в ядре.
В других экспериментах исследовалась возможность образования, скопления личинок из «ядра», искусственно созданного на периферии системы. В зависимости от числа личинок в начальном ядре возникают различные ситуации (рис. С).
Если число личинок в ядре мало по сравнению с общим числом личинок, то скопление не образовывалось (рис. D) Если же число. личинок в ядре велико, то скопление растет (рис. Е). При среднем 244
Рис. Е. Рост начальных ядер из 20 (О---О) и 30 (•———•) личинок. Общая численность популяции в каждом эксперименте 80 личинок.
(не слишком большом и не слишком малом) числе личинок в ядре могут возникать структуры новых типов: появляются и сосуществуют два, три или четыре новых скопления с временем жизни не меньшим, чем время наблюдения (рис. F и G).
В экспериментах с однородными начальными условиями такие многокластерные структуры никогда не наблюдались. По-видимому, на бифуркационной диаграмме они соответствуют устойчивым состояниям при допустимых значениях параметров, характеризующих систему, недостижимых из однородных начальных условий. Затравочное ядро выполняет функцию своего рода возмущения, которым необходимо воздействовать на систему для того, чтобы возбудить ее и перевести в область бифуркационной диаграммы, соответствующей" семействам многокластерных распределений.
Постройка гнезда (термитника) термитами — одна из тех когерентных активностей, которые дали некоторым ученым повод для умозрительных утверждений о «коллективном разуме» в сообществах насекомых. Проявляется этот «коллективный разум» довольно необычным способом: для участия в постройке такого огромного и сложного сооружения, как термитник, термитам
245
Рис. F. Многокластерные распределения. В начальном ядре 15 личинок. Общая численность популяции 80 личинок.
необходимо очень мало информации. Первая стадия строительной активности (закладка основания), как показал Грассе, является результатом внешне беспорядочного поведения термитов. На этой стадии они приносят и беспорядочно разбрасывают комочки земли, но каждый комочек пропитывают гормоном, привлекающим других термитов. Ситуацию можно представить следующим образом: начальной «флуктуацией» является несколько большая концентрация комочков земли, которая рано или поздно возникнет в какой-то точке области обитания термитов. Возросшая плотность термитов в окрестности этой точки, привлеченных несколько большей концентрацией гормона, приводит к нарастанию флуктуации. Поскольку число термитов в окрестности точки увеличивается, постольку вероятность сбрасывания ими комочков земли в этой окрестности возрастает, что в свою очередь приводит к увеличению концентрации гормона-аттрактанта. Так воздвигаются «опоры». Расстояние между ними определяется радиусом распро-
246
Рис. G. Рост скопления, искусственно созданного на периферии (верхний рис.), индуцирует образование еще одного небольшого скопления (нижний рис.).
247
странения гормона. Недавно были описаны и другие аналогичные примеры.
Хотя принцип порядка Больцмана позволяет описывать химические или биологические процессы, в которых неоднородности выравниваются, а начальные условия забываются, он не может объяснить ситуации, подобные только что описанным, где несколько «решений», принятых в условиях потери устойчивости, могут направить развитие системы, состоящей из большого числа взаимодействующих единиц, к некоторой глобальной структуре.
Когда новая структура возникает в результате конечного возмущения, флуктуация, приводящая к смене режимов, не может сразу «одолеть» начальное состояние. Она должна сначала установиться в некоторой конечной области и лишь затем распространиться и «заполнить» все пространство. Иначе говоря, существует механизм нуклеации. В зависимости от того, лежат ли размеры начальной области флуктуации ниже или выше критического значения (в случае химических диссипативных структур этот порог зависит, в частности, от кинетических констант и коэффициента диффузии), флуктуация либо затухает, либо распространяется на всю систему. Явления нуклеации хорошо известны из классической теории фазового перехода: например, в газе непрестанно образуются и затем испаряются капельки конденсата. Когда же температура и давление достигают точки, в которой становится устойчивым жидкое состояние, может образоваться капля критических размеров (тем меньших, чем ниже температура и чем выше давление), Если размеры капли превышают порог нуклеации, газ почти мгновенно превращается в жидкость.
Как показывают теоретические исследования и численное моделирование, критические размеры ядра возрастают с эффективностью механизмов диффузии, связывающих между собой все области системы. Иначе говоря, чем быстрее передается сигнал по «каналам связи» внутри системы, тем выше процент безрезультатных флуктуаций и, следовательно, тем устойчивее система. Этот аспект проблемы критического размера означает, что в подобных ситуациях «внешний мир», т. е. все, что окружает флуктуирующую область, всегда стремится погасить флуктуации. Затухнут ли флуктуации или усилятся, зависит от эффективности «канала связи» между
248
флуктуирующей областью и внешним миром. Таким образом, критические размеры определяются конкуренцией между «интегративной силой» системы и химическими механизмами, приводящими к усилению флуктуации.
Описанная нами модель применима, в частности, к результатам, полученным в последнее время in vitro при экспериментальных исследованиях зарождения раковых опухолей5. В этих исследованиях отдельная ракова
Рис. 19. Нуклеация капли жидкости в перенасыщенном паре. а) капля меньше критического размера; b) капля больше критического размера. Существование порога для диссипативных структур подтверждено экспериментально.
клетка рассматривается как флуктуация, способная спонтанно и непрестанно появляться и размножаться посредством репликации. Возникнув, раковая клетка сталкивается с популяцией цитотоксических клеток и либо погибает, либо выживает. В зависимости от значений различных параметров, характеризующих процессы репликации и гибели раковых клеток, мы можем предсказывать либо регресс, либо разрастание опухоли. Такого рода кинетические исследования привели к открытию неожиданных свойств взаимодействия цитотоксических клеток и опухоли: было установлено, что цитотоксические клетки могут принимать мертвые опухолевые клетки за живые. Такие ошибки существенно затрудняют разрушение опухоли.
Вопрос о пределах сложности системы поднимался довольно часто. Действительно, чем сложнее система, тем более многочисленны типы флуктуаций, угрожающих
249
ее устойчивости. Позволительно, однако, спросить, как же в таком случае существуют такие сложные системы, какими является экологическая или социальная структура человеческого общества? Каким образом им удается избежать перманентного хаоса? Частичным ответом на подобные вопросы может быть ссылка на стабилизирующее влияние связи между частями систем, процессов диффузии. В сложных системах, где отдельные виды растений, животных и индивиды вступают между собой в многочисленные и разнообразные взаимодействия, связь между различными частями системы не может не быть достаточно эффективной. Между устойчивостью, обеспечиваемой связью, и неустойчивостью из-за флуктуации имеется конкуренция. От исхода этой конкуренции зависит порог устойчивости.
4. Структурная устойчивость
В каких случаях мы начинаем говорить об эволюции в ее собственном смысле? Как известно, диссипативные структуры требуют сильно неравновесных условий. Тем не менее уравнения реакций с диффузией содержат параметры, допускающие сдвиг в слабо неравновесную область. На бифуркационной диаграмме система может эволюционировать и приближаясь к равновесию, и удаляясь от него, подобно тому как жидкость может переходить от ламинарного течения к турбулентному и возвращаться к ламинарному. Сколько-нибудь жесткой и определенной схемы эволюции не существует.
С совершенно иной ситуацией мы встречаемся в моделях, в которых размеры системы входят в качестве параметра бифуркации: рост, происходящий необратимо во времени, приводит к необратимой эволюции. Однако такой тип развития является достаточно узким частным случаем, хотя вполне возможно, что он имеет некоторое отношение к морфогенетическому развитию.
Ни в биологической, ни в экологической или социальной эволюции мы не можем считать заданным определенное множество взаимодействующих единиц или определенное множество преобразований этих единиц. Это означает, что определение системы необходимо модифицировать в ходе эволюции. Простейший из примеров такого рода эволюции связан с понятием структурной
250
устойчивости. Речь идет о реакции заданной системы на введение новых единиц, способных размножаться и вовлекать во взаимодействие различные процессы, протекающие в системе.
Проблема устойчивости системы относительно изменений такого типа сводится к следующему. Вводимые в небольшом количестве в систему новые составляющие приводят к возникновению новой сети реакций между ее компонентами. Новая сеть реакций начинает конкурировать со старым способом функционирования системы. Если система структурно устойчива относительно вторжения новых единиц, то новый режим функционирования не устанавливается, а сами новые единицы («инноваторы») погибают. Но если структурные флуктуации успешно «приживаются» (например, если новые единицы размножаются достаточно быстро и успевают «захватить» систему до того, как погибнут), то вся система перестраивается на новый режим функционирования: ее активность подчиняется новому «синтаксису»6.
Простейшим примером такого рода может служить популяция макромолекул, образующихся в результате-полимеризации внутри системы, в которую поступают мономеры А и В. Предположим, что процесс полимеризации автокаталитический, т. е. синтезированный полимер используется в качестве образца для образования цепи с той же последовательностью структурных единиц. Такого рода синтез протекает гораздо быстрее, чем синтез в отсутствие образца для копирования. Каждый тип полимеров, отличающийся от других последовательностью расположения в цепи молекул А и В, может быть описан набором параметров, задающих скорость катализируемого синтеза копии, точность процесса копирования и среднее время жизни самой макромолекулы. Можно показать, что при определенных условиях в популяции доминирует полимер какого-то одного типа, например АВАВАВА..., а остальные полимеры могут рассматриваться как «флуктуации» относительно него. Возникающая всякий раз проблема структурной устойчивости обусловлена тем, что в результате «ошибки» при' копировании эталонного образца в системе возникает полимер нового типа, характеризуемый ранее не встречавшейся последовательностью мономеров А и В и новым набором параметров, который начинает размножаться, конкурируя с доминантными видам и за обла-
251
дание мономерами А и В. Перед нами простейший вариант классической дарвиновской идеи о «выживании
наиболее приспособленного».
Аналогичные идеи положены в основу модели предбиотической эволюции, разработанной Эйгеном и его сотрудниками. Подробности теории Эйгена можно найти в многочисленных статьях и книжных публикациях7, поэтому мы ограничимся лишь изложением самой сути. Эйген и его сотрудники показали, что только система
одного типа обладает способностью сопротивляться «ошибкам», постоянно совершаемым автокаталитическими популяциями, — а именно полимерная система, структурно устойчивая относительно появления любого полимера-«мутанта». Такая система состоит из двух множеств полимерных молекул. Молекулы первого множества выполняют функцию «нуклеиновых кислот». Каждая молекула обладает способностью к самовоспроизведению и действует как катализатор при синтезе молекул второго множества, выполняющих функцию
«протеинов». Каждая молекула второго множества катализирует самовоспроизведение молекул первого множества. Такая кросс-каталитическая связь между молекулами двух множеств может превращаться в цикл (каждая «нуклеиновая кислота» воспроизводит себя с помощью «протеина»). Этот цикл обеспечивает устойчивое выживание «нуклеиновых кислот» и «протеинов», защищенных от постоянно возникающих с высоким коэффициентом воспроизводства новых полимеров: ничто не может вмешиваться в самовоспроизводящийся цикл, образуемый «нуклеиновыми кислотами» и «протеинами». Таким образом, эволюция нового типа начинает расти на прочном фундаменте, предвосхищающем появление генетического кода.
Подход, предложенный Эйгеном, несомненно, представляет большой интерес. В среде с ограниченным запасом питательных веществ дарвиновский отбор имеет важное значение для точного самовоспроизведения. Но нам хотелось бы думать, что это не единственный аспект предбиотической эволюции. Не менее важное значение имеют сильно неравновесные условия, связанные с критическими, пороговыми значениями потоков энергии и вещества. По-видимому, разумно предположить, что некоторые из первых стадий эволюции к жизни были связаны с возникновением механизмов, способных погло-
252
щать и трансформировать химическую энергию, как бы выталкивая систему в сильно неравновесные условия. На этой стадии жизнь, или «преджизнь», была редким
событием и дарвиновский отбор не играл такой существенной роли, как на более поздних стадиях.
В нашей книге отношению между микроскопическим и макроскопическим уделяется немало внимания. Одной из наиболее важных проблем в эволюционной теории является возникающая в итоге обратная связь между макроскопическими структурами и микроскопическими событиями: макроскопические структуры, возникая из микроскопических событий, должны были бы в свою очередь приводить к изменениям в микроскопических механизмах. Как ни странно, но в настоящее время наиболее понятные случаи относятся к ситуациям, возникающим в человеческом обществе. Когда мы прокладываем дорогу или строим мост, мы можем предсказать, как это скажется на поведении окрестного населения, а оно в свою очередь определяет изменения в характере и способах связи внутри региона. Такие взаимосвязанные процессы порождают очень сложные ситуации, и это обстоятельство необходимо сознавать, приступая к их моделированию. Именно поэтому мы ограничимся описанием лишь четырех наиболее простых случаев.
5. Логистическая эволюци
Понятие структурной устойчивости находит широкое применение в социальных проблемах. Следует, однако, подчеркнуть, что всякий раз речь идет о сильном упрощении реальной ситуации, описываемой в терминах конкуренции между процессами саморепликации в среде с ограниченными пищевыми ресурсами.
В экологии классическое уравнение, описывающее такую проблему, называется логистическим уравнением. Оно описывает, как эволюционирует популяция из N особей с учетом рождаемости, смертности и количества ресурсов, доступных популяции. Логистическое уравнение можно представить в виде dN/dt=rN(K—N)—mN, где r и m — характерные постоянные рождаемости и смертности, К — «несущая способность» окружающей среды. При любом начальном значении N система со временем выходит на стационарное значение N=K—m/r, завися-
253
Рис. 20. Эволюция популяции N как функция времени t, описываемая логистической кривой. Стационарное состояние N=0 неустойчиво, а стационарное состояние N=K—т/r устойчиво относительно флуктуации величины N.
щее от разности между несущей способностью среды и отношением постоянных смертности и рождаемости. При достижении этого стационарного значения наступает насыщение: в каждый момент времени рождается столько индивидов, сколько их погибает.
Кажущаяся простота логистического уравнения до некоторой степени скрывает сложность механизмов, участвующих в процессе. Мы уже упоминали о внешнем шуме. В случае логистического уравнения он имеет особенно простой смысл. Ясно, что при учете одних лишь климатических флуктуаций коэффициенты К, т и r нельзя считать постоянными: как хорошо известно, такие флуктуации могут разрушить экологическое равновесие и даже обречь популяцию на полное вымирание. Разумеется, в системе начинаются новые процессы, такие, как создание запасов пищи и образование новых колоний, которые заходят в своем развитии настолько далеко, что позволяют в какой-то мере избежать воздействия внешних флуктуации.
Есть в логистической модели и другие тонкости. Вместо того чтобы записывать логистическое уравнение в непрерывном времени, будем сравнивать состояние популяции через заданные интервалы времени (с интервалом, например, в год). Такое дискретное логистическое
254
уравнение представимо в виде Nt+1=Nt(l+r[1—Nt/K]), где Nt и Nt+1 — популяции с интервалом в один год (членом, учитывающим смертность, мы пренебрегаем). Р. Мэй8 обратил внимание на одну замечательную особенность таких уравнений: несмотря на их простоту, они допускают необычайно много решений. При значениях параметра 0?r?2 в дискретном случае так же, как и в непрерывном, наблюдается монотонное приближение к равновесию. При значениях параметра 2
До сих пор мы рассматривали все со статической точки зрения. Обратимся теперь к механизмам, позволяющим варьировать параметры К, r и m в ходе биологической или экологической эволюции.
Следует ожидать, что в процессе эволюции значения экологических параметров К, r и m будут изменяться (так же как и многих других параметров и переменных независимо от того, допускают ли они квантификацию или не допускают). Живые сообщества непрестанно изыскивают новые способы эксплуатации существующих ресурсов или открытия новых (увеличивая тем самым значение параметра К), продления жизни или более быстрого размножения. Каждое экологическое равновесие, определяемое логистическим уравнением, носит лишь временный характер, и логистически заданная экологическая ниша последовательно заполняется серией видов, каждый из которых вытесняет предшествующие, когда его «способность» к использованию ниши, измеряемая величиной К—m/r, становится больше, чем у
255
Рис. 21. Эволюция всей популяции Х как функция времени. Популяция состоит из видов X1, Х2 и Х3, возникающих последовательно и соответствующих возрастающим значениям К—т/r (пояснения см. в тексте).
них (см. рис. 21). Таким образом, логистическое уравнение описывает весьма простую ситуацию, позволяющую количественно сформулировать дарвиновскую идею о выживании «наиболее приспособленного»: наиболее приспособленным считается тот вид, у которого в данный момент времени величина К—т/r больше.
Сколь ни ограниченна задача, описываемая логистическим уравнением, однако и она приводит к некоторым поистине замечательным примерам изобретательности природы.
Возьмем хотя бы гусениц, которые должны оставаться незамеченными, поскольку они движутся слишком медленно, чтобы успеть скрыться от врага.
Выработанные в процессе эволюции стратегии, включающие использование ядов, едких веществ, раздражающих волосков и игл, оказываются высокоэффективными при отпугивании птиц и других потенциальных хищников. Но ни одна из этих стратегий не обладает универсальной эффективностью, способной надежно защитить гусеницу от любого хищника в любое время, в особенности если хищник голоден. Идеальная стратегия состоит в том, чтобы быть как можно более незаметной. Некоторые гусеницы близки к этому идеалу, а при виде разнообразия и изощренности стратегий, используемые сотнями видов чешуекрылых, чтобы остаться незамечен-
256
ными, невольно вспоминаются слова выдающегося натуралиста XIX в. Жан Луи Агассиса: «Экстравагантность настолько глубоко отражает самую возможность существования, что вряд ли найдется какая-нибудь концепция, которую Природа не реализовала бы как слишком экстраординарную»10.
Мы не можем удержаться от искушения привести пример, заимствованный у Милтона Лава11. Трематод (плоский червь), паразитирующий в печени овцы, проходит путь от муравья до овцы, где наконец происходит самовоспроизведение. Вероятность того, что овца проглотит инфицированного муравья, сама по себе очень мала, но поведение такого муравья изменяется самым удивительным образом, и вероятность, по-прежнему оставаясь малой, становится максимальной. Можно с полным основанием сказать, что трематод «завладевает» телом своего хозяина. Он проникает в мозг муравья и вынуждает свою жертву вести себя самоубийственным образом: порабощенный муравей вместо того, чтобы оставаться на земле, взбирается по стеблю растения и, замерев на самом кончике листа, поджидает овцу. Это — поистине «остроумное» решение проблемы для паразита. Остается загадкой, как оно было отобрано.
Модели, аналогичные логистическому уравнению, позволяют исследовать и другие ситуации, возникающие в ходе биологической эволюции. Например, такие модели помогают определить условия межвидовой конкуренции, при которой определенной части популяции выгодно специализироваться на «военной», непроизводительной деятельности (таковы, например, «солдаты» у общественных насекомых). Можно также указать, в какой среде специализированный вид с ограниченным диапазоном пищевых ресурсов имеет более высокую вероятность, выжить, чем неспециализированный вид, потребляющий более разнообразные пищевые ресурсы12. Но здесь мы сталкиваемся с некоторыми весьма различными проблемами организации внутренне дифференцированных популяций. Во избежание путаницы и недоразумений необходимо установить четкие «демаркационные линии». В популяциях, где отдельные особи различимы, где каждая особь наделена памятью, обладает своим характером и опытом и призвана играть свою особую роль, применимость логистического уравнения или, более общо, простого аналога дарвиновских идей
257
становится весьма относительной. В дальнейшем мы еще вернемся к этой проблеме.
Интересно отметить, что кривая на рис. 21, показывающая, как последовательно сменяются при увеличении параметра К—т/r периоды роста и пики семейства решений логистического уравнения, может также описывать размножение некоторых технологических продедур или продуктов. Открытие или технологическое новшество, появление нового продукта нарушает сложившееся социальное, технологическое или экономическое равновесие. Такое равновесие соответствует максимуму кривой роста техники или продуктов производства, с которыми новшеству приходится вступать в конкуренцию (в ситуации, описываемой логистическим уравнением, они играют аналогичную роль13). Приведем лишь один пример. Распространение пароходов привело не только к почти полному исчезновению парусного флота, но и за счет снижения транспортных расходов и повышения скорости перевозок способствовало увеличению спроса на морской транспорт (т. е. увеличению параметра K), что в свою очередь повлекло за собой увеличение численности транспортных судов. Разумеется, ситуация, о которой мы говорим здесь, предельно упрощена и, по предположению, подчиняется чисто экономической логике: технологическое новшество в данном случае лишь удовлетворяет (хотя и иным путем) ранее существовавшую потребность, которая остается неизменной. Но в экологии и человеческом обществе имеется немало примеров инноваций, оказавшихся успешными, несмотря на отсутствие предварительной «ниши».
6. Эволюционная обратная связь
Мы сделаем первый шаг к объяснению эволюционной обратном связи, если будем считать «несущую способность» системы не постоянной, как это было до сих пор, а функцией того, как используется система.
Такое расширение модели позволит нам учесть некоторые дополнительные аспекты экономической деятельности, и в частности некоторые «эффекты усиления». Например, мы получаем возможность описать самоускоряющиеся свойства системы и пространственную дифференциацию различных уровней активности.
258
Географы уже построили модель, коррелирующую эти процессы. Мы имеем в виду модель Кристаллера, определяющую оптимальное пространственное распределение центров экономическом деятельности. Крупные городские центры располагаются в узлах шестиугольной решетки. Каждый из центров окружен кольцом городов следующего по величине масштаба, те в свою очередь окружены тяготеющими к ним еще меньшими населенными пунктами и т. д. Ясно, что в действительности такое геометрически правильное строго иерархическое распределение встречается очень редко: немало исторических, политических и географических факторов нарушают пространственную симметрию. Но рассматриваемая модель нереалистична и по другим причинам. Даже если бы мы исключили все наиболее важные источники асимметричного развития и начали с однородного в экономическом и географическом отношениям пространства, моделирование генезиса распределения, по Кристаллеру, приводит к выводу о том, что описываемая моделью статическая оптимизация является возможным, но маловероятным исходом эволюции.
Рассматриваемая модель14 использует лишь минимальный набор переменных, входящих в вычисления, аналогичные, проведенным Кристаллером. Построено несколько уравнений, обобщающих логистическое. При выводе их авторы исходят из основного предположения о том, что способность населения мигрировать есть функция локальных уровней экономической активности, определяющих своего рода локальную «несущую способность», которая в данном случае сводится к занятости населения. Но местное население есть в то же время потенциальный потребитель товаров, производимых местной промышленностью. Таким образом, в локальном. развитии существует двойная положительная обратная связь, называемая «городским мультипликатором»: и локальное население, и экономическая структура, сложившаяся при уже достигнутом уровне активности, способствуют дальнейшему его повышению. Вместе с тем каждый локальный уровень активности определяется конкуренцией с аналогичными центрами экономической активности, расположенными в других местах. Сбыт произведенных продуктов или оказываемых услуг зависит от стоимости транспортировки их к потребителю и масштабов «предприятия». Расширение любых пред-
259
Рис. 22. Возможная история «урбанизации». Черный кружок — населенные пункты только с одной функцией 1; крупный черный кружок — населенные пункты с функциями 1 и 2; крупный: черный кружок в треугольнике — населенные пункты
А
260
с функциями 1, 2 и 3; крупный черный кружок в треугольнике и квадрате — крупнейшие населенные пункты с функциями 1, 2, 3 и 4. При t=0 (не показано на рисунках) численность населения всех пунктов одинакова и составляет 67 единиц.
В
На рис. С численность населения крупнейшего центра достигает максимума (152 единицы), после чего начинается «расползание» города с образованием городов-спутников.
С
262
Аналогичный процесс наблюдается и в окрестностях главного центра.
D
263
264
приятий определяется спросом на то (товар или услугу), производству чего оно способствует и за производство чего данные предприятия конкурируют с другими. Таким образом, между относительным ростом населения и производительной деятельностью или сферой услуг существуют сильная обратная связь и нелинейные зависимости.
За исходное состояние в рассматриваемой модели приняты гипотетические начальные условия, при которых в различных точках наблюдается (сельскохозяйственная) активность «уровня 1». Модель позволяет проследить возникновение иерархически упорядоченной активности, соответствующей более высоким уровням иерархии по Кристаллеру, т. е. подразумевающей экспорт произведенной продукции в более широкую область. Модель показывает, что даже если начальное состояние совершенно однородно, то одной лишь игры случайных (т. е. не контролируемых моделью) факторов, таких, как место и время закладки различных предприятий, достаточно для нарушения симметрии — появления зон с высокой концентрацией активности и одновременным спадом экономической активности в других областях и оттоком из них населения. Проигрывание модели на ЭВМ позволяет наблюдать расцвет и упадок, подчинение одного экономического центра другому и соответственно доминирование одних центров над другими, периоды, благоприятные для развития альтернативных направлений, и сменяющие их периоды «замораживания» уже существующих структур.
В то время как симметричное распределение Кристаллера игнорирует «историю», изложенный выше сценарий учитывает ее (по крайней мере самым минимальным образом) как взаимодействие «законов», имеющих в данном случае чисто экономическую природу, и «случая», управляющего последовательностью, в которой возникают предприятия.
7. Моделирование сложности
Несмотря на свою простоту, наша модель довольно точно передает некоторые особенности эволюции сложных систем. В частности, она проливает свет на природу трудностей «управления» развитием, зависящим от