Опис файлу

clear all;

close all;

global s;

global K;

global r;

global k;

global D;

global h;

K=500;

r=0.895;

s=0.7;

D=3;

k=1;

h=1;

function y=lotkifun4(x);

global s;

global K;

global r;

global k;

global D;

global h;

y(1)=x(1).*(r.*(1-(x(1)/K))-((k.*x(2))/(x(1)+D)));

y(2)=x(2).*(s.*(1-(h.*x(2))/x(1)));

endfunction

t=[1:0.5:39];

x0=[10,10];

Z=lsode("lotkifun4", x0,t);

[X,Y]=meshgrid (5:5.5:45, 5:5.5:45); %tworzenie siatki

F=X.*(r.*(1-X./K)-k.*Y./(X+D)); %tworzenie wektrów

G=Y.*(s.*(1-h.*Y./X));

n=sqrt(F.^2+G.^2); %obliczanie długoœci wektrów

g=G./n;

f=F./n; %normalizacja wektrów

for i=1:20

j=i*4;

subplot(2,1,1) %tworzenie górnego rysunku

plot(t,Z(:,2),t,Z(:,1),'r') %tworzenie wykresu drapieżników i ofiar

hold on

legend('drapieżnicy','ofiary')

legend('location','north')

plot(t(j),Z(j,1),'o',"markersize",10) %tworzenie punktów

plot(t(j),Z(j,2),'o',"markersize",10)

xlabel('czas')

ylabel('liczebnoœć populacji')

subplot(2,1,2) %tworzeni dolnego rysunku

plot(Z(:,2),Z(:,1),'r') %rysowanie portretu fazowego

hold on

h=quiver(X,Y,f,g,1); %tworzenie pola wektorowego '1' oznacza długoœć wektoru

set (h,"maxheadsize",0.35) %ustalanie wielkoœci strzałki

plot(Z(j,2),Z(j,1),'o',"markersize",10) %rysowanie punktów

xlabel('populacja drapieżników')

ylabel('populacja ofiar')

hold off

nazwa= sprintf ( "NAZWA%04d.png",i)

print(nazwa,'-dpng','-S500,620')

close all

endfor

Ліцензування

Я, власник авторських прав на цей твір, добровільно публікую його на умовах такої ліцензії:

Цей файл доступний на умовах Creative Commons CC0 1.0 Universal Public Domain Dedication.

Особа, що пов'язала роботу з даною дією, передала роботу у суспільне надбання шляхом відмови від усіх своїх прав на роботу по всьому світу по закону про авторське право, включаючи всі пов'язані і суміжні права, в тій мірі, що допускається законом. Ви можете копіювати, змінювати, розповсюджувати і виконувати роботу, навіть на комерційній основі, не питаючи дозволу. http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse