
El procesamiento adaptativo espacio-temporal ( STAP ) es una técnica de procesamiento de señales muy utilizada en sistemas de radar . Consiste en algoritmos de procesamiento adaptativo de la matriz para facilitar la detección de objetivos. El procesamiento de señales de radar se beneficia del STAP en áreas donde la interferencia es un problema (por ejemplo, ruido de fondo , interferencias , etc.). Mediante una aplicación cuidadosa del STAP, es posible lograr mejoras de sensibilidad de varios órdenes de magnitud en la detección de objetivos.
STAP emplea una técnica de filtrado bidimensional mediante una antena de matriz en fase con múltiples canales espaciales. El acoplamiento de estos múltiples canales espaciales con formas de onda Doppler pulsadas da origen al término "espacio-tiempo". Aplicando las estadísticas del entorno de interferencia, se genera un vector de ponderación STAP adaptativo. Este vector se aplica a las muestras coherentes recibidas por el radar.
Historia
La teoría STAP fue publicada por primera vez por Lawrence E. Brennan e Irving S. Reed a principios de la década de 1970. En el momento de la publicación, tanto Brennan como Reed trabajaban en Technology Service Corporation (TSC). Si bien se introdujo formalmente en 1973, [ 1 ] sus raíces teóricas se remontan a 1959. [ 2 ]
Motivación y aplicaciones
Para el radar terrestre, los retornos desordenados tienden a estar en CC, lo que los hace fácilmente discriminados por la Indicación de Objetivo Móvil (MTI) . [ 3 ] Por lo tanto, se puede utilizar un filtro de muesca en el bin de Doppler cero. [ 2 ] Las plataformas aéreas con movimiento propio experimentan un movimiento relativo del desorden terrestre que depende del ángulo, lo que resulta en un acoplamiento ángulo-Doppler en la entrada. [ 2 ] En este caso, el filtrado 1D no es suficiente, ya que el desorden puede superponerse al Doppler del objetivo deseado desde múltiples direcciones. [ 2 ] La interferencia resultante se denomina típicamente "cresta de desorden", ya que forma una línea en el dominio ángulo-Doppler. [ 2 ] Las señales de interferencia de banda estrecha también son una fuente de interferencia y exhiben una correlación espacial significativa. [ 4 ] Por lo tanto, se debe considerar el ruido del receptor y la interferencia, y los procesadores de detección deben intentar maximizar la relación señal-a-interferencia y ruido (SINR) .
Aunque se desarrollaron principalmente para radares, las técnicas STAP tienen aplicaciones en sistemas de comunicaciones. [ 5 ]
Teoría básica

STAP es esencialmente un filtrado en el dominio espacio-temporal. [ 2 ] Esto significa que estamos filtrando sobre múltiples dimensiones, y se deben emplear técnicas de procesamiento de señales multidimensionales. [ 6 ] El objetivo es encontrar los pesos espacio-temporales óptimos enespacio -dimensional, dondees el número de elementos de antena (nuestros grados de libertad espaciales) yes el número de coeficientes del intervalo de repetición de pulsos (PRI) (nuestros grados de libertad temporales), para maximizar la relación señal-a-interferencia-ruido (SINR) . [ 2 ] Por lo tanto, el objetivo es suprimir el ruido, el desorden, los inhibidores, etc., manteniendo al mismo tiempo el retorno de radar deseado. Puede pensarse como un filtro de respuesta de impulso finito (FIR) 2D , con un filtro FIR 1D estándar para cada canal (canales espaciales dirigidos desde una matriz dirigida electrónicamente o elementos individuales), y los coeficientes de estos filtros FIR 1D corresponden a múltiples retornos (espaciados en el tiempo PRI). [ 1 ] Tener grados de libertad tanto en el dominio espacial como en el dominio temporal es crucial, ya que el desorden puede estar correlacionado en el tiempo y el espacio, mientras que los inhibidores tienden a estar correlacionados espacialmente (a lo largo de una dirección específica). [ 1 ]
En la primera figura se muestra un ejemplo sencillo y trivial de STAP, paraEste es un ejemplo idealizado de un patrón de dirección, donde la respuesta del conjunto se ha dirigido a la respuesta objetivo ideal.[ 2 ] Desafortunadamente , en la práctica, esto es una simplificación excesiva, ya que la interferencia que se debe superar al dirigir los nulos mostrados no es determinista, sino de naturaleza estadística. [ 2 ] Esto es lo que requiere que STAP sea una técnica adaptativa. Nótese que incluso en este ejemplo idealizado, en general, debemos dirigirnos sobre el plano ángulo-Doppler 2D en puntos discretos para detectar objetivos potenciales (moviendo la ubicación del lóbulo principal sinc 2D que se muestra en la figura), y hacerlo para cada uno de los intervalos de alcance en nuestro sistema.
El diagrama funcional básico se muestra a la derecha. Para cada antena, normalmente se completa un paso de conversión descendente y conversión analógica-digital. Luego, se utiliza un filtro FIR 1-D con elementos de retardo de longitud PRI para cada canal de antena orientado. Los pesos ordenados lexicográficamenteason los grados de libertad que se deben resolver en el problema STAP. Es decir, STAP tiene como objetivo encontrar los pesos óptimos para el conjunto de antenas. Se puede demostrar que, para un dadomatriz de covarianza de interferencia,Los pesos óptimos que maximizan el SINR se calculan como:
dóndees un escalar que no afecta al SINR. [ 2 ] La entrada óptima del detector viene dada por:
dóndees una instantánea espacio-temporal de los datos de entrada. La principal dificultad de STAP radica en resolver e invertir la matriz de covarianza de interferencia, que suele ser desconocida.[ 1 ] Surgen otras dificultades cuando la matriz de covarianza de interferencia está mal condicionada, lo que hace que la inversión sea numéricamente inestable. [ 5 ] En general, este filtrado adaptativo debe realizarse para cada uno de los intervalos de rango inequívocos del sistema, para cada objetivo de interés (coordenadas ángulo-Doppler), lo que supone una enorme carga computacional. [ 4 ] Pueden producirse pérdidas de dirección cuando los retornos reales del objetivo no caen exactamente en uno de los puntos de nuestro plano ángulo-Doppler 2D que hemos muestreado con nuestro vector de dirección.. [ 1 ]
Aproches
Los distintos enfoques se pueden desglosar por taxonomía de procesamiento, [ 7 ] o por simplificación del espacio de datos / fuentes de datos. [ 2 ]
Métodos directos
La solución óptima consiste en utilizar todos los grados de libertad procesando el filtro adaptativo en los elementos de la antena. Para los métodos directos adaptativos, la inversión de matriz de muestra (SMI) utiliza la matriz de covarianza de interferencia estimada (de muestra) en lugar de la matriz de covarianza de interferencia real. [ 8 ] Esto se debe a que la matriz de covarianza de interferencia real no se conoce en la práctica. [ 1 ] Si se conoce por algún medio, no es necesario estimarla y los pesos óptimos se fijan. Esto a veces se denomina variación independiente de los datos. La variación dependiente de los datos estima la matriz de covarianza de interferencia a partir de los datos. En los sistemas de comunicaciones MIMO, esto se puede hacer mediante una secuencia de entrenamiento. [ 5 ] El detector clarividente se obtiene cuando la matriz de covarianza se conoce perfectamente y se define como:
dóndees la estadística de instantánea espacio-temporal para elcelda de rango bajo la hipótesis de interferencia únicamente,. [ 1 ] Para SMI, la matriz de covarianza de interferencia para laLa celda de rango que consiste en las estadísticas de ruido interferente, desorden y bloqueadores se estima de la siguiente manera: [ 4 ]
dóndeson los datos de entrenamiento obtenidos del procesador de entrada para elcelda de rango. Por lo tanto, se promedian las instantáneas espacio-temporales que rodean la celda de rango deseada. Nótese que la instantánea espacio-temporal de la celda de rango deseada se excluye normalmente (así como varias celdas adicionales, o "celdas de guarda") para evitar el blanqueamiento de las estadísticas. [ 1 ]
El principal problema de los métodos directos es la gran complejidad computacional asociada con la estimación e inversión de matrices formadas a partir de muchos grados de libertad (gran número de elementos y/o pulsos). [ 1 ] Además, para métodos dondedebe estimarse utilizando muestras de datos, el número de muestras necesarias para lograr un error particular depende en gran medida de la dimensionalidad de la matriz de covarianza de interferencia. [ 4 ] Como resultado, para sistemas de alta dimensión, esto puede requerir un número inalcanzable de celdas de rango inequívocas. [ 1 ] Además, estas celdas de datos adyacentes deben contener estadísticas estacionarias en función del rango, lo cual rara vez es una buena suposición para el gran número de celdas requeridas (para una degradación de SINR de 3 dB desde STAP clarividente óptimo). [ 2 ] [ 1 ]
Métodos de rango reducido
Los métodos de rango reducido buscan superar las cargas computacionales del método directo al reducir la dimensionalidad de los datos o el rango de la matriz de covarianza de interferencia. [ 2 ] Esto se puede lograr formando haces y realizando STAP en el espacio de haces. [ 7 ] Se pueden usar métodos pre y post Doppler en el espacio de haces. Los métodos post Doppler también se pueden usar en la entrada completa del elemento de antena para reducir los datos solo en esta dimensión. Un ejemplo popular es la antena de centro de fase desplazado (DPCA), que es una forma de STAP independiente de datos en el espacio de haces, pre-Doppler. [ 7 ] El objetivo es realizar la formación de haces de tal manera que el haz parezca estacionario mientras el radar aerotransportado está en movimiento durante períodos de tiempo discretos para que el ruido aparezca sin Doppler. [ 2 ] Sin embargo, los errores de fase pueden causar una degradación significativa ya que el algoritmo no es adaptativo a los datos devueltos. [ 2 ] Se pueden utilizar muchos otros métodos para reducir el rango de la matriz de covarianza de interferencia, por lo que todos los métodos en la categoría de rango reducido pueden considerarse como una simplificación de la matriz de covarianza que se va a invertir:
Los métodos post-Doppler descomponen el problema STAP a partir de unproblema de filtrado adaptativo afiltros adaptativos individuales de longitud(unproblema del filtro adaptativo). [ 2 ] Al realizar un procesamiento Doppler fijo, los filtros adaptativos se vuelven solo espaciales. [ 2 ] Dado que la respuesta objetivo ya está dirigida a una ubicación angular-Doppler específica, la dimensionalidad se puede reducir mediante el preprocesamiento de múltiples bins Doppler y ángulos que rodean este punto. [ 4 ] Además de reducir la dimensionalidad del procesador adaptativo, esto a su vez reduce la cantidad de fotogramas de datos de entrenamiento necesarios al estimar la matriz de covarianza de interferencia, ya que esta cantidad depende de la dimensión. [ 2 ]
Dado que estos métodos reducen la dimensionalidad de los datos, son inherentemente subóptimos. [ 1 ] Existen varias técnicas para comparar el rendimiento de los métodos de rango reducido y los métodos directos estimados con el STAP clarividente (directo con conocimiento perfecto de la matriz de covarianza de interferencia y el vector de dirección objetivo), principalmente basadas en la pérdida SINR. [ 1 ] Un ejemplo de ello es
donde hemos tomado la relación del SINR evaluado con los pesos subóptimosy el SINR evaluado con los pesos óptimos[ 1 ] Nótese que , en general, esta cantidad es estadística y se debe tomar la esperanza para hallar la pérdida SINR promedio. La pérdida SINR clarividente también se puede calcular tomando la relación entre la SINR óptima y la SNR del sistema, lo que indica la pérdida debida a la interferencia. [ 1 ]
Métodos basados en modelos
También existen métodos basados en modelos que intentan forzar o explotar la estructura de la matriz de interferencia de covarianza. El método más aplicable de forma general es la estructura de matriz de atenuación de covarianza. [ 2 ] El objetivo es modelar la interferencia de forma compacta, momento en el que puede procesarse utilizando técnicas de componentes principales o SMI de carga diagonal (donde se añade una matriz diagonal aleatoria de pequeña magnitud para intentar estabilizar la matriz antes de invertirla). [ 2 ] Este modelado tiene el beneficio adicional de decorrelacionar la fuga del subespacio de interferencia (ISL) y es resistente al movimiento de desorden interno (ICM). [ 2 ] El método de componentes principales primero aplica el análisis de componentes principales para estimar los autovalores y autovectores dominantes, y luego aplica una atenuación de covarianza y añade un nivel de ruido estimado:
dóndees elvalor propio estimado mediante PCA,es el asociadovector propio estimado mediante PCA,implica la multiplicación elemento por elemento de matricesy,es la matriz de covarianza estimada yes el nivel de ruido estimado. [ 2 ] La estimación de la atenuación de covarianzaPuede resultar complicado, dependiendo de la complejidad del modelo subyacente que intenta emular el entorno de interferencia. Se recomienda al lector consultar [ 2 ] para obtener más información sobre este tema en particular. Una vez que esta atenuación esté suficientemente modelada, también puede aplicarse a la adaptación SMI más simple de CMT de la siguiente manera:
dóndees la matriz estimada SMI típica que se ve en el método directo aproximado,es el factor de carga diagonal, yes la matriz identidad del tamaño apropiado. Cabe señalar que esto tiene como objetivo mejorar el método SMI estándar, donde SMI utiliza un número menor de intervalos de rango en su promedio que la técnica SMI estándar. Dado que se utilizan menos muestras en los datos de entrenamiento, la matriz a menudo requiere estabilización en forma de carga diagonal. [ 2 ]
Ejemplos más restrictivos implican modelar la interferencia para forzar estructuras de Toeplitz, y pueden reducir considerablemente la complejidad computacional asociada al procesamiento al explotar esta estructura. [ 2 ] Sin embargo, estos métodos pueden sufrir debido a la falta de coincidencia del modelo, o el ahorro computacional puede verse anulado por el problema del ajuste del modelo (como el problema no lineal del ajuste a una matriz de Toeplitz o de bloques de Toeplitz) y la estimación del orden. [ 2 ]
Aplicaciones modernas
A pesar de tener casi 40 años de existencia, STAP cuenta con aplicaciones modernas.
comunicaciones MIMO
Para canales dispersivos, las comunicaciones de entrada múltiple y salida múltiple pueden formular soluciones STAP. La compensación de canal selectiva en frecuencia se puede utilizar para extender las técnicas de ecualización tradicionales para sistemas SISO que utilizan STAP. [ 5 ] Para estimar la señal transmitidaEn un receptor MIMO, podemos ponderar linealmente nuestra entrada espacio-temporal.con matriz de ponderacióncomo sigue
para minimizar el error cuadrático medio (ECM) . [ 5 ] Usando STAP con una secuencia de entrenamiento, la matriz de ponderación óptima estimada (coeficientes STAP) viene dada por: [ 5 ]
radar MIMO
STAP se ha extendido para radar MIMO para mejorar la resolución espacial para el ruido, utilizando técnicas STAP de radar SIMO modificadas. [ 9 ] Se requieren nuevos algoritmos y formulaciones que se apartan de la técnica estándar debido al gran rango del subespacio de interferencia-ruido creado por los arreglos virtuales de radar MIMO, [ 9 ] que normalmente implica explotar la estructura diagonal por bloques de la matriz de covarianza de interferencia MIMO para dividir el gran problema de inversión de matrices en problemas más pequeños. En comparación con los sistemas de radar SIMO, que tendrántransmitir grados de libertad yreciben grados de libertad, para un total deLos sistemas de radar MIMO tienengrados de libertad, lo que permite una resolución espacial adaptativa mucho mayor para la mitigación de la confusión. [ 9 ]
Véase también
Referencias
- 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Melvin, WL, Descripción general de A STAP , Revista IEEE Aerospace and Electronic Systems – Número especial de tutoriales, vol. 19 , n.° 1, enero de 2004, págs.19–35.
- 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 Guerci, JR, Procesamiento adaptativo espacio-temporal para radar , Artech House Publishers, 2003. ISBN 1-58053-377-9.
- ↑ Richards, MA, Scheer, JA y Holm, WA, Principios del radar moderno , SciTech Publishing, 2010. ISBN 1-89112-152-9.
- 1 2 3 4 5 Richards, MA, Fundamentos del procesamiento de señales de radar , McGraw-Hill Education, 2014. ISBN 0-07179-832-3.
- 1 2 3 4 5 6 Bliss, DW y Govindasamy, S., Comunicaciones inalámbricas adaptativas: canales y redes MIMO , Cambridge University Press, 2013. ISBN 1-10703-320-9.
- ↑ Dudgeon, DE y Mersereau, RM, Procesamiento de señales digitales multidimensionales , Serie de procesamiento de señales de Prentice-Hall, 1984. ISBN 0-13604-959-1.
- 1 2 3 Ward, J., Procesamiento adaptativo espacio-temporal para radar aerotransportado, Coloquio del IEE sobre procesamiento adaptativo espacio-temporal (Ref. No. 1998/241), abril de 1998, págs. 2/1–2/6.
- ↑ Van Trees, HL, Optimum Array Processing, Wiley, NY, 2002.
- 1 2 3 Li, J. y Stoica, P., Procesamiento de señales de radar MIMO , John Wiley & Sons, 2009. ISBN 0-47017-898-1.
Lecturas adicionales
- Brennan, LE y IS Reed, Teoría del radar adaptativo , IEEE AES-9, págs. 237–252, 1973
- Guerci, JR, Procesamiento adaptativo espacio-temporal para radar , Artech House Publishers, 2003. ISBN 1-58053-377-9.
- Klemm, Richard, Principios del procesamiento adaptativo espacio-temporal , IEE Publishing, 2002. ISBN 0-85296-172-3.
- Klemm, Richard, Aplicaciones del procesamiento adaptativo espacio-temporal , IEE Publishing, 2004. ISBN 0-85296-924-4.
- Melvin, WL, Una visión general de STAP, Revista IEEE Aerospace and Electronic Systems – Número especial de tutoriales, Vol. 19, No. 1, enero de 2004, págs. 19–35.
- Michael Parker, Fundamentos del radar – Parte 4: Procesamiento adaptativo espacio-temporal , EETimes, 28/06/2011
- Procesamiento de señales de radar