JP3208424B2 - Surface array multi-channel adaptive phasing method and apparatus - Google Patents
Surface array multi-channel adaptive phasing method and apparatusInfo
- Publication number
- JP3208424B2 JP3208424B2 JP03267999A JP3267999A JP3208424B2 JP 3208424 B2 JP3208424 B2 JP 3208424B2 JP 03267999 A JP03267999 A JP 03267999A JP 3267999 A JP3267999 A JP 3267999A JP 3208424 B2 JP3208424 B2 JP 3208424B2
- Authority
- JP
- Japan
- Prior art keywords
- sensor
- sensors
- surface array
- channel adaptive
- array multi
- 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
Links
- 238000000034 method Methods 0.000 title claims description 43
- 230000003044 adaptive effect Effects 0.000 title claims description 37
- 238000004364 calculation method Methods 0.000 claims description 46
- 230000006978 adaptation Effects 0.000 claims description 32
- 230000035945 sensitivity Effects 0.000 claims description 13
- 230000003111 delayed effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 6
- 230000005855 radiation Effects 0.000 description 6
- 238000005070 sampling Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Description
【0001】[0001]
【発明の属する技術分野】本発明は、水中又は空中にお
いて、音波又は振動を発信し目標からの反響音又は振動
を受信するか或いは目標からの放射音又は振動を受信す
る音波又は振動のセンサー群を有するソーナー等の装置
に適用してサイドビーム方向から到来する雑音を低減可
能な面アレイマルチチャンネル適合整相方法及び装置に
関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a group of sensors for transmitting sound waves or vibrations and receiving reverberation sounds or vibrations from a target or receiving radiation sounds or vibrations from a target in water or in the air. The present invention relates to a surface array multi-channel adaptive phasing method and apparatus capable of reducing noise arriving from a side beam direction by being applied to a device such as a sonar having a phase difference.
【0002】[0002]
【従来の技術】水中又は空中において、音波を発信し目
標からの反響音を受信する受波器出力或いは目標からの
放射音を受信する受波器出力は周囲に存在する雑音に覆
われている。このため、目標からの受信信号を見つける
ために、受波器出力を整相処理して雑音レベルを低減さ
せる必要がある。2. Description of the Related Art Underwater or in the air, the output of a receiver for transmitting sound waves and receiving reverberation sound from a target or the output of a receiver for receiving sound radiated from a target is covered with noise existing in the surroundings. . Therefore, in order to find a received signal from a target, it is necessary to perform a phasing process on the output of the receiver to reduce the noise level.
【0003】水中又は空中において、音波を発信し目標
からの反響音を受信する受波器出力或いは目標からの放
射音を受信する受波器出力の従来の整相方法及び装置
は、各受波器出力にシェーディング係数と呼ばれる固定
の数値を掛けそれらを加算していた。A conventional phasing method and apparatus for a receiver output for transmitting a sound wave and receiving a reverberation sound from a target or a receiver output for receiving a radiation sound from a target in water or in the air is described below. The output of the container was multiplied by a fixed value called a shading coefficient and added.
【0004】一方、LMS(Least Mean Square)アル
ゴリズムは以下の計算により出力を計算する。On the other hand, the LMS (Least Mean Square) algorithm calculates the output by the following calculation.
【0005】入力ベクトルをXk、タップ重みをWkと
すると、出力ykは yk=Wk TXk (但し、k:サンプリング時期) となる。瞬時誤差信号εkは希望応答dkと出力ykの
差を取り εk=dk−yk=dk−Wk TXk で表す。タップ重みは Wk+1=Wk+2μεkXk で計算する。μはステップパラメータである。Assuming that the input vector is X k and the tap weight is W k , the output y k is y k = W k T X k (where k: sampling time). Instantaneous error signal epsilon k expressed by taking the difference between the output y k and the desired response d k ε k = d k -y k = d k -W k T X k. The tap weight is calculated by W k + 1 = W k +2 με k X k . μ is a step parameter.
【0006】しかし、LMS(Least Mean Square)ア
ルゴリズムは、信号出力と雑音出力との間に相関がある
場合には信号出力も削減されるため、現在のところ我が
国のソーナーの整相処理には用いられていない。However, since the LMS (Least Mean Square) algorithm reduces the signal output when there is a correlation between the signal output and the noise output, the LMS algorithm is currently used for phasing of sonar in Japan. Not been.
【0007】また、フロストの適合処理{O. L. Frost,
“An Algorithm for LinearyConstraind Adaptive Arra
y Processing”, Proc. IEEE, vol 60, No8, August,92
6-935(1972)}においては、Xkをk番目のサンプル、
Wkをk番目のサンプルに対するタップ重みとし、k+
1番目のタップ重みを次式で求めて順次出力y kを計算
する。 yk=WkXk Wk+1/2=Wk+2μykXk [e1 k+1/2…eL k+1/2]=[f−[1
1 … 1]Wk+1/2]/D[0007] Frost adaptation processing {O. L. Frost,
“An Algorithm for LinearyConstraind Adaptive Arra
y Processing ”, Proc. IEEE, vol 60, No8, August, 92
In 6-935 (1972) X, XkIs the k-th sample,
WkIs the tap weight for the k-th sample, and k +
The first tap weight is obtained by the following equation and sequentially output y kCalculate
I do. yk= WkXk Wk + 1/2= Wk+ 2μykXk [E1 k + 1/2… EL k + 1/2] = [F− [1
1 ... 1] Wk + 1/2] / D
【数1】 Wk+1=Wk+2μykXk+Ek+1/2 yk+1=Wk+1Xk+1 ここで、μはステップパラメータ、Dは素子数、Lはタ
ップ数、fはフィルターである。(Equation 1) W k + 1 = W k +2 μy k X k + E k + 1/2 y k + 1 = W k + 1 X k + 1 where μ is a step parameter, D is the number of elements, L is the number of taps, and f is a filter.
【0008】[0008]
【発明が解決しようとする課題】ところで、シェーディ
ング係数の選択の方法には、メインビーム幅を狭くする
方法、サイドビームの感度を下げる方法等種々の方法が
あり、周囲の雑音を考慮して係数を選択している。しか
し、これらの方法はいずれも雑音源は十分遠距離に存在
することが前提になっていた。このため、近距離にある
雑音源、及び存在方位や位置が時間とともに変動する雑
音源に対しては最適なシェーディング係数になっていな
いため、最適な整相処理になっていなかった。There are various methods for selecting a shading coefficient, such as a method for narrowing the width of the main beam and a method for decreasing the sensitivity of the side beam. Is selected. However, these methods all assume that the noise source is located at a sufficiently long distance. For this reason, the optimal shading coefficient has not been obtained for a noise source located at a short distance and a noise source whose existence azimuth and position fluctuates with time, so that the optimal phasing process has not been achieved.
【0009】本発明は、上記の点に鑑み、中央の受波器
と両端の2つの受波器からなる組を複数個配置すること
で、各組における中央の受波器と両端の受波器との指向
性の差及び位相差を利用し、これに適合処理計算を適用
して、雑音源がサイドビーム方向であれば、いずれの距
離に存在しても、いずれの方向に存在しても、常に雑音
レベルが最小になるような最適なシェーディング係数の
選択が可能な面アレイマルチチャンネル適合整相方法及
び装置を提供することを目的とする。In view of the above, the present invention provides a plurality of sets each including a central receiver and two receivers at both ends, so that the center receiver and the receivers at both ends in each set are arranged. Using the directivity difference and phase difference with the detector, and applying the adaptive processing calculation to this, if the noise source is in the side beam direction, it exists in any direction regardless of the distance. Another object of the present invention is to provide a surface array multi-channel adaptive phasing method and apparatus capable of always selecting an optimal shading coefficient so as to minimize the noise level.
【0010】本発明のその他の目的や新規な特徴は後述
の実施の形態において明らかにする。[0010] Other objects and novel features of the present invention will be clarified in embodiments described later.
【0011】[0011]
【課題を解決するための手段】上記目的を達成するため
に、本願請求項1の発明は、音波又は振動を発信し目標
からの反響音又は振動を受信するか或いは目標からの放
射音又は振動を受信する音波又は振動のセンサー出力の
整相方法において、直線上に等間隔に配置された3つの
センサーからなるセンサーの組を複数配置し、複数配置
された前記センサーの組がすべて中央のセンサーを共有
しており、各組における両端の2つのセンサーの加算出
力と中央のセンサー出力の受信指向性の差並びに位相の
差に適合処理計算を適用し、サイドローブ方向から到来
する雑音を低減させることを特徴とする。In order to achieve the above object, the invention according to claim 1 of the present invention transmits a sound wave or a vibration and receives a reverberation sound or a vibration from a target, or emits a sound or a vibration from a target. In the phasing method of the sound wave or vibration sensor output that receives a plurality of sets of three sensors arranged at equal intervals on a straight line, a plurality of sets of sensors are arranged.
Pairs of shared sensors share a central sensor
In addition, the adaptive processing calculation is applied to the difference between the reception directivity and the phase difference between the sum output of the two sensors at both ends and the center sensor output in each set to reduce the noise coming from the sidelobe direction. Features.
【0012】[0012]
【0013】本願請求項2の発明は、前記請求項1にお
いて、複数の前記センサーの組がM(1≦M)本の直線
上に配置されているものである。[0013] The present application claims 2 invention, our in claim 1
And a plurality of sets of the sensors are arranged on M (1 ≦ M) straight lines.
【0014】本願請求項3の発明は、前記請求項2にお
いて、前記M本の直線が面状配置であり、各直線同士の
なす角を等しい角度としたものである。According to a third aspect of the present invention, in the second aspect , the M straight lines have a planar arrangement, and angles formed by the straight lines are equal to each other.
【0015】本願請求項4の発明は、前記請求項2にお
いて、前記M本の直線が面状配置であり、各直線同士の
なす角の一部又は全部を異なる角度としたものである。According to a fourth aspect of the present invention, in the second aspect , the M straight lines have a planar arrangement, and some or all of the angles formed by the straight lines have different angles.
【0016】本願請求項5の発明は、前記請求項2にお
いて、i(1≦i≦M)本目の直線上に、前記センサー
の組をNi(1≦Ni)個配置したものである。According to a fifth aspect of the present invention, in the second aspect , N i (1 ≦ N i ) sets of the sensors are arranged on the i- th (1 ≦ i ≦ M) straight line. .
【0017】本願請求項6の発明は、前記請求項5にお
いて、N1=N2=…=Ni=…=NMとしたものであ
る。[0017] the invention of claim 6 is in the claim 5, in which the N 1 = N 2 = ... = N i = ... = N M.
【0018】本願請求項7の発明は、前記請求項5にお
いて、一部又は全部のNiを異なる数にしたものであ
る。According to a seventh aspect of the present invention, in the fifth aspect , some or all of the Nis are different.
【0019】本願請求項8の発明は、前記請求項5にお
いて、i本目の直線上の1組目の両端のセンサー間隔を
di1、2組目の両端のセンサー間隔をdi2、j組目
の両端のセンサー間隔をdij、Ni組目の両端のセン
サー間隔をdiNiとしたとき、di2=2×di1、
di3=3×di1、…、dij=j×di1、…、d
iNi=Ni×di1のように1組目の両端のセンサー
間隔の整数倍としたものである。According to an eighth aspect of the present invention, in the fifth aspect , the sensor interval between both ends of the first set on the i-th straight line is di1 , the sensor interval between both ends of the second set is di2 , and the jth set. sensor interval d ij of both ends of, when the sensor distance across the N i th set was d iNi, d i2 = 2 × d i1,
d i3 = 3 × d i1 ,..., d ij = j × d i1,.
It is an integral multiple of the sensor interval between both ends of the first set, such as iNi = N i × d i1 .
【0020】本願請求項9の発明は、前記請求項5にお
いて、i本目の直線上のj組目の両端のセンサー間隔を
dij(1≦j≦Ni)としたとき、前記センサー間隔
dijの一部又は全部を1組目の両端のセンサー間隔の
整数倍とは異なる間隔にしたものである。According to a ninth aspect of the present invention, in the fifth aspect , when the sensor interval between the two ends of the j-th set on the i-th straight line is d ij (1 ≦ j ≦ N i ), the sensor interval d A part or all of ij is different from an integral multiple of the sensor interval at both ends of the first set.
【0021】本願請求項10の発明は、前記請求項5に
おいて、1本目の直線の両端のセンサー間隔をd11、
2本目の直線の両端のセンサー間隔をd21、i本目の
直線の両端のセンサー間隔をdi1、M本目の直線の両
端のセンサー間隔をdM1としたとき、d11=d21
=…=di1=…=dM1としたものである。According to a tenth aspect of the present invention, in the fifth aspect , the distance between the sensors at both ends of the first straight line is d 11 ,
If the sensor interval at both ends of the second straight line is d 21 , the sensor interval at both ends of the i-th straight line is d i1 , and the sensor interval at both ends of the M-th straight line is d M1 , d 11 = d 21
=... = D i1 =... = D M1 .
【0022】本願請求項11の発明は、前記請求項5に
おいて、i本目の直線の両端のセンサー間隔をd
i1(1≦i≦M)としたとき、前記センサー間隔d
i1の一部又は全部を異なる間隔にしたものである。According to an eleventh aspect of the present invention, in the fifth aspect , the distance between the sensors at both ends of the i-th straight line is set to d.
i1 (1 ≦ i ≦ M), the sensor interval d
A part or all of i1 is provided at different intervals.
【0023】本願請求項12の発明は、前記請求項1か
ら11のいずれかにおいて、前記センサーの組における
両端の2つのセンサー出力の和の0.5倍の値と中央の
センサー出力との和が最小になるように、適合処理計算
を行うものである。The twelfth aspect of the present invention is directed to the first aspect.
11. In any one of 11., the adaptive processing calculation is performed so that the sum of the value of 0.5 times the sum of the outputs of the two sensors at both ends in the pair of sensors and the output of the central sensor is minimized. .
【0024】本願請求項13の発明は、前記請求項12
において、前記適合処理計算を、LMSアルゴリズムを
用いて行うものである。The invention of claim 13 of the present application is the invention of claim 12
In the above, the calculation of the adaptation processing is performed using an LMS algorithm.
【0025】本願請求項14の発明は、前記請求項13
において、前記適合処理計算を、フロストの適合処理を
用いて行うものである。The invention of claim 14 of the present application is the invention of claim 13
In the above, the calculation of the adaptation processing is performed by using the frost adaptation processing.
【0026】本願請求項15の発明は、前記請求項13
において、前記適合処理計算を、LMSアルゴリズムを
用いて行う場合に、前記センサーの組における両端の2
つのセンサーの加算出力に対するタップ重みと中央のセ
ンサー出力に対するタップ重みの和を常に1に保つこと
によって、目標からの受信信号に対する利得を常に1に
保つものである。The invention of claim 15 of the present application is the invention of claim 13
In the above, when the adaptation processing calculation is performed using the LMS algorithm, two ends at both ends in the sensor set are used.
By keeping the sum of the tap weights for the added outputs of the two sensors and the tap weights for the central sensor output at one, the gain for the received signal from the target is always kept at one.
【0027】本願請求項16の発明は、前記請求項15
において、前記センサーの組における両端の2つのセン
サーの加算出力に対するタップ重みと中央のセンサー出
力に対するタップ重みの和を常に1に保つために、適合
処理計算毎にタップ重みの和を求め、それぞれのタップ
重みをこの和で割り戻してタップ重みの合計を常に1に
保つものである。[0027] The invention of claim 16 of the present application is the invention of claim 15.
In order to always keep the sum of the tap weight for the added output of the two sensors at both ends in the set of sensors and the tap weight for the central sensor output at 1, the sum of the tap weights is calculated for each adaptation process calculation, The tap weights are divided by this sum to keep the total tap weights at one.
【0028】本願請求項17の発明は、前記請求項1
3,15又は16において、前記LMSアルゴリズムの
タップ重みの計算における前記タップ重みとして正負い
ずれの数値も選択し、前記タップ重みが負になった場合
に当該タップ重みの絶対値の上限値を定めるものであ
る。[0028] The invention of claim 17 of the present application is the invention of claim 1.
In 3, 15 or 16 , the tap weights in the calculation of tap weights of the LMS algorithm are both positive and negative, and when the tap weights become negative, the upper limit value of the absolute value of the tap weights is determined. It is.
【0029】本願請求項18の発明は、前記請求項1
3,15又は16において、前記LMSアルゴリズムの
タップ重みの計算における前記タップ重みが負になった
ら、その絶対値を取り、正の数値のみを選択するもので
ある。[0029] The invention of claim 18 of the present application is the above-described claim 1.
In 3, 15, or 16 , when the tap weight in the calculation of the tap weight of the LMS algorithm becomes negative, its absolute value is taken and only a positive numerical value is selected.
【0030】本願請求項19の発明は、音波又は振動を
発信し目標からの反響音又は振動を受信するか或いは目
標からの放射音又は振動を受信する音波又は振動のセン
サー出力の整相装置において、直線上に等間隔に配置さ
れた3つのセンサーからなるセンサーの組が複数配置さ
れかつ複数の前記センサーの組がすべて中央のセンサー
を共有しているとともに、各センサー出力をそれぞれ遅
延させる遅延器と、各センサーの組における両端の2つ
のセンサーの遅延出力をそれぞれ加算する加算器と、各
加算器出力に0.5をそれぞれ乗じる乗算器と、各セン
サーの組における中央の受波器の遅延出力と各乗算器出
力とを受けて適合処理計算を行う適合処理器とを備え、
該適合処理器の適合処理計算によりサイドローブの受信
感度を下げて雑音を低減させることを特徴とするもので
ある。According to a nineteenth aspect of the present invention, there is provided a phasing device for transmitting a sound wave or vibration and receiving a reverberation sound or vibration from a target or a sound wave or vibration sensor output for receiving a radiation sound or vibration from a target. , Placed at equal intervals on a straight line
Multiple sensor sets consisting of three sensors
And the set of sensors is a central sensor
And a delay unit that delays the output of each sensor, an adder that adds the delay outputs of the two sensors at both ends of each sensor set, and multiplies each adder output by 0.5. A multiplier and a matching processor that receives the delay output of the central receiver and the output of each multiplier in each sensor set and performs a matching process calculation,
It is characterized in that the reception sensitivity of the side lobe is reduced by the adaptive processing calculation of the adaptive processor to reduce noise.
【0031】[0031]
【発明の実施の形態】以下、本発明に係る面アレイマル
チチャンネル適合整相方法及び装置の実施の形態を図面
に従って説明する。BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a block diagram showing an embodiment of a method and an apparatus for phasing a surface array multi-channel adaptation according to the present invention.
【0032】図1は本発明の実施の形態の全体構成図、
図2は実施の形態における受波器配置を示す斜視図、図
3は音源の配置例を示す斜視図、図4は実施の形態での
適合処理計算による受信感度の計算結果の特性図であ
る。FIG. 1 is an overall configuration diagram of an embodiment of the present invention.
FIG. 2 is a perspective view showing an arrangement of the receivers in the embodiment, FIG. 3 is a perspective view showing an example of the arrangement of the sound sources, and FIG. 4 is a characteristic diagram of a calculation result of the receiving sensitivity by the adaptation processing calculation in the embodiment. .
【0033】図1及び図2において、水中又は空中にて
音波を発信し目標からの反響音を受信するか或いは目標
からの放射音を受信する音波センサーとしての受波器1
から13はすべてXY平面上にあり(つまり面状配置で
あり)、図2のように受波器2、1、3は第1組目、受
波器4、1、5は第2組目、受波器6、1、7は第3組
目、受波器8、1、9は第4組目、受波器10、1、1
1は第5組目、受波器12、1、13は第6組目の受波
器の組(3つの受波器からなる3チャンネル適合整相を
行う受波器の組)をそれぞれ示す。受波器1は各組共通
の受波器であり、各組の両端の受波器の中央に位置す
る。In FIG. 1 and FIG. 2, a receiver 1 as a sound wave sensor for transmitting a sound wave underwater or in the air and receiving a reverberation sound from a target or receiving a radiation sound from a target.
To 13 are all on the XY plane (that is, in a planar arrangement), as shown in FIG. 2, receivers 2, 1, and 3 are a first set, and receivers 4, 1, and 5 are a second set. , Receivers 6, 1, 7 are a third set, receivers 8, 1, 9 are a fourth set, and receivers 10, 1, 1,
Reference numeral 1 denotes a fifth set, and receivers 12, 1 and 13 denote a sixth set of receivers (a set of three receivers that perform three-channel adaptive phasing). . The receiver 1 is a receiver common to each set, and is located at the center of the receiver at both ends of each set.
【0034】図2は、説明を簡単にするため第1組と第
4組、第2組と第5組、第3組と第6組の受波器がそれ
ぞれ3本の直線上にある場合について例示したが、一般
的にはM(1≦M)本の直線があり、i(1≦i≦M)
本目の直線上にNi(1≦N i)個の受波器がある場合
を考える。FIG. 2 shows the first set and the first set for the sake of simplicity.
4 sets, 2nd set and 5th set, 3rd set and 6th set of receivers
Although the case where they are on three straight lines is illustrated,
Typically, there are M (1 ≦ M) straight lines, and i (1 ≦ i ≦ M)
N on the straight linei(1 ≦ N i) When there are receivers
think of.
【0035】線14は受波器2と受波器3の間隔、線1
5は受波器4と受波器5の間隔、線16は受波器6と受
波器7の間隔、線17は受波器8と受波器9の間隔、線
18は受波器10と受波器11の間隔、線19は受波器
12と受波器13の間隔を示す。Line 14 is the distance between receiver 2 and receiver 3, line 1
5 is the distance between the receiver 4 and the receiver 5, line 16 is the distance between the receiver 6 and the receiver 7, line 17 is the distance between the receiver 8 and the receiver 9, and line 18 is the receiver A line 19 indicates an interval between the receiver 10 and the receiver 11, and a line 19 indicates an interval between the receiver 12 and the receiver 13.
【0036】通常は直線同士のなす角21は全て等しく
するが、例えば、縦と横の間隔が異なる格子状に配置さ
れたアレイのように、等しくできない場合及び特に目的
があって等しくしない場合でも支障はない。M本の直線
同士のなす角を全て等しくする場合に直線同士のなす角
はπ/Mとなる。Normally, the angles 21 formed by the straight lines are all equal. However, even if the angles cannot be made equal to each other, for example, an array arranged in a lattice having different vertical and horizontal intervals, and even if there is a particular purpose, they are not equal. No problem. When the angles formed by the M straight lines are all equal, the angle formed by the straight lines is π / M.
【0037】i本目の直線上に前記受波器の組がNi個
のある場合、i本目の直線上の1組目の両端の受波器間
隔をdi1、2組目の両端の受波器間隔をdi2、j組
目の両端の受波器間隔をdij、Ni組目の両端の受波
器間隔をdiNiとすれば、di2=2×di1、d
i3=3×di1、…、dij=j×di1、…、di
Ni=Ni×di1のように1組目の受波器間隔の整数
倍とするが、整数倍にできない場合及び特に目的があっ
て整数倍にしない場合でも支障はない。[0037] If the i-th set of straight line on the wave receiver is a N i number, the receivers spacing at both ends of the first set on the straight line of the i-th across the d i1, 2 pair in receiving if the duplexer spacing d i2, j-th set receivers interval d ij of both ends of the receivers spacing across the N i th set as d iNi, d i2 = 2 × d i1, d
i3 = 3 × d i1, ... , d ij = j × d i1, ..., d i
Although it is set to an integral multiple of the receiver interval of the first set such as Ni = N i × d i1 , there is no problem even when the integral multiple cannot be used and when there is a special purpose and the integral multiple is not used.
【0038】M本の直線上の前記受波器の組において、
1本目の直線の両端の受波器間隔をd11、2本目の直
線の両端の受波器間隔をd21、i本目の直線の両端の
受波器間隔をdi1、M本目の直線の両端の受波器間隔
をdM1としたとき、通常はd11=d21=…=d
i1=…=dM1とするが、等しくできない場合及び特
に目的があって等しくしない場合でも支障はない。In the set of receivers on M straight lines,
The receiver interval at both ends of the first straight line is d 11 , the receiver interval at both ends of the second straight line is d 21 , the receiver interval at both ends of the i-th straight line is di 1 , and the M-th straight line is Assuming that the receiver spacing at both ends is d M1 , usually d 11 = d 21 =... = D
i1 = ... = but the d M1, trouble even if you do not equally if there is a case and, in particular, the purpose can not be equal not.
【0039】M本の直線上の前記受波器の組の個数Ni
は、通常はN1=N2=…=Ni=…=NMとするが、
等しくできない場合及び特に目的があって等しくしない
場合でも支障はない。The number N i of sets of the receivers on M straight lines
Is usually N 1 = N 2 = ... = N i = ... = but the N M,
There is no problem even if they cannot be made equal, or even if there is a particular purpose and they are not made equal.
【0040】各受波器は曲面上に配置することもある
が、この場合には遅延器で遅延を行い各受波器を平面上
の配置に置き換える。Each receiver may be arranged on a curved surface. In this case, the delay is delayed by a delay unit, and each receiver is replaced with an arrangement on a plane.
【0041】図1の全体構成図において、受波器1から
13で受波した音波を遅延器31から43で遅延を行い
所要の方向にメインローブ(メインビーム)を向ける
が、ここでは説明を簡単にするため各受波器は全て平面
上にあり、各遅延器の遅延量は零とすれば、本実施の形
態の装置は図2のZ軸方向にメインローブを持ち、その
他の方向はサイドローブ(サイドビーム)を持つ。第1
組目の受波器の組における両端の受波器2、3に対応し
た遅延器32、33の出力を加算器44で、第2組目の
受波器の組における両端の受波器4、5に対応した遅延
器34、35の出力を加算器45で、第3組目の受波器
の組における両端の受波器6、7に対応した遅延器3
6、37の出力を加算器46で、第4組目の受波器の組
における両端の受波器8、9に対応した遅延器38、3
9の出力を加算器47で、第5組目の受波器の組におけ
る両端の受波器10、11に対応した遅延器40、41
の出力を加算器48で、第6組目の受波器の組における
両端の受波器12、13に対応した遅延器42、43の
出力を加算器49でそれぞれ加算した後、それぞれ乗算
器50、51、52、53、54、55で0.5を乗ず
る。そして、各受波器の組の中央の受波器1に対応した
遅延器31の出力と、乗算器50、51、52、53、
54、55の出力を適合処理器56に入力し、適合処理
計算を行う(3チャンネルの受波器を持つ各組における
両端の2つの受波器の加算出力と中央の受波器出力の受
信指向性の差並びに位相の差に適合処理計算を適用す
る)。In the overall configuration diagram of FIG. 1, the sound waves received by the receivers 1 to 13 are delayed by the delay units 31 to 43 and the main lobe (main beam) is directed in a required direction. For the sake of simplicity, all the receivers are on a plane, and the delay amount of each delay unit is assumed to be zero. As a result, the apparatus of this embodiment has a main lobe in the Z-axis direction in FIG. Has side lobes (side beams). First
The outputs of the delay units 32 and 33 corresponding to the receivers 2 and 3 at both ends of the second set of receivers are added by the adder 44 to the receivers 4 at both ends of the second set of receivers. The outputs of the delay units 34 and 35 corresponding to the receivers 5 and 5 are added by an adder 45 to the delay units 3 corresponding to the receivers 6 and 7 at both ends in the third set of receivers.
The outputs of 6 and 37 are added by an adder 46 to delayers 38 and 3 corresponding to the receivers 8 and 9 at both ends in the fourth set of receivers.
The output of 9 is added by an adder 47 to delay units 40 and 41 corresponding to the receivers 10 and 11 at both ends in the fifth set of receivers.
The outputs of the delay units 42 and 43 corresponding to the receivers 12 and 13 at both ends of the sixth set of receivers are added by an adder 49, respectively. Multiply 0.5 by 50, 51, 52, 53, 54, 55. Then, the output of the delay unit 31 corresponding to the center receiver 1 of each set of receivers and the multipliers 50, 51, 52, 53,
The outputs of 54 and 55 are input to the adaptation processor 56, and an adaptation processing calculation is performed. (Reception of the sum output of the two receivers at both ends and the output of the center receiver in each set having a 3-channel receiver) The adaptation calculation is applied to the difference in directivity and the difference in phase).
【0042】図3において、球61、62、63はXY
平面上、YZ平面上、ZX平面上にある音源である。角
67、68、69は入射角度である。In FIG. 3, balls 61, 62 and 63 are XY
The sound source is on a plane, YZ plane, or ZX plane. Angles 67, 68 and 69 are incident angles.
【0043】図4は、Mが2で直線のなす角21がπ/
2、Nが1、音源と受波器1の距離22が両端の受波器
間隔dの約1.3倍、周波数が[音速/(2×d)]及
び[音速/d]の場合の本実施の形態の装置の入射角度
別の受信感度の計算結果である。グラフ71、72は音
源が61にある場合、グラフ73、74は音源が62に
ある場合、グラフ75、76は音源が63にある場合で
ある。本実施の形態の装置はサイドローブ方向の受信感
度を低下させることが分かる。FIG. 4 shows that when M is 2 and the angle 21 formed by the straight line is π /
2, N = 1, the distance 22 between the sound source and the receiver 1 is about 1.3 times the receiver distance d at both ends, and the frequency is [sound speed / (2 × d)] and [sound speed / d]. 9 is a calculation result of reception sensitivity for each incident angle of the device of the present embodiment. Graphs 71 and 72 show the case where the sound source is at 61, graphs 73 and 74 show the case where the sound source is at 62, and graphs 75 and 76 show the case where the sound source is at 63. It can be seen that the apparatus of the present embodiment reduces the reception sensitivity in the side lobe direction.
【0044】なお、Mが2、Nが1で受波器1、2、
3、4及び5で構成する装置における遅延器31、3
2、33、34及び35の出力x1k,x2k,
x3k,x4k,x5kは、下記の式(1)から式(5)で示
される。 x1k = a1k cos(ωt+ψ1k) …(1) x2k = a2k cos(ωt+ψ2k) …(2) x3k = a3k cos(ωt+ψ3k) …(3) x4k = a4k cos(ωt+ψ4k) …(4) x5k = a5k cos(ωt+ψ5k) …(5) (但し、k:k番目のサンプリング値、a1k〜
a5k:振幅成分、ω:音波の角周波数、ψ1k〜ψ
5k:位相成分)遅延器32と33の出力を加算器44
で加算し、乗算器50で0.5を乗じた値は式(6)で表
される。 x23k = 0.5×(x2k+x3k) …(6) 遅延器34と35の出力を加算器45で加算し、乗算器
51で0.5を乗じた値は式(7)で表される。 x45k = 0.5×(x4k+x5k) …(7) そして、遅延器31、乗算器50、51の出力値を適合
処理器56に入力し適合処理計算を行って上記の図4の
ような結果を得ることができる。Note that M is 2, N is 1, and the receivers 1, 2,.
Delay devices 31, 3 in the device composed of 3, 4 and 5
2 , 33, 34 and 35 outputs x 1k , x 2k ,
x 3k , x 4k , and x 5k are represented by the following equations (1) to (5). x1k = a1kcos (ωt + ψ1k ) (1) x2k = a2kcos (ωt + ψ2k ) (2) x3k = a3kcos (ωt + ψ3k ) ... (3) x4k = a4kcos (ωt) 4k ) ... (4) x5k = a5kcos (ωt + ψ5k ) ... (5) (where k: k-th sampling value, a1k to
a 5k : amplitude component, ω: angular frequency of sound wave, ψ 1k to ψ
5k : phase component) The outputs of the delay units 32 and 33 are added to an adder 44.
The value obtained by multiplying by 0.5 and multiplying by 0.5 by the multiplier 50 is expressed by Expression (6). x23k = 0.5 × ( x2k + x3k ) (6) The value obtained by adding the outputs of the delay units 34 and 35 by the adder 45 and multiplying by 0.5 by the multiplier 51 is represented by the following equation (7). . x 45k = 0.5 × (x 4k + x 5k ) (7) Then, the output values of the delay unit 31 and the multipliers 50 and 51 are input to the adaptation processing unit 56 and the adaptation processing calculation is performed to execute the processing shown in FIG. The result can be obtained.
【0045】前記適合処理計算は、請求項16,17及
び18又は19による計算方法やフロストの適合処理計
算を用いて行うことができる。The adaptation processing calculation can be performed using a calculation method according to claims 16, 17 and 18 or 19 or a frost adaptation processing calculation.
【0046】前記適合処理計算を、請求項16,17及
び18又は19による計算方法を用いて行う場合につい
て、以下に述べる。3チャンネルの受波器からなる前記
受波器の組における両端の2つの受波器の加算出力に
0.5を乗算した結果をy1とし、中央の受波器出力を
y2とする。適合処理器における各々の受波器の組につ
いての適合処理出力zは z = W1y1+W2y2 …(8) ここで、W1とW2は適合処理アルゴリズム(LMSア
ルゴリズムを用いる)におけるタップ重みであり、W1
とW2の和が1となるように選択する。つまり、両端の
2つの受波器の加算出力に0.5を乗算した結果(=y
1)に対するタップ重みW1と中央の受波器出力(=y
2)に対するタップ重みW2の和を常に1に保つことに
よって、目標からの受信信号に対する利得を常に1に保
つようにする。The case where the adaptation processing calculation is performed by using the calculation method according to claim 16, 17 and 18 or 19 will be described below. 3 channels of the result of multiplying 0.5 to the sum output of the two receivers at both ends in the receivers of the sets of wave receiver and y 1, the central receivers output y 2. The adaptation processing output z for each set of receivers in the adaptation processor is z = W 1 y 1 + W 2 y 2 (8) where W 1 and W 2 are adaptation processing algorithms (using the LMS algorithm) , And W 1
The sum of the W 2 is selected to be 1. That is, the result of multiplying the added output of the two receivers at both ends by 0.5 (= y
Tap weight W 1 and the center of the receiving transducer output to 1) (= y
By always keeping the 1 the sum of tap weights W 2 for 2), to keep the gain for the received signal from the target always 1.
【0047】適合処理アルゴリズムにおける誤差信号ε
を適合処理計算の出力zが最小になるように設定して式
(9)とすれば、 ε = −z …(9) タップ重みW1,W2は次式(10a)又は式(10b)で計算す
ることができる。 Wk+1 j =|Wk j+2μεkyk j| …(10a) Wk+1 j = Wk j+2μεkyk j …(10b) 但し、式(10a)はタップ重みを全て正の数を採用するよ
うにした場合、式(10b)はタップ重みを正負いずれの場
合でも採用する場合で、タップ重みW1を求めるときj
=1、タップ重みW2を求めるときj=2である。な
お、式(10b)のタップ重みを正負いずれの場合でも採用
する場合で、当該タップ重みが負になったときは当該タ
ップ重みの絶対値の上限値を定め、タップ重みが発散し
ないように設定するとよい。The error signal ε in the adaptation processing algorithm
Is set so that the output z of the adaptive processing calculation is minimized, and the equation
Assuming (9), ε = −z (9) The tap weights W 1 and W 2 can be calculated by the following equation (10a) or (10b). W k + 1 j = | W k j + 2με k y k j | ... (10a) W k + 1 j = W k j + 2με k y k j ... (10b) where the formula (10a) adopts the number of all the tap weights positive If you like to, if equation (10b) is employed in any case the tap weights positive and negative, when determining the tap weights W 1 j
= 1, and j = 2 when determining the tap weights W 2. In addition, when the tap weight of the equation (10b) is adopted in both positive and negative cases, when the tap weight becomes negative, the upper limit of the absolute value of the tap weight is determined, and the tap weight is set so as not to diverge. Good to do.
【0048】このように計算したタップ重みW1,W2
(LMSアルゴリズムにより求められたそれぞれの収束
値)に対して、音波入射方向がメインローブ方向となっ
て到来した目標からの受波信号が低減されないように、
W1とW2の和を常に1にするための拘束条件として次
式(11)を用いて再度タップ重みを計算し直す。つまり、
適合処理計算毎にタップ重みの和を求め、それぞれのタ
ップ重みをこの和で割り戻してタップ重みの合計を常に
1に保つ計算を式(11)によって行う。 W'k+1 j = Wk+1 j/(Wk+1 1+ Wk+1 2) …(11) ここで、W'k+1 jは割り戻し後の値、 W
k+1 jは割り戻し前のである(但し、j=1,
2)。The tap weights W 1 and W 2 thus calculated
With respect to (the respective convergence values obtained by the LMS algorithm), the received signal from the target arriving with the sound wave incident direction being the main lobe direction is not reduced.
The tap weight is calculated again using the following equation (11) as a constraint condition for always making the sum of W 1 and W 2 equal to 1. That is,
The sum of the tap weights is calculated for each adaptation process calculation, and each tap weight is divided by the sum to keep the total of the tap weights at 1 by using equation (11). W 'k + 1 j = W k + 1 j / (W k + 1 1 + W k + 1 2) ... (11) Here, W' k + 1 j are values after rebate, W
k + 1 j is before the rebate (where j = 1,
2).
【0049】このように、サイドローブ方向の雑音を低
減するように最適値に設定したタップ重みW1,W2を
用いて適合処理器56で各受波器の組を適合整相するこ
とで、従来の加算整相に比べて、サイドローブを低減で
きる。As described above, by using the tap weights W 1 and W 2 set to the optimal values so as to reduce the noise in the side lobe direction, the adaptive processor 56 adaptively phasing each receiver set. Thus, side lobes can be reduced as compared with the conventional addition phasing.
【0050】なお、サイドローブ方向からの音波が定常
的なものであれば、タップ重みは最適値に設定後一定と
なるが、定常的でなければタップ重みは最適値を求めて
変化する。If the sound wave from the side lobe direction is stationary, the tap weight becomes constant after being set to the optimum value, but if it is not stationary, the tap weight changes by finding the optimum value.
【0051】この実施の形態によれば、次の通りの効果
を得ることができる。According to this embodiment, the following effects can be obtained.
【0052】(1) 水中又は空中において、音波を発信
し目標からの反響音を受信するか或いは目標からの放射
音を受信する受波器出力を整相する場合において、等間
隔に配置された両端及び中央の合計3個の受波器からな
る組を複数設け、各受波器の組における両端の2つの受
波器の加算出力と中央の受波器の出力の受信指向性の差
並びに位相の差に適合処理計算を適用し、サイドローブ
方向から到来する雑音を低減可能である。(1) Underwater or in the air, when sound waves are transmitted and a reverberation sound from a target is received or when a receiver output for receiving a radiation sound from a target is phased, the receivers are arranged at equal intervals. A plurality of sets each having a total of three receivers at both ends and a center are provided, and a difference in reception directivity between the sum output of the two receivers at both ends and the output of the center receiver in each set of receivers; By applying the adaptation processing calculation to the phase difference, it is possible to reduce noise coming from the side lobe direction.
【0053】(2) 前記受波器の組を複数本の直線上に
面状配置することで、多様な方向のサイドローブに対し
て抑圧機能を持たせることができる。(2) By arranging the sets of the receivers in a plane on a plurality of straight lines, a function of suppressing side lobes in various directions can be provided.
【0054】(3) 前記受波器の組における両端の2つ
の受波器出力の和の0.5倍の値と中央の受波器出力と
の和が最小になるように、適合処理計算を行えばよく、
構成が複雑化することがない。(3) Adaptation processing calculation so that the sum of the value of 0.5 times the sum of the outputs of the two receivers at both ends in the set of receivers and the output of the center receiver is minimized. Should be done,
The configuration is not complicated.
【0055】(4) 前記適合処理計算は、請求項16,
17及び18又は19による計算方法やフロストの適合
処理を用いて行うことができる。(4) The adaptation processing calculation is performed according to claim 16,
The calculation can be performed by using the calculation method 17 or 18 or 19 or the frost adaptation processing.
【0056】(5) 前記適合処理計算を、請求項16,
17及び18又は19による計算方法を用いて行う場合
に、前記受波器の組における両端の2つの受波器の加算
出力に対するタップ重みと中央の受波器出力に対するタ
ップ重みの和を常に1に保つことによって、目標からの
受信信号に対する利得を常に1に保つことができる。具
体的には、前記受波器の組における両端の2つの受波器
の加算出力に対するタップ重みと中央の受波器出力に対
するタップ重みの和を常に1に保つために、適合処理計
算毎にタップ重みの和を求め、それぞれのタップ重みを
この和で割り戻してタップ重みの合計を常に1に保つ。(5) The adaptation processing calculation is performed according to claim 16,
When using the calculation method of 17 and 18 or 19, the sum of the tap weights for the added outputs of the two receivers at both ends of the set of receivers and the tap weights for the center receiver output is always 1 , The gain for the received signal from the target can always be kept at 1. Specifically, in order to always keep the sum of the tap weight for the added output of the two receivers at both ends in the set of receivers and the tap weight for the center receiver output at 1, The sum of the tap weights is obtained, and each tap weight is divided by the sum to keep the total of the tap weights at one.
【0057】なお、本発明は音波の反射又は放射を受波
するソーナー装置等の整相処理に有用であり、魚群探知
機にも適用可能である。The present invention is useful for a phasing process of a sonar device or the like that receives reflection or radiation of a sound wave, and is applicable to a fish finder.
【0058】また、受波器の一部又は全部を振動センサ
ーに置き換え、伝達関数(振動を音響に変換するための
もの)による処理を行い、適合処理計算を行う構成とす
ることも可能である。It is also possible to replace a part or the whole of the receiver with a vibration sensor, perform a process using a transfer function (for converting vibration into sound), and perform an adaptive process calculation. .
【0059】また、受波器の全部を振動センサーに置き
換え、振動に対して適合処理計算を行う構成とすること
も可能である。It is also possible to replace all of the receivers with a vibration sensor, and to perform an adaptive processing calculation on the vibration.
【0060】以上本発明の実施の形態について説明して
きたが、本発明はこれに限定されることなく請求項の記
載の範囲内において各種の変形、変更が可能なことは当
業者には自明であろう。Although the embodiments of the present invention have been described above, it is obvious to those skilled in the art that the present invention is not limited to the embodiments and various modifications and changes can be made within the scope of the claims. There will be.
【0061】[0061]
【発明の効果】以上説明したように、本発明に係る面ア
レイマルチチャンネル適合整相方法及び装置によれば、
中央の受波器と両端の2つの受波器からなる受波器の組
を複数配置し、複数配置された前記センサーの組がすべ
て中央のセンサーを共有するようにし、前記受波器の組
における中央の受波器と両端の受波器との指向性の差及
び位相差を利用し、これに適合処理計算を適用すること
により、雑音源がサイドビーム方向であれば、いずれの
距離に存在しても、いずれの方向に存在しても、常に雑
音レベルが最小になるような最適なシェーディング係数
の選択を行うことができる。この結果、雑音源の抑圧効
果を向上させることが可能となる。As described above, according to the surface array multi-channel adaptive phasing method and apparatus according to the present invention,
A plurality of sets of receivers including a center receiver and two receivers at both ends are arranged, and a plurality of sets of the sensors are arranged.
Using a difference in directivity and a phase difference between the center receiver and the receivers at both ends in the set of receivers, and applying an adaptive processing calculation to this. Therefore, if the noise source is in the side beam direction, it is possible to always select the optimal shading coefficient that minimizes the noise level regardless of the distance or the direction in which the noise source exists. . As a result, it is possible to improve the noise source suppression effect.
【図1】本発明に係る面アレイマルチチャンネル適合整
相方法及び装置の実施の形態の全体構成図(但しM=
3、N=2の場合)である。FIG. 1 is an overall configuration diagram of an embodiment of a surface array multi-channel adaptive phasing method and apparatus according to the present invention (where M =
3, N = 2).
【図2】実施の形態における受波器配置を説明するため
の斜視図(但しM=3、N=2の場合)である。FIG. 2 is a perspective view for explaining a receiver arrangement in the embodiment (provided that M = 3 and N = 2).
【図3】実施の形態における受波器に入射する音波の音
源の位置を説明するための斜視図である。FIG. 3 is a perspective view for explaining a position of a sound source of a sound wave incident on the receiver in the embodiment.
【図4】実施の形態において、M=2、N=1の場合に
図3に示す音源から音波が入射したときの受信感度の計
算結果の1例を説明するための特性図である。FIG. 4 is a characteristic diagram illustrating an example of a calculation result of reception sensitivity when a sound wave enters from the sound source illustrated in FIG. 3 when M = 2 and N = 1 in the embodiment.
1〜13 受波器 14 1組目の両端の受波器の間隔 15 2組目の両端の受波器の間隔 16 3組目の両端の受波器の間隔 17 4組目の両端の受波器の間隔 18 5組目の両端の受波器の間隔 19 6組目の両端の受波器の間隔 21 第1本目の直線と第2本目の直線のなす角 22 中心受波器と音源との距離 23 XY平面上への音源の投影位置と中心受波器の作
る直線とX軸のなす角 24 音源と中心受波器の作る直線とZ軸のなす角 31〜43 遅延器 44〜49 加算器 50〜55 番乗算器 56 適合処理器 61 XY平面上にある音源 62 YZ平面上にある音源 63 ZX平面上にある音源 64 受波器1と音源61との距離 65 受波器1と音源62との距離 66 受波器1と音源63との距離 67 X軸と線64のなす角 68 Z軸と線65のなす角 69 Z軸と線66のなす角 71 音源がXY平面にある場合の入射角度別受信感度 72 音源がXY平面にある場合の入射角度別受信感度 73 音源がYZ平面にある場合の入射角度別受信感度 74 音源がYZ平面にある場合の入射角度別受信感度 75 音源がZX平面にある場合の入射角度別受信感度 76 音源がZX平面にある場合の入射角度別受信感度1-13 Receiver 14 Spacing between receivers at both ends of first set 15 Spacing between receivers at both ends of second set 16 Spacing between receivers at both ends of third set 17 Receiver at both ends of fourth set Spacing between receivers 18 Spacing between receivers at both ends of fifth set 19 Spacing between receivers at both ends of sixth set 21 Angle between first straight line and second straight line 22 Center receiver and sound source 23 Angle between the projection position of the sound source on the XY plane and the straight line formed by the center receiver and the X axis 24 The angle formed by the sound source and the straight line formed by the center receiver and the Z axis 31 to 43 Delay device 44 to Reference Signs List 49 adder 50-55 multiplier 56 adaptation processor 61 sound source on XY plane 62 sound source on YZ plane 63 sound source on ZX plane 64 distance between receiver 1 and sound source 61 65 receiver 1 66 The distance between the receiver 1 and the sound source 63 67 The angle between the X axis and the line 64 68 The Z axis Angle formed by 65 69 Angle formed between Z axis and line 66 71 Reception sensitivity by incident angle when sound source is on XY plane 72 Reception sensitivity by incidence angle when sound source is on XY plane 73 When reception source is on YZ plane Receiving sensitivity by incident angle 74 Receiving sensitivity by incident angle when the sound source is on the YZ plane 75 Receiving sensitivity by incident angle when the sound source is on the ZX plane 76 Receiving sensitivity by incident angle when the sound source is on the ZX plane
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 平6−18663(JP,A) 特開 平10−207490(JP,A) 特許2723865(JP,B2) 特許3005680(JP,B2) 特許2720845(JP,B2) (58)調査した分野(Int.Cl.7,DB名) G01S 3/80 - 3/86 G01S 5/18 - 5/30 G01S 7/52 - 7/64 G01S 15/00 - 15/96 ──────────────────────────────────────────────────続 き Continuation of the front page (56) References JP-A-6-18663 (JP, A) JP-A-10-207490 (JP, A) Patent 2723865 (JP, B2) Patent 30000568 (JP, B2) Patent 2720845 (JP, B2) (58) Fields surveyed (Int. Cl. 7 , DB name) G01S 3/80-3/86 G01S 5/18-5/30 G01S 7 /52-7/64 G01S 15/00- 15/96
Claims (19)
又は振動を受信するか或いは目標からの放射音又は振動
を受信する音波又は振動のセンサー出力の整相方法にお
いて、 直線上に等間隔に配置された3つのセンサーからなるセ
ンサーの組を複数配置し、複数配置された前記センサー
の組がすべて中央のセンサーを共有しており、各組にお
ける両端の2つのセンサーの加算出力と中央のセンサー
出力の受信指向性の差並びに位相の差に適合処理計算を
適用し、サイドローブ方向から到来する雑音を低減させ
ることを特徴とする面アレイマルチチャンネル適合整相
方法。1. A method for phasing a sensor output of a sound wave or a vibration for transmitting a sound wave or a vibration and receiving a reverberation sound or a vibration from a target or receiving a sound or a vibration radiated from a target, comprising the steps of: A plurality of sets of sensors consisting of three sensors arranged in a plurality of sensors, wherein the plurality of sensors are arranged
All of the sets share the central sensor, apply the adaptive processing calculation to the difference between the reception directivity and the phase difference of the two sensors at both ends and the center sensor output in each set, and apply the sidelobe direction. Surface array multi-channel adaptive phasing method, characterized in that noise coming from the array is reduced.
本の直線上に配置されている請求項1記載の面アレイマ
ルチチャンネル適合整相方法。2. The method according to claim 1, wherein a set of the plurality of sensors is M (1 ≦ M).
2. The surface array multi-channel adaptive phasing method according to claim 1, wherein the phasing method is arranged on a straight line.
線同士のなす角を等しい角度とした請求項2記載の面ア
レイマルチチャンネル適合整相方法。3. The surface array multi-channel adaptive phasing method according to claim 2, wherein the M straight lines are arranged in a plane, and angles formed by the straight lines are equal to each other.
線同士のなす角の一部又は全部を異なる角度とした請求
項2記載の面アレイマルチチャンネル適合整相方法。4. The surface array multi-channel adaptive phasing method according to claim 2, wherein the M straight lines have a planar arrangement, and some or all of the angles formed by the straight lines are different.
センサーの組をNi(1≦Ni)個配置した請求項2記
載の面アレイマルチチャンネル適合整相方法。5. The surface array multi-channel adaptive phasing method according to claim 2 , wherein N i (1 ≦ N i ) sets of sensors are arranged on the i-th (1 ≦ i ≦ M) straight line.
請求項5記載の面アレイマルチチャンネル適合整相方
法。6. N 1 = N 2 = ... = N i = ... = N M a surface array multichannel adapted phasing method according to claim 5, wherein the.
求項5記載の面アレイマルチチャンネル適合整相方法。7. The surface array multi-channel adaptive phasing method according to claim 5, wherein some or all of the Nis are different numbers.
ー間隔をdi1、2組目の両端のセンサー間隔を
di2、j組目の両端のセンサー間隔をdij、Ni組
目の両端のセンサー間隔をdiNiとしたとき、di2
=2×di1、di3=3×di1、…、dij=j×
di1、…、diNi=Ni×di1のように1組目の
両端のセンサー間隔の整数倍とした請求項5記載の面ア
レイマルチチャンネル適合整相方法。8. i-th straight line of the first set of sensor spacing d i1 ends, the second set of sensors intervals across the sensor spacing d i2, j-th group at both ends d ij, N i sets when the sensor interval of both ends of the eye was d iNi, d i2
= 2 × d i1 , d i3 = 3 × d i1 ,..., D ij = j ×
6. The surface array multi-channel adaptive phasing method according to claim 5 , wherein di 1 ,..., di Ni = N i × d i1 is an integral multiple of the sensor interval between both ends of the first set.
ー間隔をdij(1≦j≦Ni)としたとき、前記セン
サー間隔dijの一部又は全部を1組目の両端のセンサ
ー間隔の整数倍とは異なる間隔にした請求項5記載の面
アレイマルチチャンネル適合整相方法。9. When a sensor interval between both ends of a j-th set on an i-th straight line is d ij (1 ≦ j ≦ N i ), a part or all of the sensor interval d ij is set to both ends of the first set. 6. The surface array multi-channel adaptive phasing method according to claim 5, wherein the interval is different from an integral multiple of the sensor interval.
d11、2本目の直線の両端のセンサー間隔をd21、
i本目の直線の両端のセンサー間隔をdi1、M本目の
直線の両端のセンサー間隔をdM1としたとき、d11
=d21=…=di1=…=dM1とした請求項5記載
の面アレイマルチチャンネル適合整相方法。10. The sensor interval at both ends of the first straight line is d 11 , the sensor interval at both ends of the second straight line is d 21 ,
When the sensor interval at both ends of the i-th straight line is d i1 and the sensor interval at both ends of the M-th straight line is d M1 , d 11
= D 21 = ... = d i1 = ... = d M1 to the surface array multichannel adapted phasing method according to claim 5, wherein the.
di1(1≦i≦M)としたとき、前記センサー間隔d
i1の一部又は全部を異なる間隔にした請求項5記載の
面アレイマルチチャンネル適合整相方法。11. When the sensor interval at both ends of the i-th straight line is d i1 (1 ≦ i ≦ M), the sensor interval d
6. The surface array multi-channel adaptive phasing method according to claim 5 , wherein a part or the whole of i1 is set at different intervals.
のセンサー出力の和の0.5倍の値と中央のセンサー出
力との和が最小になるように、適合処理計算を行う請求
項1から11のいずれかに記載の面アレイ3チャンネル
適合整相方法。12. As the sum of 0.5 times the value and the center of the sensor outputs of the two sum of the sensor output at both ends in the set of the sensor is minimized, claim 1 for adaptation process calculation 11 The phase-matching method for three-channel array according to any one of the above.
ズムを用いて行う請求項12記載の面アレイマルチチャ
ンネル適合整相方法。13. The surface array multi-channel adaptive phasing method according to claim 12 , wherein the adaptive processing calculation is performed using an LMS algorithm.
処理を用いて行う請求項13記載の面アレイマルチチャ
ンネル適合整相方法。14. The surface array multi-channel adaptive phasing method according to claim 13 , wherein the adaptive processing calculation is performed using frost adaptive processing.
ズムを用いて行う場合に、前記センサーの組における両
端の2つのセンサーの加算出力に対するタップ重みと中
央のセンサー出力に対するタップ重みの和を常に1に保
つことによって、目標からの受信信号に対する利得を常
に1に保つ請求項13記載の面アレイマルチチャンネル
適合整相方法。15. When the adaptation processing calculation is performed using an LMS algorithm, the sum of the tap weight for the added output of the two sensors at both ends in the set of sensors and the tap weight for the central sensor output is always set to 1. 14. The surface array multi-channel adaptive phasing method according to claim 13 , wherein the gain for the received signal from the target is always kept at 1 by keeping.
のセンサーの加算出力に対するタップ重みと中央のセン
サー出力に対するタップ重みの和を常に1に保つため
に、適合処理計算毎にタップ重みの和を求め、それぞれ
のタップ重みをこの和で割り戻してタップ重みの合計を
常に1に保つ請求項15記載の面アレイマルチチャンネ
ル適合整相方法。16. In order to always keep the sum of the tap weights for the added outputs of the two sensors at both ends in the set of sensors and the tap weights for the central sensor output at 1, the sum of the tap weights is calculated for each adaptation process calculation. 16. The surface array multi-channel adaptive phasing method according to claim 15 , wherein each tap weight is divided by the sum to keep a total of tap weights at one.
の計算において、前記タップ重みとして正負いずれの数
値も選択し、前記タップ重みが負になった場合に当該タ
ップ重みの絶対値の上限値を定める請求項13,15又
は16記載の面アレイマルチチャンネル適合整相方法。17. The calculation of tap weights of the LMS algorithm, wherein both positive and negative numerical values are selected as the tap weights, and when the tap weights become negative, an upper limit value of the absolute value of the tap weights is determined. 13, 15 or
Is a surface array multi-channel adaptive phasing method according to item 16 .
の計算において、前記タップ重みが負になったら、その
絶対値を取り、正の数値のみを選択する請求項13,1
5又は16記載の面アレイマルチチャンネル適合整相方
法。18. The calculation of the tap weights of the LMS algorithm, if the tap weight is negative, claim that takes the absolute value, selects only positive numbers 13, 1
17. The surface array multi-channel adaptive phasing method according to 5 or 16 .
音又は振動を受信するか或いは目標からの放射音又は振
動を受信する音波又は振動のセンサー出力の整相装置に
おいて、直線上に等間隔に配置された3つのセンサーからなるセ
ンサーの組が複数配置されかつ複数の前記センサーの組
がすべて中央のセンサーを共有しているとともに、 各センサー出力をそれぞれ遅延させる遅延器と、各セン
サーの組における両端の2つのセンサーの遅延出力をそ
れぞれ加算する加算器と、各加算器出力に0.5をそれ
ぞれ乗じる乗算器と、各センサーの組における中央の受
波器の遅延出力と各乗算器出力とを受けて適合処理計算
を行う適合処理器とを備え、 該適合処理器の適合処理計算によりサイドローブの受信
感度を下げて雑音を低減させることを特徴とする面アレ
イマルチチャンネル適合整相装置。19. The phasing device sensor output of acoustic or vibration receiving the radiated sound or vibration from or targeted to receive the reflected sound or vibration from the ultrasonic or outgoing Goal vibration, equally spaced on a straight line Consisting of three sensors
A plurality of sensor sets and a plurality of the sensor sets
Share the central sensor, delay the outputs of each sensor, add the delay outputs of the two sensors at each end of each sensor set, and add 0 to each adder output. .5, and an adaptation processor that receives the delayed output of the central receiver and the output of each multiplier in each set of sensors and performs an adaptation processing calculation. A plane array multi-channel adaptive phasing apparatus characterized in that the reception sensitivity of side lobes is reduced by calculation to reduce noise.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP03267999A JP3208424B2 (en) | 1999-02-10 | 1999-02-10 | Surface array multi-channel adaptive phasing method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP03267999A JP3208424B2 (en) | 1999-02-10 | 1999-02-10 | Surface array multi-channel adaptive phasing method and apparatus |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2000230973A JP2000230973A (en) | 2000-08-22 |
| JP3208424B2 true JP3208424B2 (en) | 2001-09-10 |
Family
ID=12365571
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP03267999A Expired - Lifetime JP3208424B2 (en) | 1999-02-10 | 1999-02-10 | Surface array multi-channel adaptive phasing method and apparatus |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3208424B2 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6807125B2 (en) * | 2002-08-22 | 2004-10-19 | International Business Machines Corporation | Circuit and method for reading data transfers that are sent with a source synchronous clock signal |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2720845B2 (en) | 1994-09-01 | 1998-03-04 | 日本電気株式会社 | Adaptive array device |
| JP2723865B2 (en) | 1995-11-20 | 1998-03-09 | 防衛庁技術研究本部長 | Spatial filter and azimuth analysis apparatus and method |
-
1999
- 1999-02-10 JP JP03267999A patent/JP3208424B2/en not_active Expired - Lifetime
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2720845B2 (en) | 1994-09-01 | 1998-03-04 | 日本電気株式会社 | Adaptive array device |
| JP2723865B2 (en) | 1995-11-20 | 1998-03-09 | 防衛庁技術研究本部長 | Spatial filter and azimuth analysis apparatus and method |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2000230973A (en) | 2000-08-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ahmad et al. | Synthetic aperture beamformer for imaging through a dielectric wall | |
| US8098842B2 (en) | Enhanced beamforming for arrays of directional microphones | |
| JP6737464B2 (en) | Test system, waveform simulator device, test method and program | |
| Stergiopoulos | Implementation of adaptive and synthetic-aperture processing schemes in integrated active-passive sonar systems | |
| US20060114148A1 (en) | Robust optimal shading scheme for adaptive beamforming with missing sensor elements | |
| CN101779140A (en) | Method for measuring angle of arrival of coherent source using spatial smoothing on arbitrary sensor network | |
| CN106855622A (en) | A kind of angle-measuring method of phased array at subarray level radar | |
| JPWO2018181201A1 (en) | Transmission device, reception device, transmission method and reception method | |
| US6058075A (en) | System for canceling interferers from broadband active sonar signals using adaptive beamforming methods | |
| Hao et al. | An optimization method for frequency-invariant beamforming with arbitrary sensor arrays | |
| JP5161474B2 (en) | Unwanted wave suppression device | |
| JP3396798B2 (en) | Target position localization method | |
| US20050088337A1 (en) | Vertically stacked turnstile array | |
| JP3208424B2 (en) | Surface array multi-channel adaptive phasing method and apparatus | |
| JP4072149B2 (en) | Distributed aperture antenna device | |
| Bi et al. | Design of a robust steerable differential beamformer with linear acoustic vector sensor arrays | |
| Tolstoy | Applications of matched-field processing to inverse problems in underwater acoustics | |
| Bertilone et al. | Array gain for a cylindrical array with baffle scatter effects | |
| Edelson et al. | Limitations on the overlap-correlator method imposed by noise and signal characteristics | |
| CN1875292A (en) | Resolving directional information in sonar arrays | |
| JP6088165B2 (en) | Detection device, detection method, and detection program | |
| Pal et al. | Frequency invariant MVDR beamforming without filters and implementation using MIMO radar | |
| Huang et al. | Robust adaptive beamforming for large-scale arrays | |
| JP3005680B1 (en) | Three-channel adaptive phasing method and apparatus | |
| Rashida et al. | High resolution wideband acoustic beamforming and underwater target localization using 64-element linear hydrophone array |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |