JP2021135860A - Frequency characteristic estimation device and frequency characteristic estimation method - Google Patents
Frequency characteristic estimation device and frequency characteristic estimation method Download PDFInfo
- Publication number
- JP2021135860A JP2021135860A JP2020032761A JP2020032761A JP2021135860A JP 2021135860 A JP2021135860 A JP 2021135860A JP 2020032761 A JP2020032761 A JP 2020032761A JP 2020032761 A JP2020032761 A JP 2020032761A JP 2021135860 A JP2021135860 A JP 2021135860A
- Authority
- JP
- Japan
- Prior art keywords
- input
- output
- signal
- value generation
- control target
- 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
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000004044 response Effects 0.000 claims abstract description 45
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000010586 diagram Methods 0.000 description 26
- 230000005540 biological transmission Effects 0.000 description 15
- 238000012546 transfer Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000005094 computer simulation Methods 0.000 description 2
- 239000013078 crystal Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000004886 process control Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
Description
本発明は、周波数特性推定装置及び周波数特性推定方法に関する。 The present invention relates to a frequency characteristic estimation device and a frequency characteristic estimation method.
制御器の構造パラメータを設定し、性能のよい制御系を設計するためには、制御対象の応答特性を把握することが求められる。特に、安定性、応答特性を向上させるためには、制御対象の周波数特性を把握することが求められる。制御対象の周波数特性は、種々の周波数の信号を制御対象に入力して、その出力信号と入力信号とを計測して推定される。また、制御対象の周波数特性を理解し易くするため、推定された周波数特性からボード線図が作成される(例えば、特許文献1)。 In order to set the structural parameters of the controller and design a control system with good performance, it is necessary to understand the response characteristics of the controlled object. In particular, in order to improve stability and response characteristics, it is necessary to understand the frequency characteristics of the controlled object. The frequency characteristics of the controlled object are estimated by inputting signals of various frequencies to the controlled object and measuring the output signal and the input signal. Further, in order to make it easier to understand the frequency characteristics of the controlled object, a Bode diagram is created from the estimated frequency characteristics (for example, Patent Document 1).
特許文献1では、閉ループ回路に種々の周波数の正弦波信号を入力して、その出力信号と入力信号との関係から周波数特性を推定し、ボード線図を作成することとしている。
In
精度よく制御対象の周波数特性を推定するためには、種々の周波数の正弦波信号に対する入出力応答の取得実験を要するので、同定に要する時間、コストが増大する。 In order to accurately estimate the frequency characteristics of the controlled object, it is necessary to perform an input / output response acquisition experiment for sinusoidal signals of various frequencies, so that the time and cost required for identification increase.
本発明は、上述の事情に鑑みてなされたものであり、短時間、低コストで制御対象の周波数特性を推定することができる周波数特性推定装置及び周波数特性推定方法を提供することを目的とする。 The present invention has been made in view of the above circumstances, and an object of the present invention is to provide a frequency characteristic estimation device and a frequency characteristic estimation method capable of estimating the frequency characteristics of a controlled object in a short time and at low cost. ..
上記目的を達成するために、この発明の第1の観点に係る周波数特性推定装置は、
制御対象に所定の入力信号を入力して得られる、少なくとも一組の入出力応答を取得する開ループデータ取得部と、
前記所定の入力信号を入力して正弦波信号を出力する推定値生成フィルタと、前記制御対象とに基づいて、複数の異なる周波数の正弦波信号に対する前記制御対象の入出力信号を生成する演算部と、を備え、
前記演算部は、
前記推定値生成フィルタが出力する正弦波信号を所定の周波数に設定し、
前記開ループデータ取得部で取得された入出力応答の入力信号を、前記推定値生成フィルタに入力して得られる第1出力信号を生成し、
前記開ループデータ取得部で取得された前記制御対象の出力信号を、前記推定値生成フィルタに入力して得られる第2出力信号を生成し、
前記第1出力信号及び前記第2出力信号に基づいて前記制御対象の入出力信号を生成する。
In order to achieve the above object, the frequency characteristic estimation device according to the first aspect of the present invention is
An open-loop data acquisition unit that acquires at least one set of input / output responses obtained by inputting a predetermined input signal to the control target, and
An estimation value generation filter that inputs the predetermined input signal and outputs a sine wave signal, and a calculation unit that generates an input / output signal of the control target for a plurality of sine wave signals of different frequencies based on the control target. And with
The calculation unit
The sine wave signal output by the estimated value generation filter is set to a predetermined frequency, and the frequency is set.
The input signal of the input / output response acquired by the open loop data acquisition unit is input to the estimated value generation filter to generate a first output signal.
The output signal of the control target acquired by the open-loop data acquisition unit is input to the estimated value generation filter to generate a second output signal.
The input / output signal to be controlled is generated based on the first output signal and the second output signal.
また、前記所定の入力信号は、
ステップ入力又はインパルス入力である、
こととしてもよい。
In addition, the predetermined input signal is
Step input or impulse input,
It may be that.
この発明の第2の観点に係る周波数特性推定方法では、
制御対象に所定の入力信号を入力して得られる、少なくとも一組の入出力応答を取得する開ループデータ取得工程と、
前記所定の入力信号を入力して正弦波信号を出力する推定値生成フィルタと、前記制御対象とに基づいて、複数の異なる周波数の正弦波信号に対する前記制御対象の入出力信号を生成するデータ生成工程と、を含み、
前記データ生成工程では、
前記推定値生成フィルタが出力する正弦波信号を所定の周波数に設定し、
前記開ループデータ取得工程で取得された入出力応答の入力信号を、前記推定値生成フィルタに入力して得られる第1出力信号を生成し、
前記開ループデータ取得工程で取得された前記制御対象の出力信号を、前記推定値生成フィルタに入力して得られる第2出力信号を生成し、
前記第1出力信号及び前記第2出力信号に基づいて前記制御対象の入出力信号を生成する。
In the frequency characteristic estimation method according to the second aspect of the present invention,
An open-loop data acquisition process that acquires at least one set of input / output responses obtained by inputting a predetermined input signal to the control target, and
Data generation that generates input / output signals of the control target for a plurality of different frequency sine wave signals based on the estimation value generation filter that inputs the predetermined input signal and outputs the sine wave signal and the control target. Including the process
In the data generation step,
The sine wave signal output by the estimated value generation filter is set to a predetermined frequency, and the frequency is set.
The input signal of the input / output response acquired in the open loop data acquisition step is input to the estimated value generation filter to generate a first output signal.
The output signal of the controlled object acquired in the open-loop data acquisition step is input to the estimated value generation filter to generate a second output signal.
The input / output signal to be controlled is generated based on the first output signal and the second output signal.
また、前記所定の入力信号は、
ステップ入力又はインパルス入力である、
こととしてもよい。
In addition, the predetermined input signal is
Step input or impulse input,
It may be that.
本発明の周波数特性推定装置及び周波数特性推定方法によれば、少なくとも一組の入出力応答に基づいて、複数の周波数応答を生成することができるので、短時間、低コストで制御対象の周波数特性を推定することが可能である。 According to the frequency characteristic estimation device and the frequency characteristic estimation method of the present invention, a plurality of frequency responses can be generated based on at least one set of input / output responses, so that the frequency characteristics to be controlled can be controlled in a short time and at low cost. Can be estimated.
以下、図を参照しつつ、本発明の実施の形態に係る周波数特性推定装置1について説明する。
Hereinafter, the frequency
周波数特性推定装置1は、例えばコンピュータ装置であり、図1のブロック図に示すように、制御部21、記憶部22、表示部23、入力部24を備える。
The frequency
制御部21は、CPU(Central Processing Unit)、ROM(Read Only Memory)、RAM(Random Access Memory)、水晶発振器等から構成されており、周波数特性推定装置1の動作を制御するとともに、外部から取得された制御対象52の入出力応答に基づいて、制御対象52の周波数特性を推定する。制御部21は、制御部21のROM、記憶部22等に記憶されている各種動作プログラム及びデータをRAMに読み込んでCPUを動作させることにより、図1に示す制御部21の各機能を実現させる。これにより、制御部21は、開ループデータ取得部211、演算部212、ボード線図作成部213として動作する。
The
開ループデータ取得部211は、図1に示すように、制御対象52を含む制御システム30、外部記憶装置31、ネットワークを介したデータベース32等と入出力インターフェースを介して選択的に接続されて、制御対象52の開ループデータを取得する。取得される開ループデータは、制御対象52に所定の入力信号を入力した場合の、入出力応答データである。制御対象52に入力される所定の入力信号は、ステップ入力、インパルス入力等である。
As shown in FIG. 1, the open-loop
演算部212は、開ループデータ取得部211で取得された開ループデータ、推定値生成フィルタ51に基づいて、複数の周波数の周波数応答を生成する。推定値生成フィルタ51は、開ループデータで用いられている入力信号を入力して、正弦波信号を出力する特性を有するフィルタである。推定値生成フィルタ51の詳細な構成については後述する。
The
ボード線図作成部213は、演算部212で生成された周波数応答に基づいて、制御対象52の周波数特性を表すボード線図を作成する。
The Bode
記憶部22は、ハードディスク、フラッシュメモリ等の不揮発性メモリであり、開ループデータの入力信号に対応する推定値生成フィルタ51、生成された周波数応答、ボード線図等を記憶する。
The
表示部23は、コンピュータ装置である周波数特性推定装置1に備えられた表示用デバイスであり、例えば液晶パネルである。表示部23は、開ループデータ取得部211で取得された開ループデータ、ボード線図作成部213で作成されたボード線図等を表示する。
The
入力部24は、開ループデータの取得指示、周波数応答を算出するための設定周波数等を入力するための入力デバイスである。入力部24は、周波数特性推定装置1に備えられたキーボード、マウス等である。
The
続いて、周波数特性推定装置1を用いた周波数特性推定方法について、図2のフローチャートを参照しつつ、具体的に説明する。本実施の形態では、開ループデータとして予め取得された制御対象52のステップ応答を、外部記憶装置31から取得して、制御対象52の周波特性を推定する場合を例として説明する。
Subsequently, the frequency characteristic estimation method using the frequency
図2に示すように、開ループデータ取得工程として、制御対象52のステップ応答に係る入出力データを取得する(ステップS11)。本ステップS11で取得される開ループデータは、図3に示すように、実際の制御対象52にステップ入力である入力U0(s)を与え、制御対象52からの出力Y0(s)を測定することによって得られる入出力データである。本実施の形態に係るステップ入力は単位ステップ入力である。
As shown in FIG. 2, as an open-loop data acquisition step, input / output data related to the step response of the
また、本実施の形態に係る開ループデータは、予め実験データとして取得されて外部記憶装置31に記憶されており、開ループデータ取得部211は、外部記憶装置31と接続されて開ループデータを読み込むことにより、開ループデータを取得する。
Further, the open loop data according to the present embodiment is acquired in advance as experimental data and stored in the
続いて、データ生成工程として、ステップS11で取得した開ループデータに基づいて、複数の周波数の正弦波信号に対する制御対象52の入出力データを生成する。
Subsequently, as a data generation step, input / output data of the
ここで、データ生成工程における周波数応答の生成方法について説明する。図4(A)に示すように、システム50は、周波数特性を推定する対象である制御対象52と、推定値生成フィルタ51とを含むこととする。推定値生成フィルタ51は、開ループデータの入力信号である所定の入力信号を入力して、正弦波信号を出力する伝達特性を有するフィルタとして定義される。
Here, a method of generating a frequency response in the data generation step will be described. As shown in FIG. 4A, the
本実施の形態に係る推定値生成フィルタ51は、以下の式に示すように、単位ステップ信号を入力して、正弦波信号を出力する伝達特性X(s)を有する。
X(s)=ωs/(s2+ω2) ・・・(1)
As shown in the following equation, the estimated
X (s) = ωs / (s 2 + ω 2 ) ・ ・ ・ (1)
図4(A)は、推定値生成フィルタ51の出力側に制御対象52を接続したシステム50のブロック線図であり、図4(B)は、推定値生成フィルタ51の入力側に制御対象52を接続したシステム50’のブロック線図である。上記の2つのブロック線図に係る全体の伝達特性は、等価であるので、同じ信号の入力に対して、同じ信号を出力する。
FIG. 4A is a block diagram of the
データ生成工程で、生成したい周波数応答は、図4(A)に示す制御対象52の入力U(s)と出力Y(s)との関係である。すなわち、入力U(s)として種々の周波数の正弦波信号を与え、これらに対する出力Y(s)を生成することにより、制御対象52の周波数応答を生成する。
The frequency response to be generated in the data generation step is the relationship between the input U (s) and the output Y (s) of the
具体的には、推定値生成フィルタ51の伝達特性は式(1)のとおりであるので、推定値生成フィルタ51の周波数特性パラメータωを調整することにより、単位ステップ入力である入力U0(s)から種々の周波数の正弦波信号を生成し、制御対象52へ入力することができる。
Specifically, since the transmission characteristics of the estimated
また、図4(A)のシステム50と図4(B)のシステム50’との伝達特性は等価であるので、図4(A)の推定値生成フィルタ51の伝達特性と図4(B)の推定値生成フィルタ51の伝達特性とが等価であれば、等価な入力U0(s)に対するシステム50、50’の出力Y(s)は等価となる。
Further, since the transmission characteristics of the
図4(B)のシステム50’における制御対象52の入出力関係は、ステップS11で取得された開ループデータと置き換えることができるので、周波数特性パラメータωを所定の周波数に設定した推定値生成フィルタ51へ、取得された開ループデータの出力Y0(s)を入力することにより、オフラインで出力Y(s)を算出することができる。
Since the input / output relationship of the
算出された出力Y(s)は、上述したように図4(A)のシステム50の出力Y(s)と等価であるので、推定値生成フィルタ51の出力として算出可能なU(s)と、算出されたY(s)との関係から制御対象52の入出力データ、すなわち周波数応答を得ることができる。以上のように、実験データとして取得された開ループデータと、推定値生成フィルタ51の特性とに基づいて、種々の周波数における制御対象52の周波数応答をオフラインの周波数特性推定装置1で生成し、制御対象52の周波数特性を推定することができる。
Since the calculated output Y (s) is equivalent to the output Y (s) of the
図2のフローチャートに戻り、データ生成工程では、まず、推定値生成フィルタ51の伝達特性を決定する(ステップS12)。推定値生成フィルタ51は、単位ステップ入力から正弦波信号を出力するものであり、式(1)の伝達特性を有する。
Returning to the flowchart of FIG. 2, in the data generation step, first, the transmission characteristic of the estimated
式(1)及び図4(A)に示されるように、推定値生成フィルタ51は、開ループデータにおける所定の入力信号(単位ステップ信号)をキャンセルし、正弦波信号を生成する伝達特性を有する。予め所定の入力信号が設定されていない場合、開ループデータ取得部211で取得された開ループデータの入力信号に対応する推定値生成フィルタ51として、記憶部22に記憶されている推定値生成フィルタ51を、演算部212が読み込んで使用することとすればよい。
As shown in the equation (1) and FIG. 4 (A), the estimated
例えば、開ループデータ取得部211で取得された開ループデータの入力信号が、単位インパルス信号であった場合、演算部212は、予め記憶部22に記憶されている、インパルス入力に対応する推定値生成フィルタ51を読み込んで用いる。この場合の推定値生成フィルタ51の伝達特性は、X(s)=ω/(s2+ω2)である。
For example, when the input signal of the open loop data acquired by the open loop
続いてステップS13として、推定値生成フィルタ51の周波数特性を設定する。具体的には、推定値生成フィルタ51から出力される正弦波の周波数を決定する式(1)の周波数特性パラメータωを所定の値、例えば0.01rad/secに設定する。
Subsequently, in step S13, the frequency characteristics of the estimated
周波数特性パラメータωを設定した推定値生成フィルタ51に、ステップS11で取得された開ループデータの入力U0(s)を入力して、第1出力信号^U(s)を生成する(ステップS14)。推定値生成フィルタ51の伝達特性は、ステップ入力から正弦波信号を出力するように設定されているので、第1出力信号^U(s)は正弦波信号となる。
The input U 0 (s) of the open loop data acquired in step S11 is input to the estimated
また、周波数特性パラメータωを設定した推定値生成フィルタ51に、ステップS11で取得された開ループデータの出力Y0(s)を入力して、第2出力信号^Y(s)を算出する(ステップS15)。
Further, the output Y 0 (s) of the open loop data acquired in step S11 is input to the estimated
図4(A)、(B)に示すように、ステップS14、S15で算出されたシステム50の第1出力信号^U(s)及びシステム50’の第2出力信号^Y(s)は、制御対象52の入力及び出力に相当する。演算部212は、制御対象52の入出力信号に相当する第1出力信号^U(s)と第2出力信号^Y(s)との関係から、周波数応答を生成して、制御対象52のゲイン、位相を算出し、記憶部22に記憶させる(ステップS16)。
As shown in FIGS. 4A and 4B, the first output signal ^ U (s) of the
以下、制御対象52の周波数応答の生成を繰り返し、制御対象52の周波数特性を把握するために必要となるデータの取得を行う。より具体的には、推定値生成フィルタ51の周波数特性パラメータωが予め設定された上限値、例えば100rad/secに達するまで(ステップS17のNO)、所定の刻み幅、例えばωが10n≦ω<10n+1[rad/s](n:整数)の場合、ωを10n[rad/s]増加させて(ステップS18)、ステップS14からステップS16を繰り返す。
Hereinafter, the generation of the frequency response of the
推定値生成フィルタ51の周波数特性パラメータωが予め設定された上限値、例えば100rad/secに達した場合(ステップS17のYES)、制御対象52の周波数特性を把握するために十分なデータが取得されたものと判断し、周波数応答の生成を終了し、次の処理に進む。
When the frequency characteristic parameter ω of the estimated
ボード線図作成部213は、各周波数で算出されたゲイン、位相を記憶部22から読み出し、ボード線図を作成する(ステップS19)。また、ボード線図作成部213は、作成したボード線図を表示部23に表示させるととともに、記憶部22に記憶させ、周波数特性推定処理は終了する。
The Bode
以上説明したように、本実施の形態に係る周波数特性推定装置及び周波数特性推定方法によれば、一組の入出力応答である開ループデータに基づいて、オフラインの周波数特性推定装置1で種々の周波数における周波数応答を生成し、制御対象52の周波数特性を推定することができるので、短時間、低コストで制御対象52の周波数特性を推定することが可能である。
As described above, according to the frequency characteristic estimation device and the frequency characteristic estimation method according to the present embodiment, the offline frequency
また、本実施の形態に係る周波数特性推定方法は、一組の入出力応答に基づいて、数値計算によって、種々の周波数応答を生成して周波数特性を推定するので、制御対象52の伝達特性、すなわち伝達関数の次数等が不明な場合でも、容易にボード線図を作成し、周波数特性を推定することが可能である。
Further, in the frequency characteristic estimation method according to the present embodiment, various frequency responses are generated by numerical calculation based on a set of input / output responses to estimate the frequency characteristics. Therefore, the transfer characteristics of the
本実施の形態では、開ループデータは、一組の入出力応答であることとしたが、これに限られず、二組以上の入出力応答に基づいて周波数特性を推定することとしてもよい。この場合、周波数特性推定装置1は、各組の入出力応答に対して制御対象52の周波数特性を推定し、これらを並べて表示部23に表示させてもよい。また、演算部212は、推定された周波数特性の平均値を、制御対象52の周波数特性として推定することとしてもよい。
In the present embodiment, the open-loop data is a set of input / output responses, but the present invention is not limited to this, and the frequency characteristics may be estimated based on two or more sets of input / output responses. In this case, the frequency
また、本実施の形態に係る推定値生成フィルタ51の伝達特性は、X(s)=ωs/(s2+ω2)としたが、これに限られず、周波数応答を生成可能な特性であればよい。例えば、余弦関数を表す伝達特性であるX(s)=s2/(s2+ω2)としてもよい。
Further, the transmission characteristic of the estimated
以下、本実施の形態に係る周波数特性推定方法を用いて推定を行った場合の数値例について説明する。本例では、制御対象52の伝達特性G(s)を、以下の伝達関数で示される伝達特性として、コンピュータシミュレーションによる推定を行った。
Hereinafter, a numerical example in the case of estimation using the frequency characteristic estimation method according to the present embodiment will be described. In this example, the transfer characteristic G (s) of the
制御対象52に単位ステップ入力U0(s)を入力した場合の開ループデータである入力U0(t)、出力Y0(t)は、図5に示すとおりである。また、本発明の方法により、推定値生成フィルタ51の周波数特性パラメータωを10−2rad/secから102rad/secまで変化させて、生成された周波数応答に基づいて作成したボード線図を図6に示す。図6では、比較のため、コンピュータ上で作成した上記の制御対象52のボード線図も示している。
The input U 0 (t) and the output Y 0 (t), which are open-loop data when the unit step input U 0 (s) is input to the
図6に示すように、ゲイン特性は30rad/sec以下の範囲、位相特性は20rad/sec以下の範囲でコンピュータシミュレーションの結果とよく一致(ゲイン特性±10dB以下、位相特性±10deg以下)しており、良好な推定結果を得られていることがわかる。 As shown in FIG. 6, the gain characteristic is in the range of 30 rad / sec or less and the phase characteristic is in the range of 20 rad / sec or less, which is in good agreement with the result of the computer simulation (gain characteristic ± 10 dB or less, phase characteristic ± 10 deg or less). It can be seen that good estimation results have been obtained.
本発明によれば、上記の数値例のように、同定しにくいシステム、すなわち、非最小位相系で、2つの共振ピークを持つ振動系であり、相対次数が3であるシステムであっても、容易に周波数特性を推定することが可能である。 According to the present invention, even in a system that is difficult to identify, that is, a non-minimum phase system, an oscillating system having two resonance peaks, and a relative order of 3, as in the above numerical example. It is possible to easily estimate the frequency characteristics.
本発明は、制御対象の構造が不明である制御系の周波数特性推定に好適である。特に、入出力応答の取得に時間を要し、操業データの取得が難しいプロセス系の制御システムに好適である。 The present invention is suitable for estimating the frequency characteristics of a control system whose structure to be controlled is unknown. In particular, it is suitable for process control systems in which it takes time to acquire input / output responses and it is difficult to acquire operation data.
1 周波数特性推定装置、21 制御部、211 開ループデータ取得部、212 演算部、213 ボード線図作成部、22 記憶部、23 表示部、24 入力部、30 制御システム、31 外部記憶装置、32 データベース、50,50’ システム、51 推定値生成フィルタ、52 制御対象 1 Frequency characteristic estimation device, 21 control unit, 211 open loop data acquisition unit, 212 calculation unit, 213 board diagram creation unit, 22 storage unit, 23 display unit, 24 input unit, 30 control system, 31 external storage device, 32 Database, 50, 50'system, 51 estimate generation filter, 52 controlled object
Claims (4)
前記所定の入力信号を入力して正弦波信号を出力する推定値生成フィルタと、前記制御対象とに基づいて、複数の異なる周波数の正弦波信号に対する前記制御対象の入出力信号を生成する演算部と、を備え、
前記演算部は、
前記推定値生成フィルタが出力する正弦波信号を所定の周波数に設定し、
前記開ループデータ取得部で取得された入出力応答の入力信号を、前記推定値生成フィルタに入力して得られる第1出力信号を生成し、
前記開ループデータ取得部で取得された前記制御対象の出力信号を、前記推定値生成フィルタに入力して得られる第2出力信号を生成し、
前記第1出力信号及び前記第2出力信号に基づいて前記制御対象の入出力信号を生成する、
周波数特性推定装置。 An open-loop data acquisition unit that acquires at least one set of input / output responses obtained by inputting a predetermined input signal to the control target, and
An estimation value generation filter that inputs the predetermined input signal and outputs a sine wave signal, and a calculation unit that generates an input / output signal of the control target for a plurality of sine wave signals of different frequencies based on the control target. And with
The calculation unit
The sine wave signal output by the estimated value generation filter is set to a predetermined frequency, and the frequency is set.
The input signal of the input / output response acquired by the open loop data acquisition unit is input to the estimated value generation filter to generate a first output signal.
The output signal of the control target acquired by the open-loop data acquisition unit is input to the estimated value generation filter to generate a second output signal.
An input / output signal to be controlled is generated based on the first output signal and the second output signal.
Frequency characteristic estimator.
ステップ入力又はインパルス入力である、
請求項1に記載の周波数特性推定装置。 The predetermined input signal is
Step input or impulse input,
The frequency characteristic estimation device according to claim 1.
前記所定の入力信号を入力して正弦波信号を出力する推定値生成フィルタと、前記制御対象とに基づいて、複数の異なる周波数の正弦波信号に対する前記制御対象の入出力信号を生成するデータ生成工程と、を含み、
前記データ生成工程では、
前記推定値生成フィルタが出力する正弦波信号を所定の周波数に設定し、
前記開ループデータ取得工程で取得された入出力応答の入力信号を、前記推定値生成フィルタに入力して得られる第1出力信号を生成し、
前記開ループデータ取得工程で取得された前記制御対象の出力信号を、前記推定値生成フィルタに入力して得られる第2出力信号を生成し、
前記第1出力信号及び前記第2出力信号に基づいて前記制御対象の入出力信号を生成する、
周波数特性推定方法。 An open-loop data acquisition process that acquires at least one set of input / output responses obtained by inputting a predetermined input signal to the control target, and
Data generation that generates input / output signals of the control target for a plurality of different frequency sine wave signals based on the estimation value generation filter that inputs the predetermined input signal and outputs the sine wave signal and the control target. Including the process
In the data generation step,
The sine wave signal output by the estimated value generation filter is set to a predetermined frequency, and the frequency is set.
The input signal of the input / output response acquired in the open loop data acquisition step is input to the estimated value generation filter to generate a first output signal.
The output signal of the controlled object acquired in the open-loop data acquisition step is input to the estimated value generation filter to generate a second output signal.
An input / output signal to be controlled is generated based on the first output signal and the second output signal.
Frequency characteristic estimation method.
ステップ入力又はインパルス入力である、
請求項3に記載の周波数特性推定方法。 The predetermined input signal is
Step input or impulse input,
The frequency characteristic estimation method according to claim 3.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020032761A JP7294662B2 (en) | 2020-02-28 | 2020-02-28 | Frequency characteristic estimation device and frequency characteristic estimation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020032761A JP7294662B2 (en) | 2020-02-28 | 2020-02-28 | Frequency characteristic estimation device and frequency characteristic estimation method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021135860A true JP2021135860A (en) | 2021-09-13 |
JP7294662B2 JP7294662B2 (en) | 2023-06-20 |
Family
ID=77662377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020032761A Active JP7294662B2 (en) | 2020-02-28 | 2020-02-28 | Frequency characteristic estimation device and frequency characteristic estimation method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7294662B2 (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005050058A (en) * | 2003-07-31 | 2005-02-24 | Fuji Electric Holdings Co Ltd | Transfer function estimation device for characteristic-unknown system |
JP2006010596A (en) * | 2004-06-29 | 2006-01-12 | Fuji Electric Systems Co Ltd | Electronic transformer with oscillation determination function |
-
2020
- 2020-02-28 JP JP2020032761A patent/JP7294662B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005050058A (en) * | 2003-07-31 | 2005-02-24 | Fuji Electric Holdings Co Ltd | Transfer function estimation device for characteristic-unknown system |
JP2006010596A (en) * | 2004-06-29 | 2006-01-12 | Fuji Electric Systems Co Ltd | Electronic transformer with oscillation determination function |
Also Published As
Publication number | Publication date |
---|---|
JP7294662B2 (en) | 2023-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schoukens et al. | Identification of Wiener–Hammerstein systems by a nonparametric separation of the best linear approximation | |
KR101021895B1 (en) | Method and system for processing a sound field representation | |
Liao et al. | Mesh-robustness of an energy stable BDF2 scheme with variable steps for the Cahn–Hilliard model | |
JP7040189B2 (en) | Electromagnetic wave measurement point calculation device and radiation interference wave measurement device | |
Altoè et al. | Transmission line cochlear models: improved accuracy and efficiency | |
Jin et al. | A new non-parametric stationarity test of time series in the time domain | |
JP2009175917A (en) | Controlled parameter adjustment method and controlled parameter adjustment program | |
CN112560974B (en) | Information fusion and vehicle information acquisition method and device | |
TWI774919B (en) | Information processing device, program, process execution device and information processing system | |
JP2017017782A (en) | Servo controller with function for optimizing control gain in accordance with evaluation function in online manner | |
KR20180136492A (en) | Motor control device | |
Worden et al. | On Gaussian process NARX models and their higher-order frequency response functions | |
Varona et al. | Automatic adaptive sampling in parametric model order reduction by matrix interpolation | |
Giordano et al. | Consistency aspects of Wiener-Hammerstein model identification in presence of process noise | |
JP2021135860A (en) | Frequency characteristic estimation device and frequency characteristic estimation method | |
Lee et al. | A Novel PID Controller Tuning for Desired Time Response Using a Pulse Response | |
Meher et al. | Relay‐based identification of Wiener model | |
Martín et al. | An enhanced identification test monitoring procedure for MIMO systems relying on uncertainty estimates | |
Kraszewski et al. | Nonlinear Kalman filtering in the presence of additive noise | |
Yuan et al. | Sequential B-spline surface construction using multiresolution data clouds | |
Ježek et al. | New algorithm for spectral factorization and its practical application | |
Cohen Tenoudji et al. | Optimal filtering: Wiener and Kalman filters | |
Khakimzyanov et al. | Analyzing Nonlinear Behavior Sequences Through ASM | |
Isaev et al. | Analytical Representation of the Complex Frequency Response of a Hydrophone | |
WO2022183468A1 (en) | Control method for turntable with six-degrees-of-freedom adjustment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220616 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20230530 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20230531 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230601 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7294662 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |