




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、系统工程作业本次作业在matlab中的代码如下: 求可达矩阵M:在matlab中输入代码如下:A=0 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 1 0 0 0 0 0 0 1 0 0;1 0 0 1 0 0 0 0 0 0 0 0 0;1 1 1 0 0 0 0 1 1 0 0 0 0;1 0 0 1 0 0 0 0 0 0 0 0 0;1 0 0 1 0 0 0 0 0 0 0 0 0;1 0 0 0 0 0 0 0 0
2、 0 0 1 0;1 0 0 0 1 1 0 0 0 1 0 0 0;1 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 0 1 1 0 1 1 1 0 0 0;I=eye(13);M=(A+I)13;M(M=0)=1;M 求综合影响矩阵T:在matlab中输入代码如下:A=0 0 0 0 0 0 0 0 0 0 0 0 0;1 0 0 0 0 0 0 0 0 0 0 0 0 ;1 0 0 0 0 0 0 0 0 0 0 0 0 ;1 0 0 0 0 0 0 0 0 0 0 0 0 ;1 0 0 1 0 0 0 0 0 0 1 0 0 ;1 0 0 1 0 0 0 0 0 0 0
3、0 0 ;1 1 1 0 0 0 0 1 1 0 0 0 0 ;1 0 0 1 0 0 0 0 0 0 0 0 0 ;1 0 0 1 0 0 0 0 0 0 0 0 0 ;1 0 0 0 0 0 0 0 0 0 0 1 0 ;1 0 0 0 1 1 0 0 0 1 0 0 0 ;1 0 0 0 0 0 0 0 0 0 0 0 0 ;1 0 0 0 1 1 0 1 1 1 0 0 0 ;d=0;s=0;for i=1:13for j=1:13s=s+A(i,j);if dsd=s;endends=0;endG=A/d;I=eye(13);T=G/(I-G)T 求综合影响矩阵中的行和、列和、中心度
4、、原因度:在matlab中代码如下:x=0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;for i=1:13for j=1:13x(i,1)=x(i,1)+T(i,j);endendx(行和)y=0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0;for i=1:13for j=1:13y(1,i)=y(1,i)+T(j,i);endendy(列和)z=0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;for i=1:13z(i,1)
5、=x(i,1)+y(1,i);endz(中心度)c=0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0;for i=1:13c(i,1)=x(i,1)-y(1,i);endc(原因度)食物链的例子代码基本相同,就是矩阵输入不同,再次则不以此写出,只需将A矩阵替换为 0 1 1 1 1 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 1 1;0 0 0 0 0 0 0 0 0 1 1 1; 0 0 0 0 0 0 0 0 0 0 1 1;0 0 0 0 0 1 1 1 1 0 0 0;0 0 0 0 0 0 1 0 1 0
6、0 0;0 0 0 0 0 0 0 1 1 0 0 0;0 0 0 0 0 0 0 0 0 1 0 0;0 0 0 0 0 0 0 0 0 0 1 1;0 0 0 0 0 0 0 0 0 0 0 1;0 0 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 0 0 0 0 0 0;即可。其他代码一致。具体题目及其答案如下第二题:一个系统的邻接矩阵为A,求A的可达矩阵,并对其进行分解。求得A的可达矩阵M如下:计算系统的可达性集合和先行集合,以及二者的共同集合,如下表所示:iL(si)F(si)L(si)F(si)111-13121,22,7231,33,7341,44-9,11,1
7、3451,4-6,10-125,11,135,1161,4,65,6,11,13671-4,7-97781,4,87,8,13891,4,97,9,139101,10,125,10,11,1310111,4-6,10-125,11,135,11121,125,10-1312131,4-6,8-131313由上表可知: T=7,13,由于L(s3)L(s7),所以s3、s7属于同一个区域中,即M不可分区。 s1是第一级要素。划去s1对应的行列,所得表如下:iL(si)F(si)L(si)F(si)222,72333,73444-9,11,13454-6,10-125,11,135,1164,65
8、,6,11,13672-4,7-97784,87,8,13894,97,9,1391010,125,10,11,1310114-6,10-125,11,135,1112125,10-1312134-6,8-131313由上表可知:s2、s3、s4、s12是第二级要素。同理可得:s6、s8、s9、s10是第三级要素;s5、s7、s11是第四级要素;s13是第五级要素。按照上面所划分的要素等级,对可达矩阵M进行重新排列,得到如下:第三题:对于上述问题研究的系统,利用DEMATEL方法判断各要素之间的综合影响关系,并求出各要素的原因度和中心度。计算规范化直接影响矩阵:继续算得综合影响矩阵各要素之间的
9、综合影响矩阵要 素s1s2s3s4s5s6s7s8s9s10s11s12s13s10000000000000s20.1667000000000000s30.1667000000000000s40.1667000000000000s50.2397000.17620.02860.02860000.02860.17140.00480s60.1944000.1667000000000s70.28700.16670.16670.0556000016670.16670000s80.1944000.1667000000000s90.1944000.1677000000000s100.194400000000
10、000.16670s110.2714000.05710.17140.17140000.17140.02860.02860s120.1667000000000000s130.3362000.11270.17140.171400.16670.16670.17140.02860,02860各系统要素的原因度以及中心度要素行和列和原因度中心度s102.5787-2.57872.5787s20.16670.166700.3334s30.16670.166700.3334s40.16670.9027-0.7361.0694s50.67790.37140.30651.0493s60.36110.3714-0
11、.01030.7325s71.009401.00941.0094s80.36110.33340.02770.6945s90.36210.33340.02870.6955s100.36110.3714-0.01030.7325s110.89990.22860.67131.1285s120.16670.2287-0.0620.3954s131.353701.35371.3537从上表的结果中可以得到:各要素在系统中的重要程度依次是:s1、s13、s11、s4、s5、s7、s6、s10、s9、s8、s12、s2、s3;系统中原因要素是:s1、s4、s6、s10、s12;系统中结果要素是:s5、s7、
12、s8、s9、s11、s13在该系统中,s2、s3的原因度最低,其中心度也最低,可以考虑删除该要素,达到减少要素的目的;第四题:求解问题二中的系统结构模型。观察矩阵,不难发现,s5与s11的相应的行列元素完全相同,可以把二者当成一个系统看待,从而可缩减相应的行和列,得到缩减矩阵如下:在矩阵中,先找出一级与二级之间的关系,再找出二级与三级之间的关系,直到把各级找完为止。从中知=1,说明节点s2与处于第一级的节点s1有关,即s2s1,以此类推,则有:第一级:s1第二级:s2s1,s3s1,s4s1,s12s1;第三级:s6s4,s8s4,s9s4,s10s12;第四级:s5s6,s5s10,s7s2
13、,s7s3,s7s8,s7s9;第五级:s13s5。根据上述画出结构矩阵如下:从而,可以绘制出系统的多级递阶有向结构图,如下图所示:第五题:探讨ISM方法与DEMATEL方法的区别与联系。区别:DEMATEL方法是由美国学者提出的一种运用图论和矩阵论原来进行系统因素分析的方法,它借助系统中各因素之间的逻辑关系构建直接影响矩阵,计算各因素对其他因素的影响程度以及被影响程度,从而计算各因素的中心度和原因度。根据因素所对应的中心度和原因度,得出该因素所属种类(原因因素还是结果因素),也可根据中心度和原因度的取值掉正系统的结构图,使系统的结构更加合理。ISM方法是现代系统工程学中广泛应用的一种分析和揭
14、示系统结构的方法。ISM作为一种分析系统结构的方法,将系统要素之间复杂、凌乱的关系分解成清晰的、多层级的结构形式。以ISM法得到的系统结构是一种宏观的定性结构,它揭示系统的几何学的定性结构,而不是对其结构做出精确的代数描述,或给出数量上、统计上的性质。也就是说。ISM法在分析系统结构时,是以系统的各要素为研究对象,以各要素间是否存在某种关系(从属关系、并列关系、因果关系等)为构造模型的唯一依据,即如果两要素间存在着某种关系,在模型中它们之间就会有连线相接,否则两个要素就是独立的。可见,DEMATEL方法以要素间的影响程度为依据,以分析各要素的影响度、被影响度、中心度和原因度为核心,它不仅反映了
15、系统各要素间的相互影响关系及相应的影响程度,而且反映了各要素在系统中的重要程度;ISM方法是以分析系统要素间的联系(形成关系)为基础,以分析系统层级化结构为核心,并可以进一步进行系统的序列化、聚类分析等,因此它不需要具体因素间的数量化关心。联系:从数学模型可以看出,DEMATEL方法和ISM方法并不是两个彼此相互独立的方法,两个基本矩阵的转化是这两种方法相结合进行系统分析的基础。DEMATEL方法中的直接影响矩阵可以转化为ISM方法中的邻接矩阵。DEMATEL方法中的综合影响矩阵可以转化为ISM方法中的可达矩阵。DEMATEL方法与ISM方法各有侧重,但是研究的每一方面都是系统分析的重要组成部
16、分。第六题:食物链的分析求得A的可达矩阵M如下:计算系统的可达性集合和先行集合,以及二者的共同集合,如下表所示:iL(si)F(si)L(si)F(si)11-121122,11,121,2233,10-121,3344,11,121,4455-121,5566-121,5,6677-121,5-7788,10,121,5-8899,11,121,5-7,991010,121,3,5-8,101011111-7,9,111112121-10,1212由上表可知: T=1,由于L(s1)L(si),所以s1与其他要素属于同一个区域中,即M不可分区。 s11、s12是第一级要素。划去s12对应的行
17、列,所得表如下:iL(si)F(si)L(si)F(si)11-1011221,2233,101,33441,4455-101,5566-101,5,6677-101,5-7788,101,5-88991,5-7,9910101,3,5-8,1010由上表可知:s2、s4、s9、s10是第二级要素。同理可得:s3、s8是第三级要素; s7是第四级要素;s6是第五级要素;s5是第六级要素;s1是第七级要素。按照上面所划分的要素等级,对可达矩阵M进行重新排列,得到如下:对于食物链的例子,用DEMATEL方法,算的综合影响矩阵T如下:要 素s1s2s3s4s5s6s7s8s9s10s11s12s10
18、0.25000.25000.25000.25000.06250.07810.08200.09770.08300.21190.2327s200000000000.25000.2500s30000000000.25000.25000.3125s400000000000.25000.2500s5000000.25000.31250.32810.39060.08200.09770.1182s60000000.25000.06250.31250.01560.07810.0820s700000000.25000.25000.06250.06250.0781s80000000000.250000.0625s
19、900000000000.25000.2500s10000000000000.2500s11000000000000s12000000000000各系统要素的原因度以及中心度要素行和列和原因度中心度s11.8479 0.0000 1.8479 1.8479 s20.5000 0.2500 0.2500 0.7500 s30.8125 0.2500 0.5625 1.0625 s40.5000 0.2500 0.2500 0.7500 s51.5791 0.2500 1.3291 1.8291 s60.8008 0.3125 0.4883 1.1133 s70.7031 0.6406 0.0625 1.3438 s80.3125 0.7227 -0.4102
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 金矿尾矿处理与资源化利用技术考核试卷
- 酿造食品企业的法律法规遵守与合规考核试卷
- 慢性阻塞性肺疾病疾病查房
- 急救仪器使用与维护指南
- 急性呼吸窘迫综合征护理要点
- 呼吸机脱机指征标准
- Cladosporide-C-生命科学试剂-MCE
- 2025年新高考数学一轮复习讲义(学生版)
- 食品饮料行业2025年包装废弃物处理与资源化利用研究报告
- 2025年睡眠医疗市场趋势预测:诊疗服务模式创新与行业可持续发展路径
- 耳穴治疗学试题及答案
- 2024版压力容器设计审核机考题库-简答题3-1
- 抗肿瘤临床应用管理办法
- 施工现场脚手架搭设的示例图解
- 2024年甘肃兰州中考满分作文《向内扎根向阳而生》
- 肝性脑病的临床观察与护理
- 2025(统编版)语文五年级下册第八单元解析+任务目标+大单元教学设计
- 《责任和担当》课件
- 涉外合同审查培训
- 2025年度医疗健康咨询服务兼职医生聘用合同
- 2025年外卖平台合作配送合同范本
评论
0/150
提交评论