next up previous contents
Next: 5.2.2 Cylindrical potential Up: 5.2 Pattern association algorithm Previous: 5.2 Pattern association algorithm


5.2.1 Spherical potential

First, we consider a spherical potential field in feature space. Here, it is not necessary to compute the principal components. All we need is the mean of the data in the feature space,

$\displaystyle \bf\bar{\boldsymbol{\varphi}}$ = $\displaystyle {\frac{{1}}{{n}}}$$\displaystyle \sum_{{i=1}}^{n}$$\displaystyle \varphi$($\displaystyle \bf x_{i}^{}$) . (5.1)

The potential of a point $ \bf z$ in the original space is the squared distance from its mapping $ \varphi$($ \bf z$) to the center $ \bf\bar{\boldsymbol{\varphi}}$,

ES($\displaystyle \bf z$) = ||$\displaystyle \varphi$($\displaystyle \bf z$) - $\displaystyle \bf\bar{\boldsymbol{\varphi}}$||2 = $\displaystyle \varphi$($\displaystyle \bf z$)T$\displaystyle \varphi$($\displaystyle \bf z$) - 2 $\displaystyle \varphi$($\displaystyle \bf z$)T$\displaystyle \bf\bar{\boldsymbol{\varphi}}$ + $\displaystyle \bf\bar{\boldsymbol{\varphi}}^{T}_{}$$\displaystyle \bf\bar{\boldsymbol{\varphi}}$ . (5.2)

Using (5.1), the scalar products can be replaced by the kernel function k (see section 2.4),

ES($\displaystyle \bf z$) = k($\displaystyle \bf z$,$\displaystyle \bf z$) - $\displaystyle {\frac{{2}}{{n}}}$$\displaystyle \sum_{{i=1}}^{n}$k($\displaystyle \bf z$,$\displaystyle \bf x_{i}^{}$) + $\displaystyle {\frac{{1}}{{n^2}}}$$\displaystyle \sum_{{i,j=1}}^{n}$k($\displaystyle \bf x_{i}^{}$,$\displaystyle \bf x_{j}^{}$) . (5.3)

All parts of this equation are known. The last term is constant, and can therefore be omitted. For radial basis functions (section 2.4.3), the first term is also constant, and the potential can be simplified to:

ES($\displaystyle \bf z$) = - $\displaystyle {\frac{{2}}{{n}}}$$\displaystyle \sum_{{i=1}}^{n}$k($\displaystyle \bf z$,$\displaystyle \bf x_{i}^{}$) . (5.4)

This function is proportional to the Parzen window density estimator (Parzen, 1962).


next up previous contents
Next: 5.2.2 Cylindrical potential Up: 5.2 Pattern association algorithm Previous: 5.2 Pattern association algorithm
Heiko Hoffmann
2005-03-22