(57)【特許請求の範囲】
【請求項１】
複数の時系列信号の混合信号を表すデータを取得するデータ取得部と、
前記データを用いて、前記複数の時系列信号のうちの一の時系列信号をモデリングするためのモデルパラメータであって、各ファクタの状態が初期状態である初期状態確率と、各ファクタの状態が遷移する遷移確率と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の平均値を求めるのに用いられる、各ファクタの各状態に固有の固有波形と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の分散とを含むモデルパラメータを有する確率生成モデルの各ファクタの状態にいる状態確率を求める状態推定を行う状態推定部と、
前記状態確率を用いて、前記確率生成モデルの学習を行うモデル学習部と
を備え、
前記モデル学習部は、
前記一の時系列信号に特有の制約の下で、前記固有波形を求める波形分離学習を行う波形分離学習部と、
前記分散を求める分散学習を行う分散学習部と、
前記初期状態確率と、前記遷移確率とを求める状態変動学習を行う状態変動学習部と
を有する
データ処理装置。
【請求項２】
前記確率生成モデルは、FHMM(Factorial Hidden Markov Model)である
請求項１に記載のデータ処理装置。
【請求項３】
前記モデル学習部は、各ファクタが３以上の状態を有する前記FHMMの学習を行う
請求項２に記載のデータ処理装置。
【請求項４】
前記ファクタの３以上の状態は、前記一の時系列信号に関する３以上の信号波形に対応する
請求項３に記載のデータ処理装置。
【請求項５】
前記複数の時系列信号の混合信号に関する前記FHMMを記憶するモデル記憶部をさらに備え、
前記状態推定部は、前記モデル記憶部から前記FHMMの１以上のモデルパラメータを取得して状態推定を行う
請求項２ないし４のいずれかに記載のデータ処理装置。
【請求項６】
前記状態推定部は、1時刻に状態が遷移するファクタの数を制限する状態遷移制限の下で、前記状態確率を求める
請求項１ないし５のいずれかに記載のデータ処理装置。
【請求項７】
前記状態推定部は、
前記平均値、及び、前記分散を用いて、各ファクタの状態の組み合わせにおいて、前記データが観測される観測確率を求め、
前記観測確率、及び、前記遷移確率を用いて、前記データの系列Y1，Y2，・・・，YTについて、前記データY1，Y2，・・・，Ytを観測し、時刻tに、各ファクタの状態の組み合わせzにいる前向き確率αt,zと、時刻tに、各ファクタの状態の組み合わせzにいて、その後、前記データYt，Yt+1，・・・，YTを観測する後ろ向き確率βt,zとを求め、
前記前向き確率αt,z、及び、前記後ろ向き確率βt,zを用いて、時刻tに、各ファクタの状態の組み合わせzにいる事後確率γt,zを求め、
前記事後確率γt,zを周辺化することにより、前記状態確率を求める
場合において、
各ファクタの状態の組み合わせを、パーティクルとして、
1時刻後のパーティクルを、前記状態遷移制限の下で予測し、前記前向き確率αt,zに基づいて、所定数のパーティクルをサンプリングすることを繰り返しながら、前記パーティクルとしての状態の組み合わせzの前記前向き確率αt,zを求め、
1時刻前のパーティクルを、前記状態遷移制限の下で予測し、前記後ろ向き確率βt,zに基づいて、所定数のパーティクルをサンプリングすることを繰り返しながら、前記パーティクルとしての状態の組み合わせzの前記後ろ向き確率βt,zを求める
請求項６に記載のデータ処理装置。
【請求項８】
前記分散学習部は、前記ファクタごとに個別の分散、又は、各ファクタの状態ごとに個別の分散を求める
請求項１ないし７のいずれかに記載のデータ処理装置。
【請求項９】
前記波形分離学習部は、前記固有波形を用いて求められる前記時系列信号の振幅が負の値にならない負荷制約の下で、前記固有波形を求める
請求項１ないし８のいずれかに記載のデータ処理装置。
【請求項１０】
前記波形分離学習部は、前記固有波形が、前記一の時系列信号について用意された複数の基底波形の１以上の組み合わせで表現される基底波形制約の下で、前記固有波形を求める
請求項１ないし９のいずれかに記載のデータ処理装置。
【請求項１１】
前記FHMMを使用した前記データ処理において、少なくとも一部の前記データの処理を、クラウドコンピューティングで行う
請求項２ないし５のいずれかに記載のデータ処理装置。
【請求項１２】
複数の時系列信号の混合信号を表すデータを取得するデータ取得部と、
前記データを用いて、前記複数の時系列信号のうちの一の時系列信号をモデリングするためのモデルパラメータであって、各ファクタの状態が初期状態である初期状態確率と、各ファクタの状態が遷移する遷移確率と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の平均値を求めるのに用いられる、各ファクタの各状態に固有の固有波形と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の分散とを含むモデルパラメータを有する確率生成モデルの各ファクタの状態にいる状態確率を求める状態推定を行う状態推定部と、
前記状態確率を用いて、前記確率生成モデルの学習を行うモデル学習部と
を備え、
前記モデル学習部は、
前記一の時系列信号に特有の制約の下で、前記固有波形を求める波形分離学習を行う波形分離学習部と、
前記分散を求める分散学習を行う分散学習部と、
前記初期状態確率と、前記遷移確率とを求める状態変動学習を行う状態変動学習部と
を有するデータ処理装置のデータ処理方法であって、
前記データ取得部が前記データを取得するデータ取得ステップと、
前記状態推定部が前記状態推定を行う状態推定ステップと、
前記モデル学習部が前記確率生成モデルの学習を行うモデル学習ステップと
を有し、
前記モデル学習ステップは、
前記波形分離学習部が前記波形分離学習を行う波形分離学習ステップと、
前記分散学習部が前記分散学習を行う分散学習ステップと、
前記状態変動学習部が前記状態変動学習を行う状態変動学習ステップと
を有するデータ処理方法。
【請求項１３】
複数の時系列信号の混合信号を表すデータを取得するデータ取得部と、
前記データを用いて、前記複数の時系列信号のうちの一の時系列信号をモデリングするためのモデルパラメータであって、各ファクタの状態が初期状態である初期状態確率と、各ファクタの状態が遷移する遷移確率と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の平均値を求めるのに用いられる、各ファクタの各状態に固有の固有波形と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の分散とを含むモデルパラメータを有する確率生成モデルの各ファクタの状態にいる状態確率を求める状態推定を行う状態推定部と、
前記状態確率を用いて、前記確率生成モデルの学習を行うモデル学習部であって、
前記一の時系列信号に特有の制約の下で、前記固有波形を求める波形分離学習を行う波形分離学習部と、
前記分散を求める分散学習を行う分散学習部と、
前記初期状態確率と、前記遷移確率とを求める状態変動学習を行う状態変動学習部と
を有するモデル学習部と
して、コンピュータを機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本技術は、データ処理装置、データ処理方法、及び、プログラムに関し、特に、例えば、家庭にある複数の電気機器それぞれの消費電力等を、容易かつ正確に求めることができるようにするデータ処理装置、データ処理方法、及び、プログラムに関する。
【背景技術】
【０００２】
例えば、家庭内等にある家電（家庭用電化製品）等の電気機器の個々の消費電力や消費電流を、家庭内のユーザに提示する、消費電力等のいわゆる「見える化」を実現する方法としては、例えば、各コンセントに、スマートタップを取り付ける方法がある。
【０００３】
スマートタップは、そのスマートタップが取り付けられたコンセント（に繋がっている家電）で消費されている電力を計測する計測機能と、外部との通信機能とを有する。
【０００４】
スマートタップにおいて、計測機能で計測した（消費）電力を、通信機能で、ディスプレイ等に送信し、ディスプレイにおいて、スマートタップからの消費電力を表示することにより、各家電の消費電力の「見える化」を実現することができる。
【０００５】
しかしながら、家庭内のすべてのコンセントに、スマートタップを取り付けることは、コストの面で、容易ではない。
【０００６】
また、いわゆるビルトインのエアコン（エアコンディショナ）等の、家に備え付けの家電は、コンセントを介さずに、電力線に、直接接続されている場合があり、そのような家電については、スマートタップを用いることは、困難である。
【０００７】
そこで、例えば、家庭等において、配電盤（分電盤）で計測した電流の情報から、その先に繋がっている家庭内の個々の家電の消費電力等を求めるNILM(Non-Intrusive Load Monitoring)と呼ばれる技術が注目されている。
【０００８】
NILMでは、例えば、一箇所で計測した電流等を用いて、そこから先に繋がっている個々の家電（負荷）の消費電力が、個別に計測することなく求められる。
【０００９】
例えば、特許文献１には、一箇所で計測した電流と電圧から、有効電力及び無効電力を算出し、それぞれの変化量をクラスタリングすることで、電気機器の同定を行うNILMの技術が開示されている。
【００１０】
特許文献１に記載の技術では、家電のオン時とオフ時の有効電力及び無効電力が変化することを利用するために、有効電力及び無効電力の変化点が検出される。このため、変化点の検出に失敗すると、家電を正確に同定することが困難になる。
【００１１】
また、最近の家電については、稼働状態をオンとオフの2状態で表現することが困難であり、オン時とオフ時の有効電力及び無効電力の変化を利用するだけでは、家電を正確に同定することが困難になっている。
【００１２】
そこで、例えば、特許文献２及び３には、SVM(Support Vector Machine)等のLMC（Large Margin Classifier）を、家電の識別モデル（判別モデル、Classification）として用いるNILMの技術が記載されている。
【００１３】
しかしながら、識別モデルを用いたNILMでは、HMM(Hidden Markov Model)等の生成モデルとは違い、各家電について、既知の学習データを用意し、その学習データを用いた識別モデルの学習を、事前に済ませておく必要がある。
【００１４】
このため、識別モデルを用いたNILMでは、既知の学習データを用いて識別モデルの学習を行っていない家電については、対応することが困難である。
【００１５】
そこで、例えば、非特許文献１及び２には、既知の学習データを用いた事前の学習が必要な識別モデルではない、生成モデルであるHMMを用いたNILMの技術が記載されている。
【００１６】
しかしながら、単純なHMMを用いたNILMでは、家電の数が増加すると、HMMの状態の数が膨大となり、実装することが困難となる。
【００１７】
具体的には、例えば、個々の家電の稼働状態が、オンとオフの2状態である場合に、M個の家電の稼働状態（の組み合わせ）を表現するのに必要なHMMの状態の数は、2M個になり、状態の遷移確率の数は、状態の数の2乗である(2M)2個になる。
【００１８】
したがって、家庭内の家電の数Mが、最近では多いとはいえない、例えば、20個であったとしても、HMMの状態の数は、220=1,048,576個という膨大な数になり、遷移確率の数は、その2乗の1,099,511,627,776個というテラオーダのさらに膨大な数になる。
【先行技術文献】
【特許文献】
【００１９】
【特許文献１】米国特許第4858141号明細書
【特許文献２】特開2001-330630号公報
【特許文献３】国際公開第01/077696号パンフレット
【非特許文献】
【００２０】
【非特許文献１】Bons M., Deville Y., Schang D. 1994. Non-intrusive electrical load monitoring using Hidden Markov Models.Third international Energy Efficiency and DSM Conference, October 31, Vancouver, Canada. , p. 7
【非特許文献２】中村久栄、伊藤公一、鈴木達也、「隠れマルコフモデルに基づいた電気機器の稼働状況モニタリングシステム」、電気学会論文誌Ｂ、Vol. 126、No. 12、pp. 12231229、2006年
【発明の概要】
【発明が解決しようとする課題】
【００２１】
現在においては、稼働状態が、オンとオフの2状態だけではない家電、すなわち、モードや設定等によって消費電力（電流）が変動するエアコン等の家電（変動負荷家電）のような電気機器の個々の消費電力等を、容易かつ正確に求めることができるNILMの技術の提案が要請されている。
【００２２】
本技術は、このような状況に鑑みてなされたものであり、電気機器の個々の消費電力等を、容易かつ正確に求めることができるようにするものである。
【課題を解決するための手段】
【００２３】
本技術の一側面のデータ処理装置、又は、プログラムは、複数の時系列信号の混合信号を表すデータを取得するデータ取得部と、前記データを用いて、前記複数の時系列信号のうちの一の時系列信号をモデリングするためのモデルパラメータであって、各ファクタの状態が初期状態である初期状態確率と、各ファクタの状態が遷移する遷移確率と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の平均値を求めるのに用いられる、各ファクタの各状態に固有の固有波形と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の分散とを含むモデルパラメータを有する確率生成モデルの各ファクタの状態にいる状態確率を求める状態推定を行う状態推定部と、前記状態確率を用いて、前記確率生成モデルの学習を行うモデル学習部とを備え、前記モデル学習部は、前記一の時系列信号に特有の制約の下で、前記固有波形を求める波形分離学習を行う波形分離学習部と、前記分散を求める分散学習を行う分散学習部と、前記初期状態確率と、前記遷移確率とを求める状態変動学習を行う状態変動学習部とを有するデータ処理装置、又は、データ処理装置として、コンピュータを機能させるためのプログラムである。
【００２４】
本技術の一側面のデータ処理方法は、複数の時系列信号の混合信号を表すデータを取得するデータ取得部と、前記データを用いて、前記複数の時系列信号のうちの一の時系列信号をモデリングするためのモデルパラメータであって、各ファクタの状態が初期状態である初期状態確率と、各ファクタの状態が遷移する遷移確率と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の平均値を求めるのに用いられる、各ファクタの各状態に固有の固有波形と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の分散とを含むモデルパラメータを有する確率生成モデルの各ファクタの状態にいる状態確率を求める状態推定を行う状態推定部と、前記状態確率を用いて、前記確率生成モデルの学習を行うモデル学習部とを備え、前記モデル学習部は、前記一の時系列信号に特有の制約の下で、前記固有波形を求める波形分離学習を行う波形分離学習部と、前記分散を求める分散学習を行う分散学習部と、前記初期状態確率と、前記遷移確率とを求める状態変動学習を行う状態変動学習部とを有するデータ処理装置のデータ処理方法であって、前記データ取得部が前記データを取得するデータ取得ステップと、前記状態推定部が前記状態推定を行う状態推定ステップと、前記モデル学習部が前記確率生成モデルの学習を行うモデル学習ステップとを有し、前記モデル学習ステップは、前記波形分離学習部が前記波形分離学習を行う波形分離学習ステップと、前記分散学習部が前記分散学習を行う分散学習ステップと、前記状態変動学習部が前記状態変動学習を行う状態変動学習ステップとを有するデータ処理方法である。
【００２５】
以上のような一側面においては、複数の時系列信号の混合信号を表すデータが取得され、前記データを用いて、前記複数の時系列信号のうちの一の時系列信号をモデリングするためのモデルパラメータであって、各ファクタの状態が初期状態である初期状態確率と、各ファクタの状態が遷移する遷移確率と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の平均値を求めるのに用いられる、各ファクタの各状態に固有の固有波形と、各ファクタの状態の組み合わせにおいて観測される前記データの観測値の分散とを含むモデルパラメータを有する確率生成モデルの各ファクタの状態にいる状態確率を求める状態推定が行われる。そして、前記状態確率を用いて、前記確率生成モデルの学習が行われる。前記確率生成モデルの学習では、前記一の時系列信号に特有の制約の下で、前記固有波形を求める波形分離学習と、前記分散を求める分散学習と、前記初期状態確率と、前記遷移確率とを求める状態変動学習とが行われる。
【００２６】
なお、データ処理装置は、独立した装置であっても良いし、１つの装置を構成している内部ブロックであっても良い。
【００２７】
また、プログラムは、伝送媒体を介して伝送することにより、又は、記録媒体に記録して、提供することができる。
【発明の効果】
【００２８】
本技術の一側面によれば、電気機器の個々の消費電力等を、容易かつ正確に求めることができる。
【図面の簡単な説明】
【００２９】
【図１】本技術のデータ処理装置を適用した監視システムの一実施の形態の概要を説明する図である。
【図２】家電分離で行われる波形分離学習の概要を説明する図である。
【図３】本技術を適用した監視システムの第１実施の形態の構成例を示すブロックである。
【図４】FHMMを説明する図である。
【図５】FHMMによる家電分離の定式化の概要を説明する図である。
【図６】監視システムが行う、EMアルゴリズムに従ったFHMMの学習の処理（学習処理）を説明するフローチャートである。
【図７】監視システムがステップＳ１３で行うEステップの処理を説明するフローチャートである。
【図８】FHMMの前向き確率αt,z、及び、後ろ向き確率βt,zと、HMMの前向き確率αt,i、及び、後ろ向き確率βt,jとの関係を説明する図である。
【図９】監視システムがステップＳ１４で行うMステップの処理を説明するフローチャートである。
【図１０】監視システムが行う、家電#mの情報を提示する情報提示処理を説明するフローチャートである。
【図１１】情報提示処理で行われる消費電力U(m)の表示例を示す図である。
【図１２】本技術を適用した監視システムの第２実施の形態の構成例を示すブロックである。
【図１３】監視システムがステップＳ１３で行うEステップの処理を説明するフローチャートである。
【図１４】監視システムがステップＳ１４で行うMステップの処理を説明するフローチャートである。
【図１５】本技術を適用した監視システムの第３実施の形態の構成例を示すブロックである。
【図１６】FHMMの状態の組み合わせzに、パーティクルフィルタを適用して、前向き確率αt,pを求める方法を説明する図である。
【図１７】FHMMの状態の組み合わせzに、パーティクルフィルタを適用して、後ろ向き確率βt,pを求める方法を説明する図である。
【図１８】FHMMの状態の組み合わせzに、パーティクルフィルタを適用して、事後確率γt,pを求める方法を説明する図である。
【図１９】監視システムがステップＳ１３で行うEステップの処理を説明するフローチャートである。
【図２０】監視システムがステップＳ１３で行うEステップの処理を説明するフローチャートである。
【図２１】本技術を適用した監視システムの第４実施の形態の構成例を示すブロックである。
【図２２】監視システムが負荷制約を課して行うステップＳ１４のMステップの処理を説明するフローチャートである。
【図２３】負荷制約を説明する図である。
【図２４】基底波形制約を説明する図である。
【図２５】監視システムが基底波形制約を課して行うステップＳ１４のMステップの処理を説明するフローチャートである。
【図２６】本技術を適用した監視システムの第５実施の形態の構成例を示すブロックである。
【図２７】FHMMの学習を行う監視システムによる話者分離の概要を説明する図である。
【図２８】本技術を適用した監視システムの第６実施の形態の構成例を示すブロックである。
【図２９】監視システムが行うモデル学習の処理（学習処理）を説明するフローチャートである。
【図３０】本技術を適用したコンピュータの一実施の形態の構成例を示すブロック図である。
【発明を実施するための形態】
【００３０】
［本技術の概要］
【００３１】
図１は、本技術のデータ処理装置を適用した監視システムの一実施の形態の概要を説明する図である。
【００３２】
各家庭において、電力会社から提供される電力は、配電盤（分電盤）に引き込まれ、配電盤から、家庭内の（コンセント等に接続された）家電等の電気機器に供給される。
【００３３】
本技術を適用した監視システムは、配電盤の、家庭内に電力を供給している、いわば根元等の1箇所において、家庭内の１以上の家電で消費されている電流の総和を計測し、その電流の総和の系列（電流波形系列）から、家庭内の、例えば、エアコンや掃除機等の個々の家電で消費されている電力（電流）を分離する家電分離を行う。
【００３４】
なお、監視システムに入力する入力データとしては、各家電で消費されている電流の総和そのものの他、各家電で消費されている電流の総和に関する総和データを採用することができる。
【００３５】
総和データとしては、加算が成立する値の総和を採用することができる。具体的には、総和データとしては、各家電で消費されている電流の総和そのものの他、例えば、各家電で消費されている電力の総和や、各家電で消費されている電流の波形をFFT(Fast Fourier Transform)等することによって得られる周波数成分の総和等を採用することができる。
【００３６】
また、家電分離では、総和データから、個々の家電で消費されている電力の他、個々の家電で消費されている電流に関する情報を分離することができる。具体的には、家電分離では、総和データから、例えば、個々の家電で消費されている電流や、その周波数成分を分離することができる。
【００３７】
以下の説明では、総和データとして、例えば、各家電で消費されている電流の総和を採用するとともに、家電分離において、例えば、総和データとしての電流の総和の波形から、個々の家電で消費されている電流の波形を分離することとする。
【００３８】
図２は、家電分離で行われる波形分離学習の概要を説明する図である。
【００３９】
波形分離学習では、各時刻tの総和データとしての電流波形Ytが、各家電#mで消費されている電流の電流波形W(m)の加算値（総和）であるとして、電流波形Ytから、個々の家電#mで消費されている電流波形W(m)が求められる。
【００４０】
図２では、家庭内に、5個の家電#1ないし#5があり、その5個の家電#1ないし#5のうちの、家電#1,#2,#4、及び、#5がオン状態（電力を消費する状態）になっており、家電#3がオフ状態（電力を消費しない状態）になっている。
【００４１】
そのため、図２では、総和データとしての電流波形Ytが、家電#1,#2,#4、及び、#5それぞれの消費電流W(1),W(2),W(4)、及び、W(5)の加算値（総和）になっている。
【００４２】
［本技術を適用した監視システムの第１実施の形態］
【００４３】
図３は、本技術を適用した監視システムの第１実施の形態の構成例を示すブロックである。
【００４４】
図３において、監視システムは、データ取得部１１、状態推定部１２、モデル記憶部１３、モデル学習部１４、ラベル取得部１５、及び、データ出力部１６を有する。
【００４５】
データ取得部１１は、総和データとしての電流波形Yの時系列（電流時系列）、及び、その電流波形Yに対応する電圧の波形（電圧波形）Vの時系列（電圧時系列）を取得し、状態推定部１２、モデル学習部１４、及び、データ出力部１６に供給する。
【００４６】
すなわち、データ取得部１１は、例えば、電流と電圧を計測する計測器（センサ）で構成される。
【００４７】
データ取得部１１は、例えば、配電盤等において、監視システムが設置されている家庭内の各家電で消費されている電流の総和としての電流波形Yを計測するとともに、対応する電圧波形Vを計測し、状態推定部１２、モデル学習部１４、及び、データ出力部１６に供給する。
【００４８】
状態推定部１２は、データ取得部１１からの電流波形Yと、モデル記憶部１３に記憶された、監視システムが設置されている家庭内の家電全体のモデルである全体モデル（のモデルパラメータ）φとを用いて、各家電の稼働状態を推定する状態推定を行う。そして、状態推定部１２は、状態推定の推定結果としての各家電の稼働状態Γを、モデル学習部１４、ラベル取得部１５、及び、データ出力部１６に供給する。
【００４９】
すなわち、図３において、状態推定部１２は、評価部２１、及び、推定部２２を有する。
【００５０】
評価部２１は、モデル記憶部１３に記憶された全体モデルφを構成する複数の家電モデル#1ないし#Mそれぞれの状態の各組み合わせにおいて、データ取得部１１から（状態推定部１２に）供給される電流波形Yが観測される程度を評価した評価値Eを求め、推定部２２に供給する。
【００５１】
推定部２２は、評価部２１から供給される評価値Eを用い、モデル記憶部１３に記憶された全体モデルφを構成する複数の家電モデル#1ないし#Mそれぞれにおいている状態、すなわち、家電モデル#mが表す家電（家電モデル#mによってモデル化された家電）の稼働状態Γを推定し、モデル学習部１４、ラベル取得部１５、及び、データ出力部１６に供給する。
【００５２】
モデル記憶部１３は、複数の家電全体のモデルである全体モデル（のモデルパラメータ）φを記憶する。
【００５３】
全体モデルφは、複数であるM個の、家電の（消費電流を表現する）モデルである家電モデル#1ないし#Mから構成される。
【００５４】
全体モデルのパラメータφには、家電モデル#mが表す家電の稼働状態ごとの消費電流を表す電流波形パラメータが含まれる。
【００５５】
全体モデルのパラメータφには、その他、例えば、家電モデル#mが表す家電の稼働状態の遷移（変動）を表す状態変動パラメータ、家電モデル#mが表す家電の稼働状態の初期状態を表す初期状態パラメータ、及び、全体モデルにおいて観測（生成）される電流波形Yの観測値の分散に関する分散パラメータを含めることができる。
【００５６】
モデル記憶部１３に記憶された全体モデルのモデルパラメータφは、状態推定部１２の評価部２１と推定部２２、ラベル取得部１５、及び、データ出力部１６で参照され、モデル学習部１４の後述する波形分離学習部３１、分散学習部３２、及び、状態変動学習部３３によって更新される。
【００５７】
モデル学習部１４は、データ取得部１１から供給される電流波形Yと、状態推定部１２（の推定部２２）から供給される状態推定の推定結果（各家電の稼働状態）Γとを用いて、モデル記憶部１３に記憶された全体モデルのモデルパラメータφを更新するモデル学習を行う。
【００５８】
すなわち、図３において、モデル学習部１４は、波形分離学習部３１、分散学習部３２、及び、状態変動学習部３３を有する。
【００５９】
波形分離学習部３１は、データ取得部１１から（モデル学習部１４に）供給される電流波形Yと、状態推定部１２（の推定部２２）から供給される各家電の稼働状態Γとを用いて、モデルパラメータφとしての電流波形パラメータを求める（更新する）波形分離学習を行い、その波形分離学習によって得られる電流波形パラメータによって、モデル記憶部１３に記憶された電流波形パラメータを更新する。
【００６０】
分散学習部３２は、データ取得部１１から（モデル学習部１４に）供給される電流波形Yと、状態推定部１２（の推定部２２）から供給される各家電の稼働状態Γとを用いて、モデルパラメータφとしての分散パラメータを求める（更新する）分散学習を行い、その分散学習によって得られる分散パラメータによって、モデル記憶部１３に記憶された分散パラメータを更新する。
【００６１】
状態変動学習部３３は、状態推定部１２（の推定部２２）から供給される各家電の稼働状態Γを用いて、モデルパラメータφとしての初期状態パラメータ、及び、状態変動パラメータを求める（更新する）状態変動学習を行い、その状態変動学習によって得られる初期状態パラメータ、及び、状態変動パラメータによって、モデル記憶部１３に記憶された初期状態パラメータ、及び、状態変動パラメータをそれぞれ更新する。
【００６２】
ラベル取得部１５は、状態推定部１２（の推定部２２）から供給される各家電の稼働状態Γ、モデル記憶部１３に記憶された全体モデルφ、及び、データ出力部１６で得られる、各家電モデル#mが表す家電の消費電力U(m)等を、必要に応じて用いて、各家電モデル#mが表す家電を識別するための家電ラベルL(m)を取得し、データ出力部１６に供給する。
【００６３】
データ出力部１６は、データ取得部１１から供給される電圧波形V、状態推定部１２（の推定部２２）から供給される各家電の稼働状態Γφ、及び、モデル記憶部１３に記憶された全体モデルを用いて、各家電モデル#mが表す家電の消費電力U(m)を求め、図示せぬディスプレイに表示すること等によって、ユーザに提示する。
【００６４】
なお、データ出力部１６において、各家電モデル#mが表す家電の消費電力U(m)は、ラベル取得部１５から供給される家電ラベルL(m)とともに、ユーザに提示することができる。
【００６５】
以上のように構成される監視システムにおいて、モデル記憶部１３に記憶される全体モデルとしては、例えば、FHMM(Factorial Hidden Markov Model)を採用することができる。
【００６６】
図４は、FHMMを説明する図である。
【００６７】
すなわち、図４Ａは、通常のHMMのグラフィカルモデルを示しており、図４Ｂは、FHMMのグラフィカルモデルを示している。
【００６８】
通常のHMMでは、時刻tに、その時刻tにいる１つの状態Stにおいて、１つの観測値Ytが観測される。
【００６９】
一方、FHMMでは、時刻tに、その時刻tにいる複数の状態S(1)t，S(2)t，・・・，S(M)tの組み合わせにおいて、１つの観測値Ytが観測される。
【００７０】
FHMMは、Zoubin Ghahramaniらが提案した確率生成モデルであり、その詳細については、例えば、Zoubin Ghahramani, and Michael I. Jordan, Factorial Hidden Markov Models’, Machine Learning Volume 29, Issue 2-3 ,Nov./Dec. 1997（以下、文献Ａともいう）に記載されている。
【００７１】
図５は、FHMMによる家電分離の定式化の概要を説明する図である。
【００７２】
ここで、FHMMは、複数のHMMを有する構成になっている。FHMMが有する各HMMを、ファクタともいい、m番目のファクタを、ファクタ#mとも記載する。
【００７３】
FHMMにおいて、時刻tにいる複数の状態S(1)tないしS(M)tの組み合わせとは、時刻tにいる各ファクタ#mの状態の組み合わせ（ファクタ#1の状態と、ファクタ#2の状態と、・・・、ファクタ#Mの状態とのセット）である。
【００７４】
図５は、ファクタの数Mが3のFHMMを示している。
【００７５】
家電分離においては、例えば、１つのファクタが１つの家電に対応する（１つのファクタと１つの家電とを対応させる）。図５では、ファクタ#mが、家電#mに対応している。
【００７６】
FHMMにおいて、ファクタを構成する状態の数は、ファクタごとに任意であるが、図５では、3つのファクタ#1,#2,#3の状態の数は、いずれも、4個になっている。
【００７７】
図５では、時刻t=t0において、ファクタ#1は、4つの状態#11,#12,#13,#14のうちの、状態#14（太線の丸印で示す）になっており（いて）、ファクタ#2は、4つの状態#21,#22,#23,#24のうちの、状態#21（太線の丸印で示す）になっている。また、時刻t=t0において、ファクタ#3は、4つの状態#31,#32,#33,#34のうちの、状態#33（太線の丸印で示す）になっている。
【００７８】
家電分離においては、ファクタ#mの状態は、そのファクタ#mが対応する家電#mの稼働状態に対応する。
【００７９】
例えば、家電#1に対応するファクタ#1において、状態#11は、家電#1のオフ状態に対応し、状態#14は、家電#1の、いわゆる通常モードのオン状態に対応する。また、例えば、家電#1に対応するファクタ#1において、状態#12は、家電#1の、いわゆるスリープモードのオン状態に対応し、状態#13は、家電#1の、いわゆる省エネモードのオン状態に対応する。
【００８０】
FHMMにおいて、ファクタ#mの状態#miでは、各ファクタの状態ごとに固有の波形である固有波形W(m)#miが観測（生成）される。
【００８１】
図５において、ファクタ#1では、時刻t=t0にいる状態#14において、固有波形W(1)#14が観測され、ファクタ#2では、時刻t=t0にいる状態#21において、固有波形W(2)#21が観測される。さらに、ファクタ#3では、時刻t=t0にいる状態#33において、固有波形W(3)#33が観測される。
【００８２】
FHMMでは、各ファクタにいる状態において観測される固有波形を合成した合成波形が、FHMMにおいて観測される観測値として生成される。
【００８３】
ここで、固有波形の合成としては、例えば、固有波形の総和（加算）を採用することができる。その他、固有波形の合成としては、例えば、固有波形の重み付け加算や、固有波形の論理和（固有波形の値が0と1の場合）等を採用することができるが、家電分離では、固有波形の総和が採用される。
【００８４】
FHMMの学習では、FHMMにおいて、各時刻t＝・・・,t0,t1,・・・の総和データとしての電流波形・・・，Yt0，Yt0+1，・・・が観測されるように、FHMMのモデルパラメータが求められる（更新される）。
【００８５】
モデル記憶部１３（図３）に記憶される全体モデルφとして、以上のようなFHMMを採用する場合、全体モデルφを構成する家電モデル#mは、ファクタ#mに相当する。
【００８６】
なお、FHMMのファクタの数Mとしては、家庭内に存在すると想定される家電の最大数よりもマージンとしての所定数だけ大きい値が採用される。
【００８７】
また、全体モデルφとしてのFHMMとしては、各ファクタが3以上の状態を有するFHMMが採用される。
【００８８】
ファクタが2状態しか有しない場合には、そのファクタに対応する家電の稼働状態として、例えば、オフ状態とオン状態との２状態しか表現することができず、モードや設定等によって消費電力（電流）が変動するエアコン等の家電（以下、変動負荷家電ともいう）について、正確な消費電力等を求めることが困難になるからである。
【００８９】
すなわち、全体モデルφとしてのFHMMとして、各ファクタが3以上の状態を有するFHMMを採用することにより、変動負荷家電について、正確な消費電力等を求めることが可能となる。
【００９０】
全体モデルφとして、FHMMを採用する場合、FHMMにおいて観測される電流波形Ytの系列と、各ファクタ#mの状態S(m)tの組み合わせStの系列との同時分布P({St,Yt})は、マルコフ性を仮定することにより、式（１）で計算される。
【００９１】
【数１】
0006020880.tif 000002
・・・（１）
【００９２】
ここで、同時分布P({St,Yt})は、時刻tにおいて、各ファクタ#mの状態S(m)tの組み合わせ（M個のファクタそれぞれの状態の組み合わせ）Stにおいて、電流波形Ytが観測される確率を表す。
【００９３】
P(S１)は、最初の時刻t=1において、各ファクタ#mの状態S(m)1の組み合わせS1にいる初期状態確率を表す。
【００９４】
P(St｜St-1)は、時刻t-1に、状態の組み合わせSt-1にいて、時刻tに、状態の組み合わせStに遷移する遷移確率を表す。
【００９５】
P(Yt｜St)は、時刻tに、状態の組み合わせStにおいて、電流波形Ytが観測される観測確率を表す。
【００９６】
時刻tにおける状態の組み合わせStは、M個のファクタ#1ないし#Mそれぞれの、時刻tにおける状態S(1)t，S(2)t，・・・，S(M)tの組み合わせであり、式St＝{S(1)t，S(2)t，・・・，S(M)t}で表現される。
【００９７】
なお、家電#mの稼働状態は、他の家電#m'とは独立に変動すると仮定し、ファクタ#mの状態S(m)tは、他のファクタ#m'の状態S(m')tとは独立に遷移することとする。
【００９８】
また、FHMMのファクタ#mとしてのHMMの状態の数K(m)としては、他のファクタ#m'としてのHMMの状態の数K(m')とは独立な数を採用することができる。但し、ここでは、説明を簡単にするため、ファクタ#1ないし#Mの状態の数K(1)ないしK(M)は、式K(1)＝K(2)＝・・・＝K(M)＝Kで表される通り、同一の数Kであることとする。
【００９９】
FHMMにおいて、式（１）の同時分布P({St,Yt})の計算に必要な初期状態確率P(S１)、遷移確率P(St｜St-1)、及び、観測確率P(Yt｜St)は、以下のように計算することができる。
【０１００】
すなわち、初期状態確率P(S１)は、式（２）に従って計算することができる。
【０１０１】
【数２】
0006020880.tif 000003
・・・（２）
【０１０２】
ここで、P(S(m)1)は、ファクタ#mの状態S(m)1が、最初の時刻t=1にいる状態（初期状態）である初期状態確率を表す。
【０１０３】
初期状態確率P(S(m)1)は、例えば、ファクタ#mのk番目（k＝1,2,・・・,K）の状態の初期状態確率を、k行目のコンポーネントとするK行の列ベクトル（K行1列の行列）である。
【０１０４】
遷移確率P(St｜St-1)は、式（３）に従って計算することができる。
【０１０５】
【数３】
0006020880.tif 000004
・・・（３）
【０１０６】
ここで、P(S(m)t|S(m)t-1)は、ファクタ#mにおいて、時刻t-1に状態S(m)t-1にいて、時刻tに状態S(m)tに遷移する遷移確率を表す。
【０１０７】
遷移確率P(S(m)t|S(m)t-1)は、例えば、ファクタ#mのk番目の状態#kから、k'番目（k'＝1,2,・・・,K）の状態#k'に遷移する遷移確率を、k行k'列目のコンポーネントとするK行K列の行列（正方行列）である。
【０１０８】
観測確率P(Yt｜St)は、式（４）に従って計算することができる。
【０１０９】
【数４】
0006020880.tif 000005
・・・（４）
【０１１０】
ここで、ダッシュ（’）は転置を表し、上付きの−１は、逆数（逆行列）を表す。また、｜C｜は、Cの絶対値（行列式）（determinant演算）を表す。
【０１１１】
また、Dは、観測値Ytの次元を表す。
【０１１２】
例えば、図３のデータ取得部１１は、電圧が、負の値から正の値に変化するときのゼロ交差のあるタイミングを、電流の位相が0になっているタイミングとして、所定のサンプリング間隔で、１周期分（日本では、1/50又は1/60秒）の電流をサンプリングし、そのサンプリング値をコンポーネントとする列ベクトルを、１時刻分の電流波形Ytとして出力する。
【０１１３】
データ取得部１１が１周期分の電流をサンプリングするサンプリングの回数がD回であるとすると、電流波形Ytは、D行の列ベクトルである。
【０１１４】
式（４）の観測確率P(Yt｜St)によれば、観測値Ytは、平均値（平均ベクトル）がμtで、分散（共分散行列）がCの正規分布に従う。
【０１１５】
平均値μtは、例えば、電流波形Ytと同様のD行の列ベクトルであり、分散Cは、例えば、D行D列の行列（対角成分が分散の行列）である。
【０１１６】
平均値μtは、図５で説明した固有波形W(m)を用い、式（５）で表される。
【０１１７】
【数５】
0006020880.tif 000006
・・・（５）
【０１１８】
ここで、ファクタ#mの状態#kの固有波形を、W(m)kと表すこととすると、ファクタ#mの状態#kの固有波形W(m)kは、例えば、電流波形Ytと同様のD行の列ベクトルである。
【０１１９】
また、固有波形W(m)は、ファクタ#mの各状態#1,#2,・・・,#Kの固有波形W(m)1,W(m)2,・・・,W(m)Kの集まりであり、ファクタ#mの状態#kの固有波形W(m)kである列ベクトルを、k列目のコンポーネントとするD行K列の行列である。
【０１２０】
さらに、S*(m)tは、時刻tにいるファクタ#mの状態を表し、以下、S*(m)tを、時刻tのファクタ#mの現在状態ともいう。時刻tのファクタ#mの現在状態S*(m)tは、例えば、式（６）に示すように、K行のうちの1行のみのコンポーネントが0で、他の行のコンポーネントが0のK行の列ベクトルである。
【０１２１】
【数６】
0006020880.tif 000007
・・・（６）
【０１２２】
時刻tにいるファクタ#mの状態が、状態#kである場合、時刻tのファクタ#mの現在状態S*(m)tとしてのK行の列ベクトルS*(m)tは、k行目のコンポーネントだけが1とされ、他のコンポーネントは0とされる。
【０１２３】
式（５）によれば、時刻tにいる各ファクタ#mの状態#kの固有波形W(m)kの総和が、時刻tの電流波形Ytの平均値μtとして求められる。
【０１２４】
FHMMのモデルパラメータφは、式（２）の初期状態確率P(S(m)1)、式（３）の遷移確率P(S(m)t|S(m)t-1)、式（４）の分散C、及び、式（５）の固有波形W(m)（＝W(m)1,W(m)2,・・・,W(m)K）であり、図３のモデル学習部１４では、これらのFHMMのモデルパラメータφが求められる。
【０１２５】
すなわち、波形分離学習部３１は、波形分離学習によって、固有波形W(m)を、電流波形パラメータとして求める。分散学習部３２では、分散学習によって、分散Cを、分散パラメータとして求める。状態変動学習部３３は、状態変動学習によって、初期状態確率P(S(m)1)、及び、遷移確率P(S(m)t|S(m)t-1)を、それぞれ、初期状態パラメータ、及び、状態変動パラメータとして求める。
【０１２６】
ここで、前述したように、例えば、個々の家電の稼働状態が、オンとオフの2状態であっても、通常のHMMによって、20個の家電の稼働状態（の組み合わせ）を表現する場合には、HMMの状態の数は、220=1,048,576個になり、遷移確率の数は、その2乗の1,099,511,627,776個になる。
【０１２７】
一方、FHMMによれば、稼働状態がオンとオフの2状態しかないM個の家電は、各ファクタが2状態の、M個のファクタによって表現することができる。したがって、各ファクタにおいて、状態の数は、2個であり、遷移確率の数は、その2乗の4個であるから、FHMMによって、M=20個の家電（ファクタ）の稼働状態を表現する場合には、FHMMの状態の数（総数）は、40＝2×20個という少ない数で済み、遷移確率の数も、80＝4×20個という少ない数で済む。
【０１２８】
FHMMの学習、すなわち、FHMMのモデルパラメータφとしての初期状態確率P(S(m)1)、遷移確率P(S(m)t|S(m)t-1)、分散C、及び、固有波形W(m)の更新は、例えば、文献Ａに記載されているように、EM(Expectation-Maximization)アルゴリズムに従って行うことができる。
【０１２９】
EMアルゴリズムによるFHMMの学習では、式（７）の条件付き完全データ対数尤度の期待値Q(φnew|φ)を最大化するために、Eステップの処理とMステップの処理とが、交互に繰り返される。
【０１３０】
【数７】
0006020880.tif 000008
・・・（７）
【０１３１】
ここで、条件付き完全データ対数尤度の期待値Q(φnew|φ)とは、モデルパラメータφの下で、完全データ{St,Yt}が観測される場合において、新たなモデルパラメータφnewの下で、完全データ{St,Yt}が観測される対数尤度log(P({St,Yt}|φnew))の期待値を意味する。
【０１３２】
EMアルゴリズムのEステップの処理では、式（７）の条件付き完全データ対数尤度の期待値Q(φnew|φ)（に相当する値）が求められ、EMアルゴリズムのMステップの処理では、Eステップの処理で求められた期待値Q(φnew|φ)を、より大にする新たなモデルパラメータφnewが求められ、モデルパラメータφが、（期待値Q(φnew|φ)をより大にする）新たなモデルパラメータφnewに更新される。
【０１３３】
図６は、監視システム（図３）が行う、EMアルゴリズムに従ったFHMMの学習の処理（学習処理）を説明するフローチャートである。
【０１３４】
ステップＳ１１において、モデル学習部１４は、モデル記憶部１３に記憶されたFHMMのモデルパラメータφとしての初期状態確率P(S(m)1)、遷移確率P(S(m)t|S(m)t-1)、分散C、及び、固有波形W(m)を初期化し、処理は、ステップＳ１２に進む。
【０１３５】
ここで、初期状態確率P(S(m)1)としてのK行の列ベクトルのk行目のコンポーネント、すなわち、ファクタ#mのk番目の初期状態確率π(m)kは、例えば、1/Kに初期化される。
【０１３６】
遷移確率P(S(m)t|S(m)t-1)としてのK行K列の行列のi行j列目のコンポーネント（i,j＝1,2,・・・,K）、すなわち、ファクタ#mにおいて、i番目の状態#iからj番目の状態#jに遷移する遷移確率P(m)i,jは、例えば、乱数を用いて、式P(m)i,1＋P(m)i,2＋・・・＋P(m)i,K＝1を満たすように、初期化される。
【０１３７】
分散CとしてのD行D列の行列は、例えば、乱数を用いて、対角成分を設定し、他の成分（コンポーネント）を0としたD行D列の対角行列に初期化される。
【０１３８】
固有波形W(m)としてのD行K列の行列のk列目の列ベクトル、すなわち、ファクタ#mの状態#kの固有波形W(m)kとしてのD行の列ベクトルの各コンポーネントは、例えば、乱数を用いて初期化される。
【０１３９】
ステップＳ１２では、データ取得部１１が、所定の時間T分の電流波形を取得し、各時刻t=1,2・・・,Tの電流波形を、計測波形Y1,Y2,・・・,YTとして、状態推定部１２、及び、モデル学習部１４に供給して、処理は、ステップＳ１３に進む。
【０１４０】
ここで、データ取得部１１は、時刻t=1,2・・・,Tの電流波形とともに、電圧波形も取得する。データ取得部１１は、時刻t=1,2・・・,Tの電圧波形を、データ出力部１６に供給する。
【０１４１】
データ出力部１６では、データ取得部１１からの電圧波形が、消費電力の算出に用いられる。
【０１４２】
ステップＳ１３では、状態推定部１２が、データ取得部１１からの計測波形Y1ないしYTを用いて、Eステップの処理を行い、処理は、ステップＳ１４に進む。
【０１４３】
すなわち、ステップＳ１３では、状態推定部１２は、データ取得部１１からの計測波形Y1ないしYTを用いて、モデル記憶部１３に記憶されたFHMMの各ファクタ#mの各状態にいる状態確率等を求める状態推定を行い、その状態推定の推定結果を、モデル学習部１４、及び、データ出力部１６に供給する。
【０１４４】
ここで、図５で説明したように、家電分離においては、ファクタ#mの状態は、そのファクタ#mが対応する家電#mの稼働状態に対応する。したがって、FHMMのファクタ#mの状態#kにいる状態確率は、家電#mの稼働状態が状態#kである程度を表すので、そのような状態確率を求める状態推定は、家電の稼働状態を求めている（推定している）、ということができる。
【０１４５】
ステップＳ１４では、モデル学習部１４は、データ取得部１１からの計測波形Y1ないしYT、及び、状態推定部１２からの状態推定の推定結果を用いて、Mステップの処理を行い、処理は、ステップＳ１５に進む。
【０１４６】
すなわち、ステップＳ１４では、モデル学習部１４は、データ取得部１１からの計測波形Y1ないしYT、及び、状態推定部１２からの状態推定の推定結果を用いて、モデル記憶部１３に記憶された、各ファクタが３以上の状態を有するFHMMの学習を行うことにより、モデル記憶部１３に記憶されたFHMMのモデルパラメータφとしての初期状態確率π(m)k、遷移確率P(m)i,j、分散C、及び、固有波形W(m)を更新する。
【０１４７】
ステップＳ１５では、モデル学習部１４は、モデルパラメータφの収束条件が満たされているかどうかを判定する。
【０１４８】
ここで、モデルパラメータφの収束条件としては、例えば、Eステップ及びMステップの処理があらかじめ設定された所定の回数だけ繰り返されたことや、FHMMにおいて、計測波形Y1ないしYTが観測される尤度の、モデルパラメータφの更新前と更新後との間の変化量が、あらかじめ設定された閾値以内であること等を採用することができる。
【０１４９】
ステップＳ１５で、モデルパラメータφの収束条件が満たされていないと判定された場合、処理は、ステップＳ１３に戻り、以下、同様の処理が繰り返される。
【０１５０】
また、ステップＳ１５において、モデルパラメータφの収束条件が満たされると判定された場合、学習処理は、終了する。
【０１５１】
図７は、図３の監視システムが、図６のステップＳ１３で行うEステップの処理を説明するフローチャートである。
【０１５２】
ステップＳ２１において、評価部２１は、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの分散C、及び、固有波形W(m)、並びに、データ取得部１１からの計測波形Yt＝{Y1,Y2,・・・,YT}を用い、各時刻t＝{1,2,・・・,T}の、状態の各組み合わせStについて、式（４）の観測確率P(Yt|St)を、評価値Eとして求め、推定部２２に供給して、処理は、ステップＳ２２に進む。
【０１５３】
ステップＳ２２では、推定部２２が、評価部２１からの観測確率P(Yt|St)、及び、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの遷移確率P(m)i,j（及び初期状態確率π(m)）を用いて、計測波形Y1,Y2,・・・,Ytを観測し、時刻tに、状態の組み合わせ(時刻tにいるファクタ#1の状態と、ファクタ#2の状態と、・・・、ファクタ#Mの状態との組み合わせ)zにいる前向き確率αt,zを求め、処理は、ステップＳ２３に進む。
【０１５４】
ここで、HMMの前向き確率の求め方については、例えば、C.M. ビショップ, 「パターン認識と機械学習 下 ベイズ理論による統計的予測」, シュプリンガー・ジャパン, 2008（以下、文献Ｂともいう）の336ページに記載されている。
【０１５５】
前向き確率αt,zは、例えば、1時刻前の前向き確率αt-1,wを用いた漸化式αt,z＝Σαt-1,wP(z|w)P(Yt|z)に従って求めることができる。
【０１５６】
漸化式αt,z＝Σαt-1,wP(z|w)P(Yt|z)において、Σは、wを、FHMMの状態の組み合わせのすべてに変えてとるサメーションを表す。
【０１５７】
また、漸化式αt,z＝Σαt-1,wP(z|w)P(Yt|z)において、wは、1時刻前の時刻t-1にいる状態の組み合わせを表す。P(z|w)は、時刻t-1に、状態の組み合わせwにいて、時刻tに状態の組み合わせzに遷移する遷移確率を表し、P(Yt|z)は、時刻tに、状態の組み合わせzにおいて、計測波形Ytを観測する観測確率を表す。
【０１５８】
なお、前向き確率αt,zの初期値、すなわち、時刻t=1のときの前向き確率α1,zとしては、状態の組み合わせzを構成する各ファクタ#mの状態#kの初期状態確率π(m)kの積が採用される。
【０１５９】
ステップＳ２３では、推定部２２が、評価部２１からの観測確率P(Yt|St)、及び、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの遷移確率P(m)i,jを用いて、時刻tに、状態の組み合わせzにいて、その後、計測波形Yt,Yt+1,・・・,YTを観測する後ろ向き確率βt,zを求め、処理は、ステップＳ２４に進む。
【０１６０】
ここで、HMMの後ろ向き確率の求め方については、例えば、上述の文献Ｂの336ページに記載されている。
【０１６１】
後ろ向き確率βt,zは、例えば、1時刻後の後ろ向き確率βt+1,wを用いた漸化式βt,z＝ΣP(Yt|z)P(w|z)βt+1,wに従って求めることができる。
【０１６２】
漸化式βt,z＝ΣP(Yt|z)P(w|z)βt+1,wにおいて、Σは、wを、FHMMの状態の組み合わせのすべてに変えてとるサメーションを表す。
【０１６３】
また、漸化式βt,z＝ΣP(Yt|z)P(w|z)βt+1,wにおいて、wは、1時刻後の時刻t+1にいる状態の組み合わせを表す。P(w|z)は、時刻tに、状態の組み合わせzにいて、時刻t+1に状態の組み合わせwに遷移する遷移確率を表し、P(Yt|z)は、時刻tに、状態の組み合わせzにおいて、計測波形Ytを観測する観測確率を表す。
【０１６４】
なお、後ろ向き確率βt,zの初期値、すなわち、時刻t=Tのときの後ろ向き確率βT,zとしては、1が採用される。
【０１６５】
ステップＳ２４では、推定部２２は、前向き確率αt,zと後ろ向き確率βt,zを用いて、式（８）に従い、全体モデルφとしてのFHMMにおいて、時刻tに、状態の組み合わせzにいる事後確率γt,zを求め、処理は、ステップＳ２５に進む。
【０１６６】
【数８】
0006020880.tif 000009
・・・（８）
【０１６７】
ここで、式（８）の右辺の分母のΣは、wを、時刻tに取り得る状態の組み合わせStのすべてに変えてとるサメーションを表す。
【０１６８】
式（８）によれば、事後確率γt,zは、前向き確率αt,zと後ろ向き確率βt,zとの積αt,zβt,zを、その積αt,wβt,wの、FHMMが取り得る状態の組み合わせw∈Stについての総和Σαt,wβt,wで正規化することにより求められる。
【０１６９】
ステップＳ２５では、推定部２２は、事後確率γt,zを用い、ファクタ#mにおいて、時刻tに、状態S(m)tにいる事後確率<S(m)t>と、時刻tに、ファクタ#mにおいて、状態S(m)tにいて、他のファクタ#nにおいて、状態S(n)tにいる事後確率<S(m)tS(n)t'>とを求め、処理は、ステップＳ２６に進む。
【０１７０】
ここで、事後確率<S(m)t>は、式（９）に従って求められる。
【０１７１】
【数９】
0006020880.tif 000010
・・・（９）
【０１７２】
式（９）によれば、ファクタ#mにおいて、時刻tに、状態S(m)tにいる事後確率<S(m)t>は、時刻tに、状態の組み合わせzにいる事後確率γt,zを、ファクタ#mの状態を含まない状態の組み合わせzについて周辺化することにより求められる。
【０１７３】
なお、事後確率<S(m)t>は、例えば、時刻tに、ファクタ#mのK個の状態のうちのk番目の状態いる状態確率（事後確率）を、k行目のコンポーネントとするK行の列ベクトルであることとする。
【０１７４】
事後確率<S(m)tS(n)t'>は、式（１０）に従って求められる。
【０１７５】
【数１０】
0006020880.tif 000011
・・・（１０）
【０１７６】
式（１０）によれば、時刻tに、ファクタ#mにおいて、状態S(m)tにいて、他のファクタ#nにおいて、状態S(n)tにいる事後確率<S(m)tS(n)t'>は、時刻tに、状態の組み合わせzにいる事後確率γt,zを、ファクタ#mの状態と、ファクタ#nの状態との両方を含まない状態の組み合わせzについて周辺化することにより求められる。
【０１７７】
なお、事後確率<S(m)tS(n)t'>は、例えば、時刻tに、ファクタ#mの状態#kと、他のファクタ#nの状態#k'とにいる状態確率（事後確率）を、k行k'列目のコンポーネントとするK行K列の行列であることとする。
【０１７８】
ステップＳ２６では、推定部２２は、前向き確率αt,z、後ろ向き確率βt,z、遷移確率P(z|w)、及び、評価部２１からの観測確率P(Yt|St)を用いて、ファクタ#mにおいて、時刻t-1に、状態S(m)t-1にいて、次の時刻tに、状態S(m)tにいる事後確率<S(m)t-1S(m)t'>を求める。
【０１７９】
そして、推定部２２は、事後確率<S(m)t>，<S(m)tS(n)t'>、及び、<S(m)t-1S(m)t'>を、状態推定の推定結果として、モデル学習部１４、ラベル取得部１５、及び、データ出力部１６に供給し、Eステップの処理からリターンする。
【０１８０】
ここで、事後確率<S(m)t-1S(m)t'>は、式（１１）に従って求められる。
【０１８１】
【数１１】
0006020880.tif 000012
・・・（１１）
【０１８２】
式（１１）の事後確率<S(m)t-1S(m)t'>の算出にあたり、状態の組み合わせwから状態の組み合わせzに遷移する遷移確率P(z|w)は、式（３）に従い、状態の組み合わせwを構成するファクタ#1の状態#i(1)から状態の組み合わせzを構成するファクタ#1の状態#j(1)への遷移確率P(1)i(1),j(1)、状態の組み合わせwを構成するファクタ#2の状態#i(2)から状態の組み合わせzを構成するファクタ#2の状態#j(2)への遷移確率P(2)i(2),j(2)、・・・、及び、状態の組み合わせwを構成するファクタ#Mの状態#i(M)から状態の組み合わせzを構成するファクタ#Mの状態#j(M)への遷移確率P(M)i(M),j(M)の積P(1)i(1),j(1)×P(2)i(2),j(2)×・・・×P(M)i(M),j(M)として求められる。
【０１８３】
なお、事後確率<S(m)t-1S(m)t'>は、例えば、ファクタ#mにおいて、時刻t-1に状態#iにいて、次の時刻tに状態jにいる状態確率（事後確率）を、i行j列目のコンポーネントとするK行K列の行列であることとする。
【０１８４】
図８は、FHMMの前向き確率αt,z、及び、後ろ向き確率βt,zと、（通常の）HMMの前向き確率αt,i、及び、後ろ向き確率βt,jとの関係を説明する図である。
【０１８５】
FHMMについては、そのFHMMと等価なHMMを構成することができる。
【０１８６】
あるFHMMと等価なHMMは、そのFHMMの各ファクタの状態の組み合わせzに相当する状態を有する。
【０１８７】
そして、FHMMの前向き確率αt,z、及び、後ろ向き確率βt,zは、そのFHMMと等価なHMMの前向き確率αt,i、及び、後ろ向き確率βt,jに一致する。
【０１８８】
図８Ａは、2つの状態#1及び#2を有するファクタ#1と#2からなるFHMMを示している。
【０１８９】
図８ＡのFHMMでは、ファクタ#1の状態#kとファクタ#2の状態#k'との組み合わせz＝[k,k']として、ファクタ#1の状態#1とファクタ#2の状態#1との組み合わせ[1,1]、ファクタ#1の状態#1とファクタ#2の状態#2との組み合わせ[1,2]、ファクタ#1の状態#2とファクタ#2の状態#1との組み合わせ[2,1]、及び、ファクタ#1の状態#2とファクタ#2の状態#2との組み合わせ[2,2]の4通りがある。
【０１９０】
図８Ｂは、図８ＡのFHMMと等価なHMMを示している。
【０１９１】
図８ＢのHMMは、図８ＡのFHMMの状態の4つの組み合わせ[1,1]，[1,2]，[2,1]、及び、[2,2]それぞれに相当する4つの状態#(1,1)，#(1,2)，#(2,1)、及び、#(2,2)を有する。
【０１９２】
そして、図８ＡのFHMMの前向き確率αt,z＝{αt,[1,1]，αt,[1,2]，αt,[2,1]，αt,[2,2]}は、図８ＢのHMMの前向き確率αt,i＝{αt,(1,1)，αt,(1,2)，αt,(2,1)，αt,(2,2)}と一致する。
【０１９３】
同様に、図８ＡのFHMMの後ろ向き確率βt,z＝{βt,[1,1]，βt,[1,2]，βt,[2,1]，βt,[2,2]}は、図８ＢのHMMの前向き確率βt,i＝{βt,(1,1)，βt,(1,2)，βt,(2,1)，βt,(2,2)}と一致する。
【０１９４】
例えば、上述の式（８）の右辺の分母、つまり、積αt,wβt,wの、FHMMが取り得る状態の組み合わせw∈Stについての総和Σαt,wβt,wは、図８ＡのFHMMについては、式Σαt,wβt,w＝αt,[1,1]βt,[1,1]＋αt,[1,2]βt,[1,2]＋αt,[2,1]βt,[2,1]＋αt,[2,2]βt,[2,2]で表される。
【０１９５】
図９は、図３の監視システムが、図６のステップＳ１４で行うMステップの処理を説明するフローチャートである。
【０１９６】
ステップＳ３１において、波形分離学習部３１は、データ取得部１１からの計測波形Yt、並びに、推定部２２からの事後確率<S(m)t>、及び、<S(m)tS(n)t'>を用いて、波形分離学習を行うことにより、固有波形W(m)の更新値W(m)newを求め、その更新値W(m)newによって、モデル記憶部１３に記憶された固有波形W(m)を更新して、処理は、ステップＳ３２に進む。
【０１９７】
すなわち、波形分離学習部３１は、波形分離学習として、式（１２）を計算することにより、固有波形W(m)の更新値W(m)newを求める。
【０１９８】
【数１２】
0006020880.tif 000013
・・・（１２）
【０１９９】
ここで、Wnewは、D行K列の行列であるファクタ#mの固有波形W(m)の更新値W(m)newが、ファクタ（のインデクス）#m順に、左から右方向に並んだD行K×M列の行列である。D行K×M列の行列である固有波形（の更新値）Wnewの(m-1)K+k列の列ベクトルは、ファクタ#mの状態#kの固有波形W(m)k（の更新値）になっている。
【０２００】
<St'>は、K行の列ベクトルである事後確率<S(m)t>が、ファクタ#m順に、上から下方向に並んだK×M行の列ベクトルを転置したK×M列の行ベクトルである。K×M列の行ベクトルである事後確率<St'>の(m-1)K+k列目のコンポーネントは、時刻tに、ファクタ#mの状態#kにいる状態確率になっている。
【０２０１】
<StSt'>は、K行K列の行列である事後確率<S(m)tS(n)t'>が、ファクタ#m順に、上から下方向に並び、かつ、ファクタ#n順に、左から右方向に並んだK×M行K×M列の行列である。K×M行K×M列の行列である事後確率<StSt'>の(m-1)K+k行(n-1)K+k'列目のコンポーネントは、時刻tに、ファクタ#mの状態#kと、他のファクタ#nの状態#k'とにいる状態確率になっている。
【０２０２】
上付のアスタリスク(*)は、逆行列、又は、擬似逆行列を表す。
【０２０３】
式（１２）を計算する波形分離学習によれば、計測波形Ytと式（５）の平均値μt＝ΣW(m)S*(m)tとの誤差が、なるべく小さくなるように、計測波形Ytが、固有波形W(m)に分離される。
【０２０４】
ステップＳ３２では、分散学習部３２が、データ取得部１１からの計測波形Yt、推定部２２からの事後確率<S(m)t>、及び、及び、モデル記憶部１３に記憶された固有波形W(m)を用いて、分散学習を行うことにより、分散Cの更新値Cnewを求め、モデル記憶部１３に記憶された分散Cを更新して、処理は、ステップＳ３３に進む。
【０２０５】
すなわち、分散学習部３２は、分散学習として、式（１３）を計算することにより、分散Cの更新値Cnewを求める。
【０２０６】
【数１３】
0006020880.tif 000014
・・・（１３）
【０２０７】
ステップＳ３３では、状態変動学習部３３が、推定部２２からの事後確率<S(m)t>、及び、<S(m)t-1S(m)t'>を用いて、状態変動学習を行うことにより、遷移確率P(m)i,jの更新値P(m)i,jnew、及び、初期状態確率π(m)の更新値π(m)newを求め、その更新値P(m)i,jnew、及び、π(m)newによって、モデル記憶部１３に記憶された遷移確率P(m)i,j、及び、初期状態確率π(m)を更新して、Mステップの処理からリターンする。
【０２０８】
すなわち、状態変動学習部３３は、状態変動学習として、式（１４）、及び、式（１５）を計算することにより、それぞれ、遷移確率P(m)i,jの更新値P(m)i,jnew、及び、初期状態確率π(m)の更新値π(m)newを求める。
【０２０９】
【数１４】
0006020880.tif 000015
・・・（１４）
【０２１０】
【数１５】
0006020880.tif 000016
・・・（１５）
【０２１１】
ここで、<S(m)t-1,iS(m)t,j>は、K行K列の行列である事後確率<S(m)t-1S(m)t'>のi行j列目のコンポーネントであり、ファクタ#mにおいて、時刻t-1に状態#iにいて、次の時刻tに状態#jにいる状態確率を表す。
【０２１２】
<S(m)t-1,i>は、K行の列ベクトル事後確率<S(m)t-1>のi行目のコンポーネントであり、時刻t-1に、ファクタ#mの状態#iいる状態確率を表す。
【０２１３】
π(m)（π(m)new）は、ファクタ#mの状態#kの初期状態確率π(m)k（の更新値π(m)knew）を、k行目のコンポーネントとするK行の列ベクトルである。
【０２１４】
図１０は、監視システム（図３）が行う、家電#mの情報を提示する情報提示処理を説明するフローチャートである。
【０２１５】
ステップＳ４１において、データ出力部１６は、データ取得部１１からの電圧波形（電流波形である計測波形Ytに対応する電圧波形）Vt、状態推定部１２からの状態推定の推定結果としての事後確率<S(m)t>、及び、モデル記憶部１３に記憶された固有波形W(m)を用いて、各ファクタ#mの消費電力U(m)を求め、処理は、ステップＳ４２に進む。
【０２１６】
ここで、データ出力部１６は、時刻tの電圧波形Vtと、時刻tのファクタ#mに対応する家電#mの消費電流Atとを用いて、時刻tのファクタ#mに対応する家電#mの消費電力U(m)を求める。
【０２１７】
データ出力部１６において、時刻tのファクタ#mに対応する家電#mの消費電流Atは、以下のようにして求められる。
【０２１８】
すなわち、データ出力部１６は、例えば、ファクタ#mにおいて、事後確率<S(m)t>が最大の状態#kの固有波形W(m)を、時刻tのファクタ#mに対応する家電#mの消費電流Atとして求める。
【０２１９】
また、データ出力部１６は、例えば、K行の列ベクトルである事後確率<S(m)t>のコンポーネントである、時刻tのファクタ#mの各状態の状態確率を重みとして用いた、ファクタ#mの各状態の固有波形W(m)1，W(m)2，・・・，W(m)Kの重み付け加算値を、時刻tのファクタ#mに対応する家電#mの消費電流Atとして求める。
【０２２０】
なお、FHMMの学習が進行し、ファクタ#mが、家電#mを適切に表現する家電モデルとなると、時刻tのファクタ#mの各状態の状態確率は、時刻tの家電#mの稼働状態に対応する状態の状態確率がほぼ1になり、残りのK-1個の状態の状態確率はほぼ0となる。
【０２２１】
その結果、ファクタ#mにおいて、事後確率<S(m)t>が最大の状態#kの固有波形W(m)と、時刻tのファクタ#mの各状態の状態確率を重みとして用いた、ファクタ#mの各状態の固有波形W(m)1，W(m)2，・・・，W(m)Kの重み付け加算値とは、ほぼ同一となる。
【０２２２】
ステップＳ４２では、ラベル取得部１５が、各家電モデル#mが表す家電#m、すなわち、FHMMの各ファクタ#mに対応する家電#mを識別する家電ラベルL(m)を取得し、データ出力部１６に供給して、処理は、ステップＳ４３に進む。
【０２２３】
ここで、ラベル取得部１５では、例えば、データ出力部１６において求められる各ファクタ#mに対応する家電#mの消費電流Atや、消費電力U(m)、その消費電力U(m)から認識される家電#mの使用時間帯を表示し、その消費電流Atや、消費電力U(m)、使用時間帯にあてはまる家電の名称を、ユーザに入力してもらうことで、ユーザにより入力された家電の名称を、家電ラベルL(m)として取得することができる。
【０２２４】
また、ラベル取得部１５では、例えば、様々な家電について、その消費電力や、電流波形（消費電流）、使用時間帯等の属性と、家電の名称とを対応付けて登録したデータベースを事前に用意しておき、データベースにおいて、データ出力部１６で求められる各ファクタ#mに対応する家電#mの消費電流Atや、消費電力U(m)、その消費電力U(m)から認識される家電#mの使用時間帯に対応付けられている家電の名称を、家電ラベルL(m)として取得することができる。
【０２２５】
なお、ラベル取得部１５において、既に、家電ラベルL(m)が取得され、データ出力部１６に供給されている家電#mに対応するファクタ#mについては、ステップＳ４２の処理をスキップすることができる。
【０２２６】
ステップＳ４３では、データ出力部１６が、各ファクタ#mの消費電力U(m)を、そのファクタ#mの家電ラベルL(m)とともに、図示せぬディスプレイに表示すること等により、ユーザに提示して、情報提示処理は終了する。
【０２２７】
図１１は、図１０の情報提示処理で行われる消費電力U(m)の表示例を示す図である。
【０２２８】
データ出力部１６では、例えば、図１１に示すように、各ファクタ#mに対応する家電#mの消費電力U(m)の時系列が、その家電#mの名称等の家電ラベルL(m)とともに、図示せぬディスプレイに表示される。
【０２２９】
以上のように、監視システムでは、各ファクタが３以上の状態を有するFHMMによって、各家電の稼働状態をモデル化するFHMMの学習が行われるので、モードや設定等によって消費電力（電流）が変動するエアコン等の変動負荷家電について、正確な消費電力等を求めることができる。
【０２３０】
また、監視システムでは、配電盤等の1箇所で、家庭内の各家電で消費されている電流の総和を計測することにより、家庭内の各家電の消費電力を求めることができるので、各コンセントにスマートタップを取り付ける場合に比較して、コスト、及び、労力のいずれの面からも、容易に、家庭内の各家電の消費電力の「見える化」を実現することができる。
【０２３１】
そして、家庭内の各家電の消費電力の「見える化」によれば、例えば、家庭での節電の意識を高めることができる。
【０２３２】
なお、監視システムで得られる家庭内の各家電の消費電力については、例えば、サーバ等で収集し、各家庭の家電の消費電力から、家電の使用時間帯、ひいていは、生活パターンを推定し、マーケティング等に役立てることができる。
【０２３３】
［本技術を適用した監視システムの第２実施の形態］
【０２３４】
図１２は、本技術を適用した監視システムの第２実施の形態の構成例を示すブロックである。
【０２３５】
なお、図中、図３の場合と対応する部分については、同一の符号を付してあり、以下では、その説明は、適宜省略する。
【０２３６】
図１２の監視システムは、データ取得部１１ないしデータ出力部１６を有する点で、図３の監視システムと共通する。
【０２３７】
但し、図１２の監視システムは、モデル学習部１４の分散学習部３２に代えて、個別分散学習部５２が設けられている点で、図３の監視システムと相違する。
【０２３８】
図３の分散学習部３２では、分散学習において、全体モデルとしてのFHMMについて、1の分散Cが求められるが、図１２の個別分散学習部５２では、分散学習において、全体モデルとしてのFHMMについて、ファクタ#mごと、又は、各ファクタ#mの状態#kごとに、個別の分散（以下、個別分散ともいう）σ(m)が求められる。
【０２３９】
ここで、個別分散σ(m)が、各ファクタ#mの状態#kごとに個別の分散である場合、個別分散σ(m)は、例えば、ファクタ#mの状態#kの分散σ(m)kを、k列目のコンポーネントとするK列の行ベクトルである。
【０２４０】
なお、以下では、説明を簡単にするため、個別分散σ(m)が、ファクタ#mごとに個別の分散である場合、個別分散σ(m)としてのK列の行ベクトルのすべてのコンポーネントσ(m)1，σ(m)2，・・・，σ(m)Kを、ファクタ#mに共通の同一の値（ファクタ#mについての分散としての同一のスカラ値、又は、（共分散）行列）とすることで、ファクタ#mごとに個別の個別分散σ(m)を、各ファクタ#mの状態#kごとに個別の個別分散σ(m)と同一に扱うこととする。
【０２４１】
個別分散σ(m)としてのK列の行ベクトルのk列目のコンポーネントであるファクタ#mの状態#kの分散σ(m)kは、分散としてのスカラ値、又は、共分散行列としてのD行D列の行列である。
【０２４２】
個別分散学習部５２において、個別分散σ(m)が求められる場合、評価部２１では、その個別分散σ(m)を用い、式（４）に代えて、式（１６）に従って、観測確率P(Yt｜St)が求められる。
【０２４３】
【数１６】
0006020880.tif 000017
・・・（１６）
【０２４４】
ここで、式（１６）のΣtは、式（１７）に従って求められる。
【０２４５】
【数１７】
0006020880.tif 000018
・・・（１７）
【０２４６】
S*(m)tは、式（６）に示したように、時刻tにいるファクタ#mの状態を表し、K行のうちの1行のみのコンポーネントが0で、他の行のコンポーネントが0のK行の列ベクトルである。
【０２４７】
個別分散学習部５２において、個別分散σ(m)は、以下のようにして求められる。
【０２４８】
すなわち、実際に計測（観測）される計測波形Ytに対する全体モデルφとしてのFHMMの分散σの期待値<σ>は、式（１８）に示すように、計測波形Ytとしての電流波形と、全体モデルφとしてのFHMMにおいて生成（観測）される生成波形Y^tとしての電流波形との自乗誤差|Yt-Y^t|2の期待値<|Yt-Y^t|2>に等しい。
【０２４９】
【数１８】
0006020880.tif 000019
・・・（１８）
【０２５０】
分散σの期待値<σ>は、各ファクタ#mの状態#kごとに個別の個別分散σ(m)の期待値σ(m)<S(m)t>のファクタ#mについての総和σ(1)<S(1)t>+σ(2)<S(2)t>+・・・+σ(M)<S(M)t>に等価である。
【０２５１】
また、計測波形Ytと生成波形Y^tとの自乗誤差|Yt-Y^t|2の期待値<|Yt-Y^t|2>は、計測波形Ytと、個別波形W(m)の期待値W(m)<S(m)t>のファクタ#mについての総和W(1)<S(1)t>+W(2)<S(2)t>+・・・+W(M)<S(M)t>との自乗誤差|Yt-（W(1)<S(1)t>+W(2)<S(2)t>+・・・+W(M)<S(M)t>）|2に等価である。
【０２５２】
したがって、式（１８）は、式（１９）と等価である。
【０２５３】
【数１９】
0006020880.tif 000020
・・・（１９）
【０２５４】
式（１９）を満たす個別分散σ(m)＝σ(m)newは、個別分散σ(m)が負ではないという制約がある制約付き２次計画法により、式（２０）に従って求めることができる。
【０２５５】
【数２０】
0006020880.tif 000021
・・・（２０）
【０２５６】
個別分散学習部５２は、式（２０）に従って、個別分散σ(m)の更新値σ(m)newを求め、その更新値σ(m)newによって、モデル記憶部１３に記憶された個別分散σ(m)を更新する。
【０２５７】
以上のように、個別分散学習部５２において、ファクタ#mごと、又は、各ファクタ#mの状態#kごとの個別分散σ(m)を求める場合には、1の分散Cを用いる場合よりも、全体モデルφとしてのFHMMの、家電を表現する表現能力が向上するとともに、推定部２２での状態推定の精度が向上するので、より正確な消費電力等を求めることができる。
【０２５８】
特に、例えば、吸引の状態によって消費電流がばらつく掃除機等の変動負荷家電について、個別分散σ(m)を用いることにより、より正確な消費電力を求めることができる。
【０２５９】
図１２の監視システムでは、分散Cに代えて、個別分散σ(m)が用いられることを除き、図３の監視システムと同様の処理（図６の学習処理や、図１０の情報提示処理）が行われる。
【０２６０】
図１３は、図１２の監視システムが、図６のステップＳ１３で行うEステップの処理を説明するフローチャートである。
【０２６１】
ステップＳ６１において、評価部２１は、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの個別分散σ(m)、及び、固有波形W(m)、並びに、データ取得部１１からの計測波形Yt＝{Y1,Y2,・・・,YT}を用い、各時刻t＝{1,2,・・・,T}の、状態の各組み合わせStについて、（式（４）ではなく）式（１６）及び式（１７）に従って、観測確率P(Yt|St)を求め、推定部２２に供給して、処理は、ステップＳ６２に進む。
【０２６２】
以下、ステップＳ６２ないしＳ６６では、図７のステップＳ２２ないしＳ２６とそれぞれ同様の処理が行われる。
【０２６３】
図１４は、図１２の監視システムが、図６のステップＳ１４で行うMステップの処理を説明するフローチャートである。
【０２６４】
ステップＳ７１において、波形分離学習部３１は、図９のステップＳ３１と同様に、波形分離学習によって、固有波形W(m)の更新値W(m)newを求め、その更新値W(m)newによって、モデル記憶部１３に記憶された固有波形W(m)を更新して、処理は、ステップＳ７２に進む。
【０２６５】
ステップＳ７２では、個別分散学習部５２が、データ取得部１１からの計測波形Yt、推定部２２からの事後確率<S(m)t>、及び、及び、モデル記憶部１３に記憶された固有波形W(m)を用いて、分散学習としての式（２０）に従って、個別分散σ(m)の更新値σ(m)newを求める。
【０２６６】
そして、個別分散学習部５２は、、個別分散σ(m)の更新値σ(m)newによって、モデル記憶部１３に記憶された個別分散σ(m)を更新し、処理は、ステップＳ７２からステップＳ７３に進む。
【０２６７】
ステップＳ７３では、状態変動学習部３３は、図９のステップＳ３３と同様に、状態変動学習によって、遷移確率P(m)i,jの更新値P(m)i,jnew、及び、初期状態確率π(m)の更新値π(m)newを求め、その更新値P(m)i,jnew、及び、π(m)newによって、モデル記憶部１３に記憶された遷移確率P(m)i,j、及び、初期状態確率π(m)を更新して、Mステップの処理からリターンする。
【０２６８】
［本技術を適用した監視システムの第３実施の形態］
【０２６９】
図１５は、本技術を適用した監視システムの第３実施の形態の構成例を示すブロックである。
【０２７０】
なお、図中、図３の場合と対応する部分については、同一の符号を付してあり、以下では、その説明は、適宜省略する。
【０２７１】
図１５の監視システムは、データ取得部１１ないしデータ出力部１６を有する点で、図３の監視システムと共通する。
【０２７２】
但し、図１５の監視システムは、状態推定部１２の推定部２２に代えて、近似推定部４２が設けられている点で、図３の監視システムと相違する。
【０２７３】
近似推定部４２は、1時刻に状態が遷移するファクタの数を制限する状態遷移制限の下で、事後確率（状態確率）<S(m)t>等を求める。
【０２７４】
ここで、例えば、事後確率<S(m)t>は、式（９）に従い、事後確率γt,zを周辺化することで求められる。
【０２７５】
事後確率γt,zは、式（８）に従い、前向き確率αt,zと後ろ向き確率βt,zを用い、FHMMの各ファクタの状態の組み合わせzすべてについて求められる。
【０２７６】
ところで、FHMMの各ファクタの状態の組み合わせzの数は、KM通りだけ存在し、ファクタの数Mの増加に対して、指数のオーダで増加する。
【０２７７】
したがって、FHMMが、ある程度の数のファクタを有する場合に、各ファクタの状態の組み合わせzすべてについて、前向き確率αt,z、及び、後ろ向き確率βt,z、さらには、事後確率γt,zを、文献Ａに記載されているように厳密に計算するのでは、計算量が膨大になる。
【０２７８】
そこで、FHMMの状態推定については、例えば、Gibbs Samplingや、Completely Factorized Variational Approximation、Structured Variational Approximation等の近似推定手法が文献Ａにおいて提案されている。しかしながら、これらの近似推定手法では、計算量が依然として大きいことや、近似による精度の低下が著しいことがある。
【０２７９】
ところで、家庭内のすべての家電の稼働状態が、各時刻に変化（変動）する可能性は、停電等の異常なケースを除けば、極めて低い。
【０２８０】
そこで、近似推定部４２は、1時刻に状態が遷移するファクタの数を制限する状態遷移制限の下で、事後確率<S(m)t>等を求める。
【０２８１】
状態遷移制限としては、1時刻に状態が遷移するファクタの数を、例えば、1個や数個等以下に制限することを採用することができる。
【０２８２】
状態遷移制限によれば、事後確率<S(m)t>等を求めるのに用いられる事後確率γt,z、ひいては、前向き確率αt,z、及び、後ろ向き確率βt,zを、厳密に計算する必要がある状態の組み合わせzの数が、大幅に削減され、さらに、その削減される状態の組み合わせは、生じる可能性が極めて低い組み合わせであるため、事後確率（状態確率）<S(m)t>等の精度を大きく損なうことなく、計算量を大幅に少なくすることができる。
【０２８３】
近似推定部４２において、1時刻に状態が遷移するファクタの数を制限する状態遷移制限の下で、事後確率<S(m)t>等を求める方法としては、Eステップの処理に、パーティクルフィルタを適用する方法、すなわち、FHMMの状態の組み合わせzに、パーティクルフィルタを適用する方法がある。
【０２８４】
ここで、パーティクルフィルタについては、例えば、上述の文献Ｂの364ページに記載されている。
【０２８５】
また、本実施の形態では、状態遷移制限として、1時刻に状態が遷移するファクタの数を、例えば、1個以下に制限することを採用することとする。
【０２８６】
FHMMの状態の組み合わせzに、パーティクルフィルタを適用する場合、パーティクルフィルタのp番目のパーティクルSpは、FHMMの状態のある組み合わせを表現する。
【０２８７】
ここで、パーティクルSpは、そのパーティクルSpが表現する、ファクタ#1の状態S(1)p、ファクタ#2の状態S(2)p、・・・、ファクタ#Mの状態S(M)pの組み合わせであり、式Sp＝{S(1)p，S(2)p，・・・，S(M)p}で表現される。
【０２８８】
また、FHMMの状態の組み合わせzに、パーティクルフィルタを適用する場合、状態の組み合わせStにおいて、計測波形Ytが観測される式（４）の観測確率P(Yt｜St)に代えて、パーティクルSp（が表現する状態の組み合わせ）において、計測波形Ytが観測される観測確率P(Yt｜Sp)が用いられる。
【０２８９】
観測確率P(Yt｜Sp)は、例えば、式（２１）に従って計算することができる。
【０２９０】
【数２１】
0006020880.tif 000022
・・・（２１）
【０２９１】
式（２１）の観測確率P(Yt｜Sp)は、平均値（平均ベクトル）が、μtではなく、μpである点で、式（４）の観測確率P(Yt｜St)と異なっている。
【０２９２】
平均値μpは、平均値μtと同様に、電流波形Ytと同様のD行の列ベクトルであり、固有波形W(m)を用いて、式（２２）で表される。
【０２９３】
【数２２】
0006020880.tif 000023
・・・（２２）
【０２９４】
ここで、S*(m)pは、パーティクルSpとしての状態の組み合わせを構成する状態のうちの、ファクタ#mの状態を表し、以下、パーティクルSpのファクタ#mの状態S*(m)pともいう。パーティクルSpのファクタ#mの状態S*(m)pは、例えば、式（２３）に示すように、K行のうちの1行のみのコンポーネントが0で、他の行のコンポーネントが0のK行の列ベクトルである。
【０２９５】
【数２３】
0006020880.tif 000024
・・・（２３）
【０２９６】
パーティクルSpとしての状態の組み合わせを構成するファクタ#mの状態が、状態#kである場合、パーティクルSpのファクタ#mの状態S*(m)pとしてのK行の列ベクトルS*(m)pは、k行目のコンポーネントだけが1とされ、他のコンポーネントは0とされる。
【０２９７】
式（２２）によれば、パーティクルSpとしての状態の組み合わせを構成する各ファクタ#mの状態#kの固有波形W(m)kの総和が、時刻tの電流波形Ytの平均値μpとして求められる。
【０２９８】
図１６は、FHMMの状態の組み合わせzに、パーティクルフィルタを適用して、前向き確率αt,pを求める方法を説明する図である。
【０２９９】
なお、図１６では（後述する図１７及び図１８でも同様）、例えば、FHMMのファクタの数が4個であり、各ファクタが、オン状態とオフ状態とを稼働状態として表す2個の状態を有することとする。
【０３００】
評価部２１は、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの分散C、及び、固有波形W(m)、並びに、データ取得部１１からの計測波形Yt＝{Y1,Y2,・・・,YT}を用い、各時刻t＝{1,2,・・・,T}の、各パーティクルSp＝{S1,S2,・・・,SR}としての状態の組み合わせについて、式（２１）の観測確率P(Yt|Sp)を求め、近似推定部４２に供給する。
【０３０１】
そして、近似推定部４２は、評価部２１からの観測確率P(Yt|Sp)、及び、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの遷移確率P(m)i,j（及び初期状態確率π(m)）を用いて、計測波形Y1,Y2,・・・,Ytを観測し、時刻tに、パーティクルSpとしての状態の組み合わせにいる前向き確率αt,pを求める。
【０３０２】
前向き確率αt,pは、例えば、1時刻前の前向き確率αt-1,rを用いた漸化式αt,p＝Σαt-1,rP(Sp|r)P(Yt|Sp)に従って求めることができる。
【０３０３】
漸化式αt,p＝Σαt-1,rP(Sp|r)P(Yt|Sp)において、rは、1時刻に状態が遷移するファクタの数を1状態以下に制限する状態遷移制限の下で、時刻tのパーティクルSpとしての状態の組み合わせに遷移することができる時刻t-1のパーティクルとしての状態の組み合わせrを表す。
【０３０４】
したがって、状態の組み合わせrとなり得るのは、パーティクルSpとしての状態の組み合わせと同一の状態の組み合わせと、パーティクルSpとしての状態の組み合わせとは1個だけ状態が異なる状態の組み合わせとなる。
【０３０５】
漸化式αt,p＝Σαt-1,rP(Sp|r)P(Yt|Sp)において、Σは、状態の組み合わせrすべてについてとるサメーションを表す。
【０３０６】
また、漸化式αt,p＝Σαt-1,rP(Sp|r)P(Yt|Sp)において、P(Sp|r)は、時刻t-1にパーティクルとしての状態の組み合わせrにいて、時刻tにパーティクルSpとしての状態の組み合わせに遷移する遷移確率を表し、全体モデルφとしてのFHMMの遷移確率P(m)i,jを用い、式（３）に従って求められる。
【０３０７】
なお、前向き確率αt,pの初期値としては、例えば、状態の組み合わせすべての中から、パーティクルSpとして、ランダムに選択される状態の組み合わせを構成する各ファクタ#mの状態#kの初期状態確率π(m)kの積が採用される。
【０３０８】
近似推定部４２は、FHMMが取り得る状態の組み合わせのすべてではなく、パーティクルSpとしての状態の組み合わせについてだけ、前向き確率αt,pを求める。
【０３０９】
時刻tに存在するパーティクルSpの数をR個とすると、R個の前向き確率αt,pが求められる。
【０３１０】
近似推定部４２は、時刻tのR個のパーティクルSp＝{S1,S2,・・・,SR}について、前向き確率αt,pを求めると、その前向き確率αt,pに基づいて、R個のパーティクルSp＝{S1,S2,・・・,SR}の中からのパーティクルのサンプリングを行う。
【０３１１】
すなわち、近似推定部４２は、R個のパーティクルSp＝{S1,S2,・・・,SR}の中から、前向き確率αt,pの大きい順に、所定数であるP個のパーティクルSp＝{S1,S2,・・・,SP}をサンプリングし、そのP個のパーティクルSpだけを、時刻tのパーティクルとして残す。
【０３１２】
なお、P≧Rの場合は、R個のパーティクルSpすべてがサンプリングされる。
【０３１３】
また、パーティクルSpをサンプリングする個数Pは、FHMMが取り得る状態の組み合わせの数よりも小さい値であり、例えば、監視システムに許容される計算コスト等に基づいて設定される。
【０３１４】
近似推定部４２は、時刻tのP個のパーティクルSp＝{S1,S2,・・・,SP}をサンプリングすると、そのP個のパーティクルSp＝{S1,S2,・・・,SP}それぞれについて、1時刻後である時刻t+1のパーティクルを、状態遷移制限の下で予測する。
【０３１５】
1時刻に状態が遷移するファクタの数を1状態以下に制限する状態遷移制限の下では、時刻t+1のパーティクルSqとしては、パーティクルSpとしての状態の組み合わせと同一の状態の組み合わせと、パーティクルSpとしての状態の組み合わせとは1個だけ状態が異なる状態の組み合わせとが予測される。
【０３１６】
上述したように、FHMMのファクタの数が4個であり、各ファクタが、オン状態とオフ状態とを稼働状態として表す2個の状態を有する場合には、時刻tのあるパーティクルSpについては、その時刻tのパーティクルSpとしての状態の組み合わせと同一の1通りの状態の組み合わせと、時刻tのパーティクルSpとしての状態の組み合わせとは1個だけ状態が異なる4通りの状態の組み合わせとの、合計で5通りの組み合わせが、時刻t+1のパーティクルSqとして予測される。
【０３１７】
近似推定部４２は、時刻t+1のパーティクルSq＝{S1,S2,・・・,SQ}の予測後、その時刻t+1のパーティクルSqについて、時刻tのパーティクルSpの場合と同様に、前向き確率αt+1,qを求め、以下、同様にして、各時刻t=1,2,・・・,Tの前向き確率を求める。
【０３１８】
以上のように、近似推定部４２は、各ファクタ#mの状態の組み合わせを、パーティクルとして、1時刻後のパーティクルを、状態遷移制限の下で予測し、前向き確率αt,zに基づいて、所定数のパーティクルをサンプリングすることを繰り返しながら、各時刻tについて、パーティクルとしての状態の組み合わせzの前向き確率αt,zを求める。
【０３１９】
図１７は、FHMMの状態の組み合わせzに、パーティクルフィルタを適用して、後ろ向き確率βt,pを求める方法を説明する図である。
【０３２０】
近似推定部４２は、評価部２１からの観測確率P(Yt|St)、及び、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの遷移確率P(m)i,jを用いて、時刻tに、パーティクルSpとしての状態の組み合わせにいて、その後、計測波形Yt,Yt+1,・・・,YTを観測する後ろ向き確率βt,pを求める。
【０３２１】
後ろ向き確率βt,pは、例えば、1時刻後の後ろ向き確率βt+1,rを用いた漸化式βt,p＝ΣP(Yt|Sp)P(r|Sp)βt+1,rに従って求めることができる。
【０３２２】
漸化式βt,p＝ΣP(Yt|Sp)P(r|Sp)βt+1,rにおいて、rは、1時刻に状態が遷移するファクタの数を1状態以下に制限する状態遷移制限の下で、時刻tのパーティクルSpとしての状態の組み合わせから遷移することができる時刻t+1のパーティクルとしての状態の組み合わせrを表す。
【０３２３】
したがって、状態の組み合わせrとなり得るのは、パーティクルSpとしての状態の組み合わせと同一の状態の組み合わせと、パーティクルSpとしての状態の組み合わせとは1個だけ状態が異なる状態の組み合わせとなる。
【０３２４】
漸化式βt,p＝ΣP(Yt|Sp)P(r|Sp)βt+1,rにおいて、Σは、状態の組み合わせrすべてについてとるサメーションを表す。
【０３２５】
また、漸化式βt,p＝ΣP(Yt|Sp)P(r|Sp)βt+1,rにおいて、P(r|Sp)は、時刻tにパーティクルSpとしての状態の組み合わせにいて、時刻t+1にパーティクルとしての状態の組み合わせrに遷移する遷移確率を表し、全体モデルφとしてのFHMMの遷移確率P(m)i,jを用い、式（３）に従って求められる。
【０３２６】
なお、後ろ向き確率βt,pの初期値βT,pとしては、1が採用される。
【０３２７】
近似推定部４２は、前向き確率αt,pと同様に、FHMMが取り得る状態の組み合わせのすべてではなく、パーティクルSpとしての状態の組み合わせについてだけ、後ろ向き確率βt,pを求める。
【０３２８】
時刻tに存在するパーティクルSpの数をR個とすると、R個の後ろ向き確率βt,pが求められる。
【０３２９】
近似推定部４２は、時刻tのR個のパーティクルSp＝{S1,S2,・・・,SR}について、後ろ向き確率βt,pを求めると、その後ろ向き確率βt,pに基づいて、R個のパーティクルSp＝{S1,S2,・・・,SR}の中からのパーティクルのサンプリングを行う。
【０３３０】
すなわち、近似推定部４２は、前向き確率αt,pの場合と同様に、R個のパーティクルSp＝{S1,S2,・・・,SR}の中から、後ろ向き確率βt,pの大きい順に、所定数であるP個のパーティクルSp＝{S1,S2,・・・,SP}をサンプリングし、そのP個のパーティクルSpだけを、時刻tのパーティクルとして残す。
【０３３１】
近似推定部４２は、時刻tのP個のパーティクルSp＝{S1,S2,・・・,SP}をサンプリングすると、そのP個のパーティクルSp＝{S1,S2,・・・,SP}それぞれについて、1時刻前である時刻t-1のパーティクルを、状態遷移制限の下で予測する。
【０３３２】
1時刻に状態が遷移するファクタの数を1状態以下に制限する状態遷移制限の下では、時刻t-1のパーティクルSqとしては、パーティクルSpとしての状態の組み合わせと同一の状態の組み合わせと、時刻tのパーティクルSpとしての状態の組み合わせとは1個だけ状態が異なる状態の組み合わせとが予測される。
【０３３３】
上述したように、FHMMのファクタの数が4個であり、各ファクタが、オン状態とオフ状態とを稼働状態として表す2個の状態を有する場合には、時刻tのあるパーティクルSpについては、その時刻tのパーティクルSpとしての状態の組み合わせと同一の1通りの状態の組み合わせと、時刻tのパーティクルSpとしての状態の組み合わせとは1個だけ状態が異なる4通りの状態の組み合わせとの、合計で5通りの組み合わせが、時刻t-1のパーティクルSqとして予測される。
【０３３４】
近似推定部４２は、時刻t-1のパーティクルSq＝{S1,S2,・・・,SQ}について、時刻tのパーティクルSpの場合と同様に、後ろ向き確率βt+1,qを求め、以下、同様にして、各時刻t=T,T-1,・・・,1の後ろ向き確率を求める。
【０３３５】
以上のように、近似推定部４２は、各ファクタ#mの状態の組み合わせを、パーティクルとして、1時刻前のパーティクルを、状態遷移制限の下で予測し、後ろ向き確率βt,zに基づいて、所定数のパーティクルをサンプリングすることを繰り返しながら、各時刻tについて、パーティクルとしての状態の組み合わせzの後ろ向き確率βt,zを求める。
【０３３６】
図１８は、FHMMの状態の組み合わせzに、パーティクルフィルタを適用して、事後確率γt,pを求める方法を説明する図である。
【０３３７】
近似推定部４２は、各時刻tにおいて、前向き確率αt,pが求められたパーティクル（以下、前向きパーティクルともいう）Spとしての状態の組み合わせと、後ろ向き確率βt,p'が求められたパーティクル（以下、後ろ向きパーティクルともいう）Sp'としての状態の組み合わせとが、同一の状態の組み合わせzである場合、その状態の組み合わせzについて、前向き確率αt,zと後ろ向き確率βt,zを用いて、式（２３）に従い、全体モデルφとしてのFHMMにおいて、時刻tに、状態の組み合わせzにいる事後確率γt,zを求める。
【０３３８】
【数２４】
0006020880.tif 000025
・・・（２４）
【０３３９】
ここで、式（２４）の右辺の分母のΣは、wを、時刻tの前向きパーティクルSpと後ろ向きパーティクルSp'との両方に共通に残っている状態の組み合わせSp∩Sp'のすべてに変えてとるサメーションを表す。
【０３４０】
なお、時刻tに前向きパーティクルSpとして残っていない状態の組み合わせzの前向き確率αt,z、及び、時刻tに後ろ向きパーティクルSp'として残っていない状態の組み合わせzの後ろ向き確率βt,zは、例えば、いずれも、0とされる。
【０３４１】
したがって、時刻tの前向きパーティクルSpと後ろ向きパーティクルSp'との両方に共通に残っている状態の組み合わせSp∩Sp'以外の状態の組み合わせの事後確率γt,zは、0となる。
【０３４２】
図１５の監視システムでは、FHMMの状態の組み合わせzに、パーティクルフィルタを適用することにより、状態遷移制限の下で、事後確率<S(m)t>等を求めることを除き、図３の監視システムと同様の処理が行われる。
【０３４３】
図１９、及び、図２０は、図１５の監視システムが、図６のステップＳ１３で行うEステップの処理を説明するフローチャートである。
【０３４４】
なお、図２０は、図１９に続くフローチャートである。
【０３４５】
ステップＳ８１において、近似推定部４２は、時刻（を表す変数）tを、前向き確率についての初期値としての1に初期化して、処理は、ステップＳ８２に進む。
【０３４６】
ステップＳ８２では、近似推定部４２は、FHMMが取り得る状態の組み合わせの中から、所定の数の状態の組み合わせを、例えば、ランダムに、パーティクルSpとして選択（サンプリング）する。なお、ステップＳ８２では、FHMMが取り得る状態の組み合わせすべてを、パーティクルSpとして選択することができる。
【０３４７】
近似推定部４２は、パーティクルSpとしての状態の組み合わせについて、前向き確率αt,pの初期値α1,pを求め、処理は、ステップＳ８２からステップＳ８３に進む。
【０３４８】
ステップＳ８３では、時刻tが、計測波形Yt＝{Y1,Y2,・・・,YT}の系列の最後の時刻（系列長）Tに等しいかどうかを判定する。
【０３４９】
ステップＳ８３において、時刻tが、計測波形Ytの系列の最後の時刻Tに等しくないと判定された場合、すなわち、時刻tが時刻T未満である場合、処理は、ステップＳ８４に進み、近似推定部４２は、所定の条件としての状態遷移制限の下で、パーティクルSpとしての状態の組み合わせについて、1時刻後である時刻t+1のパーティクルを予測し、処理は、ステップＳ８５に進む。
【０３５０】
ステップＳ８５では、近似推定部４２は、時刻tを1だけインクリメントして、処理は、ステップＳ８６に進む。
【０３５１】
ステップＳ８６では、評価部２１は、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの分散C、及び、固有波形W(m)、並びに、データ取得部１１からの計測波形Yt＝{Y1,Y2,・・・,YT}を用い、時刻tの、各パーティクルSp＝{S1,S2,・・・,SR}としての状態の組み合わせについて、式（２１）の観測確率P(Yt|Sp)を求める。評価部２１は、観測確率P(Yt|Sp)を、近似推定部４２に供給して、処理は、ステップＳ８６からステップＳ８７に進む。
【０３５２】
ステップＳ８７では、近似推定部４２は、評価部２１からの観測確率P(Yt|Sp)、及び、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの遷移確率P(m)i,jを用いて、計測波形Y1,Y2,・・・,Ytを観測し、時刻tに、パーティクルSpとしての状態の組み合わせにいる前向き確率αt,pを求め、処理は、ステップＳ８８に進む。
【０３５３】
ステップＳ８８では、近似推定部４２は、前向き確率αt,pに基づき、パーティクルSpの中から、前向き確率αt,pの大きい順に、所定数のパーティクルをサンプリングし、時刻tのパーティクルとして残す。
【０３５４】
その後、処理は、ステップＳ８８からステップＳ８３に戻り、以下、同様の処理が繰り返される。
【０３５５】
そして、ステップＳ８３において、時刻tが、計測波形Ytの系列の最後の時刻Tに等しいと判定された場合、処理は、図２０のステップＳ９１に進む。
【０３５６】
ステップＳ９１では、近似推定部４２は、時刻tを、後ろ向き確率についての初期値としてのTに初期化して、処理は、ステップＳ９２に進む。
【０３５７】
ステップＳ９２では、近似推定部４２は、FHMMが取り得る状態の組み合わせの中から、所定の数の状態の組み合わせを、例えば、ランダムに、パーティクルSpとして選択する。なお、ステップＳ９２では、FHMMが取り得る状態の組み合わせすべてを、パーティクルSpとして選択することができる。
【０３５８】
近似推定部４２は、パーティクルSpとしての状態の組み合わせについて、後ろ向き確率βt,pの初期値αT,pを求め、処理は、ステップＳ９２からステップＳ９３に進む。
【０３５９】
ステップＳ９３では、時刻tが、計測波形Yt＝{Y1,Y2,・・・,YT}の系列の最初の時刻1に等しいかどうかを判定する。
【０３６０】
ステップＳ９３において、時刻tが、計測波形Ytの系列の最初の時刻1に等しくないと判定された場合、すなわち、時刻tが時刻1より大である場合、処理は、ステップＳ９４に進み、近似推定部４２は、所定の条件としての状態遷移制限の下で、パーティクルSpとしての状態の組み合わせについて、1時刻前である時刻t-1のパーティクルを予測し、処理は、ステップＳ９５に進む。
【０３６１】
ステップＳ９５では、近似推定部４２は、時刻tを1だけデクリメントして、処理は、ステップＳ９６に進む。
【０３６２】
ステップＳ９６では、評価部２１は、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの分散C、及び、固有波形W(m)、並びに、データ取得部１１からの計測波形Yt＝{Y1,Y2,・・・,YT}を用い、時刻tの、各パーティクルSp＝{S1,S2,・・・,SR}としての状態の組み合わせについて、式（２１）の観測確率P(Yt|Sp)を求める。評価部２１は、観測確率P(Yt|Sp)を、近似推定部４２に供給して、処理は、ステップＳ９６からステップＳ９７に進む。
【０３６３】
ステップＳ９７では、近似推定部４２は、評価部２１からの観測確率P(Yt|St)、及び、モデル記憶部１３に記憶された全体モデルφとしてのFHMMの遷移確率P(m)i,jを用いて、時刻tに、パーティクルSpとしての状態の組み合わせにいて、その後、計測波形Yt,Yt+1,・・・,YTを観測する後ろ向き確率βt,pを求め、処理は、ステップＳ９８に進む。
【０３６４】
ステップＳ９８では、近似推定部４２は、後ろ向き確率βt,pに基づき、パーティクルSpの中から、後ろ向き確率βt,pの大きい順に、所定数のパーティクルをサンプリングし、時刻tのパーティクルとして残す。
【０３６５】
その後、処理は、ステップＳ９８からステップＳ９３に戻り、以下、同様の処理が繰り返される。
【０３６６】
そして、ステップＳ９３において、時刻tが、計測波形Ytの系列の最初の時刻1に等しいと判定された場合、処理は、ステップＳ９９に進む。
【０３６７】
ステップＳ９９では、近似推定部４２は、前向き確率αt,z（＝αt,p)と後ろ向き確率βt,z（＝βt,p)を用いて、式（２４）に従い、全体モデルφとしてのFHMMにおいて、時刻tに、状態の組み合わせzにいる事後確率γt,zを求め、処理は、ステップＳ１００に進む。
【０３６８】
以下、ステップＳ１００及びＳ１０１では、図７のステップＳ２５及びＳ２６とそれぞれ同様の処理が行われ、事後確率<S(m)t>，<S(m)tS(n)t'>、及び、<S(m)t-1S(m)t>が求められる。
る。
【０３６９】
以上のように、近似推定部４２において、FHMMの状態の組み合わせzに、パーティクルフィルタを適用することで、状態遷移制限の下で、事後確率<S(m)t>等を求めることにより、可能性が低い状態の組み合わせについての計算を省略し、事後確率<S(m)t>等の計算の効率化を図ることができる。
【０３７０】
［本技術を適用した監視システムの第４実施の形態］
【０３７１】
図２１は、本技術を適用した監視システムの第４実施の形態の構成例を示すブロックである。
【０３７２】
なお、図中、図３の場合と対応する部分については、同一の符号を付してあり、以下では、その説明は、適宜省略する。
【０３７３】
図２１の監視システムは、データ取得部１１ないしデータ出力部１６を有する点で、図３の監視システムと共通する。
【０３７４】
但し、図２１の監視システムは、モデル学習部１４の波形分離学習部３１に代えて、制約付き分離学習部５１が設けられている点で、図３の監視システムと相違する。
【０３７５】
図３の分離学習部３１では、特に制約なしで、固有波形W(m)が求められるが、制約付き波形分離学習部５１は、家電（家電分離）に特有の所定の制約の下で、固有波形W(m)を求める。
【０３７６】
ここで、図３の分離学習部３１では、Mステップにおいて、固有波形W(m)の算出にあたり、特に、制約は課されていない。
【０３７７】
そのため、図３の分離学習部３１では、計測波形Ytが、複数の電流波形を重畳した波形であることのみを考慮して、波形分離学習により、固有波形W(m)が算出される。
【０３７８】
すなわち、図３の分離学習部３１の波形分離学習では、計測波形Ytが、複数の電流波形を重畳した波形であるとして、その複数の電流波形としての固有波形W(m)に分離される。
【０３７９】
固有波形W(m)の算出に、何らの制約もない場合、波形分離学習の解の自由度（冗長性）、すなわち、計測波形Ytから分離することができる複数の電流波形の自由度は、非常に高ため、家電の電流波形としては不適切な電流波形としての固有波形W(m)が求められることがあり得る。
【０３８０】
そこで、制約付き波形分離学習部５１では、家電に特有の所定の制約の下で、波形分離学習を行うことで、家電の電流波形としては不適切な電流波形としての固有波形W(m)が求められることを防止し、家電の正確な電流波形としての固有波形W(m)を求める。
【０３８１】
ここで、所定の制約の下で行われる波形分離学習を、制約付き波形分離学習ともいう。
【０３８２】
図２１の監視システムでは、制約付き波形分離学習が行われることを除き、図３の監視システムと同様の処理が行われる。
【０３８３】
制約付き波形分離学習において、家電に特有の制約としては、例えば、負荷制約や基底波形制約がある。
【０３８４】
負荷制約は、家電#mの電流波形としての固有波形W(m)と、家電#mにかかる電圧の電圧波形、すなわち、電流波形である計測波形Ytに対応する電圧波形Vtとの乗算により求められる家電の消費電力U(m)が、負の値にならない（家電#mが発電をしない）という制約である。
【０３８５】
基底波形制約は、家電#mの各稼働状態の消費電流の電流波形である固有波形W(m)は、その家電#mについて、基底波形としてあらかじめ用意された複数の波形の1個以上の組み合わせで表現されるという制約である。
【０３８６】
図２２は、図２１の監視システムが、負荷制約を課して行う図６のステップＳ１４のMステップの処理を説明するフローチャートである。
【０３８７】
ステップＳ１２１において、制約付き波形分離学習部５１は、データ取得部１１からの計測波形Yt、並びに、推定部２２からの事後確率<S(m)t>、及び、<S(m)tS(n)t'>を用い、負荷制約の下で、波形分離学習を行うことにより、固有波形W(m)の更新値W(m)newを求め、その更新値W(m)newによって、モデル記憶部１３に記憶された固有波形W(m)を更新して、処理は、ステップＳ１２２に進む。
【０３８８】
すなわち、制約付き波形分離学習部５１は、負荷制約の下での波形分離学習として、式（２５）を、制約付きの２次計画法を用いて解くことにより、固有波形W(m)の更新値W(m)newを求める。
【０３８９】
【数２５】
0006020880.tif 000026
・・・（２５）
【０３９０】
ここで、式（２５）では、固有波形W(m)の更新値W(m)newとなる固有波形Wに、負荷制約としての式0≦V'Wを満たす制約が課されている。
【０３９１】
Vは、計測波形Ytとしての電流波形に対応する電圧波形Vtを表すD行の列ベクトルであり、V'は、列ベクトルVを転置した行ベクトルである。
【０３９２】
なお、式（２５）の固有波形の更新値Wnew、事後確率<St'>、及び、<StSt'>は、式（１２）で説明したそれと一致する。
【０３９３】
すなわち、固有波形の更新値Wnewは、D行K×M列の行列であり、その(m-1)K+k列の列ベクトルは、ファクタ#mの状態#kの固有波形W(m)kの更新値になっている。
【０３９４】
事後確率<St'>は、K×M列の行ベクトルであり、その(m-1)K+k列目のコンポーネントは、時刻tに、ファクタ#mの状態#kにいる状態確率になっている。
【０３９５】
事後確率<StSt'>は、K×M行K×M列の行列であり、その(m-1)K+k行(n-1)K+k'列目のコンポーネントは、時刻tに、ファクタ#mの状態#kと、他のファクタ#nの状態#k'とにいる状態確率になっている。
【０３９６】
なお、式（２５）の固有波形Wは、固有波形の更新値Wnewと同様のD行K×M列の行列であり、式（２５）の右辺のargmin{X}のXを最小にする固有波形Wが、固有波形の更新値Wnewとして、２次計画法により求められる。
【０３９７】
ステップＳ１２２及びＳ１２３では、図９のステップＳ３２及びＳ３３と同様の処理が行われる。
【０３９８】
図２３は、負荷制約を説明する図である。
【０３９９】
すなわち、図２３Ａは、負荷制約のない波形分離学習により求められる電流波形としての固有波形W(m)を示しており、図２３Ｂは、負荷制約の下で行われた波形分離学習により求められる電流波形としての固有波形W(m)を示している。
【０４００】
なお、図２３では、ファクタの数Mを2個として、ファクタ#1のある状態#kの固有波形W(1)kと、ファクタ#2のある状態#k'の固有波形W(2)k'とを示してある。
【０４０１】
図２３Ａ及び図２３Ｂでは、ファクタ#1の固有波形W(1)kである電流波形は、いずれも、電圧波形Vtと同相になっている。
【０４０２】
しかしながら、図２３Ａでは、負荷制約がされていないため、ファクタ#2の固有波形W(2)k'である電流波形は、電圧波形Vtと逆相になっている。
【０４０３】
一方、図２３Ｂでは、負荷制約が課された結果、ファクタ#2の固有波形W(2)k'である電流波形は、電圧波形Vtと同相になっている。
【０４０４】
計測波形Ytは、図２３Ａに示したファクタ#1の固有波形W(1)kと、ファクタ#2の固有波形W(2)k'とに分離することもできるし、図２３Ｂに示したファクタ#1の固有波形W(1)kと、ファクタ#2の固有波形W(2)k'とに分離することもできる。
【０４０５】
しかしながら、計測波形Ytを、図２３Ａに示したファクタ#1の固有波形W(1)kと、ファクタ#2の固有波形W(2)k'とに分離した場合には、固有波形W(2)k'が電圧波形Vtと逆相になっているファクタ#2に対応する家電#2は、発電を行っていることになり、適切な家電分離を行うことが困難となるおそれがある。
【０４０６】
一方、負荷制約を課した場合には、計測波形Ytは、図２３Ｂに示したファクタ#1の固有波形W(1)kと、ファクタ#2の固有波形W(2)k'とに分離される。
【０４０７】
図２３Ｂの固有波形W(1)k、及び、固有波形W(2)k'は、いずれも、電圧波形Vtと同相になっており、したがって、負荷制約によれば、ファクタ#1に対応する家電#1、及び、ファクタ#2に対応する家電#2が、いずれも、電力を消費する負荷であるとして、適切な家電分離を行うことができる。
【０４０８】
図２４は、基底波形制約を説明する図である。
【０４０９】
図２４において、Yは、D行の列ベクトルである計測波形Ytが、時刻t順に、左から右方向に並んだD行T列の行列である。D行T列の行列である計測波形Yのt列目の列ベクトルは、時刻tの計測波形Ytになっている。
【０４１０】
図２４において、Wは、D行K列の行列であるファクタ#mの固有波形W(m)が、ファクタ#m順に、左から右方向に並んだD行K×M列の行列である。D行K×M列の行列である固有波形Wの(m-1)K+k列の列ベクトルは、ファクタ#mの状態#kの固有波形W(m)kになっている。
【０４１１】
図２４において、Fは、K×Mを表す。
【０４１２】
図２４において、<St>は、K行の列ベクトルである時刻tの事後確率<S(m)t>がファクタ#m順に上から下方向に並んだF=K×M行の列ベクトルが、時刻t順に、左から右方向に並んだをK×M行T列の行列である。K×M行T列の行列である事後確率<St>の(m-1)K+k行t列目のコンポーネントは、時刻tに、ファクタ#mの状態#kにいる状態確率になっている。
【０４１３】
波形分離学習では、FHMMにおいて、図２４に示すように、固有波形Wと事後確率<St>との積W×<St>が、計測波形Yとして観測されることとして、固有波形Wが求められる。
【０４１４】
基底波形制約は、上述したように、家電#mの各稼働状態の消費電流の電流波形である固有波形W(m)が、その家電#mについて、基底波形としてあらかじめ用意された複数の波形の1個以上の組み合わせで表現されるという制約であるため、基底波形制約を課した波形分離学習では、図２４に示すように、固有波形Wが、所定の数Nの基底波形Bと、所定の係数Aとの積B×Aで表されることとして、固有波形Wが求められる。
【０４１５】
ここで、所定の数Nの基底波形Bのうちの、n番目の基底波形Bを、Bnと表すこととすると、Bnは、例えば、波形のサンプル値をコンポーネントとするD行の列ベクトルであり、基底波形Bは、D行の列ベクトルである基底波形Bnが、インデクスn順に、左から右方向に並んだD行N列の行列である。
【０４１６】
係数Aは、N行K×M列の行列であり、n行(m-1)K+k列目のコンポーネントは、固有波形W(m)kを、N個の基底波形B1,B2,・・・,BNの組み合わせ（重畳）で表現するのに、n番目の基底波形Bnに乗算される係数になっている。
【０４１７】
ここで、固有波形W(m)kを、N個の基底波形B1ないしBNの組み合わせで表現するのに、N個の基底波形B1ないしBNと乗算される係数としての、例えば、N行の列ベクトルを、係数A(m)kと表し、その係数A(m)kを、ファクタ#mのK個の状態#1ないし#Kの順に、左から右方向に並べたN行K列の行列を、係数A(m)と表すこととすると、係数Aは、係数A(m)を、ファクタ#m順に、左から右方向に並べた行列である。
【０４１８】
基底波形Bは、例えば、画像処理等で用いられるICA(Independent Component Analysis)やNMF(Non-negative Matrix Factorization)等の基底抽出を、計測波形Yを対象に行うことにより用意（取得）することができる。
【０４１９】
また、メーカ等が、家電の基底波形を、サイト等で公開している場合には、基底波形Bは、サイトにアクセスすることにより用意することができる。
【０４２０】
図２５は、図２１の監視システムが、基底波形制約を課して行う図６のステップＳ１４のMステップの処理を説明するフローチャートである。
【０４２１】
ステップＳ１３１において、制約付き波形分離学習部５１は、データ取得部１１からの計測波形Yt、並びに、推定部２２からの事後確率<S(m)t>、及び、<S(m)tS(n)t'>を用い、基底波形制約の下で、波形分離学習を行うことにより、固有波形W(m)の更新値W(m)newを求め、その更新値W(m)newによって、モデル記憶部１３に記憶された固有波形W(m)を更新して、処理は、ステップＳ１３２に進む。
【０４２２】
すなわち、制約付き波形分離学習部５１は、基底波形制約の下での波形分離学習として、式（２６）を、制約付きの２次計画法を用いて解くことにより、固有波形Wを、基底波形Bの組み合わせで表現するための係数Anewを求める。
【０４２３】
【数２６】
0006020880.tif 000027
・・・（２６）
【０４２４】
ここで、式（２６）では、固有波形Wが、所定の数Nの基底波形Bと、所定の係数Aとの積B×Aで表されること、及び、係数Aを最小化すること(minA)が、基底波形制約として課されている。
【０４２５】
係数Aを最小化すること(minA)とは、係数AであるN行K×M列の行列の各コンポーネントの値（大きさ）を最小化することを意味し、ひいては、係数AであるN行K×M列の行列を、なるべくスパースにすることを意味する。
【０４２６】
係数AであるN行K×M列の行列を、なるべくスパースにすることにより、固有波形W(m)kは、なるべく少ない数の基底波形Bnの組み合わせで表されることとなる。
【０４２７】
制約付き波形分離学習部５１は、式（２６）に従って、係数Anewを求めると、式（２７）に従い、係数Anewと基底波形Bとを用いて、固有波形W(m)の更新値W(m)newを求める。
【０４２８】
【数２７】
0006020880.tif 000028
・・・（２７）
【０４２９】
ステップＳ１３２及びＳ１３３では、図９のステップＳ３２及びＳ３３と同様の処理が行われる。
【０４３０】
家電で消費される消費電流が、基底波形の1以上の組み合わせの波形になる場合には、基底波形制約を課すことにより、家電で消費される消費電流としてあり得ない電流の波形が、固有波形として求められることを防止して、家電に適切な固有波形を求めることができる。
【０４３１】
なお、以上においては、負荷制約と、基底波形制約とを、別個に課すこととしたが、負荷制約と、基底波形制約とは、同時に課すことができる。
【０４３２】
負荷制約と基底波形制約とを同時に課す場合には、制約付き波形分離学習部５１は、Mステップにおいて、負荷制約、及び、基底波形制約の下での波形分離学習として、式（２８）を、制約付きの２次計画法を用いて解くことにより、固有波形Wを、基底波形Bの組み合わせで表現するための係数Anewであって、固有波形Wを用いて求められる家電の消費電力を負の値にしない係数Anewを求める。
【０４３３】
【数２８】
0006020880.tif 000029
・・・（２８）
【０４３４】
そして、制約付き波形分離学習部５１は、式（２８）に従って求められた係数Anewと、基底波形Bとを用い、式（２９）に従って、固有波形W(m)の更新値W(m)newを求める。
【０４３５】
【数２９】
0006020880.tif 000030
・・・（２９）
【０４３６】
［本技術を適用した監視システムの第５実施の形態］
【０４３７】
図２６は、本技術を適用した監視システムの第５実施の形態の構成例を示すブロックである。
【０４３８】
なお、図中、図３の場合と対応する部分については、同一の符号を付してあり、以下では、その説明は、適宜省略する。
【０４３９】
図２６の監視システムは、データ取得部１１ないしデータ出力部１６を有する点で、図３の監視システムと共通する。
【０４４０】
但し、第１に、図２６の監視システムは、モデル学習部１４の分散学習部３２に代えて、図１２の個別分散学習部５２が設けられている点で、図３の監視システムと相違する。
【０４４１】
第２に、図２６の監視システムは、状態推定部１２の推定部２２に代えて、図１５の近似推定部４２が設けられている点で、図３の監視システムと相違する。
【０４４２】
第３に、図２６の監視システムは、モデル学習部１４の波形分離学習部３１に代えて、図２１の制約付き分離学習部５１が設けられている点で、図３の監視システムと相違する。
【０４４３】
したがって、図２６の監視システムでは、ファクタ#mごと、又は、各ファクタ#mの状態#kごとに、個別の分散（個別分散）σ(m)を求める個別分散学習が行われる。
【０４４４】
さらに、図２６の監視システムでは、1時刻に状態が遷移するファクタの数を制限する状態遷移制限の下で、事後確率（状態確率）<S(m)t>等を求める近似推定が行われる。
【０４４５】
また、図２６の監視システムでは、制約付き波形分離学習により、固有波形W(m)が求められる。
【０４４６】
なお、監視システムでは、個別分散学習、近似推定、及び、制約付き波形分離学習のそれぞれを単独で行うことや、個別分散学習、近似推定、及び、制約付き波形分離学習のすべてを行うことの他、個別分散学習、近似推定、及び、制約付き波形分離学習のうちの任意の2つを行うことが可能である。
【０４４７】
［家電分離以外への監視システムの適用］
【０４４８】
以上、FHMMの学習を行う監視システムにおいて、総和データとしての電流波形を監視（モニタ）して、家電分離を行う場合について説明したが、FHMMの学習を行う監視システムは、その他、1以上の信号が重畳されている重畳信号を監視して、その重畳信号に重畳されている信号を分離する任意のアプリケーションに適用することができる。
【０４４９】
図２７は、FHMMの学習を行う監視システムによる話者分離の概要を説明する図である。
【０４５０】
FHMMの学習を行う監視システムによれば、家電分離で計測波形Ytとして用いた電流波形に代えて、複数の話者の発話が重畳されている音声信号を用いることで、その複数の話者の発話が重畳されている音声信号から、各話者の音声を分離する話者分離を行うことができる。
【０４５１】
FHMMの学習を行う監視システムにおいて、話者分離を行う場合にも、家電分離の場合と同様に、個別分散学習や、近似推定、制約付き波形分離学習を行うことができる。
【０４５２】
話者分離において、ファクタ#mごと、又は、各ファクタ#mの状態#kごとの個別分散σ(m)を用いる場合には、1の分散Cを用いる場合よりも、FHMMの表現能力が向上するので、話者分離の精度を向上させることができる。
【０４５３】
また、複数の話者のすべての発話状態が、各時刻に変化する可能性は極めて低いことから、1時刻に、発話状態が変化する話者を、１人や数人以下に制限しても、話者分離の精度に、大きな影響はない。
【０４５４】
そこで、話者分離では、1時刻に状態が遷移するファクタの数を、例えば、1個以下に制限する（発話状態が変化する話者を１人以下に制限する）状態遷移制限の下で、事後確率<S(m)t>等を求める近似推定を行うことができる。そして、状態遷移制限の下で、事後確率<S(m)t>等を求める近似推定によれば、事後確率<S(m)t>等を求めるのに用いられる事後確率γt,zを、厳密に計算する必要がある状態の組み合わせzの数が、大幅に削減されるので、計算量を大幅に少なくすることができる。
【０４５５】
さらに、話者分離において、個別波形W(m)は、人の音声の波形であるから、その周波数成分は、人の音声が取り得る周波数帯域に存在する。したがって、話者分離で行う制約付き波形分離学習では、個別波形W(m)の周波数成分を、人の音声が取り得る周波数帯域内の周波数成分に制限するという、人の音声に特有の制約を採用することができる。この場合、人の音声の波形として、適切な固有波形W(m)を求めることができる。
【０４５６】
なお、話者分離において、制約付き波形分離学習で課す制約としては、その他、例えば、家電分離の場合と同様に、基底波形制約を採用することができる。
【０４５７】
［本技術を適用した監視システムの第６実施の形態］
【０４５８】
図２８は、本技術を適用した監視システムの第６実施の形態の構成例を示すブロックである。
【０４５９】
なお、図中、図３の場合と対応する部分については、同一の符号を付してあり、以下では、その説明は、適宜省略する。
【０４６０】
図２８の監視システムは、データ取得部１１ないしデータ出力部１６を有する点で、図３の監視システムと共通する。
【０４６１】
但し、第１に、図２８の監視システムは、状態推定部１２において、評価部２１及び推定部２２に代えて、評価部７１及び推定部７２が、それぞれ設けられている点で、図３の監視システムと相違する。
【０４６２】
第２に、図２８の監視システムは、モデル学習部１４において、波形分離学習部３１に代えて、制約付き分離学習部８１が設けられており、分散学習部３２及び状態変動学習部３３が設けられていない点で、図３の監視システムと相違する。
【０４６３】
図３の監視システムでは、FHMMを、全体モデルφとして用いて、家電分離を行うが、図２８の監視システムでは、FHMMではないモデル、すなわち、例えば、固有波形W(m)のみをモデルパラメータとして有するモデル（以下、波形モデルともいう）を、全体モデルφとして用いて、家電分離を行う。
【０４６４】
したがって、図２８では、モデル記憶部１３には、波形モデルが、全体モデルφとして記憶される。ここで、波形モデルは、固有波形W(m)をモデルパラメータとして有するが、かかる波形モデルにおいては、1つの固有波形W(m)が、家電モデル#mに相当する。
【０４６５】
状態推定部１２は、評価部７１及び推定部７２を有し、データ取得部１１から供給される総和データとしての計測波形Yt、及び、モデル記憶部１３に記憶された波形モデルを用いて、複数としてのM個の家電#1,#2,・・・,#mの稼働状態を求める状態推定を行う。
【０４６６】
すなわち、評価部７１は、モデル記憶部１３に記憶された全体モデルφとしての波形モデルを構成する各家電モデル#mにおいて、データ取得部１１から供給される電流波形Yが観測される程度を評価した評価値Eを求め、推定部７２に供給する。
【０４６７】
推定部７２は、評価部７１から供給される評価値Eを用い、各家電モデル#mが表す家電の各時刻tの稼働状態C(m)t,kを、例えば、整数計画法によって推定し、モデル学習部１４、ラベル取得部１５、及び、データ出力部１６に供給する。
【０４６８】
ここで、推定部７２は、例えば、式（３０）の整数計画問題を、整数計画法によって解くことにおり、家電#mの稼働状態C(m)t,kを推定する。
【０４６９】
【数３０】
0006020880.tif 000031
・・・（３０）
【０４７０】
ここで、式（３０）において、Eは、計測波形Ytと、全体モデルφとしての波形モデルにおいて観測される総和データとしての電流波形ΣΣW(m)kC(m)t,kとの誤差を表し、推定部７２では、この誤差Eを最小化する(minimize)稼働状態C(m)t,kが求められる。
【０４７１】
また、式（３０）において、固有波形W(m)kは、家電#mの稼働状態C(m)t,kに固有の電流の波形である固有波形を表し、計測波形Ytと同様に、D行の列ベクトルである。
【０４７２】
さらに、式（３０）において、K(m)は、家電#mの稼働状態C(m)t,kの数（種類数）を表す。
【０４７３】
稼働状態C(m)t,kは、0以上の整数で、スカラ値であり、時刻tの家電#mの稼働状態を表す。稼働状態C(m)t,kは、例えば、家電#mのモード（設定）に相当し、家電#mにおいて、オンとなることができるモードが、例えば、1個以下に制限されていることとする。
【０４７４】
家電#mにおいて、オンとなることができるモードが1個以下に制限されていることは、式（３１）で表される。
【０４７５】
【数３１】
0006020880.tif 000032
・・・（３１）
【０４７６】
式（３１）によれば、0以上の整数である稼働状態C(m)t,kが取り得る値は、0又は1となる。
【０４７７】
なお、家電の稼働状態を、整数計画法によって推定することについては、例えば、「電気機器の非侵入型稼働状態モニタリングシステム」、稲垣伸吉、江上司、鈴木達也(名古屋大学)、中村久栄、伊藤公一(（株）トーエネック)、第42回計測自動制御学会離散事象システム研究会，pp.33-38，Dec 20，2008，大阪大学に記載されている。
【０４７８】
以上のように、推定部７２が、式（３０）を解くことにより、家電#mの稼働状態C(m)t,kを推定する場合、評価部７１は、データ取得部１１から供給される計測波形Yt、及び、モデル記憶部１３に記憶された波形モデルの固有波形W(m)kを用いて、式（３０）の誤差Eを求め、評価値Eとして、推定部７２に供給する。
【０４７９】
モデル学習部１４において、制約付き波形分離学習部８１は、図２１の制約付き波形分離学習部５１と同様に、家電に特有の所定の制約の下での波形分離学習（制約付き波形分離学習）を行うことで、家電の電流波形としては不適切な電流波形としての固有波形W(m)kが求められることを防止し、家電の正確な電流波形としての固有波形W(m)kを求める。
【０４８０】
すなわち、制約付き波形分離学習部８１は、データ取得部１１から供給される計測波形Yt、並びに、推定部７２から供給される家電#mの稼働状態C(m)t,kを用い、例えば、負荷制約の下で、波形分離学習を行うことにより、固有波形W(m)kの更新値W(m)newkを求め、その更新値W(m)newkによって、モデル記憶部１３に記憶された固有波形W(m)kを更新する。
【０４８１】
すなわち、制約付き波形分離学習部８１は、負荷制約の下での波形分離学習として、例えば、式（３２）の２次計画問題を、２次計画法によって解くことにより、固有波形W(m)kの更新値W(m)newkを求める。
【０４８２】
【数３２】
0006020880.tif 000033
・・・（３２）
【０４８３】
ここで、式（３２）では、固有波形W(m)kの更新値W(m)newkとなる固有波形Wに、負荷制約としての式0≦V'Wを満たす制約が課されている。
【０４８４】
式（３２）において、Vは、計測波形Ytとしての電流波形に対応する電圧波形Vtを表すD行の列ベクトルであり、V'は、列ベクトルVを転置した行ベクトルである。
【０４８５】
また、稼働状態C(m)t,kがオン、オフを表す家電#mのモードを、モード#kということとすると、式（３２）において、Wは、(m-1)K+k列の列ベクトルが、家電#mのモード#k（がオンのとき）の固有波形になっているD行K×M列の行列である。なお、Kは、ここでは、例えば、K(1),K(2),・・・,K(M)の中の最大値を表すこととする。
【０４８６】
固有波形の更新値Wnewは、D行K×M列の行列であり、その(m-1)K+k列の列ベクトルは、ファクタ#mの状態#kの固有波形W(m)kの更新値になっている。
【０４８７】
また、式（３２）において、Ctは、(m-1)K+k行のコンポーネントが、稼働状態C(m)t,kになっているK×M行の列ベクトルである。
【０４８８】
なお、制約付き波形分離学習部８１では、その他の制約を課すことができる。すなわち、制約付き波形分離学習部８１では、例えば、図２１の制約付き波形分離学習部５１と同様に、負荷制約に代えて、基底波形制約を課すことや、負荷制約及び基底波形制約の両方を課すことができる。
【０４８９】
図２９は、図２８の監視システムが行う波形モデルの学習の処理（学習処理）を説明するフローチャートである。
【０４９０】
ステップＳ１５１において、モデル学習部１４は、モデル記憶部１３に記憶された全体モデルのモデルパラメータφとしての固有波形W(m)kを初期化し、処理は、ステップＳ１５２に進む。
【０４９１】
ここで、固有波形W(m)kとしてのD行の列ベクトルの各コンポーネントは、例えば、乱数を用いて初期化される。
【０４９２】
ステップＳ１５２では、データ取得部１１が、所定の時間T分の電流波形を取得し、各時刻t=1,2・・・,Tの電流波形を、計測波形Y1,Y2,・・・,YTとして、状態推定部１２、及び、モデル学習部１４に供給して、処理は、ステップＳ１５３に進む。
【０４９３】
ここで、データ取得部１１は、時刻t=1,2・・・,Tの電流波形とともに、電圧波形も取得する。データ取得部１１は、時刻t=1,2・・・,Tの電圧波形を、データ出力部１６に供給する。
【０４９４】
データ出力部１６では、情報提示処理（図１０）において、データ取得部１１からの電圧波形が、消費電力の算出に用いられる。
【０４９５】
ステップＳ１５３では、状態推定部１２の評価部７１が、データ取得部１１からの計測波形Y1ないしYT、及び、モデル記憶部１３に記憶された波形モデルの固有波形W(m)kを用いて、各家電#mの稼働状態C(m)t,kを求めるための式（３０）の評価値Eとしての誤差Eを求める。
【０４９６】
さらに、評価部７１は、誤差Eを、推定部７２に供給し、処理は、ステップＳ１５３からステップＳ１５４に進む。
【０４９７】
ステップＳ１５４では、推定部７２は、評価部７１からの式（３０）の誤差Eを最小化することにより、各家電#mの稼働状態C(m)t,kを推定し、モデル学習部１４、ラベル取得部１５、及び、データ出力部１６に供給して、処理は、ステップＳ１５５に進む。
【０４９８】
ステップＳ１５５では、モデル学習部１４の制約付き波形分離学習部８１は、データ取得部１１から供給される計測波形Yt、並びに、推定部７２から供給される家電#mの稼働状態C(m)t,kを用い、負荷制約等の所定の制約の下で、波形分離学習を行うことにより、固有波形W(m)kの更新値W(m)newkを求める。
【０４９９】
すなわち、制約付き波形分離学習部８１は、例えば、負荷制約の下での波形分離学習として、式（３２）を解くことにより、固有波形W(m)kの更新値W(m)newkを求める。
【０５００】
そして、制約付き波形分離学習部８１は、固有波形W(m)kの更新値W(m)newkによって、モデル記憶部１３に記憶された固有波形W(m)kを更新し、処理は、ステップＳ１５５からステップＳ１５６に進む。
【０５０１】
ステップＳ１５６では、モデル学習部１４は、モデルパラメータφの収束条件が満たされているかどうかを判定する。
【０５０２】
ここで、モデルパラメータφの収束条件としては、例えば、ステップＳ１５４の稼働状態C(m)t,kの推定と、ステップＳ１５５の制約付き波形分離学習による固有波形W(m)kの更新とがあらかじめ設定された所定の回数だけ繰り返されたことや、式（３０）の誤差Eの、モデルパラメータφの更新前と更新後との間の変化量が、あらかじめ設定された閾値以内であること等を採用することができる。
【０５０３】
ステップＳ１５６で、モデルパラメータφの収束条件が満たされていないと判定された場合、処理は、ステップＳ１５３に戻り、以下、ステップＳ１５３ないしＳ１５６の処理が繰り返される。
【０５０４】
ステップＳ１５３ないしＳ１５６の処理が繰り返し行われることにより、ステップＳ１５４の稼働状態C(m)t,kの推定と、ステップＳ１５５の制約付き波形分離学習による固有波形W(m)kの更新とが交互に繰り返され、ステップＳ１５４で求められる稼働状態C(m)t,kと、ステップＳ１５５で求められる固有波形W(m)k（の更新値）とは、精度が増していく。
【０５０５】
その後、ステップＳ１５６において、モデルパラメータφの収束条件が満たされると判定された場合、学習処理は、終了する。
【０５０６】
なお、状態推定部１２では、整数計画法や、FHMMを用いない任意のその他の方法で、家電の稼動状態を推定することができる。
【０５０７】
［本技術を適用したコンピュータの説明］
【０５０８】
次に、上述した一連の処理は、ハードウェアにより行うこともできるし、ソフトウェアにより行うこともできる。一連の処理をソフトウェアによって行う場合には、そのソフトウェアを構成するプログラムが、汎用のコンピュータ等にインストールされる。
【０５０９】
そこで、図３０は、上述した一連の処理を実行するプログラムがインストールされるコンピュータの一実施の形態の構成例を示している。
【０５１０】
プログラムは、コンピュータに内蔵されている記録媒体としてのハードディスク１０５やROM１０３に予め記録しておくことができる。
【０５１１】
あるいはまた、プログラムは、リムーバブル記録媒体１１１に格納（記録）しておくことができる。このようなリムーバブル記録媒体１１１は、いわゆるパッケージソフトウエアとして提供することができる。ここで、リムーバブル記録媒体１１１としては、例えば、フレキシブルディスク、CD-ROM(Compact Disc Read Only Memory)，MO(Magneto Optical)ディスク，DVD(Digital Versatile Disc)、磁気ディスク、半導体メモリ等がある。
【０５１２】
なお、プログラムは、上述したようなリムーバブル記録媒体１１１からコンピュータにインストールする他、通信網や放送網を介して、コンピュータにダウンロードし、内蔵するハードディスク１０５にインストールすることができる。すなわち、プログラムは、例えば、ダウンロードサイトから、ディジタル衛星放送用の人工衛星を介して、コンピュータに無線で転送したり、LAN(Local Area Network)、インターネットといったネットワークを介して、コンピュータに有線で転送することができる。
【０５１３】
コンピュータは、CPU(Central Processing Unit)１０２を内蔵しており、CPU１０２には、バス１０１を介して、入出力インタフェース１１０が接続されている。
【０５１４】
CPU１０２は、入出力インタフェース１１０を介して、ユーザによって、入力部１０７が操作等されることにより指令が入力されると、それに従って、ROM(Read Only Memory)１０３に格納されているプログラムを実行する。あるいは、CPU１０２は、ハードディスク１０５に格納されたプログラムを、RAM(Random Access Memory)１０４にロードして実行する。
【０５１５】
これにより、CPU１０２は、上述したフローチャートにしたがった処理、あるいは上述したブロック図の構成により行われる処理を行う。そして、CPU１０２は、その処理結果を、必要に応じて、例えば、入出力インタフェース１１０を介して、出力部１０６から出力、あるいは、通信部１０８から送信、さらには、ハードディスク１０５に記録等させる。
【０５１６】
なお、入力部１０７は、キーボードや、マウス、マイク等で構成される。また、出力部１０６は、LCD(Liquid Crystal Display)やスピーカ等で構成される。
【０５１７】
ここで、本明細書において、コンピュータがプログラムに従って行う処理は、必ずしもフローチャートとして記載された順序に沿って時系列に行われる必要はない。すなわち、コンピュータがプログラムに従って行う処理は、並列的あるいは個別に実行される処理（例えば、並列処理あるいはオブジェクトによる処理）も含む。
【０５１８】
また、プログラムは、１のコンピュータ（プロセッサ）により処理されるものであっても良いし、複数のコンピュータによって分散処理されるものであっても良い。さらに、プログラムは、遠方のコンピュータに転送されて実行されるものであっても良い。
【０５１９】
さらに、本明細書において、システムとは、複数の構成要素（装置、モジュール（部品）等）の集合を意味し、すべての構成要素が同一筐体中にあるか否かは問わない。したがって、別個の筐体に収納され、ネットワークを介して接続されている複数の装置、及び、１つの筐体の中に複数のモジュールが収納されている１つの装置は、いずれも、システムである。
【０５２０】
なお、本技術の実施の形態は、上述した実施の形態に限定されるものではなく、本技術の要旨を逸脱しない範囲において種々の変更が可能である。
【０５２１】
例えば、本技術は、１つの機能をネットワークを介して複数の装置で分担、共同して処理するクラウドコンピューティングの構成をとることができる。
【０５２２】
また、上述のフローチャートで説明した各ステップは、１つの装置で実行する他、複数の装置で分担して実行することができる。
【０５２３】
さらに、１つのステップに複数の処理が含まれる場合には、その１つのステップに含まれる複数の処理は、１つの装置で実行する他、複数の装置で分担して実行することができる。
【０５２４】
なお、本技術は、以下のような構成をとることができる。
【０５２５】
［１］
複数の電気機器で消費される電流の総和に関する総和データを用いて、FHMM(Factorial Hidden Markov Model)の各ファクタの状態にいる状態確率を求める状態推定を行う状態推定部と、
前記状態確率を用いて、各ファクタが３以上の状態を有する前記FHMMの学習を行うモデル学習部と
を備えるデータ処理装置。
［２］
前記FHMMは、
各ファクタの状態の組み合わせにおいて観測される前記総和データの観測値の平均値を求めるのに用いられる、各ファクタの各状態に固有の固有波形と、
各ファクタの状態の組み合わせにおいて観測される前記総和データの観測値の分散と、
各ファクタの状態が初期状態である初期状態確率と、
各ファクタの状態が遷移する遷移確率と
を、モデルパラメータとして有し、
前記モデル学習部は、
前記固有波形を求める波形分離学習を行う波形分離学習部と、
前記分散を求める分散学習を行う分散学習部と、
前記初期状態確率と、前記遷移確率とを求める状態変動学習を行う状態変動学習部と
を有する
［１］に記載のデータ処理装置。
［３］
前記分散学習部は、前記ファクタごとに個別の分散、又は、各ファクタの状態ごとに個別の分散を求める
［２］に記載のデータ処理装置。
［４］
前記状態推定部は、1時刻に状態が遷移するファクタの数を制限する状態遷移制限の下で、前記状態確率を求める
［２］又は［３］に記載のデータ処理装置。
［５］
前記状態推定部は、
前記平均値、及び、前記分散を用いて、各ファクタの状態の組み合わせにおいて、前記総和データが観測される観測確率を求め、
前記観測確率、及び、前記遷移確率を用いて、前記総和データの系列Y1，Y2，・・・，YTについて、前記総和データY1，Y2，・・・，Ytを観測し、時刻tに、各ファクタの状態の組み合わせzにいる前向き確率αt,zと、時刻tに、各ファクタの状態の組み合わせzにいて、その後、前記総和データYt，Yt+1，・・・，YTを観測する後ろ向き確率βt,zとを求め、
前記前向き確率αt,z、及び、前記後ろ向き確率βt,zを用いて、時刻tに、各ファクタの状態の組み合わせzにいる事後確率γt,zを求め、
前記事後確率γt,zを周辺化することにより、前記状態確率を求める
場合において、
各ファクタの状態の組み合わせを、パーティクルとして、
1時刻後のパーティクルを、前記状態遷移制限の下で予測し、前記前向き確率αt,zに基づいて、所定数のパーティクルをサンプリングすることを繰り返しながら、前記パーティクルとしての状態の組み合わせzの前記前向き確率αt,zを求め、
1時刻前のパーティクルを、前記状態遷移制限の下で予測し、前記後ろ向き確率βt,zに基づいて、所定数のパーティクルをサンプリングすることを繰り返しながら、前記パーティクルとしての状態の組み合わせzの前記後ろ向き確率βt,zを求める
［４］に記載のデータ処理装置。
［６］
前記波形分離学習部は、前記電気機器に特有の制約の下で、前記固有波形を求める
［２］ないし［５］のいずれかに記載のデータ処理装置。
［７］
前記波形分離学習部は、前記固有波形を用いて求められる前記電気機器の消費電力が負の値にならない負荷制約の下で、前記固有波形を求める
［６］に記載のデータ処理装置。
［８］
前記波形分離学習部は、前記固有波形が、前記電気機器について用意された複数の基底波形の１以上の組み合わせで表現される基底波形制約の下で、前記固有波形を求める
［６］に記載のデータ処理装置。
［９］
複数の電気機器で消費される電流の総和に関する総和データを用いて、FHMM(Factorial Hidden Markov Model)の各ファクタの状態にいる状態確率を求める状態推定を行い、
前記状態確率を用いて、各ファクタが３以上の状態を有する前記FHMMの学習を行う
ステップを含むデータ処理方法。
［１０］
複数の電気機器で消費される電流の総和に関する総和データを用いて、FHMM(Factorial Hidden Markov Model)の各ファクタの状態にいる状態確率を求める状態推定を行う状態推定部と、
前記状態確率を用いて、各ファクタが３以上の状態を有する前記FHMMの学習を行うモデル学習部と
して、コンピュータを機能させるためのプログラム。
【符号の説明】
【０５２６】
１１  データ取得部，  １２  状態推定部，  １３  モデル記憶部，  １４  モデル学習部，  １５  ラベル取得部，  １６  データ出力部，  ２１  評価部，  ２２  推定部，  ３１  波形分離学習部，  ３２  分散学習部，  ３３  状態変動学習部，  ４２  推定部，  ５１  制約付き波形分離学習部，  ５２  個別分散学習部，  ７１  評価部，  ７２  推定部，  ８１  制約付き波形分離学習部，  １０１  バス，  １０２  CPU，  １０３  ROM，  １０４  RAM，  １０５  ハードディスク，  １０６  出力部，  １０７  入力部，  １０８  通信部，  １０９  ドライブ，  １１０  入出力インタフェース，  １１１  リムーバブル記録媒体
【図１】
0006020880.tif 000034
【図２】
0006020880.tif 000035
【図３】
0006020880.tif 000036
【図４】
0006020880.tif 000037
【図５】
0006020880.tif 000038
【図６】
0006020880.tif 000039
【図７】
0006020880.tif 000040
【図８】
0006020880.tif 000041
【図９】
0006020880.tif 000042
【図１０】
0006020880.tif 000043
【図１１】
0006020880.tif 000044
【図１２】
0006020880.tif 000045
【図１３】
0006020880.tif 000046
【図１４】
0006020880.tif 000047
【図１５】
0006020880.tif 000048
【図１６】
0006020880.tif 000049
【図１７】
0006020880.tif 000050
【図１８】
0006020880.tif 000051
【図１９】
0006020880.tif 000052
【図２０】
0006020880.tif 000053
【図２１】
0006020880.tif 000054
【図２２】
0006020880.tif 000055
【図２３】
0006020880.tif 000056
【図２４】
0006020880.tif 000057
【図２５】
0006020880.tif 000058
【図２６】
0006020880.tif 000059
【図２７】
0006020880.tif 000060
【図２８】
0006020880.tif 000061
【図２９】
0006020880.tif 000062
【図３０】
0006020880.tif 000063