shiyan33.m

支持向量数据描述.里面有对应的图,一看就明白

%3.3图
clc
clear all
load('E:\智能信息处理\程序\data\banana\shiyan1.mat');
 plot(X(:,1),X(:,2),'b.',Xe1(:,1),Xe1(:,2),'r*',Xe2(:,1),Xe2(:,2),'kx');
 legend('训练样本','测试目标样本','测试非目标样本');
 hold on
n = size(X,1); 
 epsilon =1e-8;
 C=1;
 p1=0.3;
  H = zeros(n,n); 
    for i=1:n
       for j=1:n
            H(i,j)=exp(-(X(i,:)-X(j,:))*(X(i,:)-X(j,:))'/(2*p1^2));
       end;
    end;
    f = -ones(n,1); 
    H = H+1e-10*eye(size(H));
    lb = zeros(n,1);    
    ub = C*ones(n,1);   
    Aeq = ones(1,n);, beq = 1;    
    st = cputime;  
    [alpha,FV,EF] = quadprog(2*H, f,[],[], Aeq, beq, lb, ub);
    svi = find( alpha > epsilon);   
    nsv = length(svi);
    R=0;
    for k=1:nsv
     for i=1:n
      R=R+2*alpha(i)*exp(-(X(svi(k),:)-X(i,:))*(X(svi(k),:)-X(i,:))'/(2*p1^2));
     end;
   end;
   R=1+alpha'*H*alpha-R/nsv;
 
n=size(X,1);
x=sym('x');
y=sym('y'); 
p1=0.3;
f=0;
for i=1:n
  f=f+2*alpha(i)*exp(-([x,y]-X(i,:))*([x,y]-X(i,:))'/(2*p1^2));
end;    
f=1+alpha'*H*alpha-f-R;
ezplot(f);