spb/traite-evaluation-qc

1#!/usr/bin/env python32"""Figure ch05 — sensibilité de la capitalisation d'un écart de taxe foncière : perte de valeur (en % de la valeur) selon3l'écart de taux (points de base de la valeur) pour quatre taux de capitalisation et deux degrés de capitalisation.4Calculs : ch05_rente.capitalisation_taxe (formule exacte) ; aucune donnée observée.5Sortie : assets/figures/ch05/ch05-fig-capitalisation-taxe.pdf"""6import sys, pathlib7ROOT = pathlib.Path(__file__).resolve().parents[3]8sys.path.insert(0, str(ROOT / "analysis/python")); sys.path.insert(0, str(ROOT / "analysis/python/cases"))9import numpy as np10import matplotlib.pyplot as plt11from traite.style import setup, palette, savefig, LARGEUR12import ch05_rente as c51314def main():15    setup()16    pb = np.linspace(0, 100, 101); tau0 = 0.007017    fig, ax = plt.subplots(figsize=(LARGEUR, 3.2))18    for k, i in enumerate([0.04, 0.05, 0.06, 0.07]):19        exact = np.array([c5.capitalisation_taxe(1.0, tau0, tau0 + p / 1e4, i)["dV_pct_exact"] for p in pb]) * 10020        ax.plot(pb, exact, color=palette.cat[k], lw=1.7, label=f"$i$ = {i*100:.0f} %, capitalisation complète (exacte)")21        if i == 0.05:22            ax.plot(pb, 0.6 * exact, color=palette.cat[k], lw=1.2, ls="--", label="$i$ = 5 %, capitalisation à 60 %")23            ax.plot(pb, -pb / 1e4 / i * 100, color=palette.cat[k], lw=0.9, ls=":", label="$i$ = 5 %, approximation $-\\Delta T / i$")24    for x, lab in [(10, "10 pb"), (22.4, "écart nominal\nQuébec–Gatineau 2026"), (50, "50 pb")]:25        ax.axvline(x, color=palette.grid, lw=0.8, ls="--"); ax.text(x + 1, -0.6, lab, fontsize=6.8, color=palette.ink2, va="top")26    ax.set_xlim(0, 100); ax.set_ylim(-25, 0)27    ax.set_xlabel("Écart de taux de taxe foncière (points de base de la valeur marchande)")28    ax.set_ylabel("Variation de la valeur (%)")29    ax.legend(fontsize=6.8, loc="lower left")30    out = savefig(fig, "ch05/ch05-fig-capitalisation-taxe.pdf"); print("OK", out)3132if __name__ == "__main__":33    main()34