JPH0425563B2 - - Google Patents

Info

Publication number
JPH0425563B2
JPH0425563B2 JP57180192A JP18019282A JPH0425563B2 JP H0425563 B2 JPH0425563 B2 JP H0425563B2 JP 57180192 A JP57180192 A JP 57180192A JP 18019282 A JP18019282 A JP 18019282A JP H0425563 B2 JPH0425563 B2 JP H0425563B2
Authority
JP
Japan
Prior art keywords
value
control
control system
loop
determined
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.)
Expired - Lifetime
Application number
JP57180192A
Other languages
Japanese (ja)
Other versions
JPS5969813A (en
Inventor
Sakae Tezuka
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.)
JFE Steel Corp
Original Assignee
Kawasaki Steel 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 Kawasaki Steel Corp filed Critical Kawasaki Steel Corp
Priority to JP57180192A priority Critical patent/JPS5969813A/en
Publication of JPS5969813A publication Critical patent/JPS5969813A/en
Publication of JPH0425563B2 publication Critical patent/JPH0425563B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • G05B23/0205Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
    • G05B23/0218Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
    • G05B23/0224Process history based detection method, e.g. whereby history implies the availability of large amounts of data
    • G05B23/024Quantitative history assessment, e.g. mathematical relationships between available data; Functions therefor; Principal component analysis [PCA]; Partial least square [PLS]; Statistical classifiers, e.g. Bayesian networks, linear regression or correlation analysis; Neural networks

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Measurement Of Resistance Or Impedance (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Description

【発明の詳細な説明】[Detailed description of the invention] 【産業上の利用分野】[Industrial application field]

本発明は、制御系の診断方法及び装置に係り、
特に、自動制御装置等の状態が正常であるか異常
であるかを診断する際に用いるのに好適な、制御
系の診断方法及び装置に関する。
The present invention relates to a method and device for diagnosing a control system,
In particular, the present invention relates to a control system diagnostic method and device suitable for use in diagnosing whether the state of an automatic control device or the like is normal or abnormal.

【従来の技術】[Conventional technology]

近年、製造設備等は、大規模且つ高級となり、
各機器の健全性が、それらの機器単独の生産性の
みならず、膨大な機器よりなる製造設備全体の生
産性をも左右する傾向にある。従つて、製造設備
全体及び各構成機器の健全性を確認し、異常の前
兆を早急に把握して、適切な処置を施す必要があ
る。 このような目的で、制御系の設定信号を階段状
に変え、それに対応する制御量の変化を観測す
る、いわゆるステツプ応答を調べる方法や、設定
信号として正弦波を印加し、その周波数を順次変
えて、それに対応した制御量の変化を観測する、
いわゆる周波数応答を調べる方法が実用化されて
いる。
In recent years, manufacturing equipment has become larger and more sophisticated.
The health of each piece of equipment tends to affect not only the productivity of each piece of equipment alone, but also the productivity of the entire manufacturing facility, which is made up of a huge number of pieces of equipment. Therefore, it is necessary to confirm the health of the entire manufacturing facility and each component, to quickly identify signs of abnormality, and to take appropriate measures. For this purpose, there is a method to examine the so-called step response, in which the setting signal of the control system is changed stepwise and the corresponding change in the controlled variable is observed, or a sine wave is applied as the setting signal and its frequency is sequentially changed. and observe the corresponding change in the control amount.
A method of examining so-called frequency response has been put into practical use.

【発明が解決しようとする課題】[Problem to be solved by the invention]

しかしながら、これらのステツプ応答或いは周
波数応答を調べる方法は、いずれも、 (1) 操業中の信号とは全く別のステツプ信号又は
正弦波信号を用いる必要があるため、操業中に
実施することは不可能である。 (2) これらの入力信号は、操業中の信号とは、そ
の形状、振幅等が異なることから、操業中の制
御系の状態を診断する手段としては不十分であ
る。 (3) 多重の制御ループよりなる制御系において
は、被観測ループ以外の制御ループとの干渉に
より、被観測ループの状態を正確に診断するこ
とが不可能である。 等の欠点を有していた。 即ち、例えば第1図に示すような、メインルー
プの調節器Aと、マイナーループの調節器Bと、
マイナーループの制御対象Cと、メインループの
制御対象D及びEと、マイナーループの帰還要素
Fと、メインループの帰還要素Gと、並列ループ
への影響要素Hと、並列ループの調節器Iと、並
列ループの制御対象Jと、並列ループからの影響
要素Kと、並列ループの制御対象Lと、並列ルー
プの帰還要素Mとを有する制御系を考えると、要
素B、C、Fよりなるマイナーループは、それ自
体で一つの閉じた制御ループであると同時に、
B・C/(1+B・C・F)というブロツクとし
て、メインループの制御対象の一つになつてい
る。従つて、マイナーループの各要素B、C、F
の特性は、マイナーループの制御性に影響するの
みならず、メインループの制御性にも重大な影響
を及ぼすことになる。 第1図において、aは、メインループの目標
値、bは、メインループの制御偏差(=a−n)、
cは、メインループの操作量(=マイナーループ
の目標値)、dは、メインループの外乱、eは、
マイナーループの制御偏差(=c+d−m)、f
は、マイナーループの操作量、gは、マイナール
ープの外乱、h=f+g、iは、マイナーループ
の制御量、j=i・D、k=j・E、1は、メイ
ンループの制御量、mは、マイナーループの帰還
量、nは、メインループの帰還量、oは、メイン
ループと並列で且つ互いに干渉し合うループ(以
下並列ループと称する)への干渉外乱、pは、並
列ループの目標値、qは、並列ループの制御偏差
(=p−w)、rは、並列ループの操作量、s=
r・J、tは、並列ループからの干渉外乱、u=
s・L、vは、並列ループの制御量、wは、並列
ループの帰還量である。 ここで、メインループの制御性を診断するため
に、メインループの目標値aを変化させた時、そ
れに対する応答は、メインループの主な要素であ
るA、D、E、Gのみでなく、マイナーループの
要素B、C、Fの影響も受ける。又、マイナール
ープの制御性を診断するために、マイナーループ
の目標値cを変化させ、それに対応するマイナー
ループの帰還量mを観測しようとする時、マイナ
ーループの制御量iは、メインループの要素D、
E、G、Aを経由して、マイナーループの目標値
cを変化させてしまう。更に、メインループの制
御量1は、メインループの内部状態量kの他に、
メインループの内部状態量jにH・I・J・K・
M/(1+I・J・L・M)が乗ぜられた信号t
を含むことから、I、J、L、Mよりなる並列ル
ープにも影響され、メインループの主な要素であ
るA、D、E、Gのみでは定まらない。 一方、制御系の状態を診断するための他の方法
として、或る要素の入力に白色雑音信号又は可変
周波数の正弦波信号を印加し、対応する出力の振
幅と位相角の変化を伝達関数として測定する方法
がある。 しかしながら、この方法は、対象とするが要素
が線形で1入力1出力の場合にしか適用できず、
実際の制御系に、この方法を適用すると、線形の
1入力1出力系に近似したことによる誤差が発生
する。 即ち、例えば前出第1図に示したような、2入
力2出力系においては、メインループの制御量
(出力)1は、メインループの目標値(入力)a
のみならず、並列ループの目標値(他の入力)p
の影響も受け、更に、外乱d、gにも影響され
る。同様に、並列ループの制御量(出力)vは、
2入力a,pと2外乱d,gに影響される。従つ
て、このような系では、或る特定の入力と出力と
の間の関係のみを抽出しようとしても、他の系の
影響を排除し得ないために、抽出不可能である。 又、公知のように、伝達関数法はラブラス変換
を基礎としているので、非線形系に対しては適用
できない。実際の制御系の伝達関数の変化を解析
する上で、線形要素の特性の変化によるものか、
非線形性の変化によるものかを区別しなければな
らないが、従来の伝達関数法では不可能である。 このように、従来の方法は、いずれも重大な欠
点があり、制御系の正常性を適確に診断すること
は困難であつた。 本発明は、前記従来の欠点を解消するべくなさ
れたもので、制御系の正常性を、容易に、且つ、
適確に診断することができる制御系の診断方法及
び装置を提供することを目的とする。
However, all of these methods of examining the step response or frequency response require the use of step signals or sine wave signals that are completely different from the signals used during operation; therefore, they cannot be carried out during operation. It is possible. (2) Since these input signals differ in shape, amplitude, etc. from signals during operation, they are insufficient as a means for diagnosing the state of the control system during operation. (3) In a control system consisting of multiple control loops, it is impossible to accurately diagnose the state of the observed loop due to interference with control loops other than the observed loop. It had the following drawbacks. That is, for example, as shown in FIG. 1, a main loop regulator A, a minor loop regulator B,
Controlled object C of the minor loop, controlled objects D and E of the main loop, feedback element F of the minor loop, feedback element G of the main loop, element H influencing the parallel loop, and regulator I of the parallel loop. , a control system having a parallel loop controlled object J, an influence element K from the parallel loop, a parallel loop controlled object L, and a parallel loop feedback element M, a minor consisting of elements B, C, and F is considered. The loop is itself a closed control loop, and at the same time
The block B.C/(1+B.C.F) is one of the control targets of the main loop. Therefore, each element B, C, F of the minor loop
The characteristics of not only affect the controllability of the minor loop but also have a significant effect on the controllability of the main loop. In FIG. 1, a is the target value of the main loop, b is the control deviation of the main loop (=an),
c is the operation amount of the main loop (=target value of the minor loop), d is the disturbance of the main loop, and e is the
Minor loop control deviation (=c+d-m), f
is the operating amount of the minor loop, g is the disturbance of the minor loop, h=f+g, i is the controlled amount of the minor loop, j=i・D, k=j・E, 1 is the controlled amount of the main loop, m is the amount of feedback of the minor loop, n is the amount of feedback of the main loop, o is the interference disturbance to the loop that is parallel to the main loop and interferes with each other (hereinafter referred to as parallel loop), p is the amount of feedback of the parallel loop. The target value, q is the control deviation of the parallel loop (=p-w), r is the manipulated variable of the parallel loop, s=
rJ, t is interference disturbance from parallel loop, u=
s·L, v is the control amount of the parallel loop, and w is the feedback amount of the parallel loop. Here, in order to diagnose the controllability of the main loop, when the target value a of the main loop is changed, the response to it is not only A, D, E, and G, which are the main elements of the main loop, but also It is also affected by minor loop elements B, C, and F. Furthermore, in order to diagnose the controllability of the minor loop, when changing the target value c of the minor loop and observing the corresponding feedback amount m of the minor loop, the control amount i of the minor loop is the same as that of the main loop. element D,
The target value c of the minor loop changes via E, G, and A. Furthermore, the control amount 1 of the main loop includes, in addition to the internal state amount k of the main loop,
The internal state quantity j of the main loop is H・I・J・K・
Signal t multiplied by M/(1+I・J・L・M)
, it is also affected by the parallel loop consisting of I, J, L, and M, and cannot be determined only by the main elements of the main loop, A, D, E, and G. On the other hand, as another method for diagnosing the state of a control system, a white noise signal or a variable frequency sine wave signal is applied to the input of a certain element, and changes in the amplitude and phase angle of the corresponding output are used as a transfer function. There is a way to measure it. However, this method can only be applied when the target element is linear and has one input and one output.
When this method is applied to an actual control system, errors occur due to approximation to a linear one-input, one-output system. That is, for example, in a two-input, two-output system as shown in FIG.
as well as the target value (other input) of the parallel loop p
It is also affected by disturbances d and g. Similarly, the control amount (output) v of the parallel loop is
It is affected by two inputs a and p and two disturbances d and g. Therefore, in such a system, even if an attempt is made to extract only the relationship between a certain specific input and output, it is impossible to extract it because the influence of other systems cannot be excluded. Furthermore, as is well known, the transfer function method is based on the Labras transformation, and therefore cannot be applied to nonlinear systems. When analyzing changes in the transfer function of an actual control system, it is important to check whether it is due to changes in the characteristics of linear elements.
It is necessary to distinguish whether this is due to a change in nonlinearity, but this is not possible with the conventional transfer function method. As described above, all of the conventional methods have serious drawbacks, and it has been difficult to accurately diagnose the normality of the control system. The present invention was made in order to eliminate the above-mentioned conventional drawbacks, and it is possible to easily maintain the normality of the control system.
It is an object of the present invention to provide a control system diagnostic method and device that can accurately diagnose a control system.

【課題を解決するための手段】[Means to solve the problem]

本発明は、制御系の診断方法において、まず、
制御系の正常性を診断するに際して、制御系の制
御偏差の2乗積分値を演算し、該2乗積分値が、
予め設定された基準値より大である時に、制御ル
ープとしての制御性不良と判定し、制御系の各要
素について、該要素の入力の出力に対する寄与度
を表わすコヒーレンス関数を演算し、該コヒーレ
ンス関数の値が、予め設定された基準値より大で
ある場合には、該要素の入力に対する出力の関係
を伝達関数として求め、該伝達関数のゲイン余裕
や位相余裕が、予め設定されたそれぞれの基準範
囲外である時に、該要素の制御特性不良と判定
し、一方、前記コヒーレンス関数の値が、前記基
準値以下である場合には、該要素の制御特性の良
否を判定しないことによつて、前記目的を達成し
たものである。 又、本発明は、前記目的を達成する装置を、制
御系の制御偏差の2乗積分値を演算する演算手段
と、前記2乗積分値を、予め設定された基準値と
比較する第1の比較手段と、制御系の各要素につ
いて、該要素の入力の出力に対する寄与度を表わ
すコヒーレンス関数及び該要素の入力に対する出
力の関係を表わす伝達関数を求める周波数解析手
段と、前記コヒーレンス関数の値を、予め設定さ
れた基準値と比較する第2の比較手段と、前記伝
達関数を、予め設定されたそれぞれの基準値と比
較して、ゲイン余裕及び位相余裕を求め第3の比
較手段と、前記各比較手段の出力によつて、前記
2乗積分値がその基準値より大である時は、制御
ループとしての制御性不良と判定し、前記コヒー
レンス関数の値がその基準値より大である場合
に、前記伝達関数のゲイン余裕と位相余裕の値
が、予め設定されたそれぞれの基準範囲外である
時は、該要素の制御特性不良と判定し、前記コヒ
ーレンス関数の値が、前記基準値以下である時
は、該要素の制御特性の良否を判定しない論理回
路とを用いて構成したものである。
The present invention provides a method for diagnosing a control system, first of all,
When diagnosing the normality of the control system, the square integral value of the control deviation of the control system is calculated, and the square integral value is
When the value is larger than a preset reference value, it is determined that the controllability of the control loop is poor, and for each element of the control system, a coherence function representing the contribution of the input of the element to the output is calculated, and the coherence function is calculated. If the value of is larger than a preset reference value, the relationship between the input and output of the element is determined as a transfer function, and the gain margin and phase margin of the transfer function are determined according to each preset reference value. When the value of the coherence function is outside the range, it is determined that the control characteristics of the element are poor; on the other hand, when the value of the coherence function is less than or equal to the reference value, the control characteristics of the element are not determined to be good or bad; The above objective has been achieved. Further, the present invention provides an apparatus for achieving the above object, comprising: a calculation means for calculating a square integral value of a control deviation of a control system; and a first calculating means for calculating a square integral value of a control deviation of a control system; a comparison means, a frequency analysis means for determining, for each element of the control system, a coherence function representing the degree of contribution of the input of the element to the output and a transfer function representing the relationship between the input and the output of the element; , a second comparison means for comparing with a preset reference value; a third comparison means for comparing the transfer function with each preset reference value to obtain a gain margin and a phase margin; According to the output of each comparison means, when the square integral value is larger than the reference value, it is determined that the controllability of the control loop is poor, and when the value of the coherence function is larger than the reference value. When the values of the gain margin and phase margin of the transfer function are outside the respective preset reference ranges, it is determined that the control characteristics of the element is poor, and the value of the coherence function is determined to be less than or equal to the reference value. When this is the case, it is constructed using a logic circuit that does not judge whether the control characteristics of the element are good or bad.

【作用】[Effect]

本発明は、 (1) 線形性や入出力数に関係なく、制御偏差の2
乗積分値で制御性が評価できること、 (2) 伝達関数の信頼性をコヒーレンス関係で評価
できること、 に着目してなされたものである。 以下、前出第1図に示したような、メインレー
プの中にマイナーループを含み、更に、他に並列
ループを有する制御ループを例にとつて、本発明
の方法を説明する。 診断しようとするメインループは、その目標値
がa、制御量が1、制御偏差がb、操作量がcで
あり、操作量cに対する制御量1の伝達関数をN
とすると、第2図に示す如く、調節器Aと、等価
制御対象Nと、帰還要素Gとからなる制御ループ
として表現することができる。ここで、等価制御
対象Nは、1入力1出力の単純な要素ではなく、
入力cの他に、外乱d、g、pが加わるとみなす
ことができる。 この制御ループ全体を閉ループ10とすると、
この閉ループ10の入出力a、l、状態量b、
c、n、外乱d、g、pを、第2図に示す如く、
診断装置12に入力する。この診断装置12は、
例えば第3図に示す如く、制御系の制御偏差bの
2乗積分値ISEを演算する演算器14と、前記2
乗積分値ISEを基準ISE0と比較する第1の比較器
16と、制御系の各要素のコヒーレンス関数CF
と伝達関数を求める2次元周波数解析器18と、
前記コヒーレンス関数CFの値を基準値CF0と比
較する第2の比較器20と、前記伝達関数を基準
値と比較して、ゲイン余裕及び位相余裕を求める
第3の比較器22と、前記各比較器16,20,
22の出力によつて制御系の正常性を判定するた
めの論理回路24とから構成されている。 前記のような診断装置12における制御系の正
常性の診断は、次のようにして行われる。 即ち、まず、制御ループとしての制御性を、制
御偏差bの2乗積分値ISEとして評価する。 ISE=∫〓0b2dt……(1) ここで、積分時間τは、制御ループの時定数に
対応して設定されている。 制御が緩慢な場合、及び、過敏な場合のいずれ
の場合にも、制御偏差bの絶対値は大きく、更
に、操業面に対する影響の度合から、制御偏差b
の2乗で評価するのが適している。なお、過渡的
な応答を評価する際には、その過渡現象の経過時
間にわたつて、制御偏差bの二乗値を積分して評
価するのが適している。 2乗積分値ISEが、そのループの基準値ISE0
りも小さい場合は、制御ループとしての制御性に
問題はないが、基準値ISE0以上となつた場合は、
各要素の特性を把握し、修理や取替えを含む必要
な調整を実施して、正常な状態に復旧する。 制御ループとしての制御性が把握できても、こ
の中の各要素の特性を正確に把握することが、正
しい運転のために必要となる。このために以下に
示す伝達関数によつて評価する。この際、伝達関
数の信頼度を高めるため、コヒーレンス関数で信
頼の置ける角周波数域において、伝達関数による
評価を行うこともできる。 具体的には、まず、目標値aと制御量1とを2
次元周波数解析器18で2次元周波数解析し、両
者間のコヒーレンス関数CFalを求める。このコ
ヒーレンス関数CFalは、入力aの出力1に対す
る寄与度を表わす関数であり、数学的には次式で
求められる。 CFa1=|C(ω)|2/{I(ω)・0(ω)2}……
(2) ここではI(ω)、0(ω)は、それぞれ入力信
号a、出力信号1のオートパワースペクトル、C
(ω)は、入力と出力のクロスパワースペクトル、
ωは角周波数である。 前出(2)式によつて求められたコヒーレンス関数
CFalが大である程、制御量1は目標値aに強く
影響されていることを示し、一方、コヒーレンス
関数CFalが小さい場合は、制御量1は、目標値
a以外の入力(又は外乱)に強く影響されている
か、又は、目標値aから制御量1に至る間の要素
が強い非線形性を有していることを示している。 従つて、コヒーレンス関数CFalが基準値CFal0
より大きい場合、即ち、制御量1が目標値aに強
く対応している場合には、目標値a(入力)に対
する制御量l(出力)の関係を伝達関数として求
める。一方、コヒーレンス関数CFalが基準値
CFal0以下でいる場合には、制御量1が目標値a
以外の他の要因に支配されていることになるの
で、その伝達関数の信頼性が低い。従つて、該要
素(目標値a)の制御特性の良否について判定す
ることなく、例えば、制御量lに強い影響を与え
ている他人の入力、外乱又は要素の制御特性を調
査することができる。 コヒーレンス関数CFalが大きく、伝達関数を
求める時は、制御ループの余裕で制御性を評価す
る。 例えば、操業を停止させて診断する場合には、
目標値印加点に模擬信号として白色雑音、又は、
可変周波数の正弦波信号を印加すると共に、2次
元周波数解析器18の一方の入力端子に接続し、
更に、模擬信号に対応する制御量1を2次元周波
数解析器18の他方の入力端子に接続して、次式
で数学的に表わされる、2入力間の伝達関数と、
前記2乗積分値を、予め設定された基準値と比較
する第1の比較手段と、制御系の各要素につい
て、該要素の入力の出力に対する寄与度を表わす
コヒーレンス関数及び該要素の入力に対する出力
の関係を表わす伝達関数を求める周波数解析手段
と、前記コヒーレンス関数の値を、予め設定され
た基準値と比較する第2の比較手段と、前記伝達
関数を、予め設定されたそれぞれの基準値と比較
して、ゲイン余裕及び位相余裕を求める第3の比
較手段と、前記各比較手段の出力によつて、前記
2乗積分値がその基準値より大である時は、制御
ループとしての制御性不良と判定し、前記コヒー
レンス関数の値がその基準値より大である場合
に、前記伝達関数のゲイン余裕と位相余裕の値
が、予め設定されたそれぞれの基準範囲外である
時は、該要素の制御特性不良と判定し、前記コヒ
ーレンス関数の値が、前記基準値以下である時
は、該要素の制御特性の良否を判定しない論理回
路とを用いてを求める。 G=C(ω)/I(ω)……(3) 一方、操業中に操業に外乱を与えずに伝達関数
を測定する場合は、模擬信号を印加せず、測定し
ようとする要素の入力と出力を2次元周波数解析
器18の2入力端子に接続し求める。 このようにして求められた伝達関数から、制御
性を診断するために、次の2つの方法がある。 (1) ゲイン|G|が0dBとなる周波数f1における
位相θf1の−180°に対する余裕、即ち、 位相余裕Me=θf1+180°と、位相θが−180°
になる周波数f2におけるゲイン|G|f2dBに対
する余裕、即ち、ゲイン余裕MG=|G|f2が、
それぞれ基準の値M〓O、MGOに対して、どれだ
け離れているかによつて判定する方法。 この方法は、制御系の特性を正確に診断する
ために有効であるが、広い周波数領域にわたつ
て有意な信号を得る必要があることから、操業
中に用いることは不適切である。 (2) ゲイン|G|が基準の低周波数域の値|G0
|から3dB低下した時の周波数f−3dBと、位
相θが基準の低周波数域の値θ0から90°遅れた
時の周波数f90°とが、それぞれの基準値(f−
3dB)0、(f90°)0からどれだけずれているかを診
断する方法。 この方法、精度的には多少不安定であるが、
比較的狭い周波数領域で十分評価し得る。 前記のような方法で、各要素の特性を定量的に
評価し、装置としての正常性を診断することがで
きる。第4図A〜Cに、本発明による、伝達関
数、コヒーレンス関数、クロスパワースペクトル
及びオートパワースペクトルの測定波形の例を示
す。 なお、前記コヒーレンス関数や伝達関数は、前
出第2図に示す如く、ブラウン管26に表示して
観察することや、或いは、ドツトプリンタやデジ
タルブロツク等のハードコピー機28でコピーす
ることが可能である。又、2乗積分値ISEは、数
字で表示したり、コピーしたりすることが可能で
あり、或いは、図示しないデジタル/アナログ変
換器と指示記録計等によつてアナログ的に指示ま
たは記録することも可能である。更に、電子計算
機を活用することにより、多数の測定点を順次捜
査し、それぞれの要素の特性を自動的に測定する
ことも可能である。
The present invention has the following features: (1) Regardless of linearity or number of inputs and outputs, the control deviation is
This was done by focusing on the following points: (2) the reliability of the transfer function can be evaluated using the coherence relationship; and (2) the reliability of the transfer function can be evaluated using the coherence relationship. The method of the present invention will be explained below by taking as an example a control loop that includes a minor loop in the main loop and further has parallel loops, as shown in FIG. 1 above. The main loop to be diagnosed has a target value a, a controlled variable 1, a control deviation b, and a manipulated variable c, and the transfer function of the controlled variable 1 to the manipulated variable c is N.
Then, as shown in FIG. 2, it can be expressed as a control loop consisting of a regulator A, an equivalent controlled object N, and a feedback element G. Here, the equivalent controlled object N is not a simple element with one input and one output, but
In addition to the input c, it can be considered that disturbances d, g, and p are added. If this entire control loop is defined as a closed loop 10, then
The input/output a, l of this closed loop 10, the state quantity b,
c, n, disturbances d, g, p, as shown in Figure 2,
input to the diagnostic device 12; This diagnostic device 12 is
For example, as shown in FIG.
A first comparator 16 that compares the product integral value ISE with the reference ISE 0 , and a coherence function CF of each element of the control system.
and a two-dimensional frequency analyzer 18 for determining the transfer function.
a second comparator 20 that compares the value of the coherence function CF with a reference value CF 0 ; a third comparator 22 that compares the transfer function with the reference value to obtain a gain margin and a phase margin; Comparators 16, 20,
22, and a logic circuit 24 for determining the normality of the control system based on the output of the control system 22. Diagnosis of the normality of the control system in the diagnostic device 12 as described above is performed as follows. That is, first, the controllability as a control loop is evaluated as the square integral value ISE of the control deviation b. ISE=∫〓 0 b 2 dt...(1) Here, the integration time τ is set corresponding to the time constant of the control loop. The absolute value of the control deviation b is large in both cases of slow and sensitive control, and furthermore, due to the degree of influence on the operation, the control deviation b
It is appropriate to evaluate using the square of . Note that when evaluating a transient response, it is appropriate to integrate and evaluate the square value of the control deviation b over the elapsed time of the transient phenomenon. If the square integral value ISE is smaller than the reference value ISE 0 of the loop, there is no problem with the controllability of the control loop, but if it becomes greater than the reference value ISE 0 ,
Understand the characteristics of each element and make necessary adjustments, including repairs and replacements, to restore normal conditions. Even if the controllability of the control loop can be understood, it is necessary to accurately understand the characteristics of each element within it for correct operation. For this purpose, evaluation is performed using the transfer function shown below. At this time, in order to increase the reliability of the transfer function, it is also possible to perform evaluation using the transfer function in an angular frequency range where the coherence function is reliable. Specifically, first, the target value a and the control amount 1 are set to 2.
A two-dimensional frequency analyzer 18 performs two-dimensional frequency analysis to obtain a coherence function CFal between the two. This coherence function CFal is a function representing the degree of contribution of input a to output 1, and is mathematically determined by the following equation. CFa1=|C(ω)| 2 /{I(ω)・0(ω) 2 }...
(2) Here, I(ω) and 0(ω) are the auto power spectra of input signal a and output signal 1, respectively, and C
(ω) is the input and output cross power spectrum,
ω is the angular frequency. Coherence function determined by equation (2) above
The larger CFal is, the more the controlled variable 1 is influenced by the target value a. On the other hand, when the coherence function CFal is small, the controlled variable 1 is influenced by inputs (or disturbances) other than the target value a. This indicates that the control amount is strongly influenced, or that the elements between the target value a and the control amount 1 have strong nonlinearity. Therefore, the coherence function CFal is the reference value CFal 0
If the value is larger, that is, if the controlled variable 1 strongly corresponds to the target value a, the relationship between the controlled variable l (output) and the target value a (input) is determined as a transfer function. On the other hand, the coherence function CFal is the reference value
If CFal is less than 0 , the control amount 1 is the target value a
Since the transfer function is dominated by other factors, the reliability of the transfer function is low. Therefore, without determining whether the control characteristics of the element (target value a) are good or bad, it is possible to investigate, for example, inputs from others, disturbances, or control characteristics of the element that have a strong influence on the control amount l. When the coherence function CFal is large and the transfer function is determined, controllability is evaluated using the margin of the control loop. For example, when diagnosing after stopping operations,
White noise as a simulated signal at the target value application point, or
Applying a variable frequency sine wave signal and connecting to one input terminal of the two-dimensional frequency analyzer 18,
Furthermore, the control amount 1 corresponding to the simulated signal is connected to the other input terminal of the two-dimensional frequency analyzer 18, and a transfer function between the two inputs is mathematically expressed by the following equation,
a first comparing means for comparing the square integral value with a preset reference value; and a coherence function representing the degree of contribution of the input of the element to the output for each element of the control system, and an output for the input of the element. frequency analysis means for determining a transfer function representing the relationship between the two; second comparison means for comparing the value of the coherence function with a preset reference value; and a second comparison means for comparing the value of the coherence function with a preset reference value. A third comparison means for comparing the gain margin and a phase margin, and the output of each of the comparison means, when the square integral value is larger than the reference value, the controllability as a control loop is determined. When the value of the coherence function is determined to be defective and the value of the coherence function is larger than its reference value, and the value of the gain margin and phase margin of the transfer function are outside the respective preset reference ranges, the element When the control characteristic of the element is determined to be poor and the value of the coherence function is less than or equal to the reference value, the control characteristic of the element is determined using a logic circuit that does not determine whether the control characteristic of the element is good or bad. G=C(ω)/I(ω)……(3) On the other hand, when measuring the transfer function without causing any disturbance to the operation during operation, do not apply a simulated signal and input the element to be measured. and the output are connected to the two input terminals of the two-dimensional frequency analyzer 18. There are two methods for diagnosing controllability from the transfer function obtained in this way. (1) The margin for the phase θf 1 of -180° at the frequency f 1 where the gain |G| becomes 0 dB, that is, the phase margin Me = θ f1 +180°, and the phase θ is -180°.
The margin for the gain |G|f 2 dB at the frequency f 2 that becomes, that is, the gain margin M G = |G|f 2 is
A method of judging based on how far away it is from the reference values M〓 O and M GO , respectively. Although this method is effective for accurately diagnosing the characteristics of a control system, it is inappropriate to use it during operation because it requires obtaining significant signals over a wide frequency range. (2) Gain |G | is the value in the low frequency range |G 0
The frequency f-3dB when the frequency has decreased by 3dB from
3dB) 0 , (f90°) How to diagnose how far it deviates from 0 . Although this method is somewhat unstable in terms of accuracy,
It can be sufficiently evaluated in a relatively narrow frequency range. By the method described above, the characteristics of each element can be quantitatively evaluated and the normality of the device can be diagnosed. 4A to 4C show examples of measured waveforms of a transfer function, coherence function, cross power spectrum, and auto power spectrum according to the present invention. The coherence function and transfer function can be displayed on a cathode ray tube 26 for observation, as shown in FIG. 2, or can be copied using a hard copy machine 28 such as a dot printer or a digital block. . In addition, the square integral value ISE can be displayed numerically or copied, or can be indicated or recorded in an analog manner using a digital/analog converter and an indicator recorder (not shown). is also possible. Furthermore, by utilizing a computer, it is also possible to sequentially investigate a large number of measurement points and automatically measure the characteristics of each element.

【実施例】【Example】

第5図に、本発明に係る診断方法の実施例の流
れ図を示す。この実施例においては、伝達関数の
評価方法として、余裕Mを用いた例を示している
が、勿論、f−3dB、f90°で評価することも可能
である。更に、第5図においては、それぞれの判
定において異常と判定された場合に実施すべき調
整等についても記載している。 以下、第5図に示した流れ図を詳細に説明す
る。 まず、ステツプ101において、閉ループ10
の入出力a、l、状態量b、c、n、外乱d、
g、pを読込む。次いで、ステツプ102に進み、
診断過程を表わすカウンタα、β、γを、すべて
リセツトする。次いで、ステツプ103に進み、前
出(1)式を用いて算出した制御偏差bの2乗積分値
ISEが、基準値ISE0未満であるか否かを判定す
る。判定結果が正である場合には、制御系は正常
であると判断して、診断を終了する。 一方、前出ステツプ103における判定結果が否
である場合には、ステツプ104に進み、カウンタ
αの計数値が1以下であるか否かを判定する。判
定結果が否である場合には、ステツプ105に進み、
他の要因を調査して、診断を終了する。 一方、前出ステツプ104における判定結果が正
である場合には、ステツプ106に進み、前出(2)式
により算出されるコヒーレンス関数CFalが、そ
の基準値CFaloより大であるか否かを判定する。
判定結果が正である場合には、ステツプ107に進
み、余裕Malが、基準値Malo±△Mal内にある
か否かを判定する。判定結果が正である場合に
は、ステツプ108に進み、基準値ISE0、CFalo、
Malo、△Malを見直す。次いで、ステツプ109に
進み、他の要因を調査して、診断を終了する。 一方、前出ステツプ107における判定結果が否
である場合には、ステツプ110に進み、カウンタ
γの計数値が1未満であるか否かを判定する。判
定結果が正である場合には、ステツプ111に進み、
制御量1と帰還量n間のコヒーレンス関数CFln
が、その基準値CFlnoより大であり、且つ、両者
間の余裕Mlnが、基準値Mlno±△Mln内にある
か否かを判定する。判定結果が否である場合に
は、ステツプ112に進み、帰還要素Gを調整する。
ステツプ112終了後、或いは、前出ステツプ111に
おける判定結果が正である場合には、ステツプ
113に進み、操作量cと制御量1間のコヒーレン
ス関数CFclが、その基準値CFcloより大であり、
且つ、両者間の余裕Mclが、基準値Mclo±△
Mcl内にあるか否かを判定する。判定結果が否か
である場合には、ステツプ114に進み、等価制御
対象Nの伝達関数を調整する。ステツプ114終了
後、或いは、前出ステツプ113における判定結果
が正である場合には、ステツプ115に進み、制御
偏差bと操作量c間のコヒーレンス関数CFbcが、
その基準値CFbcoより大であり、且つ、両者間の
余裕Mbcが、基準値Mbco±△Mbc内にあるか否
かを判定する。判定結果が否かである場合には、
ステツプ116に進み、調節器Aを調整する。ステ
ツプ116終了後、或いは、前出ステツプ115におけ
る判定結果が正である場合には、ステツプ117に
進み、カウンタγを1だけカウントアツプして、
前出ステツプ107に戻る。 一方、前出ステツプ110における判定結果が否
である場合には、ステツプ118に進み、基準値
CFlno、Mlno、△Mln、CFclo、Mclo、△Mcl、
CFbco、Mbco、△Mbcを見直し、更に、ステツ
プ119で、他の要因を調査して、前出ステツプ106
に戻る。 又、前出ステツプ106における判定結果が否で
ある場合には、ステツプ120に進み、カウンタβ
の計数値が1未満であるか否かを判定する。判定
結果が否である場合には、ステツプ121に進み、
他の要因を調査した後、ステツプ122で、カウン
タαを1だけカウントアツプして前出ステツプ
103に戻る。 一方、前出ステツプ120における判定結果が正
である場合には、ステツプ123に進み、外乱dと
制御量1間のコヒーレンス関数CFdl、外乱gと
制御量1間のコヒーレンス関数CFgl、外乱pと
制御量1間のコヒールレンス関数CFplが、それ
ぞれの基準値CFdlo、CFglo、CFploより大であ
るか否かを判定する。判定結果が正である場合に
は、ステツプ124に進み、外乱d、g、pを低減
させる。一方、前出ステツプ123の判定結果が否
である場合には、ステツプ125に進み、帰還要素
G、等価制御対象N、調節器Aを調節する。該ス
テツプ125、或いは、前出ステツプ124終了後、ス
テツプ126に進み、カウンタβを1だけカウント
アツプして、前出ステツプ106に戻る。 なお前記実施例においては、メインループを診
断する場合について説明しているが、メインルー
プの診断の一環として、又は、全くの別の目的か
らマイナーループや並列ループを診断する場合に
も、メインループについて説明したのと同様の手
順で診断すればよい。 次に、厚板圧延機の油圧式自動板厚制御系統
(以下油圧AGC系統と称する)において、制御状
態が過敏になつた時の診断例について説明する。 油圧AGC系統は、アナログ制御系よりなり、
油圧シリンダ位置制御系の中に、油圧サーボ弁開
度制御系を含み、更に、油圧シリンダ位置制御系
と並列に自動板厚制御系が設けられている。この
制御系において、油圧シリンダ位置制御系は、自
動板厚制御系と干渉し合つているため、コヒーレ
ンス関数が小さい。油圧サーボ弁開度制御系は、
単一の制御系であり、更に、殆んど線形系とみな
し得ることから、コヒーレンス関数が大きい。 従つて、診断に際しては、まず、制御系の制御
偏差を周期的に測定して、その測定値を2乗演算
器で平行演算し、それを加算器で加算した。加算
器は、設定加算回数(この時は50回)毎にリセツ
トされる。この加算値をデジタルプリンタで繰返
し印字した。すると、この加算値(ISE)が
7857/50回となり、基準値1000/50回よりも大で
あつたため、制御系の目標値aと制御量1とのコ
ヒーレンス関数CFalを2次元周波数解析器で求
めたところ、常用の周波数範囲0〜50Hzにおい
て、コヒーレンス関数CFalが基準値0.9以上であ
つた。従つて、目標値aと制御量1、操作量cと
の間の伝達関数を、2次元周波数解析器で求めた
ところ、調節器のゲイン不足(基準値の70%)が
判明した。 これに対して、調節器のゲインを基準値の102
%に調整したところ、制御性は飛躍的に向上し、
ISEは0.048/50回と正常値に戻つた。
FIG. 5 shows a flowchart of an embodiment of the diagnostic method according to the present invention. In this embodiment, an example is shown in which the margin M is used as a method for evaluating the transfer function, but it is of course also possible to evaluate using f-3 dB and f90°. Further, FIG. 5 also describes adjustments, etc. that should be performed when each determination is determined to be abnormal. The flowchart shown in FIG. 5 will be explained in detail below. First, in step 101, the closed loop 10
input/output a, l, state quantities b, c, n, disturbance d,
Read g and p. Next, proceed to step 102,
All counters α, β, and γ representing the diagnostic process are reset. Next, the process proceeds to step 103, where the square integral value of the control deviation b calculated using the above equation (1) is calculated.
It is determined whether ISE is less than the reference value ISE0 . If the determination result is positive, it is determined that the control system is normal, and the diagnosis is terminated. On the other hand, if the determination result in step 103 is negative, the process proceeds to step 104, where it is determined whether the count value of counter α is 1 or less. If the determination result is negative, proceed to step 105.
Investigate other factors to finalize the diagnosis. On the other hand, if the determination result in step 104 is positive, the process proceeds to step 106, where it is determined whether the coherence function CFal calculated by equation (2) is greater than its reference value CFalo. do.
If the determination result is positive, the process proceeds to step 107, where it is determined whether the margin Mal is within the reference value Malo±ΔMal. If the judgment result is positive, proceed to step 108 and set the reference values ISE 0 , CFalo,
Review Malo, △Mal. Next, the process proceeds to step 109 to investigate other factors and complete the diagnosis. On the other hand, if the determination result in step 107 is negative, the process proceeds to step 110, where it is determined whether the count value of counter γ is less than 1. If the judgment result is positive, proceed to step 111.
Coherence function CFLn between control amount 1 and feedback amount n
is larger than the reference value CFLno, and whether or not the margin Mln between the two is within the reference value Mlno±ΔMln. If the determination result is negative, the process proceeds to step 112 and the feedback element G is adjusted.
After step 112 is completed, or if the judgment result in step 111 is positive, step
Proceeding to step 113, the coherence function CFcl between the manipulated variable c and the controlled variable 1 is larger than its reference value CFclo,
Moreover, the margin Mcl between the two is the reference value Mcl±△
Determine whether it is within Mcl. If the result of the determination is negative, the process proceeds to step 114, where the transfer function of the equivalent controlled object N is adjusted. After step 114, or if the judgment result in step 113 is positive, the process proceeds to step 115, where the coherence function CFbc between the control deviation b and the manipulated variable c is
It is determined whether the reference value CFbco is greater than the reference value CFbco and the margin Mbc between the two is within the reference value Mbco±ΔMbc. If the judgment result is negative,
Proceed to step 116 and adjust regulator A. After step 116, or if the judgment result in step 115 is positive, proceed to step 117, increment the counter γ by 1, and
Return to step 107 above. On the other hand, if the judgment result in step 110 is negative, the process proceeds to step 118 and the reference value is
CFlno, Mlno, △Mln, CFclo, Mclo, △Mcl,
Review CFbco, Mbco, and △Mbc, and further investigate other factors in step 119, and repeat step 106 above.
Return to If the judgment result in step 106 is negative, the process proceeds to step 120, where the counter β is
It is determined whether the count value of is less than 1. If the determination result is negative, proceed to step 121.
After investigating other factors, in step 122, count up the counter α by 1 and repeat the previous step.
Return to 103. On the other hand, if the judgment result in step 120 is positive, the process proceeds to step 123, where the coherence function CFdl between the disturbance d and the control amount 1, the coherence function CFgl between the disturbance g and the control amount 1, the coherence function CFgl between the disturbance p and the control It is determined whether the coherence function CFpl between the quantities 1 is larger than the respective reference values CFdlo, CFglo, and CFplo. If the determination result is positive, the process proceeds to step 124 and the disturbances d, g, and p are reduced. On the other hand, if the determination result in step 123 is negative, the process proceeds to step 125, where the feedback element G, equivalent control object N, and regulator A are adjusted. After step 125 or step 124 described above is completed, the process proceeds to step 126, increments the counter β by 1, and returns to step 106 described above. Although the above embodiment describes the case of diagnosing the main loop, the main loop can also be used as part of diagnosing the main loop or when diagnosing minor loops or parallel loops for completely different purposes. The diagnosis can be made using the same procedure as explained above. Next, an example of diagnosis when the control state becomes oversensitive in a hydraulic automatic plate thickness control system (hereinafter referred to as a hydraulic AGC system) of a thick plate rolling mill will be described. The hydraulic AGC system consists of an analog control system.
The hydraulic cylinder position control system includes a hydraulic servo valve opening control system, and is further provided with an automatic plate thickness control system in parallel with the hydraulic cylinder position control system. In this control system, the hydraulic cylinder position control system interferes with the automatic plate thickness control system, so the coherence function is small. The hydraulic servo valve opening control system is
Since it is a single control system and can be regarded as an almost linear system, it has a large coherence function. Therefore, in the diagnosis, first, the control deviation of the control system was periodically measured, the measured values were subjected to parallel calculations using a square calculator, and the results were added using an adder. The adder is reset every set number of additions (50 times in this case). This added value was repeatedly printed using a digital printer. Then, this added value (ISE) becomes
The result was 7857/50 times, which was larger than the reference value of 1000/50 times, so when the coherence function CFal between the target value a of the control system and the controlled variable 1 was determined using a two-dimensional frequency analyzer, it was found that the commonly used frequency range was 0. At ~50Hz, the coherence function CFal was above the reference value of 0.9. Therefore, when the transfer function between the target value a, the controlled variable 1, and the manipulated variable c was obtained using a two-dimensional frequency analyzer, it was found that the gain of the regulator was insufficient (70% of the reference value). In contrast, the gain of the regulator is set to the reference value of 102
When adjusted to %, controllability improved dramatically,
ISE returned to normal value at 0.048/50 times.

【発明の効果】【Effect of the invention】

以上説明した通り、本発明によれば、制御系な
異常を、容易に、且つ、適確に診断することがで
きるという優れた効果を有する。
As described above, the present invention has the excellent effect of being able to easily and accurately diagnose abnormalities in the control system.

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

第1図は、本発明に係る制御系の診断方法が適
用される、複雑な制御ループの1例を示すブロツ
ク線図、第2図は、本発明に係る診断装置を用い
て、前出第1図に示したような制御対象のメイン
ループを診断している状態を示すブロツク線図、
第3図は、前記診断装置の構成を示すブロツク線
図、第4図A〜Cは、本発明による、伝達関数、
コヒーレンス関数、クロスパワースペクトル及び
オートパワースペクトルの測定波形の例を示す線
図、第5図は、本発明に係る診断方法の実施例を
示す流れ図である。 10……閉ループ、A……調節器、N……等価
制御対象、G……帰還要素、12……診断装置、
14……演算器、16,20,22……比較器、
18……2次元周波数解析器、24……論理回
路。
FIG. 1 is a block diagram showing an example of a complicated control loop to which the control system diagnosis method according to the present invention is applied, and FIG. A block diagram showing the state in which the main loop of the controlled object as shown in Figure 1 is being diagnosed,
FIG. 3 is a block diagram showing the configuration of the diagnostic device, and FIGS. 4A to 4C show transfer functions according to the present invention,
FIG. 5 is a diagram showing examples of measurement waveforms of a coherence function, a cross power spectrum, and an auto power spectrum. FIG. 5 is a flowchart showing an embodiment of the diagnostic method according to the present invention. 10...Closed loop, A...Adjuster, N...Equivalent control object, G...Feedback element, 12...Diagnostic device,
14... Arithmetic unit, 16, 20, 22... Comparator,
18...Two-dimensional frequency analyzer, 24...Logic circuit.

Claims (1)

【特許請求の範囲】 1 制御系の正常性を診断するに際して、 制御系の制御偏差の2乗積分値を演算し、 該2乗積分値が、予め設定された基準値より大
である時に、制御ループとしての制御性不良と判
定し、 制御系の各要素について、該要素の入力の出力
に対する寄与度を表わすコヒーレンス関数を演算
し、 該コヒーレンス関数の値が、予め設定された基
準値より大である場合には、該要素の入力に対す
る出力の関係を伝達関数として求め、 該伝達関数のゲイン余裕や位相余裕が、予め設
定されたそれぞれの基準範囲外である時に、該要
素の制御特性不良と判定し、 一方、前記コヒーレンス関数の値が、前記基準
値以下である場合には、該要素の制御特性の良否
を判定しないことを特徴とする制御系の診断方
法。 2 制御系の制御偏差の2乗積分値を演算する演
算手段と、 前記2乗積分値を、予め設定された基準値と比
較する第1の比較手段と、 制御系の各要素について、該要素の入力の出力
に対する寄与度を表わすコヒーレンス関数及び該
要素の入力に対する出力の関係を表わす伝達関数
を求める周波数解析手段と、 前記コヒーレンス関数の値を、予め設定された
基準値と比較する第2の比較手段と、 前記伝達関数を、予め設定されたそれぞれの基
準値と比較して、ゲイン余裕及び位相余裕を求め
る第3の比較手段と、 前記各比較手段の出力によつて、前記2乗積分
値がその基準値より大である時は、制御ループと
しての制御性不良と判定し、前記コヒーレンス関
数の値がその基準値より大である場合に、前記伝
達関数のゲイン余裕と位相余裕の値が、予め設定
されたそれぞれの基準範囲外である時は、該要素
の制御特性不良と判定し、前記コヒーレンス関数
の値が、前記基準値以下である時は、該要素の制
御特性の良否を判定しない論理回路と、 を備えたことを特徴とする制御系の診断装置。
[Claims] 1. When diagnosing the normality of the control system, calculate the square integral value of the control deviation of the control system, and when the square integral value is larger than a preset reference value, It is determined that the controllability of the control loop is poor, and for each element of the control system, a coherence function representing the contribution of the input of the element to the output is calculated, and if the value of the coherence function is greater than a preset reference value. If so, the relationship between the input and output of the element is determined as a transfer function, and when the gain margin and phase margin of the transfer function are outside the respective preset reference ranges, it is determined that the control characteristics of the element are defective. A method for diagnosing a control system, characterized in that, on the other hand, when the value of the coherence function is less than or equal to the reference value, the control characteristics of the element are not determined. 2. Calculating means for calculating the square integral value of the control deviation of the control system; first comparing means for comparing the square integral value with a preset reference value; and for each element of the control system, the element frequency analysis means for determining a coherence function representing the contribution of the input to the output of the element and a transfer function representing the relationship between the input and the output of the element; and a second frequency analysis means for comparing the value of the coherence function with a preset reference value. a comparison means; a third comparison means for determining a gain margin and a phase margin by comparing the transfer function with respective preset reference values; When the value is larger than the reference value, it is determined that the controllability as a control loop is poor, and when the value of the coherence function is larger than the reference value, the values of the gain margin and phase margin of the transfer function are determined. is outside the respective preset reference ranges, it is determined that the control characteristics of the element are poor, and when the value of the coherence function is below the reference value, the control characteristics of the element are determined to be good or bad. A control system diagnostic device comprising: a logic circuit that does not perform judgment;
JP57180192A 1982-10-14 1982-10-14 Method and apparatus for diagnosing control system Granted JPS5969813A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP57180192A JPS5969813A (en) 1982-10-14 1982-10-14 Method and apparatus for diagnosing control system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP57180192A JPS5969813A (en) 1982-10-14 1982-10-14 Method and apparatus for diagnosing control system

Publications (2)

Publication Number Publication Date
JPS5969813A JPS5969813A (en) 1984-04-20
JPH0425563B2 true JPH0425563B2 (en) 1992-05-01

Family

ID=16078999

Family Applications (1)

Application Number Title Priority Date Filing Date
JP57180192A Granted JPS5969813A (en) 1982-10-14 1982-10-14 Method and apparatus for diagnosing control system

Country Status (1)

Country Link
JP (1) JPS5969813A (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2628757B1 (en) * 1988-03-17 1992-01-17 Atochem METHOD FOR REGULATING AN ELECTROLYSIS CELL, ITS APPLICATION TO THE PRODUCTION OF CHLORINE AND SODA BY ELECTROLYSIS OF SODIUM CHLORIDE IN SOLUTION
JP3431390B2 (en) * 1996-03-18 2003-07-28 三菱電機株式会社 Cable insulation diagnosis method
US6804618B2 (en) * 1997-09-29 2004-10-12 Fisher Controls International, Llc Detection and discrimination of instabilities in process control loops
JP2016157313A (en) * 2015-02-25 2016-09-01 東芝三菱電機産業システム株式会社 Trend monitoring system of steel plant

Also Published As

Publication number Publication date
JPS5969813A (en) 1984-04-20

Similar Documents

Publication Publication Date Title
JPS6348249Y2 (en)
CN115950609B (en) Bridge deflection anomaly detection method combining correlation analysis and neural network
CN109990712A (en) A kind of width gage on-line calibration method
JPH0425563B2 (en)
CN115855165B (en) Multi-dimensional precision point inspection method and system for thermal power equipment
CN117630800A (en) Fault diagnosis method and system for automatic calibrating device of electric energy meter
KR101946631B1 (en) System and method for monitoring state of structure based compressive sensing
CN114050386B (en) Debugging method based on cavity filter scattering parameter and screw mapping extraction
JP3035089B2 (en) Process characteristic estimation method and apparatus, process monitoring method and process control method using the estimation method
US11846539B2 (en) Method and system for measuring interference in a checkweighing device
CN111366281B (en) XRD method residual stress detection accuracy testing and shear stress calculation method
JP3201882B2 (en) Error Correction Method for Rockwell Testing Machine
RU2222865C1 (en) Computer-aided diagnostic complex
CN111854929B (en) Method and system for evaluating uncertainty of sound power level of extra-high voltage main equipment
Yadav Measurement systems analysis and a study of anova method
JPH07253370A (en) Measuring system of six component forces
JP2002131439A (en) Deterioration diagnosing device for semiconductor detector
JPH04169849A (en) Threshold setting method for discriminating data
JPH0731076B2 (en) Abnormality diagnosis method for rotating machinery
JPS6358508A (en) Method and device for diagnosis of plant abnormality
JPS58148971A (en) Decision device for insulation and conduction of signal line of substrate
JP3071875B2 (en) IC test equipment
KR0118990B1 (en) Method and device for the control of sheet thickness
JPH06309580A (en) Monitor and diagnosis device for plant
CN113671435A (en) Method and device for evaluating step response characteristics of direct current transformer