两种基于方差的全局灵敏度分析W指标改进算法

doi:10.7527/S1000-6893.2016.0090

固体力学与飞行器总体设计

本期目录 | 过刊浏览 | 高级检索

前一篇 | 后一篇

两种基于方差的全局灵敏度分析W指标改进算法

巩祥瑞, 吕震宙, 左健巍

西北工业大学航空学院, 西安 710072

收稿日期:2015-04-10 修回日期:2016-03-22 出版日期:2016-06-15 发布日期:2016-03-25
通讯作者: 吕震宙女,博士,教授,博士生导师。主要研究方向:飞行器设计及可靠性工程。Tel:029-88460480 E-mail:zhenzhoulu@nwpu.edu.cn E-mail:zhenzhoulu@nwpu.edu.cn
作者简介:巩祥瑞男,硕士研究生。主要研究方向:飞行器设计及可靠性工程。Tel:029-88460480 E-mail:gxrui1991@163.com;左健巍男,硕士研究生。主要研究方向:飞行器设计及可靠性工程。Tel:029-88460480 E-mail:jianweizuo@mail.nwpu.edu.cn
基金资助:
国家自然科学基金(51475370);中央高校基本科研业务费专项资金(3102015BJ(II)CG009)

Two improved methods for variance-based global sensitivity analysis' W-indices

GONG Xiangrui, LYU Zhenzhou, ZUO Jianwei

School of Aeronautics, Northwestern Polytechnical University, Xi'an 710072, China

Received:2015-04-10 Revised:2016-03-22 Online:2016-06-15 Published:2016-03-25
Supported by:
National Natural Science Foundation of China (51475370);Fundamental Research Funds for the Central University (3102015BJ(II)CG009)

摘要/Abstract

摘要：

在全局灵敏度分析(SA)中,基于方差的灵敏度分析指标(包括Sobol指标和W指标)应用广泛。其中Sobol指标是将输入随机变量固定于特定点时,求得其对输出响应量的平均影响;而W指标求解当输入随机变量在各自分布区间上缩减变化时,输入变量对输出响应量的影响程度。相比Sobol指标,W指标所反映的信息更加全面。但目前对W指标的求解方法还比较欠缺,双层重复抽样蒙特卡罗(DLRS MC)方法和双层一次抽样蒙特卡罗(DLSS MC)方法是两种传统的求解方法。针对W指标的求解问题,提出了两种新算法:改进的蒙特卡罗模拟(AMCS)和基于稀疏网格积分(SGI)的方法。AMCS只需抽取一组样本便可计算出所有变量的各阶W指标,由于该方法是通过筛选策略来计算条件区间上的方差,避免了DLSS MC法中由于小数取整带来的计数误差,从而提高了计算W指标的精度。基于SGI的方法则利用稀疏网格积分来计算三重矩进而得到W指标,由于该方法继承了稀疏网格积分的高效性,因而进一步提高了W指标的计算效率。最后,给出了两个数值算例和一个工程算例,用于验证所提方法求解W指标的准确性和高效性。

关键词: 全局灵敏度分析, Sobol指标, W指标, 改进的蒙特卡罗模拟法, 稀疏网格积分法

Abstract:

In the global sensitivity analysis (SA), the variance-based sensitivity indices, such as Sobol's indices and W-indices, are used widely. Sobol's indices estimate the average variation of model output when input variables are fixed in their points. W-indices measure the average reduction of model output if input variables are reduced in their distributions. Compared with Sobol's indices, W-indices include more information. The double loop repeated set Monte Carlo (DLRS MC) and double loop single set Monte Carlo (DLSS MC) are two traditional methods, but these available methods for solving W-indices are still defective. In order to calculate W-indices efficiently, two new methods are presented. They are advanced Monte Carlo simulation (AMCS) and sparse grid integration (SGI)-based method. The AMCS only needs one set of samples to estimate all W-indices. Since screening method is used to estimate the variance in the conditional interval, and the count error induced by taking an integer in DLSS MC can be avoided, the accuracy of AMCS is higher than that of DLSS MC. The SGI-based method estimates W-indices by evaluating threefold statistic moment by SGI, in which the efficiency of the SGI is inherited. Finally, two numerical examples and an engineering example are employed to demonstrate the reasonability and efficiency of the presented methods.

Key words: global sensitivity analysis, Sobol's indices, W-indices, advanced Monte Carlo simulation method, sparse grid integration

中图分类号:

巩祥瑞, 吕震宙, 左健巍. 两种基于方差的全局灵敏度分析W指标改进算法[J]. 航空学报, 2016, 37(6): 1888-1898.

GONG Xiangrui, LYU Zhenzhou, ZUO Jianwei. Two improved methods for variance-based global sensitivity analysis' W-indices[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2016, 37(6): 1888-1898.

参考文献

[1] ZIO E, PODOFILLINI L. Accounting for components interactions in the differential importance measure[J]. Reliability Engineering & System Safety, 2006, 91(10-11):1163-1174.
[2] DO VAN P, BARROS A, BERENGUER C. Differential importance measure of Markov models using perturbation analysis[C]//European Safety and Reliability Conference (ESREL 2009). Abingdon, UK:Taylor & Francis, 2009:981-987.
[3] BORGONOVO E, APOSTOLAKIS G E. Comparison of global sensitivity analysis techniques and importance measures in PSA[J]. Reliability Engineering & System Safety, 2003, 79(2):175-185.
[4] VON GRIENSVEN A, MEIXNER T, GRUNWALD S, et al. A global sensitivity analysis tool for the parameters of multi-variable catchment models[J]. Journal of Hydrology, 2006, 324(1-4):10-23.
[5] SALTELLI A. Sensitivity analysis for importance assessment[J]. Risk Analysis, 2002, 22(3):579-590.
[6] HELTON J C, DAVIS F J. Latin hypercube sampling and the propagation of uncertainty in analyses of complex systems[J]. Reliability Engineering & System Safety, 2003, 81(1):23-69.
[7] HELTON J C. Uncertainty and sensitivity analysis techniques for use in performance assessment for radioactive waste disposal[J]. Reliability Engineering & System Safety, 1993, 42(2):327-367.
[8] SALTELLI A, MARIVOET J. Non-parametric statistics in sensitivity analysis for model output:A comparison of selected techniques[J]. Reliability Engineering & System Safety, 1990, 28(2):229-253.
[9] CHUN M H, HAN S J, TAK N I L. An uncertainty importance measure using a distance metric for the change in a cumulative distribution function[J]. Reliability Engineering & System Safety, 2000, 70(3):313-321.
[10] LIU H, CHEN W, SUDJIANTO A. Relative entropy based method for probabilistic sensitivity analysis in engineering design[J]. Journal of Mechanical Design, 2006, 128(2):326-336.
[11] BORGONOVO E. A new uncertainty importance measure[J]. Reliability Engineering & System Safety, 2007, 92(6):771-784.
[12] CUI L J, LU Z Z, ZHAO X P. Moment-independent importance measure of basic random variable and its probability density evolution solution[J]. Science China Technological Sciences, 2010, 53(4):1138-1145.
[13] LI L Y, LU Z Z, FENG J, et al. Moment-independent importance measure of basic variable and its state dependent parameter solution[J]. Structural Safety, 2012, 38:40-47.
[14] SOBOL' I M. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates[J]. Mathematics and Computers in Simulation, 2001, 55(1):271-280.
[15] SALTELLI A, ANNONI P, AZZINI I, et al. Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index[J]. Computer Physics Communications, 2010, 181(2):259-270.
[16] ARCHER G E B, SALTELLI A, SOBOL I M. Sensitivity measures, ANOVA-like techniques and the use of bootstrap[J]. Journal of Statistical Computation and Simulation, 1997, 58(2):99-120.
[17] SOBOL' I M. On sensitivity estimation for nonlinear mathematical models[J].Matematicheskoe Modelirovanie, 1990, 2(1):112-118.
[18] SALTELLI A, TARANTOLA S. On the relative importance of input factors in mathematical models:Safety assessment for nuclear waste disposal[J]. Journal of the American Statistical Association, 2002, 97(459):702-709.
[19] WEI P, LU Z Z, SONG J. A new variance-based global sensitivity analysis technique[J]. Computer Physics Communications, 2013, 184(11):2540-2551.
[20] TARANTOLA S, KOPUSTINSKAS V, BOLADO-LAVIN R, et al. Sensitivity analysis using contribution to sample variance plot:Application to a water hammer model[J]. Reliability Engineering & System Safety, 2012, 99:62-73.
[21] RATTO M, PAGANO A, YOUNG P. State dependent parameter metamodelling and sensitivity analysis[J]. Computer Physics Communications, 2007, 177(11):863-876.
[22] BARTHELMANN V, NOVAK E, RITTER K. High dimensional polynomial interpolation on sparse grids[J]. Advances in Computational Mathematics, 2000, 12(4):273-288.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

主管单位：中国科学技术协会主办单位：中国航空学会北京航空航天大学

两种基于方差的全局灵敏度分析W指标改进算法

Two improved methods for variance-based global sensitivity analysis' W-indices

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 3

编辑推荐

Metrics

本文评价

[1]	冯凯旋, 吕震宙, 蒋献. 基于偏导数的全局灵敏度指标的高效求解方法[J]. 航空学报, 2018, 39(3): 221699-221699.
[2]	刘付超, 魏鹏飞, 周长聪, 张政, 岳珠峰. 含旋转铰间隙平面运动机构可靠性灵敏度分析[J]. 航空学报, 2018, 39(11): 422133-422141.
[3]	阮文斌, 刘洋, 熊磊. 基于全局灵敏度分析的侧向气动导数不确定性对侧向飞行载荷的影响[J]. 航空学报, 2016, 37(6): 1827-1832.