JP2019074400A - 位置推定装置、空調システム、位置推定方法及びプログラム - Google Patents

Abstract

【課題】空間内の音源の位置を推定する位置推定装置を提供する。【解決手段】位置推定装置は、空間に設けられた複数のマイクが、空間において発された音を検出した各々の時刻を取得し、取得した各々の前記時刻に基づいて、複数のマイクのうち基準となる１つのマイクが音を検出した時刻と他の各々のマイクが検出した時刻との差である到達時間差を算出し、予め定められた到達時間差および所定の音源位置の関係を定めたルックアップテーブルと、実際に算出した到達時間差とに基づいて、音の発信位置を推定する。【選択図】図２

Description

本発明は、位置推定装置、空調システム、位置推定方法及びプログラムに関する。

オフィス等の空間を空調する場合、その空間の室温等が均一になるように制御することが多い。これに対し、例えば、ユーザが存在するエリアを対象として室温等の空調制御を行う方法が提供されている。このような制御を行う場合、空調システムは、ユーザが存在する位置を正確に推定する必要がある。従来から位置の推定方法には、様々な方法が提供されている。

例えば、位置が既知である４つのＧＰＳ衛星から受信した電波の到達時間によって、受信機の位置を推定する方法が知られている。ここで、４つの衛星の座標をそれぞれ（Ａ１、Ｂ１、Ｃ１）、（Ａ２、Ｂ２、Ｃ２）、（Ａ３、Ｂ３、Ｃ３）、（Ａ４、Ｂ４、Ｃ４）とし、受信機の位置の座標を（ｘ、ｙ、ｚ）とし、受信機の時間のずれをｄ、電波の速度を光速ｃとすると、以下の連立方程式が得られる。
ｆ１＝（ｘ−Ａ１）^２＋（ｙ−Ｂ１）^２＋（ｚ−Ｃ１）^２−（ｃ（ｔ１−ｄ））^２＝０
ｆ２＝（ｘ−Ａ２）^２＋（ｙ−Ｂ２）^２＋（ｚ−Ｃ２）^２−（ｃ（ｔ２−ｄ））^２＝０
ｆ３＝（ｘ−Ａ３）^２＋（ｙ−Ｂ３）^２＋（ｚ−Ｃ３）^２−（ｃ（ｔ３−ｄ））^２＝０
ｆ４＝（ｘ−Ａ４）^２＋（ｙ−Ｂ４）^２＋（ｚ−Ｃ４）^２−（ｃ（ｔ４−ｄ））^２＝０
受信機の位置はこの方程式を解くことにより得られるが、一般にこの方程式の解は、ニュートン−ラプソン法等の近似計算を用いて求めることが多い。この近似計算は、処理負荷が高く、同様の計算を空調機の制御に適用しようとしても、空調機に搭載できるコンピュータの性能を超える可能性がある。また、近似計算によって解を導出する場合、探索の初期値によっては、局所解に陥ってしまうことがある。

なお、特許文献１には、人型ロボットが人間と会話によるコミュニケーションを行うために音声が発せられた方向を推定する方法が記載されている。

特表２０１６−５３７６２２号公報

空調機に搭載されるコンピュータでも無理なく処理できる計算負荷で、局所解に陥ることなく、ユーザが存在する位置を推定できる方法が求められている。

そこでこの発明は、上述の課題を解決することのできる位置推定装置、空調システム、位置推定方法及びプログラムを提供することを目的としている。

本発明の一態様は、空間に設けられた複数のマイクが、前記空間において発された音を検出した各々の時刻を取得する時刻情報取得部と、取得した各々の前記時刻に基づいて、複数の前記マイクのうち基準となる１つの前記マイクが前記音を検出した時刻と他の各々のマイクが検出した時刻との差である到達時間差を算出する到達時間差算出部と、予め定められた前記到達時間差および所定の音源位置の関係を定めたルックアップテーブルと、実際に前記到達時間差算出部が算出した前記到達時間差とに基づいて、前記音の発信位置を推定する位置推定部と、を備える位置推定装置である。

本発明の一態様における前記位置推定部は、前記到達時間差算出部が算出した前記到達時間差と、前記ルックアップテーブルに基づく前記到達時間差との差が最小となるときの前記ルックアップテーブルに基づく前記音源位置を、前記音の発信位置として推定する。

本発明の一態様における前記ルックアップテーブルには、複数の前記マイクが配置された位置の重心を原点とする座標系で表された前記音源位置と前記到達時間差との関係が記録されている。

本発明の一態様における前記ルックアップテーブルには、前記マイクが配置された空間の環境条件に応じた前記到達時間差および前記音源位置の関係が記録されている。

本発明の一態様における前記空間の環境条件は、温度および湿度および気圧のうち少なくとも一つである。

本発明の一態様における前記ルックアップテーブルには、所定の温度における前記到達時間差および前記音源位置の関係が記録され、前記位置推定部は、前記ルックアップテーブルに基づく前記到達時間差を実際の温度に基づいて補正し、補正後の前記ルックアップテーブルと前記到達時間差算出部が算出した前記到達時間差とに基づいて、前記音の発信位置を推定する。

本発明の一態様における前記位置推定装置は、前記空間の内部を移動しながら発された前記音について前記位置推定部が推定した前記音の発信位置に基づいて、前記空間の形状を示す地図情報を作成する地図情報作成部、をさらに備える。

本発明の一態様は、上記の位置推定装置を備え、前記位置推定装置が推定した前記音の発信位置を対象として空調を行う、空調システムである。

本発明の一態様は、空間に設けられた複数のマイクが、前記空間において発された音を検出した各々の時刻を取得するステップと、取得した各々の前記時刻に基づいて、複数の前記マイクのうち基準となる１つの前記マイクが前記音を検出した時刻と他の各々のマイクが検出した時刻との差である到達時間差を算出するステップと、予め定められた前記到達時間差および所定の音源位置の関係を定めたルックアップテーブルと、実際に前記到達時間差を算出するステップで算出した前記到達時間差とに基づいて、前記音の発信位置を推定するステップと、を有する位置推定方法である。

本発明の一態様は、コンピュータを、空間に設けられた複数のマイクが、前記空間において発された音を検出した各々の時刻を取得する手段、取得した各々の前記時刻に基づいて、複数の前記マイクのうち基準となる１つの前記マイクが前記音を検出した時刻と他の各々のマイクが検出した時刻との差である到達時間差を算出する手段、予め定められた前記到達時間差および所定の音源位置の関係を定めたルックアップテーブルと、実際に前記到達時間差を算出する手段が算出した前記到達時間差とに基づいて、前記音の発信位置を推定する手段、として機能させるためのプログラムである。

本発明によれば、少ない計算負荷で、正確に位置を推定することができる。

本発明の一実施形態における空調システムの一例を示す図である。 本発明の一実施形態における位置推定装置の一例を示すブロック図である。 本発明の一実施形態におけるマイクの配置例を示す図である。 本発明の一実施形態における到達時間差を説明する第１の図である。 本発明の一実施形態における到達時間差を説明する第２の図である。 本発明の一実施形態における到達時間差を説明する第３の図である。 本発明の一実施形態におけるＬＵＴの一例を示す図である。 本発明の一実施形態における位置推定結果の一例を示す第１の図である。 本発明の一実施形態における位置推定結果の一例を示す第２の図である。 本発明の一実施形態における位置推定処理の一例を示すフローチャートである。 本発明の一実施形態における地図情報の一例を示す図である。 本発明の一実施形態における地図情報の作成処理の一例を示すフローチャートである。 本発明の一実施形態における位置推定装置のハードウェア構成の一例を示す図である。

＜実施形態＞
以下、本発明の一実施形態による冷媒回路を図１〜図１３を参照して説明する。
図１は、本発明の一実施形態における空調システムの一例を示す図である。
図示するように空調システムは、空調機１０と、端末装置４０と、マイクＭ１〜Ｍ４とを含む。空調機１０は、制御装置２０を備える。制御装置２０は、空調機１０の運転を制御する。制御装置２０は、空間ｗの一部のエリアを対象として、空調を行うよう空調機１０を制御する機能を有している。制御装置２０は、位置推定装置３０を備える。位置推定装置３０は、端末装置４０が発した音ｓに基づいて、音ｓが発された位置を推定する。制御装置２０は、例えば、位置推定装置３０が推定した位置を対象として空調制御を行う。端末装置４０は、所定の周波数の音ｓを出力する機能を有している。また、端末装置４０は、端末装置４０が存在する空間ｗにおけるエリアを指定する情報ｋを、空調機１０（制御装置２０）へ送信する機能を有している。端末装置４０は、例えば、スマートフォン等の情報処理装置である。あるいは、端末装置４０は、音を出力する装置と、情報ｋを出力する装置の２台で構成されていてもよい。マイクＭ１〜Ｍ４は、端末装置４０が発した音を検出する。

位置推定装置３０は、端末装置４０が発した音ｓをマイクＭ１〜Ｍ４の各々が検出した時刻の差に基づいて、その音源（端末装置４０）の位置を推定し、制御装置２０は、位置推定装置３０が推定した位置を対象として空調制御を行う。

また、位置推定装置３０は、音源の位置推定を利用して空間ｗの地図情報を作成する機能を有している。制御装置２０は、位置推定装置３０が作成した地図情報に対して、端末装置４０を介してユーザが指定した空間ｗにおけるエリアを指定する情報ｋを取得し、情報ｋが指定するエリアを対象とする空調制御を行う。次に位置推定装置３０の機能について説明する。

図２は、本発明の一実施形態における位置推定装置の一例を示すブロック図である。
位置推定装置３０は、音情報取得部３１と、到達時間差算出部３２と、位置推定部３３と、センサ情報取得部３４と、地図情報作成部３５と、記憶部３６を備える。
音情報取得部３１は、マイクＭ１〜Ｍ４が音ｓを検出して電気信号に変換した音信号をマイクＭ１〜Ｍ４別に取得する。音情報取得部３１は、マイクＭ１〜Ｍ４と各マイクが音ｓを検出した時刻（音情報取得部３１が音信号を取得した時刻）とを対応付けて記憶部３６に記録する。

到達時間差算出部３２は、マイクＭ１〜Ｍ４のうちの１つのマイク（例えばマイクＭ１）を基準として、基準となるマイクＭ１が音ｓを検出した時刻と、他の３つのマイクＭ２〜Ｍ４が音ｓを検出した時刻との差、つまり、音ｓがマイクＭ１に到達する時刻と、マイクＭ２〜Ｍ４に到達する時刻との時間差（到達時間差）を算出する。

位置推定部３３は、到達時間差算出部３２が算出した到達時間差と記憶部３６が記憶するルックアップテーブル３７（以下、ＬＵＴ３７と記載する）とに基づいて、音ｓが発された位置を推定する。ＬＵＴ３７は、マイクＭ１〜Ｍ４のうち２つのマイクの間での到達時間差と空間ｗにおける所定の音源位置との対応関係を規定したテーブルを含む。

センサ情報取得部３４は、温度センサ、湿度センサ、気圧センサから、空間ｗの温度、湿度、気圧の情報を取得する。温度センサ、湿度センサ、気圧センサは、空調機１０が備えるものでもよいし、空間ｗに設けられたものでもよい。また、これらのセンサは、全て設けられている必要は無く、何れか１つでよい。

地図情報作成部３５は、位置推定部３３が推定した位置情報を用いて、空間ｗの地図情報を作成する。より具体的には、地図情報作成部３５は、端末装置４０を空間ｗの形状に沿って移動させながら出力した音ｓについて位置推定部３３が推定した位置情報を用いて、空間ｗの地図情報を作成する。地図情報作成部３５は、位置推定装置３０とは独立した装置として設けてもよい。

記憶部３６は、ＬＵＴ３７、位置推定部３３が推定した端末装置４０の位置情報、地図情報作成部３５が作成した地図情報などを記憶する。記憶部３６は、位置推定装置３０が備える記憶装置でもよいし、外部（例えば、データセンタ等）の記憶装置であってもよい。

次に図３〜図９を用いて本実施形態における位置推定方法を説明する。
図３は、本発明の一実施形態におけるマイクの配置例を示す図である。
図３にマイクＭ１〜Ｍ４の配置例を示す。マイクＭ１〜Ｍ４は、空間ｗの天井に、例えばマイクＭ１〜Ｍの各々が、空調機１０の吹き出し口の中心を重心とする正方形の頂点をなすように配置される。

図４は、本発明の一実施形態における到達時間差を説明する第１の図である。
図４に到達時間差を算出するための座標系の例を示す。まず、４つのマイクＭ１〜Ｍ４の重心を原点ｏとする。すると、空間ｗの位置Ｐ１は、原点からの距離ｒと、Ｘ軸の正方向を０度としたときの角度θを用いて、極座標（ｒ、θ）で表すことができる。ある時刻に位置Ｐ１（音源Ｐ１と記載する場合がある）にて端末装置４０が音ｓを発してから、マイクＭ１、Ｍ２、Ｍ３、Ｍ４の各々が音ｓを検出するまでの時間をそれぞれＴｍ１、Ｔｍ２、Ｔｍ３、Ｔｍ４とする。音源Ｐ１の位置が変化すると、音源Ｐ１とマイクＭ１〜Ｍ４の距離が変化するので、Ｔｍ１〜Ｔｍ４も変化する。Ｔｍ１〜Ｔｍ４は、ｒおよびθの関数として表すことができる。図５にこの様子を示す。

図５は、本発明の一実施形態における到達時間差を説明する第２の図である。
図５は、距離ｒ別に角度θとＴｍ１〜Ｔｍ４との関係を示したグラフである。例えば、左上端のグラフは、マイクＭ１〜Ｍ４の重心からの距離ｒが０．５ｍの場合の角度θとＴｍ１〜Ｔｍ４の関係を示したグラフである。到達時間Ｔｍ１〜Ｔｍ４は、角度θの変化に伴うマイクＭ１〜Ｍ４と音源Ｐ１の距離の変化に応じて図示するように変化する。他のグラフについても同様である。例えば、右下端のグラフは、マイクＭ１〜Ｍ４の重心からの距離ｒが３ｍの場合の角度θとＴｍ１〜Ｔｍ４の関係を示したグラフである。左上端のグラフ（距離ｒ＝０．５ｍ）と比べ、距離が離れた分だけ、全体的に到達時間が長くなっている。

また、マイクＭ１を基準のマイクとし、マイクＭ１が音ｃを検出する時刻を基準（時刻＝０）としたときに、他のマイクＭ２〜Ｍ４の各々が音ｃを検出する時刻（到達時間差）は以下である。つまり、マイクＭ２の到達時間差（ΔＴ２１）はＴｍ２−Ｔｍ１、マイクＭ３の到達時間差（ΔＴ３１）はＴｍ３−Ｔｍ１、マイクＭ４の到達時間差（ΔＴ４１）はＴｍ４−Ｔｍ１である。音源Ｐ１の位置が変化すると、音源Ｐ１とマイクＭ１〜Ｍ４の各々の距離が変化するので、ΔＴ２１、ΔＴ３１、ΔＴ４１も変化する。つまり、ΔＴ２１、ΔＴ３１、ΔＴ４１は、ｒおよびθの関数として表すことができる。図６にこの様子を示す。

図６は、本発明の一実施形態における到達時間差を説明する第３の図である。
図６は、角度θと到達時間差との関係を距離ｒ毎に示したグラフを、到達時間差ΔＴ２１、ΔＴ３１、ΔＴ４１別に示したものである。例えば、左端のグラフは、到達時間差ΔＴ２１について、距離ｒ＝０．５ｍ、１ｍ、２ｍの場合の到達時間差と角度θの関係を示している。音ｓの音源Ｐ１の位置とマイクＭ１およびマイクＭ２の位置関係の変化に応じて、到達時間は図示するように変化する。他のグラフについても同様である。

図５、図６で例示したグラフは、例えば、音速と距離と到達時間の関係（ｄ＝ｃｔ、ｄは距離、ｃは音速、ｔは到達時間）から計算により作成されたものである。なお、音速は温度によって変化することから、図５、図６のグラフを所定の温度ごとに用意しておいてもよい。図６で例示した到達時間差ΔＴ２１、ΔＴ３１、ΔＴ４１のグラフが予め用意されており、実際に空間ｗのある位置Ｐ１にて音ｓが発せられると、その音ｓについて計測したΔＴ２１ｍ、ΔＴ３１ｍ、ΔＴ４１ｍに基づいて、音源Ｐ１の位置が推定できる。例えば、到達時間差ΔＴ２１ｍが２０００μｓであれば、音源Ｐ１の位置情報は、左端のグラフより、ｒ＝０．５ｍ、θ＝２４０度で示される位置であると推定できる。また、到達時間差ΔＴ２１ｍが０μｓであれば、音源Ｐ１の位置は、ｒ＝０．５ｍ，１ｍ，２ｍの何れか、θ＝０度または１８０度で示される位置である。候補が複数あるので、残りの２つのグラフを、ΔＴ３１ｍ、ΔＴ４１ｍに基づいて逆引きして、音源Ｐ１の位置を推定することができる。つまり、ΔＴ２１ｍ、ΔＴ３１ｍ、ΔＴ４１ｍについて、対応するグラフを逆引きして得られた位置情報（座標情報）の候補うち一致するものが音源Ｐ１の推定位置である。
また、到達時間差Δ２１等のグラフが温度ごとに用意されているならば、空間ｗの室温に応じた到達時間差Δ２１等のグラフを参照して逆引きを行うことで、より正確な音源の位置を推定することができる。

なお、説明の便宜上、音源Ｐ１の位置を平面の極座標（図４）で表すこととして説明を行ったが、マイクの数が４つ以上で、到達時間差のグラフが３つ以上用意できれば、音源Ｐ１の３次元における位置を推定することができる。また、例えば、オフィス等の、着席して業務を行う空間ｗを想定し、空調対象の高さ方向を固定することができるならば、２次元の位置を推定するだけで良い。その場合は、マイクの数は最低限３つ存在すればよく、到達時間差のグラフは２つ以上用意できれば、音源Ｐ１の位置推定は可能である。

図６で例示したグラフをデータテーブルの形式で表したものがルックアップテーブルである。図７に音源位置を２次元座標で表したときのＬＵＴ３７の一例を示す。
図７は、本発明の一実施形態におけるＬＵＴの一例を示す図である。
図６で例示したものは、音源Ｐ１の位置を極座標で表したグラフであるのに対し、図７のＬＵＴ３７は、音源Ｐ１の位置を直交座標で表したときのＸ，Ｙ座標と到達時間差Δ２１との関係を規定したテーブルである。ＬＵＴ３７の各列は図４で例示したマイクＭ１〜Ｍ４の重心を原点とする直交座標におけるＸ座標の値が対応し、各行はＹ座標の値が対応する。テーブルの各セルには、対応する行と列が示す座標位置に音源Ｐ１が存在するときの到達時間Δｔ２１の値が記載されている。このテーブルは次のような手順で作成することができる。（１）例えば、Ｘ＝［−１０、１０］、Ｙ＝［−１０、１０］［ｍ］の範囲で空間を格子状に分割する。格子の間隔は、空間ｗの大きさに応じて任意に決定する。（２）次に各格子点（ｘ、ｙ）から音を発したときの各マイクＭｉ（ｉ＝１〜４）への到達時間ｔｉを計算する。ここで、格子点と各マイクＭｉの距離をｄ_ｉ、音速をｃとすると、到達時間ｔ_ｉ＝ｄ_ｉ／ｃである。（３）次にマイクＭ１を基準としたΔＴ２１、ΔＴ３１、ΔＴ４１を計算し、ΔＴ２１（ｘ、ｙ）、ΔＴ３１（ｘ、ｙ）、ΔＴ４１（ｘ、ｙ）のマップを作成する。図７のＬＵＴはこのようにして作成されたテーブルの一例である。

図６のグラフと同様、Ｚ軸（高さ方向）の項目をテーブルに含め、音源位置の３次元座標と到達時間差Δｔ２１との関係を定義したＬＵＴ３７を用いることができる。図７のテーブルは記憶部３６に記録されている。図示したものと同様、到達時間差ΔＴ３１、到達時間差ΔＴ４１に対するＬＵＴ３７も記憶部３６に記録されている。また、温度毎にそれぞれの温度における到達時間差ΔＴ２１〜ΔＴ４１についてのＬＵＴ３７を記憶部３６に記録しておいてもよい。湿度、気圧についても同様である。

例えば、計測した到達時間差Δｔ２１ｍの値が２１３１．６の場合、位置推定部３３は、ＬＵＴ３７を参照して、音源の位置を（１０、１０）と推定する。位置推定部３３は、同様に計測した到達時間差ΔＴ３１ｍの値と到達時間差ΔＴ３１に関するＬＵＴ３７、計測した到達時間差ΔＴ４１ｍの値と到達時間差ΔＴ４１に関するＬＵＴ３７を用いて、音源Ｐ１の座標情報を算出する。ＬＵＴ３７に登録されていない到達時間差については、線形補間、補外などにより、対応する座標情報を算出する。位置推定部３３は、各マイク間の到達時間差に対応するＬＵＴ３７を用いて算出した複数の座標情報が一致すれば、その座標が示す位置を音源Ｐ１の位置として推定する。あるいは、位置推定部３３は、ＬＵＴ３７に基づいて得られた複数の座標情報の平均を算出して、その座標が示す位置を音ｓの音源Ｐ１の位置としてもよい。

また、例えば、位置推定部３３は、次式（１）によって定義される到達時間差の実測値とＬＵＴ３７に登録された到達時間差の誤差量ε_ＲＭＳが最小となるような座標（ｘ、ｙ）を算出し、その座標（ｘ、ｙ）を音源Ｐ１の位置情報として推定してもよい。

ここで、ΔＴ２１ｍ、ΔＴ３１ｍ、ΔＴ４１ｍは到達時間差の実測値、ΔＴ２１（ｘ、ｙ）、ΔＴ３１（ｘ、ｙ）、ΔＴ４１（ｘ、ｙ）は、ＬＵＴ３７に登録された到達時間差の値または補間、補外により算出された到達時間差の値である。式（１）の最小値は、局所探索法、最急降下法、最小二乗法などにより効率的に算出することができる。
なお、ＬＵＴ３７では、極座標と到達時間差との対応関係が規定されていてもよい。

次に本実施形態の位置推定方法により音源Ｐ１の位置を推定した結果の例を示す。
図８は、本発明の一実施形態における位置推定結果の一例を示す第１の図である。
音源Ｐ１の真の位置のＸＹＺ座標を（１，３，３）とする。図には高さ方向（Ｚ軸）を省略したマイクＭ１〜Ｍ４と音源Ｐ１の位置関係が示してある。本実施形態の位置推定方法による音源Ｐ１の座標の推定結果は（１，３，３）である。一方、従来の４つの連立方程式の解をニュートンラプソン法によって算出する方法を用いると、座標の推定結果は（１，３，３）である。この推定結果により、本実施形態の到達時間差ΔＴ２１〜ΔＴ４１とＬＵＴ３７とを用いた位置推定を行うことで、計算負荷を低減しつつ、音源Ｐ１の真の位置を正確に推定できることが確認できた。

図９は、本発明の一実施形態における位置推定結果の一例を示す第２の図である。
音源Ｐ１の真の位置のＸＹＺ座標を（０，５，３）とする。本実施形態の位置推定方法による音源Ｐ１の座標の推定結果は（０，５，３）である。一方、従来の４つの連立方程式の解をニュートンラプソン法によって算出する方法を用いると、座標の推定結果は（０，−６．９，−２．６）である。従来法による解を図中に丸印で示す。この例のように２つのマイク（Ｍ１およびＭ３、Ｍ２およびＭ４）から等距離にあるような位置に音源Ｐ１がある場合、音源Ｐ１の位置は特異点となり、従来法を用いると正しく位置推定を行うことができない。これに対し、本実施形態の位置推定方法によれば、２つのマイクから等距離にあるような位置に音源Ｐ１がある場合であっても、位置を正確に推定できることが確認できた。

次に本実施形態の位置推定方法の流れについて説明する。
図１０は、本発明の一実施形態における位置推定方法の一例を示すフローチャートである。
前提として、空間ｗの天井にはマイクＭ１〜Ｍ４が設けられている。また、記憶部３６には、到達時間差ΔＴ２１、ΔＴ３１、ΔＴ４１と音源Ｐ１の位置（座標情報）との対応テーブル、つまり図７で例示したＬＵＴ３７が記録されている。
ユーザは、空間ｗにおける所望の位置（空調対象とするエリア）に移動し、端末装置４０を操作して、所定の周波数の音ｓを出力させる。すると、マイクＭ１〜Ｍ４が音ｓを検出する（ステップＳ１１）。音情報取得部３１は、マイクＭ１〜Ｍ４の各々から音信号を取得する。音情報取得部３１は、マイクＭ１〜Ｍ４の各々から音信号を取得した時刻と各マイクの識別情報を対応付けて記憶部３６に記録する。次に到達時間差算出部３２が、マイクＭ１による音信号の取得時刻を基準とする他のマイクＭ２、Ｍ３、Ｍ４の到達時間差ΔＴ２１ｍ、ΔＴ３１ｍ、ΔＴ４１ｍを算出する（ステップＳ１２）。次に到達時間差算出部３２は、算出した実際の到達時間差ΔＴ２１ｍ、ΔＴ３１ｍ、ΔＴ４１ｍを位置推定部３３へ出力する。

次に到達時間差算出部３２は、図７で例示したＬＵＴ３７を参照して（ステップＳ１３）、各到達時間差ΔＴ２１ｍ、ΔＴ３１ｍ、ΔＴ４１ｍと対応する（最も近い）音源Ｐ１の座標情報の候補値を取得する。このとき、ＬＵＴ３７に到達時間差ΔＴ２１ｍ等と等しい値が登録されていなければ、補間、補外などを行って音源Ｐ１の座標情報の候補値を取得する。次に位置推定部３３は、ＬＵＴ３７から取得した音源Ｐ１の座標情報の候補値に基づいて、音源Ｐ１の位置を推定する（ステップＳ１４）。例えば、到達時間差算出部３２は、座標情報の候補値を上記した式（１）に適用して、さらに誤差量ε_ＲＭＳが最小となるような座標（ｘ、ｙ、ｚ）を算出し、その座標（ｘ、ｙ、ｚ）が示す位置を音源Ｐ１の位置として推定する。あるいは、到達時間差算出部３２は、到達時間差ΔＴ２１ｍに基づいてＬＵＴ３７から取得した座標情報１と、到達時間差ΔＴ３１ｍに基づいてＬＵＴ３７から取得した座標情報２と、到達時間差ΔＴ４１ｍに基づいてＬＵＴ３７から取得した座標情報３の平均を算出し、その値を音源Ｐ１の位置情報として推定してもよい。位置推定部３３は、推定した位置情報（座標情報）を記憶部３６に記録する。

ＬＵＴ３７には、温度および湿度および気圧のうち少なくとも一つごとに到達時間差ΔＴ２１〜ΔＴ４１と音源位置の座標情報とが定められていてもよい。その場合、位置推定部３３は、温度および湿度および気圧のうち少なくとも一つと、実際の到達時間差ΔＴ２１ｍ等に基づいて、ＬＵＴ３７からその温度等に応じたテーブルを選択する。位置推定部３３は、選択したテーブルを参照して音源Ｐ１の座標情報の候補値を取得する。例えば、センサ情報取得部３４は、温度センサから空間ｗの室温を取得し、位置推定部３３に出力する。位置推定部３３は、その温度と到達時間差ΔＴ２１ｍに基づいてＬＵＴ３７のうち最も近い温度について用意されたテーブルを参照して座標情報を取得する。到達時間差ΔＴ３１ｍ、ΔＴ４１ｍについても同様である。そして、位置推定部３３は、式（１）を用いて音源Ｐ１の位置を推定する。

また、上記のように音速ｃ［ｍ／ｓ］、音源Ｐ１とマイクＭ１の間の距離をｄ［ｍ］、到達時間ｔ［ｓ］の間には次の関係がある。
ｄ＝ｃ（Ｔ）×ｔ
ここで、音速ｃ（Ｔ）は温度Ｔの関数で、経験的に以下が成り立つ。
ｃ（Ｔ）＝３３１．５＋０．６×Ｔ［ｍ／ｓ］ （Ｔ［ｄｅｇＣ］はセ氏温度）
従って、様々な温度ごとにＬＵＴ３７を用意しておかなくても、例えば、ある温度（２５ｄｅｇＣ）におけるＬＵＴ３７を用意しておけば、センサ情報取得部３４が取得した温度を用いて、上記した（１）〜（３）の手順によって、温度の影響分を加算して補正した音速に基づく、補正後の到達時間差を算出し、任意の温度に応じたテーブルを作成することができる。

また、温度以外にも音速に影響する環境条件として、気圧や湿度に対するシミュレーションなどを行って、気圧や湿度に応じた音速を算出できるようにし、気圧や湿度に応じたテーブルを作製できるようにしてもよい。

本実施形態によれば、マイクＭ１〜Ｍ４に音源Ｐ１から発された音が到着する時刻の差と、空間ｗ内の様々な位置から音を発したときにマイクＭ１〜Ｍ４にその音が到着する時刻の差と音を発した位置とを対応付けたＬＵＴ３７とに基づいて音源Ｐ１の正確な位置を推定することができる。また、本実施形態の到着時間差とＬＵＴ３７とに基づく位置推定方法は、計算負荷が比較的低く、局所解に陥る等の問題もなく、正確な位置を推定できることから、空調機１０が備えるマイコン等の位置推定装置３０にも好適である。
本実施形態の位置推定装置３０を備えることにより、制御装置２０は、ユーザが所望するエリアを対象とする空調制御を行うことができる。

また、端末装置４０をユーザが所持するようにし、所定の時間間隔で音ｓを出力するように設定しておくことで、ユーザが空間ｗ内を移動しても、制御装置２０は、その移動に追従して空調対象とするエリアを変更することができるのでユーザの満足度が向上する。また、音源Ｐ１の３次元座標の推定を行うと、空間ｗが複数のフロアに跨るような場合でも、例えば、ユーザが１階から２階へ移動すれば、空調対象エリアを２階に自動的に変更することができ、利便性が向上する。

また、位置推定装置３０が推定した位置を蓄積するようにし、一定期間の推定結果に基づいて空調制御を行うようにしてもよい。例えば、空間ｗのユーザが、ある位置（例えば、ソファがある場所）にいる時間が長ければ、制御装置２０は、端末装置４０から音ｓの出力がなくても当該位置を対象とした空調制御を行うようにしてもよい。また例えば、端末装置４０から断続的に音を出力するような設定にしておき、ユーザの移動に伴う位置情報を蓄積することによってユーザの動線を解析し、空間ｗの出入り口の位置を推定する。そして、出入り口の推定結果に基づいて、制御装置２０は、推定した出入り口の方向には、送風を行わないような制御を行ってもよい。これにより省エネの効果が期待できる。

なお、上記の実施形態では、マイクＭ１〜Ｍ４の全てを同時に用いることとしたが、音源Ｐ１の平面での位置の推定は、マイクが３つあれば可能である。平面上の位置推定を行う場合、マイクＭ１〜Ｍ３だけを使用して位置推定を行い、何れかのマイクが故障した場合に残りのマイクを使用するようにしてもよい。あるいは、４つのマイクのうち、任意の３つのマイクを毎回選択して、選択したマイクによって位置推定を行うようにしてもよい。

または、平面での位置推定を行う場合、１回の音源Ｐ１の位置推定について、マイクＭ１、Ｍ２、Ｍ３を用いた推定、マイクＭ１、Ｍ３、Ｍ４を用いた推定、マイクＭ１、Ｍ２、Ｍ４を用いた推定、マイクＭ２、Ｍ３、Ｍ４を用いた推定の４回の位置推定処理を行って、４回の推定結果の平均を、最終的な音源Ｐ１の位置として決定してもよい。

また、上記の実施形態では、音ｓを端末装置４０から出力することとしたが、ユーザの声を予め登録しておき、位置推定装置３０が、ユーザの声が発された位置を推定して、制御装置２０が、その推定位置に対して空調制御を行うようにしてもよい。

＜地図情報の作成＞
次にこれまでに説明した方法で推定した位置情報を用いて地図情報を作成し、その地図情報を空調制御に利用する方法について説明する。例えば、空調機１０の据え付け時に据え付け作業を行うスタッフが、図示しない端末装置４０ａから音ｓａを出力しつつ、空間ｗの内壁沿いに移動する。位置推定装置３０は、所定の周期でマイクＭ１〜Ｍ４が出力した音信号を取得し、その位置情報を推定する。地図情報作成部３５は、これらの位置情報が示す点を結んで、空間ｗの形状を表した地図情報を作成する。地図情報は、空間ｗの内側形状を模した２次元または３次元の図である。図１１に地図情報の一例を示す。

図１１は、本発明の一実施形態における地図情報の一例を示す図である。
地図情報作成部３５は、記憶部３６に記録された座標情報を用いて、図１１で例示する地図情報Ｇを作成する。さらに地図情報作成部３５は、例えば、地図情報Ｇの内側の領域を、空調機１０が空調制御を区別して行うことが可能な大きさの単位領域で分割してもよい。図１１の各セル（Ｅ１１など）が分割後の単位領域である。地図情報作成部３５は、地図情報Ｇを記憶部３６に記録する。

地図情報Ｇは、例えば、制御装置２０によって、空間ｗのユーザが所有する端末装置４０へ送信され、端末装置４０の表示画面に表示される。すると、ユーザは、端末装置４０の表示画面から所望のエリアを選択する（図１１のＰ２）。すると、端末装置４０は、Ｐ２が含まれるエリアの識別情報Ｅ３３（情報ｋ）を、制御装置２０へ送信する。すると、制御装置２０は、Ｅ３３が示すエリアを対象とする空調制御を行う。

図１２は、本発明の一実施形態における地図情報の作成処理の一例を示すフローチャートである。
前提として、据え付けスタッフは、空間ｗを壁沿いに一方方向に移動するとする。まず、空調機１０の据え付け作業を行うスタッフが、空間ｗの内の所定の位置に移動して、端末装置４０ａを操作し、端末装置４０ａから音ｓａを出力させる（ステップＳ２１）。マイクＭ１〜Ｍ４は、音ｓａを検出し（ステップＳ２２）、音信号を位置推定装置３０へ出力する。位置推定装置３０では、音情報取得部３１が音信号を取得する。また、位置推定部３３が、図１０のフローチャートで説明した処理によって、端末装置４０ａが音ｓａを発した音源の位置を推定する（ステップＳ２３）。位置推定部３３は、推定した音源の位置情報を、音ｓａを検出した時刻と対応付けて記憶部３６に記録する（ステップＳ２４）。次に地図情報作成部３５が、全ての位置について位置情報を記録済みかどうかを判定する（ステップＳ２５）。例えば、据え付けスタッフが、空間ｗの計測が終了したことを示す操作を端末装置４０ａに対して行うと、端末装置４０ａが測定完了信号を位置推定装置３０へ送信する。地図情報作成部３５は、測定完了信号を取得すると、全ての位置について位置情報を記録済みであると判定する。全ての位置情報が記録済みではない場合（ステップＳ２５；Ｎｏ）、ステップＳ２１からの処理を繰り返す。全ての位置情報が記録済みの場合（ステップＳ２５；Ｙｅｓ）、地図情報作成部３５は、記憶部３６に記録された位置情報を、位置情報と対応して記録された時刻順に結んで空間ｗの形状を表す地図情報を作成する（ステップＳ２６）。

本実施形態の位置推定方法を活用すると、空間ｗの地図情報を作成することができる、この地図情報は、空調制御の対象とするエリアの選択に用いることができる。なお、地図情報作成部３５が作成する地図情報は、空間ｗの立体形状を示した３次元の地図情報であってもよい。

図１３は、本発明の一実施形態における位置推定装置のハードウェア構成の一例を示す図である。
コンピュータ９００は、ＣＰＵ９０１、主記憶装置９０２、補助記憶装置９０３、インタフェース９０４を備える。
上述の位置推定装置３０は、コンピュータ９００を備える。そして、上述した各処理部の動作は、プログラムの形式で補助記憶装置９０３に記憶されている。ＣＰＵ９０１は、プログラムを補助記憶装置９０３から読み出して主記憶装置９０２に展開し、当該プログラムに従って上記処理を実行する。また、ＣＰＵ９０１は、プログラムに従って、上述した記憶部３６に対応する記憶領域を主記憶装置９０２に確保する。

少なくとも１つの実施形態において、補助記憶装置９０３は、一時的でない有形の媒体の一例である。一時的でない有形の媒体の他の例としては、インタフェース９０４を介して接続される磁気ディスク、光磁気ディスク、光ディスク、半導体メモリ等が挙げられる。また、このプログラムが通信回線によってコンピュータ９００に配信される場合、配信を受けたコンピュータ９００が当該プログラムを主記憶装置９０２に展開し、上記処理を実行してもよい。

上述した位置推定装置３０における各処理の過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムを位置推定装置３０のコンピュータが読み出して実行することによって、上記処理が行われる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしてもよい。

また、上記プログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であってもよい。
位置推定装置３０は、１台のコンピュータで構成されていても良いし、通信可能に接続された複数のコンピュータで構成されていてもよい。

その他、本発明の趣旨を逸脱しない範囲で、上記した実施の形態における構成要素を周知の構成要素に置き換えることは適宜可能である。また、この発明の技術範囲は上記の実施形態に限られるものではなく、本発明の趣旨を逸脱しない範囲において種々の変更を加えることが可能である。音情報取得部３１は時刻情報取得部の一例である。音源Ｐ１の位置は発信位置の一例である。

１０・・・空調機
２０・・・制御装置
３０・・・位置推定装置
３１・・・音情報取得部
３２・・・到達時間差算出部
３３・・・位置推定部
３４・・・センサ情報取得部
３５・・・地図情報作成部
３６・・・記憶部
３７・・・ＬＵＴ（ルックアップテーブル）
４０・・・端末装置
Ｍ１、Ｍ２、Ｍ３、Ｍ４・・・マイク
ｗ・・・空間
ｓ・・・音
９００・・・コンピュータ
９０１・・・ＣＰＵ
９０２・・・主記憶装置
９０３・・・補助記憶装置
９０４・・・インタフェース

Claims (10)

1. 空間に設けられた複数のマイクが、前記空間において発された音を検出した各々の時刻を取得する時刻情報取得部と、
   取得した各々の前記時刻に基づいて、複数の前記マイクのうち基準となる１つの前記マイクが前記音を検出した時刻と他の各々のマイクが検出した時刻との差である到達時間差を算出する到達時間差算出部と、
   予め定められた前記到達時間差および所定の音源位置の関係を定めたルックアップテーブルと、実際に前記到達時間差算出部が算出した前記到達時間差とに基づいて、前記音の発信位置を推定する位置推定部と、
   を備える位置推定装置。
2. 前記位置推定部は、
   前記到達時間差算出部が算出した前記到達時間差と、前記ルックアップテーブルに基づく前記到達時間差との差が最小となるときの前記ルックアップテーブルに基づく前記音源位置を、前記音の発信位置として推定する、
   請求項１に記載の位置推定装置。
3. 前記ルックアップテーブルには、複数の前記マイクが配置された位置の重心を原点とする座標系で表された前記音源位置と前記到達時間差との関係が記録されている、
   請求項１または請求項２に記載の位置推定装置。
4. 前記ルックアップテーブルには、前記マイクが配置された空間の環境条件に応じた前記到達時間差および前記音源位置の関係が記録されている、
   請求項１から請求項３の何れか１項に記載の位置推定装置。
5. 前記空間の環境条件は、温度および湿度および気圧のうち少なくとも一つである、
   請求項４に記載の位置推定装置。
6. 前記ルックアップテーブルには、所定の温度における前記到達時間差および前記音源位置の関係が記録され、
   前記位置推定部は、前記ルックアップテーブルに基づく前記到達時間差を実際の温度に基づいて補正し、補正後の前記ルックアップテーブルと前記到達時間差算出部が算出した前記到達時間差とに基づいて、前記音の発信位置を推定する、
   請求項１から請求項３の何れか１項に記載の位置推定装置。
7. 前記空間の内部を移動しながら発された前記音について前記位置推定部が推定した前記音の発信位置に基づいて、前記空間の形状を示す地図情報を作成する地図情報作成部、
   をさらに備える請求項１から請求項６の何れか１項に記載の位置推定装置。
8. 請求項１から請求項７の何れか１項に記載の位置推定装置を備え、
   前記位置推定装置が推定した前記音の発信位置を対象として空調を行う、空調システム。
9. 空間に設けられた複数のマイクが、前記空間において発された音を検出した各々の時刻を取得するステップと、
   取得した各々の前記時刻に基づいて、複数の前記マイクのうち基準となる１つの前記マイクが前記音を検出した時刻と他の各々のマイクが検出した時刻との差である到達時間差を算出するステップと、
   予め定められた前記到達時間差および所定の音源位置の関係を定めたルックアップテーブルと、実際に前記到達時間差を算出するステップで算出した前記到達時間差とに基づいて、前記音の発信位置を推定するステップと、
   を有する位置推定方法。
10. コンピュータを、
    空間に設けられた複数のマイクが、前記空間において発された音を検出した各々の時刻を取得する手段、
    取得した各々の前記時刻に基づいて、複数の前記マイクのうち基準となる１つの前記マイクが前記音を検出した時刻と他の各々のマイクが検出した時刻との差である到達時間差を算出する手段、
    予め定められた前記到達時間差および所定の音源位置の関係を定めたルックアップテーブルと、実際に前記到達時間差を算出する手段が算出した前記到達時間差とに基づいて、前記音の発信位置を推定する手段、
    として機能させるためのプログラム。

Images