CN112884348B - 基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法 - Google Patents
基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法 Download PDFInfo
- Publication number
- CN112884348B CN112884348B CN202110268715.5A CN202110268715A CN112884348B CN 112884348 B CN112884348 B CN 112884348B CN 202110268715 A CN202110268715 A CN 202110268715A CN 112884348 B CN112884348 B CN 112884348B
- Authority
- CN
- China
- Prior art keywords
- deviation
- node
- nodes
- source
- network
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06395—Quality analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Abstract
本发明公开了基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法,主要涉及航天火工品生产领域;包括步骤:S1、提取偏差源、观测特征,确定输入、输出节点;S2、分析变量间依赖关系,确定初始网络结构;S3、确定批次状态转移概率,建立动态贝叶斯网络;S4、多次随机抽样,扩充数据集;S5、构建偏差敏感度矩阵;S6、对初始网络进行结构学习和参数学习并完成测点优化;本发明改进了传统诊断方法中经济成本高、偏差源难以定位、过度依赖专家经验的不足,提高偏差诊断的精度和可靠性,减少人力物力投入。
Description
技术领域
本发明涉及航天火工品生产领域,具体是基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法。
背景技术
随着我国国民经济水平的不断提升和航天事业的发展,航天火工品的需求也不断增加,更加高效、精确的火工品生产加工成为了亟待解决的问题。由于航天火工品制造的特殊性,火工品制造过程在保证稳定的同时,在安全性、精度等方面都高于其它行业的要求。在航天火工品起爆器的生产过程中,零件制造偏差、定位偏差、装配误差以及操作缺陷等误差源都会对起爆器制造偏差造成影响,进而对整个火工品质量乃至下游产品的功能产生很大的影响,因此起爆器质量控制技术受到学界与业界的高度重视。目前,国内航天起爆器的生产加工仍然存在一些问题:
1)、起爆器的生产处于半自动化状态,部分加工步骤需要由人工来完成生产,过多依赖于主观经验,因此产品质量的一致性难以得到保证。
2)、起爆器产品有小批量、多批次的特点,生产加工过程中参数众多、数据量大、规格复杂,观测数据存在波动性、传递性和耦合性,对产品偏差的定位和识别存在一定困难。
针对上述存在的问题,我们迫切需要找到一种分析方法使得产品在生产加工中质量的一致性能够得到保障。
发明内容
本发明的目的在于解决现有技术中存在的问题,提供一种基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法,它不依赖于具体偏差解析模型的偏差源诊断策略,其中应用了动态贝叶斯网络方法,可被应用于生产制造过程的偏差源诊断和质量一致性分析;改进了传统诊断方法中的经济成本高、偏差源难以定位、过度依赖专家经验的不足,提高了偏差诊断的精度和可靠性,减少了人力物力投入。
本发明为实现上述目的,通过以下技术方案实现:
基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法,包括步骤:
S1、提取偏差源、观测特征,确定输入、输出节点;
S2、分析变量间依赖关系,确定初始网络结构;
S3、确定批次状态转移概率,建立动态贝叶斯网络;
S4、多次随机抽样,扩充数据集;
S5、构建偏差敏感度矩阵;
S6、对初始网络进行结构学习和参数学习并完成测点优化。
进一步的,步骤S1中,确定输入节点包括零件偏差S0、定位偏差Vf、装配偏差Vw以及其他偏差源Vq,输出节点为关键产品特征观测点S1,S2…,Sn。
进一步的,步骤S2中,用网络节点表示偏差源节点和观测节点,用有向边表示偏差源节点对观测节点的因果关系,用概率表示输入与输出之间的关系强度,给出条件概率表用来量化上层节点对下层节点的不确定性影响,由此确定模型初始结构。
进一步的,步骤S3中,定义Pi、Pj用以表示依赖关系,Pi表示前一批次成功的前提下下一批次成功的概率,Pj表示前一批次失败的前提下下一批次失败的概率,由此构建生产制造偏差的动态贝叶斯网络。
进一步的,步骤S4中,收集检测数据,基于标准正态分布N(0,σ2),对比检测数据与均值状态之间的距离大小来实现检测数据的离散化,针对小数据集D,多次随机抽样取数据子集Dn=0.5D。
进一步的,步骤S5中,定义产品质量偏差与各级偏差之间的敏感度矩阵S,描绘偏差源向量{Vu}与总成偏差{Vw}之间的线性映射关系,基于区分度指标法选取高敏感度特征参数,计算节点间的互信息大小,由此分析节点之间的依赖关系。
进一步的,步骤S6中,将先验概率与样本数据融合,推理出偏差诊断的后验概率,并依据敏感度矩阵和互信息计算结果,对初始网络进行结构学习和参数学习并完成测点优化。
进一步的,在初始网络基础上,基于条件独立性检验,对动态贝叶斯网络进行更新的算法具体工作如下:
步骤①、设定学习循环次数上限u以及判定阈值ε,给出初始网络结构X,定义观测节点Sj的潜在偏差源节点集
步骤②、计算偏差源节点Fi与第j个观测节点Sj之间的互信息大小:
I(Fi;Sj)=H(Fi)-H(Fi|Sj);
步骤③、将前k个与Sj具有最大互信息的偏差源节点作为Sj的潜在父节点,
步骤④、条件判定,如果或当k=u时,确认添加从π(Sj)到Sj的有向边,删除初始结构中/>到节点Sj的有向边,同时更新网络参数;
步骤⑤、如果不满足条件,则k←k+1,进入下一轮学习;当步骤④中两个条件均无法满足时,增加Sj的潜在父节点集中偏差源节点数目,从步骤②开始重新检验,直到Sj与其他偏差源节点条件独立或者父节点数目达到u时,停止学习;
步骤⑥、依据以上流程依次对所有观测节点进行条件独立性检验,实现对网络结构和参数的更新,完成测点优化。
对比现有技术,本发明的有益效果在于:
1、本发明基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法,旨在改进当前起爆器生产加工偏差诊断过程过度依赖专家经验的不足,具有工程意义;
2、本发明针对火工品生产小数据集、多批次的特点,基于变量间依赖关系和状态转移概率分布,设计偏差源诊断方法,其在实际运用中有更高的精度和可靠性;
3、本发明针对实际系统偏差的波动性、传递性和耦合性,基于多源信息融合分析,在实际应用中具有不依赖具体、精确的偏差解析数学模型的优势。
附图说明
附图1是本发明的偏差源诊断方法实施流程图;
附图2是本发明的动态贝叶斯网络更新流程图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所限定的范围。
实施例:本发明公开了一种基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法,包括基于变量间依赖关系的模型建立,各参数耦合强度的表示,基于批次状态转移概率分布的构造,将质量诊断问题转化为依据历史数据和当前观测的后验估计。
具体的,如附图1所示,包括步骤:
S1、提取偏差源、观测特征,确定输入、输出节点;
S2、分析变量间依赖关系,确定初始网络结构;
S3、确定批次状态转移概率,建立动态贝叶斯网络;
S4、多次随机抽样,扩充数据集;
S5、构建偏差敏感度矩阵;
S6、对初始网络进行结构学习和参数学习并完成测点优化。
进一步的,步骤S1中,针对航天火工品起爆器生产制造过程中的工艺步骤和质量控制点,确定输入节点包括零件偏差S0、定位偏差Vf、装配偏差Vw以及其他偏差源Vq,输出节点为关键产品特征观测点S1,S2…,Sn。
进一步的,步骤S2中,用网络节点表示偏差源节点和观测节点,用有向边表示偏差源节点对观测节点的因果关系,用概率表示输入与输出之间的关系强度,给出条件概率表用来量化上层节点对下层节点的不确定性影响,由此确定模型初始结构。
进一步的,步骤S3中,确定批次间的状态转移概率,定义Pi、Pj用以表示依赖关系,Pi表示前一批次成功的前提下下一批次成功的概率,Pj表示前一批次失败的前提下下一批次失败的概率,由此构建生产制造偏差的动态贝叶斯网络。
进一步的,步骤S4中,收集检测数据,基于标准正态分布N(0,σ2),对比检测数据与均值状态之间的距离大小来实现检测数据的离散化,针对小数据集D,多次随机抽样取数据子集Dn=0.5D,以此扩充数据集的容量、增加数据的多样性,从而提高学习的精度。
进一步的,步骤S5中,定义产品质量偏差与各级偏差之间的敏感度矩阵S,描绘偏差源向量{Vu}与总成偏差{Vw}之间的线性映射关系,基于区分度指标法选取高敏感度特征参数,计算节点间的互信息大小,由此分析节点之间的依赖关系。
进一步的,步骤S6中,将先验概率与样本数据融合,推理出偏差诊断的后验概率,并依据敏感度矩阵和互信息计算结果,对初始网络进行结构学习和参数学习并完成测点优化。
进一步的,在初始网络基础上,基于条件独立性检验,对动态贝叶斯网络进行更新的算法具体工作如下:
步骤①、设定学习循环次数上限u以及判定阈值ε,给出初始网络结构X,定义观测节点Sj的潜在偏差源节点集
步骤②、计算偏差源节点Fi与第j个观测节点Sj之间的互信息大小:
I(Fi;Sj)=H(Fi)-H(Fi|Sj);
步骤③、将前k个与Sj具有最大互信息的偏差源节点作为Sj的潜在父节点,
步骤④、条件判定,如果或当k=u时,确认添加从π(Sj)到Sj的有向边,删除初始结构中/>到节点Sj的有向边,同时更新网络参数;
步骤⑤、如果不满足条件,则k←k+1,进入下一轮学习;当步骤④中两个条件均无法满足时,增加Sj的潜在父节点集中偏差源节点数目,从步骤②开始重新检验,直到Sj与其他偏差源节点条件独立或者父节点数目达到u时,停止学习;
步骤⑥、依据以上流程依次对所有观测节点进行条件独立性检验,实现对网络结构和参数的更新,完成测点优化。
本发明改进了传统诊断方法中的经济成本高、偏差源难以定位、过度依赖专家经验的不足,提高了偏差诊断的精度和可靠性,减少了人力物力投入。
Claims (3)
1.基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法,其特征在于,包括步骤:
S1、提取偏差源和观测特征,确定输入和输出节点;
S2、分析变量间依赖关系,确定初始网络结构;步骤S2中,用网络节点表示偏差源节点和观测节点,用有向边表示偏差源节点对观测节点的因果关系,用概率表示输入与输出之间的关系强度,给出条件概率表用来量化上层节点对下层节点的不确定性影响,由此确定模型初始结构;
S3、确定批次状态转移概率,建立动态贝叶斯网络;步骤S3中,定义Pi、Pj用以表示依赖关系,Pi表示前一批次成功的前提下下一批次成功的概率,Pj表示前一批次失败的前提下下一批次失败的概率,由此构建生产制造偏差的动态贝叶斯网络;
S4、多次随机抽样,扩充数据集;
S5、构建偏差敏感度矩阵;步骤S5中,定义产品质量偏差与各级偏差之间的敏感度矩阵S,描绘偏差源向量{Vu}与总成偏差{Vw}之间的线性映射关系,基于区分度指标法选取高敏感度特征参数,计算节点间的互信息大小,由此分析节点之间的依赖关系;
S6、对初始网络进行结构学习和参数学习并完成测点优化;步骤S6中,将先验概率与样本数据融合,推理出偏差诊断的后验概率,并依据敏感度矩阵和互信息计算结果,对初始网络进行结构学习和参数学习并完成测点优化;
在初始网络基础上,基于条件独立性检验,对动态贝叶斯网络进行更新的算法具体工作如下:
步骤①、设定学习循环次数上限u以及判定阈值ε,给出初始网络结构X,定义观测节点Sj的潜在偏差源节点集
步骤②、计算偏差源节点Fi与第j个观测节点Sj之间的互信息大小:
I(Fi;Sj)=H(Fi)-H(FiSj);
步骤③、将前k个与Sj具有最大互信息的偏差源节点作为Sj的潜在父节点,
步骤④、条件判定,如果或当k=u时,确认添加从π(Sj)到Sj的有向边,删除初始结构中/>到节点Sj的有向边,同时更新网络参数;
步骤⑤、如果不满足条件,则k←k+1,进入下一轮学习;当步骤④中两个条件均无法满足时,增加Sj的潜在父节点集中偏差源节点数目,从步骤②开始重新检验,直到Sj与其他偏差源节点条件独立或者父节点数目达到u时,停止学习;
步骤⑥、依据以上流程依次对所有观测节点进行条件独立性检验,实现对网络结构和参数的更新,完成测点优化。
2.根据权利要求1所述的基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法,其特征在于:步骤S1中,确定输入节点包括零件偏差S0、定位偏差Vf、装配偏差Vw以及其他偏差源Vq,输出节点为关键产品特征观测点S1,S2…,Sn。
3.根据权利要求1所述的基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法,其特征在于:步骤S4中,收集检测数据,基于标准正态分布N(0,σ2),对比检测数据与均值状态之间的距离大小来实现检测数据的离散化,针对小数据集D,多次随机抽样取数据子集Dn=0.5D。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110268715.5A CN112884348B (zh) | 2021-03-12 | 2021-03-12 | 基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110268715.5A CN112884348B (zh) | 2021-03-12 | 2021-03-12 | 基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112884348A CN112884348A (zh) | 2021-06-01 |
CN112884348B true CN112884348B (zh) | 2023-09-15 |
Family
ID=76041179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110268715.5A Active CN112884348B (zh) | 2021-03-12 | 2021-03-12 | 基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112884348B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113393146B (zh) * | 2021-06-25 | 2022-12-20 | 重庆大学 | 基于贝叶斯网络的多阶段柔性测试系统可靠性评估方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108134680A (zh) * | 2016-11-30 | 2018-06-08 | 中国科学院沈阳自动化研究所 | 一种基于贝叶斯网络的系统测量节点优化配置方法 |
CN111368888A (zh) * | 2020-02-25 | 2020-07-03 | 重庆邮电大学 | 基于深度动态贝叶斯网络的服务功能链故障诊断方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10110976B2 (en) * | 2016-10-10 | 2018-10-23 | Utilidata, Inc. | Systems and methods for system measurements integrity determination |
-
2021
- 2021-03-12 CN CN202110268715.5A patent/CN112884348B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108134680A (zh) * | 2016-11-30 | 2018-06-08 | 中国科学院沈阳自动化研究所 | 一种基于贝叶斯网络的系统测量节点优化配置方法 |
CN111368888A (zh) * | 2020-02-25 | 2020-07-03 | 重庆邮电大学 | 基于深度动态贝叶斯网络的服务功能链故障诊断方法 |
Non-Patent Citations (1)
Title |
---|
基于小数据集贝叶斯网络建模的偏差源诊断方法;刘银华等;《上海交通大学学报》;第46卷(第5期);第701-705页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112884348A (zh) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109255173B (zh) | 考虑区间不确定性的结构失效概率区间计算方法 | |
CN108664700B (zh) | 基于不确定数据包络分析的加速退化信息融合建模方法 | |
CN109992921B (zh) | 一种燃煤电厂锅炉热效率的在线软测量方法及系统 | |
TWI584134B (zh) | 製程異因分析方法與製程異因分析系統 | |
CN108875771A (zh) | 一种基于稀疏高斯伯努利受限玻尔兹曼机和循环神经网络的故障分类模型及方法 | |
CN109490072B (zh) | 一种土木工程建筑用检测系统及其检测方法 | |
CN111950627A (zh) | 一种多源信息融合方法及其应用 | |
CN112884348B (zh) | 基于动态贝叶斯网络的航天起爆器生产偏差源诊断方法 | |
CN116167640B (zh) | 一种lcp薄膜生产质量检测数据分析方法及系统 | |
CN114970688A (zh) | 基于LSTMAD算法和Hermite插值法的滑坡监测数据预处理方法 | |
CN106815428B (zh) | 一种基于智能优化算法的风洞天平校准数据处理方法 | |
CN108537249B (zh) | 一种密度峰值聚类的工业过程数据聚类方法 | |
CN108563201B (zh) | 一种dmaic驱动的小批量零件加工工艺过程质量改进方法 | |
TWI709054B (zh) | 預測模型的建立裝置、建立方法與產品品質監控系統 | |
CN110222814B (zh) | 基于嵌入式dcnn的乙烯裂解炉炉管重管识别方法 | |
CN110197199B (zh) | 嵌入式dcnn和边缘计算的裂解炉重管温度识别方法 | |
CN112329868A (zh) | 基于clara聚类的制造加工设备群能效状态评价方法 | |
TWI801973B (zh) | 錯誤因素之推定裝置及推定方法 | |
CN115081485B (zh) | 一种基于ai的漏磁内检测数据自动分析方法 | |
CN116451170A (zh) | 一种基于多源异构数据的配电网检测方法 | |
CN113312587A (zh) | 基于arima预测和回归预测的传感器采集数据缺失值处理方法 | |
CN114861810A (zh) | 煤气化装置工艺诊断方法及系统 | |
CN111507003A (zh) | 基于质量相关动态特性提取的脱丁烷塔关键变量预测方法 | |
CN117556366B (zh) | 基于数据筛选的数据异常检测系统及方法 | |
CN114295095B (zh) | 一种自由曲面检测最佳测点数确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |