JP2010266967A - Device and program for adjustment of pid - Google Patents

Device and program for adjustment of pid Download PDF

Info

Publication number
JP2010266967A
JP2010266967A JP2009116015A JP2009116015A JP2010266967A JP 2010266967 A JP2010266967 A JP 2010266967A JP 2009116015 A JP2009116015 A JP 2009116015A JP 2009116015 A JP2009116015 A JP 2009116015A JP 2010266967 A JP2010266967 A JP 2010266967A
Authority
JP
Japan
Prior art keywords
pid
transfer function
parameter
control amount
obtaining
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2009116015A
Other languages
Japanese (ja)
Other versions
JP5549112B2 (en
Inventor
Chikashi Nakazawa
親志 中沢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuji Electric Co Ltd
Original Assignee
Fuji Electric Systems Co Ltd
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 Fuji Electric Systems Co Ltd filed Critical Fuji Electric Systems Co Ltd
Priority to JP2009116015A priority Critical patent/JP5549112B2/en
Publication of JP2010266967A publication Critical patent/JP2010266967A/en
Application granted granted Critical
Publication of JP5549112B2 publication Critical patent/JP5549112B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a PID (Proportional Integral Derivate) adjustment device 30 or the like allowing widely spread of an adjustment method or the like of a PID adjuster to on-site engineers, and allowing a significant reduction of adjustment labor and a time. <P>SOLUTION: A transfer function parameter acquisition unit 32 obtains a parameter of a plurality of transfer functions representing a control target 10 based on time series data of an operation amount MV input to the control target 10 and time series data of a control amount PV output from the control target 10 according to the operation amount MV. An adaptability acquisition unit 34 obtains an adaptation rate between an estimation control amount and an actual control amount output by each transfer function having the parameter. An optimal transfer function selection unit 36 selects the transfer function having optimal adaptability based on the adaptation rate. An optimal PID parameter search unit 38 searches for a PID parameter of the PID adjuster 20 such that a predetermined evaluation index is satisfied from results of a disturbance response to disturbance d and a setting value response to a setting value SV with the selected transfer function (11) or the like as a target. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

本発明は、制御対象を閉ループにて制御するPID調節器のPIDパラメータを求めるPID調整装置等に関する。   The present invention relates to a PID adjusting device for obtaining a PID parameter of a PID adjuster that controls a controlled object in a closed loop.

PID(Proportional Integral Derivate)制御は、石油、化学、鉄鋼、紙パルプ、水環境等のいわゆるプロセス産業における制御方式において、90%以上を占めるといわれている基盤的制御技術であり、PID制御は、実用上不可欠な制御方式となっている。従来、制御対象の同定方法とその同定方法から得られたモデルである伝達関数をもとにPID調節器の制御パラメータを決定するPID調節器の調整方法、調整に用いるプログラムおよびPID調節装置に関しては、シミュレーションにより試行錯誤的に実施する方法が提案されてきた(特許文献1参照)。特許文献1に記載されているように、制御器シミュレーション部は作業者との対話処理に従って制御パラメータ(PIDパラメータ)の値を設定し、作業者は当該PIDパラメータが適切なものである場合、制御器シミュレーション部へシミュレーションの終了指示を発行していた。   PID (Proportional Integral Derivate) control is a fundamental control technology that is said to occupy over 90% in so-called process industries such as petroleum, chemical, steel, paper pulp, and water environment. It is an indispensable control method for practical use. Conventionally, regarding a method for identifying a control object, a method for adjusting a PID regulator that determines a control parameter of the PID regulator based on a transfer function that is a model obtained from the identification method, a program used for adjustment, and a PID regulator A method of performing trial and error by simulation has been proposed (see Patent Document 1). As described in Patent Document 1, the controller simulation unit sets the value of the control parameter (PID parameter) according to the dialogue processing with the worker, and the worker controls the PID parameter when the PID parameter is appropriate. A simulation end instruction was issued to the vessel simulator.

上述のように、シミュレーションにより試行錯誤的にPIDパラメータを調整する場合、作業者は制御対象に関する経験・ノウハウを必要とする。このため、PID調節器の調整方法およびPID調節装置等に関して、シミュレーションにより試行錯誤的に実施する方法は広く現場のエンジニアに普及させるには困難であるという問題があった。さらに、シミュレーションにより試行錯誤的に実施する方法は、調整労力と時間とを極めて要するという問題があった。   As described above, when adjusting the PID parameters by trial and error by simulation, the operator needs experience and know-how regarding the controlled object. For this reason, there has been a problem that the method of performing the trial and error by simulation with respect to the adjustment method of the PID controller and the PID adjustment device is difficult to spread widely to engineers in the field. Furthermore, the method of performing trial and error by simulation has a problem that adjustment labor and time are extremely required.

そこで、本発明の目的は、上記問題を解決するためになされたものであり、PID調節器の調整方法およびPID調節装置等に関して、広く現場のエンジニアに普及させることが可能であり、且つ調整労力と時間とを著しく短縮することができるPID調整装置等を提供することにある。   Therefore, an object of the present invention is to solve the above problems, and it is possible to widely disseminate the adjustment method of the PID adjuster, the PID adjuster and the like to engineers on the site, and the adjustment effort. It is an object of the present invention to provide a PID adjusting device and the like that can significantly reduce the time.

この発明のPID調整装置は、制御対象を閉ループにて制御するPID調節器のPIDパラメータを求めるPID調整装置であって、記録部に記録された制御対象に入力する操作量の時系列データと該操作量に応じて制御対象から出力される制御量の時系列データとに基づき、制御対象を表現する複数の伝達関数のパラメータを求める伝達関数パラメータ取得手段と、前記伝達関数パラメータ取得手段により求められたパラメータを有する各伝達関数が出力した推定制御量と実際の制御量との適合性を求める適合性取得手段と、前記適合性取得手段により求められた適合性に基づき、最良の適合性を有する伝達関数を選択する又は選択させる最良伝達関数選択手段と、前記最良伝達関数選択手段により選択された又は選択させた伝達関数を対象として、設定値応答及び外乱応答の結果から予め決めておいた評価指標を満たすべくPID調節器のPIDパラメータを探索する最適PIDパラメータ探索手段とを備えたことを特徴とする。   A PID adjusting device according to the present invention is a PID adjusting device that obtains a PID parameter of a PID adjuster that controls a controlled object in a closed loop, the time series data of the operation amount input to the controlled object recorded in a recording unit, and the PID adjusting device Based on time series data of the controlled variable output from the controlled object according to the manipulated variable, transfer function parameter acquiring means for determining a plurality of transfer function parameters representing the controlled object, and obtained by the transfer function parameter acquiring means A fitness acquisition means for determining the compatibility between the estimated control amount output by each transfer function having the parameters and the actual control amount, and the best compatibility based on the compatibility obtained by the compatibility acquisition means. The best transfer function selection means for selecting or selecting a transfer function is paired with the transfer function selected or selected by the best transfer function selection means. As, characterized in that a optimum PID parameter search means for searching the PID parameters of the PID controller to meet the evaluation index determined in advance from results of the set value response and disturbance response.

ここで、この発明のPID調整装置において、前記最適PIDパラメータ探索手段が用いる設定値応答及び外乱応答の結果は、制御量に白色ノイズを注入した状態の設定値応答及び外乱応答の結果とすることができる。   Here, in the PID adjusting device of the present invention, the result of the set value response and the disturbance response used by the optimum PID parameter search means is the result of the set value response and the disturbance response in a state where white noise is injected into the control amount. Can do.

ここで、この発明のPID調整装置において、前記最適PIDパラメータ探索手段が用いる評価指標は、制御量とステップ状の設定値との差の2乗を時間で積分した指標とすることができる。   Here, in the PID adjustment apparatus of the present invention, the evaluation index used by the optimum PID parameter search means can be an index obtained by integrating the square of the difference between the control amount and the step-like set value with time.

ここで、この発明のPID調整装置において、前記適合性取得手段が用いる適合性は、推定制御量、実際の制御量及び実際の制御量の平均値に基づく所定の適合率とすることができる。   Here, in the PID adjustment apparatus of the present invention, the suitability used by the suitability acquisition unit can be a predetermined suitability ratio based on an estimated control amount, an actual control amount, and an average value of the actual control amounts.

この発明のPID調整プログラムは、制御対象を閉ループにて制御するPID調節器のPIDパラメータを求めるPID調整プログラムであって、PID調整御プログラムを実行するPID調整装置のコンピュータに、記録部に記録された制御対象に入力する操作量の時系列データと該操作量に応じて制御対象から出力される制御量の時系列データとに基づき、制御対象を表現する複数の伝達関数のパラメータを求める伝達関数パラメータ取得ステップ、前記伝達関数パラメータ取得ステップで求められたパラメータを有する各伝達関数が出力した推定制御量と実際の制御量との適合性を求める適合性取得ステップ、前記適合性取得ステップで求められた適合性に基づき、最良の適合性を有する伝達関数を選択する又は選択させる最良伝達関数選択ステップ、前記最良伝達関数選択ステップで選択された又は選択させた伝達関数を対象として、設定値応答及び外乱応答の結果から予め決めておいた評価指標を満たすべくPID調節器のPIDパラメータを探索する最適PIDパラメータ探索ステップを実行させるためのPID調整プログラムである。   The PID adjustment program according to the present invention is a PID adjustment program for obtaining a PID parameter of a PID adjuster that controls a controlled object in a closed loop, and is recorded in a recording unit in a computer of the PID adjustment apparatus that executes the PID adjustment control program. A transfer function for obtaining parameters of a plurality of transfer functions expressing the control object based on the time series data of the operation amount input to the control object and the time series data of the control amount output from the control object according to the operation amount A parameter obtaining step, a suitability obtaining step for obtaining suitability between an estimated control amount output by each transfer function having the parameters obtained in the transfer function parameter obtaining step and an actual control amount, and obtained in the suitability obtaining step. The best transfer function that selects or selects the transfer function with the best fit based on the best fit Searching for the PID parameter of the PID controller to satisfy the evaluation index determined in advance from the result of the set value response and the disturbance response for the transfer function selected or selected in the selection step and the best transfer function selection step It is a PID adjustment program for executing the optimal PID parameter search step.

ここで、この発明のPID調整プログラムにおいて、前記最適PIDパラメータ探索ステップで用いる設定値応答及び外乱応答の結果は、制御量に白色ノイズを注入した状態の設定値応答及び外乱応答の結果とすることができる。   Here, in the PID adjustment program of the present invention, the result of the set value response and the disturbance response used in the optimum PID parameter search step is the result of the set value response and the disturbance response in a state where white noise is injected into the control amount. Can do.

本発明のPID調整装置等によれば、PID調整装置の最適PIDパラメータ探索部が、最良伝達関数選択部により選択されたまたは選択させた伝達関数を対象として、設定値SVに対する設定値応答および外乱dに対する外乱応答の結果から予め決めておいた評価指標を満たすように、PID調節器のPIDパラメータを探索することができる。この結果、作業者が制御対象に関する経験・ノウハウに基づき、PIDパラメータの決定を試行錯誤的に行う必要がなくなるため、PID調節器の調整方法およびPID調節装置等に関して、広く現場のエンジニアに普及させることが可能となり、加えて調整労力と時間とを著しく短縮することができるという効果がある。   According to the PID adjustment device or the like of the present invention, the optimum PID parameter search unit of the PID adjustment device targets the transfer function selected or selected by the best transfer function selection unit, and the set value response and disturbance to the set value SV. The PID parameter of the PID adjuster can be searched so as to satisfy a predetermined evaluation index from the result of the disturbance response to d. As a result, it is not necessary for the operator to determine the PID parameters based on experience and know-how related to the control target, so that the PID controller adjustment method, PID adjustment device, and the like are widely disseminated to on-site engineers. In addition, the adjustment labor and time can be remarkably shortened.

本発明の実施例1におけるPID調整装置を示す図である。It is a figure which shows the PID adjustment apparatus in Example 1 of this invention. 伝達関数パラメータ取得部32の構成を示す図である。4 is a diagram illustrating a configuration of a transfer function parameter acquisition unit 32. FIG. 最良伝達関数選択部36が適合性FIT(%)に基づき、最良の適合性を有する伝達関数を選択するまたは選択させる状態を示す図である。It is a figure which shows the state in which the best transfer function selection part 36 selects or makes it select the transfer function which has the best fitness based on fitness FIT (%). 最適PIDパラメータ探索部38の機能を説明するための図である。4 is a diagram for explaining a function of an optimum PID parameter search unit 38. FIG. 図4で説明した設定値SV、外乱dおよび制御量PVのタイムチャートである。5 is a time chart of the set value SV, disturbance d, and control amount PV described in FIG. PID調整装置30の各機能をすべてソフトウェアにより実現した場合の処理の流れを示すフローチャートである。It is a flowchart which shows the flow of a process when each function of the PID adjustment apparatus 30 is implement | achieved by software. 本発明の実施例2におけるPID調整装置を示す図である。It is a figure which shows the PID adjustment apparatus in Example 2 of this invention. 実施例2における最適PIDパラメータ探索部38の機能を説明するための図である。It is a figure for demonstrating the function of the optimal PID parameter search part 38 in Example 2. FIG. 図8で説明した設定値SV、外乱dおよび制御量PVのタイムチャートを示す図である。It is a figure which shows the time chart of setting value SV demonstrated in FIG. 8, disturbance d, and control amount PV. 本発明のPID調整プログラムを実行するPID調整装置30内のコンピュータの内部回路40を示すブロック図である。It is a block diagram which shows the internal circuit 40 of the computer in the PID adjustment apparatus 30 which executes the PID adjustment program of this invention.

以下、各実施例について図面を参照して詳細に説明する。   Hereinafter, each embodiment will be described in detail with reference to the drawings.

図1は、本発明の実施例1におけるPID調整装置を示す。図1において、符号10は制御対象、20は制御対象10を閉ループにて制御するPID調節器、30はPID調節器20のPIDパラメータを求めるPID調整装置である。図1に示されるように、PID調節器20から制御対象10へ操作量MVが入力され、当該操作量MVと外乱dとに応じて制御対象10から制御量PVが出力される。PID調整装置30は、図1に示されるように、伝達関数パラメータ取得部(伝達関数パラメータ取得手段)32、適合性取得部(適合性取得手段)34、最良伝達関数選択部(最良伝達関数選択手段)36および最適PIDパラメータ探索部(最適PIDパラメータ探索手段)38の各機能を有している。以下、図1および他の図面を適宜参照しながら説明する。   FIG. 1 shows a PID adjustment apparatus according to Embodiment 1 of the present invention. In FIG. 1, reference numeral 10 denotes a control object, 20 denotes a PID adjuster that controls the control object 10 in a closed loop, and 30 denotes a PID adjustment device that obtains a PID parameter of the PID adjuster 20. As illustrated in FIG. 1, the operation amount MV is input from the PID adjuster 20 to the control target 10, and the control amount PV is output from the control target 10 according to the operation amount MV and the disturbance d. As shown in FIG. 1, the PID adjustment device 30 includes a transfer function parameter acquisition unit (transfer function parameter acquisition unit) 32, a suitability acquisition unit (compatibility acquisition unit) 34, and a best transfer function selection unit (best transfer function selection). Means) 36 and an optimum PID parameter search unit (optimum PID parameter search means) 38. Hereinafter, description will be made with reference to FIG. 1 and other drawings as appropriate.

図1に示されるように、伝達関数パラメータ取得部32は、記録部(後述するディスク等の記録装置48)に記録された制御対象10に入力する操作量MVの時系列データ(点線で示される方のMV)とこの操作量MVに応じて制御対象10から出力される制御量PVの時系列データ(点線で示される方のPV)とに基づき、制御対象10を表現する複数の伝達関数のパラメータを求める。   As shown in FIG. 1, the transfer function parameter acquisition unit 32 includes time series data (indicated by dotted lines) of the operation amount MV input to the control object 10 recorded in the recording unit (recording device 48 such as a disk to be described later). MV) and the time series data (PV indicated by a dotted line) of the control amount PV output from the control object 10 in accordance with the manipulated variable MV, a plurality of transfer functions representing the control object 10 Find the parameters.

図2は、伝達関数パラメータ取得部32の構成を示す。図2において、符号11ないし17はオフラインにて予め構造を決めておいた制御対象10を表現する複数の伝達関数である。図2に示されるように、伝達関数11は一次系(または一次おくれ要素)+むだ時間(またはむだ時間要素)、伝達関数12は二次系(または二次おくれ要素)+むだ時間、伝達関数13は一次系+積分系(積分要素)+むだ時間、伝達関数14は二次系+積分系、伝達関数15は二次系(振動系:二次おくれ要素の減衰係数ζが0<ζ<1であり、ステップ応答が振動的)+むだ時間、伝達関数16は積分系+むだ時間、伝達関数17は非最小位相系である。以上のように、図2では7種類の伝達関数11ないし17の例を示したが、他の伝達関数を用いてもよいことは勿論である。伝達関数パラメータ取得部32は、操作量MVの時系列データと操作量MVに応じて制御対象10から出力される制御量PVの時系列データとから、各伝達関数11ないし17のパラメータである、ゲイン定数G、時定数T、T 、T 、T 、むだ時間L、減衰係数ζ、固有角周波数ωを求める。これらの各伝達関数11ないし17のパラメータは、以下の式1で示される目的関数Jの最小値を求めることにより得ることができる。 FIG. 2 shows the configuration of the transfer function parameter acquisition unit 32. In FIG. 2, reference numerals 11 to 17 denote a plurality of transfer functions that represent the control object 10 whose structure is previously determined offline. As shown in FIG. 2, the transfer function 11 is a first-order system (or first-order element) + dead time (or time-delay element), the transfer function 12 is a second-order system (or second-order element) + dead time, and a transfer function. 13 is a primary system + integration system (integration element) + dead time, transfer function 14 is a secondary system + integration system, transfer function 15 is a secondary system (vibration system: damping coefficient ζ of the secondary sway element is 0 <ζ < 1 and step response is oscillatory) + dead time, transfer function 16 is integral system + dead time, and transfer function 17 is a non-minimum phase system. As described above, FIG. 2 shows examples of seven types of transfer functions 11 to 17, but other transfer functions may be used as a matter of course. The transfer function parameter acquisition unit 32 is a parameter of each transfer function 11 to 17 from the time series data of the operation amount MV and the time series data of the control amount PV output from the control target 10 according to the operation amount MV. gain constant G *, the time constant T *, T 1 *, T 2 *, T 3 *, the dead time L *, attenuation coefficient zeta *, obtains the natural angular frequency omega *. These parameters of the transfer functions 11 to 17 can be obtained by obtaining the minimum value of the objective function J shown by the following expression 1.

Figure 2010266967
Figure 2010266967

ここで、PVは推定された制御量である。式1は(無制約)非線形最適化問題となるため、例えばニュートン法、準ニュートン法、最急降下法等を適用することにより式1の最小値Jを求めることができる。 Here, PV * is an estimated control amount. Since Equation 1 is an (unconstrained) nonlinear optimization problem, the minimum value J of Equation 1 can be obtained by applying, for example, the Newton method, the quasi-Newton method, the steepest descent method, or the like.

図1に示されるように、適合性取得部34は、伝達関数パラメータ取得部32により求められたパラメータを有する各伝達関数が出力した推定制御量(点線で示されるPV)と実際の制御量(点線で示される方のPV)との適合性を求める。伝達関数パラメータ取得部32により求められたk番目の推定制御量をPV(k)、k番目の実際の制御量をPV(k)、制御量PV(k)の平均値を、 As shown in FIG. 1, the suitability acquisition unit 34 includes an estimated control amount (PV * indicated by a dotted line) output by each transfer function having parameters obtained by the transfer function parameter acquisition unit 32 and an actual control amount. The compatibility with the PV indicated by the dotted line is obtained. The k-th estimated control amount obtained by the transfer function parameter acquisition unit 32 is PV * (k), the k-th actual control amount is PV (k), and the average value of the control amounts PV (k) is

Figure 2010266967
Figure 2010266967

とすると、適合性は、以下の式2で表される適合率FIT(%)により求めることができる。 Then, the conformity can be obtained from the conformity rate FIT (%) expressed by the following formula 2.

Figure 2010266967
Figure 2010266967

ここで、Nは時系列データのデータ数を表す。本明細書では、推定制御量PVと実際の制御量PVとの2つのデータ比較に上記適合率FIT(%)を採用した。しかし、適合率FIT(%)の代わりに2つのデータの相関を用いることも可能である。 Here, N represents the number of data of time series data. In the present specification, the precision FIT (%) is used for comparing two data of the estimated control amount PV * and the actual control amount PV. However, it is also possible to use the correlation of two data instead of the precision FIT (%).

図1に示されるように、最良伝達関数選択部36は、適合性取得部34により求められた適合性FIT(%)(点線で示されるFIT)に基づき、最良の適合性を有する伝達関数を選択するまたは選択させる。図3は、最良伝達関数選択部36が適合性FIT(%)に基づき、最良の適合性を有する伝達関数を選択するまたは選択させる状態を示す。図3で図2と同じ符号を付した箇所は同じ要素を示すため説明は省略する。図3に表形式で、伝達関数欄、適合率欄および選択欄と、各伝達関数11ないし17とが示されている。上記選択するまたは選択させる状態は、当該表形式により後述するディスプレイ等の表示装置44に表示することができる。最良伝達関数選択部36は、図3に示すように、各伝達関数11等に対して計算された適合率FIT(%)を各伝達関数欄に対応する適合率欄に表示し、最も適合率FIT(%)が良い伝達関数を選択欄にチェック等を表示することにより選択する。図3では、伝達関数11を選択した状態が表示されている。選択は最良伝達関数選択部36が適合率FIT(%)の値により自動的に選択するが、手動にて任意の伝達関数11等を選択することも可能である。この場合、後述するマウス等の入力操作部51を用いて選択欄中の所望の矩形をクリックすることにより選択することができる。   As shown in FIG. 1, the best transfer function selection unit 36 selects a transfer function having the best fitness based on the fitness FIT (%) (FIT indicated by a dotted line) obtained by the fitness acquisition unit 34. Select or have it selected. FIG. 3 shows a state where the best transfer function selection unit 36 selects or selects a transfer function having the best fitness based on the fitness FIT (%). In FIG. 3, the same reference numerals as those in FIG. FIG. 3 shows, in tabular form, the transfer function column, the matching rate column, the selection column, and the transfer functions 11 to 17. The state to be selected or selected can be displayed on the display device 44 such as a display described later in the table format. As shown in FIG. 3, the best transfer function selection unit 36 displays the relevance ratio FIT (%) calculated for each transfer function 11 and the like in the relevance ratio column corresponding to each transfer function field, A transfer function with a good FIT (%) is selected by displaying a check or the like in the selection column. In FIG. 3, a state in which the transfer function 11 is selected is displayed. The selection is automatically selected by the best transfer function selection unit 36 based on the value of the precision FIT (%), but any transfer function 11 or the like can be manually selected. In this case, selection can be made by clicking a desired rectangle in the selection column using an input operation unit 51 such as a mouse described later.

図1に示されるように、最適PIDパラメータ探索部38は、最良伝達関数選択部36により選択されたまたは選択させた伝達関数11等(点線で示される11等)を対象として、設定値(後述するSV)に対する設定値応答および外乱(後述するd)に対する外乱応答の結果から予め決めておいた評価指標を満たすように、PID調節器20のPIDパラメータを探索する。   As shown in FIG. 1, the optimum PID parameter search unit 38 sets a set value (described later) for the transfer function 11 or the like (11 indicated by a dotted line) selected or selected by the best transfer function selection unit 36. The PID parameter of the PID adjuster 20 is searched so as to satisfy a predetermined evaluation index from the result of the set value response to SV) and the disturbance response result to disturbance (d described later).

図4は、最適PIDパラメータ探索部38の機能を説明するための図である。図4で図1と同じ符号を付した箇所は同じ要素を示すため説明は省略する。図4に示されるように、初めに設定値SVにステップ状の信号を入力し、ステップ応答が収束した後、外乱dにステップ上の信号を入力し、出力である制御量PVから評価指標を計算する。図5は、図4で説明した設定値SV、外乱dおよび制御量PVのタイムチャートを示す。図5(A)は設定値SV、図5(B)は外乱d、図5(C)は制御量PV、横軸はいずれも時間であり、時間tは操作量MVに外乱dが加えられた時間を示す。図4に示される評価指標の計算では、以下の式3で表されるISE(Integral of Squared Error)を用いて評価指標を計算する。 FIG. 4 is a diagram for explaining the function of the optimum PID parameter search unit 38. In FIG. 4, portions denoted by the same reference numerals as those in FIG. As shown in FIG. 4, first, a step-like signal is input to the set value SV, and after the step response has converged, a signal on the step is input to the disturbance d, and the evaluation index is calculated from the control amount PV as the output. calculate. FIG. 5 shows a time chart of the set value SV, the disturbance d, and the control amount PV described in FIG. 5A is the set value SV, FIG. 5B is the disturbance d, FIG. 5C is the control amount PV, the horizontal axis is time, and the time t d is the addition of the disturbance d to the manipulated variable MV. Time given. In the calculation of the evaluation index shown in FIG. 4, the evaluation index is calculated using ISE (Integral of Squared Error) expressed by the following Expression 3.

Figure 2010266967
Figure 2010266967

ここで、e(t)=PV−SVである。即ち、評価指標は制御量PVとステップ状の設定値SVとの差の2乗を時間で0から無限大まで積分した指標である。図5では、e(t)につき、制御量PVの立ち上がり時間の部分と外乱応答が加えられた部分とを示してある。評価指標である式3は、例えば、パーティクル・スォーム・オプティマイゼイション(Particle Swarm Optimization:PSO)等の探索手法により、PID調節器20のパラメータである、比例ゲインK、積分時間T、微分時間Tを求めることができる。 Here, e (t) = PV-SV. That is, the evaluation index is an index obtained by integrating the square of the difference between the control amount PV and the step-like set value SV from 0 to infinity over time. FIG. 5 shows a rise time portion of the control amount PV and a portion to which a disturbance response is added for e (t). Equation 3, which is an evaluation index, is, for example, a proportional gain K P , integration time T i , differentiation, which is a parameter of the PID controller 20 by a search method such as Particle Swarm Optimization (PSO). You can determine the time T D.

PSOは、簡単化した社会モデルのシミュレーションを通して開発されたモダン・ヒューリスティック(Modern heuristic : MH)手法の1つであり、鳥の群れの動きを連続変数の2次元空間で表現することを通して開発された。   PSO is one of the modern heuristic (MH) methods developed through the simulation of a simplified social model, and was developed through expressing the movement of a flock of birds in a two-dimensional space of continuous variables. .

PSOについては、J.Kennedy and R.Eberhartによる“Partic1e Swarm Optimization”(Proc. of IEEE International Conference on Neural Networks, Vol. IV,pp.1942-1948, Perth, Australia, 1995.)または吉田・福山他による「電圧信頼度を考慮したParticle Swarm Optimizationによる電圧無効電力制御方式の検討」(電気学会論文誌B, 119巻12号, 1999年12月)、特開2000−116003「電圧無効電力制御方法」、特開2002−51464「配電系統における状態推定法」等を参照されたい。   Regarding PSO, “Partic1e Swarm Optimization” by J. Kennedy and R. Eberhart (Proc. Of IEEE International Conference on Neural Networks, Vol. IV, pp.1942-1948, Perth, Australia, 1995.), Yoshida, Fukuyama, etc. "Study of Voltage Reactive Power Control Method by Particle Swarm Optimization Considering Voltage Reliability" (The Institute of Electrical Engineers of Japan B, Vol.119, No.12, December 1999), JP 2000-11603 "Voltage Reactive Power Control Method" JP 2002-51464 “State estimation method in distribution system” and the like.

PSOでは、各エージェント(PSOでは、群れの個々をエージェントと呼ぶ。)の位置(状態量)をx、y座標で表し、各エージェントの速度をvx(x方向の速度)、vy(y方向の速度)で表現する。各エージェントの位置および速度の情報から、次時点における各エージェントの位置を更新することができる。   In PSO, the position (state quantity) of each agent (in PSO, each group is called an agent) is represented by x and y coordinates, and the velocity of each agent is represented by vx (velocity in x direction) and vy (in y direction). (Speed). From the position and speed information of each agent, the position of each agent at the next time point can be updated.

上記エージェントの概念に基づき、鳥の群れ全体が何らかの目的関数を最適化するような行動をとると考えると、以下のような最適化が考えられる。即ち、各エージェントは、各々の探索における目的関数の従前における最も評価の高い最良値pbestと、その時の位置(状態量)を示すx、y座標とを覚えている。各エージェントは、各エージェントが保有する最良値pbestの集団の中で最も評価の高い値、即ち、当該集団の従前の目的関数の最良値gbestを共有している。各エージェントは、現在の位置(x、y座標)および現在の速度(vx、vy)を持っている。各エージェントは、現在の位置(x、y座標)および速度(vx、vy)と、最良値pbest、gbestとの間の距離に応じて、最良値pbest、gbestの存在する位置へ進むように方向を変更しようとする。   Based on the above-mentioned agent concept, the following optimization can be considered if the whole flock of birds takes an action that optimizes some objective function. That is, each agent remembers the best value pbest with the highest evaluation before the objective function in each search and the x and y coordinates indicating the position (state quantity) at that time. Each agent shares the highest value among the group of best values pbest held by each agent, that is, the best value gbest of the previous objective function of the group. Each agent has a current position (x, y coordinates) and a current velocity (vx, vy). Each agent moves to the position where the best values pbest and gbest exist according to the distance between the current position (x and y coordinates) and velocity (vx and vy) and the best values pbest and gbest. Try to change.

上記変更しようとする行動は速度で表現される。各エージェントの速度は、現在の速度と最良値pbest、gbestとを用いて、式4のように修正することができる。   The action to be changed is expressed by speed. The speed of each agent can be corrected as shown in Equation 4 using the current speed and the best values pbest and gbest.

k+1=w×V +c×rand()×(pbest−s
+c×rand()×(gbest−s ) (4)
V i k + 1 = w × V i k + c 1 × rand () × (pbest i −s i k )
+ C 2 × rand () × (gbest i −s i k ) (4)

但し、V は、エージェントiの速度、rand()は、0〜1までの一様乱数、skは、エージェントiの探索k回目の位置(探索点)、pbestは、エージェントiの最良値p b e s t、wはエージェント速度に対する重み関数、c1、は、各項に対する重み係数である。 Where V i k is the speed of agent i, rand () is a uniform random number from 0 to 1, s i k is the position (search point) of agent i in search k, and pbest i is agent i best value pbest, w is the weighting function for the agent velocity, c 1, c 2 are weighting coefficients for each term.

上記式4を用いることにより、各エージェントiの従前の最良値pbestおよび集団の最良値gbestに確率的に近づくような速度V k+1を求めることができ、速度V k+1により、各エージェントiの現在の位置(探索点)s を式5のように修正することができる。 By using the above equation 4, it is possible to obtain a speed V i k + 1 that probabilistically approaches the previous best value pbest i of each agent i and the best value gbest of the group, and each agent i is determined by the speed V i k + 1. The current position (search point) s i k can be corrected as shown in Equation 5.

k+1=s +V k+1 (5) s i k + 1 = s i k + V i k + 1 (5)

具休的には、PSOのエージェントの状態変数は、PIDパラメータである比例ゲインK、積分時間T、微分時間Tとする。 In practice, the state variables of the PSO agent are the proportional gain K P , the integration time T i , and the differentiation time T D which are PID parameters.

次に、エージェントの初期値として、エージェントの位置sおよび速度vを設定する。なお、エージェントの数iは事前に設定しておく。 Next, the agent position s i and the velocity v i are set as initial values of the agent. Note that the number i of agents is set in advance.

続いて、式4および式5の計算を繰り返すことにより、最適解を探索する。ここで、繰り返し回数kは事前に設定しておく。最適解を探索する際に、繰返しの各ステップ毎に、各エージェントiの評価と、各探索点の最良値pbestおよび集団の最良値gbestの更新を行う。 Subsequently, the optimum solution is searched by repeating the calculations of Equations 4 and 5. Here, the number of repetitions k is set in advance. When searching for the optimum solution, the evaluation of each agent i and the update of the best value pbest i of each search point and the best value gbest of the group are performed at each step of the iteration.

ここで、各エージェントiの評価は、以下の手順で行うことができる。まず、エージェントiの状態変数sを、PIDパラメータである比例ゲインK、積分時間T、微分時間Tに設定する、次に、設定値SVに対する設定値応答および外乱dに対する外乱応答の結果から、評価指標(式3)を評価する。すべてのエージェントの評価結果に基づいて、各探索点の最良値pbestおよび集団の最良値gbestを決定する。指定された繰り返し回数kが終了した後、集団の最良値gbestの状態量をPIDパラメータである比例ゲインK、積分時間T、微分時間Tとしたものを、最適なPIDパラメータとする。以上により、PID調節器20の最適なPIDパラメータを求めることができる。 Here, the evaluation of each agent i can be performed by the following procedure. First, the state variable s i of agent i, proportional gain K P is a PID parameter, integral time T i, set the derivative time T D, then the disturbance response to setpoint response and the disturbance d for the set value SV From the result, the evaluation index (Formula 3) is evaluated. Based on the evaluation results of all agents, the best value pbest i of each search point and the best value gbest of the group are determined. After the designated number of repetitions k has ended, the optimal PID parameter is obtained by setting the state quantity of the best value gbest of the group as a proportional gain K P , integration time T i , and differentiation time T D as PID parameters. As described above, the optimum PID parameter of the PID controller 20 can be obtained.

図1に示されるPID調整装置30の伝達関数パラメータ取得部32等の各機能は、ソフトウェアまたはハードウェアを一部または全部に用いて実現することができる。図6は、PID調整装置30の各機能をすべてソフトウェアにより実現した場合の処理の流れをフローチャートで示す。ソフトウェアまたはハードウェアを一部または全部に用いて実現した場合も同様のフローチャートにより説明することができる。   Each function of the transfer function parameter acquisition unit 32 and the like of the PID adjustment device 30 shown in FIG. 1 can be realized by using software or hardware partly or entirely. FIG. 6 is a flowchart showing the flow of processing when all the functions of the PID adjustment apparatus 30 are realized by software. A case where the software or hardware is partially or wholly implemented can be described by the same flowchart.

図6に示されるように、制御対象10を閉ループにて制御するPID調節器20のPIDパラメータを求めるPID調整プログラムは、当該PID調整御プログラムを実行するPID調整装置のコンピュータ(後述するCPU41)に、以下の各ステップを実行させる。まず、記録部(後述する記録装置48等)に記録された制御対象10に入力する操作量MVの時系列データと操作量MVに応じて制御対象10から出力される制御量PVの時系列データとに基づき、制御対象10を表現する複数の伝達関数11等のパラメータを求める(伝達関数パラメータ取得ステップ。ステップS10)。次に、伝達関数パラメータ取得ステップ(ステップS10)で求められたパラメータを有する各伝達関数11等が出力した推定制御量PVと実際の制御量PVとの適合性(適合率FIT(%))を求める(適合性取得ステップ。ステップS12)。続いて、適合性取得ステップ(ステップS12)で求められた適合性に基づき、最良の適合性を有する伝達関数を選択するまたは選択させる(最良伝達関数選択ステップ。ステップS14)。最後に、最良伝達関数選択ステップ(ステップS14)で選択されたまたは選択させた伝達関数を対象として、設定値SVに対する設定応答および外乱dに対する外乱応答の結果から予め決めておいた評価指標(式3)を満たすように、PID調節器20のPIDパラメータを探索する(最適PIDパラメータ探索ステップ。ステップS16)。 As shown in FIG. 6, the PID adjustment program for obtaining the PID parameter of the PID adjuster 20 that controls the control target 10 in a closed loop is stored in a computer (CPU 41 described later) of the PID adjustment apparatus that executes the PID adjustment control program. The following steps are executed. First, the time series data of the operation amount MV input to the control object 10 recorded in the recording unit (recording device 48 or the like described later) and the time series data of the control amount PV output from the control object 10 according to the operation amount MV. Based on the above, parameters such as a plurality of transfer functions 11 expressing the control object 10 are obtained (transfer function parameter acquisition step, step S10). Next, the suitability (fitness rate FIT (%)) between the estimated control amount PV * output by each transfer function 11 having the parameters obtained in the transfer function parameter acquisition step (step S10) and the actual control amount PV. (Adaptability acquisition step. Step S12). Subsequently, based on the suitability obtained in the suitability acquisition step (step S12), the transfer function having the best suitability is selected or selected (best transfer function selection step; step S14). Finally, with respect to the transfer function selected or selected in the best transfer function selection step (step S14), an evaluation index (formula) determined in advance from the result of the set response to the set value SV and the disturbance response to the disturbance d The PID parameter of the PID adjuster 20 is searched so as to satisfy 3) (optimum PID parameter searching step, step S16).

以上より、本発明の実施例1によれば、PID調整装置30は、伝達関数パラメータ取得部32、適合性取得部34、最良伝達関数選択部36および最適PIDパラメータ探索部38の各機能を有している。伝達関数パラメータ取得部32は、記録部(後述する記録装置48等)に記録された制御対象10に入力する操作量MVの時系列データとこの操作量MVに応じて制御対象10から出力される制御量PVの時系列データとに基づき、制御対象10を表現する複数の伝達関数のパラメータ(ゲイン定数G、時定数T、T 、T 、T 、むだ時間L、減衰係数ζ、固有角周波数ω)を求める。伝達関数は、オフラインにて予め構造を決めておいた制御対象10を表現する複数(例えば図2に示される7つの伝達関数11等)を用いることが好適であるが、これらの伝達関数11等に限定されるものではない。各伝達関数11ないし17のパラメータは、式1で示される目的関数Jの最小値を求めることにより得ることができる。適合性取得部34は、伝達関数パラメータ取得部32により求められたパラメータを有する各伝達関数が出力した推定制御量(点線で示されるPV)と実際の制御量(点線で示される方のPV)との適合性を求める。適合性は、式2で表される適合率FIT(%)により求めることができる。適合率FIT(%)の代わりに2つのデータの相関を用いることも可能である。最良伝達関数選択部36は、適合性取得部34により求められた適合性FIT(%)に基づき、最良の適合性を有する伝達関数を選択するまたは選択させる。上記選択するまたは選択させる状態は後述するディスプレイ等の表示装置44に表示することができる。選択は最良伝達関数選択部36が適合率FIT(%)の値により自動的に選択するが、手動にて任意の伝達関数11等を選択することも可能である。最適PIDパラメータ探索部38は、最良伝達関数選択部36により選択されたまたは選択させた伝達関数11等を対象として、設定値SVに対する設定値応答および外乱dに対する外乱応答の結果から予め決めておいた評価指標を満たすように、PID調節器20のPIDパラメータを探索する。評価指標の計算は、式3で表されるISEを用いて行うことができる。評価指標である式3は、例えば、PSO等の探索手法を有する探索手段により、PID調節器20のパラメータである、比例ゲインK、積分時間T、微分時間Tを求めることができる。 As described above, according to the first embodiment of the present invention, the PID adjustment device 30 has the functions of the transfer function parameter acquisition unit 32, the suitability acquisition unit 34, the best transfer function selection unit 36, and the optimum PID parameter search unit 38. is doing. The transfer function parameter acquisition unit 32 is output from the control target 10 according to the time series data of the operation amount MV input to the control target 10 recorded in the recording unit (a recording device 48 described later) and the operation amount MV. Based on the time series data of the controlled variable PV, parameters of a plurality of transfer functions (gain constant G * , time constant T * , T 1 * , T 2 * , T 3 * , dead time L * representing the controlled object 10 are represented . , Attenuation coefficient ζ * , natural angular frequency ω * ). It is preferable to use a plurality of transfer functions (for example, the seven transfer functions 11 shown in FIG. 2) representing the control object 10 whose structure has been previously determined offline. It is not limited to. The parameters of the transfer functions 11 to 17 can be obtained by obtaining the minimum value of the objective function J expressed by the equation 1. The conformity acquisition unit 34 outputs the estimated control amount (PV * indicated by the dotted line) output by each transfer function having the parameter obtained by the transfer function parameter acquisition unit 32 and the actual control amount (PV indicated by the dotted line). ). The suitability can be obtained from the fit rate FIT (%) expressed by Equation 2. It is also possible to use the correlation of two data instead of the precision FIT (%). The best transfer function selection unit 36 selects or selects a transfer function having the best fitness based on the fitness FIT (%) obtained by the fitness acquisition unit 34. The selection or selection state can be displayed on a display device 44 such as a display described later. The selection is automatically selected by the best transfer function selection unit 36 based on the value of the precision FIT (%), but any transfer function 11 or the like can be manually selected. The optimum PID parameter search unit 38 is determined in advance from the result of the set value response to the set value SV and the disturbance response to the disturbance d for the transfer function 11 or the like selected or selected by the best transfer function selection unit 36. The PID parameter of the PID controller 20 is searched so as to satisfy the evaluation index. The calculation of the evaluation index can be performed using the ISE represented by Expression 3. Equation 3, which is an evaluation index, can determine, for example, the proportional gain K P , integration time T i , and differentiation time T D that are parameters of the PID adjuster 20 by search means having a search method such as PSO.

以上のように本発明の実施例1によれば、PID調整装置30の最適PIDパラメータ探索部38が、最良伝達関数選択部36により選択されたまたは選択させた伝達関数11等を対象として、設定値SVに対する設定値応答および外乱dに対する外乱応答の結果から予め決めておいた評価指標を満たすように、PID調節器20のPIDパラメータを探索することができる。この結果、作業者が制御対象10に関する経験・ノウハウに基づき、PIDパラメータの決定を試行錯誤的に行う必要がなくなるため、PID調節器の調整方法およびPID調節装置等に関して、広く現場のエンジニアに普及させることが可能となり、加えて調整労力と時間とを著しく短縮することができる。   As described above, according to the first embodiment of the present invention, the optimum PID parameter searching unit 38 of the PID adjusting device 30 sets the transfer function 11 or the like selected or selected by the best transfer function selecting unit 36 as a target. The PID parameter of the PID adjuster 20 can be searched so as to satisfy a predetermined evaluation index from the results of the set value response to the value SV and the disturbance response to the disturbance d. As a result, it is not necessary for the operator to determine the PID parameters based on experience and know-how related to the controlled object 10, so that the method for adjusting the PID adjuster, the PID adjuster, etc. are widely used by engineers on site. In addition, the adjustment labor and time can be significantly reduced.

図7は、本発明の実施例2におけるPID調整装置を示す。図7で図1と同じ符号を付した箇所は同じ要素を示すため説明は省略する。実施例2が実施例1と異なる点は、図7に示されるように、制御量PVに白色ノイズw(上記式4の重み関数とは異なる。)が注入した状態となっている点にある。   FIG. 7 shows a PID adjustment apparatus according to Embodiment 2 of the present invention. In FIG. 7, the same reference numerals as those in FIG. The difference between the second embodiment and the first embodiment is that, as shown in FIG. 7, white noise w (different from the weighting function of Equation 4 above) is injected into the control amount PV. .

図8は、実施例2における最適PIDパラメータ探索部38の機能を説明するための図である。図8で図7と同じ符号を付した箇所は同じ要素を示すため説明は省略する。図8に示されるように、制御量PVに白色ノイズwを注入した状態で、初めに設定値SVにステップ状の信号を入力し、ステップ応答が収束した後、外乱dにステップ上の信号を入力し、出力である制御量PVから評価指標を計算する。図9は、図8で説明した設定値SV、外乱dおよび制御量PVのタイムチャートを示す。図9(A)は設定値SV、図9(B)は外乱d、図9(C)は制御量PV、横軸はいずれも時間であり、時間tdは操作量MVに外乱dが加えられた(+)時間を示す。図9(C)に示されるように、制御量PVに白色ノイズwが注入された状態がわかる。制御量PVに白色ノイズwを注入することで、PIDパラメータのハイゲイン化(最適PIDパラメータのゲインが大きくなる現象)を防ぐことができる。以上のように、実施例2では、PID調整装置30の最適PIDパラメータ探索部38が用いる設定値SVに対する設定値応答および外乱dに対する外乱応答の結果は、制御量PVに白色ノイズwを注入した状態の設定値応答および外乱応答の結果である。最適PIDパラメータ探索部38は、制御量PVに白色ノイズwを注入した状態の設定値応答および外乱応答の結果から、予め決めておいた評価指標を満たすように、PID調節器20のパラメータを探索する。図8に示される評価指標の計算では、実施例1と同様に、式3で表されるISEを用いて評価指標を計算する。図9では、実施例1と同様に、e(t)につき、制御量PVの立ち上がり時間の部分と外乱応答が加えられた部分とを示してある。評価指標である式3は、実施例1と同様に、PSO等の探索手法により、PID調節器20のパラメータである、比例ゲインK、積分時間T、微分時間Tを求めることができる。 FIG. 8 is a diagram for explaining the function of the optimum PID parameter search unit 38 in the second embodiment. In FIG. 8, the same reference numerals as those in FIG. As shown in FIG. 8, in a state where white noise w is injected into the control amount PV, a step-like signal is first input to the set value SV, and after the step response has converged, the signal on the step is applied to the disturbance d. The evaluation index is calculated from the control amount PV that is input and output. FIG. 9 is a time chart of the set value SV, the disturbance d, and the control amount PV described in FIG. 9A is the set value SV, FIG. 9B is the disturbance d, FIG. 9C is the control amount PV, the horizontal axis is time, and the time td is the disturbance d added to the manipulated variable MV. (+) Time. As shown in FIG. 9C, it can be seen that the white noise w is injected into the control amount PV. By injecting the white noise w into the control amount PV, it is possible to prevent a high gain of the PID parameter (a phenomenon in which the gain of the optimal PID parameter increases). As described above, in the second embodiment, the white noise w is injected into the control amount PV as the result of the set value response to the set value SV and the disturbance response to the disturbance d used by the optimum PID parameter search unit 38 of the PID adjusting device 30. It is a result of a set value response of a state and a disturbance response. The optimum PID parameter search unit 38 searches for the parameters of the PID adjuster 20 so as to satisfy a predetermined evaluation index from the result of the set value response and the disturbance response in the state where the white noise w is injected into the control amount PV. To do. In the calculation of the evaluation index shown in FIG. 8, the evaluation index is calculated using the ISE represented by Expression 3, as in the first embodiment. In FIG. 9, as in the first embodiment, the rise time portion of the control amount PV and the portion to which the disturbance response is added are shown for e (t). As in the first embodiment, Expression 3, which is an evaluation index, can determine the proportional gain K P , integration time T i , and differentiation time T D that are parameters of the PID controller 20 by a search method such as PSO. .

実施例2におけるPID調整装置30の他の機能については実施例1と同様であるため、説明は省略する。PID調整装置30の各機能は、実施例1と同様に、ソフトウェアまたはハードウェアを一部または全部に用いて実現することができる。PID調整プログラムの機能についても、最適PIDパラメータ探索ステップ(ステップS16)で用いる設定値応答および外乱応答の結果は、制御量PVに白色ノイズwを注入した状態の設定値応答および外乱応答の結果である点を除き実施例1と同様であるため、説明は省略する。   Since the other functions of the PID adjusting device 30 in the second embodiment are the same as those in the first embodiment, description thereof is omitted. Each function of the PID adjustment device 30 can be realized by using a part or all of software or hardware as in the first embodiment. Regarding the function of the PID adjustment program, the result of the set value response and the disturbance response used in the optimum PID parameter search step (step S16) is the result of the set value response and the disturbance response in the state where the white noise w is injected into the control amount PV. Since it is the same as that of Example 1 except for a certain point, description is abbreviate | omitted.

以上より、本発明の実施例2によれば、制御量PVに白色ノイズwを注入した状態で、初めに設定値SVにステップ状の信号を入力し、ステップ応答が収束した後、外乱dにステップ上の信号を入力し、出力である制御量PVから評価指標を計算する。最適PIDパラメータ探索部38は、制御量PVに白色ノイズwを注入した状態の設定値応答および外乱応答の結果から、予め決めておいた評価指標を満たすように、PID調節器20のパラメータを探索する。以上により、実施例1の効果に加えて、制御量PVに白色ノイズwを注入することで、PIDパラメータのハイゲイン化(最適PIDパラメータのゲインが大きくなる現象)を防ぐことができるという効果を得ることができる。   As described above, according to the second embodiment of the present invention, in the state where the white noise w is injected into the control amount PV, a step-like signal is first input to the set value SV, and after the step response has converged, the disturbance d A signal on the step is input, and an evaluation index is calculated from the control amount PV that is an output. The optimum PID parameter search unit 38 searches for the parameters of the PID adjuster 20 so as to satisfy a predetermined evaluation index from the result of the set value response and the disturbance response in the state where the white noise w is injected into the control amount PV. To do. As described above, in addition to the effect of the first embodiment, by injecting the white noise w into the control amount PV, it is possible to prevent the high gain of the PID parameter (a phenomenon in which the gain of the optimal PID parameter is increased). be able to.

上述したように、本発明のPID調整装置30は全体をハードウェアとして実現することもでき、あるいは一部または全体をソフトウェアとして実現することも可能である。図10は、本発明のPID調整プログラムを実行するPID調整装置30内のコンピュータの内部回路40を示すブロック図である。図10に示されるように、CPU41、ROM42、RAM43、画像制御部46、コントローラ47、入力制御部50および外部I/F部52はバス53に接続されている。図10において、上述の本発明のPID調整プログラムは、ROM42、ディスク48またはDVD若しくはCD−ROM49等の記録媒体(脱着可能な記録媒体を含む。)に記録されている。ディスク48(記録部)には、伝達関数パラメータ取得部32が用いる操作量MVの時系列データファイル61および制御量の時系列データファイル62等を記録しておくことができる。これらの時系列データファイル61等はDVD若しくはCD−ROM49等の記録媒体49に記録されていてもよい。伝達関数パラメータ取得部32が求めた複数の伝達関数のパラメータ、適合性取得部34が求めた適合率FIT(%)、最良伝達関数選択部36が選択したまたは選択させた伝達関数等は一時的にRAM43に記録しておくことができる。PID調整プログラムは、ROM42からバス53を介し、あるいはディスク48またはDVD若しくはCD−ROM49等の記録媒体からコントローラ47を経由してバス53を介しRAM43へロードされる。画像制御部46は、所望の画像データ、例えば図3に示される最良伝達関数選択部36が適合性FIT(%)に基づき、最良の適合性を有する伝達関数を選択するまたは選択させる状態を示す画像データをVRAM45へ送出する。ディスプレイ等の表示装置44はVRAM45から送出された上記画像データ等を表示する。VRAM45は表示装置44の一画面分のデータ容量に相当する容量を有している画像メモリである。入力操作部51はコンピュータに入力を行うためのマウス、テンキー等の入力装置であり、入力制御部50は入力操作部51と接続され入力制御等を行う。外部I/F部52はコンピュータの外部と接続する際のインタフェース機能を有している。   As described above, the PID adjustment apparatus 30 of the present invention can be realized entirely as hardware, or part or all can be realized as software. FIG. 10 is a block diagram showing an internal circuit 40 of a computer in the PID adjustment apparatus 30 that executes the PID adjustment program of the present invention. As shown in FIG. 10, the CPU 41, ROM 42, RAM 43, image control unit 46, controller 47, input control unit 50, and external I / F unit 52 are connected to a bus 53. In FIG. 10, the above-described PID adjustment program of the present invention is recorded on a recording medium (including a removable recording medium) such as a ROM 42, a disk 48, a DVD or a CD-ROM 49. On the disk 48 (recording unit), the time series data file 61 of the operation amount MV used by the transfer function parameter acquisition unit 32, the time series data file 62 of the control amount, and the like can be recorded. These time series data files 61 and the like may be recorded on a recording medium 49 such as a DVD or a CD-ROM 49. The parameters of a plurality of transfer functions obtained by the transfer function parameter acquisition unit 32, the relevance ratio FIT (%) obtained by the suitability acquisition unit 34, the transfer function selected or selected by the best transfer function selection unit 36, etc. are temporary. Can be recorded in the RAM 43. The PID adjustment program is loaded into the RAM 43 from the ROM 42 via the bus 53 or from a recording medium such as the disk 48 or DVD or CD-ROM 49 via the controller 47 via the bus 53. The image control unit 46 indicates a state in which desired image data, for example, the best transfer function selection unit 36 shown in FIG. 3 selects or selects a transfer function having the best fitness based on the fitness FIT (%). The image data is sent to the VRAM 45. A display device 44 such as a display displays the image data sent from the VRAM 45. The VRAM 45 is an image memory having a capacity corresponding to the data capacity of one screen of the display device 44. The input operation unit 51 is an input device such as a mouse or a numeric keypad for inputting to the computer. The input control unit 50 is connected to the input operation unit 51 and performs input control. The external I / F unit 52 has an interface function when connecting to the outside of the computer.

上述のようにCPU41が本発明のPID調整プログラムを実行することにより、本発明の目的を達成することができる。PID調整プログラムは上述のようにDVD若しくはCD−ROM49等の記録媒体の形態でコンピュータCPU41に供給することができ、PID調整プログラムを記録したDVD若しくはCD−ROM49等の記録媒体も同様に本発明を構成することになる。PID調整プログラムを記録した記録媒体としては上述された記録媒体の他に、例えばメモリ・カード、メモリ・スティック、光ディスク、FD等を用いることができる。   As described above, the CPU 41 executes the PID adjustment program of the present invention, whereby the object of the present invention can be achieved. As described above, the PID adjustment program can be supplied to the computer CPU 41 in the form of a recording medium such as a DVD or a CD-ROM 49, and a recording medium such as a DVD or a CD-ROM 49 on which the PID adjustment program is recorded similarly applies the present invention. Will be composed. As a recording medium on which the PID adjustment program is recorded, for example, a memory card, a memory stick, an optical disc, an FD, or the like can be used in addition to the recording medium described above.

本発明の活用例として、石油、化学、鉄鋼、紙パルプ、水環境等のいわゆるプロセス産業におけるPID制御方式に適用することができる。   As an application example of the present invention, it can be applied to a PID control system in a so-called process industry such as petroleum, chemical, steel, paper pulp, and water environment.

10 制御対象、 11、12、13、14、15、16、17 伝達関数、 20 PID調節器、 30 PID調整装置、 32 伝達関数パラメータ取得部、 34適合性取得部、 36 最良伝達関数選択部、 38 最適PIDパラメータ探索部、 40 内部回路、 41 CPU、 42 ROM、 43 RAM、 44 表示装置、 45 VRAM、 46 画像制御部、 47 コントローラ、 49 DVD、CD−ROM、 50 入力制御部、 51 入力操作部、 52 外部I/F部、 53 バス、 61 操作量MVの時系列データファイル、 62 制御量の時系列データファイル。   10 control object 11, 12, 13, 14, 15, 16, 17 transfer function, 20 PID adjuster, 30 PID adjustment device, 32 transfer function parameter acquisition unit, 34 compatibility acquisition unit, 36 best transfer function selection unit, 38 Optimal PID parameter search unit, 40 internal circuit, 41 CPU, 42 ROM, 43 RAM, 44 display device, 45 VRAM, 46 image control unit, 47 controller, 49 DVD, CD-ROM, 50 input control unit, 51 input operation Part, 52 external I / F part, 53 bus, 61 time series data file of manipulated variable MV, 62 time series data file of controlled variable.

特開2004−38428号公報JP 2004-38428 A

Claims (6)

制御対象を閉ループにて制御するPID調節器のPIDパラメータを求めるPID調整装置であって、
記録部に記録された制御対象に入力する操作量の時系列データと該操作量に応じて制御対象から出力される制御量の時系列データとに基づき、制御対象を表現する複数の伝達関数のパラメータを求める伝達関数パラメータ取得手段と、
前記伝達関数パラメータ取得手段により求められたパラメータを有する各伝達関数が出力した推定制御量と実際の制御量との適合性を求める適合性取得手段と、
前記適合性取得手段により求められた適合性に基づき、最良の適合性を有する伝達関数を選択する又は選択させる最良伝達関数選択手段と、
前記最良伝達関数選択手段により選択された又は選択させた伝達関数を対象として、設定値応答及び外乱応答の結果から予め決めておいた評価指標を満たすべくPID調節器のPIDパラメータを探索する最適PIDパラメータ探索手段とを備えたことを特徴とするPID調整装置。
A PID adjusting device for obtaining a PID parameter of a PID adjuster that controls a control target in a closed loop,
Based on the time-series data of the operation amount input to the control object recorded in the recording unit and the time-series data of the control amount output from the control object according to the operation amount, a plurality of transfer functions representing the control object A transfer function parameter obtaining means for obtaining a parameter;
Suitability obtaining means for obtaining suitability between the estimated control amount output by each transfer function having the parameter obtained by the transfer function parameter obtaining means and the actual control amount;
A best transfer function selecting means for selecting or selecting a transfer function having the best suitability based on the suitability obtained by the suitability obtaining means;
Optimal PID for searching for a PID parameter of a PID controller so as to satisfy a predetermined evaluation index from results of set value response and disturbance response for the transfer function selected or selected by the best transfer function selection means A PID adjusting device comprising a parameter search means.
請求項1記載のPID調整装置において、前記最適PIDパラメータ探索手段が用いる設定値応答及び外乱応答の結果は、制御量に白色ノイズを注入した状態の設定値応答及び外乱応答の結果であることを特徴とするPID調整装置。   The PID adjustment apparatus according to claim 1, wherein the result of the set value response and the disturbance response used by the optimum PID parameter search means is a result of the set value response and the disturbance response in a state where white noise is injected into the control amount. A characteristic PID adjusting device. 請求項1又は2記載のPID調整装置において、前記最適PIDパラメータ探索手段が用いる評価指標は、制御量とステップ状の設定値との差の2乗を時間で積分した指標であることを特徴とするPID調整装置。   3. The PID adjustment apparatus according to claim 1, wherein the evaluation index used by the optimum PID parameter search means is an index obtained by integrating the square of the difference between the control amount and the step-like set value with time. PID adjustment device. 請求項1乃至3のいずれかに記載のPID調整装置において、前記適合性取得手段が用いる適合性は、推定制御量、実際の制御量及び実際の制御量の平均値に基づく所定の適合率であることを特徴とするPID調整装置。   4. The PID adjustment apparatus according to claim 1, wherein the compatibility used by the compatibility acquisition unit is a predetermined adaptation rate based on an estimated control amount, an actual control amount, and an average value of the actual control amount. There is a PID adjusting device. 制御対象を閉ループにて制御するPID調節器のPIDパラメータを求めるPID調整プログラムであって、PID調整御プログラムを実行するPID調整装置のコンピュータに、
記録部に記録された制御対象に入力する操作量の時系列データと該操作量に応じて制御対象から出力される制御量の時系列データとに基づき、制御対象を表現する複数の伝達関数のパラメータを求める伝達関数パラメータ取得ステップ、
前記伝達関数パラメータ取得ステップで求められたパラメータを有する各伝達関数が出力した推定制御量と実際の制御量との適合性を求める適合性取得ステップ、
前記適合性取得ステップで求められた適合性に基づき、最良の適合性を有する伝達関数を選択する又は選択させる最良伝達関数選択ステップ、
前記最良伝達関数選択ステップで選択された又は選択させた伝達関数を対象として、設定値応答及び外乱応答の結果から予め決めておいた評価指標を満たすべくPID調節器のPIDパラメータを探索する最適PIDパラメータ探索ステップを実行させるためのPID調整プログラム。
A PID adjustment program for obtaining a PID parameter of a PID adjuster that controls a controlled object in a closed loop, the computer of the PID adjustment device that executes the PID adjustment control program,
Based on the time-series data of the operation amount input to the control object recorded in the recording unit and the time-series data of the control amount output from the control object according to the operation amount, a plurality of transfer functions representing the control object A transfer function parameter obtaining step for obtaining a parameter;
A suitability obtaining step for obtaining suitability between the estimated control amount output by each transfer function having the parameters obtained in the transfer function parameter obtaining step and the actual control amount;
A best transfer function selecting step for selecting or selecting a transfer function having the best suitability based on the suitability obtained in the suitability obtaining step;
Optimum PID for searching the PID parameter of the PID controller to satisfy the evaluation index determined in advance from the set value response and disturbance response results for the transfer function selected or selected in the best transfer function selection step A PID adjustment program for executing a parameter search step.
請求項5記載のPID調整プログラムにおいて、前記最適PIDパラメータ探索ステップで用いる設定値応答及び外乱応答の結果は、制御量に白色ノイズを注入した状態の設定値応答及び外乱応答の結果であることを特徴とするPID調整プログラム。   6. The PID adjustment program according to claim 5, wherein the result of the set value response and the disturbance response used in the optimum PID parameter search step is a result of the set value response and the disturbance response in a state where white noise is injected into the control amount. A featured PID adjustment program.
JP2009116015A 2009-05-12 2009-05-12 PID adjustment device and PID adjustment program Active JP5549112B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009116015A JP5549112B2 (en) 2009-05-12 2009-05-12 PID adjustment device and PID adjustment program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009116015A JP5549112B2 (en) 2009-05-12 2009-05-12 PID adjustment device and PID adjustment program

Publications (2)

Publication Number Publication Date
JP2010266967A true JP2010266967A (en) 2010-11-25
JP5549112B2 JP5549112B2 (en) 2014-07-16

Family

ID=43363913

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009116015A Active JP5549112B2 (en) 2009-05-12 2009-05-12 PID adjustment device and PID adjustment program

Country Status (1)

Country Link
JP (1) JP5549112B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3173881A1 (en) 2015-11-30 2017-05-31 Omron Corporation Control parameter tuning device, control parameter tuning method, control parameter tuning program
KR20210013971A (en) * 2019-07-29 2021-02-08 김푸름 Automatic learning tuning system of motor controller using PSO
CN113204876A (en) * 2021-04-30 2021-08-03 广东电网有限责任公司电力科学研究院 Noise gain calculation method, device, equipment and medium of PD controller
WO2023026552A1 (en) * 2021-08-26 2023-03-02 オムロン株式会社 Temperature control unit and temperature control method
CN116281422A (en) * 2023-05-18 2023-06-23 青岛宏大纺织机械有限责任公司 Control method and system of tension controller of bobbin winder

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2714567C1 (en) * 2019-05-06 2020-02-18 федеральное государственное бюджетное образовательное учреждение высшего образования "Ивановский государственный энергетический университет имени В.И. Ленина" (ИГЭУ) Automatic adjustment method of regulator

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0723404A (en) * 1993-06-30 1995-01-24 Nippon Dempa Kogyo Co Ltd Video camera equipment
JPH0816205A (en) * 1994-06-30 1996-01-19 Fujitsu Ltd Servo control system and its tuning method
JP2001273005A (en) * 2000-03-27 2001-10-05 Seiko Epson Corp Servo control system, electronic instrument and method for adjusting gain
JP2004038428A (en) * 2002-07-02 2004-02-05 Yamatake Corp Method for generating model to be controlled, method for adjusting control parameter, program for generating the model, and program for adjusting the parameter
JP2007287153A (en) * 2006-04-13 2007-11-01 Fisher Rosemount Syst Inc Robust process model identification method and system in model based control technique

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0723404A (en) * 1993-06-30 1995-01-24 Nippon Dempa Kogyo Co Ltd Video camera equipment
JPH0816205A (en) * 1994-06-30 1996-01-19 Fujitsu Ltd Servo control system and its tuning method
JP2001273005A (en) * 2000-03-27 2001-10-05 Seiko Epson Corp Servo control system, electronic instrument and method for adjusting gain
JP2004038428A (en) * 2002-07-02 2004-02-05 Yamatake Corp Method for generating model to be controlled, method for adjusting control parameter, program for generating the model, and program for adjusting the parameter
JP2007287153A (en) * 2006-04-13 2007-11-01 Fisher Rosemount Syst Inc Robust process model identification method and system in model based control technique

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3173881A1 (en) 2015-11-30 2017-05-31 Omron Corporation Control parameter tuning device, control parameter tuning method, control parameter tuning program
KR20210013971A (en) * 2019-07-29 2021-02-08 김푸름 Automatic learning tuning system of motor controller using PSO
CN113204876A (en) * 2021-04-30 2021-08-03 广东电网有限责任公司电力科学研究院 Noise gain calculation method, device, equipment and medium of PD controller
WO2023026552A1 (en) * 2021-08-26 2023-03-02 オムロン株式会社 Temperature control unit and temperature control method
CN116281422A (en) * 2023-05-18 2023-06-23 青岛宏大纺织机械有限责任公司 Control method and system of tension controller of bobbin winder
CN116281422B (en) * 2023-05-18 2023-08-18 青岛宏大纺织机械有限责任公司 Control method and system of tension controller of bobbin winder

Also Published As

Publication number Publication date
JP5549112B2 (en) 2014-07-16

Similar Documents

Publication Publication Date Title
JP5549112B2 (en) PID adjustment device and PID adjustment program
JP4223894B2 (en) PID parameter adjustment device
JP2009116515A (en) Pid parameter adjustment support device
JP5768834B2 (en) Plant model management apparatus and method
TWI743986B (en) Motor control device and motor control method
JP2012208789A (en) Reinforcement learning device, control device, and reinforcement learning method
KR20050010985A (en) Control object model generation method, program thereof, control parameter adjustment method, and program thereof
EP3899797A1 (en) Multi-agent reinforcement learning with matchmaking policies
CN101390024A (en) Operation control method, operation control device and operation control system
Norgaard et al. NNSYSID and NNCTRL tools for system identification and control with neural networks
Konakoglu et al. 2D coordinate transformation using artificial neural networks
Vrabie Online adaptive optimal control for continuous-time systems
JP4698578B2 (en) Methods and articles for detecting, verifying, and repairing collinearity
JP6652398B2 (en) Parameter selection method, parameter selection program and parameter selection device
CN110210548A (en) A kind of picture dynamic self-adapting compression method based on intensified learning
CN117540203A (en) Multi-directional course learning training method and device for cooperative navigation of clustered robots
Liu et al. Efficient reinforcement learning control for continuum robots based on inexplicit prior knowledge
CN111095133A (en) Method and apparatus for providing an adaptive self-learning control program for deployment on a target field device
WO2021186500A1 (en) Learning device, learning method, and recording medium
JP5582487B2 (en) Process state prediction method
Fakhari et al. Quantum inspired reinforcement learning in changing environment
CN112583322A (en) Active-disturbance-rejection control method for linear motor and related device
JP6850313B2 (en) Simulation equipment, simulation method and simulation program
JP6848710B2 (en) Plant control and adjustment equipment and methods
CN115972216A (en) Parallel robot positive motion solving method, control method, equipment and storage medium

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20110422

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120313

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130411

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130514

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130712

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20131112

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140108

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20140422

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140505

R150 Certificate of patent or registration of utility model

Ref document number: 5549112

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250