lidit

import numpy as np

import pandas as pd

import matplotlib.pyplot as plt

from datetime import datetime, timedelta

 

# 1. Initiële Parameters (29 jan 2026)

start_date = datetime(2026, 1, 29)

idelta_sunset_base = -0.042

rc = np.sqrt(72)

 

# 2. Simulatie van de Maanfase (7-daagse trend)

# We benaderen de maancyclus voor de komende week (richting Volle Maan)

days = np.arange(0, 8)

dates = [start_date + timedelta(days=int(d)) for d in days]

# Maanverlichting stijgt van ~85% naar 100% (Volle Maan rond 1 feb) en neemt daarna af

moon_phases = [0.85, 0.92, 0.98, 1.00, 0.97, 0.92, 0.85, 0.78]

 

# 3. Datalogger Logica: Berekening i-delta per nacht

idelta_log = []

for phase in moon_phases:

# Formule: idelta_Lidit = idelta_sunset + (M_phase * 0.012)

val = idelta_sunset_base + (phase * 0.012)

idelta_log.append(round(val, 4))

 

# 4. Dataframe voor Validatie

df_lidit = pd.DataFrame({

'Datum': [d.strftime('%d-%m-%Y') for d in dates],

'Maanfase (%)': [f"{p*100:.0f}%" for p in moon_phases],

'i-delta Lidit': idelta_log

})

 

# 5. Visualisatie van de trend

plt.style.use('dark_background')

fig, ax1 = plt.subplots(figsize=(10, 5))

 

color = '#00CCFF'

ax1.set_xlabel('Datum')

ax1.set_ylabel('i-delta Stabiliteit', color=color)

ax1.plot(df_lidit['Datum'], df_lidit['i-delta Lidit'], color=color, marker='o', linewidth=2, label='i-delta Trend')

ax1.tick_params(axis='y', labelcolor=color)

ax1.grid(alpha=0.2, linestyle='--')

 

ax2 = ax1.twinx()

color = '#FFCC00'

ax2.set_ylabel('Maanverlichting (%)', color=color)

ax2.bar(df_lidit['Datum'], moon_phases, alpha=0.2, color=color, label='Maanfase')

ax2.tick_params(axis='y', labelcolor=color)

 

plt.title("LIDIT-FASE VOORSPELLING: 7-DAAGSE LUNAIRE DRIFT")

fig.tight_layout()

plt.show()

 

print(df_lidit.to_string(index=False))

Done