Das halte ich für falsch!hikaru hat geschrieben:Die Idee ist algorithmisch gut, aber um ehrlich zu sein finde ich den Weg als Lösung der Aufgabe hässlich.
Den Effekt gibt es, aber er ist hier nicht sonderlich groß.rvb hat geschrieben:05.01.2019 11:50:55Ein einfaches Beispiel: angenommen, du hättest einen perfekten Zufallsgenerator, der Zahlen zwischen 0 und 8 absolut gleichverteilt liefert, aber möchtest eine Zufallszahl im Intervall [0,5] haben. Wenn du dazu einfach das Ergebnis des Zufallsgenerators modulo 6 nimmst, dann kriegst du zwar eine Zufallszahl im Intervall [0,5], aber die landet mit 2/3 Wahrscheinlichkeit im Intervall [0,2] und nur mit 1/3 Wahrscheinlichkeit im Intervall [3,5].
Das beste /dev/urandom hilft also nichts, wenn man hinterher ein Modulo dranhängt.
Ist richtig.
Der ist nur ±1, wenn der passende Winkel dazukommt. Dass bei Winkel nahe 0 viele x nahe Radius sind, muss so sein, weil das auf ein großes Stück in y-Richtung verteilt wird.hikaru hat geschrieben:05.01.2019 12:33:45Der Radius den heinz mit 2000 definiert hat, wird hier auf der Hochachse auf ±1 normiert und die Werte sollten durch den Zufallsgenerator gleichmäßig in diesem Wertebereich verteilt sein.
Bei A und B läuft der Kreisbogen parallel zur y-Achse. Bei C und D läuft er π/4 dazu. Das erklärt die Wurzel 2. Auf gleich langem Kreisbogen müssen gleich viele Punkte landen.hikaru hat geschrieben:05.01.2019 12:33:45Um ein halbwegs anschauliches Beispiel zu bringen, schauen wir uns mal 4 Punkte an, die alle auf dem Rarius r liegen:
A: φ=0°; B: φ=1°; C: φ=45°; D: φ=46°
Der Kosinus von A und B unterscheidet sich erst bei 10^-4, während C und D sich schon bei 10^-2 unterscheiden. Der Abstand zwischen A und B ist nahe 1, da sie kartesisch fast die selbe x-Koordinate haben. C und D unterscheiden sich dagegen in x unyd y um jeweils 1, was dazu führt, dass ihr Abstand nicht 1 sondern Wurzel 2 ist.
Genauso ist es, und da man ja sonst nichts zu tun hat...smutbert hat geschrieben:Ich liebe das Debianforum für solche Threads!!
40575Code: Alles auswählen
c++ z2z.cc -o z2z -I/usr/include/SDL -lSDL_gfx -lSDL -WallNein, wenn ich die Punkte, die du nach nopaste kopiert hast, mit meinem Testprogramm auftrage, kommt das raus:
Danke fuers ausprobieren!MSfree hat geschrieben:06.01.2019 19:33:27Nein, wenn ich die Punkte, die du nach nopaste kopiert hast, mit meinem Testprogramm auftrage, kommt das raus:
Das Kreuz ist also in den Daten vorhanden. In der Folge muß natürlich auch povray ein Kreuz anzeigen.
40576Code: Alles auswählen
translate <sin(W)*(-200+400*rand(Z1)),cos(W)*(-200+400*rand(Z1)),0>Code: Alles auswählen
#declare Z1=seed(0);Das ist es leider auch nicht...MSfree hat geschrieben:06.01.2019 21:34:11Es ist denkbar, daß mit seed(0) das Kreuz besonders auffällig ist, mit anderen seeds aber nicht.
40569 nochmal genauer angesehen.Schade.
Die Funktion malt Punkte und sie ist für den Allgemeinfall ausgelegt.Welchen Zweck hat Deine Funktion void DrawDots? Wozu skalierst Du da was...
Die Skalierung, die sich in meinem Fall ergibt,ist 1.0, denn die Punkte haben einen Wertebereich von -200 bis +200 und die Zeichenfläche ist 400 x 400 Pixel. Wie gesagt, die Skalierung ist in diesem Fall unnötig, aber für den Allgemeinfall ausgelegt.Wenn Du die Koordinaten in x(i),y(i) direckt "Zeichnen" wuerdest, wuede da das gleiche Bild herauskommen? Also ohne Kreuz?
Ich verstehe noch nicht, warum das zu der Kreuzbildung führt. Wenn man zwei verschiedene Zufallszahlen für den Radius in x- und y-Richtung wählt, liegt der Punkt zwar auf einer Ellipse statt auf einem Kreis, aber gemittelt über 10000 Punkte würde ich trotzdem eine gleichmäßige Verteilung* erwarten.heinz hat geschrieben:06.01.2019 22:15:20Ich habe den Radius nicht einmal berechnet und dann fuer sin und cos gemeinsam benutzt, sondern es so gemacht.
(Ich hatte vor auf diese Art auch einfach Ellipsen zu erstellen...)
for(int z=0 ; z<10000 ; z++)
{
float winkel=MyRandom(360.0);
float x = sin( winkel * (Pi / 180.0) ) * MyRandom(200.0);
float y = cos( winkel * (Pi / 180.0) ) * MyRandom(200.0);
}
Zuerst kommt der Winkel. Daraus entsteht ein Punkt auf dem Kreis mit Radius 200, und das gleichverteilt. Beachte nun das Rechteck, das den Ursprung und diesen Punkt gegenüber hat. In dieses Rechteck wird ein Punkt gesetzt, und das gleichverteilt..hikaru hat geschrieben:07.01.2019 10:36:28Ich verstehe noch nicht, warum das zu der Kreuzbildung führt. Wenn man zwei verschiedene Zufallszahlen für den Radius in x- und y-Richtung wählt, liegt der Punkt zwar auf einer Ellipse statt auf einem Kreis, aber gemittelt über 10000 Punkte würde ich trotzdem eine gleichmäßige Verteilung* erwarten.
Fuer mich leider nicht...
Code: Alles auswählen
Pseudocode:
1 while ( Count < 10000 )
2 {
3 Wxy=rand()%360;
4 Wz=rand()%360;
5 Rx=rand()%400-200;
6 Ry=rand()%200-100;
7 Rz=rand()%100-50;
8 X = sin( Wz * (pi/180) ) * cos( Wxy * (pi/180) ) * Rx ;
9 Y = sin( Wz * (pi/180) ) * sin( Wxy * (pi/180) ) * Ry ;
10 Z = cos( Wz * (pi/180) ) * Rz ;
11 Count++;
12 }
Code: Alles auswählen
9 Y = ( ( sin( Wz * (pi/180) ) * sin( Wxy * (pi/180) ) ) / 2 ) * Rx ;
Code: Alles auswählen
9 Y = ( ( sin( Wz * (pi/180) ) * sin( Wxy * (pi/180) ) ) / ( Rx / Ry ) ) * Rx ;
Du musst präzise formulieren, was du haben willst. Dann hast du auch schon die Lösung gefunden.heinz hat geschrieben:07.01.2019 23:18:54Lassen wir mal die haeufung der Punkte in der Mitte ausser acht.
Wie muesste denn die Formel aussehen um fuer x, y und z verschiedene Radien angeben zu koennen, ohne das ein Kreuz entsteht?
Code: Alles auswählen
int radX=2000;
int radY=1000;
int radZ=500;
double rad_faktorY=(double)radX/(double)radY;
double rad_faktorZ=(double)radX/(double)radZ;
while(count < 10000)
{
double winkel=(double)(rand()%36000000)/100000.0;
double winkelZ=(double)(rand()%36000000)/100000.0;
double radius=rand()%radX;
double XX=sin(winkelZ*(Pi/180.0)) * cos(winkel*(Pi/180.0)) * radius;
double YY=( ( sin(winkelZ * (Pi/180.0)) * sin(winkel * (Pi/180.0)) ) / rad_faktorY ) * radius;
double ZZ=( cos(winkelZ * (Pi/180.0)) / rad_faktorZ ) * radius;
}
leider nicht.reox hat geschrieben:11.01.2019 08:56:27Kann es sein, dass das nur ein Artefakt ist, weil du die Achsen unterschiedlich skalierst?
Code: Alles auswählen
int radX=2000;
int radY=2000;
int radZ=2000;
Code: Alles auswählen
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
N = 10000
r = 1
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
theta = np.pi * np.random.rand(N)
phi = 2 * np.pi * np.random.rand(N)
radius = r * np.random.rand(N)
x = np.sin(theta) * np.cos(phi) * radius;
y = np.sin(theta) * np.sin(phi) * radius;
z = np.cos(theta) * radius
ax.scatter(x, y, z)
plt.show()
Also wenn ich die Z-Koordinaten aus der erzeugten POV-Datei mal nach der Haeufigkeit aufliste, erkennt man sehr gutreox hat geschrieben:11.01.2019 16:49:02Eine Frage wäre jetzt: Wie kann man zB ein Histogramm erzeugen aus den Koordinaten um zu schauen ob die Verteilung passt?
Code: Alles auswählen
grep '^translate' ./zufall.pov |cut -f3 -d','|tr -d '>'|cut -f1 -d'.'|sort -g|uniq -c >./z-h.txt