Let’s talk about ambiguity and correlation. The ambiguity function is a core component of radar signal processing practice and theory. The autocorrelation function and the cyclic autocorrelation function, are key elements of generic signal processing and cyclostationary signal processing, respectively. Ambiguity and correlation both apply a quadratic functional to the data or signal of interest, and they both weight that quadratic functional by a complex exponential (sine wave) prior to integration or summation.
Are they the same thing? Well, my answer is both yes and no.
The Ambiguity Function
The common definition of the ambiguity function is
where is the delay or lag parameter and
is a frequency parameter that we will interpret below. Notice that the integral is over all time. This implies that
should be an energy signal. An energy signal has finite energy, so that the integral
exists and is greater than zero. Because it exists, it is also finite. Any bounded integrable (or summable in discrete time) signal that has finite temporal support is an energy signal; think of a rectangle, triangle, rapidly decaying exponential pulse, etc. An energy signal doesn’t have to have finite support, but that is a common property of energy signals, especially those used in radar.
The cross ambiguity function replaces one of the terms in (1) with another energy signal
:
Consider a radar-like signal consisting of a pulse with finite temporal support, such as a burst of a sine wave or a burst of a chirp waveform. The signal is transmitted, strikes a reflective target, experiences the Doppler effect, and a portion of its energy propagates back to a receiver (often co-located with the transmitter). So the received reflected pulse is
where is a complex-valued constant,
is the round-trip delay from transmitter to target to receiver, and
is the Doppler shift frequency, which is a good approximation to the Doppler effect in many cases of practical interest.
If you substitute from (4) and
into (3), and ignore the additive noise
, you can see that the ambiguity magnitude will peak when
and
. For this reason, the
parameter is called the Doppler frequency or Doppler parameter.
The “auto” ambiguity in (1) is useful for studying the radar-processing properties of candidate radar waveforms, whereas the cross version in (3) is more directly related to radar signal processing itself.
When using the ambiguity functions, the range of the Doppler frequency is limited to those Doppler shifts that are physically possible given the possible relative motions between the target and the transmitter/receiver. This is a (thankfully, for computational cost reasons) small range. For example, for a stationary transmitter/receiver, a target moving at mph relative to the transmitter/receiver, and an operating frequency of
GHz, the Doppler shift is about
kHz.
I think the ambiguity function is most often called just that: ambiguity or ‘the ambiguity function.’ But the cross version is often called the CAF, and some people call the auto ambiguity function the ‘complex ambiguity function,’ or CAF. So that’s two CAFs.
The Cyclic Autocorrelation Function
Back to our wheelhouse. The cyclic autocorrelation function is defined by an infinite-time average:
Notice here that if is an energy signal, the limit will be zero. On the other hand, if
is a power signal, the limit can be finite and non-zero, depending on the nature of the signal(s) in
, the value of the delay
, and crucially on the value of the cycle frequency
. The power of a signal is the temporal average of the energy:
The cyclic cross correlation function is defined by
In the correlation functions, the frequency parameter–the cycle frequency–that gives rise to non-zero results is a function of the modulation types and parameters of the signals that are present in and
. Cycle frequencies can be quite small, such as those for narrowband digital signals operating at
kHz, or quite large, such as those for wideband digital signals such as CDMA (
MHz), WCDMA (
MHz), WiFi DSSS BPSK (
MHz), etc. The range of
is typically much larger than the range of
.
The cyclic autocorrelation function is also called the CAF. So that’s three CAFs, increasing the potential for confusion. Which CAF are we talking about? If they are all the same, there is no problem. If they aren’t, then we need to know how they differ.
Estimators
To answer the post-title question, we need to introduce estimators of the ambiguity and cyclic autocorrelation functions. The two functions involve an integral over all time (a sum for discrete-time versions), so that natural estimators are simply the infinite-time definitions with the infinite intervals replaced by finte intervals. For the cyclic autocorrelation function, this is the estimator that follows simply from the definition in the fraction-of-time probability framework (The Literature [R67, R8, R131]), and is also the familiar estimator that arises in the context of a cycloergodic cyclostationary stochastic process.
The cyclic correlation estimators:
For the ambiguity, we have the estimators:
The Question of Equivalence for Infinite-Time Mathematical Objects
When we think of ambiguity and cyclic correlation in terms of infinite-time mathematical operations ((1), (3), (5), and (7)), they are not equivalent since they produce very different answers for identical inputs. Ambiguity in this context applies to energy signals (think: signals of finite duration) and cyclic correlation applies to power signals (think: persistent random processes or signals). Ambiguity will blow up if applied to a power signal, and cyclic correlation will be zero for an energy signal.
Conclusion: They are not the same thing.
The Question of Equivalence for Finite-Time Estimators
When we think of ambiguity and cyclic correlation in terms of finite-time estimators, there is an equivalence in that one function is a scaled version of the other. Let’s look at the ambiguity for . First, impose symmetry with respect to
:
or
We recognize the integral in (13) as the estimate of the cyclic autocorrelation multiplied by ,
So the estimates differ by a scaling factor as well as a phase shift of
. Since the user of the estimator selects
and
, this phase factor is easily compensated. Therefore, the estimates of ambiguity and cyclic correlation contain the same information.
Conclusion: They are the same thing.
Related Topics
In the theory of complex-valued cyclostationary signals, the involved random variables and random signals are not proper. This means that their conjugate correlation functions can be non-zero. The normal correlation between random variables and
is
The conjugate correlation is defined as the correlation between and the conjugate of
, or
:
As we’ve explained elsewhere on the CSP Blog, this impropriety leads to the conjugate cyclic autocorrelation and conjugate cross correlation functions
and
It is unclear whether there is any useful analog of these conjugate correlation functions in the world of ambiguity processing. So in this sense, the conclusion is that the cyclic autocorrelation functions are not the same as the ambiguity function.
Let me know what you think in the Comments section below.
So in fourier transformed space do the peaks of lagged correlations of frequency show up the doppler shift, and of phase show up the delay?
Sounds similar to camera autofocus and automatic sub-pixel image registration.
The convolution version would be the correlation after reversing one time-series; which leaves frequency unaffected but flips its phase.
I confess I do not grasp the meaning of this question. Can you rephrase and use mathematical notation? In your comments, you can use LaTeX to format equations and symbols. Start a mathematical object (equation, symbol, etc.) with a dollar sign, type the word ‘latex’, include the latex commands, and end with a dollar sign.
“* latex \alpha *” produces
, where * should be replaced by $
On re-reading I’ve realised I had missed your explanation of the doppler parameter. Is the additive frequency shift an approximation?
Up to constants and scaling the fourier transform A of a is defined
and the Convolution operator in time domain is equivalent to pointwise multiplication in the frequency domain
If signal
then
,
a scaling of frequencies.
I would write
and the convolution theorem for Fourier transforms says
But, yes, I agree that convolution implies that each frequency component in the signal
is scaled by the corresponding frequency component in the signal
. In the context of filtering, identify
as the “signal” and
as the “impulse response of the filter” (linear time-invariant transformation).
I think by “additive frequency shift” you mean a simple frequency shift like those implied by the multiplication of a signal by a complex sine wave,
. This kind of simple frequency shift is a good approximation to the full Doppler effect (which affects different frequency components of the signal differently in general) under mild assumptions on the bandwidth, center frequency, and involved speeds.
Hi Chad,
great article, thank you! Helped me to clarify my understanding of both, the ambiguity and the cyclic autocorrelation function.
Jonas
Thank you Jonas!