diff --git a/noisy_deriv.m b/noisy_deriv.m new file mode 100644 index 0000000000000000000000000000000000000000..3495fc868671c5c33c844d87f448681f50e6b15b --- /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 +