Wie können sortierte, gleichmäßig verteilte Werte in einem Intervall effizient generiert werden?

Angenommen, ich möchte aus dem Intervall eine Reihe von Zufallszahlen generieren (a, b). Die generierte Sequenz sollte auch die Eigenschaft haben, dass sie sortiert ist. Ich kann mir zwei Möglichkeiten vorstellen, um dies zu erreichen.

Sei ndie Länge der zu erzeugenden Sequenz.

1. Algorithmus:

Let `offset = floor((b - a) / n)`
for i = 1 up to n:
   generate a random number r_i from (a, a+offset)
   a = a + offset
   add r_i to the sequence r

2. Algorithmus:

for i = 1 up to n:
    generate a random number s_i from (a, b)
    add s_i to the sequence s
sort(r)

Meine Frage ist, erzeugt Algorithmus 1 Sequenzen, die so gut sind wie die von Algorithmus 2 erzeugten?

Der erste Algorithmus schlägt aus zwei Gründen schlecht fehl :

1. Wenn Sie den Boden von kann dies drastisch reduziert werden. In der Tat, wenn b - a < n ist , ist es Null, was Ihnen eine Menge gibt, deren Werte alle gleich sind!$(a-b)/n$$b-a \lt n$

2. Wenn Sie nicht das Wort ergreifen, sind die resultierenden Werte zu gleichmäßig verteilt. Zum Beispiel gibt es in jeder einfachen Zufallsstichprobe von iid einheitlichen Variablen (etwa zwischen a = 0 und b = 1 ) eine ( 1 - 1 / n ) n ≈ 1 / e ≈ 37 % Wahrscheinlichkeit, dass die größte nicht sein wird im oberen Intervall von 1 - 1 / n bis 1 . Mit Algorithmus 1 gibt es eine $n$$a=0$$b=1$$(1-1/n)^n\approx 1/e\approx 37\%$$1-1/n$$1$$100\%$Chance, dass das Maximum in diesem Intervall liegt. Für einige Zwecke ist diese Supergleichmäßigkeit gut, aber im Allgemeinen ist es ein schrecklicher Fehler, weil (a) viele Statistiken ruiniert werden, aber (b) es sehr schwierig sein kann, festzustellen, warum.

3. Wenn Sie das Sortieren vermeiden möchten, generieren Sie stattdessen unabhängige exponentiell verteilte Variablen. Normalisieren Sie ihre kumulative Summe auf den Bereich ( 0 , 1 ), indem Sie durch die Summe dividieren. Löschen Sie den größten Wert (der immer 1 sein wird ). Skalieren Sie auf den Bereich ( a , b ) .$n+1$$(0,1)$$1$$(a,b)$

Histogramme aller drei Algorithmen werden angezeigt. (Jedes zeigt die kumulativen Ergebnisse von unabhängigen Sätzen mit jeweils n = 100 Werten.) Das Fehlen einer sichtbaren Variation im Histogramm für Algorithmus 1 zeigt das Problem dort. Die Variation der beiden anderen Algorithmen ist genau das, was zu erwarten ist - und was Sie$1000$$n=100$ von einem Zufallszahlengenerator benötigen .

Weitere (amüsante) Möglichkeiten zum Simulieren unabhängiger gleichmäßiger Variablen finden Sie unter Simulieren von Zeichnungen aus einer gleichmäßigen Verteilung mithilfe von Zeichnungen aus einer Normalverteilung .

Hier ist der RCode, der die Figur erzeugt hat.

b <- 1
a <- 0
n <- 100
n.iter <- 1e3

offset <- (b-a)/n
as <- seq(a, by=offset, length.out=n)
sim.1 <- matrix(runif(n.iter*n, as, as+offset), nrow=n)
sim.2 <- apply(matrix(runif(n.iter*n, a, b), nrow=n), 2, sort)
sim.3 <- apply(matrix(rexp(n.iter*(n+1)), nrow=n+1), 2, function(x) {
  a + (b-a) * cumsum(x)[-(n+1)] / sum(x)
})

par(mfrow=c(1,3))
hist(sim.1, main="Algorithm 1")
hist(sim.2, main="Algorithm 2")
hist(sim.3, main="Exponential")

Der erste Algorithmus erzeugt zu gleichmäßig verteilte Zahlen

Siehe auch Reihen mit geringer Diskrepanz .

$[0;1]$

(Wie bereits ausgeführt, ist dies eine gewünschte Eigenschaft zB für Schichtung sein kann. Low-Diskrepanz Serien wie Halton und Sobel haben ihre Fälle verwenden.)

Ein richtiger, aber teurer Ansatz (für echte Werte)

... soll Beta-verteilte Zufallszahlen verwenden. Die Rangordnungsstatistik der Gleichverteilung ist Beta-verteilt. Sie können dies verwenden, um zufällig die kleinste , dann die zweitkleinste, ... Wiederholung zu zeichnen .

$[0;1]$$\text{Beta}[1,n]$$n$$1-X\sim \text{Beta}[n, 1]$$-\ln (1-X)\sim \text{Exponential}[n]$$\frac{-\ln(U[0;1])}{n}$

Was den folgenden Algorithmus ergibt:

x = a
for i in range(n, 0, -1):
    x += (b-x) * (1 - pow(rand(), 1. / i))
    result.append(x) 

Es kann numerische Instabilitäten geben, und das Berechnen powund Teilen für jedes Objekt kann sich als langsamer als das Sortieren herausstellen.

Für ganzzahlige Werte müssen Sie möglicherweise eine andere Verteilung verwenden.

Das Sortieren ist unglaublich billig, verwenden Sie es also einfach

$O(n \log n)$

Es hängt auch davon ab, was Sie mit den Zufallszahlen machen. Bei numerischen Integrationsproblemen würde Methode 1 (wenn sie durch Entfernen des Bodenoperators korrigiert wird) eine überlegene Punktmenge erzeugen. Was Sie tun, ist eine Form der geschichteten Probenahme und hat den Vorteil, dass Verklumpungen vermieden werden. Es ist beispielsweise unmöglich, alle Ihre Werte im Bereich von 0- (ba) / n zu erhalten. Für andere Anwendungen kann dies jedoch sehr schlecht sein. Dies hängt davon ab, was Sie damit tun möchten.