Текущая версия на 16:12, 6 декабря 2016

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

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

Алгоритм

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

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

Архитектурная акустика/Рей-трейсинг: различия между версиями

Текущая версия на 16:12, 6 декабря 2016

Алгоритм

Свойства алгоритма

Навигация

Архитектурная акустика/Рей-трейсинг: различия между версиями

Текущая версия на 16:12, 6 декабря 2016

Алгоритм

Свойства алгоритма

Навигация

Поиск