(57)【要約】【課題】時系列の連続データ将来値を高精度で予測できるようにするとともに、上記予測結果の品質評価を行なうことができるようにする。【解決手段】時系列信号値を基にして遅延ベクトルを生成する遅延ベクトル生成処理と、上記遅延ベクトル生成処理によって生成した遅延ベクトルが描く再構成アトラクタ上で、現在時刻の近傍点の挙動から将来値を予測する将来値予測処理と、上記将来値予測処理によって予測した将来値の信頼性を、上記近傍点の軌道の平行度合いに基づいて評価する信頼性評価処理とを行なうようにすることにより、時系列データの法則性の存在程度に基づいて将来値を高精度に予測できるようにするとともに、上記予測した将来値の品質評価を高精度に行なうことができるようにする。【選択図】        図１
2004310500.tif 000002
【特許請求の範囲】
【請求項１】
時系列信号値を基にして遅延ベクトルを生成する遅延ベクトル生成処理と、
上記遅延ベクトル生成処理によって生成した遅延ベクトルが描く再構成アトラクタ上で、現在時刻の近傍点の挙動から将来値を予測する将来値予測処理と、
上記将来値予測処理によって予測した将来値の信頼性を、上記近傍点の軌道の平行度合いに基づいて評価する信頼性評価処理とを行なうことを特徴とする時系列連続データの将来予測方法。
【請求項２】
上記将来値予測処理においては、上記再構成アトラクタ上の近傍ベクトルの位置から内挿問題に相当する解に基づいて将来値を予測することを特徴とする請求項１に記載の時系列連続データの将来予測方法。
【請求項３】
上記信頼性評価処理においては、再構成されたアトラクタ軌道が周囲の近傍点の軌道に対して平行に走る度合いを定常性評価指標として定義するとともに、上記定常性評価指標に基づいて上記将来値予測処理によって予測した将来値の信頼性を評価するようにしたことを特徴とする請求項１または２に記載の時系列連続データの将来予測方法。
【請求項４】
過去の時系列データを正規化する時系列データ正規化処理と、
上記時系列データ正規化処理によって正規化された過去の時系列データに基づいて遅延ベクトルを生成する遅延ベクトル生成処理と、
上記遅延ベクトル生成処理によって生成された遅延ベクトルを多次元の位相空間に埋め込んで再構成アトラクタを生成する再構成アトラクタ生成処理と、
上記再構成アトラクタ生成処理によって生成された再構成アトラクタ上で現在時刻における複数の近傍点を探索する近傍点探索処理と、
上記近傍点探索処理によって探索された複数の近傍点の挙動に基づいて将来値を予測する将来値予測処理と、
上記将来値予測処理によって予測された将来値の信頼性を上記近傍点の軌道平行度合いに基づいて評価する信頼性評価処理とを有することを特徴とする時系列連続データの将来予測方法。
【請求項５】
上記時系列データは地域冷暖房システムにおける日最大冷水負荷量であり、エネルギー需給の将来予測を行なうことを特徴とする請求項１〜４の何れか１項に記載の時系列連続データの将来予測方法。
【請求項６】
時系列信号値を基にして遅延ベクトルを生成する遅延ベクトル生成手段と、
上記遅延ベクトル生成手段によって生成した遅延ベクトルが描く再構成アトラクタ上で、現在時刻の近傍点の挙動から将来値を予測する将来値予測手段と、
上記将来値予測手段によって予測した将来値の信頼性を、上記近傍点の軌道の平行度合いに基づいて評価する信頼性評価手段とを行なうことを特徴とする時系列連続データの将来予測装置。
【請求項７】
上記将来値予測手段においては、上記再構成アトラクタ上の近傍ベクトルの位置から内挿問題に相当する解に基づいて将来値を予測することを特徴とする請求項６に記載の時系列連続データの将来予測装置。
【請求項８】
上記信頼性評価手段においては、再構成されたアトラクタ軌道が周囲の近傍点の軌道に対して平行に走る度合いを定常性評価指標として定義するとともに、上記定常性評価指標に基づいて上記将来値予測手段によって予測した将来値の信頼性を評価するようにしたことを特徴とする請求項６または７に記載の時系列連続データの将来予測装置。
【請求項９】
過去の時系列データを正規化する時系列データ正規化手段と、
上記時系列データ正規化手段によって正規化された過去の時系列データに基づいて遅延ベクトルを生成する遅延ベクトル生成手段と、
上記遅延ベクトル生成手段によって生成された遅延ベクトルを多次元の位相空間に埋め込んで再構成アトラクタを生成する再構成アトラクタ生成手段と、
上記再構成アトラクタ生成手段によって生成された再構成アトラクタ上で現在時刻における複数の近傍点を探索する近傍点探索手段と、
上記近傍点探索手段によって探索された複数の近傍点の挙動に基づいて将来値を予測する将来値予測手段と、
上記将来値予測手段によって予測された将来値の信頼性を上記近傍点の軌道平行度合いに基づいて評価する再構成アトラクタ生成手段とを有することを特徴とする時系列連続データの将来予測装置。
【請求項１０】
上記時系列データは地域冷暖房システムにおける日最大冷水負荷量であり、エネルギー需給の将来予測を行なうことを特徴とする請求項６〜９の何れか１項に記載の時系列連続データの将来予測装置。
【請求項１１】
時系列信号値を基にして遅延ベクトルを生成する遅延ベクトル生成処理と、
上記遅延ベクトル生成処理によって生成した遅延ベクトルが描く再構成アトラクタ上で、現在時刻の近傍点の挙動から将来値を予測する将来値予測処理と、
上記将来値予測処理によって予測した将来値の信頼性を、上記近傍点の軌道の平行度合いに基づいて評価する信頼性評価処理とをコンピュータに実行させることを特徴とするコンピュータプログラム。
【請求項１２】
過去の時系列データを正規化する時系列データ正規化処理と、
上記時系列データ正規化処理によって正規化された過去の時系列データに基づいて遅延ベクトルを生成する遅延ベクトル生成処理と、
上記遅延ベクトル生成処理によって生成された遅延ベクトルを多次元の位相空間に埋め込んで再構成アトラクタを生成する再構成アトラクタ生成処理と、
上記再構成アトラクタ生成処理によって生成された再構成アトラクタ上で現在時刻における複数の近傍点を探索する近傍点探索処理と、
上記近傍点探索処理によって探索された複数の近傍点の挙動に基づいて将来値を予測する将来値予測処理と、
上記将来値予測処理によって予測された将来値の信頼性を上記近傍点の軌道平行度合いに基づいて評価する再構成アトラクタ生成処理とをコンピュータに実行させることを特徴とするコンピュータプログラム。
【請求項１３】
上記請求項１１または１２に記載のコンピュータプログラムを記録したことを特徴とするコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は時系列連続データの将来予測方法、装置、コンピュータプログラム及び記録媒体に関し、特に、時系列の連続データの測定時点における実測値に基づいて上記連続データの将来の将来値を高精度に予測するために用いて好適なものである。
【０００２】
【従来の技術】
従来、時系列連続データを使用して将来の値を予測する方法が種々提案されている。例えば、エネルギー需給の将来予測を行なう分野においては、地域冷暖房の集中管理システムの日最大冷水負荷量を解析対象とし、１年前の日最大冷水負荷量の時系列データを教師データとし、翌年の日最大冷水負荷量の予測を行なうようにしている。
【０００３】
従来、このような予測を行なう方法として、例えば、過去に測定された時系列データの実測値の推移をいくつものパターンに分類しておき、これらのパターンを参照することによって予測する方法や、過去に測定されたデータの実測値から最小二乗法を用いて予測する方法が用いられていた。また、統計的解析による因果関係を用いた方法、すなわち時間のずれを取り入れた多重回帰分析による予測方法も知られている。
【０００４】
しかしながら、上述したようなパターン化による予測方法では、パターン化に必然的な限界があり、また、最小二乗法を用いる予測方法や多重回帰分析による予測方法にも予測精度に限界があった。そのため、従来の予測方法においては人間が経験や勘によって予測を補助することも多く見受けられる現状があった。
【０００５】
このような不都合を型決するために、例えば、特許文献１や特許文献２において、数学的な手法を用いて、過去の時系列データから将来の時系列データを予測する方法が提案された。
【０００６】
【特許文献１】
特開平７−２８８９３１号公報
【特許文献２】
特開平１１−８５２１３号公報
【０００７】
【発明が解決しようとする課題】
上記特許文献１によれば、いわゆる「アトラクタ」と呼ばれる軌道を構築して予測計算を行なっている。また、上記特許文献２においては、時系列の連続データの測定時点及び上記測定時点における実測値に基づいて未来の時点のデータを予測している。したがって、これらの予測方法によれば、従来と比較して高精度な予測を行なうことが可能となる。
【０００８】
しかしながら、上記特許文献１及び２に記載されている予測方法を用いた場合も、予測誤差が発生することは避けられない。また、これらの文献１及び２に記載されている予測方法の場合、予測結果に対する精度がどれほどのものであるのか不明である問題があった。
【０００９】
本発明は、上述のような問題を解決するために成されたものであり、時系列の連続データ将来値を高精度で予測できるようにするとともに、上記予測結果の品質評価を行なうことができるようにすることを目的とする。
【００１０】
【課題を解決するための手段】
本発明の時系列連続データの将来予測方法は、時系列信号値を基にして遅延ベクトルを生成する遅延ベクトル生成処理と、上記遅延ベクトル生成処理によって生成した遅延ベクトルが描く再構成アトラクタ上で、現在時刻の近傍点の挙動から将来値を予測する将来値予測処理と、上記将来値予測処理によって予測した将来値の信頼性を、上記近傍点の軌道の平行度合いに基づいて評価する信頼性評価処理とを行なうことを特徴としている。
また、本発明の他の特徴とするところは、上記将来値予測処理においては、上記再構成アトラクタ上の近傍ベクトルの位置から内挿問題に相当する解に基づいて将来値を予測することを特徴としている。
また、本発明のその他の特徴とするところは、上記信頼性評価処理においては、再構成されたアトラクタ軌道が周囲の近傍点の軌道に対して平行に走る度合いを定常性評価指標として定義するとともに、上記定常性評価指標に基づいて上記将来値予測処理によって予測した将来値の信頼性を評価するようにしたことを特徴としている。
また、本発明のその他の特徴とするところは、過去の時系列データを正規化する時系列データ正規化処理と、上記時系列データ正規化処理によって正規化された過去の時系列データに基づいて遅延ベクトルを生成する遅延ベクトル生成処理と、上記遅延ベクトル生成処理によって生成された遅延ベクトルを多次元の位相空間に埋め込んで再構成アトラクタを生成する再構成アトラクタ生成処理と、上記再構成アトラクタ生成処理によって生成された再構成アトラクタ上で現在時刻における複数の近傍点を探索する近傍点探索処理と、上記近傍点探索処理によって探索された複数の近傍点の挙動に基づいて将来値を予測する将来値予測処理と、上記将来値予測処理によって予測された将来値の信頼性を上記近傍点の軌道平行度合いに基づいて評価する信頼性評価処理とを有することを特徴としている。
また、本発明のその他の特徴とするところは、上記時系列データは地域冷暖房システムにおける日最大冷水負荷量であり、エネルギー需給の将来予測を行なうことを特徴としている。
【００１１】
本発明の時系列連続データの将来予測装置は、時系列信号値を基にして遅延ベクトルを生成する遅延ベクトル生成手段と、上記遅延ベクトル生成手段によって生成した遅延ベクトルが描く再構成アトラクタ上で、現在時刻の近傍点の挙動から将来値を予測する将来値予測手段と、上記将来値予測手段によって予測した将来値の信頼性を、上記近傍点の軌道の平行度合いに基づいて評価する信頼性評価手段とを行なうことを特徴としている。
また、本発明の他の特徴とするところは、上記将来値予測手段においては、上記再構成アトラクタ上の近傍ベクトルの位置から内挿問題に相当する解に基づいて将来値を予測することを特徴としている。
また、本発明のその他の特徴とするところは、上記信頼性評価手段においては、再構成されたアトラクタ軌道が周囲の近傍点の軌道に対して平行に走る度合いを定常性評価指標として定義するとともに、上記定常性評価指標に基づいて上記将来値予測手段によって予測した将来値の信頼性を評価するようにしたことを特徴としている。
また、本発明のその他の特徴とするところは、過去の時系列データを正規化する時系列データ正規化手段と、上記時系列データ正規化手段によって正規化された過去の時系列データに基づいて遅延ベクトルを生成する遅延ベクトル生成手段と、上記遅延ベクトル生成手段によって生成された遅延ベクトルを多次元の位相空間に埋め込んで再構成アトラクタを生成する再構成アトラクタ生成手段と、上記再構成アトラクタ生成手段によって生成された再構成アトラクタ上で現在時刻における複数の近傍点を探索する近傍点探索手段と、上記近傍点探索手段によって探索された複数の近傍点の挙動に基づいて将来値を予測する将来値予測手段と、上記将来値予測手段によって予測された将来値の信頼性を上記近傍点の軌道平行度合いに基づいて評価する再構成アトラクタ生成手段とを有することを特徴としている。
また、本発明のその他の特徴とするところは、上記時系列データは地域冷暖房システムにおける日最大冷水負荷量であり、エネルギー需給の将来予測を行なうことを特徴としている。
【００１２】
本発明のコンピュータプログラムは、時系列信号値を基にして遅延ベクトルを生成する遅延ベクトル生成処理と、上記遅延ベクトル生成処理によって生成した遅延ベクトルが描く再構成アトラクタ上で、現在時刻の近傍点の挙動から将来値を予測する将来値予測処理と、上記将来値予測処理によって予測した将来値の信頼性を、上記近傍点の軌道の平行度合いに基づいて評価する信頼性評価処理とをコンピュータに実行させることを特徴としている。
また、本発明の他の特徴とするところは、過去の時系列データを正規化する時系列データ正規化処理と、上記時系列データ正規化処理によって正規化された過去の時系列データに基づいて遅延ベクトルを生成する遅延ベクトル生成処理と、上記遅延ベクトル生成処理によって生成された遅延ベクトルを多次元の位相空間に埋め込んで再構成アトラクタを生成する再構成アトラクタ生成処理と、上記再構成アトラクタ生成処理によって生成された再構成アトラクタ上で現在時刻における複数の近傍点を探索する近傍点探索処理と、上記近傍点探索処理によって探索された複数の近傍点の挙動に基づいて将来値を予測する将来値予測処理と、上記将来値予測処理によって予測された将来値の信頼性を上記近傍点の軌道平行度合いに基づいて評価する再構成アトラクタ生成処理とをコンピュータに実行させることを特徴としている。
【００１３】
本発明の記録媒体は、上記に記載のコンピュータプログラムをコンピュータ読み取り可能に記録したことを特徴としている。
【００１４】
【発明の実施の形態】
次に、添付図面を参照しながら本発明の時系列連続データの将来予測方法、装置、コンピュータプログラム及び記録媒体の実施の形態について説明する。
【００１５】
先ず、最初に本発明の時系列連続データの将来予測方法を説明する。
予測方法の概念を図１に示し、解析の手順に従って予測方法の概要を以下に説明する。
（イ）時系列データｘ（ｔ）を使って遅延ベクトルｖ（ｔ）を作成し、上記遅延ベクトルｖ（ｔ）による再構成アトラクタを構成する。
【００１６】
【数１】
2004310500.tif 000003
【００１７】
ここで、ｍは埋め込み次元、τは遅れ時間を表している。
現在時刻をｔ０とし、再構成アトラクタ上で、ｖ（ｔ０）のＮ個の近傍ベクトルｖｉＮＮ（ｔ０）を探索決定する。上記添字ｉは、ｉ番目の大きさを持つ近傍ベクトルを示し、ｉ＝１，２…、Ｎである。上記近傍ベクトルｖｉＮＮ（ｔ０）の探索には、以下に示すユークリッドノルムＬｉを採用した。
【００１８】
【数２】
2004310500.tif 000004
【００１９】
（ロ）上で決定した近傍ベクトルの再構成アトラクタ上におけるｒ時刻後の位置ｖｉＮＮ（ｔ０）より、ベクトルｖ（ｔ０）のｒ時刻後の位置を推定し、ｘ（ｔ０＋ｒ）を決定する。
（ハ）ベクトルｖ（ｔ０）のｒ時刻後の位置を推定するため、動径基底関数（Ｒａｄｉａｌ  Ｂａｓｉｓ    Ｆｕｎｃｔｉｏｎ、以下ＲＢＦ）を使い、以下のように定式化する。
近傍ベクトルｖｉＮＮ（ｔ０）を以下の（３）式で定義する。
【００２０】
【数３】
2004310500.tif 000005
【００２１】
近傍ベクトルｖｉＮＮ（ｔ０）及びｘ（ｔ０＋ｒ）の値を（４）式に代入して構成したＮ元連立方程式を解き、重み係数γｋ及びバイアスβを決定する。
【００２２】
【数４】
2004310500.tif 000006
【００２３】
ここで、動径基底関数（ＲＢＦ）φは、（５）式に示すガウス関数により定義する。ｂは動径基底関数の影響範囲の指標となる係数である。
【００２４】
【数５】
2004310500.tif 000007
【００２５】
上記（４）式より決定した重み係数γｋ及びバイアスβを使い、ｘ（ｔ０＋ｒ）の推定値ｘ´（ｔ０＋ｒ）を下記の（６）式より計算する。
【００２６】
【数６】
2004310500.tif 000008
【００２７】
次に、予測結果の品質評価の考え方にについて説明する。
本実施の形態の予測結果の品質評価法は、上述したように、動径基底関数を使って定式化した近傍ベクトルの動きに基づき将来予測を行なうもので、近傍点の動きに或る法則性が存在することを前提にしていることから決定論的予測手法と呼ばれる。
【００２８】
そして、再構成アトラクタ上の近傍ベクトルｖｉＮＮ（ｔ０＋ｒ）の位置から推定値ｖ（ｔ０＋ｒ）の位置を推定し、推定値ｘ´（ｔ０＋ｒ）を決定する、いわゆる内挿問題となっている。近傍点のｖｉＮＮ（ｔ０）から近傍ベクトルｖｉＮＮ（ｔ０＋ｒ）への動きに強い法則性が存在すれば、ｖ（ｔ０）もそれと同伴して動く確率が高くなり、推定値ｘ´（ｔ０＋ｒ）を高い精度で予測することができる。
【００２９】
一方、当手法と対極の位置にある確率論的予測手法では、時系列ｘ（ｔ）を確率変数として扱い、ＡＲモデルの例を（７）式に示すように、現在時刻ｔ０までの時系列データｘ（ｔ０−１）、ｘ（ｔ０−２）、…、ｘ（ｔ０−ｎ）を使った回帰式に基づく時間外挿により将来値ｘ´（ｔ０＋ｒ）を予測する外挿問題となっている。
【００３０】
【数７】
2004310500.tif 000009
【００３１】
ここで、ｃは定数、φｉは回帰係数、ｅｔ０は誤差項を表す。
近傍点ｖｉＮＮ（ｔ０）からｖｉＮＮ（ｔ０＋ｒ）への動きに強い法則性が存在すれば、予測点ｘ´（ｔ０＋ｒ）の周囲に近傍点ｖｉＮＮ（ｔ０＋ｒ）が適切に配置されるので、外挿問題である確率論的予測手法より内挿問題である当手法の方が、精度の高い推定が可能となる。
【００３２】
したがって、近傍点ｖｉＮＮ（ｔ０）の動きの法則性を使って、予測結果の信頼性を評価することができる。そして、近傍点ｖｉＮＮ（ｔ０）の動きに基づく法則性の定量記述に、前章で定義した定常性評価指標を使用する。定常性評価指標の定義を以下に再記述する。
【００３３】
すなわち、状態の挙動を「ｒ」の時間スケールで観測したときに、時間発展の様子が決定論的、すなわち、或る法則性に支配されて推移するように見えるということは、図２に示すように、再構成された軌道群の近接した部分が「ｒ」時刻後に近接した部分に移されることを意味する。これは、アトラクタ軌道が周囲の近傍点の軌道に対して平行に走ることを意味している。
【００３４】
この考えに基づき、状態が或る法則下で安定した定常状態を維持する度合いを表す定常性評価指標（Ｓｔｅａｄｙ  Ｓｔａｔｅ  Ｉｎｄｅｘ）を次式で定義する。
定常性評価指標（ＳＳ１）＝（ｒ時刻後に生き残った近傍点数）／（時刻ｔ０における近傍点数）・・・（８）
【００３５】
この指標は、時系列変化の法則性依存度を示し、法則性依存度が大きく定常状態に近いほど「１」に近づき、完全なランダム状態である白色ノイズでは「０」となる。法則性依存度が大きいということは、状態が或る法則下で安定した状態を維持することを表し、同数値指標が状態の安定性の評価指標にもなっている。
【００３６】
＜予測結果の実施例＞
図３（ａ）は、１年前の日最大冷水負荷量を１年間の同データの分散値で除算することにより正規化したものを示し、本解析の教師データとなるものである。また、図３（ｂ）は、図３（ａ）の翌年の最大冷水負荷量を示し、本解析の予測対象となるものであり、図３（ａ）と同様に、１年間の同データの分散値で除算して正規化している。
【００３７】
図３（ｂ）の時系列データの再構成アトラクタ上にある各時刻における近傍点を、図３（ａ）の時系列データの再構成アトラクタ上の点から決定し、上述した方法に基づいて、翌日の最大冷水負荷量を推定した。推定結果を、図３（ｂ）に示す。
図３（ｃ）には、図３（ｂ）の実績値と予測値との誤差を示した。誤差の計算は、下記（９）式を使用した。
〔予測誤差〕＝〔予測値〕−〔実績値〕        ・・・（９）
【００３８】
また、図３（ｃ）に、定常性評価指標を併せて表示した。図３（ｃ）より、定常性評価指標と予測誤差との間に対応関係があり、定常性評価指標が低下すると予測誤差が大きくなる傾向にあるのが判る。
【００３９】
図４は、図３（ｃ）の各時刻における予測誤差と定常性評価指標の関係をプロットした図である。ここでの予測誤差は、絶対値表示をしている。
図４より、定常性評価指標が「０．７」以上のときは、予測誤差の絶対値が０〜０．４の範囲内に安定的に分布しているのが判る。
【００４０】
一方、定常性評価指標が「０．７」未満になると、定常性評価指標値の低下に伴い、「０．５」以上の大きな予測誤差が発生する確率が大きくなっていることが判る。したがって、このケースの場合、定常性評価指標が「０．７」以上のときは、本手法による予測結果に対し品質保証が可能であるが、同指標が「０．７」未満になる領域は、大きな予測誤差が発生する確率が大きくなるため、予測結果に対し他の手法との併用等を行なうようにするケアが必要であるという情報を顧客に提示することが可能となる。
【００４１】
上述したように、本実施の形態の時系列連続データの将来予測方法によれば、カオス時系列解析手法をエネルギー需給予測に適用して、時系列データの法則性の存在程度を定常性評価指標で表現し、同指標の値が高位にあるときは高精度の予測結果が得られるとともに、同指標により予測精度の品質評価を行なうことができることを示した。
【００４２】
図５は、本発明の時系列連続データの将来予測方法の処理手順の一例を説明するフローチャートである。
図５に示したように、処理が開始されると、最初のステップＳ５１において時系列データの正規化処理が行われる。これは、過去の時系列データを正規化する処理であり、上記実施例の場合には１年前の日最大冷水負荷量を１年間の同データの分散値で除算する処理であり、本解析の教師データを得る処理である。
【００４３】
次に、ステップＳ５２に進み、遅延ベクトル生成処理を行なう。これは、上記ステップＳ５１における時系列データ正規化処理によって正規化された過去の時系列データｘ（ｔ）を使って遅延ベクトルｖ（ｔ）を生成する処理である。
【００４４】
次に、ステップＳ５３に進み、再構成アトラクタを生成する処理を行なう。これは、上記ステップＳ５２における遅延ベクトル生成処理によって生成された遅延ベクトルを多次元の位相空間に埋め込む処理、すなわち、多次元の位相空間上で軌道を構築する処理である。
【００４５】
次に、ステップＳ５４に進み、近傍点探索処理を行なう。これは、上記再構成アトラクタ生成処理によって生成された再構成アトラクタ上で現在時刻における複数の近傍点を探索する処理である。
【００４６】
次に、ステップＳ５５に進み、将来値予測処理を行なう。この処理は、上記近ステップＳ５４における傍点探索処理によって探索された複数の近傍点の挙動に基づいて将来値を予測する処理である。
【００４７】
次に、ステップＳ５６に進み、信頼性評価処理を行なう。これは、上記将来値予測処理によって予測された将来値の信頼性を上記近傍点の軌道平行度合いに基づいて評価する処理である。詳述しないが、上述したように信頼性が低い場合には大きな予測誤差が発生する確率が大きくなるため、予測結果に対して他の手法との併用等を行なうことを顧客に提示する。
【００４８】
次に、上述した各処理を実行可能なコンピュータシステムの一例を図６に示す。
図６において、６００はコンピュータＰＣである。ＰＣ６００は、ＣＰＵ６０１を備え、ＲＯＭ６０２またはハードディスク（ＨＤ）６１１に記憶された、あるいはフレキシブルディスクドライブ（ＦＤ）６１２より供給されるデバイス制御ソフトウェアを実行し、システムバス６０４に接続される各デバイスを総括的に制御する。
【００４９】
前記ＰＣ６００のＣＰＵ６０１，ＲＯＭ６０２またはハードディスク（ＨＤ）６１１に記憶されたプログラムにより、本実施形態の各機能手段が構成される。
【００５０】
６０３はＲＡＭで、ＣＰＵ６０１の主メモリ、ワークエリア等として機能する。６０５はキーボードコントローラ（ＫＢＣ）であり、キーボード（ＫＢ）６０９から入力される信号をシステム本体内に入力する制御を行なう。６０６は表示コントローラ（ＣＲＴＣ）であり、表示装置（ＣＲＴ）６１０上の表示制御を行なう。６０７はディスクコントローラ（ＤＫＣ）で、ブートプログラム（起動プログラム：パソコンのハードやソフトの実行（動作）を開始するプログラム）、複数のアプリケーション、編集ファイル、ユーザファイルそしてネットワーク管理プログラム等を記憶するハードディスク（ＨＤ）６１１、及びフレキシブルディスク（ＦＤ）６１２とのアクセスを制御する。
【００５１】
６０８はネットワークインタフェースカード（ＮＩＣ）で、ＬＡＮ６２０を介して、ネットワークプリンタ、他のネットワーク機器、あるいは他のＰＣと双方向のデータのやり取りを行なう。
【００５２】
また、上記コンピュータシステムにより、本発明の時系列連続データの将来予測装置がプログラム構成される。
図７に、時系列連続データの将来予測装置７００の一例を示す。この時系列連続データの将来予測装置７００は、記憶手段７０１、読み出し手段７０２、時系列データ正規化手段７０３、遅延ベクトル生成手段７０４、再構成アトラクタ生成手段７０５、近傍点探索手段７０６、将来値予測手段７０７、信頼性評価手段７０８等を有していて、外部から供給される過去の時系列信号値７０は記憶手段７０１に記憶される。
【００５３】
そして、上記記憶手段７０１に記憶された過去の時系列データ７０は読み出し手段７０２によって読み出され、時系列データ正規化手段７０３に与えられる。上記時系列データ正規化手段７０３は、過去の時系列データを正規化する手段であり、与えられた過去の時系列データ７０から教師データを生成する処理を行なう。
【００５４】
上記遅延ベクトル生成手段７０４は、遅延ベクトルを生成する手段であり、上記時系列データ正規化手段７０３によって正規化された過去の時系列データｘ（ｔ）を使って遅延ベクトルｖ（ｔ）を生成する。
【００５５】
上記再構成アトラクタ生成手段７０５は、上記遅延ベクトル生成手段７０４によって生成された遅延ベクトルを多次元の位相空間に埋め込む処理、すなわち、多次元の位相空間上で軌道を構築する処理を行なう。
【００５６】
上記近傍点探索手段７０６は、上記再構成アトラクタ生成手段７０５によって生成された再構成アトラクタ上で、現在時刻における複数の近傍点を探索する。また、上記将来値予測手段７０７は、上記近傍点探索手段７０６によって探索された複数の近傍点の挙動に基づいて将来値を予測する。上記信頼性評価手段７０８は、上記将来値予測手段７０７によって予測された将来値の信頼性を上記近傍点の軌道の平行度合いに基づいて評価する。
【００５７】
（本発明の他の実施の形態）
本発明は複数の機器から構成されるシステムに適用しても１つの機器からなる装置に適用しても良い。
【００５８】
また、前述した実施の形態の機能を実現するように各種のデバイスを動作させるように、前記各種デバイスと接続された装置あるいはシステム内のコンピュータに対し、記録媒体から、またはインターネット等の伝送媒体を介して前記実施の形態の機能を実現するためのソフトウェアのプログラムコードを供給し、そのシステムあるいは装置のコンピュータ（ＣＰＵあるいはＭＰＵ）に格納されたプログラムに従って前記各種デバイスを動作させることによって実施したものも、本発明の範疇に含まれる。
【００５９】
また、この場合、前記ソフトウェアのプログラムコード自体が前述した実施の形態の機能を実現することになり、そのプログラムコード自体、およびそのプログラムコードをコンピュータに供給するための手段、例えばかかるプログラムコードを格納した記録媒体は本発明を構成する。かかるプログラムコードを記憶する記録媒体としては、例えばフレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
【００６０】
また、コンピュータが供給されたプログラムコードを実行することにより、前述の実施の形態で説明した機能が実現されるだけでなく、そのプログラムコードがコンピュータにおいて稼働しているＯＳ（オペレーティングシステム）あるいは他のアプリケーションソフト等と共同して前述の実施の形態で示した機能が実現される場合にもかかるプログラムコードは本発明の実施の形態に含まれることは言うまでもない。
【００６１】
さらに、供給されたプログラムコードがコンピュータの機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに格納された後、そのプログラムコードの指示に基づいてその機能拡張ボードや機能拡張ユニットに備わるＣＰＵ等が実際の処理の一部または全部を行い、その処理によって前述した実施の形態の機能が実現される場合にも本発明に含まれる。
【００６２】
【発明の効果】
以上説明してきたように、本発明によれば、時系列信号値を基にして遅延ベクトルを生成し、上記生成した遅延ベクトルが描く再構成アトラクタ上で、現在時刻の近傍点の挙動から将来値を予測するとともに、上記予測した将来値の信頼性を、上記近傍点の軌道の平行度合いに基づいて評価するようにしたので、時系列連続データの将来値を高精度に予測することができる。また、将来値において大きな誤差が発生する確率が大きいと予測される場合には、予測結果に対して他の手法との併用等を行なうことが必要であるという情報をユーザーに提示することができる。
【００６３】
また、本発明のその他の特徴によれば、上記時系列データの法則性の存在程度を定常性評価指標で表現し、上記定常性評価指標に基づいて上記予測された将来値の信頼性を評価するようにしたので、予測精度の品質評価を高精度に行なうことができる。
【図面の簡単な説明】
【図１】本発明の実施の形態を示し、予測方法の概念を示す図である。
【図２】定常性評価指標の概念を説明する図である。
【図３】本発明の実施例を示し、地域冷暖房システムの日最大冷水負荷量の予測結果例を示す図である。
【図４】予測誤差と定常性評価指標との関係を示す図である。
【図５】時系列連続データの将来予測方法の処理手順の一例を示すフローチャートである。
【図６】時系列連続データの将来予測方法を実行可能なコンピュータシステムの一例を示すブロック図である。
【図７】コンピュータシステムにより構成される時系列連続データの将来予測装置の一例を示す機能構成図である。
【符号の説明】
７００  時系列連続データの将来予測装置
７０１  記憶手段
７０２  読み出し手段
７０３  時系列データ正規化手段
７０４  遅延ベクトル生成手段
７０５  再構成アトラクタ生成手段
７０６  近傍点探索手段
７０７  将来値予測手段
７０８  信頼性評価手段
【図１】
2004310500.tif 000010
【図２】
2004310500.tif 000011
【図３】
2004310500.tif 000012
【図４】
2004310500.tif 000013
【図５】
2004310500.tif 000014
【図６】
2004310500.tif 000015
【図７】
2004310500.tif 000016