La modélisation de la diffusion par des matériaux hétérogènes rugueux est un sujet complexe tant au niveau théorique que numérique. Dans notre travail, les dimensions des différents éléments constitutifs du milieu sont dans le domaine résonant et leurs propriétés sont caractéristiques du domaine optique ou infrarouge. Nous chercherons à calculer numériquement la Fonction de Distribution de la Réflectance Bidirectionnelle pour une géométrie 2D. Nous développons une approche de Monte Carlo basée sur la combinaison d'une méthode des différences finies dans le domaine temporel (FDTD) pour le calcul électromagnétique et de méthodes de génération de milieux hétérogènes rugueux (Metropolis). <br /><br />Nous appliquons cette approche pour étudier la diffusion par des surfaces rugueuses de fonctions d'auto-corrélation gaussienne ou exponentielle. Les résultats des calculs sont comparés avec ceux donnés par la Méthode des Moments. Ensuite, nous proposons une méthode originale de détermination d'un indice effectif relié à la propagation de l'énergie cohérente. Nous l'utilisons pour déterminer les domaines de validité de modèles approchés : Maxwell-Garnett, Bruggeman, Foldy-Twersky et Keller. Nous utilisons ces résultats pour étudier la diffusion en champ lointain par un ensemble bidisperse de particules, et notamment la possibilité de remplacer l'influence des petites particules par le milieu homogénéisé décrit précédemment. Enfin, nous étudions le problème du couplage surface/volume. Nous vérifions une hypothèse de découplage basée sur l'utilisation de l'indice effectif défini précédemment pour la diffusion surfacique notamment dans le cas où la rugosité est directement corrélée à la position des particules.