JP2001290514A - Simulation method of control simulation device - Google Patents
Simulation method of control simulation deviceInfo
- Publication number
- JP2001290514A JP2001290514A JP2000105094A JP2000105094A JP2001290514A JP 2001290514 A JP2001290514 A JP 2001290514A JP 2000105094 A JP2000105094 A JP 2000105094A JP 2000105094 A JP2000105094 A JP 2000105094A JP 2001290514 A JP2001290514 A JP 2001290514A
- Authority
- JP
- Japan
- Prior art keywords
- module
- simulation
- equation
- state variable
- calculating
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000004088 simulation Methods 0.000 title claims abstract description 38
- 239000011159 matrix material Substances 0.000 claims abstract description 21
- 230000010354 integration Effects 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims description 22
- 238000005094 computer simulation Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
Landscapes
- Complex Calculations (AREA)
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は、制御系のシミュレ
ーションを行う装置、特に高速なシミュレーションが要
求される制御系のシミュレーション装置におけるシミュ
レーション方法に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an apparatus for simulating a control system, and more particularly to a simulation method for a simulation apparatus for a control system which requires a high-speed simulation.
【0002】[0002]
【従来の技術】自然界の現象の中には、数学的に常微分
方程式で記述できるものが多く、制御の分野において
も、制御対象となるプラントおよびそれを制御するコン
トローラなどを常微分方程式でモデル化(記述)し、そ
の常微分方程式を解析することによりプラントおよびコ
ントローラを含めた系全体の挙動を把握する。そのため
にはモデル化された常微分方程式を解く必要があるが、
しかし解析的に解ける常微分方程式は非常に限られ、現
実的には計算機を用いて数値的に解くという手法(数値
積分法)が採られるのが一般的である。2. Description of the Related Art Many phenomena in the natural world can be mathematically described by ordinary differential equations. In the field of control, a plant to be controlled and a controller for controlling the same are modeled by ordinary differential equations. Then, the behavior of the whole system including the plant and the controller is grasped by analyzing (or describing) the ordinary differential equation. To do this, we need to solve the modeled ODE,
However, the ordinary differential equations that can be solved analytically are extremely limited, and in practice, a method of numerically solving using a computer (numerical integration method) is generally adopted.
【0003】そのような数値積分法としてこれまで様々
な手法が提案されており、シミュレーションの対象や目
的などによってそれら手法の中から適切な手法を選択し
用いることとなる。主な手法としてはルンゲクッタ法に
代表される陽的積分法、台形法に代表される陰的積分法
などがある。Various methods have been proposed as such a numerical integration method, and an appropriate method is selected from these methods depending on the object or purpose of the simulation. The main methods include an explicit integration method represented by the Runge-Kutta method and an implicit integration method represented by the trapezoidal method.
【0004】[0004]
【発明が解決しようとする課題】シミュレーションの対
象期間を長く設定する場合、特にシミュレーション対象
が大規模複雑となると、シミュレーションの各時間ステ
ップに要する計算量は膨大なものとなり、シミュレーシ
ョンに要する時間が長くなる。この計算量を減らすため
には時間ステップを大きく設定することが考えられる
が、陽的積分法は、時間ステップを大きく設定すると原
理上計算誤差が大きくなる傾向があり、シミュレーショ
ン結果の信頼性が損なわれる。一方、陰的積分法は、一
般的に時間ステップを大きく設定しても計算誤差が大き
くなりにくいという性質を有しているが、各時間ステッ
プに要する計算量が陽的積分法よりも大きくなり、結果
的にシミュレーションの速度が向上しない場合があると
いう問題を有している。When the simulation target period is set to be long, especially when the simulation target is large-scale and complicated, the amount of calculation required for each time step of the simulation is enormous, and the time required for the simulation is long. Become. To reduce the amount of calculation, it is conceivable to set a large time step.However, in the explicit integration method, if the time step is set large, the calculation error tends to increase in principle, and the reliability of the simulation result is impaired. It is. On the other hand, the implicit integration method generally has the property that the calculation error is unlikely to increase even if a large time step is set, but the amount of calculation required for each time step is larger than the explicit integration method. As a result, there is a problem that the simulation speed may not be improved.
【0005】[0005]
【課題を解決するための手段】上記課題を解決するため
に成された本発明では、制御系を構成する複数のモジュ
ールをつなぎ合わせることによって1つの制御系を構成
し、各モジュールの動特性を支配する常微分方程式で表
された状態方程式を、非線型連立方程式の求解にニュー
トン法を用いた陰的数値積分法により時間ステップ毎に
解くことによって、毎時間ステップにおける各モジュー
ルの出力値を計算する制御シミュレーション装置のシミ
ュレーション方法において、陰的数値積分法における非
線型連立方程式をニュートン法により求解するに当た
り、各状態変数を微小に変動させ、その変動に対する各
状態変数の時間微係数の変動量を計算し、それらの値か
らニュートン法のヤコビアン行列を算出する際に、各モ
ジュールを、その微小に変動させる状態を有するモジュ
ールの出力を直接入力とするモジュールと、入力を時間
遅れなく出力とするモジュールを1個以上介して入力と
するモジュールと、入力としないモジュールとに、系を
構成するモジュールの接続関係および各モジュールの出
力方程式から選別することで、状態変数の時間微係数の
変動量が0とならない状態変数を特定し、その変動量の
みを計算によって求めることによって、ヤコビアン行列
の各要素の値を計算する際の計算量を減らし、よって数
値積分に要する時間を短縮し、シミュレーションを高速
化する。According to the present invention which has been made to solve the above-mentioned problems, one control system is formed by connecting a plurality of modules constituting a control system, and the dynamic characteristics of each module are determined. Calculate the output value of each module at each time step by solving the equation of state expressed by the governing ordinary differential equation at each time step by solving the nonlinear simultaneous equations using the implicit numerical integration method using Newton's method In solving the nonlinear simultaneous equations in the implicit numerical integration method by the Newton method in the simulation method of the control simulation device, each state variable is minutely changed, and the amount of change of the time derivative of each state variable with respect to the change is calculated. When calculating and calculating the Newton method Jacobian matrix from those values, each module is The modules that constitute a system include a module that directly receives an output of a module that has a state that fluctuates, a module that receives an input via at least one module that outputs an input without time delay, and a module that does not receive an input. By determining from the connection relation of the module and the output equation of each module, the state variable in which the variation of the time derivative of the state variable does not become 0 is specified, and only the variation is obtained by calculation, thereby calculating each element of the Jacobian matrix. Reduces the amount of calculation when calculating the value of, thereby shortening the time required for numerical integration and speeding up the simulation.
【0006】[0006]
【発明の実施の形態】本実施例の制御シミュレーション
装置を実現する構成を図1に示す。図1で、101は中
央演算装置(CPU)、102はプログラムの格納およ
び演算結果の格納に用いられる常時書き込み読み出し可
能記憶装置(RAM)、103は計算結果を表示するた
めのディスプレイ装置、104はプログラムおよびデー
タを格納するためのハードディスク装置、105は10
1から104でを結合するバスである。FIG. 1 shows a configuration for realizing a control simulation apparatus according to this embodiment. In FIG. 1, 101 is a central processing unit (CPU), 102 is a constantly writable and readable storage device (RAM) used for storing programs and calculation results, 103 is a display device for displaying calculation results, and 104 is a display device. Hard disk drive for storing programs and data, 105 is 10
A bus connecting 1 to 104.
【0007】この構成により、図2のフローチャートに
従って数値積分が行われる。以下図3のようなシミュレ
ーション対象を例にし、陰的積分法として台形法を用い
た場合について説明する。なお、図3の301は出力を
y1 、状態変数をx1 ,x2、状態方程式d(x1 ,x
2 )/dt=f1 (x1 ,x2 )、出力方程式をy1=
h1 (x1 ,x2 )を動特性として持つモジュール、3
02は入力をu2 、状態変数をx3 、出力をy2 、状態
方程式をdx3 /dt=f2 (x3 )+g2 (u2 )、
出力方程式をy2 =h2 (x3 )として持つモジュー
ル、303は入力をu3 、状態変数をx4 ,x5 、状態
方程式をd(x4 ,x5 )/dt=f3 (x4 ,x5 )
+g3 (u3 )として持つモジュール、304は入力を
u4 、出力をy4 、状態変数をx6 、状態方程式をdx
6 /dt=f4 (x6 )+g4 (u 4 )、出力方程式を
y4 =h4 (x6 )+k4 (u4 )として持つモジュー
ル、305は入力をu5 、状態変数をx7 、状態方程式
をdx7 /dt=f5 (x7)+g5 (u5 )として持
つモジュールである。ここで、304は、出力方程式が
その入力u4 の関数となっているため、入力u4 が時間
遅れなく出力y4 に現れるという性質(直達成分)を持
っているモジュールであり、その他のモジュールはその
ような性質は持っていない。このような各モジュールの
直達成分の有無は、直達成分フラグによって各モジュー
ル毎に記憶される。また、モジュール間の接続関係は、
102に示すように、各モジュールの次モジュールへの
ポインタによってRAM上に記憶される。With this configuration, the flowchart of FIG.
Therefore, numerical integration is performed. A simulation as shown in Figure 3 below
Using the trapezoidal method as an implicit integration method
The following describes the case where The output 301 in FIG.
y1, And state variables x1, XTwo, The state equation d (x1, X
Two) / Dt = f1(X1, XTwo), The output equation is y1=
h1(X1, XTwo) As a dynamic characteristic, 3
02 is input uTwo, And state variables xThree, Output yTwo,Status
Dx equationThree/ Dt = fTwo(XThree) + GTwo(UTwo),
Output equation is yTwo= HTwo(XThree) Module
303, input uThree, And state variables xFour, XFive,Status
The equation is d (xFour, XFive) / Dt = fThree(XFour, XFive)
+ GThree(UThree)), 304 is an input
uFour, Output yFour, And state variables x6And the state equation is dx
6/ Dt = fFour(X6) + GFour(U Four), The output equation
yFour= HFour(X6) + KFour(UFour) Module
305 input uFive, And state variables x7, Equation of state
To dx7/ Dt = fFive(X7) + GFive(UFive)
One module. Where 304 is the output equation
That input uFourInput uFourIs time
Output y without delayFourHas the property of appearing in
The other modules are
It does not have such properties. Of each such module
The direct achievement flag is set for each module by the direct achievement flag.
Is stored for each file. The connection relationship between modules is
As shown in 102, each module is assigned to the next module.
It is stored on the RAM by the pointer.
【0008】またシミュレーションの現在時刻をt、時
間ステップ幅をh、状態変数をまとめたベクトルをX
(=(x1 ,x2 ,x3 ,x4 ,x5 ,x6 ,
x7 ))、関数F(X)をF(X)=(f1 (x1 ,x
2 )f2 (x3 )f3 (x4 ,x5 )f 4 (x6 )f5
(x7 ))、時刻tにおけるXをXt 、時刻t+hにお
けるXをXt+h と表す。[0008] Also, the current time of the simulation is t, hour
H is the interval step width, and X is a vector summarizing the state variables.
(= (X1, XTwo, XThree, XFour, XFive, X6,
x7)) And the function F (X) is expressed as F (X) = (f1(X1, X
Two) FTwo(XThree) FThree(XFour, XFive) F Four(X6) FFive
(X7)), X at time t is XtAt time t + h
X to Xt + hIt expresses.
【0009】まず、S201において初期状態X0 を0
に設定し、S202において現在時刻tを0に設定す
る。また時間ステップhは操作者が任意に与えるものと
する。これらは102に示すようにRAMに記憶され
る。S203において時刻tから時間ステップh後の状
態Xt+h を非線型連立方程式 (Xt +Xt+h )/h=(F(Xt )+F(Xt+h ))/2 (1) を満たすようなXt+h を図4に示すステップに従ってニ
ュートン法を用いて求める。そしてS204において現
在時刻をhだけ進め、S205において終了判断が行わ
れ、まだ終了していない場合はS203に戻る。First, in step S201, the initial state X 0 is set to 0.
, And the current time t is set to 0 in S202. The time step h is arbitrarily given by the operator. These are stored in RAM as shown at 102. In S203, the state Xt + h after the time step h from the time t is converted into the nonlinear simultaneous equation ( Xt + Xt + h ) / h = (F ( Xt ) + F ( Xt + h )) / 2 (1) Xt + h that satisfies is satisfied by using the Newton method according to the steps shown in FIG. Then, in S204, the current time is advanced by h. In S205, an end determination is made, and if it has not ended, the process returns to S203.
【0010】次に図4のニュートン法について述べる。
式(1)を満たすXt+h は図4に示すステップに従って
ニュートン法を用いて求める。S401では、式(1)
の解Xt+h の暫定解Xt+h (0) を適当に定め、S402
では、そのi行j列(i,j=17)の要素が J(i,j)=δij+0.5×h×(Δ(dxi /dt)/Δxj ) (2) で定義されるヤコビアン行列(以下J行列と呼ぶ)を計
算する。ここで、式(2)におけるδijはクロネッカー
のデルタ δij=1(i=jの場合) 0(i≠jの場合) (3) を表し、Δ(dxi /dt)/Δxj は、状態変数xj
がΔxj だけ微小変動したときの状態変数xi の時間微
係数の微小変動量を表しており、この値は各モジュール
の状態方程式および出力方程式から計算される。例え
ば、j=1,i=3の場合は、x1 がΔx1 だけ微小変
動したときのモジュール301の出力y1(=u2 )を
計算し、u2 からモジュール302の状態変数x3 の時
間微係数の値を求め、予め計算しておいたx1 が微小変
動する前のx3 の時間微係数の値からの差分を計算する
ことで求められる。次に、S403でJ行列を用いて解
の更新方向dを線形連立方程式 Jd=−F(Xt+h ) (4) を解き求める。次に、S404において Xt+h (n+1) ←Xt+h (n) +μd (5) と更新する。そしてS405において解の収束を判定
し、収束していない場合S402に戻り、収束している
場合は終了する。Next, the Newton method of FIG. 4 will be described.
X t + h that satisfies the equation (1) is obtained by using the Newton method according to the steps shown in FIG. In S401, equation (1)
Define the solution X t + h of interim solution X t + h (0) appropriately, S402
So the elements of the row i and column j (i, j = 17) is defined by the J (i, j) = δ ij + 0.5 × h × (Δ (dx i / dt) / Δx j) (2) A Jacobian matrix (hereinafter referred to as a J matrix) is calculated. Here, [delta] ij in equation (2) (in the case of i = j) delta [delta] ij = 1 Kronecker (for i ≠ j) 0 represents a (3), Δ (dx i / dt) / Δx j is , State variable x j
Represents the minute variation of the time derivative of the state variable x i when the value fluctuates by Δx j , and this value is calculated from the state equation and output equation of each module. For example, when j = 1 and i = 3, the output y 1 (= u 2 ) of the module 301 when x 1 slightly fluctuates by Δx 1 is calculated, and the state variable x 3 of the module 302 is calculated from u 2 . obtains the value of the time derivative, obtained by x 1 which had been previously calculated to calculate the difference from the value of the time derivative of the previous x 3 to slight change. Next, in S403, the solution update direction d is solved by using the J matrix to solve the linear simultaneous equation Jd = −F (X t + h ) (4). Next, in S404, Xt + h (n + 1) ← Xt + h (n) + μd (5) is updated. Then, in step S405, the convergence of the solution is determined. If the convergence has not occurred, the process returns to step S402, and if the convergence has occurred, the process ends.
【0011】S402では、J行列を計算する必要があ
るが、従来技術では全てのi,jについて式(2)を計
算しJ行列を計算していたが、シミュレーション対象が
大規模複雑になるに従い状態数が膨大になり、式(2)
の計算量が膨大になり、結果シミュレーションの速度が
著しく低くなる。例えば図3では、状態変数が7個ある
ので式(2)の計算を72 =49回行う必要がある。In S402, it is necessary to calculate the J matrix. In the prior art, the formula (2) is calculated for all i and j to calculate the J matrix. The number of states becomes enormous, and equation (2)
The amount of calculation becomes huge, and as a result, the speed of the simulation becomes extremely low. For example, in FIG. 3, since there are seven state variables, the calculation of equation (2) needs to be performed 7 2 = 49 times.
【0012】そこで本発明では、J行列を計算する際に
全てのi,jについて式(2)を計算するのではなく、
J行列のうち0となる要素を前もって特定しておき、そ
の要素に関しては式(2)の計算を省く事でJ行列を計
算する際の全体の計算量を低減させる。一般的に、シミ
ュレーション装置におけるJ行列は0要素が多いため、
大幅な計算量の低減が期待される。Therefore, in the present invention, when calculating the J matrix, instead of calculating equation (2) for all i and j,
An element that becomes 0 in the J matrix is specified in advance, and the calculation amount of the J matrix is reduced by omitting the calculation of Expression (2) for the element, thereby reducing the total amount of calculation when calculating the J matrix. In general, the J matrix in a simulation device has many 0 elements,
A large reduction in the amount of calculation is expected.
【0013】本発明でのJ行列の計算は図5に示すステ
ップに従って計算する。S501において、変動させる
状態変数xj (j=1〜7)を任意に選択する。そして
S502においてそのxj の変動幅Δxj をパラメータ
により適当に設定する。S503においてxj の変動の
影響を求める状態変数xi を任意に選択する。The calculation of the J matrix in the present invention is performed according to the steps shown in FIG. In S501, a state variable x j (j = 1 to 7) to be changed is arbitrarily selected. And appropriately setting the variation range [Delta] x j of the x j by the parameter in S502. Arbitrarily select a state variable x i to determine the effects of changes in x j in S503.
【0014】S504においてxj を有するモジュール
の出力がxi を有するモジュールの入力に直接接続され
ているかを、102のRAM上に記憶した各モジュール
の次モジュールポインタを順次辿って判断し、直接接続
されていない場合S505で途中1つ以上のモジュール
を経由して間接的に接続しているかを、同じく102の
RAM上に記憶された各モジュールの次モジュールへの
ポインタを辿ることで判断する。直接接続されている場
合、および間接的に接続されかつ途中経由しているモジ
ュールがすべて直達成分(入力信号が微分要素もしくは
積分要素を経由せず出力信号に現れる成分)を有してい
る場合、S507において式(2)に従ってJ行列のi
行j列の成分を計算する。In step S504, it is determined whether the output of the module having x j is directly connected to the input of the module having x i by sequentially following the next module pointer of each module stored in the RAM of 102, and is directly connected. If not, it is determined in step S505 whether the connection is indirectly via one or more modules in the middle by tracing the pointer to the next module of each module stored in the RAM 102. When directly connected, and when all the modules connected indirectly and passing therethrough have direct achievements (components in which the input signal does not pass through the differential element or the integral element and appears in the output signal), In S507, i of the J matrix according to the equation (2)
Calculate the components in row j column.
【0015】ここで、直達成分の有無は102のRAM
上に記憶された各モジュールの直達成分存在フラグによ
って判断する。それ以外の場合は式(2)の計算を行う
までもなくJ行列のi行j列の成分は0となるため、式
(2)の計算を省略する事ができ、計算量を削減する事
ができる。また、xi とxj が同じモジュールに含まれ
ている場合もJ行列のi行j列の成分は0となり、i=
jの場合は式(2)より1となる。例えば、j=1の場
合、状態変数xj を有しているモジュールは301であ
り、301の出力を直接入力としているのは302と3
04であり、また305は直達項をもつモジュール30
4を経由して301に接続しているので、状態変数xj
が微少変動する際に計算すべき時間微係数は、Δ(dx
3 /dt)/Δx1 、Δ(dx6 /dt)/Δx1 、Δ
(dx7 /dt)/Δx1 であり、これらの計算結果と
式(2)からJ(3,1),J(6,1),J(7,
1)を計算する。また、状態変数x2 はx1 と同じモジ
ュールに含まれているのでJ(2,1)=0となり、ま
たJ(1,1)は式(2)から1となる。状態変数
x 4 ,x5 を有するモジュールは、301の出力を直接
入力としておらず、また直達成分をもつモジュールを経
由しても接続されていないためJ(4,1),J(5,
1)は0となる。次に、S509において全てのxj お
よび全てのxj について行ったか判断を行い、行った場
合終了し、それ以外の場合次のxi もしくはxj を選択
し、同様の処理を全てのxi およびxj について繰り返
す。Here, the presence or absence of the direct achievement is determined by 102 RAMs.
According to the direct achievement presence flag of each module stored above
I judge. Otherwise, calculate equation (2)
Since the component at the i-th row and the j-th column of the J matrix becomes 0 soon, the expression
The calculation of (2) can be omitted, and the calculation amount can be reduced.
Can be. Also, xiAnd xjAre included in the same module
Also, the component of the i-th row and the j-th column of the J matrix becomes 0, and i =
In the case of j, it becomes 1 according to equation (2). For example, if j = 1
The state variable xjThe module having
The reason why the output of 301 is directly input is 302 and 3
04 and a module 305 having a direct term
4, the state variable xj
Is small when Δ fluctuates slightly, Δ (dx
Three/ Dt) / Δx1, Δ (dx6/ Dt) / Δx1, Δ
(Dx7/ Dt) / Δx1And these calculation results and
From equation (2), J (3,1), J (6,1), J (7,
Calculate 1). Also, the state variable xTwoIs x1Same moji as
J (2,1) = 0 because it is included in the
J (1,1) becomes 1 from Expression (2). State variables
x Four, XFiveModule directly outputs the output of 301
The module that is not input and has direct achievement
J (4,1), J (5,
1) becomes 0. Next, in step S509, all xjYou
And all xjJudge whether you went about
Ends, otherwise the next xiOr xjchoose
And the same processing is performed for all xiAnd xjRepeat about
You.
【0016】以上の方法により、図3で示されるシミュ
レーション対象では、J(1,1),J(2,2),J
(3,3),J(4,4),J(5,5),J(6,
6),J(7,7)が1、J(3,1),J(6,
1),J(7,1),J(3,2),J(6,2),J
(7,2),J(4,3),J(5,3),J(7,
6)について式(2)を用いて値を計算、残りの残り3
2要素は式(2)を計算するまでもなく0と特定でき、
式(2)の計算回数を16回に低減する事ができる。According to the above method, in the simulation object shown in FIG. 3, J (1,1), J (2,2), J (1,2)
(3,3), J (4,4), J (5,5), J (6,
6), J (7,7) is 1, J (3,1), J (6,
1), J (7, 1), J (3, 2), J (6, 2), J
(7,2), J (4,3), J (5,3), J (7,
For 6), calculate the value using equation (2), and calculate the remaining 3
Two elements can be specified as 0 without calculating equation (2),
The number of calculations of the equation (2) can be reduced to 16 times.
【0017】[0017]
【発明の効果】本発明により、シミュレーション対象で
あるモデルの特性を用いて、陰的数値積分法を用いた制
御シミュレーションの実行に要する計算量を低減するこ
とができ、これまで知られている制御シミュレーション
装置よりも高速な制御シミュレーションが可能となるこ
とから、特に大規模複雑な制御シミュレーションを高速
に実行する用途に有用される。According to the present invention, the amount of calculation required to execute a control simulation using an implicit numerical integration method can be reduced using characteristics of a model to be simulated. Since the control simulation can be performed at a higher speed than that of the simulation device, it is particularly useful for the purpose of executing a large-scale and complicated control simulation at a high speed.
【図1】本実施例の装置構成を示す図である。FIG. 1 is a diagram showing a device configuration of the present embodiment.
【図2】本実施例の制御シミュレーションの手順を示す
図である。FIG. 2 is a diagram illustrating a procedure of a control simulation according to the present embodiment.
【図3】本実施例の制御シミュレーションのモデル例を
表す図である。FIG. 3 is a diagram illustrating a model example of a control simulation according to the present embodiment.
【図4】本実施例の制御シミュレーションにおけるニュ
ートン法を用いた非線型連立方程式の求解の手順を示す
図である。FIG. 4 is a diagram showing a procedure for solving a non-linear system of equations using the Newton method in the control simulation of the present embodiment.
【図5】本実施例におけるニュートン法のヤコビアン行
列を計算する手順を示す図である。FIG. 5 is a diagram showing a procedure for calculating a Jacobian matrix of the Newton method in the embodiment.
101…CPU 102…RAM 103…CRT 104…HDD 105…バス S201〜S205…本実施例の制御シミュレーション
の手順である。 S301〜S305…本実施例の制御シミュレーション
のモデル例における各モジュールである。 S401〜S405…本実施例の制御シミュレーション
におけるニュートン法の手順である。 S501〜S509…本実施例におけるニュートン法の
ヤコビアン行列を計算する手順である。101: CPU 102: RAM 103: CRT 104: HDD 105: Buses S201 to S205: Control simulation procedures of the present embodiment. S301 to S305 are modules in the control simulation model example of the present embodiment. S401 to S405 are the procedures of the Newton method in the control simulation of the present embodiment. S501 to S509 are the procedures for calculating the Jacobian matrix of the Newton method in this embodiment.
Claims (1)
なぎ合わせることによって1つの制御系を構成し、各モ
ジュールの動特性を支配する常微分方程式で表された状
態方程式を、非線型連立方程式の求解にニュートン法を
用いた陰的数値積分法により時間ステップ毎に解くこと
によって、毎時間ステップにおける各モジュールの出力
値を計算する制御シミュレーション装置におけるシミュ
レーション方法において、前記ニュートン法により求解
するに当たり、前記状態変数を微小に変動させ、その変
動に対する各状態変数の時間微係数の変動量を計算し、
それらの値からニュートン法のヤコビアン行列を算出す
る際に、各モジュールから、微小に変動させる状態変数
を有するモジュールの出力を直接入力とするモジュール
と、入力を時間遅れなく出力とするモジュールを1つ以
上介して入力とするモジュールとを微少に変動させる状
態変数毎に選別して、前記選別したモジュールに含まれ
る状態変数だけの時間微係数の変動量を計算することに
より、ヤコビアン行列の各要素の値を計算する際の計算
量を減らし、シミュレーションを高速化することを特徴
とする制御シミュレーション装置におけるシミュレーシ
ョン方法。1. A control system is constructed by connecting a plurality of modules constituting a control system, and a state equation expressed by an ordinary differential equation governing dynamic characteristics of each module is converted into a nonlinear simultaneous equation. By solving for each time step by an implicit numerical integration method using Newton's method for solving, in a simulation method in a control simulation device for calculating an output value of each module in each time step, when solving by the Newton's method, The state variable is fluctuated minutely, and the fluctuation amount of the time derivative of each state variable with respect to the fluctuation is calculated,
When calculating the Jacobian matrix of the Newton's method from these values, one module directly receives the output of a module having a state variable that fluctuates minutely, and one module outputs the input without time delay from each module. The module to be input through the above is selected for each state variable that fluctuates minutely, and the amount of change of the time derivative of only the state variable included in the selected module is calculated, whereby each element of the Jacobian matrix is calculated. A simulation method in a control simulation device, characterized in that a calculation amount for calculating a value is reduced and a simulation is speeded up.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000105094A JP2001290514A (en) | 2000-04-06 | 2000-04-06 | Simulation method of control simulation device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000105094A JP2001290514A (en) | 2000-04-06 | 2000-04-06 | Simulation method of control simulation device |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2001290514A true JP2001290514A (en) | 2001-10-19 |
Family
ID=18618520
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000105094A Pending JP2001290514A (en) | 2000-04-06 | 2000-04-06 | Simulation method of control simulation device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2001290514A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101071070B1 (en) | 2009-03-17 | 2011-10-10 | 군산대학교산학협력단 | Method for solving equations and recording medium storing the method |
CN103793595A (en) * | 2013-12-30 | 2014-05-14 | 中山职业技术学院 | Numerical simulation method of finned tube heat exchanger for variable frequency CO2 heat-pump water heaters |
-
2000
- 2000-04-06 JP JP2000105094A patent/JP2001290514A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101071070B1 (en) | 2009-03-17 | 2011-10-10 | 군산대학교산학협력단 | Method for solving equations and recording medium storing the method |
CN103793595A (en) * | 2013-12-30 | 2014-05-14 | 中山职业技术学院 | Numerical simulation method of finned tube heat exchanger for variable frequency CO2 heat-pump water heaters |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2008521138A (en) | Method for analysis of control systems | |
WO2018152319A1 (en) | Spreadsheet recalculation algorithm for directed acyclic graph processing | |
CN107133190A (en) | The training method and training system of a kind of machine learning system | |
US20130226530A1 (en) | Mesh generation system | |
JP6290029B2 (en) | Production control support device, production control support method and program | |
EP2128778A1 (en) | Simulation control program, recording medium, simulator and simulation control method | |
CN103942380B (en) | Graphical control system design and simulation tool | |
JP2001290514A (en) | Simulation method of control simulation device | |
CN111651922A (en) | Method, system, device and storage medium for analyzing simulated manual calculation of force method calculation | |
Miskov-Zivanov et al. | THiMED: Time in hierarchical model extraction and design | |
Feng et al. | Unifying relationships between complexity and stability in mutualistic ecological communities | |
JP2807383B2 (en) | Operation limit value monitoring device | |
US20230206094A1 (en) | Inferring device, inferring method and non-transitory computer readable medium | |
CN105354385A (en) | Implementing method of grid sequencing in semiconductor process simulation | |
EP3663992A1 (en) | Learning-finished model integration method, device, program, ic chip, and system | |
Jafari et al. | Decoupled Interpolation-based Method for Numerical Simulation of Digital Controllers | |
JPH02226186A (en) | Simulation method | |
US20240119204A1 (en) | Modeling simulation system for plasma control | |
JPH01217562A (en) | System for generating computing grid | |
CN107194083B (en) | Lyapunov time lag stability criterion optimization method based on redundancy item discrimination | |
CN104899183A (en) | Rapid Gauss-Jordan elimination method for symbolic linear system | |
JPS59117660A (en) | Processing system for restoration of simulation state | |
US20120232866A1 (en) | Method for real time computation of the state variables of a hybrid differential-algebraic process model | |
JPH11154170A (en) | Circuit simulator | |
Busch | Continuous Simulation with Ordinary Differential Equations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060906 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090126 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090203 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090330 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20100316 |