В построении математических моделей сложных процессов выделяются следующие этапы.
1. Прежде всего, те реальные явления, которые хотят смоделировать, должны быть тщательно изучены: выявлены главные компоненты и установлены законы, определяющие характер взаимодействия между ними. Если неясно, как связаны между собой реальные объекты, построение адекватной модели невозможно. На этом этапе должны быть сформулированы те вопросы, ответ на которые должна дать модель. Прежде чем строить математическую модель природного явления, надо иметь гипотезу о его течении.
2. Разрабатывается математическая теория, описывающая изучаемые процессы с необходимой детальностью. На ее основе строится модель в виде системы абстрактных взаимодействий. Установленные законы должны быть облечены в точную математическую форму. Конкретные модели могут быть представлены в аналитической форме (системой аналитических уравнений) или в виде логической схемы машинной программы. Модель природного явления есть строгое математическое выражение сформулированной гипотезы.
3. Проверка модели – расчет на основе модели и сличение результатов с действительностью. При этом проверяется правильность сформулированной гипотезы. При значительном расхождении сведений модель отвергают или совершенствуют. При согласованности результатов модели используют для прогноза, вводя в них различные исходные параметры.
Следует, однако, отметить, что сама по себе математическая модель не может служить абсолютным доказательством правильности той или иной гипотезы, так как может оказаться, что разные гипотезы приводят к сходным результатам, но она служит одним из путей анализа реальности.
Расчетные методы в случае правильно построенной модели помогают увидеть то, что трудно или невозможно проверить в эксперименте, позволяют воспроизводить такие процессы, наблюдение которых в природе потребовало бы много сил и больших промежутков времени. В математических моделях можно «проигрывать» разные варианты – вычленять разные связи, комбинировать отдельные факторы, упрощать или усложнять структуру систем, менять последовательность и силу воздействий – все это дает возможность лучше понять механизмы, действующие в природных условиях.
Моделируют различные по характеру процессы, происходящие в реальной среде, как, например, отдельные типы экологических взаимодействий хищник – жертва, паразит – хозяин, конкурентные отношения, мутуализм и др. Математическими моделями описываются и проверяются разные варианты динамики численности, популяций, продукционные процессы в экосистемах, условия стабилизации сообществ, ход восстановления систем при разных формах нарушений и многие другие явления. Сами методы математического моделирования биологических систем развиваются, совершенствуются и разнообразятся.
Например, одну из простейших математических моделей для системы паразит – хозяин в динамике численности насекомых разработал в 1925 г. статистик А. Лотка, который вывел следующие уравнения:
где N1– численность популяции хозяина; N2 – численность популяции паразита; r1– удельная скорость увеличения популяции хозяина; d2 – удельная скорость гибели популяции паразита; p1 и р2– константы. График процесса паразитической инвазии, построенный по таким уравнениям, обнаруживает, что в результате взаимодействия двух видов должны возникать осцилляции (колебания) с постоянной амплитудой, которая зависит от соотношения между скоростями увеличения численности двух видов.
В это же время математик В. Вольтерра выявил сходные закономерности для системы хищник – жертва, обрабатывая статистические данные рыбного промысла. Один из выведенных им законов – «закон периодического цикла» – гласит, что процесс уничтожения одного вида другим может привести к периодическим колебаниям численности популяций обоих видов, зависящих только от коэффициентов роста популяций хищника и жертвы и от исходной относительной численности.
В период, когда были сделаны эти расчеты, экологи вели поиск причин циклических колебаний численности, которые были обнаружены к тому времени у ряда видов. Делались попытки отыскать внешние факторы (космические, солнечные, атмосферные), ответственные за периодические изменения популяций. Модели А. Лотки и В. Вольтерра позволили выдвинуть идею, что периодический колебательный режим в популяциях может возникнуть в результате межвидовых отношений и без внешнего периодического воздействия. Эта идея оказалась плодотворной для дальнейшего развития теории динамики численности популяций. Однако сама модель являлась не адекватной, т. е. не описывала действительность, так как в природе практически не обнаруживаются подобные непрерывные осцилляции с постоянной амплитудой у пар видов, связанных по типу хищник – жертва или паразит – хозяин.
Уравнения А. Лотки и В. Вольтерра были чрезвычайно упрощенными, так как исходили из целого ряда нереальных допущений: что изменение численности популяции одного вида немедленно вызывает ответную реакцию популяции другого вида, что «аппетиты» хищника беспредельны, поиски жертв случайны, что плодовитость хищников пропорциональна численности всей популяции жертв.
Как показал Г. Ф. Гаузе (1934, 1935), даже в условиях упрощенного эксперимента с простейшими трудно добиться соблюдения этих допущений. В его опытах с инфузориями удалось получить лишь два цикла хищник – жертва, после чего система пришла к разрушению. В природе колебания численностей имеют более сложный характер. Во взаимодействиях хищника и жертвы широко распространен эффект «запаздывания» из-за разницы в скоростях размножения, играют роль такие показатели, как степень насыщения («функциональная реакция») хищников, время, затрачиваемое ими на поиск и поимку добычи, способность переключаться на другую пищу, защитные приспособления жертв, размещение их в пространстве и территориальное поведение, возрастная и половая структура популяций и многое другое. Кроме того, рост численности популяций может сдерживаться и другими причинами, в том числе внутривидовыми взаимоотношениями.
В 1933 г. А. Никольсон, несколько усложнив математическую модель Лотки и введя в систему дополнительных хозяев и паразитов, показал, что это ослабляет осцилляции. В 1936 г. А. Н. Колмогоров разработал новые подходы и описал также возможности устойчивого стационарного состояния системы взаимодействующих через трофические связи видов. Позднее для систем хищник – жертва, паразит – хозяин было предложено множество других моделей. С введением в модели дополнительных параметров сильно усложняется математический аппарат и техника расчетов. Многие из этих ограничений позволило снять использование электронно-вычислительных машин.
В экологии сначала преобладали математические модели, основанные на предположениях о существовании в природе четких причинно-следственных зависимостей между популяциями в сообществах (так называемый детерминистский подход). В настоящее время меняется сам подход к математическому моделированию в экологии. Разработаны так называемые имитационные модели, основное внимание в которых уделяется именно разнообразию внутренней структуры популяций и сообществ. Вместо отбрасывания «несущественных» связей математики пытаются определить роль внутреннего разнообразия в поддержании существования надорганизменных систем.
Математическое моделирование широко применяется при решении экологических проблем, связанных с антропогенными воздействиями на природную среду. В современных математических моделях выделяют тактические и стратегические модели. Тактические модели экосистем и популяций служат для экологического прогнозирования их состояния, в том числе при разного рода экзогенных воздействиях. Стратегические модели строят в основном с исследовательскими целями, для вскрытия общих законов функционирования биологических систем, таких, как стабильность, разнообразие, устойчивость к воздействиям, способность возвращаться в исходное состояние. В задачи стратегических моделей входит изучение с помощью ЭВМ последствий разных стратегий управления экосистемами, чтобы иметь возможность выбрать оптимальную.