一种基于特征提取和面积度量的动态仿真结果验证方法技术

本发明专利技术公开了一种基于特征提取和面积度量的动态仿真结果验证方法，属于计算机仿真技术领域。具体的，采用离散正交多项式提取时间序列的均值、斜率、曲率等特征，进而计算各个特征的概率分布面积差异，最后综合得到验证结果。本发明专利技术能解决具有不确定性的动态输出的仿真结果验证问题。同时，特征系数具有明确的物理含义，可根据正交多项式拟合效果和动态输出的实际物理含义选择待提取的特征系数。的实际物理含义选择待提取的特征系数。的实际物理含义选择待提取的特征系数。

【技术实现步骤摘要】
一种基于特征提取和面积度量的动态仿真结果验证方法


[0001]本专利技术涉及一种基于特征提取和面积度量的动态仿真结果验证方法，属于计算机仿真


技术介绍

[0002]仿真技术已经成为继理论研究和实验研究之后第三种认识、改造客观世界的重要手段，且由于其具有经济性、安全性、可重复性、无破坏性、不受气候条件限制、不受场地空间限制等优点已经广泛的应用于工业、农业、医学、经济、航天、航空等领域。随着仿真技术的广泛应用和快速发展，仿真结果的可信性也越来越受人们的关注。二十世纪五十年代末期，仿真模型可信度的问题首次被提出，并随之开展了一定的研究，出现了图灵测试、 TIC法、假设检验等一系列经典的仿真结果验证方法。进入二十一世纪，随着仿真系统复杂程度的增加，仿真输出也呈现多元性、不确定性、相关性等特点，以Oberkampf、Mahadevan 等为代表的学者相继提出了概率分布面积差异、区间估计、贝叶斯假设检验等仿真结果验证方法。
[0003]针对仿真系统的动态输出，此方面的研究相对较少。有学者将时间序列中每一时间点作为一个变量，将动态时间序列的仿真结果验证问题转化为多元静态变量的仿真结果验证问题，但由于时间序列的自相关性和高维度，此种方法的计算量较大且结果不够准确；此外，动态输出的仿真结果验证方法主要集中于采用特征提取方法提取时间序列的特征，进而对特征进行一致性分析，如基于概率主成分分析和贝叶斯因子相结合的仿真结果验证方法、基于KL变换(Karhunen
‑
Loe`ve Transformation)的动态输出验证方法等，但已有的方法中采用主成分分析或概率主成分分析等提取的特征物理含义不明确；国内外相关学者还提出了基于小波变换的模型验证方法、时间域上的可靠性准则等。

技术实现思路

[0004]本专利技术的目的是提出一种基于特征提取和面积度量的动态仿真结果验证方法，以解决已有的方法中采用主成分分析或概率主成分分析等提取的特征物理含义不明确、计算量较大且结果不够准确等问题。
[0005]一种基于特征提取和面积度量的动态仿真结果验证方法，假设动态输出的仿真时间序列表示为Y＝{y
i
}，i＝1,2,...,M，参考时间序列表示为R＝{r
i
}，i＝1,2,...,M，动态输出的仿真数据和参考数据包含多个时间序列。此动态仿真结果验证方法分为单验证点和多验证点两个部分，
[0006]在单验证点，动态输出仿真结果验证包括以下步骤：
[0007]S100、在单验证点x，得到仿真模型的动态输出为Y
l
，l＝1,2,...,N
s
，实际观测数据为R
k
， k＝1,2,...,N
r
；
[0008]S200、采用离散切比雪夫多项式拟合时间序列，得到仿真数据和参考数据的正交展开系数分别为{α
sl
}
q
和{α
rk
}
q
，q＝1,...,Q，其中Q为正交展开系数的个数，由此得到{α
sl
}
q
的累积分布函数为F
sq
(α)，{α
rk
}
q
的经验累积分布函数为F
rq
(α)；
[0009]S300、采用面积度量计算F
sq
(α)和F
rq
(α)的面积差异为d
q
，即并对面积差异归一化得到各个正交展开系数仿真数据和参考数据的一致性c
q
，面积差异的归一化计算公式为：其中x
max
和x
min
分别为对应的正交展开系数的仿真和参考数据中的最大值和最小值；
[0010]S400、综合各个正交展开系数的一致性，计算得到动态输出的验证结果其中ω
q
为各个正交展开系数的权重且
[0011]在多验证点，动态输出仿真结果验证包括以下步骤：
[0012]S500、在验证点{x1,x2,...,x
P
}，得到仿真模型的动态输出为实际观测数据为其中l＝1,2,...,N
s
，k＝1,2,...,N
r
，j＝1,2,...,P；
[0013]S600、采用离散切比雪夫多项式拟合时间序列，得到仿真数据和参考数据的正交展开系数分别为和q＝1,...,Q，其中Q为正交展开系数的个数；
[0014]S700、计算仿真数据的正交展开系数的累积分布函数{F
j
}
q
，根据概率积分变换，各个系数的累积分布函数{F
j
}
q
转换为均匀分布U(0,1)，同时，参考数据的正交展开系数根据对应系数的{F
j
}
q
转换为一系列的u值，转换关系为
[0015]S800、计算各个系数的与均匀分布U(0,1)的面积差异其中F(u)
q
为的经验累积分布函数，并对面积差异归一化得到各个正交展开系数仿真数据和参考数据的一致性c
q
，面积差异的归一化计算公式为：c
q
＝1
‑
2*d
q
；
[0016]S900、综合各个正交展开系数的一致性计算的得到动态输出的验证结果其中ω
q
为各个正交展开系数的权重且
[0017]本专利技术的有益效果为：本专利技术提供一种基于特征提取和面积度量的动态仿真结果验证方法。针对动态输出的仿真结果验证问题，首先采用正交多项式提取时间序列的均值、斜率、曲率等特征系数；其次，由于时间序列各个特征系数间相互独立，采用概率分布面积差异法度量各个特征系数的差异；最后，综合特征系数的差异得到动态输出的验证结果。此方法能解决具有不确定性的动态输出的仿真结果验证问题。同时，特征系数具有明确的物理含义，可根据正交多项式拟合效果和动态输出的实际物理含义选择待提取的特征系数。
附图说明
[0018]图1为单验证点动态输出仿真结果验证方法的流程图；
[0019]图2为多验证点动态输出仿真结果验证方法的流程图；
[0020]图3为视线角示意图，其中，图3(a)为视线角；图3(b)为2阶正交多项式拟合视线角；
[0021]图4为候选模型的面积度量示意图，其中，图4(a)为特征系数1的面积度量，(b) 为特征系数2的面积度量，(c)为特征系数3的面积度量；
[0022]图5为候选模型的U
‑
pooling度量示意图，其中，图5(a)为特征系数1的U
‑
pooling 度量；图5(b)为特征系数2的U
‑
pooling度量；图5(c)为特征系数3的U
‑
pooling度量。
具体实施方式
[0023]下面将结合本专利技术本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于特征提取和面积度量的动态仿真结果验证方法，假设动态输出的仿真时间序列表示为Y＝{y
i
}，i＝1,2,...,M，参考时间序列表示为R＝{r
i
}，i＝1,2,...,M，动态输出的仿真数据和参考数据包含多个时间序列，此动态仿真结果验证方法分为单验证点和多验证点两个部分，在单验证点，动态输出仿真结果验证包括以下步骤：S100、在单验证点x，得到仿真模型的动态输出为Y
l
，l＝1,2,...,N
s
，实际观测数据为R
k
，k＝1,2,...,N
r
；S200、采用离散切比雪夫多项式拟合时间序列，得到仿真数据和参考数据的正交展开系数分别为{α
sl
}
q
和{α
rk
}
q
，q＝1,...,Q，其中Q为正交展开系数的个数，由此得到{α
sl
}
q
的累积分布函数为F
sq
(α)，{α
rk
}
q
的经验累积分布函数为F
rq
(α)；S300、采用面积度量计算F
sq
(α)和F
rq
(α)的面积差异为d
q
，即并对面积差异归一化得到各个正交展开系数仿真数据和参考数据的一致性c
q
，面积差异的归一化计算公式为：其中x
max
和x
min<...

【专利技术属性】
技术研发人员：张欢，李伟，杨明，马萍，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：