CN106202694A - 基于组合预测方法的组合Kriging模型构建方法 - Google Patents
基于组合预测方法的组合Kriging模型构建方法 Download PDFInfo
- Publication number
- CN106202694A CN106202694A CN201610530706.8A CN201610530706A CN106202694A CN 106202694 A CN106202694 A CN 106202694A CN 201610530706 A CN201610530706 A CN 201610530706A CN 106202694 A CN106202694 A CN 106202694A
- Authority
- CN
- China
- Prior art keywords
- function
- kriging
- combination
- model
- kriging model
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于组合预测方法的组合Kriging模型构建方法,先使用7种不同类型相关函数分别构建7个单一相关函数Kriging近似模型,然后以平均绝对相对误差最小为优化目标,构建每个单一相关函数Kriging模型的权重值的优化模型,获得最优权重值,最后利用最优权重值和基于7种不同类型相关函数构建的7个单一相关函数Kriging近似模型,构建组合Kriging模型,本发明提出的组合Kriging模型具有预测精度高、预测性能更加稳定的特点,在基于计算机数值仿真计算的工程优化设计领域具有较高的实用价值。
Description
技术领域
本发明属于工程设计优化技术领域,尤其涉及一种基于组合预测方法的组合Kriging模型构建方法。
技术背景
近年来,高精度的数值仿真技术在工程设计及优化工作中得到了广泛的应用,采用数值仿真技术可以精确的模拟工程结构在不同受力状态和工况环境下的力学特性与工作性能,为工程设计人员提供设计依据,在提高工程设计效率的同时降低设计成本。为了获得更加精确的仿真计算结果,通常需要构建更精细的仿真模型,从而使得数值仿真模型计算规模增大,导致数值仿真计算耗时增加,降低了工程设计效率。为了缩小计算规模,提高计算效率,近似模型(Surrogate model)技术在高精度数值仿真计算中开始得到应用。近似模型能够提供高精度的预测功能,减少工程优化设计中的仿真计算时间,有效改善设计效率。
Kriging模型是一种典型的近似模型,具有计算效率高、非线性拟合能力强的特点,能够适应高维非线性数据的拟合与预测。因此,其在工程优化设计领域得到了广泛的应用。相关函数(Spatial correlation function)是Kriging模型获得样本信息的重要工具,在构建近似模型过程中可以基于不同的类型的相关函数构建不同类型的Kriging模型。典型的Kriging相关函数包括三次函数(Cubic)、指数函数(Exp)、幂函数(Expg)、高斯函数(Gaussian)、线性函数(Linear)、球函数(Spherical)和样条函数(Spline)7种。对同一问题进行预测时,基于不同相关函数构建的Kriging模型的预测精度是不同的。从组合预测的思想来看,使用单一的、不合适的相关函数构建Kriging模型,会使得模型获得的样本信息具有片面性,导致部分样本信息丢失、样本信息不完整的问题,不利于保证Kriging模型预测精度。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种基于组合预测方法的组合Kriging模型构建方法,构建的Kriging模型具有预测精度更高、稳定性更好的特点。
为了达到上述目的,本发明采取的技术方案为:
一种基于组合预测方法的组合Kriging模型构建方法,包括以下步骤:
1)设已知样本点设计变量X={x1,x2,…,xn},其对应的响应量为Y={y(x1),y(x2),…,y(xn)};
2)构建单一相关函数Kriging模型:以样本点数据(X,Y)为对象,分别基于7个不同相关函数构建7个单一Kriging模型7个不同相关函数为三次函数(Cubic)、指数函数(Exp)、幂函数(Expg)、高斯函数(Gaussian)、线性函数(Linear)、球函数(Spherical)和样条函数(Spline);
3)确定各单一相关函数Kriging模型的最优权重值ωi:对于单一相关函数Kriging模型要求其对应的权重值具有如下特性:a)预测精度最高的单一相关函数Kriging模型权重值最大;b)所有权重值之和
以平均绝对相对误差(Mean absolute percentage error,MAPE)最小为优化目标,构建权重值ωi的优化模型,对7个权重值进行优化,获得权重值ω1,ω2,…,ω7;
4)构造组合Kriging模型:在获得7个最优权重值之后,基于组合预测方法的Kriging模型表述为
本发明的有益效果:本发明引入组合预测方法对Kriging模型样本信息的完整性进行改善,综合考虑不同类型相关函数描述获得的样本信息,将单一相关函数的Kriging模型视作为组合模型的子模型,以减少由于单一、不合适相关函数选择而带来的预测误差。与现有的单一相关函数Kriging模型相比,本发明提出的组合Kriging模型的样本信息由多个相关函数描述,可以较好的避免样本信息丢失的问题,能够提高样本信息的完整性,其预测精度更高、预测性能更加稳定,在工程优化设计领域中具有较高的实际应用价值。本发明减少工程优化问题的计算成本,提高工程优化设计效率。
附图说明
图1为本发明方法的流程图。
图2是实施例二维函数原始响应面。
图3是实施例基于7个相关函数构建的单一相关函数Kriging模型预测结果。
图4是实施例组合Kriging模型预测结果。
具体实施方式
以下结合附图及实施例对本发明作进一步说明。
参照图1,一种基于组合预测方法的组合Kriging模型构建方法,包括以下步骤:
1)设已知样本点设计变量X={x1,x2,…,xn},其对应的响应量为Y={y(x1),y(x2),…,y(xn)};
考虑一个如式(1)所示的二维函数,其对应的一组样本点对应的响应值为其原始响应面如图2所示,
2)构建单一相关函数Kriging模型:不同相关函数获得的样本点数据信息是不同的,为了获得所有7个相关函数获得的样本点数据信息,以样本点数据(X,Y)为对象,分别基于7个不同相关函数构建7个单一Kriging模型 参照图3,7个不同相关函数为三次函数(Cubic)、指数函数(Exp)、幂函数(Expg)、高斯函数(Gaussian)、线性函数(Linear)、球函数(Spherical)和样条函数(Spline);
3)确定各单一相关函数Kriging模型的最优权重值ωi:为了充分利用不同相关函数获得的样本点数据信息,需要合理的确定各单一相关函数Kriging模型的权重值;对于单一相关函数Kriging模型要求其对应的权重值具有如下特性:a)预测精度最高的单一相关函数Kriging模型权重值最大;b)所有权重值之和
从样本数据集合中抽取m个样本点,以平均绝对相对误差(Mean absolutepercentage error,MAPE)最小为优化目标,构建如式(2)所示的权重值ωi优化模型,获得优化的各单一Kriging模型权重值;
Find:ωi
ωi≥0
其中,y(xi)是第i个样本点的实际响应值;是利用组合Kriging模型获得的第i个样本点的预测值;
4)构造组合Kriging模型:在获得7个最优权重值之后,构建如式(3)所示的该二维函数组合Kriging模型,
参照图4,图4给出了组合Kriging模型的预测结果,表1为该组合Kriging模型预测精度与其他单一相关函数Kriging模型预测精度对比。
表1组合Kriging模型与其他单一相关函数Kriging模型预测精度对比
从表1可知,组合Kriging模型比其他7种单一相关函数Kriging模型的预测精度及稳定性更高。
Claims (1)
1.一种基于组合预测方法的组合Kriging模型构建方法,其特征在于,包括以下步骤:
1)设已知样本点设计变量X={x1,x2,…,xn},其对应的响应量为Y={y(x1),y(x2),…,y(xn)};
2)构建单一相关函数Kriging模型:以样本点数据(X,Y)为对象,分别基于7个不同相关函数构建7个单一Kriging模型7个不同相关函数为三次函数(Cubic)、指数函数(Exp)、幂函数(Expg)、高斯函数(Gaussian)、线性函数(Linear)、球函数(Spherical)和样条函数(Spline);
3)确定各单一相关函数Kriging模型的最优权重值ωi:对于单一相关函数Kriging模型要求其对应的权重值具有如下特性:a)预测精度最高的单一相关函数Kriging模型权重值最大;b)所有权重值之和
以平均绝对相对误差(Mean absolute percentage error,MAPE)最小为优化目标,构建权重值ωi的优化模型,对7个权重值进行优化,获得权重值ω1,ω2,…,ω7;
4)构造组合Kriging模型:在获得7个最优权重值之后,基于组合预测方法的Kriging模型表述为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610530706.8A CN106202694A (zh) | 2016-07-06 | 2016-07-06 | 基于组合预测方法的组合Kriging模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610530706.8A CN106202694A (zh) | 2016-07-06 | 2016-07-06 | 基于组合预测方法的组合Kriging模型构建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106202694A true CN106202694A (zh) | 2016-12-07 |
Family
ID=57473594
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610530706.8A Pending CN106202694A (zh) | 2016-07-06 | 2016-07-06 | 基于组合预测方法的组合Kriging模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106202694A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107704949A (zh) * | 2017-09-12 | 2018-02-16 | 中国船舶重工集团公司第七0九研究所 | 基于历史数据驱动的航班进离场调度优化方法及系统 |
CN109752632A (zh) * | 2019-01-26 | 2019-05-14 | 云南电网有限责任公司迪庆供电局 | 一种高海拔输电线路金具起晕电压海拔校正方法 |
CN111735686A (zh) * | 2020-05-16 | 2020-10-02 | 厦门大学 | 一种用于钢骨混凝土结构性能分析的基于备选模型的Kriging预测方法 |
-
2016
- 2016-07-06 CN CN201610530706.8A patent/CN106202694A/zh active Pending
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107704949A (zh) * | 2017-09-12 | 2018-02-16 | 中国船舶重工集团公司第七0九研究所 | 基于历史数据驱动的航班进离场调度优化方法及系统 |
CN107704949B (zh) * | 2017-09-12 | 2021-09-03 | 中国船舶重工集团公司第七0九研究所 | 基于历史数据驱动的航班进离场调度优化方法及系统 |
CN109752632A (zh) * | 2019-01-26 | 2019-05-14 | 云南电网有限责任公司迪庆供电局 | 一种高海拔输电线路金具起晕电压海拔校正方法 |
CN109752632B (zh) * | 2019-01-26 | 2021-08-10 | 云南电网有限责任公司迪庆供电局 | 一种高海拔输电线路金具起晕电压海拔校正方法 |
CN111735686A (zh) * | 2020-05-16 | 2020-10-02 | 厦门大学 | 一种用于钢骨混凝土结构性能分析的基于备选模型的Kriging预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104636563B (zh) | 高速压力机上横梁可靠性设计方法 | |
CN105843073B (zh) | 一种基于气动力不确定降阶的机翼结构气动弹性稳定性分析方法 | |
CN110619432B (zh) | 一种基于深度学习的特征提取水文预报的方法 | |
Ma et al. | Ultra-short-term wind generation forecast based on multivariate empirical dynamic modeling | |
CN105203106B (zh) | 一种基于模拟退火粒子群算法的wMPS网络布局优化方法 | |
CN112819962B (zh) | 数字图像相关中非均匀网格划分及局部网格疏密方法 | |
CN105930562A (zh) | 一种非概率条件下的结构性能优化设计方法 | |
CN105260607A (zh) | 一种串并联耦合的多模型水文预报方法 | |
CN105868853A (zh) | 一种短期风电功率组合概率预测方法 | |
CN106202694A (zh) | 基于组合预测方法的组合Kriging模型构建方法 | |
CN103985000A (zh) | 基于函数型非参数回归的中长期典型日负荷曲线预测方法 | |
CN103106331B (zh) | 基于降维和增量式极限学习机的光刻线宽智能预测方法 | |
CN112329349B (zh) | 一种边坡可靠度评估方法 | |
CN110069854A (zh) | 多重分辨tweno格式对可压流场问题的模拟方法 | |
CN102419827A (zh) | 基于rbf神经网络的沸腾换热预测方法 | |
CN109871605A (zh) | 折弯成形方法、系统、装置及介质 | |
CN102254184B (zh) | 一种多物理域特征信息融合方法 | |
CN113849943A (zh) | 一种耦合压力先验信息的供水管网节点需水量校核方法 | |
CN116595827B (zh) | 无限维度条带喷丸成形工艺规划方法和系统 | |
CN104008258A (zh) | 基于位移形态的钢结构火灾温度场逆向推定方法 | |
CN113722961B (zh) | 一种基于广义协同高斯过程模型的结构不确定性量化解析方法 | |
CN106709167A (zh) | 一种基于相似度理论的装配有效性评估方法 | |
Storcz et al. | Applicability of Multivariate Linear Regression in Building Energy Demand Estimation. | |
CN114297582A (zh) | 基于多探针局部敏感哈希负二项式回归模型的离散计数数据的建模方法 | |
CN110673088B (zh) | 混合视距和非视距环境中基于到达时间的目标定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161207 |