Стохастическая модель

Форум пользователей пакета Mathematica

Модератор: Admin

MaximRi
Сообщения: 2
Зарегистрирован: Ср ноя 14, 2012 4:58 pm

Стохастическая модель

Сообщение MaximRi » Ср ноя 14, 2012 5:05 pm

Уважаемые коллеги помогите пожалуйста написать код, который бы рассчитывал стохастическую модель детерминистической модели см ниже:
(* -> A -> B -> C -> *)
(*C inbibitor of reaction -> A*)

(*Parameters*)
V0 = 0.1; Vmax1 = 0.05; Km1 = 5; Vmax2 = 0.02; Km2 = 6; kd = 0.001;

(*ODE*)
s = NDSolve[{A'[t] == V0*1/(1 + (C[t]/2)^2) - Vmax1*A[t]/(Km1 + A[t]),B'[t] == Vmax1*A[t]/(Km1 + A[t]) - Vmax2*B[t]/(Km2 + B[t]),C'[t] == Vmax2*B[t]/(Km2 + B[t]) - kd*C[t], A[0] == 0, B[0] == 0, C[0] == 0}, {A[t], B[t], C[t]}, {t, 0, 20000}]

(*fig*)
p1 = Plot[A[t] /. s, {t, 0, 20000}, PlotRange -> All,
PlotStyle -> Hue[.0]];(*Red*)
p2 = Plot[B[t] /. s, {t, 0, 20000}, PlotRange -> All,
PlotStyle -> Hue[.3]]; (*Green*)
p3 = Plot[C[t] /. s, {t, 0, 20000}, PlotRange -> All,
PlotStyle -> Hue[.6]];(*Blue*)

Заранее спасибо. Работа будет вознаграждена в денежном эквиваленте.

С уважением, Максим

MaximRi
Сообщения: 2
Зарегистрирован: Ср ноя 14, 2012 4:58 pm

Re: Стохастическая модель

Сообщение MaximRi » Вт дек 04, 2012 12:19 pm

Спасибо. Задача решена
MaximRi писал(а):Уважаемые коллеги помогите пожалуйста написать код, который бы рассчитывал стохастическую модель детерминистической модели см ниже:
(* -> A -> B -> C -> *)
(*C inbibitor of reaction -> A*)

(*Parameters*)
V0 = 0.1; Vmax1 = 0.05; Km1 = 5; Vmax2 = 0.02; Km2 = 6; kd = 0.001;

(*ODE*)
s = NDSolve[{A'[t] == V0*1/(1 + (C[t]/2)^2) - Vmax1*A[t]/(Km1 + A[t]),B'[t] == Vmax1*A[t]/(Km1 + A[t]) - Vmax2*B[t]/(Km2 + B[t]),C'[t] == Vmax2*B[t]/(Km2 + B[t]) - kd*C[t], A[0] == 0, B[0] == 0, C[0] == 0}, {A[t], B[t], C[t]}, {t, 0, 20000}]

(*fig*)
p1 = Plot[A[t] /. s, {t, 0, 20000}, PlotRange -> All,
PlotStyle -> Hue[.0]];(*Red*)
p2 = Plot[B[t] /. s, {t, 0, 20000}, PlotRange -> All,
PlotStyle -> Hue[.3]]; (*Green*)
p3 = Plot[C[t] /. s, {t, 0, 20000}, PlotRange -> All,
PlotStyle -> Hue[.6]];(*Blue*)

Заранее спасибо. Работа будет вознаграждена в денежном эквиваленте.

С уважением, Максим