Eine Verteilungsfunktion für rationale Zahlen im Einheitsintervall berechnen 

Bringen Sie die rationalen Zahlen zwischen Null und eins in eine Reihenfolge entsprechend ihrer Binärdarstellung.

Die gesamte Wolfram-Language Eingabe zeigenEingabe verbergen

In[1]:=

X RankedRational[n_Integer?Positive] := Block[{dig = IntegerDigits[n, 2], avec, cf}, avec = Reverse[Pick[Range[Length[dig] - 1, 0, -1], dig, 1]]; cf = ConstantArray[0, Length[avec] + 1]; cf[[3 ;; -1]] = Differences[avec]; cf[[-1]] += 1; cf[[2]] = avec[[1]] + 1; cf // FromContinuedFraction ]

In[2]:=

X Table[RankedRational[n], {n, 1, 16}]

Out[2]=

Die kumulative Verteilungsfunktion rationaler Zahlen in dieser Folge ist die MinkowskiQuestionMark-Funktion.

In[3]:=

X cdf[x_] = CDF[EmpiricalDistribution[Table[RankedRational[n], {n, 1024}]], x];

Die gesamte Wolfram-Language Eingabe zeigenEingabe verbergen

In[4]:=

X Plot[{MinkowskiQuestionMark[x], cdf[x] }, {x, 0, 1}, PlotStyle -> {Red, Directive[Thick, Black, Dashing[.04]]}, PlotTheme -> "Detailed", AspectRatio -> 1, ImageSize -> Medium]

Out[4]=