BSD

Bonjour,

J'ai un projet à rendre assez urgemment, en physique informatique, et je ne sais tout simplement pas comment utiliser l'algorythme donné par le prof:


L'énoncé: Utiliser la méthode DMC (Diffusion Monte Carlo) pour calculer la fonction d'onde de l'état fondamental de l'atome d'hydrogène.


L'algorythme: 1. Eparpiller N marcheurs aléatoires. On attribue à chacun un poids W.

2. A chaque intervalle de temps dt,
a) pour chaque marcheur, R'=R+€*(racine de 2*D*dt)=R+€*(racine de dt) lorsque la densité de probabilité est g(€)=exp(-€^2/2dt)/((2*pi*dt)^(d/2));
b) son poids devient W'=W*exp(dV(R)*dt), avec dV=combien le poids a augmenté ou diminué.

3. Si W>1, si nc=int(W+z)>=2, créer nc copies du marcheur au point R' avec un poids W'=W-1;
si W<1, effacer le marcheur avec une probabilité de 1-W

[ici, D=1/2; d est le nombre de dimensions; "int" est la partie entière; et je ne sais pas ce que représente z]

J'aimerais faire tout ça dans Matlab, mais si la personne qui peut m'aider préfère un autre programme, ça peut se faire aussi.

Merci d'avance à qui passera du temps à m'expliquer tout ça