sobol灵敏度分析matlab_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 sequence
做的事情。这边我们生成
的矩阵为:
5
、
将矩阵的前
D
列设置为矩阵
A
,后
D
列设置为
B
列,在我们的例子中就是矩阵
m
的前
3
列设置为
矩阵
A
,后
3
列设置为矩阵
B
。
0.5
0.5
0.5
0.75
0.25
0.25
0.25
0.75
0.75
0.375
0.375
0.625
A
0.5
0.5
0.5
0.25
0.75
0.75
0.75
0.25
0.25
0.875
0.375
0.125
B
构造
N*D
的矩阵
AB
i
(i = 1,2,…,D),即用矩阵
B
中的第
i
列替换矩阵
A
的第
i
列,以本体为例:
2
4
1
2
3
1
Y = sin(
) + 7
(sin(
))
+ 0.1
sin(
)
x
x
x
x
0.5
0.5
0.5
0.5
0.5
0.5
0.75
0.25
0.25
0.25
0.75
0.75
0.25
0.75
0.75
0.75
0.25
0.25
0.375
0.375
0.625
0.875
0.375
0.12
5
m
sobol灵敏度分析matlab_sobol全局灵敏性分析相关推荐
- 数模笔记_多变量最优化的拉格朗日乘子方法中的灵敏性分析和影子价格
Date: 2_21 Name: Guo Yehao Theme: Sensitivity Analysis and Shadow Price in Optimality with multiple ...
- 数学建模灵敏度分析_数学建模中的灵敏度分析,到底在分析什么?
-2019/6/13- 概念 在数学建模中使用最优化方法时,我们常常会忽略对模型进行灵敏度分析,若缺少这一个步骤,会使得模型的可靠度受到质疑,那到底什么是灵敏度分析? 参考百度:灵敏度分析,是研究与分 ...
- 基于Autoware分析op_global_planner全局路径规划模块重规划
基于Autoware分析op_global_planner全局路径规划模块重规划 本次带来的是Autoware中的路网路径规划op_global_planner的重规划问题,在实际应用中我们肯定不仅仅 ...
- Linux内核 eBPF基础:kprobe原理源码分析:源码分析
Linux内核 eBPF基础 kprobe原理源码分析:源码分析 荣涛 2021年5月11日 在 <Linux内核 eBPF基础:kprobe原理源码分析:基本介绍与使用>中已经介绍了kp ...
- 对《GitHub服务中断24小时11分钟事故分析报告》的分析
对<GitHub服务中断24小时11分钟事故分析报告>的分析 声明 本文是根据公众号"高效开发运维"翻译整理的GitHub服务中断分析报告的信息进行分析.可能存在某些观 ...
- 仪器分析和计算机,化学分析和仪器分析的区别
化学分析和仪器分析的区别: 化学分析是指利用化学反应和它的计量关系来确定被测物质的组成和含量的一类分析方法.测定时需使用化学试剂.天平和一些玻璃器皿. 仪器分析(近代分析法或物理分析法):是基于与物质 ...
- R语言|clusterprofile超几何分布富集分析 GO,KEGG富集分析,循环Fisher‘s test
超几何分布富集分析 GO,KEGG富集分析,循环Fisher's test ID转换 GO KEGG 把KEGG里的geneid转回名称(readable table) 超几何分布以及生成data f ...
- eds能谱图分析实例_成分分析的四大神器—XRF、ICP、EDX和WDX
成分分析技术主要用于对未知物.未知成分等进行分析,通过成分分析技术可以快速确定目标样品中的各种组成成分是什么,帮助实验人员对样品进行定性定量分析,鉴别等.今天,小析姐就给大家介绍四种成分分析的常见设备 ...
- AI人工智能分析-人脸识别和分析(人脸检测跟踪、获取特征长度、提取用于人脸特征、比较相似度)
AI人工智能分析-人脸识别和分析(人脸检测跟踪.获取特征长度.提取用于人脸特征.比较相似度) 人工智能(Artificial Intelligence),英文缩写为AI.它是研究.开 ...
最新文章
- 网线的制作方法和千兆网线的施工注意事项
- linux ipv6 组播地址,ipv6组播地址是什么意思
- Airtable(数据管理)
- 借助JVM生日的时机,说说关于JVM你所不知道的那些事
- android ztree,ztree实现权限功能(横向显示)
- DDD战略设计相关核心概念的理解
- A_Pancers团队项目设计完善编码测试
- django会话跟踪技术
- 谷歌修复安卓System 组件中的多个 RCE 漏洞
- gcc动态链接库基本知识 [转]
- 2020年20种最佳Android应用程序模板
- sir节点matlab模型,SIR模型实现(matlab)
- 傻瓜也能看懂的C语言病毒代码(无害)
- mediasoup server 启动失败排查
- Erphp loggedin 异地IP登录自动禁封用户 WordPress插件
- Javascript中Promise对象
- ACK Acknowledgement 确认 AES Advanced Encryption Standard 高级加密标准 ATM Asynchronous Transfer Mode异步传输模式
- HLK-W801-LVGL8之触摸驱动
- for循环与range方法
- 魂断“中国百慕大”?GIS斩祸根!