On 19 Set, 11:04, davide.dei..._at_gmail.com wrote:
> On 18 Set, 17:10, Henry <enrico.ce..._at_gmail.com> wrote:
>
> > Prima di rispondere qualche domanda,
>
> > Che generatore di numeri pseudo-casuali stai usando ?
>
> > Sarebbe interessante vedere la porzione di codice in cui applici
> > l'algoritmo di Metropolis.
>
> > Per determinare la bont� di un generatore esistono diversi test,
> > quello proposto nell' esercitazione consiste nel vedere con quale
> > frequenza (altezza dell'istogramma) un certo numero viene generato.
> > Una caratteristica fondamentale di un generatore � che i numeri che
> > vengono generati siano uniformemente distribuiti (tipicamente nell'
> > intervallo [0,1]) e che abbia un periodo molto grande.
>
> > Henry
>
> il generatore di numeri e' il ran1 del fortran 77 che genera numeri
> compresi tra 0,1 a distr.uniforme.
> non capisco fondamentalmente una cosa: io la sequenza degli xi la
> genero con il metropolis; questa sequenza devo generarla ulteriormente
> per esempio m volte??? perche per un certo N hai una lunghezza della
> sequenza...se poi devo vedere quante volte un numero e' stato generato
> dovrei ripetere il tutto altre m volte....
Scusa ma non capisco..... :(
Un conto sono i numeri generati "casualmente" altra cosa � il criterio
con il quale accetti o meno un risultato con l'algoritmo di Metropolis
(come detto prima, sarebbe interessante vedere il codice o anche lo
pseudo codice sul quale stai lavorando).
Nota sul generatore di numeri pseudo casuali:
Fai molta attenzione al generatore che usi, nel caso in esame
probabilmente la casualit� dei numeri generati *forse* non � cos�
critica, ma se in futuro dovrai affrontare problemi pi� articolati
sappi che dovrai controllare accuratamente il generatore. Come ti
dicevo nel mio precedente post esistono alcuni test per controllare
l'efficacia dei generatori, tra questi il test di Kolmogorov-Smirnov e
un test "spettrale" che tu stesso puoi provare, potrai osservare che
alcuni dei pi� "famosi" generatori di numeri casuali, non sono poi
cos� casuali. Tra questi il generatore del fortran, quello di Matlab,
Mathematica e Derive sono da evitare per simulazioni "serie". Ad oggi
che io sappia uno dei migliori generatori (presenti anche in sorgente
per Fortran) � il Mersenne Twister di cui ti lascio il link...... Buon
divertimento :)
http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/VERSIONS/FORTRAN/fortran.html
Ciao
Henry
Received on Wed Sep 19 2007 - 22:15:28 CEST