matlab实现感知机算法--统计学习小灶
1 clear all; 2 clc; 3 %% 4 %算法 5 %输入:训练数据集T = {(x1,y1),(x2,y2),...,(xn,yn)};学习率η 6 %输出:w,b;感知机模型f(x) = sign(w*x+b) 7 %(1)选取初值w0,b0 8 %(2)在训练集中选取数据(xi,yi) 9 %(3)如果yi(w*xi+b)<=0 10 % w = w+η*yi*xi 11 % b = b+ηyi 12 %(4)转至(2),直至训练集中没有误分类点 13 %% 14 15 %初始化 16 X = [3 3 1;4 3 1;1 1 -1];%训练集 17 [sn,fn] = size(X); 18 y = X(:,fn);%标签 19 X(:,fn) = []; 20 fn = fn-1; 21 w = zeros(2,1); 22 b = 0; 23 yta = 1;%参数η 24 ok = 0; 25 %% 26 %循环体 27 while(ok==0)%%%不断的扫描直到所有的实例都满足无误分类点才结束 28 for i = 1:sn 29 if (y(i)*(w'*X(i,:)'+b)<=0) 30 31 w = w + yta*y(i)*X(i,:)' 32 b = b + yta*y(i) 33 else 34 if(i == sn) ok=1; 35 break 36 end 37 end 38 end 39 end 40 41 %% 42 %图示 43 figure(1) 44 plot(X(1,:),X(2,:),'r.') 45 hold on 46 plot(X(3,:),'*') %样本 47 xmin=min(X(:,1)); 48 xmax=max(X(:,1)); 49 % ymin=min(X(:,2)); 50 % ymax=max(X(:,2)); 51 xindex=xmin-1:(xmax-xmin)/100:xmax+1; 52 yindex=w(1)*xindex+w(2)*xindex+b; 53 plot(xindex,yindex)%分类面