simulation acoustique par la methode des sources images

fonction « croisement.m » qui calcule les coordonnées du point d’intersection de deuxsegments en fonction des coordonnées de leurs extrémités, renvoie un 1 si le pointd’intersection appartient bien au segment et 0 sinon. Ainsi, il est possible de dire si unrayon acoustique se réfléchit bien à l’intérieur d’une paroi pour garantir la visibilitéd’une source.- Si une paroi potentiellement masquante se trouve en travers du chemin d’un rayonalors, le programme détecte cette intersection et la source devient invisible.Ici, la fonction « croisement.m » est au centre de l’algorithme et mériterait un travaild’optimisation pour accélérer les calculs. Elle utilise la sous-fonction « intersection.m » quirenvoie les coordonnées du point d’intersection de deux droites à partir de leurs coefficientshomogènes. Le procédé utilisé est un solveur d’équation linéaire (« linsolve.m ») proposé parMatLab qui résoud le système d’équation définit par ces droites. C’est précisément cetalgorithme qui pourrait faire l’objet d’une attention particulière lors d’une étude ultérieure.Tout comme pour la gestion de l’absorption par les parois de la salle, la fonction« sourcevu.m » doit remonter l’historique des sources pour retracer les rayons point par point(voir II.1.3 pour la méthode de construction des rayons). Elle utilise encore le fait que laposition de la source mère est indiquée à chaque ligne du tableau « mater ». Il est ainsipossible de retrouver la suite de réflexions ayant engendrée une source.Ici, nous avons pu exposer les problématiques ainsi que les outils mathématiquespermettant de traiter des questions de critères de sources dans le cas général. Afin de bienexpliquer ces notions, il est essentiel de comprendre le passage de la constellation au tracé desrayons acoustiques qui, en définitive, ne sont que deux manières distinctes de représenter unmême phénomène.73