楊繼坤，徐廷學，趙 亮，陳 卓

(1.海軍航空工程學院 研究生管理大隊，山東 煙臺264001;2.海軍航空工程學院 兵器科學與技術系，山東 煙臺264001;3.中國人民解放軍92076 部隊，北京102202)

0 引 言

對導彈武器系統(tǒng)可用度的評估，通常是根據裝備某個時期內統(tǒng)計的各項時間因素值來進行計算，由于統(tǒng)計時間周期長，影響因素眾多，這種方法不能很好地反映裝備真實可用度水平［1］。目前，研究導彈武器系統(tǒng)等可修系統(tǒng)的主要數學工具是隨機過程理論。文獻［2］評定了出現(xiàn)保障延誤時間的可修系統(tǒng)穩(wěn)態(tài)可用度。文獻［3-4］利用馬爾可夫更新過程分析了偵察系統(tǒng)的可靠性和裝甲裝備的使用可用度。文獻［5］基于更新過程進行了導彈武器系統(tǒng)的可用度分析。

在對武器系統(tǒng)進行可用度分析時，由于經濟上或者試驗環(huán)境難以模擬等原因，并不能確切知道影響可用度的時間分布和參數值，而只能根據已有的實驗和專家的經驗確定這些分布類型和參數區(qū)間，將其描述為隨機模糊變量是更客觀、更合理的。目前，隨機模糊理論在可靠性領域有一定的應用。文獻［6-7］基于模糊隨機理論建立了結構可靠性分析模型。文獻［8］基于模糊理論進行了導彈系統(tǒng)貯存可靠性仿真方法研究。

本文首先基于馬爾可夫更新過程和拉普拉斯變換工具，建立考慮技術陣地設備部分戰(zhàn)損情況下和具有保障延誤的串聯(lián)可修導彈武器系統(tǒng)可用度分析模型。針對可用度分析中存在模糊因素的特點，將導彈視作模糊系統(tǒng)，引入隨機模糊的概念，采用三角模糊數量化模糊信息，估算可用度水平。

1 馬爾可夫更新過程

設隨機變量Zn取值在E={0，1，…，K}中，隨機變量Tn取值在［0，∞)中，且有0=T0≤T1≤T2≤…≤Tn，n=0，1，…，如果對所有j ∈E，t ≥0，有P{Zn+1=j，Tn+1-Tn≤t| Z0，Z1，…，Zn，T0，T1，…，Tn}=P{Zn+1=j，Tn+1-Tn≤t| Zn}，則隨機過程(Z，T)={Zn，Tn，n=0，1，…}稱為狀態(tài)空間E 上的馬爾可夫更新過程。經簡化求得馬爾可夫更新方程為［3-4］

即:

其中，Ai(t)為系統(tǒng)可用度;gi(t)為定義在［0，∞)上的非負、在任何有限區(qū)間上的有界函數，在計算可用度時，通常這個函數是不同裝備服從任意分布的壽命、維修、保障延誤的時間。

在本文中，G(t)的Laplace 變換為

Laplace-Stieltjes 變換為

R(s)表示復變數s 的實部，Gk(t)表示G(t)的k 重卷積［9］，即，且G(0)(t)=1，t ≥0。

2 隨機模糊理論

設(Θ，P(Θ)，Cr)為可信性空間，其中Θ 為非空集合，P(Θ)為Θ 的冪集，Cr 為定義在P(Θ)上的可信性測度［10-12］。

定義1:假設ξ 為一從可信性空間(Θ，P(Θ)，Cr)到實數集R 上的函數，則ξ 為模糊變量。

定義2:令ξ 為可信性空間(Θ，P(Θ)，Cr)上的模糊變量，α ∈(0，1］，稱和分別為ξ 的α-悲觀值和α-樂觀值。

定義3:設ξ 為模糊變量，它的期望值E［ξ］是有限的，則。

定義4:若a=(al，am，an)，0 ＜al≤am≤al，稱a 為一個三角模糊數，如果其隸屬度函數可以表示為

定義5:設a 為三角模糊數，對?λ ∈(0，1)，a的截集aλ為有界閉區(qū)間［(am-al)λ+al，(am-an)λ+an］。

3 導彈武器系統(tǒng)可用度分析模型

在導彈武器系統(tǒng)的可用度研究中，大多數文獻考慮的系統(tǒng)都是假定部件故障后立即得到修理。但實際情況并非如此，而是常常有一段隨機的等待時間(稱這段時間為保障延誤時間)，例如申請等待備件等。而且在現(xiàn)實戰(zhàn)場環(huán)境下，敵方主要打擊的目標是導彈技術測試陣地，對技術陣地的打擊會損毀技術陣地的測試和維修設備(包括測試中的導彈)，使其部分戰(zhàn)損，需要對其進行更換(或修理)后再工作，這都將影響到導彈技術測試、對接任務的進行，使技術陣地不能及時提供裝配好的導彈供作戰(zhàn)使用。因此，在系統(tǒng)可用度分析中，考慮技術陣地測試設備部分戰(zhàn)損和保障延誤的可修系統(tǒng)十分必要，而且更有重要的理論意義和實用價值［13-14］。

3.1 模型構建步驟

1)明確系統(tǒng)組成。設系統(tǒng)為由n 個分系統(tǒng)串聯(lián)而成的可修系統(tǒng);

2)確定隨機模糊變量。時間分布未知的分系統(tǒng)和修理設備，要通過專家的知識與經驗來描述其分布類型和參數情況，這些量即為模糊變量;

3)建立基于馬爾可夫更新過程的可用度模型，得出可用度的解析模型(本章下面兩小節(jié)主要對這方面進行研究)。

4)運用三角模糊數估算可用度水平。將隨機模糊變量用三角模糊數來表示，根據隨機模糊相關定理，估算其α-悲觀值和α-樂觀值，再根據定義3 求解其穩(wěn)態(tài)可用度。

3.2 基本條件和假設

1)系統(tǒng)的第i 個部件壽命Xi服從指數分布Fi(t)=1-e-λit，t ≥0，λi＞0;

2)系統(tǒng)發(fā)生故障，如果不能及時得到維修，視為不可用狀態(tài)，產生的保障延誤時間Wi服從一般分布Hi(t)，其平均保障延誤時間為

3)進行故障維修時，修理時間Yi服從一般的分布Gi(t)，其平均修理時間為

修復后系統(tǒng)立即轉為工作狀態(tài);

4)當1 個部件故障等待修理或者正在修理時，其他部件停止工作，不再發(fā)生故障，此時系統(tǒng)處于故障狀態(tài)。部件修理完恢復如新，由指數分布的性質，知道在同一時刻不會有多個部件故障;

5)技術陣地修理設備在執(zhí)行任務過程中，會因部分戰(zhàn)損而出現(xiàn)失效，其服務期U 服從指數分布U(t)=1-e-αt，t ≥0，α ＞0，修理設備失效后立即進行更換，其更換時間V 統(tǒng)一服從一般分布V(t)，且平均更換時間為，修理設備更換如新;

6)Xi，Yi，Wi，U，V 相互獨立。

3.3 構建系統(tǒng)可用度模型

定義系統(tǒng)的廣義修理時間γ 為:在導彈技術陣地，測試和修理設備從開始測試和修理故障系統(tǒng)的時刻起，直到該部件技術準備完成的時間長度，其中包括測試與修理設備可能的戰(zhàn)損而對測試和修理設備進行更換的時間。令，則

其LS 變換和平均廣義修理時間為

A(t)=P{X(t)=1 | X(0)=1}。

若令τk，δk，γk分別表示系統(tǒng)的第k 個工作壽命、第k 個保障延誤時間和第k 個廣義修理時間，不難證明，{τk+δk+γk，k=1，2，…}是個更新過程。其分布函數為

對式(6)2 端作LS(Laplace-Stieltjes)變換得

式(8)右端第1 項中，當t ＜Xj時，時刻t 系 統(tǒng)必處于工作狀態(tài)，自然有X(t)=1;第2 項中，當Xj≤t ≤Xj+Wj+γj，即系統(tǒng)處于修理狀態(tài)，不可能與X(t)=1 同時發(fā)生，故第2 項為0;第3 項如下式推導所示。

從而

4 算例分析

假設某型導彈武器系統(tǒng)由n 個分系統(tǒng)串聯(lián)而成，通過可靠性分析，建立系統(tǒng)可靠性框圖如圖1所示。

圖1 系統(tǒng)可靠性框圖Fig.1 Reliablity block diagram of system

根據已有的實驗和專家經驗確定，其分系統(tǒng)和測試設備壽命服從指數分布，保障延誤時間和維修時間服從一般分布，應用三角模糊數來評定其隨機模糊參數如表1所示。

表1 不確定參數的三角模糊表示Tab.1 The uncertain parameters of triangular fuzzy representation

根據三角模糊數的相關定理，可以得出如下不同參數的α-悲觀值和α-樂觀值(此處用x 表示α):

其中x ∈(0，1］。

根據定義2 可知，可用度的悲觀值和樂觀值可分別表示為

由定義3 可得出

由算例分析結果可知，系統(tǒng)的可用度為0.5811，如果不考慮測試和維修設備故障問題，則β=0，系統(tǒng)的可用度為0.618 0，與經過前期MC 仿真得到的0.592 3 比較接近。結果表明，所建模型能夠很好地反映裝備實際組成和運行情況，該模型和算法適用于導彈武器系統(tǒng)可用度的分析。

5 結 語

本文在隨機模糊參數的背景下，研究了導彈武器系統(tǒng)可用度評估模型，通過分析可得到如下結論:

1)建立的馬爾可夫更新過程模型考慮到技術陣地設備部分戰(zhàn)損和保障延誤時間等因素，符合實際使用維修和作戰(zhàn)要求。

2)所建可用度模型中存在模糊因素，引入隨機模糊理論，采用三角模糊數量化模糊信息(分布類型和參數)能合理、有效地解決這一問題。

3)為復雜的導彈武器系統(tǒng)可用度研究提供了一條技術途徑，對改進武器系統(tǒng)，制導武器系統(tǒng)的維護和管理有促進作用。

［1］黃奇.兩不同部件串聯(lián)可修系統(tǒng)的可用度的一個新的計算方法［J］.數學理論與應用，2004，24(2):123-126.

［2］張幗奮，武洪萍.出現(xiàn)保障延誤時間的可修系統(tǒng)穩(wěn)態(tài)可用度的評定［J］.浙江大學學報，2006，33(3):268-271.

［3］沈吉鋒，張永志，等.基于馬爾科夫更新過程的偵查系統(tǒng)可靠性分析［J］.兵工自動化，2010，29(3):7-10.

［4］劉福勝，吳緯，單志偉.基于馬爾科夫更新過程的裝甲裝備使用可用度模型［J］.裝甲兵工程學院學報，2010，24(5):15-17.

［5］楊繼坤，徐廷學.基于更新過程的導彈武器系統(tǒng)可用度分析［J］.兵工自動化，2012，31(2):14-16.

［6］田艷芳，師義民，李鳳.基于隨機模糊理論的結構可靠性分析［J］.模糊系統(tǒng)與數學，2008，22(1):163-166.

［7］王磊，劉文.基于隨機模糊參數的結構模糊可靠性分析模型［J］.北京航空航天大學學報，2005，31(4):412-415.

［8］陳海建，滕克難，李波.基于模糊理論的導彈系統(tǒng)貯存可靠性仿真方法研究［J］.彈箭與制導學報，2011，31(2):33-36.

［9］黃卓，郭波.更新理論推論過程及應用［J］.數學的實踐與認識，2006，36(9):200-204.

［10］楊倫標，高英儀.模糊數學［M］.廣州:華南理工大學出版社，1995.

［11］LIU Y B，QIAO Z，WANG G Y.Fuzzy random reliability of structures based on fuzzy random variables［J］.Fuzzy Sets and Systems，1997，86:345-355.

［12］郭書祥，呂震宇.概率模型含有模糊參數時的模糊失效概率計算方法［J］.機械強度，2003，25(5):527-529.

［13］MACHERE Y，KOEHN P，SPARROW D.Improving reliability and operational availability of military systems［A］.IEEE Aerospace Conference［C］，2005.3489-3957.

［14］劉鯖潔，陳桂明，劉小方.BP 神經網絡在導彈武器系統(tǒng)使用可用度預測中的應用［J］.戰(zhàn)術導彈技術，2011(1):26-28.