特許 7438086 データの系列間の差異を判定する方法、装置およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-02-15

(45)【発行日】2024-02-26

(54)【発明の名称】データの系列間の差異を判定する方法、装置およびプログラム

(51)【国際特許分類】

   G06F 18/22 20230101AFI20240216BHJP

   G06F 11/34 20060101ALI20240216BHJP

   G01M 17/007 20060101ALN20240216BHJP

【ＦＩ】

G06F18/22

G06F11/34 176

G01M17/007 Z

【請求項の数】 8

(21)【出願番号】P 2020191919

(22)【出願日】2020-11-18

(65)【公開番号】P2022080707

(43)【公開日】2022-05-30

【審査請求日】2022-12-16

(73)【特許権者】

【識別番号】000233055

【氏名又は名称】株式会社日立ソリューションズ

(74)【代理人】

【識別番号】110002572

【氏名又は名称】弁理士法人平木国際特許事務所

(72)【発明者】

【氏名】南波 優

(72)【発明者】

【氏名】齋藤 信

(72)【発明者】

【氏名】佐藤 圭太

【審査官】山本 俊介

(56)【参考文献】

【文献】特開２００６－２８８９８３（ＪＰ，Ａ）

【文献】特開２０１４－２１７６９１（ＪＰ，Ａ）

【文献】特開２０２０－０４３６５０（ＪＰ，Ａ）

【文献】国際公開第２０１４／１８４９２８（ＷＯ，Ａ１）

【文献】特開平０７－１４７３０９（ＪＰ，Ａ）

【文献】特開２０２０－１４１３１６（ＪＰ，Ａ）

【文献】@katsuki104(イノ カツキ)，時系列データのタイムラグを検知したい！！（相互相関関数），Qiita [online]，2019年02月25日，https://qiita.com/katsuki104/items/d0d2a38efb91557e031c，[2023年11月13日検索]

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｆ １８／２２

Ｇ０６Ｆ １１／３４

Ｇ０１Ｍ １７／００７

(57)【特許請求の範囲】

【請求項1】

コンピュータによってデータの系列間の差異を判定する方法であって、
前記コンピュータが、第１系列および第２系列を取得するステップと、
前記コンピュータが、前記第１系列と前記第２系列との相関係数がより大きくなるように前記第２系列をずらすことにより、第３系列を取得する、ずらしステップと、
前記コンピュータが、前記第１系列と前記第３系列との差分の絶対値に基づいて、前記第１系列と前記第３系列との差異を判定する、判定ステップと、
前記コンピュータが、判定の結果を出力する、出力ステップと、
を含み、
前記判定ステップは、
前記コンピュータが、前記差分の絶対値が所定閾値以上である区間を取得するステップと、
前記コンピュータが、各前記区間における前記差分の絶対値の平均値を算出するステップと、を含む、方法。

【請求項2】

前記判定ステップは、さらに、前記コンピュータが、各前記区間について、前記平均値に基づき、差異のグレードを決定するステップを含む、請求項１に記載の方法。

【請求項3】

前記出力ステップにおいて、前記コンピュータは、各前記区間について、当該区間を表す情報を、当該区間に係る前記平均値に基づく順位に従って出力する、請求項１に記載の方法。

【請求項4】

前記第１系列、前記第２系列および前記第３系列はいずれも時系列のデータであり、
前記ずらしステップにおいて、前記コンピュータは、前記第２系列を、時間についてずらす、請求項１に記載の方法。

【請求項5】

前記ずらしステップにおいて、前記コンピュータは、前記第２系列を、前記相関係数が最大となるようにずらす、請求項１に記載の方法。

【請求項6】

さらに、前記コンピュータが、前記判定ステップにおいて用いられる基準値の入力を、ＧＵＩを介して受け付けるステップを含む、請求項１に記載の方法。

【請求項7】

請求項１～６のいずれか一項に記載の方法を実行することにより、データの系列間の差異を判定する装置。

【請求項8】

請求項１～６のいずれか一項に記載の方法をコンピュータに実行させるプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、データの系列間の差異を判定する方法、装置およびプログラムに関する。

【背景技術】

【0002】

車両の試験等のために、データを測定して評価する方法が周知である。このような方法の例は、特許文献１に記載される。

【0003】

データの評価において、異なる測定行為で得られた複数のデータ系列間の差異の評価が行われる場合がある。このような評価に用いられるプロットデータ可視化システムの構成例を図１３に示す。従来のプロットデータ可視化システムは、入力データ処理機能および表示機能を発揮し、記憶領域に前機種のプロットデータ（データの系列）および新機種のプロットデータを格納する。

【0004】

図１４に、従来技術に係るプロットデータ可視化方法のフローチャートの例を示す。プロットデータ可視化システムは、まず入力処理機能によってユーザによる旧機種のプロットデータと新機種のプロットデータの入力を受け付け、次に表示機能によってこれらのプロットデータを表示する。

【先行技術文献】

【特許文献】

【0005】

【文献】特開２０１７－５８３５７号公報

【発明の概要】

【発明が解決しようとする課題】

【0006】

しかしながら、従来の技術では、異なる測定行為で得られた複数のデータ系列間の差異の評価が困難であるという問題があった。

【0007】

たとえば、自動車製品の評価において、前機種と新機種とで想定外の差異が発生していないかを確認するために、前機種のデータと、新機種のデータとが比較され、差異が評価される場合がある。そのような場合に、正確な評価を行うためには、前機種と新機種のデータ系列を整合させる（たとえば測定開始時刻を一致させる）ことが必要となる。

【0008】

従来の技術では、この整合作業を作業者が目視によって行っており、時間がかかるので迅速な処理が困難である。また、目視による整合作業では必ずしも適切な整合が実現されず、差異の評価精度を向上させるのが困難である。

【0009】

本発明はこのような課題を解決するためになされたものであり、データの系列間の差異を判定する方法、装置およびプログラムにおいて、複数のデータ系列間の差異の評価をより容易にするものを提供することを目的とする。

【課題を解決するための手段】

【0010】

本発明に係る方法の一例は、
データの系列間の差異を判定する方法であって、
第１系列および第２系列を取得するステップと、
前記第１系列と前記第２系列との相関係数がより大きくなるように前記第２系列をずらすことにより、第３系列を取得する、ずらしステップと、
前記第１系列と前記第３系列との差異を判定する、判定ステップと、
判定の結果を出力する、出力ステップと、
を備える。

【0011】

本発明に係る装置の一例は、上述の方法を実行することにより、データの系列間の差異を判定する。

【0012】

本発明に係るプログラムの一例は、上述の方法をコンピュータに実行させる。

【発明の効果】

【0013】

本発明に係る、データの系列間の差異を判定する方法、装置およびプログラムによれば、複数のデータ系列間の差異の評価がより容易となる。

【図面の簡単な説明】

【0014】

【図1】発明の実施例１に係る差異判定装置の構成例を示す図。

【図2】図１の差異判定装置の機能ブロック表現例を示す図。

【図3】データの系列の構成例を示す図。

【図4】図１の差異判定装置によって実行される方法を表すフローチャートの例。

【図5】図４のステップＳ１で取得される２つの系列の具体例を示す図。

【図6】図４のステップＳ１およびＳ２において用いることができるＧＵＩの画面構成例。

【図7】相関係数を算出するための式の例。

【図8】図４のステップＳ３で取得される第３系列の具体例を示す図。

【図9】差異のグレード付け処理の具体例を説明する図。

【図10】判定結果の出力例を示す図。

【図11】図４のステップＳ３におけるより具体的な処理手順を表すフローチャートの例。

【図12】図４のステップＳ５におけるより具体的な処理手順を表すフローチャートの例。

【図13】従来のプロットデータ可視化システムの構成例を示す図。

【図14】従来のプロットデータ可視化方法のフローチャートの例。

【発明を実施するための形態】

【0015】

以下、本発明の実施の形態を添付図面に基づいて説明する。
［実施例１］
図１に、本発明の実施例１に係る差異判定装置１０の構成例を示す。差異判定装置１０は、本明細書に記載される方法を実行することにより、データの系列間の差異を判定する装置である。

【0016】

差異判定装置１０は、たとえば公知のコンピュータを用いて構成することができる。差異判定装置１０は、たとえば、入力部１１と、表示部１２と、通信部１３と、ＣＰＵ１４と、記憶部１５とを備える。

【0017】

入力部１１は、たとえばキーボードおよびマウスを含む。表示部１２は、たとえば液晶ディスプレイを含む。通信部１３は、たとえばネットワークインタフェースを含む。ＣＰＵ１４は演算装置として機能する。記憶部１５はたとえば半導体メモリおよび磁気ディスク装置を含む。記憶部１５は、複数のデータの系列を記憶することができる。たとえば、第１系列のデータファイル１６と、第２系列のデータファイル１７とを記憶する。

【0018】

また、記憶部１５は図示しないプログラムを記憶してもよい。ＣＰＵ１４がこのプログラムを実行することにより、差異判定装置１０は本明細書において説明される機能を実現してもよい。すなわち、このプログラムは、本明細書に記載される方法をコンピュータに実行させるプログラムである。

【0019】

図２に、差異判定装置１０の機能ブロック表現例を示す。差異判定装置１０は、入力データ処理機能２０と、演算結果表示機能２１と、相関係数演算機能２２と、時間ズレ補正機能２３と、問題区間判定機能２４と、グレード付け機能２５と、記憶機能２６とを発揮する。記憶機能２６は、第１系列の格納領域２７と、第２系列の格納領域２８とを実現する。これらの機能は、図１に示す入力部１１、表示部１２、通信部１３、ＣＰＵ１４、および記憶部１５が、互いに協働することにより発揮されてもよい。

【0020】

図３に、データの系列の構成例を示す。この構成例は、たとえば第１系列に係る構成を示すが、第２系列も同様の形式で構成することができる。データの系列において、複数のデータ要素が互いに関連付けられる。図３の例では、系列は複数の行および６つの列を含み、各行において、時刻を表すデータ要素と、速度を表すデータ要素と、加速度を表すデータ要素と、ヨーを表すデータ要素と、ピッチを表すデータ要素と、ロールを表すデータ要素とが、互いに関連付けられている。

【0021】

なお、図３の例では６列のデータ要素を含むが、本実施例では系列は２列のデータから構成される。たとえば、図３の例は、時刻および速度からなる系列と、時刻および加速度からなる系列と、時刻およびヨーからなる系列と、時刻およびピッチからなる系列と、時刻およびロールからなる系列とを含む。

【0022】

図４に、差異判定装置１０によって実行される方法を表すフローチャートの例を示す。このフローチャートは、データの系列間の差異を判定する方法を表す。なお、差異判定装置１０の動作の説明において、差異判定装置１０内部で各処理を実行する手段を具体的に特定しない場合があるが、当業者は各処理を適宜設計することができる（たとえば画面上でデータを入力する処理は、入力部１１、表示部１２、ＣＰＵ１４および記憶部１５が協働して、入力データ処理機能２０を実現するという処理として、周知技術に基づき設計可能である）。

【0023】

図４の処理において、まず差異判定装置１０は、第１系列および第２系列を取得する（ステップＳ１）。たとえば、第１系列は前機種のプロットデータであり、第２系列は新機種のプロットデータである。本実施例では、各系列は２つのデータ要素（たとえば時刻および速度）からなる組の集合である。とくに、本実施例では、各系列は時系列のデータである。ステップＳ１の処理は、たとえば入力データ処理機能２０を介して実現される。

【0024】

図５に、ステップＳ１で取得される２つの系列の具体例を示す。横軸は時刻（たとえばｍｓ）を表し、縦軸は速度（たとえばｍ／ｓ）を表す。第１系列Ｄ１を実線で、第２系列Ｄ２を破線で表す。この例では、第１系列Ｄ１および第２系列Ｄ２の測定処理の時間的ズレの結果として、第１系列Ｄ１に対して第２系列Ｄ２の変化が見かけ上遅れており、すなわち右側にずれている。このため、これらの系列を直接的に比較することは不適切であると言える。

【0025】

図４に戻り、ステップＳ１の後に、差異判定装置１０は、判定基準の入力を受け付ける（ステップＳ２）。判定基準は、系列間に差異があるか否か、または差異がどの程度であるかを判定する際の基準となる情報である。

【0026】

図６に、ステップＳ１およびＳ２（および後述のステップＳ６）において用いることができるＧＵＩ（グラフィックユーザインタフェース）の画面構成例を示す。ＧＵＩは、第１系列を含むデータのファイル名を入力するファイル名入力欄３１と、第２系列を含むデータのファイル名を入力するファイル名入力欄３２と、相関係数を計算すべき項目名を入力する項目名入力欄３３と、判定に用いられる基準値を入力する基準値入力欄３４とを含む。基準値は、ステップＳ４に関連して後述する相関係数の閾値と、ステップＳ５に関連して後述する許容誤差とを含んでもよい。

【0027】

このように、差異判定装置１０は、ＧＵＩを介して基準値を受け付けることができるので、基準値を容易に変更することができ、様々な値を有する測定データに関する判定を行うことができる。

【0028】

また、ＧＵＩは、第１系列を図示する第１系列表示領域３５と、第２系列を図示する第２系列表示領域３６と、第１系列と第２系列との差異を図示する差異表示領域３７とを含んでもよい。

【0029】

図４に戻り、ステップＳ２の後に、差異判定装置１０は、第１系列と第２系列との相関係数がより大きくなるように第２系列をずらすことにより、第３系列を取得する（ステップＳ３、ずらしステップ）。とくに、本実施例では、相関係数が最大となるように第２系列をずらす。ここで、「系列をずらす」とは、ある列のデータ要素をすべて同じ値だけ増加または減少させることをいい、たとえば２次元プロット図上で波形を平行移動させることに対応する。本実施例では第２系列は時間についてずらされ、すなわち波形が時間軸方向に平行移動される。ここで、第２系列は時系列のデータであるので、第３系列もまた時系列のデータとなる。

【0030】

相関係数の定義および算出方法は、公知文献に基づいて当業者が適宜設計可能であるが、たとえば図７に示す式に基づいて相関係数を算出することができる。図７において、ｒは系列ｘと系列ｙとの相関係数を表す。系列ｘおよびｙはそれぞれｎ個のデータ要素ｘ_ｉおよび要素ｙ_ｉ（ただし１≦ｉ≦ｎ）を含み、ｘおよびｙの上にバー記号を付した記号はそれぞれ系列ｘおよびｙにおけるデータ要素の平均値を表す。相関係数ｒは、－１以上１以下の値となり、値が１に近いほど正の相関が高いことを示す。

【0031】

ステップＳ３の処理は、たとえば相関係数演算機能２２および時間ズレ補正機能２３を介して実現可能である。また、ステップＳ３におけるより具体的な処理手順の例は、図１１を参照して後述する。

【0032】

このように、相関係数に基づいて自動的に第２系列をずらすことができるので、作業者の目視による手作業が不要となる。また、正確な演算に基づく適切な系列の整合が実現される。とくに、相関係数が最大となるようにずらすことにより、最適な整合が実現される。

【0033】

図８に、ステップＳ３で取得される第３系列の具体例を示す。各軸は図５と同様である。第１系列Ｄ１を実線で、第３系列Ｄ３を破線で表す。図５と比較すると、第１系列Ｄ１と第２系列Ｄ２とは時間的にあまりよく整合していないが、第１系列Ｄ１と第３系列Ｄ３とは時間的によく整合しており、第１系列Ｄ１と第３系列Ｄ３とをより正確に比較することができると言える。

【0034】

ステップＳ３の後、差異判定装置１０は、ステップＳ４およびＳ５を実行することにより、第１系列Ｄ１と第３系列Ｄ３との差異を判定する（判定ステップ）。

【0035】

まず、差異判定装置１０は、第１系列Ｄ１と第３系列Ｄ３との相関係数に基づいて差異を判定する（ステップＳ４、第１判定ステップ）。差異の判定処理は、たとえば相関係数と所定閾値との比較に基づいて行われる。この閾値は、ステップＳ２で判定基準の一部として入力することができる（たとえば図６の基準値入力欄３４）。相関係数が閾値以下であれば、第１系列Ｄ１と第３系列Ｄ３とに差異があると判定され、相関係数が閾値を超えていれば、第１系列Ｄ１と第３系列Ｄ３とに差異がないと判定される。

【0036】

閾値はたとえば０．９とすることができる。この場合には、たとえば相関係数が０．９９であれば差異がないと判定され、たとえば相関係数が０．８９であれば差異があると判定される。

【0037】

次に、差異判定装置１０は、第１系列Ｄ１と第３系列Ｄ３との差分の絶対値に基づいて差異を判定する（ステップＳ５、第２判定ステップ）。ここで、差分とは、ずらされたデータ要素以外のデータ要素の差分をいい、実施例１では速度のデータ要素の差分をいう。たとえば、各時刻における第１系列Ｄ１の速度と、第３系列Ｄ３の速度とについて、差分の絶対値が計算される。

【0038】

差異の判定処理は、たとえば差分の絶対値の最大値（最大誤差）と所定閾値（許容誤差）との比較に基づいて行われる。許容誤差は、ステップＳ２で判定基準の一部として入力することができる（たとえば図６の基準値入力欄３４）。最大誤差が許容誤差未満であれば、第１系列Ｄ１と第３系列Ｄ３とに差異がないと判定され、最大誤差が許容誤差以上であれば、第１系列Ｄ１と第３系列Ｄ３とに差異があると判定される。

【0039】

図８の例では、許容誤差はたとえば０．０２とすることができる。この場合には、たとえば最大誤差が０．００５であれば差異がないと判定され、たとえば最大誤差が０．０５５であれば差異があると判定される。

【0040】

また、ステップＳ５において差異があると判定された場合には、差異判定装置１０は、差異のグレード付けを行う。この場合、差異判定装置１０は、まず差分の絶対値が許容誤差以上である区間（以下「問題区間」と呼ぶ）を取得し、各問題区間における差分の絶対値の平均値を算出する。そして、各問題区間について、この平均値に基づき、差異のグレードを決定する。これらの処理は、たとえば問題区間判定機能２４およびグレード付け機能２５を介して実現可能である。

【0041】

図９を用いて、差異のグレード付け処理の具体例を説明する。なお、説明の便宜上、図９に示す各系列の内容は図５および図８のものとは異なる。

【0042】

差異判定装置１０は、差分の絶対値が許容誤差以上である区間を取得する。図９の例では、差分の絶対値が許容誤差０．０２以上である区間として、問題区間Ｐ１～Ｐ３が取得される。問題区間Ｐ１は時刻３７秒において開始し、差分の絶対値の平均値は０．０３７である。問題区間Ｐ２は時刻５２秒において開始し、差分の絶対値の平均値は０．０２７である。問題区間Ｐ３は時刻８１秒において開始し、差分の絶対値の平均値は０．０２０である。

【0043】

各問題区間における差異のグレードは、この平均値が大きいほど高くなるように決定される。たとえば、平均値が許容誤差の２倍未満である場合には、グレードは最低のＡグレードとなり、平均値が許容誤差の２倍以上３倍未満である場合には、グレードはＡグレードより高いＢグレードとなり、平均値が許容誤差の３倍以上４倍未満である場合には、グレードはＢグレードより高いＣグレードとなり、平均値が許容誤差の４倍以上である場合には、グレードは最高のＤグレードとなる。このような基準によれば、図９の例では問題区間Ｐ１～Ｐ３いずれもグレードＡとなる。

【0044】

以上のステップＳ４およびＳ５が判定ステップを構成するが、ステップＳ４において相関係数に基づく差異がないと判定され、かつ、ステップＳ５において差分の絶対値に基づく差異もないと判定された場合には、第１系列Ｄ１と第３系列Ｄ３との差異はないと判定される。一方、ステップＳ４およびＳ５のいずれかまたは双方で差異があると判定された場合には、第１系列Ｄ１と第３系列Ｄ３との差異があると判定される。

【0045】

このように、複数の判定を組み合わせて最終的な判定を行うことにより、多面的な判定が実現でき、柔軟な判定基準の設定が可能となる。

【0046】

ステップＳ５の後、差異判定装置１０は判定の結果を出力する（ステップＳ６、出力ステップ）。ステップＳ６の処理は、たとえば演算結果表示機能２１を介して、表示によって実現される。

【0047】

図１０に、判定結果の出力例を示す。この出力例は図９の例に対応するものである。第１系列Ｄ１と第３系列Ｄ３とに差異があるので、分析結果欄に「問題あり」と表示されている。第１系列Ｄ１と第３系列Ｄ３とに差異がない場合には、分析結果欄に「問題なし」と表示される。

【0048】

また、詳細欄に、各問題区間に関する情報が出力される。この例では、１行目が問題区間Ｐ１に対応し、２行目が問題区間Ｐ２に対応し、３行目が問題区間Ｐ３に対応する。各問題区間を特定する情報として、各問題区間の開始時刻が表示される（たとえば問題区間Ｐ１については「３７秒付近」と表示される）。

【0049】

各問題区間を表す情報は、特定の順位に従って出力されてもよい。図１０の例では、各問題区間について、当該問題区間を表す情報が、当該問題区間に係る平均値に基づく順位に従って出力されている（上述のように、問題区間Ｐ１、Ｐ２、Ｐ３における平均値はそれぞれ、０．０３７、０．０２７、０．０２０である）。また、各問題区間についてグレードが出力されてもよい。

【0050】

このように、差異判定装置１０は、平均値が高くなった問題区間から順に評価結果として出力するので、とくに差異が大きい区間を容易に認識することができ、優先して対応すべき区間を迅速に決定することができる。

【0051】

また、差異判定装置１０は、判定の結果を、各系列の図示とともに、差異表示領域３７（図６）に表示してもよい。図６の例では、問題区間が着色されて他の区間と区別され、容易に認識できるようになっている。

【0052】

なお、図１０の例では、ステップＳ５の判定結果のみを示しており、ステップＳ４の判定結果については示されないが、図１０の表示内容に加えて、またはこれに代えて、ステップＳ４の判定結果を表示してもよい。たとえば、ステップＳ４において相関係数が閾値を超えている場合に、「前機種と新機種のプロットデータが大きく異なっています」というメッセージを表示してもよい。

【0053】

また、判定結果の出力において、判定処理に係る他の情報を出力してもよい。たとえば、ステップＳ３において第２系列がずらされた時間、ステップＳ４で算出された相関係数、ステップＳ５において差異があると判定されたか否か、ステップＳ５で算出された最大誤差、ステップＳ５で算出された平均値、ステップＳ５で取得された問題区間の持続時間または終了時刻、等を出力してもよい。

【0054】

また、判定結果の出力は、図６のＧＵＩにおいて行われてもよい。たとえば、図９に示すグラフが、差異表示領域３７に表示されてもよい。

【0055】

以上説明するように、本発明の実施例１に係る差異判定装置１０は、相関係数に基づいて自動的かつ適切に第２系列をずらすことができるので、作業者の目視による手作業が不要となる。また、正確な演算に基づく適切な系列の整合が実現される。このようにして、複数のデータ系列間の差異の評価をより容易にすることができる。

【0056】

図１１に、図４のステップＳ３におけるより具体的な処理手順を表すフローチャートの例を示す。
まず、差異判定装置１０は変数を初期化する（ステップＳ１１）。すなわち、最も大きい相関係数の値を０とし、最も相関係数が大きくなった時の時間ズレ補正値を０とし、第２系列をずらす時間を０とする。

【0057】

次に、差異判定装置１０は、以下のステップＳ１２～Ｓ１６を、自動ズレ補正の時間についてループして実行する。ここで、「自動ズレ補正の時間」とは、第２系列をずらす時間をいい、０から所定の最大補正範囲（時間）まで、一定の周期を単位として増加する。この最大補正範囲は、事前に決定される値（たとえば５秒）であり、たとえば図６のＧＵＩを介して入力できるように構成されてもよい（とくに個別の入力欄は図示しない）。

【0058】

ループにおいて、まず差異判定装置１０は、第２系列をずらす時間に、１周期分の時間を加算する（ステップＳ１２）。この周期は、たとえば第２系列の測定周期であってもよく、任意に指定可能な値であってもよい。たとえば図６のＧＵＩを介して周期を入力できるように構成されてもよい。

【0059】

次に、差異判定装置１０は、第２系列をずらすことにより第３系列を取得し（ステップＳ１３）、第１系列と第３系列との相関係数を算出する（ステップＳ１４）。そして、処理中のループのステップＳ１４において算出された相関係数が、最も大きい相関係数の値より大きいか否かを判定する（ステップＳ１５）。

【0060】

ステップＳ１５においてＹＥＳと判定された場合には、最も大きい相関係数の値を表す変数に、処理中のループのステップＳ１４において算出された相関係数を代入するとともに、最も相関係数が大きくなった時の時間ズレ補正値を表す変数に、処理中のループの第２系列をずらす時間を代入する（ステップＳ１６）。ステップＳ１５の判定においてＮＯと判定された場合には、ステップＳ１６は実行されない。その後、処理は次のループに入る。

【0061】

以上のような処理手順により、当業者は図４のステップＳ３を実行することができる。

【0062】

図１２に、図４のステップＳ５におけるより具体的な処理手順を表すフローチャートの例を示す。
まず、差異判定装置１０は変数を初期化する（ステップＳ２１）。図１２のステップＳ２１に示す各変数の定義は次の通りであり、［ ］記号は配列を表す。
Ｃｏｕｎｔ：ループカウント
Ｃｏｕｎｔ＿ｄｉｆｆ：許容誤差以上の差異が連続で発生している回数
Ｎｕｍｂｅｒ：許容誤差以上の差異が発生している区間（問題区間）の数
Ｉｎｓｔａｎｔ＿ｄｉｆｆ：許容誤差以上の差異の値
Ｔｏｔａｌ＿ｄｉｆｆ：処理中の問題区間中の差異の合計値
Ａｖｅｒａｇｅ＿ｄｉｆｆ［ ］：処理中の問題区間中の差異の平均値
Ｓｔａｒｔ［ ］：許容誤差以上の差異が発生し始める時間（問題区間の開始時刻）
Ｅｎｄ［ ］：許容誤差以上の差異が連続で発生している時間帯の終了時間（問題区間の終了時刻）

【0063】

次に、差異判定装置１０は、以下のステップＳ２２～Ｓ３２を、第３系列の配列の要素についてループして実行する。この配列は、ずらされたデータ要素（たとえば時刻）以外のデータ要素（たとえば速度）を配列として表したものである。以下、第１系列および第３系列の配列を、それぞれ第１系列［ｉ］および第３系列［ｉ］と表す（ただしｉはインデックスを表す整数である）。

【0064】

ループにおいて、差異判定装置１０は、まずＣｏｕｎｔの値を１だけ増加させる（ステップＳ２２）。次に、Ｉｎｓｔａｎｔ＿ｄｉｆｆに、第１系列［Ｃｏｕｎｔ］－第３系列［Ｃｏｕｎｔ］の絶対値を代入する（ステップＳ２３）。そして、Ｉｎｓｔａｎｔ＿ｄｉｆｆが許容誤差以上であるか否かを判定する（ステップＳ２４）。

【0065】

ステップＳ２４において、Ｉｎｓｔａｎｔ＿ｄｉｆｆが許容誤差以上であると判定された（すなわち問題区間である）場合には、以下に説明するステップＳ２５～Ｓ３１を実行する。そうでない場合（すなわち問題区間ではない）場合には、Ｃｏｕｎｔ＿ｄｉｆｆおよびＴｏｔａｌ＿ｄｉｆｆを０とし（ステップＳ３２）、処理は次のループに入る。

【0066】

Ｉｎｓｔａｎｔ＿ｄｉｆｆが許容誤差以上であると判定された（すなわち問題区間である）場合には、まず差異判定装置１０は、Ｃｏｕｎｔ＿ｄｉｆｆを１だけ増加させ（ステップＳ２５）、Ｔｏｔａｌ＿ｄｉｆｆにＩｎｓｔａｎｔ＿ｄｉｆｆを加算する（ステップＳ２６）。

【0067】

次に、差異判定装置１０は、Ｃｏｕｎｔ＿ｄｉｆｆが１に等しいか否かを判定する（ステップＳ２７）。Ｃｏｕｎｔ＿ｄｉｆｆが１に等しい場合（すなわち問題区間が開始された直後）には、Ｎｕｍｂｅｒを１だけ増加させ（ステップＳ２８）、それまでに検出された問題区間の数をカウントするとともに、Ｓｔａｒｔ［Ｎｕｍｂｅｒ］に第３系列［Ｃｏｕｎｔ］を代入して問題区間の開始時刻を記憶する（ステップＳ２９）。

【0068】

ステップＳ２９の後に、または、ステップＳ２７においてＣｏｕｎｔ＿ｄｉｆｆが１に等しくない場合（すなわち問題区間が継続している場合）には、Ｅｎｄ［Ｎｕｍｂｅｒ］に第３系列［Ｃｏｕｎｔ］を代入する（ステップＳ３０）。問題区間が継続している場合にはステップＳ２８およびＳ２９は実行されないので、Ｓｔａｒｔ［Ｎｕｍｂｅｒ］の値は維持され、Ｅｎｄ［Ｎｕｍｂｅｒ］のみが増加することになる。

【0069】

次に、差異判定装置１０は、Ｔｏｔａｌ＿ｄｉｆｆをＣｏｕｎｔ＿ｄｉｆｆで除算し、商をＡｖｅｒａｇｅ＿ｄｉｆｆ［Ｎｕｍｂｅｒ］に代入する（ステップＳ３１）。このようにして、問題区間における差分の絶対値の平均値が算出される。その後、処理は次のループに入る。

【0070】

以上のような処理手順により、当業者は図４のステップＳ５を実行することができる。

【0071】

上述の実施例１において、次のような変形を施すことができる。実施例１の差異判定装置１０は、２つの系列における単一のデータ種類（たとえば速度）の差異を判定するが、複数のデータ種類の差異を判定してもよい。たとえば、図３のようなデータに基づき、速度の差異と、加速度の差異とを、それぞれ判定してもよい。その場合には、結果の出力において、２通りの判定結果を整列して表示してもよい。具体例として、横軸を共通として時間を表し、２つのグラフを上下に配置してもよい。このような構成によれば、複数種類のデータ間の関係を容易に認識することができる。

【0072】

実施例１では各系列は時系列であり、時間についてずらされるが、系列は時間に関する要素を含まないものであってもよく、時間以外のデータ要素についてずらされてもよい。

【0073】

実施例１では差異判定装置１０は単一のコンピュータによって構成されるが、複数のコンピュータを含むシステムによって構成されてもよい。

【符号の説明】

【0074】

１０…差異判定装置（差異を判定する装置）
３１，３２…ファイル名入力欄
３３…項目名入力欄
３４…基準値入力欄
３５…第１系列表示領域
３６…第２系列表示領域
３７…差異表示領域
Ｄ１…第１系列
Ｄ２…第２系列
Ｄ３…第３系列

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】