plant.m

用MATLAB编写的支持向量机

x1=[1.4 0.2;1.4 0.2;1.3 0.2;1.5 0.2;1.4 0.2;...
    1.7 0.4;1.4 0.3;1.5 0.2;1.4 0.2;1.5 0.1;...
    1.5 0.2;1.6 0.2;1.4 0.1;1.1 0.1;1.2 0.2;...
    1.5 0.4;1.3 0.4;1.4 0.3;1.7 0.3;1.5 0.3;...
    1.7 0.2;1.5 0.4;1.0 0.2;1.7 0.5;1.9 0.2;...
    1.6 0.2;1.6 0.4;1.5 0.2;1.4 0.2;1.6 0.2;...
    1.6 0.2;1.5 0.4;1.5 0.1;1.4 0.2;1.5 0.1;...
    1.2 0.2;1.3 0.2;1.5 0.1;1.3 0.2;1.5 0.2;...
    1.3 0.3;1.3 0.3;1.3 0.2;1.6 0.6;1.9 0.4;...
    1.4 0.3;1.6 0.2;1.4 0.2;1.5 0.2;1.4 0.2];
x2=[4.7 1.4;4.5 1.5;4.9 1.5;4.0 1.3;4.6 1.5;...
    4.5 1.3;4.7 1.6;3.3 1.0;4.6 1.3;3.9 1.4;...
    3.5 1.0;4.2 1.5;4.0 1.0;4.7 1.4;3.6 1.3;...
    4.4 1.4;4.5 1.5;4.1 1.0;4.5 1.5;3.9 1.1;...
    4.8 1.8;4.0 1.3;4.9 1.5;4.7 1.2;4.3 1.3;...
    4.4 1.4;4.8 1.4;5.0 1.7;4.5 1.5;3.5 1.0;...
    3.8 1.0;3.7 1.0;3.9 1.2;5.1 1.6;4.5 1.5;...
    4.5 1.6;4.7 1.5;4.4 1.3;4.1 1.3;4.0 1.3;...
    4.4 1.2;4.6 1.4;4.0 1.2;3.3 1.0;4.2 1.3;...
    4.2 1.2;4.2 1.3;4.3 1.3;3.0 1.1;4.1 1.3];
x3=[6.0 2.5;5.1 1.9;5.9 2.1;5.6 1.8;5.8 2.2;...
    6.6 2.1;4.5 1.7;6.3 1.8;5.8 1.8;6.1 2.5;...
    5.1 2.0;5.3 1.9;5.5 2.1;5.0 2.0;5.1 2.4;...
    5.3 2.3;5.5 1.8;6.7 2.2;6.9 2.3;5.0 1.5;...
    5.7 2.3;4.9 2.0;6.7 2.0;4.9 1.8;5.7 2.1;...
    6.0 1.8;4.8 1.8;4.9 1.8;5.6 2.1;5.8 1.6;...
    6.1 1.9;6.4 2.0;5.6 2.2;5.1 1.5;5.6 1.4;...
    6.3 2.3;5.6 2.4;5.5 1.8;4.8 1.8;5.4 2.1;...
    5.6 2.4;5.1 2.3;5.1 1.9;5.9 2.3;5.7 2.5;...
    5.2 2.3;5.0 1.9;5.2 2.0;5.4 2.3;5.1 1.8];
N1=size(x1,1);
N2=size(x2,1);
N3=size(x3,1);
y1=ones(1,N1);
y2=-ones(1,N2);
y3=-ones(1,N3);

figure;
plot(x1(:,1),x1(:,2),'o',x2(:,1),x2(:,2),'+',x3(:,1),x3(:,2),'d');
axis([0 8 -1 4]);
hold on;

X=[x1;x2;x3];
y=[y1,y2,y3];
x=X;
C=Inf;
[g,z,a,k,d]=matlabsvm1(C,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),'rp') 
    end
end

pp=-3:0.1:8;
A=sum(a.*(y').*X(:,2));
B=sum(a.*(y').*(X*(X(k(1),:))'))-y(k(1))-sum(a.*(y').*X(:,1))*pp;
E=B+1;
F=B-1;
t1=B/A;
t2=E/A;
t3=F/A;
plot(pp,t1,'r',pp,t2,'--',pp,t3,'--')