matlab代做highspeedlogic★基于matlab编程的最小距离分类判别

a=[2 2 2 3;0 2 4 3];   
b=[0 -2 -1 1 3;3 2 -1 -2 -1];
m1=mean(a,2);
m2=mean(b,2);
p=(m1+m2)/2;
x=0.2:2.25;
y=x;%两类中点连线方程.
plot(x,y);
hold on;
plot(m1(1),m1(2),'h');
plot(m2(1),m2(2),'h');
x=-4:6;
y=-x+2.45;%分类函数;y=kx+c,两类中点连线的斜率为1,故k=-1(垂直),直线过点p(1.225 1.225),易求c=2.45.
plot(x,y,'c');
for i=1:4
plot(a(1,i),a(2,i),'ro')
end
for i=1:5
plot(b(1,i),b(2,i),'bo')
end
grid on;
axis equal

w1=[[2,0];[2,2];[2,4];[3,3]];
w1(:,1);
w2=[[0,3];[-2,2];[-1,-1];[1,-2];[3,-1]];
n1=length(w1);
n2=length(w2);
w1_t=0;
for i=1:n1
 w1_t=w1_t+w1(i,:);
end
w1_e=w1_t./n1;
w2_t=0;
for i=1:n1
 w2_t=w2_t+w2(i,:);
end
w2_e=w2_t./n1;
w12mid=(w1_e+w2_e)./2;
w12mid(1);
%slop=-(w2_e(2)-w1_e(2))./(w2_e(1)-w1_e(1));
%t=10;
%a=[w12mid(1)-t,w12mid(1),w12mid(1)+t]';
%b=[w12mid(2)-t.*slop,w12mid(1),w12mid(1)+t.*slop]';
%figure(3)
%line(a,b);
%hold on
plot(w1(:,1)',w1(:,2)','ro');
hold on
plot(w1_e(:,1),w1_e(:,2),'k*');
hold on
plot(w2(:,1)',w2(:,2)','gx');
hold on
plot(w2_e(:,1),w2_e(:,2),'k*');
hold on
slop=-(w2_e(2)-w1_e(2))./(w2_e(1)-w1_e(1));
t=10;
a=[w12mid(1)-t,w12mid(1),w12mid(1)+t]';
b=[w12mid(2)-t.*slop,w12mid(1),w12mid(1)+t.*slop]';
line(a,b);
hold on
hold off