Matched Filter - Derivation of The Matched Filter

Derivation of The Matched Filter

The following section derives the matched filter for a discrete-time system. The derivation for a continuous-time system is similar, with summations replaced with integrals.

The matched filter is the linear filter, that maximizes the output signal-to-noise ratio.

Though we most often express filters as the impulse response of convolution systems, as above (see LTI system theory), it is easiest to think of the matched filter in the context of the inner product, which we will see shortly.

We can derive the linear filter that maximizes output signal-to-noise ratio by invoking a geometric argument. The intuition behind the matched filter relies on correlating the received signal (a vector) with a filter (another vector) that is parallel with the signal, maximizing the inner product. This enhances the signal. When we consider the additive stochastic noise, we have the additional challenge of minimizing the output due to noise by choosing a filter that is orthogonal to the noise.

Let us formally define the problem. We seek a filter, such that we maximize the output signal-to-noise ratio, where the output is the inner product of the filter and the observed signal .

Our observed signal consists of the desirable signal and additive noise :

Let us define the covariance matrix of the noise, reminding ourselves that this matrix has Hermitian symmetry, a property that will become useful in the derivation:

where denotes the Hermitian (conjugate) transpose of, and denotes expectation. Let us call our output, the inner product of our filter and the observed signal such that

We now define the signal-to-noise ratio, which is our objective function, to be the ratio of the power of the output due to the desired signal to the power of the output due to the noise:

We rewrite the above:

We wish to maximize this quantity by choosing . Expanding the denominator of our objective function, we have

Now, our becomes

We will rewrite this expression with some matrix manipulation. The reason for this seemingly counterproductive measure will become evident shortly. Exploiting the Hermitian symmetry of the covariance matrix, we can write

\mathrm{SNR} = \frac{ | {(R_v^{1/2}h)}^\mathrm{H} (R_v^{-1/2}s) |^2 } { {(R_v^{1/2}h)}^\mathrm{H} (R_v^{1/2}h) },

We would like to find an upper bound on this expression. To do so, we first recognize a form of the Cauchy-Schwarz inequality:

which is to say that the square of the inner product of two vectors can only be as large as the product of the individual inner products of the vectors. This concept returns to the intuition behind the matched filter: this upper bound is achieved when the two vectors and are parallel. We resume our derivation by expressing the upper bound on our in light of the geometric inequality above:

\mathrm{SNR} = \frac{ | {(R_v^{1/2}h)}^\mathrm{H} (R_v^{-1/2}s) |^2 } { {(R_v^{1/2}h)}^\mathrm{H} (R_v^{1/2}h) } \leq \frac{ \left[ 			{(R_v^{1/2}h)}^\mathrm{H} (R_v^{1/2}h) 		\right] 		\left[ 			{(R_v^{-1/2}s)}^\mathrm{H} (R_v^{-1/2}s) 		\right] } { {(R_v^{1/2}h)}^\mathrm{H} (R_v^{1/2}h) }.

Our valiant matrix manipulation has now paid off. We see that the expression for our upper bound can be greatly simplified:

\mathrm{SNR} = \frac{ | {(R_v^{1/2}h)}^\mathrm{H} (R_v^{-1/2}s) |^2 } { {(R_v^{1/2}h)}^\mathrm{H} (R_v^{1/2}h) } \leq s^\mathrm{H} R_v^{-1} s.

We can achieve this upper bound if we choose,

where is an arbitrary real number. To verify this, we plug into our expression for the output :

\mathrm{SNR} = \frac{ | {(R_v^{1/2}h)}^\mathrm{H} (R_v^{-1/2}s) |^2 } { {(R_v^{1/2}h)}^\mathrm{H} (R_v^{1/2}h) } = \frac{ \alpha^2 | {(R_v^{-1/2}s)}^\mathrm{H} (R_v^{-1/2}s) |^2 } { \alpha^2 {(R_v^{-1/2}s)}^\mathrm{H} (R_v^{-1/2}s) } = \frac{ | s^\mathrm{H} R_v^{-1} s |^2 } { s^\mathrm{H} R_v^{-1} s } = s^\mathrm{H} R_v^{-1} s.

Thus, our optimal matched filter is

We often choose to normalize the expected value of the power of the filter output due to the noise to unity. That is, we constrain

This constraint implies a value of, for which we can solve:

yielding

giving us our normalized filter,

If we care to write the impulse response of the filter for the convolution system, it is simply the complex conjugate time reversal of .

Though we have derived the matched filter in discrete time, we can extend the concept to continuous-time systems if we replace with the continuous-time autocorrelation function of the noise, assuming a continuous signal, continuous noise, and a continuous filter .

Read more about this topic:  Matched Filter

Famous quotes containing the word matched:

    The sun set; but set not his hope:
    Stars rose; his faith was earlier up:
    Fixed on the enormous galaxy,
    Deeper and older seemed his eye:
    And matched his sufferance sublime
    The taciturnity of time.
    Ralph Waldo Emerson (1803–1882)