再生核希尔伯特空间(RKHS)
中级网管
2021年06月04日 19:21

1.由函数到向量

首先函数是可以被看做向量的,例如定义在%5Ba%2Cb%5D上的函数f,如果取%5CDelta%20x为间隔得到点a%2C%20x_1%2Cx_2%2C...%2Cx_n%2Cb,那么就可以把函数变为一个向量(f(a)%2Cf(x_1)%2C...%2Cf(x_n)%2Cf(b)),如果%5CDelta%20x趋向于0那么显然可以把该函数作为一个无穷维向量看待。同理当x为向量时该结论也成立。

这样我们就可以定义函数内积%3Cf%2Cg%3E%3D%20%5Clim_%7B%5CDelta%20x%5Cto%200%7D%5Cdisplaystyle%5Csum_%7Bi%7Df(x_i)g(x_i)%5CDelta%20x%3D%5Cint%20f(x)g(x)dx

注意到这里跟向量内积相比多加了%5CDelta%20x项,其原因在于对于向量来说,其维数是离散的,而对于函数来说是连续的,因此多加了%5CDelta%20x来标准化。感觉也可以理解为在%5CDelta%20x半径的邻域上函数f,g取值相同。

2.向量基到函数基

线性空间中所有向量都能由其基的线性组合表示,设R%5Ep维线性空间上的一组基为%5C%7Be_i%5C%7D,则对任意向量都有a%3D%5Cdisplaystyle%5Csum_%7Bi%3D1%7D%5Ep%5Clambda_ie_i,由%3Ca%2Ce_i%3E%3D%3C%5Cdisplaystyle%5Csum_%7Bj%3D1%7D%5Ep%5Clambda_je_j%2Ce_i%3E%3D%5Clambda_i%3Ce_i%2Ce_i%3E%3D%5Clambda_i%7Ce_i%7C%5E2,%5Clambda_i%3D%5Cfrac%7B1%7D%7B%7Ce_i%7C%5E2%7D%3Ca%2C%20e_i%3E.

同理如果取得一组函数基%5C%7Bf_i(x)%5C%7D,则对任意函数都有

g(x)%3D%5Cdisplaystyle%5Csum_i%5Cmu_if_i(x).

下面介绍傅里叶例子。

引入一组基函数%5C%7Bh_p%3De%5E%7Bi2%5Cpi%20px%2FT%7D%5C%7D_%7Bp%3D-%5Cinfty%7D%5E%7B%5Cinfty%7D,定义域为%5B0%2CT%5D,其中p为整数,这组基构成一个函数空间,且任意定义在%5B0%2CT%5D上的函数可以由其线性表示

f%3D%5Cdisplaystyle%5Csum_pc_ph_p.

易证基正交性

%3Ch_p%2Ch_q%3E%3D%5Cint_0%5ETh_p%5Cbar%7Bh%7D_qdx%3D%5Cint_0%5ETe%5E%7Bi2%5Cpi%20px%2FT%7De%5E%7B-i2%5Cpi%20qx%2FT%7D%3D0. 

注意到这里因为是复领域,函数内积定义有所变化。

基的大小

%7Ch_p%7C%5E2%3D%3Ch_p%2Ch_p%3E%3DT.

故

c_p%3D%5Cfrac%7B1%7D%7BT%7D%3Cf%2Ch_p%3E%3D%5Cfrac%7B1%7D%7BT%7D%5Cint_0%5ETf(x)e%5E%7B-i2%5Cpi%20px%2FT%7Ddx.

3.核函数

定义函数K(x%2Cy),若满足:

K(x%2Cy)%3DK(y%2Cx)

%5Cint%20f(x)K(x%2Cy)f(y)dxdy%5Cge0,f为任意函数

则称其为核函数。

(条件2另一种表述为K对应的Gram矩阵为半正定的。)

可以把核函数看成是一个二维无穷大矩阵,类比于矩阵的特征分解,可以找到一组特征值及其对应的特征函数,这些特征函数作为一组正交基可以构成一个函数空间。

设%5Clambda为一特征值,%5Cpsi(x)为其对应的特征方程,则有如下关系成立

%5Cint%20K(x%2Cy)%5Cpsi(x)dx%3D%5Clambda%20%5Cpsi(y)

对于任意特征值%5Clambda_1%5Cneq%5Clambda_2及其对应的特征函数%5Cpsi_1%2C%5Cpsi_2有

%5Cbegin%7Baligned%7D%0A%5Cint%5Clambda_1%5Cpsi_1(x)%5Cpsi_2(x)dx%20%26%20%3D%5Cint%5Cint%20K(y%2Cx)%5Cpsi_1(y)dy%5Cpsi_2(x)dx%20%5C%5C%26%3D%20%5Cint%5Cint%20K(x%2Cy)%5Cpsi_2(x)dx%5Cpsi_1(y)dy%20%5C%5C%26%3D%5Cint%20%5Clambda_2%5Cpsi_2(y)%5Cpsi_1(y)dy%0A%5Cend%7Baligned%7D

故

%3C%5Cpsi_1%2C%5Cpsi_2%3E%3D%5Cint%5Cpsi_1(x)%5Cpsi_2(x)dx%3D0.

其正交性得证.

与矩阵特征分解类似地有

K(x%2Cy)%3D%5Cdisplaystyle%5Csum_i%5Clambda_i%5Cpsi_i(x)%5Cpsi_i(y).

这里注意下i代表的是第几个基函数,%5Cpsi_i应被看作函数转化的无穷向量,x代表的是该函数在x处的取值,

这叫做Mercer理论. 这无穷个正交的特征函数可以作为一组基构成一个函数空间.

4.再生核希尔伯特空间

上面介绍的一组基可以构成一个再生核希尔伯特空间H,这个空间需要注意的是一个完备的定义了内积的线性空间,并且它的核具有再生性。下面主要介绍核的再生性。

以%5C%7B%5Csqrt%7B%5Clambda_i%7D%5Cpsi_i%5C%7D_%7Bi%3D1%7D%5E%5Cinfty作为一组正交基构成一个希尔伯特空间。那么对空间中的任意函数f都有

f%3D%5Cdisplaystyle%5Csum_if_i%5Csqrt%7B%5Clambda_i%7D%5Cpsi_i

因此f可以表示为向量f%3D(f_1%2Cf_2%2C...)_H%5ET.

由上面3中的式子有

K(x%2C%5Ccdot)%3D%5Cdisplaystyle%5Csum_i%5Clambda_i%5Cpsi_i(x)%5Cpsi_i

右边在空间H中可以照上面f的表示一样表示为一个向量

K(x%2C%5Ccdot)%3D(%5Csqrt%7B%5Clambda_1%7D%5Cpsi_1(x)%2C%5Csqrt%7B%5Clambda_2%7D%5Cpsi_2(x)%2C...)%5ET.

因此有

%3CK(x%2C%5Ccdot)%2CK(y%2C%5Ccdot)%3E_H%3D%5Cdisplaystyle%5Csum_i%5Clambda_i%5Cpsi_i(x)%5Cpsi_i(y)%3DK(x%2Cy).

这就是核的再生性。

5.它的作用

如果我们定义一种映射%5Cphi(x)%3DK(x%2C%5Ccdot)%3D(%5Csqrt%7B%5Clambda_1%7D%5Cpsi_1(x)%2C%5Csqrt%7B%5Clambda_2%7D%5Cpsi_2(x)%2C...)%5ET.

那么有%3C%5Cphi(x)%2C%5Cphi(y)%3E_H%3DK(x%2Cy).

由这个结果可以看出我们并不需要清楚地知道这个映射%5Cphi和特征空间的基,对于一个对称半正定函数K,必定存在一种映射和特征空间使得如上式子成立,这也叫做kernel-trick.

对于分类问题来说非线性问题需要进行高维转换,优化过程中对偶表示产生了内积,这种情况下就可以运用此法.

6.补充

上面介绍的是由核函数推出相应的再生核希尔伯特空间.

实际上标准定义为

设有一集合X及一希尔伯特空间H,它的元素是X上的函数f%3AX%5Cto%20R. 我们称H为再生核希尔伯特空间:若存在一函数k%3AX%5Cto%20R%20满足:

f(x)%3D%3Cf%2Ck(%5Ccdot%2Cx)%3E_H%2C%5Cforall%20f%5Cin%20H.

且

H%3Dspan%5C%7Bk(%5Ccdot%2Cx)%3Ax%5Cin%20X%5C%7D.

上面推导并不严密,关于Merce理论方面还得再看看更具体的数学推导。

参考及推荐

[1]http://songcy.net/posts/story-of-basis-and-kernel-part-1/

[2]http://songcy.net/posts/story-of-basis-and-kernel-part-2/

[3]https://qddmj.cn/rkhs.htm

[4]https://zhuanlan.zhihu.com/p/63979860

[5]https://blog.csdn.net/ChangHengyi/article/details/80577318