myprsvm.m

模式识别,对iris数据进行分类,svm算法

function SVMlabel=MyPRSVM(SV,Alpha,Thsvm,Testdata,Kerneltype,Sigma)
dim=size(SV,2);
Y=SV(:,1);
SV=SV(:,2:dim);
switch Kerneltype
    case 'linear'
        K=SV*Testdata';
    case 'poly'
        K=(SV*Testdata'+1).^2;
    case 'gauss'
        n=length(Y);
        tmp=zeros(n,1);
        for i=1:n
            t=SV(i,:)-Testdata;
            tmp(i,1)=-0.5*(t*t');
        end        
        K=exp(tmp./(Sigma.^2));
    otherwise
        K=0;
end
SVMlabel=Alpha.*Y.*K;
SVMlabel=sum(SVMlabel)+Thsvm;