CA2273530C - Radio environment analysis apparatus - Google Patents
Radio environment analysis apparatus Download PDFInfo
- Publication number
- CA2273530C CA2273530C CA002273530A CA2273530A CA2273530C CA 2273530 C CA2273530 C CA 2273530C CA 002273530 A CA002273530 A CA 002273530A CA 2273530 A CA2273530 A CA 2273530A CA 2273530 C CA2273530 C CA 2273530C
- Authority
- CA
- Canada
- Prior art keywords
- signal
- incoming
- signals
- receiver
- array antenna
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/74—Multi-channel systems specially adapted for direction-finding, i.e. having a single antenna system capable of giving simultaneous indications of the directions of different signals
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01Q—ANTENNAS, i.e. RADIO AERIALS
- H01Q3/00—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system
- H01Q3/26—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture
Abstract
A radio environment analysis apparatus includes a receiver and signal processor. The receiver receives a signal modulated by phase shift keying (PSK) using a pseudo noise (PN) code sequence by an array antenna at a plurality of fixed points on the plane, converts the received signals into intermediate frequencies or demodulated signals, and outputs them.
The signal processor processes a signal output from the receiver to analyze the arrival angle, delay time, and relative power of the signal incoming to the receiver.
The signal processor has a plurality of normalization units, first estimation unit, arithmetic unit, and second estimation unit. The normalization units individually normalize signals from the receiver and output them as a normalized signal group. The first estimation unit calculates a covariance matrix, and estimates the arrival angle of each incoming signal.
The arithmetic unit calculates a weight of the array antenna for suppressing incoming signals except for an incoming signal having the estimated arrival angle, and calculates the product of the weight of the array antenna and the normalized signal group. The second estimation unit calculates a covariance matrix and estimates the delay time and relative power.
The signal processor processes a signal output from the receiver to analyze the arrival angle, delay time, and relative power of the signal incoming to the receiver.
The signal processor has a plurality of normalization units, first estimation unit, arithmetic unit, and second estimation unit. The normalization units individually normalize signals from the receiver and output them as a normalized signal group. The first estimation unit calculates a covariance matrix, and estimates the arrival angle of each incoming signal.
The arithmetic unit calculates a weight of the array antenna for suppressing incoming signals except for an incoming signal having the estimated arrival angle, and calculates the product of the weight of the array antenna and the normalized signal group. The second estimation unit calculates a covariance matrix and estimates the delay time and relative power.
Description
~~34~b_h Specification Title of the Invention Radio Environment Analysis Apparatus Rac- grn, nd of the Invention The present invention relates to a radio environment analysis apparatus and, more particularly, to a radio environment analysis apparatus for measuring the arrival angle, delay time, and relative power of a signal incoming to the antenna in a multipath environment such as an urban or indoor space.
To realize high-speed digital communication, frequency resources must be effectively used. Many studies have been made for techniques such as a method using the milliwave band where a wide band is essentially ensured from a specific band, a method of reducing a repetitively used frequency with a decrease in cell area, and dynamic zone control effectively using adjacent cell base stations in accordance with the number of users.
As the communication speed increases, interference between multipath codes occurs to decrease the communication quality in an environment such as an urban or indoor space. LVhen the cell area decreases, interference in the same channel occurs.
One of effective means for removing interference between multipath codes and interference in the same channel is to introduce an adaptive antenna effective for removing a long delayed wave and an adaptive equalizer effective for removing a short delayed wave. In this case, to effectively operate the adaptive antenna and adaptive equalizer, the arrival angle and delay time of a signal incoming to the antenna must be measured. Also in dynamic zone control, particularly the arrival angle of a signal incoming to the antenna must be measured for beam control of the antenna.
Known examples of an algorithm of measuring the arrival angle, delay time, and relative power of a signal incoming to the antenna in the multipath environment are a 2D-MUSIC method, FFT-MUSIC method, and 2D-Unitary ESPRIT method.
Evaluation of the indoor propagation environment by the 2D-MUSIC method is reported in "High-Resolution Analysis of Indoor Multipath Propagation Structure", IEICE TRANS. COMMUN., VOL. E78-B, NO. 11, pp. 1450 - 1457, NOVEMBER 1995 (reference 1).
The principle of the FFT-MUSIC method is described in "Estimation of Propagation Delay Time and Direction of Arrival of Indoor Quasi-Millimeter Multipath Waves Using FFT-MUSIC with Triangular Antenna Array", TECHNICAL REPORT OF IEICE, A~P95-120, pp. 79 -84, February 1996 (reference 2).
Evaluation of propagation characteristics in the urban space using the 2D-Unitary ESPRIT method is described in "High-Resolution 3-D Direction-of-Arrival Determination for Urban Mobile Radio", IEEE TRANSACTIONS
ON ANTENNAS AND PROPAGATION, VOL. 45, N0. 4, APRIL 1997 (reference 3).
The principle of simultaneously estimating a combination of the arrival angle and delay time of a signal incoming to the antenna is described in "Paired Estimation of Propagation Delay Times and Directions of Arrival of Multipath Waves Using 2D Unitary ESPRIT", TECHNICAL REPORT OF IEICE, A~P97-78, JULY 1997 (reference 4).
However, these conventional methods of estimating the arrival angle and delay time of a signal incoming to the antenna suffer the following problems.
First, the 2D-MUSIC method in reference 1 and the 2D-Unitary ESPRIT method in reference 4 can estimate only the azimuth plane of the arrival angle and cannot estimate the elevation plane. At present, since a mobile telephone base station in operation controls the elevation plane directivity of the antenna in accordance with the situation, the azimuth and elevation planes of the arrival angle must be simultaneously grasped.
Second, the position of a reception antenna element is mechanically moved for measurement. This is also the same in the 2D-Unitary ESPRIT method in reference 3. The mechanical movement of the reception antenna increases the measurement time. Considering mobile communication, the antenna and equalizer must be controlled in real time, so these methods are difficult to apply to mobile communication.
Third, since a reception signal is normalized using a transmission signal as a reference, the transmission signal must be supplied to the receiver via a cable. Connecting communication stations by the cable essentially contradicts mobile communication.
To the contrary, the FFT-MUSIC method in reference 2 can simultaneously estimate the arrival angles of azimuth and elevation planes and the delay time without requiring any mechanical movement of the reception antenna and any connection cable between communication stations. However, this method does not resolve the arrival angle by a superresolution technique, resulting in low angle resolution and low precision. In reference 3, although the arrival angles of elevation and azimuth planes are estimated by the 2D-Unitary ESPRIT method, the delay time is estimated using a conventional channel counter, resulting in low time resolution and low precision.
smmmary of the Invention It is an object of the present invention to provide a radio environment analysis apparatus capable of simultaneously estimating a combination of the arrival angles of elevation and azimuth planes and the delay time with high resolution and high precision.
It is another object of the present invention to provide a radio environment measurement apparatus applicable to even mobile communication.
To achieve the above objects, according to the present invention, there is provided a radio environment analysis apparatus comprising a receiver for receiving a signal modulated by phase shift keying (PSK) using a pseudo noise (PN) code sequence by an array antenna at a plurality of fixed points on a plane, converting the plurality of received signals into intermediate frequencies or demodulated signals, and outputting the resultant signals, and a signal processor for processing a signal output from the receiver to analyze an arrival angle, delay time, and relative power of the signal incoming to the receiver, the signal processor having a plurality of normalization means for individually normalizing respective signals from the receiver and outputting the normalized signals as a normalized signal group, first estimation means for calculating a covariance matrix from the signal group normalized by the normalization means, and estimating an arrival angle of each incoming signal using one of a 2D-MUSIC method and 2D-Unitary ESPRIT method, arithmetic means for calculating a weight of the array antenna for - suppressing incoming signals except for an incoming signal having the arrival angle estimated by the first estimation means, and calculating a product of the calculated weight of the array antenna and the normalized signal group, and second estimation means for calculating a covariance matrix from a product signal group output from the arithmetic means, and estimating a delay time and relative power using one of the MUSIC
method and ESPRIT method.
Fig. 1 is a block diagram showing a radio environment analysis apparatus according to the first embodiment of the present invention;
Fig. 2 is a block diagram showing signal processing operation of a DSP shown in Fig. 1;
Fig. 3 is an explanatory view showing moving average calculation processing by the 2D-Unitary ESPRIT
method;
Fig. 4 is an explanatory view showing a moving average calculation method using the MUSIC method;
Fig. 5 is a block diagram showing a radio environment analysis apparatus according to the second embodiment of the present invention;
Fig. 6 is a block diagram showing signal processing operation of a DSP shown in Fig. 5; and Fig. 7 is a block diagram showing a radio environment analysis apparatus according to the third embodiment of the present invention.
To realize high-speed digital communication, frequency resources must be effectively used. Many studies have been made for techniques such as a method using the milliwave band where a wide band is essentially ensured from a specific band, a method of reducing a repetitively used frequency with a decrease in cell area, and dynamic zone control effectively using adjacent cell base stations in accordance with the number of users.
As the communication speed increases, interference between multipath codes occurs to decrease the communication quality in an environment such as an urban or indoor space. LVhen the cell area decreases, interference in the same channel occurs.
One of effective means for removing interference between multipath codes and interference in the same channel is to introduce an adaptive antenna effective for removing a long delayed wave and an adaptive equalizer effective for removing a short delayed wave. In this case, to effectively operate the adaptive antenna and adaptive equalizer, the arrival angle and delay time of a signal incoming to the antenna must be measured. Also in dynamic zone control, particularly the arrival angle of a signal incoming to the antenna must be measured for beam control of the antenna.
Known examples of an algorithm of measuring the arrival angle, delay time, and relative power of a signal incoming to the antenna in the multipath environment are a 2D-MUSIC method, FFT-MUSIC method, and 2D-Unitary ESPRIT method.
Evaluation of the indoor propagation environment by the 2D-MUSIC method is reported in "High-Resolution Analysis of Indoor Multipath Propagation Structure", IEICE TRANS. COMMUN., VOL. E78-B, NO. 11, pp. 1450 - 1457, NOVEMBER 1995 (reference 1).
The principle of the FFT-MUSIC method is described in "Estimation of Propagation Delay Time and Direction of Arrival of Indoor Quasi-Millimeter Multipath Waves Using FFT-MUSIC with Triangular Antenna Array", TECHNICAL REPORT OF IEICE, A~P95-120, pp. 79 -84, February 1996 (reference 2).
Evaluation of propagation characteristics in the urban space using the 2D-Unitary ESPRIT method is described in "High-Resolution 3-D Direction-of-Arrival Determination for Urban Mobile Radio", IEEE TRANSACTIONS
ON ANTENNAS AND PROPAGATION, VOL. 45, N0. 4, APRIL 1997 (reference 3).
The principle of simultaneously estimating a combination of the arrival angle and delay time of a signal incoming to the antenna is described in "Paired Estimation of Propagation Delay Times and Directions of Arrival of Multipath Waves Using 2D Unitary ESPRIT", TECHNICAL REPORT OF IEICE, A~P97-78, JULY 1997 (reference 4).
However, these conventional methods of estimating the arrival angle and delay time of a signal incoming to the antenna suffer the following problems.
First, the 2D-MUSIC method in reference 1 and the 2D-Unitary ESPRIT method in reference 4 can estimate only the azimuth plane of the arrival angle and cannot estimate the elevation plane. At present, since a mobile telephone base station in operation controls the elevation plane directivity of the antenna in accordance with the situation, the azimuth and elevation planes of the arrival angle must be simultaneously grasped.
Second, the position of a reception antenna element is mechanically moved for measurement. This is also the same in the 2D-Unitary ESPRIT method in reference 3. The mechanical movement of the reception antenna increases the measurement time. Considering mobile communication, the antenna and equalizer must be controlled in real time, so these methods are difficult to apply to mobile communication.
Third, since a reception signal is normalized using a transmission signal as a reference, the transmission signal must be supplied to the receiver via a cable. Connecting communication stations by the cable essentially contradicts mobile communication.
To the contrary, the FFT-MUSIC method in reference 2 can simultaneously estimate the arrival angles of azimuth and elevation planes and the delay time without requiring any mechanical movement of the reception antenna and any connection cable between communication stations. However, this method does not resolve the arrival angle by a superresolution technique, resulting in low angle resolution and low precision. In reference 3, although the arrival angles of elevation and azimuth planes are estimated by the 2D-Unitary ESPRIT method, the delay time is estimated using a conventional channel counter, resulting in low time resolution and low precision.
smmmary of the Invention It is an object of the present invention to provide a radio environment analysis apparatus capable of simultaneously estimating a combination of the arrival angles of elevation and azimuth planes and the delay time with high resolution and high precision.
It is another object of the present invention to provide a radio environment measurement apparatus applicable to even mobile communication.
To achieve the above objects, according to the present invention, there is provided a radio environment analysis apparatus comprising a receiver for receiving a signal modulated by phase shift keying (PSK) using a pseudo noise (PN) code sequence by an array antenna at a plurality of fixed points on a plane, converting the plurality of received signals into intermediate frequencies or demodulated signals, and outputting the resultant signals, and a signal processor for processing a signal output from the receiver to analyze an arrival angle, delay time, and relative power of the signal incoming to the receiver, the signal processor having a plurality of normalization means for individually normalizing respective signals from the receiver and outputting the normalized signals as a normalized signal group, first estimation means for calculating a covariance matrix from the signal group normalized by the normalization means, and estimating an arrival angle of each incoming signal using one of a 2D-MUSIC method and 2D-Unitary ESPRIT method, arithmetic means for calculating a weight of the array antenna for - suppressing incoming signals except for an incoming signal having the arrival angle estimated by the first estimation means, and calculating a product of the calculated weight of the array antenna and the normalized signal group, and second estimation means for calculating a covariance matrix from a product signal group output from the arithmetic means, and estimating a delay time and relative power using one of the MUSIC
method and ESPRIT method.
Fig. 1 is a block diagram showing a radio environment analysis apparatus according to the first embodiment of the present invention;
Fig. 2 is a block diagram showing signal processing operation of a DSP shown in Fig. 1;
Fig. 3 is an explanatory view showing moving average calculation processing by the 2D-Unitary ESPRIT
method;
Fig. 4 is an explanatory view showing a moving average calculation method using the MUSIC method;
Fig. 5 is a block diagram showing a radio environment analysis apparatus according to the second embodiment of the present invention;
Fig. 6 is a block diagram showing signal processing operation of a DSP shown in Fig. 5; and Fig. 7 is a block diagram showing a radio environment analysis apparatus according to the third embodiment of the present invention.
The present invention will be described in detail below with reference to the accompanying drawings.
Fig. 1 shows a radio environment analysis apparatus according to the first embodiment of the present invention. In Fig. 1, the radio environment analysis apparatus comprises a transmitter 1 for transmitting a signal modulated by PSK (Phase Shift Keying) using a PN (Pseudo Noise) code sequence, a receiver 2 for outputting signals received by a plurality of antenna elements arrayed on the plane after conversion into intermediate frequencies or demodulated signals, and a signal processor 3 for simultaneously estimating a combination of the arrival angles (elevation and azimuth angles), delay time, and relative power of a reception signal on the basis of an output from the receiver 2.
The transmitter 1 is constituted by a PN code generator 11 for generating a PN code, a BPSK (Binary PSK) modulator 12 for phase-modulating the PN code output from the PN code generator 11, a frequency converter/amplifier 13 for converting the modulated signal output from the BPSK modulator 12 into a predetermined measurement frequency, and amplifying and outputting it, and a transmission antenna 14 for radiating the signal output from the frequency converter/amplifier 13 into the space.
The receiver 2 is constituted by a planar-array antenna 15 having M X N (M and N are integers equal to or larger than 2) antenna elements 1611 to 16MN arrayed on the plane, and a plurality of frequency converters 1711 to 17MN for outputting signals received by the antenna elements 1611 and l6MNafter conversion into intermediate frequencies or basebands.
The signal processor 3 is constituted by a plurality of A/D (Analog-to-Digital) converters (A-D) 1811 to 18MN for digitally converting the frequency-converted reception signals from the frequency converters 1711 to 17MN, and a digital signal processor (DSP) 19 for simultaneously estimating a combination of the arrival angles (elevation and azimuth angles), delay time, and relative power of a signal incoming to the planar-array antenna 15 on the basis of the plurality of digital reception signals from the A/D converters 1811 to 18MN and a criterion signal.
Signal processing operation of the DSP 19 will be explained with reference to Fig. 2.
In Fig. 2, a plurality of signals that are received by the antenna elements 1611 to 16MN of the planar-array antenna 15 and digitally converted are transformed by fast Fourier transformers 10111 to 101MN.
Then, the time-series signals are converted into signals in the frequency range. To resolve a plurality of correlative signals such as multipath signals for the - g _ delay time, a moving average must be calculated in the frequency range. In calculating the moving average, the spectrum must be flat in a target frequency range.
The group delay devices of the transmitter and receiver must be excluded from processing targets for waiting for reception.
For this purpose, in an FFT 103, a transmission signal is directly input to each channel of the receiver and subjected to fast Fourier transform.
In normalization units 10211 to 102MN, reception signals converted into the frequency range are normalized by a calibration (criterion) signal from the FFT 103 that is converted into the frequency range.
The arrival angles (elevation and azimuth angles) can be estimated from a reception signal having a frequency unique to the planar-array antenna. In a covariance matrix estimation unit 104, a moving-averaged covariance matrix is estimated using the reception signals normalized to the center frequency of the band from the normalization units 10211 to 102MN.
Fig. 3 shows moving average calculation processing by the 2D-Unitary ESPRIT method.
In Fig. 3, a complex signal received by the planar-array antenna is divided into (M-m+1) x (N-n+1), m x n minor matrices Xpq. In general, m is (M+1)/2, and n is (N+1)/2. Each divided reception signal matrix Xpq is vector-transformed as given by - g -xPq = vec (XPq) . . . ( 1 ) to calculate a correlation matrix R defined by __ H
RP4 xP9x P9 . . . ( 2 ) This operation is done for each divided minor matrix XPq to average a covariance matrix RPq by M-m+1 N-n+1 R (M - m + 1) (N - n + 1) R Pq . . . ( 3 ) P_-1 9-1 The averaged covariance matrix R is unitarily converted by QHmn - RQmn . . . ( 4 ) Where QHmn QHm ~ QHn Qmn - Qm ~ Qn Qm: mth-order unitary matrix Qn: nth-order unitary matrix ~: Kronecker's operator That is, if m(n) is an even number (m(n) - 2k:
k is an integer), Ik J Ik . ( 4A) Q2k - ' ' IIk - j IIk If m(n) is an odd number (m(n) - 2k+1), Ik 0 J Ik Q2k+1 1 ~ ~T . . .
IIk 0 j IIk where 0 is a 0 vector, Ik is the kth-order unit matrix, and IIk is defined by IIk - ' To increase resolution and precision, an ensemble average is calculated in a real matrix generation unit 105 by Ry = Re [ E [ QHmnRQmn ] ] . . . ( 5 ) where E[] is the ensemble average.
The real part of the covariance matrix is extracted.
An arrival angle estimation unit 106 estimates a combination of the arrival elevation and azimuth angles by the following procedure.
First, Ry is decomposed into eigenvalues, and the obtained eigenvalues ~.i (i = 11 to mn) are arranged in the decreasing order to be expressed by where az is noise power.
From equation (6), the number L of arrival signals can be estimated from the distribution of the eigenvalues.
Then, a partial signal space Es to be defined is calculated by ES = [elez...eL] ... (7) where ei is an eigenvector corresponding to the eigenvalue Matrices E" and E" to be defined are calculated by E _ K~,~ Es , E _ K"' Es . . . ( 8 ) N v K~z Es K~z Es where Ku 1 = IM ~ Kl, Ku z = IM ~ Kz Kvi = K3 ~ IN, K"z = K9 ~ IN
Kl = Re [ QHN_1 J2QN ] , Kz = Im [ QHN_1 JzQN ]
K3 Re [QHM_lJzQM] , Kq = Im[QHM_1JZQM]
IM: Mth-order unit matrix IN: Nth-order unit matrix Jz: selection matrix for extracting the Mth or Nth order from the second row of the matrix A matrix EXY to be defined is calculated by EXY = [KuIES~K~zEs] ... (9) Further, a matrix:
lO EXYHEXY - (KN~ES)x [Ka~Es ~ K~,zEs] . . . (10) (K,,zEs) is calculated and decomposed into eigenvalues.
A matrix E is calculated by arranging the obtained eigenvectors ei (i = 1 to 2L) by E = [elez...ezL] ... (11) The matrix E is divided into four L x L matrices by E - En E~z . . . ( 12 ) Ez~ Ez2 A matrix '~If a is calculated by . (13) ~ a = -Eiz [ Ezz ] 1 . .
Subsequently, this processing using equations (9) to (13) is repetitively executed for E" to define a matrix 'iIr" corresponding to the matrix '~',~ . 1IJ'~ + j'~I~'" is decomposed into eigenvalues. Letting the obtained eigenvalue be wu(i) + jw"(i) (i = 1 to L), an arrival elevation angle 6i and azimuth angle chi can be obtained by solving simultaneous equations:
(,JN(1) _ ~ OX COS 81 COS l~l ~ (i) _ ~ Dy cos B; sin r~l . . . ( 14 ) where O x and O y are array element intervals and ~, is the wavelength of a transmission signal.
In this way, the arrival azimuth and elevation angles of L arrival signals can be obtained.
After that, the delay time and relative power are obtained.
An arrival signal designation unit 107 designates arbitrary one of the L detected arrival signals. An adaptive weight calculation unit 108 stores the designated signal to calculate the weight of the adaptive antenna for suppressing other unnecessary signals.
Letting S" be the covariance matrix and a(c~, 8) be a steering vector indicating the direction of the stored signal when an unnecessary signal is received by the array antenna, the weight of the adaptive antenna is calculated by W = a Sua* ( ~~ a ) . . . (15) An array calculation unit 109 calculates the product of the array reception data and weight by y = W"X . . . ( 16 ) At this time, the product is calculated for a complex reception signal having a plurality of frequencies normalized in the frequency range. If J
data are used in the frequency range, array outputs y(1) to y(J) are obtained. Using the array output data y(1) to y(J), an estimation moving averaging unit 110 calculates a correlation matrix S by S = YY" ... (17) where Y = [Y(1)...y(J)1T
To resolve correlative signals for the delay time, moving averaging shown in Fig. 4 is done for the calculated correlation matrix. Fig. 4 shows a moving average calculation method using the MUSIC method for J
- 8. The correlation matrix S after moving average calculation is decomposed into eigenvalues, and the obtained eigenvalues are arranged in the decreasing order to be expressed by ~,, ? ~,z ~ . . . ? ~.~ » ~.~+~ _ ~.t = oz . .. (18) From equation (18), the number d of arrival waves can be estimated.
Then, J-d eigenvectors corresponding to eigenvalues equal to noise power are extracted to calculate a MUSIC spectrum P(t) by P(t) ° ~ H 1 H . . . ( 19 ) i=d+la (~ ) eiei a (~ ) Where a ( i ) - a j 1 , . . . , e-~wJt ~ T
w 1. . . c~ J: frequencies in use In this evaluation function, a point exhibiting a sharp peak is an estimate of the delay time.
The relative power is given by P = ( A*A ) -lA* ( S - a2I ) A ( A*A ) -1 . . . ( 2 0 ) where A = [a(i,),..., a(i~)]
I: unit matrix Finally, a delay time/relative power estimation unit 111 calculates an estimate of power from the diagonal term of the matrix P. Note that processing from the arrival signal designation unit 107 to the delay time/relative power estimation unit 111 is repetitively executed for all arrival angles of sequentially extracted signals having arrival angles estimated by the arrival angle estimation unit 106.
In the first embodiment, processing of the arrival angle estimation unit 106 and processing of the relative power estimation unit 111 respectively employ an algorithm based on the 2D-Unitary ESPRIT method and an algorithm based on the MUSIC method. Instead, the arrival angle estimation unit 106 and relative power estimation unit 111 may respectively employ the 2D-MUSIC
method and ESPRIT method.
Fig. 5 shows a radio environment analysis apparatus according to the second embodiment of the present invention.
In Fig. 5, a transmitter 1 comprises a frequency sweep signal generator 21 for generating a frequency sweep signal as a transmission signal. The transmission frequency from the transmitter 1 changes within a predetermined band by the frequency sweep signal generator 21 under the control of a signal processor 3.
Part of the transmission signal is extracted by a coupler 22 and sent to a receiver 2 via a cable 23.
The receiver 2 additionally comprises a frequency converter 17MN+1 for converting the transmission signal input via the cable 23 into an intermediate frequency or demodulated signal, similar to an array antenna reception signal, and outputting the converted signal to the signal processor. The signal processor 3 comprises a control signal generator 24 controlled by a DSP 19 in order to change the transmission signal frequency of the transmitter 1 within a predetermined band.
Signal processing operation of the DSP 19 shown in Fig. 6 will be explained.
Transmission and array antenna reception signals frequency-converted by the receiver 2 are transformed by discrete Fourier transformers (DFTs) 11211 to 112MN to obtain the amplitudes and phases. The amplitudes and phases of the transmission and array antenna reception signals are normalized in normalization units 10211 to 102MN.
The subsequent processing procedure is basically the same as in the first embodiment (see Fig. 2). Note that estimating the delay time requires array antenna reception signals having a plurality of frequencies. Hence, array outputs having a plurality of frequencies are calculated in place of the transmission frequency in order to perform processing of an array calculation unit 109.
The second embodiment requires the connection cable 23 between the transmitter 1 and receiver 2, and thus is applicable to a static or quasi-static system such as an indoor radio LAN (Local Area Network).
A radio environment analysis apparatus according to the third embodiment of the present.
invention will be described with reference to Fig. 7.
In the third embodiment, the present invention is directly applied to hardware used in the above-described references 1 and 4. That is, a receiver 2 comprises a scanner 25 and network analyzer 26, whereas a signal processor 3 comprises a personal computer 27.
As a transmission signal, a frequency sweep signal generated by a frequency sweep signal generator 21 is used. The transmission frequency sequentially changes within a predetermined band under the control of the personal computer 27 in order to estimate the delay time. Part of the transmission signal is extracted by a coupler 22, input to the network analyzer 26 via a cable 23, and used as a reference signal.
A single antenna element 16 is attached to the scanner 25 and controlled by the personal computer 27 so as to sequentially move to alignment positions of the array antenna. The third embodiment obtains an array antenna reception signal necessary for estimating the arrival elevation and azimuth angles by spatial movement of the antenna element 16.
The subsequent processing procedure is the same as in the second embodiment (see Fig. 6). Note that since the third embodiment adopts the network analyzer 26, the DFTs 11211 to 112MN and normalization units 10211 to 102MN used to normalize the reception signal can be eliminated. The delay time is estimated using a measurement value prepared by calibrating the frequency characteristics of the network analyzer 26.
As has been described above, the present invention can estimate not only the arrival azimuth angle but also the arrival elevation angle, the delay time, or the propagation parameter of relative power at the same time with high resolution and high precision.
Further, the present invention can provide a radio environment analysis apparatus applicable to even mobile communication without requiring any connection cable between communication stations and any mechanical movement of the antenna element.
The basic processing of the present invention can also be applied to hardware requiring the connection cable between communication stations and mechanical movement of the antenna element.
Fig. 1 shows a radio environment analysis apparatus according to the first embodiment of the present invention. In Fig. 1, the radio environment analysis apparatus comprises a transmitter 1 for transmitting a signal modulated by PSK (Phase Shift Keying) using a PN (Pseudo Noise) code sequence, a receiver 2 for outputting signals received by a plurality of antenna elements arrayed on the plane after conversion into intermediate frequencies or demodulated signals, and a signal processor 3 for simultaneously estimating a combination of the arrival angles (elevation and azimuth angles), delay time, and relative power of a reception signal on the basis of an output from the receiver 2.
The transmitter 1 is constituted by a PN code generator 11 for generating a PN code, a BPSK (Binary PSK) modulator 12 for phase-modulating the PN code output from the PN code generator 11, a frequency converter/amplifier 13 for converting the modulated signal output from the BPSK modulator 12 into a predetermined measurement frequency, and amplifying and outputting it, and a transmission antenna 14 for radiating the signal output from the frequency converter/amplifier 13 into the space.
The receiver 2 is constituted by a planar-array antenna 15 having M X N (M and N are integers equal to or larger than 2) antenna elements 1611 to 16MN arrayed on the plane, and a plurality of frequency converters 1711 to 17MN for outputting signals received by the antenna elements 1611 and l6MNafter conversion into intermediate frequencies or basebands.
The signal processor 3 is constituted by a plurality of A/D (Analog-to-Digital) converters (A-D) 1811 to 18MN for digitally converting the frequency-converted reception signals from the frequency converters 1711 to 17MN, and a digital signal processor (DSP) 19 for simultaneously estimating a combination of the arrival angles (elevation and azimuth angles), delay time, and relative power of a signal incoming to the planar-array antenna 15 on the basis of the plurality of digital reception signals from the A/D converters 1811 to 18MN and a criterion signal.
Signal processing operation of the DSP 19 will be explained with reference to Fig. 2.
In Fig. 2, a plurality of signals that are received by the antenna elements 1611 to 16MN of the planar-array antenna 15 and digitally converted are transformed by fast Fourier transformers 10111 to 101MN.
Then, the time-series signals are converted into signals in the frequency range. To resolve a plurality of correlative signals such as multipath signals for the - g _ delay time, a moving average must be calculated in the frequency range. In calculating the moving average, the spectrum must be flat in a target frequency range.
The group delay devices of the transmitter and receiver must be excluded from processing targets for waiting for reception.
For this purpose, in an FFT 103, a transmission signal is directly input to each channel of the receiver and subjected to fast Fourier transform.
In normalization units 10211 to 102MN, reception signals converted into the frequency range are normalized by a calibration (criterion) signal from the FFT 103 that is converted into the frequency range.
The arrival angles (elevation and azimuth angles) can be estimated from a reception signal having a frequency unique to the planar-array antenna. In a covariance matrix estimation unit 104, a moving-averaged covariance matrix is estimated using the reception signals normalized to the center frequency of the band from the normalization units 10211 to 102MN.
Fig. 3 shows moving average calculation processing by the 2D-Unitary ESPRIT method.
In Fig. 3, a complex signal received by the planar-array antenna is divided into (M-m+1) x (N-n+1), m x n minor matrices Xpq. In general, m is (M+1)/2, and n is (N+1)/2. Each divided reception signal matrix Xpq is vector-transformed as given by - g -xPq = vec (XPq) . . . ( 1 ) to calculate a correlation matrix R defined by __ H
RP4 xP9x P9 . . . ( 2 ) This operation is done for each divided minor matrix XPq to average a covariance matrix RPq by M-m+1 N-n+1 R (M - m + 1) (N - n + 1) R Pq . . . ( 3 ) P_-1 9-1 The averaged covariance matrix R is unitarily converted by QHmn - RQmn . . . ( 4 ) Where QHmn QHm ~ QHn Qmn - Qm ~ Qn Qm: mth-order unitary matrix Qn: nth-order unitary matrix ~: Kronecker's operator That is, if m(n) is an even number (m(n) - 2k:
k is an integer), Ik J Ik . ( 4A) Q2k - ' ' IIk - j IIk If m(n) is an odd number (m(n) - 2k+1), Ik 0 J Ik Q2k+1 1 ~ ~T . . .
IIk 0 j IIk where 0 is a 0 vector, Ik is the kth-order unit matrix, and IIk is defined by IIk - ' To increase resolution and precision, an ensemble average is calculated in a real matrix generation unit 105 by Ry = Re [ E [ QHmnRQmn ] ] . . . ( 5 ) where E[] is the ensemble average.
The real part of the covariance matrix is extracted.
An arrival angle estimation unit 106 estimates a combination of the arrival elevation and azimuth angles by the following procedure.
First, Ry is decomposed into eigenvalues, and the obtained eigenvalues ~.i (i = 11 to mn) are arranged in the decreasing order to be expressed by where az is noise power.
From equation (6), the number L of arrival signals can be estimated from the distribution of the eigenvalues.
Then, a partial signal space Es to be defined is calculated by ES = [elez...eL] ... (7) where ei is an eigenvector corresponding to the eigenvalue Matrices E" and E" to be defined are calculated by E _ K~,~ Es , E _ K"' Es . . . ( 8 ) N v K~z Es K~z Es where Ku 1 = IM ~ Kl, Ku z = IM ~ Kz Kvi = K3 ~ IN, K"z = K9 ~ IN
Kl = Re [ QHN_1 J2QN ] , Kz = Im [ QHN_1 JzQN ]
K3 Re [QHM_lJzQM] , Kq = Im[QHM_1JZQM]
IM: Mth-order unit matrix IN: Nth-order unit matrix Jz: selection matrix for extracting the Mth or Nth order from the second row of the matrix A matrix EXY to be defined is calculated by EXY = [KuIES~K~zEs] ... (9) Further, a matrix:
lO EXYHEXY - (KN~ES)x [Ka~Es ~ K~,zEs] . . . (10) (K,,zEs) is calculated and decomposed into eigenvalues.
A matrix E is calculated by arranging the obtained eigenvectors ei (i = 1 to 2L) by E = [elez...ezL] ... (11) The matrix E is divided into four L x L matrices by E - En E~z . . . ( 12 ) Ez~ Ez2 A matrix '~If a is calculated by . (13) ~ a = -Eiz [ Ezz ] 1 . .
Subsequently, this processing using equations (9) to (13) is repetitively executed for E" to define a matrix 'iIr" corresponding to the matrix '~',~ . 1IJ'~ + j'~I~'" is decomposed into eigenvalues. Letting the obtained eigenvalue be wu(i) + jw"(i) (i = 1 to L), an arrival elevation angle 6i and azimuth angle chi can be obtained by solving simultaneous equations:
(,JN(1) _ ~ OX COS 81 COS l~l ~ (i) _ ~ Dy cos B; sin r~l . . . ( 14 ) where O x and O y are array element intervals and ~, is the wavelength of a transmission signal.
In this way, the arrival azimuth and elevation angles of L arrival signals can be obtained.
After that, the delay time and relative power are obtained.
An arrival signal designation unit 107 designates arbitrary one of the L detected arrival signals. An adaptive weight calculation unit 108 stores the designated signal to calculate the weight of the adaptive antenna for suppressing other unnecessary signals.
Letting S" be the covariance matrix and a(c~, 8) be a steering vector indicating the direction of the stored signal when an unnecessary signal is received by the array antenna, the weight of the adaptive antenna is calculated by W = a Sua* ( ~~ a ) . . . (15) An array calculation unit 109 calculates the product of the array reception data and weight by y = W"X . . . ( 16 ) At this time, the product is calculated for a complex reception signal having a plurality of frequencies normalized in the frequency range. If J
data are used in the frequency range, array outputs y(1) to y(J) are obtained. Using the array output data y(1) to y(J), an estimation moving averaging unit 110 calculates a correlation matrix S by S = YY" ... (17) where Y = [Y(1)...y(J)1T
To resolve correlative signals for the delay time, moving averaging shown in Fig. 4 is done for the calculated correlation matrix. Fig. 4 shows a moving average calculation method using the MUSIC method for J
- 8. The correlation matrix S after moving average calculation is decomposed into eigenvalues, and the obtained eigenvalues are arranged in the decreasing order to be expressed by ~,, ? ~,z ~ . . . ? ~.~ » ~.~+~ _ ~.t = oz . .. (18) From equation (18), the number d of arrival waves can be estimated.
Then, J-d eigenvectors corresponding to eigenvalues equal to noise power are extracted to calculate a MUSIC spectrum P(t) by P(t) ° ~ H 1 H . . . ( 19 ) i=d+la (~ ) eiei a (~ ) Where a ( i ) - a j 1 , . . . , e-~wJt ~ T
w 1. . . c~ J: frequencies in use In this evaluation function, a point exhibiting a sharp peak is an estimate of the delay time.
The relative power is given by P = ( A*A ) -lA* ( S - a2I ) A ( A*A ) -1 . . . ( 2 0 ) where A = [a(i,),..., a(i~)]
I: unit matrix Finally, a delay time/relative power estimation unit 111 calculates an estimate of power from the diagonal term of the matrix P. Note that processing from the arrival signal designation unit 107 to the delay time/relative power estimation unit 111 is repetitively executed for all arrival angles of sequentially extracted signals having arrival angles estimated by the arrival angle estimation unit 106.
In the first embodiment, processing of the arrival angle estimation unit 106 and processing of the relative power estimation unit 111 respectively employ an algorithm based on the 2D-Unitary ESPRIT method and an algorithm based on the MUSIC method. Instead, the arrival angle estimation unit 106 and relative power estimation unit 111 may respectively employ the 2D-MUSIC
method and ESPRIT method.
Fig. 5 shows a radio environment analysis apparatus according to the second embodiment of the present invention.
In Fig. 5, a transmitter 1 comprises a frequency sweep signal generator 21 for generating a frequency sweep signal as a transmission signal. The transmission frequency from the transmitter 1 changes within a predetermined band by the frequency sweep signal generator 21 under the control of a signal processor 3.
Part of the transmission signal is extracted by a coupler 22 and sent to a receiver 2 via a cable 23.
The receiver 2 additionally comprises a frequency converter 17MN+1 for converting the transmission signal input via the cable 23 into an intermediate frequency or demodulated signal, similar to an array antenna reception signal, and outputting the converted signal to the signal processor. The signal processor 3 comprises a control signal generator 24 controlled by a DSP 19 in order to change the transmission signal frequency of the transmitter 1 within a predetermined band.
Signal processing operation of the DSP 19 shown in Fig. 6 will be explained.
Transmission and array antenna reception signals frequency-converted by the receiver 2 are transformed by discrete Fourier transformers (DFTs) 11211 to 112MN to obtain the amplitudes and phases. The amplitudes and phases of the transmission and array antenna reception signals are normalized in normalization units 10211 to 102MN.
The subsequent processing procedure is basically the same as in the first embodiment (see Fig. 2). Note that estimating the delay time requires array antenna reception signals having a plurality of frequencies. Hence, array outputs having a plurality of frequencies are calculated in place of the transmission frequency in order to perform processing of an array calculation unit 109.
The second embodiment requires the connection cable 23 between the transmitter 1 and receiver 2, and thus is applicable to a static or quasi-static system such as an indoor radio LAN (Local Area Network).
A radio environment analysis apparatus according to the third embodiment of the present.
invention will be described with reference to Fig. 7.
In the third embodiment, the present invention is directly applied to hardware used in the above-described references 1 and 4. That is, a receiver 2 comprises a scanner 25 and network analyzer 26, whereas a signal processor 3 comprises a personal computer 27.
As a transmission signal, a frequency sweep signal generated by a frequency sweep signal generator 21 is used. The transmission frequency sequentially changes within a predetermined band under the control of the personal computer 27 in order to estimate the delay time. Part of the transmission signal is extracted by a coupler 22, input to the network analyzer 26 via a cable 23, and used as a reference signal.
A single antenna element 16 is attached to the scanner 25 and controlled by the personal computer 27 so as to sequentially move to alignment positions of the array antenna. The third embodiment obtains an array antenna reception signal necessary for estimating the arrival elevation and azimuth angles by spatial movement of the antenna element 16.
The subsequent processing procedure is the same as in the second embodiment (see Fig. 6). Note that since the third embodiment adopts the network analyzer 26, the DFTs 11211 to 112MN and normalization units 10211 to 102MN used to normalize the reception signal can be eliminated. The delay time is estimated using a measurement value prepared by calibrating the frequency characteristics of the network analyzer 26.
As has been described above, the present invention can estimate not only the arrival azimuth angle but also the arrival elevation angle, the delay time, or the propagation parameter of relative power at the same time with high resolution and high precision.
Further, the present invention can provide a radio environment analysis apparatus applicable to even mobile communication without requiring any connection cable between communication stations and any mechanical movement of the antenna element.
The basic processing of the present invention can also be applied to hardware requiring the connection cable between communication stations and mechanical movement of the antenna element.
Claims (7)
1. A radio environment analysis apparatus comprising:
a receiver for receiving a signal modulated by phase shift keying (PSK) using a pseudo noise (PN) code sequence by an array antenna at a plurality of fixed points on a plane, converting the plurality of received signals into intermediate frequencies or demodulated signals, and outputting the resultant signals; and a signal processor for processing a signal output from said receiver to analyze an arrival angle, delay time, and relative power of the signal incoming to said receiver, said signal processor having:
a plurality of normalization means for individually normalizing respective signals from said receiver and outputting the normalized signals as a normalized signal group;
first estimation means for calculating a covariance matrix from the signal group normalized by said normalization means, and estimating an arrival angle of each incoming signal using one of a 2D-MUSIC
method and 2D-Unitary ESPRIT method;
arithmetic means for calculating a weight of said array antenna for suppressing incoming signals except for an incoming signal having the arrival angle estimated by said first estimation means, and calculating a product of the calculated weight of said array antenna and the normalized signal group; and second estimation means for calculating a covariance matrix from a product signal group output from said arithmetic means, and estimating a delay time and relative power using one of the MUSIC method and ESPRIT method.
a receiver for receiving a signal modulated by phase shift keying (PSK) using a pseudo noise (PN) code sequence by an array antenna at a plurality of fixed points on a plane, converting the plurality of received signals into intermediate frequencies or demodulated signals, and outputting the resultant signals; and a signal processor for processing a signal output from said receiver to analyze an arrival angle, delay time, and relative power of the signal incoming to said receiver, said signal processor having:
a plurality of normalization means for individually normalizing respective signals from said receiver and outputting the normalized signals as a normalized signal group;
first estimation means for calculating a covariance matrix from the signal group normalized by said normalization means, and estimating an arrival angle of each incoming signal using one of a 2D-MUSIC
method and 2D-Unitary ESPRIT method;
arithmetic means for calculating a weight of said array antenna for suppressing incoming signals except for an incoming signal having the arrival angle estimated by said first estimation means, and calculating a product of the calculated weight of said array antenna and the normalized signal group; and second estimation means for calculating a covariance matrix from a product signal group output from said arithmetic means, and estimating a delay time and relative power using one of the MUSIC method and ESPRIT method.
2. An apparatus according to claim 1, wherein said array antenna comprises a plurality of antenna elements arrayed on the plane, and said normalization means normalizes respective outputs from said plurality of antenna elements.
3. An apparatus according to claim 1, wherein said array antenna comprises one antenna element moving through fixed points on the plane, and said normalization means normalizes respective outputs from said antenna element at the fixed points.
4. An apparatus according to claim 1, wherein said normalization means performs Fourier transform for a signal from said receiver and normalizes the signal with the same spectrum as a transmission signal.
5. An apparatus according to claim 1, wherein said apparatus further comprises a transmitter for transmitting a signal PSK-modulated using a PN code sequence, and said normalization means normalizes the signal from said receiver by a frequency sweep signal supplied from said transmitter.
6. An apparatus according to claim 1, wherein said first estimation means estimates a combination of arrival elevation and azimuth angles as an arrival angle of an incoming signal, and said arithmetic means stores information of the incoming signal having the arrival elevation and azimuth angles estimated by said first estimation means, and calculates a weight of said array antenna for suppressing other incoming signals in accordance with the storage result.
7. An apparatus according to claim 1, wherein said signal processor stores each incoming signal and calculates a weight for suppressing other incoming signals to calculate a normalized signal group and product for all incoming signals having estimated arrival angles, and calculates a delay time and relative power from a covariance matrix of the obtained product signal group to estimate a combination of the arrival angles, delay times, and relative powers of all signals incoming to said receiver
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP152629/98 | 1998-06-02 | ||
| JP10152629A JPH11344517A (en) | 1998-06-02 | 1998-06-02 | Environmental analyzer for electric wave |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2273530A1 CA2273530A1 (en) | 1999-12-02 |
| CA2273530C true CA2273530C (en) | 2002-10-15 |
Family
ID=15544573
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002273530A Expired - Fee Related CA2273530C (en) | 1998-06-02 | 1999-06-01 | Radio environment analysis apparatus |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US6084928A (en) |
| EP (1) | EP0963003B1 (en) |
| JP (1) | JPH11344517A (en) |
| KR (1) | KR100340832B1 (en) |
| AU (1) | AU748803B2 (en) |
| CA (1) | CA2273530C (en) |
| DE (1) | DE69932772T2 (en) |
Families Citing this family (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19754031A1 (en) * | 1997-12-05 | 1999-06-17 | Siemens Ag | Method and measuring arrangement for measuring the properties of radio channels |
| JP3233088B2 (en) * | 1998-01-22 | 2001-11-26 | 松下電器産業株式会社 | Directivity control antenna device |
| US6311043B1 (en) * | 1998-10-27 | 2001-10-30 | Siemens Aktiengesellschaft | Method and measurement configuration for measuring the characteristics of radio channels |
| KR100487243B1 (en) * | 1998-12-17 | 2005-08-31 | 삼성전자주식회사 | Device and method for estimating the position of terminal in mobile communication system |
| US6351238B1 (en) * | 1999-02-23 | 2002-02-26 | Matsushita Electric Industrial Co., Ltd. | Direction of arrival estimation apparatus and variable directional signal receiving and transmitting apparatus using the same |
| DE19925580B4 (en) * | 1999-06-04 | 2010-09-09 | Eads Deutschland Gmbh | Level and bearing angle-based detection and segmentation of programs in the HF range |
| JP4169884B2 (en) * | 1999-09-24 | 2008-10-22 | 富士通株式会社 | Communication device using adaptive antenna |
| JP3315955B2 (en) * | 1999-09-29 | 2002-08-19 | 松下電器産業株式会社 | Base station apparatus and wireless communication method |
| JP3331595B2 (en) * | 2000-06-05 | 2002-10-07 | 日東紡績株式会社 | Indoor environment design system, indoor environment evaluation system, indoor environment design method, and indoor environment evaluation method |
| JP4306098B2 (en) * | 2000-06-30 | 2009-07-29 | 株式会社デンソー | Communication device |
| EP1387181A1 (en) * | 2000-12-12 | 2004-02-04 | Matsushita Electric Industrial Co., Ltd. | Radio-wave arrival-direction estimating apparatus and directional variable transceiver |
| GB2376567B (en) * | 2001-06-12 | 2005-07-20 | Mobisphere Ltd | Improvements in or relating to smart antenna arrays |
| US6600446B2 (en) * | 2001-06-29 | 2003-07-29 | Lockheed Martin Corporation | Cascadable architecture for digital beamformer |
| US7181167B2 (en) * | 2001-11-21 | 2007-02-20 | Texas Instruments Incorporated | High data rate closed loop MIMO scheme combining transmit diversity and data multiplexing |
| US7684525B2 (en) * | 2002-01-09 | 2010-03-23 | Nokia Corporation | Method and receiver for reception of a composite signal |
| US6980614B2 (en) * | 2002-01-14 | 2005-12-27 | Raytheon Company | System and method for subband beamforming using adaptive weight normalization |
| KR100526299B1 (en) * | 2002-10-04 | 2005-11-03 | 송기홍 | An apparatus for converting the propagation signal according to the propagation and its method |
| JP4808984B2 (en) | 2005-03-31 | 2011-11-02 | 株式会社デンソーアイティーラボラトリ | Arrival wave direction estimation device |
| DE102005023174B4 (en) * | 2005-05-19 | 2016-03-03 | Robert Bosch Gmbh | A method for interface adaptation of a hardware baseband receiver in satellite communication systems, interface adapter for hardware baseband receivers and a corresponding computer program and a corresponding computer-readable storage medium |
| JP4757629B2 (en) | 2005-12-28 | 2011-08-24 | 株式会社デンソーアイティーラボラトリ | Arrival direction estimation device |
| JP4576515B2 (en) * | 2007-03-06 | 2010-11-10 | 学校法人慶應義塾 | Event detection device |
| JPWO2009028080A1 (en) * | 2007-08-30 | 2010-11-25 | 富士通株式会社 | Advertisement distribution system, call control device, and advertisement distribution method |
| KR101330082B1 (en) * | 2007-10-15 | 2013-11-18 | 삼성전자주식회사 | Apparatus and method for coherent sources direction of arrival using polarization in multipath system |
| WO2009081981A1 (en) * | 2007-12-25 | 2009-07-02 | Honda Elesys Co., Ltd. | Electronic scanning radar apparatus, received wave direction estimating method, and received wave direction estimating program |
| JP5554018B2 (en) * | 2009-06-12 | 2014-07-23 | 三菱電機株式会社 | Angle measuring device, monopulse angle measuring device, monopulse radar, multi-static radar |
| JP5645211B2 (en) * | 2009-08-28 | 2014-12-24 | 一般財団法人電力中央研究所 | Method and apparatus for detecting direction of arrival of pulsed radio noise |
| JP5494567B2 (en) | 2011-05-17 | 2014-05-14 | 株式会社デンソー | Radar apparatus, inspection system, and inspection method |
| GB2502108B (en) * | 2012-05-16 | 2014-10-15 | Canon Kk | Reception quality assessment |
| US20220209875A1 (en) * | 2019-04-11 | 2022-06-30 | Nec Corporation | Spatial data creating apparatus, spatial data creating method, and program |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4965732A (en) * | 1985-11-06 | 1990-10-23 | The Board Of Trustees Of The Leland Stanford Junior University | Methods and arrangements for signal reception and parameter estimation |
| JPH0476466A (en) * | 1990-07-18 | 1992-03-11 | Sanki Denshi Kogyo Kk | Electromagnetic wave measuring device |
| JP2901170B2 (en) * | 1993-05-27 | 1999-06-07 | ケイディディ株式会社 | Satellite / land mobile communication system integration method |
| JPH0865743A (en) * | 1994-08-23 | 1996-03-08 | N T T Ido Tsushinmo Kk | Mobile communication system |
| JP2937068B2 (en) * | 1995-04-07 | 1999-08-23 | 三菱電機株式会社 | Direction detection device |
| US5930243A (en) * | 1996-10-11 | 1999-07-27 | Arraycomm, Inc. | Method and apparatus for estimating parameters of a communication system using antenna arrays and spatial processing |
-
1998
- 1998-06-02 JP JP10152629A patent/JPH11344517A/en active Pending
-
1999
- 1999-06-01 AU AU32365/99A patent/AU748803B2/en not_active Ceased
- 1999-06-01 CA CA002273530A patent/CA2273530C/en not_active Expired - Fee Related
- 1999-06-01 EP EP99109358A patent/EP0963003B1/en not_active Expired - Lifetime
- 1999-06-01 US US09/323,782 patent/US6084928A/en not_active Expired - Lifetime
- 1999-06-01 DE DE69932772T patent/DE69932772T2/en not_active Expired - Lifetime
- 1999-06-02 KR KR1019990020159A patent/KR100340832B1/en not_active IP Right Cessation
Also Published As
| Publication number | Publication date |
|---|---|
| EP0963003A3 (en) | 2001-02-28 |
| EP0963003B1 (en) | 2006-08-16 |
| JPH11344517A (en) | 1999-12-14 |
| KR20000005826A (en) | 2000-01-25 |
| CA2273530A1 (en) | 1999-12-02 |
| DE69932772D1 (en) | 2006-09-28 |
| DE69932772T2 (en) | 2007-08-16 |
| AU748803B2 (en) | 2002-06-13 |
| EP0963003A2 (en) | 1999-12-08 |
| KR100340832B1 (en) | 2002-06-15 |
| AU3236599A (en) | 1999-12-09 |
| US6084928A (en) | 2000-07-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2273530C (en) | Radio environment analysis apparatus | |
| US7057555B2 (en) | Wireless LAN with distributed access points for space management | |
| US6489923B1 (en) | Position location method and apparatus for a mobile telecommunications system | |
| JP6278961B2 (en) | Radar equipment | |
| US6642888B2 (en) | Radio-wave arrival-direction estimating apparatus and directional variable transceiver | |
| US6529745B1 (en) | Radio wave arrival direction estimating antenna apparatus | |
| US20190369223A1 (en) | Radar apparatus | |
| Kumari et al. | A low-resolution ADC proof-of-concept development for a fully-digital millimeter-wave joint communication-radar | |
| Lin et al. | 3D wideband mmWave localization for 5G massive MIMO systems | |
| US11754658B2 (en) | Radio station for client localization in multipath indoor environment | |
| WO2019068305A1 (en) | Receiving device and methods thereof | |
| Saleem et al. | Investigation of massive MIMO channel spatial characteristics for indoor subway tunnel environment | |
| Shahid et al. | Comparative analysis between direction of arrival algorithms | |
| CA2293940A1 (en) | Method for estimating arrival direction of desired wave | |
| CA2244369A1 (en) | Multiplex radio communication apparatus | |
| EP3632000B1 (en) | Determination of electrical phase relation in a communications network | |
| JP2000138520A (en) | Antenna system | |
| JP2004061468A (en) | Method for estimating arrival direction of multiplex wave using spatial characteristic and receiving beam formation device using the same | |
| Wang et al. | MUSIC Algorithm for IRS-Assisted AOA Estimation | |
| Ullah et al. | Power scalable angle of arrival estimation using pilot design with orthogonal subsequences | |
| El Arja et al. | Joint TOA/DOA measurements for UWB indoor propagation channel using MUSIC algorithm | |
| Neesen et al. | A novel doppler-based direction-of-arrival estimation scheme | |
| Shaikh et al. | Complexity reduced direction finding massive MIMO system using EM lens | |
| Arai et al. | Hardware implementation of smart antenna systems for high speed wireless communication | |
| Jami et al. | Improved method for estimating angle of arrival in multipath conditions using the'MUSIC'algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| MKLA | Lapsed |
Effective date: 20130603 |
|
| MKLA | Lapsed |
Effective date: 20130603 |