ES2250675T3 - Modulo discriminador de trayectos multiples para un sistema de navegacion. - Google Patents

Modulo discriminador de trayectos multiples para un sistema de navegacion.

Info

Publication number
ES2250675T3
ES2250675T3 ES02748961T ES02748961T ES2250675T3 ES 2250675 T3 ES2250675 T3 ES 2250675T3 ES 02748961 T ES02748961 T ES 02748961T ES 02748961 T ES02748961 T ES 02748961T ES 2250675 T3 ES2250675 T3 ES 2250675T3
Authority
ES
Spain
Prior art keywords
cdot
signal
module
otimes
code
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
ES02748961T
Other languages
English (en)
Inventor
Robert Schweikert
Thomas Worz
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Agence Spatiale Europeenne
Original Assignee
Agence Spatiale Europeenne
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Agence Spatiale Europeenne filed Critical Agence Spatiale Europeenne
Application granted granted Critical
Publication of ES2250675T3 publication Critical patent/ES2250675T3/es
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7097Interference-related aspects
    • H04B1/711Interference-related aspects the interference being multi-path interference
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/22Multipath-related issues

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Instructional Devices (AREA)
  • Noise Elimination (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

Módulo discriminador de trayectos múltiples para un sistema de comunicación y/o de navegación que emplea una modulación con espectro repartido, módulo que presenta una entrada capaz de recibir señales de navegación, un muestreador para suministrar señales muestreadas k a una frecuencia doble de la frecuencia aparente fc del código de dichas señales, así como un submódulo (DISCR) de cálculo de una señal de error ek a partir de dichas señales muestreadas k y del késimo valor del código de reparto generado localmente, el cual es designado mediante CIKIL, caracterizado porque el submódulo (DISCR) emplea el cálculo de una relación entre valores de correlación que incluye dos valores de correlación (Z-K, Z-K) adelantados con relación al valor de fase real, procediendo dichos valores de una correlación entre unas muestras de la señal recibida y dicha señal de código de reparto generada localmente, según una de las siguientes fórmulas: ek = Kâ u (Zk -- / Zk - + Sâ) ek = Kâ (Z--k / Z-k + Sâ) siendono nulo e incluido entre -1 y +1 ek = Kâ (2 Z--k / Z-k + Z+k + Sâ) donde { Z+K = [r k+0, 5 uCKL]h b/k Z-K = [r k-t2 uCKL]h b/k Z--K = [r k-t1 uCKL]h b/k o { Z+K = [rk uCk + 1/2L] h b/k Z-K = [rk uC (k-t2) L] h b/k Z--K = [rk uC(k-t1) L] h b/k K = constante hbk designa la frecuencia impulsional de un filtro de paso bajo designa la parte real Sâ = - g(-t1Tc) / g(-t2Tc) donde g (aTc) = SenÑaCosÑâa / ða[1-(2âa)] designa el factor de atenuación de la señal SRC.

Description

Módulo discriminador de trayectos múltiples para un sistema de navegación.
La presente invención tiene por objeto un módulo discriminador de trayectos múltiples para un sistema de navegación, así como un sistema de navegación que incluya dicho módulo.
En los últimos diez años, los sistemas de modulación con espectro repartido y secuencia directa (DS/SS) han adquirido una importancia creciente.
Actualmente, esta técnica se emplea no sólo en los sistemas de navegación por satélite, GPS y GLONASS, sino que se ha introducido también en los sistemas de comunicación terrestre y por satélite, como por ejemplo, US-Standard IS-95, GLOBALSTAR y, más recientemente, en la tercera generación de teléfonos móviles según el estándar UMTS, así como en el sistema europeo de navegación por satélite GALILEO.
El concepto de la modulación DS/SS, por ejemplo, la modulación bifásica BPSK, emplea la introducción de un código de ruidos pseudo-aleatorios PRN, que tiene como consecuencia que la señal modulada resultante presenta una pasabanda más ancha que una señal que transmitiera únicamente los símbolos de datos. Es en este sentido que se "reparte" la densidad espectral de la señal.
En el receptor, una réplica generada localmente del código PRN transmitido se alinea con la fase del código de la señal recibida. En particular, en los receptores de navegación, un alineamiento de la fase del código es esencial para determinar de forma precisa el tiempo de llegada TOA que se utiliza para la determinación de la distancia geométrica entre el transmisor y el receptor. Una vez realizado el alineamiento, es posible llevar a cabo la estimación de la fase de la portadora y la determinación de los símbolos de los datos transmitidos.
Dicho alineamiento se realiza de forma clásica en un receptor, por medio de un bucle cerrado en fase DLL ("Delay-Locked-Loop"), del que se describe un ejemplo en el artículo de M. SIMON et al. publicado en el libro "Spread Spectrum Communications Handbook" publicado por McGRAWHILL Incorporated, 2ª edición 1994.
Dicho alineamiento utiliza el resultado de una correlación entre la señal recibida y unas versiones anticipadas E y retardadas L de una señal de código de referencia generado localmente para calcular una señal de error que es proporcional al error de fase del código (diferencia entre la fase del código estimada y la fase recibida).
Dicha señal de error debe indicar en qué dirección debe estar desfasada la fase de la señal de referencia (adelanto o retraso) para estar en sincronización con la señal recibida. Un espacio entre los códigos adelantados E y retrasados L es generalmente de un bit de una secuencia de pseudo-ruido ("chip").
Las señales del tipo SRC (que tienen como espectro una cosinusoide sobreelevada) están definidas en el estándar UMTS. Este tipo de señal SRC puede adoptarse asimismo en el sistema GALILEO mencionado anteriormente. Se describe un receptor digital que empleaduchas señales, en el artículo de R. de GAUDENZI et al. titulado "A Digital Chip Timing Recovery Loop for Band-Limited Direct-Sequence Spread-Spectrum Signals", publicado en IEEE Trans. On Comm., vol. 41, nº 11, páginas 1760-1769, noviembre de 1993.
La precisión de la medición del tiempo de llegada es perturbada de forma negativa por la presencia de una distorsión debida a trayectos múltiples, de lo que resulta, en el caso de la telemetría, una disminución de la precisión de la determinación de la posición y, en el caso de la transmisión de datos, un incremento de la tasa de error de bit o de trama. Esto ocurre especialmente en el caso en que la distorsión multi-trayecto se representa especialmente por medio de una
reflexión única que procede de un punto que se sitúa en el entorno inmediato del receptor, con una dinámica escasa.
La superposición de las señales directas y reflejadas puede, de este modo, causar distorsión fortuita ("jitter") que afecta a las mediciones del tiempo de llegada TOA del bucle DLL.
De ello resulta que unas técnicas que permitan reducir el impacto de dichos trayectos múltiples sobre la determinación de la fase del código, tienen un interés primordial, especialmente en el ámbito de la navegación.
Hasta ahora, los procedimientos de compensación de trayectos múltiples destinados a la telemetría se han desarrollado especialmente en el marco de receptores GPS.
Por ello, un gran número de estos algoritmos utiliza el hecho de que el caudal aparente ("Chip Rate") del código C/A disponible públicamente es mucho más escaso que la pasabanda de la transmisión. En este caso, es ventajoso reducir el desfase temporal entre las señales de código de referencia adelantadas E y retrasadas L, hasta que tengan un valor que es inferior a un bit una secuencia de pseudo-ruidos ("One Chip Duration") para disminuir el error inducido por los trayectos multi-haces.
En el marco de los sistemas del tipo SRC, dado que el espectro de frecuencias está estrictamente limitado a 1 + \beta veces el caudal aparente del código (\beta designa el factor de atenuación de los impulsos SRC), los anteriores métodos no son eficaces para compensar los trayectos multi-haces.
Además, en el artículo de Phillip G. MATTOS titulado "Multipath elimination for the low-cost consumer GPS" publicado en el Acta de la conferencia ION GPS 1996, de Kansas City, páginas 665-671, se sugiere la sustitución de los puntos de correlación adelantados E y retrasados L mediante dos puntos de correlación adelantados. Por el contrario, este artículo no propone medios de implementación de esta técnica.
La presente invención tiene por objeto un módulo discriminador que convenga para su uso en un sistema de comunicación o de navegación con espectro repartido y, más concretamente, según una modulación del tipo SRC.
Según una primera variante, la invención se refiere a un módulo discriminador de trayectos múltiples para un sistema de comunicación y/o de navegación, que emplea una modulación con espectro repartido, módulo que presenta una entrada capaz de recibir señales de navegación, un muestreador para suministrar señales muestreadoras con una frecuencia el doble de la frecuencia aparente fc del código de dichas señales, así como un submódulo de cálculo de una señal de error e_{k} a partir de dichas señales muestreadas y de un código de reparto CIVIL generado localmente, caracterizado porque:
e_{k} = K_{\beta} \cdot \Re (Z_{k} - - / Z_{k} - + S_{\beta})
\vskip1.000000\baselineskip
donde \{Z^{-}{}_{k} = [r_{k-t2} \cdot C_{|K|L}] \otimes h \ b/k
\hskip1cm
Z^{- -}{}_{k} = [r_{k-t1} \cdot C_{|K|L}] \otimes h \ b/k
\vskip1.000000\baselineskip
\vskip1.000000\baselineskip
\hskip0.5cm
o \{Z^{-}{}_{k} = [r_{k} \cdot C_{|K - t2|L}] \otimes h \ b/k
\hskip1cm
Z^{- -}{}_{k} = [r_{k} \cdot C_{|K - t1|L}] \otimes h \ b/k
K\beta = constante
y
S_{\beta} = - g(-t1Tc) / g(-t2Tc)
donde
g (aTc) = Sen \Pi a \ Cos \Pi \beta a / \pi a[1-(2 \beta a)]
designando \beta el factor de atenuación de la señal SRC.
Según una segunda variante, la invención se refiere a un módulo discriminador de trayectos múltiples para un sistema de comunicación y/o de navegación que emplea una modulación con espectro repartido, módulo que presenta una entrada capaz de recibir señales de navegación, un muestreador para suministrar señales muestreadoras a una frecuencia el doble de la frecuencia aparente fc del código de dichas señales, así como un submódulo de cálculo de una señal de error e_{k} a partir de dichas señales muestreadas y de un código de reparto CIVIL generado localmente, caracterizado porque:
e_{k} = K_{\beta} (Z_{k} - - / Z_{k} - + S_{\beta})
\vskip1.000000\baselineskip
donde \{Z^{-}{}_{k} = [r_{k - t2} \cdot C_{|K|L}] \otimes h \ b/k
\hskip1cm
Z^{--}{}_{k} = [r_{k - t1} \cdot C_{|K|L}] \otimes h \ b/k
\vskip1.000000\baselineskip
\vskip1.000000\baselineskip
\hskip0.5cm
o \{Z^{-}{}_{k} = [r_{k} \cdot C_{|K - t2|L}] \otimes h \ b/k
\hskip1cm
Z^{--}{}_{k} = [r_{k} \cdot C_{|K-t1|L}] \otimes h \ b/k
K\beta = constante
y
S_{\beta} = g(-t1Tc) / g(-t2Tc)
donde
g (aTc) = Sen \Pi a \ Cos \Pi \beta a / \pi a[1 - (2 \beta a)]
designando \beta el factor de atenuación de la señal SRC.
Según una tercera variante, la invención se refiere a un módulo discriminador de trayectos múltiples para un sistema de comunicación y/o de navegación que emplea una modulación con espectro repartido, módulo que presenta una entrada capaz de recibir señales de navegación, un muestreador para suministrar señales muestreadoras a una frecuencia el doble de la frecuencia aparente fc del código de dichas señales, así como un submódulo de cálculo de una señal de error e_{k} a partir de dichas señales muestreadas y de un código de reparto CIVIL generado localmente, caracterizado porque:
e_{k} = K_{\beta} \Re (2 Z^{- -}{}_{k} / Z^{-}{}_{k} + Z^{+}{}_{k} + S_{\beta})
\vskip1.000000\baselineskip
donde \{Z^{+}{}_{K} = [r \ _{k + 1/2} \cdot C|K|L] \otimes h \ b/k
\hskip1.3cm
Z^{+}{}_{K} = [r \ k - t2 \cdot C|K|L] \otimes h \ b/k
\hskip1.3cm
Z^{-}{}_{K} = [r \ k - t1 \cdot C|K|L] \otimes h \ b/k
\vskip1.000000\baselineskip
\vskip1.000000\baselineskip
\hskip1cm
o \{Z^{+}{}_{K} = [rk \cdot C (k + 1 / 2)L] \otimes h \ b/k
\hskip1cm
Z^{-}{}_{K} = [rk \cdot C |k-t2|L] \otimes h \ b/k
\hskip1cm
Z^{--}{}_{K} = [rk \cdot C|k-t1|L] \otimes h \ b/k
\beta = constante
y
S_{\beta} = g(-t1Tc) / g(-t2Tc)
donde
g (aTc) = Sen\Pi a Cos\Pi \beta a / \pi a[1-(2\beta a)]
designando \beta el factor de atenuación de la señal SRC.
En cada unos de los casos anteriores, h^{b}_{k} designa la respuesta impulsional de un filtro de paso bajo.
El módulo discriminador puede caracterizarse porque:
K\beta = {1 / d/d\varepsilon (g((\varepsilon - t1)Tc) / g((\varepsilon - t2)Tc))_{\varepsilon = 0}}
Para cada una de las tres variantes anteriores, se puede obtener t1 = 1,5 y t2 = 0,5
La invención se refiere asimismo a un sistema de navegación caracterizado porque presenta un módulo discriminador como el definido anteriormente.
Finalmente, la invención se refiere a un sistema de navegación caracterizado porque presenta un módulo discriminador que genera una señal de error e_{k} que asegura, de manera conocida en sí, una corrección en bucle cerrado a partir de señales muestreadas Z^{+}_{K} y Z^{-}_{K}, al que está asociado un módulo discriminador como el definido anteriormente, para generar en bucle abierto a partir de una señal de error e_{k}, una señal de corrección de la salida de fase del código.
Otras características y ventajas de la invención aparecerán mejor mediante la lectura de la siguiente descripción, realizada a título de ejemplo no limitativo, con relación a los dibujos, en los cuales:
- la figura 1 ilustra la influencia de una señal multi-trayectos que incluye un trayecto directo y por lo menos un trayecto reflejado;
- la figura 2 muestra un bloque diagrama de un bucle DLL según el artículo de R. de GAUDENZI et al. indicado anteriormente;
\newpage
- las figuras 3a y 3b representan una curva de la función de autocorrelación g(t), respectivamente para un módulo según la figura 1 y según la invención;
- la figura 4 representa un bloque diagrama de un bucle DLL que incorpora un discriminador según la invención; y
- la figura 5 muestra un bloque diagrama de un bucle DLL que incorpora un discriminador clásico, así como un módulo discriminador según la presente invención.
Una señal s_{T}(t) en banda de base con espectro repartido en secuencia directa (DS-SS) se representa mediante la fórmula 3-1, con la palabra de código repartido de una longitud de 2 bits de secuencia de pseudo-ruido (o "chips") y presentando cada símbolo de dato dp,q,i M/L palabras de código repartidas
(3-1)ST(t) = \surd a \cdot Ps \cdot \sum_{i = -\infty}{}^{\infty} (d \ p, \lfloor i \rfloor_{M} \cdot CP,|i|_{L} + j \cdot b \cdot dQ,|i|_{M} \cdot CQ,|i|_{L}) \cdot gT(t - iTc)
Los factores a y b están en los siguientes casos.
a = 1,0 y b = 0 BPSK DS - SS
a = 0,5 y b = 1, con dp_{\cdot \ \lfloor \ i \rfloor \ M} = dQ_{\cdot \lfloor \ i \rfloor \ M} QPN DS - SS
a = 0,5 y b = 1, con dp_{\cdot \ \lfloor \ i \rfloor \ M} = dQ_{\cdot \lfloor \ i \rfloor \ M} Bi-BPSKDS-SS
donde,
P_{s}
potencia transmitida
d_{P/Q,i}
símbolos de datos
\quad
(d_{P/Q,i}\varepsilon \ [-1,1])
C_{P/Q,i} bit ("chip") de una palabra de código repartido de longitud
\quad
L (c_{P/Q,i}\varepsilon \ [-1,1])
T_{c} = 1/f_{c}
duración de un bit de una secuencia de pseudo-ruido (o duración de un "chip")
g_{T}(t)
forma del impulso de un bit o "chip", por ejemplo SRC
M
longitud de los símbolos de datos en la duración de un bit o "chip"
\lfloor i\rfloor_{M} = int(i/M)
|i|_{M}= i \ mod \ M
Tras la transmisión por un canal con adición de ruido blanco Gausiano (AWGN) que presenta una densidad espectral simétrica No/2, la señal recibida filtrada r(t) se muestra mediante la fórmula (3-2)
\quad
r(t) = S_{\tau}(t - \tau) \otimes g_{r}(t)
\quad
= \surd a \cdot P_{s} \cdot \sum_{i = -\infty}{}^{\infty} (d_{P} \cdot \lfloor i\rfloor_{M} \cdot C_{P} \cdot |i|_{L} + j \cdot b \cdot d_{Q} \cdot \lfloor i \rfloor_{M} \cdot C_{Q} \cdot |i|_{L})
\cdot g(t - \tau - iT_{c}) \cdot exp(j(\Delta \omega(t) + \pi(t))) + n(t)
(3-2)
g(t) = gT(t) \otimes gT(t) designa la forma del impulso de un bit o "chip" tras el filtrado y \otimes designa una convolución.
g(t) constituye pues la función de autocorrelación de gT(t).
\Delta\omega(t) designa el desfase residual de la frecuencia de portadora tras la mezcla de la señal en la banda de base.
\phi(t) designa la fase de la portadora de la señal recibida.
Debido a la similitud entre los componentes en fase y en cuadratura de la señal según la fórmula 3-2, se opta por conservar únicamente la componente en fase para la continuación. Por ello, los símbolos p,q que designan estas dos situaciones no se han conservado para aligerar la anotación.
Por lo tanto, la fórmula 3-2 se convierte en:
(3-3)r(t) = \surd P_{s} \cdot \sum_{i=-\infty}{}^{\infty} d\lfloor i\rfloor_{M} \cdot c|i|_{L} \cdot g(t - \tau - iT_{c}) \cdot exp(j(\Delta\omega (t) + \phi (t))) + n(t)
para un sistema del tipo BPSK DS-SS, por ejemplo.
La figura 1 es un modelo de trayecto múltiple utilizable para mostrar el efecto de los trayectos múltiples en el bucle DLL. Además de la señal directa procedente del satélite, la antena recibe una segunda versión retrasada de esta misma señal, denominada componente "multi-trayectos", debida a una reflexión, procediendo el retraso de esta segunda señal del trayecto más largo que ha debido recorrer.
La suma de estas dos señales recibidas por la antena puede expresarse mediante la siguiente fórmula:
(3-4)S'_{T}(t) = S_{T}(t - \tau) + \alpha \cdot S_{T} (t - \tau - \Delta \tau) \cdot exp (j \phi)
donde
a atenuación de la señal reflejada con relación a la señal directa
\Delta \tau retraso de la señal reflejada con relación a la señal directa
\varphi = 2\pi \cdot \Delta \tau \cdot f_{c}/c desfase de fase de la portadora de la señal reflejada con relación a la señal directa (c: velocidad de la luz, fc: frecuencia de la portadora)
Tras el filtrado (véase las fórmulas 3-2 y 3-3), se obtiene:
(3-5)r(t) = r(t) + \alpha - r(t - \Delta \tau) - exp(j\phi)
Se ha dado la explicación teniendo en cuenta una única reflexión. En la práctica, existen varios componentes que corresponden a trayectos reflejados múltiples que se superponen a la señal directa.
En el documento de R. de GAUDENZI et al. ya mencionado, se emplea un bucle DLL cuya arquitectura se recuerda a continuación, con relación a la figura 2.
Se observa:
-
que el convertidor analógico-digital podría estar dispuesto en otro lugar;
-
que el desfase temporal destinado a alinear la señal de reflexión local sobre la señal recibida puede aplicarse, bien a la señal de referencia local, bien a la señal recibida.
La señal en banda de base r(t) obtenida mediante filtrado con un filtro de características GT(f) se muestrea a dos veces la frecuencia de bit o "chip", es decir, a 2fc. Las muestras correspondientes a instantes "semi-enteros" (k + 0,5) Tc + \tau son dirigidas hacia el bucle DLL, mientras que las demás muestras correspondientes a instantes "enteros" kTc + \tau son dirigidas hacia la correlación puntual y el seguimiento de la fase de la portadora y la desmodulación de los datos (circuito NCO).
Las muestras correspondientes a instantes semi-enteros se dan mediante:
(3-6)rk + 1/2 = \surd P_{S} \cdot exp(j\phi) \cdot \sum_{i=\infty}{}^{\infty} d\lfloor i\rfloor M \cdot c|i|L \cdot g((\varepsilon k + k + 1/2 - i)Tc) \cdot + nk + 1/2
Estas muestras r_{k+1/2} son dirigidas según dos ramas. En la rama superior (en la figura 2), las muestras son retrasadas de un bit o "chip" Tc antes de ser multiplicadas por el késimo valor del código de reparto C_{|k|L} generado localmente por el generador de código SCGEN. Esta multiplicación es seguida, en cada una de las ramas, por un filtrado de paso bajo H_{b}(z).
Se obtienen así unas muestras:
Z_{k}{}^{+} = [rk + 1/2 \cdot c|k|_{L}] \otimes h^{b}{}_{k}
(3-7)Z_{k}{}^{-} = [rk - 1/2 \cdot c|k|_{L}] \otimes h^{b}{}_{k}
h^{b}_{k} designa la respuesta impulsional del filtro de paso bajo. La pasabanda de dicho filtro está limitada en la práctica hacia abajo por:
\ding{226} el ritmo f_{s} (f_{s}=f_{c}/M) de los símbolos de datos, en cuya ausencia se pierde energía útil, o
\newpage
\ding{226} la dinámica entre el transmisor y el receptor; se supone que la distancia entre el transmisor y el receptor es constante frente a la pasabanda de H^{b}(z).
La señal de error e_{k}se genera de la siguiente manera:
(3-8)e_{k} = |Z^{-}{}_{k}|^{2} - |Z^{+}{}_{k}|^{2}
Para obtener una señal de error directamente explotable por el circuito operativo NCO, se filtra generalmente la señal e_{k} mediante otro filtro digital, el filtro de bucle, cuya función de transferencia es H^{d}(z). Dado que e_{k} es independiente de los símbolos de datos, las características de H^{d}(z) se determinan principalmente mediante la respuesta deseada del bucle DLL a la dinámica entre el transmisor y el receptor, debiendo el bucle de estimación de fase ser capaz de seguir una distancia linear creciente o decreciente entre el transmisor y el receptor sin error residual de seguimiento.
La característica n(\varepsilon) indica la dependencia de la señal de error ek en función del error de fase de código \varepsilon.
(3-9)n(\varepsilon) = E[e_{k}|\varepsilon_{k} = \varepsilon \forall k]
E[\bullet] designa la probabilidad
Tenemos:
(3-10)n(\varepsilon) = g^{2}[(\varepsilon - 0,5)T_{c}] - g^{2} [(\varepsilon + 0,5)T_{c}]
con una codificación SRC de factor de atenuación \beta, g(\varepsilonTc) que vale:
g(\varepsilon T_{c}) = sen(\Pi \varepsilon) / \pi \varepsilon \cdot cos(\phi \beta \varepsilon) / 1 - (2 \beta \varepsilon)^{2}
La función de autocorrelación g(\varepsilonT_{c}) donde (3=0,35 se muestra en la figura 3a. Las muestras adelantadas E y retrasadas L según (3,8) se indican para \varepsilon=0 (sin error de fase de código).
Según la fórmula (3,10), la curva en S resultante se ilustra en la figura 3a. Mediante esta curva, el bucle DLL comanda el interpolador, de manera que la señal recibida es alineada con la señal generada localmente (es decir, \varepsilon=0).
El discriminador de la invención (figura 4) emplea dos muestras adelantadas E1 (Z^{-}_{k}) y E2 (Z^{-}_{k}) en instantes (\varepsilon-t_{1})T_{c} y (\varepsilon-t_{2})T_{c}, que son generados según la fórmula:
ek = K_{\beta} \cdot \Re (Zk - - / Zk- + S_{\beta})
\bullet
\Re (\bullet) designa la parte real
\bullet
K_{\beta} es una constante
\bullet
S_{\beta} es un desfase
\bullet
Las dos muestras E2 y E1 se definen (véase (3,7)) mediante:
Z^{-}{}_{k} = _{L} r_{k - t2} \cdot C_{P,|K|L}\rfloor \otimes hb/k
\hskip0.5cm
y
\hskip0.5cm
Z^{-}{}_{k} = _{L} r_{k - t1} \cdot C_{P,|K|L}\rfloor \otimes hb/k
donde, por ejemplo, t_{1}=1,5 y t_{2}=0,5
También se pueden generar estas dos muestras desfasando la réplica local del código según la fórmula:
Z^{-}{}_{k} = [rk \cdot C_{P,|(k-t2)|L}] \otimes hb/k
y
Z^{--}{}_{k} = [rk \cdot C_{P,|(k-t1)|L}] \otimes hb/k
\newpage
Se debe elegir el desfase para que el valor esperado de la señal de error E[e_{k}] sea nulo cuando el error de la fase del código \varepsilon_{k} es nulo.
Aplicando la fórmula (3-12), resulta:
(3-13)E[\Re (Z^{--}{}_{K} / Z^{-}{}_{K} + S_{\beta})|\varepsilon k=0]\equiv 0 \Rightarrow S_{\beta} = E[-\Re(Z^{--}{}_{K} / Z^{-}{}_{K})|\varepsilon k=0] = -g(-t1Tc)/g(-t2Tc)
El factor de pendiente K\beta se elige preferiblemente para que el valor de la pendiente
E [d / d\varepsilon \ e_{k} |\varepsilon_{k} = 0] \ sea \ igual \ a \ 1 \ cuando \varepsilon_{K} = 0
De lo que resulta:
(3-14)K_{\beta} = 1 / E[d/d\varepsilon\Re(Z^{--}{}_{k}/Z^{-}{}_{k} + S\beta)|\varepsilon k=0] = 1/ d/d\varepsilon(g((\varepsilon - t1)T_{c}) /g((\varepsilon - t2)T_{c}))\varepsilon = 0
donde, por ejemplo, t_{1}=1,5 y t_{2}=0,5.
Como se muestra en la fórmula (3-12), el discriminador ek de la invención es independiente de la fase \Pi de la portadora. Esta dependencia se suprime gracias al hecho de que se genera la relación S_{K}^{-}/Z_{K}^{-}.
La función de autocorrelación g(t) para \beta = 0,35 se representa en la figura 3b.
Se puede utilizar el módulo discriminador de dos maneras:
-
puede sustituir directamente un discriminador conocido,
-
puede integrarse en un discriminador conocido para suministrar señales de corrección de bucle abierto para reducir el error multitrayectos inducido.
La figura 4 representa el bloque diagrama de una rama DLL que incluye un módulo discriminador de la invención.
Comparando las figuras 2 y 4, se observa que:
La rama que genera Z_{k} en la figura 2 es sustituida por una rama que genera Z^{-}_{k}con un elemento de retraso de una cantidad 2T_{c};
El cálculo de la señal de error e_{K} en función de Z_{k} se realiza en el circuito DISCO mediante la fórmula 3-12 correspondiente al módulo de la invención.
En el caso de una señal del tipo SRC, la curva en S n(\varepsilon) se obtiene mediante:
(3-15)n(\varepsilon) = K\beta \Re (g((\varepsilon - t1)T_{c}) / g((\varepsilon - t2)T_{c}) + S\beta)
donde
g(\varepsilon T_{c}) = sen(\pi \varepsilon) / \pi \varepsilon \cdot cos(\pi \beta \varepsilon) /1 - (2\beta \varepsilon)^{2}
donde, por ejemplo, t1=1,5 y t2=0,5.
El módulo discriminador corresponde con una buena aproximadamente al comportamiento buscado, es decir, que la salida de la curva en S es proporcional a la entrada e (e=K.\varepsilon) para -0,5\varepsilon \leq 0,5.
Eligiendo K_{\beta} según la fórmula, se obtiene 3-14, K = 1 y e = \varepsilon.
Se puede utilizar el módulo discriminador para una estimación en bucle abierto para corregir la salida de la fase de código de un módulo conocido, como se ilustra en la figura 4.
Con relación al módulo de la figura 4, se dispone una rama adicional con un elemento de retraso 2Tc para generar la señal Z^{-}_{k}.
Se puede emplear un filtro digital H^{d}_{com}(z) como filtro de paso bajo de la salida del nuevo módulo discriminador. El funcionamiento de la rama DLL no cambia con relación al caso de la figura 2.
La fase del código se corrige mediante la salida del nuevo módulo discriminador que es alimentada por medio del filtro de paso bajo H^{d}_{com}(z). Dado que el nuevo módulo está menos influenciado por los trayectos múltiples, el desfase contenido en la fase estimada del código de la rama DLL puede corregirse en gran parte.
Como variante, la fórmula puede sustituirse por una función de amplitud:
e_{k} = K_{\beta}(\lambda|Z^{--}{}_{k}| / |Z^{-}{}_{k}| + S_{\beta})
siendo \lambda un número no nulo incluido entre -1 y +1.
Según otra variante, la señal de error e_{k} puede determinarse mediante:
ek = K\beta \cdot \Re(2.z^{--}{}_{k} / z^{-}{}_{k} + z^{+}{}_{k} + S_{\beta})
Esta expresión se adapta especialmente a una corrección en rama abierta (figura 5). Debido al efecto de media de salidas de los filtros de paso bajo (Z^{-}_{k}+Z^{+}_{k})/2, la potencia del ruido en la variable resultante está, en efecto, dividida entre dos, lo que conduce a un menor ruido en la señal e_{k}.
Se observa que los instantes de muestreo t_{1}= 1,5 y t_{2} = 0,5 pueden tener otros valores, y que el desfase t_{1}-t_{2} entre dichos instantes de muestreo puede ser distinto de un bito o "chip".

Claims (7)

1. Módulo discriminador de trayectos múltiples para un sistema de comunicación y/o de navegación que emplea una modulación con espectro repartido, módulo que presenta una entrada capaz de recibir señales de navegación, un muestreador para suministrar señales muestreadas \Gamma_{k} a una frecuencia doble de la frecuencia aparente fc del código de dichas señales, así como un submódulo (DISCR) de cálculo de una señal de error e_{k} a partir de dichas señales muestreadas \Gamma_{k} y del késimo valor del código de reparto generado localmente, el cual es designado mediante CIKIL, caracterizado porque el submódulo (DISCR) emplea el cálculo de una relación entre valores de correlación que incluye dos valores de correlación (Z^{-}_{K}, Z^{-}_{K}) adelantados con relación al valor de fase real, procediendo dichos valores de una correlación entre unas muestras de la señal recibida y dicha señal de código de reparto generada localmente, según una de las siguientes fórmulas:
e_{k} = K_{\beta} \cdot \Re (Z_{k} - - / Z_{k} - + S_{\beta})
\vskip1.000000\baselineskip
e_{k} = K_{\beta} (\lambda Z^{--}_{k} / Z^{-}_{k} + S_{\beta})
siendo \lambda no nulo e incluido entre -1 y +1
e_{k} = K_{\beta} \Re (2 \ Z^{--}{}_{k} / Z^{-}{}_{k} + Z^{+}{}_{k} + S_{\beta})
\vskip1.000000\baselineskip
donde \{Z^{+}{}_{K} = [r_{k + 0.5} \cdot C|K|L] \otimes \ h \ b/k
\hskip1.3cm
Z^{-}{}_{K} = [r \ k - t2 \cdot C|K|L] \otimes \ h \ b/k
\hskip1.3cm
Z^{--}{}_{K} = [r \ k - t1 \cdot C|K|L] \otimes \ h \ b/k
\vskip1.000000\baselineskip
\vskip1.000000\baselineskip
\hskip1cm
o \{Z^{+}{}_{K} = [rk \cdot C |k + 1 / 2|L] \otimes h \ b/k
\hskip1cm
Z^{-}{}_{K} = [rk \cdot C (k-t2)L] \otimes \ h \ b/k
\hskip1cm
Z^{--}{}_{K} = [rk \cdot C(k-t1)L] \otimes \ h \ b/k
K_{\beta} = constante
h^{b}_{k} designa la frecuencia impulsional de un filtro de paso bajo
\Re designa la parte real
S_{\beta} = - g(-t1Tc) / g(-t2Tc)
donde
g (aTc) = Sen \Pi a \ Cos\Pi \beta a / \pi a[1 - (2\beta a)]
\beta designa el factor de atenuación de la señal SRC.
2. Módulo discriminador, según la reivindicación 1, caracterizado porque
K\beta = 1 / d/d\varepsilon (g((\varepsilon - t1)Tc) / g((\varepsilon - t2)Tc))_{\varepsilon=0}
\varepsilon designa el error de fase del código.
3. Módulo discriminador, según una de las reivindicaciones anteriores, caracterizado porque t1 - t2 = 1.
4. Módulo discriminador, según la reivindicación 3, caracterizado porque t1 = 1,5 y t2 = 0,5.
5. Módulo discriminador, según una de las reivindicaciones anteriores, caracterizado porque \lambda = -1.
6. Módulo discriminador, según una de las reivindicaciones anteriores, caracterizado porque \lambda = 1.
7. Sistema de navegación que presenta un módulo discriminador que genera una señal de error e'k que asegura, de manera conocida en sí, una corrección en bucle cerrado a partir de señales muestreadas Z^{+}_{K} y Z^{-}_{K}, caracterizado porque lleva asociado un módulo discriminador, según una de las reivindicaciones anteriores, para generar en bucle abierto una señal de corrección de la salida de fase del código.
ES02748961T 2001-06-25 2002-06-19 Modulo discriminador de trayectos multiples para un sistema de navegacion. Expired - Lifetime ES2250675T3 (es)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0108363 2001-06-25
FR0108363A FR2826528B1 (fr) 2001-06-25 2001-06-25 Module discriminateur de trajets multiples pour un systeme de navigation

Publications (1)

Publication Number Publication Date
ES2250675T3 true ES2250675T3 (es) 2006-04-16

Family

ID=8864743

Family Applications (1)

Application Number Title Priority Date Filing Date
ES02748961T Expired - Lifetime ES2250675T3 (es) 2001-06-25 2002-06-19 Modulo discriminador de trayectos multiples para un sistema de navegacion.

Country Status (8)

Country Link
US (1) US8050309B2 (es)
EP (1) EP1400028B1 (es)
JP (1) JP4133808B2 (es)
CA (1) CA2450940C (es)
DE (1) DE60206501T2 (es)
ES (1) ES2250675T3 (es)
FR (1) FR2826528B1 (es)
WO (1) WO2003003603A1 (es)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7995683B2 (en) * 2007-10-24 2011-08-09 Sirf Technology Inc. Noise floor independent delay-locked loop discriminator
JP5104295B2 (ja) * 2007-12-26 2012-12-19 セイコーエプソン株式会社 コード位相算出方法、プログラム及びコード位相算出装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19712751A1 (de) * 1997-03-26 1998-10-08 Deutsch Zentr Luft & Raumfahrt Empfänger zum Empfangen von Signalen eines Satellitennavigationssystems
US6735242B1 (en) * 1999-08-30 2004-05-11 Nokia Corporation Time tracking loop for pilot aided direct sequence spread spectrum systems
US6925108B1 (en) * 2000-05-26 2005-08-02 Freescale Semiconductor, Inc. Ultrawide bandwidth system and method for fast synchronization
US7010073B2 (en) * 2001-01-19 2006-03-07 Qualcomm, Incorporated Delay lock loops for wireless communication systems

Also Published As

Publication number Publication date
FR2826528A1 (fr) 2002-12-27
EP1400028B1 (fr) 2005-10-05
US8050309B2 (en) 2011-11-01
JP2004531178A (ja) 2004-10-07
CA2450940A1 (en) 2003-01-09
WO2003003603A1 (fr) 2003-01-09
DE60206501T2 (de) 2006-07-27
JP4133808B2 (ja) 2008-08-13
DE60206501D1 (de) 2005-11-10
FR2826528B1 (fr) 2003-10-31
US20040240530A1 (en) 2004-12-02
CA2450940C (en) 2012-05-29
EP1400028A1 (fr) 2004-03-24

Similar Documents

Publication Publication Date Title
US6282231B1 (en) Strong signal cancellation to enhance processing of weak spread spectrum signal
US8170085B2 (en) Multipath error estimation in satellite navigation receivers
ES2292048T3 (es) Metodo y aparato para la adquisicion y seguimiento de una señal portadora de desplazamiento binario (boc).
ES2298584T3 (es) Arquitectura de hardware para procesar señales portadoras de desviacion binaria alternativa galileo (altboc).
US8000381B2 (en) Unbiased code phase discriminator
CA2568553C (en) Improved discriminator function for gps code alignment
ES2305824T3 (es) Metodo y dispositivo para demodular señales portadoras de desplazamiento binario alterno de galileo (altboc).
AU724006B2 (en) Mitigation of multipath effects in global positioning system receivers
US6414987B1 (en) Code multipath estimation for weighted or modified tracking
US8000378B2 (en) Narrow correlator technique for multipath mitigation
JPH10505664A (ja) Pコードで変調された信号の再生及び追従のための全地球測位システム(gps)受信器
Yang et al. Acquisition of the GPS L5 signal using coherent combining of I5 and Q5
JP2009535640A (ja) 衛星ナビゲーション受信機用適応符号発生器
Avila-Rodriguez et al. CBOC: an implementation of MBOC
KR101133398B1 (ko) 무선 내비게이션 수신기 내의 멀티패스 완화 방법
US6236355B1 (en) Fast recovery of lock point information for GPS RTK receiver
EP2188907B1 (en) Method and device for multipath mitigation
ES2250675T3 (es) Modulo discriminador de trayectos multiples para un sistema de navegacion.
US6345068B1 (en) Hierarchical delay lock loop code tracking system with multipath correction
US6687316B1 (en) High resolution correlator technique for spread spectrum ranging system code and carrier multipath mitigation
US7298776B2 (en) Acquisition of a gated pilot signal with coherent and noncoherent integration
Julien et al. 1-bit processing of composite BOC (CBOC) signals
US7010021B2 (en) Method of generating a time shifted signal
ES2322973T3 (es) Sistema y metodo correlador altamente paralelo para gps.
Tavasoli Applications of spread spectrum communication in global positioning system