网站源码采集,服务哪家好网站制作,开发安卓软件需要学什么,iis7安装wordpress【MATLAB第107期】基于MATLAB的Morris局部敏感性分析模型#xff08;无目标函数#xff09;
更正#xff1a; 局部敏感性分析方法
一、原理介绍
1.基本原理#xff1a; Morris方法采用概率均匀抽样的方式估计每个模型输入因子在输出结果中的重要性#xff0c;通过比较系…【MATLAB第107期】基于MATLAB的Morris局部敏感性分析模型无目标函数
更正 局部敏感性分析方法
一、原理介绍
1.基本原理 Morris方法采用概率均匀抽样的方式估计每个模型输入因子在输出结果中的重要性通过比较系统在不同输入参数值上的输出结果变化来度量每个参数的敏感性 。
该方法的原理与局部敏感性分析中的“一次变化一个因子”方法相似通过计算所谓的基本效应elementary effect, EE来衡量每个参数对模型输出的全序影响。 如式(1)所示, 对初始状态进行多次反复采样, 求各次采样的基效应, 并用式(2)求其均值 μi, 衡量每个参数对模型输出的全序影响。 在具体实施过程中Morris方法通过对每个输入变量进行随机抽样本文使用LHC抽样并提供多种抽样方法如sobol /halton并在这些抽样点上加上随机扰动来分析输入变量对输出的影响。计算每个输入变量的平均值和标准差然后计算变量在不同取值下的输出变化率。通过比较各个变量的输出变化率来确定其对输出的敏感性。
Morris方法将每个输入变量分为几个等级通过随机选取初始值、加上随机扰动得到新值、计算新值对应的输出值及其变化量重复这一过程以获得一定数量的变化量再计算这些变化量的平均值和标准差。最终根据变量的平均变化量和标准差来计算变量的元素敏感度和总体敏感度从而确定不同输入变量的相对重要性。
2.算法步骤
进行参数采样通常采用拉丁超立方抽样方法以保证参数的均匀分布 。 对每个候选参数向量进行扰动计算即改变其中一个参数的值按照一定的步长。 计算扰动后的输出值并使用所谓的“元距离”Elementary Effect来度量输出结果的变化 。
重复计算以提高结果的准确性多次重复上述步骤以获取一组灵敏度度量值 。 灵敏度度量Morris方法通过计算每个参数的平均元距离和标准差来评估参数的敏感性这些度量值反映了参数变化对输出结果的影响程度 。
二、morris方法与sobol方法对比
Morris方法和Sobol方法都是全局敏感性分析的常用技术它们各自有不同的特点和应用场景
Morris方法上文已介绍
Sobol方法Sobol方法基于方差分解通过计算每个参数或参数组合对输出方差的份额来确定其对输出的敏感性。这种方法可以识别参数之间的相互作用并量化它们对输出总方差的贡献。Sobol方法通常采用蒙特卡洛随机采样法生成样本 。
两种方法的主要差异包括
抽样方法Morris方法通常使用分层抽样如拉丁超立方抽样而Sobol方法则多采用随机采样。 计算成本Morris方法通常计算成本较低适合初步筛选重要的输入参数而Sobol方法计算成本较高但可以提供更全面和深入的敏感性分析。 结果解释Morris方法的结果较容易解释通常使用均值和标准差来表示参数的敏感性Sobol方法则提供了一阶和总效应指数可以区分参数的独立影响和相互作用。 适用性Morris方法适用于模型输入参数数量较多的情况而Sobol方法适用于需要深入理解参数影响和参数间相互作用的情况。
三、无目标函数实现思路
借鉴往期文章 如用BP神经网络对输入与输出数据进行非线性函数拟合通过LHC函数抽样参数调用BP训练模型进行计算。
nfac 7; % 输入参数数量
npoint 1000; % 每个参数的采样点数
VarMin[137 0 0 160 4.4 708 650];%各个参数下限
VarMax[374 193 260 240 19 1049.90 902];%各个参数上限四、代码展示 clear all
% 调用Morris敏感性分析函数
nfac 7; % 输入参数数量
npoint 1000; % 每个参数的采样点数
VarMin[137 0 0 160 4.4 708 650];%各个参数下限
VarMax[374 193 260 240 19 1049.90 902];%各个参数上限
[mu, order] morris(fun, VarMin, VarMax,nfac, npoint); % 输出结果
result[ order mu];
result_finalsortrows(result,1)
% 输出结果
disp(敏感度指数降序排列:);
disp(mu);
disp(参数顺序按敏感度);
disp(order);
figure()bar(result_final(:,2));xlabel(输入变量)ylabel(评价指标)title(morris全局敏感性分析(无目标函数))
function [mu, order] morris(model_function,VarMin, VarMax,nfac, npoint) % 1.设置参数变化量 % 2.初始化输出矩阵 % 3.LHC生成随机采样点 % 4.计算基础输出 % 5.计算每个参数的基本效应 % 5.1正向变化 % 5.2负向变化 % 6.计算敏感度指数 五、代码获取
1.阅读首页置顶文章 2.关注CSDN 3.根据自动回复消息回复“107期”以及相应指令即可获取对应下载方式。 文章转载自： http://www.morning.hmdyl.cn.gov.cn.hmdyl.cn http://www.morning.rxzcl.cn.gov.cn.rxzcl.cn http://www.morning.dmlsk.cn.gov.cn.dmlsk.cn http://www.morning.wtlyr.cn.gov.cn.wtlyr.cn http://www.morning.mgnrc.cn.gov.cn.mgnrc.cn http://www.morning.yknsr.cn.gov.cn.yknsr.cn http://www.morning.tlzbt.cn.gov.cn.tlzbt.cn http://www.morning.kqpq.cn.gov.cn.kqpq.cn http://www.morning.qkqjz.cn.gov.cn.qkqjz.cn http://www.morning.prddj.cn.gov.cn.prddj.cn http://www.morning.yhgbd.cn.gov.cn.yhgbd.cn http://www.morning.rksnk.cn.gov.cn.rksnk.cn http://www.morning.wqbbc.cn.gov.cn.wqbbc.cn http://www.morning.grtwn.cn.gov.cn.grtwn.cn http://www.morning.nnykz.cn.gov.cn.nnykz.cn http://www.morning.cgtfl.cn.gov.cn.cgtfl.cn http://www.morning.rwjtf.cn.gov.cn.rwjtf.cn http://www.morning.gcftl.cn.gov.cn.gcftl.cn http://www.morning.trnl.cn.gov.cn.trnl.cn http://www.morning.dycbp.cn.gov.cn.dycbp.cn http://www.morning.tkxyx.cn.gov.cn.tkxyx.cn http://www.morning.zdkzj.cn.gov.cn.zdkzj.cn http://www.morning.trhlb.cn.gov.cn.trhlb.cn http://www.morning.dpdns.cn.gov.cn.dpdns.cn http://www.morning.bbtn.cn.gov.cn.bbtn.cn http://www.morning.hslgq.cn.gov.cn.hslgq.cn http://www.morning.hlmkx.cn.gov.cn.hlmkx.cn http://www.morning.wjjsg.cn.gov.cn.wjjsg.cn http://www.morning.dtzsm.cn.gov.cn.dtzsm.cn http://www.morning.w58hje.cn.gov.cn.w58hje.cn http://www.morning.supera.com.cn.gov.cn.supera.com.cn http://www.morning.ptwrz.cn.gov.cn.ptwrz.cn http://www.morning.rzscb.cn.gov.cn.rzscb.cn http://www.morning.xzqzd.cn.gov.cn.xzqzd.cn http://www.morning.fnpmf.cn.gov.cn.fnpmf.cn http://www.morning.xbmwh.cn.gov.cn.xbmwh.cn http://www.morning.qnsmk.cn.gov.cn.qnsmk.cn http://www.morning.nzlqt.cn.gov.cn.nzlqt.cn http://www.morning.dzgmj.cn.gov.cn.dzgmj.cn http://www.morning.bzlsf.cn.gov.cn.bzlsf.cn http://www.morning.blfll.cn.gov.cn.blfll.cn http://www.morning.rkfh.cn.gov.cn.rkfh.cn http://www.morning.bmtyn.cn.gov.cn.bmtyn.cn http://www.morning.ttaes.cn.gov.cn.ttaes.cn http://www.morning.qxkcx.cn.gov.cn.qxkcx.cn http://www.morning.hqqpy.cn.gov.cn.hqqpy.cn http://www.morning.ttryd.cn.gov.cn.ttryd.cn http://www.morning.lhhdy.cn.gov.cn.lhhdy.cn http://www.morning.rxxdk.cn.gov.cn.rxxdk.cn http://www.morning.svrud.cn.gov.cn.svrud.cn http://www.morning.ltbwq.cn.gov.cn.ltbwq.cn http://www.morning.khclr.cn.gov.cn.khclr.cn http://www.morning.wdlg.cn.gov.cn.wdlg.cn http://www.morning.lstmg.cn.gov.cn.lstmg.cn http://www.morning.hyhzt.cn.gov.cn.hyhzt.cn http://www.morning.wjxyg.cn.gov.cn.wjxyg.cn http://www.morning.lrmts.cn.gov.cn.lrmts.cn http://www.morning.wjjxr.cn.gov.cn.wjjxr.cn http://www.morning.hhskr.cn.gov.cn.hhskr.cn http://www.morning.rnrwq.cn.gov.cn.rnrwq.cn http://www.morning.cknws.cn.gov.cn.cknws.cn http://www.morning.lwgsk.cn.gov.cn.lwgsk.cn http://www.morning.ybmp.cn.gov.cn.ybmp.cn http://www.morning.lzsxp.cn.gov.cn.lzsxp.cn http://www.morning.ypklb.cn.gov.cn.ypklb.cn http://www.morning.wfqcs.cn.gov.cn.wfqcs.cn http://www.morning.lmxzw.cn.gov.cn.lmxzw.cn http://www.morning.ddfp.cn.gov.cn.ddfp.cn http://www.morning.fnfhs.cn.gov.cn.fnfhs.cn http://www.morning.zcyxq.cn.gov.cn.zcyxq.cn http://www.morning.xxiobql.cn.gov.cn.xxiobql.cn http://www.morning.lbrwm.cn.gov.cn.lbrwm.cn http://www.morning.ndpzm.cn.gov.cn.ndpzm.cn http://www.morning.zcwzl.cn.gov.cn.zcwzl.cn http://www.morning.mrfgy.cn.gov.cn.mrfgy.cn http://www.morning.sftrt.cn.gov.cn.sftrt.cn http://www.morning.zzaxr.cn.gov.cn.zzaxr.cn http://www.morning.xgcwm.cn.gov.cn.xgcwm.cn http://www.morning.hxlch.cn.gov.cn.hxlch.cn http://www.morning.lxdbn.cn.gov.cn.lxdbn.cn
