Skip to content

Tuorui "v1ncent19" Peng

En voyage dans l'espace de Hilbert.

Mathematics & Statistics4 min read中文

RKHS in a Dual Space delta function Perspective

Recently I was introduced a new point of view of RKHS by my advisor. I think it's quite interesting and worth sharing. It provides a new understanding of how RKHS is constructed and, it's surprisingly related to the Dirac delta function in physics (and it does provide me some more 'maths' understanding of the delta function).

The original note could be found in Scattered Data Approximation by Holger Wendland. Thus here I directly write my understanding of the note in Chinese, with some of my own understanding added.

For some preliminary knowledge about RKHS see my note in RKHS theory.

Understanding of Dirac δ\delta

从前我们(在物理的视角中)对 Dirac δ\delta function 的理解是作为一个单点位置的激发函数,say (如果不带角标地版本默认为是 0 处的 δ\delta 函数)

δx0()δ(x0),s.t.δx0(x)={,x=x00,xx0,δx0(x)dx=1\begin{align} \delta_{x_0}(\, \cdot \, )\equiv \delta(\, \cdot \, -x_0),\quad s.t. \delta_{x_0}(x)=\begin{cases} \infty, & x=x_0 \\ 0, & x\neq x_0 \end{cases},\quad \int_{-\infty}^{\infty} \delta_{x_0}(x)dx=1 \end{align}

Well 如果要更严格化地讲 Dirac δ\delta 所谓一个“函数”定义,我们一般会用一个函数列逼近,比如一些典型的例子有

δn(x)=n1[1n,1n](x)δn(x)=nπsinnxnx\begin{align} \delta _n(x) = n\mathbb{1}_{[-\frac{1}{n},\frac{1}{n}]}(x) \\ \delta _n(x) = \dfrac{ n }{ \pi } \dfrac{ \sin nx }{ nx } \end{align}

Dirac δ\delta as a dual vector

一种更直接的方式是将 Dirac δx0\delta_{x_0} 看作是一个函数空间对偶空间中的一个元素,它对函数空间中的元素 ff 的作用效果是

δx0,f=f(x0)\begin{align} \langle \delta_{x_0},f \rangle = f(x_0) \end{align}

这里的函数空间(和它的对偶空间)只要仿照广义相对论里面的张量逆变矢量空间(以及对偶的协变矢量空间)的概念理解就可以了。

  • Review:

    逆变矢量空间我们记为VV,里面的元素vVv\in V就是我们平时熟悉使用的各种物理量,比如位移矢量r\vec{r},速度矢量v\vec{v},电场强度矢量E\vec{E},等等都是逆变的,在广相里我们会用抽象指标记号,say vμv^\mu 表示它。在VV的基础上我们可以将其对偶空间VV^*定义为所有VV上(线性)函数的集合,即

    V:={ω:VR}\begin{align} V^*:= \{ \omega : V\mapsto \mathbb{R} \} \end{align}

    也就是VV^*中的元素将vVv\in V映射成实数,这也就是广相里协变矢量的定义,say ων\omega_\nu。如果进一步有了度规gμνg_{\mu \nu }作为同构映射之后,我们可以将VV中的元素vμv^\mu映射到VV^*中的元素ων\omega_\nu,并将之记为它的对偶元素

    vνgμνvμ\begin{align} v_\nu \equiv g_{\mu \nu }v^\mu \end{align}

    当然在没有度规的情况下这种同构映射是 arbitrary 的,度规的存在才引进了这样一种“独特”的对应关系(同构映射),也让我们在广相里可以方便地使用某个量的“一对”矢量和对偶矢量。

这个时候将上面的矢量空间换成函数空间 V={f}V=\{f\} 就能发现 δx0\delta _{x_0} 更好的一种理解方式是:它是 VV^* 中的元素,将任意 fVf\in V 映射到 f(x0)f(x_0)

RKHS in a dual space delta function PoV

可以简单回顾一下所谓的 Reproducing Kernel K(,)K(\, \cdot \, ,\, \cdot \, ): 一个对称非负的二元函数,再生表现为

K(,x0),fHK=f(x0)\begin{align} \langle K(\, \cdot \, ,x_0),f \rangle _{\mathcal{H}_K} = f(x_0) \end{align}

对比一下 δx0\delta _{x_0} (此时作为对偶矢量)的效果立刻就能看出相似之处:

δx0[f]=f(x0)\begin{align} \delta _{x_0}[f] = f(x_0) \end{align}

这里的一个直接解释是 Riesz 表示定理(这里我直接给出内容了,不会证)

  • 对于矢量空间 VV,对于任意的 ωV \omega \in V^*,都存在唯一的 ϕωV \phi _\omega \in V w.r.t. ϕωV=ωV\Vert \phi_\omega \Vert_{V} = \Vert \omega \Vert_{V^*},使得 vV\forall v\in V ω(v)=ϕω,v\begin{align} \omega (v) = \langle \phi_\omega ,v \rangle \end{align}

所以这里的意思完全就是,在上述的 Hilbert Space HK\mathcal{H}_K 和对应的 Kernel K(,)K(\, \cdot \, ,\, \cdot \, ) 下,我们有

δx0=K(,x0),HK\begin{align} \delta _{x_0} = \langle K(\, \cdot \, ,x_0),\, \cdot \, \rangle _{\mathcal{H}_K} \end{align}

(这里一个方便的观点是,given x0x_0, K(,x0)HKK(\, \cdot \, ,x_0) \in \mathcal{H}_K)anyway i.e. 这里 Reproducing Kernel 的作用也是通过 Riesz 表示定理来构造了一个F:HKHKF: \mathcal{H}_K^*\mapsto \mathcal{H}_KF(δx0)=K(,x0)F(\delta _{x_0}) = K(\, \cdot \, ,x_0)。借用这个映射我们可以将内积也映射过去

δx0,δx1HK=K(,x0),K(,x1)HK\begin{align} \langle \delta _{x_0},\delta _{x_1} \rangle _{\mathcal{H}_K^*} = \langle K(\, \cdot \, ,x_0),K(\, \cdot \, ,x_1) \rangle _{\mathcal{H}_K} \end{align}

How is it useful

这种观点能够方便地帮助我们理解 RKHS 的一些性质,比如说

  • 正定性。n,{αi,xi}i=1n\forall n, \{\alpha_i,x_i \}_{i=1}^n i,j=1nαiαjK(xi,xj)=i,j=1nαiαjK(,xi),K(,xj)HK=i,j=1nαiαjδxi,δxjHK=i=1nαiδxi,j=1nαjδxjHK=i=1nαiδxiHK20\begin{align} \sum_{i,j=1}^n\alpha _i\alpha _jK(x_i,x_j)=& \sum_{i,j=1}^n \alpha _i\alpha _j \langle K(\, \cdot \, ,x_i),K(\, \cdot \, ,x_j) \rangle _{\mathcal{H}_K} \\ =& \sum_{i,j=1}^n \alpha _i\alpha _j \langle \delta _{x_i},\delta _{x_j} \rangle _{\mathcal{H}_K^*} \\ =& \langle \sum_{i=1}^n \alpha _i \delta _{x_i},\sum_{j=1}^n \alpha _j \delta _{x_j} \rangle _{\mathcal{H}_K^*} \\ =& \Vert \sum_{i=1}^n \alpha _i \delta _{x_i} \Vert _{\mathcal{H}_K^*}^2 \geq 0 \end{align}
  • 核再生 K(x,),K(,y)HK=δy(K(x,))=K(x,y)\begin{align} \langle K(x, \, \cdot \, ) , K(\, \cdot \, ,y) \rangle _{\mathcal{H}_K} =& \delta _y(K(x,\, \cdot \, )) = K(x,y)\\ \end{align}