JPH0694763A - 電波源可視化装置 - Google Patents

電波源可視化装置

Info

Publication number
JPH0694763A
JPH0694763A JP24676392A JP24676392A JPH0694763A JP H0694763 A JPH0694763 A JP H0694763A JP 24676392 A JP24676392 A JP 24676392A JP 24676392 A JP24676392 A JP 24676392A JP H0694763 A JPH0694763 A JP H0694763A
Authority
JP
Japan
Prior art keywords
radio wave
antenna
output
radio
image
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
Application number
JP24676392A
Other languages
English (en)
Other versions
JP3286970B2 (ja
Inventor
Hitoshi Kitayoshi
均 北吉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Advantest Corp
Original Assignee
Advantest Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Advantest Corp filed Critical Advantest Corp
Priority to JP24676392A priority Critical patent/JP3286970B2/ja
Publication of JPH0694763A publication Critical patent/JPH0694763A/ja
Application granted granted Critical
Publication of JP3286970B2 publication Critical patent/JP3286970B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

(57)【要約】 (修正有) 【目的】 電波源と観測面の距離zや観測電波波長が電
波源間隔1よりも大きい場合でもたやすく電波源像を分
離することができることを可能とする。 【構成】 観測面13のx,y座標の各点(x,y)に
ついて水平偏波受信アンテナ22、垂直偏波受信アンテ
ナ23を移動させながら受信出力を求め、これら出力を
それぞれ離散的フーリエ変換器38,39で離散的フー
リエ変換し、また観測面13の近くに設けた固定アンテ
ナ21の出力を離散的フーリエ変換し、その固定アンテ
ナ出力のフーリエ変換に対する水平アンテナ出力のフー
リエ変換、垂直アンテナ出力のフーリエ変換の各離散的
フーリエの比を求める。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】この発明は、例えばEMC(エレ
クトロ マグネティック コンパティビリティ)の測定
や、電波センシング、即ち反射、屈折、吸収などのパタ
ーン解析に利用され、電波源の位置の可視化を行う電波
源可視化装置に関する。
【0002】
【従来の技術】図5に従来の電波源可視化装置を示す。
密着させて測定することができない電波発生源を含む被
測定体11に対し、これより距離zだけ離れ、その被測
定体11と平行に、サーチコイル12により線状走査、
あるいは二次元走査を手動または自動的に行い、そのサ
ーチコイル12の出力をデジタル信号にAD変換器15
により変換してパーソナルコンピュータよりなる制御表
示器14に取り込んで受信電波強度の分布図を作ってい
た。
【0003】
【発明が解決しようとする課題】つまり従来においては
電波強度計やサーチコイルとスペクトラムアナライザの
組み合わせにより前記観測面を走査して受信電波強度分
布を得ている。しかし電波源から電界強度計に入射する
電波は直接波のみならず天井や床などからの反射波があ
り特に反射波の存在を知ることができない。
【0004】また複数の電波源が存在する場合、その複
数の電波間で干渉が生じるため電波源の間隔lが電波源
の電波の波長λより大きくても、被測定体11、つまり
電波源と観測面13との間隔zが電波波長λより大きい
場合は電波源の像を得ることができない。さらに複数の
電波波面が合成されるため電波源の間隔lが電波波長λ
より小さい場合や距離zが電波波長より小さい場合にお
いても電波源の間隔lが距離zより小さい場合は電波源
の像の分離ができない。
【0005】また電波強度計を走査中に電波源より発生
する電波強度が変動すると、得られる電波像が歪んでし
まうなどの欠点があった。
【0006】
【課題を解決するための手段】請求項1の発明によれ
ば、観測面近くに固定アンテナが設けられ、また観測面
を垂直偏波受信アンテナと水平偏波受信アンテナとを組
みにして移動させられ、その垂直偏波受信アンテナ及び
水平偏波受信アンテナの固定アンテナに対する各出力比
が取られる。これらの比について各受信点についての空
間伝搬関数との畳込み積分を行って受動ホログラム像の
再生演算を行う。これにより観測像から電波源像つまり
電波源分布が得られる。
【0007】請求項2の発明によれば、垂直偏波受信ア
ンテナ及び水平偏波受信アンテナを組みとして観測面が
走査され、その両アンテナの各測定受信位置における出
力振幅の比のアークタンジェントを求めて合成電波入射
角を演算し、この合成電波入射角のラプラシアン、つま
り空間二階微分をおこない、これを像として電波源分布
を得る。
【0008】さらに請求項1または2の発明において各
測定点における垂直偏波受信アンテナ、水平偏波受信ア
ンテナの各出力の位相角を演算し、その位相角の差の状
況分布を求め、これにより直接波と反射波とのそんざい
割合の分布を求める。また請求項1または2の発明にお
いて、さらに受動ホログラム像を周波数を変えて二つ作
り、その二つの像の同一点(位置)における位相差を垂
直偏波成分、水平偏波成分について求め、これにより反
射波の直接波に対する遅れ時間が求まる。
【0009】
【実施例】図1に請求項1の発明の実施例を示す。観測
面13の近くに固定アンテナ21が設けられる。この観
測面13の各点を移動走査することができるように水平
偏波受信アンテナ22、垂直偏波受信アンテナ23の組
が設けられる。つまりアンテナ22,23は移動アンテ
ナであり、この移動アンテナ22,23により観測面1
3を面走させて、その各点における受信電界強度を測定
する。これらアンテナ21,22,23の受信出力はそ
れぞれ帯域通過フィルタ24,25,26に通されて所
望帯域のみ取り出され、さらにそれぞれ増幅器27,2
8,29で増幅され、その増幅出力は局部発振器30か
らの信号により混合器31,32,33でそれぞれ周波
数変換され、それらの信号出力は低域通過ろ波器34,
35,36に通され、そのろ波出力は離散的フーリエ変
換器37,38,39によりそれぞれ離散的フーリエ変
換され、そのフーリエ変換された結果SR ,SH ,SV
がそれぞれ制御演算表示部41へ供給される。
【0010】制御演算表示部41は移動アンテナ22,
23をそれぞれ同時にその姿勢を保持した状態で観測面
13の各点を順次移動させて観測面13を走査し、その
各点においてそれぞれ得られた信号、つまり観測面13
の座標位置(x,y)における離散的フーリエ変換の結
果SR (x,y),SH (x,y),SV (x,y)に
ついて固定アンテナ出力フーリエ変換SR (x,y)に
対する水平偏波受信アンテナ22の出力フーリエ変換S
H (x,y)の比hH (x,y)と垂直偏波受信アンテ
ナ23の出力フーリエ変換SV (x,y)の比h
V (x,y)を演算する。
【0011】 hV (x,y)=SV (x,y)/SR (x,y) hH (x,y)=SH (x,y)/SR (x,y) …(1) この受信出力比の測定面13における状態は受動ホログ
ラム、つまり固定アンテナの出力に対する空間的な位相
分布、言い替えれば固定アンテナ出力と移動アンテナ出
力とを干渉させた相関スペクトルを表す。
【0012】また各点におけるこの固定アンテナ出力に
対する比hH (x,y)の位相角とhV (x,y)の位
相角の差θH-V (x,y)をそれぞれ各測定点(x,
y)について演算する。 θH-V (x,y)=∠hH (x,y)−∠hV (x,y) …(2) この位相角度の差θH-V (x,y)は観測面13におけ
る直接波と反射波の存在の割合の分布を示すことにな
る。つまり反射波の有無を知ることができ、またどこの
部分が反射波が多いかを知ることができる。
【0013】さらに各点(x,y)の固定アンテナ出力
に対する移動アンテナ出力の振幅の比のアークタンジェ
ントを求めることによって合成電波の入射角θ
H/V (x,y)の状態を得る。 θH/V (x,y)=tan-1(‖hH (x,y)‖/‖hV (x,y)‖) …(3) また各点における固定アンテナに対する出力比の強度M
(x,y)を求める。
【0014】 M(x,y)=√(‖hH (x,y)‖2 +‖hV (x,y)‖2 ) …(4) これにより各点(x,y)における電波強度が分かり、
観測面13における電波強度分布を知ることができる。
さらに電波源像を求めるため各点の合成電波入射角θ
H/V (x,y)のラプラシアン、つまり空間二階微分Δ
θH/V (x,y)を求める。この場合(4)式で得た電
波強度分布でその各ラプラシアン値に重み付けを行い、
受信電界強度が小さいところでの電波源像に対するθ
H/V (x,y)のラプラシアン値を小さくすることによ
り雑音の影響を小さくすることが望ましい。
【0015】さらに電波源間隔lが電波波長λよりも大
きい場合や電波源と観測面との間の距離zが電波波長λ
より大きい場合、電波干渉像となるため固定アンテナ出
力に対する移動アンテナ出力の比する比hV (x,y)
及びhH (x,y)に対してホログラム再生演算を行い
V (x,y)、IH (x,y)を得る。 I(X′,Y′)=∬〔h(x,y)×{exp(−j2πR/λ)}/R〕 dxdy R=√{z2 +(x−x′)2 +(y−y′)2 } このようにホログラム再生演算を行った結果のI
v (x,y)とIH (x,y)を(2)式、(3)式、
(4)式におけるhV (x,y)、hH (x,y)の代
わりに用いてそれぞれ位相角差や合成入射角及び強度な
どを求める。
【0016】また上述において、測定されたhV (x,
y)、hH (x,y)や演算されたΔθH/V (x,y)
に雑音が多く含まれるような場合には、それぞれに対し
て二次元空間低域通過フィルタを演算処理により作用さ
せるとよい。つまりf(x,y)に対する二次元空間低
域通過フィルタ作用は次式で表される。 F(x,y)=∬f{(x,y)×sin(λ-1×R)λ-1・R}dx ,dy R =(√{(x−x′)2 +(y−y′)2 })×2π λは光速を空間遮断周波数で割った値 ここでΔθH/V (x,y)の各点における値を求め、そ
の測定観測面上における状態を見ると、電波源像がz>
lでも得られることについて述べると、観測面で一点に
入射する電波の方向及び強さは複数の電波源から放射さ
れた電波のベクトル和である。従って電波源の間隔lに
対して観測面と電波源との距離zが大きくなればなるほ
ど合成電波の波面が滑らかに変化し、電波源の分離が難
しくなる。しかし観測面に対して電波源の方向が垂直に
なる位置、つまり電波源より観測面に対して垂線を降ろ
した位置ではこの両者間の距離zに関係なく合成電波の
入射角度が変曲し、即ちその角度の見る方向が逆とな
る。従って合成電波の入射角のアンテナ移動方向におけ
る二階微分ΔθH/V (x,y)のピークは観測面から垂
直に見た電波源の位置を与えることになる。
【0017】この関係を図2及び図3に電波源間距離を
3.3メートルとし、移動アンテナ22,23を観測面
13に沿って10メートル移動させた場合における受信
アンテナの出力の合成振幅を曲線42として合成電波の
入射角を曲線43でそのラプラシアンを曲線44で示
す。図2Aは電波源と観測面との間の距離zが1メート
ル、図2Bはzが3メートル、図3Aはzが6メート
ル、図3Bはzが12メートルの場合である。これら曲
線からzが電波源間隔lより小さければ図2Aに示すよ
うに強度分布42、入射角分布43さらにそのラプラシ
アンの分布44においてもそれぞれ電波源位置において
大きなピークとなり、その電波源像を容易に求めること
ができる。しかし図2Bに示すようにzとlとがほぼ近
くなると強度分布42による電波源像を得ることは困難
であるが入射角分布43及びそのラプラシアン分布44
はまだ十分電波源像を得ることができる。図3Aに示す
ようにzが電波源間隔より倍程度大きくなると入射角の
分布43によっても像を知ることは困難であるが、その
ラプラシアン分布44においては明瞭なピークで電波源
を知ることができる。さらに図3Bに示すように電波源
間隔に対して距離zが相当大きくなっても入射角のラプ
ラシアン分布44のそのピーク値は小さくなるが、その
ピークは明らかに検出され、電波源の位置を知ることが
できる。
【0018】上述したように、ホログラムを再生しその
強度分布あるいは入射角のラプラシアン分布を求めても
これらが直接波により得られているのか、反射波により
得られているのかを分離することは困難である。これが
直接波か反射波によりその電波源が得られているかを知
るには(2)式で求めた位相角差の分布状態を求めれば
これにより区別することができる。
【0019】発振器30の発振周波数f0 を変えてそれ
ぞれについて受動ホログラム像、つまり式(1)のhV
(x,y)、hH (x,y)の分布状態をその二つの再
生像の位相差を求めることによって反射波の直接波に対
する遅れ時間を求めることができる。即ち周波数f1
より求めたIV 1 (x,y)と周波数f2 により求め
たIV 2 (x,y)との位相差を求め、また同様にI
H 1 (x,y)とI H 2 (x,y)の位相差を求め
る。このようにすることによって水平偏波成分における
反射と垂直偏波成分における反射の仕方が異なり、直接
波に対する反射波の遅れを知ることができる。
【0020】上述において、離散的フーリエ変換器3
7,38,39の代わりに高速フーリエ変換器を用いて
もよい。さらに移動アンテナとしては垂直偏波成分、水
平偏波成分の2方向成分ではなく3方向成分をそれぞれ
各別に検出する3つのアンテナを用いてもよい。また一
つの電波源像を得る為に電波源に対する観測面の角度を
変えて2回測定して、電波像を立体的に得ることもでき
る。移動アンテナ22,23の出力を切り替えて一つの
スペクトル解析部に供給するようにしてもよい。
【0021】請求項2の発明においては、固定アンテナ
21を省略して移動アンテナ22,23の出力のフーリ
エ変換SV (x,y),SH (x,y)をそれぞれ
(2)式、(3)式、(4)式におけるhV (x,
y),hH (x,y)の代わりに用いて上述したことを
実行するものである。次に測定データを示す。電波源間
隔l=10cmにおいて、測定周波数を200MHzと
し、観測面積を20cm×18cmにおいて、(4)式
で得た強度分布と入射角のラプラシアンとの状況を、z
を5、10、15、20、25、30、35センチにそ
れぞれ変えたばあいについてその大きさの値の分布状態
を観測面状に対応して図4に示した。この場合において
も合成振幅から分離不可能であった二つの電波源が入射
角のラプラシアンによって二つ斜線領域に明確に分離で
きるが、特にzが1より大においても斜線領域が二つあ
り、分離できることが分かる。また合成振幅と入射角の
ラプラシアンはともに距離zが大きくなるにつれて左側
にパターンが移動しているが、これは観測平面と電波面
平面の平衡ずれのためと考えられる。また入射角ラプラ
シアン像においても床、天井からの反射波イメージが確
認できることが分かる。
【0022】
【発明の効果】以上述べた様に、この発明によれば電波
源間隔1が電波波長λよりも大きく、また観測面の距離
zが電波波長λより大きい場合やこれらが波長よりも小
さくても電波源間隔1よりも観測面間距離zが大きい場
合のいずれの場合においても入射角のラプラシアンを求
めることによって電波像を得ることができる。さらに各
垂直偏波受信アンテナと水平偏波受信アンテナとの出力
位相角の差の分布を求めることによって直接波と反射波
とを区別することができる。また電波源出力が変動して
も請求項1の発明では固定アンテナ出力で移動アンテナ
出力を基準化しているため影響を受けず電波像が歪む恐
れはない。
【図面の簡単な説明】
【図1】この発明の実施例を示すブロック図。
【図2】強度分布、入射角分布、そのラプラシアン分布
を示す図。
【図3】強度分布、入射角分布、そのラプラシアン分布
を示す図。
【図4】電波源から観測面までの距離を各種の状態とし
た場合に置ける強度の二次元像及びラプラシアンの二次
元像を示す図。
【図5】従来の電波像観測装置を示す図。

Claims (4)

    【特許請求の範囲】
  1. 【請求項1】 観測平面近くにおかれた固定アンテナ
    と、 上記観測平面を移動可能な垂直偏波受信アンテナ及び水
    平偏波受信アンテナと、 上記固定アンテナの出力に対する上記垂直偏波受信アン
    テナの出力比及び上記水平偏波受信アンテナの出力比を
    取る手段と、 これら両比に対し受動ホログラム像の再生演算を行う手
    段と、 を具備する電波源可視化装置。
  2. 【請求項2】 観測平面を移動可能に設けられた垂直偏
    波受信アンテナ及び水平偏波受信アンテナと、 上記観測平面内における各同一位置において、上記両ア
    ンテナよりの出力の振幅比のアークタンジェントを計算
    して合成電波の入射角を演算する手段と、 その入射角のラプラシアンを演算して上記観測面上にお
    けるラプラシアン像を得る手段と、 を具備する電波源可視化装置。
  3. 【請求項3】 同一地点における上記垂直偏波受信アン
    テナ出力の位相角と上記水平偏波受信アンテナの出力の
    位相角との分布を上記観測平面について演算する手段を
    具備することを特徴とする請求項1または請求項2記載
    の電波源可視化装置。
  4. 【請求項4】 複素受動ホログラム再生像を周波数を変
    えて二つ作り、その二つの像の位相差を求めその位相差
    像を得る手段を具備することを特徴とする請求項1また
    は2記載の電波源可視化装置。
JP24676392A 1992-09-16 1992-09-16 電波源可視化装置 Expired - Fee Related JP3286970B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP24676392A JP3286970B2 (ja) 1992-09-16 1992-09-16 電波源可視化装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP24676392A JP3286970B2 (ja) 1992-09-16 1992-09-16 電波源可視化装置

Publications (2)

Publication Number Publication Date
JPH0694763A true JPH0694763A (ja) 1994-04-08
JP3286970B2 JP3286970B2 (ja) 2002-05-27

Family

ID=17153304

Family Applications (1)

Application Number Title Priority Date Filing Date
JP24676392A Expired - Fee Related JP3286970B2 (ja) 1992-09-16 1992-09-16 電波源可視化装置

Country Status (1)

Country Link
JP (1) JP3286970B2 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0829468A (ja) * 1994-07-20 1996-02-02 Nec Gumma Ltd 電磁界強度測定装置
US6147482A (en) * 1997-06-26 2000-11-14 Fuji Xerox Co., Ltd. Method for and apparatus of detecting and displaying radiated noise
JP2018525617A (ja) * 2015-07-09 2018-09-06 コグニティヴ システムズ コーポレイション 高周波カメラシステム

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0829468A (ja) * 1994-07-20 1996-02-02 Nec Gumma Ltd 電磁界強度測定装置
US6147482A (en) * 1997-06-26 2000-11-14 Fuji Xerox Co., Ltd. Method for and apparatus of detecting and displaying radiated noise
JP2018525617A (ja) * 2015-07-09 2018-09-06 コグニティヴ システムズ コーポレイション 高周波カメラシステム

Also Published As

Publication number Publication date
JP3286970B2 (ja) 2002-05-27

Similar Documents

Publication Publication Date Title
JP3144688B2 (ja) 空気現象の検出および測定のための方法および装置ならびにそのような装置に使用する送信機および受信機
Amin et al. Wideband synthetic aperture beamforming for through-the-wall imaging [lecture notes]
Yuan et al. 3D Fourier imaging based on 2D heterodyne detection at THz frequencies
US8060339B2 (en) Multistatic concealed object detection
Gao et al. Implementation of the phase shift migration in MIMO-sidelooking imaging at terahertz band
Gennarelli et al. Radar imaging through cinderblock walls: Achievable performance by a model-corrected linear inverse scattering approach
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
Hayashi et al. $ F $–$ k $ Filter Designs to Suppress Direct Waves for Bistatic Ground Penetrating Radar
JP3286970B2 (ja) 電波源可視化装置
Lay-Ekuakille et al. Optimizing and post processing of a smart beamformer for obstacle retrieval
JP3621988B2 (ja) イオノゾンデ装置
CN105629229B (zh) 飞机无损检测系统和方法
JPH11326480A (ja) 特定域電波可視化方法及び装置
Popescu et al. Point spread function estimation for a terahertz imaging system
JP2013210268A (ja) 2次元像再構成方法
Nilsson et al. Radar interferometric measurements with a planar patch antenna array
WO2020061294A1 (en) Crack detection using synthetic aperture radar microwave reflectometer
CN117555038A (zh) 毫米波成像设备、安检设备以及安检方法
JPH10232252A (ja) 3次元波源分布のホログラム観測方法
JP3034675B2 (ja) マイクロ波による内部温度分布の計測装置
RU2541682C1 (ru) Способ анализа акустического поля сложного источника
CN107315180A (zh) 一种两频点太赫兹全息成像方法及系统
CN222636301U (zh) 毫米波成像设备、安检设备
RU2478982C2 (ru) Способ определения звукового давления движущегося протяженного источника акустического поля
Jandt et al. Radar imaging based on multiple non-coherent antennas

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20020129

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080315

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090315

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees