JP2022148844A - 計測データ同期方法、試験方法、及びコンピュータプログラム - Google Patents

Abstract

【課題】複数の計測器から出力される時系列データを精度良く同期できる計測データ同期方法を提供すること。【解決手段】計測データ同期方法は、第１～第５計測器から出力される時系列をデータ解析装置によって同期する方法であって、定常区間と立ち上がり区間とを含む基準信号を同時に第１～第５計測器へ入力したときに、第１～第５計測器から出力される第１～第５基準データを取得するステップＳＴ１と、第１～第５基準データに対し所定の時間幅内の中央値を抽出するメディアンフィルタ処理を施すことにより、第１～第５フィルタ処理後データを生成するステップＳＴ３～ＳＴ５と、第１～第５フィルタ処理後データの各々から立ち上がりタイミング情報として第１～第５立ち上がり時刻を算出するステップＳＴ６と、第１～第５立ち上がり時刻に基づいて第１～第５計測器から出力される第１～第５計測データを同期するステップＳＴ７と、を備える。【選択図】図３

Description

本発明は、計測データ同期方法、試験方法、及びコンピュータプログラムに関する。より詳しくは、複数の計測器から出力される時系列データをコンピュータによって同期する計測データ同期方法、試験方法、及びコンピュータプログラムに関する。

近年、二次電池式電気自動車、燃料電池自動車、ハイブリッド自動車、及びプラグインハイブリッド自動車等、電気エネルギを利用して走行する電動車両の開発が盛んである。電動車両におけるエネルギの流れは、走行中に目まぐるしく変化する。このため電動車両の性能を評価するための試験は、内燃機関で発生した動力のみによって走行する内燃機関車両の試験と比較して、より多くの計測器を必要としかつ手順も複雑である。例えばハイブリッド車両の試験では、車両の環境情報、速度、走行距離、電流、電圧、電力量、及び排ガス関連データ（排ガス濃度、排ガス流量等）等、様々なデータを取得するための複数の計測器が組み合わせて用いられる。

電動車両の試験では、試験の内容に応じて各計測器を作動させるタイミングが異なる。また一般的に用いられる計測器には、他の計測器と同期して作動させる機能が無いものも多い。このため各計測器によって得られた複数の計測データをコンピュータによって解析するためには、複数の計測データを同期させるための前処理を予め行っておく必要がある。従来、このような計測データの同期はオペレータによる手作業で行っていたため、作業に時間がかかり、また試験結果にばらつきが生じるおそれもあった。

特開２０１５－１９０８６１号公報

ところで１次元の時系列データである各計測データから、例えば図８に示すような信号の立ち上がり点を抽出することができれば、この立ち上がり点を基準として各計測データを同期させることができる場合がある。

図８には、ステップ状の基準信号（細一点鎖線参照）を計測器に入力した場合における、計測器の出力（細実線参照）の一例を示す。また図８には、細実線で示す計測データに対し、特許文献２に示す処理を施して得られる計測データを太破線で示す。

図８において細実線で示すように、一般的に計測器の計測データには、高調波のノイズやオフセットが含まれているため、立ち上がり点を精度良く抽出することができない。そこで特許文献１に示すように、計測データに対し微分処理、積分処理、及びローパスフィルタを用いた高調波除去処理を施すことにより、計測データから高調波のノイズやオフセットを除去することが考えられる。しかしながら図８において太破線で示すように、計測データに対しローパスフィルタを適用し、高周波成分を除去すると、立ち上がりが遅れてしまい、立ち上がり点を精度良く抽出することができず、ひいては計測データを精度良く同期させることができなくなってしまう。

本発明は、複数の計測器から出力される時系列データを精度良く同期できる計測データ同期方法、試験方法、及びコンピュータプログラムを提供することを目的とする。

（１）本発明に係る計測データ同期方法は、複数の計測器（例えば、後述の計測器３１～３５）から出力される時系列データをコンピュータ（例えば、後述のデータ解析装置５）によって同期する方法であって、定常区間と立ち上がり区間とを含む基準信号を同時に前記複数の計測器へ入力したときに、各計測器から出力される時系列データである基準データを取得するステップと、前記複数の基準データに対し所定の時間幅内の中央値を抽出するメディアンフィルタ処理を施すことにより、複数のフィルタ処理後データを生成するステップと、前記複数のフィルタ処理後データの各々から立ち上がりタイミング情報を取得するステップと、前記複数の立ち上がりタイミング情報に基づいて前記複数の計測器から出力される時系列データを同期するステップと、を備えることを特徴とする。

（２）この場合、前記フィルタ処理後データを生成するステップでは、前記フィルタ処理後データのうち前記定常区間に含まれる信号値が評価時間にわたり所定範囲内に収まるように前記時間幅を調整することが好ましい。

（３）この場合、前記基準データを取得するステップでは、前記立ち上がり区間においてランプ状に変化する信号を前記基準信号として用い、前記フィルタ処理後データを生成するステップでは、前記基準データに対し微分処理を施した後、前記メディアンフィルタ処理を施すことが好ましい。

（４）この場合、前記立ち上がりタイミング情報を取得するステップでは、前記フィルタ処理後データにおいて信号値が所定の閾値を超える時刻を立ち上がり時刻として取得することが好ましい。

（５）この場合、前記基準データを取得するステップでは、前記立ち上がり区間においてランプ状に変化する信号を前記基準信号として用い、前記立ち上がりタイミング情報を取得するステップでは、前記フィルタ処理後データのうち前記立ち上がり区間内に含まれる複数の信号値に基づいて生成される一次関数の傾き及び切片を算出し、前記傾き及び前記切片に基づいて前記フィルタ処理後データの立ち上がり箇所の遅れ時間を算出することが好ましい。

（６）本発明に係る試験方法は、複数の計測器（例えば、後述の計測器３１～３５）から出力される時系列データをコンピュータ（例えば、後述のデータ解析装置５）で処理することによって供試体の性能を評価する方法であって、定常区間と立ち上がり区間とを含む基準信号を同時に前記複数の計測器の各々の基準信号用チャンネルへ入力したときに、各計測器から出力される時系列データである基準データと、前記供試体の状態と相関のある信号を前記複数の計測器の各々の計測用チャンネルへ入力したときに、各計測器から出力される時系列データである計測データと、を取得するステップと、前記複数の基準データに対し所定の時間幅内の中央値を抽出するメディアンフィルタ処理を施すことにより、複数のフィルタ処理後データを生成するステップと、前記複数のフィルタ処理後データの各々から立ち上がりタイミング情報取得するステップと、前記複数の立ち上がりタイミング情報に基づいて前記複数の計測データを同期するステップと、前記同期した複数の計測データに基づいて前記供試体の性能を評価するステップと、を備えることを特徴とする。

（１）本発明では、定常区間と立ち上がり区間とを含む基準信号を同時に複数の計測器へ入力したときに各計測器から出力される時系列データである基準データを取得し、複数の基準データに対し所定の時間幅内の中央値を抽出するメディアンフィルタ処理を施すことによって複数のフィルタ処理後データを生成し、これらフィルタ処理後データの各々から立ち上がりタイミング情報を取得し、これら立ち上がりタイミング情報に基づいて、複数の計測器から出力される時系列データを同期する。特に本発明では、各基準データに対しメディアンフィルタ処理を施すことにより、基準データに含まれる信号の立ち上がりに遅れを生じさせることなく高周波ノイズを除去することができる。よって本発明によれば、メディアンフィルタ処理を経て生成されたフィルタ処理後データを用いることにより、精度良く立ち上がりタイミング情報を取得できるので、各計測器から出力される時系列データを精度良く同期させることができる。

（２）本発明では、フィルタ処理後データのうち定常区間に含まれる信号値が評価時間にわたり所定範囲内に収まるように、メディアンフィルタ処理における時間幅を調整することにより、立ち上がりに遅れを生じさせることなく、インパルス性ノイズのような突発的なノイズを適切に除去することができるので、結果として計測器の時系列データを精度良く同期させることができる、

（３）本発明では、立ち上がり区間においてランプ状に変化する信号を基準信号として用い、基準データに対し微分処理を施した後、メディアンフィルタ処理を施すことによってフィルタ処理後データを生成する。これにより、各計測器の計測データにおける立ち上がりに遅れを生じさせることなく、これら計測データに含まれるオフセット成分を除去することができるので、結果として計測器の時系列データを精度良く同期させることができる。

（４）本発明では、フィルタ処理後データにおいて信号値が所定の閾値を超える時刻を立ち上がり時刻として取得することにより、簡易な演算によって立ち上がり時刻を取得することができる。

（５）本発明では、立ち上がり区間においてランプ状に変化する信号を基準信号として用い、フィルタ処理後データのうち立ち上がり区間内に含まれる複数の信号値に基づいて生成される一次関数の傾き及び切片を算出し、これら傾き及び切片に基づいてフィルタ処理後データの立ち上がり箇所の遅れ時間を算出する。本発明によれば、このようにして算出された立ち上がり箇所の遅れ時間を用いることにより、各計測器の時系列データをさらに精度良く同期させることができる。

（６）本発明では、定常区間と立ち上がり区間とを含む基準信号を同時に複数の計測器の各々の基準信号用チャンネルへ入力したときに各計測器から出力される時系列データである基準データと、供試体の状態と相関のある信号を複数の計測器の各々の計測用チャンネルへ入力したときに各計測器から出力される時系列データである計測データと、を取得する。また本発明では、複数の基準データに対しメディアンフィルタ処理を施すことによって複数のフィルタ処理後データを生成し、これらフィルタ処理後データの各々から立ち上がりタイミング情報を取得し、これら立ち上がりタイミング情報に基づいて、複数の計測器から出力される時系列データを同期し、これら同期した複数の計測データに基づいて供試体の供試体の性能を評価する。特に本発明では、各基準データに対しメディアンフィルタ処理を施すことにより、基準データに含まれる信号の立ち上がりに遅れを生じさせることなく高周波ノイズを除去することができる。よって本発明によれば、メディアンフィルタ処理を経て生成されたフィルタ処理後データを用いることにより、精度良く立ち上がりタイミング情報を取得できるので、各計測器から出力される計測データを精度良く同期させることができ、ひいては精度良く供試体の性能を評価することができる。

本発明の第１実施形態に係る試験システムの構成を示す図である。 基準信号の一例を示す図である。 基準信号の一例を示す図である。 計測データ同期方法及び試験方法の具体的な手順を示すフローチャートである。 メディアンフィルタ処理の効果を説明するための図である。 第１～第５フィルタ処理後データを同一の時間軸上にプロットしたものの一例である。 本発明の第２実施形態に係る計測データ同期方法及び試験方法の具体的な手順を示すフローチャートである。 第１及び第２フィルタ処理後データの立ち上がり箇所の遅れ時間を算出する手順を説明するための図である。 ステップ状の基準信号を計測器に入力した場合における計測器の出力の一例を示す図である。

＜第１実施形態＞
以下、本発明の第１実施形態について、図面を参照しながら詳細に説明する。
図１は、本実施形態に係る計測データ同期方法及び試験方法が適用された試験システムＳの構成を示す図である。

試験システムＳは、二次電池式電気自動車、燃料電池自動車、ハイブリッド自動車、及びプラグインハイブリッド自動車等、電気エネルギを利用して走行する電動車両Ｖを供試体とし、この電動車両Ｖの性能を評価する装置である。

試験システムＳは、電動車両Ｖが載置されるシャシダイナモメータシステム１と、電動車両Ｖの排気管から排出される排気が通流する希釈トンネル２と、電動車両Ｖ及び電動車両Ｖの近傍に取り付けられた複数のセンサの検出信号やシャシダイナモメータシステム１から送信される信号等に基づいて各種データを生成する計測器群３と、計測器群３によって生成されたデータを用いて各種演算を行うデータ解析装置５と、を備える。

シャシダイナモメータシステム１は、その外周面において電動車両Ｖの駆動輪と接するように電動車両Ｖの下方に設けられた複数のローラ１１と、これらローラ１１に接続された複数のダイナモメータ１２と、ローラ１１上を走行する電動車両Ｖに適切な負荷が付与されるようにダイナモメータ１２を制御する制御装置１３と、を備える。

電動車両Ｖの運転者は、図示しないドライバーズエイドに表示される運転パターン（例えば、車速パターン）に従って電動車両Ｖを運転する。

計測器群３は、走行中の電動車両Ｖの状態と相関のある様々な物理量を計測する複数の計測器３１，３２，３３，３４，３５と、これら計測器３１～３５によって計測されるデータを同期させるための基準信号を生成する基準信号生成装置３８と、を備える。これら計測器３１～３５は、それぞれ、少なくとも計測用チャンネル及び基準信号用チャンネルを備える多チャンネルである。基準信号生成装置３８は、各計測器３１～３５の基準信号用チャンネルに接続されている。基準信号生成装置３８において生成した基準信号は、各計測器３１～３５の基準信号用チャンネルに入力される。

第１計測器３１は、主に電動車両Ｖの環境情報（例えば、室温、湿度、風速、及び大気圧等）を計測するための装置である。第１計測器３１の計測用チャンネルには、上記環境情報を計測するために必要な各種センサが接続されている。第１計測器３１は、基準信号用チャンネルに入力される基準信号の時系列データである第１基準データと、計測用チャンネルへの入力信号に基づいて生成した環境情報の時系列データでありかつ第１基準データと同期する第１計測データと、を含む第１データファイル３１ｄを出力する。

第２計測器３２は、主に電動車両Ｖ及び外部電源の電力情報（例えば、電流、電圧、及び電力量等）を計測するための装置である。第２計測器３２の計測用チャンネルには、上記電力情報を計測するために必要な各種センサが接続されている。第２計測器３２は、基準信号用チャンネルに入力される基準信号の時系列データである第２基準データと、計測用チャンネルへの入力信号に基づいて生成した電力情報の時系列データでありかつ第２基準データと同期する第２計測データと、を含む第２データファイル３２ｄを出力する。

第３計測器３３は、主に電動車両Ｖの走行情報（例えば、速度や走行距離等）や排ガス情報を計測するための装置である。第３計測器３３の計測用チャンネルには、ダイナモメータシステム１の制御装置１３や排ガス情報を計測するために必要な各種センサが接続されている。第３計測器３３は、基準信号用チャンネルに入力される基準信号の時系列データである第３基準データと、計測用チャンネルへの入力信号に基づいて生成した走行情報や排ガス情報の時系列データでありかつ第３基準データと同期する第３計測データと、を含む第３データファイル３３ｄを出力する。

第４計測器３４は、主に電動車両Ｖの車載コンピュータによって監視されるモニタ情報（例えば、バッテリのＳＯＣ等）を計測するための装置である。第４計測器３４の計測用チャンネルには、電動車両Ｖの車載コンピュータに接続されている。第４計測器３４は、基準信号用チャンネルに入力される基準信号の時系列データである第４基準データと、計測用チャンネルへの入力信号に基づいて生成したモニタ情報の時系列データでありかつ第４基準データと同期する第４計測データと、を含む第４データファイル３４ｄを出力する。

第５計測器３５は、主に電動車両Ｖの排ガス情報を計測するための装置である。第５計測器３５の計測用チャンネルには、希釈トンネル２に設けられたセンサが接続されている。第５計測器３５は、基準信号用チャンネルに入力される基準信号の時系列データである第５基準データと、計測用チャンネルへの入力信号に基づいて生成した排ガス情報の時系列データでありかつ第５基準データと同期する第５計測データと、を含む第５データファイル３５ｄを出力する。

基準信号生成装置３８は、オペレータによる操作に応じて、少なくとも定常区間と立ち上がり区間とを含む基準信号を生成し、各計測器３１～３５の基準信号用チャネルへ同時に入力する。本実施形態では、基準信号生成装置３８は、オペレータによる手動操作に応じて基準信号を生成する場合について説明するが、本発明はこれに限らない。基準信号生成装置３８は、ドライバーズエイドの開始信号と連動して基準信号を生成してもよい。

図２Ａ及び図２Ｂは、基準信号の一例を示す図である。図２Ａは、立ち上がり区間においてステップ状に変化する例を示し、図２Ｂは、立ち上がり区間においてランプ状に変化する例を示す。

図１に戻り、データ解析装置５は、ＣＰＵ（Central Processing Unit）等の演算処理手段、データファイル３１ｄ～３５ｄや各種プログラムを記憶するＨＤＤ（Hard Disk Drive）やＳＳＤ（Solid State Drive）等の補助記憶手段、演算処理手段がプログラムを実行する上で一時的に必要とされるデータを格納するためのＲＡＭ（Random Access Memory）といった主記憶手段、オペレータが各種操作を行うキーボードといった操作手段、及びオペレータに各種情報を表示するディスプレイといった表示手段等のハードウェアによって構成されるコンピュータである。

図３は、計測データ同期方法及び試験方法の具体的な手順を示すフローチャートである。図３のフローチャートに示す各ステップは、データ解析装置にインストールされたコンピュータプログラムを実行することによって実現される。

始めにステップＳＴ１では、データ解析装置は、第１～第５計測器から出力される第１～第５データファイルを取得し、ステップＳＴ２へ移る。上述のように各データファイルは、基準信号を同時に第１～第５計測器の各々の基準信号用チャンネルへ入力したときに、各計測器から出力される基準データと、電動車両の状態と相関のある信号を第１～第５計測器の各々の計測用チャンネルへ入力したときに、各計測器から出力される計測データと、を含む。なお以下では、図２Ｂに示すように立ち上がり区間においてランプ状に変化する信号を基準信号として用いた場合について説明する。

次にステップＳＴ２では、データ解析装置は、ステップＳＴ１で取得した第１～第５基準データに対し微分処理を施し、ステップＳＴ３へ移る。これにより第１～第５基準データの定常区間におけるオフセットが除去される。

次にステップＳＴ３では、データ解析装置は、ステップＳＴ２における微分処理を経た第１～第５基準データに対し、所定の時間幅（２Ｋ＋１）内の中央値を抽出するメディアンフィルタ処理を施すことによって、第１～第５フィルタ処理後データを生成し、ステップＳＴ４に移る。

ここでメディアンフィルタ処理とは、１次元の時系列データであるフィルタ入力ｘ（ｎ）に対し（ｎは整数）、下記式（１）に示すように、時間幅２Ｋ＋１（Ｋは、整数）内の中央値（ｘ（ｎ－Ｋ），…，ｘ（ｎ＋Ｋ）のうちＫ＋１番目に大きい値）を抽出し、１次元の時系列データであるフィルタ出力ｙ（ｎ）を出力する処理である。
ｙ（ｎ）＝
ＭＥＤ［ｘ（ｎ－Ｋ），ｘ（ｎ－Ｋ＋１），…，ｘ（ｎ），…，ｘ（ｎ＋Ｋ）］ （１）

次にステップＳＴ４では、データ解析装置は、ステップＳＴ３において所定の時間幅（２Ｋ＋１）の下で生成した第１～第５フィルタ処理後データのうち定常区間に含まれる信号値が所定の評価時間にわたり値０を中心とする正の上限値と負の下限値との間の範囲内（後述の図４参照）に収まるか否かを判定する。データ解析装置は、ステップＳＴ４における判定結果がＮＯである場合には、整数Ｋの値をカウントアップすることによってメディアンフィルタ処理における時間幅を広げた後（ステップＳＴ５）、ステップＳＴ３に戻る。またデータ解析装置は、ステップＳＴ４における判定結果がＹＥＳである場合には、ステップＳＴ６に移る。

図４は、ステップＳＴ３～ＳＴ５におけるメディアンフィルタ処理の効果を説明するための図である。図４には、図８と同様に、ステップ状の基準信号（細一点鎖線参照）を計測器の基準信号用チャンネルに入力した場合に、この計測器から出力される基準データ（細実線参照）の一例を示す。また図４には、細実線で示す基準データに上記ステップＳＴ３～ＳＴ５に示す処理を施すことによって得られるフィルタ処理後データを太破線で示す。

図４において細実線で示すように、第１～第５計測器から出力される第１～第５基準データには、高調波のノイズが含まれている。このため、ステップＳＴ３～ＳＴ５に示すように、定常区間に含まれる信号値が値０を中心とする正の上限値と負の下限値との間の所定範囲内に収まるようにメディアンフィルタ処理の時間幅（２Ｋ＋１）を調整することにより、定常区間における信号値の変動を小さくすることができる。また図４と図８とを比較して明らかなように、本発明では基準データに対しメディアンフィルタ処理を施すことによってフィルタ処理後データを生成することにより、高調波ノイズを除去しながら、元の基準データに対し遅れを生じさせることなくフィルタ処理後データを生成することができる。

図３に戻り、ステップＳＴ６では、データ解析装置は、ステップＳＴ４において算出された第１～第５フィルタ処理後データの各々から立ち上がりタイミング情報を算出し、ステップＳＴ７に移る。ここで立ち上がりタイミング情報とは、各フィルタ処理後データにおいて信号値が立ち上がるタイミングに関する情報である。

図５は、第１～第５フィルタ処理後データを同一の時間軸上にプロットしたものの一例である。なお理解を容易にするため、図５に示す例では、各フィルタ処理後データに含まれる高調波ノイズの図示を省略する。上述のように第１～第５計測器は同期して作動させていないため、各フィルタ処理後データにおける信号値の立ち上がりタイミングは不揃いである。

ステップＳＴ６では、データ解析装置は、第１、第２、第３、第４、及び第５フィルタ処理後データにおいて信号値が図５において破線で示す閾値を超える時刻である第１立ち上がり時刻ｔ１、第２立ち上がり時刻ｔ２、第３立ち上がり時刻ｔ３、第４立ち上がり時刻ｔ４、及び第５立ち上がり時刻ｔ５を立ち上がりタイミング情報として算出する。従って図５に示す例では、第２基準データを時間（ｔ２－ｔ１）だけ進めることにより、第２基準データと第１基準データとを同期させることができ、第３基準データを時間（ｔ１－ｔ３）だけ遅らせることにより、第３基準データと第１基準データとを同期させることができ、第４基準データを時間（ｔ４－ｔ１）だけ進めることにより、第４基準データと第１基準データとを同期させることができ、第５基準データを時間（ｔ１－ｔ５）だけ遅らせることにより、第５基準データと第１基準データとを同期させることができる。

図３に戻り、ステップＳＴ７では、データ解析装置は、ステップＳＴ６において算出した第１～第５立ち上がり時刻に基づいて、ステップＳＴ１において取得した第１～第５計測データを同期する。上述のように第１～第５基準データは、それぞれ第１～第５計測データと同期していることから、第１～第５立ち上がり時刻を用いることにより、これら第１～第５計測データを同期させることができる。

次にステップＳＴ８では、データ解析装置は、ステップＳＴ７において同期した第１～第５計測データを用いて所定の演算を行うことにより、電動車両の性能を評価する。

本実施形態に係る計測データ同期方法（図３のステップＳＴ１～ＳＴ７）及び試験方法（図３のステップＳＴ１～ＳＴ８）によれば、以下の効果を奏する。
（１）本実施形態では、定常区間と立ち上がり区間とを含む基準信号を同時に第１～第５計測器へ入力したときに各計測器から出力される時系列データである第１～第５基準データを取得し、これら第１～第５基準データに対しメディアンフィルタ処理を施すことによって第１～第５フィルタ処理後データを生成し、これら第１～第５フィルタ処理後データの各々から第１～第５立ち上がり時刻を算出し、これら第１～第５立ち上がり時刻に基づいて、第１～第５計測器から出力される第１～第５計測データを同期する。特に本実施形態では、各基準データに対しメディアンフィルタ処理を施すことにより、各基準データに含まれる信号の立ち上がりに遅れを生じさせることなく高周波ノイズを除去することができる。よって本実施形態によれば、メディアンフィルタ処理を経て生成された第１～第５フィルタ処理後データを用いることにより、精度良く第１～第５立ち上がり時刻を取得できるので、第１～第５計測器から出力される第１～第５計測データを精度良く同期させることができる。

（２）本実施形態では、第１～第５フィルタ処理後データのうち定常区間に含まれる信号値が評価時間にわたり所定範囲内に収まるように、メディアンフィルタ処理における時間幅を調整することにより、信号値の立ち上がりに遅れを生じさせることなく、インパルス性ノイズのような突発的なノイズを適切に除去することができるので、結果として第１～第５計測データを精度良く同期させることができる。

（３）本実施形態では、立ち上がり区間においてランプ状に変化する信号を基準信号として用い、第１～第５基準データに対し微分処理を施した後、メディアンフィルタ処理を施すことによって第１～第５フィルタ処理後データを生成する。これにより、第１～第５計測器の第１～第５計測データにおける立ち上がりに遅れを生じさせることなく、これら第１～第５計測データに含まれるオフセット成分を除去することができるので、結果として第１～第５計測データを精度良く同期させることができる。

（４）本実施形態では、第１～第５フィルタ処理後データにおいて信号値が所定の閾値を超える時刻を立ち上がり時刻として取得することにより、簡易な演算によって立ち上がり時刻を取得することができる。

（５）本発明では、定常区間と立ち上がり区間とを含む基準信号を同時に複数の計測器の各々の基準信号用チャンネルへ入力したときに各計測器から出力される時系列データである基準データと、供試体の状態と相関のある信号を複数の計測器の各々の計測用チャンネルへ入力したときに各計測器から出力される時系列データである計測データと、を取得する。また本発明では、複数の基準データに対しメディアンフィルタ処理を施すことによって複数のフィルタ処理後データを生成し、これらフィルタ処理後データの各々から立ち上がりタイミング情報を取得し、これら立ち上がりタイミング情報に基づいて、複数の計測器から出力される時系列データを同期し、これら同期した複数の計測データに基づいて供試体の供試体の性能を評価する。特に本発明では、各基準データに対しメディアンフィルタ処理を施すことにより、基準データに含まれる信号の立ち上がりに遅れを生じさせることなく高周波ノイズを除去することができる。よって本発明によれば、メディアンフィルタ処理を経て生成されたフィルタ処理後データを用いることにより、精度良く立ち上がりタイミング情報を取得できるので、各計測器から出力される計測データを精度良く同期させることができ、ひいては精度良く供試体の性能を評価することができる。

＜第２実施形態＞
次に、本発明の第２実施形態について、図面を参照しながら詳細に説明する。
本実施形態は、第１実施形態とデータ解析装置における処理の手順が異なる。

図６は、本実施形態に係る計測データ同期方法及び試験方法の具体的な手順を示すフローチャートである。なお図６におけるステップＳＴ１１，ＳＴ１７は、図３におけるステップＳＴ１，ＳＴ８と同じであるので、詳細な説明を省略する。なお以下では、図２Ｂに示すように立ち上がり区間においてランプ状に変化する信号を基準信号として用いた場合について説明する。

ステップＳＴ１２では、データ解析装置は、ステップＳＴ１１で取得した第１～第５基準データに対し、式（１）を参照して説明したメディアンフィルタ処理を施すことによって、第１～第５フィルタ処理後データを生成し、ステップＳＴ１３に移る。

次にステップＳＴ１３では、データ解析装置は、ステップＳＴ１２において所定の時間幅（２Ｋ＋１）の下で生成した第１～第５フィルタ処理後データのうち定常区間に含まれる信号値が所定の評価時間にわたり値０を中心とする正の上限値と負の下限値との間の範囲内（上述の図４参照）に収まるか否かを判定する。データ解析装置は、ステップＳＴ１３における判定結果がＮＯである場合には、整数Ｋの値をカウントアップすることによってメディアンフィルタ処理における時間幅を広げた後（ステップＳＴ１４）、ステップＳＴ１２に戻る。またデータ解析装置は、ステップＳＴ１３における判定結果がＹＥＳである場合には、ステップＳＴ１５に移る。すなわち本実施形態は、第１～第５基準データに対し、微分処理を施さずにメディアンフィルタ処理を施すことによって第１～第５フィルタ処理後データを生成する点において、第１実施形態と異なる。

次にステップＳＴ１５では、データ解析装置は、ステップＳＴ１３において算出された第１～第５フィルタ処理後データの各々から立ち上がりタイミング情報を算出し、ステップＳＴ１６に移る。より具体的には、データ解析装置は、第１～第５フィルタ処理後データから、各々の所定の基準に対する立ち上がり箇所の遅れ時間を立ち上がりタイミング情報として算出する。なお以下では、第１～第５フィルタ処理後データにおける立ち上がり箇所の遅れ時間を、それぞれ第１～第５遅れ時間という。

図７は、第１及び第２フィルタ処理後データの立ち上がり箇所の遅れ時間Δｔ１，Δｔ２を算出する手順を説明するための図である。なお以下では、第１フィルタ処理後データを基準とした場合、すなわち第１遅れ時間Δｔ１を０とした場合について説明する。

まず、データ解析装置は、各フィルタ処理後データから、立ち上がり区間内の信号値（図７の例では、０よりもやや大きな値に定められた下限値とこの下限値よりも大きな値に定められた上限値との間に含まれる信号値）を抽出する。次にデータ解析装置は、各フィルタ処理後データの立ち上がり区間内の関数形を２つの係数（傾き及び切片）によって特徴付けられる一次関数によって近似するとともに、抽出した複数の信号値を用いて最小二乗法によって各一次関数の傾き及び切片を算出する。なお下記式（２－１）及び（２－２）には、第１及び第２フィルタ処理後データの立ち上がり区間の信号値（Ｙ１，Ｙ２）を時間（Ｘ１ｉ，Ｘ２ｉ）の一次関数として近似した例を示す。なお下記式（２－１）及び（２－２）において係数（ａ１，ａ２）は傾きであり、係数（ｂ１，ｂ２）は切片である。
Ｙ１＝ａ１・Ｘ１ｉ＋ｂ１ （２－１）
Ｙ２＝ａ２・Ｘ２ｉ＋ｂ２ （２－２）

次にデータ解析装置は、各フィルタ処理後データから算出した傾き及び切片を用いることによって、所定の基準（本実施形態では、第１フィルタ処理後データ）に対する立ち上がり箇所の遅れ時間を算出する。より具体的には、第１フィルタ処理後データを基準とした場合、第２遅れ時間Δｔ２は、先に算出した傾き（ａ１，ａ２）及び切片（ｂ１，ｂ２）を用いて下記式（３）によって算出される。なお下記式（３）では、第２フィルタ処理後データの立ち上がり区間内における一次関数の傾き（ａ２）は、第１フィルタ処理後データの立ち上がり区間内における一次関数の傾き（ａ１）と等しいと仮定したが、本発明はこれに限らない。例えば各フィルタ処理後データの立ち上がり区間内における一次関数の傾きは、上述のように最小二乗法によって算出した各フィルタ処理後データの傾きの平均値としてもよい。
Δｔ２＝（ｂ２－ｂ１）／ａ１ （３）

図６に戻り、ステップＳＴ１６では、データ解析装置は、ステップＳＴ１５において算出した第１～第５遅れ時間に基づいて、ステップＳＴ１１において取得した第１～第５計測データを同期する。上述のように第１～第５基準データは、それぞれ第１～第５計測データと同期していることから、第１～第５遅れ時間を用いることにより、これら第１～第５計測データを同期させることができる。すなわち、上述のように第１基準データを基準とした場合（第１遅れ時間を０とした場合）、第２～第５計測データを、それぞれ第２～第５遅れ時間だけ進ませることにより、これら第１～第５計測データを同期させることができる。

本実施形態に係る計測データ同期方法（図６のステップＳＴ１１～ＳＴ１６）及び試験方法（図６のステップＳＴ１１～ＳＴ１７）によれば、以下の効果を奏する。
（６）本実施形態では、立ち上がり区間においてランプ状に変化する信号を基準信号として用い、フィルタ処理後データのうち立ち上がり区間内に含まれる複数の信号値に基づいて生成される一次関数の傾き及び切片を算出し、これら傾き及び切片に基づいてフィルタ処理後データの立ち上がり箇所の遅れ時間を算出する。本実施形態によれば、このようにして算出された立ち上がり箇所の遅れ時間を用いることにより、各計測器の時系列データをさらに精度良く同期させることができる。

以上、本発明の一実施形態について説明したが、本発明はこれに限らない。本発明の趣旨の範囲内で、細部の構成を適宜変更してもよい。

例えば、上記実施形態では、基準信号生成装置３８によって生成した基準信号を各計測器３１～３５に入力した場合について説明したが、本発明はこれに限らない。例えば、車速検出信号やドライバーズエイドから出力される車速指令信号を基準信号として各計測器３１～３５に入力してもよい。

Ｓ…試験システム
Ｖ…電動車両（供試体）
１…シャシダイナモメータシステム
３…計測器群
３１…第１計測器
３１ｄ…第１データファイル
３２…第２計測器
３２ｄ…第２データファイル
３３…第３計測器
３３ｄ…第３データファイル
３４…第４計測器
３４ｄ…第４データファイル
３５…第５計測器
３５ｄ…第５データファイル
３８…基準信号生成装置
５…データ解析装置（コンピュータ）

Claims (7)

1. 複数の計測器から出力される時系列データをコンピュータによって同期する計測データ同期方法であって、
   定常区間と立ち上がり区間とを含む基準信号を同時に前記複数の計測器へ入力したときに、各計測器から出力される時系列データである基準データを取得するステップと、
   前記複数の基準データに対し所定の時間幅内の中央値を抽出するメディアンフィルタ処理を施すことにより、複数のフィルタ処理後データを生成するステップと、
   前記複数のフィルタ処理後データの各々から立ち上がりタイミング情報を取得するステップと、
   前記複数の立ち上がりタイミング情報に基づいて前記複数の計測器から出力される時系列データを同期するステップと、を備えることを特徴とする計測データ同期方法。
2. 前記フィルタ処理後データを生成するステップでは、前記フィルタ処理後データのうち前記定常区間に含まれる信号値が評価時間にわたり所定範囲内に収まるように前記時間幅を調整することを特徴とする請求項１に記載の計測データ同期方法。
3. 前記基準データを取得するステップでは、前記立ち上がり区間においてランプ状に変化する信号を前記基準信号として用い、
   前記フィルタ処理後データを生成するステップでは、前記基準データに対し微分処理を施した後、前記メディアンフィルタ処理を施すことを特徴とする請求項１又は２に記載の計測データ同期方法。
4. 前記立ち上がりタイミング情報を取得するステップでは、前記フィルタ処理後データにおいて信号値が所定の閾値を超える時刻を立ち上がり時刻として取得することを特徴とする請求項１から３の何れかに記載の計測データ同期方法。
5. 前記基準データを取得するステップでは、前記立ち上がり区間においてランプ状に変化する信号を前記基準信号として用い、
   前記立ち上がりタイミング情報を取得するステップでは、前記フィルタ処理後データのうち前記立ち上がり区間内に含まれる複数の信号値に基づいて生成される一次関数の傾き及び切片を算出し、前記傾き及び前記切片に基づいて前記フィルタ処理後データの立ち上がり箇所の遅れ時間を算出することを特徴とする請求項１又は２に記載の計測データ同期方法。
6. 複数の計測器から出力される時系列データをコンピュータによって同期させるためのコンピュータプログラムであって、
   定常区間と立ち上がり区間とを含む基準信号を同時に前記複数の計測器へ入力したときに、各計測器から出力される時系列データである基準データを取得するステップと、
   定常区間と立ち上がり区間とを含む基準信号を同時に前記複数の計測器へ入力した場合に、各々から出力される時系列データを基準データとして取得するステップと、
   前記複数の基準データに対し所定の時間幅内の中央値を抽出するメディアンフィルタ処理を施すことにより、複数のフィルタ処理後データを生成するステップと、
   前記複数のフィルタ処理後データの各々から立ち上がりタイミング情報を取得するステップと、
   前記複数の立ち上がりタイミング情報に基づいて前記複数の計測器から出力される時系列データを同期するステップと、を前記コンピュータに実行させるためのコンピュータプログラム。
7. 複数の計測器から出力される時系列データをコンピュータで処理することによって供試体の性能を評価する試験方法であって、
   定常区間と立ち上がり区間とを含む基準信号を同時に前記複数の計測器の各々の基準信号用チャンネルへ入力したときに、各計測器から出力される時系列データである基準データと、前記供試体の状態と相関のある信号を前記複数の計測器の各々の計測用チャンネルへ入力したときに、各計測器から出力される時系列データである計測データと、を取得するステップと、
   前記複数の基準データに対し所定の時間幅内の中央値を抽出するメディアンフィルタ処理を施すことにより、複数のフィルタ処理後データを生成するステップと、
   前記複数のフィルタ処理後データの各々から立ち上がりタイミング情報取得するステップと、
   前記複数の立ち上がりタイミング情報に基づいて前記複数の計測データを同期するステップと、
   前記同期した複数の計測データに基づいて前記供試体の性能を評価するステップと、を備えることを特徴とする試験方法。
   

Images