sobol全局灵敏性分析.docxVIP

Sobol全局灵敏性分析 最近在研究全局敏感性分析方法中的Sobol方法，看了一些国内的论文，发现一个通病，就是公式一挂就可以得出结果了，真心觉得这种论文很“恶心”，主要原因是自己看不太懂。直到在维基百科上面找到了这种方法的详细解释，今天我们就根据网页上的步骤用一个例子来走一遍。?1.假设现在有一个函数:，从公式中可以看到有x1、x2、x3三个自变量对应变量Y有影响。 2.然后一般会给这三个参数一个取值范围，这里假设三个自变量的取值范围都设为[0,1]。 敏感性分析的目的就是求取这三个参数对于Y值得贡献。当然我们这边可能有人一下子就可以分析出那个参数对于Y值影响最大，但是在解决实际问题时，这个函数一般都是未知，我们只能将其视作一个黑盒子，只有输入和输出，这时我们对其进行敏感性分析就很有必要了。经过敏感性分析我们就能找出对结果影响较大的参数。这样对于调整结果是很有帮助的。?3.接着上面的例子，首先我们得根据三个自变量的范围进行采样，这边采样的方法一般都是蒙特卡洛采样以及一系列基于蒙特卡洛采样的变种，这个例子中我们采用了 Sobol sequence，具体的采样原理在这就不说了，大家可以自行谷歌。为了方便讲解例子我们设置采样的样本数为4(N = 4)，自变量数目为3(D = 3)。我们按照上述网页的步骤。?4.生成N * 2D（即4行6列）的样本矩阵。这个就是我们Sobol se