(57)【特許請求の範囲】
【請求項１】
入力時系列データに基づいて、前記入力時系列データの支配方程式の項として想定される微分項を含む項別の時系列データを作成する項別時系列データ作成部と、
前記入力時系列データの時刻ごとに各項と係数との積の線形和の正値を生成する線形和生成部と、
前記時系列データが対象とする時間区間または空間区間内の前記線形和の正値の合計値が設定値以下になるかまたは０に近づくように前記各項の係数を決定する係数決定部とを備えるデータ分析装置。
【請求項２】
前記係数決定部で決定された係数と各項との積の線形和を、前記入力時系列データの支配方程式として出力する支配方程式出力部をさらに備える請求項１に記載のデータ分析装置。
【請求項３】
前記項は、時間または空間に対する微分項の他、時間、空間または物理量のべき乗項、時間、空間または物理量の関数項、あるいは前記項同士の積をとった積項を含む請求項１に記載のデータ分析装置。
【請求項４】
前記線形和の正値は、前記線形和の自乗値または絶対値である請求項１に記載のデータ分析装置。
【請求項５】
前記係数決定部は、ガウスザイデル法に基づいて各項の係数を決定する請求項１に記載のデータ分析装置。
【請求項６】
前記支配方程式を離散化して解いた解と、前記入力時系列データとの比較結果に基づいて、前記支配方程式を評価する支配方程式評価部をさらに備える請求項１に記載のデータ分析装置。
【請求項７】
入力時系列データに基づいて、電力系統の潮流方程式に係る項別の時系列データを作成する項別時系列データ作成部と、
前記入力時系列データの時刻ごとに各項と係数との積の線形和の正値を生成する線形和生成部と、
前記時系列データが対象とする時間区間または空間区間内の前記線形和の正値の合計値が、設定値以下になるかまたは０に近づくように前記各項の係数を決定する係数決定部とを備える電力潮流解析装置。
【請求項８】
前記潮流方程式に係る各項の係数を前記電力系統の各線路のアドミタンスとして、前記電力系統の各線路のインピーダンスマップを作成する請求項７に記載の電力潮流解析装置。
【請求項９】
入力時系列データに基づいて、前記入力時系列データの支配方程式の項として想定される微分項を含む項別の時系列データを作成し、
前記入力時系列データの時刻ごとに各項と係数との積の線形和の正値を生成し、
前記時系列データが対象とする時間区間または空間区間内の前記線形和の正値の合計値が設定値以下になるかまたは０に近づくように前記各項の係数を決定する、
ことをコンピュータが実行するデータ分析方法。
【請求項１０】
前記支配方程式の項は、時間または空間に対する微分項の他、時間、空間または物理量のべき乗項、時間、空間または物理量の関数項、あるいは前記項同士の積をとった積項を含む請求項９に記載のデータ分析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、時系列データから支配方程式を求めることが可能なデータ分析装置、電力潮流解析装置およびデータ分析方法に関する。
【背景技術】
【０００２】
時系列データを近似的に関数表現する手法として、関数を時間のべき乗和として与え、時系列データと時間のべき乗和との誤差を最小化するように、時間の各べき乗項の係数を決める最小自乗法がある。この手法によれば、時系列データを近似関数として表すことが可能である。この手法は、時系列データを近似的に関数で表すことができ、分析や制御への応用で一定の役割を果たしてきた。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】Abdulrahman Baqais, Generic Algorithm for function approximation: An experimental investigation, International Journal of Artificial Intelligence and Applications (IJAIA), Vol. 7, No. 3, May 2016.
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、時系列データを関数表現するために、関数を時間のべき乗和として与える手法は、時系列データを近似関数あるいはデジタルデータとして評価するものであり、時系列データの背後にある支配方程式を求めることはできなかった。【０００５】
本発明は、上記事情に鑑みなされたものであり、その目的は、時系列データの背後にある支配方程式を求めることが可能なデータ分析装置、電力潮流解析装置およびデータ分析方法を提供することにある。
【課題を解決するための手段】
【０００６】
上記目的を達成するため、第１の観点に係るデータ分析装置は、入力時系列データに基づいて、前記入力時系列データの支配方程式の項として想定される微分項を含む項別の時系列データを作成する項別時系列データ作成部と、前記入力時系列データの時刻ごとに各項と係数との積の線形和の正値を生成する線形和生成部と、前記時系列データが対象とする時間区間または空間区間内の前記線形和の正値の合計値が設定値以下になるかまたは０に近づくように前記各項の係数を決定する係数決定部とを備える。
【発明の効果】
【０００７】
本発明によれば、時系列データの背後にある支配方程式を求めることができる。
【図面の簡単な説明】
【０００８】
【図１】図１は、第１実施形態に係るデータ分析装置の構成を示すブロック図である。
【図２】図２は、第２実施形態に係るデータ分析装置の構成を示すブロック図である。
【図３】図３は、第３実施形態に係るデータ分析装置の構成を示すブロック図である。
【図４】図４は、第４実施形態に係るデータ分析装置の構成を示すブロック図である。
【図５】図５は、第５実施形態に係る電力潮流解析装置の構成を示すブロック図である。
【図６】図６は、第６実施形態に係るデータ分析装置に適用可能なハードウェア構成を示すブロック図である。
【図７】図７は、比較例に係るデータ分析装置の構成を示すブロック図である。
【発明を実施するための形態】
【０００９】
実施形態について、図面を参照して説明する。なお、以下に説明する実施形態は特許請求の範囲に係る発明を限定するものではなく、また、実施形態の中で説明されている諸要素及びその組み合わせの全てが発明の解決手段に必須であるとは限らない。
【００１０】
図１は、第１実施形態に係るデータ分析装置の構成を示すブロック図である。
図１において、データ分析装置には、端末３Ａが設けられている。端末３Ａは、データ記憶媒体１、外部出力部２および表示部４に接続されている。端末３Ａには、時系列データ取得部１１、項別時系列データ作成部１２、線形和生成部１３、係数決定部１４および支配方程式出力部１５が設けられている。
【００１１】
時系列データ取得部１１は、時系列データ２１を取得する。この時系列データ２１は、時刻ｔごとの関数ｆ（ｔ）で表される波形２０をデジタル化したデジタルデータで表すことができる。時系列データ取得部１１は、データ記憶媒体１から過去の計測データを時系列データ２１として取得するようにしてもよいし、外部出力部２からリアルタイムの計測データを時系列データ２１として取得するようにしてもよい。
【００１２】
項別時系列データ作成部１２は、時系列データ２１に基づいて、時系列データ２１の支配方程式の項として想定される微分項を含む項別の時系列データ２２を作成する。時系列データ２２の項は、時間または空間に対する微分項の他、時間、空間または物理量のべき乗項、時間、空間または物理量の関数項、あるいは項同士の積をとった積項を含んでいてもよい。線形和生成部１３は、時系列データ２２の時刻ｔごとに各項と係数Ａ１〜Ａｍ（ｍは正の整数）との積の線形和の正値２３を生成する。なお、線形和の正値２３は、線形和の自乗をとった値でもよいし、線形和の絶対値をとった値でもよい。係数決定部１４は、時系列データ２１が対象とする時間区間または空間区間内の線形和の正値２３の合計値が設定値以下になるかまたは０に近づくように各項の係数Ａ１〜Ａｍを決定する。支配方程式出力部１５は、係数決定部１４で決定された係数Ａ１〜Ａｍと各項との積の線形和を、時系列データ２１の支配方程式として出力する。
【００１３】
時系列データ取得部１１で取得された時系列データ２１は、項別時系列データ作成部１２に送られる。項別時系列データ作成部１２は、時系列データから時間に対する一次微分、２次微分、時間のべき乗、物理量のべき乗、或いは各項の積の時系列データ２２を表形式で作成する。
【００１４】
例えば、項別時系列データ作成部１２は、以下のような項を予め準備することができる。
（Ａ１）時系列の一次関数、二次関数、三次関数および四次関数
（Ａ２）関数fのべき乗式f, f2 ,f3, ・・・, f10
（Ａ３）点xのべき乗式x, x2 ,x3, ・・・, x10
（Ａ４）点xの関数：sin(Cx), cos(Cx), eCx, ・・・
（Ａ５）上記（Ａ１）〜（Ａ４）の各項の積をとった積項
なお、（Ａ４）中のＣ(定数)は、遺伝的アルゴリズムで振り、最適値を出すことができる。この時、時系列データ２１の支配方程式の項として想定される全ての項を予め準備することができる。これらの項は、ライブラリとして予め保持することができる。
【００１５】
そして、各時刻に対して以下の時系列データ２１が与えられたものとする。
時刻        ：(t0, t1, t2, t3, ・・・, tN)
時系列データ：(f0, f1, f2, f3, ・・・, fN)
ただし、Ｎは２以上の整数である。ここで、時間間隔をΔとすると、ｉ番目の時刻における一次微分は（数式１）、二次微分は（数式２）で与えることができる。
【００１６】
【数１】
0006980105.tif 000002
【００１７】
【数２】
0006980105.tif 000003
【００１８】
三次以上の微分についても同様に、時間間隔Δとデータｆｉを用いて表すことができる。
ここで、時系列データ(f0, f1, f2, f3, ・・・, fN)が一次微分、二次微分、・・に依存する微分方程式を満足するものとする。この時、線形和生成部１３は、一次微分を基準にすると、時系列データ２２の各項と係数Ａ１〜Ａｍとの積の線形和は、以下の（数式３）で与えることができる。
【００１９】
【数３】
0006980105.tif 000004
【００２０】
次に、線形和生成部１３は、例えば、（数式３）で与えられる線形和を自乗することにより線形和を正値に変換し、誤差ｅｉを算出する。誤差ｅｉは、以下の（数式４）で与えることができる。
【００２１】
【数４】
0006980105.tif 000005
【００２２】
係数決定部１４は、時間(t0, t1, t2, t3, ・・・, tN)を通じて（数式４）で与えられる誤差ｅｉの合計値を算出することにより、誤差Ｅを算出する。誤差Ｅは、以下の（数式５）で与えることができる。
【００２３】
【数５】
0006980105.tif 000006
【００２４】
そして、係数決定部１４は、（数式５）で与えられる誤差Ｅを、設定値以下になるかまたは０に近づくように誤差Ｅの各項の係数Ａ１〜Ａｍを決定する。支配方程式出力部１５は、これらの係数Ａ１〜Ａｍと各項との積の線形和を、時系列データ２１の支配方程式として出力する。
【００２５】
これにより、時系列データ２１から、その支配方程式を求めることが可能となり、時系列データ２１が表す物理現象を解釈することが可能となる。また、異なる複数の時系列データ２１が得られた時に、それらの時系列データ２１に同一の物理現象が反映されているのか、異なる物理現象が反映されているのかを解釈することができ、異なる時系列データ２１間の関連性を容易に判断することが可能となる。
【００２６】
図２は、第２実施形態に係るデータ分析装置の構成を示すブロック図である。
図２の構成では、図１の端末３Ａの代わりに端末３Ｂが設けられている。端末３Ｂには、図１の係数決定部１４としてガウスザイデル法適用部１６が設けられている。ガウスザイデル法適用部１６は、ガウスザイデル法に基づいて連立方程式を数値的に解くことにより、（数式５）の誤差Ｅの各項の係数Ａ１〜Ａｍを決定する。（数式５）の誤差Ｅを最小にするには、以下の（数式６）を満足する係数Ａ１〜Ａｍを決定すればよい。
【００２７】
【数６】
0006980105.tif 000007
【００２８】
例えば、係数Ａ１の場合、（数式６）から以下の（数式７）および（数式８）を求めることができる。
【００２９】
【数７】
0006980105.tif 000008
【００３０】
【数８】
0006980105.tif 000009
【００３１】
係数Ａ２以下についても同様に求めると、係数Ａ１〜Ａｍについての一連の連立一次方程式を得ることができる。この連立一次方程式を解くことにより、（数式５）の誤差Ｅを最小にする係数Ａ１〜Ａｍを決定することができ、時系列データ２１の背後にある支配方程式を求めることができる。
【００３２】
一般に、ｎ次連立方程式は、以下の（数式９）で与えることができる。
【００３３】
【数９】
0006980105.tif 000010
【００３４】
ただし、ａｉｊはｎ×ｎ行列である。求めるべき係数はｘｊである。（数式９）を数値的に解くには、ガウスザイデル法を用いることができる。ガウスザイデル法は収束計算であり、一般に、以下の（数式１０）の漸化式で与えることができる。
【００３５】
【数１０】
0006980105.tif 000011
【００３６】
簡単な事例として３次元連立方程式の場合は、以下の（数式１１）で与えることができる。
【００３７】
【数１１】
0006980105.tif 000012
【００３８】
この場合、各ｘ１、ｘ２、ｘ３は、以下の（数式１２）、（数式１３）および（数式１４）の漸化式で与えることができる。
【００３９】
【数１２】
0006980105.tif 000013
【００４０】
【数１３】
0006980105.tif 000014
【００４１】
【数１４】
0006980105.tif 000015
【００４２】
収束計算では、許容誤差をｘ（ｉ＋１）−ｘ（ｉ）＜１０−６などと設定し、その許容誤差を満足するまで（数式１２）、（数式１３）および（数式１４）の計算を繰り返すことができる。
【００４３】
ここで、時系列データ２１の支配方程式の係数Ａ１〜Ａｍの算出にガウスザイデル法を用いることにより、支配方程式を精度よく求めることが可能となる。
【００４４】
なお、上述した第２実施形態では、係数Ａ１〜Ａｍの算出にガウスザイデル法を用いる方法について説明したが、ガウスザイデル法以外の連立方程式を数値的に解く方法で係数Ａ１〜Ａｍを算出するようにしてもよい。
【００４５】
図３は、第３実施形態に係るデータ分析装置の構成を示すブロック図である。
図３の構成では、図１の端末３Ａの代わりに端末３Ｃが設けられている。端末３Ｃには、図１の係数決定部１４として遺伝的アルゴリズム適用部１７が設けられている。遺伝的アルゴリズム適用部１７は、係数Ａ１〜Ａｍの候補を遺伝子で表現した個体を複数用意し、適応度の高い個体を優先的に選択して組み換えおよび突然変異などの操作を繰り返しながら係数Ａ１〜Ａｍを探索する。
【００４６】
遺伝的アルゴリズムは、概ね以下の手順で行われる。
・STEP-A1：適切な目的関数を設定する。本実施形態では、（数式５）の誤差Ｅの項の数となる。
・STEP-A2：「現世代」および「次世代」と呼ばれる二つの集合を準備する。各集合にはＭ（Ｍは２以上の整数）個の個体が準備される。
・STEP-A3：「現世代」集合に、初期値として乱数にてＭ個の個体を与える。
・STEP-A4：「現世代」の各個体の目的関数に対する適応度を評価する。
・STEP-A5：以下の（Ｂ１）〜（Ｂ３）の操作をある確率で行い、その結果を「次世代」として保存する。
（Ｂ１）複製：現状の個体の遺伝子をそのまま用いる。
（Ｂ２）交叉：二つの個体を選択し、遺伝子の組み換えを行う。
（Ｂ３）突然変異：選択した一つの個体の遺伝子に変化を与える。
・STEP-A6：「次世代」集合の個体数がＭ個生成されるまで、（Ｂ１）〜（Ｂ３）の操作を繰り返す。個体数がＭ個になった時点で「次世代」を新しい「現世代」として入れ替える。
・STEP-A7：上記の「現世代」→「次世代」→「現世代」→・・・の操作を最大世帯数設定値回数まで繰り返す。最終的な「現世代」の中で、目的関数に対する適応度が最も高いものを解とする。
【００４７】
第１実施形態では、項の係数Ａ１〜Ａｍを相対的に比較し、係数Ａ１〜Ａｍの小さいものを消去していた。一方、この遺伝的アルゴリズムを用いると、（数式５）で与えられる誤差Ｅを最小とし、支配方程式中の項の個数を最小とするように目的関数を設定することにより、誤差Ｅを最小化するととともに、相対的に係数Ａ１〜Ａｍの小さい項を自動的に削除することができる。これにより、時系列データ２１から求められる支配方程式に含まれる項の数を最小化でき、よりシンプルな解を求めることができる。
【００４８】
なお、項の係数Ａ１〜Ａｍのうち、その値が他の係数に比べて著しく小さい場合（例えば、１／１００以下の場合）、遺伝的アルゴリズム以外の最適化商法により、その係数の係る項を省略するようにしてもよい。
【００４９】
第１実施形態から第３実施形態を組み合わせた場合、以下のアルゴリズムにより時系列データ２１の背後にある支配方程式を求めることができる。
（Ｃ１）時系列データ：時間：(t0, t1, t2, t3, ・・・, tN)、データ：(f0, f1, f2, f3, ・・・, fN)を与える。
（Ｃ２）時系列データから一次微分、二次微分、三次微分、・・・、f、f2など、時系列データの支配方程式に想定される項を作成する。
（Ｃ３）遺伝的アルゴリズムにより、（Ｃ２）から想定される項の組み合わせを作成する。
（Ｃ４）各項に係数Ａ１〜Ａｍを割り振り、ガウスザイデル法で係数Ａ１〜Ａｍを求める。
（Ｃ５）遺伝的アルゴリズムでトータル誤差が最小になる項の組合せを見つける。
（Ｃ６）求まった遺伝的アルゴリズムの結果を元に、微分方程式の数値計算をし、確認をする。
なお、遺伝的アルゴリズムは、多目的とし、目的関数を「誤差最小」および「項数最小」とする。
【００５０】
図４は、第４実施形態に係るデータ分析装置の構成を示すブロック図である。
図４の構成では、図１の端末３Ａの代わりに端末３Ｄが設けられている。端末３Ｄには、支配方程式離散化部１８および支配方程式評価部１９が図１の構成に追加されている。支配方程式離散化部１８は、支配方程式出力部１５から出力された支配方程式を離散化する。支配方程式評価部１９は、支配方程式を離散化して解いた解と、時系列データ２１との比較結果に基づいて支配方程式を評価する。
以下、支配方程式を数値的に離散化する方法を具体的に説明する。一般的には、支配方程式のモデル式は、以下の（数式１５）で与えることができる。ただし、Ｃ１、Ｃ２、・・・は定数である。
【００５１】
【数１５】
0006980105.tif 000016
【００５２】
（数式１５）は、以下の（数式１６）のまとめることができる。
【００５３】
【数１６】
0006980105.tif 000017
【００５４】
本アルゴリズムでは、誤差ｅの合計を最小にするように各項の定数を決定する。この時、決定した項および定数の妥当性を検証するには、逆に微分方程式を適切な初期条件の元に解いて、時系列データ２１と比較検証する必要がある。
以下、その解き方について示す。
微分の点iにおける一次微分の離散式を、以下の（数式１７）で定義する。微分の点iにおける二次微分の離散式を、以下の（数式１８）で定義する。
【００５５】
【数１７】
0006980105.tif 000018
【００５６】
【数１８】
0006980105.tif 000019
【００５７】
三次以上の微分の離散式についても同様に求める。これを踏まえて、具体的なアルゴリズムを以下に示す。
・STEP-B1：決定した項の内の微分の最大次数を求める。
・STEP-B2：最大次数以下の項で（数式１６）の誤差ｅを表す。
例えば、微分方程式が一次微分のみの場合は、（数式１６）の誤差ｅは以下の（数式１９）となる。
【００５８】
【数１９】
0006980105.tif 000020
【００５９】
（数式１７）の一次微分の離散式を（数式１９）に代入することで、以下の（数式２０）が得られる。
【００６０】
【数２０】
0006980105.tif 000021
【００６１】
（数式２０）を変形することで、以下の（数式２１）が得られる。
【００６２】
【数２１】
0006980105.tif 000022
【００６３】
この時、初期条件が１つ与えられれば、（数式２１）を解くことができる。支配方程式評価部１９は、（数式２１）の解を時系列データ２１と比較することにより、支配方程式を評価することができる。より次数の高い微分式の場合も同様の手順で解くことができる。
【００６４】
図５は、第５実施形態に係る電力潮流解析装置の構成を示すブロック図である。
図５の構成では、図４の端末３Ｄの代わりに端末３Ｅが設けられている。端末３Ｅには、図１の項別時系列データ作成部１２として潮流解析用項別時系列データ作成部３２が設けられている。潮流解析用項別時系列データ作成部３２は、時系列データ２１に基づいて、電力系統の潮流方程式に係る項別の時系列データ２２を作成する。
【００６５】
時系列データ取得部１１は、外部出力部２を介し、電力系統の模式図１３を取得することができる。この時、時系列データ取得部１１は、電力系統の模式図１３から、ノードｉから系統に流入する電流値およびノードｉの電圧値を取得することができる。
【００６６】
電力系統の潮流方程式において、一般にノードｉにおける電力方程式は、以下の（数式２２）で与えることができる。また、ノードｉから系統に流入する電流は、以下の（数式２３）で与えることができる。ただし、ドットは複素数、バーは複素共役を表す。
【００６７】
【数２２】
0006980105.tif 000023
【００６８】
【数２３】
0006980105.tif 000024
【００６９】
ただし、Ｐｉは有効電力（ノードｉから系統に流入する向きを正）、Ｑｉは無効電力（ノードｉから系統に流入する向きを正）、Ｖｍドットは、ノードｉに接続される線路ｍの電圧、Ｙｉｍドットは、ノードｉに接続される線路ｍのアドミタンスである。（数式２２）および（数式２３）から、以下の（数式２４）を得ることができる。
【００７０】
【数２４】
0006980105.tif 000025
【００７１】
ここで、時系列データ取得部１１は、各時刻ｔにおけるノードｉから系統に流入する電流値およびノードｉの電圧値を時系列データ２１として取得する。潮流解析用項別時系列データ作成部３２は、各時刻ｔにおけるノードｉから系統に流入する電流値およびノードｉの電圧値を項別の時系列データ２２とする。
【００７２】
次に、線形和生成部１３は、（数式２４）を正値に変換し、ノードｉから系統に流入する電流値およびノードｉの電圧値を（数式２４）に代入することで誤差ｅｉを求める。次に、係数決定部１４は、線形和生成部１３で求めた誤差ｅｉの合計値が設定値以下になるかまたは０に近づくように各線路ｍのアドミタンスを決定する。
【００７３】
この第５実施形態では、第１実施形態から第４実施形態と異なり、電力系統の潮流方程式が予め与えられている。そして、支配方程式出力部１５は、係数決定部１４で算出された係数を各線路ｍのアドミッタンスの値として出力する。通常、電力の潮流計算では、アドミタンスは系統から既知の情報として与えられるが、その詳細な値を入手するのが困難であることも多い。この第５実施形態によれば、各ノードｉでの電力潮流が与えられれば、電力系統の模式図３１の各線路のアドミタンスを自動的に決定することができ、系統解析に必須な系統における各線路のインピーダンスマップを自動作成することができる。
【００７４】
図６は、第６実施形態に係るデータ分析装置に適用可能なハードウェア構成を示すブロック図である。
図６において、データ分析装置１００には、プロセッサ１０１、通信制御デバイス１０２、通信インターフェース１０３、主記憶デバイス１０４および外部記憶デバイス１０５が設けられている。プロセッサ１０１、通信制御デバイス１０２、通信インターフェース１０３、主記憶デバイス１０４および外部記憶デバイス１０５は、内部バス１０６を介して相互に接続されている。主記憶デバイス１０４および外部記憶デバイス１０５は、プロセッサ１０１からアクセス可能である。
【００７５】
また、データ分析装置１００の外部には、センサ１２０が設けられている。センサ１２０は、入出力インターフェース１０７を介して内部バス１０６に接続されている。
【００７６】
プロセッサ１０１は、データ分析装置１００全体の動作制御を司るハードウェアである。主記憶デバイス１０４は、例えば、ＳＲＡＭまたはＤＲＡＭなどの半導体メモリから構成することができる。主記憶デバイス１０４には、プロセッサ１０１が実行中のプログラムを格納したり、プロセッサ１０１がプログラムを実行するためのワークエリアを設けたりすることができる。
【００７７】
外部記憶デバイス１０５は、大容量の記憶容量を有する記憶デバイスであり、例えば、ハードディスク装置やＳＳＤ（Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ）である。外部記憶デバイス１０５は、各種プログラムの実行ファイルやプログラムの実行に用いられるデータを保持することができる。外部記憶デバイス１０５には、データ分析プログラム１０５Ａを格納することができる。データ分析プログラム１０５Ａは、データ分析装置１００にインストール可能なソフトウェアであってもよいし、データ分析装置１００にファームウェアとして組み込まれていてもよい。
【００７８】
通信制御デバイス１０２は、外部との通信を制御する機能を有するハードウェアである。通信制御デバイス１０２は、通信インターフェース１０３を介してネットワーク１０９に接続される。ネットワーク１０９は、インターネットなどのＷＡＮ（Ｗｉｄｅ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）であってもよいし、ＷｉＦｉなどのＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）であってもよいし、ＷＡＮとＬＡＮが混在していてもよい。
【００７９】
入出力インターフェース１０７は、センサ１２０から入力されるセンサデータをプロセッサ１０１が処理可能なデータ形式に変換する。入出力インターフェース１０７には、ＡＤコンバータを設けるようにしてもよい。
【００８０】
データ分析装置１００は、ネットワーク１０９またはセンサ１２０から時系列データ２１を取得し、外部記憶デバイス１０５に格納することができる。プロセッサ１０１がデータ分析プログラム１０５Ａを主記憶デバイス１０４に読み出し、データ分析プログラム１０５Ａを実行することにより、時系列データ２１の背後にある支配方程式を求めることができる。
【００８１】
この時、データ分析プログラム１０５Ａは、図１の項別時系列データ作成部１２、線形和生成部１３、係数決定部１４および支配方程式出力部１５の機能を実現することができる。
なお、データ分析プログラム１０５Ａの実行は、複数のプロセッサやコンピュータに分担させてもよい。あるいは、プロセッサ１０１は、ネットワーク１０９を介してクラウドコンピュータなどにデータ分析プログラム１０５Ａの全部または一部の実行を指示し、その実行結果を受け取るようにしてもよい。
【００８２】
図７は、比較例に係るデータ分析装置の構成を示すブロック図である。
図７の構成では、図１の端末３Ａの代わりに端末３Ｆが設けられている。端末３Ｆには、時系列データ取得部１１、べき乗関数和設定部４１、係数決定部４２および係数出力部４３が設けられている。べき乗関数和設定部４１は、時間ｔのべき乗関数和を設定する。このべき乗関数和は、べき関数による近似関数で時系列データ２１を表現する。このべき乗関数和は、係数をＡｉ、Ｃｉ、・・・とすると、以下の（数式２５）で与えることができる。
【００８３】
【数２５】
0006980105.tif 000026
【００８４】
係数決定部４２は、べき乗線形関数和と時系列データ２１の差を最小にするように係数Ａｉ、Ｃｉ、・・・を求める。この時、係数決定部４２は、遺伝的アルゴリズムを用いることができる。係数出力部４３は、係数決定部４２で求められた係数Ａｉ、Ｃｉ、・・・を出力する。このべき乗線形関数和は、べき関数による近似関数で時系列データ２１を表現するものであり、時系列データ２１の背後にある支配方程式ではない。
【符号の説明】
【００８５】
１  データ記憶媒体、２  外部出力部、３Ａ〜３Ｅ  端末、４  表示部、１１  時系列データ取得部、１２  項別時系列データ作成部、１３  線形和生成部、１４  係数決定部、１５  支配方程式出力部
【図１】
0006980105.tif 000027
【図２】
0006980105.tif 000028
【図３】
0006980105.tif 000029
【図４】
0006980105.tif 000030
【図５】
0006980105.tif 000031
【図６】
0006980105.tif 000032
【図７】
0006980105.tif 000033