JP2011118755A - 学習装置、学習方法、およびプログラム - Google Patents
学習装置、学習方法、およびプログラム Download PDFInfo
- Publication number
- JP2011118755A JP2011118755A JP2009276691A JP2009276691A JP2011118755A JP 2011118755 A JP2011118755 A JP 2011118755A JP 2009276691 A JP2009276691 A JP 2009276691A JP 2009276691 A JP2009276691 A JP 2009276691A JP 2011118755 A JP2011118755 A JP 2011118755A
- Authority
- JP
- Japan
- Prior art keywords
- data
- likelihood
- state
- time
- hidden markov
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
- G08G1/0133—Traffic data processing for classifying traffic situation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/29—Graphical models, e.g. Bayesian networks
- G06F18/295—Markov models or related models, e.g. semi-Markov models; Markov random fields; Networks embedding Markov models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0108—Measuring and analyzing of parameters relative to traffic conditions based on the source of data
- G08G1/0112—Measuring and analyzing of parameters relative to traffic conditions based on the source of data from the vehicle, e.g. floating car data [FCD]
Abstract
【課題】欠落した時系列データから、より正確に隠れマルコフモデルを学習する。
【解決手段】データ補間部31は、時系列データにおいて欠落しているデータを補間する。状態尤度計算部32乃至バックワード尤度計算部35は、正常位置データに対する状態尤度と補間位置データに対する状態尤度とを異なる条件で計算して、データが補間された時系列データに対する隠れマルコフモデルの尤度を計算する。初期確率遷移確率推定部36乃至観測確率推定部38は、時系列データから隠れマルコフモデルを推定するとともに、状態尤度計算部32乃至バックワード尤度計算部35により計算される尤度が大きくなるように、隠れマルコフモデルを更新する。本発明は、例えば、ユーザの活動モデルを学習する学習装置に適用できる。
【選択図】図4
【解決手段】データ補間部31は、時系列データにおいて欠落しているデータを補間する。状態尤度計算部32乃至バックワード尤度計算部35は、正常位置データに対する状態尤度と補間位置データに対する状態尤度とを異なる条件で計算して、データが補間された時系列データに対する隠れマルコフモデルの尤度を計算する。初期確率遷移確率推定部36乃至観測確率推定部38は、時系列データから隠れマルコフモデルを推定するとともに、状態尤度計算部32乃至バックワード尤度計算部35により計算される尤度が大きくなるように、隠れマルコフモデルを更新する。本発明は、例えば、ユーザの活動モデルを学習する学習装置に適用できる。
【選択図】図4
Description
本発明は、学習装置、学習方法、およびプログラムに関し、特に、時系列データからユーザの活動状態を隠れマルコフモデルとして学習する学習装置、学習方法、およびプログラムに関する。
近年、ユーザが身に着けられるセンサであるウェアラブルセンサから得られる時系列データを用いてユーザの状態をモデル化して学習し、学習により得られたモデルを用いてユーザの現在の状態を認識する研究が盛んである(例えば、特許文献1,2、非特許文献1)。
また、本出願人は、未来の所望の時刻におけるユーザの活動状態の複数の可能性を確率的に予測する方法を、特願2009−180780(以下、先願1という)として先に提案している。先願1の方法では、時系列データからユーザの活動状態を確率的状態遷移モデルとして学習し、学習した確率的状態遷移モデルを用いて現在の活動状態を認識し、「所定時間後」のユーザの活動状態を確率的に予測することができる。そして、「所定時間後」のユーザの活動状態の予測の一例として、ユーザの現在の位置を認識し、所定時間後のユーザの行き先(場所)を予測する例が示されている。
さらに、本出願人は、複数の目的地への到達確率、経路、到達時間を予測する方法、および、目的地の候補を自動的に検出する方法を、特願2009−208064(以下、先願2という)として先に提案している。また、先願2では、確率的状態遷移モデルを構成する状態ノードの中から、目的地に対応する状態ノードを見つける方法が提案されている。
‘‘Life Patterns: structure from wearable sensors’’,Brian Patrick Clarkson, Doctor Thesis, MIT, 2002
先願1および2では、GPSセンサにより測定された位置(緯度経度)の時系列データから、ユーザの活動状態を隠れマルコフモデルとして学習する例が示されている。
この場合の隠れマルコフモデルは、等間隔にセンシングされた時系列データから、その背後の「隠れた」変数のダイナミクスを状態の確率的遷移として表現する確率モデルである。このように、時系列データから隠れマルコフモデルを学習する場合、使用する時系列データは、等間隔でサンプリングされたデータであることが望ましい。
一方、GPSセンサは、人工衛星からの信号を捕捉して、現在地の緯度経度を測定するセンサである。従って、GPSセンサは、人工衛星からの信号の届かない地下や建物内などでは、現在地の緯度経度を測定することができない。そして、GPSセンサが緯度経度を測定できなくなってから測定を再開するまでの間、位置(緯度経度)の時系列データに欠落期間が生じる。
このような欠落期間が生じた時系列データを用いても、隠れマルコフモデルの学習は可能である。しかし、この場合に得られるモデルには、欠落期間の前に対応する状態ノードと欠落期間の後に対応する状態ノードとを短絡する状態遷移が現れる。
例えば、図1は、先願2に示される学習方法により、実験者が鉄道を利用して自宅から会社に出勤する際にGPSセンサにより位置を10日間にわたって測定した時系列データを用いて、通勤経路を隠れマルコフモデルとして学習した結果を示している。なお、図中、小さい丸は、学習した隠れマルコフモデルの各状態ノードの位置を示している。また、各状態ノード間を結ぶ直線は、状態ノード間の状態遷移を示している。また、図1の横軸は経度を表し、縦軸は緯度を表している。
また、図2は、通勤に利用した路線の各駅の位置関係が分かるように、図1の上に各駅の位置と駅名を示した図である。
この路線は、長原駅付近から荏原中延駅付近まで地下を通るため、その付近で時系列データの欠落が生じる。そのため、例えば、洗足池駅付近と旗の台駅付近、旗の台駅付近と戸越銀座駅付近、洗足池駅付近と荏原中延駅付近、洗足池駅付近と大崎広小路駅付近をそれぞれ短絡したようなデータが時系列データに出現する。従って、この時系列データを用いて学習処理を行うと、洗足池駅付近と旗の台駅付近、旗の台駅付近と戸越銀座駅付近、洗足池駅付近と荏原中延駅付近、洗足池駅付近と大崎広小路駅付近の状態ノードをそれぞれ直接結ぶ状態遷移(すなわち経路)を含んだモデルが学習されてしまう。
例えば、このように学習されたモデルを用いて、目的地への到達時間を予測した場合、目的地までの経路にワープしたような区間が出現する。そのため、目的地に到達するまでに遷移する状態ノードの数(以下、ステップ数という)が、本来あるべき値より少なくなり、目的地への到達時間の予測精度が悪化する。
また、同じ経路に対して、正常な経路と短絡された経路が存在することになり、それぞれ別の分岐した経路として認識されてしまう。さらに、時系列データが欠落する場所が一定でない場合、欠落する場所に応じて異なる短絡経路が学習され、その結果、学習したモデルに、大量の分岐経路が現れる。その結果、目的地までの経路を探索する場合に、本来存在しない多数の分岐経路を探索する必要が生じ、処理にかかる負荷が増大し、処理時間が長くなったり、処理結果の信頼度が低下したりする。
本発明は、このような状況を鑑みてなされたものであり、欠落した時系列データから、より正確に隠れマルコフモデルを学習できるようにするものである。
本発明の一側面の学習装置は、
時系列データにおいて欠落しているデータを補間する補間手段と、
前記時系列データから隠れマルコフモデルを推定する推定手段と、
推定された前記隠れマルコフモデルの尤度を計算する尤度計算手段と
を備え、
前記尤度計算手段は、欠落していないデータである正常データに対する尤度と補間されたデータである補間データに対する尤度とを異なる条件で計算して、データが補間された前記時系列データに対する前記隠れマルコフモデルの尤度を計算し、
前記推定手段は、前記尤度計算手段により計算される尤度が大きくなるように、前記隠れマルコフモデルを更新する。
時系列データにおいて欠落しているデータを補間する補間手段と、
前記時系列データから隠れマルコフモデルを推定する推定手段と、
推定された前記隠れマルコフモデルの尤度を計算する尤度計算手段と
を備え、
前記尤度計算手段は、欠落していないデータである正常データに対する尤度と補間されたデータである補間データに対する尤度とを異なる条件で計算して、データが補間された前記時系列データに対する前記隠れマルコフモデルの尤度を計算し、
前記推定手段は、前記尤度計算手段により計算される尤度が大きくなるように、前記隠れマルコフモデルを更新する。
本発明の一側面の学習方法は、
時系列データから隠れマルコフモデルを学習する学習装置が、
前記時系列データにおいて欠落しているデータを補間し、
前記時系列データから前記隠れマルコフモデルを推定し、
推定された前記隠れマルコフモデルの尤度を計算するとともに、
欠落していないデータである正常データに対する尤度と補間されたデータである補間データに対する尤度とを異なる条件で計算して、データが補間された前記時系列データに対する前記隠れマルコフモデルの尤度を計算し、
計算される尤度が大きくなるように、前記隠れマルコフモデルを更新する
ステップを含む。
時系列データから隠れマルコフモデルを学習する学習装置が、
前記時系列データにおいて欠落しているデータを補間し、
前記時系列データから前記隠れマルコフモデルを推定し、
推定された前記隠れマルコフモデルの尤度を計算するとともに、
欠落していないデータである正常データに対する尤度と補間されたデータである補間データに対する尤度とを異なる条件で計算して、データが補間された前記時系列データに対する前記隠れマルコフモデルの尤度を計算し、
計算される尤度が大きくなるように、前記隠れマルコフモデルを更新する
ステップを含む。
本発明の一側面のプログラムは、
時系列データにおいて欠落しているデータを補間し、
前記時系列データから隠れマルコフモデルを推定し、
推定された前記隠れマルコフモデルの尤度を計算するとともに、
欠落していないデータである正常データに対する尤度と補間されたデータである補間データに対する尤度とを異なる条件で計算して、データが補間された前記時系列データに対する前記隠れマルコフモデルの尤度を計算し、
計算される尤度が大きくなるように、前記隠れマルコフモデルを更新する
ステップを含む処理をコンピュータに実行させる。
時系列データにおいて欠落しているデータを補間し、
前記時系列データから隠れマルコフモデルを推定し、
推定された前記隠れマルコフモデルの尤度を計算するとともに、
欠落していないデータである正常データに対する尤度と補間されたデータである補間データに対する尤度とを異なる条件で計算して、データが補間された前記時系列データに対する前記隠れマルコフモデルの尤度を計算し、
計算される尤度が大きくなるように、前記隠れマルコフモデルを更新する
ステップを含む処理をコンピュータに実行させる。
本発明の一側面においては、時系列データにおいて欠落しているデータが補間され、前記時系列データから隠れマルコフモデルが推定され、推定された前記隠れマルコフモデルの尤度が計算されるとともに、欠落していないデータである正常データに対する尤度と補間されたデータである補間データに対する尤度とが異なる条件で計算され、データが補間された前記時系列データに対する前記隠れマルコフモデルの尤度が計算され、計算される尤度が大きくなるように、前記隠れマルコフモデルが更新される。
本発明の一側面によれば、欠落した時系列データから、より正確に隠れマルコフモデルを学習することができる。
以下、本発明を実施するための形態(以下、実施の形態という)について説明する。なお、説明は以下の順序で行う。
1.実施の形態
2.変形例
1.実施の形態
2.変形例
<1.実施の形態>
[予測システムの実施の形態のブロック図]
図3は、本発明を適用した予測システムの一実施の形態の構成例を示すブロック図である。
[予測システムの実施の形態のブロック図]
図3は、本発明を適用した予測システムの一実施の形態の構成例を示すブロック図である。
予測システム1は、GPSセンサ11、時系列データ記憶部12、行動学習部13、行動認識部14、行動予測部15、到達時間予測部16、操作部17、および表示部18により構成される。
予測システム1は、GPSセンサ11により取得される現在地および現在時刻を示す時系列データから、ユーザの活動状態(行動・活動パターンを表した状態)を確率的状態遷移モデルとして学習する学習処理を行う。また、予測システム1は、学習処理により得られたパラメータで表される確率的状態遷移モデル(ユーザ活動モデル)を用いて、ユーザが指定した目的地までの経路および時間を予測する予測処理も行う。
図3において、点線の矢印は、学習処理におけるデータの流れを示しており、実線の矢印は、予測処理におけるデータの流れを示している。
GPSセンサ11は、一定時間間隔(たとえば、15秒間隔)で人工衛星からの信号を受信し、自身の現在地の緯度経度を測定する。そして、GPSセンサ11は、学習処理においては、測定した位置(緯度経度)と測定した時刻(以下、測定時刻ともいう)のデータ(以下、GPSデータという)を時系列データ記憶部12に供給する。また、GPSセンサ11は、予測処理においては、GPSデータを行動認識部14に供給する。
時系列データ記憶部12は、GPSセンサ11により連続して取得されたGPSデータ、すなわち、位置と測定時刻の時系列データを記憶する。ユーザの行動・活動パターンを学習するので、例えば、数日分程度など、ある程度の期間について蓄積された時系列データが必要である。
行動学習部13は、時系列データ記憶部12に記憶されている時系列データに基づいて、GPSセンサ11が組み込まれた機器を携行するユーザの活動状態を、確率的状態遷移モデルとして学習する。時系列データはユーザの位置と測定時刻を示すデータであるので、確率的状態遷移モデルとして学習されるユーザの活動状態は、ユーザの現在地の時系列変化、すなわちユーザの移動経路を表す状態となる。学習に使用される確率的状態遷移モデルとしては、例えば、エルゴディックHMM(Hidden Markov Model,隠れマルコフモデル)などの、隠れ状態を含む確率的状態遷移モデルを採用することができる。本実施の形態では、確率的状態遷移モデルとして、エルゴディックHMMにスパース制約を与えたものを採用する。なお、スパース制約を与えたエルゴディックHMM、エルゴディックHMMのパラメータの算出方法等については、図7乃至図9を参照して後述する。
また、上述したように、GPSセンサ11が現在地を測定することができない場合、時系列データに欠落期間が生じる。このような場合、行動学習部13は、図13を参照して後述するように、時系列データにおいて欠落したデータを補間して、ユーザの活動状態を学習する。
行動学習部13は、学習結果を示すデータを表示部18に供給し、表示させる。また、行動学習部13は、学習処理により得られた確率的状態遷移モデルのパラメータを行動認識部14および行動予測部15に供給する。
行動認識部14は、学習により得られたパラメータの確率的状態遷移モデルを用いて、GPSセンサ11からリアルタイムに供給される位置の時系列データから、ユーザの現在の活動状態、すなわち、ユーザの現在地を認識する。行動認識部14は、ユーザの現在の状態ノードのノード番号を行動予測部15に供給する。
行動予測部15は、学習により得られたパラメータの確率的状態遷移モデルを用いて、行動認識部14から供給される状態ノードのノード番号が示すユーザの現在地から、ユーザが取りうる経路を過不足なく探索(予測)する。また、行動予測部15は、探索された経路ごとの生起確率を計算することにより、探索された経路が選択される確率である選択確率を予測する。
到達時間予測部16には、行動予測部15から、ユーザが取りうる経路と、その選択確率が供給される。また、到達時間予測部16には、操作部17から、ユーザが指定した目的地を示す情報が供給される。
到達時間予測部16は、探索結果としてのユーザが取りうる経路から、ユーザが指定した目的地を含む経路を抽出し、抽出された各経路について目的地までの到達時間を予測する。また、到達時間予測部16は、目的地へ到達する確率である到達確率を予測する。到達時間予測部16は、目的地について複数の経路が存在する場合、その複数の経路の選択確率の和を、目的地の到達確率として算出する。目的地への経路が1つのみである場合には、その経路の選択確率が、そのまま目的地の到達確率となる。そして、到達時間予測部16は、予測結果を示す情報を表示部18に供給し、表示させる。
操作部17は、ユーザが入力した目的地についての情報を受け付け、到達時間予測部16に供給する。表示部18は、行動学習部13または到達時間予測部16から供給される情報を表示する。
[行動学習部の機能の構成の例]
図4は、予測システム1の行動学習部13の機能の構成の例を示すブロック図である。
図4は、予測システム1の行動学習部13の機能の構成の例を示すブロック図である。
行動学習部13は、データ補間部31、状態尤度計算部32、状態尤度修正部33、フォワード尤度計算部34、バックワード尤度計算部35、初期確率遷移確率推定部36、観測確率推定前処理部37、および、観測確率推定部38により構成される。
データ補間部31は、時系列データ記憶部12に記憶されている時系列データを取得し、図13を参照して後述するように、時系列データの欠落期間のデータを補間する。データ補間部31は、補間した時系列データを状態尤度計算部32に供給する。
状態尤度計算部32は、図12を参照して後述するように、ユーザの活動モデルを表すHMMの各状態の時系列データに対する尤度(以下、状態尤度ともいう)を計算する。状態尤度計算部32は、計算した状態尤度を示すデータ、および、時系列データを状態尤度修正部33に供給する。
状態尤度修正部33は、図12を参照して後述するように、時系列データの欠落期間に対する状態尤度を修正し、修正した状態尤度を示すデータ、および、時系列データをフォワード尤度計算部34に供給する。
フォワード尤度計算部34は、図12を参照して後述するように、ユーザの活動モデルを表すHMMの時系列データに対するフォワード尤度を計算する。フォワード尤度計算部34は、状態尤度修正部33から取得したデータに、計算したフォワード尤度を示すデータを追加して、バックワード尤度計算部35に供給する。
バックワード尤度計算部35は、図12を参照して後述するように、ユーザの活動モデルを表すHMMの時系列データに対するバックワード尤度を計算する。バックワード尤度計算部35は、フォワード尤度計算部34から取得したデータに、計算したバックワード尤度を示すデータを追加して、初期確率遷移確率推定部36に供給する。
初期確率遷移確率推定部36は、図12を参照して後述するように、ユーザの活動モデルを表すHMMの初期確率および状態遷移確率を推定する。初期確率遷移確率推定部36は、バックワード尤度計算部35から取得したデータに、計算した初期確率および状態遷移確率を示すデータを追加して、観測確率推定前処理部37に供給する。
観測確率推定前処理部37は、図12を参照して後述するように、ユーザの活動モデルを表すHMMの観測確率を推定する前に、時系列データの欠落期間に対するフォワード尤度およびバックワード尤度を変更する。観測確率推定前処理部37は、初期確率遷移確率推定部36から取得したデータに対して、フォワード尤度およびバックワード尤度を変更したデータを観測確率推定部38に供給する。
観測確率推定部38は、図12を参照して後述するように、ユーザの活動モデルを表すHMMの初期確率を推定する。観測確率推定部38は、HMMのパラメータ(初期確率、状態遷移確率、観測確率)を示すデータを、行動認識部14および行動予測部15、または、状態尤度計算部32に供給する。
[予測システムのハードウエア構成例]
以上のように構成される予測システム1は、例えば、図5に示されるハードウエア構成を採用することができる。即ち、図5は、予測システム1のハードウエア構成例を示すブロック図である。
以上のように構成される予測システム1は、例えば、図5に示されるハードウエア構成を採用することができる。即ち、図5は、予測システム1のハードウエア構成例を示すブロック図である。
図5において、予測システム1は、3台のモバイル端末51−1乃至51−3とサーバ52とにより構成されている。モバイル端末51−1乃至51−3は、同一機能を有する同型のモバイル端末51であるが、モバイル端末51−1乃至51−3では、それを所有するユーザが異なる。従って、図5では、3台のモバイル端末51−1乃至51−3のみが示されているが、実際には、ユーザ数に応じた数のモバイル端末51が存在する。
モバイル端末51は、無線通信及びインターネット等のネットワークを介した通信により、サーバ52とデータの授受を行うことができる。サーバ52は、モバイル端末51から送信されてくるデータを受信し、受信したデータに対し所定の処理を行う。そして、サーバ52は、データ処理の処理結果を無線通信等によりモバイル端末51に送信する。
従って、モバイル端末51とサーバ52は、無線または有線による通信を行う通信部を少なくとも有する。
さらに、モバイル端末51が、図3のGPSセンサ11、操作部17、および表示部18を備え、サーバ52が、図3の時系列データ記憶部12、行動学習部13、行動認識部14、行動予測部15、および到達時間予測部16を備える構成を採用することができる。
この構成が採用される場合、学習処理において、モバイル端末51が、GPSセンサ11により取得された時系列データを送信する。サーバ52は、受信した学習用の時系列データに基づき、ユーザの活動状態を確率的状態遷移モデルにより学習する。そして、予測処理において、モバイル端末51が、操作部17を介してユーザによって指定された目的地を送信するとともに、GPSセンサ11によりリアルタイムに取得されるGPSデータを送信する。サーバ52は、学習により得られたパラメータを用いて、ユーザの現在の活動状態、すなわち、ユーザの現在地を認識し、さらに、指定された目的地までの経路および時間を処理結果としてモバイル端末51に送信する。モバイル端末51は、サーバ52から送信されてきた処理結果を表示部18に表示する。
また例えば、モバイル端末51が、図3のGPSセンサ11、行動認識部14、行動予測部15、到達時間予測部16、操作部17、および表示部18を備え、サーバ52が、図3の時系列データ記憶部12および行動学習部13を備える構成を採用することができる。
この構成が採用される場合、学習処理において、モバイル端末51が、GPSセンサ11により取得された時系列データを送信する。サーバ52は、受信した学習用の時系列データに基づき、ユーザの活動状態を確率的状態遷移モデルにより学習し、学習により得られたパラメータをモバイル端末51に送信する。そして、予測処理において、モバイル端末51が、GPSセンサ11によりリアルタイムに取得されるGPSデータを、サーバ52から受信したパラメータを用いて、ユーザの現在地を認識し、さらに、指定された目的地までの経路および時間を演算する。そして、モバイル端末51は、演算結果としての目的地までの経路および時間を表示部18に表示する。
以上のようなモバイル端末51とサーバ52との間の役割分担は、それぞれのデータ処理装置としての処理能力や通信環境に応じて決定することができる。
学習処理は、処理に要する1回あたりの時間は非常に長いが、それほど頻繁に処理する必要はない。従って、一般的には、携行可能なモバイル端末51よりもサーバ52の方が処理能力が高いので、サーバ52に、一日に一回程度蓄積された時系列データに基づいて学習処理(パラメータの更新)を行わせるようにすることができる。
一方、予測処理は、時々刻々とリアルタイムに更新される位置データに対応させて迅速に処理し、表示することが望ましいので、モバイル端末51で処理を行う方が望ましい。通信環境がリッチであれば、上述したようにサーバ52に予測処理も行わせ、予測結果のみをサーバ52から受信する方が、携行可能な小型化が要求されるモバイル端末51の負荷が軽減され、望ましい。
また、モバイル端末51単独で、データ処理装置として学習処理および予測処理を高速に行うことが可能である場合には、図3の予測システム1の構成すべてをモバイル端末51が備えるようにすることも勿論可能である。
[入力される時系列データの例]
図6は、予測システム1で取得された位置の時系列データの例を示している。図6において、横軸は経度を表し、縦軸は緯度を表している。なお、図6では、図を分かりやすくするために、各位置における時刻を示していないが、実際には、各位置における時刻も時系列データに含まれる。
図6は、予測システム1で取得された位置の時系列データの例を示している。図6において、横軸は経度を表し、縦軸は緯度を表している。なお、図6では、図を分かりやすくするために、各位置における時刻を示していないが、実際には、各位置における時刻も時系列データに含まれる。
図6に示される時系列データは、実験者の1ヶ月半程度の期間に蓄積された時系列データを示している。図6に示されるように、時系列データは、主に、自宅周辺と、勤務先などの4か所の外出先を移動したデータとなっている。なお、この時系列データには、人工衛星を捕捉できず、位置データが欠落しているデータも含まれている。
[エルゴディックHMMについて]
次に、予測システム1が、学習モデルとして採用するエルゴディックHMMについて説明する。
次に、予測システム1が、学習モデルとして採用するエルゴディックHMMについて説明する。
図7は、HMMの例を示している。
HMMは、状態と状態間遷移とを有する状態遷移モデルである。
図7は、3状態のHMMの例を示している。
図7において(以降の図においても同様)、丸印は、状態を表し、矢印は、状態遷移を表す。なお、状態は、上述のユーザの活動状態に対応し、状態ノードと同義である。
また、図7において、si(図7では、i=1,2,3)は、状態(ノード)を表し、aijは、状態siから状態sjへの状態遷移確率を表す。さらに、bj(x)は、状態sjへの状態遷移時に、観測値xが観測される出力確率密度関数(観測確率)を表し、πiは、状態siが初期状態である初期確率を表す。
なお、出力確率密度関数bj(x)としては、例えば、混合正規確率分布等が用いられる。
ここで、HMM(連続HMM)は、状態遷移確率aij、出力確率密度関数bj(x)、及び初期確率πiによって定義される。これらの状態遷移確率aij、出力確率密度関数bj(x)、及び初期確率πiを、HMMのパラメータλ={aij,bj(x), πi,i=1,2,・・・,M,j=1,2,・・・,M}という。Mは、HMMの状態数を表す。
HMMのパラメータλを推定する方法としては、Baum-Welchの最尤推定法が広く利用されている。Baum-Welchの最尤推定法は、EMアルゴリズム(EM(Expectation-Maximization) algorithm)に基づくパラメータの推定方法である。
Baum-Welchの最尤推定法によれば、観測される時系列データx=x1,x2,・・・,xTに基づき、その時系列データが観測(生起)される確率である生起確率から求まる尤度を最大化するように、HMMのパラメータλの推定が行われる。ここで、xtは、時刻tに観測される信号(サンプル値)を表し、Tは、時系列データの長さ(サンプル数)を表す。
なお、Baum-Welchの最尤推定法を用いたパラメータλの推定方法の詳細については、図12を参照して後述する。
また、Baum-Welchの最尤推定法については、例えば、‘‘パターン認識と機械学習(下)’’,C.M.ビショップ著,P. 333(英語原書:‘‘Pattern Recognition and Machine Learning (Information Science and Statistics) ’’,Christopher M. BishopSpringer, New York, 2006.)(以下、文献Aと称する)に記載されている。
なお、Baum-Welchの最尤推定法は、尤度最大化に基づくパラメータ推定方法ではあるが、最適性を保証するものではなく、HMMの構造やパラメータλの初期値によっては、局所解に収束することがある。
HMMは、音声認識で広く利用されているが、音声認識で利用されるHMMでは、一般に、状態の数や状態遷移の仕方等はあらかじめ決定される。
図8は、音声認識で利用されるHMMの例を示している。
図8のHMMは、left-to-right型と呼ばれる。
図8では、状態数は3になっており、状態遷移は、自己遷移(状態siから状態siへの状態遷移)と、左から右隣の状態への状態遷移とのみを許す構造に制約されている。
図8のHMMのように、状態遷移に制約があるHMMに対して、図7に示した、状態遷移に制約がないHMM、すなわち、任意の状態siから任意の状態sjへの状態遷移が可能なHMMは、エルゴディック(Ergodic)HMMと呼ばれる。
エルゴディックHMMは、構造としては最も自由度の高いHMMであるが、状態数が多くなると、パラメータλの推定が困難となる。
例えば、エルゴディックHMMの状態数が、1000である場合、状態遷移の数は、100万(=1000×1000)となる。
したがって、この場合、パラメータλのうちの、例えば、状態遷移確率aijについては、100万個の状態遷移確率aijを推定することが必要となる。
そこで、状態に対して設定する状態遷移には、例えば、スパース(Sparse)な構造であるという制約(スパース制約)をかけることができる。
ここで、スパースな構造とは、任意の状態から任意の状態への状態遷移が可能なエルゴディックHMMのような密な状態遷移ではなく、ある状態から状態遷移することができる状態が非常に限定されている構造である。なお、ここでは、スパースな構造であっても、他の状態への状態遷移は、少なくとも1つ存在し、また、自己遷移は存在することとする。
図9は、スパース制約を与えたHMMを示している。
ここで、図9では、2つの状態を結ぶ双方向の矢印は、その2つの状態の一方から他方への状態遷移と、他方から一方への状態遷移とを表す。また、図9において、各状態は、自己遷移が可能であり、その自己遷移を表す矢印の図示は、省略されている。
図9では、16個の状態が、2次元空間上に格子状に配置されている。すなわち、図9では、横方向に、4個の状態が配置され、縦方向にも、4個の状態が配置されている。
いま、横方向に隣接する状態どうしの距離、及び、縦方向に隣接する状態どうしの距離を、いずれも1とすると、図9Aは、距離が1以下の状態への状態遷移は可能とし、他の状態への状態遷移はできないというスパース制約を与えたHMMを示している。
また、図9Bは、距離が√2以下の状態への状態遷移は可能とし、他の状態への状態遷移はできないというスパース制約を与えたHMMを示している。
本実施の形態では、GPSセンサ11が取得したGPSデータが、時系列データとして、時系列データ記憶部12に供給される。行動学習部13は、時系列データ記憶部12に記憶されているGPSデータの時系列データのうちの位置の時系列データx=x1,x2,・・・,xTを用い、ユーザ活動モデルを表すHMMのパラメータλを推定する。
即ち、ユーザの移動軌跡を表す各時刻の位置(緯度経度)のデータが、HMMの状態sjのいずれかに対応する地図上の一点から、所定の分散値の広がりを持って正規分布した確率変数の観測データであると考える。行動学習部13は、各状態sjに対応する地図上の一点とその分散値、および状態遷移確率aijを最適化する。
行動認識部14は、学習により得られたユーザ活動モデル(HMM)に対して、ビタビ法を適用し、GPSセンサ11からの位置データx=x1,x2,・・・,xTが観測される尤度を最大にする状態遷移の過程(状態の系列)(パス)(以下、最尤パスともいう)を求める。これにより、ユーザの現在の活動状態、即ち、ユーザの現在地に対応する状態siが認識される。
ここで、ビタビ法とは、各状態siを始点とする状態遷移のパスの中で、時刻tに、状態siから状態sjに状態遷移する状態遷移確率aijと、その状態遷移において、位置データx=x1,x2,・・・,xTのうちの時刻tのサンプル値xtが観測される確率(出力確率密度関数bj(x)から求められる出力確率)とを、処理後時系列データxの長さTに亘って累積した値(生起確率)を最大にするパス(最尤パス)を決定するアルゴリズムである。ビタビ法の詳細については上述の文献AのP.347に記載されている。
[行動予測部15による経路の探索処理]
次に、行動予測部15による経路の探索処理について説明する。
次に、行動予測部15による経路の探索処理について説明する。
学習により得られたHMMの各状態siは、地図上の所定の点(位置)を表し、状態siと状態sjが結ばれているとき、状態siから状態sjを移動する経路を表していると考えることができる。
この場合、状態siに対応する各点は、端点、通過点、分岐点、ループのいずれかに分類することができる。端点とは、自己遷移以外の確率が極めて小さく(自己遷移以外の確率が所定の値以下であり)、次に移動可能な点がない点である。通過点とは、自己遷移以外に有意な遷移が一つある、換言すれば、次に移動可能な点が一つある点である。分岐点とは、自己遷移以外に有意な遷移が二つ以上ある、換言すれば、次に移動可能な点が二つ以上ある点である。ループとは、これまで通過した経路上のどれかと一致する点である。
目的地への経路を探索する場合、異なる経路がある場合には、それぞれの経路について必要時間等の情報を提示することが望まれる。そこで、可能な経路を過不足なく探索するために、次の条件を設定する。
(1)一度分岐した経路は再度合流した場合でも、別の経路とみなす。
(2)経路内に端点か、これまで通過した経路内に含まれる点が現れた場合、その経路の探索を終了する。
(1)一度分岐した経路は再度合流した場合でも、別の経路とみなす。
(2)経路内に端点か、これまで通過した経路内に含まれる点が現れた場合、その経路の探索を終了する。
行動予測部15は、行動認識部14により認識されたユーザの現在の活動状態、即ち、ユーザの現在の点を出発点として、次の移動先としての状態遷移が可能な点を端点、通過点、分岐点、ループのいずれかに分類することを(2)の終了条件まで繰り返す。
行動予測部15は、現在の点が端点であると分類された場合、現在の点をここまでの経路に接続してから、この経路の探索を終了する。
一方、現在の点が通過点であると分類された場合、行動予測部15は、現在の点をここまでの経路に接続してから、次の点へ移動する。
また、現在の点が分岐点であると分類された場合、行動予測部15は、現在の点をこれまでの経路に接続して、さらに分岐の数だけ、これまでの経路を複製し、分岐点と接続する。そして、行動予測部15は、分岐点の1つを次の点として移動する。
現在の点がループであると分類された場合、行動予測部15は、現在の点をこれまでの経路に接続せずに、この経路の探索を終了する。なお、現在の点から、1つ前の点に経路を逆戻りする場合はループに含まれるため考慮しない。
[探索処理の例]
図10は、行動予測部15による経路の探索処理の簡単な例を示している。
図10は、行動予測部15による経路の探索処理の簡単な例を示している。
図10の例において、状態s1が現在地である場合、最終的に3通りの経路が探索されることになる。1つめの経路は、状態s1から状態s5,状態s6等を経由して状態s10までの経路(以下、経路Aともいう。)である。2つめの経路は、状態s1から状態s5,状態s11,状態s14,状態s23等を経由して状態s29までの経路(以下、経路Bともいう。)である。3つめの経路は、状態s1から状態s5,状態s11,状態s19,状態s23等を経由して状態s29までの経路(以下、経路Cともいう。)である。
行動予測部15は、探索された各経路が選択される確率(経路の選択確率)を計算する。経路の選択確率は、経路を構成する状態間の遷移確率を順次乗算することで求められる。ただし、次の状態に遷移する場合のみを考慮し、その場所に滞留する場合は考慮する必要がないので、学習により求められた各状態の状態遷移確率aijから、自己遷移確率を除いて規格化された遷移確率[aij]を用いて、経路の選択確率が求められる。
したがって、例えば、図10の状態s5の状態遷移確率aijが、自己遷移確率a5,5=0.5,遷移確率a5,6=0.2,遷移確率a5,11=0.3である場合、状態s5から状態s6または状態s11に分岐する場合の遷移確率[a5,6]および遷移確率[a5,11]は、それぞれ、0.4,0.6となる。
探索された経路の状態siのノード番号iが、(y1,y2,・・・,yn)であるとき、この経路の選択確率は、規格化された遷移確率[aij]を用いて、次式(2)で表すことができる。
図10の例では、経路Aの選択確率は、0.4である。また、経路Bの選択確率は、0.24=0.6×0.4である。経路Cの選択確率は、0.36=0.6×0.6である。そして、計算された経路の選択確率の総和は1=0.4+0.24+0.36であり、過不足ない探索を実現することができることがわかる。
以上のように、現在地に基づいて探索された各経路とその選択確率が、行動予測部15から到達時間予測部16に供給される。
到達時間予測部16は、行動予測部15によって探索された経路から、ユーザが指定した目的地を含む経路を抽出し、抽出された各経路について目的地までの時間を予測する。
例えば、図10の例では、探索された3つの経路A乃至Cのうち、目的地である状態s28を含む経路は経路Bと経路Cである。到達時間予測部16は、経路Bまたは経路Cを通って、目的地である状態s28に到達するまでの時間を予測する。
なお、目的地を含む経路が多数あり、全ての経路を表示すると見づらくなる場合や経路の提示数が所定数に設定されている場合には、目的地を含む全ての経路のなかから、表示部18に表示させる経路(以下、適宜、表示経路ともいう。)を決定する必要がある。そのような場合、行動予測部15では、各経路について選択確率が算出されているので、到達時間予測部16は、選択確率の高い順に、所定数の経路を表示経路として決定することができる。
現在時刻t1の現在地が状態sy1であり、時刻(t1,t2,・・・,tg)における決定された経路が(sy1,sy2,・・・,syg)であるとする。換言すれば、決定された経路の状態siのノード番号iが(y1,y2,・・・,yg)であるとする。以下、簡単のため、位置に相当する状態siを、単に、そのノード番号iで表わす場合もある。
現在時刻t1での現在地y1は、行動認識部14の認識により確定しているので、現在時刻t1の現在地がy1である確率Py1(t1)は、
Py1(t1)=1
である。また、現在時刻t1にy1以外の他の状態にいる確率は0である。
Py1(t1)=1
である。また、現在時刻t1にy1以外の他の状態にいる確率は0である。
一方、所定の時刻tnにノード番号ynにいる確率Pyn(tn)は、
で表すことができる。式(3)の右辺第一項は、もともとその位置ynにいて、自己遷移した場合の確率を表し、右辺第二項は、1つ前の位置yn-1から位置ynに遷移してきた場合の確率を表している。式(3)では、経路の選択確率の計算とは異なり、学習により得られた状態遷移確率aijがそのまま利用される。
即ち、予測値<tg>は、現在時刻から、「その直前の時刻tg-1に状態sygの1つ前の状態syg-1にいて、時刻tgに状態sygに移動するとき」までの時間の期待値で表される。
[背景技術]の欄で提示した先願1の方法により、目的地への到達時間の予測値を求める場合、「所定時間後」の目的地に相当する状態の状態遷移確率aijを時間tに関して積分する必要があるが、この場合、積分区間をどれくらいに設定するのかが問題となる。先願1の方法では、ループする経路を通って目的地に到達する場合を区別できない。そのため、目的地までの途中の経路にループが存在している場合、積分区間を長めに設定すると、ループを介して二度目、三度目に目的地に到達する場合も含まれるため、目的地までの到達時間を正しく算出することができない。
式(4)で表される、本発明による目的地への到達時間の予測値の演算では、同じく、時刻tについて積分(Σ)する必要はある。しかし、探索経路において、ループする経路を通って目的地に到達する場合が除外されているため、期待値の積分区間として十分長い区間を設定することが可能である。式(4)における積分区間は、例えば、学習した経路のなかで移動時間の最大の移動時間の1倍や2倍程度とすることができる。
[ユーザ活動モデル学習処理]
次に、図11のフローチャートを参照して、ユーザの移動経路を、ユーザの活動状態を表す確率的状態遷移モデルとして学習するユーザ活動モデル学習処理について説明する。
次に、図11のフローチャートを参照して、ユーザの移動経路を、ユーザの活動状態を表す確率的状態遷移モデルとして学習するユーザ活動モデル学習処理について説明する。
初めに、ステップS1において、GPSセンサ11は、GPSデータを取得する。即ち、GPSセンサ11は、現在地の緯度経度を測定し、測定した位置(緯度経度)と測定時刻を示すGPSデータを時系列データ記憶部12に供給する。
ステップS2において、時系列データ記憶部12は、GPSセンサ11により連続して取得されたGPSデータ、すなわち、位置と測定時刻の時系列データを記憶する。
ステップS3において、行動学習部13は、パラメータ推定処理を行い、ユーザ活動モデル学習処理は終了する。
[パラメータ推定処理]
ここで、図12を参照して、図11のステップS3のパラメータ推定処理の詳細について説明する。
ここで、図12を参照して、図11のステップS3のパラメータ推定処理の詳細について説明する。
ステップS21において、データ補間部31は、時系列データ記憶部12に記憶されているGPSデータの時系列データ、即ち、位置と測定時刻の時系列データを取得する。
ステップS22において、データ補間部31は、データ補間処理を実行する。
[データ補間処理]
ここで、図13のフローチャートを参照して、図12のステップS22のデータ補間処理の詳細について説明する。
ここで、図13のフローチャートを参照して、図12のステップS22のデータ補間処理の詳細について説明する。
ステップS51において、データ補間部31は、時系列データの欠落期間を探索する。上述したように、GPSセンサ11が人工衛星からの信号を捕捉できている間は、GPSデータは一定の時間間隔(以下、サンプリング周期ともいう)で取得されている。従って、データ補間部31は、時系列データの各データの測定時刻に基づいて、GPSデータが取得された間隔がサンプリング周期より長い期間を探索し、検出した期間を欠落期間とする。
ステップS52において、データ補間部31は、時刻データを補間する。即ち、データ補間部31は、GPSデータのサンプリング周期に基づいて、時系列データの欠落期間における測定時刻を補間する。
ステップS53において、データ補間部31は、補間した測定時刻に対応する位置データを補間する。
例えば、位置データを補間する方法として、GPSデータの欠落期間における位置を、欠落期間の直前に測定された位置(以下、消失点という)に固定する方法が考えられる。この方法は容易ではあるが、欠落期間の直後に測定された位置(以下、復活点という)に関係なく、欠落期間における位置が固定されてしまうため、補間した位置データの信頼度は低くなる。
そこで、データ補間部31は、消失点と復活点の間を線形補間することにより位置データを補間する。即ち、データ補間部31は、消失点と復活点を結ぶ直線上に、欠落期間のサンプル数と同じ数の点を等間隔に配置し、配置した点の位置を示すデータを、補間した各測定時刻に対応する位置データとして補間する。
なお、直線補間以外の方法、例えば、一次以上の多項式補間、対数補間、シンク関数を用いた補間などを用いて、位置データを補間するようにしてもよい。
ステップS54において、データ補間部31は、位置データの信頼度を設定する。具体的には、補間した位置データ(以下、補間位置データという)は、欠落期間以外の期間(以下、正常期間という)の元々存在する位置データ(以下、正常位置データともいう)と比較して、データの信頼度が低い。そこで、例えば、データ補間部31は、正常位置データと補間位置データとを後で区別できるように、2つのデータを区別するパラメータを時系列データに追加する。これにより、後の処理で正常位置データと補間位置データとを区別し、補間位置データを信頼度が低いデータとして扱うことが可能になる。
また、例えば、データ補間部31は、各位置データの値が、実際にその時刻に存在していた位置(以下、単に実際の位置ともいう)にどの程度近いかを表す信頼度として、各位置データの値を中心とする分散(あるいは標準偏差)を設定する。例えば、データ補間部31は、正常位置データの分散を0または所定の定数に設定し、補間位置データの分散を正常位置データの分散より大きい値の定数に設定する。
あるいは、補間位置データの分散を定数ではなく、その位置に応じた値に設定するようにしてもよい。具体的には、補間位置データに示される位置は、欠落期間の初めまたは終わりに近いデータほど、実際の位置との差が小さくなり、欠落期間の初めおよび終わりから遠いデータほど、実際の位置との差が大きくなると考えられる。そこで、例えば、データ補間部31は、欠落期間の初めまたは終わりに時間的に近い補間位置データほど分散を小さい値に設定し(即ち、信頼度を高く設定し)、欠落期間の初めおよび終わりから時間的に遠い補間位置データほど分散を大きい値に設定する(即ち、信頼度を低く設定する)。
なお、この信頼度として設定された分散(あるいは標準偏差)は、補間位置データに示される位置は必ずしも正しい値ではないが、その位置の周辺に、設定された分散(あるいは標準偏差)で、正しい位置が確率的に見つかることを表す。
ステップS55において、データ補間部31は、補間した時系列データを状態尤度計算部32に供給する。即ち、データ補間部31は、欠落期間の測定時刻と位置データを補間し、各位置データの信頼度を設定した時系列データを状態尤度計算部32に供給する。その後、データ補間処理は終了する。
[パラメータ推定処理の続き]
図12に戻り、ステップS23において、状態尤度計算部32は、位置データに対する各状態の尤度を計算する。具体的には、状態尤度計算部32は、ユーザ活動モデルを表すHMMの状態siへの遷移時に、時系列データの時刻tの位置データxtが出力されると仮定した場合の状態尤度P(si|xt)を、次式(5)により計算する。
図12に戻り、ステップS23において、状態尤度計算部32は、位置データに対する各状態の尤度を計算する。具体的には、状態尤度計算部32は、ユーザ活動モデルを表すHMMの状態siへの遷移時に、時系列データの時刻tの位置データxtが出力されると仮定した場合の状態尤度P(si|xt)を、次式(5)により計算する。
また、式(5)のDは、位置データの次元数を示している。いまの場合、位置データは緯度と経度の2次元なので、D=2となる。そして、xt(1)、xt(2)が、それぞれ位置データxtの緯度、経度を表すものとする。また、状態siへの遷移時に出力される位置データの緯度および経度の出力確率密度関数がそれぞれ単一正規分布に従い、μsi(1)、σsi(1)は、緯度の出力確率密度関数の中心値および標準偏差を表し、μsi(2)、σsi(2)は、経度の出力確率密度関数の中心値および標準偏差を表すものとする。
なお、式(5)は、Baum-Welchの最尤推定法で一般的に用いられる式である。
状態尤度計算部32は、全ての状態siと位置データxtの組み合わせについて、式(5)により状態尤度P(si|xt)を計算し、計算した状態尤度P(si|xt)を示すデータおよび時系列データを状態尤度修正部33に供給する。
ステップS24において、状態尤度修正部33は、補間した位置データに対する各状態の尤度を修正する。例えば、状態尤度修正部33は、補間した位置データxtに対する状態尤度P(si|xt)を、次式(6)により修正する。
式(6)は、欠落期間内の時刻tの位置データxtが出力された際に各状態siに遷移する確率が全て等しいと仮定した場合の状態尤度P(si|xt)を示しており、補間位置データに対する状態尤度P(si|xt)が全て同じ値となる。
あるいは、例えば、状態尤度修正部33は、補間した位置データxtに対する状態尤度P(si|xt)を、次式(7)により修正する。
これにより、状態尤度P(si|xt)に補間位置データの信頼度が反映され、補間位置データに対する状態尤度P(si|xt)が、正常位置データに対する状態尤度P(si|xt)より相対的に低い値に設定される。即ち、同じ値の補間位置データと正常位置データに対して、補間位置データに対する状態尤度P(si|xt)が、正常位置データに対する状態尤度P(si|xt)より低い値に設定される。また、補間位置データに対する状態尤度P(si|xt)は、補間位置データの信頼度が高いほど(分散σt(d)2が小さいほど)高くなり、補間位置データの信頼度が低いほど(分散σt(d)2が大きいほど)低くなる。
状態尤度修正部33は、全ての状態siと全ての補間した位置データxtの組み合わせについて、式(6)または式(7)により状態尤度P(si|xt)を修正する。状態尤度修正部33は、修正した状態尤度P(si|xt)を含む全ての状態尤度P(si|xt)のデータおよび時系列データをフォワード尤度計算部34に供給する。
このステップS24の処理により、結果的に、補間位置データに対する状態尤度P(si|xt)と正常位置データに対する状態尤度P(si|xt)が異なる条件で計算されることになる。なお、ステップS23において、補間位置データに対する状態尤度P(si|xt)を計算しないようにしてもよい。
ステップS25において、フォワード尤度計算部34は、フォワード尤度を計算する。即ち、フォワード尤度計算部34は、次式(8)および(9)を用いて、時刻tにおける状態siのフォワード尤度αt(si)を時刻1から最終の時刻Tまで順番に計算する。
なお、式(8)および(9)は、Baum-Welchの最尤推定法のフォワードアルゴリズムにおいて一般的に用いられる式である。
フォワード尤度計算部34は、各時刻tにおける全ての状態siのフォワード尤度αt(si)を計算する。そして、フォワード尤度計算部34は、状態尤度修正部33から取得したデータに、計算したフォワード尤度αt(si)を示すデータを追加して、バックワード尤度計算部35に供給する。
ステップS26において、バックワード尤度計算部35は、バックワード尤度を計算する。即ち、バックワード尤度計算部35は、次式(10)および(11)を用いて、時刻tにおける状態siのバックワード尤度βt(si)を、最終の時刻Tから時刻1まで逆順に計算する。
なお、式(10)および(11)は、Baum-Welchの最尤推定法のバックワードアルゴリズムにおいて一般的に用いられる式である。
バックワード尤度計算部35は、各時刻tにおける全ての状態siのバックワード尤度βt(si)を計算する。そして、バックワード尤度計算部35は、フォワード尤度計算部34から取得したデータに、計算したバックワード尤度βt(si)を示すデータを追加して、初期確率遷移確率推定部36に供給する。
このように、ステップS23乃至S26の処理により、データを補間した時系列データに対する隠れマルコフモデルの各種の尤度が計算される。
ステップS27において、初期確率遷移確率推定部36は、初期確率、状態遷移確率を更新する。即ち、初期確率遷移確率推定部36は、各状態siの初期確率πsi、各状態間の状態遷移確率aijを、次の式(12)および式(13)で求まる初期確率πsi’、状態遷移確率aij’にそれぞれ更新する。
初期確率遷移確率推定部36は、バックワード尤度計算部35から取得したデータに、更新した初期確率πsiおよび状態遷移確率aijを示すデータを追加して、観測確率推定前処理部37に供給する。
ステップS28において、観測確率推定前処理部37は、補間した位置データの観測確率への寄与度を低減する。具体的には、観測確率推定前処理部37は、欠落期間内の時刻tにおけるフォワード尤度αt(si)およびバックワード尤度βt(si)を全て0に設定する。観測確率推定前処理部37は、初期確率遷移確率推定部36から取得したデータに対して、このようにフォワード尤度αt(si)およびバックワード尤度βt(si)を変更したデータを観測確率推定部38に供給する。
ステップS29において、観測確率推定部38は、観測確率を更新する。即ち、観測確率推定部38は、各状態siの出力確率密度関数の中心値μsi(d)、分散σsi(d)2を、次の式(14)および式(15)で求まる中心値μsi(d)’、分散σsi(d)’2にそれぞれ更新する。
なお、式(14)および式(15)は、Baum-Welchの最尤推定法で一般的に用いられる式である。
ここで、ステップS28において、欠落期間内の時刻tにおけるフォワード尤度αt(si)およびバックワード尤度βt(si)が全て0に設定されている。従って、特別な計算式を用いることなく、各状態siの観測確率(中心値μsi(d)、分散σsi(d)2)への補間位置データの寄与度が、正常位置データの寄与度より低くなり、ほぼ0になる。これにより、観測確率の推定に、信頼度が低い補間位置データがほとんど反映されなくなり、観測確率の推定精度が向上する。
ステップS30において、観測確率推定部38は、パラメータの更新を終了するか否かを判定する。観測確率推定部38は、後述する所定の更新終了条件がまだ成立していない場合、パラメータの更新を終了しないと判定し、処理はステップS31に進む。
ステップS31において、観測確率推定部38は、更新されたパラメータを状態尤度計算部32に供給する。即ち、観測確率推定部38は、ステップS27乃至S29の処理で更新された、各状態siの初期確率πsi、中心値μsi(d)および分散σsi(d)2、並びに、各状態間の状態遷移確率aijを示すデータを状態尤度計算部32に供給する。
その後、ステップS30において、パラメータの更新を終了すると判定されるまで、ステップS23乃至S31の処理が繰り返し実行される。これにより、状態siの系列の各種の尤度、すなわち、状態尤度P(si|xt)、フォワード尤度αt(si)、バックワード尤度βt(si))が次第に増加し、最終的に最大になるように、HMMのパラメータの更新が行われる。
一方、ステップS30において、観測確率推定部38は、例えば、各尤度の増分が所定の値以下となり、パラメータの更新の収束条件を満たした場合、パラメータの更新を終了すると判定し、処理はステップS32に進む。あるいは、観測確率推定部38は、例えば、ステップS23乃至S31の更新処理が所定の回数実行された場合、パラメータの更新を終了すると判定し、処理はステップS32に進む。
ステップS32において、観測確率推定部38は、最終的なパラメータを状態尤度計算部32に供給する。即ち、観測確率推定部38は、ステップS23乃至S31の更新処理で最終的に求められた、各状態siの初期確率πsi、中心値μsi(d)および分散σsi(d)2、並びに、各状態間の状態遷移確率aijを示すデータを、行動認識部14および行動予測部15に供給する。
その後、パラメータ推定処理は終了する。
[到達時間予測処理]
次に、図11のユーザ活動モデル学習処理により得られたユーザ活動モデルを表す確率的状態遷移モデルのパラメータを用いて、目的地までの経路を探索し、到達時間を算出してユーザに提示する到達時間予測処理について説明する。
次に、図11のユーザ活動モデル学習処理により得られたユーザ活動モデルを表す確率的状態遷移モデルのパラメータを用いて、目的地までの経路を探索し、到達時間を算出してユーザに提示する到達時間予測処理について説明する。
図14は、到達時間予測処理のフローチャートを示している。なお、この例では、図14の処理の前に、目的地が予め決められているものとするが、図14の処理の途中で目的地を入力するようにしてもよい。
初めに、ステップS101において、GPSセンサ11は、位置の時系列データを取得し、行動認識部14に供給する。行動認識部14には、所定のサンプル数の位置の時系列データが一時的に記憶される。
ステップS102において、行動認識部14は、学習により得られたパラメータに基づくユーザ活動モデルから、ユーザの現在の活動状態を認識する。即ち、行動認識部14は、ユーザの現在地を認識する。そして、行動認識部14は、ユーザの現在の状態ノードのノード番号を行動予測部15に供給する。
ステップS103において、行動予測部15は、現在探索している状態ノード(以下、適宜、現在の状態ノードともいう。)に対応する点が、端点、通過点、分岐点、またはループのいずれであるかを判定する。ステップS102の処理直後は、ユーザの現在地に対応する状態ノードが現在の状態ノードとなる。
ステップS103で、現在の状態ノードに対応する点が端点であると判定された場合、処理はステップS104に進み、行動予測部15は、現在の状態ノードをここまでの経路に接続し、この経路の探索を終了し、ステップS111に進む。なお、現在の状態ノードが現在地に対応する状態ノードである場合、ここまでの経路は存在しないため、接続する処理は行われない。ステップS105,S107、S110でも同様である。
ステップS103で、現在の状態ノードに対応する点が通過点であると判定された場合、処理はステップS105に進み、行動予測部15は、現在の状態ノードをここまでの経路に接続する。そして、ステップS106において、行動予測部15は、次の状態ノードを現在の状態ノードとし、移動する。ステップS106の処理後、処理はステップS103に戻る。
ステップS103で、現在の状態ノードに対応する点が分岐点であると判定された場合、処理はステップS107に進み、行動予測部15は、現在の状態ノードをここまでの経路と接続する。そして、ステップS108において、行動予測部15は、ここまでの経路を分岐の数だけ複製し、分岐先の状態ノードと接続する。さらにステップS109において、行動予測部15は、複製した経路の1つを選択し、その選択した経路の先の状態ノードを現在の状態ノードとし移動する。ステップS109の処理後、処理はステップS103に戻る。
一方、ステップS103で、現在の状態ノードに対応する点がループであると判定された場合、処理はステップS110に進み、行動予測部15は、現在の状態ノードをここまでの経路と接続せずに、この経路の探索を終了し、ステップS111に進む。
ステップS111において、行動予測部15は、未探索の経路があるかを判定する。ステップS111で、未探索の経路があると判定された場合、処理はステップS112に進み、行動予測部15は、現在地の状態ノードに戻り、未探索の経路の次の状態ノードを現在の状態ノードとし移動する。ステップS112の処理後、処理はステップS103に戻る。これにより、未探索の経路について、端点またはループにより探索が終了するまで経路の探索が実行される。
ステップS111で、未探索の経路がないと判定された場合、処理はステップS113に進み、行動予測部15は、探索された各経路の選択確率(生起確率)を計算する。行動予測部15は、各経路と、その選択確率を到達時間予測部16に供給する。
ステップS114において、到達時間予測部16は、行動予測部15によって探索された経路から、入力された目的地を含む経路を抽出し、目的地の到達確率を算出する。具体的には、到達時間予測部16は、目的地に複数の経路が存在する場合、その複数の経路の選択確率の和を、目的地の到達確率として計算する。目的地への経路が1つしかない場合、経路の選択確率が、そのまま、目的地の到達確率とされる。
ステップS115において、到達時間予測部16は、抽出された経路の数が提示数として予め設定された所定個数より多いかを判定する。
ステップS115で、抽出された経路の数が所定個数より多いと判定された場合、処理はステップS116に進み、到達時間予測部16は、表示部18に表示する所定個数の経路を決定する。例えば、到達時間予測部16は、選択される可能性の高い順に、所定個数の経路を決定することができる。
一方、ステップS115で、抽出された経路の数が所定個数以下であると判定された場合、ステップS116の処理がスキップされる。即ち、この場合、目的地に到達するための全ての経路が表示部18に表示される。
ステップS117において、到達時間予測部16は、表示部18に表示すると決定された各経路の到達時間を計算する。そして、到達時間予測部16は、目的地の到達確率および目的地までの経路と到達時間を表示する画像の信号を表示部18に供給する。
ステップS118において、表示部18は、到達時間予測部16から供給される画像の信号に基づいて、目的地の到達確率および目的地までの経路と到達時間を表示して、処理を終了する。
[本発明の効果の説明]
以上のように、予測システム1では、GPSセンサ11により取得される位置と測定時刻の時系列データの欠落期間のデータを補間して、補間した時系列データを用いて、ユーザの活動状態を隠れマルコフモデルとして学習する学習処理を行う。従って、時系列データに欠落期間が存在していても、正確にユーザの活動状態を学習することができる。
以上のように、予測システム1では、GPSセンサ11により取得される位置と測定時刻の時系列データの欠落期間のデータを補間して、補間した時系列データを用いて、ユーザの活動状態を隠れマルコフモデルとして学習する学習処理を行う。従って、時系列データに欠落期間が存在していても、正確にユーザの活動状態を学習することができる。
ここで、図15を参照して、予測システム1の学習処理の効果の具体例を説明する。
図15は、図1および図2を参照して上述した実験と同じ時系列データを用いて、予測システム1により学習処理を行った結果を示す図である。なお、このときの学習処理では、図12のステップS24における状態尤度P(si|xt)の修正に、式(6)を用いている。
図15と図1を比較すると、図15では、距離が離れた状態ノードを短絡する状態遷移の数が減っていることが分かる。
なお、式(6)により、欠落期間の補間位置データに対する状態尤度P(si|xt)を全て同じ値にしても、フォワード尤度αt(si)およびバックワード尤度βt(si)を計算する際に、補間した位置データに対して、その近くにある状態ノードが自動的に割り当てられるため、状態ノードを短絡する状態遷移が学習されることが防止される。
また、状態尤度P(si|xt)の修正に上述した式(7)を用いても、ほぼ同じ学習結果を得ることができる。ただし、式(7)を用いた方が、より状態尤度P(si|xt)の信頼度が向上するため、より正確にユーザ活動モデルを学習することができる。
なお、参考までに、図16乃至図18を参照して、同じ時系列データを用いて、予測システム1とは異なる対策を行った場合の学習結果の例について説明する。
図16は、時系列データの欠落期間において、消失点と同じ位置のデータを補間した時系列データを用いて学習処理を行った結果を示している。即ち、欠落期間中、消失点と同じ位置でホールドされた時系列データを用いて学習処理を行った結果を示している。この方法を用いれば、ホールドされた位置付近に複数の状態ノードが割り当てられたり、ホールドされた位置付近の状態ノードの自己遷移確率が高くなったりする。従って、目的地までのステップ数の不足の問題は解消され、目的地への到達時間の予測精度は向上する。しかし、図1と比較して、短絡経路の数を削減することはできず、処理負荷が増大する問題は解消されない。
図17は、欠落期間が所定の時間以上になったら、時系列データを分割し、別の系統の時系列データとするようにした場合の学習結果を示している。ちなみに、図17は、欠落期間が2分以上になったら時系列データを分割するようにした場合の例を示している。この場合、学習に用いる時系列データに欠落期間がなくなるため、図17に示されるように、通勤経路が一本の経路として学習されることが分かる。ただし、この場合、時系列データを分割する条件となる時間を適正に設定する必要がある。
例えば、図18は、欠落期間が1分以上になったら、時系列データを分割するようにした場合の学習結果を示している。この場合、時系列データを分割するタイミングが多くなりすぎて、その結果、通勤経路が分断され、1本につながっていない。このように経路が分断されてしまうと、例えば、目的地への到達時間を予測する場合に、目的地に到達する経路を1つも見つけられずに処理が終わってしまう可能性がある。
これに対して、予測システム1では、各種のパラメータの調整を行う必要がなく、同じ時系列データに対して、常に同じ学習結果を得ることができる。
また、予測システム1は、学習処理により得られたパラメータで表される確率的状態遷移モデルを用いて、入力された目的地の到達確率および目的地までの経路と到達時間を予測して、ユーザに提示する。
したがって、予測システム1によれば、ユーザが指定した目的地の到達確率および目的地までの経路とその到達時間を予測して、ユーザに提示することができる。また、上述したように、予測システム1では、時系列データが欠落していても、正確にユーザの活動モデルが学習されるため、それらの予測精度が向上する。
<2.変形例>
以上の説明では、1系統の時系列データに対する学習処理を中心に説明したが、本発明は、複数系統の時系列データに対する学習処理にも適用することができる。
以上の説明では、1系統の時系列データに対する学習処理を中心に説明したが、本発明は、複数系統の時系列データに対する学習処理にも適用することができる。
また、以上の説明では、ユーザの活動モデルの学習に、位置の時系列データのみを用いる例を示したが、測定時刻の時系列データも用いるようにしてもよい。ただし、測定時刻は、位置データと異なり、補間した値に誤差が生じない。従って、例えば、補間位置データに対する状態尤度P(si|xt)の計算式(7)を、次式(16)に変更するようにすればよい。
式(16)を用いることにより、補間した緯度および経度の信頼度のみが状態尤度P(si|xt)に反映されるようになる。
なお、式(16)の代わりに、式(6)を用いることも可能である。
さらに、以上の説明では、位置の時系列データを用いて、ユーザの活動モデルを表すHMMを学習する例を示したが、本発明は、データの種類を問わず、時系列データからHMMを学習する場合に適用できる。
また、学習に用いる時系列データの欠落期間以外の期間(正常期間)において、データのサンプリング間隔が一定でない場合、サンプリング間隔が一定になるように、再サンプリング処理等を行うようにしてもよい。
上述した一連の処理は、ハードウエアにより実行することもできるし、ソフトウエアにより実行することもできる。一連の処理をソフトウエアにより実行する場合には、そのソフトウエアを構成するプログラムが、コンピュータにインストールされる。ここで、コンピュータには、専用のハードウエアに組み込まれているコンピュータや、各種のプログラムをインストールすることで、各種の機能を実行することが可能な、例えば汎用のパーソナルコンピュータなどが含まれる。
図19は、上述した一連の処理をプログラムにより実行するコンピュータのハードウエアの構成例を示すブロック図である。
コンピュータにおいて、CPU(Central Processing Unit)201,ROM(Read Only Memory)202,RAM(Random Access Memory)203は、バス204により相互に接続されている。
バス204には、さらに、入出力インタフェース205が接続されている。入出力インタフェース205には、入力部206、出力部207、記憶部208、通信部209、ドライブ210、およびGPSセンサ211が接続されている。
入力部206は、キーボード、マウス、マイクロホンなどよりなる。出力部207は、ディスプレイ、スピーカなどよりなる。記憶部208は、ハードディスクや不揮発性のメモリなどよりなる。通信部209は、ネットワークインタフェースなどよりなる。ドライブ210は、磁気ディスク、光ディスク、光磁気ディスク、或いは半導体メモリなどのリムーバブル記録媒体212を駆動する。GPSセンサ211は、図1のGPSセンサ11に対応する。
以上のように構成されるコンピュータでは、CPU201が、例えば、記憶部208に記憶されているプログラムを、入出力インタフェース205及びバス204を介して、RAM203にロードして実行することにより、上述した一連の処理が行われる。
コンピュータ(CPU201)が実行するプログラムは、例えば、パッケージメディア等としてのリムーバブル記録媒体212に記録して提供することができる。また、プログラムは、ローカルエリアネットワーク、インターネット、デジタル衛星放送といった、有線または無線の伝送媒体を介して提供することができる。
コンピュータでは、プログラムは、リムーバブル記録媒体212をドライブ210に装着することにより、入出力インタフェース205を介して、記憶部208にインストールすることができる。また、プログラムは、有線または無線の伝送媒体を介して、通信部209で受信し、記憶部208にインストールすることができる。その他、プログラムは、ROM202や記憶部208に、あらかじめインストールしておくことができる。
なお、コンピュータが実行するプログラムは、本明細書で説明する順序に沿って時系列に処理が行われるプログラムであっても良いし、並列に、あるいは呼び出しが行われたとき等の必要なタイミングで処理が行われるプログラムであっても良い。
なお、本明細書において、フローチャートに記述されたステップは、記載された順序に沿って時系列的に行われる場合はもちろん、必ずしも時系列的に処理されなくとも、並列に、あるいは呼び出しが行われたとき等の必要なタイミングで実行されてもよい。
なお、本明細書において、システムとは、複数の装置により構成される装置全体を表すものである。
本発明の実施の形態は、上述した実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において種々の変更が可能である。
1 予測システム, 13 行動学習部, 14 行動認識部, 15 行動予測部, 16 到達時間予測部, 31 データ補間部, 32 状態尤度計算部, 33 状態尤度修正部, 34 フォワード尤度計算部, 35 バックワード尤度計算部, 36 初期確率遷移確率推定部, 37 観測確率推定前処理部, 38 観測確率推定部
Claims (7)
- 時系列データにおいて欠落しているデータを補間する補間手段と、
前記時系列データから隠れマルコフモデルを推定する推定手段と、
推定された前記隠れマルコフモデルの尤度を計算する尤度計算手段と
を備え、
前記尤度計算手段は、欠落していないデータである正常データに対する尤度と補間されたデータである補間データに対する尤度とを異なる条件で計算して、データが補間された前記時系列データに対する前記隠れマルコフモデルの尤度を計算し、
前記推定手段は、前記尤度計算手段により計算される尤度が大きくなるように、前記隠れマルコフモデルを更新する
学習装置。 - 前記尤度計算手段は、前記補間データに対する尤度を前記正常データに対する尤度より低くする
請求項1に記載の学習装置。 - 前記補間手段は、前記時系列データの欠落期間の直前の第1のデータと直後の第2のデータに基づいて、前記欠落期間におけるデータを補間するとともに、補間したデータに対する信頼度を、前記欠落期間の初めまたは終わりに近いデータほど高く設定し、前記欠落期間の初めおよび終わりから遠いデータほど低く設定し、
前記尤度計算手段は、前記補間データの前記信頼度が低いほど、前記補間データに対する尤度を低くする
請求項2に記載の学習装置。 - 前記尤度計算手段は、前記補間データに対する前記隠れマルコフモデルの各状態の尤度を全て同じ値にする
請求項1に記載の学習装置。 - 前記推定手段は、前記隠れマルコフモデルの観測確率への前記補間データの寄与度を前記正常データの寄与度より低くする
請求項1乃至4のいずれかに記載の学習装置。 - 時系列データから隠れマルコフモデルを学習する学習装置が、
前記時系列データにおいて欠落しているデータを補間し、
前記時系列データから前記隠れマルコフモデルを推定し、
推定された前記隠れマルコフモデルの尤度を計算するとともに、
欠落していないデータである正常データに対する尤度と補間されたデータである補間データに対する尤度とを異なる条件で計算して、データが補間された前記時系列データに対する前記隠れマルコフモデルの尤度を計算し、
計算される尤度が大きくなるように、前記隠れマルコフモデルを更新する
ステップを含む学習方法。 - 時系列データにおいて欠落しているデータを補間し、
前記時系列データから隠れマルコフモデルを推定し、
推定された前記隠れマルコフモデルの尤度を計算するとともに、
欠落していないデータである正常データに対する尤度と補間されたデータである補間データに対する尤度とを異なる条件で計算して、データが補間された前記時系列データに対する前記隠れマルコフモデルの尤度を計算し、
計算される尤度が大きくなるように、前記隠れマルコフモデルを更新する
ステップを含む処理をコンピュータに実行させるためのプログラム。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009276691A JP2011118755A (ja) | 2009-12-04 | 2009-12-04 | 学習装置、学習方法、およびプログラム |
US12/917,853 US20110137831A1 (en) | 2009-12-04 | 2010-11-02 | Learning apparatus, learning method and program |
CN201010570580XA CN102087709A (zh) | 2009-12-04 | 2010-11-26 | 学习设备、学习方法以及程序 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009276691A JP2011118755A (ja) | 2009-12-04 | 2009-12-04 | 学習装置、学習方法、およびプログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2011118755A true JP2011118755A (ja) | 2011-06-16 |
Family
ID=44082979
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009276691A Withdrawn JP2011118755A (ja) | 2009-12-04 | 2009-12-04 | 学習装置、学習方法、およびプログラム |
Country Status (3)
Country | Link |
---|---|
US (1) | US20110137831A1 (ja) |
JP (1) | JP2011118755A (ja) |
CN (1) | CN102087709A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015060237A (ja) * | 2013-09-17 | 2015-03-30 | 日本電気株式会社 | 予測モデル学習装置、予測モデル学習方法およびコンピュータプログラム |
JPWO2018131311A1 (ja) * | 2017-01-10 | 2019-11-07 | 日本電気株式会社 | センシングシステム、センサノード装置、センサ測定値処理方法及びプログラム |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013192533A2 (en) | 2012-06-22 | 2013-12-27 | Google, Inc. | Contextual traffic or transit alerts |
KR101630389B1 (ko) | 2012-06-22 | 2016-06-15 | 구글 인코포레이티드 | 현재 위치 또는 시간에 관한 정보의 제시 |
US20140278838A1 (en) | 2013-03-14 | 2014-09-18 | Uber Technologies, Inc. | Determining an amount for a toll based on location data points provided by a computing device |
US9965521B1 (en) * | 2014-02-05 | 2018-05-08 | Google Llc | Determining a transition probability from one or more past activity indications to one or more subsequent activity indications |
CN106233208B (zh) * | 2014-05-12 | 2017-09-29 | 三菱电机株式会社 | 时间序列数据处理装置 |
US10424036B2 (en) | 2014-06-02 | 2019-09-24 | Uber Technologies, Inc. | Maintaining data for use with a transport service during connectivity loss between systems |
CN105279170A (zh) * | 2014-06-27 | 2016-01-27 | 华为技术有限公司 | 一种活动识别方法及系统 |
CN105446995A (zh) * | 2014-07-14 | 2016-03-30 | 富士通株式会社 | 估计数据更新时刻的方法和设备、数据集成方法和设备 |
US9503516B2 (en) | 2014-08-06 | 2016-11-22 | Google Technology Holdings LLC | Context-based contact notification |
WO2018227389A1 (en) | 2017-06-13 | 2018-12-20 | Beijing Didi Infinity Technology And Development Co., Ltd. | Systems and methods for determining estimated time of arrival |
KR102481794B1 (ko) * | 2018-04-18 | 2022-12-28 | 한국전자통신연구원 | 사용자 행동의 시계열 분석을 활용한 위험 상황 인지 방법 및 장치 |
KR20190136578A (ko) * | 2018-05-31 | 2019-12-10 | 삼성전자주식회사 | 음성 인식 방법 및 장치 |
KR102592830B1 (ko) * | 2018-12-05 | 2023-10-23 | 현대자동차주식회사 | 차량용 센서퓨전 타겟 예측 장치 및 그의 센서 퓨전 타겟 예측 방법과 그를 포함하는 차량 |
KR102501530B1 (ko) * | 2018-12-31 | 2023-02-21 | 한국전자통신연구원 | 시계열 데이터 처리 장치 및 이의 동작 방법 |
CN112580837A (zh) * | 2020-12-29 | 2021-03-30 | 广州航天海特系统工程有限公司 | 一种设备检修方法、装置、设备及存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19709518C5 (de) * | 1997-03-10 | 2006-05-04 | Harman Becker Automotive Systems Gmbh | Verfahren und Vorrichtung zur Spracheingabe einer Zieladresse in ein Zielführungssystem im Echtzeitbetrieb |
CN101549498B (zh) * | 2009-04-23 | 2010-12-29 | 上海交通大学 | 智能帮扶式助行机器人的自主跟踪与导航系统 |
-
2009
- 2009-12-04 JP JP2009276691A patent/JP2011118755A/ja not_active Withdrawn
-
2010
- 2010-11-02 US US12/917,853 patent/US20110137831A1/en not_active Abandoned
- 2010-11-26 CN CN201010570580XA patent/CN102087709A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015060237A (ja) * | 2013-09-17 | 2015-03-30 | 日本電気株式会社 | 予測モデル学習装置、予測モデル学習方法およびコンピュータプログラム |
JPWO2018131311A1 (ja) * | 2017-01-10 | 2019-11-07 | 日本電気株式会社 | センシングシステム、センサノード装置、センサ測定値処理方法及びプログラム |
US11514277B2 (en) | 2017-01-10 | 2022-11-29 | Nec Corporation | Sensing system, sensor node device, sensor measurement value processing method, and program |
JP7206915B2 (ja) | 2017-01-10 | 2023-01-18 | 日本電気株式会社 | センシングシステム、センサノード装置、センサ測定値処理方法及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
US20110137831A1 (en) | 2011-06-09 |
CN102087709A (zh) | 2011-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2011118755A (ja) | 学習装置、学習方法、およびプログラム | |
US9589082B2 (en) | Data processing device that calculates an arrival probability for a destination using a user's movement history including a missing portion | |
US8572008B2 (en) | Learning apparatus and method, prediction apparatus and method, and program | |
US20110137833A1 (en) | Data processing apparatus, data processing method and program | |
JP5495014B2 (ja) | データ処理装置、データ処理方法、およびプログラム | |
US9285235B2 (en) | Information processing device, information processing method, and program | |
JP2011252844A (ja) | データ処理装置、データ処理方法、およびプログラム | |
US20110313957A1 (en) | Data processing apparatus, data processing method and program | |
CN106714110B (zh) | 一种Wi-Fi位置指纹地图自动构建方法及系统 | |
US8665743B2 (en) | User behavior model and statistical transition map to assist advanced WLAN applications | |
WO2017031856A1 (zh) | 信息预测的方法和装置 | |
US8457880B1 (en) | Telematics using personal mobile devices | |
JP2012003494A (ja) | 情報処理装置、情報処理方法及びプログラム | |
KR101058296B1 (ko) | 모바일 디바이스에 대한 실시간 위치 추적 방법 및 시스템 | |
JP2012008771A (ja) | 情報処理装置、情報処理システム、情報処理方法およびプログラム | |
JP2013205170A (ja) | 情報処理装置、情報処理方法、およびプログラム | |
JP2011059924A5 (ja) | ||
CN105120433A (zh) | 基于连续采样及模糊聚类处理的wlan室内定位方法 | |
JP2013205171A (ja) | 情報処理装置、情報処理方法、およびプログラム | |
JP5845604B2 (ja) | 情報処理装置、情報処理方法、およびプログラム | |
JP5664398B2 (ja) | 情報処理装置、情報処理方法、およびプログラム | |
JP2008294593A (ja) | 位置推定システム、位置推定装置及びそれらに用いる位置推定方法並びにそのプログラム | |
CN114646313A (zh) | 一种用户轨迹定位方法、电子设备及计算机存储介质 | |
US20230114585A1 (en) | Ranging-type positioning system and ranging-type positioning method based on crowdsourced calibration | |
KR101203200B1 (ko) | 무선 센서 네트워크에서의 노드 위치 인식 방법 및 그 방법을 이용한 노드 위치 인식 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20130205 |