03.06.2013 Views

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

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

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

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Efecto Doppler<br />

Navegación por posicionamiento<br />

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

GPS: Otros conceptos<br />

Ecuación <strong>de</strong>l efecto Doppler: fR = fT<br />

GPS: segmentos<br />

Cálculo <strong>de</strong> posición. Errores.<br />

Cálculo <strong>de</strong> velocidad.<br />

<br />

1 − v r · a<br />

c<br />

fR es la frecuencia recibida.<br />

fT es la frecuencia transmitida (conocida).<br />

<br />

, don<strong>de</strong>:<br />

v r = ˙s − ˙u es la velocidad relativa satélite-usuario.<br />

a = s−u<br />

s−u es el vector unitario en la dirección satélite-usuario.<br />

Si ya hemos obtenido la posición siguiendo los métodos<br />

anteriormente <strong>de</strong>scritos se pue<strong>de</strong> consi<strong>de</strong>rar a conocido.<br />

Por tanto la diferencia <strong>de</strong> fase ∆f vendrá dada por:<br />

v r · a<br />

∆f = fT − fR = −fT .<br />

c<br />

Por otro lado el observable no es directamente ∆f , porque el<br />

reloj <strong>de</strong>l segmento <strong>de</strong> usuario no tiene la suficiente precisión e<br />

introduce errores <strong>de</strong> medida <strong>de</strong> frecuencia <strong>de</strong> la siguiente<br />

forma: fM = fR − fM ˙ tu, don<strong>de</strong> fM es la frecuencia medida.<br />

Navegación por posicionamiento<br />

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

GPS: Otros conceptos<br />

Deriva <strong>de</strong>l reloj <strong>de</strong> usuario<br />

GPS: segmentos<br />

Cálculo <strong>de</strong> posición. Errores.<br />

Cálculo <strong>de</strong> velocidad.<br />

Para enten<strong>de</strong>r la ecuación fM = fR − fM ˙<br />

tu, imaginemos que el<br />

usuario mi<strong>de</strong> una señal dada por y = sin(a · τ), don<strong>de</strong> τ es el<br />

tiempo <strong>de</strong>l receptor.<br />

El receptor <strong>de</strong>l usuario <strong>de</strong>duce que tiene una señal <strong>de</strong> a<br />

Hercios. Por tanto fM = a.<br />

Pero si τ = t, don<strong>de</strong> t es el tiempo GPS, se introduce un<br />

error. Este error es τ = t + tu, don<strong>de</strong> tu es la <strong>de</strong>riva <strong>de</strong>l reloj<br />

<strong>de</strong>l usuario. Imaginemos que tu c1 + c2t. Luego ˙<br />

tu = c2.<br />

Entonces realmente y = sin(at + ac1 + ac2t), lo que es una<br />

señal <strong>de</strong> a + ac2 Hercios (ac1 es un <strong>de</strong>sfase y no influye en la<br />

frecuencia <strong>de</strong> la señal). Luego fR es igual a a + ac2.<br />

En efecto, se verifica: a = a + ac2 − ac2.<br />

33 / 48<br />

Aunque ˙<br />

tu pue<strong>de</strong> ser muy pequeño, tiene un efecto muy<br />

significativo en el resultado real, ya que estará multiplicado<br />

por c. Por tanto una <strong>de</strong>riva <strong>de</strong> 1 microsegundo por segundo<br />

(10 −6 ) daría errores <strong>de</strong>l or<strong>de</strong>n <strong>de</strong> 300 m/s! 34 / 48<br />

Navegación por posicionamiento<br />

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

GPS: Otros conceptos<br />

Algoritmo <strong>de</strong> cálculo <strong>de</strong> velocidad I<br />

GPS: segmentos<br />

Cálculo <strong>de</strong> posición. Errores.<br />

Cálculo <strong>de</strong> velocidad.<br />

Por tanto tenemos las dos ecuaciones: fT − fR = −fT<br />

fM = fR − fM ˙<br />

tu.<br />

Eliminando fR, se llega a:<br />

c fT − fM<br />

fT<br />

= −v r · a − c fM<br />

Puesto que v r = ˙s − ˙u, escribimos:<br />

c fT − fM<br />

fT<br />

Llamemos d = c fT −fM<br />

fT<br />

˙<br />

tu<br />

fT<br />

+ ˙s · a = ˙u · a − c fM<br />

˙<br />

tu<br />

fT<br />

+ ˙s · a; es un vector conocido en<br />

v r · a<br />

c<br />

función <strong>de</strong> los datos, la medida <strong>de</strong> fase <strong>de</strong> la portadora, el<br />

cálculo orbital, y la estimación anterior <strong>de</strong> la posición.<br />

Para cada satélite (un mínimo como ya vimos <strong>de</strong> 4) se<br />

tendrá una ecuación: d i = ˙u · ai − c fM i<br />

tu ˙<br />

fT i<br />

Navegación por posicionamiento<br />

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

GPS: Otros conceptos<br />

Algoritmo <strong>de</strong> cálculo <strong>de</strong> velocidad II<br />

GPS: segmentos<br />

Cálculo <strong>de</strong> posición. Errores.<br />

Cálculo <strong>de</strong> velocidad.<br />

Obsérvese que las componentes <strong>de</strong> ai a partir <strong>de</strong> una posición<br />

anteriormente estimada û son:<br />

axi =<br />

(xi −ˆxu)<br />

ayi =<br />

azi =<br />

√<br />

(xi −ˆxu) 2 +(yi −ˆyu) 2 ,<br />

+(zi −ˆzu) 2<br />

(yi −ˆyu)<br />

√<br />

(xi −ˆxu) 2 +(yi −ˆyu) 2 ,<br />

+(zi −ˆzu) 2<br />

(zi −ˆzu)<br />

√<br />

(xi −ˆxu) 2 +(yi −ˆyu) 2 .<br />

+(zi −ˆzu) 2<br />

Estos valores ya se habían calculado anteriormente en la<br />

estimación <strong>de</strong> posición! Luego son conocidos. Llegamos a:<br />

d i = axi ˙ux + ayi ˙uy + azi ˙uz − c fM i<br />

tu ˙<br />

fT i<br />

Aproximamos fM i 1. Por tanto llegamos a la ecuación<br />

fT i<br />

d = Hg, don<strong>de</strong> H es la matriz que se usó para estimar la<br />

posición, d es un vector con las medidas y datos, y<br />

g = [ ˙ux ˙uy ˙uz − c ˙tu] T que hay que calcular.<br />

Resolvemos el problema por mínimos cuadrados como antes.<br />

y<br />

35 / 48<br />

36 / 48

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

Saved successfully!

Ooh no, something went wrong!