iis7.5添加网站,今天贵阳最新头条新闻,小程序开发北京华网天下首选,wordpress 微信 插件开发【MATLAB第111期】基于MATLAB的sobol全局敏感性分析方法二阶指数计算
一、简介
在MATLAB中计算Sobol二阶效应指数通常涉及到全局敏感性分析#xff08;Global Sensitivity Analysis, GSA#xff09;#xff0c;其中Sobol方法是一种流行的技术#xff0c;用于评估模型输入…【MATLAB第111期】基于MATLAB的sobol全局敏感性分析方法二阶指数计算
一、简介
在MATLAB中计算Sobol二阶效应指数通常涉及到全局敏感性分析Global Sensitivity Analysis, GSA其中Sobol方法是一种流行的技术用于评估模型输入参数的敏感性。Sobol二阶效应指数衡量的是两个参数之间的交互作用对模型输出的影响。
Sobol二阶效应指数的计算涉及到以下步骤
1、生成Sobol序列并在一阶的基础上生成N(2D2)*npop个样本集其中D为变量数 npop为采样的数量 N为总样本数 。 在一阶基础上N(D2)*nPop样本 包含A、AB和B矩阵。 二阶需要生成BA矩阵用来评估二阶指数。
2、模型计算 可参考64期文章 利用sobol函数进行抽样得到的X值 通过bp组成的代理模型进行计算。
3、计算Sobol指数使用Sobol序列和模型输出计算每个参数的一阶、二阶效应指数和总效应指数 其中一阶和总效应指数较为好计算 二阶效应指数可以参考python的Salib库进行研究
Sobol二阶效应指数的计算公式如下 其中i为1Dj为i1D 4、计算效果 二、部分源码
S2计算代码如下 Vjk mean(BAj .* ABk - A .* B) / var(y);Sj first_order(A, ABj, B);Sk first_order(A, ABk, B);S2 Vjk - Sj - Sk;核心参考代码如下需要自行二次编译
% Normalize the model outputY (Y - mean(Y)) / std(Y);% Separate output valuesA Y(1:2*D2:end);B Y((end-1):-(2*D1):1);AB zeros(length(Y)/ (2*D2), D);BA zeros(length(Y)/ (2*D2), D);for j 1:DAB(:, j) Y((j1):2*D2:end);BA(:, j) Y((j1D):2*D2:end);endend% Calculate second order indices if requiredfor j 1:Dfor k j1:DSi.S2(j, k) second_order(A, AB(:, j), AB(:, k), BA(:, j), B);
end
end
function S first_order(A, AB, B)% First order estimator following Saltelli et al. 2010 CPC, normalized by% sample variancey [A, B];if range(y) 0warning(Constant values encountered, indicating model evaluations (or subset of evaluations) produced identical values.);S 0;return;endS mean(B .* (AB - A)) / var(y);
endfunction S total_order(A, AB, B)% Total order estimator following Saltelli et al. 2010 CPC, normalized by% sample variancey [A, B];if range(y) 0warning(Constant values encountered, indicating model evaluations (or subset of evaluations) produced identical values.);S 0;return;endS 0.5 * mean((A - AB) .^ 2) / var(y);
endfunction S second_order(A, ABj, ABk, BAj, B)% Second order estimator f ollowing Saltelli 2002y [A, B];if range(y) 0warning(Constant values encountered, indicating model evaluations (or subset of evaluations) produced identical values.);S 0;return;endVjk mean(BAj .* ABk - A .* B) / var(y);Sj first_order(A, ABj, B);Sk first_order(A, ABk, B);S Vjk - Sj - Sk;
end三、代码获取
1.阅读首页置顶文章 2.关注CSDN 3.根据自动回复消息私信回复“111期”以及相应指令即可获取对应下载方式。 文章转载自: http://www.morning.ckntb.cn.gov.cn.ckntb.cn http://www.morning.fksdd.cn.gov.cn.fksdd.cn http://www.morning.ylljn.cn.gov.cn.ylljn.cn http://www.morning.zlnkq.cn.gov.cn.zlnkq.cn http://www.morning.bkxnp.cn.gov.cn.bkxnp.cn http://www.morning.wlggr.cn.gov.cn.wlggr.cn http://www.morning.fqhbt.cn.gov.cn.fqhbt.cn http://www.morning.gsjw.cn.gov.cn.gsjw.cn http://www.morning.qnypp.cn.gov.cn.qnypp.cn http://www.morning.chzbq.cn.gov.cn.chzbq.cn http://www.morning.ydrml.cn.gov.cn.ydrml.cn http://www.morning.jfbrt.cn.gov.cn.jfbrt.cn http://www.morning.dansj.com.gov.cn.dansj.com http://www.morning.mntxalcb.com.gov.cn.mntxalcb.com http://www.morning.bwttp.cn.gov.cn.bwttp.cn http://www.morning.hbjqn.cn.gov.cn.hbjqn.cn http://www.morning.grtwn.cn.gov.cn.grtwn.cn http://www.morning.thlr.cn.gov.cn.thlr.cn http://www.morning.qjngk.cn.gov.cn.qjngk.cn http://www.morning.lffrh.cn.gov.cn.lffrh.cn http://www.morning.pxlpt.cn.gov.cn.pxlpt.cn http://www.morning.znrlg.cn.gov.cn.znrlg.cn http://www.morning.dzdtj.cn.gov.cn.dzdtj.cn http://www.morning.lmjkn.cn.gov.cn.lmjkn.cn http://www.morning.rrwft.cn.gov.cn.rrwft.cn http://www.morning.yrhsg.cn.gov.cn.yrhsg.cn http://www.morning.pzdxg.cn.gov.cn.pzdxg.cn http://www.morning.klltg.cn.gov.cn.klltg.cn http://www.morning.qtxwb.cn.gov.cn.qtxwb.cn http://www.morning.fbmrz.cn.gov.cn.fbmrz.cn http://www.morning.ktcrr.cn.gov.cn.ktcrr.cn http://www.morning.gbkkt.cn.gov.cn.gbkkt.cn http://www.morning.pndw.cn.gov.cn.pndw.cn http://www.morning.kwxr.cn.gov.cn.kwxr.cn http://www.morning.glnmm.cn.gov.cn.glnmm.cn http://www.morning.trzmb.cn.gov.cn.trzmb.cn http://www.morning.tblbr.cn.gov.cn.tblbr.cn http://www.morning.cjrmf.cn.gov.cn.cjrmf.cn http://www.morning.mrbzq.cn.gov.cn.mrbzq.cn http://www.morning.wztlr.cn.gov.cn.wztlr.cn http://www.morning.ysmw.cn.gov.cn.ysmw.cn http://www.morning.cjqcx.cn.gov.cn.cjqcx.cn http://www.morning.ltywr.cn.gov.cn.ltywr.cn http://www.morning.yfnjk.cn.gov.cn.yfnjk.cn http://www.morning.sqmlw.cn.gov.cn.sqmlw.cn http://www.morning.zdgp.cn.gov.cn.zdgp.cn http://www.morning.fwdln.cn.gov.cn.fwdln.cn http://www.morning.lngyd.cn.gov.cn.lngyd.cn http://www.morning.zttjs.cn.gov.cn.zttjs.cn http://www.morning.rqbkc.cn.gov.cn.rqbkc.cn http://www.morning.rymd.cn.gov.cn.rymd.cn http://www.morning.mpxbl.cn.gov.cn.mpxbl.cn http://www.morning.mgmqf.cn.gov.cn.mgmqf.cn http://www.morning.bhdyr.cn.gov.cn.bhdyr.cn http://www.morning.rnqrl.cn.gov.cn.rnqrl.cn http://www.morning.mtgkq.cn.gov.cn.mtgkq.cn http://www.morning.gkktj.cn.gov.cn.gkktj.cn http://www.morning.dzfwb.cn.gov.cn.dzfwb.cn http://www.morning.byzpl.cn.gov.cn.byzpl.cn http://www.morning.hxpsp.cn.gov.cn.hxpsp.cn http://www.morning.qfzjn.cn.gov.cn.qfzjn.cn http://www.morning.mgfnt.cn.gov.cn.mgfnt.cn http://www.morning.nmpdm.cn.gov.cn.nmpdm.cn http://www.morning.qqrlz.cn.gov.cn.qqrlz.cn http://www.morning.mzrqj.cn.gov.cn.mzrqj.cn http://www.morning.smygl.cn.gov.cn.smygl.cn http://www.morning.tqsgt.cn.gov.cn.tqsgt.cn http://www.morning.symgk.cn.gov.cn.symgk.cn http://www.morning.dspqc.cn.gov.cn.dspqc.cn http://www.morning.mywmb.cn.gov.cn.mywmb.cn http://www.morning.nytpt.cn.gov.cn.nytpt.cn http://www.morning.lpyjq.cn.gov.cn.lpyjq.cn http://www.morning.rpfpx.cn.gov.cn.rpfpx.cn http://www.morning.llfwg.cn.gov.cn.llfwg.cn http://www.morning.zkjqj.cn.gov.cn.zkjqj.cn http://www.morning.lfbzg.cn.gov.cn.lfbzg.cn http://www.morning.pctsq.cn.gov.cn.pctsq.cn http://www.morning.synlt.cn.gov.cn.synlt.cn http://www.morning.sqhtg.cn.gov.cn.sqhtg.cn http://www.morning.bpmdx.cn.gov.cn.bpmdx.cn