Распространение радиоволн ВЧ/Реализация/ControlPointSet::EncompassingAperture

encompassing_aperture_t EncompassingAperture(const spacial_vector<double>& ptAsObservedFrom) const; //1
encompassing_aperture_t EncompassingAperture() const; //2

Для заданной точки строит сектор encompassing_aperture_t сферы с единичным радиусом, описанной вокруг этой точки, так, что все точки текущего множества контрольных точек при проекции на сферу будут лежать внутри сектора.

1. Построение ведется для заданной точки.
2. То же, но построение ведется для начала мировой системы координат.

Параметры

[in]

ptAsObservedFrom

Точка, вокруг которой строится сфера.

Возвращаемое значение:

Результирующий объект-сектор.

См. также

Source::OptimizeRadiationVectors	Строит сектор сферы encompassing_aperture_t, область которого ограничивает множество, реализуемое системой классов radiation_set_t, элементов поля, излучаемого источником.
ReflectingObject::EncompassingAperture	Для заданной точки строит сектор encompassing_aperture_t сферы с единичным радиусом, описанной вокруг этой точки, так, что отражающий объект ReflectingObject при проекции на сферу будет лежать внутри сектора.

Пусть $V={\begin{pmatrix}{\vec {v}}_{0},\cdots ,{\vec {v}}_{n-1}\end{pmatrix}}$ - непустое входное множество геометрических точек, и ${\vec {P}}$ - некоторая заданная точка, вокруг которой требуется построить сектор $\Omega =\left\langle {\vec {c}},\alpha \right\rangle$ сферы с единичным радиусом, где ${\vec {c}}$ - вектор, ассоциированный с направленным отрезком от ${\vec {P}}$ до центральной точки сектора на поверхности сферы, и $\alpha$ - ангулярный радиус сектора, как это реализуется классом encompassing_aperture_t. Проекции всех точек из $V$ на сферу должны принадлежать $\Omega$.

Пусть также все точки из $V$ лежат по одну сторону от $P$, т.е. $P$ не находится геометрически в облаке, образованном $V$. Если это не так, то множество $V$ должно быть разбито на подмножества точек, лежащих по одну сторону от $P$, а затем сектора, образованные применением описываемого алгоритма над $P$ и подмножествами, должны быть объединены с помощью функции unify.

Фактически, это означает, что все углы, образованные всеми возможными парами отрезков от $P$ до $\forall {\vec {v}}_{i}\in V$ должны быть меньше $\pi$. Это имеет место, когда $V$ образовано, например, вершинами одной отражающей поверхности ReflectingObject или контрольными точками одной плоскости вывода результатов PlainControlPointSet. См. также Source::OptimizeRadiationVectors.

Поиск центральной точки сферы и ангулярного радиуса при конструировании объекта encompassing_aperture_t на основе позиции точки ${\vec {P}}$, вокруг которой строится сфера, и двух точек, которые ограничивают площадь сектора. Вектора ${\vec {v}}_{1}$ и ${\vec {v}}_{2}$ на рисунке соответствуют векторам ${\vec {v'}}_{i}$ и ${\vec {v'}}_{j}$ в статье.

Если $\left|V\right|\equiv 1$, то $\Omega =\left\langle {\frac {{\vec {v}}_{0}-{\vec {P}}}{\left|{\vec {v}}_{0}-{\vec {P}}\right|}},0\right\rangle$. Далее по тексту $\left|V\right|\geq 2$.

Создадим множество $V'=\bigcup \limits _{i=0}^{n-1}\left({\vec {v}}_{i}-{\vec {P}}\right)$ и выберем пару ${\vec {v'}}_{i},{\vec {v'}}_{j}\in V'$ тех точек, отрезки из ${\vec {P}}$ до которых образуют между собой наибольший угол $2\alpha$.

Если угол равен нулю, то $\Omega =\left\langle {\frac {{\vec {v'}}_{i}}{\left|{\vec {v'}}_{i}\right|}},0\right\rangle$.

Поэтому далее по тексту рассматривается оставшийся случай, когда $\left|V\right|\geq 2$, $0<2\alpha <\pi$ и $0<\alpha <{\frac {\pi }{2}}$.

Тогда задача нахождения $\Omega$ будет сводиться к нахождению биссектрисы между Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \vec{v'}_i} и ${\vec {v'}}_{j}$, лежащий на ней вектор Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \vec{c}} единичной длины и угол $\alpha ={\frac {1}{2}}\arccos {\frac {{\vec {v'}}_{i}\cdot {\vec {v'}}_{j}}{\left|{\vec {v'}}_{i}\right|\cdot \left|{\vec {v'}}_{j}\right|}}$.

Для нахождения вектора Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \vec{c}} соединим прямой точки ${\vec {v'}}_{i}$ и Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \vec{v'}_j} и обозначим как $\tau {\vec {c}}$ вектор, с которым ассоциирован направленный отрезок, имеющий начало в точке ${\vec {P}}$, направленный вдоль биссектриссы между Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \vec{v'}_i} и ${\vec {v'}}_{j}$ до точки пересечения с прямой. Обозначим как Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \vec{v}_c} вектор, с которым ассоциирован направленный отрезок с началом в конце отрезка ${\vec {v'}}_{i}$ и концом в точке пересечения. То есть

Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \vec{c}=\vec{v}_i+\vec{v}_c} .

Согласно теореме о биссектрисе эта точка делит отрезок между ${\vec {v'}}_{i}$ и Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \vec{v'}_j} в пропорции, равной отношению длин $\left|{\vec {v'}}_{i}\right|$ и Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \left|\vec{v'}_j\right|} . Поэтому

${\vec {v}}_{c}=\left|{\vec {v}}_{c}\right|{\frac {{\vec {v'}}_{j}-{\vec {v'}}_{i}}{\left|{\vec {v'}}_{j}-{\vec {v'}}_{i}\right|}}$;

Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \frac{\left|\vec{v}_c\right|}{\left|\vec{v'}_j-\vec{v'}_i\right|-\left|\vec{v}_c\right|}=\frac{\left|\vec{v}_i\right|}{\left|\vec{v}_j\right|}\Rightarrow\vec{v}_c=\frac{\left|\vec{v}_i\right|\cdot\left|\vec{v}_j-\vec{v}_i\right|}{\left|\vec{v}_i\right| + \left|\vec{v}_j\right|}} .

Поэтому

$\tau {\vec {c}}={\frac {1}{\left|{\vec {v}}_{i}\right|\cdot \left|{\vec {v}}_{j}\right|}}\left(\left|{\vec {v}}_{j}\right|{\vec {v}}_{i}+\left|{\vec {v}}_{i}\right|{\vec {v}}_{j}\right)$.

Пусть Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \tau'=\frac{\tau}{\left|\vec{v}_i\right|\cdot\left|\vec{v}_j\right|}} и

${\vec {c'}}={\begin{pmatrix}{\vec {v}}_{i}&{\vec {v}}_{j}\end{pmatrix}}{\begin{pmatrix}\left|{\vec {v}}_{j}\right|\\\left|{\vec {v}}_{i}\right|\end{pmatrix}}$.

Тогда Невозможно разобрать выражение (MathML с запасными SVG или PNG (рекомендуется для современных браузеров и инструментов повышения доступности): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «https://wikimedia.org/api/rest_v1/»:): {\displaystyle \vec{c}=\frac{\tau'\vec{c'}}{\tau'\left|\vec{c'}\right|}=\frac{\vec{c'}}{\left|\vec{c'}\right|}} и $\alpha ={\frac {1}{2}}\arccos {\frac {{\vec {v'}}_{i}\cdot {\vec {v'}}_{j}}{\left|{\vec {v'}}_{i}\right|\cdot \left|{\vec {v'}}_{j}\right|}}$.