Problem Set 7 for Random Processes | ECE 534, Assignments of Electrical and Electronics Engineering

Material Type: Assignment; Class: Random Processes; Subject: Electrical and Computer Engr; University: University of Illinois - Urbana-Champaign; Term: Fall 2004;

Typology: Assignments

Pre 2010

Uploaded on 02/24/2010

koofers-user-s6g
koofers-user-s6g 🇺🇸

10 documents

1 / 2

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
ECE 534 RANDOM PROCESSES SPRING 2004
PROBLEM SET 7 Due Wednesday, December 8
Wiener Filtering of Random Processes
Assigned Reading: Chapter 7 of the notes.
1. A quadratic predictor
Suppose Xis a mean zero, stationary discrete-time random process and that nis an integer with
n1. Consider estimating Xn+1 by a nonlinear one-step predictor of the form
b
Xn+1 =h0+
n
X
k=1
h1(k)Xk+
n
X
j=1
j
X
k=1
h2(j, k)XjXk
(a) Find equations in term of the moments (second and higher, if needed) of Xfor the triple
(h0, h1, h2) to minimize the one step prediction error: E[(Xn+1 b
Xn+1)2].
(b) Explain how your answer to part (a) simplifies if Xis a Gaussian random process.
2. Estimation of a filtered narrowband random process in noise
Suppose Xis a mean zero real-valued stationary Gaussian random process with the spectral density
shown.
1
f
8 Hz
S (2 f)π
X
8 Hz
10 Hz
4
10 Hz
4
(a) Explain how Xcan be simulated on a computer using a pseudo-random number generator that
generates standard normal random variables. Try to use the minimum number per unit time. How
many normal random variables does your construction require per simulated unit time?
(b) Suppose Xis passed through a linear time-invariant system with approximate transfer function
H(2πf ) = 107/(107+f2). Find an approximate numerical value for the power of the output.
(c) Let Zt=Xt+Wtwhere Wis a Gaussian white noise random process, independent of X, with
RW(τ) = δ(τ). Find hto minimize the mean square error E[(Xtb
Xt)2], where b
X=hZ.
(d) Find the mean square error for the estimator of part (c).
3. Interpolating a Gauss Markov process
Let Xbe a real-valued, mean zero stationary Gaussian process with RX(τ) = e−|τ|. Let a > 0.
Suppose X0is estimated by b
X0=c1Xa+c2Xawhere the constants c1and c2are chosen to
minimize the mean square error (MSE).
(a) Use the orthogonality principle to find c1,c2, and the resulting minimum MSE, E[(X0b
X0)2].
(Your answers should depend only on a.)
(b) Use the orthogonality principle again to show that b
X0as defined above is the minimum MSE
estimator of X0given (Xs:|s| a). (This implies that Xhas a two-sided Markov property.)
4. Filtering a WSS signal plus noise
Suppose Xand Nare jointly WSS, mean zero, continuous time random processes with RXN 0.
The processes are the inputs to a system with the block diagram shown, for some transfer functions
K1(ω) and K2(ω):
1
pf2

Partial preview of the text

Download Problem Set 7 for Random Processes | ECE 534 and more Assignments Electrical and Electronics Engineering in PDF only on Docsity!

ECE 534 RANDOM PROCESSES SPRING 2004

PROBLEM SET 7 Due Wednesday, December 8

Wiener Filtering of Random Processes

Assigned Reading: Chapter 7 of the notes.

  1. A quadratic predictor Suppose X is a mean zero, stationary discrete-time random process and that n is an integer with n ≥ 1. Consider estimating Xn+1 by a nonlinear one-step predictor of the form

X̂ n+1 = h 0 +

∑^ n

k=

h 1 (k)Xk +

∑^ n

j=

∑^ j

k=

h 2 (j, k)Xj Xk

(a) Find equations in term of the moments (second and higher, if needed) of X for the triple (h 0 , h 1 , h 2 ) to minimize the one step prediction error: E[(Xn+1 − X̂ n+1)^2 ]. (b) Explain how your answer to part (a) simplifies if X is a Gaussian random process.

  1. Estimation of a filtered narrowband random process in noise Suppose X is a mean zero real-valued stationary Gaussian random process with the spectral density shown.

1 f

8 Hz

S (^) X (2 πf)

8 Hz

10 4 Hz 10 4 Hz

(a) Explain how X can be simulated on a computer using a pseudo-random number generator that generates standard normal random variables. Try to use the minimum number per unit time. How many normal random variables does your construction require per simulated unit time? (b) Suppose X is passed through a linear time-invariant system with approximate transfer function H(2πf ) = 10^7 /(10^7 + f 2 ). Find an approximate numerical value for the power of the output. (c) Let Zt = Xt + Wt where W is a Gaussian white noise random process, independent of X, with RW (τ ) = δ(τ ). Find h to minimize the mean square error E[(Xt − X̂ t)^2 ], where X̂ = h ∗ Z. (d) Find the mean square error for the estimator of part (c).

  1. Interpolating a Gauss Markov process Let X be a real-valued, mean zero stationary Gaussian process with RX (τ ) = e−|τ^ |. Let a > 0. Suppose X 0 is estimated by X̂ 0 = c 1 X−a + c 2 Xa where the constants c 1 and c 2 are chosen to minimize the mean square error (MSE). (a) Use the orthogonality principle to find c 1 , c 2 , and the resulting minimum MSE, E[(X 0 − X̂ 0 )^2 ]. (Your answers should depend only on a.) (b) Use the orthogonality principle again to show that X̂ 0 as defined above is the minimum MSE estimator of X 0 given (Xs : |s| ≥ a). (This implies that X has a two-sided Markov property.)
  2. Filtering a WSS signal plus noise Suppose X and N are jointly WSS, mean zero, continuous time random processes with RXN ≡ 0. The processes are the inputs to a system with the block diagram shown, for some transfer functions K 1 (ω) and K 2 (ω):

N

K 1 + K 2

X Y=X^ out+Nout

Suppose that for every value of ω, Ki(ω) 6 = 0 for i = 1 and i = 2. Because the two subsystems are linear, we can view the output process Y as the sum of two processes, Xout, due to the input X, plus Nout, due to the input N. Your answers to the first four parts should be expressed in terms of K 1 , K 2 , and the power spectral densities SX and SN. (a) What is the power spectral density SY? (b) Find the signal-to-noise ratio at the output (the power of Xout divided by the power of Nout). (c) Suppose Y is passed into a linear system with transfer function H, designed so that the output at time t is X̂ t, the best linear estimator of Xt given (Ys : s ∈ IR). Find H. (d) Find the resulting minimum mean square error. (e) The correct answer to part (d) (the minimum MSE) does not depend on the filter K 2. Why?

  1. Properties of a particular Gaussian process Let X be a zero-mean, wide-sense stationary Gaussian random process in continuous time with autocorrelation function RX (τ ) = (1 + |τ |)e−|τ^ |^ and power spectral density SX (ω) = (2/(1 + ω^2 ))^2. Answer the following questions, being sure to provide justification. (a) Is X mean ergodic in the m.s. sense? (b) Is X a Markov process? (c) Is X differentiable in the m.s. sense? (d) Find the causal, minimum phase filter h (or its transform H) such that if white noise with autocorrelation function δ(τ ) is filtered using h then the output autocorrelation function is RX. (e) Express X as the solution of a stochastic differential equation driven by white noise.
  2. Spectral decomposition and factorization (a) Let x be the signal with Fourier transform given by ̂x(2πf ) =

[ sinc(100f )ej^2 πf T^

]

. Find the energy of x for all real values of the constant T. (b) Find the spectral factorization of the power spectral density S(ω) = (^) ω (^4) +16^1 ω (^2) +100. (Hint: 1 + 3j is a pole of S.)

  1. Causal estimation of a channel input process Let X = (Xt : t ∈ IR) and N = (Nt : t ∈ IR) denote WSS random processes with RX (τ ) = 32 e−|τ^ | and RN (τ ) = δ(τ ). Think of X as an input signal and N as noise, and suppose X and N are orthogonal to each other. Let k denote the impulse response function given by k(τ ) = 2e−^3 τ^ I{τ ≥ 0 }, and suppose an output process Y is generated according to the block diagram shown:

k + X

N

Y

That is, Y = X ∗ k + N. Suppose Xt is to be estimated by passing Y through a causal filter with impulse response function h, and transfer function H. Find the choice of H and h in order to minimize the mean square error.