用于实时计算过程模型状态量的方法和模拟器.pdf

描述和示出了一种用于在过程计算机(1)上实时计算过程模型的状态量(xM)的方法，所述过程计算机具有第一计算单元(2)和处理接口(3)，其中所述过程计算机(1)被设置为能够由所述过程计算机(1)通过所述处理接口(3)采集物理过程(4)的至少一个状态量(xP)，和/或能够由所述过程计算机(1)输出至少一个输出量(y)以影响所述物理过程(4)，其中由所述第一计算单元(2)利用显式积分法计算所述过程模型，以确定状态量(xM，e)。本发明的目的是至少部分地避免在用于实时计算过程模型状态量的已知方法中的缺点，尤其是在最大程度确保计算的数值精度和稳定性的前提下能够保证计算的实时性。根据本发明，该任务这样来解决：提供了第二计算单元(5)，基本上与利用所述第一计算单元(2)对过程模型进行计算同步地，利用这个第二计算单元(5)通过隐式积分法执行该过程模型，为了校正在第一计算单元(2)上所执行的对过程模型的计算，调用利用第二计算单元(5)和隐式积分法所计算的过程模型的状态量(xM，i)，并将以这种方式利用显式积分法求得的状态量(xM，e)用作该过程模型的状态量(xM)。

1.  一种用于在过程计算机(1)上实时计算过程模型的状态量(x_M)的方法，所述过程计算机具有第一计算单元(2)和处理接口(3)，其中所述过程计算机(1)被设置为能够由所述过程计算机(1)通过所述处理接口(3)采集物理过程(4)的至少一个状态量(x_P)，和/或能够由所述过程计算机(1)输出至少一个输出量(y)以影响所述物理过程(4)，其中由所述第一计算单元(2)利用显式积分法计算所述过程模型，以确定状态量(x_M，e)，其特征在于，
提供了第二计算单元(5)，基本上与利用所述第一计算单元(2)对过程模型进行计算同步地，利用这个第二计算单元(5)通过隐式积分法执行该过程模型，为了校正在第一计算单元(2)上所执行的对过程模型的计算，调用利用第二计算单元(5)和隐式积分法所计算的过程模型的状态量(x_M，i)，并将以这种方式利用显式积分法求得的状态量(x_M，e)用作该过程模型的状态量(x_M)。

2.  根据权利要求1的方法，其特征在于，为了用第一计算单元(2)计算当前的状态量(x_M，e)，使用由第二计算单元(5)所计算的已经过去的状态量(x_M)——如果它存在的话。

3.  一种用于在过程计算机(1)上实时计算过程模型的状态量(x_M)的方法，所述过程计算机具有第一计算单元(2)和处理接口(3)，其中所述过程计算机(1)被设置为能够由所述过程计算机(1)通过所述处理接口(3)采集物理过程(4)的至少一个状态量(x_P)，和/或能够由所述过程计算机(1)输出至少一个输出量(y)以影响所述物理过程(4)，其中由所述第一计算单元(2)利用显式积分法计算所述过程模型，以确定状态量(x_M，e)，其特征在于，
提供了第二计算单元(5)，基本上与利用所述第一计算单元(2)对过程模型进行计算同步地，利用这个第二计算单元(5)通过隐式积分法计算该过程模型，只要利用第二计算单元(5)和隐式积分法所计算的状态量(x_M，i)是实时的，就将利用第二计算单元(5)和隐式积分法所计算的状态量(x_M，i)用作该过程模型的状态量(x_M)，否则的话，将利用显式积分法所求得的状态量(x_M，e)用作该过程模型的状态量(x_M)。

4.  根据权利要求1至3中任一项的方法，其特征在于，在第二计算单元(5)上实现的隐式积分法至少部分地用比在第一计算单元(2)上实现的显式积分法更大的采样步宽来工作，尤其是用可变的采样步宽来工作。

5.  根据权利要求4的方法，其特征在于，在第二计算单元(5)上，位于两个计算出的采样时刻之间的状态量(x_M，i)通过内插来计算，尤其是通过过程模型或过程模型一部分的多项式拟合函数来计算，特别是通过密集输出法来计算。

6.  根据权利要求1至5中任一项的方法，其特征在于，在第二计算单元(5)上，为计算过程模型的将来的状态量(x_M，i^k+1)所需的该物理过程的未来状态量(x_P^k+1)通过由该物理过程的已经过去的状态量的外插来得到。

7.  根据权利要求1至6中任一项的方法，其特征在于，当第一计算单元(2)和第二计算单元(5)对过程模型的状态量(x_M)的计算延迟结束的情况下，该过程模型的状态量(x_M)通过由已经求得的该过程模型的状态量(x_M)的先前的值的外插来计算。

8.  根据权利要求1至7中任一项的方法，其特征在于，过程模型利用所述第二计算单元(5)和至少一个另外的第三计算单元上的隐式积分法并行处理，并且在所述第二计算单元(5)和所述至少一个另外的第三计算单元上并行地共同计算该过程模型的状态量(x_M，i)。

9.  一种模拟器，尤其是在环模拟器，用于实时计算过程模型的状态量(x_M)，具有第一计算单元(2)、至少一个第二计算单元(5)和处理接口(3)，其中能够由所述模拟器通过所述处理接口(3)采集物理过程(4)的至少一个状态量(x_P)，并且能够由所述模拟器输出至少一个状态量(y)以影响所述物理过程，其中所述模拟器被设置为能够利用它来实现如权利要求1至8中任一项所述的方法。

用于实时计算过程模型状态量的方法和模拟器
技术领域
本发明涉及一种用于在过程计算机上实时计算过程模型状态量的方法和模拟器，所述过程计算机具有第一计算单元和处理接口，其中所述过程计算机被设置为能够由所述过程计算机通过所述处理接口采集物理过程的至少一个状态量，和/或能够由所述过程计算机输出至少一个输出量以影响所述物理过程，其中由所述第一计算单元利用显式积分法计算过程模型，以确定状态量。另外，本发明还涉及另外一种该类型的方法和用于实现这些方法的模拟器。
背景技术
实践中，先前所给出的用于实时计算过程模型状态量的方法以及实现这类方法的模拟器尤其是不能离开机电和过程自动化领域来考虑。通常，过程计算机或模拟器通过其处理接口与真实的物理过程相联系，使得物理过程的状态量能够通过测量由过程计算机或模拟器所采集，并且过程计算机或模拟器能够经由该处理接口通过输出量的输出对真实的物理过程施加影响。
对于过程计算机或模拟器有大量的应用场合，其例如可以根据在过程计算机或模拟器上执行的是何种类型的过程模型来进行分类。
当在过程计算机上主要实现的是控制算法的时候，过程计算机通常是控制设备，其中过程模型可以是所实现的控制算法的一部分，例如在从线性或者还有非线性的多个量的系统去耦时就是这种情况。一种典型的应用场合在这里是调节器设计，其例如根据所采用的开发工具、例如在所谓的“快速控制原型(RCP)”的框架内实现。
在其它经常出现的应用场合下，由过程计算机计算一个工艺流程的数学映像作为过程模型，如汽车发动机或传动机构的数学模型。这样的过程模型通过过程计算机的处理接口与一个“真实的”技术物理过程相关联，在上述汽车应用场合的情况下，往往与能够测试其功能的控制设备相关联；目前它可能是发动机控制设备。由过程模型的输出量对该控制设备施加相应的信号(例如发动机转数、凸轮轴的姿态等)，并且该控制设备以特定的状态量、即物理过程的状态量对其做出反应，其至少部分地通过处理接口反作用于过程计算机。在这种应用场合下，过程计算机经常被称为“模拟器”。模拟器也可以作为开放的控制回路“开环”工作，其中物理过程的状态量并不反作用于模拟器。
借助前面的例子可知，“过程计算机”和“模拟器”的概念是指等同的数字采样及计算工具，利用它能够计算出过程模型。下文中提到“过程计算机”或“模拟器”时，只是考虑到针对特定应用的常用语言习惯，但在功能上这些概念在下文中通常是可互换的。如果利用模拟器似乎使测试对象(这里是真实的过程)构成了一个工作循环，则经常用“在环(in-the-Loop)”模拟器来表示，在本例中是硬件在环模拟器。
在过程计算机上计算出的过程模型是通过微分方程对物理技术系统的数学描述。这些微分方程是过程模型的与时间有关的状态量，其中所述微分方程通过公共的状态量彼此关联。为了数值求解过程模型，即为了由微分方程求出状态量，采用了已知的数值积分法。该积分法可按照不同的观点来分类，例如可以分成显式积分法和隐式积分法，以及单步积分法和多步积分法。
首先，这些积分法的共同之处在于，仅仅在某些特定的离散时间点计算一个过程模型的状态量x_M，结果，这些状态量并不是作为时间连续的函数x_M(t)，而是仅仅位于离散的时间支持点x_M(t_k)。针对该过程模型的状态量x_M利用积分法在一个支持点求出的值在方法的限制下始终是有误差的，因此只能得到对于状态量x_M的实际值的近似值。对于各种数值积分均存在以下任务：即要求出状态量x_M在时刻t_k+1的将来值，即状态量x_M，k+1，其中一般地来说有下列公式：
x‾M,k=1=x‾M,k+f‾(x·M,k+1,xM,k,x·M,k,xM,k-1,...)---(1)]]>
显式积分法规定：在时刻t_k+1要新计算的状态量x_M的值、即x_M，k+1仅根据时间上位于前面的值x_M，k和等等来计算，即函数f可以显式地通过该值来表述。通过显式积分法对状态量x_M，e的计算在位于前面的离散时间点进行，即根据已经过去的状态量的值及其一次或多次求导递归地进行。根据状态量的值需要回溯多少个时间步骤，将其称为单步积分法或者多步积分法；这同样适用于隐式积分法。
隐式积分法规定：如公式1所示，函数f本身与当前所求得的状态量或者其一次或多次求导相关。已知的是，这种相关性在数值上只能迭代求解，因此这意味着：在数值积分步骤内必须执行一次迭代循环，而这在显式积分法中则不需要。通常当连续计算出的近似值彼此相差小于一个固定设置的最大值、即处于特定误差范围以下时结束为求取用隐式积分法所计算的状态量x_M，i的迭代。
显式积分法所具有的优点是：为计算当前状态量x_M，e所需的计算步骤数目是已知的，并且近似保持相等，从而能够确保状态量x_M的当前值位于特定的、可预见的时间点。这对于所有实时应用来说都是十分重要的，因为否则的话就不能实现时间离散的采样系统了。但显式积分法的缺点是：计算精度具有不确定性，尤其是在刚性系统、即具有显著不同的特征值的系统的情况下，当显式数值积分法的步宽被选择为足够小时，只能得到可接受的计算结果。对于刚性系统，显式积分法的步宽可能必须要缩小，这样就不能再在实时性条件下进行计算了。
与此相反，隐式积分法的优点是通过预先给定误差范围可以非常精确地预定计算精度，但相应的缺点是不知道在一个积分步骤内为了计算出所要求的精度需要多少次迭代。因此隐式积分法尤其是在刚性系统的情况下只能非常有限地适用于实时应用。
发明内容
因此，本发明的目的是至少部分地避免在用于实时计算过程模型状态量的已知方法中所提到的上述缺点，尤其是在最大程度确保计算的数值精度和稳定性的前提下能够保证计算的实时性。
首先根据本发明的第一个技术方案，本发明所述的用于实时计算与真实物理过程相关联的过程模型的状态量的方法的主要特征在于，提供了第二计算单元，基本上与利用第一计算单元计算过程模型同步地，利用所述第二计算单元通过隐式积分法执行该过程模型，为了校正在第一计算单元上所执行的对过程模型的计算，调用利用第二计算单元和隐式积分法所计算的过程模型状态量，并将以这种方式通过显式积分法所求得的状态量用作该过程模型的状态量。
根据本发明的方法利用了以下原理：过程计算机通过硬件配置提供了多个计算单元，其可能是多处理器系统的形式、即具有多个作为部件的可区分的处理器，或者可能是多内核处理器(Multicore)，其中一个处理器包含多个在很大程度上彼此独立工作的内核，它们实际上可以同时处理不同的进程，这些进程不是仅仅串行地先后执行，而是以单个线程的形式被处理，因而不会产生时间增益。
当说到“基本上”与利用第一计算单元计算过程模型同步地，利用所述第二计算单元通过隐式积分法执行该过程模型时，这意味着：所述两种积分法彼此同时地执行，但允许硬件和方法步骤上的差别。例如可能是这样的情况：显式和隐式的积分法在不同的、即彼此稍微有所偏差的计算时刻对状态量x_M，e及x_M，i执行计算，此时对过程模型的两种计算在本发明的意义上仍然可以认为是“基本上”同步执行的，因为这两者均涉及相同的物理时间，并具有共同的物理时基。
通过将利用显式积分法所求得的状态量x_M，e用作过程模型的状态量x_M，可以保持计算过程模型状态量时的实时条件。通过以利用隐式积分法所求得的状态量x_M，i来校正利用显式积分法所求得的状态量x_M，e，可以附带地保证显式积分法的误差不会随便累加。其中尽管尽可能经常采取这种校正是具有优点的，但并不要求在显式积分法的每一个计算步骤中都进行这种校正。
根据本发明的一个特别有利的实施例，为了用第一计算单元计算当前的状态量x_M，e，使用由第二计算单元所计算的已经过去的状态量x_M，i——如果它存在的话。在这里，显式积分法(只要可能的话)用隐式积分法的更为精确的数据来执行。当用隐式积分法所计算的状态量x_M，i目前不完整时(例如由于隐式积分法没有能够及时结束)，由显式积分法得到的肯定是带有错误的状态信息。
根据本发明的另一独立技术方案，开始提到的以及所描述的任务这样来解决：提供了第二计算单元，基本上与利用第一计算单元计算过程模型同步地，利用所述第二计算单元通过隐式积分法计算该过程模型，只要利用第二计算单元和隐式积分法所计算的状态量是实时的，就将利用第二计算单元和隐式积分法所计算的状态量用作该过程模型的状态量，否则的话，将利用显式积分法所求得的状态量用作该过程模型的状态量。
在这种本发明所述的方法中，也就是说，当隐式积分法没有及时结束的情况下，即在实时性方面没能得到所希望的结果时，在第一计算单元上所执行的显式积分法作为对该过程模型进行计算的返回点。对于在第二计算单元上执行的隐式积分法违背了实时性条件的情况，采用在第一计算单元上执行的显式积分法的计算结果作为替代。很显然，可以把隐式积分法的状态量x_M，i提供给第一计算单元上的显式积分法，从而可以保证显式积分法也能得到良好的精度。
在本发明所述方法的一种有利的实施例中，在第二计算单元上实现的隐式积分法至少部分地用比在第一计算单元上实现的显式积分法更大的采样步宽来工作，其中隐式积分法的采样步宽尤其优选的是可变的，即根据状态量x_M或者根据这些状态量随时间的变化而进行适配。在过程模型的状态量x_M必须位于固定的时间扫描栅内、即必须提供恒定采样率的应用场合下，隐式积分法在可变计算步宽下所提供的状态量x_M，i可以通过下述方法来同步，即在第二计算单元上，位于两个计算出的采样时刻之间的状态量通过内插来计算，尤其是通过过程模型或过程模型一部分的多项式拟合函数来计算，尤其是通过密集输出(Dense-Output)法来计算。当隐式积分法的采样步宽与显式积分法的采样步宽相差不太多时，在求取两个计算步骤之间过程模型的状态量时，简单的、即线性的内插已足以得到可用的结果。
当隐式积分法的采样步宽与显式积分法的采样步宽相差较大或者对隐式积分法的精度有较高要求时，在采样时间点之间对过程模型进行多项式拟合。通过采用已知的密集输出法，实践中能够以高精度在位于两个采样时间点之间的任意时刻计算出该过程模型的所求状态量。
采用隐式积分法的一个特别的困难在于，为了计算将来的状态量、即时刻t_k+1的状态量，还需要该物理过程在该时刻出现的状态量x_p，但它很显然并不是已知的，因为时刻t_k+1是未来的时间点。根据本发明，在这种情况下建议：在第二计算单元上，为计算过程模型的将来的状态量x_M，i所需的该物理过程的未来状态量x_p可以通过由该物理过程的已经过去的状态量的外插来得到。
对于第一计算单元上的显式积分法以及第二计算单元上的隐式积分法均延迟结束的情况，根据本发明的另一个有利的实施例建议：为了计算过程模型的状态量x_M，该状态量通过由已经求得的该过程模型的状态量x_M的先前的值的外插来计算。由此保证了在本发明的方法中始终能够提供在所要求的采样时间点的状态量。
根据本发明的另一个具有优点的实施例，过程模型利用第二计算单元和至少一个另外的第三计算单元上的隐式积分法并行处理，并且在第二计算单元和至少一个另外的第三计算单元上并行地共同计算该过程模型的状态量x_M，i。例如在隐式积分法中求解相关联的方程系的任务可以相对较好地并行处理，其中在现有技术中已知有可以自动执行这种并行处理的方法。通过并行处理实现了以下可能性，即在有疑惑的情况下所花的时间比同步处理的显式积分法多得多的隐式积分法能够以很高的精度执行。
此外，开始所提到和所描述的任务通过一种模拟器来实现，尤其是在环模拟器，用于实时计算过程模型的状态量，具有第一计算单元、至少一个第二计算单元和处理接口，其中能够由所述模拟器通过所述处理接口采集物理过程的至少一个状态量，并且可以由所述模拟器输出至少一个状态量，以影响所述物理过程，其中所述模拟器整体上被设置为能够利用它来实现前面所描述的方法。
附图说明
现在逐一给出不同的可能方案，其能够实现并进一步形成本发明所述的方法和本发明所述的模拟器。为此参照从属于权利要求1和3的从属权利要求以及与附图相结合对优选实施例所进行的描述。如图所示：
图1示出与物理过程相关联的过程计算机以及与之相应的模拟器，
图2示出关于第一计算单元和第二计算单元上的显式积分法和隐式积分法的可能的处理时间的示例性示意图，以及所计算的状态量的应用，
图3示出关于第一计算单元和第二计算单元上的显式积分法和隐式积分法的可能的处理时间的另一示例性示意图，以及所计算的状态量的替代应用，并且
图4示出通过外插计算该物理过程的将来的状态量的示意图。
具体实施方式
图1在示意图中示出了具有第一计算单元2和处理接口3的过程计算机1或模拟器1。关于过程计算机1的其它必要细节、例如存储器和将第一计算单元2和存储器相连的数据总线在这里没有示出，其对于本发明来说并不是很重要。
过程计算机1被设置为由过程计算机通过处理接口3采集物理过程4的至少一个状态量x_p，并可以由过程计算机1输出至少一个用于影响物理过程4的输出量y，其中该过程模型由第一计算单元2利用显式积分法进行计算，以求得状态量x_M，e。
如开始处详细描述的，尤其是在计算刚性过程模型的状态量时，显式积分法的应用在实时性条件方面是存在问题的，因为为了确保特定精度，步宽必须选择得很小，从而可能不能再在实时性条件下执行计算。
解决这一问题的第一种方法在图2中示出，其特征在于，提供了第二计算单元5，利用这个第二计算单元5，基本上与利用第一计算单元2对过程模型进行计算同步地，利用隐式积分法执行该过程模型，为了校正在第一计算单元2上所执行的对过程模型的计算，调用利用第二计算单元5和隐式积分法所计算的过程模型的状态量x_M，i，并将以这种方式利用显式积分法求得的状态量x_M，e用作该过程模型的状态量x_M。
在图1中所示的过程计算机中，第一计算单元2和第二计算单元5是一个多内核处理器的两个独立的内核；在另外的在这里没有示出的实施例中，第一计算单元2和第二计算单元5也可以是一个多处理器系统中的两个独立的处理器。第一计算单元2和第二计算单元5对一个共用的存储器进行访问，通过该存储器可以交换数据，例如所计算出的状态量。这在图1中没有详细示出。
在图1中所示的实施例中，这样来实现对利用显式积分法所求得的状态量x_M，e的校正：为了用第一计算单元2计算当前的状态量x_M，e，使用由第二计算单元5所计算的已经过去的状态量x_M，i——如果它存在的话。这种方法在图2中详细示出。时间轴被简要地划分成用k-2、k-1、k、k+1和k+2来表示的积分时间步骤。图中所示的是相应地由第一计算单元2所计算的状态量x_M，e和由第二计算单元所计算的状态量x_M，i。
正如借助所示箭头所看到的，为了在第一计算单元2上利用显式积分法计算将来的状态量，使用了隐式积分法所计算的状态量x_M，i。这只有当在第二计算单元上所执行的隐式积分法没有及时结束时才不能实现，如图2中在时间点k-2和k-1之间的时间间隔就是这种情况。在这种情况下，为了计算状态x_M，e，k-1，利用显式积分法来计算由显式积分法所得到的状态量x_M，e，k-2。原理上，显式积分法总是得到用作该过程模型的状态量x_M的状态量x_M，e。这样，整体上以较大的可靠性保证了计算的实时性。
图3中示出了用于实时计算过程模型的状态量x_M的另一种方法，其首先具有同样的特征，即提供了第二计算单元5。利用这个第二计算单元5，基本上与利用第一计算单元2对过程模型进行计算同步地，利用隐式积分法计算该过程模型，只要利用第二计算单元5和隐式积分法所计算的状态量x_M，i是实时的，就将利用第二计算单元5和隐式积分法所计算的状态量x_M，i用作该过程模型的状态量x_M。否则的话，将利用显式积分法所求得的状态量x_M，e用作该过程模型的状态量x_M。
这两个计算单元中的哪一个为该过程模型提供计算出的状态数据这里仅取决于在第二计算单元5上所执行的隐式积分法是否及时结束。在图3中可以看到，在时间点k-2、k、k+1和k+2之前的计算间隔内，分别由第二计算单元提供通过隐式积分法所计算的状态量x_M，i作为该过程模型的状态量x_M，只有在时间点k-1之前的时间间隔内，才使用由第一计算单元上的显式积分法所求得的状态量x_M，e作为该过程模型的状态量x_M。
在图2和图3所示的实施例中还可以看到，在第一计算单元上的显式积分法用通过隐式积分法所求得的状态量x_M，i来执行，只要它存在的话，即隐式积分法在先前的计算步骤中及时结束。
在图2和图3中示出了在第二计算单元上实现的隐式积分法用和显式积分法同样的计算步宽来执行。在另一种这里没有示出的方法中，隐式积分法的步宽被选为大于显式积分法的步宽。优选的是，计算步宽是可变的，并与状态量随时间的发展相适应，即在瞬时过程的情况下选择得小，而在状态变量只有很小的改变的调谐状态下选择得大。
在导言部分已经说明了隐式积分法的特点是过程模型的状态量的计算在最近的感兴趣的时间步骤内不仅与已经过去的状态量有关，而且还与感兴趣的时间点本身的状态量或其导数有关。因此隐式积分法需要迭代求解。由于这种特性，不仅涉及该状态量本身，还可能涉及相关物理过程的状态量。这意味着，为了计算感兴趣的将来的状态量x_M，i，k+1，还可能需要该物理过程的将来的状态量x_P，k+1，它当然还不存在。
在图4中示出了在这种情况下，为计算所需的物理过程的将来的状态量x_P，k+1，通过外插由该物理过程的已经过去的状态量x_P，k获得；通过外插，分别得到估计值x‾~P,k-2,x‾~P,k-1,x‾~P,k,x‾~P,k+1,x‾~P,k+2,]]>它们分别通过起伏的箭头来表示，用来计算状态量x_M，i。