▶ 株式会社デンソーＥＭＣエンジニアリングサービスの特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2022134585
(43)【公開日】2022-09-15
(54)【発明の名称】信号源推定装置
(51)【国際特許分類】
G01R 29/08 20060101AFI20220908BHJP
G01R 29/10 20060101ALI20220908BHJP
H04B 17/309 20150101ALI20220908BHJP
【ＦＩ】
G01R29/08 D
G01R29/10 Z
H04B17/309
【審査請求】有
【請求項の数】5
【出願形態】ＯＬ
(21)【出願番号】P 2021033796
(22)【出願日】2021-03-03
(71)【出願人】
【識別番号】502135679
【氏名又は名称】株式会社デンソーＥＭＣエンジニアリングサービス
(74)【代理人】
【識別番号】100106149
【弁理士】
【氏名又は名称】矢作 和行
(74)【代理人】
【識別番号】100121991
【弁理士】
【氏名又は名称】野々部 泰平
(74)【代理人】
【識別番号】100145595
【弁理士】
【氏名又は名称】久保 貴則
(72)【発明者】
【氏名】赤堀 一郎
(72)【発明者】
【氏名】下起 洋司
(72)【発明者】
【氏名】杉浦 大地
(57)【要約】
【課題】ターゲット信号に影響する信号の信号源を容易に推定することができる信号源推定装置を提供する。
【解決手段】解析対象の信号であるターゲット信号ｙと、プローブにより検出された信号であるプローブ信号とのコヒーレンス関数γＺＹ ２を算出するコヒーレンス関数算出部８０と、プローブ信号ｚのパワースペクトルＷＺＺを算出するパワースペクトル算出部５２と、コヒーレンス関数γＺＹ ２と、プローブ信号スペクトルＺのパワースペクトルＷＺＺを平均化した値との積ＷＶＺを算出し、算出した値を出力する乗算部９２とを備える信号源推定装置。
【選択図】図２
【解決手段】解析対象の信号であるターゲット信号ｙと、プローブにより検出された信号であるプローブ信号とのコヒーレンス関数γＺＹ ２を算出するコヒーレンス関数算出部８０と、プローブ信号ｚのパワースペクトルＷＺＺを算出するパワースペクトル算出部５２と、コヒーレンス関数γＺＹ ２と、プローブ信号スペクトルＺのパワースペクトルＷＺＺを平均化した値との積ＷＶＺを算出し、算出した値を出力する乗算部９２とを備える信号源推定装置。
【選択図】図２
【特許請求の範囲】
【請求項1】
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ
２）を算出するコヒーレンス関数算出部（８０、１９０）と、
前記検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
前記コヒーレンス関数と、平均化した前記検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、
を備える信号源推定装置。
前記検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
前記コヒーレンス関数と、平均化した前記検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、
を備える信号源推定装置。
【請求項2】
請求項１に記載の信号源推定装置であって、
前記検出子の位置を検出する位置検出部（１１０）と、
前記乗算部が算出した前記積を、前記位置検出部が検出した前記検出子の位置に対応付けて、表示装置（３０）に表示する表示処理部（２２）と、
を備える、信号源推定装置。
前記検出子の位置を検出する位置検出部（１１０）と、
前記乗算部が算出した前記積を、前記位置検出部が検出した前記検出子の位置に対応付けて、表示装置（３０）に表示する表示処理部（２２）と、
を備える、信号源推定装置。
【請求項3】
請求項２に記載の信号源推定装置であって、
前記検出子により信号源の位置を推定する範囲を撮影するカメラ（１１）を備え、
前記表示処理部は、前記表示装置に、前記カメラが撮影した画像を視認可能にしつつ、前記カメラが撮影した画像に、前記検出子の二次元位置と前記積との関係を示す画像を重畳表示する、信号源推定装置。
前記検出子により信号源の位置を推定する範囲を撮影するカメラ（１１）を備え、
前記表示処理部は、前記表示装置に、前記カメラが撮影した画像を視認可能にしつつ、前記カメラが撮影した画像に、前記検出子の二次元位置と前記積との関係を示す画像を重畳表示する、信号源推定装置。
【請求項4】
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ
２）を算出するコヒーレンス関数算出部（８０、１９０）と、
前記検出子の位置を検出する位置検出部（１１０）と、
前記コヒーレンス関数を、前記位置検出部が検出した前記検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、
を備える信号源推定装置。
前記検出子の位置を検出する位置検出部（１１０）と、
前記コヒーレンス関数を、前記位置検出部が検出した前記検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、
を備える信号源推定装置。
【請求項5】
請求項４に記載の信号源推定装置であって、
前記検出子により信号源の位置を推定する範囲を撮影するカメラ（１１）を備え、
前記表示処理部は、前記表示装置に、前記カメラが撮影した画像を視認可能にしつつ、前記カメラが撮影した画像に、前記検出子の二次元位置と前記コヒーレンス関数との関係を示す画像を重畳表示する、信号源推定装置。
前記検出子により信号源の位置を推定する範囲を撮影するカメラ（１１）を備え、
前記表示処理部は、前記表示装置に、前記カメラが撮影した画像を視認可能にしつつ、前記カメラが撮影した画像に、前記検出子の二次元位置と前記コヒーレンス関数との関係を示す画像を重畳表示する、信号源推定装置。
【発明の詳細な説明】
【技術分野】
【0001】
信号源の位置を推定する信号源推定装置に関する。
【背景技術】
【0002】
特許文献１には、被測定対象物からどのような距離でどのような信号が発生しているかを表示する装置が開示されている。特許文献１に開示されている装置は、被測定対象物から発生する電磁界をプローブ等により検出し、プローブ等の位置に対応する画面の区画を、強度に応じた色で表示する。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特開２０１３－２３８５８１号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
プローブ等により検出した電磁界の強度には、解析対象の信号（以下、ターゲット信号）に影響しない信号源の信号も影響する。したがって、特許文献１に開示された装置は、ターゲット信号に影響する信号源を推定することが必ずしも容易ではない。
【0005】
本開示は、この事情に基づいて成されたものであり、本開示の目的は、ターゲット信号に影響する信号の信号源を容易に推定することができる信号源推定装置を提供することにある。
【課題を解決するための手段】
【0006】
上記目的は独立請求項に記載の特徴の組み合わせにより達成され、また、下位請求項は更なる有利な具体例を規定する。特許請求の範囲に記載した括弧内の符号は、一つの態様として後述する実施形態に記載の具体的手段との対応関係を示すものであって、開示した技術的範囲を限定するものではない。
【0007】
上記目的を達成するための１つの開示は、
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
コヒーレンス関数と、平均化した検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、を備える信号源推定装置である。
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
コヒーレンス関数と、平均化した検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、を備える信号源推定装置である。
【0008】
ターゲット信号と検出子信号とのコヒーレンス関数は、詳しくは後述するが、ターゲット信号のパワースペクトルに占める、検出子が検出した検出子信号に対応する信号源信号のパワースペクトルの割合を示す。コヒーレンス関数と検出子信号のパワースペクトルの積が大きな値になる場合、ターゲット信号のパワースペクトルに占める信号源信号のパワースペクトルの割合が高く、かつ、検出子信号のパワースペクトルも大きいことになる。よって、この積が大きな値になる検出子信号を検出する位置が、ターゲット信号に影響を与えている信号源の位置を示していると言える。したがって、この積を用いることで、ターゲット信号に影響を与えている信号源の位置を容易に推定できる。
【0009】
上記目的を達成するための他の１つの開示は、
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子の位置を検出する位置検出部（１１０）と、
コヒーレンス関数を、位置検出部が検出した検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、を備える信号源推定装置である。
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子の位置を検出する位置検出部（１１０）と、
コヒーレンス関数を、位置検出部が検出した検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、を備える信号源推定装置である。
【0010】
ターゲット信号と検出子信号とのコヒーレンス関数は、ターゲット信号のパワースペクトルに占める、検出子が検出した検出子信号に対応する信号源信号のパワースペクトルの割合を示す。信号源推定装置は、このコヒーレンス関数を、検出子の位置に対応付けて表示装置に表示する。よって、操作者は、表示装置を見ることで、ターゲット信号に影響を与えている信号源の位置を容易に推定できる。
【図面の簡単な説明】
【0011】
【図1】信号源推定システム１０の全体構成図である。
【図2】処理装置２０の信号処理部２１が実行する機能を示すブロック図である。
【図3】一例とする測定系をモデル化した図である。
【図4】測定系を一般化した図である。
【図5】処理装置２０の表示処理部２２が実行する機能を示すブロック図である。
【図6】オーバーレイ画像を例示する図である。
【図7】オーバーレイ後画像を例示する図である。
【図8】積ＷＶＺの特性を説明するためのシミュレーションを示す図である。
【図10】第２実施形態の信号処理部２２１の構成を示すブロック図である。
【発明を実施するための形態】
【0012】
＜第１実施形態＞
以下、実施形態を図面に基づいて説明する。図１は、信号源推定システム１０の全体構成図である。信号源推定システム１０は、カメラ１１、ターゲット信号検出線１２、プローブ１４、処理装置２０、表示装置３０を備えている。
以下、実施形態を図面に基づいて説明する。図１は、信号源推定システム１０の全体構成図である。信号源推定システム１０は、カメラ１１、ターゲット信号検出線１２、プローブ１４、処理装置２０、表示装置３０を備えている。
【0013】
カメラ１１は、プローブ１４により信号源の位置を推定する範囲を撮影するために、所定の位置に固定されている。信号源は、ターゲット信号ｙに含まれる信号源信号ｘを発生させている部分である。ターゲット信号ｙは解析対象の信号である。解析する内容は、ターゲット信号ｙに含まれる信号源信号ｘの発生位置を推定することである。信号源信号ｘがターゲット信号ｙにおいて不要な信号であれば、信号源信号ｘはノイズであることになり、信号源はノイズ源である。
【0014】
ターゲット信号検出線１２は、ターゲット信号ｙを検出しようとする箇所に接続されている。ターゲット信号ｙは、ターゲット信号検出線１２を用いて検出する以外に、電界プローブあるいは磁界プローブを用いて検出したり、アンテナにより検出したりすることができる。
【0015】
プローブ１４は、検出子信号であるプローブ信号ｚを検出するための検出子である。プローブ信号ｚは、空間を伝播する電磁信号である。プローブ信号ｚは、信号源を推定するために検出する信号であるが、プローブ１４の伝達関数Ｈによる影響だけ、信号源信号ｘとは相違する。
【0016】
プローブ１４は、操作者が手に持って移動させる。操作者は、プローブ１４を、ターゲット信号ｙに影響する信号源を推定する対象装置に接近させて、その信号源の可能性がある範囲を順次移動させる。図１では、ＥＣＵ１６が備える複数の電子機器１８がそれぞれ信号源である。ただし、全部の電子機器１８が、ターゲット信号ｙに影響を与えている信号源であるとは限らない。そこで、この信号源推定システム１０により、ターゲット信号ｙに影響を与えている信号源がどれであるかを推定する。なお、プローブ信号ｚおよびターゲット信号ｙを検出する対象装置は、ＥＣＵ１６のような小型のものに限られず、車両などであってもよい。
【0017】
処理装置２０は、信号源推定装置に相当するものであり、少なくとも１つのプロセッサを備えた構成により実現できる。たとえば、処理装置２０は、ＣＰＵ、ＲＯＭ、ＲＡＭ、Ｉ／Ｏ、およびこれらの構成を接続するバスラインなどを備えたコンピュータにより実現できる。ＲＯＭには、汎用的なコンピュータを処理装置２０として機能させるためのプログラムが格納されている。ＣＰＵが、ＲＡＭの一時記憶機能を利用しつつ、ＲＯＭに記憶されたプログラムを実行することで、処理装置２０は、信号処理部２１および表示処理部２２として機能する。処理装置２０が、信号処理部２１および表示処理部２２として機能することは、上記プログラムに対応する方法が実行されることを意味する。
【0018】
処理装置２０が実行する処理の詳細は後述するが、処理装置２０は、カメラ１１が撮影した対象装置の画像（以下、カメラ画像）に、オーバーレイ画像を重畳させたオーバーレイ後画像の画像データ（以下、オーバーレイ後画像データ）を作成する。オーバーレイ画像の一例は、コヒーレンス関数γＺＹ
２とプローブ１４の二次元位置との関係を示す画像である。オーバーレイ画像の他の例は、コヒーレンス関数γＺＹ
２とプローブ信号スペクトルＺのパワースペクトルＷＺＺとの積ＷＶＺとプローブ１４の二次元位置との関係を示す画像である。処理装置２０は、オーバーレイ後画像を表示装置３０に表示させる。
【0019】
表示装置３０は、二次元表示およびカラー表示が可能な表示装置である。表示装置３０は、オーバーレイ後画像などを操作者に向けて表示する。本実施形態の表示装置３０は、表示画面に、可視化表示領域３１とグラフ表示領域３２を表示できる。
【0020】
［信号処理部２１の説明］
信号処理部２１は、図２に詳細を示す各部を備える。図２をもとに、信号処理部２１が実行する機能を説明する。信号処理部２１は、周波数変換部４０、パワースペクトル算出部５０、クロススペクトル算出部６０、平均化処理部７０、コヒーレンス関数算出部８０、乗算部９０、乗算部９２を備える。
信号処理部２１は、図２に詳細を示す各部を備える。図２をもとに、信号処理部２１が実行する機能を説明する。信号処理部２１は、周波数変換部４０、パワースペクトル算出部５０、クロススペクトル算出部６０、平均化処理部７０、コヒーレンス関数算出部８０、乗算部９０、乗算部９２を備える。
【0021】
周波数変換部４０は、２つのフレーム切り出し部４１、４２と、２つの離散フーリエ変換部（以下、ＤＦＴ部）４３、４４とを備える。一方のフレーム切り出し部４１には、図示しないＡＤ変換回路によりデジタル値に変換されたターゲット信号ｙが入力される。フレーム切り出し部４１は、逐次入力されるターゲット信号ｙから、ＤＦＴ部４３での一度の処理に用いる時間（すなわち窓）分の信号を切り出す。そして、切り出した信号をＤＦＴ部４３に入力する。
【0022】
他方のフレーム切り出し部４２には、図示しないＡＤ変換回路によりデジタル値に変換されたプローブ信号ｚが入力される。フレーム切り出し部４２は、逐次入力されるプローブ信号ｚから、ＤＦＴ部４４での一度の処理に用いる時間分の信号を切り出す。そして、切り出した信号をＤＦＴ部４４に入力する。
【0023】
ＤＦＴ部４３は、入力されたターゲット信号ｙを離散フーリエ変換する。フーリエ変換後の信号は、ターゲット信号ｙが周波数スペクトルに変換された信号である。変換後の信号をターゲット信号スペクトルＹとする。
【0024】
ＤＦＴ部４４は、入力されたプローブ信号ｚを離散フーリエ変換する。フーリエ変換後の信号は、プローブ信号ｚが周波数スペクトルに変換された信号である。変換後の信号をプローブ信号スペクトルＺとする。
【0025】
パワースペクトル算出部５０は、ターゲット信号スペクトルＹが入力される第１パワースペクトル算出部５１とプローブ信号スペクトルＺが入力される第２パワースペクトル算出部５２とを備える。第１パワースペクトル算出部５１は、ターゲット信号スペクトルＹのパワースペクトル、すなわち、｜Ｙ｜２を算出する。第２パワースペクトル算出部５２は、プローブ信号スペクトルＺのパワースペクトル、すなわち、｜Ｚ｜２を算出する。
【0026】
クロススペクトル算出部６０は、プローブ信号スペクトルＺとターゲット信号スペクトルＹのクロススペクトルＷＺＹを算出する。本実施形態では、具体的には、プローブ信号スペクトルＺの共役複素数Ｚ＊とターゲット信号スペクトルＹとの積を算出する。
【0027】
平均化処理部７０は、第１平均化処理部７１、第２平均化処理部７２、第３平均化処理部７３を備える。第１平均化処理部７１は、ターゲット信号スペクトルＹのパワースペクトルを平均化する。ターゲット信号スペクトルＹのパワースペクトルの平均値をＷＹＹとする。第２平均化処理部７２は、プローブ信号スペクトルＺのパワースペクトルを平均化する。プローブ信号スペクトルＺのパワースペクトルの平均値をＷＺＺとする。第３平均化処理部７３はクロススペクトルを平均化する。クロススペクトルの平均値をＷＺＹとする。なお、以降、記号Ｗで示す値は、平均であることを表記していなくても平均化処理された後の値を意味するものとし、「平均」は省略する。
【0028】
これら平均化処理部７０においてスペクトルを平均化する時間は任意に設定可能である。スペクトルを平均化する時間は、たとえば、数十ミリ秒程度とすることができる。また、スペクトルを平均化する時間は、数秒程度としてもよい。３つの平均化処理部７１、７２、７３の平均化時間は互いに同一とすることができる。平均化処理の手法は、ｎ項移動平均、指数移動平均など、種々の平均化手法を採用することができる。
【0029】
コヒーレンス関数算出部８０は、（１）式に示すコヒーレンス関数γＺＹ
２を算出する。コヒーレンス関数算出部８０は、算出したコヒーレンス関数γＺＹ
２を表示処理部２２に送る。コヒーレンス関数γＺＹ
２は、後述するように、ターゲット信号ｙに占める信号源信号ｘに起因するパワースペクトルの割合を示す。
【0030】
【数1】
【0031】
乗算部９２は、プローブ信号スペクトルＺのパワースペクトルＷＺＺとコヒーレンス関数γＺＹ
２との積ＷＶＺを算出し、その積ＷＶＺを表示処理部２２に送る。なお、平均化処理部７１が算出したターゲット信号スペクトルＹのパワースペクトルＷＹＹおよび平均化処理部７２が算出したプローブ信号スペクトルＺのパワースペクトルＷＺＺも、表示処理部２２に送られる。
【0032】
［数学的説明］
次に、コヒーレンス関数γＺＹ ２がターゲット信号ｙに占める信号源信号ｘに起因するパワースペクトルの割合を示すことを理論的に説明する。
次に、コヒーレンス関数γＺＹ ２がターゲット信号ｙに占める信号源信号ｘに起因するパワースペクトルの割合を示すことを理論的に説明する。
【0033】
図３に、一例とする測定系をモデル化した図を示す。図３に示す測定系では、３つの信号源Ａ、Ｂ、Ｃからの信号ｘＡ、ｘＢ、ｘＣが、それぞれ、伝達ゲインＧＡ、ＧＢ、ＧＣで伝達された後、合成されてターゲット信号ｙとなっている。
【0034】
各信号ｘＡ、ｘＢ、ｘＣは相互に無相関であるとする。このとき、ある１つの信号源に着目したとき、他の信号源はノイズｎとして１つにまとめることができる。したがって、図３に例示する測定系を含め、どのような複雑な測定系でも、測定系は図４に一般化することができる。
【0035】
図４より、ターゲット信号ｙの周波数スペクトルであるターゲット信号スペクトルＹは（２）式で表すことができる。また、プローブ信号ｚの周波数スペクトルであるプローブ信号スペクトルＺは（３）式で表すことができる。なお、Ｇは信号源信号ｘがターゲット信号ｙの観測点まで伝達されるときの伝達関数、Ｈはプローブ１４の伝達関数である。Ｘは信号源信号ｘの周波数スペクトル、Ｎはノイズの周波数スペクトルである。
【0036】
Ｙ＝ＧＸ＋Ｎ （２）
Ｚ＝ＨＸ （３）
（２）式、（３）式より、ＷＺＺは（４）式の変形ができ、ＷＺＹは（５）式の変形ができる。
Ｚ＝ＨＸ （３）
（２）式、（３）式より、ＷＺＺは（４）式の変形ができ、ＷＺＹは（５）式の変形ができる。
【数4】
【0037】
【数5】
【数6】
【0038】
コヒーレンス関数γＺＹ
２は、プローブ信号ｚを用いて算出しており、プローブ信号ｚと信号源信号ｘの間には、伝達関数Ｈの影響分の相違がある。しかし、上記式を用いた証明により、コヒーレンス関数γＺＹ
２には、伝達関数Ｈは影響しないことが示された。
【0039】
さらに、Ｖ＝ＧＸであるので、（７）式が成立する。なお、Ｖは、信号ｖの周波数スペクトルであり、信号ｖは信号源信号ｘが伝達関数Ｇによりターゲット信号ｙの観測点まで伝達された信号である。
【0040】
【数7】
（７）式を変形すると（８）式が得られる。
【0041】
【数8】
【0042】
［表示処理部２２の説明］
表示処理部２２は、図５に詳細を示す各部を備える。次に、図５をもとに、表示処理部２２が実行する機能を説明する。表示処理部２２は、データ読み込み部１００、位置検出部１１０、メモリ１２０、グラフデータ作成部１３０、オーバーレイデータ作成部１４０、オーバーレイ処理部１５０を備えている。
表示処理部２２は、図５に詳細を示す各部を備える。次に、図５をもとに、表示処理部２２が実行する機能を説明する。表示処理部２２は、データ読み込み部１００、位置検出部１１０、メモリ１２０、グラフデータ作成部１３０、オーバーレイデータ作成部１４０、オーバーレイ処理部１５０を備えている。
【0043】
データ読み込み部１００は、ＷＶＶ、ＷＹＹ、ＷＺＺ、ＷＶＺ、γＺＹ
２を信号処理部２１から読み込む。以下、データ読み込み部１００が信号処理部２１から読み込んだデータを総称して読み込みデータとする。
【0044】
位置検出部１１０は、カメラ画像を取得し、そのカメラ画像をもとにプローブ１４の二次元位置を逐次検出する。プローブ１４によりプローブ信号ｚを検出している状態では、カメラ画像にプローブ１４が写っている。したがって、カメラ画像をもとにプローブ１４の二次元位置を検出できる。
【0045】
メモリ１２０は、データ読み込み部１００が読み込んだ読み込みデータと、位置検出部１１０が検出したプローブ１４の位置とを対応付けて記憶する。このメモリ１２０には、処理装置２０が備えるＲＡＭを用いることができる。メモリ１２０に記憶したデータは、ハードディスクなどの不揮発性の記憶装置に保存してもよい。
【0046】
グラフデータ作成部１３０は、表示装置３０のグラフ表示領域３２にグラフを描画するためのグラフ描画データを作成する。グラフ表示領域３２に描画するグラフは、たとえば、横軸を周波数とする。縦軸は、選択した区画のターゲット信号スペクトルＹのパワースペクトルＷＹＹ、プローブ信号スペクトルＺのパワースペクトルＷＺＺ、コヒーレンス関数γＺＹ
２、プローブ信号スペクトルＺのパワースペクトルＷＺＺとコヒーレンス関数γＺＹ
２との積ＷＶＺとする。また、縦軸をパワースペクトルＷＶＶとしたグラフをグラフ表示領域３２に表示してもよい。
【0047】
グラフの横軸に含まれる周波数帯は、操作者の入力操作によりグラフデータ作成部１３０に入力される周波数範囲により変更可能である。また、区画は、可視化表示領域３１の１つの区画であり、操作者の入力操作により定まる。
【0048】
オーバーレイデータ作成部１４０は、カメラ画像に重畳するオーバーレイデータを作成する。オーバーレイデータは、オーバーレイ画像を示すデータである。オーバーレイ画像は、プローブ１４の二次元位置と読み込みデータが示す値の大きさとの対応を示す二次元画像である。プローブ１４の二次元位置は、可視化表示領域３１の各区画に対応した離散化した位置である。読み込みデータが示す値の大きさは、色もしくは濃淡で表す。オーバーレイデータ作成部１４０は、可視化表示領域３１の各区画の読み込みデータをメモリ１２０から取得し、取得したデータをもとに、オーバーレイデータを作成する。
【0049】
オーバーレイデータ作成部１４０にも、操作者の入力操作により周波数範囲が入力される。オーバーレイデータ作成部１４０に入力される周波数範囲は、グラフデータ作成部１３０に入力される周波数範囲とは異なる範囲を設定できる。読み込みデータは周波数により数値が異なる。各区画に対応する読み込みデータの値は、周波数範囲の最大値、平均値など、解析目的に応じて定まる値とする。
【0050】
図６に、オーバーレイ画像を例示する。図６において、左上に示しているカメラ画像以外の３つの画像が、オーバーレイ画像である。図６の右上の画像は、可視化表示領域３１の区画ごとにプローブ信号スペクトルＺのパワースペクトルＷＺＺを示す画像である。図６の左下の画像は、可視化表示領域３１の区画ごとにコヒーレンス関数γＺＹ
２を示す画像である。図６の右下の画像は、可視化表示領域３１の区画ごとに積ＷＶＺを示す画像である。
【0051】
オーバーレイ処理部１５０は、表示装置３０の可視化表示領域３１に表示させるオーバーレイ後画像を示すデータであるオーバーレイ後画像データを作成する。そして、作成したオーバーレイ後画像データを表示装置３０に出力する。これにより、表示装置３０の可視化表示領域３１にオーバーレイ後画像が表示される。
【0052】
図７に、オーバーレイ後画像を例示する。図７において、左上に示しているカメラ画像以外の３つの画像がオーバーレイ後画像である。図７の右上の画像は、カメラ画像に図６の右上の画像を重畳した画像である。図７の左下の画像は、カメラ画像に図６の左下の画像を重畳した画像である。図７の右下の画像は、カメラ画像に図６の右下の画像を重畳した画像である。いずれのオーバーレイ後画像も、カメラ画像を視認可能にしつつも、カメラ画像に、オーバーレイ画像が重畳表示された画像である。
【0053】
オーバーレイ後画像が表示されることで、ノイズ源（すなわち信号源）とノイズの大きさとの関係を容易に理解できる。図７には、４つの画像が並んで示されているが、可視化表示領域３１には、選択的に１つの画像を表示する。ただし、これとは別の態様として、事前に設定された１つあるいは複数の画像が可視化表示領域３１に表示されるようになっていてもよい。
【0054】
［３つのオーバーレイ後画像の違い］
図７に示す３つのオーバーレイ後画像の違いを説明する。プローブ信号ｚのパワースペクトルＷＺＺを示す画像には、３つのノイズ源１、２、３が観測できる。しかし、これら３つのノイズ源１、２、３がターゲット信号ｙに影響を与えているかどうかは、この画像からは分からない。
図７に示す３つのオーバーレイ後画像の違いを説明する。プローブ信号ｚのパワースペクトルＷＺＺを示す画像には、３つのノイズ源１、２、３が観測できる。しかし、これら３つのノイズ源１、２、３がターゲット信号ｙに影響を与えているかどうかは、この画像からは分からない。
【0055】
コヒーレンス関数γＺＹ
２を示す画像を見ると、左下の１つの領域のみが、コヒーレンス関数γＺＹ
２が高いことが分かる。この画像を見ることで、ターゲット信号ｙに影響を与えているノイズ源は、ノイズ源２であることが分かる。さらに、積ＷＶＺを示す画像を見ると、ターゲット信号ｙに影響を与えているノイズ源が、ノイズ源２であることが、より明確に分かる。
【0056】
コヒーレンス関数γＺＹ
２と積ＷＶＺとの違いを示すために、図８に示すシミュレーションを説明する。図８に示すシミュレーションでは、２個のノイズ源１、２をｘ軸上の＋３０ｍｍおよび－３０ｍｍの位置に配置する。２個のノイズ源１、２の周波数差は周波数分解能より小さい。ターゲット信号ｙに対してはノイズ源１のみが寄与している。高さを固定してプローブ１４を走査させて、プローブ信号ｚを検出する。このプローブ信号ｚをもとに、コヒーレンス関数γＺＹ
２と、積ＷＶＺを算出する。
【0057】
図９に、プローブ信号スペクトルＺのパワースペクトルＷＺＺ、コヒーレンス関数γＺＹ
２、積ＷＶＺの位置変化を示す。図９上段のグラフに示されるように、プローブ信号スペクトルＺのパワースペクトルＷＺＺは、－３０ｍｍと３０ｍｍの位置にピークがある。図９中段のグラフに示されるように、コヒーレンス関数γＺＹ
２は、ノイズ源２の位置では０になっている。これに対して、ノイズ源１の位置では１になっている。ただし、ノイズ源１の位置以外でも、コヒーレンス関数γＺＹ
２は１など高い値になっている。しかし、ノイズ源の位置は分かっているのが通常である。図８の例では、ノイズ源の位置は－３０ｍｍおよび３０ｍｍであることは分かっている。したがって、３０ｍｍの位置のコヒーレンス関数γＺＹ
２が０であることが分かれば、ターゲット信号ｙに影響を与えているノイズ源はノイズ源１であることが分かる。このように、位置とコヒーレンス関数γＺＹ
２の関係から、ターゲット信号ｙに影響を与えているノイズ源を推定することができる。
【0058】
【0059】
［第１実施形態のまとめ］
以上、説明した第１実施形態の信号源推定システム１０によれば、ターゲット信号ｙとプローブ信号ｚを取得し、それらの信号から、コヒーレンス関数γＺＹ ２を算出する。プローブ信号ｚに伝達関数Ｈの影響があるにも関わらず、コヒーレンス関数γＺＹ ２は、上述した数学的説明が示すように、ターゲット信号ｙのパワースペクトルに占める信号源信号ｘのパワースペクトルの割合を表している。
以上、説明した第１実施形態の信号源推定システム１０によれば、ターゲット信号ｙとプローブ信号ｚを取得し、それらの信号から、コヒーレンス関数γＺＹ ２を算出する。プローブ信号ｚに伝達関数Ｈの影響があるにも関わらず、コヒーレンス関数γＺＹ ２は、上述した数学的説明が示すように、ターゲット信号ｙのパワースペクトルに占める信号源信号ｘのパワースペクトルの割合を表している。
【0060】
信号源推定システム１０は、このコヒーレンス関数γＺＹ
２を、プローブ１４の位置に対応付けて表示装置３０に表示する。よって、信号源推定システム１０を操作する操作者は、表示装置３０を見ることで、ターゲット信号ｙに影響を与えている信号源の位置を容易に推定できる。
【0061】
特に、本実施形態では、コヒーレンス関数γＺＹ
２をプローブ１４の位置に対応付けた画像として、オーバーレイ後画像を表示装置３０に表示する。コヒーレンス関数γＺＹ
２のオーバーレイ後画像は、カメラ画像を視認可能にしつつ、カメラ画像に、コヒーレンス関数γＺＹ
２とプローブ１４の二次元位置とを対応付けた画像が重畳された画像である（図７左下画像）。したがって、どの信号源が、ターゲット信号ｙに影響を与えている信号源であるかを特に推定しやすい。
【0062】
また、本実施形態では、プローブ信号スペクトルＺのパワースペクトルＷＺＺと、コヒーレンス関数γＺＹ
２の積ＷＶＺを算出する。この積ＷＶＺが大きな値になる場合、ターゲット信号ｙのパワースペクトルＷＹＹに占める信号源信号ｘのパワースペクトルＷＸＸの割合が高く、かつ、プローブ信号ｚのパワースペクトルＷＺＺも大きいことになる。よって、この積ＷＶＺが大きな値になるプローブ信号ｚを検出する位置が、ターゲット信号ｙに影響を与えている信号源の位置であると推定できる。
【0063】
プローブ１４の二次元位置と積ＷＶＺとの関係を示す画像も、カメラ画像を視認可能にしつつ、カメラ画像に重畳表示されている（図７右下画像）。操作者は、この画像を見ることで、ターゲット信号ｙに影響を与えている信号源の位置を明確に推定できる。
【0064】
＜第２実施形態＞
次に、第２実施形態を説明する。この第２実施形態以下の説明において、それまでに使用した符号と同一番号の符号を有する要素は、特に言及する場合を除き、それ以前の実施形態における同一符号の要素と同一である。また、構成の一部のみを説明している場合、構成の他の部分については先に説明した実施形態を適用できる。
次に、第２実施形態を説明する。この第２実施形態以下の説明において、それまでに使用した符号と同一番号の符号を有する要素は、特に言及する場合を除き、それ以前の実施形態における同一符号の要素と同一である。また、構成の一部のみを説明している場合、構成の他の部分については先に説明した実施形態を適用できる。
【0065】
図１０に第２実施形態の信号処理部２２１を示す。信号処理部２２１は、第１実施形態の信号処理部２１に代えて用いることができるものである。
【0066】
信号処理部２２１は、信号処理部２１が備えていたコヒーレンス関数算出部８０、乗算部９０を備えていない。代わりに、信号処理部２２１は、伝達関数推定部１８０、２乗部１８２、乗算部１８４、コヒーレンス関数算出部１９０を備える。信号処理部２２１の他の構成は、信号処理部２１と同じである。
【0067】
伝達関数推定部１８０は、プローブ信号ｚからターゲット信号ｙへの伝達系に、クロススペクトル法による伝達関数推定を適用する。すなわち、平均化されたクロススペクトルＷＺＹを、プローブ信号スペクトルＺのパワースペクトルＷＺＺで割る。伝達関数推定部１８０は、ＷＺＹをＷＺＺで割った値を２乗部１８２に出力する。
【0068】
２乗部１８２は、伝達関数推定部１８０から取得した値を２乗する。乗算部１８４は、２乗部１８２が出力した値と、プローブ信号スペクトルＺのパワースペクトルＷＺＺの積を算出し、算出した値を表示処理部２２に出力する。算出した値は、次に説明するように、ターゲット信号ｙのパワースペクトルに含まれている信号源信号ｘに起因する成分のパワースペクトルＷＶＶである。
【0069】
次に、乗算部１８４が算出した値がパワースペクトルＷＶＶであることを説明する。下記（９）式は、プローブ信号スペクトルＺとターゲット信号スペクトルＹに対してクロススペクトル法を適用した場合の式変形であり、クロススペクトル法に、さらに、前述した（２）式、（３）式を適用している。
（９）式を２乗し、かつＷＺＺを乗じると、（１０）式に示す式変形ができる。（１０）式において、１行目は（９）式の両辺を２乗して、かつＷＺＺを乗じた式であり、２行目は、ＷＺＺによる約分をし、かつ、分母分子にＷＹＹを乗じた式である。
【数9】
【0070】
【数10】
【0071】
コヒーレンス関数算出部１９０は、乗算部１８４が算出した値を、第１平均化処理部７１にて平均化されたターゲット信号スペクトルＹのパワースペクトルＷＹＹで割る。上述したように、乗算部１８４が算出した値はパワースペクトルＷＶＶを意味している。これをパワースペクトルＷＹＹで割った値は、前述した（７）式より、コヒーレンス関数γＺＹ
２になる。コヒーレンス関数算出部１９０は、算出した値を表示処理部２２に出力する
この第２実施形態でも、計算手順は異なるものの、第１実施形態と同じ各種の値が算出できる。よって、第１実施形態と同じ効果が得られる。
この第２実施形態でも、計算手順は異なるものの、第１実施形態と同じ各種の値が算出できる。よって、第１実施形態と同じ効果が得られる。
【0072】
以上、実施形態を説明したが、開示した技術は上述の実施形態に限定されるものではなく、次の変形例も開示した範囲に含まれ、さらに、下記以外にも要旨を逸脱しない範囲内で種々変更して実施できる。
【0073】
＜変形例１＞
実施形態では、コヒーレンス関数γＺＹ ２あるいは積ＷＶＺをプローブ１４の位置に対応付けた画像として、オーバーレイ画像をカメラ画像に重畳したオーバーレイ後画像（図７左下画像、図７右下画像）を表示装置３０に表示していた。しかし、コヒーレンス関数γＺＹ ２あるいは積ＷＶＺをプローブ１４の位置に対応付けた画像として、オーバーレイ画像をカメラ画像に重畳せず、オーバーレイ画像をそのまま表示装置３０に表示してもよい。オーバーレイ画像は、図６に例示した画像である。図６に示す画像からでも、ターゲット信号ｙに影響を与える信号源を容易に推定できる。
実施形態では、コヒーレンス関数γＺＹ ２あるいは積ＷＶＺをプローブ１４の位置に対応付けた画像として、オーバーレイ画像をカメラ画像に重畳したオーバーレイ後画像（図７左下画像、図７右下画像）を表示装置３０に表示していた。しかし、コヒーレンス関数γＺＹ ２あるいは積ＷＶＺをプローブ１４の位置に対応付けた画像として、オーバーレイ画像をカメラ画像に重畳せず、オーバーレイ画像をそのまま表示装置３０に表示してもよい。オーバーレイ画像は、図６に例示した画像である。図６に示す画像からでも、ターゲット信号ｙに影響を与える信号源を容易に推定できる。
【符号の説明】
【0074】
１０：信号源推定システム １１：カメラ １２：ターゲット信号検出線 １４：プローブ １６：ＥＣＵ １８：電子機器 ２０：処理装置 ２１：信号処理部 ２２：表示処理部 ３０：表示装置 ３１：可視化表示領域 ３２：グラフ表示領域 ４０：周波数変換部 ５０：パワースペクトル算出部 ６０：クロススペクトル算出部 ７０：平均化処理部 ８０：コヒーレンス関数算出部 ９０：乗算部 ９２：乗算部 １００：データ読み込み部 １１０：位置検出部 １２０：メモリ １３０：グラフデータ作成部 １４０：オーバーレイデータ作成部 １５０：オーバーレイ処理部 １８０：伝達関数推定部 １８４：乗算部 １９０：コヒーレンス関数算出部 ２２１：信号処理部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【手続補正書】
【提出日】2022-06-06
【手続補正1】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項1】
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
前記検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
前記コヒーレンス関数と、平均化した前記検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、
前記検出子の位置を検出する位置検出部（１１０）と、
前記乗算部が算出した前記積を、前記位置検出部が検出した前記検出子の位置に対応付けて、表示装置（３０）に表示する表示処理部（２２）と、
を備える信号源推定装置。
前記検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
前記コヒーレンス関数と、平均化した前記検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、
前記検出子の位置を検出する位置検出部（１１０）と、
前記乗算部が算出した前記積を、前記位置検出部が検出した前記検出子の位置に対応付けて、表示装置（３０）に表示する表示処理部（２２）と、
を備える信号源推定装置。
【請求項2】
請求項１に記載の信号源推定装置であって、
前記検出子により信号源の位置を推定する範囲を撮影するカメラ（１１）を備え、
前記表示処理部は、前記表示装置に、前記カメラが撮影した画像を視認可能にしつつ、前記カメラが撮影した画像に、前記検出子の二次元位置と前記積との関係を示す画像を重畳表示する、信号源推定装置。
前記検出子により信号源の位置を推定する範囲を撮影するカメラ（１１）を備え、
前記表示処理部は、前記表示装置に、前記カメラが撮影した画像を視認可能にしつつ、前記カメラが撮影した画像に、前記検出子の二次元位置と前記積との関係を示す画像を重畳表示する、信号源推定装置。
【請求項3】
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ
２
）を算出するコヒーレンス関数算出部（８０、１９０）と、
前記検出子の位置を検出する位置検出部（１１０）と、
前記コヒーレンス関数を、前記位置検出部が検出した前記検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、
前記検出子により信号源の位置を推定する範囲を撮影するカメラ（１１）と、
を備え、
前記表示処理部は、前記表示装置に、前記カメラが撮影した画像を視認可能にしつつ、前記カメラが撮影した画像に、前記検出子の二次元位置と前記コヒーレンス関数との関係を示す画像を重畳表示する、信号源推定装置。
前記検出子の位置を検出する位置検出部（１１０）と、
前記コヒーレンス関数を、前記位置検出部が検出した前記検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、
前記検出子により信号源の位置を推定する範囲を撮影するカメラ（１１）と、
を備え、
前記表示処理部は、前記表示装置に、前記カメラが撮影した画像を視認可能にしつつ、前記カメラが撮影した画像に、前記検出子の二次元位置と前記コヒーレンス関数との関係を示す画像を重畳表示する、信号源推定装置。
【請求項4】
解析対象の信号であるターゲット信号（ｙ）と、それぞれ所定位置に存在する複数の信号源を備える対象装置に接近させて順次移動させた検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
前記検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
前記コヒーレンス関数と、平均化した前記検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、
を備える信号源推定装置。
前記検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
前記コヒーレンス関数と、平均化した前記検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、
を備える信号源推定装置。
【請求項5】
解析対象の信号であるターゲット信号（ｙ）と、それぞれ所定位置に存在する複数の信号源を備える対象装置に接近させて順次移動させた検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
前記検出子の位置を検出する位置検出部（１１０）と、
前記コヒーレンス関数を、前記位置検出部が検出した前記検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、
を備える信号源推定装置。
前記検出子の位置を検出する位置検出部（１１０）と、
前記コヒーレンス関数を、前記位置検出部が検出した前記検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、
を備える信号源推定装置。
【手続補正2】
【補正対象書類名】明細書
【補正対象項目名】０００７
【補正方法】変更
【補正の内容】
【0007】
上記目的を達成するための１つの開示は、
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
コヒーレンス関数と、平均化した検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、
検出子の位置を検出する位置検出部（１１０）と、
乗算部が算出した積を、位置検出部が検出した検出子の位置に対応付けて、表示装置（３０）に表示する表示処理部（２２）と、を備える信号源推定装置である。
上記目的を達成するための１つの開示は、
解析対象の信号であるターゲット信号（ｙ）と、それぞれ所定位置に存在する複数の信号源を備える対象装置に接近させて順次移動させた検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ ２ ）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
コヒーレンス関数と、平均化した検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、を備える信号源推定装置である。
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
コヒーレンス関数と、平均化した検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、
検出子の位置を検出する位置検出部（１１０）と、
乗算部が算出した積を、位置検出部が検出した検出子の位置に対応付けて、表示装置（３０）に表示する表示処理部（２２）と、を備える信号源推定装置である。
上記目的を達成するための１つの開示は、
解析対象の信号であるターゲット信号（ｙ）と、それぞれ所定位置に存在する複数の信号源を備える対象装置に接近させて順次移動させた検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ ２ ）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子信号のパワースペクトル（ＷＺＺ）を算出するパワースペクトル算出部（５２）と、
コヒーレンス関数と、平均化した検出子信号のパワースペクトルとの積（ＷＶＺ）を算出し、算出した値を出力する乗算部（９２）と、を備える信号源推定装置である。
【手続補正3】
【補正対象書類名】明細書
【補正対象項目名】０００９
【補正方法】変更
【補正の内容】
【0009】
上記目的を達成するための他の１つの開示は、
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子の位置を検出する位置検出部（１１０）と、
コヒーレンス関数を、位置検出部が検出した検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、検出子により信号源の位置を推定する範囲を撮影するカメラ（１１）と、を備え、
表示処理部は、表示装置に、カメラが撮影した画像を視認可能にしつつ、カメラが撮影した画像に、検出子の二次元位置とコヒーレンス関数との関係を示す画像を重畳表示する、信号源推定装置である。
上記目的を達成するための１つの開示は、
解析対象の信号であるターゲット信号（ｙ）と、それぞれ所定位置に存在する複数の信号源を備える対象装置に接近させて順次移動させた検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ ２ ）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子の位置を検出する位置検出部（１１０）と、
コヒーレンス関数を、位置検出部が検出した検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、を備える信号源推定装置である。
解析対象の信号であるターゲット信号（ｙ）と、検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ２）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子の位置を検出する位置検出部（１１０）と、
コヒーレンス関数を、位置検出部が検出した検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、検出子により信号源の位置を推定する範囲を撮影するカメラ（１１）と、を備え、
表示処理部は、表示装置に、カメラが撮影した画像を視認可能にしつつ、カメラが撮影した画像に、検出子の二次元位置とコヒーレンス関数との関係を示す画像を重畳表示する、信号源推定装置である。
上記目的を達成するための１つの開示は、
解析対象の信号であるターゲット信号（ｙ）と、それぞれ所定位置に存在する複数の信号源を備える対象装置に接近させて順次移動させた検出子（１４）により検出された信号である検出子信号（ｚ）とのコヒーレンス関数（γＺＹ ２ ）を算出するコヒーレンス関数算出部（８０、１９０）と、
検出子の位置を検出する位置検出部（１１０）と、
コヒーレンス関数を、位置検出部が検出した検出子の位置に対応付けて表示装置（３０）に表示する表示処理部（２２）と、を備える信号源推定装置である。