ES2340672T3 - Metodo y sistemas de desconvolucion para periodos irregulares de repeticion de impulsos, tratamiento doppler y de ecos parasitos. - Google Patents

Metodo y sistemas de desconvolucion para periodos irregulares de repeticion de impulsos, tratamiento doppler y de ecos parasitos. Download PDF

Info

Publication number
ES2340672T3
ES2340672T3 ES03101597T ES03101597T ES2340672T3 ES 2340672 T3 ES2340672 T3 ES 2340672T3 ES 03101597 T ES03101597 T ES 03101597T ES 03101597 T ES03101597 T ES 03101597T ES 2340672 T3 ES2340672 T3 ES 2340672T3
Authority
ES
Spain
Prior art keywords
deconvolution
spectra
dft
echoes
clutter
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
ES03101597T
Other languages
English (en)
Inventor
Radmila Erkocevic-Pribic
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.)
Thales Nederland BV
Original Assignee
Thales Nederland BV
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 Thales Nederland BV filed Critical Thales Nederland BV
Application granted granted Critical
Publication of ES2340672T3 publication Critical patent/ES2340672T3/es
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/414Discriminating targets with respect to background clutter
    • 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/08Systems for measuring distance only
    • G01S13/10Systems for measuring distance only using transmission of interrupted, pulse modulated waves
    • G01S13/22Systems for measuring distance only using transmission of interrupted, pulse modulated waves using irregular pulse repetition frequency
    • G01S13/225Systems for measuring distance only using transmission of interrupted, pulse modulated waves using irregular pulse repetition frequency with cyclic repetition of a non-uniform pulse sequence, e.g. staggered PRF
    • 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/52Discriminating between fixed and moving objects or between objects moving at different speeds
    • G01S13/522Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves
    • G01S13/524Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/288Coherent receivers

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Polysaccharides And Polysaccharide Derivatives (AREA)
  • Saccharide Compounds (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Detergent Compositions (AREA)

Abstract

Método de desconvolución de una señal muestreada en períodos irregulares de repetición de impulsos, x(tm), que comprende las siguientes operaciones: - (S1) conversión de las muestras x(tm) irregulares en muestras r(iTε) regulares; - (S2) cálculo del espectro dft(r) de estas muestras regulares; - (S3) aislamiento de los espectros de los ecos parásitos dft(r) suponiendo que los ecos parásitos se extienden por más de unos pocos intervalos de alcance; - (S4) estimación de las líneas espectrales de los ecos parásitos a partir de la media y de la anchura de los espectros aislados de los ecos parásitos; - (S5) resta de los espectros de los ecos parásitos estimados del espectro total dft(r); - (S6) desconvolución de los espectros restantes; en el que, en las operaciones de resta (S5) y de desconvolución (S6), los cálculos se ajustan al tipo de ecos parásitos mediante los espectros de los ecos parásitos dados por la operación (S3).

Description

Método y sistemas de desconvolución para períodos irregulares de repetición de impulsos, tratamiento Doppler y de ecos parásitos.
Este invento se refiere al tratamiento de señales de radar. En particular, este invento está relacionado con el tratamiento Doppler y el filtrado de ecos parásitos en una señal muestreada en períodos irregulares de repetición de impulsos (PRT).
\vskip1.000000\baselineskip
Descripción de la técnica relacionada
En un radar usual, el tratamiento Doppler hace uso de una transformada discreta de Fourier (DFT) sobre señales muestreadas regularmente porque el período de repetición de impulsos es constante dentro de un tren de impulsos. Por ello, dicho tratamiento usual que emplea la transformada discreta de Fourier (DFT) no es aplicable en el caso de una señal muestreada en períodos irregulares de repetición de impulsos (PRT).
La expresión "período irregular de repetición de impulsos" (PRT) se utiliza, usualmente, para impedir que a velocidades ciegas (ambigüedad de la velocidad), los inhibidores bloqueen el período de repetición de impulsos (PRT) y, también, para resolver ambigüedades Doppler.
Los actuales trabajos sobre muestreo irregular están motivados, en su mayor parte, por el tratamiento de la imagen. Se han desarrollado algoritmos eficientes basados en operaciones de tramas como en el trabajo "Métodos numéricos eficientes en la teoría de muestreo no uniforme", escrito por Feichtinger, H.G. y otros en Numerische Mathematik 69 (NUHAG). La transformada en óndulas en conjuntos de filtros digitales como en "Sistemas multifrecuencia y conjuntos de filtro", escrito por Vaidyanathan, P.P., Prentice-Hall Publishers (1993), ilustra eficientes aplicaciones 1D con muestras irregulares obtenidas mediante decimación en enteros o expansión de muestras uniformes.
En el tratamiento Doppler de señales de radar, se ha estudiado el escalonamiento del período de repetición de impulsos (PRT), pero haciendo énfasis en el análisis de Fourier, en vez de en resolver al problema del muestreo irregular.
Una excepción la constituye el algoritmo del NSSL propuesto en la patente norteamericana 6.081.221 relacionada con el radar Doppler atmosférico. Este algoritmo del NSSL ofrece una solución no ambigua en la velocidad Doppler y solamente realiza filtrado de ecos parásitos terrestres. Tales ecos parásitos terrestres tienen una frecuencia media igual a cero, por cuanto no se desplazan. De este modo, los ecos parásitos que tengan una frecuencia media como los ecos parásitos marinos o los debidos a la lluvia..., no pueden ser filtrados con este algoritmo propuesto por el NSSL.
El documento "filtrado de ecos parásitos y estimación del momento espectral para radares Doppler atmosféricos utilizando períodos escalonados de repetición de impulsos (PRT)", publicado en Marzo de 2000 por M. Sachidananda, parece ser la técnica anterior más próxima.
\vskip1.000000\baselineskip
Sumario del invento
Este invento resuelve los inconvenientes antes mencionados, en particular los relacionados con la ambigüedad de la velocidad y el filtrado de cualquier tipo de eco parásito, proporcionando un método de desconvolución que filtra cualquier clase de eco parásito, incluso ecos parásitos variables como los originados por el mar, la lluvia,...
Los objetos y realizaciones del presente invento se describen en lo que sigue en las reivindicaciones 1 a 11.
\vskip1.000000\baselineskip
Breve descripción de los dibujos
Otras características y ventajas del invento resultarán evidentes a partir de la siguiente descripción de ejemplos de realizaciones del invento con referencia a los dibujos, que muestran detalles esenciales para el invento, y a partir de las reivindicaciones. Los detalles individuales pueden incorporarse en una realización del invento por separado o conjuntamente formando cualquier combinación.
- Las figuras 1a, 1b, 1c y 1d muestran espectros obtenidos a partir de una señal muestreada en períodos irregulares de repetición de impulsos (PRT) mediante cuatro técnicas diferentes de acuerdo con la tecnología anterior; los métodos NSSL, NUHAG, de óndulas y LSP, respectivamente.
- Las figuras 2a y 2b representan ecos parásitos de, respectivamente, la tierra y el mar, obtenidos por el NSSL.
- La figura 3 ilustra una gráfica de proceso de pasos significativos del método para filtrar ecos parásitos de acuerdo con el invento.
- Las figuras 4a, 4b, 4c, 4d y 4e ilustran el método de acuerdo con el invento, respectivamente, mediante espectros del muestreo irregular utilizado, la señal simulada, los ecos parásitos marinos simulados, la entrada al filtro de ecos parásitos y la salida, para una señal exponencial compleja libre de ruido, a la frecuencia de 1,2 y un eco parásito marino de configuración Gaussiana con amplitud media y un ancho de banda correspondiente a un estado 5 del mar cuando PRF y RF fuesen, respectivamente, 1 kHz y 1 GHz.
Las figuras 5a, 5b, 5c, 5d y 5e ilustran el método de acuerdo con el invento aplicado al eco de blancos y a ecos parásitos marinos de la figura 4 con ruido Gaussiano añadido y amplitud con distribución de Weibull, respectivamente, mediante espectros del muestreo irregular utilizado, la señal simulada, los ecos parásitos marinos simulados, la entrada al filtro de ecos parásitos y la salida para una señal exponencial compleja a la frecuencia de 1,28, un eco parásito marino con configuración Gaussiana con amplitud media y ancho de banda correspondiente a un estado 5 del mar, cuando PRF y RF fuesen, respectivamente, 1 kHz y 1 GHz y una frecuencia de blanco de 1,28.
\vskip1.000000\baselineskip
Descripción más detallada
El muestreo irregular no ha tenido un uso tan extendido como el muestreo uniforme debido a que los análisis de tiempo y de frecuencia de las muestras irregulares tienen una implicación bastante elevada.
El muestreo aleatorio implica, usualmente, números aleatorios añadidos a intervalos regulares de muestreo. El espectro de energía basado en la DFT de muestras aleatorias consiste en el espectro de energía de la señal más ruido aditivo no correlacionado. Por ejemplo, incluso con una señal libre de ruido muestreada a la frecuencia de muestreo de Nyquist, la relación entre señal de salida y ruido podría ser de, solamente, 1.
El muestreo irregular determinístico supone repetir periódicamente la misma secuencia de K intervalos irregulares con el intervalo de muestreo medio T_{S}. Usualmente, se dice que el muestreo es entrelazado cuando K=2 y multifrecuencia o en racimo cuando K tiene una longitud arbitraria. Como un conjunto de muestreo de esta clase contiene K conjuntos de muestreo regulares, la DFT da K picos (para una componente de frecuencia de una señal) dentro del intervalo limitado por la frecuencia de muestreo 1/ T_{S}.
Para cualquier conjunto de muestreo regular {nT}, tal que T\leq1/(2B), una señal real x(t) limitada a una banda B de frecuencia, se escribe como x(t) = \sum_{n}x(nT).sen c[2B(t-nT)]. Las muestras irregulares se comprenden perfectamente en teoría, pero su método resulta ser, usualmente, demasiado complicado. Para un conjunto de muestreo irregular {t_{n}}, la reconstrucción con bases bi-ortogonales {sen c[2B(t-nT)]} y {\Psi_{n}(t)}, se escribe como: x(t) = \sum_{n} c_{n}.sen c[2B(t- t_{n})] = \sum_{n}x(t_{n}) \Psi_{n}(t), donde c_{n} es el producto interno de x(t) con {\Psi_{n}(t)}. Si el conjunto {t_{n}} está limitado como |t-t_{n}|<1/(8B), \Psi_{n}(t) es una función de interpolación de Lagrange. Las funciones base también pueden ser tramas, es decir, bases cuya ortogonalidad no se requiera. Las condiciones de trama son mucho más débiles y más útiles con fines
prácticos.
Las siguientes técnicas pertenecientes a la tecnología anterior ilustran el estado de la técnica en relación con el análisis espectral de muestras irregulares.
El periodograma de Lomb-Scargle (LSP) es el clásico periodograma basado en la DFT corregido preservando el comportamiento estadístico y la invariancia por traslación en el tiempo. Si las muestras son irregulares, el periodograma de Lomb-Scargle (LSP) sobre muestras multifrecuencia revela el patrón determinístico como se muestra en la figura 1d.
La transformada en óndulas de x(t) viene dada como X(\alpha,t) con una óndula desplazada en \tau y escalada en \alpha. El espectro se representa mediante |X(\alpha,t)|^{2}, denominado escalograma. Los corrientes algoritmos eficientes de transformada en óndulas soportan el muestreo uniforme con ajuste de escala diádico y desplazamiento (\alpha=2^{m} y \tau=n2^{-m}). En este estudio, cualquier muestreo irregular se lleva a cabo en una óndula Morlet existente. La transformada en óndulas es más robusta, pero su resolución de frecuencia es muy mala, como se muestra en la figura 1c. En general, la transformada en óndulas no es aplicable debido a que los ecos de radar son señales periódicas y, así, las óndulas difícilmente pueden servir como sus funciones base convenientes.
En la descomposición de tramas NUHAG, primero se re-escribe una señal de banda limitada, compleja, x(t), de energía finita, dada por M muestras irregulares, como polinomios trigonométricos p(t) de período 1 y grado K, siendo K<M/2. La solución para los coeficientes de Fourier a_{k}, |k|\leqK, x(\tau)=\sum_{k}a_{k}.exp(j2\pif_{k}.\tau), |f_{k}|\leqB, se basa en las propiedades del operador de tramas: \sum_{m}p(t_{m}).D_{K}(t-t_{m}) donde D_{K}(t)= \sum_{k}exp(j2\pik.t), representa la trama. El método NUHAG ilustrado por la figura 1b es aplicable para cualquier tipo de muestreo pero solamente hasta la frecuencia de muestreo lo que, desafortunadamente, no es de aplicación en el campo del radar.
La desconvolución de magnitud del NSSL se utiliza para el tratamiento mediante radar Doppler terrestre de un esquema de muestreo entrelazado. En general, un conjunto de muestreo NSSL {t_{m}} es multifrecuencia con una tasa K y el intervalo medio T_{K}. Un intervalo de tiempo NSSL (t_{m+1} - t_{m}), así como toda la secuencia KT_{S} son múltiplos enteros del mayor intervalo de tiempo común T_{\varepsilon}, KT_{S} =LT_{\varepsilon}, de modo que el conjunto regular más pequeño, {iT_{\varepsilon}} puede contener {t_{m}}, {t_{m}}\subset{iT_{\varepsilon}}.
Así, las muestras irregulares x(t_{m}) son convertidas en muestras rellenas con ceros r(iT_{\varepsilon}) que son producto de un esquema de muestreo c_{i}, c_{i}=\delta(iT_{\varepsilon}-t_{m}), y las muestras regulares x(iT_{\varepsilon}). Basándose en esta relación: r=diag(c).x (en forma vectorial), el espectro de x puede derivarse como sigue:
1
donde C es una matriz de Toeplitz cuyos vectores de fila son desplazados cíclicamente dft(c). Como C es singular y, por ello, no es invertible, la idea del NSSL es utilizar las magnitudes en su lugar. Solamente es aplicable si no existen adiciones complejas en el producto C.dft(x). Esta condición implica el ancho de banda x, pero ésta no es una restricción para la mayoría de los radares.
El espectro |dft(r)| contiene L réplicas del espectro |dft(x)| que son ponderadas por un coeficiente a partir de |dft(c)|, en un intervalo de frecuencias que es L/K veces más amplio que la frecuencia de muestreo media f_{s}. La desconvolución proporciona la réplica más intensa, es decir, el espectro de señal |dft(x)|.
Como se muestra en la figura 1a, el método del NSSL es comparable con el método NUHAG (figura 1b) para muestras multifrecuencia. Además, el método del NSSL permite resolver los problemas del radar de ambigüedad Doppler y filtrado de ecos parásitos terrestres. Desafortunadamente, el método del NSSL solamente es aplicable en el campo de los radares terrestres.
El método del NSSL solamente soporta frecuencias superiores a las frecuencias de Nyquist, a saber, hasta L/K veces la frecuencia de muestreo.
En los radares de seguimiento, en los que los Dopplers esperados son razonablemente conocidos, pueden bastar los espectros de señales soportados por el NSSL. A saber, el uso de las magnitudes en
2
implica que el espectro de la señal x no puede ser más amplio que N/L líneas espectrales, donde N y L son el número de muestras regulares (es decir, la longitud de los vectores r y x) y la longitud de la secuencia irregular básica, respectivamente. También puede funcionar en los casos en que el espectro de la señal sea más amplio que N/L, pero la distancia entre las líneas espectrales no puede ser un múltiplo entero de N/L.
El intervalo de frecuencias no ambiguas puede expandirse en forma ilimitada, pero el filtrado de los ecos parásitos interviene más al aumentar la complejidad del muestreo.
De acuerdo con el invento, si se conoce el ancho de banda B_{cl} de los ecos parásitos, pueden estimarse las amplitudes de los ecos parásitos complejos y filtrarse de |dft(r)| antes de la desconvolución. Así, el filtro de acuerdo con el invento está destinado al tipo de ecos parásitos particular del entorno en el que se transmite la señal, por ejemplo, ecos parásitos marinos como se muestra en las figuras 2 y 3.
Puede estimarse cualquier clase de ecos parásitos. Luego, se centra un conjunto {i} en un entero correspondiente a la frecuencia media Doppler como se muestra en las figuras 3a y 3b para, respectivamente, ecos parásitos terrestres y ecos parásitos marinos.
En general, las amplitudes g_{cl}(i, k) de la componente espectral de orden i de los ecos parásitos, f_{i} \in B_{cl}, i \in [-I, J] repetidas en la línea distinta de cero de orden k del espectro de código, k = 1, L, puede estimarse como sigue:
3
donde los vectores d y a_{L} contienen L componentes espectrales distintas de cero de dft(c), d_{1}={dft(c)}_{1+(l-1)N/L}, |dft(c)|=1, a_{L}\equiv[d_{2}d_{3}...d_{L}d_{1}] y z=dft(r).
\newpage
En general, el intervalo entero i \in [-I, J] es desconocido. Para estimar las líneas espectrales de ecos parásitos, f_{i} \in B_{cl} utilizando la correlación, se supone que los ecos parásitos tienen una configuración Gaussiana y que las señales recibidas son exponenciales complejas.
La configuración Gaussiana implica que la anchura del espectro, \sigma_{cl}, puede estimarse a partir de los coeficientes de correlación \rho(\tau) en los dos primeros retardos, como sigue:
4
\vskip1.000000\baselineskip
El espectro \sigma_{cl} de los ecos parásitos de tres veces \sigma_{f} desde la frecuencia media hacia la izquierda y hacia la derecha implica el 99,7% del espectro. Un \sigma_{cl} óptimo es un compromiso entre tomar el máximo posible de los espectros de los ecos parásitos y dejar todo lo posible de los espectros de señales y, como tal, debe tratarse como adaptable a un eco de radar particular.
La frecuencia Doppler media \overline{f_{d}} puede estimarse a partir del argumento de la autocorrelación en el primer retardo r(\tau),r(\tau)=\rho(\tau)exp(j\omega_{d}\tau), como sigue:
5
\vskip1.000000\baselineskip
El método de desconvolución de acuerdo con el invento se ilustra en la figura 3. Este método puede resumirse en las siguientes operaciones:
[S1: operación de conversión] Las muestras x(t_{m}) irregulares se convierten en muestras r(iT_{\varepsilon})regulares;
[S2: operación de DFT] Se calcula el espectro dft(r) de estas muestras regulares;
[S3: operación de aislamiento] Se aíslan los espectros de los ecos parásitos suponiendo que éstos se extienden por más de unos pocos intervalos de alcance;
[S4: estimación] Se estiman las líneas espectrales de los ecos parásitos a partir de la media y de la anchura de los espectros aislados de los ecos parásitos;
[S5: resta] Se restan los espectros de los ecos parásitos del espectro total dft(r);
[S6: desconvolución] Se realiza la desconvolución de los espectros restantes.
En cuyas operaciones S5 y S6 se ajustan los cálculos al tipo de ecos parásitos mediante los espectros de ecos parásitos obtenidos en la operación S3.
En una segunda realización, también puede calcularse el espectro dft(c) y sus L componentes distintas de cero hallados en la operación S2. En una tercera realización, en la operación S4 pueden estimarse las amplitudes de los espectros de los ecos parásitos. Así, pueden reducirse las operaciones de convolución, basándose en las L componentes espectrales distintas de cero entre las operaciones S5 y S6. En otras realizaciones, pueden combinarse estas realizaciones segunda y tercera.
Además, las líneas espectrales de los ecos parásitos pueden estimarse a partir de la media en \overline{f_{d}} = arg[r(\tau_{1})]/(2\pi\tau_{1}) y/o la anchura en
6
\newpage
de los espectros aislados [S4]. Además, las amplitudes de los espectros de los ecos parásitos pueden estimarse mediante
7
\vskip1.000000\baselineskip
La operación de convolución puede reducirse en
8
antes de que serealice la desconvolución [S6] de los espectros restantes en
9
El método de desconvolución de acuerdo con el invento se ilustra en las figuras 4a, 4b, 4c, 4d y 4e cuando el muestreo multifrecuencia contiene cinco componentes espectrales distintas de cero como se muestra en la figura 4a, que modulan la señal de entrada de la figura 4d que contiene el eco del blanco libre de ruidos de la figura 4b y los ecos parásitos marinos no aleatorios de la figura 4c. Todos estos espectros se dan en dB en función de la frecuencia Doppler normalizada. Se supone que la densidad espectral de la energía (psd) de los ecos parásitos marinos tiene configuración Gaussiana. El filtrado de los ecos parásitos y la desconvolución tienen como resultado el contenido de frecuencia del eco del blanco como se muestra en la figura 4e.
Este método de desconvolución, en cualquiera de sus realizaciones, se aplica en un sistema de radar mediante un sistema de desconvolución, que comprende:
-
Medios para convertir las muestras irregulares x(t_{m}) en muestras regulares r(iT_{\varepsilon});
-
Medios para calcular el espectro dft(r) de estas muestras regulares;
-
Medios para aislar los espectros de los ecos parásitos en dft(r) suponiendo que los ecos parásitos se extienden por más de unos pocos intervalos de alcance;
-
Medios para estimar las líneas espectrales de los ecos parásitos a partir de la media y la anchura de los espectros parásitos aislados;
-
Medios para restar los espectros de los ecos parásitos estimados del espectro total dft(r);
-
Medios para realizar la desconvolución de los espectros restantes.
En las figuras 5a, 5b, 5c, 5d y 5e, el método de desconvolución de acuerdo con el invento se aplica a un escenario más real. El ruido Gaussiano se suma al eco del blanco, como se muestra en la figura 5b, y los ecos parásitos marinos de la figura 5c se hacen aleatorios. Los ecos parásitos se distribuyen, normalmente, en una señal muestreada en un período de repetición de impulsos irregular (PRT) que tiene la correlación determinada por la densidad espectral de energía con configuración Gaussiana, y con una amplitud por alcance con una distribución, supuestamente, según Weibull. El mismo muestreo (figura 5a) modula la señal de entrada de la figura 5d que contiene un eco de blanco con ruidos (figura 5b) y ecos parásitos marinos aleatorios (figura 5c). Todos estos espectros se dan en dB en función de la frecuencia Doppler normalizada. El método de desconvolución de acuerdo con el invento funciona mejor si la frecuencia Doppler normalizada se ajusta a un valor de 1,28. Esto puede conseguirse cambiando la frecuencia de muestreo. La estimación y la resta de los ecos parásitos, seguidas por la desconvolución, dan como resultado, de nuevo, el contenido de frecuencia del eco del blanco, como se muestra en la figura 5e.
Pueden utilizarse otras secuencias irregulares que pudieran ser óptimas para la extensión del alcance Doppler, pero también para el filtrado de los ecos parásitos.
Una ventaja del método de desconvolución de acuerdo con el invento es que trabaja por encima de las frecuencias de Nyquist y que ofrece un método de filtrado para cualquier clase de ecos parásitos: cualesquiera ecos parásitos "variables" como los generados por el mar o por la lluvia, pero también ecos parásitos no variables, como los ecos parásitos terrestres.
Así, otra aplicación para el filtro y el método propuestos puede ser el filtrado de ecos parásitos de "lluvia" o, incluso, cualquier filtrado de ecos parásitos (para ecos parásitos "variables" y "no variables"). De este modo, permite el uso de este filtro de ecos parásitos en cualquier clase de radar y en cualquier entorno. Ello se debe a que el filtrado de los ecos parásitos de acuerdo con el invento es posible para ecos parásitos cuyas líneas espectrales no estén, necesariamente, centradas en cero.
Más generalmente, un sistema de desconvolución de esta clase puede utilizarse para realizar la desconvolución de cualquier tipo de señal muestreada irregular y no solamente una de radar.

Claims (11)

1. Método de desconvolución de una señal muestreada en períodos irregulares de repetición de impulsos, x(t_{m}), que comprende las siguientes operaciones:
-
(S1) conversión de las muestras x(t_{m}) irregulares en muestras r(iT_{\varepsilon}) regulares;
-
(S2) cálculo del espectro dft(r) de estas muestras regulares;
-
(S3) aislamiento de los espectros de los ecos parásitos dft(r) suponiendo que los ecos parásitos se extienden por más de unos pocos intervalos de alcance;
-
(S4) estimación de las líneas espectrales de los ecos parásitos a partir de la media y de la anchura de los espectros aislados de los ecos parásitos;
-
(S5) resta de los espectros de los ecos parásitos estimados del espectro total dft(r);
-
(S6) desconvolución de los espectros restantes;
en el que, en las operaciones de resta (S5) y de desconvolución (S6), los cálculos se ajustan al tipo de ecos parásitos mediante los espectros de los ecos parásitos dados por la operación (S3).
\vskip1.000000\baselineskip
2. Método de desconvolución de acuerdo con cualquiera de las reivindicaciones precedentes, en el que la operación (S2) comprende, también, las siguientes operaciones secundarias:
-
cálculo del espectro del esquema de muestreo, dft(c)
-
búsqueda de sus L componentes distintas de cero.
\vskip1.000000\baselineskip
3. Método de desconvolución de acuerdo con la reivindicación precedente, en el que los pasos de convolución en las operaciones de resta (S5) y de desconvolución (S6) se reducen basándose en las L componentes espectrales distintas de cero.
4. Método de desconvolución de acuerdo con cualquiera de las reivindicaciones precedentes, en el que la operación (S4) abarca también, las amplitudes de la estimación de los espectros de los ecos parásitos.
5. Sistema de desconvolución de una señal muestreada en períodos irregulares de repetición de impulsos, x(t_{m}), que comprende:
-
medios para convertir las muestras irregulares x(t_{m}) en muestras regulares r(iT_{\varepsilon});
-
medios para calcular el espectro dft(r) de estas muestras regulares;
-
medios para aislar los espectros de los ecos parásitos en dft(r) suponiendo que los ecos parásitos se extienden por más de unos pocos intervalos de alcance;
-
medios para estimar las líneas espectrales de los ecos parásitos a partir de la media y la anchura de los espectros de los ecos parásitos aislados;
-
medios para restar los espectros de los ecos parásitos estimados del espectro total dft(r);
-
medios para realizar la desconvolución de los espectros restantes;
en el que el cálculo de los medios para restar y de los medios para realizar la desconvolución, se ajustan al tipo de ecos parásitos mediante los espectros de los ecos parásitos dados por los medios para aislar.
\vskip1.000000\baselineskip
6. Sistema de desconvolución de acuerdo con cualquiera de las reivindicaciones 6 o 7, en el que los medios para calcular:
-
calculan el espectro dft(c) del esquema de muestreo, y
-
buscan sus L componentes distintas de cero.
\vskip1.000000\baselineskip
7. Sistema de desconvolución de acuerdo con la reivindicación precedente, en el que las operaciones de convolución de los medios para restar y de los medios para realizar la desconvolución, se reducen basándose en la L componentes espectrales distintas de cero.
8. Sistema de desconvolución de acuerdo con cualquiera de las reivindicaciones 6 a 8, en el que los medios para estimar también estiman las amplitudes de los espectros de los ecos parásitos.
9. Uso del método de desconvolución de acuerdo con la reivindicación 1 en un sistema de radar.
10. Uso del método de desconvolución de acuerdo con la reivindicación 1 como filtrado de ecos parásitos.
11. Uso del método de desconvolución de acuerdo con la reivindicación precedente como filtrado de ecos parásitos marinos.
ES03101597T 2002-06-05 2003-06-02 Metodo y sistemas de desconvolucion para periodos irregulares de repeticion de impulsos, tratamiento doppler y de ecos parasitos. Expired - Lifetime ES2340672T3 (es)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
NL1020763A NL1020763C2 (nl) 2002-06-05 2002-06-05 Werkwijze en -systeem voor het deconvolueren van een met een onregelmatige PRT bemonsterd signaal, en de toepassingsmogelijkheden ervan.
NL1020763 2002-06-05

Publications (1)

Publication Number Publication Date
ES2340672T3 true ES2340672T3 (es) 2010-06-08

Family

ID=29546434

Family Applications (1)

Application Number Title Priority Date Filing Date
ES03101597T Expired - Lifetime ES2340672T3 (es) 2002-06-05 2003-06-02 Metodo y sistemas de desconvolucion para periodos irregulares de repeticion de impulsos, tratamiento doppler y de ecos parasitos.

Country Status (11)

Country Link
US (1) US6803876B2 (es)
EP (1) EP1369703B1 (es)
JP (1) JP2004037457A (es)
KR (1) KR100966060B1 (es)
AT (1) ATE458204T1 (es)
CA (1) CA2431094C (es)
DE (1) DE60331277D1 (es)
ES (1) ES2340672T3 (es)
IL (1) IL156293A (es)
NL (1) NL1020763C2 (es)
ZA (1) ZA200304397B (es)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL1020763C2 (nl) * 2002-06-05 2003-12-08 Thales Nederland Bv Werkwijze en -systeem voor het deconvolueren van een met een onregelmatige PRT bemonsterd signaal, en de toepassingsmogelijkheden ervan.
US6937185B1 (en) * 2003-08-14 2005-08-30 Lockheed Martin Corporation Rain versus target discrimination for doppler radars
US6861974B1 (en) * 2003-10-16 2005-03-01 Lockheed Martin Corporation Clutter normalization by adaptation of transmit waveform
US7589666B2 (en) * 2004-12-30 2009-09-15 Vaisala, Inc. System and method for processing data in weather radar
JP2006250934A (ja) * 2005-03-08 2006-09-21 Synthesys Research Inc デューティ・サイクルを測定するための方法及び装置
US7583222B2 (en) 2005-08-26 2009-09-01 Vaisala Oyj Method for using pulse compression in weather radar
WO2012032333A1 (en) * 2010-09-08 2012-03-15 Bae Systems Plc Radar system and methods of controlling a radar system and the provision of a range-doppler target image signature
JP2013142775A (ja) * 2012-01-11 2013-07-22 Sony Corp 表示装置、電子機器、表示方法、並びにプログラム
CN102608597B (zh) * 2012-03-19 2014-07-23 西安电子科技大学 基于非完全数据解卷积的实孔径前视成像方法
FR3008191B1 (fr) * 2013-07-05 2015-07-17 Thales Sa Procede de detection d'au moins une cible par radar doppler a impulsions avec mesure non ambigue de la vitesse radiale et radar doppler a impulsions pour la mise en oeuvre d'un tel procede
EP3144695B1 (de) * 2015-09-17 2019-04-10 HENSOLDT Sensors GmbH Vorrichtung und verfahren zur messung von position und geschwindigkeit eines objekts
US11346929B2 (en) * 2017-01-16 2022-05-31 Koninklijke Philips N.V. Systems and methods for ultrafast ultrasound imaging
CN109001703A (zh) * 2018-08-10 2018-12-14 南京信息工程大学 一种基于小波包多阈值处理的海杂波去噪办法
CN109597060B (zh) * 2018-12-07 2021-06-25 北京敏视达雷达有限公司 一种雷达测速方法及装置
CN110146891A (zh) * 2019-01-17 2019-08-20 西安电子科技大学 基于盲区数据恢复的单通道参差合成孔径雷达成像方法
CN110161491B (zh) * 2019-06-28 2021-01-12 电子科技大学 一种针对微弱生命体的测距和呼吸频率估计方法
CN115327504B (zh) * 2022-07-04 2024-05-10 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于测量条件参数的海杂波幅度分布无类型化预测方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4400700A (en) * 1981-06-08 1983-08-23 The United States Of America As Represented By The Secretary Of The Army Doppler frequency analysis of radar signals
US4737788A (en) * 1985-04-04 1988-04-12 Motorola, Inc. Helicopter obstacle detector
US5376940A (en) * 1985-12-19 1994-12-27 Hughes Aircraft Company Helicopter recognition radar processor
US6081221A (en) * 1998-08-27 2000-06-27 Zrnic; Dusan S. Signal processing and systematic phase coding for mitigation of range and velocity ambiguities in Doppler weather radars
NL1020763C2 (nl) * 2002-06-05 2003-12-08 Thales Nederland Bv Werkwijze en -systeem voor het deconvolueren van een met een onregelmatige PRT bemonsterd signaal, en de toepassingsmogelijkheden ervan.

Also Published As

Publication number Publication date
IL156293A (en) 2008-12-29
KR100966060B1 (ko) 2010-06-28
US20040056794A1 (en) 2004-03-25
JP2004037457A (ja) 2004-02-05
ZA200304397B (en) 2004-02-13
KR20030094116A (ko) 2003-12-11
DE60331277D1 (de) 2010-04-01
CA2431094C (en) 2011-07-26
NL1020763C2 (nl) 2003-12-08
EP1369703A2 (en) 2003-12-10
ATE458204T1 (de) 2010-03-15
US6803876B2 (en) 2004-10-12
CA2431094A1 (en) 2003-12-05
IL156293A0 (en) 2004-01-04
EP1369703A3 (en) 2005-01-05
EP1369703B1 (en) 2010-02-17

Similar Documents

Publication Publication Date Title
ES2340672T3 (es) Metodo y sistemas de desconvolucion para periodos irregulares de repeticion de impulsos, tratamiento doppler y de ecos parasitos.
Knill et al. High range and Doppler resolution by application of compressed sensing using low baseband bandwidth OFDM radar
JP6177467B1 (ja) ドップラー能力を高めたレーダー動作
Angelosante et al. Estimating multiple frequency-hopping signal parameters via sparse linear regression
Wang et al. Interleaved OFDM radar signals for simultaneous polarimetric measurements
JP2005530164A (ja) 物体検出用システムにおける妨害を抑制するための方法
CN113640752B (zh) 一种基于脉间相位频谱双捷变的波形设计方法
Li et al. Fast adaptive pulse compression based on matched filter outputs
Mohr et al. Optimized complementary waveform subsets within an FM noise radar CPI
Barbarossa et al. Parameter estimation of 2D multi-component polynomial phase signals: An application to SAR imaging of moving targets
Ojowu et al. RFI suppression for synchronous impulse reconstruction UWB radar using RELAX
Akhtar et al. Formation of range-doppler maps based on sparse reconstruction
Cao et al. Wavelet-based waveform for effective sidelobe suppression in radar signal
Levanon Noncoherent radar pulse compression based on complementary sequences
Dwyer Range and Doppler information from fourth-order spectra
Zhou et al. High-resolution UWB ranging based on phase-only correlator
WO2014184760A1 (en) Coherent radar
RU2740782C1 (ru) Способ радиолокационной съёмки Земли и околоземного пространства радиолокатором с синтезированной апертурой антенны в неоднозначной по дальности полосе с селекцией движущихся целей на фоне отражений от подстилающей поверхности и радиолокатор с синтезированной апертурой антенны для его реализации
Yang et al. Parameter identifiability of monostatic MIMO chaotic radar using compressed sensing
CA2441328A1 (en) Agile prt deconvolution method and systems, and its uses
Atzeni et al. Optical signal-processing through dual-channel ultrasonic light modulators
Skinner et al. Matched FSK/PSK radar
RU2107926C1 (ru) Способ импульсной радиолокации системой фазоманипулированных сигналов
McEliece et al. Some combinatorial aspects of spread-spectrum communication systems
Davies et al. Wideband acoustic communications dispelling narrowband myths