JP4173882B2 - Arrival direction estimation method - Google Patents

Description

  The present invention relates to a method for estimating the direction of an incoming wave, and more particularly to a technique effective for improving accuracy when estimating the direction of arrival using a MUSIC algorithm.

Conventionally, a large-scale system using a base station facility and a radio wave measuring vehicle has been used to detect illegal radio wave sources that adversely affect systems such as mobile phone systems.
Recently, however, radio wave measurement vehicles are not convenient in downtown areas where there are many sources, so people are directly carrying and using portable spectrum analyzers and directional antennas (mainly Yagi antennas). .
The method of detection using a spectrum analyzer and a Yagi antenna is to manually scan the Yagi antenna 360 degrees in the horizontal plane with a human hand while looking at the screen after sweeping the frequency with the spectrum analyzer, and visually observing the change on the screen. Record that there is a source in the direction of change. Move to another location and perform the same measurement. This process is repeated to narrow down the source points mainly by three-point surveying and finally identify them. Therefore, this method has a problem that a lot of time is required.
Therefore, in recent years, an apparatus for estimating the arrival direction of an incoming wave using a MUSIC (MUltiple Slgnal Classification) algorithm that is an arrival direction estimation algorithm is being used. An apparatus for estimating the direction of arrival of an incoming wave using this MUSIC algorithm can specify a transmission source with high accuracy and in a short time. (See Non-Patent Document 1 below)

As prior art documents related to the invention of the present application, there are the following.
Masayuki Nakano, Masafumi Hirono, Yukio Sato, Takehiro Miyamoto, Masatsugu Inogome, Akimichi Hamada, Koichi Ichige, Hiroyuki Arai, “Direction of Arrival System: Overview”, Shingaku Sodai, September 2005. Nobuyoshi Kikuma, "Adaptive signal processing by array antenna", Science and Technology Publishing, 1998

An apparatus that estimates the arrival direction of an incoming wave using the MUSIC algorithm described above uses a plurality of antenna elements and estimates the arrival direction using the MUSIC algorithm.
In this MUSIC algorithm, the direction of arrival is estimated by null steering for the received signal, and the result is displayed as a spectrum called a MUSIC spectrum. For the calculation of the spectrum, an ideal antenna pattern is usually used.
However, the conventional apparatus has a problem that the estimation accuracy of the arrival direction is deteriorated due to the difference in pattern between an ideal antenna pattern and an actual antenna used for measurement.
The present invention has been made to solve the problems of the prior art, and an object of the present invention is to improve the arrival direction estimation accuracy in the arrival direction estimation method using the MUSIC algorithm. It is to provide a technology that becomes possible.

Of the inventions disclosed in this application, the outline of typical ones will be briefly described as follows.
In order to solve the above-described problems, the present invention has a plurality of antenna elements, an antenna unit that receives radio waves from a transmission source, and an RF that amplifies an output signal from the antenna unit and converts it into an IF signal. And an AD conversion unit that converts the output signal from the RF unit into a digital signal, and a calculation unit that performs signal processing on the digital signal from the AD conversion unit and estimates the arrival direction of the incoming wave An arrival direction estimation method of an incoming wave in an arrival direction estimation apparatus, wherein the calculation unit calculates in advance a mode vector corresponding to an actual antenna pattern, and apparatus calibration for the received signal Step 2 for correcting the mutual coupling between the antenna elements, Step 3 for creating a correlation matrix using the signal subjected to the two corrections in Step 2, and the step Step 4 for applying the FB spatial average to the correlation matrix created in Step 3, Step 5 for obtaining the eigenvector of the matrix subjected to the FB spatial average in Step 4, and the actual antenna pattern calculated in Step 1 A step 6 for calculating a MUSIC spectrum using the mode vector and a vector corresponding to the noise subspace of the eigenvector obtained in step 5, wherein the step 1 includes apparatus calibration for the reference signal; Step 11 for correcting mutual coupling between antenna elements, Step 12 for creating a correlation matrix using the reference signal subjected to two corrections in Step 11, and FB spatial averaging for the correlation matrix created in Step 12 13 and the eigenvector of the matrix subjected to the FB space average in step 13 Look, the vector corresponding to the signal subspace of the eigen vectors, characterized by a step 14 of extracting as a mode vectors corresponding to the actual antenna patterns.

The effects obtained by the representative ones of the inventions disclosed in the present application will be briefly described as follows.
According to the present invention, it is possible to improve the direction of arrival estimation accuracy in the method of estimating the direction of arrival of incoming waves using the MUSIC algorithm.

Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
In all the drawings for explaining the embodiments, parts having the same functions are given the same reference numerals, and repeated explanation thereof is omitted.
[Example]
FIG. 1 is a block diagram showing a schematic configuration of an arrival direction estimation apparatus to which an arrival direction estimation method of an arrival wave according to an embodiment of the present invention is applied.
The arrival direction estimation device shown in FIG. 1A has a plurality of antenna elements, receives an radio wave from a transmission source, amplifies an output signal from the antenna unit 10, and converts it to an IF signal. An RF unit 11, an AD conversion unit 12 that converts an output signal from the RF unit 11 into a digital signal, and a personal computer that performs signal processing on the digital signal from the AD conversion unit 12 and estimates an arrival direction of an incoming wave 13. Here, the personal computer 13 constitutes a calculation unit.
In FIG. 1A, the antenna unit 10, the RF unit 11, and the AD conversion unit 12 within the broken line can be integrally formed.
As illustrated in FIG. 2, the antenna unit 10 illustrated in FIG. 1A includes four antenna elements (R1 to R4).
The RF unit 11 amplifies output signals output from the four antenna elements (R1 to R4), and converts them into 4-channel IF signals. The four-channel IF signal is converted into a digital signal by the AD conversion unit 12 and then converted into a digital baseband signal by the calculation unit of the personal computer 13.
When the antenna unit 10 has four or more antenna elements, as shown in FIG. 1B, the high frequency switch 15 is disposed between the antenna unit 10 and the RF unit 11, and the high frequency switch 15 The output signals from the four antenna elements may be selected.

Prior to describing the arrival direction estimation method of the arrival wave of the present embodiment, a conventional arrival direction estimation method of the arrival wave will be described with reference to FIG. FIG. 9 is a flowchart for explaining the procedure of a conventional arrival direction estimation method for incoming waves.
In the conventional arrival direction estimation method of the incoming wave, the personal computer 13 estimates the arrival direction of the incoming wave according to the following procedure.
(1) Sampling data of received signals of four channels of incoming waves received by the four antenna elements (R1 to R4) of the antenna unit 10 is acquired. (Step 200 in FIG. 9)
(2) Device system calibration is performed on the acquired sampling data. (Step 201 in FIG. 9)
(3) Correction of mutual coupling between antenna elements is performed on the acquired sampling data. (Step 202 in FIG. 9)
(4) Create a correlation matrix of sampling data subjected to device system calibration and correction of mutual coupling between antenna elements. (Step 203 in FIG. 9)
(5) FB spatial averaging is performed on the obtained correlation matrix. (Step 204 in FIG. 9)
(6) The eigenvector of the matrix subjected to FB space averaging is obtained. (Step 205 in FIG. 9)
(5) The MUSIC spectrum is calculated using the mode vector corresponding to the ideal (omni) antenna pattern and the vector corresponding to the noise subspace of the eigenvector obtained in the above step. (Step 206 in FIG. 9)
However, in the conventional method for estimating the arrival direction of an incoming wave, the arrival direction estimation accuracy deteriorates due to the difference between the ideal antenna pattern and the actual antenna used in the measurement.

Next, the arrival direction estimation method of the arrival wave of the present embodiment will be described with reference to FIG. FIG. 3A is a flowchart for explaining the procedure of the arrival direction estimation method of the incoming wave according to the embodiment of the present invention.
In the arrival direction estimation method of the arrival wave of this embodiment, the personal computer 13 estimates the arrival direction of the arrival wave by the following procedure.
(1) Sampling data of received signals of four channels of incoming waves received by the four antenna elements (R1 to R4) of the antenna unit 10 is acquired. (Step 100 in FIG. 3A)
(2) Device system calibration is performed on the acquired sampling data. (Step 101 in FIG. 3A)
(3) Correction of mutual coupling between antenna elements is performed on the acquired sampling data. (Step 102 in FIG. 3A)
(4) Create a correlation matrix of sampling data subjected to device system calibration and correction of mutual coupling between antenna elements. (Step 103 in FIG. 3A)
(5) FB spatial averaging is performed on the obtained correlation matrix. (Step 104 in FIG. 3A)
(6) The eigenvector of the matrix subjected to FB space averaging is obtained. (Step 105 in FIG. 3A)
(7) The MUSIC spectrum is calculated using the mode vector corresponding to the actual antenna pattern calculated in advance and the vector corresponding to the noise subspace of the eigenvector obtained in the above step. (Step 106 in FIG. 3A)
In this embodiment, device system calibration, correction of mutual coupling between antenna elements, and a mode vector corresponding to an actual antenna pattern are used in order to improve the direction of arrival estimation accuracy.

Each item will be described below.
(A) Device system calibration, i.e., channel-to-channel phase difference due to cable manufacturing error between antenna unit 10 and RF unit 11, channel-to-channel phase within RF unit 11, channel-to-channel variation due to amplitude difference, etc. The phase / amplitude difference is corrected by the following method.
A reference signal having the same amplitude and phase is input to each channel, and the phase and amplitude from the reference channel are obtained based on the signal. By applying the inverse characteristics of the obtained phase difference and amplitude difference to the received signal of each channel, the phase difference and amplitude difference are corrected.
FIG. 4 is a graph showing an example of the phase difference between the channels measured by the arrival direction estimation apparatus shown in FIG. In FIG. 4 and FIGS. 5 and 6 described later, the reference signal is input to the antenna unit 10 with the arrival direction (θ) shifted by 0.5 degrees and the phase difference from the reference channel is measured. The solid line indicates the theoretical value, and the broken line indicates the actual measurement value.
Moreover, in FIG. 4 and FIGS. 5 and 6 described later, the vertical axis indicates the phase difference (φ), and the horizontal axis indicates the plotted points plotted by 0.5 degrees. That is, the horizontal axis represents the angle (θ) of the arrival direction.
Assuming that the inter-channel phase difference and the amplitude difference are values shown in the following table (1), FIG. 5 shows the result of correcting the inter-channel phase / amplitude difference for the received signal of each channel shown in FIG. Show.
Comparing FIG. 4 and FIG. 5, it can be seen that the actual measurement value is closer to the theoretical value.

(B) Correction of mutual coupling between antenna elements is performed by the following method.
In general, mutual coupling between antenna elements refers to a phenomenon in which, when antenna elements are arranged close to each other, unlike when one antenna element is arranged, the antenna pattern is distorted and varies due to the influence of the antenna elements. .
The correction of mutual coupling between antenna elements is performed by obtaining the mutual coupling between antenna elements using the result of measurement with a single transmission source, and giving the inverse characteristics of the mutual coupling between antenna elements to the received signal of each channel. .
FIG. 6 shows the result of correcting the mutual coupling between the antenna elements for the reception signals of the respective channels shown in FIG. When FIG. 5 is compared with FIG. 6, it can be seen that the actually measured value is closer to the theoretical value, but there are still places where the actually measured value is deviated from the theoretical value. This shift causes the arrival direction estimation accuracy to deteriorate.

(C) The extraction of the mode vector corresponding to the actual antenna pattern is executed according to the procedure shown in FIG.
FIG. 3B is a flowchart showing a procedure for extracting a mode vector corresponding to an actual antenna pattern in the arrival direction estimation method for an incoming wave according to the embodiment of the present invention.
(1) Sampling data of reference signals of four channels of reference waves received by the four antenna elements (R1 to R4) of the antenna unit 10 are acquired. (Step 110 in FIG. 3B)
(2) Device system calibration is performed on the acquired sampling data. (Step 111 in FIG. 3B)
(3) Correction of mutual coupling between antenna elements is performed on the acquired sampling data. (Step 112 in FIG. 3B)
(4) Create a correlation matrix of sampling data subjected to device system calibration and correction of mutual coupling between antenna elements. (Step 113 in FIG. 3B)
(5) FB spatial averaging is performed on the obtained correlation matrix. (Step 114 in FIG. 3B)
(6) The eigenvector of the matrix subjected to FB space averaging is obtained. (Step 115 in FIG. 3B)
(7) A vector corresponding to the signal subspace of the obtained eigenvector is extracted as a mode vector corresponding to the actual antenna pattern. (Step 116 in FIG. 3B)
Note that mode vectors other than the reference wave are obtained by interpolation.

・・・・・・・・・・・・・・ （３）

(3)

FIG. 7 is a graph showing an example of the MUSIC spectrum obtained from this example.
In addition, FIG. 8 performs two corrections, ie, apparatus calibration and correction of mutual coupling between antenna elements, and uses a mode vector corresponding to an ideal antenna pattern and a vector corresponding to a noise subspace of an eigenvector. It is a graph which shows an example of the calculated | required MUSIC spectrum.
When FIG. 7 is compared with FIG. 8, it can be seen that in this embodiment, the dynamic range (difference between the incoming wave and the noise portion) is improved by 20 dB or more.
Further, the estimation error of the arrival direction estimated by the conventional method is 2 degrees, but the estimation error of the arrival direction estimated by the method of this embodiment is 0.5 degrees.
Thus, in this embodiment, it is possible to improve the estimation accuracy of the arrival direction. The same effect was also obtained in the case of two correlated waves.
As mentioned above, the invention made by the present inventor has been specifically described based on the above embodiments. However, the present invention is not limited to the above embodiments, and various modifications can be made without departing from the scope of the invention. Of course.

It is a block diagram which shows schematic structure of the arrival direction estimation apparatus with which the arrival direction estimation method of the arrival wave of the Example of this invention is applied. It is a figure which shows arrangement | positioning of the antenna element of the antenna part shown to Fig.1 (a). It is a flowchart for demonstrating the procedure of the arrival direction estimation method of the incoming wave of the Example of this invention. It is a flowchart which shows the procedure for extracting the mode vector corresponding to an actual antenna pattern in the arrival direction estimation method of the incoming wave of the Example of this invention. It is a graph which shows an example of the phase difference between each channel measured with the arrival direction estimation apparatus shown in FIG. It is a graph which shows the result of having corrected the phase and amplitude difference between channels to the received signal of each channel shown in FIG. It is a graph which shows the result of having corrected the mutual coupling between antenna elements to the received signal of each channel shown in FIG. It is a graph which shows an example of the MUSIC spectrum calculated | required by the arrival direction estimation method of the incoming wave of the Example of this invention. It is a graph which shows an example of the MUSIC spectrum calculated | required by the arrival direction estimation method of the conventional incoming wave. It is a flowchart for demonstrating the procedure of the conventional arrival direction estimation method of an incoming wave.

Explanation of symbols

DESCRIPTION OF SYMBOLS 10 Antenna part 11 RF part 12 AD conversion part 13 Personal computer 15 High frequency switch R1-R4 Antenna element

Claims (1)

An antenna unit having a plurality of antenna elements and receiving radio waves from a transmission source;
An RF unit for amplifying an output signal from the antenna unit and converting it to an IF signal;
An AD converter that converts an output signal from the RF unit into a digital signal;
An arrival direction estimation method for an arrival wave in an arrival direction estimation device configured by a signal processing unit for performing digital signal processing on the digital signal from the AD conversion unit and estimating an arrival direction of the arrival wave,
The calculation unit calculates in advance a mode vector corresponding to an actual antenna pattern;
Step 2 for performing device system calibration and correction of mutual coupling between antenna elements on the received signal;
Creating a correlation matrix using the signal subjected to the two corrections in step 2, and
Step 4 for performing FB spatial averaging on the correlation matrix created in Step 3;
Obtaining an eigenvector of the matrix subjected to the FB space averaging in the step 4, and
Calculating a MUSIC spectrum using the mode vector corresponding to the actual antenna pattern calculated in step 1 and the vector corresponding to the noise subspace of the eigenvector determined in step 5;
The step 1 includes a step 11 of performing device system calibration and correction of mutual coupling between antenna elements on the reference signal;
Creating a correlation matrix using the reference signal subjected to the two corrections in step 11, and
Applying an FB spatial average to the correlation matrix created in step 12;
And a step 14 of obtaining an eigenvector of the matrix subjected to the FB space averaging in the step 13 and extracting a vector corresponding to the signal subspace of the eigenvector as a mode vector corresponding to an actual antenna pattern. Arrival direction estimation method for incoming waves.

Images