原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
模糊数学期末论文
模糊数学期末论文
针对北京、天津、河北、山西、内蒙古、辽宁、吉林、黑龙江，我国北方的8个省或直辖市省级行政区农村人口的食品消费数据进行模糊聚类分析。
第一步：数据标准化
数据矩阵
设论域表示我国这8个省市的样本数据，每个样本用6个指标表示其性状，即
原始数据如下【1】：
样本编号 行政区 粮食 蔬菜 猪牛羊肉 家禽 蛋类 水产品 1 北京 99.37 96.25 18.2 3.29 10.88 5.65 2 天津 140.92 62.87 12.6 1.23 11.6 9.62 3 河北 180.15 53.66 7.21 0.89 7.22 2.61 4 山西 177.65 70.56 5.71 0.69 6.05 0.8 5 内蒙古 187.04 70.9 27.96 2.79 5.36 2.06 6 辽宁 179.36 167.92 18.52 2.36 2.06 1.01 7 吉林 153.26 99.25 11.32 3.55 8.56 4.04 8 黑龙江 156.23 83.67 10.1 3.41 6.22 4.29 ??数据标准化
令,消除量纲的影响，利用matlab程序（见附件一）得到样本标准化矩阵：
其中，
第二步：标定——建立模糊相似矩阵
采用欧式距离法建立相似矩阵，令
为了使，选取（在这里，我首先计算了各样本间的欧式距离，发现最大的距离为2.7893，所以选择了），利用matlab程序（见附件二）得模糊相似矩阵如下：
第三步：聚类
用平方法求传递闭包,经过4次平方法得到传递闭包矩阵，即模糊矩阵。（matlab源程序见附件三）：
取，分为8类
取，分为7类，
取，分为6类，
取，分为5类，
取，分为4类，
取，分为3类，
取，分为2类，
取，分为1类，
针对以上8个省级行政区，从实际意义上说，将其分为3-4种情况较好，即和较好。
以上聚类对于我国北方地区调整居民食品供给具有十分重要的意义，能给政府以及市场的经济调整提供较好的作用。
附件一：
function biaozhunhua(A)
[m,n]=size(A);
B=zeros(1,n);
S=zeros(1,n);
E=zeros(m,n);
for i=1:n
sum=0;
for j=1:m
sum=sum+A(j,i);
end
B(i)=sum/m;
var1=0;
for j=1:m
var1=var1+(A(j,i)-B(i))^2;
end
S(i)=(var1/m)^1/2;
end
for i=1:m
for j=1:n
E(i,j)=(A(i,j)-B(j))/S(j);
end
end
E
附件二
function xiangsijuzhen(A)
[m,n]=size(A);
r=zeros(m);
for i=1:m
for j=1:m
s=0;
for k=1:n
s=(A(i,k)-A(j,k))^2;
end
r(i,j)=1-1/3*s^1/2;
r(j,i)=1-1/3*s^1/2;
end
end
r
附件三
function chuandi(R)
a=size(R);
B=zeros(a);
flag=0;
l=0;
while flag==0
for i= 1: a
for j= 1: a
for k=1:a
B( i , j ) = max(min( R( i , k) , R( k, j) ) , B( i , j ) ) ;%R与R内积，先取小再取大
end
end
end
if B==R
flag=1;
else
R=B;%循环计算R传递闭包
end
l=l+1;
end
R
l
文档评论（0）