Распространение радиоволн ВЧ/Метод зеркальных отображений

Алгоритм

Основное течение

1. $\forall ~ \rho_k \in$ G.Множество контрольных точек()
   1. $\forall ~ s_j \in$ G.Множество первичных источников()
      1. $s'\leftarrow s_j$
      2. $stack_{t}' \leftarrow t0_{mn}$
      3. $stack_{\rho}' \leftarrow \rho_k$
      4. $I \leftarrow 0$
      5. Переход 1.1.8.1.4
      6. $stack_{t}' \leftarrow stack_{t}$
      7. $stack_{\rho}' \leftarrow stack_{\rho}$
      8. $\forall ~ f_m\in$ G.Множество отражающих объектов()
         1. $\forall ~ t_{mn} \in f_m$.Множество отражающих поверхностей()
            1. $\rho' \leftarrow$ $\rho'$.Построить зеркальное отображение($t_{mn}$.Плоскость грани())
            2. $stack_{t}' \leftarrow t_{mn}$
            3. $stack_{\rho}' \leftarrow \rho'$
            4. $I' \leftarrow I$
            5. Position $\leftarrow$ $s_j$.Антенна().Позиция()
            6. $\rho' \leftarrow stack_{\rho}'$
            7. $t_{mn} \leftarrow stack_{t}'$
            8. Ray $\leftarrow$ Луч().Создать(Position, Вектор(Position,$\rho'$.Позиция(), 1))
            9. Если $I' > 0$
               1. $t_{mn} \leftarrow stack_{t}'$
               2. $[cross, P] \leftarrow$ Ray.Пересечение($t_{mn}$.Плоскость грани())
               3. Если $cross$
                  1. Если $t_{mn}$.Принадлежность($P$) $\wedge$ Расстояние(Position,$\rho'$.Позиция()) $>$ Расстояние(Position,$P$)
                     1. Переход 1.1.8.1.10
                  2. Иначе
                     1. Переход 1.1.8.1.13
               4. Иначе
                  1. Переход 1.1.8.1.13
            10. $\forall ~ f_m'\in$ G.Множество отражающих объектов()
                1. $\forall ~ t_{mn}' \in f_m'$.Множество отражающих поверхностей()
                   1. $[cross', P'] \leftarrow$ Ray.Пересечение($t_{mn}'$.Плоскость грани())
                   2. Если $cross'$
                      1. Если $t_{mn}'$.Принадлежность($P'$) $\wedge$ Расстояние(Position,$P$) $>$ Расстояние(Position,$P'$)
                         1. Переход 1.1.8.1.13
            11. Если $I' > 0$
                1. Angle $\leftarrow$ Ray.Угол пересечения($t_{mn}$.Плоскость грани())
                2. $s_j \leftarrow$ Вторичный источник при рейтрейсинге.Создать($s_j$.Напряженность$(\omega_n,~\theta,~\varphi$, Расстояние(Position,$P$), G.Среда распространения()$),$ Angle, $P,~t_{mn}$)
                3. $\theta, \varphi \leftarrow 0$
                4. $I' \leftarrow I'-1$
                5. Переход 1.1.8.5
            12. $\rho_k$.Зарегистрировать($s_j$.Напряженность$(\omega_n,~\theta,~\varphi$, Расстояние(Position, $\rho_k$.Позиция()), G.Среда распространения()$)$)
            13. Если $E > E_{end}$
                1. $I \leftarrow I+1$
                2. $stack_t \leftarrow t_{m,n}$
                3. $stack_{\rho} \leftarrow \rho'$
                4. Переход 1.1.6
            14. $I \leftarrow I-1$
            15. Если $I < 0$
                1. Переход 1.1.9
            16. $t_{m,n} \leftarrow stack_t$
            17. $\rho' \leftarrow stack_{\rho}$
      9. $s_j \leftarrow s'$