01.05.2017 Views

563489578934

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Sec. 6–6 The Gaussian Random Process 449<br />

or<br />

y Hx<br />

(6–115)<br />

where the elements of the N × N matrix H are related to the impulse response of the linear<br />

network by<br />

h ij = [h(t i - l j )] ¢l<br />

(6–116)<br />

We will now show that y is described by an N-dimensional Gaussian PDF when x is described<br />

by an N-dimensional Gaussian PDF. This may be accomplished by using the theory for a multivariate<br />

functional transformation, as given in Appendix B. From Eq. (B–99), the PDF of y is<br />

The Jacobian is<br />

f y (y) = f x(x)<br />

|J(y/x)|<br />

`<br />

x = H -1 y<br />

dy 1 dy 1<br />

Á dy 1<br />

dx 1 dx 2 dx N<br />

Ja y dy 2 dy 2<br />

Á dy 2<br />

x b = Det G dx 1 dx 2 dx N W = Det [H] ! K<br />

o o o<br />

dy N dy N<br />

Á dy N<br />

dx 1 dx 2 dx N<br />

(6–117)<br />

(6–118)<br />

where K is a constant. In this problem, J (y/x) is a constant (not a function of x) because<br />

y = Hx is a linear transformation. Thus,<br />

f y (y) = 1<br />

|K| f x(H -1 y)<br />

or<br />

1<br />

f y (y) =<br />

(6–119)<br />

(2p) N/2 1/2<br />

|K| |Det C x |<br />

e-(1/2) [(H-1 y -m x ) T -1<br />

C x (H -1 y-m x )]<br />

where the subscript x has been appended to the quantities that are associated with x(t). But we<br />

know that m y = Hm x , and from matrix theory, we have the property [AB] T = B T A T , so that<br />

the exponent of Eq. (6–119) becomes<br />

where<br />

or<br />

- 1 2 [(y - m y) T (H -1 ) T ]C x -1 [ H -1 (y - m y )] = - 1 2 [(y - m y) T C ȳ 1 (y - m y )]<br />

C ȳ 1 = (H -1 ) T C x -1 H -1<br />

C y = HC x H T<br />

(6–120)<br />

(6–121)<br />

(6–122)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!