【MATLAB第109期】基于MATLAB的带置信区间的RSA区域敏感性分析方法,无目标函数

【MATLAB第108期】基于MATLAB的带置信区间的RSA区域敏感性分析方法,无目标函数

参考第64期文章【MATLAB第64期】【保姆级教程】基于MATLAB的SOBOL全局敏感性分析模型运用(含无目标函数,考虑代理模型)

创新点:

1、采用两种方法计算待置信区间的敏感性值

(1)通过对输出值设置阈值条件,绘制不同变量X与Y之间满足和不满足条件样本的累积分布函数(CDF),检查参数变化范围,并根据CDF之间的最大垂直距离绘制敏感度值。并使用自举法重复收敛分析以得出置信区间。

2、 将输入样本拆分为与输出等距范围相对应的"n组"数据集(Wagener等人于2001年首次提出),通过CDF之间最大垂直距离的统计量(最大值或中值)来评估不同数据集中X的CDF之间的距离,计算区域敏感度值。

1、建立代理模型

(1)数据设置:常用的案例数据 ,103*8 ,前7列代表输入变量, 最后1列代表因变量。

(2)数据变量上下限

X1-X7的最小值为:137 0 0 160 4.4 708 650

X1-X7的最大值为:374 193 260 240 19 1049.90 902

(3)代理模型选择:BP

2、数据抽样

matlab 复制代码
SampStrategy = 'lhs' ; % 拉丁超立方抽样
N = 400 ; % 样本数
Y = model_evaluation(myfun,X) ;

3、阈值设置

绘制输出样本(Y)与输入向量样本的一个分量(即X的一列)的关系,并对输出值设置阈值。

如阈值设置为20,如下图所示 :

X1-X7变量, 横坐标对应的是不同变量的数值范围,纵坐标代表Y的数值范围。

(1)绘制数据集X中满足给定阈值条件和不满足条件的CDF图


(2)通过平行坐标图检查行为参数化的范围

4a、绘制灵敏度指数(带阈值)

(1)参数CDF之间的最大垂直距离绘制

(2)并通过bootstrapping自举法评估鲁棒性,使用越来越多的样本重复计算,以评估收敛性,重复收敛分析以得出置信区间

4b、绘制灵敏度指数(分组法)

将数据集X中的样本拆分为与Y的'ngroup'等距值对应的'nggroups'子集。

然后,它通过CDF之间最大垂直距离的统计量(最大值或中值)来评估不同数据集中X的CDF之间的距离。

(1)将输入样本拆分为与输出等距范围相对应的"n组"数据集

(2)绘制参数CDF,自定义标签和图例:

5、代码获取

1.阅读首页置顶文章

2.关注CSDN

3.根据自动回复消息,私信回复"109期"以及相应指令,即可获取对应下载方式。

相关推荐
Evand J5 小时前
【MATLAB例程介绍】(17)永磁同步电机(PMSM)基于自适应无迹卡尔曼滤波(AUKF)的电流噪声在线估计。测量噪声协方差自适应与转速估计对比仿真
开发语言·matlab·卡尔曼滤波·无迹卡尔曼滤波·pmsm·电机滤波
严肃小白兔1 天前
MATLAB 绘图与科研论文插图制作:图形设计、统计表达与矢量导出
开发语言·matlab·科研绘图
严肃小白兔3 天前
MATLAB 数值计算中的精度与误差问题:浮点运算、误差传播与数值稳定性
人工智能·算法·matlab
Evand J3 天前
【PMSM例程】(14)永磁同步电机无迹卡尔曼滤波(UKF)定子电阻在线辨识与补偿。电阻变化跟踪及预测误差对比仿真。附完整代码的下载链接
matlab·电机·滤波·卡尔曼滤波·永磁同步电机·无迹卡尔曼滤波·pmsm
爱喝雪碧的可乐3 天前
MATLAB
开发语言·算法·matlab
严肃小白兔4 天前
MATLAB 优秀开源项目与工具箱推荐:涵盖通信、图像处理、脑机接口等九大领域
图像处理·matlab·信息与通信·开源项目
Evand J6 天前
【MATLAB例程12】永磁同步电机基于增广无迹卡尔曼滤波的未知负载转矩在线观测
开发语言·matlab
分享资料6 天前
MATLAB 中分号 (;)的作用
matlab
Evand J7 天前
【PMSM例程11】永磁同步电机(PMSM)基于UKF的转速与转子位置估计。附下载链接
算法·matlab·卡尔曼滤波·永磁同步电机·ukf·pmsm·位置估计
严肃小白兔8 天前
优秀 MATLAB 作者、书籍推荐:国外经典教材与技术资源整理
matlab