1. 程式人生 > >svm核函式的理解和選擇

svm核函式的理解和選擇

特徵空間的隱式對映:核函式

    咱們首先給出核函式的來頭:在上文中,我們已經瞭解到了SVM處理線性可分的情況,而對於非線性的情況,SVM 的處理方法是選擇一個核函式 κ(⋅,⋅) ,通過將資料對映到高維空間,來解決在原始空間中線性不可分的問題。

    此外,因為訓練樣例一般是不會獨立出現的,它們總是以成對樣例的內積形式出現,而用對偶形式表示學習器的優勢在為在該表示中可調引數的個數不依賴輸入屬性的個數,通過使用恰當的核函式來替代內積,可以隱式得將非線性的訓練資料對映到高維空間,而不增加可調引數的個數(當然,前提是核函式能夠計算對應著兩個輸入特徵向量的內積)。

    線上性不可分的情況下,支援向量機首先在低維空間中完成計算,然後通過核函式將輸入空間對映到高維特徵空間,最終在高維特徵空間中構造出最優分離超平面,從而把平面上本身不好分的非線性資料分開。如圖7-7

所示,一堆資料在二維空間無法劃分,從而對映到三維空間裡劃分

    而在我們遇到核函式之前,如果用原始的方法,那麼在用線性學習器學習一個非線性關係,需要選擇一個非線性特徵集,並且將資料寫成新的表達形式,這等價於應用一個固定的非線性對映,將資料對映到特徵空間,在特徵空間中使用線性學習器,因此,考慮的假設集是這種型別的函式:
    這裡ϕ:X->F是從輸入空間到某個特徵空間的對映,這意味著建立非線性學習器分為兩步:
  1. 首先使用一個非線性對映將資料變換到一個特徵空間F,
  2. 然後在特徵空間使用線性學習器分類。
    而由於對偶形式就是線性學習器的一個重要性質,這意味著假設可以表達為訓練點的線性組合,因此決策規則可以用測試點和訓練點的內積來表示:
    如果有一種方式可以在特徵空間中直接計算內積φ(xi · φ(x),就像在原始輸入點的函式中一樣,就有可能將兩個步驟融合到一起建立一個非線性的學習器,這樣直接計演算法的方法稱為核函式方法:     核是一個函式K,對所有x,z(-X,滿足,這裡φ是從X到內積特徵空間F的對映。

核函式:如何處理非線性資料

    來看個核函式的例子。如下圖所示的兩類資料,分別分佈為兩個圓圈的形狀,這樣的資料本身就是線性不可分的,此時咱們該如何把這兩類資料分開呢(下文將會有一個相應的三維空間圖)?

 

    事實上,上圖所述的這個資料集,是用兩個半徑不同的圓圈加上了少量的噪音生成得到的,所以,一個理想的分界應該是一個“圓圈”而不是一條線(超平面)。如果用 X

1 和 X2 來表示這個二維平面的兩個座標的話,我們知道一條二次曲線(圓圈是二次曲線的一種特殊情況)的方程可以寫作這樣的形式:

    注意上面的形式,如果我們構造另外一個五維的空間,其中五個座標的值分別為 Z1=X1Z2=X21Z3=X2Z4=X22Z5=X1X2,那麼顯然,上面的方程在新的座標系下可以寫作:

    關於新的座標 Z ,這正是一個 hyper plane 的方程!也就是說,如果我們做一個對映 ϕ:R2R5 ,將 X 按照上面的規則對映為 Z ,那麼在新的空間中原來的資料將變成線性可分的,從而使用之前我們推導的線性分類演算法就可以進行處理了。這正是Kernel 方法處理非線性問題的基本思想。

    再進一步描述 Kernel 的細節之前,不妨再來看看這個例子對映過後的直觀例子。當然,你我可能無法把 5 維空間畫出來,不過由於我這裡生成資料的時候就是用了特殊的情形,具體來說,我這裡的超平面實際的方程是這個樣子(圓心在 X2 軸上的一個正圓):

    因此我只需要把它對映到 Z1=X21Z2=X22Z3=X2 這樣一個三維空間中即可,下圖即是對映之後的結果,將座標軸經過適當的旋轉,就可以很明顯地看出,資料是可以通過一個平面來分開的(pluskid:下面的gif 動畫,先用 Matlab 畫出一張張圖片,再用 Imagemagick 拼貼成):

    核函式相當於把原來的分類函式

    對映成:

    而其中的可以通過求解如下 dual 問題而得到的:

    這樣一來問題就解決了嗎?似乎是的:拿到非線性資料,就找一個對映 ,然後一股腦把原來的資料對映到新空間中,再做線性 SVM 即可。不過事實上沒有這麼簡單!其實剛才的方法稍想一下就會發現有問題:在最初的例子裡,我們對一個二維空間做對映,選擇的新空間是原始空間的所有一階和二階的組合,得到了五個維度;如果原始空間是三維,那麼我們會得到 19 維的新空間,這個數目是呈爆炸性增長的,這給 的計算帶來了非常大的困難,而且如果遇到無窮維的情況,就根本無從計算了。所以就需要 Kernel 出馬了。

    不妨還是從最開始的簡單例子出發,設兩個向量,而即是到前面說的五維空間的對映,因此對映過後的內積為:

(公式說明:上面的這兩個推導過程中,所說的前面的五維空間的對映,這裡說的前面便是文中2.2.1節的所述的對映方式,回顧下之前的對映規則,再看那第一個推導,其實就是計算x1,x2各自的內積,然後相乘相加即可,第二個推導則是直接平方,去掉括號,也很容易推出來

    另外,我們又注意到:

二者有很多相似的地方,實際上,我們只要把某幾個維度線性縮放一下,然後再加上一個常數維度,具體來說,上面這個式子的計算結果實際上和對映

     之後的內積的結果是相等的,那麼區別在於什麼地方呢?

  1. 一個是對映到高維空間中,然後再根據內積的公式進行計算;
  2. 而另一個則直接在原來的低維空間中進行計算,而不需要顯式地寫出對映後的結果

    (公式說明:上面之中,最後的兩個式子,第一個算式,是帶內積的完全平方式,可以拆開,然後,通過湊一個得到,第二個算式,也是根據第一個算式湊出來的

    回憶剛才提到的對映的維度爆炸,在前一種方法已經無法計算的情況下,後一種方法卻依舊能從容處理,甚至是無窮維度的情況也沒有問題。

    我們把這裡的計算兩個向量在隱式對映過後的空間中的內積的函式叫做核函式 (Kernel Function) ,例如,在剛才的例子中,我們的核函式為:

核函式能簡化對映空間中的內積運算——剛好“碰巧”的是,在我們的 SVM 裡需要計算的地方資料向量總是以內積的形式出現的。對比剛才我們上面寫出來的式子,現在我們的分類函式為:

    其中 由如下 dual 問題計算而得:

    這樣一來計算的問題就算解決了,避開了直接在高維空間中進行計算,而結果卻是等價的!當然,因為我們這裡的例子非常簡單,所以我可以手工構造出對應於的核函數出來,如果對於任意一個對映,想要構造出對應的核函式就很困難了。

幾個核函式

    通常人們會從一些常用的核函式中選擇(根據問題和資料的不同,選擇不同的引數,實際上就是得到了不同的核函式),例如:

  • 多項式核,顯然剛才我們舉的例子是這裡多項式核的一個特例(R = 1,d = 2。雖然比較麻煩,而且沒有必要,不過這個核所對應的對映實際上是可以寫出來的,該空間的維度是,其中 是原始空間的維度。
  • 高斯核,這個核就是最開始提到過的會將原始空間對映為無窮維空間的那個傢伙。不過,如果選得很大的話,高次特徵上的權重實際上衰減得非常快,所以實際上(數值上近似一下)相當於一個低維的子空間;反過來,如果選得很小,則可以將任意的資料對映為線性可分——當然,這並不一定是好事,因為隨之而來的可能是非常嚴重的過擬合問題。不過,總的來說,通過調控引數,高斯核實際上具有相當高的靈活性,也是使用最廣泛的核函式之一。下圖所示的例子便是把低維線性不可分的資料通過高斯核函式對映到了高維空間:
  • 線性核,這實際上就是原始空間中的內積。這個核存在的主要目的是使得“對映後空間中的問題”和“對映前空間中的問題”兩者在形式上統一起來了(意思是說,咱們有的時候,寫程式碼,或寫公式的時候,只要寫個模板或通用表示式,然後再代入不同的核,便可以了,於此,便在形式上統一了起來,不用再分別寫一個線性的,和一個非線性的)

核函式的本質

        上面說了這麼一大堆,讀者可能還是沒明白核函式到底是個什麼東西?我再簡要概括下,即以下三點:
  1. 實際中,我們會經常遇到線性不可分的樣例,此時,我們的常用做法是把樣例特徵對映到高維空間中去(如上文2.2節最開始的那幅圖