HW3 Answers

From SeedWiki

Jump to: navigation, search

Contents

Part 1 (50 points)

A

1.
\begin{array}{l}
E(X)=0, E(Y)=0\\
V(X)=V(Y)=V(Z+XX)=V(Z)+V(XX)=\frac{20^2}{12}+1=34.\dot 3\\
Cov(X,Y)=E((X-E(X)(Y-E(Y))=E(XY)=E((Z+XX)(Z+YY))=\\
E(Z^2+Z\cdot XX+Z\cdot YY+XX\cdot YY) = E(Z^2)=V(Z)+(E(Z))^2=\frac{20^2}{12}=33.\dot 3\\
\end{array}
Not independent (Cov(X,Y) \ne 0)

2.
\begin{array}{l}
E(X)=0, E(Y)=0\\
V(X)=V(Y)=\frac{20^2}{12}=33.\dot 3\\
Cov(X,Y)=E((X-E(X)(Y-E(Y))=E(XY)=0\\
\end{array}
Not independent (P(Y=3) \ne p(Y=3|X=1))

3.
\begin{array}{l}
E(X)=0, E(Y)=0\\
V(X)=V(Y)=E(X^2)-E(X)^2=16\\
Cov(X,Y)=0\\
\end{array}
Independent

4.
\begin{array}{l}
E(X)=0, E(Y)=0\\
V(X)=V(Y)=E(X^2)-E(X)^2=E(25(\cos^2T))=25E(\cos^2T)
=25\cdot 4\int_0^{\frac{\pi}{2}}\frac{1}{2\pi}\cos^2x dx\\
=\frac{50}{\pi}\int_0^{\frac{\pi}{2}}\cos^2x dx
=\left[\frac{50}{\pi}\left( \frac{1}{2}\cos x\sin x +\frac{1}{2}x \right) \right]_0^{\frac{\pi}{2}}
=\frac{50}{\pi}(\frac{\pi}{4})=12.5\\
Cov(X,Y)=0\\
\end{array}
Not independent

5.
\begin{array}{l}
E(X)=0, E(Y)=0\\
V(X)=25, V(Y)=16\\
Cov(X,Y)=0\\
\end{array}
Independent

B

S1=2
S2=5
S3=1
S4=4
S5=3


Part 2 (20 points)

1

\begin{array}{l}
E(X^2+Y^2)=E(X^2)+E(Y^2)=16+16=32\\
E(\max(X,Y))=\frac{3}{4}(4)+\frac{1}{4}(-4)=3-1=2\\
\end{array}

2

\begin{array}{l}
E(X^2+Y^2)=E(X^2)+E(Y^2)=2E(X^2)=2\int_0^{10}\frac{1}{10}x^2dx \\
=2\left[ \frac{x^3}{3}\cdot \frac{1}{10}\right]_0^{10}=\frac{200}{3}=66.\dot 6\\
E(\max(X,Y))=\frac{3}{4}E(|X|)+\frac{1}{4}E(-|X|)
=\frac{3}{4}E(|X|)-\frac{1}{4}E(|X|)
=\frac{1}{2}E(|X|)=\frac{1}{2}\cdot 5=2.5\\
\end{array}

3

\begin{array}{l}
E(X^2+Y^2)=2E(X^2)=2(V(X)-(E(X))^2)=2V(X)=2\cdot 12.5=25\\
E(\max(X,Y))=2\int_{\frac{1}{4}\pi}^{\frac{3}{4}\pi}5\cdot\frac{1}{2\pi}\sin x dx
= \frac{5}{\pi}\left[\cos x\right]_{\frac{1}{4}\pi}^{\frac{3}{4}\pi}
=\frac{5}{\pi}\sqrt{2}=2.2508\\
\end{array}

Part 3 (50 points)

Using the MATLAB normcdf function, we find that

P(N1>1) = 1 - normcdf(1,0,1) = 0.1587
P(N2>1) = 1 - normcdf(1,0,2) = 0.3085

Running mean(G>1) gives us 0.1610 . Thus, we conclude that G was created using distribution N1.

Sanov bound

p = 0.1587 q = 0.1610 m = size(G,2)

sanov = exp(exp(-m*(q*log(q/p)+(1-q)*log((1-q)/(1-p))))

1.0214e-025

Hoeffding bound

hoeff = exp(-2*(p-q)^2*m)

1.2392e-019

Chernoff bound

chernoff = exp(-m/3*(p-q)^2/p)

6.1213e-011

Personal tools
projects