JP4155508B2 - Control parameter adjustment method for control device - Google Patents

Control parameter adjustment method for control device Download PDF

Info

Publication number
JP4155508B2
JP4155508B2 JP2003020706A JP2003020706A JP4155508B2 JP 4155508 B2 JP4155508 B2 JP 4155508B2 JP 2003020706 A JP2003020706 A JP 2003020706A JP 2003020706 A JP2003020706 A JP 2003020706A JP 4155508 B2 JP4155508 B2 JP 4155508B2
Authority
JP
Japan
Prior art keywords
control
test
control parameter
output
parameter adjustment
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 - Fee Related
Application number
JP2003020706A
Other languages
Japanese (ja)
Other versions
JP2004234219A (en
Inventor
政志 中本
史郎 日野
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.)
Toshiba Corp
Original Assignee
Toshiba Corp
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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP2003020706A priority Critical patent/JP4155508B2/en
Publication of JP2004234219A publication Critical patent/JP2004234219A/en
Application granted granted Critical
Publication of JP4155508B2 publication Critical patent/JP4155508B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は、制御装置の制御パラメータ調整方法に係わり、特に制御系の時系列の応答データを用いて制御パラメータを調整する制御装置の制御パラメータ調整方法に関する。
【0002】
【従来の技術】
制御装置により制御対象を制御する場合、制御対象の特性に合わせて制御パラメータを調整する必要がある。制御パラメータの調整には、制御対象の応答データから制御パラメータを算出する方法が一般に行われている。
【0003】
以下の非特許文献1には、従来の制御パラメータの調整方法として実用化されている合計6種の方法が掲載されている。
また、以下の特許文献1には、非特許文献1に未掲載のニューラルネットを応用した調整方法が記載されている。
【0004】
ところで、これら2つの先行技術文献に記載されている方法は、原則的にテストを一度行って制御パラメータを決定する方法と、テストを繰り返し行ってパラメータを更新していく方法の2通りの方法に大別できる。テストの繰り返しを行わない方法の一つとして、制御対象の入出力からプラントの数式モデルを求め、この数式モデルからモデルマッチングあるいはジーグラ・ニコルス法で制御パラメータを求める方法が特許文献2のほか前記非特許文献1にも掲載されている。
【0005】
以下、これら従来の制御パラメータ調整方法について、その概要を説明する。図9は従来のパラメータ調整方法の一例を示すもので、モデル化してパラメータを計算する方法を説明するための構成図である。図9に示す制御システムの制御系は、制御対象1を制御する制御装置2と、この制御装置2の制御パラメータを調整するパラメータ計算部10と、制御装置2の出力操作量u(t)にテスト信号w(t)を重畳するテスト信号発生部13と、制御対象1の入力信号v(t)と出力信号y(t)とから制御対象1のモデルを作成するモデル化部9とから構成されている。なお、4−1は目標値r(t)と制御量y(t)とを加算し、偏差e(t)を求める信号加算部、4−2は操作量u(t)にテスト信号w(t)を重畳する信号加算部、また、4−3は制御対象1に雑音n(t)が加わることを表す加算部である。
【0006】
制御対象1のモデル化には制御対象1の伝達関数などのパラメトリックなモデルと、むだ時間、立ち上がり時間、減衰率またはオーバーシュート量などの特徴量をモデル化する方法などがある。そして、モデル化部9でモデル化を行う場合には、制御対象1の入出力信号v(t)、y(t)に多くの情報量を含める必要があり、このため一般には入力信号v(t)に重畳するテスト信号w(t)を変化させるようにしている。
【0007】
パラメータ計算部10は、モデル化部9からモデル化した制御対象1の情報(伝達関数のパラメータ、特徴量など)を入力して制御パラメータを計算により求める。この制御パラメータ計算部10における制御パラメータの計算には、一般にモデルマッチングあるいはジーグラ・ニコルス法が用いられるが、ファジー推論を用いる方法もある(前記非特許文献1の129頁〜132頁、図6.13〜図6.16を参照)。
【0008】
テストを繰り返し行って制御パラメータを更新する方法としては、前記特許文献1に記載のニューラルネットを応用した方法がある。以下、このニューラルネットを応用した制御パラメータ更新方法の概要について図10を参照して説明する。
【0009】
図10はこの制御パラメータ更新方法を説明するための構成図であり、この制御システムの制御系は、制御対象1を制御する制御装置2と、学習用データ部11と、調整ルール学習システム14と、制御装置2の制御パラメータを算出するパラメータ算出部12と、前記制御装置2から出力された操作量u(t)にテスト信号w(t)を重畳して前記制御対象1に入力するテスト信号発生部13とを備えている。
【0010】
そしてこの図10による制御パラメータ更新方法は、操作量u(t)に重畳するテスト信号w(t)をあるいは目標信号r(t)を変更してテストを行い、このテストの繰り返しによって制御対象1の入出力信号v(t)、y(t)および制御偏差e(t)などの情報をパラメータ算出部12に入力し、調整ルール学習システム14から出力される調整ルールに従って制御パラメータの更新を行うようにしたものである。
【0011】
この制御パラメータの更新ルールはニューラルネットを利用した調整ルール学習システム14から与えられ、この調整ルール学習システム14は学習用データ11によりルールが決められる。
【0012】
【非特許文献1】
システム制御情報学会編、須田信英著作代表「PID制御」朝倉書店出版、1992年7月20日、p.107−117、118−172
【特許文献1】
特許第2862308号公報(第1−18頁、第1図−第28図)
【特許文献2】
特許第2755644号公報(第1−8頁、第1図−第18図(b))
【0013】
【発明が解決しようとする課題】
一般に、雑音n(t)や制御対象の非線型性などによるモデル化には誤差が生じる。同様に特徴量を抽出する場合も誤差が生じる。このため、テストの繰り返しを前提としない方法では、テストにより得られた情報から、必ずしも正しい制御パラメータが決定されるか否かは明確ではない。
さらにファジー推論を用いる方法は、パラメータを算出するルールが恣意的で客観性がないため、ルールの構築に時間と労力を必要とする。
【0014】
一方、テストを繰り返して制御パラメータを更新するニューラルネットを用いた従来の方法では、調整ルールを定めるルールの学習に多大な時間と多くのテストデータを必要とする。
また、いずれの制御パラメータ調整方法も、選択された制御パラメータが妥当であるかどうかを定量的に評価・判断する方法がない。
【0015】
さらに、いずれの方法も制御対象の情報を明確に得るためには、制御装置を制御対象から切り離した、いわゆる開ループでテストを行うことが望ましいが、その場合には不安定な制御対象には適用が不可能である。
【0016】
そこで、本発明の目的は、短時間で制御パラメータを決定でき、かつ選ばれた制御パラメータの定量的な評価を同時に行うことのできる制御装置の制御パラメータ調整方法を提供することにある。
【0017】
【課題を解決するための手段】
上記の目的を達成するために、請求項1に係わる制御装置の制御パラメータ調整方法の発明は、制御量を目標値に一致させるように制御対象を制御する制御装置の制御パラメータ調整方法において、制御系を制御対象と制御装置とを組み合わせた閉ループの状態で1回目のテスト信号を加え、そのときの応答データを用いて2回目のテスト信号を作り、この2回目のテスト信号を加えたときの応答データを用いて制御パラメータの変化量を決め、これら2回一組としたテストを繰り返し、設定した評価関数を最小とする制御パラメータを決定することを特徴とする。
【0018】
【発明の実施の形態】
(第1の実施の形態)
図1は本発明による制御パラメータ調整方法を適用した制御システムの1実施の形態を示す図である。
【0019】
本実施の形態による制御システムの制御系は、制御対象1を制御する制御装置2と、この制御装置2の制御パラメータを調整するとともに、テスト信号 w(t) を発生させる制御パラメータ調整システム3と、前記制御装置2から出力された操作量u(t)を偏微分して、操作量u(t)の、操作パラメータxについての変化率を求め、この変化率を前記制御パラメータ調整システム3に入力する信号変換部5と、制御対象1から出力される制御量y(t)を偏微分して、制御量y(t)の、操作パラメータxについての変化率を求め、この変化率を前記制御パラメータ調整システム3に出力する信号変換部7とから構成されている。
【0020】
次に、制御パラメータ調整システム3による制御パラメータの調整方法について説明する。なお、事前の条件として、制御対象1は制御装置2により必ずしも最適ではないが安定に制御されているものとする。
【0021】
まず、制御性能を測る制御指標を設定する。この制御指標としては目標値r(t)と出力y(t)との偏差e(t)をはじめとして操作量u(t)、操作量の変化量u(t)−u(t-1)などを使って1回目のテストからN回目のテストまでのテスト期間の総和を考える。具体的な制御指標として、以下の式(1)から(3)などがある。
【0022】
【数1】

Figure 0004155508
【0023】
上記の(1)式〜(3)式において、tはサンプル、λ(t)はサンプルtの重み係数、uは操作量の定常値とし、総和はテスト期間t=1からt=Nまでの加算平均とする。
【0024】
これら(1)式〜(3)式に示すいずれかの評価指標Jを選んでテストを行い、このテストの時系列応答データ(以下、単に応答データという)から決まる評価指標の値を最小にする制御パラメータxを求めることで、制御装置2のパラメータ調整を行う。
【0025】
次に、評価指標Jを最小にする調整方法について上記(1)式を例にして説明する。評価指標Jは制御装置2のパラメータxの関数であり、制御パラメータを変数xと書くと、評価指標Jは制御パラメータの関数J(x)として書くことができる。
【0026】
このパラメータの関数J(x)が制御パラメータxについて最小になる条件は∂J(x)/∂x=0となる場合である。(1)式を制御パラメータxで偏微分すると、(4)式に示す関数の感度すなわち勾配が求まる。そして、感度(勾配)を最小すなわち∂J(x)/∂x=0とする制御パラメータxをニュートン法により、逐次更新すると(5)式に示す制御パラメータの値が求まる。
【0027】
【数2】
Figure 0004155508
【0028】
なお、(5)式において、xkはk回目の更新により求められた制御パラメ−タの値とする。Rkは正定値行列であり、添字kはk回目の更新時の値を示す。
上記(5)式による制御パラメータの更新には(4)式の値が必要であり、(4)式中の偏差e(t)、操作量u(t)はテストを行って応答データから直接求める。一方、(4)式中の制御量y(t)、操作量u(t)の制御パラメータxについての変化率∂y(t)/∂x、∂u(t)/∂xはデータから直接得ることはできないが、2回のテストを行って応答データから求める。
【0029】
図1における制御対象1と制御装置2とからなる閉ループ系について、目標値r(t)と出力である制御量y(t)との関係、目標値r(t)と操作量u(t)との関係はそれぞれ次のようになる。
【0030】
【数3】
Figure 0004155508
【0031】
上式から、∂y(t)/∂x、∂u(t)/∂xはそれぞれ、以下の(6)式、(7)式から求めることができる。
【0032】
【数4】
Figure 0004155508
なお、ここで(6)式および(7)式の∂C(x)/∂xは制御装置2の制御演算式から計算することができる。また、1/1+PC(x)・e(t)およびP/1+PC(x)・e(t)は、それぞれ1回目のテストの結果得られた偏差e(t)すなわちr(t)−y(t)を、2回目のテストの際、テスト信号w(t)とした場合の操作量u(t)および出力y(t)である。
【0033】
このように、本実施の形態は2回のテストを一組としてテストを行うことにより(4)式の計算に必要なデータを応答データとして求め、さらにこの応答データを(5)式に従って計算することにより、制御パラメータを逐次更新することができる。
【0034】
次に、以上説明した制御パラメータの更新の考え方に基づいて、図1における制御パラメータ調整システム3の動作を説明する。
この方法は(4)式の中の制御量y(t)、操作量u(t)の、制御パラメータxについての変化率∂y(t)/∂x、∂u(t)/∂xを求めるために、(6)式、(7)式を用いる方法である。
【0035】
▲1▼.1回目のテスト:
図3の動作説明図で示すように、制御パラメータ調整システム3は、発生させたテスト信号w(t)を目標値r(t)として用い、目標値そのものをr(t)からr1(t)に変化させる。そしてこのときの偏差e(t)および操作量u(t)を計測し、応答データe1(t)およびu1(t)を前記制御パラメータ調整システム3に入力する。この制御パラメータ調整システム3ではこの入力したデータe1(t)、u1(t)を1回目のテストの応答データとして保存する。
【0036】
▲2▼.2回目のテスト:
図4の動作説明図で示すように、保存しておいた1回目のテストの応答データである偏差e1(t)をテスト信号w2(t)とし、これを信号加算部4−2において操作量u2(t)に重畳させてv2(t)を制御対象1に加えることにより、2回目のテストを行い、操作量u2(t)、出力y2(t)を測定する。そして操作量u2(t)、出力y2(t)についてそれぞれ信号変換部5、7で偏微分により∂C(x)/∂xを求め、それぞれ∂y(t)/∂x、∂u(t)/∂xを生成し、制御パラメータ調整システム3に入力する。
【0037】
▲3▼.∂J(x)/∂xの計算:
制御パラメータ調整システム3は、取込んだ∂y(t)/∂x、∂u(t)/∂xを(4)式に代入して∂J(x)/∂xを計算する。(4)式で用いるe(t)、u(t)は1回目のテストの応答データe1(t)、u1(t)を用い、∂y(t)/∂x、∂u(t)/∂xは2回目のテストから得た応答データから生成した値を用いる。
【0038】
▲4▼.パラメータの更新:
制御パラメータ調整システム3は、(5)式に従って制御パラメータxを逐次更新することにより、xk+1を得る。λは適当な進み幅である。
【0039】
なお、(5)式において、Rkは単位行列でもよいし、∂y(t)/∂xおよび∂u(t)/∂xを使ってガウス・ニュートン法による次式(8)から求めるようにしてもよい。
【0040】
【数5】
Figure 0004155508
【0041】
以上述べたように、第1の実施の形態によれば、制御指標を定量的に表す評価関数を最小にする制御パラメータを求めることができる。また、制御パラメータの更新に評価関数の制御パラメータについての感度(勾配)を用いるため、最適値への収束を早くすることができる。さらに、評価関数の制御パラメータについての感度(勾配)を求めるため、制御パラメータの収束状況を定量的に把握することができる。
また、テストが閉ループで実施されるため、不安定な制御対象にも適用できる。
【0042】
さらに、制御量y(t)、操作量u(t)の、制御パラメータxについての変化率∂y(t)/∂xおよび∂u(t)/∂xを計算する際、それぞれ(6)式、(7)式を用いているため,計算が簡単になるという特長を有する。
【0043】
(第2の実施の形態)
第2の実施の形態について図2に示す制御システム構成図を参照して説明する。本実施の形態と、第1の実施の形態との主な相違点は、(4)式の中の∂y(t)/∂x、∂u(t)/∂xを求めるために、本実施の形態では第1の実施の形態における信号変換部5および7をそれぞれ信号変換部6および8に替えて(9)、(10)式を演算するようにした点である。
【0044】
すなわち、図2における制御対象1と制御装置2とからなる閉ループ系について、目標値と出力である制御量y(t)との関係、目標値と操作量u(t)との関係はそれぞれ次のようになる。
【0045】
【数6】
Figure 0004155508
【0046】
上式から、∂y(t)/∂xおよび∂u(t)/∂xはそれぞれ、(9)式および(10)式として求めることができる。
【0047】
【数7】
Figure 0004155508
【0048】
なお、ここで、(9)および(10)式の1/C(x)・∂C(x)/∂xは制御装置2の制御演算式から計算することができる。また、C(x)/1+PC(x)・e(t)およびPC(x)/1+PC(x)・e(t)は、それぞれ1回目のテストで得られた偏差e(t)すなわちr(t)-y(t)を、2回目のテストの際、目標値r(t)とした場合の操作量u(t)および出力y(t)である。
【0049】
このように、第2の実施の形態の場合も第1の実施の形態と同様、2回のテストを1組としてテストを行うことにより(4)式の計算に必要なデータを応答データとして求め、さらにこの応答データを(5)式に従って計算することにより、制御パラメータを逐次更新することができる。
以上説明した制御パラメータの更新の考え方に基づいて、図2における制御パラメータ調整システム3の動作を説明する。
【0050】
▲1▼1回目のテスト:
図3の動作説明図において、1回目のテストを行う。このテスト方法は第1の実施の形態による方法と同一である。
【0051】
▲2▼.2回目のテスト:
図5の動作説明図で示すように、制御パラメータ調整システム3は保存しておいた1回目のテストの応答データである偏差e1(t)をテスト信号w2(t)として出力する。そしてこのテスト信号を目標値r(t)にして2回目のテストを行い、操作量u2(t)、出力y2(t)を測定する。さらに操作量u2(t)、出力y2(t)についてそれぞれ信号変換部6、8で偏微分によりそれぞれ1/C(x)・∂C(x)/∂xを求め、それぞれ∂y(t)/∂xおよび∂u(t)/∂xを生成し、制御パラメータ調整システム3に入力する。
【0052】
▲3▼.∂J(x)/∂xの計算:
制御パラメータ調整システム3は、取込んだ∂y(t)/∂x、∂u(t)/∂xを(4)式に代入して∂J(x)/∂xを計算する。(4)式で用いるe(t)、u(t)は1回目のテスト結果e1(t)、u1(t)を用い、∂y(t)/∂xおよび∂u(t)/∂xは2回目のテストから得た応答データから生成した値を用いる。
【0053】
▲4▼.パラメータの更新:
制御パラメータ調整システム3は、(5)式に従って制御パラメータxを逐次更新することにより、xk+1を得る。λは適当な進み幅である。
【0054】
なお、(5)式において、Rkは単位行列でもよいし、∂y(t)/∂xおよび∂u(t)/∂xを使ってガウス・ニュートン法による前記(8)式から求めるようにしてもよい。
【0055】
以上述べたように、第2の実施の形態によれば(9)式、(10)式を計算するので、(6)式、(7)式を計算する第1の実施の形態と比較して計算自体は複雑になるが,2回のテストが共に目標値信号の変更で実施できる。このためコントローラと制御パラメータ調整システム3との信号のインターフェイスが簡単になる。特に、多くの汎用のコントローラは標準的に目標値信号を外部から入力できるようになっているため、コントローラのインターフェイスや、制御プログラムの変更が不要である。このため、現場に実装済みの制御装置2にハードウェアやソフトウェアの変更なしで、そのまま適用することができる。
【0056】
(第3の実施の形態)
本実施の形態は、第1および第2の実施の形態によるパラメータの調整方法において、評価指標を(2)式とした場合である。評価指標を(2)式とする場合は、以上説明した操作量の評価をu(t)からu(t)-u(t-1)に変更すればよい。
【0057】
なお、制御対象が定位系すなわち制御量の定常値を変化させるために操作量も定常的に変化するような系の場合、テスト信号で目標値をステップ状に変更させると、操作量u(t)は時間が経過してもゼロにならない。偏差e(t)はゼロに向かうため、前述した(1)式の評価関数では、重み係数λ(t)がゼロ以外では、操作量u(t)が大きく評価され、結果的に操作量を小さくする制御パラメータが選択される。このため定位系では望ましい応答が得られない場合がある。
【0058】
しかしながら、この第3の実施の形態では、前記(2)式の評価関数を採用したので、制御対象の特性にかかわらず、目標値のステップ状の変化に対して時間の経過にしたがい偏差e(t)と操作量の差分u(t)-u(t-1)は共にゼロに向かうため、偏差と操作量の評価が重み係数λ(t)で調整でき、定位系における前記の不具合を解消することができる。
【0059】
(第4の実施の形態)
本実施の形態は、第1および第2の実施の形態で述べたパラメータの調整方法において、評価指標を(3)式とした場合である。評価指標を(3)式とする場合は、以上説明した操作量の評価をu(t)からu(t)-uに変更すればよい。
【0060】
第4の実施の形態によれば、評価関数の中の、偏差e(t)と操作量の変化分u(t)-uは時間の経過に従ってゼロに向かうため、前述した第3の実施の形態と同様の効果を得ることができる。
【0061】
(第5の実施の形態)
次に、第5の実施の形態による調整方法を図2に示す。この実施の形態と第2の実施の形態との相違点は、第2の実施の形態の場合、制御パラメータ調整システム3が1回目のテストで目標値r(t)を変化させるのに対して、本実施の形態の場合は1回目のテストで操作量の出口のテスト信号w(t)を変化させるようにした点にあり、2回目のテストの結果は図4と同じである。原理は第2の実施の形態の場合と同様である。
【0062】
この第5の実施の形態は、(4)式の中の∂y(t)/∂xおよび∂u(t)/∂xを求めるために(6)、(7)式を用いる方法では制御パラメータの更新は次のようになる。
【0063】
▲1▼.1回目のテスト:
図6の動作説明図において、制御パラメータ調整システム3からテスト信号w(t)を発生させ、目標値を変化させる。このときの偏差e(t)および操作量u(t)を計測して、e1(t)、u1(t)として制御パラメータ調整システム3に保存する。
【0064】
▲2▼.2回目のテスト:
図7の操作説明図において、1回目のテストの結果保存されている偏差e1(t)をテスト信号w(t)として用いて2回目のテストを行い、操作量u(t)、出力y(t)を測定する。そしてこれらの測定量を信号変換部5および7で偏微分してそれぞれ∂Cx/∂xの計算を行い、生成した∂y(t)/∂xおよび∂u(t)/∂xを制御パラメータ調整システム3に入力する。
【0065】
▲3▼.∂J(x)/∂xの計算:
制御パラメータ調整システム3は、(4)式に従って、∂J(x)/∂xを計算する。(4)式におけるe(t)、u(t)は1回目のテスト結果e1(t)、u1(t)を用い、また、∂y(t)/∂x、∂u(t)/∂xは2回目のテストから得た応答データを使って生成した値を用いる。
【0066】
▲4▼.パラメータの更新:
制御パラメータ調整システム3は、(4)式で求めた値を(5)式に代入し、制御パラメータxを求め、逐次更新する。(5)式のλは適当な進み幅である。なお、Rkは単位行列でもよいし、∂y(t)/∂x、∂u(t)/∂xを使ってガウス・ニュートン法による前記(8)式を用いてもよい。
【0067】
以上述べた第5の実施の形態では、計算について第1の実施の形態と同様の効果のほか、さらに次のような効果を得ることができる。すなわち、第1の実施の形態から第4の実施の形態の場合、1回目のテストで目標値を変更するため、これにしたがって制御量が目標値に一致するように変化する。このため制御対象の運転状態が変化してしまう。
【0068】
しかしながら、この第5の実施の形態では、操作量にテスト信号を与えて一時的に制御量を変化させてはいるが、目標値が一定であるため制御量は元の値に戻る。このため、制御対象の運転状態に与える影響を小さくできるといった効果を得ることができる。
【0069】
(第6の実施の形態)
第6の実施の形態について図2に示す制御システム構成図を参照して説明する。本実施の形態は、(4)式の中の∂y(t)/∂x、∂u(t)/∂xを求めるために(9)(10)式を用いて制御パラメータを調整するようにした方法に関するものである。
【0070】
▲1▼.1回目のテスト:
第1回目のテストは、第5の実施の形態と同一であり、また、動作説明図は図6のとおりである。
【0071】
▲2▼.2回目のテスト:
図8の動作説明図で示すように、目標値r(t)を1回目のテストの目標値と出力との偏差e1(t)としてテストを行い、操作量u(t)、出力y(t)を測定し、信号変換部6および8でそれぞれ1/C(x)・∂C(x)/∂xの計算を行い、∂y(t)/∂xおよび∂u(t)/∂xを生成し、制御パラメータ調整システム3に入力し保存する。
【0072】
▲3▼.∂J(x)/∂xの計算:
制御パラメータ調整システム3は、取込んだ∂y(t)/∂x、∂u(t)/∂xを(4)式に代入して∂J(x)/∂xを計算する。(4)式で用いるe(t)、u(t)は1回目のテストの応答データe1(t)、u1(t)を用い、∂y(t)/∂xおよび∂u(t)/∂xは2回目のテストから得た応答データから生成した値を用いる。
【0073】
▲4▼.パラメータの更新:
制御パラメータ調整システム3は、(5)式に従って制御パラメータxを逐次更新することにより、xk+1を得る。λは適当な進み幅である。
なお、(5)式において、Rkは単位行列でもよいし、∂y(t)/∂xおよび∂u(t)/∂xを使ってガウス・ニュートン法による前記(8)式から求めるようにしてもよい。
【0074】
以上述べたように、第6の実施の形態の場合、2回目のテストは目標値を変化するが、1回目のテストで制御量がゼロに向かうため、2回目のテストの目標値もゼロに向かう信号となる。このため、運転状態に与える影響が小さくなり、第5の実施の形態と同様の効果を得ることができる。
【0075】
(第7の実施の形態)
本実施の形態は、以上述べたパラメータの調整方法において、評価指標を(2)式とした場合である。評価指標を(2)式とする場合は、以上説明した操作量の評価をu(t)からu(t)-u(t-1)に変更すればよい。
この第7の実施の形態によれば、第3の実施の形態と第5の実施の形態あるいは第6の実施の形態を合わせた効果を得ることができる。
【0076】
(第8の実施の形態)
本実施の形態によるパラメータの調整方法は、評価指標を(3)式とする場合のパラメータの調整方法に関するものである。この場合、操作量の評価をu(t)からu(t)-uに変更すればよい。
この第8の実施の形態によれば、第4、第5あるいは第6の実施の形態を合わせた効果が得られる。
【0077】
【発明の効果】
以上述べたように本発明によれば、制御指標を定量的に表す評価関数を最小にする制御パラメータを求めることができる。
また、制御パラメータの更新に評価関数の制御パラメータについての感度(勾配)を用いるため、最適値への収束を早くすることができる。
【0078】
さらに、評価関数の制御パラメータについての感度(勾配)が求められるため、制御パラメータの収束状況が定量的に把握することができる。
またさらに、テストを閉ループで実施するため、不安定な制御対象についても適用することができる。
【図面の簡単な説明】
【図1】本発明を説明するためのシステム構成例の一例を示す図。
【図2】本発明を説明するためのシステム構成例の他の例を示す図。
【図3】発明の第1の実施の形態におけるパラメータ調整の第1回目のテスト時の作用を示す図。
【図4】発明の第1の実施の形態におけるパラメータ調整の第2回目のテスト時の作用を示す図。
【図5】本発明の第2の実施の形態におけるパラメータ調整の第2回目のテスト時の作用を示す図。
【図6】本発明の第5および第6の実施の形態におけるパラメータ調整の第1回目のテスト時の作用を示す図。
【図7】本発明の第5の実施の形態におけるパラメータ調整の第2回目のテスト時の作用を示す図。
【図8】本発明の第6の実施の形態におけるパラメータ調整の第2回目のテスト時の作用を示す図。
【図9】従来の制御パラメータ更新方法の一例を説明するための構成図。
【図10】従来の制御パラメータ更新方法の他の例を説明するための構成図。
【符号の説明】
1…制御対象、2…制御装置、3…制御パラメータ調整システム、4−1〜4−3…信号加算部、5〜8…信号変換部。[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a control parameter adjustment method for a control device, and more particularly to a control parameter adjustment method for a control device that adjusts control parameters using time-series response data of a control system.
[0002]
[Prior art]
When a control target is controlled by a control device, it is necessary to adjust control parameters in accordance with the characteristics of the control target. For adjusting the control parameter, a method of calculating the control parameter from response data to be controlled is generally performed.
[0003]
Non-Patent Document 1 below lists a total of six methods that have been put to practical use as conventional control parameter adjustment methods.
Also, the following Patent Document 1 describes an adjustment method using a neural network that is not yet published in Non-Patent Document 1.
[0004]
By the way, the methods described in these two prior art documents are basically divided into two methods: a method in which a test is performed once to determine a control parameter, and a method in which a test is repeatedly performed to update a parameter. It can be divided roughly. As one of the methods that do not repeat the test, there is a method of obtaining a plant mathematical model from input / output of a controlled object, and obtaining a control parameter from this mathematical model by model matching or Ziegler-Nichols method in addition to the above-mentioned non-patent documents. It is also published in Patent Document 1.
[0005]
The outline of these conventional control parameter adjustment methods will be described below. FIG. 9 shows an example of a conventional parameter adjustment method, and is a configuration diagram for explaining a method of calculating parameters by modeling. The control system of the control system shown in FIG. 9 includes a control device 2 that controls the controlled object 1, a parameter calculation unit 10 that adjusts the control parameter of the control device 2, and an output operation amount u (t) of the control device 2. The test signal generator 13 for superimposing the test signal w (t), and the modeling unit 9 for creating a model of the control target 1 from the input signal v (t) and the output signal y (t) of the control target 1 Has been. Note that 4-1 is a signal adder that adds the target value r (t) and the control amount y (t) to obtain the deviation e (t), and 4-2 is the test signal w ( A signal addition unit that superimposes t), and 4-3 is an addition unit that indicates that noise n (t) is added to the control target 1.
[0006]
The modeling of the controlled object 1 includes a parametric model such as a transfer function of the controlled object 1 and a method of modeling a feature quantity such as a dead time, a rise time, an attenuation rate, or an overshoot amount. When modeling is performed by the modeling unit 9, it is necessary to include a large amount of information in the input / output signals v (t) and y (t) of the controlled object 1. The test signal w (t) superimposed on t) is changed.
[0007]
The parameter calculation unit 10 receives information (transfer function parameters, feature amounts, etc.) of the controlled object 1 modeled from the modeling unit 9 and obtains control parameters by calculation. In general, model matching or Ziegler-Nichols method is used for calculation of the control parameter in the control parameter calculation unit 10, but there is a method using fuzzy inference (see Non-Patent Document 1, pages 129 to 132, FIG. 6). 13 to FIG. 6.16).
[0008]
As a method for updating a control parameter by repeatedly performing a test, there is a method using a neural network described in Patent Document 1. Hereinafter, an outline of a control parameter update method using this neural network will be described with reference to FIG.
[0009]
FIG. 10 is a block diagram for explaining this control parameter update method. The control system of this control system includes a control device 2 that controls the control object 1, a learning data unit 11, an adjustment rule learning system 14, and the like. A parameter calculation unit 12 that calculates a control parameter of the control device 2 and a test signal that is input to the control target 1 by superimposing a test signal w (t) on the operation amount u (t) output from the control device 2 And a generator 13.
[0010]
In the control parameter updating method according to FIG. 10, the test is performed by changing the test signal w (t) superimposed on the operation amount u (t) or the target signal r (t), and the control object 1 is obtained by repeating this test. Information such as input / output signals v (t), y (t) and control deviation e (t) is input to the parameter calculator 12 and the control parameters are updated in accordance with the adjustment rules output from the adjustment rule learning system 14. It is what I did.
[0011]
The update rule for the control parameter is given from the adjustment rule learning system 14 using a neural network, and the rule is determined by the learning data 11 in the adjustment rule learning system 14.
[0012]
[Non-Patent Document 1]
Edited by System Control Information Society, Nobuhide Suda's representative “PID Control”, Asakura Shoten Publishing, July 20, 1992, p. 107-117, 118-172
[Patent Document 1]
Japanese Patent No. 2862308 (page 1-18, FIGS. 1 to 28)
[Patent Document 2]
Japanese Patent No. 2755644 (pages 1-8, FIGS. 1-18 (b))
[0013]
[Problems to be solved by the invention]
In general, an error occurs in modeling due to noise n (t), nonlinearity of a control target, and the like. Similarly, an error occurs when extracting a feature amount. For this reason, it is not clear whether the correct control parameter is determined from the information obtained by the test in the method that does not assume the repetition of the test.
Furthermore, the method using fuzzy inference requires time and labor to construct rules because the rules for calculating parameters are arbitrary and not objective.
[0014]
On the other hand, in a conventional method using a neural network in which a test parameter is updated by repeating a test, it takes a lot of time and a lot of test data to learn a rule for determining an adjustment rule.
In addition, none of the control parameter adjustment methods has a method for quantitatively evaluating and judging whether or not the selected control parameter is valid.
[0015]
Furthermore, in order to clearly obtain information about the controlled object in either method, it is desirable to perform a so-called open loop test in which the control device is separated from the controlled object. It is impossible to apply.
[0016]
Accordingly, an object of the present invention is to provide a control parameter adjustment method for a control device that can determine a control parameter in a short time and can simultaneously perform quantitative evaluation of selected control parameters.
[0017]
[Means for Solving the Problems]
In order to achieve the above object, a control parameter adjustment method for a control device according to a first aspect of the present invention provides a control parameter adjustment method for a control device that controls a control target so that a control amount matches a target value. When the first test signal is added in a closed loop state in which the control target and the control device are combined in the system, the second test signal is created using the response data at that time, and this second test signal is added. It is characterized in that the amount of change in the control parameter is determined using the response data, a test that is a set of two times is repeated, and the control parameter that minimizes the set evaluation function is determined.
[0018]
DETAILED DESCRIPTION OF THE INVENTION
(First embodiment)
FIG. 1 is a diagram showing an embodiment of a control system to which a control parameter adjusting method according to the present invention is applied.
[0019]
The control system of the control system according to the present embodiment includes a control device 2 that controls the controlled object 1, a control parameter adjustment system 3 that adjusts the control parameters of the control device 2 and generates a test signal w (t). The operation amount u (t) output from the control device 2 is partially differentiated to obtain a change rate of the operation amount u (t) with respect to the operation parameter x, and this change rate is given to the control parameter adjustment system 3. The input signal converter 5 and the controlled variable y (t) output from the controlled object 1 are partially differentiated to determine the rate of change of the controlled variable y (t) with respect to the operation parameter x. And a signal conversion unit 7 for outputting to the control parameter adjustment system 3.
[0020]
Next, a control parameter adjustment method by the control parameter adjustment system 3 will be described. As a precondition, it is assumed that the control object 1 is not necessarily optimal but is stably controlled by the control device 2.
[0021]
First, a control index for measuring control performance is set. This control index includes the deviation e (t) between the target value r (t) and the output y (t), the manipulated variable u (t), and the manipulated variable change u (t) -u (t-1) The total test period from the first test to the Nth test is considered. Specific control indexes include the following formulas (1) to (3).
[0022]
[Expression 1]
Figure 0004155508
[0023]
In the above equations (1) to (3), t is a sample, λ (t) is a weighting factor of sample t, u is a steady value of the manipulated variable, and the total is from test period t = 1 to t = N Addition average of.
[0024]
A test is performed by selecting one of the evaluation indexes J shown in the equations (1) to (3), and the value of the evaluation index determined from the time series response data (hereinafter simply referred to as response data) of this test is minimized. By determining the control parameter x, the parameter of the control device 2 is adjusted.
[0025]
Next, an adjustment method for minimizing the evaluation index J will be described using the above equation (1) as an example. The evaluation index J is a function of the parameter x of the control device 2. When the control parameter is written as a variable x, the evaluation index J can be written as a function J (x) of the control parameter.
[0026]
The condition that the function J (x) of this parameter is minimized with respect to the control parameter x is when ∂J (x) / ∂x = 0. When the equation (1) is partially differentiated by the control parameter x, the sensitivity, that is, the gradient of the function shown in the equation (4) is obtained. Then, when the control parameter x that minimizes the sensitivity (gradient), that is, ∂J (x) / ∂x = 0, is sequentially updated by the Newton method, the value of the control parameter shown in equation (5) is obtained.
[0027]
[Expression 2]
Figure 0004155508
[0028]
In equation (5), xk is the value of the control parameter obtained by the k-th update. Rk is a positive definite matrix, and the subscript k indicates the value at the k-th update.
The value of equation (4) is required to update the control parameter using equation (5) above, and the deviation e (t) and manipulated variable u (t) in equation (4) are tested directly from the response data. Ask. On the other hand, the change rates ∂y (t) / ∂x and ∂u (t) / ∂x for the control parameter x of the control amount y (t) and the manipulated variable u (t) in the equation (4) are directly calculated from the data. Although it cannot be obtained, it is obtained from response data after two tests.
[0029]
In the closed loop system composed of the control object 1 and the control device 2 in FIG. 1, the relationship between the target value r (t) and the output control amount y (t), the target value r (t) and the manipulated variable u (t). The relationship is as follows.
[0030]
[Equation 3]
Figure 0004155508
[0031]
From the above equation, ∂y (t) / ∂x and ∂u (t) / ∂x can be obtained from the following equations (6) and (7), respectively.
[0032]
[Expression 4]
Figure 0004155508
Here, ∂C (x) / ∂x in the equations (6) and (7) can be calculated from the control arithmetic expression of the control device 2. Also, 1/1 + PC (x) · e (t) and P / 1 + PC (x) · e (t) are the deviation e (t) obtained as a result of the first test, that is, r (t ) −y (t) are the operation amount u (t) and the output y (t) when the test signal w (t) is used in the second test.
[0033]
As described above, in the present embodiment, the data required for the calculation of the equation (4) is obtained as response data by performing the test with two tests as a set, and the response data is further calculated according to the equation (5). Thus, the control parameters can be updated sequentially.
[0034]
Next, the operation of the control parameter adjustment system 3 in FIG. 1 will be described based on the concept of updating the control parameters described above.
In this method, the change rates ∂y (t) / ∂x and ∂u (t) / ∂x with respect to the control parameter x of the control amount y (t) and the manipulated variable u (t) in the equation (4) are calculated. In order to obtain this, it is a method using the equations (6) and (7).
[0035]
(1). First test:
As shown in the operation explanatory diagram of FIG. 3, the control parameter adjustment system 3 uses the generated test signal w (t) as the target value r (t) and converts the target value itself from r (t) to r 1 (t ). The deviation e (t) and the manipulated variable u (t) at this time are measured, and response data e 1 (t) and u 1 (t) are input to the control parameter adjustment system 3. The control parameter adjustment system 3 stores the input data e 1 (t) and u 1 (t) as response data for the first test.
[0036]
(2). Second test:
As shown in the operation explanatory diagram of FIG. 4, the deviation e 1 (t), which is the response data of the first test stored, is used as the test signal w 2 (t), and this is added to the signal adder 4-2. by superimposed on the manipulated variable u 2 (t) is added v 2 a (t) to the control object 1, it performs a second test, the operation amount u 2 (t), measuring an output y 2 (t). Then, for the manipulated variable u 2 (t) and output y 2 (t), 変 換 C (x) / ∂x is obtained by partial differentiation in the signal converters 5 and 7, respectively, and ∂y (t) / ∂x and ∂u, respectively. (t) / ∂x is generated and input to the control parameter adjustment system 3.
[0037]
(3). Calculation of ∂J (x) / ∂x:
The control parameter adjustment system 3 calculates ∂J (x) / ∂x by substituting the acquired ∂y (t) / ∂x and ∂u (t) / ∂x into the equation (4). The e (t) and u (t) used in the equation (4) are the response data e 1 (t) and u 1 (t) of the first test, and ∂y (t) / ∂x and ∂u (t ) / ∂x uses a value generated from response data obtained from the second test.
[0038]
(4). Update parameters:
The control parameter adjustment system 3 obtains x k + 1 by sequentially updating the control parameter x according to the equation (5). λ is an appropriate advance width.
[0039]
In Eq. (5), Rk may be a unit matrix, or ∂y (t) / ∂x and ∂u (t) / ∂x may be used to obtain from the following equation (8) using the Gauss-Newton method. May be.
[0040]
[Equation 5]
Figure 0004155508
[0041]
As described above, according to the first embodiment, the control parameter that minimizes the evaluation function that quantitatively represents the control index can be obtained. Moreover, since the sensitivity (gradient) about the control parameter of the evaluation function is used to update the control parameter, the convergence to the optimum value can be accelerated. Furthermore, since the sensitivity (gradient) for the control parameter of the evaluation function is obtained, the convergence state of the control parameter can be quantitatively grasped.
Further, since the test is performed in a closed loop, it can be applied to an unstable control target.
[0042]
Further, when calculating the change rates ∂y (t) / ∂x and ∂u (t) / ∂x for the control parameter x of the control amount y (t) and the manipulated variable u (t), respectively (6) Since the equations (7) and (7) are used, the calculation is easy.
[0043]
(Second Embodiment)
A second embodiment will be described with reference to the control system configuration diagram shown in FIG. The main difference between this embodiment and the first embodiment is that in order to obtain ∂y (t) / ∂x and ∂u (t) / ∂x in the equation (4), In the embodiment, the signal converters 5 and 7 in the first embodiment are replaced with the signal converters 6 and 8, respectively, and the equations (9) and (10) are calculated.
[0044]
That is, the relationship between the target value and the controlled variable y (t) and the relationship between the target value and the manipulated variable u (t) are as follows for the closed loop system composed of the control object 1 and the control device 2 in FIG. become that way.
[0045]
[Formula 6]
Figure 0004155508
[0046]
From the above equation, ∂y (t) / ∂x and ∂u (t) / ∂x can be obtained as equations (9) and (10), respectively.
[0047]
[Expression 7]
Figure 0004155508
[0048]
Here, 1 / C (x) · ∂C (x) / ∂x in the equations (9) and (10) can be calculated from the control arithmetic expression of the control device 2. In addition, C (x) / 1 + PC (x) · e (t) and PC (x) / 1 + PC (x) · e (t) are the deviations e (t ), That is, r (t) -y (t) is the manipulated variable u (t) and output y (t) when the target value r (t) is used in the second test.
[0049]
Thus, in the case of the second embodiment as well, as in the first embodiment, the data required for the calculation of the equation (4) is obtained as response data by performing the test with two tests as one set. Further, by calculating this response data according to the equation (5), the control parameters can be updated sequentially.
Based on the control parameter update concept described above, the operation of the control parameter adjustment system 3 in FIG. 2 will be described.
[0050]
(1) First test:
In the operation explanatory diagram of FIG. 3, the first test is performed. This test method is the same as the method according to the first embodiment.
[0051]
(2). Second test:
As shown in the operation explanatory diagram of FIG. 5, the control parameter adjustment system 3 outputs the stored deviation e 1 (t), which is the response data of the first test, as the test signal w 2 (t). Then, the second test is performed with this test signal as the target value r (t), and the manipulated variable u 2 (t) and the output y 2 (t) are measured. Further, 1 / C (x) and ∂C (x) / ∂x are obtained by partial differentiation in the signal converters 6 and 8 for the manipulated variable u 2 (t) and the output y 2 (t), respectively, and ∂y ( t) / ∂x and ∂u (t) / ∂x are generated and input to the control parameter adjustment system 3.
[0052]
(3). Calculation of ∂J (x) / ∂x:
The control parameter adjustment system 3 calculates ∂J (x) / ∂x by substituting the acquired ∂y (t) / ∂x and ∂u (t) / ∂x into the equation (4). The e (t) and u (t) used in the equation (4) use the first test results e 1 (t) and u 1 (t), and ∂y (t) / ∂x and ∂u (t) / ∂x uses a value generated from response data obtained from the second test.
[0053]
(4). Update parameters:
The control parameter adjustment system 3 obtains x k + 1 by sequentially updating the control parameter x according to the equation (5). λ is an appropriate advance width.
[0054]
In Eq. (5), Rk may be a unit matrix or may be obtained from Eq. (8) using the Gauss-Newton method using ∂y (t) / ∂x and ∂u (t) / ∂x. May be.
[0055]
As described above, according to the second embodiment, the expressions (9) and (10) are calculated. Compared with the first embodiment that calculates the expressions (6) and (7), Although the calculation itself becomes complicated, both tests can be performed by changing the target value signal. This simplifies the signal interface between the controller and the control parameter adjustment system 3. In particular, since many general-purpose controllers can normally input a target value signal from the outside, it is not necessary to change the controller interface or the control program. For this reason, it can apply as it is, without the change of hardware or software, to the control apparatus 2 already mounted in the field.
[0056]
(Third embodiment)
The present embodiment is a case where the evaluation index is the expression (2) in the parameter adjustment method according to the first and second embodiments. When the evaluation index is the expression (2), the operation amount evaluation described above may be changed from u (t) to u (t) -u (t-1).
[0057]
If the control target is a localization system, i.e., a system in which the manipulated variable also changes steadily in order to change the steady value of the controlled variable, the manipulated variable u (t ) Will not be zero over time. Since the deviation e (t) goes to zero, the manipulated variable u (t) is greatly evaluated when the weighting factor λ (t) is not zero in the evaluation function of the above-described equation (1). The control parameter to be reduced is selected. For this reason, a desired response may not be obtained in a localization system.
[0058]
However, in the third embodiment, since the evaluation function of the above equation (2) is adopted, the deviation e (as the time elapses with respect to the step-like change of the target value regardless of the characteristics of the controlled object. Since the difference between t) and the manipulated variable u (t) -u (t-1) is both zero, the deviation and manipulated variable evaluation can be adjusted with the weighting factor λ (t), eliminating the above-mentioned problem in the localization system. can do.
[0059]
(Fourth embodiment)
This embodiment is a case where the evaluation index is set to the expression (3) in the parameter adjustment method described in the first and second embodiments. When the evaluation index is the expression (3), the evaluation of the manipulated variable described above may be changed from u (t) to u (t) -u .
[0060]
According to the fourth embodiment, the deviation e (t) and the change u (t) -u of the manipulated variable in the evaluation function go to zero with the passage of time. The same effect as that of the embodiment can be obtained.
[0061]
(Fifth embodiment)
Next, FIG. 2 shows an adjustment method according to the fifth embodiment. The difference between this embodiment and the second embodiment is that, in the case of the second embodiment, the control parameter adjustment system 3 changes the target value r (t) in the first test. In the case of the present embodiment, the test signal w (t) at the exit of the manipulated variable is changed in the first test, and the result of the second test is the same as in FIG. The principle is the same as that in the second embodiment.
[0062]
In the fifth embodiment, the method using the equations (6) and (7) to obtain た め y (t) / ∂x and ∂u (t) / ∂x in the equation (4) is controlled. The parameter update is as follows.
[0063]
(1). First test:
In the operation explanatory diagram of FIG. 6, a test signal w (t) is generated from the control parameter adjustment system 3, and the target value is changed. The deviation e (t) and the manipulated variable u (t) at this time are measured and stored in the control parameter adjustment system 3 as e 1 (t) and u 1 (t).
[0064]
(2). Second test:
In the operation explanatory diagram of FIG. 7, the second test is performed using the deviation e 1 (t) stored as a result of the first test as the test signal w (t), and the operation amount u (t) and the output y Measure (t). These measured quantities are partially differentiated by the signal converters 5 and 7 to calculate ∂Cx / ∂x, respectively, and the generated ∂y (t) / ∂x and ∂u (t) / ∂x are used as control parameters. Input to the adjustment system 3.
[0065]
(3). Calculation of ∂J (x) / ∂x:
The control parameter adjustment system 3 calculates ∂J (x) / ∂x according to the equation (4). In the equation (4), e (t) and u (t) use the first test results e 1 (t) and u 1 (t), and ∂y (t) / ∂x and ∂u (t) / ∂x uses the value generated using the response data obtained from the second test.
[0066]
(4). Update parameters:
The control parameter adjustment system 3 substitutes the value obtained by the equation (4) into the equation (5), obtains the control parameter x, and sequentially updates it. In equation (5), λ is an appropriate advance width. Rk may be a unit matrix, or the above equation (8) based on the Gauss-Newton method using て y (t) / ∂x and ∂u (t) / ∂x.
[0067]
In the fifth embodiment described above, the following effects can be obtained in addition to the same effects as those of the first embodiment regarding the calculation. That is, in the case of the first embodiment to the fourth embodiment, the target value is changed in the first test, and accordingly, the control amount is changed to match the target value. For this reason, the operation state of the controlled object changes.
[0068]
However, in the fifth embodiment, although the control amount is temporarily changed by giving a test signal to the operation amount, the control amount returns to the original value because the target value is constant. For this reason, the effect that the influence which it has on the driving | running state of a control object can be made small can be acquired.
[0069]
(Sixth embodiment)
A sixth embodiment will be described with reference to the control system configuration diagram shown in FIG. In the present embodiment, control parameters are adjusted using equations (9) and (10) in order to obtain ∂y (t) / ∂x and ∂u (t) / ∂x in equation (4). It is related to the method.
[0070]
(1). First test:
The first test is the same as in the fifth embodiment, and the operation explanatory diagram is as shown in FIG.
[0071]
(2). Second test:
As shown in the operation explanatory diagram of FIG. 8, the test is performed with the target value r (t) as a deviation e 1 (t) between the target value of the first test and the output, and the manipulated variable u (t) and output y ( t) and the signal converters 6 and 8 calculate 1 / C (x) and ∂C (x) / ∂x respectively, and ∂y (t) / ∂x and ∂u (t) / ∂ x is generated, input to the control parameter adjustment system 3, and stored.
[0072]
(3). Calculation of ∂J (x) / ∂x:
The control parameter adjustment system 3 calculates ∂J (x) / ∂x by substituting the acquired ∂y (t) / ∂x and ∂u (t) / ∂x into the equation (4). The e (t) and u (t) used in the equation (4) are the response data e 1 (t) and u 1 (t) of the first test, and ∂y (t) / ∂x and ∂u (t ) / ∂x uses a value generated from response data obtained from the second test.
[0073]
(4). Update parameters:
The control parameter adjustment system 3 obtains x k + 1 by sequentially updating the control parameter x according to the equation (5). λ is an appropriate advance width.
In Eq. (5), Rk may be a unit matrix or may be obtained from Eq. (8) using the Gauss-Newton method using ∂y (t) / ∂x and ∂u (t) / ∂x. May be.
[0074]
As described above, in the case of the sixth embodiment, the target value changes in the second test, but since the control amount goes to zero in the first test, the target value in the second test also becomes zero. It becomes a signal to go. For this reason, the influence which it has on a driving | running state becomes small, and the effect similar to 5th Embodiment can be acquired.
[0075]
(Seventh embodiment)
The present embodiment is a case where, in the parameter adjustment method described above, the evaluation index is represented by equation (2). When the evaluation index is the expression (2), the operation amount evaluation described above may be changed from u (t) to u (t) -u (t-1).
According to the seventh embodiment, it is possible to obtain an effect obtained by combining the third embodiment and the fifth embodiment or the sixth embodiment.
[0076]
(Eighth embodiment)
The parameter adjustment method according to the present embodiment relates to a parameter adjustment method when the evaluation index is expressed by equation (3). In this case, may be changed evaluation of the operation amount from u (t) to u (t) -u ∞.
According to the eighth embodiment, an effect obtained by combining the fourth, fifth, or sixth embodiment can be obtained.
[0077]
【The invention's effect】
As described above, according to the present invention, a control parameter that minimizes an evaluation function that quantitatively represents a control index can be obtained.
Moreover, since the sensitivity (gradient) about the control parameter of the evaluation function is used to update the control parameter, the convergence to the optimum value can be accelerated.
[0078]
Furthermore, since the sensitivity (gradient) for the control parameter of the evaluation function is obtained, the convergence state of the control parameter can be grasped quantitatively.
Furthermore, since the test is performed in a closed loop, it can be applied to an unstable control target.
[Brief description of the drawings]
FIG. 1 is a diagram showing an example of a system configuration example for explaining the present invention.
FIG. 2 is a diagram showing another example of a system configuration example for explaining the present invention.
FIG. 3 is a diagram showing an operation during a first test of parameter adjustment in the first embodiment of the invention.
FIG. 4 is a diagram showing an action during a second test of parameter adjustment in the first embodiment of the invention.
FIG. 5 is a diagram illustrating an operation during a second test of parameter adjustment according to the second embodiment of the present invention.
FIG. 6 is a diagram showing an effect during a first test of parameter adjustment in the fifth and sixth embodiments of the present invention.
FIG. 7 is a diagram illustrating an operation during a second test of parameter adjustment according to the fifth embodiment of the present invention.
FIG. 8 is a diagram illustrating an operation during a second test of parameter adjustment according to a sixth embodiment of the present invention.
FIG. 9 is a configuration diagram for explaining an example of a conventional control parameter update method.
FIG. 10 is a configuration diagram for explaining another example of a conventional control parameter update method.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 1 ... Control object, 2 ... Control apparatus, 3 ... Control parameter adjustment system, 4-1 to 4-3 ... Signal addition part, 5-8 ... Signal conversion part.

Claims (7)

制御対象を制御する制御装置と、この制御装置の制御パラメータを調整するとともにテスト信号を発生させる制御パラメータ調整システムと、前記制御装置から出力された操作量を偏微分して、操作量の操作パラメータについての変化率を求め、この変化率を前記制御パラメータ調整システムに入力する信号変換部と、前記制御対象から出力される制御量を偏微分して、制御量の、操作パラメータについての変化率を求め、この変化率を前記制御パラメータ調整システムに出力する信号変換部とから制御系を構成し、制御量を目標値に一致させるように制御対象を制御する制御装置の制御パラメータ調整方法において、
前記制御系を制御対象と制御装置とを組み合わせた閉ループの状態で1回目のテスト信号を加え、そのときの応答データを用いて2回目のテスト信号を作り、この2回目のテスト信号を加えたときの応答データを用いて制御パラメータの変化量を決め、これら2回一組としたテストを繰り返し、設定した評価関数を最小とする制御パラメータを決定することを特徴とする制御装置の制御パラメータ調整方法。
A control device for controlling a control target, a control parameter adjustment system for adjusting a control parameter of the control device and generating a test signal, and an operation parameter for an operation amount by partially differentiating an operation amount output from the control device The rate of change of the control amount with respect to the operation parameter is obtained by partially differentiating the control amount output from the control target and the signal conversion unit that inputs the rate of change to the control parameter adjustment system. In the control parameter adjustment method of the control device for controlling the control object so as to match the control amount with the target value , comprising a control system from the signal conversion unit that calculates the rate of change to the control parameter adjustment system .
A first test signal is added in a closed loop state in which the control system and the control device are combined in the control system, a second test signal is created using the response data at that time, and the second test signal is added. The control parameter adjustment of the control device is characterized in that the amount of change of the control parameter is determined using the response data at the time, the test is repeated twice as a set, and the control parameter that minimizes the set evaluation function is determined. Method.
前記評価関数は、目標値に対する前記制御対象の出力の望ましい変化と実際に得られる制御対象の出力との偏差の2乗と操作量の変化分の2乗に重みをかけた量の和をテスト期間についての加算平均値とすることを特徴とする請求項1記載の制御装置の制御パラメータ調整方法。The evaluation function tests a sum of amounts obtained by weighting the square of the deviation between the desired change in the output of the control target with respect to the target value and the actual output of the control target and the square of the change in the manipulated variable. control parameter adjusting method of a control apparatus according to claim 1, wherein the average value and to Rukoto for periods. 前記評価関数は、目標値に対する前記制御対象の出力の望ましい変化と実際に得られる制御対象の出力との偏差の2乗と、操作量と操作量の将来値との差の2乗に重みをかけた量の和をテスト期間についての加算平均値とすることを特徴とする請求項1 記載の制御装置の制御パラメータ調整方法。The evaluation function gives a weight to the square of the difference between the desired change in the output of the controlled object relative to the target value and the actual output of the controlled object, and the square of the difference between the manipulated variable and the future value of the manipulated variable. control parameter adjusting method of a control apparatus according to claim 1, wherein the average value and to Rukoto for the test period the sum of the amount applied. 前記1回目のテストの際、目標値を変化することを特徴とする請求項1記載の制御装置の制御パラメータ調整方法。  2. The control parameter adjustment method for a control device according to claim 1, wherein the target value is changed during the first test. 前記1回目のテストの際、操作量にテスト信号を加えることを特徴とする請求項1の制御装置の制御パラメータ調整方法。  2. The control parameter adjustment method for a control apparatus according to claim 1, wherein a test signal is added to the operation amount in the first test. 前記2回目のテストの際、1回目のテスト目標値に対する前記制御対象の出力の望ましい変化と実際に得られる制御対象の出力との偏差を、新しい目標値とすることを特徴とする請求項1の制御装置の制御パラメータ調整方法。  The difference between a desired change in the output of the controlled object with respect to the first test target value and the actually obtained output of the controlled object is set as a new target value in the second test. Method for adjusting control parameters of the control apparatus. 前記2回目のテストの際、1回目のテスト目標値に対する前記制御対象の出力の望ましい変化と実際に得られる制御対象の出力との偏差を、操作量に加えることを特徴とする請求項1記載の制御装置の制御パラメータ調整方法。  The difference between a desired change in the output of the controlled object with respect to the first test target value and an actually obtained output of the controlled object is added to the manipulated variable during the second test. Method for adjusting control parameters of the control apparatus.
JP2003020706A 2003-01-29 2003-01-29 Control parameter adjustment method for control device Expired - Fee Related JP4155508B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2003020706A JP4155508B2 (en) 2003-01-29 2003-01-29 Control parameter adjustment method for control device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003020706A JP4155508B2 (en) 2003-01-29 2003-01-29 Control parameter adjustment method for control device

Publications (2)

Publication Number Publication Date
JP2004234219A JP2004234219A (en) 2004-08-19
JP4155508B2 true JP4155508B2 (en) 2008-09-24

Family

ID=32950260

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003020706A Expired - Fee Related JP4155508B2 (en) 2003-01-29 2003-01-29 Control parameter adjustment method for control device

Country Status (1)

Country Link
JP (1) JP4155508B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5467299B2 (en) * 2007-11-30 2014-04-09 エフイーアイ カンパニー Iterative feedback adjustment in a scanning probe microscope
CN112578674B (en) * 2020-12-29 2022-05-17 润电能源科学技术有限公司 Excitation signal adjusting method and device of control system and readable storage medium

Also Published As

Publication number Publication date
JP2004234219A (en) 2004-08-19

Similar Documents

Publication Publication Date Title
Xia An inverse model of MR damper using optimal neural network and system identification
JP4918207B2 (en) Computer system and method for constraining nonlinear approximation of empirical processes.
JP4908433B2 (en) Control parameter adjustment method and control parameter adjustment program
US5553195A (en) Dynamic neural net
Guzmán et al. Interactive teaching of constrained generalized predictive control
CN108682023A (en) Close coupling Unscented kalman tracking filter algorithm based on Elman neural networks
JP2007500379A (en) Software optimization device for intelligent control system
Fischer et al. Predictive control based on local linear fuzzy models
JP2019021186A (en) Data processor, control system, data processing method, and program
US7123971B2 (en) Non-linear model with disturbance rejection
Fu et al. Physics-data combined machine learning for parametric reduced-order modelling of nonlinear dynamical systems in small-data regimes
Tsokanas et al. Adaptive model predictive control for actuation dynamics compensation in real-time hybrid simulation
JP4155508B2 (en) Control parameter adjustment method for control device
JP3214876B2 (en) Neural network configuration method and neural network construction support system
Du et al. Reinforcement learning
JPH06332506A (en) Nonlinear controller
Euler-Rolle et al. Automated generation of feedforward control using feedback linearization of local model networks
CN110909492A (en) Sewage treatment process soft measurement method based on extreme gradient lifting algorithm
JP3740382B2 (en) Process model parameter adjusting apparatus and method
JP2674167B2 (en) Process control method using neural network model
CN114139442A (en) Method for optimizing steam turbine valve flow characteristic function based on K-means-LSTM
RADBAKHSH et al. Physics-Informed Neural Network for Analyzing Elastic Beam Behavior
Bal et al. JMASM 55: MATLAB Algorithms and Source Codes of'cbnet'Function for Univariate Time Series Modeling with Neural Networks (MATLAB)
CN117807403B (en) Steel transformation control method and device based on behavior cloning, medium and computer equipment
Escudero et al. Adaptive linearization for nonlinear systems using continuous Neural Networks

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20050404

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20060829

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20070222

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20070620

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070703

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070830

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080205

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080331

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: 20080702

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20080704

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110718

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110718

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120718

Year of fee payment: 4

LAPS Cancellation because of no payment of annual fees