Moltiplico max per un random tra 0 ed 1, non tra 0 e rand_max. Il metodo accept/reject dice che prima estraggo un numero random x, e poi estraggo un secondo numero random r2 tra 0 e max (max > f(x) per ogni x). Se r2 < f(x) accetto x, se r2 > f(x) scarto x e riprovo. Quando finalmente un x viene accettato, aggiungo questo nuovo x al set di numeri random che stavo cercando. AB
praticamente: -estraggo un r1 tra 0 e 1 -estraggo un r2 tra 0 e 2 -se r2 < exp(-4r1*r1) allora accetto r1 altrimenti lo scarto. Ma se estraggo tipo lo zero(r1) e poi 2(r2),lo zero viene eliminato anche se in realtà è quello con frequenza maggiore!quindi non ho capito un ca**o...
No, hai capito benissimo, la tecnica è quella che tu descrivi. Lo zero può venir eliminato, ma può anche farcela a passare (dipende da r2), e questo succede più facilmente a lui che ad un r1 vicino ad 1. Prova. AB
Ma se per caso lo scarto, poi nella mia estrazione mi manca lo zero. Ma in realtà questo non avviene perchè faccio tante estrazioni e quindi prima o poi anche lo zero verrà riconsiderato,no? E sicuramente verrà riconsiderato molte più volte rispetto agli altri numeri dato che il valore della funzione per lui è molto alto. Evviva!forse ho capito(spero),grazie!
Esatto. AB PS: se vai a vedere su wikipedia (inglese) le voci: montecarlo method montecarlo option model montecarlo integration puoi farti un'idea del 10 percento delle applicazioni di questa tecnica.
non ho capito il procedimento:perchè pone un val max=2 e poi lo moltilpca per un valore casuale tra 0 e rand max per ottenere r2?
RispondiEliminaMoltiplico max per un random tra 0 ed 1, non tra 0 e rand_max.
RispondiEliminaIl metodo accept/reject dice che prima estraggo un numero random x, e poi estraggo un secondo numero random r2 tra 0 e max (max > f(x) per ogni x). Se r2 < f(x) accetto x, se r2 > f(x) scarto x e riprovo.
Quando finalmente un x viene accettato, aggiungo questo nuovo x al set di numeri random che stavo cercando. AB
L'esercizio consiste nell'estrarre numeri secondo la frequenza che voglio io?
RispondiEliminaSi. Io ho usato una funzione gaussiana, ma è chiaro che se uno riesce a farlo con una funzione riesce a farlo con tutte. AB
RispondiEliminapraticamente:
RispondiElimina-estraggo un r1 tra 0 e 1
-estraggo un r2 tra 0 e 2
-se r2 < exp(-4r1*r1) allora accetto r1 altrimenti lo scarto.
Ma se estraggo tipo lo zero(r1) e poi 2(r2),lo zero viene eliminato anche se in realtà è quello con frequenza maggiore!quindi non ho capito un ca**o...
No, hai capito benissimo, la tecnica è quella che tu descrivi. Lo zero può venir eliminato, ma può anche farcela a passare (dipende da r2), e questo succede più facilmente a lui che ad un r1 vicino ad 1. Prova. AB
RispondiEliminaMa se per caso lo scarto, poi nella mia estrazione mi manca lo zero. Ma in realtà questo non avviene perchè faccio tante estrazioni e quindi prima o poi anche lo zero verrà riconsiderato,no?
RispondiEliminaE sicuramente verrà riconsiderato molte più volte rispetto agli altri numeri dato che il valore della funzione per lui è molto alto.
Evviva!forse ho capito(spero),grazie!
Esatto. AB PS: se vai a vedere su wikipedia (inglese) le voci:
RispondiEliminamontecarlo method
montecarlo option model
montecarlo integration
puoi farti un'idea del 10 percento delle applicazioni di questa tecnica.