火灾自动报警系统模拟与可靠性解算作品说明书

1、火灾自动报警系统模拟与可靠性解算摘要本作品由以火灾自动报警系统为主题的虚拟公寓模型和可靠性解算程序两部分组成。前者以游戏软件“我的世界”为平台实现模型搭建与情景模拟，具有真实还原度高，可操作性强的特点。后者通过对火灾自动预警系统正常、失效（误报、漏报）状态的模拟与分析，应用系统安全及可靠性理论，确定了火灾自动报警系统应用效能评价指标体系，并根据层次分析法确定了各指标的量化方法，最终建立火灾自动报警系统可靠度计算模型。另外，我小组以Visual Basic为工具，编写了系统各部分对整体可靠性影响的定量解算软件，希望能对今后火灾报警系统可靠性能的改善与提升提供指导意义。为验证该模型的合理性与软件的

2、适用性，在统计分析及专家评判的基础上，以模型为背景，该系统确定了6个目标层指标和3个选择层指标，分别赋予权重并进行解算和一致性检验，解算结果如表1所示。表1 数据输出表探测器报警控制器灭火装置量化结果0.44020.15900.3657因此：探测器对整个系统的可靠性影响最大本作品的创新点在于：1）虚拟模型的外观设计新意；2）利用游戏软件创建模型、进行情景模拟的实用新意；3）火灾自动报警系统可靠度定量解算软件的编写新意。说明书技术领域本作品涉及一种将决策者对复杂系统的决策思维过程模型化、数量化的过程，适用范围较广，不仅可用于火灾自动报警系统可靠性能的定量解算，还可为其他问题最佳方案的选择提供依据

3、。创新背景火灾是无情的，它吞噬生命、毁坏财富，给人们带来的损失是无法估量的，据最新调查，当前全球每年约发生火灾600至700万次，约有6至7万人死于火灾，其中吸入有毒烟气致死的人数占死亡总人数的70%75%。因此，如何在火灾初期阶段及时有效地提醒人员尽快逃离现场对于减少火灾中的人员伤亡具有极大的实用价值和社会效益。随着科学技术发展，火灾自动报警系统发展如火如荼，并取得了较好的防范效果。火灾自动报警系统是由触发器件、火报警装置、火灾警报装置以及具有其他辅助功能装置组成的火灾报警系统。它是为早发现和通报火情并及时采取有效措施，控制和扑灭火灾而设置在建筑中或其他场所的一种自动消防设施。火灾初期，系统

4、里各种自动探测装置，如感温探测器、感烟探测器等，都及时地给出火灾报警电信号，传送到火灾报警控制器开启声光报警来警示相关工作人员，同时控制主机上还显示出火灾发生的部位及时问等信息，工作人员确认后便可以进行联动控制以及其他辅助功能的控制。倒如启动消火栓泵、喷淋泵、防火卷帘等消防用的灭火、减灾设备，及时扑灭火灾，防止火势蔓延，最大限度地减少生命和财产损失。可见，火灾自动报警系统是同火灾做斗争的有力工具，在消防安全保卫工作中发挥着重要作用。但是，由于火灾自动报警系统长期处于备战状态，而且其探测装置对环境的变化及其他干扰又十分敏感，所以导致系统的误报率较高。为了明晰系统的运行情况，保障系统的正常运行，合

5、理对系统进行维护保养并提高其性能，这就要求对实际应用中的火灾自动报警系统可靠性进行分析。由此可见，系统的可靠性分析具有重大意义。创新内容本作品要解决的问题，是采用一种方法使模型操作简单方便，模拟效果直观明朗，并编写一个软件用以解算火灾预警系统各组成部分的可靠度，呈现客观、准确、实时、可靠的定量分析结果。为解决上述问题，本作品采用的方法步骤如下：1、虚拟模型部分1）利用“我的世界”游戏软件。“我的世界”是一款沙盒游戏，以让每一个玩家在三维空间中自由地创造和破坏不同种类的方块为主体。这款游戏最大的特色在于其颇具自由性，使之更加适合于用作模型的搭建与情景的模拟；2）本作品中，社区及公寓部分均由方块组

6、合搭建，并通过安装一些特殊模组实现了火灾预警系统模拟模型的高真实度、高还原度。其模拟性之强，模拟效果之好，可操作性之高，使界面颇具身临其境之感；3）分别对火灾自动预警系统正常、失效（误报、漏报）状态的模拟与分析。2、解算软件部分1）通过应用系统安全及可靠性理论，运用数学建模思想，确定火灾自动报警系统各部分可靠度计算模型；2）根据层次分析法确定了各指标的量化方法，最终建立火灾自动报警系统可靠度计算模型；3）以Visual Basic为工具，编写了系统各部分对整体可靠性影响的定量解算软件；4）在统计分析及专家评判的基础上，以模型为背景，确定了6个目标层指标和3个选择层指标，分别赋予权重进行解算和一

7、致性检验。具体实施方式1、虚拟模型部分1）准备过程：下载“我的世界”游戏软件和模型创建所需要的各个模组；2）创建过程：通过组合搭建各种方块实现社区及公寓的结构设计，另外。通过特殊模组中的元素对模型加以修饰与完善；3）情景模拟与录制：情景一：火灾发生时，火灾自动报警系统正常工作；情景二：火灾发生时，火灾自动报警系统延迟报警；视频录制采用KK录像机。2、解算软件部分1）理论准备过程（1）建模方法选择通过对本问题的分析与查找文献，最后确定采用层次分析法对该问题进行建模，以下对其进行简单介绍。层次分析法，简称AHP方法，是一种定性与定量相结合的决策分析方法。它是一种将决策者对复杂系统的决策思维过程模型

8、化、数量化的过程。应用这种方法，决策者通过将复杂问题分解为若干层次和若干因素，在各因素之间进行简单的比较和计算，就可以得出不同方案的权重，为最佳方案的选择提供依据。层次分析法（AHP）基本原理：AHP法首先把问题层次化，按问题性质和总目标将此问题分解成不同层次，构成一个多层次的分析结构模型，分为最低层（供决策的方案、措施等），相对于最高层（总目标）的相对重要性权值的确定或相对优劣次序的排序问题。（2）分析火灾自动预警及响应系统各部分组成火灾自动报警系统一般由触发器件、火灾报警装置、火灾警报装置和电源组成，复杂系统还包括消防控制设备。可将系统功能组件归类和简化为四类单元：a.触发单元(T)，包括

9、手动报警按钮、火灾探测器、状态输入接口等，作为输入信号源；b.处理与控制单元(C)，主要为火灾报警控制器 (含区域报警控制器)； c.警报与联动控制单元(A)，包括警报装置、火灾应急照明、火警电话、应急广播、联动控制接口等；d.供电单元(P)，包括双电源互投箱、交直流转 化单元、备用电源；（3）确定指标体系选择层（方案层）：探测器、报警控制器、灭火装置目标层：灵敏度、误动作率（误报率）、动作有效性、环境稳定性、可维修性、寿命时长如下图所示：图1 指标体系结构图（4）建立判断矩阵在层次分析法中，为了使判断定量化，关键在于设法使任意两个方案对于某一准则的相对优越程度得到定量描述。一般对单一准则来说

10、，两个方案进行比较总能判断出优劣，层次分析法采用1-9标度方法，对不同情况的评比给出数量标度。如表2所示。表2 标度表标度定义与说明1两个元素对某个属性具有同样重要性3两个元素比较，一元素比另一元素稍微重要5两个元素比较，一元素比另一元素明显重要7两个元素比较，一元素比另一元素重要得多9两个元素比较，一元素比另一元素极端重要2、4、6、8表示需要在上述两个标准之间拆衷时的标度1/bij两个元素的反比较（5）用和积法计算最大特征向量将判断矩阵的每一列元素作归一化处理，其元素的一般项为：将每一列经归一化处理后的判断矩阵按行相加为：对向量W=（W1，W2Wn)归一化处理：W=（W1，W2Wn)即为所

11、求的特征向量的近似解。（6）计算判断矩阵最大特征根（7）计算判断矩阵一致性指标（8）一致性检验当C.R0.10时，便认为判断矩阵具有可以接受的一致性。当C.R0.10时，就需要调整和修正判断矩阵，使其满足C.R.= 0.1 Then MsgBox (请调整和修正判断矩阵) Else For i = 1 To n w0(i) = Format(w0(i), 0.000) Next i End If 选择层部分Dim K As IntegerDim h() As DoubleReDim h(1 To n, 1 To m, 1 To m) For K = 1 To n For i = 1 To m

12、For j = 1 To m h(K, i, j) = xlssheet2.Cells(4 * (K - 1) + i + 1, j + 1) Next j Next i Next K 输入各个H矩阵 Dim hsum() As Double ReDim hsum(1 To n, 1 To m) For K = 1 To n h矩阵个数 For j = 1 To m hsum(K, j) = 0 For i = 1 To m hsum(K, j) = h(K, i, j) + hsum(K, j) Next i Next j Next K Dim h1() As Double ReDim h1

13、(1 To n, 1 To m, 1 To m) For K = 1 To n For j = 1 To m For i = 1 To m h1(K, i, j) = h(K, i, j) / hsum(K, j) Next i Next j Next K 归一化后的h1矩阵 Dim hsall() As Double ReDim hsall(1 To n) Dim hsum2() As Double ReDim hsum2(1 To n, 1 To m) For K = 1 To n hsall(K) = 0 For i = 1 To m hsum2(K, i) = 0 For j = 1

14、To m hsum2(K, i) = h1(K, i, j) + hsum2(K, i) Next j hsall(K) = hsum2(K, i) + hsall(K) Next i Next K Dim hw() As Double ReDim hw(1 To n, 1 To m) For K = 1 To n For i = 1 To m hw(K, i) = hsum2(K, i) / hsall(K) Next i Next K h矩阵的最大特征向量求出 Dim hjzj() As Double ReDim hjzj(1 To n, 1 To m) hjzj指h矩阵的矩阵积 For

15、K = 1 To n For i = 1 To m hjzj(K, i) = 0 For j = 1 To m hjzj(K, i) = hjzj(K, i) + h1(K, i, j) * hw(K, j) Next j Next i Next K 矩阵的积 Dim hf() As Double ReDim hf(1 To n) Dim hfci() As Double ReDim hfci(1 To n) Dim hfcr() As Double ReDim hfcr(1 To n) For K = 1 To n hf(K) = 0 For i = 1 To m hf(K) = hf(K)

16、 + hjzj(K, i) / (n * hw(K, i) Next i hfci(K) = (hf(K) - n) / (n - 1) If n = 15 Then ri = 1.59 Else If n = 14 Then ri = 1.58 Else If n = 13 Then ri = 1.56 Else If n = 12 Then ri = 1.54 Else If n = 11 Then ri = 1.52 Else If n = 10 Then ri = 1.49 Else If n = 9 Then ri = 1.46 Else If n = 8 Then ri = 1.4

17、1 Else If n = 7 Then ri = 1.32 Else If n = 6 Then ri = 1.24 Else If n = 5 Then ri = 1.12 Else If n = 4 Then ri = 0.9 Else If n = 3 Then ri = 0.58 Else If n = 2 Or n = 1 Then ri = 0 Else End If hfcr(K) = hfci(K) / ri Next K 计算矩阵最大特征根、随机一致性指标、随机一致性比率 For K = 1 To n If hfcr(K) = 0.1 Then MsgBox (请调整和修正

18、判断矩阵) Else For i = 1 To m hw(K, i) = Format(hw(K, i), 0.000) Next i End If Next K 计算各个目标得分Dim zf() As Double zf表总分ReDim zf(1 To m) For j = 1 To m zf(j) = 0 For i = 1 To n zf(j) = zf(j) + w0(i) * hw(i, j) Next i zf(j) = Format(zf(j), 0.0000) Next j 求最大值并记录位置Dim p, q As IntegerMax = zf(1)q = 1 For p =

19、 2 To m If zf(p) Max Then Max = zf(p) q = p End If Next p Picture6.Print 最高得分为第; q; 个目标，得分为; Max xlbook.Close (True)xlbook2.Close (True)Set xlssheet = NothingSet xlbook = NothingSet VBExcel = NothingSet xlssheet2 = NothingSet xlbook2 = NothingSet VBExcel = NothingEnd SubPrivate Sub Command2_Click()

20、Dim VBExcel As Excel.Application Dim VBExcel2 As Excel.Application Dim xlbook, xlbook2 As Excel.Workbook Dim xlssheet, xlssheet2 As Excel.Worksheet Set VBExcel = CreateObject(Excel.Application) Set VBExcel2 = CreateObject(Excel.Application) Set xlbook = VBExcel.Workbooks.Open(C:UserszjDesktop创新比赛数据输

21、入表1.xlsx) 指定你的路径 Set xlssheet = xlbook.Worksheets(1) 指定工作表，也可以使用sheet1 Set xlbook2 = VBExcel.Workbooks.Open(C:UserszjDesktop创新比赛数据输入表2.xlsx) 指定你的路径 Set xlssheet2 = xlbook2.Worksheets(1) 指定工作表，也可以使用sheet1 VBExcel.Visible = True m = InputBox(请输入选择层因子个数, 输入数据) n = InputBox(请输入目标层因子个数, 输入数据) Dim p1() Re

22、Dim p1(1 To n, 1 To n) Dim p2() ReDim p2(1 To n, 1 To n) p1为原始数据矩阵 p2为归一化后数据矩阵 For i = 1 To n For j = 1 To n p1(i, j) = xlssheet.Cells(i + 1, j + 1) Next j Next i 输入P1矩阵 Dim sum() As Double ReDim sum(1 To n) For j = 1 To n sum(j) = 0 For i = 1 To n sum(j) = p1(i, j) + sum(j) Next i Next j For j = 1

23、To n For i = 1 To n p2(i, j) = p1(i, j) / sum(j) Next i Next j 归一化后的P2矩阵 Dim sum2() As Double ReDim sum2(1 To n) sall = 0 For i = 1 To n sum2(i) = 0 For j = 1 To n sum2(i) = p2(i, j) + sum2(i) Next j sall = sum2(i) + sall Next i Dim w0() As Single ReDim w0(1 To n) For i = 1 To n w0(i) = sum2(i) / sa

24、ll Next i P矩阵的最大特征向量求出 Dim jzj() As Double ReDim jzj(1 To n) For i = 1 To n jzj(i) = 0 For j = 1 To n jzj(i) = jzj(i) + p1(i, j) * w0(j) Next j Next i 矩阵的积 Dim f As Double f = 0 For i = 1 To n f = f + jzj(i) / (n * w0(i) Next i 计算矩阵最大特征根 fci = (f - n) / (n - 1) If n = 15 Then ri = 1.59 Else If n = 1

25、4 Then ri = 1.58 Else If n = 13 Then ri = 1.56 Else If n = 12 Then ri = 1.54 Else If n = 11 Then ri = 1.52 Else If n = 10 Then ri = 1.49 Else If n = 9 Then ri = 1.46 Else If n = 8 Then ri = 1.41 Else If n = 7 Then ri = 1.32 Else If n = 6 Then ri = 1.24 Else If n = 5 Then ri = 1.12 Else If n = 4 Then

26、 ri = 0.9 Else If n = 3 Then ri = 0.58 Else If n = 2 Or n = 1 Then ri = 0 Else End If fcr = fci / ri 计算随机一致性比率 If fcr = 0.1 Then MsgBox (请调整和修正判断矩阵) Else Picture1.Print 目标层矩阵最大特征根为： For i = 1 To n w0(i) = Format(w0(i), 0.000) Picture1.Print w0(i) Next i End If 选择层部分Dim K As IntegerDim h() As DoubleR

27、eDim h(1 To n, 1 To m, 1 To m) For K = 1 To n For i = 1 To m For j = 1 To m h(K, i, j) = xlssheet2.Cells(4 * (K - 1) + i + 1, j + 1) Next j Next i Next K 输入各个H矩阵 Dim hsum() As Double ReDim hsum(1 To n, 1 To m) For K = 1 To n h矩阵个数 For j = 1 To m hsum(K, j) = 0 For i = 1 To m hsum(K, j) = h(K, i, j)

28、+ hsum(K, j) Next i Next j Next K Dim h1() As Double ReDim h1(1 To n, 1 To m, 1 To m) For K = 1 To n For j = 1 To m For i = 1 To m h1(K, i, j) = h(K, i, j) / hsum(K, j) Next i Next j Next K 归一化后的h1矩阵 Dim hsall() As Double ReDim hsall(1 To n) Dim hsum2() As Double ReDim hsum2(1 To n, 1 To m) For K =

29、1 To n hsall(K) = 0 For i = 1 To m hsum2(K, i) = 0 For j = 1 To m hsum2(K, i) = h1(K, i, j) + hsum2(K, i) Next j hsall(K) = hsum2(K, i) + hsall(K) Next i Next K Dim hw() As Double ReDim hw(1 To n, 1 To m) For K = 1 To n For i = 1 To m hw(K, i) = hsum2(K, i) / hsall(K) Next i Next K h矩阵的最大特征向量求出 Dim

30、hjzj() As Double ReDim hjzj(1 To n, 1 To m) hjzj指h矩阵的矩阵积 For K = 1 To n For i = 1 To m hjzj(K, i) = 0 For j = 1 To m hjzj(K, i) = hjzj(K, i) + h1(K, i, j) * hw(K, j) Next j Next i Next K 矩阵的积 Dim hf() As Double ReDim hf(1 To n) Dim hfci() As Double ReDim hfci(1 To n) Dim hfcr() As Double ReDim hfcr(

31、1 To n) For K = 1 To n hf(K) = 0 For i = 1 To m hf(K) = hf(K) + hjzj(K, i) / (n * hw(K, i) Next i hfci(K) = (hf(K) - n) / (n - 1) If n = 15 Then ri = 1.59 Else If n = 14 Then ri = 1.58 Else If n = 13 Then ri = 1.56 Else If n = 12 Then ri = 1.54 Else If n = 11 Then ri = 1.52 Else If n = 10 Then ri =

32、1.49 Else If n = 9 Then ri = 1.46 Else If n = 8 Then ri = 1.41 Else If n = 7 Then ri = 1.32 Else If n = 6 Then ri = 1.24 Else If n = 5 Then ri = 1.12 Else If n = 4 Then ri = 0.9 Else If n = 3 Then ri = 0.58 Else If n = 2 Or n = 1 Then ri = 0 Else End If hfcr(K) = hfci(K) / ri Next K 计算矩阵最大特征根、随机一致性指

33、标、随机一致性比率 For K = 1 To n Picture2.Print 第; K; 个目标层矩阵的最大特征根为： If hfcr(K) = 0.1 Then MsgBox (请调整和修正判断矩阵) Else For i = 1 To m hw(K, i) = Format(hw(K, i), 0.000) Picture2.Print hw(K, i); Next i End If Picture2.Print Next K 计算各个目标得分Dim zf() As Double zf表总分ReDim zf(1 To m) For j = 1 To m zf(j) = 0 For i =

34、 1 To n zf(j) = zf(j) + w0(i) * hw(i, j) Next i zf(j) = Format(zf(j), 0.0000) Picture3.Print 第; j; 个目标层矩阵总得分为:; zf(j) Next j 比较大小排序用Dim temp As Double temp为临时交换用For i = 1 To m - 1 For j = i + 1 To m If zf(i) ; Next K Picture4.Print zf(K); xlbook.Close (True)xlbook2.Close (True)Set xlssheet = Nothing

35、Set xlbook = NothingSet VBExcel = NothingSet xlssheet2 = NothingSet xlbook2 = NothingSet VBExcel = NothingEnd Sub附件21、验算原始数据：判断矩阵A表示为下：Hp1p2p3p4p5p6p111/21/3112p2211/2213p3321233p411/21/211/21p5111/3211p61/21/31/3111探测器、报警控制器、灭火装置分别对每个目标的层性赋予标度：（1）灵敏度P1：H1ABCA177B1/711C1/711（2）误动作率（误报率）P2：H2ABCA134B1/311/2C1/421（3）动作有效性P3：H1ABCA111/2B111/3C231（4）环境稳定性P4：H1ABCA1