模式识别作业第三章1.docVIP

第三章作业 1.设有10个两维模式, ;若选择阈值的系数为,使用最大最小算法进行聚类。 解：先定义几个函数： 文件名为caldistance.m函数 %求其他样本点与选定点的第m个的最大距离，并找到这个最大距离点 Function [distance,i]=caldistance(x,m) [M,N]=size(x); distance=0; index=0; for i=1:N index=index+1; temp=sum((x(:,m)-x(:,i)).^2); distance(index)=sqrt(temp); v(:,i)=distance(index); end end 文件名为jiaoxiaojuli.m函数 %其他样本分别与第m、n个样本之间的欧氏距离距离，取较小的 function [vv i]=jiaoxiaojuli(x,m,n) [M,N]=size(x); ddistance=0; for i=1:N temp=sum((x(:,m)-x(:,i)).^2); ddistance(1)=sqrt(temp); temp1=sum((x(:,n)-x(:,i)).^2); ddistance(2)=sqrt(temp1); ddistance=min(ddistance(1),ddistance(2)); vv(:,i)=ddistance; end end 文件名为fenlei.m函数 %其他各样本点被分到离z1、z2近的类中 function [w1 w]=fenlei(x,m,n) [M,N]=size(x); for i=1:N a1=sqrt(sum((x(:,m)-x(:,i)).^2)); a2=sqrt(sum((x(:,n)-x(:,i)).^2)); if a1a2 w1(1,i)=i; w(1,i)=0; else w1(1,i)=0; w(1,i)=i; end end end 文件名为z3jiaoxiao.m函数 %第m、n和q个样本与其他样本之间的欧氏距离取较小的一个 function [s]=z3jiaoxiao(x,m,n,q) [M,N]=size(x); ddistance=0; for i=1:N temp=sum((x(:,m)-x(:,i)).^2); ddistance(1)=sqrt(temp); temp=sum((x(:,n)-x(:,i)).^2); ddistance(2)=sqrt(temp); temp=sum((x(:,q)-x(:,i)).^2); ddistance(3)=sqrt(temp); ddistance=min(min(ddistance(1),ddistance(2)),ddistance(3)); s(:,i)=ddistance; end end 文件名为z3fenlei.m函数 %将x中的各样本点分类到离第m、n、p较近的类中 function [w1 w2 w3]=z3fenlei(x,m,n,p) [M,N]=size(x); for i=1:N a1=sqrt(sum((x(:,m)-x(:,i)).^2)); a2=sqrt(sum((x(:,n)-x(:,i)).^2)); a3=sqrt(sum((x(:,p)-x(:,i)).^2)); if a1=a2a1=a3 w1(1,i)=i; w2(1,i)=0; w3(1,i)=0; elseif a2=a1a2=a3 w1(1,i)=0; w2(1,i)=i; w3(1,i)=0; elseif a3a1a3a2 w1(1,i)=0; w2(1,i)=0; w3(1,i)=i; end end end 总程序为： x1=[0 0];x2=[1 1];x3=[2 2];x4=[3 7];x5=[3 6];x6=[4 6];x7=[5 7];x8=[6 3];x9=[7 3];x10=[7 4]; X=[x1 x2 x3 x4 x5 x6 x7 x8 x9 x10]; z1=x2; [distance,i]=