Исследование математической модели численности популяции

Отрицательный коэффициент хищничества. Завышенный коэффициент хищничества. Упрощенный вариант действительности, используемый для изучения ее ключевых свойств. Выводы по влиянию коэффициента хищничества на математическую модель численности популяции.

Рубрика Экономико-математическое моделирование
Вид курсовая работа
Язык русский
Дата добавления 26.12.2014
Размер файла 1,3 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

СОДЕРЖАНИЕ

Аннотация

Введение

1. Описание модели.

2. Исследование модели.

2.1 Имеется математическое решение, но нет физического решения.

2.1.1 Отрицательный коэффициент хищничества (d=-1).

2.1.2 Коэффициент хищничества равный нулю.

2.1.3 Завышенный коэффициент хищничества.

2.2 Имеются как математическое решение так физическое решение.

2.2.1 Количество жертв значительно превышает количество хищников.

2.2.2 Количество жертв не значительно превышает количество хищников.

2.2.3 Количество жертв и количество хищников примерно равны.

2.2.4 Количество хищников не значительно превышает количество жертв.

2.2.5 Количество хищников значительно превышает количество жертв.

Заключение.

Список использованной литературы.

Аннотация

В данной курсовой работе будет проводиться исследование математической модели численности популяции («хищник-жертва») в зависимости от коэффициента хищничества a (скорость, с которой встречи хищников с жертвами позволяют хищникам прибавлять численность своей популяции).

Введение

Актуальность темы.

Впервые математическая модель «хищник-жертва» была получена А. Лоткой (1925 г.), который использовал для описания динамики взаимодействующих биологических популяций. Чуть позже и независимо от Лотки аналогичные (и более сложные) модели были разработаны итальянским математиком В. Вольтерра (1926 г.), глубокие исследования которого в области экологических проблем заложили фундамент математической теории биологических сообществ или так называемой математической экологии.

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

Цель работы.

Целью курсовой работы является закрепление теоретических знаний, полученных по данной и смежным дисциплинам, освоение современных методов и средств моделирования систем. А также выявление зависимости математической модели «хищник-жертва» от коэффициента хищничества.

Основные методы исследования.

В данной работе в качестве метода исследования используется исследование с помощью математического моделирования.

Математическая модель - это упрощенный вариант действительности, используемый для изучения ее ключевых свойств. Чарльз Лейв и Джеймс Марч дают такое определение модели: “Модель - это упрощенная картина реального мира. Она обладает некоторыми, но не всеми свойствами реального мира. Она представляет собой множество взаимосвязанных предположений о мире. Как и любая картина, модель проще тех явлений, которые она по замыслу отображает или объясняет”.

Вклад автора курсовой работы в рассматриваемую проблему.

Моим вкладом в рассматриваемую проблему является исследование влияния коэффициента хищничества на математическую модель «хищник-жертва».

Практическая (теоретическая) ценность курсовой работы.

Ценность данной курсовой работы заключается в выявленной зависимости влияния коэффициента хищничества на математическую модель «хищник-жертва».

Реализация и внедрение курсовой работы.

Данную курсовую работу можно использовать при моделировании долговременных отношений между видами хищника и жертвы в какой-либо экосистеме. К примеру: караси и щуки, зайцы и волки, мыши и лисы, микробы и антитела и т. д.

1. Описание модели

Иногда простая математическая модель хорошо описывает сложную биологическую систему. Примером этого служат долговременные отношения между видами хищника и жертвы в какой-либо экосистеме. Математические расчеты роста популяции отдельно взятого вида показывают, что пределы плотности популяции можно описать простыми уравнениями. Кривая численности популяции растет экспоненциально, пока она небольшая, а затем выравнивается, когда она достигает пределов возможности экосистемы поддерживать ее. Простое продолжение этой концепции позволяет нам понять экосистему, в которой взаимодействуют два вида -- хищник и жертва.

Итак, если число растительноядных жертв x, а число плотоядных хищников y, то вероятность, что хищник встретится с травоядным, пропорциональна произведению xy. Другими словами, чем выше численность одного из видов, тем выше вероятность таких встреч. В отсутствие хищников популяция жертвы будет расти экспоненциально (по крайне мере вначале), а в отсутствие жертв популяция хищника сократится до нуля -- либо из-за голода, либо в результате миграции. Теперь, если dx -- изменение популяции растительноядных за время dt, а dy изменение популяции плотоядных за тот же интервал времени, то две популяции описываются уравнениями:

Здесь a -- скорость роста численности травоядных в отсутствие хищников, а c -- скорость сокращения численности плотоядных в отсутствие травоядных. Постоянные b и d -- скорость, с которой встречи хищников с жертвами удаляют травоядных из популяции, и скорость, с которой эти встречи позволяют хищникам прибавлять численность своей популяции. Знак минус в первом уравнении показывает, что встречи сокращают популяцию жертвы, а знак плюс во втором говорит о том, что встречи увеличивают популяцию хищника. Как видите, любое изменение численности травоядных влияет на численность плотоядных, и наоборот. Две популяции необходимо рассматривать вместе.

Решение этих уравнений показывает, что обе популяции развиваются циклически. Если популяция травоядных увеличивается, вероятность встреч хищник--жертва возрастает, и, соответственно (после некоторой временной задержки), растет популяция хищников. Но рост популяции хищников приводит к сокращению популяции травоядных (также после некоторой задержки), что ведет к снижению численности потомства хищников, а это повышает число травоядных и так далее.

2. Исследование модели

Разделим исследование модели на следующие этапы:

§ Имеется математическое решение, но нет физического решения;

o при отрицательном коэффициенте хищничества;

o при нулевом коэффициенте хищничества;

o при завышенном коэффициенте хищничества.

§ Имеются как математическое решение так физическое решение;

o Количество жертв значительно превышает количество хищников;

o Количество жертв не значительно превышает количество хищников;

o Количество жертв и количество хищников примерно равны;

o Количество хищников не значительно превышает количество жертв;

o Количество хищников значительно превышает количество жертв;

2.1 Имеется математическое решение, но нет физического решения

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

· при отрицательном коэффициенте хищничества;

· при нулевом коэффициенте хищничества;

· при завышенном коэффициенте хищничества.

2.1.1 Отрицательный коэффициент хищничества (d=-1)

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

Пример в Maple 9 при d=-1:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2-1*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..2.5);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..1);

> with(plots):odeplot(V,[t,y(t)],0..0.8);

> odeplot(V,[x(t),y(t)],0..2.6);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2-1*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=

array([0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.25,2.5]));

2.1.2 Коэффициент хищничества равный нулю

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

Пример в Maple 9 при d=0:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+0*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],7..10);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..1);

> with(plots):odeplot(V,[t,y(t)],0..0.8);

> odeplot(V,[x(t),y(t)],0..2.6);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+0*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=

array([0,0.5,1,1.5,2,2.5,3,3.5,4,4.5,5]));

2.1.3 Завышенный коэффициент хищничества

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

Пример в Maple 9 при d=5:

Здесь я принял решение, что не будет физического решения в связи с тем, что количество жертв доходит до очень низкого значения близкого к нулю (у>0,002). Получается что в данной модели количество жертв очень мало и скорее всего с чем меньше количества жертв те меньше вероятность увеличения их популяции.

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+5*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> with(plots):odeplot(V,[t,x(t)],0.5..2);

> odeplot(V,[x(t),y(t)],0..3.6);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+5*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=

array([0,0.5,1,1.5,2,2.5,3,3.5,4,4.5,5]));

Пример в Maple 9 при d=19.76:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+19.76*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..21);

> odeplot(V,[x(t),y(t)],0..8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+19.76*x(t))*y(t),x(0)=3,y(0)=1};

>res1:=dsolve(eq1,type=numeric,output=array([0,1,2,3,4,5,6,7,8,9,10,11,12]));

Пример в Maple 9 при d=50:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+50*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],7..10);

> odeplot(V,[x(t),y(t)],0..10);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+50*x(t))*y(t),x(0)=3,y(0)=1};

>

res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,2,3,4,5,6,7,8,9,10,11,12]));

Пример в Maple 9 при d=100:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+100*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..50);

> odeplot(V,[x(t),y(t)],0..10);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+100*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,2,3,4,5,33,34]));

2.2 Имеются как математическое решение так физическое решение

Т.е. складывается такая ситуация когда имеются как математическое решение, так и физическое решению, иными словами при таких параметрах математической модели возможно физическое решение. Произведём следующее разделение:

· количество жертв значительно превышает количество хищников;

· количество жертв не значительно превышает количество хищников;

· количество жертв и хищников примерно равны;

· количество хищников не значительно превышает количество жертв;

· количество хищников значительно превышает количество жертв;

2.2.1 Количество жертв значительно превышает количество хищников

В этом случае количество хищников очень мало, а количество жертв очень велико это связано с тем, что коэффициент хищничества, т.е. скорость, с которой встречи хищников с жертвами позволяют хищникам прибавлять численность своей популяции, очень мал (d в интервале от 0,01 до 1).

Пример в Maple 9 при d=0,01:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+0.01*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> with(plots):odeplot(V,[t,y(t)],0..5);

> odeplot(V,[x(t),y(t)],0..4);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+0.01*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,1,2,3,4,5,6,7,8]));

Пример в Maple 9 при d=0,1:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+0.1*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+0.1*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

Пример в Maple 9 при d=0,5:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+0.5*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+0.5*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

Пример в Maple 9 при d=0.73(минимум численности жертв совпадает с максимумом численности хищников):

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+0.73*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..5);

> odeplot(V,[x(t),y(t)],0..4);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+0.73*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

Пример в Maple 9 при d=1:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+1*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+1*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

2.2.2 Количество жертв не значительно превышает количество хищников

В этом случае количество хищников мало, а количество жертв велико это связано с тем, что коэффициент хищничества, т.е. скорость, с которой встречи хищников с жертвами позволяют хищникам прибавлять численность своей популяции, мал (d в интервале от 1,2 до 1,6).

Пример в Maple 9 при d=1,2:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+1.2*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+1.2*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

Пример в Maple 9 при d=1,4:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+1.4*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+1.4*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

Пример в Maple 9 при d=1,6:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+1.6*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+1.6*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

2.2.3 Количество жертв и количество хищников примерно равны

В этом случае количество хищников и количество жертв примерно равны, это связано с тем, что коэффициент хищничества, т.е. скорость, с которой встречи хищников с жертвами позволяют хищникам прибавлять численность своей популяции имеет среднее значение (d в интервале от 1,8 до 2,1).

Пример в Maple 9 при d=1,8:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+1.8*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+1.8*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

Пример в Maple 9 при d=2:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+2*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+2*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

Пример в Maple 9 при d=2,1:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+2.1*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+2.1*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

2.2.4 Количество хищников не значительно превышает количество жертв

В этом случае количество жертв мало, а количество хищников велико это связано с тем, что коэффициент хищничества, т.е. скорость, с которой встречи хищников с жертвами позволяют хищникам прибавлять численность своей популяции, высок (d в интервале от 2,2 до 3,5).

Пример в Maple 9 при d=2,2:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+2.2*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.8);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+2.2*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

Пример в Maple 9 при d=3:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+3*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..2.9);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+3*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

Пример в Maple 9 при d=3,5:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+3.5*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> odeplot(V,[x(t),y(t)],0..3);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+3.5*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

2.2.5 Количество хищников значительно превышает количество жертв

В этом случае количество жертв мало, а количество хищников велико это связано с тем, что коэффициент хищничества, т.е. скорость, с которой встречи хищников с жертвами позволяют хищникам прибавлять численность своей популяции, высок (d около 4).

Пример в Maple 9 при d=4:

> V:=dsolve({diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+4*x(t))*y(t),x(0)=3,y(0)=1},{x(t),y(t)},numeric);

> with(plots):odeplot(V,[[t,x(t)],[t,y(t)]],0..10);

> with(plots):odeplot(V,[t,x(t)],0.5..2);

> odeplot(V,[x(t),y(t)],0..3.6);

> eq1:={diff(x(t),t)=(4-3*y(t))*x(t),diff(y(t),t)=(-2+4*x(t))*y(t),x(0)=3,y(0)=1};

> res1:=dsolve(eq1,type=numeric,output=array([0,0.5,1,1.5,2,2.5,3,3.5,4]));

Заключение

коэффициент хищничество математический модель

В заключении можно сказать что, исследуя математическую модель численности популяции («хищник-жертва») в зависимости от коэффициент хищничества a (скорость, с которой встречи хищников с жертвами позволяют хищникам прибавлять численность своей популяции) я сделал следующие выводы по влиянию коэффициента хищничества на данную модель:

при коэффициенте хищничества находящемся в интервале (-?;0Размещено на http://www.allbest.ru/

] физического решения нет, хотя и есть математическое решение;

при коэффициенте хищничества находящемся в интервале (0,01;1) колРазмещено на http://www.allbest.ru/

ичество жертв будет значительно превышать количества хищников;

при коэффициенте хищничества находящемся в интервале (1,2;1,6) кРазмещено на http://www.allbest.ru/

оличество жертв будет не значительно превышать количества хищников;

при коэффициенте хищничества находящемся в интервале (1,8;2,1) кРазмещено на http://www.allbest.ru/

оличество жертв будет примерно равно количеству хищников;

при коэффициенте хищничества находящемся в интервале (2,2;3,5) кРазмещено на http://www.allbest.ru/

оличество хищников будет не значительно превышать количества жертв;

при коэффициенте хищничества имеющем значение близкое 4Размещено на http://www.allbest.ru/

количество хищников будет значительно превышать количества жертв;

при коэффициенте хищничества находящемся в интервале (5;Размещено на http://www.allbest.ru/

?) физического решения нет, хотя и есть математическое решение;

с увеличением коэффициента хищничества количество хищников увелиРазмещено на http://www.allbest.ru/

чивается, а количество жертв уменьшается и наоборот, с уменьшением коэффициента хищничества количество хищников уменьшается, а количество жертв увеличивается, т.е. количество хищников прямо пропорционально коэффициенту хищничества, а количество жертв обратно пропорционально;

Как видно в перечисленных выводах нет чёткой границы в интервалах, это связанно с тем, что невозможно поставить чёткую границу когда количество одних превышает количества других.

В качестве хищника и жертвы могут быть взяты караси и щуки, зайцы и волки, мыши и лисы, микробы и антитела, рыба и планктон и т. д. Другими словами в качестве хищников взяты те, кто питается только жертвами, тем самым, увеличивая свою популяцию. А в качестве жертв те, кем питаются хищники, причем в отсутствии хищников жертвы должны увеличивать свою популяцию.

Список использованной литературы

1) Ризниченко Г. Ю., Рубин А. Б. Математические модели биологических продукционных процессов, Изд. МГУ, Москва, 1993.

2) Семакин И.Г., Хеннер Е.К. Информатика. Задачник-практикум в 2 т.: ЛБЗ, 1999 г., с. 188-189.

3) Вольтерра В. Математическая теория борьбы за существование. Наука, Москва, 1976.

4) http://www.math.rsu.ru/mexmat/kvm/MME/dsarch/hishnik.html

5) http://elementy.ru/trefil/83

6) http://www.dmb.biophys.msu.ru/registry?article=9821

Размещено на Allbest.ru

...

Подобные документы

  • Множественная линейная регрессия: спецификация модели, оценка параметров. Отбор факторов на основе качественного теоретико-экономического анализа. Коэффициент регрессии при фиктивной переменной. Проблемы верификации модели. Коэффициент детерминации.

    контрольная работа [88,0 K], добавлен 08.09.2014

  • Построение математической модели выбранного экономического явления методами регрессионного анализа. Линейная регрессионная модель. Выборочный коэффициент корреляции. Метод наименьших квадратов для модели множественной регрессии, статистические гипотезы.

    курсовая работа [1,1 M], добавлен 22.05.2015

  • Построение эконометрической модели. Описания, анализ и прогнозирование явлений и процессов в экономике. Использование регрессионных моделей. Построение корреляционной матрицы. Коэффициент множественной детерминации. Значение статистики Дарбина-Уотсона.

    курсовая работа [61,0 K], добавлен 10.03.2013

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

    курсовая работа [298,3 K], добавлен 20.03.2011

  • Модель взаимодействия хищников и жертв на плоскости. Построение модели, которая даёт визуальное представление о происходящем взаимодействии хищников и жертв. Стабилизация численности жертв к равновесной численности в условиях отсутствия хищников.

    реферат [538,9 K], добавлен 02.03.2011

  • Определение временных и пространственных данных в эконометрике. Коэффициент детерминации и средняя ошибка аппроксимации как показатели качества однофакторной модели в эконометрике. Особенности построения множественной регрессивной модели. Временные ряды.

    контрольная работа [804,3 K], добавлен 15.11.2012

  • Построение модели парной регрессии и расчет индекса парной корреляции. Построение производственной функции Кобба-Дугласа, коэффициент детерминации . Зависимость среднедушевого потребления от размера дохода и цен. Расчет параметров структурной модели.

    контрольная работа [1,6 M], добавлен 05.01.2012

  • Оценка адекватности эконометрических моделей статистическим данным. Построение доверительных зон регрессий спроса и предложения. Вычисление коэффициента регрессии. Построение производственной мультипликативной регрессии, оценка ее главных параметров.

    контрольная работа [1,2 M], добавлен 25.04.2010

  • Роль экономико-математических методов в оптимизации экономических решений. Этапы построения математической модели и решение общей задачи симплекс-методом. Составление экономико-математической модели предприятия по производству хлебобулочных изделий.

    курсовая работа [1,3 M], добавлен 09.07.2015

  • Линеаризация математической модели регулирования. Исследование динамических характеристик объекта управления по математической модели. Исследование устойчивости замкнутой системы управления линейной системы. Определение устойчивости системы управления.

    курсовая работа [1,6 M], добавлен 07.08.2013

  • Группировка предприятий по среднегодовой стоимости производственных фондов. Сглаживание скользящей средней и ее центрирование. Определение коэффициента линейной регрессионной модели и показателей детерминации. Коэффициенты эластичности и их интерпретация.

    контрольная работа [493,4 K], добавлен 06.05.2015

  • Сущность экономико-математической модели, ее идентификация и определение достаточной структуры для моделирования. Построение уравнения регрессии. Синтез и построение модели с учетом ее особенностей и математической спецификации. Верификация модели.

    контрольная работа [73,9 K], добавлен 23.01.2009

  • Основные понятия математической статистики. Нахождение коэффициента эластичности модели. Проведение экономического анализа, составление прогноза и построение доверительной области. Вычисление зависимости показателя от фактора. Проверка созданной модели.

    контрольная работа [173,9 K], добавлен 19.06.2009

  • Модели распределения доходов. Количественный подход к анализу полезности и спроса. Кривые безразличия, решение задачи об оптимальном выборе потребителя. Функции спроса и коэффициент эластичности. Предельная полезность и предельная норма замещения.

    презентация [470,8 K], добавлен 28.04.2013

  • Модели распределения доходов. Количественный подход к анализу полезности и спроса. Отношение предпочтения и функция полезности. Кривые безразличия, решение задачи оптимального выбора потребителя. Функции спроса, изменение цен и коэффициент эластичности.

    курсовая работа [412,7 K], добавлен 11.02.2011

  • Разработка оптимального режима процесса получения максимального выхода химического вещества. Получение математической модели процесса с применением метода центральных композиционных ортогональных планов второго порядка. Исследование поверхности отклика.

    курсовая работа [104,3 K], добавлен 20.07.2012

  • Оценить влияние определенных факторов на изучаемый показатель и друг на друга с помощью коэффициентов линейной корреляции. Среднее квадратическое отклонение фактора. Коэффициент линейной корреляции. Линейные регрессионные модели изучаемого показателя.

    контрольная работа [381,3 K], добавлен 21.04.2010

  • Расчет параметров уравнения регрессии, среднего коэффициента эластичности и средней ошибки аппроксимации по рынку вторичного жилья. Определение идентификации моделей денежного и товарного рынков, выбор метода оценки параметров модели, оценка его качества.

    контрольная работа [133,1 K], добавлен 23.06.2010

  • Понятие и особенности прогнозирования. Стандартная ошибка предсказываемого среднего значения. Прогнозирование при наличии авторегрессии ошибок. Точечное и интервальное прогнозирование, основанное на модели линейной регрессии, коэффициент ее детерминации.

    контрольная работа [827,9 K], добавлен 08.01.2016

  • Взаимосвязи экономических переменных. Понятие эконометрической модели. Коэффициент корреляции и его свойства. Линейная парная регрессия. Метод наименьших квадратов. Основные предпосылки и принципы регрессионного анализа. Статистика Дарбина-Уотсона.

    шпаргалка [142,4 K], добавлен 22.12.2011

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.