JPH0312710B2 - - Google Patents

Info

Publication number
JPH0312710B2
JPH0312710B2 JP58064547A JP6454783A JPH0312710B2 JP H0312710 B2 JPH0312710 B2 JP H0312710B2 JP 58064547 A JP58064547 A JP 58064547A JP 6454783 A JP6454783 A JP 6454783A JP H0312710 B2 JPH0312710 B2 JP H0312710B2
Authority
JP
Japan
Prior art keywords
output
spatial frequency
cosine
maximum point
beams
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.)
Expired - Lifetime
Application number
JP58064547A
Other languages
English (en)
Other versions
JPS59190676A (ja
Inventor
Hiroshi Nishimura
Masao Igarashi
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.)
Oki Electric Industry Co Ltd
Original Assignee
Oki Electric Industry Co Ltd
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 Oki Electric Industry Co Ltd filed Critical Oki Electric Industry Co Ltd
Priority to JP58064547A priority Critical patent/JPS59190676A/ja
Publication of JPS59190676A publication Critical patent/JPS59190676A/ja
Publication of JPH0312710B2 publication Critical patent/JPH0312710B2/ja
Granted legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Description

【発明の詳細な説明】 (技術分野) 本発明は、同一の指向性パターンを持つ受波素
子の配列からなるアレイを用いてマルチビームを
形成し、該マルチビームの出力から信号源位置の
方向余弦を推定するにあたつて、前記受波素子の
指向性パターンによつて生ずる方向余弦の推定誤
差を除去する方向余弦推定方式に関するものであ
る。
(背景技術) ソーナーや音響測位及びレーダーにおいては、
空間上に直線(1次元)、平面(2次元)あるい
は3次元的に素子を配列したアレイを用いてマル
チビームを形成し、このマルチビームの離散出力
の最大点を求め、該最大点近傍のマルチビームの
出力に対し空間周波数領域において補間操作を施
し、補間後の出力の最大点を求め、該最大点に対
応する空間周波数により信号源位置の方向余弦を
求める方向余弦推定方式が広く用いられている。
ところで、従来の方式では、前記素子の指向性
は無指向性、すなわち方向に対する素子の受波感
度は一定であるという仮定を前提としている。と
ころが、前記素子は通常ある程度の指向性を持つ
ており、従つて従来の方式によつて方向余弦を推
定すると、素子の指向性のために誤差を生ずると
いう欠点があつた。
第1図は、アレイが直線(1次元)アレイの場
合の幾何学的説明図であり、111,112,…,
11i,…,11Mは各々受波素子、12は該受波
素子が配列される基準軸、13は信号源方向を示
す直線、θxは前記基準軸12に対する前記信号源
方向を示す直線13の方向余弦角である。
第2図は、各受波素子の指向性パターンが無指
向性であると仮定した従来の方向余弦推定方式の
機能ブロツク図であり、211,212,…,21
i,…,21Mは各々増幅器、221,222,…,
22i,…,22Mは各々帯域通過フイルタ、23
はビームフオーマ、y1,y2,…,yo,yNはマルチ
ビームの出力、24は第1の最大点検出器、kn
はマルチビームの出力が最大となる空間周波数の
離散値、25は補間器、y(k)は補間後のビー
ム出力、26は第2の最大点検出器、k^は前記補
間後のビーム出力y(k)が最大となる空間周波
数、27は変換器、cosθ^xは前記k^に対応した方向
余弦推定値、28は出力端子である。
信号源が直線(1次元)アレイの開口に比べ充
分遠方にある場合を考える。信号源方向13の方
向余弦角がθxであり、基準軸12の基準位置座標
での受波信号Sp(t)が Sp(t)=A・ej(2〓ft+〓) (1) で与えられるとする。ただしAは振幅、fは信号
周波数、φは初期位相である。すると第i番目素
子11iへの入力信号Si(t)は、θx方向から到来
する信号波面の基準点と第i番目素子11iとの
到来時間差を考えると で与えられる。ただしAiは振幅、ζiは第i番目1
1iの位置座標、cは伝搬速度である。ここで空
間周波数(spatial frequency)k〔文献:S.
Haykin,Editor J.H.Justice N.L.Owsley J.L.
Yen A.C.Kak: “Array Signal
Processing”,PRENTICE−HALL SIGNAL
PROCESSING SERIES,P119,1985〕を定義
する。
kΔ= 2πcos θx/λ (3) ただしλは信号の波長である。
空間周波数kの物理的な意味は、基準軸上の位
置座標ζiにおける、信号源方向の方向余弦角θx方
向に計測した時の基準位置座標との位相差がkζi
(radians)になることを示す。前記空間周波数k
とf/c=1/λの関係を用い、キヤリア成分及
び初期位相成分を省略し(2)式を書き直すと Si(t)=Ai・e-jk〓i (4) となる。
ビームフオーマ23は、空間周波数kのN個の
離散値k1,k2,…,ko,…,kNに対応する方向
余弦角θx1,θx2,…,θxo,θxN方向にN個のマル
チビームを形成し、該マルチビームの出力y1,
y2,…,yo,…,yNを出力する。ビームフオーマ
出力yoを得るため、前記ビームフオーマ23は、
第i番目素子11iからの信号に対し位相の補正
ejkn 〓i,ko=2πcos θxo/λを行う。前記位相補正
後 のすべてのM個の信号の和をとり、その積分時間
T2−T1でのパワーを次式yoとして出力する。
yo=|∫T2 T1 M/Σ/i=1 Si(t)・ejkn 〓i
dt
|2 (5) 第1の最大点検出器24はy1,y2,…,yo,
…,yNの最大値ynを求め、そのときの空間周波
数knを出力する。
補間器25は前記空間周波数knの近傍のマル
チビーム出力に補間操作を施し、空間周波数kの
連続関数としてビーム出力y(k)を出力する。
第2の最大点検出器26はy(k)の最大値
max y(k)を求め、そのときの空間周波数k^を
出力する。
変換器27は式(3)に基づいて前記k^から方向余
弦 cosθ^x=λk^/2π (6) を出力する。第3図は、素子111,112,11
i,…,11Mが無指向性、すなわち素子の規格化
した指向性パターンが全ての空間周波数kに対し
てb(k)=1が成立するときの、第2図で示される
従来の方式の説明図である。b(k)=1である場合
には、従来方式で求められる前期k¨は真値と一致
し、従つて前記cosθ^xも真値と一致する。
第4図は素子111,112,…,11i,…,1
1Mが無指向性でない場合に、従来方式をそのま
ま適用した場合の説明図である。
文献「R.J.Urick:“Principles of
Underwater Sound”,McGraw−Hill,P.57,
1967」で明らかにされているように、アレイを構
成する各素子111,112,…,11i,…,11
Mが同一の指向性パターンb(k)を持つときのマルチ
ビームの指向性パターンは、前記各素子が無指向
性であるときのマルチビームの指向性パターンと
前記素子の指向性パターンの積で与えられるか
ら、この場合の前記ビームフオーマの出力y′1,
y′2,…,y′o,…,y′Nは各素子が無指向性である
ときのビームフオーマの出力y1,y2,…,yo,
…,yNを用いて、 y′o=b(ko)・yo;n=1,2,…,n,…,
N (7) で表現できるる。ここでb(ko)は空間周波数の
離散値koにおける指向性パターン値、すなわち素
子感度である。
従つてy′1,y′2,…,y′o,…,y′Nに補間操作を
施し、その最大点から求めた空間周波数k^′は一般
に前記k^と一致しない。すなわち、従来の方向余
弦推定方式では、アレイを構成する素子が指向性
をもつ場合、推定誤差Δcosθx=λΔk/2π(ただし、 Δk=k^′−k^)を生ずることになる。
(発明の課題) 本発明の目的はこの欠点を除去するため、ビー
ムフオーマの出力に素子感度の逆数を掛け、等価
的に素子が無指向性である場合と同じビーム出力
を得、無指向性でない素子を使い従来方式を適用
することによつて生ずる空間周波数の推定誤差を
なくしたものであり、以下詳細に説明する。
(発明の構成および作用) 第5図は、本発明の実施例であつて511,5
12,…,51o,…,51Nは各々乗算器、52
1,522,…,52o,…,52Nは各々レジスタ
である。レジスタ521,522,…,52o,…,
52Nには素子の指向性パターンの空間周波数k1,
k2,…,ko,…,kNにおける値の逆数、すなわ
ち素子感度の逆数1/b(k1),1/b(k2),…,
1/b(ko),…,1/b(kN)の値が定数として
記憶されており、乗算器511,512,…,51
o,…,51Nは該定数と前記ビームフオーマ23
の出力y′1,y′2,…,y′o,…,y′Nとの積y′1/b
(k1),y′2/b(k2),…,y′o/b(ko),…,y
′N/
b(kN)を算出し、前記第1の最大点検出器24
と補間器25に出力する。
前記の式(7)より、乗算器511,512,…,5
1o,…,51Nの出力は、 y′o/b(ko)=b(ko) ・yo/b(ko)=yo ; n=1,2,…,N
(8) で与えられ、従つて第1の最大点検出器24によ
る最大点knの検出及び補間器25による補間操
作は、素子の指向性パターンとは無関係に、素子
が無指向性である場合と同じビーム出力y1,y2,
…,yo,…,yNに対してなされることになる。
以上説明したように本発明の実施例では、同一
の指向性をもつ素子を使用する方向余弦推定方式
においてはビームフオーマ出力に受波感度の逆数
を掛けることにより、(8)式で示されるように第1
の最大点検出器による最大点knの検出及び補間
器による補間操作は、素子の指向性パターンとは
無関係に、素子が無指向性である場合と同じビー
ム出力に対してなされることになり、従来方式の
適用によつて生ずる空間周波数の推定誤差Δkを
無くすことができ、従つて方向余弦の推定誤差
Δcosθxを無くすことができるという利点がある。
(発明の効果) 本発明は、アレイを構成する素子に指向性があ
る場合にも、マルチビームの出力を用いて信号源
位置の方向余弦を正確に推定することができると
いう利点があり、ソーナー、音響測位、レーダー
における方向余弦推定方式に利用することができ
る。
【図面の簡単な説明】
第1図は、アレイが直線(1次元)アレイの場
合の幾何学的説明図、第2図は、各受波素子の指
向性が無指向性であると仮定した場合の従来の方
向余弦推定方式の機能ブロツク図、第3図は各素
子が無指向性の場合の従来方式の説明図、第4図
は各素子が無指向性でなく、同一の指向性パター
ンをもつ場合に従来方式をそのまま適用した場合
の説明図、第5図は本発明の実施例の機能ブロツ
ク図である。 符号の説明(第2図、第5図)、211,212,
…,21i,…,21Mは各々増幅器、221,2
22,…,22i,…,22Mは各々帯域通過フイ
ルタ、23はビームフオーマ、y1,y2,…,yo,
…,yNはマルチビームの出力、24は第1の最大
点検出器、knはマルチビームの出力が最大とな
る空間周波数の離散値、25は補間器、y(k)は補
間後のビーム出力、26は第2の最大点検出器、
k^は前記補間後のビーム出力y(k)が最大となる空間
周波数、27は変換器、cosθ^xは前記k^に対応した
方向余弦推定値、28は出力端子である。511,
512,…,51o,…,51Nは乗算器、521,
522,…,52o,…,52Nはレジスタである。

Claims (1)

    【特許請求の範囲】
  1. 1 同一の指向性パターンを持つ受波素子の配列
    からなるアレイを用いて、空間周波数の複数の離
    散値に対応する方向にビーム(以下マルチビーム
    という)を形成し、該マルチビームの出力から信
    号源位置の方向余弦を推定する方向余弦推定装置
    において、前記空間周波数の離散値に対応する方
    向の素子感度の逆数と該方向のマルチビーム出力
    との積を算出する乗算器と、該乗算器の出力に対
    して補間を行なう補間器と、該補間器の出力が最
    大となる空間周波数を求める最大点検出器を有
    し、該空間周波数に対応する方向余弦を信号源位
    置の方向余弦の推定値とすることを特徴とする方
    向余弦推定方式。
JP58064547A 1983-04-14 1983-04-14 方向余弦推定方式 Granted JPS59190676A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP58064547A JPS59190676A (ja) 1983-04-14 1983-04-14 方向余弦推定方式

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP58064547A JPS59190676A (ja) 1983-04-14 1983-04-14 方向余弦推定方式

Publications (2)

Publication Number Publication Date
JPS59190676A JPS59190676A (ja) 1984-10-29
JPH0312710B2 true JPH0312710B2 (ja) 1991-02-20

Family

ID=13261352

Family Applications (1)

Application Number Title Priority Date Filing Date
JP58064547A Granted JPS59190676A (ja) 1983-04-14 1983-04-14 方向余弦推定方式

Country Status (1)

Country Link
JP (1) JPS59190676A (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0367303U (ja) * 1989-11-02 1991-07-01
US5227803A (en) * 1992-07-22 1993-07-13 Hughes Aircraft Company Transponder location and tracking system and method

Also Published As

Publication number Publication date
JPS59190676A (ja) 1984-10-29

Similar Documents

Publication Publication Date Title
AU2016227509B2 (en) Adaptive beamformer for sonar imaging
US4244026A (en) Velocity measuring correlation sonar
US6482160B1 (en) High resolution 3D ultrasound imaging system deploying a multidimensional array of sensors and method for multidimensional beamforming sensor signals
US9479867B2 (en) Method and circuitry for direction of arrival estimation using microphone array with a sharp null
EP2063292B1 (en) Calibrating a multibeam sonar apparatus
CN109283536A (zh) 一种多波束测深声呐水体成像波束形成算法
JP2002511591A (ja) センサアレイを用いての電波源の方向決定
JP7302293B2 (ja) 方位推定装置、方位推定方法、及びプログラム
RU2515179C1 (ru) Способ определения направления на гидроакустический маяк-ответчик в условиях многолучевого распространения навигационного сигнала
CN105022050A (zh) 一种多传感器阵列的水声信道离散噪声源抑制方法
CN115494485B (zh) 一种二维低旁瓣恒定分辨率前视声呐成像方法及系统
CN118818507B (zh) 一种宽幅角扫描成像装置、方法、系统
US7792313B2 (en) High precision beamsteerer based on fixed beamforming approach beampatterns
Ma et al. Spatiotemporal two-dimensional deconvolution beam imaging technology
US8824699B2 (en) Method of, and apparatus for, planar audio tracking
Blomberg et al. Adaptive beamforming applied to a cylindrical sonar array using an interpolated array transformation
JPH0511912B2 (ja)
Greening et al. Adaptive beamforming for nonstationary arrays
CN116593965A (zh) 一种基于远场机会声源的水平阵阵型估计方法
Xu et al. Coherent source direction estimation for three-row bathymetric sidescan sonars
JP3583908B2 (ja) 目標計測装置
JP2690606B2 (ja) 音源数決定方法
JP4712305B2 (ja) 電波方向探知装置
JP2001074835A (ja) バイスタティックソーナーの左右判別方法
JP2886330B2 (ja) 信号方位測定装置