1、第三章第三章 系统模型系统模型系统模型概述结构模型层次分析法模型概念及特征系统模型的分类建模原则及常用方法结构模型概念及特征解析结构模型的建立应 用 案 例第1节 系统模型概述一切客观存在的事物及其运动形态称为“实体”（即原型）。为便于实验、分析和预测，总是先把所需研究的系统结构型态或运动形态变成易于考察的形式，即转化为“模型”。1 1 系统模型定义系统模型定义系统模型是对现实系统（实体）的特征及其变化规律的一种模仿、抽象或描述。模仿、抽象或描述。一、系统模型及其特征一、系统模型及其特征系统的属性是多方面的，系统模型只是系统某一方面本质属性的描述，所以同一系统或试题，模型不是唯一的；模型建立是
2、以模型与原型之间的相似性为基础的，这里的相似可以是外表的相似，内部结构的相似或仅为功能的相似。模型可以是定量的，也可以是定性的，或是两者的结合模型。2 2 系统模型的特征系统模型的特征它是现实系统的抽象或模仿；它是由反映系统本质或特征的主要因素构成；它集中体现这些主要因素之间的关系。说明说明3 3使用系统模型的必要性使用系统模型的必要性系统开发的需要。在开发一个新系统时，系统尚未建立，无法直接实验；经济性考虑。大型复杂系统直接实验价格昂贵；安全性考虑。有些系统直接实验是很危险的，有时根本不允许；时间上考虑。社会、经济、生态系统，惯性大，反应周期长；系统模型易操作，分析结果易于理解。这里介绍按与
3、实体的关系分：1 1 形象模型（实体与比例模型）形象模型（实体与比例模型）这种模型保留着实体的外形特征，仅在尺度上成比例的改变。2 2 模拟模型模拟模型根据相似系统原理，利用一种系统代替或近似描述另一种系统，前者为后者的模拟模型。3 3 数学模型数学模型用各种数学符号、数值描述工程、技术、管理、经济等有关因素及它们之间数量关系的模型。包括网络模型、图表模型、逻辑模型和解析模型。二、系统模型的分类二、系统模型的分类现实系统形象模型模拟模型文字模型网络模型图表模型逻辑模型解析模型物理模型数学模型研究的速度变化方便性抽象性现实性建模费用建模时间增加减少系统模型的分类及特征比较系统模型的分类及特征比较
4、1 1 建模的原则建模的原则（1）现实性：把本质的东西和关系反映进去，非本质的东西去掉，而又不影响反映现实的真实程度。（2）简明性：模型既要精确，又要简明。（3）适应性：在运算分析方面、适应问题的变化、操作方面等具有适应性。（4）完整性（5）规范性：尽量借鉴标准形式。一般处理原则：力求达到真实性，在真实的基础上达到简明性，最后尽可能达到适应性要求。三、建模的原则及常用方法三、建模的原则及常用方法2 建模一般过程（1 1）明确建模目的和要求；）明确建模目的和要求；（2 2）弄清系统或子系统中的主要因素及）弄清系统或子系统中的主要因素及其相互关系其相互关系 ；（3 3）选择模型方法；）选择模型方法
5、；（4 4）确定模型结构；）确定模型结构；（5 5）估计模型参数；）估计模型参数；（6 6）模型试运行；）模型试运行；（7 7）对模型进行实验研究；）对模型进行实验研究；（8 8）对模型进行必要修正。）对模型进行必要修正。3 3 常用建模方法常用建模方法直接分析法直接分析法：对内部结构和特征已经清楚的系统，可利用已知的定律和定理，经过一定的分析和推理，得到系统模型。如线性规划模型、存储模型。数学关系式表达数学关系式表达实验法实验法：对内部结构和特征不清楚或不很清楚的系统，如能进行实验观察，可通过实验方法测量其输入和输出，再根据一定的辨识方法，得到系统模型。模拟法模拟法统计分析法：统计分析法：系
6、统结构不很清楚，且不允许直接进行实验的系统，可以采用数据收集和统计分析的方法建立系统模型。类似法类似法：建立系统的类似模型。拟合法拟合法启发性思考法。启发性思考法。设一个质量为m，长度为l的摆，其偏离中心线的角度为(很小)，(t)st:方程的解是以 为周期的简谐震动。022mgmldtdGLT2建立单摆简谐运动的类似模型建立单摆简谐运动的类似模型mglL-C电路，电路中q(t)st:解是以 为周期的简谐震动。一一对应模拟。CLLCTqLLCdtqd20122)()(1ttqglLCL-C电路图蒙特卡罗的特点是在所研究系统的模型中模拟随机事件，即对于所求的值应该设定什么样的概率过程为题进行求解的
7、技术方法。启发性思考法启发性思考法蒙特卡罗法计算值蒙特卡罗法计算值11在边长为1的的正方形中任意打N个点，并将n个点置于扇形部分，如使点数N足够大，则认为近似等于正方形和扇形面积之比，即：N/n=12/(12 1/4)即：4n/N与概率现象本身没有任何关系的问题，也可用概率的方法来解决，是一种“想法的转换”，即启发性思考方法。第第2 2节节 结构模型结构模型（Structure Model）在开发和改造一个系统时，首先需要了解系统中各要素间存在怎样的关系，即了解和掌握系统的结构，即建立系统的结构模型。1 1 结构模型结构模型就是用有向连接图来描述系统各要素间的关系，以表示一个作为要素集合体的系
8、统模型。一结构模型的概念及原理一结构模型的概念及原理（1）结构模型是一种几何模型：节点表示系统的要素，有向边表示要素间的关系。（2）结构模型是以定性分析为主的模型。（3）结构模型可以用矩阵形式描述，进行定性与定量分析。结构模型的建模方法很多如关联树（如问题树、关联树（如问题树、目标树、决策树）法、解释结构模型化（目标树、决策树）法、解释结构模型化（ISMISM）方法、系统动力学（方法、系统动力学（SDSD）结构模型化方法等。结构模型化方法等。本部分要求大家主要学习和掌握本部分要求大家主要学习和掌握解析结构模型ISM ISM（Interpret Structure Model）方法方法（规范方法
9、、实用化方法）。（规范方法、实用化方法）。S3S1S2S4S5基本性质基本性质ISM是美国华费尔特教授于1973年作为分析复杂的社会经济系统有关问题而开发的一种方法。其特点是把复杂的系统分解为若干子系统或要素，利用人们的实践经验和知识，以及计算机的帮助，最终将系统构造成多级递阶的结构模型。ISM的程序为：组织构造ISM小组（10人左右）设定问题选择系统要素，制定系统明细表。构思有向图，建立连接矩阵和可达矩阵。对可达矩阵进行分解，建立结构模型。由结构模型转化为解析结构模型。2 2 解析结构模型原理解析结构模型原理设定设定问题问题、形、形成意成意识模识模型型找出找出 影响影响 要素要素要素要素关系
10、关系分析分析（关（关系图系图）建立可建立可达矩阵达矩阵(M)和缩和缩减减 矩阵矩阵（M/）矩阵矩阵层次层次化处化处理理(ML/)绘制绘制 多级多级 递阶递阶 有向有向 图图建立建立 解释解释 结构结构 模型模型分析分析 报告报告比较比较/F 学习学习初步分析初步分析规范分析规范分析综合分析综合分析ISMISM方法原理图方法原理图1 1有向连接图有向连接图由若干节点和有向边连接而成的图象，即为节点和有向边的集合。表示为：G=S，E2 2邻接矩阵邻接矩阵A A描述图中节点两两之间的直接关系。A中元素3 3可达矩阵可达矩阵R R用矩阵形式反映有向连接图各节点之间通过一定路径可以到达的程度。Si经若干
11、路径到达Sj 否则二、解析结构模型的建立二、解析结构模型的建立jijiijsRsRssa,0,101ijr可达矩阵=邻接矩阵A+单位矩阵I，并经过一定的运算后求得。即有 A1=A+I再设 A2=(A+I)2(用布尔代数运算规则)一般地，通过依此运算后，可得：A1 A2 An-1=An则有R=An-1=(A+I)n-1R-可达矩阵，它表明各节点间经过长度不大于(n-1)条通道可以到达的程度。对于节点数n为个的图，最长的通路长度肯定不超过(n-1).例：例：现有如下图所示7个要素组成的系统，试建立它的关系，并求邻接矩阵和可达矩阵。有向连接图7154632由此可得邻接矩阵AA的元素全为零的行所对应的
12、节点为汇点。A的元素全为零的列所对应的节点为源泉。对应每一节点的行中元素值为1的数量，是离开该节点的有向边数。对应每一节点的列中元素值为1的数量，是进入该节点的有向边数。0000010000100000000000110000000100000000010000000A矩阵矩阵A A的特性的特性建立可以矩阵R。经计算后得：(A+I)1 (A+I)2=(A+I)3 R=(A+I)2 1000011011100000100000111000011110000000110000001R#布尔代数运算规则：0+0=0，0+1=1，1+0=1，1+1=1，0 0=0，0 1=0，1 0=0，1 1=14
13、 4 可达矩阵的分解可达矩阵的分解区域分解1(S)将要素分成区域，不同区域的要素相互间是没有关系的。首先将R中的元素划分为可达集和先行集（1）要素Si的可达集R(Si)R中第Si行矩阵元素为1对应的列要素的集合。即：(N为节点集合，rij=1表示 Si 与Sj关联)1)(ijjirNSSR区域分解区域分解（2）要素Sj的先行集A(Sj)R中第Sj 列矩阵元素为1所对应的行要素的集合。即：（3）共同集合T可达集R(Si)与先行集A(Sj)的交集等于先行集A(Sj)的要素集合，即：1)(ijijrNSSA)()()(jjiiSASASRNST（4）确立不同区域 任取属于共同集的两要素Su,Sv，若
14、 ，则Su,Sv属同一区域；若 ，则Su,Sv属于不同区域。这样运算后的集合称区域分解，可写成：其中M为区域数。)()(vuSRSR)()(vuSRSRPmPPS,)(21级间分解2(P)将系统中的所有要素，以可达矩阵为准则划分不同层次。在一个多级结构中，它的最上层要素Si的R(Si)，只能由Si自身和Si的强连通要素组成；同时Si的先行集只能由由Si自身和结构中的下一级可能到达的要素以及Si的强连通要素组成。若Si是最上层单元，需满足：找出最高一级要素后，将其从可达矩阵中划去相应的行与列，在从剩下的可达矩阵中寻找新的最高级要素，依此类推。级间分解级间分解)()()(jiiSASRSR级间划分
15、可用下式表示：，其中K为级次若定义：L0=，则：其中：分别是由 要素组成的子图求得的可达集和先行集。强连通划分3(L)：级间分解后，每级要素中可能有强连通要素，一般构成一个回路，只需选择一个要素即可。kLLLP,)(212)()()(111110ikjkikkikSRSASRLLLPSL)(),(11jkikSASR110kLLLP强连通划分强连通划分接例接例 可达矩阵分解（区域划可达矩阵分解（区域划分）分）I=(j)R(Si)A(Sj)R(Si)A(Sj)T=A(Sj)R(Si)A(Sj)=R(Si)123456711,23,4,5,64,5,654,5,61,2,71,2,72,733,4
16、,63,4,5,63,4,671234,654,673715因为：R(3)R(3)A(7)=A(7)=，则S3，S7分属不同区域，所以，区域划分为：7216543211,)(sssssssPPP因为：S1，S5满足：所以，S1，S5分属两区域的最高层次。即；L1=S1，S5再有N-L0 L1进行第二级分解。接例接例 可达矩阵分解（级间分可达矩阵分解（级间分解）解）)()()(jiiSASRSRi=(j)R(Si)A(Sj)R(Si)A(Sj)2346723，4，64，64，62，72，733，4，63，4，67234，64，67该表的最高级，即为可达矩阵的第二级要素L2=2,4,6由N-L0-
17、L1-L2,得：i=(j)R(Si)A(Sj)R(Si)A(Sj)37373737该表的最高级，即为可达矩阵的第三级要素为：L3=3,7这样，经过三级划分，将R中的7个单元划分成三层次，即2 2(P)=L(P)=L1 1，L L2 2，L L3 3(强连通划分强连通划分)44，6 6 属强连通块。属强连通块。S1S2S7S3S4S5S6第第1级级 第第2级级 第第3级级区域区域划分划分级位级位划分划分强连接强连接 要素要素 缩减缩减剔除剔除超级超级关系关系去掉去掉自身自身关系关系绘图绘图结束结束以可达矩阵以可达矩阵M M为基础，以矩阵变换为主线的递阶结构模型的建立过程：为基础，以矩阵变换为主线
18、的递阶结构模型的建立过程：M MM M（P P）M M（L L）M M（L L）M M（L)A DL)A D（AA）原例的递阶结构模型：原例的递阶结构模型：案例：人口系统影响总人口增长问题案例：人口系统影响总人口增长问题新中国成立以来，人们的期望寿命有了较大提高，相对死亡率降低了，国民收入的不断增长，生活水平不断提高，计划生育政策贯彻不力等等，导致我国人口速度增长过快。为此，成立了各方面人员参加的研究小组对人口增长问题进行了研究，主要任务为：应用应用ISMISM讨论和确定我国总人口增长的影响因素；讨论和确定我国总人口增长的影响因素；根据经验和对话建立可达矩阵，解析结构模型；根据经验和对话建立可
19、达矩阵，解析结构模型；通过模型中各因素分析，为制定有关人口政策、通过模型中各因素分析，为制定有关人口政策、控制人口等政策提供依据。控制人口等政策提供依据。经ISM小组讨论后，认为主要影响因素有11个，并经多次讨论后确定它们之间的关系。三、解析结构模型的应用三、解析结构模型的应用 S1期望寿命S2保健水平S3生育能力S4计生政策S5思想风格S6营养S7环境污染S8国民收入S9国民素质S10出生率S11死亡率S12总人口：S Si i 与与S Sj j互互有关系；有关系；：S Sj j与与S Si i有有关系；关系；：S Si i 与与S Sj j有有关系关系1111111111111111111
20、111111111111111111111111111根据以上对话过程，建立可达矩阵根据以上对话过程，建立可达矩阵1234567891011121,11,121,2,3,11,123,10,124,10,121,3,6,10,11,121,7,11,121,3,4,8,10,11,124,9,10,1210,1211,12121,2,6,7,822,3,6,84,8,967893,4,6,8,9,101,2,6,7,8,101-12123456789101112I=jR(Si)A(Sj)R(Si)A(Sj)R(2)R(6)R(7)R(8)R(9)共同集合不存在空集，所以没有区域之分。首先找出R
21、（12）=R（12）A（12）所以第一层次为要素12第二层次为要素10，11第三层次为要素1，3，4第四层次为要素2，6，7，8，9总人口出生率死亡率生育能力思想风俗期望寿命计生政策保健水平营养水平国民收入环境污染国民素质人口系统解析结构模型人口系统解析结构模型第第3 3节节 层次分层次分析法析法它是由美国著名数学家Satty1972年提出的，这是一种以定性分析与定量分析相结合的系统分析方法，主要特点为：分析思路清晰，可将分析人员的思维过程系统化、数学化和模型化；分析时所需的定量数据较少，但要求对问题的本质、包含的因素及其内在的关系分析清晰；多用于多准则、多目标问题以及其他各种不同类型的系统分
22、析与评价。层次分析法首先把复杂系统所包含的因素及相互关系形成多层次结构模型，一般可用解析结构模型建立。多层次结构模型一般有三类：1 1 完全相关性结构完全相关性结构特点是上一级的每一要素与下一级的全部要素相关，即上一级每一要素都作为下一级的评价项目而起作用。2 2 完全独立性结构完全独立性结构特点是上一级要素都各自有独立的、完全不相同的下级要素。3 3 混合结构混合结构是一种即非完全相关又非完全独立的结构。一、建立层次结构模型一、建立层次结构模型1 1 建立判断矩阵建立判断矩阵判断矩阵是以上一层的某一要素ak作为准则，对本级的要素进行两两比较而确定的矩阵。如A层元素ak与下一层次元素B1，B2
23、，Bn有关系，则其判断矩阵为二、判断矩阵二、判断矩阵nnnjnninijiinjnjnjniknnijbbbbbbbbbbbbbbbbBBBBBBBBabB21212222211112112121计算要素相对重计算要素相对重要度的基本信息要度的基本信息bij表示对ak而言，Bj元素与Bi元素相对重要性的比值，即 bij=wi/wjbij 值是根据资料数据，专家意见和系统分析人员的经验反复研究后确定的。由于是对单一准则的两两比较，一般不难给出评分数据。2 2 判断矩阵标度判断矩阵标度(1-9)(1-9)及其定义及其定义 标度标度 定义定义1 对ak而言，Bi与Bj同样重要3 对ak而言，Bi与B
24、j稍重要一些5 对ak而言，Bi与Bj明显重要7 对ak而言，Bi与Bj非常重要9 对ak而言，Bi与Bj特别重要2、4、6、8其重要性介于上述两相邻标度的中间倒数 Bi 与Bj 的比值为bij，Bj与Bi的比值为bji=1/bij 一般对n阶判断矩阵，需对n(n-1)/2个矩阵元素给出数值。根据判断矩阵，计算对上一层次某元素，本层次各元素相对重要次序。1判断矩阵的一致性及其检验基本原理判断矩阵的一致性及其检验基本原理#（1 1）比较称重法的分析矩阵）比较称重法的分析矩阵设n个物体，重量分别为W1，W2，Wn，且，两两比较n个物体得矩阵：三、层次单排序及一致性检验三、层次单排序及一致性检验)1
25、()(212112111nnnnniiinnnijwwwwwwwwwwwwwwwwwwaA)2(,21TnwwwW)3(321nwnwnwnwAWv定义重量向量：定义重量向量：v向量向量W W右乘右乘A A得：得：v由矩阵理论知：由矩阵理论知：n n恰为矩阵恰为矩阵A A唯一非零的特征唯一非零的特征根，且为最大特征根，而根，且为最大特征根，而W W则为其所对应的特征则为其所对应的特征向量。向量。对于给定的一组物体，设想计算出它们在总重量中占份额，其方法有两种：一是直接称出每一物体，可得各自的相对重量；二是不称或无法称出每一物体重量的条件下，设法对每一物体的相对重量作出判断。对于大量的社会问题，
26、在没有统一度量标尺可用之时便可能由这种层层比较分析中估量出各因素的重要程度，这恰是比较分析判断法的优越性所在。（介绍的正是该方法）（2 2）一致性问题及检验思想）一致性问题及检验思想 式(3)要求wi/wj(i,j=1,2,n)的比值，这在物理问题中不难。但在社会性的决策问题中，尤其是在征求众多人的判断过程中，实际给出的aij判断值与理想的wi/wj有一定偏差。由矩阵理论可知，若有数1,2,n满足A的特征方程：AW=W则称1,2,n为矩阵A的特征值。且有如下性质：Aaaannnn21221121由(1)式得:已知方程(3)AW=nW的最大特征根为:max=n 作为判断矩阵作为判断矩阵A A一致
27、性的充分必要条件。一致性的充分必要条件。当矩阵A具有完全一致性时，当矩阵A不具有完全一致性，nnnniiniimax2max12max1,0,naaaniinn12211,1(3)(3)一致性指标的定义及计算一致性指标的定义及计算层次分析法用除去判断矩阵最大特征根以外的其余特征根的平均值作为判断矩阵偏离一致性的指标,即:同时引入判断矩阵的平均随机一致性指标RI值(表6-13)随机一致性比率随机一致性比率(CR)=CR)=判断矩阵的一致性指标判断矩阵的一致性指标(CI)/CI)/同阶平均随机一致性指标同阶平均随机一致性指标(RI)RI)作为一致性检验标准。当CR0.1CR0.1时,判断矩阵的一致
28、性是可以接受的。1maxnnCI2 2 层次单排序及一致性检验层次单排序及一致性检验(1)判断矩阵的特征值的近似算法将判断矩阵每一列正规化每一列经正规化的判断矩阵按行相加对向量 正规化所得的 即为所求特征向量.njibbbnkkjijij,2,1,1njbWnjiji,2,1,1,21niWWWWniWWWnjji,2,1,1TnWWWW,21(2)计算判断矩阵的最大特征根计算判断矩阵的最大特征根 表示向量（AW）的第i个分量。(3)(3)一致性检验指标一致性检验指标iniiiAWWAWn)()(11max1.01maxRICICRnnCI利用单排序结果，自上而下地求出每一级要素相对系统总体的
29、权重。四、综合排序四、综合排序njjnmninnjmjijjmimimiminjiiwwwwWwwwwwwwwwwwwwwwwwwwwBBBBCCCCWB212121222221111211212121CjWjimijiijWWW1如果C层次某些因素对于Bi单排序的一致性指标CIi相应的平均随机一致性指标为RIi，则C层次总排序随即一致性比率为：当CR0.1时，层次总排序结果具有满意的一致性。miiimiiiRIwCIwCR11减少交通事故的损失防止事故发生防止事故发生的损失促进恢复提高司机安全操作责任感提高车辆操作性能改善道路设施提高车辆安全保障功能加强十字路口交通管理充实急救医疗体制健全医
30、疗体制充实残疾人治疗培训体制完全独立性结构完全独立性结构购置一台满意的设备功 能价 格维护性方案1方案2方案3A AB BC C目标层目标层准则层准则层方案层方案层例：例：建建立立如如图图的的层层次次结结构构模模型型完全相关性结构完全相关性结构现在三种备选的方案中：C1 性能较好 价格一般 维护需一般水平C2 性能最好 价格较贵 维护需一般水平C3 性能差 价格便宜 维护容易假定该企业在设备的使用上首先要求功能强，其次要求易维护，最后考虑价格便宜。AB1 B2 B3B1B2B3 1 5 3 1/5 1 1/3 1/3 3 1B1C1 C2 C3C1C2C3 1 1/4 2 4 1 8 1/2
31、1/8 1B2C1 C2 C3C1C2C3 1 4 1/3 1/4 1 1/8 3 8 1B3C1 C2 C3C1C2C3 1 1 1/3 1 1 1/3 3 3 1一、经讨论确定各判定矩阵一、经讨论确定各判定矩阵二、层次单排序及一致性检验二、层次单排序及一致性检验B1C1 C2 C3C1C2C3 1 1/4 2 4 1 8 1/2 1/8 1和5.5 1.375 11B1C1 C2 C3C1C2C3.1818 .1818 .1818.7272 .7272 .7272.0910 .0910 .0910和ijbiiwwW/.54542.1816.27300.18180.72720.09103.0
32、1.0即判断矩阵的特征向量:W=(0.1818,.7272,0.0910)T1 1计算判断矩计算判断矩 阵特征向量阵特征向量099.30910.02728.07272.01824.21818.05456.031)(2728.01824.25454.00910.07272.01818.018/12/181424/1111maxniiinWAWAW10.00853.058.00495.00495.0133099.31maxRICICRnnCI2 2 计算最大特征根计算最大特征根3 3 一致性检验一致性检验B2C1 C2 C3C1C2C3 1 4 1/3 1/4 1 1/8 3 8 1W0.2572
33、0.07380.6690max=3.018CI=0.009CR=0.01550.10同理计算同理计算B3C1 C2 C3C1C2C3 1 1 1/3 1 1 1/3 3 3 1W0.15910.15910.6818max=3.085CI=0.0425CR=0.07330.10AB1 B2 B3B1B2B3 1 5 3 1/5 1 1/3 1/3 3 1W0.63330.10350.2532max=3.039CI=0.0195CR=0.03360.100.1818 0.2572 0.15910.7272 0.0738 0.15910.0910 0.6690 0.6818C1C2C3 0.6333
34、 0.1035 0.2532Wj0.18200.5085*0.2995 B1 B2 B3BiCjWiWij01.00741.058.0043.0043.00425.02532.0009.01035.00495.06333.011RICICRRIWRICIWCIniiiniii三、层次总排序及一致性检验三、层次总排序及一致性检验体系安全性驾驶员系统车辆系统道路系统交通环境系统管理系统性能特性技能特性视角特性身心特性行为特性教育系统救护系统监控系统管理水平夜间行车交通饱和度车辆混合率恶劣天气安全装置线形组合道路线型路面系统制动系统操稳系统防撞系统动力系统车辆悬挂系统高速公路交通安全保障体系 交 通 安 全 意 识 人 驾 驶 员 培 训、考 核 事 故 急 救 等 汽 车 安 全 设 计、制 造 交 通 车 辆 安 全 检 测、汽 车 维 保 使 用 情 况 交 通 流、交 通 速 度 特 性 交 通 安 全 体 系 道 路 设 计、修 建 道 路 道 路 养 护 道 路 安 全 保 护 体 系 法 规（公 路 法、交 通 法 规 等）交 通 环 境 自 然 环 境（气 候：风 雨 雪 雾、时 间 因 素 等）管 理 体 系、机 构 管 理 管 理 者 素 质 现 代 管 理 手 段 图2 道 路 交 通 安 全 保 障 体 系 感谢下感谢下载载
