zufallszahlen (Autodesk/Lisp)

Ja ... nach welcher statistischen Verteilung (gleichverteilt / Normalverteilt )denn ?

Ist zwar standardmäßig nicht in Autolisp integriert, kann man sich aber proggen
Z.B. für Gleichverteilte

;;; - ------------------------------------------------------------------------------ - ;
;;; - Funktionsname : DT:RANDOM-U01 - ;
;;; - Beschreibung : Erzeugung einer Liste mit gleichverteilten Zufallszahlen - ;
;;; - zwischen 0 und 1 der angegebenen Listlänge U(0,1) - ;
;;; - Parameter : ANZAHL [INT] - Anzahl der zu ermittelnen Zufallszahlen - ;
;;; - X0 [INT] - Startwert für Generator - ;
;;; - bei nil wird ein Startwert generiert - ;
;;; - Rückgabe : RANDOM [LIST(REAL)] - Liste mit Pseudozufallszahlen - ;
;;; - Beispiel : (DT:RANDOM 100 nil) - ;
;;; - ------------------------------------------------------------------------------ - ;
;;; - Info : Ein solcher Zufallsgenerator heißt linearer Kongruenzgenerator,- ;
;;; - wobei man im Fall C=0 auch von einem multiplikativen Generator,- ;
;;; - im Falle C/=0 von einem gemischten Kongruenzgenerator spricht - ;
;;; - * Es gilt : die Zufallszahlen treten periodisch auf - ;
;;; - * Die Periode hängt stark von den gewählten Parametern ab - ;
;;; - * Ein linearer kongruenter Generator kann max. (m-1) - ;
;;; - verschiedene Zufallszahlen liefern (vollständige Periode) - ;
;;; - Vollständige Perioden erhält man wenn folgende Theoreme erfüllt sind - ;
;;; - 1) c und m sind relative Primzahlen, haben also nur 1 als gemeinsamen Teiler - ;
;;; - 2) Jede Primzahl, die m dividiert, dividiert auch a-1 - ;
;;; - 3) a-1 ist durch 4 teilbar, wenn m es auch ist - ;
;;; - Da m die Periodenlänge bestimmt, ist für m eine große Primzahl zu wählen - ;
;;; - Für a sollte gelten 0.01*m<a<0.99*m und die Ziffern von a sollten kein - ;
;;; - einfaches Muster aufweisen. nach Fishmann/Greenb. a = sqrt(m-6c/m*(1-c/m)) - ;
;;; - ------------------------------------------------------------------------------ - ;
(defun DT:RANDOM-U01(ANZAHL X0 / M A C X RANDOM)
(DT:RANDOMINIT X0)
(if(=(type ANZAHL)'INT)
(progn
(setq M 2147483647) ;2^31-1 ;_Konstante für Modul
(setq A 65539) ;_Konstante für Multiplikator
(setq C 2957) ;_Konstante für Inkrement
(setq X *RANDOMSEED*) ;_Startwert
(repeat ANZAHL
(setq X (rem(+(* (float A) X)C)M))
(setq RANDOM(cons (/ X (float M)) RANDOM))
)
(setq *RANDOMSEED* X) ;_setzen de SEED-Variable, um bei weitern Aufrufen
) ;_den Generator nicht neu initialisieren zu müssen
)
RANDOM
)

------------------
- Thomas -
"Bei 99% aller Probleme ist die umfassende Beschreibung des Problems bereits mehr als die Hälfte der Lösung desselben."

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

RTOPO Hydrology CAD for topography CAD APP für GIS, Vermessung und Kartografie