zl程序教程

您现在的位置是:首页 >  后端

当前栏目

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

MATLAB方法代码 实现 分析 模型 用于
2023-09-14 09:05:23 时间

 👨‍🎓个人主页:研学社的博客  

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。

⛳️座右铭:行百里者,半于九十。

📋📋📋本文目录如下:🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现


💥1 概述

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

📚2 运行结果

部分代码:

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);

部分代码:

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);

🎉3 参考文献

部分理论来源于网络,如有侵权请联系删除。

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

🌈4 Matlab代码实现