感知机例题 matlab

D=[3 3 1;4 3 1;1 1 -1];
X=D(:,[1,2]); %取1，2列为X数据
Y=D(:,3);%取第三列为标签
plotData2(X,Y);%将样本点画出来 
axis([-5 5 -5 5]);%横纵坐标范围
hold on
[x1,y1] = ML(X,Y);
plot(x1,y1,'-r','Linewidth',3);

function [x1,y1] = ML(X,Y)
%第一步，初始化参数
x1=-5:0.01:5; %x1坐标轴
m=length(Y);%样本大小
w=[0;0]; %训练数据是二维的，故w也是二维的
b=0;
r=0.5;  %学习率
% t=0;%迭代次数
error=1;
while error>0 
    error=0;
    for i=1:m
          if (((w'*X(i,:)'+b)*Y(i))<=0)
              error=error+1;
               w=w+r*Y(i)*X(i,:)';
               b=b+r*Y(i);
               y1=(-w(1)*x1-b)/w(2);
               plot(x1,y1,'-b');
               pause(0.5);
          end 
     end
end
end

function plotData2(X,Y)
  figure;hold on;
  pos=find(Y==1);%找到Y矩阵中等于1的数存于变量pos正中
  n