beiyesi.m

来自「几个模式识别的作业程序」· M 代码 · 共 16 行

s1=[0 0;2 0;2 2; 0 2];
s2=[4 4;6 4; 6 6;4 6];
m1=mean(s1);
m2=mean(s2);
n1=zeros(2,2);n2=zeros(2,2);
for i=1:4
    n1=n1+s1(i,:)'*s1(i,:);
end
c1=n1/4-m1'*m1;
for i=1:4
    n2=n2+s2(i,:)'*s2(i,:);
end
c1=n2/4-m2'*m2; 
% 则判别函数：d1=m1*c1*x+m1*c1*m1';d2=m2*c2*x+m2*c2*m2';d=d1-d2;
plot(s1(:,1),s1(:,2),'ro',s2(:,1),s2(:,2),'bo');
xlabel('x1');ylabel('x2');grid;