The Ishigami function¶

The Ishigami function of Ishigami & Homma (1990) is recurrent test case for sensitivity analysis methods and uncertainty.

Let $a=7$ and $b=0.1$ (see Crestaux et al. (2007) and Marrel et al. (2009)). We consider the function

$g(X_1,X_2,X_3) = \sin(X_1)+a \sin (X_2)^2 + b X_3^4 \sin(X_1)$

for any $X_1,X_2,X_3\in[-\pi,\pi]$

We assume that the random variables $X_1,X_2,X_3$ are independent and have the uniform marginal distribution in the interval from $-\pi$ to $\pi$ :

$X_1,X_2,X_3\sim \mathcal{U}(-\pi,\pi).$

Analysis¶

The expectation and the variance of $Y$ are

$E(Y) = \frac{a}{2}$

and

$V(Y) = \frac{1}{2} + \frac{a^2}{8} + \frac{b^2 \pi^8}{18} + \frac{b\pi^4}{5}.$

The Sobol’ decomposition variances are

$V_1 = \frac{1}{2} \left(1 + b\frac{\pi^4}{5} \right)^2, \qquad V_2 = \frac{a^2}{8}, \qquad V_{1,3} = b^2 \pi^8 \frac{8}{225}$

and $V_3=V_{1,2} = V_{2,3}=V_{1,2,3} = 0$ .

This leads to the following first order Sobol’ indices:

$S_1 = \frac{V_1}{V(Y)}, \qquad S_2 = \frac{V_2}{V(Y)}, \qquad S_3 = 0,$

and the following total order indices:

$ST_1 = \frac{V_1+V_{1,3}}{V(Y)}, \qquad ST_2 = S_2, \qquad S_3 = \frac{V_{1,3}}{V(Y)}.$

The third variable $X_3$ has no effect at first order (because $X_3^4$ it is multiplied by $\sin(X_1)$ ), but has a total effet because of the interactions with $X_1$ . On the other hand, the second variable $X_2$ has no interactions which implies that the first order indice is equal to the total order indice for this input variable.

References¶

Ishigami, T., & Homma, T. (1990, December). An importance quantification technique in uncertainty analysis for computer models. In Uncertainty Modeling and Analysis, 1990. Proceedings., First International Symposium on (pp. 398-403). IEEE.

Sobol’, I. M., & Levitan, Y. L. (1999). On the use of variance reducing multipliers in Monte Carlo computations of a global sensitivity index. Computer Physics Communications, 117(1), 52-61.

Saltelli, A., Chan, K., & Scott, E. M. (Eds.). (2000). Sensitivity analysis (Vol. 134). New York: Wiley.

Crestaux, T., Martinez, J.-M., Le Maitre, O., & Lafitte, O. (2007). Polynomial chaos expansion for uncertainties quantification and sensitivity analysis. SAMO 2007, http://samo2007.chem.elte.hu/lectures/Crestaux.pdf.

Load the use case¶

We can load this classical model from the use cases module as follows :

>>> from openturns.usecases import ishigami_function as ishigami_function
>>> # Load the Ishigami use case
>>> im = ishigami_function.IshigamiModel()

API documentation¶

See IshigamiModel.

OpenTURNS

An Open source initiative for the Treatment of Uncertainties, Risks'N Statistics

Table of Contents

Previous topic

Next topic

This Page

The Ishigami function¶

Analysis¶

References¶

Load the use case¶

API documentation¶

Examples based on this use case¶