Phase-comparison monopulse

Phase-comparison monopulse, a simple form of phase-interferometry, is a technique used in radio frequency (RF) applications such as radar and direction finding to accurately estimate the direction of arrival of a signal from the phase difference of the signal measured on two (or more) separated antennas or more typically from displaced phase centers of an array antenna. Phase comparison monopulse differs from amplitude comparison monopulse in that the former uses displaced phase centers with a common beam pointing direction, while the latter uses a common phase center and displaced beam pointing directions.

In phase-comparison monopulse, typically an array is subdivided into sub-arrays, and then a "sum" and a "difference" or "del" channel are formed. For a linear array, these subarrays would each be half of the elements, divided in the middle. For a planar array, these sub-arrays would be the four quadrants of the array, each with 1/4 of the array's elements. In a linear array, the output of each sub-array is summed to form the "sum" channel, and the same outputs are subtracted to form the "del" channel. The monopulse ratio is formed by dividing the imaginary part of the del channel by the real part of the sum channel. This ratio gives an error signal that indicates to a high degree of accuracy the actual target angle as compared to the center of the beam. For a planar array, one sum channel is formed as the sum of the outputs of all four quadrants, but two del channels are formed, one for the elevation dimension and one for the orthogonal azimuth dimension. Two monopulse ratios are formed just as with a linear array, each one indicating the deviation angle in one dimension from the center of the beam.

There are some common misconceptions about phase comparison monopulse. First, only one beam is formed. Monopulse processing is done entirely with the received signal in the array manifold (beam forming network). Speaking in terms of only one dimension for clarity, such as with a linear array, the signal is received by the array and summed into each of two subarrays with displaced phase centers. The sum channel is formed simply by adding these two subarray outputs, and the result is exactly the same as if the entire array was initially summed in one step. The del channel is formed simply by subtracting these same subarray outputs. Second, phase-comparison monopulse doesn't technically actually do a phase comparison, but rather simply divides the del channel by the sum channel to arrive at a ratio wherein the angle information is encoded. The following mathematical derivation should make it clear why this is so.

Mathematics

We can define the beam pattern (array factor) of a uniform linear array (ULA) with N elements, as:

$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 }$, where ${\vec {v}}_{\theta }$ is the array manifold vector and ${\vec {w}}$ is a vector of complex weights representing amplitude and phase adjustments applied to each antenna element. The manifold vector, ${\vec {v}}_{\theta }$, fully encapsulates all of the spatial properties of the array. $d$ is the distance between elements of the array, and $\theta$ is the angle of arrival of an incident plane wave, defined from end-fire, i.e., $\theta =90^{\circ }$ is a signal from array broadside.

File:Uniform linear array (ULA) geometry.png

Uniform linear array (ULA) geometry

It is common to perform a variable substitution to $\psi$-space, where $\psi ={\frac {2\pi }{\lambda }}dcos\theta$, and therefore we have:

$B_{\psi }\left(\psi \right)=\sum _{n=0}^{N-1}w_{n}^{*}e^{j\left(n-{\frac {N-1}{2}}\right)\psi }$

and we can more easily see that $\psi$ is simply the phase shift between adjacent elements. The ${\frac {N-1}{2}}$ term simply references the absolute phase to the physical center of the array.

The weight vector is a combination of a steering vector that steers the beam in a target direction, $\psi _{T}$, using phase adjustments and an amplitude taper that is often applied to reduce sidelobes. Thus, $\left[{\vec {w}}\right]_{n}=a_{n}e^{j\left(n-{\frac {N-1}{2}}\right)\psi _{T}}$, and

$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 }}$, where $\psi _{\Delta }=\psi _{T}-\psi$.

We can clearly see now that the beam pattern, in $\psi$-space, is the spatial equivalent of the discrete time Fourier transform (DTFT) of the array amplitude tapering vector times a linear phase term. The advantage of $\psi$-space is that the beam shape is identical no matter where it is steered, and is only a function of the deviation of the desired target phase from the actual target phase.

Let us now assume an un-tapered, normalized array with $a_{n}={\frac {1}{N}}$. The beam pattern can be easily shown to be the familiar aliased sinc (asinc) function:

$B_{\psi }\left(\psi _{\Delta }\right)={\frac {1}{N}}{\frac {sin\left(N{\frac {\psi _{\Delta }}{2}}\right)}{sin{\frac {\psi _{\Delta }}{2}}}}$

This pattern is also known, for monopulse purposes, as the "sum" pattern, as it was obtained by summing all of the elements together. Going forward we will suppress the $\Delta$ subscript and instead use only $\psi$ with the understanding that it represents the deviation of the steered target phase and the actual target phase.

Difference Pattern

Let us now develop the monopulse "difference" or "del" pattern by dividing the array into two equal halves called subarrays. We could have just as easily derived the sum pattern by first determining the pattern of each subarray individually and adding these two results together. In monopulse practice, this is what is actually done. The reader is left to show that ${\vec {v}}_{\psi }\left(\psi \right)$ is conjugate symmetric, so it can be re-written in terms of only its first half, ${\vec {v}}_{\psi _{1}}\left(\psi \right)$ using an exchange matrix, ${\textbf {J}}$, that "flips" this vector. Note that ${\textbf {J}}\cdot {\textbf {J}}={\textbf {I}}$. Assuming that N is even (we could just as easily develop this using an odd N),

${\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}}$

If we assume that the weight matrix is also conjugate symmetric (a good assumption), then

${\vec {w}}={\begin{bmatrix}{\vec {w}}_{1}\\\cdots \\{\textbf {J}}{\vec {w}}_{1}^{*}\end{bmatrix}}$

and the sum beam pattern can be rewritten as:

$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]$

The difference or "del" pattern can easily be inferred from the sum pattern simply by flipping the sign of the weights for the second half of the array:

$\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]$

Again assuming that $a_{n}={\frac {1}{N}}$, the del pattern can be shown to reduce to:

$\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}}}}$

Mono-pulse sum and difference (del) patterns

Monopulse Ratio

The monopulse ratio is formed as:

${\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)$

Mono-pulse ratio within 1 beamwidth of the main response axis

Therefore, a monopulse processor functions by measuring the target signal on the "del" channel, dividing the imaginary part of this result by the real part of the same target signal measured on the "sum" channel, then converting this ratio to a deviation angle using the relationships:

$\psi ={\frac {4}{N}}atan\left({\frac {\Delta _{\psi }}{\Sigma _{\psi }}}\right)$

and

$\theta =acos\left({\frac {\psi \lambda }{2\pi d}}\right)$

Note that this deviation angle is added to the beam pointing angle to arrive at the more accurate estimate of the actual target bearing angle.

See also

• Very-long-baseline interferometry
• Amplitude-comparison monopulse
• Monopulse radar