支持向量機：最大邊際分割超平面?

本案例繪制了二分類時追求最大邊際的分割超平面，使用的算法是線性核下的支持向量機分類器。

輸入：

print(__doc__)

import numpy as np
import matplotlib.pyplot as plt
from sklearn import svm
from sklearn.datasets import make_blobs


# 我們創建40個用來分割的數據點
X, y = make_blobs(n_samples=40, centers=2, random_state=6)

# 擬合模型，并且為了展示作用，并不進行標準化
clf = svm.SVC(kernel='linear', C=1000)
clf.fit(X, y)

plt.scatter(X[:, 0], X[:, 1], c=y, s=30, cmap=plt.cm.Paired)

# 繪制decision function的結果
ax = plt.gca()
xlim = ax.get_xlim()
ylim = ax.get_ylim()

# 創造網格來評估模型
xx = np.linspace(xlim[0], xlim[1], 30)
yy = np.linspace(ylim[0], ylim[1], 30)
YY, XX = np.meshgrid(yy, xx)
xy = np.vstack([XX.ravel(), YY.ravel()]).T
Z = clf.decision_function(xy).reshape(XX.shape)

# 繪制決策邊界和邊際
ax.contour(XX, YY, Z, colors='k', levels=[-1, 0, 1], alpha=0.5,
           linestyles=['--', '-', '--'])
# 繪制支持向量
ax.scatter(clf.support_vectors_[:, 0], clf.support_vectors_[:, 1], s=100,
           linewidth=1, facecolors='none', edgecolors='k')
plt.show()

腳本的總運行時間：（0分鐘0.071秒）