JP3571952B2 - 2-DOF control device - Google Patents
2-DOF control device Download PDFInfo
- Publication number
- JP3571952B2 JP3571952B2 JP01602299A JP1602299A JP3571952B2 JP 3571952 B2 JP3571952 B2 JP 3571952B2 JP 01602299 A JP01602299 A JP 01602299A JP 1602299 A JP1602299 A JP 1602299A JP 3571952 B2 JP3571952 B2 JP 3571952B2
- Authority
- JP
- Japan
- Prior art keywords
- degree
- freedom
- control
- parameter
- control device
- 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.)
- Expired - Lifetime
Links
Images
Landscapes
- Feedback Control In General (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、石油,化学,鉄鋼,電力,製紙などの各種プラント計装分野やサーボ制御分野などに利用される2自由度制御装置に係り、特にプラント各箇所の流量,圧力,液位,組成,位置,速度,回転数などのプロセス変量をフィードバック制御する2自由度制御装置に関する。
【0002】
【従来の技術】
外乱入力時および目標値変更時の何れでも制御量を良好にフィードバック制御するためには、従来から2自由度PID制御装置が用いられている。
【0003】
ここで、2自由度PID制御とは、外乱に対する抑制特性と目標値変更に対する追従特性との双方を同時に最適化するPID制御である。
【0004】
図20は従来の一般的な2自由度PIDまたは2自由度I−PD制御装置のブロック構成図である。
【0005】
この制御装置は2自由度PID(以下、2自由度I−PDを含む)制御器100および制御対象101によって構成され、この2自由PID度制御器100は、比例系、積分系および微分系からなっている。
【0006】
前記比例系は、目標値Rに係数要素111の比例項用2自由度化パラメータbを乗算し、この乗算値から減算要素112にて制御量Yを減算し、加算要素113に導入する。一方、積分系は、比較要素114において目標値Rと制御量Yとの制御偏差Eを求め、得られた制御偏差Eについて積分時間TIをもつ積分要素 115で積分し、加算要素113に導入する。さらに、微分系は、目標値Rに係数要素116の微分項用2自由度化パラメータcを乗算して減算要素117に導き、ここで乗算値から制御量Yを減算し、得られた減算値について微分時間TD をもつ微分要素118で微分した後、前記加算要素113に導入する。
【0007】
このようにして得られた減算要素112出力、積分要素115出力および微分要素118出力は加算要素113で加算された後、この加算値に比例ゲイン係数要素119の比例ゲインKpを乗算することによって操作量MVを求めた後、外 乱量とともに制御対象101に印加する。
【0008】
また、従来、図21に示すように、目標値変更を検出し、その影響を打ち消すために予測先行的に操作量変化を判断し、先回りして操作を実行するフィードフォワード制御方式を採用した2自由度PID制御器が用いられている。
【0009】
この2自由度PID制御器100は、目標値Rを目標値フィードフォワード要素121に入力し、ここで比例項用2自由度化パラメータbおよび微分項用2自由度化パラメータcを用いて、(b+cTDs)なる演算を実行し、得られた出 力を発信する。また、比較要素114において目標値Rと制御量Yとの制御偏差Eを求めた後、この制御偏差Eについて積分時間TIをもつ積分要素115で積 分し、減算要素122に導入する。一方、微分系では,制御量Yに比例・微分要素123の比例・微分演算式(1+TDs)を用いて微分し、減算要素122に 導入する。ここで、積分要素115の出力から比例・微分要素123の出力を減算し、この減算結果に比例ゲイン係数要素119の比例ゲインKPを乗算し、こ の乗算値と目標値フィードフォワード要素121の出力とを加算要素124で加算し、制御対象101(図示せず)に印加する操作量MVを求める構成である。
【0010】
さらに、従来、図13に示すように、目標値フィルタ要素131を付加した2自由度PID制御装置が用いられている。
【0011】
この2自由度PID制御器100は、目標値Rを目標値フィルタ要素131の2自由度化パラメータb,cを用いて、下記演算式により演算を実行し、最適な目標値追従特性となる演算目標値を算出する。
【0012】
(1+bTIs+cTITDs2)/(1+TIs) ……(1)
以上のようにして求めた演算目標値と制御量Yとを比較要素114により制御偏差Eを求めた後、この制御偏差Eに比例・積分要素132にて{1+1/(TIs)}を演算する。また、制御量Yに微分要素133にて微分演算を実行する 。そして、減算要素134において比例・積分要素132の出力から微分要素133の出力を減算し、この得られた減算出力に比例ゲイン係数要素119の比例ゲインKpを乗算することによって操作量MVを求め、外乱量とともに制御対象 101(図示せず)に印加する。
【0013】
【発明が解決しようとする課題】
ところで、以上のような各種の2自由度PID制御器では、何れも外乱に対する制御応答が最適となるようなPIDパラメータ(比例ゲインKP,積分時間TI,微分時間TD)を調整した後、目標値変更に対する追従特性を最適とするため の2自由度化パラメータb,cを試行錯誤的に調整する必要がある。また、予め制御対象モデルを構築し、オフラインにてある評価関数の下にシミュレーションを実行し、最適値探索により伝達関数の各パラメータ,つまりPIDパラメータおよび2自由度化パラメータb,cを設計することが行われている。
【0014】
その結果、外乱および目標値変更に対して最適な制御応答となるようなPIDパラメータおよび2自由度化パラメータb,cを得るために、熟練者による煩雑な調整作業が必要であり、また調整作業に多くの時間を必要とする問題がある。
【0015】
本発明は上記事情にかんがみてなされたもので、目標値変更に対する最適な制御応答を得るための2自由度化パラメータを迅速、かつ、熟練を要せずに容易に発生する2自由度制御装置を提供することにある。
【0016】
【課題を解決するための手段】
上記課題を解決するために、本発明は、制御パラメータおよび2自由度化パラメータを用いて、制御対象に対して外乱抑制性および目標値追従性が共に良好に制御する2自由度制御装置において、外乱入力制御応答に関する立上り時間σD と、目標値変更制御応答に関する立上り時間σRと、前記制御パラメータである 積分時間TIあるいは当該積分時間TIおよび微分時間TDとを用いて、前記目標 値変更制御応答を規定する前記2自由度化パラメータを発生するパラメータ発生機構を設けた構成である。
【0017】
なお、以上の入力項目の他に、外乱入力制御応答規範モデルの2次項係数α2Dおよび目標値変更制御応答規範モデルの2次項係数α2Rを用い、所定の演算式により得られる値に比例する値を、前記目標値変更制御応答を規定する2自由度化パラメータとしてもよい。
【0018】
この発明は、以上のような手段を講じたことにより、予め定める所定の入力項目、つまり外乱入力制御応答に関する立上り時間σD、目標値変更制御応答に関 する立上り時間σR、前記制御パラメータである積分時間TIあるいは当該積分時間TIおよび微分時間TDを用い、さらに外乱入力制御応答規範モデルの2次項係数をα2Dおよび目標値変更制御応答規範モデルの2次項係数をα2Rを加え、実制御系の伝達関数と規範モデルの伝達関数とから導かれる所定の演算を実行することにより、目標値変更制御応答を規定する2自由度化パラメータを発生でき、従来のような試行錯誤的に調整する作業を省略できる。
【0019】
また、別の発明は、制御対象の伝達特性タイプ、むだ時間および限界周期、または制御対象の伝達特性タイプおよび基準化むだ時間を用いて、目標値変更制御応答を規定する前記2自由度化パラメータを発生するパラメータ発生機構を設けたものである。
【0020】
この発明は、以上のような手段を講じたことにより、予め制御対象の伝達特性タイプ、むだ時間および限界周期、または制御対象の伝達特性タイプおよび基準化むだ時間をもとにグラフ関数化すれば、迅速、かつ、容易に2自由度化パラメータを発生でき、上記発明と同様の作用を奏するものである。
【0021】
さらに、別の発明は、外乱抑制用I−PD制御系における第1の規範モデルの属性値と目標値追従性を規定する第2の規範モデルの属性値とが入力され、これら両属性値の関数として係数値を発生する関数要素と、この関数要素の係数値、第1の規範モデルの立上り時間および積分時間を用いて目標値追従性を良好にする2自由度化パラメータを発生する演算要素とを有するパラメータ発生機構を設けたものである。
【0022】
この発明は、外乱抑制用I−PD制御系における第1の規範モデルの属性値と目標値追従性を規定する第2の規範モデルの属性値との関数として係数値を発生し、この係数値、第1の規範モデルの立上り時間および積分時間を用いて、目標値追従性を良好にする2自由度化パラメータを発生するので、試行錯誤的な調整作業を行わずに自動的に2自由度化パラメータを発生できる。
【0023】
【発明の実施の形態】
以下、本発明の実施の形態について図面を参照して説明する。
(第1の実施の形態)
本発明に係わる2自由度制御装置の一実施の形態を説明するに先立ち、本装置を実現するための基本的な考え方について図1を参照して説明する。
(1) 先ず、目標値Rおよび制御量Yの入力のもとに、従来と同様に比例ゲインKP、積分時間TI、微分時間TDの他、比例項用2自由度化パラメータb、微分項用2自由度化パラメータcを用いて演算し、操作量MVを求める2自由度PID制御器1を設けた構成とする。なお、2自由度PID制御器1は例えば図11の2自由度PID制御器100に相当する。
(2) 次に、かかる2自由度PID制御器1において、予め外乱に対する制御応答が最適となるようなPIDパラメータである比例ゲインKP、積分時間TI 、微分時間TDを調整し、2自由度PID制御器1に設定する。
(3) さらに、新たにパラメータ発生機構2を設け、ある入力項目Aのもとに自動的に2自由度化パラメータb,cを算出する。
(4) このようにしてパラメータ発生機構2により算出される2自由度化パラメータb,cを発生し、2自由度PID制御器1に設定するものとする。
【0024】
以上の考え方を具体的に実現するためには、いかなる入力項目Aを用いるか、さらにパラメータ発生機構2は入力項目Aのもとにいかなる演算を行ってパラメータb,cを求めるかが重要になってくる。
【0025】
そこで、以下、入力項目Aおよび演算式を導き出す手順について説明する。
【0026】
従来の一般的な2自由度制御系(図20参照)に関して伝達関数をもって表わせば、図2(a)に示すような構成となる。
【0027】
同図において1は2自由度PID制御器、3は制御対象、4は加算要素、Rは目標値、Yは制御量、Dは外乱である。
【0028】
すなわち、2自由度PID制御器1は、目標値Rと制御量Yが入力され、外乱に対する抑制特性のための制御パラメータである比例ゲインKP、積分時間TI、微分時間TD、比例項用2自由度化パラメータbおよび微分項用2自由度化パラ メータcを用い、以下のような演算式によって操作量MVを求めることができる。
【0029】
KP{bR−Y+(1/TIs)(R−Y)+TDs(cR−Y)}……(2)
以上の演算式によって求められた操作量MVは外乱Dとともに制御対象3に印加する。
【0030】
ところで、図2(a)に示す2自由度PID制御系は等価ブロック変換すると、仮想的に図2(b)のような制御系で表わすことができる。5はフィードフォワード制御要素、6は加減算要素、7はフィードバック制御要素である。
【0031】
この制御系においてフィードフォワード制御要素5は、入力される目標値Rに対し、{KP(b+1/TIs+cTDs)}なる伝達関数を演算し、出力Wを発信する。また、フィードバック制御要素7は、制御対象3の出力である制御量Yに対し、{KP(1+1/TIs+TD)}なる伝達関数を演算し、出力Vを発信す る。そして、加減算要素6においてフィードフォワード制御要素5の出力Wからフィードバック制御要素7の出力Vを減算し、得られた減算出力を制御対象3に印加し、制御量Yを発生させる。このとき、外乱Dは加減算要素6に加算的に加わり、制御対象3に入力される。
【0032】
また、制御対象3の伝達関数をGP(s)とすると、フィードバック制御系に おける外乱Dから制御量Yへの伝達関数WD(s)は次式により表わせる。
【0033】
他方、外乱Dによる制御量Yの望ましい制御応答を規定する規範モデルは既に知られているように次式とする(PID制御、p14(2.11)式またはp31(2.35)式、著者:須田信英、編者:システム制御情報学会、発行所:株式会社 朝倉書店参照)。
【0034】
MD(s)={s/(KP/TI)}・[1/(1+σDs+α2DσD2s2+α3DσD3s3)] ……(4)
ここで、外乱入力に対する制御応答に関する規範モデルの係数α2D,α3Dは 速応性とロバスト性との兼ね合いから、例えば表1に示す値とする。この表1は前述する文献PID制御、p14〜p15参照。
【0035】
【表1】
【0036】
ところで、以上のような(3)式の実制御系と(4)式の規範モデルとは等しい関係になければならない。よって、WD(s)=MD(s)について解けば、外乱入力に対する制御応答に関する立上り時間σDおよびPIDパラメータKP、TI、TD が決定できる。そこで、決定されたパラメータを設定したPIDフィードバック制御系に対し、制御対象の入力側にステップ状外乱を与えたとすると、図7(a)に示すような制御応答が得られる。
【0037】
次に、図2(b)に示す制御系に対し、決定された制御パラメータKP,TI,TDを設定し、このときの目標値Rから制御量Yへの伝達関数WR(s)をもって表わすと、次式のようになる。
【0038】
WR(s)=[KP{b+1/(TIs)+cTDs}]・MD(s) ……(5)ここで、前段の[ ]内はフィードフォワード制御要素5の伝達関数、後段のMD(s)は前記(4)式の規範モデルの伝達関数である。
【0039】
従って、この(5)式は下記する(6)式で表わせる。
【0040】
WR(s)=(1+bTIs+cTITDs2)/(1+σDs+α2DσD2s2+α3DσD3s3) ……(6)
他方、目標値Rによる制御量Yの望ましい制御応答を規定する規範モデルは次式で表わすものとする(前述文献 PID制御、p120(6.2)式参照)。
【0041】
MR(s)=1/(σRs+α2RσR2s2+α3RσR3s3) …(7)
ここで、目標値変化に対する制御応答に関する規範モデルの係数α2R、α3Rは速応性とロバスト性との兼ね合いから、例えば前記表1に示す値とする。
【0042】
そこで、WR(s)=MR(s)について解くと、目標値変化に対する制御応答に関する立上り時間σRは下式の3次方程式の最小実根から求められる。
【0043】
(1−2α2R+α3R)σR3−(1−α2R)σDσR2+α2DσD2σR−α3DσD3=0 ……(8)
ここで、(6)式のWR(s)=B/Cとし、また(7)式のMR(s)=1/Zと置けば、WR(s)=MR(s)の関係から、
B/C=1/Z → Z=C/B ……(9)
が導き出せる。
【0044】
そこで、C/Bについて演算し、1次項および2次項まで解を導き出し、前記Zの1次項および2次項とそれぞれ等しい関係に置いて計算すれば、比例項用2自由度化パラメータbは、
b=(σD−σR)/TI ……(10)
なる演算式から決定できる。また、微分項用2自由度化パラメータcは、
なる演算式から決定できる。
【0045】
なお、(11)式においてα2D=α2R=0.5の場合には、
c=(σD−σR)2/(2TITD) ……(12)となり、2自由度化パラメータb,cは(σD−σR)の関数であるとも言える。
【0046】
従って、以上の説明から明らかなように、2自由度化パラメータb,cを求めるための入力項目Aおよび演算式を得ることができたので、以下、その具体例について説明する。
【0047】
図3は本発明に係る2自由度PID制御装置の一実施の形態を示す構成図である。
【0048】
同図(a)に示す制御装置は、2自由度PID制御器1にパラメータ発生機構2が接続され、このパラメータ発生機構2には、入力項目Aとして、積分時間TI、微分時間TD、外乱入力制御応答を規定する立上り時間σD、目標値変更制御 応答を規定する立上り時間σRが入力され、後記する図4に示す演算手段により 2自由度化パラメータb,cを算出し、2自由度PID制御器1に設定する構成である。
【0049】
同図(b)に示す制御装置は、同図(a)に示す入力項目の他に、さらに
外乱入力制御応答規範モデルの2次項係数α2D、目標値変更制御応答規範モデ ルの2次項係数α2Rも追加入力され、後記する図5に示す演算手段により2自 由度化パラメータb,cを算出し、2自由度PID制御器1に設定する構成である。
【0050】
すなわち、図4(a)に示すパラメータ発生機構2は比較要素11および除算要素12からなり、比較要素11は入力される外乱入力制御応答を規定する立上り時間σDから同じく入力される目標値変更制御応答を規定する立上り時間σRを減算し、その減算出力を除算要素12に導入する。この除算要素12は、減算結果(σD−σR)を積分時間TIで除算し、前記(10)式の比例項用2自由度化 パラメータbを算出し、2自由度PID制御器1に設定する。
【0051】
なお、外乱入力制御応答を規定する立上り時間σDはPIDパラメータを部分 的モデルマッチング法で設計する際に自動的に決定でき、或いは他の方法を用いてPIDパラメータを決定する場合にはステップ外乱が入った時に制御量Yがピーク値に達する時間にほぼ等しいことから、この時間からも決定可能である。
【0052】
また、目標値変更制御応答を規定する立上り時間σRは前述したように演算に より一義的に決定でき、或いは目標値ステップ変化時の制御量がほぼ50%に到達する時間として任意に設定できる。
【0053】
次に、図4(b)に示すパラメータ発生機構2は、比較要素13、除算要素14、乗算要素15、除算要素16および係数要素17が設けられている。この比較要素13および除算要素14は、図4(a)と同様に比例項用2自由度化パラメータbを算出する。さらに、乗算要素15において比較要素13の出力と除算要素14の出力とを乗算し、この乗算結果を除算要素16に導入し、ここで乗算結果を微分時間TDで除した後、係数要素17で1/2を乗算すれば、前記(1 2)式に基づく微分項用2自由度化パラメータc=(σD−σR)2/(2TITD)を算出でき、2自由度PID制御器1に発信する。
【0054】
さらに、図5は図3(b)のパラメータ発生機構2の内部構成を示す図である。
【0055】
このパラメータ発生機構2は、第1の演算要素21および第2の演算要素22からなり、第1の演算要素21は(10)式の演算により比例項用2自由度化パラメータbを算出し、2自由度PID制御器1に発信する。第2の演算要素22は(11)式の演算により微分項用2自由度化パラメータcを算出し、同様に2自由度PID制御器1に発信する。
【0056】
次に、パラメータ発生機構2の他の例について図6を参照して説明する。
【0057】
図6(a)は外乱入力および目標値変更による応答制御を規定する立上り時間σD、σRの差信号(σD−σR)を除算要素12に導き、ここで差信号を積分時間TIで除算することによって比例項用2自由度化パラメータbを算出する例であ る。
【0058】
図6(b)は同図(a)と同様に両立上り時間σD、σRの差信号(σD−σR)を入力する点を除けば、図4(b)と同様な構成により2自由度化パラメータb,cを算出し、2自由度PID制御器1に発信する例である。
【0059】
なお、以上のようにして決定された入力項目Aのもとに2自由度化パラメータb,cを設定した場合の目標値ステップ変更時の制御応答例は図7(b)に示す通りである。また、比較のため、従来のI−PD性御時(b=c=0)の制御応答は図7(c)に、基本PID制御時(b=c=1)の制御応答は図7(d)に示す通りである。
【0060】
これらの制御応答図から明らかなように、b=c=0または1に限定されることなく、b,cを任意の所定値とすることにより、目標値変更時に大きなオーバシュートを発生させることなく、追従性の優れた2自由度PID制御を実現できる。
【0061】
なお、σD、σR、TI、TDは以上記述した値に限られるものではなく、任意の定めた値を用いて、(10)式、(11)式の演算式により、2自由度化パラメータb,c決定し設定してもよい。
【0062】
従って、以上のような実施の形態によれば、予め定められた入力項目をパラメータ発生機構2に入力すれば、前述する導出された所定の演算式に基づいて2自由度化パラメータb,cを求めて2自由度PID制御器1に設定可能であり、従来のように長時間をかけて試行錯誤的に調整しつつ最適パラメータを決定するといった手間を省くことができ、調整作業の迅速化および熟練者によらずにパラメータを決定できる。
【0063】
(第2の実施の形態)
第1の実施の形態では、(1)式〜(11)式の導出過程に基づいて入力項目A(σD、σR 、T I 、T D 、σ2D 、σ2R )を選定し、これら入力項目Aのもとに所定の演算式により2自由度化パラメータb,cを決定したが、本実施の形態では、グラフ関数を用いて2自由度化パラメータb,cを発信させる構成とするものである。
【0064】
すなわち、本実施の形態においては、遅れ次数をもつ各種の制御対象に関し、前述する各式で求められたパラメータb,cをグラフで表わせば、図8(a)、(b)に示すようになるので、予めパラメータ発生機構2に図8に示すような関係にあるグラフを関数として設定しておき、外部から入力される制御対象の伝達関数(テスト・バッチ)タイプ、むだ時間Lおよび限界周期Tc、或いはLとTcとの代わりに基準化むだ時間LN(=L/Tc)を受けたとき、関数としての 2自由度化パラメータb,cを発信するようにしてもよい。
【0065】
図9および図10は本発明に係る2自由度PID制御装置の実施の形態を示す構成図である。
【0066】
この制御装置は図3と同様に2自由度PID制御器1とパラメータ発生機構2とで構成されているが、図9はパラメータ発生機構2に図8に示すようなグラフ関数が設定され、入力される制御対象タイプ、むだ時間Lおよび限界周期Tcに基づいて2自由度化パラメータb,cを発生させる例である。
【0067】
図10は同じくパラメータ発生機構2に図8に示すグラフ関数が設定され、入力される制御対象タイプおよび基準化むだ時間LN(L/Tc)に基づいて同じ く2自由度化パラメータb,cを発信し、2自由度PID制御器1に設定する例である。
【0068】
従って、この実施の形態によれば、各種の制御対象に関し、第1の実施の形態により決定されたパラメータb、cをグラフ関数として予めパラメータ発生機構2に設定しておけば、制御対象タイプ、むだ時間などの入力項目を入力するだけで、迅速に制御対象タイプに応じて最適な2自由度化パラメータb,cを発生することができ、煩雑な調整作業を必要とせずに目標値変更に対して追従性の優れた2自由度PID制御装置を実現できる。
(第3の実施の形態)
この実施の形態は、2自由度PID制御装置に代えて2自由度I−PD制御装置に適用した例である。I−PD制御は比例動作および微分動作とも制御量だけによって依存して動作する制御である。
【0069】
この2自由度I−PD制御系は、基本的には図20ないし図22と同様な構成であることから、2自由度I−PD制御装置においても概念的には図1と同様な構成である図11で表わすことができる。
【0070】
すなわち、2自由度I−PD制御装置は、2自由度I−PD制御器1′にパラメータ発生機構2が接続され、このパラメータ発生機構2は入力される入力項目Aに基づいて比例項用および微分項用2自由度化パラメータb,cを発生する機構である。
【0071】
この2自由度PID制御系は、図20の従来構成から図12(a)に示すようなブロック図で表すことができる。
【0072】
図12において31は2自由度I−PD制御器(2自由度PID制御器)、32は制御対象、33は加算要素、Rは制御量、Yは制御量、Dは外乱である。
【0073】
この2自由度I−PD制御器31は、目標値Rと制御量Yが入力され、外乱に対する抑制特性のためのPIDパラメータである比例ゲインKP、積分時間TI、微分時間TD、比例項用2自由度化パラメータbおよび微分項用2自由度化パラ メータcを用い、以下のような演算式により操作量MVを求める。
【0074】
KP{bR−Y+(1/TIs)(R−Y)+TDs(cR−Y)}……(13)
このようにして求めた操作量MVは外乱Dとともに制御対象32に印加する。
【0075】
ところで、図12(a)に示す制御系は等価ブロック変換すると、仮想的に図12(b)のようなフィードフォワード制御要素34とフィードバック制御要素35とで構成される。32は制御対象、36は加減算要素である。
【0076】
この制御系においてフィードフォワード制御要素34は、入力される目標値Rに対し、制御要素34の伝達関数を乗算し出力Wを発信する。また、フィードバック制御要素35は、制御対象32の出力である制御量Yに対し、制御要素35の伝達関数を演算し出力Vを発信する。そして、加減算要素36はフィードフォワード制御要素34の出力Wからフィードバック制御要素35の出力Vを減算し、得られた減算出力を制御対象32に印加し、制御量Yを発生させる。このとき、外乱Dは加減算要素36に加算的に加わり、制御対象32に入力される。
【0077】
また、制御対象32の伝達関数をGP(s)とすると、フィードバック制御系 における外乱Dから制御量Yへの伝達関数WD(s)は次式で表わせる。
【0078】
他方、外乱Dによる制御量Yの望ましい制御応答を規定する規範モデルは既に知られているように次式とする(PID制御、p14(2.11)式またはp31(2.35)式、著者:須田信英、編者:システム制御情報学会、発行所:株式会社 朝倉書店)。
【0079】
MD(s)={(TI/KP)s}/(1+σDs+α2DσD2s2+α3DσD3s3+α4DσD4s4) ……(15)
ここで、外乱入力に対する制御応答に関する規範モデルの係数α2D,α3D, α4Dは速応性とロバスト性との兼ね合いから、例えば表2に示す値とする。
【0080】
【表2】
【0081】
ここで、前述と同様にWD(s)=MD(s)として解くと、外乱入力に対する制御応答に関する立上り時間σDおよびPIDパラメータKP、TI、TDが決定できる。この決定されたパラメータを設定したI−PDフィードバック制御系に対し、制御対象の入力側にステップ状外乱を与えたとすると、図17に示すような 制御応答が得られる。
【0082】
次に、図12(b)に示す制御系に対し、決定された制御パラメータKP,TI,TDを設定し、このときの目標値Rから制御量Yへの伝達関数WR(s)をもって表わすと、次式のようになる。
【0083】
WR(s)=[KP{b+1/(TIs)+cTDs}] ・MD(s)……(16)
ここで、前段の[ ]内はフィードフォワード制御要素34の伝達関数、後段のMD(s)は前記(15)式の規範モデルの伝達関数である。
【0084】
従って、この(16)式は下記する(17)式で表わせる。
【0085】
WR(s)=(1+bTIs+cTITDs2)/(1+σDs+α2DσD2s2+α3DσD3s3+α4DσD4s4)……(17)
他方、目標値Rによる制御量Yの望ましい制御応答を規定する規範モデルは次式で表わすものとする(前述する文献 PID制御 、p120(6.2)式参照)。
【0086】
MR(s)=1/(1+σRs+α2RσR2s2+α3RσR3s3)…(18)
ここで、目標値変化に対する制御応答に関する規範モデルの係数α2R、α3Rは速応性とロバスト性との兼ね合いから、例えば前記表2に示す値とする。
【0087】
そこで、WR(s)=MR(s)とし、高次項を無視して解くと、目標値変更に対する制御応答に関する立上り時間σRと外乱入力に対する制御応答に関する立 上り時間σDとの立上り時間比σR/σDは下式の3次方程式の最小実根から求ま る。
【0088】
(1−2α2 R+α3 R)(σR/σD)3−(1−α2 R)(σR/σD)2+α2D(σR/σD)−α3 D=0 ……(19)
また、微分時間を用いない2自由度I−P制御では、立上り時間比σR/σDは次の2次方程式の最小実根として求まる。
【0089】
(1−α2R)(σR/σD)2−(σR/σD)+α2D=0 ……(20)
次に、比例項用2自由度化パラメータb算出用係数βは、
β=1−(σR/σD) ……(21)
また、微分項用2自由度化パラメータc算出用係数γは、
γ=(1−α2R)(σR/σD)2−(σR/σD)+α2D ……(22)
とする。
【0090】
よって、これらの値を用いると、比例項用2自由度化パラメータbは、
b=βσD/TI ……(23)
と決定でき、また微分項用2自由度化パラメータcは、
c=γσD2/(TITD) ……(24)
と決定できる。
【0091】
そこで、これら係数βおよび係数γを発生させる関数要素が必要となるが、これら係数を発生させる関数要素は以下のように作成する。
【0092】
すなわち、ある補間係数λDから定まる(15)式の規範モデルMD(s)と、ある補間係数λRから定まる前記(18)式の規範モデルMR(s)とを用いて、前記(19)式の方程式または前記(20)式の方程式から立上り時間比σR/ σDを求め、前記(21)式、(22)式から係数β、γを求める。
【0093】
そして、これら補間係数λD,λRの種々の値に対して係数β、γを求め、x軸をλD、y軸をλR、z軸をβおよびγとした行列データテーブルを作成し、補間係数λD,λRを入力すれば、係数β、γを発生する関数要素を作成できる。
【0094】
なお、補間係数λと最大感度Msとの間には図15に示す関係があるので、λD,λRの代わりに最大感度MSD、MSRを用いてもよい。また、2次の規範モデルを用いた場合には、λRの代わりとして、規範モデルの分母2次項の係数α2Rを用いてもよい。
【0095】
さらに、関数要素例として、例えば図16(a)に示すように、x軸に3次規範モデルのλD、y軸に2次規範モデルのα2Rとすれば、2自由度I−P制御用の係数βを求めるための関数が得られる。また、図16(b)、(c)に示すように、x軸に4次規範モデルのλD、y軸に3次規範モデルのλRとすれば、2自由度I−PD制御用の係数β、γを求めるための関数が得られる。
【0096】
図13は以上のような説明を踏まえつつ具体化したパラメータ発生機構の一例を示す構成図である。
【0097】
このパラメータ発生機構2は関数要素41と演算要素42とで構成されている。この関数要素41は、外乱抑制用I−PD制御系の規範モデルの属性値XD( λD、MSDなど)と目標値追従性を規定する規範モデルの属性値XR(γR、MSR 、α2Rなど)との種々の値の組みに対して係数βの値を蓄えた行列データテーブルであって、入力される属性値XDと属性値XRとに基づいて、行列データテーブルから蓄えた係数βを発生させる。なお、行列データテーブルの代わりに、前記(19)式または(20)式と(21)式との演算により、係数βを発生させる要素であってもよい。
【0098】
前記演算要素42は乗算要素43および除算要素44からなり、そのうち乗算要素43は関数要素41から発生される係数βに外乱抑制用I−PD制御系の規範モデルの立上り時間σDを乗算し、除算要素44に導入する。この除算要素4 4は、乗算要素43によって得られる乗算値を積分時間TIで除して比例項用2 自由度化パラメータbを算出し、2自由度I−PD制御器1′、31に設定する。
【0099】
図14は2自由度化パラメータ発生機構の他の例を示す構成図である。
【0100】
このパラメータ発生機構2は、2つの関数要素45、46と、これら関数要素45,46にそれぞれ対応する2つの演算要素47,48とによって構成されている。
【0101】
これら関数要素45、46は、それぞれ規範モデルの属性値XDとXRとの種々の値の組みに対して係数βの値を蓄えた行列データテーブルであって、入力される属性値XDと属性値XRとに基づいて、それぞれ行列データテーブルから蓄えた係数β,γを発生させる。なお、関数要素45である行列データテーブルの代わりに、前記(19)式または(20)式と(21)式との演算により係数βを発生させ、また関数要素46である行列データテーブルの代わりに、前記(19)式と(22)式の演算により係数γを発生させる要素であってもよい。
【0102】
前記演算要素47は図13に示す演算要素42と同様に乗算要素49および除算要素50からなり、乗算要素49は関数要素45から発生される係数βに外乱抑制用I−PD制御系の規範モデルの立上り時間σDを乗算し、除算要素50に 導入する。この除算要素50は、乗算要素49にて得られた乗算値を積分時間TIで除して比例項用2自由度化パラメータbを算出する。
【0103】
一方の演算要素48は、第1の乗算要素51、第2の乗算要素52、除算要素53および乗算要素54からなり、前記(24)式の演算により微分項用2自由度化パラメータcを算出する。具体的には、第2の乗算要素52にて得られた立上り時間σDの乗算値σD2と係数γとを乗算要素51にて乗算し、除算要素53に導入する。この除算要素53は乗算要素51からの乗算値を乗算要素54により得られる積分時間TIと微分時間TDとの乗算値で除算し、(24)式により微分項用2自由度化パラメータcを算出する。そして、各演算要素47,48で得られたパラメータb,cを例えば図20に示す2自由度I−PD制御器1′に発信する。
【0104】
なお、2自由度I−PD制御器1′としては、図20に示したものに限定されず、これを等価ブロック変換して得られる図21に示す目標値フィードフォワード型や図22に示す目標値フィルタ型、さらにはループ補償型やフィードバック補償型などであってもよい(PID制御、p75、著者:須田信英、朝倉書店)。
【0105】
また、図20に示す微分要素118として完全微分型で表記したが、1次遅れ要素が付加された不完全微分型であってもよい。
【0106】
次に、パラメータ発生機構への入力について説明する。外乱を最も良好に抑制できるように規範モデルを用いてPID値を設計した場合、この時に使用した規範モデルで補間係数λD、最大感度MSD、分母2次項の係数α2R、立上り時間σDは一義的に決まるので、これらを使用する。規範モデルを用いずにPID値を調整した場合、操作端にステップ入力を加えたときの補間係数λDを種々変えて描 いた制御応答(図17)と比較し、形状がもっとも近い波形の補間係数λDから 、規範モデルおよび立上り時間σD(制御応答がほぼピークに達する時間)を同 定して使用する。
【0107】
なお、図17(a)は4次の規範モデルの制御応答図であって、微分動作を用いた2自由度I−PD制御系を用いる場合に使用される。また、図17(b)は3次の規範モデルの制御応答図であって、微分動作を用いない2自由度I−P制御系に使用される。
【0108】
また、目標値変更の過度時に発生する操作量の変化を小さくする場合やロバスト性を増やしたい場合には立上り時間σDとして前記により定まる値よりも大き い値を、また応答を早めたい場合には小さい値を入力する。
【0109】
また、例えばI−PD制御系の規範モデルを自動的に同定する機構を付加し、I−PD制御器で用いられている積分時間TI、微分時間TDとともに同定結果の入力項目を本発明のパラメータ発生機構に自動的に入力すれば、オート(セルフ)チューニング2自由度I−PD制御系を構成することができる。
【0110】
さらに、目標値変更の追従性を規定する規範モデルとしては、補間係数λRや 分母2次項の係数α2Rを種々変えて描いた制御応答波形(図18)を参照し、オーバシュート量などから、その制御系で最も適切な形状を与える補間係数λRや 分母2次項の係数α2R、つまり規範モデルの属性値を決定して入力する。なお、図18(a)は目標値単位ステップ変更時の補間係数λRを種々変えた場合の3 次規範モデルの制御応答図であって、2自由度I−PD制御に用いられ、図18(b)は目標値単位ステップ変更時の分母2次項の係数を種々変えた場合の2次の規範モデルの制御応答図であって、2自由度I−P制御時に用いられる。
【0111】
図19は、本発明によるパラメータ発生機構2から2自由度化パラメータを発生し、2自由度I−P(D)制御器に設定し、制御対象を制御したときの制御応答の一例である。なお、比較のために、1自由度I−P(D)制御器およびPI(D)制御器で制御した場合についてその制御応答を示している。
【0112】
この図から明らかなように、本発明によるI−PD制御の場合、目標値変更時に大きなオーバシュートを発生させることなく、しかも速応性をもたせて応答させることができ、目標値追従性を大幅に向上させることができる。
【0113】
【発明の効果】
以上説明したように本発明によれば、所定の入力項目のもとに目標値変更時の制御応答を規定する2自由度化パラメータを自動的に発生でき、従来のように長時間かけて試行錯誤的に調整する作業を省略でき、また熟練者に頼らずに短時間に2自由度化パラメータを発生できる。
【0114】
よって、外乱に対しても目標値変更に対しても、迅速に最適な制御性を実現できる2自由度制御装置を提供できる。
【図面の簡単な説明】
【図1】本発明に係わる2自由度制御装置の第1および第2の実施の形態を説明する概念構成図。
【図2】本発明の原理を説明する2自由度PID制御系のブロック線図。
【図3】本発明に係わる2自由度制御装置における入力項目を含む全体のブロック構成図。
【図4】図3に示すパラメータ発生機構の具体的な構成例を示す図。
【図5】図3に示すパラメータ発生機構の具体的な構成例を示す図。
【図6】図3に示すパラメータ発生機構の具体的な構成例を示す図。
【図7】本発明に係わる2自由度制御装置を用いた場合の制御応答例を示す図。
【図8】図3に示すパラメータ発生機構に設定する関数を説明する図。
【図9】本発明に係わる2自由度制御装置の第2の実施形態を示すブロック構成図。
【図10】本発明に係わる2自由度制御装置の第2の実施形態の他の例を示すブロック構成図。
【図11】本発明に係わる2自由度制御装置の第3の実施の形態を説明する概念構成図。
【図12】本発明の原理を説明する2自由度I−PD制御系のブロック線図。
【図13】パラメータ発生機構の具体的な構成例を示す図。
【図14】パラメータ発生機構の具体的な構成例を示す図。
【図15】規範モデルの補間係数と最大感度との関係を示す図。
【図16】関数要素の関数特性を示す図。
【図17】PID値設計あるいはI−PD制御系の特性同定に使用するステップ外乱入力時の規範モデルの制御応答を示す図。
【図18】2自由度I−P(D)制御設計に使用する目標値単位ステップ変更時の規範モデルの制御応答を示す図。
【図19】本発明に係わる2自由度制御装置を用いた場合の制御応答例を示す図。
【図20】従来公知の一般的な2自由度PIDまたは2自由度I−PD制御系の概略構成を示す図。
【図21】従来公知のフィードフォワード要素をもつ2自由度制御装置の構成図。
【図22】従来公知の目標値フィルタ要素をもつ2自由度制御装置の構成図。
【符号の説明】
1…2自由度PID制御器
1′…2自由度I−PD制御器
2…パラメータ発生機構
21、22…演算要素
σD…外乱入力制御応答に関する立上り時間
σR…目標値変更制御応答に関する立上り時間
TI…積分時間
TD…微分時間
α2D…外乱入力制御応答規範モデルの2次項係数
α2R…目標値変更制御応答規範モデルの2次項係数
41,45,46…関数要素
42,47,48…演算要素[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a two-degree-of-freedom controller used in various instrumentation fields and servo control fields of petroleum, chemical, steel, electric power, papermaking, and the like, and particularly relates to a flow rate, a pressure, a liquid level, and a composition of each part of a plant. The present invention relates to a two-degree-of-freedom control device that performs feedback control of process variables such as position, speed, and rotation speed.
[0002]
[Prior art]
Conventionally, a two-degree-of-freedom PID control device has been used to satisfactorily perform feedback control of a control amount both when a disturbance is input and when a target value is changed.
[0003]
Here, the two-degree-of-freedom PID control is PID control for simultaneously optimizing both the suppression characteristic for disturbance and the follow-up characteristic for target value change.
[0004]
FIG. 20 is a block diagram of a conventional general 2-DOF PID or 2-DOF I-PD control device.
[0005]
This control device includes a two-degree-of-freedom PID (hereinafter, including two-degree-of-freedom I-PD)
[0006]
The proportional system multiplies the target value R by the two-degree-of-freedom parameter b for the proportional term of the
[0007]
The output of the
[0008]
Conventionally, as shown in FIG. 21, a feedforward control method is adopted in which a change in a target value is detected, a change in an operation amount is determined in advance in order to cancel the influence thereof, and an operation is executed ahead of schedule. A degree of freedom PID controller is used.
[0009]
The two-degree-of-
[0010]
Further, conventionally, as shown in FIG. 13, a two-degree-of-freedom PID control device to which a target
[0011]
The two-degree-of-
[0012]
(1 + bTIs + cTITDs2) / (1 + T)Is) …… (1)
The control target E and the control amount Y obtained as described above are used to determine the control deviation E by the
[0013]
[Problems to be solved by the invention]
By the way, in each of the various PID controllers described above, the PID parameter (proportional gain KP, Integration time TI, Differentiation time TD), It is necessary to adjust the two-degree-of-freedom parameters b and c for optimizing the follow-up characteristic to the change of the target value by trial and error. In addition, a control target model is constructed in advance, a simulation is executed under a certain evaluation function offline, and each parameter of the transfer function, that is, a PID parameter and two-degree-of-freedom parameters b and c are designed by searching for an optimum value. Has been done.
[0014]
As a result, in order to obtain the PID parameter and the two-degree-of-freedom parameters b and c that provide an optimal control response to disturbance and a change in the target value, a complicated adjustment work by a skilled person is necessary, and the adjustment work is required. There is a problem that requires a lot of time.
[0015]
SUMMARY OF THE INVENTION The present invention has been made in view of the above circumstances, and provides a two-degree-of-freedom control device that quickly and easily generates a two-degree-of-freedom parameter for obtaining an optimal control response to a change in a target value without skill. Is to provide.
[0016]
[Means for Solving the Problems]
In order to solve the above-described problems, the present invention provides a two-degree-of-freedom control device that uses a control parameter and a two-degree-of-freedom parameter to control both a disturbance suppression property and a target value tracking property with respect to a control target. A rise time σD relating to a disturbance input control response, a rise time σR relating to a target value change control response, and an integral time T which is the control parameter.IAlternatively, the integration time TIAnd derivative time TDAnd a parameter generating mechanism for generating the two-degree-of-freedom parameter defining the target value change control response.
[0017]
In addition, in addition to the above input items, a value proportional to a value obtained by a predetermined arithmetic expression using a second-order coefficient α2D of the disturbance input control response reference model and a second-order coefficient α2R of the target value change control response reference model. , A two-degree-of-freedom parameter defining the target value change control response.
[0018]
According to the present invention, by taking the above-described means, predetermined predetermined input items, that is, a rise time σD relating to a disturbance input control response, a rise time σR relating to a target value change control response, and an integral Time TIAlternatively, the integration time TIAnd derivative time TDAnd α2D as a second-order coefficient of the disturbance input control response reference model and α2R as a second-order coefficient of the target value changing control response reference model, and a predetermined value derived from the transfer function of the actual control system and the transfer function of the reference model. By executing the above calculation, the two-degree-of-freedom parameter defining the target value change control response can be generated, and the work of adjusting by trial and error as in the related art can be omitted.
[0019]
Another invention provides the two-degree-of-freedom parameter defining a target value change control response using a transfer characteristic type of a controlled object, a dead time and a limit period, or a transfer characteristic type and a standardized dead time of a controlled object. Is provided with a parameter generating mechanism for generating the
[0020]
According to the present invention, by taking the above-described means, it is possible to previously form a graph function based on the transfer characteristic type of the controlled object, the dead time and the limit period, or the transfer characteristic type and the standardized dead time of the controlled object. It is possible to generate the two-degree-of-freedom parameter quickly and easily, and to achieve the same effect as the above-mentioned invention.
[0021]
Further, in another invention, the attribute value of the first reference model and the attribute value of the second reference model that defines the target value followability in the disturbance suppression I-PD control system are input, and both of these attribute values are input. A function element for generating a coefficient value as a function, and a calculation element for generating a two-degree-of-freedom parameter for improving target value followability using the coefficient value of the function element, the rise time and the integration time of the first reference model And a parameter generating mechanism having the following.
[0022]
According to the present invention, a coefficient value is generated as a function of an attribute value of a first reference model and an attribute value of a second reference model defining target value followability in a disturbance suppression I-PD control system. , Using the rise time and the integration time of the first reference model to generate the two-degree-of-freedom parameter for improving the followability of the target value, the two-degree-of-freedom parameter is automatically adjusted without performing trial and error adjustment work. Parameter can be generated.
[0023]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, embodiments of the present invention will be described with reference to the drawings.
(First Embodiment)
Prior to describing an embodiment of the two-degree-of-freedom control apparatus according to the present invention, a basic concept for realizing the present apparatus will be described with reference to FIG.
(1) First, based on the input of the target value R and the control amount Y, the proportional gain KP, Integration time TI, Differentiation time TDIn addition to the above, a two-degree-of-
(2) Next, in the two-degree-of-
(3) Further, a new
(4) It is assumed that the two-degree-of-freedom parameters b and c calculated by the
[0024]
In order to concretely realize the above concept, it is important what kind of input item A is used, and what kind of calculation the
[0025]
Therefore, a procedure for deriving the input item A and the arithmetic expression will be described below.
[0026]
If a conventional general two-degree-of-freedom control system (see FIG. 20) is represented by a transfer function, a configuration as shown in FIG. 2A is obtained.
[0027]
In the figure, 1 is a two-degree-of-freedom PID controller, 3 is a control target, 4 is an addition element, R is a target value, Y is a control amount, and D is a disturbance.
[0028]
That is, the two-degree-of-
[0029]
KP{BR-Y + (1 / TIs) (RY) + TDs (cR-Y)} (2)
The manipulated variable MV obtained by the above equation is applied to the
[0030]
Incidentally, the two-degree-of-freedom PID control system shown in FIG. 2A can be virtually represented by a control system as shown in FIG. 5 is a feedforward control element, 6 is an addition / subtraction element, and 7 is a feedback control element.
[0031]
In this control system, the
[0032]
Further, the transfer function of the
[0033]
On the other hand, a reference model that defines a desired control response of the control amount Y due to the disturbance D is expressed by the following equation (PID control, equation p14 (2.11) or equation p31 (2.35), author). : Nobuhide Suda, Editor: The Society of System Control and Information Engineers, Publisher: Asakura Shoten Co., Ltd.).
[0034]
MD(S) = {s / (KP/ TI)} · [1 / (1 + σDs + α2DσD)2s2+ Α3DσD3s3)] …… (4)
Here, the coefficients α2D and α3D of the reference model regarding the control response to the disturbance input are set to, for example, the values shown in Table 1 in consideration of the balance between the quick response and the robustness. In Table 1, see the above-mentioned document PID control, p14 to p15.
[0035]
[Table 1]
[0036]
Incidentally, the actual control system of the above equation (3) and the reference model of the equation (4) must be in the same relation. Therefore, WD(S) = MDSolving for (s) gives rise time σ for control response to disturbance inputDAnd PID parameter KP, TI, TD ButCan decide. Therefore, if a step-like disturbance is applied to the input side of the control target in the PID feedback control system in which the determined parameters are set, a control response as shown in FIG. 7A is obtained.
[0037]
Next, for the control system shown in FIG.P, TI, TD, and the transfer function W from the target value R to the control amount Y at this time is set.RWhen expressed by (s), the following equation is obtained.
[0038]
WR(S) = [KP{B + 1 / (TIs) + cTDs}] ・ MD(5) Here, [] in the preceding stage is the transfer function of the
[0039]
Therefore, this equation (5) can be expressed by the following equation (6).
[0040]
WR(S) = (1 + bT)Is + cTITDs2) / (1 + σDs + α2DσD)2s2+ Α3DσD3s3) …… (6)
On the other hand, a reference model that defines a desirable control response of the control amount Y based on the target value R is represented by the following equation (see the above-mentioned document PID control, p120 (6.2) equation).
[0041]
MR(S) = 1 / (σRs + α2RσR2s2+ Α3RσR3s3…… (7)
Here, the coefficients α2R and α3R of the reference model relating to the control response to the change in the target value are set to the values shown in Table 1 above, for example, in view of the balance between quick response and robustness.
[0042]
Then, WR(S) = MRWhen solving for (s), the rise time σR relating to the control response to the change in the target value is obtained from the minimum real root of the following cubic equation.
[0043]
(1-2α2R + α3R) σR3− (1-α2R) σDσR2+ Α2DσD2σR-α3DσD3= 0 (8)
Here, W in equation (6)R(S) = B / C, and M in equation (7)RIf (s) = 1 / Z, WR(S) = MRFrom the relationship (s),
B / C = 1 / Z → Z = C / B (9)
Can be derived.
[0044]
Therefore, by performing an operation on C / B, deriving a solution up to the first-order and second-order terms, and calculating the Z and B terms in the same relation as the first-order and second-order terms, respectively, the two-degree-of-freedom parameter b for the proportional term becomes
b = (σD−σR) / TI …… (10)
It can be determined from the following arithmetic expression. The two-degree-of-freedom parameter c for the differential term is
It can be determined from the following arithmetic expression.
[0045]
In the case of α2D = α2R = 0.5 in the equation (11),
c = (σD-σR)2/ (2TITD) (12) and the two-degree-of-freedom parameters b and c are (σD−σR).
[0046]
Therefore, as is apparent from the above description, the input item A and the arithmetic expression for obtaining the two-degree-of-freedom parameters b and c were obtained, and a specific example will be described below.
[0047]
FIG. 3 is a configuration diagram showing one embodiment of a two-degree-of-freedom PID control device according to the present invention.
[0048]
In the control device shown in FIG. 1A, a
[0049]
The control device shown in FIG. 2B further includes, in addition to the input items shown in FIG.
A quadratic term coefficient α2D of the disturbance input control response reference model and a quadratic term coefficient α2R of the target value change control response reference model are additionally input, and the calculation means shown in FIG. And the two-degree-of-
[0050]
That is, the
[0051]
Note that the rise time σD that defines the disturbance input control response can be automatically determined when the PID parameter is designed by the partial model matching method, or when the PID parameter is determined by using another method, the step disturbance can be determined. Since the time when the control amount Y reaches the peak value upon entering, it can be determined from this time.
[0052]
Further, the rise time σR defining the target value change control response can be uniquely determined by the calculation as described above, or can be arbitrarily set as the time required for the control amount at the time of the target value step change to reach approximately 50%.
[0053]
Next, the
[0054]
FIG. 5 is a diagram showing the internal configuration of the
[0055]
The
[0056]
Next, another example of the
[0057]
FIG. 6A shows a difference signal (σD−σR) between the rise times σD and σR that defines the response control by disturbance input and target value change to the
[0058]
FIG. 6 (b) has two degrees of freedom by a configuration similar to that of FIG. 4 (b), except that a difference signal (σD−σR) between the rise times σD and σR is input as in FIG. This is an example in which parameters b and c are calculated and transmitted to the two-degree-of-
[0059]
FIG. 7B shows an example of a control response when the target value step is changed when the two-degree-of-freedom parameters b and c are set based on the input item A determined as described above. . For comparison, FIG. 7C shows a control response in the conventional I-PD control (b = c = 0), and FIG. 7C shows a control response in the basic PID control (b = c = 1). As shown in d).
[0060]
As is apparent from these control response diagrams, b and c are not limited to 0 or 1, and by setting b and c to arbitrary predetermined values, a large overshoot does not occur when the target value is changed. Thus, two-degree-of-freedom PID control with excellent followability can be realized.
[0061]
Note that σD, σR, TI, TDIs not limited to the values described above, and the two-degree-of-freedom parameters b and c may be determined and set according to the arithmetic expressions of Expressions (10) and (11) using any predetermined values. .
[0062]
Therefore, according to the above-described embodiment, if a predetermined input item is input to the
[0063]
(Second embodiment)
In the first embodiment, based on the derivation process of the equations (1) to (11), the input item A (σD, ΣR , T I , T D ,σ2D ,σ2R )Is selected, and the two-degree-of-freedom parameters b and c are determined by a predetermined arithmetic expression based on these input items A. In the present embodiment, the two-degree-of-freedom parameters b and c are determined using a graph function. It is configured to transmit.
[0064]
That is, in the present embodiment, regarding the various control objects having the delay order, the parameters b and c obtained by the above-described equations can be represented by graphs as shown in FIGS. Therefore, a graph having a relationship as shown in FIG. 8 is set in advance in the
[0065]
9 and 10 are configuration diagrams showing an embodiment of a two-degree-of-freedom PID control device according to the present invention.
[0066]
This control device comprises a two-degree-of-
[0067]
In FIG. 10, the graph function shown in FIG. 8 is similarly set in the
[0068]
Therefore, according to this embodiment, if the parameters b and c determined according to the first embodiment are set in advance in the
(Third embodiment)
This embodiment is an example in which the present invention is applied to a 2-DOF I-PD control device instead of a 2-DOF PID control device. The I-PD control is a control that operates depending on only the control amount in both the proportional operation and the differential operation.
[0069]
Since the two-degree-of-freedom I-PD control system has basically the same configuration as in FIGS. 20 to 22, the two-degree-of-freedom I-PD control device also has a conceptually similar configuration to FIG. This can be represented in FIG.
[0070]
That is, in the two-degree-of-freedom I-PD control device, a
[0071]
This two-degree-of-freedom PID control system can be represented by a block diagram as shown in FIG. 12A from the conventional configuration of FIG.
[0072]
In FIG. 12, 31 is a two-degree-of-freedom I-PD controller (two-degree-of-freedom PID controller), 32 is a control target, 33 is an addition element, R is a control amount, Y is a control amount, and D is a disturbance.
[0073]
The two-degree-of-freedom I-
[0074]
KP{BR-Y + (1 / TIs) (RY) + TDs (cR-Y)} (13)
The operation amount MV obtained in this way is applied to the
[0075]
By the way, the control system shown in FIG. 12A is virtually composed of a
[0076]
In this control system, the
[0077]
Further, the transfer function of the controlled
[0078]
On the other hand, a reference model that defines a desired control response of the control amount Y due to the disturbance D is expressed by the following equation (PID control, equation p14 (2.11) or equation p31 (2.35); : Nobuhide Suda, Editor: The Society of Systems, Control and Information Engineers, Publisher: Asakura Shoten Co., Ltd.)
[0079]
MD(S) = {(TI/ KP) S} / (1 + σDs + α2DσD)2s2+ Α3DσD3s3+ Α4DσD4s4) ...... (15)
Here, the coefficients α2D, α3D, and α4D of the reference model relating to the control response to the disturbance input are set to, for example, the values shown in Table 2 in consideration of the balance between the quick response and the robustness.
[0080]
[Table 2]
[0081]
Here, as described above, WD(S) = MD(S), the rise time .sigma.D and the PID parameter K for the control response to the disturbance inputP, TI, TDCan be determined. Assuming that a step-like disturbance is applied to the input side of the control target in the I-PD feedback control system in which the determined parameters are set, a control response as shown in FIG. 17 is obtained.
[0082]
Next, with respect to the control system shown in FIG.P, TI, TDIs set, and a transfer function W from the target value R to the control amount Y at this time is set.RWhen expressed by (s), the following equation is obtained.
[0083]
WR(S) = [KP{B + 1 / (TIs) + cTDs}] · MD (s) ... (16)
Here, [] in the first stage is the transfer function of the
[0084]
Therefore, this equation (16) can be expressed by the following equation (17).
[0085]
WR(S) = (1 + bT)Is + cTITDs2) / (1 + σDs + α2DσD)2s2+ Α3DσD3s3+ Α4DσD4s4) ... (17)
On the other hand, a reference model that defines a desired control response of the control amount Y based on the target value R is represented by the following equation (see the above-mentioned document PID control, p120 (6.2) equation).
[0086]
MR(S) = 1 / (1 + σRs + α2RσR2s2+ Α3RσR3s3)… (18)
Here, the coefficients α2R and α3R of the reference model relating to the control response to the change in the target value are set to the values shown in Table 2 above, for example, in view of the balance between quick response and robustness.
[0087]
Then, WR(S) = MR(S) and solving with ignoring higher order terms, the rise time ratio σR / σD of the rise time σR relating to the control response to the change of the target value and the rise time σD relating to the control response to the disturbance input is represented by the following cubic equation: From the smallest real root of.
[0088]
(1-2α2 R+ Α3 R) (ΣR/ ΣD)3− (1-α2 R) (ΣR/ ΣD)2+ Α2D(ΣR/ ΣD) -Α3 D= 0 (19)
In the two-degree-of-freedom IP control without using the differential time, the rise time ratio σR / σD is obtained as the minimum real root of the following quadratic equation.
[0089]
(1-α2R) (σR / σD)2− (ΣR / σD) + α2D = 0 (20)
Next, the coefficient β for calculating the two-degree-of-freedom parameter b for the proportional term is
β = 1− (σR / σD) (21)
The coefficient γ for calculating the two-degree-of-freedom parameter c for the differential term is
γ = (1−α2R) (σR / σD)2− (ΣR / σD) + α2D (22)
And
[0090]
Therefore, using these values, the two-degree-of-freedom parameter b for the proportional term is
b = βσD / TI ...... (23)
And the two-degree-of-freedom parameter c for the differential term is
c = γσD2/ (TITD) ...... (24)
Can be determined.
[0091]
Therefore, a function element for generating these coefficients β and γ is required. The function elements for generating these coefficients are created as follows.
[0092]
That is, using the reference model MD (s) of the equation (15) determined from a certain interpolation coefficient λD and the reference model MR (s) of the above equation (18) determined from a certain interpolation coefficient λR, The rise time ratio σR / σD is obtained from the equation or the equation (20), and the coefficients β and γ are obtained from the equations (21) and (22).
[0093]
Then, coefficients β and γ are obtained for various values of the interpolation coefficients λD and λR, and a matrix data table is created in which the x-axis is λD, the y-axis is λR, and the z-axis is β and γ. , ΛR, a function element for generating the coefficients β and γ can be created.
[0094]
Since there is a relationship shown in FIG. 15 between the interpolation coefficient λ and the maximum sensitivity Ms, the maximum sensitivities MSD and MSR may be used instead of λD and λR. When a second-order normative model is used, a coefficient α2R of a second-order denominator of the normative model may be used instead of λR.
[0095]
Further, as an example of a function element, as shown in FIG. 16A, if the x-axis is λD of the third-order norm model and the y-axis is α2R of the second-order norm model, for two-degree-of-freedom IP control, A function for obtaining the coefficient β is obtained. Further, as shown in FIGS. 16B and 16C, if the x-axis is λD of the fourth-order norm model and the y-axis is λR of the third-order norm model, the coefficient β for two degrees of freedom I-PD control , Γ are obtained.
[0096]
FIG. 13 is a configuration diagram showing an example of a parameter generation mechanism embodied based on the above description.
[0097]
This
[0098]
The
[0099]
FIG. 14 is a configuration diagram showing another example of the two-degree-of-freedom parameter generation mechanism.
[0100]
The
[0101]
These
[0102]
The
[0103]
One
[0104]
Note that the two-degree-of-freedom I-PD controller 1 'is not limited to the one shown in FIG. 20, but a target value feed-forward type shown in FIG. It may be a value filter type, or a loop compensation type or a feedback compensation type (PID control, p75, author: Nobuhide Suda, Asakura Shoten).
[0105]
Although the
[0106]
Next, the input to the parameter generating mechanism will be described. When the PID value is designed using the reference model so that disturbance can be suppressed most satisfactorily, the interpolation coefficient λD, the maximum sensitivity MSD, the coefficient α2R of the second-order denominator, and the rise time σD are uniquely determined in the reference model used at this time. Use these as they are decided. When the PID value is adjusted without using the reference model, it is compared with the control response (Fig. 17) drawn by changing the interpolation coefficient λD when a step input is applied to the operation terminal (Fig. 17), and the interpolation coefficient of the waveform having the closest shape is obtained. From λD, the reference model and the rise time σD (the time when the control response almost reaches a peak) are specified and used.
[0107]
FIG. 17A is a control response diagram of a fourth-order reference model, and is used when a two-degree-of-freedom I-PD control system using a differential operation is used. FIG. 17B is a control response diagram of a third-order reference model, which is used in a two-degree-of-freedom IP control system that does not use a differential operation.
[0108]
In addition, when the change in the manipulated variable that occurs when the target value change is excessive is reduced or when the robustness is desired to be increased, the rise time σD is set to a value larger than the value determined as described above. Enter a smaller value.
[0109]
Further, for example, a mechanism for automatically identifying a reference model of the I-PD control system is added, and an integration time T used in the I-PD controller is used.I, Differentiation time TDIf the input items of the identification result are automatically input to the parameter generating mechanism of the present invention, an auto (self) tuning two-degree-of-freedom I-PD control system can be configured.
[0110]
Further, as a reference model for defining the followability of the target value change, a control response waveform (FIG. 18) drawn by variously changing the interpolation coefficient λR and the coefficient α2R of the denominator quadratic term is referred to. The interpolation coefficient λR that gives the most appropriate shape in the control system and the coefficient α2R of the denominator quadratic term, that is, the attribute values of the reference model are determined and input. FIG. 18A is a control response diagram of the third-order model when the interpolation coefficient λR is variously changed when the target value unit step is changed, and is used for two-degree-of-freedom I-PD control. b) is a control response diagram of the second-order normative model when the coefficient of the denominator quadratic term is variously changed at the time of changing the target value unit step, and is used in two-degree-of-freedom IP control.
[0111]
FIG. 19 is an example of a control response when a two-degree-of-freedom parameter is generated from the
[0112]
As is clear from this figure, in the case of the I-PD control according to the present invention, it is possible to respond with a quick response without generating a large overshoot at the time of changing the target value, and it is possible to greatly improve the target value following ability. Can be improved.
[0113]
【The invention's effect】
As described above, according to the present invention, a two-degree-of-freedom parameter that defines a control response when a target value is changed can be automatically generated based on a predetermined input item. The work of adjusting by mistake can be omitted, and the two-degree-of-freedom parameter can be generated in a short time without relying on a skilled person.
[0114]
Therefore, it is possible to provide a two-degree-of-freedom control device capable of quickly realizing optimum controllability both for disturbance and target value change.
[Brief description of the drawings]
FIG. 1 is a conceptual configuration diagram illustrating first and second embodiments of a two-degree-of-freedom control device according to the present invention.
FIG. 2 is a block diagram of a two-degree-of-freedom PID control system for explaining the principle of the present invention.
FIG. 3 is an overall block diagram including input items in a two-degree-of-freedom control device according to the present invention.
FIG. 4 is a diagram showing a specific configuration example of a parameter generation mechanism shown in FIG. 3;
FIG. 5 is a diagram showing a specific configuration example of a parameter generation mechanism shown in FIG. 3;
FIG. 6 is a diagram showing a specific configuration example of a parameter generation mechanism shown in FIG. 3;
FIG. 7 is a diagram showing an example of a control response when the two-degree-of-freedom control device according to the present invention is used.
FIG. 8 is a view for explaining functions set in the parameter generation mechanism shown in FIG. 3;
FIG. 9 is a block diagram showing a two-degree-of-freedom control device according to a second embodiment of the present invention.
FIG. 10 is a block diagram showing another example of the two-degree-of-freedom control apparatus according to the second embodiment of the present invention.
FIG. 11 is a conceptual diagram illustrating a two-degree-of-freedom control apparatus according to a third embodiment of the present invention.
FIG. 12 is a block diagram of a two-degree-of-freedom I-PD control system illustrating the principle of the present invention.
FIG. 13 is a diagram showing a specific configuration example of a parameter generation mechanism.
FIG. 14 is a diagram showing a specific configuration example of a parameter generation mechanism.
FIG. 15 is a diagram showing a relationship between an interpolation coefficient of a reference model and a maximum sensitivity.
FIG. 16 is a diagram showing function characteristics of function elements.
FIG. 17 is a diagram showing a control response of a reference model at the time of a step disturbance input used for PID value design or characteristic identification of an I-PD control system.
FIG. 18 is a diagram showing a control response of a reference model when changing a target value unit step used in two-degree-of-freedom IP (D) control design.
FIG. 19 is a diagram showing an example of a control response when the two-degree-of-freedom control device according to the present invention is used.
FIG. 20 is a diagram showing a schematic configuration of a conventionally known general two-degree-of-freedom PID or two-degree-of-freedom I-PD control system;
FIG. 21 is a configuration diagram of a two-degree-of-freedom control device having a conventionally known feedforward element.
FIG. 22 is a configuration diagram of a two-degree-of-freedom control device having a conventionally known target value filter element.
[Explanation of symbols]
1 ... 2 degrees of freedom PID controller
1 '... 2-DOF I-PD controller
2. Parameter generation mechanism
21, 22 ... operation element
σD: Rise time for disturbance input control response
σR: Rise time for target value change control response
TI… Integration time
TD… Derivative time
α2D: second order coefficient of disturbance input control response reference model
α2R: second order coefficient of target value change control response reference model
41, 45, 46 ... function elements
42, 47, 48 ... arithmetic elements
Claims (10)
外乱入力制御応答に関する立上り時間と、目標値変更制御応答に関する立上り時間と、前記制御パラメータである積分時間あるいは当該積分時間および微分時間とを用いて、前記目標値変更制御応答を規定する前記2自由度化パラメータを発生するパラメータ発生機構を備えたことを特徴とする2自由度制御装置。In a two-degree-of-freedom control device that controls both disturbance suppression and target value tracking with respect to a controlled object using a control parameter and a two-degree-of-freedom parameter,
Using the rise time related to the disturbance input control response, the rise time related to the target value change control response, and the integration time or the integration time and the derivative time that are the control parameters, the two freedoms that define the target value change control response are used. A two-degree-of-freedom control device comprising a parameter generation mechanism for generating a degree-of-degree parameter.
外乱入力制御応答に関する立上り時間をσD、外乱入力制御応答規範モデルの 2次項係数をα2D、目標値変更制御応答に関する立上り時間をσR、目標値変更 制御応答規範モデルの2次項係数をα2R、積分時間をTI、微分時間をTDとしたとき、
{(1−α2R)・σR2−σD・σR+α2D・σD2}/TI・TD
なる演算式により得られる値に比例する値を、前記目標値変更制御応答を規定する2自由度化パラメータとする演算手段を設けたことを特徴とする2自由度制御装置。The parameter generating mechanism according to claim 1,
The rise time relating to the disturbance input control response is σD, the second order coefficient of the disturbance input control response reference model is α2D, the rise time relating to the target value change control response is σR, the target value change The second order coefficient of the control response reference model is α2R, and the integration time Is T I and the differentiation time is T D ,
{(1−α2R) · σR 2 −σD · σR + α2D · σD 2 } / T I · T D
A two-degree-of-freedom control device, comprising: a calculating means for setting a value proportional to a value obtained by an arithmetic expression as a two-degree-of-freedom parameter defining the target value change control response.
制御対象の伝達特性タイプ、むだ時間および限界周期、または制御対象の伝達特性タイプおよび基準化むだ時間を用いて、目標値変更制御応答を規定する前記2自由度化パラメータを発生するパラメータ発生機構を備えたことを特徴とする2自由度制御装置。In a two-degree-of-freedom control device that controls both disturbance suppression and target value tracking with respect to a controlled object using a control parameter and a two-degree-of-freedom parameter,
A parameter generation mechanism that generates the two-degree-of-freedom parameter that defines a target value change control response using the transfer characteristic type of the controlled object, the dead time and the limit period, or the transfer characteristic type and the standardized dead time of the controlled object. A two-degree-of-freedom control device, comprising:
各種の制御対象に関し、請求項1ないし請求項3において得られる目標値変更制御応答を規定する2自由度化パラメータをグラフ関数として設定してなることを特徴とする2自由度制御装置。The parameter generating mechanism according to claim 3,
4. A two-degree-of-freedom control device, wherein a two-degree-of-freedom parameter defining a target value change control response obtained in any one of claims 1 to 3 is set as a graph function for various control targets.
外乱抑制用I−PD制御系における第1の規範モデルの属性値と目標値追従性を規定する第2の規範モデルの属性値とが入力され、これら両属性値の関数として係数値を発生する関数要素と、この関数要素の係数値、第1の規範モデルの立上り時間および積分時間を用いて目標値追従性を良好にする2自由度化パラメータを発生する演算要素とを有するパラメータ発生機構を備えたことを特徴とする2自由度制御装置。In a two-degree-of-freedom control device that controls both disturbance suppression and target value tracking with respect to a controlled object using a control parameter and a two-degree-of-freedom parameter,
The attribute value of the first reference model and the attribute value of the second reference model that defines the target value followability in the disturbance suppression I-PD control system are input, and a coefficient value is generated as a function of these two attribute values. A parameter generation mechanism having a function element and an operation element for generating a two-degree-of-freedom parameter for improving target value followability by using a coefficient value of the function element, a rise time and an integration time of the first reference model. A two-degree-of-freedom control device, comprising:
前記関数要素は、前記第1の規範モデルの属性値と前記第2の規範モデルの属性値の種々の値の組に対する係数値を保持する行列データテーブルを設け、入力される前記両属性値から係数値を発生することを特徴とする2自由度制御装置。The two-degree-of-freedom control device according to claim 5,
The function element is provided with a matrix data table that holds coefficient values for various sets of attribute values of the first reference model and attribute values of the second reference model. A two-degree-of-freedom control device for generating a coefficient value.
前記第1の規範モデルは、外乱を良好に抑制するPIDパラメータを設計するために使用するモデル、あるいはI−PD制御系の特性を同定して得られるモデルであることを特徴とする2自由度制御装置。In the two-degree-of-freedom control device according to claim 5 or 6,
The first reference model is a model used for designing a PID parameter that favorably suppresses a disturbance, or a model obtained by identifying characteristics of an I-PD control system. Control device.
前記規範モデルの属性値は、前記両規範モデルの補間係数、該両規範モデルの最大感度および該両規範モデルの伝達関数係数値のうち、何れか1つを用いることを特徴とする2自由度制御装置。The two-degree-of-freedom control device according to any one of claims 5 to 7,
The attribute value of the reference model uses any one of an interpolation coefficient of the reference models, a maximum sensitivity of the reference models, and a transfer function coefficient value of the reference models. Control device.
前記演算要素は、前記関数要素から発生される係数値に前記第1の規範モデルの立上り時間を乗じ、この乗算値を積分時間で除することにより2自由度化パラメータを発生することを特徴とする2自由度制御装置。The two-degree-of-freedom control device according to claim 5,
The arithmetic element generates a two-degree-of-freedom parameter by multiplying a coefficient value generated from the function element by a rise time of the first reference model, and dividing the multiplied value by an integration time. To control two degrees of freedom.
前記演算要素は、前記関数要素から発生される係数値に前記第1の規範モデルの立上り時間の2乗を乗じ、この乗算値を積分時間と微分時間との乗算値で除算することにより2自由度化パラメータを発生することを特徴とする2自由度制御装置。The two-degree-of-freedom control device according to claim 5,
The arithmetic element multiplies the coefficient value generated from the function element by the square of the rise time of the first reference model, and divides the multiplied value by the multiplied value of the integral time and the derivative time to obtain two free elements. A two-degree-of-freedom control device for generating a degree of freedom parameter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP01602299A JP3571952B2 (en) | 1999-01-25 | 1999-01-25 | 2-DOF control device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP01602299A JP3571952B2 (en) | 1999-01-25 | 1999-01-25 | 2-DOF control device |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2000214902A JP2000214902A (en) | 2000-08-04 |
JP3571952B2 true JP3571952B2 (en) | 2004-09-29 |
Family
ID=11904950
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP01602299A Expired - Lifetime JP3571952B2 (en) | 1999-01-25 | 1999-01-25 | 2-DOF control device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3571952B2 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5136863B2 (en) * | 2010-02-05 | 2013-02-06 | 横河電機株式会社 | Operation monitoring device |
JP6683270B1 (en) * | 2019-02-04 | 2020-04-15 | オムロン株式会社 | Disturbance suppression device, disturbance suppression method, and program |
-
1999
- 1999-01-25 JP JP01602299A patent/JP3571952B2/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2000214902A (en) | 2000-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102033492A (en) | Linear neuron on-line learning adaptive control method and controller for passive system | |
JP4166637B2 (en) | Method for adjusting process control apparatus and adjustment tool therefor | |
JP2772106B2 (en) | 2-DOF adjustment device | |
Tseng et al. | Digital shadow identification from feed drive structures for virtual process planning | |
CN112327604B (en) | Feedforward compensation expected dynamic PI and PID control method | |
JP3571952B2 (en) | 2-DOF control device | |
Zhmud et al. | The two methods of reverse overshoot suppression in automation systems | |
Kersten et al. | Application-based discussion of verified simulations of interval enclosure techniques | |
Jorge et al. | Novel characterization and tuning methods for integrating processes | |
Kurien et al. | Overview of different approaches of pid controller tuning | |
JP2005284828A (en) | Controller, method for tuning control parameter, program, and recording medium | |
López et al. | A tuning procedure for the cascade proportional integral retarded controller | |
JPWO2008018496A1 (en) | Control method and control apparatus | |
JPH09146610A (en) | Multivariable nonlinear process controller | |
JP3641181B2 (en) | PID controller | |
JP4272229B2 (en) | Process control device | |
US6847851B1 (en) | Apparatus for improved general-purpose PID and non-PID controllers | |
CN112947048B (en) | Control method, system and medium of multivariable coupling control system | |
JPH0272404A (en) | Deciding method for membership function | |
Pathiran et al. | A comparative study on the well-known feedback controller design methods | |
Lü et al. | An improved model-free adaptive control with G function fuzzy reasoning regulation design and its applications | |
RU2302028C1 (en) | Method for controlling dynamic objects | |
JP2844289B2 (en) | IMC controller | |
JPS629405A (en) | Process controller | |
López et al. | Two delay-based nonlinear controllers and their tuning methods |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040609 |
|
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: 20040622 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040625 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080702 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090702 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090702 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100702 Year of fee payment: 6 |