9-2-SVM-demo.py 1.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142
  1. import numpy as np
  2. import matplotlib.pyplot as plt
  3. from sklearn.svm import SVC
  4. def main():
  5. # 随机40个点,符号正态分布
  6. np.random.seed(0)
  7. x = np.r_[np.random.randn(20, 2) - [2, 2], np.random.randn(20, 2) + [2, 2]]
  8. y = [0] * 20 + [1] * 20
  9. # 导入svm并训练
  10. clf = SVC(kernel='linear')
  11. clf.fit(x, y)
  12. # 构造超平面
  13. w = clf.coef_[0] # 得到w
  14. a = -w[0] / w[1] # 找到斜率
  15. xx = np.linspace(-5, 5) # -5,5返回均匀间隔的数字
  16. yy = a * xx - (clf.intercept_[0]) / w[1] # clf.intercept_[0]#用来获得截距
  17. b = clf.support_vectors_[0] # 求出过切线的点
  18. yy_down = a * xx + (b[1] - a * b[0]) # 下边界
  19. b = clf.support_vectors_[-1]
  20. yy_up = a * xx + (b[1] - a * b[0]) # 上边界
  21. # matplotlib画图
  22. # 超平面
  23. plt.plot(xx, yy, 'k-')
  24. plt.plot(xx, yy_down, 'k--')
  25. plt.plot(xx, yy_up, 'k--')
  26. # 离得最近的向量
  27. plt.scatter(clf.support_vectors_[:, 0], clf.support_vectors_[:, 1], s=80, facecolors='red')
  28. # 散点
  29. plt.scatter(x[:, 0], x[:, 1], c=y, cmap=plt.cm.Paired)
  30. plt.axis('tight')
  31. plt.show()
  32. if __name__ == "__main__":
  33. main()