wie erzeugt integer-Zufallszahlen in fortran 90 im Bereich [0,5]?

Ich bin Art von neuen in die fortran-proramming.
Kann jemand bitte helfen Sie mir mit der Lösung.
ich habe ein problem in der Erzeugung von ganzzahligen Zufallszahlen
im Bereich [0,5] in fortran random number using
random_seed und rand

  • rand ist eine GNU-Erweiterung. Der Fortran 95-Kombination ist random_seed und random_number.
  • Ist eine Folge von genau 5, für Sie akzeptabel? Oder meinten Sie bitten, für den Bereich [0,5) ? Die Vorschläge, die unten nicht geben Ihnen genau den Endpunkt Wert (außer vielleicht durch Rundung.)
  • Ich bin mir nicht ganz sicher, was du meinst... Könntest du das näher?
  • Man verwendet den Begriff "[0,5]". Dies ist die "Intervall-notation", das heißt, der Wert kann überall von exakt 0, um genau 5. Auf die überprüfung, dass vielleicht, was Sie wollen. francescalus' Antwort weiter unten ausführlich erklärt.
InformationsquelleAutor user3531410 | 2014-04-14
  1. 4

    Was: 

    program rand_test
  use,intrinsic :: ISO_Fortran_env
  real(REAL32)  :: r(5)
  integer       :: i(5)

  ! call init_random_seed() would go here

  call random_number(r)

  ! Uniform distribution requires floor: Thanks to @francescalus 
  i = floor( r*6._REAL32 )

  print *, i
end program
    • Verwenden real :: ... und 6. dann. Aber bitte, tun Sie sich einen gefallen und wechseln Sie zu expliziten Definitionen!
    • i = floor( r*3. ) - 1
    • das ist ok, aber wird es geben, gleichmäßige Verteilung .....?\
    • wird es geben, gleichmäßige Verteilung für -1,0,1?
    • Ich sehe nicht, warum nicht - du bist nur eine Verschiebung der Grenzen.
    • YA sicher. ein weiteres problem ist es, die mit floating-point-Fehler M6101:MATH. WARUM IST ES PASSIERT?
    • Was ich fand, dass, wenn ich bin deviding einige integer-Zahl 580 mit der Nummer 1.97 e-7, dann ist die division geben einen negativen integer-Zahl. warum ist das so?
    • derzeit werden eine Menge für die Unterstützung, die mir mit der random Nummer..Ya ich bin nur eine Verschiebung der Grenzen.Dank

    InformationsquelleAutor Alexander Vogt
  2. 4

    Zur Unterstützung der Antwort von Alexander Vogt, ich werde verallgemeinern.

    Die intrinsische random_number(u) gibt eine reelle Zahl u (oder ein array von solchen) aus der Gleichverteilung über dem Intervall [0,1). [Das heißt, es enthält 0 aber nicht 1.]

    Haben eine diskrete Gleichverteilung auf den ganzen zahlen {n, n+1, ..., m-1, m} schnitzen die kontinuierliche Verbreitung bis in die m+1-n gleich große Stücke, die Abbildung jeder chunk eine ganze Zahl. Eine Möglichkeit könnte sein:

    call random_number(u)
j = n + FLOOR((m+1-n)*u)  ! We want to choose one from m-n+1 integers

    Wie Sie sehen können, für die erste Frage {0, 1, 2, 3, 4, 5} diese reduziert sich auf

    call random_number(u)
j = FLOOR(6*u)            ! n=0 and m=5

    und für den anderen Fall in Ihrem Kommentar { -1, 0, 1}

    call random_number(u)
j = -1 + FLOOR(3*u)       ! n=-1 and m=1

    Natürlich auch andere Transformationen erforderlich, Sätze von nicht-zusammenhängenden ganzen zahlen, und man sollte darauf achten, zu numerischen Problemen.

    • um zu vermeiden, neigen aufgrund der begrenzten Anzahl von bits in real Mantisse (53 für IEEE 754 double precision, 24 für single-precision), Sie möglicherweise benötigen, rufen Sie random_number mehrere Male, z.B., siehe wie _randbelow(n) erfolgt über random() in Python
    InformationsquelleAutor francescalus
Schreibe einen Kommentar