跳到正文
格致开物MATHWIKI

局部敏感性分析

局部敏感性分析(local sensitivity analysis)研究模型在指定基点附近,输入的微小变化会引起多大的输出变化。对可微模型,它用偏导数或雅可比矩阵表示一阶响应。这个响应依赖基点、保持不变的条件和变量单位;它既不是整个输入范围内的重要性排名,也不自动给出观测数据中的因果关系。[1]

串联电路中的一次变化

取一个理想直流电路,电源电压为 V,两个固定温度下的欧姆电阻 R1,R2>0 串联,忽略电源内阻。电流模型为 I(V,R1,R2)=VR1+R2. 以下用 V=12V、R1=10Ω、R2=20Ω 作为教学基点,得到 I0=0.4A。要判断第一只电阻增加一点时的变化,先保持电压和第二只电阻不变: ∂I∂R1=−V(R1+R2)2=−175A/Ω. 负号表示电阻增加会使电流下降,数值表示在该基点附近,每增加一欧姆,电流约下降 0.0133333A。这是切线斜率的解释,不是对任意增量都精确成立的比例律。

例如把 R1 从十欧姆改为十一欧姆,一阶预测为 Ilin=0.4−175×1≈0.3866667A. 重新计算模型则为 Iexact=12/31≈0.3870968A,二者相差约 0.0004301A。下图在同一基点比较原曲线和切线;离开基点越远,曲率造成的差别越明显。

电流随第一只电阻变化的曲线与在十欧姆处的切线,十一欧姆处两者数值分别约0.38710和0.38667安培
电压固定为十二伏,第二只电阻固定为二十欧姆。标出的两个近邻值对应正文的一欧姆扰动。

导数、线性化与余项

一般地,令输入向量为 p=(p1,…,pd),标量输出为 f(p),基点 p0 位于定义域内部。若 f 在该点可微,则 f(p0+δp)=f(p0)+∑i=1dSiδpi+o(‖δp‖),Si=∂f∂pi(p0). 这里的小 o 表示余项除以扰动长度后趋于零。只存在各个偏导数并不足以保证这个联合近似;可微性才保证所有足够小方向上的线性化。输出是多维向量时,将各输出对各输入的偏导数排成雅可比矩阵,一阶变化写成 Jδp。

在电路例子中,三个绝对敏感度是 SV=130A/V,SR1=SR2=−175A/Ω. 因此多个小变化同时发生时,预测电流增量为 δI≈δV30−δR1+δR275, 其中电压和电阻的数值分别按伏和欧姆代入。系数有不同单位,不能只比较它们的数字大小就判断哪个输入“更重要”。

这个例子还能精确算出一种余项。设总电阻原为 R>0,电压固定,总电阻改变 δ,且 R+δ>0。同分母展开得 VR+δ−(VR−VδR2)=Vδ2R2(R+δ). 将右边与实际输出 V/(R+δ) 相除,在 V≠0 时,线性近似的绝对相对误差恰为 δ2/R2。前例为 1/900,约百分之零点一一。误差为二次量,解释了足够小扰动时一阶近似的作用;当总电阻接近零时,原模型本身也显露奇异性。

相对敏感度与单位变换

若基点的 pi≠0 且 f(p0)≠0,定义相对敏感度或弹性 Ei=pif(p0)Si. 它没有单位,将一阶关系改写成 δff(p0)≈∑iEiδpipi. 当参数和输出均为正时,也可写成 Ei=∂ln⁡f/∂ln⁡pi。这不是任意有限百分比变化的精确等式。

电路中的相对敏感度为 EV=1,ER1=−13,ER2=−23. 所以电压单独增加百分之一,电流约增加百分之一;第一、第二只电阻各自单独增加百分之一,电流分别约下降百分之三分之一和百分之三分之二。两只电阻对“增加同样一欧姆”的响应相等,但对“增加同样百分比”的响应不同,因为原始电阻并不相等。

三个有正负号的相对敏感度柱,电压为1,第一只电阻为负三分之一,第二只电阻为负三分之二
柱长比较相同相对输入变化的一阶效应;正负号保留响应方向。

把电阻的数值单位从欧姆改为千欧姆,绝对导数会变为原来的千倍,而参数数值变为原来的千分之一,乘积 piSi 不变。因此相对敏感度对这种乘法单位换算不变。若输出或参数基值为零,前面的除法没有定义,应报告绝对导数,或明确另选的参考尺度;不能把无定义的比值填写为零。

方向、约束与不确定度

偏导数固定其他输入。若实际变化受到约束,必须沿允许的方向计算。例如保持总电阻三十欧姆,令 R2=30−R1,则 dIdR1=∂I∂R1+∂I∂R2dR2dR1=−175+175=0. 两种结论并不冲突:一欧姆的单独增加与两只电阻之间的一欧姆转移,是不同的实验。

若输入在基点附近随机变化,协方差矩阵为 Σ,则在线性近似有效、相关二阶矩存在的条件下,有常用的不确定度传播近似 Var⁡(f)≈JΣJ𝖳. 这是把线性化代入方差公式所得,忽略了非线性余项;不能将其当作任意大输入不确定度下的精确式。协方差与线性近似的适用性也是不确定度传播中需要保留的条件。[2]例如只让两只电阻有独立的小波动,标准差各为 0.3Ω,电流标准差近似为 0.32+0.32752≈0.005657A. 若两电阻波动始终等大反向,总电阻不变,此时协方差交叉项恰好抵消;漏掉它们就会错误地预测电流波动。

最后,局部导数为零只说明一阶项为零。例如 f(x)=x2 在零处导数为零,但有限改变 x 会产生 x2 的输出。对较宽范围、非单调关系或相互作用,需要全局敏感性分析。方差贡献的概率定义另见索博尔指数。这些方法回答的问题不同,不能直接用局部导数替代。

参考文献

  1. ↑ European Commission, Joint Research Centre,Methods,Global Sensitivity Analysis 部分,对局部导数与全局变化范围的区分。
  2. ↑ NIST,NIST TN 1297: 5. Combined Standard Uncertainty,第5.1及5.4节,关于协方差、标准差组合与线性近似的说明。