主從多機可靠性建模仿真 摘要：本文使用Matlab的中的Simulink編制主 ...

1、主從多機可靠性建模仿真摘要：本文使用matlab的中的simulink編制主從多機的模型，進(jìn)行仿真。系統(tǒng)部件的每一個元件使用stateflow建立狀態(tài)轉(zhuǎn)移機模型。然后根據(jù)元件狀態(tài)推定部件狀態(tài)，根據(jù)部件狀態(tài)推定子系統(tǒng)和系統(tǒng)狀態(tài)。此方法簡單有效，直觀易懂，適合復(fù)雜系統(tǒng)的建模和仿真。關(guān)鍵詞：主從多機，可靠性，建模，仿真a simulation for the reliability of a master/slave systemabstract: a simulation of system reliability based on simulink is proposed in this pap

2、er. the state of each element of the system is based on a fsm model. the state of each part is based on the states of its elements. and the state of the whole system is based on its parts. this method has more virtues than usual methods and it can be applied to any complex system.key words: master/s

3、lave system, reliability, model, simulation1 引言matlab提供的simulink是一個用來對動態(tài)系統(tǒng)進(jìn)行建模、仿真和分析的軟件包，它支持連續(xù)、離散及兩者混合的線性和非線性系統(tǒng)。simulink為用戶提供了很方便的圖形化功能模塊以便連接成一個模擬系統(tǒng)，簡化了設(shè)計的流程，減輕了設(shè)計負(fù)擔(dān)。在sinmlink環(huán)境中，利用鼠標(biāo)就可以在模型窗口中直觀地“畫” 出系統(tǒng)模型，然后直接進(jìn)行仿真。1.1 研究對象的物理模型元件部件的基本構(gòu)成單位，也是模型中系統(tǒng)的最小組成單位。 所有元件連續(xù)無故障運行時間 的概率密度分布都遵從負(fù)指數(shù)分布部件由一個或多個元件構(gòu)成。在模型

4、中構(gòu)成一個部件的元件組合方式具有“單純性”，即只能是以下三種情形之一：串聯(lián)組合、并聯(lián)組合、k-out-of-n組合。多個部件可以組合為具有一定功能的子系統(tǒng)。子系統(tǒng)由若干部件構(gòu)成，能完成一組特定功能。為討論方便，模型中整個系統(tǒng)被劃分為5個子系統(tǒng)：通信主機子系統(tǒng)、集線器子系統(tǒng)、三個相互獨立的通信從機子系統(tǒng)。串聯(lián)組合參與組合的所有元（部）件之一失效，則整體失效。并聯(lián)組合參與組合的所有元（部）件全部失效，則整體才失效。k-out-of-n組合參與組合的所有n個元（部）件中，若至少有k個正常，則整體正常；換言之，多于(n-k)個同時失效，則整體才失效。系統(tǒng)的故障系統(tǒng)內(nèi)若存在有處于失效狀態(tài)的部件，則認(rèn)為系

5、統(tǒng)有故障。系統(tǒng)的重大故障系統(tǒng)運行中的下述情形，列為重大故障。這些情形可能同時出現(xiàn)。l l通信主機子系統(tǒng)失效或無法通信（雙機熱備時，雙機同時失效或無法通信）。l l集線器發(fā)生故障。l l因某個通信從機的某類特定故障，引發(fā)通信總線阻塞，整個通信網(wǎng)絡(luò)癱瘓。l l三個通信從機子系統(tǒng)中，有兩個或兩個以上同時失效或無法通信。系統(tǒng)的一般故障除列為重大故障之外的其他故障。（通信）總線阻塞系統(tǒng)中的通信總線是各通信機間的共享信道。在同一時刻，至多只能有一臺通信機的接口電路工作于發(fā)送信息狀態(tài)，并占用信道。當(dāng)某臺通信機出現(xiàn)特定故障，使其通信接口電路無法退出發(fā)送狀態(tài)，會引發(fā)通信總線阻塞。此時，因信道被占用，其他處于正常

6、狀態(tài)的通信機間也無法進(jìn)行通信。雙機熱備份通信主機子系統(tǒng)可配置為雙機熱備份，實現(xiàn)冗余容錯。結(jié)構(gòu)相同的雙機同步運行。正常情況下，一臺為主用機，承擔(dān)完整的信息輸入處理和控制信號輸出職責(zé)；一臺為備用機，僅有信息輸入，控制信號輸出通路被切斷。當(dāng)主用機發(fā)生故障時，主備職責(zé)自動倒換，可保持子系統(tǒng)整體不失效。值得注意的是，由于是熱備份，所以處于備用狀態(tài)的主機也在實際運行，運行中也可能出現(xiàn)故障狀態(tài)，甚至需要人工修復(fù)。圖1-1 基本配置下系統(tǒng)的組成被研究的系統(tǒng)為一個“1主3從”通信系統(tǒng)。在基本配置下，其組成示意圖如圖1-1。系統(tǒng)中有1臺通信主機和3臺通信從機。所有通信主/從機使用共享信道建立物理連接。典型的實例如

7、rs485標(biāo)準(zhǔn)的二線制總線結(jié)構(gòu)。在該通信系統(tǒng)的傳輸協(xié)議中有以下規(guī)定：l 每臺通信機擁有唯一的地址標(biāo)識l 僅主機擁有發(fā)起對話的權(quán)限，主機與從機之間可進(jìn)行“主問從答”式對話l 從機與從機間不進(jìn)行對話l 必須保證在同一時刻，至多只有一臺通信機占用信道發(fā)送信息，其余通信機處于接收監(jiān)聽狀態(tài)每臺通信機內(nèi)部有兩塊主要電路板：控制電路板和接口電路板?？刂齐娐钒迳嫌形⑻幚砥?，運行相應(yīng)的軟件程序。接口電路板負(fù)責(zé)與通信總線聯(lián)系，它的某類故障，會引發(fā)總線阻塞。光電隔離器完成電-光-電信號轉(zhuǎn)換，通過電氣隔離增加安全性。其內(nèi)部采用了提高可靠性為目的的元件冗余設(shè)計。光電隔離器無法阻止接口電路故障引發(fā)的總線阻塞。集線器是一個

8、通信線纜的匯接裝置。12 元件和部件的模型2.1 時鐘信號發(fā)生器生成周期為1個仿真單位，占空比為50%的方波，每一個上升沿表示一天的到來。如圖2-1所示：圖2-1 時鐘信號發(fā)生器2.2 元件子系統(tǒng)模型2.2.1 構(gòu)成控制硬部件、接口部件、集線器部件的元件模型構(gòu)成控制硬部件、接口部件、集線器部件的所有元件特性滿足獨立同分布，對應(yīng)參數(shù)：元件發(fā)生故障的概率為： = 0.00047988利用stateflow建立如圖2-2所示的狀態(tài)轉(zhuǎn)移機模型。為了簡便算法，假定任何元件從故障狀態(tài)恢復(fù)后的第一天，一定處于無故障狀態(tài)，故設(shè)立recover狀態(tài)，表示故障恢復(fù)后第一天的狀態(tài)。圖2-2所示的狀態(tài)轉(zhuǎn)移機模型有一個

9、輸入端口random，用于讀取隨機數(shù)發(fā)生器產(chǎn)生的隨機數(shù)；有一個輸出端口state，用于輸出當(dāng)前元件的狀態(tài)；還有一個上升沿觸發(fā)端口day，用于引入時鐘信號。其外部結(jié)構(gòu)如圖2-3所示。由于simulink自帶的uniform number generater模塊產(chǎn)生的隨機數(shù)與所給的seed有關(guān)，如果seed不變，每次產(chǎn)生的隨機序列都相同。所以隨機數(shù)發(fā)生器采用了matlab函數(shù)rand()，產(chǎn)生0至1之間的平均分布的隨機數(shù)。以便進(jìn)行多次仿真。將圖2-3所示的元件模型封裝為元件子系統(tǒng)模塊，便于部件的搭建。如圖2-4所示。圖2-2 構(gòu)成控制硬部件、接口部件、集線器部件的元件狀態(tài)轉(zhuǎn)移機模型圖2-3 狀態(tài)轉(zhuǎn)

10、移機外部接口 圖2-4 封裝后的元件子系統(tǒng)模塊2.2.2 構(gòu)成控制軟部件的元件模型構(gòu)成控制軟部件的單元件，對應(yīng)分布參數(shù)： 元件發(fā)生故障的概率為： = 0.004788未采取措施1的情況下，構(gòu)成控制軟部件的單元件的狀態(tài)轉(zhuǎn)移機模型與構(gòu)成控制硬部件、接口部件、集線器部件元件狀態(tài)轉(zhuǎn)移機模型大致相同，此處不再贅述。采取措施1的情況下，元件發(fā)生故障后可自動恢復(fù)的概率為：0.95*p2 = 0.0045，需要人工修復(fù)的概率為：0.05*p2 = 0.0002394。利用stateflow建立如圖2-5所示的狀態(tài)轉(zhuǎn)移機模型，并將其與外部模塊一起封裝成類似于圖2-4所示的元件子系統(tǒng)。由于控制軟部件由單元件組成，

11、因此該元件子系統(tǒng)模塊就是控制軟部件模塊。圖2-5構(gòu)成控制軟部件的元件狀態(tài)轉(zhuǎn)移機模型2.2.3 構(gòu)成光電隔離器部件的元件模型構(gòu)成光隔部件的所有元件特性滿足獨立同分布，對應(yīng)參數(shù):光隔部件組成元件的故障可以自行恢復(fù)，故障狀態(tài)只持續(xù)1日。利用stateflow建立如圖2-6所示的狀態(tài)轉(zhuǎn)移機模型。并將其與外部模塊一起封裝成類似于圖2-4所示的光電隔離器元件子系統(tǒng)。圖2-6構(gòu)成光電隔離器部件的元件狀態(tài)轉(zhuǎn)移機模型2.3 系統(tǒng)部件模型2.3.1控制硬部件、接口部件、集線器部件模型控制硬部件、接口部件、集線器部件都可以看做由20個統(tǒng)計特性獨立元件構(gòu)成的部件，并適用串聯(lián)組合。因此，這些部件的狀態(tài)應(yīng)該是構(gòu)成部件的所

12、有元件狀態(tài)的邏輯與運算結(jié)果。圖2-7控制硬部件、接口部件、集線器部件模型2.3.2控制軟部件模型與控制軟部件的元件模型一致。2.3.3光電隔離器部件模型光電隔離器可看作由5個統(tǒng)計特性獨立元件構(gòu)成的部件，其中任意3個元件無故障即能使該部件正常發(fā)揮效能，適用k-out-of-n組合，k=3, n=5。使用simulink中的sum 模塊，將5個元件的輸出狀態(tài)相加，再使用interval test模塊，測試sum模塊的輸出，如果sum模塊的輸出大于2，說明有2個以上的模塊正常工作，interval test模塊輸出1，光電隔離器工作正常；否則interval test模塊輸出0，光電隔離器出現(xiàn)故障。

13、如圖2-8所示。圖2-8 光電隔離器部件模型3 故障的判定和記錄3.1 系統(tǒng)狀態(tài)記錄在一些元件和部件的狀態(tài)輸出端連接一個to workspace模塊，在每一個時鐘脈沖高電平到來的時候記錄一次，順序記錄到數(shù)列中，1表示工作正常，0表示出現(xiàn)故障。仿真結(jié)束后賦值給matlab工作區(qū)內(nèi)的變量，以便統(tǒng)計仿真結(jié)果。3.2、3.3、3.4所述的三種情況中的to workspace模塊互相獨立，互不干擾。例如：系統(tǒng)運行10天，某元件的狀態(tài)數(shù)列為：1 1 0 0 0 1 1 1 1 1，表示該元件在第2天時發(fā)生故障，故障持續(xù)了3天，其余時刻工作正常。采用數(shù)列和矩陣的方式記錄仿真結(jié)果，在使用matlab處理結(jié)果時

14、，帶來了很大的方便。本仿真系統(tǒng)的運行時間為8年，所以數(shù)列的每一行有365*8=2920個元素。3.2元件故障為了統(tǒng)計人工修理的次數(shù)，需要在所有的會出現(xiàn)“需要人工修復(fù)的故障”的元件的狀態(tài)輸出端外連接to workspace模塊。經(jīng)統(tǒng)計，主機子系統(tǒng)不采用“雙機熱備”時，一共需要統(tǒng)計184個元件的狀態(tài)；主機子系統(tǒng)采用“雙機熱備”時，需要統(tǒng)計225個元件的狀態(tài)。下面以主機子系統(tǒng)不采用“雙機熱備”時為例：仿真結(jié)束后，得到184個元件狀態(tài)數(shù)列，每一個數(shù)列有2920個元素，代表每個元件在2920天中的每一天的狀態(tài)。將這些數(shù)列組合起來，得到一個184行、2920列的矩陣a。要計算一次仿真時人工修理的次數(shù)，只要

15、計算矩陣a內(nèi)每一行“連0序列段”的個數(shù)。我們使用的算法是：在矩陣a第一列前加上一列“全1列”，然后計算矩陣a內(nèi)滿足以下條件的元素個數(shù)： 該元素的值為0，但是與它同一行前一個元素的值為1。3.3系統(tǒng)故障系統(tǒng)故障的定義：系統(tǒng)內(nèi)若存在有失效狀態(tài)的部件，則認(rèn)為系統(tǒng)有故障。與3.2所述類似，在所有部件的狀態(tài)輸出端外連接to workspace模塊。經(jīng)統(tǒng)計，主機子系統(tǒng)不采用“雙機熱備”時，一共需要統(tǒng)計15個部件的狀態(tài)；主機子系統(tǒng)采用“雙機熱備”時，需要統(tǒng)計18個部件的狀態(tài)。下面以主機子系統(tǒng)不采用“雙機熱備”時為例：仿真結(jié)束后，得到15個部件狀態(tài)數(shù)列，每一個數(shù)列有2920個元素，代表每個部件在2920天中的

16、每一天的狀態(tài)。將這些數(shù)列組合起來，得到一個15行、2920列的矩陣b。課題要求中規(guī)定：原因不同的兩個或多個故障，如果首尾相連或者重疊發(fā)生，視作同一次故障。因此計算系統(tǒng)故障發(fā)生次數(shù)時，先對矩陣b進(jìn)行處理，把矩陣b每一列的所有元素做邏輯與運算，得到一個1行、2920列的矩陣s。矩陣s中“連0序列段”的個數(shù)，就是一次仿真時系統(tǒng)故障的發(fā)生次數(shù)。3.4系統(tǒng)重大故障根據(jù)“系統(tǒng)重大故障”的定義，在特定的部件狀態(tài)輸出端外連接to workspace模塊。這里的特定部件為：1、 主機子系統(tǒng)的三個部件的狀態(tài)做邏輯與運算后的狀態(tài)輸出。采用“雙機熱備”時，對兩個主機子系統(tǒng)的狀態(tài)做邏輯或運算后輸出。2、 集線器部件的狀

17、態(tài)輸出。3、 在三個通信從機的接口電路旁分別附加一個有10%概率輸出0的模塊，將其輸出與接口電路狀態(tài)輸出做邏輯或運算后連接to workspace模塊。如圖3-1所示。4、 將3個通信子系統(tǒng)的狀態(tài)按照k-out-of-n的形式連接后輸出，當(dāng)k2時，輸出為0。圖3-1 10%概率輸出0模塊不采取措施1、2、3時，需要統(tǒng)計6處狀態(tài)，采取措施2時，需要統(tǒng)計3處狀態(tài)，采取措施1、2、3時，需要統(tǒng)計3處狀態(tài)。計算系統(tǒng)重大故障發(fā)生次數(shù)的算法與3.3中所述的類似。先得到狀態(tài)矩陣c，把矩陣c每一列的所有元素做邏輯與運算，得到一個1行、2920列的矩陣f。矩陣f中“連0序列段”的個數(shù)，就是一次仿真時系統(tǒng)故障的發(fā)

18、生次數(shù)。3.5 無故障運行天數(shù)計算3.3中得到的矩陣s中的“1”的個數(shù)，得到的結(jié)果就是一次仿真時系統(tǒng)的無故障運行天數(shù)。3.6 連續(xù)無故障運行天數(shù)計算3.3中得到的矩陣s中的“連1序列段”的個數(shù)和每段的長度。我們將結(jié)果存入矩陣中。例如：矩陣3 5 7表示有三個“連1序列段”，長度分別為3、5、7，表示系統(tǒng)有三段連續(xù)無故障運行時間，分別是3天、5天和7天。等到100次仿真做完，再匯總計算均值。3.7 連續(xù)無重大故障運行天數(shù)計算方法與3.6中所述類似，對矩陣f進(jìn)行統(tǒng)計即可，此處不再贅述。4 仿真的控制和結(jié)果統(tǒng)計建立5個模型，分別對應(yīng)課題要求上的5種情況，然后編寫5個m文件控制一次仿真的進(jìn)行并返回結(jié)果

19、，最后編寫一個主程序?qū)γ糠N情況循環(huán)執(zhí)行100次并統(tǒng)計結(jié)果，輸出到文本文件。5 仿真結(jié)果分析每個模型仿真運行100次后，得到的結(jié)果如表5-1所示：表5-1 仿真結(jié)果匯總?cè)斯ば蘩淼钠骄螖?shù)系統(tǒng)故障的平均次數(shù)系統(tǒng)重大故障平均次數(shù)平均無故障運行天數(shù)平均連續(xù)無故障運行天數(shù)平均連續(xù)無重大故障運行天數(shù)不采取措施1、2、3308.7200214.7500119.34002124.65009.864721.5077單獨采取措施1257.6800189.2500106.35002241.210011.796524.4923單獨采取措施2305.2200211.8900103.60002134.960010.043

20、624.9534單獨采取措施3375.8500239.520070.02001983.85008.256138.9686同時采取措施1、2、3309.6600213.530044.78002125.73009.922261.2814通過仔細(xì)分析仿真結(jié)果后我們發(fā)現(xiàn)：單獨采取措施1后，“人工修理平均次數(shù)”、“系統(tǒng)故障平均次數(shù)”顯著減少，這是由于采用了看門狗機制，使得絕大部分的軟件故障能夠自己恢復(fù)，減少了“人工修理平均次數(shù)”和“系統(tǒng)故障平均次數(shù)”。單獨采取措施2后，“系統(tǒng)重大故障平均次數(shù)”與不采取措施相比，只是略有減少，這可能是由于由三臺從機引起的“系統(tǒng)重大故障”在“系統(tǒng)重大故障”總數(shù)中所占比例不是

21、很大的關(guān)系。也有可能是樣本數(shù)據(jù)發(fā)生了很大的偏移導(dǎo)致。單從仿真數(shù)據(jù)上來看，單獨采取措施2，對系統(tǒng)可靠性的提升不是很顯著。單獨采取措施3后，“人工修理平均次數(shù)”和“系統(tǒng)故障平均次數(shù)”顯著增加，這是因為增加了一臺通信主機子系統(tǒng)，系統(tǒng)元件數(shù)增加，“人工修理平均次數(shù)”和“系統(tǒng)故障平均次數(shù)”也增加了?！跋到y(tǒng)重大故障平均次數(shù)”與不采取措施相比顯著減少，這再次說明由通信主機引起的“系統(tǒng)重大故障”在“系統(tǒng)重大故障”總數(shù)中所占比例較高。從仿真數(shù)據(jù)上看，采取措施3以后，系統(tǒng)在面對“系統(tǒng)重大故障”時，可靠性有了顯著提高，但這個提高是以系統(tǒng)普通故障數(shù)增加為代價換來的。同時采取措施1、2、3后，“人工修理平均次數(shù)” 和“

22、系統(tǒng)故障平均次數(shù)”與不采取任何措施時相近，這是因為采取措施1的正面效果與采取措施3的副作用相抵消的緣故。“系統(tǒng)重大故障的平均次數(shù)”和“平均連續(xù)無重大故障運行天數(shù)”均達(dá)到了最高。從仿真結(jié)果上看：同時采取措施1、2、3與不采取措施相比，只是顯著提高了系統(tǒng)在面對“系統(tǒng)重大故障”時的可靠性。綜上所述，采取措施1可以提升系統(tǒng)在面對普通故障時的可靠性，而且沒有副作用；采取措施3可以提升系統(tǒng)在面對重大故障時的可靠性，但是會降低系統(tǒng)在面對普通故障時的可靠性。如果措施1、3同時使用，與不采取措施相比，措施1的正面效果被措施3完全抵消，僅能提升系統(tǒng)在面對重大故障時的可靠性。從仿真數(shù)據(jù)來看，措施2對系統(tǒng)可靠性的提升

23、影響不大。6 參考文獻(xiàn)1 “主從多機通信系統(tǒng)可靠性建模研究”課題要求v37 附錄7.1 程序代碼主程序：function main% 這是主程序，控制仿真的進(jìn)行，統(tǒng)計仿真的數(shù)據(jù)times = 100; %仿真進(jìn)行的次數(shù)，也就是數(shù)據(jù)樣本的個數(shù)fid = fopen(result.txt,wt); %存放結(jié)果的文件for t = 1:5%-計算仿真模型的數(shù)據(jù)-temp1 = 0;temp2 = 0;temp3 = 0;temp4 = 0;temp5 = ;temp6 = ;a1 = 0;a2 = 0;a3 = 0;a4 = 0;a5 = 0;a51 = 0;a52 = 0;a6 = 0;a61 =

24、 0;a62 = 0;for n = 1:timesnteval(temp1, temp2, temp3, temp4, temp5, temp6 = dowork,num2str(t)a1 = a1 + temp1;a2 = a2 + temp2;a3 = a3 + temp3;a4 = a4 + temp4;a51 = a51 + sum(temp5);a52 = a52 + length(temp5);a61 = a61 + sum(temp6);a62 = a62 + length(temp6);end;a1 = a1/times; %人工修理的平均次數(shù)a2 = a2/times; %

25、系統(tǒng)故障的平均發(fā)生次數(shù)a3 = a3/times; %系統(tǒng)重大故障的平均發(fā)生次數(shù)a4 = a4/times; %平均無故障運行天數(shù)a5 = a51/a52; %平均連續(xù)無故障運行天數(shù)a6 = a61/a62; %平均連續(xù)無重大故障運行天數(shù)output = a1 a2 a3 a4 a5 a6;fprintf(fid,%6.4f %6.4f %6.4f %6.4f %6.4f %6.4fn,output);%-仿真模型的數(shù)據(jù)計算結(jié)束-end;fclose(fid);運行第一個仿真模型的程序dowork1，其他的類似function count1, count2, count3, count4, c

26、ount5, count6 = dowork1% work1用于運行第一個仿真模型（不采取措施1、2、3）并統(tǒng)計結(jié)果% data1 - data184 是每個元件的狀態(tài)數(shù)列% part1 - part15 是每個部件件的狀態(tài)數(shù)列% fatal1 - fatal5 是系統(tǒng)的重大故障數(shù)列，fatal6是根據(jù)slave1,slave2,slave3推測出來的重大故障數(shù)列% slave1,slave2,slave3是表示子系統(tǒng)工作狀態(tài)的數(shù)列sim(work1,0,2920); %開始仿真，如果時間為8年，那么應(yīng)該這樣設(shè)定:0,2920d1 = ; %初始化d1，d1是所有元件的狀態(tài)矩陣，它的行數(shù)為元件

27、數(shù)，列為每一天d2 = ; %初始化d2，d2是所有部件的狀態(tài)矩陣，表示系統(tǒng)故障。它的行數(shù)為系統(tǒng)的部件個數(shù)，列為每一天d3 = ; %初始化d3，d3是系統(tǒng)的重大故障狀態(tài)矩陣，表示重大系統(tǒng)故障。s = ; %初始化s，s表示整個系統(tǒng)的狀態(tài)矩陣f = ; %初始化f，f表示整個系統(tǒng)的重大故障狀態(tài)矩陣%-以下是計算人工修理次數(shù)的代碼-count1 = 0; %count1是人工修理的次數(shù)for i = 1:184eval(d1 = d1,data,num2str(i),;); %將所有的狀態(tài)數(shù)列整合成矩陣end;d1 = d1;d1(:,1) = 1; %初始狀態(tài)設(shè)定為1m, n = size(d

28、1); %矩陣d1有m行，n列for i = 1:mfor j = 1:nif (d1(i,j) = 0) & (d1(i,j-1) = 0)count1 = count1 + 1;end;end;end;%-以下是計算系統(tǒng)故障數(shù)的代碼-count2 = 0; %系統(tǒng)故障的發(fā)生次數(shù)for i = 1:15eval(d2 = d2,part,num2str(i),;); %將所有的狀態(tài)數(shù)列整合成矩陣end;d2 = d2;d2(:,1) = 1; %初始狀態(tài)設(shè)定為1s = d2(1,:); %將d2的第一行賦值給sm1,n1 = size(d2); %矩陣d2有m1行，n1列for i = 1:

29、m1for j = 1:n1s(1,j) = s(1,j) & d2(i,j); %把s分別與d2的每一行的相應(yīng)元素做邏輯與運算end;end;for j = 1:n1if (s(1,j) = 0) & (s(1,j-1) = 0)count2 = count2 +1;end;end;%-以下是計算系統(tǒng)重大故障發(fā)生次數(shù)的代碼-count3 = 0; %系統(tǒng)重大故障發(fā)生次數(shù)fatal6 = slave1 + slave2 + slave3; %把三個子系統(tǒng)的狀態(tài)相加。m2,n2 = size(fatal6); %矩陣fatal6有m2行，n2列for i = 1:m2if fatal6(i,1) 2fatal6(i,1)