Oui en effet, lorsqu'on cherche à simuler des variables aléatoires indépendantes, il y a beaucoup d'algorithmes disponibles pour toutes les lois imaginables. Ce polycopié traite assez bien le cas "indépendant":
[url=http://]http://samos.univ-paris1.fr/archives/ftp/cours/millet/montecarlo.pdf[/url]
Mais pour simuler la corrélation, il semble y avoir moins de recherche sur le sujet.
En ce qui concerne, mon problème, j'ai trouvé une première solution mais qui ne fonctionne que pour des corrélations faibles (et positives, mais c'est le cas qui m'intéresse). Je la poste ici: peut-être les commentaires aideront à l'améliorer et pourquoi pas trouver une solution pour une corrélation forte.
On définit X,Y,Z trois lois de Poisson indépendantes telles que:
et on pose:
et
ce qui permet bien d'avoir 2 lois de Poisson (puisque la somme de 2 lois de Poisson indépendantes est encore une loi de Poisson) avec un coefficient de corrélation de
. Seulement il faut respecter