CN109034225B - 一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法 - Google Patents

一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法 Download PDF

Info

Publication number
CN109034225B
CN109034225B CN201810778887.5A CN201810778887A CN109034225B CN 109034225 B CN109034225 B CN 109034225B CN 201810778887 A CN201810778887 A CN 201810778887A CN 109034225 B CN109034225 B CN 109034225B
Authority
CN
China
Prior art keywords
interval
sample
gray
random
parameter
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
Application number
CN201810778887.5A
Other languages
English (en)
Other versions
CN109034225A (zh
Inventor
方圣恩
陈杉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuzhou University
Original Assignee
Fuzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fuzhou University filed Critical Fuzhou University
Priority to CN201810778887.5A priority Critical patent/CN109034225B/zh
Publication of CN109034225A publication Critical patent/CN109034225A/zh
Application granted granted Critical
Publication of CN109034225B publication Critical patent/CN109034225B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/29—Graphical models, e.g. Bayesian networks
    • G06F18/295—Markov models or related models, e.g. semi-Markov models; Markov random fields; Networks embedding Markov models

Landscapes

  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analogue/Digital Conversion (AREA)
  • Complex Calculations (AREA)

Abstract

本发明涉及一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法，首先将结构随机参数表示为连续区间灰数，并求得该区间灰数的核与扩散度；接着基于区间灰数的核与扩散度构造随机参数的灰色先验分布，不断进行样本抽取与判别，获得初始灰色区间随机样本集；最后对步骤S2获得的样本集进行不断扰动和更新，并将所得样本的上下界置为更新后的区间灰数上下界；当区间宽度小于预设值时，终止迭代，以该区间灰数作为结构不确定性参数修正和估计的依据。本发明能够有效有效提高随机参数的估计精度。

Description

一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估 计方法
技术领域
本发明涉及工程结构健康监测中的模型修正技术领域，特别是一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法。
背景技术
结构健康监测中通常假定损伤会导致结构局部刚度发生变化，因此利用模型修正方法估计结构当前刚度参数是损伤识别的重要途径之一。近年来随着计算机技术的不断发展，模型修正技术[1,2]在结构健康监测领域得到了不断发展和应用。传统的确定性模型修正技术难以考虑实际结构边界条件未知、材料离散性与观测噪声等不确定因素的影响，实用性不佳，使得考虑参数不确定性的模型修正方法逐渐受到重视。
当前不确定性参数识别方法主要基于贝叶斯模型修正、随机有限元模型修正和统计模式识别[3]。相比于另外两种方法，贝叶斯模型修正能够充分利用结构的历史数据和当前实测数据信息，不断修正和更新待识别参数的先验概率分布，具备在线监测的独特优势[4]。贝叶斯模型修正时通常将随机参数的先验和后验概率假定为某一常见分布如正态分布和均匀分布[5-7]，然而实际工程中由于测试数据往往不充分，很多时候只能知道参数的大致取值范围，其概率分布形式不得而知，使得基于概率统计方法分析得到的参数后验概率分布估计值得商榷。因此，将结构随机参数考虑为存在一定取值范围的灰数，并将此类灰色随机变量纳入已有的贝叶斯模型修正过程更符合实际情况。
贝叶斯模型修正首先根据工程师的经验和历史数据大致估计结构随机参数的先验概率分布情况，然后在获取当前实测数据的基础上，通过贝叶斯公式不断修正更新先验分布，最终得到随机参数的后验概率分布估计。可见贝叶斯模型修正的本质是实现随机参数先验分布向后验分布的转化，但当前贝叶斯推断过程随机参数的选取是人为假定的(比如假定随机参数服从某一典型正态分布)，若真实分布与假定的不同，可能会得到错误的修正结果[8]。从事物认知的角度来看，表征结构损伤的随机参数本质上可看作部分已知、部分未知的灰色量。一方面可以根据结构的工作状态大致推知损伤导致的随机参数变动范围，另一方面由于主观认知不足及实测数据欠缺，需要收集更多信息进一步判断参数值的大小(损伤程度)。这种对随机参数存在一定程度认识但又无法确知参数大小的情况与灰色系统的概念类似，同时贝叶斯模型修正过程与灰色系统的白化过程也有相似之处。因此，将结构随机参数认为是灰性的，即将灰色随机变量与贝叶斯模型修正相结合，能更加客观地反映实际工程结构随机参数的演化情况。
其中参考文献如下：
[1]M.I.Friswell,J.E.Mottershead,Finite Element Model Updating inStructural Dynamics[M].Kluwer Academic Press,Dordrecht,1995。
[2]吴晓菊.结构有限元模型修正综述[J].《特种结构》,2009,26(1):39-45。
[3]宗周红,牛杰,王浩.基于模型确认的结构概率损伤识别方法研究进展[J].《土木工程学报》,2012(8):121-130。
[4]R.Rocchetta,M.Broggi,Q.Huchet,et al.On-line Bayesian modelupdating for structural health monitoring.Mechanical Systems&SignalProcessing,2018,103:174-195。
[5]J.L.Beck,L.S.Katafygiotis.Updating Models and TheirUncertainties.I:Bayesian Statistical Framework:Journal of EngineeringMechanics:(ASCE).Journal of Engineering Mechanics,1998,124(4):455-461。
[6]L.S.Katafygiotis,J.L.Beck.Updating Models and TheirUncertainties.II:Model Identifiability.Journal of Engineering Mechanics,1998,124(4):463-467。
[7]颜王吉,曹诗泽,任伟新.结构系统识别不确定性分析的Bayes方法及其进展[J].《应用数学和力学》,2017,38(1):44-59。
[8]张建新.基于贝叶斯方法的有限元模型修正研究[D].重庆大学,2014。
发明内容
有鉴于此，本发明的目的是提出一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法，能够有效提高随机参数的估计精度。
本发明采用以下方案实现：一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法，具体包括以下步骤：
步骤S1：将结构随机参数表示为连续区间灰数，并求得该区间灰数的核与扩散度；
步骤S2：基于区间灰数的核与扩散度构造随机参数的灰色先验分布，不断进行样本抽取与判别，获得初始灰色区间随机样本集；
步骤S3：对步骤S2获得的样本集进行不断扰动和更新，并将所得样本的上下界置为更新后的区间灰数上下界；当区间宽度小于预设值时，终止迭代，以该区间灰数作为结构不确定性参数修正和估计的依据。
其中本发明采用群体蒙特卡洛抽样生成所需的参数样本。
本发明在贝叶斯模型修正方法中引入灰色随机变量，在假定结构随机参数为区间灰数的基础上，以其真实值服从正态分布作为补充信息，然后基于群体蒙特卡洛抽样完成灰色随机参数的后验概率分布估计。该发明首先在贝叶斯模型修正框架内考虑随机参数的灰色性质，然后从参数取值分布随机性的角度考虑区间灰数的特点，使用正态分布随机函数完成区间灰数的数值覆盖。与传统贝叶斯模型修正方法所采用的马尔科夫链蒙特卡洛抽样不同，本发明采用群体蒙特卡洛抽样生成所需的参数样本，然后估计参数的灰色区间，能够有效提高随机参数的估计精度。
本发明结构随机参数的区间灰数表示；使用正态分布随机函数完成区间灰数的数值覆盖，将区间灰数转化为等效的实数集合；运用群体蒙特卡洛抽样生成所需的参数样本，最终实现参数灰色区间估计。具体阐述如下：
对于结构随机参数θ，给定命题P，p(x)表示命题信息。在无法确知命题真实取值的情况下，可以给定该命题的可能取值区间I。若θ在区间I内的取值连续，那么命题P可理解为p(x)意义下的连续型灰数
Figure GDA0002668913770000041
这里命题信息p(x)表示灰数
Figure GDA0002668913770000042
在x处取值的可能性。令参数在该命题下的真值为x0，可见灰数
Figure GDA0002668913770000043
描述了一个含有命题真值x0的数值区间I，而且
Figure GDA0002668913770000044
的灰色程度(即不确定程度)依赖于命题信息p(x)。
若θ取值区间
Figure GDA0002668913770000045
其中
Figure GDA0002668913770000046
记连续型区间灰数为
Figure GDA0002668913770000047
p；若参数真值x0的取值可能性p(x)服从正态分布，则称
Figure GDA0002668913770000048
p为正态分布区间灰数，简记为
Figure GDA0002668913770000049
本发明讨论的均为正态分布区间灰数。
基于正态分布的“3σ”原则定义区间灰数的核与信息扩散度。区间灰数的核：
Figure GDA0002668913770000051
区间灰数的信息扩散度：
Figure GDA0002668913770000052
于是有：
Figure GDA0002668913770000053
基于上述假定，区间灰数
Figure GDA0002668913770000054
真值x0的取值可能性服从正态分布
Figure GDA0002668913770000055
那么x0落在
Figure GDA0002668913770000056
内的概率为99.73％，基本实现了对灰数
Figure GDA0002668913770000057
的数值覆盖。
由于区间灰数本身无法做一般概率统计意义下的处理，因此需要将其转化为信息等效的实数集合。由前述可知，使用蒙特卡洛抽样产生的随机样本有99.73％都将落在以
Figure GDA0002668913770000058
为中心的3σ区间内，可以满足信息等效条件。另外，随机抽样获得的样本最大与最小值能够对应区间灰数的上下界信息，也满足了数据完整性原则。目前贝叶斯模型修正通常使用马尔科夫链蒙特卡洛抽样来完成参数收敛过程，但该方法存在陷入低概率区域抽样的风险，为此本发明采用群体蒙特卡洛抽样方法，对真值服从正态分布的区间灰数进行大量样本模拟，并将该过程与贝叶斯模型修正有机结合实现参数概率分布的更新。
进一步地，步骤S1具体包括以下步骤：
步骤S11：生成初始灰色区间随机样本集X1；
步骤S12：设结构随机参数θ的初始值为连续型区间灰数
Figure GDA0002668913770000059
求得该区间灰数的核与扩散度：
Figure GDA00026689137700000510
Figure GDA0002668913770000061
式中，
Figure GDA0002668913770000062
为区间灰数的核的初始值，σ0为区间灰数的扩散度的初始值。
进一步地，步骤S2具体包括以下步骤：
步骤S21：定义正态分布
Figure GDA0002668913770000063
为参数初始概率分布π(θ)，从π(θ)中随机抽取一个参数样本
Figure GDA0002668913770000064
将
Figure GDA0002668913770000065
代入结构数值模型求得仿真样本
Figure GDA0002668913770000066
步骤S22：构造并计算仿真样本
Figure GDA0002668913770000067
与实测样本D1之间的误差函数
Figure GDA0002668913770000068
对于给定容差ε1，若
Figure GDA0002668913770000069
则重新抽取
Figure GDA00026689137700000610
否则保存该样本并进行下一个样本的抽取和判别，直到生成足够数目的参数样本。
进一步地，步骤S3具体包括以下步骤：
步骤S31：将参数样本的最小值a 1和最大值
Figure GDA00026689137700000611
置为更新后的区间灰数上下界，即
Figure GDA00026689137700000612
并给所有样本设置相同的抽样权重w1＝1/N，其中N为样本数；计算当前迭代步t＝1时该样本集的方差
Figure GDA00026689137700000613
样本方差
Figure GDA00026689137700000614
定义为当前迭代步样本集方差的2倍：
Figure GDA00026689137700000615
步骤S32：后续迭代步t≥2时，对样本集Xt中的样本进行扰动更新：从上一迭代步t-1的样本集Xt-1里结合权重随机抽取样本
Figure GDA00026689137700000616
并用上一迭代步t-1的样本集方差
Figure GDA00026689137700000617
扰动
Figure GDA00026689137700000618
使其成为
Figure GDA00026689137700000619
具体操作上可以通过从
Figure GDA00026689137700000620
中抽样来实现；
步骤S33：将
Figure GDA00026689137700000621
代入结构数值模型求得仿真样本
Figure GDA00026689137700000622
并判断仿真样本与实测样本之间的误差是否小于预设的容差，若是，保存该样本
Figure GDA00026689137700000623
并计算其权重值，进入步骤S34，否则返回步骤S32；其中样本权重由当前样本的概率密度函数、上一迭代步得到的各样本的权重以及一个高斯概率转移核联合求得；
步骤S34：在样本集中所有样本都更新之后结束当前迭代步；每一步迭代之后计算样本集的方差，获得该迭代步内参数样本的上下界a t和
Figure GDA0002668913770000071
将其置换为更新后的区间灰数上下界，即
Figure GDA0002668913770000072
步骤S35：判断区间宽度是否小于预设的阈值，若是，则终止迭代，得到更新完毕的参数区间灰数
Figure GDA0002668913770000073
否则，返回步骤S32。
考虑到实际结构边界条件未知、材料离散性与观测噪声等不确定性因素的影响，模型修正和参数估计过程中考虑参数的不确定性很有必要。传统的贝叶斯模型修正方法通常将随机参数的先验和后验概率假定为某一常见分布，然而由于实测数据的匮乏，很多时候只能知道参数的大致取值范围，想要基于统计数据获得参数的概率分布函数并不现实，那么基于概率统计方法分析得到的参数后验概率分布估计值得商榷。因此，在“少数据、贫信息”的情况下，将结构随机参数考虑为存在一定取值范围的灰数，并将此类灰色随机变量纳入已有的贝叶斯模型修正过程更符合实际情况。
与现有技术相比，本发明有以下有益效果：
1、将灰色随机变量纳入贝叶斯模型修正框架内，能够客观反映实际工程结构随机参数的灰色性质；
2、使用正态分布随机函数完成区间灰数的数值覆盖，实现了从灰数域向实数域的等信息转换；
3、基于群体蒙特卡洛抽样实现参数后验分布样本的生成并进行参数的灰色区间估计，相较于传统的马尔科夫链蒙特卡洛抽样能够有效提高随机参数的估计精度。
附图说明
图1为本发明的方法流程示意图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出，以下详细说明都是示例性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
如图1所示，本实施例提供了一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法，具体包括以下步骤：
步骤S1：将结构随机参数表示为连续区间灰数，并求得该区间灰数的核与扩散度；
步骤S2：基于区间灰数的核与扩散度构造随机参数的灰色先验分布，不断进行样本抽取与判别，获得初始灰色区间随机样本集；
步骤S3：对步骤S2获得的样本集进行不断扰动和更新，并将所得样本的上下界置为更新后的区间灰数上下界；当区间宽度小于预设值时，终止迭代，以该区间灰数作为结构不确定性参数修正和估计的依据。
其中本实施例采用群体蒙特卡洛抽样生成所需的参数样本。
本实施例在贝叶斯模型修正方法中引入灰色随机变量，在假定结构随机参数为区间灰数的基础上，以其真实值服从正态分布作为补充信息，然后基于群体蒙特卡洛抽样完成灰色随机参数的后验概率分布估计。该发明首先在贝叶斯模型修正框架内考虑随机参数的灰色性质，然后从参数取值分布随机性的角度考虑区间灰数的特点，使用正态分布随机函数完成区间灰数的数值覆盖。与传统贝叶斯模型修正方法所采用的马尔科夫链蒙特卡洛抽样不同，本发明采用群体蒙特卡洛抽样生成所需的参数样本，然后估计参数的灰色区间，能够有效提高随机参数的估计精度。
本实施例结构随机参数的区间灰数表示；使用正态分布随机函数完成区间灰数的数值覆盖，将区间灰数转化为等效的实数集合；运用群体蒙特卡洛抽样生成所需的参数样本，最终实现参数灰色区间估计。具体阐述如下：
对于结构随机参数θ，给定命题P，p(x)表示命题信息。在无法确知命题真实取值的情况下，可以给定该命题的可能取值区间I。若θ在区间I内的取值连续，那么命题P可理解为p(x)意义下的连续型灰数
Figure GDA0002668913770000091
这里命题信息p(x)表示灰数
Figure GDA0002668913770000092
在x处取值的可能性。令参数在该命题下的真值为x0，可见灰数
Figure GDA0002668913770000093
描述了一个含有命题真值x0的数值区间I，而且
Figure GDA0002668913770000094
的灰色程度(即不确定程度)依赖于命题信息p(x)。
若θ取值区间
Figure GDA0002668913770000095
其中
Figure GDA0002668913770000096
记连续型区间灰数为
Figure GDA0002668913770000097
p；若参数真值x0的取值可能性p(x)服从正态分布，则称
Figure GDA0002668913770000098
p为正态分布区间灰数，简记为
Figure GDA0002668913770000101
本发明讨论的均为正态分布区间灰数。
基于正态分布的“3σ”原则定义区间灰数的核与信息扩散度。区间灰数的核：
Figure GDA0002668913770000102
区间灰数的信息扩散度：
Figure GDA0002668913770000103
于是有：
Figure GDA0002668913770000104
基于上述假定，区间灰数
Figure GDA0002668913770000105
真值x0的取值可能性服从正态分布
Figure GDA0002668913770000106
那么x0落在
Figure GDA0002668913770000107
内的概率为99.73％，基本实现了对灰数
Figure GDA0002668913770000108
的数值覆盖。
由于区间灰数本身无法做一般概率统计意义下的处理，因此需要将其转化为信息等效的实数集合。由前述可知，使用蒙特卡洛抽样产生的随机样本有99.73％都将落在以
Figure GDA0002668913770000109
为中心的3σ区间内，可以满足信息等效条件。另外，随机抽样获得的样本最大与最小值能够对应区间灰数的上下界信息，也满足了数据完整性原则。目前贝叶斯模型修正通常使用马尔科夫链蒙特卡洛抽样来完成参数收敛过程，但该方法存在陷入低概率区域抽样的风险，为此本发明采用群体蒙特卡洛抽样方法，对真值服从正态分布的区间灰数进行大量样本模拟，并将该过程与贝叶斯模型修正有机结合实现参数概率分布的更新。
在本实施例中，步骤S1具体包括以下步骤：
步骤S11：生成初始灰色区间随机样本集X1；
步骤S12：设结构随机参数θ的初始值为连续型区间灰数
Figure GDA00026689137700001010
求得该区间灰数的核与扩散度：
Figure GDA0002668913770000111
Figure GDA0002668913770000112
式中，
Figure GDA0002668913770000113
为区间灰数的核的初始值，σ0为区间灰数的扩散度的初始值。
在本实施例中，步骤S2具体包括以下步骤：
步骤S21：定义正态分布
Figure GDA0002668913770000114
为参数初始概率分布π(θ)，从π(θ)中随机抽取一个参数样本
Figure GDA0002668913770000115
将
Figure GDA0002668913770000116
代入结构数值模型求得仿真样本
Figure GDA0002668913770000117
步骤S22：构造并计算仿真样本
Figure GDA0002668913770000118
与实测样本D1之间的误差函数
Figure GDA0002668913770000119
对于给定容差ε1，若
Figure GDA00026689137700001110
则重新抽取
Figure GDA00026689137700001111
否则保存该样本并进行下一个样本的抽取和判别，直到生成足够数目的参数样本。
在本实施例中，步骤S3具体包括以下步骤：
步骤S31：将参数样本的最小值a 1和最大值
Figure GDA00026689137700001112
置为更新后的区间灰数上下界，即
Figure GDA00026689137700001113
并给所有样本设置相同的抽样权重w1＝1/N，其中N为样本数；计算当前迭代步t＝1时该样本集的方差
Figure GDA00026689137700001114
样本方差
Figure GDA00026689137700001115
定义为当前迭代步样本集方差的2倍：
Figure GDA00026689137700001116
步骤S32：后续迭代步t≥2时，对样本集Xt中的样本进行扰动更新：从上一迭代步t-1的样本集Xt-1里结合权重随机抽取样本
Figure GDA00026689137700001117
并用上一迭代步t-1的样本集方差
Figure GDA00026689137700001118
扰动
Figure GDA00026689137700001119
使其成为
Figure GDA00026689137700001120
具体操作上可以通过从
Figure GDA00026689137700001121
中抽样来实现；
步骤S33：将
Figure GDA00026689137700001122
代入结构数值模型求得仿真样本
Figure GDA00026689137700001123
并判断仿真样本与实测样本之间的误差是否小于预设的容差，若是，保存该样本
Figure GDA00026689137700001124
并计算其权重值，进入步骤S34，否则返回步骤S32；其中样本权重由当前样本的概率密度函数、上一迭代步得到的各样本的权重以及一个高斯概率转移核联合求得；
步骤S34：在样本集中所有样本都更新之后结束当前迭代步；每一步迭代之后计算样本集的方差，获得该迭代步内参数样本的上下界a t和
Figure GDA0002668913770000121
将其置换为更新后的区间灰数上下界，即
Figure GDA0002668913770000122
步骤S35：判断区间宽度是否小于预设的阈值，若是，则终止迭代，得到更新完毕的参数区间灰数
Figure GDA0002668913770000123
否则，返回步骤S32。
考虑到实际结构边界条件未知、材料离散性与观测噪声等不确定性因素的影响，模型修正和参数估计过程中考虑参数的不确定性很有必要。传统的贝叶斯模型修正方法通常将随机参数的先验和后验概率假定为某一常见分布，然而由于实测数据的匮乏，很多时候只能知道参数的大致取值范围，想要基于统计数据获得参数的概率分布函数并不现实，那么基于概率统计方法分析得到的参数后验概率分布估计值得商榷。因此，在“少数据、贫信息”的情况下，将结构随机参数考虑为存在一定取值范围的灰数，并将此类灰色随机变量纳入已有的贝叶斯模型修正过程更符合实际情况。
以上所述仅为本发明的较佳实施例，凡依本发明申请专利范围所做的均等变化与修饰，皆应属本发明的涵盖范围。

Claims (1)

1.一种结合随机变量灰性和贝叶斯模型修正的工程结构不确定性参数估计方法，其特征在于：包括以下步骤：
步骤S1：将结构随机参数表示为连续区间灰数，并求得该区间灰数的核与扩散度；
步骤S2：基于区间灰数的核与扩散度构造结构随机参数的灰色先验分布，不断进行样本抽取与判别，获得初始灰色区间随机样本集；
步骤S3：对步骤S2获得的样本集进行不断扰动和更新，并将所得样本的上下界置为更新后的区间灰数上下界；当区间宽度小于预设值时，终止迭代，以该区间灰数作为该结构随机参数修正和估计的依据；
其中，步骤S1具体包括以下步骤：
步骤S11：生成初始灰色区间随机样本集X1；
步骤S12：设结构随机参数θ的初始值为连续型区间灰数
Figure FDA0002668913760000011
求得该区间灰数的核与扩散度：
Figure FDA0002668913760000012
Figure FDA0002668913760000013
式中，
Figure FDA0002668913760000014
为区间灰数的核的初始值，σ0为区间灰数的扩散度的初始值；
其中，步骤S2具体包括以下步骤：
步骤S21：定义正态分布
Figure FDA0002668913760000015
为参数初始概率分布π(θ)，从π(θ)中随机抽取一个参数样本
Figure FDA0002668913760000016
将
Figure FDA0002668913760000017
代入结构数值模型求得仿真样本
Figure FDA0002668913760000018
步骤S22：构造并计算仿真样本
Figure FDA0002668913760000021
与实测样本D1之间的误差函数
Figure FDA0002668913760000022
对于给定容差ε1，若
Figure FDA0002668913760000023
则重新抽取
Figure FDA0002668913760000024
否则保存该样本并进行下一个样本的抽取和判别，直到生成足够数目的参数样本；
其中，步骤S3具体包括以下步骤：
步骤S31：将参数样本的最小值a 1和最大值
Figure FDA0002668913760000025
置为更新后的区间灰数上下界，即
Figure FDA0002668913760000026
并给所有样本设置相同的抽样权重w1＝1/N，其中N为样本数；计算当前迭代步t＝1时该样本集的方差
Figure FDA0002668913760000027
样本方差
Figure FDA0002668913760000028
定义为当前迭代步样本集方差的2倍：
Figure FDA0002668913760000029
步骤S32：后续迭代步t≥2时，对样本集Xt中的样本进行扰动更新：从上一迭代步t-1的样本集Xt-1里结合权重随机抽取样本
Figure FDA00026689137600000210
并用上一迭代步t-1的样本集方差
Figure FDA00026689137600000211
扰动
Figure FDA00026689137600000212
使其成为
Figure FDA00026689137600000213
步骤S33：将
Figure FDA00026689137600000214
代入结构数值模型求得仿真样本
Figure FDA00026689137600000215
并判断仿真样本与实测样本之间的误差是否小于预设的容差，若是，保存该样本
Figure FDA00026689137600000216
并计算其权重值，进入步骤S34，否则返回步骤S32；
步骤S34：在样本集中所有样本都更新之后结束当前迭代步；每一步迭代之后计算样本集的方差，获得该迭代步内参数样本的上下界a t和
Figure FDA00026689137600000217
将其置换为更新后的区间灰数上下界，即
Figure FDA00026689137600000218
步骤S35：判断区间宽度是否小于预设的阈值，若是，则终止迭代，得到更新完毕的参数区间灰数
Figure FDA00026689137600000219
否则，返回步骤S32。
CN201810778887.5A 2018-07-16 2018-07-16 一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法 Active CN109034225B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810778887.5A CN109034225B (zh) 2018-07-16 2018-07-16 一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810778887.5A CN109034225B (zh) 2018-07-16 2018-07-16 一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法

Publications (2)

Publication Number Publication Date
CN109034225A CN109034225A (zh) 2018-12-18
CN109034225B true CN109034225B (zh) 2021-01-29

Family

ID=64643093

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810778887.5A Active CN109034225B (zh) 2018-07-16 2018-07-16 一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法

Country Status (1)

Country Link
CN (1) CN109034225B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109885893B (zh) * 2019-01-24 2021-07-09 北京航空航天大学 基于不对称信息的贝叶斯加速验收抽样方案的设计方法
CN110008568A (zh) * 2019-03-29 2019-07-12 中国人民解放军战略支援部队航天工程大学 小样本下装备平均修复时间的非统计估计模型
CN110009255B (zh) * 2019-04-16 2021-07-20 西南交通大学 基于三参数区间泛灰数的铁路车站能力表示方法
CN110263465B (zh) * 2019-06-26 2023-10-17 上海工程技术大学 一种用于二阶阶梯杆结构不确定性的改进灰度方法
CN114117759A (zh) * 2021-11-12 2022-03-01 大连海事大学 基于非参数模型的大型船舶轴系固有频率不确定性分析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101504717A (zh) * 2008-07-28 2009-08-12 上海高德威智能交通系统有限公司 特征区域的定位方法、车身深浅色与车身颜色的识别方法
CN101799369A (zh) * 2010-01-27 2010-08-11 北京信息科技大学 一种动力设备故障监测预报中的灰色建模方法
CN107730084A (zh) * 2017-09-18 2018-02-23 杭州安脉盛智能技术有限公司 基于灰色预测和风险评估的变压器检修决策方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8201087B2 (en) * 2007-02-01 2012-06-12 Tegic Communications, Inc. Spell-check for a keyboard system with automatic correction

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101504717A (zh) * 2008-07-28 2009-08-12 上海高德威智能交通系统有限公司 特征区域的定位方法、车身深浅色与车身颜色的识别方法
CN101799369A (zh) * 2010-01-27 2010-08-11 北京信息科技大学 一种动力设备故障监测预报中的灰色建模方法
CN107730084A (zh) * 2017-09-18 2018-02-23 杭州安脉盛智能技术有限公司 基于灰色预测和风险评估的变压器检修决策方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Solution of nonlinear programming problem with interval grey parameters;Shu-Hua Mao et al;《2015 IEEE International Conference on Grey Systems and Intelligent Services (GSIS)》;20151026;91-96页 *
随机信息中系统参数的灰色贝叶斯区间估计研究;李勇等;《西南师范大学学报（自然科学版）》;20111031;第36卷(第5期);59-61页 *

Also Published As

Publication number Publication date
CN109034225A (zh) 2018-12-18

Similar Documents

Publication Publication Date Title
CN109034225B (zh) 一种结合随机变量灰性和贝叶斯模型修正的不确定性参数估计方法
CN108304679A (zh) 一种自适应可靠性分析方法
Li et al. Uncertainty analysis in ecological studies
CN110427654B (zh) 一种基于敏感状态的滑坡预测模型构建方法及系统
CN113239022B (zh) 医疗诊断缺失数据补全方法及补全装置、电子设备、介质
CN109145399B (zh) 一种基于改进的粒子滤波算法的疲劳裂纹扩展预测方法
CN115935834A (zh) 一种基于深度自回归网络和持续学习策略的历史拟合方法
CN115270239A (zh) 基于动力特性和智能算法响应面法的桥梁可靠性预测方法
CN115982141A (zh) 一种针对时序数据预测的特征优化方法
CN115455793A (zh) 基于多尺度模型修正的高层结构复杂构件受力分析方法
Zhang et al. Predicting growth and interaction of multiple cracks in structural systems using Dynamic Bayesian Networks
CN116187153B (zh) 基于层次贝叶斯的水工结构数字孪生模型更新方法
CN116629010A (zh) 一种基于随机过程的退化模型确认及其试验设计方法
Galioto et al. Likelihood-based generalization of Markov parameter estimation and multiple shooting objectives in system identification
Liu Leave-group-out cross-validation for latent Gaussian models
KR101549349B1 (ko) 데이터 기반 모델의 건물 에너지 모델링 방법과 이를 이용한 에너지 최적관리와 제어 방법
CN114580151A (zh) 一种基于灰色线性回归-马尔科夫链模型的需水预测方法
He et al. Singular-Value Decomposition Feature-Extraction Method for Cost-Performance Prediction
Li Sensitivity analysis and uncertainty integration for system diagnosis and prognosis
CN115048868B (zh) 基于时序神经网络的动态测量系统不确定度的评定方法
Wagener et al. Stochastic formulation of a conceptual hydrological model
CN117131465B (zh) 单管塔损伤识别方法、装置、电子设备及可读存储介质
Vojinovic A complementary modelling approach to manage uncertainty of computationally expensive models
WO2020071430A1 (ja) 情報処理装置、情報処理システム、情報処理方法及びプログラムが格納された非一時的なコンピュータ可読媒体
Du Combining statistical methods with dynamical insight to improve nonlinear estimation

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