1. Field of the Invention
The present invention relates to control units for engine systems, in particular control units having a separate arithmetic unit for evaluating data-based function models, for example, Gaussian process models.
2. Description of the Related Art
Up to now, function models in control units, i.e., the control path and system models, have been implemented through the specification of characteristic maps, characteristic curves, or functions emulating a physical system. These models are adapted by a user by adjusting model parameters to the conditions of the physical system.
The use of non-parametric data-based function models constitutes an alternative, with which the functions of physical systems may be essentially emulated without parameter specifications. For example, a Gaussian process model, which is essentially defined using hyperparameters and nodes, may be used as a data-based function model. The data-based function models are created based on training data which may be ascertained in a testing system. The nodes for the Gaussian process model may correspond to the training data, may be selected from these data, or may be generated from the training data.
In particular, local effects may not be properly mapped by the created data-based function model under certain circumstances. If a data-based function model has already been determined based on training data of an initial training data record, it is difficult to take into consideration training data of a subsequently ascertained training data record in a proper manner in the data-based function model which has already been created. However, simply merging the training data records with or without variation of the hyperparameters of the data-based model allows local effects to be properly taken into consideration only if the training data of the subsequently added training data record have a sufficient number of measuring points. Thus, the additional measuring points may have a sufficient weight in relation to the measuring points of the training data of the initial training data record. Furthermore, it is required that the measuring points of the subsequently added training data record do not conflict with already existing measuring points of the initial training data record, i.e., have a relatively large deviation from them. Otherwise, data-based function models are obtained having high measurement noise and an accordingly high modeling error for the function values in the range of the local effect.
Control units having a microcontroller and a separate model calculation unit for calculating data-based models in a control unit are known from the related art. Thus, for example, a control unit having an additional logic circuit is known from Published German patent application document DE 10 2010 028 259 A1, which is designed for calculating exponential functions in order to support the execution of Bayesian regression methods, which are required in particular for the calculation of Gaussian process models.
Furthermore, an additional method for adding measuring points of an additional training data record to an existing Gaussian process model is known from C. Plagemann, K. Kersting, W. Burgard, “Nonstationary Gaussian Process Regression Using Point Estimates of Local Smoothness,” ICML Proceedings, pp. 204-2116, 2006. However, this method is inefficient, in particular since the parameter optimization is difficult.
According to a first aspect of the present invention, a method for creating a data-based function model is provided. The method includes the following steps:
One idea of the aforementioned method for creating a data-based function model is to additively interconnect multiple data-based partial models in order to implement a control path or system model in a control unit. This enables a proper representation of local effects, which are modeled as a separate data-based partial model. The aforementioned method provides for forming an additive chain of a number of data-based partial models in which, based on a first data-based partial model, each additional data-based partial model specifies an additive deviation from the first data-based partial model or the sum of the first data-based partial model and the additional data-based partial models already taken into consideration. In this way, local effects may be well modeled, in particular based on the breadth of the data-based model determined by the hyperparameters, without substantially impairing the first data-based partial model.
The additive merging of multiple data-based partial models has the additional advantage that they may be calculated in a simple manner in a control unit using an additional model calculation unit.
By providing the model calculation unit, which calculates multiple data-based partial models in succession in an independent manner and adds up the resulting partial results, a simple ascertainment of a function value based on the aforementioned additive function model is made possible. By avoiding the arithmetic unit having to prepare the model calculation unit for the calculation of each individual partial model, the calculation of a corresponding function value may be considerably accelerated.
Furthermore, it may be provided that each of the additional training data records contains training data which are associated with a local effect using a classification or clustering method.
In particular, the additional data-based partial model may be ascertained from the difference training data record in such a way that it drops to a constant value of zero in the extrapolation range.
Furthermore, the data-based partial models may correspond to Gaussian process models, the one or the multiple additional data-based partial models being ascertained from the corresponding difference training data record based on an average value function which is constantly zero.
According to one additional aspect, a method for calculating a function value of a data-based function model is provided in a control unit which has an arithmetic unit and a separate model calculation unit, the data-based function model having multiple partial models. The method includes the following steps:
The combination of the first-mentioned method for creating a data-based function model in connection with the calculation in a model calculation unit with the aid of a DMA unit enables the execution of a linked calculation method in the model calculation unit, which in turn the calculations to be carried out rapidly without loading the arithmetic unit.
Furthermore, the additional configuration data may be retrieved from a memory area which is a function of the calculation address, in particular by applying a predefined address offset to the calculation address.
According to one specific embodiment, the stop condition may be specified by the retrieved additional configuration data.
It may be provided that the stop condition specifies that the arithmetic unit identifies the value stored in the accumulator as a function value of a data-based function model.
In particular, the retrieval of the additional configuration data from the memory unit, the transmission of the relevant additional configuration data to the model calculation unit, and the starting of the ascertainment of the additional partial function value of the additional partial model may be carried out by the DMA unit.
According to one additional aspect, a device, in particular an arithmetic unit, is provided for creating a data-based function model, the device being designed to:
According to one additional aspect, a control unit for calculating a function value of a data-based function model is provided, the data-based function model having multiple partial models, including:
According to one additional aspect, a computer program is provided which is configured to carry out all steps of the aforementioned method for creating a data-based function model.
A method is initially described based on the flow chart of
In step S1, a first data-based partial model ffirst
In step S2, a second training data record
(xi, yi), i=1, . . . N
is furthermore provided, where xi represents the p-dimensional measuring points and yi represents the scalar output values.
In step S3, the deviations {tilde over (y)}i between the model predictions ffirst
{tilde over (y)}
i
=y
i
−f
first
partial
model(xi)
In step S4, an additional, second data-based model fsecond
In step S5, the first, initially provided data-based partial model and the second, additional data-based partial model may be additively connected:
f(x)=ffirst
where f(x) corresponds to the function value of the overall data-based function model.
This concept may be generalized in a simple manner by allowing any number of additional partial models. For example, measuring points of an additional training data record, which, for example, were associated with a certain local effect via a classification or clustering method, may be modeled in an additional data-based partial model. The additional data-based partial model is created from a difference training data record whose output values correspond to the differences between the output values of the additional training data record to be taken into consideration and the corresponding function values of a data-based model formed from an addition of the previously ascertained data-based partial models at the measuring points of the additional training data record. The additional data-based partial model is then ascertained based on the difference training data record. This approach may be repeated for all additional training data records.
For a quantity I of additional data-based partial models which were ascertained in the above-described manner, the following applies:
For each individual result of a partial model, a model-dependent scaling factor an (where n=0, . . . , I) exists, with which the weight of the relevant partial model
may be taken into consideration. In addition, a global model offset also exists, which is defined as c0 and may be configured as a part of hyperparameters. These parameters belong to the model parameters and describe the output normalization of the function model to be ascertained.
Control unit 1 receives external input signals, for example, from sensors or other control units such as other modules of engine system 1, via input interface 11. The received signals include input values and may also indicate user specifications in addition to sensor signals.
Control unit 1 is used to carry out the control unit functions in order to calculate one or multiple output values as a function of the input values which are received via input interface 11, and to output them via output interface 13. Control unit 1 must calculate control path and system models for carrying out the functions, and uses a separate model calculation unit 15 for this which is designed for calculating data-based function models. In particular, model calculation unit 15 is specialized in the calculation of exponential functions, which are required for calculating Gaussian process models.
Data are stored in memory unit 12 which were determined for certain output values offline, i.e., in test measurements, and, for example, include hyperparameters of one or multiple data-based function models and nodes (node data), which contain a part of or all training data for the relevant data-based model.
In addition, parameters and variables which were received or calculated by arithmetic unit 10 may also be stored in memory unit 12.
Furthermore, control unit 1 includes DMA unit 14, which makes it possible for model calculation unit 15 to access memory unit 12 directly, in addition to arithmetic unit 10. Furthermore, DMA unit 14 may start a calculation of a data-based model by retrieving corresponding parameters and node data from memory unit 12 and forwarding them to model calculation unit 15.
For calculating data-based models, it may be provided that arithmetic unit 10 transmits the hyperparameters and node data (training data) directly to model calculation unit 15, or communicates, in the form of a pointer, to a memory area of memory unit 12 of model calculation unit 15, from which area of memory unit 12 the hyperparameters required for calculating the data-based model and node data (training data) may be retrieved as configuration data. Model calculation unit 15 may then access memory unit 12 and retrieve the parameters and node data (training data) with the aid of DMA unit 14. Furthermore, DMA unit 14 may be designed to automatically read in a next record of hyperparameters and node data as configuration data into the model calculation unit, after a function value of a data-based partial model has been calculated. The consecutive calculations are terminated if a stop condition is fulfilled.
As schematically depicted in
First configuration folder 21 represents a memory area in memory unit 12 in which the parameters and node data for the calculation of the first data-based partial model are stored in configuration registers 22.
After the start prompt, DMA unit 14 retrieves the data from first configuration folder 21 in step S11 and transmits these data in step S12 to model calculation unit 15.
After the hyperparameters and node data for the first partial model have been transferred into model calculation unit 15, an accumulator 31 of model calculation unit 15 is reset to the aforementioned model offset value c0, which may be received as a part of the hyperparameters, and in step S13, a first partial function value of the first data-based partial model is calculated, and the first partial function value is added to the value stored in accumulator 31 of model calculation unit 15.
In step S14, the completion of the execution of the calculation of the first partial function value is communicated to DMA unit 14.
Thereupon, in step S15, DMA unit 14 requests the parameters and node data for a second partial model of a start address of a second configuration folder 21 from memory unit 14. The start address of second configuration folder 21 is a function of calculation data address BA and is ascertained through addition of a predefined address offset O, which indicates the memory size of configuration folder 21. The parameters and node data for the second partial model are transferred to model calculation unit 15 in order to calculate the second partial model.
The second partial model is calculated in step S16 and the resulting second partial function value is subsequently added to the value stored in the accumulator. If a corresponding parameter has been stored in the relevant configuration folder 21, the second partial function value may be weighted beforehand with the aid of a weighting factor ai, in particular applied multiplicatively.
In step S17, the completion of the execution of the calculation of the second partial function value is communicated to DMA unit 14.
Steps S15 through S17 are now repeated as often as corresponds to the number of configuration folders 21 which are stored in memory unit 12 for the calculation of a corresponding partial function value.
Configuration folder 21, which is used for calculating the last partial result, may be identified accordingly. Alternatively, configuration folder 21 following the last calculation may be correspondingly configured as a stop command. In this way, after the calculation of the nth partial result and the additive application to the value stored in the accumulator, in step S18, a calculation end-interrupt may be transferred directly or indirectly via DMA unit 14 to arithmetic unit 10, so that arithmetic unit 10 retrieves the final result from the accumulator as a function value of the model calculation by the model calculation unit 15.
The last configuration register 22 of each configuration folder 21 may predefine a corresponding configuration for the start of the calculation in model calculation unit 15. In this way, the configuration process of model calculation unit 15 and the start of the calculation process may be managed and initiated via DMA unit 14.
After the calculation of the last partial function value and after the ascertainment of the overall function value in the accumulator, in order for model calculation unit 15 to provide this value for immediate retrieval by arithmetic unit 10, configuration folder 21 of the nth partial model may set an interrupt setting on arithmetic unit 10. Alternatively, the interrupt setting may also be set on DMA unit 14, which in turn forwards the interrupt to arithmetic unit 10.
During the transmission of the parameters and node data to model calculation unit 15, a DMA function “DMA—Repeat—Copy” known per se for DMA units may be used, which may be triggered for the first partial model by arithmetic unit 10 and subsequently by the calculation end-interrupt of model calculation unit 15.
Number | Date | Country | Kind |
---|---|---|---|
10 2013 206 292.3 | Apr 2013 | DE | national |