diff --git a/lotkavolterra_antibiotics.m b/lotkavolterra_antibiotics.m index 7f1583df261e3b3c1d3527a2227b1ef063cea84d..a71829e43890af69f8d9987e5564f06b6078dec7 100644 --- a/lotkavolterra_antibiotics.m +++ b/lotkavolterra_antibiotics.m @@ -78,8 +78,13 @@ function dxdt = lotka_volterra_with_antibiotics(t, x, r, K, A, antibiotics) % Növekedési ráták csökkentése az időfüggő antibiotikum hatással r_mod = r_mod .* (1 - antibiotic_effect * sensitivity_vector'); - resistance_factor = 1 + 0.01 * t; % Rezisztencia fokozatos növekedése + + scaled_time = t - 30; % Átskálázott idő: antibiotikum hatás kezdete legyen a 0 + resistance_growth_rate = 0.01; + resistance_factor = 1 + resistance_growth_rate * scaled_time; r_mod = r_mod .* (1 ./ resistance_factor); + max_increase = 1.5; % Maximális eltartóképesség növekedési faktor (pl. 1.5 = 50%-os növekedés) + K_mod = K_mod .* (1 + (max_increase - 1) * (1 - exp(-resistance_growth_rate * scaled_time))); % Az antibiotikumok antagonisztikus hatásának figyelembevétele interaction_modifier = ones(size(r_mod)); % Kezdetben nincs hatás