JPH09120303A - プロセス同定方法およびその装置 - Google Patents
プロセス同定方法およびその装置Info
- Publication number
- JPH09120303A JPH09120303A JP29932195A JP29932195A JPH09120303A JP H09120303 A JPH09120303 A JP H09120303A JP 29932195 A JP29932195 A JP 29932195A JP 29932195 A JP29932195 A JP 29932195A JP H09120303 A JPH09120303 A JP H09120303A
- Authority
- JP
- Japan
- Prior art keywords
- impulse response
- wavelet transform
- response function
- equation
- estimated
- 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 68
- 238000005316 response function Methods 0.000 claims abstract description 16
- 238000006243 chemical reaction Methods 0.000 claims abstract description 3
- 238000012067 mathematical method Methods 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims 1
- 230000009466 transformation Effects 0.000 abstract 2
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 abstract 1
- 238000004886 process control Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 238000001914 filtration Methods 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 238000012369 In process control Methods 0.000 description 1
- 238000010965 in-process control Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Feedback Control In General (AREA)
Abstract
や統計的な方法の複雑さや勘と経験への依存をなくし、
精度の高いプロセス同定方法を提供する。 【解決手段】 プロセスの入出力信号をセンサーを用い
て検出し、計測された信号をA−D変換によりデジタル
信号に変換し、さらにインパルス応答関数のウェーブレ
ット変換を数学的な手法を用いて推定し、推定したイン
パルス応答関数のウェーブレット変換をウェーブレット
逆変換することによりインパルス応答関数を推定して、
プロセスを同定する。
Description
けるプロセス同定方法に関するものである。
は、物理的、化学的な法則から数式モデルを導く方法や
自己回帰モデルに代表される統計的な方法が用いられて
いた。
的、化学的な法則から数式モデルを導く方法において
は、プロセスが複雑すぎて定式化が困難であったり、仮
に定式化できたとしても式中の物理量が計測できず実用
的でない場合もあった。
な方法においては、採取したデータの質がモデルの精度
に大きく影響する。このため、通常採取したデータから
解析に必要なプロセスノイズと有害な観測ノイズを分離
するために、フィルタリング等の前処理を行うこと多
い。しかし、フィルタリング等の前処理は、解析者の勘
と経験に依ることが多く、解析者によって同定結果が大
きく異なる場合があった。
であり、解析者の勘と経験に依ることなく、精度の高い
プロセス同定方法を提供することにある。
出力信号を計測し、計測された信号からインパルス応答
のウェーブレット変換を数学的な手法を用いて推定し、
推定したインパルス応答のウェーブレット変換をウェー
ブレット逆変換することでインパルス応答を推定するこ
とを特徴とするプロセス同定方法である。
スに入力される入力信号と出力される出力信号をデジタ
ル信号に変換するA/D変換手段と、デジタル変換した
入出力信号に基づいてインパルス応答関数のウェーブレ
ット変換を推定するインパルス応答のウェーブレット変
換推定手段と、推定したインパルス応答関数のウェーブ
レット変換をウェーブレット逆変換してインパルス応答
関数を推定する手段を備えたことを特徴とするプロセス
同定装置である。
する。プロセス同定方法として、物理的、化学的な法則
から数式モデルを導く方法は、プロセスが複雑すぎて定
式化が困難であったり、仮に定式化できたとしても式中
の物理量が計測できず実用的でない場合もあるという問
題点があった。
な方法は、採取したデータの質がモデルの精度に大きく
影響する。このため、通常採取したデータから解析に必
要なプロセスノイズと有害な観測ノイズを分離するため
に、フィルタリング等の前処理を行うこと多い。しか
し、フィルタリング等の前処理は、解析者の勘と経験に
依ることが多く、解析者によって同定結果が大きく異な
る場合があるという問題点があった。
法として近年注目されているウェーブレット変換に着眼
した。ウェーブレット変換は、次式で定義される。
ーリングパラメータ、bはシフトパラメータ、h(t)
は基本ウェーブレット、*は複素共役を示す。スケーリ
ングパラメータaは、周波数特性を変化させる働きを持
ち、シフトパラメータbは、時間軸を変化させる働きを
持つ。
トと呼ばれる関数を時間軸に対してスケーリングした
り、シフトしたものを積分核とした積分変換である。こ
のウェーブレット変換を用いることで信号を時間−周波
数領域に展開することができ、従来、時間領域や周波数
領域で比較的困難とされていた信号処理が可能となる。
そこで、ウェーブレット変換をプロセス同定に適用する
ことにより、解析に必要なプロセスノイズと有害な観測
ノイズを分離することができる。このことから、ウェー
ブレット変換を用いてプロセスを同定する有益な方法を
見出した。またその方法を具体的に実施する装置を提供
する。
図である。処理の流れとしては、まず、プロセスの入出
力に対応する信号をセンサを用いて検出し、計測された
信号(アナログ信号)をAーD変換器によってデジタル
信号に変換する。さらに、インパルス応答のウェーブレ
ット変換の最小二乗推定値wg(a,b)を次式によっ
て求める。
u(t)はプロセスの入力、h(t)は、基本ウェーブ
レット、aはスケーリングパラメータ、bはシフトパラ
メータを示す。
ット変換の推定値wg(a,b)を次式に代入すること
によって、インパルス応答関数g(t)を推定する。
むだ時間や時定数等が正確に把握でき、プロセスの動特
性にあわせた制御系設計、改造が可能となる。これによ
り、制御性の向上が図られ、安定運転、省エネルギーに
つながる。また、工程解析やプロセスシミュレーション
にも活用することができる。さらに、モデル予測制御装
置等に組み込むことにより、プロセスの動特性の変動に
あわせて自動的に追従するような制御システムを構築す
ることも可能である。
以下に説明する。シミュレーションでは、プロセスモデ
ルとして、むだ時間+一次遅れを仮定した。プロセスモ
デルを以下に示す。
た。図2は、シミュレーションに用いた入力データとそ
の時間応答をプロットしたものである。ここで横軸の時
間は、データのサンプリング間隔に対応した時間を表す
ものである。(以下も同様とする。)本実施例では、基
本ウェーブレットとして次式を用いた。ここでは、ω=
2.0、γ=5.336とした。
ーブレット変換の最小二乗推定値wg(a,b)を示し
たものである。また、図4は、プロセスモデルのインパ
ルス応答(上段)と、推定したインパルス応答(下段)
を示したものである。この図から、本発明の方法によっ
て、インパルス応答が正確に推定されることがわかる。
ーブレット変換の最小二乗推定値であることを、以下に
説明する。
(t)、インパルス応答をg(t)とすると(6)式が
成り立つ。
ーブレット変換をwg(a,b)とするとウェーブレッ
ト逆変換により(4)式が成り立つ。ここで、h(t)
は、基本ウェーブレットであり、aはスケーリングパラ
メータ、bはシフトパラメータである。
と、
力u(t)にある種のフィルターを通したものと考える
ことができる。(9)式をa,bに関して適当に離散化
すると、
乗推定値は、残差の二乗和Qを最小とするような値であ
る。いま、
て0とおく。すなわち、
ロセスの入出力データを計測し、計測したデータからイ
ンパルス応答を精度良く推定することにより、制御性の
向上を図ることが容易に可能となる。これは、ウェーブ
レット変換をプロセス同定に適用することにより、解析
に必要なプロセスノイズと有害な観測ノイズを分離で
き、精度良く同定できたことによるものである。
ロック図である。
るためのもので、むだ時間を10時間、時定数を10時
間とした場合のプロセスの入力特性図、出力特性図であ
る。
るためのもので、プロセスのインパルス応答のウェーブ
レット変換を推定した結果を示す説明図である。
るためのもので、モデルのインパルス応答と推定したイ
ンパルス応答を示す特性図である。
換推定装置、5 インパルス応答同定装置。
Claims (2)
- 【請求項1】 プロセスの入出力信号を計測し、計測さ
れた信号からインパルス応答関数のウェーブレット変換
を数学的な手法を用いて推定し、推定したインパルス応
答関数のウェーブレット変換をウェーブレット逆変換す
ることでインパルス応答関数を推定することを特徴とす
るプロセス同定方法。 - 【請求項2】 プロセスに入力される入力信号と出力さ
れる出力信号をデジタル信号に変換するA/D変換手段
と、デジタル変換した入出力信号に基づいてインパルス
応答関数のウェーブレット変換を推定するインパルス応
答のウェーブレット変換推定手段と、推定したインパル
ス応答関数のウェーブレット変換をウェーブレット逆変
換してインパルス応答関数を推定する手段を備えたこと
を特徴とするプロセス同定装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP29932195A JP3502208B2 (ja) | 1995-10-25 | 1995-10-25 | プロセス同定方法およびその装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP29932195A JP3502208B2 (ja) | 1995-10-25 | 1995-10-25 | プロセス同定方法およびその装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH09120303A true JPH09120303A (ja) | 1997-05-06 |
JP3502208B2 JP3502208B2 (ja) | 2004-03-02 |
Family
ID=17871034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP29932195A Expired - Fee Related JP3502208B2 (ja) | 1995-10-25 | 1995-10-25 | プロセス同定方法およびその装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3502208B2 (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999022280A1 (de) * | 1997-10-27 | 1999-05-06 | Siemens Aktiengesellschaft | Verfahren zur bestimmung von stationären zudständen eines technischen prozesses |
WO2000002103A1 (de) * | 1998-07-07 | 2000-01-13 | Siemens Aktiengesellschaft | Verfahren zum betreiben einer regelungseinrichtung und vorrichtung zur durchführung des verfahrens |
JP2006520038A (ja) * | 2003-02-14 | 2006-08-31 | ドレッサ、インク | オンライン弁診断を実行するための方法、システムおよび記憶媒体 |
WO2007139088A1 (ja) * | 2006-05-29 | 2007-12-06 | Kyoto University | 変数決定方法、変数決定装置、プログラム及び記録媒体 |
JP2009104659A (ja) * | 2001-12-18 | 2009-05-14 | Mts Systems Corp | 制御システムのための制御パラメータを確定する方法 |
US7890216B2 (en) | 2007-01-05 | 2011-02-15 | Dresser, Inc. | Control valve and positioner diagnostics |
-
1995
- 1995-10-25 JP JP29932195A patent/JP3502208B2/ja not_active Expired - Fee Related
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999022280A1 (de) * | 1997-10-27 | 1999-05-06 | Siemens Aktiengesellschaft | Verfahren zur bestimmung von stationären zudständen eines technischen prozesses |
WO2000002103A1 (de) * | 1998-07-07 | 2000-01-13 | Siemens Aktiengesellschaft | Verfahren zum betreiben einer regelungseinrichtung und vorrichtung zur durchführung des verfahrens |
US6480750B2 (en) | 1998-07-07 | 2002-11-12 | Siemens Aktiengesellschaft | Controlling system and method for operating a controlling system |
JP2009104659A (ja) * | 2001-12-18 | 2009-05-14 | Mts Systems Corp | 制御システムのための制御パラメータを確定する方法 |
JP2006520038A (ja) * | 2003-02-14 | 2006-08-31 | ドレッサ、インク | オンライン弁診断を実行するための方法、システムおよび記憶媒体 |
WO2007139088A1 (ja) * | 2006-05-29 | 2007-12-06 | Kyoto University | 変数決定方法、変数決定装置、プログラム及び記録媒体 |
US7890216B2 (en) | 2007-01-05 | 2011-02-15 | Dresser, Inc. | Control valve and positioner diagnostics |
Also Published As
Publication number | Publication date |
---|---|
JP3502208B2 (ja) | 2004-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107329932B (zh) | 基于非线性调频分量分解的时频域模态参数辨识方法 | |
CN107607065A (zh) | 一种基于变分模态分解的冲击回波信号分析方法 | |
CN1441906A (zh) | 波形检测装置和使用该装置的状态监视系统 | |
Ferrero et al. | Modeling and processing measurement uncertainty within the theory of evidence: Mathematics of random–fuzzy variables | |
CN109798920A (zh) | 基于改进emd的mems陀螺随机误差建模滤波方法 | |
JP3502208B2 (ja) | プロセス同定方法およびその装置 | |
Gillberg | Methods for frequency domain estimation of continuous-time models | |
Xianmin | A new method with high confidence for validation of computer simulation models of flight systems | |
Malek | Blind compensation of memoryless nonlinear distortions in sparse signals | |
CN113010038B (zh) | 一种基于超分辨率重构的超声兰姆波触控载荷识别方法 | |
Schimmack et al. | Scaling-based least squares methods with implemented Kalman filter approach for nano-parameters identification | |
Ulapane et al. | System Identification of Static Nonlinear Elements: A Unified Approach of Active Learning, Over-fit Avoidance, and Model Structure Determination | |
Vergara et al. | On estimating the center frequency of ultrasonic pulses | |
JPH0677211B2 (ja) | プラント・モデリング装置 | |
Lahalle et al. | Continuous ARMA spectral estimation from irregularly sampled observations | |
Shenshu et al. | Adaptive filtering of color noise using the Kalman filter algorithm | |
JP2843991B2 (ja) | 音紋パターンの変換方法 | |
Asadi et al. | Optimal number of neurons for a two layer neural network model of a process | |
CN113297908B (zh) | 一种脉冲激励下基于模式搜索的非线性刚度识别方法 | |
He et al. | Delay Estimation of Dynamic System Based on Correlation Coefficient | |
CN110243562B (zh) | 一种适用于工程结构振动监测信号的在线递推谱分析方法 | |
Rafajłowicz | Time-domain optimization of input signals for distributed-parameter systems identification | |
Robertsson et al. | Analyzing bacteriological growth using wavelet transform | |
JP4156177B2 (ja) | ピエゾスキャナの歪み補正方法 | |
JPS60135772A (ja) | 信号分析装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Effective date: 20031202 Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20031204 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071212 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 5 Free format text: PAYMENT UNTIL: 20081212 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 5 Free format text: PAYMENT UNTIL: 20081212 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091212 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091212 Year of fee payment: 6 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091212 Year of fee payment: 6 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 7 Free format text: PAYMENT UNTIL: 20101212 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 8 Free format text: PAYMENT UNTIL: 20111212 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 9 Free format text: PAYMENT UNTIL: 20121212 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 9 Free format text: PAYMENT UNTIL: 20121212 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Year of fee payment: 10 Free format text: PAYMENT UNTIL: 20131212 |
|
LAPS | Cancellation because of no payment of annual fees |