JPH071222B2

JPH071222B2 - プロセス状態解析方式

Info

Publication number: JPH071222B2
Application number: JP3036113A
Authority: JP
Inventors: 聡大石; 幸一川口; 和夫小山; 賢二中村
Original assignee: Hitachi Ltd; Tokyo Gas Co Ltd
Current assignee: Hitachi Ltd; Tokyo Gas Co Ltd
Priority date: 1991-03-01
Filing date: 1991-03-01
Publication date: 1995-01-11
Anticipated expiration: 2010-01-11
Also published as: JPH04274725A

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、大量の時系列データを
利用して、プロセス状態解析を行なう解析手法に係り、
特に波形状態の一種であるハンチングの有無とその発散
・収束の情報を得るに好適なプロセス状態解析方式に関
する。

【０００２】

【従来の技術】従来のプロセス状態解析方式において
は、図１７に示すように、プロセスデータを入出力する
プロセスデータ入出力装置１と、プロセスデータを記憶
しプラントを制御及び監視する情報を演算する電子計算
機２０と、この情報とプロセスデータとを表示する表示
手段（ＣＲＴ）２１と、表示手段２１を経由して演算に
用いる情報を入力する入力手段（マンマシン機構）３と
よりなり、数値解析手法であるＦＦＴ（Ｆast Ｆouri−
er Ｔransform）を用ることにより、時系列データをグ
ラフに表現した時の波形の主成分周波数と、その周波数
における振幅とを解析結果として得ていた。

【０００３】知識処理システムやインテリジェントアラ
ームシステムを実現しようとする場合、次のような問題
があげられる。

【０００４】従来から人間（オペレータ）の直感や経験
によって判断していた状態把握、例えば、オペレータが
記録計のトレンドチャート紙の印字状況や、計算機シス
テムのＣＲＴ等に表示されるトレンドグラフの波形を参
照して、現在や過去の状況がどのような状態であるかを
判断する。例えば「このトレンドは上昇している」、
「このトレンドは急降下している」、「波形がハンチン
グしている」、「波形が発散状態でハンチングしてい
る」等と言うように技術用語や専門用語で表現される状
態変数に置き換えて判断結果を出すのである。

【０００５】このように、知識処理システムやインテリ
ジェントアラームシステムを実現しようとする場合、計
算機システム内に蓄えられている時系列データ（数値デ
ータ）を上記のような状態変数に置換する処理、特にハ
ンチングを理解する必要については、最適な処理方法が
なかった。

【０００６】

【発明が解決しようとする課題】従来のプロセス状態解
析方式にあっては、知識処理システムやインテリジェン
トアラームシステムを実現しようとする場合、オペレー
タがトレンドグラフを経験的によって判断し、状態把握
を行なっており、プロセス状態のハンチングを理解する
最適な処理方法がないという問題点があった。

【０００７】本発明の目的は、特に「ハンチング」と呼
ばれる周期的に変動または不規則に動揺する現象にある
時系列な数値データを、数値解析処理によってその特徴
を抽出することのできるプロセス状態解析方式を提供す
ることにある。

【０００８】

【課題を解決するための手段】前記の目的を達成するた
め、本発明に係るプロセス状態解析方式は、プラントの
プロセスデータを入出するプロセスデータ入出力装置
と、プロセスデータを記憶しかつ演算してプラントを制
御及び監視する電子計算機と、この演算結果とプロセス
データとを表示する表示手段と、演算に用いる情報を入
力する入力手段とよりなるプロセス状態解析方式におい
て、電子計算機は、プロセスデータの時系列データを重
複時間帯を持たせて複数分割し、分割したそれぞれの時
系列データをＦＦＴ解析し、その解析結果によりハンチ
ングの状態を判断する解析手段を備えている構成とす
る。

【０００９】そして解析手段は、分割数Ｎに１を加算し
たＮ＋１回のＦＦＴ解析を行ない、第１回のＦＦＴ解析
をすべての時系列データに行なってハンチングの有無を
判断する構成でもよい。

【００１０】また、解析手段は、それぞれのＦＦＴ解析
結果を一次回帰処理し、ハンチングの状態を発散傾向、
収束傾向又は安定状態の状態変数に置き換える構成であ
る。

【００１１】さらに表示手段は、それぞれの状態変数
を、表示又は印字で出力する手段であるものとする。

【００１２】そして、解析手段は、ＦＦＴ解析を周期的
に実行し、それぞれのＦＦＴ解析結果を一次回帰処理
し、将来のハンチングの状態を予測する構成でもよい。

【００１３】また表示手段は、マンマシン機構とＣＲＴ
よりなり、マンマシン機構を経由して状態変数がＣＲＴ
に表示される構成でもよい。

【００１４】

【作用】本発明のプロセス状態解析方式によれば、電子
計算機に設けたヒストリカルデータベースは、プロセス
データをはじめとする様々なデータを時系列的に時々刻
々と収集しデータベースに蓄えておき、データベースア
クセス機構は、そのデータベースへのアクセス管理を行
なう。

【００１５】マンマシン機構（入力手段）や他プログラ
ムから、プロセス解析機構を起動する情報（データ名
称、解析データ収集開始時刻、解析データ個数、ハンチ
ング条件各種規定値等）を与えることによって、指定さ
れた条件のもとで解析が行なわれる。

【００１６】また、解析結果は状態変数としてハンチン
グ「有／無」有りの場合は、「発散／収束／安定」とい
う文字が日本語表示され、状態変数を管理コードでも出
力される。

【００１７】プロセス状態解析機構（解析手段）は、
ヒストリカルデータ収集処理、ＦＴＴ実行処理、時
系列データ分割処理、一次回帰処理、状態変数変換
処理を行なう。

【００１８】ヒストリカルデータ収集処理起動情報として与えられた解析データ収集開始時刻と、
データ数よりターゲットになるヒストリカルデータとを
データベースアクセス機構を通して収集する。解析デー
タ収集開始時刻は、絶対時刻と相対時刻のどちらでも指
定可能である。

【００１９】ＦＴＴ実行処理本処理は、によって収集した時系列データを周波数系
列データに変換する。

【００２０】時系列データ分割処理本処理は、ハンチングの詳細情報を得るための前処理と
して、によって収集した時系列データを分割する。

【００２１】分割方法はランニングスペクトルとなるよ
うに行ない、ＦＴＴ実行後の最大振幅値に時間的要素が
付加される。

【００２２】一次回帰処理分割した時系列データを用い、それぞれＦＦＴを実施
し、それぞれの結果を求める。この時分割数分の最大振
幅値を時系列的に並べ、一次回帰処理を行なう。

【００２３】状態変数変換処理で得られる近似直線の傾きにより、オペレータが用い
ていた状態変数に置き換える処理を行なう。

【００２４】

【実施例】本発明を用いたそれぞれの実施例を図面を参
照しながら説明する。

【００２５】実施例１；インテリジェントアラーム本実施例は、プロセス制御システムにおいてプロセスア
ラームを従来の計器レベルのＨigh Ｅrror，Ｌow Ｅrro
r，偏差Ｅrr−or及び変化率Ｅrrorに加え、プロセス監
視処理によるプロセス傾向アラーム及びプロセス予測ア
ラームが該当する。

【００２６】プロセス傾向アラームは、本発明であるハ
ンチング状態解析処理やトレンド状態解析処理、相関状
態解析処理、ステップ状態解析処理及びパルス状態解析
処理等を用い、プロセス監視処理のパトロールによっ
て、予め指定されていたプロセス変数（Ｔag）に対し
て、「ハンチングが発生しています」、「トレンド状態
が急上昇で推移しています」、「トレンド状態が緩やか
な下降で推移しています」等、瞬時値レベルのアラーム
警報から過去データを加味したアラーム警報の出力が可
能になる。

【００２７】本実施例は、図１に示すようにプロセスデ
ータ入出力装置１と、データ収集機構２と、マンマシン
機構３と、プロセス制御処理機構４と、プラントデータ
ベース５と、ヒストリカルデータベース６と、データベ
ースアクセス機構７と、パトロール監視機構８と、プロ
セス状態解析機構９とより構成される。

【００２８】対象プラントにおける各種センサは、プロ
セス入出力装置１を介してデータ収集機構２によってデ
ータ収集され、プラントデータベース５に書き込まれ
る。プラントデータベース５のデータは、周期的にヒス
トリカルデータベース６に移され、ヒストリカルデータ
として長期間保存される。プラントデータベース５の内
容やヒストリカルデータベース６の内容は、マンマシン
機構３を経由して、容易に参照することができる。プロ
セス制御処理機構４は、収集されたプロセスデータに演
算処理を行ないデータ加工を行なう。また、データによ
ってはプロセス入出力装置１を経由して、対象プラント
に出力されるものもある。プロセス制御処理機構４で
は、収集したプロセスデータの合理性チェック等も実施
しており、この時、プロセスデータが測定レンジ範囲内
にあるか、プロセスデータがアラーム警報出力設定値を
越えているかをチェックするため、プロセスデータと前
回収集したプロセスデータとの偏差や変化率を算出し、
それぞれがアラーム警報出力設定値を越えているか等、
様々なチェックの実施結果によってアラーム警報をマン
マシン機構３に出力しＣＲＴ２１に表示する。

【００２９】プロセス状態解析機構（解析手段）には前
述の各種解析処理が格納されている。パトロール監視
機構８は、予め登録されているプロセス変数（Ｔag）を
管理し、周期的にプロセス状態解析機構９に実行要求を
出すスケジューリングを行ない、各種状態解析処理を実
行させ、解析結果を得て、その結果をアラーム警報とし
てマンマシン機構３の１つであるグラフィックＣＲＴ２
１に表示出力する。

【００３０】これらの機構が動作することにより、従来
のアラームにはない情報をオペレータに対して、警報と
いう形で情報提供できる。

【００３１】次にプロセス状態解析機構（解析手段）に
格納されているそれぞれの解析処理について図２〜図８
を参照しながら説明する。

【００３２】プロセスの「ハンチング」の状態を判断す
るには、与えられた時系列データが周期的要素を持って
変動しているか、さらにはその変動の状態が発散傾向、
収束傾向又は安定傾向を持っているかについて解析する
必要がある。

【００３３】周期的要素の変動については、与えられた
時系列データをＦＦＴ処理し数値解析処理により判断す
ることが可能である。

【００３４】図２にｆ（ｔ）で示されるＮ個の時系列デ
ータｆ（０），ｆ（１），ｆ（２）……ｆ（Ｎ−１）に
対するフーリエ変換は、（１）式で表わされる。

【００３５】

【数１】

【００３６】この値をｃｏｏｌｅｙ−Ｔｕｋｅｙのアル
ゴリズムによって求めることによって、時系列データｆ
（ｔ）は周波数系列データＡ（ｆ）に変換される。この
処理結果に基づいてＡ（₀）を直流成分、Ａ（ｆ）の最
大値（ｍａｘ（Ａ（ｆ）))を最大振幅値、ｍａｘ（Ａ
（ｆ））である最大振幅値となるｆを卓越周波数とす
る。直流成分と最大振幅値との比率ｍａｘ（Ａ（ｆ))／
Ａ₀が規定値より大きい場合を「ハンチング有」状態と
する。

【００３７】

【数２】

【００３８】Ｃ＝規定値云い換えると、与えられた時系列データは様々な周波数
形の合成波と考えることができ、この合成波のベースと
なる直流成分値と合成波のうち、最も振幅の大きい波形
を比較し、この比較結果により「ハンチング有」の状態
か否かを決定する。

【００３９】次に「ハンチング」の変動状態の解析は、
与えられた時系列データをＮ分割し、この分割後のＮ個
の時系列データをそれぞれＦＦＴ処理し、分割ＦＦＴを
行なう。ここでは、分割数Ｎは２以上と規定し、分割
される時系列データは、図３〜図５に示すように分割す
るものとする。分割手段は同図の例に限らないが、基本
的には時系列データが他の時系列データと、重複する部
分（重複時間帯）を有することが重要である。分割数Ｎ
の決定は任意であるが、分割数Ｎが増大すると解析結果
の精度が向上する。しかし、精度の向上は、解析結果の
出力までの計算処理時間に比例するため、計算処理時間
を確保しかつある程度の精度を得るようにする必要があ
り、本実施例では分割数Ｎ＝３で実施した。

【００４０】時系列データをＮ分割した後、それぞれの
時系列データに対してＦＦＴ処理を実行する。得られた
計算結果を｛Ａ₁（₀₁）₁、ｍａｘ（Ａ₁(ｆ₁))、ｆ₁｝、
{Ａ₂（₀₂）、ｍａｘ（Ａ₂（ｆ₂))、ｆ₂｝、
｛Ａ₃（₀₃）、ｍａｘ（Ａ₃（ｆ₃))、ｆ₃｝…とする。

【００４１】Ａ（₀ｘ）はＸ番目の時系列データのＦＦ
Ｔ処理結果の直流成分であり、ｍａｘ（Ａ（ｆｘ))は、
Ｘ番目の時系列データのＦＦＴ処理結果の最大振幅値、
ｆｘは、Ｘ番目の時系列データのＦＦＴ処理結果のｍａ
ｘ（Ａ（ｆ₀))、ｆ₀を得る卓越周波数（最大周波数）で
ある。

【００４２】このＮ個の解析結果は、これらの並びもま
た時系列に変化するデータと考えることができる。

【００４３】ただし、この｛Ａ（₀ｘ）、ｍａｘ(Ａ(ｆ
ｘ))、ｆｘ｝の決定方法は、分割前に行った処理結果
を｛Ａ（₀₀）、ｍａｘ(Ａ(ｆ₀))、ｆ₀｝としたとき、
分割後の個々の処理結果は、｛Ａ₁（₀₀）、ｍａｘ(Ａ
₁(ｆ₀))、ｆ₀｝、｛Ａ₂(₀₀）、ｍａｘ(Ａ₂(ｆ₀))、
ｆ₀｝、｛Ａ₃(₀₀）、ｍａｘ（Ａ₃(ｆ₀))、ｆ₀｝となる
ようなデータを求めるものである。

【００４４】第２図において、の結果は次のように
なる。

【００４５】分割前に行なった結果を｛Ａ₀，Ａｍ，
ｆｍ｝としたとき、分割後の個々の処理結果は、
｛(Ａ₁(ｆ₀))、ａ₁，ｆｍ｝、｛(Ａ₂(ｆ₀))、ａ₂，ｆ
ｍ｝、｛(Ａ₃(ｆ₀))、ａ₃，ｆｍ｝となる。

【００４６】このことは、で求めた最大振幅値を有す
る卓越周波数に注目し、分割後のＮ個の時系列データに
おける卓越周波数に対する最大振幅値の推移を見ていく
ものである。

【００４７】この最大振幅値を時系列的に並べ一次回帰
処理を行なうことにより、横軸ｔに対する最大振幅値の
変化を直線Ｌで近似し、「変動している場合、その状態
についての詳細情報」を求めることが可能となる。

【００４８】一般的に、近似した直線Ｌの傾きαが、０
近傍、０近傍を除く正の値、又は０近傍を除く負の値に
区分することができ、０近傍……安定状態０近傍を除く正の値……発散傾向０近傍を除く負の値……収束傾向にそれぞれ対応することができる。

【００４９】０近傍などの具体的かつ正確な数値につい
ては、対象プロセス及びデータ種別などに応じて、設計
知識などを加味して決定する必要がある。

【００５０】本処理は、時系列データの分割数Ｎにさら
に１回を加えたＮ＋１回のＦＦＴ処理と、１回の１次回
帰処理とから成り立っている。すなわち、第１回目の分
割前のすべての時系列データに対するＦＦＴ処理で「ハ
ンチングの有無」を判断し、残りのＮ回のＦＦＴ処理
と、１回の１次回帰処理によって「ハンチングの詳細情
報（発散、収束、安定）」を判断する。

【００５１】発散傾向とは図６に示すように、時間的経
過に伴って波形の振幅が大きくなっていく状態をさす。

【００５２】収束傾向とは図７に示すように、時間的経
過に伴って波形の振幅が小さくなっていく状態をさす。

【００５３】安定状態とは図８に示すように、時間的経
過に伴って波形の振幅が変動せず、安定した波形が出力
される状態をさす。（プロセスにおいて、この状態を安
定であるとは云わないが、ここでは、安定という用語を
用いた。）このように、波形の形状を状態変数である
「発散」、「収束」、「安定」という用語に置き換える
ことにより、オペレータが通常使用している用語とし
て、メッセージ、オペレーションガイド又は電報等で出
力する計算機システムが実現可能となる。

【００５４】図９は、前記の解析処理のフローを示すフ
ローチャートである。

【００５５】実施例２；知識処理システム本実施例は、プロセス制御システムに知識処理機構を付
加したシステムにおいて、知識処理機構内で実行される
推論内で、本発明であるプロセス状態解析機構９によっ
て得られる結果を推論データとして使用できるようにし
たものである、推論にてオペレータが使用している言語
（状態変数）をそのまま利用することが可能になる。

【００５６】実際のルール記述においては、と言う表現が可能になるため、オペレータが通常考えて
いる表現をそのまま用いることができ、知識整理、ルー
ル入力及び保守が、従来の表現に比べて効率的になる。

【００５７】本実施例は、図１０に示すように、プロセ
スデータ入出力装置１と、データ収集機構２と、マンマ
シン機構３と、プロセス制御処理機構４と、プラントデ
ータベース５と、ヒストリカルデータベース６と、デー
タベース機構７と、推論機構１０と、知識データベース
１１と、プロセス状態解析機構９とより構成される。

【００５８】対象プラントにおける各種センサは、プロ
セス入出力装置１を経由してデータ収集機構２によって
データ収集され、プラントデータベース５に書き込まれ
る。また、これ以外には、他計算機システムから図示し
ないデータ通信装置を経由して収集したデータを、プラ
ントデータベース５に書き込む場合も同様である。プラ
ントデータベース５のデータは周期的にヒストリカルデ
ータベース６に移され、ヒストリカルデータとして長期
保存される。プラントデータベース５の内容やヒストリ
カルデータベース６の内容は、マンマシン機構３を経由
して、容易に参照することができる。プロセス制御処理
機構４は、収集されたプロセスデータに演算処理を行な
いデータ加工を行なう。またデータによってはプロセス
入出力装置１を経由して、対象プラントに出力されるも
のもある。知識データベース１１にはＩＦ……ＴＨＥＮ
……形式で記述されたプロダクションルールが登録され
ており、推論機構１０によって知識データベース１１に
登録されたプロダクションルールが実行される。ルール
が実行された場合、ルール中のＴＨＥＮ部にプロセス状
態解析機構９への実行要求マクロ名称が記述されている
時、ルールが実行されることによってプロセス状態解析
機構９へ実行要求がかかる。この要求によって、プロセ
ス状態解析機構９が各種状態解析結果を行なう。解析処
理結果は、知識データベース１１の事実形知識を取扱う
フレームに反映される。反映した状態において、前述の
状態変数で表現したルールが実行されることになる。

【００５９】マンマシン機構３では、推論で実行された
状態解析処理結果を参照することができる。

【００６０】実施例３；プロセス予測図１に示す構成を有するシステムにおいて、図９に示す
処理を拡張することにより、予測処理が可能となり、オ
ペレータに将来のプロセスデータの推移（予測結果）を
知らせることができる。

【００６１】すなわち、図９における処理結果であるα
を図１に示すヒストリカルデータベース６に逐次蓄える
ことにより、ハンチング傾向状態の時系列データが生成
される。この時系列データをさらに一次回帰処理し、図
１１に示す近似直線Ｍを求めることにより、長期的なハ
ンチング傾向遷移をとらえることができる。すなわち、
図１１に示すｔ₁を現在時刻とすると、一次回帰処理結
果の近似直線Ｍを延長することにより、将来の時刻にお
けるαを予測することが可能となる。

【００６２】実施例４；ＣＲＴ表示画面前記処理結果のＣＲＴ表示画面の画面例を、図１２〜図
１６を参照しながら説明する。

【００６３】図１２の画面例１は、プロセス状態解析結
果をオペレータに知らせるアラームメッセージ出力を示
し、一般画面出力エリアの上端にアラームメッセージ出
力エリアが設けられている。

【００６４】図１３の画面例３は、プロセス状態解析実
行画面を示す。プロセス状態解析実行をＣＲＴから要求
する画面であって、解析する時間と種別とを入力し、実
行後は画面例３に展開する。画面例３は、プロセス状態
解析結果表示画面を示す。プロセス状態解析の実行結果
として発散、収束、安定を表示し、さらに用いた時系列
データをグラフとして表示する。画面例４はマルチウイ
ンドによる画面２，３を組合せた表示を示し、プロセス
状態解析で使用される内部データを詳細情報として表示
する。

【００６５】図１４の画面例５は、プロセス状態解析機
構を用いたプロセス予測アラーム（インテリジェントア
ラーム）画面を示す。プロセス状態解析を常時起動し、
プロセス状態が上昇中であることを表示するが、ハンチ
ングがある場合、発散、収束、安定を表示する。

【００６６】図１５の画面例６は、プロセス状態解析機
構を用いたｎ分後警報（インテリジェントアラーム）出
力Ｔａｇ一覧画面を示す。プロセス状態解析結果より、
ｎ分経過後の状態を表示するが、トレンド状態が上昇傾
向にあって、例えば５分後に上限警報値を越えてしまう
プロセスデータ名を表示している。

【００６７】図１６の画面例７は、プロセス状態解析機
構を用いた相関チェック結果（インテリジェントアラー
ム）出力画面を示す。画面例６と同様であるが、プロセ
ス状態解析に相関処理を適用した場合、プロセスデータ
名を表示している。

【００６８】

【発明の効果】本発明のプロセス状態解析方式によれ
ば、プロセス状態把握を計算機による自動処理化が図
れ、オペレータが記憶している表現と同じ表現で表示が
可能であるため理解されやすい。

【００６９】そして、プロセスアラームが検出できるた
め、プロセス状態把握する時、将来の挙動変化の予測が
できる。

【００７０】また知識処理に状態把握機能を持たせるこ
とができるため、長期的な挙動状況が把握可能となり、
ハンチングのゆっくりとした挙動もとらえることができ
る。

【図面の簡単な説明】

【図１】本発明の実施例１を示す構成図である。

【図２】本発明のハンチング状態解析処理を説明する図
である。

【図３】本発明のハンチング状態解析処理を説明する図
である。

【図４】本発明のハンチング状態解析処理を説明する図
である。

【図５】本発明のハンチング状態解析処理を説明する図
である。

【図６】本発明のハンチング状態解析処理を説明する図
である。

【図７】本発明のハンチング状態解析処理を説明する図
である。

【図８】本発明のハンチング状態解析処理を説明する図
である。

【図９】第２図〜第９図の解析処理を示すフローチャー
トである。

【図１０】本発明の実施例２を示す構成図である。

【図１１】本発明の実施例３を示す構成図である。

【図１２】本発明の実施例４を示す画面例の図である。

【図１３】本発明の実施例４を示す画面例の図である。

【図１４】本発明の実施例４を示す画面例の図である。

【図１５】本発明の実施例４を示す画面例の図である。

【図１６】本発明の実施例４を示す画面例の図である。

【図１７】従来の技術を示す図である。

【符号の説明】

１プロセス入出力装置３マンマシン機構（入力手段）９プロセス状態解析機構（解析手段）２０電子計算機２１ＣＲＴ（表示手段）

───────────────────────────────────────────────────── フロントページの続き (72)発明者小山和夫神奈川県横浜市鶴見区末広町１−７−７東京瓦斯株式会社生産技術センター内 (72)発明者中村賢二神奈川県横浜市鶴見区末広町１−７−７東京瓦斯株式会社生産技術センター内

Claims

【特許請求の範囲】

【請求項１】プラントのプロセスデータを入出するプ
ロセスデータ入出力装置と、前記プロセスデータを記憶
しかつ演算して前記プラントを制御及び監視する電子計
算機と、この演算結果と前記プロセスデータとを表示す
る表示手段と、前記演算に用いる情報を入力する入力手
段とよりなるプロセス状態解析方式において、前記電子
計算機は、前記プロセスデータの時系列データを重複時
間帯を持たせて複数分割し、分割したそれぞれの時系列
データをＦＦＴ解析し、その解析結果によりハンチング
の状態を判断する解析手段を備えていることを特徴とす
るプロセス状態解析方式。
【請求項２】解析手段は、分割数Ｎに１を加算したＮ
＋１回のＦＦＴ解析を行ない、第１回のＦＦＴ解析をす
べての時系列データに行なってハンチングの有無を判断
することを特徴とする請求項１記載のプロセス状態解析
方式。
【請求項３】解析手段は、それぞれのＦＦＴ解析結果
を一次回帰処理し、ハンチングの状態を発散傾向、収束
傾向又は安定状態の状態変数に置き換えることを特徴と
する請求項１又は２記載のプロセス状態解析方式。
【請求項４】表示手段は、それぞれの状態変数を、表
示又は印字で出力する手段であることを特徴とする請求
項１記載のプロセス状態解析方式。
【請求項５】解析手段は、ＦＦＴ解析を周期的に実行
し、それぞれのＦＦＴ解析結果を一次回帰処理し、将来
のハンチングの状態を予測することを特徴とする請求項
１，２又は３記載のプロセス状態解析方式。
【請求項６】表示手段は、マンマシン機構とＣＲＴよ
りなり、該マンマシン機構を経由して状態変数がＣＲＴ
に表示されることを特徴とする請求項１又は４記載のプ
ロセス状態解析方式。