模式作业设计.doc

1.1题目给出K-均值算法的程序框图，编写程序，自选一组分别属于三类的三维模式样本（图2.13），并给它们进行聚类分析。1.2算法原理K-均值算法也称C-均值算法，是根据函数准则进行分类的聚类算法，基于使聚类准则最小化。此处所用的聚类准则函数是聚类集中每一个样本点到该类聚类中心的距离平方和，对于第j个聚类集，准则函数定义为式中，Sj表示第j个聚类集，也称聚类域，其聚类中心为；为第j个聚类集中所包含的样本个数。（1）任选K个初始聚类中心Z1(1)，Z2(1)，ZK(1)，。括号内的序号代表寻找聚类中心的迭代运算的次序号。一般可选择样本集中前K个样本作为初始聚类中心。（2）按最小距离原则将其余样本分配到K个聚类中心中的某一个，即：若，则。式中，k代表迭代运算次序号；K代表聚类中心的个数。（3）计算各个聚类中心的新向量值，即以均值向量作为新的聚类中心。这一步要分别计算K个聚类中心的样本均值向量，故该算法被称为K-均值算法。（4）如果，则回到（2），将模式样本逐个重新分类，并重复迭代计算；如果，算法收敛，计算完毕。1.3程序流程图1.4 MATLAB程序代码clear all;clc;data=input(请输入样本数据矩阵：);X=data(:,1);Y=data(:,2);figure(1);plot(X,Y,r*,LineWidth,3);axis(0 9 0 8)xlabel(x);ylabel(y);hold on;grid on;m=size(data,1);n=size(data,2);counter=0;k=input(请输入聚类数目：);if km disp(输入的聚类数目过大，请输入正确的k值); k=input(请输入聚类数目：);endM=cell(1,m);for i=1:k M1,i=zeros(1,n);endMold=cell(1,m);for i=1:k Mold1,i=zeros(1,n);end%随机选取k个样本作为初始聚类中心%第一次聚类,使用初始聚类中心 p=randperm(m);%产生m个不同的随机数for i=1:k M1,i=data(p(i),:);endwhile true counter=counter+1; disp(第); disp(counter); disp(次迭代);count=zeros(1,k);%初始化聚类CC=cell(1,k);for i=1:k C1,i=zeros(m,n);end%聚类for i=1:m gap=zeros(1,k); for d=1:k for j=1:n gap(d)=gap(d)+(M1,d(j)-data(i,j)2; endend y,l=min(sqrt(gap); count(l)=count(l)+1; C1,l(count(l),:)=data(i,:);end Mold=M; disp(聚类中心为：); for i=1:k disp(M1,i); end disp(聚类结果为：); for i=1:k disp(C1,i); end sumvar=0; for i=1:k E=0; disp(单个误差平方和为：); for j=1:count(i) for h=1:n E=E+(M1,i(h)-C1,i(j,h)2; end end disp(E); sumvar=sumvar+E; end disp(总体误差平方和为：); disp(sumvar); %计算新的聚类中心,更新M,并保存旧的聚类中心 for i=1:k M1,i=sum(C1,i)/count(i); end %检查前后两次聚类中心是否变化，若变化则继续迭代；否则算法停止； tally=0; for i=1:k if abs(Mold1,i-M1,i)0; CountError=0; for i=1:ro Temp=W*ExpandData(i,:); if TempP2) disp(属于第一类); else disp(属于第二类); end;end3.5结果与分析程序没有完整的考虑两类问题中的每个情况，只解答了两类样本的协方差矩阵相等，即的情况。而且绘制判别函数的语句可读性差。4.1题目自己选择一幅彩色图像，选择一种分类算法，使用颜色特征对彩色图像中的像素进行分类。最终输出彩色图像按照颜色的分割结果。4.4程序代码clear all;clc;file_name=input(输入文件图像路径);I_rgb=imread(file_name);figure();image(I_rgb);title(原始图像);whos file_name;m,n,k=size(file_name);new=double(reshape(file_name,m*n,k);whos new;startdata=input(背景和前背景初始凝聚点);idClass=kmeans(new,2,Start,startdata);idfile=(idClass=1);result=