IP Force 特許公報掲載プロジェクト 2022.1.31 β版

ホーム > 特許ランキング > ＢＩＰＲＯＧＹ株式会社

このエントリーをはてなブックマークに追加

知財求人 - 知財ポータルサイト「IP Force」

▶ 日本ユニシス株式会社の特許一覧

(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024019830
(43)【公開日】2024-02-14
(54)【発明の名称】時系列データ解析システム
(51)【国際特許分類】
   G06F 17/10 20060101AFI20240206BHJP
   G05B 23/02 20060101ALI20240206BHJP
   G01D 3/032 20060101ALI20240206BHJP
   G01K 3/08 20060101ALI20240206BHJP
   G01D 1/10 20060101ALI20240206BHJP
【ＦＩ】
G06F17/10 V
G05B23/02 Z
G01D3/032
G01K3/08
G01D1/10
【審査請求】未請求
【請求項の数】16
【出願形態】ＯＬ
(21)【出願番号】P 2022122546
(22)【出願日】2022-08-01
(71)【出願人】
【識別番号】591030237
【氏名又は名称】ＢＩＰＲＯＧＹ株式会社
(74)【代理人】
【識別番号】100094569
【弁理士】
【氏名又は名称】田中 伸一郎
(74)【代理人】
【識別番号】100103610
【弁理士】
【氏名又は名称】▲吉▼田 和彦
(74)【代理人】
【識別番号】100109070
【弁理士】
【氏名又は名称】須田 洋之
(74)【代理人】
【識別番号】100067013
【弁理士】
【氏名又は名称】大塚 文昭
(74)【代理人】
【識別番号】100086771
【弁理士】
【氏名又は名称】西島 孝喜
(74)【代理人】
【識別番号】100109335
【弁理士】
【氏名又は名称】上杉 浩
(74)【代理人】
【識別番号】100120525
【弁理士】
【氏名又は名称】近藤 直樹
(74)【代理人】
【識別番号】100139712
【弁理士】
【氏名又は名称】那須 威夫
(72)【発明者】
【氏名】大塚 洋
(72)【発明者】
【氏名】吉本 昌平
【テーマコード（参考）】
2F075
3C223
5B056
【Ｆターム（参考）】
2F075AA04
2F075AA05
2F075EE14
2F075EE18
2F075FF02
3C223AA11
3C223BA01
3C223CC01
3C223DD01
3C223EB01
3C223FF04
3C223FF06
3C223FF08
3C223GG01
3C223HH03
3C223HH28
5B056BB03
5B056GG03
5B056HH07
(57)【要約】
【課題】ノイズが多く含まれるデータに対してもより有効にデータの解析を行うことを可能とするデータ解析システム及び時系列データ解析方法を提供する。
【解決手段】時系列データ解析システムは、バンドパスフィルタ処理部１６、移動平均処理部１８、微分処理部２０、積分階差処理２２、減算処理部２４、パルス生成処理部２６、閾値処理部２８、分割点探索処理部３０、範囲設定処理部３２を含んでいる。これらの各処理部において、時系列データに対して各種の処理をさまざまに適用し実行することによって、時系列データの中の屈折点、すなわちトレンドが変化する点を高い精度で見いだすことができる。範囲設定処理部３２は、例えば、ユーザが時間軸に直接範囲を設定する処理を実行する。さらに、分割点探索処理部によって、コンピュータに対する負荷の少ない処理により時系列データによくフィットする近似曲線を求める。
【選択図】図１
【特許請求の範囲】
【請求項1】
時系列データに対して時間微分を実行し、必要に応じて時間微分後の時系列データの波形図を表示する時間微分処理手段と、
処理対象となる時系列データに対し、前記時間微分処理部による時間微分処理を実行するときに時間範囲を時間軸上で設定する範囲設定処理部と、
前記範囲設定処理部によって設定された時間範囲において、処理対象となる時系列データがピークとなる位置を検出するピーク位置検出部と、
を備える、時系列データ解析システム。
【請求項2】
前記範囲設定処理部は、時系列データの時間軸上に、前記時間範囲の少なくとも一方の境界をユーザが直接設定するものである、請求項１に記載の時系列データ解析システム。
【請求項3】
前記範囲設定処理部は、ユーザが縦軸の特定の値を設定し、それに対応する時間軸上の値を前記時間範囲の少なくとも一方の境界とするものである、請求項１に記載の時系列データ解析システム。
【請求項4】
前記範囲設定処理部は、ピーク位置検出部によって検出された時系列データの縦軸のピーク値を検索し、当該ピーク値に対応する時間軸上の値を、前記時間範囲の一方の境界とするものである、請求項１に記載の時系列データ解析システム。
【請求項5】
さらに、ノイズ除去手段として、時系列データに対して移動平均処理を実行し、必要に応じて移動平均処理後の時系列データの波形図を表示する移動平均処理部、及び、時系列データに対して積分階差処理を実行し、必要に応じて積分階差処理後の時系列データの波形図を表示する積分階差処理部のいずれか一方又は両方を備える、請求項１に記載の時系列データ解析システム。
【請求項6】
時系列データに対してバンドパスフィルタ処理を実行するバンドパスフィルタ処理部をさらに有する請求項１に記載の時系列データ解析システム。
【請求項7】
ユーザがユーザデータを入力するユーザ入力処理部をさらに有する請求項５に記載の時系列データ解析システム。
【請求項8】
外部から時系列データをシステム内に取り込むデータ入力部をさらに有する請求項５に記載の時系列データ解析システム。
【請求項9】
時系列データに対して時間微分を実行し、必要に応じて時間微分後の時系列データの波形図を表示する時間微分処理手段と、
時系列データに対して、その値が連続して正である時間の長さをパルス幅、連続して正である時間の長さをパルス高とする第１のパルス波形と、その値が連続して負である時間の長さをパルス幅、連続して負である時間の長さをパルス高とする第２のパルス波形とを求め、必要に応じて第１のパルス波形及び第２のパルス波形の両方又はいずれか一方の波形図を表示するパルス生成部と、
複数の時系列データの間又は複数のパルス波形の間で減算処理を実行し、必要に応じてその減算後の時系列データ又はパルス波形を表示する減算処理部と、
パルス波形の値を予め指定された閾値と比較し、当該閾値より大きいパルス又は小さいパルスを選択する閾値処理部と、
を備える時系列データ解析システム。
【請求項10】
さらに、ノイズ除去手段として、時系列データに対して移動平均処理を実行し、必要に応じて移動平均処理後の時系列データの波形図を表示する移動平均処理部、及び、時系列データに対して積分階差処理を実行し、必要に応じて積分階差処理後の時系列データの波形図を表示する積分階差処理部のいずれか一方又は両方を備える、請求項９に記載の時系列データ解析システム。
【請求項11】
時系列データに対してバンドパスフィルタ処理を実行するバンドパスフィルタ処理部をさらに有する請求項９に記載の時系列データ解析システム。
【請求項12】
ユーザがユーザデータを入力するユーザ入力処理部をさらに有する請求項９に記載の時系列データ解析システム。
【請求項13】
外部から時系列データをシステム内に取り込むデータ入力部をさらに有する請求項９に記載の時系列データ解析システム。
【請求項14】
時系列データ解析システムに取り込まれた時系列データの屈折点を求める方法であって、
前記時系列データ解析システムに含まれる微分処理部によって、取り込まれた第１の時系列データに対して３回の時間微分を行って第２の時系列データを得る第１ステップと、
前記時系列データ解析システムに含まれるパルス生成部によって、前記第２の時系列データが連続して正である時間の長さをパルス幅、連続して正である時間の長さをパルス高とする第１のパルス波形と、前記第４の時系列データが連続して負である時間の長さをパルス幅、連続して負である時間の長さをパルス高とする第２のパルス波形とを求める第２ステップと、
前記時系列データ解析システムに含まれる減算処理部によって、前記第１のパルス波形と前記第２のパルス波形との減算を行って第３のパルス波形を求める第３ステップと、
前記微分処理部によって、前記第３のパルス波形に対して時間微分を行って第４のパルス波形を求める第４ステップと、
前記時系列データ解析システムに含まれる閾値処理部によって、前記第４のパルス波形に対して閾値処理を行い、所定の閾値を超えるパルスのみを選択して、前記第１の時系列データの屈折点とする第５ステップと、
を含む、時系列データ解析方法。
【請求項15】
前記第１ステップで、２回目の時間微分を行って得られる時系列データに対して積分階差を行ってノイズを除去したあとに３回目の時間微分を実行して第２の時系列データとする、請求項１４に記載の時系列データ解析方法。
【請求項16】
時系列データの波形図の近似曲線を求めることによって時系列データを解析する時系列データ解析方法であって、
予め指定された始点と終点との間の前記時系列データ上に第１分割点から第ｎ分割点（ｎは２以上の整数）までを仮に配置するステップと、
第１分割点を始点から第２分割点まで移動させたときに始点と第１分割点とを結ぶ第１直線がその範囲で時系列データと最も良くフィットする点を第１分割点の位置と決定し、
第２分割点を第１分割点から第３分割点まで移動させたときに第１分割点と第２分割点とを結ぶ第２直線がその範囲で時系列データと最も良くフィットする点を第２分割点の位置と決定するステップと、
以下同様の手順を続行し、第ｎ分割点を第ｎ－１分割点から終点まで移動させたときに第ｎ－１分割点と第２分割点とを結ぶ第ｎ直線がその範囲で時系列データと最も良くフィットする点を第ｎ分割点の位置と決定するステップと、
前記始点と、第１分割点と、・・・第ｎ分割点と、終点とをこの順番で直線によって連結して得られる曲線を前記時系列データの近似曲線とするステップと、
を含む時系列データ解析方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、各種の測定装置や各種のセンサなどから出力される時系列データを解析するための、時系列データ解析システム及び時系列データ解析方法に関する。
【背景技術】
【0002】
製造業の生産現場では、さまざまな測定装置やセンサによって各種のデータが取得されており、製品の品質の維持・向上のためには、これらのデータを解析することが必要不可欠とされている。このようなデータ解析のための従来の方法としては、これまでのところ、ＭＴ法などの波形全体の特徴を捉える方法が主流であった。
【発明の概要】
【発明が解決しようとする課題】
【0003】
しかしながら、実際のデータには多数のノイズが含まれており、このようなノイズを多く含むデータに対しては、従来のＭＴ法などでは、波形に含まれる「特定の形」を演算で取り出すことが難しかった。
【0004】
そこで、本発明は、ノイズが多く含まれるデータに対してもより有効にデータの解析を行うことを可能とするデータ解析システム及び時系列データ解析方法を提供することを目的とする。
【課題を解決するための手段】
【0005】
上記の課題を解決するために、本発明に係る時系列データ解析システムは、
時系列データに対して時間微分を実行し、必要に応じて時間微分後の時系列データの波形図を表示する時間微分処理手段と、
処理対象となる時系列データに対し、前記時間微分処理部による時間微分処理を実行するときに時間範囲を時間軸上で設定する範囲設定処理部と、
前記範囲設定処理部によって設定された時間範囲において、処理対象となる時系列データがピークとなる位置を検出するピーク位置検出部と、
を備える（第１発明）。
【0006】
上記の第１発明において、前記範囲設定処理部は、時系列データの時間軸上に、前記時間範囲の少なくとも一方の境界をユーザが直接設定することができる（第２発明）。
【0007】
上記の第１発明において、前記範囲設定処理部は、ユーザが縦軸の特定の値を設定し、それに対応する時間軸上の値を前記時間範囲の少なくとも一方の境界とすることができる（第３発明）。
【0008】
上記の第１発明において、前記範囲設定処理部は、ピーク位置検出部によって検出された時系列データの縦軸のピーク値を検索し、当該ピーク値に対応する時間軸上の値を、前記時間範囲の一方の境界とすることができる（第４発明）。
【0009】
上記の第１発明において、さらに、ノイズ除去手段として、時系列データに対して移動平均処理を実行し、必要に応じて移動平均処理後の時系列データの波形図を表示する移動平均処理部、及び、時系列データに対して積分階差処理を実行し、必要に応じて積分階差処理後の時系列データの波形図を表示する積分階差処理部のいずれか一方又は両方を備えることができる（第５発明）。
【0010】
上記の第１発明において、時系列データに対してバンドパスフィルタ処理を実行するバンドパスフィルタ処理部をさらに有することができる（第６発明）。
【0011】
上記の第５発明において、ユーザがユーザデータを入力するユーザ入力処理部をさらに有することができる（第７発明）。
【0012】
上記の第５発明において、外部から時系列データをシステム内に取り込むデータ入力部をさらに有することができる（第８発明）。
【0013】
上記の課題を解決するために、本発明に係る時系列データ解析システムは、
時系列データに対して時間微分を実行し、必要に応じて時間微分後の時系列データの波形図を表示する時間微分処理手段と、
時系列データに対して、その値が連続して正である時間の長さをパルス幅、連続して正である時間の長さをパルス高とする第１のパルス波形と、その値が連続して負である時間の長さをパルス幅、連続して負である時間の長さをパルス高とする第２のパルス波形とを求め、必要に応じて第１のパルス波形及び第２のパルス波形の両方又はいずれか一方の波形図を表示するパルス生成部と、
複数の時系列データの間又は複数のパルス波形の間で減算処理を実行し、必要に応じてその減算後の時系列データ又はパルス波形を表示する減算処理部と、
パルス波形の値を予め指定された閾値と比較し、当該閾値より大きいパルス又は小さいパルスを選択する閾値処理部と、
を備える（第９発明）。
【0014】
上記の第９発明において、さらに、ノイズ除去手段として、時系列データに対して移動平均処理を実行し、必要に応じて移動平均処理後の時系列データの波形図を表示する移動平均処理部、及び、時系列データに対して積分階差処理を実行し、必要に応じて積分階差処理後の時系列データの波形図を表示する積分階差処理部のいずれか一方又は両方を備えることができる（第１０発明）。
【0015】
上記の第９発明において、時系列データに対してバンドパスフィルタ処理を実行するバンドパスフィルタ処理部をさらに有することができる（第１１発明）。
【0016】
上記の第９発明において、ユーザがユーザデータを入力するユーザ入力処理部をさらに有することができる（第１２発明）。
【0017】
上記の第９発明において、外部から時系列データをシステム内に取り込むデータ入力部をさらに有することができる（第１３発明）。
【0018】
上記の課題を解決するために、時系列データ解析システムに取り込まれた時系列データの屈折点を求める方法は、
前記時系列データ解析システムに含まれる微分処理部によって、取り込まれた第１の時系列データに対して３回の時間微分を行って第２の時系列データを得る第１ステップと、
前記時系列データ解析システムに含まれるパルス生成部によって、前記第２の時系列データが連続して正である時間の長さをパルス幅、連続して正である時間の長さをパルス高とする第１のパルス波形と、前記第４の時系列データが連続して負である時間の長さをパルス幅、連続して負である時間の長さをパルス高とする第２のパルス波形とを求める第２ステップと、
前記時系列データ解析システムに含まれる減算処理部によって、前記第１のパルス波形と前記第２のパルス波形との減算を行って第３のパルス波形を求める第３ステップと、
前記微分処理部によって、前記第３のパルス波形に対して時間微分を行って第４のパルス波形を求める第４ステップと、
前記時系列データ解析システムに含まれる閾値処理部によって、前記第４のパルス波形に対して閾値処理を行い、所定の閾値を超えるパルスのみを選択して、前記第１の時系列データの屈折点とする第５ステップと、
を含む（第１４発明）。
【0019】
上記の第１４発明において、前記第１ステップで、２回目の時間微分を行って得られる時系列データに対して積分階差を行ってノイズを除去したあとに３回目の時間微分を実行して第２の時系列データとするとすることができる（第１５発明）。
【0020】
上記の課題を解決するために、本発明における時系列データ解析方法は、
時系列データの波形図の近似曲線を求めることによって時系列データを解析する時系列データ解析方法であって、
予め指定された始点と終点との間の前記時系列データ上に第１分割点から第ｎ分割点（ｎは２以上の整数）までを仮に配置するステップと、
第１分割点を始点から第２分割点まで移動させたときに始点と第１分割点とを結ぶ第１直線がその範囲で時系列データと最も良くフィットする点を第１分割点の位置と決定し、
第２分割点を第１分割点から第３分割点まで移動させたときに第１分割点と第２分割点とを結ぶ第２直線がその範囲で時系列データと最も良くフィットする点を第２分割点の位置と決定するステップと、
以下同様の手順を続行し、第ｎ分割点を第ｎ－１分割点から終点まで移動させたときに第ｎ－１分割点と第２分割点とを結ぶ第ｎ直線がその範囲で時系列データと最も良くフィットする点を第ｎ分割点の位置と決定するステップと、
前記始点と、第１分割点と、・・・第ｎ分割点と、終点とをこの順番で直線によって連結して得られる曲線を前記時系列データの近似曲線とするステップと、
を含む（第１６発明）。
【図面の簡単な説明】
【0021】
【図1】本発明の実施の一形態に係る時系列データ解析システムの機能ブロックズである。
【図2】射出成形時に樹脂が流し込まれてから固化するまでの圧力の変化を測定する圧力センサの出力の時系列データである。
【図3】図２の波形図の一部を時間方向に引き延ばした波形図である。
【図4】時系列データに時間微分を行って得られる時系列データの一例を示す波形図である。
【図5】時系列データに対して移動平均をとって得られる時系列データの一例を示す波形図である。
【図6A】時系列データに対する積分階差について説明するための図である。
【図6B】時系列データに対する積分階差について説明するための図である。
【図7】ゼロを通過する時点を閾値処理で判定する場合の問題点を示す図である。
【図8】微分によって得られた波形図に基づいてパルスを生成するパルス生成部の動作を示す波形図である。
【図9】時系列データに対して各種処理を実行して、元の時系列データの屈折点を求める例を示す波形図である。
【図10】図９で求めた屈折点に対して閾値処理を実行して所定の閾値以上のパルスを選択することを示した波形図である。
【図11】ノイズが多く含まれる波形図から正しくピークの時点を求める方法を示す波形図である。
【図12】元の波形図の一つの屈折点を求める方法の一例を示した波形図である。
【図13】元の波形図の一つの屈折点を求める方法の一例を示した波形図である。
【図14】本発明の実施の一形態のシステムで実行できる処理ことをまとめた図である。
【図15】本発明の実施の一形態のシステムのユーザインターフェイスを示す図である。
【図16】時系列データ全体に対して階差処理を行った場合の問題点を説明する図である。
【図17】移動平均の度合いを強めることによってノイズを低減できるとこを説明する図である。
【図18】移動平均の点数を増やすことによる問題点を説明する図である。
【図19】オペレータが手動で時間範囲を設定する場合について説明する図である。
【図20】縦軸の値域に着目して範囲を決める方法について説明する図である。
【図21】特徴点の特性に基づいて検索する範囲を設定する方法について説明する図である。
【図22】速度データの最大点に基づいて範囲を設定する方法について説明する図である。
【図23】複数の時系列データのお互いの関係性に基づいて処理範囲を設定する処理について説明する図である。
【図24】実際に時系列データに対して処理範囲を設定する具体的な方法について説明する図である。
【図25】実際に時系列データに対して処理範囲を設定する具体的な方法について説明する図である。
【図26】実際に時系列データに対して処理範囲を設定する具体的な方法について説明する図である。
【図27】元の時系列データの近似曲線を求める方法を説明するための図である。
【図28】元の時系列データの近似曲線を求める方法を説明するための図である。
【図29】元の時系列データの近似曲線を求める方法を説明するための図である。
【発明を実施するための形態】
【0022】
図１は、本発明の実施の一形態に係る時系列データ解析システムの主要部の機能ブロック図である。図１に示す各処理部は、パーソナルコンピュータを含む一般のコンピュータにおけるプロセッサ、ＤＲＡＭ等の一時的記憶装置、ＨＤＤ装置等の非一時的記憶装置、データ入力手段としてのキーボードやマウス、ディスプレイ装置等のハードウェアと、記憶装置から読み出されてプロセッサによって実行されるコンピュータプログラム等のソフトウェアとの組み合わせによって、各種処理が行われる。
【0023】
図１において、データ入力処理部１０は、測定装置や各種のセンサなどからの時系列データをシステム内に取り込む処理を行う。ここで、時系列データとは、一定の時間間隔で取り込まれたデータ、あるいは既にシステムに取り込まれた時系列データに対して何らかの処理が行われたデータを指す。ユーザ入力処理部１２は、ユーザがキーボードやマウス等から入力したデータをシステム内に取り込む処理を行う。画像出力処理部１４は、文字、図形等を含むユーザインターフェイスや、実際の時系列の波形データ等をディスプレイや、必要に応じてプリンタ等の出力装置に出力する処理を行う。
【0024】
バンドパスフィルタ処理部１６は、システムに取り込まれた時系列データに対して、あるいは既にシステムに取り込まれて何らかの処理がなされた時系列データに対して、バンドパスフィルタ処理を実行し、必要に応じて得られた波形図を表示する。なお、バンドパスフィルタ処理部１６は、時系列データに対して直接バンドパスフィルタ処理を実行するだけでなく、必要に応じて高速フーリエ変換を行い、周波数空間でバンドパスフィルタ処理を行ってから逆高速フーリエ変換を行って時間空間に戻す処理を実行することもできる。移動平均処理部１８は、システムに取り込まれた時系列データに対して、あるいは既にシステムに取り込まれて何らかの処理がなされた時系列データに対して、隣接する所定の数の時間位置で平均をとる処理を実行する。この所定の数は、デフォルトで設定されているものをそのまま使用することも、後述するユーザインターフェイスを介してユーザが設定した値を使用することもできる。
【0025】
微分処理部２０は、システムに取り込まれた時系列データに対して、あるいは既にシステム内で何らかの処理がなされた時系列データに対して、時間微分演算を実行し、必要に応じて得られた波形図を表示する。時系列データを取り込む時間間隔（サンプリング周波数の逆数）が例えば１０００分の１秒であったとすると、この１０００分の１秒の間隔で隣り合う２つのデータの値の差を１０００分の１秒で割ることによって時間微分が得られる。ただし、この時系列ごとに１つ置き又は２つ置き等時間間隔を広げて微分処理を実行することもできる。このように有限の時間間隔で時間微分をとることを、階差とも呼ぶ。積分階差処理部２２は、システムに取り込まれた時系列データに対して、あるいは既にシステム内で何らかの処理がなされた時系列データに対して、積分演算を実行し、さらに得られた積分データに対して、階差処理を実行して、必要に応じて得られた波形図を表示する。
【0026】
減算処理部２４は、ある二つの波形データの間で減算処理を実行し、必要に応じて得られた波形図を表示する。パルス生成部２６は、時系列データに対して、その値が連続して正である時間の長さをパルス幅、連続して正である時間の長さをパルス高とする第１のパルス波形と、その値が連続して負である時間の長さをパルス幅、連続して負である時間の長さをパルス高とする第２のパルス波形とを求め、必要に応じて第１のパルス波形及び第２のパルス波形の両方又はいずれか一方の波形図を表示する。閾値処理部２８は、ユーザが閾値を設定することによって、複数のパスルの中からその閾値を超える、又はその閾値を下回るパルスを選択するという処理を実行し、必要に応じて選択されたパルスを表示する。分割点探索処理部３０は、ある時系列データのある始点から終点までの範囲を、その範囲に含まれるいくつかの点（分割点）を直線でつないだ直線からなる曲線で近似するときに、その時系列データの特性を判断するのに十分な近似ができるように、これらの分割点の位置を決めるという処理を実行する。
【0027】
処理範囲設定処理部３２は、システムに取り込まれた時系列データ、あるいは既にシステム内で何らかの処理がなされた時系列データの全体を処理対象とするのではなく、ある一定のルールに基づいて時間範囲を設定し、その設定された時間範囲の時系列データだけを処理対象とする処理を実行する。処理範囲設定処理部３２によって行う具体的な処理の内容については後に詳述する。
【0028】
図１に示す移動平均処理部１８及び積分階差処理部２２は、ノイズ除去処理を行うものであり、システムの構成として必ず必要なものというわけではなく、必要に応じてこれらのうちのいずれか一方又は両方を備えるようシステムを構成することができる。また、図１に示したすべての処理部を備えていなければ本発明が実現できないというものではない。例えば、微分処理部２０、範囲設定処理部３２、ピーク位置検出部３４だけを含んでいれば、時系列データに対する最低限必要な処理を行うシステムを構築することができるし、微分処理部２０、パルス生成部２６、減算処理部２４、閾値処理部２８だけを含んでいても、時系列データに対する最低限必要な処理を行うシステムを構築することができる。あとは、実際の処理の必要性に応じて図１に示すその他の処理部を追加してシステムを構築することができる。
【0029】
図２は、射出成形において、樹脂が型に流し込まれてから固化するまでの間に加わる圧力（縦軸）が時間（横軸）とともにどのように変化するかを示した、実際の圧力センサの出力の時系列データを示している。このようなデータの時系列データ解析システムへの取り込みは、図１のデータ入力処理部１０によって行われる。
【0030】
図２には、信号波形に重ねて４カ所に円１０ａ、１０ｂ、１０ｃ、１０ｄが示されている。図２に示す信号波形において、一番左の円１０ａは、圧力がほぼゼロの状態から急激に増大し始める時点であり、樹脂が型の内部全体に行き渡る直前の瞬間を示している。そして、１０ａの時点から圧力が急激に増大して円１０ｂのピークに到達する。このピークは、樹脂が型全体に充填され、樹脂の流し込みが停止された時点に対応する。樹脂の流し込みが停止されると圧力は急速に低下するとともに、円１０ｃの時点からしばらく圧力一定の状態が続く。この圧力一定の期間は、樹脂の温度が冷やされている期間に対応する。そして、樹脂が十分に冷やされて固化すると、円１０ｄの時点から再び急速に圧力が低下する。
【0031】
図２に示した型の内部の圧力の時系列的な変化の仕方は、型の大きさや種類、樹脂の組成等によって変わってくる。したがって、例えば図２の円１０ａ、１０ｂ、１０ｃ、１０ｄに対応する特徴的な変化が時系列の中で現れる時点（屈折点）を捉えることで、成形品の不良の発生や、射出成形作業の不具合等を知ることができる。すなわち、射出成形作業そのものや成形品の状態を、圧力センサによってモニタすることができる。
【0032】
ところで、射出成形作業や鋼板の切り出し作業を圧力センサや加重センサでモニタする場合には、大量のデータが得られる。このような大量のデータから上述のような特徴的な屈折点を正確に抽出するには、迅速な処理が必要となるので、処理の高速化が求められ、それだけ計算機に対する負荷が増大する。そこで、本実施形態では、以下のようなデータ処理を行う。
【0033】
上述の射出成形の例で、樹脂が型の内部全体に行き渡る瞬間（図２の円１０ａ）の前後の変化を、時間を大きく引き延ばすと、図３（ａ）のようになると考えられる。図３（ａ）に円１１ａで示す、樹脂が型の内部全体に行き渡る時点から圧力が上昇し始める。すなわち、この時点が屈折点となる。屈折点を知りたい場合に、まず図１に示した微分処理部２０により、図３（ａ）に示したグラフに対して時間微分をとると、図３（ｂ）のような時系列データの曲線が得られる。図３（ａ）の変化を変位（displacement）と考えれば、図３（ｂ）の波形は速度（velocity）に対応する。そして図３（ｂ）の曲線に対し、微分処理部２０によりさらに微分をとると、図３（ｃ）に示す曲線が得られる。これは加速度（acceleration）に対応する。この図３（ｃ）の曲線において、円１１ｂのような極大点が見つかれば、その時点が、樹脂が型の内部全体に行き渡った瞬間の屈折点であると判断することができる。
【0034】
そこで、図３（ｃ）と同じ図４（ａ）の曲線に対し、微分処理部２０によりさらにもう１回時間微分を行って、図４（ｂ）の曲線を得る。この図４（ｂ）の曲線がゼロを通過する時点（上向き矢印で示す）が分かれば、この時点が図４（ａ）の曲線の極大点、すなわち樹脂が型の内部全体に行き渡った時点であると判断することができる。しかしながら、図３（ａ）に示した曲線は、樹脂が型の内部全体に行き渡った場合の圧力の変化を想定した曲線に過ぎず、実際の圧力センサの信号波形はこのような滑らかなグラフにはならず、これにかなりのノイズが載った状態で出力される。
【0035】
図５は、ノイズが多く含まれる時系列データから極大点を求めようとした場合の問題点を示している。図５（ａ）は、図４（ａ）に対応する、極大点を求めようとする信号波形の一例を示しており、図４（ａ）の曲線と比較するとかなりのノイズが含まれている。そこで、まず、図１の移動平均処理部１８により移動平均をとる。この図５（ａ）の信号波形に対して、４９点移動平均をとったものが図５（ｂ）の波形であり、９９点移動平均をとったものが図５（ｃ）の波形であり、１９９点移動平均をとったものが図５（ｄ）の波形である。ここで、複数点移動平均は、階差すなわち微分を行ったことに相当するが、一般には、移動平均の点数が多くなればそれだけノイズは減少する。
【0036】
図５の（ｂ）、（ｃ）、（ｄ）の波形は、異なる点数で移動平均をとったものであるが、（ｂ）、（ｃ）、（ｄ）の波形がゼロを通過するそれぞれの時点ｔb、ｔc、ｔdは一致してはおらず、前後に多少ずれている。このことは、データの処理の仕方によって、すなわちこの場合だと、何点で移動平均処理を行うかによって、得られた波形がゼロを通過する時点が変わってしまうことを示している。したがって、上記図４に示した、加速度に対して微分（階差）をとるという方法では、ノイズが含まれる出力信号から、トレンドが変わる時点を正確に求めることはできない。
【0037】
そこで、上記の問題を改善するために、本実施形態では、積分階差という処理を実行する。すなわち、ゼロを通過する時点を求めたい波形に対して、ある等間隔の時間で積分を行い、この積分値に対して差分をとる。具体的には、
という計算を実行する。
【0038】
上記の積分階差の式について、図６Ａ及び図６Ｂを参照しながら説明する。まず、図６Ａ及び図６Ｂの（ａ）の波形（図５（ａ）の波形と同じである）において、いま注目している時点をｉとし、積分範囲をｗとすると、上記の式の右辺第１項は、波形ａの値をｉからｉ＋ｗまで積算したもの、すなわちｉからｉ＋ｗまでのａの曲線の下側の面積（ハッチングで示す）に対応する。同様に右辺第２項は、波形ａの値をｉ－ｗからｉまで積算したもの、すなわちｉ－ｗからｉまでの曲線の下側の面積（梨子地で示す）に対応する。したがって、積分階差は、ｉの右側の面積とｉの左側の面積との差を表している。ここで、ｗが９９であれば９９点積分階差、ｗが１９９であれば１９９点積分階差という。
【0039】
図６Ａ及び図６Ｂにおいて、（ｂ）の波形は、（ａ）の波形ａに対して、９９点積分階差を実行して得られた波形図を示している。そして、図６Ａでは、積分階差の値がゼロとなる場合に対応し、図６Ｂでは、右辺第１項と第２の差がほぼ極小となっている場合に対応している。図６Ａ及び図６Ｂにおける（ｂ）の波形を見ると分かるように、元の波形の特徴を大きく失うことなく、極大点、すなわちゼロを通過する時点を特定するのに支障がない範囲でノイズが除去されていることが分かる。図６Ａの（ａ）の波形を見ると、ｉの時点よりもやや左側にこの波形が実際にピークとなる時点がある。人がこのピークを見れば、このピークがノイズに起因するものだと判断することができるが、このような判断をコンピュータに実行させようとすることは非常に難しい。これに対して、上記の積分階差処理を実行して得られる図６Ａの（ｂ）の波形をみると、この波形がゼロを通過する時点が、図６Ａの（ａ）の波形曲線の全体的な極大点と一致していることが分かる。
【0040】
次に、図６Ａの（ｂ）に示すような波形がゼロを通過する時点をコンピュータによって判定する具体的な方法について説明する。図６Ａの（ｂ）に示す波形も、実際には飛び飛びの時間（横軸）に対してプロットし、これを曲線で結んだ時系列のデータである。図７（ａ）、（ｂ）は、上の波形図が、図６Ａの（ｂ）に対応するものであり、下の波形図は上の波形図を時間微分したものである。したがって、下の波形図がゼロになる時点が、上の波形図がゼロを通過する時点となる。しかし下の波形図も上の波形図と同様に飛び飛びの時間に対してプロットされるものである。このように時間解像度が有限であるため、図７（ａ）、（ｂ）の下に示す時系列データに完全なゼロが含まれるとは稀である。
【0041】
データがゼロとなる点を探索する場合には、図７（ａ）、（ｂ）の下の波形図に示すように、縦軸のゼロを中心とした一定幅のしきい値（横軸の時間軸を中心とする一定幅の帯状の部分）を設定し、この範囲に含まれるデータをゼロと判定することが考えられる。図７（ａ）、（ｂ）の下の図で、閾値の範囲にあるデータについて、分かりやすいように縦軸に平行な線で示している。この場合、閾値の幅が（ａ）のように狭いと、曲線が大きく変化する範囲に該当するデータが存在しない場合があったり、また、閾値の幅が（ｂ）のように広いと、該当するデータが多数存在して本来見つけたいデータを見つけることができないといった事態が生じる。このように、アプリケーションごとに適切な閾値を設定することは容易ではなく、効率も悪いという問題がある。別の方法として、時系列データを１点ずつ探索し、符号がプラスからマイナスへ、あるいはマイナスからプラスへ変わる時点を見つける方法も考えられる。しかしながら、この方法では処理がスカラー処理となるため、効率が悪い。
【0042】
そこで本実施形態では、次のようにして波形図がゼロを通過する時点を見いだす。図８（ａ）は、図７（ａ）、（ｂ）の上の波形図に対応し、図８（ｂ）は、図７（ａ）、（ｂ）の下の波形図に対応するもので、図８（ｂ）は図８（ａ）の波形を時間微分した波形である。この図８（ｂ）の波形について、図１に示したパルス生成部２６によって、図８（ｃ）、（ｄ）のようなパルス波形を導く。図８（ｃ）のパルス波形は、図８（ｂ）の波形において、連続して正である時間の長さをパルス幅、連続して正である時間の長さをパルス高として導いたものである。例えば、図８（ｂ）の中央の大きく変化する部分でＴ1の時間幅に対応する図８（ｃ）の部分にそのＴ1の値に対応する高さのパルスが示されている。一方、図８（ｄ）のパルス波形は、図８（ｂ）の波形において、連続して負である時間の長さをパルス幅、連続して負である時間の長さをパルス高として導いたものである。例えば、図８（ｂ）の中央の大きく変化する部分でＴ2の時間幅に対応する図８（ｄ）の部分にそのＴ2の値に対応する高さのパルスが示されている。図８の（ｃ）と（ｄ）は、縦軸の値は別として、時間軸で見ると互い違いになっている。
【0043】
次に、図９に示すように（図９（ａ）～（ｄ）は、図８（ａ）～（ｄ）に対応する）、図１の減算処理部２４によって、図９（ｄ）と図９（ｃ）を減算して差分をとると、図９（ｅ）の波形図が得られ、さらに図９（ｅ）の波形図に対して時間微分を行うと、図９（ｆ）に示すような、正から負へ切り換わった時点と負から正に切り換わった時点のみにパルスが残る波形図が得られる。そしてこのパルスが正か負かによって、正から負に切り換わったのか、負から正に切り換わったのかが分かる。図９（ｆ）の波形図の各ピークは、極大点の時間的な位置に対応し、各ピークの高さは各階差が同じ符号で連続する長さを反映している。なお、ここまでの処理はすべてベクトル処理であり、データを一つひとつ検索する必要がないため、計算処理の負荷が大幅に軽減される。
【0044】
さらに、図９（ｆ）の波形図に対して、前述の積分階差を実行すると、図１０（ｈ）の波形図が得られる（図１０の（ａ）は、型の内部の圧力を検出する圧力センサの出力波形であり、（ｂ）～（ｇ）は、図９の（ａ）～（ｆ）に対応する）。ここで、図１の閾値処理部２６により、ある一定の閾値を設定しておいて、これを超えるパルスのみを選択すれば、この時点が図１０（ａ）の時系列データの屈折点、すなわちトレンドが変化した時点であることが分かる。
【0045】
次に、ノイズの含まれるデータから実際のピークとなる時点を探索する例について説明する。図１１において、（ａ）は実際の型の内部の圧力を検出する圧力センサの出力データであり、このかなりのノイズが含まれるデータから圧力の最大点を求める。図１１において、横軸は時間である。この図１１（ａ）の出力データに対して、高速フーリエ変換（ＦＦＴ）を実行し、得られたスペクトルに対してバンドパスフィルタで不要な成分を除去し、さら逆ＦＦＴを実行して得られたのが、図１１（ｂ）に示す波形である。なお、図１１（ｂ）において、「０－４０Ｈｚ」は、４０Ｈｚより大きい周波数の波形をカットして０～４０Ｈｚの周波数の波形だけを取り出すことを意味する。
【0046】
そして、図１１（ｂ）の波形に対して、９点の移動平均をとったものが図１１（ｃ）に示した波形である。この図１１（ｃ）の波形の最大点が、圧力の最大時点と考えることができる。ところで、図１１（ａ）の波形と図１１（ｃ）の波形を比較すると、図１１（ａ）の実際のピークの位置は、図１１（ｃ）のピークの位置よりも左側にずれている。これは、ノイズの影響によるものと考えられる。このようにノイズによってピークの位置がずれるような場合であっても、図１１の（ａ）→（ｂ）→（ｃ）の処理を行うことによって、ノイズの影響を排除して圧力が最大となる時点を求めることができる。この図１１の（ａ）→（ｂ）→（ｃ）の処理を行う部分は、機能的には図１に示すピーク位置検出部３４に対応する。
【0047】
次に、図２の圧力センサの出力波形の円１０ｃに当たる時点、すなわち樹脂が充填され圧力がピークになった直後に、圧力が急激に低下したあとにほぼ一定の安定な期間が始まる時点を見いだす処理について説明する。図１２において、（ａ）は図２に示した圧力センサの出力波形を、横軸の時間軸の方向に伸ばした状態を示している。図１２の（ａ）の波形に対して、図１のバンドパスフィルタ処理部１６によって、図１１の（ｂ）と同様のバンドパスフィルタ処理を施して得られたのが、図１２の（ｂ）に示す波形である。
【0048】
そして、（ｂ）の波形に対して、図１の移動平均処理部１８によって９点の移動平均をとって得られたのが（ｃ）に示す波形である。そしてさらに、（ｃ）の波形に対して微分処理部２０によって時間微分を行って得られたのが（ｄ）に示す波形であり、これにさらに時間微分を行って得られたのが（ｅ）に示す波形である。この（ｅ）に示す波形には、横軸の値が６００のやや手前にピークが見られるが、（ａ）の波形と対照させると、この時点が、安定な期間の開始する時点であることが分かる。
【0049】
続いて、図２の圧力センサの出力波形の円１０ｄに当たる時点、すなわち樹脂が十分に冷やされて固化して圧力が低下する時点を見いだす処理について説明する。図１３において、（ａ）は図２に示した圧力センサの出力波形を、横軸の時間軸の方向に伸ばした状態を示しており、（ａ）の円で示した部分が、樹脂が十分に冷やされて固化して圧力が低下する時点に該当する。図１３（ａ）の波形に対して、図１のバンドパスフィルタ処理部１６によって、図１１の（ｂ）と同様のバンドパスフィルタ処理を施して得られたのが、図１３の（ｂ）に示す波形である。
【0050】
そして、図１３の（ｂ）の波形に対して、図１の移動平均処理部１８によって９点の移動平均をとって得られたのが（ｃ）に示す波形である。そしてさらに、（ｃ）の波形に対して微分処理部２０によって時間微分を行って得られたのが（ｄ）に示す波形である。ここでは、この（ｄ）に示す波形から、あるほぼ一定の正の値から一時的に負になる時点が、図１３（ａ）の円で示した部分に対応することが分かる。なお、この例では、（ｄ）の波形に対してさらに時間微分を行って得られる（ｅ）の波形、すなわち加速度に対応する波形では、下降の動きが緩やかであるため、この時点をうまく検出することは難しい。
【0051】
これまで説明してきたように、図１に示した各処理部により、時系列データに対して各種の処理をさまざまに適用し実行することによって、時系列データの中の屈折点、すなわちトレンドが変化する点を高い精度で見いだすことができる。図１４は、本実施形態のシステムでできることをまとめた図である。データ処理としてできる基本的な処理は、バンドパスフィルタ処理、移動平均処理、階差積分処理、微分処理の４つで、これらを組み合わせることによって各種のデータ処理が可能となる。そして、処理がなされたデータに対して、最大値、最大値の位置（時間）、最小値の位置（時間）、閾値処理、積分のうちの１つ又はいくつかを求めることによって屈折点、すなわちトレンドの変化点を特定することができる。
【0052】
図１５は、図１に示した時系列データ解析システムのユーザインターフェイスの一例を示している。このようなユーザインターフェイスを利用することにより、ユーザが実際に時系列データに対する処理を実行する際の利便性が向上する。このユーザインターフェイスを介して、任意の時系列データに対して上述の各処理を実行し、その処理結果を直ちに視覚的に見ることができる。この画面の左上の「データソース」欄には、処理対象となる時系列データの名称が記載され、このすぐ右側の「Ｒａｗ ｄａｔａ」欄に、圧力センサや加重センサの出力データなどの、処理の対象となるコンピュータに入力された実際の時系列データが表示される。
【0053】
さらに、その右側には「ＦＦＴ ｓｐｅｃｔｒｕｍ」欄が設けられ、ここには「Ｒａｗ ｄａｔａ」欄に表示されている波形図に対して高速フーリエ変換（ＦＦＴ）が行われたＦＦＴスペクトルが表示される。このＦＦＴ処理は、「Ｒａｗ ｄａｔａ」欄に表示される入力データすべてに対して実行することもできるし、ユーザがＦＦＴ処理を実行するかどうかを選択するようにもできる。なお、この「ＦＦＴ ｓｐｅｃｔｒｕｍ」欄では、バンドパスの設計が可能である。図１５の「ＦＦＴ ｓｐｅｃｔｒｕｍ」欄の波形を見ると、横軸の５０Ｈｚ前後にピークが集中していることが分かる。そこで、ここではバンドパスの上限として２５Ｈｚが指定されており（２５Ｈｚより高い周波数部分のグレーの網かけがそのことを示している）、２５Ｈｚ以下の波形だけが抽出されている。
【0054】
また、図１５において、「データソース」欄の下には、サンプリングレートを選択する欄が設けられている。このサンプリングデータは、時系列データの隣り合う個々のデータの時間間隔に対応するものであり、ここでは「１００Ｈｚ（１００分の１秒間隔）」、「１０００Ｈｚ（１０００分の１秒間隔）」、「１００ｋＨｚ（１０万分の１秒間隔）」のうちのいずれかを選択することができるようになっている（図１５では１０００Ｈｚが選択されている）。その下の「検出モード」欄では、時系列データに対して変位（displacement）、速度（velocity）、加速度（acceleration）のいずれかを選択して表示させることができる。さらにその下の「バンドパス下限（Ｈｚ）」及び「バンドパス上限（Ｈｚ）」では、バンドパス処理の上限値及び下限値を、画面上でグラフィック的に設定できるようになっている。
【0055】
その下の「移動平均範囲」欄は、移動平均をとる範囲、すなわち平均をとる連続する時系列データの数を設定するものであり、１から１９までの範囲で設定できる。その下の「ｃｅｎｔｅｒ」欄は、移動平均の中央を選択するスイッチである。その下の「屈折点検出」欄は、屈折点検出を実行するかどうかをチェック欄にチェックするかしないかで設定する。その下の「ピーク方向」欄は、ピークを検出する場合に、上方向のピーク及び下方向のピークのいずれか又は両方を検出することを選択する欄である。「ピーク幅」欄は、ピークを検出する場合に、ここで設定した幅以下の幅のピークだけ検出するように設定する欄である。「ピーク閾値（最大ピークに対する割合）」は、ここで設定した閾値を超える値を有するもののみをピークとして検出するように設定する欄である。そして、図１５の右下で最も広い領域には、これまでの各欄で設定した内容に基づいて入力された時系列データを処理して得られる波形図が示されている。
【0056】
以上の説明から理解されるように、時系列データが特徴的な変化を示す屈折点は、変位、速度、加速度のいずれかのピークとして見いだすことができる。そして、変位から速度を求める方法、及び、速度から加速度を求める実際の方法は、実装上は時間微分、すなわち階差をとることによる。しかしながら、入力された時系列データの全体に対して階差をとる処理を行うと、その時系列データにおける長期のトレンドが失われることになる。長期のトレンドが失われることは、短期の変動、すなわちノイズを強調することになって、特徴点の抽出の妨げとなる。図１６はこのことを説明する図であり、（ａ）は、ある荷重の変化を示す時系列データ、（ｂ）は、（ａ）の時系列データを時間微分して得られる速度データ、（ｃ）は（ｂ）の速度データを時間微分して得られた加速度データを示している。図１６の（ｃ）に示すように、（ａ）の時系列データ全体に対して時間微分を求める処理を連続して行うと、ノイズが強調されて、（ａ）に円で示したトレンドが変化した特徴点を見つけることが困難となる。
【0057】
また、図１７は、移動平均の度合いを強めることによってノイズを低減できるとこを示している。図１７の（ａ）は、速度及び加速度に対して５点移動平均をとったもの、（ｂ）は９点移動平均をとったもの、（ｃ）は１９点移動平均をとったものを示している。これらの図で、図の下に示した上向きの矢印は、特徴点と判定された時点を示しているが、（ａ）では、荷重が低下し始める特徴点ではなくそれよりも右側にずれた時点を特徴点として抽出している。（ｂ）では、（ａ）に比べてノイズは低減されており、特徴点を検出できているが、たまたま検出できた可能性が高く、堅牢性に問題がありうる。（ｃ）では、ノイズはかなり除去されており、ノイズはさらに低減されて特徴件が検出されている。
【0058】
しかしながら、移動平均の点数を増やすと、図１８に示すような問題が生じうる。図１８において、実線が荷重の変動を示す実際の曲線を、点線は元の荷重のデータに対して１９点移動平均を取った曲線を、一点鎖線は元の荷重のデータに対して３９点移動平均を取った曲線をそれぞれ示している。図１８において、例えば荷重のかかり始めを示す曲線の立ち上がりを知りたいときに、３９点移動平均を取ってしまうと、実際の立ち上がりの時点から左側（時間の早い側）に大きくずれてします。このように、移動平均を取る点数を増やしてゆくと、元の曲線から大きく波形が変化してしまい、単純に移動平均の点数を増やしてゆくだけでは正確な特徴点が難しくなる。
【0059】
図１９の一番上のグラフは、射出成形において、樹脂が型に流し込まれてから固化するまでの間に加わる圧力（縦軸）が時間（横軸）とともにどのように変化するかを示した圧力センサの出力の時系列データのグラフを示しており、２番目のグラフは一番上の時系列データに対し、一回の微分を行って速度を求めた時系列データのグラフであり、一番下のグラフは、さらにもう一回微分して加速度を求めた時系列データを示している。一番上のグラフでは、横軸の値が２５０の辺り（ｔ1で示す）で縦軸の値が急激に上昇し、３００を少し過ぎた辺り（ｔ2で示す）でピークに達して低下し始め、４００を少し超えた辺り（ｔ3で示す）で低下速度が急速に減少してからしばらくほぼ一定の値が続いている。このように、横軸が２００から５００までの間に大きく３つのトレンドの変化が認められ、それぞれに意味が異なっている。
【0060】
圧力センサからの出力が、図１９の一番上のグラフのような挙動を示すことを分かっていれば、例えばｔ3の特徴点を抽出しようとするオペレータは、横軸が４００を少し超えた辺りの前後の適当な範囲に限定して、図３から図１３に示したような処理を実行するようにすれば、特徴点の抽出処理を効率化することができるとともに、精度も向上する。そこで、本実施形態では、図１９のような時系列データに対して特徴点の抽出処理を実行するときに、図１９の時系列データを表示させた状態で、横軸が４００を少し超えた辺りの前後にグレーの時間幅を手動で直接設定することによって、その設定範囲についてのみ図３から図１３に示したような処理を実行するようシステムに指示することができる。
【0061】
図２０は、横軸の時間軸ではなく、縦軸の値域に着目して範囲を決める方法を示している。図２０に示すグラフは、図１９に示したグラフと同じであるが、２番目に示した速度のグラフにおいて、縦軸の値がマイナス０．１以下になる範囲を設定し、この設定範囲についてのみ図３から図１３に示したような処理を実行するようシステムに指示する。このようにすることによって、横軸が３００前後になる加速度が大きく変化する部分を、特徴点の探索処理から排除することができる。このように、縦軸の値に基づいて閾値を指定し、それよりも元の値が小さい範囲又は大きい範囲を指定することで、本当に見つけたい特徴点を確実に見つけることができる。
【0062】
図２１は、図１９、図２０と同じグラフであり、ここでは、図１９のｔ3に対応する時点の前後で加速度が最大となる時点を見つけようとしている。そこで、時刻ｔ3の前後のある範囲（ここでは１５０ポイントの範囲）に限定し図３から図１３に示したような処理を実行するようシステムに指示する。このようにして特徴点の特性に基づいて検索する範囲を設定することによっても、特徴点の抽出処理を効率化することができ、精度も向上する。
【0063】
図２２は、図１９～図２１に示したセンサセンサの時間的変化を示すグラフのうち、加速度最小点と速度最大点の間を検索範囲とし、この範囲で加速度が最大となる点を検索するたに、図３から図１３に示したような処理を実行するようシステムに指示する例である。
【0064】
次に、複数の時系列データのお互いの関係性に基づいて処理範囲を設定する処理について説明する。図２３の上のグラフには、同時に取り込まれた３つの時系列データが示されている。これらは、射出成形の型の３カ所に設けられた圧力センサからの出力が、樹脂が型に重点される前後で圧力が時間的にどのように変化するかを示したグラフである。３つの圧力センサから同時に３チャネルのデータが発生し、システムに入力される。図２３の上に示すグラフから分かるように、同じ一つの型の内部でも、場所によって圧力の変化が同一ではないが、全体としての変化の傾向は類似していることが分かる。この３つセンサからのデータの関係性をみて、範囲を設定することについて説明する。
【0065】
図２３の上に示した３チャネルの時系列データから、圧力上昇が開始される点を見つける場合を考える。図２３の下に示すグラフは、上に示す３つのデータについての分散（３つの時系列データの各時点における平均値とそれぞれのデータとの差分をとり、さらに各差分を２乗して合計したもの）を求めたものである。図２３の下のグラフを見ると、樹脂が充填される前はほとんど変化がない。圧力センサが設けられた位置は互いに異なるため、圧力が上昇し始める時刻が異なるため、３つの圧力が上昇を始める直前に分散がなり最大となり、そしていったん樹脂が充填されると、圧力センサの出力は再び類似した挙動に戻る。そこで、図２３の下に示す例では、分散が最初にピークとなる時点の前後の－５０～＋１０ポイントの範囲を検索対象として設定すれば、３つの圧力センサの立ち上がり点を検索範囲とすることができる。
【0066】
次に、図２４～２６を参照して、ＤＳＬを用いて、実際に時系列データに対して処理範囲を設定する具体的な方法について説明する。図２４は、下部に図１９に示した時系列データのグラフを示し、上部に、下の時系列データに対してオペレータが処理範囲を指定するためのデータ入力のためのＤＳＬの一例を示している。この例では、オペレータが、図の下部の時系列データを見ながら「HORIZONTAL_LIMIT = 」の右側の括弧内に「３８０」及び「４５０」と入力した状態を示している。この操作により、時系列データには、横軸（時間軸）が３５０から４５０までの範囲にグレーの帯が表示され、この範囲が処理範囲として指定される。このグレーの帯の幅を見て適切でないと判断したときは、「HORIZONTAL_LIMIT = 」の右側の括弧内の数字を変更して適切な範囲を再度設定することができる。なお、「ROLLING_MEAN = 5」は、移動平均の点数が「５」であることを示しており、「VERTICAL_LIMIT = [None, None]」は、垂直方向の範囲指定をしないことを示しており、「TARGET = IDXMIN(ACC)」は、加速度（ＡＣＣ）が最小になる位置を検索することを指示したことを示している
【0067】
図２５は、図２０に示した時系列データに対して、縦軸の値域に着目してＤＳＬを用いて処理範囲を設定する具体的な方法を示している。図２５では、「HORIZONTAL_LIMIT = [None, None]」となっており、横軸について範囲設定は行わないが、代わりに「VERTIAL_LIMIT = [None, -0.1，VCT]」と設定されて、２番目の速度の縦軸の値が－０．１以下の範囲について検索することを指示したこと示している。なお、「ROLLING_MEAN = ５」は、移動平均の点数が「５」であることを示しており、「TARGET = IDXMAN(ACC)」は、加速度が最大にある位置を検索することを指示したことを示している。
【0068】
図２６は、の下部には、図２２に示した時系列データに対して、加速度最大点を検索するための、ＤＳＬを用いた処理範囲の設定の仕方を示している。図２６において、「HORIZONTAL_LIMIT = [IDXMAX(DST), IDXMAX(DST)+150]」とあるのは、元波形（ＤＳＴ）の値域が最大となる点を検索範囲の左端とし、元波形の縦軸の値が＋１５０である点を検索範囲の右端として設定したことを示している。また、「TARGET = IDXMAX(ACC)」とあるのは、加速度が最大となる位置を検索することを指示したことを示している。
【0069】
以上のように、オペレータが時系列データの特性を知っている場合には、その特性に合わせて容易に範囲を設定して、時系列データの全体ではなく設定された一部の範囲だけをシステム検索させることができる。このようにすることによって、検索処理効率及び検索精度を向上させることができ、合わせて図１８に関連して説明した問題点を回避することができる。
【0070】
次に、時系列データに対し少ない数の直線を分割点で連結して時系列データにより良くフィットする近似曲線を求める実施形態につい手説明する。図２７は、プレス加工の際に使用する歪みゲージというセンサからの出力データの波形を示している。図２７の出力データの波形に矢印ｘ1、ｘ2で示した２つのピークの間の波形の様子を見ることによって、プレス加工がうまくいったか、あるいは何か不都合が生じたか等を知ることができることが経験的に知られている。そこで、プレス加工がうまくいったかどうかを自動的に判定するために、２つのピークであるｘ1とｘ2の間をなるべく少ない数で分割し、その分割点を直線でつなぐことによって最も良く出力データの波形にフィットする近似曲線を求め、この近似曲線に基づいて上記のような判断を行う。そのために、本実施形態では、図１に示した時系列データ解析システムの分割点探索処理部３０により、以下のような処理を実行する。
【0071】
図２８に示した波形は、図２７に示した波形と同じではないが、大体の傾向は一致しているので、図２８を用いて、分割点探索処理について説明する。図２８において、始点ｘ1と終点ｘ2という２つのピークの間にｐ1、ｐ2、ｐ3という３つの分割点を波形データの上に仮に置き、その間を直線でつなぐ。分割点ｐ1、ｐ2、ｐ3の仮の位置は、分割点探索処理部２４により、例えばｘ1とｘ2の間に均等に配置することもできるが、それ以外の仮の配置をユーザが選択して設定することもできる。
【0072】
ｐ1、ｐ2、ｐ3の最初の位置が決まれば、ｘ1とｐ1とをつなぐ直線ａ、ｐ1とｐ2とをつなぐ直線ｂ、ｐ2とｐ3とをつなぐ直線ｃ、ｐ3とｘ2とをつなぐ直線ｄという４つの直線ａ、ｂ、ｃ、ｄが得られる。分割点ｐ1、ｐ2、ｐ3はあくまでも仮に置かれた点なので、この時点の４つの直線ａ、ｂ、ｃ、ｄからなる曲線が、図２７の出力データの波形を最も良く表している可能性は低い。しかしながら、分割点ｐ1、ｐ2、ｐ3の数がこのように少数であっても、これらをｘ1、ｘ2の２点間のすべての点で任意に移動させて最小二乗法により直線で近似する回帰演算は、アルゴリズム自体は短いが、複数の分割点を組み合わせることで、計算量は大幅に増大する。
【0073】
そこで、３つの分割点ｐ1、ｐ2、ｐ3を移動させて、図２７の出力データの波形を最も良く表す４つの曲線ａ、ｂ、ｃ、ｄを求めることを考える。図２９は、３つの分割点ｐ1、ｐ2、ｐ3の位置を最終的にどのように決めるかを説明するための図である。まず、図２９（ａ）では、ｘ1と、仮に決められたｐ2の位置を固定する。図２９（ａ）の横軸の目盛りでは、ｘ1は２１０の辺り、ｐ2は２４０の辺りに位置している。この間でｐ1を決められた間隔、例えば横軸の目盛りで１ずつ移動させ、各点において、ｘ1とｐ1とを結ぶ直線ａが、対応する範囲の出力データの波形と最も良くフィットする点としてｐ1の位置を決定する。
【0074】
ここで、最もよくフィットする点とは、例えば直線ａと出力データ波形とを最小二乗法により近似する回帰演算を行って、回帰の誤差が最小となるｐ1の点として決定する。ここまでの処理では、横軸の２１０から２４０までの３０点で計算が必要となるが、分割点ｐ1は１つだけなので、ループ処理は１回で済み、コンピュータの負荷はそれほど大きくはない。
【0075】
次に、図２９（ｂ）に示すように、決定したｐ1から、仮に決められたｐ3までの間でｐ2を１ずつ移動させ、各点においてｐ1とｐ2とを結ぶ直線ｂが、対応する範囲の出力データの波形と最も良くフィットする点としてｐ2の位置を決定する。さらに、図２９（ｃ）に示すように、決定したｐ2からｘ2までの間でｐ3を１ずつ移動させ、各点においてｐ2とｐ3とを結ぶ直線ｃが、対応する範囲の出力データの波形と最も良くフィットする点としてｐ3の位置を決定する。以上のようにして決定したｐ1、ｐ2、ｐ3を結ぶ直線ａ、ｂ、ｃ、ｄによって得られる近似曲線は、出力データの波形と比較してかなりよく近似された曲線となっている。これにより、直線ａ、ｂ、ｃ、ｄによって得られる近似曲線に基づいて、プレス加工がうまくいったかどうかを判定することができる。しかも、上で説明した処理は、移動させる点が１つだけなので、計算負荷はそれほど高くなく、通常のパーソナルコンピュータでも十分に実行可能である。
【0076】
なお、図２８及び図２９に示した処理を一回行っただけでは十分にフィットする曲線が得られない場合には、このような処理を複数回行って、回帰誤差が最も小さくなる回の曲線を採用することもできる。そのような場合であっても、一回の処理にかかるコンピュータに対する負荷及び要する時間が大きくないので、全体としてコンピュータに対する負荷及び要する時間を削減することができる。
【0077】
例えば、射出成形では金型内部に圧力センサを設けて、実際に射出成形を行う際に内部の圧力を圧力センサで検知し、その出力データの波形を得ることができる。金型内部の圧力は、液体状の樹脂を注入している段階では高まるが、樹脂の注入が終わり、内部の樹脂が冷えて固まるに従って圧力は徐々に低下し、最終的に金型から成形品を取り出す段階では内部の圧力は急減する。射出成形の工程が適正に行われたかどうかを、このような出力センサの出力データの波形を見ることによって確認することができる。このような場合に、本実施形態のシステムを活用することができる。
【0078】
上述の例では、分割点の数はｐ1、ｐ2、ｐ3の３点だったが、このような圧力変化の特性が予め分かっていれば、それに基づいて、分割点の数を増減させてその用途に合う適切な数とすることができる。また、本実施形態を適用できるデータは、時系列のデータであれば、どのようなものにも適用できる。例えば、図２に示した圧力センサの出力データも時系列データであることら、第２実施形態における適切な数の分割点を使ってフィットする曲線を求め、それに基づいて鋼板の切断作業が適切に行われたかどうかを判定する、という用途に本実施形態のシステムを用いることができる。
【符号の説明】
【0079】
１０ データ入力部
１２ ユーザ入力処理部
１４ 画像出力処理部
１６ バンドパスフィルタ処理部
１８ 移動平均処理部
２０ 微分処理部
２２ 積分階差処理
２４ 減算処理部
２６ パルス生成部
２８ 閾値処理部
３０ 分割点探索処理部
３２ 範囲設定処理部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6A】
【図6B】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
【図28】
【図29】
Copyright © 2010-2024 Science Impact Inc. All Rights Reserved.