课程作业西瓜书机器学习课后习题：第六章

Posted 2022-09-03 海轰Pro

tags:

篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了课程作业西瓜书机器学习课后习题：第六章相关的知识，希望对你有一定的参考价值。

简介

Hello！
非常感谢您阅读海轰的文章，倘若文中有错误的地方，欢迎您指出～

ଘ(੭ˊᵕˋ)੭
昵称：海轰
标签：程序猿｜C++选手｜学生
简介：因C语言结识编程，随后转入计算机专业，获得过国家奖学金，有幸在竞赛中拿过一些国奖、省奖…已保研
学习经验：扎实基础 + 多做笔记 + 多敲代码 + 多思考 + 学好英语！

唯有努力💪

本文仅记录自己感兴趣的内容

说明

作业要求：每章从课后习题中选取3道题做就可以了

答案来源：题目的解答过程来自于网络，依据个人所学进行了一些修改、总结

仅供参考

6.1

试证明样本空间中任意点x到超平面(w，b)的距离为式(6.2)

6.2

试使用LIBSVM，在西瓜数据集3.0α上分别用线性核和高斯核训练一个SVM，并比较其支持向量的差别。

答：实验代码如下：

from sklearn import svm

X=[
    [ 0.697, 0.46 ],
    [ 0.774, 0.376],
    [ 0.634, 0.264],
    [ 0.608, 0.318],
    [ 0.556, 0.215],
    [ 0.403, 0.237],
    [ 0.481, 0.149],
    [ 0.437, 0.211],
    [ 0.666, 0.091],
    [ 0.243, 0.267],
    [ 0.245, 0.057],
    [ 0.343, 0.099],
    [ 0.639, 0.161],
    [ 0.657, 0.198],
    [ 0.36 , 0.37 ],
    [ 0.593, 0.042],
    [ 0.719, 0.103]
]
y=[1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0]

plt.figure(figsize=(18,6))

for i in range(0, len(X)):
    plt.subplot(1,3,1)
    if y[i] == 1:
        plt.scatter(X[i][0],X[i][1],c='r',marker='*')
    else:
        plt.scatter(X[i][0],X[i][1],c='g',marker='*')


print("-"*20+"线性核"+"-"*20)
clf1=svm.SVC(C = 1,kernel='linear')
print("交叉验证评分",cross_val_score(clf1,X,y,cv=5,scoring='accuracy').mean())
clf1.fit(X,y)
print("支持向量数目:",clf1.n_support_.sum())
print("支持向量:")
print(clf1.support_vectors_)

for i in X:
    res=clf1.predict(np.array(i).reshape(1, -1))
    plt.subplot(1,3,2)
    if res > 0:
        plt.scatter(i[0],i[1],c='r',marker='*')
    else :
        plt.scatter(i[0],i[1],c='g',marker='*')

print("-"*20+"高斯核"+"-"*20)
clf2=svm.SVC(C=1,kernel='rbf')
print("交叉验证评分",cross_val_score(clf2,X,y,cv=5,scoring='accuracy').mean())
clf2.fit(X,y)
print("支持向量数目",clf2.n_support_.sum())
print("支持向量:")
print(clf2.support_vectors_)



for i in X:
    res=clf2.predict(np.array(i).reshape(1, -1))
    plt.subplot(1,3,3)
    if res > 0:
        plt.scatter(i[0],i[1],c='r',marker='*')
    else :
        plt.scatter(i[0],i[1],c='g',marker='*')