WO2021214839A1 - 制御装置、制御システム、制御方法及びプログラム - Google Patents

制御装置、制御システム、制御方法及びプログラム Download PDF

Info

Publication number
WO2021214839A1
WO2021214839A1 PCT/JP2020/017093 JP2020017093W WO2021214839A1 WO 2021214839 A1 WO2021214839 A1 WO 2021214839A1 JP 2020017093 W JP2020017093 W JP 2020017093W WO 2021214839 A1 WO2021214839 A1 WO 2021214839A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
frequency
control
frequency component
control target
Prior art date
Application number
PCT/JP2020/017093
Other languages
English (en)
French (fr)
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
Application filed by 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to JP2020555933A priority Critical patent/JP6833133B1/ja
Priority to PCT/JP2020/017093 priority patent/WO2021214839A1/ja
Publication of WO2021214839A1 publication Critical patent/WO2021214839A1/ja

Links

Images

Classifications

    • 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

  • This disclosure relates to control devices, control systems, control methods and programs.
  • a control device having a tuning function that automatically determines the value of an initial control parameter or adjusts the value of a control parameter from the data under actual control is known. There is.
  • Patent Document 1 uses it for PID (Proportional Integral Differential) control from a measured value which is a control response without adding a special tuning signal to a controlled object.
  • PID controller that can be automatically changed and set by obtaining a PID constant is disclosed.
  • the PID controller described in Patent Document 1 calculates a first PID constant using a known step response method based on the target value and the measured value when the deviation between the target value and the measured value exceeds the limit value.
  • the second PID constant is calculated using a known expert method.
  • the PID controller described in Patent Document 1 changes an already set PID constant to any PID constant or uses these PID constants based on the evaluation of the degree of certainty of these PID constants. Change to a third PID constant based on it.
  • Patent Document 2 discloses a plant control device that automatically adjusts the value of a control parameter of the controller according to a mathematical model representing a dynamic characteristic of a controlled object.
  • the plant control device described in Patent Document 2 represents an input signal representing an operation amount to the control target when an identification signal generated using known information of the control target is input to the controller, and a control amount of the control target.
  • a mathematical model is calculated using a known least-squares method based on the output signal.
  • the PID controller described in Patent Document 1 when the deviation between the target value and the measured value exceeds the limit value, it is considered that the step signal, which is a dedicated input signal in the normal step response method, is input, and the first PID The constant is calculated. Therefore, the first PID constant cannot be calculated unless the deviation exceeds the limit value, and even if the deviation can be calculated exceeding the limit value, the degree of certainty is lower than that calculated by the normal step response method. There is. Further, the PID controller described in Patent Document 1 calculates the second PID constant by using the expert method, but the expert method calculates the PID constant based on the inference based on the operator's empirical rule.
  • the PID controller described in Patent Document 1 determines the PID constant to be changed based on the calculated evaluation of the degree of certainty of these PID constants, but the PID constant determined from the above-mentioned problem The accuracy itself may be low.
  • a mathematical model representing the dynamic characteristics of the controlled object is indicated by a differential equation. Therefore, when trying to calculate a mathematical model by the minimum square method based on the input signal and output signal to be controlled, it is necessary to perform differential operations on these signals, but when performing differential operations on signals containing noise, There is a problem of amplifying the high frequency component of noise. Therefore, the plant control device described in Patent Document 2 has a problem that the mathematical model cannot be calculated accurately.
  • the present disclosure has been made in view of the above circumstances, and an object of the present disclosure is to accurately calculate the values of control parameters.
  • the control device includes a control target control unit, a frequency component extraction unit, a frequency characteristic calculation unit, a model calculation unit, and a control parameter change unit.
  • the control target control unit calculates the operation amount of the control target using the value of the preset control parameter based on the deviation between the target amount and the control amount of the control target, and outputs the signal of the calculated operation amount. Output to the control target to control the control target.
  • the frequency component extraction unit extracts the signal of the first frequency component and the signal of the second frequency component from the signal of the manipulated amount, and extracts the signal of the first frequency component and the signal of the second frequency component from the signal of the controlled amount. ..
  • the frequency characteristic calculation unit is a control target including a phase delay and a gain based on the signal of the first frequency component extracted from the signal of the operation amount and the signal of the first frequency component extracted from the signal of the control amount by the frequency component extraction unit.
  • the phase delay and gain based on the signal of the second frequency component extracted from the signal of the manipulated amount and the signal of the second frequency component extracted from the signal of the controlled amount by the frequency component extractor while calculating the characteristics of the first frequency of The characteristics of the second frequency of the controlled object including the above are calculated.
  • the model calculation unit calculates the model of the control target based on the characteristics of the first frequency and the characteristics of the second frequency of the control target calculated by the frequency characteristic calculation unit.
  • the control parameter changing unit calculates the value of the control parameter based on the model to be controlled calculated by the model calculation unit and changes the setting.
  • the first frequency of the controlled object including the phase delay and the gain based on the signal of the first frequency component extracted from the signal of the manipulated amount and the signal of the first frequency component extracted from the signal of the controlled amount.
  • the characteristics of the second frequency of the controlled object including the phase delay and gain based on the signal of the second frequency component extracted from the signal of the manipulated amount and the signal of the second frequency component extracted from the signal of the controlled amount while calculating the characteristics. Is calculated, and the model to be controlled is calculated based on these calculated characteristics. Therefore, the value of the control parameter can be calculated based on the model of the controlled object calculated without performing the differential operation of the manipulated variable signal and the differential operation of the controlled variable signal. As a result, when calculating the value of the control parameter, the value of the control parameter can be calculated more accurately than the control system that performs the differential calculation of the manipulated variable signal or the differential operation of the controlled variable signal.
  • FIG. 1 shows the functional structure of the control system which concerns on embodiment of this disclosure.
  • Explanatory diagram of the method of calculating the phase lag from the extracted signal Explanatory diagram of the method of calculating the gain from the extracted signal
  • Flow chart showing the flow of control parameter change processing according to this embodiment
  • FIG. 2 for explaining the operation of the control system according to the present embodiment.
  • FIG. 3 for explaining the operation of the control system according to the present embodiment.
  • control system is a system in which a control device controls a control target while automatically adjusting the values of control parameters.
  • FIG. 1 is an overall explanatory view of the control system according to the embodiment of the present disclosure.
  • the control system 100 controls based on the control target 200, the control device 300 that controls the control target 200 based on the target amount SV, and the operation amount MV from the control device 300.
  • a detector 400 for detecting the controlled amount PV of the controlled target 200 is provided.
  • the control device 300 controls, for example, the temperature, humidity, pressure, flow rate, speed, position, etc. of the controlled object 200. Therefore, the detector 400 detects, for example, temperature, humidity, pressure, flow rate, speed, position, etc. as the controlled amount PV of the controlled object 200.
  • the control device 300 calculates the operation amount MV that eliminates the deviation E by using the subtractor 310 that calculates the deviation E between the target amount SV and the control amount PV, and the value of the held control parameter, and sets the control target 200.
  • a control target control unit 320 to be controlled is provided.
  • the control device 300 includes a frequency component extraction unit 330 that extracts the signal fs of a specific frequency component from the input signal, and an extraction result determination unit that determines whether or not the signal fs of the specific frequency component is normally extracted. It is equipped with 340.
  • control device 300 includes a frequency characteristic calculation unit 350 that calculates the characteristics of the specific frequency of the control target from the extracted signal fs of the specific frequency component, and the control device 300 is the control target from the calculated characteristics of the specific frequency of the control target.
  • a model calculation unit 360 for calculating a model is provided. Further, the control device 300 calculates and sets the control parameter value based on the change determination unit 370 that determines whether or not to change the control parameter value from the input signal, and the calculated model of the control target.
  • a control parameter changing unit 380 to be changed is provided.
  • the subtractor 310 calculates a deviation E obtained by subtracting the control amount PV from the target amount SV using the signal of the target amount SV input by the user and the signal of the control amount PV output from the detector 400, and calculates the deviation E.
  • the signal E is output to the control target control unit 320.
  • the control target control unit 320 is, for example, a controller such as a PI (Proportional Integral) controller or a PID controller.
  • the control target control unit 320 calculates an operation amount MV for eliminating the deviation E by using a preset control parameter value based on the signal of the deviation E input from the subtractor 310, and signals the operation amount MV. Is output to the control target 200 to control the control target 200. Further, the control target control unit 320 adjusts the control parameters when, for example, the target amount SV is changed, or the control amount PV changes to a value away from the target amount SV and the deviation E increases. , The signal of the manipulated variable MV is also output to the frequency component extraction unit 330.
  • the frequency component extraction unit 330 extracts the signal fsMV of a specific frequency component from the signal of the manipulated variable MV input from the control target control unit 320, and also extracts the signal fsMV of the specific frequency component from the signal of the controlled variable PV output from the detector 400.
  • the component signal fsPV is extracted.
  • the frequency component extraction unit 330 extracts signals of frequency components of the first frequency f1, the second frequency f2, and the third frequency f3. That is, the frequency component extraction unit 330 extracts the signal fsMV1 of the first frequency component, the signal fsMV2 of the second frequency component, and the signal fsMV3 of the third frequency component from the signal of the manipulated variable MV.
  • the frequency component extraction unit 330 extracts the signal fsPV1 of the first frequency component, the signal fsPV2 of the second frequency component, and the signal fsPV3 of the third frequency component from the signal of the controlled variable PV.
  • the first frequency f1, the second frequency f2, and the third frequency f3 are predetermined by the user of the control device 300.
  • the frequency component extraction unit 330 includes signals of the first frequency component fsMV1, fsPV1, signals of the second frequency component fsMV2, fsPV2, signals of the third frequency component fsMV3, from signals input using the bandpass filter and the resonance filter. Extract fsPV3.
  • the transfer functions BPF1 to BPF3 of the bandpass filter and the transfer functions RF1 to RF3 of the resonance filter are the mathematical formulas shown in (2-1) to (2-3) of the following equation 1 using the Laplace operator s. And the mathematical formulas shown in (3-1) to (3-3) of Equation 2.
  • the frequency component extraction unit 330 filters the signal of the manipulated variable MV with the bandpass filter of the transfer function BPF1 and then filters it with the resonance filter of the transfer function RF1 to obtain the first frequency component.
  • the signal fsMV1 is extracted and excited.
  • the frequency component extraction unit 330 filters the signal of the manipulated variable MV with the bandpass filter of the transfer function BPF2, and then filters the signal with the resonance filter of the transfer function RF2 to obtain the signal of the second frequency component. Extract and excite fsMV2.
  • the frequency component extraction unit 330 filters the signal of the manipulated variable MV with the bandpass filter of the transfer function BPF3, and then filters the signal with the resonance filter of the transfer function RF3 to obtain the signal of the third frequency component.
  • the frequency component extraction unit 330 filters the signal of the control amount PV with the bandpass filter of the transfer function BPF1 and then filters it with the resonance filter of the transfer function RF1 to obtain the signal of the first frequency component. Extract and excite fsPV1. Further, the frequency component extraction unit 330 filters the signal of the control amount PV with the bandpass filter of the transfer function BPF2, and then filters the signal with the resonance filter of the transfer function RF2 to obtain the signal of the second frequency component. Extract and excite fsPV2.
  • the frequency component extraction unit 330 filters the signal of the control amount PV with the bandpass filter of the transfer function BPF3, and then filters the signal with the resonance filter of the transfer function RF3 to obtain the signal of the third frequency component. Extract and excite fsPV3.
  • the transfer functions BPF1 to BPF3 of the bandpass filter are used.
  • the transfer functions RF1 to RF3 of the resonance filter are used as the transfer function of the continuous system, but the transfer function of the filter is not limited to this.
  • the transfer function of the filter may be a discrete transfer function.
  • the extraction result determination unit 340 determines whether or not the signal fs input from the frequency component extraction unit 330 is a signal that has been normally extracted, and outputs data Je indicating the extraction determination result to the frequency component extraction unit 330. ..
  • the extraction result determination unit 340 determines the frequency when the frequency of the signal extracted as the first frequency signals fsMV1 and fsPV1 is different from that of the first frequency f1, or when the amplitude of the extracted signal is smaller than the predetermined threshold value SA.
  • the component extraction unit 300 determines that the signal extraction has failed.
  • the extraction result determination unit 340 determines the frequency component when the frequency of the signal extracted as the second frequency component signals fsMV2 and fsPV2 is different from that of the second frequency f2, or when the amplitude of the extracted signal is smaller than the threshold SA. The extraction unit 300 determines that the signal extraction has failed. Further, the extraction result determination unit 340 determines the frequency component when the frequency of the signal extracted as the third frequency component signals fsMV3 and fsPV3 is different from that of the third frequency f3, or when the amplitude of the extracted signal is smaller than the threshold SA. The extraction unit 300 determines that the signal extraction has failed.
  • the frequency component extraction unit 330 was input from the control target control unit 320 when it was determined from the data Je indicating the extraction determination result input from the extraction result determination unit 340 that the extraction failed.
  • the signal fsMV of a specific frequency component is extracted again from the signal of the manipulated variable MV, and the signal fsPV of a specific frequency component is extracted again from the signal of the controlled variable PV output from the detector 400.
  • the frequency component extraction unit 330 identifies the signal that failed to be extracted from the data Je indicating the extraction determination result, and re-extracts the signal of the frequency component whose number of failed extractions exceeds a predetermined allowable number of times. Instead, extraction may be performed for a new frequency.
  • the frequency characteristic calculation unit 350 calculates the characteristics of the control target 200 based on the signal fs input from the frequency component extraction unit 330, and outputs the data if indicating the characteristics of the control target 200 to the model calculation unit 360.
  • the frequency characteristic calculation unit 350 calculates the phase delay ⁇ 1 and the gain R1 based on the signals fsMV1 and fsPV1 of the first frequency component extracted by the frequency component extraction unit 330. Further, the frequency characteristic calculation unit 350 calculates the phase delay ⁇ 2 and the gain R2 based on the signals fsMV2 and fsPV2 of the second frequency component extracted by the frequency component extraction unit 330. Further, the frequency characteristic calculation unit 350 calculates the phase delay ⁇ 3 and the gain R3 based on the signals fsMV3 and fsPV3 of the third frequency component extracted by the frequency component extraction unit 330.
  • FIG. 3 is an explanatory diagram of a method of calculating the phase lag from the extracted signal.
  • the vibration center of the signal fsMV extracted from the signal of the operation amount MV is c
  • the vibration center of the signal fsPV extracted from the signal of the control amount PV is c'
  • the signal of the operation amount MV is t p1 , t p2 , ..., T p4
  • the timing at which the maximum or minimum value of the signal fsPV extracted from the signal of the control amount PV is observed is t p1.
  • phase difference ⁇ i the average value of the phase difference
  • T c the frequencies of the extracted signals fsMV and fsPV.
  • the phase difference ⁇ i the average value ⁇ of the phase difference
  • the phase delay ⁇ are represented by the following mathematical formulas (4-1) to (4-3).
  • the frequency characteristic calculation unit 350 calculates the phase delay ⁇ 1 from the signals fsMV1 and fsPV1 of the first frequency component by using the mathematical formulas shown in (4-1) to (4-3) described above. Further, the frequency characteristic calculation unit 350 calculates the phase delay ⁇ 2 from the signals fsMV2 and fsPV2 of the second frequency component by using the mathematical formulas shown in (4-1) to (4-3) described above. Further, the frequency characteristic calculation unit 350 calculates the phase delay ⁇ 3 from the signals fsMV3 and fsPV3 of the third frequency component by using the mathematical formulas shown in (4-1) to (4-3) described above.
  • FIG. 4 is an explanatory diagram of a method of calculating the gain from the extracted signal.
  • the timing at which the signal fsMV extracted from the signal of the manipulated variable MV intersects the vibration center c is t 1 , t 2 , ..., T 5
  • the signal fsPV extracted from the signal of the controlled variable PV is 'the timing of crossing the t 1' vibration center c, t 2 ', ..., t 5', integrating 1 ⁇ i ⁇ 4 and to signal fsMV extracted from the manipulated variable MV when from t i until t (i + 1) the value of S MVi, the control amount of the extracted signal fsPV from PV 'from t (i + 1)' t i the integrated value to S PVi, the gain and R.
  • the integrated values SMVi , SMVi, and gain R are represented by the mathematical formulas shown in (5-1) to (5-3) of Equation 3 below.
  • the frequency characteristic calculation unit 350 calculates the gain R1 from the signals fsMV1 and fsPV1 of the first frequency component by using the mathematical formulas shown in (5-1) to (5-3) described above. Further, the frequency characteristic calculation unit 350 calculates the gain R2 from the signals fsMV2 and fsPV2 of the second frequency component by using the mathematical formulas shown in (5-1) to (5-3) described above. Further, the frequency characteristic calculation unit 350 calculates the gain R3 from the signals fsMV3 and fsPV3 of the third frequency component by using the mathematical formulas shown in (5-1) to (5-3) described above.
  • the data if indicating the characteristics of the control target 200 output by the frequency characteristic calculation unit 350 to the model calculation unit 360 includes the calculated phase delays ⁇ 1 to ⁇ 3 and the gains R1 to R3.
  • the frequency characteristic calculation unit 350 is as shown by the mathematical formulas shown in FIGS. 4, 5, (4-1) to (4-3), and (5-1) to (5-3).
  • the phase delay ⁇ and the gain R were calculated using the data of twice the vibration period of the extracted signals fsMV and fsPV, but the present invention is not limited to this.
  • the phase lag ⁇ and the gain R may be calculated using data for a time that is a positive integer multiple of the half cycle of the vibration cycle.
  • the model calculation unit 360 calculates the model of the control target 200 based on the phase delays ⁇ 1 to ⁇ 3 and the gains R1 to R3 included in the data if indicating the characteristics of the control target 200 input from the frequency characteristic calculation unit 350.
  • the data md indicating the calculated model is output to the control parameter changing unit 380.
  • the model of the controlled object 200 is, for example, a transfer function P (s) as a mathematical model of the controlled object 200.
  • the model calculation unit 360 calculates the transfer function P (s) of the control target 200 by calculating the coefficient of the transfer function P (s).
  • the structure of the transfer function P (s) is predetermined according to the structure of the control target control unit 320.
  • the phase delay characteristic of the PI controller is 0 [deg] or more and less than 90 [deg], so that the phase is delayed by 180 [deg] in the high frequency region.
  • a delay control target 200 is assumed.
  • the transfer function P (s) is represented by the mathematical formula shown in (6-1) of Equation 4 below, and the characteristic equation is a quadratic equation.
  • the phase lag characteristic calculated from the transfer function P (s) is ⁇ m and the gain characteristic is Rm
  • the phase lag ⁇ m and the gain Rm are the following equations 4 (6-2) and (6-3). ) Is shown by the formula.
  • the model calculation unit 360 has the phase delays ⁇ 1 to ⁇ 3 calculated by the frequency characteristic calculation unit 350, the gains R1 to R3, and the phase delays ⁇ m1 to ⁇ m3 and the gains Rm1 to Rm3 calculated from the transfer function P (s).
  • the phase delay characteristic of the PID controller may be 90 [deg] or more, so that the phase is delayed by 270 [deg] in the high frequency region 3
  • the next lag control target 200 is assumed.
  • the characteristic equation of the transfer function P (s) may be a cubic equation.
  • the change determination unit 370 determines whether or not to change the control parameter setting based on the control amount PV signal, and outputs the data Jc indicating the change determination result to the control parameter change unit 380.
  • the change determination unit 370 changes the setting of the control parameter depending on whether or not the rise time, the overshoot amount, the settling time, the attenuation ratio, etc. specified from the signal of the control amount PV satisfy the predetermined evaluation criteria. Determine whether or not to do so.
  • the control parameter change unit 380 indicates the model of the control target 200 input from the model calculation unit 360 when it is specified that the change is determined from the data Jc indicating the change determination result input from the change determination unit 370.
  • the value of a new control parameter is calculated based on the transfer function P (s) specified from the data md, and the setting of the control target control unit 320 is changed.
  • a known coefficient projection can be used as a method of calculating the value of the new control parameter from the transfer function P (s).
  • the coefficient projection is a method of adjusting the value of the control parameter of the controller based on the ratio of the coefficients of the specific polynomial which is the denominator of the transfer function P (s).
  • FIG. 2 is a block diagram showing a hardware configuration of the control device according to the present embodiment.
  • the control device 300 shown in FIG. 1 is composed of the hardware shown in FIG.
  • the control device 300 includes a control unit 51, a main storage unit 52, an external storage unit 53, an operation unit 54, a display unit 55, and a transmission / reception unit 56.
  • the control unit 51 executes the process according to the control program 59.
  • the control unit 51 includes a CPU (Central Processing Unit).
  • the control unit 51 controls the control target control unit 320, the frequency component extraction unit 330, the extraction result determination unit 340, the frequency characteristic calculation unit 350, the model calculation unit 360, the change determination unit 370, and the control device 300, according to the control program 59. It functions as a parameter change unit 380.
  • the control unit 51 performs a control target control step performed by the control target control unit 320, a frequency component extraction step performed by the frequency component extraction unit 330, an extraction result determination step performed by the extraction result determination unit 340, and a frequency characteristic calculation unit 350.
  • the frequency characteristic calculation step, the model calculation step performed by the model calculation unit 360, the change determination step performed by the change determination unit 370, and the control parameter change step performed by the control parameter change unit 380 are executed.
  • the main storage unit 52 loads the control program 59 and is used as a work area of the control unit 51.
  • the main storage unit 52 includes a RAM (Random-Access Memory).
  • the external storage unit 53 stores the control program 59 in advance.
  • the external storage unit 53 supplies the data stored in this program to the control unit 51 according to the instruction of the control unit 51, and stores the data supplied from the control unit 51.
  • the external storage unit 53 includes a non-volatile memory such as a flash memory, an HDD (Hard Disk Drive), and an SSD (Solid State Dive).
  • the operation unit 54 supplies the input information to the control unit 51.
  • the operation unit 54 includes information input components such as a keyboard, a mouse, a touch panel, and operation buttons.
  • the display unit 55 displays information input via the operation unit 54, information output by the control unit 51, and the like.
  • the display unit 55 includes display devices such as an LCD (Liquid Crystal Display) and an organic EL (Electro-Luminescence) display.
  • the transmission / reception unit 56 transmits / receives information.
  • the transmission / reception unit 56 includes information communication components such as a network termination device and a wireless communication device connected to the network.
  • the main storage unit 52, the external storage unit 53, the operation unit 54, the display unit 55, and the transmission / reception unit 56 are all connected to the control unit 51 via the internal bus 50.
  • control unit 51 uses the main storage unit 52, the external storage unit 53, the operation unit 54, the display unit 55, and the transmission / reception unit 56 as resources to extract the control target control unit 320 and the frequency component.
  • the functions of the unit 330, the extraction result determination unit 340, the frequency characteristic calculation unit 350, the model calculation unit 360, the change determination unit 370, and the control parameter change unit 380 are realized.
  • FIG. 5 is a flowchart showing the flow of the control parameter change process according to the present embodiment.
  • the control target control unit 320 determines whether or not the deviation E between the target amount SV calculated by the subtractor 310 and the control amount PV has increased (step S101). When the deviation E does not increase (step S101; N), the control target control unit 320 ends the control parameter change process.
  • the frequency component extraction unit 330 receives the signal fsMV1 of the first frequency component, the signal fsMV2 of the second frequency component, and the signal of the third frequency component from the signal of the manipulated variable MV.
  • the signal fsPV1 of the first frequency component, the signal fsPV2 of the second frequency component, and the signal fsPV3 of the third frequency component are extracted from the signal of the control amount PV (step S102).
  • the extraction result determination unit 340 determines whether or not the frequency of the signal extracted as the first frequency signal fsMV1 and fsPV1 is the same as that of the first frequency f1, and the extraction result determination unit 340 determines whether or not the frequency of the signal extracted as the second frequency component signal fsMV2 and fsPV2. It is determined whether or not the frequency is the same as the second frequency f2, and whether or not the frequency of the signals extracted as the signals fsMV3 and fsPV3 of the third frequency component is the same as the third frequency f3 (step S103).
  • the frequency component extraction unit 330 and the extraction result determination unit 340 repeat the processes of steps S102 and S103 until the signal extraction is successful. ..
  • the frequencies of the signals extracted as the first frequency signals fsMV1 and fsPV1 are the same as the first frequency f1
  • the frequencies of the signals extracted as the second frequency component signals fsMV2 and fsPV2 are the same as the second frequency f2.
  • the extraction result determination unit 340 determines the extracted signals fsMV1 to fsMV3. , It is determined whether or not the amplitude of fsPV1 to fsPV3 is larger than the threshold SA (step S104).
  • the frequency component extraction unit 330 and the extraction result determination unit 340 perform steps S102 to S104 until the signals are successfully extracted. Repeat the process of.
  • the frequency characteristic calculation unit 350 determines the phase delay ⁇ 1 based on the extracted signals fsMV1 to fsMV3 and fsPV1 to fsPV3.
  • step S105 the model calculation unit 360 calculates the transfer function P (s) based on the calculated phase delays ⁇ 1 to ⁇ 3 and the gains R1 to R3 (step S106).
  • the change determination unit 370 sets control parameters depending on whether or not the rise time, overshoot amount, settling time, attenuation ratio, etc. specified from the signal of the control amount PV satisfy predetermined evaluation criteria. It is determined whether or not it is necessary to change (step S107). When it is determined that it is not necessary to change the setting (step S107; N), the control parameter changing unit 380 ends the control parameter changing process without changing the value of the control parameter. On the other hand, when it is determined that the setting needs to be changed (step S107; Y), the control parameter changing unit 380 calculates a new control parameter value based on the calculated transfer function P (s), and controls the control target. The setting of the control unit 320 is changed (step S108), and the control parameter change process is terminated.
  • the control target control unit 320 sets the value of the preset control parameter based on the deviation E between the target quantity SV and the control quantity PV.
  • the operation amount MV for which the deviation E is eliminated is calculated by using the operation amount MV, and the signal of the operation amount MV is output to the control target 200 to control the control target 200.
  • the frequency component extraction unit 330 extracts the signal fsMV1 of the first frequency component, the signal fsMV2 of the second frequency component, and the signal fsMV3 of the third frequency component from the signal of the manipulated amount MV, and also extracts the signal fsMV3 of the third frequency component from the signal of the controlled amount PV.
  • the signal fsPV1 of the first frequency component, the signal fsPV2 of the second frequency component, and the signal fsPV3 of the third frequency component are extracted. Further, the frequency characteristic calculation unit 350 calculates the characteristics of the first frequency f1 of the controlled object 200 including the phase delay ⁇ 1 and the gain R1 based on the extracted signals fsMV1 and fsPV1 of the first frequency component. Further, the frequency characteristic calculation unit 350 calculates the characteristics of the second frequency f2 of the control target 200 including the phase delay ⁇ 2 and the gain R2 based on the extracted signals fsMV2 and fsPV2 of the second frequency component.
  • the frequency characteristic calculation unit 350 calculates the characteristics of the third frequency f3 of the control target 200 including the phase delay ⁇ 3 and the gain R3 based on the extracted signals fsMV3 and fsPV3 of the third frequency component.
  • the model calculation unit 360 calculates the model of the control target 200 based on the calculated characteristics of the first frequency f1 of the control target 200, the characteristics of the second frequency f2, and the characteristics of the third frequency f3.
  • the control parameter changing unit 380 calculates the value of the control parameter based on the calculated model of the controlled object 200 and changes the setting.
  • the value of the control parameter is set based on the model of the control target 200 calculated from the signal of the operation amount MV and the signal of the control amount PV that the control target control unit 320 actually controls the control target 200. Can be calculated.
  • the control parameter value is calculated more accurately than the control system that does not calculate the control parameter value based on the control target model calculated from the operation amount signal that actually controls the control target and the control amount signal. can.
  • the frequency f1 of the control target 200 including the phase delays ⁇ 1 to ⁇ 3 and the gains R1 to R3 based on the signals fsMV1 to fsMV3 and fsPV1 to fsPV3 extracted from the signals of the manipulated variable MV and the controlled variable PV.
  • the characteristics of ⁇ f3 are calculated, and the model of the controlled object 200 is calculated based on these calculated characteristics. Therefore, the value of the control parameter can be calculated based on the model of the control target 200 calculated without performing the differential calculation of the signal of the manipulated variable MV and the differential calculation of the signal of the controlled variable PV.
  • the value of the control parameter when calculating the value of the control parameter, the value of the control parameter can be calculated more accurately than the control system that performs the differential calculation of the manipulated variable signal or the differential operation of the controlled variable signal. Further, by doing so, the model of the control target 200 can be calculated without actually exciting the control target 200 during the calculation of the model of the control target 200. As a result, the control parameters can be adjusted without deteriorating the control performance as compared with the control system that actually excites the control target 200 during the calculation of the model of the control target 200.
  • FIG. 6 is an explanatory diagram 1 of the operation of the control system according to the present embodiment.
  • the detector 400 that actually outputs the control amount PV to the frequency component extraction unit 330 is omitted.
  • the signal path of the control device 300 includes path 1 which is a function of the frequency component extraction unit 330 to extract the signal fsMV of a specific frequency component from the signal of the operation amount MV, and the frequency component. It can be divided into a path 2 which is a function of the extraction unit 330 to extract the signal fsPV of a specific frequency component from the signal of the control amount PV.
  • FIG. 7 is an explanatory diagram 2 of the operation of the control system according to the present embodiment.
  • the frequency component extraction unit 330 controls the signal fsMV of a specific frequency component extracted from the signal of the manipulated variable MV. It is preferable to output to 200 to actually excite the controlled object 200 to acquire signals of the controlled variable PV', and then analyze these input / output signals.
  • the control target 200 is excited by outputting a signal that is generated only for calculating the model of the control target 200 and is not output to the control target 200 in actual control, such as the extracted signal fsMV, the control target 200 becomes the control target 200.
  • the plant control device described in Patent Document 2 described above generates a signal for identifying a controlled object, actually excites the controlled object, and then identifies the controlled object, and controls with the identification signal. Even if the controlled object can be identified accurately by exciting the object, the control performance may deteriorate.
  • FIG. 8 is an explanatory diagram 3 of the operation of the control system according to the present embodiment.
  • the above-mentioned path 2 can be equivalently exchanged from the relationship with the path 1.
  • the frequency component extraction unit 330 outputs the signal fsMV extracted from the signal of the manipulated variable MV to the controlled variable 200 and outputs the signal of the controlled variable fsPV when the controlled variable 200 is excited.
  • the signal of the controlled variable fsPV is equivalent to the signal of the controlled variable PV'obtained by actually exciting the controlled variable 200 with the signal fsMV extracted from the signal of the manipulated variable MV.
  • the control system 100 according to the present embodiment is controlled by using the signal fsPV equivalent to the signal of the control amount PV'obtained by actually exciting the control target 200 with the signal fsMV extracted from the signal of the operation amount MV.
  • the characteristics of the specific frequency of the target 200 can be calculated, and the model of the control target 200 can be calculated based on these calculated characteristics.
  • the control system 100 according to the present embodiment can accurately calculate the values of the control parameters.
  • the control system 100 according to the present embodiment can obtain the same signal fsMV without actually exciting the control target 200 with the signal fsMV extracted from the signal of the manipulated variable MV and acquiring the signal of the controlled variable PV'. You can get it.
  • the control system 100 according to the present embodiment can adjust the control parameters without deteriorating the control performance.
  • the frequency component extraction unit 330 extracts three types of frequency component signals fsMV1 to fsMV3 and fsPV1 to fsPV3 from the signal of the manipulated variable MV or the controlled variable PV.
  • the number is not limited to three.
  • Nth frequency component signals fsMVn and fsPVn may be added to make n types.
  • the frequency characteristic calculation unit 350 calculates the characteristics of the first frequency f1 and the characteristics of the second frequency f2 of the control target 200, and the model calculation unit 360 calculates the calculated control target 200.
  • the model of the controlled object 200 may be calculated based on the characteristics of the first frequency f1 and the characteristics of the second frequency f2.
  • the frequency characteristic calculation unit 350 calculates the characteristics of the first frequency f1 and the second frequency f2 of the controlled object 200, and ... the characteristics of the nth frequency fn, and models the model.
  • the calculation unit 360 may calculate the model of the control target 200 based on the calculated characteristics of the first frequency f1 of the control target 200, the characteristics of the second frequency f2, and the characteristics of the nth frequency fn.
  • the frequency component extraction unit 330 uses a resonance filter to obtain signals fsMV1, fsPV1 and a second frequency component of one frequency component from the signal of the manipulated variable MV or the controlled variable PV.
  • the signals fsMV2 and fsPV2 of the above and the signals fsMV3 and fsPV3 of the third frequency component are extracted, and the attenuation coefficient z of the resonance filter is 0.01, which is smaller than 0.2.
  • the influence of disturbance is more than that of the control system in which the signal of the 1st frequency component, the signal of the 2nd frequency component, and the signal of the 3rd frequency component are not extracted from the signal of the manipulated amount or the controlled amount by using the resonance filter. It is possible to accurately extract the signals fsMV1 and fsPV1 of one frequency component, the signals fsMV2 and fsPV2 of the second frequency component and the signals fsMV3 and fsPV3 of the third frequency component.
  • the influence of the disturbance is suppressed and the signals fsMV1 and fsPV1 of the first frequency component and the signals fsMV2 of the second frequency component are accurately suppressed as compared with the control system in which the attenuation coefficient of the resonance filter is 0.2 or more.
  • FsPV2 and the signals fsMV3 and fsPV3 of the third frequency component can be extracted.
  • the attenuation coefficient z of the resonance filter is set to 0.01, but the present invention is not limited to this, and is described above as long as the attenuation coefficient z of the resonance filter is an arbitrary value smaller than 0.2. It has an effect.
  • the frequency component extraction unit 330 extracts and excites the signal by using not only the resonance filter but also the bandpass filter, but the bandpass filter may be omitted.
  • the frequency component extraction unit 330 preferably extracts and excites the signal using a resonance filter having an attenuation coefficient z smaller than 0.2, but the attenuation is not limited to this.
  • the signal may be extracted and excited using a resonance filter having a coefficient z of 0.2 or more, or the signal may be extracted and excited without using a resonance filter. For example, the signal may be extracted and excited using only a bandpass filter.
  • the extraction result determination unit 340 extracts when the frequency of the signal extracted as the first frequency component signals fsMV1 and fsPV1 is different from that of the first frequency f1.
  • the frequency component extraction unit 300 determines that the signal extraction has failed.
  • the extraction result determination unit 340 determines the frequency component when the frequency of the signal extracted as the second frequency component signals fsMV2 and fsPV2 is different from that of the second frequency f2, or when the amplitude of the extracted signal is smaller than the threshold SA. The extraction unit 300 determines that the signal extraction has failed.
  • the extraction result determination unit 340 determines the frequency component when the frequency of the signal extracted as the third frequency component signals fsMV3 and fsPV3 is different from that of the third frequency f3, or when the amplitude of the extracted signal is smaller than the threshold SA.
  • the extraction unit 300 determines that the signal extraction has failed.
  • the frequency component extraction unit 330 extracts the signal again. By doing so, it is possible to calculate the characteristics of the frequencies f1 to f3 of the control target 200 using only the normally extracted signals fsMV1 to fsMV3 and fsPV1 to fsPV3 to calculate the model of the control target 200. Therefore, the model of the controlled object 200 can be calculated more systematically than the control system that does not calculate the characteristics of the frequency of the controlled object using only the normally extracted signals, and the values of the control parameters can be calculated accurately.
  • the frequency component extraction unit 330 retries all the processing of extracting the signal when it is determined that the extraction has failed, but the present invention is not limited to this. Only the extraction of the signal of the frequency component that failed to be extracted may be retried. As in the present embodiment, when it is determined that the extraction has failed, the frequency component extraction unit 330 preferably retries the process of extracting the signal, but extracts the signal. You do not have to retry the process. In this case, for example, the control parameter adjustment process may be interrupted or stopped because the signal necessary for calculating the model of the control target 200 could not be extracted.
  • the extraction result determination unit 340 determines that the frequency of the signal extracted as the signal fsMV or fsPV of the specific frequency component is different from the specific frequency, or the amplitude of the extracted signal is a threshold value.
  • SA it is preferable, but not limited to, it is determined that the signal extraction has failed.
  • the extraction result determination unit 340 does not have to determine the amplitude of the extracted signal extracted as the signals fsMV and fsPV of a specific frequency component. Further, for example, it is not necessary to omit the extraction result determination unit 340 and determine whether or not the signals extracted as the signals fsMV and fsPV of specific frequency components are normally extracted.
  • the model calculation unit 360 is calculated from the phase delays ⁇ 1 to ⁇ 3 calculated by the frequency characteristic calculation unit 350, the gains R1 to R3, and the transfer function P (s).
  • the coefficients a 0 , a 1 , and a 2 of the transfer function P (s) that minimize the errors of the phase delays ⁇ m1 to ⁇ m3 and the gains Rm1 to Rm3 are calculated.
  • the phase calculated from the phase delay and gain calculated by the frequency characteristic calculation unit 350 and the transfer function calculated from the transfer function, and the phase calculated by the control system that does not calculate the coefficient of the transfer function that minimizes the error of the gain are obtained.
  • the transfer function P (s) that systematically reproduces the delays ⁇ 1 to ⁇ 3 and the gains R1 to R3 can be calculated, and the value of the control parameter can be calculated accurately.
  • the model calculation unit 360 has the phase delays ⁇ 1 to ⁇ 3 calculated by the frequency characteristic calculation unit 350, the gains R1 to R3, and the phase delays ⁇ m1 to ⁇ m3 calculated from the transfer function P (s). It is preferable, but not limited to, the coefficients a 0 , a 1 , and a 2 of the transfer function P (s) that minimize the error of the gains Rm1 to Rm3. For example, in the model calculation unit 360, any one of the combinations ( ⁇ 1, R1) to ( ⁇ 3, R3) of the phase delay and the gain calculated by the frequency characteristic calculation unit 350 is calculated from the transfer function P (s).
  • the coefficients a 0 , a 1 , and a 2 of the transfer function P (s) that are the same as the combination of the phase delay and the gain ( ⁇ m1, Rm1) to ( ⁇ m3, Rm3) may be calculated.
  • the model calculation unit 360 calculates the model of the control target 200 by calculating the coefficients a 0 , a 1 , and a 2 of the transfer function P (s) of the control target 200.
  • the model of the target 200 is not limited to this as long as the value of the control parameter can be calculated.
  • the model of the control target 200 may be a relational expression between the transfer function P (s) of the control target 200 and the control parameters.
  • the first frequency f1, the second frequency f2, and the third frequency f3 are predetermined by the user of the control device 300.
  • the first frequency f1, the second frequency f2, and the second frequency f1 which can accurately calculate the model of the control target 200 and the value of the control parameter by reflecting the many years of experience and knowledge of a skilled engineer.
  • the user can set the three frequencies f3.
  • the first frequency f1, the second frequency f2, and the third frequency f3 are predetermined by the user of the control device 300, but are not limited thereto, and for example, the control target control unit.
  • the frequency band to be considered in the design of the control system 100 is determined from the values of the control parameters already held by the 320 and the values based on the known features of the controlled object 200, and the first frequency f1 within the determined frequency band is determined.
  • the second frequency f2 and the third frequency f3 may be determined.
  • the upper limit values of the center frequencies ⁇ B1 , ⁇ B2 , ⁇ B3 of the bandpass filter and the center frequencies ⁇ n1 , ⁇ n2 , and ⁇ n3 of the resonance filter that can be set by the sampling theorem based on the calculation cycle of the control device 300.
  • it is possible to set the first frequency f1, the second frequency f2, and the third frequency f3, which can accurately calculate the model of the control target 200 and the value of the control parameter by reflecting the information necessary for control. can.
  • the first frequency f1, the second frequency f2, and the third frequency f3 are all predetermined by the user, but the present invention is not limited to this.
  • the first frequency f1 is predetermined by the user
  • the second frequency f2 is based on the values of the control parameters already held by the control target control unit 320 as described above, and the values based on the known features of the control target 200. It is automatically determined by the control device 300 within the determined frequency band, and the third frequency f3 may be automatically determined by the control device 300 without any limitation.
  • the model of the control target 200 calculated by the model calculation unit 360 is the transfer function P (s) of the control target 200, and has a structure of the control target control unit 320.
  • the structure of the transfer function P (s) is predetermined. By doing so, it is possible to calculate a model of the control target 200 suitable for the structure of the control target control unit 320.
  • the structure of the transfer function P (s) is defined corresponding to the structure of the control target control unit 320 as in the present embodiment, but the structure of the transfer function P (s) is controlled. It does not have to correspond to the structure of the target control unit 320.
  • the characteristic equation of the transfer function P (s) does not have to be the minimum order, and may be a higher order equation having a higher order than the cubic equation.
  • the change determination unit 370 determines whether or not to change the setting of the control parameter based on the signal of the control amount PV. Then, when the change determination unit 370 determines that the control parameter setting needs to be changed, the control parameter change unit 380 determines a new control parameter value based on the model of the control target 200 calculated by the model calculation unit 360. And change the setting. By doing so, the value of the control parameter can be calculated and the setting can be changed only when it is necessary to change the setting of the control parameter. As in the present embodiment, when the change determination unit 370 determines that the control parameter setting needs to be changed, the control parameter change unit 380 calculates a new control parameter value and changes the setting. Is preferred, but not limited to. For example, when the deviation E increases and the model of the control target 200 is newly calculated, the control parameter changing unit 380 may always calculate the value of the new control parameter and change the setting. By doing so, the function of the change determination unit 370 can be omitted.
  • the coefficients a 0 , a 1 , and a 2 of the transfer function P (s) that minimize the errors of the delay ⁇ m1 to ⁇ m3 and the gains Rm1 to Rm3 have been calculated, but are not limited thereto.
  • the model calculation unit 360 a data i f indicating the characteristics of the controlled object 200, by the machine learning using learning data including the data md showing a model of the control object 200 calculated based on the characteristics of the controlled object 200
  • the data if indicating the characteristics of the control target 200 newly calculated by the frequency characteristic calculation unit 350 may be applied to the generated trained model to calculate the model of the control target 200.
  • FIG. 9 is an explanatory diagram of a method of calculating a model to be controlled using a neural network.
  • the model calculation unit 360 can use a known algorithm as the machine learning algorithm, for example, a neural network may be used.
  • the neural network according to the modified example of the present embodiment includes an input layer, an intermediate layer, and an output layer.
  • the data showing the characteristics of the first frequency f1 of the control target 200 is i f1
  • the second frequency of the control target 200 is the second frequency.
  • the model parameter A includes the coefficients a 0 , a 1 , and a 2 of the transfer function P (s), and A ⁇ ⁇ a 0 , a 1 , a 2 ⁇ holds.
  • data i f1 , if f2 , ..., Ifn indicating the characteristics of the controlled object 200 from the first frequency f1 to the nth frequency fn are input to the input layer, and the model parameter A is input from the output layer. Is output. Further, the weighting coefficient W1 in the connection between the node of the input layer and the node of the intermediate layer and the weighting coefficient W2 in the connection between the node of the intermediate layer and the node of the output layer can be set for each connection. There is.
  • a value obtained by multiplying the output value from the node of the input layer and the weighting coefficient W1 corresponding to the connection between the nodes is input to the node of the intermediate layer, and the node of the output layer is input from the node of the intermediate layer.
  • a value obtained by multiplying the output value by the weighting coefficient W2 corresponding to the connection between the nodes is input.
  • the output value from each node of the input layer, the intermediate layer, and the output layer may be the value of a linear function between the input values, and is the value of a non-linear function using the input value, for example, the value of the sigmoid function. You may.
  • the model calculation unit 360 uses data obtained by combining the data if indicating the characteristics of the controlled object 200 and the data md indicating the model of the controlled object 200 calculated based on the characteristics as learning data, and a plurality of types of learning data.
  • a neural network as a trained model can be generated. Any method can be adopted as the method for calculating the weighting coefficients W1 and W2, and for example, a known error propagation method or gradient descent method may be used.
  • the model calculation unit 360, the frequency characteristic calculating unit 350 outputs the model parameter A in accordance with the characteristics when the input to the neural network that generated the data i f indicating the characteristics of the controlled object 200 newly calculated Can be done.
  • the value of the control parameter can be set more than that of the control system that does not calculate the model to be controlled by applying the newly calculated data showing the characteristics of the controlled object to the trained model generated by machine learning. It can be calculated accurately.
  • the model calculation unit 360 generates the neural network by calculating the weighting coefficients W1 and W2 described above using the learning data, but the neural network is not limited to this, and for example, the control device 300 and the like. May generate a neural network by newly providing a different learning device.
  • the model calculation unit 360 generated a trained model using a neural network as an example of a machine learning algorithm, but the present invention is not limited to this.
  • a trained model may be generated by using a reinforcement learning algorithm such as Q-learning or TD-learning.
  • a trained model may be generated by using a known algorithm of supervised learning or unsupervised learning.
  • a trained model may be generated by using a known learning algorithm such as deep learning, genetic programming, functional logic programming, or a support vector machine.
  • the detector 400 is provided separately from the control device 300, but the present invention is not limited to this, and the detector 400 may be integrated with the control device 300, for example, the detector. 400 may be included in the control device 300.
  • the central part for processing the control device 300 including the control unit 51, the main storage unit 52, the external storage unit 53, the operation unit 54, the display unit 55, the transmission / reception unit 56, and the internal bus 50 is a dedicated system. However, it can be realized by using a normal computer system.
  • a computer program for executing the above operation is stored and distributed in a computer-readable recording medium such as a flexible disk, a DVD (Digital Versatile Disc) -ROM (Read-Only Memory), and the like.
  • a control device 300 that executes the above processing may be configured. Further, the control device 300 may be configured by storing the computer program in a storage device of the server device on the communication network and downloading it by a normal computer system.
  • control device 300 when the function of the control device 300 is realized by sharing the OS (operating system) and the application program, or when the function is realized by the cooperation between the OS and the application program, only the application program part is stored as a recording medium or storage. It may be stored in the device.
  • a computer program may be posted on a bulletin board system (BBS, Bulletin Board System) on a communication network, and the computer program may be provided via the network. Then, the above processing may be executed by starting this computer program and executing it in the same manner as other application programs under the control of the OS.
  • BSS bulletin board System
  • control unit 50 internal bus, 51 control unit, 52 main memory unit, 53 external storage unit, 54 operation unit, 55 display unit, 56 transmission / reception unit, 59 control program, 100 control system, 200 control target, 300 control device, 310 subtractor, 320 Control target control unit, 330 Frequency component extraction unit, 340 Extraction result determination unit, 350 Frequency characteristic calculation unit, 360 Model calculation unit, 370 Change determination unit, 380 Control parameter change unit, 400 Detector.

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

制御対象制御部（３２０）は、予め設定された制御パラメータの値を用いて偏差（Ｅ）が解消する操作量（ＭＶ）の信号を制御対象（２００）に出力して制御する。周波数成分抽出部（３３０）は、操作量（ＭＶ）の信号および制御量（ＰＶ）の信号から第１周波数成分の信号と第２周波数成分の信号とを抽出する。周波数特性算出部（３５０）は、抽出した第１周波数成分の信号と第２周波数成分の信号とに基づく位相遅れおよびゲインを含む制御対象（２００）の第１周波数の特性と第２周波数の特性とを算出する。モデル算出部（３６０）は、算出した制御対象（２００）の第１周波数の特性と第２周波数の特性とに基づいて制御対象（２００）のモデルを算出する。制御パラメータ変更部（３８０）は、算出した制御対象（２００）のモデルに基づいて制御パラメータの値を算出して設定を変更する。

Description

制御装置、制御システム、制御方法及びプログラム
　本開示は、制御装置、制御システム、制御方法及びプログラムに関する。
　制御対象を制御する制御システムにおいて、自動的に初期の制御パラメータの値を決定したり、実際の制御中のデータから制御パラメータの値を調整したりするチューニング機能を備えた制御装置が知られている。
　このようなチューニング機能を備えた制御装置の一例として、特許文献１には、制御対象に対して特別なチューニング信号を加えずに制御応答である測定値からＰＩＤ（Proportional Integral Differential）制御に用いられるＰＩＤ定数を求めて自動的に変更設定可能なＰＩＤ調節計が開示されている。特許文献１に記載のＰＩＤ調節計は、目標値と測定値の偏差がリミット値を超えるときに、目標値と測定値とに基づいて公知のステップ応答法を用いて第１のＰＩＤ定数を算出するとともに、公知のエキスパート法を用いて第２のＰＩＤ定数を算出する。また、特許文献１に記載のＰＩＤ調節計は、これらのＰＩＤ定数の確かさの度合いの評価に基づいて、既に設定されているＰＩＤ定数を何れかのＰＩＤ定数に変更したりこれらのＰＩＤ定数に基づく第３のＰＩＤ定数に変更したりする。
　また、このようなチューニング機能を備えた制御装置の一例として、特許文献２には、制御対象の動特性を表す数学モデルに従って制御器の制御パラメータの値を自動調整するプラント制御装置について開示されている。特許文献２に記載のプラント制御装置は、制御対象の既知情報を用いて生成された同定信号を制御器に入力したときの制御対象への操作量を表す入力信号と制御対象の制御量を表す出力信号とに基づいて公知の最小二乗法を用いて数学モデルを算出している。
特開平１０－６９３０１号公報 特開平１１－１０２２０３号公報
　特許文献１に記載のＰＩＤ調節計では、目標値と測定値の偏差がリミット値を超えたときに通常のステップ応答法における専用の入力信号であるステップ信号が入力されたとみなして第１のＰＩＤ定数を算出している。このため、偏差がリミット値を超えなければ第１のＰＩＤ定数を算出できず、偏差がリミット値を超えて算出できたとしても通常のステップ応答法で算出するよりも確かさの度合いが低い問題がある。また、特許文献１に記載のＰＩＤ調節計は、エキスパート法を用いて第２のＰＩＤ定数を算出しているが、エキスパート法は作業者の経験則に基づく推論に基づいてＰＩＤ定数を算出するため、第２のＰＩＤ定数の値が収束するまでの算出時間が長くなる問題があり、必ずしも確かさの度合いが高い値に収束しない問題もある。そして、特許文献１に記載のＰＩＤ調節計は、算出されたこれらのＰＩＤ定数の確かさの度合いの評価に基づいて変更するＰＩＤ定数を決定しているが、上述した問題から決定したＰＩＤ定数の精度自体が低い虞がある。
　また、特許文献２に記載のプラント制御装置では、制御対象の動特性を表す数学モデルが微分方程式で示される。このため、制御対象の入力信号と出力信号とに基づいて最小二乗法で数学モデルを算出しようとすると、これらの信号の微分演算を行う必要があるが、ノイズを含む信号の微分演算を行うとノイズの高周波成分を増幅させてしまう問題がある。よって、特許文献２に記載のプラント制御装置では、数学モデルを精度良く算出できない問題がある。
　本開示は、上述のような事情に鑑みてなされたものであり、制御パラメータの値を精度良く算出することを目的とする。
　上記目的を達成するため、本開示に係る制御装置は、制御対象制御部と、周波数成分抽出部と、周波数特性算出部と、モデル算出部と、制御パラメータ変更部とを備える。制御対象制御部は、目標量と制御対象の制御量との偏差に基づいて、予め設定された制御パラメータの値を用いて前記制御対象の操作量を算出し、算出した前記操作量の信号を制御対象に出力して制御対象を制御する。周波数成分抽出部は、操作量の信号から第１周波数成分の信号および第２周波数成分の信号を抽出するとともに、制御量の信号から第１周波数成分の信号および第２周波数成分の信号を抽出する。周波数特性算出部は、周波数成分抽出部が操作量の信号から抽出した第１周波数成分の信号と制御量の信号から抽出した第１周波数成分の信号とに基づく位相遅れとゲインとを含む制御対象の第１周波数の特性を算出するとともに、周波数成分抽出部が操作量の信号から抽出した第２周波数成分の信号と制御量の信号から抽出した第２周波数成分の信号とに基づく位相遅れおよびゲインを含む制御対象の第２周波数の特性を算出する。モデル算出部は、周波数特性算出部が算出した制御対象の第１周波数の特性と第２周波数の特性とに基づいて前記制御対象のモデルを算出する。制御パラメータ変更部は、モデル算出部が算出した制御対象のモデルに基づいて制御パラメータの値を算出して設定を変更する。
　本開示によれば、操作量の信号から抽出した第１周波数成分の信号と制御量の信号から抽出した第１周波数成分の信号とに基づく位相遅れとゲインとを含む制御対象の第１周波数の特性を算出するとともに、操作量の信号から抽出した第２周波数成分の信号と制御量の信号から抽出した第２周波数成分の信号とに基づく位相遅れおよびゲインを含む制御対象の第２周波数の特性を算出し、算出したこれらの特性に基づいて制御対象のモデルが算出される。このため、操作量の信号の微分演算および制御量の信号の微分演算を行わずに算出された制御対象のモデルに基づいて制御パラメータの値を算出できる。この結果、制御パラメータの値を算出するときに操作量の信号の微分演算または制御量の信号の微分演算を行う制御システムよりも制御パラメータの値を精度良く算出できる。
本開示の実施の形態に係る制御システムの機能構成を示す図 本実施の形態に係る制御装置のハードウェア構成を示すブロック図 抽出した信号から位相遅れを算出する方法の説明図 抽出した信号からゲインを算出する方法の説明図 本実施の形態に係る制御パラメータ変更処理の流れを示すフローチャート 本実施の形態に係る制御システムの作用説明図１ 本実施の形態に係る制御システムの作用説明図２ 本実施の形態に係る制御システムの作用説明図３ ニューラルネットワークを用いて制御対象のモデルを算出する方法の説明図
　以下に、本開示を実施するための形態に係る制御装置、制御システム、制御方法及びプログラムについて図面を参照して詳細に説明する。なお、図中同一または相当する部分には同じ符号を付す。
　本開示の実施の形態に係る制御システムは、制御装置が自動的に制御パラメータの値を調整しながら制御対象を制御するシステムである。
　図１は、本開示の実施の形態に係る制御システムの全体説明図である。
　本実施の形態に係る制御システム１００は、図１に示すように、制御対象２００、目標量ＳＶに基づいて制御対象２００を制御する制御装置３００、制御装置３００からの操作量ＭＶに基づいて制御された制御対象２００の制御量ＰＶを検出する検出器４００を備える。なお、制御装置３００は、例えば、制御対象２００の温度、湿度、圧力、流量、速度、位置等を制御する。このため、検出器４００は、制御対象２００の制御量ＰＶとして、例えば、温度、湿度、圧力、流量、速度、位置等を検出する。
　制御装置３００は、目標量ＳＶと制御量ＰＶとの偏差Ｅを算出する減算器３１０、保持している制御パラメータの値を用いて偏差Ｅを解消する操作量ＭＶを算出して制御対象２００を制御する制御対象制御部３２０を備える。また、制御装置３００は、入力された信号から特定の周波数成分の信号ｆｓを抽出する周波数成分抽出部３３０、特定の周波数成分の信号ｆｓが正常に抽出されたか否かを判定する抽出結果判定部３４０を備える。また、制御装置３００は、抽出された特定の周波数成分の信号ｆｓから制御対象の特定の周波数の特性を算出する周波数特性算出部３５０、算出された制御対象の特定の周波数の特性から制御対象のモデルを算出するモデル算出部３６０を備える。また、制御装置３００は、入力された信号から制御パラメータの値を変更するか否かを判定する変更判定部３７０、算出された制御対象のモデルに基づいて制御パラメータの値を算出して設定を変更する制御パラメータ変更部３８０を備える。
　減算器３１０は、使用者によって入力された目標量ＳＶの信号と検出器４００から出力された制御量ＰＶの信号とを用いて目標量ＳＶから制御量ＰＶを減算した偏差Ｅを算出し、偏差Ｅの信号を制御対象制御部３２０に出力する。
　制御対象制御部３２０は、例えば、ＰＩ（Proportional Integral）制御器、ＰＩＤ制御器等の制御器である。制御対象制御部３２０は、減算器３１０から入力された偏差Ｅの信号に基づいて、予め設定された制御パラメータの値を用いて偏差Ｅを解消する操作量ＭＶを算出し、操作量ＭＶの信号を制御対象２００に出力することで制御対象２００を制御する。また、制御対象制御部３２０は、例えば、目標量ＳＶが変更されたり、制御量ＰＶが目標量ＳＶから離れた値に変化して偏差Ｅが増大したりしたとき、制御パラメータの調整を行うため、操作量ＭＶの信号を周波数成分抽出部３３０にも出力する。
　周波数成分抽出部３３０は、制御対象制御部３２０から入力された操作量ＭＶの信号から特定の周波数成分の信号ｆｓＭＶを抽出するとともに、検出器４００から出力された制御量ＰＶの信号から特定の周波数成分の信号ｆｓＰＶを抽出する。周波数成分抽出部３３０は、第１周波数ｆ１、第２周波数ｆ２および第３周波数ｆ３の周波数成分の信号を抽出する。すなわち、周波数成分抽出部３３０は、操作量ＭＶの信号から第１周波数成分の信号ｆｓＭＶ１、第２周波数成分の信号ｆｓＭＶ２および第３周波数成分の信号ｆｓＭＶ３を抽出する。また、周波数成分抽出部３３０は、制御量ＰＶの信号から第１周波数成分の信号ｆｓＰＶ１、第２周波数成分の信号ｆｓＰＶ２および第３周波数成分の信号ｆｓＰＶ３を抽出する。なお、第１周波数ｆ１、第２周波数ｆ２、第３周波数ｆ３は、制御装置３００の使用者によって予め定められている。
　周波数成分抽出部３３０は、バンドパスフィルタと共振フィルタとを用いて入力された信号から第１周波数成分の信号ｆｓＭＶ１、ｆｓＰＶ１、第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２、第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３を抽出する。ここで、バンドパスフィルタの中心周波数ωＢ１、ωＢ２、ωＢ３と、共振フィルタの中心周波数ωｎ１、ωｎ２、ωｎ３とを以下の（１－１）～（１－３）に示す数式のように同じ値とし、共振フィルタの減衰係数ｚについて、ｚ＝０．０１が成立するものとする。
ωＢ１＝ωｎ１＝２×π×ｆ１　　　　…（１－１）
ωＢ２＝ωｎ２＝２×π×ｆ２　　　　…（１－２）
ωＢ３＝ωｎ３＝２×π×ｆ３　　　　…（１－３）
　このとき、バンドパスフィルタの伝達関数ＢＰＦ１～ＢＰＦ３と共振フィルタの伝達関数ＲＦ１～ＲＦ３とは、ラプラス演算子ｓを用いて以下の数１の（２－１）～（２－３）に示す数式と、数２の（３－１）～（３－３）に示す数式とで示される。
Figure JPOXMLDOC01-appb-M000001
Figure JPOXMLDOC01-appb-M000002
　この場合、周波数成分抽出部３３０は、操作量ＭＶの信号について、伝達関数ＢＰＦ１のバンドパスフィルタでフィルタ処理を行った後、伝達関数ＲＦ１の共振フィルタでフィルタ処理を行うことで第１周波数成分の信号ｆｓＭＶ１を抽出、励起する。また、周波数成分抽出部３３０は、操作量ＭＶの信号について、伝達関数ＢＰＦ２のバンドパスフィルタでフィルタ処理を行った後、伝達関数ＲＦ２の共振フィルタでフィルタ処理を行うことで第２周波数成分の信号ｆｓＭＶ２を抽出、励起する。また、周波数成分抽出部３３０は、操作量ＭＶの信号について、伝達関数ＢＰＦ３のバンドパスフィルタでフィルタ処理を行った後、伝達関数ＲＦ３の共振フィルタでフィルタ処理を行うことで第３周波数成分の信号ｆｓＭＶ３を抽出、励起する。また、周波数成分抽出部３３０は、制御量ＰＶの信号について、伝達関数ＢＰＦ１のバンドパスフィルタでフィルタ処理を行った後、伝達関数ＲＦ１の共振フィルタでフィルタ処理を行うことで第１周波数成分の信号ｆｓＰＶ１を抽出、励起する。また、周波数成分抽出部３３０は、制御量ＰＶの信号について、伝達関数ＢＰＦ２のバンドパスフィルタでフィルタ処理を行った後、伝達関数ＲＦ２の共振フィルタでフィルタ処理を行うことで第２周波数成分の信号ｆｓＰＶ２を抽出、励起する。そして、周波数成分抽出部３３０は、制御量ＰＶの信号について、伝達関数ＢＰＦ３のバンドパスフィルタでフィルタ処理を行った後、伝達関数ＲＦ３の共振フィルタでフィルタ処理を行うことで第３周波数成分の信号ｆｓＰＶ３を抽出、励起する。
　なお、本実施の形態では、上述した（２－１）～（２－３）、（３－１）～（３－３）に示す数式に示すように、バンドパスフィルタの伝達関数ＢＰＦ１～ＢＰＦ３と共振フィルタの伝達関数ＲＦ１～ＲＦ３とを連続系の伝達関数としたが、フィルタの伝達関数についてはこれに限定されない。例えば、フィルタの伝達関数を離散系の伝達関数としてもよい。
　抽出結果判定部３４０は、周波数成分抽出部３３０から入力された信号ｆｓが正常に抽出された信号であるか否かを判定し、抽出判定結果を示すデータＪｅを周波数成分抽出部３３０に出力する。抽出結果判定部３４０は、第１周波数の信号ｆｓＭＶ１、ｆｓＰＶ１として抽出した信号の周波数が第１周波数ｆ１と異なる場合、または、抽出した信号の振幅が予め定められた閾値ＳＡよりも小さい場合、周波数成分抽出部３００が信号の抽出に失敗したと判定する。また、抽出結果判定部３４０は、第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２として抽出した信号の周波数が第２周波数ｆ２と異なる場合、または、抽出した信号の振幅が閾値ＳＡよりも小さい場合、周波数成分抽出部３００が信号の抽出に失敗したと判定する。また、抽出結果判定部３４０は、第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３として抽出した信号の周波数が第３周波数ｆ３と異なる場合、または、抽出した信号の振幅が閾値ＳＡよりも小さい場合、周波数成分抽出部３００が信号の抽出に失敗したと判定する。
　なお、周波数成分抽出部３３０は、抽出結果判定部３４０から入力された抽出判定結果を示すデータＪｅから抽出に失敗したと判定されたことが特定されたとき、制御対象制御部３２０から入力された操作量ＭＶの信号から特定の周波数成分の信号ｆｓＭＶを再度抽出するとともに、検出器４００から出力された制御量ＰＶの信号から特定の周波数成分の信号ｆｓＰＶを再度抽出する。このとき、周波数成分抽出部３３０は、抽出判定結果を示すデータＪｅから抽出に失敗した信号を特定し、抽出に失敗した回数が予め定めた許容回数を超えた周波数成分の信号については再抽出を行わず、新たな周波数について抽出を行ってもよい。
　周波数特性算出部３５０は、周波数成分抽出部３３０から入力された信号ｆｓに基づいて制御対象２００の特性を算出し、制御対象２００の特性を示すデータiｆをモデル算出部３６０に出力する。周波数特性算出部３５０は、周波数成分抽出部３３０が抽出した第１周波数成分の信号ｆｓＭＶ１、ｆｓＰＶ１に基づく位相遅れθ１およびゲインＲ１を算出する。また、周波数特性算出部３５０は、周波数成分抽出部３３０が抽出した第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２に基づく位相遅れθ２およびゲインＲ２を算出する。また、周波数特性算出部３５０は、周波数成分抽出部３３０が抽出した第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３に基づく位相遅れθ３およびゲインＲ３を算出する。
　図３は、抽出した信号から位相遅れを算出する方法の説明図である。
　具体的には、図３に示すように、操作量ＭＶの信号から抽出した信号ｆｓＭＶの振動中心をｃ、制御量ＰＶの信号から抽出した信号ｆｓＰＶの振動中心をｃ′、操作量ＭＶの信号から抽出した信号ｆｓＭＶの極大値または極小値を観測したタイミングをｔｐ１，ｔｐ２，…，ｔｐ４、制御量ＰＶの信号から抽出した信号ｆｓＰＶの極大値または極小値を観測したタイミングをｔｐ１′，ｔｐ２′，…，ｔｐ４′、１≦ｉ≦４としたときに抽出した信号ｆｓＭＶ、ｆｓＰＶのタイミングｔｐｉ，ｔｐｉ′の位相差をφｉ、抽出した信号ｆｓＭＶ、ｆｓＰＶの位相差の平均値をφ、位相遅れをθとし、抽出した信号ｆｓＭＶ、ｆｓＰＶの周波数をＴｃする。この場合、位相差φｉ、位相差の平均値φ、位相遅れθは、以下の（４－１）～（４－３）に示す数式で示される。
φｉ＝ｔｐ１′－ｔｐ１　　　　　　　　　　 …（４－１）
φ＝（φ１＋φ２＋φ３＋φ４）／４ｐ　　　 …（４－２）
θ＝－２πφ／Ｔｃ　　　　 　　　　　   …（４－３）
　よって、周波数特性算出部３５０は、上述した（４－１）～（４－３）に示す数式を用いて第１周波数成分の信号ｆｓＭＶ１、ｆｓＰＶ１から位相遅れθ１を算出する。また、周波数特性算出部３５０は、上述した（４－１）～（４－３）に示す数式を用いて第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２から位相遅れθ２を算出する。また、周波数特性算出部３５０は、上述した（４－１）～（４－３）に示す数式を用いて第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３から位相遅れθ３を算出する。
　図４は、抽出した信号からゲインを算出する方法の説明図である。
　また、図４に示すように、操作量ＭＶの信号から抽出した信号ｆｓＭＶが振動中心ｃと交差するタイミングをｔ１，ｔ２，…，ｔ５、制御量ＰＶの信号から抽出した信号ｆｓＰＶが振動中心ｃ′と交差するタイミングをｔ１′，ｔ２′，…，ｔ５′、１≦ｉ≦４としたときに操作量ＭＶから抽出した信号ｆｓＭＶをｔｉからｔ（ｉ＋１）まで積分した値をＳＭＶｉ、制御量ＰＶから抽出した信号ｆｓＰＶをｔｉ′からｔ（ｉ＋１）′まで積分した値をＳＰＶｉ、ゲインをＲとする。この場合、積分した値ＳＭＶｉ、ＳＭＶｉおよびゲインＲは、以下の数３の（５－１）～（５－３）に示す数式で示される。
Figure JPOXMLDOC01-appb-M000003
　よって、周波数特性算出部３５０は、上述した（５－１）～（５－３）に示す数式を用いて第１周波数成分の信号ｆｓＭＶ１、ｆｓＰＶ１からゲインＲ１を算出する。また、周波数特性算出部３５０は、上述した（５－１）～（５－３）に示す数式を用いて第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２からゲインＲ２を算出する。また、周波数特性算出部３５０は、上述した（５－１）～（５－３）に示す数式を用いて第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３からゲインＲ３を算出する。なお、周波数特性算出部３５０がモデル算出部３６０に出力する制御対象２００の特性を示すデータiｆには算出した位相遅れθ１～θ３とゲインＲ１～Ｒ３とが含まれる。
　なお、本実施の形態では、周波数特性算出部３５０は、図４、図５、（４－１）～（４－３）、（５－１）～（５－３）に示す数式で示すように、抽出した信号ｆｓＭＶ、ｆｓＰＶの振動周期の２倍の時間のデータを用いて位相遅れθおよびゲインＲを算出したが、これに限定されない。例えば、振動周期の半周期の正の整数倍の時間のデータを用いて位相遅れθおよびゲインＲを算出してもよい。
　モデル算出部３６０は、周波数特性算出部３５０から入力された制御対象２００の特性を示すデータiｆに含まれる位相遅れθ１～θ３およびゲインＲ１～Ｒ３に基づいて制御対象２００のモデルを算出し、算出したモデルを示すデータｍｄを制御パラメータ変更部３８０に出力する。ここで、制御対象２００のモデルとは、例えば、制御対象２００の数式モデルとしての伝達関数Ｐ（ｓ）である。この場合、モデル算出部３６０は、伝達関数Ｐ（ｓ）の係数を算出することで、制御対象２００の伝達関数Ｐ（ｓ）を算出する。なお、伝達関数Ｐ（ｓ）の構造は、制御対象制御部３２０の構造に対応して予め定められている。例えば、制御対象制御部３２０としてＰＩ制御器を用いる場合、ＰＩ制御器の位相遅れの特性が０［ｄｅｇ］以上９０［ｄｅｇ］未満であるため、高周波領域で位相が１８０［ｄｅｇ］遅れる２次遅れの制御対象２００が想定される。この場合、伝達関数Ｐ（ｓ）は、以下の数４の（６－１）に示す数式で示され、特性方程式が二次方程式となる。また、伝達関数Ｐ（ｓ）から算出される位相遅れの特性をθｍ、ゲインの特性をＲｍとすると、位相遅れθｍ、ゲインＲｍは、以下の数４の（６－２）、（６－３）に示す数式で示される。
Figure JPOXMLDOC01-appb-M000004
　このため、伝達関数Ｐ（ｓ）から算出される第１周波数ｆ１の位相遅れθｍ１およびゲインＲｍ１は、上述した（６－２）、（６－３）に示す数式においてω＝２πｆ１とすることで算出できる。また、伝達関数Ｐ（ｓ）から算出される第２周波数ｆ２の位相遅れθｍ２およびゲインＲｍ２は、上述した（６－２）、（６－３）に示す数式においてω＝２πｆ２とすることで算出できる。そして、伝達関数Ｐ（ｓ）から算出される第３周波数ｆ３の位相遅れθｍ３およびゲインＲｍ３は、上述した（６－２）、（６－３）に示す数式においてω＝２πｆ３とすることで算出できる。
　よって、この場合、モデル算出部３６０は、周波数特性算出部３５０が算出した位相遅れθ１～θ３およびゲインＲ１～Ｒ３と伝達関数Ｐ（ｓ）から算出される位相遅れθｍ１～θｍ３およびゲインＲｍ１～Ｒｍ３の誤差が最小となる伝達関数Ｐ（ｓ）の係数ａ０，ａ１，ａ２を算出する。なお、誤差が最小となる係数ａ０，ａ１，ａ２を算出する方法としては、例えば、公知の最小二乗法を用いることが可能である。また、可能であればθ１＝θｍ１、θ２＝θｍ２、θ３＝θｍ３の連立方程式を解くことで係数ａ０，ａ１，ａ２を算出してもよい。
　なお、例えば、制御対象制御部３２０としてＰＩＤ制御器を用いる場合、ＰＩＤ制御器の位相遅れの特性が９０［ｄｅｇ］以上となる可能性があるため、高周波領域で位相が２７０［ｄｅｇ］遅れる３次遅れの制御対象２００が想定される。この場合、伝達関数Ｐ（ｓ）の特性方程式が３次方程式となってもよい。
　変更判定部３７０は、制御量ＰＶの信号に基づいて制御パラメータの設定を変更するか否かを判定し、変更判定結果を示すデータＪｃを制御パラメータ変更部３８０に出力する。変更判定部３７０は、例えば、制御量ＰＶの信号から特定された立ち上がり時間、オーバーシュート量、整定時間、減衰比等が予め定められた評価基準を満たすか否かによって、制御パラメータの設定を変更するか否かを判定する。
　制御パラメータ変更部３８０は、変更判定部３７０から入力された変更判定結果を示すデータＪｃから変更すると判定されたことが特定されたとき、モデル算出部３６０から入力された制御対象２００のモデルを示すデータｍｄから特定される伝達関数Ｐ（ｓ）に基づいて新たな制御パラメータの値を算出し、制御対象制御部３２０の設定を変更する。なお、伝達関数Ｐ（ｓ）から新たな制御パラメータの値を算出する方法としては、例えば、公知の係数図法を用いることが可能である。なお、係数図法とは、伝達関数Ｐ（ｓ）の分母である特定多項式の係数の比率に基づいて制御器の制御パラメータの値を調整する手法である。
　図２は、本実施の形態に係る制御装置のハードウェア構成を示すブロック図である。
　図１に示す制御装置３００は、図２に示すハードウェアによって構成されている。制御装置３００は、制御部５１、主記憶部５２、外部記憶部５３、操作部５４、表示部５５、送受信部５６を備える。
　制御部５１は、制御プログラム５９に従って処理を実行する。制御部５１は、ＣＰＵ（Central Processing Unit）を備える。制御部５１は、制御プログラム５９に従って、制御装置３００が備える制御対象制御部３２０、周波数成分抽出部３３０、抽出結果判定部３４０、周波数特性算出部３５０、モデル算出部３６０、変更判定部３７０、制御パラメータ変更部３８０として機能する。例えば、制御部５１は、制御対象制御部３２０が行う制御対象制御ステップ、周波数成分抽出部３３０が行う周波数成分抽出ステップ、抽出結果判定部３４０が行う抽出結果判定ステップ、周波数特性算出部３５０が行う周波数特性算出ステップ、モデル算出部３６０が行うモデル算出ステップ、変更判定部３７０が行う変更判定ステップ、制御パラメータ変更部３８０が行う制御パラメータ変更ステップを実行する。
　主記憶部５２は、制御プログラム５９をロードし、制御部５１の作業領域として用いられる。主記憶部５２は、ＲＡＭ（Random-Access Memory）を備える。
　外部記憶部５３は、制御プログラム５９を予め記憶する。外部記憶部５３は、制御部５１の指示に従って、このプログラムが記憶するデータを制御部５１に供給し、制御部５１から供給されたデータを記憶する。外部記憶部５３は、フラッシュメモリ、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Dive）等の不揮発性メモリを備える。
　操作部５４は、入力情報を制御部５１に供給する。操作部５４は、キーボード、マウス、タッチパネル、操作ボタン等の情報入力部品を備える。
　表示部５５は、操作部５４を介して入力された情報、制御部５１が出力した情報等を表示する。表示部５５は、ＬＣＤ（Liquid Crystal Display）、有機ＥＬ（Electro-Luminescence）ディスプレイ等の表示装置を備える。
　送受信部５６は、情報を送受信する。送受信部５６は、ネットワークに接続する網終端装置、無線通信装置等の情報通信部品を備える。
　制御装置３００では、主記憶部５２、外部記憶部５３、操作部５４、表示部５５及び送受信部５６はいずれも内部バス５０を介して制御部５１に接続されている。
　図１に示す制御装置３００は、制御部５１が主記憶部５２、外部記憶部５３、操作部５４、表示部５５及び送受信部５６を資源として用いることによって、制御対象制御部３２０、周波数成分抽出部３３０、抽出結果判定部３４０、周波数特性算出部３５０、モデル算出部３６０、変更判定部３７０、制御パラメータ変更部３８０の機能を実現する。
　図５は、本実施の形態に係る制御パラメータ変更処理の流れを示すフローチャートである。
　次に、図５に示すフローチャートを用いて制御装置３００が制御パラメータを自動調整する動作について説明する。先ず、図５に示すように、制御対象制御部３２０は、減算器３１０が算出した目標量ＳＶと制御量ＰＶとの偏差Ｅが増大したか否かを判定する（ステップＳ１０１）。偏差Ｅが増大しない場合（ステップＳ１０１；Ｎ）、制御対象制御部３２０は、制御パラメータ変更処理を終了する。一方、偏差Ｅが増大した場合（ステップＳ１０１；Ｙ）、周波数成分抽出部３３０は、操作量ＭＶの信号から第１周波数成分の信号ｆｓＭＶ１、第２周波数成分の信号ｆｓＭＶ２、第３周波数成分の信号ｆｓＭＶ３を抽出するとともに、制御量ＰＶの信号から第１周波数成分の信号ｆｓＰＶ１、第２周波数成分の信号ｆｓＰＶ２、第３周波数成分の信号ｆｓＰＶ３を抽出する（ステップＳ１０２）。
　信号抽出後、抽出結果判定部３４０は、第１周波数の信号ｆｓＭＶ１、ｆｓＰＶ１として抽出した信号の周波数が第１周波数ｆ１と同一か否か、第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２として抽出した信号の周波数が第２周波数ｆ２と同一か否か、第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３として抽出した信号の周波数が第３周波数ｆ３と同一か否かを判定する（ステップＳ１０３）。
　抽出した信号ｆｓＭＶ１～ｆｓＭＶ３、ｆｓＰＶ１～ｆｓＰＶ３の周波数が異なる場合（ステップＳ１０３；Ｎ）、周波数成分抽出部３３０および抽出結果判定部３４０は、信号の抽出に成功するまでステップＳ１０２、Ｓ１０３の処理を繰り返す。一方、第１周波数の信号ｆｓＭＶ１、ｆｓＰＶ１として抽出した信号の周波数が第１周波数ｆ１と同一であり、且つ、第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２として抽出した信号の周波数が第２周波数ｆ２と同一であり、且つ、第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３として抽出した信号の周波数が第３周波数ｆ３と同一である場合（ステップＳ１０３；Ｙ）、抽出結果判定部３４０は、抽出した信号ｆｓＭＶ１～ｆｓＭＶ３、ｆｓＰＶ１～ｆｓＰＶ３の振幅が閾値ＳＡよりも大きいか否かを判定する（ステップＳ１０４）。
　抽出した信号ｆｓＭＶ１～ｆｓＭＶ３、ｆｓＰＶ１～ｆｓＰＶ３の振幅が閾値ＳＡ以下となる場合（ステップＳ１０４；Ｎ）、周波数成分抽出部３３０および抽出結果判定部３４０は、信号の抽出に成功するまでステップＳ１０２～Ｓ１０４の処理を繰り返す。一方、抽出した信号ｆｓＭＶ１～ｆｓＭＶ３、ｆｓＰＶ１～ｆｓＰＶ３の振幅が閾値ＳＡよりも大きい場合（ステップＳ１０４；Ｙ）、周波数特性算出部３５０は、抽出した信号ｆｓＭＶ１～ｆｓＭＶ３、ｆｓＰＶ１～ｆｓＰＶ３に基づく位相遅れθ１～θ３およびゲインＲ１～Ｒ３を算出する（ステップＳ１０５）。特性算出後、モデル算出部３６０は、算出した位相遅れθ１～θ３およびゲインＲ１～Ｒ３に基づいて伝達関数Ｐ（ｓ）を算出する（ステップＳ１０６）。
　モデル算出後、変更判定部３７０は、制御量ＰＶの信号から特定された立ち上がり時間、オーバーシュート量、整定時間、減衰比等が予め定められた評価基準を満たすか否かによって、制御パラメータの設定を変更する必要があるか否かを判定する（ステップＳ１０７）。設定を変更する必要がないと判定した場合（ステップＳ１０７；Ｎ）、制御パラメータ変更部３８０は、制御パラメータの値を変更せずに制御パラメータ変更処理を終了する。一方、設定を変更する必要があると判定した場合（ステップＳ１０７；Ｙ）、制御パラメータ変更部３８０は、算出した伝達関数Ｐ（ｓ）に基づいて新たな制御パラメータの値を算出し、制御対象制御部３２０の設定を変更し（ステップＳ１０８）、制御パラメータ変更処理を終了する。
　以上説明したように、本実施の形態に係る制御システム１００によれば、制御対象制御部３２０は、目標量ＳＶと制御量ＰＶとの偏差Ｅに基づいて、予め設定された制御パラメータの値を用いて偏差Ｅが解消する操作量ＭＶを算出し、操作量ＭＶの信号を制御対象２００に出力して制御対象２００を制御する。また、周波数成分抽出部３３０は、操作量ＭＶの信号から第１周波数成分の信号ｆｓＭＶ１、第２周波数成分の信号ｆｓＭＶ２、第３周波数成分の信号ｆｓＭＶ３を抽出するとともに、制御量ＰＶの信号から第１周波数成分の信号ｆｓＰＶ１、第２周波数成分の信号ｆｓＰＶ２、第３周波数成分の信号ｆｓＰＶ３を抽出する。また、周波数特性算出部３５０は、抽出した第１周波数成分の信号ｆｓＭＶ１、ｆｓＰＶ１に基づく位相遅れθ１およびゲインＲ１を含む制御対象２００の第１周波数ｆ１の特性を算出する。また、周波数特性算出部３５０は、抽出した第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２に基づく位相遅れθ２およびゲインＲ２を含む制御対象２００の第２周波数ｆ２の特性を算出する。また、周波数特性算出部３５０は、抽出した第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３に基づく位相遅れθ３およびゲインＲ３を含む制御対象２００の第３周波数ｆ３の特性を算出する。また、モデル算出部３６０は、算出した制御対象２００の第１周波数ｆ１の特性と第２周波数ｆ２の特性と第３周波数ｆ３の特性とに基づいて制御対象２００のモデルを算出する。そして、制御パラメータ変更部３８０は、算出した制御対象２００のモデルに基づいて制御パラメータの値を算出して設定を変更する。
　このようにすることで、制御対象制御部３２０が実際に制御対象２００を制御している操作量ＭＶの信号と制御量ＰＶの信号から算出した制御対象２００のモデルに基づいて制御パラメータの値を算出できる。この結果、実際に制御対象を制御している操作量の信号と制御量の信号から算出した制御対象のモデルに基づいて制御パラメータの値を算出しない制御システムよりも制御パラメータの値を精度良く算出できる。また、このようにすることで、操作量ＭＶおよび制御量ＰＶの信号から抽出した信号ｆｓＭＶ１～ｆｓＭＶ３、ｆｓＰＶ１～ｆｓＰＶ３に基づく位相遅れθ１～θ３とゲインＲ１～Ｒ３とを含む制御対象２００の周波数ｆ１～ｆ３の特性を算出し、算出したこれらの特性に基づいて制御対象２００のモデルが算出される。このため、操作量ＭＶの信号の微分演算および制御量ＰＶの信号の微分演算を行わずに算出された制御対象２００のモデルに基づいて制御パラメータの値を算出できる。この結果、制御パラメータの値を算出するときに操作量の信号の微分演算または制御量の信号の微分演算を行う制御システムよりも制御パラメータの値を精度良く算出できる。さらに、このようにすることで、制御対象２００のモデルの算出中に実際に制御対象２００を励起しなくても制御対象２００のモデルを算出できる。この結果、制御対象２００のモデルの算出中に実際に制御対象２００を励起する制御システムよりも制御性能を悪化させることなく制御パラメータを調整できる。
　図６は、本実施の形態に係る制御システムの作用説明図１である。なお、説明を簡略化するため、図６では周波数成分抽出部３３０に制御量ＰＶを実際に出力している検出器４００が省略されている。
　ここで、図６に示すように、制御装置３００の信号の経路は、周波数成分抽出部３３０が操作量ＭＶの信号から特定の周波数成分の信号ｆｓＭＶを抽出する機能であるパス１と、周波数成分抽出部３３０が制御量ＰＶの信号から特定の周波数成分の信号ｆｓＰＶを抽出する機能であるパス２とに分けることができる。
　図７は、本実施の形態に係る制御システムの作用説明図２である。
　モデル算出部３６０が制御対象２００のモデルを精度良く算出するためには、図７に示すように、周波数成分抽出部３３０が操作量ＭＶの信号から抽出した特定の周波数成分の信号ｆｓＭＶを制御対象２００に出力して実際に制御対象２００を励起させて制御量ＰＶ′の信号を取得した上でこれらの入出力信号の解析を行うことが好ましい。しかしながら、抽出した信号ｆｓＭＶのように、制御対象２００のモデルを算出するためだけに生成されて実際の制御では制御対象２００に出力されない信号を出力して制御対象２００を励起すると、制御対象２００に何らかの悪影響を及ぼして制御性能が低下する虞がある。例えば、上述した特許文献２に記載のプラント制御装置は、制御対象の同定用の信号を生成して実際に制御対象を励起した上で制御対象の同定を行っており、同定用の信号で制御対象を励起したことで制御対象を精度よく同定できたとしても制御性能が低下する虞がある。
　図８は、本実施の形態に係る制御システムの作用説明図３である。
　これに対して、本実施の形態では、図８に示すように、上述したパス２はパス１との関係から等価交換が可能である。変換後のパス２では、周波数成分抽出部３３０が操作量ＭＶの信号から抽出した信号ｆｓＭＶを制御対象２００に出力して制御対象２００を励起させたときの制御量ｆｓＰＶの信号を出力する機能となる。このとき、制御量ｆｓＰＶの信号は、操作量ＭＶの信号から抽出した信号ｆｓＭＶで実際に制御対象２００を励起させて取得した制御量ＰＶ′の信号と等価である。
　したがって、本実施の形態に係る制御システム１００は、操作量ＭＶの信号から抽出した信号ｆｓＭＶで実際に制御対象２００を励起させて取得した制御量ＰＶ′の信号と同等の信号ｆｓＰＶを用いて制御対象２００の特定の周波数の特性を算出し、算出したこれらの特性に基づいて制御対象２００のモデルを算出できる。この結果、本実施の形態に係る制御システム１００は、制御パラメータの値を精度良く算出できる。また、本実施の形態に係る制御システム１００は、操作量ＭＶの信号から抽出した信号ｆｓＭＶで実際に制御対象２００を励起させて制御量ＰＶ′の信号を取得しなくても同等の信号ｆｓＭＶを取得できる。この結果、本実施の形態に係る制御システム１００は、制御性能を悪化させることなく制御パラメータを調整できる。
　なお、本実施の形態では、周波数成分抽出部３３０が操作量ＭＶまたは制御量ＰＶの信号から３種類の周波数成分の信号ｆｓＭＶ１～ｆｓＭＶ３、ｆｓＰＶ１～ｆｓＰＶ３を抽出したが、抽出する周波数成分の信号は２種類以上である限りにおいて３種類に限定されない。例えば、第１周波数成分の信号ｆｓＭＶ１、ｆｓＰＶ１と第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２のみの２種類としてもよく、第４周波数成分の信号ｆｓＭＶ４、ｆｓＰＶ４、第５周波数成分の信号ｆｓＭＶ５、ｆｓＰＶ５、…、第ｎ周波数成分の信号ｆｓＭＶｎ、ｆｓＰＶｎを追加してｎ種類としてもよい。抽出する信号を２種類とした場合、周波数特性算出部３５０は、制御対象２００の第１周波数ｆ１の特性と第２周波数ｆ２の特性とを算出し、モデル算出部３６０は、算出した制御対象２００の第１周波数ｆ１の特性と第２周波数ｆ２の特性とに基づいて制御対象２００のモデルを算出すればよい。また、抽出する信号をｎ種類とした場合、周波数特性算出部３５０は、制御対象２００の第１周波数ｆ１の特性と第２周波数ｆ２の特性と…第ｎ周波数ｆｎの特性とを算出し、モデル算出部３６０は、算出した制御対象２００の第１周波数ｆ１の特性と第２周波数ｆ２の特性と…第ｎ周波数ｆｎの特性とに基づいて制御対象２００のモデルを算出すればよい。
　また、本実施の形態に係る制御システム１００によれば、周波数成分抽出部３３０は、共振フィルタを用いて操作量ＭＶまたは制御量ＰＶの信号から１周波数成分の信号ｆｓＭＶ１、ｆｓＰＶ１と第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２と第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３とを抽出しており、共振フィルタの減衰係数ｚは、０．２よりも小さい０．０１である。
　このようにすることで、共振フィルタを用いて操作量または制御量の信号から１周波数成分の信号と第２周波数成分の信号と第３周波数成分の信号とを抽出しない制御システムよりも外乱の影響を抑制して精度良く１周波数成分の信号ｆｓＭＶ１、ｆｓＰＶ１と第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２と第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３とを抽出できる。また、このようにすることで、共振フィルタの減衰係数が０．２以上となる制御システムよりも外乱の影響を抑制して精度良く１周波数成分の信号ｆｓＭＶ１、ｆｓＰＶ１と第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２と第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３とを抽出できる。
　なお、本実施の形態では、共振フィルタの減衰係数ｚを０．０１としたが、これに限定されず、共振フィルタの減衰係数ｚが０．２よりも小さい任意の値である限りにおいて上述した作用効果を奏する。
　なお、本実施の形態では、周波数成分抽出部３３０は、共振フィルタだけでなくバンドパスフィルタも用いて信号を抽出、励起しているが、バンドパスフィルタについては省略してもよい。
　なお、本実施の形態のように、周波数成分抽出部３３０は、減衰係数ｚが０．２よりも小さい共振フィルタを用いて信号を抽出、励起することが好ましいが、これに限定されず、減衰係数ｚが０．２以上となる共振フィルタを用いて信号を抽出、励起してもよく、共振フィルタを用いずに信号を抽出、励起してもよい。例えば、バンドパスフィルタのみを用いて信号を抽出、励起してもよい。
　また、本実施の形態に係る制御システム１００によれば、抽出結果判定部３４０は、第１周波数成分の信号ｆｓＭＶ１、ｆｓＰＶ１として抽出した信号の周波数が第１周波数ｆ１と異なる場合、または、抽出した信号の振幅が閾値ＳＡよりも小さい場合、周波数成分抽出部３００が信号の抽出に失敗したと判定する。また、抽出結果判定部３４０は、第２周波数成分の信号ｆｓＭＶ２、ｆｓＰＶ２として抽出した信号の周波数が第２周波数ｆ２と異なる場合、または、抽出した信号の振幅が閾値ＳＡよりも小さい場合、周波数成分抽出部３００が信号の抽出に失敗したと判定する。また、抽出結果判定部３４０は、第３周波数成分の信号ｆｓＭＶ３、ｆｓＰＶ３として抽出した信号の周波数が第３周波数ｆ３と異なる場合、または、抽出した信号の振幅が閾値ＳＡよりも小さい場合、周波数成分抽出部３００が信号の抽出に失敗したと判定する。そして、周波数成分抽出部３３０は、抽出に失敗したと判定されたことが特定されたとき、信号の抽出を再度行う。
　このようにすることで、正常に抽出された信号ｆｓＭＶ１～ｆｓＭＶ３、ｆｓＰＶ１～ｆｓＰＶ３のみを用いて制御対象２００の周波数ｆ１～ｆ３の特性を算出して制御対象２００のモデルを算出できる。よって、正常に抽出された信号のみを用いて制御対象の周波数の特性を算出しない制御システムよりも制御対象２００のモデルを制度良く算出でき、制御パラメータの値を精度良く算出できる。
　なお、本実施の形態では、周波数成分抽出部３３０は、抽出に失敗したと判定されたことが特定されたとき、信号を抽出する処理を全て再試行しているが、これに限定されず、抽出に失敗した周波数成分の信号の抽出のみを再試行してもよい。
　なお、本実施の形態のように、周波数成分抽出部３３０は、抽出に失敗したと判定されたことが特定されたとき、信号を抽出する処理を再試行することが好ましいが、信号を抽出する処理を再試行しなくてもよい。この場合、例えば、制御対象２００のモデルを算出するために必要な信号を抽出できなかったとして制御パラメータ調整処理を中断したり、中止したりしてもよい。
　なお、本実施の形態のように、抽出結果判定部３４０は、特定の周波数成分の信号ｆｓＭＶ、ｆｓＰＶとして抽出した信号の周波数が当該特定の周波数と異なる場合、または、抽出した信号の振幅が閾値ＳＡよりも小さい場合、信号の抽出に失敗したと判定することが好ましいが、これに限定されない。例えば、抽出結果判定部３４０は、特定の周波数成分の信号ｆｓＭＶ、ｆｓＰＶとして抽出した抽出した信号の振幅について判定しなくてもよい。また、例えば、抽出結果判定部３４０を省略して特定の周波数成分の信号ｆｓＭＶ、ｆｓＰＶとして抽出した信号について正常に抽出されたか否かを判定しなくてもよい。
　また、本実施の形態に係る制御システム１００によれば、モデル算出部３６０は、周波数特性算出部３５０が算出した位相遅れθ１～θ３およびゲインＲ１～Ｒ３と伝達関数Ｐ（ｓ）から算出される位相遅れθｍ１～θｍ３およびゲインＲｍ１～Ｒｍ３の誤差が最小となる伝達関数Ｐ（ｓ）の係数ａ０，ａ１，ａ２を算出する。
　このようにすることで、周波数特性算出部３５０が算出した位相遅れおよびゲインと伝達関数から算出される位相遅れおよびゲインの誤差が最小となる伝達関数の係数を算出しない制御システムよりも算出した位相遅れθ１～θ３およびゲインＲ１～Ｒ３を制度良く再現する伝達関数Ｐ（ｓ）を算出でき、制御パラメータの値を精度良く算出できる。
　なお、本実施の形態のように、モデル算出部３６０は、周波数特性算出部３５０が算出した位相遅れθ１～θ３およびゲインＲ１～Ｒ３と伝達関数Ｐ（ｓ）から算出される位相遅れθｍ１～θｍ３およびゲインＲｍ１～Ｒｍ３の誤差が最小となる伝達関数Ｐ（ｓ）の係数ａ０，ａ１，ａ２を算出することが好ましいが、これに限定されない。例えば、モデル算出部３６０は、周波数特性算出部３５０が算出した位相遅れとゲインとの組み合わせ（θ１，Ｒ１）～（θ３，Ｒ３）のうちの何れかが伝達関数Ｐ（ｓ）から算出される位相遅れとゲインとの組み合わせ（θｍ１，Ｒｍ１）～（θｍ３，Ｒｍ３）と同一となる伝達関数Ｐ（ｓ）の係数ａ０，ａ１，ａ２を算出してもよい。
　なお、本実施の形態では、モデル算出部３６０は、制御対象２００の伝達関数Ｐ（ｓ）の係数ａ０，ａ１，ａ２を算出することで制御対象２００のモデルを算出したが、制御対象２００のモデルについては制御パラメータの値を算出できる限りにおいてこれに限定されない。例えば、制御対象２００のモデルは、制御対象２００の伝達関数Ｐ（ｓ）と制御パラメータとの関係式であってもよい。
　また、本実施の形態に係る制御システム１００によれば、第１周波数ｆ１、第２周波数ｆ２、第３周波数ｆ３は、制御装置３００の使用者によって予め定められている。
　このようにすることで、例えば、熟練の技術者の長年の経験、知見等を反映して精度良く制御対象２００のモデル、制御パラメータの値を算出できる第１周波数ｆ１、第２周波数ｆ２、第３周波数ｆ３を使用者が設定することができる。
　なお、本実施の形態では、第１周波数ｆ１、第２周波数ｆ２、第３周波数ｆ３は、制御装置３００の使用者によって予め定められているが、これに限定されず、例えば、制御対象制御部３２０が既に保持している制御パラメータの値、制御対象２００の既知の特徴に基づく値から制御システム１００の設計において考慮しなければならない周波数帯域を決定し、決定した周波数帯域内で第１周波数ｆ１、第２周波数ｆ２、第３周波数ｆ３を決定してもよい。また、例えば、制御装置３００の演算周期に基づいてサンプリング定理によって設定可能なバンドパスフィルタの中心周波数ωＢ１、ωＢ２、ωＢ３、共振フィルタの中心周波数ωｎ１、ωｎ２、ωｎ３の上限値が算出できるので、これらの中心周波数が上限値以下となる第１周波数ｆ１、第２周波数ｆ２、第３周波数ｆ３を決定してもよい。
　このようにすることで、制御に必要な情報を反映して精度良く制御対象２００のモデル、制御パラメータの値を算出できる第１周波数ｆ１、第２周波数ｆ２、第３周波数ｆ３を設定することができる。
　なお、本実施の形態では、第１周波数ｆ１、第２周波数ｆ２、第３周波数ｆ３が全て使用者によって予め定められているが、これに限定されない。例えば、第１周波数ｆ１が使用者によって予め定められ、第２周波数ｆ２が上述したように制御対象制御部３２０が既に保持している制御パラメータの値、制御対象２００の既知の特徴に基づく値から決定された周波数帯域内で制御装置３００によって自動的に決定されており、第３周波数ｆ３が何の制限もなく制御装置３００によって自動的に決定されていてもよい。
　また、本実施の形態に係る制御システム１００によれば、モデル算出部３６０が算出する制御対象２００のモデルは、制御対象２００の伝達関数Ｐ（ｓ）であり、制御対象制御部３２０の構造に対応して伝達関数Ｐ（ｓ）の構造が予め定められている。
　このようにすることで、制御対象制御部３２０の構造に適した制御対象２００のモデルを算出することができる。
　なお、本実施の形態のように、制御対象制御部３２０の構造に対応して伝達関数Ｐ（ｓ）の構造が定められていることが好ましいが、伝達関数Ｐ（ｓ）の構造は、制御対象制御部３２０の構造に対応していなくてもよい。例えば、伝達関数Ｐ（ｓ）の特性方程式が最低限の次数でなくてもよく、３次方程式よりも次数の高い高次方程式であってもよい。
　また、本実施の形態に係る制御システム１００によれば、変更判定部３７０は、制御量ＰＶの信号に基づいて制御パラメータの設定を変更するか否かを判定する。そして、制御パラメータ変更部３８０は、変更判定部３７０が制御パラメータの設定を変更する必要があると判定したとき、モデル算出部３６０が算出した制御対象２００のモデルに基づいて新たな制御パラメータの値を算出して設定を変更する。
　このようにすることで、制御パラメータの設定を変更する必要があるときのみ制御パラメータの値を算出して設定を変更することができる。
　なお、本実施の形態のように、変更判定部３７０が制御パラメータの設定を変更する必要があると判定したときに制御パラメータ変更部３８０が新たな制御パラメータの値を算出して設定を変更することが好ましいが、これに限定されない。例えば、偏差Ｅが増大して制御対象２００のモデルが新たに算出されたときには制御パラメータ変更部３８０が常に新たな制御パラメータの値を算出して設定を変更してもよい。このようにすることで、変更判定部３７０の機能を省略できる。
　なお、本実施の形態では、モデル算出部３６０が算出した制御対象２００の特性を示すデータiｆに基づいて位相遅れθ１～θ３およびゲインＲ１～Ｒ３と伝達関数Ｐ（ｓ）から算出される位相遅れθｍ１～θｍ３およびゲインＲｍ１～Ｒｍ３の誤差が最小となる伝達関数Ｐ（ｓ）の係数ａ０，ａ１，ａ２を算出したが、これに限定されない。例えば、モデル算出部３６０は、制御対象２００の特性を示すデータiｆと、制御対象２００の特性に基づいて算出した制御対象２００のモデルを示すデータｍｄとを含む学習データを用いた機械学習によって生成された学習済モデルに対して、周波数特性算出部３５０が新たに算出した制御対象２００の特性を示すデータiｆを適用し、制御対象２００のモデルを算出してもよい。
　図９は、ニューラルネットワークを用いて制御対象のモデルを算出する方法の説明図である。
　このとき、モデル算出部３６０は、機械学習のアルゴリズムとして公知のアルゴリズムを用いることができるので、例えば、ニューラルネットワークを用いてもよい。図９に示すように、本実施の形態の変更例に係るニューラルネットワークは、入力層、中間層、出力層を備えている。ここで、周波数特性算出部３５０が新たに算出した制御対象２００の特性を示すデータiｆのうち、制御対象２００の第１周波数ｆ１の特性を示すデータをiｆ１、制御対象２００の第２周波数ｆ２の特性を示すデータをiｆ２、…、制御対象２００の第ｎ周波数ｆｎの特性を示すデータをiｆｎとし、モデルパラメータをＡとする。なお、モデルパラメータＡは、伝達関数Ｐ（ｓ）の係数ａ０，ａ１，ａ２を含み、Ａ∋｛ａ０，ａ１，ａ２｝が成立する。
　変更例に係るニューラルネットワークでは、入力層に制御対象２００の第１周波数ｆ１から第ｎ周波数ｆｎまでの特性を示すデータiｆ１、iｆ２、…、iｆｎが入力され、出力層からモデルパラメータＡが出力される。また、入力層のノードと中間層のノードとの間の結合における重み係数Ｗ１、中間層のノードと出力層のノードとの間の結合における重み係数Ｗ２は、結合毎に設定できるようになっている。このため、中間層のノードには入力層のノードからの出力値とノードどうしの結合に対応する重み係数Ｗ１とが乗算された値が入力され、出力層のノードには中間層のノードからの出力値とノードどうしの結合に対応する重み係数Ｗ２とが乗算された値が入力される。なお、入力層・中間層・出力層の各ノードからの出力値は、入力値どうしの線形関数の値であってもよく、入力値を用いた非線形関数の値、例えばシグモイド関数の値であってもよい。
　よって、モデル算出部３６０は、制御対象２００の特性を示すデータiｆと当該特性に基づいて算出した制御対象２００のモデルを示すデータｍｄとを組み合わせたデータを学習データとし、複数種類の学習データを用いて上述した重み係数Ｗ１，Ｗ２を算出することで学習済モデルとしてのニューラルネットワークを生成することができる。なお、重み係数Ｗ１，Ｗ２を算出する方法としては、任意の方法を採用可能であり、例えば、公知の誤差伝搬法、勾配降下法を用いてもよい。この結果、モデル算出部３６０は、周波数特性算出部３５０が新たに算出した制御対象２００の特性を示すデータiｆを生成したニューラルネットワークに入力すれば当該特性に応じたモデルパラメータＡを出力することができる。
　このようにすることで、機械学習によって生成された学習済モデルに対して新たに算出した制御対象の特性を示すデータを適用して制御対象のモデルを算出しない制御システムよりも制御パラメータの値を精度良く算出できる。
　なお、本変更例では、モデル算出部３６０が学習データを用いて上述した重み係数Ｗ１，Ｗ２を算出することでニューラルネットワークを生成しているが、これに限定されず、例えば、制御装置３００とは異なる学習装置を新たに設けてニューラルネットワークを生成させてもよい。
　なお、本変更例では、モデル算出部３６０は、機械学習のアルゴリズムの一例としてのニューラルネットワークを用いて学習済モデルを生成したが、これに限定されない。例えば、Ｑ学習（Q-learning）、ＴＤ学習（TD-learning）等の強化学習（Reinforcement Learning）のアルゴリズムを用いて学習済モデルを生成してもよい。また、例えば、教師あり学習、教師なし学習の公知のアルゴリズムを用いて学習済モデルを生成してもよい。また、例えば、深層学習（Deep Learning）、遺伝的プログラミング、機能論理プログラミング、サポートベクターマシン等の公知の学習アルゴリズムを用いて学習済モデルを生成してもよい。
　なお、本実施の形態では、制御装置３００とは別体の検出器４００を設けているが、これに限定されず、検出器４００は制御装置３００と一体であってもよく、例えば、検出器４００が制御装置３００に含まれていてもよい。
　なお、制御部５１、主記憶部５２、外部記憶部５３、操作部５４、表示部５５、送受信部５６、内部バス５０を備える制御装置３００の処理を行う中心となる部分は、専用のシステムによらず、通常のコンピュータシステムを用いて実現可能である。例えば、前記の動作を実行するためのコンピュータプログラムを、コンピュータが読み取り可能な記録媒体、例えば、フレキシブルディスク、ＤＶＤ（Digital Versatile Disc）－ＲＯＭ（Read-Only Memory）等に格納して配布し、当該コンピュータプログラムをコンピュータにインストールすることにより、上記の処理を実行する制御装置３００を構成してもよい。また、通信ネットワーク上のサーバ装置が有する記憶装置に当該コンピュータプログラムを格納しておき、通常のコンピュータシステムがダウンロードすることで制御装置３００を構成してもよい。
　また、制御装置３００の機能を、ＯＳ（オペレーティングシステム）とアプリケーションプログラムの分担により実現する場合、または、ＯＳとアプリケーションプログラムとの協働により実現する場合には、アプリケーションプログラム部分のみを記録媒体や記憶装置に格納してもよい。
　また、搬送波にコンピュータプログラムを重畳し、通信ネットワークを介して提供することも可能である。例えば、通信ネットワーク上の掲示板（BBS, Bulletin Board System）にコンピュータプログラムを掲示し、ネットワークを介してコンピュータプログラムを提供してもよい。そして、このコンピュータプログラムを起動し、ＯＳの制御下で、他のアプリケーションプログラムと同様に実行することにより、前記の処理を実行してもよい。
　本開示は、本開示の広義の精神と範囲を逸脱することなく、様々な実施形態及び変形が可能とされるものである。また、上述した実施形態は、本開示を説明するためのものであり、本開示の範囲を限定するものではない。つまり、本開示の範囲は、実施形態ではなく、請求の範囲によって示される。そして、請求の範囲内及びそれと同等の開示の意義の範囲内で施される様々な変形が、本開示の範囲内とみなされる。
５０　内部バス、５１　制御部、５２　主記憶部、５３　外部記憶部、５４　操作部、５５　表示部、５６　送受信部、５９　制御プログラム、１００　制御システム、２００　制御対象、３００　制御装置、３１０　減算器、３２０　制御対象制御部、３３０　周波数成分抽出部、３４０　抽出結果判定部、３５０　周波数特性算出部、３６０　モデル算出部、３７０　変更判定部、３８０　制御パラメータ変更部、４００　検出器。

Claims (12)

  1. 　目標量と制御対象の制御量との偏差に基づいて、予め設定された制御パラメータの値を用いて前記制御対象の操作量を算出し、算出した前記操作量の信号を前記制御対象に出力して前記制御対象を制御する制御対象制御部と、
    　前記操作量の信号から第１周波数成分の信号および第２周波数成分の信号を抽出するとともに、前記制御量の信号から第１周波数成分の信号および第２周波数成分の信号を抽出する周波数成分抽出部と、
    　前記周波数成分抽出部が前記操作量の信号から抽出した前記第１周波数成分の信号と前記制御量の信号から抽出した前記第１周波数成分の信号とに基づく位相遅れとゲインとを含む前記制御対象の第１周波数の特性を算出するとともに、前記周波数成分抽出部が前記操作量の信号から抽出した前記第２周波数成分の信号と前記制御量の信号から抽出した前記第２周波数成分の信号とに基づく位相遅れおよびゲインを含む前記制御対象の第２周波数の特性を算出する周波数特性算出部と、
    　前記周波数特性算出部が算出した前記制御対象の前記第１周波数の特性と前記第２周波数の特性とに基づいて前記制御対象のモデルを算出するモデル算出部と、
    　前記モデル算出部が算出した前記制御対象のモデルに基づいて制御パラメータの値を算出して設定を変更する制御パラメータ変更部と、
    　を備える制御装置。
  2. 　前記周波数成分抽出部は、共振フィルタを用いて前記操作量または前記制御量の信号から第１周波数成分の信号と第２周波数成分の信号とを抽出し、
    　前記共振フィルタの減衰係数は、０．２よりも小さい、
    　請求項１に記載の制御装置。
  3. 　前記周波数成分抽出部が前記第１周波数の信号として抽出した信号の周波数が前記第１周波数と異なる場合、または、前記周波数成分抽出部が前記第２周波数成分の信号として抽出した信号の周波数が前記第２周波数と異なる場合、または、前記周波数成分抽出部が前記第１周波数の信号として抽出した信号の振幅が予め定めた閾値よりも小さい場合、または、前記周波数成分抽出部が前記第２周波数の信号として抽出した信号の振幅が前記閾値よりも小さい場合、前記周波数成分抽出部が信号の抽出に失敗したと判定する抽出結果判定部を更に備え、
    　前記周波数成分抽出部は、前記抽出結果判定部が信号の抽出に失敗したと判定したとき、信号の抽出を再度行う、
    　請求項１又は２に記載の制御装置。
  4. 　前記モデル算出部は、前記周波数特性算出部が算出した前記制御対象の前記第１周波数の特性と前記制御対象の伝達関数から算出される前記制御対象の前記第１周波数の特性との誤差が最小となり、且つ、前記周波数特性算出部が算出した前記制御対象の前記第２周波数の特性と前記伝達関数から算出される前記制御対象の前記第２周波数の特性との誤差が最小となる前記伝達関数の係数を算出することで、前記制御対象の数式モデルを算出する、
    　請求項１から３の何れか１項に記載の制御装置。
  5. 　前記モデル算出部は、前記制御対象の前記第１周波数の特性と前記第２周波数の特性とを示すデータと、前記制御対象の前記第１周波数の特性と前記第２周波数の特性とに基づいて算出した前記制御対象のモデルを示すデータとを含む学習データを用いた機械学習によって生成された学習済モデルに対して、前記周波数特性算出部が新たに算出した前記制御対象の前記第１周波数の特性と前記第２周波数の特性とを示すデータを適用し、前記制御対象のモデルを算出する、
    　請求項１から４の何れか１項に記載の制御装置。
  6. 　前記第１周波数および前記第２周波数のうちの少なくとも何れか一方の周波数は、少なくとも前記制御対象の既知の特徴を示す値と前記制御パラメータの値と前記制御装置の演算周期とのうちの何れか一つの値に基づいて予め定められている、
    　請求項１から５の何れか１項に記載の制御装置。
  7. 　前記第１周波数および前記第２周波数のうちの少なくとも何れか一方の周波数は、前記制御装置の使用者によって予め定められている、
    　請求項１から６の何れか１項に記載の制御装置。
  8. 　前記制御対象のモデルは、前記制御対象の伝達関数であり、
    　前記制御対象制御部の構造に対応して前記伝達関数の構造が予め定められている、
    　請求項１から７の何れか１項に記載の制御装置。
  9. 　前記制御量の信号に基づいて前記制御パラメータの設定を変更するか否かを判定する変更判定部をさらに備え、
    　前記制御パラメータ変更部は、前記変更判定部が前記制御パラメータの設定を変更すると判定したとき、前記モデル算出部が算出した前記制御対象のモデルに基づいて新たな制御パラメータの値を算出して設定を変更する、
    　請求項１から８の何れか１項に記載の制御装置。
  10. 　目標量と制御対象の制御量との偏差に基づいて、予め設定された制御パラメータの値を用いて前記制御対象の操作量を算出し、算出した前記操作量の信号を前記制御対象に出力して前記制御対象を制御する制御対象制御部と、
    　前記操作量の信号から第１周波数成分の信号および第２周波数成分の信号を抽出するとともに、前記制御量の信号から第１周波数成分の信号および第２周波数成分の信号を抽出する周波数成分抽出部と、
    　前記周波数成分抽出部が前記操作量の信号から抽出した前記第１周波数成分の信号と前記制御量の信号から抽出した前記第１周波数成分の信号とに基づく前記制御対象の第１周波数の特性を算出するとともに、前記周波数成分抽出部が前記操作量の信号から抽出した前記第２周波数成分の信号と前記制御量の信号から抽出した前記第２周波数成分の信号とに基づく前記制御対象の第２周波数の特性を算出する周波数特性算出部と、
    　前記周波数特性算出部が算出した前記制御対象の前記第１周波数の特性と前記第２周波数の特性とに基づいて前記制御対象のモデルを算出するモデル算出部と、
    　前記モデル算出部が算出した前記制御対象のモデルに基づいて制御パラメータの値を算出して設定を変更する制御パラメータ変更部と、
    　を備える制御システム。
  11. 　目標量と制御対象の制御量との偏差に基づいて、予め設定された制御パラメータの値を用いて前記制御対象の操作量を算出し、算出した前記操作量の信号を前記制御対象に出力して前記制御対象を制御する制御対象制御ステップと、
    　前記操作量の信号から第１周波数成分の信号および第２周波数成分の信号を抽出するとともに、前記制御量の信号から第１周波数成分の信号および第２周波数成分の信号を抽出する周波数成分抽出ステップと、
    　前記周波数成分抽出ステップで前記操作量の信号から抽出した前記第１周波数成分の信号と前記制御量の信号から抽出した前記第１周波数成分の信号とに基づく前記制御対象の第１周波数の特性を算出するとともに、前記周波数成分抽出ステップで前記操作量の信号から抽出した前記第２周波数成分の信号と前記制御量の信号から抽出した前記第２周波数成分の信号とに基づく前記制御対象の第２周波数の特性を算出する周波数特性算出ステップと、
    　前記周波数特性算出ステップで算出した前記制御対象の前記第１周波数の特性と前記第２周波数の特性とに基づいて前記制御対象のモデルを算出するモデル算出ステップと、
    　前記モデル算出ステップで算出した前記制御対象のモデルに基づいて制御パラメータの値を算出して設定を変更する制御パラメータ変更ステップと、
    　を行う制御方法。
  12. 　コンピュータを、
    　目標量と制御対象の制御量との偏差に基づいて、予め設定された制御パラメータの値を用いて前記制御対象の操作量を算出し、算出した前記操作量の信号を前記制御対象に出力して前記制御対象を制御する制御対象制御部、
    　前記操作量の信号から第１周波数成分の信号および第２周波数成分の信号を抽出するとともに、前記制御量の信号から第１周波数成分の信号および第２周波数成分の信号を抽出する周波数成分抽出部、
    　前記周波数成分抽出部が前記操作量の信号から抽出した前記第１周波数成分の信号と前記制御量の信号から抽出した前記第１周波数成分の信号とに基づく前記制御対象の第１周波数の特性を算出するとともに、前記周波数成分抽出部が前記操作量の信号から抽出した前記第２周波数成分の信号と前記制御量の信号から抽出した前記第２周波数成分の信号とに基づく前記制御対象の第２周波数の特性を算出する周波数特性算出部、
    　前記周波数特性算出部が算出した前記制御対象の前記第１周波数の特性と前記第２周波数の特性とに基づいて前記制御対象のモデルを算出するモデル算出部、
    　前記モデル算出部が算出した前記制御対象のモデルに基づいて制御パラメータの値を算出して設定を変更する制御パラメータ変更部、
    　として機能させるプログラム。
PCT/JP2020/017093 2020-04-20 2020-04-20 制御装置、制御システム、制御方法及びプログラム WO2021214839A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2020555933A JP6833133B1 (ja) 2020-04-20 2020-04-20 制御装置、制御システム、制御方法及びプログラム
PCT/JP2020/017093 WO2021214839A1 (ja) 2020-04-20 2020-04-20 制御装置、制御システム、制御方法及びプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2020/017093 WO2021214839A1 (ja) 2020-04-20 2020-04-20 制御装置、制御システム、制御方法及びプログラム

Publications (1)

Publication Number Publication Date
WO2021214839A1 true WO2021214839A1 (ja) 2021-10-28

Family

ID=74661678

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/017093 WO2021214839A1 (ja) 2020-04-20 2020-04-20 制御装置、制御システム、制御方法及びプログラム

Country Status (2)

Country Link
JP (1) JP6833133B1 (ja)
WO (1) WO2021214839A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7409343B2 (ja) * 2021-03-17 2024-01-09 横河電機株式会社 コントローラ、制御方法及び制御プログラム

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007219691A (ja) * 2006-02-15 2007-08-30 Seiko Epson Corp Ｐｉｄ制御装置および制御パラメータ更新方法
JP2017167663A (ja) * 2016-03-14 2017-09-21 株式会社東芝 プロセス診断装置、プロセス診断方法及びコンピュータプログラム
JP2018128735A (ja) * 2017-02-06 2018-08-16 ファナック株式会社 サーボ制御装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007219691A (ja) * 2006-02-15 2007-08-30 Seiko Epson Corp Ｐｉｄ制御装置および制御パラメータ更新方法
JP2017167663A (ja) * 2016-03-14 2017-09-21 株式会社東芝 プロセス診断装置、プロセス診断方法及びコンピュータプログラム
JP2018128735A (ja) * 2017-02-06 2018-08-16 ファナック株式会社 サーボ制御装置

Also Published As

Publication number Publication date
JPWO2021214839A1 (ja) 2021-10-28
JP6833133B1 (ja) 2021-02-24

Similar Documents

Publication Publication Date Title
US9720387B2 (en) Scaling and parameterizing a controller
JP4223894B2 (ja) Ｐｉｄパラメータ調整装置
US6745088B2 (en) Multi-variable matrix process control
EP3379349B1 (en) Control device, control program, and control system
JP2009175917A (ja) 制御パラメータ調整方法および制御パラメータ調整プログラム
WO2021214839A1 (ja) 制御装置、制御システム、制御方法及びプログラム
JP2021535500A (ja) マスフローコントローラおよびコントローラアルゴリズム
Saab Development of multivariable PID controller gains in presence of measurement noise
JP4815391B2 (ja) モデルパラメータ推定演算装置及び方法、モデルパラメータ推定演算処理プログラム並びにそれを記録した記録媒体
EP3575893A1 (en) Control parameter computation method and control parameter computation device
JP5125754B2 (ja) Ｐｉｄコントローラのチューニング装置、ｐｉｄコントローラのチューニング用プログラムおよびｐｉｄコントローラのチューニング方法
JP5749041B2 (ja) アクティブフィードバック制御装置及びプログラム
US20220276621A1 (en) Control device for plant and controlling method of the same
EP4063975A1 (en) Rst smith predictor
JP6848710B2 (ja) プラント制御調整装置及び方法
WO2019127974A1 (zh) 用于优化风力发电机组的载荷的方法和装置
JP4485896B2 (ja) Ｐｉｄ制御装置
US20230176532A1 (en) Control assist device, control device, and control assist method
JP2023031898A (ja) 情報処理装置
US11726470B2 (en) Method and device for determining a measure of quality continuously
JP3012713B2 (ja) 波形同定システム
Hunnekens Performance optimization of hybrid controllers for linear motion systems
CN118090102A (zh) 基于模糊神经网络的实时混合试验加载装置时滞补偿方法
JP2013109519A (ja) 設計装置及び設計プログラム
JPH0784611A (ja) コントローラ

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2020555933

Country of ref document: JP

Kind code of ref document: A

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

Ref document number: 20932678

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20932678

Country of ref document: EP

Kind code of ref document: A1