動径基底関数がカーネルであることを証明する方法は?


35

動径基底関数がカーネルであることを証明する方法は?私の知る限り、これを証明するためには、次のいずれかを証明する必要があります。k(x,y)=exp⁡(−||x−y||2)2σ2)

  1. ベクトルのセットマトリックス =は半正定です。x1,x2,...,xnK(x1,x2,...,xn)(k(xi,xj))n×n

  2. =ようなマッピングを提示できます。Φk(x,y)⟨Φ(x),Φ(y)⟩

何か助け?


1
より明確にリンクするために、この質問で機能マップについても説明します。特に、テイラー級数とRKHSとダグラスが提供する埋め込みの一般的なバージョンの両方を議論する鉱山シリーズに基づいたマーククレーセンの答え。L2
— ドゥーガル

回答:


26

禅は、ここで方法1を使用した方法2:地図中心とする球対称ガウス分布にヒルベルト空間における。これが正確に機能するためには、標準偏差と定数係数を微調整する必要があります。たとえば、ある次元では、xxL2

∫−∞∞exp⁡[−(x−z)2/(2σ2)]2πσexp⁡[−(y−z)2/(2σ2)2πσdz=exp⁡[−(x−y)2/(4σ2)]2πσ.

したがって、標準偏差を使用し、ガウス分布をスケーリングしてを取得します。この最後の再スケーリングは、正規分布のL ^ 2ノルムが一般に1ではないために発生します。 K(X、Y)=⟨Φ(X)、Φ(Y)⟩L21σ/2k(x,y)=⟨Φ(x),Φ(y)⟩L21


2
@ Zen、Douglas Zare:すばらしい回答をありがとう。公式の回答を今すぐ選択するにはどうすればよいですか?
— レオ

23

方法1を使用します。方法2を使用して、Douglas Zareの回答を確認してください。

とき、私はケースを証明する実数であるので、K (X 、Y )= EXP (- (X - Y )2 / 2 σ 2)。一般的なケースは、同じ議論から必要な変更を加えたものであり、実行する価値があります。バツ、yk(x,y)=exp⁡(−(x−y)2/2σ2)

一般性を失うことなく、その仮定。σ2=1

を書きます。ここで、h (t )= exp ( − t 2k(バツ、y)=h(バツ−y)は、分布を持つ確率変数Zの特性関数です。

h(t)=exp⁡(−t22)=E[e私tZ]
ZN(0、1)

実数および場合、 これはが正の半正定関数、つまりカーネルあることを伴います。a 1、… 、a n n ∑ j 、k = 1 a jバツ1、…、バツna1、…、an

∑j,k=1najakh(xj−xk)=∑j,k=1najakE[ei(xj−xk)Z]=E[∑j,k=1najeixjZake−ixkZ]=E[|∑j=1najeixjZ|2]≥0,
k

この結果をより一般的に理解するには、Bochnerの定理をチェックしてください:http : //en.wikipedia.org/wiki/Positive-definite_function


2
これは、2つの注意事項がある正しい方向の良い出発点です。(a)は表示された期待値と等しくありません(指数の符号を確認してください)。(b)これは、xとyはスカラーであり、ベクトルではありません。博覧会は素晴らしく、きれいであり、これらの小さなギャップをすぐに埋めることができると確信しているので、その間は賛成しました。:h(t)xy
— 枢機

1
Tks!ここで急いでいます。:
— 禅

1
申し訳ありませんが、ここで必要な変更を加えた方法はわかりません。フォームに渡す前にノルムを開発すると、製品を得て、製品と合計を交換することはできません。そして、私は単に、いい形を得るためにh形に渡した後、ノルムをどのように展開するかわかりません。私をそこに少し連れて行ってもらえますか?:)h
— Alburkerk

23

さまざまな目的のために、3番目の方法を追加します。pdカーネルを作成することが知られている一連の一般的な手順からカーネルを構築します。レッツ以下のカーネルとのドメイン表しφ特徴マップを。Xφ

  • スケーリング: 場合 PDカーネルがある、そうであるγ κ任意の定数のγ > 0で。κγκγ>0

    証明:がκの特徴マップの場合、√φκのための有効な特徴マップですγκ。γφγκ

  • 合計: 場合は及びκ 2は、PDカーネルです、そうですκ 1 + κ 2。κ1κ2κ1+κ2

    証明:連結し、特徴がマップ及びφ 2を取得するために、X ↦ [ φ 1(X )、φ 2(X ) ]。φ1φ2x↦[φ1(x)φ2(x)]

  • 限界: 場合 PDカーネルであり、κ (X 、Y ):= LIM N → ∞ κ N(X 、Yは)全てに対して存在するX 、Y、次いで、κは PDです。κ1,κ2,…κ(x,y):=limn→∞κn(x,y)x,yκ

    証明:それぞれについてとすべて{ (X I、C I)} mはI = 1 ⊆ X × R我々は持っているΣはmのiは= 1つの、C I κ nと(X I、X J)C J ≥ 0。制限をn → ∞とすると、κに対して同じプロパティが得られます。m,n≥1{(xi,ci)}i=1m⊆X×R∑i=1mciκn(xi,xj)cj≥0n→∞κ

  • 製品: 場合及びκ 2は、 PDのカーネルであり、そうであるGは(X 、Y )= κ 1(X 、Y )κ1κ2。g(x,y)=κ1(x,y)κ2(x,y)

    証明:Schurの積定理からすぐに続きますが、Schölkopfand Smola(2002)は次の素晴らしく初歩的な証明を与えます。してみましょう 独立です。したがって、 C o v(V i W i、V j W j)= C o v(V i、V j)

    (V1,…,Vm)∼N(0,[κ1(xi,xj)]ij)(W1,…,Wm)∼N(0,[κ2(xi,xj)]ij)
    共分散行列は、の共分散行列考慮ので、PSDでなければならない(V 1 、W 1、... 、V N 、W Nは)それを証明しています。
    Cov(ViWi,VjWj)=Cov(Vi,Vj)Cov(Wi,Wj)=κ1(xi,xj)κ2(xi,xj).
    (V1W1,…,VnWn)
  • パワーズ: 場合 PDカーネルが、そうであるκ N(X 、Y ):= κ (X 、Y )nは任意の正の整数を、nは。κκn(x,y):=κ(x,y)nn

    証明:「products」プロパティから直接。

  • 指数: 場合 PDカーネルが、そうであるE κ(X 、Y ):= EXP (κ (X 、Y ))。κeκ(x,y):=exp⁡(κ(バツ、y))

    eκ(x,y)=limN→∞∑n=0N1n!κ(x,y)n

  • κf:X→Rg(x,y):=f(x)κ(x,y)f(y)

    x↦f(x)φ(x)

k(x,y)=exp⁡(−12σ2‖x−y‖2)=exp⁡(−12σ2‖x‖2)exp⁡(1σ2xTy)exp⁡(−12σ2‖y‖2).
κ(x,y)=xTy1σ2x↦exp⁡(−12σ2‖x‖2)
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.