用于灵敏性分析的方法模型（Matlab代码实现）

           目录

1 概述

2 运行结果

3 参考文献

‍4 Matlab代码

---

1 概述

本代码用于灵敏性分析的方法模型。

2 运行结果

 

3 参考文献

[1]姚宏亮,苌健,王浩,李俊照.灵敏性分析下的因果网络参数的扰动学习研究[J].计算机科学与探索,2012,6(02):165-174.

‍4 Matlab代码

主函数部分代码：

clear all

n = 1000; % size of sample (n); number of simulations/runs

LHS_params = initialize_params(n);

R_0 = func_output(n, LHS_params);
R_0 = R_0';

% szy: need to get data from each params to use func_prcc()
numberParams = size(LHS_params, 2); % get number of params
LHS_matrix = [];
for i=1:numberParams
    LHS_matrix = [LHS_matrix LHS_params(i).data];
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   
% Sensitivity Analysis Results:

% PRCC (sensitivity indexes of R_0 with respect to each parameter):
PRCC_R0 = function_Prcc(R_0, LHS_matrix);

% plots:
plot_hist_paramsAndOutput(LHS_params, R_0);
plot_PRCC(PRCC_R0, numberParams);