Calcul de distances moyennes

Version imprimable

Voir 40 message(s) de cette discussion en une page

08/09/2010, 00h17
pseudocode

Citation:

Envoyé par Champialex

Finalement j'ai codé le modèle physique (villes placées initialement au hasard, de même charge, donc se repoussent, et en plus on considère une succession d'états d'équilibre : la vitesse de déplacement est nulle au début de chaque étape), ca marche pas mal, le problème c'est qu'il y a une accumulation des villes sur les bords du carré...

Il faut gérer les bords comme des particules fixes chargées => les particules sont repoussées par les bords.

1 pièce(s) jointe(s)

J'ai codé le truc, mais bon j'ai du faire une erreur dans la force qu'exerce le bord, puisque la répartition des villes n'est pas symétriques...

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
 
if(nouvelX>1/2) //nouvelX : nouvelle position de selon l'axe des x, x est l'ancienne position
{
	if(x==1)
	{
		nouvelX = nouvelX-r.nextDouble()*100000*pasTemporel*pasTemporel; //Force exercée si l'on était déja sur le bord
	}
	else
	{
		nouvelX = nouvelX - pasTemporel*pasTemporel*forceBord/Math.pow((1-x),2);
	}
}
else if (nouvelX<1/2) 
{
	if(x==0)
	{
		nouvelX = nouvelX+r.nextDouble()*100000*pasTemporel*pasTemporel;
	}
	else
	{
		nouvelX =nouvelX+pasTemporel*pasTemporel*forceBord/Math.pow(x,2);
	}
}
if(nouvelY>1/2)
{
	if(y==1)
	{
		nouvelY = nouvelY-r.nextDouble()*100000*pasTemporel*pasTemporel;
	}
	else
	{
		nouvelY = nouvelY-pasTemporel*pasTemporel*forceBord/Math.pow((1-y),2);
	}
}
else if (nouvelY<1/2) 
{
	if(y==0)
	{
		nouvelY = nouvelY+r.nextDouble()*100000*pasTemporel*pasTemporel;
	}
	else
	{
		nouvelY = nouvelY+ pasTemporel*pasTemporel*forceBord/Math.pow(y,2);
	}
}

08/09/2010, 19h50
Champialex
Le problème à été résolut grâce à l'aide d'un bienfaiteur anonyme! Il faut en fait remplacer les tests comme
Code:

1 2 if(x==1)
par :
Code:

1 2 if(1-x<tol)
où tol est un paramètre choisit de façon adaptée au problème.
12/09/2010, 14h55
Fredchkek

Citation:

Je ne vois pas trop l'intérêt de faire la boucle MC : pour qu'elle soit utile, il faudrait l'utiliser quelque part, par exemple en faisant la moyenne des distances moyennes observées...

Cette boucle est pourtant capitale ! Elle permet de faire une moyenne statistique sur la distance entre les villes. Donc elle bien utilisée !

Pour le reste des commentaires, j'ai l'impression que tu n'as pas très bien compris l'algorithme (?). Le but est bien de calculer la distance moyenne entre villes d'un voyageur qui les visites. Son choix est donc important et il a intérêt à choisir la ville la plus proche à chaque fois.

Pour le modèle physique, c'est intéressant mais cela ne permet "que" de distribuer uniformément les villes dans le pays ( au prix d'une résolution CPU-consuming j'ai l'impression).
12/09/2010, 20h13
Champialex

Citation:

Envoyé par Fredchkek

Cette boucle est pourtant capitale ! Elle permet de faire une moyenne statistique sur la distance entre les villes. Donc elle bien utilisée !

Pour le reste des commentaires, j'ai l'impression que tu n'as pas très bien compris l'algorithme (?).

Effectivement, j'avais lu trop vite... J'étais un peu trop dans mon truc :aie:

Citation:

Envoyé par Fredchkek

Pour le modèle physique, c'est intéressant mais cela ne permet "que" de distribuer uniformément les villes dans le pays ( au prix d'une résolution CPU-consuming j'ai l'impression).

Mon algo est en O(n²) , même si à mon avis le facteur devant le n² est beaucoup plus grand pour le mien :mrgreen:
Sinon, comme pour le modèle on a besoin des distances entre villes, les fournir ne demande pas trop d'effort...

Voir 40 message(s) de cette discussion en une page