(57)【要約】【課題】従来よりも確実に時系列のパラメータを推定する方法、推定したパラメータよりデータを予測する方法、パラメータ推定装置、データ予測装置、コンピュータプログラム、及び記録媒体を提供する。【解決手段】複数のデータからなる時系列を読み込み（Ｓ１）、白色雑音を生成し（Ｓ２）、白色雑音の分散σｅ２を計算し（Ｓ３）、ラグ０の白色雑音の自己相関関数をσｅ２、ラグが０以外の白色雑音の自己相関関数を０、時系列と該時系列よりも未来の白色雑音の値との相互相関関数の値を０として、計算に必要な情報を計算し（Ｓ４）、時系列を当てはめたＡＲＭＡモデルのパラメータを計算する（Ｓ５）。計算したパラメータが収束したか否かを判定し（Ｓ６）、収束してない場合は、計算したパラメータからＡＲＭＡモデルの式を用いて新たな白色雑音を生成し（Ｓ７）、パラメータが収束するまで計算を繰り返す。【選択図】　　　　図２
2004102831.tif 000002
【特許請求の範囲】
【請求項１】
複数のデータから構成される時系列、及び平均値が０で分散がσｅ２である白色雑音から、前記時系列の自己相関関数、前記白色雑音の自己相関関数、及び前記時系列と前記白色雑音との相互相関関数を用いて、前記時系列をＡＲＭＡモデルに当てはめたときのパラメータを推定する方法において、
ラグが０の前記白色雑音の自己相関関数の値をσｅ２とし、
ラグが０以外の前記白色雑音の自己相関関数の値を０とし、
前記時系列と前記白色雑音との相互相関関数のうち、前記時系列に含まれるデータと該データよりも未来の時点の白色雑音の値との相互相関関数の値を０として、
前記パラメータを計算することを特徴とするパラメータ推定方法。
【請求項２】
記憶部及び演算部を備えたコンピュータを用いて、Ｎ個のデータから構成される時系列を（ｐ，ｐ）次のＡＲＭＡモデルに当てはめたときのパラメータを推定する方法において、
Ｎ個のデータから構成される時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ｝を記憶部に記憶し、
平均値が０であるＮ個の値が時系列で並んだ白色雑音｛ｅｎ　：ｎ＝１，２，…，Ｎ｝を演算部にて生成し、
生成した白色雑音の分散σｅ２を演算部にて計算し、
前記ＡＲＭＡモデルにおけるＡＲ過程のｐ個のパラメータを行列要素とする行列ａ＝［ａ１，ａ２，…，ａｐ］Ｔ、及び前記ＡＲＭＡモデルにおけるＭＡ過程のｐ個のパラメータを行列要素とする行列ｂ＝［ｂ１，ｂ２，…，ｂｐ］Ｔを、
【数１】
2004102831.tif 000016
と演算部にて計算し、
計算した行列ａ及び行列ｂに含まれる前記ＡＲＭＡモデルのパラメータ、並びに前記時系列の値を前記ＡＲＭＡモデルの式に代入して、新たな白色雑音を演算部にて計算し、
計算した新たな白色雑音を用いて新たな行列ａ及び行列ｂを演算部にて計算し、
行列ａ及び行列ｂが収束するまで行列ａ及び行列ｂの計算を演算部にて繰り返す
ことを特徴とするパラメータ推定方法。
【請求項３】
記憶部及び演算部を備えたコンピュータを用いて、複数のデータから構成される時系列から該時系列に加わる新たなデータを予測する方法において、
時系列を記憶部に記憶し、
請求項２に記載のパラメータ推定方法を用い、前記時系列をＡＲＭＡモデルに当てはめたときのパラメータを演算部にて計算し、
計算したパラメータ、及び記憶している前記時系列の値を前記ＡＲＭＡモデルの式に代入して、前記時系列に加わる新たなデータを演算部にて計算する
ことを特徴とするデータ予測方法。
【請求項４】
Ｎ個のデータから構成される時系列を（ｐ，ｐ）次のＡＲＭＡモデルに当てはめたときのパラメータを推定する装置において、
Ｎ個のデータから構成される時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ｝を記憶する手段と、
平均値が０であるＮ個の値が時系列で並んだ白色雑音｛ｅｎ　：ｎ＝１，２，…，Ｎ｝を生成する手段と、
生成した白色雑音の分散σｅ２を計算する手段と、
前記ＡＲＭＡモデルにおけるＡＲ過程のｐ個のパラメータを行列要素とする行列ａ＝［ａ１，ａ２，…，ａｐ］Ｔ、及び前記ＡＲＭＡモデルにおけるＭＡ過程のｐ個のパラメータを行列要素とする行列ｂ＝［ｂ１，ｂ２，…，ｂｐ］Ｔを、
【数２】
2004102831.tif 000017
と計算する手段と、
計算した行列ａ及び行列ｂに含まれる前記ＡＲＭＡモデルのパラメータ、並びに前記時系列の値を前記ＡＲＭＡモデルの式に代入して、新たな白色雑音を計算する手段と、
計算した新たな白色雑音を用いて新たな行列ａ及び行列ｂを計算する手段と、
行列ａ及び行列ｂが収束するまで行列ａ及び行列ｂの計算を繰り返す手段と
を備えることを特徴とするパラメータ推定装置。
【請求項５】
複数のデータから構成される時系列から該時系列に加わる新たなデータを予測する装置において、
時系列を記憶する手段と、
請求項４に記載のパラメータ推定装置と、
該パラメータ推定装置を用いて、前記時系列から、前記時系列をＡＲＭＡモデルに当てはめたときのパラメータを計算する手段と、
計算したパラメータ、及び記憶している前記時系列の値を前記ＡＲＭＡモデルの式に代入して、前記時系列に加わる新たなデータを計算する手段と
を備えることを特徴とするデータ予測装置。
【請求項６】
コンピュータに、Ｎ個のデータから構成される時系列を（ｐ，ｐ）次のＡＲＭＡモデルに当てはめたときのパラメータを推定させるコンピュータプログラムにおいて、
コンピュータに、平均値が０であるＮ個の値が時系列で並んだ白色雑音｛ｅｎ　：ｎ＝１，２，…，Ｎ｝を生成させる手順と、
コンピュータに、生成した白色雑音の分散σｅ２を計算させる手順と、
コンピュータに、記憶している、Ｎ個のデータから構成される時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ｝を用いて、前記ＡＲＭＡモデルにおけるＡＲ過程のｐ個のパラメータを行列要素とする行列ａ＝［ａ１，ａ２，…，ａｐ］Ｔ、及び前記ＡＲＭＡモデルにおけるＭＡ過程のｐ個のパラメータを行列要素とする行列ｂ＝［ｂ１，ｂ２，…，ｂｐ］Ｔを、
【数３】
2004102831.tif 000018
と計算させる手順と、
コンピュータに、計算した行列ａ及び行列ｂに含まれる前記ＡＲＭＡモデルのパラメータ、並びに前記時系列の値を前記ＡＲＭＡモデルの式に代入して、新たな白色雑音を計算させる手順と、
コンピュータに、計算した新たな白色雑音を用いて新たな行列ａ及び行列ｂを計算させる手順と、
コンピュータに、行列ａ及び行列ｂが収束するまで行列ａ及び行列ｂの計算を繰り返させる手順と
を含むことを特徴とするコンピュータプログラム。
【請求項７】
コンピュータに、複数のデータから構成される時系列から該時系列に加わる新たなデータを予測させるコンピュータプログラムにおいて、
コンピュータに、請求項６に記載のコンピュータプログラムを用いて、記憶している時系列から、該時系列をＡＲＭＡモデルに当てはめたときのパラメータを計算させる手順と、
コンピュータに、計算したパラメータ、及び記憶している前記時系列の値を前記ＡＲＭＡモデルの式に代入して、前記時系列に加わる新たなデータを計算させる手順と
を含むことを特徴とするコンピュータプログラム。
【請求項８】
コンピュータに、Ｎ個のデータから構成される時系列を（ｐ，ｐ）次のＡＲＭＡモデルに当てはめたときのパラメータを推定させるコンピュータプログラムを記録してあるコンピュータでの読み取りが可能な記録媒体において、
コンピュータに、平均値が０であるＮ個の値が時系列で並んだ白色雑音｛ｅｎ　：ｎ＝１，２，…，Ｎ｝を生成させる手順と、
コンピュータに、生成した白色雑音の分散σｅ２を計算させる手順と、
コンピュータに、記憶している、Ｎ個のデータから構成される時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ｝を用いて、前記ＡＲＭＡモデルにおけるＡＲ過程のｐ個のパラメータを行列要素とする行列ａ＝［ａ１，ａ２，…，ａｐ］Ｔ、及び前記ＡＲＭＡモデルにおけるＭＡ過程のｐ個のパラメータを行列要素とする行列ｂ＝［ｂ１，ｂ２，…，ｂｐ］Ｔを、
【数４】
2004102831.tif 000019
と計算させる手順と、
コンピュータに、計算した行列ａ及び行列ｂに含まれる前記ＡＲＭＡモデルのパラメータ、並びに前記時系列の値を前記ＡＲＭＡモデルの式に代入して、新たな白色雑音を計算させる手順と、
コンピュータに、計算した新たな白色雑音を用いて新たな行列ａ及び行列ｂを計算させる手順と、
コンピュータに、行列ａ及び行列ｂが収束するまで行列ａ及び行列ｂの計算を繰り返させる手順と
を含むコンピュータプログラムを記録してあることを特徴とするコンピュータでの読み取りが可能な記録媒体。
【請求項９】
コンピュータに、複数のデータから構成される時系列から該時系列に加わる新たなデータを予測させるコンピュータプログラムを記録してあるコンピュータでの読み取りが可能な記録媒体において、
コンピュータに、請求項６に記載のコンピュータプログラムを用いて、記憶している時系列から、該時系列をＡＲＭＡモデルに当てはめたときのパラメータを計算させる手順と、
コンピュータに、計算したパラメータ、及び記憶している前記時系列の値を前記ＡＲＭＡモデルの式に代入して、前記時系列に加わる新たなデータを計算させる手順と
を含むコンピュータプログラムを記録してあることを特徴とするコンピュータでの読み取りが可能な記録媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、データの時系列からパラメータを推定する方法、推定されたパラメータを用いて新たなデータを予測する方法、パラメータ推定装置、データ予測方法、コンピュータをパラメータ推定装置またはデータ予測装置として実現するためのコンピュータプログラム、及び記録媒体に関する。
【０００２】
【従来の技術】
蓄積されたデータの時系列を解析し、新たなデータの値を予想する方法は、従来より様々な方法が提案されている。株価または販売量などのデータから構成される一般の経済的な時系列は、白色ノイズによって励起された系の出力とみることができる。データの時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ，…｝を、定常エルゴード的正規過程ｘ（ｔ）の標本時系列であるとしたときには、（ｐ，ｑ）次のＡＲＭＡモデルでは、｛ｘｎ　｝の要素間の関係は、以下の式で表される。
【０００３】
【数５】
2004102831.tif 000003
【０００４】
ここで、｛ｅｎ　｝は、平均値０及び分散σｅ２の正規性雑音であり、｛ａｉ　：ｉ＝１，２，…，ｐ｝及び｛ｂｊ　：ｊ＝１，２，…，ｑ｝は、
Ａ（Ｚ−１）＝１＋ａ１　Ｚ−１＋…＋ａｐ　Ｚ−ｐ
Ｂ（Ｚ−１）＝１＋ｂ１　Ｚ−１＋…＋ｂｑ　Ｚ−ｑ
とおいたときに、Ａ（Ｚ−１）及びＢ（Ｚ−１）が規約であり、更に定常条件、可逆条件および強正実条件を満たすものとする。
【０００５】
標本の時系列｛ｘｎ　｝から（１）式のパラメータ｛ａｉ　｝及び｛ｂｊ　｝を推定するための方法として、非特許文献１に開示された疑似線形回帰法がある。簡単にするためにｐ＝ｑとし、ベクトルＺｎ　及びベクトルθを次の様に定義する。
Ｚｎ　＝［−ｘｎ−１，…，−ｘｎ−ｐ，ｅｎ−１，…，ｅｎ−ｐ］Ｔ
θ＝［ａ１，ａ２，…，ａｐ，ｂ１，ｂ２，…，ｂｐ］Ｔ
但し、ｂｐ　＝０
（１）式をｘｎ　についての式に書き直したときには、（１）式は、
ｘｎ　＝θＴ　Ｚｎ　＋ｅｎ　　　…（２）
と表すことができる。Ｎ個のデータから構成される時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ｝が得られたとき、ベクトルθを時系列から最小二乗法により求めるためには、
【０００６】
【数６】
2004102831.tif 000004
【０００７】
の右辺を最小にするθを求めればよい。従って、（３）式の右辺をθで偏微分して０に等しいとすることにより、時系列から推定したベクトルθは、次の式で与えられる。
【０００８】
【数７】
2004102831.tif 000005
【０００９】
ラグがｋである｛ｘｎ　｝の自己相関関数をＲｋ　とすると、
Ｒｋ　＝Ｅ［ｘｎ　ｘｎ＋ｋ　］
である。ここで、Ｅ［ｙｎ　］は、データ列｛ｙｎ　：ｎ＝１，２，…｝の期待値である。また、
Ｒ−ｋ＝Ｒｋ
である。ラグがｋである｛ｅｎ　｝の自己相関関数をＳｋとすると、同様に、
Ｓｋ＝Ｅ［ｅｎ　ｅｎ＋ｋ　］，Ｓ−ｋ＝Ｓｋ
である。また、｛ｘｎ　｝と｛ｅｎ　｝との相互相関関数は、
Ｔｋ　＝Ｅ［ｅｎ　ｘｎ＋ｋ　］
Ｔ−ｋ＝Ｅ［ｅｎ　ｘｎ−ｋ　］
とおくと、
Ｅ［ｘｎ　ｅｎ＋ｋ　］＝Ｅ［ｅｎ　ｘｎ−ｋ　］＝Ｔ−ｋ
Ｅ［ｘｎ　ｅｎ−ｋ　］＝Ｅ［ｅｎ　ｘｎ＋ｋ　］＝Ｔｋ
となる。
【００１０】
（４）式の右辺は、｛ｘｎ　｝の自己相関関数、｛ｅｎ　｝の自己相関関数、及び｛ｘｎ　｝と｛ｅｎ　｝との相互相関関数で表すことができ、（４）式は、
【００１１】
【数８】
2004102831.tif 000006
【００１２】
となる。ベクトルａ及びｂを、
ａ＝［ａ１，ａ２，…，ａｐ］Ｔ
ｂ＝［ｂ１，ｂ２，…，ｂｐ］Ｔ
と定義し、（５）式に含まれる行列を、
【００１３】
【数９】
2004102831.tif 000007
【００１４】
と定義したときには、（５）式は、
【００１５】
【数１０】
2004102831.tif 000008
【００１６】
と表すことができる。ベクトルａ，ｂ，ｒ，ｔはｐ行１列の行列であり、行列Ｒ，Ｔ，Ｓはｐ行ｐ列の行列であるため、（７）式はベクトルａ及びｂについての連立方程式の形に書き直すことができ、
Ｒ・ａ−τ・ｂ＝−ｒ　　　　　　　　…（８）
−τＴ　・ａ＋Ｓ・ｂ＝ｔ　　　　　　…（９）
となる。（９）式より、
ｂ＝Ｓ−１（ｔ＋τＴ　・ａ）　　　　…（１０）
であり、（１０）式を（８）式へ代入することにより、
ａ＝（Ｒ−τ・Ｓ−１・τＴ　）−１（τ・Ｓ−１・ｔ−ｒ）　　　　…（１１）
が求められる。（１１）式よりベクトルａが得られ、得られたベクトルａを（１０）式に代入することにより、ベクトルｂを計算することができる。
【００１７】
疑似線形回帰法では、白色雑音｛ｅｎ　｝を乱数などを用いて生成し、生成した｛ｅｎ　｝と時系列｛ｘｎ　｝とから、（１１）式および（１０）式を用いて、（１）式のパラメータ｛ａｉ　｝及び｛ｂｊ　｝を計算し、計算した｛ａｉ　｝及び｛ｂｊ　｝から、（１）式を用いて白色雑音｛ｅｎ　｝を計算し、パラメータ｛ａｉ　｝及び｛ｂｊ　｝が収束するまで計算を繰り返し、パラメータ｛ａｉ　｝及び｛ｂｊ　｝を求める。
【００１８】
【非特許文献１】
片山徹「システム同定入門」朝倉書店　１９９４
【００１９】
【発明が解決しようとする課題】
疑似線形回帰法は、問題の性質によっては、パラメータ｛ａｉ　｝及び｛ｂｊ　｝が収束せず、パラメータの推定ができない場合があるという問題がある。また、ｐ行ｐ列の行列計算を行うため、計算時間が長いという問題があった。
【００２０】
本発明は、以上の問題を解決するためになされたものであって、その目的とするところは、計算に用いる数値のうち理論的に０となるはずの数値を０とおいて計算を行うことにより、パラメータの値の収束を容易にし、同時に計算時間の短縮を図ることができるパラメータ推定方法、パラメータ推定装置、コンピュータを該パラメータ推定装置として実現するためのコンピュータプログラム、及び記録媒体を提供することにある。
【００２１】
また、本発明の他の目的とするところは、前記パラメータ推定方法により推定したパラメータを用いてデータの推定を行うことにより、従来の手法ではデータの予測ができなかった時系列についてもデータの予測が可能となるデータ予測方法、データ予測装置、コンピュータを該データ予測装置として実現するためのコンピュータプログラム、及び記録媒体を提供することにある。
【００２２】
【課題を解決するための手段】
第１発明に係るパラメータ推定方法は、複数のデータから構成される時系列、及び平均値が０で分散がσｅ２である白色雑音から、前記時系列の自己相関関数、前記白色雑音の自己相関関数、及び前記時系列と前記白色雑音との相互相関関数を用いて、前記時系列をＡＲＭＡモデルに当てはめたときのパラメータを推定する方法において、ラグが０の前記白色雑音の自己相関関数の値をσｅ２とし、ラグが０以外の前記白色雑音の自己相関関数の値を０とし、前記時系列と前記白色雑音との相互相関関数のうち、前記時系列に含まれるデータと該データよりも未来の時点の白色雑音の値との相互相関関数の値を０として、前記パラメータを計算することを特徴とする。
【００２３】
第２発明に係るパラメータ推定方法は、記憶部及び演算部を備えたコンピュータを用いて、Ｎ個のデータから構成される時系列を（ｐ，ｐ）次のＡＲＭＡモデルに当てはめたときのパラメータを推定する方法において、Ｎ個のデータから構成される時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ｝を記憶部に記憶し、平均値が０であるＮ個の値が時系列で並んだ白色雑音｛ｅｎ　：ｎ＝１，２，…，Ｎ｝を演算部にて生成し、生成した白色雑音の分散σｅ２を演算部にて計算し、前記ＡＲＭＡモデルにおけるＡＲ過程のｐ個のパラメータを行列要素とする行列ａ＝［ａ１，ａ２，…，ａｐ］Ｔ、及び前記ＡＲＭＡモデルにおけるＭＡ過程のｐ個のパラメータを行列要素とする行列ｂ＝［ｂ１，ｂ２，…，ｂｐ］Ｔを、
【００２４】
【数１１】
2004102831.tif 000009
【００２５】
と演算部にて計算し、計算した行列ａ及び行列ｂに含まれる前記ＡＲＭＡモデルのパラメータ、並びに前記時系列の値を前記ＡＲＭＡモデルの式に代入して、新たな白色雑音を演算部にて計算し、計算した新たな白色雑音を用いて新たな行列ａ及び行列ｂを演算部にて計算し、行列ａ及び行列ｂが収束するまで行列ａ及び行列ｂの計算を演算部にて繰り返すことを特徴とする。
【００２６】
第３発明に係るデータ予測方法は、記憶部及び演算部を備えたコンピュータを用いて、複数のデータから構成される時系列から該時系列に加わる新たなデータを予測する方法において、時系列を記憶部に記憶し、第２発明に係るパラメータ推定方法を用い、前記時系列をＡＲＭＡモデルに当てはめたときのパラメータを演算部にて計算し、計算したパラメータ、及び記憶している前記時系列の値を前記ＡＲＭＡモデルの式に代入して、前記時系列に加わる新たなデータを演算部にて計算することを特徴とする。
【００２７】
第４発明に係るパラメータ推定装置は、Ｎ個のデータから構成される時系列を（ｐ，ｐ）次のＡＲＭＡモデルに当てはめたときのパラメータを推定する装置において、Ｎ個のデータから構成される時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ｝を記憶する手段と、平均値が０であるＮ個の値が時系列で並んだ白色雑音｛ｅｎ　：ｎ＝１，２，…，Ｎ｝を生成する手段と、生成した白色雑音の分散σｅ２を計算する手段と、前記ＡＲＭＡモデルにおけるＡＲ過程のｐ個のパラメータを行列要素とする行列ａ＝［ａ１，ａ２，…，ａｐ］Ｔ、及び前記ＡＲＭＡモデルにおけるＭＡ過程のｐ個のパラメータを行列要素とする行列ｂ＝［ｂ１，ｂ２，…，ｂｐ］Ｔを、
【００２８】
【数１２】
2004102831.tif 000010
【００２９】
と計算する手段と、計算した行列ａ及び行列ｂに含まれる前記ＡＲＭＡモデルのパラメータ、並びに前記時系列の値を前記ＡＲＭＡモデルの式に代入して、新たな白色雑音を計算する手段と、計算した新たな白色雑音を用いて新たな行列ａ及び行列ｂを計算する手段と、行列ａ及び行列ｂが収束するまで行列ａ及び行列ｂの計算を繰り返す手段とを備えることを特徴とする。
【００３０】
第５発明に係るデータ予測装置は、複数のデータから構成される時系列から該時系列に加わる新たなデータを予測する装置において、時系列を記憶する手段と、第４発明に係るパラメータ推定装置と、該パラメータ推定装置を用いて、前記時系列から、前記時系列をＡＲＭＡモデルに当てはめたときのパラメータを計算する手段と、計算したパラメータ、及び記憶している前記時系列の値を前記ＡＲＭＡモデルの式に代入して、前記時系列に加わる新たなデータを計算する手段とを備えることを特徴とする。
【００３１】
第６発明に係るコンピュータプログラムは、コンピュータに、Ｎ個のデータから構成される時系列を（ｐ，ｐ）次のＡＲＭＡモデルに当てはめたときのパラメータを推定させるコンピュータプログラムにおいて、コンピュータに、平均値が０であるＮ個の値が時系列で並んだ白色雑音｛ｅｎ　：ｎ＝１，２，…，Ｎ｝を生成させる手順と、コンピュータに、生成した白色雑音の分散σｅ２を計算させる手順と、コンピュータに、記憶している、Ｎ個のデータから構成される時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ｝を用いて、前記ＡＲＭＡモデルにおけるＡＲ過程のｐ個のパラメータを行列要素とする行列ａ＝［ａ１，ａ２，…，ａｐ］Ｔ、及び前記ＡＲＭＡモデルにおけるＭＡ過程のｐ個のパラメータを行列要素とする行列ｂ＝［ｂ１，ｂ２，…，ｂｐ］Ｔを、
【００３２】
【数１３】
2004102831.tif 000011
【００３３】
と計算させる手順と、コンピュータに、計算した行列ａ及び行列ｂに含まれる前記ＡＲＭＡモデルのパラメータ、並びに前記時系列の値を前記ＡＲＭＡモデルの式に代入して、新たな白色雑音を計算させる手順と、コンピュータに、計算した新たな白色雑音を用いて新たな行列ａ及び行列ｂを計算させる手順と、コンピュータに、行列ａ及び行列ｂが収束するまで行列ａ及び行列ｂの計算を繰り返させる手順とを含むことを特徴とする。
【００３４】
第７発明に係るコンピュータプログラムは、コンピュータに、複数のデータから構成される時系列から該時系列に加わる新たなデータを予測させるコンピュータプログラムにおいて、コンピュータに、第６発明に係るコンピュータプログラムを用いて、記憶している時系列から、該時系列をＡＲＭＡモデルに当てはめたときのパラメータを計算させる手順と、コンピュータに、計算したパラメータ、及び記憶している前記時系列の値を前記ＡＲＭＡモデルの式に代入して、前記時系列に加わる新たなデータを計算させる手順とを含むことを特徴とする。
【００３５】
第８発明に係るコンピュータでの読み取りが可能な記録媒体は、コンピュータに、Ｎ個のデータから構成される時系列を（ｐ，ｐ）次のＡＲＭＡモデルに当てはめたときのパラメータを推定させるコンピュータプログラムを記録してあるコンピュータでの読み取りが可能な記録媒体において、コンピュータに、平均値が０であるＮ個の値が時系列で並んだ白色雑音｛ｅｎ　：ｎ＝１，２，…，Ｎ｝を生成させる手順と、コンピュータに、生成した白色雑音の分散σｅ２を計算させる手順と、コンピュータに、記憶している、Ｎ個のデータから構成される時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ｝を用いて、前記ＡＲＭＡモデルにおけるＡＲ過程のｐ個のパラメータを行列要素とする行列ａ＝［ａ１，ａ２，…，ａｐ］Ｔ、及び前記ＡＲＭＡモデルにおけるＭＡ過程のｐ個のパラメータを行列要素とする行列ｂ＝［ｂ１，ｂ２，…，ｂｐ］Ｔを、
【００３６】
【数１４】
2004102831.tif 000012
【００３７】
と計算させる手順と、コンピュータに、計算した行列ａ及び行列ｂに含まれる前記ＡＲＭＡモデルのパラメータ、並びに前記時系列の値を前記ＡＲＭＡモデルの式に代入して、新たな白色雑音を計算させる手順と、コンピュータに、計算した新たな白色雑音を用いて新たな行列ａ及び行列ｂを計算させる手順と、コンピュータに、行列ａ及び行列ｂが収束するまで行列ａ及び行列ｂの計算を繰り返させる手順とを含むコンピュータプログラムを記録してあることを特徴とする。
【００３８】
第９発明に係るコンピュータでの読み取りが可能な記録媒体は、コンピュータに、複数のデータから構成される時系列から該時系列に加わる新たなデータを予測させるコンピュータプログラムを記録してあるコンピュータでの読み取りが可能な記録媒体において、コンピュータに、第６発明に係るコンピュータプログラムを用いて、記憶している時系列から、該時系列をＡＲＭＡモデルに当てはめたときのパラメータを計算させる手順と、コンピュータに、計算したパラメータ、及び記憶している前記時系列の値を前記ＡＲＭＡモデルの式に代入して、前記時系列に加わる新たなデータを計算させる手順とを含むコンピュータプログラムを記録してあることを特徴とする。
【００３９】
本発明者は、ＡＲＭＡモデルの式の構成から導かれるａ　ｐｒｉｏｒｉ　ｋｎｏｗｌｅｄｇｅを導入することにより、従来の疑似線形回帰法の改良を試みた。（１）式より、ＡＲＭＡモデルでは、現在のデータの値は、過去のデータの値および過去の白色雑音の値に関連づけられている。しかし、現在のデータの値は、未来の白色雑音の値とは関連性がないため、データと該データよりも未来の白色雑音との相互相関関数は、理論的に、
Ｅ［ｘｎ　ｅｎ＋ｋ　］＝０
となる。同様に、
Ｅ［ｅｎ　ｘｎ−ｋ　］＝０
である。これにより、（５）式に含まれるＴ−ｋは理論的に０とすることができる。また、白色雑音｛ｅｎ　｝の自己相関関数は、理論的に、
【００４０】
【数１５】
2004102831.tif 000013
【００４１】
となる。実際の場合は、Ｎの値が有限であるため、相互相関関数および自己相関関数を計算した値は、理論的な値とは違った値となる。理論的な値が得られる部分は、理論的な値と置き換えることにより、（５）式は、
【００４２】
【数１６】
2004102831.tif 000014
【００４３】
と表すことができる。単位行列Ｉを用いて、（１２）式は、更に、
【００４４】
【数１７】
2004102831.tif 000015
【００４５】
と表すことができる。（１３）式を、ベクトルａ及びｂについての連立方程式の形に書き直すと、
Ｒ・ａ−Ｔ・ｂ＝−ｒ　　　　　　　　…（１４）
−ＴＴ　・ａ＋σｅ２ｂ＝ｔ　　　　　　　　　　　　…（１５）
となる。（１５）式より、
ｂ＝（ｔ＋ＴＴ　・ａ）／σｅ２　　…（１６）
が得られ、（１６）式を（１４）式に代入することにより、
ａ＝｛Ｒ−（１／σｅ２）Ｔ・ＴＴ　｝−１｛（１／σｅ２）Ｔ・ｔ−ｒ｝…（１７）
が得られる。
【００４６】
第１発明においては、時系列をＡＲＭＡモデルに当てはめたときの前記ＡＲＭＡモデルのパラメータを推定する際に、白色雑音の自己相関関数、及び時系列と白色雑音との相互相関関数に理論的な値を用いることにより、確実にパラメータの値を推定することができる。
【００４７】
第２、第４、第６及び第８発明においては、白色雑音｛ｅｎ　｝を乱数などを用いて生成し、生成した白色雑音と時系列とから、（１７）式および（１６）式を用いてベクトルａ及びｂを計算し、計算したベクトルａ及びｂから、（１）式を用いて新たな白色雑音｛ｅｎ　｝を計算し、推定すべきパラメータであるベクトルａ及びｂが収束するまで計算を繰り返すことにより、従来に比べて単純化された行列を含む式を用いるため、ベクトルａ及びｂがより収束しやすくなる。
【００４８】
第３、第５、第７及び第９発明においては、従来より確実にパラメータを推定できる方法により推定したパラメータを用いてデータを予測することにより、従来ではデータの予測ができなかった時系列についてもデータを予測することができる。
【００４９】
【発明の実施の形態】
以下本発明をその実施の形態を示す図面に基づき具体的に説明する。
図１は、本発明のデータ予測装置の構成を示すブロック図である。データ予測装置１は、本発明のパラメータ推定装置の機能を兼ね備えており、コンピュータを用いて構成されている。データ予測装置１は、演算を行うＣＰＵ（演算部）１１と、演算に伴って発生する一時的な情報を記憶するＲＡＭ（記憶部）１２と、ＣＤ−ＲＯＭドライブ等の外部記憶装置１３と、ハードディスク等の内部記憶装置１４とを備えており、ＣＤ−ＲＯＭ等の本発明に係る記録媒体２から本発明に係るコンピュータプログラム２０を外部記憶装置１３にて読み取り、読み取ったコンピュータプログラム２０を内部記憶装置１４に記憶し、ＲＡＭ１２にコンピュータプログラム２０をロードし、ＣＰＵ１１は、ロードしたコンピュータプログラム２０に基づいてデータ予測装置１に必要な処理を実行する。データ予測装置１は、キーボード又はマウス等の入力装置１５と、液晶ディスプレイ又はＣＲＴディスプレイ等の出力装置１６とを備えており、データの入力を初めとするオペレータからの操作を受け付ける構成となっている。
【００５０】
なお、データ予測装置１は、図示しない通信インタフェースを備え、通信インタフェースに接続されている図示しない外部のサーバ装置から本発明に係るコンピュータプログラム２０をダウンロードし、ＣＰＵ１１にて処理を実行する形態であってもよい。
【００５１】
図２は、本発明のデータ予測装置が行うパラメータ推定の処理およびデータ予測の処理を示すフローチャートである。データ予測装置１は、オペレータの操作により入力装置１５に入力されたデータを読み込む、又は内部記憶装置１４に記憶してあるデータを読み込む等して、Ｎ個のデータから構成されるデータの時系列｛ｘｎ　：ｎ＝１，２，…，Ｎ｝をＲＡＭ１２に読み込む（Ｓ１）。ＣＰＵ１１は、コンピュータプログラム２０をＲＡＭ１２へロードし、ロードしたコンピュータプログラム２０に従って、Ｎ個の値から構成される平均値０の白色雑音｛ｅｎ　：ｎ＝１，２，…，Ｎ｝を生成する（Ｓ２）。このとき、ＣＰＵ１１は、乱数生成のルーチンを用い、時刻の値を初期値として乱数を発生させる等して、白色雑音を生成する。ＣＰＵ１１は、次に、ＲＡＭ１２にロードしたコンピュータプログラム２０に従って、生成した白色雑音｛ｅｎ　｝の分散σｅ２を計算し（Ｓ３）、データの時系列｛ｘｎ　｝及び白色雑音｛ｅｎ　｝から、（６）式で定義された行列ｒ，ｔ，Ｒと、（１３）式で定義された行列Ｔとを計算し（Ｓ４）、計算した行列ｒ，ｔ，Ｒ，Ｔ及び白色雑音の分散σｅ２を、（１７）式および（１６）式に代入して、求めるべきパラメータであるベクトルａ及びｂを計算する（Ｓ５）。
【００５２】
ＣＰＵ１１は、次に、ＲＡＭ１２にロードしたコンピュータプログラム２０に従って、計算したベクトルａ及びｂが収束したか否かを判定し（Ｓ６）、計算したベクトルａ及びｂが収束していない場合は（Ｓ６：ＮＯ）、ベクトルａ及びｂとデータの時系列｛ｘｎ　｝とを（１）式に代入して、新たな白色雑音｛ｅｎ　｝を計算し（Ｓ７）、処理をステップＳ３へ戻し、ベクトルａ及びｂが収束するまでステップＳ３〜Ｓ７の処理を繰り返す。ステップＳ６では、ＣＰＵ１１は、例えば、最新のベクトルａ，ｂと前回に計算したベクトルａ，ｂであるベクトルａ’，ｂ’とから、
Δ＝｜ａ−ａ’｜＋｜ｂ−ｂ’｜
を計算し、Δの値が所定値よりも小さい場合に、ベクトルａ及びｂが収束したと判定する。ステップＳ６においてベクトルａ及びｂが収束している場合には（Ｓ６：ＹＥＳ）、ＣＰＵ１１は、ＲＡＭ１２にロードしたコンピュータプログラム２０に従って、計算したパラメータを用いて、ＡＲＭＡモデルにより時系列｛ｘｎ　｝の予測データを計算し（Ｓ８）、処理を終了する。以上の処理により、ＡＲＭＡモデルのパラメータが推定され、時系列の新たなデータが予測できる。
【００５３】
次に、シミュレーションを行って本発明の有用性を検証した例を示す。ＡＲＡＭＡモデルのパラメータを定めて時系列を生成し、生成した時系列から本発明のパラメータ推定方法を用いてパラメータを推定するシミュレーションを行った。図３は、シミュレーションを行うためのパラメータを示した図表である。（２，２）次のＡＲＭＡモデルを用いることとし、図３に示した３種類のモデルについてシミュレーションを行った。
【００５４】
図４は、シミュレーションを行ったときの収束率を示す図表である。与えられたパラメータを用いてＮ個のデータから構成される時系列を生成し、生成した時系列を用いてパラメータの推定を行うシミュレーションを１００回行い、ベクトルａ及びｂが収束してパラメータの推定値が得られた成功率を図中に示している。従来方法を用いた場合は、モデル３の場合には確実に計算が収束するが、モデル２の場合では、Ｎが大きいときに収束せず、モデル１の場合では全く収束しない。本発明を用いた場合では、どのモデルの場合でも、３０％前後の確率で収束している。これにより、本発明は、どのような問題に対しても一定の確率でパラメータが推定できることがわかる。
【００５５】
図５は、パラメータの推定精度を示した図表である。図３に示したパラメータａ１，ａ２，ｂ１，ｂ２　と、夫々の推定値ａ１’，ａ２’，ｂ１’，ｂ２’とから、
Δａ＝｜ａ１−ａ１’｜＋｜ａ２−ａ２’｜
Δｂ＝｜ｂ１−ｂ１’｜＋｜ｂ２−ｂ２’｜
Δ＝Δａ＋Δｂ
を、収束した場合について夫々計算して平均した値を図中に示している。計算が収束してパラメータの推定値が得られた場合については、本発明と従来方法とでパラメータの推定精度は同程度となっている。
【００５６】
以上詳述した如く、本発明においては、複数のデータから構成される時系列をＡＲＭＡモデルに当てはめて該ＡＲＭＡモデルのパラメータを推定する際に、従来ではパラメータの推定ができなかった時系列についても、パラメータの推定が可能となり、パラメータが収束する確率が向上した分だけ、パラメータ推定の計算時間を短縮する事ができる。これにより、従来ではデータの予測ができなかった時系列について、データの予測が可能となり、株価または売上高などの経済的な時系列について本発明を適用して、従来よりも確実に株価または売上高などのデータの予測が可能となる。
【００５７】
なお、本実施の形態においては、本発明のデータ予測装置１は、本発明のパラメータ推定装置の機能を兼ね備えている形態を示したが、これに限るものではなく、パラメータ推定装置と、ＡＲＭＡモデルによりデータを予測するデータ予測装置とを夫々別のコンピュータにて構成し、互いに情報を入出力してパラメータの推定とデータの予測とを行う形態としても良い。
【００５８】
【発明の効果】
第１発明においては、時系列をＡＲＭＡモデルに当てはめたときの前記ＡＲＭＡモデルのパラメータを推定する際に、白色雑音の自己相関関数、及び時系列と白色雑音との相互相関関数に理論的な値を用いることにより、確実にパラメータの値を推定することができる。
【００５９】
第２、第４、第６及び第８発明においては、従来に比べて単純化された行列を含む式を用いるため、求めるべきパラメータがより収束しやすくなり、従来ではパラメータの推定ができなかった時系列についても、パラメータの推定が可能となる。また、パラメータが収束する確率が向上した分だけ、パラメータ推定の計算時間を短縮する事ができる。
【００６０】
第３、第５、第７及び第９発明においては、従来より確実にパラメータを推定できる方法により推定したパラメータを用いてデータを予測することにより、従来ではデータの予測ができなかった時系列についてもデータを予測することができる。また、株価または売上高などの経済的な時系列について本発明を適用して、従来よりも確実に株価または売上高などのデータの予測が可能となる等、本発明は優れた効果を奏する。
【図面の簡単な説明】
【図１】本発明のデータ予測装置の構成を示すブロック図である。
【図２】本発明のデータ予測装置が行うパラメータ推定の処理およびデータ予測の処理を示すフローチャートである。
【図３】シミュレーションを行うためのパラメータを示した図表である。
【図４】シミュレーションを行ったときの収束率を示す図表である。
【図５】パラメータの推定精度を示した図表である。
【符号の説明】
１　データ予測装置（パラメータ推定装置）
１１　ＣＰＵ（演算部）
１２　ＲＡＭ（記憶部）
２　記録媒体
２０　コンピュータプログラム
【図１】
2004102831.tif 000020
【図２】
2004102831.tif 000021
【図３】
2004102831.tif 000022
【図４】
2004102831.tif 000023
【図５】
2004102831.tif 000024