JPS61110203A - 数式処理制御システム - Google Patents
数式処理制御システムInfo
- Publication number
- JPS61110203A JPS61110203A JP59230185A JP23018584A JPS61110203A JP S61110203 A JPS61110203 A JP S61110203A JP 59230185 A JP59230185 A JP 59230185A JP 23018584 A JP23018584 A JP 23018584A JP S61110203 A JPS61110203 A JP S61110203A
- Authority
- JP
- Japan
- Prior art keywords
- formula
- program
- mathematical
- processing
- mathematical formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F9/00—Arrangements for program control, e.g. control units
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/0265—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion
- G05B13/028—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion using expert systems only
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S706/00—Data processing: artificial intelligence
- Y10S706/902—Application using ai with detail of the ai system
- Y10S706/903—Control
- Y10S706/906—Process plant
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S706/00—Data processing: artificial intelligence
- Y10S706/902—Application using ai with detail of the ai system
- Y10S706/919—Designing, planning, programming, CAD, CASE
- Y10S706/922—Computer program preparation
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔発明の利用分野〕
本発明は数式処理による制御に係り、特に使用条件を満
足する数式を選択しプログラムに変換制御をおこなう数
式処理制御システムを提供することにある。
足する数式を選択しプログラムに変換制御をおこなう数
式処理制御システムを提供することにある。
計算機制御システムや情報処理システムでは、対象の数
式モデルをプログラムとして表現し、対象の将来像を予
測したり、直接計測できない状態骨を推定する方法が一
般的に用いられている。
式モデルをプログラムとして表現し、対象の将来像を予
測したり、直接計測できない状態骨を推定する方法が一
般的に用いられている。
これらの数式にも使用上の条件があって、どのような場
合でも使えるということはない。このたの、細組かの数
式を用意し、条件によって使い分けるということもめず
らしくかい。一方、当初考えていないような条件がみつ
かったり、適用範囲を広げたりする場合には、新たに数
式モデルの見直しをおこなってプログラムの修正、追加
を行なうととＫなる。しかしながらこの過程で様々な間
違いが入いりこみ、システムトラブルが発生することが
多い。この間違いを防ぐために細心の注意と多大の労力
が払われるが、人手に頼る現状のやり方では完全を期す
ことは困難である。
合でも使えるということはない。このたの、細組かの数
式を用意し、条件によって使い分けるということもめず
らしくかい。一方、当初考えていないような条件がみつ
かったり、適用範囲を広げたりする場合には、新たに数
式モデルの見直しをおこなってプログラムの修正、追加
を行なうととＫなる。しかしながらこの過程で様々な間
違いが入いりこみ、システムトラブルが発生することが
多い。この間違いを防ぐために細心の注意と多大の労力
が払われるが、人手に頼る現状のやり方では完全を期す
ことは困難である。
例えば加熱炉燃焼制御では、与えられたいくつかの制御
目標を同時に満足するような制御を実現する必要がある
。例えば、抽出時の鋼片温度を目標範囲におさめること
、燃料消費量を最小にすること、炉温の変化速度を一定
値以下にすること等様々の制御目標がある。一方これら
の目標は必らずしも固定したものでかい。
目標を同時に満足するような制御を実現する必要がある
。例えば、抽出時の鋼片温度を目標範囲におさめること
、燃料消費量を最小にすること、炉温の変化速度を一定
値以下にすること等様々の制御目標がある。一方これら
の目標は必らずしも固定したものでかい。
加熱炉の炉内現象は完全にはわかっていないため、制御
目標が経験的に決められていることが多い。したがって
操業Ｐ、駒を積むにしたがい、よりよい制御目標が見つ
かり修正する必要が生じてくる。また操業方針が変更に
なったり、当初計画に含まれていない種類の偉材を加熱
するような場合にも、制御目標の修正、追加が必要にな
る。これに伴ってプログラムの修正、追加をおこなうこ
とになるが、これには次のような問題がともなう。
目標が経験的に決められていることが多い。したがって
操業Ｐ、駒を積むにしたがい、よりよい制御目標が見つ
かり修正する必要が生じてくる。また操業方針が変更に
なったり、当初計画に含まれていない種類の偉材を加熱
するような場合にも、制御目標の修正、追加が必要にな
る。これに伴ってプログラムの修正、追加をおこなうこ
とになるが、これには次のような問題がともなう。
（１）修正、追加プログラムにひそむバグによるシステ
ムトラブル。
ムトラブル。
（ＩＩ） （ＩＩｌ−正、追加に多大の労力を必要と
する。
する。
そこでこれらの間鴨点をもう少し明確にするために、従
来システムにおいて制御目標がどのような形で用いられ
ていたかを分析してみる。
来システムにおいて制御目標がどのような形で用いられ
ていたかを分析してみる。
ＩＩｌ図は圧延用鋼片を加熱する多帯式加熱炉の燃焼制
御システムである（特開昭５８−２１０１２１）第２図
において１０は加熱炉、１１は鋼片、２０は燃焼制御シ
ステＡ、３０ａ、３０ｂ、３００は比例積分型調節器、
２１は炉内鋼片の移動スケジュールを予測するスケジュ
ール予測部、２２は炉内鋼片の最適昇温曲線の予測計算
部、２３は鋼片内部の温度を炉温計算値から計算する推
定部、２４は炉温の設定値を決定する演算部を表わして
いる。
御システムである（特開昭５８−２１０１２１）第２図
において１０は加熱炉、１１は鋼片、２０は燃焼制御シ
ステＡ、３０ａ、３０ｂ、３００は比例積分型調節器、
２１は炉内鋼片の移動スケジュールを予測するスケジュ
ール予測部、２２は炉内鋼片の最適昇温曲線の予測計算
部、２３は鋼片内部の温度を炉温計算値から計算する推
定部、２４は炉温の設定値を決定する演算部を表わして
いる。
予測計算部２２は複数の制御目標を満足する鋼片昇温曲
線を計算する。次の（１）〜（５）式は代表的な制御目
標を数式で記述したものである。
線を計算する。次の（１）〜（５）式は代表的な制御目
標を数式で記述したものである。
Ｊ＝Σｖ（ｉ）Δτ−→最小 ・・・・・
・（１）７Ｌ≦′″５ｏｒｒ！≦５０
・・・・・・（２）Δθし≦Δθ。。、≦Δθυ
・・・・・・（３）０≦ｖ（１）≦ｖ”（
ｉ）（１＝１〜ｎ） ””（４）ＴＬ（ｉ）≦Ｔ
（１）≦Ｔ ｏ（１） （Ｈ＝ １〜ｎ ） ・”
・・・（５）ここでｎは現時刻から抽出までの予定時間
（これは移動スケジュール予測部２１で計算される）を
ｎ分割することを表わし、ｉはその第１番目の時間ゾー
ンを、Δτ亀は第１時刻ゾーンの時間長を意味する。従
ってＴ（１）、 ｖ（ｉ）は第１時間ゾーンにおける炉
温、および計算対集−片の加熱に要する燃料流量を表わ
している。またθｏａｔ 、Δθｏｖ〒は抽出時の鋼片
温度の平均値および均熱度（鋼片表面と内部の温度差）
を表わしている。また添字りとＵのついた項はそれぞれ
の式の下限値と上限値を表わしそいる。（１）式は現在
から抽出までの加熱に要する燃料を最小化する評価関数
、（２）、 （ａ）式は細片抽出温度に関する制約式、
（４）、 （５）式は設備能力からくる燃料流量と炉温
の制約式を表わしている。
・（１）７Ｌ≦′″５ｏｒｒ！≦５０
・・・・・・（２）Δθし≦Δθ。。、≦Δθυ
・・・・・・（３）０≦ｖ（１）≦ｖ”（
ｉ）（１＝１〜ｎ） ””（４）ＴＬ（ｉ）≦Ｔ
（１）≦Ｔ ｏ（１） （Ｈ＝ １〜ｎ ） ・”
・・・（５）ここでｎは現時刻から抽出までの予定時間
（これは移動スケジュール予測部２１で計算される）を
ｎ分割することを表わし、ｉはその第１番目の時間ゾー
ンを、Δτ亀は第１時刻ゾーンの時間長を意味する。従
ってＴ（１）、 ｖ（ｉ）は第１時間ゾーンにおける炉
温、および計算対集−片の加熱に要する燃料流量を表わ
している。またθｏａｔ 、Δθｏｖ〒は抽出時の鋼片
温度の平均値および均熱度（鋼片表面と内部の温度差）
を表わしている。また添字りとＵのついた項はそれぞれ
の式の下限値と上限値を表わしそいる。（１）式は現在
から抽出までの加熱に要する燃料を最小化する評価関数
、（２）、 （ａ）式は細片抽出温度に関する制約式、
（４）、 （５）式は設備能力からくる燃料流量と炉温
の制約式を表わしている。
燃焼制御システムの目的は（１）〜（５）の数式と同時
に満足する解を求めることにあるが、そのために引例で
は様々の工夫をおこなっている。例えば、（１）〜（５
）式のｖ（ｉ）、θ００？ と炉温Ｔ（ｉ）を導びっ
ける炉内現象モデル式の誘導、そのモデル式をベースに
して（１）〜（５）式を線形化するアリゴリズムの導出
。
に満足する解を求めることにあるが、そのために引例で
は様々の工夫をおこなっている。例えば、（１）〜（５
）式のｖ（ｉ）、θ００？ と炉温Ｔ（ｉ）を導びっ
ける炉内現象モデル式の誘導、そのモデル式をベースに
して（１）〜（５）式を線形化するアリゴリズムの導出
。
線形化された（１）〜（５）式の解を求める線形計画法
アルゴリズムの採用等がめる。最終的（それらのアルゴ
リズムは適当なプログラム言語を用いて表現され、実用
に供されることになる。すなわち、ある数式（もしくは
数式群）が与えられたとき、それを計算機で実行可能な
形にするまでには、数式間の関係を記述するモデルを導
出するステップ、数式を計算機で解くためのアルゴリズ
ムを決定するステップ、そのモデルやアルゴリズムをプ
ログラムに変換するというステップ等を踏む必要がある
。従って新しい数式が追加になったり、修正になるとア
リゴリズムの修正、プログラムの修正が必要になる。ま
たそれ以前に作ったアルゴリズムやプログラムとの整合
性も充分考慮されねばならない。
アルゴリズムの採用等がめる。最終的（それらのアルゴ
リズムは適当なプログラム言語を用いて表現され、実用
に供されることになる。すなわち、ある数式（もしくは
数式群）が与えられたとき、それを計算機で実行可能な
形にするまでには、数式間の関係を記述するモデルを導
出するステップ、数式を計算機で解くためのアルゴリズ
ムを決定するステップ、そのモデルやアルゴリズムをプ
ログラムに変換するというステップ等を踏む必要がある
。従って新しい数式が追加になったり、修正になるとア
リゴリズムの修正、プログラムの修正が必要になる。ま
たそれ以前に作ったアルゴリズムやプログラムとの整合
性も充分考慮されねばならない。
これらのアルゴリズムの検討、プログラム修正方式の決
定、プログラムの作成といった作業はシステムエンジニ
アとかプログラマとよげれる人間がおこなうことになる
が、数式からプログラムに変換する各ステップのミスを
完全に防ぐことは非常ＫＪＩＩい。
定、プログラムの作成といった作業はシステムエンジニ
アとかプログラマとよげれる人間がおこなうことになる
が、数式からプログラムに変換する各ステップのミスを
完全に防ぐことは非常ＫＪＩＩい。
なお、（１）〜（５）式はめる典形的な制御目標を表わ
したものであって、常に（１）〜（５）式を使えば良い
というものではない。例えば鋼片の材質が異なれば（１
）〜（５）の代りに別の制御目標が用いられることがあ
るし、同じ材質のものでも抽出時間の長さや鋼片の充法
等によって（１）〜（５）式の中のめる式を使わなかっ
たり、別のものでおきかえるといったことが必要（なる
。炉の操業条件に応じてどの数式をどのように組合せて
使うか（数式使用条件）は、操業Ｐ％に負うところが大
きい。このため操業経験が深まるとともに数式の使い方
も変化することもしばしば発生する。従来のように数式
をプログラムとして書いている場合には、その変化に合
せてプログラムの変更が必要になり、やはり大きな問題
を生ずる。
したものであって、常に（１）〜（５）式を使えば良い
というものではない。例えば鋼片の材質が異なれば（１
）〜（５）の代りに別の制御目標が用いられることがあ
るし、同じ材質のものでも抽出時間の長さや鋼片の充法
等によって（１）〜（５）式の中のめる式を使わなかっ
たり、別のものでおきかえるといったことが必要（なる
。炉の操業条件に応じてどの数式をどのように組合せて
使うか（数式使用条件）は、操業Ｐ％に負うところが大
きい。このため操業経験が深まるとともに数式の使い方
も変化することもしばしば発生する。従来のように数式
をプログラムとして書いている場合には、その変化に合
せてプログラムの変更が必要になり、やはり大きな問題
を生ずる。
本発明の目的は、かかる従来方式の欠点に鑑みユーザが
数式を入力すると、その数式に対応するプログラムを自
動的に生成し、作成済みの実行プログラムと自動的にリ
ンクする機能をもった数式処理システムを提供すること
にある。
数式を入力すると、その数式に対応するプログラムを自
動的に生成し、作成済みの実行プログラムと自動的にリ
ンクする機能をもった数式処理システムを提供すること
にある。
本発明は数式もしくは数式と数式使用条件の組をプログ
ラムとは独立したデータとして入力して記憶する。記憶
した数式の構文解析をおこなってそれに基づいて演算プ
ログラムを生成し、この生成したプログラムを実行して
情報処理システムにその結果を返す。このような機能を
備えることにより情報処理に任意の数式処理を埋め込ん
だのと同等の効果を与える。
ラムとは独立したデータとして入力して記憶する。記憶
した数式の構文解析をおこなってそれに基づいて演算プ
ログラムを生成し、この生成したプログラムを実行して
情報処理システムにその結果を返す。このような機能を
備えることにより情報処理に任意の数式処理を埋め込ん
だのと同等の効果を与える。
本発明の具体的実施例を第１図、第３図〜第１０図を用
いて加熱炉制御の場合について説明する。第２図におい
て１００は制御部、２００は数式処理部、１０１は移動
スケジュール予測部、１０３は鋼片温度予測部、１０２
は予測計算部、１０４は操作量予測部を表わす。ここで
１０１゜１０３および１０４は、ｔａｘ図の２１．２３
および２４と同一の機能と構成をもつものでその詳細説
明は省略する。一方制御目標に直接関係している予測計
算１０２は第１図の２２と異なり、数式処理部２００と
のインタラクションを介して最適な鋼片昇温曲線の決定
する。
いて加熱炉制御の場合について説明する。第２図におい
て１００は制御部、２００は数式処理部、１０１は移動
スケジュール予測部、１０３は鋼片温度予測部、１０２
は予測計算部、１０４は操作量予測部を表わす。ここで
１０１゜１０３および１０４は、ｔａｘ図の２１．２３
および２４と同一の機能と構成をもつものでその詳細説
明は省略する。一方制御目標に直接関係している予測計
算１０２は第１図の２２と異なり、数式処理部２００と
のインタラクションを介して最適な鋼片昇温曲線の決定
する。
すなわち、予測計算部１０２は先ず数式処理部２００に
対し、炉操業データから、どのような制御目標を使うべ
きかを問い合せる。数式処理部２００では数式使用条件
を炉操業データを比較し数式を選択する。さらに数式処
理部２００ではこの選択された数式を実行可能なプログ
ラムに自動変換し、予測計算部１０２に返す。予測計算
部１０２はこのプログラムを受は取って実行し、その結
果に基づいて昇温曲線を決定する。
対し、炉操業データから、どのような制御目標を使うべ
きかを問い合せる。数式処理部２００では数式使用条件
を炉操業データを比較し数式を選択する。さらに数式処
理部２００ではこの選択された数式を実行可能なプログ
ラムに自動変換し、予測計算部１０２に返す。予測計算
部１０２はこのプログラムを受は取って実行し、その結
果に基づいて昇温曲線を決定する。
次に予測計算部１０２と数式処理部２００の詳細を第３
図を用いて説明する。実線はプログラムフローを、点紳
はデータフローを示す。第３図において１０２は第２図
の予測計算部に対応し、１０２−ａ〜１０２−Ｊはその
処理のステップを表わしている。また３００は後述する
インタープリタを表わし、２００は第２図の２００に対
応する数式処理部を表わす。また２０１はユーザとの対
話形式により数式を入力するマンマシン処理部、２０２
は数式選択部、２０３は選択された数式を実行可能なプ
ログラム（自動変換するプログラム生成部を表わす。
図を用いて説明する。実線はプログラムフローを、点紳
はデータフローを示す。第３図において１０２は第２図
の予測計算部に対応し、１０２−ａ〜１０２−Ｊはその
処理のステップを表わしている。また３００は後述する
インタープリタを表わし、２００は第２図の２００に対
応する数式処理部を表わす。また２０１はユーザとの対
話形式により数式を入力するマンマシン処理部、２０２
は数式選択部、２０３は選択された数式を実行可能なプ
ログラム（自動変換するプログラム生成部を表わす。
またＦ１〜Ｆ９けデータファイルを表わしている。ＦＩ
Ｆｉ数式選択に必要な炉内状況データ（鋼種、板厚、在
炉時間等）を数式処理部にわたすためのインターフェー
スファイル、Ｆ２〜Ｆ４は数式処理部２００から予測計
算部１０２へ渡す実行可能プログラムの入ったファイル
、Ｆ５．・Ｆ６はマンマシン部２０１を介して取り込ん
だ数式および数式使用条件を格納するファイル、Ｆ７〜
Ｆ９は予測計算部で用いるワークファイルである。
Ｆｉ数式選択に必要な炉内状況データ（鋼種、板厚、在
炉時間等）を数式処理部にわたすためのインターフェー
スファイル、Ｆ２〜Ｆ４は数式処理部２００から予測計
算部１０２へ渡す実行可能プログラムの入ったファイル
、Ｆ５．・Ｆ６はマンマシン部２０１を介して取り込ん
だ数式および数式使用条件を格納するファイル、Ｆ７〜
Ｆ９は予測計算部で用いるワークファイルである。
予測計算部１０２の１０２−ａでは昇温予定曲線の対象
となっている鋼片（対象材と呼ぶ）の鋼種、板厚、抽出
までの時間（在炉時間）、炉温分布といった炉内状況を
表わすデータをファイルＦ１に格納する。次に１０２−
ｂで数式処理部２００を起動する。数式処理部２００の
処理はあとで詳しく説明するが、炉内状況データに基づ
いて数式（評価関数や制約式）を選択する。さらにそれ
らの数式を計算可能なプログラム（ここでは中間コード
と呼ぶ）に変換して、Ｙの結果をインターフェースファ
イルＦ２．Ｆ３．Ｆ４に格納する。
となっている鋼片（対象材と呼ぶ）の鋼種、板厚、抽出
までの時間（在炉時間）、炉温分布といった炉内状況を
表わすデータをファイルＦ１に格納する。次に１０２−
ｂで数式処理部２００を起動する。数式処理部２００の
処理はあとで詳しく説明するが、炉内状況データに基づ
いて数式（評価関数や制約式）を選択する。さらにそれ
らの数式を計算可能なプログラム（ここでは中間コード
と呼ぶ）に変換して、Ｙの結果をインターフェースファ
イルＦ２．Ｆ３．Ｆ４に格納する。
予測計算部１０２の１０２−Ｃ〜１０２−Ｈは（１）〜
（５）式で代表される数式から線形計画プログラムで使
用するデータを生成するステップであり、１０２−Ｉが
線形計画プログラムを実行するステップ、１０２−ＪＦ
ｉ線形線形計画上る計算結果が求まる最適炉温分布（Ｔ
（１）、 ’ｒ（２）、・・・・・・、Ｔか））ｏパを
求め、その炉温分布で鋼片を加熱した場合の鋼片昇温曲
線を計算するステップでおる。この中で制御目標を表わ
す数式に直接関係するのは、１０２−Ｅ、１０２−Ｇお
よび１０２−Ｈであり、ここで実行すべきプログラムは
数式処理部２００から［ｆｔタインターフェースフフイ
ルＦ２．Ｆ３゜Ｆ４に中間コードとして格納されている
（詳細は後述）、１０２−Ｅ、１０２−ｃ＊、１０２−
Ｈではこれらの中間コードを実行するために、インター
プリタ−３００を起動するだけであり、特定の制御目標
に対応する数式の処理手続きは１０２−Ｅ、１０２−Ｇ
、１０２−Ｈには−切かかれていない。従って１０２−
ａ〜１０２−Ｊの冬季続きけ特定の制御目標に依存しな
い標準的プログラムである。
（５）式で代表される数式から線形計画プログラムで使
用するデータを生成するステップであり、１０２−Ｉが
線形計画プログラムを実行するステップ、１０２−ＪＦ
ｉ線形線形計画上る計算結果が求まる最適炉温分布（Ｔ
（１）、 ’ｒ（２）、・・・・・・、Ｔか））ｏパを
求め、その炉温分布で鋼片を加熱した場合の鋼片昇温曲
線を計算するステップでおる。この中で制御目標を表わ
す数式に直接関係するのは、１０２−Ｅ、１０２−Ｇお
よび１０２−Ｈであり、ここで実行すべきプログラムは
数式処理部２００から［ｆｔタインターフェースフフイ
ルＦ２．Ｆ３゜Ｆ４に中間コードとして格納されている
（詳細は後述）、１０２−Ｅ、１０２−ｃ＊、１０２−
Ｈではこれらの中間コードを実行するために、インター
プリタ−３００を起動するだけであり、特定の制御目標
に対応する数式の処理手続きは１０２−Ｅ、１０２−Ｇ
、１０２−Ｈには−切かかれていない。従って１０２−
ａ〜１０２−Ｊの冬季続きけ特定の制御目標に依存しな
い標準的プログラムである。
なお以下では数式処理部２００によって制御目標として
（１）、 （２）、 （４）、 （５）式が選択された
ものとして説明を進める。また、説明が必要以上に？！
雑になるのを避けるため（２）、 （４）、 （５）式
の上下限値は定数で、時間分割数ｎと時間長Δτｌの値
も決っているものとする。
（１）、 （２）、 （４）、 （５）式が選択された
ものとして説明を進める。また、説明が必要以上に？！
雑になるのを避けるため（２）、 （４）、 （５）式
の上下限値は定数で、時間分割数ｎと時間長Δτｌの値
も決っているものとする。
予測計算部１０２の処理の内容をもう少し詳しく説明す
る。予測計算部の１０２−Ｃでは適当な炉温分布（基準
炉温分布） （Ｔ（１）、 Ｔ（２）、・・・・・・。
る。予測計算部の１０２−Ｃでは適当な炉温分布（基準
炉温分布） （Ｔ（１）、 Ｔ（２）、・・・・・・。
Ｔ（ｎ））を決定し、１０２−Ｄの部分で炉内現象モデ
ルを用いて鋼片の内部温度分布と燃料流量をシミュレー
ションで決定する。炉内現象モデルとしては熱伝導方程
式、燃料流量式を用いれば良い。
ルを用いて鋼片の内部温度分布と燃料流量をシミュレー
ションで決定する。炉内現象モデルとしては熱伝導方程
式、燃料流量式を用いれば良い。
この炉内現象モデルで計算される諸量は炉内の状況を推
定するための基本的な物理景であり、どのような制御目
標を用いるかとは直接関係ない。
定するための基本的な物理景であり、どのような制御目
標を用いるかとは直接関係ない。
１０２−Ｄは、この炉内現象を計算するプログラムその
ものを表わしている。
ものを表わしている。
ここでは計算で求まる鋼片の内部温度分布を第４図に示
したようにθ（１）、・・・・・・θ（５）で表わし、
各時間ゾーンの燃料流量の計算値をｖ（１）、ｖ（２）
、・・・・・・、ｖ（ｎ）で表わす。それらの計算結果
はプロセスデータ格納テーブルＦ７の所定の格納場所に
格納される。
したようにθ（１）、・・・・・・θ（５）で表わし、
各時間ゾーンの燃料流量の計算値をｖ（１）、ｖ（２）
、・・・・・・、ｖ（ｎ）で表わす。それらの計算結果
はプロセスデータ格納テーブルＦ７の所定の格納場所に
格納される。
制御目標の中には炉内現象モデルで直接計算していない
ような物理！（例えば（２）式のθｏｖｔが用いられて
いることがある。その値を決定するＫは、それらの物理
量を炉内現象モデルで計算した基本的物理骨と結びつけ
る関係式が必要である。θｏａｔの場合には、次のよう
な関係式が（２）式の定義と一緒におこなわれているも
のとする（マンマシン処理部２０１０所で説明する）。
ような物理！（例えば（２）式のθｏｖｔが用いられて
いることがある。その値を決定するＫは、それらの物理
量を炉内現象モデルで計算した基本的物理骨と結びつけ
る関係式が必要である。θｏａｔの場合には、次のよう
な関係式が（２）式の定義と一緒におこなわれているも
のとする（マンマシン処理部２０１０所で説明する）。
この関係式に対応するプログラムは数式処理２００によ
って自動生成されインター７エースフアイルＦ２に格納
されている。一方このプログラムの実行に用いられるデ
ータは、１０２−Ｄの所で説明したようにプロセスデー
タ格納テーブル・Ｆ７に格納されている。ステップ１０
２−Ｅからインタープリタ３００を起動すると、インタ
ープリタ−はファイルＦ２から中間コードと取り出し実
行する。この計算によりθｏｒ１テの値が決まり、ワー
クテーブルＦ８の゛所定の格納場所に格納される。
って自動生成されインター７エースフアイルＦ２に格納
されている。一方このプログラムの実行に用いられるデ
ータは、１０２−Ｄの所で説明したようにプロセスデー
タ格納テーブル・Ｆ７に格納されている。ステップ１０
２−Ｅからインタープリタ３００を起動すると、インタ
ープリタ−はファイルＦ２から中間コードと取り出し実
行する。この計算によりθｏｒ１テの値が決まり、ワー
クテーブルＦ８の゛所定の格納場所に格納される。
以上の１０２−Ｃ〜１０２−Ｈの計算は異ったいくつか
の炉温分布に関して繰返しおこなわれる。
の炉温分布に関して繰返しおこなわれる。
すなわち、先ず線形中心となる炉温分布（Ｔｏ（１）。
・・−・−・、 ’ｒｅ （ｎ） １に関し１０２−Ｃ
〜１０２−Ｅの計算をおこなってθｏＴ１↑やｖ （ｉ
）を決定する（これを００υ？、Ｏ’１０ （１）と表
わすことにする）。次に第１時刻ゾーンの炉温だけを’
ｒｏ（Ｊ）＋ΔＴ（、（Ｔは定数）として、同様の計算
を繰返す（ｊ＝１〜ｎ）。
〜１０２−Ｅの計算をおこなってθｏＴ１↑やｖ （ｉ
）を決定する（これを００υ？、Ｏ’１０ （１）と表
わすことにする）。次に第１時刻ゾーンの炉温だけを’
ｒｏ（Ｊ）＋ΔＴ（、（Ｔは定数）として、同様の計算
を繰返す（ｊ＝１〜ｎ）。
その結果をここではθｏｎ〒１Ｊｌｖｌ（１）とする。
結果はいずれもワークテーブルＦ８に格納される。
これらの結果から各時間ゾーンの炉温変化が各物理量に
及ぼす影響の大きさを次式によって計算できる。
及ぼす影響の大きさを次式によって計算できる。
数式処理部２００は（１）〜（５）式の中から００υ〒
やｖ（ｉ）のような物理量に関し、それらの線形係数を
計算するプログラム（７）、 （８）式に対応）を自動
的に生成し、インターフェースファイルＦ３に格納シて
おく。従って第３図の１０２−Ｇからインタープリタ３
００を起動すると、そのプログラムをインタープリタ３
００が実行する。これにより線形係数αｊ＋ｒｌ＋ｊの
値が求まる。
やｖ（ｉ）のような物理量に関し、それらの線形係数を
計算するプログラム（７）、 （８）式に対応）を自動
的に生成し、インターフェースファイルＦ３に格納シて
おく。従って第３図の１０２−Ｇからインタープリタ３
００を起動すると、そのプログラムをインタープリタ３
００が実行する。これにより線形係数αｊ＋ｒｌ＋ｊの
値が求まる。
ところで、制御目標の中に現われる物理量θｏａ’ｙと
、ｖ（１）は以上より次のように表わされる。
、ｖ（１）は以上より次のように表わされる。
ｖ（ｉ）＝ｖｏ（ｏ−ｔ−ΣｒＩＪＸ＋
・”ｉ（１（ｉ＝１〜ｎ） ここでＸｉは第１時間ゾーンの炉温微少変化量を表わし
ており、本実施例では、このＸｉを独立変数とみなして
いる。（９）、６０式ｔｌ−（１）、 （２）、 （４
）、 （５）式に代入して整理すると次のようになる。
・”ｉ（１（ｉ＝１〜ｎ） ここでＸｉは第１時間ゾーンの炉温微少変化量を表わし
ており、本実施例では、このＸｉを独立変数とみなして
いる。（９）、６０式ｔｌ−（１）、 （２）、 （４
）、 （５）式に代入して整理すると次のようになる。
Ｊ＝Ｃ，＋ｃｘ→最小 ・・・・・・ａ
ｏｂ Ｌ＜？ −ｘ ≦ｂ ａ−（Ｊ２ ここで、 ｃｏ ＝ Σ Ｗｅ（すΔＴ ＋
−０３・・・（ｌ！１９ ０３〜ａ６１式の各要素の値が決まれば線形計画法の欅
準プログラムにより（１υ、（Ｉ２を満す解が求められ
る。（１１〜（（１式には各要素の計算式が与えられて
いるが、これらの計算式は制御目標としてどの数式が選
ばれるかによって異ったものになる。数式処理部２００
では、数式が選択されたあと、その数式からＣ，、ｃ、
ｂＬ、ｂＵ、Ｈの各要素を計算するプログラムを生成し
、ファイルＦ４に格納しておく。＠３図の１０２−Ｈか
らインタープリタを起動し、ファイル４Ｐのプログラム
を実行すると、線形計画プログラムで必要なデータが生
成される（詳細は数式処理２００で述べる）。
ｏｂ Ｌ＜？ −ｘ ≦ｂ ａ−（Ｊ２ ここで、 ｃｏ ＝ Σ Ｗｅ（すΔＴ ＋
−０３・・・（ｌ！１９ ０３〜ａ６１式の各要素の値が決まれば線形計画法の欅
準プログラムにより（１υ、（Ｉ２を満す解が求められ
る。（１１〜（（１式には各要素の計算式が与えられて
いるが、これらの計算式は制御目標としてどの数式が選
ばれるかによって異ったものになる。数式処理部２００
では、数式が選択されたあと、その数式からＣ，、ｃ、
ｂＬ、ｂＵ、Ｈの各要素を計算するプログラムを生成し
、ファイルＦ４に格納しておく。＠３図の１０２−Ｈか
らインタープリタを起動し、ファイル４Ｐのプログラム
を実行すると、線形計画プログラムで必要なデータが生
成される（詳細は数式処理２００で述べる）。
この計算に引き続き線形計画性計算プログラムを１０２
−Ｉで起動すること（より標価関数ａυ式を最小（する
炉温の微小変化の最適値（Ｘ、 、ＯＦ〒。
−Ｉで起動すること（より標価関数ａυ式を最小（する
炉温の微小変化の最適値（Ｘ、 、ＯＦ〒。
Ｘ□ＯＦ？、・・・・・・Ｉ Ｘｓ＠ＯＦ！）が求まる
。このとき最適炉温分布（Ｔ（１）、 Ｔ（２）、 ・
・・−、Ｔ（ｎ） ） ＯＦ？は（Ｔ、（１）＋ ＸＩ
、ＯＦ？、 Ｔ６（２）十Ｘｌ０Ｐ？、 ・＋１
＋、Ｔｏ（ｎ）＋ｘｍ、ａｐｔ）で与えられる。第３図
の１０２−Ｊではこの最適炉温で加熱した場合の鋼片温
度変化を熱伝導方程式により求め、それを最適な鋼片昇
温曲線とする。
。このとき最適炉温分布（Ｔ（１）、 Ｔ（２）、 ・
・・−、Ｔ（ｎ） ） ＯＦ？は（Ｔ、（１）＋ ＸＩ
、ＯＦ？、 Ｔ６（２）十Ｘｌ０Ｐ？、 ・＋１
＋、Ｔｏ（ｎ）＋ｘｍ、ａｐｔ）で与えられる。第３図
の１０２−Ｊではこの最適炉温で加熱した場合の鋼片温
度変化を熱伝導方程式により求め、それを最適な鋼片昇
温曲線とする。
次に数式処理部２００の機能と構成を説明する。
第３図のマンマシン処理部２０１はユーザが数式の使用
条件および数式を入力するのをサポートシ、。
条件および数式を入力するのをサポートシ、。
その入力データをファイルＦ５．Ｆ６に格納する。
第５図（ａ） ＫファイルＦ５中の数式使用条件、同Φ
）にファイルＦ６中の数式の定義式の例を示す。また第
５図（Ｃ）は数式中に表われる物理量を炉内現象モデル
で計算される基本的物理量に結びつける式でファイルＦ
５に格納されている。第５図（ａ）に示したように数式
使用条件はｒＩＦ−ＴＨＥＮ・・・」の形で入力する。
）にファイルＦ６中の数式の定義式の例を示す。また第
５図（Ｃ）は数式中に表われる物理量を炉内現象モデル
で計算される基本的物理量に結びつける式でファイルＦ
５に格納されている。第５図（ａ）に示したように数式
使用条件はｒＩＦ−ＴＨＥＮ・・・」の形で入力する。
「〜」はある数式を使用する上で成り立つべき条件を表
わしており（前提部と呼ぶ）、それが成立するときに用
いられる数式の種類と名前が「・・・」に示されている
「・・・」を結論部と呼ぶ）。例えばルール１は［鋼片
の種類（Ｋｉｎｄ）が低炭素鋼（ｌｏｗ Ｃａｒｂｏ
ｒ）であれば、制約式（Ｃｏｎｄｉｔｉｏｎ） とし
て抽出温度平均値制約ａＶＥＱを使え」ということを意
味している。第５図（ｂ）の第１式はａｖＥｑの定義式
を表わしており、この中に登場するｔｈｅｔａ Ｑｕｔ
の定義式が第５図（１））の下方に示されている。これ
らの数式使用条件や数式、物理量の定義式はファイルＦ
５．Ｆ６に文字列として記憶されている。従って新しい
数式の追加や修正は、これらの文字列ファイルの編集を
意味する。結局マンマシン処理部２０１はこれらの文字
列ファイルの橢隼を対話的におこなうのをサポートする
。
わしており（前提部と呼ぶ）、それが成立するときに用
いられる数式の種類と名前が「・・・」に示されている
「・・・」を結論部と呼ぶ）。例えばルール１は［鋼片
の種類（Ｋｉｎｄ）が低炭素鋼（ｌｏｗ Ｃａｒｂｏ
ｒ）であれば、制約式（Ｃｏｎｄｉｔｉｏｎ） とし
て抽出温度平均値制約ａＶＥＱを使え」ということを意
味している。第５図（ｂ）の第１式はａｖＥｑの定義式
を表わしており、この中に登場するｔｈｅｔａ Ｑｕｔ
の定義式が第５図（１））の下方に示されている。これ
らの数式使用条件や数式、物理量の定義式はファイルＦ
５．Ｆ６に文字列として記憶されている。従って新しい
数式の追加や修正は、これらの文字列ファイルの編集を
意味する。結局マンマシン処理部２０１はこれらの文字
列ファイルの橢隼を対話的におこなうのをサポートする
。
第３図の数式選択部２０２は、予測計算部からインター
フェースファイルＦ１を渡されたとき、このファイルの
中にあるデータ（鋼材のｍｍ、板厚、在炉時間等）とフ
ァイルＦ５中の数式使用条件を比較し、どの数式を使う
べきかを決定する。
フェースファイルＦ１を渡されたとき、このファイルの
中にあるデータ（鋼材のｍｍ、板厚、在炉時間等）とフ
ァイルＦ５中の数式使用条件を比較し、どの数式を使う
べきかを決定する。
例えば一種（［１ｎｄ）として低炭素鋼（ＬＯＷＣａｒ
ｂｏｎ）というデータがインター７エースフアイルの中
にあれば、ファイルＦ５中のルール１により式ａｖｌｉ
：ｑを使うべきことが決定される。同様にしてインター
７エースフアイルＦ１中の全てのデータをファイルＦ５
中の全ルールに適用し、昇温曲線の計算に必要な全ての
数式を選び出す。
ｂｏｎ）というデータがインター７エースフアイルの中
にあれば、ファイルＦ５中のルール１により式ａｖｌｉ
：ｑを使うべきことが決定される。同様にしてインター
７エースフアイルＦ１中の全てのデータをファイルＦ５
中の全ルールに適用し、昇温曲線の計算に必要な全ての
数式を選び出す。
この数式選択部２０２の動作はファイルＦ１中の数式使
用条件の内容に関係なく、一定の処理をおこなう。すな
わち■ファイルＦ５中の条件を１つとり出し、■ファイ
ル１の中のデータの中に、とり出した条件の前提部と一
致するものがあるかどうか調べ、■一致すれば、その結
論部を図示していないファイルＦ、１０に格納し、■不
一致でめれは別の条件に対し■〜■を行なう。この処理
は使用条件が全てチェックされるまで続く。この説明か
らもわかるように数式選択部は条件が成立するかどうか
の判定（すなわち推論）の処理フローの制御を担当して
いる。使用条件の内容は全てファイルＦ５の中にあユ、
数式選択部からみると処理の対象となるデータにしかす
ぎない。
用条件の内容に関係なく、一定の処理をおこなう。すな
わち■ファイルＦ５中の条件を１つとり出し、■ファイ
ル１の中のデータの中に、とり出した条件の前提部と一
致するものがあるかどうか調べ、■一致すれば、その結
論部を図示していないファイルＦ、１０に格納し、■不
一致でめれは別の条件に対し■〜■を行なう。この処理
は使用条件が全てチェックされるまで続く。この説明か
らもわかるように数式選択部は条件が成立するかどうか
の判定（すなわち推論）の処理フローの制御を担当して
いる。使用条件の内容は全てファイルＦ５の中にあユ、
数式選択部からみると処理の対象となるデータにしかす
ぎない。
従って数式の使用条件４しくけそれに関連した数式が変
更になっても数式選択部２０２の処理手続き（すなわち
プログラム）を変更する必要はない。
更になっても数式選択部２０２の処理手続き（すなわち
プログラム）を変更する必要はない。
次に中間コード生成部２０３の処理を第６図を用いて説
明する。中間コード生成部２０３の役割け、ファイルＦ
ＩＯ中の選択済の数式（文字列として表わされている）
を実行可能なプログラムに変換することにある。ここで
実行可能プログラムは、演算の種類を表わすコードと、
各演算に用いられるデータ（もしくけデータ格納番地）
の組からカる。第７図に代表的演算コードを示す。
明する。中間コード生成部２０３の役割け、ファイルＦ
ＩＯ中の選択済の数式（文字列として表わされている）
を実行可能なプログラムに変換することにある。ここで
実行可能プログラムは、演算の種類を表わすコードと、
各演算に用いられるデータ（もしくけデータ格納番地）
の組からカる。第７図に代表的演算コードを示す。
’！、６図のれ環フローに従い、以下では次の例を用い
て中間コード生成部の処理を説明する。
て中間コード生成部の処理を説明する。
選択された数式
％式％
数式中の物理量の定義式
ｔ ｈ ｅ ｔ ａ Ｏｕ ｔ＝ （ＴＨＩ）＋２ Ｆ
［２）＋ ２ Ｔ［３）＋２Ｔ）［４）＋−Ｔ庄５））
／８ ・・・（１ｇステップ１（２０３−ａ） 選択された数式群の中から１つずつ数式をとり出し、数
式を構成要素に分解し、各要素が演算。
［２）＋ ２ Ｔ［３）＋２Ｔ）［４）＋−Ｔ庄５））
／８ ・・・（１ｇステップ１（２０３−ａ） 選択された数式群の中から１つずつ数式をとり出し、数
式を構成要素に分解し、各要素が演算。
変数、定数のいずれに該当しているかを調べる。
この手法自体は構文解析として知られており、その手法
自体は本特許の本質とけ直接関係ないので省略する。弐
〇Ｆｊを２０３−ａで処理した結果を第８図（ａ）に示
す。ここで、数式を構成する要素を示しており骨はその
要素が変数であることを示す。
自体は本特許の本質とけ直接関係ないので省略する。弐
〇Ｆｊを２０３−ａで処理した結果を第８図（ａ）に示
す。ここで、数式を構成する要素を示しており骨はその
要素が変数であることを示す。
第８図（ａ）のような表現はリスト表現と呼ばれること
がある。
がある。
ステップ２（２０３−ｂ）
次に第８図（ａ）のような数式の中に炉内現職モデルで
直接計算しているもの以外の変数（以下では未定義車数
と呼ぶ）があるかどうか調べる。もし未定Ｉ！ＬＩ変数
が含まれていれば、それの定義式をファイルＦ６からと
抄出してくる。第８図（ａ）の場合には変数ｔｈｅｔａ
Ｑｕｔが含まれているので、ファイルＦ６から０９式を
とり出してきて、ステップ１と同様にこの定義式を演算
順序を考慮したリストに変換する（第８図（ｂ））。
直接計算しているもの以外の変数（以下では未定義車数
と呼ぶ）があるかどうか調べる。もし未定Ｉ！ＬＩ変数
が含まれていれば、それの定義式をファイルＦ６からと
抄出してくる。第８図（ａ）の場合には変数ｔｈｅｔａ
Ｑｕｔが含まれているので、ファイルＦ６から０９式を
とり出してきて、ステップ１と同様にこの定義式を演算
順序を考慮したリストに変換する（第８図（ｂ））。
ステップ３（２０３−ｃ）
第８図（ｂ）のリストを実行可能プログラムに変換する
。この変換はカッコ（・・・）のレベルが深い方から浅
い方に向って計算が実行されるようにプログラム化され
る。第８図（Ｃ）は０１式の実行可能プログラムを示す
。ｍ８図（ｃ）の各行の第１要素は演算コードを、ｙ４
２．第３の要素は引数、第４要素は演算結果の格納場所
を示している。また定数以外の要素は全てデータとアク
セスする記憶場所を示している。第８図（ｂ）中のＴＨ
Ｉ、ＴＨ２，・・・・・・。
。この変換はカッコ（・・・）のレベルが深い方から浅
い方に向って計算が実行されるようにプログラム化され
る。第８図（Ｃ）は０１式の実行可能プログラムを示す
。ｍ８図（ｃ）の各行の第１要素は演算コードを、ｙ４
２．第３の要素は引数、第４要素は演算結果の格納場所
を示している。また定数以外の要素は全てデータとアク
セスする記憶場所を示している。第８図（ｂ）中のＴＨ
Ｉ、ＴＨ２，・・・・・・。
ＴＨ５は、第９図に示した変数変換テーブル（基づいて
、−片内部渾度の計算値が実際に格納されるＰＤ（４３
，ＦＤ（５）、・・・・・・、ＦＤ（８）に変換されて
い乙。第８図（Ｃ）中の１Ｗｉｌは、ｔｈｅｔａｏｕｔ
ｃ７）計算値を格納する記憶場所を示すものでファイル
Ｆ８の空きエリアに確保される。第３図１０２−Ｅの所
で述べたようにｔｈｅｔａＱｕｔの値は（ｎ＋１）回計
算されるので、記憶場所の大きさｔ（ｎ＋１）ケとられ
る。’ＩＷＴαｌ＃におけるＷＴ４（Ｉｎはその先頭の
アドレス、チは複数個記憶場所が存在することを示す。
、−片内部渾度の計算値が実際に格納されるＰＤ（４３
，ＦＤ（５）、・・・・・・、ＦＤ（８）に変換されて
い乙。第８図（Ｃ）中の１Ｗｉｌは、ｔｈｅｔａｏｕｔ
ｃ７）計算値を格納する記憶場所を示すものでファイル
Ｆ８の空きエリアに確保される。第３図１０２−Ｅの所
で述べたようにｔｈｅｔａＱｕｔの値は（ｎ＋１）回計
算されるので、記憶場所の大きさｔ（ｎ＋１）ケとられ
る。’ＩＷＴαｌ＃におけるＷＴ４（Ｉｎはその先頭の
アドレス、チは複数個記憶場所が存在することを示す。
第８図（ｄ）Ｋ ｔｈｅｔａ Ｏｕｔすなわち００７丁
の格納場所のとりがたを示す。
の格納場所のとりがたを示す。
ステップ４（２０３−ｄ）
ステップ３で決められた各変数の格納場所に関するデー
タを用いれば線形係数の計算プログラムは容易に作り出
せる。例えば第１時間ゾーンの線形係数α、は番地ＷＴ
Ｑυに格納されたｔｈｅｔａＱｕｔから番地ＷＴ（ｌ［
）中のｔｈｅｔａＱｕｔを引いてΔＴで割ることによっ
て求まる。従ってそのプログラムは第８図（ｅ）のよう
に表わされる。このプログラムでＷＴα腸はα、の格納
場所を示しておりファイルＦ７上の空きエリアの先頭で
ある。
タを用いれば線形係数の計算プログラムは容易に作り出
せる。例えば第１時間ゾーンの線形係数α、は番地ＷＴ
Ｑυに格納されたｔｈｅｔａＱｕｔから番地ＷＴ（ｌ［
）中のｔｈｅｔａＱｕｔを引いてΔＴで割ることによっ
て求まる。従ってそのプログラムは第８図（ｅ）のよう
に表わされる。このプログラムでＷＴα腸はα、の格納
場所を示しておりファイルＦ７上の空きエリアの先頭で
ある。
ステップ５ （２０３−６）
ステップ５では制御目標を表わす数式の線形式に関する
。計算プログラムを導く。Ｕの場合、次のような処理を
おこなうプログラムを生成する。
。計算プログラムを導く。Ｕの場合、次のような処理を
おこなうプログラムを生成する。
■ 線形近似式の下限値は
［１１２０−線形中心におけるｔｈｅｔａＱｕｔ〕（＠
８図（ｄ）のＷＴに格納されている）■ 線形近似式の
上限値は （１１５０−線形中心におけるｔｈｅｔａＱｕｔ）■
第１線形係数＝α、（＠８図（ｅ）のＷＴ’ ）第２線
形係数＝α！（第８図（ｅ）のＷＴ’＋１）ｙ４０線形
係数＝α、（第８図（ｅ）のＷＴ’ ＋ｒｌ−１）この
ような処理は選択された全数式に関しておこなわれる。
８図（ｄ）のＷＴに格納されている）■ 線形近似式の
上限値は （１１５０−線形中心におけるｔｈｅｔａＱｕｔ）■
第１線形係数＝α、（＠８図（ｅ）のＷＴ’ ）第２線
形係数＝α！（第８図（ｅ）のＷＴ’＋１）ｙ４０線形
係数＝α、（第８図（ｅ）のＷＴ’ ＋ｒｌ−１）この
ような処理は選択された全数式に関しておこなわれる。
ｔｔｉ式が第２番目に処理される数式の鴫曾、上記の処
理を行なうプログラムは第８図（ｆ）のように表わされ
る。ここでＢ Ｌ（２）、 Ｂ Ｕ（２）、 ＨＨ（２
，Ｉ）（Ｉ＝１〜ｎ）は第に番目の線形式の下限’ｔ’
Ｌ＋上限値、＃形係数噴を格納する場所で、線形計画法
の標準プログラムの引数となるものである。従って第８
図（ｆ）のプログラムを実行すれば、その計算結果が線
形計画法の標準プログラムに渡されることになる。
理を行なうプログラムは第８図（ｆ）のように表わされ
る。ここでＢ Ｌ（２）、 Ｂ Ｕ（２）、 ＨＨ（２
，Ｉ）（Ｉ＝１〜ｎ）は第に番目の線形式の下限’ｔ’
Ｌ＋上限値、＃形係数噴を格納する場所で、線形計画法
の標準プログラムの引数となるものである。従って第８
図（ｆ）のプログラムを実行すれば、その計算結果が線
形計画法の標準プログラムに渡されることになる。
＠８図（Ｃ）、 （ｅ）、 （ｆ）のプログラムを中間
コードと呼び、それぞれファイルＦ２．Ｆ３．Ｆ４に格
納して予測計算部１０２に渡す。
コードと呼び、それぞれファイルＦ２．Ｆ３．Ｆ４に格
納して予測計算部１０２に渡す。
最後にインタープリタ−３００の機能について説明する
。インタープリタ−３００は中間コードの格納されたフ
ァイルを渡されると、このファイルの先頭から順にデー
タを取り出し、解釈実行する。第１０図にインタープリ
タ−の構成を示す。
。インタープリタ−３００は中間コードの格納されたフ
ァイルを渡されると、このファイルの先頭から順にデー
タを取り出し、解釈実行する。第１０図にインタープリ
タ−の構成を示す。
第１０図において、３１０はコード判定部、３２１゜３
２２、・・・・・・は演算実行部、３０１はテンポラル
データファイルを表わす。コード判定部は、中間コード
格納ファイル（Ｆ２．Ｆ３もしくはＦ４）から先頭のデ
ータを引き出し、演算コードの種別を判定する。次に演
算コードの種別で決まる引数の個数分だけ中間コード格
納ファイルからデータをとり出す。第８図（ｂ）かられ
かるように演算コードの後のデータは、その演算に使う
データの格納番地を示している場合と、定数を示してい
る場合とがある。コード判定部は、データが番地をさし
ている“場合にはファイル（Ｆ７．Ｆ８ ｏｒ Ｆ９）
の膣当番地から演算に使う値を取り出す。これらの値を
引数として、各演算コードに対応する演算奥行部の１つ
を呼び出す。各演算実行部は、演算コードの種類分だけ
用意されている。これらの演算実行部はコンパイル済の
プログラムであり、引数に対し指定された演算をおこな
い、その結果を返す。演算結果を格納する場所も中間コ
ードとして指定されている（第８図（Ｃ） 、 （６）
、 （ｆ） ）ので、コード判定部３１０は、結果を
指定の番地に格納することができる。
２２、・・・・・・は演算実行部、３０１はテンポラル
データファイルを表わす。コード判定部は、中間コード
格納ファイル（Ｆ２．Ｆ３もしくはＦ４）から先頭のデ
ータを引き出し、演算コードの種別を判定する。次に演
算コードの種別で決まる引数の個数分だけ中間コード格
納ファイルからデータをとり出す。第８図（ｂ）かられ
かるように演算コードの後のデータは、その演算に使う
データの格納番地を示している場合と、定数を示してい
る場合とがある。コード判定部は、データが番地をさし
ている“場合にはファイル（Ｆ７．Ｆ８ ｏｒ Ｆ９）
の膣当番地から演算に使う値を取り出す。これらの値を
引数として、各演算コードに対応する演算奥行部の１つ
を呼び出す。各演算実行部は、演算コードの種類分だけ
用意されている。これらの演算実行部はコンパイル済の
プログラムであり、引数に対し指定された演算をおこな
い、その結果を返す。演算結果を格納する場所も中間コ
ードとして指定されている（第８図（Ｃ） 、 （６）
、 （ｆ） ）ので、コード判定部３１０は、結果を
指定の番地に格納することができる。
インタープリタ−３００は、このようにして、中間コー
ド格納ファイルからデータをとり出しながら実行を続け
、最終的に線形計画プログラムで必要とするデータを決
定する。
ド格納ファイルからデータをとり出しながら実行を続け
、最終的に線形計画プログラムで必要とするデータを決
定する。
本実施例によれば、加熱炉煤焼制御に用いる制御目標を
、人手（よるプログラム修正なしに追加。
、人手（よるプログラム修正なしに追加。
修正することができるので、燃焼制御の機能の拡張、変
更が簡単に実現できる。
更が簡単に実現できる。
本実施例では、制御目標を表わす数式をグログラムに自
動変換するようにしており、その変換を数式選択時にお
こなうものとして説明したが、これはマンマシン処理部
で制御目標を入力する際（自動変換しておいても良い。
動変換するようにしており、その変換を数式選択時にお
こなうものとして説明したが、これはマンマシン処理部
で制御目標を入力する際（自動変換しておいても良い。
こうすれば数式選択毎に毎回変換する必要がなくな９、
予測計算部にすぐにプログラムを渡すことができる。こ
れにより処理速度を大幅に改善できる。
予測計算部にすぐにプログラムを渡すことができる。こ
れにより処理速度を大幅に改善できる。
また実施例では数式を中間コードに変換し、インタープ
リタで実行する方式を述べたが、本システムを実行する
計算機のマシン語に直接変換することもできる。この場
合には異った機種に移殖する際に変換プログラムを作り
直さなければならないというデメリットはあるが、処理
速度はやはり大幅に改善される。
リタで実行する方式を述べたが、本システムを実行する
計算機のマシン語に直接変換することもできる。この場
合には異った機種に移殖する際に変換プログラムを作り
直さなければならないというデメリットはあるが、処理
速度はやはり大幅に改善される。
また第２図の実施例では制御目標を表わす数式を数式処
理部でプログラムに変換する例を述べたが、あらかじめ
使用する数式が全てわかっていて、その使用条件だけが
変更になる可能性がある場合も考えられる。この場合に
は、数式を全てプログラムモジュールとして制御部に用
意しておき、どの数式を使用すべきかを判定したり、数
式の中で用いるデータを条件によって判定するための推
論部をｆｓｃａ図の数式処理部２００のかわりに用いれ
ば良い。第３図の数式処理と同様に推論部では制御目標
の使用条件を記憶するファイルが設けられてお９、ユー
ザとの対話によりその使用条件の追加が、変更ができる
。推論部は、第３図の数式処理と同様に制御部から送ら
れたデータに基づいてどの使用条件が成立するかを判断
するが、成立した使用条件に対応するプログラムを生成
するかわりに、制御系に内蔵したどのプログラムモジュ
ールを起動するかを決定し、その結果を制御系に知らせ
る。この方式では数式のプログラムを自動的に生成する
ことはできないが、従来のシステムに比べ制御目標の使
い方を容易（変更できるというメリットがある。
理部でプログラムに変換する例を述べたが、あらかじめ
使用する数式が全てわかっていて、その使用条件だけが
変更になる可能性がある場合も考えられる。この場合に
は、数式を全てプログラムモジュールとして制御部に用
意しておき、どの数式を使用すべきかを判定したり、数
式の中で用いるデータを条件によって判定するための推
論部をｆｓｃａ図の数式処理部２００のかわりに用いれ
ば良い。第３図の数式処理と同様に推論部では制御目標
の使用条件を記憶するファイルが設けられてお９、ユー
ザとの対話によりその使用条件の追加が、変更ができる
。推論部は、第３図の数式処理と同様に制御部から送ら
れたデータに基づいてどの使用条件が成立するかを判断
するが、成立した使用条件に対応するプログラムを生成
するかわりに、制御系に内蔵したどのプログラムモジュ
ールを起動するかを決定し、その結果を制御系に知らせ
る。この方式では数式のプログラムを自動的に生成する
ことはできないが、従来のシステムに比べ制御目標の使
い方を容易（変更できるというメリットがある。
本発明によれば、複数の制御目標を同時に満足する制御
において、直接人手によってプログラムの変更をおこな
うことなく制御目標の追加、変更ができる。このためシ
ステムの機能の拡張、変更が容易になる。
において、直接人手によってプログラムの変更をおこな
うことなく制御目標の追加、変更ができる。このためシ
ステムの機能の拡張、変更が容易になる。
第１図は本発明の実施例、第２図は従来の加熱炉燃焼制
御の例、第３図は予測計算部と数式処理部の詳細説明、
第４図は一片内部温度の説明図、第５図は制御目標を表
わす数式使用条件と数式の具体例、第６図は中間コード
生成部の処理フロー、第７図は演算コード、ＷＪ８図は
プログラム生成の過程を説明するための図、ｆｉｌｃ９
図は変数変換テーブル、第１０図はインタープリタ−の
構成。 １００・・・制御部、１０１・・・玖動スケジュール予
測部、１０１・・・予測計算部、１０３・・・鋼片温度
予測部、１０４・・・操作量予測部、２００・・・数式
処理部。
御の例、第３図は予測計算部と数式処理部の詳細説明、
第４図は一片内部温度の説明図、第５図は制御目標を表
わす数式使用条件と数式の具体例、第６図は中間コード
生成部の処理フロー、第７図は演算コード、ＷＪ８図は
プログラム生成の過程を説明するための図、ｆｉｌｃ９
図は変数変換テーブル、第１０図はインタープリタ−の
構成。 １００・・・制御部、１０１・・・玖動スケジュール予
測部、１０１・・・予測計算部、１０３・・・鋼片温度
予測部、１０４・・・操作量予測部、２００・・・数式
処理部。
Claims (1)
- 【特許請求の範囲】 １、数式モデル処理に基づいて制御をおこなう数式処理
制御システムにおいて、複数の数式と各数式の使用条件
を入力して記憶する記憶手段と該制御システムから該数
式使用を判断するためのデータを与えそのデータと該使
用条件が一致する数式群を選択する手段と該選択された
数式群を実行可能プログラムに自動的に変換する手段と
、該プログラムを実行した結果を該処理システムに帰還
する手段とを備えたことを特徴とする数式処理制御シス
テム。 ２、前記特許請求の範囲第１項記載において、該記憶手
段に数式とその数式の使用条件とを組として記憶し、該
数式の使用条件を満足する数式を選択実行することを特
徴とする数式処理制御システム。 ３、前記特許請求の範囲第１項記載において、該数式の
使用条件はその制御系の制御目標であることを特徴とす
る数式処理制御システム。 ４、前記特許請求の範囲第３項記載において、該帰還手
段からのデータに基づいて、 記憶済みの各制御目標の使用条件が成立するかどうかを
判断し、該条件が成立した数式をプログラムに変換し、
該プログラムにより制御をおこなうことを特徴とする数
式処理制御システム。 ５、前記特許請求の範囲第４項において、制御目標の使
用条件を１つずつ記憶ファイルから取り出して該条件と
該データの比較し、その条件が成立すると判断された場
合には、該条件の組となつている数式をファイルに登録
し、条件が成立しない場合には別の使用条件に関し同様
の処理をおこなう推論をおこなうことを特徴とする数式
処理制御システム。 ６、前記特許請求の範囲第３、第４、第５項において数
式をプログラムに変換するにあたり、数式を演算子、変
数、定数を要素とするリストに分解し、各変数に記憶場
所を割り当て、該演算子の優先順位に基づく処理手順を
決定し該演算子とその引数となる変数と定数の組を該処
理手順に配列してそれをプログラムとすることを特徴と
すの数式処理制御システム。 ７、前記特許請求の範囲第１、第２、第３、第４、第５
、第６項において、該数式処理により得られたプログラ
ムを解釈実行するインタープリタを設けたことを特徴と
する数式処理制御システム。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59230185A JPH0789283B2 (ja) | 1984-11-02 | 1984-11-02 | 数式処理制御システム |
| KR1019850008138A KR920003499B1 (ko) | 1984-11-02 | 1985-11-01 | 지식처리를 응용한 계산기 제어 시스템 |
| US06/793,773 US4740886A (en) | 1984-11-02 | 1985-11-01 | Computer control system utilizing knowledge processing |
| CN85108971A CN1011267B (zh) | 1984-11-02 | 1985-11-01 | 利用智能处理的计算机控制系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59230185A JPH0789283B2 (ja) | 1984-11-02 | 1984-11-02 | 数式処理制御システム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS61110203A true JPS61110203A (ja) | 1986-05-28 |
| JPH0789283B2 JPH0789283B2 (ja) | 1995-09-27 |
Family
ID=16903925
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59230185A Expired - Lifetime JPH0789283B2 (ja) | 1984-11-02 | 1984-11-02 | 数式処理制御システム |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4740886A (ja) |
| JP (1) | JPH0789283B2 (ja) |
| KR (1) | KR920003499B1 (ja) |
| CN (1) | CN1011267B (ja) |
Families Citing this family (47)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62270712A (ja) * | 1986-05-20 | 1987-11-25 | Nippon Kokan Kk <Nkk> | 高炉状況検出方法 |
| US4910691A (en) * | 1987-09-30 | 1990-03-20 | E.I. Du Pont De Nemours & Co. | Process control system with multiple module sequence options |
| US4907167A (en) * | 1987-09-30 | 1990-03-06 | E. I. Du Pont De Nemours And Company | Process control system with action logging |
| US5218557A (en) * | 1987-10-08 | 1993-06-08 | Digital Equipment Corporation | Expert system for assisting in the design of a complex system |
| US4989166A (en) * | 1988-07-12 | 1991-01-29 | Hitachi, Ltd. | Method for synthesizing analysis model and flow analysis system |
| US5260883A (en) * | 1988-07-21 | 1993-11-09 | Daniel Flow Products, Inc. | Interactive data entry system and method for the design, engineering and manufacture of meter tubes |
| US5121496A (en) * | 1988-07-25 | 1992-06-09 | Westinghouse Electric Corp. | Method for creating, maintaining and using an expert system by recursively modifying calibration file and merging with standard file |
| KR920011084B1 (ko) * | 1988-08-04 | 1992-12-26 | 미쓰비시전기 주식회사 | 엘리베이터 시험장치 |
| JPH0245831A (ja) * | 1988-08-08 | 1990-02-15 | Ricoh Co Ltd | 知識推諭処理装置 |
| JPH07100262B2 (ja) * | 1988-10-07 | 1995-11-01 | 三菱電機株式会社 | 放電加工終了判定方法及びその装置 |
| US5058043A (en) * | 1989-04-05 | 1991-10-15 | E. I. Du Pont De Nemours & Co. (Inc.) | Batch process control using expert systems |
| US4975865A (en) * | 1989-05-31 | 1990-12-04 | Mitech Corporation | Method and apparatus for real-time control |
| JP2779208B2 (ja) * | 1989-06-12 | 1998-07-23 | 津田駒工業株式会社 | 織機回転数の制御方法および装置 |
| US5195026A (en) * | 1989-09-12 | 1993-03-16 | Hitachi, Ltd. | Method and apparatus for computer controlled nonlinear optimization |
| JPH03127122A (ja) * | 1989-10-12 | 1991-05-30 | Fujitsu Ltd | データ処理装置における演算処理方式 |
| US5159685A (en) * | 1989-12-06 | 1992-10-27 | Racal Data Communications Inc. | Expert system for communications network |
| US5126933A (en) * | 1990-02-07 | 1992-06-30 | White Iii Charles A | Self-learning memory unit for process controller and self-updating function generator |
| US5140537A (en) * | 1990-05-21 | 1992-08-18 | Hewlett-Packard Company | Modeling a factory with human operators and validating the model |
| CZ293613B6 (cs) * | 1992-01-17 | 2004-06-16 | Westinghouse Electric Corporation | Způsob monitorování chodu zařízení pomocí CPU |
| US5428525A (en) * | 1992-07-01 | 1995-06-27 | Cappelaere; Patrice G. | Computer system and method for signal control prioritizing and scheduling |
| JPH0772767A (ja) * | 1993-06-15 | 1995-03-17 | Xerox Corp | 対話型ユーザ支援システム |
| US6948070B1 (en) * | 1995-02-13 | 2005-09-20 | Intertrust Technologies Corporation | Systems and methods for secure transaction management and electronic rights protection |
| US7095854B1 (en) * | 1995-02-13 | 2006-08-22 | Intertrust Technologies Corp. | Systems and methods for secure transaction management and electronic rights protection |
| US7133845B1 (en) * | 1995-02-13 | 2006-11-07 | Intertrust Technologies Corp. | System and methods for secure transaction management and electronic rights protection |
| WO1996027155A2 (en) * | 1995-02-13 | 1996-09-06 | Electronic Publishing Resources, Inc. | Systems and methods for secure transaction management and electronic rights protection |
| EP0886813B1 (en) * | 1996-03-12 | 2003-01-08 | TYCO Electronics Corporation | Electrical heating systems |
| US5930798A (en) * | 1996-08-15 | 1999-07-27 | Predicate Logic, Inc. | Universal data measurement, analysis and control system |
| US6112126A (en) * | 1997-02-21 | 2000-08-29 | Baker Hughes Incorporated | Adaptive object-oriented optimization software system |
| US6230062B1 (en) | 1999-01-08 | 2001-05-08 | Voyan Technology | Adaptation to unmeasured variables |
| US6278898B1 (en) * | 1999-04-07 | 2001-08-21 | Voyan Technology | Model error bounds for identification of stochastic models for control design |
| CA2481432A1 (en) | 2002-04-19 | 2003-10-30 | Ronald Cass | Processing mixed numeric and/or non-numeric data |
| US7483868B2 (en) * | 2002-04-19 | 2009-01-27 | Computer Associates Think, Inc. | Automatic neural-net model generation and maintenance |
| US7777743B2 (en) * | 2002-04-19 | 2010-08-17 | Computer Associates Think, Inc. | Viewing multi-dimensional data through hierarchical visualization |
| US7444310B2 (en) * | 2002-04-19 | 2008-10-28 | Computer Associates Think, Inc. | Automatic model maintenance through local nets |
| GB0209543D0 (en) * | 2002-04-26 | 2002-06-05 | Rolls Royce Plc | The automation and optimisation of the design of a component |
| US7313279B2 (en) * | 2003-07-08 | 2007-12-25 | Computer Associates Think, Inc. | Hierarchical determination of feature relevancy |
| EP1636738A2 (en) * | 2003-05-23 | 2006-03-22 | Computer Associates Think, Inc. | Adaptive learning enhancement to auotmated model maintenance |
| US8589869B2 (en) | 2006-09-07 | 2013-11-19 | Wolfram Alpha Llc | Methods and systems for determining a formula |
| US9213768B1 (en) | 2009-05-15 | 2015-12-15 | Wolfram Alpha Llc | Assumption mechanism for queries |
| US8601015B1 (en) | 2009-05-15 | 2013-12-03 | Wolfram Alpha Llc | Dynamic example generation for queries |
| US8484015B1 (en) | 2010-05-14 | 2013-07-09 | Wolfram Alpha Llc | Entity pages |
| US8812298B1 (en) | 2010-07-28 | 2014-08-19 | Wolfram Alpha Llc | Macro replacement of natural language input |
| US9058029B2 (en) * | 2011-03-31 | 2015-06-16 | Brad Radl | System and method for creating a graphical control programming environment |
| US9069814B2 (en) | 2011-07-27 | 2015-06-30 | Wolfram Alpha Llc | Method and system for using natural language to generate widgets |
| US9734252B2 (en) | 2011-09-08 | 2017-08-15 | Wolfram Alpha Llc | Method and system for analyzing data using a query answering system |
| US9851950B2 (en) | 2011-11-15 | 2017-12-26 | Wolfram Alpha Llc | Programming in a precise syntax using natural language |
| US9405424B2 (en) | 2012-08-29 | 2016-08-02 | Wolfram Alpha, Llc | Method and system for distributing and displaying graphical items |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5833709A (ja) * | 1981-08-24 | 1983-02-28 | Nippon Steel Corp | 数式制御モデルのオンライン変更方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| ZA801538B (en) * | 1979-03-19 | 1981-04-29 | United Kingdom Government | Method and apparatus for estimating slant visibility in fog |
| JPS58105373A (ja) * | 1980-01-12 | 1983-06-23 | メゾガツダサジ・エス・エレルミスゼリパリ・スゼルヴエゾ・ヴアララ−ト | 電子装置 |
| GB2079492B (en) * | 1980-06-25 | 1984-06-06 | Yokogawa Electric Works Ltd | Programmable controller |
| US4472770A (en) * | 1981-06-22 | 1984-09-18 | Li Chou H | Self-optimizing method and machine |
| US4459259A (en) * | 1982-06-29 | 1984-07-10 | The United States Of America As Represented By The United States Department Of Energy | Digital computer operation of a nuclear reactor |
| JPH0736123B2 (ja) * | 1983-05-09 | 1995-04-19 | 株式会社日立製作所 | 設備群制御方法 |
| DE3390539C2 (de) * | 1983-10-14 | 1990-03-08 | Ford Werke Ag | Verfahren zum Betreiben eines adaptiv gesteuerten technischen Prozesses |
-
1984
- 1984-11-02 JP JP59230185A patent/JPH0789283B2/ja not_active Expired - Lifetime
-
1985
- 1985-11-01 KR KR1019850008138A patent/KR920003499B1/ko not_active IP Right Cessation
- 1985-11-01 CN CN85108971A patent/CN1011267B/zh not_active Expired
- 1985-11-01 US US06/793,773 patent/US4740886A/en not_active Expired - Lifetime
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5833709A (ja) * | 1981-08-24 | 1983-02-28 | Nippon Steel Corp | 数式制御モデルのオンライン変更方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1011267B (zh) | 1991-01-16 |
| KR920003499B1 (ko) | 1992-05-01 |
| CN85108971A (zh) | 1986-11-05 |
| US4740886A (en) | 1988-04-26 |
| KR860004354A (ko) | 1986-06-20 |
| JPH0789283B2 (ja) | 1995-09-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPS61110203A (ja) | 数式処理制御システム | |
| Yalaoui et al. | Integrated production planning and preventive maintenance in deteriorating production systems | |
| US7444312B2 (en) | Method and arrangement for designing a technical system | |
| TW200413869A (en) | Method and apparatus for adaptively determining weight factors within the context of an objective function | |
| MXPA06000105A (es) | Metodo y aparato para evaluar una solucion propuesta para un problema de restriccion. | |
| CN110826695B (zh) | 数据处理方法、装置和计算机可读存储介质 | |
| Spirin et al. | Software for the raw-materials management system in blast-furnace smelting | |
| Lavrov et al. | Software for decision-making support in blast-furnace operation | |
| EP3920113A1 (en) | Optimization apparatus, optimization method, and optimization program | |
| Kim et al. | Strategy to coordinate actions through a plant parameter prediction model during startup operation of a nuclear power plant | |
| Jouse et al. | Safety control of nuclear power operations using self-programming neural networks | |
| Vescovi et al. | Modeling and simulation of a complex industrial process | |
| Penny | An analysis, both theoretical and by simulation, of a time-shared computer system | |
| Bazargan-Lari et al. | Intra-cell and inter-cell layout designs for cellular manufacturing | |
| Kobysh et al. | Situational control of hot blast stoves group based on decision tree | |
| Long et al. | Cycle 1 Development of Initial HCI Engineering Design Principles for Domestic Energy Planning and Control | |
| JP2579996B2 (ja) | 炉心状態予測装置 | |
| Deguchi et al. | Operation support system using dynamic simulation for a combined batch/continuous plant | |
| JPH01257300A (ja) | 計画作成支援方法およびシステム | |
| Gladkov et al. | Hybrid Approach to Solving the Problems of Operational Production Planning | |
| Kellett et al. | Concurrent METATEM as a coordination language | |
| Castillo-Durán et al. | A stability assessment of optimum Fuel Reload Patterns for a BWR | |
| Yang et al. | Computationally efficient nonlinear model predictive control | |
| Vannucci et al. | A simplified approach to the simulation of direct reduction of iron ore | |
| JP2907621B2 (ja) | 操業知識自動獲得方法および装置 |