《1 概述》
1 概述
通常在进行卫星故障检测和诊断时, 人们往往想到确定性推理理论。然而在卫星故障检测和诊断的许多方面, 如症状信息获取、决策信息系统生成和约简等方面都不可避免地要运用合情推理和容错能力。不确定性推理理论可以实现合情推理和容错能力, 而确定性推理理论却不能。因而仅仅应用确定性推理理论是难以对某些卫星故障进行自动检测和诊断的。目前, 人们正在研究不确定性推理理论在卫星故障检测和诊断中的应用, 相应的系统也正在开发。不确定性推理理论包括若干理论。笔者研究的卫星故障检测和诊断的三种新方法, 分别基于不确定性推理理论中的包含度理论、粗糙集理论和证据推理理论
《2 三种新方法简介》
2 三种新方法简介
《2.1包含度理论方法》
2.1包含度理论方法
设
·X为对象集合, Ai⊂X (i≤k) 为X的分划, 即满足Ai∩Aj=Φ (i≠j) , 且
·Ax 是X的一个分划, Ax={ Ai; i≤k }。
·D是X上分划的全体。
·X上的两个分划Ax={ Ai; i≤k } 和Bx={Bj; j≤l}。
若有Ai⊂Bj, 称Bx 依赖于Ax (Ax⊂Bx) 。 如果Ai∈Ax, Bj∈Bx 且Ai⊂Bj, 称Ai 影响到Bj。
D (Bj / Ai) 和D (Bx / Ax) 分别是X和D上的包含度。 包含度定义如下:
若
如果是 Ai, 则为 Bji 。
《2.2粗糙集理论方法》
2.2粗糙集理论方法
定义1 决策信息系统可表示为S=〈U, R, V, f〉。U为对象集, R=C∪D为属性集 (C为条件属性, D为决策属性) 。
IND (B) 为一等价关系, 且
定义2 U为对象集, P是定义在U上的一个等价关系簇, R∈P。如果IND (PR}) =IND (P) , R对P是绝对不必要的, 否则, 是绝对必要的。
定义3 U为对象集, P是定义在U上的一个等价关系簇, R∈P。如果每一个R∈P对P都是绝对必要的, 则P是独立的, 否则, 是互相依赖的。
定义4 U为对象集, P是定义在U上的一个等价关系簇。P簇所有绝对必要关系组成的集合为关系簇P的绝对核 (CORE (P) ) 。
定义5 U为对象集, P和Q是定义在U上的两个等价关系簇且Q⊂P。如果 IND (Q) =IND (P) , 且Q是独立的, 则Q是P的一个绝对约简。
定义6 对X和不分明关系B, 包含于X的最大可定义集B- (X) 和最小可定义集B- (X) 都能够根据B确定。B- (X) 和B- (X) 分别称为X的B下逼近和上逼近:
定义7 U为对象集, P和Q是定义在U上的两个等价关系簇, Q的P正域为
定义8 U为对象集, P和Q是定义在U上的两个等价关系簇。如果POSP (Q) =POS (Pr}) (Q) , 则r是P中相对于Q可省略的, 否则, 是不可省略的。
定义9 U为对象集, P和Q是定义在U上的两个等价关系簇。如果POSP (Q) ≠POS (Pr}) (Q) , P为相对于Q独立的。
定义10 U为对象集, P和Q是定义在U上的两个等价关系簇, 且S⊂P。如果 POSS (Q) =POSP (Q) , 则S为P的Q约简 (REDQ (P) ) 。
定义11 U为对象集, P和Q是定义在U上的两个等价关系簇, P的所有Q不可省略关系簇为P的Q核 (COREQ (P) ) 。
定义12 U为对象集, P和Q是定义在U上的两个等价关系簇, 如果POSP (Q) =U, 则U是P上相对于Q一致的。
定义13 决策信息系统条件属性集合C的相对约简C′是C相对于决策属性D的最大独立子集。
《2.3证据推理理论方法》
2.3证据推理理论方法
证据推理理论的关键概念是mass函数。mass函数的值是由专家依据他们的态度、经验、知识等确定的。对于同一决策问题, 各位专家给出的值往往不同。如果仅应用确定性推理理论, 依据专家们给出的mass函数值进行决策是困难的。若应用证据推理理论, 决策就比较容易, 因为证据推理理论可以合理地合并这些mass函数。证据推理理论的这一特点对于卫星故障检测和诊断是非常有用的。在卫星故障检测和诊断实际中, 经常会有这样的情况:对于同一个现象, 专家们根据其经验和知识给出的结果是不同的, 其原因是定位卫星故障和确定导致卫星故障的因素都是困难的。此时, 就可以应用证据推理理论合成各位专家的不同观点, 得到合理的结果
mass函数合成方法如下 (H为假设空间) :
其中
设m是X上的mass函数, 则:
近似精度μ=B (A) /L (A) 。通过求μmax, 确定导致卫星故障的原因。
应用证据推理理论时, 要注意以下五点:
·导致卫星故障的因素是不相关的。
· 可以对各位专家的mass函数给出不同的权重, 导致不同因素故障的权重也可以不同。
·一种故障通常可能由多种因素同时引起, 应启发专家在确定mass函数值时考虑多种因素。
·如果不同因素的μ比较接近, 则故障可能是由这几种因素共同作用的结果。
·要考虑非显著因素的综合效果。若非显著因素的综合效果大于某个显著因素的效果, 则应该用另一种方法进行故障诊断, 以验证证据推理理论的诊断结果。
《3 三种新方法的实际应用》
3 三种新方法的实际应用
《3.1包含度理论的应用》
3.1包含度理论的应用
·卫星是否穿越地球北极;
·是否实施轨道控制;
·星上主份陀螺是否在线。
根据包含度理论, 卫星姿态控制计算机故障信息系统的原始形式如表1。
Table 1 Information system of the fault of satellite attitude control computer (original form)
《表1》
故障 对象 |
症状 |
决策 | ||
穿越地球北极 |
轨道控制 | 主份陀螺在线 | ||
F1 |
是 | 是 | 是 | D1 |
F2 |
是 | 否 | 否 | D2 |
F3 |
否 | 是 | 否 | D3 |
F4 |
否 | 是 | 否 | D3 |
F5 |
是 | 是 | 是 | D1 |
F6 |
是 | 是 | 是 | D1 |
F7 |
是 | 否 | 否 | D2 |
F8 |
否 | 是 | 否 | D3 |
定义 F表示故障, C表示故障症状属性, D表示故障决策属性。
症状属性中的“是”用1表示, “否”用0表示;决策属性中的“D1”用1表示, “D2”用2表示, “D3”用3表示。为表述方便, 用PNP表示“穿越地球北极”, 用OM表示“轨道控制”, 用MGO表示“主份陀螺在线”, 则卫星姿态控制计算机故障信息系统的形式化形式如表2。
Table 2 Information system of the fault of satellite attitude control computer (formulary form)
《表2》
F |
C |
D | ||
PNP |
OM | MGO | ||
F1 |
1 | 1 | 1 | 1 |
F2 |
1 | 0 | 0 | 2 |
F3 |
0 | 1 | 0 | 3 |
F4 |
0 | 1 | 0 | 3 |
F5 |
1 | 1 | 1 | 1 |
F6 |
1 | 1 | 1 | 1 |
F7 |
1 | 0 | 0 | 2 |
F8 |
0 | 1 | 0 | 3 |
分别按PNP, OM和MGO划分F:
则PNP对F的划分 (PNP*) 为
则OM对F的划分 (OM*) 为
则MGO对F的划分 (MGO*) 为
于是, 得到关于F依不同属性的所有划分:
按D划分F:
则D对F的划分 (D*) 为
由于PNP*>D*, OM*>D*, MGO*>D*, 所以, 依赖PNP*, OM*, MGO*其中之一不能形成决策规则。但是, PNP*×OM*≤D*, OM*×MGO*≤D*, 表明可以用PNP与OM或OM与MOG划分D。
例如:
则有
因而得到以下决策规则:
规则1:
如果 (PNP=1, OM=1) , 则 D=D1;
如果 (PNP=1, OM=0) , 则 D=D2;
如果 (PNP=0, OM=1) , 则 D=D3。
规则2:
如果 (OM=1, MGO=1) , 则 D=D1;
如果 (OM=0, MGO=0) , 则 D=D2;
如果 (OM=1, MGO=0) , 则 D=D3。
《3.2粗糙集理论的应用》
3.2粗糙集理论的应用
设引起卫星姿态故障的姿态控制部件包括推力器、陀螺、动量轮、磁力矩器。为表述方便, “推力器”用Thruster表示, “陀螺”用Gyro表示、“动量轮”用MW表示, “磁力矩器”用MT表示
Table 3 Symptom information
《表3》
部件 |
名称 | 症状信息 | ||
P1 |
Thruster | TH1 | TH2 | TH3 |
P2 |
Gyro | G1 | G2 | G3 |
P3 |
MW | M1 | M2 | |
P4 |
MT | TO1 | TO2 |
为表述方便, 用TH1, TH2, TH3分别表示Thruster遥测数据的三种取值;用G1, G2, G3分别表示Gyro遥测数据的三种取值;用M1, M2分别表示MW遥测数据的两种取值;用TO1, TO2分别表示MT遥测数据的两种取值。
根据卫星历史测量数据得到的初始决策信息系统如表4。
Table 4 Decision information system
《表4》
部件属性值 | 故障状态 | |||
TH1 |
G3 | M2 | TO2 | 无 |
TH1 |
G3 | M2 | TO1 | 无 |
TH2 |
G3 | M2 | TO2 | 有 |
TH3 |
G2 | M2 | TO2 | 有 |
TH3 |
G1 | M1 | TO2 | 有 |
TH3 |
G1 | M1 | TO1 | 无 |
TH2 |
G1 | M1 | TO1 | 有 |
TH1 |
G2 | M2 | TO2 | 无 |
TH1 |
G1 | M1 | TO2 | 有 |
TH3 |
G2 | M1 | TO2 | 有 |
TH1 |
G2 | M1 | TO1 | 有 |
TH2 |
G2 | M2 | TO1 | 有 |
定义C为条件属性集 (部件属性) :
C={推力器 (a1) , 陀螺 (a2) , 动量轮 (a3) , 磁力矩器 (a4) }。
设D为决策属性集:
则有如下“不分明关系 (IND) ”和“正域 (POS) ”:
IND (C) ={{1}, {2}, {3}, {4}, {5}, {6}, {7}, {8}, {9}, {10}, {11}, {12}};
IND (D) ={{1, 2, 6, 8}, {3, 4, 5, 7, 9, 10, 11, 12}};
POSC (D) =U;
IND (C\a1) ={{1, 3}, {2}, {4, 8}, {5, 9}, {6, 7}, {10}, {11}, {12}};
IND (C\a2) ={{1, 8}, {2}, {3}, {4}, {5, 10}, {6}, {7}, {9}, {11}, {12}};
IND (C\a3) ={{1}, {2}, {3}, {4, 10}, {5}, {6}, {7}, {8}, {9}, {11}, {12}};
IND (C\a4) ={{1, 2}, {3}, {4}, {5, 6}, {7}, {8}, {9}, {11}, {12}}。
由此可得
因此, 在判断该故障时, 陀螺或动量轮的因素是不必要的。但是:
所以, 条件属性集{推力器 (a1) 、陀螺 (a2) 、磁力矩器 (a4) }和{推力器 (a1) 、动量轮 (a3) 、磁力矩器 (a4) }能够决定决策属性集{Decision (d) }, C的所有D约简为
C的D核为
根据基于辨识矩阵和逻辑运算的决策表属性约简算法, 得到决策信息系统的两个约简A和B, 如表5, 表6。
Table 5 Reduction A of the decision information system (Thruster, Gyro, MT)
《表5》
部件属性值 | 故障状态 | ||
TH1 |
G3 | TO2 | 无 |
TH1 |
G3 | TO1 | 无 |
TH2 |
G3 | TO2 | 有 |
TH3 |
G2 | TO2 | 有 |
TH3 |
G1 | TO2 | 有 |
TH3 |
G1 | TO1 | 无 |
TH2 |
G1 | TO1 | 有 |
TH1 |
G2 | TO2 | 无 |
TH1 |
G1 | TO2 | 有 |
TH1 |
G2 | TO1 | 有 |
TH2 |
G2 | TO1 | 有 |
Table 6 Reduction B of the decision information system (Thruster, MW, MT)
《表6》
部件属性值 | 故障状态 | ||
TH1 |
M2 | TO2 | 无 |
TH1 |
M2 | TO1 | 无 |
TH2 |
M2 | TO2 | 有 |
TH3 |
M2 | TO2 | 有 |
TH3 |
M1 | TO2 | 有 |
TH3 |
M1 | TO1 | 无 |
TH2 |
M1 | TO1 | 有 |
TH1 |
M1 | TO2 | 有 |
TH1 |
M1 | TO1 | 有 |
TH2 |
M2 | TO1 | 有 |
《3.3证据推理理论的应用》
3.3证据推理理论的应用
设导致卫星姿态故障的部件有动量轮故障 (P1) 、推力器故障 (P2) 、姿态敏感器故障 (P3) , 专家给出以下判断 (YPj表示是该部件导致姿态故障, NPj表示不是该部件导致姿态故障)
mass 函数合并:
1) 合并m1 和m2, m = m1⊕m2
2) 合并m9 和 m6
3) 合并m3 和 m5
4) 合并m7 和 m10
5) 合并m4 和 m8
6) 最后
求μmax
诊断结果为“推力器故障导致姿态故障”。
《4 结论》
4 结论
笔者首先分析了仅使用确定性推理理论进行卫星故障检测和诊断的缺陷, 然后研究了应用不确定性推理理论中的包含度理论、粗糙集理论和证据推理理论进行卫星故障检测和诊断的方法, 最后给出了分别应用不确定性推理理论中的上述三种理论进行卫星故障检测和诊断的三个实例。通过这些分析和研究, 得出以下结论:
1) 同时应用确定性推理理论和不确定性推理理论, 可使卫星故障检测和诊断得到合理解决。
2) 由于人们经常遇到只应用确定性推理理论无法对某些卫星故障进行检测和诊断的困境, 因而应用以不确定性推理理论进行卫星故障检测和诊断的方法是十分有效和实用的。