JP7257992B2 - 情報処理装置、情報処理方法およびプログラム - Google Patents
情報処理装置、情報処理方法およびプログラム Download PDFInfo
- Publication number
- JP7257992B2 JP7257992B2 JP2020111794A JP2020111794A JP7257992B2 JP 7257992 B2 JP7257992 B2 JP 7257992B2 JP 2020111794 A JP2020111794 A JP 2020111794A JP 2020111794 A JP2020111794 A JP 2020111794A JP 7257992 B2 JP7257992 B2 JP 7257992B2
- Authority
- JP
- Japan
- Prior art keywords
- information processing
- power distribution
- correlation coefficient
- calculation unit
- estimation
- 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.)
- Active
Links
Images
Landscapes
- Mobile Radio Communication Systems (AREA)
Description
本発明は、情報処理装置、情報処理方法およびプログラムに関する。
特許文献１に記載された受信電力推定装置では、空間的に配置された電波センサーでの受信電力と、遮蔽物を考慮した伝搬損失推定値から、無線局の送信アンテナ利得（ＥＩＲＰの指向性）を推定することが行われる。また、当該受信電力推定装置では、送信アンテナ利得の推定値と、第２の伝搬損失推定値から、空間的な受信電力を推定することが行われる（特許文献１の段落００６７～００７６、図６～１０などを参照。）。
しかしながら、上述のようなアンテナ指向性を推定する手法では、電波センサーがない方向における送信アンテナ利得を線形補間で推定するが、実際に使用されているアンテナ指向性パターンが複雑である場合、電波センサーがない方向で送信アンテナ利得の推定値に誤差が生じる場合があった。このため、送信アンテナ利得の推定値の誤差により、電波センサーがない方向における場所の推定受信電力にも誤差が生じる場合があった。
また、上述のようなアンテナ指向性を推定する手法では、無線局の正確な位置に基づいて伝搬損失推定値を求めることが行われるが、例えば、推定により求められた無線局の位置が用いられるときには、当該位置に誤差があるため、伝搬損失推定値にも誤差が生じる場合があった。
また、上述のようなアンテナ指向性を推定する手法では、無線局の正確な位置に基づいて伝搬損失推定値を求めることが行われるが、例えば、推定により求められた無線局の位置が用いられるときには、当該位置に誤差があるため、伝搬損失推定値にも誤差が生じる場合があった。
本発明は、このような事情を考慮してなされたもので、アンテナ指向性のパターンの方向を精度良く推定することができる情報処理装置、情報処理方法およびプログラムを提供することを課題とする。
一構成例として、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定する推定演算部を備える、情報処理装置である。
一構成例として、推定演算部が、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定する、情報処理方法である。
一構成例として、推定演算部が、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定するステップを、コンピュータに実行させるためのプログラムである。
本発明に係る情報処理装置、情報処理方法およびプログラムによれば、アンテナ指向性のパターンの方向を精度良く推定することができる。
以下、図面を参照し、本発明の実施形態について説明する。
（第１実施形態）
＜情報処理システム＞
図１は、本発明の実施形態に係る情報処理システム１の概略的な構成を示す図である。
情報処理システム１は、情報処理装置１１と、複数の電波センサー２１－１１～２１－ＭＮと、ネットワーク３１を備える。
情報処理装置１１と、それぞれの電波センサー２１－１１～２１－ＭＮは、ネットワーク３１と接続されている。そして、情報処理装置１１と、それぞれの電波センサー２１－１１～２１－ＭＮとは、ネットワーク３１を介して、通信することが可能である。
ネットワーク３１は、有線のネットワークであってもよく、無線のネットワークであってもよく、あるいは、有線と無線の両方を含むネットワークであってもよい。
＜情報処理システム＞
図１は、本発明の実施形態に係る情報処理システム１の概略的な構成を示す図である。
情報処理システム１は、情報処理装置１１と、複数の電波センサー２１－１１～２１－ＭＮと、ネットワーク３１を備える。
情報処理装置１１と、それぞれの電波センサー２１－１１～２１－ＭＮは、ネットワーク３１と接続されている。そして、情報処理装置１１と、それぞれの電波センサー２１－１１～２１－ＭＮとは、ネットワーク３１を介して、通信することが可能である。
ネットワーク３１は、有線のネットワークであってもよく、無線のネットワークであってもよく、あるいは、有線と無線の両方を含むネットワークであってもよい。
それぞれの電波センサー２１－１１～２１－ＭＮは、空間的に分散されて設置されている。そして、それぞれの電波センサー２１－１１～２１－ＭＮは、所定の周波数について電波の受信電力を検出（測定）し、その検出結果を、ネットワーク３１を介して、情報処理装置１１に送信する。
ここで、本実施形態では、ＭおよびＮをそれぞれ２以上の整数として、（Ｍ×Ｎ）個の電波センサー２１－１１～２１－ＭＮを示す。本実施形態では、説明の便宜上、（Ｍ×Ｎ）個の電波センサー２１－１１～２１－ＭＮは、所定の方向に等間隔にＭ行分が配置されており、各行ごとに当該所定の方向とは直交する方向に等間隔にＮ列分が配置されているとする。
なお、複数の電波センサーの数および配置は、それぞれ、任意であってもよい。ＭとＮとは、例えば、同じ値であってもよく、あるいは、異なる値であってもよい。
ここで、本実施形態では、ＭおよびＮをそれぞれ２以上の整数として、（Ｍ×Ｎ）個の電波センサー２１－１１～２１－ＭＮを示す。本実施形態では、説明の便宜上、（Ｍ×Ｎ）個の電波センサー２１－１１～２１－ＭＮは、所定の方向に等間隔にＭ行分が配置されており、各行ごとに当該所定の方向とは直交する方向に等間隔にＮ列分が配置されているとする。
なお、複数の電波センサーの数および配置は、それぞれ、任意であってもよい。ＭとＮとは、例えば、同じ値であってもよく、あるいは、異なる値であってもよい。
＜情報処理装置＞
情報処理装置１１は、入力部１１１と、出力部１１２と、通信部１１３と、記憶部１１４と、制御部１１５を備える。
制御部１１５は、取得部１３１と、推定演算部１３２と、出力制御部１３３を備える。
情報処理装置１１は、入力部１１１と、出力部１１２と、通信部１１３と、記憶部１１４と、制御部１１５を備える。
制御部１１５は、取得部１３１と、推定演算部１３２と、出力制御部１３３を備える。
入力部１１１は、例えば、ユーザ（人）によって操作される操作部と、外部の装置から直接的に情報を入力する機能部などを備える。
出力部１１２は、例えば、情報を表示出力する画面と、情報を音出力するスピーカと、外部の装置に直接的に出力する機能部などを備える。
通信部１１３は、ネットワーク３１を介して外部の装置（本実施形態では、電波センサー２１－１１～２１－ＭＮ）と通信する。
記憶部１１４は、情報を記憶する。
出力部１１２は、例えば、情報を表示出力する画面と、情報を音出力するスピーカと、外部の装置に直接的に出力する機能部などを備える。
通信部１１３は、ネットワーク３１を介して外部の装置（本実施形態では、電波センサー２１－１１～２１－ＭＮ）と通信する。
記憶部１１４は、情報を記憶する。
制御部１１５は、各種の処理または制御を行う。
取得部１３１は、情報を取得する。本実施形態では、取得部１３１は、通信部１１３によって電波センサー２１－１１～２１－ＭＮから受信された検出結果の情報を取得する。
推定演算部１３２は、取得部１３１によって取得された情報に基づいて、所定のアンテナの指向性を推定するための演算を行う。
出力制御部１３３は、推定演算部１３２によって行われた演算の結果を出力部１１２により出力するように制御する。
取得部１３１は、情報を取得する。本実施形態では、取得部１３１は、通信部１１３によって電波センサー２１－１１～２１－ＭＮから受信された検出結果の情報を取得する。
推定演算部１３２は、取得部１３１によって取得された情報に基づいて、所定のアンテナの指向性を推定するための演算を行う。
出力制御部１３３は、推定演算部１３２によって行われた演算の結果を出力部１１２により出力するように制御する。
＜本実施形態の前提条件＞
本実施形態では、１個の無線局のアンテナ指向性のパターンの方向を推定する。無線局としては、無線により通信を行う局が用いられ、例えば、携帯電話システムなどの基地局装置あるいは中継装置などが用いられてもよい。
また、無線局の位置（本実施形態では推定位置であり、他の例として、実際の位置であってもよい）、空中線電力、および、伝搬損失は既知であるとする。これらの情報は、例えば、既知の技術により推定などにより求められてもよく、あるいは、ユーザなどにより設定されてもよい。
また、推定対象となるアンテナ指向性のパターンは既知であるとする。当該パターンは、例えば、無線局の登録情報などから得られてもよく、あるいは、ユーザなどにより設定されてもよい。
また、アンテナ指向性のパターンは１種類のみであるとする。
また、空間として、地図上の２次元の平面を想定する。なお、本実施形態では、図示を簡易化するために、地図の情報の図示を省略するが、例えば、空間が表される場合に、２次元または３次元の地図の情報が組み合わされて表示等されてもよい。
本実施形態では、１個の無線局のアンテナ指向性のパターンの方向を推定する。無線局としては、無線により通信を行う局が用いられ、例えば、携帯電話システムなどの基地局装置あるいは中継装置などが用いられてもよい。
また、無線局の位置（本実施形態では推定位置であり、他の例として、実際の位置であってもよい）、空中線電力、および、伝搬損失は既知であるとする。これらの情報は、例えば、既知の技術により推定などにより求められてもよく、あるいは、ユーザなどにより設定されてもよい。
また、推定対象となるアンテナ指向性のパターンは既知であるとする。当該パターンは、例えば、無線局の登録情報などから得られてもよく、あるいは、ユーザなどにより設定されてもよい。
また、アンテナ指向性のパターンは１種類のみであるとする。
また、空間として、地図上の２次元の平面を想定する。なお、本実施形態では、図示を簡易化するために、地図の情報の図示を省略するが、例えば、空間が表される場合に、２次元または３次元の地図の情報が組み合わされて表示等されてもよい。
ここで、本実施形態では、無線局の位置が推定により取得されている場合を想定する。無線局の推定位置は、当該無線局の実際の位置とは異なり得る。
無線局の位置を推定する手法としては、任意の手法が用いられてもよく、例えば、既存のアルゴリズムを使用する手法が用いられてもよい。無線局の位置が推定される場合、推定誤差が発生し得る。位置の推定の精度は、無線局の位置を推定する手法に応じて決まり、例えば、当該推定誤差の最大値（最大推定誤差）は当該手法に応じて決まる。
なお、本実施形態に係る無線局のアンテナ指向性のパターンの推定は、無線局の位置が推定される場合に好適であるが、無線局の実際の位置が把握されている場合に適用されてもよい。
無線局の位置を推定する手法としては、任意の手法が用いられてもよく、例えば、既存のアルゴリズムを使用する手法が用いられてもよい。無線局の位置が推定される場合、推定誤差が発生し得る。位置の推定の精度は、無線局の位置を推定する手法に応じて決まり、例えば、当該推定誤差の最大値（最大推定誤差）は当該手法に応じて決まる。
なお、本実施形態に係る無線局のアンテナ指向性のパターンの推定は、無線局の位置が推定される場合に好適であるが、無線局の実際の位置が把握されている場合に適用されてもよい。
情報処理装置１１において、推定演算部１３２によって行われる［処理Ｔ１］～［処理Ｔ５］を説明する。
［処理Ｔ１］
推定演算部１３２は、複数の電波センサー２１－１１～２１－ＭＮから取得された検出結果に基づいて、空間的な電力分布を取得する。
当該電力分布は、無線局２１１が実際の位置にあるときに、無線局２１１から送信される電波に基づいて得られる電力分布である。
［処理Ｔ１］
推定演算部１３２は、複数の電波センサー２１－１１～２１－ＭＮから取得された検出結果に基づいて、空間的な電力分布を取得する。
当該電力分布は、無線局２１１が実際の位置にあるときに、無線局２１１から送信される電波に基づいて得られる電力分布である。
図２は、本発明の実施形態に係る無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布を示す図である。
ただし、図２の例では、説明の便宜上、無線局２１１が示される位置は、推定された位置となっている。当該電力分布は、無線局２１１が実際の位置にあるときに、無線局２１１から送信される電波に基づいて得られる電力分布である。
図２の例では、アンテナ指向性の推定対象となる１個の無線局２１１と、その周囲に存在する電波センサー２１－１１～２１－ＭＮを示してある。なお、図２の例では、推定演算に使用される電波センサー２１－１１～２１－ＭＮを示してあるが、他の電波センサーが情報処理システム１に備えられていてもよい。
また、図２の例では、図示を簡易化するために、それぞれの電波センサー２１－１１～２１－ＭＮによる検出結果（受信電力の強度）の概略を示してあるが、当該検出結果は受信電力の数値［ｄＢ］であり、当該数値が情報処理装置１１により取得される。
ただし、図２の例では、説明の便宜上、無線局２１１が示される位置は、推定された位置となっている。当該電力分布は、無線局２１１が実際の位置にあるときに、無線局２１１から送信される電波に基づいて得られる電力分布である。
図２の例では、アンテナ指向性の推定対象となる１個の無線局２１１と、その周囲に存在する電波センサー２１－１１～２１－ＭＮを示してある。なお、図２の例では、推定演算に使用される電波センサー２１－１１～２１－ＭＮを示してあるが、他の電波センサーが情報処理システム１に備えられていてもよい。
また、図２の例では、図示を簡易化するために、それぞれの電波センサー２１－１１～２１－ＭＮによる検出結果（受信電力の強度）の概略を示してあるが、当該検出結果は受信電力の数値［ｄＢ］であり、当該数値が情報処理装置１１により取得される。
また、図２の例では、図を見やすくするために、無線局２１１の位置を四角形を用いて表すが、当該四角形の中心が当該位置であるとする。同様に、図２の例では、図を見やすくするために、それぞれの電波センサー２１－１１～２１－ＭＮの位置を丸を用いて表すが、当該丸の中心が当該位置であるとする。これらについては、以降の図（図３～図５、図７～図８）においても同様である。
［処理Ｔ２］
推定演算部１３２は、無線局２１１のアンテナ指向性のパターンを仮の方向に配置した場合について、電波センサー２１－１１～２１－ＭＮによる検出結果（空間的な電力分布）を推定する。推定演算部１３２は、例えば、無線局２１１の空中線電力（送信機出力電力）、無線局２１１から電波センサー２１－１１～２１－ＭＮへの方向におけるアンテナ利得、伝搬損失に基づいて、電波センサー２１－１１～２１－ＭＮによる受信電力を推定する。
ここでは、無線局２１１の位置として、推定された位置が用いられる。
推定演算部１３２は、無線局２１１のアンテナ指向性のパターンを仮の方向に配置した場合について、電波センサー２１－１１～２１－ＭＮによる検出結果（空間的な電力分布）を推定する。推定演算部１３２は、例えば、無線局２１１の空中線電力（送信機出力電力）、無線局２１１から電波センサー２１－１１～２１－ＭＮへの方向におけるアンテナ利得、伝搬損失に基づいて、電波センサー２１－１１～２１－ＭＮによる受信電力を推定する。
ここでは、無線局２１１の位置として、推定された位置が用いられる。
図３は、本発明の実施形態に係る無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および仮のアンテナ指向性のパターン３１１の方向における電力分布を示す図である。
ここで、図３の例では、無線局２１１が示される位置は、推定された位置となっている。当該電力分布は、無線局２１１が推定された位置にあると仮定ときに得られると推定される電力分布である。
図３には、仮のアンテナ指向性のパターン３１１について、中心方向３３１と、角度幅３５１を示してある。
なお、図３の例では、図２の例と同様に、図示を簡易化するために、それぞれの電波センサー２１－１１～２１－ＭＮによる検出結果（受信電力の強度）の概略を示してあるが、当該検出結果は受信電力の数値［ｄＢ］である。
また、情報処理装置１１では、パターン３１１について、中心方向３３１および角度幅３５１は、例えば、パターン３１１の情報として、記憶部１１４に記憶され、制御部１１５によって処理される。
ここで、図３の例では、無線局２１１が示される位置は、推定された位置となっている。当該電力分布は、無線局２１１が推定された位置にあると仮定ときに得られると推定される電力分布である。
図３には、仮のアンテナ指向性のパターン３１１について、中心方向３３１と、角度幅３５１を示してある。
なお、図３の例では、図２の例と同様に、図示を簡易化するために、それぞれの電波センサー２１－１１～２１－ＭＮによる検出結果（受信電力の強度）の概略を示してあるが、当該検出結果は受信電力の数値［ｄＢ］である。
また、情報処理装置１１では、パターン３１１について、中心方向３３１および角度幅３５１は、例えば、パターン３１１の情報として、記憶部１１４に記憶され、制御部１１５によって処理される。
［処理Ｔ３］
推定演算部１３２は、電波センサー２１－１１～２１－ＭＮによる検出結果に基づいて得られた空間的な電力分布（測定値に基づく電力分布）と、無線局２１１について仮のアンテナ指向性のパターン３１１の方向を想定した場合に電波センサー２１－１１～２１－ＭＮにより得られると想定される空間的な電力分布（想定値に基づく電力分布）との相関係数を演算する。
推定演算部１３２は、電波センサー２１－１１～２１－ＭＮによる検出結果に基づいて得られた空間的な電力分布（測定値に基づく電力分布）と、無線局２１１について仮のアンテナ指向性のパターン３１１の方向を想定した場合に電波センサー２１－１１～２１－ＭＮにより得られると想定される空間的な電力分布（想定値に基づく電力分布）との相関係数を演算する。
ここで、本実施形態に係る情報処理装置１１では、図２および図３に示される電力分布から、無線局２１１の推定位置を含む所定の領域にある電力分布を除外して、無線局２１１のアンテナ指向性のパターンを推定する。換言すると、本実施形態に係る情報処理装置１１では、図２および図３に示される電力分布のうち、無線局２１１の推定位置を含む所定の領域の外側にある電力分布を抽出して、無線局２１１のアンテナ指向性のパターンを推定する。
本実施形態では、当該所定の領域として、無線局２１１の推定位置を中心とする半径ｍ（ｍは、０よりも大きい値）の円の内側の領域が用いられる。なお、当該所定の領域の境界線については、当該所定の領域の内側に含められてもよく、あるいは、含められなくてもよい。
本実施形態では、当該所定の領域として、無線局２１１の推定位置を中心とする半径ｍ（ｍは、０よりも大きい値）の円の内側の領域が用いられる。なお、当該所定の領域の境界線については、当該所定の領域の内側に含められてもよく、あるいは、含められなくてもよい。
本実施形態では、当該所定の領域は、固定的に設定されている。
ここで、円の半径ｍとしては、例えば、無線局の位置を推定するアルゴリズムにおける位置推定精度に応じた値が用いられてもよく、具体例として、当該位置推定精度に相当する値が用いられてもよく、あるいは、当該位置推定精度に所定の係数が乗じられた結果に相当する値（または、当該位置推定精度に相当する値に所定の係数が乗じられた結果）が用いられてもよい。
位置推定精度に相当する値としては、例えば、最大推定誤差の値が用いられてもよい。
なお、円の半径ｍとしては、位置推定精度に相当する値、または、それよりも大きい値とすることが、位置の推定精度の観点で、好ましい。つまり、無線局の推定位置を含む除外領域に、当該無線局の実際の位置が含まれることが好ましい。
ここで、円の半径ｍとしては、例えば、無線局の位置を推定するアルゴリズムにおける位置推定精度に応じた値が用いられてもよく、具体例として、当該位置推定精度に相当する値が用いられてもよく、あるいは、当該位置推定精度に所定の係数が乗じられた結果に相当する値（または、当該位置推定精度に相当する値に所定の係数が乗じられた結果）が用いられてもよい。
位置推定精度に相当する値としては、例えば、最大推定誤差の値が用いられてもよい。
なお、円の半径ｍとしては、位置推定精度に相当する値、または、それよりも大きい値とすることが、位置の推定精度の観点で、好ましい。つまり、無線局の推定位置を含む除外領域に、当該無線局の実際の位置が含まれることが好ましい。
図４は、本発明の実施形態（第１実施形態）に係る所定の領域が除外された場合における無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布を示す図である。
図４の例では、図２に示される無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布と比べて、所定の領域における電力分布が除かれた様子を示してある。
図４の例では、当該所定の領域は、無線局２１１の推定位置を中心とする半径ｍの円の領域である円領域４１１の内側の領域である。
図４の例では、図２に示される無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布と比べて、所定の領域における電力分布が除かれた様子を示してある。
図４の例では、当該所定の領域は、無線局２１１の推定位置を中心とする半径ｍの円の領域である円領域４１１の内側の領域である。
図５は、本発明の実施形態（第１実施形態）に係る所定の領域が除外された場合における無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および仮のアンテナ指向性のパターン３１１の方向における電力分布を示す図である。
図５の例では、図３に示される無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布と比べて、所定の領域における電力分布が除かれた様子を示してある。
図５の例では、当該所定の領域は、無線局２１１の推定位置を中心とする半径ｍの円の領域である円領域４２１の内側の領域である。
ここで、図４に示される円領域４１１と図５に示される円領域４２１とは同じである。
図５の例では、図３に示される無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布と比べて、所定の領域における電力分布が除かれた様子を示してある。
図５の例では、当該所定の領域は、無線局２１１の推定位置を中心とする半径ｍの円の領域である円領域４２１の内側の領域である。
ここで、図４に示される円領域４１１と図５に示される円領域４２１とは同じである。
本実施形態に係る情報処理装置１１では、説明の便宜上、電波センサー２１－１１～２１－ＭＮのうちで円領域４１１（あるいは、円領域４２１でも同様）の外側に位置する電波センサーを、有効な電波センサーと呼んで説明する。
換言すると、有効な電波センサーは、電波センサー２１－１１～２１－ＭＮから、円領域４１１（あるいは、円領域４２１でも同様）の内側に位置する電波センサーを除いた電波センサーである。
換言すると、有効な電波センサーは、電波センサー２１－１１～２１－ＭＮから、円領域４１１（あるいは、円領域４２１でも同様）の内側に位置する電波センサーを除いた電波センサーである。
ここで、本実施形態では、それぞれの電波センサー２１－１１～２１－ＭＮの位置は既知であり、無線局２１１の推定位置は既知であり、これらの位置を表す情報が情報処理装置１１に設定される。情報処理装置１１では、当該情報に基づいて、無線局２１１の推定位置とそれぞれの電波センサー２１－１１～２１－ＭＮの位置との距離を計算することが可能である。なお、位置を表す情報としては、任意の情報が用いられてもよく、例えば、緯度および経度の情報が用いられてもよい。
本実施形態では、式（１）に示される相関係数ｒが用いられる。
式（１）において、ｔ（ｔ＝１～Ｔ）は変数であり、Ｔは電波センサー２１－１１～２１－ＭＮのうちで有効な電波センサーの数（本実施形態では、Ｍ×Ｎよりも小さい数）を表しており、Ｐ０，ｔはｔ番目の電波センサーによる受信電力（測定値）を表しており、Ｐ１，ｔは仮のアンテナ指向性のパターン３１１の方向を想定した場合におけるｔ番目の電波センサーによる受信電力（想定値）を表している。また、Ｐ０，ｔの上に横線が引かれたものはＰ０，ｔの平均値を表しており、Ｐ１，ｔの上に横線が引かれたものはＰ１，ｔの平均値を表している。
ここで、ｔ＝１～Ｔのそれぞれに対応する複数の電波センサー２１－１１～２１－ＭＮの順序は任意であってもよい。
なお、相関係数としては、例えば、他のパラメータを含む演算式が用いられてもよく、あるいは、同じパラメータであっても他の演算式による相関係数が用いられてもよい。
式（１）において、ｔ（ｔ＝１～Ｔ）は変数であり、Ｔは電波センサー２１－１１～２１－ＭＮのうちで有効な電波センサーの数（本実施形態では、Ｍ×Ｎよりも小さい数）を表しており、Ｐ０，ｔはｔ番目の電波センサーによる受信電力（測定値）を表しており、Ｐ１，ｔは仮のアンテナ指向性のパターン３１１の方向を想定した場合におけるｔ番目の電波センサーによる受信電力（想定値）を表している。また、Ｐ０，ｔの上に横線が引かれたものはＰ０，ｔの平均値を表しており、Ｐ１，ｔの上に横線が引かれたものはＰ１，ｔの平均値を表している。
ここで、ｔ＝１～Ｔのそれぞれに対応する複数の電波センサー２１－１１～２１－ＭＮの順序は任意であってもよい。
なお、相関係数としては、例えば、他のパラメータを含む演算式が用いられてもよく、あるいは、同じパラメータであっても他の演算式による相関係数が用いられてもよい。
このように、本実施形態では、相関係数を求めるときに、演算対象とする電波センサー２１－１１～２１－ＭＮの範囲が一部の電波センサーに設定されている。
具体的には、電力分布の抽出は、それぞれの電波センサー２１－１１～２１－ＭＮを単位として、電波センサー２１－１１～２１－ＭＮを抽出することにより行われる。つまり、本実施形態では、所定の領域の外側にある電波センサーを相関係数の演算対象とし、当該所定の領域の内側にある電波センサーを相関係数の演算対象外とする。
具体的には、電力分布の抽出は、それぞれの電波センサー２１－１１～２１－ＭＮを単位として、電波センサー２１－１１～２１－ＭＮを抽出することにより行われる。つまり、本実施形態では、所定の領域の外側にある電波センサーを相関係数の演算対象とし、当該所定の領域の内側にある電波センサーを相関係数の演算対象外とする。
ここで、電波センサー２１－１１～２１－ＭＮによる受信電力（測定値と想定値の一方または両方）に重み付けが行われてもよい。
例えば、相関係数の演算において、無線局２１１からの距離に応じた重み付けが、電波センサー２１－１１～２１－ＭＮによる受信電力（測定値と想定値の一方または両方）に行われてもよい。具体例として、無線局２１１に近いほど重み付けの値を高くし、無線局２１１から遠いほど重み付け値を低くしてもよく、これにより、無線局２１１に対して遠方に存在する別の無線局の影響を低減し、また、無線局２１１からの電力が低い領域でのノイズの影響を低減することができる。逆に、無線局２１１に近いほど重み付けの値を低くし、無線局２１１から遠いほど重み付け値を高くしてもよい。
例えば、相関係数の演算において、無線局２１１からの距離に応じた重み付けが、電波センサー２１－１１～２１－ＭＮによる受信電力（測定値と想定値の一方または両方）に行われてもよい。具体例として、無線局２１１に近いほど重み付けの値を高くし、無線局２１１から遠いほど重み付け値を低くしてもよく、これにより、無線局２１１に対して遠方に存在する別の無線局の影響を低減し、また、無線局２１１からの電力が低い領域でのノイズの影響を低減することができる。逆に、無線局２１１に近いほど重み付けの値を低くし、無線局２１１から遠いほど重み付け値を高くしてもよい。
［処理Ｔ４］
推定演算部１３２は、仮のアンテナ指向性のパターン３１１の方向を初期値から変化させながら、それぞれの仮定の方向について相関係数を演算する。
推定演算部１３２は、仮のアンテナ指向性のパターン３１１の方向を初期値から変化させながら、それぞれの仮定の方向について相関係数を演算する。
［処理Ｔ５］
推定演算部１３２は、算出された相関係数の値が最も高くなったアンテナ指向性のパターン３１１の方向を、推定結果とする。
推定演算部１３２は、算出された相関係数の値が最も高くなったアンテナ指向性のパターン３１１の方向を、推定結果とする。
ここで、相関係数を演算する仮定の方向の数としては、任意に設定されてもよい。
また、仮定の方向を変化させる手法としては、任意の手法が用いられてもよい。
また、仮定の方向を変化させる手法としては、任意の手法が用いられてもよい。
一例として、均等なステップで仮定の方向を変化させる手法が用いられてもよい。
図６は、本発明の実施形態に係る情報処理装置１１において行われる処理の手順の一例を示す図である。
図６の例では、無線局２１１（１点とみなす）を中心として全方向を３６０度の角度として、これをＬ（Ｌは２以上の整数）で等分割したステップ（一定の角度ステップ）が用いられている。また、ｎ（ｎ＝１～Ｌ）を変数とする。Ｌは、例えば、３６などであってもよい。
図６は、本発明の実施形態に係る情報処理装置１１において行われる処理の手順の一例を示す図である。
図６の例では、無線局２１１（１点とみなす）を中心として全方向を３６０度の角度として、これをＬ（Ｌは２以上の整数）で等分割したステップ（一定の角度ステップ）が用いられている。また、ｎ（ｎ＝１～Ｌ）を変数とする。Ｌは、例えば、３６などであってもよい。
（ステップＳ１）
推定演算部１３２は、ｎ＝０に設定する。そして、推定演算部１３２は、ステップＳ２の処理へ移行する。
推定演算部１３２は、ｎ＝０に設定する。そして、推定演算部１３２は、ステップＳ２の処理へ移行する。
（ステップＳ２）
推定演算部１３２は、相関係数ｒ（θｎ）を計算する。そして、推定演算部１３２は、ステップＳ３の処理へ移行する。ここで、相関係数ｒ（θｎ）は、仮定の方向を角度θｎとした場合における相関係数である。
推定演算部１３２は、相関係数ｒ（θｎ）を計算する。そして、推定演算部１３２は、ステップＳ３の処理へ移行する。ここで、相関係数ｒ（θｎ）は、仮定の方向を角度θｎとした場合における相関係数である。
（ステップＳ３）
推定演算部１３２は、ｎに１を追加（加算）する。そして、推定演算部１３２は、ステップＳ４の処理へ移行する。
推定演算部１３２は、ｎに１を追加（加算）する。そして、推定演算部１３２は、ステップＳ４の処理へ移行する。
（ステップＳ４）
推定演算部１３２は、ｎがＬよりも小さいか否かを判定する。
この判定の結果、推定演算部１３２は、ｎがＬよりも小さいと判定した場合（ステップＳ４：ＹＥＳ）、ステップＳ２の処理へ移行する。
一方、この判定の結果、推定演算部１３２は、ｎがＬよりも小さくないと判定した場合（ステップＳ４：ＮＯ）、ステップＳ５の処理へ移行する。
推定演算部１３２は、ｎがＬよりも小さいか否かを判定する。
この判定の結果、推定演算部１３２は、ｎがＬよりも小さいと判定した場合（ステップＳ４：ＹＥＳ）、ステップＳ２の処理へ移行する。
一方、この判定の結果、推定演算部１３２は、ｎがＬよりも小さくないと判定した場合（ステップＳ４：ＮＯ）、ステップＳ５の処理へ移行する。
（ステップＳ５）
推定演算部１３２は、算出されたＬ個の相関係数のうち最大となる相関係数ＭＡＸ（ｒ（θｎ））が求められた角度θｎをアンテナ指向性のパターン３１１の方向の推定結果と判定して、当該推定結果の情報を出力する。そして、推定演算部１３２は、本フローの処理を終了する。
推定演算部１３２は、算出されたＬ個の相関係数のうち最大となる相関係数ＭＡＸ（ｒ（θｎ））が求められた角度θｎをアンテナ指向性のパターン３１１の方向の推定結果と判定して、当該推定結果の情報を出力する。そして、推定演算部１３２は、本フローの処理を終了する。
他の例として、相関係数が高い領域の角度分割を繰り返す手法が用いられてもよい。
例えば、アンテナ指向性のパターン３１１の方向の角度を定義する。一例として、北の方向を０度（°）とし、右回り（時計回り）に３６０度の角度を定義する。この場合、東の方向が９０度となり、南の方向が１８０度となり、西の方向が２７０度となる。
推定演算部１３２は、初期の領域から始めて、対象となる領域（対象領域）を２つ以上の分割領域に分けて、それぞれの分割領域ごとに１つの仮定の方向を設定して相関係数を計算し、相関係数が最も高い仮定の方向を含む分割領域を次の対象領域とする。そして、推定演算部１３２は、最終的に絞り込まれた分割領域に含まれる仮定の方向を採用する。初期の領域としては、例えば、３６０度の領域が用いられてもよい。
例えば、アンテナ指向性のパターン３１１の方向の角度を定義する。一例として、北の方向を０度（°）とし、右回り（時計回り）に３６０度の角度を定義する。この場合、東の方向が９０度となり、南の方向が１８０度となり、西の方向が２７０度となる。
推定演算部１３２は、初期の領域から始めて、対象となる領域（対象領域）を２つ以上の分割領域に分けて、それぞれの分割領域ごとに１つの仮定の方向を設定して相関係数を計算し、相関係数が最も高い仮定の方向を含む分割領域を次の対象領域とする。そして、推定演算部１３２は、最終的に絞り込まれた分割領域に含まれる仮定の方向を採用する。初期の領域としては、例えば、３６０度の領域が用いられてもよい。
以上のように、本実施形態に係る情報処理システム１における情報処理装置１１では、無線システムにおいて無線送信を行う無線局２１１のアンテナ指向性のパターン３１１を推定する。
情報処理装置１１では、所定の領域を除外して、複数の電波センサー２１－１１～２１－ＭＮを用いて取得された空間的な電力分布と、無線局２１１のアンテナ指向性のパターン３１１の方向を仮の方向としたときに電波センサー２１－１１～２１－ＭＮで得られると想定される空間的な電力分布との相関係数を求める。そして、情報処理装置１１では、仮の方向を変化させて、相関係数の値が最も高くなった仮の方向を、推定結果の方向とする。
情報処理装置１１では、所定の領域を除外して、複数の電波センサー２１－１１～２１－ＭＮを用いて取得された空間的な電力分布と、無線局２１１のアンテナ指向性のパターン３１１の方向を仮の方向としたときに電波センサー２１－１１～２１－ＭＮで得られると想定される空間的な電力分布との相関係数を求める。そして、情報処理装置１１では、仮の方向を変化させて、相関係数の値が最も高くなった仮の方向を、推定結果の方向とする。
したがって、情報処理装置１１では、無線局２１１の推定位置が用いられる場合においても、アンテナ指向性のパターン３１１の方向を高精度に推定することが可能となる。そして、情報処理装置１１では、例えば、電波センサー２１－１１～２１－ＭＮがない方向の場所の受信電力の推定誤差を低減して、当該受信電力を高精度に推定することができる。
なお、アンテナ指向性のパターン３１１の推定、あるいは、受信電力の推定においては、例えば、無線局２１１の空中線電力（送信機出力電力）、無線局２１１のアンテナから見た対象となる場所の方向のアンテナ利得、伝搬損失などの情報が用いられてもよく、これにより、推定の精度を向上させることができる。
なお、アンテナ指向性のパターン３１１の推定、あるいは、受信電力の推定においては、例えば、無線局２１１の空中線電力（送信機出力電力）、無線局２１１のアンテナから見た対象となる場所の方向のアンテナ利得、伝搬損失などの情報が用いられてもよく、これにより、推定の精度を向上させることができる。
ここで、無線局２１１の推定位置の近傍を相関計算の範囲から除外することの効果について説明する。
無線局２１１の推定位置に誤差がある場合に、アンテナ指向性のパターンの推定精度に影響を与え得る。
例えば、無線局２１１の実際の位置の近傍の範囲では電力が高いため、当該範囲を除外した方が、相関係数の精度が良いと考えられる。つまり、相関係数を計算する範囲に、無線局２１１の推定位置と実際の位置との間にある電力が高い部分が含まれると、電力が高い部分での相関にずれが出て、相関係数の精度に影響が生じると考えられる。
無線局２１１の推定位置に誤差がある場合に、アンテナ指向性のパターンの推定精度に影響を与え得る。
例えば、無線局２１１の実際の位置の近傍の範囲では電力が高いため、当該範囲を除外した方が、相関係数の精度が良いと考えられる。つまり、相関係数を計算する範囲に、無線局２１１の推定位置と実際の位置との間にある電力が高い部分が含まれると、電力が高い部分での相関にずれが出て、相関係数の精度に影響が生じると考えられる。
一例として、無線局２１１の推定位置が、実際の位置に対して、アンテナの指向方向にずれた場合に、（本実施形態が適用されない場合における）相関係数の精度への影響が大きいと考えられ、本実施形態により当該影響を低減する効果が顕著である。
一方、無線局２１１の推定位置が、実際の位置に対して、アンテナの指向方向とは逆方向にずれた場合には、相関係数の精度の劣化は小さいと考えられるが、本実施形態を適用しても、相関係数の精度が十分に得られると考えられる。
また、無線局２１１の推定位置が、実際の位置に対して、他の方向にずれた場合には、当該他の方向に応じて、（本実施形態が適用されない場合における）相関係数の精度への影響が生じると考えられ、本実施形態の効果が得られる。
このようなことから、無線局２１１の推定位置のずれ方向に関わらず、本実施形態が有効である。
一方、無線局２１１の推定位置が、実際の位置に対して、アンテナの指向方向とは逆方向にずれた場合には、相関係数の精度の劣化は小さいと考えられるが、本実施形態を適用しても、相関係数の精度が十分に得られると考えられる。
また、無線局２１１の推定位置が、実際の位置に対して、他の方向にずれた場合には、当該他の方向に応じて、（本実施形態が適用されない場合における）相関係数の精度への影響が生じると考えられ、本実施形態の効果が得られる。
このようなことから、無線局２１１の推定位置のずれ方向に関わらず、本実施形態が有効である。
（第２実施形態）
本実施形態では、第１実施形態とは異なる点について詳しく説明し、同様な点については説明を省略する。
本実施形態では、概略的には、第１実施形態と比べて、アンテナ指向性のパターンを推定する際に、［処理Ｔ３］において、図４および図５に示される所定の領域として異なる領域を用いる。
本実施形態では、第１実施形態とは異なる点について詳しく説明し、同様な点については説明を省略する。
本実施形態では、概略的には、第１実施形態と比べて、アンテナ指向性のパターンを推定する際に、［処理Ｔ３］において、図４および図５に示される所定の領域として異なる領域を用いる。
ここで、本実施形態に係る情報処理装置１１では、図２および図３に示される電力分布から、無線局２１１の推定位置を含む所定の第１領域にある電力分布を除外するとともに、無線局２１１の推定位置を含まない所定の第２領域にある電力分布を除外して、無線局２１１のアンテナ指向性のパターンを推定する。換言すると、本実施形態に係る情報処理装置１１では、図２および図３に示される電力分布のうち、無線局２１１の推定位置を含む所定の第１領域の外側にあり、かつ、無線局２１１の推定位置を含まない所定の第２領域にはない電力分布を抽出して、無線局２１１のアンテナ指向性のパターンを推定する。
本実施形態では、これらの除外される領域を所定の領域とする。
本実施形態では、当該所定の第１領域として、無線局２１１の推定位置を中心とする半径ｍ１（ｍ１は、０よりも大きい値）の円（第１円）の内側の領域が用いられる。なお、当該所定の第１領域の境界線については、当該所定の第１領域の内側に含められてもよく、あるいは、含められなくてもよい。
また、本実施形態では、当該所定の第２領域として、無線局２１１の推定位置を中心とする半径ｍ２（ｍ２は、ｍ１よりも大きい値）の円（第２円）の外側の領域が用いられる。なお、当該所定の第２領域の境界線については、当該所定の第２領域の内側に含められてもよく、あるいは、含められなくてもよい。
そして、本実施形態では、当該所定の領域として、当該所定の第１領域と当該所定の第２領域とを合わせた領域が用いられる。
本実施形態では、当該所定の第１領域として、無線局２１１の推定位置を中心とする半径ｍ１（ｍ１は、０よりも大きい値）の円（第１円）の内側の領域が用いられる。なお、当該所定の第１領域の境界線については、当該所定の第１領域の内側に含められてもよく、あるいは、含められなくてもよい。
また、本実施形態では、当該所定の第２領域として、無線局２１１の推定位置を中心とする半径ｍ２（ｍ２は、ｍ１よりも大きい値）の円（第２円）の外側の領域が用いられる。なお、当該所定の第２領域の境界線については、当該所定の第２領域の内側に含められてもよく、あるいは、含められなくてもよい。
そして、本実施形態では、当該所定の領域として、当該所定の第１領域と当該所定の第２領域とを合わせた領域が用いられる。
ここで、第１円の半径ｍ１としては、例えば、無線局の位置を推定するアルゴリズムにおける位置推定精度に応じた値が用いられてもよく、具体例として、当該位置推定精度に相当する値が用いられてもよく、あるいは、当該位置推定精度に所定の係数が乗じられた結果に相当する値（または、当該位置推定精度に相当する値に所定の係数が乗じられた結果）が用いられてもよい。
位置推定精度に相当する値としては、例えば、最大推定誤差の値が用いられてもよい。
なお、第１円の半径ｍ１としては、位置推定精度に相当する値、または、それよりも大きい値とすることが、位置の推定精度の観点で、好ましい。つまり、無線局の推定位置を含む除外領域に、当該無線局の実際の位置が含まれることが好ましい。
位置推定精度に相当する値としては、例えば、最大推定誤差の値が用いられてもよい。
なお、第１円の半径ｍ１としては、位置推定精度に相当する値、または、それよりも大きい値とすることが、位置の推定精度の観点で、好ましい。つまり、無線局の推定位置を含む除外領域に、当該無線局の実際の位置が含まれることが好ましい。
また、第２円の半径ｍ２としては、例えば、所定の幅ｄを定めて、（ｍ２＝ｍ１＋ｄ）により設定されてもよい。
なお、第２円の半径ｍ２としては、第２円の内側では受信電力がノイズに埋もれず、第２円の外側では受信電力がノイズに埋もれるような値とすることが、位置の推定精度の観点で、好ましい。つまり、受信電力がノイズに埋もれる領域は、除外されることが好ましい。受信電力がノイズに埋もれる領域は、例えば、受信電力のレベルと、ノイズのレベルに基づいて設定されてもよく、一例として、無線局からの送信電力と、伝搬損失に基づいて設定されてもよい。
なお、第２円の半径ｍ２としては、第２円の内側では受信電力がノイズに埋もれず、第２円の外側では受信電力がノイズに埋もれるような値とすることが、位置の推定精度の観点で、好ましい。つまり、受信電力がノイズに埋もれる領域は、除外されることが好ましい。受信電力がノイズに埋もれる領域は、例えば、受信電力のレベルと、ノイズのレベルに基づいて設定されてもよく、一例として、無線局からの送信電力と、伝搬損失に基づいて設定されてもよい。
図７は、本発明の実施形態（第２実施形態）に係る所定の領域が除外された場合における無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布を示す図である。
図７の例では、図２に示される無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布と比べて、所定の領域における電力分布が除かれた様子を示してある。
図７の例では、当該所定の領域は、無線局２１１の推定位置を中心とする半径ｍ１の円の領域である第１円領域５１１の内側の第１領域と、無線局２１１の推定位置を中心とする半径ｍ２の円の領域である第２円領域５１２の外側の第２領域と、を合わせた領域である。
図７の例では、図２に示される無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布と比べて、所定の領域における電力分布が除かれた様子を示してある。
図７の例では、当該所定の領域は、無線局２１１の推定位置を中心とする半径ｍ１の円の領域である第１円領域５１１の内側の第１領域と、無線局２１１の推定位置を中心とする半径ｍ２の円の領域である第２円領域５１２の外側の第２領域と、を合わせた領域である。
図８は、本発明の実施形態（第２実施形態）に係る所定の領域が除外された場合における無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および仮のアンテナ指向性のパターン３１１の方向における電力分布を示す図である。
図８の例では、図２に示される無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布と比べて、所定の領域における電力分布が除かれた様子を示してある。
図８の例では、当該所定の領域は、無線局２１１の推定位置を中心とする半径ｍ１の円の領域である第１円領域５２１の内側の第１領域と、無線局２１１の推定位置を中心とする半径ｍ２の円の領域である第２円領域５２２の外側の第２領域と、を合わせた領域である。
ここで、図７に示される第１円領域５１１と図８に示される第１円領域５２１とは同じである。
また、図７に示される第２円領域５１２と図８に示される第２円領域５２２とは同じである。
図８の例では、図２に示される無線局２１１と電波センサー２１－１１～２１－ＭＮの配置および電力分布と比べて、所定の領域における電力分布が除かれた様子を示してある。
図８の例では、当該所定の領域は、無線局２１１の推定位置を中心とする半径ｍ１の円の領域である第１円領域５２１の内側の第１領域と、無線局２１１の推定位置を中心とする半径ｍ２の円の領域である第２円領域５２２の外側の第２領域と、を合わせた領域である。
ここで、図７に示される第１円領域５１１と図８に示される第１円領域５２１とは同じである。
また、図７に示される第２円領域５１２と図８に示される第２円領域５２２とは同じである。
本実施形態に係る情報処理装置１１では、電波センサー２１－１１～２１－ＭＮのうちで所定の領域の外側に位置する電波センサーを有効な電波センサーとして、有効な電波センサーについて相関係数を計算する。
換言すると、有効な電波センサーは、電波センサー２１－１１～２１－ＭＮから、所定の領域の内側に位置する電波センサーを除いた電波センサーである。
具体的には、有効な電波センサーは、無線局２１１の推定位置を中心とする半径ｍ１の円の外側であり、かつ、無線局２１１の推定位置を中心とする半径ｍ２の円の内側である領域である。
換言すると、有効な電波センサーは、電波センサー２１－１１～２１－ＭＮから、所定の領域の内側に位置する電波センサーを除いた電波センサーである。
具体的には、有効な電波センサーは、無線局２１１の推定位置を中心とする半径ｍ１の円の外側であり、かつ、無線局２１１の推定位置を中心とする半径ｍ２の円の内側である領域である。
このように、本実施形態では、相関係数を求めるときに、演算対象とする電波センサー２１－１１～２１－ＭＮの範囲が一部の電波センサーに設定されている。
具体的には、電力分布の抽出は、それぞれの電波センサー２１－１１～２１－ＭＮを単位として、電波センサー２１－１１～２１－ＭＮを抽出することにより行われる。つまり、本実施形態では、所定の領域の外側にある電波センサーを相関係数の演算対象とし、当該所定の領域の内側にある電波センサーを相関係数の演算対象外とする。
具体的には、電力分布の抽出は、それぞれの電波センサー２１－１１～２１－ＭＮを単位として、電波センサー２１－１１～２１－ＭＮを抽出することにより行われる。つまり、本実施形態では、所定の領域の外側にある電波センサーを相関係数の演算対象とし、当該所定の領域の内側にある電波センサーを相関係数の演算対象外とする。
以上のように、本実施形態に係る情報処理システム１における情報処理装置１１では、所定の領域を除外した領域（本実施形態では、同心円で挟まれた領域）について相関係数を計算することで、無線システムにおいて無線送信を行う無線局２１１のアンテナ指向性のパターン３１１を推定する。
したがって、情報処理装置１１では、第１実施形態と同様な効果を得ることができ、さらに、無線局２１１の推定位置から遠方におけるノイズの影響を低減して、相関係数の精度への影響を低減することができる。
したがって、情報処理装置１１では、第１実施形態と同様な効果を得ることができ、さらに、無線局２１１の推定位置から遠方におけるノイズの影響を低減して、相関係数の精度への影響を低減することができる。
（第３実施形態）
本実施形態では、第１実施形態～第２実施形態とは異なる点について詳しく説明し、同様な点については説明を省略する。
本実施形態では、概略的には、第１実施形態～第２実施形態と比べて、アンテナ指向性のパターンを推定する際に、アンテナ指向性のパターンの仮定の方向を変化させるとともに、相関係数の演算対象外とする所定の領域を変化させる。つまり、本実施形態では、アンテナ指向性のパターンの仮定の方向および相関係数の演算対象外とする所定の領域を変化させて、相関係数の最大値を求める。
本実施形態では、第１実施形態～第２実施形態とは異なる点について詳しく説明し、同様な点については説明を省略する。
本実施形態では、概略的には、第１実施形態～第２実施形態と比べて、アンテナ指向性のパターンを推定する際に、アンテナ指向性のパターンの仮定の方向を変化させるとともに、相関係数の演算対象外とする所定の領域を変化させる。つまり、本実施形態では、アンテナ指向性のパターンの仮定の方向および相関係数の演算対象外とする所定の領域を変化させて、相関係数の最大値を求める。
図９は、本発明の実施形態（第３実施形態）に係る情報処理装置１１において行われる処理の手順の一例を示す図である。
図９に示されるフローの処理は、例えば、図６に示されるステップＳ２の処理において、相関関数ｒ（θｎ）を計算する際に行われる。
なお、図６に示されるフローと図９に示されるフローとの組み合わせの仕方は、例えば、同様な結果が得られれば、任意の組み合わせの処理手順が用いられてもよい。
図９に示されるフローの処理は、例えば、図６に示されるステップＳ２の処理において、相関関数ｒ（θｎ）を計算する際に行われる。
なお、図６に示されるフローと図９に示されるフローとの組み合わせの仕方は、例えば、同様な結果が得られれば、任意の組み合わせの処理手順が用いられてもよい。
図９の例において、第１実施形態における図４および図５の例に適用される場合には、円の半径は、円領域４１１の半径ｍ（円領域４２１の半径ｍも同じ）に相当する。
また、図９の例において、第２実施形態における図７および図８の例に適用される場合には、円の半径は、第１円領域５１１の半径ｍ１（第１円領域５２１の半径ｍ１も同じ）と、第２円領域５１２の半径ｍ２（第２円領域５２２の半径ｍ２も同じ）との一方または両方である。つまり、図９の例において、第２実施形態における図７および図８の例に適用される場合には、２つの円のうちの任意の一方に適用されてもよく、あるいは、両方に適用されてもよい。
また、図９の例において、第２実施形態における図７および図８の例に適用される場合には、円の半径は、第１円領域５１１の半径ｍ１（第１円領域５２１の半径ｍ１も同じ）と、第２円領域５１２の半径ｍ２（第２円領域５２２の半径ｍ２も同じ）との一方または両方である。つまり、図９の例において、第２実施形態における図７および図８の例に適用される場合には、２つの円のうちの任意の一方に適用されてもよく、あるいは、両方に適用されてもよい。
（ステップＳ２１）
推定演算部１３２は、円の半径を初期値に設定する。そして、推定演算部１３２は、ステップＳ２２の処理へ移行する。
ここで、本実施形態では、円の半径の初期値、および、当該円の半径を変更する範囲が、あらかじめ設定されている。この設定内容は、例えば、情報処理装置１１にあらかじめ記憶されていてもよく、あるいは、ユーザの操作によって情報処理装置１１に可変に設定されてもよい。
推定演算部１３２は、円の半径を初期値に設定する。そして、推定演算部１３２は、ステップＳ２２の処理へ移行する。
ここで、本実施形態では、円の半径の初期値、および、当該円の半径を変更する範囲が、あらかじめ設定されている。この設定内容は、例えば、情報処理装置１１にあらかじめ記憶されていてもよく、あるいは、ユーザの操作によって情報処理装置１１に可変に設定されてもよい。
（ステップＳ２２）
推定演算部１３２は、相関係数ｒ（θｎ）を計算する。そして、推定演算部１３２は、ステップＳ２３の処理へ移行する。ここで、相関係数ｒ（θｎ）は、仮定の方向を角度θｎとし、かつ、円の半径が現在設定されている値である場合における相関係数である。
推定演算部１３２は、相関係数ｒ（θｎ）を計算する。そして、推定演算部１３２は、ステップＳ２３の処理へ移行する。ここで、相関係数ｒ（θｎ）は、仮定の方向を角度θｎとし、かつ、円の半径が現在設定されている値である場合における相関係数である。
（ステップＳ２３）
推定演算部１３２は、円の半径を変更する。そして、推定演算部１３２は、ステップＳ２４の処理へ移行する。
ここで、円の半径を変更する態様は、例えば、あらかじめ定められている。
一例として、円の半径を変更する態様として、現在の値に所定値（ここでは、正の値とする）を加算した結果へ変更することで、円の半径を次第に大きくしていく態様が用いられてもよい。
他の例として、円の半径を変更する態様として、現在の値から所定値（ここでは、正の値とする）を減算した結果へ変更することで、円の半径を次第に小さくしていく態様が用いられてもよい。
推定演算部１３２は、円の半径を変更する。そして、推定演算部１３２は、ステップＳ２４の処理へ移行する。
ここで、円の半径を変更する態様は、例えば、あらかじめ定められている。
一例として、円の半径を変更する態様として、現在の値に所定値（ここでは、正の値とする）を加算した結果へ変更することで、円の半径を次第に大きくしていく態様が用いられてもよい。
他の例として、円の半径を変更する態様として、現在の値から所定値（ここでは、正の値とする）を減算した結果へ変更することで、円の半径を次第に小さくしていく態様が用いられてもよい。
なお、第２実施形態に係る図７および図８の例において、２つの円の半径ｍ１、ｍ２の両方を変更する場合には、例えば、両方の円の半径ｍ１、ｍ２を同じ態様で変更してもよく、あるいは、それぞれの円の半径ｍ１、ｍ２ごとに異なる態様で変更してもよい。
また、第２実施形態に係る図７および図８の例において、２つの円の半径ｍ１、ｍ２の両方を変更する場合には、例えば、２つの円の半径ｍ１、ｍ２のうちの一方を固定して他方を変更する処理手順および当該他方を固定して当該一方を変更する処理手順の組み合わせが用いられてもよく、あるいは、２つの円の半径ｍ１、ｍ２の両方を同時に変更する処理手順が用いられてもよい。
また、第２実施形態に係る図７および図８の例において、２つの円の半径ｍ１、ｍ２の両方を変更する場合には、例えば、２つの円の半径ｍ１、ｍ２のうちの一方を固定して他方を変更する処理手順および当該他方を固定して当該一方を変更する処理手順の組み合わせが用いられてもよく、あるいは、２つの円の半径ｍ１、ｍ２の両方を同時に変更する処理手順が用いられてもよい。
（ステップＳ２４）
推定演算部１３２は、変更後の円の半径が、当該円の半径を変更する範囲を超えたか否かを判定する。
この判定の結果、推定演算部１３２は、変更後の円の半径が、当該円の半径を変更する範囲を超えたと判定したときには、当該円の半径の変更は終了したと判定し（ステップＳ２４：ＹＥＳ）、ステップＳ２５の処理へ移行する。
一方、この判定の結果、推定演算部１３２は、変更後の円の半径が、当該円の半径を変更する範囲を超えていないと判定したときには、当該円の半径の変更は終了していないと判定し（ステップＳ２４：ＮＯ）、ステップＳ２２の処理へ移行する。
推定演算部１３２は、変更後の円の半径が、当該円の半径を変更する範囲を超えたか否かを判定する。
この判定の結果、推定演算部１３２は、変更後の円の半径が、当該円の半径を変更する範囲を超えたと判定したときには、当該円の半径の変更は終了したと判定し（ステップＳ２４：ＹＥＳ）、ステップＳ２５の処理へ移行する。
一方、この判定の結果、推定演算部１３２は、変更後の円の半径が、当該円の半径を変更する範囲を超えていないと判定したときには、当該円の半径の変更は終了していないと判定し（ステップＳ２４：ＮＯ）、ステップＳ２２の処理へ移行する。
（ステップＳ２５）
推定演算部１３２は、算出された複数の相関係数のうち最大となる相関係数を、相関係数ｒ（θｎ）の算出値として、選択する。そして、推定演算部１３２は、本フローの処理を終了する。
推定演算部１３２は、算出された複数の相関係数のうち最大となる相関係数を、相関係数ｒ（θｎ）の算出値として、選択する。そして、推定演算部１３２は、本フローの処理を終了する。
以上のように、本実施形態に係る情報処理システム１における情報処理装置１１では、相関係数の計算において除外する所定の領域を変化させて、無線システムにおいて無線送信を行う無線局２１１のアンテナ指向性のパターン３１１を推定する。
したがって、情報処理装置１１では、第１実施形態あるいは第２実施形態よりも高精度で、推定を行うことができる。
したがって、情報処理装置１１では、第１実施形態あるいは第２実施形態よりも高精度で、推定を行うことができる。
（以上の実施形態の変形例）
＜変形例１＞
以上の実施形態では、円を用いて各種の領域を設定したが、円の代わりに、他の任意の図形が用いられてもよく、例えば、三角形、四角形、五角形、他の多角形、楕円、などが用いられてもよい。
また、以上の実施形態では、無線局２１１の推定位置を中心として図形（以上の実施形態では、円）を配置したが、他の配置が用いられてもよい。
また、以上の実施形態では、１つの図形または２つの図形を用いて所定の領域を設定したが、他の例として、３つ以上の図形を用いて所定の領域を定義してもよい。
＜変形例１＞
以上の実施形態では、円を用いて各種の領域を設定したが、円の代わりに、他の任意の図形が用いられてもよく、例えば、三角形、四角形、五角形、他の多角形、楕円、などが用いられてもよい。
また、以上の実施形態では、無線局２１１の推定位置を中心として図形（以上の実施形態では、円）を配置したが、他の配置が用いられてもよい。
また、以上の実施形態では、１つの図形または２つの図形を用いて所定の領域を設定したが、他の例として、３つ以上の図形を用いて所定の領域を定義してもよい。
＜変形例２＞
以上の実施形態では、アンテナ指向性のパターンとして、既知の１種類が想定されていた。これに対して、無線局２１１の登録情報などに複数のアンテナが登録されている場合などには、複数種類のアンテナ指向性のパターンが使用され得る。
そこで、推定演算部１３２は、複数種類のアンテナ指向性のパターンを順番に用いて以上の実施形態の場合と同様に各パターンについて相関係数の最大値を求め、複数種類のパターンのなかで相関係数の最大値が最大となるパターンを推定結果とする。そして、推定演算部１３２は、当該パターンにおいて、相関係数の最大値が得られた方向を推定結果とする。これにより、推定演算部１３２は、アンテナ指向性のパターンおよび方向の推定結果を得る。
以上の実施形態では、アンテナ指向性のパターンとして、既知の１種類が想定されていた。これに対して、無線局２１１の登録情報などに複数のアンテナが登録されている場合などには、複数種類のアンテナ指向性のパターンが使用され得る。
そこで、推定演算部１３２は、複数種類のアンテナ指向性のパターンを順番に用いて以上の実施形態の場合と同様に各パターンについて相関係数の最大値を求め、複数種類のパターンのなかで相関係数の最大値が最大となるパターンを推定結果とする。そして、推定演算部１３２は、当該パターンにおいて、相関係数の最大値が得られた方向を推定結果とする。これにより、推定演算部１３２は、アンテナ指向性のパターンおよび方向の推定結果を得る。
＜変形例３＞
以上の実施形態では、２次元の平面（２次元平面）を考慮したが、他の例として、推定演算部１３２は、アンテナ指向性のパターンの方向を推定する際に、３次元の空間（３次元空間）を考慮した演算を行う。
例えば、第１実施形態における［処理Ｔ４］では、推定演算部１３２は、仮定するアンテナ指向性のパターンの方向を変化させて相関係数を求める際、アンテナ指向性のパターンの方向を均等なステップで変化させる場合に、仮の方向を２次元の水平方向で変化させる。つまり、第１実施形態では、無線局２１１のアンテナと電波センサー２１－１１～２１－ＭＮとで設置高が同じであり、かつ、無線局２１１のアンテナが大地に対して垂直に設置されている場合を想定している。この場合、電波センサー２１－１１～２１－ＭＮは無線局２１１のアンテナから見て水平方向に存在するため、アンテナ指向性のパターンは、アンテナの水平面における指向性のパターンを用いることで十分である。
なお、アンテナ指向性のパターンの方向を推定する空間的な範囲と比べて、無線局２１１および無線局２１１のアンテナは、１点であるとみなしている。
以上の実施形態では、２次元の平面（２次元平面）を考慮したが、他の例として、推定演算部１３２は、アンテナ指向性のパターンの方向を推定する際に、３次元の空間（３次元空間）を考慮した演算を行う。
例えば、第１実施形態における［処理Ｔ４］では、推定演算部１３２は、仮定するアンテナ指向性のパターンの方向を変化させて相関係数を求める際、アンテナ指向性のパターンの方向を均等なステップで変化させる場合に、仮の方向を２次元の水平方向で変化させる。つまり、第１実施形態では、無線局２１１のアンテナと電波センサー２１－１１～２１－ＭＮとで設置高が同じであり、かつ、無線局２１１のアンテナが大地に対して垂直に設置されている場合を想定している。この場合、電波センサー２１－１１～２１－ＭＮは無線局２１１のアンテナから見て水平方向に存在するため、アンテナ指向性のパターンは、アンテナの水平面における指向性のパターンを用いることで十分である。
なお、アンテナ指向性のパターンの方向を推定する空間的な範囲と比べて、無線局２１１および無線局２１１のアンテナは、１点であるとみなしている。
これに対して、無線局２１１のアンテナと電波センサー２１－１１～２１－ＭＮとで設置高が異なる場合、あるいは、無線局２１１のアンテナが大地に対して傾けて設置される場合などに対応する。
例えば、推定演算部１３２は、第１実施形態における［処理Ｔ４］において、第１実施形態の場合の代わりに、アンテナ指向性の方向（水平方向）とアンテナの傾きを変化させて、相関係数を求める。
また、推定演算部１３２は、電波センサー２１－１１～２１－ＭＮでの受信電力を、例えば、無線局２１１の空中線電力（送信機出力電力）、無線局２１１のアンテナから見たその場所（電波センサー２１－１１～２１－ＭＮの場所）の指向方向（水平方向）およびアンテナから見た水平方向に対する垂直方向の角度から求まるアンテナ利得、伝搬損失、から演算により求める。また、当該受信電力の演算において、例えば、ビルなどの障害物の情報が用いられてもよい。なお、これらの値は、それぞれ、例えば、既知の手法によって求められてもよく、あるいは、ユーザなどによって設定されてもよい。
無線局２１１のアンテナから電波センサー２１－１１を見た指向方向（水平方向）と、無線局２１１のアンテナから電波センサー２１－１１を見た水平方向に対する垂直方向の角度に基づいて、無線局２１１のアンテナの利得は垂直方向の利得と水平方向の利得との積（ｄＢの単位では、値の加算）となる。
例えば、推定演算部１３２は、第１実施形態における［処理Ｔ４］において、第１実施形態の場合の代わりに、アンテナ指向性の方向（水平方向）とアンテナの傾きを変化させて、相関係数を求める。
また、推定演算部１３２は、電波センサー２１－１１～２１－ＭＮでの受信電力を、例えば、無線局２１１の空中線電力（送信機出力電力）、無線局２１１のアンテナから見たその場所（電波センサー２１－１１～２１－ＭＮの場所）の指向方向（水平方向）およびアンテナから見た水平方向に対する垂直方向の角度から求まるアンテナ利得、伝搬損失、から演算により求める。また、当該受信電力の演算において、例えば、ビルなどの障害物の情報が用いられてもよい。なお、これらの値は、それぞれ、例えば、既知の手法によって求められてもよく、あるいは、ユーザなどによって設定されてもよい。
無線局２１１のアンテナから電波センサー２１－１１を見た指向方向（水平方向）と、無線局２１１のアンテナから電波センサー２１－１１を見た水平方向に対する垂直方向の角度に基づいて、無線局２１１のアンテナの利得は垂直方向の利得と水平方向の利得との積（ｄＢの単位では、値の加算）となる。
具体例として、２次元の平面における半径ｍの円の内側、外側の代わりに、３次元の空間における半径ｍの球の内側、外側が用いられてもよい。
また、２次元の平面における半径ｍ１の円の内側、外側、および、半径ｍ２の円の内側、外側の代わりに、３次元の空間における半径ｍ１の球の内側、外側、および、半径ｍ２の球の内側、外側が用いられてもよい。
また、２次元の平面における半径ｍ１の円の内側、外側、および、半径ｍ２の円の内側、外側の代わりに、３次元の空間における半径ｍ１の球の内側、外側、および、半径ｍ２の球の内側、外側が用いられてもよい。
＜変形例４＞
以上の実施形態では、１つの無線局２１１を考慮したが、複数の無線局を考慮してもよい。
例えば、推定演算部１３２は、［処理Ｔ１］～［処理Ｔ５］を行うことで得られた１局目の無線局２１１のアンテナ指向性のパターンの方向の推定結果を用いて、１局目の無線局２１１のアンテナ指向性のパターンによる電波センサー２１－１１～２１－ＭＮでの受信電力を推定する。
そして、推定演算部１３２は、無線局２１１の推定結果の電力分布を、電波センサー２１－１１～２１－ＭＮで取得された電力分布から差し引いて、その結果として空間的な電力分布を求める。推定演算部１３２は、このようにして求められた空間的な電力分布を、仮想的に、電波センサー２１－１１～２１－ＭＮによる検出結果から得られる電力分布とみなして、２局目の無線局２１１についてアンテナ指向性のパターンの方向を推定する。
推定演算部１３２は、３局目以降の無線局についても同様に、直前までの無線局のアンテナによる電波センサー２１－１１～２１－ＭＮでの推定受信電力を、電波センサー２１－１１～２１－ＭＮで取得された受信電力から差し引いて得られる空間的な電力分布を用いて、アンテナ指向性のパターンの方向の推定を行う。
このような一連の処理を、すべての無線局に対して順番に行うことで、各無線局のアンテナ指向性のパターンの方向を推定することができる。
以上の実施形態では、１つの無線局２１１を考慮したが、複数の無線局を考慮してもよい。
例えば、推定演算部１３２は、［処理Ｔ１］～［処理Ｔ５］を行うことで得られた１局目の無線局２１１のアンテナ指向性のパターンの方向の推定結果を用いて、１局目の無線局２１１のアンテナ指向性のパターンによる電波センサー２１－１１～２１－ＭＮでの受信電力を推定する。
そして、推定演算部１３２は、無線局２１１の推定結果の電力分布を、電波センサー２１－１１～２１－ＭＮで取得された電力分布から差し引いて、その結果として空間的な電力分布を求める。推定演算部１３２は、このようにして求められた空間的な電力分布を、仮想的に、電波センサー２１－１１～２１－ＭＮによる検出結果から得られる電力分布とみなして、２局目の無線局２１１についてアンテナ指向性のパターンの方向を推定する。
推定演算部１３２は、３局目以降の無線局についても同様に、直前までの無線局のアンテナによる電波センサー２１－１１～２１－ＭＮでの推定受信電力を、電波センサー２１－１１～２１－ＭＮで取得された受信電力から差し引いて得られる空間的な電力分布を用いて、アンテナ指向性のパターンの方向の推定を行う。
このような一連の処理を、すべての無線局に対して順番に行うことで、各無線局のアンテナ指向性のパターンの方向を推定することができる。
一構成例として、情報処理装置１１では、推定演算部１３２は、電力の測定結果に基づく空間的な電力分布から対象無線局以外の無線局のアンテナ指向性のパターンについて推定された電力分布を差し引いた結果の電力分布と、前記対象無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定する。
＜変形例５＞
以上の実施形態では、１つのアンテナ指向性のパターンごとについて推定を行ったが、他の例として、複数のアンテナ指向性のパターンをまとめて推定を行ってもよい。
例えば、推定演算部１３２は、無線局２１１の複数のアンテナ指向性のパターンのそれぞれを仮の方向に配置した場合について、電波センサー２１－１１～２１－ＭＮによる検出結果（空間的な電力分布）を推定する。この推定結果は、これら複数のアンテナ指向性のパターンによる電力分布の合成結果に相当する。
この場合、例えば、推定演算部１３２は、複数のアンテナ指向性のパターンのそれぞれについて、別々または同時に、仮の方向を変化させて、推定を行う。
なお、無線局２１１で使用されるアンテナの本数が既知である場合には既知の本数が用いられればよいが、他の例として、無線局２１１で使用されるアンテナの本数が不明である場合には、当該本数が推定されてもよい。
ここで、１つの無線局における複数のアンテナ指向性のパターンをまとめて推定を行う代わりに、複数の無線局のそれぞれにおけるアンテナ指向性のパターンをまとめて推定を行ってもよい。この場合、それぞれの無線局におけるアンテナ指向性のパターンの数は、１以上の任意の値であってもよい。
以上の実施形態では、１つのアンテナ指向性のパターンごとについて推定を行ったが、他の例として、複数のアンテナ指向性のパターンをまとめて推定を行ってもよい。
例えば、推定演算部１３２は、無線局２１１の複数のアンテナ指向性のパターンのそれぞれを仮の方向に配置した場合について、電波センサー２１－１１～２１－ＭＮによる検出結果（空間的な電力分布）を推定する。この推定結果は、これら複数のアンテナ指向性のパターンによる電力分布の合成結果に相当する。
この場合、例えば、推定演算部１３２は、複数のアンテナ指向性のパターンのそれぞれについて、別々または同時に、仮の方向を変化させて、推定を行う。
なお、無線局２１１で使用されるアンテナの本数が既知である場合には既知の本数が用いられればよいが、他の例として、無線局２１１で使用されるアンテナの本数が不明である場合には、当該本数が推定されてもよい。
ここで、１つの無線局における複数のアンテナ指向性のパターンをまとめて推定を行う代わりに、複数の無線局のそれぞれにおけるアンテナ指向性のパターンをまとめて推定を行ってもよい。この場合、それぞれの無線局におけるアンテナ指向性のパターンの数は、１以上の任意の値であってもよい。
一構成例として情報処理装置１１では、推定演算部１３２は、電力の測定結果に基づく空間的な電力分布と、複数のアンテナ指向性のパターンのそれぞれの方向をそれぞれの仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向の組み合わせのそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向の組み合わせを前記複数の前記アンテナ指向性のパターンの方向の推定結果と判定する。
＜構成例＞
一構成例として、情報処理装置１１では、推定演算部１３２によって、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定する。
一構成例として、情報処理方法（以上の実施形態では、情報処理装置１１において行われる方法）では、推定演算部１３２によって、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定する。
一構成例として、プログラム（以上の実施形態では、情報処理装置１１においてプロセッサーにより実行されるプログラム）は、推定演算部１３２によって、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定するステップを、コンピュータ（以上の実施形態では、情報処理装置１１を構成するコンピュータ）に実行させるためのプログラムである。
一構成例として、情報処理装置１１では、推定演算部１３２によって、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定する。
一構成例として、情報処理方法（以上の実施形態では、情報処理装置１１において行われる方法）では、推定演算部１３２によって、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定する。
一構成例として、プログラム（以上の実施形態では、情報処理装置１１においてプロセッサーにより実行されるプログラム）は、推定演算部１３２によって、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定するステップを、コンピュータ（以上の実施形態では、情報処理装置１１を構成するコンピュータ）に実行させるためのプログラムである。
なお、以上に説明した情報処理装置１１などの任意の装置における任意の構成部の機能を実現するためのプログラムを、コンピュータ読み取り可能な記録媒体に記録し、そのプログラムをコンピュータシステムに読み込ませて実行するようにしてもよい。なお、ここでいう「コンピュータシステム」とは、オペレーティングシステム（ＯＳ：Ｏｐｅｒａｔｉｎｇ Ｓｙｓｔｅｍ）あるいは周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ（Ｃｏｍｐａｃｔ Ｄｉｓｃ）－ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークあるいは電話回線等の通信回線を介してプログラムが送信された場合のサーバーやクライアントとなるコンピュータシステム内部の揮発性メモリ（ＲＡＭ）のように、一定時間プログラムを保持しているものも含むものとする。
また、上記のプログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワークあるいは電話回線等の通信回線のように情報を伝送する機能を有する媒体のことをいう。
また、上記のプログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、上記のプログラムは、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイルであってもよい。差分ファイルは、差分プログラムと呼ばれてもよい。
また、上記のプログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、上記のプログラムは、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイルであってもよい。差分ファイルは、差分プログラムと呼ばれてもよい。
以上に説明した任意の装置における任意の構成部の機能は、プロセッサーにより実現されてもよい。例えば、実施形態における各処理は、プログラム等の情報に基づき動作するプロセッサーと、プログラム等の情報を記憶するコンピュータ読み取り可能な記録媒体により実現されてもよい。ここで、プロセッサーは、例えば、各部の機能が個別のハードウェアで実現されてもよく、あるいは、各部の機能が一体のハードウェアで実現されてもよい。例えば、プロセッサーはハードウェアを含み、当該ハードウェアは、デジタル信号を処理する回路およびアナログ信号を処理する回路のうちの少なくとも一方を含んでもよい。例えば、プロセッサーは、回路基板に実装された１または複数の回路装置、あるいは、１または複数の回路素子のうちの一方または両方を用いて、構成されてもよい。回路装置としてはＩＣ（Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）などが用いられてもよく、回路素子としては抵抗あるいはキャパシターなどが用いられてもよい。
ここで、プロセッサーは、例えば、ＣＰＵであってもよい。ただし、プロセッサーは、ＣＰＵに限定されるものではなく、例えば、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、あるいは、ＤＳＰ（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒ）等のような、各種のプロセッサーが用いられてもよい。また、プロセッサーは、例えば、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）によるハードウェア回路であってもよい。また、プロセッサーは、例えば、複数のＣＰＵにより構成されていてもよく、あるいは、複数のＡＳＩＣによるハードウェア回路により構成されていてもよい。また、プロセッサーは、例えば、複数のＣＰＵと、複数のＡＳＩＣによるハードウェア回路と、の組み合わせにより構成されていてもよい。また、プロセッサーは、例えば、アナログ信号を処理するアンプ回路あるいはフィルター回路等のうちの１以上を含んでもよい。
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
１…情報処理システム、１１…情報処理装置、２１－１１～２１－ＭＮ…電波センサー、３１…ネットワーク、１１１…入力部、１１２…出力部、１１３…通信部、１１４…記憶部、１１５…制御部、１３１…取得部、１３２…推定演算部、１３３…出力制御部、２１１…無線局、３１１…パターン、３３１…中心方向、３５１…角度幅、４１１、４２１、５１１～５１２、５２１～５２２…円領域
Claims (12)
- 所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定する推定演算部を備える、
情報処理装置。 - 前記推定演算部は、複数の前記パターンのそれぞれについて前記相関係数を演算し、最大の相関係数が得られた前記パターンおよび前記仮の方向を前記アンテナ指向性のパターンおよび方向の推定結果と判定する、
請求項１に記載の情報処理装置。 - 前記推定演算部は、２次元平面における前記パターンの情報を用いて前記相関係数を演算する、
請求項１または請求項２に記載の情報処理装置。 - 前記推定演算部は、３次元空間における前記パターンの情報を用いて前記相関係数を演算する、
請求項１または請求項２に記載の情報処理装置。 - 前記推定演算部は、電力の前記測定結果に基づく空間的な前記電力分布と、前記無線局の前記アンテナ指向性の前記パターンの方向を前記仮の方向としたときに想定される空間的な前記電力分布との一方または両方について所定の重み付けを行い、前記重み付けが行われた結果に基づいて前記相関係数を演算する、
請求項１から請求項４のいずれか１項に記載の情報処理装置。 - 前記所定の領域は、所定の円の内側の領域を含む、
請求項１から請求項５のいずれか１項に記載の情報処理装置。 - 前記所定の領域は、さらに、前記所定の円よりも半径が大きい所定の他の円の外側の領域を含む、
請求項６に記載の情報処理装置。 - 前記推定演算部は、前記所定の領域を変化させて前記相関係数を演算し、最大の相関係数が得られた前記パターンおよび前記仮の方向を前記アンテナ指向性のパターンおよび方向の推定結果と判定する、
請求項１から請求項７のいずれか１項に記載の情報処理装置。 - 前記推定演算部は、前記無線局の推定位置を中心として設定される前記所定の領域を用いる、
請求項１から請求項８のいずれか１項に記載の情報処理装置。 - 前記推定演算部は、前記無線局の位置推定精度に応じた前記所定の領域を用いる、
請求項１から請求項８のいずれか１項に記載の情報処理装置。 - 推定演算部が、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定する、
情報処理方法。 - 推定演算部が、所定の領域を除外して、電力の測定結果に基づく空間的な電力分布と、無線局のアンテナ指向性のパターンの方向を仮の方向としたときに想定される空間的な電力分布との相関係数を演算し、複数の前記仮の方向のそれぞれについて得られた前記相関係数のうちで最大の相関係数が得られた前記仮の方向を前記アンテナ指向性のパターンの方向の推定結果と判定するステップを、
コンピュータに実行させるためのプログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020111794A JP7257992B2 (ja) | 2020-06-29 | 2020-06-29 | 情報処理装置、情報処理方法およびプログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020111794A JP7257992B2 (ja) | 2020-06-29 | 2020-06-29 | 情報処理装置、情報処理方法およびプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2022010966A JP2022010966A (ja) | 2022-01-17 |
| JP7257992B2 true JP7257992B2 (ja) | 2023-04-14 |
Family
ID=80147854
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020111794A Active JP7257992B2 (ja) | 2020-06-29 | 2020-06-29 | 情報処理装置、情報処理方法およびプログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7257992B2 (ja) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005195532A (ja) | 2004-01-09 | 2005-07-21 | National Institute Of Information & Communication Technology | 電子機器から漏洩する電波の波源位置の推定装置 |
| JP2007129697A (ja) | 2005-10-05 | 2007-05-24 | Matsushita Electric Ind Co Ltd | 無線通信装置 |
| JP2013214896A (ja) | 2012-04-03 | 2013-10-17 | Nippon Telegr & Teleph Corp <Ntt> | 基地局装置、無線通信方法、及び無線通信システム |
| JP2017208679A (ja) | 2016-05-18 | 2017-11-24 | 日本電気株式会社 | 受信電力推定装置、受信電力推定方法および受信電力推定プログラム |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH088814A (ja) * | 1994-06-23 | 1996-01-12 | Nippon Motorola Ltd | 移動無線通信方式 |
| JPH10232252A (ja) * | 1997-02-20 | 1998-09-02 | Advantest Corp | ３次元波源分布のホログラム観測方法 |
| JP6510358B2 (ja) * | 2015-08-07 | 2019-05-08 | 日本電信電話株式会社 | 無線通信システム及び無線通信方法 |
-
2020
- 2020-06-29 JP JP2020111794A patent/JP7257992B2/ja active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005195532A (ja) | 2004-01-09 | 2005-07-21 | National Institute Of Information & Communication Technology | 電子機器から漏洩する電波の波源位置の推定装置 |
| JP2007129697A (ja) | 2005-10-05 | 2007-05-24 | Matsushita Electric Ind Co Ltd | 無線通信装置 |
| JP2013214896A (ja) | 2012-04-03 | 2013-10-17 | Nippon Telegr & Teleph Corp <Ntt> | 基地局装置、無線通信方法、及び無線通信システム |
| JP2017208679A (ja) | 2016-05-18 | 2017-11-24 | 日本電気株式会社 | 受信電力推定装置、受信電力推定方法および受信電力推定プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2022010966A (ja) | 2022-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | CSI-based fingerprinting for indoor localization: A deep learning approach | |
| JP5741223B2 (ja) | 情報処理装置、補正方法、及び補正プログラム | |
| US8102314B2 (en) | Method and apparatus for determining the location of a mobile object | |
| KR101325579B1 (ko) | 이동 단말 장치의 방향을 추정하는 장치 및 방법 | |
| Okşar | A Bluetooth signal strength based indoor localization method | |
| US20150338502A1 (en) | Determining a position of a terminal | |
| KR101163335B1 (ko) | 수신신호강도의 거리 추정 방식에 의거한 실내 무선 측위 방법 및 이 방법을 실시하기 위한 프로그램이 기록된 기록매체 | |
| WO2003102622A1 (en) | Probabilistic model for a positioning technique | |
| JP2007013500A (ja) | 無線端末位置推定システム、無線端末位置推定システムの位置推定方法及びデータ処理装置 | |
| US9660740B2 (en) | Signal strength distribution establishing method and wireless positioning system | |
| KR20120115895A (ko) | 수신신호강도의 거리추정방식에 의거하여 4개 이상의 앵커노드를 이용한 실내 무선 측위 방법 및 이 방법을 실시하기 위한 프로그램이 기록된 기록매체 | |
| US20160091341A1 (en) | Method and apparatus for object localization | |
| Wang et al. | A near-optimal least squares solution to received signal strength difference based geolocation | |
| KR100863917B1 (ko) | 하나 이상의 기지국에 대한 무선 이동 장치의 위치 추정 | |
| JP7257992B2 (ja) | 情報処理装置、情報処理方法およびプログラム | |
| JP2016520801A (ja) | 位置判定を改善する送信機の方向性プルーニング | |
| Cimdins et al. | Sundew: Design and evaluation of a model-based device-free localization system | |
| US20160091349A1 (en) | Method and apparatus for object localizing | |
| US20210239785A1 (en) | Interference source searching method and interference source searching apparatus | |
| JP6331072B2 (ja) | ホワイトスペース検出装置、ホワイトスペース検出方法、及びプログラム | |
| KR101694521B1 (ko) | 전파지문지도 생성 장치 및 방법 | |
| CN115840192A (zh) | 一种基于空间估计谱置信度估计的室内定位方法 | |
| JP7325358B2 (ja) | 情報処理装置、情報処理方法およびプログラム | |
| US20190219664A1 (en) | System and method for range and angle determination to an array of radio receivers | |
| JP2021135250A (ja) | 情報処理装置、情報処理方法およびプログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220531 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20230214 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20230307 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230404 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7257992 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |