版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、第五章聚类分析安庆师范学院 胡云峰习题5.3.下表给出我国历年职工人数(单位:万人),请用有序样品的fisher法聚类。年份全民所有制集体所有制195215802319541881121195624235541958453266219605044925196233031012196434651136196639391264196841701334197047921424197256101524197460071644197668601813197874512048198080192425解 第一步数据标准化后计算直径d程序:x=1580 23;1881 121;2423 554;4532 66
2、2;5044 925;3303 1012;3465 1136;. 3939 1264;4170 1334;4792 1424;5610 1524;6007 1644;6860 1813;. 7451 2048;8019 2425;stdr=std(x);n,m=size(x);x=x./stdr(ones(n,1),:); n p=size(x);d=zeros(n,n);for i=1:1:n; for j=1:1:n; if ij t=i:1:j; xgjunzhi=(1/(j-i+1)*sum(x(t,:); y=zeros(1,j-i+1); for s=i:1:j y(s)=(x(s
3、,:)-xgjunzhi)*(x(s,:)-xgjunzhi); end s=i:1:j; d(i,j)=sum(y); else d(i,j)=0; end endendd=d输出结果 矩阵太大,所以用excel处理了一下d=00000.0000000000.448980.24578002.06321.39810.6000003.92562.6511.18020.110004.50223.00911.42380.569530.40005.1793.43531.66480.825760.538310.0006.08234.0211.9761.0230.633430.127810.0007.03
4、114.65022.32551.23130.7550.263410.112750.008.33225.57622.90941.60451.05310.606190.338810.131220.010.3127.10344.01172.41261.77721.37930.923140.526640.315410.12.6968.99725.44223.51142.75482.35531.6691.04570.654960.256320.03671000016.29111.9987.86885.50384.56864.11933.10322.14681.47070.771220.308580.12
5、76200021.11716.12811.3218.42987.23166.64875.21163.83122.77931.68770.88810.460160.10709002822.16716.52812.97811.38610.5468.55966.6275.07163.45392.17481.34430.598320.199510我们只看下三角所有元素,其它元素理解为空第二步我们计算损失函数矩阵l程序:%设计一个把样品分为两类的程序,以及对应最后一类分割点d=d;l=zeros(n-1,n-1);alp=zeros(n-1,n-1);for m=2:n; s=zeros(1,m-1);
6、 for j=2:m s(1,j-1)=d(1,j-1)+d(j,m); end l(m-1,1)=min(s(1,1:m-1); for j=1:m-1 if l(m-1,1)=s(1,j); alp(m-1,1)=j+1; end endend%分为k类for k=3:n; for m=k:n s=zeros(1,m-k+1); for j=k:m; s(1,j-k+1)=l(j-2,k-2)+d(j,m); end l(m-1,k-1)=min(s(1,1:m-k+1); for j=1:m-k+1 if l(m-1,k-1)=s(1,j); alp(m-1,k-1)=j+k-1; en
7、d end endend输出结果 这里由于表太大,用excel处理一下l=0000.000000000.448980.00000000.559960.133550.0000001.01850.559960.133550.000001.27470.58040.153990.0430070.0001.4720.687770.261360.150380.0430070.001.68030.823370.396960.166440.0554640.0328970.002.05351.16620.711620.285210.166440.0554640.0328970.02.86161.77970.92
8、2770.496360.265840.154860.0554640.0328970.3.96041.93661.07970.653280.321920.203150.0921740.0554640.0328970.0124600005.95282.36211.47471.02020.593790.321920.203150.0921740.0554640.03290.0124560008.71882.94162.04371.18680.760370.429010.310240.199270.0921740.055460.0328970.01245600alp=20033044444550000000000466660000000004666670000000046668880000000466888890000004688100041111111000041121212120004111111113130014140101255在这里解释一下这两个矩阵表示分为k类,k从2到15;列表示样本数m,m从2到15接下来我们根据结果
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- GB/T 26668-2024电子电气产品材料声明
- 昆明第一中学2025届高三年级第二次联考 地理试卷(含答案)
- 初中生物实验苏教版全攻略
- 新版苏教五年级音乐教学案
- 探索人教版观沧海课件秘密
- 小学数学解方程教材
- 北师大版英语教材解读
- 初二数学解方程解题策略
- 北师大版校园教学设计点评
- 四年级下册苏教版语文习作题目解析分享
- 华星光电机电包述标文件中电四
- 大学学生宿舍维修管理细则
- “阿里巴巴”并购“饿了么”案例分析
- 罪犯违规违纪问题分析
- 混凝土抗冻性检测1快冻法
- 医院食堂经营方案
- DB34T 4397-2023 城市公共交通首末站设置规范
- 下肢深静脉血栓的预防课件
- 项目建设领导小组工作方案三篇
- 成品检验规范
- 擦窗机施工组织设计
评论
0/150
提交评论