Радиальная базисная функция

редактировать

A радиальная базисная функция (RBF ) - функция с действительным знаком φ {\ textstyle \ varphi}{\ textstyle \ varphi} , значение которого зависит только от расстояния между вводом и некоторой фиксированной точкой, либо исходной точкой, так что φ (x) = φ (‖ Икс ‖) {\ textstyle \ varphi (\ mathbf {x}) = \ varphi (\ left \ | \ mathbf {x} \ right \ |)}{\ textstyle \ varphi (\ mathbf {x}) = \ varphi (\ left \ | \ mathbf {x} \ right \ |)} или другой фиксированный точка c {\ textstyle \ mathbf {c}}{\ textstyle \ mathbf {c}} , называемая центром, так что φ (x) = φ (‖ x - c ‖) {\ textstyle \ varphi (\ mathbf {x}) = \ varphi (\ left \ | \ mathbf {x} - \ mathbf {c} \ right \ |)}{\ textstyle \ varphi (\ mathbf {x}) = \ varphi (\ left \ | \ mathbf {x} - \ mathbf {c} \ right \ |)} . Любая функция φ {\ textstyle \ varphi}{\ textstyle \ varphi} , удовлетворяющая свойству φ (x) = φ (‖ x ‖) {\ textstyle \ varphi (\ mathbf {x}) = \ varphi (\ left \ | \ mathbf {x} \ right \ |)}{\ textstyle \ varphi (\ mathbf {x}) = \ varphi (\ left \ | \ mathbf {x} \ right \ |)} - это радиальная функция. Расстояние обычно составляет евклидово расстояние, хотя иногда используются другие метрики. Они часто используются как набор {φ k} k {\ displaystyle \ {\ varphi _ {k} \} _ {k}}{\ displaystyle \ {\ varphi _ {k } \} _ {k}} , который образует основу для некоторое интересующее функциональное пространство, отсюда и название.

Суммы радиальных базисных функций обычно используются для аппроксимации заданных функций. Этот процесс аппроксимации также можно интерпретировать как простую разновидность нейронной сети ; Это был контекст, в котором они первоначально применялись к машинному обучению в работе Дэвида Брумхеда и Дэвида Лоу в 1988 году, которая была основана на основополагающем исследовании Майкла Дж. Д. Пауэлла 1977 года. RBF также используются в качестве ядра в классификации векторов поддержки. Этот метод оказался достаточно эффективным и гибким, поэтому радиальные базисные функции теперь применяются во множестве инженерных приложений.

Содержание

  • 1 Определение
    • 1.1 Примеры
  • 2 Аппроксимация
  • 3 Сеть RBF
  • 4 См. Также
  • 5 Ссылки
  • 6 Дополнительная литература

Определение

Радиальная функция - это функция φ: [0, ∞) → R {\ textstyle \ varphi: [0, \ infty) \ в \ mathbb {R}}{\ textstyle \ varphi: [0, \ infty) \ to \ mathbb {R}} . В паре с метрикой в ​​векторном пространстве ‖ ⋅ ‖: V → [0, ∞) {\ textstyle \ | \ cdot \ |: V \ to [0, \ infty)}{\ textstyle \ | \ cdot \ |: V \ to [0, \ infty)} a функция φ c = φ (‖ x - c ‖) {\ textstyle \ varphi _ {\ mathbf {c}} = \ varphi (\ | \ mathbf {x} - \ mathbf {c} \ |)}{\ textstyle \ varphi _ {\ mathbf {c}} = \ varphi (\ | \ mathbf {x} - \ mathbf {c} \ |)} называется радиальным ядром с центром в c {\ textstyle \ mathbf {c}}{\ textstyle \ mathbf {c}} . Радиальная функция и связанные с ней радиальные ядра называются радиальными базисными функциями, если для любого набора узлов {xk} k = 1 n {\ displaystyle \ {\ mathbf {x} _ {k} \} _ { k = 1} ^ {n}}{\ displaystyle \ {\ mathbf {x} _ {k} \} _ {k = 1} ^ {n}}

  • Ядра φ x 1, φ x 2,…, φ xn {\ displaystyle \ varphi _ {\ mathbf {x} _ {1}}, \ varphi _ {\ mathbf {x} _ {2}}, \ dots, \ varphi _ {\ mathbf {x} _ {n}}}{\ displaystyle \ varphi _ { \ mathbf {x} _ {1}}, \ varphi _ {\ mathbf {x} _ {2}}, \ dots, \ varphi _ {\ mathbf {x} _ {n}}} линейно независимы (например, φ (r) = r 2 {\ displaystyle \ varphi (r) = r ^ {2}}{\ displaystyle \ varphi (r) = r ^ {2}} in V = R {\ displaystyle V = \ mathbb {R}}{\ displaystyle V = \ mathbb {R}} не является радиальная базисная функция)
[φ (‖ x 1 - x 1 ‖) φ (‖ x 2 - x 1 ‖)… φ (‖ xn - x 1 ‖) φ (‖ x 1 - x 2 ‖) φ (‖ x 2 - x 2 ‖)… φ (‖ xn - x 2 ‖) ⋮ ⋮ ⋱ ⋮ φ (‖ x 1 - xn ‖) φ (‖ x 2 - xn ‖)… φ (‖ xn - xn ‖)] { \ displaystyle {\ begin {bmatrix} \ varphi (\ | \ mathbf {x} _ {1} - \ mathbf {x} _ {1} \ |) \ varphi (\ | \ mathbf {x} _ {2} - \ mathbf {x } _ {1} \ |) \ точки \ varphi (\ | \ mathbf {x} _ {n} - \ mathbf {x} _ {1} \ |) \\\ varphi (\ | \ mathbf {x } _ {1} - \ mathbf {x} _ {2} \ |) \ varphi (\ | \ mathbf {x} _ {2} - \ mathbf {x} _ {2} \ |) \ точки \ varphi (\ | \ mathbf {x} _ {n} - \ mathbf {x} _ {2} \ |) \\\ vdots \ vdots \ ddots \ vdots \\\ varphi (\ | \ mathbf { x} _ {1} - \ mathbf {x} _ {n} \ |) \ varphi (\ | \ mathbf {x} _ {2} - \ mathbf {x} _ {n} \ |) \ точки \ varphi (\ | \ mathbf {x} _ {n} - \ mathbf {x} _ {n} \ |) \\\ end {bmatrix}}}{\ displaystyle {\ begin {bmatrix} \ varphi (\ | \ mathbf {x} _ {1} - \ mathbf {x} _ {1} \ |) \ varphi (\ | \ mathbf {x} _ {2} - \ mathbf {x} _ {1} \ |) \ dots \ varphi (\ | \ mathbf {x} _ {n} - \ mathbf {x} _ {1} \ |) \\\ varphi (\ | \ mathbf {x} _ {1} - \ mathbf {x} _ {2} \ |) \ varphi (\ | \ mathbf {x} _ {2} - \ mathbf {x} _ {2} \ |) \ dots \ varphi (\ | \ mathbf {x} _ {n} - \ mathbf {x} _ {2} \ |) \\ \ vdots \ vdots \ ddots \ vdots \\\ varphi (\ | \ mathbf {x} _ {1} - \ mathbf {x} _ {n} \ |) \ varphi (\ | \ mathbf {x } _ {2} - \ mathbf {x} _ {n} \ |) \ dots \ varphi (\ | \ mathbf {x} _ {n} - \ mathbf {x} _ {n} \ |) \ \\ конец {bmatrix}}} неособое число.

Примеры

Обычно используемые типы радиальных базисных функций включают (запись r = ‖ x - xi ‖ {\ textstyle r = \ left \ | \ mathbf {x} - \ mathbf {x} _ {i} \ right \ |}{\ textstyle r = \ left \ | \ mathbf {x} - \ mathbf {x} _ {i} \ right \ |} и использование ε {\ textstyle \ varepsilon}{\ textstyle \ varepsilon} для указания параметра формы, который может использоваться для масштабирования входных данных радиального ядра):

  • Бесконечно гладкие RBF

Эти радиальные базисные функции взяты из C ∞ (R) {\ displaystyle C ^ {\ infty} (\ mathbb {R})}{\ displaystyle C ^ {\ infty} (\ mathbb {R})} и являются строго положительно определенные функции, требующие настройки параметра формы ε {\ displaystyle \ varepsilon}\ varepsilon

  • Gaussian :φ (r) = e - (ε r) 2 {\ displaystyle \ varphi (r) = e ^ {- (\ varepsilon r) ^ {2}}}{\ displaystyle \ varphi (r) = e ^ {- (\ varepsilon r) ^ {2}}}
A функция Гаусса для нескольких вариантов ε {\ displaystyle \ varepsilon}\ varepsilon . График масштабированного Функция рельефности с несколькими вариантами ε {\ displaystyle \ varepsilon}\ varepsilon .
  • :φ (r) = 1 + (ε r) 2 {\ displaystyle \ varphi (r) = {\ sqrt {1 + (\ varepsilon r) ^ {2}}}}{\ displaystyle \ varphi (r) = {\ sqrt {1+ (\ varepsilon r) ^ {2}}}}
  • :φ (r) = 1 1 + (ε r) 2 {\ displaystyle \ varphi (r) = {\ dfrac {1} {1 + (\ varepsilon r) ^ {2}}}}{\ displaystyle \ varphi (r) = {\ dfrac {1} {1 + (\ varepsilon r) ^ {2}}}}
  • :φ (r) = 1 1 + (ε r) 2 {\ displaystyle \ varphi (r) = {\ dfrac {1} {\ sqrt {1 + (\ varepsilon r) ^ {2}}}}}{\ displaystyle \ varphi (r) = {\ dfrac {1} {\ sqrt {1 + (\ varepsilon r) ^ {2}}}}}
  • Полигармонический сплайн :φ (r) = rk, k = 1, 3, 5, … Φ (г) знак равно рк пер ⁡ (г), к = 2, 4, 6,… {\ displaystyle {\ begin {align} \ varphi (r) = r ^ {k}, k = 1,3, 5, \ dotsc \\\ varphi (r) = r ^ {k} \ ln (r), k = 2,4,6, \ dotsc \ end {выравнивается}}}{\ displayst yle {\ begin {align} \ varphi (r) = r ^ {k}, k = 1,3,5, \ dotsc \\\ varphi (r) = r ^ {k} \ ln (r), k = 2,4,6, \ dotsc \ end {align}}} * Для полигармонической четной степени шлицы (k = 2, 4, 6,…) {\ displaystyle (k = 2,4,6, \ dotsc)}{\ displaystyle (k = 2,4,6, \ dotsc)} , чтобы избежать числовых проблем при r = 0 { \ displaystyle r = 0}r = 0 где ln ⁡ (0) = - ∞ {\ displaystyle \ ln (0) = - \ infty}{\ displaystyle \ ln (0) = - \ infty } , вычислительная реализация часто записывается как φ (r) = rk - 1 ln ⁡ (rr) {\ displaystyle \ varphi (r) = r ^ {k-1} \ ln (r ^ {r})}{\ displaystyle \ varphi (r) = r ^ {k-1} \ ln (r ^ {r})} .

Эти RBF компактно поддерживаются и поэтому не -нуль только в пределах радиуса 1 / ε {\ displaystyle 1 / \ varepsilon}1 / \ varepsilon и, следовательно, имеют разреженные матрицы дифференцирования

φ (r) = {exp ⁡ (- 1 1 - (ε r) 2) для r < 1 ε 0 otherwise {\displaystyle \varphi (r)={\begin{cases}\exp \left(-{\frac {1}{1-(\varepsilon r)^{2}}}\right){\t_dv{ for }}r<{\frac {1}{\varepsilon }}\\0{\t_dv{ otherwise}}\end{cases}}}{\ displaystyle \ varphi (r) = {\ begin {cases} \ exp \ left (- {\ frac {1} {1 - (\ varepsilon r) ^ {2}}} \ right) {\ t_dv {for}} r <{\ frac {1 } {\ varepsilon}} \\ 0 {\ t_dv {else}} \ end {cases}}}

Аппроксимация

Радиальные базисные функции обычно используются для построения аппроксимации функций вида

y (x) Знак равно ∑ я знак равно 1 N вес φ (‖ Икс - хи ‖), {\ Displaystyle у (\ mathbf {x}) = \ сумма _ {я = 1} ^ {N} w_ {я} \, \ varphi (\ left \ | \ mathbf {x} - \ mathbf {x} _ {i} \ right \ |),}{\ displaystyle y (\ mathbf {x}) = \ sum _ {i = 1} ^ {N} w_ {i} \, \ varphi (\ left \ | \ mathbf {x} - \ mathbf {x} _ {i} \ right \ |),}

где аппроксимирующая функция y (x) {\ textstyle y (\ mathbf {x}) }{\ textstyle y (\ mathbf {x})} представлен как сумма N {\ displaystyle N}N радиальных базисных функций, каждая из которых связана с другим центром xi {\ textstyle \ mathbf {x} _ {i}}{\ textstyle \ mathbf {x } _ {я}} , и взвешенные с помощью соответствующего коэффициента wi. {\ textstyle w_ {i}.}{\ textstyle w_ {i}.} Веса wi {\ textstyle w_ {i}}{\ textstyle w_ {i }} можно оценить с помощью матричных методов линейных наименьших квадратов, поскольку аппроксимирующая функция линейна по весам wi {\ textstyle w_ {i}}{\ textstyle w_ {i }} .

Схемы аппроксимации такого типа особенно используются в прогнозировании временных рядов и управление из нелинейных систем, демонстрирующих достаточно простое хаотическое поведение и трехмерную реконструкцию в компьютерной графике (например, иерархический RBF и Деформация пространства позы ).

RBF Network

Две ненормализованные радиальные базисные функции Гаусса в одном входном измерении. Базовые функциональные центры расположены в x 1 = 0,75 {\ textstyle x_ {1} = 0,75}{\ textstyle x_ {1} = 0,75} и x 2 = 3,25 {\ textstyle x_ {2} = 3,25} <211.>Сумма

y (x) = ∑ i = 1 N wi φ (‖ x - xi ‖), {\ displaystyle y (\ mathbf {x}) = \ sum _ {i = 1} ^ {N} w_ {i} \, \ varphi (\ left \ | \ mathbf {x} - \ mathbf {x} _ {i} \ right \ |),}
{\ displaystyle y (\ mathbf {x}) = \ sum _ {i = 1} ^ {N} w_ {i} \, \ varphi (\ left \ | \ mathbf {x} - \ mathbf {x} _ {i} \ right \ |),} также можно интерпретировать как довольно простой однослойный тип искусственной нейронной сети, называемой сетью с радиальными базисными функциями, при этом радиальные базисные функции принимают на себя роль функций активации сети. Можно показать, что любую непрерывную функцию на компактном интервале в принципе можно интерполировать с произвольной точностью суммой этой формы, если достаточно большое число N {\ textstyle N}{\ textstyle N} радиальных базисных функций.

Аппроксимация y (x) {\ textstyle y (\ mathbf {x})}{\ textstyle y (\ mathbf {x})} дифференцируема по весам wi {\ textstyle w_ {i}}{\ textstyle w_ {i }} . Таким образом, веса можно было узнать, используя любой из стандартных итерационных методов для нейронных сетей.

Использование радиальных базисных функций таким образом приводит к разумному подходу к интерполяции при условии, что набор подгонки был выбран таким, что он систематически охватывает весь диапазон (идеально расположены эквидистантные точки данных). Однако без полиномиального члена, который ортогонален радиальным базисным функциям, оценки за пределами подгоночного набора, как правило, работают плохо.

См. Также

Ссылки

Дополнительная литература

  • Харди, Р.Л. (1971). «Мультиквадрические уравнения топографии и других нерегулярных поверхностей». Журнал геофизических исследований. 76 (8): 1905–1915. Bibcode : 1971JGR.... 76.1905H. doi : 10.1029 / jb076i008p01905.
  • Харди, Р.Л. (1990). "Теория и приложения многоквадратично-бигармонического метода, 20 лет открытий, 1968 1988". Комп. Math Applic. 19 (8/9): 163–208. doi : 10.1016 / 0898-1221 (90) 90272-l.
  • Press, WH; Теукольский С.А.; Феттерлинг, штат Вашингтон; Фланнери, Б.П. (2007), «Раздел 3.7.1. Интерполяция радиальной базисной функции», Численные рецепты: Искусство научных вычислений (3-е изд.), Нью-Йорк: Cambridge University Press, ISBN 978-0-521-88068-8
  • Сираяноне, С., 1988, Сравнительные исследования кригинга, многоквадратично-бигармонических и других методов решения проблем минеральных ресурсов, доктор философии. Диссертация, кафедра наук о Земле, Государственный университет Айовы, Эймс, Айова.
  • Sirayanone, S.; Харди, Р.Л. (1995). «Многоквадратично-бигармонический метод, используемый для минеральных ресурсов, метеорологических и других приложений». Журнал прикладных наук и вычислений. 1 : 437–475.
Последняя правка сделана 2021-06-03 05:53:51
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте