基于仿真的自动驾驶可靠性估计(二)
前 言
基于仿真的自动驾驶可靠性估计(一)中已经介绍,使用定步长泛化、朴素蒙特卡罗等方法生成验证自动驾驶系统的仿真场景难以在可以接受的成本内精确估计被测试系统在指定逻辑场景即测试空间内的失效概率。本文将由此出发,介绍若干可以用来估计罕见事件发生概率的可靠性分析方法。
01. 问题定义
我们将一个拥有D个可泛化参数的逻辑场景等价于 维参数空间 ,其中 是一组随机变量,其分布函数由逻辑场景决定,而 的一组具体取值 即为此逻辑场景下的一个具体场景。
被测试自动驾驶系统或车辆在某一具体场景 中失效,指该系统及车辆的某项评价指标(kpi)处于失效域中。我们用 表示某项kpi在场景 上的取值,不失一般性地用 表示此项kpi处于失效域中。故被测试系统或车辆在某一逻辑场景中的失效概率等价于对应测试空间中,参数组合落入失效域的概率,即
其中为此逻辑场景中参数的联合分布。
本文剩余部分将介绍若干精确高效估计此概率的可靠性分析方法并通过数值实现说明其效率。一般我们会在标准正态空间中使用这些可靠性分析方法,故我们提前使用Rosenblatt变换将测试空间转化为标准正态分布空间,之后的讨论均基于标准正态分布空间。
02.一阶可靠性分析方法(FORM)
一种经典的可靠性分析算法是将 在失效临界面(即 的区域)上的一个点 利用泰勒展开进行线性近似,即 再通过 来估计 。如下图所示,右上角的区域为失效域,FORM用线性展开即黄色线以上的区域来代替实际失效区域。
显然失效临界的形状及 的选取对近似的结果有很大的影响,直观想来我们会选取失效临界中最可能的点(实际是似然函数最大的点)MPP做为展开的中心 ,即 此时利用正态分布的性质及 的约束可以证明 其中 为标准正态分布函数的概率分布函数, 为原点到MPP的距离。可以看出FORM的精度取决于失效临界与线性边界的差距,而FORM算法的效率只取决于寻找MPP算法的速度。最早使用Rackwith-Fiessler算法来快速寻找MPP点,然而此算法并不能保证收敛,故在Rackwith-Fiessler算法不收敛时可以使用更为复杂的优化算法来实现,例如NLPQL。
03. 方向采样算法
将测试空间考虑成一个多维球体,若在每个方向上存在最多一个失效临界时,可以通过方向采样算法估计失效概率 。其基本原理是
1. 在半径为 的层球面进行均匀采样,并记录其中的失效场景;
2. 对失效场景所在的方向插值求解此方向的失效临界场景;
3. 注意测试空间为 维标准正态分布空间,故其中任意一点与原点距离的平方服从自由度为的卡方分布( ),故可对失效概率估计如下
其中 是总采样数, 是失效方向数, 是第 个失效方向的临界场景与原点间的距离。
通过上面的介绍可以看出,方向采样算法无法观察到距离原点超过R的失效场景。在标准正态分布空间中无法观测的失效场景概率占比小于等于 ,故往往 取比较大的值例如 至 ,如此漏掉的失效场景占比仅在 至 ,是完全可以接受的。影响此算法精度的另一个因素是在外层球面采样的均匀程度,我们使用球面拉丁超立方采样代替传统的方法获取高维球面上的超均匀样本,用相同数量的样本实现对球面更好的覆盖。
三维球面的均匀采样,左侧为传统方法,右侧为球面拉丁超立方采样,红色为寻找到的危险场景
04. 重要性采样算法
重要性采样算法。其中蓝色点为安全场景,黄色点为本轮采样中较危险场景,红色点为失效场景,紫色点为失效临界场景。其中左上和右小角的是小边界为线性,左下和右上角的失效边界为非线性。
05. 数值实验
下面将用2个数值实验,证明不同可靠性分析算法的有效性。
数值实验一:
可靠性分析方法 | FORM |
测试空间 | |
kpi函数 | |
理论失败概率 |
由于kpi函数的连续性及凸凹性较好,使用Rackwith-Fiessler算法快速迭代寻找MPP,得到如下实验结果
失败概率估计 | |
相对偏差 | 16% |
相对误差 | 0,FORM方法不使用随机采样 |
kpi函数调用次数 | 5次,5次迭代寻找MPP |
可以看当线性展开可以比较好的逼近失效临界是,FORM方法可以很好的估计失效概率。此例中虽然无法在全局使用线性展开逼近失效临界,但在MPP点附近逼近效果较好,即在对失效概率贡献最大的区域逼近效果较好,故可以得到较好的估计效果。
数值实验二
实验二中我们使用 测试空间中的Katsuki函数作为kpi函数,即 Katsuki函数通过是具备了不同形状的失效临界,其理论失效概率为 。使用方向采样及重要性采样对其进行估计,下面将展示两种算法的结果
此实验中方向采样算法的球面采样数为100,临界面寻找的最大迭代次数为10;重要性采样每轮采样数为500。因为测试空间维度较低,方向采样有着非常优秀的表现,同时两者对比蒙特卡罗方法都体现出了超高的效率。需要注意的是,虽然重要性采样使用的kpi函数平均次数更多,但在并行仿真下实际只有4轮仿真;而kpi调用次数较少的方向采样因为缺乏并行仿真支撑能力,耗时反而更长。
最后对不同的可靠性分析算法能够适应的情景做出分析,需要注意的是蒙特卡罗方法或拉丁超立方/sobol方法虽然不受测试空间的限制,但其效率仅能应对概率在 及更大的情况。
-
汽车测试网V课堂
-
微信公众号
-
汽车测试网手机站
最新资讯
-
直播|中国汽研:重型车温室气体企业平均方
2024-11-23 17:42
-
美国FCC正式划定5.9GHz频段用于C-V2X技术
2024-11-23 17:41
-
标准立项 | 《自主代客泊车 预期功能安全测
2024-11-23 16:13
-
长三角(盐城)智能网联汽车试验场气象多级
2024-11-23 16:11
-
CSAE汽车可靠性技术分会第二届委员会第二次
2024-11-23 16:09