感知机对偶形式

a=[3,3;4,3;1,1];
y=[1;1;-1];
scatter(a(:,1),a(:,2))
ok=0;
alpha=zeros(size(a,1), 1);  
bias=0;  
eta=1;
gram=a*a';
while ok==0
for i=1:size(a,1)
    temp=0;
    for j=1:size(a,1)
         temp=temp+alpha(j,1)*y(j,1)*gram(i,j);
    end
    if y(i,1)*(temp+bias)<=0
        alpha(i,1)=alpha(i,1)+eta;
        bias=bias+y(i,1)*eta;
        ok=0;
        break
    else
        if i==size(a,1)
            ok=1;
            break
        end
    end
end
end
w=a'*(alpha.*y);
b=bias;