python如何做敏感度分析,使用SALib工具箱从测量数据进行Python敏感性分析

I would like to understand, how to use the SALib python toolbox to perform a Sobol sensitivity analysis (to study parameters and crossed parameters influence)

From the original example I'm supposed to proceed this way:

from SALib.sample import saltelli

from SALib.analyze import sobol

from SALib.test_functions import Ishigami

import numpy as np

problem = {

'num_vars': 3,

'names': ['x1', 'x2', 'x3'],

'bounds': [[-np.pi, np.pi]]*3

}

# Generate samples

param_values = saltelli.sample(problem, 1000)

# Run model (example)

Y = Ishigami.evaluate(param_values)

# Perform analysis

Si = sobol.analyze(problem, Y, print_to_console=True)

# Returns a dictionary with keys 'S1', 'S1_c