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

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

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

Алгоритм

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$ и если не в_тени$(t,p,G)$
                     1. $\rho_i.$ Зарегистрировать${\frac {I'}{|t-p|^{2}}},p,\omega ))$
               7. Для Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): \forall ~ r' \in G.Множество отражающих объектов ()
                  1. Для Невозможно разобрать выражение (Ошибка преобразования. Сервер («https://wikimedia.org/api/rest_») сообщил: «Cannot get mml. Server problem.»): \forall ~ f' \in r' G.Множество плоскостей ()
                     1. Если Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): \exists t'=\beta\vec{D'}\cap f'
                     2. Если не В_тени (t,t',G)
                        1. Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): D' \leftarrow Отразить Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): (D',f')
                        2. Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle I' \leftarrow \frac{I'(1-f'.absorptance(\omega)}{|t-t'|^2}}
                        3. Невозможно разобрать выражение (Ошибка преобразования. Сервер («https://wikimedia.org/api/rest_») сообщил: «Cannot get mml. Server problem.»): t \leftarrow t'
                        4. Переход на шаг 1.2.1.1.1.5