[PDF]
ECE 460: Communication and Information Theory
Prof. B.-P. Paris
Homework 8
Solution
![∫ 1 ∫ 1 ∫ 1∫ 1
fX,Y(x,y)dxdy = K (x + y)dxdy
0 0 [0∫ 10∫ 1 ∫ 1∫ 1 ]
= K xdxdy + ydxdy
[ 0 0 0 0]
1 2||1 1 1 2||1 1
= K --x ||y |0 + -y ||x|0
2 0 2 0
= K](hwsol_460140x.png)
2)

3) By exploiting the symmetry of fX,Y and the fact that it has to integrate to 1, one immediately sees that the answer to this question is 1/2. The “mechanical” solution is:

4)

The region over which we integrate in order to find p(X > Y,X + 2Y > 1) is marked with an A in the following figure.
Thus
![∫ 1 ∫ x
p(X > Y,X + 2Y > 1) = (x + y)dxdy
13 1-2x
∫ 1 [ 1 - x 1 2 1 - x 2 ]
= 1 x(x - -----) + --(x - (-----) ) dx
∫ 3 ( 2 2) 2
= 1 15-x2 - 1x - 1- dx
13 8 4 8
49
= ----
∫108∫
p(X + 2Y > 1) = 1 1 (x + y)dxdy
0 1-2x
∫ 1 [ 1 - x 1 1 - x ]
= x(1 - -----) + --(1 - (------)2) dx
∫0 ( 2 )2 2
1 3- 2 3- 3-
= 0 8 x + 4x + 8 dx
3 1 ||1 3 1 ||1 3 ||1
= --× -x3 || + --× --x2|| + -x||
8 3 0 4 2 0 8 0
7-
= 8](hwsol_460145x.png)
5) When X = Y the volume under integration has measure zero and thus

6) Conditioned on the fact that X = Y , the new p.d.f of X is

In words, we re-normalize fX,Y (x,y) so that it integrates to 1 on the region characterized by X = Y . The result depends only on x. Then p(X > 1 2|X = Y ) = ∫ 1∕21f X|X=Y (x)dx = 3∕4.
7)

8) FX(x|X + 2Y > 1) = p(X ≤ x,X + 2Y > 1)∕p(X + 2Y > 1)
![∫ ∫
p(X ≤ x,X + 2Y > 1) = x 1 (v + y)dvdy
0 1-2v
∫ x[3 3 3]
= --v2 + -v + -- dv
0 8 4 8
1- 3 3- 2 3-
= 8x + 8 x + 8x](hwsol_460149x.png)

![∫ 1
E [X |X + 2Y > 1] = xfX (x|X + 2Y > 1)dx
∫01( )
= 3x3 + 6-x2 + 3x
0 7 7 7
3 1 4||1 6 1 3||1 3 1 2||1 17
= --× -x || + --× -x || + --× --x || = ---
7 4 0 7 3 0 7 2 0 28](hwsol_460151x.png)
1) fX,Y (x,y) is a PDF, hence its integral over the supporting region of x, and y is 1.

2)

3)

4) If x < y then fX|Y (x|y) = 0. If x ≥ y, then with u = x - y ≥ 0 we obtain

5)
![∫ ∞ ∫ ∞
E [X |Y = y] = xe -x+ydx = ey xe- xdx
y[ y ]
||∞ ∫ ∞
= ey - xe-x|| + e- xdx
y y
= ey(ye -y + e-y) = y + 1](hwsol_460156x.png)
6) In this part of the problem we will use extensively the following definite integral

![∫ ∞ ∫ ∞ -x-y ∫ ∞ - y ∫ ∞ - x
E [XY ] = 0 y xy2e dxdy = 0 2ye y xe dxdy
∫ ∞ ∫ ∞ ∫ ∞
= 2ye -y(ye-y + e-y)dy = 2 y2e-2ydy + 2 ye-2ydy
0 0 0
= 2 1-2! + 2 1-1! = 1
23 22](hwsol_460158x.png)
![∫ ∞ - x - x ∫ ∞ -x ∫ ∞ -2x
E [X ] = 2 xe (1 - e )dx = 2 xe dx - 2 xe dx
0 0 0
= 2 - 2 1--= 3-
∫ 22 2
∞ -2y 1-- 1-
E [Y ] = 2 0 ye dy = 222 = 2
2 ∫ ∞ 2 -x - x ∫ ∞ 2 -x ∫ ∞ 2 -2x
E [X ] = 2 x e (1 - e )dx = 2 x e dx - 2 x e dx
0 0 0
= 2 ⋅ 2! - 2-12! = 7
∫ 23 2
2 ∞ 2 -2y 1-- 1-
E [Y ] = 2 0 y e dy = 2 232! = 2](hwsol_460159x.png)
![3 1 1
COV (X, Y ) = E [XY ] - E [X ]E [Y ] = 1 - --⋅--= --
2 2 4](hwsol_460160x.png)
and
![COV (X, Y ) 1
ρX,Y = -----2----------2-1∕2-----2----------2-1∕2 = √---
(E [X ] - (E [X ]) ) (E [Y ] - (E[Y ]) ) 5](hwsol_460161x.png)
____________________________________________________________________________________________________
| REAL*8 | x,t,a,q,pi,p,b1,b2,b3,b4,b5 |
| PARAMETER | (p=.2316419d+00, b1=.31981530d+00, |
| + | b2=-.356563782d+00, b3=1.781477937d+00, |
| + | b4=-1.821255978d+00, b5=1.330274429d+00) |
| C- |
| pi=4.*atan(1.) |
| C-INPUT |
| PRINT*, | ’Enter -x-’ |
| READ*, | x |
| C- |
| t=1./(1.+p*x) |
| a=b1*t + b2*t**2. + b3*t**3. + b4*t**4. + b5*t**5. |
| q=(exp(-x**2./2.)/sqrt(2.*pi))*a |
| C-OUTPUT |
| PRINT*, | q |
| C- |
| STOP |
| END |
The results of this approximation along with the actual values of Q(x) (taken from text Table 4.1) are tabulated in the following table. As it is observed a very good approximation is achieved.
| x | Q(x) | Approximation |
| 1. | 1.59 × 10-1 | 1.587 × 10-1 |
| 1.5 | 6.68 × 10-2 | 6.685 × 10-2 |
| 2. | 2.28 × 10-2 | 2.276 × 10-2 |
| 2.5 | 6.21 × 10-3 | 6.214 × 10-3 |
| 3. | 1.35 × 10-3 | 1.351 × 10-3 |
| 3.5 | 2.33 × 10-4 | 2.328 × 10-4 |
| 4. | 3.17 × 10-5 | 3.171 × 10-5 |
| 4.5 | 3.40 × 10-6 | 3.404 × 10-6 |
| 5. | 2.87 × 10-7 | 2.874 × 10-7 |
![∫ ∞ ∫ ∞
-∞ -∞ fX,Y(x,y)dxdy
∫ 0 ∫ 0 x2+y2 ∫ ∞ ∫ ∞ x2+y2
= K-e---2--dxdy + K-e- -2--dxdy
-∞∫ -∞ π ∫ 0 0∫ π ∫
K-- 0 - x2 0 - y2 K-- ∞ - x2- ∞ - y2
= π -∞ e 2 dx -∞ e 2 dx + π 0 e 2 dx 0 e 2 dx
K [ 1√ --- ]
= --- 2(-- 2π)2 = K
π 2](hwsol_460162x.png)
2) If x < 0 then


e-x2
2 which implies that fX(x) is a zero-mean
Gaussian random variable with variance 1. Since fX,Y (x,y) is symmetric to
its arguments and the same is true for the region of integration we conclude
that fY (y) is a zero-mean Gaussian random variable of variance
1.
3) fX,Y (x,y) has not the same form as a binormal distribution. For xy < 0, fX,Y (x,y) = 0 but a binormal distribution is strictly positive for every x, y.
4) The random variables X and Y are not independent for if xy < 0 then
fX(x)fY (y)
0 whereas fX,Y (x,y) = 0.
5)
![1 ∫ 0 ∫ 0 x2+y2 1 ∫ ∞ ∫ ∞ x2+y2
E[XY ] = -- XY e---2--dxdy + -- e---2--dxdy
π ∫- ∞ -∞ ∫ π 0∫ 0 ∫
1- 0 - x22 0 - y22- 1- ∞ - x22 ∞ - y22
= π - ∞ Xe dx -∞ Y e dy + π 0 Xe dx 0 Y e dy
1 1 2
= --(- 1 )(- 1) +--= --
π π π](hwsol_460167x.png)
0 and
E[X] = E[Y ] = 0, so that E[XY ] - E[X]E[Y ]
0.
6) In general fX|Y (x,y) = fX,Y (x,y) fY (y) . If y > 0, then

If y ≤ 0, then

Thus

which is not a Gaussian distribution.