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

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

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

Алгоритм

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. Если $I' > G$.Минимальный уровень звука()
                  1. Для $\forall ~ \rho_i \in$ G.Множество плоскостей ()
                     1. Если $\exists p=\alpha\vec{D'}\cap \rho_i, 0 < \alpha < \infty$ и если не в_тени$(t,p,G)$
                        1. $\rho_i.$ Зарегистрировать$( \frac{I'}{|t-p|^2},p,\omega)$
                  2. Для $\forall ~ r' \in$ G.Множество отражающих объектов ()
                     1. Для $\forall ~ f' \in r'$ .Множество поверхностей()
                        1. Если $\exists t'=\beta\vec{D'}\cap f', 0 < \beta < \infty$
                           1. Если не В_тени (t,t',G)
                              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.6