WO2017018227A1 - Dispositif de contrôle optimal, procédé de contrôle optimal, programme informatique et système de contrôle optimal - Google Patents

Dispositif de contrôle optimal, procédé de contrôle optimal, programme informatique et système de contrôle optimal Download PDF

Info

Publication number
WO2017018227A1
WO2017018227A1 PCT/JP2016/070821 JP2016070821W WO2017018227A1 WO 2017018227 A1 WO2017018227 A1 WO 2017018227A1 JP 2016070821 W JP2016070821 W JP 2016070821W WO 2017018227 A1 WO2017018227 A1 WO 2017018227A1
Authority
WO
WIPO (PCT)
Prior art keywords
control
extreme value
parameter
value control
amount
Prior art date
Application number
PCT/JP2016/070821
Other languages
English (en)
Japanese (ja)
Inventor
理 山中
由紀夫 平岡
Original Assignee
株式会社東芝
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
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=57885542&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=WO2017018227(A1) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by 株式会社東芝 filed Critical 株式会社東芝
Priority to CN201680044271.4A priority Critical patent/CN107924162B/zh
Publication of WO2017018227A1 publication Critical patent/WO2017018227A1/fr

Links

Images

Classifications

    • C—CHEMISTRY; METALLURGY
    • C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F1/00—Treatment of water, waste water, or sewage
    • C—CHEMISTRY; METALLURGY
    • C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F3/00—Biological treatment of water, waste water, or sewage
    • C02F3/34—Biological treatment of water, waste water, or sewage characterised by the microorganisms used
    • 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

Definitions

  • Embodiments described herein relate generally to an optimal control device, an optimal control method, a computer program, and an optimal control system.
  • Extreme value control is a model-free real-time optimal control technique that does not use a complex model of a plant.
  • the outline of extreme value control is to search for an operation amount that optimizes an evaluation amount based on a control amount of a process to be controlled by forcibly changing the operation amount.
  • control parameters various parameters related to extreme value control
  • JP 2012-215575 A Japanese Patent Publication No. 8-23332 Japanese Patent Publication No. 8-23329 Japanese Patent Publication No. 6-60594 JP 2009-258068 A Japanese Patent No. 4286880 Japanese Patent No. 4309326 Japanese Patent No. 5300827 Japanese Patent Laid-Open No. 2004-171531
  • the problem to be solved by the present invention is to provide an optimum control device, an optimum control method, a computer program, and an optimum control system capable of executing extreme value control with a control parameter corresponding to the characteristics of a process to be controlled. .
  • the optimal control system of the embodiment includes a control target parameter determination unit, an extreme value control parameter determination unit, and an extreme value control unit.
  • the control target parameter determination unit is configured to control the control target process based on an operation amount of the control target process and an evaluation amount indicating an index related to optimization of the control target process based on the control amount that changes according to the operation amount.
  • the control target parameter indicating the characteristic is determined.
  • the extreme value control parameter determining unit determines an extreme value control parameter for executing the extreme value control based on the controlled object parameter determined by the controlled object parameter determining unit.
  • the extreme value control unit executes extreme value control of the process to be controlled using the extreme value control parameter determined by the extreme value control parameter determination unit.
  • the extreme value control unit executes extreme value control for changing the manipulated variable so that the evaluation value is directed to an optimum value using the extreme value control parameter determined by the extreme value control parameter determining unit.
  • the figure explaining the concept of extreme value control The block diagram which shows the structural example of an extreme value controller.
  • the functional block diagram which shows the function structure of the optimal control apparatus of embodiment.
  • FIG. 1 is a diagram for explaining the concept of extreme value control.
  • Extreme value control is a control method that adaptively searches for the optimum value of the evaluation amount based on the change in the evaluation amount according to the change in the manipulated variable.
  • the evaluation amount is a value serving as an optimization index for a process to be controlled (hereinafter referred to as “control target process”).
  • the evaluation amount is an index value determined based on the control amount of the process to be controlled, and the relationship between the evaluation amount and the control amount is represented by a predetermined evaluation function.
  • This evaluation function may be set based on an arbitrary evaluation criterion as long as it is based on the control amount.
  • the evaluation amount may be the control amount itself.
  • this evaluation function is an unknown function with respect to the manipulated variable.
  • the manipulated variable is changed by a periodic signal called a dither signal. Usually, this dither signal is often given as a sine wave.
  • the manipulated variable is continuously vibrated by a dither signal, and the change (increase / decrease) in the evaluation value is observed. Then, based on the change in the evaluation amount with respect to the change in the operation amount, the operation amount is changed in such a direction that the evaluation amount approaches the optimum value of the evaluation function. The optimum value of the evaluation function is searched by repeating such changes in the operation amount.
  • 1A is an unknown evaluation function with respect to the manipulated variable.
  • an unknown evaluation function is assumed as a downward convex quadratic function.
  • the evaluated quantity changes as shown in FIG.
  • the evaluation amount decreases with respect to the increase in the operation amount, it can be seen that the operating point has changed on the left side of the minimum value 101 of the evaluation function curve 100.
  • the evaluation amount is changed as shown in FIG. In this case, since the evaluation amount increases with respect to the increase in the operation amount, it can be seen that the operating point has changed on the right side from the minimum value 101.
  • PID control Proportional-Integral-Derivative Control
  • PID control Proportional-Integral-Derivative Control
  • extreme value control is an optimal value search type control method that searches for an operation amount that optimizes the evaluation amount, and therefore represents the relationship between the operation amount and the control amount as in PID control.
  • No process model is required in advance. Therefore, extreme value control is an effective control method even for a process to be controlled in which a target value cannot be set in advance, and has the potential to be widely used in the future.
  • An extreme value controller that performs extreme value control based on this principle can be realized with a relatively simple configuration.
  • FIG. 2 is a block diagram showing a configuration example of the extreme value controller.
  • 2 includes a high-pass filter 21 (LPF: Low-Pass Filter), a dither signal output unit 22, a low-pass filter 23 (HPF: High-Pass Filter), and an integrator 24.
  • LPF Low-Pass Filter
  • HPF High-Pass Filter
  • integrator 24 an integrator 24.
  • the configuration of the extreme value controller 2 is as complex as the conventional PID controller. Therefore, the extreme value controller 2 can be easily mounted using hardware such as a PLC (Programmable Logic Controller), like the PID controller.
  • PLC Programmable Logic Controller
  • the extreme value controller 2 forcibly changes the operation amount of the process 200 to be controlled by applying a dither signal (M) having a periodic change.
  • M dither signal
  • modulation the operation amount of the control target process 200 periodically changes.
  • the control target process 200 outputs a control amount in accordance with the input of the modulated operation amount.
  • the control target process 200 acquires and outputs the evaluation amount that has changed in accordance with the change in the control amount.
  • the evaluation amount output from the control target process 200 is fed back to the extreme value controller 2.
  • the change (response) of the evaluation amount with respect to the change of the operation amount appears with a certain time delay.
  • extreme value control is a control method for searching for an extreme value of an unknown evaluation function with respect to an operation amount. Therefore, the evaluation function of the control target process 200 is premised on having a minimum value, but the value is unknown.
  • the high-pass filter 21 removes a constant bias corresponding to the unknown minimum value from the feedback evaluation amount. In other words, this process is a process for always adjusting the unknown local minimum value to zero, and is a preprocess necessary for determining the direction of change (increase or decrease) that the integrator 24 gives to the manipulated variable. .
  • the dither signal output unit 22 causes the dither signal (D) to act on the evaluation amount adjusted in this way. Thereby, the same frequency component as that of the dither signal (M) is extracted from the evaluation amount changed by the modulation.
  • this operation is referred to as demodulation (demodulation).
  • demodulation demodulation
  • the role of demodulation is as follows.
  • the evaluation function for the operation amount of the control target process 200 is unknown. Therefore, the evaluation function may include a non-linear element.
  • the evaluation function is assumed to be a non-linear function convex downward (convex upward in the case of local maximum search). Due to such a non-linear element, there is a high possibility that a harmonic component or a subharmonic component corresponding to the frequency ⁇ of the dither signal (M) appears in the evaluation amount.
  • Demodulation is a process for removing the influence of such harmonics and subharmonics. By this demodulation, a component having the same frequency ⁇ as that of the dither signal (M) in which the evaluation amount is changed is extracted from the components included in the evaluation amount.
  • the demodulated evaluation amount is input to the low-pass filter 23.
  • the low-pass filter 23 extracts a steady component (low frequency component) from the evaluation amount.
  • the steady component is considered to indicate whether the evaluation amount has changed in the increasing direction or the decreasing direction due to the application of the dither signal (M).
  • the integrator 24 functions as an estimator that estimates the direction (gradient) of the operation amount to be moved to bring the evaluation amount close to the minimum value by integrating the steady component extracted by the low-pass filter 23.
  • This kind of gradient estimation is based on the most basic gradient method as an adaptive control system estimation method.
  • the dither signal (M) is adjusted by the gradient estimated by the integrator 24 so as to give a change that brings the evaluation amount close to the minimum value with respect to the operation amount.
  • the maximum value search can be realized, for example, by inverting the sign of the integrator 24.
  • FIG. 3 is a functional block diagram illustrating a functional configuration of the optimum control apparatus according to the embodiment.
  • the plant 300 and the evaluation amount acquisition unit 400 correspond to the control target process 200 of FIG.
  • the evaluation amount acquisition unit 400 acquires the evaluation amount J based on the control amount Y output from the plant 300.
  • the optimal control device 1 performs extreme value control of the plant 300 based on the evaluation amount J output from the evaluation amount acquisition unit 400.
  • the optimal control device 1 includes a CPU (Central Processing Unit) connected via a bus, a memory, an auxiliary storage device, and the like, and executes an optimal control program.
  • the optimal control device 1 functions as a device including a control target parameter determination unit 11, an extreme value control parameter adjustment unit 12, and an extreme value control controller 13 by executing an optimal control program. All or some of the functions of the optimal control device 1 may be realized by using hardware such as ASIC (Application Specific Integrated Circuit), PLD (Programmable Logic Device), and FPGA (Field Programmable Gate Array). .
  • the optimal control program may be recorded on a computer-readable recording medium.
  • the computer-readable recording medium is, for example, a portable medium such as a flexible disk, a magneto-optical disk, a ROM, a CD-ROM, or a storage device such as a hard disk built in the computer system.
  • the optimal control program may be transmitted via a telecommunication line.
  • the control target parameter determination unit 11 determines a control target parameter indicating the characteristics of the plant 300.
  • the control target parameter is a parameter such as a time constant or dead time of the plant 300.
  • the control target parameter determination unit 11 determines the control target parameter based on the operation amount U input to the plant 300 and the evaluation amount J output from the evaluation amount acquisition unit 400.
  • the extreme value control parameter adjustment unit 12 determines an extreme value control parameter for executing extreme value control on the plant 300.
  • the extreme value control parameters are parameters such as the frequency of the low-pass filter and the high-pass filter, the frequency and amplitude of the dither signal, and the gain of the integrator.
  • the extreme value control parameter adjustment unit 12 determines an extreme value control parameter based on the control target parameter determined by the control target parameter determination unit 11.
  • the extreme value controller 13 executes extreme value control of the process to be controlled based on the extreme value control parameter determined by the extreme value control parameter adjustment unit 12. Details of the configuration of the extreme value controller 13 are the same as those of the extreme value controller 2 shown in FIG.
  • the process to be controlled that is the control target of the optimal control device 1 is not limited to the plant.
  • the controlled process may be an arbitrary process having an evaluation amount to be optimized.
  • the function of the optimal control apparatus 1 will be described in detail by taking a water treatment plant that realizes a biological wastewater treatment process as an example of the plant 300 as an example.
  • FIG. 4 is a diagram showing an outline of a water treatment plant.
  • the water treatment plant 500 of FIG. 4 includes an anaerobic tank 510, an oxygen-free tank 520, an aerobic tank 530, and a final sedimentation basin 540.
  • Anaerobic tank 510 is equipment for activating microorganisms.
  • the anoxic tank 520 is equipment for removing nitrogen.
  • the aerobic tank 530 is equipment for decomposing organic substances, removing phosphorus, and nitrifying ammonia.
  • the final sedimentation basin 540 is equipment for precipitating activated sludge.
  • the water treatment plant 500 is equipped with facilities such as a pump for conveying water and sludge between the above facilities, a blower for supplying air into the tank, and a sensor for measuring the concentration of substances in the air or water.
  • the chemical injection pump 511 is a pump that supplies chemicals such as a carbon source that activates microorganisms to the anaerobic tank 510.
  • the circulation pump 531 is a pump that controls the circulation amount of the water to be treated that circulates between the aerobic tank 530 and the anoxic tank 520.
  • the blower 532 controls the amount of aeration by supplying air to the aerobic tank 530.
  • the return sludge pump 541 is a pump that returns the sludge from the final sedimentation basin 540 to the anoxic tank 520.
  • the excess sludge extraction pump 542 is a pump that extracts excess sludge from the final sedimentation basin 540.
  • the sensor 512 and the sensor 543 measure the quality of discharged water in the anaerobic tank 510 and the final sedimentation basin 540, respectively.
  • the manipulated variable is the return rate of the returned sludge
  • the controlled variable is the concentration of nitrogen and phosphorus contained in the discharged water (hereinafter referred to as “discharged nitrogen concentration” and “released phosphorus, respectively”). Called “concentration”).
  • the return rate is obtained by dividing the discharge flow rate of the return sludge pump 541 by the inflow amount.
  • the released nitrogen concentration and the released phosphorus concentration are acquired by the sensor 512 and the sensor 543.
  • the control amount may be the amount of nitrogen and phosphorus contained in the discharge water (hereinafter referred to as “discharge nitrogen amount” and “discharge phosphorus amount”, respectively).
  • the released nitrogen amount and the released phosphorus amount can be obtained by multiplying the released nitrogen concentration and the released phosphorus concentration by the released flow rate, respectively.
  • an evaluation function for acquiring the evaluation amount based on the control amount output from the water treatment plant 500 is set in advance.
  • the evaluation function here defines an unknown evaluation function for the manipulated variable as a function of the control amount.
  • the evaluation function is a function representing the relationship between the released nitrogen concentration and released phosphorus concentration and the evaluation amount. This evaluation function needs to be set so as to take an extreme value between the control amount at the upper limit of the operation amount (return rate) and the control amount at the lower limit of the operation amount.
  • the evaluation amount is expressed as the sum of the water quality cost based on the concept of drainage levy and the power cost of the return sludge pump 541 (hereinafter referred to as “total cost”). Can be considered.
  • the power cost of the return sludge pump 541 can be calculated from the return sludge flow rate, the rated power of the return sludge pump 541, and the like.
  • water quality cost is expressed by the following equation.
  • COD means chemical oxygen demand
  • BOD means biochemical oxygen demand
  • TN means released nitrogen
  • TP means released phosphorus.
  • the conversion factor for each cost may be determined based on the actual drainage levy, or may be determined by other methods. In general, it is known that, among COD, BOD, TN, and TP, those that greatly change by changing the return rate are TN and TP. Therefore, here, the water quality cost is expressed by the following equation (2).
  • the evaluation function may be set based only on the water quality cost.
  • the evaluation function is expressed as the total cost including the operation cost (electric power cost). (Return rate) It is set to take an extreme value between the control amount at the upper limit and the control amount at the lower limit of the operation amount.
  • a function that directly represents the evaluation of water quality may be set as the evaluation function.
  • the evaluation amount may be calculated as in the following formula (3).
  • TN lim and TP lim are parameters representing a threshold level corresponding to the regulation value and management value of the discharged water quality.
  • the evaluation amount rapidly increases when the threshold level is exceeded. Therefore, it can be expected that the extreme value control functions so as to suppress the evaluation amount within the threshold level.
  • the evaluation function setting method necessary for extreme value control has been described using the water treatment plant 500 as shown in FIG. 4 as an example.
  • setting of the evaluation function is not necessary.
  • An example of this is the control of wind turbine blades in a wind power plant.
  • the evaluation amount is the power generation amount
  • the operation amount is the rotation angle of the windmill blade.
  • the evaluation amount acquisition unit 400 may not be provided.
  • extreme value control may be applicable by acquiring the evaluation amount.
  • the control target parameter determination unit 11 determines the control target parameter based on the control amount and the evaluation amount acquired as described above. Hereinafter, control target parameters determined by the control target parameter determination unit 11 will be described.
  • time constant One of the control target parameters to be determined by the control target parameter determination unit 11 is a time constant of the control target process.
  • This time constant is normally used as a control parameter for PID control.
  • the time constant may be simply set based on the above definition, but the evaluation amount includes a plurality of control amounts as in the water treatment plant 500 described above. When acquired based on (TN concentration and TP concentration), the larger time constant (that is, the slower response speed) is set as a representative value.
  • the time constant may be identified by actually changing the manipulated variable and measuring the response time, it can be identified using any open-loop (open-loop) or closed-loop (closed-loop) system identification method. Also good.
  • One of the control target parameters to be determined by the control target parameter determination unit 11 is a dead time (until the response of the control amount is obtained from the input of the operation amount or until the evaluation amount is acquired from the input of the operation amount ( Delay time).
  • the dead time can be identified by any identification method, but the simplest method is to collect time-series data for manipulated variables and time-series data for controlled variables and calculate the correlation coefficient while shifting the time. It is. In this case, the time difference until the value of the correlation coefficient becomes maximum is the dead time.
  • the value of the second order differential value G is not always constant. Therefore, in such a case, a representative value such as an average value or a median value is used, or a second-order differential value G is determined for each of several operating points.
  • parameter identification is performed using the process model represented by the following equations (4) and (5). It is also possible.
  • Equation (4) is a model representing “first-order lag + dead time process” that is also commonly used in PID control.
  • Expression (5) represents the evaluation function of the model represented by Expression (4) as a quadratic function having the simplest convex shape.
  • Y (t) in Expression (4) represents a control amount at time t.
  • u represents an operation amount.
  • Ku represents the process gain.
  • T represents a time constant, and s represents a Laplace operator.
  • L represents the dead time.
  • J (t) in Equation (5) represents the evaluation amount at time t.
  • K y represents the gain of the evaluation function.
  • K arg represents the control amount y giving the optimum value
  • K min represents the optimum value (minimum value) of the evaluation function.
  • K arg and K min in equation (5) are often not known.
  • the process model of Expression (4) can be identified using a normal identification method if the control amount can be measured.
  • the K y of formula (5) is calculated directly from the relational expression between the control amount and the evaluation value, or be estimated by measuring the variation of the evaluation value J when changing a control amount it can. If K u and K y can be estimated in this way, the second-order differential value G (in the steady state) can be calculated by a simple calculation formula as shown in the following formula (6). it can.
  • One of the control target parameters to be determined by the control target parameter determination unit 11 is an upper limit value and a lower limit value of the operation amount.
  • the upper and lower limit values of the operation amount implemented in the plant monitoring control system or the like may be used as they are, or the maximum value and the minimum value in the past operation data of the plant are respectively set as upper limit values. And it may be used as a lower limit.
  • control target parameter determination unit 11 may automatically determine the control target parameter using the plant measurement data, or input a parameter identified by performing a necessary test such as a step response test. It may be configured to accept. Moreover, when identification using plant data is difficult, it may be configured to accept input of assumed values for some or all of the control target parameters.
  • the extreme value control parameter adjustment unit 12 determines an extreme value control parameter based on the control target parameter determined by the control target parameter determination unit 11.
  • the extreme value control parameters determined by the extreme value control parameter adjusting unit 12 will be described.
  • Extreme control parameter adjustment unit 12 calculates the fluctuation range U R of the operation amount from the upper limit value and the lower limit value of the manipulated variable.
  • the extreme control parameter adjusting section 12 with respect to the variation range U R of the manipulated variables, in advance set the parameters k 1 indicating whether to allow the oscillating at how much variation width dither signal.
  • the extreme value control parameter adjusting unit 12 determines the amplitude a of the dither signal based on the set parameter k 1 as shown in the following equation (7).
  • the parameter k 1 is the ratio of the amplitude a to the variation range U R of the manipulated variables.
  • the parameter k 1 may be set to about 0.01 to 0.1.
  • the amplitude a of the dither signal is determined based on the parameter k 1 .
  • the amplitude a of the dither signal is a signal that is forcibly added for modulation of the manipulated variable, and the signal fed back is input to the low-pass filter. It is assumed that the signal output from the low-pass filter hardly oscillates, and the amplitude of the steady operation amount vibration is assumed to be approximately the same as the amplitude a of the dither signal.
  • the extreme value control parameter adjusting unit 12 does not adversely affect the normal control of the process to be controlled, and the extreme value The amplitude can be determined so as to produce a change in the evaluation quantity that is necessary for control.
  • the value of the parameter k 1 it is preferable that the default value is set.
  • the extreme value control parameter adjustment unit 12 determines the frequency of the dither signal based on two control target parameters, a time constant and a dead time. First, the extreme value control parameter adjustment unit 12 acquires a time constant T mod as a comprehensive time constant of the process to be controlled.
  • the time constant T mod is calculated by the following equation (8).
  • T mod in Expression (8) corresponds to a time constant in the case where “first order delay + dead time process” expressed by Expression (4) is approximated by Padé.
  • the extreme value control parameter adjustment unit 12 acquires the bandwidth of the process to be controlled based on T mod acquired by Expression (8).
  • the bandwidth ⁇ b of the process to be controlled is calculated by the following equation (9).
  • the extreme value control parameter adjusting unit 12 determines the frequency of the dither signal based on the bandwidth ⁇ b of the process to be controlled acquired by Expression (9).
  • the frequency ⁇ of the dither signal is calculated by the following equation (10).
  • k 2 is a parameter that takes values in the range of 5-10.
  • this setting makes it possible to regard the process to be controlled as a static process that is stationary as seen from the time scale of the scanning signal for extreme value search (that is, the dither signal).
  • k 2 has a default value as in k 1 .
  • the dynamics of the controlled process and the time scale of the scanning of the extremum search set 10 so as to be completely separated as the default value for k 2.
  • the extreme value control parameter adjusting unit 12 determines the extreme value control parameters related to the low-pass filter and the high-pass filter using the parameters determined above.
  • a method for determining control parameters related to the low-pass filter and the high-pass filter will be described.
  • the extreme value control parameter adjustment unit 12 determines the frequency ⁇ 1 of the low-pass filter by the following equation (11).
  • the extreme value controller 13 may be configured not to include a low-pass filter.
  • ⁇ 1 may be set to ⁇ (infinity).
  • the high-pass filter is used for the purpose of removing the bias component in order to make the minimum value of the evaluation function zero.
  • the high pass filter needs to pass the dither signal. Therefore, the extreme value control parameter adjusting unit 12 determines the frequency ⁇ 2 of the high-pass filter by the following equation (12).
  • ⁇ 2 may be an arbitrary value as long as it is sufficiently smaller than 0.1 ⁇ ⁇ .
  • ⁇ 2 may be a sufficiently small fluctuation value that satisfies ⁇ 2 ⁇ 0.1 ⁇ ⁇ .
  • the extreme value control parameter adjustment unit 12 determines the gain of the integrator (hereinafter referred to as “integral gain”) based on the control target parameter and the extreme value control parameter determined as described above.
  • the integral gain KI is expressed by the following equations (13) and (14) using the frequency ⁇ of the dither signal.
  • k 3 is preferably set to 5 to 10.
  • the above KI 0 is a factor that greatly affects the performance of extreme value control.
  • a method for setting the KI 0 will be described.
  • Equation (14) is a mathematical formula derived based on an average system (average system) used for stability analysis of the extreme value control system.
  • An average system is a system in which when a periodic input is applied to a certain system, the behavior (output) of the system can be expressed by an average value (average) in that period.
  • the process to be controlled is a static process having no dynamics
  • the average system of the extreme value control system is expressed by the following equation (15).
  • D J represents a gradient related to the periodic average (xx * ) of the input of the evaluation function J.
  • x * is the equilibrium point of x.
  • ⁇ is a time function scaled by the frequency ⁇ of the dither signal. It is a value represented by the following formula (16).
  • the average system of Equation (15) represents the dynamics related to the convergence of extreme value control. Specifically, it represents how fast the evaluation amount converges to the minimum value (minimum value) with respect to the operation amount given periodic vibration by the dither signal.
  • the slope of D J of formula (15) can be expressed from equation (5) as in the following equation (17).
  • the time constant of the average system indicated by Equation (19) is the time constant on the time axis ⁇ .
  • FIG. 5 summarizes the method for determining the control target parameter, the method for adjusting the extreme value control parameter, and the conventional design guideline in the embodiment described above.
  • the optimal control apparatus 1 of the embodiment configured as described above is capable of determining an extreme value control parameter and adjusting an extreme value control parameter by only setting five control object parameters. With part 12. Therefore, by using the optimum control device 1 of the embodiment, it is possible to realize extreme value control of the process to be controlled on a scale that can be implemented in the PLC.
  • the above-described optimal control device 1 may be configured to acquire process data during operation and automatically update the extreme value control parameters. With this configuration, the optimal control device 1 can perform control with extreme value control parameters corresponding to the state of the process to be controlled that changes over time, and improve the performance of searching for the optimum value by extreme value control. Can be made.
  • the optimal control device 1 displays the prediction information indicating the prediction of how the behavior of the control target process changes by setting the control target parameter and the extreme value control parameter determined by the above method. , Not shown).
  • FIG. 6 is a diagram illustrating a specific example of display of prediction information.
  • the display screen 600 includes a mode display area 610, a control target parameter display area 620, a pre-adjustment extreme value control parameter display area 630, a post-adjustment extreme value control parameter display area 640, a pre-adjustment control information display area 650, and a post-adjustment control.
  • An information display area 660 is provided.
  • the mode display area 610 is an area in which an extreme value control execution mode is displayed.
  • the example of FIG. 6 shows a case where manual control is selected from two modes of automatic control and manual control.
  • the control target parameter display area 620 is an area in which the value of the control target parameter determined by the optimal control device 1 is displayed.
  • the pre-adjustment extreme value control parameter display area 630 is an area in which the value of the current extreme value control parameter before change is displayed.
  • the adjusted extreme value control parameter display area 640 is an area in which the value of the newly determined extreme value control parameter is displayed.
  • the pre-adjustment control information display area 650 is an area in which a control result with the current extreme value control parameter is displayed.
  • the post-adjustment control information display area 660 is an area in which a prediction of a control result when extreme value control is performed using a newly determined extreme value control parameter is displayed.
  • the optimal control device 1 includes a prediction unit (not shown) that simulates the behavior of the control target process based on the process model of the control target process, and a display control unit (not shown) that generates the display screen.
  • the display control unit generates a display screen that displays the value of the extreme value control parameter before and after the change by combining the prediction result simulated by the prediction unit and the control result of the current extreme value control parameter.
  • the optimal control device 1 stores a plurality of parameter sets of the control target parameters and the extreme value control parameters, and uses a predetermined switching criterion. It may be configured as a gain scheduling type control device that switches the extreme value control parameter accordingly.
  • the optimal control apparatus 1 may be configured to periodically identify the control target parameter and update the extreme value control parameter when the difference from the current control target parameter value exceeds a predetermined threshold value.
  • the threshold value for determining the difference between the parameters to be controlled may be set simply based on the absolute value of the difference, or the ratio of the change based on the parameter after the change ((parameter value before adjustment ⁇ parameter after adjustment) (Value) ⁇ adjusted parameter value ⁇ 100 (%)). In the latter case, for example, the threshold value may be ⁇ (%), and the extreme value control parameter may be updated when the rate of change becomes ⁇ ⁇ % or more.
  • the prediction unit of the first modification is provided, the extreme value is obtained when the similarity between the actual control result and the prediction result is evaluated by the sum of squares of errors and the response speed is assumed to change greatly.
  • the control parameter may be updated.
  • the optimum control device 1 updates the extreme value control parameter as a device that executes extreme value control. However, the optimum control device 1 determines whether the current control state is valid based on the criteria for updating the parameter. A function as a diagnostic device for diagnosing sex may be provided. In this case, for example, when the update of the extreme value control parameter is recommended, the optimal control device 1 may include a notification unit that notifies an operator such as a plant to that effect.
  • control target parameter determination unit that determines the control target parameter of the control target process based on the operation amount and the evaluation amount, and the control target determined by the control target parameter determination unit And an extreme value control parameter determination unit for determining an extreme value control parameter for executing extreme value control based on the parameter, thereby enabling extreme values with control parameters appropriately set according to the characteristics of the process to be controlled. Value control can be performed.

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Organic Chemistry (AREA)
  • Hydrology & Water Resources (AREA)
  • Microbiology (AREA)
  • Chemical & Material Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Environmental & Geological Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Artificial Intelligence (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Software Systems (AREA)
  • Medical Informatics (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Feedback Control In General (AREA)
  • Purification Treatments By Anaerobic Or Anaerobic And Aerobic Bacteria Or Animals (AREA)

Abstract

Selon un mode de réalisation, la présente invention concerne un système de contrôle optimal comportant une unité de détermination de paramètres contrôlés, une unité de détermination de paramètres de contrôle de valeur extrême, une unité de contrôle de valeur extrême. L'unité de détermination de paramètres contrôlés détermine un paramètre contrôlé indiquant une caractéristique d'un processus à contrôler, cette détermination étant effectuée sur la base d'une quantité de fonctionnement et d'une quantité d'évaluation représentant un indicateur concernant l'optimisation du processus à contrôler sur la base d'une quantité de contrôle qui varie en fonction de la quantité de fonctionnement. L'unité de détermination de paramètre de contrôle de valeur extrême détermine un paramètre de valeur extrême sur la base du paramètre contrôlé déterminé par l'unité de détermination de paramètres contrôlés. L'unité de contrôle de valeur extrême exécute le contrôle de valeur extrême au moyen du paramètre de contrôle de valeur extrême déterminé par l'unité de détermination de paramètre de contrôle de valeur extrême.
PCT/JP2016/070821 2015-07-29 2016-07-14 Dispositif de contrôle optimal, procédé de contrôle optimal, programme informatique et système de contrôle optimal WO2017018227A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201680044271.4A CN107924162B (zh) 2015-07-29 2016-07-14 最优控制装置、最优控制方法、记录介质及最优控制系统

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2015-149830 2015-07-29
JP2015149830A JP6523854B2 (ja) 2015-07-29 2015-07-29 最適制御装置、最適制御方法、コンピュータプログラム及び最適制御システム

Publications (1)

Publication Number Publication Date
WO2017018227A1 true WO2017018227A1 (fr) 2017-02-02

Family

ID=57885542

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2016/070821 WO2017018227A1 (fr) 2015-07-29 2016-07-14 Dispositif de contrôle optimal, procédé de contrôle optimal, programme informatique et système de contrôle optimal

Country Status (3)

Country Link
JP (1) JP6523854B2 (fr)
CN (1) CN107924162B (fr)
WO (1) WO2017018227A1 (fr)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6744145B2 (ja) * 2016-06-15 2020-08-19 株式会社東芝 制御装置、制御方法及びコンピュータプログラム
JP6763831B2 (ja) * 2017-07-03 2020-09-30 横河電機株式会社 制御システム及び制御方法
JP7154774B2 (ja) * 2018-02-27 2022-10-18 株式会社東芝 最適制御装置、制御方法及びコンピュータプログラム
JP7082900B2 (ja) * 2018-04-24 2022-06-09 三菱重工業株式会社 制御装置及び燃焼システム
JP7132016B2 (ja) * 2018-07-27 2022-09-06 株式会社東芝 制御装置、制御方法及びコンピュータプログラム
JP7267779B2 (ja) * 2019-03-04 2023-05-02 株式会社東芝 最適制御装置、最適制御方法及びコンピュータプログラム
WO2020241657A1 (fr) * 2019-05-29 2020-12-03 東芝インフラシステムズ株式会社 Dispositif de commande optimale, procédé de commande optimale et programme informatique

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09325801A (ja) * 1996-06-06 1997-12-16 Nitto Kogyo Kk 極値制御回路
US6098010A (en) * 1997-11-20 2000-08-01 The Regents Of The University Of California Method and apparatus for predicting and stabilizing compressor stall
JP2015513706A (ja) * 2012-05-10 2015-05-14 三菱電機株式会社 システムの動作を制御するコントローラー及び方法

Family Cites Families (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01276304A (ja) * 1988-04-28 1989-11-06 Toshiba Corp 最適化制御装置
US5404423A (en) * 1992-06-16 1995-04-04 Nippon Telegraph And Telephone Corporation Method and apparatus for indetification, forecast, and control of a non-linear flow on a physical system network using a neural network
FR2700080B1 (fr) * 1992-12-30 1995-01-27 Unite Hermetique Sa Alimentation optimale d'un moteur électrique.
US6627309B2 (en) * 2001-05-08 2003-09-30 3M Innovative Properties Company Adhesive detackification
JP4309326B2 (ja) * 2004-10-06 2009-08-05 本田技研工業株式会社 プラントの制御装置
JP2007287063A (ja) * 2006-04-20 2007-11-01 Hitachi Ltd 最適制御方法、最適制御システム、統括制御装置およびローカル制御装置
JP4639166B2 (ja) * 2006-05-18 2011-02-23 本田技研工業株式会社 制御装置
GB2463827B (en) * 2007-07-17 2012-09-05 Johnson Controls Tech Co Extremum seeking control with actuator saturation control
GB2463218B (en) * 2007-07-17 2012-12-05 Johnson Controls Tech Co Extremum seeking control with reset control
WO2010004611A1 (fr) * 2008-07-07 2010-01-14 本田技研工業株式会社 Régulateur
CN101408752B (zh) * 2008-10-21 2014-03-26 中国人民解放军海军航空工程学院 基于混沌退火和参数扰动的神经网络极值控制方法及系统
CN101825867A (zh) * 2009-03-06 2010-09-08 株式会社东芝 设备最优化运行系统,最优运转点计算方法以及最优运转点计算程序
US8412357B2 (en) * 2010-05-10 2013-04-02 Johnson Controls Technology Company Process control systems and methods having learning features
CN101859106B (zh) * 2010-06-23 2011-12-14 浙江大学 一种发酵生产过程控制方法及应用
CN102176120B (zh) * 2011-02-24 2013-08-14 同济大学 随机地震激励系统磁流变阻尼最优控制的方法
CN103259488B (zh) * 2012-02-17 2016-04-27 通用电气公司 电机控制方法、控制系统及控制装置
CN102880182B (zh) * 2012-09-12 2015-01-14 北京航空航天大学 一种存在网络随机延迟的微小型无人飞行器控制方法
JP6216112B2 (ja) * 2012-11-30 2017-10-18 アズビル株式会社 多変数制御装置および方法
US9448546B2 (en) * 2013-03-15 2016-09-20 Rockwell Automation Technologies, Inc. Deterministic optimization based control system and method for linear and non-linear systems
CN103365213B (zh) * 2013-07-15 2015-10-21 温州大学 用于兆瓦级逆变系统的极值优化自整定数字pid控制方法
CN103926947B (zh) * 2014-04-16 2017-03-22 安阳师范学院 大跨度索桥结构非线性系统半主动振动控制方法
CN104122795A (zh) * 2014-07-15 2014-10-29 河海大学常州校区 基于新型极值函数指标的智能自整定pid室温控制算法
CN104361204A (zh) * 2014-10-20 2015-02-18 上海电机学院 利用自组织极值优化处理进行控制优化的方法
CN104554251A (zh) * 2014-12-09 2015-04-29 河南理工大学 基于道路坡度信息的混合动力汽车节能预测控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09325801A (ja) * 1996-06-06 1997-12-16 Nitto Kogyo Kk 極値制御回路
US6098010A (en) * 1997-11-20 2000-08-01 The Regents Of The University Of California Method and apparatus for predicting and stabilizing compressor stall
JP2015513706A (ja) * 2012-05-10 2015-05-14 三菱電機株式会社 システムの動作を制御するコントローラー及び方法

Also Published As

Publication number Publication date
CN107924162A (zh) 2018-04-17
JP2017033104A (ja) 2017-02-09
CN107924162B (zh) 2021-03-05
JP6523854B2 (ja) 2019-06-05

Similar Documents

Publication Publication Date Title
WO2017018227A1 (fr) Dispositif de contrôle optimal, procédé de contrôle optimal, programme informatique et système de contrôle optimal
Yang et al. Active disturbance rejection control on first-order plant
JP6147687B2 (ja) 制御パラメータ調整装置、および制御パラメータ調整プログラム
Roy et al. Level control of two tank system by fractional order integral state feedback controller tuned by PSO with experimental validation
JP2019027212A (ja) 水道の維持管理支援装置及び水道の維持管理支援システム
Jämsä-Jounela et al. Evaluation of control performance: methods, monitoring tool and applications in a flotation plant
WO2020241657A1 (fr) Dispositif de commande optimale, procédé de commande optimale et programme informatique
KR101125318B1 (ko) 수처리공정에 이용되는 오존주입률 결정장치 및 방법
Adgar et al. Enhancement of coagulation control using the streaming current detector
CN107250929B (zh) 控制系统、控制装置、控制方法以及非暂时性存储介质
WO2019167998A1 (fr) Dispositif de commande optimale, procédé de commande et programme informatique
JP6744145B2 (ja) 制御装置、制御方法及びコンピュータプログラム
JP2021186761A (ja) 学習モデル生成装置、推測装置および散気量制御装置
RU2368934C2 (ru) Адаптивная система управления
Gerov et al. Parameter estimation method for the unstable time delay process
JP7132016B2 (ja) 制御装置、制御方法及びコンピュータプログラム
di Capaci et al. Identification and estimation of valve stiction by the use of a smoothed model
RU2339988C1 (ru) Адаптивная система управления
RU2343524C1 (ru) Адаптивная система управления
JP2024017830A (ja) 最適制御装置、最適制御方法およびコンピュータプログラム
RU2612084C1 (ru) Способ адаптивного управления плохо формализуемым объектом
Aquino et al. An emotional controller plc implementation for an industrial fan system
Vrána et al. PID contoller autotuning based on nonlinear tuning rules
Jiao et al. Research on the weak signal detection method based on adaptive vibrational resonance
Ucun et al. Towards optimizing canal system operations: Extremum-seeking controller design via a frequency-based control approach

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16830332

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 22.06.2018)

122 Ep: pct application non-entry in european phase

Ref document number: 16830332

Country of ref document: EP

Kind code of ref document: A1