JP6838454B2 - Control method, control device and control program - Google Patents

Control method, control device and control program Download PDF

Info

Publication number
JP6838454B2
JP6838454B2 JP2017058723A JP2017058723A JP6838454B2 JP 6838454 B2 JP6838454 B2 JP 6838454B2 JP 2017058723 A JP2017058723 A JP 2017058723A JP 2017058723 A JP2017058723 A JP 2017058723A JP 6838454 B2 JP6838454 B2 JP 6838454B2
Authority
JP
Japan
Prior art keywords
equation
matrix
lmi
linear
control
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.)
Active
Application number
JP2017058723A
Other languages
Japanese (ja)
Other versions
JP2018163395A (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.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2017058723A priority Critical patent/JP6838454B2/en
Publication of JP2018163395A publication Critical patent/JP2018163395A/en
Application granted granted Critical
Publication of JP6838454B2 publication Critical patent/JP6838454B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Description

本発明は、フィードバック制御技術に関する。 The present invention relates to a feedback control technique.

動的システムとは、出力が現在の入力だけでなく過去の入力にも依存するシステムのことである。特に以下のような、離散時間の線形状態方程式(Liner State Equation)(1)および線形出力方程式(Linear Output Equation)(2)で表される動的システムは、離散時間線形時不変システムと呼ばれる。 A dynamic system is one in which the output depends not only on the current input but also on the past input. In particular, the dynamic system represented by the Linear State Equation (1) and the Linear Output Equation (2) of discrete time as follows is called a linear time-invariant system of discrete time.

Figure 0006838454
Figure 0006838454

Figure 0006838454
Figure 0006838454

図1に、上記の離散時間線形時不変システムのブロック線図を示す。Ax[k]とBu[k]との加算によりx[k+1]が求められる。また、y[k]がCx[k]として求められる。q-1はx[k]のkを1減らす作用素である。 FIG. 1 shows a block diagram of the above-mentioned discrete-time linear time-invariant system. X [k + 1] can be obtained by adding Ax [k] and Bu [k]. Further, y [k] is obtained as Cx [k]. q -1 is an operator that reduces k of x [k] by 1.

但し、動的システムの状態を必ずしも直接観測できるとは限らず、直接観測できない場合には観測可能な出力から状態が推定される。例えばフィードバック制御においては、観測可能な出力に基づき観測或いは制御の対象(以下、制御対象と呼ぶ)の内部状態が推定され、推定された内部状態に応じて入力が決定される。 However, it is not always possible to directly observe the state of the dynamic system, and if it cannot be observed directly, the state is estimated from the observable output. For example, in feedback control, the internal state of an observation or control target (hereinafter referred to as a control target) is estimated based on the observable output, and the input is determined according to the estimated internal state.

可検出である線形システムについて、時間の経過とともに収束する状態の値を計算する装置はオブザーバと呼ばれる。オブザーバによる推定値は、十分に時間が経過すれば真の値に収束することが保証されている。従来の技術においては、推定値が真の値に収束したと見做すことができる十分な時間が経過した後の推定値が制御に利用されている。 For a detectable linear system, a device that calculates the value of a state that converges over time is called an observer. Observer estimates are guaranteed to converge to true values over time. In the conventional technique, the estimated value after a sufficient time that can be regarded as having converged to the true value is used for control.

しかし、制御対象が安定ではない場合、推定値が真の値に収束することを待っていると、内部状態が変化し制御が間に合わないことがある。図2及び図3に、時間と状態との関係を示す。図2において、縦軸は状態x1の値を表し、横軸は時間を表す。図3において、縦軸は状態x2の値を表し、横軸は時間を表す。真の値は実線でプロットされており、推定値は破線でプロットされている。図2及び図3のいずれにおいても、推定値が真の値に収束するまでに状態が変化しているので、収束後の推定値に基づく制御を実行する前に適切ではない状態が発生する可能性がある。一方で、推定値が真の値に収束する前においては、真の状態が満たす制約(たとえば、線形の等式関係)が満たされていないために、推定値に基づく制御を行うことができない場合がある。 However, when the controlled object is not stable and waiting for the estimated value to converge to the true value, the internal state may change and the control may not be in time. 2 and 3 show the relationship between time and state. In FIG. 2, the vertical axis represents the value of the state x 1 and the horizontal axis represents the time. In FIG. 3, the vertical axis represents the value of the state x 2 and the horizontal axis represents the time. True values are plotted with solid lines and estimates are plotted with dashed lines. In both FIGS. 2 and 3, the state changes before the estimated value converges to the true value, so an inappropriate state may occur before performing control based on the estimated value after convergence. There is sex. On the other hand, before the estimated value converges to the true value, the constraint that the true state satisfies (for example, the linear equation relation) is not satisfied, so that the control based on the estimated value cannot be performed. There is.

また、例えば射影操作により推定値を真の状態と同じ制約を満たすように変換して制御に利用する場合、推定値が得られる度に変換を行う必要があるため、演算時間の増加を招く。従って、安定ではない状況(例えば、制御に利用するための値を生成するための時間的な余裕が無い状況)においては、不正確な値に基づく制御を行うことになり、結果として不安定な制御につながる可能性がある。図4に、演算時間が増加した場合における状態x1の値と時間との関係を示し、図5に、演算時間が増加した場合における状態x2の値と時間との関係を示す。図4及び図5のいずれにおいても、真の値が振動的に変動しており、制御が不安定になっている。 Further, for example, when the estimated value is converted so as to satisfy the same constraint as the true state by a projection operation and used for control, it is necessary to perform the conversion every time the estimated value is obtained, which causes an increase in calculation time. Therefore, in an unstable situation (for example, a situation where there is no time to generate a value to be used for control), control based on an inaccurate value is performed, and as a result, it is unstable. May lead to control. FIG. 4 shows the relationship between the value of the state x 1 and the time when the calculation time increases, and FIG. 5 shows the relationship between the value of the state x 2 and the time when the calculation time increases. In both FIGS. 4 and 5, the true value fluctuates oscillatingly, and the control becomes unstable.

特開平10−54279号公報Japanese Unexamined Patent Publication No. 10-54279

R.F.Stengel, "Optimal Control and Estimation", Dover edition, Dover Publications, 1994.R.F.Stengel, "Optimal Control and Optimization", Dover edition, Dover Publications, 1994. Gene F. Franklin, J. David Powell and Michael Workman, "Digital Control of Dynamic Systems", third edition, Ellis-Kagle Press, 1998.Gene F. Franklin, J. David Powell and Michael Workman, "Digital Control of Dynamic Systems", third edition, Ellis-Kagle Press, 1998. Mustapha Ait Rami, Fernand Tadeo and Uwe Helmke, "Positive observers for linear positive systems, and their implications", International Journal of Control, vol.84, no.4, pp.716-725, 2011.Mustapha Ait Rami, Fernand Tadeo and Uwe Helmke, "Positive observers for linear positive systems, and their implications", International Journal of Control, vol.84, no.4, pp.716-725, 2011.

本発明の目的は、1つの側面では、必ずしも真の値に収束しているとは限らない推定値を用いて、制御対象の制御を適切に行う技術を提供することである。 An object of the present invention is to provide a technique for appropriately controlling a controlled object by using an estimated value that does not necessarily converge to a true value in one aspect.

一態様に係る制御方法は、制御対象の内部状態を表す変数が満たす条件に対応する方程式に基づき、制御対象の内部状態の推定に用いられるゲイン行列を生成し、生成されたゲイン行列を用いるオブザーバにより生成された、制御対象の内部状態の推定結果に基づき、制御対象を制御する処理を含む。 The control method according to one aspect generates a gain matrix used for estimating the internal state of the controlled object based on an equation corresponding to a condition satisfied by a variable representing the internal state of the controlled object, and an observer using the generated gain matrix. Includes processing to control the control target based on the estimation result of the internal state of the control target generated by.

1つの側面では、必ずしも真の値に収束しているとは限らない推定値を用いて、制御対象の制御を適切に行えるようになる。 In one aspect, it is possible to appropriately control the controlled object by using an estimated value that does not necessarily converge to the true value.

図1は、離散時間線形時不変システムのブロック線図である。FIG. 1 is a block diagram of a discrete-time linear time-invariant system. 図2は、時間と状態との関係を示す図である。FIG. 2 is a diagram showing the relationship between time and state. 図3は、時間と状態との関係を示す図である。FIG. 3 is a diagram showing the relationship between time and state. 図4は、時間と状態との関係を示す図である。FIG. 4 is a diagram showing the relationship between time and state. 図5は、時間と状態との関係を示す図である。FIG. 5 is a diagram showing the relationship between time and state. 図6は、状態空間における平面と値とを示す図である。FIG. 6 is a diagram showing a plane and a value in the state space. 図7は、オブザーバのブロック線図である。FIG. 7 is a block diagram of the observer. 図8は、状態空間における平面と値とを示す図である。FIG. 8 is a diagram showing a plane and a value in the state space. 図9は、第1の実施の形態の動的システムの概要を示す図である。FIG. 9 is a diagram showing an outline of the dynamic system of the first embodiment. 図10は、制御装置の機能ブロック図である。FIG. 10 is a functional block diagram of the control device. 図11は、制御対象装置の機能ブロック図である。FIG. 11 is a functional block diagram of the controlled device. 図12は、設計部が実行する処理の処理フローを示す図である。FIG. 12 is a diagram showing a processing flow of processing executed by the design unit. 図13は、状態推定部が実行する処理の処理フローを示す図である。FIG. 13 is a diagram showing a processing flow of processing executed by the state estimation unit. 図14は、通信部が実行する処理の処理フローを示す図である。FIG. 14 is a diagram showing a processing flow of processing executed by the communication unit. 図15は、通信部が実行する処理の処理フローを示す図である。FIG. 15 is a diagram showing a processing flow of processing executed by the communication unit. 図16は、第2の実施の形態の動的システムの概要を示す図である。FIG. 16 is a diagram showing an outline of the dynamic system of the second embodiment. 図17は、コンピュータの機能ブロック図である。FIG. 17 is a functional block diagram of the computer.

[実施の形態1]
本実施の形態においては、可検出であり且つ状態が以下のような線形等式(3)を満たす離散時間線形時不変システムを制御対象とする。
[Embodiment 1]
In the present embodiment, the control target is a discrete-time linear time-invariant system that is detectable and whose state satisfies the linear equation (3) as follows.

Figure 0006838454
Figure 0006838454

Figure 0006838454
Figure 0006838454

線形状態方程式(Liner State Equation)(1)および線形出力方程式(Linear Output Equation)(2)は上に示したとおりであり、A、B、C、Wおよびυは既知であるとする。線形等式(3)は、図6に示すように、状態の値が状態空間内の平面上に有る(すなわち、保存則が成立している)ことを意味する。 It is assumed that the Linear State Equation (1) and the Linear Output Equation (2) are as shown above, and A, B, C, W and υ are known. The linear equation (3) means that the value of the state is on a plane in the state space (that is, the conservation law holds), as shown in FIG.

一般に、オブザーバによる推定値^x[k+1]は以下のように表される。 In general, the observer's estimate ^ x [k + 1] is expressed as follows.

Figure 0006838454
Figure 0006838454

オブザーバのブロック線図は図7に示すとおりである。 The block diagram of the observer is as shown in FIG.

行列Lはゲイン行列(オブザーバゲインとも呼ばれる)である。ゲイン行列は、例えばLMI(Linear Matrix Inequality)可解問題を解くことで求められる。LMI可解問題は効率的な数値計算手法が知られた問題であり、ソルバ(例えばSDP(SemiDefinite Program)ソルバ)も存在するので、存在判定法及び構成法として適切である。 The matrix L is a gain matrix (also called an observer gain). The gain matrix can be obtained by solving, for example, an LMI (Linear Matrix Inequality) solvable problem. The LMI solvable problem is a problem for which an efficient numerical calculation method is known, and since a solver (for example, an SDP (SemiDefinite Program) solver) also exists, it is appropriate as an existence determination method and a construction method.

まず、入力されたA、B及びCに基づき以下のLMI可解問題を解くことで、行列X及び行列Yが求められる。 First, the matrix X and the matrix Y are obtained by solving the following LMI solvable problem based on the input A, B, and C.

Figure 0006838454
Figure 0006838454

この式は、左辺の行列が正定であることを表す。そして、ゲイン行列はL=X-1Yとして計算される。これにより推定値^x[k+1]を求めることができる。 This equation indicates that the matrix on the left side is positive. Then, the gain matrix is calculated as L = X -1 Y. This makes it possible to obtain the estimated value ^ x [k + 1].

但し、背景技術の欄で述べたように、このような方法で求めた推定値は等式制約を満たさないことがある。このことは、図8に示すように、推定値が状態空間内の平面上に無い場合があることを意味する。また、推定値が状態空間内の平面上に存在するように変換を行うと、演算時間が長くなり制御が不安定になることがある。 However, as described in the background technology section, the estimated values obtained by such a method may not satisfy the equation constraint. This means that, as shown in FIG. 8, the estimated value may not be on a plane in the state space. Further, if the conversion is performed so that the estimated value exists on a plane in the state space, the calculation time becomes long and the control may become unstable.

そこで以下では、通常のオブザーバと比較して計算量を増加させることなく(条件A−1)、真の値と同じ線形等式を満たしつつ(条件A−2)且つ時間の経過とともに真の値に収束する推定値を計算するオブザーバが存在する(条件A−3)か判定する方法を検討する。そして、そのようなオブザーバが存在する場合にはそのようなオブザーバを具体的に実現する(条件B)ことが可能であり且つ効率的な数値計算手法に基づいた(条件C)方法を検討する。 Therefore, in the following, the true value is satisfied with the passage of time while satisfying the same linear equation as the true value (condition A-1) without increasing the amount of calculation as compared with a normal observer (condition A-1). Consider a method for determining whether there is an observer (condition A-3) that calculates an estimated value that converges on. Then, when such an observer exists, a method (condition C) based on an efficient numerical calculation method capable of concretely realizing such an observer (condition B) will be examined.

条件A−1に関して、通常のオブザーバと比較して計算量が増加しないための必要十分条件は、計算式が通常のオブザーバと同じ形式であることである。すなわち、本実施の形態の計算式が式(5)の形式であることである。 Regarding condition A-1, the necessary and sufficient condition for the calculation amount not to increase as compared with the normal observer is that the calculation formula has the same format as the normal observer. That is, the calculation formula of the present embodiment is in the form of the formula (5).

条件A−2に関して、推定値が真の値と同じ線形等式を満たすための必要十分条件は、W^x[0]=υかつWL=0である。証明は以下のとおりである。まず、以下が成立する。 For condition A-2, the necessary and sufficient conditions for the estimated value to satisfy the same linear equation as the true value are W ^ x [0] = υ and WL = 0. The proof is as follows. First, the following holds.

Figure 0006838454
Figure 0006838454

(十分性)
WL=0であるとき、以下が成立する。
(Sufficiency)
When WL = 0, the following holds.

Figure 0006838454
Figure 0006838454

つまり、W^x[k]はWx[k]と同じ初期値を有し且つWx[k]と同じ漸化式に従うので、以下が成立する。 That is, since W ^ x [k] has the same initial value as Wx [k] and follows the same recurrence formula as Wx [k], the following holds.

Figure 0006838454
Figure 0006838454

(必要性)
逆にWL≠0であるときは、以下の式が成立しない場合が存在する。
(necessity)
On the contrary, when WL ≠ 0, the following equation may not hold.

Figure 0006838454
Figure 0006838454

従って、W^x[k]はWx[k]と異なる漸化式に従うので、式(9)が成立しない。以上で条件A−2についての必要十分条件の証明を終了する。 Therefore, since W ^ x [k] follows a recurrence formula different from Wx [k], the formula (9) does not hold. This completes the proof of the necessary and sufficient conditions for condition A-2.

条件A−3に関して、以下を誤差系と呼ぶ。 Regarding condition A-3, the following is called an error system.

Figure 0006838454
Figure 0006838454

時間の経過とともに推定値が真の値に収束することを「誤差系が安定である」という。誤差系を安定にするような推定値の計算方法が存在する(すなわち、そのようなゲイン行列が存在する)ための必要十分条件(すなわち、条件A−3の必要十分条件)は、式(6)のLMI可解問題の解が存在することである。 The fact that the estimated value converges to the true value with the passage of time is called "the error system is stable". The necessary and sufficient condition (that is, the necessary and sufficient condition of condition A-3) for the existence of the calculation method of the estimated value that stabilizes the error system (that is, the existence of such a gain matrix) is given by the equation (6). ) Is that there is a solution to the LMI solvable problem.

条件Bに関して、誤差系を安定にするオブザーバを実現するためには、式(5)のゲイン行列をL=X-1Yとして求めればよい。 Regarding the condition B, in order to realize an observer that stabilizes the error system, the gain matrix of the equation (5) may be obtained as L = X -1 Y.

以上から、条件A−3の必要十分条件である式(6)と条件A−2の必要十分条件であるWL=0とを連立させ、さらにL=X-1Yを代入した以下の問題を解くことができれば、条件A−1、条件A−2、条件A−3および条件Bが満たされる。 From the above, the following problem is solved by combining equation (6), which is a necessary and sufficient condition of condition A-3, with WL = 0, which is a necessary and sufficient condition of condition A-2, and substituting L = X -1 Y. If it can be solved, the conditions A-1, A-2, A-3 and B are satisfied.

Figure 0006838454
Figure 0006838454

しかし、式(12)は行列を変数とする非線形問題であるので、効率的な数値計算手法を適用することができず、条件Cが満たされない。 However, since Eq. (12) is a non-linear problem with a matrix as a variable, an efficient numerical calculation method cannot be applied, and the condition C is not satisfied.

そこで本実施の形態においては、推定値が線形等式を満たすための必要条件をゲイン行列に関する等式で表し、この等式を用いて条件A−1、条件A−2、条件A−3及び条件BをLMI可解問題に帰着させる。 Therefore, in the present embodiment, the necessary conditions for the estimated value to satisfy the linear equation are expressed by an equation related to the gain matrix, and these equations are used to express conditions A-1, condition A-2, condition A-3, and conditions A-3. Condition B is reduced to the LMI solvable problem.

まず、線形等式の係数行列Wを以下のように特異値分解(Singular Value Decomposition)する。U及びVは直交行列である。なお、一般の等式(f(x[k])=0)ではなく線形等式に限定しているため特異値分解を実行することができる。 First, the coefficient matrix W of the linear equation is singular value decomposed (Singular Value Decomposition) as follows. U and V are orthogonal matrices. It should be noted that the singular value decomposition can be executed because it is limited to the linear equation instead of the general equation (f (x [k]) = 0).

Figure 0006838454
Figure 0006838454

以下の条件(a)及び条件(b)を満たすゲイン行列Lが存在するか判定し、存在した場合にはゲイン行列Lを計算する。 It is determined whether or not a gain matrix L that satisfies the following conditions (a) and (b) exists, and if so, the gain matrix L is calculated.

条件(a)は、ゲイン行列LがΛ∈R(n-r)×pと直交行列Vとを用いて以下のように表されることである。 The condition (a) is that the gain matrix L is expressed as follows using Λ ∈ R (nr) × p and the orthogonal matrix V.

Figure 0006838454
Figure 0006838454

なお、条件(a)はWL=0であるための必要十分条件である。この証明については付録を参照されたい。このように、条件(A−2)の必要十分条件を式(14)の形で表現することで、LMI可解問題に帰着させることが可能になる。 The condition (a) is a necessary and sufficient condition for WL = 0. See the appendix for this proof. In this way, by expressing the necessary and sufficient condition of the condition (A-2) in the form of the equation (14), it becomes possible to reduce it to the LMI solvable problem.

条件(b)は、A−LCの全ての固有値複素平面における原点を中心とする半径1の円の内部に存在することである。 The condition (b) is that it exists inside a circle having a radius of 1 centered on the origin in all the eigenvalue complex planes of A-LC.

条件(a)及び条件(b)を満たすゲイン行列Lが存在する場合には、ゲイン行列Lを計算し、計算したゲイン行列Lに基づく式(5)に従った状態推定を行うようにオブザーバを設定する。以下では、このように設定されたオブザーバを線形等式保存オブザーバと呼ぶ。 When the gain matrix L satisfying the condition (a) and the condition (b) exists, the observer calculates the gain matrix L and estimates the state according to the equation (5) based on the calculated gain matrix L. Set. In the following, the observer set in this way is referred to as a linear equation preservation observer.

条件(a)及び条件(b)を満たすゲイン行列Lが存在するか判定し且つゲイン行列Lが存在する場合にゲイン行列Lを計算する手順は、以下のとおりである。 The procedure for determining whether or not the gain matrix L satisfying the condition (a) and the condition (b) exists and calculating the gain matrix L when the gain matrix L exists is as follows.

まず、以下のLMI可解問題を解く。 First, the following LMI solvable problem is solved.

Figure 0006838454
Figure 0006838454

式(15)が可解である場合(すなわち、解X及びYが存在する場合)、条件(a)及び条件(b)を満たすゲイン行列Lが存在すると判定される。式(15)が可解ではない場合、条件(a)及び条件(b)を満たすゲイン行列Lが存在しないと判定される。式(15)が可解である場合、以下の式に基づきゲイン行列Lが計算される。 When the equation (15) is solvable (that is, when the solutions X and Y exist), it is determined that the gain matrix L satisfying the condition (a) and the condition (b) exists. If the equation (15) is not solvable, it is determined that there is no gain matrix L that satisfies the condition (a) and the condition (b). If equation (15) is solvable, the gain matrix L is calculated based on the following equation.

Figure 0006838454
Figure 0006838454

式(15)によってゲイン行列Lの存在を判定できること及び式(16)によってゲイン行列を計算できることの証明を以下に示す。 The proof that the existence of the gain matrix L can be determined by the equation (15) and that the gain matrix can be calculated by the equation (16) is shown below.

式(11)で表された誤差系が安定である(n→∞,x[k]−^x[k]→0)ための必要十分条件は、A−LCの全ての固有値が複素平面の原点を中心とする半径1の円内に存在することであり、それは以下の行列不等式(すなわち、離散時間線形時不変システムに対するリアプノフ不等式)を満たす行列Pおよびゲイン行列Lが存在することと等価である。この点については、「蛯原義雄,「LMIによるシステム制御」,森北出版,2012」を参照されたい。 The necessary and sufficient condition for the error system represented by the equation (11) to be stable (n → ∞, x [k] − ^ x [k] → 0) is that all the eigenvalues of A-LC are in the complex plane. It exists in a circle with a radius of 1 centered on the origin, which is equivalent to the existence of a matrix P and a gain matrix L that satisfy the following matrix inequality (that is, the Riapnov inequality for a discrete-time linear time-invariant system). is there. For this point, refer to "Yoshio Ebihara," System Control by LMI ", Morikita Publishing, 2012".

Figure 0006838454
Figure 0006838454

WL=0と等価な条件である式(14)をこの行列不等式に代入すると、以下のようになる。 Substituting equation (14), which is a condition equivalent to WL = 0, into this matrix inequality gives the following.

Figure 0006838454
Figure 0006838454

この行列不等式を満たす行列P及び行列Λが存在することは、誤差系を安定とする式(14)の形式のゲイン行列Lが存在するための必要十分条件である。第1式にVVT=Iを挿入して整理し且つ第2式の両辺左からVTをかけ右からVをかけると、以下のようになる。 The existence of the matrix P and the matrix Λ that satisfy this matrix inequality is a necessary and sufficient condition for the existence of the gain matrix L of the form of the equation (14) that stabilizes the error system. When VV T = I is inserted into the first equation and arranged, and V T is applied from the left on both sides of the second equation and V is applied from the right, the result is as follows.

Figure 0006838454
Figure 0006838454

この式を満たす行列P及び行列Λが存在することは、誤差系を安定とする式(14)の形式のゲイン行列Lが存在するための必要十分条件である。 The existence of the matrix P and the matrix Λ that satisfy this equation is a necessary and sufficient condition for the existence of the gain matrix L of the form of the equation (14) that stabilizes the error system.

ここでX:=VTPVとおくと、以下のようになる。 Here, if X: = VT PV, it becomes as follows.

Figure 0006838454
Figure 0006838454

この式を満たす行列X及び行列Λが存在することは、誤差系を安定とする式(14)の形式のゲイン行列Lが存在するための必要十分条件である。 The existence of the matrix X and the matrix Λ that satisfy this equation is a necessary and sufficient condition for the existence of the gain matrix L of the form of the equation (14) that stabilizes the error system.

この式も離散時間線形時不変システムに対するリアプノフ不等式であり、以下の式と等価である。この点については、「小原敦美,「行列不等式アプローチによる制御系設計」,コロナ社,2016」を参照されたい。 This equation is also a Lyapunov inequality for a discrete-time linear time-invariant system and is equivalent to the following equation. For this point, refer to "Atsumi Ohara," Control System Design by Matrix Inequality Approach ", Corona Publishing Co., Ltd., 2016".

Figure 0006838454
Figure 0006838454

よって、この式を満たす行列X及び行列Λが存在することは、誤差系を安定とする式(14)の形式のゲイン行列Lが存在するための必要十分条件である。 Therefore, the existence of the matrix X and the matrix Λ that satisfy this equation is a necessary and sufficient condition for the existence of the gain matrix L of the form of the equation (14) that stabilizes the error system.

さらにY:=XΛとおくと以下の線形行列不等式を得ることができる。 Further, if Y: = XΛ, the following linear matrix inequality can be obtained.

Figure 0006838454
Figure 0006838454

よって、この線形行列不等式が解X>0及び解Yを持つときに限り誤差系を安定とする式(14)の形式のゲイン行列Lが存在する。そして、具体的なゲイン行列LはYの定義から式(16)のように求められる。 Therefore, there is a gain matrix L of the form of equation (14) that stabilizes the error system only when this linear matrix inequality has solution X> 0 and solution Y. Then, the specific gain matrix L can be obtained from the definition of Y as in Eq. (16).

以下では、以上の議論に基づく本実施の形態の具体的な態様を説明する。 Hereinafter, a specific embodiment of the present embodiment based on the above discussion will be described.

図9に、第1の実施の形態の動的システムの概要を示す。制御装置1及び制御装置1により制御される制御対象装置3は、例えばLAN(Local Area Network)又はインターネット等のネットワーク2に接続される。 FIG. 9 shows an outline of the dynamic system of the first embodiment. The control device 1 and the control target device 3 controlled by the control device 1 are connected to a network 2 such as a LAN (Local Area Network) or the Internet.

制御装置1は、例えば、パーソナルコンピュータ、サーバ、或いはマイクロコントローラである。図10に、制御装置1の機能ブロック図を示す。制御装置1は、設計部101と、オブザーバとしての処理を実行する状態推定部103と、通信部105と、出力データ格納部111と、推定値格納部113と、ワークデータ格納部115とを含む。 The control device 1 is, for example, a personal computer, a server, or a microcontroller. FIG. 10 shows a functional block diagram of the control device 1. The control device 1 includes a design unit 101, a state estimation unit 103 that executes processing as an observer, a communication unit 105, an output data storage unit 111, an estimated value storage unit 113, and a work data storage unit 115. ..

設計部101、状態推定部103及び通信部105は、プログラムが図17におけるメモリ2501にロードされCPU(Central Processing Unit)2503により実行されることで実現される。出力データ格納部111、推定値格納部113及びワークデータ格納部115は、図17におけるメモリ2501又はHDD(Hard Disk Drive)2505上に設けられる。 The design unit 101, the state estimation unit 103, and the communication unit 105 are realized by loading the program into the memory 2501 in FIG. 17 and executing it by the CPU (Central Processing Unit) 2503. The output data storage unit 111, the estimated value storage unit 113, and the work data storage unit 115 are provided on the memory 2501 or the HDD (Hard Disk Drive) 2505 in FIG.

設計部101は、ワークデータ格納部115に格納されているパラメータに基づきゲイン行列Lを計算し、計算したゲイン行列をワークデータ格納部115に格納する。状態推定部103は、ワークデータ格納部115に格納されているデータ及び出力データ格納部111に格納されているデータに基づき制御対象装置3の状態の推定値を計算し、計算した推定値を推定値格納部113に格納する。通信部105は、制御対象装置3の出力を制御対象装置3から受信し、受信した出力を出力データ格納部111に格納する。また、通信部105は、制御対象装置3への入力を制御対象装置3へ送信する。 The design unit 101 calculates the gain matrix L based on the parameters stored in the work data storage unit 115, and stores the calculated gain matrix in the work data storage unit 115. The state estimation unit 103 calculates an estimated value of the state of the controlled device 3 based on the data stored in the work data storage unit 115 and the data stored in the output data storage unit 111, and estimates the calculated estimated value. It is stored in the value storage unit 113. The communication unit 105 receives the output of the control target device 3 from the control target device 3, and stores the received output in the output data storage unit 111. Further, the communication unit 105 transmits the input to the control target device 3 to the control target device 3.

図11に、制御対象装置3の機能ブロック図を示す。制御対象装置3は、物理センサであるセンサ301と、被制御部302と、通信部303とを含む。 FIG. 11 shows a functional block diagram of the controlled device 3. The controlled target device 3 includes a sensor 301 which is a physical sensor, a controlled unit 302, and a communication unit 303.

センサ301は、観測可能な出力yを計測する。被制御部302は、制御装置1から受信した入力uに従って動作が制御される。通信部303は、制御装置1から入力uを受信する。また、通信部303は、センサ301により計測された出力yを制御装置1へ送信する。 The sensor 301 measures the observable output y. The operation of the controlled unit 302 is controlled according to the input u received from the control device 1. The communication unit 303 receives the input u from the control device 1. Further, the communication unit 303 transmits the output y measured by the sensor 301 to the control device 1.

例えば制御対象装置3が車両である場合、被制御部302は例えばモータであり、入力uはモータの回転数であり、出力yは例えば制御対象装置3の位置情報であり、状態xは例えば位置情報と制御対象装置3の速度情報との組合せである。また、例えば制御対象装置3が複数の貯水タンクである場合、被制御部302は例えば複数の貯水タンクの水量を調節するポンプであり、入力uは複数の貯水タンク間での水の輸送量であり、出力yは例えば複数の貯水タンクのうちいずれかの貯水タンクの水量であり、状態xは例えば複数の貯水タンクそれぞれの水量である。但し、実施形態はこのような例に限られるわけではない。例えば、電気回路、化学プラント、溶鉱炉、空調機、ロボット、サーバ、サーバルームおよびデータセンタ等の動的システムを対象としてもよい。 For example, when the controlled target device 3 is a vehicle, the controlled unit 302 is, for example, a motor, the input u is the rotation speed of the motor, the output y is, for example, the position information of the controlled target device 3, and the state x is, for example, the position. It is a combination of the information and the speed information of the control target device 3. Further, for example, when the controlled device 3 is a plurality of water storage tanks, the controlled unit 302 is, for example, a pump that adjusts the amount of water in the plurality of water storage tanks, and the input u is the amount of water transported between the plurality of water storage tanks. Yes, the output y is, for example, the amount of water in any of the plurality of water storage tanks, and the state x is, for example, the amount of water in each of the plurality of water storage tanks. However, the embodiment is not limited to such an example. For example, dynamic systems such as electric circuits, chemical plants, blast furnaces, air conditioners, robots, servers, server rooms and data centers may be targeted.

図12を用いて、制御装置1の設計部101が実行する処理を説明する。 The process executed by the design unit 101 of the control device 1 will be described with reference to FIG.

設計部101は、パラメータ(ここでは、式(15)の行列A、B、C及びW)の入力を受け付け(図12:ステップS1)、ワークデータ格納部115に格納する。なお、ステップS1においては、記憶領域に予め格納されたパラメータを設計部101が読み出すようにしてもよい。 The design unit 101 receives the input of the parameters (here, the matrices A, B, C and W of the equation (15)) (FIG. 12: step S1) and stores them in the work data storage unit 115. In step S1, the design unit 101 may read out the parameters stored in advance in the storage area.

設計部101は、ステップS1においてワークデータ格納部115に格納されたパラメータを用いて、式(15)で表されるLMI可解問題の解X及びYを求める処理を実行する(ステップS3)。 The design unit 101 executes a process of obtaining solutions X and Y of the LMI solvable problem represented by the equation (15) using the parameters stored in the work data storage unit 115 in step S1 (step S3).

設計部101は、ステップS3において解が存在するか判定する(ステップS5)。解が存在しない場合(ステップS5:Noルート)、設計部101は、以下の処理を実行する。具体的には、設計部101は、通常の方法で状態推定部103が推定値を計算するように設定を行う(ステップS7)。そして処理は終了する。通常の方法とは、例えば、通常のオブザーバで推定値を求めた後、等式制約を満たすように推定値を変換する処理である。ステップS7の処理は本実施の形態の主要な処理ではないので、詳細な説明を省略する。 The design unit 101 determines in step S3 whether a solution exists (step S5). When the solution does not exist (step S5: No route), the design unit 101 executes the following processing. Specifically, the design unit 101 sets the state estimation unit 103 to calculate the estimated value by a normal method (step S7). And the process ends. The usual method is, for example, a process of obtaining an estimated value by a normal observer and then converting the estimated value so as to satisfy the equation constraint. Since the process of step S7 is not the main process of the present embodiment, detailed description thereof will be omitted.

一方、解が存在する場合(ステップS5:Yesルート)、設計部101は、以下の処理を実行する。具体的には、設計部101は、求められた解X及びYを用いて、式(16)でゲイン行列Lを計算する(ステップS9)。 On the other hand, when a solution exists (step S5: Yes route), the design unit 101 executes the following processing. Specifically, the design unit 101 calculates the gain matrix L by the equation (16) using the obtained solutions X and Y (step S9).

設計部101は、ステップS9において計算されたゲイン行列Lを、ワークデータ格納部115における、状態推定部103を設定するための領域に格納する(ステップS11)。そして処理は終了する。 The design unit 101 stores the gain matrix L calculated in step S9 in the area for setting the state estimation unit 103 in the work data storage unit 115 (step S11). And the process ends.

以上のような処理を実行すれば、等式制約を満たす推定値を算出することが可能なゲイン行列Lを、状態推定部103に使用させることができるようになる。 By executing the above processing, the state estimation unit 103 can use the gain matrix L capable of calculating the estimated value satisfying the equation constraint.

図13を用いて、制御装置1の状態推定部103が実行する処理を説明する。 The process executed by the state estimation unit 103 of the control device 1 will be described with reference to FIG.

状態推定部103は、ステップS1において入力され且つワークデータ格納部115に格納されたパラメータ及びステップS11においてワークデータ格納部115に格納されたゲイン行列Lを読み出す(図13:ステップS13)。 The state estimation unit 103 reads out the parameters input in step S1 and stored in the work data storage unit 115 and the gain matrix L stored in the work data storage unit 115 in step S11 (FIG. 13: step S13).

状態推定部103は、推定値^xを初期化する(ステップS15)。 The state estimation unit 103 initializes the estimated value ^ x (step S15).

状態推定部103は、前時刻の出力yを出力データ格納部111から取得し且つ前時刻の入力uを取得する(ステップS17)。 The state estimation unit 103 acquires the output y of the previous time from the output data storage unit 111 and the input u of the previous time (step S17).

状態推定部103は、ステップS13において読み出されたパラメータ及びゲイン行列と、ステップS17において取得された入力u及び出力yと、算出又は初期化された推定値^xとを用いて、^x←A^x+Bu+L(y−C^x)により推定値^xを更新する(ステップS19)。状態推定部103は、更新後の推定値^xを推定値格納部113に格納する。 The state estimation unit 103 uses the parameter and gain matrix read in step S13, the inputs u and output y acquired in step S17, and the calculated or initialized estimated value ^ x to ^ x ←. The estimated value ^ x is updated by A ^ x + Bu + L (y−C ^ x) (step S19). The state estimation unit 103 stores the updated estimated value ^ x in the estimated value storage unit 113.

状態推定部103は、処理を終了するか判定する(ステップS21)。例えば処理の終了指示を受け付けた場合(ステップS21:Yesルート)には処理は終了し、処理の終了指示を受け付けていない場合(ステップS21:Noルート)には処理はステップS17に戻る。 The state estimation unit 103 determines whether to end the process (step S21). For example, when the process end instruction is received (step S21: Yes route), the process ends, and when the process end instruction is not accepted (step S21: No route), the process returns to step S17.

以上のように、本実施の形態の線形等式保存オブザーバであれば、推定値が等式制約を満たすことが保証されるので、必ずしも真の値に収束していない推定値を用いて制御対象を制御することができるようになる。 As described above, in the case of the linear equation preservation observer of the present embodiment, since the estimated value is guaranteed to satisfy the equation constraint, the controlled object is controlled by using the estimated value that does not necessarily converge to the true value. Will be able to control.

図14及び図15を用いて、制御装置1の通信部105が実行する処理を説明する。 A process executed by the communication unit 105 of the control device 1 will be described with reference to FIGS. 14 and 15.

まず、通信部105が入力uを送信する処理を説明する。通信部105は、入力uを取得する(図14:ステップS31)。入力uは、推定値^xが既に算出された場合には推定値^xに基づき設計部101により算出される。推定値^xが算出されていない場合には、入力uとして初期値が設定される。 First, the process of transmitting the input u by the communication unit 105 will be described. The communication unit 105 acquires the input u (FIG. 14: step S31). The input u is calculated by the design unit 101 based on the estimated value ^ x when the estimated value ^ x has already been calculated. If the estimated value ^ x has not been calculated, the initial value is set as the input u.

通信部105は、ステップS31にて取得された入力uを制御対象装置3に送信する(ステップS33)。そして処理は終了する。 The communication unit 105 transmits the input u acquired in step S31 to the control target device 3 (step S33). And the process ends.

次に、通信部105が出力yを受信する処理を説明する。通信部105は、出力yを制御対象装置3から受信する(図15:ステップS41)。 Next, the process of receiving the output y by the communication unit 105 will be described. The communication unit 105 receives the output y from the control target device 3 (FIG. 15: step S41).

通信部105は、ステップS41にて受信した出力yを出力データ格納部111に格納する(ステップS43)。そして処理は終了する。 The communication unit 105 stores the output y received in step S41 in the output data storage unit 111 (step S43). And the process ends.

以上のような処理を実行すれば、制御対象装置3の制御を適切に行うことができるようになる。 By executing the above processing, the control target device 3 can be appropriately controlled.

[実施の形態2]
第1の実施の形態においては制御装置1がゲイン行列Lを求める処理および推定値を求める処理の両方を実行するが、第2の実施の形態においては両処理が別々の装置において実行される。具体的には、図16に示すように、ネットワーク2に接続された設計装置4がゲイン行列Lを求める処理を実行し、ネットワーク2に接続された状態推定装置5が推定値を算出する処理を実行する。このようなシステム形態であっても、制御対象装置3を適切に制御することができる。
[Embodiment 2]
In the first embodiment, the control device 1 executes both the process of obtaining the gain matrix L and the process of obtaining the estimated value, but in the second embodiment, both processes are executed in separate devices. Specifically, as shown in FIG. 16, the design device 4 connected to the network 2 executes the process of obtaining the gain matrix L, and the state estimation device 5 connected to the network 2 calculates the estimated value. Execute. Even in such a system form, the control target device 3 can be appropriately controlled.

以上本発明の一実施の形態を説明したが、本発明はこれに限定されるものではない。例えば、上で説明した制御装置1及び制御対象装置3の機能ブロック構成は実際のプログラムモジュール構成に一致しない場合もある。 Although one embodiment of the present invention has been described above, the present invention is not limited thereto. For example, the functional block configuration of the control device 1 and the control target device 3 described above may not match the actual program module configuration.

また、処理フローにおいても、処理結果が変わらなければ処理の順番を入れ替えることも可能である。さらに、並列に実行させるようにしても良い。 Also, in the processing flow, it is possible to change the order of processing if the processing result does not change. Further, it may be executed in parallel.

また、制御装置1、設計装置4及び状態推定装置5はネットワーク2に有線で接続されてもよい。 Further, the control device 1, the design device 4, and the state estimation device 5 may be connected to the network 2 by wire.

[付録]
本付録では、条件(a)がWL=0であるための必要十分条件であることの証明を示す。
[appendix]
This appendix shows the proof that the condition (a) is a necessary and sufficient condition for WL = 0.

以下が成立する。 The following holds.

Figure 0006838454
Figure 0006838454

よって、以下が成立する。 Therefore, the following holds.

Figure 0006838454
Figure 0006838454

Uは直交行列であるので、この式はさらに以下の式と等価である。 Since U is an orthogonal matrix, this equation is further equivalent to the following equation.

Figure 0006838454
Figure 0006838454

式(25)の下線部分とVTLとの積が0となるのは、VTLがΛ∈R(n-r)×pを用いて以下のように表されるときに限られる。 The product of the underlined and V T L of formula (25) becomes 0 is limited to when the V T L is expressed as follows using the Λ∈R (nr) × p.

Figure 0006838454
Figure 0006838454

Vは直交行列であるので、これは式(14)と等価である。 Since V is an orthogonal matrix, this is equivalent to equation (14).

以上で付録を終了する。 This is the end of the appendix.

なお、上で述べた制御装置1、設計装置4及び状態推定装置5は、コンピュータ装置であって、図17に示すように、メモリ2501とCPU2503とハードディスク・ドライブ(HDD)2505と表示装置2509に接続される表示制御部2507とリムーバブル・ディスク2511用のドライブ装置2513と入力装置2515とネットワークに接続するための通信制御部2517とがバス2519で接続されている。オペレーティング・システム(OS:Operating System)及び本実施例における処理を実施するためのアプリケーション・プログラムは、HDD2505に格納されており、CPU2503により実行される際にはHDD2505からメモリ2501に読み出される。CPU2503は、アプリケーション・プログラムの処理内容に応じて表示制御部2507、通信制御部2517、ドライブ装置2513を制御して、所定の動作を行わせる。また、処理途中のデータについては、主としてメモリ2501に格納されるが、HDD2505に格納されるようにしてもよい。本発明の実施例では、上で述べた処理を実施するためのアプリケーション・プログラムはコンピュータ読み取り可能なリムーバブル・ディスク2511に格納されて頒布され、ドライブ装置2513からHDD2505にインストールされる。インターネットなどのネットワーク及び通信制御部2517を経由して、HDD2505にインストールされる場合もある。このようなコンピュータ装置は、上で述べたCPU2503、メモリ2501などのハードウエアとOS及びアプリケーション・プログラムなどのプログラムとが有機的に協働することにより、上で述べたような各種機能を実現する。 The control device 1, the design device 4, and the state estimation device 5 described above are computer devices, and as shown in FIG. 17, the memory 2501, the CPU 2503, the hard disk drive (HDD) 2505, and the display device 2509 The display control unit 2507 to be connected, the drive device 2513 for the removable disk 2511, the input device 2515, and the communication control unit 2517 for connecting to the network are connected by a bus 2519. The operating system (OS: Operating System) and the application program for executing the processing in this embodiment are stored in the HDD 2505, and are read from the HDD 2505 to the memory 2501 when executed by the CPU 2503. The CPU 2503 controls the display control unit 2507, the communication control unit 2517, and the drive device 2513 according to the processing contents of the application program to perform a predetermined operation. Further, the data in the process of processing is mainly stored in the memory 2501, but may be stored in the HDD 2505. In the embodiment of the present invention, the application program for performing the above-described processing is stored and distributed on the computer-readable removable disk 2511 and installed from the drive device 2513 to the HDD 2505. It may be installed on the HDD 2505 via a network such as the Internet and a communication control unit 2517. Such a computer device realizes various functions as described above by organically collaborating with the hardware such as the CPU 2503 and the memory 2501 described above and the program such as the OS and the application program. ..

以上述べた本発明の実施の形態をまとめると、以下のようになる。 The embodiments of the present invention described above can be summarized as follows.

本実施の形態の第1の態様に係る制御方法は、(A)制御対象の内部状態を表す変数が満たす条件に対応する方程式に基づき、制御対象の内部状態の推定に用いられるゲイン行列を生成し、(B)生成されたゲイン行列を用いるオブザーバにより生成された、制御対象の内部状態の推定結果に基づき、制御対象を制御する処理を含む。 The control method according to the first aspect of the present embodiment generates a gain matrix used for estimating the internal state of the controlled object based on (A) an equation corresponding to a condition satisfied by a variable representing the internal state of the controlled object. (B) The process of controlling the controlled object based on the estimation result of the internal state of the controlled object generated by the observer using the generated gain matrix is included.

上記方程式に基づくゲイン行列を用いれば推定結果が条件を満たすので、必ずしも真の値に収束しているとは限らない推定結果を用いて、制御対象の制御を適切に行えるようになる。 If the gain matrix based on the above equation is used, the estimation result satisfies the condition, so that the control target can be appropriately controlled by using the estimation result that does not necessarily converge to the true value.

また、上記方程式は、ゲイン行列と、制御対象の内部状態を表す変数が満たす等式制約の係数行列に対する特異値分解の結果に含まれる直交行列とを含んでもよい。このような方程式を利用すれば、LMI可解問題を適用できるようになる。 Further, the above equation may include a gain matrix and an orthogonal matrix included in the result of singular value decomposition for the coefficient matrix of the equation constraint satisfied by the variable representing the internal state of the controlled object. By using such an equation, the LMI solvable problem can be applied.

また、ゲイン行列を生成する処理において、(a1)入力されたデータに基づくLMI(Linear Matrix Inequality)可解問題の解が存在するか判定し、(a2)LMI可解問題の解が存在する場合、LMI可解問題の解と方程式とを用いてゲイン行列を生成してもよい。 Further, in the process of generating the gain matrix, (a1) it is determined whether there is a solution of the LMI (Linear Matrix Inequality) solvable problem based on the input data, and (a2) there is a solution of the LMI solvable problem. , The gain matrix may be generated using the solution of the LMI solvable problem and the equation.

本実施の形態の第2の態様に係る制御装置は、(C)制御対象の内部状態を表す変数が満たす条件に対応する方程式に基づき、制御対象の内部状態の推定に用いられるゲイン行列を生成する生成部(実施の形態における設計部101は生成部の一例である)と、(D)生成部により生成されたゲイン行列を用いるオブザーバ(実施の形態における状態推定部103はオブザーバの一例である)により生成された、制御対象の内部状態の推定結果に基づき、制御対象を制御する制御部(実施の形態における通信部105は制御部の一例である)とを有する。 The control device according to the second aspect of the present embodiment generates a gain matrix used for estimating the internal state of the controlled object based on (C) an equation corresponding to a condition satisfied by a variable representing the internal state of the controlled object. The generation unit (the design unit 101 in the embodiment is an example of the generation unit) and the observer using the gain matrix generated by the (D) generation unit (the state estimation unit 103 in the embodiment is an example of the observer). ), Based on the estimation result of the internal state of the controlled object, has a control unit (the communication unit 105 in the embodiment is an example of the control unit) that controls the controlled object.

なお、上記方法による処理をコンピュータに実行させるためのプログラムを作成することができ、当該プログラムは、例えばフレキシブルディスク、CD−ROM、光磁気ディスク、半導体メモリ、ハードディスク等のコンピュータ読み取り可能な記憶媒体又は記憶装置に格納される。尚、中間的な処理結果はメインメモリ等の記憶装置に一時保管される。 A program for causing a computer to execute the processing by the above method can be created, and the program can be a computer-readable storage medium such as a flexible disk, a CD-ROM, a magneto-optical disk, a semiconductor memory, or a hard disk. Stored in storage. The intermediate processing result is temporarily stored in a storage device such as a main memory.

以上の実施例を含む実施形態に関し、さらに以下の付記を開示する。 The following additional notes will be further disclosed with respect to the embodiments including the above embodiments.

(付記1)
コンピュータに、
制御対象の内部状態を表す変数が満たす条件に対応する方程式に基づき、前記制御対象の内部状態の推定に用いられるゲイン行列を生成し、
生成された前記ゲイン行列を用いるオブザーバにより生成された、前記制御対象の内部状態の推定結果に基づき、前記制御対象を制御する、
処理を実行させる制御プログラム。
(Appendix 1)
On the computer
Based on the equation corresponding to the condition satisfied by the variable representing the internal state of the controlled object, a gain matrix used for estimating the internal state of the controlled object is generated.
The control target is controlled based on the estimation result of the internal state of the control target generated by the observer using the generated gain matrix.
A control program that executes processing.

(付記2)
前記方程式は、前記ゲイン行列と、前記制御対象の内部状態を表す変数が満たす等式制約の係数行列に対する特異値分解の結果に含まれる直交行列とを含む、
付記1記載の制御プログラム。
(Appendix 2)
The equation includes the gain matrix and an orthogonal matrix contained in the result of singular value decomposition on the coefficient matrix of the equation constraint satisfied by the variable representing the internal state of the controlled object.
The control program described in Appendix 1.

(付記3)
前記ゲイン行列を生成する処理において、
入力されたデータに基づくLMI(Linear Matrix Inequality)可解問題の解が存在するか判定し、
前記LMI可解問題の解が存在する場合、前記LMI可解問題の解と前記方程式とを用いて前記ゲイン行列を生成する、
付記1又は2記載の制御プログラム。
(Appendix 3)
In the process of generating the gain matrix,
Determine if there is a solution to the LMI (Linear Matrix Inequality) solvable problem based on the input data.
If a solution to the LMI solvable problem exists, the gain matrix is generated using the solution of the LMI solvable problem and the equation.
The control program according to Appendix 1 or 2.

(付記4)
コンピュータが、
制御対象の内部状態を表す変数が満たす条件に対応する方程式に基づき、前記制御対象の内部状態の推定に用いられるゲイン行列を生成し、
生成された前記ゲイン行列を用いるオブザーバにより生成された、前記制御対象の内部状態の推定結果に基づき、前記制御対象を制御する、
処理を実行する制御方法。
(Appendix 4)
The computer
Based on the equation corresponding to the condition satisfied by the variable representing the internal state of the controlled object, a gain matrix used for estimating the internal state of the controlled object is generated.
The control target is controlled based on the estimation result of the internal state of the control target generated by the observer using the generated gain matrix.
A control method that executes processing.

(付記5)
制御対象の内部状態を表す変数が満たす条件に対応する方程式に基づき、前記制御対象の内部状態の推定に用いられるゲイン行列を生成する生成部と、
前記生成部により生成された前記ゲイン行列を用いるオブザーバにより生成された、前記制御対象の内部状態の推定結果に基づき、前記制御対象を制御する制御部と、
を有する制御装置。
(Appendix 5)
A generator that generates a gain matrix used to estimate the internal state of the controlled object based on the equation corresponding to the condition satisfied by the variable representing the internal state of the controlled object.
A control unit that controls the control target based on an estimation result of an internal state of the control target generated by an observer using the gain matrix generated by the generation unit.
Control device with.

1 制御装置 2 ネットワーク
3 制御対象装置 4 設計装置
5 状態推定装置
101 設計部 103 状態推定部
105 通信部 111 出力データ格納部
113 推定値格納部 115 ワークデータ格納部
301 センサ 302 被制御部
303 通信部
1 Control device 2 Network 3 Control target device 4 Design device 5 State estimation device 101 Design unit 103 State estimation unit 105 Communication unit 111 Output data storage unit 113 Estimated value storage unit 115 Work data storage unit 301 Sensor 302 Controlled unit 303 Communication unit

Claims (3)

コンピュータに、
制御対象の線形状態方程式及び線形出力方程式に含まれる係数行列と前記制御対象の内部状態を表す変数が満たす線形等式制約の係数行列に対する特異値分解の結果に含まれる直交行列とに基づくLMI(Linear Matrix Inequality)可解問題の解が存在するか否かを判定し、
前記LMI可解問題の解が存在する場合、前記LMI可解問題の解から前記制御対象の内部状態の推定に用いられるゲイン行列を生成し、
生成された前記ゲイン行列を用いるオブザーバにより生成された、前記制御対象の内部状態の推定結果に基づき、前記制御対象を制御する、
処理を実行させる制御プログラム。
On the computer
LMI (LMI) based on the coefficient matrix included in the linear state equation and linear output equation of the controlled object and the orthogonal matrix included in the result of singularity decomposition for the coefficient matrix of the linear equation constraint satisfied by the variable representing the internal state of the controlled object. Linear Matrix Inequality) Determines if a solution to the solvable problem exists and
When the solution of the LMI solvable problem exists, a gain matrix used for estimating the internal state of the controlled object is generated from the solution of the LMI solvable problem.
The control target is controlled based on the estimation result of the internal state of the control target generated by the observer using the generated gain matrix.
A control program that executes processing.
コンピュータが、
制御対象の線形状態方程式及び線形出力方程式に含まれる係数行列と前記制御対象の内部状態を表す変数が満たす線形等式制約の係数行列に対する特異値分解の結果に含まれる直交行列とに基づくLMI(Linear Matrix Inequality)可解問題の解が存在するか否かを判定し、
前記LMI可解問題の解が存在する場合、前記LMI可解問題の解から前記制御対象の内部状態の推定に用いられるゲイン行列を生成し、
生成された前記ゲイン行列を用いるオブザーバにより生成された、前記制御対象の内部状態の推定結果に基づき、前記制御対象を制御する、
処理を実行する制御方法。
The computer
LMI (LMI) based on the coefficient matrix included in the linear state equation and linear output equation of the controlled object and the orthogonal matrix included in the result of singularity decomposition for the coefficient matrix of the linear equation constraint satisfied by the variable representing the internal state of the controlled object. Linear Matrix Inequality) Determines if a solution to the solvable problem exists and
When the solution of the LMI solvable problem exists, a gain matrix used for estimating the internal state of the controlled object is generated from the solution of the LMI solvable problem.
The control target is controlled based on the estimation result of the internal state of the control target generated by the observer using the generated gain matrix.
A control method that executes processing.
制御対象の線形状態方程式及び線形出力方程式に含まれる係数行列と前記制御対象の内部状態を表す変数が満たす線形等式制約の係数行列に対する特異値分解の結果に含まれる直交行列とに基づくLMI(Linear Matrix Inequality)可解問題の解が存在するか否かを判定し、前記LMI可解問題の解が存在する場合、前記LMI可解問題の解から前記制御対象の内部状態の推定に用いられるゲイン行列を生成する生成部と、
前記生成部により生成された前記ゲイン行列を用いるオブザーバにより生成された、前記制御対象の内部状態の推定結果に基づき、前記制御対象を制御する制御部と、
を有する制御装置。
LMI (LMI) based on the coefficient matrix included in the linear state equation and linear output equation of the controlled object and the orthogonal matrix included in the result of singularity decomposition for the coefficient matrix of the linear equation constraint satisfied by the variable representing the internal state of the controlled object. Linear Matrix Inequality) Determines if there is a solution to the solvable problem, and if there is a solution to the LMI solvable problem, it is used to estimate the internal state of the controlled object from the solution of the LMI solvable problem. A generator that generates a gain matrix and
A control unit that controls the control target based on an estimation result of an internal state of the control target generated by an observer using the gain matrix generated by the generation unit.
Control device with.
JP2017058723A 2017-03-24 2017-03-24 Control method, control device and control program Active JP6838454B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017058723A JP6838454B2 (en) 2017-03-24 2017-03-24 Control method, control device and control program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017058723A JP6838454B2 (en) 2017-03-24 2017-03-24 Control method, control device and control program

Publications (2)

Publication Number Publication Date
JP2018163395A JP2018163395A (en) 2018-10-18
JP6838454B2 true JP6838454B2 (en) 2021-03-03

Family

ID=63860230

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017058723A Active JP6838454B2 (en) 2017-03-24 2017-03-24 Control method, control device and control program

Country Status (1)

Country Link
JP (1) JP6838454B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116224099B (en) * 2023-05-06 2023-07-21 力高(山东)新能源技术股份有限公司 Method for dynamically and adaptively estimating battery SOC

Also Published As

Publication number Publication date
JP2018163395A (en) 2018-10-18

Similar Documents

Publication Publication Date Title
Stojanovic et al. Joint state and parameter robust estimation of stochastic nonlinear systems
WO2021092980A1 (en) Longitudinal federated learning optimization method, apparatus and device, and storage medium
JP2015181041A (en) Method and system for empirical modeling of time-varying, parameter-varying, and nonlinear systems via iterative linear subspace computation
Van et al. A robust fault diagnosis and accommodation scheme for robot manipulators
EP2818947A1 (en) Method for automatically setting controller bandwidth
JP2014215822A (en) State estimating apparatus, method, and program
WO2016129078A1 (en) Route selection device and route selection program
JP6838454B2 (en) Control method, control device and control program
KR20230008005A (en) Method and apparatus for determining information
US11088960B2 (en) Information processing apparatus and verification system
JP6617771B2 (en) Linear parameter variation model estimation system, method and program
US20200326952A1 (en) Modification procedure generation device, modification procedure generation method and storage medium for storing modification procedure generation program
US11526690B2 (en) Learning device, learning method, and computer program product
JP6640615B2 (en) Orbit calculation device and orbit calculation program
JP5634355B2 (en) Target tracking system and program and method thereof, angle tracking device and program and method thereof, target tracking device and program and method thereof
WO2018008391A1 (en) Linear parameter varying model estimation system, method, and program
US20220172091A1 (en) Learning parameters of bayesian network using uncertain evidence
JP2008165412A (en) Performance-computing device
JP2019101712A (en) Abnormality estimation device, abnormality estimation method and program
JP7111159B2 (en) Information sharing device, information sharing method and information sharing program
JP2022167093A (en) Information processing device, information processing method, and program
JP2018084493A (en) Sensor network analyzer, sensor network analysis method and sensor network analysis program
JP2018040601A (en) Bias error estimation device, bias error estimation system, bias error estimation method and bias error estimation program
Amato et al. On the region of asymptotic stability of nonlinear quadratic systems
JP2015075830A (en) Parallel processing management program, parallel processing management method, and parallel processing management device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191212

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20201014

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20201020

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201218

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210125

R150 Certificate of patent or registration of utility model

Ref document number: 6838454

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150