Articulo de referencia

Monopulso de comparación de fase

El monopulso de comparación de fase es una técnica utilizada en aplicaciones de radiofrecuencia (RF), como radares y localización de dirección , para estimar con precisión la di...

El monopulso de comparación de fase es una técnica utilizada en aplicaciones de radiofrecuencia (RF), como radares y localización de dirección , para estimar con precisión la dirección de llegada de una señal a partir de la diferencia de fase de la señal medida en dos (o más) antenas separadas [ 1 ] o, más comúnmente, a partir de los centros de fase desplazados de una antena de arreglo. El monopulso de comparación de fase se diferencia del monopulso de comparación de amplitud en que el primero utiliza centros de fase desplazados con una dirección de apuntamiento del haz común, mientras que el segundo utiliza un centro de fase común y direcciones de apuntamiento del haz desplazadas. [ 2 ]

En la técnica de comparación de fase monopulso, un arreglo se subdivide en sub-arreglos, formando así un canal de "suma" y un canal de "diferencia" o "del". En un arreglo lineal, cada sub-arreglo consta de la mitad de los elementos, divididos por la mitad. En un arreglo planar, estos sub-arreglos corresponden a los cuatro cuadrantes del arreglo, cada uno con 1/4 de sus elementos. En un arreglo lineal, la salida de cada sub-arreglo se suma para formar el canal de "suma", y las mismas salidas se restan para formar el canal de "del". La relación monopulso se obtiene dividiendo la parte imaginaria del canal del entre la parte real del canal de suma. Esta relación proporciona una señal de error que indica con alta precisión el ángulo real del objetivo con respecto al centro del haz. Para una matriz planar, se forma un canal suma como la suma de las salidas de los cuatro cuadrantes, pero se forman dos canales delta, uno para la dimensión de elevación y otro para la dimensión de azimut ortogonal. Se forman dos relaciones monopulso, al igual que con una matriz lineal, cada una indicando el ángulo de desviación en una dimensión desde el centro del haz. [ 3 ]

Existen algunos conceptos erróneos comunes sobre el monopulso de comparación de fase. Primero, solo se forma un haz. El procesamiento monopulso se realiza completamente con la señal recibida en la variedad de la matriz y la red de formación de haces. Hablando en términos de una sola dimensión para mayor claridad, como con una matriz lineal, la señal es recibida por la matriz y se suma en cada una de las dos submatrices con centros de fase desplazados. El canal de suma se forma simplemente sumando las salidas de estas dos submatrices, y el resultado es exactamente el mismo que si toda la matriz se hubiera sumado inicialmente en un solo paso. El canal de eliminación se forma simplemente restando estas mismas salidas de submatrices. Segundo, el monopulso de comparación de fase técnicamente no realiza una comparación de fase, sino que simplemente divide el canal de eliminación por el canal de suma para llegar a una relación donde se codifica la información angular. [ 4 ] La siguiente derivación matemática debería aclarar por qué esto es así.

Matemáticas

Patrón de suma

Podemos definir el patrón de haz ( factor de matriz ) de una matriz lineal uniforme (ULA) con N elementos, como: [ 5 ]

Bθ(θ)=wHvθ(θ)=norte=0norte1wnorte[vθ(θ)]norte=norte=0norte1wnortemij(nortenorte12)2πλddoosθ{\displaystyle B_{\theta }\left(\theta \right)={\vec {w}}^{H}{\vec {v}}_{\theta }\left(\theta \right)=\sum _{n=0}^{N-1}w_{n}^{*}\left[{\vec {v}}_{\theta }\left(\theta \right)\right]_{n}=\sum _{n=0}^{N-1}w_{n}^{*}e^{j\left(n-{\frac {N-1}{2}}\right){\frac {2\pi }{\lambda }}dcos\theta }}, dóndevθ{\displaystyle {\vec {v}}_{\theta }}es el vector de variedad de matriz yw{\displaystyle {\vec {w}}}es un vector de pesos complejos que representan ajustes de amplitud y fase aplicados a cada elemento de la antena. El vector de variedad,vθ{\displaystyle {\vec {v}}_{\theta }}, encapsula completamente todas las propiedades espaciales de la matriz. d{\displaystyle d}es la distancia entre los elementos del arreglo, yθ{\displaystyle \theta }es el ángulo de llegada de una onda plana incidente , definido a partir del frente de onda, es decir,θ=90{\displaystyle \theta =90^{\circ }}es una señal de la matriz en dirección perpendicular.

Es común realizar una sustitución de variables paraψ{\displaystyle \psi }-espacio, dondeψ=2πλddoosθ{\displaystyle \psi ={\frac {2\pi }{\lambda }}dcos\theta }y por lo tanto tenemos:

Bψ(ψ)=norte=0norte1wnortemij(nortenorte12)ψ{\displaystyle B_{\psi }\left(\psi \right)=\sum _{n=0}^{N-1}w_{n}^{*}e^{j\left(n-{\frac {N-1}{2}}\right)\psi }}

y podemos ver más fácilmente queψ{\displaystyle \psi }es simplemente el desfase entre elementos adyacentes.norte12{\displaystyle {\frac {N-1}{2}}}El término simplemente hace referencia a la fase absoluta con respecto al centro físico del conjunto.

Observe que este resultado es el mismo si, en lugar de eso, primero sumamos cada mitad de la matriz y luego sumamos esos resultados.

Bψ(ψ)=norte=0norte21wnortemij(nortenorte12)ψ+norte=norte2norte1wnortemij(nortenorte12)ψ{\displaystyle B_{\psi }\left(\psi \right)=\sum _{n=0}^{{\frac {N}{2}}-1}w_{n}^{*}e^{j\left(n-{\frac {N-1}{2}}\right)\psi }+\sum _{n={\frac {N}{2}}}^{N-1}w_{n}^{*}e^{j\left(n-{\frac {N-1}{2}}\right)\psi }}

El vector de peso es una combinación de un vector de dirección que dirige el haz en una dirección determinada,ψS{\displaystyle \psi _{S}}, utilizando ajustes de fase y una atenuación de amplitud que se aplica a menudo para reducir los lóbulos laterales . Por lo tanto,[w]norte=anortemij(nortenorte12)ψS{\displaystyle \left[{\vec {w}}\right]_{n}=a_{n}e^{j\left(n-{\frac {N-1}{2}}\right)\psi _{S}}}, y

Bψ(ψΔ)=mij(norte12)ψΔnorte=0norte1anortemijnorteψΔ{\displaystyle B_{\psi }\left(\psi _{\Delta }\right)=e^{j\left({\frac {N-1}{2}}\right)\psi _{\Delta }}\sum _{n=0}^{N-1}a_{n}e^{-jn\psi _{\Delta }}}, dóndeψΔ=ψSψ{\displaystyle \psi _{\Delta }=\psi _{S}-\psi }.

Ahora podemos ver claramente que el patrón del haz, enψ{\displaystyle \psi }-espacio, es el equivalente espacial de la transformada discreta de Fourier (DTFT) del vector de atenuación de amplitud de la matriz multiplicado por un término de fase lineal . La ventaja deψ{\displaystyle \psi }-La ventaja es que la forma del haz es idéntica sin importar hacia dónde se dirija, y es solo una función de la desviación de la fase objetivo deseada con respecto a la fase objetivo real.

Supongamos ahora una matriz normalizada sin atenuación conanorte=1norte{\displaystyle a_{n}={\frac {1}{N}}}. Se puede demostrar fácilmente que el patrón del haz es la conocida función sinc (asinc) con alias:

Bψ(ψΔ)=1nortesinorte(norteψΔ2)sinorteψΔ2{\displaystyle B_{\psi }\left(\psi _{\Delta }\right)={\frac {1}{N}}{\frac {sin\left(N{\frac {\psi _{\Delta }}{2}}\right)}{sin{\frac {\psi _{\Delta }}{2}}}}}

Este patrón también se conoce, para fines de monopulso, como el patrón "suma", ya que se obtuvo sumando todos los elementos. En adelante suprimiremos elΔ{\displaystyle \Delta }subíndice y en su lugar utilice soloψ{\displaystyle \psi }con el entendimiento de que representa la desviación entre la fase objetivo dirigida y la fase objetivo real.

Patrón de diferencia

Desarrollemos ahora el patrón de "diferencia" o "del" monopulso dividiendo la matriz en dos mitades iguales llamadas submatrices. Podríamos haber derivado fácilmente el patrón de suma determinando primero el patrón de cada submatriz individualmente y sumando estos dos resultados. En la práctica monopulso, esto es lo que realmente se hace. Se deja al lector la tarea de demostrar quevψ(ψ){\displaystyle {\vec {v}}_{\psi }\left(\psi \right)}es simétrica conjugada, por lo que puede reescribirse en términos de solo su primera mitad,vψ1(ψ){\displaystyle {\vec {v}}_{\psi _{1}}\left(\psi \right)}utilizando una matriz de intercambio,J{\displaystyle {\textbf {J}}}, que "invierte" este vector.

J=[001100100]{\displaystyle {\textbf {J}}={\begin{bmatrix}0&\cdots &0&1\\\vdots &\ddots &1&0\\0&\cdot ^{\cdot ^{\cdot }}&\ddots &\vdots \\1&0&\cdots &0\end{bmatrix}}}

Tenga en cuenta queJJ=I{\displaystyle {\textbf {J}}\cdot {\textbf {J}}={\textbf {I}}}. Suponiendo que N es par (podríamos desarrollarlo igualmente usando un N impar), [ 6 ]

vψ(ψ)=[vψ1(ψ)Jvψ1(ψ)]{\displaystyle {\vec {v}}_{\psi }\left(\psi \right)={\begin{bmatrix}{\vec {v}}_{\psi _{1}}\left(\psi \right)\\\cdots \\{\textbf {J}}{\vec {v}}_{\psi _{1}}^{*}\left(\psi \right)\end{bmatrix}}}

Si asumimos que la matriz de pesos también es simétrica conjugada (una buena suposición), entonces

w=[w1Jw1]{\displaystyle {\vec {w}}={\begin{bmatrix}{\vec {w}}_{1}\\\cdots \\{\textbf {J}}{\vec {w}}_{1}^{*}\end{bmatrix}}}

y el patrón de haz suma se puede reescribir como: [ 7 ]

Bψ(ψ)=Σψ(ψ)=wHvψ(ψ)=[w1Hw1TJ][vψ1(ψ)Jvψ1(ψ)]=w1Hvψ1(ψ)+w1Tvψ1(ψ)=2Rmi[w1Hvψ1(ψ)]{\displaystyle B_{\psi }\left(\psi \right)=\Sigma _{\psi }\left(\psi \right)={\vec {w}}^{H}{\vec {v}}_{\psi }\left(\psi \right)={\begin{bmatrix}{\vec {w}}_{1}^{H}&\vdots &{\vec {w}}_{1}^{T}{\textbf {J}}\end{bmatrix}}{\begin{bmatrix}{\vec {v}}_{\psi _{1}}\left(\psi \right)\\\cdots \\{\textbf {J}}{\vec {v}}_{\psi _{1}}^{*}\left(\psi \right)\end{bmatrix}}={\vec {w}}_{1}^{H}{\vec {v}}_{\psi _{1}}\left(\psi \right)+{\vec {w}}_{1}^{T}{\vec {v}}_{\psi _{1}}^{*}\left(\psi \right)=2Re\left[{\vec {w}}_{1}^{H}{\vec {v}}_{\psi _{1}}\left(\psi \right)\right]}

El patrón de diferencia o "del" se puede inferir fácilmente del patrón de suma simplemente invirtiendo el signo de los pesos para la segunda mitad de la matriz:

Δψ(ψ)=[w1Hw1TJ][vψ1(ψ)Jvψ1(ψ)]=w1Hvψ1(ψ)w1Tvψ1(ψ)=2Imetro[w1Hvψ1(ψ)]{\displaystyle \Delta _{\psi }\left(\psi \right)={\begin{bmatrix}{\vec {w}}_{1}^{H}&\vdots &-{\vec {w}}_{1}^{T}{\textbf {J}}\end{bmatrix}}{\begin{bmatrix}{\vec {v}}_{\psi _{1}}\left(\psi \right)\\\cdots \\{\textbf {J}}{\vec {v}}_{\psi _{1}}^{*}\left(\psi \right)\end{bmatrix}}={\vec {w}}_{1}^{H}{\vec {v}}_{\psi _{1}}\left(\psi \right)-{\vec {w}}_{1}^{T}{\vec {v}}_{\psi _{1}}^{*}\left(\psi \right)=2Im\left[{\vec {w}}_{1}^{H}{\vec {v}}_{\psi _{1}}\left(\psi \right)\right]}

Suponiendo nuevamente queanorte=1norte{\displaystyle a_{n}={\frac {1}{N}}}, se puede demostrar que el patrón del se reduce a:

Δψ(ψ)=2norteImetro[norte=0norte21mij(nortenorte12)ψ]=2nortesinorte2(norteψ4)sinorteψ2{\displaystyle \Delta _{\psi }\left(\psi \right)={\frac {2}{N}}Im\left[\sum _{n=0}^{{\frac {N}{2}}-1}e^{-j\left(n-{\frac {N-1}{2}}\right)\psi }\right]={\frac {2}{N}}{\frac {sin^{2}\left(N{\frac {\psi }{4}}\right)}{sin{\frac {\psi }{2}}}}}
Patrones de suma y diferencia (del) monopulso

Relación monopulso

La relación monopulso se forma de la siguiente manera:

ΔψΣψ=2nortesinorte2(norteψ4)sinorteψ21nortesinorte(norteψ2)sinorteψ2=2sinorte2(norteψ4)sinorte(norteψ2)=1doos(norteψ2)sinorte(norteψ2)=tanorte(norteψ4){\displaystyle {\frac {\Delta _{\psi }}{\Sigma _{\psi }}}={\frac {{\frac {2}{N}}{\frac {sin^{2}\left(N{\frac {\psi }{4}}\right)}{sin{\frac {\psi }{2}}}}}{{\frac {1}{N}}{\frac {sin\left(N{\frac {\psi }{2}}\right)}{sin{\frac {\psi }{2}}}}}}={\frac {2sin^{2}\left(N{\frac {\psi }{4}}\right)}{sin\left(N{\frac {\psi }{2}}\right)}}={\frac {1-cos\left(N{\frac {\psi }{2}}\right)}{sin\left(N{\frac {\psi }{2}}\right)}}=tan\left(N{\frac {\psi }{4}}\right)}

Se observa que, dentro del ancho de haz de 3 dB del sistema, la relación monopulso es casi lineal. De hecho, para muchos sistemas, una aproximación lineal es suficiente. También se puede apreciar que la relación monopulso es continua dentro del ancho de haz entre nulos, pero presenta asíntotas en los nulos del haz. Por lo tanto, la relación monopulso solo es precisa para medir el ángulo de desviación de un objetivo dentro del lóbulo principal del sistema. Sin embargo, los objetivos detectados en las zonas periféricas del sistema, si no se corrigen, producirán resultados erróneos.

Relación monopulso dentro de 1 ancho de haz del eje de respuesta principal

Concepto de operaciones

Antes de realizar el procesamiento monopulso, un sistema debe detectar primero un objetivo, lo cual hace normalmente utilizando el canal suma. Todas las mediciones típicas que realiza un sistema no monopulso se llevan a cabo utilizando el canal suma, por ejemplo, distancia, Doppler y ángulo. Sin embargo, la medición del ángulo es limitada, ya que el objetivo podría estar en cualquier punto dentro del ancho del haz suma, por lo que el sistema solo puede asumir que la dirección de apuntamiento del haz coincide con el ángulo real del objetivo. En realidad, por supuesto, el ángulo real del objetivo y el ángulo de direccionamiento del haz serán diferentes.

Por lo tanto, un procesador monopulso funciona detectando y midiendo primero la señal objetivo en el canal suma. Luego, solo cuando sea necesario para los objetivos detectados, mide la misma señal en el canal "del", dividiendo la parte imaginaria de este resultado por la parte real del canal "suma", y luego convirtiendo esta relación en un ángulo de desviación utilizando las siguientes relaciones:

ψΔ=ψSψ=4norteardotanorte(ΔψΣψ){\displaystyle \psi _{\Delta }=\psi _{S}-\psi ={\frac {4}{N}}arctan\left({\frac {\Delta _{\psi }}{\Sigma _{\psi }}}\right)}

y

θ=ardodoos((ψSψΔ)λ2πd)=ardodoos(λ2πd(2πλddoosθS4norteardotanorte(ΔψΣψ)))=ardodoos(doosθS2λnorteπdardotanorte(ΔψΣψ)){\displaystyle \theta =arccos\left({\frac {\left(\psi _{S}-\psi _{\Delta }\right)\lambda }{2\pi d}}\right)=arccos\left({\frac {\lambda }{2\pi d}}\left({\frac {2\pi }{\lambda }}dcos\theta _{S}-{\frac {4}{N}}arctan\left({\frac {\Delta _{\psi }}{\Sigma _{\psi }}}\right)\right)\right)=arccos\left(cos\theta _{S}-{\frac {2\lambda }{N\pi d}}arctan\left({\frac {\Delta _{\psi }}{\Sigma _{\psi }}}\right)\right)}

Este ángulo de desviación, que puede ser positivo o negativo, se suma al ángulo de apuntamiento del haz para obtener una estimación más precisa del ángulo de rumbo real del objetivo. Por supuesto, si el conjunto es bidimensional, como un conjunto planar, existen dos canales del, uno para la elevación y otro para el acimut, y por lo tanto se forman dos relaciones monopulso.

Véase también

Referencias

  1. Mahafza, Bassem R. (1998). Introducción al análisis de radar; Procesamiento de señales de radar en ingeniería eléctrica . CRC Press . pág.  251. ISBN 0-8493-1879-3.
  2. Sherman, Samuel M. Principios y técnicas de monopulso, 2.ª edición . Artech House. pág. 72. 
  3. Sherman, Samuel M. Principios y técnicas de monopulso, 2.ª edición . Artech House.
  4. Sherman, Samuel M. Principios y técnicas de monopulso, 2.ª edición . Artech House. págs. 70–74 . 
  5. Van Trees, HL (2002). Procesamiento óptimo de matrices, Parte IV de la teoría de detección, estimación y modulación . John Wiley & Sons, Inc. pág. 39. 
  6. Van Trees, HL (2002). Procesamiento óptimo de matrices, Parte IV de la teoría de detección, estimación y modulación . John Wiley & Sons, Inc. pág. 40. 
  7. Van Trees, HL (2002). Procesamiento óptimo de matrices, Parte IV de la teoría de detección, estimación y modulación . John Wiley & Sons, Inc. pág. 40.