Active sonar likelihood function
We will find that the natural logarithm of the measurement likelihood ratio simplifies the detection expression considerably:
The determinant of
is
It is convenient to partition
into sub-matrices compatible with r.
The determinant of
is
Which becomes
Using
Sylvester’s determinant theorem :
We obtain:
Which equals
Now we partition the observed ping history y into
so that:
The
Woodbury matrix identity states:
So that
And hence
So that the log likelihood function becomes
Putting the noise variance inside the observation vector yields for the log likelihood function
The log-likelihood ratio shows that the assumed echo shape,
, the assumed energy of the signal
, the ambient noise density
, and the expected location in time of the echo needed to select
are the parameters required to evaluate the log-likelihood function. The assumed echo shape and energy can vary by hypothesis h, but the noise properties
have been assumed to be the same for each decision hypothesis h.
Signal to noise ratio of the detector
The magnitude squared term,
, is a matched filter, where the observations are cross-correlated with the signal template. The observations are normalized by the noise variance before cross correlation, which is a form of pre-whitening. The term
is the Energy to Noise Density Ratio (ENR) of the detection problem. Note that the signal to noise ratio of the matched filter output can be written as:
This result is a general result for matched filters. A matched filter’s SNR is the Energy to Noise Density Ratio of the problem. The energy of a signal being detected is related to the average amplitude and the duration of the signal. The SNR output of the matched filter is independent of the details of the waveform being detected, only the signal energy and the noise spectral density determine the matched filter response.
The likelihood function is dependent only the ENR as well.
A-priori assumptions
Before the ping history is received, we assess the probability of each hypothesis,
and
. The a-priori information may have come from previous pings, or are probabilities assigned by the sonar system to begin a target search.
Using the logarithm of probability density function ratios simplifies the expressions:
Using the likelihood ratio notation,
Once we compute
, we can declare a target is present, with confidence
by computing:
Target Present if:
Because the target hypothesis space contains many hypotheses, this detection problem is can be considered a composite hypothesis test.
An alternative approach to detection of a target with an unknown range is solved by finding the target range hypothesis with the greatest measurement likelihood as the detection statistic [Kay]. This is referred to as Generalized Maximum Likelihood Ratio Testing (GLRT).
Target Present if:
The GLRT approach is often easier to implement than the Bayes detection approach, because one avoids the integration/summation over a-priori probabilities.