EP2791827A2 - Method for evaluating the solution to a multicriteria optimization problem - Google Patents

Method for evaluating the solution to a multicriteria optimization problem

Info

Publication number
EP2791827A2
EP2791827A2 EP12795756.1A EP12795756A EP2791827A2 EP 2791827 A2 EP2791827 A2 EP 2791827A2 EP 12795756 A EP12795756 A EP 12795756A EP 2791827 A2 EP2791827 A2 EP 2791827A2
Authority
EP
European Patent Office
Prior art keywords
space
solution
variation
model
variables
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.)
Ceased
Application number
EP12795756.1A
Other languages
German (de)
French (fr)
Inventor
Klemens WALLNER
Alejandra GARCIA
Adnand DRAGOTI
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
AVL List GmbH
Original Assignee
AVL List GmbH
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by AVL List GmbH filed Critical AVL List GmbH
Publication of EP2791827A2 publication Critical patent/EP2791827A2/en
Ceased legal-status Critical Current

Links

Classifications

    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02D—CONTROLLING COMBUSTION ENGINES
    • F02D41/00—Electrical control of supply of combustible mixture or its constituents
    • F02D41/02—Circuit arrangements for generating control signals
    • F02D41/14—Introducing closed-loop corrections
    • F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
    • F02D41/1406—Introducing closed-loop corrections characterised by the control or regulation method with use of a optimisation method, e.g. iteration
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/0205—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
    • G05B13/024—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system in which a parameter or coefficient is automatically adjusted to optimise the performance
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/0265—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/20—Drawing from basic elements, e.g. lines or circles
    • G06T11/206—Drawing of charts or graphs
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design

Definitions

  • the present invention relates to a method for evaluating the solution of a multi-criteria optimization problem.
  • Such multicriteric optimization problems generally exist if target values of several objective functions are to be optimized simultaneously depending on several restrictions (such as boundary conditions, physical limits, etc.). However, if the objective functions are in conflict with each other, the simultaneous optimization of all objective functions is often a problem.
  • multicriteria optimization problems usually can not find unique solutions but only a set of possible multidimensional solution points in a multi-dimensional space a multi-dimensional surface, the so-called Pareto front, which all represent an optimal compromise of the multi-criteria optimization problem. Individual points of this Pareto front therefore represent different, but in each case optimal, compromises between the objective functions.
  • Such multi-criteria optimization problems are known per se and there are also a number of mathematical methods for solving such problems.
  • a method is known with which an ECU is optimized during operation by means of a multi-criteria optimization problem with respect to an exhaust gas soot consumption compromise.
  • a single aggregate objective function (AOF) is used, which combines the weighted objective functions into a functional one common solution is the linearly weighted summation of the objective functions described here.
  • Each objective function is assigned a weighting factor. tor, from which a scalar objective function is derived as the sum.
  • the actual optimization can be carried out using conventional approaches, for example by means of sequential quadratic programming (SQP), an effective, iterative method for nonlinear limited optimization, which is necessary for the desired reduction of the required computing power.
  • SQL sequential quadratic programming
  • US Pat. No. 7,921,371 B1 describes a method for visualizing the solution of a multi-dimensional multi-criteria optimization problem. All target values are displayed on parallel, adjacent axes, whereby also the smallest and largest value of the found optimal solution is visualized, which represent the entire possible range of the optimal solutions for a target size. Therefore, the pareto front is represented in the form of parallel axes for the number of target variables. For each target size, there is a fixed target value, also shown in the diagram, with the parallel axes shifted vertically so that the target values all lie on a horizontal line.
  • WO 01/67395 A1 discloses a method in which, in a representation of all possible optimal solutions, those are marked (here by a different color), which satisfy a certain user-specified criterion, such as a constraint.
  • the solutions are projected onto two or three-dimensional surfaces in order to be able to recognize connections.
  • This object is achieved according to the invention by representing the set of optimal solutions of the multicriteria optimization problem in a model space as a two- or three-dimensional diagram of the objective functions and simultaneously displaying at least one of the objective functions as a function of at least one variation variable in a variation space
  • Model space and the variation space are interactively interconnected by marked for each selected solution in the model space, the solution size underlying the solution in the variation space.
  • a mathematical model can be used as the objective function, which is determined from a number of measurements of the objective function as a function of the variables. This makes it possible to apply the inventive method practically to any optimization problems.
  • MOP Multi-Criteria Optimization Problems
  • the variation quantities x are e.g. during a calibration, the settings, e.g. at the test bench where the calibrator carries out his measurements. Since these variation quantities x represent valid points in the variation space, the calibrator knows that he can adjust the variation quantities x in this range. For this reason, in this variation space, around the variation quantities x, an envelope is laid, which is referred to as design space.
  • the design room thus contains all the variables x valid for the respective application.
  • evolutionary algorithms for multicriteria optimization do not require any weighting or a priori information, so that these methods have been increasingly used in recent years and in particular have proven to be effective and robust methods.
  • genetic algorithms - based on selection, recombination and mutation - were used to achieve a continuous approach to a desired target. They are easy to apply to a wide range of problems and are very robust in the search for global optima, even if there are a variety of local optima.
  • conflicting requirements are searched for in a set of compromise solutions that approximate the best possible solutions.
  • the quality of an approximation can be quantified by the volume dominated by it in the target space, the S metric. Maximizing the S-metric is a desirable goal and, at the same time, sufficient scalar replacement of the original objective function.
  • a genetic algorithm uses these within the selection and thus achieves excellent results, especially when more than three goals have to be optimized where other multicriteria genetic algorithms fail.
  • the currently most popular genetic algorithm for determining the Pareto front is the Non-Dominated Sorting Genetic Algorithm (NSGA-II), which has proven to be extremely efficient, especially for the determination of a global optimum.
  • NSGA-II is a high-performance, elitist algorithm that prioritizes non-dominating individuals and preserves the diversity of solutions. The algorithm creates an initial population within the vector of variational magnitudes x and approaches iteratively in an iterative process that is based on selection,
  • the peculiarity of the visual presentation lies in a split representation of model space 1 and variation space 2, as shown in FIG.
  • f 2 (x) "NOX” and f 3 (x) "Fuel Consumption” The illustration of the design space game of FIG. Three-dimensional space, the fi through the objective functions (x) "Smoke”, is spanned. the valid range 3 in which the solutions can move is contained in this model space 1.
  • the Pareto front 4 contains the solutions found for the multi-criteria optimization problem within this valid range 3.
  • the Model space 1 can also be represented by several two- or three-dimensional representations of the k dimensions of the model space. These dimensions of the k dimensions in the two- or three-dimensional representations can be made dependent on the multi-criteria optimization problem and the preference of the user is represented by a number of two or three-dimensional representations of objective functions f j (x) and variation quantities x In Place.
  • the objective function is f 3 (x) each represented as a function of the three varying sizes x- ⁇ "exhaust gas temperature" x 2 "EGR rate” and x 3 "rail pressure".
  • any combinations of objective functions f j (x) and variation quantities x are conceivable here.
  • the objective functions f j (x) can be known functions of the variables x.
  • an objective function f j (x) is a mathematical model that is determined from measurements or experiments. In this case, measurements are carried out on the object of the multi-criteria optimization problem, eg on an internal combustion engine, a drive train, a transmission, a vehicle, etc., for example on corresponding test benches or in the course of test drives.
  • the desired objective functions f j are measured as a function of the variables x and possibly other variables. Mathematical models of the objective functions f j are then created from these measured variables.
  • Mathematical models of the objective functions f j are then created from these measured variables.
  • Possible models include a polynomial regression model, a fast neural network or an intelligent neural network. Because of this approach, additional measurements, ie real measurements, will not necessarily be 100% accurate on this model.
  • the methods for determining the models therefore also provide a model confidence interval indicating the bandwidth in which further measurements are likely to move. This means that a model with a slim model confidence range fits relatively well on the measurements made and has accordingly good explanatory power. The closer the model confidence interval, the better the models fit the measurements, and the more likely the model-determined solutions to the multi-criteria optimization problem will actually yield the sought-after values.
  • well-known objective functions ( ⁇ ) can also have a model confidence range, which in turn indicates how exactly an objective function f j agrees with a real measurement.
  • the model confidence range is thus a measure of the accuracy of the model or a target function with respect to real measurements.
  • the model confidence ranges 5 can also be shown in the individual diagrams, for example in the form of an upper and lower limit as shown in FIG.
  • the special feature of this type of representation is that the objective functions ( ⁇ ) in the model space 1 and the Pareto front 4 are represented as a set of possible optimal solutions of the multi-criteria optimization problem together with the variables x and therefore can also be analyzed together.
  • the representation in the variation space 2 is interactively adapted to a selection of a point in the model space 1.
  • a crosshair 6 is provided. With the example, an interesting point 7 of the Pareto front 4, or the valid area 3, is selected.
  • the crosshair 6 automatically marks the variation quantities x for this punk 7 in the model space 1.
  • the respective values of the variables x at this point can also be indicated, as indicated in FIG. 1.
  • the model trust area 5 can be represented so that the user additionally receives information about how trustworthy the underlying objective function fj (or mathematical model) is on this point.
  • this type of representation also allows the analysis of the effects of changes in the specification of the constraints g, and / or the range x min , x ma x of the variables x.
  • changes lead to other solutions that can then be easily compared directly.
  • it can be provided, for example, to change the limits of the range of variation variables x, for example by means of a slider in the variation space 2, whereby the representation of the solution in the model space 1 can change at the same time.
  • a number of measurements are performed on the internal combustion engine, wherein the target variables of the target functions f j (x) NOx, soot and consumption depending on the variables x, eg exhaust gas temperature, EGR rate, rail pressure, are measured.
  • the number and the amount of measurements can be predetermined, eg by a given Design of Experiment.
  • the measurements are used to determine mathematical models and model confidence regions 5 for the objective functions f j (x).
  • the multicriteria optimization problem for optimizing the objective functions f j (x) can be solved and the solution in the split representation of model space 1 and variation space 2 can be analyzed.
  • the calibrator can examine various optimal solutions of the pareto-front 4 with regard to the underlying variation variables x and the model confidence region 5. From these possible optimal solutions, the calibrator then determines one of the solutions as the best possible compromise.
  • the experience of the calibrator plays a major role.
  • the model confidence ranges and the Dependencies of the variables x are also the values of additional model channels that have not been optimized as target functions, as well as the robustness of the settings, eg whether the models in the vicinity of the optimum changes greatly, low influenceability due to component tolerances, etc., are taken into account.
  • This can be repeated for all operating points required for the calibration (eg speed, torque, load) of the internal combustion engine.
  • a predetermined number of operating points eg ten to twenty operating points, are generally required.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Automation & Control Theory (AREA)
  • Medical Informatics (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Image Processing (AREA)

Abstract

Solving a multidimensional multicriteria optimization problem is difficult because the correlations and dependencies between solutions, target functions and variation variables can be detected only with difficulty. In order to facilitate this, it is proposed that a model space (1) and a variation space (2) be displayed simultaneously and in interactively linked fashion.

Description

Verfahren zur Auswertung der Lösung eines multikriteriellen Optimierungsproblems  Method for evaluating the solution of a multi-criteria optimization problem
Die gegenständliche Erfindung betrifft ein Verfahren zur Auswertung der Lösung eines multikriteriellen Optimierungsproblems. The present invention relates to a method for evaluating the solution of a multi-criteria optimization problem.
Bei der Kalibrierung von Fahrzeugsteuergeräten (xCU), wie z.B. ein Motorsteuergerät ECU oder ein Getriebesteuergerät TCU, auf Fahrzeug- und Komponentenprüfständen, aber auch bei Erprobungsfahrten, wird der Kalibrieringenieur regelmäßig vor die Aufgabe gestellt, einen Kompromiss zwischen unterschiedlichen Zielgrößen zu optimieren. Probleme, bei denen Zielgrößen mehrerer Zielfunktionen gleichzeitig in Abhängigkeit mehrerer Einschränkungen zu optimieren sind, sind in der Praxis und insbesondere in der Kalibrierung weit verbreitet. Beispiele dafür sind z.B. der NOx-Ruß Kompromiss, oder der NOx-Ruß-Verbrauch Kompromiss bei klassischen verbrennungsmotorischen Antrieben, oder der Batterieschädigung- Verbrauch Kompromiss bei Hybridanwendungen, oder der Kompromiss zwischen sportlicher aber trotzdem komfortabler Abstimmung der Schaltvorgänge bei Getrieben. Generell ergeben sich bei Fahrzeugen durch Antriebskonzepte mit Automatikgetrieben und/oder Hybridi- sierung mehr und mehr widersprüchliche Anforderungen an die einzelnen Komponenten, die in Form eines Kompromiss dargestellt und optimiert werden müssen. Solche multikriteriellen Optimierungsprobleme bestehen allgemein, wenn Zielgrößen mehrerer Zielfunktionen gleichzeitig in Abhängigkeit mehrerer Einschränkungen (wie Randbedingungen, physikalische Grenzen, etc.) zu optimieren sind. Wenn die Zielfunktionen in einem Zielkonflikt zueinander stehen stellt sich das gleichzeitige Optimieren aller Zielfunktionen aber häufig als Problem dar. Insbesondere können bei solchen multikriteriellen Optimierungsproblemen in der Regel keine eindeutigen Lösungen gefunden werden, sondern nur eine Menge von möglichen mehrdimensionalen Lösungspunkten in einem mehrdimensionalen Raum, also eine mehrdimensionale Fläche, die sogenannte Pareto-Front, die alle einen optimalen Kompromiss des multikriteriellen Optimierungsproblems darstellen. Einzelne Punkte dieser Pareto-Front stellen daher unterschiedliche, aber jeweils optimale Kompromisse zwischen den Zielfunktionen dar. Solche multikriteriellen Optimierungsprobleme sind an sich bekannt und es gibt auch eine Reihe von mathematischen Methoden zur Lösung solcher Probleme. Aus der EP 2 192 294 A1 ist z.B. ein Verfahren bekannt, mit dem eine ECU im laufenden Betrieb mittels eines multikriteriellen Optimierungsproblemes hinsichtlich eines Abgas-Ruß- Verbrauch Kompromisses optimiert wird. Dabei wird eine einzelne Aggregate Objective Function (AOF) verwendet, die die gewichteten Zielfunktionen in einem Funktional zusammenfasse Eine gebräuchliche Lösung stellt die hier beschriebene linear gewichtete Auf- summierung der Zielfunktionen dar. Dabei wird jede Zielfunktion mit einem Gewichtungsfak- tor versehen, woraus sich als Summe eine skalare Zielfunktion ableitet. Hierbei kann die eigentliche Optimierung mit herkömmlichen Ansätzen durchgeführt werden, z.B. mittels Sequential Quadratic Programming (SQP), einer effektiven, iterativen Methode für nichtlineare beschränkte Optimierung, was für die gewünschte Reduktion der benötigten Rechen- leistung notwendig ist. Allerdings hängt die Aussagekraft einer solchen Optimierung in starkem Maße von der Wahl der Gewichtungsfaktoren ab, die aber in vielen Fällen nicht zuverlässig im Vorfeld festgelegt werden können, womit die Ergebnisse einer solchen Optimierung nicht immer zufriedenstellend sind bzw. diese Methode überhaupt nur auf eine begrenzte Gruppe vom Optimierungsproblemen anwendbar ist. Ein weiteres Problem bei multikriteriellen Optimierungsproblemen ist, die gefundene Lösung zu visualisieren und so darzustellen, dass eine einfache, aussagekräftige Analyse der Lösung möglich ist. Zwei- und dreidimensionale Zusammenhänge sind vom Menschen noch erfassbar. Allerdings handelt es sich in der Regel um mehrdimensionale Zusammenhänge, womit für die Auswertung der Ergebnisse der multikriteriellen Optimierung Wege gefunden werden müssen, die eine einfache aber trotzdem aussagekräftige Auswertung ermöglichen. Insbesondere stellt sich die Aufgabe, aus den gefundenen Kompromissen einen bestimmten Kompromiss als letztendliche Lösung des Optimierungsproblems auszuwählen. Hier ist es insbesondere schwierig die Zusammenhänge der einzelnen Zielfunktionen und der Variationsgrößen in der gefundenen Lösung zu erfassen. Auch dazu gibt es schon bekannte Ver- fahren und Ansätze, um dieses Problem zu lösen. When calibrating vehicle control units (xCU), such as an engine control unit ECU or a TCU, on vehicle and component test benches, but also during test drives, the calibration engineer is regularly faced with the task of optimizing a compromise between different target values. Problems in which target values of several objective functions are simultaneously to be optimized depending on a number of restrictions are widespread in practice and in particular in calibration. Examples include the NOx-soot compromise, or the NOx-soot consumption compromise in classic internal combustion engine drives, or the battery damage- consumption compromise in hybrid applications, or the tradeoff between sporty but still comfortable tuning of gearshift operations. In general, in vehicles with drive concepts with automatic transmissions and / or hybridization, more and more conflicting demands are made on the individual components, which must be represented and optimized in the form of a compromise. Such multicriteric optimization problems generally exist if target values of several objective functions are to be optimized simultaneously depending on several restrictions (such as boundary conditions, physical limits, etc.). However, if the objective functions are in conflict with each other, the simultaneous optimization of all objective functions is often a problem. In particular, such multicriteria optimization problems usually can not find unique solutions but only a set of possible multidimensional solution points in a multi-dimensional space a multi-dimensional surface, the so-called Pareto front, which all represent an optimal compromise of the multi-criteria optimization problem. Individual points of this Pareto front therefore represent different, but in each case optimal, compromises between the objective functions. Such multi-criteria optimization problems are known per se and there are also a number of mathematical methods for solving such problems. From EP 2 192 294 A1, for example, a method is known with which an ECU is optimized during operation by means of a multi-criteria optimization problem with respect to an exhaust gas soot consumption compromise. In this case, a single aggregate objective function (AOF) is used, which combines the weighted objective functions into a functional one common solution is the linearly weighted summation of the objective functions described here. Each objective function is assigned a weighting factor. tor, from which a scalar objective function is derived as the sum. In this case, the actual optimization can be carried out using conventional approaches, for example by means of sequential quadratic programming (SQP), an effective, iterative method for nonlinear limited optimization, which is necessary for the desired reduction of the required computing power. However, the meaningfulness of such optimization largely depends on the choice of weighting factors, which in many cases can not be determined reliably in advance, which means that the results of such an optimization are not always satisfactory or even limited to a limited group from the optimization problem is applicable. Another problem with multi-criteria optimization problems is to visualize the solution found and present it in such a way that a simple, meaningful analysis of the solution is possible. Two- and three-dimensional relationships are still detectable by humans. However, as a rule, these are multi-dimensional correlations, with which methods for the evaluation of the results of the multi-criteria optimization must be found, which allow a simple but nevertheless meaningful evaluation. In particular, the task is to select from the compromises found a certain compromise as the ultimate solution of the optimization problem. Here, it is particularly difficult to grasp the relationships between the individual objective functions and the variables in the solution found. There are also known methods and approaches to solve this problem.
In US 7 921 371 B1 ist ein Verfahren zur Visualisierung der Lösung eines mehrdimensionalen multikriteriellen Optimierungsproblems beschrieben. Dabei werden alle Zielgrößen auf parallelen, nebeneinander liegenden Achsen dargestellt, wobei auch jeweils der kleinste und größte Wert der gefundenen optimalen Lösung visualisiert ist, die den gesamten möglichen Bereich der optimalen Lösungen für eine Zielgröße repräsentieren. Es wird daher die Pareto- Front in Form von parallelen Achsen für die Anzahl der Zielgrößen dargestellt. Für jede Zielgröße gibt es einen festgelegten Zielwert, der ebenfalls im Diagramm dargestellt ist, wobei die parallelen Achsen vertikal so verschoben sind, dass die Zielwerte alle auf einer horizontalen Linie zu liegen kommen. Die oberen und unteren Grenzen der einzelnen Zielgrößen können nun durch einen Benutzer variiert werden, wobei gleichzeitig auch alle anderen Achsen beeinflusst werden, also nur mehr die Lösungen dargestellt werden, die die vom Benutzer gewählte Grenzen einer Zielgröße erfüllen. Der Benutzer erhält damit die Möglichkeit die Zusammenhänge zwischen den einzelnen Zielgrößen zu analysieren und dadurch einen Kompromiss aus der Menge aller möglichen Kompromisse (Pareto-Front) auszuwählen. Aus der WO 01/67395 A1 ist wiederum ein Verfahren bekannt, bei dem in einer Darstellung aller möglichen optimalen Lösungen jene markiert (hier durch eine andere Farbe) werden, die ein bestimmtes, vom Benutzer vorgegebenes Kriterium, wie z.B. eine Randbedingung, erfüllen. Dabei werden die Lösungen auf zwei oder dreidimensionale Flächen projiziert, um Zusammenhänge erkennen zu können. Dargestellt werden sowohl Zielgrößen (abhängige Größen), als auch Variationsgrößen (unabhängige Größen). Damit kann die Lösungsmenge gezielt eingeschränkt werden und visuell dargestellt werden, was es ebenfalls erlaubt, Zusammenhänge zu analysieren und zu erkennen. Allerdings ermöglicht diese Methode nur sehr pauschale Aussagen zu treffen und erlaubt keine detaillierte Analyse der Lösung eines multikriteriellen Optimierungsproblems. US Pat. No. 7,921,371 B1 describes a method for visualizing the solution of a multi-dimensional multi-criteria optimization problem. All target values are displayed on parallel, adjacent axes, whereby also the smallest and largest value of the found optimal solution is visualized, which represent the entire possible range of the optimal solutions for a target size. Therefore, the pareto front is represented in the form of parallel axes for the number of target variables. For each target size, there is a fixed target value, also shown in the diagram, with the parallel axes shifted vertically so that the target values all lie on a horizontal line. The upper and lower limits of the individual target variables can now be varied by one user, whereby at the same time all the other axes are also influenced, ie only the solutions that meet the user-selected limits of a target variable are displayed. The user thus has the opportunity to analyze the relationships between the individual target values and thereby to choose a compromise from the set of all possible compromises (Pareto front). WO 01/67395 A1, in turn, discloses a method in which, in a representation of all possible optimal solutions, those are marked (here by a different color), which satisfy a certain user-specified criterion, such as a constraint. The solutions are projected onto two or three-dimensional surfaces in order to be able to recognize connections. Both target variables (dependent variables) and variables (independent variables) are shown. Thus, the amount of solution can be specifically limited and visually displayed, which also allows to analyze and recognize relationships. However, this method only allows to make very general statements and does not allow a detailed analysis of the solution of a multi-criteria optimization problem.
Es ist daher eine Aufgabe der gegenständlichen Erfindung ein Verfahren zur Visualisierung und Analyse der Lösung eines multikriteriellen Optimierungsproblems anzugeben, das eine einfache aber trotzdem detaillierte Auswertung der gefundenen Lösungen und damit die Auswahl einer der gefundenen, optimalen Lösungen als Lösung des multikriteriellen Optimierungsproblems ermöglicht. It is therefore an object of the present invention to provide a method for visualizing and analyzing the solution of a multi-criteria optimization problem, which allows a simple but nevertheless detailed evaluation of the solutions found and thus the selection of one of the found, optimal solutions as a solution of the multi-criteria optimization problem.
Diese Aufgabe wird erfindungsgemäß dadurch gelöst, indem die Menge der optimalen Lö- sungen des multikriteriellen Optimierungsproblems in einem Modellraum als zwei- oder dreidimensionales Diagramm der Zielfunktionen dargestellt sind und gleichzeitig in einem Variationsraum zumindest eine der Zielfunktionen in Abhängigkeit von zumindest einer Variationsgröße dargestellt ist und der Modellraum und der Variationsraum interaktiv miteinander verbunden sind, indem für jede selektierte Lösung im Modellraum die der Lösung zugrunde- liegende Variationsgröße im Variationsraum markiert wird. Diese Art der Aufbereitung und Darstellung des multikriteriellen Optimierungsproblems, indem die optimalen Lösungen des multikriteriellen Optimierungsproblems zusammen mit den Variationsgrößen dargestellt sind, erlaubt eine einfache und aussagekräftige grafische Analyse der Lösungen des multikriteriellen Optimierungsproblems mit dem Ziel eine Lösung als bestmöglicher Kompromiss auszuwählen. Durch die Koppelung der grafischen Ansichten von Variationsraum und Modellraum ergibt sich eine besonders hohe Aussagefähigkeit bei der Analyse multikriterieller Optimierungsprobleme. So können einzelne Variationsgrößen hinsichtlich Ihrer Optimalität und Abhängigkeit zu anderen Variationsgrößen analysiert werden. This object is achieved according to the invention by representing the set of optimal solutions of the multicriteria optimization problem in a model space as a two- or three-dimensional diagram of the objective functions and simultaneously displaying at least one of the objective functions as a function of at least one variation variable in a variation space Model space and the variation space are interactively interconnected by marked for each selected solution in the model space, the solution size underlying the solution in the variation space. This way of processing and presenting the multi-criteria optimization problem by presenting the optimal solutions of the multi-criteria optimization problem together with the variation sizes allows a simple and meaningful graphic analysis of the solutions of the multi-criteria optimization problem with the aim of choosing a solution as the best possible compromise. The coupling of the graphical views of the variation space and the model space results in a particularly high informative value in the analysis of multi-criteria optimization problems. Thus, individual variation quantities can be analyzed with regard to their optimality and dependence on other variables.
In einer bevorzugten Ausgestaltung kann als Zielfunktion ein mathematisches Modell ver- wendet werden, das aus einer Anzahl von Messungen der Zielfunktion in Abhängigkeit von den Variationsgrößen ermittelt wird. Das ermöglicht es, das erfindungsgemäße Verfahren praktisch auf beliebige Optimierungsprobleme anzuwenden. In a preferred embodiment, a mathematical model can be used as the objective function, which is determined from a number of measurements of the objective function as a function of the variables. This makes it possible to apply the inventive method practically to any optimization problems.
Ganz besonders vorteilhaft ist, im Variationsraum für die Zielfunktion zusätzlich einen Modellvertrauensbereich dazustehen, womit die Lösungen gleichzeitig in Bezug auf den jeweils zugehörigen Modellvertrauensbereich bewertet werden können. Insgesamt lassen sich auf diese Weise bessere Kalibrierergebnisse erzeugen und darüber hinaus lässt sich der Kalibriervorgang teilautomatisiert durchführen, d. h. in kosten- bzw. zeiteffektiver Weise und nicht zuletzt komfortabel und reproduzierbar. It is also particularly advantageous to additionally provide a model confidence range in the variation space for the objective function, with which the solutions can be evaluated simultaneously with respect to the respectively associated model confidence range. Overall, it is possible to produce better calibration results in this way, and moreover, the calibration process can be carried out semi-automatically, ie in a cost-effective and time-effective manner, and not least comfortable and reproducible.
Die gegenständliche Erfindung wird nachfolgend unter Bezugnahme auf die Figur 1 erläutert, die schematisch und beispielhaft eine vorteilhafte Ausgestaltung der Erfindung zeigt. The subject invention will be explained below with reference to Figure 1, which shows schematically and by way of example an advantageous embodiment of the invention.
Das Problem Zielgrößen zu finden, die mehrere Zielfunktionen gleichzeitig in Abhängigkeit mehrerer Einschränkungen optimieren, in der Regel minimieren, werden multikriterielle Optimierungsprobleme (MOP) genannt. Die mathematische Beschreibung des Problems lautet wie folgt: min/ max (f i(x),f 2(x),...,fk(x)) The problem of finding targets that typically optimize multiple target functions concurrently with multiple constraints is called Multi-Criteria Optimization Problems (MOP). The mathematical description of the problem is as follows: min / max (f i (x), f 2 (x), ..., f k (x))
x  x
(MOP) gi (x) < 0 i = l,...,m (MOP) gi (x) <0 i = l, ..., m
X mm■ < X < X max  X mm ■ <X <X max
Wobei fj(x), mit j=1 ,...,k, die Zielfunktion darstellt, die minimiert bzw. maximiert werden soll. Diese Optimierung soll sowohl unter Beachtung von Nebenbedingungen g, als auch für einen beschränkten Bereich xmin, xmax der Variationsgrößen x={x1,...xn} erfolgen. Where f j (x), where j = 1, ..., k, represents the objective function that should be minimized or maximized. This optimization should take place both under consideration of constraints g, as well as for a limited range x min , x max of the variables x = {x 1 , ... x n }.
Die Variationsgrößen x liegen im Variationsraum, unter dem man den n-dimensionalen Raum (n = Anzahl der Variationsgrößen) versteht, der durch die Variationsgrößen x aufgespannt wird. Die Variationsgrößen x sind z.B. bei einer Kalibrierung die Einstellungen, z.B. am Prüfstand, an denen der Kalibrateur seine Messungen vornimmt. Da diese Variationsgrößen x gültige Punkte im Variationsraum darstellen weiß der Kalibrateur, dass er die Variationsgrößen x in diesem Bereich verstellen kann. Aus diesem Grund wird in diesem Variati- onsraum um die Variationsgrößen x eine Hülle gelegt, die als Designraum bezeichnet wird. Der Designraum enthält somit alle für den jeweiligen Anwendungsfall gültigen Variationsgrößen x. The variation quantities x lie in the space of variation, under which one understands the n-dimensional space (n = number of variation quantities), which is spanned by the variation quantities x. The variation quantities x are e.g. during a calibration, the settings, e.g. at the test bench where the calibrator carries out his measurements. Since these variation quantities x represent valid points in the variation space, the calibrator knows that he can adjust the variation quantities x in this range. For this reason, in this variation space, around the variation quantities x, an envelope is laid, which is referred to as design space. The design room thus contains all the variables x valid for the respective application.
Wenn die Zielfunktionen , wie meistens der Fall, in einem Zielkonflikt zueinander stehen, ist es schwierig alle Zielfunktionen fj gleichzeitig zu minimieren. Deshalb wurde das Konzept der nichtdominierenden Individuen eingeführt. Ein Punkt im Zielfunktionsraum X* wird als nicht- dominierendes Individuum bezeichnet, sofern es keine gültige Lösung innerhalb des Designraumes gibt für die gilt: a) , = 1 k / f i (x) ^ f i (x*) b) 3 j e {1 , 2,..., k} / f j (x) < f j (x) In anderen Worten sofern es keine Möglichkeit gibt, den Punkt in einer der Zielfunktionen zu verbessen, ohne dabei eine der anderen zu verschlechtern. Üblicherweise und bekanntermaßen existiert für ein multikriterielles Optimierungsproblem nicht nur eine solche Lösung sondern eine Menge von möglichen Lösungen, die als Pareto-Front bezeichnet wird. Es existieren viele theoretische Ansätze zur Lösung solcher multikriterieller Optimierungsprobleme, jedoch liefern nicht alle bekannten Methoden gute Lösungen bei der Annäherung an die Pareto-Front. If the objective functions, such as most often the case, are in a conflict with each other, it is difficult to minimize all objective functions f j simultaneously. Therefore, the concept of non-dominant individuals was introduced. A point in the objective function space X * is called a non-dominating individual, unless there is a valid solution within the design space for which: a), = 1 k / fi (x) ^ fi (x * ) b) 3 each {1 , 2, ..., k} / fj (x) <fj (x) In other words, if there is no way to improve the point in one of the objective functions without worsening one of the others. Usually and as is known, for a multi-criteria optimization problem, not only such a solution exists, but a set of possible solutions called the Pareto front. There are many theoretical approaches to solving such multicriteria optimization problems, but not all known methods provide good approaches to the Pareto front.
Eine einzelne Aggregate Objective Function (AOF), wie z.B. in der EP 2 192 294 A1 beschrieben, ist vermutlich der intuitivste Ansatz zum Lösen multikriterieller Optimierungs- probleme. Evolutionärer Algorithmen zur multikriteriellen Optimierung benötigen hingegen keine Wichtungen und auch keine a priori Information, so dass diese Verfahren in den letzten Jahren vermehrt zum Einsatz gekommen sind und haben sich insbesondere als effektive und robuste Verfahren erwiesen. Speziell genetischen Algorithmen - basierend auf Selektion, Rekombination und Mutation - kamen hierbei zum Einsatz, die eine kontinuierliche Annähe- rung an ein gewünschtes Ziel realisieren. Sie sind leicht auf verschiedenste Probleme anzuwenden und sind sehr robust bei der Suche nach globalen Optima, auch wenn eine Vielzahl an lokalen Optima existiert. Bei der multikriteriellen Optimierung wird zu zueinander im Widerspruch stehenden Anforderungen eine Menge von Kompromisslösungen gesucht, welche die bestmöglichen Lösungen approximieren. Die Güte einer Approximation lässt sich durch das von ihr dominierte Volumen im Zielraum, der S Metrik, quantifizieren. Eine Maximierung der S-Metrik ist erstrebenswertes Ziel und gleichzeitig hinreichender skalarer Ersatz der ursprünglichen Zielfunktion. Ein genetischer Algorithmus setzt diese innerhalb der Selektion ein und erreicht dadurch hervorragende Ergebnisse, insbesondere wenn mehr als drei Ziele zu optimieren sind, wo andere multikriterielle genetische Algorithmen versagen. Der derzeit populärste genetische Algorithmus zur Ermittlung der Pareto-Front ist der Non-Dominated Sorting Genetic Algorithm (NSGA-II), der sich insbesondere zur Ermittlung eines möglichst globalen Optimums als äußerst effizient erwiesen hat. NSGA-II ist ein performanter, elitisti- scher Algorithmus welcher nichtdominierende Individuen priorisiert und die Diversität der Lösungen bewahrt. Der Algorithmus erstellt eine Initialpopulation innerhalb des Vektors von Variationsgrößen x und nähert sich in einem iterativen Prozess, welcher auf Selektions-,A single Aggregate Objective Function (AOF), such as described in EP 2 192 294 A1, is probably the most intuitive approach to solve multi-criteria optimization problems. By contrast, evolutionary algorithms for multicriteria optimization do not require any weighting or a priori information, so that these methods have been increasingly used in recent years and in particular have proven to be effective and robust methods. Specifically genetic algorithms - based on selection, recombination and mutation - were used to achieve a continuous approach to a desired target. They are easy to apply to a wide range of problems and are very robust in the search for global optima, even if there are a variety of local optima. In multicriteria optimization, conflicting requirements are searched for in a set of compromise solutions that approximate the best possible solutions. The quality of an approximation can be quantified by the volume dominated by it in the target space, the S metric. Maximizing the S-metric is a desirable goal and, at the same time, sufficient scalar replacement of the original objective function. A genetic algorithm uses these within the selection and thus achieves excellent results, especially when more than three goals have to be optimized where other multicriteria genetic algorithms fail. The currently most popular genetic algorithm for determining the Pareto front is the Non-Dominated Sorting Genetic Algorithm (NSGA-II), which has proven to be extremely efficient, especially for the determination of a global optimum. NSGA-II is a high-performance, elitist algorithm that prioritizes non-dominating individuals and preserves the diversity of solutions. The algorithm creates an initial population within the vector of variational magnitudes x and approaches iteratively in an iterative process that is based on selection,
Kreuzungs- und Mutationsoperationen basiert, mittels Favorisierung von nichtdominierenden Individuen (Elitismus) mit jeder neuen Generation näher an, bis ein Lösungskriterium erfüllt ist. Solche genetischen Algorithmen sind naturgemäß verhältnismäßig rechenintensiv, jedoch auch leicht für heutige Mehrkernprozessoren und verteilte Rechnerstrukturen zu paral- lelisieren. Da diese Algorithmen an sich bekannt sind, wird hier nicht näher darauf eingegangen, insbesondere da das gewählte Verfahren zur Lösung des mulitkriteriellen Optimie- rungsproblems keinen Einfluss auf die erfindungsgemäße Aufbereitung und Analyse der Ergebnisse hat. Crossing and mutation operations are closer to each new generation by favoring non-dominant individuals (elitism) until a solution criterion is met. Such genetic algorithms are by nature relatively computationally intensive, but also easy to parallelize for today's multi-core processors and distributed computer structures. Since these algorithms are known per se, they will not be discussed in more detail here, especially since the chosen method for the solution of the multi-criteria optimization has no influence on the preparation and analysis of the results according to the invention.
Um das Ergebnis der multikriteriellen Optimierung (unabhängig vom gewählten Lösungsalgorithmus) auswerten zu können, wird eine spezielle Art der Aufbereitung und Darstellung der verwendeten Größen genutzt. Diese Auswertung auf Basis einer grafischen Analyse wird nachfolgend beschrieben. In order to be able to evaluate the result of the multi-criteria optimization (independent of the selected solution algorithm), a special type of processing and representation of the variables used is used. This evaluation based on a graphical analysis is described below.
Die Besonderheit an der visuellen Aufbereitung liegt in einer geteilten Darstellung von Modellraum 1 und Variationsraum 2, wie in Fig. 1 dargestellt. Der Modellraum 1 ist dabei der k- dimensionale Raum, der durch die k Zielfunktionen (χ), j=1 ,...,k aufgespannt wird. Im Bei- spiel nach Fig. 1 z.B. der dreidimensionale Raum, der durch die Zielfunktionen f-i(x)„Smoke", f2(x)„NOX" und f3(x)„Fuel Consumption" aufgespannt wird. Die Abbildung des Designraumes in diesen Modellraum 1 ist der gültige Bereich 3, in dem sich die Lösungen bewegen können. Die Pareto-Front 4 enthält die gefundenen Lösungen des multikriteriellen Optimierungsproblems innerhalb dieses gültigen Bereichs 3. Falls der Modellraum 1 eine höhere Dimension als drei aufweist, so kann der Modellraum 1 auch durch mehrere zwei- oder dreidimensionale Darstellungen der k Dimensionen des Modellraumes dargestellt werden. Welche der k Dimensionen in den verschiedenen zwei- oder dreidimensionale Darstellungen zusammengefasst werden, kann vom multikriteriellen Optimierungsproblem und von der Präferenz des Benutzers abhängig gemacht werden. Der Variationsraum 2 wird durch eine Anzahl zwei oder dreidimensionalen Darstellungen von Zielfunktionen fj(x) und Variationsgrößen x repräsentiert. Im Beispiel nach Fig.1 ist die Zielfunktion f3(x) jeweils in Abhängigkeit der drei Variationsgrößen x-ι„Abgastemperatur", x2 „EGR Rate" und x3„Raildruck" dargestellt. Hier sind aber beliebige Kombinationen von Zielfunktionen fj(x) und Variationsgrößen x denkbar. Die Zielfunktionen fj(x) können dabei bekannte Funktionen der Variationsgrößen x sein. Es ist aber auch denkbar, dass eine Zielfunktion fj(x) ein mathematisches Modell ist, das aus Messungen oder Versuchen ermittelt wird. Dabei werden am Gegenstand des multikriteriellen Optimierungsproblems, also z.B. an einem Verbrennungsmotor, einem Antriebsstrang, einem Getriebe, einem Fahrzeug, etc., Messungen durchgeführt, z.B. auf entspre- chenden Prüfständen oder im Zuge von Testfahrten. Dabei werden die gewünschten Zielfunktionen fj in Abhängigkeit von den Variationsgrößen x und eventuell anderen Größen gemessen. Aus diesen Messgrößen werden dann mathematische Modelle der Zielfunktionen fj erstellt. Dazu gibt es ebenso eine Reihe von bekannten Methoden, mit denen eine möglichst gute Abdeckung mit möglichst wenigen Messungen erreicht werden können, um möglichst gute mathematischen Modelle zu erhalten. Mögliche Modelle sind z.B. ein polynomisches Regressionsmodell, ein Fast Neural Network oder ein Intelligent Neural Network. Aufgrund dieser Vorgehensweise werden zusätzliche Messungen, also reale Messwerte, nicht unbedingt 100% exakt auf diesem Modell liegen. Die Methoden zur Ermittlung der Modelle liefern daher auch einen Modellvertrauensbereich, der die Bandbreite angibt, in der sich weitere Messungen wahrscheinlich bewegen. Das bedeutet, dass ein Modell mit einem schlanken Modellvertrauensbereich relativ gut auf die gemachten Messungen passt und eine dementsprechend gute Aussagekraft hat. Je enger der Modellvertrauensbereich desto besser passen die Modelle auf die Messungen und desto wahrscheinlicher ergeben die über die Model- le ermittelten Lösungen des multikriteriellen Optimierungsproblems tatsächliche die gesuchten Werte. Selbstverständlich können auch bekannte Zielfunktionen (χ) einen Modellvertrauensbereich aufweisen, der wiederum angibt, wie genau eine Zielfunktion fj mit einer realen Messung übereinstimmt. Der Modellvertrauensbereich ist somit ein Maß für die Genauigkeit des Modells bzw. einer Zielfunktion bezogen auf reale Messungen. Im Variationsraum 2 können daher in den einzelnen Diagrammen auch die Modellvertrauensbereiche 5 dargestellt sein, z.B. in Form einer oberen und unteren Grenze wie in Fig.1 ersichtlich. The peculiarity of the visual presentation lies in a split representation of model space 1 and variation space 2, as shown in FIG. The model space 1 is the k-dimensional space, which is spanned by the k target functions (χ), j = 1, ..., k. In examples 1, for example f 2 (x) "NOX" and f 3 (x) "Fuel Consumption" The illustration of the design space game of FIG. Three-dimensional space, the fi through the objective functions (x) "Smoke", is spanned. the valid range 3 in which the solutions can move is contained in this model space 1. The Pareto front 4 contains the solutions found for the multi-criteria optimization problem within this valid range 3. If the model space 1 has a dimension larger than three, then the Model space 1 can also be represented by several two- or three-dimensional representations of the k dimensions of the model space.These dimensions of the k dimensions in the two- or three-dimensional representations can be made dependent on the multi-criteria optimization problem and the preference of the user is represented by a number of two or three-dimensional representations of objective functions f j (x) and variation quantities x In Place. In the example of Figure 1, the objective function is f 3 (x) each represented as a function of the three varying sizes x-ι "exhaust gas temperature" x 2 "EGR rate" and x 3 "rail pressure". However, any combinations of objective functions f j (x) and variation quantities x are conceivable here. The objective functions f j (x) can be known functions of the variables x. However, it is also conceivable that an objective function f j (x) is a mathematical model that is determined from measurements or experiments. In this case, measurements are carried out on the object of the multi-criteria optimization problem, eg on an internal combustion engine, a drive train, a transmission, a vehicle, etc., for example on corresponding test benches or in the course of test drives. In this case, the desired objective functions f j are measured as a function of the variables x and possibly other variables. Mathematical models of the objective functions f j are then created from these measured variables. There are also a number of known methods with which the best possible coverage can be achieved with as few measurements as possible to get good mathematical models. Possible models include a polynomial regression model, a fast neural network or an intelligent neural network. Because of this approach, additional measurements, ie real measurements, will not necessarily be 100% accurate on this model. The methods for determining the models therefore also provide a model confidence interval indicating the bandwidth in which further measurements are likely to move. This means that a model with a slim model confidence range fits relatively well on the measurements made and has accordingly good explanatory power. The closer the model confidence interval, the better the models fit the measurements, and the more likely the model-determined solutions to the multi-criteria optimization problem will actually yield the sought-after values. Of course, well-known objective functions (χ) can also have a model confidence range, which in turn indicates how exactly an objective function f j agrees with a real measurement. The model confidence range is thus a measure of the accuracy of the model or a target function with respect to real measurements. In the variation space 2, therefore, the model confidence ranges 5 can also be shown in the individual diagrams, for example in the form of an upper and lower limit as shown in FIG.
Das besondere an dieser Art der Darstellung ist, dass dadurch die Zielfunktionen (χ) im Modellraum 1 bzw. die Pareto-Front 4 als Menge der möglichen optimalen Lösungen des multikriteriellen Optimierungsproblems zusammen mit den Variationsgrößen x dargestellt sind und daher auch gemeinsam analysiert werden können. Dabei wird die Darstellung im Variationsraum 2 interaktiv an eine Auswahl eines Punktes im Modellraum 1 angepasst. Dazu ist z.B. ein Fadenkreuz 6 vorgesehen. Mit dem z.B. ein interessanter Punkt 7 der Pareto- Front 4, oder des gültigen Bereichs 3, selektiert wird. Im Variationsraum 2 markiert das Fa- denkreuz 6 automatisch die Variationsgrößen x für diesen Punk 7 im Modellraum 1 . Gleichzeitig können auch die jeweiligen Werte der Variationsgrößen x in diesem Punkt angegeben werden, wie in Fig. 1 angedeutet. Ebenso kann der Modellvertrauensbereich 5 dargestellt sein, sodass der Benutzer zusätzlich Information erhält, wie vertrauenswürdig die zugrundeliegende Zielfunktion fj (bzw. mathematische Modell) in diesem Punkt ist. Einmal ermittelte Modelle können natürlich auch für spätere Aufgaben eingesetzt werden. The special feature of this type of representation is that the objective functions (χ) in the model space 1 and the Pareto front 4 are represented as a set of possible optimal solutions of the multi-criteria optimization problem together with the variables x and therefore can also be analyzed together. The representation in the variation space 2 is interactively adapted to a selection of a point in the model space 1. For this example, a crosshair 6 is provided. With the example, an interesting point 7 of the Pareto front 4, or the valid area 3, is selected. In the variation space 2, the crosshair 6 automatically marks the variation quantities x for this punk 7 in the model space 1. At the same time, the respective values of the variables x at this point can also be indicated, as indicated in FIG. 1. Likewise, the model trust area 5 can be represented so that the user additionally receives information about how trustworthy the underlying objective function fj (or mathematical model) is on this point. Once determined models can of course be used for later tasks.
Durch die Koppelung der grafischen Ansichten von Variationsraum 2 und Modellraum 1 ergibt sich eine besonders hohe Aussagefähigkeit bei der Analyse multikriterieller Optimierungsprobleme. So können einzelne Variationsgrößen x nicht nur hinsichtlich Ihrer Optimali- tät und Abhängigkeit zu anderen Variationsgrößen x, sondern gleichzeitig in Bezug auf den jeweils zugehörigen Modellvertrauensbereich 5 bewertet werden. Somit wurde die Möglich- keit geschaffen, die Lösungsmenge im Detail zu analysieren, um den tatsächlich bestmöglichen Kompromiss zwischen den Zielfunktionen (χ) zu ermitteln. The coupling of the graphical views of Variationsraum 2 and Modellraum 1 results in a particularly high informative value in the analysis of multi-criteria optimization problems. Thus, individual variation quantities x can be evaluated not only with regard to their optimality and dependence on other variables x, but also simultaneously with respect to the respectively associated model confidence range 5. Thus, the possible to analyze the solution set in detail in order to find the best possible compromise between the objective functions (χ).
Diese Art der Darstellung erlaubt aber auch die Analyse der Einflüsse von Änderungen in der Vorgabe der Nebenbedingungen g, und/oder des Bereichs xmin, xmax der Variationsgrößen x. Solche Änderungen führen zu anderen Lösungen, die dann einfach direkt miteinander verglichen werden können. Dazu kann z.B. vorgesehen sein, die Grenzen des Bereichs der Variationsgrößen x z.B. mittels eines Schiebers im Variationsraum 2 zu ändern, womit sich gleichzeitig die Darstellung der Lösung im Modellraum 1 ändern kann. Ebenso könnte vorgesehen sein, die Lösungen im Modellraum 1 so zu filtern, dass nurm ehr solche angezeigt werden, die bestimmte Bereiche einer oder mehrerer Variationsgrößen x erfüllen. Durch die gekoppelte Darstellung und die Interaktivität der Darstellung sind solche Einflüsse einfacher zu erfassen. However, this type of representation also allows the analysis of the effects of changes in the specification of the constraints g, and / or the range x min , x ma x of the variables x. Such changes lead to other solutions that can then be easily compared directly. For this purpose, it can be provided, for example, to change the limits of the range of variation variables x, for example by means of a slider in the variation space 2, whereby the representation of the solution in the model space 1 can change at the same time. It could also be provided to filter the solutions in the model space 1 such that only those that fulfill certain areas of one or more variables x are displayed. The coupled presentation and the interactivity of the presentation make such influences easier to grasp.
Aufgrund der Komplexität der in diesem Zusammenhang notwendigen Berechnungen ist es vorteilhaft, die Softwarearchitektur der erfindungsgemäßen Methode so zu gliedern, dass sich komplexe Aufgaben parallelisieren lassen. Dadurch können einzelne Aufgaben parallel auf verschiedenen Prozessoren oder auch auf verschiedenen Rechnern ausgeführt werden. Diese verteilte multikriterielle Optimierung erlaubt dem Anwender seine Optimierungsaufgaben performant und hochqualitativ in einem beliebig skalierbaren, verteilten System durchzuführen. Diese Steigerung der Performance ermöglicht darüber hinaus mit der gekoppelten Visualisierung von Variationsraum 2 und Modellraum 1 eine deutlich gesteigerte Aussagefähigkeit bei der Analyse multikriterieller Optimierungsprobleme. Due to the complexity of the calculations required in this context, it is advantageous to structure the software architecture of the method according to the invention so that complex tasks can be parallelized. As a result, individual tasks can be executed in parallel on different processors or on different computers. This distributed multi-criteria optimization allows the user to carry out his optimization tasks in a high-quality, high-performance manner in an arbitrarily scalable, distributed system. In addition, with the coupled visualization of variation space 2 and model space 1, this increase in performance enables significantly increased informative value in the analysis of multi-criteria optimization problems.
Ein möglicher Ablauf einer Kalibrierung eines Motorsteuergeräts ECU eines Verbrennungsmotors hinsichtlich des NOx-Ruß-Verbrauchs wird nachfolgend als Beispiel einer A possible sequence of a calibration of an engine control unit ECU of an internal combustion engine with respect to the NOx-soot consumption will be described below as an example
multikriteriellen Optimierung beschrieben. Zuerst werden eine Anzahl von Messungen am Verbrennungsmotor durchgeführt, wobei die Zielgrößen der Zielfunktionen fj(x) NOx, Ruß und Verbrauch in Abhängigkeit der Variationsgrößen x, z.B. Abgastemperatur, EGR Rate, Raildruck, gemessen werden. Die Anzahl und die Menge der Messungen kann vorgegeben sein, z.B. durch ein vorgegebenes Design of Experiment. Anhand der Messungen werden mathematische Modelle und Modellvertrauensbereiche 5 für die Zielfunktionen fj(x) ermittelt. Danach kann das multikriterielle Optimierungsproblem zur Optimierung der Zielfunktionen fj(x) gelöst und die Lösung in der geteilten Darstellung von Modellraum 1 und Variationsraum 2 analysiert werden. Dabei kann der Kalibrateur verschiedene optimale Lösungen der Pare- to-Front 4 hinsichtlich der zugrundeliegenden Variationsgrößen x und des Modellvertrauensbereichs 5 untersuchen. Aus diesen möglichen optimalen Lösungen bestimmt der Kalibrateur dann eine der Lösungen als bestmöglichen Kompromiss. Hierbei spielt die Erfahrung des Kalibrateurs eine große Rolle. Dazu können neben den Modellvertrauensbereichen und den Abhängigkeiten der Variationsgrößen x auch die Werte zusätzlicher Modellkanäle die nicht als Zielfunktionen optimiert wurden, sowie die Robustheit der Einstellungen, z.B. ob sich die Modelle in der Nähe des Optimums stark ändert, geringe Beeinflussbarkeit durch Bauteiltoleranzen, etc., berücksichtigt werden. Das kann für alle für die Kalibrierung benötigten Be- triebspunkte (z.B. Drehzahl, Drehmoment, Last) des Verbrennungsmotors wiederholt werden. Für eine Kalibrierung werden in der Regel eine vorgegebene Anzahl von Betriebspunkten, z.B. zehn bis zwanzig Betriebspunkte, benötigt. described multi-criteria optimization. First, a number of measurements are performed on the internal combustion engine, wherein the target variables of the target functions f j (x) NOx, soot and consumption depending on the variables x, eg exhaust gas temperature, EGR rate, rail pressure, are measured. The number and the amount of measurements can be predetermined, eg by a given Design of Experiment. The measurements are used to determine mathematical models and model confidence regions 5 for the objective functions f j (x). Then the multicriteria optimization problem for optimizing the objective functions f j (x) can be solved and the solution in the split representation of model space 1 and variation space 2 can be analyzed. In this case, the calibrator can examine various optimal solutions of the pareto-front 4 with regard to the underlying variation variables x and the model confidence region 5. From these possible optimal solutions, the calibrator then determines one of the solutions as the best possible compromise. Here, the experience of the calibrator plays a major role. In addition to the model confidence ranges and the Dependencies of the variables x are also the values of additional model channels that have not been optimized as target functions, as well as the robustness of the settings, eg whether the models in the vicinity of the optimum changes greatly, low influenceability due to component tolerances, etc., are taken into account. This can be repeated for all operating points required for the calibration (eg speed, torque, load) of the internal combustion engine. For a calibration, a predetermined number of operating points, eg ten to twenty operating points, are generally required.

Claims

Patentansprüche claims
1 . Verfahren zur Auswertung der Lösung eines multikriteriellen Optimierungsproblems für j=1 ,...,k von Variationsgrößen x abhängigen Zielfunktionen (χ), wobei die Menge der opti- malen Lösungen des multikriteriellen Optimierungsproblems in einem Modellraum (1 ) als zwei- oder dreidimensionales Diagramm der j=1 ,...,k Zielfunktionen fj(x) dargestellt sind und gleichzeitig in einem Variationsraum (2) zumindest eine der Zielfunktionen fj(x) in Abhängigkeit von zumindest einer Variationsgröße (x) dargestellt ist und der Modellraum (1 ) und der Variationsraum (2) interaktiv miteinander verbunden sind, indem für jede selektierte Lösung im Modellraum (1 ) die der Lösung zugrundeliegende Variationsgröße (x) im Variationsraum (2) markiert wird. 1 . Method for evaluating the solution of a multi-criteria optimization problem for j = 1,..., K of variance variables x-dependent objective functions (χ), where the set of optimal solutions of the multi-criteria optimization problem in a model space (1) is a two- or three-dimensional diagram the j = 1, ..., k target functions f j (x) are shown and at the same time in a variation space (2) at least one of the objective functions f j (x) is represented as a function of at least one variation variable (x) and the model space ( 1) and the variation space (2) are interactively connected to one another by marking for each selected solution in the model space (1) the solution variable (x) on which the solution is based in the variation space (2).
2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, dass als eine Zielfunktion fj(x) ein mathematisches Modell verwendet wird, das aus einer Anzahl von Messungen der Zielfunktion fj(x) in Abhängigkeit von den Variationsgrößen (x) ermittelt wird. 2. The method according to claim 1, characterized in that as a target function f j (x) a mathematical model is used which is determined from a number of measurements of the objective function f j (x) as a function of the variables (x).
3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass im Variationsraum (2) für die Zielfunktion fj(x) zusätzlich ein Modellvertrauensbereich (5) dargestellt ist. 3. The method according to claim 1 or 2, characterized in that in the variation space (2) for the objective function f j (x) additionally a model confidence region (5) is shown.
EP12795756.1A 2011-12-12 2012-11-08 Method for evaluating the solution to a multicriteria optimization problem Ceased EP2791827A2 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
AT0182011A AT510328A2 (en) 2011-12-12 2011-12-12 METHOD FOR EVALUATING THE SOLUTION OF A MULTICRITERIAL OPTIMIZATION PROBLEM
PCT/EP2012/072165 WO2013087307A2 (en) 2011-12-12 2012-11-08 Method for evaluating the solution to a multicriteria optimization problem

Publications (1)

Publication Number Publication Date
EP2791827A2 true EP2791827A2 (en) 2014-10-22

Family

ID=45724213

Family Applications (1)

Application Number Title Priority Date Filing Date
EP12795756.1A Ceased EP2791827A2 (en) 2011-12-12 2012-11-08 Method for evaluating the solution to a multicriteria optimization problem

Country Status (5)

Country Link
US (1) US9760532B2 (en)
EP (1) EP2791827A2 (en)
JP (1) JP5940667B2 (en)
AT (1) AT510328A2 (en)
WO (1) WO2013087307A2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4235319A1 (en) 2022-02-25 2023-08-30 Rheinisch-Westfälische Technische Hochschule (RWTH) Aachen Automated feature calibration

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150202464A1 (en) * 2014-01-23 2015-07-23 Mitsubis Multi-Criteria Optimization in Particle Beam Dose Optimization
JP6237321B2 (en) * 2014-02-21 2017-11-29 富士通株式会社 Operation plan creation program, operation plan creation method, and operation plan creation device
CN107002576B (en) * 2014-11-17 2020-10-23 大众汽车有限公司 Control apparatus for internal combustion engine
DE102015202135A1 (en) * 2015-02-06 2016-08-11 Volkswagen Aktiengesellschaft Control method and control for an internal combustion engine
DE102016106976B4 (en) * 2016-04-15 2018-10-31 Iav Gmbh Ingenieurgesellschaft Auto Und Verkehr Method for determining a model of a technical system
DE102016208238A1 (en) * 2016-05-12 2017-11-16 Volkswagen Aktiengesellschaft Control method for a hybrid drive, control unit and hybrid drive
DE102016208236A1 (en) * 2016-05-12 2017-11-30 Volkswagen Ag Control method for an internal combustion engine, control unit and internal combustion engine
AT518676B1 (en) * 2016-05-17 2018-02-15 Avl List Gmbh Method for calibrating a technical system
DE102017215251B4 (en) * 2017-08-31 2019-04-18 Volkswagen Aktiengesellschaft Method and control unit for emission control of an internal combustion engine
US10884721B2 (en) 2018-05-08 2021-01-05 Autodesk, Inc. Branch objects for dependent optimization problems
AT521422B1 (en) * 2018-06-21 2020-09-15 Avl List Gmbh Procedure for calibrating a technical system
EP3612011A1 (en) * 2018-08-14 2020-02-19 ABB Schweiz AG Method of controlling cooling in a data centre
AT521927B1 (en) 2018-12-10 2020-10-15 Avl List Gmbh Procedure for the calibration of a technical system
CN110533263A (en) * 2019-09-11 2019-12-03 绍兴建元电力集团有限公司 A kind of hot integrated system Multipurpose Optimal Method of electric-gas-based on improvement NSGA-II algorithm
CN111813669B (en) * 2020-07-04 2023-10-13 毛澄映 Adaptive random test case generation method based on multi-target group intelligence
CN113325896B (en) * 2021-05-31 2022-03-01 浙江工业大学 Multi-target temperature optimization control method of intelligent retail machine
IT202100020744A1 (en) * 2021-08-02 2023-02-02 Fpt Motorenforschung Ag Method of modeling a powertrain and controlling the modeled powertrain

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6417852B2 (en) 2000-03-09 2002-07-09 Vladimir Sevastyanov Method of visualization and graphical analysis for multidimensional functions
US6954908B2 (en) * 2002-12-10 2005-10-11 Cadence Design Systems, Inc. Circuit design point selection method and apparatus
US7921371B1 (en) 2006-03-22 2011-04-05 Versata Development Group, Inc. System and method of interactive, multi-objective visualization
JP5003499B2 (en) * 2008-01-14 2012-08-15 富士通株式会社 Multi-objective optimization design support apparatus, method, and program
JP5151733B2 (en) * 2008-06-27 2013-02-27 富士通株式会社 Design improvement support apparatus, method, and program for multi-objective optimization design
EP2192294B1 (en) 2008-11-27 2013-02-13 Continental Automotive GmbH Method for operating auto ignition combustion engine
JP5375429B2 (en) * 2009-02-12 2013-12-25 富士通株式会社 Design support apparatus, method, and program
JP2011253477A (en) * 2010-06-04 2011-12-15 Hitachi Ltd Design support device and design support method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
KOLLAT ET AL: "A framework for Visually Interactive Decision-making and Design using Evolutionary Multi-objective Optimization (VIDEO)", ENVIRONMENTAL MODELLING & SOFTWARE, ELSEVIER, AMSTERDAM, NL, vol. 22, no. 12, 8 August 2007 (2007-08-08), pages 1691 - 1704, XP022191725, ISSN: 1364-8152, DOI: 10.1016/J.ENVSOFT.2007.02.001 *
See also references of WO2013087307A2 *
T. KRUSE ET AL: "Modern Statistical Modeling and Evolutionary Optimization Methods for the Broad Use in ECU Calibration", PROCEEDINGS OF THE 17TH WORLD CONGRESS THE INTERNATIONAL FEDERATION OF AUTOMATIC CONTROL; SEOUL, KOREA; JULY 6-11, 2008., vol. 43, 12 July 2010 (2010-07-12), Red Hook, NY, pages 739 - 743, XP055464923, ISSN: 1474-6670, ISBN: 978-1-123-47890-7, DOI: 10.3182/20100712-3-DE-2013.00031 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4235319A1 (en) 2022-02-25 2023-08-30 Rheinisch-Westfälische Technische Hochschule (RWTH) Aachen Automated feature calibration
DE102022104648A1 (en) 2022-02-25 2023-08-31 Rheinisch-Westfälische Technische Hochschule (RWTH) Aachen, Körperschaft des öffentlichen Rechts AUTOMATED FUNCTIONAL CALIBRATION

Also Published As

Publication number Publication date
JP5940667B2 (en) 2016-06-29
AT510328A2 (en) 2012-03-15
US9760532B2 (en) 2017-09-12
WO2013087307A3 (en) 2013-09-06
US20140344320A1 (en) 2014-11-20
WO2013087307A2 (en) 2013-06-20
JP2014533387A (en) 2014-12-11

Similar Documents

Publication Publication Date Title
EP2791827A2 (en) Method for evaluating the solution to a multicriteria optimization problem
EP2999998B1 (en) Methods for ascertaining a model of an output variable of a technical system
DE202010017612U1 (en) Improved genetic algorithm with multiple targets for component design optimization
DE102008048478A1 (en) Sampling strategy using genetic algorithms in the optimization of a technical design
DE112006000821T5 (en) Tax system and tax procedure
EP1626317A2 (en) Method for developing a nonlinear, stationary or dynamical model of a control variable of a construction machine
DE102020202335A1 (en) Method and device for creating a model of a technical system from measurements
WO2020057937A1 (en) Method and assistance system for parameterizing an anomaly detection method
WO2016146528A1 (en) Method for generating a model ensemble for calibrating a control device
EP3458699B1 (en) Method for calibrating a technical system
DE10209146A1 (en) Method and system for the automatic planning of experiments
WO2016198047A1 (en) Method for establishing a simulation model for mapping at least one functional process of a drive train component
DE102009021781A1 (en) Engine-operating method for calculating an engine-operating map for a vehicle&#39;s control device creates a map with a specified number of nodes while measuring data points to calculate a map value
WO2009053137A2 (en) Method for the computer-aided ascertainment of at least one property of a hair colouration
DE102014225039A1 (en) Method and apparatus for providing sparse Gaussian process models for calculation in an engine control unit
DE102019128655A1 (en) Method for providing a computer-aided control for a technical system
Wei et al. Visual diagnostics of parallel performance in training large-scale dnn models
EP3942372B1 (en) Method for validating system parameters of an energy system, method for operating an energy system, and energy management system for an energy system
DE102020003746A1 (en) CREATION OF CHARACTERISTICS FOR A CONTROL OF A COMBUSTION ENGINE
DE102016106976B4 (en) Method for determining a model of a technical system
EP2996001A1 (en) Method for computer-based analysis of a data set from observations
DE602004011732T2 (en) Control method for achieving a linear shear response
AT521865B1 (en) Method and system for simulating a number of particles emitted by an internal combustion engine
EP4341876A1 (en) Computer-implemented method and system for determining optimized system parameters of a technical system using a cost function
DE102022104648A1 (en) AUTOMATED FUNCTIONAL CALIBRATION

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20140129

AK Designated contracting states

Kind code of ref document: A2

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

RIN1 Information on inventor provided before grant (corrected)

Inventor name: WALLNER, KLEMENS

Inventor name: DRAGOTI, ADNAND

Inventor name: GARCIA, ALEJANDRA

DAX Request for extension of the european patent (deleted)
17Q First examination report despatched

Effective date: 20151216

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

APBK Appeal reference recorded

Free format text: ORIGINAL CODE: EPIDOSNREFNE

APBN Date of receipt of notice of appeal recorded

Free format text: ORIGINAL CODE: EPIDOSNNOA2E

APBR Date of receipt of statement of grounds of appeal recorded

Free format text: ORIGINAL CODE: EPIDOSNNOA3E

APAF Appeal reference modified

Free format text: ORIGINAL CODE: EPIDOSCREFNE

REG Reference to a national code

Ref country code: DE

Ref legal event code: R003

APBT Appeal procedure closed

Free format text: ORIGINAL CODE: EPIDOSNNOA9E

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION HAS BEEN REFUSED

18R Application refused

Effective date: 20220628