系统工程课程设计运用ISM优化大城市公共交通发展Word文档格式.docx
《系统工程课程设计运用ISM优化大城市公共交通发展Word文档格式.docx》由会员分享，可在线阅读，更多相关《系统工程课程设计运用ISM优化大城市公共交通发展Word文档格式.docx（11页珍藏版）》请在冰点文库上搜索。
但是公共交通作为一个系统工程,其优化方法和措施很多,很难直观地区分措施的重要程度,故在进行公交优化时确定优化措施的主次及实施先后等问题存在一定的难度。
为了在短时间内更有效地采取措施,分清主次,使公交发展更切实可行,促进公共交通的良好发展,在此次课设中采用解释结构模型来优化大城市的公共交通。
2.2分析问题
影响南京市公共交通发展的因素很多,根据实际情况和参考资料进行相应的分析,对优化措施进行归纳和总结,其构成要素见表2.1。
2.3该问题的调查问卷
通过调差问卷的形式,可以使问题现实化,问题结论更有可信度。
在调差问卷的过程中能掌握实际生活中的实际的问题,在对实际问题的实际调查与研究过程中,运用具体的方法解决具体的问题,是具体问题具体化,最终找到最优的解决方案。
调查问卷见附录
（一）。
2.4ISM的建立
1.系统中这12个要素是有机的联系在一起的,而这些要素之间又是相互影响,相互作用的,将这种影响及其作用关系用矩阵、及邻接矩阵来表示出来。
矩阵的元素aij1表示要素Ai对Aj有直接影响,否则aij0。
在对本问题的系统分析中,建立邻接矩阵如表2.2。
表2.1系统的构成要素
要素编号要素名称要素定义
A1票价体系各站点区间内票价的构成体系
A2公交运营成本公共交通在运营中产生的成本
A3公交站点优化使公共交通站点合理布局的过程
A4快速公交发展规划使公共交通快捷、迅速的一系列发展计划
A5公交专用道只允许公交通过的线路
A6公交投资力度对公交投资多与少的一个判断
A7公交换乘枢纽乘客换乘公交的大型节点
A8公交优先信号控制对公交优先通过的一种信号的控制
A9公交运营车辆技术水平公共交通运营车辆的技术水平
A10限制私家车发展政策使私家车合理发展的相关政策
A11公交司乘人员素质公交司机和乘客的素质
A12公交优先法律体系关于公交优先的法律体系
表2.2邻接矩阵
A1A2A3A4A5A6A7A8A9A10A11A12
A1000000000000
A2000000000000
A3110000000000
A4000010010000
A5000000000000
A6000000001000
A7000010000000
A8000000000000
A9000100000000
A10001000100000
A11000100000000
A12000000000010
2.5解决问题
在此设计过程中,为了使复杂问题简单化,明晰化,我们运用解释结构模型法（InterpretativeStructuralModelingMethod,简称ISM方法）解决问题。
下面对此种方法做以全面的介绍。
2.5.1ISM解释结构模型叙述
解释结构模型法（InterpretativeStructuralModelingMethod,简称ISM方法）是现代系统工程中广泛应用的一种分析方法,能够利用系统要素之间已知的零乱关系,用于分析复杂系统要素间关联结构,揭示出系统内部结构。
核心思想:
把复杂系统分解为若干子系统要素,利用人机交互,将系统构造成一个多级递阶的结构模型,如图2.1所示。
图2.1递阶层次结构
ISM的应用:
ISM特别适用于变量众多、关系复杂而结构不清晰的系统分析,也可用于方案的排序。
ISM的应用十分广泛,从能源问题到地区经济开发、企事业甚至个人范围的问题,都可用ISM来建立结构模型,并据此进行系统分析。
物流领域:
质量工程项目、业务流程再造、制造企业ERP影响因素分析等。
1.解释结构模型的工作程序如下:
（1）建立系统要素关系表;
（2）根据系统要素关系表,作相应有向图,并建立邻接矩阵;
（3）通过矩阵运算求出该系统的可达矩阵M;
（4）对可达矩阵M进行区域分解和级间分解;
（5）建立系统解释结构模型。
2.系统结构的矩阵表达:
1邻接矩阵:
表示系统要素间基本二元关系或直接联系情况的矩阵。
2可达矩阵:
表示系统要素间任意次传递性二元关系或有向图上两个节点之间通过任意长的路径可以到达的情况。
图2.2有向图
图2.3可达矩阵图
3.可达矩阵的计算:
（1）邻接矩阵+单位矩阵新矩阵即A+IA+I
（2）依次运算:
A+I1≠A+I2≠A+I3≠?
?
≠A+Ir-1A+IrM即当A+Ir-1A+Ir时,矩阵A+Ir-1就是可达矩阵其中运算中用到的布尔代数法则为:
0+00,0+11,1+11
0×
00,1×
11
4.建立递阶结构模型的规范方法:
建立反映系统问题要素间层次关系的递阶结构模型,在可达矩阵的基础上进行,一般要经过区域划分、级位划分、骨架矩阵提取和多级递阶有向图绘制等四个阶段。
2.5.2ISM建模过程
1.区域划分
区域划分即将系统的构成要素集合,分割成关于给定二元关系的相互独立的区域的过程。
首先以可达矩阵M为基础,划分与要素Si（i1,2,…,n）相关联的系统要素的类型（如可达集、先行集等）,并找出在整个系统（所有要素集合S）中有明显特征的要素。
有关要素集合的定义如下:
①达集RSi:
在可达矩阵或有向图中,由Si可到达的诸要素所构成的集合,记为R（Si）。
②先行集ASi:
在可达矩阵或有向图中,可到达Si的诸要素所构成的集合,记为A（Si）。
③共同集CSi:
可达集和先行集的共同部分,即交集,记为CSi;
系统要素Si的可达集RSi、先行集ASi、共同集CSi之间的关系如图2.1所示:
图2.1关系图
④起始集BS和终止集ES:
起始集:
是在S中只到达其他要素而不被其他要素到达的要素所构成的集合,记为B（S）。
B（S）中的要素在有向图中只有箭线流出,而无箭线流入,是系统的输入要素。
判断方法:
当CSiASi时,Si即是起始集的元素。
终止集:
当CSiRSi时,Si即是终止集的元素。
得到以上特征集后判断系统要素集合S是否可分割方法有两种:
（1）判断起始集BS中的要素及其可达集RSi要素能否分割;
（2）判断终止集ES中的要素及其先行集ASi要素能否分割;
重点介绍利用起始集进行判断的方法:
利用起始集B（S）判断区域能否划分的规则如下:
在B（S）中任取两个要素bu、bv:
①如果R（bu）∩R（bv）≠ψ,则bu、bv及R（bu）、R（bv）中的要素属同一区域。
若对所有u和v均有此结果（均不为空集）,则区域不可分。
②如果R（bu）∩R（bv）ψ,则bu、bv及R（bu）、R（bv）中的要素不属同一区域,系统要素集合S至少可被划分为两个相对独立的区域。
区域划分的结果可记为:
∏（S）P1,P2,…,Pk,…,Pm。
其中Pk为第k个相对独立区域的要素集合。
相应的经过区域划分后的可达矩阵变为块对角矩阵,记作M（P）。
2.级位划分
区域内的级位划分,即确定某区域内各要素所处层次地位的过程。
这是建立多级递阶结构模型的关键工作。
设P是由区域划分得到的某区域要素集合,若用L1,L2,…,L表示从高到低的各级要素集合（其中为最大级位数）,则级位划分的结果可写出:
∏（P）L1,L2,…,L。
级位划分的基本做法是:
找出整个系统要素集合的最高级要素（终止集要素）后,可将它们去掉,再求剩余要素集合（形成部分图）的最高级要素,依次类推,直到确定出最低一级要素集合（即L）。
即找到共同集等于可达集的要素,CSiRSi
3.提取骨架矩阵
提取骨架矩阵,是通过对可达矩阵M（L）的缩约和检出,建立起M（L）的最小实现矩阵,即骨架矩阵A′。
这里的骨架矩阵,也即为M的最小实现多级递阶结构矩阵。
对经过区域和级位划分后的可达矩阵M（L）的缩检共分三步,即:
（1）检查各层次中的强连接要素,建立可达矩阵M（L）的缩减矩阵M′L
（2）去掉M′L中已具有邻接二元关系的要素间的超级二元关系,得到经进一步简化后的新矩阵M?
L。
（3）进一步去掉M?
（L）中自身到达的二元关系,即减去单位矩阵,将M?
（L）主对角线上的“1”全变为“0”,得到经简化后具有最小二元关系个数的骨架矩阵A′。
4.绘制多节递阶有向图
根据骨架矩阵A′,绘制出多级递阶有向图D（A′）,即建立系统要素的递阶结构模型。
绘图一般分为如下三步:
1.分区域从上到下逐级排列系统构成要素。
2.同级加入被删除的与某要素有强连接关系的要素（如例1中与S4强连接的S6）,及表征它们相互关系的有向弧。
3.按A′所示的邻接二元关系,用级间有向弧连接成有向图DA′
以可达矩阵M为基础,以矩阵变换为主线的递阶结构模型的建立过程:
2.5.3ISM方法建模过程?
规范方法
在系统结构不十分复杂的情况下,可以采用简便的方法来建模。
主要过程:
1.判定二元关系,建立可达矩阵及其缩减矩阵
已知一有向图如下图2.2所示:
图2.2有向图由有向图2.2得方格图2.3,如下图所示:
图2.3方格图其中:
A?
列要素对行要素有直接影响;
V?
行要素对列要素有直接影响;
X?
行列两要素相互影响;
（）?
逻辑推断递推关系。
在问题确定后,相关人员根据问题要素绘制方格图,直观地确定各要素之间的二元关系,并在两要素交汇处用不同符号加以标示。
再由逻辑推断出要素间各次递推的二元关系,用加括号的符号注在图上。
由方格图可以得到邻接矩阵如下:
经计算有A+I1≠A+I2A+I3,所以可达矩阵就是A+I2,常用M表示
2.对可达矩阵的缩减矩阵进行层次化处理
删除强连接的要素,得到缩减矩阵;
对缩减矩阵进行重排,按每行“1”元素的多少,由少到多顺序排列,调整行和列;
在新矩阵中,从左上角到右下角,依次分解出最大阶数的单位矩阵,并加注方框。
每个方框就表示一个层次。
3.绘制多级递阶有向图
先把所有要素按已有层次排列,然后按照M′中两方框交汇处的“1”元素,画出不同层次要素间直接联系的有向弧,得到多级递阶有向图。
根据系统要素建立的邻接矩阵,使用MATLAB编程求出可达矩阵和级别划分矩阵。
程序见附录,运行截图见图2.4。
图2.4程序运行截图
源程序见附录。
3.根据可达矩阵得到可达集、先行集和共同集如表2.3所示。
表2.3可达集、先行集和共同集列表
Ai可达集RAi先行集AAi共同集CAi
A111,3,101
A222,3,102
A31,2,33,103
A44,5,84,6,9,11,124
A554,5,6,7,9,10,11,125
A64,5,6,8,966
A75,77,107
A884,6,8,9,11,128
A94,5,8,96,99
A101,2,3,5,7,101010
A114,5,8,1111,1211
A124,5,8,11,121212
根据级别划分矩阵绘制多级制递阶有向图,见图2.5。
图2.5多级制递阶有向图
由图2.5中的多级制递阶有向图可知,公交结构最优系统是一个具有四层的多级递阶系统。
影响因素层次分析如表2.4所示。
公交结构优化的最直接影响因素是:
A1-票价体系,A2-公交运营成本,A5-公交专用道,A8-公交优先信号控制;
影响第二层的因素是:
A3-公交站点优化,A7-公交换乘枢纽,A4-快速公交发展规划;
影响第三层的因素为:
A10-限制私家车发展政策,A11-公交司乘人员素质,A9-公交运营车辆技术水平;
影响第四层的因素为:
A12-公交优先法律体系,A6-公交投资力度。
表2.4影响因素层次分析
优化措施的层次优化措施
第一层A1、A2、A5、A8
第二层A3、A7、A4
第三层A10、A11、A9
第四层A12、A6
我们可以清楚的看到该系统是一个有四级的递阶结构模型。
由此可以分析出,大城市公交优化的相关优势如下:
通过建立完善的公交优先的法律体系来保障相关政策和措施的制定,加大投资力度。
与此同时,提高公交运营车辆的技术水平,制定限制私家车辆发展政策,提高司乘人员素质。
随之进行快速公交发展规划,进行公交站点优化,建立公交换乘枢纽。
在此基础上采取相应的具体方法来优化公交系统,如:
设立公交专用车道,进行公交优先信号控制,建立合理的票价体系,进行公交运营成本核算等。
3.结论
通过运用解释结构模型（InterpretativeStructuralModelingMethod,简称ISM方法）来优化大城市公共交通的发展,确定影响公共交通发展的各种因素的各级层次通
过逐层逐次的详细分析和计算,明确地知道,公共交通发展最直接的影响因素是票价体系、公交运营陈本、公交专用道、公交优先信号控制。
明确了影响因素的层次并针对公交优先模型的多级递阶结构提出了相应的公交优先措施,为衡量公交优化措施的主次,在短时间内有效地、有侧重地实施公交优化措施提供了参考,为大城市公共交通的进一步优化发展提供了有力的依据。
只是,本次课设只对影响因素进行了定性分析,没有体现影响因素对公共交通发展的影响程度的量化,即没有定量分析。
因此,作为二十一世纪的大学生,在今后的学习或工作中需结合其他算法,如模糊综合评价法、网络分析法等进行定量的分析与评价,对所考察的问题进行更深层次的的研究和讨论,以此来提升自身在专业知识方面的扩充。
同时,通过本次的课程设计,作为只在理论的知识里生活的大学生,我们更进一步掌握了实践定性分析的方法与步骤,熟悉了理论与实践相结合的具体分析过程,初步学会运用理论知识进行实际案例的评价与决策,并且能够结合实际了解系统评价指标体系的相关结构组成。
同时在各小组团结配合共同完成实验任务时让我们明白团队合作的关键之所在,以及明白了众人的力量是强大的。
4.参考文献
[1]汪应洛.系统工程[M].北京:
机械工业出版社,2004:
69-81.
[2]姜启源,谢金星,叶俊.数学模型[M].高等教育出版社,2011.
[3]刘思峰,党耀国.预测方法与技术[M].高等教育出版社,2005.
[4]田志友,王浣成.解释结构模型在服务蓝图设计中的应用[J].工业工程与管理,2004,
[5]穆尔,高会生,刘童娜,李聪聪.MATLAB实用教程[M].电子工业出版社,2010.
附录:
使用MATLAB编程求出可达矩阵和级别划分矩阵的程序如下。
%ISM解释型结构模型根据邻接矩阵求可达矩阵,进行级划分的算法:
%step1求可达矩阵
%A是邻接矩阵,可以由用户输入
A[100000000000010000000000111000000000000110010000000010000000000001001000000010100000000000010000000100001000001000100100000100000010000000000011];
NsizeA,1;
%N是矩阵的阶数,是所有要素的数目
r[N:
N];
rA;
forn1:
1:
Nfori1:
Nforj1:
Nsum0.0;
fork1:
Nsumsum+ri,k*Ak,j;
%此处是采用普通的矩阵乘法,但是下面可以根据sum的值是否大于或等于1%来判断连通型,从而等价得到ri,j的实际值。
endifsum1ri,j1;
elseri,j0;
endendendendRr;
%R是可达矩阵
%step2级别划分L[N:
%L是二维数组,存储级划分的结果,下面首先初始化为0fori1:
NLi,j0;
endend
forp1:
N%p是存储层次级数l的变量,也控制了总循环的次数k1;
%k是记录每一级p内要素,在L内存储下标的变量,形式为Lp,k%下面是利用二重循环,求解p级内的要素fori1:
Nsign0;
%是标志变量,初始为0,如果对某要素考察后,其值仍为0,则表明该要素是顶点要素。
否则,不是顶点要素sum0;
%sum是计数器,用于判断当前考察的要素所在的矩阵行向量是否是值全为0的向量。
因为本程序的算法是这样的,%即如果已经发现某要素是某级的顶点,则在求下一级顶点时,需要去除上面各级的要素。
在本程序是通过变通的方法实现%相同的目的,即所有已经是顶点的要素的所在行和列的值全部重新置为0。
这样就等价于删除了这些上级顶点。
从而简化了%程序的算法。
forj1:
Nsumsum+ri,j;
ifri,j1&
ri,j~rj,i%算法的关键,RiRi交Ai,只需要判断要素i所在的行中,%所有值为1的矩阵元素Ri,j,其对称矩阵元素Rj,i的值如果也为1,则说明RiRi交Ai成立。
否则%只要有一个元素不满足该条件（此时标志变量的值赋为1）,就说明该要素i不是顶点。
进行下一要素i+1的考察。
sign1breakendendifsum~0%sum不为0,说明该要素不是已经求出的上级的顶点,可能是新级别的顶点,如果sum为0,则说明是已经%求出的顶点ifsign0%sign为0,说明是新顶点Lp,ki;
%在L内记录新顶点,p为级别kk+1;
end%已经对p级别的顶点计算完毕,下面的程序是把p级的顶点所在的行和列的矩阵元素全部置为0,以在后面的计算表示老顶点forg1:
NifLp,g~0forh1:
NrLp,g,h0;
rh,Lp,g0;
endendend%进行下一次循环,找出p+1级别的顶点元素endclc%清命令窗口%输出可达矩阵,级别划分矩阵。
RL
%对于同级内部的强连通集的最大回路集合的划分,和不连通集的划分,比较容易。
程序不再给出算法,自行判断。
升级会员 