From 761ac7bcabe308b4d5eb7aa1cb20ba384bafaae4 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Rivny=C3=A1k=20T=C3=ADmea?= <rivnyak.timea@hallgato.ppke.hu> Date: Tue, 1 Jun 2021 21:02:54 +0000 Subject: [PATCH] noisy_deriv for SIRC with noise --- noisy_deriv.m | 13 +++++++++++++ 1 file changed, 13 insertions(+) create mode 100644 noisy_deriv.m diff --git a/noisy_deriv.m b/noisy_deriv.m new file mode 100644 index 0000000..3495fc8 --- /dev/null +++ b/noisy_deriv.m @@ -0,0 +1,13 @@ +function dydt = noisy_deriv(y, nu, mu, epsilon, beta, gamma, Gamma, q) + S = y(1); + I = y(2); + R = y(3); + C = y(4); + noise = 25 * randn(4); + dSdt = nu - (beta*I + epsilon*beta*C)*S - mu*S + noise(1); + dIdt = (beta*I + epsilon*beta*C)*S - gamma*I - mu*I + noise(2); + dRdt = gamma*(1-q)*I + Gamma*C - mu*R + noise(3); + dCdt = gamma*q*I - Gamma*C - mu*C + noise(4); + dydt = [dSdt; dIdt; dRdt; dCdt]; +end + -- GitLab