Fusjonert MEO- og LEO-posisjonering, navigasjon og tidssynkronisering
Tre scenariotyper gir posisjonering, navigasjon og tidssynkronisering (PNT) med satellitter i lav jordbane (LEO) alene eller fusjonert med de globale satellittnavigasjonssystemene (GNSS) i middels høy jordbane (MEO): leo-pvt (posisjon, hastighet og tid), leo-ppp (presis punktposisjonering) og ntn-positioning (et 5G ikke-terrestrisk nettverk). Kilde: src/leo_fusion/.
Motoren er systemuavhengig. Hver konstellasjon er Walker-skall, eksplisitte baneelementer eller en GNSS-forhåndsinnstilling i typen constellation-design; hvert signal er en bærebølge, en brikketakt og en omhyllingskurve for bærebølge-til-støy-tetthet (C/N0); hvert feilbudsjett er et eksplisitt énsigma. Navngitte systemer er valgfrie forhåndsinnstillinger, én fil hver, og hver funksjon, test og medfølgende scenario kjører uten noen bestemt av dem.
Systemer og forhåndsinnstillinger#
Et scenario lister systemene sine som [[system]]-tabeller. Dette fullstendige scenarioet slår sammen GPS-grunnlinjekonstellasjonen med et generisk LEO-lag og kjører slik det står (kshana <file>.toml; motor 0.40.0 gir en median posisjonsfortynning av presisjon på 2,43 for bare GPS og 1,17 fusjonert, og en rotmiddelkvadrat 3D-feil på 2,06 m og 0,47 m):
kind = "leo-pvt"
mode = "joint"
seed = 3
duration_s = 600.0
step_s = 60.0
[user]
lat_deg = 40.42
lon_deg = -3.70
height_m = 650.0
[[system]]
name = "GPS"
preset = "gps-baseline"
cn0_dbhz = [38.0, 48.0]
sisre_m = 0.6
[[system]]
name = "Any LEO layer"
carrier_hz = 1.5e9 # any band
chip_rate_hz = 10.23e6 # omit with doppler_only = true for a signal of opportunity
cn0_dbhz = [45.0, 55.0] # at the elevation mask and at the zenith
sisre_m = 0.3 # signal-in-space range error of orbit and clock, one sigma
clock = "estimated" # its own receiver clock (inter-system bias solved) or "known"
isb_ns = 40.0 # true inter-system bias used to simulate its pseudoranges
[[system.shell]]
total = 240
planes = 12
phasing = 1
altitude_km = 1000.0
inclination_deg = 60.0
preset = "gps-baseline" (Global Positioning System, GPS; eller galileo, beidou-meo, glonass) tar en publisert GNSS-konstellasjon; leo_preset = "..." fyller ut det tabellen utelater fra en LEO-forhåndsinnstilling. C/N0 løper lineært i sinus til elevasjonen fra masken til senit (et modelleringsvalg: en virkelig passeringsomhyllingskurve avhenger av sendeantennens mønster). Pseudoavstandens énsigma er motorens termiske støy i delay-lock-loop (DLL) (navsignal::dll_code_jitter_chips) ved den C/N0-en kombinert med avstandsfeilen i signal i rommet (SISRE), med mindre sigma_pr_m fastsetter den; avstandsratens énsigma er bølgelengden ganger sigma_doppler_hz. Fordi sigma er en inndata, kan en egen signalmodell mate den uten at posisjoneringskoden vet hvor den kom fra.
| Forhåndsinnstilling | Hva det omfatter | Kilder |
|---|---|---|
xona-pulsar |
X1 ved 1593,3225 MHz, 1,023 Mchip/s, mottatt effekt −148,2 til −139,1 dBW; X5 ved 1190,51625 MHz, 10,23 Mchip/s, −144,9 til −136,2 dBW; 258 satellitter i 18 baneplan på 1080 km, 53 og 97 grader | PUBLIC: arXiv:2509,19551. Representativt: fordelingen mellom de to banehellingene, faseforskyvningen, avstandsfeilen |
iridium-stl |
Iridium Satellite Time and Location (STL): 66 satellitter, seks nesten polare baneplan på 780 km og 86,4 grader; bursts ved 1616 til 1626 MHz, kvadraturfaseskiftmodulasjon (QPSK) ved 25 000 symboler/s; mottatt effekt 300 til 2400 ganger GPS-minimum | PUBLIC: Resilient Navigation and Timing Foundation (RNTF), National Institute of Standards and Technology (NIST), konstellasjon. AVLEDET: effekten i dBW. Representativt: bærebølge i båndets senter, avstandsfeil |
starlink-sop |
Bare Doppler: 11,325 GHz, 240 MHz fyr med ortogonal frekvensdelt multipleksing (OFDM), C/N0 opptil 57,07 dB-Hz, Doppler-feil 30 Hz; første skall på 1584 satellitter på 550 km og 53 grader | PUBLIC: Kozhaya, Saroufim og Kassas 2025, Federal Communications Commission (FCC) 21-48. Representativt: faseforskyvning, C/N0 ved horisonten |
centispace |
Binær faseskiftmodulasjon (BPSK) ved 2,046 Mchip/s nær L1 og L5; banebestemmelse på centimeternivå | PUBLIC: innledende vurdering, GPS Solutions. Representativt: bærebølger, effekt, hele den operative konstellasjonen |
generic-c-band |
Et C-båndssystem i tildelingen 5010 til 5030 MHz for radionavigasjonssatellitttjenesten (RNSS), for systemer der parametrene ikke er offentlige | PUBLIC: tildelingen i International Telecommunication Union (ITU) Radio Regulations. Alt annet er representativt |
atomic-zero-clock |
En efemeride-og-klokke-modell: efemeride med polynom av sjette orden, 30 s oppdatering, klokke styrt mot GNSS-tid; SISRE 0,26 m | PUBLIC: InsideGNSS. AVLEDET: SISRE. Kombiner med en hvilken som helst konstellasjon gjennom ephemeris_preset |
celeste-iod |
Celeste-demonstrasjonen i bane (IOD) fra Den europeiske romorganisasjonen (ESA), én konfigurasjon av en parametrisert signalutforming | PUBLIC bane (ESA, Institute of Navigation (ION) abstract 16907); WORKSHOP signalparametere, presentert på ESA Navigation Innovation and Support Programme (NAVISP) LEO-PNT-workshopen, 2026 |
Hvert tall som stammer fra workshopen ligger i modulen fusion i src/celeste_iod.rs, brukt av scenarios/celeste-iod-fused-pvt.toml (tests/workshop_preset_isolation.rs håndhever hvor de kan forekomme). Slik holdes forhåndsinnstillingen tilbake: hvert tall fra workshopen i motoren ligger i én fil, src/celeste_iod.rs, som build.rs bare kompilerer inn når den finnes, og i scenarioene scenarios/*celeste-iod*.toml som bare finnes i repositoriet. Når den filen og de scenarioene slettes, holdes forhåndsinnstillingen tilbake uten kildeendring; alle andre tester, uavhengige referanser og scenarioer kjører uten den, og antallet scenariofiler i README teller fortsatt de tilbakeholdte filene.
leo-pvt#
Én type, fire moduser (mode = "doppler" | "joint" | "polar" | "timing").
Doppler. En mottaker som måler Doppler f_D ved bølgelengde λ observerer avstandsraten ρ̇ = −λ f_D = u·(v_s − v_u) + ḋ, med u siktlinjen og ḋ mottakerklokkens drift. En batch-tilpasning med Gauss-Newton løser posisjonen, driften og eventuelt en konstant hastighet, med analytiske partiellderiverte ∂ρ̇/∂r = −(I − u uᵀ)(v_s − v_u)/ρ: bare satellitthastigheten på tvers av siktlinjen bærer posisjonsinformasjon, så geometrien er hvordan siktlinjene dreier seg i løpet av vinduet. Rapporten gir feilen mot vindulengde, omhyllingskurven for Doppler, Doppler-rate og rykk for hvert system, og en tabell for én passering. Én passering lokaliserer brukeren langs sporet (tidspunktet for null Doppler) og på tvers av det (hvor bratt Doppler-kurven er); feilen på tvers av sporet vokser uten grense når passeringen går rett over, og speilbildet på tvers av bakkesporet passer nesten like godt (bare jordens rotasjon skiller dem). Et system med bare Doppler (doppler_only = true, forhåndsinnstillingen starlink-sop) kjører i denne modusen.
Felles. For hver epoke, vektet minste kvadraters metode over pseudoavstandene til bare GNSS-systemene, bare LEO-systemene og alle sammen, én mottakerklokke per system (forskyvningen mellom systemer, ISB, estimeres) eller en kjent kringkastet forskyvning på referansetidsskalaen. Den rapporterer den geometriske, posisjons-, horisontale, vertikale og tidsmessige fortynningen av presisjon (GDOP, PDOP, HDOP, VDOP, TDOP), feilene, de estimerte forskyvningene, og DOP mot antallet LEO-satellitter lagt til GNSS-geometrien.
Polar. Satellitter i sikte, median DOP og tilgjengelighet mot breddegrad for MEO GNSS, LEO og begge.
Tidssynkronisering. På en kjent posisjon er hver LEO-pseudoavstand en klokkemåling. Mottakeren kjører et Kalman-filter med to tilstander (fase, frekvens) der prosesstøyen er oscillatorklassens hvite og tilfeldig-gange-frekvensstøy (slot-timing og klokkemodellene), holder tiden mellom passeringer, og omregner systemtid til koordinert universaltid (UTC) med IS-GPS-200-uttrykket Δt_UTC = Δt_LS + A0 + A1 (t_E − t_ot + 604800 (WN − WN_t)); usikkerheten i den kringkastede forskyvningen er en skjevhet per kjøring.
leo-ppp#
Et float-EKF (utvidet Kalman-filter) for presis punktposisjonering (PPP) på ionosfærefri kode og fase: statiske koordinater, en hvit mottakerklokke, én forskyvning mellom systemer per ekstra system, en zenitvåt forsinkelse som følger tilfeldig gange, og en float-flertydighet per satellittbue. Feilen i den presise banen og klokken er felles for en satellitts kode og fase, så hver satellitt er én korrelert oppdatering med to rader, i en Joseph-stabilisert form som holder kovariansen positiv definitt over tusenvis av faseoppdateringer på millimeternivå. Målingene simuleres fra samme modell, så den normaliserte estimeringsfeilen i kvadrat (NEES) for posisjonen må følge en khikvadratfordeling med tre frihetsgrader; testen kontrollerer det over Monte Carlo-seed. Konvergens er første epoke etter hvilken de horisontale og vertikale feilene holder seg under terskelverdiene (10 cm) til kjøringen er slutt.
En LEO-satellitt krysser himmelen på omtrent ti minutter, en MEO-satellitt flytter seg noen få grader, så LEO-buer skiller flertydighetene fra posisjonen mye raskere.
ntn-positioning#
En 5G-nedlenke i et ikke-terrestrisk nettverk (NTN) i S-båndet for mobil satellitttjeneste (MSS), 3rd Generation Partnership Project (3GPP)-bånd n256 (2170 til 2200 MHz). Cramér-Rao-grensen (CRB) for ankomsttid er c / (2π β sqrt(2 (C/N0) T)), med β rotmiddelkvadrat-båndbredden (Gabor) (B/√12 for et flatt OFDM-spekter), og for frekvensen til en tone sqrt(3 / (2π² (C/N0) T³)). Typen sammenligner kanaler (som standard en New Radio-kanal på 5 MHz og en smalbåndskanal på 200 kHz) gjennom posisjonsbestemmelser fra nedlenkens ankomsttid med en ukjent mottakerklokke og en Doppler-posisjonsbestemmelse med én satellitt over en passering.
Medfølgende scenarioer#
| Scenario | Type, modus | Hva den viser (med seed, deterministisk) |
|---|---|---|
leo-doppler-positioning |
leo-pvt, doppler |
Xona Pulsar X1-laget (258 satellitter), 1 Hz Doppler-støy: 3D-feil 7,5 m etter ti minutter (26 satellitter); på én enkelt passering er énsigma på tvers av sporet 303 m når passeringen går nesten rett over (25 km fra sporet) mot 43 m ved 250 km fra sporet |
starlink-sop-doppler-positioning |
leo-pvt, doppler |
Mulighetssignaler med bare Doppler, 30 Hz Doppler-feil ved 10 Hz: formelt horisontalt énsigma omtrent 10 m etter 20 s med 12 satellitter og 6 m etter 60 s |
meo-leo-fused-pvt |
leo-pvt, joint |
GPS og Galileo median PDOP 1,47 og rotmiddelkvadrat (RMS) 3D-feil 1,06 m; Xona X5 alene 1,73 og 0,49 m; fusjonert 1,06 og 0,41 m; LEO-skjevheten på 40 ns (11,99 m) estimert som 12,06 m |
leo-ppp-convergence |
leo-ppp |
Konvergens 7,4 min (bare MEO GNSS), 4,8, 3,2, 2,7 og 2,3 min med 60, 96, 192 og 288 LEO-satellitter |
ntn-5g-positioning |
ntn-positioning |
Avstandsgrense 0,26 til 1,03 m for 4,5 MHz og 6,5 til 26 m for 180 kHz; posisjonsbestemmelser nær 2 m mot 26 m |
polar-arctic-leo-coverage |
leo-pvt, polar |
GPS og Galileo VDOP 1,20 ved ekvator, 1,51 ved 89,9 grader; med Iridium 1,14 ved polen |
leo-timing-utc |
leo-pvt, timing |
Iridium-tidsoverføring over 6 timer ved Boulder: 9,9 til 11,4 ns RMS ved nominell C/N0 på tvers av fire oscillatorklasser, 15,1 til 25,0 ns 30 dB lavere |
celeste-iod-fused-pvt (valgfri) |
leo-pvt, joint |
GNSS sammen med de to IOD-satellittene på referansetidsskalaen (bare i repositoriet, tilbakeholdt sammen med forhåndsinnstillingen) |
Typene kjører også inne i andre scenarioer: leo-pvt er medlem av sveipet leo-focus-fused-pnt-sisre og av kampanjene leo-vertical-autonomous-vehicle, leo-vertical-polar-arctic, leo-vertical-critical-infrastructure-timing og leo-vertical-5g-network-timing; leo-ppp av sveipet leo-focus-ppp-altitude; ntn-positioning av sveipet leo-focus-ntn-bandwidth; og leo-pvt og leo-ppp er trinnene for fusjonert posisjonsbestemmelse og PPP i typen leo-pnt-chain. Resultatene deres står på LEO-PNT.md.
Publiserte tall og merkelapper#
| Tall | Kilde | Kshana | Merkelapp |
|---|---|---|---|
| Iridium-Doppler opptil ±36 kHz | RNTF | 35,9 kHz, maksimum over et døgn med passeringer, innenfor 5% | VALIDATED (tests/leo_doppler_reference.rs) |
| Xona Pulsar X1 maksimal Doppler 32 til 34 kHz | arXiv:2509,19551 | 33,6 kHz fra én satellitt i hvert publiserte skall (1080 km, 53 og 97 grader), innenfor intervallet uten utvidelse; skallet på 53 grader alene har topp på 31,8 kHz, siden jordens rotasjon bremser en prograd bane i forhold til bakken | VALIDATED (samme test) |
| Xona Pulsar X1 rykk ±1,26 Hz/s² | arXiv:2509,19551 | omtrent 1,0 Hz/s² på en passering rett over | rapportert, ikke validert |
| PPP-konvergens 9,6 → 7,0, 3,2, 2,1, 1,3 min med 60/96/192/288 LEO | Li et al., Journal of Geodesy 93:749 (2019), doi 10.1007/s00190-018-1195-2 | 7,4 → 4,8, 3,2, 2,7, 2,3 min, monoton | MODELLED konsistens i trenden (konstellasjoner, støy og stasjoner er ikke gjengitt) |
| Starlink omtrent 2 m på 20 s med tre satellitter | Kozhaya, Saroufim og Kassas 2025 | formelt horisontalt énsigma omtrent 10 m etter 20 s med 12 satellitter, med en perfekt efemeride og all feil samlet i et Doppler-sigma på 30 Hz ved 10 Hz | MODELLED sammenligning (mottaker og målerate er ikke gjengitt) |
| Iridium-tidssynkronisering under 40 ns fra UTC(NIST) | NIST | omtrent 10 ns RMS modellert | MODELLED sammenligning |
Alt annet er MODELLED med begrunnelsen i verifikasjonsmatrisen (docs/VERIFICATION-MATRIX.md): DOP kontrolleres mot et håndutledet tilfelle med fire satellitter, posisjonsbestemmelsene mot kovariansen sin, Doppler-modellen mot identiteten for avstandsakselerasjon, PPP-filteret mot khikvadrat-NEES-båndet, og grensene mot kvadratur.
Begrensninger, oppgitt uten omsvøp#
- Totrekkslegemebaner med den sekulære driften fra J2, jordens andre sonale harmoniske; ingen tilpasningsfeil i kringkastet efemeride utover SISRE-énsigma, ingen satellittantennemønster.
- Ingen ionosfærisk eller troposfærisk forsinkelse i
leo-pvt(PPP-filteret har en tilstand for våt forsinkelse; dets ionosfærefrie kombinasjon fjerner førsteordens ionosfære), ingen scintillasjon, flerveis, syklusglipp eller løsning av heltallsflertydighet. - NTN-nøyaktigheten er en grense, ikke et oppnådd resultat.
- Forhåndsinnstillingsverdier merket representative er illustrerende og beskriver ingen operatørs utforming.