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

Работа добавлена на сайт bukvasha.net: 2015-10-28

Поможем написать учебную работу

Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.

Предоплата всего

от 25%

Подписываем

договор

Имя

Выберите тип работы:

Принимаю Политику конфиденциальности

Скидка 25% при заказе до 21.9.2024

СОДЕРЖАНИЕ:
Аннотация................................................................................................................................... 2

Введение..................................................................................................................................... 3

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

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

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

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

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

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

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

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

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

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

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

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

Заключение............................................................................................................................... 29

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

Аннотация

В данной курсовой работе будет проводиться исследование математической модели численности популяции («хищник-жертва») в зависимости от коэффициент хищничества
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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, 8.58856138564497228, .367879347816813306], [1., 44.7734235824848526, .135335242733586658], [1.50000000000000000, ...

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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, .155084361788828830e-1, 4.45671058992297198], [1., .162619865502825509e-2, 1.65800249409337620], [1.5000000000000...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [1., .641972708661369706e-12, 3.77213844452433688], [2., .263123108874330629e-12, .510503200528171885], [3., .740929043476753777e-11, ....

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, .603215359693454228e-8, 21.8621457137488946], [1., -.337113169365853714e-12, 8.04263240339943408], [2., -.5694963...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [1., 43.5493026821559610, .154006065970763939], [2., 1219.60129559745747, 1.32789904118668423], [3., 2.80441607694117501, 1.46320115037...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, 7.91489347278273581, .466627549090192716], [1., 32.3258803157660353, .406246377208466182], [1.50000000000000000, ...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, 5.27390888445715778, 1.00878867815856799], [1., 5.66246113174579868, 1.63650355497690447], [1.50000000000000000, ...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, 3.96850746224258800, 1.38800788137017128], [1., 2.53822515113379943, 1.73763031795569045], [1.50000000000000000, ...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, 2.75053639270340612, 1.83379104028982453], [1., 1.20715844761316404, 1.69168687616629308], [1.50000000000000000, ...

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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, 2.06820835949503224, 2.13779988426323486], [1., .758431481622431348, 1.63720395671045216], [1.50000000000000000, ...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, 1.54727176479670602, 2.40857274912695329], [1., .498879166411486442, 1.58688150947265694], [1.50000000000000000, ...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, 1.15663712180585288, 2.64424219146750650], [1., .338592824265573244, 1.54550548242685726], [1.50000000000000000, ...

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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, .866470907752631914, 2.84733551695322528], [1., .234918132894744920, 1.51373805168522946], [1.50000000000000000, ...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, .651658526387901005, 3.02239398022668615], [1., .165595576421101426, 1.49086785740149020], [1.50000000000000000, ...

Пример в 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]));
res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, .566156566506158420, 3.10098549309420424], [1., .139673095235192502, 1.48243051530274618], [1.50000000000000000, ...

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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, .492515810078810845, 3.17443269341173638], [1., .118106293892748640, 1.47580177094093746], [1.50000000000000000, ...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, .169453862514519270, 3.63578884404674430], [1., .328013105726006408e-1, 1.47298561235965364], [1.5000000000000000...

Пример в 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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, .904162546287280850e-1, 3.86118248818633480], [1., .152175130452786318e-1, 1.50221336948716289], [1.5000000000000...

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]));

res1 := matrix([[vector([t, x(t), y(t)])], [matrix([[0., 3., 1.], [.500000000000000000, .493835355362726853e-1, 4.06679384444323678], [1., .715851428912346426e-2, 1.54570987154564942], [1.5000000000000...

Заключение.

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

при коэффициенте хищничества находящемся в интервале (-∞;0] физического решения нет, хотя и есть математическое решение;

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

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

при коэффициенте хищничества находящемся в интервале (1,8;2,1) количество жертв будет примерно равно количеству хищников;

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

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

при коэффициенте хищничества находящемся в интервале (5; ∞) физического решения нет, хотя и есть математическое решение;

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

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

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

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

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

Bukvasha