JP2012113676A - Configuring method, system, and program for controller - Google Patents

Configuring method, system, and program for controller Download PDF

Info

Publication number
JP2012113676A
JP2012113676A JP2010264660A JP2010264660A JP2012113676A JP 2012113676 A JP2012113676 A JP 2012113676A JP 2010264660 A JP2010264660 A JP 2010264660A JP 2010264660 A JP2010264660 A JP 2010264660A JP 2012113676 A JP2012113676 A JP 2012113676A
Authority
JP
Japan
Prior art keywords
max
mid
candidate
controller
plant
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2010264660A
Other languages
Japanese (ja)
Other versions
JP5713338B2 (en
Inventor
Junya Shimizu
淳也 清水
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.)
International Business Machines Corp
Original Assignee
International Business Machines 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 International Business Machines Corp filed Critical International Business Machines Corp
Priority to JP2010264660A priority Critical patent/JP5713338B2/en
Publication of JP2012113676A publication Critical patent/JP2012113676A/en
Application granted granted Critical
Publication of JP5713338B2 publication Critical patent/JP5713338B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a technique for automatically configuring a controller with respect to each candidate for a scheduling parameter vector.SOLUTION: An LPV model of a plant is described, the number of scheduling parameters is given, and a maximum value and a minimum value found from a high-order request of restriction conditions are given to each element of a scheduling parameter vector comprising a plurality of elements. Each candidate has a ν gap calculated among Min, Max, and an intermediate value between them respectively, respective candidates are determined as scheduling parameter candidates in the decreasing order of differences between ν gaps. The number of parameters given first is selected and converted into an LMI, the controller is found based on an Hbasis for all combination end points associated with the maximum value and minimum value of the selected scheduling parameter, and a controller for a generalized plant is constituted by bilinear interpolation weighted with the ν gap calculated during the parameter selection.

Description

この発明は、自動車のプラントなどの制御器を構成する技法に関し、より詳しくは、線形パラメータ変動モデルによってプラントを記述し、それを線形行列不等式に変換することでゲインスケジュール制御を行うための技法に関する。   The present invention relates to a technique for configuring a controller such as an automobile plant, and more particularly to a technique for performing gain schedule control by describing a plant by a linear parameter variation model and converting it into a linear matrix inequality. .

従来より、動作中に構成パラメータが変化するプラントに対する制御として、多くの場合、データマップによりオブザーバのゲインを調整して状態変数を推定しながら制御を行う、可変構造制御の方法が取られている。   Conventionally, as a control for a plant whose configuration parameters change during operation, in many cases, a variable structure control method has been used in which control is performed while adjusting an observer gain using a data map and estimating a state variable. .

しかし、データマップに全ての可能性を盛り込むことは多大な労力を要し、また、実験的にマップを作成する場合には、プラントのパラメータ変更への対応が難しい。   However, incorporating all the possibilities into the data map requires a great deal of labor, and when creating a map experimentally, it is difficult to respond to plant parameter changes.

そこで、線形パラメータ変動(LPV: Linear Parameter-Varying)モデルによってプラントを記述し、変動範囲制約を線形行列不等式(LMI: Linear Matrix Inequalities)に変換することで全ての場合を網羅するゲインスケジューリング制御が一つの解として検討されてきた。   Therefore, a plant is described using a linear parameter-variable (LPV) model, and gain scheduling control covering all cases is achieved by converting the fluctuation range constraint into linear matrix inequalities (LMI). Has been considered as a solution.

特開平10−220475号公報は、強いジャイロ効果が作用するLPVシステムに対し、全ての回転数に渡り連続して最適に安定な制御の可能な磁気軸受装置を提供するために、LMIベースで設計され、凸最適化を満足する、異なる回転数を対象とした2台のH補償器K(θmax)とK(θmin)を配設し、この各H補償器は、それぞれ内挿補間により求めた係数α1 及びα2 を乗じ、並列に接続する。両H補償器の出力は加算し、係数α1 及びα2 はその総和が1となるように、また回転速度センサ16で計測した回転数ωz に依存しつつ連続的に変化させ、半径方向位置センサ8により計測した半径方向位置と指令値r間の誤差を算出し、その誤差に基づきLMIゲインスケジュールH演算を行う技法を開示する。 Japanese Patent Laid-Open No. 10-220475 is designed on an LMI basis to provide a magnetic bearing device capable of optimally and continuously controlling over all rotation speeds for an LPV system in which a strong gyro effect acts. Two H compensators K (θ max ) and K (θ min ) for different rotational speeds satisfying the convex optimization are arranged, and each H compensator is interpolated. Multiply by the coefficients α1 and α2 found by interpolation and connect them in parallel. The outputs of both H compensators are added, and the coefficients α1 and α2 are continuously changed so that the sum thereof is 1 and depending on the rotational speed ωz measured by the rotational speed sensor 16, and the radial position sensor A technique for calculating an error between the radial position measured by 8 and the command value r and performing an LMI gain schedule H∞ calculation based on the error is disclosed.

特開平11−110003号公報は、パラメータ変動の上限値、下限値により定まる範囲の中において、目標応答を入力(指定)することにより、コントローラを導出し、目標応答を満足できる制御性能を得るビークル制御系の設計方法を提供することを目的とし、(A)ビークルの諸元と、外乱条件と、パラメータ変動範囲と、目標応答を設定し、(B)ビークル1の力学モデルを作成し、(C)ビークル1の力学モデルに基づきLPVモデル化を行い、(D)前記設定条件を逆LMI問題として定式化して、LMI設計法によりコントローラを導出し、(E)目標応答の標準であるステップ応答を制約条件とする制御系設計アルゴリズムをLMIステップ制約サーボ系の構成で作成し、(F)シミュレータによるコントローラ制御性能の検証を行い、(G)検証結果が不十分であれば、目標応答設定の調整に戻り、検証結果が十分であれば、制御系を完成させることを開示する。   Japanese Patent Application Laid-Open No. 11-110003 discloses a vehicle that obtains control performance that satisfies a target response by deriving a controller by inputting (designating) a target response within a range determined by an upper limit value and a lower limit value of parameter fluctuations. Aiming at providing a control system design method, (A) setting the specifications of the vehicle, disturbance conditions, parameter fluctuation range, and target response, (B) creating a dynamic model of the vehicle 1 ( C) LPV modeling is performed based on the dynamic model of the vehicle 1, (D) the setting condition is formulated as an inverse LMI problem, a controller is derived by the LMI design method, and (E) a step response that is a standard of target response A control system design algorithm with a constraint condition is created with an LMI step-constrained servo system configuration. Performed testimony, if insufficient (G) verification result back to the adjustment of the target response setting, the verification result if sufficient, disclose that complete the control system.

しかし、上記従来技術においては、十分な性能を得るためのスケジューリング・パラメータ候補の選択が難しく、その選択は、特別な計測をしたり、ボード線図を眺めながら主観的に決定するなど、人手に頼っていたので、時間がかかるとともに、なかなか、最適なパラメータ候補の選択の選択に至らなかった。   However, in the above prior art, it is difficult to select a scheduling parameter candidate for obtaining sufficient performance, and the selection is made manually, such as by making a special measurement or subjectively determining while looking at the Bode diagram. Since it relied on, it took time and it was difficult to select an optimal parameter candidate.

ところで、S.Aubouet et al.,”H-infinity/LPV observer for an industrial semi-active suspension,” IEEE Proc. Conf. on Control Applications, pp.756-763, 2009には、LPVモデルと各制約条件を考慮したパラメータの変動幅を表現し、LMIに変換してH基準に基づく制御則を各パラメータの最大値・最小値の組合せ端点において求める技法が記述されている。 By the way, S. Aubouet et al., “H-infinity / LPV observer for an industrial semi-active suspension,” IEEE Proc. Conf. On Control Applications, pp.756-763, 2009 includes the LPV model and each constraint. Describes a technique for expressing a control variation based on the H criterion at a combination end point of the maximum value / minimum value of each parameter by expressing the fluctuation range of the parameter in consideration of.

また、G.Vinnicombe, “Frequency Domain Uncertainty and the Graph Topology,” IEEE Trans. Automatic Control, Vol.38, No.9, pp.1371-1383, 1993には、プラントパラメータベクトルの各候補について、νギャップを計算する技法が記述されている。   Also, G. Vinnicombe, “Frequency Domain Uncertainty and the Graph Topology,” IEEE Trans. Automatic Control, Vol. 38, No. 9, pp.1371-1383, 1993 describes the ν gap A technique for computing is described.

特開平10−220475号公報Japanese Patent Laid-Open No. 10-220475 特開平11−110003号公報Japanese Patent Laid-Open No. 11-110003

S.Aubouet et al.,”H-infinity/LPV observer for an industrial semi-active suspension,” IEEE Proc. Conf. on Control Applications, pp.756-763, 2009S. Aubouet et al., “H-infinity / LPV observer for an industrial semi-active suspension,” IEEE Proc. Conf. On Control Applications, pp.756-763, 2009 G.Vinnicombe, “Frequency Domain Uncertainty and the Graph Topology,” IEEE Trans. Automatic Control, Vol.38, No.9, pp.1371-1383, 1993G. Vinnicombe, “Frequency Domain Uncertainty and the Graph Topology,” IEEE Trans. Automatic Control, Vol.38, No.9, pp.1371-1383, 1993

この発明の目的は、スケジューリング・パラメータ・ベクトルの各候補について計算されたνギャップに基づき、制御器を自動的に構成するための技法を提供することにある。   It is an object of the present invention to provide a technique for automatically configuring a controller based on the ν gap calculated for each candidate scheduling parameter vector.

本発明によれば、ユーザーが先ず、プラントのLPVモデルを記述し、スケジューリング・パラメータ(θの要素)とその数を与える。複数の要素からなるスケジューリング・パラメータ・ベクトルの各々の要素には、制約条件の上位要求から出てくる最大値と最小値が与えられる。記述されたLPVモデルと、スケジューリング・パラメータ・ベクトル、及びその要素の最大値と最小値は、本発明の処理を実行するシステムに保存される。   According to the present invention, the user first describes the plant LPV model and gives scheduling parameters (elements of θ) and their numbers. Each element of the scheduling parameter vector composed of a plurality of elements is given a maximum value and a minimum value that come out from the upper request of the constraint condition. The described LPV model, scheduling parameter vector, and maximum and minimum values of its elements are stored in a system that performs the processing of the present invention.

システムは、スケジューリングパラメータベクトルの各候補について、そのMin, Max、さらにそれらの中間値の各々の間について、νギャップを計算する。   For each candidate scheduling parameter vector, the system calculates a ν gap for each of its Min, Max, and their intermediate values.

システムは次に、各候補についてνギャップ間の差が大きいものから順番にスケジューリング・パラメータ候補とし、最初に与えられたパラメータ数を選ぶ。   Next, the system selects scheduling parameters candidates in descending order of the difference between the ν gaps for each candidate, and selects the number of parameters given first.

システムは次に、LPVモデルをLMIに変換し、選定したスケジューリング・パラメータの最大値と最小値に関する全ての組み合わせ端点についてH基準によって制御器を求める。 The system then converts the LPV model to LMI and determines the controller according to the H∞ criterion for all combined endpoints for the maximum and minimum values of the selected scheduling parameter.

システムは次に、パラメータ選定時に計算したνギャップによって重み付けした双一次補間によって一般化プラントへの制御器を構成する。   The system then configures the controller to the generalized plant by bilinear interpolation weighted by the ν gap calculated during parameter selection.

この発明によれば、多くの制約条件が存在する場合に、モデル設計者に依存してきたスケジューリングパラメータを客観的基準で自動選定し、かつ、それを利用した従来よりも適切な制御器設計が可能となるため、変動範囲が漠然としか捉えられない場合にもマップ作成の労力を低減することが可能となる。   According to the present invention, when there are many constraints, scheduling parameters that have been dependent on the model designer are automatically selected based on objective criteria, and more appropriate controller design than before can be made using this. Therefore, it is possible to reduce the effort for creating the map even when the fluctuation range is only vague.

また、νギャップによるパラメータ空間特性の把握は、SysML等による要求分析との対応付けを容易にする可能性があり、事後的に制約を追加する際に制御ストラテジーを見直すべきか否かの自動判別を可能とするためモデリング設計の柔軟性を向上するために用いることができる。   In addition, grasping the parameter space characteristics by the ν gap may facilitate the association with the requirement analysis by SysML or the like, and automatically discriminates whether or not the control strategy should be reconsidered when constraints are added later. Can be used to improve the flexibility of modeling design.

さらに、線形補間でなく、双一次補間を用いることは、非ユークリッド幾何的測地距離を考慮することになり、これにより、フィードバック制御の誤差を小さくすることにつながる。   Furthermore, using bilinear interpolation instead of linear interpolation takes into account non-Euclidean geometric geodetic distances, thereby reducing feedback control errors.

本発明で使用される一例の制御コンピュータの概要ブロック図である。It is a general | schematic block diagram of an example control computer used by this invention. 本発明に従い制御器を決定する処理を実行するためのプログラム等の機能ブロック図である。It is a functional block diagram of a program or the like for executing processing for determining a controller according to the present invention. 本発明に従い制御器を決定する処理を実行するためのフローチャートを示す図である。It is a figure which shows the flowchart for performing the process which determines a controller according to this invention. LMIへの変換によるH基準制御器導出を説明するための図である。It is a figure for demonstrating Hinfinity reference | standard controller derivation | leading-out by conversion to LMI. 一般プラントからLMI問題への変換を説明するための図である。It is a figure for demonstrating conversion from a general plant to a LMI problem.

以下、図面に従って、本発明の実施例を説明する。これらの実施例は、本発明の好適な態様を説明するためのものであり、発明の範囲をここで示すものに限定する意図はないことを理解されたい。また、以下の図を通して、特に断わらない限り、同一符号は、同一の対象を指すものとする。   Embodiments of the present invention will be described below with reference to the drawings. It should be understood that these examples are for the purpose of illustrating preferred embodiments of the invention and are not intended to limit the scope of the invention to what is shown here. Further, throughout the following drawings, the same reference numerals denote the same objects unless otherwise specified.

図1を参照すると、本発明に従い制御器118を制御するための制御コンピュータのハードウェアのブロック図が示されている。図1において、システム・バス102には、CPU104と、主記憶(RAM)106と、ハードディスク・ドライブ(HDD)108と、キーボード110と、マウス112と、ディスプレイ114が接続されている。CPU104は、好適には、32ビットまたは64ビットのアーキテクチャに基づくものであり、例えば、インテル社のPentium(商標)4、インテル社のCore(商標) 2 DUO、AMD社のAthlon(商標)などを使用することができる。主記憶106は、好適には、2GB以上の容量、より好ましくは、4GB以上の容量をもつものである。   Referring to FIG. 1, a block diagram of the control computer hardware for controlling the controller 118 in accordance with the present invention is shown. In FIG. 1, a CPU 104, a main memory (RAM) 106, a hard disk drive (HDD) 108, a keyboard 110, a mouse 112, and a display 114 are connected to the system bus 102. The CPU 104 is preferably based on a 32-bit or 64-bit architecture, such as Intel Pentium (trademark) 4, Intel Core (trademark) 2 DUO, AMD Athlon (trademark), etc. Can be used. The main memory 106 preferably has a capacity of 2 GB or more, more preferably a capacity of 4 GB or more.

ハードディスク・ドライブ108には、個々に図示しないが、オペレーティング・システム及び本発明に係る処理プログラムなどが、予め格納されている。オペレーティング・システムは、Linux(商標)、マイクロソフト社のWindows Vista、Windows XP(商標)、Windows(商標)2000、アップルコンピュータのMac OS(商標)などの、CPU104に適合する任意のものでよい。   The hard disk drive 108 stores in advance an operating system, a processing program according to the present invention, and the like (not shown). The operating system may be any compatible with the CPU 104, such as Linux (trademark), Microsoft Windows Vista, Windows XP (trademark), Windows (trademark) 2000, Apple Mac OS (trademark).

キーボード110及びマウス112は、オペレーティング・システムが提供するグラフィック・ユーザ・インターフェースに従い、ディスプレイ114に表示されたアイコン、タスクバー、ウインドウなどのグラフィック・オブジェクトを操作するために使用される。キーボード110及びマウス112はまた、後述するデータは記録用プログラムを操作するためにも使用される。   The keyboard 110 and the mouse 112 are used to operate graphic objects such as icons, taskbars, and windows displayed on the display 114 in accordance with a graphic user interface provided by the operating system. The keyboard 110 and the mouse 112 are also used to operate a recording program with data described later.

ディスプレイ114は、これには限定されないが、好適には、1024×768以上の解像度をもち、32ビットtrue colorのLCDモニタである。ディスプレイ114は、制御器118によるプラントの動作の波形などを表示するために使用される。   The display 114 is preferably, but not limited to, a 32-bit true color LCD monitor with a resolution of 1024 × 768 or higher. The display 114 is used for displaying a waveform or the like of the operation of the plant by the controller 118.

ハードディスク・ドライブ108にはさらに、後述するLPVモデル204、LPVモデル204のためのスケジューリング・パラメータ206、νギャップ計算モジュール208、LMI変換モジュール210、制御器構成モジュール212及び、全体の処理を統合するメイン・プログラム202が格納されている。メイン・プログラム202、νギャップ計算モジュール208、LMI変換モジュール210、及び制御器構成モジュール212は、C、C++、C#、Java(R)などの既存の任意のプログラミング言語で作成することができる。これらのモジュールは、オペレーティング・システムの働きで適宜主記憶106にロードされて実行される。図示しないが、メイン・プログラム202は、適当なGUIにより、オペレータが操作するためのウィンドウをディスプレイ114に表示し、オペレータがキーボード110及びマウス112などを使って処理を開始し、あるいは停止することができる。   The hard disk drive 108 further includes an LPV model 204 to be described later, a scheduling parameter 206 for the LPV model 204, a ν gap calculation module 208, an LMI conversion module 210, a controller configuration module 212, and a main unit that integrates the overall processing. A program 202 is stored. The main program 202, the ν gap calculation module 208, the LMI conversion module 210, and the controller configuration module 212 can be created in any existing programming language such as C, C ++, C #, Java®. These modules are appropriately loaded into the main memory 106 and executed by the operation of the operating system. Although not shown, the main program 202 displays a window for the operator to operate on the display 114 with an appropriate GUI, and the operator can start or stop the process using the keyboard 110, the mouse 112, or the like. it can.

バス102にはさらに、PCI、USBなどの適当なインターフェース・ボード116を介して、図4に示すプラント402を制御するための制御器118が接続されている。制御器118は、インターフェース・ボード116を通じて、計算の結果得られた制御行列Kを格納するレジスタをもつ。   Further, a controller 118 for controlling the plant 402 shown in FIG. 4 is connected to the bus 102 via an appropriate interface board 116 such as PCI or USB. The controller 118 has a register for storing a control matrix K obtained as a result of calculation through the interface board 116.

プラント402は例えば、エンジン、ブレーキ、エアコンなど制御対象の機構装置である。   The plant 402 is a mechanism device to be controlled such as an engine, a brake, and an air conditioner.

プラント402は、実機であってもよいし、Simulink(R)、Modelicaなどのシミュレーション・モデリング・ツールで作成されたソフトウェア的なモデルであってもよい。プラント402がソフトウェア的なモデルである場合、ハードディスク・ドライブ108に格納されたモジュールとなり、すると、制御器118もハードディスク・ドライブ108に格納されたモジュールとなるので、インターフェース・ボード116を介したインターフェースは不要である。   The plant 402 may be a real machine or a software model created by a simulation modeling tool such as Simulink (R) or Modelica. When the plant 402 is a software model, the module is stored in the hard disk drive 108, and the controller 118 is also a module stored in the hard disk drive 108. Therefore, the interface via the interface board 116 is It is unnecessary.

次に、図2の機能ブロック図を参照して、本発明に従い制御器を決定する処理を実行するためのプログラム等の論理的構成について説明する。   Next, with reference to the functional block diagram of FIG. 2, a logical configuration of a program or the like for executing a process for determining a controller according to the present invention will be described.

メイン・プログラム202は、ユーザとのインターフェースを含め、全体の動作を制御する機能を果たす。   The main program 202 functions to control the entire operation including the interface with the user.

線形パラメータ変動(LPV)モデル204は、下記のような線形方程式により、プラント402を記述するモデルであり、力学方程式などを考慮して、予めユーザにより記述される。

Figure 2012113676

ここで、A、B、C、Dは、状態空間モデルの係数行列である。 The linear parameter variation (LPV) model 204 is a model that describes the plant 402 by the following linear equation, and is described in advance by the user in consideration of the dynamic equation and the like.
Figure 2012113676

Here, A, B, C, and D are coefficient matrices of the state space model.

スケジューリング・パラメータ206は、LPVモデル204で使用される時変パラメータθ(t)であり、次の式で示すように、r次元のベクトルである。

Figure 2012113676

θ(t)の各要素θj(t) (j = 1,2,...,r)は、システム要件に従い、予めユーザにより記述される。 The scheduling parameter 206 is a time-varying parameter θ (t) used in the LPV model 204, and is an r-dimensional vector as shown by the following equation.
Figure 2012113676

Each element θ j (t) (j = 1, 2,..., r) of θ (t) is described by the user in advance according to the system requirements.

また、θ(t)の各要素θj(t) (j = 1,2,...,r)には、下記の式に示すように、制約条件の上位要求から出てくる最大値θi maxと、最小値θi minが予めユーザにより決定され、スケジューリング・パラメータ206としてハードティスク・ドライブ108に保存されている。

Figure 2012113676
In addition, each element θ j (t) (j = 1,2, ..., r) of θ (t) has the maximum value i max and the minimum value θ i min are determined in advance by the user and stored in the hard disk drive 108 as the scheduling parameter 206.
Figure 2012113676

νギャップ計算モジュール208は、スケジューリング・パラメータ206の最大値θj max、最小値θj min及びそれらの間の中間値θj midとから、G.Vinnicombe, “Frequency Domain Uncertainty and the Graph Topology,” IEEE Trans. Automatic Control, Vol.38, No.9, pp.1371-1383, 1993に記述されているような方法でνギャップを計算する機能をもつ。なお、中間値θj midは、次のようにして計算される。

Figure 2012113676
The ν gap calculation module 208 calculates G. Vinnicombe, “Frequency Domain Uncertainty and the Graph Topology,” from the maximum value θ j max , the minimum value θ j min and the intermediate value θ j mid between them. It has a function to calculate the ν gap by the method described in IEEE Trans. Automatic Control, Vol.38, No.9, pp.1371-1383, 1993. The intermediate value θ j mid is calculated as follows.
Figure 2012113676

線形行列不等式(LMI)変換モジュール210は、LPVモデル204をLMIに変換し、選定されたスケジューリング・パラメータの最大値と最小値を用いて、H基準により、制御器を求める機能をもつ。 The linear matrix inequality (LMI) conversion module 210 has a function of converting the LPV model 204 to LMI and obtaining a controller based on the H∞ criterion using the maximum and minimum values of the selected scheduling parameter.

制御器構成モジュール212は、パラメータ選定時に計算したνギャップによって重み付けした双一次補間によって、一般化プラントへの制御器118を構成する。この際、制御器118の構成は、図1に示すインターフェース・ボード116を介して行われる。   The controller configuration module 212 configures the controller 118 to the generalized plant by bilinear interpolation weighted by the ν gap calculated at the time of parameter selection. At this time, the controller 118 is configured via the interface board 116 shown in FIG.

次に、図3のフローチャートを参照して、本発明に従い制御器を決定する処理について説明する。図3において、ステップ302では、ユーザが、LPVモデルを記述し、またスケジューリング・パラメータとその要素数を与える。その際、スケジューリング・パラメータの各要素の最大値と最小値も、仕様に従い入力される。それらはそれぞれ、メイン・プログラム202によって、図2のLPVモデル204及びスケジューリング・パラメータ206として、ハードディスク・ドライブ108に保存される。   Next, processing for determining a controller according to the present invention will be described with reference to the flowchart of FIG. In FIG. 3, in step 302, a user describes an LPV model and gives scheduling parameters and the number of elements thereof. At this time, the maximum value and the minimum value of each element of the scheduling parameter are also input according to the specification. Each of them is stored in the hard disk drive 108 by the main program 202 as the LPV model 204 and scheduling parameter 206 of FIG.

ステップ304では、メイン・プログラム202がνギャップ計算モジュール206を呼び出して、制約条件の上位要求から出てくるスケジューリング・パラメータ・ベクトルの各候補に関して、Min、Max、さらにその中間値の各々について、例えば、G.Vinnicombe, “Frequency Domain Uncertainty and the Graph Topology,” IEEE Trans. Automatic Control, Vol.38, No.9, pp.1371-1383, 1993に記述されている技法により、νギャップが計算される。   In step 304, the main program 202 calls the ν gap calculation module 206, and for each candidate of the scheduling parameter vector coming out of the upper request of the constraint condition, for each of Min, Max and its intermediate value, for example, , G. Vinnicombe, “Frequency Domain Uncertainty and the Graph Topology,” IEEE Trans. Automatic Control, Vol.38, No.9, pp.1371-1383, 1993 is used to calculate the ν gap. .

先ず、一般的なνギャップの計算について説明する。2つのプラント・モデルの伝達関数をP1(jω)、P2(jω)とすると、その2つの伝達関数の間のνギャップは、次の式で定義される。ここで、ωは角周波数である。

Figure 2012113676
First, general calculation of the ν gap will be described. If the transfer functions of two plant models are P 1 (jω) and P 2 (jω), the ν gap between the two transfer functions is defined by the following equation. Here, ω is an angular frequency.
Figure 2012113676

ここで、κ(X,Y)は次のように定義されている。

Figure 2012113676

また、
Figure 2012113676

は、最大特異値をあらわす。 Here, κ (X, Y) is defined as follows.
Figure 2012113676

Also,
Figure 2012113676

Represents the maximum singular value.

1入力1出力系であれば、νギャップは次のようになる。

Figure 2012113676
In the case of a 1-input 1-output system, the ν gap is as follows.
Figure 2012113676

ここで、θのi番目(i = 1,...,r)の要素について最大値を代入した場合の伝達関数をTi max、最小値を代入した場合の伝達関数をTi minとする。また、その中間値を代入した場合の伝達関数をTi midとする。 Here, T i max is the transfer function when the maximum value is substituted for the i-th (i = 1, ..., r) element of θ, and T i min is the transfer function when the minimum value is substituted. . Also, let T i mid be the transfer function when the intermediate value is substituted.

そこで、メイン・プログラム202は、上記νギャップの計算式を用いて、νギャップ計算モジュール206を呼び出して、次値を計算する。
νi_mid_max = δν(Ti mid,Ti max)
νi_min_mid = δν(Ti min,Ti mid)
νi_min_max = δν(Ti min,Ti max)
Therefore, the main program 202 calls the ν gap calculation module 206 using the above ν gap calculation formula and calculates the next value.
ν i _mid_max = δ ν (T i mid , T i max )
ν i _min_mid = δ ν (T i min , T i mid )
ν i _min_max = δ ν (T i min , T i max )

メイン・プログラム202は、ステップ306で、
(|νi_min_mid - νi_mid_max| + νi_min_max) (ただしi = 1,...,r)の値が大きいiから順番にスケジューリング・パラメータ候補とし、最初に与えられたパラメータ数を選ぶ。
The main program 202 is step 306,
(| ν i _min_mid-ν i _mid_max | + ν i _min_max) (where i = 1, ..., r), starting from i with the largest value, select the scheduling parameter candidates and select the number of parameters given first. .

このとき、第1の実施例では、最初に与えられたパラメータのみが選ばれ、第2の実施例では、次に与えられたパラメータも選ばれる。この第1の実施例と第2の実施例について、後で説明する。   At this time, in the first embodiment, only the first given parameter is selected, and in the second embodiment, the next given parameter is also selected. The first embodiment and the second embodiment will be described later.

ステップ308では、メイン・プログラム202は、LMI変換モジュール210を呼び出して、LPVモデル204をLMIに変換し、H基準制御器を導出する。 In step 308, the main program 202 calls the LMI conversion module 210 to convert the LPV model 204 to LMI and derives an H∞ reference controller.

そこでまず、その一般論を説明する。すなわち、LPVモデル204の変動部分を外生信号として捉え、さらにそれを大きなプラントを考え、それを図4に示すように、一般化プラント404と定義する。外生信号の内部フィルタや外生信号そのものは、LPVモデル204の、どのパラメータをスケジューリング・パラメータとして選ぶかが決定されれば、その周波数特性から自動的に決定することができる。   First, the general theory will be explained. That is, the fluctuation part of the LPV model 204 is regarded as an exogenous signal, and further considered as a large plant, which is defined as a generalized plant 404 as shown in FIG. The internal filter of the exogenous signal and the exogenous signal itself can be automatically determined from their frequency characteristics if it is determined which parameter of the LPV model 204 is selected as the scheduling parameter.

そこで、

Figure 2012113676

と表現できる一般化プラント404において、外生信号wに対して制御量zをなるべく小さく抑えることが制御目的なので、この間の伝達関数Gzwを小さくする制御器Kを設計することになる。ここで、
Figure 2012113676

である。 Therefore,
Figure 2012113676

In the generalized plant 404 that can be expressed as follows, since the control purpose is to keep the control amount z as small as possible with respect to the exogenous signal w, a controller K that reduces the transfer function G zw during this period is designed. here,
Figure 2012113676

It is.

このGzwの大きさの尺度としてHノルムを使い、これがH基準となる。その定義は、次のとおりである。

Figure 2012113676
The H∞ norm is used as a measure of the magnitude of G zw , and this is the H∞ standard. The definition is as follows.
Figure 2012113676

ここで

Figure 2012113676

は、最大特異値をあらわし、ωは角周波数をあらわす。この定義により、H制御問題は、一般化プラントGに対し、u = Kyのフィードバック制御により、閉ループ系を内部安定化し、かつ、与えられた正数γに大して、
Figure 2012113676

を満たす制御器Kを求めることと定式化できる。γは任意に与えてもよいが、自動決定する方法も従来技術として存在する。これを解くための一つのアプローチが、LMIによるものである。 here
Figure 2012113676

Represents the maximum singular value, and ω represents the angular frequency. By this definition, the H control problem is to stabilize the closed loop system internally by feedback control of u = Ky for the generalized plant G, and to increase the given positive number γ,
Figure 2012113676

It is possible to formulate a controller K that satisfies the above. γ may be arbitrarily given, but there is a conventional method for automatically determining γ. One approach to solve this is by LMI.

すなわち、zとwの部分だけを取り出して

Figure 2012113676

としたとき、ある定理によれば、次の2つの条件は等価となる。
(i) 行列Aは安定で、図5に示すブロックにおいて、
Figure 2012113676

が成立する。
(ii) 次式を満たすX > 0が存在する。
Figure 2012113676
That is, take out only the parts z and w
Figure 2012113676

Then, according to a theorem, the following two conditions are equivalent.
(i) The matrix A is stable and in the block shown in FIG.
Figure 2012113676

Is established.
(ii) There exists X> 0 that satisfies the following formula.
Figure 2012113676

つまり、条件(i)は、上記H基準であり、それを満たす行列Xが制御器であり、条件(ii)のXについての線形行列不等式(LMI)を満たすものを求めると、H基準の制御器が求められるという次第である。 That is, the condition (i) is the above-mentioned H criterion, and the matrix X that satisfies the condition is a controller, and when the one that satisfies the linear matrix inequality (LMI) for X in the condition (ii) is obtained, the H criterion The controller is required.

一般に、フィードバック制御の場合、そのまま代入していくと、Xについて線形にならないので、変数変換などして上記の形式にもっていくことになる。この技術的詳細に関しては、適宜S.Aubouet et al.,”H-infinity/LPV observer for an industrial semi-active suspension,” IEEE Proc. Conf. on Control Applications, pp.756-763, 2009も参照されたい。   Generally, in the case of feedback control, if it is substituted as it is, it will not become linear with respect to X, so it will be converted into the above form by variable conversion. For technical details, see also S. Aubouet et al., “H-infinity / LPV observer for an industrial semi-active suspension,” IEEE Proc. Conf. On Control Applications, pp. 756-763, 2009. I want.

しかし、全てのスケジューリング・パラメータについてwを入力することはできないので、メイン・プログラム202は、ステップ308で、LMI変換モジュール210によって、端点において上記LMIを用いたH基準による制御器を計算し、次にステップ310で、パラメータ選定時に計算したνギャップによって重み付けした双一次補間によって、一般化プラントへの制御器を構成する。 However, since w cannot be input for all scheduling parameters, the main program 202 calculates a controller based on the H∞ criterion using the LMI at the end point by the LMI conversion module 210 in step 308, Next, in step 310, a controller for the generalized plant is configured by bilinear interpolation weighted by the ν gap calculated at the time of parameter selection.

なお、従来技術では、一般化プラントへの制御器Kは、下記のように、アフィン変換で表現できると暗黙的に仮定されていた。

Figure 2012113676

ここで、K(i)は、i番目の制御器である。 In the prior art, it has been implicitly assumed that the controller K for the generalized plant can be expressed by affine transformation as follows.
Figure 2012113676

Here, K (i) is the i-th controller.

これを、次の式のような、パラメータ選定時に計算したνギャップによって重み付けした双一次補間に置き換えたのが、本発明の1つの特徴である。

Figure 2012113676
One feature of the present invention is that this is replaced by bilinear interpolation weighted by the ν gap calculated at the time of parameter selection as in the following equation.
Figure 2012113676

ここで、Nとして何個を選ぶのかで、実施例が分かれる。第1の実施例では、最初に与えられたパラメータのみが選ばれ、すなわち、そのパラメータについて、K(1)をminに固定されたときの制御器、K(2)を、maxに固定されたときの制御器とし、minとmidのそれぞれに固定したときの伝達関数間のνギャップをν1、midとmaxのそれぞれに固定したときの伝達関数間のνギャップをν2とし、p1とp2を下記の式で与えることによって、双一次補間により、制御器が構成される。

Figure 2012113676
Here, the embodiment is divided according to how many N are selected. In the first embodiment, only the first given parameter is selected, that is, for that parameter, the controller when K (1) is fixed to min and K (2) is fixed to max. the controller of time, [nu gap between the transfer function when the [nu gap between the transfer function of the time fixed to the respective min and mid is fixed to each of the [nu 1, mid and max and [nu 2, and p 1 The controller is constructed by bilinear interpolation by giving p 2 by the following equation.
Figure 2012113676

第2の実施例では、最初に与えられたパラメータと、その次のパラメータの2つが選ばれ、すなわち、それらのパラメータについて、{p1_min, p2_min}の端点で設計された制御器をK(1)、{p1_min, p2_max}の端点で設計された制御器をK(2)、{p1_max, p2_min}の端点で設計された制御器をK(3)、{p1_max, p2_max}の端点で設計された制御器をK(4)とし、p1_minとp2_minのそれぞれに固定したときの伝達関数間のνギャップをν1、p1_minとp2_maxのそれぞれに固定したときの伝達関数間のνギャップをν2、p1_maxとp2_minのそれぞれに固定したときの伝達関数間のνギャップをν3、p1_maxとp2_maxのそれぞれに固定したときの伝達関数間のνギャップをν4とし、p1、p2、p3、p4を下記の式で与えることによって、双一次補間により、制御器が構成される。

Figure 2012113676
In the second embodiment, two parameters, the first given parameter and the next parameter, are selected, that is, a controller designed at the end point of {p1_min, p2_min} for these parameters is represented by K (1) The controller designed at the end point of {p1_min, p2_max} is K (2) , the controller designed at the end point of {p1_max, p2_min} is designed at the end point of K (3) , {p1_max, p2_max} When the controller is K (4) , the ν gap between the transfer functions when fixed at p1_min and p2_min is ν 1 , and the ν gap between the transfer functions when fixed at p1_min and p2_max is ν 2 , p1_max And p2_min, the ν gap between the transfer functions when fixed to ν 3 , and the ν gap between the transfer functions when fixed to p1_max and p2_max, respectively, as ν 4, and p 1 , p 2 , p 3 , p The controller is constructed by bilinear interpolation by giving 4 by the following equation.
Figure 2012113676

本発明の特徴に従い、線形補間でなく、νギャップによる重み付けを用いることは、非ユークリッド幾何的測地距離を考慮することになり、これにより、フィードバック制御の誤差を小さくすることにつながる。   In accordance with the features of the present invention, using weighting by ν gaps rather than linear interpolation will take into account non-Euclidean geometric geodesic distances, thereby reducing feedback control errors.

なお、3つ以上のパラメータを選んで制御器を構成する実施例も考えられるが、あまり選ぶパラメータの数を増やすと、制御器が複雑になって機能しない場合も生じるので、妥当な数に留めることが望ましい。   An embodiment in which the controller is configured by selecting three or more parameters is also conceivable, but if the number of parameters to be selected is increased too much, the controller may become complicated and may not function. It is desirable.

以上本発明の実施例を説明してきたが、本発明は、特定のハードウェア、ソフトウェアなどのプラット・フォームに限定されず、任意のプラット・フォーム上で実現可能であることを、この分野の当業者は理解するであろう。   Although the embodiments of the present invention have been described above, the present invention is not limited to specific hardware and software platforms, and can be implemented on any platform. The merchant will understand.

102 システム・バス
104 CPU
106 主記憶
108 ハードティスク
108 ハードディスク・ドライブ
114 ディスプレイ
116 インターフェース・ボード
118 制御器
202 メイン・プログラム
204 LPVモデル
206 ギャップ計算モジュール
206 スケジューリング・パラメータ
208 ギャップ計算モジュール
210 変換モジュール
212 制御器構成モジュール
402 プラント
404 一般化プラント
102 System bus 104 CPU
106 Main memory 108 Hard disk 108 Hard disk drive 114 Display 116 Interface board 118 Controller 202 Main program 204 LPV model 206 Gap calculation module 206 Scheduling parameter 208 Gap calculation module 210 Conversion module 212 Controller configuration module 402 Plant 404 Generalized plant

Claims (12)

コンピュータの処理により、スケジューリング・パラメータが与えられた、LPVモデルで記述されたプラントから、制御器を構成する方法であって、
制約条件の要求に基づく、前記LPVモデルのプラント・パラメータ・ベクトルの各候補のMin、Max及び中間値の各々について、νギャップを計算するステップと、
前記νギャップ間の差の大きさに基づき、前記プラント・パラメータ・ベクトルの候補を選定するステップと、
前記LPVモデルをLMIに変換し、前記選定したプラント・パラメータ・ベクトルの最大値と最小値に関する全ての組み合わせ端点についてH基準により、制御器を求めるステップと、
前記選定されたプラント・パラメータ・ベクトルにより計算された前記νギャップによって重み付けした双一次補間によって、前記プラントへの制御器を構成するステップを有する、
制御器の構成方法。
A method for constructing a controller from a plant described in an LPV model given scheduling parameters by computer processing, comprising:
Calculating a ν gap for each Min, Max and intermediate value of each candidate plant parameter vector of the LPV model based on constraint requirements;
Selecting the plant parameter vector candidates based on the magnitude of the difference between the ν gaps;
Converting the LPV model into an LMI and obtaining a controller according to the H criterion for all combination endpoints for the maximum and minimum values of the selected plant parameter vector;
Configuring a controller to the plant by bilinear interpolation weighted by the ν gap calculated by the selected plant parameter vector.
How to configure the controller.
前記νギャップを計算するステップが、スケジューリング・パラメータのi番目の要素の最大値を代入した場合の伝達関数をTi max、最小値を代入した場合の伝達関数をTi min、その中間値を代入した場合の伝達関数をTi midとし、δν()をνギャップを計算する関数であるとしたとき、
νi_mid_max = δν(Ti mid,Ti max)
νi_min_mid = δν(Ti min,Ti mid)
νi_min_max = δν(Ti min,Ti max)
を計算し、
前記プラント・パラメータ・ベクトルの候補を選定するステップが、
(|νi_min_mid - νi_mid_max| + νi_min_max)の値に基づき、前記プラント・パラメータ・ベクトルの候補を選定する、請求項1に記載の方法。
In the step of calculating the ν gap, the transfer function when the maximum value of the i-th element of the scheduling parameter is substituted is T i max , the transfer function when the minimum value is substituted is T i min , and the intermediate value is When the substitution transfer function is T i mid and δ ν () is a function that calculates the ν gap,
ν i _mid_max = δ ν (T i mid , T i max )
ν i _min_mid = δ ν (T i min , T i mid )
ν i _min_max = δ ν (T i min , T i max )
Calculate
Selecting the plant parameter vector candidates;
The method according to claim 1, wherein the plant parameter vector candidate is selected based on a value of (| ν i _min_mid−ν i _mid_max | + ν i _min_max).
前記プラント・パラメータ・ベクトルの候補を選定するステップが、
(|νi_min_mid - νi_mid_max| + νi_min_max)の値が最大の候補のみ1つを選定し、前記双一次補間は、選定されたパラメータのminとmaxを固定した2点で行われる、請求項2に記載の方法。
Selecting the plant parameter vector candidates;
Only one candidate with the maximum value of (| ν i _min_mid-ν i _mid_max | + ν i _min_max) is selected, and the bilinear interpolation is performed at two points with fixed min and max of the selected parameters. The method according to claim 2.
前記プラント・パラメータ・ベクトルの候補を選定するステップが、
(|νi_min_mid - νi_mid_max| + νi_min_max)の値が最大の候補と、その次に最大の候補の2つを選定し、前記双一次補間は、最大の候補をp1、次に最大の候補をp2としたとき、p1_minとp2_min、p1_minとp2_max、p1_maxとp2_min、p1_maxとp2_maxをそれぞれ固定した4点で行われる、請求項2に記載の方法。
Selecting the plant parameter vector candidates;
The candidate with the largest value of (| ν i _min_mid-ν i _mid_max | + ν i _min_max) and the next largest candidate are selected, and the bilinear interpolation uses p1 as the largest candidate, The method according to claim 2, wherein when the maximum candidate is p2, p1_min and p2_min, p1_min and p2_max, p1_max and p2_min, and p1_max and p2_max are fixed at four points, respectively.
コンピュータの処理により、スケジューリング・パラメータが与えられた、LPVモデルで記述されたプラントから、制御器を構成するプログラムであって、
前記コンピュータに、
制約条件の要求に基づく、前記LPVモデルのプラント・パラメータ・ベクトルの各候補のMin、Max及び中間値の各々について、νギャップを計算するステップと、
前記νギャップ間の差の大きさに基づき、前記プラント・パラメータ・ベクトルの候補を選定するステップと、
前記LPVモデルをLMIに変換し、前記選定したプラント・パラメータ・ベクトルの最大値と最小値に関する全ての組み合わせ端点についてH基準により、制御器を求めるステップと、
前記選定されたプラント・パラメータ・ベクトルにより計算された前記νギャップによって重み付けした双一次補間によって、前記プラントへの制御器を構成するステップを実行させる、
制御器の構成プログラム。
A program that configures a controller from a plant described in an LPV model given scheduling parameters by computer processing,
In the computer,
Calculating a ν gap for each Min, Max and intermediate value of each candidate plant parameter vector of the LPV model based on constraint requirements;
Selecting the plant parameter vector candidates based on the magnitude of the difference between the ν gaps;
Converting the LPV model into an LMI and obtaining a controller according to the H criterion for all combination endpoints for the maximum and minimum values of the selected plant parameter vector;
Causing the step of configuring a controller to the plant by bilinear interpolation weighted by the v gap calculated by the selected plant parameter vector;
Controller configuration program.
前記νギャップを計算するステップが、スケジューリング・パラメータのi番目の要素の最大値を代入した場合の伝達関数をTi max、最小値を代入した場合の伝達関数をTi min、その中間値を代入した場合の伝達関数をTi midとし、δν()をνギャップを計算する関数であるとしたとき、
νi_mid_max = δν(Ti mid,Ti max)
νi_min_mid = δν(Ti min,Ti mid)
νi_min_max = δν(Ti min,Ti max)
を計算し、
前記プラント・パラメータ・ベクトルの候補を選定するステップが、
(|νi_min_mid - νi_mid_max| + νi_min_max)の値に基づき、前記プラント・パラメータ・ベクトルの候補を選定する、請求項5に記載のプログラム。
In the step of calculating the ν gap, the transfer function when the maximum value of the i-th element of the scheduling parameter is substituted is T i max , the transfer function when the minimum value is substituted is T i min , and the intermediate value is When the substitution transfer function is T i mid and δ ν () is a function that calculates the ν gap,
ν i _mid_max = δ ν (T i mid , T i max )
ν i _min_mid = δ ν (T i min , T i mid )
ν i _min_max = δ ν (T i min , T i max )
Calculate
Selecting the plant parameter vector candidates;
6. The program according to claim 5, wherein the plant parameter vector candidate is selected based on a value of (| ν i _min_mid−ν i _mid_max | + ν i _min_max).
前記プラント・パラメータ・ベクトルの候補を選定するステップが、
(|νi_min_mid - νi_mid_max| + νi_min_max)の値が最大の候補のみ1つを選定し、前記双一次補間は、選定されたパラメータのminとmaxを固定した2点で行われる、請求項6に記載のプログラム。
Selecting the plant parameter vector candidates;
Only one candidate with the maximum value of (| ν i _min_mid-ν i _mid_max | + ν i _min_max) is selected, and the bilinear interpolation is performed at two points with fixed min and max of the selected parameters. The program according to claim 6.
前記プラント・パラメータ・ベクトルの候補を選定するステップが、
(|νi_min_mid - νi_mid_max| + νi_min_max)の値が最大の候補と、その次に最大の候補の2つを選定し、前記双一次補間は、最大の候補をp1、次に最大の候補をp2としたとき、p1_minとp2_min、p1_minとp2_max、p1_maxとp2_min、p1_maxとp2_maxをそれぞれ固定した4点で行われる、請求項6に記載のプログラム。
Selecting the plant parameter vector candidates;
The candidate with the largest value of (| ν i _min_mid-ν i _mid_max | + ν i _min_max) and the next largest candidate are selected, and the bilinear interpolation uses p1 as the largest candidate, The program according to claim 6, wherein the program is executed at four points where p1_min and p2_min, p1_min and p2_max, p1_max and p2_min, and p1_max and p2_max are fixed, respectively, when the maximum candidate is p2.
コンピュータの処理により、スケジューリング・パラメータが与えられた、LPVモデルで記述されたプラントから、制御器を構成するシステムであって、
記憶手段と、
前記記憶手段に格納された前記LPVモデル及びそのスケジューリング・パラメータと、
制約条件の要求に基づく、前記LPVモデルのプラント・パラメータ・ベクトルの各候補のMin、Max及び中間値の各々について、νギャップを計算する手段と、
前記νギャップ間の差の大きさに基づき、前記プラント・パラメータ・ベクトルの候補を選定する手段と、
前記LPVモデルをLMIに変換し、前記選定したプラント・パラメータ・ベクトルの最大値と最小値に関する全ての組み合わせ端点についてH基準により、制御器を求める手段と、
前記選定されたプラント・パラメータ・ベクトルにより計算された前記νギャップによって重み付けした双一次補間によって、前記プラントへの制御器を構成する手段を有する、
制御器の構成システム。
A system that configures a controller from a plant described in an LPV model given scheduling parameters by computer processing,
Storage means;
The LPV model and its scheduling parameters stored in the storage means;
Means for calculating a ν gap for each of the Min, Max and intermediate values of each candidate plant parameter vector of the LPV model, based on constraint requirements;
Means for selecting candidate plant parameter vectors based on the magnitude of the difference between the ν gaps;
Means for converting the LPV model into LMI and obtaining a controller according to H criteria for all combinational endpoints relating to the maximum and minimum values of the selected plant parameter vector;
Means for configuring a controller to the plant by bilinear interpolation weighted by the v gap calculated by the selected plant parameter vector;
Controller configuration system.
前記νギャップを計算する手段が、スケジューリング・パラメータのi番目の要素の最大値を代入した場合の伝達関数をTi max、最小値を代入した場合の伝達関数をTi min、その中間値を代入した場合の伝達関数をTi midとし、δν()をνギャップを計算する関数であるとしたとき、
νi_mid_max = δν(Ti mid,Ti max)
νi_min_mid = δν(Ti min,Ti mid)
νi_min_max = δν(Ti min,Ti max)
を計算し、
前記プラント・パラメータ・ベクトルの候補を選定する手段が、
(|νi_min_mid - νi_mid_max| + νi_min_max)の値に基づき、前記プラント・パラメータ・ベクトルの候補を選定する、請求項9に記載の制御器の構成システム。
The means for calculating the ν gap is T i max for the transfer function when the maximum value of the i-th element of the scheduling parameter is substituted, T i min for the transfer function when the minimum value is substituted, and the intermediate value thereof When the substitution transfer function is T i mid and δ ν () is a function that calculates the ν gap,
ν i _mid_max = δ ν (T i mid , T i max )
ν i _min_mid = δ ν (T i min , T i mid )
ν i _min_max = δ ν (T i min , T i max )
Calculate
Means for selecting candidate plant parameter vectors
The controller configuration system according to claim 9, wherein candidates for the plant parameter vector are selected based on a value of (| ν i _min_mid-ν i _mid_max | + ν i _min_max).
前記プラント・パラメータ・ベクトルの候補を選定する手段が、
(|νi_min_mid - νi_mid_max| + νi_min_max)の値が最大の候補のみ1つを選定し、前記双一次補間は、選定されたパラメータのminとmaxを固定した2点で行われる、請求項10に記載の制御器の構成システム。
Means for selecting candidate plant parameter vectors
Only one candidate with the maximum value of (| ν i _min_mid-ν i _mid_max | + ν i _min_max) is selected, and the bilinear interpolation is performed at two points with fixed min and max of the selected parameters. The system for configuring a controller according to claim 10.
前記プラント・パラメータ・ベクトルの候補を選定する手段が、
(|νi_min_mid - νi_mid_max| + νi_min_max)の値が最大の候補と、その次に最大の候補の2つを選定し、前記双一次補間は、最大の候補をp1、次に最大の候補をp2としたとき、p1_minとp2_min、p1_minとp2_max、p1_maxとp2_min、p1_maxとp2_maxをそれぞれ固定した4点で行われる、請求項10に記載の制御器の構成システム。
Means for selecting candidate plant parameter vectors
The candidate with the largest value of (| ν i _min_mid-ν i _mid_max | + ν i _min_max) and the next largest candidate are selected, and the bilinear interpolation uses p1 as the largest candidate, The controller configuration system according to claim 10, wherein p1_min and p2_min, p1_min and p2_max, p1_max and p2_min, and p1_max and p2_max are fixed at four points, respectively, when the maximum candidate is p2.
JP2010264660A 2010-11-29 2010-11-29 Controller configuration method, system, and program Expired - Fee Related JP5713338B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010264660A JP5713338B2 (en) 2010-11-29 2010-11-29 Controller configuration method, system, and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010264660A JP5713338B2 (en) 2010-11-29 2010-11-29 Controller configuration method, system, and program

Publications (2)

Publication Number Publication Date
JP2012113676A true JP2012113676A (en) 2012-06-14
JP5713338B2 JP5713338B2 (en) 2015-05-07

Family

ID=46497784

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010264660A Expired - Fee Related JP5713338B2 (en) 2010-11-29 2010-11-29 Controller configuration method, system, and program

Country Status (1)

Country Link
JP (1) JP5713338B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018008391A1 (en) * 2016-07-07 2018-01-11 日本電気株式会社 Linear parameter varying model estimation system, method, and program
WO2018016628A1 (en) * 2016-07-22 2018-01-25 株式会社明電舎 Shaft torque control device
CN110501294A (en) * 2019-08-07 2019-11-26 西安文理学院 A kind of multivariate calibration methods based on information fusion
WO2020070963A1 (en) * 2018-10-02 2020-04-09 株式会社明電舎 Control apparatus design method, control apparatus, and axial torque control apparatus
CN112631553A (en) * 2020-12-07 2021-04-09 北京无线电测量研究所 Radar servo system model establishing method, computer equipment and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10220475A (en) * 1997-02-06 1998-08-21 Seiko Seiki Co Ltd Magnetic bearing device using lmi-based gain schedule control
JPH11110003A (en) * 1997-10-02 1999-04-23 Mitsubishi Heavy Ind Ltd Design method for vehicle control system
JP2003345403A (en) * 2002-05-28 2003-12-05 Osaka Industrial Promotion Organization Method and device for designing, method, device and system for controlling, computer program and recording medium
JP2005354774A (en) * 2004-06-09 2005-12-22 Nissan Motor Co Ltd Drive force control device
JP2008045687A (en) * 2006-08-18 2008-02-28 Jtekt Corp Magnetic bearing device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10220475A (en) * 1997-02-06 1998-08-21 Seiko Seiki Co Ltd Magnetic bearing device using lmi-based gain schedule control
JPH11110003A (en) * 1997-10-02 1999-04-23 Mitsubishi Heavy Ind Ltd Design method for vehicle control system
JP2003345403A (en) * 2002-05-28 2003-12-05 Osaka Industrial Promotion Organization Method and device for designing, method, device and system for controlling, computer program and recording medium
JP2005354774A (en) * 2004-06-09 2005-12-22 Nissan Motor Co Ltd Drive force control device
JP2008045687A (en) * 2006-08-18 2008-02-28 Jtekt Corp Magnetic bearing device

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
G.VINNICOMBE: "Frequency Domain Uncertainty and the Graph Topology", IEEE TRANSACTIONS ON AUTOMATIC CONTROL, vol. 38, no. 9, JPN6014027824, September 1993 (1993-09-01), pages 1371 - 1383, ISSN: 0003009296 *
P.APKARIAN ET AL: "Self-scheduled H∞ Control of Linear Parameter-varying Systems: a Design Example", AUTOMATICA, vol. 31, no. 9, JPN6014027822, 1995, pages 1251 - 1261, ISSN: 0003009295 *
S.AUBOUET ET AL: "H∞/LPV observer for an industrial semi-active suspension", IEEE PROC. CONF. ON CONTROL APPLICATIONS, JPN6014027827, 2009, pages 756 - 763, XP031542885, ISSN: 0003009297 *
藤森 篤 他: "ゲインスケジューリングを用いた無人宇宙往還機の飛行制御系設計(v-gapを用いた動作点選択と局所多目的ゲイ", 日本航空宇宙学会論文集, vol. 51, no. 594, JPN6014027820, July 2003 (2003-07-01), JP, pages 339 - 345, ISSN: 0003009294 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2018008391A1 (en) * 2016-07-07 2019-04-25 日本電気株式会社 Linear parameter variation model estimation system, method and program
WO2018008391A1 (en) * 2016-07-07 2018-01-11 日本電気株式会社 Linear parameter varying model estimation system, method, and program
US10558179B2 (en) 2016-07-22 2020-02-11 National University Corporation Chiba University Shaft torque control device
WO2018016628A1 (en) * 2016-07-22 2018-01-25 株式会社明電舎 Shaft torque control device
JPWO2018016628A1 (en) * 2016-07-22 2019-02-28 株式会社明電舎 Shaft torque control device
WO2020070963A1 (en) * 2018-10-02 2020-04-09 株式会社明電舎 Control apparatus design method, control apparatus, and axial torque control apparatus
JP2020057183A (en) * 2018-10-02 2020-04-09 株式会社明電舎 Controller designing method, controller and axial torque controller
JP7046774B2 (en) 2018-10-02 2022-04-04 株式会社明電舎 Controller design method
US11422064B2 (en) 2018-10-02 2022-08-23 Meidensha Corporation Control apparatus design method, control apparatus, and axial torque control apparatus
CN110501294A (en) * 2019-08-07 2019-11-26 西安文理学院 A kind of multivariate calibration methods based on information fusion
CN110501294B (en) * 2019-08-07 2021-09-28 西安文理学院 Multivariate correction method based on information fusion
CN112631553A (en) * 2020-12-07 2021-04-09 北京无线电测量研究所 Radar servo system model establishing method, computer equipment and storage medium
CN112631553B (en) * 2020-12-07 2023-11-28 北京无线电测量研究所 Radar servo system model building method, computer equipment and storage medium

Also Published As

Publication number Publication date
JP5713338B2 (en) 2015-05-07

Similar Documents

Publication Publication Date Title
JP7171192B2 (en) Interactive system for modulating feedback loops
JP5713338B2 (en) Controller configuration method, system, and program
JP5723582B2 (en) Method, program and system for configuring controller in sliding mode control system
JP2008521138A (en) Method for analysis of control systems
CN109270946A (en) The attitude control method of flexible spacecraft
JP7090734B2 (en) Control system, control method and storage medium
JP5255924B2 (en) Power system reduction model creation device, power system reduction model creation method, and power system reduction model creation program
JP2019505889A (en) Cost function design system, cost function design method, and cost function design program
JP4815391B2 (en) Model parameter estimation calculation apparatus and method, model parameter estimation calculation processing program, and recording medium recording the same
CN115755629A (en) Robust controller design method and device suitable for active vibration isolation and electronic equipment
Zou et al. Finite-time adaptive interval type-2 fuzzy tracking control for Mecanum-wheel mobile robots
da Silva et al. Computer-aided integrated design for machines with varying dynamics
CN114310874B (en) Cooperative robot joint moment control method, system and computer equipment
JP6984597B2 (en) Linear parameter variation model estimation system, method and program
JP4612989B2 (en) Mechatronics system simulation method
JP5851941B2 (en) Control system design support apparatus, method, and program
CN115146416A (en) Engineering machinery hydraulic system model parameter identification method and device and excavator
CN114147710B (en) Robot control method and device, robot and storage medium
JP5888013B2 (en) Neural network design method, program, and digital analog fitting method
JP2020035325A (en) Design system, learned model generation method, and design program
US9329584B2 (en) Method, system and program for constructing a controller
WO2013151633A1 (en) Thermal deflection analysis
CN113937783A (en) Wind power plant reactive voltage coordination control method and device based on model prediction control
JPH10323070A (en) Motor controller
CN117331311B (en) Robot dynamics parameter estimation method based on acceleration-free recursive filtering regression

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130903

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20140604

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140708

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20140828

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140828

RD12 Notification of acceptance of power of sub attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7432

Effective date: 20140828

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20140829

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150123

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20150126

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150126

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

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20150220

RD14 Notification of resignation of power of sub attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7434

Effective date: 20150220

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20150304

R150 Certificate of patent or registration of utility model

Ref document number: 5713338

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees