JPH09134113A - 波源像可視化方法及び装置 - Google Patents
波源像可視化方法及び装置Info
- Publication number
- JPH09134113A JPH09134113A JP28984895A JP28984895A JPH09134113A JP H09134113 A JPH09134113 A JP H09134113A JP 28984895 A JP28984895 A JP 28984895A JP 28984895 A JP28984895 A JP 28984895A JP H09134113 A JPH09134113 A JP H09134113A
- Authority
- JP
- Japan
- Prior art keywords
- spectrum
- source image
- signal
- observation
- wave
- 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.)
- Granted
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Holo Graphy (AREA)
Abstract
(57)【要約】
【課題】音波や電波ホログラムの原理に従い、２次元干
渉測定を行って波源像の可視化を行う場合に、コヒーレ
ンシに関する事前の知識を必要とせず、部分コヒーレン
トな波源分布の評価が可能であって、スペクトルの速い
変化に追従できるようにする。 【解決手段】固定センサ１４からの信号ａ(ｔ)と走査セ
ンサ１３からの信号ｂ(ｔ)をそれぞれフーリエ変換して
スペクトルＳa(ω),Ｓa(ω)を得る。クロススペクトル
の平均値〈Ｓa *(ω)Ｓb(ω)〉や各パワースペクトルの
平方根の平均値〈｜Ｓa(ω)｜〉,〈｜Ｓb(ω)｜〉を算
出し、これらに基づいてコヒーレンス関数γ2(ω)を定
める。クロススペクトルを｜Ｓa(ω)｜で規格化し、さ
らにコヒーレンス関数γ2(ω)に応じて観測周波数帯域
で積分して相関値Ｃab(ｘ,ｙ)を求める。観測対象１１
への距離ｚに応じて、波源像の再生処理の演算方法を変
える。
渉測定を行って波源像の可視化を行う場合に、コヒーレ
ンシに関する事前の知識を必要とせず、部分コヒーレン
トな波源分布の評価が可能であって、スペクトルの速い
変化に追従できるようにする。 【解決手段】固定センサ１４からの信号ａ(ｔ)と走査セ
ンサ１３からの信号ｂ(ｔ)をそれぞれフーリエ変換して
スペクトルＳa(ω),Ｓa(ω)を得る。クロススペクトル
の平均値〈Ｓa *(ω)Ｓb(ω)〉や各パワースペクトルの
平方根の平均値〈｜Ｓa(ω)｜〉,〈｜Ｓb(ω)｜〉を算
出し、これらに基づいてコヒーレンス関数γ2(ω)を定
める。クロススペクトルを｜Ｓa(ω)｜で規格化し、さ
らにコヒーレンス関数γ2(ω)に応じて観測周波数帯域
で積分して相関値Ｃab(ｘ,ｙ)を求める。観測対象１１
への距離ｚに応じて、波源像の再生処理の演算方法を変
える。
Description
【０００１】
【発明の属する技術分野】本発明は、ホログラフィの原
理に基づいて波源像を可視化する方法及び装置に関す
る。
理に基づいて波源像を可視化する方法及び装置に関す
る。
【０００２】
【従来の技術】電波ホログラムや音波ホログラムを用い
て波源像を可視化する技術が実用化され、例えば不要電
磁波輻射の低減や騒音低減などのために使用されてい
る。観測対象の物体の電波ホログラムを得て解析するこ
とにより、その物体のどこからどの程度の不要輻射が生
じているかを知ることができ（輻射源分布の可視化）、
いわゆるＥＭＩ(Electromagnetic Interference)対策を
進める上で有効な情報を得ることができる。音波ホログ
ラフィによる音源分布の可視化については、例えば、上
羽 貞行、”音波ホログラフィと騒音源探査”、計測と
制御、第１６巻第５号、４２７〜４３３頁、１９７７年
に記載がある。
て波源像を可視化する技術が実用化され、例えば不要電
磁波輻射の低減や騒音低減などのために使用されてい
る。観測対象の物体の電波ホログラムを得て解析するこ
とにより、その物体のどこからどの程度の不要輻射が生
じているかを知ることができ（輻射源分布の可視化）、
いわゆるＥＭＩ(Electromagnetic Interference)対策を
進める上で有効な情報を得ることができる。音波ホログ
ラフィによる音源分布の可視化については、例えば、上
羽 貞行、”音波ホログラフィと騒音源探査”、計測と
制御、第１６巻第５号、４２７〜４３３頁、１９７７年
に記載がある。
【０００３】図２は、２次元干渉観測を行って２次元複
素インタフェログラム（複素ホログラム）を求めこのイ
ンタフェログラムを再生することによって波源分布の表
示を行う従来の波源像可視化装置の構成を示すブロック
図であり、上述の参考文献に示されているものである。
素インタフェログラム（複素ホログラム）を求めこのイ
ンタフェログラムを再生することによって波源分布の表
示を行う従来の波源像可視化装置の構成を示すブロック
図であり、上述の参考文献に示されているものである。
【０００４】観測対象５１に対して距離ｚ0だけ離れた
場所に走査観測面５２を設定し、走査観測面５２内を移
動する走査センサ５３を用い、走査観測面５２の各点に
おいて、観測対象５１からの電磁波や音波を検出して信
号ｂ(ｔ)を得る。また、走査センサ５３とは別に固定セ
ンサ５４を設け、固定センサ５４でも観測対象５１から
の電磁波や音波を検出して信号ａ(ｔ)を得る。なお、ｔ
は時刻を表わし、ｘ,ｙは走査観測面５２内での２次元
座標を表わし、ｘ',ｙ'は観測対象５１での２次元座標
を表わしている。
場所に走査観測面５２を設定し、走査観測面５２内を移
動する走査センサ５３を用い、走査観測面５２の各点に
おいて、観測対象５１からの電磁波や音波を検出して信
号ｂ(ｔ)を得る。また、走査センサ５３とは別に固定セ
ンサ５４を設け、固定センサ５４でも観測対象５１から
の電磁波や音波を検出して信号ａ(ｔ)を得る。なお、ｔ
は時刻を表わし、ｘ,ｙは走査観測面５２内での２次元
座標を表わし、ｘ',ｙ'は観測対象５１での２次元座標
を表わしている。
【０００５】各センサ５３,５４からの信号ｂ(ｔ),ａ
(ｔ)は、バンドパスフィルタ５５,５６で所定の周波数
帯域に制限され、信号ｂ(ｔ)はロックインアンプ５８,
５９に入力し、信号ａ(ｔ)は位相シフタ５７に入力す
る。位相シフタ５７からは、位相遅れのない信号（０
°）が一方のロックインアンプ５８の参照信号として出
力し、π／４だけ遅れた信号（９０°）が他方のロック
インアンプ５８,５９の参照信号として出力する。した
がって、一方のロックインアンプ５８から、走査観測面
５２上の点(ｘ,ｙ)での相関値の実数成分ＪRe(ｘ,ｙ)が
出力され、他方のロックインアンプ５９からは点(ｘ,
ｙ)での相関値の虚数成分ＪIm(ｘ,ｙ)が出力されること
になる。走査観測面５２の各点での相関値を得ることに
より、２次元複素インタフェログラムが得られることに
なる。
(ｔ)は、バンドパスフィルタ５５,５６で所定の周波数
帯域に制限され、信号ｂ(ｔ)はロックインアンプ５８,
５９に入力し、信号ａ(ｔ)は位相シフタ５７に入力す
る。位相シフタ５７からは、位相遅れのない信号（０
°）が一方のロックインアンプ５８の参照信号として出
力し、π／４だけ遅れた信号（９０°）が他方のロック
インアンプ５８,５９の参照信号として出力する。した
がって、一方のロックインアンプ５８から、走査観測面
５２上の点(ｘ,ｙ)での相関値の実数成分ＪRe(ｘ,ｙ)が
出力され、他方のロックインアンプ５９からは点(ｘ,
ｙ)での相関値の虚数成分ＪIm(ｘ,ｙ)が出力されること
になる。走査観測面５２の各点での相関値を得ることに
より、２次元複素インタフェログラムが得られることに
なる。
【０００６】このようにして得られた各点での相関値に
基づき波源像再生処理部６０により波形像を再生し、こ
の波形像を表示部６１に表示する。具体的には、観測対
象５１上の点(ｘ',ｙ')における輻射の強度をＩ(ｘ',
ｙ')で表わすと、この強度Ｉ(ｘ',ｙ')は、以下のよう
に表わされる。ｊは虚数単位である。
基づき波源像再生処理部６０により波形像を再生し、こ
の波形像を表示部６１に表示する。具体的には、観測対
象５１上の点(ｘ',ｙ')における輻射の強度をＩ(ｘ',
ｙ')で表わすと、この強度Ｉ(ｘ',ｙ')は、以下のよう
に表わされる。ｊは虚数単位である。
【０００７】
【数１】 波源像再生処理部６０は式(1)に表わされる計算を実行
する。
する。
【０００８】ところで、式(1)の右辺は、波源がコヒー
レントである場合には振幅を表わすが、波源がインコヒ
ーレントである場合には振幅の平方を表わす。したがっ
て、表示部６１において観測対象５１における振幅分布
を表示させる場合、波源がコヒーレントである場合には
Ｉ(ｘ',ｙ')を表示部６１に対して出力し、波源がイン
コヒーレントの場合には
レントである場合には振幅を表わすが、波源がインコヒ
ーレントである場合には振幅の平方を表わす。したがっ
て、表示部６１において観測対象５１における振幅分布
を表示させる場合、波源がコヒーレントである場合には
Ｉ(ｘ',ｙ')を表示部６１に対して出力し、波源がイン
コヒーレントの場合には
【０００９】
【外１】 を表示部６１に対して出力する必要がある。
【００１０】
【発明が解決しようとする課題】上述した従来の波源像
可視化装置には、ロックインアンプの応答性の制限のた
めに、速いスペクトル変動に追従できないという問題点
がある。また、上述したように、波源がコヒーレントか
インコヒーレントであるかによって再生像の解釈が異な
り、事前にコヒーレント波源であるかインコヒーレント
波源であるかを知らないと、再生像の解釈を行うことが
できないこともある。電子装置からの不要輻射の可視化
の場合、時間、空間的に、コヒーレント、インコヒーレ
ント、部分コヒーレントな成分が混在し、しかも事前に
その割合等を予測することが困難なことが多いので、結
局、従来の装置によっては定量的な解釈を行うことが難
しい。さらに、従来の装置には、強度分布のみが求めら
れるので、波源状態の時間、空間的なコヒーレンシを評
価できないという問題点もある。
可視化装置には、ロックインアンプの応答性の制限のた
めに、速いスペクトル変動に追従できないという問題点
がある。また、上述したように、波源がコヒーレントか
インコヒーレントであるかによって再生像の解釈が異な
り、事前にコヒーレント波源であるかインコヒーレント
波源であるかを知らないと、再生像の解釈を行うことが
できないこともある。電子装置からの不要輻射の可視化
の場合、時間、空間的に、コヒーレント、インコヒーレ
ント、部分コヒーレントな成分が混在し、しかも事前に
その割合等を予測することが困難なことが多いので、結
局、従来の装置によっては定量的な解釈を行うことが難
しい。さらに、従来の装置には、強度分布のみが求めら
れるので、波源状態の時間、空間的なコヒーレンシを評
価できないという問題点もある。
【００１１】本発明の目的は、コヒーレンシに関する事
前の知識を必要とせず、部分コヒーレントな波源分布の
評価が可能であって、スペクトルの速い変化に追従でき
る波源可視化方法及び装置を提供することにある。
前の知識を必要とせず、部分コヒーレントな波源分布の
評価が可能であって、スペクトルの速い変化に追従でき
る波源可視化方法及び装置を提供することにある。
【００１２】
【課題を解決するための手段】本発明の波源可視化方法
は、観測対象からの放射波動に対して２次元干渉測定を
実施し、観測対象における波源像を可視化する方法であ
って、２次元干渉測定の観測面に対する基準となる位置
での放射波動を検出して第１の信号として出力する第１
のセンサと、観測面内の任意の位置での放射波動を検出
して第２の信号として出力する第２のセンサとを使用
し、第１の信号及び第２の信号をそれぞれフーリエ変換
して第１のスペクトル及び第２のスペクトルを獲得し、
第１のスペクトルのパワースペクトルと、第２のスペク
トルのパワースペクトルと、第１のスペクトルと第２の
スペクトルの間のクロススペクトルとを求め、クロスス
ペクトルの平均値と各パワースペクトルから算出された
平均値とに基づきコヒーレンス関数を決定し、観測面内
の点ごとに、クロススペクトルを得てこれを第１のスペ
クトルのパワースペクトルの平方根で規格化し、規格化
されたクロススペクトルに対して観測周波数帯域内でコ
ヒーレンス関数の値に応じた積分演算を実行して相関値
を算出し、観測対象と観測面との距離に応じて波源像の
再生のための演算処理方法を選択し、各点ごとの相関値
に基づいて波源像を再生する。
は、観測対象からの放射波動に対して２次元干渉測定を
実施し、観測対象における波源像を可視化する方法であ
って、２次元干渉測定の観測面に対する基準となる位置
での放射波動を検出して第１の信号として出力する第１
のセンサと、観測面内の任意の位置での放射波動を検出
して第２の信号として出力する第２のセンサとを使用
し、第１の信号及び第２の信号をそれぞれフーリエ変換
して第１のスペクトル及び第２のスペクトルを獲得し、
第１のスペクトルのパワースペクトルと、第２のスペク
トルのパワースペクトルと、第１のスペクトルと第２の
スペクトルの間のクロススペクトルとを求め、クロスス
ペクトルの平均値と各パワースペクトルから算出された
平均値とに基づきコヒーレンス関数を決定し、観測面内
の点ごとに、クロススペクトルを得てこれを第１のスペ
クトルのパワースペクトルの平方根で規格化し、規格化
されたクロススペクトルに対して観測周波数帯域内でコ
ヒーレンス関数の値に応じた積分演算を実行して相関値
を算出し、観測対象と観測面との距離に応じて波源像の
再生のための演算処理方法を選択し、各点ごとの相関値
に基づいて波源像を再生する。
【００１３】本発明において、放射波動としては、音波
や電波など、２次元干渉測定が可能な任意の波動を用い
ることができる。マイクロ波帯の電波の波源像を求める
場合には、いったん、受信信号を中間周波数帯の信号に
周波数変換してからフーリエ変換を行うことも可能であ
る。
や電波など、２次元干渉測定が可能な任意の波動を用い
ることができる。マイクロ波帯の電波の波源像を求める
場合には、いったん、受信信号を中間周波数帯の信号に
周波数変換してからフーリエ変換を行うことも可能であ
る。
【００１４】本発明の方法において、相関値の算出に際
しては、例えば、コヒーレンス関数の値が所定のしきい
値を越えている周波数範囲内で積分演算を実行する。ま
た、波源像再生のための演算処理方法としては、例え
ば、波動干渉における近傍界領域に対応するもの、フレ
ネル領域に対応するもの、フラウンホーファー領域に対
応するものなどを用意する。さらに、本発明の方法にお
いては、再生された波源像の位相情報と振幅情報とを同
時に表示することができ、この場合には、異なるトリガ
条件で測定を実行して位相情報と振幅情報を比較するこ
とによって、時間的、空間的コヒーレンスを評価するこ
とが可能になる。
しては、例えば、コヒーレンス関数の値が所定のしきい
値を越えている周波数範囲内で積分演算を実行する。ま
た、波源像再生のための演算処理方法としては、例え
ば、波動干渉における近傍界領域に対応するもの、フレ
ネル領域に対応するもの、フラウンホーファー領域に対
応するものなどを用意する。さらに、本発明の方法にお
いては、再生された波源像の位相情報と振幅情報とを同
時に表示することができ、この場合には、異なるトリガ
条件で測定を実行して位相情報と振幅情報を比較するこ
とによって、時間的、空間的コヒーレンスを評価するこ
とが可能になる。
【００１５】本発明の波源可視化装置は、観測対象から
の放射波動に対して２次元干渉測定を実施し、観測対象
における波源像を可視化する装置であって、２次元干渉
測定の観測面に対する基準となる位置での放射波動を検
出して第１の信号として出力する第１のセンサと、観測
面内の任意の位置での放射波動を検出して第２の信号と
して出力する第２のセンサと、第１の信号及び第２の信
号をそれぞれフーリエ変換して第１のスペクトル及び第
２のスペクトルを算出するフーリエ変換手段と、第１の
スペクトルのパワースペクトルの平均値と、第２のスペ
クトルのパワースペクトルの平均値と、第１のスペクト
ルと第２のスペクトルの間のクロススペクトルの平均値
とを求める平均値計算手段と、クロススペクトルの平均
値と各パワースペクトルから算出された平均値とに基づ
きコヒーレンス関数を決定するコヒーレンス関数計算手
段と、観測面内の点ごとに、クロススペクトルを得てこ
れを第１のスペクトルのパワースペクトルの平方根で規
格化し、規格化されたクロススペクトルに対して観測周
波数帯域内でコヒーレンス関数の値に応じた積分演算を
実行して相関値を算出する相関値算出手段と、観測対象
と観測面との距離に応じて波源像の再生のための演算処
理方法を選択し、各点ごとの相関値に基づいて波源像を
再生する波源像再生処理手段と、を有する。
の放射波動に対して２次元干渉測定を実施し、観測対象
における波源像を可視化する装置であって、２次元干渉
測定の観測面に対する基準となる位置での放射波動を検
出して第１の信号として出力する第１のセンサと、観測
面内の任意の位置での放射波動を検出して第２の信号と
して出力する第２のセンサと、第１の信号及び第２の信
号をそれぞれフーリエ変換して第１のスペクトル及び第
２のスペクトルを算出するフーリエ変換手段と、第１の
スペクトルのパワースペクトルの平均値と、第２のスペ
クトルのパワースペクトルの平均値と、第１のスペクト
ルと第２のスペクトルの間のクロススペクトルの平均値
とを求める平均値計算手段と、クロススペクトルの平均
値と各パワースペクトルから算出された平均値とに基づ
きコヒーレンス関数を決定するコヒーレンス関数計算手
段と、観測面内の点ごとに、クロススペクトルを得てこ
れを第１のスペクトルのパワースペクトルの平方根で規
格化し、規格化されたクロススペクトルに対して観測周
波数帯域内でコヒーレンス関数の値に応じた積分演算を
実行して相関値を算出する相関値算出手段と、観測対象
と観測面との距離に応じて波源像の再生のための演算処
理方法を選択し、各点ごとの相関値に基づいて波源像を
再生する波源像再生処理手段と、を有する。
【００１６】本発明の装置においては、第１の信号を蓄
積してフーリエ変換手段に出力する第１のデータメモリ
と、第２の信号を蓄積してフーリエ変換手段に出力する
第２のデータメモリと、第１のデータメモリ及び第２の
データメモリに対してデータ蓄積のためのトリガ信号を
発生するトリガ回路とをさらに設けることができる。ま
た、波形像再生処理手段で再生された波源像について振
幅情報と位相情報とを同時に表示する表示手段をさらに
設けることができる。
積してフーリエ変換手段に出力する第１のデータメモリ
と、第２の信号を蓄積してフーリエ変換手段に出力する
第２のデータメモリと、第１のデータメモリ及び第２の
データメモリに対してデータ蓄積のためのトリガ信号を
発生するトリガ回路とをさらに設けることができる。ま
た、波形像再生処理手段で再生された波源像について振
幅情報と位相情報とを同時に表示する表示手段をさらに
設けることができる。
【００１７】
【発明の実施の形態】次に、本発明の実施の形態につい
て、図面を参照して説明する。図１は本発明の実施の一
形態の波源像可視化装置の構成を示すブロック図であ
り、この波源像可視化装置は、本発明の方法に基づいて
可視化波源像を表示するものである。
て、図面を参照して説明する。図１は本発明の実施の一
形態の波源像可視化装置の構成を示すブロック図であ
り、この波源像可視化装置は、本発明の方法に基づいて
可視化波源像を表示するものである。
【００１８】複数の波源集合を有する観測対象１１に対
して、距離ｚだけ離れた位置に走査観測面１２を設け、
走査観測面１２内を移動する走査センサ１３を配置す
る。また、走査センサ１３とは別に、移動しない固定セ
ンサ１４を配置する。走査センサ１３は第２のセンサに
対応し、固定センサ１４は第１のセンサに対応する。走
査センサ１３と固定センサ１４は、観測波動が音波であ
るときには例えばマイクロホンであり、電波であるとき
にはアンテナである。以下の説明では、走査観測面１２
での２次元座標をｘ,ｙで表わし、観測対象１１での２
次元座標をｘ',ｙ'で表わす。また本発明では、観測対
象１１と走査センサ１３や固定センサ１４との距離ｒ
が、可視化スペクトルの帯域幅Δω（観測周波数帯域の
角周波数表示による帯域幅）に比べて十分小さい、すな
わち、ｒ≪ｖ／Δω（ｖは波動の速度）であって、観測
波動を平均波長λ0で扱えることを前提とする。
して、距離ｚだけ離れた位置に走査観測面１２を設け、
走査観測面１２内を移動する走査センサ１３を配置す
る。また、走査センサ１３とは別に、移動しない固定セ
ンサ１４を配置する。走査センサ１３は第２のセンサに
対応し、固定センサ１４は第１のセンサに対応する。走
査センサ１３と固定センサ１４は、観測波動が音波であ
るときには例えばマイクロホンであり、電波であるとき
にはアンテナである。以下の説明では、走査観測面１２
での２次元座標をｘ,ｙで表わし、観測対象１１での２
次元座標をｘ',ｙ'で表わす。また本発明では、観測対
象１１と走査センサ１３や固定センサ１４との距離ｒ
が、可視化スペクトルの帯域幅Δω（観測周波数帯域の
角周波数表示による帯域幅）に比べて十分小さい、すな
わち、ｒ≪ｖ／Δω（ｖは波動の速度）であって、観測
波動を平均波長λ0で扱えることを前提とする。
【００１９】走査センサ１３及び固定センサ１４の出力
側には、それぞれ、走査センサ１３及び固定センサ１４
からの観測信号ｂ(ｔ),ａ(ｔ)のうちの所定の周波数帯
域の成分のみを帯域通過させて中間周波数ＩＦの信号に
変換する前変換部１５,１６が設けられている。各前変
換部１５,１６には、基準周波数ｆrefが供給されてい
る。前変換部１５,１６は、例えばＲＦスペクトラムア
ナライザをゼロスパンモードで基準周波数ｆrefに位相
ロックさせることで実現できる。音波や比較的低い周波
数の電波を観測波動とするような場合には、デジタル信
号処理によって帯域制限や周波数変換処理を行う前変換
部を使用してもよい。前変換部１５,１６の出力側に
は、前変換部１５,１６の出力をサンプリングして記憶
するデータメモリ１８,１９が配置している。外部同期
入力あるいは前変換部１５の出力に応じてトリガ信号を
発生するトリガ回路１７が設けられており、各データメ
モリ１８,１９はトリガ信号に応じてデータのサンプリ
ングを開始するように構成されている。また、各データ
メモリ１８,１９にはサンプリング周波数ｆSが供給され
ている。
側には、それぞれ、走査センサ１３及び固定センサ１４
からの観測信号ｂ(ｔ),ａ(ｔ)のうちの所定の周波数帯
域の成分のみを帯域通過させて中間周波数ＩＦの信号に
変換する前変換部１５,１６が設けられている。各前変
換部１５,１６には、基準周波数ｆrefが供給されてい
る。前変換部１５,１６は、例えばＲＦスペクトラムア
ナライザをゼロスパンモードで基準周波数ｆrefに位相
ロックさせることで実現できる。音波や比較的低い周波
数の電波を観測波動とするような場合には、デジタル信
号処理によって帯域制限や周波数変換処理を行う前変換
部を使用してもよい。前変換部１５,１６の出力側に
は、前変換部１５,１６の出力をサンプリングして記憶
するデータメモリ１８,１９が配置している。外部同期
入力あるいは前変換部１５の出力に応じてトリガ信号を
発生するトリガ回路１７が設けられており、各データメ
モリ１８,１９はトリガ信号に応じてデータのサンプリ
ングを開始するように構成されている。また、各データ
メモリ１８,１９にはサンプリング周波数ｆSが供給され
ている。
【００２０】フーリエ変換部２０は、サンプリングされ
てデータメモリ１８,１９に格納されたデータをそれぞ
れフーリエ変換してスペクトルＳa(ω),Ｓb(ω)を算出
するためのものであり、フーリエ変換部２０の出力側に
は、スペクトルＳa(ω),Ｓb(ω)それぞれのパワースペ
クトルの平均値とこれらスペクトルＳa(ω),Ｓb(ω)の
クロススペクトルの平均値を算出する平均値計算部２１
が設けられている。ωは角周波数である。さらに、平均
値計算部２１での算出結果に基づいてコヒーレンス関数
γ2(ω)を算出する第１の演算部２２と、コヒーレンス
関数γ2(ω)から定まる積分操作関数ｆ(ω)と平均値計
算部２１での算出結果に基づいて相関値Ｃab(ｘ,ｙ)を
算出する第２の演算部２３と、相関値Ｃabに基づいて波
形像を再生する波形像再生処理部２４と、再生された波
形像を表示する表示部２５とが設けられている。第１の
演算部２２及び第２の演算部２３は、それぞれ、コヒー
レンス関数計算手段、相関値算出手段に対応する。後述
するように、距離ｚと観測波動の波長λ0との関係によ
って波形像の再生方法が異なるため、波形像再生処理部
２４には距離ｚも入力している。上述のトリガ回路１
７、各データメモリ１８,１９及びフーリエ変換部２０
は、例えば、ＦＦＴスペクトラムアナライザ（デジタル
スペクトラムアナライザ）を用いて実現できる。
てデータメモリ１８,１９に格納されたデータをそれぞ
れフーリエ変換してスペクトルＳa(ω),Ｓb(ω)を算出
するためのものであり、フーリエ変換部２０の出力側に
は、スペクトルＳa(ω),Ｓb(ω)それぞれのパワースペ
クトルの平均値とこれらスペクトルＳa(ω),Ｓb(ω)の
クロススペクトルの平均値を算出する平均値計算部２１
が設けられている。ωは角周波数である。さらに、平均
値計算部２１での算出結果に基づいてコヒーレンス関数
γ2(ω)を算出する第１の演算部２２と、コヒーレンス
関数γ2(ω)から定まる積分操作関数ｆ(ω)と平均値計
算部２１での算出結果に基づいて相関値Ｃab(ｘ,ｙ)を
算出する第２の演算部２３と、相関値Ｃabに基づいて波
形像を再生する波形像再生処理部２４と、再生された波
形像を表示する表示部２５とが設けられている。第１の
演算部２２及び第２の演算部２３は、それぞれ、コヒー
レンス関数計算手段、相関値算出手段に対応する。後述
するように、距離ｚと観測波動の波長λ0との関係によ
って波形像の再生方法が異なるため、波形像再生処理部
２４には距離ｚも入力している。上述のトリガ回路１
７、各データメモリ１８,１９及びフーリエ変換部２０
は、例えば、ＦＦＴスペクトラムアナライザ（デジタル
スペクトラムアナライザ）を用いて実現できる。
【００２１】次に、この装置を用いた波源像の可視化に
ついて説明する。
ついて説明する。
【００２２】走査センサ１３を移動させることにより、
走査観測面１２内の任意の点(ｘ,ｙ)において観測対象
１１からの波動を受信して信号ｂ(ｔ)を取得し、同時
に、固定センサ１４によっても波動を受信して信号ａ
(ｔ)を取得する。そして、前変換部１５,１６により、
所定の周波数帯域に帯域制限してから、信号ｂ(ｔ),ａ
(ｔ)を中間周波数ＩＦの信号に周波数変換し、データメ
モリ１８,１９により、観測周波数帯域幅より十分に高
いサンプリングレートによって、Ｔ秒間にわたり、これ
らの信号をサンプリングする。そして、サンプリングさ
れたデータをフーリエ変換部２０によってフーリエ変換
することにより、信号ａ(ｔ)に対するスペクトルＳ
a(ω)と信号ｂ(ｔ)に対するスペクトルＳb(ω)とを得
る。スペクトルＳa(ω),Ｓb(ω)は平均値計算部２１に
入力し、各スペクトルＳa(ω),Ｓb(ω)のパワースペク
トルの平方根の平均値〈｜Ｓa(ω)｜〉,〈｜Ｓb(ω)
｜〉や、スペクトルＳa(ω),Ｓb(ω)間のクロススペク
トルのパワースペクトル〈Ｓa *(ω)Ｓb(ω)〉が算出さ
れる。ここで、〈・〉はアンサンブル平均を示し、*は複
素共役を示している。
走査観測面１２内の任意の点(ｘ,ｙ)において観測対象
１１からの波動を受信して信号ｂ(ｔ)を取得し、同時
に、固定センサ１４によっても波動を受信して信号ａ
(ｔ)を取得する。そして、前変換部１５,１６により、
所定の周波数帯域に帯域制限してから、信号ｂ(ｔ),ａ
(ｔ)を中間周波数ＩＦの信号に周波数変換し、データメ
モリ１８,１９により、観測周波数帯域幅より十分に高
いサンプリングレートによって、Ｔ秒間にわたり、これ
らの信号をサンプリングする。そして、サンプリングさ
れたデータをフーリエ変換部２０によってフーリエ変換
することにより、信号ａ(ｔ)に対するスペクトルＳ
a(ω)と信号ｂ(ｔ)に対するスペクトルＳb(ω)とを得
る。スペクトルＳa(ω),Ｓb(ω)は平均値計算部２１に
入力し、各スペクトルＳa(ω),Ｓb(ω)のパワースペク
トルの平方根の平均値〈｜Ｓa(ω)｜〉,〈｜Ｓb(ω)
｜〉や、スペクトルＳa(ω),Ｓb(ω)間のクロススペク
トルのパワースペクトル〈Ｓa *(ω)Ｓb(ω)〉が算出さ
れる。ここで、〈・〉はアンサンブル平均を示し、*は複
素共役を示している。
【００２３】第１の演算部２２では、平均値計算部２１
での計算結果を利用して、下記式に基づき、所定の観測
周波数帯域におけるコヒーレンス関数γ2(ω)が計算さ
れる。
での計算結果を利用して、下記式に基づき、所定の観測
周波数帯域におけるコヒーレンス関数γ2(ω)が計算さ
れる。
【００２４】
【数２】 そして、このコヒーレンス関数関数γ2(ω)に対してし
きい値αを設定することにより、積分操作関数ｆ(ω)が
得られる。
きい値αを設定することにより、積分操作関数ｆ(ω)が
得られる。
【００２５】
【数３】 第２の演算部２３では、下記式に基づき、走査観測面１
２上での走査センサ１３の位置(ｘ,ｙ)ごとに、｜Ｓ
a(ω)｜で規格化された複素相関値Ｃab(ｘ,ｙ)が計算さ
れる。ω0は観測周波数帯域の中心角周波数である。
２上での走査センサ１３の位置(ｘ,ｙ)ごとに、｜Ｓ
a(ω)｜で規格化された複素相関値Ｃab(ｘ,ｙ)が計算さ
れる。ω0は観測周波数帯域の中心角周波数である。
【００２６】
【数４】 この相関値Ｃab(ｘ,ｙ)は、遅延時間τが０であるとき
の信号Ｓa(ｔ)と信号Ｓb(ｔ)との相互相関に相当する値
となる。このとき、積分操作関数ｆ(ω)を得るために使
用するコヒーレンス関数γ2(ω)は、走査センサ１３を
移動する度に新しい値に更新してもよいし、走査観測面
１２上の代表する点(ｘ0,ｙ0)で得た値を共通値として
他の点で使用してもよい。また、相関値Ｃab(ｘ,ｙ)を
得るために使用されるクロススペクトルと固定センサ１
４の受信信号ａ(ｔ)のパワースペクトルの平方根とは、
平均値操作を行わずに用いてもよい。これに対し、コヒ
ーレンス関数γ2(ω)の演算には、必ず、平均化された
クロススペクトルと平均化されたパワースペクトル平方
根を使用しなければならない。
の信号Ｓa(ｔ)と信号Ｓb(ｔ)との相互相関に相当する値
となる。このとき、積分操作関数ｆ(ω)を得るために使
用するコヒーレンス関数γ2(ω)は、走査センサ１３を
移動する度に新しい値に更新してもよいし、走査観測面
１２上の代表する点(ｘ0,ｙ0)で得た値を共通値として
他の点で使用してもよい。また、相関値Ｃab(ｘ,ｙ)を
得るために使用されるクロススペクトルと固定センサ１
４の受信信号ａ(ｔ)のパワースペクトルの平方根とは、
平均値操作を行わずに用いてもよい。これに対し、コヒ
ーレンス関数γ2(ω)の演算には、必ず、平均化された
クロススペクトルと平均化されたパワースペクトル平方
根を使用しなければならない。
【００２７】走査観測面１２上の各点について相関値Ｃ
ab(ｘ,ｙ)を求めたら、波形像再生処理部２４におい
て、求めた相関値に基づき波形像を再生し、例えば鳥瞰
図表示によって表示部２５に表示する。波形像の再生
は、観測対象１１と走査観測面１２との距離ｚに応じ
て、近傍界領域での観測、フレネル(Fresnel)領域での
観測、フラウンホーファー(Fraunhofer)領域での観測の
３通りの方法で実行される。
ab(ｘ,ｙ)を求めたら、波形像再生処理部２４におい
て、求めた相関値に基づき波形像を再生し、例えば鳥瞰
図表示によって表示部２５に表示する。波形像の再生
は、観測対象１１と走査観測面１２との距離ｚに応じ
て、近傍界領域での観測、フレネル(Fresnel)領域での
観測、フラウンホーファー(Fraunhofer)領域での観測の
３通りの方法で実行される。
【００２８】観測波動の平均波長λ0に対してｚが著し
く小さい場合（ｚ≪λ0）に代表される近傍界領域の観
測では、相関値Ｃab(ｘ,ｙ)をそのままで表示するよう
にする。
く小さい場合（ｚ≪λ0）に代表される近傍界領域の観
測では、相関値Ｃab(ｘ,ｙ)をそのままで表示するよう
にする。
【００２９】近傍界領域とフラウンホーファー領域の中
間の領域であるフレネル領域の観測では、 Ｆ-1[Ｆ[Ｃab(ｘ,ｙ)]・Ｆ[Ｐ(ｘ,ｙ;ｚ)]] を表示するようにする。ここでＦ[・]はフーリエ変換、
Ｆ-1[・]はフーリエ逆変換を示し、Ｐ(ｘ,ｙ;ｚ)は下記
式で表わされる。
間の領域であるフレネル領域の観測では、 Ｆ-1[Ｆ[Ｃab(ｘ,ｙ)]・Ｆ[Ｐ(ｘ,ｙ;ｚ)]] を表示するようにする。ここでＦ[・]はフーリエ変換、
Ｆ-1[・]はフーリエ逆変換を示し、Ｐ(ｘ,ｙ;ｚ)は下記
式で表わされる。
【００３０】
【数５】 一方、ｚ≫λ0に代表されるフラウンホーファー領域の
観測では、
観測では、
【００３１】
【外２】 で表わされる値に応じて表示が行われるようにする。こ
こでｋ0＝２π／λ0であり、ξとηはそれぞれ走査観測
面１２から見た方位角と仰角であって、ξ＝ｘ'／ｚ、
η＝ｙ'／ｚである。
こでｋ0＝２π／λ0であり、ξとηはそれぞれ走査観測
面１２から見た方位角と仰角であって、ξ＝ｘ'／ｚ、
η＝ｙ'／ｚである。
【００３２】このように再生される振幅値は、波源の空
間コヒーレンシに関わりなく、ω0±(Δω／２)の範囲
での絶対平均振幅値を与える。実際には、観測対象１１
の大きさと観測波動の平均波長λ0との関係に応じて、
距離ｚを定めるようにすることが望ましい。すなわち、
平均波長λ0が観測面積に対して大きい場合には、ｚ＜
λ0として近傍界領域で観測が行われるようにする。一
方、観測対象１１の面積が走査観測面１２の面積よりも
大きい場合には、ｚ≫λ0としてフラウンホーファー領
域で観測が行われるようにし、これらの中間である場合
にはフレネル領域での観測が行われるようにする。
間コヒーレンシに関わりなく、ω0±(Δω／２)の範囲
での絶対平均振幅値を与える。実際には、観測対象１１
の大きさと観測波動の平均波長λ0との関係に応じて、
距離ｚを定めるようにすることが望ましい。すなわち、
平均波長λ0が観測面積に対して大きい場合には、ｚ＜
λ0として近傍界領域で観測が行われるようにする。一
方、観測対象１１の面積が走査観測面１２の面積よりも
大きい場合には、ｚ≫λ0としてフラウンホーファー領
域で観測が行われるようにし、これらの中間である場合
にはフレネル領域での観測が行われるようにする。
【００３３】表示部２５において各再生位置での振幅と
位相を同時に示すことで、空間・時間的な波源のコヒー
レンスを知ることができる。観測データを得るためのト
リガ条件を変えながら観測を行った場合、トリガ条件ご
との振幅の変化の様子は時間コヒーレンシを示し、位相
の変化の様子は空間コヒーレンシを示す。振幅と位相の
同時表示は、例えば鳥瞰図表示を行う場合には、山の高
さで振幅を表わし、色や明るさで位相を示すことによっ
て行える。あるいは、振幅を明るさ、位相を色として２
次元的に表現してもよいし、等高線図により振幅を等高
線で表わし位相を色または明るさで表現してもよい。
位相を同時に示すことで、空間・時間的な波源のコヒー
レンスを知ることができる。観測データを得るためのト
リガ条件を変えながら観測を行った場合、トリガ条件ご
との振幅の変化の様子は時間コヒーレンシを示し、位相
の変化の様子は空間コヒーレンシを示す。振幅と位相の
同時表示は、例えば鳥瞰図表示を行う場合には、山の高
さで振幅を表わし、色や明るさで位相を示すことによっ
て行える。あるいは、振幅を明るさ、位相を色として２
次元的に表現してもよいし、等高線図により振幅を等高
線で表わし位相を色または明るさで表現してもよい。
【００３４】ここで、本実施の形態における波源像の可
視化に関する制約条件について説明する。観測周波数帯
域幅Δωについては、観測対象１１と観測位置との距離
をｒとして、この距離ｒにおいて帯域幅Δω内での周波
数変動に対して位相がほぼ一定に見えるように設定する
必要がある。また、コヒーレンス関数γ2(ω)での評価
を行うために、走査センサ１３と固定センサ１４で同じ
信号成分を受信し、走査センサ１３側の遅延時間と固定
センサ１４側の遅延時間の差に比べてサンプリング時間
Ｔが十分大きく、かつ、このサンプリング時間Ｔによっ
て定まる周波数分解能によって、走査センサ１３からの
スペクトルＳa(ω)と固定センサ１４からのスペクトル
Ｓb(ω)とがスペクトル領域で識別できることが必要で
ある。
視化に関する制約条件について説明する。観測周波数帯
域幅Δωについては、観測対象１１と観測位置との距離
をｒとして、この距離ｒにおいて帯域幅Δω内での周波
数変動に対して位相がほぼ一定に見えるように設定する
必要がある。また、コヒーレンス関数γ2(ω)での評価
を行うために、走査センサ１３と固定センサ１４で同じ
信号成分を受信し、走査センサ１３側の遅延時間と固定
センサ１４側の遅延時間の差に比べてサンプリング時間
Ｔが十分大きく、かつ、このサンプリング時間Ｔによっ
て定まる周波数分解能によって、走査センサ１３からの
スペクトルＳa(ω)と固定センサ１４からのスペクトル
Ｓb(ω)とがスペクトル領域で識別できることが必要で
ある。
【００３５】以上の観点からすると、サンプリング周波
数ｆSをできるだけ低く設定するためには、前変換部１
５，１６での帯域制限の帯域幅と観測周波数帯域幅Δω
とを等しくすればよい。しかしながら、帯域制限にはア
ナログ回路によるバンドパスフィルタを使用するのが一
般的であり、前変換部１５と前変換部１６とでこのバン
ドパスフィルタの特性が異なる可能性がある。そこで、
バンドパスフィルタの帯域幅を観測周波数帯域幅Δωよ
りも大きくし、バンドパスフィルタの周波数特性が比較
的理想に近い部分で観測周波数帯域を設定するようにす
ればよい。
数ｆSをできるだけ低く設定するためには、前変換部１
５，１６での帯域制限の帯域幅と観測周波数帯域幅Δω
とを等しくすればよい。しかしながら、帯域制限にはア
ナログ回路によるバンドパスフィルタを使用するのが一
般的であり、前変換部１５と前変換部１６とでこのバン
ドパスフィルタの特性が異なる可能性がある。そこで、
バンドパスフィルタの帯域幅を観測周波数帯域幅Δωよ
りも大きくし、バンドパスフィルタの周波数特性が比較
的理想に近い部分で観測周波数帯域を設定するようにす
ればよい。
【００３６】以上、本発明の実施の形態について説明し
たが、本発明は走査観測面内を１個の走査センサを移動
させて測定を行うことに限定されるものではない。例え
ば、走査観測面内に２次元的に複数の固定したセンサを
配列し、これらセンサによって観測波動を同時に受信す
るようにしてもよい。この場合には、これら２次元的に
配列したセンサの１つを上述の固定センサとしてもよい
し、いくつかのセンサからの受信信号を加算して上述の
信号ａ(ｔ)（固定センサからの信号に相当）を得るよう
にしてもよい。
たが、本発明は走査観測面内を１個の走査センサを移動
させて測定を行うことに限定されるものではない。例え
ば、走査観測面内に２次元的に複数の固定したセンサを
配列し、これらセンサによって観測波動を同時に受信す
るようにしてもよい。この場合には、これら２次元的に
配列したセンサの１つを上述の固定センサとしてもよい
し、いくつかのセンサからの受信信号を加算して上述の
信号ａ(ｔ)（固定センサからの信号に相当）を得るよう
にしてもよい。
【００３７】
【発明の効果】以上説明したように本発明は、固定セン
サからの受信信号と走査センサからの受信信号のクロス
スペクトルやそれぞれのパワースペクトルを算出し、コ
ヒーレンスの度合を示すコヒーレンス関数を計算し、ク
ロススペクトルを得てこれを固定センサからの信号のパ
ワースペクトルの平方根で規格化し、規格化されたクロ
ススペクトルを所定の観測周波数帯域内でコヒーレンス
関数の値に応じて積分演算を実行して相関値を算出する
とともに、観測対象と走査観測面との距離に応じて波源
像の再生のための演算処理方法を選択し、各点ごとの相
関値に基づいて波源像を再生することにより、コヒーレ
ンシに関する事前の知識を必要とせず、部分コヒーレン
トな波源分布の評価が可能であって、スペクトルの速い
変化に追従できるようになるという効果がある。
サからの受信信号と走査センサからの受信信号のクロス
スペクトルやそれぞれのパワースペクトルを算出し、コ
ヒーレンスの度合を示すコヒーレンス関数を計算し、ク
ロススペクトルを得てこれを固定センサからの信号のパ
ワースペクトルの平方根で規格化し、規格化されたクロ
ススペクトルを所定の観測周波数帯域内でコヒーレンス
関数の値に応じて積分演算を実行して相関値を算出する
とともに、観測対象と走査観測面との距離に応じて波源
像の再生のための演算処理方法を選択し、各点ごとの相
関値に基づいて波源像を再生することにより、コヒーレ
ンシに関する事前の知識を必要とせず、部分コヒーレン
トな波源分布の評価が可能であって、スペクトルの速い
変化に追従できるようになるという効果がある。
【図１】本発明の実施の一形態の波源可視化装置の構成
を示すブロック図である。
を示すブロック図である。
【図２】従来の波源可視化装置の構成を示すブロック図
である。
である。
１１ 観測対象 １２ 走査観測面 １３ 走査センサ １４ 固定センサ １５,１６ 前変換部 １７ トリガ回路 １８,１９ データメモリ ２０ フーリエ変換部 ２１ 平均値計算部 ２２,２３ 演算部 ２４ 波形像再生処理部 ２５ 表示部
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 ＦＩ 技術表示箇所 Ｇ０１Ｓ 15/89 Ｇ０１Ｓ 15/89 Ｂ
Claims (7)
- 【請求項１】 観測対象からの放射波動に対して２次元
干渉測定を実施し、前記観測対象における波源像を可視
化する方法であって、 前記２次元干渉測定の観測面に対する基準となる位置で
の放射波動を検出して第１の信号として出力する第１の
センサと、前記観測面内の任意の位置での放射波動を検
出して第２の信号として出力する第２のセンサとを使用
し、 前記第１の信号及び前記第２の信号をそれぞれフーリエ
変換して第１のスペクトル及び第２のスペクトルを獲得
し、 前記第１のスペクトルのパワースペクトルと、前記第２
のスペクトルのパワースペクトルと、前記第１のスペク
トルと前記第２のスペクトルの間のクロススペクトルと
を求め、 前記クロススペクトルの平均値と前記各パワースペクト
ルから算出された平均値とに基づきコヒーレンス関数を
決定し、 前記観測面内の点ごとに、前記クロススペクトルを得て
これを前記第１のスペクトルのパワースペクトルの平方
根で規格化し、前記規格化されたクロススペクトルに対
して観測周波数帯域内で前記コヒーレンス関数の値に応
じた積分演算を実行して相関値を算出し、 前記観測対象と前記観測面との距離に応じて波源像の再
生のための演算処理方法を選択し、前記各点ごとの相関
値に基づいて波源像を再生する波源像可視化方法。 - 【請求項２】 前記相関値の算出に際し、前記コヒーレ
ンス関数の値が所定のしきい値を越えている周波数範囲
で前記積分演算を実行する請求項１に記載の波源像可視
化方法。 - 【請求項３】 再生された波源像の位相情報と振幅情報
とを同時に表示する請求項１または２に記載の波源像可
視化方法。 - 【請求項４】 異なるトリガ条件で測定を実行して位相
情報と振幅情報を比較し、時間的、空間的コヒーレンス
を評価する請求項３に記載の波源像可視化方法。 - 【請求項５】 観測対象からの放射波動に対して２次元
干渉測定を実施し、前記観測対象における波源像を可視
化する装置であって、 前記２次元干渉測定の観測面に対する基準となる位置で
の放射波動を検出して第１の信号として出力する第１の
センサと、 前記観測面内の任意の位置での放射波動を検出して第２
の信号として出力する第２のセンサと、 前記第１の信号及び前記第２の信号をそれぞれフーリエ
変換して第１のスペクトル及び第２のスペクトルを算出
するフーリエ変換手段と、 前記第１のスペクトルのパワースペクトルの平均値と、
前記第２のスペクトルのパワースペクトルの平均値と、
前記第１のスペクトルと前記第２のスペクトルの間のク
ロススペクトルの平均値とを求める平均値計算手段と、 前記クロススペクトルの平均値と前記各パワースペクト
ルから算出された平均値とに基づきコヒーレンス関数を
決定するコヒーレンス関数計算手段と、 前記観測面内の点ごとに、前記クロススペクトルを得て
これを前記第１のスペクトルのパワースペクトルの平方
根で規格化し、前記規格化されたクロススペクトルに対
して観測周波数帯域内で前記コヒーレンス関数の値に応
じた積分演算を実行して相関値を算出する相関値算出手
段と、 前記観測対象と前記観測面との距離に応じて波源像の再
生のための演算処理を選択し、前記各点ごとの相関値に
基づいて波源像を再生する波源像再生処理手段と、を有
する波源像可視化方法。 - 【請求項６】 前記第１の信号を蓄積して前記フーリエ
変換手段に出力する第１のデータメモリと、前記第２の
信号を蓄積して前記フーリエ変換手段に出力する第２の
データメモリと、前記第１のデータメモリ及び前記第２
のデータメモリに対してデータ蓄積のためのトリガ信号
を発生するトリガ回路と、を有する請求項５に記載の波
源像可視化装置。 - 【請求項７】 前記波形像再生処理手段で再生された波
源像について振幅情報と位相情報とを同時に表示する表
示手段をさらに有する請求項５または６に記載の波源像
可視化装置。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28984895A JP3642844B2 (ja) | 1995-11-08 | 1995-11-08 | 波源像可視化方法及び装置 |
| DE19645953A DE19645953B4 (de) | 1995-11-08 | 1996-11-07 | Verfahren zum Bestimmen einer Korrelation und Gerät zur Durchführung des Verfahrens |
| US08/746,293 US5748314A (en) | 1995-11-08 | 1996-11-07 | Correlation function measurement method and apparatus, and wave source image visualization method and apparatus based on correlation function measurement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28984895A JP3642844B2 (ja) | 1995-11-08 | 1995-11-08 | 波源像可視化方法及び装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH09134113A true JPH09134113A (ja) | 1997-05-20 |
| JP3642844B2 JP3642844B2 (ja) | 2005-04-27 |
Family
ID=17748556
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP28984895A Expired - Fee Related JP3642844B2 (ja) | 1995-11-08 | 1995-11-08 | 波源像可視化方法及び装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3642844B2 (ja) |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6188356B1 (en) | 1998-05-12 | 2001-02-13 | Advantest Corporation | Radio monitoring method and apparatus |
| US6198539B1 (en) | 1997-08-22 | 2001-03-06 | Advantest Corporation | Hologram observation method and hologram observation apparatus |
| US6594605B2 (en) | 2000-02-03 | 2003-07-15 | Advantest Corp. | Correlation function measuring method and apparatus |
| JP2005207847A (ja) * | 2004-01-22 | 2005-08-04 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波ホログラフィ法とｍｕｓｉｃ法を用い電波源の特定に用いられる電波監視装置 |
| JP2005227205A (ja) * | 2004-02-16 | 2005-08-25 | Toshiba Corp | 目標距離測定装置 |
| JP2006284356A (ja) * | 2005-03-31 | 2006-10-19 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波ホログラフィ電波源探査装置 |
| JP2006337280A (ja) * | 2005-06-03 | 2006-12-14 | Toshiba Corp | 電波発生源可視化装置及び電波発生源可視化方法 |
| JP2006337281A (ja) * | 2005-06-03 | 2006-12-14 | Toshiba Corp | 電波発生源可視化装置及び電波発生源可視化方法 |
| JP2006349631A (ja) * | 2005-06-20 | 2006-12-28 | Toshiba Corp | 電波監視装置 |
| JP2007017393A (ja) * | 2005-07-11 | 2007-01-25 | Toshiba Corp | 電波発射源表示装置 |
| JP2007212648A (ja) * | 2006-02-08 | 2007-08-23 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波発射源可視化装置及びその方法 |
| JP2007212228A (ja) * | 2006-02-08 | 2007-08-23 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波発射源可視化装置及びその方法 |
| JP2007212229A (ja) * | 2006-02-08 | 2007-08-23 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波発射源可視化装置及びその方法 |
| JP2007212227A (ja) * | 2006-02-08 | 2007-08-23 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波発射源可視化装置 |
| JP2018525617A (ja) * | 2015-07-09 | 2018-09-06 | コグニティヴ システムズ コーポレイション | 高周波カメラシステム |
| WO2021117269A1 (ja) * | 2019-12-12 | 2021-06-17 | 国立研究開発法人宇宙航空研究開発機構 | フェーズドアレイアンテナ装置 |
| JP2021518569A (ja) * | 2018-03-19 | 2021-08-02 | セブン ベル ゲーエムベーハーＳｅｖｅｎ Ｂｅｌ ＧｍｂＨ | 音源の空間的位置の特定のための装置、システムおよび方法 |
| WO2021176872A1 (ja) * | 2020-03-03 | 2021-09-10 | ソニーグループ株式会社 | 情報表示制御方法と情報表示制御装置およびプログラム |
| JP2022134585A (ja) * | 2021-03-03 | 2022-09-15 | 株式会社デンソーＥｍｃエンジニアリングサービス | 信号源推定装置 |
-
1995
- 1995-11-08 JP JP28984895A patent/JP3642844B2/ja not_active Expired - Fee Related
Cited By (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6198539B1 (en) | 1997-08-22 | 2001-03-06 | Advantest Corporation | Hologram observation method and hologram observation apparatus |
| US6347060B1 (en) | 1997-08-22 | 2002-02-12 | Advantest Corp. | Hologram observation method and hologram observation apparatus and multi-dimensional hologram data processing method and multi-dimensional hologram data processing apparatus |
| US6188356B1 (en) | 1998-05-12 | 2001-02-13 | Advantest Corporation | Radio monitoring method and apparatus |
| US6594605B2 (en) | 2000-02-03 | 2003-07-15 | Advantest Corp. | Correlation function measuring method and apparatus |
| DE10103822B4 (de) * | 2000-02-03 | 2009-06-04 | Advantest Corp. | Korrelationsfunktions-Messverfahren und -Vorrichtung |
| JP2005207847A (ja) * | 2004-01-22 | 2005-08-04 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波ホログラフィ法とｍｕｓｉｃ法を用い電波源の特定に用いられる電波監視装置 |
| JP4544869B2 (ja) * | 2004-01-22 | 2010-09-15 | 総務大臣 | 電波ホログラフィ法とｍｕｓｉｃ法を用い電波源の特定に用いられる電波監視装置 |
| JP2005227205A (ja) * | 2004-02-16 | 2005-08-25 | Toshiba Corp | 目標距離測定装置 |
| JP4519478B2 (ja) * | 2004-02-16 | 2010-08-04 | 株式会社東芝 | 目標距離測定装置 |
| JP2006284356A (ja) * | 2005-03-31 | 2006-10-19 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波ホログラフィ電波源探査装置 |
| JP4723916B2 (ja) * | 2005-06-03 | 2011-07-13 | 株式会社東芝 | 電波発生源可視化装置及び電波発生源可視化方法 |
| JP2006337281A (ja) * | 2005-06-03 | 2006-12-14 | Toshiba Corp | 電波発生源可視化装置及び電波発生源可視化方法 |
| JP2006337280A (ja) * | 2005-06-03 | 2006-12-14 | Toshiba Corp | 電波発生源可視化装置及び電波発生源可視化方法 |
| JP4686265B2 (ja) * | 2005-06-20 | 2011-05-25 | 株式会社東芝 | 電波監視装置 |
| JP2006349631A (ja) * | 2005-06-20 | 2006-12-28 | Toshiba Corp | 電波監視装置 |
| JP2007017393A (ja) * | 2005-07-11 | 2007-01-25 | Toshiba Corp | 電波発射源表示装置 |
| JP2007212227A (ja) * | 2006-02-08 | 2007-08-23 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波発射源可視化装置 |
| JP2007212229A (ja) * | 2006-02-08 | 2007-08-23 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波発射源可視化装置及びその方法 |
| JP2007212228A (ja) * | 2006-02-08 | 2007-08-23 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波発射源可視化装置及びその方法 |
| JP2007212648A (ja) * | 2006-02-08 | 2007-08-23 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波発射源可視化装置及びその方法 |
| JP2018525617A (ja) * | 2015-07-09 | 2018-09-06 | コグニティヴ システムズ コーポレイション | 高周波カメラシステム |
| JP2021518569A (ja) * | 2018-03-19 | 2021-08-02 | セブン ベル ゲーエムベーハーＳｅｖｅｎ Ｂｅｌ ＧｍｂＨ | 音源の空間的位置の特定のための装置、システムおよび方法 |
| WO2021117269A1 (ja) * | 2019-12-12 | 2021-06-17 | 国立研究開発法人宇宙航空研究開発機構 | フェーズドアレイアンテナ装置 |
| WO2021176872A1 (ja) * | 2020-03-03 | 2021-09-10 | ソニーグループ株式会社 | 情報表示制御方法と情報表示制御装置およびプログラム |
| JP2022134585A (ja) * | 2021-03-03 | 2022-09-15 | 株式会社デンソーＥｍｃエンジニアリングサービス | 信号源推定装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3642844B2 (ja) | 2005-04-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3642844B2 (ja) | 波源像可視化方法及び装置 | |
| US5748314A (en) | Correlation function measurement method and apparatus, and wave source image visualization method and apparatus based on correlation function measurement | |
| US10386339B2 (en) | Modal vibration analysis system | |
| US9225365B2 (en) | Sampling method, apparatus, probe, reception beamforming apparatus, and medical imaging system performing the sampling method | |
| JPH0782053B2 (ja) | アパーチャ合成ラジオメータ | |
| Osumi et al. | Microwave holographic imaging method with improved resolution | |
| Liu et al. | About the application of the van Cittert-Zernike theorem in ultrasonic imaging | |
| US20160135694A1 (en) | Medical radar method and system | |
| JP6261839B1 (ja) | 合成開口レーダ信号処理装置 | |
| JP4388113B2 (ja) | 三次元形状測定装置 | |
| JP4734057B2 (ja) | 移動物体検出装置 | |
| US6140960A (en) | Hologram observation method for three-dimensional wave source distribution, and stereoscopic directivity estimation method of antenna and wave distribution observation method based on hologram observation | |
| WO2017141352A1 (ja) | 物体検出装置及びセンサ装置 | |
| CN113671495B (zh) | 基于Zynq平台的太赫兹雷达探测系统以及方法 | |
| CN114894793B (zh) | 基于消除伪影的成像方法、成像系统及服务器 | |
| EP1198757B1 (en) | Correlation analysis in the phase domain | |
| Ozsoy et al. | Pencil back-projection method for SAR imaging | |
| JPH10239364A (ja) | 波動分布観測方法及び波動分布表示方法 | |
| JP6257072B2 (ja) | 白色干渉計装置による表面形状の測定方法 | |
| JPH10232252A (ja) | ３次元波源分布のホログラム観測方法 | |
| JP2823045B2 (ja) | 目標運動解析方法 | |
| KR101641129B1 (ko) | x 밴드 레이더를 이용한 해류속도 및 연안 수심 측정 시스템 및 그 시스템을 이용한 해류속도 및 연안 수심 측정방법 | |
| JPS61751A (ja) | 超音波媒体の特性測定装置 | |
| JP3446148B2 (ja) | 超音波ドプラ診断装置 | |
| Comesana et al. | Assessing vehicle exterior noise using a virtual phased array (VPA) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040401 |
|
| 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: 20050112 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050126 |
|
| R150 | Certificate of patent (=grant) or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080204 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090204 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100204 Year of fee payment: 5 |
|
| LAPS | Cancellation because of no payment of annual fees |