Исследование системы уравнений Тьюринга

[Yodo]

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


x'=A+x^2*y - (B+1)x
y'=Bx-x^2*y

Установить соотношение параметров A и В, определяющее точку бифуркации. По результатам моделирования построить графики зависимости частоты осцилляций динамических переменных от параметраов А и В
Нашел каим-то чудом коэффициенты ,при которых начинается осцилляция, но дальше вообще не понимаю что надо делать, как найти эту частоту осцилляций и определить точку бифуркации?
[MATLAB]
function dy = turing(t,y,a,b)
dy = zeros(2,1);
dy(1)=a+y(1)*y(1)*y(2)-(b+1)*y(1);
dy(2)=b*y(1)-y(1)*y(1)*y(2);
[/MATLAB]

[MATLAB]
clear all; clc;
t=[0 100];
y0=[1 1];
a=1;
b=2;
[T,Y]=ode45(@turing,t,y0,[],a,b);
plot(T,Y(:,1),'-',T,Y(:,2),'-.')
grid on
title('Численное решение');
legend('X','Y');
xlabel('t,сек');
disp('Вывод значений функций')
[Y(:,1),Y(:,2)][/MATLAB]
вот график осцилляций

[Yodo]

никто не знает?

[s-andriano]

Думаю, этот вопрос лучше бы задавать на форуме химиков.

Если в численном решении наблюдаются осцилляции, можно измерить их перилд и из него определить частоту.