JPH10229306A - Method and device for controlling adaptive array antenna - Google Patents

Method and device for controlling adaptive array antenna

Info

Publication number
JPH10229306A
JPH10229306A JP9042876A JP4287697A JPH10229306A JP H10229306 A JPH10229306 A JP H10229306A JP 9042876 A JP9042876 A JP 9042876A JP 4287697 A JP4287697 A JP 4287697A JP H10229306 A JPH10229306 A JP H10229306A
Authority
JP
Japan
Prior art keywords
calculating
processing
calculation
digital
signal
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
JP9042876A
Other languages
Japanese (ja)
Other versions
JP3438513B2 (en
Inventor
Yukihiro Kamiya
幸宏 神谷
Shinichi Nomoto
真一 野本
Hideo Kobayashi
英雄 小林
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
Kokusai Denshin Denwa KK
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 Kokusai Denshin Denwa KK filed Critical Kokusai Denshin Denwa KK
Priority to JP04287697A priority Critical patent/JP3438513B2/en
Publication of JPH10229306A publication Critical patent/JPH10229306A/en
Application granted granted Critical
Publication of JP3438513B2 publication Critical patent/JP3438513B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Abstract

PROBLEM TO BE SOLVED: To improve transmission speed by placing priority of at least part of processing in the calculation processing that adaptively calculates a weight coefficient of a signal related to each antenna element to be lower than the priority of other calculation processing for which high-speed signal processing is a requirement so as to reduce the calculation amount. SOLUTION: The priority of at least a part of processing in the calculation processing that adaptively calculates a weight coefficient of a signal related to each of antenna elements 30-1 to 30-m is placed lower than the priority of other calculation processing, for which high speed signal processing is a requirement. Then the calculation processing of a weight coefficient calculation circuit 32, consisting of a covariance matrix calculation section, an inverse matrix calculation section and a weight coefficient calculation section is conducted for a period CT5 which is longer than a sampling period TS, and the calculation processing by multipliers 33-1 to 33-m and an adder 34 is conducted for the sampling period TS. That is, since a plurality of processes whose output frequency of the calculation result differs are processed by a multi-task method, the calculation amount per the sampling period is reduced.

Description

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

【0001】[0001]

【発明の属する技術分野】本発明は、各種無線通信に使
用される適応型アレーアンテナの制御方法及び装置に関
する。
[0001] 1. Field of the Invention [0002] The present invention relates to a method and an apparatus for controlling an adaptive array antenna used for various wireless communications.

【0002】[0002]

【従来の技術】適応型アレーアンテナは、各アンテナ素
子で受信した信号に適当な重み係数を掛け合わせて合成
することにより、ある方向に強い指向性を持たせたり、
逆にある方向から到来する信号に対してのみ、感度をゼ
ロとして受信しないようにすることができるアンテナで
ある。この適応型アレーアンテナは、電波環境の変化に
応じて適切な重み係数を供給することにより、最適なビ
ーム形成を行って空間フィルタリングを実現し、他の無
線通信システムとの干渉を避けるために利用される。
2. Description of the Related Art An adaptive array antenna has a strong directivity in a certain direction by multiplying a signal received by each antenna element by an appropriate weighting coefficient and combining the signals.
Conversely, it is an antenna capable of preventing sensitivity only to a signal arriving from a certain direction with zero sensitivity. This adaptive array antenna provides optimal weighting according to changes in the radio wave environment, performs optimal beamforming, realizes spatial filtering, and is used to avoid interference with other wireless communication systems. Is done.

【0003】このようなアンテナが、電波環境に適応し
たビーム形成を行うためには、重み係数を逐次的に計算
し、更新して行くことが必要となる。この重み係数の適
応的制御アルゴリズムとしては、従来よりさまざまな方
法が提案されている。また、これら制御法の実現方法と
して、アナログ信号をディジタル信号に変換し、計算機
を用いてディジタル信号処理する方法が広く検討されて
いる。
In order for such an antenna to perform beam forming adapted to the radio wave environment, it is necessary to sequentially calculate and update weighting coefficients. Various methods have been conventionally proposed as adaptive control algorithms for the weight coefficient. As a method of implementing these control methods, a method of converting an analog signal into a digital signal and performing digital signal processing using a computer has been widely studied.

【0004】図1は、特開平6−196921号公報に
開示されている従来のフェーズドアレー装置を示す。
FIG. 1 shows a conventional phased array device disclosed in Japanese Patent Application Laid-Open No. 6-196921.

【0005】同図に示すように、複数のアンテナ素子1
0−1〜10−Mからなるアレーアンテナ10により受
信された信号は、送受信モジュールRM−1〜RM−M
内のアナログ/ディジタル変換器11によりディジタル
信号に変換される。変換されたディジタル信号R1 〜R
M はマルチビーム形成回路12において処理されて各ビ
ームに対応する出力E1 〜EN となり、ビーム選択回路
13により、そのうちの電力の大きいビームSE1 〜S
N が選択され、同相分配回路14により分配される。
この出力は分岐した後、その分岐の一方であるSEB1
〜SEBN が適応制御プロセッサ15に入力され、ディ
ジタル信号処理により重み係数W1 〜WN が計算され
る。これら重み係数W1 〜WN は逐次的に更新されて、
可変利得増幅器(乗算器)16−1〜16−Nに供給さ
れ、同相分配回路14の出力のもう一方の分岐である出
力SEA1 〜SEAN と掛け合わされる。その出力が同
相合成器17により合成され、受信ベースバンド信号が
得られる。
[0005] As shown in FIG.
Signals received by the array antenna 10 composed of 0-1 to 10-M are transmitted and received by the transmitting / receiving modules RM-1 to RM-M.
Is converted into a digital signal by the analog / digital converter 11 in the inside. The converted digital signals R 1 to R
M multibeam are processed in forming circuit 12 corresponding to the beam power E 1 to E N, and the by the beam selecting circuit 13, a large beam SE 1 to S of the power of which
E N is selected and distributed by the in-phase distribution circuit 14.
This output is branched and then one of the branches, SEB 1
~SEB N is input to the adaptive control processor 15, the weighting factor W 1 to W-N is calculated by the digital signal processing. These weighting factors W 1 to W N are sequentially updated,
The signals are supplied to variable gain amplifiers (multipliers) 16-1 to 16 -N and multiplied by outputs SEA 1 to SEAN which are the other branches of the output of the in-phase distribution circuit 14. The outputs are combined by the in-phase combiner 17 to obtain a received baseband signal.

【0006】送信においては、適応制御プロセッサ15
の出力W1 〜WN に応じて位相演算プロセッサ18で算
出した位相量DP1 〜DPM が位相器19に供給され、
送信ベースバンド信号を同相分配器20で同相分配した
送信信号F1 〜FM の位相制御が施され送信ビームが形
成される。
In transmission, the adaptive control processor 15
Phase amount DP 1 to DP M calculated by the phase calculating processor 18 in response to the output W 1 to W-N of is supplied to the phaser 19,
The transmission signals F 1 to F M obtained by in-phase distribution of the transmission baseband signal by the in-phase distributor 20 are subjected to phase control to form a transmission beam.

【0007】[0007]

【発明が解決しようとする課題】ディジタル信号処理を
するこの種の制御装置は、装置の小型化などのメリット
がある反面、重み係数を算出する計算機へ多大な負荷が
かかる問題があり、伝送速度の面からその適用範囲が制
限されている。即ち、このような適応型アレーアンテナ
制御装置は、その動作速度が、駆動する計算機の計算速
度に依存するため、結果として信号伝送速度も計算機の
能力に制限されるので、その適用が限られてしまうので
ある。
This type of control device for performing digital signal processing has the advantage of downsizing the device and the like, but has the problem that a heavy load is imposed on the computer for calculating the weighting factor, and the transmission speed is high. Therefore, its application range is limited. That is, since the operation speed of such an adaptive array antenna control device depends on the calculation speed of the computer to be driven, the signal transmission speed is also limited by the capability of the computer, so its application is limited. It will be lost.

【0008】例えば、図1に示した従来例では、ディジ
タル信号処理を行う各部の処理周期は基本的に同一であ
る。即ち、各部の処理優先度は互いに等しい。このよう
なディジタル信号処理システムの伝送速度は、アナログ
信号をディジタル信号に変換する際のサンプリング周波
数により、サンプリング定理の関係から決定されるが、
従来例では、各部の処理優先度が互いに等しいので、こ
のシステムを用いてリアルタイムに適応制御を行おうと
すると、各サンプル値に対して全ての計算プロセスを実
行してから次のサンプル値を得ることとなる。そのた
め、計算プロセスが全て実行される時間より短いサンプ
リング周期を実現することができなくなり、結果とし
て、計算機の処理速度によって伝送速度が規定されてし
まうこととなる。
For example, in the conventional example shown in FIG. 1, the processing cycle of each unit for performing digital signal processing is basically the same. That is, the processing priorities of the respective units are equal to each other. The transmission speed of such a digital signal processing system is determined by the sampling frequency when converting an analog signal to a digital signal from the relationship of the sampling theorem.
In the conventional example, since the processing priorities of the respective units are equal to each other, when performing an adaptive control in real time using this system, it is necessary to execute all calculation processes for each sample value and then obtain the next sample value. Becomes For this reason, it is impossible to realize a sampling period shorter than the time during which the entire calculation process is executed, and as a result, the transmission speed is defined by the processing speed of the computer.

【0009】一方、電波環境の変化速度(例えば、信号
波源及び干渉波源の移動速度、移動体端末の移動速度)
は、一般に、ディジタル信号処理のサンプリング速度に
対して十分遅い。このことを考慮すると、前述した従来
例におけるビーム選択手段及び適応制御プロセッサの動
作速度、即ち重み係数の更新速度は、それほど高速であ
る必要はなく、サンプリング間隔で重みを更新する必要
はないこととなる。従って、同一の計算優先度(処理優
先度)で全ての処理を行うことは、不必要な計算負荷を
ディジタル信号処理部に与え、結果としてサンプリング
速度を低下させ、伝送速度を低下させることになる。
On the other hand, the changing speed of the radio wave environment (for example, the moving speed of the signal wave source and the interference wave source, the moving speed of the mobile terminal)
Is generally slow enough for the sampling rate of digital signal processing. In consideration of this, the operation speed of the beam selecting means and the adaptive control processor in the conventional example described above, that is, the update speed of the weight coefficient does not need to be so high, and it is not necessary to update the weight at the sampling interval. Become. Therefore, performing all the processes with the same calculation priority (processing priority) gives an unnecessary calculation load to the digital signal processing unit, and as a result, lowers the sampling speed and lowers the transmission speed. .

【0010】このように従来の適応型アレーアンテナ重
み係数の計算方法によると、ディジタル信号処理装置へ
の計算負荷によりそのサンプリング周期が制限され、結
果として伝送速度が制限されるという問題が生じてい
た。
As described above, according to the conventional method for calculating the weight coefficient of the adaptive array antenna, there is a problem that the sampling period is limited by the calculation load on the digital signal processing device, and as a result, the transmission speed is limited. .

【0011】従って本発明の目的は、現状の限られたデ
ィジタル信号処理を行う計算機の計算能力をより効率的
に用いて、より高速な伝送速度に対して適応するアレー
アンテナ制御方法及び装置を提供することにある。
Accordingly, it is an object of the present invention to provide an array antenna control method and apparatus which adapts to a higher transmission rate by more efficiently using the computing power of a computer which performs current limited digital signal processing. Is to do.

【0012】[0012]

【課題を解決するための手段】本発明によれば、複数の
アンテナ素子を有する適応型アレーアンテナの制御をデ
ィジタル信号処理を用いて行う場合に、各アンテナ素子
に関する信号の重み係数を適応的に計算する計算処理部
における少なくとも一部の処理の優先度を、信号の高速
処理が要求される他の計算処理の優先度より低く設定し
た適応型アレーアンテナの制御方法及び装置が提供され
る。
According to the present invention, when control of an adaptive array antenna having a plurality of antenna elements is performed using digital signal processing, the weight coefficient of a signal relating to each antenna element is adaptively adjusted. There is provided an adaptive array antenna control method and apparatus in which the priority of at least a part of processing in a calculation processing unit for calculation is set lower than the priority of other calculation processing requiring high-speed signal processing.

【0013】通常、アレーアンテナが適用される環境に
おいて、信号源及び干渉源並びに受信点の移動速度等の
電波環境の幾何学的変化速度は、ディジタル信号処理の
サンプリング速度と比較すると非常に低速である。従っ
て、重み係数は、サンプリング間隔毎に更新される必要
はなく、ある一定の比較的長い周期で重み係数の保持と
更新を繰り返せば、実用上十分である。ただし、漸近的
に重み係数を求めるアルゴリズムにおいて、重み係数は
サンプル値毎に計算するが、目標値は一定の比較的長い
周期で保持と更新を繰り返す。このように、電波環境の
変化の速度を考慮すると、アレーアンテナ制御装置のう
ち、計算結果の高速な更新が必要ない部分を分割して実
行することにより、サンプリング周期あたりの計算量を
軽減し、結果として同じ処理能力のプロセッサを用いた
場合に、より高速な伝送速度の実現が可能となる。
Normally, in an environment to which an array antenna is applied, the geometric change speed of the radio wave environment such as the moving speed of a signal source and an interference source and a receiving point is very low as compared with the sampling speed of digital signal processing. is there. Therefore, the weight coefficient does not need to be updated at each sampling interval, and it is practically sufficient to repeat the holding and updating of the weight coefficient at a certain relatively long period. However, in the algorithm for asymptotically obtaining the weighting coefficient, the weighting coefficient is calculated for each sample value, but the target value is repeatedly held and updated at a constant relatively long cycle. As described above, in consideration of the speed of the change of the radio wave environment, by dividing and executing a portion of the array antenna control device that does not require a high-speed update of the calculation result, the calculation amount per sampling cycle is reduced, As a result, when processors having the same processing capacity are used, a higher transmission speed can be realized.

【0014】図2は、例として時間軸上に示した計算の
更新頻度の異なる2つの計算プロセスを用いて本発明の
基本原理を説明するための図である。この例では、プロ
セス1とプロセス2との2つの計算プロセスがあり、前
者は後者よりも優先度が高い、即ち、プロセス2よりプ
ロセス1の方が、計算結果の高速な更新が必要であるも
のとする。
FIG. 2 is a diagram for explaining the basic principle of the present invention using two calculation processes having different update frequencies of the calculation shown on the time axis as an example. In this example, there are two calculation processes, a process 1 and a process 2, and the former has a higher priority than the latter, that is, the process 1 requires a faster update of the calculation result than the process 2. And

【0015】同図(a)は、従来技術の信号処理により
行われてきたシングルタスク法を示しており、全ての計
算プロセスを同一の優先度で行っている。このため、時
間的に優先度の高いプロセス1もプロセス2も同じ更新
頻度で計算が進められ、各サンプル時間間隔で全ての計
算を行ってから次のサンプル値を得るので、プロセス1
の計算を高速化するためには、計算機の性能を改善する
しかない。
FIG. 1A shows a single task method which has been performed by the conventional signal processing, and all calculation processes are performed with the same priority. Therefore, the process 1 and the process 2 having higher temporal priority are calculated at the same update frequency, and the next sample value is obtained after performing all calculations at each sample time interval.
The only way to speed up the calculation is to improve the performance of the computer.

【0016】一方、同図(b)は、本発明の基本原理に
関するマルチタスク法を示しており、計算プロセスがタ
スク1及びタスク2の2つに分割されている。高速計算
が必要なプロセス1は各サンプリング時間間隔以内に計
算を完了し、各サンプリング間隔で結果を算出する必要
がないプロセス2については、各々のサンプリング間隔
でプロセス1の計算を行った残りの時間を使って行うこ
ととし、プロセス2の計算が途中でもサンプリング間隔
時間が終了すると、計算を中断し、次のサンプル値を取
り込んでプロセス1の計算を新たに開始する。プロセス
1の計算が完了した後、再び前のサンプリング時間間隔
で途中であったプロセス2の計算を再開する。以上のよ
うに、高速な計算結果出力が必要ない部分の計算処理
を、分割して実行することにより、サンプリング周期当
たりの計算量を軽減し、結果として、同じ処理能力のプ
ロセッサを用いた場合に伝送速度の向上が可能となる。
FIG. 1B shows a multitask method according to the basic principle of the present invention, in which the calculation process is divided into two tasks, task 1 and task 2. Process 1 which requires high-speed calculation completes the calculation within each sampling time interval, and process 2 which does not need to calculate the result at each sampling interval, the remaining time for calculating process 1 at each sampling interval. When the sampling interval time ends even during the calculation of the process 2, the calculation is interrupted, the next sample value is fetched, and the calculation of the process 1 is newly started. After the calculation of the process 1 is completed, the calculation of the process 2 which was halfway at the previous sampling time interval is restarted. As described above, by dividing and executing the calculation processing of the portion that does not require high-speed calculation result output, the calculation amount per sampling cycle is reduced, and as a result, when a processor having the same processing capacity is used. Transmission speed can be improved.

【0017】本発明に関する方法が、アンテナ素子で受
信した複数のアナログ信号をディジタル信号に変換し、
変換された複数のディジタル信号に基づいて複数のアン
テナ素子対応の重み係数を適応的に計算し、計算した複
数の重み係数と複数のディジタル信号とをそれぞれ乗算
した後に合成する処理を含んでいる場合に、高速処理が
要求される他の計算処理が、ディジタル信号への変換処
理、乗算処理及び合成処理を少なくとも含んでいること
が好ましい。
A method according to the invention comprises converting a plurality of analog signals received by an antenna element into digital signals,
In the case of including a process of adaptively calculating weighting factors corresponding to a plurality of antenna elements based on a plurality of converted digital signals, multiplying the calculated plurality of weighting factors by a plurality of digital signals, and then combining them. In addition, it is preferable that other calculation processing requiring high-speed processing includes at least conversion processing into a digital signal, multiplication processing, and synthesis processing.

【0018】高速処理が要求される他の計算処理が、デ
ィジタル信号への変換処理のサンプリング周期内に完了
されることが好ましい。
It is preferable that other calculation processing requiring high-speed processing be completed within the sampling period of the conversion processing into a digital signal.

【0019】重み係数の計算処理が、ディジタル信号へ
の変換処理のサンプリング周期より長い周期で完了され
ることが好ましい。
It is preferable that the calculation processing of the weight coefficient is completed in a cycle longer than the sampling cycle of the conversion processing into the digital signal.

【0020】複数のアンテナ素子を有する適応型アレー
アンテナの制御をディジタル信号処理を用いて行う本発
明に関する装置は、アンテナ素子で受信した複数のアナ
ログ信号をディジタル信号に変換するアナログ/ディジ
タル変換手段と、変換された複数のディジタル信号に基
づいて複数のアンテナ素子対応の重み係数を適応的に計
算する重み係数計算手段と、重み係数計算手段から出力
される複数の重み係数と前記変換手段から得られる複数
のディジタル信号とをそれぞれ乗算する乗算手段と、乗
算手段の出力を合成する合成手段とを備えており、重み
係数計算手段における計算処理における少なくとも一部
の処理の優先度を、高速処理が要求されるアナログ/デ
ィジタル変換手段、乗算手段及び合成手段の計算処理の
優先度より低く設定することが好ましい。
An apparatus according to the present invention for controlling an adaptive array antenna having a plurality of antenna elements by using digital signal processing comprises an analog / digital conversion means for converting a plurality of analog signals received by the antenna elements into digital signals. Weight coefficient calculating means for adaptively calculating a weight coefficient corresponding to a plurality of antenna elements based on a plurality of converted digital signals; a plurality of weight coefficients output from the weight coefficient calculating means; A multiplying means for multiplying each of the plurality of digital signals; and a synthesizing means for synthesizing an output of the multiplying means. High-speed processing requires priority of at least a part of processing in the weight coefficient calculating means. Lower than the priority of the calculation processing of the analog / digital conversion means, the multiplication means and the synthesis means. It is preferable to.

【0021】本発明の一実施形態においては、重み係数
計算手段が、アナログ/ディジタル変換手段から得られ
た複数のサンプル値の共分散行列の時間平均値を、アナ
ログ/ディジタル変換手段のサンプリング周期より長い
周期で算出する第1の計算手段と、第1の計算手段から
得られた共分散行列の逆行列を、サンプリング周期より
長い周期で算出する第2の計算手段と、第2の計算手段
により得られた信号行列の逆行列からアレーアンテナの
重み係数を算出し、サンプリング周期より長い所定時間
の間、重み係数を保持すると共に、所定時間を周期とし
て重み係数を更新する第3の計算手段とを含んでいる。
In one embodiment of the present invention, the weight coefficient calculating means calculates a time average value of a covariance matrix of a plurality of sample values obtained from the analog / digital converting means based on a sampling period of the analog / digital converting means. A first calculating means for calculating in a long cycle, a second calculating means for calculating an inverse matrix of a covariance matrix obtained from the first calculating means in a cycle longer than a sampling cycle, and a second calculating means. A third calculating means for calculating a weighting factor of the array antenna from the inverse matrix of the obtained signal matrix, holding the weighting factor for a predetermined time longer than the sampling period, and updating the weighting factor at a predetermined time period; Contains.

【0022】本発明の他の実施形態においては、重み係
数計算手段が、アナログ/ディジタル変換手段のサンプ
リング周期毎にアナログ/ディジタル変換手段から入力
される複数のサンプル値の共分散行列の時間平均値を求
めて信号波到来角度をサンプリング周期より長い周期で
算出する第1の計算手段と、第1の計算手段により得ら
れた信号波到来角度に基づいてアレーアンテナの重み係
数を算出し、サンプリング周期より長い所定時間の間、
重み係数を保持すると共に、所定時間を周期として重み
係数を更新する第2の計算手段とを含んでいる。
In another embodiment of the present invention, the weighting factor calculating means includes a time average value of a covariance matrix of a plurality of sample values input from the analog / digital converting means for each sampling cycle of the analog / digital converting means. Calculating a signal wave arrival angle at a period longer than the sampling period, and calculating a weighting factor of the array antenna based on the signal wave arrival angle obtained by the first calculation unit. For a longer period of time,
A second calculating unit that holds the weighting coefficient and updates the weighting coefficient with a predetermined time as a cycle.

【0023】この場合、第1の計算手段が、アナログ/
ディジタル変換手段から得られた複数のサンプル値の共
分散行列の時間平均値を、アナログ/ディジタル変換手
段のサンプリング周期より長い周期で算出する第3の計
算手段と、第3の計算手段から得られた共分散行列の固
有値及び固有ベクトルを、サンプリング周期より長い周
期で算出する第4の計算手段と、第4の計算手段により
得られた信号行列の固有値及び固有ベクトルから信号到
来角度を求める第5の計算手段とを含んでいるかもしれ
ない。
In this case, the first calculating means uses analog /
A third calculating means for calculating a time average value of a covariance matrix of a plurality of sample values obtained from the digital converting means at a cycle longer than a sampling cycle of the analog / digital converting means; and a third calculating means. Fourth calculating means for calculating the eigenvalues and eigenvectors of the covariance matrix at a period longer than the sampling period, and fifth calculation for obtaining the signal arrival angle from the eigenvalues and eigenvectors of the signal matrix obtained by the fourth calculating means. Means.

【0024】本発明のさらに他の実施形態においては、
重み係数計算手段が、アナログ/ディジタル変換手段の
サンプリング周期毎にアナログ/ディジタル変換手段か
ら入力される複数のサンプル値の共分散行列の時間平均
値を求めて信号波到来角度をサンプリング周期より長い
周期で算出し保持する第1の計算手段と、第1の計算手
段により得られた信号波到来角度と合成手段からフィー
ドバックされる受信信号と各素子出力のサンプル値とか
ら、アレーアンテナの重み係数を、サンプリング周期で
更新する第2の計算手段とを含んでいる。
In still another embodiment of the present invention,
The weighting factor calculating means calculates a time average value of a covariance matrix of a plurality of sample values input from the analog / digital conversion means for each sampling cycle of the analog / digital conversion means, and determines a signal wave arrival angle in a cycle longer than the sampling cycle. A first calculating means which calculates and holds the weighting coefficient of the array antenna based on a signal wave arrival angle obtained by the first calculating means, a received signal fed back from the synthesizing means, and a sample value of each element output. , A second calculating means for updating at a sampling period.

【0025】[0025]

【発明の実施の形態】図3は、本発明の一実施形態とし
て、適応型アレーアンテナの制御装置を概略的に示すブ
ロック図である。同図において、30は複数のアンテナ
素子30−1〜30−mを有するアレーアンテナ、31
−1〜31−mはアンテナ素子30−1〜30−mから
のアナログ出力をディジタル信号に変換するアナログ/
ディジタル変換器、32はアナログ/ディジタル変換器
31−1〜31−mの出力及び足し算器34のフィード
バック出力から重み係数を算出する重み係数計算回路、
33−1〜33−mはアナログ/ディジタル変換器31
−1〜31−mの出力に重み係数を乗算する乗算器、3
4は乗算器33−1〜33−mの出力を合成する足し算
器をそれぞれ示している。
FIG. 3 is a block diagram schematically showing a control device of an adaptive array antenna according to an embodiment of the present invention. In the figure, reference numeral 30 denotes an array antenna having a plurality of antenna elements 30-1 to 30-m;
-1 to 31-m are analog / digital converters for converting analog outputs from the antenna elements 30-1 to 30-m into digital signals.
A digital converter 32, a weight coefficient calculating circuit for calculating a weight coefficient from the outputs of the analog / digital converters 31-1 to 31-m and the feedback output of the adder 34;
33-1 to 33-m are analog / digital converters 31
A multiplier for multiplying the outputs of −1 to 31-m by a weighting factor, 3
Reference numeral 4 denotes an adder for combining the outputs of the multipliers 33-1 to 33-m.

【0026】アンテナ素子30−1〜30−mからのア
ナログ出力は、各素子に対応するアナログ/ディジタル
変換器31−1〜31−mにより、所定のサンプリング
周期TS でサンプリングされ、ディジタル信号に変換さ
れる。変換されたディジタル信号は、対応する重み係数
を乗算器33−1〜33−mで乗算され、足し算器34
により各アンテナ素子に対応する出力が合成される。こ
れらの計算は、リアルタイムの適応制御を実現するため
に、所定のサンプリング周期TS 以内に結果が算出され
なければならない。この結果は、同様に、サンプリング
周期毎に重み係数を出力する重み係数計算回路32へフ
ィードバックされ、重み係数が計算されて乗算器33−
1〜33−mに供給される。即ち、本発明における信号
の高速処理が要求される計算処理は、本実施形態では、
アナログ/ディジタル変換器31−1〜31−m、乗算
器33−1〜33−m及び足し算器34の計算処理に対
応している。
The analog outputs from the antenna elements 30-1 to 30-m are sampled at predetermined sampling periods T S by analog / digital converters 31-1 to 31-m corresponding to each element, and are converted into digital signals. Is converted. The converted digital signal is multiplied by a corresponding weight coefficient in multipliers 33-1 to 33-m, and an adder 34
As a result, outputs corresponding to the respective antenna elements are combined. The results of these calculations must be calculated within a predetermined sampling period T S in order to realize real-time adaptive control. The result is similarly fed back to a weighting factor calculation circuit 32 that outputs a weighting factor for each sampling period, a weighting factor is calculated, and a multiplier 33-
1-33-m. That is, the calculation processing that requires high-speed processing of the signal in the present invention is, in the present embodiment,
It corresponds to the calculation processing of the analog / digital converters 31-1 to 31-m, the multipliers 33-1 to 33-m, and the adder 34.

【0027】重み係数計算回路32は、サンプリング周
期TS より長い周期CTS (ただし、Cは1より大きい
定数)で重み係数を算出して更新し、この周期CTS
では、前回に算出した値を保持して出力する。なお、図
3に記載されているTS 及びCTS は、その周期で信号
が入出力されることを表わしている。
The weighting coefficient calculating circuit 32, the sampling period T S is longer than the period CT S (however, C is a constant greater than 1) and update calculates a weight coefficient, within this period CT S, was calculated last time Output while retaining the value. It should be noted that T S and C S S shown in FIG. 3 indicate that signals are input and output in that cycle.

【0028】図4は、図3の実施形態における重み係数
計算回路32の構成例を示すブロック図である。この構
成例では、入力信号行列の逆行列を用いて重み係数を求
めている。
FIG. 4 is a block diagram showing an example of the configuration of the weight coefficient calculation circuit 32 in the embodiment of FIG. In this configuration example, the weight coefficient is obtained using the inverse matrix of the input signal matrix.

【0029】入力信号の共分散行列の逆行列を用いる方
法としては、例えばSMI(Sample Matri
x Inverse)アルゴリズムがある。この制御ア
ルゴリズムは、例えば「SMI法に基づくアダプティブ
アレーを用いた高速ディジタル陸上移動通信の多重波抑
圧」、小川・川端・伊藤、電子情報通信学会論文誌、B
−II Vol.J75−B−II No.11 19
92年11月に詳しく説明されている。
As a method of using an inverse matrix of a covariance matrix of an input signal, for example, an SMI (Sample Matrix)
x Inverse) algorithm. This control algorithm is described in, for example, "Multiple Wave Suppression in High-Speed Digital Land Mobile Communication Using an Adaptive Array Based on the SMI Method", Ogawa / Kawabata / Ito, Transactions of the Institute of Electronics, Information and Communication Engineers, B
-II Vol. J75-B-II No. 11 19
This is described in detail in November 1992.

【0030】このSMIアルゴリズムによれば、最適重
み係数Wopt は、次のように表わされる。 Wopt=Rxx -1xr ここで、Rxx及びVxrは、xk (i)をk番目のアンテ
ナ素子の時刻iにおけるサンプリング値、r(i)を各
バーストの先頭に置かれるトレーニング信号系、X
(i)=[x1 (i)x2 (i)・・・xm (i)]T
とし、T が転置、* が複素共役、E[ ]が期待値(ア
ンサンブル平均)を表わすものとすると、 Rxx=E[X* (i)XT (i)] Vxr=E[X* (i)r(i)] である。
According to the SMI algorithm, the optimum weight coefficient W opt is expressed as follows. W opt = R xx -1 V xr where R xx and V xr are trainings in which x k (i) is a sampling value of the k-th antenna element at time i, and r (i) is placed at the head of each burst. Signal system, X
(I) = [x 1 ( i) x 2 (i) ··· x m (i)] T
Let T be the transpose, * be the complex conjugate, and E [] be the expected value (ensemble average), R xx = E [X * (i) X T (i)] V xr = E [X * (I) r (i)].

【0031】図3に示す複数のアンテナ素子30−1〜
30−mからの入力が、図4に示すアナログ/ディジタ
ル変換器31−1〜31−mにおいて所定のサンプリン
グ周期TS でディジタル信号に変換されたあと、各サン
プル値x1 (i)〜xm (i)は入力信号の共分散行列
計算部32aに入力され、共分散行列の時間平均Rxx
びVxrが計算される。ここで計算された共分散行列の時
間平均Rxx及びVxrは、逆行列計算部32bに入力され
る。ただし、共分散行列の時間平均の出力は、サンプリ
ング周期TS より長い周期CTS (ただし、Cは1より
大きい定数)で行われる。また、逆行列計算部32bに
おいて共分散行列の時間平均値を取り込むのは、共分散
行列計算部32aの計算結果の出力周期に等しい周期C
S においてであり、同様に逆行列Rxx -1の算出もサン
プリング周期TS より長い周期CTS で行われる。共分
散行列計算部32aで計算された相関ベクトルVxr及び
逆行列計算部32bで計算された逆行列Rxx -1は、重み
係数計算部32cに供給され、最適重み係数Wopt が算
出される。重み係数計算部32cは逆行列計算部部32
bの計算周期に等しい周期、即ちサンプリング周期TS
より長い周期CTSで計算を行って重み係数の更新を行
うが、重み係数の乗算器33−1〜33−mへの供給
は、サンプリング周期TS に等しい周期で供給する。即
ち、周期CTS内では、前の計算結果を保持している。
The plurality of antenna elements 30-1 to 30-1 shown in FIG.
After the input from 30-m is converted into a digital signal at a predetermined sampling period T S by analog / digital converters 31-1 to 31-m shown in FIG. 4, each sample value x 1 (i) to x m (i) is input to the input signal covariance matrix calculator 32a, and the time average Rxx and Vxr of the covariance matrix are calculated. The time averages R xx and V xr of the covariance matrix calculated here are input to the inverse matrix calculation unit 32b. However, the output of the time average of the covariance matrix, the sampling period T S is longer than the period CT S (however, C is a constant greater than 1) takes place in. The inverse matrix calculator 32b takes in the time average value of the covariance matrix because the period C is equal to the output cycle of the calculation result of the covariance matrix calculator 32a.
And at the T S, likewise calculates the inverse matrix R xx -1 is also performed in a long period CT S than the sampling period T S. The correlation vector V xr calculated by the covariance matrix calculation unit 32a and the inverse matrix R xx -1 calculated by the inverse matrix calculation unit 32b are supplied to the weight coefficient calculation unit 32c, and the optimum weight coefficient W opt is calculated. . The weight coefficient calculation unit 32c is an inverse matrix calculation unit 32
b, ie, the sampling period T S
Updating the weight coefficients by performing calculations in a longer period CT S, but supplied to the multiplier 33-1 to 33-m of the weighting factors, supplied with a period equal to the sampling period T S. That is, in the period CT S, holds the previous calculation results.

【0032】以上述べたように本実施形態では、共分散
行列計算部32a、逆行列計算部32b及び重み係数計
算部32cから構成される重み係数計算回路32の計算
処理をサンプリング周期TS より長い周期CTS で行
い、乗算器33−1〜33−m及び足し算器34の計算
処理をサンプリング周期TS で行っている。即ち、計算
結果の出力頻度の違う複数のプロセスをマルチタスクで
処理しているので、サンプリング周期当たりの計算量を
軽減し、結果として、同じ処理能力のプロセッサを用い
た場合に伝送速度の向上が可能となる。
As described above, in the present embodiment, the calculation processing of the weight coefficient calculation circuit 32 composed of the covariance matrix calculation section 32a, the inverse matrix calculation section 32b and the weight coefficient calculation section 32c is longer than the sampling period T S. performed in the period CT S, it is performed calculation processing of the multipliers 33-1 to 33-m and the adder 34 with a sampling period T S. That is, since a plurality of processes with different output frequencies of the calculation results are processed by multitasking, the amount of calculation per sampling cycle is reduced, and as a result, the transmission speed is improved when processors having the same processing capacity are used. It becomes possible.

【0033】図5は、本発明の他の実施形態として、適
応型アレーアンテナの制御装置を概略的に示すブロック
図である。同図において、51−1〜51−mは図示し
てないが図3のアンテナ素子30−1〜30−mと同様
のアンテナ素子からのアナログ出力をディジタル信号に
変換するアナログ/ディジタル変換器、52はアナログ
/ディジタル変換器51−1〜51−mの出力から重み
係数を算出する重み係数計算回路、53−1〜53−m
はアナログ/ディジタル変換器51−1〜51−mの出
力に重み係数を乗算する乗算器、54は乗算器53−1
〜53−mの出力を合成する足し算器をそれぞれ示して
いる。重み係数計算回路52は、到来角度計算部52a
及び重み係数推定部52bから構成されている。
FIG. 5 is a block diagram schematically showing a control device of an adaptive array antenna according to another embodiment of the present invention. In the figure, 51-1 to 51-m are analog / digital converters for converting analog outputs from antenna elements similar to the antenna elements 30-1 to 30-m in FIG. Reference numeral 52 denotes a weight coefficient calculation circuit for calculating a weight coefficient from the outputs of the analog / digital converters 51-1 to 51-m, and 53-1 to 53-m.
Is a multiplier for multiplying the output of the analog / digital converters 51-1 to 51-m by a weight coefficient, and 54 is a multiplier 53-1.
5 shows adders for combining outputs of .about.53-m. The weight coefficient calculation circuit 52 includes an arrival angle calculation unit 52a.
And a weight coefficient estimating unit 52b.

【0034】本実施形態は、複数のアンテナ素子からな
るアレーアンテナを制御するための制御アルゴリズムで
あるHA(Howells−Applebaum)アル
ゴリズムを、ディジタル信号処理を用いて実現した構成
例である。このHAアルゴリズムは、逐次的に到来角度
推定部が到来角度を推定して更新し、最急降下法などの
漸近的方法又は逆行列計算をともなう直接解法により、
重み係数を計算する重み係数推定部に供給する。ただ
し、希望波の信号到来角度が目標値として必要である。
このHAアルゴリズムは、例えば「アダプティブアンテ
ナ理論の進展と今後の展望」、小川・菊間、電子情報通
信学会論文誌、B−II Vol.J75−B−II
No.11 1992年11月に詳しく説明されてい
る。
This embodiment is an example of a configuration in which an HA (Howells-Applebaum) algorithm, which is a control algorithm for controlling an array antenna composed of a plurality of antenna elements, is implemented using digital signal processing. In this HA algorithm, the arrival angle estimator sequentially estimates and updates the arrival angle, and uses an asymptotic method such as a steepest descent method or a direct solution method with an inverse matrix calculation.
The weight coefficient is supplied to a weight coefficient estimating unit that calculates the weight coefficient. However, the signal arrival angle of the desired wave is required as the target value.
This HA algorithm is described in, for example, "Progress and Future Prospects of Adaptive Antenna Theory", Ogawa and Kikuma, IEICE Transactions, B-II Vol. J75-B-II
No. 11 Detailed description in November 1992.

【0035】HAアルゴリズムによれば、最適重み係数
opt は、次のように表わされる。 Wopt=Rxx -1* 又は Wopt=Ruu -1* ここで、希望波をXs (t)、干渉波をXi (t)、雑
音をXn (t)とし、それぞれ列ベクトルで表わすもの
とすると、 Rxx=E[{Xs (t)+Xi (t)+Xn (t)}*
{Xs (t)+Xi (t)+Xn (t)}T ] Ruu=E[{Xi (t)+Xn (t)}* {Xi (t)
+Xn (t)}T ] と表わされ、これらは信号波共分散行列である。また、
Sはステアリングベクトルといわれ、基準素子に対する
各素子の受信信号の相対位相を示し、例えば、等間隔の
1次元アレーの場合では、 S=[s12 ・・・ smTn =2(m−1)πd sinθ/λ と表わされる。ここで、dは素子間隔、λは受信波の波
長、θは到来角度である。さらに、T が転置、* が複素
共役、E[ ]が期待値(アンサンブル平均)を表わ
す。
According to the HA algorithm, the optimum weight coefficient W opt is expressed as follows. W opt = R xx −1 S * or W opt = R uu −1 S * where X s (t) is a desired wave, X i (t) is an interference wave, and X n (t) is noise. If it is represented by a column vector, R xx = E [{X s (t) + X i (t) + X n (t)} *
{X s (t) + X i (t) + X n (t)} T] R uu = E [{X i (t) + X n (t)} * {X i (t)
+ X n (t)} T ], which are signal wave covariance matrices. Also,
S is said to steering vector indicates the relative phase of the received signals of the respective elements relative to the reference element, for example, in the case of equally spaced one-dimensional array, S = [s 1 s 2 ··· s m] T s n = 2 (m-1) πd sin θ / λ. Here, d is the element interval, λ is the wavelength of the received wave, and θ is the arrival angle. Further, T represents transpose, * represents complex conjugate, and E [] represents an expected value (ensemble average).

【0036】上述したようにHAアルゴリズムには、こ
のような逆行列を含む計算を行って直接最適重みを求め
る方法と、最急降下法などを用いた漸近的解法がある
が、アレーアンテナのリアルタイム制御を行うには、直
接解法が望ましい。これは、漸近的な方法では、その収
束速度が電波環境の影響を受けるが、直接解法によれ
ば、そのような問題が避けられるからである。しかし、
逆行列計算に伴う、計算機への計算負荷が従来は問題と
なっていたが、本発明のごとく、重み係数を適応的に計
算する計算処理における少なくとも一部の処理の優先度
を、高速処理が要求される他の計算処理の優先度より低
く設定することによってこれは解消できることとなる。
As described above, the HA algorithm includes a method of directly calculating an optimum weight by performing a calculation including such an inverse matrix and an asymptotic solution method using a steepest descent method. To do this, a direct solution is desirable. This is because the convergence speed is affected by the radio wave environment in the asymptotic method, but such a problem can be avoided by the direct solution. But,
Conventionally, the calculation load on the computer involved in the inverse matrix calculation has been a problem. However, as in the present invention, the priority of at least a part of the calculation processing for adaptively calculating the weight coefficient is determined by the high-speed processing. This can be solved by setting the priority lower than the required priority of other calculation processing.

【0037】図5の実施形態において、複数のアンテナ
素子からの入力が、アナログ/ディジタル変換器51−
1〜51−mにおいて所定のサンプリング周期TS でデ
ィジタル信号に変換される。それらサンプル値x1
(i)〜xm (i)は、サンプリング周期TS 毎に到来
角度計算部52a及び重み係数推定部52bに入力され
る。到来角度計算部52aでは、サンプリング周期TS
で入力があるのに対し、このサンプリング周期TS より
長い周期CTS (ただし、Cは1より大きい定数)で到
来角度を計算しこの周期CTS で出力する。
In the embodiment shown in FIG. 5, inputs from a plurality of antenna elements are converted to analog / digital converters 51-.
From 1 to 51-m, the signal is converted into a digital signal at a predetermined sampling period T S. These sample values x 1
(I) to x m (i) are input to the arrival angle calculator 52a and the weighting coefficient estimator 52b for each sampling period T S. In the arrival angle calculator 52a, the sampling period T S
, The arrival angle is calculated with a period CT S longer than the sampling period T S (where C is a constant greater than 1), and output at this period CT S.

【0038】到来角度計算部52aの構成例が図6に示
されている。この例では、複数のアンテナ素子からなる
アレーアンテナの複数出力を用いて到来角度推定が行わ
れる。一般に、アレーアンテナの制御アルゴリズムにお
いて、信号の到来角度を用いて重み係数の適応制御を行
うアルゴリズムは少なくない。このため、信号到来角度
推定を精度よく行う種々のアルゴリズムが提案されてい
るが、その中にMUSICがある。このMUSICアル
ゴリズムは、Rxxの最小の固有値に対応する固有ベクト
ルがS* と直交することを利用して信号到来角度の推定
を行うものであり、例えば「High−Resolut
ion Techniques inSignal P
rocessing Antenna」、Y.Ogaw
a,N.Nikuma、IEICE Transact
ions on Commun. vol.E78−B
No.11 Nov.1995に詳しく説明されてい
る。
FIG. 6 shows an example of the structure of the arrival angle calculator 52a. In this example, the arrival angle estimation is performed using a plurality of outputs of an array antenna including a plurality of antenna elements. In general, in the control algorithm of the array antenna, there are many algorithms for adaptively controlling the weight coefficient using the arrival angle of a signal. For this reason, various algorithms for accurately estimating the angle of arrival of a signal have been proposed. Among them, there is MUSIC. The MUSIC algorithm, which eigenvector corresponding to the smallest eigenvalue of R xx is an estimate of the signal arrival angles by utilizing the fact that perpendicular to the S *, for example, "High-Resolut
ion Technologies inSignal P
processing Antenna "; Ogaw
a, N .; Nikuma, IEICE Transact
ions on Commun. vol. E78-B
No. 11 Nov. 1995.

【0039】図6において、複数のアンテナ素子からの
入力が、アナログ/ディジタル変換器51−1〜51−
mにおいて所定のサンプリング周期TS でディジタル信
号に変換されたあと、そのサンプル値x1 (i)〜xm
(i)はサンプリング周期TS 毎に共分散行列計算部5
2a1 に入力され、共分散行列の時間平均値Rxx又はR
uuが計算される。ここで計算された共分散行列の時間平
均は、固有値及び固有ベクトル計算部52a2 に入力さ
れる。ただし、共分散行列の時間平均の出力は、サンプ
リング周期TS より長い周期CTS (ただし、Cは1よ
り大きい定数)で行われる。また、固有値及び固有ベク
トル計算部52a2 において共分散行列の時間平均値を
取り込むのは、共分散行列計算部52a1 の計算結果の
出力周期に等しい周期CTS においてであり、同様に固
有値及び固有ベクトル計算の結果の算出もサンプリング
周期TS より長い周期CTS で行われる。固有値及び固
有ベクトル計算部52a2 で計算された固有値及び固有
ベクトルは、信号波到来角度推定部(MUSICアルゴ
リズム)52a3 に供給され、サンプリング周期TS
り長い周期CTS で到来角度が計算されて出力される。
到来角度推定部52a3 で求められた到来角度は、重み
係数推定部52bに供給され、最適重み係数Wopt が算
出される。重み係数推定部52bは到来角度推定部52
3 の計算周期に等しい周期、即ちサンプリング周期T
S より長い周期CTS で計算を行って重み係数の更新を
行うが、重み係数の乗算器53−1〜53−mへの供給
は、サンプリング周期TS に等しい周期で供給する。即
ち、周期CTS 内では、前の計算結果を保持している。
In FIG. 6, inputs from a plurality of antenna elements are analog-to-digital converters 51-1 to 51-.
m, the signal is converted into a digital signal at a predetermined sampling period T S , and the sampled values x 1 (i) to x m
(I) is a covariance matrix calculation unit 5 for each sampling period T S.
2a 1, and the time average value R xx or R of the covariance matrix
uu is calculated. Time average here is calculated covariance matrix is input to the eigenvalues and eigenvectors calculating section 52a 2. However, the output of the time average of the covariance matrix, the sampling period T S is longer than the period CT S (however, C is a constant greater than 1) takes place in. Also, capture the time-averaged value of the covariance matrix in the eigenvalue and eigenvector calculation section 52a 2 is a the period equal CT S to the output period of the covariance matrix calculating unit 52a 1 of the calculation results are also eigenvalues and eigenvectors calculated Is also calculated at a period CT S longer than the sampling period T S. Calculated eigenvalues and eigenvectors eigenvalues and eigenvectors calculating section 52a 2 is supplied to the DOA estimating unit (MUSIC algorithm) 52a 3, the arrival angle is output is calculated in a long period CT S than the sampling period T S You.
Arrival angle determined by the arrival angle estimation unit 52a 3 is supplied to the weighting coefficient estimator 52 b, the optimal weight coefficient W opt is calculated. The weighting coefficient estimating unit 52b
a period equal to the calculation period of a 3 , that is, the sampling period T
The weight coefficient is updated by performing the calculation at a cycle CT S longer than S , but the weight coefficient is supplied to the multipliers 53-1 to 53-m at a cycle equal to the sampling cycle T S. That is, in the period CT S, holds the previous calculation results.

【0040】乗算器53−1〜53−mの出力は、足し
算器54に印加され各アンテナ素子に対応する出力が合
成される。乗算器53−1〜53−m及び足し算器54
の計算は、リアルタイムの適応制御を実現するために、
所定のサンプリング周期TS以内に結果が算出されなけ
ればならない。即ち、本発明における信号の高速処理が
要求される計算処理は、本実施形態では、アナログ/デ
ィジタル変換器51−1〜51−m、乗算器53−1〜
53−m及び足し算器54の計算処理に対応している。
The outputs of the multipliers 53-1 to 53-m are applied to an adder 54, and the outputs corresponding to the respective antenna elements are synthesized. Multipliers 53-1 to 53-m and adder 54
Is calculated to achieve real-time adaptive control.
The result must be calculated within a predetermined sampling period T S. That is, in the present embodiment, the calculation processing that requires high-speed processing of signals in the present invention is performed by the analog / digital converters 51-1 to 51-m and the multipliers 53-1 to 53-1.
53-m and the calculation processing of the adder 54.

【0041】以上述べたように本実施形態では、到来角
度計算部52a及び重み係数推定部52bから構成され
る重み係数計算回路52の計算処理をサンプリング周期
Sより長い周期CTS で行い、乗算器53−1〜53
−m、足し算器54及びディジタル/アナログ変換器5
5の計算処理をサンプリング周期TS で行っている。即
ち、計算結果の出力頻度の違う複数のプロセスをマルチ
タスクで処理しているので、サンプリング周期当たりの
計算量を軽減し、結果として、同じ処理能力のプロセッ
サを用いた場合に伝送速度の向上が可能となる。
[0041] In this embodiment as described above, performs a calculation process consists weighting factor calculator circuit 52 from the incoming angle calculation unit 52a and the weight coefficient estimator 52b at the sampling period T S is longer than the period CT S, multiplication Vessels 53-1 to 53
−m, adder 54 and digital / analog converter 5
The calculation processing of No. 5 is performed at the sampling period T S. That is, since a plurality of processes with different output frequencies of the calculation results are processed by multitasking, the amount of calculation per sampling cycle is reduced, and as a result, the transmission speed is improved when processors having the same processing capacity are used. It becomes possible.

【0042】図7は、本発明のさらに他の実施形態とし
て、適応型アレーアンテナの制御装置を概略的に示すブ
ロック図である。同図において、71−1〜71−mは
図示してないが図3のアンテナ素子30−1〜30−m
と同様のアンテナ素子からのアナログ出力をディジタル
信号に変換するアナログ/ディジタル変換器、72はア
ナログ/ディジタル変換器71−1〜71−mの出力及
び足し算器74のフィードバック出力から重み係数を算
出する重み係数計算回路、73−1〜73−mはアナロ
グ/ディジタル変換器71−1〜71−mの出力に重み
係数を乗算する乗算器、74は乗算器73−1〜73−
mの出力を合成する足し算器をそれぞれ示している。重
み係数計算回路72は、到来角度計算部72a及び重み
係数推定部72bから構成されている。
FIG. 7 is a block diagram schematically showing a control device of an adaptive array antenna according to still another embodiment of the present invention. 3, 71-1 to 71-m are not shown, but the antenna elements 30-1 to 30-m of FIG.
An analog / digital converter 72 for converting an analog output from an antenna element into a digital signal in the same manner as described above, and a weight coefficient 72 is calculated from the outputs of the analog / digital converters 71-1 to 71-m and the feedback output of the adder 74. Weight coefficient calculation circuits, 73-1 to 73-m are multipliers for multiplying the outputs of analog / digital converters 71-1 to 71-m by weight coefficients, and 74 is multipliers 73-1 to 73-m.
An adder for combining the outputs of m is shown. The weight coefficient calculation circuit 72 includes an arrival angle calculation section 72a and a weight coefficient estimation section 72b.

【0043】本実施形態も図5の実施形態と同様に、複
数のアンテナ素子からなるアレーアンテナを制御するた
めの制御アルゴリズムであるHAアルゴリズムを、ディ
ジタル信号処理を用いて実現した構成例である。ただ
し、本実施形態では、重み係数を、HAループを介して
漸近的な方法で収束させている。このため、到来角度計
算部72aでは、サンプリング周期TS より長い周期C
S (ただし、Cは1より大きい定数)で到来角度を計
算しこの周期CTS で出力するが、重み係数推定部72
bでは、サンプリング周期TS で重み係数を算出して更
新しその値を出力する。本実施形態のその他の構成、ア
ルゴリズム、作用効果等は、図5の実施形態の場合と全
く同様である。本発明における信号の高速処理が要求さ
れる計算処理は、本実施形態では、アナログ/ディジタ
ル変換器71−1〜71−m、乗算器73−1〜73−
m、足し算器74及び重み係数推定部72bの計算処理
に対応している。
As in the embodiment of FIG. 5, the present embodiment is an example of a configuration in which an HA algorithm, which is a control algorithm for controlling an array antenna composed of a plurality of antenna elements, is implemented using digital signal processing. However, in the present embodiment, the weight coefficients are converged in an asymptotic manner via the HA loop. For this reason, the arrival angle calculator 72a sets the cycle C longer than the sampling cycle T S.
T S (however, C is a constant greater than 1), but calculates and outputs the arrival angle in this period CT S, weighting coefficient estimator 72
In b, the weight coefficient is calculated and updated at the sampling period T S , and the updated value is output. Other configurations, algorithms, operational effects, and the like of this embodiment are exactly the same as those of the embodiment of FIG. In the present embodiment, the calculation processing that requires high-speed processing of signals in the present invention is performed by analog / digital converters 71-1 to 71-m and multipliers 73-1 to 73-.
m, the calculation processing of the adder 74 and the weight coefficient estimation unit 72b.

【0044】以上述べた実施形態は全て本発明を例示的
に示すものであって限定的に示すものではなく、本発明
は他の種々の変形態様及び変更態様で実施することがで
きる。従って本発明の範囲は特許請求の範囲及びその均
等範囲によってのみ規定されるものである。
The embodiments described above all show the present invention by way of example and not by way of limitation, and the present invention can be embodied in other various modifications and alterations. Therefore, the scope of the present invention is defined only by the appended claims and their equivalents.

【0045】[0045]

【発明の効果】以上詳細に説明したように本発明では、
複数のアンテナ素子を有する適応型アレーアンテナの制
御をディジタル信号処理を用いて行う場合に、各アンテ
ナ素子に関する信号の重み係数を適応的に計算する計算
処理における少なくとも一部の処理の優先度を、信号の
高速処理が要求される他の計算処理の優先度より低く設
定している。このように、高速な計算結果出力が必要な
い部分の計算処理を、分割して実行することにより、サ
ンプリング周期当たりの計算量を軽減し、結果として、
同じ処理能力のプロセッサを用いた場合に伝送速度の向
上が可能となる。即ち、現状の限られたディジタル信号
処理機構の計算能力をより効率的に用いて、より高速な
伝送速度に対して適用することができる。
As described in detail above, in the present invention,
When performing the control of the adaptive array antenna having a plurality of antenna elements using digital signal processing, the priority of at least a part of the processing in the calculation processing of adaptively calculating the weight coefficient of the signal for each antenna element, The priority is set lower than the priority of other calculation processing requiring high-speed signal processing. In this way, by dividing and executing the calculation processing of the portion that does not require high-speed calculation result output, the calculation amount per sampling cycle is reduced, and as a result,
When processors having the same processing capacity are used, the transmission speed can be improved. That is, the present invention can be applied to higher transmission speeds by more efficiently using the calculation capability of the currently limited digital signal processing mechanism.

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

【図1】従来のフェーズドアレー装置を示すブロック図
である。
FIG. 1 is a block diagram showing a conventional phased array device.

【図2】本発明の基本原理を説明するための図である。FIG. 2 is a diagram for explaining a basic principle of the present invention.

【図3】本発明の一実施形態として適応型アレーアンテ
ナの制御装置を概略的に示すブロック図である。
FIG. 3 is a block diagram schematically showing a control device of an adaptive array antenna according to an embodiment of the present invention.

【図4】図3の実施形態における重み係数計算回路の構
成例を示すブロック図である。
FIG. 4 is a block diagram illustrating a configuration example of a weight coefficient calculation circuit in the embodiment of FIG. 3;

【図5】本発明の他の実施形態として適応型アレーアン
テナの制御装置を概略的に示すブロック図である。
FIG. 5 is a block diagram schematically showing a control device of an adaptive array antenna according to another embodiment of the present invention.

【図6】図5の実施形態における到来角度計算部の構成
例を示すブロック図である。
FIG. 6 is a block diagram illustrating a configuration example of an arrival angle calculator in the embodiment of FIG. 5;

【図7】本発明のさらに他の実施形態として適応型アレ
ーアンテナの制御装置を概略的に示すブロック図であ
る。
FIG. 7 is a block diagram schematically showing a control device of an adaptive array antenna according to still another embodiment of the present invention.

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

30 アレーアンテナ 30−1〜30−m アンテナ素子 31−1〜31−m、51−1〜51−m、71−1〜
71−m アナログ/ディジタル変換器 32、52、72 重み係数計算回路 32a、52a1 共分散行列計算部 32b、72b 逆行列計算部 32c 重み係数計算部 33−1〜33−m、53−1〜53−m 乗算器 34、54、74 足し算器 52a、72a 到来角度計算部 52a2 固有値及び固有ベクトル計算部 52a3 信号到来角度推定部 52b 重み係数推定部
30 array antenna 30-1 to 30-m antenna element 31-1 to 31-m, 51-1 to 51-m, 71-1
71-m analog / digital converter 32, 52, 72 weight coefficient calculation circuit 32a, 52a 1 covariance matrix calculation section 32b, 72b inverse matrix calculation section 32c weight coefficient calculation section 33-1 to 33-m, 53-1 53-m multipliers 34, 54, 74 Adders 52a, 72a Arrival angle calculator 52a 2 Eigenvalue and eigenvector calculator 52a 3 Signal arrival angle estimator 52b Weight coefficient estimator

Claims (10)

【特許請求の範囲】[Claims] 【請求項1】 複数のアンテナ素子を有する適応型アレ
ーアンテナの制御をディジタル信号処理を用いて行う方
法であって、該各アンテナ素子に関する信号の重み係数
を適応的に計算する計算処理における少なくとも一部の
処理の優先度を、前記信号の高速処理が要求される他の
計算処理の優先度より低く設定したことを特徴とする適
応型アレーアンテナの制御方法。
1. A method for controlling an adaptive array antenna having a plurality of antenna elements using digital signal processing, wherein at least one of the calculation processes for adaptively calculating a weight coefficient of a signal related to each antenna element is performed. A method of controlling an adaptive array antenna, wherein the priority of the processing of the unit is set lower than the priority of another calculation processing requiring high-speed processing of the signal.
【請求項2】 前記アンテナ素子で受信した複数のアナ
ログ信号をディジタル信号に変換し、該変換された複数
のディジタル信号に基づいて前記複数のアンテナ素子対
応の重み係数を適応的に計算し、該計算した複数の重み
係数と前記複数のディジタル信号とをそれぞれ乗算した
後に合成する処理を含む方法であって、前記高速処理が
要求される他の計算処理が、前記ディジタル信号への変
換処理、前記乗算処理及び前記合成処理を少なくとも含
んでいることを特徴とする請求項1に記載の方法。
2. Converting a plurality of analog signals received by the antenna element into digital signals, adaptively calculating weighting factors corresponding to the plurality of antenna elements based on the converted plurality of digital signals, A method comprising a process of multiplying each of the calculated plurality of weighting coefficients and the plurality of digital signals and then synthesizing each of the plurality of weighted coefficients, wherein another calculation process requiring the high-speed processing is a conversion process to the digital signal, The method of claim 1, comprising at least a multiplication operation and the combining operation.
【請求項3】 前記高速処理が要求される他の計算処理
が、前記ディジタル信号への変換処理のサンプリング周
期内に完了されることを特徴とする請求項2に記載の方
法。
3. The method according to claim 2, wherein the other calculation processing requiring the high-speed processing is completed within a sampling period of the conversion processing into the digital signal.
【請求項4】 前記重み係数の計算処理が、前記ディジ
タル信号への変換処理のサンプリング周期より長い周期
で完了されることを特徴とする請求項2又は3に記載の
方法。
4. The method according to claim 2, wherein the calculation process of the weight coefficient is completed in a cycle longer than a sampling cycle of the conversion process into the digital signal.
【請求項5】 複数のアンテナ素子を有する適応型アレ
ーアンテナの制御をディジタル信号処理を用いて行う装
置であって、該各アンテナ素子に関する信号の重み係数
を適応的に計算する計算処理における少なくとも一部の
処理の優先度を、前記信号の高速処理が要求される他の
計算処理の優先度より低く設定したことを特徴とする適
応型アレーアンテナの制御装置。
5. An apparatus for controlling an adaptive array antenna having a plurality of antenna elements by using digital signal processing, wherein at least one of the calculation processing for adaptively calculating a weight coefficient of a signal relating to each of the antenna elements is performed. A control device for an adaptive array antenna, wherein the priority of the processing of the unit is set lower than the priority of another calculation processing requiring high-speed processing of the signal.
【請求項6】 複数のアンテナ素子を有する適応型アレ
ーアンテナの制御をディジタル信号処理を用いて行う装
置であって、前記アンテナ素子で受信した複数のアナロ
グ信号をディジタル信号に変換するアナログ/ディジタ
ル変換手段と、該変換された複数のディジタル信号に基
づいて前記複数のアンテナ素子対応の重み係数を適応的
に計算する重み係数計算手段と、該重み係数計算手段か
ら出力される複数の重み係数と前記変換手段から得られ
る複数のディジタル信号とをそれぞれ乗算する乗算手段
と、該乗算手段の出力を合成する合成手段とを備えてお
り、前記重み係数計算手段における計算処理における少
なくとも一部の処理の優先度を、高速処理が要求される
前記アナログ/ディジタル変換手段、前記乗算手段及び
前記合成手段の計算処理の優先度より低く設定したこと
を特徴とする適応型アレーアンテナの制御装置。
6. An apparatus for controlling an adaptive array antenna having a plurality of antenna elements by using digital signal processing, wherein the plurality of analog signals received by the antenna elements are converted into digital signals. Means, weight coefficient calculating means for adaptively calculating weight coefficients corresponding to the plurality of antenna elements based on the converted plurality of digital signals, a plurality of weight coefficients output from the weight coefficient calculating means, A multiplying means for multiplying each of the plurality of digital signals obtained from the converting means; and a synthesizing means for synthesizing an output of the multiplying means. Calculation of the analog / digital converting means, the multiplying means and the synthesizing means for which high-speed processing is required. A control device for an adaptive array antenna, wherein the control is set lower than the processing priority.
【請求項7】 前記重み係数計算手段が、前記アナログ
/ディジタル変換手段から得られた複数のサンプル値の
共分散行列の時間平均値を、該アナログ/ディジタル変
換手段のサンプリング周期より長い周期で算出する第1
の計算手段と、該第1の計算手段から得られた共分散行
列の逆行列を、前記サンプリング周期より長い周期で算
出する第2の計算手段と、該第2の計算手段により得ら
れた信号行列の逆行列を用いてアレーアンテナの重み係
数を算出し、前記サンプリング周期より長い所定時間の
間、該重み係数を保持すると共に、該所定時間を周期と
して重み係数を更新する第3の計算手段とを含んでいる
ことを特徴とする請求項6に記載の装置。
7. The weight coefficient calculating means calculates a time average value of a covariance matrix of a plurality of sample values obtained from the analog / digital converting means at a cycle longer than a sampling cycle of the analog / digital converting means. First
Calculating means, a second calculating means for calculating an inverse matrix of the covariance matrix obtained from the first calculating means at a cycle longer than the sampling cycle, and a signal obtained by the second calculating means. A third calculating means for calculating a weighting factor of the array antenna by using an inverse matrix of the matrix, holding the weighting factor for a predetermined time longer than the sampling period, and updating the weighting factor with the predetermined time as a period 7. The apparatus of claim 6, comprising:
【請求項8】 前記重み係数計算手段が、前記アナログ
/ディジタル変換手段のサンプリング周期毎に該アナロ
グ/ディジタル変換手段から入力される複数のサンプル
値の共分散行列の時間平均値を求めて信号波到来角度を
前記サンプリング周期より長い周期で算出する第1の計
算手段と、該第1の計算手段により得られた信号波到来
角度に基づいてアレーアンテナの重み係数を算出し、前
記サンプリング周期より長い所定時間の間、該重み係数
を保持すると共に、該所定時間を周期として重み係数を
更新する第2の計算手段とを含んでいることを特徴とす
る請求項6に記載の装置。
8. The weighting coefficient calculating means obtains a time average value of a covariance matrix of a plurality of sample values input from the analog / digital conversion means for each sampling cycle of the analog / digital conversion means, and obtains a signal wave. First calculating means for calculating the arrival angle at a period longer than the sampling period, and calculating a weighting factor of the array antenna based on the signal wave arrival angle obtained by the first calculating unit; 7. The apparatus according to claim 6, further comprising: a second calculating unit that holds the weighting factor for a predetermined time and updates the weighting factor with the predetermined time as a cycle.
【請求項9】 前記第1の計算手段が、前記アナログ/
ディジタル変換手段から得られた複数のサンプル値の共
分散行列の時間平均値を、該アナログ/ディジタル変換
手段のサンプリング周期より長い周期で算出する第3の
計算手段と、該第3の計算手段から得られた共分散行列
の固有値及び固有ベクトルを、前記サンプリング周期よ
り長い周期で算出する第4の計算手段と、該第4の計算
手段により得られた信号行列の固有値及び固有ベクトル
から信号到来角度を求める第5の計算手段とを含んでい
ることを特徴とする請求項8に記載の装置。
9. The method according to claim 8, wherein the first calculating means calculates the analog /
A third calculating means for calculating a time average value of a covariance matrix of a plurality of sample values obtained from the digital converting means in a cycle longer than a sampling cycle of the analog / digital converting means; and Fourth calculating means for calculating the eigenvalues and eigenvectors of the obtained covariance matrix at a period longer than the sampling period, and obtaining the signal arrival angle from the eigenvalues and eigenvectors of the signal matrix obtained by the fourth calculating means 9. Apparatus according to claim 8, including fifth calculating means.
【請求項10】 前記重み係数計算手段が、前記アナロ
グ/ディジタル変換手段のサンプリング周期毎に該アナ
ログ/ディジタル変換手段から入力される複数のサンプ
ル値の共分散行列の時間平均値を求めて信号波到来角度
を前記サンプリング周期より長い周期で算出し保持する
第1の計算手段と、該第1の計算手段により得られた信
号波到来角度と前記合成手段からフィードバックされる
受信信号と各素子出力のサンプル値とから、アレーアン
テナの重み係数を、前記サンプリング周期で算出して更
新する第2の計算手段とを含んでいることを特徴とする
請求項6に記載の装置。
10. The signal weighting means according to claim 1, wherein said weight coefficient calculating means obtains a time average value of a covariance matrix of a plurality of sample values inputted from said analog / digital conversion means at each sampling cycle of said analog / digital conversion means. First calculating means for calculating and holding the angle of arrival at a period longer than the sampling period; a signal wave arrival angle obtained by the first calculating means, a received signal fed back from the synthesizing means, and an output of each element; 7. The apparatus according to claim 6, further comprising: second calculating means for calculating and updating a weight coefficient of the array antenna at the sampling period from the sample value.
JP04287697A 1997-02-13 1997-02-13 Method and apparatus for controlling adaptive array antenna Expired - Fee Related JP3438513B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP04287697A JP3438513B2 (en) 1997-02-13 1997-02-13 Method and apparatus for controlling adaptive array antenna

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP04287697A JP3438513B2 (en) 1997-02-13 1997-02-13 Method and apparatus for controlling adaptive array antenna

Publications (2)

Publication Number Publication Date
JPH10229306A true JPH10229306A (en) 1998-08-25
JP3438513B2 JP3438513B2 (en) 2003-08-18

Family

ID=12648253

Family Applications (1)

Application Number Title Priority Date Filing Date
JP04287697A Expired - Fee Related JP3438513B2 (en) 1997-02-13 1997-02-13 Method and apparatus for controlling adaptive array antenna

Country Status (1)

Country Link
JP (1) JP3438513B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002124900A (en) * 2000-08-03 2002-04-26 Lucent Technol Inc Method for adaptive antenna array and device to be used in multiplex connection wireless communication system

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102624431B (en) * 2012-03-13 2014-12-17 中兴通讯股份有限公司 Method and system for improving communication capacity of multiple antennae at vehicle-mounted terminal

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002124900A (en) * 2000-08-03 2002-04-26 Lucent Technol Inc Method for adaptive antenna array and device to be used in multiplex connection wireless communication system

Also Published As

Publication number Publication date
JP3438513B2 (en) 2003-08-18

Similar Documents

Publication Publication Date Title
CN108462521B (en) Anti-interference realization method of self-adaptive array antenna
EP1043801B1 (en) Adaptive array antenna system
JP4574266B2 (en) Method and apparatus for receiving digital television signals using spatial diversity and beamforming
EP0670608B1 (en) Apparatus and method for adaptively controlling array antenna comprising adaptive control means with improved initial value setting arrangement
US6738018B2 (en) All digital phased array using space/time cascaded processing
US5854612A (en) Method and apparatus for controlling array antenna
JPH11298345A (en) Cdma adaptive antenna receiver and communication system
JP2002077011A (en) Adaptive antenna receiver
US20100060523A1 (en) Adaptive array antenna apparatus and adaptive control method therefor
JP3767799B2 (en) Null direction control method and apparatus for array antenna
Tabra et al. Hybrid MVDR-LMS beamforming for massive MIMO
JP3438513B2 (en) Method and apparatus for controlling adaptive array antenna
CN114204284B (en) Anti-interference method and system for phased array antenna
JP3438527B2 (en) Signal wave arrival angle estimation device and array antenna control device
JP2000353997A (en) Adaptive array antenna device
JP4219490B2 (en) Receiving method and receiving apparatus
JPH11308130A (en) Interference wave suppressing device
JP4576742B2 (en) Transmission / reception frequency division multiplexing radio equipment
CN113395098B (en) Method and device for multi-antenna signal combination and transmitting signal shaping
JP3405156B2 (en) Array antenna control method and device
JP2000209018A (en) Array antenna control method, device therefor and recording medium recording array antenna control program
Kaur et al. Comparative study of different algorithms to implement smart antenna array-a review
Godara Adaptive postbeamformer interference canceler with improved performance in the presence of broadband directional sources
Verma et al. Digital Beamforming Using RLS-QRD Algorithm
JP3959073B2 (en) Array antenna communication device

Legal Events

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20030513

LAPS Cancellation because of no payment of annual fees