层次分析数学建模案例教学文稿.doc

层次分析数学建模案例精品文档收集于网络，如有侵权请联系管理员删除 基于层次分析法的护岸框架最优方案选择【摘要】长期以来，四面六边透水框架在河道整治等工程中，因其取材方便、自身稳定性、透水性、阻水性好、适合地形变化等特性优点而被广泛的应用。但是，在抛投和使用过程中，存在被水流冲击而翻滚移位、结构强度的不足、难以合理互相钩连的问题，使框架群不能达到理想的堆砌效果。本文主要探讨如何合理设计改进现有护岸框架，以最大程度减少框架群被水流冲击翻滚移位的情况，增加框架群在使用过程中互相钩连程度和结构强度，达到减速促淤效果。针对问题，我们结合四面六边透水框架本身的优势特性，在原有框架的基础上进行改进设计，根据三角形稳定性的特性，通过应用机理分析，进行物理图形构造，设计出三种供选方案。模型一：构建四面六边带触脚框架模型（图5.2），该模型在四面六边透水框架的基础上，运用触脚设计，较好的融合增强四面六边透水框架本身的优点特性，使框架达到不易翻滚，并与其他的框架自然地相互钩连。模型二：构建六面九边带触脚框架模型（图5.6），该模型是对模型一的改进，综合模型一和原型模型的结构，不仅具备良好的亲水性、阻水性和稳定性，而且触脚比模型一更多，使框架更加稳定，不易翻滚、框架群之间也更容易钩连；同时，模型二施工简单，更容易构造，也更加节约经济造价成本。模型三：构建双四面六边护岸框架模型（图5.12），该模型设计内外双层四面六边透水框架体，旨在增加护岸框架结构强度和稳定性及阻水性。运用内外双层结构设计，形成内外双层保障。由三角形的稳定性可以得知该模型结构强度高、稳定性强。模型四：应用层次分析法对如何科学、合理地进行选择护岸框架，进行系统的分析。选取施工时架空率易接近4到6、结构强度、不易翻滚程度、框架群间易钩连程度、生产成本及易生产、施工简易度六个因素指标为准则层，选取原有护岸框架和本文设计的三个框架模型作为方案层，运用Matlab软件计算比较，最后得出结论为：模型二（六面九边带触脚框架模型）为最优护岸框架模型。【关键词】护岸框架 层次分析法 立体图形 触脚设计 Matlab一、 问题重述在江河中，堤岸、江心洲的迎水区域被水流长期冲刷侵蚀。在河道整治工程中，需要在受侵蚀严重的部位设置一些人工设施，以减弱水流的冲刷，促进该处泥沙的淤积，以保护河岸形态的稳定。现在常用的设施包括四面六边透水框架等。这是一种由钢筋混泥土框杆相互焊接而成的正四面体结构，常见的尺寸为边长约1m，框杆截面约0.10.1m，将一定数量的框架投入水中，在水中形成框杆群，可以使水流消能减速，达到减弱冲击，防冲促淤的效果。对四面六边透水框架在抛投时和在使用过程中，可能被水流冲击而翻滚移位，使框架群不能达到理想的堆砌效果，对功能有不利影响。为了使框架在水中互相钩连，需要设计新的形状。但已有的多数设计方案都存在问题，主要集中在两个方面：结构强度不足，以及虽然原则上能够互相钩连，但依然不清楚最终堆砌而成的形状是否合理。请你建立合理的数学模型，设计一个良好的框架结构。二、 问题的分析近10年来，四面六边透水框架作为一种新型的护岸防冲方式，在河道整治等工程中，得到了广泛的应用。这种新型的护岸方式，不仅可以节省工程材料，而且具有良好的减速促淤效果，在投放区域起到了防止冲刷的作用，但是，在实际使用过程中，四面六边透水框架可能被水流冲击翻滚移位，使框架群不能达到理想的推砌效果，对功能有不利的影响。主要表现为：结构强度的不足，以及虽然原则上能够互相钩连，但是不清楚最终形状是否合理。针对所给问题，我们要设计出一种新型的护岸框架结构，以解决以下实际生活中护岸框架存在的问题：1) 结构强度的不足；2) 科学、合理的相互钩连；3) 减轻框架群在使用过程中的翻滚移位率。我们知道混泥土四面六边透水框架群作为新型的江河护岸工程技术1，和其他传统的护岸框架（丁坝、抛石护脚等）相比，拥有取材方便、自身稳定性好、透水性好、阻水性佳、基础不易被冲刷、适合地形变化等特性优点，特别是框架的尺寸、架空率和铺设长度的设计，更是有效的综合了透水性、稳定性和阻水性。因此，我们考虑在设计新型护岸框架群时，应结合混泥土四面六边透水框架群在护岸技术上的优势，在混泥土四面六边透水框架群的基础上优化设计，发挥四面六边透水框架群的优势，并尽量弥补四面六边透水框架群在结构强度、易钩连程度、翻滚移位程度上的不足，并综合考虑设计后的框架结构在架空程度、经济生产成本、施工的难易程度等指标，通过机理分析，确定出参数关系，从而设计出四面六边带触脚框架模型（模型一）、六面九边带触脚框架模型（模型二）和双四面六边透水框架群（模型三）然后，我们利用软件2，建立框架群层次分析模型3（模型四）通过建立目标层、决策层和方案层，可以选取施工时架空率接近4-6的程度、结构强度、易翻滚程度、易钩连程度、生产成本、施工简易度六个指标对模型一、模型二、模型三所设计的改价护岸框架和四面六边透水框架群原型进行综合分析评价，以确立出最优的新型护岸框架方案。三、 模型假设1. 护岸框架焊接牢固。2. 护岸框架材料均匀，规格一致。3. 设计的各类框架选材和四面六边透水框架一样。4. 设计的各类框架在施工过程中，不出现偷工减料等现象。5. 不考虑不同市场间的材料的差价。四、 符号说明符号含义单个四面六边透水框架的尺寸（边长与边的横截面的边长之比）单个四面六边透水框架的边长四面六边透水框架每一边的横截面（正方形）的边长护岸框架的触脚长成对比较矩阵（准则层对目标层）成对比较矩阵（方案层会目标层）矩阵中的元素矩阵的归一化特征向量矩阵的阶数矩阵的最大特征值一致性指标随机一致性指标一致性比率矩阵中的元素五、 模型的建立及求解5.1 模型一（四面六边带触脚框架模型）本模型通过机理分析，在四面六边框架（如图5.1）的基础上分析设计。5.1.1 现有四面六边框架分析图5.1 四面六边框架四面六边框架群作为近10年来新型的江河护岸工程技术，具有以下优缺点（如表5.1所示）表5.1 四面六边框架群优缺点优点本身可以近似的看作正三棱锥每一个面都为正三角形，根据物理学三角形具有固定性，四面六边框架群本身就具备良好的稳定性。框架群合理架空率取值分配，拥有恰当的空隙结构区别于传统的护岸框架结构，基础不易被冲刷、适合地形变化，具有良好的透水性和阻水性缺点结构强度的不足虽然原则上能够互相钩连，但依然不清楚最终堆砌而成的形状是否合理被水流冲击而翻滚移位在抛投时和在使用过程中，可能被水流冲击而翻滚移位，使框架群不能达到理想的堆砌效果，对功能有不利影响。针对四面六边框架群在实际护岸中存在的问题和不足，我们为了保证四面六边框架群区别于传统护岸框架的优点，在四面六边框架群的框架结构上进行优化设计。5.1.2 第一种优化设计方案四面六边带触脚框架立体效果图（如图5.2）图5.2 四面六边带触脚框架立体效果图四面六边带触脚框架平面结构构造（如图5.3所示）图5.3 四面六边带触脚框架正视图、左视图、俯视图5.1.3 四面六边带触脚框架优点分析四面六边带触脚框架和四面六边框架相比，具有以下优点（表5.2所示）表5.2 四面六边带触脚框架优点分析表1不易翻滚：四面六边带触脚框架和四面六边框架相比，在每个边上新增了触脚，能有效地阻止水流冲击而导致的翻滚移位；护岸框架一般设置在浅水区域，触脚设计能更好地运用地形优势，加大触脚与淤泥的结合度，在一定程度上也减缓了易翻滚度。2框架群之间更容易钩连：四面六边带触脚框架的触脚设计，容易与框架群自然地相互钩连，这种自然的相互钩连方式相比而言，更加合理。5.1.4 四面六边带触脚框架设计步骤方法步骤一由于带触脚的四面六边框架，相对于四面六边框架有表5.2的优势，所以，我们在四面六边框架的基础上，沿着一边延伸的棱长。（如图5.4所示）图5.4 四面六边透水框架延长一边棱长步骤二确立参数尺寸。如图5.5为图5.4延长的触脚横截面。图5.5 延长的棱长横截面如图5.5所示，为框架每一边的横截面（正方形）的边长，通过机理分析文献资料和实验数据1，我们知道尺寸并不是越大或者越小就越好，的最优取值为15至18。为了方便计算，我们取中间值16，使：现实中，触脚并不是越长越好，实际中，触脚越长越容易折断，触脚越短，效果也就不明显，为了方便计算，我们约定：触脚长为。则：即触脚长为边长的1/5。步骤三按照步骤一到步骤二的方法，同时延长其他几个边，到达如图5.2所示的四面六边带触脚框架。四面六边带触脚框架在四面六边框架的基础上，新增了触脚，使框架结构不尽具备四面六边框架的基本优点，而且不易翻滚，容易钩连，克服了四面六边框架本身的不易钩连。5.2 模型二（六面九边带触脚框架模型）模型二(六面九边框架带触脚框架模型)是模型一的改进模型。5.2.1 第二种优化设计方案：六面九边框架优化设计方案模型二的立体效果图（如图5.6所示）图5.6 六面九边带触脚框架立体图六面九边框架设计平面结构构造（如图5.7）图5.7 六面九边框架带触脚正视图、左视图、俯视图模型二由上下两个四面六边框架构成，形成六面九边框架，每个面都是三角形的六面九边框架，具有更强的结构强度，稳定性虽然相比四面六边框架有所降低，但是与传统的护岸框架稳定性具有一定的优势。综合模型一带触脚的特性，使模型二组成的框架群更加容易相互钩连，触脚容易扎入泥土，使框架群更加稳定，不易翻滚。5.2.2 六面九边带触脚框架在水流中的示例图图5.8 六面九边透水框架旋转如图5.8所示的框架结构，投到岸边，触脚插入泥土中的效果如图5.9所示，图5.9 六面九边框架触脚插入泥沙示意图此时为框架放在岸边无水状态下的示意图，当有水流时，框架如图5.10所示。图5.10 六面九边框架放入水流中示意图5.2.3 六面九边带触脚框架的优点分析六面九边带触脚框架和模型一相比，拥有以下优点（如表5.3所示）表5.3 六面九边带触脚框架优势分析表1更加不易翻滚：如5.2.2六面九边带触脚框架在水流中的示意图（图5.10）六面九边带触脚框架在结构上，由上下两个四面六边框架构成，而且每条边延伸的触脚长，容易与框架群自然地相互钩连，触脚设计也更容易扎深泥沙层，使框架结构更加不易翻滚2框架群之间更容易钩连：模型二的触脚设计，容易使框架群自热相互钩连，而不是为了钩连而钩连，真正达到科学、合理，而且相比模型一而言，模型二触脚更多，钩连效果更为明显。3施工简单，由两个四面六边透水框架上下组成，更容易构造。4节约经济成本，由两个四面六边透水框架上下组成的结构，使六面九边框架在成本计算上，每一个比同等条件下两个四面六边透水框架节约3根混合泥土的造价成本。5.2.4 六面九边带触脚框架设计步骤步骤一参数确定设置。如图5.1所示，我们在四面六边透水框架的基础上设计，取每一个框架的一根边的的边长为，横截面边长为，每一跟延伸的触脚长为，参数计算方法同模型一的参数计算。步骤二在图5.1的基础上，焊接如图5.11所示的六面九边框架结构，由上下两个四面六边透水框架构成。图5.11 无触脚的六面九边框架步骤三如模型一的步骤一到步骤三，在框架结构的边上延伸触脚长，形成最终效果图（如图5.6所示）。5.3 模型三（双四面六边护岸框架模型）模型三的设计旨在不大范围改变四面六边透水框架的整体基本特性的基础上，增加四面六边透水框架结构强度和稳定性，设计内外双层四面六边透水框架体，以达到更好的防冲促淤的效果。5.3.1 第三种优化方案：双四面六边护岸框架模型优化设计结构方案模型三的立体效果图（如图5.12所示）图5.12 双四面六边护岸框架立体效果图双四面六边护岸结构内外层框架结构构造模型三外部框架结构为（如图5.13）图5.13 双四面六边护岸外层框架结构如图5.13所示，双四面六边护岸外层框架中，底部新增一条辅助轴，以达到固定内部框架的目的。模型三内部框架结构为（如图5.14所示），图5.14 双四面六边护岸内层框架结构双四面六边护岸内外层整体构造（如图5.15）图5.15 双四面六边护岸内外层框架结构构造双四面六边护岸框架设计平面结构构造（如图5.16所示），图5.16 双四面六边护岸框架正视图、左视图、俯视图5.3.2 双四面六边护岸框架模型优点分析双四面六边护岸框架和四面六边框架相比，拥有以下优点（如表5.4所示）表5.4 双四面六边护岸框架优势分析表1结构强度：双四面六边护岸框架在结构上由内外两层四面六边框架构成，结构复杂，而且结构强度比模型一、模型二和四面六边框架更好2稳定性强：结构设计上，由于都是四面六边框架，保留了四面六边框架本身在透水性、阻水性上的优势，比模型一稳定性更强3双层结构设计，由两个四面六边透水框架内外两层组成，很大程度上相互保护，外层四面六边透水框架损坏，也还有内层结构4独特框架设计，内外双层框架结构同时又把整体框架划分为小型四面六边透水框架，根据三角形的稳定性，小型的框架结构在很大程度上提高了结构本身稳定性。5.3.3 双四面六边护岸框架设计步骤方法步骤一参数确定设置。外框架的边长设置仍然为，横截面同样边长为，外框架层辅助轴长度为，如图5.12所示内框架层棱长为，内框架层横截面边长为。步骤二制作框架时，有意识的在外框架层内制作如图所示的内框架层结构，内框架层结构和外框架结构方向刚好相反，使总体框架形成内外双层结构，并形成4个小型四面六边框架。步骤三在外框架底层焊接长度的辅助轴，以达到更好的固定内框架层的目的，效果如图5.12所示。5.4 四面六边透水框架框架原型模型为了方便比较设计后的护岸框架模型，我们选取四面六边透水框架，作为原型与之比较，效果图如图5.1所示。四面六边透水框架设计平面结构构造（如图5.17所示），图5.17 四面六边透水框架正视图、左视图、俯视图原型的尺寸为：每边长度为，横截面边长为。原型的优缺点如表5.1所示。5.5 模型四（层次分析法确定优化设计方案选择）5.5.1问题分析四种护岸框架模型，哪种护岸框架效果最好，我们这里采用在原有方案和以上设计的三种方案基础上，建立框架群层次分析模型3的决策方法，首先，可以选取施工时架空率接近4至6的程度、结构强度、易翻滚程度、易钩连程度、生产成本、施工简易度六个指标对所设计的三种改进护岸框架和四面六边透水框架原型进行综合分析评价，确立出最优的新型护岸框架模型方案。5.5.2建立层次结构模型目标是选择一个最优护岸框架，选取准则为施工时架空率易接近4至6、结构强度、不易翻滚程度、框架群间易钩连程度、生产成本及易生产、施工简易度，方案层为原型模型、模型一、模型二、模型三4种方案。准则说明：(1)施工时架空率易接近4至6:架空率表示框架群内部空隙的相对大小。我们根据文献资料和大量的实验数据1得知，设计的护岸框架结构架空率在4至6的时候，阻水性和亲水性效果最佳，到达最好的减速促淤效果。(2)结构强度：指在河岸上的护岸框架自身的稳定性，结构强度越大则自身的稳定越强，构造也越复杂，结构的强度同时也是四面六边透水框架在实际运用中存在，需要克服的问题之一。(3)不易翻滚程度：指护岸框架在实际投入运用过程中，受水流冲击而不易翻滚移位的程度，一个好的护岸框架设计应尽量避免框架结构受水流而翻滚移位，所以不易翻滚程度同样也是评价设计好坏的重要指标因素之一。(4)框架群间易钩连程度：指护岸框架与其他框架的相互钩连程度，护岸框架结构越容易与其他的框架钩连，就易形成框架群，在实际投入使用过程中，也就越稳定，越能形成良好的减速促淤效果。(5)生产成本及易生产：在实际生活投入使用过程中，同样经济生产成本也是重要指标之一，结构强度、稳定性、阻水性、透水性再好，如果生产成本太大，现实中也是不会投入使用的。(6)施工简易度：即生产施工框架的便易程度，越容易生产的框架结构，就越容易生产出来，生产效率也就越高，但是就不能保证其减速促淤效果以及框架结构本身的阻水性、透水性和稳定性，所以施工简易度虽然也是一项影响因素，但和另外因素相比，比重相对较低。5.5.3 构建层次分析图我们运用AHP对其进行系统的分析，得出如图5.18所示的结构模型图图5.18 护岸框架形状决策问题的层次结构图方案层中的框架A是四面六边透水框架框架原型模型，框架B是模型一（四面六边带触脚框架模型），框架C是模型二（六面九边带触脚框架模型），框架D是模型三（双四面六边护岸框架模型）其中，重要性的标度如表5.5所示表5.5 数值重要性标度表数值标度含义与同样重要比稍微重要比比较重要比十分重要比绝对重要上述相邻判断的中间值5.5.4 构建判断矩阵我们利用两两比较的方法，由各准则对比，得到正反矩阵如下矩阵中，表示施工时架空率易接近4到6比结构强度稍微重要，因为当架空率接近4到6时框架的减速消能效果最好。达到设计的初衷：减速消能。同理，施工时架空率易接近4到6因素比框架群的易钩连程度稍微重要；比生产成本及易生产程度明显比较重要，与施工的简易相比，比较重要，所以根据判断依据，。如果护岸框架容易翻滚，那么就算施工时，架空率接近4到6，随着护岸框架的翻滚，架空率也会很快变化，所以不易翻滚程度与施工时架空率易接近4到6同等重要。所以。矩阵中表示结构强度与不易翻滚程度和框架群易钩连程度同等重要，因为这些因素都是影响框架消能减速效果的重要因素；结构强度比生产成本及易生产、施工简易更为重要，因为我们的目的是设计出一个较之有效的护岸框架，如果说生产成本高或者有些难施工也可以接受的。所以。矩阵中表示不易翻滚程度比框架群易钩连程度微小重要，因为如果框架群易翻滚，无论钩连程度好坏也会散开；不易翻滚程度比生产成本和施工简易更为重要，因为我们需要的是能够到达减速促淤效果的护岸框架。所以，。同理，框架群的易钩连程度与生产成本及易生产和施工简易相比是比较重要。所以，。而如果生产成本低而且容易生产，但施工起来很困难也是不可行的，所以生产成本与施工简易相比，后者比前者稍微重要。所以。利用软件2对矩阵A求最大特征值（代码详见附录），得：=6.3303，归一化特征向量为：一致性指标有资料3可查出随机一致性指标的数值表（如表5.6所示），表5.6随机一致性指标的数值123456789000.580.91.121.241.321.411.45由上表知，当时，所以一致性比列即一致性检验通过，上述可作为权向量。在护岸框架形状决策中，已经得到了第2层（准则层）对第1层（目标层）的权向量。用同样的方法构造第3层（方案层）对第2层的每一个准则的成对比较阵，分别如下 这里矩阵中的元素是方案（护岸框架形状）形状与形状对于准则因素的优越性比较尺度。由第3层的陈对比较阵通过Matlab2计算出权向量（代码详见附录），最大特征根，一致性指标和一致性比率结果列入下表（表5.7）。表5.7 护岸框架形状决策问题第3层的计算结果1234560.09100.10930.08030.10510.44760.19960.28710.10930.27530.30780.28290.14040.50000.20900.53040.49440.16360.49500.12200.57250.11400.09270.10590.16504.09614.00424.17424.03254.07104.06060.03200.00140.05800.01080.02370.02020.03560.00150.06450.01200.02630.0225由上表可知，所有矩阵的一致检验均通过。根据以上数据分别计算出框架A、框架B、框架C、框架D在目标中的组合权重列表（如表5.8）表5.8 各框架在目标中的组合权重计算结果0.33310.16970.23120.16710.03920.0596总排序框架A0.09100.10930.08030.10510.44760.19960.1144框架B0.28710.10930.27530.30780.28290.14040.2487框架C0.50000.20900.53040.49440.16360.49500.4432框架D0.12200.57250.11400.09270.10590.16500.1936层次分析法中方案层对目标层的组合一致性比率为其中为第2层对第1层的权向量的第个元素。护岸框架形状决策问题中方案层对目标层的组合一致性比率为所以，组合一致性检验通过检验。结果表明，在护岸框架形状中，框架C（六面九边带触脚框架）的综合护岸减速效果最为优越。六、 结论根据护岸框架所需求的亲水性、稳定性、阻水性、结构强度、互相钩连程度、易翻滚程度和成本综合考量，我们设计出了三种框架模型和原型：模型一：构建四面六边带触脚框架模型，运用触脚设计，很好的结合四面六边透水框架本身的优点特性，在构造上，使框架达到不易翻滚，并与其他的框架自然、合理的相互钩连；模型二：构建六面九边带触脚框架模型，不仅具备良好的亲水性、阻水性和稳定性，而且触脚比模型一更多，使框架更加不易翻滚、框架群之间也更容易钩连；同时，模型二施工简单，更容易构造，也更加节约经济造价成本。模型三：构建双四面六边护岸框架模型，该模型增加护岸框架结构强度和稳定性。运用内外双层结构设计，形成内外双层保障，而且独特框架设计把整体框架划分为4个小型四面六边透水框架，利用三角形的稳定性，该模型结构强度高、稳定性强。四面六边透水框架原型模型：具有稳定性、透水性、阻水性好的特点。三个模型各有优缺点，我们通过模型四（层次分析模型）运用AHP对其进行系统的分析。通过建立判断矩阵，运用一致性检验，通过Matlab软件计算、比较，最后得出结论为：模型二（六面九边带触脚框架模型）最为优越。所以，我们最终选取我们设计的模型二：六面九边带触脚框架模型方案。七、 模型的评价6.1模型的优点（1）本模型运用层次分析法（AHP），通过构建层次结构和判断矩阵，结合占权比重大的指标准则，科学、系统、合理地选择优化方案；（2）模型设计过程中，运用机理分析确立参数的选择；（3）独特触脚设计结构，使框架结构达到不易翻滚的目的，并与其他的框架自然、合理的相互钩连；（4）模型设计运用三角形的稳定性定理；（5）立体水流示意模拟，让框架结构更加形象；（6）运用Matlab软件进行层次分析的计算特征值、特征向量和一致性检验；（7）立体图形结合平面构图，直观、具体地反应出框架设计结构，使人一目了然。6.2模型的缺点模型的缺点层次分析模型（AHP）模型中，准则层的指标选择虽然是本问题中按护岸框架影响因素权重最大的六种因素，但是还是存在一定的主观性。八、 模型的改进与推广7.1模型的改进在建模过程中，因时间仓促，所构建的设计框架模型有限，这样使方案层的选择方式也有限，这些问题留在以后进行改进。7.2模型的推广模型可以推广到其他框架结构的设计；立体图形的构造；物理模型。层次分析法可以推广于人力资源管理模型、人口统计模型以及旅游选择、资源分配等选择性、决策性和最优的问题。九、 参考文献1 张文捷，王玢，麻夏，王南海，江河护岸新技术，北京：中国水利水电出版社，2002.2 赵海滨，MATLAB应用大全，北京：清华大学出版社，2012.3 姜启源,谢金星，数学模型，北京：高等教育出版社，2011.十、 附录计算成对比较的最大特征值、权向量、一致性比率的Matlab代码clear all;close all;a=1 3 1 3 6 5;1/3 1 1 1 4 4;1 1 1 2 4 4;1/3 1 1/2 1 5 5;1/6 1/4 1/4 1/5 1 1/3;1/5 1/4 1/4 1/5 3 1;x,y=eig(a);eigenvalue=diag(y);lamda_a=eigenvalue(1)cia=(lamda_a-6)/5;cra=cia/1.24wa=x(:,1)/sum(x(:,1)b1=1 1/3 1/4 1/2;3 1 1/2 3;4 2 1 5;2 1/3 1/5 1;x,y=eig(b1);eigenvalue=diag(y);lamda_b1=eigenvalue(1)cib1=(lamda_b1-4)/3;crb1=cib1/0.9wb1=x(:,1)/sum(x(:,1)b2=1 1 1/2 1/5;1 1 1/2 1/5;2 2 1 1/3;5 5 3 1;x,y=eig(b2);eigenvalue=diag(y);lamda_b2=eigenvalue(1)cib2=(lamda_b2-4)/3;crb2=cib2/0.9wb2=x(:,1)/sum(x(:,1)b3=1 1/3 1/5 1/2;3 1 1/3 4;5 3 1 4;2 1/4 1/4 1;x,y=eig(b3);eigenvalue=diag(y);lamda_b3=eigenvalue(1)cib3=(lamda_b3-4)/3;crb3=cib3/0.9wb3=x(:,1)/sum(x(:,1)b4=1 1/3 1/4 1;3 1 1/2 4;4 2 1 5;1 1/4 1/5 1;x,y=eig(b4);eigenvalue=diag(y);lamda_b4=eigenvalue(1)cib4=(lamda_b4-4)/3;crb4=cib4/0.9wb4=x(:,1)/sum(x(:,1)b5=1 2 3 3;1/2 1 2 3;1/3 1/2 1 2;1/3 1/3 1/2 1;x,y=eig(b5);eigenvalue=diag(y);lamda_b5=eigenvalue(1)cib5=(lamda_b5-4)/3;crb5=cib5/0.9wb5=x(:,1)/sum(x(:,1)b6=1 2 1/3 1;1/2 1 1/3 1;3 3 1 3;1 1 1/3 1;x,y=eig(b6);eigenvalue=diag(y);lamda_b6=eigenvalue(1)cib6=(lamda_b6-4)/3;crb6=cib6/0.9wb6=x(:,1)/sum(x(:,1)w_sum=wb1,wb2,wb3,wb4,wb5,wb6*waci=cib1,cib2,cib3,cib4,cib5,cib6;cr=ci*wa/sum(0.9*wa)