JPH1115362A - Simulation method, and simulator - Google Patents

Simulation method, and simulator

Info

Publication number
JPH1115362A
JPH1115362A JP16260697A JP16260697A JPH1115362A JP H1115362 A JPH1115362 A JP H1115362A JP 16260697 A JP16260697 A JP 16260697A JP 16260697 A JP16260697 A JP 16260697A JP H1115362 A JPH1115362 A JP H1115362A
Authority
JP
Japan
Prior art keywords
simulation
time step
value
simulation time
input value
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
Application number
JP16260697A
Other languages
Japanese (ja)
Inventor
Taiji Tsutsui
泰治 筒井
Yasunori Sueyoshi
康則 末吉
Takeshi Toshi
武史 都志
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.)
Kubota Corp
Original Assignee
Kubota 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 Kubota Corp filed Critical Kubota Corp
Priority to JP16260697A priority Critical patent/JPH1115362A/en
Publication of JPH1115362A publication Critical patent/JPH1115362A/en
Pending legal-status Critical Current

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

PROBLEM TO BE SOLVED: To enable the follow-up to the rapid input change with excellent accuracy by performing the time increment determination procedure to determine the simulation time increment immediately after an event related to the prescribed input value to a simulator is generated before the simulation is performed. SOLUTION: The time increment determination procedure to determine the simulation time increment for the simulation to be newly performed immediately after generation of an event related to one or two or more input values to a process simulator 4, is first performed according to the change in the input value from a controller 2 immediately after a time increment determining means 5 completes the communication. Then, the process simulator 4 performs the simulation with the simulation time increment obtained by performing the time increment determination procedure. Each simulation time increment of the subsequently continuous simulations is performed by the process simulator 4 with the simulation time increment determined in the immediately previous simulation.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、シミュレータの出
力値の変化量に応じてシミュレーション時間刻みを変化
させるシミュレーション方法及びシミュレータに関し、
特に、実時間でシミュレーションを実行する必要のある
プラント施設の運転者訓練用の運転訓練シミュレータ及
びそのシミュレーション方法に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a simulation method and a simulator for changing a simulation time interval according to a change amount of an output value of a simulator.
In particular, the present invention relates to an operation training simulator for operator training of a plant facility that needs to execute a simulation in real time and a simulation method thereof.

【0002】[0002]

【従来の技術】従来、この種のシミュレーション方法及
びシミュレータでは、シミュレーションの対象が、解析
的に厳密に解法できない数学モデルによって表現される
場合は、一般に多変数の入力と多変数の出力と多変数の
内部状態を与えて、その各変数間の関係を状態方程式と
出力方程式で表した数式モデルを計算機で数値的に解法
してシミュレーションを実行する。この状態方程式をオ
イラー法やルンゲクッタ法等の解法手法によって解法す
る際のシミュレーション時間刻みとして、例えば、0.
1μ秒〜10秒といった範囲を設定し、これらの解法手
法によって算出された出力値とシミュレーション前の旧
出力値との変化量を所定の評価基準に基づいて評価し、
評価結果としての変化量が小さければ、シミュレーショ
ン時間刻みを長くし、その変化量が大きければ、シミュ
レーション時間刻みを短くして次ステップのシミュレー
ションを逐次実行する手法が採られ、シミュレーション
の精度と実行時間の短縮が図られている。尚、シミュレ
ータの出力値は、シミュレーションを実行して得られた
結果を意味し、具体的には、前記状態方程式を解法した
場合は、内部状態変数を意味する。
2. Description of the Related Art Conventionally, in a simulation method and a simulator of this kind, when a simulation target is expressed by a mathematical model that cannot be analytically rigorously solved, generally a multivariable input, a multivariable output and a multivariable Is given, and a simulation is executed by numerically solving a mathematical model in which the relation between each variable is expressed by a state equation and an output equation by a computer. As a simulation time step when solving this state equation by a solution method such as the Euler method or the Runge-Kutta method, for example, 0.
A range such as 1 μs to 10 seconds is set, and the amount of change between the output value calculated by these solution methods and the old output value before simulation is evaluated based on a predetermined evaluation criterion,
If the amount of change as the evaluation result is small, the simulation time step is lengthened, and if the amount of change is large, the simulation time step is shortened and the simulation of the next step is executed sequentially. Is shortened. The output value of the simulator means a result obtained by executing the simulation. More specifically, when the state equation is solved, it means an internal state variable.

【0003】[0003]

【発明が解決しようとする課題】しかしながら、上述し
た従来のシミュレーション方法及びシミュレータでは、
入力値が急激に変化した場合などに、シミュレーション
時間刻みが不適切となり、また、その変更が間に合わ
ず、シミュレーション結果が実際の動作と齟齬を来した
り、シミュレーション自体が発散してしまう虞がある。
一方、これらの問題を解消するために、シミュレーショ
ン時間刻みをシミュレーション精度を確保できる小さな
値に予め固定しておく方法も実用化されているが、この
場合は、シミュレーションの実行時間が長くなり、実時
間でシミュレーションを実行する必要がある上記の運転
訓練シミュレータ等のシミュレータへの適用が困難とな
る。また、数学モデルの簡素化によってシミュレーショ
ンの高速化を図っても、数学モデルの簡素化によってシ
ミュレーション精度を悪化させる虞も生じる。本発明
は、かかる実情に鑑みてなされたものであり、その目的
は、上述の問題点を解消し、シミュレーションが高速に
実行可能で、然も、急激な入力変化に対して精度よく追
従可能なシミュレーション方法及びシミュレータを提供
する点にある。
However, in the conventional simulation method and simulator described above,
When the input value changes suddenly, the simulation time interval becomes inappropriate, the change cannot be made in time, the simulation result may be inconsistent with the actual operation, or the simulation itself may diverge. .
On the other hand, in order to solve these problems, a method of fixing the simulation time interval to a small value that can secure the simulation accuracy in advance has been put to practical use, but in this case, the simulation execution time becomes longer and the actual execution time becomes longer. It becomes difficult to apply to a simulator such as the above-mentioned driving training simulator which needs to execute a simulation in time. Further, even if the speed of the simulation is increased by simplifying the mathematical model, there is a possibility that the accuracy of the simulation is deteriorated by the simplified mathematical model. The present invention has been made in view of the above circumstances, and an object of the present invention is to solve the above-described problems, to execute a simulation at high speed, and to accurately follow a sudden input change. It is to provide a simulation method and a simulator.

【0004】[0004]

【課題を解決するための手段】この目的を達成するため
の本発明の第一の特徴構成は、シミュレータの出力値の
変化量に応じてシミュレーション時間刻みを変化させる
シミュレーション方法に関し、特許請求の範囲の欄の請
求項1に記載した通り、前記シミュレータへの1または
2以上の入力値に係るイベント発生直後のシミュレーシ
ョン時間刻みを決定する時間刻み決定手順をシミュレー
ション実行前に実行する点にある。
A first feature of the present invention for achieving this object is a simulation method for changing a simulation time interval according to a change amount of an output value of a simulator. As described in claim 1 in the column of (1), a time step determining procedure for determining a simulation time step immediately after the occurrence of an event relating to one or more input values to the simulator is performed before executing the simulation.

【0005】同第二の特徴構成は、上述の第一の特徴構
成のシミュレーション方法に関し、特許請求の範囲の欄
の請求項2に記載した通り、前記イベント発生に係る各
入力値の予め想定される最大変化量に対してシミュレー
ションが発散しない限界シミュレーション時間刻みのテ
ーブルを予め作成しておき、前記時間刻み決定手順にお
いて、イベント発生直前のシミュレーション実行時に設
定された補正前シミュレーション時間刻みを前記テーブ
ルの前記限界シミュレーション時間刻みに基づいて補正
してイベント発生直後のシミュレーション時間刻みとす
る点にある。
The second characteristic configuration relates to a simulation method of the above-mentioned first characteristic configuration, and as described in claim 2 of the claims, the input values relating to the event occurrence are assumed in advance. A table of the limit simulation time interval at which the simulation does not diverge for the maximum change amount is created in advance, and in the time interval determination procedure, the pre-correction simulation time interval set at the time of executing the simulation immediately before the event occurs is stored in the table. The point is that the correction is made based on the limit simulation time interval to make the simulation time interval immediately after the event occurs.

【0006】同第三の特徴構成は、上述の第二の特徴構
成のシミュレーション方法に関し、特許請求の範囲の欄
の請求項3に記載した通り、前記時間刻み決定手順にお
いて、イベント発生時に入力値が変化した入力別に、前
記限界シミュレーション時間刻みと前記補正前シミュレ
ーション時間刻みの差を実際の入力値変化量と予め想定
される最大変化量の比で按分したもの、または、それに
所定の安全係数を乗じたものを前記補正前シミュレーシ
ョン時間刻みに加算して補正後シミュレーション時間刻
みを演算し、イベント発生時に単一の入力で入力値が変
化した場合は、前記補正後シミュレーション時間刻みを
イベント発生直後のシミュレーション時間刻みとし、イ
ベント発生時に複数の入力で入力値が変化した場合は、
入力値が変化した入力別の前記補正後シミュレーション
時間刻みの最小値、または、その最小値に所定の安全係
数を乗じたものを算出してイベント発生直後のシミュレ
ーション時間刻みとする点にある。
The third characteristic configuration relates to a simulation method of the above-mentioned second characteristic configuration. As described in claim 3 of the claims, in the time step determining procedure, an input value is set when an event occurs. For each input that has changed, the difference between the limit simulation time step and the pre-correction simulation time step is proportionally divided by the ratio between the actual input value change amount and the maximum change amount assumed in advance, or a predetermined safety coefficient. The multiplied value is added to the pre-correction simulation time step to calculate the post-correction simulation time step.If the input value changes with a single input when an event occurs, the post-correction simulation time step is used immediately after the event occurrence. If the input value changes with multiple inputs at the time of the event,
The point is that the minimum value of the post-correction simulation time increment for each input whose input value has changed or the minimum value multiplied by a predetermined safety coefficient is calculated and used as the simulation time increment immediately after the event occurrence.

【0007】同第四の特徴構成は、上述の第一、第二、
または、第三の特徴構成のシミュレーション方法に関
し、特許請求の範囲の欄の請求項4に記載した通り、前
記イベントが入力値の変化の有無にかかわらず一定時間
間隔で発生する点にある。
[0007] The fourth characteristic configuration is the above-mentioned first, second, and third features.
Alternatively, with regard to the simulation method of the third characteristic configuration, as described in claim 4 in the claims section, the event occurs at a fixed time interval regardless of the presence or absence of a change in the input value.

【0008】同第五の特徴構成は、シミュレータの出力
値の変化量に応じてシミュレーション時間刻みを変化さ
せるシミュレーション時間刻み可変型のシミュレータに
関し、特許請求の範囲の欄の請求項5に記載した通り、
1または2以上の入力値を受信した直後のシミュレーシ
ョン時間刻みをシミュレーション実行前に決定する時間
刻み決定手段を備えている点にある。
The fifth characteristic configuration relates to a simulation time step variable type simulator that changes the simulation time step according to the amount of change in the output value of the simulator, as described in claim 5 of the claims. ,
The present invention is characterized in that a time interval determining means for determining a simulation time interval immediately after receiving one or more input values before executing the simulation is provided.

【0009】同第六の特徴構成は、上述の第五の特徴構
成のシミュレータに関し、特許請求の範囲の欄の請求項
6に記載した通り、入力値受信時の各入力値の予め想定
される最大変化量に対してシミュレーションが発散しな
い限界シミュレーション時間刻みのテーブルを予め記憶
した記憶手段を備え、前記時間刻み決定手段が、入力値
受信直前のシミュレーション実行時に設定された補正前
シミュレーション時間刻み一時的に記憶し、前記テーブ
ルの前記限界シミュレーション時間刻みに基づいてその
補正前シミュレーション時間刻みを補正して入力値受信
直後のシミュレーション時間刻みとして出力する点にあ
る。
The sixth characteristic configuration relates to the simulator having the above-mentioned fifth characteristic configuration. As described in claim 6 of the claims, each input value at the time of receiving the input value is assumed in advance. Storage means for preliminarily storing a table of a limit simulation time interval at which the simulation does not diverge with respect to the maximum change amount, wherein the time interval determination means temporarily sets a pre-correction simulation time interval set at the time of the simulation immediately before receiving the input value; In that the simulation time step before correction is corrected based on the limit simulation time step in the table, and is output as a simulation time step immediately after receiving the input value.

【0010】同第七の特徴構成は、上述の第六の特徴構
成のシミュレータに関し、特許請求の範囲の欄の請求項
7に記載した通り、前記時間刻み決定手段が、各入力値
毎に前記限界シミュレーション時間刻みと前記補正前シ
ミュレーション時間刻みの差を実際の入力値変化量と予
め想定される最大変化量の比で按分したもの、または、
それに所定の安全係数を乗じたものを前記補正前シミュ
レーション時間刻みに加算して補正後シミュレーション
時間刻みを演算する第1補正後時間刻み演算手段と、入
力値が複数の入力において変化した場合に、前記補正後
時間刻み演算手段の各入力値毎の演算結果の内の最小
値、または、その最小値に所定の安全係数を乗じて補正
後シミュレーション時間刻みとする第2補正後時間刻み
演算手段を備えてなる点にある。
The seventh feature configuration relates to the simulator having the sixth feature configuration described above, and as described in claim 7 in the claims section, the time step determining means sets the time interval for each input value. The difference between the limit simulation time step and the pre-correction simulation time step is proportionally divided by the ratio of the actual input value change amount and the previously assumed maximum change amount, or
A first corrected time step calculating means for calculating a post-correction simulation time step by adding a value obtained by multiplying the multiplied value by a predetermined safety coefficient to the pre-correction simulation time step, and when an input value changes in a plurality of inputs, A second corrected time step calculating means for calculating a corrected simulation time step by multiplying a minimum value or a minimum value of the calculation results for each input value of the post-correction time step calculating means by a predetermined safety coefficient; The point is to prepare.

【0011】同第八の特徴構成は、上述の第五、第六、
または、第七の特徴構成のシミュレータに関し、特許請
求の範囲の欄の請求項8に記載した通り、入力値の変化
の有無にかかわらず一定時間間隔で入力値を受信する点
にある。
The eighth characteristic configuration is the fifth, sixth,
Alternatively, regarding the simulator having the seventh characteristic configuration, as described in claim 8 in the claims section, the input value is received at fixed time intervals regardless of the presence or absence of a change in the input value.

【0012】以下に作用を説明する。第一または第五の
特徴構成によれば、入力値の変化の前後において、入力
値変化直前のシミュレーションの実行によって設定され
た入力値変化直後のシミュレーションに使用されるシミ
ュレーション時間刻みが、そのまま使用されるのではな
く、前記入力値の変化に応じて、そのシミュレーション
時間刻みが適切に補正されるか新たに適正値が与えられ
るために、極度のシミュレーション誤差の発生やシミュ
レーションの発散が回避でき、高精度のシミュレーショ
ンができるのである。然も、シミュレーション時間刻み
がシミュレータの出力値の変化量に応じて変化するシミ
ュレーション時間刻み可変型のシミュレータを維持でき
るため、シミュレーション実行時間が不必要に長くなる
のも回避でき、高精度と高速実行性の両立が図れるので
ある。
The operation will be described below. According to the first or fifth feature configuration, before and after the change of the input value, the simulation time interval used for the simulation immediately after the change of the input value set by executing the simulation immediately before the change of the input value is used as it is. Instead, the simulation time interval is appropriately corrected or a new appropriate value is given in accordance with the change in the input value, so that the occurrence of an extreme simulation error and the divergence of the simulation can be avoided. Simulation of accuracy can be performed. Of course, since the simulation time step can be maintained in a variable simulation time step type simulator in which the simulation time step changes according to the amount of change in the output value of the simulator, it is possible to avoid unnecessary increase in the simulation execution time, and to achieve high accuracy and high speed execution. Sex can be balanced.

【0013】第二または第六の特徴構成によれば、前記
テーブルには、各入力値の最大変化量に対応する限界シ
ミュレーション時間刻みがテーブル化されているため、
その限界シミュレーション時間刻みと、既に設定されて
いる入力値変化のない場合のシミュレーション時間刻み
とに基づいて、実際の入力値の変化量に応じた最適なシ
ミュレーション時間刻みを、入力値の変化量から設定す
ることができるのである。
According to the second or sixth characteristic configuration, the table stores the limit simulation time interval corresponding to the maximum change amount of each input value.
Based on the limit simulation time step and the previously set simulation time step when there is no change in the input value, the optimum simulation time step according to the actual input value change amount is calculated from the input value change amount. It can be set.

【0014】更に、第三または第七の特徴構成によれ
ば、第二または第六の特徴構成による実際の入力値の変
化量に応じた最適なシミュレーション時間刻みを簡単に
演算できるのである。従って、シミュレーション時間刻
みの補正に要する計算時間を短縮でき、総合的な実行処
理時間の短縮が図れるのである。
Further, according to the third or seventh feature configuration, it is possible to easily calculate the optimum simulation time interval according to the amount of change of the actual input value by the second or sixth feature configuration. Therefore, the calculation time required for correcting the simulation time increment can be shortened, and the overall execution processing time can be shortened.

【0015】第四または第八の特徴構成によれば、実シ
ステムにおいてシミュレーションの対象である特定のプ
ロセスが実システムにおける特定のコントローラの制御
対象であって、一定時間間隔で制御される場合に、前記
コントローラの制御量をシミュレーションの入力値とし
て、当該システムのシミュレーションに適用できるので
ある。例えば、実時間でシミュレーションを実行する必
要のあるプラント施設の運転者訓練用の運転訓練シミュ
レータを、プラント施設を制御するコントローラと、そ
のコントローラの制御対象をモデル化したシミュレータ
で構成するに当たり、本特徴構成が適用できるのであ
る。
According to the fourth or eighth aspect, when a specific process to be simulated in a real system is a control target of a specific controller in the real system and is controlled at a fixed time interval, The control amount of the controller can be applied to the simulation of the system as an input value of the simulation. For example, an operation training simulator for operator training of a plant facility that needs to execute a simulation in real time is composed of a controller that controls the plant facility and a simulator that models the control target of the controller. The configuration can be applied.

【0016】[0016]

【発明の実施の形態】以下に本発明の実施の形態を図面
に基づいて説明する。
Embodiments of the present invention will be described below with reference to the drawings.

【0017】図1に示すように、プラント施設の運転者
訓練用の運転訓練シミュレーションシステムは、本発明
に係るシミュレーション方法を実行するシミュレータで
ある運転訓練シミュレータ1とコントローラ2とこれら
を接続するシステムバス3で構成される。前記運転訓練
シミュレータ1は、実際のプラント施設で発生する各種
プロセスをシミュレーションするプロセスシミュレータ
4とそのプロセスシミュレータ4にシミュレーション時
間刻みの初期値を供給する時間刻み決定手順S2を実行
する時間刻み決定手段5で構成される。具体的には、前
記運転訓練シミュレータ1、並びに、前記プロセスシミ
ュレータ4及び前記時間刻み決定手段5は、本発明に係
るシミュレーション方法を実行するプログラムを実行可
能なエンジニアリングワークステーション等の汎用コン
ピュータ6で構成される。前記コントローラ2は実際の
プラント施設で使用されるコントローラと同じものが使
用されるが、前記汎用コンピュータ6で等価なものを構
成しても構わない。この場合、前記システムバス3は、
実質的に、前記汎用コンピュータ6内に形成される。
As shown in FIG. 1, an operation training simulation system for operator training of a plant facility includes an operation training simulator 1 which is a simulator for executing the simulation method according to the present invention, a controller 2, and a system bus connecting these. 3 The operation training simulator 1 includes a process simulator 4 for simulating various processes generated in an actual plant facility and a time step determining means 5 for executing a time step determining step S2 for supplying an initial value of a simulation time step to the process simulator 4. It consists of. Specifically, the driving training simulator 1, the process simulator 4 and the time step determining means 5 are configured by a general-purpose computer 6 such as an engineering workstation capable of executing a program for executing the simulation method according to the present invention. Is done. The controller 2 is the same as the controller used in the actual plant facility, but the general-purpose computer 6 may be equivalent. In this case, the system bus 3
Substantially, it is formed in the general-purpose computer 6.

【0018】前記システムバス3を介して、前記運転訓
練シミュレータ1と前記コントローラ2が1秒間隔で通
信を行う。この通信間隔は実システムに合わせて任意に
設定変更可能である。具体的には、図2に示すように、
前記通信間隔で設定される単位時間の一部分を前記運転
訓練シミュレータ1と前記コントローラ2間の通信時間
に割り当て、その通信時間内に前記プロセスシミュレー
タ4の出力値でありシミュレーション結果であるプロセ
ス値を前記コントローラ2へ出力し、反対に、前記コン
トローラ2はそのプロセス値に応じた操作量を前記プロ
セスシミュレータ4の入力値として出力し、前記プロセ
スシミュレータ4に対してイベント発生S1させる。ま
た、前記入力値は前記時間刻み決定手段5へも入力され
る。
The driving training simulator 1 and the controller 2 communicate with each other via the system bus 3 at intervals of one second. The communication interval can be arbitrarily changed according to the actual system. Specifically, as shown in FIG.
A part of the unit time set by the communication interval is allocated to a communication time between the driving training simulator 1 and the controller 2, and a process value which is an output value of the process simulator 4 and a simulation result is obtained within the communication time. The controller 2 outputs an operation amount corresponding to the process value as an input value of the process simulator 4, and causes the process simulator 4 to generate an event S1. Further, the input value is also input to the time step determining means 5.

【0019】図2に示すように、前記単位時間の内の前
記通信時間以降の部分は、主として、前記運転訓練シミ
ュレータ1がシミュレーションの実行に使用する。具体
的には、前記時間刻み決定手段5が、前記通信終了直後
に、前記コントローラ2からの入力値の変化量に応じ
て、前記イベント発生S1直後に新たに実行すべきシミ
ュレーション用のシミュレーション時間刻みを決定する
時間刻み決定手順S2を先ず実行し、引き続き、その時
間刻み決定手順S2を実行して得られたシミュレーショ
ン時間刻みで、前記プロセスシミュレータ4がシミュレ
ーションの実行S3を開始する。前記単位時間内でそれ
以降継続するシミュレーションの各シミュレーション時
間刻みは、従来通り、前記プロセスシミュレータ4が直
前のシミュレーション時に決定されたシミュレーション
時間刻みでシミュレーションを実行し、その出力値の変
化量に応じて次のシミュレーションのシミュレーション
時間刻みを決定する手順を繰り返し実行し、前記プロセ
スシミュレータ4の一連のシミュレーション実行S3が
終了する。尚、従来技術によるシミュレーション方法ま
たはシミュレータの場合の前記単位時間の使用区分を、
参考のために、図3に示す。
As shown in FIG. 2, the portion after the communication time in the unit time is mainly used by the driving training simulator 1 for executing a simulation. More specifically, the time step determining means 5 performs a simulation time step for a simulation to be newly executed immediately after the event occurrence S1 in accordance with the change amount of the input value from the controller 2 immediately after the end of the communication. Is executed first, and then the process simulator 4 starts the execution S3 of the simulation at the simulation time interval obtained by executing the time interval determination procedure S2. For each simulation time step of the simulation that continues thereafter within the unit time, the process simulator 4 executes the simulation at the simulation time step determined at the time of the immediately preceding simulation as in the related art, and according to the amount of change in the output value. The procedure for determining the simulation time interval of the next simulation is repeatedly executed, and a series of simulation execution S3 of the process simulator 4 is completed. In the case of the simulation method or the simulator according to the prior art, the use division of the unit time is
It is shown in FIG. 3 for reference.

【0020】以下、前記プロセスシミュレータ4による
シミュレーション時間刻みの決定方法につき説明する。
前記プロセスシミュレータ4のシミュレーション対象の
プロセスが数式1に示す状態方程式と数式2に示す出力
方程式で構成される数式モデルで表現される。また、こ
の状態方程式の解法手法としてオイラー法を採用する。
Hereinafter, a method of determining the simulation time interval by the process simulator 4 will be described.
The process to be simulated by the process simulator 4 is represented by a mathematical model composed of a state equation represented by Equation 1 and an output equation represented by Equation 2. The Euler method is adopted as a method of solving the equation of state.

【0021】[0021]

【数1】dX=A*X+B*U## EQU1 ## dX = A * X + B * U

【0022】[0022]

【数2】Y=C*X+D*U## EQU2 ## Y = C * X + D * U

【0023】但し、Xはi番目の状態変数xi を成分と
する内部状態ベクトル、dXはi番目の状態変数の時間
微分値dxi を成分とする内部状態微分値ベクトル、U
はi番目の入力変数ui を成分とする入力ベクトル、Y
はi番目の出力変数yi を成分とする出力ベクトル、A
はシステム行列、Bは制御行列、Cは観測行列、Dは直
達行列を示す。オイラー法の場合、前記微分dxi が前
進差分(xi (t+h)−xi (t))/hで置換され
る。このときのhがシミュレーション時間刻みである。
[0023] However, X is i-th internal state vector with state variables x i component, dX internal state differential value vector with i-th time derivative dx i components of the state variable, U
Input vector to the i-th input variable u i component, Y
Is an output vector having the i-th output variable y i as a component, A
Denotes a system matrix, B denotes a control matrix, C denotes an observation matrix, and D denotes a direct matrix. For Euler, the differential dx i is replaced by forward difference (x i (t + h) -x i (t)) / h. H at this time is a simulation time interval.

【0024】次のシミュレーションのシミュレーション
時間刻みhNEXTは、現シミュレーション結果の前記内部
状態微分値ベクトルdXの各状態変数の微分値dxi
つき、数式3に示す評価値ei を全て算出し、数式4の
判定処理を実行し、数式4を満足した場合は、現シミュ
レーションのシミュレーション時間刻みhを半分に短縮
し、数式4を満足せず数式5を満足した場合は、前記シ
ミュレーション時間刻みhを2倍にして導出される。数
式4及び数式5を満足しない場合は、シミュレーション
時間刻みhNEXTは現シミュレーションのシミュレーショ
ン時間刻みhを維持する。尚、下記の数式3乃至5にお
いて、dxi-OLD は前回のシミュレーション時の微分値
dxi であり、Tは事前に設定される定数で、シミュレ
ーション精度を調整する機能を有する。また、max
(E)は評価値ei を成分とする評価値ベクトルEの成
分の内の最大値を選択する関数値である。
[0024] The next simulation simulation time step h NEXT, per differential value dx i of each state variable of the internal state differential value vector dX of the current simulation results to calculate all the evaluation value e i shown in equation 3, equation When the determination process of Equation 4 is performed and Expression 4 is satisfied, the simulation time interval h of the current simulation is reduced to half, and when Expression 4 is not satisfied and Expression 5 is satisfied, the simulation time interval h is set to 2 It is derived by doubling. If Equations 4 and 5 are not satisfied, the simulation time step h NEXT maintains the simulation time step h of the current simulation. Note that in Equation 3 to 5 below, dx i-OLD is the differential value dx i at the previous simulation, T is a constant that is set in advance, has a function of adjusting the simulation accuracy. Also, max
(E) is a function value to select the maximum value of the components of the evaluation value vector E to component evaluation value e i.

【0025】[0025]

【数3】ei =(dxi −dxi-OLD )/(|dx
i-OLD |+T)
[Equation 3] e i = (dx i −dx i-OLD ) / (| dx
i-OLD | + T)

【0026】[0026]

【数4】h*max(E)>10*T## EQU4 ## h * max (E)> 10 * T

【0027】[0027]

【数5】h*max(E)<TH * max (E) <T

【0028】次に、前記時間刻み決定手段5によるシミ
ュレーション時間刻みの決定方法につき説明する。前記
時間刻み決定手段5は、前記コントローラ2からの各入
力値ui の最大変化量ΔuMAXiと、当該最大変化時にシ
ミュレーションが発散しない限界の限界シミュレーショ
ン時間刻みhXiを各別に予め測定しておき、その限界シ
ミュレーション時間刻みhXiをテーブル化して、前記汎
用コンピュータ6に設けられた記憶手段7に格納されて
いる。
Next, a method of determining the simulation time step by the time step determining means 5 will be described. Said time increment determining means 5, and the maximum variation Delta] u MAXi of each input value u i from the controller 2, the limit simulation time step h Xi limit the simulation when the maximum change does not diverge measured in advance to each other The limit simulation time increment h Xi is stored in a table in the storage means 7 provided in the general-purpose computer 6.

【0029】前記入力値ui の内、前記コントローラ2
との通信時に実際に変化したものについてのみ、各別に
数式6を演算して入力値ui 別のシミュレーション時間
刻みhi を算出する。但し、数式6において、h0iは前
記コントローラ2との通信直前のシミュレーションで数
式1乃至5の手順で決定されたシミュレーション時間刻
みであり、Δui は前記コントローラ2との通信時にお
ける実際の入力値uiの変化量であり、a1 は安全係数
であり1以上の値である。
Of the input values u i , the controller 2
Only those actually changed during communication with, and calculates an input value u i by the simulation time step h i to each other by calculating Equation 6. Here, in Expression 6, h 0i is a simulation time interval determined by the procedure of Expressions 1 to 5 in a simulation immediately before communication with the controller 2, and Δu i is an actual input value at the time of communication with the controller 2. a variation of u i, a 1 is one or more valued in the safety factor.

【0030】[0030]

【数6】hi =h0i+((hXi−h0i)*(Δui /Δ
MAXi))*a1
H i = h 0i + ((h Xi −h 0i ) * (Δu i / Δ
u MAXi )) * a 1

【0031】前記入力値ui の内、前記コントローラ2
との通信時に実際に変化したものが1入力の場合は、数
式6のhi をシミュレーション時間刻みhとする。2入
力以上が変化した場合は、更に、数式7に基づいて,シ
ミュレーション時間刻みhを決定する。但し、数式7に
おいて、min(H)は入力値ui 別のシミュレーショ
ン時間刻みhi を成分とするシミュレーション時間刻み
ベクトルHの成分の内の最小値を選択する関数値であ
り、a2 は安全係数であり1以下の値である。。
[0031] Of the input values u i, the controller 2
If the input that has actually changed at the time of communication with is one input, let h i in Equation 6 be the simulation time interval h. If two or more inputs change, the simulation time increment h is further determined based on Equation 7. However, in Equation 7, min (H) is a function value for selecting the minimum value among the components of the simulation time step vector H having the simulation time step h i as a component for each input value u i , and a 2 is a safe value. The coefficient is a value of 1 or less. .

【0032】[0032]

【数7】h=min(H)*a2 H = min (H) * a 2

【0033】尚、数式6または数式7において、前記安
全係数a1 またはa2 が1に固定される場合は夫々の乗
算は不要である。
Note that in Equations 6 and 7, when the safety coefficient a 1 or a 2 is fixed at 1, each multiplication is unnecessary.

【0034】以下に、別実施形態を説明する。 〈1〉上記の実施の形態では、前記通信間隔は一定時間
間隔で発生し、必ずしも前記コントローラ2が出力する
前記プロセスシミュレータ4の入力値は変化しない構成
となっているが、前記入力値の少なくとも一つが変化し
た場合のみ前記運転訓練シミュレータ1と前記コントロ
ーラ2間の通信を行うように不定間隔でイベントを発生
させても構わない。 〈2〉前記プロセスシミュレータ4に組み込まれた数式
モデルは、必ずしも数式1及び2の形式で表現されなく
ても構わない。また、状態方程式の解法はオイラー法以
外でも構わない。例えば、ルンゲクッタ法を使用するの
も好ましい。 〈3〉前記コントローラ2からの各入力値ui の最大変
化量ΔuMAXiに対応する限界シミュレーション時間刻み
Xiと、通信直前のシミュレーションで決定されたシミ
ュレーション時間刻みh0iからシミュレーション時間刻
みhを決定する手法は、必ずしも数式6と数式7によら
なくても構わない。例えば、前記入力値u i の特性に応
じて、数式6のような線型補間を施さずに、非線形の補
間を行っても構わない。また、前記安全係数a1 を全て
の入力値ui に対して一定値とせず、前記入力値ui
に変化させても構わない。更に、前記安全係数a2 も実
際に変化した入力値ui または入力数に対応して変化さ
せても構わない。 〈4〉前記時間刻み決定手順S2は、必ずしも、各入力
値ui の最大変化量Δu MAXi別の前記限界シミュレーシ
ョン時間刻みhXiをテーブル化して使用しなくても構わ
ない。例えば、各入力値ui の最大変化時の入力パター
ン別に限界シミュレーション時間刻みを測定して予めテ
ーブル化しても構わない。また、各入力値ui の変化量
から直接シミュレーション時間刻みhをニューラルネッ
トワーク等を応用して決定できる構成としても構わな
い。
Hereinafter, another embodiment will be described. <1> In the above embodiment, the communication interval is a fixed time
Occurs at intervals and necessarily output by the controller 2
Configuration in which the input value of the process simulator 4 does not change
But at least one of the input values changes
The driving training simulator 1 and the controller
Event is generated at irregular intervals so that communication between
You may let it. <2> Formulas incorporated in the process simulator 4
The model is not necessarily expressed in the form of Equations 1 and 2.
It does not matter. The solution of the equation of state is less than Euler's method.
It can be outside. For example, using the Runge-Kutta method
Is also preferred. <3> Each input value u from the controller 2iMaximum change
Amount ΔuMAXiLimit simulation time step corresponding to
hXiAnd the stain determined in the simulation immediately before communication.
Hour time increment0iFrom simulation time
The method for determining h is not necessarily based on Equations 6 and 7.
You don't have to. For example, the input value u iAccording to the characteristics of
Therefore, without performing the linear interpolation as in Equation 6, the nonlinear interpolation is performed.
You can go between them. Also, the safety coefficient a1All
Input value uiTo the input value uiAnother
It may be changed to. Further, the safety factor aTwoAlso real
Input value u that has changediOr change according to the number of inputs
You can do it. <4> The time step determination step S2 is not necessarily performed by each input.
Value uiMaximum change amount Δu MAXiAnother said limit simulation
Time intervalXiYou don't have to use
Absent. For example, each input value uiInput pattern when the maximum change of
Measure the limit simulation time interval for each
May be used. Also, each input value uiAmount of change
Simulation time increment h directly from the neural network
It may be a configuration that can be determined by applying
No.

【0035】〈5〉上記の実施形態では、前記プロセス
シミュレータ4及び前記時間刻み決定手段5は、本発明
に係るシミュレーション方法を実行するプログラムを実
行可能なエンジニアリングワークステーション等の汎用
コンピュータ6で構成されずに、専用プロセッサで各別
にまたは一体で構成されても構わない。
<5> In the above embodiment, the process simulator 4 and the time step determining means 5 are constituted by a general-purpose computer 6 such as an engineering workstation capable of executing a program for executing the simulation method according to the present invention. Instead, a dedicated processor may be used separately or integrally.

【0036】〈6〉上記の実施形態では、プラント施設
の運転者訓練用の運転訓練シミュレーションシステムを
例に、前記コントローラ2が存在するシステムを想定し
たが、必ずしも前記コントローラ2が存在する必要がな
く、前記数式モデル自体にフィードバック系統が含まれ
てあっても構わない。
<6> In the above-described embodiment, a system in which the controller 2 exists is assumed by taking an operation training simulation system for operator training of a plant facility as an example. However, the controller 2 does not necessarily need to exist. Alternatively, the mathematical expression model itself may include a feedback system.

【0037】[0037]

【発明の効果】以上説明したように、本発明によれば、
シミュレータに対する入力値が急激に変化しても、前記
時間刻み決定手段または前記時間刻み決定手順の実行に
よりシミュレーション時間刻みが適切に設定されシミュ
レーションを実行できるため、それ以降の入力値が変化
しないシミュレーションでは、従来通りに、シミュレー
タの出力値の変化量に応じてシミュレーション時間刻み
を変化させるシミュレーション方法が使用できた。図4
に、本実施形態によるシミュレーション結果と図3に示
す従来手法でのシミュレーション結果との出力値の差を
例示する。結果として、本発明によれば、シミュレーシ
ョンが高速に実行可能で、然も、急激な入力変化に対し
て精度よく追従可能なシミュレーション方法及びシミュ
レータを提供することができた。更に、イベント発生の
時間間隔が一定時間に固定されている場合は、シミュレ
ーション時間の短縮効果によって、シミュレーション実
行後の空き時間を、シミュレーションの実行にかかるコ
ンピュータを他の機能の実行に有効利用できるという効
果も期待できるのである。
As described above, according to the present invention,
Even if the input value to the simulator changes abruptly, the simulation can be executed with the simulation time interval appropriately set by the execution of the time interval determination means or the time interval determination procedure. As in the past, a simulation method that changes the simulation time interval according to the amount of change in the output value of the simulator could be used. FIG.
FIG. 3 illustrates the difference between the output value of the simulation result according to the present embodiment and the simulation result obtained by the conventional method shown in FIG. As a result, according to the present invention, it is possible to provide a simulation method and a simulator that can execute a simulation at a high speed and can accurately follow a sudden input change. Furthermore, when the time interval of the event occurrence is fixed to a fixed time, the effect of shortening the simulation time is that the idle time after the execution of the simulation can be effectively used by the computer for executing the simulation for executing other functions. The effect can be expected.

【0038】尚、特許請求の範囲の項に、図面との対照
を便利にするために符号を記すが、該記入により本発明
は添付図面の構成に限定されるものではない。
In the claims, reference numerals are provided for convenience of comparison with the drawings, but the present invention is not limited to the configuration shown in the attached drawings.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明に係るシミュレータのブロック構成図FIG. 1 is a block diagram of a simulator according to the present invention.

【図2】 本発明に係る時間刻み決定手順とシミュレーシ
ョン実行とイベント発生の時間関係を示すタイミング図
FIG. 2 Time step determination procedure and simulation according to the present invention
Diagram showing the time relationship between the action execution and the event occurrence

【図3】従来技術におけるシミュレーション実行とイベ
ント発生の時間関係を示すタイミング図
FIG. 3 is a timing chart showing a time relationship between a simulation execution and an event occurrence in the related art.

【図4】本発明に係るシミュレーション方法と従来のシ
ミュレーション方法とのシミュレーション結果の比較図
FIG. 4 is a comparison diagram of simulation results between a simulation method according to the present invention and a conventional simulation method.

【符号の説明】[Explanation of symbols]

1 運転訓練シミュレータ 2 コントローラ 3 システムバス 4 プロセスシミュレータ 5 時間刻み決定手段 6 汎用コンピュータ 7 記憶手段 S1 イベント発生 S2 時間刻み決定手順 S3 シミュレーション実行 DESCRIPTION OF SYMBOLS 1 Driving training simulator 2 Controller 3 System bus 4 Process simulator 5 Time step determining means 6 General-purpose computer 7 Storage means S1 Event occurrence S2 Time step determining procedure S3 Simulation execution

Claims (8)

【特許請求の範囲】[Claims] 【請求項1】 シミュレータの出力値の変化量に応じて
シミュレーション時間刻みを変化させるシミュレーショ
ン方法において、前記シミュレータへの1または2以上
の入力値に係るイベント発生(S1)直後のシミュレー
ション時間刻みを決定する時間刻み決定手順(S2)を
シミュレーション実行(S3)前に実行することを特徴
とするシミュレーション方法。
1. A simulation method for changing a simulation time step according to a change amount of an output value of a simulator, wherein the simulation time step is determined immediately after an event (S1) relating to one or more input values to the simulator. A simulation method characterized in that a time step determining step (S2) to be performed is executed before executing a simulation (S3).
【請求項2】 前記イベント発生(S1)に係る各入力
値の予め想定される最大変化量に対してシミュレーショ
ンが発散しない限界シミュレーション時間刻みのテーブ
ルを予め作成しておき、前記時間刻み決定手順(S2)
において、イベント発生(S1)直前のシミュレーショ
ン実行(S3)時に設定された補正前シミュレーション
時間刻みを前記テーブルの前記限界シミュレーション時
間刻みに基づいて補正してイベント発生(S1)直後の
シミュレーション時間刻みとする請求項1記載のシミュ
レーション方法。
2. A limit simulation time interval table in which the simulation does not diverge for a previously assumed maximum change amount of each input value related to the event occurrence (S1) is created in advance, and the time interval determination procedure ( S2)
In the above, the simulation time interval before correction set at the time of the simulation execution (S3) immediately before the event occurrence (S1) is corrected based on the limit simulation time interval of the table to be the simulation time interval immediately after the event occurrence (S1). The simulation method according to claim 1.
【請求項3】 前記時間刻み決定手順(S2)におい
て、イベント発生(S1)時に入力値が変化した入力別
に、前記限界シミュレーション時間刻みと前記補正前シ
ミュレーション時間刻みの差を実際の入力値変化量と予
め想定される最大変化量の比で按分したもの、または、
それに所定の安全係数を乗じたものを前記補正前シミュ
レーション時間刻みに加算して補正後シミュレーション
時間刻みを演算し 前記イベント発生(S1)時に単一
の入力で入力値が変化した場合は、前記補正後シミュレ
ーション時間刻みを前記イベント発生(S1)直後のシ
ミュレーション時間刻みとし、 前記イベント発生(S1)時に複数の入力で入力値が変
化した場合は、入力値が変化した入力別の前記補正後シ
ミュレーション時間刻みの最小値、または、その最小値
に所定の安全係数を乗じたものを算出して前記イベント
発生(S1)直後のシミュレーション時間刻みとする請
求項2記載のシミュレーション方法。
3. In the time step determining step (S2), a difference between the limit simulation time step and the pre-correction simulation time step is calculated for each input whose input value has changed at the time of the event occurrence (S1). And proportionally divided by the ratio of the maximum change assumed in advance, or
A value obtained by multiplying the calculated value by a predetermined safety coefficient is added to the simulation time step before correction to calculate the simulation time step after correction. If the input value changes with a single input when the event occurs (S1), the correction is performed. The post-simulation time step is a simulation time step immediately after the event occurrence (S1). If the input value changes by a plurality of inputs at the time of the event occurrence (S1), the post-correction simulation time for each input whose input value has changed 3. The simulation method according to claim 2, wherein a minimum value of the step or a value obtained by multiplying the minimum value by a predetermined safety coefficient is calculated, and the calculated value is used as a simulation time step immediately after the occurrence of the event (S1).
【請求項4】 前記イベントが入力値の変化の有無にか
かわらず一定時間間隔で発生する請求項1、2または3
記載のシミュレーション方法。
4. The event according to claim 1, wherein the event occurs at a fixed time interval regardless of whether or not the input value has changed.
The described simulation method.
【請求項5】 シミュレータの出力値の変化量に応じて
シミュレーション時間刻みを変化させるシミュレーショ
ン時間刻み可変型のシミュレータにおいて、1または2
以上の入力値を受信した直後のシミュレーション時間刻
みをシミュレーション実行(S3)前に決定する時間刻
み決定手段(5)を備えたシミュレータ。
5. A simulation time step variable type simulator that changes a simulation time step according to a change amount of an output value of the simulator.
A simulator provided with time step determining means (5) for determining a simulation time step immediately before receiving the above input value before executing the simulation (S3).
【請求項6】 入力値受信時の各入力値の予め想定され
る最大変化量に対してシミュレーションが発散しない限
界シミュレーション時間刻みのテーブルを予め記憶した
記憶手段(7)を備え、前記時間刻み決定手段(5)
が、入力値受信直前のシミュレーション実行(S3)時
に設定された補正前シミュレーション時間刻み一時的に
記憶し、前記テーブルの前記限界シミュレーション時間
刻みに基づいてその補正前シミュレーション時間刻みを
補正して入力値受信直後のシミュレーション時間刻みと
して出力する請求項5記載のシミュレータ。
6. A storage means (7) for storing in advance a limit simulation time step table in which the simulation does not diverge with respect to a previously assumed maximum change amount of each input value at the time of receiving the input value, and the time step determination is performed. Means (5)
Temporarily stores the pre-correction simulation time step set at the time of the simulation execution (S3) immediately before receiving the input value, and corrects the pre-correction simulation time step based on the limit simulation time step in the table to input the input value. The simulator according to claim 5, wherein the simulator outputs the simulation time interval immediately after reception.
【請求項7】 前記時間刻み決定手段(5)が、各入力
値毎に前記限界シミュレーション時間刻みと前記補正前
シミュレーション時間刻みの差を実際の入力値変化量と
予め想定される最大変化量の比で按分したもの、また
は、それに所定の安全係数を乗じたものを前記補正前シ
ミュレーション時間刻みに加算して補正後シミュレーシ
ョン時間刻みを演算する第1補正後時間刻み演算手段
と、入力値が複数の入力において変化した場合に、前記
補正後時間刻み演算手段の各入力値毎の演算結果の内の
最小値、または、その最小値に所定の安全係数を乗じて
補正後シミュレーション時間刻みとする第2補正後時間
刻み演算手段を備えてなる請求項6記載のシミュレー
タ。
7. The time step determining means (5) calculates, for each input value, a difference between the limit simulation time step and the pre-correction simulation time step between an actual input value change amount and a previously assumed maximum change amount. A first corrected time step calculating means for calculating the corrected simulation time step by adding a value proportionally divided by the ratio or a value obtained by multiplying the ratio by a predetermined safety coefficient to the pre-correction simulation time step; In the case of a change in the input of (1), a minimum value of the calculation results for each input value of the post-correction time step calculation means or a minimum safety value is multiplied by a predetermined safety coefficient to obtain a post-correction simulation time step. 7. The simulator according to claim 6, comprising a time step calculating means after two corrections.
【請求項8】 入力値の変化の有無にかかわらず一定時
間間隔で入力値を受信する請求項5、6または7記載の
シミュレータ。
8. The simulator according to claim 5, wherein the input value is received at a fixed time interval regardless of whether the input value has changed.
JP16260697A 1997-06-19 1997-06-19 Simulation method, and simulator Pending JPH1115362A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP16260697A JPH1115362A (en) 1997-06-19 1997-06-19 Simulation method, and simulator

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP16260697A JPH1115362A (en) 1997-06-19 1997-06-19 Simulation method, and simulator

Publications (1)

Publication Number Publication Date
JPH1115362A true JPH1115362A (en) 1999-01-22

Family

ID=15757798

Family Applications (1)

Application Number Title Priority Date Filing Date
JP16260697A Pending JPH1115362A (en) 1997-06-19 1997-06-19 Simulation method, and simulator

Country Status (1)

Country Link
JP (1) JPH1115362A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7392109B2 (en) 2000-07-31 2008-06-24 Kabushiki Kaisha Toyota Chuokenkyusho System for integrally generating NC data
US8990786B2 (en) 2010-12-15 2015-03-24 International Business Machines Corporation Program optimizing apparatus, program optimizing method, and program optimizing article of manufacture
JP2016162270A (en) * 2015-03-03 2016-09-05 三菱電機株式会社 Proximity detection device and proximity detection method

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7392109B2 (en) 2000-07-31 2008-06-24 Kabushiki Kaisha Toyota Chuokenkyusho System for integrally generating NC data
US8990786B2 (en) 2010-12-15 2015-03-24 International Business Machines Corporation Program optimizing apparatus, program optimizing method, and program optimizing article of manufacture
JP2016162270A (en) * 2015-03-03 2016-09-05 三菱電機株式会社 Proximity detection device and proximity detection method

Similar Documents

Publication Publication Date Title
US7346403B2 (en) Pid parameter adjustment device
US4918643A (en) Method and apparatus for substantially improving the throughput of circuit simulators
CN109100990B (en) Control device, recording medium, and control system
JP5227254B2 (en) Real-time calculation method and simulator of state quantity of process model
EP3385796A1 (en) Control device
CN108170029B (en) Parameter self-tuning method of MIMO full-format model-free controller based on partial derivative information
KR20060111400A (en) Control object model generation apparatus and generation method
JP2019507405A (en) How to configure collaborative simulation for integrated systems
JPH05127706A (en) Neural network type simulator
JPH1115362A (en) Simulation method, and simulator
Kersten et al. Application-based discussion of verified simulations of interval enclosure techniques
JPH0756608A (en) Signal processing method
CN110705159A (en) Heat source model parameter solving method, device, equipment and storage medium
JPH08110802A (en) Pid controller
Joseph et al. Homotopy based error estimator for adaptive modified euler in transient stability simulations
ATE532050T1 (en) METHOD FOR DETERMINING STRESS-STRAIN CURVES USING SPLINE INTERPOLATION BASED ON CHARACTERISTIC POINTS AND USING NEURONAL NETWORKS
JP3340923B2 (en) SAC controller
CN108008634B (en) Parameter self-tuning method of MISO partial-format model-free controller based on partial derivative information
Elnagar Optimal control computation for integro‐differential aerodynamic equations
JPH117314A (en) Plant simulator
Thielemans et al. Offline Learning Control to Improve the Accuracy of Real-time Explicit Co-simulation
JP3254618B2 (en) controller
CN118052183A (en) Numerical integration method and system for nonlinear circuit electromagnetic transient simulation
Rios-Avila f able: Estimation of marginal effects for models with alternative variable transformations
JP2660976B2 (en) Digital integration module for sampling controller.