Как определить хаотическое поведение в следующей системе ODE? - PullRequest
Новая
       26

Как определить хаотическое поведение в следующей системе ODE?

0 голосов
/ 22 июня 2019

Я пытаюсь выяснить, имеет ли следующая система ODE хаотическую траекторию. Для этого я написал следующий код в MATLAB. 

% x(1) is y and x(2) is k

ODEfcn = @(t,x,a) [0.25*(a*x(1)*(1-x(1)) - 0.0015*x(2) - 0.3*x(1)); a*x(1)*(1-x(1)) - 0.0515*x(2)];
ic = [1 1];
tspan = [0 50];
a = 2;                                                     % Substitute Correct Value
[t,y] = ode45(@(t,x)ODEfcn(t,x,a), tspan, ic);
figure
plot(t, y)
grid

for k = 1:numel(t)
    ddt(k,:) = ODEfcn(t(k),y(k,:),a);                       % Calculate Derivatives From Solved Values & ‘ODEfcn’
end
figure
plot(ddt(:,1), ddt(:,2))
grid
axis equal
xlabel('$\frac{dy(t)}{dt}$', 'Interpreter','latex')
ylabel('$\frac{dk(t)}{dt}$', 'Interpreter','latex')
...