一种液体火箭发动机启动阶段的Petri网量子贝叶斯故障诊断方法技术

本发明专利技术公开了一种液体火箭发动机启动阶段的Petri网量子贝叶斯故障诊断方法，该方法通过引入量子干涉，使得量子贝叶斯网络可以极大地影响概率推断，特别是当网络的不确定性水平非常高时(没有观察到证据)。通过建立关于不可观故障变迁的量子贝叶斯网络模型，利用可观节点触发关系和不可观节点的触发概率，手动选择量子参数计算故障节点触发概率，判断系统故障状态。本方法能够针对部分可观的离散动态系统进行故障诊断，其优点是能够对不可观的部分也进行准确的故障诊断。

【技术实现步骤摘要】
一种液体火箭发动机启动阶段的Petri网量子贝叶斯故障诊断方法
本专利技术涉及离散动态系统系统故障诊断领域，尤其涉及一种液体火箭发动机启动阶段的Petri网量子贝叶斯故障诊断方法。
技术介绍
推力器故障约占航天器、卫星姿态和轨道控制系统故障的四分之一，推力器故障会极大地改变整个航天器的行为，甚至可能导致整个任务的失败。因此，异常组件行为必须及早检测。航天器的故障检测和诊断往往依赖硬件和传感器冗余，但由于重量和成本等限制，给航天器加传感器并不总是成功的。为克服观测信息不全的问题，同时满足实时性要求，诊断系统一般采用基于定性模型的诊断方法，Petri网最早应用于故障识别与诊断，应用在电力系统、通信系统等。目前现有的故障诊断技术是Petri网，其最早应用于故障识别与诊断，应用在电力系统、通信系统等。通过运用一种基于增量算法，从当前可能状态及其概率的集合出发，利用概率模型提出了一种评估未来故障概率的方法。而部分可观随机Petri网，计算与给定的时间观测序列相一致的时间和非时间标记轨迹的概率，根据故障概率进行诊断。本专利技术创造提出一种液体火箭发动机启动阶段的Petri网量子贝叶斯故障诊断方法，由于量子干涉的引入，量子贝叶斯网络可以极大地影响概率推断，特别是当网络的不确定性水平非常高时(没有观察到证据)。通过建立关于不可观故障变迁的量子贝叶斯网络模型，利用可观节点触发关系和不可观节点的触发概率，手动选择量子参数计算故障节点触发概率，判断系统故障状态。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种液体火箭发动机启动阶段的Petri网量子贝叶斯故障诊断方法，本专利技术能够大大提升液体火箭发动机启动阶段故障诊断过程中的效率以及对故障诊断的准确性。为实现上述的目的，本专利技术提供一种液体火箭发动机启动阶段的Petri网量子贝叶斯故障诊断方法，其特征在于，包括如下步骤：步骤S1、根据液体火箭发动机启动阶段的工作状态建立液体火箭发动机启动阶段的Petri网模型；步骤S2、构建Petri网模型的SCG图，遍历SCG图中满足可观变迁触发序列的路径，并且判断是否包含故障变迁，估计液体火箭发动机的故障状态；步骤S3、建立QBPN模型，通过QBPN模型估计液体火箭发动机的故障概率。进一步的，所述步骤3具体包括：步骤S301、根据前向路径中的可观变迁触发状态计算故障变迁触发概率；步骤S302、利用QB函数根据后向路径可观变迁触发状态修正故障变迁触发概率；步骤S303、取所有故障变迁中最大触发概率作为液体火箭发动机的故障概率。本专利技术的有益效果是：本专利技术能够针对部分可观的离散动态系统进行故障诊断，其优点是能够对不可观的部分也进行准确的故障诊断。因为引入了量子干涉，故通过建立关于不可观故障变迁的量子贝叶斯网络模型，利用可观节点触发关系和不可观节点的触发概率，手动选择量子参数计算故障节点触发概率，就能判断系统故障状态。附图说明图1是实施例1中液体火箭发动机启动阶段工作过程中各关键节点的Petri网模型。图2是实施例1中故障变迁t10对应的QBPN模型。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。需要说明的是，在正式说明本专利技术实施例之前，有两点需要提前说明，第一本专利技术是对于液体火箭发动机启动阶段的可能出现的故障进行预测，液体火箭发动机启动阶段工作原理为高压氦气通入推进剂储箱挤压推进剂进入下游燃烧室混合点火产生推力，电磁阀未打开或燃料加注阀故障，则说明液体火箭发动机在启动阶段发生了故障。第二，介绍与本专利技术相关的两种算法，第一种算法是：1、基于SCG的状态估计算法输入：POPN模型及可能故障变迁和可观变迁序列σo；输出：系统的故障状步骤1、初始化根节点C0标记node(C0)＝E，对应标识M0，ti∈A(M0)；步骤2、当存在一个标记node(Ck)＝E的节点则选择一新节点Ck标记node(Ck)＝E，遍历任意ti∈A(M0)，ti在标识Mk触发Mq＝Mk+C(·,ti)建立新节点Cq标识为Mq，如果存在和Cq相同节点，则Cq节点标记node(Ck)＝Z否则Cq节点标记node(Ck)＝E；步骤3、用标记为E的节点构建SCG图；步骤4、令所有满足S(σo)路径放入集合；步骤5、遍历i＝1,…,r，如果则如果则如果则其中，为每个集合分配一个故障状态；所有满足可观变迁序列σo的路径放在集合S(σ0)；Pre∈(N)m×n、Post∈(N)m×n分别为变迁的前向弧矩阵和后向弧矩阵，C＝Post-Pre为POPN的关联矩阵，维数为m×n(N为非负整数集)，M0为初始标识；ψ为量子贝叶斯表；ti∈A(M0)表示从M0状态到其他状态的路径集合。在该算法中，步骤1，步骤2，步骤3为构建SCG图的过程，步骤4为遍历所有满足可观变迁序列σo的路径放在集合S(σ0)中，步骤5为根据集合S(σ0)中路径是否包含故障变迁判断系统故障状态。2、基于QBPN模型的故障诊断算法输入：POPN模型、量子贝叶斯概率表ψ和可能故障变迁输出：系统故障概率步骤1、初始化步骤2、遍历当则i＜j≤k,tj∈Tu，复制tk与之间所有变迁tj、库所P和弧到QBPN，当则k＜j≤i,tj∈Tu，复制与tk之间所有变迁tj、库所P和弧到QBPN，步骤3、如果tk没有发生权则删除tk及其输入输出弧；步骤4、遍历ti∈TQB，根据量子概率表ψ对ti设置P(ti)和量子概率振幅根据前向路径可观变迁状态计算步骤5、遍历令计算系统故障概率其中，变迁ti在标识M处点火，记为M[ti>，σ＝t0t1t2…th∈T是变迁序列集合，对应标识为M0[t0>M1[t1>…[th>Mh+1，简记为M0[σ>Mh+1，th为最后触发的变迁；T是n维的变迁集，To是可观变迁集，Tu是不可观变迁集；量子贝叶斯Petri网定义为其中为QBPN中可观变迁集，为QBPN中不可观变迁集，且ψ是ti∈TQB的量子概率表。在该算法中：步骤1，步骤2，步骤3和步骤4为在POPN模型基础上构建QBPN模型的过程，其中步骤2为根据故障变迁所在序列分别向前向后构建含有故障变迁的关系网络，步骤3为删除没有发生权的变迁，步骤4为给定变迁触发的量子概率幅值，并计算故障变迁的先验概率，步骤5为根据后向路径可观变迁触发状态修正故障变迁触发概率。QBPN模型诊断系统故障状态，首先根据前向路径中可观变迁触发状态计算故障变迁触发概率，再利用QB函数根据后向路径可本文档来自技高网...

【技术保护点】
1.一种液体火箭发动机启动阶段的Petri网量子贝叶斯故障诊断方法，其特征在于，包括如下步骤：/n步骤S1、根据液体火箭发动机启动阶段的工作状态建立液体火箭发动机启动阶段的Petri网模型；/n步骤S2、构建Petri网模型的SCG图，遍历SCG图中满足可观变迁触发序列的路径，并且判断是否包含故障变迁，初步估计液体火箭发动机的故障状态；/n步骤S3、建立QBPN模型，通过QBPN模型估计液体火箭发动机的故障概率。/n

【技术特征摘要】
1.一种液体火箭发动机启动阶段的Petri网量子贝叶斯故障诊断方法，其特征在于，包括如下步骤：
步骤S1、根据液体火箭发动机启动阶段的工作状态建立液体火箭发动机启动阶段的Petri网模型；
步骤S2、构建Petri网模型的SCG图，遍历SCG图中满足可观变迁触发序列的路径，并且判断是否包含故障变迁，初步估计液体火箭发动机的故障状态；
步骤S3、建立QBPN模型，通过QBPN模型估计液体...

【专利技术属性】
技术研发人员：张信哲，刘久富，汪恒宇，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏;32