Trochę spóźniłem się na imprezę, ale musiałem wdrożyć ogólne rozwiązanie i okazało się, że żadne z rozwiązań nie jest w stanie zaspokoić moich potrzeb.
Przyjęte rozwiązanie jest dobre dla małych zasięgów; jednakże maximum - minimum
może to być nieskończoność dla dużych zakresów. Tak więc poprawioną wersją może być ta wersja:
public static double NextDoubleLinear(this Random random, double minValue, double maxValue)
{
// TODO: some validation here...
double sample = random.NextDouble();
return (maxValue * sample) + (minValue * (1d - sample));
}
To ładnie generuje liczby losowe, nawet między double.MinValue
a double.MaxValue
. Ale to wprowadza kolejny „problem”, który ładnie przedstawiono w tym poście : jeśli użyjemy tak dużych zakresów, wartości mogą wydawać się zbyt „nienaturalne”. Na przykład po wygenerowaniu 10 000 losowych podwójnych wartości z zakresu od 0 do double.MaxValue
wszystkich wartości znajdowały się w przedziale od 2,9579E + 304 do 1,7976E + 308.
Stworzyłem więc również inną wersję, która generuje liczby w skali logarytmicznej:
public static double NextDoubleLogarithmic(this Random random, double minValue, double maxValue)
{
// TODO: some validation here...
bool posAndNeg = minValue < 0d && maxValue > 0d;
double minAbs = Math.Min(Math.Abs(minValue), Math.Abs(maxValue));
double maxAbs = Math.Max(Math.Abs(minValue), Math.Abs(maxValue));
int sign;
if (!posAndNeg)
sign = minValue < 0d ? -1 : 1;
else
{
// if both negative and positive results are expected we select the sign based on the size of the ranges
double sample = random.NextDouble();
var rate = minAbs / maxAbs;
var absMinValue = Math.Abs(minValue);
bool isNeg = absMinValue <= maxValue ? rate / 2d > sample : rate / 2d < sample;
sign = isNeg ? -1 : 1;
// now adjusting the limits for 0..[selected range]
minAbs = 0d;
maxAbs = isNeg ? absMinValue : Math.Abs(maxValue);
}
// Possible double exponents are -1022..1023 but we don't generate too small exponents for big ranges because
// that would cause too many almost zero results, which are much smaller than the original NextDouble values.
double minExponent = minAbs == 0d ? -16d : Math.Log(minAbs, 2d);
double maxExponent = Math.Log(maxAbs, 2d);
if (minExponent == maxExponent)
return minValue;
// We decrease exponents only if the given range is already small. Even lower than -1022 is no problem, the result may be 0
if (maxExponent < minExponent)
minExponent = maxExponent - 4;
double result = sign * Math.Pow(2d, NextDoubleLinear(random, minExponent, maxExponent));
// protecting ourselves against inaccurate calculations; however, in practice result is always in range.
return result < minValue ? minValue : (result > maxValue ? maxValue : result);
}
Niektóre testy:
Oto posortowane wyniki generowania 10000 losowych liczb podwójnych z przedziału od 0 i Double.MaxValue
przy obu strategiach. Wyniki wyświetlane są w skali logarytmicznej:
Chociaż liniowe wartości losowe wydają się na pierwszy rzut oka błędne, statystyki pokazują, że żadna z nich nie jest „lepsza” od drugiej: nawet strategia liniowa ma równomierny rozkład, a średnia różnica między wartościami jest prawie taka sama w przypadku obu strategii .
Gra z różnymi zakresami pokazała mi, że strategia liniowa staje się „rozsądna” w zakresie od 0 ushort.MaxValue
do „rozsądnej” wartości minimalnej 10,78294704 (dla ulong
zakresu minimalna wartość to 3,03518E + 15 int
;: 353341). Oto te same wyniki obu strategii wyświetlane w różnych skalach:
Edytować:
Niedawno uczyniłem moje biblioteki open source, zapraszam do obejrzenia RandomExtensions.NextDouble
metody z pełną walidacją.