Difference between revisions of "Theory of Stochastic Signals/Uniformly Distributed Random Variables"

From LNTwww
Line 12: Line 12:
  
 
Die Grafik zeigt  
 
Die Grafik zeigt  
*links die Wahrscheinlichkeitsdichtefunktion $f_{x}(x)$ und
+
*links die Wahrscheinlichkeitsdichtefunktion $f_{x}(x)$,
 
*rechts die Verteilungsfunktion $F_{x}(r)$
 
*rechts die Verteilungsfunktion $F_{x}(r)$
  
Line 20: Line 20:
  
 
Aus der Grafik und der Definition können folgende Eigenschaften abgeleitet werden:  
 
Aus der Grafik und der Definition können folgende Eigenschaften abgeleitet werden:  
*Die Wahrscheinlichkeitsdichtefunktion (WDF) $f_{x}(x)$ besitzt im Bereich von $x_{\rm min}$ bis $x_{\rm max}$ den konstanten Wert $1/(x_{\rm max} - x_{\rm min})$, wobei an den beiden Bereichsgrenzen für $f_{x}(x)$ jeweils nur der halbe Wert – also der Mittelwert zwischen links- und rechtsseitigem Grenzwert – zu setzen ist.  
+
*Die Wahrscheinlichkeitsdichtefunktion ('''WDF''') besitzt im Bereich von $x_{\rm min}$ bis $x_{\rm max}$ den konstanten Wert $1/(x_{\rm max} - x_{\rm min})$. An den Bereichsgrenzen ist für $f_{x}(x)$ jeweils nur der halbe Wert – also der Mittelwert zwischen links- und rechtsseitigem Grenzwert – zu setzen.  
*Die Verteilungsfunktion (VTF) $F_{x}(r)$ steigt im Bereich von $x_{\rm min}$ bis $x_{\rm max}$ linear von $0$ auf $1$ linear an.   
+
*Die Verteilungsfunktion ('''VTF''') steigt im Bereich von $x_{\rm min}$ bis $x_{\rm max}$ linear von $0$ auf $1$ linear an.   
 
*Mittelwert und Streuung haben bei der Gleichverteilung die folgenden Werte:
 
*Mittelwert und Streuung haben bei der Gleichverteilung die folgenden Werte:
 
:$$m_{\rm 1} = \frac{\it x_ {\rm max} \rm + \it x_{\rm min}}{2},\hspace{0.5cm}
 
:$$m_{\rm 1} = \frac{\it x_ {\rm max} \rm + \it x_{\rm min}}{2},\hspace{0.5cm}
Line 32: Line 32:
 
Die Grafik zeigt zwei Signalverläufe mit gleichförmiger Amplitudenverteilung.  
 
Die Grafik zeigt zwei Signalverläufe mit gleichförmiger Amplitudenverteilung.  
 
*Links ist statistische Unabhängigkeit der einzelnen Abtastwerte vorausgesetzt, das heißt, $x_ν$ kann alle Werte zwischen $x_{\rm min}$ und $x_{\rm max}$ mit gleicher Wahrscheinlichkeit annehmen, und zwar unabhängig von der Vergangenheit $(x_{ν–1}, x_{ν–2}, \hspace{0.1cm}\text{...}).$  
 
*Links ist statistische Unabhängigkeit der einzelnen Abtastwerte vorausgesetzt, das heißt, $x_ν$ kann alle Werte zwischen $x_{\rm min}$ und $x_{\rm max}$ mit gleicher Wahrscheinlichkeit annehmen, und zwar unabhängig von der Vergangenheit $(x_{ν–1}, x_{ν–2}, \hspace{0.1cm}\text{...}).$  
*Beim rechten Signal $y(t)$ ist diese Unabhängigkeit aufeinanderfolgender Signalwerte nicht mehr gegeben. Vielmehr stellt dieses Sägezahnsignal ein deterministisches Signal dar.  
+
*Beim rechten Signal $y(t)$ ist diese Unabhängigkeit aufeinanderfolgender Signalwerte nicht mehr gegeben. Vielmehr stellt dieses Sägezahnsignal ein deterministisches Signal dar.
  
 
[[File:P_ID618__Sto_T_3_4_S2_neu100.png |center|frame|Beispiele gleichverteilter Signale]] }}
 
[[File:P_ID618__Sto_T_3_4_S2_neu100.png |center|frame|Beispiele gleichverteilter Signale]] }}
Line 49: Line 49:
 
'''(2)'''   In der ''Bildverarbeitung & Bildcodierung'' wird häufig vereinfachend mit der Gleichverteilung anstelle der tatsächlichen, meist sehr viel komplizierteren Verteilung des Originalbildes gerechnet, da der Unterschied des Informationsgehaltes zwischen einem ''natürlichen Bild'' und dem auf der Gleichverteilung basierenden Modell relativ gering ist.  
 
'''(2)'''   In der ''Bildverarbeitung & Bildcodierung'' wird häufig vereinfachend mit der Gleichverteilung anstelle der tatsächlichen, meist sehr viel komplizierteren Verteilung des Originalbildes gerechnet, da der Unterschied des Informationsgehaltes zwischen einem ''natürlichen Bild'' und dem auf der Gleichverteilung basierenden Modell relativ gering ist.  
  
'''(3)'''   Für die Modellierung übertragungstechnischer Systeme sind gleichverteilte Zufallsgrößen dagegen die Ausnahme. Ein Beispiel für eine tatsächlich (nahezu) gleichverteilte Zufallsgröße ist die Phase bei kreissymmetrischen Störungen, wie sie beispielsweise bei ''Quadraturamplitudenmodulationsverfahren'' auftreten.  
+
'''(3)'''   Für die Modellierung übertragungstechnischer Systeme sind gleichverteilte Zufallsgrößen dagegen die Ausnahme. Ein Beispiel für eine tatsächlich (nahezu) gleichverteilte Zufallsgröße ist die Phase bei kreissymmetrischen Störungen, wie sie beispielsweise bei ''Quadratur–Amplitudenmodulationsverfahren'' (QAM) auftreten.  
  
 
+
Das interaktive Applet [[Applets:WDF_VTF|WDF, VTF und Momente spezieller Verteilungen]]  berechnet unter anderem
Das interaktive Applet [[Applets:WDF_VTF|WDF, VTF und Momente spezieller Verteilungen]]  berechnet unter Anderem
+
*alle Kenngrößen der Gleichverteilung (diese wird hier als „Rechteck” bezeichnet)  
*alle Kenngrößen der Gleichverteilung (diese ist hier als „Rechteck” bezeichnet)  
 
 
*für beliebige Parameter $x_{\rm min}$ und $x_{\rm max}$.  
 
*für beliebige Parameter $x_{\rm min}$ und $x_{\rm max}$.  
  
Line 59: Line 58:
 
==Erzeugung einer Gleichverteilung mit Pseudo–Noise–Generatoren==
 
==Erzeugung einer Gleichverteilung mit Pseudo–Noise–Generatoren==
 
<br>
 
<br>
Die heute verwendeten Zufallsgeneratoren sind meist '''pseudozufällig''':  
+
{{BlaueBox|TEXT= 
 +
$\text{Definition}$&nbsp; Die heute verwendeten Zufallsgeneratoren sind meist '''pseudozufällig''':  
 
*Das bedeutet, dass die erzeugte Folge als das Ergebnis eines festen Algorithmuses eigentlich deterministisch ist, für den Anwender jedoch aufgrund der großen Periodenlänge $P$ als stochastisch erscheint.  
 
*Das bedeutet, dass die erzeugte Folge als das Ergebnis eines festen Algorithmuses eigentlich deterministisch ist, für den Anwender jedoch aufgrund der großen Periodenlänge $P$ als stochastisch erscheint.  
*Mehr hierzu im Kapitel [[Stochastische_Signaltheorie/Erzeugung_von_diskreten_Zufallsgrößen|Erzeugung von diskreten Zufallsgrößen]].
+
*Mehr hierzu im Kapitel [[Stochastische_Signaltheorie/Erzeugung_von_diskreten_Zufallsgrößen|Erzeugung von diskreten Zufallsgrößen]]. }}
  
  
Für die Systemsimulation haben Pseudo&ndash;Noise (PN)&ndash;Generatoren gegenüber echten Zufallsgeneratoren den entscheidenden Vorteil, dass die erzeugten Zufallsfolgen ohne Speicherung reproduzierbar sind, was  
+
Für die Systemsimulation haben Pseudo&ndash;Noise ('''PN''')&ndash;Generatoren gegenüber echten Zufallsgeneratoren den entscheidenden Vorteil, dass die erzeugten Zufallsfolgen ohne Speicherung reproduzierbar sind, was  
 
*zum einen den Vergleich verschiedener Systemmodelle ermöglicht, und  
 
*zum einen den Vergleich verschiedener Systemmodelle ermöglicht, und  
 
*auch die Fehlersuche wesentlich erleichtert.  
 
*auch die Fehlersuche wesentlich erleichtert.  
Line 71: Line 71:
 
Ein Zufallsgenerator sollte dabei folgende Kriterien erfüllen:  
 
Ein Zufallsgenerator sollte dabei folgende Kriterien erfüllen:  
  
'''(1)''' &nbsp; Die Zufallsgrößen $x_ν$ einer generierten Folge sollten mit sehr guter Näherung gleichverteilt sein. Bei wertdiskreter Darstellung an einem Rechner erfordert dies unter Anderem eine hinreichend ''hohe Bitauflösung,'' zum Beispiel mit 32 oder 64 Bit pro Abtastwert.  
+
'''(1)''' &nbsp; Die Zufallsgrößen $x_ν$ einer generierten Folge sollten mit sehr guter Näherung gleichverteilt sein. Bei wertdiskreter Darstellung an einem Rechner erfordert dies unter anderem eine hinreichend ''hohe Bitauflösung,'' zum Beispiel mit 32 oder 64 Bit pro Abtastwert.
 +
 
 +
'''(2)''' &nbsp; Bildet man aus der sequentiellen Zufallsfolge $〈x_ν〉$ jeweils nichtüberlappende Paare von Zufallsgrößen, beispielsweise $(x_ν, x_{ν+1})$, $(x_{ν+2}$, $x_{ν+3})$, ... , so sollten diese ''Tupel''&nbsp; in einer zweidimensionalen Darstellung innerhalb eines Quadrates ebenfalls gleichverteilt sein.  
  
'''(2)''' &nbsp; Bildet man aus der sequentiellen Zufallsfolge $〈x_ν〉$ jeweils nichtüberlappende Paare von Zufallsgrößen, beispielsweise $(x_ν, x_{ν+1}), (x_{ν+2}, x_{ν+3})$, ... , so sollten diese ''Tupel'' in einer zweidimensionalen Darstellung innerhalb eines Quadrates ebenfalls gleichverteilt sein.  
+
'''(3)''' &nbsp; Bildet man aus der sequentiellen Folge $〈x_ν〉$ nicht überlappende $n$&ndash;''Tupel'' von Zufallsgrößen &nbsp; ⇒ &nbsp; $(x_ν$, ... , $x_{ν+n–1})$, $(x_{ν+n}$, ... , $x_{ν+2n–1})$ usw., so sollten auch diese innerhalb eines $n$&ndash;dimensionalen Würfels möglichst die Gleichverteilung ergeben.  
  
'''(3)''' &nbsp; Bildet man aus der sequentiellen Folge $〈x_ν〉$ nicht überlappende $n$-''Tupel'' von Zufallsgrößen  ⇒  $(x_ν$, ... , $x_{ν+n–1}), (x_{ν+n}$, ... , $x_{ν+2n–1})$ usw., so sollten auch diese innerhalb eines $n$&ndash;dimensionalen Würfels möglichst die Gleichverteilung ergeben.
 
  
 
''Anmerkung:''
 
''Anmerkung:''
*Die erste Forderung bezieht sich ausschließlich auf die ''Amplitudenverteilung'' (WDF) und ist im Allgemeinen leichter zu erfüllen.  
+
*Die erste Forderung bezieht sich ausschließlich auf die ''Amplitudenverteilung''&nbsp; (WDF) und ist im Allgemeinen leichter zu erfüllen.  
 
*Die beiden weiteren Forderungen sollen eine „ausreichende Zufälligkeit” der Folge gewährleisten.  
 
*Die beiden weiteren Forderungen sollen eine „ausreichende Zufälligkeit” der Folge gewährleisten.  
 
*Sie betreffen die statistische Unabhängigkeit aufeinander folgender Zufallswerte.  
 
*Sie betreffen die statistische Unabhängigkeit aufeinander folgender Zufallswerte.  
Line 89: Line 90:
 
'''(1)''' &nbsp;  Diese Zufallsgeneratoren basieren auf der sukzessiven Manipulation einer Integervariablen $k$. Geschieht die Zahlendarstellung im Rechner mit $L$ Bit, so nimmt diese Variable bei geeigneter Behandlung des Vorzeichenbits alle Werte zwischen $1$ und $2^{L − 1}$ jeweils genau einmal an.  
 
'''(1)''' &nbsp;  Diese Zufallsgeneratoren basieren auf der sukzessiven Manipulation einer Integervariablen $k$. Geschieht die Zahlendarstellung im Rechner mit $L$ Bit, so nimmt diese Variable bei geeigneter Behandlung des Vorzeichenbits alle Werte zwischen $1$ und $2^{L − 1}$ jeweils genau einmal an.  
  
'''(2)''' &nbsp;  Die hieraus abgeleitete Zufallsgröße  
+
'''(2)''' &nbsp;  Die hieraus abgeleitete Zufallsgröße $x={k}/{\rm 2^{\it L - \rm 1}}$ ist ebenfalls diskret (mit Stufenzahl $M = 2^{L– 1})$:
:$$x={k}/({\rm 2^{\it L - \rm 1}}) = k\cdot \Delta x \in  \{\Delta x, \hspace{0.05cm}2\cdot \Delta x,\hspace{0.05cm}\text{ ...}\hspace{0.05cm} , \hspace{0.05cm}1-\Delta x,\hspace{0.05cm} 1\}$$
+
:$$x={k}/{\rm 2^{\it L - \rm 1}} = k\cdot \Delta x \in  \{\Delta x, \hspace{0.05cm}2\cdot \Delta x,\hspace{0.05cm}\text{ ...}\hspace{0.05cm} , \hspace{0.05cm}1-\Delta x,\hspace{0.05cm} 1\}.$$
ist demnach ebenfalls diskret (mit Stufenzahl $M = 2^{L– 1})$. Ist die Bitanzahl $L$ hinreichend groß, so ist der Abstand $Δx = 1/(2^{L– 1})$ zwischen zwei möglichen Werten sehr klein, und man kann $x$ im Rahmen der Simulationsgenauigkeit durchaus als eine kontinuierliche Zufallsgröße interpretieren.  
+
Ist die Bitanzahl $L$ hinreichend groß, so ist der Abstand $Δx = 1/2^{L– 1}$ zwischen zwei möglichen Werten sehr klein, und man kann $x$ im Rahmen der Simulationsgenauigkeit durchaus als eine kontinuierliche Zufallsgröße interpretieren.  
  
 
'''(3)''' &nbsp; Die rekursive Generierungsvorschrift eines solchen ''Multiplicative Congruential Generators'' lautet:  
 
'''(3)''' &nbsp; Die rekursive Generierungsvorschrift eines solchen ''Multiplicative Congruential Generators'' lautet:  
Line 98: Line 99:
 
'''(4)''' &nbsp; Die statistischen Eigenschaften der Folge hängen entscheidend von den Parametern $a$ und $m$ ab. Der Startwert $k_0$ hat dagegen für die Statistik eine eher untergeordnete Bedeutung.  
 
'''(4)''' &nbsp; Die statistischen Eigenschaften der Folge hängen entscheidend von den Parametern $a$ und $m$ ab. Der Startwert $k_0$ hat dagegen für die Statistik eine eher untergeordnete Bedeutung.  
  
'''(5)''' &nbsp; Die besten Ergebnisse erzielt man mit der Basis $m =2^l–1$, wobei $l$ eine beliebige natürliche Zahl angibt. Weit verbreitet ist bei Rechnern mit 32 Bit-Architektur und einem Vorzeichenbit die Basis $m = 2^{31} – 1 = 2'147'483'647.$ Ein entsprechender Algorithmus lautet:  
+
'''(5)''' &nbsp; Die besten Ergebnisse erzielt man mit der Basis $m =2\hspace{0.05cm}^l–1$, wobei $l$ eine beliebige natürliche Zahl angibt. Weit verbreitet ist bei Rechnern mit 32 Bit-Architektur und einem Vorzeichenbit die Basis $m = 2^{31} – 1 = 2\hspace{0.08cm}147\hspace{0.08cm}483\hspace{0.08cm}647.$ Ein entsprechender Algorithmus lautet:  
 
:$$k_\nu=(16807\cdot k_{\nu-1})\hspace{0.1cm} \rm mod\hspace{0.1cm}(2^{31}-1).$$
 
:$$k_\nu=(16807\cdot k_{\nu-1})\hspace{0.1cm} \rm mod\hspace{0.1cm}(2^{31}-1).$$
  
'''(6)''' &nbsp; Für einen solchen Generator ist nur der Startwert $k_0 = 0$ nicht erlaubt. Für alle anderen Startwerte beträgt die Periodendauer $P = 2^{31} − 2.$
+
'''(6)''' &nbsp; Für einen solchen Generator ist nur der Startwert $k_0 = 0$ nicht erlaubt. Für $k_0 \ne 0$ beträgt die Periodendauer $P = 2^{31} − 2.$
 
 
  
 
[[File:P_ID920__Sto_T_3_4_S5.png |right|frame| Multiplicative Congruental Generator (C-Programm)]]
 
[[File:P_ID920__Sto_T_3_4_S5.png |right|frame| Multiplicative Congruental Generator (C-Programm)]]
 
{{GraueBox|TEXT=   
 
{{GraueBox|TEXT=   
$\text{Beispiel 2:}$&nbsp;  
+
$\text{Beispiel 2:}$&nbsp; Wir analysieren den oben beschriebenen ''Multiplicative Congruental Generator'' genauer:
Der hier angegebene Algorithmus kann auf einem 32 Bit-Rechner allerdings nicht direkt implementiert werden, da das Ergebnis der Multiplikation bis zu 46 Bit beansprucht.  
+
*Den Algorithmus kann man allerdings auf einem 32 Bit&ndash;Rechner nicht direkt implementieren, da das Multiplikationsergebnis bis zu 46 Bit benötigt.  
 
+
*Er kann aber so abgewandelt werden, dass zu keinem Zeitpunkt der Berechnung der 32 Bit&ndash;Integerzahlenbereich überschritten wird.  
 
 
*Er kann aber so abgewandelt werden, dass zu keinem Zeitpunkt der Berechnung der Integerzahlenbereich von 32 Bit überschritten wird.  
 
 
 
 
 
 
*Das so modifizierte C-Programm $\text{uniform( )}$ ist rechts angegeben.}}
 
*Das so modifizierte C-Programm $\text{uniform( )}$ ist rechts angegeben.}}
  

Revision as of 07:16, 9 August 2018

Allgemeine Beschreibung und Definition


WDF und VTF der Gleichverteilung

$\text{Definition:}$  Eine Zufallsgröße $x$ bezeichnet man als gleichverteilt, wenn sie nur Werte im Bereich von $x_{\rm min}$ bis $x_{\rm max}$ annehmen kann, und zwar mit gleicher Wahrscheinlichkeit.

Die Grafik zeigt

  • links die Wahrscheinlichkeitsdichtefunktion $f_{x}(x)$,
  • rechts die Verteilungsfunktion $F_{x}(r)$


einer gleichverteilten Zufallsgröße $x$.


Aus der Grafik und der Definition können folgende Eigenschaften abgeleitet werden:

  • Die Wahrscheinlichkeitsdichtefunktion (WDF) besitzt im Bereich von $x_{\rm min}$ bis $x_{\rm max}$ den konstanten Wert $1/(x_{\rm max} - x_{\rm min})$. An den Bereichsgrenzen ist für $f_{x}(x)$ jeweils nur der halbe Wert – also der Mittelwert zwischen links- und rechtsseitigem Grenzwert – zu setzen.
  • Die Verteilungsfunktion (VTF) steigt im Bereich von $x_{\rm min}$ bis $x_{\rm max}$ linear von $0$ auf $1$ linear an.
  • Mittelwert und Streuung haben bei der Gleichverteilung die folgenden Werte:
$$m_{\rm 1} = \frac{\it x_ {\rm max} \rm + \it x_{\rm min}}{2},\hspace{0.5cm} \sigma = \frac{\it x_{\rm max} - \it x_{\rm min}}{2 \sqrt{3}}.$$
  • Bei symmetrischer WDF   ⇒   $x_{\rm min} = -x_{\rm max}$ erhält man als Sonderfall den Mittelwert $m_1 = 0$ und die Varianz $σ^2 = x_{\rm max}^2/3.$


$\text{Beispiel 1:}$  Die Grafik zeigt zwei Signalverläufe mit gleichförmiger Amplitudenverteilung.

  • Links ist statistische Unabhängigkeit der einzelnen Abtastwerte vorausgesetzt, das heißt, $x_ν$ kann alle Werte zwischen $x_{\rm min}$ und $x_{\rm max}$ mit gleicher Wahrscheinlichkeit annehmen, und zwar unabhängig von der Vergangenheit $(x_{ν–1}, x_{ν–2}, \hspace{0.1cm}\text{...}).$
  • Beim rechten Signal $y(t)$ ist diese Unabhängigkeit aufeinanderfolgender Signalwerte nicht mehr gegeben. Vielmehr stellt dieses Sägezahnsignal ein deterministisches Signal dar.
Beispiele gleichverteilter Signale

Bedeutung der Gleichverteilung für die Nachrichtentechnik


Die Bedeutung gleichverteilter Zufallsgrößen für die Informations– und Kommunikationstechnik ist darauf zurückzuführen, dass diese WDF–Form aus Sicht der Informationstheorie unter der Nebenbedingung Spitzenwertbegrenzung ein Optimum darstellt:

  • Mit keiner anderen Verteilung als der Gleichverteilung erreicht man unter dieser Voraussetzung eine größere differentielle Entropie.
  • Mit dieser Thematik beschäftigt sich das Kapitel Differentielle Entropie im Buch „Informationstheorie”.


Daneben sind unter Anderem noch folgende Punkte zu nennen:

(1)   Die Bedeutung der Gleichverteiltung für die Simulation nachrichtentechnischer Systeme ist darauf zurückzuführen, dass man entsprechende „Pseudo–Zufallsgeneratoren” relativ einfach realisieren kann, und dass sich daraus andere Verteilungen wie zum Beispiel die Gaußverteilung und die Exponentialverteilung leicht ableiten lassen.

(2)   In der Bildverarbeitung & Bildcodierung wird häufig vereinfachend mit der Gleichverteilung anstelle der tatsächlichen, meist sehr viel komplizierteren Verteilung des Originalbildes gerechnet, da der Unterschied des Informationsgehaltes zwischen einem natürlichen Bild und dem auf der Gleichverteilung basierenden Modell relativ gering ist.

(3)   Für die Modellierung übertragungstechnischer Systeme sind gleichverteilte Zufallsgrößen dagegen die Ausnahme. Ein Beispiel für eine tatsächlich (nahezu) gleichverteilte Zufallsgröße ist die Phase bei kreissymmetrischen Störungen, wie sie beispielsweise bei Quadratur–Amplitudenmodulationsverfahren (QAM) auftreten.

Das interaktive Applet WDF, VTF und Momente spezieller Verteilungen berechnet unter anderem

  • alle Kenngrößen der Gleichverteilung (diese wird hier als „Rechteck” bezeichnet)
  • für beliebige Parameter $x_{\rm min}$ und $x_{\rm max}$.


Erzeugung einer Gleichverteilung mit Pseudo–Noise–Generatoren


$\text{Definition}$  Die heute verwendeten Zufallsgeneratoren sind meist pseudozufällig:

  • Das bedeutet, dass die erzeugte Folge als das Ergebnis eines festen Algorithmuses eigentlich deterministisch ist, für den Anwender jedoch aufgrund der großen Periodenlänge $P$ als stochastisch erscheint.
  • Mehr hierzu im Kapitel Erzeugung von diskreten Zufallsgrößen.


Für die Systemsimulation haben Pseudo–Noise (PN)–Generatoren gegenüber echten Zufallsgeneratoren den entscheidenden Vorteil, dass die erzeugten Zufallsfolgen ohne Speicherung reproduzierbar sind, was

  • zum einen den Vergleich verschiedener Systemmodelle ermöglicht, und
  • auch die Fehlersuche wesentlich erleichtert.


Ein Zufallsgenerator sollte dabei folgende Kriterien erfüllen:

(1)   Die Zufallsgrößen $x_ν$ einer generierten Folge sollten mit sehr guter Näherung gleichverteilt sein. Bei wertdiskreter Darstellung an einem Rechner erfordert dies unter anderem eine hinreichend hohe Bitauflösung, zum Beispiel mit 32 oder 64 Bit pro Abtastwert.

(2)   Bildet man aus der sequentiellen Zufallsfolge $〈x_ν〉$ jeweils nichtüberlappende Paare von Zufallsgrößen, beispielsweise $(x_ν, x_{ν+1})$, $(x_{ν+2}$, $x_{ν+3})$, ... , so sollten diese Tupel  in einer zweidimensionalen Darstellung innerhalb eines Quadrates ebenfalls gleichverteilt sein.

(3)   Bildet man aus der sequentiellen Folge $〈x_ν〉$ nicht überlappende $n$–Tupel von Zufallsgrößen   ⇒   $(x_ν$, ... , $x_{ν+n–1})$, $(x_{ν+n}$, ... , $x_{ν+2n–1})$ usw., so sollten auch diese innerhalb eines $n$–dimensionalen Würfels möglichst die Gleichverteilung ergeben.


Anmerkung:

  • Die erste Forderung bezieht sich ausschließlich auf die Amplitudenverteilung  (WDF) und ist im Allgemeinen leichter zu erfüllen.
  • Die beiden weiteren Forderungen sollen eine „ausreichende Zufälligkeit” der Folge gewährleisten.
  • Sie betreffen die statistische Unabhängigkeit aufeinander folgender Zufallswerte.


Multiplicative Congruental Generator


$\text{Multiplicative Congruental Generator}$ ist das bekannteste Verfahren zur Erzeugung einer Folge $〈 x_\nu 〉$ mit gleichverteilten Werten $ x_\nu$ zwischen $0$ und $1$. Diese Methode wird hier stichpunktartig angegeben:

(1)   Diese Zufallsgeneratoren basieren auf der sukzessiven Manipulation einer Integervariablen $k$. Geschieht die Zahlendarstellung im Rechner mit $L$ Bit, so nimmt diese Variable bei geeigneter Behandlung des Vorzeichenbits alle Werte zwischen $1$ und $2^{L − 1}$ jeweils genau einmal an.

(2)   Die hieraus abgeleitete Zufallsgröße $x={k}/{\rm 2^{\it L - \rm 1}}$ ist ebenfalls diskret (mit Stufenzahl $M = 2^{L– 1})$:

$$x={k}/{\rm 2^{\it L - \rm 1}} = k\cdot \Delta x \in \{\Delta x, \hspace{0.05cm}2\cdot \Delta x,\hspace{0.05cm}\text{ ...}\hspace{0.05cm} , \hspace{0.05cm}1-\Delta x,\hspace{0.05cm} 1\}.$$

Ist die Bitanzahl $L$ hinreichend groß, so ist der Abstand $Δx = 1/2^{L– 1}$ zwischen zwei möglichen Werten sehr klein, und man kann $x$ im Rahmen der Simulationsgenauigkeit durchaus als eine kontinuierliche Zufallsgröße interpretieren.

(3)   Die rekursive Generierungsvorschrift eines solchen Multiplicative Congruential Generators lautet:

$$k_\nu=(a\cdot k_{\nu-1})\hspace{0.1cm} \rm mod \hspace{0.1cm} \it m.$$

(4)   Die statistischen Eigenschaften der Folge hängen entscheidend von den Parametern $a$ und $m$ ab. Der Startwert $k_0$ hat dagegen für die Statistik eine eher untergeordnete Bedeutung.

(5)   Die besten Ergebnisse erzielt man mit der Basis $m =2\hspace{0.05cm}^l–1$, wobei $l$ eine beliebige natürliche Zahl angibt. Weit verbreitet ist bei Rechnern mit 32 Bit-Architektur und einem Vorzeichenbit die Basis $m = 2^{31} – 1 = 2\hspace{0.08cm}147\hspace{0.08cm}483\hspace{0.08cm}647.$ Ein entsprechender Algorithmus lautet:

$$k_\nu=(16807\cdot k_{\nu-1})\hspace{0.1cm} \rm mod\hspace{0.1cm}(2^{31}-1).$$

(6)   Für einen solchen Generator ist nur der Startwert $k_0 = 0$ nicht erlaubt. Für $k_0 \ne 0$ beträgt die Periodendauer $P = 2^{31} − 2.$

Multiplicative Congruental Generator (C-Programm)

$\text{Beispiel 2:}$  Wir analysieren den oben beschriebenen Multiplicative Congruental Generator genauer:

  • Den Algorithmus kann man allerdings auf einem 32 Bit–Rechner nicht direkt implementieren, da das Multiplikationsergebnis bis zu 46 Bit benötigt.
  • Er kann aber so abgewandelt werden, dass zu keinem Zeitpunkt der Berechnung der 32 Bit–Integerzahlenbereich überschritten wird.
  • Das so modifizierte C-Programm $\text{uniform( )}$ ist rechts angegeben.

Aufgaben zum Kapitel


Aufgabe 3.5: Dreieck– und Trapezsignal

Aufgabe 3.5Z: Antennengebiete