音波パルスを用いて探索側端末より被探索側端末の位置関係を推定する際の直接波抽出方法であって、前記位置関係の推定を開始する旨の情報と、前記音波パルスの所定の送信回数及び送信間隔の情報を前記探索端末より前記被探索端末へ無線通信で通知した後、当該通知を前記被探索側端末が受信した旨を前記探索側端末に返信する推定開始確認段階と、前記推定開始確認段階の後、前記被探索側端末が録音を開始する録音開始段階と、前記推定開始確認段階の後、前記探索側端末が前記所定の送信回数及び送信間隔を設けた音波パルスを送信した後、前記被探索側端末に音波パルス送信終了を無線通信で通知するパルス送信段階と、前記音波パルス送信終了の通知を受信して前記被探索側端末が前記録音を終了し、該録音に含まれる音波パルスのうち、受信間隔が前記所定の送信間隔と所定基準内で等しいと判定され、且つ前記所定の送信回数に等しい個数の音波パルスを直接波として抽出し、該直接波の各々の受信時刻を前記探索側端末に無線通信で通知する直接波抽出段階と、前記通知された直接波の受信時刻に基づいて前記探索側端末が前記被探索側端末の位置関係を推定する推定段階とを備えることを特徴とする直接波抽出方法。前記録音開始段階に先だって、前記推定に用いる音波パルスの各々に共通な所定の波形情報を予め定める段階をさらに備え、前記パルス送信段階は前記所定の波形情報で形成される音波パルスを送信し、前記直接波抽出段階は前記録音と前記所定の波形情報との相関関数の値が所定基準を満たすピーク箇所を前記含まれる音波パルスとすることを特徴とする請求項1に記載の直接波抽出方法。前記直接波抽出段階にてさらに、前記受信間隔と前記所定の送信間隔との差が以下の式で与えられる許容値以内に収まる音波パルスを直接波として抽出することを特徴とする請求項1または2に記載の直接波抽出方法。=×÷前記探索側端末及び前記被探索側端末が通話機能を備えた携帯情報端末であって、前記被探索側端末が前記録音を行うための受波器が1個のみであり且つ前記通話機能も担い、前記探索側端末が前記音波パルスを送信する送波器が前記通話機能も担うことを特徴とする請求項1ないし3のいずれかに記載の直接波抽出方法。請求項1ないし4のいずれかに記載の直接波抽出方法で用いられる探索側端末。請求項1ないし4のいずれかに記載の直接波抽出方法で用いられる被探索側端末。上記のような、送信端末と受信端末との各種の位置関係の推定を、携帯端末のような利用可能なハードウェア及び該ハードウェアの性能に制限のある端末にて送受信される音波を用いて、負荷の軽い処理で実現できれば、様々な有用な利用が可能となる。そしてここでは、音波の中から直接波のみを抽出することが必須である。しかしながら、特許文献1に開示された技術による直接波の抽出では、受信側に複数のマイクを必要とし、一般に単数のマイクのみを備える携帯端末に適用しようとすると、不足マイクの追加によって追加ハードウェアが必要となる。また仮に追加したとしても、マイク数分の音波データ列の処理が必要となり、処理が重くなる。本発明の目的は、上記従来技術の課題を解決し、通常のハードウェアを備える端末において軽い処理により実現可能な、パルス状の音波を用いて探索側端末より被探索側端末の位置関係を推定するに際して直接波のみを抽出する直接波抽出方法並びに当該方法にて用いられる探索側端末及び被探索側端末を提供することにある。上記の目的を達成するため、本発明は、音波パルスを用いて探索側端末より被探索側端末の位置関係を推定する際の直接波抽出方法であって、前記位置関係の推定を開始する旨の情報と、前記音波パルスの所定の送信回数及び送信間隔の情報を前記探索端末より前記被探索端末へ無線通信で通知した後、当該通知を前記被探索側端末が受信した旨を前記探索側端末に返信する推定開始確認段階と、前記推定開始確認段階の後、前記被探索側端末が録音を開始する録音開始段階と、前記推定開始確認段階の後、前記探索側端末が前記所定の送信回数及び送信間隔を設けた音波パルスを送信した後、前記被探索側端末に音波パルス送信終了を無線通信で通知するパルス送信段階と、前記音波パルス送信終了の通知を受信して前記被探索側端末が前記録音を終了し、該録音に含まれる音波パルスのうち、受信間隔が前記所定の送信間隔と所定基準内で等しいと判定され、且つ前記所定の送信回数に等しい個数の音波パルスを直接波として抽出し、該直接波の各々の受信時刻を前記探索側端末に無線通信で通知する直接波抽出段階と、前記通知された直接波の受信時刻に基づいて前記探索側端末が前記被探索側端末の位置関係を推定する推定段階とを備えることを特徴とする。前記特徴によれば、受信録音に含まれる音波パルスの中から、送信した音波パルスの送信回数及び送信間隔と合致するものを直接波によるものと判定して抽出するので、探索側端末及び被探索側端末がハードウェア等に制限のある端末であっても正確に、端末間の位置関係を推定するのに必要な直接波が抽出できる。以下、図面を参照して本発明を詳細に説明する。本発明では、予め所定の送信回数及び送信間隔を設定しておき、当該設定に従うパルス状の音波を送信端末がスピーカ等にて送信し、受信端末がマイク等にて受信し録音する。受信端末は前記設定における情報を利用して、録音の中から送信端末が送信した直接波によるパルスを選択し、当該選択結果を送信端末に伝える。送信端末は当該選択結果を用いて、受信端末との位置関係を推定する。推定される位置関係としては、端末間の距離や端末同士の相対的な向きなどがある。以下、上記のような位置関係の推定を行う側と受ける側という意味で、音波の送信端末を探索側、音波の受信端末を被探索側と呼ぶこととする。図1は本発明による直接波のパルス抽出の原理を説明する図である。前提として、探索側と被探索側で、予め用いる音波パルスの時間軸上の波形情報を共有しているものとする。は単一のパルスのみを送信し受信して検出する場合の例である。被探索側は受信録音と予め共有して保持している音波パルス波形との、に示すような相互相関関数を算出する。そして相関値が所定の閾値を超えて且つピーク値となる時刻t1を直接波受信時刻と判定する。ではそのような条件を満たす箇所は時刻t1の1箇所のみであるので、当該t1を直接波受信時刻と判定すればよい。しかし、単一パルスを用いると、反射波や干渉波あるいはその他の雑音によって、に示すような例えば時刻t20及びt21といったような、2箇所以上で所定の閾値を超えるピークが得られた場合に、いずれが直接波によるものであるかが判定できない。の例では時刻t20の方が相関値が大きいが、干渉波などの可能性もあり、相関値最大を選出して直接波の受信時刻を判定することはできない。そこで本発明では、に示すように、所定回数及び所定間隔を設けた複数の音波パルスを利用する。の例では、100msの間隔を設けた同一波形の2個の音波パルスを送信し、受信録音に対して相互相関関数を算出した例である。相関値が所定の閾値を超えるピーク箇所がt30、t31、t32及びt33の4箇所得られるが、2個のパルスを100msの間隔を設けて送信したという情報を利用して、これらの中から直接波受信時刻として、2個のパルスであって間隔が100msとなるt32及びt33を選択する。なお、における回数及び間隔を予め乱数などで設定することによって、仮に同様の波形情報を発する他の音源が存在したとしても区別することができる。また、のように回数及び間隔を設定したとしても、直接波の他に反射波が同様の回数及び間隔で現れうることも考えられる。しかしながら、探索側と被探索側とが少しでも移動すると、反射波の経路は大きく変動するため、このように反射波と直接波とが区別できなくなる可能性は低い。特に本発明では探索側及び被探索側の両者に携帯情報端末を想定しているので、このような可能性は低い。なお特に本発明の利点として、以上の図1ののような直接波抽出処理が、単一マイクによる受信録音で簡素な処理によって実施できる、という点がある。図2は本発明に係る探索側端末又は被探索側端末として機能することのできる、一般的な携帯情報端末の構成の一例を示した機能ブロック図である。携帯情報端末100は、制御部101、記憶部102、入力部103、表示部104、通信部105、マイク107及びスピーカ108を含む音声入出力部106、計時部109、付加機能部110並びにバス200を備える。制御部101は各部および各部間の制御を行い、記憶部102は当該制御部等の動作に必要な一時的又は固定的な情報を記憶する。入力部103はキーパッドやタッチパネルなどを含みユーザ入力を受け取る。表示部104はディスプレイなどを含み、各種情報をユーザに向けて表示する。通信部105はアンテナ等を含み各種の形態における無線通信を担う。音声入出力部106はマイク107とスピーカ108とを含み、携帯情報端末100が通話機能を備える場合であれば通話時に利用され、その他、本発明で実施される音声データの再生や録音などにも用いられる。計時部109は時計の役割を果たす。付加機能部110はその他各種のユーザ用途に応じた機能を提供する。バス200は各部間のデータ転送を担う。なお、付加機能部110は制御部101に含まれていてもよい。以下に説明するように、本発明の探索側端末及び被探索側端末はこのような一般的な携帯情報端末100を用いて部分的機能として、または追加的機能として実現できる。また以下の説明より明らかなように、本発明の探索側端末及び被探索側端末の両端末機能を一端末において同時に備え、探索側/被探索側の各立場となる場合毎に機能を使い分ける端末も一般的な携帯情報端末100にて実現できる。なおまた、携帯情報端末100は、同等の機能を実現できる一般の情報端末であってもよい。図3は本発明に係る探索側端末と被探索側端末との、位置関係の推定を行う際の関係を含んだ機能ブロック図である。探索側端末10は相対方向推定部11、向き取得部12、入力部13、無線通信部14、音波設定送信部15、距離推定部16及びこれら各部間のデータ転送を担うバス17を備える。被探索側端末20は無線通信部21、音波受信識別部22、入力部23及びこれら各部間のデータ転送を担うバス24を備える。当該機能ブロック図において、推定する位置関係として、探索側端末から見た被探索側端末の相対方向を推定する場合は距離推定部16を利用せずに相対方向推定11及び向き取得部12を利用する。距離を推定する場合には逆に、相対方向推定11及び向き取得部12を利用せず、距離推定部16を利用する。以下、相対方向の推定を主に説明し、適宜、距離の推定を行う場合も追加で説明することとする。図4は相対方向推定の概要を説明する図である。図4では屋外の待ち合わせ広場といったような、音波到達可能な例えば数十メートル程度の平面的な場所領域にて本発明に係る探索側の端末から被探索側の端末の方向を推定する。なお、ここで音波到達可能とは、広場にいる人などの障害物による減衰などによる制約、図2の携帯情報端末100におけるマイク107の制約、当該携帯情報端末100で利用可能なバッテリによる制約、当該携帯情報端末100を利用するユーザが不快にならないような範囲での音を用いることによる制約、などによって定まるものである。図4においてまず、探索側と被探索側は無線通信可能であり、且つ当該無線により今からからに対して相対方向推定を行うということを互いに認識している状態にあるとする。次にに示すように、例えば探索側端末のユーザが探索側端末を手に持って肘の回りに回転させる等して探索側端末をほぼ水平面上にスウィングさせる。このとき、図示しているように回転の時間経過と共に、例えば5回といったような複数回のパルス状の音波S~Sを探索側端末より送信する。そして探索側端末は各音波を発生した各時刻と各時刻における探索側端末の方向とを記録する。被探索側の端末では当該方向推定の開始にあたって、にて探索側の送信する音波パルスの所定の送信回数および間隔の情報を探索側より受け取り、適宜録音を開始し、さらに音波送信を終えた旨の通知を探索側より受け取って録音を終了する。被探索側は前記情報を用いて図1で説明した手法により、録音から対応する直接波によるパルスを識別し、各音波S~Sに対応する受信時刻を求め、探索側に当該時刻を無線送信する。このような被探索側における直接波パルスの識別の結果として、別のペアとして方向推定を行う等している通行人等の発生する音、すなわち図Aで前述の干渉波が除外される。また図4には示していないが、図Aで前述の反射波やその他の雑音も除外される。探索側は受信時刻を受け取って受信間隔を算出し、記録してある送信時刻より算出される送信間隔と比較し、相対方向を推定する。図5に当該推定に用いる各音波送信時の端末向き、送信間隔及び受信間隔の一例を示す。図5では音波S、Sを送信した間の区間では送信間隔のほうが受信間隔より長く、それ以降の区間では逆の関係となっている。図5の例における上記結果の得られるメカニズムは次の通りである。すなわち、図4に示すように、探索側と被探索側との相対方向および端末のスウィング領域の関係においてドップラー効果と同様の原理により、送信間隔より受信間隔が狭い区間音波S、S送信区間と送信間隔より受信間隔が広い区間音波S、S、S送信区間との区別が生じるためである。探索側は送信間隔と受信間隔とのこのような原理による差異を抽出して、受信間隔と送信間隔との大小関係の逆転が生ずる方向である、音波SとSとを送信した間における端末の方向を被探索側の存在する相対方向であると推定する。なおまた、図4及び図5を用いて相対方向推定の場合を説明したが、距離を推定する場合には、図4のようなスウィングを加える必要はなく、より単純に次のようにすればよい。すなわち、距離推定を開始する旨の情報を探索側から被探索側へ無線で伝えて同期を取り、被探索側は録音を開始する。ここで特に、距離推定においては単に通信のための同期を取るだけではなく、両端末において所定のプロトコルを用いるなどして時刻を正確に同期させる。探索側は所定回数及び間隔の音波パルスを送信し、当該各パルスの送信時刻を記録する。ここで距離推定においては、音波パルス送信時の探索側端末は静止させておけばよく、相対方向推定の場合のようにスウィング動作は不要である。むしろ受信間隔に誤差が生じるので、音波パルス送信時に探索側端末は極力動かさない方がよい。被探索側は録音より前記所定回数及び間隔を用いて、送信パルスの直接波に相当する受信パルスを抽出し、録音にて記録されているその各受信時刻を探索側に返信する。探索側は、送信/受信の順番が互いに対応する送信パルスと受信パルスとの記録時刻の差に音速を乗じることによって、被探索側への距離を推定する。当該乗ずる音速としては、例えば340msのように所定値を設けておけばよい。図3に戻り、全体として上記図4及び図5のような相対方向推定又は距離推定を行う探索側端末10及び被探索側端末20の各機能ブロックについて説明する。相対方向推定部11は各種情報を用いて被探索側端末20の存在している相対的な向きを推定する。距離推定部16も各種情報を用いて被探索側端末20との距離を推定する。相対方向推定部11及び距離推定部16は例えば付加機能部110に含めることができる。向き取得部12は探索側端末10の向きを計時部109が取得する時刻と連動して取得し、例えば付加機能部110に含めることができる。一実施形態では向き取得部12は地磁気センサを用いて向きを取得する。入力部13探索側端末10側はボタンなどからなり、ユーザが探索側端末10を起動して相対方向推定又は距離推定を開始させる旨や、推定の各段階の完了等をボタン押下などにより探索側端末10に通知する。入力部13には携帯端末100の入力部103を利用することができる。探索側端末10の無線通信部14および被探索側端末20の無線通信部21はBluetooth登録商標、無線LAN、3G、Zigbee登録商標等により空間中やネットワーク等を介して互いに無線通信し、音波の送信時刻送信間隔や受信時刻受信間隔といったような相対方向及び距離推定に用いるデータや、推定フローにおける各段階での端末間同期信号を端末間で送受信する。無線通信部14および無線通信部21には、携帯端末100の通信部105を利用することができる。音波設定送信部15はパルス状音波の送信回数及び間隔を設定し、スピーカなどにより当該音波を設定通りに再生して送信する。音波設定送信部15にはスピーカ108、付加機能部110及び計時部109などを利用することができる。実際に送信する音波には、例えば、一般的な携帯端末100の備えるスピーカ108により送信可能である信号長が50ms、周波数帯が26kHzであるチャープ信号を用いることで、フェーディングの影響を低減できる。スピーカマイクが対応していれば、超音波を用いることもできる。音波受信識別部22はマイクなどにより音波を受信、録音し、無線通信部21を介してあらかじめ受け取っている送信回数間隔などの情報を用いて当該録音の中から音波設定送信部15より送信された直接波を識別する。音波受信識別部22にはマイク107、付加機能部110及び計時部109などを利用することができる。入力部23被探索側端末20側は、位置関係推定の各段階における探索側端末10からの同期確認への応答などからなる被探索側端末20のユーザの入力を受け付けるが、当該入力部23を用いず被探索側端末20が自動応答してもよい。入力部23にも携帯端末100の入力部103を利用することができる。図6に、図4で説明したような相対方向及び距離の両方の場合を含む位置関係の推定を行うための、探索側端末10と被探索側端末20との間のやりとりの手順を示す。また当該手順において探索側端末10で実施される各ステップのフローチャートを図7に、被探索側端末20で実施される各ステップのフローチャートを図8に示す。図6に示す各端末自身および端末間の手順には手順~と番号を付しているが、対応するステップの説明においては図7および図8においても同様にブロック内に手順の番号i1~14を付して参照している。当該手順の説明により、図3の各機能ブロックの処理の説明も行う。手順において、探索側端末10を持っているユーザが入力部13に含まれるボタンを押下することで、端末上にインストールされている位置関係推定アプリケーションが立ち上がるなどして、位置関係推定が開始される。手順にて、被探索側端末20に対して無線通信を介して相対方向推定開始を通知し、被探索側端末20においても所定のアプリケーションを立ち上げるなど準備した後、承認の返信を受ける。手順にて、音波設定送信部15がパルス音波の送信回数及び間隔を決定する。相対方向の推定に際しては、前述のようにパルス送信間隔と受信間隔とで1ms程度の差が生ずることを積極的に利用する必要があるので、例えば次のように決定する。n個のパルス状音波Sを送信する際の個の送信間隔Itを不定期間隔とし、当該送信間隔を次式を用いて決定する。Itrand×50ms...ここで、randはm以下の自然数をランダムに発生させる変数である。例えばm3として、順次Itを定め、当該間隔の総和が所定値を超えた時点を最後の音波パルスとして個数nを定める。このようにして、図5の例のような送信回数および間隔を定めることができる。あるいは同様に、何通りかこのような送信回数及び間隔を所定のものとして求めておき、2通り以上定めている場合はその中からランダムに選ぶようにしてもよい。式のように乱数を用いずに、100ms間隔で5回といったように所定間隔及び所定回数を定めておいてもよい。距離推定の場合も、相対方向推定における上記のような手法で回数及び間隔を定めることができる。なお、各音波パルス自体の情報は、図6の各手順自体の開始前に、あらかじめ探索側端末10及び被探索側端末20にて波形情報として既知であり、上記n個のパルスは全て同じ波形のパルスであるとする。なおまた、当該パルスに用いる波形情報を複数種類設けておいて、n回同じ波形でパルスを送る際に、そのうちいずれの種類を用いるかを上記手順にて決定して通知してもよい。手順にて、当該決定した送信回数及び間隔の情報を、被探索側端末20に無線送信する。なお、上記のように所定の波形情報を有する音波パルスに対して、毎回所定の範囲内でランダムに回数及び間隔を定めることにより、直接波の抽出が効果的に行える。被探索側端末20は受信確認メッセージを返信すると共に、手順に進み音波受信識別部22がマイクを起動して音波受信すなわち録音を開始する。次の手順~は、まず相対方向推定の場合から説明する。探索側端末10は受信確認を受けると手順にてユーザが入力部13のボタンを押下して方向推定のスウィング動作を開始する旨を探索側端末10に伝え、押下と同時にスウィング動作を加える。当該スウィング動作と共に手順にて音波設定送信部15によって、所定の波形情報を所定回数及び間隔を設けて再生した音波が発信されると共に、同時に手順にて当該音波発生の各時点においてスウィング動作により変化している探索側端末10の向きを向き取得部12が取得する。当該手順後のスウィング開始と同時に手順のボタン押下により起動される手順及び手順を行うことにより、探索側端末10において各音波パルス送信時刻と当該各時刻における端末の向きとが得られる。スウィングを終えると手順にてユーザが入力部13のボタンを押下して探索側端末10にスウィング動作の終了を伝える。当該手順は送信音波を全て送信した時点又はその所定時間後に探索側端末10が自動で終了判断するようにしてもよい。以上、相対方向推定の場合について説明した手順~につき、距離推定の場合にはスウィング動作を省略して、好ましくは静止状態で同様にして実行すればよい。推定精度の面からは好ましくないが、後述のように、例えば探索側端末10を持つユーザが歩くことによって探索側端末10が移動中に実行してもよい。距離推定の場合は、上記のボタン押下が探索側端末10に所定のパルス音波を発生させる合図となる。なお、距離推定の場合には前述の手順の際に、ネットワークタイムプロトコルなどを用いて、探索側端末10及び被探索側端末20の時計を正確に合わせておく。手順にて無線通信により探索側端末10より被探索側端末20に位置関係推定の終了を通知し、被探索側端末20は確認応答した後、手順にてマイクを停止して録音音波受信を終了する。さらに被探索側端末20は手順にて、手順にて受信した送信回数及び間隔の情報、予め既知の各音波パルスの波形情報を利用することにより、自身での手順における録音データから、音波受信識別部22によって直接波パルスに相当する箇所の時刻情報を抽出する。当該抽出処理の詳細については後述する。手順にて、当該抽出された直接波受信時刻は探索側端末10に無線送信され、確認応答メッセージを返信する。手順にて探索側端末10は送信した音波パルスの間隔あるいは時刻と、無線受信した直接波受信時刻と、相対方向推定においては各送信時点における端末の向きとを用いて被探索側端末20の位置関係を推定する。当該推定は、相対方向推定であれば相対方向推定部11が、距離推定であれば距離推定部16が行う。以上の手順、、及びにおける無線通信は、探索側端末10の無線通信部14及び被探索側端末20の無線通信部21を介して、必要に応じてその間にさらに所定のネットワークを介して行われる。なお、本発明においては、上述のように探索側端末10と被探索側端末20とで位置関係推定を開始する旨を伝える同期処理を行うことによって、音波送信を行う時間区間を含みその他の余分な時間を極力含まない区間として手順~間の録音が行われるため、被探索端末20側で音波受信識別部22が録音から送信音波に対応する音を識別抽出するのが容易になるという効果がある。また以上図6の手順における探索側端末10の各ステップのフローチャートである図7を説明する。ステップS0にて相対方向推定開始するとステップS1手順にて入力部13のボタンをユーザが押下する。ステップS2手順では被探索側端末20に対して相対方向推定の開始を無線通知し、処理の同期を促す。ステップS3手順では送信回数及び間隔を音波設定送信部15が決定し、ステップS4手順では当該決定情報送信間隔などを被探索側端末20に無線通知する。ステップS5手順に進みユーザによるスウィング動作開始の合図であるボタン押下を入力部13が検知する。ステップS5と同時に、相対方向推定の場合、ユーザが探索側端末10にスウィング動作を加えると共に、ステップS6手順において設定通りの音波送信および端末向き取得を行う。相対方向推定の場合、スウィング動作を終えたユーザはステップS7手順にて入力部13のボタンを押下することによりスウィング動作終了を探索側端末10に伝える。距離推定であれば単に音波送信終了を探索側端末10に伝える。前述のように当該ステップS7はユーザによるボタン押下ではなく音波送信終了後自動で行うようにしてもよい。ステップS8手順では無線通信により位置関係推定の終了を被探索側端末20に伝える。ステップS9手順では被探索側端末20が抽出した直接波時刻情報を無線にて送信し、探索側端末10が受信する。ステップS10手順にて送信時刻、直接波受信時刻、取得向きに基づいて被探索側端末20の存在する相対方向又は該端末との距離を推定し、ステップS11にてフローは終了する。またステップS5、S7及びS9にてボタン押下が検出されなかった場合や、被探索側端末20から必要な情報が送信されてこなかった場合はフロー中にてエラーが発生したと判定し、再度ステップS1よりやり直す。なお、図6および図7で探索側端末10に対するユーザのボタン押下として説明した部分は、入力部13が例えばタッチパネル等のその他のインタフェースであってもよく、同様に各種の開始合図などの入力検知として用いられる。相対方向推定の場合であれば、インタフェースとしては手順のユーザ入力検知後に探索側端末10に余分な動作を加えることなく迅速且つスムーズにスウィング動作を開始でき、スウィング終了時も手順にて探索側端末10に余分な動作を加えることなく迅速且つスムーズに端末にユーザ入力を伝えることのできるインタフェースが好ましく、ボタン押下はその一例である。さらに以上図6の手順における被探索側端末20で実行される各ステップのフローチャートである図8を説明する。ステップS20で相対方向推定を開始すると、まずステップS21手順にて探索側端末10からの位置関係推定開始通知すなわち同期リクエストを受信し、位置関係推定に応じられるよう準備を行い応答返信する。ステップS22手順では探索側端末10より送信回数及び間隔の情報を無線にて受信し、当該情報は音波受信識別部22に後の直接波抽出処理のために渡される。その後ステップS23手順にて音波受信識別部22がマイクを起動し音波受信すなわち録音を開始する。ステップS24手順にて探索側端末10より音波送信の終了通知を受信し、応答すると共に状態を同期させる。すなわちステップS25手順に進んで音波受信識別部22において起動したマイクを停止し、録音を終了する。ステップS26手順に進み、音波受信識別部22はステップS22にて受信した情報を用いて録音データより直接波を抽出し、各直接波の受信時刻を算出する。ステップS27手順に進み、当該算出結果を探索側端末10に通知するとステップS28にてフローは終了する。また、当該フローにおいてステップS22、S24などにて、探索側端末10より推定フローが順調に進んでいる旨の情報を所定時間待っても受信できない場合にはエラーが発生しているので、再度ステップS21に戻る。なお、被探索側端末20における、探索側端末10からの各種同期要求に対する応答の処理ステップS21手順のACK、S22手順のACK及びS24手順のACKは、全て自動で行うようにしてもよい。また、当該位置関係推定の利用用途に応じて、最初のステップS21を被探索側端末20を持っているユーザによる承認の入力によって行うようにしてもよい。次に、音波受信識別部22において直接波の抽出を行い、各直接波の受信時刻が算出できている前提で、相対方向推定部11が探索側端末10からみて被探索側端末20が存在する相対方向を推定する詳細につき説明する。なお、同様の前提の元で距離推定部16が両端末間の距離を推定する場合には、あるパルスの送信時刻と対応する直接波受信時刻との差を求め、音速を乗ずればよい。図4で概略的に説明したように、音波設定送信部15が設定した音波を送信している間、同時に探索側端末10はユーザが当該端末を手に持つなどして概ね水平に且つ当該ユーザの腕の長さなどに対応する概ね一定の半径の円状に、当該ユーザの位置する静止位置を中心にスウィング動作を受ける。そして当該スウィング動作中の各パルス音波送信時における探索側端末10の向きが向き取得部12により時刻と共に取得される。なお、当該スウィング動作が所定の半径を伴い端末の向き変化と共に位置移動があることは本発明が波源の移動によるドップラー効果に対応する原理を利用していることから必須である。たとえば端末を固定位置で回転させるだけでは位置移動がほとんどないので相対方向推定は困難となる。当該スウィング動作に伴う音波送信終了後、被探索側端末20において上述の音波受信識別部22の処理によって直接波受信時刻が算出され、無線通信を介して相対方向推定部11に渡される。相対方向推定部11は当該受信時刻より直接波受信間隔を求め、予め設定し音波受信識別部22より送信した音波の送信間隔、および各パルス送信時刻に取得された向きを用いて、ドップラー原理に対応するメカニズムを利用して被探索側端末20の相対方向を推定する。すなわち、スウィング中に探索側が被探索側に近づいている区間に送信した音波の受信間隔は送信間隔より狭くなり、反対に遠ざかっている区間に送信した音波の受信間隔は送信間隔より広くなる。これより、両区間の境界となる音波を特定し、その音波を送信した際の端末向きを相対方向と推定できる。具体的には次のように推定する。n個の音波Sを送信した際の送信間隔をIt、受信間隔をIr、端末向きをDとし、相対方向RDは次式で求められる。RDDIrかつIt≦Irを満たす...ここで、被探索側端末20が探索側端末10をスウィングした方向の後ろ側に位置する場合、式の条件を満たす端末向きは得られない。その場合は次式により、受信間隔が送信間隔より広くなる区間と、受信間隔が送信間隔より広く区間の境界となる音波を特定し、その音波を送信した際の端末向きの真後ろを相対方向RDと推定する。RDDnを想定している。当該抽出処理のフローを図9に示す。当該抽出処理は、候補箇所m個から直接波受信箇所n個を見つけ出すために、m箇所からn箇所を選ぶ全ての組み合わせmCn通りの場合につき、受信間隔が対応する箇所における所定の送信間隔とほぼ等しいかを逐一確かめ、誤差が最も小さいものを選ぶという処理である。以下、開始ステップS30以降の図9の各ステップを説明する。m個の音波信号から任意のn個の音波信号の組み合わせを求める。当該組み合わせの各々を識別するカウント変数をkとし、k1から開始して全ての組み合わせkにつき以下ステップS32~S37を実施する。組み合わせkの音波信号n個の受信時刻rkを求める。組み合わせkのn個の音波信号の各受信間隔rintervalkを次式で計算する。rintervalkrkrki番目の送信間隔tintervalと、対応するi番目の受信間隔rintervalkとの絶対誤差errorkを次式で計算する。なお、当該各送信間隔tintervalは前述の通り、直接波送信回数nと共に図6の手順にて被探索端末10へ通知されている。errorkrintervalktinterval絶対誤差errorkが1つでも許容値thresholdを超過した場合、絶対誤差の累積値errorsumk∞とし、へ移動する。当該許容値については後述するが、許容値を超過する場合、当該組み合わせkは直接波の受信時刻を表していないとで自動判定するために累積値を上記のように設定する。絶対誤差errorkが全て許容値threshold以内に収まっていれば、へ進み累積値errorsumkを計算する。絶対誤差の累積値errorsumkを次式で計算する。errorsumkerrorkerrork...errork全ての組み合わせkにつき以上の処理を行ったか確認する。全て行い、kmCnとなる場合、へ移動する。そうでない場合、kk1とし、未処理の次の組み合わせに対して処理を継続するため、に戻る。絶対誤差の累積値errorsumkが最小となる組み合わせkを求め、当該組み合わせkにおける音波信号n個が直接波を受信したものであると判定し、受信時刻を求める。図10を用いて反射波を除去し直接波を抽出する、図9のフローチャートによる処理適用例を説明する。図10では送信回数n2、送信間隔interval100msとしており、録音に対する相関関数の算出を行ってパルス箇所を求めた結果、4個の音波信号s、s、s、sを受信した場合を、縦軸を録音音波における受信電力として示している。以下、上記手順を実行した例の一部を記載する。6通りの組み合わせS、S、SSSSを求める。組み合わせ1Sについて以下を実行する。rinterval150error150thresholdvtinterval100などであるので、や以外の斜め方向における送信間隔と受信間隔との差を考慮する際も、AB方向成分のみの効果を考えればよい。以上のように、距離推定の場合では又はが最大の差を生ずる方向であるので、当該vtintervalの分の許容値を設定すればよい。音速を定数として340msとする。また、探索側端末10を持つユーザの音波送信中の歩速が例えば毎秒1m以下であるとする。この場合、送信間隔tintervalmsの間に受信間隔と生じる差はms以下となるので、当該値を許容値とすればよい。より一般には次のように求める。×÷上記設定において、はにおける各i番目ごとに定めてもよいし、許容値にマージンを持たせる意味で全n1送信間隔のうちの最大のものを選んでもよい。また、は例えば秒速1mあるいは2mのように、所定値を設けておけばよい。も温度などによる変動もあるが例えば340msのように定数として設けておけばよい。また、被探索側端末20の位置Bは静止していると仮定したが、被探索をユーザが持って移動している等の場合も、上記所定速度の調整により許容値を設定することができる。次に図11の相対方向推定を行う場合の許容値設定を説明する。音波送信中、探索側端末10を持つユーザが位置Aに、被探索側端末20は位置Bに存在し、ユーザが手でスウィングさせることで探索側端末10はAを中心に半径rの円上を概ね動くものとする。この場合、各音波送信時点での速度方向はに描くように変化するが、許容値を設定するに際してはの例と同様に、マージンを見込んで最も大きな差を生ずる又はのAB方向の場合を利用する。においてもAB間の距離をLとすると、本発明の通常の用途においてLは10m程度以上、rは50cm程度、1つの送信間隔における端末の移動量Δxもrと同程度である。よってLΔxの値は20程度以上となるので、移動量のうちAB方向成分のみを考慮すればよい。よっての場合と同様に、の相対方向推定の場合には許容値を次のように設定すればよい。×÷ここで×であり、上記はラジアン単位で次のようにして求めることができる。における送信間隔tintervalの両端の時刻において探索側端末10で取得されている向きの情報を、図6の手順において追加情報として被探索側端末20へ通知し、当該両端の時刻における向きの差の角度÷送信間隔tintervalとして角速度を求める。この場合、向きの情報は全てのiにつき通知を受けるものとする。あるいは、各送信間隔tintervalにおける角速度を求めて、その中から最大のものを用いてもよいし、全体の平均の角速度を用いてもよい。の相対方向推定を行う場合でさらにA及びBの位置自体が歩きなどにより移動している場合を想定して許容値を設定する場合、間隔に最も大きな差が生ずる場合を考慮して、上記で求めた許容値とで求めた許容値とを足し合わせてもよい。なお、以上の各許容値設定に際して、のとのとを一括して、探索側端末に想定される所定速度、としてもよい。すなわち、の両者を含む許容値を次のように与えることができる。×÷以上、本発明によれば音波パルスを用いた端末間の位置推定に際して必要となる直接波の抽出を、両端末が携帯情報端末などで機能的な制限がある場合でも、確実に行うようにすることができる。また本発明によれば上記のような制限を課せられた端末間で利用可能な音声入出力デバイスを用いて、位置推定の一種として有用な端末間の相対方向の推定を行うことができる。特に当該方向推定は、特殊なハードウェア等を必要とせず、通常の携帯端末で通話用に利用される単一のマイク及び通話用に利用される通常のスピーカで実施可能である。