Задача звучит следующим образом:
Численно решить систему уравнений Тьюринга, описывающих динамику изменения концентрации x и y двух веществ в рез-те специфической химической реакции.
x'=A+x^2*y - (B+1)x
y'=Bx-x^2*y
Установить соотношение параметров A и В, определяющее точку бифуркации. По результатам моделирования построить графики зависимости частоты осцилляций динамических переменных от параметров А и В.
Гугл меня ни на что не вывел. Еле-еле нашел два значения параметров,когда начинаются осцилляции. В каких пределах брать A и B ? Должны ли они быть дискретными? Периоды x и y идентичны или отличаются друг от друга.
вот кривой код на матлаб:
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)]
и график http://i076.radikal.ru/1205/90/3ac37d9d447b.jpg
где вообще можно что-нибудь отыскать по этой теме?