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

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

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

Алгоритм

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