今天看D.M.J. Tax的Data description toolbox (dd_tool)手冊,裡面有一段簡單的測試代碼:
- x=oc_set(gendatb([50,10]),'1')
- scatterd(x,'legend')
- w=svdd(target_class(x),0.1,8)
- plotc(w)
- w=svdd(x,0.1,8)
- plotc(w)
一開始會出現上篇《PRtools模式識別工具箱關於version運行出錯》 問題,得以解決。
最終運行圖像如下:
- x=oc_set(gendatb([500,100]),'1')
- scatterd(x,'legend')
-
- % svdd mapping
- % FRACREJ:棄真(第I類)錯誤10%
- % SIGMA:8
- w=svdd(x,0.1,8)
- plotc(w,'r*')
-
- % FRACREJ:棄真(第I類)錯誤10%
- % SIGMA:4
- w=svdd(x,0.1,4)
- plotc(w,'r--')
-
- % gauss mapping
- % FRACREJ:棄真(第I類)錯誤10%
- w=gauss_dd(x,0.1)
- plotc(w,'b--')
結果圖: