安全系統(tǒng)工程安全預(yù)測(cè)實(shí)用教案

1、16 6 系統(tǒng)安全預(yù)測(cè)系統(tǒng)安全預(yù)測(cè)(yc)(yc)6.1 安全預(yù)測(cè)概述安全預(yù)測(cè)概述6.2 德?tīng)柗祁A(yù)測(cè)法德?tīng)柗祁A(yù)測(cè)法6.3 回歸回歸(hugu)分析法分析法預(yù)測(cè)預(yù)測(cè)6.4 馬爾柯夫鏈預(yù)測(cè)法馬爾柯夫鏈預(yù)測(cè)法6.5 灰色系統(tǒng)預(yù)測(cè)法灰色系統(tǒng)預(yù)測(cè)法第1頁(yè)/共69頁(yè)第一頁(yè)，共70頁(yè)。26.1.1 安全預(yù)測(cè)(yc)的定義第2頁(yè)/共69頁(yè)第二頁(yè)，共70頁(yè)。36.1.1 安全預(yù)測(cè)(yc)的定義第3頁(yè)/共69頁(yè)第三頁(yè)，共70頁(yè)。46.1.2 安全預(yù)測(cè)(yc)的分類（1）按預(yù)測(cè)對(duì)象的范圍分1）宏觀預(yù)測(cè)。是指對(duì)整個(gè)(zhngg)行業(yè)、一個(gè)地區(qū)、一個(gè)企業(yè)的安全狀況的預(yù)測(cè)。2）微觀預(yù)測(cè)。是指對(duì)一個(gè)生產(chǎn)單位的生產(chǎn)系統(tǒng)或其子

2、系統(tǒng)的安全狀況的預(yù)測(cè)。（2）按時(shí)間長(zhǎng)短分1）長(zhǎng)期預(yù)測(cè)。2）中期預(yù)測(cè)。3）短期預(yù)測(cè)。第4頁(yè)/共69頁(yè)第四頁(yè)，共70頁(yè)。5（3）按所應(yīng)用的原理分1）白色理論預(yù)測(cè)。用于預(yù)測(cè)的問(wèn)題與所受影響因素已十分清楚的情況。2）灰色理論預(yù)測(cè)。也稱為灰色系統(tǒng)預(yù)測(cè)，灰色系統(tǒng)既包含有已知信息(xnx)又含有未知信息(xnx)的系統(tǒng)。安全生產(chǎn)活動(dòng)本身就是個(gè)灰色系統(tǒng)。3）黑色理論預(yù)測(cè)。也稱為黑箱系統(tǒng)或黑色系統(tǒng)預(yù)測(cè)。這種系統(tǒng)中所含的信息(xnx)多為非確定的。常用的安全預(yù)測(cè)方法：德?tīng)柗祁A(yù)測(cè)法、回歸分析法預(yù)測(cè)、馬爾柯夫鏈預(yù)測(cè)法、灰色系統(tǒng)預(yù)測(cè)法。6.1.2 安全預(yù)測(cè)(yc)的分類第5頁(yè)/共69頁(yè)第五頁(yè)，共70頁(yè)。66.1.3 安

3、全(nqun)預(yù)測(cè)的步驟*一般說(shuō)來(lái)，預(yù)測(cè)可分為(fn wi)4個(gè)階段10個(gè)步驟。確定預(yù)測(cè)目標(biāo)和任務(wù)1.確定預(yù)測(cè)目的2.制定預(yù)測(cè)計(jì)劃3.確定預(yù)測(cè)時(shí)間輸入信息4.收集預(yù)測(cè)資料5.檢驗(yàn)現(xiàn)有資料預(yù)測(cè)處理6.選擇預(yù)測(cè)方法7.建立預(yù)測(cè)模型8.進(jìn)行推測(cè)或計(jì)算輸出結(jié)果9.預(yù)測(cè)結(jié)果的鑒定10.修正預(yù)測(cè)結(jié)果階段(jidun) 步驟 第6頁(yè)/共69頁(yè)第六頁(yè)，共70頁(yè)。76 6 系統(tǒng)安全預(yù)測(cè)系統(tǒng)安全預(yù)測(cè)(yc)(yc)6.1 安全預(yù)測(cè)概述安全預(yù)測(cè)概述(i sh)6.2 德?tīng)柗祁A(yù)測(cè)法德?tīng)柗祁A(yù)測(cè)法6.3 回歸分析法預(yù)測(cè)回歸分析法預(yù)測(cè)6.4 馬爾柯夫鏈預(yù)測(cè)法馬爾柯夫鏈預(yù)測(cè)法6.5 灰色系統(tǒng)預(yù)測(cè)法灰色系統(tǒng)預(yù)測(cè)法第7頁(yè)/共69

4、頁(yè)第七頁(yè)，共70頁(yè)。86.2 德 爾 菲 預(yù) 測(cè) 法 德?tīng)柗祁A(yù)測(cè)法是美國(guó)蘭德公司于20世紀(jì)40年代首先提出并用于技術(shù)預(yù)測(cè)的一種直觀預(yù)測(cè)法。德?tīng)柗祁A(yù)測(cè)法是一種廣為適用的預(yù)測(cè)方法。它既可用于科技預(yù)測(cè)，也可用于社會(huì)、經(jīng)濟(jì)預(yù)測(cè)；既可用于短期預(yù)測(cè)，也可用于長(zhǎng)期預(yù)測(cè)。德?tīng)柗祁A(yù)測(cè)法的優(yōu)點(diǎn)：能對(duì)大量非技術(shù)性的無(wú)法定量分析的因素作出概率估計(jì)，并將概率估計(jì)結(jié)果(ji gu)告訴專家，充分發(fā)揮信息反饋和信息控制的作用，使分散的評(píng)估意見(jiàn)逐次收斂，最后集中在協(xié)調(diào)一致的評(píng)估結(jié)果(ji gu)上。第8頁(yè)/共69頁(yè)第八頁(yè)，共70頁(yè)。96.2.1 德?tīng)柗祁A(yù)測(cè)法的基本(jbn)程序 德?tīng)柗祁A(yù)測(cè)法的實(shí)質(zhì)是利用專家的知識(shí)、經(jīng)驗(yàn)、智慧

5、等無(wú)法數(shù)量化而帶來(lái)很大模糊性的信息。通過(guò)通信的方式進(jìn)行信息交換，逐步地取得較一致的意見(jiàn)，達(dá)到預(yù)測(cè)的目的。德?tīng)柗祁A(yù)測(cè)法的基本程序?yàn)椋海?）確定預(yù)測(cè)目標(biāo)目標(biāo)選擇應(yīng)是本系統(tǒng)或本專業(yè)中對(duì)發(fā)展規(guī)劃有重大影響而且意見(jiàn)分歧較大的課題，預(yù)測(cè)期限以中、遠(yuǎn)期(yun q)為宜。如工礦企業(yè)傷亡事故發(fā)展趨勢(shì)預(yù)測(cè)。第9頁(yè)/共69頁(yè)第九頁(yè)，共70頁(yè)。10（2）成立管理小組 人數(shù)為二人到十幾人，要求具備(jbi)必要的專業(yè)知識(shí)和統(tǒng)計(jì)學(xué)、數(shù)據(jù)處理等方面的知識(shí)。其任務(wù)為：對(duì)預(yù)測(cè)的工作過(guò)程進(jìn)行設(shè)計(jì)，提出可供選擇的專家名單，搞好專家征詢和輪間信息反饋工作，整理預(yù)測(cè)結(jié)果和寫(xiě)出預(yù)測(cè)報(bào)告書(shū)。（3）選擇專家 專家選擇是預(yù)測(cè)成敗的關(guān)鍵，主要

6、要求為： 1）專家總體的權(quán)威程度較高。 2）專家的代表面應(yīng)廣泛，技術(shù)、管理、情報(bào)6.2.1 德?tīng)柗祁A(yù)測(cè)法的基本(jbn)程序 第10頁(yè)/共69頁(yè)第十頁(yè)，共70頁(yè)。11 專家和高層決策人員。 3）嚴(yán)格專家的推薦和審定程序。 4）專家人數(shù)要適當(dāng)。一般為2050人，大型預(yù)測(cè)可達(dá)100人左右。（4）設(shè)計(jì)評(píng)估意見(jiàn)征詢表 沒(méi)有統(tǒng)一的格式，但要求： 1）表格的每一欄目要緊扣預(yù)測(cè)目標(biāo)。 2）表格簡(jiǎn)明扼要。專家思考時(shí)間長(zhǎng)，填表時(shí)間短，24小時(shí)(xiosh)。 3）填表方式簡(jiǎn)單。用數(shù)字及字母表示結(jié)果。6.2.1 德?tīng)柗祁A(yù)測(cè)法的基本(jbn)程序 第11頁(yè)/共69頁(yè)第十一頁(yè)，共70頁(yè)。12（5）專家征詢和輪間信息反

7、饋一般分34輪征詢。在第一輪征詢中，預(yù)先擬訂一個(gè)預(yù)測(cè)事件一覽表，并允許專家對(duì)表格補(bǔ)充和修改。管理小組應(yīng)將與課題有關(guān)的大量技術(shù)政策和經(jīng)濟(jì)條件等背景材料提供給專家。在征詢表中，最常見(jiàn)的問(wèn)題是要求(yoqi)專家對(duì)某項(xiàng)技術(shù)實(shí)現(xiàn)的日期作出預(yù)言。6.2.1 德?tīng)柗祁A(yù)測(cè)(yc)法的基本程序 第12頁(yè)/共69頁(yè)第十二頁(yè)，共70頁(yè)。13德?tīng)柗祁A(yù)測(cè)(yc)法的特點(diǎn) 德?tīng)柗品ㄊ且粋€(gè)可控制的組織集體思想(sxing)交流的過(guò)程，使得由各方面的專家組成的集體能作為一個(gè)整體來(lái)解答某個(gè)復(fù)雜問(wèn)題。（1）匿名性。專家能充分發(fā)揮自己的智慧和知識(shí)。（2）反饋性。管理小組對(duì)每一輪的預(yù)測(cè)結(jié)果指出統(tǒng)計(jì)、匯總，提供有關(guān)專家的論證依據(jù)和資

8、料，作為反饋材料發(fā)給每一位專家，供下一輪預(yù)測(cè)時(shí)參考。（3）預(yù)測(cè)結(jié)果的統(tǒng)計(jì)特性。 第13頁(yè)/共69頁(yè)第十三頁(yè)，共70頁(yè)。146.2.2 專家(zhunji)意見(jiàn)的統(tǒng)計(jì)處理 （1）數(shù)量和時(shí)間答案的處理常用中位數(shù)和上、下四分點(diǎn)的方法(fngf)，處理專家們的答案，求出預(yù)測(cè)的期望值和時(shí)間。將專家的回答按從小到大的順序排列。如：當(dāng)有n 個(gè)專家時(shí)，共有n個(gè)答數(shù)：o其中(qzhng)位數(shù)按下式計(jì)算 nnxxxx121)(22)(1211偶數(shù)奇數(shù)knxxknxxkkk第14頁(yè)/共69頁(yè)第十四頁(yè)，共70頁(yè)。15（1）數(shù)量和時(shí)間(shjin)答案的處理 上四分(s fn)點(diǎn)按下式計(jì)算 為偶數(shù)，為奇數(shù)，為偶數(shù)，為奇

9、數(shù)，）（）（上kknxxkknxkknxxkknxxkkkkkk22212212123231321223233321第15頁(yè)/共69頁(yè)第十五頁(yè)，共70頁(yè)。16下四分(s fn)點(diǎn)按下式計(jì)算 為偶數(shù)，為奇數(shù)，為偶數(shù)，為奇數(shù)，下kknxxkknxkknxxkknxxkkkkkk222122121222112221（1）數(shù)量和時(shí)間答案(d n)的處理 第16頁(yè)/共69頁(yè)第十六頁(yè)，共70頁(yè)。17數(shù)量和時(shí)間答案(d n)的處理舉例例：某企業(yè)邀請(qǐng)例：某企業(yè)邀請(qǐng)16位專家對(duì)該企業(yè)某事故發(fā)生概率進(jìn)行位專家對(duì)該企業(yè)某事故發(fā)生概率進(jìn)行預(yù)測(cè)預(yù)測(cè)(yc)，得到，得到16個(gè)數(shù)據(jù)，即個(gè)數(shù)據(jù)，即 n=16，n = 2k，k

10、 = 8為偶數(shù)。數(shù)據(jù)由小到大為：為偶數(shù)。數(shù)據(jù)由小到大為： n12345678事故發(fā)生概率事故發(fā)生概率p1031.35 1.38 1.40 1.40 1.40 1.45 1.47 1.50n910111213141516事故發(fā)生概率事故發(fā)生概率p1031.50 1.50 1.50 1.53 1.55 1.60 1.60 1.65解：解： n=16為偶數(shù)為偶數(shù)(u sh)，k = 8，則中位數(shù)為：，則中位數(shù)為：50. 1)(212981xxxxxkk第17頁(yè)/共69頁(yè)第十七頁(yè)，共70頁(yè)。18數(shù)量和時(shí)間答案的處理(chl)舉例因此，處理后的預(yù)測(cè)結(jié)果為：事故因此，處理后的預(yù)測(cè)結(jié)果為：事故(shg)發(fā)生

11、概率的發(fā)生概率的期望值為期望值為 p = 1.5010-3，上限，上限p上上 = 1.5410-3，下限下限p下下 = 1.4010-3。 k = 8為偶數(shù)(u sh)，則上中位數(shù)為：54. 1212131212323）（上x(chóng)xxxxkk k = 8為偶數(shù)，則下中位數(shù)為：40. 121254122）（下xxxxxkk第18頁(yè)/共69頁(yè)第十八頁(yè)，共70頁(yè)。19（2）等級(jí)比較答案的處理在邀請(qǐng)專家進(jìn)行安全預(yù)測(cè)時(shí)，常有對(duì)某些項(xiàng)目的重要性進(jìn)行排序的要求，如預(yù)防措施的選擇、事故(shg)原因的確定等。對(duì)這類問(wèn)題，可采用評(píng)分法處理。當(dāng)對(duì)n 個(gè)項(xiàng)目排序時(shí)，第1位得 n 分，第2位得 n-1 分，然后計(jì)算：mj

12、NBsniiij, 2 , 1,16.2.2 專家(zhunji)意見(jiàn)的統(tǒng)計(jì)處理 第19頁(yè)/共69頁(yè)第十九頁(yè)，共70頁(yè)。20（2）等級(jí)比較答案(d n)的處理 式中，m 參加比較的項(xiàng)目個(gè)數(shù)； sj 第j個(gè)項(xiàng)目的總得分； kj 第j個(gè)項(xiàng)目的得分比重； n 要求排序的項(xiàng)目個(gè)數(shù)； Bi 排在第i位項(xiàng)目的得分； Ni 將某項(xiàng)目排在第i位的專家(zhunji)人數(shù)； M對(duì)問(wèn)題作出回答的專家(zhunji)總?cè)藬?shù)。nijjiMsk1第20頁(yè)/共69頁(yè)第二十頁(yè)，共70頁(yè)。21等級(jí)(dngj)比較答案的處理舉例 例：某礦井例：某礦井(kungjng)下發(fā)生了火災(zāi)，大量煤炭正在下發(fā)生了火災(zāi)，大量煤炭正在燃燒，不

13、僅造成大量經(jīng)濟(jì)損失，而且對(duì)礦井燃燒，不僅造成大量經(jīng)濟(jì)損失，而且對(duì)礦井(kungjng)安全生產(chǎn)也構(gòu)成了威脅。為消滅火災(zāi)，安全生產(chǎn)也構(gòu)成了威脅。為消滅火災(zāi)，共提出了六個(gè)方案，請(qǐng)共提出了六個(gè)方案，請(qǐng)93位專家對(duì)從中選出三個(gè)方位專家對(duì)從中選出三個(gè)方案并對(duì)其排序。案并對(duì)其排序。 方案序號(hào)方案序號(hào)方方 案案 內(nèi)內(nèi) 容容123456 密閉火區(qū)，利用風(fēng)壓平衡法控制漏風(fēng)密閉火區(qū)，利用風(fēng)壓平衡法控制漏風(fēng) 密閉火區(qū)，向火區(qū)內(nèi)注水密閉火區(qū)，向火區(qū)內(nèi)注水 密閉火區(qū)，向火區(qū)內(nèi)注泥漿密閉火區(qū)，向火區(qū)內(nèi)注泥漿 密閉火區(qū)，向火區(qū)內(nèi)注爐煙密閉火區(qū)，向火區(qū)內(nèi)注爐煙 密閉火區(qū)，向火區(qū)內(nèi)注液氮密閉火區(qū)，向火區(qū)內(nèi)注液氮 密閉火區(qū)，利用

14、風(fēng)壓平衡法控制漏風(fēng)并向密閉火區(qū)，利用風(fēng)壓平衡法控制漏風(fēng)并向火區(qū)內(nèi)注泥漿火區(qū)內(nèi)注泥漿第21頁(yè)/共69頁(yè)第二十一頁(yè)，共70頁(yè)。22解：要求解：要求(yoqi)從六個(gè)方案中選三個(gè)，故從六個(gè)方案中選三個(gè)，故n = 3，排在第，排在第1、2、3位的得分為：位的得分為：B1 = 3， B2 = 2， B3 = 1。將第。將第1方案排在第方案排在第1、2、3位的專家人數(shù)分別為位的專家人數(shù)分別為71、15、2，則則N1 = 71， N2 = 15， N3 = 2。于是得：。于是得：24521152713311iiiNBs44. 0) 321 (93245111niiMsk等級(jí)(dngj)比較答案的處理舉例 第

15、22頁(yè)/共69頁(yè)第二十二頁(yè)，共70頁(yè)。23o用同樣方法處理其他(qt)項(xiàng)目，所得結(jié)果如下表。方案序號(hào)方案序號(hào)123456各方案得分各方案得分2453665531168各方案得分比重各方案得分比重0.440.070.120.010.060.30各方案排序各方案排序o通過(guò)比較(bjio)各方案的得分比重，認(rèn)為應(yīng)采取的措施及其順序?yàn)椋旱?、第6、第3方案，即都需要密閉火區(qū)，綜合采取風(fēng)壓平衡法控制漏風(fēng)并向火區(qū)內(nèi)注泥漿。等級(jí)比較(bjio)答案的處理舉例 第23頁(yè)/共69頁(yè)第二十三頁(yè)，共70頁(yè)。246 6 系統(tǒng)安全預(yù)測(cè)系統(tǒng)安全預(yù)測(cè)(yc)(yc)6.1 安全預(yù)測(cè)安全預(yù)測(cè)(yc)概述概述6.2 德?tīng)柗祁A(yù)測(cè)

16、德?tīng)柗祁A(yù)測(cè)(yc)法法6.3 回歸分析法預(yù)測(cè)回歸分析法預(yù)測(cè)(yc)6.4 馬爾柯夫鏈預(yù)測(cè)馬爾柯夫鏈預(yù)測(cè)(yc)法法6.5 灰色系統(tǒng)預(yù)測(cè)灰色系統(tǒng)預(yù)測(cè)(yc)法法第24頁(yè)/共69頁(yè)第二十四頁(yè)，共70頁(yè)。256.3 回歸(hugu)分析法預(yù)測(cè)&企業(yè)或部門(mén)的安全狀況(zhungkung)與影響它的各種因素是一個(gè)密切聯(lián)系的整體。&而這個(gè)整體又具有相對(duì)穩(wěn)定性和持續(xù)性，即時(shí)間序列平穩(wěn)性。這樣就可以拋開(kāi)對(duì)逐個(gè)因素的分析，就其整體利用慣性原理，對(duì)企業(yè)或部門(mén)的安全狀況(zhungkung)進(jìn)行預(yù)測(cè)提供了可能。 第25頁(yè)/共69頁(yè)第二十五頁(yè)，共70頁(yè)。266.3.1 回歸(hugu)分析法概述&am

17、p;企業(yè)或部門(mén)的安全狀況可以用一定時(shí)期內(nèi)的傷亡人次數(shù)、千人死亡率、千人負(fù)傷率、百萬(wàn)噸產(chǎn)品死亡率等指標(biāo)來(lái)表示。所有這些指標(biāo)，都可以通過(guò)預(yù)測(cè)，對(duì)其未來(lái)的變化作出估計(jì)。&回歸分析法具有預(yù)測(cè)結(jié)果比較接近實(shí)際、易于表示數(shù)據(jù)的離散性并給出預(yù)測(cè)區(qū)間等優(yōu)點(diǎn)，在工礦企業(yè)傷亡事故趨勢(shì)預(yù)測(cè)中已得到了廣泛的應(yīng)用。 &回歸分析法是研究相關(guān)關(guān)系的一種數(shù)理統(tǒng)計(jì)方法(fngf)，分為一元回歸和多元回歸。 第26頁(yè)/共69頁(yè)第二十六頁(yè)，共70頁(yè)。27回歸分析法預(yù)測(cè)(yc)的步驟（1）根據(jù)實(shí)驗(yàn)或觀察數(shù)據(jù)，繪制散點(diǎn)圖，大體確定變量之間的相關(guān)關(guān)系；（2）根據(jù)散點(diǎn)圖初步確定相關(guān)關(guān)系方程表達(dá)式的類型(lixng)，建立經(jīng)

18、驗(yàn)回歸方程，從而對(duì)變量之間的關(guān)系程度進(jìn)行精確的計(jì)算與分析。散點(diǎn)圖是利用有對(duì)應(yīng)關(guān)系的兩個(gè)變量分別作為坐標(biāo)，且將這兩個(gè)變量的統(tǒng)計(jì)值標(biāo)在該坐標(biāo)系中所成的圖形上。在繪制散點(diǎn)圖之前，應(yīng)先根據(jù)實(shí)驗(yàn)或觀察取得一組互相對(duì)應(yīng)的數(shù)據(jù)編制成數(shù)據(jù)表，然后根據(jù)散點(diǎn)圖進(jìn)行計(jì)算和分析。 第27頁(yè)/共69頁(yè)第二十七頁(yè)，共70頁(yè)。28相 關(guān) 關(guān) 系 的 分 類（1）從相關(guān)的性質(zhì)分為正相關(guān)和負(fù)相關(guān)；（2）從影響因素的多少分為單相關(guān)和復(fù)相關(guān)。單相關(guān)是兩個(gè)現(xiàn)象之間的關(guān)系；（3）從相關(guān)的表現(xiàn)形式分為直線相關(guān)和曲線相關(guān)。（4）從相關(guān)緊密程度分為完全相關(guān)、不完全相關(guān)和不相關(guān)。完全相關(guān)為函數(shù)關(guān)系，兩現(xiàn)象各自(gz)獨(dú)立，毫無(wú)關(guān)系則為不相關(guān)，

19、在完全相關(guān)和不相關(guān)之間則為不完全相關(guān)(安全系統(tǒng)多屬于此)。第28頁(yè)/共69頁(yè)第二十八頁(yè)，共70頁(yè)。296.3.2 一元線性回歸(hugu)法預(yù)測(cè)第29頁(yè)/共69頁(yè)第二十九頁(yè)，共70頁(yè)。30回歸(hugu)直線方程的確定&最小二乘法要求(yoqi)y的修均值和實(shí)際值的離差平方和為最小，即&設(shè)有n對(duì)x與y的數(shù)值(shz)，若y的修均值值以a+ bx代入，則此離差平方和成為a與b的函數(shù)，用 W(a,b) 表示，即最小值niiyy12)(2)(),(bxaybaW第30頁(yè)/共69頁(yè)第三十頁(yè)，共70頁(yè)。31回歸直線方程(fngchng)的確定&為了(wi le)使W(a,b) 成

20、為最小，可分別求W(a,b)對(duì)a及b的偏導(dǎo)且令其等于0，整理后得： &由上述方程可求得參數(shù)(cnsh)a、b分別為：2xbxaxyxbnay 222)(xnxyxxyxa22)(xnxxynyxb第31頁(yè)/共69頁(yè)第三十一頁(yè)，共70頁(yè)。32回歸(hugu)分析法預(yù)測(cè)的舉例&一元線性回歸用于企業(yè)事故趨勢(shì)分析時(shí)，方程式中各變量代表的具體意義為：x時(shí)間順序號(hào)；y事故數(shù)據(jù)；n事故數(shù)據(jù)總數(shù)。&例：某企業(yè)近10年來(lái)的事故傷亡人數(shù)如下頁(yè)表所示，現(xiàn)用一元線性回歸法預(yù)測(cè)(yc)事故的發(fā)展趨勢(shì)。 &解 首先，根據(jù)傷亡人數(shù)的統(tǒng)計(jì)數(shù)值繪制散點(diǎn)圖，得出傷亡人數(shù)與時(shí)間的關(guān)系為直線關(guān)系。然后

21、，求出參數(shù)a、b： 2 .80385105558810291355)(2222 xnxyxxyxa第32頁(yè)/共69頁(yè)第三十二頁(yè)，共70頁(yè)。33回歸(hugu)直線方程的確定舉例時(shí)間順序號(hào)時(shí)間順序號(hào)(x)傷亡人數(shù)傷亡人數(shù)(y)x2xyy21234567891066757394683435325358149162536496481100661502193763402042452564775804356562553298836462411561225102428093364第33頁(yè)/共69頁(yè)第三十三頁(yè)，共70頁(yè)。34回歸直線方程的確定(qudng)舉例89. 3385105529131058855)

22、(222xnxxynyxb38348,2913,385,588,5522yxyxyx第34頁(yè)/共69頁(yè)第三十四頁(yè)，共70頁(yè)。35（2）回歸(hugu)分析法的相關(guān)關(guān)系&因此回歸(hugu)直線方程為：y = 80.2 3.89 x 在回歸分析中，還應(yīng)研究計(jì)算得到(d do)的回歸直線是否符合實(shí)際數(shù)據(jù)變化的趨勢(shì)。為此引入相關(guān)系數(shù)r 的概念，其計(jì)算公式為： yyxxxyLLLr ,1yxnxyLxy22)(1xnxLxx22)(1ynyLyy第35頁(yè)/共69頁(yè)第三十五頁(yè)，共70頁(yè)。36相關(guān)系數(shù)取不同的數(shù)值(shz)時(shí)，分別表示實(shí)際數(shù)據(jù)和回歸直線之間的不同符合情況。 1）r = 0時(shí)，表示回

23、歸直線不符合實(shí)際數(shù)據(jù)的變化情況； （2）回歸分析法的相關(guān)(xinggun)關(guān)系第36頁(yè)/共69頁(yè)第三十六頁(yè)，共70頁(yè)。37 2）0 |r| 1時(shí)，表示回歸直線在一定程度上符合實(shí)際(shj)數(shù)據(jù)的變化趨勢(shì)。|r|越大，說(shuō)明回歸直線與實(shí)際(shj)數(shù)據(jù)變化趨勢(shì)的符合程度越大；|r|越小，則符合程度越小。 （2）回歸(hugu)分析法的相關(guān)關(guān)系第37頁(yè)/共69頁(yè)第三十七頁(yè)，共70頁(yè)。38 3）|r| =1時(shí)，表示回歸(hugu)直線完全符合實(shí)際數(shù)據(jù)的變化情況。 （2）回歸分析法的相關(guān)(xinggun)關(guān)系第38頁(yè)/共69頁(yè)第三十八頁(yè)，共70頁(yè)。39（3）回歸分析法的預(yù)測(cè)(yc)與控制2)(2nyyS

24、y在 回 歸 分 析 中 ， 還 應(yīng) 根 據(jù) 回 歸 方 程 來(lái) 預(yù) 測(cè)(yc) y 的取值范圍。當(dāng) n 較大時(shí)，y 的剩余均方差為： 當(dāng)x = x0 時(shí)，相應(yīng)的 y0 服從正態(tài)分布，則 y0 落在y02Sy0區(qū)間上的概率(gil)為0.9545，因此可得到y(tǒng) 的預(yù)測(cè)區(qū)間為y0-2Sy0， y0+2Sy0，也可求出 y 在某區(qū)間內(nèi)取值，相應(yīng) x 在什么范圍。第39頁(yè)/共69頁(yè)第三十九頁(yè)，共70頁(yè)。40（3）回歸(hugu)分析法的預(yù)測(cè)與控制xy89. 32 .80仍以前面的傷亡數(shù)據(jù)為例,說(shuō)明預(yù)測(cè)區(qū)間的求法。 該例的回歸直線(zhxin)過(guò)程為： 可以將前面表中xi 的分別代入方程求得y的均值，

25、再以實(shí)際 yi 分別減去y的均值，即可計(jì)算(j sun)得到8 .1721063.25242)(2nyySy 據(jù)此可預(yù)測(cè)y的取值范圍。例如，設(shè) x = 11，則 y = 80.2 3.8911 = 37.41第40頁(yè)/共69頁(yè)第四十頁(yè)，共70頁(yè)。41 相應(yīng)預(yù)測(cè)區(qū)間為：37.41217.8 = 37.41 35.6，即1.81,73.01，y值落在1.8173.01內(nèi)的概率為0.9545?；蛘?huzh)37.4117.8 ，即19.61,55.21。 我們可以作出預(yù)測(cè)帶如圖所示。（3）回歸(hugu)分析法的預(yù)測(cè)與控制第41頁(yè)/共69頁(yè)第四十一頁(yè)，共70頁(yè)。426.3.3 一元非線性回歸(hu

26、gu)預(yù)測(cè)法 非線性回歸分析方法是通過(guò)一定的變換，將非線性問(wèn)題轉(zhuǎn)化為線性問(wèn)題，然后利用線性回歸的方法進(jìn)行回歸分析。 非線性回歸曲線有很多種，選用哪一種曲線作為(zuwi)回歸分析則要根據(jù)實(shí)際數(shù)據(jù)在坐標(biāo)系中的變化分布形狀，也可根據(jù)專業(yè)知識(shí)確定分析曲線。 常用的非線性回歸曲線有以下幾種。 （1）雙曲線xbay1則有，令xxyy1,1xbay第42頁(yè)/共69頁(yè)第四十二頁(yè)，共70頁(yè)。43（2）冪函數(shù)baxy 則有，令aaxxyylglg,lgxbay（3）指數(shù)函數(shù)(zh sh hn sh)bxaey ) 1則有，令aayyln,lnbxayxbaey/)2則有，令aaxxyyln,1,lnxbay6.

27、3.3 一元非線性回歸(hugu)預(yù)測(cè)法第43頁(yè)/共69頁(yè)第四十三頁(yè)，共70頁(yè)。44（4）對(duì)數(shù)函數(shù)(du sh hn sh)xbaylg則有，令xxlgxbay下面以指數(shù)函數(shù) y = aebx 為例，說(shuō)明非線性曲線的回歸方法例 某鋼廠上一年的工傷人數(shù)(rn sh)的統(tǒng)計(jì)數(shù)據(jù)如下頁(yè)表所示，試用指數(shù)函數(shù) y = aebx 進(jìn)行回歸分析。則有，令aayyln,lnbxay6.3.3 一元(y yun)非線性回歸預(yù)測(cè)法第44頁(yè)/共69頁(yè)第四十四頁(yè)，共70頁(yè)。月份月份時(shí)間順序號(hào)時(shí)間順序號(hào)(x) 工傷人數(shù)工傷人數(shù)(y)yx2xyy21234567891011121234567891011121512764

28、56744212.7082.4851.9461.7921.3861.6091.7921.9461.3861.3860.69301491625364964811001211442.7084.9705.8387.1686.9309.65412.5415.5712.4713.867.62307.3336.1753.7873.2111.9312.5893.2113.787.01.9210.4800第45頁(yè)/共69頁(yè)第四十五頁(yè)，共70頁(yè)。46一元(y yun)非線性回歸曲線第46頁(yè)/共69頁(yè)第四十六頁(yè)，共70頁(yè)。47一元非線性回歸預(yù)測(cè)(yc)法舉例 故指數(shù)回歸(hugu)曲線方程為:37. 265012

29、78129.19650337.9978)(2222 xnxyxyxxa175. 06501278337.9912129.1978)(222xnxyxnyxb34.15ln73. 2eeaaaa得：由xey175. 034.15n計(jì)算計(jì)算(j sun)相關(guān)系數(shù)：相關(guān)系數(shù)：25129.1978121337.991yxnyxLyx第47頁(yè)/共69頁(yè)第四十七頁(yè)，共70頁(yè)。48一元非線性回歸(hugu)預(yù)測(cè)法舉(1222xnxLxx84. 5129.19121336.36)(1222ynyLyy 0 |r| 1，說(shuō)明用指數(shù)曲線進(jìn)行回歸分析，在一定程度上反映了該廠實(shí)際工傷人數(shù)的變

30、化趨勢(shì)。 還可進(jìn)一步計(jì)算(j sun)出 y 的預(yù)測(cè)區(qū)間。87. 084. 514325yyxxyxLLLr第48頁(yè)/共69頁(yè)第四十八頁(yè)，共70頁(yè)。496 6 系統(tǒng)安全預(yù)測(cè)系統(tǒng)安全預(yù)測(cè)(yc)(yc)6.1 安全安全(nqun)預(yù)測(cè)預(yù)測(cè)概述概述6.2 德?tīng)柗祁A(yù)測(cè)法德?tīng)柗祁A(yù)測(cè)法6.3 回歸分析法預(yù)測(cè)回歸分析法預(yù)測(cè)6.4 馬爾柯夫鏈預(yù)測(cè)法馬爾柯夫鏈預(yù)測(cè)法6.5 灰色系統(tǒng)預(yù)測(cè)法灰色系統(tǒng)預(yù)測(cè)法第49頁(yè)/共69頁(yè)第四十九頁(yè)，共70頁(yè)。506.4 馬爾柯夫鏈預(yù)測(cè)(yc)法 o如果事物的發(fā)展過(guò)程及狀態(tài)只與事物當(dāng)時(shí)的狀態(tài)有關(guān)，而與以前狀態(tài)無(wú)關(guān)(wgun)時(shí)，則此事物的發(fā)展變化稱為馬爾柯夫鏈。o如果系統(tǒng)的安全狀

31、況具有馬爾柯夫性質(zhì)，且一種狀態(tài)轉(zhuǎn)變?yōu)榱硪环N狀態(tài)的規(guī)律又是可知的，那么可以利用馬爾柯夫鏈的概念進(jìn)行計(jì)算和分析，來(lái)預(yù)測(cè)未來(lái)特定時(shí)刻的系統(tǒng)安全狀態(tài)。o馬爾柯夫鏈?zhǔn)潜碚饕粋€(gè)系統(tǒng)在變化過(guò)程中的特性狀態(tài)，可用一組隨時(shí)間進(jìn)程而變化的變量來(lái)描述。第50頁(yè)/共69頁(yè)第五十頁(yè)，共70頁(yè)。51馬爾柯夫鏈預(yù)測(cè)(yc)法 o如果系統(tǒng)在任何時(shí)刻上的狀態(tài)是隨機(jī)性的，則變化過(guò)程是一個(gè)隨機(jī)過(guò)程，當(dāng)時(shí)刻 t 變到時(shí)刻 t+1，狀態(tài)變量從某個(gè)取值變到另一個(gè)取值，系統(tǒng)就實(shí)現(xiàn)了狀態(tài)轉(zhuǎn)移(zhuny)。o系統(tǒng)從某種狀態(tài)轉(zhuǎn)移(zhuny)到各種狀態(tài)的可能性大小，可用轉(zhuǎn)移(zhuny)概率來(lái)描述。o假定系統(tǒng)的初始狀態(tài)可用狀態(tài)向量表示為： ,

32、)0()0(3)0(2)0(1)0(nsssss狀態(tài)轉(zhuǎn)移概率(gil)矩陣為: 第51頁(yè)/共69頁(yè)第五十一頁(yè)，共70頁(yè)。526.4.2 狀態(tài)(zhungti)轉(zhuǎn)移概率矩陣及其性質(zhì) o狀態(tài)(zhungti)轉(zhuǎn)移矩陣是一個(gè) n 階方陣，滿足概率矩陣的一般性質(zhì)，即有 0 pij 1 且 。nnnnnnpppppppppp212222111211o也就是說(shuō)，狀態(tài)轉(zhuǎn)移矩陣的所有行向量都是概率(gil)向量。11njijp第52頁(yè)/共69頁(yè)第五十二頁(yè)，共70頁(yè)。53o一次轉(zhuǎn)移(zhuny)向量 s(1) 為:pss)0()1(o二次轉(zhuǎn)移(zhuny)向量 s(2) 為:1)0()1(kkpsso類似(li

33、 s)地2)0()1()2(pspss6.4.3 安全預(yù)測(cè) 第53頁(yè)/共69頁(yè)第五十三頁(yè)，共70頁(yè)。54馬爾柯夫鏈預(yù)測(cè)(yc)法舉例 例：某單位對(duì)例：某單位對(duì)1250名接觸矽塵人員進(jìn)行健康檢名接觸矽塵人員進(jìn)行健康檢查時(shí)，發(fā)現(xiàn)查時(shí)，發(fā)現(xiàn)(fxin)職工的健康狀況分布如表職工的健康狀況分布如表所示。所示。 健健 康康 狀狀 況況 健健 康康疑疑 似似 矽矽 肺肺 矽矽 肺肺代代 表表 符符 號(hào)號(hào) s1(0)s2(0)s3(0)人人 數(shù)數(shù)100020050o根據(jù)統(tǒng)計(jì)資料，一年后接塵人員的健康變化規(guī)律為：健康人員繼續(xù)保持健康者剩70%，有20%變?yōu)橐伤?y s)矽肺，10%的人被定為矽肺，即10. 0

34、,20. 0,70. 0131211ppp第54頁(yè)/共69頁(yè)第五十四頁(yè)，共70頁(yè)。55馬爾柯夫鏈預(yù)測(cè)(yc)法舉例 o原有疑似矽肺者一般不可能恢復(fù)(huf)為健康者，仍保持原狀者為80%，有20%被正式定為矽肺，即 2 . 0, 8 . 0, 0232221ppp 因此(ync)，狀態(tài)轉(zhuǎn)移矩陣為 o矽肺患者一般不可能恢復(fù)為健康或返回疑似矽肺，即 1, 0, 0333231ppp0 . 1002 . 08 . 001 . 02 . 07 . 0p第55頁(yè)/共69頁(yè)第五十五頁(yè)，共70頁(yè)。56馬爾柯夫鏈預(yù)測(cè)(yc)法舉例 預(yù)測(cè)(yc)一年后接塵人員的健康狀況為： 333231232221131211

35、)0(3)0(2)0(1)0()1(pppppppppssspss0 . 1002 . 08 . 001 . 02 . 07 . 0502001000190360700即一年后，仍然健康者為700人，疑似矽肺者360人，被定為矽肺者190人。預(yù)測(cè)表明，該單位矽肺發(fā)展速度很快，必須加強(qiáng)(jiqing)防塵工作和醫(yī)療衛(wèi)生工作。 第56頁(yè)/共69頁(yè)第五十六頁(yè)，共70頁(yè)。576 6 系統(tǒng)安全預(yù)測(cè)系統(tǒng)安全預(yù)測(cè)(yc)(yc)6.1 安全預(yù)測(cè)概述安全預(yù)測(cè)概述(i sh)6.2 德?tīng)柗祁A(yù)測(cè)法德?tīng)柗祁A(yù)測(cè)法6.3 回歸分析法預(yù)測(cè)回歸分析法預(yù)測(cè)6.4 馬爾柯夫鏈預(yù)測(cè)法馬爾柯夫鏈預(yù)測(cè)法6.5 灰色系統(tǒng)預(yù)測(cè)法灰色系統(tǒng)

36、預(yù)測(cè)法第57頁(yè)/共69頁(yè)第五十七頁(yè)，共70頁(yè)。586.5 灰色(hus)系統(tǒng)預(yù)測(cè)法 o灰色系統(tǒng)理論預(yù)測(cè)通過(guò)一系列數(shù)據(jù)生成方法，如直接累加法、移動(dòng)平均法、加權(quán)累加法、自適應(yīng)性累加法等，將根本沒(méi)有(mi yu)規(guī)律的、雜亂無(wú)章的或規(guī)律性不強(qiáng)的一組原始數(shù)據(jù)變得具有明顯的規(guī)律性，解決了數(shù)學(xué)界一直認(rèn)為不能解決的微分方程建模問(wèn)題。o灰色系統(tǒng)預(yù)測(cè)是從灰色系統(tǒng)的建模、關(guān)聯(lián)度及殘差辨識(shí)的思想出發(fā)，所獲得的關(guān)于預(yù)測(cè)的概念、觀點(diǎn)與方法。將灰色系統(tǒng)理論用于廠礦企業(yè)事故預(yù)測(cè)，一般選用 GM(1,1) 模型，是一階的一個(gè)變量的微分方程模型。 第58頁(yè)/共69頁(yè)第五十八頁(yè)，共70頁(yè)。596.5.1 灰色系統(tǒng)預(yù)測(cè)(yc)建模

37、方法 設(shè)原始離散(lsn)數(shù)據(jù)序列為：o其中 N 為序列長(zhǎng)度，對(duì)其進(jìn)行(jnxng)一次累加生成處理后得,)0()0(3)0(2)0(1)0(NxxxxxNkxxkjjk, 2 , 11)0()1(o則以生成序列,)1()1(3)1(2)1(1)1(Nxxxxx 為基礎(chǔ)建立的灰色生成模型為uxdtdx)1()1(第59頁(yè)/共69頁(yè)第五十九頁(yè)，共70頁(yè)。60o稱為(chn wi)一階灰色微分方程，記為 GM(1,1) 。o 式中，、u 為待辨識(shí)參數(shù)。 o設(shè)參數(shù)(cnsh)向量TNNTxxxyu,)0()0(3)0(2 則由下式求得方程(fngchng)的最小二乘解： 12/ )(12/ )()1

38、(1)1()1(1)1(2NNxxxxB6.5.1 灰色系統(tǒng)預(yù)測(cè)建模方法 第60頁(yè)/共69頁(yè)第六十頁(yè)，共70頁(yè)。61的值作累減還原(hun yun)，即得到原始數(shù)據(jù)的估計(jì)值： o時(shí)間(shjin)響應(yīng)方程為：NTTyBBB1)(/)/()1(1)1(ueuxxtto離散響應(yīng)(xingyng)方程為：/)/()1(1)1(1ueuxxkk 式中 )0(1)1(1xx將)1(1kx)1()1(1)0(1kkkxxx6.5.1 灰色系統(tǒng)預(yù)測(cè)建模方法 第61頁(yè)/共69頁(yè)第六十一頁(yè)，共70頁(yè)。62o用上述方法(fngf)建立累加殘差生成模型： oGM(1,1) 模型的擬合殘差中往往還有一部分動(dòng)態(tài)有效信息，可以通過(guò)建立(jinl)殘差 GM(1,1)模型對(duì)原模型進(jìn)行修正。 。一般NNN,)1()1(2)1(1)1( 式中，1和u1為殘差模型(mxng)參數(shù)。累減后得(1)的還原估計(jì)值： o記殘差)1()1()1(kkkxx組成的序列為 1111)1(1)1(1/)/(1ueukk6.5.1 灰色系統(tǒng)預(yù)測(cè)建模方法 第62頁(yè)/共69頁(yè)第六十二頁(yè)，共70頁(yè)。63o若殘差模型是對(duì)第 m 個(gè)殘差開(kāi)始進(jìn)行擬合的，則修正(xizhng)后