20100373苏浩辉实验2.docVIP

重庆大学学生上机实验报告 上机课程名称：信息智能分析与预测 上机实验名称：实验2模糊数学与模糊聚类 开课学院及实验室：经济与工商管理学院 2013年 11 月 26 日 学 院 经管学院 年级、专业、班 2010级信管1班 姓 名 苏浩辉 学 号 指导教师 陈铿 成 绩 教师评语 教师签名： 【实验目的】 借助MATLAB工具，试对重庆各区县进行模糊聚类。 【实验内容】 1.按照模糊聚类的基本计算步骤处理重庆40区县的数据，并展示聚类结果 按照步骤编写代码如下： （1）主程序代码 C = xlsread(D:\CQDATA.xls,1,B5:V44);%打开原分析数据 C2 = normalization(C); %数据标准化函数-平移极差变换 C3 = set_similar_matrix(C2); %建立相似矩阵-相关系数法 C4 = pf_package(C3); %创建模糊等价矩阵-平方法求传递包 Alamd = fuzzy_lamd(C4,7); （2）输出结果 [M,N]= f_ dynamic_clustering(C4); （3）各函数代码 function y = normalization(x)%数据标准化函数-平移极差变换 [m, n] = size(x); y = zeros(m, n); for i = 1: m for j = 1: n ma = max(x(:, j)); mi = min(x(:, j)); y(i,j) = (x(i,j)-mi)/(ma - mi); end end function R = set_similar_matrix(X)%建立相似矩阵-相关系数法 [m,n] = size(X); R = [m,m]; for i=1:m for j=1:m xi = mean(X(i,:)); xj = mean(X(j,:)); xis=0;xjs=0; for k=1:n xis=xis+(X(i,k)-xi)^2;xjs=xjs+(X(j,k)-xj)^2; end s=sqrt(xis*xjs);R(i,j)=0; for k=1:n R(i,j)=R(i,j)+abs((X(i,k)-xi)*(X(j,k)-xj)); end R(i,j)=R(i,j)/s; end end function R2 = pf_ package(R)%创建模糊等价矩阵-平方法求传递包 R2 = Max_Min(R,R); k = 1; while 1 fprintf(k = %d \n, k); R2 = Max_Min(R,R); k = k + 1; if R2 == R break; else R = R2; end end function C = Max_Min(A,B) %模糊矩阵的合成，先取小后取大 [m,s1]=size(A); [s2,n]=size(B); C=[m, n]; if(s1~=s2) disp(A的列不等于B的列); else for i=1:m for j=1:n C(i,j)=0; for k=1:s1 if(A(i,k)B(k,j)) x=A(i,k); else x=B(k,j); end if(C(i,j)x) C(i,j)=x; end end end end end function Alamd = fuzzy_lamd(A,q)%求矩阵的lamd截距阵并输出聚类结果 [m,n]=size(A); p=m*n; C=A(1:p); D=sort(C); for i=1:p for j=i+1:p if(D(i)==D(j)) D(j)=0; end