agreg-activite-bts-ciel/lib/SimuGaz/SimuGaz.cpp

70 lines
2.1 KiB
C++
Raw Normal View History

#include <SimuGaz.h>
#include <Preferences.h>
#define INITIAL_CONSUMPTION 1594.87f
#define INITIAL_SECONDS 11462400 // date initiale = 13/05/2026 16:00
#define AVERAGE_DAILY_CONSUMPTION 2.66f
// Source : opendata.grdf.fr/explore/dataset/correction_climatique_grdf/
// Profil résidentiel chauffage (P011), moyenne 2019-2023, normalisée sur 12 mois
static const float GRDF_MONTHLY_COEFF[12] = {
2.21f, // Janvier
2.05f, // Février
1.62f, // Mars
0.97f, // Avril
0.55f, // Mai
0.22f, // Juin
0.15f, // Juillet
0.17f, // Août
0.38f, // Septembre
0.85f, // Octobre
1.52f, // Novembre
2.11f, // Décembre
};
Preferences prefs;
// Base annuelle : 972 m³/an pour 100m² (CRE)
// - base mensuelle moyenne = 81 m³/mois
// - base journalière moyenne = 2.66 m³/j
float consumptionM3;
uint32_t totalSeconds;
void get_stored_values(){
prefs.begin("simugaz");
consumptionM3 = prefs.getFloat("consumption_m3", INITIAL_CONSUMPTION);
totalSeconds = prefs.getUInt("seconds", INITIAL_SECONDS);
prefs.end();
}
void store_values(){
prefs.begin("simugaz");
prefs.putFloat("consumption_m3", consumptionM3);
prefs.putUInt("seconds", totalSeconds);
prefs.end();
}
uint32_t simulateGasPulses(uint32_t intervalSeconds)
{
get_stored_values();
totalSeconds += intervalSeconds;
uint32_t month = (totalSeconds / (86400UL * 30UL)) % 12;
float dailyBase = AVERAGE_DAILY_CONSUMPTION * GRDF_MONTHLY_COEFF[month];
// Variabilité journalière ±15%
float dailyVariation = 1.0f + 0.30f * ((float)(esp_random() % 1000) / 1000.0f - 0.5f);
// Résolution compteur Gazpar : 0.010 m³ par impulsion
float incrementM3 = dailyBase * intervalSeconds / 86400.0f * dailyVariation;
consumptionM3 += incrementM3;
uint32_t pulseCount = (uint32_t)(consumptionM3 / 0.010f);
uint32_t incrementPulses = (uint32_t)(incrementM3 / 0.010f);
store_values();
Serial.printf("[GAZ] Incrément de %lu impulsions - Total : %lu impulsions (%.3f m³)\n", incrementPulses, pulseCount, pulseCount * 0.010f);
return pulseCount;
}