Архитектурная акустика/Рей-трейсинг

Архитектурная акустика\Рей-трейсинг

Модель распространения звука в помещении, в которой каждая поверхность зеркально отражает падающий на нее луч. Таким образом, каждая отражающая помещения является вторичным источником в заданном направлении, которое определяется через угол отражения, равный углу падения.

Алгоритм

1. Для $\forall ~ s_i \in$ G.Множество первичных источников()
   1. Для $\forall ~ \rho_i \in$ G.Множество плоскостей вывода результатов()
      1. Для $\forall ~ p \in \rho_i$
         1. $D \leftarrow p - s_i.$Позиция()
         2. Для $\forall ~ r \in$ G.Множество отражающих объектов()
            1. Для $\forall ~ f \in r.$Множество поверхностей()
               1. Если $\exists t=\vec{D}\cap f$ И расстояние $(s_i,t)< |\vec{D}|$
                  1. Перейти на шаг 1.1.1
         3. Для $\forall ~ k \in$[$0,s_i.$Ширина полосы частот())
            1. $\omega \leftarrow s_i.$Частота($k$)
            2. $\rho_i.$ Зарегистрировать звук ($s_i.$ Интенсивность (Азимут $\vec{D},$ Зенит$\vec{D},\omega),0,\omega)$
   2. Для $\forall ~ r \in$ G.Множество отражающих объектов()
      1. Для $\forall ~ f \in r.$ G.Множество поверхностей()
         1. Для $\forall ~ t \in f.$ G.Множество отражающих элементов(Шаг дискретизации)
            1. Для $\forall ~ k \in$[$0,s_i.$Ширина полосы частот())
               1. $\omega \leftarrow s_i.$Частота($k$)
               2. $\vec{D} =\vec{t}-s_i.$Позиция()
               3. $I\leftarrow s_i.$Интенсивность(Азимут$\vec{D},$ Зенит$\vec{D},\omega$)
               4. $\vec{D'} \leftarrow$Отразить$\vec{D'},t$
               5. $I' \leftarrow (1-f.$коэффициент отражения$(\omega))\cdot I$
               6. Для $\forall ~ \rho_i \in$ G.Множество плоскостей ()
                  1. Если $\exists p=\alpha\vec{D'}\cap \rho_i$
                     1. $\rho_i.$ Зарегистрировать${\frac {I'}{|t-p|^{2}}},p,\omega )$
               7. Для $\forall ~ r' \in$ G.Множество отражающих объектов ()
                  1. Для $\forall ~ f' \in r'$ G.Множество плоскостей ()
                     1. Если $\exists t'=\beta\vec{D'}\cap f'$
                        1. $D' \leftarrow$ Отразить $(D',f')$
                        2. $I'\leftarrow {\frac {I'(1-f'.absorptance(\omega )}{|t-t'|^{2}}}$
                        3. $t \leftarrow t'$
                        4. Переход на шаг 1.2.1.1.1.5