JP2000323915A - Interference wave suppressing device - Google Patents
Interference wave suppressing deviceInfo
- Publication number
- JP2000323915A JP2000323915A JP11181767A JP18176799A JP2000323915A JP 2000323915 A JP2000323915 A JP 2000323915A JP 11181767 A JP11181767 A JP 11181767A JP 18176799 A JP18176799 A JP 18176799A JP 2000323915 A JP2000323915 A JP 2000323915A
- Authority
- JP
- Japan
- Prior art keywords
- antenna
- equation
- signal
- interference wave
- auxiliary
- 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.)
- Pending
Links
Landscapes
- Noise Elimination (AREA)
- Radio Transmission System (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は、干渉波抑圧装置に
関し、特にアンテナのサイドローブから入射する干渉波
を抑圧する干渉波抑圧装置に関するものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an interference wave suppressing device, and more particularly to an interference wave suppressing device for suppressing an interference wave incident from a side lobe of an antenna.
【0002】[0002]
【従来の技術】アダプティブアンテナは干渉波の入射方
向にアンテナパターンの零点を向けるため、所望信号の
受信に支障となる干渉信号の抑圧に有効である。このア
ダプティブアンテナの荷重計算に用いる適応アルゴリズ
ムには、信号処理演算量が少なく荷重係数ができるだけ
速く最適値に収束することが望まれる。2. Description of the Related Art An adaptive antenna directs a zero point of an antenna pattern in an incident direction of an interference wave, and thus is effective for suppressing an interference signal which hinders reception of a desired signal. It is desired that the adaptive algorithm used for the load calculation of the adaptive antenna has a small amount of signal processing operation and the weight coefficient converges to the optimum value as quickly as possible.
【0003】SMI(Sample Matrix Inversion)法
は、受信信号をディジタル信号に変換して処理を行うこ
とを前提にして、アレーアンテナにおける各素子アンテ
ナ受信信号(ディジタル同相・直交信号)から相関行列
と相互相関ベクトルを計算して連立1次方程式を構成
し、それを直接解くことにより荷重係数を求める。少な
いサンプル数で良好な解を得ることができる。以下の説
明においては各信号は等価低域表現された離散時間信号
とする。The SMI (Sample Matrix Inversion) method is based on the premise that a received signal is converted into a digital signal for processing, and a correlation matrix and a cross-correlation matrix are obtained from each element antenna received signal (digital in-phase / quadrature signal) in an array antenna. A simultaneous linear equation is constructed by calculating a correlation vector, and a weighting factor is obtained by directly solving it. A good solution can be obtained with a small number of samples. In the following description, each signal is a discrete-time signal expressed in an equivalent low-frequency range.
【0004】SMI法は、受信信号をディジタル信号に
変換して処理するのであればサイドローブキャンセラに
対しても適用可能である。図6はサイドローブキャンセ
ラにSMI法を適用して荷重係数を求める場合の構成図
である。ここでは補助アンテナ数は2としている。図6
の機能と構成について説明する。d(k)は入射方向が
既知の所望波、s1(k)、s2(k)は入射方向が未知
の干渉波である。これらは互いに相関がないものとす
る。[0004] The SMI method can be applied to a side lobe canceller as long as the received signal is converted into a digital signal and processed. FIG. 6 is a configuration diagram in a case where a load coefficient is obtained by applying the SMI method to the side lobe canceller. Here, the number of auxiliary antennas is two. FIG.
The function and configuration of will be described. d (k) is a desired wave whose incident direction is known, and s 1 (k) and s 2 (k) are interference waves whose incident direction is unknown. These are assumed to be uncorrelated.
【0005】100は所望波方向にビームを向けた主ア
ンテナ、101a,101bは指向性が等方性に近い補
助アンテナ、1a,1b,1cは受信機、2a,2b,
2cはA/D変換器である。A/D変換器2a,2b,
2cの出力信号をそれぞれx 0(k),x1(k),x2
(k)とする。これは実部が同相成分,虚部が直交成分
のディジタル複素信号とする。[0005] Reference numeral 100 denotes a main beam for directing a beam in a desired wave direction.
The antennas 101a and 101b have complementary directivities close to isotropic.
Auxiliary antennas, 1a, 1b, 1c are receivers, 2a, 2b,
2c is an A / D converter. A / D converters 2a, 2b,
The output signals of 2c are x 0(K), x1(K), xTwo
(K). This means that the real part is the in-phase component and the imaginary part is the quadrature component
Digital complex signal.
【0006】補助アンテナ出力に対応するA/D変換器
2b,2cの出力信号x1(k),x2(k)に、荷重係
数重み付け手段である荷重係数器3a,3bにより値が
それぞれw1,w2の重み付けを行いそれらを加算する。
加算した結果得られた信号yA(k)は式(1)のよう
になる。荷重係数w1,w2が適切に決められているな
ら、yA(k)は主アンテナで受信された干渉信号のレ
プリカになっている。従って、yA(k)を主アンテナ
出力に対応するA/D変換器2aの出力信号x0(k)
から差し引くことによって、サイドローブキャンセラ出
力信号z(k)は干渉信号が抑圧された信号となる。The output signals x 1 (k) and x 2 (k) of the A / D converters 2 b and 2 c corresponding to the outputs of the auxiliary antennas have values w respectively by weighting means 3 a and 3 b as weighting means. 1 and w 2 are weighted and added.
The signal y A (k) obtained as a result of the addition is as shown in equation (1). If the weighting factors w 1 and w 2 are properly determined, y A (k) is a replica of the interference signal received by the main antenna. Therefore, y A (k) is the output signal x 0 (k) of the A / D converter 2a corresponding to the main antenna output.
, The side lobe canceller output signal z (k) becomes a signal in which the interference signal is suppressed.
【0007】[0007]
【数1】 (Equation 1)
【0008】次に,図6におけるSMI法による荷重計
算法を説明する。簡単のため電波環境は定常とする。S
MI法は、例えば従来技術文献1 小川、菊間著:“ア
ダプティブアンテナ理論の進展と今後の展望”(電子情
報通信学会論文誌B−II,vol. J75-B-II,No. 11,pp.
721-732,1992年11月)に記述がある。相関行列・相互
相関ベクトル計算手段10aでは式(2)の相関行列R
と式(3)の相互相関ベクトルpを計算する。Next, a load calculation method by the SMI method in FIG. 6 will be described. For simplicity, the radio wave environment is assumed to be stationary. S
The MI method is described in, for example, Prior Art Document 1 by Ogawa and Kikuma: "Progress and Future Prospects of Adaptive Antenna Theory" (Transactions of the Institute of Electronics, Information and Communication Engineers, B-II, vol.
721-732, November 1992). In the correlation matrix / cross-correlation vector calculation means 10a, the correlation matrix R
And the cross-correlation vector p of equation (3) is calculated.
【0009】[0009]
【数2】 (Equation 2)
【0010】[0010]
【数3】 (Equation 3)
【0011】式(2)の相関行列Rと式(3)の相互相
関ベクトルpの各要素は以下のように計算する。ここで
KはSMI法で使用するサンプル数、肩文字*は複素共
役を表す。相関行列Rはエルミート行列である。Each element of the correlation matrix R of the equation (2) and the cross-correlation vector p of the equation (3) is calculated as follows. Here, K represents the number of samples used in the SMI method, and the superscript * represents a complex conjugate. The correlation matrix R is a Hermitian matrix.
【0012】[0012]
【数4】 (Equation 4)
【0013】[0013]
【数5】 (Equation 5)
【0014】荷重係数ベクトルWを式(6)のように定
義すると、SMI法で解くべき方程式は式(7)のよう
になる。If the load coefficient vector W is defined as in equation (6), the equation to be solved by the SMI method is as in equation (7).
【0015】[0015]
【数6】 (Equation 6)
【0016】[0016]
【数7】 (Equation 7)
【0017】式(7)の解(荷重係数)は、式(8)の
ようになる。式(7)は、ガウス消去法などで数値的に
解くことができるが、補助アンテナ数が3以下の場合は
容易に逆行列を直接解析的に計算できる。その場合は数
値算法により方程式を解く必要はなくなる。補助アンテ
ナ数が2の場合の解は、式(9)のようになる。荷重計
算手段15は、相関行列・相互相関ベクトル計算手段1
0aで計算した行列とベクトルから、式(8)もしくは
式(9)を計算して荷重を求め、その値を荷重係数器3
a,3bに出力する。The solution (load coefficient) of equation (7) is as shown in equation (8). Equation (7) can be numerically solved by a Gaussian elimination method or the like, but when the number of auxiliary antennas is three or less, the inverse matrix can be easily calculated directly and analytically. In such a case, there is no need to solve the equation by a numerical algorithm. The solution when the number of auxiliary antennas is 2 is as shown in Expression (9). The load calculating means 15 includes a correlation matrix / cross-correlation vector calculating means 1
Equation (8) or Equation (9) is calculated from the matrix and the vector calculated at 0a to obtain the load, and the value is calculated by the load coefficient unit 3
a and 3b.
【0018】[0018]
【数8】 (Equation 8)
【0019】[0019]
【数9】 (Equation 9)
【0020】干渉波数が補助アンテナ数と一致する場合
は、式(8)、式(9)から計算された荷重を持つ図6
のサイドローブキャンセラは、良好な干渉波抑圧性能を
持つ。ところが干渉波数が補助アンテナ数より少ない場
合、あるいは干渉波数が補助アンテナ数と一致しても干
渉波数が補助アンテナ数より少ないとみなされる場合
(補助アンテナ101aに対する補助アンテナ101b
における干渉波s1(k)の受信位相差と干渉波s
2(k)の受信位相差との2つの位相差の差が2πの整
数倍の場合)には、式(8)、式(9)から計算された
荷重を持つ図6のサイドローブキャンセラは、干渉波を
抑圧できてもサイドローブキャンセラ出力における信号
対雑音電力比が低下してしまうことがある。それは、1
つには主アンテナで受信した所望波を抑圧してしまうこ
と、2つには荷重係数ベクトルWのノルム||W||が大き
くなってしまい、補助アンテナ101a,101bに接
続された受信機で発生する雑音が荷重係数w1,w2によ
り増幅されてしまうからである。特に荷重係数ベクトル
Wのノルム||W||が大きくなる理由について記されてい
る従来技術文献はあまりないと思われるので、上記の現
象について補助アンテナ数が2の場合について説明す
る。When the number of interference waves matches the number of auxiliary antennas, FIG. 6 having the weight calculated from equations (8) and (9)
Has good interference wave suppression performance. However, when the number of interference waves is smaller than the number of auxiliary antennas, or when the number of interference waves is considered to be smaller than the number of auxiliary antennas even when the number of interference waves matches the number of auxiliary antennas (the auxiliary antenna 101b with respect to the auxiliary antenna 101a)
Phase difference of interference wave s 1 (k) and interference wave s
In the case where the difference between the two phase differences from the 2 (k) reception phase difference is an integral multiple of 2π), the side lobe canceller in FIG. 6 having the weight calculated from Expressions (8) and (9) is However, even if the interference wave can be suppressed, the signal-to-noise power ratio at the output of the side lobe canceller may decrease. It is 1
One is that the desired wave received by the main antenna is suppressed, and the other is that the norm || W || of the weighting factor vector W becomes large, and the receiver connected to the auxiliary antennas 101a and 101b This is because the generated noise is amplified by the weight coefficients w 1 and w 2 . In particular, it is considered that there are not many prior art documents describing the reason why the norm || W || of the load coefficient vector W becomes large. Therefore, the above phenomenon will be described in the case where the number of auxiliary antennas is two.
【0021】定常状態の場合を考える。定常状態の場
合、式(2)の相関行列Rと式(3)の相互相関ベクト
ルpの各要素は和ではなく、式(10)と式(11)の
ような集合平均で表現される(i,j=1,2)。この
とき、式(8)はサイドローブキャンセラ出力信号平均
電力を最小にする意味で最適な解である。今後それを最
適解と呼ぶことにし、Woptと表記する。Consider a steady state case. In the steady state, each element of the correlation matrix R of Expression (2) and the cross-correlation vector p of Expression (3) is not a sum but is represented by a set average as in Expressions (10) and (11) ( i, j = 1,2). At this time, equation (8) is an optimal solution in the sense of minimizing the average power of the side lobe canceller output signal. In the future, it will be called the optimal solution, and will be denoted as W opt .
【0022】[0022]
【数10】 (Equation 10)
【0023】[0023]
【数11】 [Equation 11]
【0024】主アンテナの位置を位相中心として、補助
アンテナ1,2における所望信号の位相中心に対する位
相差ををα0,β0、干渉波iの位相中心に対する位相差
をα i,βiとする(i=1,2)。各補助アンテナの指
向特性は等方性で値は1とする。主アンテナ、補助アン
テナ1,2における受信信号x0(k),x1(k),x
2(k)は、それぞれ式(12)〜(14)のようにな
る。n0(k),n1(k),n2(k)は、それぞれ主
アンテナ,補助アンテナ1,補助アンテナ2における受
信機雑音で互いに無相関かつ入射波とも無相関であり、
G0,G1,G2は、それぞれ所望波,干渉波1,干渉波
2に対する主アンテナにおける複素利得である。Using the position of the main antenna as the phase center,
Position with respect to phase center of desired signal in antennas 1 and 2
Α0, Β0, The phase difference of the interference wave i with respect to the phase center
Is α i, Βi(I = 1, 2). Finger of each auxiliary antenna
The direction characteristic is isotropic and the value is 1. Main antenna, auxiliary antenna
Received signal x at teners 1 and 20(K), x1(K), x
Two(K) are as shown in equations (12) to (14), respectively.
You. n0(K), n1(K), nTwo(K) is the main
Antennas, auxiliary antennas 1 and 2
It is uncorrelated with each other and uncorrelated with the incident wave due to the transmitter noise,
G0, G1, GTwoAre the desired wave, the interference wave 1, and the interference wave, respectively.
2 is the complex gain at the main antenna.
【0025】[0025]
【数12】 (Equation 12)
【0026】[0026]
【数13】 (Equation 13)
【0027】[0027]
【数14】 [Equation 14]
【0028】相関行列Rの要素rijと相互相関ベクトル
pの要素piは、それぞれ式(15)〜(18)のよう
になる。Pd=E[|d(k)|2]、Psi=E[|si
(k)|2](i=1,2)、PN=E[|ni(k)
|2](i=0,1,2)で、各信号の電力である。The element r ij of the correlation matrix R and the element p i of the cross-correlation vector p are represented by the following equations (15) to (18). P d = E [| d (k) | 2 ], P si = E [| s i
(K) | 2 ] (i = 1, 2), P N = E [| n i (k)
| 2 ] (i = 0, 1, 2), which is the power of each signal.
【0029】[0029]
【数15】 (Equation 15)
【0030】[0030]
【数16】 (Equation 16)
【0031】[0031]
【数17】 [Equation 17]
【0032】[0032]
【数18】 (Equation 18)
【0033】式(7)の解、すなわち最適解Woptを式
(19)のように表す。式(19)右辺の各要素は式
(20)〜(22)で、さらに相関行列固有値λは式
(23)で表される。The solution of the equation (7), that is, the optimum solution W opt is represented as the equation (19). The elements on the right side of Expression (19) are expressed by Expressions (20) to (22), and the correlation matrix eigenvalue λ is expressed by Expression (23).
【0034】[0034]
【数19】 [Equation 19]
【0035】[0035]
【数20】 (Equation 20)
【0036】[0036]
【数21】 (Equation 21)
【0037】[0037]
【数22】 (Equation 22)
【0038】[0038]
【数23】 (Equation 23)
【0039】(β2−α2)と(β1−α1)の差が2πの
整数倍のとき、式(22)の第4項は0になり、第2項
と第3項は式(24)のようにまとめられる。式(2
3)の根号内も同様にまとめられる。つまり、入射角が
異なるにも関わらず、電力Ps1+Ps2の1つの干渉波が
入射するのと同等となる。When the difference between (β 2 -α 2 ) and (β 1 -α 1 ) is an integral multiple of 2π, the fourth term of equation (22) is 0, and the second and third terms are (24). Equation (2
The same applies to 3). In other words, this is equivalent to the incidence of one interference wave of power P s1 + P s2 , despite the different incident angles.
【0040】[0040]
【数24】 (Equation 24)
【0041】干渉波が1波の場合、Ps2=0とおいて式
(20)〜(22)は簡単にできる。それを式(19)
に代入して干渉波が1波の場合の最適解(式(25))
を得る。式(25)の分母が相関行列Rの行列式であ
る。さらに固有値λは式(26)のようになる。In the case where the number of interference waves is one, equations (20) to (22) can be simplified with P s2 = 0. Equation (19)
Optimal solution in the case of one interference wave by substituting into (Expression (25))
Get. The denominator of Expression (25) is the determinant of the correlation matrix R. Further, the eigenvalue λ is as shown in Expression (26).
【0042】[0042]
【数25】 (Equation 25)
【0043】[0043]
【数26】 (Equation 26)
【0044】補助アンテナ1素子あたりの受信所望信号
対雑音電力比SNRin、受信干渉信号対雑音電力比IN
R1inを式(27)のように定義する。以降、それぞれ
入力SNR,入力INRと呼ぶことがある。最適解であ
る式(25)は、式(28)のように入力SNR,入力
INRで表現できる。さらに、式(28)は、INR
1in>>SNRinおよびINR1in>>1という条件下で
は式(29)のように近似できる。Desired reception signal per auxiliary antenna element
Noise power ratio SNRin,Received interference signal to noise power ratio IN
R1inIs defined as in equation (27). Hereafter, each
They may be called input SNR and input INR. The optimal solution
Equation (25) is expressed by the input SNR and input
It can be represented by INR. Further, equation (28) is obtained by calculating INR
1in>>> SNRinAnd INR1in>> Under the condition of 1
Can be approximated as in equation (29).
【0045】[0045]
【数27】 [Equation 27]
【0046】[0046]
【数28】 [Equation 28]
【0047】[0047]
【数29】 (Equation 29)
【0048】補助アンテナ数が2、干渉波数が1の場合
のサイドローブキャンセラ出力信号z(k)は式(3
0)のようになる。従って、干渉波が1波の場合のサイ
ドローブキャンセラ出力信号電力Poutは式(31)の
ようになる。式(31)の一番下の式で、第1項は所望
信号電力、第2項は干渉信号電力、第3項は雑音電力で
ある。荷重係数ノルム2乗値が大きくなると第3項の雑
音電力が大きくなる。式(25)の最適解はサイドロー
ブキャンセラ出力信号電力Poutを最小にするものであ
る。When the number of auxiliary antennas is 2 and the number of interference waves is 1, the side lobe canceller output signal z (k) is given by the following equation (3).
0). Therefore, the side lobe canceller output signal power P out in the case of one interference wave is as shown in Expression (31). In the lowermost expression of Expression (31), the first term is the desired signal power, the second term is the interference signal power, and the third term is the noise power. The noise power of the third term increases as the weight coefficient norm squared value increases. The optimal solution of the equation (25) is to minimize the side lobe canceller output signal power P out .
【0049】[0049]
【数30】 [Equation 30]
【0050】[0050]
【数31】 (Equation 31)
【0051】式(31)の一番下の式の第1項を主アン
テナ受信所望信号電力|G0|2Pdで割ると、サイドロ
ーブキャンセラ出力において所望信号電力がどのくらい
保たれているかを表す量となる。式(32)で定義され
るその量を所望信号電力増加比DRと呼ぶことにする。
この量が小さいほどサイドローブキャンセラにおいて所
望信号は抑圧されるということである。INR1in>>
SNRinおよびINR1 in>>1という条件下で最適解
(式(29))に対する所望信号電力増加比DR optを
求めると式(33)のようになる。The first term of the lowermost equation of the equation (31) is replaced by
Tenor reception desired signal power | G0|TwoPdDivide by
What is the desired signal power at the output of the
It is a quantity indicating whether it is kept. Defined by equation (32)
This amount is called a desired signal power increase ratio DR.
The smaller this amount is, the more
The wanted signal is suppressed. INR1in>>
SNRinAnd INR1 in>> Optimal solution under the condition of 1
(Equation (29)) desired signal power increase ratio DR optTo
When it is found, it becomes like Expression (33).
【0052】[0052]
【数32】 (Equation 32)
【0053】[0053]
【数33】 [Equation 33]
【0054】例えば、G0=100、G1=1、SNRin
=0dBとすると、式(33)絶対値内の第3項はせい
ぜい0.01程度である。第2項は所望波や干渉波の入
射角にもよるが、0から2/3程度の値である。第2項
が0に近くなるのはα0−α1−β0+β1が2πの整数倍
に近い場合であり、それ以外はDRoptは1より小さく
なる。入力SNRであるSNRinがかなり小さくならな
い限り、ほとんどの場合についてSLC出力における所
望信号電力は低下する(所望信号が抑圧される)ことに
なる。For example, G 0 = 100, G 1 = 1, SNR in
Assuming that = 0 dB, the third term in the absolute value of Expression (33) is at most about 0.01. The second term is a value of about 0 to 2/3, depending on the incident angle of the desired wave or the interference wave. The second term is close to 0 when α 0 −α 1 −β 0 + β 1 is close to an integral multiple of 2π, and otherwise DR opt is smaller than 1. Unless the input SNR, SNR in, is significantly reduced, in most cases the desired signal power at the SLC output will decrease (the desired signal will be suppressed).
【0055】一方、式(31)の一番下の式の第3項を
主アンテナ受信機雑音電力PNで割ると雑音電力がどの
くらい増えたかを表す量となる。式(34)で定義され
るその量を雑音電力増加比NRと呼ぶことにする。On the other hand, a quantity representing whether increased how much the three term divided by the main antenna receiver noise power P N the noise power of the formula at the bottom of the formula (31). The amount defined by the equation (34) is referred to as a noise power increase ratio NR.
【0056】[0056]
【数34】 (Equation 34)
【0057】INR1in>>SNRinおよびINR1in>
>1という条件下で最適解(式(29))に対する雑音
電力増加比NRoptを計算する。結果は複雑になるた
め、最適解に対する雑音電力増加比NRoptを式(3
5)のようにおく。式(35)の第2項のNw (i)(i=
1,2,・・・,8),Dwは、それぞれ式(36)〜
(44)のようになる。INR 1in >> SNR in and INR 1in >>
The noise power increase ratio NR opt for the optimal solution (Equation (29)) is calculated under the condition of> 1. Since the result becomes complicated, the noise power increase ratio NR opt for the optimal solution is calculated by the equation (3)
Place as in 5). N w (i ) of the second term of equation (35) (i =
1, 2,..., 8) and D w are expressed by equations (36) to
It becomes like (44).
【0058】[0058]
【数35】 (Equation 35)
【0059】[0059]
【数36】 [Equation 36]
【0060】[0060]
【数37】 (37)
【0061】[0061]
【数38】 (38)
【0062】[0062]
【数39】 [Equation 39]
【0063】[0063]
【数40】 (Equation 40)
【0064】[0064]
【数41】 [Equation 41]
【0065】[0065]
【数42】 (Equation 42)
【0066】[0066]
【数43】 [Equation 43]
【0067】[0067]
【数44】 [Equation 44]
【0068】最適解に対する雑音電力増加比NRoptの
値を見積もってみる。G0=100、G1=1、SNRin
=0dB、INR1in≧20dBとする。このとき、式
(44)のDwの値は最小で4、最大で36程度であ
る。式(36)〜(43)のNw (i )(i=1,2,…,
8)のオーダはおよそ以下の通りである。Estimate the value of the noise power increase ratio NR opt for the optimal solution. G 0 = 100, G 1 = 1, SNR in
= 0 dB and INR 1in ≧ 20 dB. At this time, the value of D w in Expression (44) is about 4 at the minimum and about 36 at the maximum. N w (i ) (i = 1, 2,...,
The order of 8) is approximately as follows.
【0069】[0069]
【表1】 [Table 1]
【0070】このように見ると、一番支配的な項はG0
の2乗の項があり、しかもINR1inで割られることの
ないNw (2)である。これをDwの最大値約36で割って
も雑音電力増加比NRoptの値は2000を越えてしま
う。入力SNRであるSNRinがG0の逆数より小さく
なるか、α0−α1−β0+β1が2πの整数倍にならない
と雑音電力増加比NRoptの値は1に近づかない。しか
し式(37)の中かっこ内が非常に小さい値になる角度
は限られている。In this way, the most dominant term is G 0
N w (2) which has a squared term and is not divided by INR 1in . Even if this is divided by the maximum value of D w about 36, the value of the noise power increase ratio NR opt exceeds 2000. The value of the noise power increase ratio NR opt does not approach 1 unless the input SNR SNR in becomes smaller than the reciprocal of G 0 or α 0 −α 1 −β 0 + β 1 does not become an integral multiple of 2π. However, the angle at which the value inside the curly braces of Expression (37) becomes very small is limited.
【0071】以上説明したように、最適解で与えられる
荷重に対しては、補助アンテナ受信所望信号対雑音電力
比がかなり小さくならないと所望信号は抑圧されてしま
うし、また、サイドローブキャンセラ出力における雑音
電力は大きく増加してしまう。つまり、信号対雑音電力
比が大きく低下してしまう。これは補助アンテナに指向
性を持たせて干渉波方向に補助アンテナのビームを向け
た場合も同様である。As described above, with respect to the load given by the optimal solution, the desired signal is suppressed unless the ratio of the desired signal to the noise of the auxiliary antenna is considerably reduced, and the output of the side lobe canceller is reduced. The noise power increases greatly. That is, the signal-to-noise power ratio is greatly reduced. The same applies to the case where the auxiliary antenna has directivity and the beam of the auxiliary antenna is directed in the interference wave direction.
【0072】[0072]
【発明が解決しようとする課題】本発明は上述の問題点
を解決するためになされたもので、干渉波数が補助アン
テナ数より少ない場合、あるいは異なる角度から入射す
る干渉波数が補助アンテナ数と一致しても干渉波数が補
助アンテナ数より少ないとみなされる場合でも、荷重係
数ベクトルWのノルム||W||が大きくならず、干渉波を
十分抑圧し、かつ、サイドローブキャンセラ出力におけ
る信号対雑音電力比が低下することのないような干渉波
抑圧装置を提供することを目的とする。SUMMARY OF THE INVENTION The present invention has been made in order to solve the above-mentioned problems. When the number of interference waves is smaller than the number of auxiliary antennas, or when the number of interference waves incident from different angles is equal to the number of auxiliary antennas. Even when the number of interference waves is considered to be smaller than the number of auxiliary antennas, the norm || W || of the weighting factor vector W does not increase, the interference waves are sufficiently suppressed, and the signal-to-noise at the sidelobe canceller output is reduced. It is an object of the present invention to provide an interference wave suppression device that does not reduce the power ratio.
【0073】[0073]
【課題を解決するための手段】この発明に係る干渉波抑
圧装置は、所望の入射信号の入射方向にビームを向けた
主アンテナおよび主アンテナに近接して配置された複数
の補助アンテナヘの入射信号と同時に入射する干渉波信
号を抑圧する干渉波抑圧装置において、主アンテナおよ
び補助アンテナの各々に設けられ、主アンテナおよび補
助アンテナ受信信号を入力してディジタル同相・直交信
号に変換するディジタル同相・直交信号変換手段と、複
数の補助アンテナの受信信号に対応するディジタル同相
・直交信号に対して荷重係数により重み付けする荷重係
数重み付け手段と、重み付けされた複数のディジタル同
相・直交信号を加算する加算手段と、主アンテナの受信
信号に対応するディジタル同相・直交信号から加算され
たディジタル同相・直交信号を減算する減算手段と、補
助アンテナの受信信号に対応する複数のディジタル同相
・直交信号の相関行列を計算する相関行列計算手段と、
主アンテナの受信信号に対応するディジタル同相・直交
信号と補助アンテナの受信信号に対応する複数のディジ
タル同相・直交信号との間の相互相関ベクトルを計算す
る相互相関ベクトル計算手段と、相関行列の固有値を計
算する相関行列固有値計算手段と、固有値および相関行
列および相互相関ベクトルに基づいて荷重係数を計算
し、荷重係数を荷重係数重み付け手段に出力する荷重計
算手段とを備えている。According to the present invention, there is provided an interference wave suppressing apparatus comprising: a main antenna for directing a beam in a direction of incidence of a desired incident signal; and an incident signal to a plurality of auxiliary antennas arranged in proximity to the main antenna. In an interference wave suppression device for simultaneously suppressing an incident interference signal, a digital in-phase / quadrature is provided for each of a main antenna and an auxiliary antenna, and receives the main antenna and the auxiliary antenna and converts them into digital in-phase / quadrature signals. Signal conversion means, weighting means for weighting digital in-phase / quadrature signals corresponding to received signals of a plurality of auxiliary antennas with weight coefficients, and addition means for adding a plurality of weighted digital in-phase / quadrature signals. , The digital in-phase added from the digital in-phase and quadrature signals corresponding to the received signal of the main antenna Subtracting means for subtracting the quadrature signal, and correlation matrix calculating means for calculating a correlation matrix of a plurality of digital phase and quadrature signals corresponding to the received signal of the auxiliary antenna,
A cross-correlation vector calculation means for calculating a cross-correlation vector between a digital in-phase / quadrature signal corresponding to the received signal of the main antenna and a plurality of digital in-phase / quadrature signals corresponding to the received signal of the auxiliary antenna; And a load calculating means for calculating a weighting factor based on the eigenvalue, the correlation matrix and the cross-correlation vector, and outputting the weighting factor to the weighting factor weighting means.
【0074】また、主アンテナは、複数の素子アンテナ
から構成されるアレーアンテナである。The main antenna is an array antenna composed of a plurality of element antennas.
【0075】また、補助アンテナは、アレーアンテナを
構成する素子アンテナの一部の素子アンテナである。The auxiliary antenna is a part of the element antenna constituting the array antenna.
【0076】また、補助アンテナは、指向性アンテナで
あり、干渉波の到来方向を計測する干渉波到来方向計測
手段と、干渉波到来方向に補助アンテナの主ビーム方向
を向けるアンテナ指向方向調整手段とをさらに備えてい
る。The auxiliary antenna is a directional antenna, and includes an interference wave arrival direction measuring means for measuring the arrival direction of the interference wave, and an antenna directivity adjustment means for directing the main beam direction of the auxiliary antenna to the interference wave arrival direction. Is further provided.
【0077】また、補助アンテナは、複数の素子アンテ
ナから構成されたアレーアンテナであり、干渉波の到来
方向を計測する干渉波到来方向計測手段と、干渉波到来
方向に補助アンテナの主ビーム方向を向けるように補助
アンテナを構成する各素子アンテナに対して重み付けし
て合成する重み付け信号合成手段とをさらに備えてい
る。The auxiliary antenna is an array antenna composed of a plurality of element antennas, and includes an interference wave arrival direction measuring means for measuring an arrival direction of the interference wave, and a main beam direction of the auxiliary antenna for the interference wave arrival direction. Weighted signal combining means for weighting and combining each element antenna constituting the auxiliary antenna so as to direct the antenna is further provided.
【0078】さらに、補助アンテナは、主アンテナを構
成する素子アンテナの少なくとも一部で構成され、干渉
波の到来方向を計測する干渉波到来方向計測手段と、干
渉波到来方向に補助アンテナの主ビーム方向を向けるよ
うに補助アンテナを構成する各素子アンテナに対して重
み付けして合成する重み付け信号合成手段とをさらに備
えている。Further, the auxiliary antenna is constituted by at least a part of an element antenna constituting the main antenna, and includes an interference wave arrival direction measuring means for measuring an arrival direction of the interference wave, and a main beam of the auxiliary antenna in the interference wave arrival direction. Weighted signal combining means for weighting and combining each element antenna constituting the auxiliary antenna so as to direct the direction is further provided.
【0079】[0079]
【発明の実施の形態】実施の形態1.補助アンテナが2
素子の場合について説明する。最大2干渉波を抑圧でき
る。図1は本発明にかかわる干渉波抑圧装置の構成図で
ある。DESCRIPTION OF THE PREFERRED EMBODIMENTS Embodiment 1 2 auxiliary antennas
The case of an element will be described. Up to two interference waves can be suppressed. FIG. 1 is a configuration diagram of an interference wave suppressing apparatus according to the present invention.
【0080】図1において、100は所望波方向にビー
ムを向けた主アンテナ、101a,101bは指向性が
等方性に近い補助アンテナ、1a,1b,1cは受信
機、2a,2b,2cはA/D変換器である。受信機1
a〜1cおよびA/D変換器2a〜2cは、受信信号を
ディジタル同相・直交信号に変換するディジタル同相・
直交信号変換手段を構成している。In FIG. 1, 100 is a main antenna for directing a beam in a desired wave direction, 101a and 101b are auxiliary antennas whose directivity is close to isotropic, 1a, 1b and 1c are receivers, 2a, 2b and 2c are receivers. An A / D converter. Receiver 1
a to 1c and A / D converters 2a to 2c are digital in-phase converters for converting received signals into digital in-phase and quadrature signals.
It constitutes orthogonal signal conversion means.
【0081】3a,3bは、補助アンテナ受信信号に対
応する複数の上記ディジタル同相・直交信号を荷重係数
により重み付けする荷重係数重み付け手段である荷重係
数器である。4は重み付けされた信号を加算する加算手
段である加算器である。5は主アンテナ受信信号に対応
するディジタル同相・直交信号から上記重み付けされた
信号を減算する減算手段である減算器である。Reference numerals 3a and 3b denote weighting factors which are weighting means for weighting a plurality of the digital in-phase / quadrature signals corresponding to the auxiliary antenna reception signal by weighting factors. Reference numeral 4 denotes an adder which is an adding means for adding the weighted signals. Reference numeral 5 denotes a subtractor as subtraction means for subtracting the weighted signal from the digital in-phase / quadrature signal corresponding to the main antenna reception signal.
【0082】10aは、補助アンテナ受信信号に対応す
る複数の上記ディジタル同相・直交信号の相関行列を計
算する相関行列計算手段であると共に、主アンテナ受信
信号に対応するディジタル同相・直交信号と補助アンテ
ナ受信信号に対応する複数のディジタル同相・直交信号
との間の相互相関ベクトルを計算する相互相関ベクトル
計算手段である相関行列・相互相関ベクトル計算手段で
ある。相関行列・相互相関ベクトル計算手段10aで
は、式(2)の相関行列Rと式(3)の相互相関ベクト
ルpを計算する。Reference numeral 10a denotes a correlation matrix calculating means for calculating a correlation matrix of a plurality of digital in-phase / quadrature signals corresponding to the auxiliary antenna reception signal, and a digital in-phase / quadrature signal corresponding to the main antenna reception signal and the auxiliary antenna. Correlation matrix / cross-correlation vector calculation means, which is a cross-correlation vector calculation means for calculating a cross-correlation vector between a plurality of digital in-phase / quadrature signals corresponding to a received signal. The correlation matrix / cross-correlation vector calculation means 10a calculates the correlation matrix R of equation (2) and the cross-correlation vector p of equation (3).
【0083】11aは相関行列の固有値を計算する相関
行列固有値計算手段である。13aは固有値に応じて相
関行列と相関関数ベクトルの各要素から荷重係数を計算
する荷重計算手段である。Reference numeral 11a denotes a correlation matrix eigenvalue calculating means for calculating the eigenvalue of the correlation matrix. Reference numeral 13a denotes a load calculating unit that calculates a load coefficient from each element of the correlation matrix and the correlation function vector according to the eigenvalue.
【0084】相関行列固有値計算手段11aでは、相関
行列Rの固有値を計算する。相関行列Rの固有値を調べ
ることで、実質的な干渉波の数を知ることができる。
「実質的」という意味は、上述のように干渉波が2波入
射しても信号処理上は1波入射とみなされる場合がある
からである。相関行列Rは2×2の行列なので2つの固
有値は解析的に計算できて、式(15)を考慮すると式
(45)のようになる。The correlation matrix eigenvalue calculating means 11a calculates the eigenvalue of the correlation matrix R. By examining the eigenvalues of the correlation matrix R, it is possible to know the actual number of interference waves.
The meaning of “substantially” is because even if two interference waves are incident as described above, it may be regarded as one wave incident in signal processing. Since the correlation matrix R is a 2 × 2 matrix, two eigenvalues can be analytically calculated, and as shown in Expression (45) when Expression (15) is considered.
【0085】[0085]
【数45】 [Equation 45]
【0086】定常状態では相関行列Rの固有値は干渉波
が2波の場合は式(23)、1波の場合は式(26)の
ようになる。固有値λのオーダは、干渉波がない場合は
小さい方がPN、大きい方が2Pd+PN程度となる。干
渉波が1波の場合は大きい方が2Ps1、小さい方がPN
から2Pd+PN程度となる。干渉波が2波の場合は、大
きい方が2(Ps1+Ps2)、小さい方は入射角に依存し
て、最小でPNから2Pd+PN程度、最大でPs1とPs2
の小さい方の2倍程度となる。干渉波が2波の場合で小
さい方の固有値がPNから2Pd+PN程度となるのは、
実質的に干渉波が1波と見なされるような角度から2つ
の干渉波が入射する場合である。従って、小さい方の固
有値の大きさから干渉波が1波(干渉波が2波入射して
も実質的に1波入射とみなされる場合も含む。以降この
注釈を付けない)か2波かを判断することができる。有
限サンプル数で相関行列Rの要素を計算することを考慮
すると、雑音電力の2倍+想定する受信信号電力最大値
の2倍をしきい値として、2つの固有値がともにしきい
値以下なら干渉波数は0、2つの固有値のうち一方がし
きい値を越えて小さい方の固有値がこのしきい値以下な
ら干渉波数は1、2つの固有値ともにしきい値を越えれ
ば干渉波数は2と判定する。In the steady state, the eigenvalue of the correlation matrix R is as shown in Equation (23) when the number of interference waves is two, and as Equation (26) when the number of interference waves is one. Order of the eigenvalues λ is better when there is no interference wave has a small P N, the larger becomes about 2P d + P N. When the interference wave is one wave, the larger one is 2P s1 and the smaller one is P N
From this, it becomes about 2P d + P N. In the case of two interference waves, the larger one is 2 (P s1 + P s2 ), the smaller one depends on the angle of incidence and is at least PN to 2P d + P N at the minimum, and P s1 and P s2 at the maximum.
Is about twice as small as The reason why the smaller eigenvalue is about 2P d + P N from P N when there are two interference waves is
This is a case where two interference waves enter from an angle at which the interference wave is regarded as substantially one wave. Therefore, from the magnitude of the smaller eigenvalue, it is determined whether the interference wave is one wave (including the case where two interference waves are incident, but it is regarded as substantially one wave incident. You can judge. Considering that the elements of the correlation matrix R are calculated with a finite number of samples, if the threshold value is twice the noise power + twice the assumed maximum value of the received signal power, if both eigenvalues are equal to or less than the threshold value, interference occurs. The wave number is 0, if one of the two eigenvalues exceeds the threshold and the smaller eigenvalue is less than or equal to this threshold, the number of interference waves is determined to be 1. If both of the two eigenvalues exceed the threshold, the number of interference waves is determined to be 2. .
【0087】相関行列固有値計算手段11aでは相関行
列Rの固有値を式(45)により計算し、判定手段12
aでは計算された固有値に基づき、上述のように干渉波
数が1か2か判定する。荷重計算手段13aでは、相関
行列Rと相互相関ベクトルpから荷重係数を計算し、そ
の結果を荷重係数器3a,3bに出力するが、判定手段
12aでの判定結果に応じて荷重計算法が異なる。以
下、それについて説明する。The correlation matrix eigenvalue calculation means 11a calculates the eigenvalue of the correlation matrix R according to equation (45),
In a, it is determined whether the number of interference waves is 1 or 2 based on the calculated eigenvalue as described above. The load calculating unit 13a calculates a load coefficient from the correlation matrix R and the cross-correlation vector p, and outputs the result to the load coefficient units 3a and 3b. The load calculating method differs depending on the result of the determination by the determining unit 12a. . Hereinafter, this will be described.
【0088】干渉波数が2と判定された場合は式(9)
により荷重を計算し、その結果を荷重係数器3a,3b
に出力する。干渉波数が0と判定された場合は、W1=
W2=0として荷重係数器3a,3bに出力する。If the number of interference waves is determined to be 2, the equation (9)
And calculate the load by using the load coefficient units 3a and 3b.
Output to If the number of interference waves is determined to be 0, W 1 =
W 2 = 0 is output to the load coefficient units 3a and 3b.
【0089】干渉波数が1と判定した場合は、以下の3
つの荷重計算方法のうちのどれかにする。今まで行って
きた数値計算例によれば、干渉信号+雑音の抑圧性能に
最も優れるのは3つめの荷重計算方法3であったのでそ
れを推奨する。If the number of interference waves is determined to be 1, the following 3
One of the two load calculation methods. According to the numerical calculation examples performed so far, the third load calculation method 3 is the most excellent in the suppression performance of the interference signal + noise, and therefore it is recommended.
【0090】(荷重計算方法1)強制的にwi=0とす
る(i=1あるいは2)。これは補助アンテナ101a
あるいは101bからの信号をカットすることに等し
い。このとき荷重係数値は式(46)もしくは式(4
7)である。サイドローブキャンセラ出力における雑音
電力をなるべく小さくする観点からは、荷重係数ベクト
ルのノルム||W||の小さい方がよい。(Load Calculation Method 1) Forcibly set w i = 0 (i = 1 or 2). This is the auxiliary antenna 101a
Alternatively, it is equivalent to cutting the signal from 101b. At this time, the load coefficient value is calculated by the equation (46) or the equation (4).
7). From the viewpoint of minimizing the noise power in the sidelobe canceller output, the smaller the norm || W || of the weighting coefficient vector, the better.
【0091】[0091]
【数46】 [Equation 46]
【0092】[0092]
【数47】 [Equation 47]
【0093】(荷重計算方法2)式(2)の相関行列R
は干渉波数が1のとき、干渉信号電力が補助アンテナ受
信機雑音電力に比べて十分に大きいと特異に近い。それ
は、上述のように大きい方の固有値が干渉信号電力の2
倍で小さい方の固有値が受信機雑音電力と所望信号電力
の2倍の和程度で、両者の大きさの違いが大きいからで
ある。そこで、いっそのこと相関行列Rは特異(ランク
が1)であると仮定して、式(7)の最小二乗最小ノル
ム解を荷重係数とする。このとき荷重係数ノルムは小さ
く抑えられるため、サイドローブキャンセラ出力におけ
る雑音電力増加は最小限となる。また、荷重計算方法3
を説明した後で示すが、この最小二乗最小ノルム解は、
干渉波電力が補助アンテナ受信機雑音電力に比べて十分
に大きいと所望信号の抑圧も抑えられる。(Load calculation method 2) Correlation matrix R of equation (2)
Is near singular when the interference signal power is sufficiently larger than the auxiliary antenna receiver noise power when the number of interference waves is 1. That is, as described above, the larger eigenvalue is 2 of the interference signal power.
This is because the smaller and larger eigenvalue is about twice the sum of the receiver noise power and the desired signal power, and the difference between the two is large. Therefore, it is further assumed that the correlation matrix R is singular (rank 1), and the least-squares least-norm solution of the equation (7) is used as the weighting factor. At this time, the load coefficient norm is kept small, so that an increase in noise power at the output of the side lobe canceller is minimized. In addition, load calculation method 3
After explaining, this least-squares least-norm solution is
If the interference wave power is sufficiently larger than the auxiliary antenna receiver noise power, suppression of the desired signal can be suppressed.
【0094】式(7)の最小二乗最小ノルム解は補助ア
ンテナ数が多い場合は数値的に特異値分解を行って求め
なければならないが、補助アンテナ数が2の場合には特
異値分解を行う必要はない。以下、これについて説明す
る。相関行列Rが特異であるなら解は式(48)を満足
する。ただしその解は無数に存在する。そのなかで荷重
係数ベクトルWのノルム||W||を最小にするものが図1
のサイドローブキャンセラの荷重として望ましい。なぜ
なら、||W||2が小さいほどサイドローブキャンセラ出
力における雑音電力が小さくなるからである。When the number of auxiliary antennas is large, the least-squares least-norm solution of equation (7) must be numerically obtained by performing singular value decomposition, but when the number of auxiliary antennas is 2, singular value decomposition is performed. No need. Hereinafter, this will be described. If the correlation matrix R is singular, the solution satisfies equation (48). However, there are countless solutions. Among them, the one that minimizes the norm || W || of the load coefficient vector W is shown in FIG.
It is desirable as the load of the side lobe canceller. This is because the smaller || W || 2 , the smaller the noise power at the sidelobe canceller output.
【0095】[0095]
【数48】 [Equation 48]
【0096】最小ノルム解は以下のように求められる。
記述を簡単にするために、a=r11,b=r12,c=p
1とおく。式(48)からw2を消去すると、解は式(4
9)のようになる。このとき、荷重係数ベクトルWのノ
ルムの2乗||W||2は式(50)のようになる。The minimum norm solution is obtained as follows.
To simplify the description, a = r 11 , b = r 12 , c = p
Set to 1 . Eliminating w 2 from equation (48), the solution becomes equation (4)
It becomes like 9). At this time, the square of the norm of the load coefficient vector W || W || 2 is as shown in Expression (50).
【0097】[0097]
【数49】 [Equation 49]
【0098】[0098]
【数50】 [Equation 50]
【0099】||W||2を最小にするw1は式(50)をw
1に関して偏微分して0とおくことにより得られる。式
(50)の偏微分は式(51)のようになる。ここで、
∂w 1/∂w1=1,∂w1 */∂w1=0であることを用
いた。これより荷重係数ベクトルのノルムが最小となる
解Wmin.norm1、すなわち荷重係数は、式(52)によ
り相関行列Rと相互相関ベクトルpの要素より解析的に
簡単に計算できる。相関行列Rは実際には特異でないた
め、このように求められた荷重係数は最適解ではない
が、それに近い干渉信号抑圧性能を有していることは多
くの数値例で確認した。|| W ||TwoW to minimize1Is the equation (50)
1Is obtained by partially differentiating with respect to. formula
The partial differential of (50) is as shown in equation (51). here,
∂w 1/ ∂w1= 1, ∂w1 */ ∂w1= 0
Was. This minimizes the norm of the load coefficient vector
Solution Wmin.norm1That is, the load coefficient is calculated by the equation (52).
From the elements of the correlation matrix R and the cross-correlation vector p
Easy to calculate. The correlation matrix R was not actually singular
Is not the optimal solution
However, they often have interference signal suppression performance close to that.
This was confirmed by a number of numerical examples.
【0100】[0100]
【数51】 (Equation 51)
【0101】[0101]
【数52】 (Equation 52)
【0102】(荷重計算方法3)荷重計算方法2と同じ
く、式(2)の相関行列Rが特異であるとして式(7)
の最小二乗最小ノルム解を荷重係数として適用するが、
解の導き方が異なる。特異値分解ではなく、相関行列R
の大きい方の固有値に対応する固有ベクトルを用いて最
小二乗最小ノルム解を導ける。従来技術文献2(内田、
鷹尾著:“ビームスペース・パーシャリアダプティブア
レイに関する検討” 電子情報通信学会技術報告 A・
P88−52,1988)では最小二乗最小ノルム解と
は全く記していないが、そこで得られている解は相関行
列が特異であるとして得られる最小二乗最小ノルム解と
等しい。補助アンテナ数が2、干渉波数が1の場合は以
下のようになる。結果として、式(2)の相関行列Rの
固有ベクトルを求める必要はない。(Load calculation method 3) Similarly to the load calculation method 2, it is assumed that the correlation matrix R of the equation (2) is singular and the equation (7)
Apply the least-squares least-norm solution of as the weighting factor,
The way to derive the solution is different. Rather than singular value decomposition, the correlation matrix R
The least-squares least-norm solution can be derived using the eigenvector corresponding to the larger eigenvalue of. Prior art document 2 (Uchida,
By Takao: "A study on beam space / perspective adaptive array" IEICE technical report
P88-52, 1988) does not describe any least-squares least-norm solution, but the solution obtained therefrom is equal to the least-square least-norm solution obtained assuming that the correlation matrix is singular. When the number of auxiliary antennas is 2 and the number of interference waves is 1, the following is performed. As a result, there is no need to find the eigenvectors of the correlation matrix R in equation (2).
【0103】固有値と固有ベクトルを用いると、最小二
乗最小ノルム解Wmin.norm2(式(52)と区別するた
めに添字に2を付けた)は式(53)のようになる。式
(2)の相関行列Rの大きい方の固有値λ1(式(4
5)の大きい方)に対応する固有ベクトルe=[e1,
e2]Tを用いると、式(53)の行列R+は式(54)
で定義される一種の擬似的な逆行列である。Using the eigenvalues and eigenvectors, the least-squares least-norm solution W min.norm2 (subscripted with 2 to distinguish it from equation (52)) is as shown in equation (53). The larger eigenvalue λ 1 of the correlation matrix R in Equation (2) (Equation (4)
The eigenvector e = [e 1 ,
e 2 ] T , the matrix R + in equation (53) is given by equation (54)
Is a kind of pseudo inverse matrix defined by
【0104】[0104]
【数53】 (Equation 53)
【0105】[0105]
【数54】 (Equation 54)
【0106】式(53)を式(4)の相関行列Rと式
(5)の相互相関ベクトルpの要素で表す。||e||2=
1とする。Re=λ1eであることより式(55)が成
り立つ。これより式(56)が得られる。さらに、特性
多項式が0に等しいことを利用すると式(57)が得ら
れる。式(57)と(55)より、大きい方の固有値に
対応する固有ベクトルは式(58)のようになる。φは
任意の実数である。Equation (53) is represented by the elements of the correlation matrix R of equation (4) and the cross-correlation vector p of equation (5). || e || 2 =
Let it be 1. Equation (55) holds because Re = λ 1 e. From this, equation (56) is obtained. Further, by utilizing that the characteristic polynomial is equal to 0, Expression (57) is obtained. From equations (57) and (55), the eigenvector corresponding to the larger eigenvalue is as shown in equation (58). φ is an arbitrary real number.
【0107】[0107]
【数55】 [Equation 55]
【0108】[0108]
【数56】 [Equation 56]
【0109】[0109]
【数57】 [Equation 57]
【0110】[0110]
【数58】 [Equation 58]
【0111】式(58)を式(53)に代入すると、式
(59)のように最小二乗最小ノルム解Wmin.norm2を
得ることができる。すなわち、結果的に固有ベクトルを
求めることなしに荷重係数ベクトルを計算できる。相関
行列Rは実際には特異でないため、このように求められ
た荷重係数は最適解ではないが、それに近い干渉信号抑
圧性能を有していることは多くの数値例で確認した。By substituting equation (58) into equation (53), a least squares minimum norm solution W min.norm2 can be obtained as in equation (59). That is, as a result, the load coefficient vector can be calculated without obtaining the eigenvector. Since the correlation matrix R is not actually singular, the weighting factor obtained in this manner is not an optimal solution, but it has been confirmed in many numerical examples that it has interference signal suppression performance close to it.
【0112】[0112]
【数59】 [Equation 59]
【0113】[0113]
【数60】 [Equation 60]
【0114】式(52)と式(59)の違いだが、相関
行列Rが完全に特異であれば両者は一致する。しかし、
現実には受信機雑音があるから相関行列Rは特異に近い
とはいえ、特異にはならない。そのため、式(52)と
式(59)は実際には等しくない。多くの数値計算例に
よれば、式(59)の方が干渉信号+雑音の抑圧性能が
高いようである。また、定常状態における式(59)の
干渉信号抑圧性能は、多くの数値計算例によれば最適解
(式(9))のそれに近いものであった。つまり、最適
解ではないにも関わらず干渉信号抑圧能力は高い。The difference between the equations (52) and (59) is that if the correlation matrix R is completely singular, they coincide. But,
In reality, since there is receiver noise, the correlation matrix R is not singular although it is close to singular. Therefore, equation (52) and equation (59) are not actually equal. According to many numerical calculation examples, it seems that Expression (59) has higher interference signal + noise suppression performance. Further, the interference signal suppression performance of Equation (59) in the steady state is close to that of the optimal solution (Equation (9)) according to many numerical calculation examples. In other words, the interference signal suppressing ability is high even though it is not the optimal solution.
【0115】次に、式(52)と式(59)の最小二乗
最小ノルム解が、所望信号の抑圧も起こさず、また、サ
イドローブキャンセラ出力における雑音電力の増大も最
適解に比べて小さいことを示す。まず、式(52)の最
小ノルム解Wmin.norm1について、式(32)の所望信
号電力増加比と式(34)の雑音電力増加比を計算す
る。式(52)を各入射信号の電力あるいは入力SN
R,入力INRで表現すると式(61)のようになる。
式(61)は、INR1in>>SNRinおよびINR1in
>>1という条件下では式(62)のように近似でき
る。Next, the least-squares least-norm solution of Equations (52) and (59) does not cause suppression of the desired signal, and the increase in noise power at the sidelobe canceller output is smaller than that of the optimal solution. Is shown. First, for the minimum norm solution W min.norm1 in Expression (52), the desired signal power increase ratio in Expression (32) and the noise power increase ratio in Expression (34) are calculated. Equation (52) is used to calculate the power or input SN of each incident signal.
When expressed by R and input INR, it is as shown in equation (61).
Equation (61) is given by: INR 1in >> SNR in and INR 1in
Under the condition of >> 1, it can be approximated as in equation (62).
【0116】[0116]
【数61】 [Equation 61]
【0117】[0117]
【数62】 (Equation 62)
【0118】INR1in>>SNRinおよびINR1in>
>1という条件下で、式(62)に対する所望信号電力
増加比DRmin.norm1と雑音電力増加比NRmin.norm1を
求める。これは式(62)を式(32)と式(34)に
代入して得られる。その結果はそれぞれ式(63)と式
(64)のようになる。INR 1in >> SNR in and INR 1in >>
Under the condition of> 1, the desired signal power increase ratio DR min.norm1 and the noise power increase ratio NR min.norm1 for Expression (62) are obtained. This is obtained by substituting equation (62) into equations (32) and (34). The results are as shown in equations (63) and (64), respectively.
【0119】[0119]
【数63】 [Equation 63]
【0120】[0120]
【数64】 [Equation 64]
【0121】式(63)の所望信号電力増加比DR
min.norm1と式(64)の雑音電力増加比NRmin.norm1
の値を見積もってみる。上述と同様、G0=100,G1
=1,SNRin=0dB,INR1in≧20dBとす
る。式(63)の絶対値内の第2項は0.01以下、第
3項も約0.01である。従って式(63)の所望信号
電力増加比DRmin.norm1はほぼ1である。つまり、所
望信号の抑圧はINR1in>>SNRinおよびINR1in
>>1、さらに主アンテナのサイドローブレベルが十分
に低い(G1/G0の値が十分小さい)限りほとんど生じ
ないことになる。The desired signal power increase ratio DR of the equation (63)
min.norm1 and the noise power increase ratio NR min.norm1 in equation (64)
Try to estimate the value of. As above, G 0 = 100, G 1
= 1, SNR in = 0 dB, INR 1in ≧ 20 dB. The second term in the absolute value of the equation (63) is 0.01 or less, and the third term is also about 0.01. Therefore, the desired signal power increase ratio DR min.norm1 in the equation (63) is almost 1. That is, the suppression of the desired signal is represented by INR 1in >> SNR in and INR 1in
>> 1 and hardly occurs as long as the side lobe level of the main antenna is sufficiently low (the value of G 1 / G 0 is sufficiently small).
【0122】式(64)の第1項は約1、第2項は約
0.5、第3項+第4項は1.5である。雑音電力増加
比NRmin.norm1の値はせいぜい3程度であり、最適解
のような大きな値とはならない。つまり、サイドローブ
キャンセラ出力における雑音電力は最適解ほど増加しな
い。以上まとめると、式(52)の荷重係数はサイドロ
ーブキャンセラ出力において所望信号電力の抑圧を起こ
さないし、雑音電力の増加も大きくない。従って信号対
雑音電力比は低下しない。In the equation (64), the first term is about 1, the second term is about 0.5, and the third and fourth terms are 1.5. The value of the noise power increase ratio NR min.norm1 is at most about 3 and does not become a large value like the optimal solution. That is, the noise power at the output of the sidelobe canceller does not increase as much as the optimal solution. In summary, the weighting factor in equation (52) does not suppress the desired signal power at the sidelobe canceller output, and the increase in noise power is not large. Therefore, the signal to noise power ratio does not decrease.
【0123】同様に、式(59)の最小ノルム解W
min.norm2について、式(32)の所望信号電力増加比
と式(34)の雑音電力増加比を計算する。式(59)
を各入射信号の電力あるいは入力SNR,入力INRで
表現する。式(65)のようになる。式(65)は、I
NR1in>>SNRinおよびINR1in>>1という条件
下では式(66)のように近似できる。Similarly, the minimum norm solution W of equation (59)
For min.norm2 , the desired signal power increase ratio of Expression (32) and the noise power increase ratio of Expression (34) are calculated. Equation (59)
Is represented by the power of each incident signal or the input SNR and input INR. Equation (65) is obtained. Equation (65) gives
Under the condition of NR 1in >> SNR in and INR 1in >> 1, approximation can be made as in equation (66).
【0124】[0124]
【数65】 [Equation 65]
【0125】[0125]
【数66】 [Equation 66]
【0126】INR1in>>SNRinおよびINR1in>
>1という条件下で、式(66)に対する所望信号電力
増加比DRmin.norm2と雑音電力増加比NRmin.norm2を
求める。これは式(66)を式(32)と式(34)に
代入して得られる。その結果はそれぞれ式(67)と式
(68)のようになる。INR 1in >> SNR in and INR 1in >>
Under the condition of> 1, a desired signal power increase ratio DR min.norm2 and a noise power increase ratio NR min.norm2 for Expression (66) are obtained. This is obtained by substituting equation (66) into equations (32) and (34). The results are as shown in equations (67) and (68), respectively.
【0127】[0127]
【数67】 [Equation 67]
【0128】[0128]
【数68】 [Equation 68]
【0129】式(67)の所望信号電力増加比DR
min.norm2と式(68)の雑音電力増加比NRmin.norm2
の値を見積もってみる。上述と同様、G0=100、G1
=1、SNRin=0dB、INR1in≧20dBとす
る。式(67)の絶対値内の第2項は0.01以下、第
3項も約0.01である。従って式(67)の所望信号
電力増加比DRmin.norm2はほぼ1である。つまり、所
望信号の抑圧はINR1in>>SNRinおよびINR1in
>>1、さらに主アンテナのサイドローブレベルが十分
に低い限りほとんど生じないことになる。The desired signal power increase ratio DR in equation (67)
min.norm2 and the noise power increase ratio NR min.norm2 in equation (68)
Try to estimate the value of. As above, G 0 = 100, G 1
= 1, SNR in = 0 dB, INR 1in ≧ 20 dB. The second term in the absolute value of the equation (67) is 0.01 or less, and the third term is also about 0.01. Therefore, the desired signal power increase ratio DR min.norm2 in equation (67) is approximately 1. That is, the suppression of the desired signal is represented by INR 1in >> SNR in and INR 1in
>> 1 and hardly any as long as the side lobe level of the main antenna is sufficiently low.
【0130】式(68)の第1項は約1、第2項は約
0.5、第3項+第4項は1.5である。雑音電力増加
比NRmin.norm2の値はせいぜい3程度であり、最適解
のような大きな値とはならない。つまり、サイドローブ
キャンセラ出力における雑音電力は最適解ほど増加しな
い。以上まとめると、式(52)の荷重係数はサイドロ
ーブキャンセラ出力において所望信号電力の抑圧を起こ
さないし、雑音電力の増加も大きくない。従って信号対
雑音電力比は低下しない。In the equation (68), the first term is about 1, the second term is about 0.5, and the third and fourth terms are 1.5. The value of the noise power increase ratio NR min.norm2 is at most about 3 and does not become a large value like the optimal solution. That is, the noise power at the output of the sidelobe canceller does not increase as much as the optimal solution. In summary, the weighting factor in equation (52) does not suppress the desired signal power at the sidelobe canceller output, and the increase in noise power is not large. Therefore, the signal to noise power ratio does not decrease.
【0131】以上示したように、実施の形態1によれ
ば、サイドローブキャンセラの補助アンテナ数が2の場
合において、干渉波数が1、もしくは干渉波数が2でも
実質的に1とみなされる場合でも、干渉波を十分抑圧
し、かつ、サイドローブキャンセラ出力における信号対
雑音電力比が低下することのないような荷重係数を極め
て簡単に解析的に計算できる干渉波抑圧装置を提供でき
る。As described above, according to Embodiment 1, when the number of auxiliary antennas of the side lobe canceller is 2, even if the number of interference waves is 1, or even if the number of interference waves is 2, it is regarded as substantially 1 In addition, it is possible to provide an interference wave suppressing apparatus which can sufficiently suppress an interference wave and which can extremely easily and analytically calculate a weighting factor such that a signal-to-noise power ratio in a sidelobe canceller output does not decrease.
【0132】実施の形態2.補助アンテナが3素子の場
合について説明する。最大3干渉波を抑圧できる。図2
は本発明にかかわる干渉波抑圧装置の構成図である。本
構成の図1と異なる点は、補助アンテナが1つ増えたこ
と、相関行列・相互相関ベクトル計算手段10bで計算
すべき行列とベクトルのサイズが大きくなっているこ
と、相関行列固有値計算手段11bでの固有値計算方法
が相関行列固有値計算手段11aと異なること、荷重計
算手段13bでの荷重計算方法が荷重計算手段13aと
異なる点にある。それ以外は同じであるので、図2にお
ける相関行列固有値計算方法と荷重計算方法についての
み説明する。Embodiment 2 The case where the auxiliary antenna has three elements will be described. Up to three interference waves can be suppressed. FIG.
1 is a configuration diagram of an interference wave suppression device according to the present invention. This configuration differs from FIG. 1 in that the number of auxiliary antennas is increased by one, the size of the matrix and vector to be calculated by the correlation matrix / cross-correlation vector calculation means 10b is large, and the correlation matrix eigenvalue calculation means 11b Is different from the correlation matrix eigenvalue calculation means 11a in that the load calculation means 13b is different from the load calculation means 13a. Since the rest is the same, only the correlation matrix eigenvalue calculation method and the load calculation method in FIG. 2 will be described.
【0133】相関行列・相互相関ベクトル計算手段10
bでは、主アンテナ受信信号x0(k)と補助アンテナ
受信信号xi(k)(i=1,2,3)から式(69)
の相関行列Rと式(70)の相互相関ベクトルpを計算
する。それらの各要素は式(4)と式(5)で計算でき
る。ただし、i,j=1,2,3とする。SMI法で解
くべき方程式は式(71)で、その荷重係数ベクトルW
は式(72)で定義される。Correlation matrix / cross-correlation vector calculation means 10
In b, from the main antenna received signal x 0 (k) and the auxiliary antenna received signal x i (k) (i = 1, 2, 3), the equation (69) is used.
And the cross-correlation vector p of equation (70) are calculated. Each of those elements can be calculated by equations (4) and (5). However, i, j = 1, 2, and 3. The equation to be solved by the SMI method is given by equation (71), and the weight coefficient vector W
Is defined by equation (72).
【0134】[0134]
【数69】 [Equation 69]
【0135】[0135]
【数70】 [Equation 70]
【0136】[0136]
【数71】 [Equation 71]
【0137】[0137]
【数72】 [Equation 72]
【0138】相関行列固有値計算手段11bでは相関行
列Rの固有値を式(73)〜式(78)により計算す
る。相関行列Rの特性多項式は固有値λの3次方程式で
あり、式(73)となる。式(73)の3次方程式の解
は公知の解の公式により式(74)〜式(77)の順に
定数を定義して計算すると固有値である式(78)を得
る。The correlation matrix eigenvalue calculation means 11b calculates the eigenvalues of the correlation matrix R according to equations (73) to (78). The characteristic polynomial of the correlation matrix R is a cubic equation of the eigenvalue λ, and is represented by Expression (73). The solution of the cubic equation of equation (73) is defined by a constant in the order of equations (74) to (77) according to a known solution formula, and calculation yields equation (78), which is an eigenvalue.
【0139】[0139]
【数73】 [Equation 73]
【0140】[0140]
【数74】 [Equation 74]
【0141】[0141]
【数75】 [Equation 75]
【0142】[0142]
【数76】 [Equation 76]
【0143】[0143]
【数77】 [Equation 77]
【0144】[0144]
【数78】 [Equation 78]
【0145】判定手段12bでは、計算された式(7
8)の固有値に基づき、相関行列Rが特異に近いかそう
でないか判定する。判定手段12bは3つの固有値の中
にしきい値以下のものがいくつあるかを判断する。補助
アンテナ数からそのような固有値の数を取り除いたもの
が実質的な干渉波の数である。しきい値は、多くの計算
例から無難と思われる雑音電力の2倍+想定される所望
信号最大電力の2倍とする。In the judgment means 12b, the calculated expression (7
Based on the eigenvalue of 8), it is determined whether the correlation matrix R is close to singular or not. The determining means 12b determines how many of the three eigenvalues are equal to or smaller than the threshold value. The actual number of interference waves is obtained by removing the number of such eigenvalues from the number of auxiliary antennas. The threshold value is set to twice the noise power considered safe from many calculation examples + twice the assumed desired signal maximum power.
【0146】荷重計算手段13bでは、相関行列Rと相
互相関ベクトルpから荷重係数を計算するが、判定手段
12bでの判定結果に応じて荷重計算手段13bにおけ
る荷重計算法が異なる。以下、それについて説明する。The load calculating means 13b calculates a load coefficient from the correlation matrix R and the cross-correlation vector p, but the load calculating method in the load calculating means 13b differs depending on the result of the judgment by the judging means 12b. Hereinafter, this will be described.
【0147】(1)3つの固有値がしきい値以上の場合 式(79)、式(80)により荷重を計算し、その結果
を荷重係数器3a,3b,3cに出力する。(1) When the three eigenvalues are equal to or larger than the threshold value: The load is calculated by the equations (79) and (80), and the result is output to the load coefficient units 3a, 3b and 3c.
【0148】[0148]
【数79】 [Expression 79]
【0149】[0149]
【数80】 [Equation 80]
【0150】(2)しきい値以下の固有値の数が2の場
合 干渉波数は1と判断される。相関行列Rは特異に近い。
Rのランクは1に近い。(A)と(B)の2つの荷重計
算方法があるが、荷重係数ベクトルのノルムをなるべく
小さくする観点から(B)を推奨する。(B)はRのラ
ンクが1であると仮定してノルムが最小となる解(荷
重)の計算である。(2) When the number of eigenvalues equal to or less than the threshold value is 2 The number of interference waves is determined to be 1. The correlation matrix R is close to singular.
The rank of R is close to one. There are two load calculation methods (A) and (B), but (B) is recommended from the viewpoint of minimizing the norm of the load coefficient vector. (B) is the calculation of the solution (load) that minimizes the norm, assuming that the rank of R is 1.
【0151】(荷重計算方法A)強制的にwi=wj=0
とする(i,j=1,2,3;i≠j)。これは3つの
補助アンテナ101a,101b,101cのうちの2
つの補助アンテナからの信号をカットすることに等し
い。このとき解(荷重)は式(81)〜式(83)の3
通りのうちのどれかとする。サイドローブキャンセラ出
力における雑音電力をなるべく小さくする観点からは、
荷重係数ベクトルのノルム||W||が小さい方がよい。(Load Calculation Method A) Forcibly w i = w j = 0
(I, j = 1, 2, 3; i ≠ j). This is two of the three auxiliary antennas 101a, 101b, 101c.
Equivalent to cutting the signal from two auxiliary antennas. At this time, the solution (load) is expressed by Equations (81) to (83).
Any of the streets. From the viewpoint of minimizing the noise power in the sidelobe canceller output,
The smaller the norm || W || of the load coefficient vector, the better.
【0152】[0152]
【数81】 [Equation 81]
【0153】[0153]
【数82】 (Equation 82)
【0154】[0154]
【数83】 [Equation 83]
【0155】(荷重計算方法B)式(69)の相関行列
Rのランクが1であると仮定したときにノルムが最小と
なる解を荷重とする。補助アンテナ数が3の場合も特異
値分解を行う必要はない。以下、これについて説明す
る。rank(R)=1なら、式(71)の解は式(8
4)を満足するw1,w2,w3で、これは無数に存在す
る。そのなかで荷重係数ベクトルWのノルム||W||を最
小にするものが図2のサイドローブキャンセラの荷重と
して最も望ましい。なぜなら、||W||2の値が小さいほ
どサイドローブキャンセラ出力における雑音電力が小さ
くなるからである。(Load Calculation Method B) Assuming that the rank of the correlation matrix R in equation (69) is 1, the solution with the minimum norm is taken as the load. Even when the number of auxiliary antennas is 3, it is not necessary to perform singular value decomposition. Hereinafter, this will be described. If rank (R) = 1, the solution of equation (71) is given by equation (8)
There are countless w 1 , w 2 , and w 3 that satisfy 4). Among them, the one that minimizes the norm || W || of the load coefficient vector W is most desirable as the load of the side lobe canceller in FIG. This is because the smaller the value of || W || 2, the smaller the noise power in the sidelobe canceller output.
【0156】[0156]
【数84】 [Equation 84]
【0157】そのような解は以下のように求められる。
記述を簡単にするために、a=r11,b=r12,c=r
13とおく。これらのa,b,cは実施の形態1とは定義
が異なる。式(84)からw3を消去すると解は式(8
5)となる。このとき、式(85)の荷重係数ベクトル
Wのノルムの2乗||W||2は式(86)となる。Such a solution is obtained as follows.
To simplify the description, a = r 11 , b = r 12 , c = r
Enter 13 . These a, b, and c have different definitions from the first embodiment. Eliminating w 3 from equation (84) gives a solution to equation (8)
5). At this time, the square || W || 2 of the norm of the load coefficient vector W in the equation (85) becomes the equation (86).
【0158】[0158]
【数85】 [Equation 85]
【0159】[0159]
【数86】 [Equation 86]
【0160】式(86)の荷重係数ベクトルWのノルム
の2乗||W||2を最小にするw1,w 2は、式(86)を
w1,w2に関して偏微分して0とおくことにより得られ
る。偏微分した結果を式(87)、式(88)に示す。The norm of the load coefficient vector W in equation (86)
Square of || W ||TwoW to minimize1, W TwoReplaces equation (86) with
w1, WTwoIs obtained by partially differentiating with respect to
You. Expressions (87) and (88) show the results of partial differentiation.
【0161】[0161]
【数87】 [Equation 87]
【0162】[0162]
【数88】 [Equation 88]
【0163】式(87)、(88)よりw1,w2が求め
られ、さらに式(84)よりw3が求められる。荷重係
数ベクトルのノルムが最小となる解(荷重係数ベクト
ル)は式(89)の通りとなる。W 1 and w 2 are obtained from equations (87) and (88), and w 3 is obtained from equation (84). The solution (load coefficient vector) that minimizes the norm of the load coefficient vector is as shown in equation (89).
【0164】[0164]
【数89】 [Equation 89]
【0165】(3)しきい値以下の固有値の数が1の場
合 相関行列Rは特異に近い。干渉波数は2と判断される。
Rのランクは2に近い。(C)と(D)の2つの荷重計
算方法が考えられる。(D)はRのランクが2であると
仮定してノルムが最小となる荷重の計算である。(3) When the number of eigenvalues equal to or smaller than the threshold value is 1 The correlation matrix R is almost singular. The number of interference waves is determined to be 2.
The rank of R is close to two. Two load calculation methods (C) and (D) can be considered. (D) is the calculation of the load that minimizes the norm, assuming that the rank of R is 2.
【0166】(荷重計算方法C)強制的にwi=0とす
る。i=1,2,3のどれかである。これは補助アンテ
ナ101a,101b,101cからのどれか1つの信
号をカットすることに等しい。このとき解くべき方程式
は2元1次連立方程式となる。式(9)と同様にして荷
重を解析的に得ることができる。(Load Calculation Method C) Forcibly set w i = 0. i = 1, 2, or 3. This is equivalent to cutting any one signal from the auxiliary antennas 101a, 101b, and 101c. The equation to be solved at this time is a two-dimensional linear simultaneous equation. The load can be obtained analytically in the same manner as in equation (9).
【0167】(荷重計算方法D)Rのランクが2である
と仮定したときにノルムが最小となる解を荷重とする。
特異値分解を行う必要はない。以下、これについて説明
する。まず、3元1次連立方程式(式(70))の中か
ら独立な2つの方程式を選び出す。それを式(90)の
ようにおく。ここで、[a,b,c]=[ri1,ri2,
ri3]、pA=pi(i=1,2のどれか)、[d,e,
f]=[rj1,rj2,rj3]、pB=pj (j=2,3
のどれか、ただし、j≠i)である。定数a,b,cは
荷重計算方法Bと意味が異なる。(Load Calculation Method D) When it is assumed that the rank of R is 2, the solution that minimizes the norm is defined as the load.
There is no need to perform singular value decomposition. Hereinafter, this will be described. First, two independent equations are selected from the ternary linear equations (Equation (70)). It is set as in equation (90). Here, [a, b, c] = [ ri1 , ri2 ,
r i3 ], p A = p i (i = 1 or 2), [d, e,
f] = [r j1 , r j2 , r j3 ], p B = p j (j = 2, 3
Where j ≠ i). The constants a, b, and c have different meanings from the load calculation method B.
【0168】[0168]
【数90】 [Equation 90]
【0169】式(90)からw2,w3をw1で表せる。
まず、式(90)の行列の2行目の式から式(91)を
得る。式(91)を式(90)の行列の1行目の式に代
入して式(92)を得る。さらに式(92)を式(9
1)に代入して式(93)を得る。From equation (90), w 2 and w 3 can be represented by w 1 .
First, Expression (91) is obtained from the expression in the second row of the matrix of Expression (90). The equation (91) is substituted into the equation of the first row of the matrix of the equation (90) to obtain the equation (92). Further, equation (92) is replaced with equation (9)
Substituting into 1) yields equation (93).
【0170】[0170]
【数91】 [Equation 91]
【0171】[0171]
【数92】 (Equation 92)
【0172】[0172]
【数93】 [Equation 93]
【0173】簡単のため、式(92)、式(93)をそ
れぞれ式(94)、式(95)のようにおく。解(荷重
係数)Wは1変数だけで式(96)のように表される。
そしてそのノルムの2乗は式(97)で与えられる。For simplicity, equations (92) and (93) are set as equations (94) and (95), respectively. The solution (load coefficient) W is represented by equation (96) using only one variable.
The square of the norm is given by equation (97).
【0174】[0174]
【数94】 [Equation 94]
【0175】[0175]
【数95】 [Equation 95]
【0176】[0176]
【数96】 [Equation 96]
【0177】[0177]
【数97】 (97)
【0178】ノルムが最小となる解は、式(97)の||
W||2を最小にするw1を求めることで得られる。式(9
7)の||W||2をw1に関して偏微分して0とおく(式
(98)).これから式(99)の通りw1を求めるこ
とができる。The solution with the minimum norm is given by ||
It is obtained by finding w 1 that minimizes W || 2 . Equation (9
|| W || 2 is partially differentiated with respect to w 1 to 0 (Equation (98)). From this, w 1 can be obtained as in equation (99).
【0179】[0179]
【数98】 [Equation 98]
【0180】[0180]
【数99】 [Equation 99]
【0181】式(99)を式(96)に代入すると、荷
重係数ベクトルは式(100)のようになる。ここで
α,β,γ,δは式(94)(95)の通りである。By substituting equation (99) into equation (96), the load coefficient vector becomes as shown in equation (100). Here, α, β, γ, and δ are as shown in equations (94) and (95).
【0182】[0182]
【数100】 [Equation 100]
【0183】(4)3つの固有値がすべてしきい値以下
の場合 干渉波は存在しないとみなし、W1=W2=W3=0とし
て、荷重係数器3a,3b,3cに出力する。[0183] regarded as (4) three eigenvalues do not exist when the interference waves below the threshold, all as W 1 = W 2 = W 3 = 0, and outputs the load coefficient units 3a, 3b, to 3c.
【0184】以上示したように、実施の形態2によれ
ば、サイドローブキャンセラの補助アンテナ数が3の場
合において、干渉波数が2以下、もしくは干渉波数が3
でも実質的に2以下とみなされる場合でも、干渉波を十
分抑圧し、かつ、サイドローブキャンセラ出力における
信号対雑音電力比が低下することのないような荷重係数
を解析的に計算できる干渉波抑圧装置を提供できる。As described above, according to Embodiment 2, when the number of auxiliary antennas of the side lobe canceller is 3, the number of interference waves is 2 or less, or the number of interference waves is 3
However, even when the signal is regarded as being substantially equal to or less than 2, the interference wave can be suppressed sufficiently, and the weighting factor can be analytically calculated so that the signal-to-noise power ratio at the side lobe canceller output does not decrease. Equipment can be provided.
【0185】実施の形態3.実施の形態1と2では主ア
ンテナの構造について特定しなかったが、複数の素子ア
ンテナからなり、所望波方向にビームを向けるように荷
重が調整されたアレーアンテナを用いることができる。
さらにこの場合、補助アンテナとしてアレーアンテナを
構成する任意の素子アンテナを用いることも可能であ
る。図3に補助アンテナとして用いる素子アンテナが2
の場合の構成図を示す。200はN個の素子アンテナ2
01a,201b,…,201nからなるアレーアンテ
ナである(素子数Nは想定する最大信号数+1以上)。Embodiment 3 Although the structure of the main antenna is not specified in the first and second embodiments, an array antenna composed of a plurality of element antennas and having a load adjusted to direct a beam in a desired wave direction can be used.
Further, in this case, any element antenna constituting the array antenna can be used as the auxiliary antenna. FIG. 3 shows two element antennas used as auxiliary antennas.
FIG. 200 is N element antennas 2
, 201n (the number of elements N is equal to or more than the assumed maximum number of signals + 1).
【0186】N個の素子アンテナのうちの適当な2素子
(図3では201cと201e)を補助アンテナとす
る。素子アンテナ201a,201b,…,201nに
は主ビームが所望波方向に向くように適切な荷重を掛け
ておく。荷重の掛け方は、アナログ処理に基づいて位相
器を用いても、あるいは各素子アンテナ受信信号をディ
ジタル同相・直交信号に変換して荷重を掛けるディジタ
ルビームフォーミングでもよい。前者ではビームフォー
ミングの後に信号をディジタル同相・直交信号に変換す
る。図3は後者の構成で、202a,202b,…,2
02nは受信機、203a,203b,…,203nは
A/D変換器、204a,204b,…,204nはビ
ームを所望波方向へ向けるための荷重係数器である。2
05は荷重係数器204a,204b,…,204n出
力信号を合成する合成器である。主アンテナとしてのア
レーアンテナ200で受信された干渉信号のレプリカを
生成するための信号はA/D変換器203a〜203n
出力のどれかからとる(図3では203cと203
e)。それ以下の構成、機能は実施の形態1と同様であ
るが、特長として、補助アンテナとしてアレーアンテナ
を構成する任意の素子アンテナを用いた場合、補助アン
テナを特に設ける必要はなくなるため、スペースファク
タに優れる。Appropriate two elements (201c and 201e in FIG. 3) of the N element antennas are used as auxiliary antennas. An appropriate load is applied to the element antennas 201a, 201b,..., 201n so that the main beam is directed to a desired wave direction. The weight may be applied by using a phase shifter based on analog processing, or by digital beam forming in which the received signal of each element antenna is converted into a digital in-phase / quadrature signal and the weight is applied. In the former, a signal is converted into a digital in-phase / quadrature signal after beamforming. FIG. 3 shows the latter configuration, in which 202a, 202b,.
.., 203n are A / D converters, and 204a, 204b,..., 204n are weighting factors for directing a beam in a desired wave direction. 2
Reference numeral 05 denotes a combiner for combining the output signals of the weighting coefficients 204a, 204b,..., 204n. Signals for generating replicas of the interference signals received by the array antenna 200 as the main antenna are A / D converters 203a to 203n.
From one of the outputs (203c and 203c in FIG. 3)
e). The configuration and functions below this are the same as those of the first embodiment. However, when an arbitrary element antenna constituting an array antenna is used as an auxiliary antenna, there is no need to provide an auxiliary antenna. Excellent.
【0187】実施の形態4.実施の形態1と2では補助
アンテナの指向特性は等方性としたが、指向性を持たせ
たアンテナを補助アンテナとして用いてもよい。ただ、
補助アンテナのビーム方向はおおまかに個々の干渉波の
方向に向ける必要があるため、干渉波到来方向計測手段
を必要とする。その計測結果に基づき、補助アンテナの
指向方向を決める。Embodiment 4 Although the directional characteristics of the auxiliary antenna are isotropic in the first and second embodiments, an antenna having directivity may be used as the auxiliary antenna. However,
Since the beam direction of the auxiliary antenna needs to be roughly directed to the direction of each interference wave, an interference wave arrival direction measuring means is required. The directional direction of the auxiliary antenna is determined based on the measurement result.
【0188】図4は干渉波到来方向計測手段を備えた補
助アンテナ数2の干渉波抑圧装置の構成図である。図1
と共通部分については説明を省く。20aは干渉波到来
方向計測手段、21は干渉波到来方向計測手段20aで
の計測結果に基づいて、補助アンテナの指向方向を干渉
波到来方向に向けるためのアンテナ指向方向調整手段で
ある。干渉波到来方向計測は、例えば主アンテナビーム
を全空間に対して走査し、電力の大きい到来波の到来方
向を検出したり、あるいは別に干渉波到来計測用のアン
テナを設けて行う。図4は前者の構成である。補助アン
テナのビーム幅が狭くなければ、精密な計測は必要な
い。アンテナ指向方向調整手段21では計測された到来
方向に2つの補助アンテナの主ビームを向けさせる。干
渉波が2波検出されれば、2つの補助アンテナ指向方向
をそれぞれの干渉波到来方向に向ける。干渉波が1波検
出されれば、2つの補助アンテナ両方の指向方向をその
到来方向に向ける。補助アンテナの指向方向調整方法は
特定しない。FIG. 4 is a block diagram of an interference wave suppressing device having two auxiliary antennas and having an interference wave arrival direction measuring means. FIG.
The description of the common parts is omitted. 20a is an interference wave arrival direction measuring means, and 21 is an antenna pointing direction adjusting means for turning the pointing direction of the auxiliary antenna to the interference wave arrival direction based on the measurement result by the interference wave arrival direction measurement means 20a. The measurement of the direction of arrival of the interference wave is performed, for example, by scanning the main antenna beam over the entire space to detect the direction of arrival of the arriving wave with high power, or by separately providing an antenna for measuring the arrival of the interference wave. FIG. 4 shows the former configuration. If the beam width of the auxiliary antenna is not narrow, precise measurement is not required. The antenna directivity adjusting means 21 directs the main beams of the two auxiliary antennas in the measured arrival direction. If two interference waves are detected, the directions of the two auxiliary antennas are directed to the respective arrival directions of the interference waves. If one interference wave is detected, the directional directions of both of the two auxiliary antennas are directed to their arrival directions. The method of adjusting the directivity of the auxiliary antenna is not specified.
【0189】干渉波到来方向計測手段20aで干渉波数
を検出しているので、実施の形態1で述べたような相関
行列固有値による干渉波数判定は必要ないように思われ
るが、主アンテナのビーム幅より近接した複数干渉波は
干渉波到来方向計測手段20aで検出できないので、相
関行列固有値による干渉波数判定を行う。その後の処理
は実施の形態1と同じである。Since the number of interference waves is detected by the interference wave arrival direction measuring means 20a, it does not seem to be necessary to determine the number of interference waves based on the correlation matrix eigenvalue as described in the first embodiment. Since a plurality of interference waves closer to each other cannot be detected by the interference wave arrival direction measuring means 20a, the number of interference waves is determined based on the correlation matrix eigenvalue. Subsequent processing is the same as in the first embodiment.
【0190】補助アンテナに指向性を持たせることで補
助アンテナ利得を高くでき、補助アンテナ出力における
雑音電力を相対的に低くすることができるため、サイド
ローブキャンセラ出力における信号対雑音電力比の点で
有利となる。By providing directivity to the auxiliary antenna, the gain of the auxiliary antenna can be increased, and the noise power at the output of the auxiliary antenna can be relatively reduced, so that the signal-to-noise power ratio at the output of the sidelobe canceller can be reduced. This is advantageous.
【0191】実施の形態5.指向性を持たせたアレーア
ンテナを補助アンテナとして用いてもよい。ただ、ビー
ム方向はおおまかに干渉波の方向に向ける必要があるた
め、干渉波到来方向計測手段を必要とする。その計測結
果に基づき、補助アンテナとしてのアレーアンテナの荷
重を干渉波方向に向けるように調整する。主アンテナを
アレーアンテナとする場合、補助アンテナ用アレーアン
テナを構成する素子アンテナとしては、主アンテナを構
成する素子アンテナの全てあるいは一部(複数)と共用
してもよいし、主アンテナとは別に設けてもよい。Embodiment 5 FIG. An array antenna having directivity may be used as an auxiliary antenna. However, since the beam direction needs to be roughly directed to the direction of the interference wave, an interference wave arrival direction measuring means is required. Based on the measurement result, the load of the array antenna as the auxiliary antenna is adjusted so as to be directed to the interference wave direction. When the main antenna is an array antenna, the element antenna constituting the auxiliary antenna array antenna may be shared with all or some (plural) of the element antennas constituting the main antenna, or separately from the main antenna. It may be provided.
【0192】図5は補助アンテナを構成するアレーアン
テナが、主アンテナを構成するアレーアンテナの一部と
共用する場合の干渉波到来方向計測手段を備えた補助ア
ンテナ数2の干渉波抑圧装置の構成図である。図1や図
3と共通部分については説明を省く。図5では2つの補
助アンテナを構成する素子アンテナが共通となっている
が、別の素子アンテナを用いてもよい。20bは干渉波
到来方向計測手段、301a,301b,…,301m
と302a,302b,…,302mは2つの補助アン
テナを構成するアレーアンテナが干渉波到来方向計測手
段20bで計測された結果に基づき干渉波到来方向にビ
ームを向けるための荷重、303a,303bはそれぞ
れ301a,301b,…,301mおよび302a,
302b,…,302mで重み付けされた信号を合成す
る重み付け信号合成手段としての合成器、22は301
a,301b,…,301mおよび302a,302
b,…,302mの荷重係数値を計算するための荷重計
算手段である。FIG. 5 shows a configuration of an interference wave suppressing device having two auxiliary antennas and having an interference wave arrival direction measuring means when the array antenna forming the auxiliary antenna is shared with a part of the array antenna forming the main antenna. FIG. 1 and 3 will not be described. In FIG. 5, the element antennas constituting the two auxiliary antennas are common, but different element antennas may be used. Reference numeral 20b denotes an interference wave arrival direction measuring means, 301a, 301b,.
, 302m are loads for the array antennas constituting the two auxiliary antennas to direct the beam in the direction of arrival of the interference wave based on the result measured by the direction-of-arrival measurement means 20b, and 303a and 303b are respectively , 301m and 302a,
A combiner as a weighted signal combining means for combining the signals weighted by 302b,.
, 301m and 302a, 302
b,..., a load calculating means for calculating a load coefficient value of 302 m.
【0193】干渉波到来方向計測手段20bは主アンテ
ナを構成する素子アンテナすべてもしくは一部(複数)
の受信信号を利用して干渉波到来方向計測を行う。図5
では全素子アンテナを用いている。計測方法は特に問わ
ないが、例えば、主アンテナビーム方向を走査したり、
MUSIC法などを使用する。補助アンテナのビーム幅
が狭くなければ精密な計測は必要ない。荷重計算手段2
2では計測された到来方向に2つの補助アンテナの主ビ
ームを向けさせるように荷重301a,301b,…,
301mと302a,302b,…,302mの係数を
計算する。干渉波が2波検出されれば、2つの補助アン
テナビーム方向をそれぞれの干渉波到来方向に向ける。
干渉波が1波検出されれば、2つの補助アンテナ両方の
ビームをその到来方向に向ける。The direction-of-arrival-of-interference measuring means 20b includes all or some (plural) element antennas constituting the main antenna.
The direction of arrival of the interference wave is measured using the received signal. FIG.
Uses all element antennas. The measurement method is not particularly limited, for example, scanning the main antenna beam direction,
The MUSIC method or the like is used. If the beam width of the auxiliary antenna is not narrow, precise measurement is not required. Load calculation means 2
2, the loads 301a, 301b,..., So as to direct the main beams of the two auxiliary antennas in the measured direction of arrival.
The coefficients of 301m and 302a, 302b, ..., 302m are calculated. If two interference waves are detected, the two auxiliary antenna beam directions are directed to the respective arrival directions of the interference waves.
If one interference wave is detected, the beams of both the two auxiliary antennas are directed in the direction of arrival.
【0194】干渉波到来方向計測手段20bで干渉波数
を数えているので、実施の形態1で述べたような相関行
列固有値による干渉波数判定は必要ないように思われる
が、干渉波到来方向の計測分解能が高くない場合は相関
行列固有値による干渉波数判定を判定手段12aで行
う。その後の処理は実施の形態1と同じである。Since the number of interference waves is counted by the interference wave arrival direction measuring means 20b, it does not seem to be necessary to determine the number of interference waves based on the correlation matrix eigenvalue as described in the first embodiment. If the resolution is not high, the determination means 12a determines the number of interference waves based on the correlation matrix eigenvalue. Subsequent processing is the same as in the first embodiment.
【0195】一方、干渉波到来方向の計測分解能が高い
場合は干渉波到来方向計測手段20bで検出した干渉波
数を判定結果としてよい。この場合は相関行列Rの固有
値を計算してそれに基づいて干渉波到来波数を判定する
必要はない。図5の波線による矢印は干渉波到来方向計
測手段20bで検出した干渉波数を判定結果とするとき
の処理の流れを示す。干渉波到来波数判定後の処理は実
施の形態1と同様である。On the other hand, when the measurement resolution of the direction of arrival of the interference wave is high, the number of interference waves detected by the direction of arrival measurement device 20b may be used as the determination result. In this case, there is no need to calculate the eigenvalue of the correlation matrix R and determine the number of incoming interference waves based on the eigenvalue. An arrow indicated by a wavy line in FIG. 5 indicates a flow of processing when the number of interference waves detected by the interference wave arrival direction measuring means 20b is used as a determination result. The processing after the determination of the number of incoming interference waves is the same as in the first embodiment.
【0196】補助アンテナにアレーアンテナを用いるこ
とで補助アンテナ利得を高くでき、補助アンテナ出力に
おける雑音電力を相対的に低くすることができるため、
サイドローブキャンセラ出力における信号対雑音電力比
の点で有利となる。また、補助アンテナを主アンテナを
構成する素子アンテナを共用にすれば、補助アンテナを
主アンテナとは別に設置する必要がないため、スペース
ファクタに優れる。By using an array antenna for the auxiliary antenna, the gain of the auxiliary antenna can be increased, and the noise power at the output of the auxiliary antenna can be relatively reduced.
This is advantageous in terms of the signal-to-noise power ratio at the side lobe canceller output. In addition, if the auxiliary antenna shares the element antenna constituting the main antenna, it is not necessary to install the auxiliary antenna separately from the main antenna, so that the space factor is excellent.
【0197】[0197]
【実施例】以下、上記の実施の形態に係わる実施例を説
明する。ここでは計算機シミュレーションにより上記実
施の形態にかかわる干渉波抑圧装置の有効性を示す。計
算機シミュレーションでは図3の構成を用いた。100
素子半波長間隔リニアアレー、補助アンテナ数2(図3
の左端の素子201aから44.5波長、49.5波長
離れている)で、所望波は0゜(1素子当たり信号対雑
音電力比=0dB、雑音電力を1に正規化)、干渉波は
1波で10゜(1素子当たり干渉信号対雑音電力比=4
0dB、主アンテナサイドローブレベルは約−30d
B)の方向から入射する。角度はリニアアレーのブロー
ドサイド方向を0゜とする。所望波と干渉波は狭帯域
で、所望波と干渉波との間に相関はない。想定する所望
信号電力最大値を雑音電力の3倍とし、干渉波数を1か
2かを判定するしきい値を2×(雑音電力)+2×(所
望信号電力最大値)=2×1+2×(3×1)=8とす
る。EXAMPLES Examples of the above embodiment will be described below. Here, the effectiveness of the interference wave suppressing apparatus according to the above embodiment is shown by computer simulation. In the computer simulation, the configuration shown in FIG. 3 was used. 100
Element half-wavelength linear array, 2 auxiliary antennas (Fig. 3
Is 44.5 wavelengths and 49.5 wavelengths away from the leftmost element 201a), the desired wave is 0 ° (signal-to-noise power ratio = 0 dB per element, noise power is normalized to 1), and the interference wave is 10 ゜ for one wave (interference signal to noise power ratio per element = 4
0dB, main antenna side lobe level is about -30d
Light is incident from the direction B). The angle is 0 ° in the broadside direction of the linear array. The desired wave and the interference wave are narrow bands, and there is no correlation between the desired wave and the interference wave. The assumed maximum value of the desired signal power is set to three times the noise power, and the threshold value for determining whether the number of interference waves is 1 or 2 is 2 × (noise power) + 2 × (maximum desired signal power) = 2 × 1 + 2 × ( 3 × 1) = 8.
【0198】SMI法で使用するサンプル数K=100
0とし、相関行列Rを式(2)、式(4)に従って計算
し、式(45)から相関行列Rの固有値を求めたとこ
ろ、20000および1.32となった。前者はしきい
値よりはるかに大きく、後者はしきい値より小さいので
干渉波数は1と判断する。荷重を式(52)で求め、1
素子当たりの信号対干渉信号+雑音電力比に対するサイ
ドローブキャンセラ出力における信号対干渉信号+雑音
電力比の改善度と、1素子当たりの信号対雑音電力比に
対するサイドローブキャンセラ出力における信号対雑音
電力比の改善度を求めたところ、前者は59.8dB、
後者は19.8dBとなった。荷重係数ベクトルのノル
ムの2乗値は0.048であり、式(34)の雑音電力
増加比は1.05である。また、式(32)の所望信号
電力増加比は1.04で、所望信号の抑圧はなかったこ
とになる。荷重を式(59)で求めた場合も式(52)
で求めた場合とほぼ同じであった。なお、荷重係数ベク
トルのノルムの2乗値は、強制的にw1=0あるいはw2
=0としたときの約1/2であった。The number of samples used in the SMI method K = 100
The correlation matrix R was calculated according to the equations (2) and (4), and the eigenvalues of the correlation matrix R were obtained from the equation (45). Since the former is much larger than the threshold and the latter is smaller than the threshold, the number of interference waves is determined to be one. The load is determined by equation (52), and 1
Improvement of the signal-to-interference signal + noise power ratio at the side lobe canceller output with respect to the signal-to-interference signal-to-noise power ratio per element, and the signal-to-noise power ratio at the side-lobe canceller output to the signal-to-noise power ratio per element When the degree of improvement was determined, the former was 59.8 dB,
The latter is 19.8 dB. The square value of the norm of the weighting coefficient vector is 0.048, and the noise power increase ratio in Expression (34) is 1.05. In addition, the desired signal power increase ratio in Expression (32) is 1.04, which means that the desired signal has not been suppressed. When the load is obtained by equation (59), equation (52)
The results were almost the same as those obtained in The square value of the norm of the load coefficient vector is forcibly set to w 1 = 0 or w 2
= 0, which is about 1/2.
【0199】一方、荷重を式(9)に従って求めた場合
は1素子当たりの信号対干渉信号+雑音電力比に対する
サイドローブキャンセラ出力における信号対干渉信号+
雑音電力比の改善度は44.6dB、1素子当たりの信
号対雑音電力比に対するサイドローブキャンセラ出力に
おける信号対雑音電力比の改善度は4.6dBとなっ
た。信号対雑音電力比の改善度は荷重を式(52)で求
めた場合よりかなり低い。そしてこれは100素子アレ
ーアンテナで得られる理論値(20dB)の1/10以
下である。荷重係数ベクトルのノルムの2乗値は19.
8であり、式(34)の雑音電力増加比は20.8であ
った。つまり、サイドローブキャンセラ出力の雑音電力
が大幅に増加したことになる。そして、式(32)の所
望信号電力増加比は0.59であった。所望信号が一部
抑圧されてしまったことになる。そのために信号対干渉
信号+雑音電力比の改善度は本発明によるものより大き
く低下した。On the other hand, when the load is obtained according to the equation (9), the signal to interference signal at the side lobe canceller output with respect to the signal-to-interference signal plus noise power ratio per element +
The improvement ratio of the noise power ratio was 44.6 dB, and the improvement ratio of the signal-to-noise power ratio at the side lobe canceller output with respect to the signal-to-noise power ratio per element was 4.6 dB. The improvement in the signal-to-noise power ratio is much lower than when the weight is determined by equation (52). This is 1/10 or less of the theoretical value (20 dB) obtained by the 100-element array antenna. The square value of the norm of the load coefficient vector is 19.
8, and the noise power increase ratio in Expression (34) was 20.8. That is, the noise power of the output of the side lobe canceller is significantly increased. Then, the desired signal power increase ratio in Expression (32) was 0.59. This means that the desired signal has been partially suppressed. As a result, the improvement of the signal-to-interference signal + noise power ratio was greatly reduced from that of the present invention.
【0200】この実施例で示すように、干渉波数が補助
アンテナ数より少ない場合は、連立方程式を直接解いて
荷重とするとサイドローブキャンセラ出力において所望
信号の抑圧と雑音電力の増加が起こることがあるのに対
して、本発明では干渉波数が補助アンテナ数より少なく
なってもそのようなことが起こらないことがわかる。As shown in this embodiment, when the number of interference waves is smaller than the number of auxiliary antennas, suppression of a desired signal and an increase in noise power may occur in the output of the sidelobe canceller if the simultaneous equations are directly solved and weighted. On the other hand, according to the present invention, even when the number of interference waves becomes smaller than the number of auxiliary antennas, such a phenomenon does not occur.
【0201】[0201]
【発明の効果】この発明に係る干渉波抑圧装置は、所望
の入射信号の入射方向にビームを向けた主アンテナおよ
び主アンテナに近接して配置された複数の補助アンテナ
ヘの入射信号と同時に入射する干渉波信号を抑圧する干
渉波抑圧装置において、主アンテナおよび補助アンテナ
の各々に設けられ、主アンテナおよび補助アンテナ受信
信号を入力してディジタル同相・直交信号に変換するデ
ィジタル同相・直交信号変換手段と、複数の補助アンテ
ナの受信信号に対応するディジタル同相・直交信号に対
して荷重係数により重み付けする荷重係数重み付け手段
と、重み付けされた複数のディジタル同相・直交信号を
加算する加算手段と、主アンテナの受信信号に対応する
ディジタル同相・直交信号から加算されたディジタル同
相・直交信号を減算する減算手段と、補助アンテナの受
信信号に対応する複数のディジタル同相・直交信号の相
関行列を計算する相関行列計算手段と、主アンテナの受
信信号に対応するディジタル同相・直交信号と補助アン
テナの受信信号に対応する複数のディジタル同相・直交
信号との間の相互相関ベクトルを計算する相互相関ベク
トル計算手段と、相関行列の固有値を計算する相関行列
固有値計算手段と、固有値および相関行列および相互相
関ベクトルに基づいて荷重係数を計算し、荷重係数を荷
重係数重み付け手段に出力する荷重計算手段とを備えて
いる。そのため、干渉波数が補助アンテナ数より少ない
場合、あるいは干渉波数が補助アンテナ数と一致しても
干渉波数が補助アンテナ数より少ないとみなされる場合
でも、荷重係数ベクトルWのノルム||W||が大きくなら
ず、干渉波を十分抑圧し、かつ、サイドローブキャンセ
ラ出力における信号対雑音電力比が低下することがな
い。The interference wave suppressing apparatus according to the present invention simultaneously enters incident signals to a main antenna and a plurality of auxiliary antennas arranged in proximity to the main antenna with a beam directed in the incident direction of a desired incident signal. In an interference wave suppression device for suppressing an interference wave signal, a digital in-phase / quadrature signal conversion unit is provided for each of a main antenna and an auxiliary antenna, and inputs a main antenna and an auxiliary antenna reception signal and converts them into digital in-phase / quadrature signals. Weighting means for weighting the digital in-phase / quadrature signals corresponding to the reception signals of the plurality of auxiliary antennas with a weighting factor; adding means for adding the weighted plurality of digital in-phase / quadrature signals; The added digital in-phase / quadrature signal is subtracted from the digital in-phase / quadrature signal corresponding to the received signal. Subtraction means, a correlation matrix calculation means for calculating a correlation matrix of a plurality of digital in-phase / quadrature signals corresponding to the reception signal of the auxiliary antenna, and reception of the digital in-phase / quadrature signal corresponding to the reception signal of the main antenna and the auxiliary antenna Cross-correlation vector calculation means for calculating a cross-correlation vector between a plurality of digital in-phase / quadrature signals corresponding to a signal, correlation-matrix eigenvalue calculation means for calculating an eigenvalue of a correlation matrix, eigenvalue and correlation matrix and cross-correlation vector And a load calculating means for calculating the load coefficient based on the weighting factor and outputting the load coefficient to the load coefficient weighting means. Therefore, even when the number of interference waves is smaller than the number of auxiliary antennas, or when the number of interference waves is considered to be smaller than the number of auxiliary antennas even when the number of interference waves matches the number of auxiliary antennas, the norm || W || The signal does not increase, the interference wave is sufficiently suppressed, and the signal-to-noise power ratio at the output of the side lobe canceller does not decrease.
【0202】また、主アンテナは、複数の素子アンテナ
から構成されるアレーアンテナである。そのため、スペ
ースファクタに優れる。The main antenna is an array antenna composed of a plurality of element antennas. Therefore, the space factor is excellent.
【0203】また、補助アンテナは、アレーアンテナを
構成する素子アンテナの一部の素子アンテナである。そ
のため、スペースファクタに優れる。The auxiliary antenna is a part of the element antenna constituting the array antenna. Therefore, the space factor is excellent.
【0204】また、補助アンテナは、指向性アンテナで
あり、干渉波の到来方向を計測する干渉波到来方向計測
手段と、干渉波到来方向に補助アンテナの主ビーム方向
を向けるアンテナ指向方向調整手段とをさらに備えてい
る。そのため、補助アンテナに指向性を持たせることで
補助アンテナ利得を高くでき、補助アンテナ出力におけ
る雑音電力を相対的に低くすることができるため、サイ
ドローブキャンセラ出力における信号対雑音電力比の点
で有利となる。The auxiliary antenna is a directional antenna, and includes an interference wave arrival direction measuring means for measuring the arrival direction of the interference wave, and an antenna directivity adjustment means for directing the main beam direction of the auxiliary antenna to the interference wave arrival direction. Is further provided. Therefore, by giving the auxiliary antenna directivity, the auxiliary antenna gain can be increased, and the noise power at the output of the auxiliary antenna can be relatively reduced, which is advantageous in terms of the signal-to-noise power ratio at the output of the sidelobe canceller. Becomes
【0205】また、補助アンテナは、複数の素子アンテ
ナから構成されたアレーアンテナであり、干渉波の到来
方向を計測する干渉波到来方向計測手段と、干渉波到来
方向に補助アンテナの主ビーム方向を向けるように補助
アンテナを構成する各素子アンテナに対して重み付けし
て合成する重み付け信号合成手段とをさらに備えてい
る。そのため、補助アンテナにアレーアンテナを用いる
ことで補助アンテナ利得を高くでき、補助アンテナ出力
における雑音電力を相対的に低くすることができるた
め、サイドローブキャンセラ出力における信号対雑音電
力比の点で有利となる。The auxiliary antenna is an array antenna composed of a plurality of element antennas, and includes an interference wave arrival direction measurement means for measuring the arrival direction of the interference wave, and a main beam direction of the auxiliary antenna for the interference wave arrival direction. Weighted signal combining means for weighting and combining each element antenna constituting the auxiliary antenna so as to direct the antenna is further provided. Therefore, by using an array antenna as the auxiliary antenna, the auxiliary antenna gain can be increased and the noise power at the output of the auxiliary antenna can be relatively reduced, which is advantageous in terms of the signal-to-noise power ratio at the output of the sidelobe canceller. Become.
【0206】さらに、補助アンテナは、主アンテナを構
成する素子アンテナの少なくとも一部で構成され、干渉
波の到来方向を計測する干渉波到来方向計測手段と、干
渉波到来方向に補助アンテナの主ビーム方向を向けるよ
うに補助アンテナを構成する各素子アンテナに対して重
み付けして合成する重み付け信号合成手段とをさらに備
えている。そのため、また、補助アンテナと主アンテナ
を構成する素子アンテナとを共用にすれば、補助アンテ
ナを主アンテナとは別に設置する必要がないため、スペ
ースファクタに優れる。Further, the auxiliary antenna is constituted by at least a part of an element antenna constituting the main antenna, and includes an interference wave arrival direction measuring means for measuring an arrival direction of the interference wave, and a main beam of the auxiliary antenna in the arrival direction of the interference wave. Weighted signal combining means for weighting and combining each element antenna constituting the auxiliary antenna so as to direct the direction is further provided. Therefore, if the auxiliary antenna and the element antenna constituting the main antenna are shared, it is not necessary to install the auxiliary antenna separately from the main antenna, so that the space factor is excellent.
【図1】 本発明の実施の形態1による干渉波抑圧装置
の構成図である。FIG. 1 is a configuration diagram of an interference wave suppression device according to a first embodiment of the present invention.
【図2】 本発明の実施の形態2による干渉波抑圧装置
の構成図である。FIG. 2 is a configuration diagram of an interference wave suppression device according to a second embodiment of the present invention.
【図3】 本発明の実施の形態3による干渉波抑圧装置
の構成図である。FIG. 3 is a configuration diagram of an interference wave suppression device according to a third embodiment of the present invention.
【図4】 本発明の実施の形態4による干渉波抑圧装置
の構成図である。FIG. 4 is a configuration diagram of an interference wave suppression device according to a fourth embodiment of the present invention.
【図5】 本発明の実施の形態5による干渉波抑圧装置
の構成図である。FIG. 5 is a configuration diagram of an interference wave suppression device according to a fifth embodiment of the present invention.
【図6】 従来の干渉波抑圧装置の構成図である。FIG. 6 is a configuration diagram of a conventional interference wave suppression device.
1a〜1d,202a〜202n 受信機(ディジタル
同相・直交信号変換手段)、2a〜2d,203a〜2
03n A/D変換器(ディジタル同相・直交信号変換
手段)、3a〜3c 荷重係数器(荷重係数重み付け手
段)、4 加算器(加算手段)、5 減算器(減算手
段)、10a,10b 相関行列・相互相関ベクトル計
算手段(相関行列計算手段,相互相関ベクトル計算手
段)、11a,11b 相関行列固有値計算手段、13
a,13b,22 荷重計算手段、21 アンテナ指向
方向調整手段、20a,20b 干渉波到来方向計測手
段、303a,303b 合成器(重み付け信号合成手
段)。1a to 1d, 202a to 202n Receiver (digital in-phase / quadrature signal conversion means), 2a to 2d, 203a to 2
03n A / D converter (digital in-phase / quadrature signal conversion means), 3a to 3c weighting coefficient unit (weighting coefficient weighting means), 4 adder (addition means), 5 subtractor (subtraction means), 10a, 10b correlation matrix Cross-correlation vector calculation means (correlation matrix calculation means, cross-correlation vector calculation means), 11a, 11b correlation matrix eigenvalue calculation means, 13
a, 13b, 22 load calculating means, 21 antenna directivity adjusting means, 20a, 20b interference wave arrival direction measuring means, 303a, 303b combiner (weighted signal combining means).
───────────────────────────────────────────────────── フロントページの続き (72)発明者 藤坂 貴彦 東京都千代田区丸の内二丁目2番3号 三 菱電機株式会社内 Fターム(参考) 5J021 AA03 AA04 AA05 AA06 AB09 DA02 DB01 EA04 FA00 FA13 FA14 FA15 FA32 GA02 GA06 JA00 5K052 AA01 BB02 DD04 EE13 FF31 GG19 GG48 5K059 CC03 DD33 DD35 EE02 ────────────────────────────────────────────────── ─── Continuation of the front page (72) Inventor Takahiko Fujisaka 2-3-2 Marunouchi, Chiyoda-ku, Tokyo F-term (reference) 5J021 AA03 AA04 AA05 AA06 AB09 DA02 DB01 EA04 FA00 FA13 FA14 FA15 FA32 GA02 GA06 JA00 5K052 AA01 BB02 DD04 EE13 FF31 GG19 GG48 5K059 CC03 DD33 DD35 EE02
Claims (6)
けた主アンテナおよび主アンテナに近接して配置された
複数の補助アンテナヘの上記入射信号と同時に入射する
干渉波信号を抑圧する干渉波抑圧装置において、 上記主アンテナおよび上記補助アンテナの各々に設けら
れ、上記主アンテナおよび上記補助アンテナ受信信号を
入力してディジタル同相・直交信号に変換するディジタ
ル同相・直交信号変換手段と、 上記複数の補助アンテナの受信信号に対応する上記ディ
ジタル同相・直交信号に対して荷重係数により重み付け
する荷重係数重み付け手段と、 上記重み付けされた複数の上記ディジタル同相・直交信
号を加算する加算手段と、 上記主アンテナの受信信号に対応するディジタル同相・
直交信号から上記加算された上記ディジタル同相・直交
信号を減算する減算手段と、 上記補助アンテナの受信信号に対応する複数の上記ディ
ジタル同相・直交信号の相関行列を計算する相関行列計
算手段と、 上記主アンテナの受信信号に対応する上記ディジタル同
相・直交信号と上記補助アンテナの受信信号に対応する
複数の上記ディジタル同相・直交信号との間の相互相関
ベクトルを計算する相互相関ベクトル計算手段と、 上記相関行列の固有値を計算する相関行列固有値計算手
段と、 上記固有値および上記相関行列および上記相互相関ベク
トルに基づいて上記荷重係数を計算し、荷重係数を上記
荷重係数重み付け手段に出力する荷重計算手段とを備え
たことを特徴とする干渉波抑圧装置。1. An interference wave suppressor for suppressing an interference wave signal simultaneously incident on said main antenna and a plurality of auxiliary antennas disposed in proximity to said main antenna and directing the beam in the incident direction of a desired incident signal. A digital in-phase / quadrature signal conversion means provided in each of the main antenna and the auxiliary antenna for receiving the main antenna and the auxiliary antenna reception signal and converting the received signal into a digital in-phase / quadrature signal; Weighting means for weighting the digital in-phase / quadrature signal corresponding to the reception signal of the antenna by a weighting factor; adding means for adding the plurality of weighted digital in-phase / quadrature signals; Digital in-phase corresponding to the received signal
Subtraction means for subtracting the added digital in-phase / quadrature signal from the quadrature signal; correlation matrix calculation means for calculating a correlation matrix of a plurality of the digital in-phase / quadrature signals corresponding to the reception signal of the auxiliary antenna; Cross-correlation vector calculation means for calculating a cross-correlation vector between the digital in-phase / quadrature signal corresponding to the reception signal of the main antenna and the plurality of digital in-phase / quadrature signals corresponding to the reception signal of the auxiliary antenna; Correlation matrix eigenvalue calculation means for calculating an eigenvalue of a correlation matrix, load calculation means for calculating the load coefficient based on the eigenvalue and the correlation matrix and the cross-correlation vector, and outputting the load coefficient to the load coefficient weighting means An interference wave suppression device comprising:
から構成されるアレーアンテナであることを特徴とする
請求項1記載の干渉波抑圧装置。2. The interference wave suppression device according to claim 1, wherein the main antenna is an array antenna including a plurality of element antennas.
ナを構成する素子アンテナの一部の素子アンテナである
ことを特徴とする請求項2記載の干渉波抑圧装置。3. The interference wave suppression device according to claim 2, wherein the auxiliary antenna is a part of an element antenna constituting the array antenna.
あり、 干渉波の到来方向を計測する干渉波到来方向計測手段
と、 上記干渉波到来方向に上記補助アンテナの主ビーム方向
を向けるアンテナ指向方向調整手段とをさらに備えたこ
とを特徴とする請求項1または2記載の干渉波抑圧装
置。4. The auxiliary antenna is a directional antenna, an interference wave arrival direction measuring means for measuring an arrival direction of an interference wave, and an antenna directivity direction for directing a main beam direction of the auxiliary antenna to the interference wave arrival direction. 3. The interference wave suppressing device according to claim 1, further comprising an adjusting unit.
ナから構成されたアレーアンテナであり、 干渉波の到来方向を計測する干渉波到来方向計測手段
と、 上記干渉波到来方向に上記補助アンテナの主ビーム方向
を向けるように上記補助アンテナを構成する各素子アン
テナに対して重み付けして合成する重み付け信号合成手
段とをさらに備えたことを特徴とする請求項1または2
記載の干渉波抑圧装置。5. The auxiliary antenna is an array antenna composed of a plurality of element antennas, an interference wave arrival direction measuring means for measuring an arrival direction of an interference wave, and a main antenna of the auxiliary antenna in the interference wave arrival direction. 3. A weighted signal synthesizing means for weighting and synthesizing each element antenna constituting the auxiliary antenna so as to direct the beam direction, and further comprising a weighted signal synthesizing means.
An interference wave suppression device as described in the above.
する素子アンテナの少なくとも一部で構成され、 干渉波の到来方向を計測する干渉波到来方向計測手段
と、 上記干渉波到来方向に上記補助アンテナの主ビーム方向
を向けるように上記補助アンテナを構成する各素子アン
テナに対して重み付けして合成する重み付け信号合成手
段とをさらに備えたことを特徴とする請求項3記載の干
渉波抑圧装置。6. The auxiliary antenna comprises at least a part of an element antenna constituting a main antenna, an interference wave arrival direction measuring means for measuring an arrival direction of an interference wave, and the auxiliary antenna in an interference wave arrival direction. 4. The interference wave suppression apparatus according to claim 3, further comprising: a weighted signal combining unit that weights and combines each element antenna constituting the auxiliary antenna so as to direct the main beam direction.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP11181767A JP2000323915A (en) | 1999-03-11 | 1999-06-28 | Interference wave suppressing device |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP6439599 | 1999-03-11 | ||
JP11-64395 | 1999-03-11 | ||
JP11181767A JP2000323915A (en) | 1999-03-11 | 1999-06-28 | Interference wave suppressing device |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2000323915A true JP2000323915A (en) | 2000-11-24 |
Family
ID=26405508
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP11181767A Pending JP2000323915A (en) | 1999-03-11 | 1999-06-28 | Interference wave suppressing device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2000323915A (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010193369A (en) * | 2009-02-20 | 2010-09-02 | Mega Chips Corp | Receiving system and receiving method |
CN103346827A (en) * | 2013-07-26 | 2013-10-09 | 四川大学 | Steady GSC (generalized sidelobe canceler) beam-former and beam optimization method thereof |
JP2014064139A (en) * | 2012-09-20 | 2014-04-10 | Softbank Mobile Corp | Interference measurement system and interference measurement method |
CN104199052A (en) * | 2014-09-22 | 2014-12-10 | 哈尔滨工程大学 | Beam sidelobe suppression method based on norm constraint |
WO2015133252A1 (en) * | 2014-03-05 | 2015-09-11 | 株式会社日立国際電気 | Noise canceller device |
JP2015167298A (en) * | 2014-03-03 | 2015-09-24 | 三菱重工業株式会社 | Signal processor and signal processing method |
JP2016017887A (en) * | 2014-07-09 | 2016-02-01 | ソフトバンク株式会社 | Interference measuring system and interference measuring method |
WO2016194519A1 (en) * | 2015-05-29 | 2016-12-08 | 株式会社日立国際電気 | Noise canceler |
JP2018029272A (en) * | 2016-08-18 | 2018-02-22 | ソフトバンク株式会社 | Base station, interference suppressing device, interference suppressing system and interference suppressing method |
CN117278063A (en) * | 2023-11-22 | 2023-12-22 | 成都通量科技有限公司 | Beam tracking circuit capable of detecting incoming wave direction and interference direction simultaneously |
-
1999
- 1999-06-28 JP JP11181767A patent/JP2000323915A/en active Pending
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010193369A (en) * | 2009-02-20 | 2010-09-02 | Mega Chips Corp | Receiving system and receiving method |
JP2014064139A (en) * | 2012-09-20 | 2014-04-10 | Softbank Mobile Corp | Interference measurement system and interference measurement method |
CN103346827A (en) * | 2013-07-26 | 2013-10-09 | 四川大学 | Steady GSC (generalized sidelobe canceler) beam-former and beam optimization method thereof |
CN103346827B (en) * | 2013-07-26 | 2016-04-06 | 四川大学 | Sane GSC Beam-former and beam optimization method thereof |
JP2015167298A (en) * | 2014-03-03 | 2015-09-24 | 三菱重工業株式会社 | Signal processor and signal processing method |
US9871542B2 (en) | 2014-03-05 | 2018-01-16 | Hitachi Kokusai Electric Inc. | Noise canceller device |
WO2015133252A1 (en) * | 2014-03-05 | 2015-09-11 | 株式会社日立国際電気 | Noise canceller device |
JP2015170877A (en) * | 2014-03-05 | 2015-09-28 | 株式会社日立国際電気 | Noise canceller apparatus |
US20160373145A1 (en) * | 2014-03-05 | 2016-12-22 | Hitachi Kokusai Electric Inc. | Noise canceller device |
JP2016017887A (en) * | 2014-07-09 | 2016-02-01 | ソフトバンク株式会社 | Interference measuring system and interference measuring method |
CN104199052A (en) * | 2014-09-22 | 2014-12-10 | 哈尔滨工程大学 | Beam sidelobe suppression method based on norm constraint |
CN104199052B (en) * | 2014-09-22 | 2016-08-17 | 哈尔滨工程大学 | A kind of beam side lobe suppressing method based on norm constraint |
WO2016194519A1 (en) * | 2015-05-29 | 2016-12-08 | 株式会社日立国際電気 | Noise canceler |
KR20170134759A (en) * | 2015-05-29 | 2017-12-06 | 가부시키가이샤 히다치 고쿠사이 덴키 | Noise canceller device |
CN107615671A (en) * | 2015-05-29 | 2018-01-19 | 株式会社日立国际电气 | Noise eliminator device |
JP6271815B2 (en) * | 2015-05-29 | 2018-01-31 | 株式会社日立国際電気 | Noise canceller device |
KR101878628B1 (en) * | 2015-05-29 | 2018-07-13 | 가부시키가이샤 히다치 고쿠사이 덴키 | Noise canceller device |
US10079614B2 (en) | 2015-05-29 | 2018-09-18 | Hitachi Kokusai Electric Inc. | Noise canceler |
CN107615671B (en) * | 2015-05-29 | 2019-06-25 | 株式会社日立国际电气 | Noise eliminator device |
JP2018029272A (en) * | 2016-08-18 | 2018-02-22 | ソフトバンク株式会社 | Base station, interference suppressing device, interference suppressing system and interference suppressing method |
CN117278063A (en) * | 2023-11-22 | 2023-12-22 | 成都通量科技有限公司 | Beam tracking circuit capable of detecting incoming wave direction and interference direction simultaneously |
CN117278063B (en) * | 2023-11-22 | 2024-02-06 | 成都通量科技有限公司 | Beam tracking circuit capable of detecting incoming wave direction and interference direction simultaneously |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3872953B2 (en) | Wireless communication device using adaptive antenna | |
JP3888189B2 (en) | Adaptive antenna base station equipment | |
JP3738705B2 (en) | Adaptive antenna device | |
US8520784B1 (en) | Coherent beam combining of independently faded signals | |
JP2006270847A (en) | Antenna device | |
JP2000323915A (en) | Interference wave suppressing device | |
Wang et al. | Design of optimum sparse array for robust MVDR beamforming against DOA mismatch | |
Al-Ardi et al. | Direction of arrival estimation in a multipath environment: An overview and a new contribution | |
JP4359778B2 (en) | Adaptive antenna receiver with good directional beam reception quality from the initial stage | |
Wang et al. | Reconfigurable adaptive linear array signal processing in GNSS applications | |
Shubair | Improved smart antenna design using displaced sensor array configuration | |
Shaikh et al. | Linear, non-linear adaptive beamforming algorithm for smart antenna system | |
JP3438527B2 (en) | Signal wave arrival angle estimation device and array antenna control device | |
Satrusallya | Evaluation of Beam Forming Capability of Linear Antenna Array for Smart Antenna System | |
RU2679486C1 (en) | Adaptive antenna system for panoramic radio receiver | |
Chen et al. | Performance analysis of super-resolution beamforming in smart antennas | |
Rashid et al. | Robust 2D beamforming with optimal loading technique | |
Mainkar et al. | Analysis of minimum variance distortionless response and least mean square beamforming algorithm for smart antenna | |
Rajiv et al. | Performance Comparison of Beam Forming Technique Using LMS and SMI algorithms | |
Singh et al. | Comparative Analysis of Conventional and Adaptive beamforming for Linear Array | |
Guruprasad et al. | Reduced complexity beam forming algorithm | |
Mohammed et al. | Estimating Angle of Arrival (AOA) for Wideband Signal by Sensor Delay Line (SDL) and Tapped Delay Line (TDL) Processors | |
Singh et al. | Beamforming Showing Effect on BER with Change in Antenna Configuration | |
SHIRVANI et al. | A new switched-beam setup for adaptive antenna array beamforming | |
Blanchard-Lapierre et al. | The JDTDOA algorithm applied to signal recovery: a performance analysis |