untitled.m

用MATLAB编写的支持向量机

clear;
N=10;
n =N/2;
randn('state',6);
x1 = randn(2,n);
y1 = ones(1,n);
x2 = 5+randn(2,n);
y2 = -ones(1,n);

figure;
plot(x1(1,:),x1(2,:),'*',x2(1,:),x2(2,:),'k.');
axis([-3 8 -3 8]);
hold on;

X1 = [x1,x2];
Y1 = [y1,y2]; 
X=X1';
y=Y1;
x=X;
C=10;
global KM 
global o 
global w
o=2;
w=0.001;
KM='Kgauss'
[gg,z,a,k,d,m]=matlabsvm2(C,KM,o,w,X,y,x)

 
a0=a';
for i=1:size(k,2)
    if a0(k(i))~=C
       plot(X(k(i),1),X(k(i),2),'s')
    elseif a0(k(i))==C
       plot(X(k(i),1),X(k(i),2),'ro') 
    end
end