Eine radiale Basisfunktion (RBF) ist eine reelle Funktion ${\displaystyle \varphi }$, deren Wert nur vom Abstand zum Ursprung abhängt, so dass ${\displaystyle \varphi (\mathbf {x} )=\varphi (\|\mathbf {x} \|)}$. Der Name kommt daher, dass die Funktion nach dieser Definition radialsymmetrisch ist und ferner diese Funktionen als Basisfunktionen einer Approximation verwendet werden. Allgemeiner kann man den Abstand zu einem Punkt c betrachten, der Zentrum genannt wird, so dass ${\displaystyle \varphi (\mathbf {x} ,\mathbf {c} )=\varphi (\|\mathbf {x} -\mathbf {c} \|)}$.

Eigenschaften

Jede Funktion ${\displaystyle \varphi (x)}$ für die ${\displaystyle \varphi (\mathbf {x} )=\varphi (\|\mathbf {x} \|)}$ gilt ist eine RBF. Als Norm ${\displaystyle \|\cdot \|}$ wird gewöhnlich die euklidische Norm gewählt, je nach Problem sind auch andere Normen möglich. Nimmt man beispielsweise die von der Lukaszyk-Karmowski-Metrik induzierte Norm, so ist es bei manchen RBFs möglich, Probleme mit schlecht konditionierten Matrizen zu umgehen, die zur Berechnung der Koeffizienten, die bei Approximation durch RBFs auftreten, erforderlich sind. Da es sich bei dem Radius um eine skalare Größe handelt, ist das Verhalten bezüglich des Aufwandes im höherdimensionalen Raum gutartig.

Typen radialer Basisfunktionen

Stückweise polynomielle RBF (${\displaystyle R_{n}}$)	${\displaystyle |r|^{n}}$, für ${\displaystyle n}$ ungerade
Thin Plate Spline (${\displaystyle TPS_{n}}$)	${\displaystyle |r|^{n}\ln(|r|)}$, für ${\displaystyle n}$ gerade
Multiquadric RBF (MQ)	${\displaystyle {\sqrt {1+(ar)^{2}}}}$
Inverse multiquadric RBF (IMQ)	${\displaystyle {\frac {1}{\sqrt {1+(ar)^{2}}}}}$
Inverse quadric RBF (IQ)	${\displaystyle {\frac {1}{1+(ar)^{2}}}}$
Polyharmonic RBF	${\displaystyle (-1)^{\left\lceil {\tfrac {n}{2}}\right\rceil }|r|^{n}\ln(|r|)}$
gaußsche RBF (GS)	${\displaystyle e^{-(ar)^{2}}}$

Approximation mit Hilfe von RBFs

Typischerweise werden Linearkombinationen von radialen Basisfunktionen zur Approximation von Funktionen genutzt. Hierbei wird die zu approximierende Funktion ${\displaystyle f(x)}$ durch eine Summe von ${\displaystyle n}$ radialen Basisfunktionen angenähert, die verschiedene Zentren ${\displaystyle \mathbf {x_{i}} }$ haben und durch die Koeffizienten ${\displaystyle \lambda _{i}}$ gewichtet sind. Weiterhin ist zu beachten, dass die Funktionswerte an den Zentren exakt eingehalten werden.

Im Allgemeinen dienen radiale Basisfunktionen zur Annäherung stetiger Funktionen, die von zwei oder mehr Variablen abhängen (multivariaten Funktionen), durch Linearkombinationen von Ausdrücken, die auf einer einzelnen Funktion einer Variablen basieren. Ist ${\displaystyle \varphi \colon \mathbb {R} ^{+}\to \mathbb {R} }$ eine radiale Basisfunktion, kann man eine multivariate Funktion wie folgt definieren:

${\displaystyle \Phi (\mathbf {x} )=\varphi (\|\mathbf {x} \|_{2})}$ für alle reellen Vektoren ${\displaystyle \mathbf {x} =(x_{1},x_{2},\ldots ,x_{n})\in \mathbb {R} ^{n}}$

Dabei ist ${\displaystyle \|\mathbf {x} \|_{2}={\sqrt {x_{1}^{2}+x_{2}^{2}+\cdots +x_{n}^{2}}}}$ der euklidische Abstand vom Koordinatenursprung. Im zweidimensionalen Fall bedeutet das, dass die Funktion ${\displaystyle \Phi }$ für jeden Punkt auf einem Kreis, dessen Mittelpunkt der Koordinatenursprung ist, denselben Wert hat. Dasselbe gilt für dreidimensionale Kugeln und höherdimensionalen Äquivalente (Hyperkugeln).

Wenn eine Interpolation auf einer Funktion ${\displaystyle f\colon \mathbb {R} ^{n}\to \mathbb {R} }$ durchgeführt wird, die an endlich vielen Punkten ${\displaystyle \mathbf {x_{i}} \in \mathbb {R} }$ bekannt ist, wird eine Menge von Basisfunktionen benötigt. Diese kann erzeugt werden, indem die radiale Basisfunktion ${\displaystyle \varphi }$ von den gegebenen Punkten verschoben wird. Die Funktion ${\displaystyle \varphi (\|\mathbf {x} -\mathbf {x_{i}} \|_{2})}$ hat dieselbe Gestalt wie die Funktion ${\displaystyle \Phi }$, aber bezogen auf den Punkt ${\displaystyle \mathbf {x_{i}} }$ statt dem Koordinatenursprung. Für jeden Punkt ${\displaystyle \mathbf {x_{i}} }$ erhält man eine neue Basisfunktion. Diese Basisfunktionen können verwendet werden, um eine Näherungsfunktion zu erzeugen:

${\displaystyle s(\mathbf {x} ):=\sum _{i=1}^{n}\lambda _{i}\cdot \varphi (\|\mathbf {x} -\mathbf {x_{i}} \|_{2})}$

Dabei ist jedes ${\displaystyle \lambda _{i}}$ eine reelle Zahl. Wenn die Bedingung ${\displaystyle s(\mathbf {x_{i}} )=f(\mathbf {x_{i}} )}$ für alle Indexe ${\displaystyle i}$ erfüllt sein soll, sind die Koeffizienten ${\displaystyle \lambda _{i}}$ die Lösungen eines linearen Gleichungssystems.

function [lam,coreFunction] = RBFcalcCoeff(a,xk,fk)

coreFunction = @(x) exp(-(a.*x).^2);   %%% definieren der RBF (phi(x))

N = length(xk);                        %%% N ausgewählte Stützstellen

A = zeros(N,N);                        %%% Initialisierung der Matrix A

for i1 = 1:N                           %%% Schleife über alle Zeilen
  for i2 = i1:N                        %%% Schleife über alle Spalten
      r = sqrt((xk(i1)-xk(i2)).^2);    %%% Bestimmung des Radius
      A(i1,i2) = coreFunction(r);      %%% Funktionswert der Matrix
      A(i2,i1) = A(i1,i2);
  end
end

lam=A\fk;                              %%% Lösung des Gleichungssystems

Positive Definitheit

Eine radiale Basisfunktion ${\displaystyle \varphi \colon \;\mathbb {R} ^{d}\to [0,\infty )}$ ist positiv definit, wenn für alle endlichen Mengen ${\displaystyle X:=\{x_{1},x_{2},\ldots ,x_{m}\}}$ aus endlich vielen Punkten ${\displaystyle x_{1},x_{2},\ldots ,x_{m}\in \mathbb {R} ^{d}}$ die symmetrischen Matrizen ${\displaystyle A_{X}:=(\varphi (\|x_{j}-x_{k}\|_{2}))_{1\leq j,k\leq m}}$ positiv definit sind.

Anwendungen

Approximationsmethoden dieser Art werden zur Modellierung nichtlinearer Systeme (mit ausreichend einfachem chaotischem Verhalten), zur 3D-Rekonstruktion in der Computergrafik (beispielsweise hierarchische RBFs) oder zur Erstellung von Antwortflächen in der Optimierung eingesetzt. Weitere Anwendungsbereiche von RBFs sind das Lösen von partiellen Differentialgleichungen (z. B. mittels gitterfreier Kollokation) oder Transformationen bei der Bildregistrierung.

Literatur

• Martin D. Buhmann: Radial Basis Functions: Theory and Implementations. Cambridge University Press, 2003
• E. Larsson, B. Fornberg: A Numerical Study of some Radial Basis Function based Solution Methods for Elliptic PDEs. In: Computers & Mathematics with Applications, Volume 46, 2003
• H. Wendland: Scattered Data Approximation. Cambridge University Press, 2005

Weblinks

• Martin Buhmann, Janin Jäger, On radial basis functions, Oberwolfach Snapshots 2019

Einzelnachweise

1. ↑ Martin Buhmann, Janin Jäger: On radial basis functions
2. ↑ R. Schaback, Georg-August-Universität Göttingen: A Practical Guide to Radial Basis Functions