C02 — Cumulative costs over a year of activity, as frequency varies

import matplotlib.pyplot as plt
import numpy as np

from cvbook.dati import carica
from cvbook.lingua import t
from cvbook.metriche import equity, rendimenti

CAPITALE = 10_000.0   # ← il capitale impegnato, in euro
                      # PROVA / TRY: il tuo capitale reale — cambia solo la
                      # scala, non le percentuali qui sotto
COSTO_GIRO = 0.0012   # ← costo tutto compreso di un giro completo (0,12%)
                      # PROVA / TRY: 0,0006 (scontato) · 0,0012 · 0,0025 (al dettaglio)
OPERAZIONI_ANNO = 52  # ← quanti giri completi fai in un anno
                      # PROVA / TRY: la tua frequenza reale, contata sull'estratto conto

speso = CAPITALE * COSTO_GIRO * OPERAZIONI_ANNO
pareggio = (1 + COSTO_GIRO) ** OPERAZIONI_ANNO - 1

print(t("in un anno paghi:        ", "in a year you pay:       ") + f"{speso:10,.2f}"
      + t(" euro", " euros"))
print(t("cioe' il:                ", "that is:                 ") + f"{speso / CAPITALE:10.2%}"
      + t(" del capitale", " of capital"))
print(t("rendimento lordo per     ", "gross return to          "))
print(t("NON perdere nulla:       ", "lose NOTHING:            ") + f"{pareggio:10.2%}")

FREQUENZE = np.array([12, 26, 52, 125, 250, 500])
COSTI = [0.0006, 0.0012, 0.0025]

with avvio.figura("schermo"):
    fig, ax = plt.subplots(figsize=(9, 5))
    for c in COSTI:
        soglia = ((1 + c) ** FREQUENZE - 1) * 100
        ax.plot(FREQUENZE, soglia, marker="o", label=t(f"costo {c:.2%} a giro", f"cost {c:.2%} per round trip"))

    # Riferimento in grigio, con l'etichetta in legenda invece che sul grafico:
    # una scritta appoggiata sulle curve e' il modo piu' rapido di rendere
    # illeggibile una figura che dice una cosa semplice.
    # Grey reference, labelled in the legend rather than on the chart: text
    # sitting on top of the curves is the fastest way to make an otherwise
    # simple figure unreadable.
    ax.axhline(
        10,
        linestyle="--",
        linewidth=1.2,
        color="#8C8C8C",
        zorder=0,
        label=t("~10% annuo: media storica di un indice azionario",
                 "~10% a year: historical average of an equity index"),
    )

    ax.set_xscale("log")
    ax.set_xticks(FREQUENZE)
    ax.set_xticklabels([str(f) for f in FREQUENZE])
    ax.set_xlabel(t("Operazioni complete all'anno", "Round trips per year"))
    ax.set_ylabel(t("Rendimento lordo necessario per pareggiare (%)", "Gross return needed to break even (%)"))
    ax.set_ylim(0, None)
    ax.legend(loc="upper left")
    plt.show()

for c in COSTI:
    print(t(f"costo {c:.2%}: ", f"cost {c:.2%}: ") + "  ".join(
        f"{f}op→{((1 + c) ** f - 1):.1%}" for f in FREQUENZE))

SERIE = "btcusdt"  # ← PROVA / TRY: "ethusdt" · "solusdt" (le tre preparate nel setup)
                   # per un'altra delle 11 serie in codice/dati/registro.json
                   # aggiungila anche a avvio.prepara([...]) qui sopra

r = rendimenti(carica(SERIE).sort("data")["chiusura"].to_numpy())
n = len(r)


def con_frequenza(rend: np.ndarray, ogni_n_giorni: int | None, costo: float) -> float:
    """Capitale finale restando sempre investiti ma rientrando ogni N giorni."""
    operazioni = np.zeros(len(rend))
    # L'ingresso del primo giorno si paga sempre, anche a chi non rientrera'
    # mai piu': e' un'operazione. Senza questa riga il compra-e-tieni era
    # l'unica riga della tabella a viaggiare gratis, e il confronto con le
    # altre frequenze partiva gia' truccato di dodici centesimi per mille euro.
    # The first day's entry is always paid, even by whoever never re-enters.
    operazioni[0] = 1.0
    if ogni_n_giorni is not None:
        operazioni[::ogni_n_giorni] = 1.0
    return float(equity(rend - operazioni * costo)[-1])


print(f"{SERIE} · {n}" + t(" giorni\n", " days\n"))
print(f"{t('frequenza', 'frequency'):>22s} " + " ".join(f"{c:>10.2%}" for c in COSTI))
for ogni, etichetta in [
    (None, t("mai (compra e tieni)", "never (buy and hold)")),
    (365, t("una volta l'anno", "once a year")),
    (30, t("una volta al mese", "once a month")),
    (7, t("una volta a settimana", "once a week")),
    (1, t("ogni giorno", "every day")),
]:
    valori = " ".join(f"{con_frequenza(r, ogni, c):10.3f}x" for c in COSTI)
    print(f"{etichetta:>22s} {valori}")