simulation acoustique par la methode des sources images

III.2.3 Programmation : les fonctions « absorption.m » et « reconstruction.m »Gestion de l’absorption par les parois.La colonne n°8 de la matrice des sources images (voir II.1.4) indique le coefficientd’atténuation dû à l’absorption par les parois : la fonction « absorption.m » se charge dedéterminer ce paramètre pour chaque source.Les sources construites de façon géométrique sont, par définition, engendré se par unesource « mère » via une paroi. Or, le tableau « mater » des sources images indique pourchaque source sa provenance, c'est-à-dire le numéro de la paroi et la position de la source« mère » qui l’ont engendrée. La nouvelle atténuation de la source est obtenue simplement enmultipliant le nouveau coefficient de réflexion de cette paroi par le coefficient d’atténuationdéjà affectée pour la source « mère ». En analysant chaque source et en posant commecondition initiale que la source primaire n’est pas atténuée, un calcul par récurrence permet dedéterminer les atténuations pour toutes les lignes du tableau.Les sources stochastiques sont générées de façon aléatoire (voir III.4), il n’est doncpas possible de retrouver leurs généalogies. En fait, le coefficient d’absorption de ces sourcesest défini de façon statistique en fonction de paramètres géométriques. Si r est la distance dela source à l’auditeur, l le libre parcours moyen dans la salle et b reflex le coefficient deréflexion moyen à considérer, alors, le coefficient d’atténuation cherché est :rlbreflexNous justifierons ce principe en détail au chapitre III.4 mais il faut noter que le coefficientb reflex à considérer est le résultat d’une analyse statistique des coefficients des sourcesgéométriques déjà calculées.La fonction « absorption.m », nous l’avons dit, est un algorithme spécifiquementdestiné à l’attribution de ces coefficients. Or certaines opérations sont redondantes avec lecalcul de la construction des sources images. Il nous semblerait donc plus simple de combinerces taches en une seule et même fonction. Cependant, le découpage en sous-bandesnécessiterait de recalculer inutilement la constellation de sources images pour chacune desbandes. En séparant la construction spatiale des sources de l’application des coefficientsd’atténuation, seul ce qui change est recalculé.92