Aufgabe 4.4Z: Höhenlinien der 2D-WDF
Aus LNTwww
Version vom 20. März 2017, 11:47 Uhr von Guenter (Diskussion | Beiträge)
Gegeben ist eine zweidimensionale Gaußsche Zufallsgröße $(x, y)$ mit dem Mittelwert $(0, 0)$ und der 2D–WDF
- $$f_{xy}(x, y) = C\cdot{\rm e}^{-(x^{\rm 2} + y^{\rm 2} +\sqrt{\rm 2}\hspace{0.05cm}\cdot \hspace{0.05cm} x \hspace{0.05cm}\cdot \hspace{0.05cm} y)}.$$
Bekannt ist weiter, dass die beiden Streuungen $\sigma_x$ und $\sigma_y$ jeweils gleich $1$ sind. In der Skizze eingetragen sind:
- eine Höhenlinie dieser WDF für $f_{xy}(x,y) =0.2$,
- die Ellipsenhauptachse (EA), und
- die Korrelationsgerade $y=K(x)$.
Hinweise:
- Die Aufgabe gehört zum Kapitel Zweidimensionale Zufallsgrößen.
- Sollte die Eingabe des Zahlenwertes „0” erforderlich sein, so geben Sie bitte „0.” ein.
- Die hier behandelte Thematik ist in zwei Lernvideos zusammengefasst:
- Gaußsche Zufallsgrößen ohne statistische Bindungen
- Gaußsche Zufallsgrößen mit statistischen Bindungen
Fragebogen
Musterlösung
- 1. Auch ohne die Angabe „σx = σy = 1” könnte man erkennen, dass die beiden Streuungen gleich sind, da im Exponenten der 2D–WDF fxy(x, y) die Koeffizienten bei x2 und y2 gleich sind. Durch Koeffizientenvergleich erhält man mit σx = σy = 1:
- $$\frac{- 2 \rho_{xy}}{\sigma_x\cdot\sigma_y} = \sqrt{2}\hspace{0.3cm}\Rightarrow\hspace{0.3cm} \rho_{xy}=\frac{-1}{\sqrt{2}} \hspace{0.15cm}\underline{\approx -0.707}.$$
- 2. Mit den unter Punkt (1) berechneten Zahlenwerten erhalten wir:
- $$C=\frac{\rm 1}{\rm 2\it\pi\cdot\sigma_x\cdot\sigma_y\cdot\sqrt{\rm 1 - \rho_{xy}^{\rm 2}}} =\frac{\rm 1}{\rm 2\pi\cdot\rm 1\cdot 1\cdot\sqrt{0.5}}=\frac{\rm 1}{\sqrt{\rm 2}\cdot \pi}\hspace{0.15cm}\underline{\approx \rm 0.225}.$$
- 3. Die allgemeine Gleichung lautet:
- $$\alpha = \frac{\rm 1}{\rm 2}\cdot \rm arctan(\rm 2 \cdot\it \rho_{xy}\cdot \frac{\sigma_x\cdot\sigma_y}{\sigma_x^{\rm 2} - \sigma_y^{\rm 2}}).$$
- Gilt σx = σy und ρxy ≠ 0, so ist der Winkel α stets ±45°. Das Vorzeichen ist abhängig vom Vorzeichen von ρxy. Im vorliegenden Fall gilt α = –45°.
- 4. Für die eingezeichnete Höhenlinie gilt:
- $$f_{xy}(x, y)=\frac{\rm 1}{\sqrt{\rm 2}\cdot \pi}\cdot \rm e^{-(\it x^{\rm 2} + \it y^{\rm 2} + \sqrt{\rm 2}\cdot\it x\cdot\it y)}=\rm 0.2$$
- $$\Rightarrow {\rm e}^{-(\it x^{\rm 2} + \it y^{\rm 2} + \sqrt{\rm 2}\cdot\it x\cdot\it y)} = \rm 0.8885 \hspace{0.5cm}\Rightarrow \hspace{0.5cm} \it x^{\rm 2} + \it y^{\rm 2} + \sqrt{\rm 2}\cdot\it x\cdot\it y = -{\rm ln(0.8885)} \approx\rm 0.118.$$
- Der Winkel der Ellipsenhauptachse ist –45°. Deshalb muss y0 = –x0 gelten. Daraus folgt weiter:
- $$x_{\rm 0}^{\rm 2} + (-x_{\rm 0})^{\rm 2} + \sqrt{\rm 2}\cdot x_{\rm 0}(-x_{\rm 0}) = 0.118$$
- $$\Rightarrow (\rm 2 - \sqrt{\rm 2})\cdot \it x_{\rm 0}^{\rm 2} = {\rm 0.118} \hspace{0.5cm}\Rightarrow \hspace{0.5cm} x_{\rm 0}^{\rm 2} \approx \frac{\rm0.118}{\rm0.585}\approx\rm 0.202; \hspace{0.5cm} x_{\rm 0}\approx\pm\rm 0.450.$$
- Die beiden Schnittpunkte der eingezeichneten Höhenlinien mit der Ellipsenhauptachse liegen somit bei (0.45, –0.45) und (–0.45, 0.45). Der Quotient x0/y0 ist in beiden Fällen –1.
- 5. Vorweg das Ergebnis: Richtig sind die Lösungsvorschläge 2 und 3.
- Mit σy = σx und dem Ergebnis aus (1) gilt für den Winkel der Korrelationsgeraden:
- $$\theta_{y\rightarrow x} = \rm arctan (\it \rho_{\it xy})=\rm arctan(-\frac{\rm 1}{\sqrt{\rm 2}})\approx -\rm 35.3^{\circ}.$$
- Das bedeutet: Die erste Aussage ist falsch und die zweite richtig. Nachfolgend der Beweis für die Richtigkeit der Aussage 3: Löst man die Ellipsengleichung (mit z = 0.118), also
- $$x^{\rm 2}+ y^{\rm 2} +\sqrt{\rm 2}\cdot \it x\cdot \it y - \it z = \rm 0 ,$$
- nach y auf, so erhält man nach Lösung einer quadratischen Gleichung
- $$y_{\rm 1/2}=\frac{\sqrt{\rm 2}}{\rm 2}\it x\pm\sqrt{\frac{x^{\rm 2}}{\rm 2}-x^{\rm 2}+\it z} \hspace{0.5cm}\Rightarrow \hspace{0.5cm} y_{\rm 1/2}=\frac{\it x}{\sqrt{\rm 2}}\pm \sqrt{\it z-\frac{x^{\rm 2}}{\rm 2}}.$$
- Die vertikale Tangente ergibt sich für den Fall, dass die beiden Lösungen y1/2 identisch sind. Das heißt: der Wurzelausdruck muss den Wert 0 ergeben. Die Lösung für positives x lautet dann:
- $$x_{\rm T}=\sqrt{\rm 2\cdot \it z}=\rm \rm 0.485.$$
- Eingesetzt in die Ellipsengleichung erhält man für den y-Wert des Tangentialpunktes:
- $$x_{\rm T}^{\rm 2} + y_{\rm T}^{\rm 2} + \sqrt{\rm 2}\cdot\it x_{\rm T} \cdot y_{\rm T} - \it z = \rm 0 \hspace{0.5cm}\Rightarrow \hspace{0.5cm}\rm 2 \it z + y_{\rm T}^{\rm 2} + \rm 2\sqrt{\it z}\cdot\it y_{\rm T} - \it z = \rm 0$$
- $$\Rightarrow y_{\rm T}^{\rm 2} + \rm 2\sqrt{\it z}\cdot\it y_{\rm T} + \it z = \rm 0 \hspace{0.5cm}\Rightarrow \hspace{0.5cm} (y_{\rm T} + \sqrt{\it z}) = \rm 0\hspace{0.5cm}\Rightarrow \hspace{0.5cm} y_{\rm T} = -\sqrt{\it z} = -0.343.$$
- Daraus ergibt sich:
- $$y_{\rm T}=-\frac{\it x_{\rm T}}{\sqrt{\rm 2}}.$$
- Das bedeutet aber auch: Der Tangentialpunkt (xT, yT) liegt exakt auf der Korrelationsgeraden:
- $$y=K(x)=-{\it x}/{\sqrt{\rm 2}}.$$