灰色关联度matlab源程序word版本.doc
《灰色关联度matlab源程序word版本.doc》由会员分享,可在线阅读,更多相关《灰色关联度matlab源程序word版本.doc(21页珍藏版)》请在咨信网上搜索。
1、灰色关联度matlab源程序精品文档灰色关联度matlab源程序 最近几天一直在写算法,其实网上可以下到这些算法的源程序的,但是为了搞懂,搞清楚,还是自己一个一个的看了,写了,作为自身的积累,而且自己的的矩阵计算类库也迅速得到补充,以后关于算法方面,基本的矩阵运算不用再重复写了,挺好的,是种积累,下面把灰关联的matlab程序与大家分享。灰色关联度分析法是将研究对象及影响因素的因子值视为一条线上的点,与待识别对象及影响因素的因子值所绘制的曲线进行比较,比较它们之间的贴近度,并分别量化,计算出研究对象与待识别对象各影响因素之间的贴近程度的关联度,通过比较各关联度的大小来判断待识别对象对研究对象的
2、影响程度。关联度计算的预处理,一般初值化或者均值化,根据我的实际需要,本程序中使用的是比较序列与参考序列组成的矩阵除以参考序列的列均值等到的,当然也可以是其他方法。%注意:由于需要,均值化方法采用各组值除以样本的各列平均值clear;clc;yangben=47.924375 25.168125 827.4105438 330.08875 1045.164375 261.37437516.3372 6.62 940.2824 709.2752 962.1284 84.87455.69666667 30.80333333 885.21 275.8066667 1052.42 435.81; %样
3、本数据fangzhen=36.27 14.59 836.15 420.41 1011.83 189.54 64.73 35.63 755.45 331.32 978.5 257.87 42.44 23.07 846 348.05 1025.4 296.69 59.34 39.7 794.31 334.63 1016.4 317.27 52.91 17.14 821.79 306.92 1141.94 122.04 4.21 4.86 1815.52 2584.68 963.61 0.00 6.01 2.43 1791.61 2338.17 1278.08 30.87 3.01 1.58 122
4、0.54 956.14 1244.75 3.91 25.65 7.42 790.17 328.88 1026.01 92.82 115.80 27 926.5 350.93 1079.49 544.38 12.63 8.75 1055.50 1379.00 875.10 1.65 ; %待判数据rows,cols=size(fangzhen);p=0.5; %分辨系数m,n=size(yangben);R=;for irow=1:rows yy=fangzhen(irow,:); data=yy;yangben; data_gyh1=mean(yangben) for i=1:m+1 for
5、j=1:n data_gyh(i,j)=data(i,j)/data_gyh1(j); end end for i=2:m+1 for j=1:n Dij(i-1,j)=abs(data_gyh(1,j)-data_gyh(i,j); end end Dijmax=max(max(Dij); Dijmin=min(min(Dij); for i=1:m for j=1:n Lij(i,j)=(Dijmin+p*Dijmax)/(Dij(i,j)+p*Dijmax); end end LijRowSum=sum(Lij); for i=1:m Rij(i)=LijRowSum(i)/n; end
6、 R=R;Rij;endRmatlab求灰色关联度矩阵源代码 2010-12-11 22:57function greyrelationaldegree(X,c)%GRAYRELATIONALDEGREE this function is used for calculating the gery%relation between squence%rememeber that the first column of the input matrix is the desicion%attribution squences.what we want to calculate is the gre
7、y ralational degree between%it and other attributions%X is the squence matrix, c is the parameter used in the function%in most of the time, the value of c is 0.5firstrow = X(1,:);reci_firstrow = 1./firstrow;reci_convert = diag(reci_firstrow);initialMIRROR = X*reci_convert;% find the initial value mi
8、rror of the sequce matrixA = initialMIRRORnrow,ncolumn = size(A);for (i=2:nrow) C = A(i,:)-A(1,:) D=abs(C); eval(B num2str(i) =D); amax = max(eval(B num2str(i) amin = min(eval(B num2str(i) maxarray(i-1)=amax minarray(i-1)=aminend %find the difference squence and the max value and min value of each s
9、quencemaxmax = max(maxarray)minmin = min(minarray)for(i=2:nrow) for(j=1:ncolumn) eval(greyrelationdegree num2str(i) (j)=(minmin+c*maxmax)/(B num2str(i) (j)+c*maxmax) endend % calculate the greyralational degree of each datafor(i=2:nrow) eval(greyrelatioanaldegree_value num2str(i) = mean (greyrelatio
10、ndegree num2str(i) ) )end基于matlab灰色关联度计算的实现 2006年07月28日 星期五 上午 11:06function r=incident_degree(x0,x1) %compute the incident degree for grey model. %Designed by NIXIUHUI,Dalian Fisher University. %17 August,2004,Last modified by NXH at 21 August,2004 %数据初值化处理 x0_initial=x0./x0(1); temp=size(x1); b=re
11、pmat(x1(:,1),1 temp(2); x1_initial=x1./b; %分辨系数选择 K=0.1; disp(The grey interconnect degree is: ); x0_ext=repmat(x0_initial,temp(1) 1); contrast_mat=abs(x0_ext-x1_initial); delta_min=min(min(contrast_mat);%delta_min在数据初值化后实际为零 delta_max=max(max(contrast_mat); a=delta_min+K*delta_max; incidence_coeffi
12、cient=a./(contrast_mat+K*delta_max);%得到关联系数 r=(sum(incidence_coefficient)/temp(2); %得到邓氏面积关联度我们根据图1的步骤和图2的数据进行编程实现,程序如下:%清除内存空间等clear;close all;clc;%载入源数据 %其实这里可以载入execl表格的n=15; %参与评价的人数m=4; %参与评价的指标个数X_0=zeros(n,m); % 数据矩阵X_2=zeros(n,m); %偏差结果的求取矩阵X_3=zeros(n,m); % 相关系数计算矩阵a1_0=13 18 17 18 17 17 18
13、 17 13 17 18 13 18 13 18;a2_0=18 18 17 17 18 13 17 13 18 13 17 13 13 17 17;a3_0=48.67 43.33 43.56 41.89 39.47 43.44 37.97 41.14 39.67 39.83 34.11 40.58 34.19 30.75 21.22;a4_0=10 10.7 3 5.4 5.4 0.7 4.2 0.5 9.3 0.85 2.9 5.45 4.2 2.7 6; %指标数X_1=a1_0,a2_0,a3_0,a4_0; %最后使用到的数据矩阵%1 寻找参考列x0=max(a1_0),max(a
14、2_0),max(a3_0),max(a4_0); %取每列的最大值(指标的最大值)%2 计算偏差结果i=1;while(i=m+1) %为什么这个地方会出问题呢 for j=1:1:n X_2(j,i)=abs(X_1(j,i)-x0(i); end; i=i+1;end%3 确定偏差的最值error_min=min(min(X_2);error_max=max(max(X_2);%4 计算相关系数i=1;p=0.5;while(i=m+1) for j=1:1:n X_3(j,i)=(error_min+p*error_max)/(X_2(j,i)+p*error_max); end; i
15、=i+1;end%X_3 %可以在此观察关联矩阵%5 计算各个学生的关连序a=zeros(1,n);for j=1:1:n for i=1:1:m a(j)=a(j)+X_3(j,i); %其实可以直接用sum end; a(j)=a(j)/m; %可以改进% 1 2 3 下一页%end%a %在此可以观测各个学生的序 %改进:如果各个指标的所占权重不一样的话,可以添加相应的权系数%6 排序b=a;c,s=sort(b);for i=1:1:n d(i)=i;endd=d;result=d b c s%7 将结果显示出来figure(1);plot(a);figure(2)bar(a); %柱
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 灰色 关联 matlab 源程序 word 版本
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。