(57)【要約】【課題】信号の受信装置近傍に信号発信源が存在する場合においても高精度に発信源の位置を推定したい。【解決手段】受信センサ１０ａ〜ｍは、所定の方向から到来する信号を受信する。信号変換部１２ａ〜ｍは、受信信号に対してアナログ／デジタル変換を行う。相関行列計算部１４は、受信信号の相関行列を計算する。積分部１６は、相関行列を積分する。固有値計算部１８は、相関行列から固有値等を計算する。関数処理部２０は、座標をパラメータとする評価関数より信号発信源の位置を推定する。強度判断部２２は、関数処理部２０によって推定された信号発信源より到来する信号の強度を計算し、その値がしきい値より小さい場合、その信号に対応する欠陥位置を欠陥位置の候補から除外する。表示部２４は、以上の推定結果を所定の形式で表示する。【選択図】　　　　図３
2004144731.tif 000002
【特許請求の範囲】
【請求項１】
受信信号の相関行列を固有値解析するステップと、
前記固有値解析の結果にもとづいて、座標をパラメータとする評価関数より信号発信源の位置を推定するステップと、
を含むことを特徴とする位置推定方法。
【請求項２】
前記評価関数は、前記パラメータの座標と受信信号の入力位置における座標のちがいをもとに定まる関数で与えられるステアリングベクトルを含むことを特徴とする請求項１に記載の位置推定方法。
【請求項３】
前記評価関数を前記パラメータの座標で掃引し、その結果から特徴点を検出し、前記信号発信源の位置を推定することを特徴とする請求項１または２のいずれかに記載の位置推定方法。
【請求項４】
位置が特定された信号発信源からの信号を受信するステップと、
前記受信した信号の相関行列を固有値解析するステップと、
前記固有値解析の結果と前記信号発生源の位置にもとづいて、座標をパラメータとする評価関数より前記信号の受信位置を推定するステップと、
を含むことを特徴とする位置推定方法。
【請求項５】
前記評価関数は、前記信号発信源の位置における座標と前記パラメータの座標のちがいをもとに定まる関数で与えられるステアリングベクトルを含むことを特徴とする請求項４に記載の位置推定方法。
【請求項６】
前記評価関数を前記パラメータの座標で掃引し、その結果から特徴点を検出し、前記信号の受信位置を推定することを特徴とする請求項４または５のいずれかに記載の位置推定方法。
【請求項７】
推定対象であり、かつパラメータとされた座標と既知の座標とのちがいをもとに定まる関数で与えられるステアリングベクトルを含んだ評価関数について、前記パラメータとされた座標の掃引による結果から特徴点を検出し、前記推定対象の座標を導出することを特徴とする位置推定方法。
【請求項８】
位置が既知の信号発信源からの信号を既知の位置で受信するステップと、
前記受信した信号の相関行列を固有値解析して、部分雑音空間の情報を導出するステップと、
前記導出した部分雑音空間の情報が、前記信号発信源の位置と前記信号を受信した位置の座標のちがいをもとに定まる関数で与えられるステアリングベクトルに略直交するように、補正値を計算するステップと、
前記受信した信号を前記補正値によって補正するステップと、
を含むことを特徴とするキャリブレーション方法。
【請求項９】
受信信号が入力される受信部と、
前記受信信号の相関行列を計算する相関行列計算部と、
前記相関行列の固有値と固有ベクトルを計算する固有値計算部と、
前記固有値と前記固有ベクトルにもとづいて、座標をパラメータとする評価関数により信号発信源の位置を推定する関数処理部と、
前記推定結果を出力する出力部と、
を含むことを特徴とする位置推定装置。
【請求項１０】
前記関数処理部における前記評価関数は、前記パラメータの座標と受信信号の入力位置における座標のちがいをもとに定まる関数で与えられるステアリングベクトルを含むことを特徴とする請求項９に記載の位置推定装置。
【請求項１１】
前記関数処理部は、前記評価関数を前記パラメータの座標で掃引し、その結果から特徴点を検出し、前記信号発信源の位置を推定することを特徴とする請求項９または１０のいずれかに記載の位置推定装置。
【請求項１２】
信号送信中において固定された周波数を有する信号を出力する送信部をさらに含むことを特徴とする請求項９から１１のいずれかに記載の位置推定装置。
【請求項１３】
推定された信号発生源からの受信信号の強度が、所定のしきい値を超えない場合に、前記受信信号に対する信号発生源を有意な検出対象から除外する処理を実行する強度判断部をさらに含むことを特徴とする請求項９から１２のいずれかに記載の位置推定装置。
【請求項１４】
位置が特定された信号発信源からの信号を受信する受信部と、
前記受信した信号の相関行列を計算する相関行列計算部と、
前記相関行列の固有値と固有ベクトルを計算する固有値計算部と、
前記固有値、前記固有ベクトルと前記信号発信源の位置にもとづいて、座標をパラメータとする評価関数により前記信号の受信位置を推定する関数処理部と、
前記推定結果を出力する出力部と、
を含むことを特徴とする位置推定装置。
【請求項１５】
前記関数処理部における前記評価関数は、前記信号発信源の位置における座標と前記パラメータの座標のちがいをもとに定まる関数で与えられるステアリングベクトルを含むことを特徴とする請求項１４に記載の位置推定装置。
【請求項１６】
前記関数処理部は、前記評価関数を前記パラメータの座標で掃引し、その結果から特徴点を検出し、前記信号の受信位置を推定することを特徴とする請求項１４または１５のいずれかに記載の位置推定装置。
【請求項１７】
位置が既知の信号発信源からの信号を既知の位置で受信する受信部と、
前記受信した信号の相関行列を固有値解析して、部分雑音空間の情報を導出する部分雑音空間計算部と、
前記導出した部分雑音空間の情報が、前記信号発信源の位置と前記信号を受信した位置の座標のちがいをもとに定まる関数で与えられるステアリングベクトルに略直交するように、補正値を計算する補正値計算部と、
前記受信した信号を前記補正値によって補正する補正部と、
を含むことを特徴とするキャリブレーション装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、信号の発信源の位置を受信した信号より推定する技術に関する。特に、受信装置近傍に信号発信源が存在する場合における発信源の位置推定技術に関する。また、材料内部の欠陥位置を受信した信号より推定する技術に関する。また、信号を受信した位置を推定する技術に関する。また、受信した信号を補正するキャリブレーション技術に関する。
【０００２】
【従来の技術】
一般に、各種構造用材料の信頼性評価、並びに確保を目的として、非破壊検査のひとつである超音波探傷試験が行われている。超音波探傷試験は、材料を壊すことなく、また人体にも無害である利点を有する。従来の超音波探傷試験における基本的な欠陥位置推定原理は以下の通りである。図１は、従来の超音波探傷試験装置を示す。被検査材料の内部に超音波信号を伝播させるために、材料表面に超音波センサである送受信センサ５０を取り付け、その送受信センサ５０から超音波パルス信号を送信する。材料内部に欠陥が存在する場合、欠陥によって超音波パルス信号が反射される。その第一反射信号は材料表面に取り付けた複数の送受信センサ５０、受信センサ５２によって電力信号として受信される。送受信センサ５０から送信された超音波パルス信号が送受信センサ５０と受信センサ５２に受信される間の往復伝播時間差により、欠陥の位置が推定される。
【０００３】
一方、受信信号より信号の到来方向を推定する高分解能位置推定技術に関する研究開発が行われている。この技術は、移動体通信やレーダのように電波伝搬環境下において、信号の到来方向や位置を推定する場合や、生活環境下における震音発生源の方向や位置を推定する場合に適用される。高分解能位置推定技術には、例えば、受信信号の相関行列を計算し、その固有値解析を行うことを基本原理した到来方向推定方法がある。到来方向推定方法のひとつはＭＵＳＩＣ（ＭＵｌｔｉｐｌｅ　ＳＩｇｎａｌ　Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）アルゴリズムである。
【０００４】
アレーアンテナにより信号の到来方向を推定する場合、ＭＵＳＩＣアルゴリズムは以下の通り適用される。図２は、１次元リニアアレイにおけるアレーアンテナ装置１２０の構成を示す。これはＭ素子のアンテナ素子６０ａ〜ｍと無線装置６２で構成される。到来する信号数をＫとすると、時刻ｔにおいてアンテナ素子６０ａ〜ｍが受信する信号をベクトルＸ（ｔ）で表すと、
【０００５】
【数１】
2004144731.tif 000003
ここで、Ｎ（ｔ）は時刻ｔにおける受信信号に含まれる白色雑音ベクトルであり、Ｓ（ｔ）およびＡは、以下の通り与えられる。
【０００６】
【数２】
2004144731.tif 000004
【数３】
2004144731.tif 000005
ここで、添え字Ｔは転置、ｓｋ（ｔ）は時刻ｔにおける目標ｋからの信号の振幅を示す。ｋは１からＫとする。ａ（θｋ）は目標ｋの到来方向を表すモードベクトルであり、Ｍ成分により構成される。
【０００７】
受信信号Ｘの相関行列をＲとすると、
【数４】
2004144731.tif 000006
である。添え字Ｈは複素共役転置を示す。相関行列Ｒを固有値解析し、その固有値をλ、固有ベクトルをｅｉとする。ここでｉは１からＭとする。到来する信号間に相関がないとすると、固有値は以下の関係を有する。
【０００８】
【数５】
2004144731.tif 000007
上式において、最小固有値に対応した雑音部分空間をＥｉとすると、次式によってＭＵＳＩＣアルゴリズムの評価関数は与えられる。
【０００９】
【数６】
2004144731.tif 000008
ここで、分母はノルムを表す。また、ａ（θ）はステアリングベクトルであり、ここでは次式で表される。
【００１０】
【数７】
2004144731.tif 000009
ＭＵＳＩＣアルゴリズムは、ステアリングベクトルの角度θを掃引し、θに対するＭＵＳＩＣ評価関数のピークを探す。ピークに対応するθが、信号の到来方向を示す。信号の到来方向に加えて信号の発信源までの距離も推定する場合、さらに使用する信号の周波数を掃引する必要がある。
【００１１】
【非特許文献１】
菊間信良著，「アレーアンテナによる適応信号処理」，科学技術出版，１９９８年，ｐ．１９５−１９９
【００１２】
【発明が解決しようとする課題】
本発明者はこうした状況下、以下の課題を認識するに至った。従来の欠陥位置推定原理を使用した超音波探傷試験において、受信センサの受信信号には、第一反射信号のみならず、材料表面での反射信号やセンサ設置面における表面波等も含まれる。これらの信号が重畳受信されるため、正確な第一反射信号の往復伝播時間差を得ることが困難である。
【００１３】
そこで本発明者は、超音波探傷試験における欠陥位置の推定精度を高めるため、従来移動体通信やレーダの分野等で信号の到来方向を高精度に推定する到来方向推定方法を超音波探傷試験へ使用することに到達した。しかし、これには以下に述べる課題がある。
【００１４】
従来の到来方向推定方法におけるＭＵＳＩＣアルゴリズムは、一般に、信号が無限遠方から到来するものとし、各受信センサに到来する信号の入射角は等しいと仮定している。そのような仮定は、従来の電波伝搬環境下においては十分成立すると考えられている。超音波探傷試験による非破壊検査においては、被検査材料の大きさは有限であるため、従来のＭＵＳＩＣアルゴリズムにおける仮定が成立せず、推定精度が低下する。さらに、電波伝搬環境下においても、受信センサの近傍より送信された信号が到来する場合に同様である。また、超音波探傷試験における超音波周波数は被検査材料の特性により定められる。そのため、信号発信源までの距離を推定する場合においても周波数を掃引できない。
【００１５】
本発明はこうした状況に鑑みなされたものであり、その目的は信号発信源が受信センサの近傍に存在する場合においても、信号発信源の位置を推定する技術の提供である。特に、信号発信源が受信センサの近傍に存在する場合においても、信号発信源の位置を推定するＭＵＳＩＣアルゴリズムの提供である。また、ＭＵＳＩＣアルゴリズムによる信号発信源までの距離推定において、周波数を掃引不要とする技術の提供である。さらに、超音波探傷試験に高精度化をもたらす技術の提供である。また、受信した信号から、当該信号を受信した位置を推定する技術の提供である。また、受信した信号から、当該信号を補正するためのキャリブレーション技術の提供である。
【００１６】
【課題を解決するための手段】
本発明においては、超音波探傷試験に対して固有値解析に基づく到来方向推定方法を使用して、欠陥位置の推定精度を向上させる。
本発明のある態様は、位置推定方法に関する。この方法は、受信信号の相関行列を固有値解析するステップと、固有値解析の結果にもとづいて、座標をパラメータとする評価関数より信号発信源の位置を推定するステップとを含む。
【００１７】
「固有値解析」とは、行列の固有値、固有ベクトルを計算すること、およびそれらに所定の処理を施すことを含む。
「座標」は、信号発信源の位置を特定するために使用される。その範囲では、直交座標系以外の座標系も含まれるものとする。
この方法によれば、評価関数のパラメータに信号発信源の位置を特定する座標が含まれるため、座標自体の特定が可能になる。また、各受信センサに到来する信号の入射角が等しいと仮定されていないため、信号発信源が受信センサの近傍に存在する場合でも、その位置を推定することが可能である。
【００１８】
評価関数は、パラメータの座標と受信信号の入力位置における座標のちがいをもとに定まる関数で与えられるステアリングベクトルを含んでもよい。また、評価関数をパラメータの座標で掃引し、その結果から特徴点を検出し、信号発信源の位置を推定してもよい。
【００１９】
「パラメータの座標と受信信号の入力位置における座標のちがいをもとに定まる関数」とは、これらの座標間の距離から定まる関数であってもよいし、その２乗などから定まる関数でもよい。
「特徴点」とは、例えば、極大値や極小値のように、何らかの定義の下で他の点と区別しうる値をいう。
【００２０】
本発明の別の態様も、位置推定方法に関する。この方法は、位置が特定された信号発信源からの信号を受信するステップと、受信した信号の相関行列を固有値解析するステップと、固有値解析の結果と信号発生源の位置にもとづいて、座標をパラメータとする評価関数より前記信号の受信位置を推定するステップとを含む。
評価関数は、信号発信源の位置における座標と前記パラメータの座標のちがいをもとに定まる関数で与えられるステアリングベクトルを含んでもよい。評価関数をパラメータの座標で掃引し、その結果から特徴点を検出し、信号の受信位置を推定してもよい。
【００２１】
本発明のさらに別の態様も、位置推定方法に関する。この方法は、推定対象であり、かつパラメータとされた座標と既知の座標とのちがいをもとに定まる関数で与えられるステアリングベクトルを含んだ評価関数について、パラメータとされた座標の掃引による結果から特徴点を検出し、推定対象の座標を導出する。
【００２２】
本発明のさらに別の態様は、キャリブレーション方法に関する。この方法は、位置が既知の信号発信源からの信号を既知の位置で受信するステップと、受信した信号の相関行列を固有値解析して、部分雑音空間の情報を導出するステップと、導出した部分雑音空間の情報が、信号発信源の位置と信号を受信した位置の座標のちがいをもとに定まる関数で与えられるステアリングベクトルに略直交するように、補正値を計算するステップと、受信した信号を前記補正値によって補正するステップとを含む。
【００２３】
本発明のさらに別の態様は、位置推定装置に関する。この装置は、受信信号が入力される受信部と、受信信号の相関行列を計算する相関行列計算部と、相関行列の固有値と固有ベクトルを計算する固有値計算部と、固有値と固有ベクトルにもとづいて、座標をパラメータとする評価関数により信号発信源の位置を推定する関数処理部と、推定結果を出力する出力部とを含む。
信号送信中において固定された周波数を有する信号を出力する送信部をさらに含んでもよい。
「信号送信中」は、信号を送信している期間を意味する。これには、信号の送信開始から送信終了までの期間だけではなく、その中の所定の期間も含むものとする。
【００２４】
本発明のさらに別の態様は、位置推定装置に関する。この装置は、位置が特定された信号発信源からの信号を受信する受信部と、受信した信号の相関行列を計算する相関行列計算部と、相関行列の固有値と固有ベクトルを計算する固有値計算部と、固有値、固有ベクトルと信号発信源の位置にもとづいて、座標をパラメータとする評価関数により信号の受信位置を推定する関数処理部と、推定結果を出力する出力部とを含む。
【００２５】
本発明のさらに別の態様は、キャリブレーション装置に関する。この装置は、位置が既知の信号発信源からの信号を既知の位置で受信する受信部と、受信した信号の相関行列を固有値解析して、部分雑音空間の情報を導出する部分雑音空間計算部と、導出した部分雑音空間の情報が、信号発信源の位置と信号を受信した位置の座標のちがいをもとに定まる関数で与えられるステアリングベクトルに略直交するように、補正値を計算する補正値計算部と、受信した信号を補正値によって補正する補正部とを含む。
なお、以上の構成要素の任意の組合せや組み替え、本発明を方法、コンピュータプログラムなどと表現したものもまた、本発明の態様として有効である。
【００２６】
【発明の実施の形態】
（実施の形態１）
以下、本発明の実施の形態を図面に基づいて説明する。図３は、実施の形態に係る位置推定装置１００の構成図である。位置推定装置１００は、受信センサ１０ａ〜ｍ、信号変換部１２ａ〜ｍ、相関行列計算部１４、積分部１６、固有値計算部１８、関数処理部２０、強度判断部２２、表示部２４を有する。
【００２７】
受信センサ１０ａ〜ｍは、所定の方向から到来する信号を受信する受信部である。信号は、電磁波、音波、超音波、水中音波などである。信号の種類に応じて、受信センサ１０ａ〜ｍの種類は異なる。受信センサ１０ａ〜ｍの種類には、超音波探傷試験用素子やアンテナ素子などがある。また、位置推定装置１００が超音波探傷試験装置の場合、受信センサ１０ａ〜ｍのいずれかが送信センサをかね、ここから信号が送信されても構わない。ここで、受信センサ１０ａ〜ｍはＭ個あるものとする。
【００２８】
信号変換部１２ａ〜ｍは、受信センサ１０ａ〜ｍで受信された信号に対して所定の変換を行う。所定の変更は、例えば、周波数変換、復調、アナログ／デジタル変換である。また、信号変換部１２ａ〜ｍは、一体で構成され、受信センサ１０ａ〜ｍで受信された信号をシリアルに処理しても構わない。
【００２９】
相関行列計算部１４は、受信信号の相関行列を計算する。受信信号は、従来技術と同様にＭ成分を有するベクトルである。また、この計算は所定のサンプリングタイミングでなされる。
積分部１６は、相関行列計算部１４により所定のサンプリングタイミングで計算された相関行列を積分する。雑音の影響を小さくするためである。積分する時間間隔は任意で構わない。
【００３０】
固有値計算部１８は、相関行列の固有値と固有ベクトルを計算する。計算方法には、例えば、べき乗法、ヤコビ法、ＱＲ分解法などがあるが、ここでは任意のもので構わない。
関数処理部２０は、固有値、固有ベクトルに到来方向推定方法を使用して、信号発信源の位置を推定する。超音波探傷試験の場合、被検査材料中の欠陥が信号発信源とみなされる。到来方向推定方法には、例えば、ＭＵＳＩＣアルゴリズムがある。詳細は後述する。
【００３１】
強度判断部２２は、関数処理部２０によって推定された信号発信源より到来する信号の強度を計算する。信号の強度とは、例えば、信号の電力である。信号電力があらかじめ設定されているしきい値より小さい場合、その信号に対応する欠陥位置を欠陥位置の候補から除外する。
表示部２４は、以上の推定結果を出力し、これを所定の形式で表示する。これは、ディスプレイ等への表示だけでなく、ファイルへの書き込みも含むものとする。さらに、これらの結果をもとにした別の信号処理も含むものとする。また、外部へ出力しても構わない。
【００３２】
この構成は、ハードウエア的には、任意のコンピュータのＣＰＵ、メモリ、その他のＬＳＩで実現でき、ソフトウエア的にはメモリのロードされた予約管理機能のあるプログラムなどによって実現されるが、ここではそれらの連携によって実現される機能ブロックを描いている。したがって、これらの機能ブロックがハードウエアのみ、ソフトウエアのみ、またはそれらの組合せによっていろいろな形で実現できることは、当業者には理解されるところである。
【００３３】
図４は、位置推定処理についての手順を示す図である。ここでは、超音波探傷試験へ適用したものとする。受信センサ１０ａ〜ｍで受信され、信号変換部１２ａ〜ｍで所定の変換が施された受信信号に対して、相関行列計算部１４は相関行列を計算する（Ｓ１０）。積分部１６は相関行列を積分し（Ｓ１２）、相関行列における雑音の影響を小さくする。総サンプリング数をＳｐとした離散サンプリングにより受信信号Ｘを得た場合、積分処理後の相関関数Ｒは、
【００３４】
【数８】
2004144731.tif 000010
と表される。固有値計算部１８は、Ｒの固有値と固有ベクトルを計算する（Ｓ１４）。従来のＭＵＳＩＣアルゴリズムと同様に、式５より到来する信号数Ｋを決定する（Ｓ１６）。本実施の形態におけるＭＵＳＩＣアルゴリズム評価関数は以下の通りである。
【００３５】
【数９】
2004144731.tif 000011
ここで、ｘ、ｙ、ｚは所定の位置を原点とした直交座標系におけるパラメータであり、これらは有限空間に限定される。また、ステアリングベクトルａ（ｘ，ｙ，ｚ）は、
【００３６】
【数１０】
2004144731.tif 000012
である。例えば、（ｘ１，ｙ１，ｚ１）はひとつの受信センサ１０ａの位置を示す。このように関数処理部２０における評価関数は、ｘ、ｙ、ｚで表されるパラメータの座標とｘ１，ｙ１，ｚ１などで表される受信信号の入力位置における座標のちがいをもとに定まる関数で与えられるステアリングベクトルを含んでいる。ステアリングベクトルはパラメータｘ、ｙ、ｚと受信センサ１０ａ〜ｍの位置の距離の関数で表される。図５は、超音波探傷試験装置に使用した場合を２次元座標として表現した一例である。
【００３７】
ＭＵＳＩＣアルゴリズム評価関数におけるパラメータｘ、ｙ、ｚを被検査材料における座標の範囲内で掃引する（Ｓ１８）。図６が掃引結果の一例を示す図である。掃引結果より、Ｋ個の極大値を選択し、Ｋ個の信号発信源を特定する（Ｓ２０）。このＫ個の信号発信源に対応する座標が被検査材料中における欠陥位置である。図６では、Ｋを５として極大値を選択している。このように関数処理部２０は、評価関数をパラメータの座標で掃引し、その結果から特徴点を検出し、信号発信源である欠陥の位置を推定する。強度判断部２２は、各欠陥位置から到来する信号と推定された信号に対する受信電力を計算する（Ｓ２２）。受信電力は、次式に示す逆行列演算を行った後、
【００３８】
【数１１】
2004144731.tif 000013
この行列Ｓの第ｋ対角成分から第ｋの到来する信号成分に対する受信電力が得られる。この受信電力はしきい値と比較される。受信信号がしきい値より大きい場合、これに対応する信号発信源を選択し（Ｓ２４）、信号発信源の座標に欠陥が存在すると決定する。一方、小さい場合、それに対応する欠陥を除外する。この処理により、一般に到来方向推定方法における誤判定を軽減できる。このように、強度判断部２２は、推定された信号発生源からの受信信号の強度が、所定のしきい値を超えない場合に、受信信号に対する信号発生源を有意な検出対象から除外する処理を実行する。表示部２４は、以上の処理結果を所定の形式で出力または表示する（Ｓ２６）。
【００３９】
なお、以上の信号発信源の位置推定処理において、例えば受信センサ１０ａで兼ねられる送信部はひとつの周波数を有する信号を出力している。送信部自体は複数の周波数を発信可能でも構わないが、ここではその中のひとつの周波数のみを使用する。
【００４０】
このように本実施の形態によれば、被検査材料の有限空間性を利用することにより、各受信センサに到来する信号の入射角が等しい仮定を不要にできるため、信号源が受信センサの近傍に存在する場合においても、信号発信源の位置をより正確に推定することが可能になる。これは、ステアリングベクトルが直交座標系における座標をパラメータとする関数で表現され、そのパラメータを掃引する処理により実現される。これに加えて、周波数掃引なしで、ＭＵＳＩＣアルゴリズムにおいて信号発信源との距離を求めることも可能になる。さらに、受信電力が小さいと推定される信号に対する信号発信源を候補から除外するため、信号発信源の誤判定を減少させることが可能である。
【００４１】
また、超音波探傷試験に関して、往復伝播時間差ではなくＭＵＳＩＣアルゴリズムを使用して欠陥位置を推定するため、第一反射信号以外の入力信号による特性劣化が軽減され、高精度な超音波探傷試験を可能にする。
【００４２】
（実施の形態２）
実施の形態１では、位置が既知の受信センサにおいて受信された信号から、信号発信源の位置を推定したが、実施の形態２では、位置が既知の信号発信源からの信号を受信した位置について、受信した信号から推定する。実施の形態２においては、工場ライン等で働く稼動体ロボット上に、信号を受信可能な受信センサを設置し、位置が既知の信号発信源からの信号にもとづいて受信センサの位置を測定する。その結果、稼動体ロボットが位置を自己認識する。
【００４３】
図８は、稼動体ロボットシステム４２の構成を示す。稼動体ロボットシステム４２は、発信源３０と位置推定装置１００を含む。
発信源３０は、位置推定装置１００に対して所定の信号を出力する。なお、この信号には、位置推定装置１００の移動についての指示が含まれてもよい。
位置推定装置１００は、稼動体ロボットの位置を検出する。ここでは、位置推定装置１００が稼動体ロボット自体であるものとし、移動のために必要なモータ等を含んで、通路４０におけるｘやｙの方向に移動する。また、受信センサ１０ａ〜ｄを含む。
【００４４】
図９は、実施の形態２に係る位置推定装置１００の構成図である。位置推定装置１００は、受信センサ１０ａ〜ｍ、センサ移動部３２、回転部３４、ケーブル３６ａ〜ｍ、信号処理部３８を含む。なお、図９は、直交座標系において、ｙを一定の値とした場合のＸ−Ｚ平面を示し、さらに図示のとおり、所定の一点が原点とされる。
【００４５】
発信源３０は、所定の信号として音波や電磁波などを発信する。また、発信源３０の座標は既知であり、これを（ｘ１，ｚ１）とする。
受信センサ１０ａ〜ｍは、発信源３０から到来する信号を受信する。これらの絶対的な座標は不明であるが、受信センサ１０ａと受信センサ１０ｂ間の距離がｃ２であるような、それぞれの受信センサ１０ａ〜ｍ間の距離は既知である。また、ひとつの受信センサ１０ａの座標をパラメータ（ｘ，ｚ）とする。
センサ移動部３２は受信センサ１０ａ〜ｍを備え、回転部３４によって、ｘ方向に移動する。
【００４６】
ケーブル３６ａ〜ｍは、それぞれ受信センサ１０ａ〜ｍと後述する信号処理部３８を接続するケーブルである。
信号処理部３８は、受信した信号をもとに所定の処理を実行する。なお、信号処理部３８は、図３における位置推定装置１００の信号変換部１２ａ〜ｍ、相関行列計算部１４、積分部１６、固有値計算部１８、関数処理部２０、強度判断部２２、表示部２４の全部あるいは一部を含むものとする。
【００４７】
実施の形態２における位置推定手段についての手順としては、図４に示したものが有効である。ここで、ステップ１８におけるＭＵＳＩＣアルゴリズムの評価関数は式９と同一である。なお、ステアリングベクトルａ（ｘ，ｙ，ｚ）は、次のとおりに示される。
【数１２】
2004144731.tif 000014
なお、（ｘ１，ｙ１，ｚ１）は発信源３０の座標、（ｘ，ｙ，ｚ）はひとつの受信センサ１０ａの座標、（ａ２，ｂ２，ｃ２）はひとつの受信センサ１０ａを基準にした第２の受信センサ１０ｂの座標、（ａＭ，ｂＭ，ｃＭ）はひとつの受信センサ１０ａを基準にした第Ｍの受信センサ１０ｍの座標である。
【００４８】
以下、図４と同様の手順において、（ｘ，ｙ，ｚ）を推定する。さらに、（ｘ，ｙ，ｚ）からその他の受信センサ１０ｂ〜ｍの座標（ｘ＋ａ２，ｙ＋ｂ２，ｚ＋ｃ２）や（ｘ＋ａＭ，ｙ＋ｂＭ，ｚ＋ｃＭ）等を計算する。これらの座標に平均等の統計処理を行って、位置推定装置１００の位置を検出する。
【００４９】
以上の構成による位置推定装置１００の動作は、以下のとおりである。センサ移動部３２は、通路４０上の任意の場所に存在する。発信源３０は音波を発信する。受信センサ１０ａ〜ｍはそれぞれ音波を受信し、電気信号に変換された音波は、ケーブル３６ａ〜ｍによって信号処理部３８に伝送される。相関行列計算部１４は、これらの相関行列を計算し、積分部１６は所定期間において相関行列を積算する。固有値計算部１８は、積算された相関行列の固有値と固有ベクトルを計算し、関数処理部２０は、ＭＵＳＩＣアルゴリズムにもとづいてひとつの受信センサ１０ａの位置を推定する。さらに、ひとつの受信センサ１０ａの位置からその他の受信センサ１０ｂ〜ｍの位置を計算する。これらの座標に平均等の統計処理を行って、位置推定装置１００の位置を検出する。さらに、センサ移動部３２を通路４０上の別の場所に移動させ、上記の処理を繰り返し実行する。
【００５０】
本実施の形態によれば、直交座標系における座標がパラメータのステアリングベクトルを含んだ評価関数にもとづいて座標が推定されるため、パラメータの座標を受信センサの位置とした場合、受信センサの位置を高精度に推定可能である。
【００５１】
（実施の形態３）
実施の形態１と２では、座標がパラメータのステアリングベクトルを含んだ評価関数によって、ＭＵＳＩＣアルゴリズムから所定の位置を推定したが、実施の形態３では、座標がパラメータのステアリングベクトルを含んだ評価関数によって、ＭＵＳＩＣアルゴリズムから複数の受信センサにおける信号の振幅と位相の偏差を推定し、さらに推定結果にもとづいて、受信センサをキャリブレーションする。すなわち、実施の形態１と２における位置推定装置の使用の前段階において、実施の形態３によって受信センサのキャリブレーションを実行する。本実施の形態において、位置が既知である校正信号の信号発信源からの受信信号を固有値解析し、その結果の雑音部分空間が、信号発信源と受信センサ間の距離の関数を含んだステアリングベクトルに直交するように上述の偏差を推定する。
【００５２】
実施の形態３に係る位置推定装置１００の構成としては、図９に示したものが有効である。
前述の式１に対応して、偏差行列Ｇを考慮すれば、受信信号ベクトルＸは以下の通り表せる。
【数１３】
2004144731.tif 000015
偏差行列Ｇを補償行列Ｇ’と推定すれば、補正された受信信号ベクトルＸ’は、以下の通り表せる。
【数１４】
2004144731.tif 000016
【００５３】
実施の形態３においては、キャリブレーション処理のために補償行列Ｇ’を推定する。図１０は、補償行列の推定手順を示す。偏差行列Ｇに対して補償行列Ｇ’を推定する場合、繰返し回数ｎにおける補償行列Ｇｎ’に対して以下に示す繰返し計算を実行する。受信信号ベクトルＸの相関行列Ｒ０を計算した後、繰返し回数ｎを１にし（Ｓ３０）、さらにＧ１’を単位行列に設定する（Ｓ３２）。行列ＲｎをＧｎ’−１Ｒ０（Ｇｎ’Ｈ）−１から計算し、Ｒｎの固有値解析によって部分雑音空間Ｅｎを求める。校正信号の信号発信源の特定された位置を（ｘｋ，ｙｋ，ｚｋ）として、逐次最小化される以下の漸化式Ｚｎを計算する（Ｓ３４）。
【数１５】
2004144731.tif 000017
ここでは、校正信号の信号発信源をＫ箇所に配置しており、さらに、ステアリングベクトルａ（ｘｋ，ｙｋ，ｚｋ）は既知の値である。さらに、Ｇｎ＋１’を以下の通りに計算する（Ｓ３６）。
【数１６】
2004144731.tif 000018
【００５４】
ここで、ｗ＝［１，０，０，・・・，０］Ｔとする。さらに、以下の条件を満たさなければ（Ｓ３８のＮ）、ｎに１を加算し（Ｓ４０）、行列Ｒｎの計算を再び実行する。
【数１７】
2004144731.tif 000019
ここで‖・‖Ｆは、フロベニウスノルムである。一方、条件を満たせば（Ｓ３８のＹ）、Ｇｎ＋１’をＧ’に代入する（Ｓ４２）。このようにして、決定したＧ’によって受信信号ベクトルＸを補正する。
【００５５】
図３や図９の位置推定装置１００の受信センサ１０ａ〜ｍにおける振幅や位相のパターンとＭＵＳＩＣに用いるステアリングベクトルに誤差がある場合や、受信センサ１０ａ〜ｍ間に振幅・位相偏差がある場合、位置推定精度が低下する。この低下を可能な限り避けるため、一般にハードウエアキャリブレーション機能が備えられている。しかし、ハードウエアキャリブレーション機能を具備するためには回路構成が複雑になり、コストの増大を招く。そのため、ソフトウエアキャリブレーション法が確立できれば、ハードウエア面とともにコスト面からも非常に有利である。本実施の形態におけるキャリブレーションはソフトウエアによっても実施可能である。
【００５６】
本実施の形態によれば、複数の受信センサにおける信号の振幅と位相のキャリブレーション処理を固有値解析にもとづいて行うため、特性の改善を可能にする。さらに、ハードウエア処理だけでなく、ソフトウエア処理も可能なので、装置規模の削減も可能にする。
【００５７】
以上、本発明を実施の形態をもとに説明した。この実施の形態は例示であり、それら各構成要素や各処理プロセスの組み合わせにいろいろな変形例が可能なこと、またそうした変形例も本発明の範囲であることは当業者に理解されるところである。そうした例を挙げる。
【００５８】
実施の形態１では、本発明が超音波探傷試験の欠陥位置推定に適用される場合を説明した。しかし、本発明はこれ以外に適用されてもよい。例えば、電波伝搬環境下における無線通信などにも適用してもよい。不要波を出力する送信装置が受信装置の近傍に存在する場合における、送信装置の座標を推定する場合である。要は、ステアリングベクトルのパラメータに対する有限空間性が成立する場合すればよい。
【００５９】
実施の形態１では、到来方向推定方法としてＭＵＳＩＣアルゴリズムを使用している。しかし、ＭＵＳＩＣアルゴリズム以外を使用してもかまわない。例えば、最小ノルム法やＥＳＰＲＩＴ（Ｅｓｔｉｍａｔｉｏｎ　ｏｆ　Ｓｉｇｎａｌ　Ｐａｒａｍｅｔｅｒｓ　ｖｉａ　Ｒｏｔａｔｉｏｎ　Ｉｎｖａｒｉａｎｃｅ　Ｔｅｃｈｎｉｑｕｅｓ）アルゴリズムなどである。要は、固有値解析にもとづいて、信号発信源が推定されればよい。
【図面の簡単な説明】
【図１】従来の技術に係る超音波探傷試験装置の構成を示す図である。
【図２】従来の技術に係るアレーアンテナの構成を示す図である。
【図３】実施の形態１に係る位置推定装置の構成を示す図である。
【図４】図３の装置による位置推定処理についての手順を示す図である。
【図５】図３の装置を超音波探傷試験装置に使用した場合の構成を示す図である。
【図６】図３の装置による評価関数の掃引結果を示す図である。
【図７】図３の装置によって推定された各座標の到来波発信源からの受信電力を示す図である。
【図８】実施の形態２に係る稼動体ロボットシステムの構成を示す図である。
【図９】図８の位置推定装置の構成を示す図である。
【図１０】実施の形態３に係る補償行列の推定手順を示す図である。
【符号の説明】
１０ａ〜ｍ　受信センサ、　１２ａ〜ｍ　信号変換部、　１４　相関行列計算部、　１６　積分部、　１８　固有値計算部、　２０　関数処理部、　２２　強度判断部、　２４　表示部、　３０　発信源、　３２　センサ移動部、　３４　回転部、　３６ａ〜ｍ　ケーブル、　３８　信号処理部、　４０　通路、　４２稼動体ロボットシステム、　５０　送受信センサ、　５２　受信センサ、　６０ａ〜ｍ　アンテナ素子、　６２　無線装置、　１００　位置推定装置、　１２０　アレーアンテナ装置。
【図１】
2004144731.tif 000020
【図２】
2004144731.tif 000021
【図３】
2004144731.tif 000022
【図４】
2004144731.tif 000023
【図５】
2004144731.tif 000024
【図６】
2004144731.tif 000025
【図７】
2004144731.tif 000026
【図８】
2004144731.tif 000027
【図９】
2004144731.tif 000028
【図１０】
2004144731.tif 000029