局部敏感性分析
局部敏感性分析(local sensitivity analysis)研究模型在指定基点附近,输入的微小变化会引起多大的输出变化。对可微模型,它用偏导数或雅可比矩阵表示一阶响应。这个响应依赖基点、保持不变的条件和变量单位;它既不是整个输入范围内的重要性排名,也不自动给出观测数据中的因果关系。[1]
串联电路中的一次变化
取一个理想直流电路,电源电压为 ,两个固定温度下的欧姆电阻 串联,忽略电源内阻。电流模型为 以下用 、、 作为教学基点,得到 。要判断第一只电阻增加一点时的变化,先保持电压和第二只电阻不变: 负号表示电阻增加会使电流下降,数值表示在该基点附近,每增加一欧姆,电流约下降 。这是切线斜率的解释,不是对任意增量都精确成立的比例律。
例如把 从十欧姆改为十一欧姆,一阶预测为 重新计算模型则为 ,二者相差约 。下图在同一基点比较原曲线和切线;离开基点越远,曲率造成的差别越明显。
导数、线性化与余项
一般地,令输入向量为 ,标量输出为 ,基点 位于定义域内部。若 在该点可微,则 这里的小 o 表示余项除以扰动长度后趋于零。只存在各个偏导数并不足以保证这个联合近似;可微性才保证所有足够小方向上的线性化。输出是多维向量时,将各输出对各输入的偏导数排成雅可比矩阵,一阶变化写成 。
在电路例子中,三个绝对敏感度是 因此多个小变化同时发生时,预测电流增量为 其中电压和电阻的数值分别按伏和欧姆代入。系数有不同单位,不能只比较它们的数字大小就判断哪个输入“更重要”。
这个例子还能精确算出一种余项。设总电阻原为 ,电压固定,总电阻改变 ,且 。同分母展开得 将右边与实际输出 相除,在 时,线性近似的绝对相对误差恰为 。前例为 ,约百分之零点一一。误差为二次量,解释了足够小扰动时一阶近似的作用;当总电阻接近零时,原模型本身也显露奇异性。
相对敏感度与单位变换
若基点的 且 ,定义相对敏感度或弹性 它没有单位,将一阶关系改写成 当参数和输出均为正时,也可写成 。这不是任意有限百分比变化的精确等式。
电路中的相对敏感度为 所以电压单独增加百分之一,电流约增加百分之一;第一、第二只电阻各自单独增加百分之一,电流分别约下降百分之三分之一和百分之三分之二。两只电阻对“增加同样一欧姆”的响应相等,但对“增加同样百分比”的响应不同,因为原始电阻并不相等。
把电阻的数值单位从欧姆改为千欧姆,绝对导数会变为原来的千倍,而参数数值变为原来的千分之一,乘积 不变。因此相对敏感度对这种乘法单位换算不变。若输出或参数基值为零,前面的除法没有定义,应报告绝对导数,或明确另选的参考尺度;不能把无定义的比值填写为零。
方向、约束与不确定度
偏导数固定其他输入。若实际变化受到约束,必须沿允许的方向计算。例如保持总电阻三十欧姆,令 ,则 两种结论并不冲突:一欧姆的单独增加与两只电阻之间的一欧姆转移,是不同的实验。
若输入在基点附近随机变化,协方差矩阵为 ,则在线性近似有效、相关二阶矩存在的条件下,有常用的不确定度传播近似 这是把线性化代入方差公式所得,忽略了非线性余项;不能将其当作任意大输入不确定度下的精确式。协方差与线性近似的适用性也是不确定度传播中需要保留的条件。[2]例如只让两只电阻有独立的小波动,标准差各为 ,电流标准差近似为 若两电阻波动始终等大反向,总电阻不变,此时协方差交叉项恰好抵消;漏掉它们就会错误地预测电流波动。
最后,局部导数为零只说明一阶项为零。例如 在零处导数为零,但有限改变 会产生 的输出。对较宽范围、非单调关系或相互作用,需要全局敏感性分析。方差贡献的概率定义另见索博尔指数。这些方法回答的问题不同,不能直接用局部导数替代。
参考文献
- ↑ European Commission, Joint Research Centre,Methods,Global Sensitivity Analysis 部分,对局部导数与全局变化范围的区分。
- ↑ NIST,NIST TN 1297: 5. Combined Standard Uncertainty,第5.1及5.4节,关于协方差、标准差组合与线性近似的说明。