(57)【要約】【課題】より簡便な手法により演算処理を低減しつつ、測定対象から発生する電磁波の波源方位およびエネルギー分布を推定することができる波源方位推定装置、波源方位推定方法及びプログラムを提供すること。【解決手段】測定対象の電子機器から離間した任意の位置において、電子機器（Ｇ）から生じる電磁波を受信するアンテナの受信信号に基づいて位置における電磁波の電磁界多成分を測定して測定値を算出する測定部（１２０）と、測定部により算出された測定値に対して最大エントロピー法を使用した信号処理を行い、信号処理の結果に基づいて電子機器から生じる電磁波の波源の方位およびエネルギー分布を推定する信号処理部（１３０）と、を備える波源方位推定装置（１００）である。【選択図】図１
2019207202.tif 000002
【特許請求の範囲】
【請求項１】
測定対象の電子機器から離間した任意の位置において、前記電子機器から生じる電磁波を受信するアンテナの受信信号に基づいて前記位置における前記電磁波の電磁界多成分を測定して測定値を算出する測定部と、
前記測定部により算出された前記測定値に対して最大エントロピー法を使用した信号処理を行い、信号処理の結果に基づいて前記電子機器から生じる前記電磁波の波源の方位およびエネルギー分布を推定する信号処理部と、を備え、
前記信号処理部は、前記信号処理において前記測定値に基づいて前記電磁波をエネルギー換算するための波動分布関数で表されたスペクトルマトリクスを作成し、前記測定値に基づく制約条件と前記スペクトルマトリクスとの関係に基づいて、エントロピーが最大となる前記波動分布関数を導出し、導出した前記波動分布関数に対して特異値分解を用いた正則化により最適化を行って前記波動分布関数を決定する、
波源方位推定装置。
【請求項２】
前記信号処理部は、前記波動分布関数と積分核との積に基づいて前記スペクトルマトリクスを設定し、
前記波動分布関数を導出する際に、前記測定値に基づいて前記スペクトルマトリクスを離散化して得られる連立方程式の解の空間的滑らかさの基準を示す関数を設定し、
前記関数に基づいて離散化した前記スペクトルマトリクスの平均２乗誤差を最小とする演算を行い前記波動分布関数を導出する、
請求項１に記載の波源方位推定装置。
【請求項３】
前記信号処理部は、前記最適化する処理において前記スペクトルマトリクスの予測値と前記測定値に基づく前記スペクトルマトリクスの値との予測誤差の２乗平均に基づいて正則化処理を行い、前記正則化処理において前記特異値分解を用いて計算回数を低減する、
請求項１または２に記載の波源方位推定装置。
【請求項４】
前記電子機器を撮像する撮像部と、
前記撮像部により撮像された前記電子機器の第１画像を表示する表示部と、
前記信号処理部により推定された前記電磁波の情報を示す第２画像を生成し、前記第２画像を前記電子機器の前記第１画像に重畳した第３画像を前記表示部に表示させる制御部と、を更に備える、
請求項１から３のうちいずれか１項に記載の波源方位推定装置。
【請求項５】
コンピュータが、
測定対象の電子機器から離間した任意の位置において、前記電子機器から生じる電磁波を受信するアンテナの受信信号に基づいて前記位置における前記電磁波の電磁界多成分を測定して測定値を算出し、
算出した前記測定値に対して最大エントロピー法を使用した信号処理を行い、
信号処理の結果に基づいて前記電子機器から生じる前記電磁波の波源の方位およびエネルギー分布を推定し、
前記信号処理において前記測定値に基づいて前記電磁波をエネルギー換算するための波動分布関数で表されたスペクトルマトリクスを作成し、
前記測定値に基づく制約条件と前記スペクトルマトリクスとの関係に基づいて、エントロピーが最大となる前記波動分布関数を導出し、
導出した前記波動分布関数に対して特異値分解を用いた正則化により最適化を行って前記波動分布関数を決定する、
波源方位推定方法。
【請求項６】
コンピュータに、
測定対象の電子機器から離間した任意の位置において、前記電子機器から生じる電磁波を受信するアンテナの受信信号に基づいて前記位置における前記電磁波の電磁界多成分を測定させて測定値を算出させ、
算出させた前記測定値に対して最大エントロピー法を使用した信号処理を行わせ、
信号処理の結果に基づいて前記電子機器から生じる前記電磁波の波源の方位およびエネルギー分布を推定させ、
前記信号処理において前記測定値に基づいて前記電磁波をエネルギー換算するための波動分布関数で表されたスペクトルマトリクスを作成させ、
前記測定値に基づく制約条件と前記スペクトルマトリクスとの関係に基づいて、エントロピーが最大となる前記波動分布関数を導出させ、
導出した前記波動分布関数に対して特異値分解を用いた正則化により最適化を行わせて前記波動分布関数を決定させる、
プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、波源方位推定装置、波源方位推定方法およびプログラムに関する。
【背景技術】
【０００２】
近年、電子機器は、電磁波を発生するものが増加している。このような電子機器から発生する電磁波は、他の電子機器に対する妨害電波の発生源となり得る。電子機器からの放射妨害波による他の電子機器の誤動作を防ぐには、電子機器から放射される妨害波の波源の位置を効率よく推定してＥＭＩ（Electromagnetic Interference）対策を検討する必要がある。
【０００３】
電子機器から発生する電磁波の波源の位置を推定する方法として、例えば、特許文献１および特許文献２に記載されたものが知られている。特許文献１に記載された技術は、測定対象の近傍の観測点で測定した磁界分布の情報を用いて、電磁波発生源の各方向成分を分離抽出して求め、測定対象内の電磁波発生源を検出し、該検出した電磁波発生源に基づいて遠方における電磁界強度を算出するものである。
【０００４】
特許文献２に記載された技術は、測定対象から離間した一点において電磁波を測定し、測定値に対してＭＵＳＩＣ法による信号処理を行うことにより妨害波の波源方位を推定している。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平１１−１４２４５３号公報
【特許文献２】特開２０１４−２０６４４０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１に記載された技術によれば、電子機器の近傍で近傍電磁波を測定する際に多数の測定点を設定して測定しなければならず、放射妨害波となる電磁波の波源方位の推定に処理時間を要する。
【０００７】
特許文献２に記載された技術によれば、MUSIC法では波源の方位しか推定できず、エネルギー分布の推定ができないという課題がある。
【０００８】
本発明は、上述した課題を解決する波源方位推定装置、波源方位推定方法及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明の第１の態様によれば、測定対象の電子機器から離間した任意の位置において、前記電子機器から生じる電磁波を受信するアンテナの受信信号に基づいて前記位置における前記電磁波の電磁界多成分を測定して測定値を算出する測定部と、前記測定部により算出された前記測定値に対して最大エントロピー法を使用した信号処理を行い、信号処理の結果に基づいて前記電子機器から生じる前記電磁波の波源の方位およびエネルギー分布を推定する信号処理部と、を備え、前記信号処理部は、前記信号処理において前記測定値に基づいて前記電磁波をエネルギー換算するための波動分布関数で表されたスペクトルマトリクスを作成し、前記測定値に基づく制約条件と前記スペクトルマトリクスとの関係に基づいて、エントロピーが最大となる前記波動分布関数を導出し、導出した前記波動分布関数に対して特異値分解を用いた正則化により最適化を行って前記波動分布関数を決定する、波源方位推定装置である。
【００１０】
また本発明は、コンピュータが、測定対象の電子機器から離間した任意の位置において、前記電子機器から生じる電磁波を受信するアンテナの受信信号に基づいて前記位置における前記電磁波の電磁界多成分を測定して測定値を算出し、算出した前記測定値に対して最大エントロピー法を使用した信号処理を行い、信号処理の結果に基づいて前記電子機器から生じる前記電磁波の波源の方位およびエネルギー分布を推定し、前記信号処理において前記測定値に基づいて前記電磁波をエネルギー換算するための波動分布関数で表されたスペクトルマトリクスを作成し、前記測定値に基づく制約条件と前記スペクトルマトリクスとの関係に基づいて、エントロピーが最大となる前記波動分布関数を導出し、導出した前記波動分布関数に対して特異値分解を用いた正則化により最適化を行って前記波動分布関数を決定する、波源方位推定方法である。
【００１１】
また本発明は、コンピュータに、測定対象の電子機器から離間した任意の位置において、前記電子機器から生じる電磁波を受信するアンテナの受信信号に基づいて前記位置における前記電磁波の電磁界多成分を測定させて測定値を算出させ、算出させた前記測定値に対して最大エントロピー法を使用した信号処理を行わせ、信号処理の結果に基づいて前記電子機器から生じる前記電磁波の波源の方位およびエネルギー分布を推定させ、前記信号処理において前記測定値に基づいて前記電磁波をエネルギー換算するための波動分布関数で表されたスペクトルマトリクスを作成させ、前記測定値に基づく制約条件と前記スペクトルマトリクスとの関係に基づいて、エントロピーが最大となる前記波動分布関数を導出させ、導出した前記波動分布関数に対して特異値分解を用いた正則化により最適化を行わせて前記波動分布関数を決定させる、プログラムである。
【発明の効果】
【００１２】
本発明によれば、より簡便な手法により演算処理を低減しつつ、測定対象から発生する電磁波の波源方位およびエネルギー分布を推定することができる。
【図面の簡単な説明】
【００１３】
【図１】第１実施形態の波源方位推定装置の構成を示すブロック図である。
【図２】測定点から見た一つの波源の波源方位を示す図である。
【図３】波源方位推定装置１００において実行される処理の流れの一例を示すフローチャートである。
【図４】第２実施形態の波源方位推定装置１００Ａの構成を示すブロック図である。
【図５】表示部１４０に表示された第３画像ＩＭの内容の一例を示す図である。
【発明を実施するための形態】
【００１４】
以下、図面を参照にしつつ、本発明の一実施形態による波源方位推定装置を説明する。
【００１５】
［第１実施形態］
図１は第１実施形態の波源方位推定装置の構成を示すブロック図である。波源方位推定装置１００は、測定対象の電子機器Ｇから発生する放射妨害波（電磁波）の波源の方位を推定するものである。波源方位推定装置１００は、例えば、アンテナ１１０と、測定部１２０と、信号処理部１３０とを備える。
【００１６】
アンテナ１１０は、空間を伝搬する電磁波（放射妨害波）を受信する。アンテナ１１０は、測定対象の電子機器Ｇから離間した一点において受信した電磁波のうち電磁界の多成分を測定する。アンテナ１１０は、例えば、測定点における電磁界多成分を検出する電磁界プローブである。電磁界多成分とは、例えば、電界の三軸方向の成分（Ｅｘ、Ｅｙ、Ｅｚ）である。測定において電子機器Ｇから発生する電磁波の波源は一つとは限らず、波源が複数個存在していてもよい。波源が複数個ある場合、アンテナ１１０で測定される電界の測定値は、複数個の波源から発生する電磁波の電界成分が合成された値となる。
【００１７】
アンテナ１１０は、電子機器Ｇから離間した位置に設置される。アンテナ１１０は、電子機器Ｇを稼働させて電磁波が発生している状態で、電子機器Ｅから離間した一点において電磁波を受信し、電磁波における三軸方向の電界成分を測定する。関連する技術においては、受信アンテナで電界または磁界の一成分のみを測定していたが、波源方位推定装置１００においては、電界の三軸方向の成分を測定値として扱う。
【００１８】
アンテナ１１０は、例えば、波源方位推定装置１００と別体に設けられている。アンテナ１１０は、例えば、波源方位推定装置１００の本体部に設けられていてもよい。アンテナ１１０は、測定部１２０と電気的に接続されている。アンテナ１１０は、例えば、測定部１２０と受信信号をＡＤ変換する増幅回路等を介して電気的に接続される。アンテナ１１０は、増幅された測定データを測定部１２０に出力する。アンテナ１１０は、測定データを数値データに変換して測定部１２０に出力してもよい。
【００１９】
測定部１２０は、例えば、アンテナ１１０から出力された測定データを取得し、取得した測定データに基づいて、電界の振幅と位相を演算し、測定値として出力する。測定部１２０は、測定値を信号処理部１３０に出力する。
【００２０】
信号処理部１３０は、例えば、測定部１２０から取得した測定値に対して信号処理を行って後述の波動分布関数を決定し、決定した波動分布関数に基づいて電子機器Ｇから発生した電磁波の波源の方位を推定する。信号処理部１３０は、例えば、最大エントロピー法を使用して波動分布関数を導出する。
【００２１】
最大エントロピー法とは、「不確かさ」を図る尺度であるエントロピーを既知の所定条件の下で最大にするよう分布を選ぶべきであるという原理に基づいて、測定値などに基づいて所定条件が与えられた場合の確率変数の分布を推定するものである。信号処理部１３０は、決定した波動分布関数に基づいて、受信した電磁波の波源方位（θ,φ）（図２参照）を推定する。波源方位の推定方法については、後に詳しく説明する。
【００２２】
［動作の説明］
次に、信号処理部１３０において実行される最大エントロピー法に基づく波源方位の推定処理について説明する。
【００２３】
図２は、測定点から見た一つの波源の波源方位を示す図である。波源方位は、測定点から見た電磁波の波源の三次元位置を角度で表したものである。波源方位は、ｘｙ平面における測定点から波源へのベクトルとのなす角を方位角φとし、ｚ軸と測定点から波源へのベクトルとのなす角度を入射角θとして表される。
【００２４】
信号処理部１３０は、測定部１２０から取得した電磁波の電界の三軸成分の測定値に基づいて、エントロピーが最大となる波動分布関数Ｆを決定する。
【００２５】
ここで、波動分布関数Ｆが未知の状態は、高いエントロピーを持っている状態であることを意味する。波動分布関数Ｆが未知の状態から波動分布関数Ｆを決定するという処理は、波動分布関数Ｆが未知の状態における高いエントロピーを下げる処理である。従って、エントロピーは、得られた測定値に基づく制約条件がある場合、制約条件下で最も効率よく下げられるような波動分布関数Ｆが決定される。すなわち、決定された波動分布関数Ｆは、制約条件下でエントロピーが最大となる分布を解として求められる。このようにして得られた解は、制約条件下で予測される最も「ありふれた」分布を表していると解釈できる。
【００２６】
信号処理部１３０は、例えば、以下の式（１）から式（６）の演算によりエントロピーを最大化する演算を行って波動分布関数を導出する。信号処理部１３０は、例えば、以下の式（７）から式（９）の演算により、導出した波動分布関数を正則化処理により最適化し、波動分布関数を決定する。
【００２７】
信号処理部１３０は、時間領域での測定値である電界の測定値に対してフーリエ変換等の信号処理を行って周波数領域のデータに変換する。時間領域とは、時間を独立変数に持つ空間である。周波数領域とは、周波数を独立変数に持つ空間である。
【００２８】
信号処理部１３０は、変換した周波数領域のデータを用いて、以下の式（１）で示されるスペクトルマトリクスＳを生成する。スペクトルマトリクスＳは、測定された電界値をエネルギー換算する行列式である。信号処理部１３０は、生成したスペクトルマトリクスＳに基づいて、測定された電磁波の周波数とエネルギーとの関係を求める。
【００２９】
ここで、波源から電磁波（平面波）が角周波数ω、入射角θ、方位角φで測定点に伝搬する場合、スペクトルマトリクスＳの各(i, j) 成分 (i, j = 1, 2, …, 6) は、例えば、波動分布関数Ｆを用いて式（1）のように表わされる。
【００３０】
【数１】
2019207202.tif 000003
【００３１】
式（１）において、aij は、単位立体角あたりのエネルギーを表す積分核を示す行列である。Ｆは後述の波動分布関数を示すＫ次縦ベクトルである。信号処理部１３０は、式（１）における電子機器Ｇの表面におけるエネルギー分布を表す波動分布関数Ｆ（ω、θ、φ）を求める。
【００３２】
図２に示されるように、電子機器Ｇの表面の方形領域を像再生領域Ｄと定義する。像再生とは、電磁波の波動領域において数値的に像を再生することである。信号処理部１３０は、像再生領域Ｄにおいて、x軸方向の辺をKx個に等分割すると共に、y軸方向の辺をKy個に等分割し、像再生領域ＤをK=KxKy個の等方形画素に分割する。ここで、分割した各画素の中央における像の値を並べて得られるＫ次縦ベクトルをＦとする。信号処理部１３０は、式（１）を数値計算により解くために、式（１）の積分核について関連した離散化を行い、式（２）を得る。
【００３３】
【数２】
2019207202.tif 000004
【００３４】
ここで、SはS(ω) の9個の実数成分Si(ω)の測定値を要素とするＭ次縦ベクトル（Ｍ＝９）である。Ａは離散化された積分核を要素とするＭ×Ｋ行列である。ｅはＳの各要素の測定誤差や積分離散化による誤差を代表するものとして導入されたＭ次縦ベクトルである。Ｍ＝９に対して通常、Ｋは非常に大きい（Ｋ＞＞Ｍ）ため、未知数の数よりも多くの連立方程式が存在し、解は一意に決まらない。
【００３５】
そこで、信号処理部１３０は、解の空間的滑らかさの基準となる関数を設定し、設定した関数に基づいて解を一意に決定する。空間的な滑らかさとは、解の関数が連続的に微分可能であることをいい、解の関数に導関数が存在し、且つ、解の関数が連続関数となることである。信号処理部１３０は、解の滑らかさを表す関数として、式（３）に示される分布関数の（ｘ，ｙ）平面における空間２階微分の像再生領域Ｄにわたる２乗積分φを用いる。
【００３６】
【数３】
2019207202.tif 000005
【００３７】
また、信号処理部１３０は、式（４）に示されるように、測定値に基づくスペクトルマトリクスと波動分布関数と離散化した積分核との積との平均２乗誤差が一定値以下という条件の下で、式（４）における解を最小化する演算を行う。
【００３８】
【数４】
2019207202.tif 000006
【００３９】
ここで、ε2は誤差限界であるが、実際の計算では、この値があらかじめ決められている必要はない。式（４）の解を得るために、式（５）の量Φ(F)を適用する。
【００４０】
【数５】
2019207202.tif 000007
【００４１】
式（５）において、γはラグランジュ未定定数の逆数であり、Ｃはラプラシアンを離散化したＫ次正方行列であり、|CF|2は空間２階微分の差分近似を用いた式（３）のφの離散表現である。Φの最小化のためにFの要素Fiで微分して、∂Φ/∂Fi=0 (i=1, 2,…, K)とおき、その解をF^( γ )とすると、式（６）が得られる。
【００４２】
【数６】
2019207202.tif 000008
【００４３】
式（６）で得られた波動分布関数F^( γ )は、エントロピーを最大化するものである。信号処理部１３０は、更に、式（６）で導出した波動分布関数F^( γ )を最適化する。信号処理部１３０は、例えば、式（７）を用いて正則化によりF^( γ )を最適化する。
【００４４】
F^( γ )から得られる再生スペクトルγは解の正則化の度合いを左右するパラメータである。信号処理部１３０は、最適なγ値を客観的に評価するため、スペクトルマトリクスSの第i要素Siを取り除き、残りの要素を使って再生されるスペクトルマトリクスSの第i要素の予測値[AF^(i)( γ )]iと、実際のSiとの予測誤差を比較する。i=1,2 …,Mにわたる予測誤差の２乗平均P(γ)は、式（７）で表される。
【００４５】
【数７】
2019207202.tif 000009
【００４６】
式（７）をそのまま適用すると計算時間が長くなるため、信号処理部１３０は、特異値分解により計算回数を低減させ解F^(γ)を最適化する。式（７）において、AC-1の特異値分解を利用すると、予測誤差G(γ)は式（８）で表される。
【００４７】
【数８】
2019207202.tif 000010
【００４８】
ここで、σk(k=1,2 …,M)は特異値、ZiはベクトルZ=UTSの第i要素、Uは特異値分解AC-1=UΣVTにおけるM次直交行列である。なお、V はK次直交行列、Σはσkを対角要素とするM×K行列であり、K≧Mが仮定されている。更に、式（６）に同じAC-1の特異値分解を利用すると、F^(γ)は式（９）に書き直される。
【００４９】
【数９】
2019207202.tif 000011
【００５０】
F^k(γ)はF^(γ)の第k要素、v’kjはC-1Vの(k,j)要素である。信号処理部１３０は、上記の式（７）から式（９）の計算処理を行うことにより、最適化した波動分布関数Ｆを決定する。式（９）で得られる波動分布関数F^k(γ)は、電子機器Ｇから発生する電磁波のエネルギー分布を示す。電子機器Ｇにおける電磁波の波源の位置は、波動分布関数Fにおける極値として表される。電子機器Ｇに波源が複数個ある場合、波源は、波動分布関数において複数の極大値として出現する位置にあると推定される。
【００５１】
信号処理部１３０は、決定された波動分布関数Ｆを用いて電子機器Ｇから発生する電磁波のエネルギー分布を推定する。また、信号処理部１３０は、決定した波動分布関数により得られた極値に基づいて、測定点から極値に対するベクトルの方位角φと入射角θを算出し、波源の位置を推定する。
【００５２】
次に、波源方位推定装置１００において実行される処理の流れについて説明する。図３は、波源方位推定装置１００において実行される処理の流れの一例を示すフローチャートである。
【００５３】
電子機器Ｇから離間した任意の位置において、アンテナ１１０により電子機器Ｇから生じる電磁波の電界の三軸成分を受信する（ステップＳ１００）。測定部１２０は、アンテナ１１０から受信した電磁波の受信信号に基づいて測定値として出力する（ステップＳ１０２）。信号処理部１３０は、測定部１２０から取得した測定値を信号処理し、式（１）により波動分布関数に基づくスペクトルマトリクスを生成する（ステップＳ１０４）。
【００５４】
信号処理部１３０は、測定値に基づいてエントロピーが最大となるように波動分布関数の導出処理を開始する（ステップＳ１０６）。信号処理部１３０は、式（１）のスペクトルマトリクスを離散化して式（２）に基づいて連立方程式を導出する（ステップＳ１０８）。信号処理部１３０は、連立方程式の解を決定するための関数を式（３）により設定する（ステップＳ１１０）。信号処理部１３０は、設定した関数に基づいて、式（４）及び式（５）により解を最小化する演算を行う（ステップＳ１１２）。
【００５５】
信号処理部１３０は、演算に基づいて、式（６）に示される波動分布関数を導出する（ステップＳ１１４）。信号処理部１３０は、式（７）を用いて正則化により導出した波動分布関数を最適化する（ステップＳ１１６）。信号処理部１３０は、最適化において式（８）に示される特異値分解により最適化の演算を行い、計算回数を低減し式（９）に示される波動分布関数を決定する（ステップＳ１１８）。信号処理部１３０は、決定した波動分布関数に基づいて、波源位置（波源方位）とエネルギー分布を推定する（ステップＳ１２０）。
【００５６】
上述した第１実施形態によれば、波源方位推定装置１００は、電子機器Ｇから離間した一点の位置において電磁波を１回測定することにより、演算処理を低減しつつ、測定対象から発生する電磁波の波源方位およびエネルギー分布を推定することができる。波源方位推定装置１００は、波動分布関数の決定において演算処理を低減することができる。
【００５７】
［第２実施形態］
第１実施形態の波源方位推定装置１００は、波動分布関数を設定して電磁波のエネルギー分布等を推定するものであった。第２実施形態の波源方位推定装置１００Ａは、電子機器を撮像した画像を表示すると共に、推定した電磁波のエネルギー分布のグラフを電子機器の画像に重畳して表示するものである。以下の説明では、第１実施形態と同様の構成については、同一の名称及び符号を用い、重複する説明については適宜省略する。
【００５８】
図４は、第２実施形態の波源方位推定装置１００Ａの構成を示すブロック図である。波源方位推定装置１００Ａは、例えば、アンテナ１１０と、測定部１２０と、信号処理部１３０と、表示部１４０と、操作部１５０と、撮像部１６０と、制御部１７０と、記憶部１８０とを備える。波源方位推定装置１００Ａは、例えば、ＰＣ（Personal Computer）、タブレット型端末等により実現される。
【００５９】
表示部１４０は、例えば、液晶ディスプレイ等の表示装置である。表示部１４０は、タッチパネルであってもよい。表示部１４０は、波源方位推定装置１００Ａの本体部と別体に設けられていてもよい。表示部１４０は、制御部１７０の制御により、画像、映像等を表示面に表示する。
【００６０】
操作部１５０は、ユーザの操作を受け付けるキーボードなどの入力装置である。操作部は、タッチパネル等により実現された表示部１４０に表示されるＧＵＩボタンにより構成されていてもよい。
【００６１】
撮像部１６０は、例えば、アンテナ１１０の位置から計測対象の電子機器Ｇを撮像するカメラ等の撮像装置である。撮像部１６０は、撮像した画像データを記憶部１８０に記憶する。撮像部１６０は、波源方位推定装置１００Ａと別体に設けられているカメラ、その他の撮像機能を備える端末装置であってもよい。
【００６２】
図５は、表示部１４０に表示された第３画像ＩＭの内容の一例を示す図である。制御部１７０は、例えば、操作部１５０により受け付けられたユーザの操作に基づいて、記憶部１８０に記憶された電子機器Ｇの第１画像を表示部１４０に表示させる。
【００６３】
制御部１７０は、例えば、操作部１５０により受け付けられたユーザの操作に基づいて、信号処理部１３０により推定された電子機器Ｇから発生した電磁波の情報を示す第２画像を生成する。エネルギー分布の情報を示す画像は、例えば、三次元で示されたエネルギー分布を示すグラフＲや、電磁波の波源Ｗの位置を示す画像や、電界強度を色別で示す画像等を含む。制御部１７０は、生成した第２画像を第１画像に重畳した第３画像ＩＭを生成し、表示部に表示させる。
【００６４】
記憶部１８０は、例えば、撮像部１６０で撮像された第１画像を記憶する。記憶部１８０は、例えば、信号処理部１３０において決定された波動分布関数を記憶する。記憶部１８０は、例えば、制御部１７０により生成された第２画像等を保存する。記憶部１８０は、例えば、フラッシュメモリやＨＤＤ（Hard Disk Drive）等の記憶装置である。
【００６５】
上述した第２実施形態によれば、表示部１４０に表示された第３画像ＩＭを確認することにより、ユーザは電子機器Ｇにおける波源Ｗの位置と、波源Ｗから発生する電磁波のエネルギー分布を直感的に把握することができる。これによりユーザは、電子機器Ｇのどの部分にＥＭＩ対策を行えばよいか、検討することができる。
【００６６】
上述の波源方位推定装置１００は内部に、コンピュータシステムを有している。そして、上述した波源方位推定する処理の過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。
【００６７】
なお、図１における処理部の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより波源方位推定する処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータシステム」は、ホームページ提供環境（あるいは表示環境）を備えたＷＷＷシステムも含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（ＲＡＭ）のように、一定時間プログラムを保持しているものも含むものとする。
【００６８】
また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【００６９】
以上、この発明の実施形態を図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。また、上記実施形態が有する１または複数のＣＰＵ等のコンピュータが実行するプログラムの一部または全部は、通信回線やコンピュータ読取可能な記録媒体を介して頒布することができる。
【符号の説明】
【００７０】
１００、１００Ａ…波源方位推定装置
１１０…アンテナ
１２０…測定部
１３０…信号処理部
１４０…表示部
１５０…操作部
１６０…撮像部
１７０…制御部
１８０…記憶部
【図１】
2019207202.tif 000012
【図２】
2019207202.tif 000013
【図３】
2019207202.tif 000014
【図４】
2019207202.tif 000015
【図５】
2019207202.tif 000016