Geodesia. Cartografía. Sistemas de referencia. Tiempos.

Geodesia. Cartografía. Sistemas de referencia. Tiempos. Geodesia. Cartografía. Sistemas de referencia. Tiempos.

03.06.2013 Views

GPS diferencial Navegación por posicionamiento GNSS: Navegación por satélite GPS: Otros conceptos Disponibilidad, integridad, continuidad Sistemas de aumento: GPS diferencial Cálculo de la actitud Para mejorar la precisión del GPS (o la integridad) se emplean las técnicas de GPS diferencial (DGPS). La idea básica es usar una o más estaciones (pseudollites), cuya posición se conoce con gran precisión, equipadas con un receptor GPS y en comunicación con el usuario (GBAS=Ground-Based Augmentation Systems). También se pueden emplear satélites extra que proporcionen medidas adicionales (SBAS=Space-Based Augmentation Systems). Por ejemplo, la red europea EGNOS. Los sistemas DGPS se clasifican como: Absolutos (ECEF) o relativos (posiciones relativas a la estación). Por zona geográfica de cobertura: Locales (10-100 km) Regionales (menos de 1000 km) Wide-area (más de 1000 km) Basados en pseudodistancias o en fases (en fases son más 8.2 Spatial and Time Correlation Characteristics of GPS Errors 381 extremely important, since they directly influence the performance achievable for any type of DGPS system. The underlying algorithms and performance of code- and carrier-based DGPS systems are presented in Sections 8.3 and 8.4, respectively. Some important DGPS message standards are introduced in Section 8.5. The final section, Section 8.6, details a number of operational and planned DGPS systems. 8.2 Spatial and Time Correlation Characteristics of GPS Errors Many of the GPS error sources discussed in Chapter 7 are highly correlated over space and time. All DGPS systems exploit these correlations to improve overall system performance. For instance, in a simple local-area DGPS system with a single reference station (see Figure 8.1), the errors in the reference station’s pseudorange and carrier-phase measurements for visible satellites are expected to be very similar to those experienced by a nearby user. If the reference station estimates the errors by leveraging its known surveyed position and provides this information in the form of corrections to the user, it is expected that the user’s position accuracy will be precisos, pueden conseguir precisión de mm.) improved as a result. This section quantifies the correlation of GPS errors between receivers separated over some distance (often referred to as the baseline, which may be interpreted as a vector) and over time. Time correlations (i.e., how rapidly the errors change with time), are also of interest, Navegación becausepor in general posicionamiento DGPS systems can- Disponibilidad, integridad, continuidad not instantaneously provide data to the end GNSS: user—even Navegación with a high-speed por satélite radio link Sistemas de aumento: GPS diferencial there is some finite delay associated with the generation, GPS: Otros transmission, conceptos reception, Cálculo de la actitud and application of the data. GPS diferencial: principios básicos de funcionamiento I 8.2.1 Satellite Clock Errors Satellite clock errors are one of the simplest GPS errors to correct. This is because a satellite clock error causes the same effect on pseudorange and carrier-phase measurements, regardless of the location of the user. For instance, if the satellite clock User Figure 8.1 Local-area DGPS concept. Satellite Reference station 45 / 48 Ejemplifiquemos el funcionamiento del DGPS con un caso simple: GBAS, absoluto, local, basado en distancias y con una sola estación. Recordemos que el usuario debe encontrar su posición u resolviendo el sistema de 4 o más ecuaciones ρi − ctu = s i − u + νu, donde νu son los errores de las señales recibidas por el usuario. Supongamos ahora que se tiene una estación (pseudollite) de posición m = [xm ym zm] T ; su distancia al satélite i es: R m i = s i − m = (xi − xm) 2 + (yi − ym) 2 + (zi − zm) 2 . Si tiene un error de reloj tm, las medidas de pseudodistancia en la estación serán: ρ m i − ctm = s i − m + νm 46 / 48 Navegación por posicionamiento GNSS: Navegación por satélite GPS: Otros conceptos Disponibilidad, integridad, continuidad Sistemas de aumento: GPS diferencial Cálculo de la actitud GPS diferencial: principios básicos de funcionamiento II La posición de la estación es fija y conocida . Por tanto conocemos la cantidad ∆ρm i = Rm i − ρm i = −ctm − νm. al receptor, y el receptor calcula La estación envía ∆ρ m i (ρi)corr = ρi + ∆ρ m i = s i − u + c(tu − tm) + (νu − νm). Si definimos tum como el error del reloj del receptor respecto al reloj de la estación, tum = tu − tm observamos que (ρi)corr = s i − u + ctum + (νu − νm). Por otro lado, ν ′ = νu − νm ≪ νu, porque νu y νm serán muy parecidos. Luego hemos conseguido reducir mucho el error. El nuevo tiempo que calculemos será con respecto a la estación. Pero la estación puede calcular su error respecto al satélite e incluirla en su mensaje de radio, de forma que tu = tum + tm. Luego recuperamos el tiempo GPS. Se consigue σUERE ≈ 0,3m + (1 − 6 cm) × (dEST −RECEP en km). 47 / 48 Navegación por posicionamiento GNSS: Navegación por satélite GPS: Otros conceptos Cálculo de la actitud mediante GPS Disponibilidad, integridad, continuidad Sistemas de aumento: GPS diferencial Cálculo de la actitud Con DGPS de precisión (basado en fases) se obtiene la actitud. Se sitúan n antenas receptoras en puntos separados de la aeronave y un único receptor. Se calculan las diferencias de posición entre antenas, asumiendo que los errores se cancelan. A cada una de estas diferencias las llamamos r i. Por ejemplo si hay 3 antenas habrá 3 medidas. Si hay 4 antenas habrá 6 medidas. En general serán n!/(2 × (n − 2)!) medidas. Evidentemente r b i es conocido. Lo que se mide es r e i . Suponiendo la posición conocida, calculamos r n i = C n e r e i . Queremos calcular C b n de las ecuaciones r b i = C b n r n i . Formemos las matrices Rb = [r b 1r b 2 . . . r b n] y Rn = [r n 1r n 2 . . . r n n]. Se tendrá que Rb = C b n Rn . Si las medidas fueran tres, se podría hacer C b n = Rb (Rn ) −1 , pero conviene mejor Rn = (C b n ) T Rb , luego C b n = (Rn (Rb ) −1 ) T . En general se usará una solución de mínimos cuadrados para más medidas. El mínimo de medidas necesarias son 3; por tanto el mínimo 48 / 48 de antenas necesario será de 3.

Sistemas de navegación integrados Filtrado óptimo de sistemas lineales: el filtro de Kalman. Navegación Aérea Tema 6: Sistemas de navegación integrados. El filtro de Kalman. Sistemas de navegación integrados Filtrado óptimo de sistemas lineales: el filtro de Kalman. Fusión de sensores. Fusión de sensores. Ejemplo: el canal vertical. INS-GPS Una aeronave actual dispone de una gran diversidad de sensores y sistemas de navegación, que pueden obtener total o parcialmente las variables de navegación PVAT. Por ejemplo hemos visto el INS, que a partir de las medidas de la IMU, el modelo de Tierra y gravedad, y una estimación inicial, nos da posición, velocidad y actitud en todo momento. También hemos visto el GPS, que igualmente es capaz de darnos todos éstos datos, o al menos (si no disponemos de múltiples antenas), la posición y la velocidad. Puede haber otros sistemas (DME-DME, etc...) Cada sistema dará una estimación diferente, sujeta a error. La idea de fusión de sensores y de los sistemas de navegación integrados, consiste en obtener una única estimación PVAT a partir de todas las anteriores, tal que el error sea el menor posible. 2 / 26 Sistemas de navegación integrados Filtrado óptimo de sistemas lineales: el filtro de Kalman. Ejemplo: el canal vertical. Fusión de sensores. Ejemplo: el canal vertical. INS-GPS Se vio en el tema 4 que el canal vertical del INS es inestable. Una forma de estabilizar el canal es usar la medida de altitud obtenida de medidas barométricas, hB. Se denomina “estimador baro-inercial de la altitud”. Recordemos que las ecuaciones del canal vertical venían dadas por: ˙ˆh = − ˆVD, ˙ˆVD = ˆρz + µe , (Re + ˆh) 2 donde ˆρz es la componente z de −(ˆω n n/e + 2ˆωn e/i )× ˆv n + â n NG . Sistemas de navegación integrados Filtrado óptimo de sistemas lineales: el filtro de Kalman. Estimador baro-inercial de la altitud I Fusión de sensores. Ejemplo: el canal vertical. INS-GPS Se modifica el canal vertical del INS de la siguiente forma, usando hB: ˙ˆh = − ˆVD − C1(ˆh − hB), µe ˙ˆVD = ˆρz + (Re + ˆ h) 2 + C2( ˆ h − hB) + C3 t donde C1, C2 y C3 son ganancias a determinar. Calculando como en el tema 4 el error de altitud y despreciando el error en el término ρz, obtenemos: δ ˙h = −δVD + C1(ˆh − hB), δ ˙VD ≈ − 2g0 δh − C2(ˆh − hB) − C3 Re t 0 0 3 / 26 ( ˆ h(τ) − hB(τ))dτ, (ˆh(τ) − hB(τ))dτ, y obsérvese que ˆh − hB = ˆh − h + h − hB = −(δh − δhB), donde δhB es el error de estimación barométrico, que suponemos aproximadamente constante. 4 / 26

GPS diferencial<br />

Navegación por posicionamiento<br />

GNSS: Navegación por satélite<br />

GPS: Otros conceptos<br />

Disponibilidad, integridad, continuidad<br />

<strong>Sistemas</strong> <strong>de</strong> aumento: GPS diferencial<br />

Cálculo <strong>de</strong> la actitud<br />

Para mejorar la precisión <strong>de</strong>l GPS (o la integridad) se emplean<br />

las técnicas <strong>de</strong> GPS diferencial (DGPS).<br />

La i<strong>de</strong>a básica es usar una o más estaciones (pseudollites),<br />

cuya posición se conoce con gran precisión, equipadas con un<br />

receptor GPS y en comunicación con el usuario<br />

(GBAS=Ground-Based Augmentation Systems).<br />

También se pue<strong>de</strong>n emplear satélites extra que proporcionen<br />

medidas adicionales (SBAS=Space-Based Augmentation<br />

Systems). Por ejemplo, la red europea EGNOS.<br />

Los sistemas DGPS se clasifican como:<br />

Absolutos (ECEF) o relativos (posiciones relativas a la<br />

estación).<br />

Por zona geográfica <strong>de</strong> cobertura:<br />

Locales (10-100 km)<br />

Regionales (menos <strong>de</strong> 1000 km)<br />

Wi<strong>de</strong>-area (más <strong>de</strong> 1000 km)<br />

Basados en pseudodistancias o en fases (en fases son más<br />

8.2 Spatial and Time Correlation Characteristics of GPS Errors 381<br />

extremely important, since they directly influence the performance achievable for<br />

any type of DGPS system. The un<strong>de</strong>rlying algorithms and performance of co<strong>de</strong>- and<br />

carrier-based DGPS systems are presented in Sections 8.3 and 8.4, respectively.<br />

Some important DGPS message standards are introduced in Section 8.5. The final<br />

section, Section 8.6, <strong>de</strong>tails a number of operational and planned DGPS systems.<br />

8.2 Spatial and Time Correlation Characteristics of GPS Errors<br />

Many of the GPS error sources discussed in Chapter 7 are highly correlated over<br />

space and time. All DGPS systems exploit these correlations to improve overall system<br />

performance. For instance, in a simple local-area DGPS system with a single reference<br />

station (see Figure 8.1), the errors in the reference station’s pseudorange and<br />

carrier-phase measurements for visible satellites are expected to be very similar to<br />

those experienced by a nearby user. If the reference station estimates the errors by<br />

leveraging its known surveyed position and provi<strong>de</strong>s this information in the form of<br />

corrections to the user, it is expected that the user’s position accuracy will be<br />

precisos, pue<strong>de</strong>n conseguir precisión <strong>de</strong> mm.)<br />

improved as a result. This section quantifies the correlation of GPS errors between<br />

receivers separated over some distance (often referred to as the baseline, which may<br />

be interpreted as a vector) and over time. Time correlations (i.e., how rapidly the<br />

errors change with time), are also of interest, Navegación becausepor in general posicionamiento<br />

DGPS systems can- Disponibilidad, integridad, continuidad<br />

not instantaneously provi<strong>de</strong> data to the end GNSS: user—even Navegación with a high-speed por satélite radio link <strong>Sistemas</strong> <strong>de</strong> aumento: GPS diferencial<br />

there is some finite <strong>de</strong>lay associated with the generation, GPS: Otros transmission, conceptos reception, Cálculo <strong>de</strong> la actitud<br />

and application of the data.<br />

GPS diferencial: principios básicos <strong>de</strong> funcionamiento I<br />

8.2.1 Satellite Clock Errors<br />

Satellite clock errors are one of the simplest GPS errors to correct. This is because a<br />

satellite clock error causes the same effect on pseudorange and carrier-phase measurements,<br />

regardless of the location of the user. For instance, if the satellite clock<br />

User<br />

Figure 8.1 Local-area DGPS concept.<br />

Satellite<br />

Reference<br />

station<br />

45 / 48<br />

Ejemplifiquemos el funcionamiento <strong>de</strong>l<br />

DGPS con un caso simple: GBAS,<br />

absoluto, local, basado en distancias y<br />

con una sola estación.<br />

Recor<strong>de</strong>mos que el usuario <strong>de</strong>be<br />

encontrar su posición u resolviendo el<br />

sistema <strong>de</strong> 4 o más ecuaciones<br />

ρi − ctu = s i − u + νu, don<strong>de</strong> νu<br />

son los errores <strong>de</strong> las señales recibidas<br />

por el usuario.<br />

Supongamos ahora que se tiene una estación (pseudollite) <strong>de</strong><br />

posición m = [xm ym zm] T ; su distancia al satélite i es:<br />

R m i = s i − m = (xi − xm) 2 + (yi − ym) 2 + (zi − zm) 2 .<br />

Si tiene un error <strong>de</strong> reloj tm, las medidas <strong>de</strong> pseudodistancia<br />

en la estación serán: ρ m i − ctm = s i − m + νm<br />

46 / 48<br />

Navegación por posicionamiento<br />

GNSS: Navegación por satélite<br />

GPS: Otros conceptos<br />

Disponibilidad, integridad, continuidad<br />

<strong>Sistemas</strong> <strong>de</strong> aumento: GPS diferencial<br />

Cálculo <strong>de</strong> la actitud<br />

GPS diferencial: principios básicos <strong>de</strong> funcionamiento II<br />

La posición <strong>de</strong> la estación es fija y conocida .<br />

Por tanto conocemos la cantidad<br />

∆ρm i = Rm i − ρm i = −ctm − νm.<br />

al receptor, y el receptor calcula<br />

La estación envía ∆ρ m i<br />

(ρi)corr = ρi + ∆ρ m i = s i − u + c(tu − tm) + (νu − νm).<br />

Si <strong>de</strong>finimos tum como el error <strong>de</strong>l reloj <strong>de</strong>l receptor respecto<br />

al reloj <strong>de</strong> la estación, tum = tu − tm observamos que<br />

(ρi)corr = s i − u + ctum + (νu − νm).<br />

Por otro lado, ν ′ = νu − νm ≪ νu, porque νu y νm serán muy<br />

parecidos. Luego hemos conseguido reducir mucho el error.<br />

El nuevo tiempo que calculemos será con respecto a la<br />

estación. Pero la estación pue<strong>de</strong> calcular su error respecto al<br />

satélite e incluirla en su mensaje <strong>de</strong> radio, <strong>de</strong> forma que<br />

tu = tum + tm. Luego recuperamos el tiempo GPS.<br />

Se consigue<br />

σUERE ≈ 0,3m + (1 − 6 cm) × (dEST −RECEP en km). 47 / 48<br />

Navegación por posicionamiento<br />

GNSS: Navegación por satélite<br />

GPS: Otros conceptos<br />

Cálculo <strong>de</strong> la actitud mediante GPS<br />

Disponibilidad, integridad, continuidad<br />

<strong>Sistemas</strong> <strong>de</strong> aumento: GPS diferencial<br />

Cálculo <strong>de</strong> la actitud<br />

Con DGPS <strong>de</strong> precisión (basado en fases) se obtiene la<br />

actitud. Se sitúan n antenas receptoras en puntos separados <strong>de</strong><br />

la aeronave y un único receptor. Se calculan las diferencias <strong>de</strong><br />

posición entre antenas, asumiendo que los errores se cancelan.<br />

A cada una <strong>de</strong> estas diferencias las llamamos r i. Por ejemplo<br />

si hay 3 antenas habrá 3 medidas. Si hay 4 antenas habrá 6<br />

medidas. En general serán n!/(2 × (n − 2)!) medidas.<br />

Evi<strong>de</strong>ntemente r b i es conocido. Lo que se mi<strong>de</strong> es r e i .<br />

Suponiendo la posición conocida, calculamos r n i = C n e r e i .<br />

Queremos calcular C b n <strong>de</strong> las ecuaciones r b i = C b n r n i .<br />

Formemos las matrices Rb = [r b 1r b 2 . . . r b n] y Rn = [r n 1r n 2 . . . r n n].<br />

Se tendrá que Rb = C b n Rn . Si las medidas fueran tres, se<br />

podría hacer C b n = Rb (Rn ) −1 , pero conviene mejor<br />

Rn = (C b n ) T Rb , luego C b n = (Rn (Rb ) −1 ) T . En general se<br />

usará una solución <strong>de</strong> mínimos cuadrados para más medidas.<br />

El mínimo <strong>de</strong> medidas necesarias son 3; por tanto el mínimo<br />

48 / 48<br />

<strong>de</strong> antenas necesario será <strong>de</strong> 3.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!