fenleixian.asv

matlab支持向量机程序SVM_SteveGu

function [w,b]=fenleixian(X,Y)
%  求分类线
%  X      - Training inputs
%  Y      - Training targets

 for i=1:n
       for j=1:n
 ker(x(i),x(j))=X(i,:)*(X(j,:))';   
 [alpha,b] = trainlssvm({X,Yc,type,gam,sig2,kernel_type,preprocess});  % 核函数为线性核函数
 w=alpha'*Y(i)*Y(j)*ker(x(i),x(j)*alpha;      
       end
 end
 epsilon = svtol(C);
% 求b 
svi = find( alpha > epsilon); % 得到 svi为支持向量;
svii=find(alpha>epsilon & alpha<(C-epsilon));
 if   length(svii)>0
   b=(1/length(svii))*sum(Y(svii)-Y(svii)*Y(svi)*ker(x(svii),x(svi)）*alpha(svi).*Y(svii));
 end