The present invention relates to a method for using and producing characteristic maps for controlling and regulating various technical devices, in particular in the area of combustion engines, fuel cells and the like.
For the purpose of modeling, calibrating and parameterizing technical devices, characteristic maps are often used, which provide an output variable as a function of input variables. Characteristic maps often do not map, or do not map completely, dependencies to be captured with the aid of physical models.
Such a characteristic map may be read out by a control unit, for example in order to obtain a model parameter, a calibration parameter or a correction parameter as an output variable as a function of operating variables and system parameters as input variables.
Such characteristic maps normally assign an associated output value of the output variable to data points of value combinations of multiple input variables, wherein, for a value combination of input variables that does not correspond to a data point, an output value of the output variables is ascertained by linear or bilinear interpolation. The distribution of the data points is normally defined offline during calibration, i.e., prior to use in the technical device, and therefore cannot be adapted afterwards to a changing behavior during the actual operation of the technical device.
German Patent Application No. DE 10 2010 040 873 A1 describes a method for ascertaining at least one output variable, which depends on a number of input variables, the output variable being described as a function of a first subset of the number of input variables by using at least two assignments, wherein the at least two assignments are formed as a function of respectively one discrete tuple of a second subset of the number of input variables, and the at least one output variable is ascertained in that for present values of the input variables of the second subset, a relationship to at least two of the discrete tuples is ascertained, and in that an interpolation is performed between the output variables of the assignments formed as a function of the at least two discrete tuples using the relationship.
U.S. Patent Application Publication No. US 2011/069336 A1 describes a method, comprising an identification of a target simplex from simplexes having points in a device-independent space, wherein each point comprises a corresponding combination of device-dependent inputs, the identification comprising: determining whether a test simplex contains a target result in the device-independent space, wherein, if the test simplex does not contain the target result, another neighboring simplex is selected as the test simplex and the determination is repeated until the target simplex is identified; and interpolating the device-dependent input of the points of the target simplex in order to identify a combination of device-dependent input for the target result.
The present invention provides a computer-implemented method for providing an output value of an output variable as a function of a value combination of input variables with the aid of a characteristic map, as well as a computer-implemented method for producing a characteristic map.
Additional developments of the present invention are disclosed herein.
According to a first aspect of the present invention, a computer-implemented method for operating a technical device with the aid of a multi-dimensional characteristic map is provided, the characteristic map being defined by data points, to each of which a characteristic map value is assigned. According to an example embodiment of the present invention, for reading out the characteristic map, an output value is determined, as a function of an input variable point to be evaluated for the technical device, with the aid of one-dimensional basis functions, which are assigned to each dimension of a data point, wherein the function values of the one-dimensional basis functions respectively have a monotone curve to a neighboring data point, at which the basis function has the function value 0, and are outside of the region between the data point and the neighboring data point 0, wherein the technical device is operated as a function of the output value.
Furthermore, for an input variable point, the function values of the one-dimensional basis functions of the data points surrounding the input variable point with respect to every dimension are multiplied in order to determine the output value.
Characteristic maps are generally used for the calibration, correction, adaptation and for the modeling of relationships that cannot be physically mapped in a complete manner. A characteristic map assigns a output variable to multiple input variables, which output variable is used in electronic control units of technical devices, in particular combustion engines, fuel cells, autonomous agents and the like.
One feature of the above method according to the present invention is to define the data points of the characteristic map with the aid of basis functions, which allow for a particularly simple production, adaptation and evaluation of the characteristic map. These basis functions may be used regardless of the input dimensionality (number of the mapping input variables of the characteristic map), it being possible to define for each data point of the characteristic map multi-dimensional basis functions as products of one-dimensional basis functions. The data points of the characteristic map correspond in this context to select value combinations of the input variables, to which a specific output value of the characteristic map is in each case directly assigned.
The basis functions are respectively assigned to one dimension of the input variables. The possibility of forming the product of the function values of the basis functions yields a simple interpolation of the output value of the output variable of the characteristic map by products of the one-dimensional basis functions and the determined output values of the output variables at the data points surrounding the queried input variable point.
According to an example embodiment of the present invention, it may be provided that for the computation of the output value for an input variable point having more than two dimensions, multiplication results of function values of the one-dimensional basis functions are stored and used repeatedly.
On account of repeated multiplications, the provision of basis functions as well as the multiplication of the function values of the basis functions for the interpolation of an operating parameter makes it possible to achieve a clear reduction of the number of required multiplications for an interpolation in comparison to conventional interpolation methods.
Furthermore, the data points of the characteristic map may form an unstructured lattice, which comprises basis units as simplexes that connect a number of directly neighboring data points to one another, which is greater by 1 than the dimensionality of the characteristic map, wherein for computing the output value, a transformation of an n-simplex surrounding the input variable point to an n+1 dimensional space is performed as a function of an input variable point and the simplex is transformed to a corresponding unit simplex, wherein the transformation is described by a multiplication with a (n+1)×(n+1) projection matrix, which results from projecting the nodes of the simplex, the output value resulting from the multiplication of the projection matrix with an input variable point complemented by a component having the value 1.
According to one specific embodiment of the present invention, an output value may be extrapolated to an input variable point lying outside of the input variable space in that characteristic map values of multiple edge data points of the characteristic map lying at the edge of the input variable space are summed in a weighted manner, the weighting depending on an angle between the straight line and the path respectively between the edge data points and the input variable point and their distance.
According to a further aspect of the present invention, a system for operating a technical device with the aid of a multi-dimensional characteristic map is provided, the characteristics map being defined by data points, to each of which a characteristic map value is assigned, wherein for reading out the characteristic map, the system is designed to determine an output value as a function of an input variable point to be evaluated for the technical device, with the aid of one-dimensional basis functions, which are assigned to each dimension of a data point, wherein the function values of the one-dimensional basis functions respectively have a monotone curve to a neighboring data point, at which the basis function has the function value 0, and are outside of the region between the data point and the neighboring data point 0, to multiply, for an input variable point, the function values of the one-dimensional basis function of the datapoints surrounding the input variable point with respect to every dimension in order to determine the output value, and to operate the technical device as a function of the output value.
According to a further aspect of the present invention, a computer-implemented method for providing a multi-dimensional characteristic map for operating a technical device is provided, wherein the characteristic map is defined by data points, to each of which a characteristic map value is assigned, an output value being determined as a function of an input variable point to be evaluated for the technical device, with the aid of one-dimensional basis functions, which are assigned to each dimension of an data point, the function values of the one-dimensional basis functions each having a monotone curve to a neighboring data point, which has the function value 0, and being outside of the neighboring data point 0, the characteristic map being calibrated or adapted using one or multiple predefined input variable points and respectively associated output values, in that the characteristic map values at adapted in such a way that the total error between the output values at the input variable points and the output values of the characteristic map is minimized for the input variable points.
It may be provided that the data points of the characteristic map form an unstructured lattice, which comprises basis units as simplexes that connect a number of directly neighboring data points to one another, which is greater by 1 than the dimensionality of the characteristic map, wherein the basis functions of the unstructured lattice are ascertained via the simplexes from selected data points, the density of the distribution of the data points being selected in such a way that the expected behavior of the output value can be mapped by linear interpolation between the data points.
According to a further aspect of the present invention, a system for providing a multi-dimensional characteristic map for operating a technical device is provided, wherein the characteristics map is defined by data points, to each of which a characteristic map value is assigned, an output value being determined as a function of an input variable point to be evaluated for the technical device, with the aid of one-dimensional basis functions, which are assigned to each dimension of an data point, the function values of the one-dimensional basis functions each having a monotone curve to a neighboring data point, at which the basis function has the function value 0, and being outside of the region between the data point and the neighboring data point 0, the system being designed to calibrate or to adapt the characteristic map using one or multiple predefined input variable points and respectively associated output values, in that the characteristic map values at adapted in such a way that the total error between the output values at the input variable points and the output values of the characteristic map is minimized for the input variable points.
Specific example embodiments of the present invention are explained in greater detail below on the basis of the figures.
For operating the technical device 2, control unit 3 provides for ascertaining an operating parameter B, which may represent a correction parameter, an adaptation parameter or a function value of a function mapping a physical behavior. For ascertaining the operating parameter B, the control unit 2 uses the characteristic map in the characteristic map memory 4 and operates the technical device 3 in accordance with the ascertained operating parameter B.
For each input variable point, a multi-dimensional basis function is defined, which is a product of the individual basis functions. An output value of the output variable may thus be computed from a characteristic map as:
where the index i takes into account each of the data points of the characteristic map lattice.
The basis functions bi are computed as products of the one-dimensional basis functions at the input value of the respective dimension of the input variable of the characteristic map.
For an individual dimension x, the basis functions, as shown in
Accordingly, the multi-dimensional basis function is then ascertained by multiplication
b
i()=bix
For training such a characteristic map, an output value of the output variable y=f′(x) is assigned to a data point. For this purpose, a learning algorithm receives an operating parameter to be learned at a specific data point x1, x2, . . . , whereby the latter can be used to improve or enter the existing learned values.
After a sufficient number of learning events, the characteristic map may indicate a correct output value of an output variable in accordance with a predefined input variable point (input variable vector). If the characteristic map is to exhibit a PT1 behavior, the output value output by the characteristic map will tend in the direction of the actual operating parameter to be learned, according to:
f′({right arrow over (x)})→({right arrow over (x)})
If an integrating behavior is to be stored, a discrete integral of the input variable point {right arrow over (x)} results as the output value of the characteristic map
where K is an integration speed parameter and τ corresponds to the preceding discrete time steps. No continuous function is available, however, as the output f′ of the characteristic map; rather, the output values for corresponding input variable points must be approximated on the basis of the characteristic map values at the data points (lattice intersections of the characteristic map or entries at the data points). Hence
where bi({right arrow over (x)}) are basis functions and yi the corresponding discrete learned characteristic map values at the data points of the characteristic map.
During an online learning step, a measurement f({right arrow over (x)}) is evaluated. First, a residual error δ is computed, which represents the error of the currently learned value. An integrator behavior corresponds to δ=f({right arrow over (x)}). For a PT1 behavior, δ=f({right arrow over (x)})−f′({right arrow over (x)}), which corresponds to the difference between the characteristic map value of the characteristic map and the output value currently to be learned at the input variable point of the measurement.
Next, the learned characteristic map values yi at the data points are modified in such a way that f′({right arrow over (x)}) corresponds better to the correct output values defined above, i.e. the residual error is compensated. This is achieve in that the basis functions are used as weights for modifying the learned characteristic map values
y
i
→y
i
+Kb
i′({right arrow over (x)})δ
where K represents a learning speed and to which K in
may correspond as the integration speed parameter.
During the offline learning, the learned characteristic map values yi are determined in such a way that the outputs f′({right arrow over (x)}) agree best with the output value of the characteristic map for the input variable point (evaluation point) .
This may be carried out by the method of the smallest squares in accordance with
where the matrix elements are given by
X
ji
=b
i({right arrow over (x)}j)
The sum operation of the equation
is performed in every line in the product x·{right arrow over (y)}.
A learned characteristic map value yi thus exists for every basis function bi(). These basis functions bi(
) are selected in order to span a multi-dimensional volume Ω, in which a learning operation is to be performed.
As seen in
The lattice points are indicated by all combinations of the points {x1}, {x2}, i.e. all gray circles in
A multidimensional basis function bi is defined for every lattice point {right arrow over (x)}1. The basis functions bi are computed as products of the one-dimensional basis functions according to every dimension of the input variables of the characteristic map. For an individual dimension x, the basis functions, as shown in
b
i()=bix
The properties given in the above definition of the basis functions may then be expanded to the higher dimensionalities
At a certain input variable point {right arrow over (x)} the basis functions, which correspond to 2N multi-dimensional data points, are unequal to 0, where N represents the number of dimensions. Thus, 2N learned values are accessed for the interpolation or modified by a learning step. The multi-dimensional data points {right arrow over (x)}1 comprise products of the one-dimensional basis functions. For every dimension, the one-dimensional basis functions of a low (index l) and an upper (index u) data point are taken into account, which enclose the input variable point {right arrow over (x)} to be evaluated. For example, in the case of three dimensions, the eight (23) multi-dimensional basis functions correspond to the eight corners of the cuboid, which enclose the input variable point {right arrow over (x)}:
where the index “l” corresponds to the lower and the index “u” corresponds to the upper data point. Since products are formed multiple times in the computation of the multi-dimensional basis functions, it is possible to use a computation tree-based approach, as illustrated in
An extrapolation of the output value on input variable points to be evaluated that lie outside of the input variable space Ω is performed by projecting the input variable point onto a limit of the input variable space Ω. Since the input variable space Ω is always convex, this projection is unambiguous.
In contrast to the above-described exemplary embodiment, characteristic maps may also be unstructured, i.e. have no hypercuboid contour. This may be expedient if the value of the input variable point (evaluation point) to be learned exists only for a non-cuboid set of data points of the input variable space. In a cubically arranged lattice, it may otherwise happen that the output values to be learned for the input variable points are not distributed over the entire input variable space and that thus some output value are never updated or accessed. This results on the one hand in a waste of resources, since the unused learned operating parameters must be stored, and on the other hand the learned values are not extrapolated in these regions during the readout, since they are not located in the extrapolation region, i.e. are not located outside of the input variable space Ω. Instead, a learned setpoint value, such as e.g. zero, is output in these regions, exactly as in the corresponding extrapolation regions.
In addition, the resolution of the learned values cannot be selected as desired with the aid of the above-described routines. With the aid of rectangular lattices, the data points can only be refined in a dimension-wise manner. The refinement in one dimension will thus be applied to all combinations of the other dimensions, regardless of whether this is necessary or not. This results in a waste of resources, since high resolutions are unnecessarily introduced in operating areas where these are not necessary. Unnecessarily high resolutions may also result in lower performance and noise suppression, since measuring noise is falsely interpreted as spatial variation.
An approach for applying unstructured characteristic maps to the learning algorithm described above is described below. The data point lattices of the characteristic maps may be selected to describe arbitrary forms and resolutions with the aid of simplexes, i.e. 1-D line segments, 2-D triangles, 3-D tetrahedrons etc., as basis units. The approach may be applied to any desired number of dimensions. In the above-described learning and evaluation approach for a cuboid data point distribution, an input variable space Ω may be spanned by the data points of the characteristic map {right arrow over (x)}i. For each data point of the characteristic map, a value yi to be learned is stored. Learning and reading out are performed with the aid of basis functions bi(). The basis functions bi(
) are defined as indicated above.
Above, the data points were defined on a rectangular characteristic map lattice, which is defined by the individual data points for every dimension. For the application of the approach described above, the data points of unstructured characteristic maps are spanned by independent data points, as shown in
The computation of the linear basis functions of unstructured characteristic map lattices may be performed efficiently with the aid of barycentric coordinates. For this purpose, a transformation of an n-simplex into an n+1-dimensional space is performed, and the simplex is transformed onto a corresponding unit simplex. For example, a 2-D triangle may be transformed onto the unit 2-simplex in three dimensions, as shown in
{right arrow over (λ)}k=Pk·{right arrow over (x)}′
Here {right arrow over (x)}′ corresponds to a (n+1)-dimensional vector as a function of the n-dimensional vector, {right arrow over (x)}, to which a component with the value 1 is appended, e.g. (x1, x2, 1). The values of Pk are obtained by projecting the nodes of the simplex e.g. for Ω1 in
P
1
·{right arrow over (x)}′
1=(1,0,0),
P
1
·{right arrow over (x)}′
2=(0,1,0),
P
1
·{right arrow over (x)}′
3=(0,0,1),
i.e., the columns of the inverse matrix P−11 correspond to the coordinates of the nodes of the simplex, to which 1 is appended.
The barycentric coordinates have the following advantages:
The basis functions in unstructured lattices may be ascertained via the simplexes from the selected data points. The data points are selected in such a way that firstly they cover the expected range of the input variable point and that secondly the density of their distribution is sufficiently high that the expected behavior of the output value may be mapped by linear interpolation between the data points.
An extrapolation from the unstructured characteristic map lattices cannot be performed in a simple manner, as described above, because the characteristic map lattice is not necessarily convex and therefore an unambiguous projection onto the limit does not always exist. Accordingly, for unstructured characteristic map lattices, it is proposed to carry out the following method in order to obtain continuous values for data points outside of the discretized input variable space Ω. This makes it possible further to avoid jumps in the output values for continuously changing input variable points.
The oriented edges Lk form the limit of the input variable space Ω, with the outwardly directed normals {right arrow over (n)}k, as illustrated in
({right arrow over (x)}−{right arrow over (x)}k)·{right arrow over (n)}k>0
where {right arrow over (x)}k is a point on the edge Lk, e.g. one of the limit nodes. For each of these edges, the edge point {right arrow over (x)}near on the edge Lk is determined, which is closest to the input variable point {right arrow over (x)} to be evaluated. This point may be on the edge or on a limit node of the edge. The corresponding output value for the extrapolation is the interpolated value at position {right arrow over (x)}near, whereby a weighting, given by
is taken into account. Here, d is the Euclidean distance between {right arrow over (x)} and the edge point {right arrow over (x)}near, and δ is the angle between the normal {right arrow over (n)} and ({right arrow over (x)}−{right arrow over (x)}near) The extrapolated output value y′ may then be computed as
Number | Date | Country | Kind |
---|---|---|---|
10 2020 208 321.5 | Jul 2020 | DE | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2021/065755 | 6/11/2021 | WO |