JP2000077926A - Controller for array antenna - Google Patents

Controller for array antenna

Info

Publication number
JP2000077926A
JP2000077926A JP10259177A JP25917798A JP2000077926A JP 2000077926 A JP2000077926 A JP 2000077926A JP 10259177 A JP10259177 A JP 10259177A JP 25917798 A JP25917798 A JP 25917798A JP 2000077926 A JP2000077926 A JP 2000077926A
Authority
JP
Japan
Prior art keywords
signal
weight coefficient
fft
output
band
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP10259177A
Other languages
Japanese (ja)
Other versions
JP3565041B2 (en
Inventor
Yukihiro Kamiya
幸宏 神谷
Yoshio Karasawa
好男 唐沢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
KDDI Corp
Original Assignee
KDD Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by KDD Corp filed Critical KDD Corp
Priority to JP25917798A priority Critical patent/JP3565041B2/en
Publication of JP2000077926A publication Critical patent/JP2000077926A/en
Application granted granted Critical
Publication of JP3565041B2 publication Critical patent/JP3565041B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Abstract

PROBLEM TO BE SOLVED: To improve the delay of convergence caused by low rating of signal processing due to band division by changing the overlap of window of FET in moving on the time base through a frequency band dividing filter means. SOLUTION: From the start of operation to the reception of an instruction from an input/output control means 9, input/output signal control means 3-1 to 3-n succeed samples buffered at every sample timing to FET calculating means 4-1 to 4-n of following steps. After the instruction is received from the input/output control means 9, on the other hand, the input/output signal control means 3-1 to 3-n dispatch sample values to the FET calculating means 4-1 to 4-n only at the timing of completely exchanging the input samples. As a result, the frequency of dispatching values to the FET calculating means is reduced into 1/band division number. Thus, the windows of FET calculating means 4-1 to 4-n buffer signals so as to overlap the samples on the time base.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、各種無線通信に使
用される、素子アンテナの重み係数を適応的に制御する
アレーアンテナの制御装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an array antenna control device for adaptively controlling a weight coefficient of an element antenna used for various wireless communications.

【0002】[0002]

【従来の技術】適応型アレーアンテナは、各素子アンテ
ナで受信した信号に、適当な重み係数を掛け合わせて合
成することにより、ある方向に強い指向性を持たせた
り、逆にある方向から到来する信号に対してのみ、感度
を0にして受信しないようにすることを可能とする。こ
のようなアンテナは、電波環境の変化に応じて適切な重
み係数を用いることにより、最適なビーム形成を行って
空間フィルタリングを実現し、他の無線通信システムと
の干渉を避けるために利用される。
2. Description of the Related Art An adaptive array antenna has a strong directivity in a certain direction or a signal coming from a certain direction by synthesizing a signal received by each element antenna by multiplying the signal by an appropriate weighting factor. It is possible to set the sensitivity to 0 only for the signal to be received and not to receive the signal. Such an antenna is used to perform optimal beamforming and realize spatial filtering by using an appropriate weighting factor according to a change in a radio wave environment, and to avoid interference with other wireless communication systems. .

【0003】近年のディジタル信号処理技術の発展、及
びそれによる装置の小型化に伴い、ディジタル信号処理
を用いた適応型アレーアンテナの移動体通信への応用が
注目を集めている。移動体通信システムへにおいて問題
となるマルチパスフェージング環境では、無数の反射
波、即ち希望信号の遅延信号が到来する。このような環
境に適応型アレーアンテナを適用することにより、理想
的には希望信号と無相関な干渉信号は空間的にキャンセ
ルする一方、希望信号の遅延信号は取り込んで希望信号
に同相合成し、出力信号のS/(N+I) を改善に貢献させる
ことが可能となる。
[0003] With the recent development of digital signal processing technology and the resulting miniaturization of devices, the application of adaptive array antennas using digital signal processing to mobile communications has attracted attention. In a multipath fading environment that poses a problem in a mobile communication system, countless reflected waves, that is, delayed signals of desired signals arrive. By applying the adaptive array antenna to such an environment, ideally, the desired signal and the uncorrelated interference signal are spatially canceled, while the delayed signal of the desired signal is captured and in-phase synthesized with the desired signal, It is possible to contribute to the improvement of S / (N + I) of the output signal.

【0004】遅延信号の同相合成に、遅延素子(TDL ; T
apped Delay Line) を用いた時間軸上の信号処理を行う
構成と、完全再構成フィルタを用いた周波数軸上の信号
処理を行う構成とが有効であることは公知である。それ
ぞれ図1及び図2aに表されている。図1は、TDL に基
づく時間領域信号処理を行う適応型アレーアンテナの構
成図である。また、図2aは、FFTを帯域分割フィル
タとして用いる周波数領域信号処理を行う適応型アレー
アンテナの構成図である。
[0004] A delay element (TDL; T
It is known that a configuration for performing signal processing on the time axis using a mapped delay line) and a configuration for performing signal processing on the frequency axis using a perfect reconstruction filter are effective. 1 and 2a, respectively. FIG. 1 is a configuration diagram of an adaptive array antenna that performs time domain signal processing based on TDL. FIG. 2A is a configuration diagram of an adaptive array antenna that performs frequency domain signal processing using FFT as a band division filter.

【0005】前述の構成を公知とする文献として、R.T.
Compton, Jr, "The RelationshipBetween Tapped Dela
y-Line FFT Processing in Adaptive Array", IEEE
Transaction on Antennas and Propagation, vol. 36 N
o. 1, January 1988 がある。該文献によれば、両構成
の性能に本質的な違いはない。
[0005] As a document that makes the above-mentioned configuration known, RT
Compton, Jr, "The RelationshipBetween Tapped Dela
y-Line FFT Processing in Adaptive Array ", IEEE
Transaction on Antennas and Propagation, vol. 36 N
o. 1, January 1988. According to the document, there is no essential difference in the performance of both configurations.

【0006】図2aの周波数領域信号処理を行う適応型
アレーアンテナは、帯域分割フィルタにFFTを採用す
るが、FFTを施すウィンドウとしてバッファを用いて
いる。このバッファは、時間軸上の信号サンプルに対し
て1サンプルづつシフトしていくように信号を取り込ん
でいく。図2bは、この様子を、時間軸上のサンプルに
対するバッファの動作として表した説明図である。
The adaptive array antenna for performing frequency domain signal processing shown in FIG. 2A employs FFT as a band division filter, but uses a buffer as a window for performing FFT. This buffer takes in the signal so as to shift by one sample with respect to the signal sample on the time axis. FIG. 2B is an explanatory diagram showing this state as the operation of the buffer for the samples on the time axis.

【0007】しかし、特に帯域分割フィルタにFFTを
用いる周波数領域信号処理を行う適応型アレーアンテナ
には、図2aに示す基本構成から派生する複数の構成が
考えられ、それぞれの性能の特徴が調査されている。こ
のような文献として、神谷、唐沢「アダプティブアレー
における時間領域信号処理と周波数領域信号処理の収束
特性に関する特徴比較」, 信学技報 AP98-19 (1998-06)
がある。該文献に述べられている周波数領域信号処理を
行う適応型アレーアンテナの構成には、分割した周波数
帯域(サブバンド)毎に信号処理を並列に行うことによ
り、重み係数更新速度を低減する構成が紹介されてい
る。図3aは、このような、FFTを帯域分割フィルタ
として用いる周波数領域信号処理を行う適応型アレーア
ンテナの構成図であり、重み係数更新速度を低減するも
のである。図3bは、図3aの構成に基づいて、時間軸
上のサンプルに対するバッファの動作を表す説明図であ
る。
However, a plurality of configurations derived from the basic configuration shown in FIG. 2A are conceivable especially for an adaptive array antenna which performs frequency domain signal processing using FFT for a band division filter, and characteristics of each performance are investigated. ing. Kamiya and Karasawa, "Characteristic Comparison of Convergence Characteristics of Time-Domain Signal Processing and Frequency-Domain Signal Processing in Adaptive Arrays", IEICE Technical Report AP98-19 (1998-06)
There is. The configuration of the adaptive array antenna that performs the frequency domain signal processing described in the document includes a configuration in which signal processing is performed in parallel for each divided frequency band (subband), thereby reducing the weight coefficient update speed. Has been introduced. FIG. 3A is a configuration diagram of such an adaptive array antenna that performs frequency domain signal processing using an FFT as a band division filter, and reduces the weight coefficient update speed. FIG. 3B is an explanatory diagram showing the operation of the buffer for the samples on the time axis based on the configuration of FIG. 3A.

【0008】図3aの構成は、帯域分割を行うFFTの
ウィンドウを、信号がオーバーラップしないようにかけ
ていくことにより、各サブバンドにおいては入力信号の
速度がFFTの次数分の1に低減されることになる。重
み係数制御アルゴリズムは、各サブバンドで独立に動作
し、帯域分割された参照信号とサブバンド信号との2乗
誤差を最小にする重み係数を逐次的に求めていく。
In the configuration of FIG. 3A, the speed of the input signal in each subband is reduced to one order of the FFT by applying the FFT window for performing the band division so that the signals do not overlap. become. The weight coefficient control algorithm operates independently for each subband, and sequentially obtains a weight coefficient that minimizes the square error between the band-divided reference signal and the subband signal.

【0009】図3aの構成では、時間領域の信号を復元
するために逆FFTの計算を必要とする。該構成では、
重み係数更新速度を低減できることから、重み係数の適
応制御アルゴリズムが計算機に与える計算負荷を大幅に
低減することができ、装置化にあたって有利な構成とな
ることが期待される。また、このことに伴って、計算負
荷が大きい高機能適応制御アルゴリズムを採用すること
が可能となる。
The configuration of FIG. 3a requires the calculation of an inverse FFT to recover the signal in the time domain. In this configuration,
Since the weight coefficient update speed can be reduced, the calculation load imposed on the computer by the adaptive control algorithm of the weight coefficient can be greatly reduced, and it is expected that the configuration will be advantageous in realizing a device. Accordingly, it is possible to adopt a high-performance adaptive control algorithm having a large calculation load.

【0010】前述のような周波数帯域分割信号処理を行
うアレーアンテナにおいて、受信信号を周波数軸上に分
割する帯域分割フィルタと、この分割した信号をもとの
スペクトルに復元する帯域復元フィルタとの組み合わせ
が必要となるが、このようなフィルタの対は完全再構成
フィルタと呼ばれている。FFT/逆FFTの組み合わ
せは、完全再構成フィルタの一つである。ディジタル信
号処理を行う際には、計算負荷の問題からも、FFT/
逆FFTの組み合わせによる完全再構成フィルタを採用
するのが有効である。
In an array antenna which performs the above-mentioned frequency band division signal processing, a combination of a band division filter for dividing a received signal on a frequency axis and a band restoration filter for restoring the divided signal to an original spectrum. , But such a filter pair is called a perfect reconstruction filter. The combination of FFT / inverse FFT is one of perfect reconstruction filters. When performing digital signal processing, the FFT /
It is effective to employ a perfect reconstruction filter based on a combination of inverse FFT.

【0011】一方、アレーアンテナの適応制御アルゴリ
ズムは、従来から種々検討されている。このうち、RL
S (Recursive Least Squares)法は、高い安定性及び高
速収束性を持つアルゴリズムとして公知である。RLS
法は、例えばSimon Haykin,"Adaptive Filter Theory",
Prentice-Hall Inc. 1996 に詳しく記載されている。
RLS法は、過去のすべてのデータを用い、アレーアン
テナ出力と参照信号との誤差の2乗値を最小とするよう
な重み係数を漸近的に求めるアルゴリズムであり、カル
マンフィルタの一形態であると解釈することもできる。
RLS法は、安定性及び収束速度の点で優れいている一
方、計算負荷が非常に大きいアルゴリズムであるため、
DSP (Digital Signal Processor) を用いる装置化を
行うにあたって、プロセッサへの計算負荷が問題とな
る。
On the other hand, various adaptive control algorithms for array antennas have been conventionally studied. Of these, RL
The S (Recursive Least Squares) method is known as an algorithm having high stability and fast convergence. RLS
The method is, for example, Simon Haykin, "Adaptive Filter Theory",
Prentice-Hall Inc. 1996.
The RLS method is an algorithm that asymptotically obtains a weight coefficient that minimizes the square value of an error between an array antenna output and a reference signal using all past data, and is interpreted as a form of a Kalman filter. You can also.
The RLS method is an algorithm that is excellent in stability and convergence speed, but has a very large calculation load.
When implementing a device using a DSP (Digital Signal Processor), the calculation load on the processor becomes a problem.

【0012】このため、従来から種々のRLS法の計算
量削減方法が提案され、高速RLS法などが提案されて
いるが、これらの計算量低減手法はアレーアンテナには
適用が不可能である。これは、高速RLS法が等化器の
ような、タップ内部の信号が時間的に1タップづつシフ
トしていく性質を利用して構成されているためであり、
アレーアンテナのように、すべてのタップの内部の値が
同時に変化していくシステムには適用が不可能である。
For this reason, various methods for reducing the amount of calculation of the RLS method have been conventionally proposed, and high-speed RLS methods and the like have been proposed, but these calculation amount reduction methods cannot be applied to array antennas. This is because the high-speed RLS method is configured using the property that the signal inside the tap shifts by one tap at a time, such as an equalizer,
It cannot be applied to a system in which the values inside all taps change simultaneously, such as an array antenna.

【0013】図4aは、高速RLS法が、適用可能な実
施形態であり、図4bは、適用不可能な実施形態を表す
説明図である。このように計算量が大きい高機能アルゴ
リズムは、低速動作をさせることができる周波数分割信
号処理型アダプティブアレーと非常に相性が良いといえ
る。
FIG. 4A shows an embodiment to which the fast RLS method can be applied, and FIG. 4B shows an embodiment to which the fast RLS method cannot be applied. It can be said that such a high-performance algorithm having a large calculation amount is very compatible with a frequency-division signal processing type adaptive array capable of operating at a low speed.

【0014】[0014]

【発明が解決しようとする課題】図3aに示した周波数
帯域分割信号処理型アダプティブアレーにRLS法を適
用した場合、分割した帯域毎に並列処理を行うことによ
りRLS法の動作を減速することができるため、このア
レーアンテナの構成とアルゴリズムの組み合わせは装置
化に当たって有効であるといえる。しかし、その一方
で、RLS法を適用したアレーアンテナの構成では、以
下のような問題が生じる。
When the RLS method is applied to the frequency band division signal processing type adaptive array shown in FIG. 3A, the operation of the RLS method can be reduced by performing parallel processing for each divided band. Therefore, it can be said that this combination of the configuration and the algorithm of the array antenna is effective in realizing the device. However, on the other hand, the following problems occur in the configuration of the array antenna to which the RLS method is applied.

【0015】図5は、S/N=10dBにおける、周波数領域信
号処理を行う適応型アレーアンテナの出力信号と参照信
号との相関係数グラフである。アレーアンテナの構成
は、図3aの構成にRLS法を重み係数制御アルゴリズ
ムとして適用したものとなっている。また、アレーアン
テナは1次元アレーで8素子のエレメントアンテナから
なる。受信信号はそれぞれ32帯域に分割される。従っ
て、必要となる重み係数は256個となる。
FIG. 5 is a graph showing a correlation coefficient between an output signal of an adaptive array antenna for performing frequency domain signal processing and a reference signal when S / N = 10 dB. The configuration of the array antenna is obtained by applying the RLS method to the configuration of FIG. 3A as a weight coefficient control algorithm. The array antenna is a one-dimensional array consisting of eight element antennas. Each received signal is divided into 32 bands. Therefore, the required weighting coefficients are 256.

【0016】図5では、横軸はステップ数、即ちRLS
法の繰り返し数を示している。該図5から、重み係数の
探索をスタートし、一度相関係数が改善された後、25
0ステップ前後で大きく劣化することがわかる。このよ
うな現象は、RLS法特有の現象であるものと考えら
れ、その原因は以下のように説明できる。
In FIG. 5, the horizontal axis is the number of steps, that is, RLS.
Shows the number of repetitions of the modulus. From FIG. 5, the search for the weighting coefficient is started, and once the correlation coefficient has been improved, 25
It can be seen that there is a significant deterioration before and after step 0. Such a phenomenon is considered to be a phenomenon peculiar to the RLS method, and the cause can be explained as follows.

【0017】RLS法は、過去のすべての入力を参照
し、その時点における最適な重み係数を求めるが、この
操作は各ステップで1つ得られる方程式を連立方程式と
して逐次的に解を求めていることに等しい。従って、連
立方程式の数が求める重み係数の数以上になるステップ
までは連立方程式の自由度が満たされないので、条件を
満たす解が無数に存在するので、最適重み係数の探索が
開始されない。このため、一度は適当な解へ収束に向か
うが、自由度が満たされるにつれ、その解では条件が満
たされない状態となり、一度誤差が増大したあと、自由
度が満たされた状態で最適重み係数の探索が開始され
る。結果として、最適重み係数の計算が、制御される重
み係数の数と同じステップが過ぎないと開始されない状
態となるため、特に帯域分割信号処理を適用した適応型
アレーアンテナでは収束が大幅に遅くなる原因となる。
The RLS method refers to all past inputs and finds an optimum weighting coefficient at that time. In this operation, an equation obtained in each step is sequentially calculated as a simultaneous equation. Is equivalent to Therefore, since the degree of freedom of the simultaneous equations is not satisfied until the number of simultaneous equations becomes equal to or greater than the number of weight coefficients to be obtained, there are countless solutions satisfying the conditions, and the search for the optimum weight coefficient is not started. For this reason, it once converges to an appropriate solution, but as the degree of freedom is satisfied, the condition is not satisfied with that solution, and once the error increases, the optimal weighting factor is increased with the degree of freedom satisfied. The search starts. As a result, since the calculation of the optimal weighting factor is not started until the same number of steps as the number of weighting factors to be controlled, convergence is significantly slowed down, especially in an adaptive array antenna to which band division signal processing is applied. Cause.

【0018】図3aの構成では、サブバンド毎に独立に
RLS法を適用できるので、自由度が満たされるまでに
必要なサンプル数は素子数に等しい8である。しかし、
重み係数更新速度を低減しているために、サンプルが帯
域分割数である32サンプリング時間に1サンプルしか
入らないため、8×32=256サンプリング時間まで
自由度が満たされないことになっている。図2に示した
周波数領域信号処理を行うアレーアンテナの構成におい
て、相関係数の収束を調べると、やはり同様に256サ
ンプリング時間で一度相関係数が劣化する。このケース
では、サンプルの入力はサンプリング周波数に等しい速
度で入力されるが、すべての重み係数を一つのRLS法
で求めようとしているために、自由度が満たされるまで
に256サンプルを必要とし、結果として図3aのケー
スと同じ時間において相関係数が劣化することとなる。
In the configuration of FIG. 3A, the RLS method can be applied independently for each sub-band, so the number of samples required until the degree of freedom is satisfied is eight, which is equal to the number of elements. But,
Since the weight coefficient update speed is reduced, only one sample is included in 32 sampling times, which is the number of band divisions, so that the degree of freedom is not satisfied up to 8 × 32 = 256 sampling times. In the configuration of the array antenna for performing the frequency domain signal processing shown in FIG. 2, when the convergence of the correlation coefficient is examined, the correlation coefficient is also degraded once in 256 sampling times. In this case, the input of the sample is input at a rate equal to the sampling frequency, but since all the weighting factors are to be obtained by one RLS method, 256 samples are required until the degree of freedom is satisfied, and the result is as follows. As a result, the correlation coefficient deteriorates at the same time as the case of FIG. 3A.

【0019】そこで、本発明は、帯域分割信号処理に基
づく適応型アレーアンテナに、重み係数制御法としてR
LS法を用いた場合、帯域分割による信号処理の低レー
ト化が原因になって生じる収束の遅れを改善するアレー
アンテナの制御装置を提供する。
Accordingly, the present invention provides an adaptive array antenna based on band division signal processing, in which R is used as a weighting factor control method.
Provided is a control device for an array antenna that improves a convergence delay caused by a reduction in signal processing rate due to band division when the LS method is used.

【0020】[0020]

【課題を解決するための手段】本発明によるアレーアン
テナ制御装置は、受信信号を周波数軸上で分割するFF
Tによる帯域分割フィルタ手段と、該手段によって分割
されたそれぞれの周波数帯域の信号に、RLS法によっ
て適応的に制御された重み係数を乗算した後、時間軸上
の信号に復元する帯域復元フィルタ手段とを有してお
り、周波数帯域分割フィルタ手段が、FFTのウィンド
ウの時間軸上の移動におけるオーバーラップを変化させ
るように構成されているものである。
SUMMARY OF THE INVENTION An array antenna control apparatus according to the present invention comprises an FF for dividing a received signal on a frequency axis.
Band dividing filter means by T, and band restoring filter means for restoring a signal on the time axis after multiplying a signal of each frequency band divided by the means by a weight coefficient adaptively controlled by the RLS method And the frequency band division filter means is configured to change the overlap in the movement of the FFT window on the time axis.

【0021】重み係数制御アルゴリズムにRLS法を用
い、FFT及び逆FFTの完全再構成フィルタに基づい
て帯域分割信号処理を行うアレーアンテナの構成では、
重み係数更新速度をFFTの適用によって低減すること
ができる。しかし、低減できることによって自由度が満
たされるまでに時間がかかるので、自由度を満たすまで
の時間は重み係数を低減せず、早期に自由度を満たし、
その後、重み係数更新速度を低減する構成とする。
In the configuration of an array antenna which performs band division signal processing based on perfect reconstruction filters of FFT and inverse FFT using the RLS method as a weight coefficient control algorithm,
The weight coefficient update speed can be reduced by applying the FFT. However, since it takes time until the degree of freedom is satisfied by being able to reduce, the time until the degree of freedom is satisfied does not reduce the weight coefficient, and the degree of freedom is satisfied early,
Thereafter, the weight coefficient update speed is reduced.

【0022】具体的には、自由度が満たされるまではバ
ッファの動作を図2bに示したようにオーバーラップし
てサンプル信号を取り込んで高速にサンプルをRLSプ
ロセッサに供給し、自由度を満たす。こうして自由度を
満たした時点で、バッファの動作を図3bに示したオー
バーラップしないようにウィンドウをかけてFFTに供
給する動作に切り替え、RLSプロセッサの動作速度を
FFTの次数(サブバンドの数)分の1に低減する。
More specifically, the operation of the buffer is overlapped as shown in FIG. 2b until the degree of freedom is satisfied, and the sample signal is fetched and the sample is supplied to the RLS processor at a high speed to satisfy the degree of freedom. When the degree of freedom is satisfied in this way, the operation of the buffer is switched to the operation of supplying the FFT with windowing so as not to overlap as shown in FIG. 3B, and the operation speed of the RLS processor is set to the order of the FFT (the number of subbands). Reduce by a factor of one.

【0023】更に、本発明によるアレーアンテナ制御装
置は、素子アンテナ出力をアナログ信号からディジタル
信号に変換するアナログ/ディジタル信号変換手段と、
該アナログ/ディジタル信号変換手段から出力された信
号を直交信号に分離する信号直交化手段と、該空間アナ
ログ/ディジタル信号変換手段で得られた各素子アンテ
ナ受信信号を、後段の帯域分割FFTに供給するための
ウィンドウを可変オーバーラップでかける信号ベクトル
化手段と、該信号ベクトル化手段で得られた各素子アン
テナ出力がベクトル化された信号を周波数領域でサブバ
ンドに周波数分割を行う帯域分割フィルタと、該FFT
計算手段により選択された所定の数のビーム出力に、後
述する重み係数計算手段から供給される重み係数を乗算
する重み係数乗算手段と、重み係数乗算手段により重み
係数を乗算された各ビーム出力を、同じサブバンド毎に
合計する加算器と、該加算器により合計された各サブバ
ンド毎の出力を処理して者と帯域を復元する信号再構成
フィルタと、前記の構成のうちの任意の点における信号
を用いて適応的に重み係数を計算して重み係数を前記重
み係数乗算器に供給する重み係数計算手段とを有するも
のである。
Further, the array antenna control apparatus according to the present invention comprises: an analog / digital signal conversion means for converting an element antenna output from an analog signal to a digital signal;
A signal orthogonalizing means for separating a signal output from the analog / digital signal converting means into an orthogonal signal, and a receiving signal for each element antenna obtained by the spatial analog / digital signal converting means to a subsequent band division FFT Signal-vectoring means for applying a window for performing variable overlap, and a band division filter for frequency-dividing a signal in which each element antenna output obtained by the signal vectorization means is vectorized into sub-bands in a frequency domain, The FFT
Weighting factor multiplying means for multiplying a predetermined number of beam outputs selected by the calculating means with a weighting factor supplied from a weighting factor calculating means to be described later; and each beam output multiplied by the weighting factor by the weighting factor multiplying means. An adder for summing for each of the same sub-bands, a signal reconstruction filter for processing the output for each of the sub-bands summed by the adder to restore a person and a band, And a weighting factor calculating means for adaptively calculating a weighting factor using the signal in (1) and supplying the weighting factor to the weighting factor multiplier.

【0024】[0024]

【発明の実施の形態】以下、図面を用いて、本発明の実
施形態について詳細に説明する。
Embodiments of the present invention will be described below in detail with reference to the drawings.

【0025】図6は、本発明によるアレーアンテナ制御
装置の一実施形態を表す構成図である。該図6には、複
数のアンテナ素子1-1 〜1-n からなるアレーアンテナ
と、これを制御するための制御装置とが記載されてい
る。
FIG. 6 is a block diagram showing an embodiment of an array antenna control device according to the present invention. FIG. 6 illustrates an array antenna including a plurality of antenna elements 1-1 to 1-n and a control device for controlling the array antenna.

【0026】アンテナ素子1-1 〜1-n で受信された信号
は、A/D 変換及び直交化手段2により所定のサンプリン
グ周波数でディジタル信号に変換され、I/Q 直交チャネ
ルに分離された後、入力信号制御手段3-1 〜3-n に入力
される。入力信号制御手段3-1 〜3-n では、後述の入出
力制御手段9から得た情報に基づいて、RLS法におけ
る自由度を満たすまでの時間においては、後段のFFT
計算手段4-1 〜4-n のウィンドウを、時間軸上でサンプ
ルがオーバーラップするように信号をバッファリングし
てFFT計算手段に信号を供給する。
The signals received by the antenna elements 1-1 to 1-n are converted into digital signals at a predetermined sampling frequency by A / D conversion and orthogonalization means 2, and are separated into I / Q orthogonal channels. Are input to the input signal control means 3-1 to 3-n. In the input signal control means 3-1 to 3-n, based on information obtained from the input / output control means 9 described later, in the time until the degree of freedom in the RLS method is satisfied, the FFT in the subsequent
The signals of the calculation means 4-1 to 4-n are buffered so that the samples overlap on the time axis, and the signals are supplied to the FFT calculation means.

【0027】図7a及び図7bは、図6の実施形態にお
ける入力信号制御手段の動作を表す説明図である。該図
は、入力信号制御手段の動作を、帯域分割数4の場合を
例として説明している。
FIGS. 7A and 7B are explanatory diagrams showing the operation of the input signal control means in the embodiment of FIG. This figure illustrates the operation of the input signal control means by taking the case of four band divisions as an example.

【0028】図7aは、動作を開始してから入出力制御
手段9から命令を受信するまでの、入力信号制御手段3-
1 〜3-n の動作を表す説明図である。即ち、毎サンプル
タイミングでバッファリングしているサンプルを後段の
FFT計算手段4-1 〜4-n に引き渡す。
FIG. 7A shows input signal control means 3 from the start of operation to the reception of a command from input / output control means 9.
It is explanatory drawing showing operation | movement of 1-3 -n. That is, the samples buffered at each sample timing are transferred to the FFT calculation means 4-1 to 4-n at the subsequent stage.

【0029】一方、図7bは、入出力制御手段9から命
令を受信した後の、入力信号制御手段3-1 〜3-n の動作
を表す説明図である。即ち、入力サンプルが完全に入れ
替わったタイミングにおいてのみ、FFT計算手段4-1
〜4-n にサンプル値を引き渡す。結果として、値をFF
T計算手段4-1 〜4-n に引き渡す頻度は、1/帯域分割
数に低減される。
FIG. 7B is an explanatory diagram showing the operation of the input signal control means 3-1 to 3-n after receiving an instruction from the input / output control means 9. That is, only at the timing when the input samples are completely replaced, the FFT calculation means 4-1
Pass sample value to ~ 4-n. As a result, the value is FF
The frequency of handing over to the T calculating means 4-1 to 4-n is reduced to 1 / band division number.

【0030】入力信号制御手段3-1 〜3-n から前述のよ
うなタイミングでサンプル値を供給されるFFT計算手
段4-1 〜4-n は、値を受け取るとFFTの演算を実行し
て出力する。ここで、FFTの次数、即ち帯域分割数は
m とする。FFT計算手段4-1 〜4-n の出力は、乗算器
5-1-1 〜5-m-n に供給され、後述するRLS法に基づく
重み係数計算手段10から供給される重み係数を信号に乗
算され、後段の加算器6-1 〜6-m において同じサブバン
ドの乗算結果毎に加算される。各加算器の出力は、IF
FT(逆FFT)計算手段7に供給される。該IFFT
計算手段7 の結果は、出力信号制御手段8に入力され、
後述の入出力制御手段9からの命令を受けて、前述の入
力信号制御手段3-1 〜3-n の動作と同期して2種類の動
作を行う。即ち、入力信号制御手段の動作が、毎サンプ
ルタイミングでFFTを行っている間は、出力信号制御
手段はIFFTの第1のポートの出力を毎サンプルタイ
ミングで取り込んで出力信号とする。従って、IFFT
の他のポートの出力は無視する。次に、入力信号制御手
段の動作が、サンプリング時間のサブバンド数倍でFF
T計算手段にサンプル値を供給している状態では、出力
信号制御手段はIFFTの全出力を採用し、それを時間
軸上の信号として出力信号に採用する。
The FFT calculating means 4-1 to 4-n, to which the sample values are supplied from the input signal control means 3-1 to 3-n at the timing as described above, execute the FFT operation upon receiving the values. Output. Here, the order of FFT, that is, the number of band divisions is
m. The outputs of the FFT calculation means 4-1 to 4-n are multipliers
The signal is multiplied by a weighting factor supplied from 5-1-1 to 5-mn and supplied from a weighting factor calculation means 10 based on the RLS method described later. It is added for each band multiplication result. The output of each adder is IF
It is supplied to FT (inverse FFT) calculation means 7. The IFFT
The result of the calculation means 7 is input to the output signal control means 8,
Upon receiving a command from the input / output control means 9 described later, two kinds of operations are performed in synchronization with the operations of the input signal control means 3-1 to 3-n. That is, while the operation of the input signal control means is performing the FFT at each sample timing, the output signal control means takes in the output of the first port of the IFFT at each sample timing and makes it an output signal. Therefore, IFFT
Ignore the output of other ports. Next, the operation of the input signal control means determines that the FF
In the state where the sample value is supplied to the T calculation means, the output signal control means adopts the entire output of the IFFT and adopts it as a signal on the time axis as the output signal.

【0031】このように、入力信号制御手段3-1 〜3-n
と出力信号制御手段8のそれぞれの2種類の動作は、入
出力制御手段9によって同期して動作するよう制御され
る。基本的には2種類の動作が切り替わるのは、RLS
法の自由度が満たされるタイミングにおいてであるが、
それより遅いタイミングで切り替えることもあり得る。
As described above, the input signal control means 3-1 to 3-n
The input and output control means 9 controls the two types of operations of the output signal control means 8 and the output signal control means 8 to operate in synchronism. Basically, the two types of operation are switched by RLS
At the timing when the freedom of law is satisfied,
Switching may be performed at a later timing.

【0032】重み係数計算手段(RLS法)10は、前
述の加算器6-1 〜6-m の出力と、対応するサブバンドに
帯域分割された参照信号11とを受信し、2乗誤差が最
小となる重み係数をRLS法に基づいて計算して乗算器
5-1-1 〜5-m-n に供給する。重み係数計算手段10も、
入出力計算手段9から命令を受信し、FFT計算手段の
動作速度と等しい速度でRLS法を実行する。また、帯
域分割された参照信号11も、入出力制御装置9から命
令を受信し、帯域分割するFFTの適用の方法を、FF
T計算手段3-1 〜3-n と同じ動作をさせて参照信号を帯
域分割する。
The weight coefficient calculating means (RLS method) 10 receives the outputs of the adders 6-1 to 6-m and the reference signal 11 divided into corresponding subbands, and the square error is calculated. Multiplier for calculating minimum weighting coefficient based on RLS method
Supply to 5-1-1 to 5-mn. The weight coefficient calculating means 10 also
An instruction is received from the input / output calculation means 9, and the RLS method is executed at a speed equal to the operation speed of the FFT calculation means. Also, the band-divided reference signal 11 receives an instruction from the input / output control device 9 and applies a method of applying FFT to perform band division.
The reference signal is band-divided by performing the same operation as the T calculation means 3-1 to 3-n.

【0033】以上詳細に説明した実施形態について、本
発明の技術思想及び見地の範囲の種々の変更、修正及び
省略は、当業者によれば容易に行うことができる。従っ
て、前述した実施形態では、あくまで例であって、何等
制約しようとするものではない。本発明は、特許請求の
範囲及びその均等物として限定するものだけに制約され
る。
Various changes, modifications, and omissions of the technical idea and scope of the present invention can be easily made by those skilled in the art with respect to the embodiment described in detail above. Therefore, the above-described embodiment is merely an example, and does not intend to limit the present invention. The invention is limited only as defined by the appended claims and equivalents thereof.

【0034】[0034]

【発明の効果】図8aから図8dは、本発明の実施形態
において、各シミュレーションステップで得られる重み
係数から計算される出力信号と参照信号との相関係数グ
ラフである。RLS法の自由度が満たされる32ステッ
プの1倍から4倍までのサンプルを図7aの動作でFF
Tに供給し、その後、図7bの動作に切り替えた場合
の、各ステップにおいて得られた重み係数を用いて計算
したアレーアンテナの出力信号と、参照信号の相関係数
の推移を示す。
FIGS. 8A to 8D are graphs of correlation coefficients between an output signal calculated from weighting coefficients obtained in each simulation step and a reference signal in the embodiment of the present invention. Samples from 1 to 4 times 32 steps satisfying the degree of freedom of the RLS method are subjected to FF by the operation of FIG.
7B shows the transition of the output signal of the array antenna calculated by using the weighting factor obtained in each step and the correlation coefficient of the reference signal when the operation is switched to the operation shown in FIG.

【0035】本発明によれば、図8aから図8dからも
明らかなように、図5に表されたように256ステップ
における相関係数の劣化が見られなくなっており、収束
特性が改善される。
According to the present invention, as is clear from FIGS. 8A to 8D, the deterioration of the correlation coefficient in 256 steps is not seen as shown in FIG. 5, and the convergence characteristic is improved. .

【図面の簡単な説明】[Brief description of the drawings]

【図1】TDL を用いる時間領域信号処理を行う適応型ア
レーアンテナの構成図である。
FIG. 1 is a configuration diagram of an adaptive array antenna that performs time domain signal processing using TDL.

【図2a】FFTを帯域分割フィルタとして用いる周波
数領域信号処理を行う適応型アレーアンテナの構成図で
ある。
FIG. 2A is a configuration diagram of an adaptive array antenna that performs frequency domain signal processing using an FFT as a band division filter.

【図2b】図2aの構成図における、時間軸上のサンプ
ルに対するバッファ動作を表す説明図である。
FIG. 2B is an explanatory diagram showing a buffer operation for a sample on a time axis in the configuration diagram of FIG. 2A.

【図3a】FFTを帯域分割フィルタとして用いる周波
数領域信号処理を行う適応型アレーアンテナにおいて重
み係数更新速度を低減する構成図である。
FIG. 3a is a configuration diagram for reducing a weight coefficient update speed in an adaptive array antenna that performs frequency domain signal processing using an FFT as a band division filter.

【図3b】図3aの構成図における、時間軸上のサンプ
ルに対するバッファ動作を表す説明図である。
FIG. 3B is an explanatory diagram showing a buffer operation for a sample on a time axis in the configuration diagram of FIG. 3A.

【図4a】高速RLS法が適用可能な実施形態の説明図
であって、タップの中を信号がシフトする構成の図であ
る。
FIG. 4A is an explanatory diagram of an embodiment to which the high-speed RLS method can be applied, in which a signal is shifted in a tap;

【図4b】高速RLS法が適用不可能な実施形態の説明
図であって、タップの中を信号が同時に変化する構成の
図である。
FIG. 4B is an explanatory diagram of an embodiment to which the high-speed RLS method cannot be applied, and is a diagram of a configuration in which a signal changes simultaneously in a tap.

【図5】周波数領域信号処理を行う適応型アレーアンテ
ナの出力信号と参照信号との相関係数グラフである。
FIG. 5 is a graph showing a correlation coefficient between an output signal of an adaptive array antenna that performs frequency domain signal processing and a reference signal.

【図6】本発明によるアレーアンテナ制御装置の実施形
態を表す構成図である。
FIG. 6 is a configuration diagram illustrating an embodiment of an array antenna control device according to the present invention.

【図7a】図6の構成において、入力信号制御手段の自
由度が満たされるまでの、毎サンプルタイミングで値を
FFT計算手段に引き渡し(帯域4分割の例)ている動
作を表す説明図である。
FIG. 7A is an explanatory diagram showing an operation of transferring a value to the FFT calculating means at each sample timing (an example of band division into four) until the degree of freedom of the input signal control means is satisfied in the configuration of FIG. 6; .

【図7b】図6の構成において、入力信号制御手段の自
由度が満たされた後の、サンプルタイミングの数倍のタ
イミングで値をFFT計算手段に引き渡し(帯域4分割
の例)ている動作を表す説明図である。
FIG. 7B shows an operation in which, after the degree of freedom of the input signal control means is satisfied, a value is passed to the FFT calculation means at a timing several times the sample timing in the configuration of FIG. 6 (an example of dividing into four bands). FIG.

【図8a】本発明による実施形態について、各シミュレ
ーションステップで得られる重み係数から計算される出
力信号と参照信号との相関係数グラフであって、8サン
プルでバッファ動作を切り替えたグラフである。
FIG. 8A is a graph showing a correlation coefficient between an output signal and a reference signal calculated from a weight coefficient obtained in each simulation step in the embodiment according to the present invention, wherein a buffer operation is switched by eight samples.

【図8b】本発明による実施形態について、各シミュレ
ーションステップで得られる重み係数から計算される出
力信号と参照信号との相関係数グラフであって、16サ
ンプルでバッファ動作を切り替えたグラフである。
FIG. 8B is a graph showing a correlation coefficient between the output signal and the reference signal calculated from the weighting coefficients obtained in each simulation step in the embodiment according to the present invention, wherein the buffer operation is switched by 16 samples.

【図8c】本発明による実施形態について、各シミュレ
ーションステップで得られる重み係数から計算される出
力信号と参照信号との相関係数グラフであって、32サ
ンプルでバッファ動作を切り替えたグラフである。
FIG. 8c is a graph showing a correlation coefficient between an output signal calculated from weighting coefficients obtained in each simulation step and a reference signal in the embodiment according to the present invention, in which a buffer operation is switched by 32 samples.

【図8d】本発明による実施形態について、各シミュレ
ーションステップで得られる重み係数から計算される出
力信号と参照信号との相関係数グラフであって、64サ
ンプルでバッファ動作を切り替えたグラフである。
FIG. 8D is a graph showing the correlation coefficient between the output signal and the reference signal calculated from the weight coefficient obtained in each simulation step in the embodiment according to the present invention, wherein the buffer operation is switched at 64 samples.

【符号の説明】[Explanation of symbols]

1-m アレーアンテナ 2 A/D 及び信号直交化手段 3-n 入力信号制御手段 4-n FFT計算手段(m次) 5-m-n 乗算器 6-m 加算器 7 逆FFT計算手段 8 出力信号制御手段 9 入出力信号制御手段 10 重み係数計算手段(RLS法) 11 帯域分割された参照信号 1-m array antenna 2 A / D and signal orthogonalization means 3-n input signal control means 4-n FFT calculation means (m order) 5-mn multiplier 6-m adder 7 inverse FFT calculation means 8 output signal control Means 9 Input / output signal control means 10 Weight coefficient calculation means (RLS method) 11 Band-divided reference signal

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】 素子アンテナの重み係数を適応的に制御
するアレーアンテナの制御装置であって、受信信号を周
波数軸上で分割するFFTによる帯域分割フィルタ手段
と、該手段によって分割されたそれぞれの周波数帯域の
信号に、RLS法によって適応的に制御された重み係数
を乗算した後、時間軸上の信号に復元する帯域復元フィ
ルタ手段とを有するアレーアンテナの制御装置におい
て、 前記周波数帯域分割フィルタ手段が、FFTのウィンド
ウの時間軸上の移動におけるオーバーラップを変化させ
るように構成されていることを特徴とするアレーアンテ
ナの制御装置。
1. An array antenna control apparatus for adaptively controlling a weight coefficient of an element antenna, comprising: a band division filter means by FFT for dividing a received signal on a frequency axis; A band restoring filter unit for multiplying a frequency band signal by a weighting coefficient adaptively controlled by an RLS method and restoring the signal on a time axis, wherein the frequency band division filter unit Is configured to change the overlap in the movement of the FFT window on the time axis.
【請求項2】 素子アンテナ出力をアナログ信号からデ
ィジタル信号に変換するアナログ/ディジタル信号変換
手段と、 前記アナログ/ディジタル信号変換手段から出力された
信号を直交信号に分離する信号直交化手段と、 前記空間アナログ/ディジタル信号変換手段で得られた
各素子アンテナ受信信号を、後段の帯域分割FFTに供
給するためのウィンドウを可変オーバーラップでかける
信号ベクトル化手段と、 前記信号ベクトル化手段で得られた各素子アンテナ出力
がベクトル化された信号を周波数領域でサブバンドに周
波数分割を行う帯域分割フィルタと、 前記FFT計算手段により選択された所定の数のビーム
出力に、後述する重み係数計算手段から供給される重み
係数を乗算する重み係数乗算手段と、 前記重み係数乗算手段により重み係数を乗算された各ビ
ーム出力を、同じサブバンド毎に合計する加算器と、 前記加算器により合計された各サブバンド毎の出力を処
理して者と帯域を復元する信号再構成フィルタと、 前記の構成のうちの任意の点における信号を用いて適応
的に重み係数を計算して重み係数を前記重み係数乗算器
に供給する重み係数計算手段とを有することを特徴とす
るアレーアンテナの制御装置。
2. An analog / digital signal converter for converting an element antenna output from an analog signal to a digital signal, a signal orthogonalizer for separating a signal output from the analog / digital signal converter into a quadrature signal, A signal vectorizing means for applying a window for variably overlapping a window for supplying each element antenna reception signal obtained by the spatial analog / digital signal conversion means to a subsequent band division FFT; A band division filter that performs frequency division of the signal in which the element antenna output is vectorized into sub-bands in the frequency domain, and a predetermined number of beam outputs selected by the FFT calculation unit are supplied from a weight coefficient calculation unit described later. Weighting factor multiplying means for multiplying the weighting factor by An adder that sums each beam output multiplied by the weight coefficient for each subband, a signal reconstruction filter that processes the output for each subband summed by the adder to restore a person and a band, Weight coefficient calculating means for adaptively calculating a weight coefficient using a signal at an arbitrary point in the configuration and supplying the weight coefficient to the weight coefficient multiplier. Control device.
JP25917798A 1998-08-31 1998-08-31 Array antenna control device Expired - Fee Related JP3565041B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP25917798A JP3565041B2 (en) 1998-08-31 1998-08-31 Array antenna control device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP25917798A JP3565041B2 (en) 1998-08-31 1998-08-31 Array antenna control device

Publications (2)

Publication Number Publication Date
JP2000077926A true JP2000077926A (en) 2000-03-14
JP3565041B2 JP3565041B2 (en) 2004-09-15

Family

ID=17330442

Family Applications (1)

Application Number Title Priority Date Filing Date
JP25917798A Expired - Fee Related JP3565041B2 (en) 1998-08-31 1998-08-31 Array antenna control device

Country Status (1)

Country Link
JP (1) JP3565041B2 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008170221A (en) * 2007-01-10 2008-07-24 Japan Radio Co Ltd Pulse compression radar device
JP2008211686A (en) * 2007-02-27 2008-09-11 Kyocera Corp Wireless communication system, wireless communication device, and wireless communication method
JP2010520720A (en) * 2007-03-02 2010-06-10 クゥアルコム・インコーポレイテッド Using adaptive antenna arrays in conjunction with online repeaters to improve signal quality
JP2012528333A (en) * 2009-05-27 2012-11-12 テレダイン アールディー インスツルメンツ System and method for determining wave characteristics from a mobile platform

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008170221A (en) * 2007-01-10 2008-07-24 Japan Radio Co Ltd Pulse compression radar device
JP2008211686A (en) * 2007-02-27 2008-09-11 Kyocera Corp Wireless communication system, wireless communication device, and wireless communication method
JP2010520720A (en) * 2007-03-02 2010-06-10 クゥアルコム・インコーポレイテッド Using adaptive antenna arrays in conjunction with online repeaters to improve signal quality
US8116239B2 (en) 2007-03-02 2012-02-14 Qualcomm Incorporated Use of a filterbank in an adaptive on-channel repeater utilizing adaptive antenna arrays
US8121535B2 (en) 2007-03-02 2012-02-21 Qualcomm Incorporated Configuration of a repeater
US8599906B2 (en) 2007-03-02 2013-12-03 Qualcomm Incorporated Closed form calculation of temporal equalizer weights used in a repeater transmitter leakage cancellation system
US8619837B2 (en) 2007-03-02 2013-12-31 Qualcomm Incorporated Use of adaptive antenna array in conjunction with an on-channel repeater to improve signal quality
JP2012528333A (en) * 2009-05-27 2012-11-12 テレダイン アールディー インスツルメンツ System and method for determining wave characteristics from a mobile platform
US8654607B2 (en) 2009-05-27 2014-02-18 Teledyne Rd Instruments, Inc. System and method for determining wave characteristics from a moving platform
US9739882B2 (en) 2009-05-27 2017-08-22 Teledyne Instruments, Inc. System and method for determining wave characteristics from a moving platform

Also Published As

Publication number Publication date
JP3565041B2 (en) 2004-09-15

Similar Documents

Publication Publication Date Title
US6590532B1 (en) Radio device
EP1540763B1 (en) Antenna array including virtual antenna elements and method
JP4162522B2 (en) Radio base apparatus, transmission directivity control method, and transmission directivity control program
US6978158B2 (en) Wide-band array antenna
JP2663820B2 (en) Decision feedback equalizer
KR101268691B1 (en) Apparatus and method for receiving data by beamforming in a smart antenna system
JPH07183840A (en) Bi-directional radiocommunication system
WO1997044908A1 (en) Adaptive antenna system and method for cellular and personal communication systems
JP2002094318A (en) Method and device for extracting signal in radio communication system
WO1999039458A1 (en) Method of directional reception using array antenna, and adaptive array antenna unit
JP4406010B2 (en) Apparatus and method for receiving signal in mobile communication system using adaptive antenna array system
US11374613B2 (en) Active interference cancellation apparatus, signal isolation control apparatus and method of actively cancelling interference
JP3565041B2 (en) Array antenna control device
JPH09261008A (en) Filter coefficient arithmetic unit of fir type digital filter for digital beam forming device, fir type digital filter for digital beam forming device and digital beam forming device
JP3473393B2 (en) Adaptive receiver
JP3597694B2 (en) Adaptive antenna device
Huang et al. An efficient subband method for wideband adaptive beamforming
KR20020049357A (en) Method for beamforming using Adaptive Array Antenna in CDMA
JP2008306662A (en) Transmitting/receiving device and communication method thereof
JP3570207B2 (en) Array antenna control method and device
CN113395098B (en) Method and device for multi-antenna signal combination and transmitting signal shaping
Ghavami Wideband beamforming using rectangular arrays without phase shifting
KR20020057594A (en) Apparatus for smart antenna using phase estimation
Li et al. An investigation into broadband smart antenna systems for wireless communication
JP2002271240A (en) Radio receiving device

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040323

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040415

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20040518

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20040531

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees