在C ++或MATLAB上实现SVM-SMO算法

时间:2016-11-07 16:12:22

标签: algorithm matlab svm

有人可以帮我找到没有以下代码汇总在MATLAB上的错误吗?请注意,k表示内核并具有自己的功能。

我尝试使用图片中的算法实现它:

enter image description here

我的实施:

while 1

    num_changed_alphas = 0;
    %step 4
     if (Y .* alphas < B)   
    [Yigi , i_WS] = max(Y .* g)
     end
     %step 5
     if (A < Y .* alphas)   
    [Yjgj , j_WS] = max(Y .* g)
     end
     %step 6
     if ( (Yigi - Yjgj ) <= eps )
         num_changed_alphas = num_changed_alphas + 1;
         break;

     end
     %step 7
     lambda = min (B(i_WS) - Y(i_WS) .* alphas(i_WS), Y(j_WS) .* alphas(j_WS) - A(j_WS) );
     lambda = min (lambda, (Yigi - Yjgj) ./ (K(i_WS,i_WS) + K(j_WS,j_WS) - 2* K(j_WS,i_WS)) );

     g = sum (g - lambda .* Y .* K (:,i_WS) + lambda .* Y .* K (:,j_WS));

     alphas = alphas + Y * lambda;

    fprintf('.');
    dots = dots + 1;
    if dots > 75
        dots = 0;
        fprintf('\n');
    end
    if exist('OCTAVE_VERSION')
        fflush(stdout);
    end
end %while

1 个答案:

答案 0 :(得分:1)

一个很大的错误就是

  

a = f(x)受g(x)

的影响

并不意味着&#34;如果g(x)则a = f(x)&#34;,而是&#34; a = f(x)只考虑这样的x&g; x)是真的&#34;。

所以

if (Y .* alphas < B)   
    [Yigi , i_WS] = max(Y .* g)
end

(以及下一个if)不是有效的操作实现

i = arg max SOMETHING   subject to CONSTRAINTS

你应该在

的行中做点什么
selected = (Y .* alphas < B) 
[Yigi , i_WS] = max(Y(selected) .* g(selected))