zoukankan      html  css  js  c++  java
  • 另一个角度来看:支持向量机(SVM)

    此文,记述数月前,看PRML第四章-线性分类模型时的思考。


    此模型得到的解析解为:


    可以看到:标签信息影响最后的参数的值,进而影响判别边界,接着会对预测产生影响。


    示例:可以看到相同的数据,使用不同的标签信息,得到不同的判别边界。


    训练完参数,使用其做预测,那么:不同的标签选择+各个类的点数Ni,就会有不同的判别边界,从而影响预测。

    如预测下图三个*点是正类<+>还是负类<o>:

    那么什么样的标签是好的呢?<对预测数据,得到的判别“合理”>

    从判别边界来看,应该让判别边界“居中”,那么新来的预测点在靠近判别边界的地方,也能得到“合理”的判断。

    由此可以看到:SVM最大边际的思想!!!

    当时,在这个地方思考了两三天,推导了部分边际的内容。觉得这个问题应该有人做过,就继续往下看了,后来看到SVM,觉得不谋而合了。

    付代码:

    clear, close all, clc
    x0 = [1 4]';
    x1 = [5 1]';
    randn('seed',1)
    N0 = 20; N1 = 40;
    X0 = repmat(x0,1,N0) + rand(2,N0);
    X1 = repmat(x1,1,N1) + rand(2,N1);
    X = [X0 X1]';
    X = [ones(N0+N1,1) X];
    
    plot(X0(1,:),X0(2,:), 'ko', 'MarkerFaceColor', 'y', 'MarkerSize', 7);
    hold on
    plot(X1(1,:),X1(2,:), 'k+','LineWidth', 2, 'MarkerSize', 7);
    axis([0 6 -3 8])
    
    %% 1
    T0 = zeros(N0,1);T1 = ones(N1,1);T = [T0;T1];
    w = inv(X'*X)*X'*T;
    k = -w(2)/w(3); b = -w(1)/w(3);
    h = refline(k,b); %已知斜率w 截距b 画直线
    set(h, 'LineWidth', 1.5),text(1,k+b-0.3, 't=\{0,+1\}')
    axis([0 6 -3 8])
    hold on
    
    %% 2
    T0 = -ones(N0,1);T1 = ones(N1,1);T = [T0;T1];
    w = inv(X'*X)*X'*T;
    k = -w(2)/w(3); b = -w(1)/w(3);
    h = refline(k,b); %已知斜率w 截距b 画直线
    set(h, 'LineWidth', 1.5, 'Color', 'r'),text(2,k*2+b, 't=\{-1,+1\}')
    axis([0 6 -3 8])
    
    %% 3
    T0 = (N0+N1)/N0*ones(N0,1);T1 = -(N0+N1)/N1*ones(N1,1);T = [T0;T1];
    w = inv(X'*X)*X'*T;
    k = -w(2)/w(3); b = -w(1)/w(3);
    h = refline(k,b); %已知斜率w 截距b 画直线
    set(h, 'LineWidth', 1.5, 'Color', 'g'),text(3,k*3+b, 't=\{+N/N0,-N/N1\}')
    axis([0 6 -3 8])
    
    hold on 
    plot(2, 2, 'r*',3,2, 'b*', 4, 1, 'c*')
    

    学习SVM的一个好的中文资料 李航著的《统计学习方法》第七章;讲得很好、很到位;列了不少参考资料。

    +维基





  • 相关阅读:
    十大经典算法总结
    十大经典算法总结
    MySQL主从同步模拟
    MySQL主从同步模拟
    高斯定理
    高斯定理
    如何修改数据决策系统登陆地址为ip
    如何修改数据决策系统登陆地址为ip
    数据库连接池问题 Max Pool Size
    数据库连接池问题 Max Pool Size
  • 原文地址:https://www.cnblogs.com/xinyuyuanm/p/3073061.html
Copyright © 2011-2022 走看看