Chapter 8

Linear filters -- applications

[10/0%]
LaTeX
Problem 8:1

Why is the assumption that R0(x)R_{0}(x) is zero outside the interval ∣x∣≤d<f0\left|x\right| \leq d < f_{0} important in (8.21)?

Context restated from the book's Section 8.3.3 (Theorem 8.4): let {X(t)}\left\{ X(t)\right\} be a narrow-banded stationary Gaussian process with mean 0 and a spectral density symmetric around ±f0\pm f_{0},

RX(f)=12R0(f+f0)+12R0(f−f0), R_{X}(f) = \frac{1}{2} R_{0}(f+f_{0}) + \frac{1}{2} R_{0}(f-f_{0}),

where R0(x)R_{0}(x) is an even function that is zero outside the interval ∣x∣≤d\left|x\right| \leq d, with d<f0d < f_{0}. Let {Y(t)}\left\{ Y(t)\right\} be the Hilbert transform of {X(t)}\left\{ X(t)\right\}, and let {X1(t)}\left\{ X_{1}(t)\right\}, {X2(t)}\left\{ X_{2}(t)\right\} be the independent Gaussian processes, each with spectral density R0(f)R_{0}(f), obtained from the complex representation Z(t)=X(t)+iY(t)=ei2πf0t(X1(t)+iX2(t))Z(t) = X(t)+iY(t) = e^{i2\pi f_{0} t}\left(X_{1}(t)+iX_{2}(t)\right). Equation (8.21) is the resulting cross-covariance function,

C(X1(t),X2(t+τ))=2∫0∞sin⁡(2π(f−f0)τ)R0(f−f0) df=0,for all τ, C\left(X_{1}(t), X_{2}(t+\tau )\right) = 2\int _{0}^{\infty } \sin \left(2\pi (f-f_{0})\tau \right) R_{0}(f-f_{0}) \, df = 0, \quad \text{for all } \tau ,

which is used to show that {X1(t)}\left\{ X_{1}(t)\right\} and {X2(t)}\left\{ X_{2}(t)\right\} are uncorrelated (hence independent, being jointly Gaussian).

?
Problem 8:2

A stationary sequence {S(t),t=0,±1,±2…}\left\{ S(t), t=0,\pm 1,\pm 2 \ldots \right\} has expectation zero and covariance function,

rS(τ)={2,τ=0,−1,∣τ∣=1,0,otherwise. r_{S}(\tau ) = \begin{cases} 2, & \tau = 0, \\ -1, & \left|\tau \right| = 1, \\ 0, & \text{otherwise.} \end{cases}

The process is disturbed by colored noise {N(t),t=0,±1,±2…}\left\{ N(t), t=0,\pm 1,\pm 2 \ldots \right\} with covariance function,

rN(τ)={2,τ=0,1,∣τ∣=1,0,otherwise. r_{N}(\tau ) = \begin{cases} 2, & \tau = 0, \\ 1, & \left|\tau \right| = 1, \\ 0, & \text{otherwise.} \end{cases}

The process and the disturbance are independent. Determine the frequency function H(f)H(f), −1/2<f≤1/2-1/2 < f \leq 1/2, and the impulse response h(t)h(t), t=0,±1,±2…t = 0,\pm 1,\pm 2 \ldots, for the optimal Wiener filter that minimizes E[(Y(t)−S(t))2]\mathbb {E}\left[\left(Y(t)-S(t)\right)^{2}\right] when Y(t)Y(t) is the output signal from the filter with X(t)=S(t)+N(t)X(t) = S(t)+N(t) as input signal.

?
Problem 8:3

Suppose that we want to measure the temperature, X(t)X(t), in a manufacturing process. Since the measurements are corrupted by additive noise, N(t)N(t), independent of the temperature, we reconstruct the true temperature with an optimal filter. Assuming rX(τ)=e−2∣τ∣r_{X}(\tau ) = e^{-2\left|\tau \right|} and rN(τ)=e−20∣τ∣r_{N}(\tau ) = e^{-20\left|\tau \right|}, find the covariance function of the output signal, {Y(t)}\left\{ Y(t)\right\}.

?
Problem 8:4

Suppose we want to transmit a bandlimited stochastic signal, S(t)S(t), on a channel disturbed by flicker noise. The spectral density of the signal is

RS(f)={1,100≤∣f∣≤1000,0,otherwise, R_{S}(f) = \begin{cases} 1, & 100 \leq \left|f\right| \leq 1000, \\ 0, & \text{otherwise,} \end{cases}

and that of the flicker noise is RN(f)=100/∣f∣R_{N}(f) = 100/\left|f\right| up to high frequencies.

?
(a)

Design an optimal noise reducing filter. Make a plot of the frequency response and compute the SNR (signal-to-noise ratio).

(b)

The following transfer function is not optimal,

G(f)={1,100≤∣f∣≤1000,0,otherwise. G(f) = \begin{cases} 1, & 100 \leq \left|f\right| \leq 1000, \\ 0, & \text{otherwise.} \end{cases}

Compute the SNR of the output signal with this filter.

Problem 8:5

An AR(2)-process {Xt}\left\{ X_{t}\right\} has a covariance function, estimated to r^X(0)=4\widehat{r}_{X}(0) = 4, r^X(±1)=2\widehat{r}_{X}(\pm 1) = 2, and r^X(±2)=−0.5\widehat{r}_{X}(\pm 2) = -0.5. Determine from these values, the parameters for the AR-process, and find the frequency function (in real-valued form) for a Wiener filter which in the best way reconstructs XtX_{t} from the measured signal Yt=Xt+NtY_{t} = X_{t}+N_{t}, when NtN_{t} is white noise with expected value 0 and variance 1. Assume that {Xt}\left\{ X_{t}\right\} and {Nt}\left\{ N_{t}\right\} are independent.

?
Problem 8:6

A deterministic signal, s(t)s(t), is corrupted by additive white noise with spectral density N0N_{0}. The signal is defined by

s(t)={Ae−bt,t≥0,0,otherwise. s(t) = \begin{cases} Ae^{-bt}, & t \geq 0, \\ 0, & \text{otherwise.} \end{cases}

The signal and the noise is processed by a linear optimal filter. We read the output of this filter at time t=T>0t = T > 0. The filter has the property that the SNR is maximized at t=Tt = T.

?
(a)

Determine the filter and the maximal SNR.

(b)

The filter used in (a) is non-causal. Make a small adjustment in order to change this, and compute the new SNR.

(c)

How large must TT be, in order for the SNR to be larger than 99 of the maximum SNR?

Problem 8:7

A sender transmits binary characters to a receiver in a digital communication system. The signal, s(t)s(t), see Figure 8.11, is used by the sender to represent a binary one and the absence of this signal is used to represent zero. We model the received signal, X(t)X(t), by

X(t)=sent signal+N(t),0≤t≤T, X(t) = \text{sent signal} + N(t), \quad 0 \leq t \leq T,

where {N(t),t∈R}\left\{ N(t), t \in \mathbb {R}\right\} is Gaussian noise with spectral density R0R_{0}. Determine the impulse response of a matched causal filter and compute the SNR (signal-to-noise ratio) of the output signal.

The signal, s(t) in Exercise 8:7.The signal, s(t) in Exercise 8:7.

?
Problem 8:8

A deterministic signal, sks_{k}, in discrete time, is given by

sk={1,k=0,1,2,3,0,otherwise. s_{k} = \begin{cases} 1, & k = 0,1,2,3, \\ 0, & \text{otherwise.} \end{cases}

We would like to determine whether this signal has been sent or not, when one receives the signal sk+Nks_{k}+N_{k} if the signal was sent and NkN_{k} if it was not. We assume that NkN_{k} is standardized white Gaussian noise.

?
(a)

We use a filter, defined by its impulse response:

h(k)={1,k=2,3,0,otherwise h(k) = \begin{cases} 1, & k = 2,3, \\ 0, & \text{otherwise} \end{cases}

and choose the decision time k=3k = 3. If the output of the filter is larger than 1 at k=3k = 3 we consider the signal to be present. What is the probability that we erroneously believe that the signal has been sent?

(b)

The filter used in (a) is not optimal. Design an optimal causal filter and specify the time of decision and the threshold level of the discriminator such that the probabilities of the two possible errors are equal. What is the probability of an error?

Problem 8:9

A weakly stationary process {St,t=0,±1,±2…}\left\{ S_{t}, t=0,\pm 1,\pm 2 \ldots \right\} has mean zero and covariance function,

rS(τ)={2,τ=0,−1,∣τ∣=1,0,otherwise. r_{S}(\tau ) = \begin{cases} 2, & \tau = 0, \\ -1, & \left|\tau \right| = 1, \\ 0, & \text{otherwise.} \end{cases}

The process is disturbed by colored noise {Nt,t=0,±1,±2…}\left\{ N_{t}, t=0,\pm 1,\pm 2 \ldots \right\} with expected value zero and covariance function

rN(τ)={2,τ=0,1,∣τ∣=1,0,otherwise. r_{N}(\tau ) = \begin{cases} 2, & \tau = 0, \\ 1, & \left|\tau \right| = 1, \\ 0, & \text{otherwise.} \end{cases}

The process and the disturbance are independent. In communication systems it is quite usual with echo, i.e., you have a signal according to

Xt=St+0.5St−1+Nt, X_{t} = S_{t}+0.5 S_{t-1}+N_{t},

where 0.5St−10.5 S_{t-1} is the damped and delayed original process realization. Determine the coefficients a,ba,b for the non-causal impulse response

h(u)={b,u=−1,a,u=0,0,otherwise. h(u) = \begin{cases} b, & u = -1, \\ a, & u = 0, \\ 0, & \text{otherwise.} \end{cases}

for an optimal filter that minimizes E[(Yt−St)2]\mathbb {E}\left[\left(Y_{t}-S_{t}\right)^{2}\right], when YtY_{t} is the output signal and XtX_{t} is the input signal to the optimal filter.

?
Problem 8:10

Find explicit expressions for the Kalman filter for an AR(1)-process observed with measurement errors,

Xt=ϕ1Xt−1+et, X_{t} = \phi _{1} X_{t-1}+e_{t}, Yt=Xt+ft, Y_{t} = X_{t}+f_{t},

when {et}\left\{ e_{t}\right\} and {ft}\left\{ f_{t}\right\} are independent normal variables with mean zero and Var⁡[et]=σe2\operatorname {Var}\left[e_{t}\right] = \sigma_{e}^{2}, Var⁡[ft]=σf2\operatorname {Var}\left[f_{t}\right] = \sigma_{f}^{2}.

?