Traduction automatique, non relue par un locuteur natif. L'anglais fait foi.

Premiers pasModifier sur GitHubSource : docs/tutorials/03-quantum-vs-classical.md

Tutoriel 3 — Résilience GNSS quantique contre classique (aboutissement)

Kinds: spoof + hybrid · Scénarios : scenarios/spoof-attack.toml, scenarios/hybrid-pnt.toml (copies pédagogiques : scenarios/security.toml, scenarios/hybrid.toml) · Difficulté : intermédiaire · ~35 min

L'aboutissement. Vous allez exécuter un scénario de détecteur de leurrage et une suite complète fusionnée de positionnement, navigation et synchronisation temporelle (PNT), et apprendre à lire ensemble la sécurité, l'intégrité et l'estime — toute la thèse : l'inertiel quantique et la synchronisation temporelle optique procurent une résilience que le refus du GNSS (Global Navigation Satellite System) retirerait autrement.

Partie A — Sécurité : détecter un leurrage temporel#

Un attaquant injecte un faux temps GNSS à rampe lente (0,1 ns/s, à partir de t = 60 s). Le récepteur recoupe le temps affirmé avec la prédiction de marche libre de sa propre horloge et signale le leurrage lorsque l'écart dépasse la borne de détection. Le détecteur est un test d'énergie χ²₁ / de Neyman–Pearson bilatéral à une probabilité de fausse alarme P_fa = 0,01 ; la figure de mérite de sécurité est 1 − P_md (un moins la probabilité de non-détection) pour une amplitude de leurrage correspondant à la spécification (20 ns).

cargo run -- scenarios/spoof-attack.toml

Summary:

scenario 2b6bd22c3b80 | spoof LinearRamp { rate_ns_per_s: 0.1 } vs 20.000 ns spec (P_fa 0.010) | quantum security 1.000 (P_md 0.000, MC 0.000) detected 70s | classical security 0.558 (P_md 0.442, MC 0.444) detected 430s
  • Quantique : security 1.000, détecté à 70 s. Le seuil de détection de l'horloge optique est une fraction de nanoseconde ; elle signale donc le leurrage presque immédiatement — bien en deçà de la spécification de 20 ns.
  • Classique : security 0.558, détecté à 430 s. Le bruit de marche libre propre à l'horloge atomique à puce (CSAC) sur la fenêtre est comparable à 20 ns ; elle ne peut donc souvent pas distinguer le leurrage de sa propre dérive — elle le manque 44 % du temps.

La référence indépendante non circulaire : P_md analytique contre Monte-Carlo#

La probabilité de non-détection du détecteur est calculée de deux façons indépendantes qui doivent concorder : une queue de χ²₁ sous forme fermée (le résultat de Neyman–Pearson) et une estimation Monte-Carlo (MC). Pour la CSAC, le P_md analytique est 0.442 et le P_md MC est 0.444 — ils concordent à 0,002 près, bien à l'intérieur de l'erreur d'échantillonnage ~few × 1/√N. Pour l'horloge optique, les deux valent 0,000. Deux calculs distincts de la même probabilité qui concordent constituent le recoupement non circulaire (Kay, Fundamentals of Statistical Signal Processing: Detection Theory ; détection NP / χ² standard). tests/tutorials.rs::tutorial3_spoof_analytic_matches_montecarlo vérifie |analytic_pmd − mc_pmd| < 0.05 pour les deux horloges ainsi que quantum security_fom > 0.9.

Partie B — PNT hybride : la suite fusionnée complète#

Le scénario d'aboutissement fusionne maintenant une horloge optique et une centrale inertielle (IMU) à atomes froids (suite quantique) face à une CSAC et une IMU de qualité navigation (suite classique), à travers la même panne GNSS de 6600 s (1,83 h), face à une spécification de synchronisation temporelle de 20 ns et à une spécification de position de 100 m. Les deux suites recalent leurs horloges toutes les 60 s sur la même liaison optique de transfert de temps inter-satellites, modélisée avec une gigue de 1 ps par recalage ([resync] dans le scénario ; src/hybrid.rs construit une liaison optical-isl pour l'une ou l'autre suite).

cargo run -- scenarios/hybrid-pnt.toml

Summary:

scenario f33d734ecc51 | quantum PNT-holdover 6600s (t 6600s/p 6600s) integrity 0.998 security n/a (no attack) | classical PNT-holdover 350s (t 6600s/p 350s) integrity 1.000 security n/a (no attack)

Lisez la répartition (t …/p …) — holdover de synchronisation temporelle contre holdover de position :

  • Quantique : PNT-holdover 6600 s — tient à la fois la synchronisation temporelle (t 6600 s) et la position (p 6600 s) pendant toute la panne.
  • Classique : PNT-holdover 350 s — la synchronisation temporelle tient les 6600 s complètes (t 6600 s), parce que le recalage par transfert de temps optique garde verrouillée même l'horloge classique. La suite est limitée en position à p 350 s : l'IMU de qualité navigation est le maillon faible. C'est la thèse de la fusion — isoler la défaillance de la suite classique sur son capteur inertiel.

integrity 0.998 pour la suite quantique est réel et n'est pas un bogue. C'est l'auto-cohérence du filtre sur un recalage bruité — et non une figure d'intégrité aéronautique (niveau de protection horizontal ou vertical, HPL/VPL). Ne surinterprétez pas ici le mot « intégrité » ; voir docs/INTEGRITY.md.

La référence indépendante non circulaire : le holdover de position de 350 s#

Le chiffre de 350 s relève de la physique de l'estime, et non d'un nombre ajusté. Un biais d'accéléromètre constant b fait croître l'erreur de position comme ½·b·T². Franchissement de la spécification de 100 m :

Navigation-grade (classical), b = 1.57e-3 m/s^2:
  T = sqrt(2 * 100 / 1.57e-3) = 357 s

Cold-atom (quantum), b = 5.88e-7 m/s^2 (Templier et al. 2022, arXiv:2209.13209):
  T = sqrt(2 * 100 / 5.88e-7) = 18,443 s >> 6600 s outage -> holds the full outage

Le moteur indique 350 s, le dernier point de la grille de 10 s avant son franchissement et 2 % en deçà des 357 s dus au seul biais : les accéléromètres du scénario portent aussi un bruit de marche aléatoire de vitesse (q_va), qui s'ajoute à la dérive du biais, de sorte que l'erreur modélisée atteint 100 m un peu plus tôt.

Loi de référence : Groves, Principles of GNSS, Inertial, and Multisensor Integrated Navigation (2e éd.), croissance de l'erreur d'estime. Ces deux valeurs de ½bT² rendent compte de la répartition 350 s / 6600 s du moteur : un contrôle physique en forme fermée, externe au filtre.

Vous pouvez voir directement le maillon faible inertiel dans le scénario d'estime (scenarios/inertial.toml, type de scénario inertial) : quantum holdover 6600s p95 41.39m | classical holdover 350s p95 30629.9m : le capteur de qualité navigation diverge jusqu'à des dizaines de kilomètres.

Comment lire ensemble les trois figures#

Chiffre Ce que cela signifie ici Ce que cela ne signifie pas
Sécurité borne analytique de détectabilité du leurrage (1 − P_md) pour une attaque configurée ; le résumé hybrid affiche n/a (no attack) parce que ce scénario n'en configure aucune pas un détecteur RAIM multi-satellites (receiver autonomous integrity monitoring)
Intégrité auto-cohérence du filtre (échantillons à l'intérieur de la borne k-σ) pas une intégrité aéronautique HPL/VPL/RAIM
PNT-holdover temps dans la spécification après la perte du GNSS, réparti entre synchronisation temporelle et position pas une précision CEP (circular error probable) ou 2DRMS (twice the distance root mean square) en 2D

La véritable intégrité autonome du récepteur — de vrais HPL/VPL avec limites d'alerte et un diagramme de Stanford — se trouve dans le paquet integrity (scenarios/integrity.toml) ; c'est l'exercice de lecture de niveau 3.

Ce que les tests fixent#

  • tutorial3_spoof_analytic_matches_montecarlo — le P_md analytique et le P_md MC concordent à moins de 0,05 pour les deux horloges ; quantum security_fom > 0.9.
  • Le résumé hybride (6600 s / 350 s) et la répartition inertielle sont fixés au moyen de leurs scénarios pédagogiques (scenarios/hybrid.toml, scenarios/inertial.toml) par teaching_scenarios_match_their_recorded_expected_summary, qui compare le résumé de chaque moteur à la ligne # expected: que consigne le fichier. La référence indépendante ½bT² ci-dessus est le contrôle physique de ces nombres.

Pièges et unités#

  • La figure de mérite (FoM) de synchronisation temporelle est en nanosecondes ; la FoM de position est en mètres et à 1 DDL (un degré de liberté : un seul axe, ni CEP ni 2DRMS).
  • integrity 0.998 est de l'auto-cohérence, et non de l'intégrité aéronautique.
  • La synchronisation temporelle de la suite classique ne survit que grâce au recalage par transfert de temps optique ; lisez la répartition (t …/p …), pas seulement le PNT-holdover affiché en tête.

Pour aller plus loin#

  • Exécutez le cas d'échec assumé : scenarios/lunanet-araim.toml (type de scénario lunar-integrity) rapporte un passage ARAIM (advanced RAIM) au pôle sud où les niveaux de protection dépassent la limite d'alerte (HPL 263–452 m > 50 m, 0 époque disponible sur 24), si bien que le système est déclaré indisponible : c'est une fonctionnalité, pas un bogue.
  • Portez le récit en orbite terrestre basse (LEO) : scenarios/leo-resilience-gnss-jammed-leo-carries.toml (type de scénario campaign) enchaîne sur une seule chronologie une phase de brouillage du GNSS, une couche de positionnement LEO qui porte l'utilisateur à travers cette phase, et un contrôle d'intégrité. Voir les capacités plus récentes pour les autres scénarios de résilience et verticaux.
  • Défendre un résultat : faites l'exercice de niveau 3 — bandes Monte-Carlo, reproductibilité et lecture d'un niveau de protection.