CN113553689A - 用于仿真技术系统的方法和设备 - Google Patents
用于仿真技术系统的方法和设备 Download PDFInfo
- Publication number
- CN113553689A CN113553689A CN202110434499.7A CN202110434499A CN113553689A CN 113553689 A CN113553689 A CN 113553689A CN 202110434499 A CN202110434499 A CN 202110434499A CN 113553689 A CN113553689 A CN 113553689A
- Authority
- CN
- China
- Prior art keywords
- following features
- simulation
- distribution function
- function
- measurement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Computational Biology (AREA)
- Operations Research (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
用于仿真技术系统的方法(10),其特征在于以下特征:借助于所述系统的仿真模型获得(11)时间序列,其中为所述仿真模型的至少一个认知参数分派能变化的值,通过对所述系统的对应测量获得(12)至少一个测量序列,对于所述认知参数的每个值,分别计算(13)针对该值获得的时间序列相对于所述测量序列的实数值误差度量(15)以及所述误差度量(15)的分布函数(16),以及为了所述仿真使用(14)以下值,在所述值的情况下所述分布函数(16)与郝维赛德函数的距离最小。
Description
技术领域
本发明涉及一种用于仿真技术系统的方法。本发明还涉及一种对应的设备、一种对应的计算机程序以及一种对应的存储介质。
背景技术
在软件技术中,用上位概念“基于模型的测试”(model-based testing,MBT)汇总了将模型用于对测试活动自动化并在测试过程中生成测试假象(Testartefakte)的用途。例如,由描述待测试系统的目标行为的模型生成测试情况是众所周知的。
特别是嵌入式系统(embedded systems)依赖于传感器的确定性输入信号,并且又通过向极为不同的执行器输出的信号来刺激其环境。因此,在这种系统的验证和前期开发阶段的过程中,对该系统的模型(model in the loop,循环中的模型,MiL)、软件(softwarein the loop,循环中的软件,SiL)、处理器(processor in the loop,循环中的处理器,PiL)或整个硬件(hardware in the loop,循环中的硬件,HiL)与环境模型一起在伺服控制中进行仿真。在车辆技术中,有时将与该原理相对应的用于根据测试阶段和测试对象检查电子控制设备的仿真器称为组件检查台、模块检查台或集成检查台。
DE10303489A1公开了这种用于测试车辆、电动工具或机器人系统的控制单元的软件的方法,其中通过测试系统至少部分地仿真可由所述控制单元控制的受控段,其方式是从所述控制单元产生输出信号并且将所述控制单元的输出信号经由第一连接传输到第一硬件模块,将第二硬件模块的信号作为输入信号经由第二连接传输到所述控制单元,其中所述输出信号作为软件中的第一控制值提供并且附加地经由通信接口相对于所述受控段实时地传输到所述测试系统。
这种仿真在各种技术领域中得到广泛应用,并且例如用于检查电动工具中的嵌入式系统、用于驱动、转向和制动系统的电机控制设备、相机系统、具有人工智能和机器学习组件的系统、机器人系统或自主车辆在它们的早期开发阶段中的适用性。然而,由于缺乏对根据现有技术的仿真模型结果的可靠性的信任,根据现有技术的仿真模型结果仅有限地包含在释放决策中。
发明内容
本发明提供了根据独立权利要求的用于仿真技术系统的方法、对应的设备、对应的计算机程序以及存储介质。
根据本发明的方案基于以下认识:验证来自仿真模型的时间信号涉及所述仿真模型的输出与来自试验的测量值之间的比较。所选择的仿真信号或测量信号在此描述了感兴趣量(quantity of interest,QOI),其可以是标量或以时间序列存在。
所提出的解决方案还考虑到以下事实:基于测量来验证仿真是最简单的验证情况。为此将两个标量或两个时间序列相互比较。然而,在实践中大多追求将多次重复试验与多次重复的仿真进行比较,以减少验证中的统计不确定性。对于标量来说,例如可以借助于所谓的可信区间或置信区间来进行的统计比较是常见的。此外,对于标量QOI来说提出了借助于所谓的概率盒(probability boxes,p-box)进行验证。这是两个累积分布函数(cumulative distribution functions,CDF)的比较。
时间序列的验证(例如,关于可能的相位误差或数量级误差)可以通过使用信号度量或验证度量从标量的验证中得出。这是数学函数,其将两个时间序列映射为一个标量并在下面称为误差度量。相关的误差度量例如是如下定义的均方误差(mean squared error,MSE):
在此,基本思想是以标量的方式表达仿真时间序列组和参考时间序列组之间的一致性或偏差。为此,基于任意误差度量将前者与后者进行成对比较。然后将这些成对比较的结果解释为仿真与参考数据之间差异的经验概率分布。
假设存在N1个不同的时间序列,它们是在所述仿真中检测的,必要时通过改变各个参数检测的。此外存在N2个时间序列,这些时间序列应当用作参考,例如来自对被仿真系统进行的测量的试验结果。现在首先考察第一仿真时间序列,并先后与所有参考时间序列进行比较,其中使用在这种情况下似乎合适的误差度量。对应地处理仿真数据组的所有其余时间序列。由此得到了所述误差度量的N1⋅N2个函数值。这些结果的分布可以通过绘制其分布函数的直方图或函数图来可视化;这种图示是普遍使用的,并且可以由专业观察者毫无问题地解释。
如果应当将这类结果用于验证,则可以对应地制定合适的质量标准作为表明仿真不合适的那些结果的最大份额。从而为了确定有效性,可能例如要求允许最大10%的比较得到大于1的误差。
根据本发明的所述方法的扩展具有利用不同类型的参数的优点:在试验中自然出现的偶然参数,以及尽管尚不确定,但是仍然应当得到确定的认知参数。基于此一方面开辟了估计所述模型的认知参数的简单的可能性,即使用参考数据组对所述模型的认知参数进行校准,不过另一方面也允许识别在所述仿真或所述参考数据组中的异常值。然后又可以将这些异常值从数据库中去除。
一般而言,可以通过这种方式扩展常规的MBT方法,以考虑所述仿真固有的不同不确定性。
通过从属权利要求中列出的措施可以对独立权利要求中说明的基本思想进行有利的扩展和改进。从而可以用于自动驾驶和其他自动化系统的领域,例如机器人。特别是在这些应用领域中,对嵌入式系统的根据本发明的确认和验证可以持久地提高所述嵌入式系统的功能可靠性。
所提出的方法在此被证明是非常灵活的,因为可以使用不同的误差度量来提供关于仿真和参考之间的平均偏差的信息,从所述平均偏差的分布中可以计算出诸如平均值、中位数和方差的统计位置和散射参数。
附图说明
在附图中示出了本发明的实施例,并且在下面的描述中详细解释了本发明的实施例。
图1示出了根据第一实施方式的方法的流程图。
图2示出了第一曲线族。
图3示出了第二曲线族。
图4示意性地示出了根据第二实施方式的控制设备。
具体实施方式
图1图解了根据本发明的方法10的基本步骤。首先借助于仿真模型获得N1个不同的时间序列(过程11)。如果改变偶然参数以及认知参数两者,则可以将仿真结果汇总在m1×n1矩阵中,在该矩阵中例如这些行包含带有固定的认知参数组的偶然参数。因此,存在m1个不同的认知参数组合—每个认知参数组合具有n1个偶然可变的参数。
然后通过对系统的对应测量,获得与时间序列相对应的测量序列(过程12),这些测量序列应当用作对所述仿真模型进行验证或校准的参考基础。所述测量序列的每个条目与每个参考时间序列的比较为所述矩阵近似添加一个其他维度。现在提供两种选择可能性来考虑以下事实:将附加维度作为偶然属性或认知属性来考察。通过这些替代的考察方式可以得到m1×(n1⋅N2)矩阵或(m1⋅N2)×n1矩阵。
下面展示了将所述附加维度视为偶然的理由。在数学上,可以通过以下方式对认知参数和偶然参数的变化建模:假设f(a,b)是取决于认知参数a和偶然参数b的仿真时间序列。因此,可以将b视为具有特定概率分布Pb的随机变量。因此,f也是随机的时间序列,其附加地取决于参数a。换句话说:b是随机的,而a表示有针对性地影响所述仿真模型的可能性。
参考g同样表示从中抽取了样本的未知分布的随机时间序列。如果给定误差度量d,则可以确定(过程13)。从上面的讲述中可以看出,这也是随机变量,更确切地说是针对认知参数a的每个值的随机变量族。来自这些变量的概率分布的样本通过上面提到的m1×(n1⋅N2)矩阵给出。
为了改善所述仿真,优选使用以下值(过程14):对于这些值来说分布函数16与郝维赛德(Heaviside)函数或单位跳跃函数,其中x<0时,x≥0时之间具有最小距离,因为对于认知参数的该值来说,所述仿真模型可以最准确地重现真实系统的特性。
在根据图2的示例中,这是以下分布函数16,其函数图设有附图标记17。例如,科学计算中常用的面积验证度量(area validation metric)适合作为用于确定该分布函数16的距离度量。
另一方面,如果将所述矩阵的附加维度分配给认知参数,则可以基于对应的图示识别出所谓的异常值。这是通过以下方式实现的,即,每个测量序列现在也都反映在自己的函数图中。因此,为非典型测量序列的误差度量15计算的分布函数16(该分布函数的函数图在根据图3的示例中设有附图标记18)即使在短暂考察时也通过其与其余曲线族的明显可识别的偏移而脱颖而出。
如图4中的示意图所示,该方法10可以例如以软件或硬件或者以软件和硬件的混合形式例如在工作站20中实现。
Claims (10)
1.一种用于仿真技术系统的方法(10),
其特征在于以下特征:
借助于所述系统的仿真模型获得(11)时间序列,其中为所述仿真模型的至少一个认知参数分派能变化的值,
通过对所述系统的对应测量获得(12)至少一个测量序列,
对于所述认知参数的每个值,分别计算(13)针对该值获得的时间序列相对于所述测量序列的实数值误差度量(15)以及所述误差度量(15)的分布函数(16),以及
为了所述仿真使用(14)以下值,在所述值的情况下所述分布函数(16)与郝维赛德函数的距离最小。
2.根据权利要求1所述的方法(10),
其特征在于以下特征:
通过所述测量获得多个彼此对应的测量序列,
对每个所述测量序列计算(13)所述误差度量(15)和所述分布函数(16),以及
从所述测量序列中删除以下测量序列,在所述测量序列的情况下所述分布函数(16)与其余分布函数(16)的距离超过预给定距离。
3.根据权利要求1或2所述的方法(10),
其特征在于以下特征:
所述距离是根据面积验证度量确定的。
4.根据权利要求1至3中任一项所述的方法(10),
其特征在于以下特征:
示出所述分布函数(16)的函数图。
5.根据权利要求1至4中任一项所述的方法(10),
其特征在于以下特征:
所述误差度量(15)是均方误差。
6.根据权利要求1至5中任一项所述的方法(10),
其特征在于以下特征:
所述系统是嵌入式系统,并且
借助于所述仿真来检查所述系统的功能可靠性。
7.根据权利要求1至6中任一项所述的方法(10),
其特征在于以下特征之一:
所述系统被设置为自主地控制车辆,或
所述系统被设置为控制制造过程。
8.一种计算机程序,其被设置为执行根据权利要求1至7中任一项所述的方法(10)。
9.一种机器可读存储介质,其上存储有根据权利要求8所述的计算机程序。
10.一种设备(20),其被设置为在集成的测试和仿真环境中执行根据权利要求1至7中任一项所述的方法(10)。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102020205131.3 | 2020-04-23 | ||
DE102020205131.3A DE102020205131A1 (de) | 2020-04-23 | 2020-04-23 | Verfahren und Vorrichtung zum Simulieren eines technischen Systems |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113553689A true CN113553689A (zh) | 2021-10-26 |
Family
ID=78101792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110434499.7A Pending CN113553689A (zh) | 2020-04-23 | 2021-04-22 | 用于仿真技术系统的方法和设备 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113553689A (zh) |
DE (1) | DE102020205131A1 (zh) |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10303489A1 (de) | 2003-01-30 | 2004-08-12 | Robert Bosch Gmbh | Verfahren und Vorrichtung zum Testen von Software einer Steuereinheit eines Fahrzeugs |
-
2020
- 2020-04-23 DE DE102020205131.3A patent/DE102020205131A1/de active Pending
-
2021
- 2021-04-22 CN CN202110434499.7A patent/CN113553689A/zh active Pending
Also Published As
Publication number | Publication date |
---|---|
DE102020205131A1 (de) | 2021-10-28 |
US20210334435A1 (en) | 2021-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107481292B (zh) | 车载摄像头的姿态误差估计方法和装置 | |
EP2715459B1 (en) | Machine-implemented method for obtaining data from a nonlinear dynamic real system during a test run | |
CN113590456A (zh) | 用于检查技术系统的方法和设备 | |
CN107972070A (zh) | 机器人性能的测试方法、测试系统及计算机可读存储介质 | |
US20200409823A1 (en) | Method and apparatus for optimal distribution of test cases among different testing platforms | |
US20210216609A1 (en) | Degradation detection system | |
Matinnejad et al. | Automated model-in-the-loop testing of continuous controllers using search | |
US11397660B2 (en) | Method and apparatus for testing a system, for selecting real tests, and for testing systems with machine learning components | |
CN111612048A (zh) | 无监督聚类异常检测的方法 | |
Zhang et al. | Neural network based uncertainty prediction for autonomous vehicle application | |
Khlamov et al. | The astronomical object recognition and its near-zero motion detection in series of images by in situ modeling | |
CN113553689A (zh) | 用于仿真技术系统的方法和设备 | |
Vaca et al. | Reverse engineering for electronic devices: temperature and soil moisture case | |
CN113704085A (zh) | 用于检查技术系统的方法和设备 | |
CN113590458A (zh) | 用于检查技术系统的方法和设备 | |
US12032880B2 (en) | Method and device for simulating a technical system | |
CN115270902A (zh) | 用于测试产品的方法 | |
CN112146890A (zh) | 用于评价和选择信号比较度量的方法和设备 | |
CN113722207A (zh) | 用于检查技术系统的方法和装置 | |
Khlamov et al. | Statistical Modeling for the Near-zero Apparent Motion Detection of Objects in Series of Images From data Stream | |
CN114651190A (zh) | 用于批准使用检测车辆的环境中的物体的传感器系统的方法、设备和计算机程序 | |
Cramer et al. | Measurement Uncertainty: Relating the uncertainties of physical and virtual measurements | |
US20210026999A1 (en) | Method and device for validating a simulation of a technical system | |
Mahajan et al. | Quantifying the sim2real gap for GPS and IMU sensors | |
CN110879530A (zh) | 用于对动态模型进行安全的训练的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |