WO2018025421A1 - 物体検知装置および物体検知方法 - Google Patents
物体検知装置および物体検知方法 Download PDFInfo
- Publication number
- WO2018025421A1 WO2018025421A1 PCT/JP2016/073205 JP2016073205W WO2018025421A1 WO 2018025421 A1 WO2018025421 A1 WO 2018025421A1 JP 2016073205 W JP2016073205 W JP 2016073205W WO 2018025421 A1 WO2018025421 A1 WO 2018025421A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- object detection
- signal
- antenna
- reception
- transmission
- Prior art date
Links
Images
Classifications
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/04—Systems determining presence of a target
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S13/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
- G01S13/343—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal using sawtooth modulation
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
Definitions
- the present invention relates to an object detection apparatus and an object detection method for detecting an object from radio waves reflected from the object or radiated from the object.
- radio waves microwaves, millimeter waves, terahertz waves, etc.
- An imaging apparatus that uses this radio wave transmission capability to image and inspect articles hidden under clothes or articles in bags has been put into practical use.
- a remote sensing technique for imaging the ground surface through a cloud from a satellite or an aircraft has been put into practical use.
- FIG. 22 is a diagram showing an object detection apparatus employing a conventional array antenna system.
- FIG. 23 is a diagram illustrating a configuration of the receiver illustrated in FIG.
- the object detection apparatus includes a transmitter 211 and a receiver 201.
- the transmitter 211 includes a transmission antenna 212.
- the receiver 201 includes receiving antennas 201 1 , 202 2 ,..., 202 N (N is the number of receiving antennas).
- the transmitter 211 irradiates an RF signal (radio wave) 213 from the transmitting antenna 212 toward the detection objects 204 1 , 204 2 ,..., 204 K (K is the number of objects).
- RF signal (radio wave) 213, the detection object 204 1, 204 2, ..., are reflected in the 204 K, the reflected wave 203 1, 203 2, ..., 203 K is generated, respectively.
- Reflected wave 203 1 generated, 203 2, ..., 203 L may receive antennas 201 1, 202 2, ..., is received at 202 N.
- the receiver 201 calculates reflected wave 203 1 receiving, 203 2, ..., based on 203 K, the detection object 204 1, 204 2, ..., a radio field intensity of the radio wave is reflected by 204 K To do. After that, the receiver 201 images the calculated radio wave intensity distribution. Thereby, the images of the detection objects 204 1 , 204 2 ,..., 204 K are obtained.
- the complex amplitudes of the incoming waves 208 1 , 208 2 ,..., 208 K are [s ( ⁇ 1 ), s ( ⁇ 2 ),..., S ( ⁇ K )].
- Receiver 201 down converter provided with a (non-shown in FIG.
- the down converter each receiving antenna 202 1, 202 2, ..., complex amplitude (baseband signal) of the RF signal received by 202 N [r 1 , r 2 ,..., r N ] are extracted.
- the complex amplitudes [r 1 , r 2 ,..., R N ] of the signals received by the receiving antennas 202 1 , 202 2 ,..., 202 N are output to the signal processing unit 205.
- Receiving antennas 202 1 , 202 2 ,..., 202 N receive signal complex amplitudes [r 1 , r 2 ,..., R N ] and incoming wave complex amplitudes [s ( ⁇ 1 ), s
- the relationship with ( ⁇ 2 ),..., s ( ⁇ K )] is given by the following equation (1).
- n (t) is a vector whose elements are noise components.
- the subscript T represents the transpose of a vector or matrix.
- d is the distance between the antennas
- lambda is incoming wave (RF signal) 208 1, 208 2, ..., a wavelength of 208 K.
- the complex amplitude r of the received signal is an amount obtained by measurement.
- the direction matrix A is an amount that can be defined (designated) in signal processing.
- the complex amplitude s of the incoming wave is an unknown, and the purpose of the incoming wave direction estimation is to determine the direction of the incoming wave s from the received signal r obtained by measurement.
- a correlation matrix R E [r ⁇ r H ] is calculated from the received signal r obtained by measurement.
- E [] represents that time-average processing is performed on the elements in parentheses, and the subscript H represents complex conjugate transpose.
- any evaluation function represented by the following equations (2) to (4) is calculated.
- E N [e K + 1 ,..., E N ] in the MUSIC method is N ⁇ (K + 1) vectors whose eigenvalue is the power of noise n (t) among the eigenvectors of the correlation matrix R It is a matrix composed of
- the evaluation function shown in Equation (2) to (4) the angle theta 1 of the incoming waves, theta 2, ⁇ ⁇ ⁇ , with a peak at theta K. Therefore, the angle of the incoming wave can be obtained by calculating the evaluation function and looking at the peak.
- the position and shape of the object can be displayed as an image from the angular distribution of the incoming waves obtained by the evaluation functions of Expressions (2) to (4).
- FIG. 24 is a diagram illustrating an example in which the beamformer method is applied in the receiver illustrated in FIG.
- phase shifters 206 1 , 206 2 ,..., 206 N and the combiner 207 of the conventional antenna array shown in FIG. 24 correspond to the signal processing unit 205 in the conventional antenna array shown in FIG. .
- Phase shifter 206 1, 206 2, ..., 206 N respectively, receive antennas 202 1, 202 2, ..., the complex amplitude 208 1 of the incoming wave received by the 202 N, 208 2, ... , to 208 N, a phase rotation ⁇ 1, ⁇ 2, ⁇ , is added [Phi N.
- Phase rotation ⁇ 1, ⁇ 2, ⁇ , ⁇ N is incoming wave 208 1 applied, 208 2, ⁇ ⁇ ⁇ , 208 N are added by the adder 207.
- the directivity component AF ( ⁇ ) obtained by removing the directivity g ( ⁇ ) of the receiving antenna 202 from the directivity E ( ⁇ ) of the array antenna is called an array factor.
- the array factor AF ( ⁇ ) represents the effect of directivity due to the formation of the array antenna.
- Reception antenna 202 n (n 1,2, ⁇ , N) signal received in a g ( ⁇ ) a n exp ( j ⁇ n).
- the signal obtained by adding in (5) is obtained as the directivity E ( ⁇ ) of equation (5).
- d is the receiving antenna 202
- lambda is arriving wave 208 1, 208 2, ..., a wavelength of 208 N.
- array factor AF theta
- Patent Documents 1 to 3 Other examples of the object detection apparatus using the array antenna method are also disclosed in Patent Documents 1 to 3.
- the object detection devices disclosed in Patent Documents 1 and 2 are receptions formed by N reception antennas by phase shifters connected to the N reception antennas incorporated in the receiver. Controls the directivity of the array antenna.
- the object detection devices disclosed in Patent Documents 1 and 2 change the directivity of the N reception array antennas formed in a beam shape, and each of the K detection objects has the reception array antenna. Direct the directional beam. Thereby, the radio wave intensity reflected by each detection object is calculated.
- the object detection device disclosed in Patent Document 3 controls the directivity of the N reception array antennas by using the frequency dependence of the N reception array antennas. Similarly to the examples of Patent Documents 1 and 2, the object detection device disclosed in Patent Document 3 also directs the directional beams of N reception array antennas to each of the K detection objects. The radio field intensity reflected by each detection object is calculated.
- FIG. 25 is a diagram showing a schematic configuration of a receiving array antenna when a conventional array antenna system is adopted.
- FIG. 26 is a diagram illustrating an object detection device that employs the Mills-Cross method.
- the object detection apparatus includes a one-dimensional array antenna 201 arranged in the vertical direction and a one-dimensional array antenna 201 arranged in the horizontal direction.
- the multiplier 221 calculates a product of signals for each set of the reception antenna in the vertical direction and the reception antenna in the horizontal direction. Therefore, a two-dimensional image can be displayed by using the calculated product.
- FIG. 27 is a diagram showing an object detection apparatus that employs a conventional synthetic aperture radar system.
- the object detection device includes a transmitter 311 and a receiver 301.
- the transmitter 311 is provided with a transmission antenna 312.
- the receiver 301 includes receiving antennas 302 1 to 302 N (N is the number of receiving antennas).
- the transmitter 311 emits an RF signal (radio wave) 313 from the transmission antenna 312 toward the detection objects 304 1 , 304 2 ,..., 304 K (K is the number of detection objects).
- RF signal (radio wave) 313, the detection object 304 1, 304 2, ..., are reflected in the 304 K, the reflected wave 303 1, 303 2, ..., 303 L are generated, respectively.
- receiver 301 from the initial position, 301 2, ..., while moving to the position of 301 N, receives the reflected wave 303 1, 303 2, ..., and 303 K at each position .
- 302 1 , 302 2 ,..., 302 N indicate reception antennas at respective positions.
- one receiving antenna functions as receiving antennas 302 1 , 302 2 ,..., 302 N. That is, in FIG. 27, one receiving antenna is a receiving array antenna (virtual array) with N antennas, like the receiving antennas 202 1 , 202 2 ,..., 202 N in the array antenna system shown in FIG. Array antenna).
- a receiving array antenna virtual array with N antennas, like the receiving antennas 202 1 , 202 2 ,..., 202 N in the array antenna system shown in FIG. Array antenna).
- the receiver 301 is based on the received reflected waves 303 1 , 303 2 ,..., 303 K , similarly to the array antenna system shown in FIG.
- the radio wave intensity reflected from the detection objects 304 1 , 304 2 ,..., 304 K is calculated. Thereafter, the receiver 301 images the calculated distribution of radio field intensity. As a result, images of the detection objects 304 1 , 304 2 ,..., 304 K are obtained.
- Patent Documents 4 to 6 disclose examples of an object detection device using a synthetic aperture radar system.
- the receiving antenna 201 1, 202 2, ..., 202 intervals of each antenna of N reflected wave 203 1 received at the receiver 201, 203 2, ..., 203 K Must be less than half of the wavelength ⁇ .
- the wavelength ⁇ is about several millimeters, so that the interval between the antennas is several millimeters or less. If this condition is not satisfied, there arises a problem that a virtual image is generated at a position where the objects 204 1 , 204 2 ,..., 204 K do not exist in the generated image.
- the resolution of the image is determined by the directional beam width ⁇ of the receiving array antenna (201 1 , 202 2 ,..., 202 N ).
- the width ⁇ of the directional beam of the receiving array antenna (201 1 , 202 2 ,..., 202 N ) is given by ⁇ to ⁇ / D.
- D is the opening size of the receiving array antenna (201 1 , 202 2 ,... 202 N ), and corresponds to the distance between the receiving antennas 202 1 and 202 N existing at both ends. That is, in order to obtain a practical resolution in imaging an article hidden under clothes or an article in a bag, the openings of the receiving array antennas (201 1 , 202 2 ,..., 202 N ).
- the size D needs to be set to about several tens of centimeters to several meters.
- the distance between the antennas of the N receiving antennas should be less than half of the wavelength ⁇ (several mm or less) and the distance between the receiving antennas existing at both ends should be at least about several tens of centimeters. From this point, the number N of antennas required per row is about several hundred.
- N reception antennas 202 are installed in each of the vertical direction and the horizontal direction.
- the total number of reception antennas required is N 2 . Therefore, in order to employ the array antenna system, the number of reception antennas necessary for the whole and the number of receivers associated therewith are approximately tens of thousands.
- the cost is extremely high in the array antenna system.
- the antenna is installed in four regions with one side of several tens of centimeters to several meters, the size and weight of the device are very large.
- the number of receiving antennas and receivers can be reduced as compared with the case where the array antenna method is adopted.
- the number of necessary reception antennas and receivers is 2N, and about several hundred reception antennas are required. Therefore, even in this case, it is difficult to solve the problems of cost, device size and weight.
- An example of an object of the present invention is an object detection device and an object detection that can solve the above-described problems and can suppress an increase in device cost, size, and weight while improving accuracy in detecting an object using radio waves. It is to provide a method.
- an object detection device for detecting an object by radio waves, A transmitter that radiates radio waves whose frequency changes continuously over time as a transmission signal; A receiving unit that acquires the transmission signal, receives the radio wave reflected by the object as a reception signal, and further mixes the acquired transmission signal with the received reception signal to generate a baseband signal.
- a data processing unit for detecting the object It is characterized by having.
- an object detection method is a method for detecting an object by radio waves, (A) radiating, as a transmission signal, a radio wave whose frequency continuously changes over time by a transmitter; (B) The transmission signal is acquired by a receiver, the radio wave reflected by the object is received as a reception signal, and the acquired transmission signal is added to the received reception signal to obtain a baseband. Generating a signal; and (C) The data processor estimates the direction of arrival of the radio wave from the measured value of the baseband signal at each sampling time, and specifies the intensity distribution of the radio wave based on the estimated direction of arrival of the radio wave, Detecting the object based on the identified intensity distribution; and It is characterized by having.
- FIG. 1 is a configuration diagram schematically showing a configuration of an object detection device according to Embodiment 1 of the present invention.
- FIG. 2 is a configuration diagram for explaining the operation principle of the object detection device according to Embodiment 1 of the present invention.
- FIG. 3 is a diagram showing a change in frequency of the radio wave transmitted in the first embodiment of the present invention.
- FIG. 4 is a diagram showing a correspondence relationship between parameters used in the conventional array antenna system and parameters used in the time virtual array system in the embodiment of the present invention.
- FIG. 5 is a diagram illustrating an operation principle of the object detection device according to Embodiment 1 of the present invention.
- FIG. 6 is a diagram showing an operation principle when the beamformer method is applied to the object detection apparatus shown in FIG. FIG.
- FIG. 7 is a characteristic diagram showing an example of the antenna gain directivity of the object detection device according to Embodiment 1 of the present invention.
- FIG. 8 is a diagram for explaining generation of a virtual image in a virtual array in one embodiment of the present invention.
- FIG. 9 is a block diagram showing an example of a specific configuration of the object detection apparatus according to Embodiment 1 of the present invention.
- FIG. 10 is a flowchart showing the operation of the object detection apparatus according to Embodiment 1 of the present invention.
- FIG. 11 is a diagram showing the configuration and operation principle of the object detection device according to the second embodiment of the present invention.
- FIG. 12 is a diagram for explaining the concept of a subarray used in the object detection device according to Embodiment 2 of the present invention.
- FIG. 13 is a flowchart showing the operation of the object detection apparatus according to the second embodiment of the present invention.
- FIG. 14 is a diagram showing the configuration and operation principle of the object detection device according to Embodiment 3 of the present invention.
- FIG. 15 is an explanatory diagram for explaining a correlation matrix calculation method for a two-dimensional frequency virtual array according to Embodiment 3 of the present invention.
- FIG. 16 is an explanatory diagram illustrating a method for calculating a correlation matrix of a two-dimensional frequency virtual array according to Embodiment 3 of the present invention.
- FIG. 17 is a flowchart showing the operation of the object detection apparatus according to the third embodiment of the present invention.
- FIG. 18 is a diagram illustrating an example of an image output from the object detection device according to Embodiment 3 of the present invention.
- FIG. 19 is a diagram illustrating a schematic configuration of the object detection device according to the fourth embodiment of the present invention.
- FIG. 20 is a block diagram specifically showing the configuration of the object detection apparatus according to Embodiment 4 of the present invention.
- FIG. 21 is a diagram showing an example of frequency control performed in the fourth embodiment of the present invention.
- FIG. 22 is a diagram showing an object detection apparatus employing a conventional array antenna system.
- FIG. 23 is a diagram illustrating a configuration of the receiver illustrated in FIG.
- FIG. 24 is a diagram illustrating an example in which the beamformer method is applied in the receiver illustrated in FIG. FIG.
- FIG. 25 is a diagram showing a schematic configuration of a receiving array antenna when a conventional array antenna system is adopted.
- FIG. 26 is a diagram illustrating an object detection device that employs the Mills-Cross method.
- FIG. 27 is a diagram showing an object detection apparatus that employs a conventional synthetic aperture radar system.
- Embodiment 1 (Embodiment 1)
- an object detection apparatus and an object detection method according to Embodiment 1 of the present invention will be described with reference to FIGS.
- FIG. 1 is a configuration diagram schematically showing a configuration of an object detection device according to Embodiment 1 of the present invention.
- the object detection apparatus 1000 includes a transmission unit 1091, a reception unit 1092, and a data processing unit 1093.
- the transmission unit 1091 radiates a radio wave whose frequency continuously changes with time as a transmission signal.
- the receiving unit 1092 acquires a transmission signal, and receives a radio wave from an object (hereinafter referred to as “target object”) 1001 to be detected as a reception signal. Further, the reception unit 1092 multiplies (mixes) the acquired transmission signal by the received reception signal to generate a baseband signal.
- the transmission unit 1091 includes a transmission antenna 1003, and the reception unit 1092 includes a reception antenna 1004.
- the number of receiving units 1092 and receiving antennas 1004 may be plural. However, in the first embodiment, the number of receiving units 1092 and receiving antennas 1004 is extremely small compared to the conventional case.
- the data processing unit 1093 estimates the arrival direction of radio waves from the measured value of the baseband signal for each sampling time. Then, the data processing unit 1093 identifies the radio wave intensity distribution based on the estimated radio wave arrival direction, and detects the object 1001 based on the identified intensity distribution.
- FIG. 2 is a diagram for explaining the operation principle of the object detection device according to Embodiment 1 of the present invention.
- FIG. 3 is a diagram showing a change in frequency of the radio wave transmitted in the first embodiment of the present invention.
- a transmitting antenna 1003 and a receiving antenna 1004 are arranged on the x-axis, and K objects 1001 1 ,... 1001 K are (x 1 , z 0 ),. , (x K , z 0 ), respectively.
- an RF signal 1010 whose carrier frequency f changes linearly is transmitted from the transmission antenna 1003.
- Carrier frequency f is between 1 chirp period (T chirp), it shall be changed from the minimum frequency f min to the maximum frequency f max.
- object 1001 1 from the transmitting antenna 1003, RF signal 1010 1 towards ⁇ ⁇ ⁇ 1001 K, ⁇ ⁇ ⁇ 1010 K is assumed to be irradiated, respectively. Furthermore, the object 1001 1, the reflected wave 1007 1 from ⁇ ⁇ ⁇ 1001 K, ⁇ ⁇ ⁇ 1007 K are intended to be received by the receiving antenna 1004.
- the combined wave of the reflected waves 1007 1 ,... 1007 K received by the receiving antenna 1004 is multiplied (mixed) with the transmission signal acquired from the transmitting unit 1091 in the receiving unit 1092 shown in FIG. As a result, a baseband signal is generated.
- the baseband signal I (t) is given by the following formula (6).
- t ′ is a time within one chirp period, and corresponds to t 0 to t M in FIG.
- the h as chirp number, 'as denoted as t-h ⁇ T chirp, t for each pass a chirp period (T chirp)' t need to take heart from t 0.
- L (x k ) is a propagation distance from the transmitting antenna to the receiving antenna via the object k.
- c is the speed of light.
- the above IF signal I (t) is an in-phase component (In-phase signal).
- a quadrature component (Quadrature signal) Q (t) is generated.
- In-phase component I (t) and quadrature component Q (t) may be generated using a quadrature modulator.
- the orthogonal component Q (t) is given by the following equation (7).
- a complex baseband signal r (t) represented by the following equation (8) is generated from the in-phase component I (t) and the quadrature component Q (t).
- the complex baseband signal r (t) is an amount that can be calculated from measured data.
- t 1 , t 2 ,..., T N are sampling times within one chirp period.
- N is the number of sampling points per chirp period.
- a vector n (t) having a noise component (random number) as an element is added to the received signal r.
- FIG. 4 is a diagram showing a correspondence relationship between parameters used in the conventional array antenna system and parameters used in the time virtual array system in the embodiment of the present invention.
- a correlation matrix R E [r ⁇ r H ] is calculated from the received signal (complex baseband signal) r defined by the equation (9) obtained by measurement, and Then, from the calculated correlation matrix R, one of the evaluation functions shown in the following equations (10) to (12) is calculated.
- Equation 9 the direction vector a (x) is the one defined by equation (9).
- E N [e K + 1 ,..., E N ] in the MUSIC method is N ⁇ (K + 1) in which the eigenvalue is the power of the noise n (t) among the eigenvectors of the correlation matrix R. ) A matrix composed of vectors.
- the conventional antenna array estimates the direction of arrival of radio waves using data received by N antennas arranged with an antenna interval d. It can be interpreted that the arrival direction of the radio wave is estimated using N pieces of received data obtained at every sampling time ⁇ t.
- data obtained at each sampling time is regarded as a virtual antenna, and a virtual antenna array (N virtual antennas arranged on the time axis ( It can be interpreted that the direction of arrival estimation is performed by constructing a temporal virtual array.
- FIG. 5 is a diagram illustrating an operation principle of the object detection device according to Embodiment 1 of the present invention.
- the receiver 1092 includes transmission signals 1003 (t 1 ), 1003 (t 2 ),..., 1003 (t N ) acquired from the transmission unit 1091 and received signals 1004 (t 1 ), 1004 (t 2 ),..., 1004 (t N ) are multiplied (mixed) to generate a received signal (complex baseband signal) r.
- the generated reception signal r is output to the signal processing unit 1095.
- the process of calculating the correlation matrix R from the received signal r, and further the process of calculating the evaluation functions of the equations (10) to (12) are performed by the signal processing unit 1095.
- the signal processing unit 1095 is included in the data processing unit 1093 in FIG.
- FIG. 6 is a diagram showing an operation principle when the beamformer method is applied to the object detection apparatus shown in FIG.
- the signal processing unit 1095 includes a phase shifter 1031 and an adder 1032. (In the following, the number of sampling points is changed from N to M.)
- the reflected waves 1007 (or their complex amplitudes) received by the virtual receiving antennas 1004 (t 1 ), 1004 (t 2 ),..., 1004 (t M ) are phase shifters 1031 (t 1 ), 1031. After receiving the phase rotations ⁇ 1 , ⁇ 2 ,..., ⁇ M at (t 2 ),..., 1031 (t M ), they are added by the adder 1032.
- the phase rotation by the phase shifters 1031 (t 1 ), 1031 (t 2 ),..., 1031 (t M ) and the addition by the adder 1032 are performed in the data processing unit 1093. It can be executed by processing, specifically by software processing using a processor.
- the position coordinates by the x-axis and the z-axis are set, the position of the transmitting unit 1091 is (0, 0), the position of the receiving unit 1092 is (x r , 0), and the position of the object 1001 is (x d , Z).
- the array factor AF (x d ) in the virtual array of the invention is calculated as the following equation (13).
- ⁇ ⁇ ⁇ t is the difference (frequency interval) of the carrier frequency f for each sampling.
- L t (x d ) is the distance between the transmission unit 1091 and the object 1001
- L r (x d ) is the distance between the reception unit 20 and the object 1001.
- c is the speed of light.
- FIG. 7 is a characteristic diagram showing an example of the antenna gain directivity of the object detection device according to the first exemplary embodiment of the present invention. Specifically, FIG. 7 shows the calculation result of the array factor AF (xd) of the virtual array using the above equations (2) and (3).
- the array factor (that is, the beam pattern) of the virtual array in such a case is shown.
- the phase rotation ⁇ m (m 1, 2,..., M) of the phase shifter 1031 (t m )
- the directivity (beam pattern) of the virtual array can be controlled.
- the beam width of the beam pattern can be calculated from the array factor AF (x d ) given by the above equations (2) and (3).
- the beam width is a factor that determines the direction of arrival estimation and the resolution of imaging (image).
- the beam width ⁇ x is given by the following equation (15).
- BW is the bandwidth of the RF carrier frequency as described above.
- BW ⁇ t ⁇ M.
- h (x r , x d , z) is a function of the position variable (x r , x d , z).
- xr xd
- h ( xr , xd , z) is given by [1+ (z / xr ) 2] 1/2.
- the beam width ⁇ x is reduced as the bandwidth BW is increased, and higher resolution performance can be obtained.
- the virtual array in the first embodiment may generate a virtual image due to the grating lobe.
- FIG. 8 is a diagram for explaining generation of a virtual image in a virtual array in one embodiment of the present invention.
- the phase amount ⁇ (x a ) is defined by the equation (16).
- the phase amount ⁇ (x a ) in the equation (16) is the phase shift of the radio wave from the transmission unit 1091 to the reception unit 1092 via the virtual image 1033 (position x a ) and the object 1001 ( This corresponds to the difference from the phase shift of the radio wave from the transmission unit 1091 to the reception unit 1092 via the position x d ).
- the phase phi (x a) is an integer multiple of 2 [pi, the same array factor is obtained in the position x d of the position x a and the object.
- ⁇ that is, a range of a position x satisfying the following conditional expression (17) can be used as a region where a virtual image is not generated (visible region).
- the visible region is expanded as the frequency interval ⁇ t is decreased, that is, as the sampling interval is decreased.
- the size (length) of the visible region is approximately inversely proportional to the frequency interval ⁇ ⁇ ⁇ t.
- the phase is controlled by the data processing unit 1093 for each measurement value for each sampling time of the baseband signal output from the receiving unit 1092. Then, by controlling the phase, the directivity of the effective antenna gain in the receiving unit 1092 is controlled, and further, the intensity distribution of the radio wave arriving at the receiving unit 1092 is measured by controlling the directivity of the antenna gain.
- the position and shape of the object 1001 can be detected. For this reason, it is not necessary to prepare a lot of receiving antennas and receivers as in the prior art.
- the first embodiment in detecting an object using radio waves, an increase in device cost, size, and weight can be suppressed while improving accuracy.
- FIG. 9 is a block diagram showing an example of a specific configuration of the object detection apparatus according to Embodiment 1 of the present invention.
- the object detection apparatus 1000 includes an output unit 1094 in addition to the transmission unit 1091, the reception unit 1092, and the data processing unit 1093.
- transmission section 1091 is configured by a transmitter
- reception section 1092 is configured by a receiver.
- the data processing unit 1093 is configured by a data processing device, that is, a computer (computer). This will be specifically described below.
- the transmission unit 1091 includes at least a power amplifier 1071, a coupler 1075, an oscillator 1103 having a frequency variable function, and a transmission control unit 1104.
- the oscillator 1103 outputs a transmission RF signal.
- the transmission RF signal output from the oscillator 1103 is amplified by the power amplifier 1071 and then transmitted as the transmission RF signal 1010 from the transmission antenna 1003.
- the transmission control unit 1104 controls the frequency of the RF signal output from the oscillator 1103.
- the RF signal output from the oscillator 1103 is output to the mixer 1042 in the receiving unit 1092 via the coupler 1075.
- the RF signal output to the mixer 1042 via the coupler 1075 is used as the LO signal of the receiving unit 1092.
- the reception unit 1092 includes a low-noise amplifier 1041, a mixer 1042, a filter 1043, an analog-digital converter 1044, and a reception control unit 1102 in addition to the reception antenna 1004. Yes.
- the receiving unit 1092 receives the radio wave (RF signal) 1007 reflected from the object 1001 by the receiving antenna 1004 provided therein.
- An RF signal 1007 received by the receiving antenna 1004 is amplified by the low noise amplifier 1041 and then input to the mixer 1042.
- the mixer 1042 mixes the received RF signal amplified by the low-noise amplifier 1041 and the RF signal (received LO signal) output from the transmission unit 1091 via the coupler 1075 to generate a baseband signal.
- a frequency signal (IF signal) is generated and output to the filter 1043.
- the filter 1043 removes noise from the baseband signal and inputs the baseband signal from which the noise has been removed to the analog-digital converter 1044.
- Analog-digital converter 1044 converts a baseband signal, which is an analog signal, into a digital baseband signal, and inputs the obtained digital baseband signal to reception control section 1102.
- the digital baseband signal obtained above corresponds to the in-phase component (In-phase signal) I (t) described in Equation (6).
- the reception control unit 1102 multiplies the in-phase component I (t) by the Hilbert transform to generate a quadrature component (Quadrature signal) Q (t). Furthermore, reception control section 1102 generates complex baseband signal r (t) from in-phase component I (t) and quadrature component Q (t) according to equation (8). The generated complex baseband signal r (t) is transferred to the data processing unit 1093. As described above, a quadrature modulator may be used instead of the mixer 1042 to generate the quadrature component Q (t).
- the data processing unit 1093 performs the processing described with reference to FIGS. 2 to 8, that is, the estimation processing of the arrival direction of the received radio wave 1007 on the delivered complex baseband signal r (t). Furthermore, the data processing unit 1106 also executes imaging processing (image generation) of the object 1001. After that, the data processing unit 1106 outputs the processing result, that is, the estimated arrival direction and the generated image to the output unit 1094.
- the output unit 1094 is a display device, for example, and displays the processing result on the screen.
- the object detection apparatus 1000 may include a plurality of transmission units 1091 and reception units 1092. Further, the data processing unit 1093 and the output unit 1094 may be incorporated in the transmission unit 1091 or the reception unit 1092.
- FIG. 7 is a flowchart showing the operation of the object detection apparatus 100 according to Embodiment 1 of the present invention.
- FIGS. 1 to 8 are referred to as appropriate.
- an object detection method is implemented by operating an object detection apparatus. Therefore, the description of the object detection method in the first embodiment is replaced with the following description of the operation of the object detection apparatus 1000.
- the transmission control unit 1104 specifies the current sampling time t m and calculates the frequency (f min + ⁇ t m ) of the RF signal transmitted by the transmission antenna 1003. (Step A1).
- the transmission control unit 1104 generates a control signal for the oscillator 1103 so that an RF signal having a frequency (f min + ⁇ t m ) is transmitted from the transmission antenna 1003, and outputs the control signal from the transmission antenna 1003.
- the RF signal having the frequency (f min + ⁇ t m ) is transmitted (step A2).
- the transmission control unit 1104 sends a control signal to the oscillator 1103 so that the output frequency of the oscillator 1103 becomes (f min + ⁇ t m ), and the oscillator 1103 has a carrier frequency of (f min + ⁇ t m ).
- RF signal is output.
- the RF signal is amplified by the power amplifier 1071 and transmitted from the transmission antenna 1003.
- the RF signal output from the oscillator 1103 is also sent to the mixer 1042 in the receiving unit 1092 via the coupler 1075.
- the receiving antenna 1004 receives the radio wave (RF signal) 1007 reflected from the object 1001 (step A 3).
- the reception control unit 1102 calculates a complex baseband signal r (t) from the in-phase component I (t) of the baseband signal obtained from the received RF signal (step A4).
- Step A4 the RF signal 1007 received by the receiving antenna 1004 is amplified by the low noise amplifier 1041, and then input to the mixer 1042.
- the mixer 1042 mixes the received RF signal amplified by the low-noise amplifier 1041 with the RF signal output from the transmission unit 1091 via the coupler 1075 as an LO signal, and generates a baseband signal (in-phase component I (t)). Is generated.
- the baseband signal (in-phase component I (t)) is input to the analog-digital converter 1044 via the filter 1043, where it is converted into a digital signal.
- the reception control unit 1102 calculates a complex baseband signal r (t) from the digitally converted baseband signal (in-phase component I (t)).
- the data processing unit 1093 estimates the arrival direction of the received radio wave 1007 using the complex baseband signal r (t), and further executes the imaging processing of the object 1001 using the estimation result ( Step A5).
- steps A1 to A5 are repeatedly executed, and the result of the repeated processing is displayed on the screen by the output unit 1094.
- the FM-CW system is adopted as a radio wave transmission and reception system. For this reason, it is not necessary to provide an oscillator in the receiving unit 1092, and the apparatus cost can be reduced also in this respect. Further, since the receiving unit 1092 does not require an oscillator, there is no need to synchronize the oscillator 1103 in the transmitting unit 1091 and the oscillator in the receiving unit 1092. As a result, a synchronization error between the transmitting unit 1091 and the receiving unit 1092 and There will be no degradation in detection accuracy.
- the process performed in the first embodiment is a process for estimating the position (particularly one-dimensional direction) of the object 1001 in the second embodiment, and the arrangement state and shape of the object 1001 in the third embodiment in a two-dimensional image. Used for display processing. These processes are also performed in the data processing unit 1093.
- Embodiment 2 Next, an object detection apparatus and an object detection method according to Embodiment 2 of the present invention will be described with reference to FIGS.
- Embodiment 2 shows an example in which the object detection apparatus 1000 shown in Embodiment 1 is used to estimate the position of an object, particularly a one-dimensional direction. Therefore, also in the second embodiment, the object detection apparatus includes the transmission unit 1091, the reception unit 1092, and the data processing unit 1093 shown in FIGS. 1 and 9. However, the second embodiment is different from the first embodiment in the number of receiving units 1092. This will be specifically described below.
- FIG. 11 is a diagram showing the configuration and operation principle of the object detection device according to the second embodiment of the present invention.
- the object detection device includes N reception units for one transmission unit.
- the object detection apparatus includes a single transmit antenna 1003, receiving antenna 1004 1 the N, ⁇ ⁇ ⁇ , 1004 n, ⁇ ⁇ ⁇ , 1004 N And.
- receiving antenna 1004 when a specific reception antenna is not indicated, it is expressed as “reception antenna 1004”.
- each receiving antenna is installed along one direction with reference to the transmitting antenna.
- the position of the transmission antenna 1003 is (d 0 , 0) in (x, z) coordinates.
- the positions of the N receiving antennas 1004 are (d x1 , 0), (d x2 , 0), ..., (d xN , 0), respectively.
- the object detection apparatus can operate even when the number N of reception antennas is 1, which is the minimum.
- N receiving antennas in order to give generality to the theory, the case of N receiving antennas is treated here.
- the positions of the transmission antenna 1003, the reception antenna 1004, and the object 1001 are fixed at the above positions.
- the data processing unit estimates the arrival direction of the radio wave from the measured value of the baseband signal received by each receiving antenna 1004. Further, the data processing unit identifies the intensity distribution of the radio wave based on the estimated arrival direction of the radio wave, and detects the position in one direction of the object 1001 based on the identified intensity distribution.
- the data processing unit constructs a time virtual array from the measured values of the baseband signal for each sampling time, and calculates a correlation matrix of the time virtual array. More specifically, the data processing unit constructs a sub-array of a time virtual array from measured values of baseband signals having different sampling times, calculates a correlation matrix for each sub-array, and calculates an average value of the correlation matrix for each sub-array. calculate. Then, the data processing unit obtains an evaluation function that reflects the position of the object 1001 based on the average value of the correlation matrix, and generates an image of the object 1001 from the obtained evaluation function.
- the object detection device will be specifically described.
- an FM-CW signal is transmitted from the transmitting antenna 1003 as in the first embodiment.
- the receiving antenna 1004 receives the reflected wave 1007 from the object 1001.
- the subscript “xn” means a signal received by the nth receiving antenna 1004 n arranged in the x-axis direction.
- the complex amplitude s xn (x d , t m ) of 1007 is an unknown number.
- the complex amplitude of the signal actually measured by the receiving antenna 1004 n is s xn ′ (t m )
- ⁇ (x d ) is an unknown number representing the reflectance of the object 1001 d .
- Exponential term of equation (21) in the right side represents the wave of the phase shift caused by the path from the transmitting antenna 1003 through object 1001 d up to the receiving antenna 1004 n.
- the size of the matrix A is MN ⁇ D
- matrix A size of n is M ⁇ D
- the size of the vector a n (x d) is the M ⁇ 1.
- the size of the matrix is expressed by the number of vertical and horizontal elements.
- the desired signal vector s is defined by the following equation (25) using the variable s 0 and ⁇ (x d ) in the right side of equation (11).
- the second embodiment it is an object to determine an evaluation function reflecting the x d dependency (that is, ⁇ (x d )) of the desired signal vector s by measurement by the receiving antenna 1004.
- the distribution and shape of the object 1001 are detected from the xd dependency of the desired signal vector s.
- the relationship of the above equation (22) can be expressed as the following equation (26) using the measurement signal vector s x , the direction matrix A, and the desired signal vector s.
- the matrix A is required to be full rank as an application condition of the MUSIC method. Adding the noise vector n (t) has an effect of effectively destroying the dependency of the column vector or row vector in the matrix A and bringing the matrix A closer to the full rank.
- the measurement signal vector s x (t) defined by Expression (23) is received by the reception antenna 1004. Then, the data processing unit calculates a correlation matrix R x represented by the following equation (27) using the received measurement signal vector s x .
- E [] in the equation (27) represents an average over the number of points (snapshot number) at time t defining the noise (random number) vector n (t).
- PN noise power
- I is a unit matrix of MN ⁇ MN order.
- H represents complex conjugate transpose.
- the sizes of the correlation matrix R x , the matrix A, and the matrix S are MN ⁇ MN order, MN ⁇ D order, and D ⁇ D order, respectively.
- the matrix A and the matrix S in the equation (28) are required to be full rank.
- Full rank means that the rank of the matrix matches the size of the matrix (the smaller of the number of rows or columns) and that all row and column vectors in the matrix are all linearly independent. Is done.
- each column vector is a function of the position xd , and each column vector is independent and has a full rank.
- the row vectors of the i-th row and the j-th row of the matrix S have the same value and are linearly dependent.
- the number of floors is lowered by 1 and is not full rank.
- a virtual array is constructed by regarding one frequency as one antenna.
- FIG. 12 is a diagram for explaining the concept of a subarray used in the object detection device according to Embodiment 2 of the present invention.
- the entire array is composed of M 0 frequency measurement data
- the measurement signal vector s xq of the subarray q in Expression (29) is given by the following Expression (30) between the direction matrix A in Expression (24) and the desired signal vector s in Expression (14). There is a relationship.
- the correlation matrix R xq of the subarray q is calculated as in the following formula (31).
- the sizes of the correlation matrix R xq , the matrix A ′, and the matrix S ′ are NM ⁇ NM order, NM ⁇ ND order, and ND ⁇ ND order, respectively.
- the relationship between the correlation matrix R x ′ of all subarray averages and the direction matrix A is calculated as in the following equation (32).
- the correlation matrix R x ′ in the equation (32) has the form A ′S ′′ A ′ H like the correlation matrix in the equation (17). Therefore, if the matrices A ′ and S ′′ are full rank,
- the evaluation function P MU (x) reflecting the x dependency (that is, ⁇ (x)) of the intensity of the desired signal vector s can be calculated.
- the matrix A ′ given by the equation (31) is also full rank.
- the row vectors of the u-th row and the v-th row of the matrix C have the same value and are linearly dependent. Disappear.
- b id is a function of the distances L 0 (x d ) and L x (x d ), and if the position x d is different, these distances have different values.
- b iu b iv (u ⁇ v) is not satisfied, and C i becomes a full rank.
- the matrix S in the equation (28) does not have a full rank because the reflectance ⁇ (x d ) can take the same value even if the position x d is different.
- the matrix S ′′ is guaranteed to be full rank because the distances L 0 (x d ) and L x (x d ) always change if the position x d changes.
- the rank of S ′′ becomes Q, and the rank of S ′′ increases by 1 whenever the number Q of subarrays is increased by one.
- each subarray is an independent signal set, and the number of subarrays is increased by one to increase the number of independent signal sets by one, so that the rank of the matrix S ′′ is also increased by one.
- the arrival direction estimation is performed by applying the MUSIC method to the correlation matrix of a general array antenna.
- the MUSIC method is applied to the correlation matrix R x ′ of all subarray averages calculated by the equation (21) (in the same manner as that applied to a general array antenna formally).
- the evaluation function P MU (x) reflecting the x dependency (that is, ⁇ (x)) of the intensity of the desired signal vector s is calculated.
- the evaluation function P MU (x) is given by the following equation (35).
- a (x) is a column vector of the direction matrix A defined by the equation (34).
- E N is given by the following equation (36).
- eigenvectors ⁇ e D + 1 , e D + 2 ,..., E MN ⁇ of (MN ⁇ D) noise spaces are used, but since at least one is required, MN ⁇ It is necessary to satisfy D ⁇ 1, that is, MN ⁇ D + 1.
- the beamformer method, the Capon method the linear method (described in Non-Patent Document 1 in the same manner as that applied to a general array antenna formally).
- the prediction method it is possible to calculate an evaluation function reflecting the x dependency (that is, ⁇ (x)) of the intensity of the desired signal vector s (t).
- the evaluation function P BF (x) based on the beam former method in the second embodiment is given by the following equation (37).
- the processing disclosed in the above-described second embodiment that is, the processing for calculating the evaluation function from the measurement result of the reflected wave and determining the position of the object from the evaluation function is executed by the data processing unit 1093 shown in FIG. Is done.
- the control by the phase shifter 1031 and the adder 1032 in the first embodiment is performed, and the received signal strength is maximized. This corresponds to the process of searching for the beam direction.
- the second embodiment it is possible to detect only the position information x d (that is, the position in the one-dimensional direction) of the coordinates (that is, the x axis) in the direction connecting the transmission unit and the reception unit.
- the object detection apparatus including the transmission unit and the reception unit has rotational symmetry about the x axis, and thus cannot be distinguished even if the coordinate value of the object 1001 other than the x axis is different. is there.
- a method for detecting position information of coordinates other than the X axis will be described later in a third embodiment.
- FIG. 13 is a flowchart showing the operation of the object detection apparatus according to the second embodiment of the present invention. Also in the second embodiment, the object detection method is implemented by operating the object detection device. Therefore, the description of the object detection method in the second embodiment is replaced with the following description of the operation of the object detection apparatus 1000.
- the transmission unit radiates radio waves while changing the frequency toward the object (step B1).
- each of the plurality of receiving units receives the reflected wave of each frequency from the object by the corresponding receiving antenna (step B2).
- Each receiving antenna is arranged in one direction as viewed from the transmitting unit.
- the data processing unit calculates the position of the object from the peak of the evaluation function (step B6).
- the calculation result is output to the output unit.
- the second embodiment it is possible to estimate a one-dimensional direction of an object without preparing a large number of receiving antennas. Also in the second embodiment, the effects described in the first embodiment can be obtained.
- Embodiment 3 shows an example in which a two-dimensional image for identifying the arrangement and shape of an object is generated based on the concept of a virtual array by the object detection apparatus 1000 shown in Embodiment 1. Therefore, also in the third embodiment, the object detection apparatus includes the transmission unit 1091, the reception unit 1092, and the data processing unit 1093 illustrated in FIGS. 1 and 9. However, the third embodiment is different from the first embodiment in the number of receiving units 1092. This will be specifically described below.
- FIG. 14 is a diagram showing the configuration and operating principle of the object detection device according to the third embodiment of the present invention.
- FIG. 14 shows the positional relationship between each antenna and the object.
- the data processing unit calculates the product of the baseband signals generated by each of the plurality of receiving units, and based on the calculated product, determines the position of the object 1001 in the N-dimensional coordinate space with the N direction as the coordinate axis. Detect.
- the transmission antenna 1003 is installed at the origin of the coordinates, and the reception antenna 1004 (x) and the reception antenna 1004 (y) of the reception unit are installed on the x-axis and the y-axis, respectively. .
- N 2.
- the direction connecting the transmitting antenna 1003 and the receiving antenna 1004 (x) and the direction connecting the transmitting antenna 1003 and the receiving antenna 1004 (y) are different from each other (not parallel). This is a desirable mode for obtaining a two-dimensional image. Note that the direction connecting the transmitting antenna 1003 and the receiving antenna 1004 (x) and the direction connecting the transmitting antenna 1003 and the receiving antenna 1004 (y) are not necessarily orthogonal to each other.
- the RF signal (radio wave) 1010 is emitted from the transmission antenna 1003 toward the object 1001 existing on the focal plane 1002. After the object 1001 is irradiated with the RF signal 1010, the reflected wave 1007 (x) and the reflected wave 1007 (y) from the object 1001 are received at the reception antenna 1004 (x) and the reception antenna 1004 (y), respectively. Received. Also in the third embodiment, as in the first and second embodiments, the carrier frequency of the RF signal 1010 output from the transmission antenna 1003 continuously changes with time.
- the third embodiment shown in FIG. 14 is intended to replace the two array antennas 201 in the Mills-cross method shown in FIG. 26 with frequency virtual arrays, respectively.
- the two array antennas 201 in the Mills-cross method shown in FIG. 27 include a virtual array composed of a combination of a transmission antenna 1003 and a reception antenna 1004 (x), a transmission antenna 1003 and a reception antenna 1004 ( It is replaced with a virtual array composed of a pair with y). Therefore, in the third embodiment, the minimum number of antennas necessary for generating a two-dimensional image is three.
- FIGS. 15 and 16 are explanatory diagrams for explaining a method of calculating a correlation matrix of a two-dimensional frequency virtual array according to Embodiment 3 of the present invention.
- 15 and 16 show calculation models for analyzing the operation of generating a two-dimensional image.
- one transmitting antenna 1003 (x 0 ) and N receiving antennas 1004 (x 1 ),. (X N ) is installed. Furthermore, in the calculation model of the third embodiment, one transmission antenna 1003 (y 0 ) and N reception antennas 1004 (y 1 ),..., 1004 (y N ) are also installed on the y axis. Has been.
- the position of the transmitting antenna 1003 (x 0 ) on the x-axis is (dx 0 , 0 , 0), and the position of the n-th receiving antenna 1004 (x n ) is (dx n , 0, 0). To do. Further, the position of the transmitting antenna 1003 (y 0 ) on the y-axis is (0, dy 0 , 0), and the position of the nth receiving antenna 1004 (y n ) is (0, dy n , 0).
- the positional relationship between the object detection device (the transmitting antenna 1003 and the receiving antenna 1004) and the object 1001 is fixed to the above positional relationship.
- the transmission antenna 1003 (x 0 ) and the transmission antenna 1003 (y 0 ) are separately arranged for the x-axis and the y-axis. This is because it has sex.
- the transmission antenna 1003 (x 0 ) and the transmission antenna 1003 (y 0 ) may be configured by a single transmission antenna, in which case this single transmission antenna transmits.
- the reception antenna on the x-axis and the reception antenna on the y-axis may be simultaneously received.
- the transmission antenna 1003 (x 0 ) and the transmission antenna 1003 (y 0 ) have M carrier frequencies ⁇ t 1 , ⁇ t 2 ,. , ⁇ t M RF signal 1010 is transmitted.
- the modulation of the RF signal 1010 is also performed in the third embodiment by the above-described FM-CW method.
- a distance L xo (x d , y d ) between the object 1001 d and the transmission antenna 1003 (x 0 ) on the x axis is given by the following equation (42).
- a distance L xn (x d , y d ) between the object 1001 d and the n-th receiving antenna 1004 (x 0 ) on the x axis is given by the following equation (43).
- the distances from the object 1001 d are similarly expressed as L yo (x d , y d ) and L yn (x d , If y d ), they are given by the following equations (44) and (45).
- a measurement signal vector s x is defined.
- the product is obtained for all combinations of the elements of the x-axis direction measurement vector s x of the above equation (50) and the elements of the y-axis direction measurement vector s y of the above equation (51). Is calculated, a direct product vector s xy shown in the following equation (52) is generated.
- the “product” here corresponds to the “product of baseband signals” described above.
- n and v are antenna numbers arranged in the x and y directions, respectively, and m and w are subscripts representing frequency numbers of signals received by antennas arranged in the x and y directions, respectively.
- the direction matrix A is defined by the following equation (53).
- the size of the directional matrix A is (MN) 2 ⁇ D
- the size of the matrix A nv is M2 ⁇ D
- the size of the vector anv (x d , y d ) is M 2 ⁇ 1.
- the matrix A nv is a directional matrix in which the n-th x-direction antenna 1004 (x n ) and the v-th y-direction antenna 1004 (y v ) are involved.
- the direction matrix A of the entire system is a collection of the direction matrices A nv of all the antenna number pairs (n, v).
- the desired signal vector s is defined by the following equation (54) using the complex amplitude s 0 and the reflectance ⁇ (x d , y d ).
- Equation (48) and (49) From Equations (48) and (49), between the measured signal vector s xy (t) in Equation (52), the direction matrix A in Equation (53), and the desired signal vector s in Equation (54), The following relational expression (55) is obtained.
- Expression (55) a vector n (t) having noise (random number) as an element is added.
- the correlation matrix R xy is calculated using the measurement signal vector s xy of Expression (52) obtained by the measurement. From the relationship of the equation (55), the relationship between the correlation matrix R xy and the direction matrix A is given by the following equation (56).
- a P N is the average power of the noise term n (t)
- I is (MN) 2 ⁇ (MN) 2-order unit matrix.
- the sizes of the correlation matrix R xy , the matrix A, and the matrix S are (MN) 2 ⁇ (MN) second order, (MN) 2 ⁇ D order, and D ⁇ D order, respectively.
- Expressions (55) and (56) are the same types as Expressions (26) and (28) in the one-dimensional direction-of-arrival estimation discussed in the second embodiment. Therefore, the evaluation function P MU (x, y) reflecting ⁇ (x d , y d ) can be calculated by applying the MUSIC method to the correlation matrix R xy in the same procedure as the one-dimensional arrival direction estimation.
- the matrix A and the matrix S in Expression (56) are full rank as an application condition of the MUSIC method.
- one subarray is constructed with M frequencies in the same procedure as the subarray method in the one-dimensional direction-of-arrival estimation described in the second embodiment, and a total of Q subarrays are constructed.
- the qth subarray signal is defined by the following equation (57).
- the qth subarray signal is obtained by shifting the subscripts m and w representing the sampling time of the component s xy (nv) (mw) of the signal vector s xy at the same time by + (q ⁇ 1).
- Correlation matrix R x q subarray q is calculated as the following equation (59).
- Equation (59) the sizes of the correlation matrix R xy q , the matrix A ′, and the matrix S ′ are (NM) 2 ⁇ (NM) second order, (NM) 2 ⁇ N 2 D order, N 2 D ⁇ N, respectively. 2 D order.
- the relationship between the correlation matrix R xy ′ of all subarray averages and the direction matrix A ′ is calculated as in the following equation (60).
- the evaluation function P MU (x, y) reflecting ⁇ (x d , y d ) is calculated by applying the MUSIC method to the correlation matrix R xy ′. it can.
- the directional matrices A 11 , A 12 ,..., A 1N ,..., A N1 ,. A 'given by) is also full rank.
- the matrix S ′′ is full rank if Q ⁇ D.
- the application condition MN ⁇ D + 1 of the MUSIC method in the one-dimensional direction-of-arrival estimation is (MN) 2 ⁇ D + 1 in the two-dimensional image generation.
- Q M 0 -M + 1 and Q ⁇ D in the subarray
- the condition of the number of sampling times (number of frequencies) M 0 required is given by the following equation (61). That is, the number of sampling time M 0 required is generally increases in proportion to the number D of to be detected position.
- the evaluation function P MU (x, y) reflecting ⁇ (x d , y d ) is calculated by applying the MUSIC method to the correlation matrix R xy ′ of all subarray averages calculated by the equation (60). To do. As a result, the evaluation function shown in the following formula (62) is obtained.
- a (x, y) is a column vector of the direction matrix A defined by the equation (42).
- E N is given by the following equation (63).
- the evaluation function of each method can also be calculated by applying the beamformer method, the Capon method, and the linear prediction method (described in Non-Patent Document 1) in the same manner as described above.
- the evaluation function P BF (x, y) based on the beam former method in the third embodiment is given by the following equation (64).
- the above-mentioned evaluation functions P BF (x, y), P CP (x, y), and P LP (x, y) are also the object 1001 d as with the evaluation function P MU (x, y) obtained by the MUSIC method.
- the process disclosed in the third embodiment that is, the process of calculating the evaluation function from the measurement result of the reflected wave and determining the position of the object from the evaluation function is also shown in FIG. 9 as in the second embodiment.
- the data processing unit 1093 executes.
- the control by the phase shifter 1031 and the adder 1032 in the first embodiment is performed, and the received signal strength is maximized. This corresponds to the process of searching for the beam direction.
- FIG. 17 is a flowchart showing the operation of the object detection apparatus according to the third embodiment of the present invention. Also in the third embodiment, the object detection method is implemented by operating the object detection device. Therefore, the description of the object detection method in the third embodiment is replaced with the following description of the operation of the object detection apparatus 1000.
- the transmission unit irradiates a target with a radio wave while changing the RF carrier frequency (step C1).
- each of the plurality of receiving units receives the reflected wave from the object by the corresponding receiving antenna (step C2).
- Each receiving antenna is arranged in two directions as viewed from the transmitting unit.
- the data processing unit calculates the position of the object from the peak of the evaluation function, and further outputs the arrangement and shape of the object to the output unit as a two-dimensional image (step C6).
- FIG. 18 is a diagram illustrating an example of an image output from the object detection device according to Embodiment 3 of the present invention.
- the object 1001 has (x, y, z) coordinate display, ( ⁇ 20 cm, ⁇ 20 cm, 100 cm), (0 cm, 0 cm, 100 cm), and (20 cm, 20 cm, 100 cm). Is arranged.
- the transmitting antenna 1003 is arranged at a position of ( ⁇ 100 cm, ⁇ 100 cm, 0 cm). It is assumed that the receiving antenna 1004 is disposed at a position (0 cm, ⁇ 100 cm, 0 cm) and a position ( ⁇ 100 cm, 0 cm, 0 cm).
- 1 chirp period (T chirp) sampling time in the number (the number of all frequencies) M 0 is 21, the number Q of sub-arrays 10, one number per sub-array (the number of frequencies) M is 12 .
- the present invention can also be applied when installed along the above.
- the receiving antenna is installed along three orthogonal directions, the position of the object in the three-dimensional space can be specified.
- FIG. 19 is a diagram illustrating a schematic configuration of the object detection device according to the fourth embodiment of the present invention.
- P indicates the number of object detection units 1202.
- the transmission unit 1091 p irradiates the objects 1201 p1 , 1201 p2 ,..., 1201 pQ with radio waves, and the reception unit 1092 p performs the object 1201 p1 , 1201 p2 ,. , 1201 receives a reflected wave from pQ .
- the states of the objects 1201 p1 , 1201 p2 ,..., 1201 pQ are detected.
- Q is the number of objects 1201.
- the object detection apparatus 1200 when the objects 1201 p1 , 1201 p2 ,..., 1201 pQ are people, clothes worn by people (1201 p1 , 1201 p2 ,..., 1201 pQ ). The presence of the article under the clothes can be detected by the radio wave transmitted through.
- the object detection apparatus 1200 detects the objects (1201 p1 , 1201 p2 ,..., 1201 pQ ) when the objects 1201 p1 , 1201 p2 ,..., 1201 pQ are objects (particularly dielectrics).
- the internal structure of the object (1201 p1 , 1201 p2 ,..., 1201 pQ ) can be detected by the transmitted radio waves.
- the object detection device 1200 can also sequentially detect the states of the objects 1201 p1 , 1201 p2 ,..., 1201 pQ by the object detection unit 1202 p when the object is the target of the flow work.
- one object detection unit 1202 p is assigned to detection or inspection of one object 1201.
- the fourth embodiment is not limited to this, with respect to the detection or inspection of one object 1201 may be assigned a plurality of object detection unit 1202 p. Further, in the fourth embodiment, to detect or inspection of a plurality of objects 1201, one object detection unit 1202 p may also be split hit.
- FIG. 20 is a block diagram specifically showing the configuration of the object detection apparatus according to Embodiment 4 of the present invention.
- the object detection apparatus 1200 of the fourth embodiment in addition to the plurality of object detection unit 1200 p, and a data control unit 1203.
- the data control unit 1203 performs the same processing as the data control unit 1093 in the first to third embodiments on each object detection unit 1202 p .
- FIG. 21 is a diagram showing an example of frequency control performed in the fourth embodiment of the present invention.
- the sampling time t is the same as in the first to third embodiments. 1 , t 2 ,..., T M , the carrier frequency (RF frequency) f continuously changes from f min to f max . Specifically, in each object detection unit 1202 p, time variation of the RF frequency is controlled to chirped.
- the data control unit 1203 controls so that the time variation of the RF frequency of the object detection unit 1202 p are shifted from each other.
- different object detection units 1202 p and object detection units 1202 r do not operate at the same RF frequency.
- each object detection unit 1202 p is controlled so that the RF frequency in each object detection unit 1202 p is different while functioning as the above-described data processing unit. .
- the data control unit 1203 executes steps A1 to A5 shown in FIG. 10, and at that time, in step A2, different ⁇ t m is set for each object detection unit.
- the synthetic aperture radar system When comparing the synthetic aperture radar system and the present embodiment, the synthetic aperture radar system has a problem that the receiver needs to be moved mechanically, which increases the time required for detecting and inspecting the object. On the other hand, in this embodiment, it is only necessary to electronically scan the reception frequency instead of the position of the receiver, so that the time for detecting and inspecting an object can be shortened as compared with the synthetic aperture radar system.
- the object detection device and the object detection method in the present embodiment the number of necessary antennas and the associated receivers can be reduced as compared with a general array antenna system, so the cost, size, There is an effect that the weight can be reduced.
- the object detection apparatus and the object detection method in the present embodiment do not require the apparatus to be moved mechanically, and thus the object detection and inspection time can be shortened. Play.
- an image of a detection target object is obtained by irradiating a detection target with radio waves having different RF frequencies at each sampling time, and detecting a radio wave reflected by the target object or a radio wave radiated from the target object. Can be generated. Therefore, according to the present embodiment, it is possible to realize image generation by high-speed scanning without reducing the number of antennas and receiving units required compared to the prior art and without having to move them.
- the present invention it is possible to suppress an increase in device cost, size, and weight while improving accuracy in detecting an object using radio waves.
- INDUSTRIAL APPLICABILITY The present invention is useful for imaging and inspecting articles hidden under clothes or articles in bags.
- Object detection device 1001, 1201 Object (object to be detected) 1002 Focal plane 1003 Transmitting antenna 1004 Receiving antenna 1007, 1010 Radio wave (RF signal) 1041 Low noise amplifier 1042 Mixer 1043 Filter 1044 Analog-digital converter 1075 Coupler 1091 Transmitter 1092 Receiver 1093 Data receiver 1094 Output unit 1103 Oscillator 1102 Reception controller 1104 Transmission controller 1202 Object detection unit 1200 Object detection device (implementation) Form 4) 1203 Data control unit
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
物体検知装置1000は、電波によって物体1001を検知するための装置である。物体検知装置1000は、時間の経過と共に周波数が連続的に変化する電波を送信信号として放射する送信部1091と、物体1001からの電波を受信信号として受信し、更に受信した受信信号に送信部1091から取得した送信信号をミキシングしてベースバンド信号を生成する受信部1092と、サンプリング時間毎のベースバンド信号の測定値から電波の到来方向を推定し、推定したサンプリング時間毎の到来方向に基づいて電波の強度分布を特定し、特定した強度分布に基づいて物体1001を検知する、データ処理部1093とを備えている。
Description
本発明は、対象物で反射又は対象物から放射された電波から対象物を検知するための、物体検知装置、及び物体検知方法に関する。
電波(マイクロ波、ミリ波、テラヘルツ波など)は、光と異なり、物体を透過する能力に優れている。この電波の透過能力を活用して、衣服の下に隠されている物品又は鞄の中の物品等を画像化して検査するイメージング装置(物体検知装置)が実用化されている。同様にして、衛星又は航空機から雲を透過して地表を画像化するリモートセンシング技術も実用化されている。
また、物体検知装置における画像化の方式としては、いくつかの方式が提案されている。一つは、アレイアンテナ方式である(例えば、非特許文献1参照)。ここで、図22及び図24を用いて、アレイアンテナ方式について説明する。図22は、従来からのアレイアンテナ方式を採用した物体検知装置を示す図である。図23は、図22に示された受信機の構成を示す図である。
図22に示すように、アレイアンテナ方式においては、物体検知装置は、送信機211と受信機201とを備えている。また、送信機211は、送信アンテナ212を備えている。受信機201は、受信アンテナ2011、2022、・・・、202Nを備えている(Nは受信アンテナの数)。
送信機211は、送信アンテナ212から、RF信号(電波)213を検知対象物2041、2042、・・・、204K(Kは対象物の数)に向けて照射する。RF信号(電波)213は、検知対象物2041、2042、・・・、204Kにおいて反射され、反射波2031、2032、・・・、203Kがそれぞれ発生する。
発生した反射波2031、2032、・・・、203Lは、受信アンテナ2011、2022、・・・、202Nにおいて受信される。受信機201は、受信した反射波2031、2032、・・・、203Kに基づいて、検知対象物2041、2042、・・・、204Kで反射された電波の電波強度を算出する。その後、受信機201は、算出した電波強度の分布を画像化する。これにより、検知対象物2041、2042、・・・、204Kそれぞれの像が得られることになる。
また、図23に示すように、アレイアンテナ方式が採用される場合、受信機201は、N本の受信アンテナ2021、2022、・・・、202Nを備えている。受信アンテナ2021、2022、・・・、202Nで、角度θk(k=1,2,・・・K)を持つK個の到来波2081、2082、・・・、208Kを受信する。到来波2081、2082、・・・、208Kの複素振幅を[s(θ1), s(θ2),・・・, s(θK)]とする。受信機201はダウンコンバータ(図23では非図示)を備えていて、前記ダウンコンバータは各受信アンテナ2021、2022、・・・、202Nで受信したRF信号の複素振幅(ベースバンド信号)[r1,r2,・・・,rN]を抽出する。受信アンテナ2021、2022、・・・、202Nで受信した信号の複素振幅[r1,r2,・・・,rN]は信号処理部205へ出力される。
受信アンテナ2021、2022、・・・、202Nにおける、受信信号の複素振幅[r1,r2,・・・,rN]と、到来波の複素振幅[s(θ1), s(θ2),・・・, s(θK)]との関係は、以下の式(1)で与えられる。
上記式(1)において、n(t)はノイズ成分を要素とするベクトルである。添字Tはベクトルないし行列の転置を表す。dはアンテナ間の距離、λは到来波(RF信号)2081、2082、・・・、208Kの波長である。
また、上記式(1)において、受信信号の複素振幅rは測定で得られる量である。方向行列Aは信号処理上で定義(指定)できる量である。到来波の複素振幅sは未知数であり、測定で得た受信信号rから到来波sの方向を決定する事が到来波方向推定の目的となる。
到来方向推定のアルゴリズムでは、測定で得た受信信号rから相関行列R=E[r・rH]を計算する。ここでE[]は括弧内の要素に時間平均の処理を施す事を表し、添字Hは複素共役転置を表す。次に、計算した相関行列Rから、以下の式(2)~(4)で示すいずれかの評価関数が計算される。
MUSIC法におけるEN=[eK+1,・・・,eN]は、相関行列Rの固有ベクトルの内、固有値がノイズn(t)の電力となるN-(K+1)個のベクトルで構成した行列である。
また、図23で示した従来型のアンテナアレイにおいて、受信信号rから相関行列Rを計算する過程、更には、式(2)~(4)の評価関数を計算する過程は、信号処理部205で実施される。
非特許文献1に記載の理論によれば、式(2)~式(4)で示した評価関数は、到来波の角度θ1,θ2,・・・,θKにおいてピークを持つ。従って、評価関数を計算してそのピークを見れば、到来波の角度を求める事ができる。式(2)~式(4)の評価関数で得た到来波の角度分布から、対象物の位置や形状を画像として表示する事ができる。
式(A2)~式(A4)で示した評価関数の内、特に式(2)のビームフォーマ法を適用する場合の信号処理部は、図24において示される。図24は、図22に示された受信機においてビームフォーマ法が適用される場合の例を示す図である。
図24で示した従来型のアンテナアレイの移相器2061、2062、・・・、206Nと合成器207が、図23で示した従来型のアンテナアレイにおける信号処理部205に対応する。移相器2061、2062、・・・、206Nは、それぞれ、受信アンテナ2021、2022、・・・、202Nで受信した到来波の複素振幅2081、2082、・・・、208Nに対し、位相回転Φ1、Φ2、・・・、ΦNを加える。位相回転Φ1、Φ2、・・・、ΦNが加えられた到来波2081、2082、・・・、208Nは、加算器207で加算される。
移相器2061、2062、・・・、206Nと、加算器207とは、アナログ回路によって実装される事もあれば、コンピュータに組み込まれたソフトウェアによって実装される事もある。また、アレイアンテナ方式では、移相器2061、2062、・・・、206Nにおける、位相回転Φ1、Φ2、・・・、ΦNの設定により、アレイアンテナの指向性が制御される。受信アンテナ202の指向性をg(θ)とし、受信アンテナ202nで受信した到来波208n(n=1,2,・・・,N)の振幅と位相とをそれぞれanおよびφnとした場合、アレイアンテナの指向性E(θ)は、以下の式(5)のように計算される。
式(5)において、アレイアンテナの指向性E(θ)から受信アンテナ202の指向性g(θ)を除去した指向性成分AF(θ)は、アレイファクターと呼ばれる。アレイファクターAF(θ)が、アレイアンテナを形成した事による指向性の効果を表す。受信アンテナ202n(n=1,2,・・・,N)で受信した信号は、g(θ)anexp(jφn)である。また、移相器206nの位相回転Φnを受けた信号g(θ)anexp(jφn)exp(jΦn)がn=1,2,・・・,Nに渡って加算器207で加算されて得られた信号が、式(5)の指向性E(θ)として得られる。
到来波2081、2082、・・・、208Nの入射角をθとした場合、到来波208nの位相φnは、-2π・n・d・sinθ/λで与えられる(n=1,2,・・・,N)。ここで、dは受信アンテナ202n(n=1,2,・・・,N)の間隔であり、λは到来波2081、2082、・・・、208Nの波長である。
上記の式(5)において、振幅anがnによらず一定とした場合、移相器206nの位相回転Φn(n=1,2,・・・,N)が到来波208nの位相φnに-1を掛けた値と等しくなるように設定すると、アレイファクターAF(θ)は角度θの方向において最大となる。このことは、即ち、移相器206nの位相回転Φnによるアレイアンテナの指向性の制御法を示している。
アレイアンテナ方式による物体検知装置の例は、その他に、特許文献1~3においても開示されている。具体的には、特許文献1及び2に開示された物体検知装置は、受信機に内蔵されたN個の受信アンテナそれぞれに接続された移相器により、N個の受信アンテナで形成される受信アレイアンテナの指向性を制御する。
そして、特許文献1及び2に開示された物体検知装置は、ビーム状に形成されたN個の受信アレイアンテナの指向性を変化させ、K個の検知対象物それぞれに対して、受信アレイアンテナの指向性ビームを向ける。これにより、各検知対象物で反射された電波強度が算出される。
また、特許文献3に開示された物体検知装置は、N個の受信アレイアンテナの周波数依存性を利用する事で、N個の受信アレイアンテナの指向性を制御している。また、特許文献3に開示された物体検知装置も、特許文献1及び2の例と同様に、K個の検知対象物それぞれに対して、N個の受信アレイアンテナの指向性ビームを向ける事で、各検知対象物で反射された電波強度を算出する。
また、実際の物体検知装置は、2次元画像を表示するため、図25で示すように、受信アンテナ202は、縦方向と横方向とにそれぞれN個ずつ設置されている。この場合、全体で必要なアンテナの数はN2個となる。図25は、従来からのアレイアンテナ方式を採用した場合の受信アレイアンテナの概略構成を示す図である。
また、2次元画像を表示するための方式としては、Mills Cross方式も知られている(例えば、非特許文献2参照)。図26は、Mills Cross方式を採用した物体検知装置を示す図である。図26に示すように、この物体検知装置は、縦方向に配列された1次元のアレイアンテナ201と、横方向に配列された1次元のアレイアンテナ201とを備えている。そして、この物体検知装置では、乗算器221は、縦方向にある受信アンテナと横方向にある受信アンテナとの組毎に、信号の積を算出する。よって、算出された積を用いることで、2次元画像を表示することが可能となる。
続いて、図27を用いて、物体検知装置における画像化の他の方式として、合成開口レーダー(SAR:Synthetic Aperture Radar)方式について説明する。図27は、従来からの合成開口レーダ方式を採用した物体検知装置を示す図である。
図27に示すように、合成開口レーダー方式において、物体検知装置は、送信機311と受信機301とを備えている。また、送信機311は、送信アンテナ312を備えている。受信機301は、受信アンテナ3021~302Nを備えている(Nは受信アンテナの数)。
送信機311は、送信アンテナ312から、RF信号(電波)313を検知対象物3041、3042、・・・、304K(Kは検知対象物の数)に向けて照射する。RF信号(電波)313は、検知対象物3041、3042、・・・、304Kにおいて反射され、反射波3031、3032、・・・、303Lがそれぞれ発生する。
この時、受信機3011は、最初の位置から、3012、・・・、301Nの位置に移動しながら、各位置において反射波3031、3032、・・・、303Kを受信する。図25において、3021、3022、・・・、302Nは、それぞれ、各位置での受信アンテナを示している。
また、これにより、1つの受信アンテナは、受信アンテナ3021、3022、・・・、302Nとして機能する。即ち、図27にいては、1つの受信アンテナが、図22で示したアレイアンテナ方式における受信アンテナ2021、2022、・・・、202Nと同じく、N本のアンテナによる受信アレイアンテナ(仮想アレイアンテナ)を形成する。
従って、図27で示した合成開口レーダー方式においても、図22で示したアレイアンテナ方式と同じく、受信機301は、受信した反射波3031、3032、・・・、303Kに基づいて、検知対象物3041、3042、・・・、304Kから反射されている電波強度を算出する。その後、受信機301は、算出した電波強度の分布を画像化する。これにより、検知対象物3041、3042、・・・、304Kそれぞれの像が得られることになる。
なお、合成開口レーダー方式による物体検知装置の例が、特許文献4~6において開示されている。
菊間信良、"アレーアンテナの基礎"、MWE2010 Digest,(2010)
B. R. Slattery,"Use of Mills cross receiving arrays in radar systems," PROC.IEE,Vol.113,No.11,NOVEMBER 1966, pp.1712-1722.
ところで、アレイアンテナ方式においては、対象物を精度良く検知しようとすると、必要となる受信アンテナの数とそれに付随する受信機の数とが非常に多くなってしまい、結果として、物体検知装置のコスト、サイズ、及び重量が大きくなるという問題がある。
上記の問題点について具体的に説明する。まず、アレイアンテナ方式の場合、受信アンテナ2011、2022、・・・、202Nの各アンテナの間隔は、受信機201において受信される反射波2031、2032、・・・、203Kの波長λの半分以下にする必要がある。例えば、反射波2031、2032、・・・、203Kがミリ波である場合は、波長λは数mm程度であるので、各アンテナの間隔は数mm以下となる。そして、この条件が満たされない場合は、生成した画像において、対象物2041、2042、・・・、204Kが存在しない位置に、虚像が発生するという問題が生じてしまう。
また、画像の分解能は受信アレイアンテナ(2011、2022、・・・、202N)の指向性ビーム幅△θで決まる。受信アレイアンテナ(2011、2022、・・・、202N)の指向性ビームの幅△θは、△θ~λ/Dにて与えられる。ここで、Dは受信アレイアンテナ(2011、2022、・・・、202N)の開口サイズであり、両端に存在する受信アンテナ2021と202Nと間の距離に相当する。つまり、衣服の下に隠されている物品又は鞄の中の物品等の画像化において実用的な分解能を得るには、受信アレイアンテナ(2011、2022、・・・、202N)の開口サイズDは数十cmから数m程度に設定されている必要がある。
上記の2つの条件、即ち、N個の受信アンテナのアンテナ間の間隔は波長λの半分以下(数mm以下)とする点と、両端に存在する受信アンテナ間の距離が少なくとも数十cm程度必要という点とから、一列あたりに必要なアンテナの数Nは数百個程度となる。
また、実際の物体検知装置では、2次元画像を表示するため、図26で示したように、受信アンテナ202は、縦方向と横方向とにそれぞれN個ずつ設置されている。この場合、全体で必要な受信アンテナの数はN2個となる。従って、アレイアンテナ方式を採用するためには、全体で必要な受信アンテナ及びそれに付随する受信機の数は数万個程度となる。
このように大量の受信アンテナと受信機とが必要となるため、上述したように、アレイアンテナ方式においては、コストは非常に高いものになる。また、一辺が数十cm~数mの四方の領域にアンテナが設置されるので、装置のサイズ及び重量は非常に大きなものとなる。
また、上述した図26に示したMills Cross方式の物体検知装置によれば、アレイアンテナ方式を採用する場合よりは、受信アンテナ及び受信機の数を減らすことは可能である。しかし、この場合であっても、必要な受信アンテナ及び受信機の数は、2N個であり、やはり数百個程度の受信アンテナが必要となる。従って、この場合であっても、コスト、装置サイズ及び重量の問題を解決することは困難である。
また、上述した図27に示した合成開口レーダー方式を採用した物体検知装置においては、受信機を機械的に動かす必要があるため、走査時間の短縮が難しいという問題がある。そして、この問題は、物体検知装置によって、物品又は人を検査する時に、単位時間当りに検査できる対象物の数が限られるという問題につながる。また、特許文献6に開示されている物体検知装置においては、受信機を動かすための機械的な機構を必要としているため、装置のサイズ及び重量が増大するという問題が発生している。
上記で議論したように、一般的な物体検知装置では、装置のコスト、サイズ、重量が非常に大きなものになる。このため、物体検知装置を実際に使用できる用途及び機会は、限定されたものになる。また、採用する方式によっては、対象物を検査する速度も限られたものになる。
本発明の目的の一例は、上記問題を解消し、電波を用いた物体の検知において、精度を向上させつつ、装置コスト、サイズ、及び重量の増大化を抑制し得る、物体検知装置及び物体検知方法を提供することにある。
上記目的を達成するため、本発明の一側面における物体検知装置は、電波によって物体を検知するための物体検知装置であって、
時間の経過と共に周波数が連続的に変化する電波を、送信信号として放射する、送信部と、
前記送信信号を取得し、前記物体で反射された前記電波を受信信号として受信し、更に、受信した前記受信信号に、取得した前記送信信号をミキシングして、ベースバンド信号を生成する、受信部と、
サンプリング時間毎の前記ベースバンド信号の測定値から、前記電波の到来方向を推定し、推定した前記電波の到来方向に基づいて、前記電波の強度分布を特定し、特定した前記強度分布に基づいて、前記物体を検知する、データ処理部と、
を備えている、ことを特徴とする。
時間の経過と共に周波数が連続的に変化する電波を、送信信号として放射する、送信部と、
前記送信信号を取得し、前記物体で反射された前記電波を受信信号として受信し、更に、受信した前記受信信号に、取得した前記送信信号をミキシングして、ベースバンド信号を生成する、受信部と、
サンプリング時間毎の前記ベースバンド信号の測定値から、前記電波の到来方向を推定し、推定した前記電波の到来方向に基づいて、前記電波の強度分布を特定し、特定した前記強度分布に基づいて、前記物体を検知する、データ処理部と、
を備えている、ことを特徴とする。
また、上記目的を達成するため、本発明の一側面における物体検知方法は、電波によって物体を検知するための方法であって、
(a)送信機によって、時間の経過と共に周波数が連続的に変化する電波を、送信信号として放射する、ステップと、
(b)受信機によって、前記送信信号を取得し、前記物体で反射された前記電波を受信信号として受信し、更に、受信した前記受信信号に、取得した前記送信信号を加算して、ベースバンド信号を生成する、ステップと、
(c)データ処理装置によって、サンプリング時間毎の前記ベースバンド信号の測定値から、前記電波の到来方向を推定し、推定した前記電波の到来方向に基づいて、前記電波の強度分布を特定し、特定した前記強度分布に基づいて、前記物体を検知する、ステップと、
を有する、ことを特徴とする。
(a)送信機によって、時間の経過と共に周波数が連続的に変化する電波を、送信信号として放射する、ステップと、
(b)受信機によって、前記送信信号を取得し、前記物体で反射された前記電波を受信信号として受信し、更に、受信した前記受信信号に、取得した前記送信信号を加算して、ベースバンド信号を生成する、ステップと、
(c)データ処理装置によって、サンプリング時間毎の前記ベースバンド信号の測定値から、前記電波の到来方向を推定し、推定した前記電波の到来方向に基づいて、前記電波の強度分布を特定し、特定した前記強度分布に基づいて、前記物体を検知する、ステップと、
を有する、ことを特徴とする。
以上のように本発明によれば、電波を用いた物体の検知において、精度を向上させつつ、装置コスト、サイズ、及び重量の増大化を抑制することができる。
(実施の形態1)
以下、本発明の実施の形態1における物体検知装置及び物体検知方法について、図1~図10を参照しながら説明する。
以下、本発明の実施の形態1における物体検知装置及び物体検知方法について、図1~図10を参照しながら説明する。
[装置構成]
最初に、図1を用いて、本実施の形態1における物体検知装置の概略構成について説明する。図1は、本発明の実施の形態1における物体検知装置の構成を概略的に示す構成図である。
最初に、図1を用いて、本実施の形態1における物体検知装置の概略構成について説明する。図1は、本発明の実施の形態1における物体検知装置の構成を概略的に示す構成図である。
図1に示す本実施の形態1における物体検知装置1000は、電波によって物体1001を検知するための装置である。図1に示すように、物体検知装置1000は、送信部1091と、受信部1092と、データ処理部1093とを備えている。
送信部1091は、時間の経過と共に周波数が連続的に変化する電波を、送信信号として放射する。受信部1092は、送信信号を取得し、検知対象となる物体(以下「対象物」と表記する。)1001からの電波を受信信号として受信する。更に、受信部1092は、受信した受信信号に、取得した送信信号を掛算(ミキシング)して、ベースバンド信号を生成する。
また、図1に示すように、本実施の形態1では、送信部1091は送信アンテナ1003を備え、受信部1092は受信アンテナ1004を備えている。図1の例では、単一の受信部1092のみが図示されているが、本実施の形態1において、受信部1092及び受信アンテナ1004の数は複数であっても良い。但し、本実施の形態1においては、受信部1092及び受信アンテナ1004の数は、従来に比べて極めて少なくなる。
データ処理部1093は、サンプリング時間毎のベースバンド信号の測定値から、電波の到来方向を推定する。そして、データ処理部1093は、推定した電波の到来方向に基づいて、電波の強度分布を特定し、特定した強度分布に基づいて、対象物1001を検知する。
ここで、物体検知装置1000の動作原理についてまず図2及び図3を用いて説明する。図2は、本発明の実施の形態1における物体検知装置の動作原理を説明する図である。図3は、本発明の実施の形態1で送信される電波の周波数の変化を示す図である。
図2で示した例では、x軸上に送信アンテナ1003と受信アンテナ1004とが配置され、K個の対象物10011、・・・1001Kが、(x1,z0),・・・, (xK,z0)の位置に、それぞれ配置されている。送信アンテナ1003からは、図3で示すように、キャリア周波数fが線形的に変化するRF信号1010が送信されるものとする。キャリア周波数fは、1チャープ周期(Tchirp)の間に、最小周波数fminから最大周波数fmaxまで変化するものとする。キャリア周波数fの帯域幅をBW(=fmax-fmin)、キャリア周波数の時間傾きをα=BW/Tchirpと定義する。
また、図2で示した例において、送信アンテナ1003から対象物10011、・・・1001Kに向けてRF信号10101、・・・1010Kが、それぞれ照射されているとする。更に、対象物10011、・・・1001Kからの反射波10071、・・・1007Kが、受信アンテナ1004で受信されるものとする。
受信アンテナ1004で受信された反射波10071、・・・1007Kの合成波は、図1で示した受信部1092において、送信部1091から取得した送信信号と掛算(ミキシング)される。これにより、ベースバンド信号が生成される。ベースバンド信号I(t)は、以下の式(6)で与えられる。
上記式(6)において、t’は1チャープ周期内の時刻であり、図3におけるt0からtMに相当する。hをチャープ番号として、t’ = t-h・Tchirpと表記されるとおり、1チャープ周期(Tchirp)を過ぎる毎にt’はt0から取り直す必要がある。σ(xk)は、対象物k (k=1,2,・・・,K)の反射率である。L(xk)は、対象物kを経由した送信アンテナから受信アンテナまでの伝搬距離である。cは光速である。
上記のIF信号I(t)は、同相成分(In-phase信号)である。同相成分I(t)にヒルベルト変換を掛けることで、直交成分(Quadrature信号)Q(t)が生成される。また、直交変調器を用いて同相成分I(t)と直交成分Q(t)とが生成されても良い。直交成分Q(t)は、以下の式(7)で与えられる。
そして、同相成分I(t)と直交成分Q(t)とから、以下の式(8)で表される複素ベースバンド信号r(t)が、生成される。
複素ベースバンド信号r(t)は、実測データから算出できる量である。ここでの目的は、実測データから得られる複素信号r(t)から、反射率σの位置xに対する依存性、特にσ(x)=0となる位置xを求める事である。σ(x)=0となる位置xが分かれば、対象物の位置や形状を決定できる。
式(8)において、σ’(x)という量が定義されている。σ(x) = 0 とσ’(x) = 0との間には同値であるという関係があるので、σ’(x) = 0となる位置xを求める事が、ここでの目的と言い直す事もできる。
上記式(8)は、以下の式(9)のように表記できる。
上記式(9)において、t1,t2,・・・,tNは、1チャープ周期内のサンプリング時間である。ここで、Nは1チャープ周期あたりのサンプリング点数となる。また、Δtはサンプリング周期であり、Δt = tn+1-tnで与えられる。式(8)から式(9)への展開にあたり、受信信号rに対し、ノイズ成分(乱数)を要素とするベクトルn(t)が付加されている。
背景技術で記した従来のアンテナアレイの動作を示す式(1)と、本実施の形態における動作を示す式(9)とを比較すると、図4で示すパラメータの対応関係(読み替え)を付ける事で、両者は同型である事が分かる。この事を利用して、従来のアンテナアレイで用いた到来方向推定アルゴリズムと同型の方式を、そのまま本実施の形態においても適用して電波の到来方向推定を実施できる。図4は、従来からのアレイアンテナ方式で用いられるパラメータと本発明の実施の形態における時間仮想アレイ方式で用いられるパラメータとの対応関係を示す図である。
すなわち、本実施の形態においては、測定で得た式(9)で定義される所の受信信号(複素ベースバンド信号)rから、相関行列R=E[r・rH]が計算され、次に、計算された相関行列Rから、以下の式(10)~(12)に示すいずれかの評価関数が計算される。
上記式(10)~(12)において、方向ベクトルa(x)は、式(9)で定義される所のものを使用する。また、上記のMUSIC法におけるEN=[eK+1,・・・,eN]は、相関行列Rの固有ベクトルの内、固有値がノイズn(t)の電力となるN-(K+1)個のベクトルで構成した行列である。
式(10)~(12)で示した評価関数は、対象物の存在位置x1,x2,・・・,xKにおいてピークを持つ。従って、評価関数を計算してそのピークを見れば、対象物の位置(存在領域)を求める事ができる。式(10)~(12)の評価関数で得られた対象物の位置分布から、対象物の位置及び形状を画像として表示する事ができる。以上の説明が本実施の形態におけるとなる。
続いて、本実施の形態における原理をより直観的に理解できるように説明を加える。ここでは特に、図4で示したパラメータの対応関係の内、従来のアンテナアレイにおけるアンテナ間隔dと、本実施の形態におけるサンプリング時間Δtとの間の対応関係に着目する。この対応関係に着目すると、従来のアンテナアレイがアンテナ間隔dを持って配置されたN本のアンテナで受信したデータを用いて電波の到来方向推定を行っているのに対し、本実施の形態ではサンプリング時間Δt毎に得られるN個の受信データで電波の到来方向の推定が行なわれていると解釈できる。言い換えると、本実施の形態では、図5に示すように、サンプリング時間毎に得られるデータを仮想的なアンテナと見なし、時間軸上に配置されたN本の仮想アンテナで仮想的なアンテナアレイ(時間仮想アレイ)を構築して到来方向推定を実施していると解釈できる。
図5は、本発明の実施の形態1における物体検知装置の動作原理を示す図である。図5に示す例では、サンプリング時間毎に、仮想的な送信アンテナ1003(t1)、1003(t2)、・・・、1003(tN)と、仮想的な受信アンテナ1004(t1)、1004(t2)、・・・、1004(tN)とによって、測定が行なわれる。
図5において、受信機1092は、送信部1091から取得した送信信号1003(t1)、1003(t2)、・・・、1003(tN)と受信信号1004(t1)、1004(t2)、・・・、1004(tN)を掛算(ミキシング)して受信信号(複素ベースバンド信号)rを生成する。生成した受信信号rは、信号処理部1095に出力される。そして、受信信号rから相関行列Rを計算する過程、更には、式(10)~(12)の評価関数を計算する過程は、信号処理部1095で実施される。なお、信号処理部1095は、図1におけるデータ処理部1093に含まれるものである。
図5で示した本実施の形態における物体検知装置の構成において、信号処理部1095をビームフォーマ法に特化した構成とすると、図6に示した構成となる。図6は、図5に示した物体検知装置にビームフォーマ法を適用した場合の動作原理を示す図である。図6において、信号処理部1095は、移相器1031と加算器1032で構成される。(以下ではサンプリング点数をNからMに差し替える。)
仮想的な受信アンテナ1004(t1)、1004(t2)、・・・、1004(tM)で受信した反射波1007(又はその複素振幅)は、移相器1031(t1)、1031(t2)、・・・、1031(tM)において位相回転Φ1、Φ2、・・・、ΦMを受けた後、加算器1032で加算される。
本実施の形態1においては、移相器1031(t1)、1031(t2)、・・・、1031(tM)による位相回転と、加算器1032による加算とは、データ処理部1093における処理、具体的には、プロセッサを用いたソフトウェアによる処理によって実行できる。
本実施の形態1における物体検知装置1000の原理は、既に述べたとおりサンプリング時間t1、t2、・・・、tMそれぞれにおける測定データで仮想的なアレイアンテナを構築し、その仮想的なアレイアンテナで到来波の方向を推定する事である。従って、図25で示した一般的なアレイアンテナと同じく、図6で示した仮想アレイにおいてもアレイファクターAF(xd)を計算できる。
ここで、x軸とz軸とによる位置座標を設定し、送信部1091の位置を(0,0)、受信部1092の位置を(xr,0)、対象物1001の位置を(xd,z)とする。仮想的な受信アンテナ21(tm)(m=1,2,・・・,M)で受信した反射波1007(tm)の振幅と位相とをそれぞれamおよびφmとした場合、本発明の仮想アレイにおけるアレイファクターAF(xd)は以下の式(13)のように計算される。
また、反射波102(tm)の位相φm(m=1,2,・・・,M)は以下の式(14)で与えられる。
ここで、式(14)において、α・Δtは、サンプリング毎のキャリア周波数fの差分(周波数間隔)である。Lt(xd)は送信部1091と対象物1001との距離である、Lr(xd)は受信部20と対象物1001との距離である。cは光速である。また、式(3)において、振幅amがmによらず一定であるとした場合、移相器1031(tm)による位相回転Φm(m=1,2,・・・,M)を反射波1007(tm)の位相φmと等しくなるように設定すると、アレイファクターAF(xd)は対象物1001(位置xd)の方向において最大となる。この事が、本実施の形態1における、移相器22(tm)の位相回転Φm(m=1,2,・・・,M)による仮想アレイの指向性の制御法を示している。
図7は、本発明の実施の形態1における物体検知装置のアンテナ利得の指向性の一例を示す特性図である。具体的には、図7は、上記の式(2)及び(3)を用いた仮想アレイのアレイファクターAF(xd)の計算結果を示している。
図7の例では、物体200の位置(xd,z)に対し、仮想アレイのビーム中心がxd=80cm、100cm、120cmの位置となるように、移相器22(tm)の位相回転Φm(m=1,2,・・・,M)が設定されている。そして、図7の例では、このような場合の、仮想アレイのアレイファクター(すなわちビームパターン)が示されている。
また、図7の例の計算において、周波数間隔α△t=250MHz、サンプリング数M=21、対象物1001のz軸座標位置z=100cm、送信部1091と受信部1092との距離xr=100cmに設定されている。
このように、図7から分かるように、本実施の形態1における仮想アレイにおいても、移相器1031(tm)の位相回転Φm(m=1,2,・・・,M)によって、仮想アレイの指向性(ビームパターン)の制御が可能である。また、上記の式(2)及び式(3)で与えられるアレイファクターAF(xd)から、ビームパターンのビーム幅を計算する事ができる。ビーム幅は、到来方向の推定とイメージング(画像)の分解能とを決定する要素である。また、本実施の形態1において、ビーム幅△xは、以下の式(15)で与えられる。
式(15)において、BWは上述したようにRFキャリア周波数の帯域幅である。周波数間隔α△tとサンプリング数Mとを用いて、BW=α△t×Mと表すことができる。また、式(4)において、h(xr,xd,z)は、位置変数(xr,xd,z)の関数である。なお、xr=xdの場合、h(xr,xd,z)は[1+(z/xr)2]1/2で与えられる。
式(15)が示すように、本実施の形態1における仮想アレイでは、帯域幅BWを広げるほど、ビーム幅△xが縮まり、より高分解能の性能が得られる。但し、一般的なアレイアンテナと同じく、本実施の形態1における仮想アレイでも、グレーティングローブによる虚像が発生する可能性がある。以下に、図8を用いて、虚像の発生について説明する。図8は、本発明の実施の形態に1おける仮想アレイでの虚像の発生を説明する図である。
また、式(16)によって、位相量φ(xa)を定義する。
また、式(16)によって、位相量φ(xa)を定義する。
式(16)における位相量φ(xa)は、図8において、虚像1033(位置xa)を経由して送信部1091から受信部1092に至るまでの電波の位相シフトと、対象物1001(位置xd)を経由して送信部1091から受信部1092に至るまでの電波の位相シフトとの差分に対応する。そして、位置xaにおいて、位相φ(xa)が2πの整数倍となる場合、位置xaと対象物の位置xdとで同じアレイファクターが得られる。即ち、位置xaに実際に対象物が存在しない場合でも、位置xaにおいて対象物1001の像(即ち、虚像1033)が発生する事となる。そのため、|φ(x)|<πを満たす領域、すなわち以下の条件式(17)を満たす位置xの範囲が、虚像の発生しない領域(可視領域)として用いる事ができる。
式(17)から、周波数間隔α△tを小さくするほど、即ち、サンプリング間隔を短くするほど、可視領域が広がる事が分かる。可視領域の大きさ(長さ)は概ね周波数間隔α・△tに反比例する。
このように、仮想アレイを用いて反射波の到来方向を推定し、その結果からイメージング処理(画像生成)を行う場合、一方向あたりの画素数は可視領域と分解能の比で与えられる。式(15)と式(17)とが示す結果から、一方向あたりの画素数=可視領域/分解能∝BW/α△t=Mという関係が得られる(BWは帯域幅、α△tは周波数間隔、Mはサンプリング数)。即ち、本実施の形態1では、必要な画素数に応じて、サンプリング数Mを設定すれば良い事となる。
本実施の形態1では、データ処理部1093によって、受信部1092が出力したベースバンド信号のサンプリング時間毎の測定値それぞれに対して、位相が制御される。そして、この位相の制御により、受信部1092における実効的なアンテナ利得の指向性が制御され、更に、アンテナ利得の指向性の制御により、受信部1092に到来する電波の強度分布が測定されるので、対象物1001の位置及び形状の検知が可能となる。このため、従来のように、大量の受信アンテナ及び受信機を用意する必要がない。本実施の形態1によれば、電波を用いた物体の検知において、精度を向上させつつ、装置コスト、サイズ、及び重量の増大化を抑制できる。
続いて、図9を用いて本実施の形態1における物体検知装置の具体的構成について説明する。図9は、本発明の実施の形態1における物体検知装置の具体的構成の一例を示すブロック図である。
図9に示すように、本実施の形態1では、物体検知装置1000は、送信部1091、受信部1092、及びデータ処理部1093に加えて、出力部1094も備えている。また、実際には、本実施の形態1では、送信部1091は送信機で構成され、受信部1092は受信機で構成されている。また、データ処理部1093は、データ処理装置、即ち、計算機(コンピュータ)で構成されている。以下、具体的に説明する。
図9に示すように、送信部1091は、送信アンテナ1003に加え、少なくとも電力増幅器1071と、カプラ1075と、周波数可変機能を持つ発振器1103と、送信制御部1104とを備えている。
送信部1091において、発振器1103は、送信RF信号を出力する。発振器1103から出力された送信RF信号は、電力増幅器1071において増幅された後、送信アンテナ1003から送信RF信号1010として送出される。
送信部1091において、送信制御部1104は、発振器1103が出力するRF信号の周波数を制御する。本実施の形態1では、発振器1103が出力するRF信号の周波数(=送信RF信号1010のキャリア周波数)は、時間の経過と共に連続的に変化するように制御される。特に、図3で示すようにRF信号の周波数を制御する事が望ましい実施の形態である。
また、発振器1103が出力したRF信号は、カプラ1075を経由して、受信部1092内のミキサ1042へと出力される。後述するように、カプラ1075を経由してミキサ1042に出力されたRF信号は、受信部1092のLO信号として使用される。
また、図9に示すように、受信部1092は、受信アンテナ1004に加え、低雑音増幅器1041と、ミキサ1042と、フィルタ1043と、アナログ-デジタル変換器1044と、受信制御部1102とを備えている。
上記図1~図6を用いて説明したとおり、受信部1092は、それに備えられた受信アンテナ1004によって、対象物1001から反射された電波(RF信号)1007を受信する。受信アンテナ1004で受信されたRF信号1007は、低雑音増幅器1041で増幅された後、ミキサ1042に入力される。
ミキサ1042は、低雑音増幅器1041で増幅された受信RF信号と、カプラ1075を経由して送信部1091から出力されてきたRF信号(受信LO信号)とをミキシングして、ベースバンド信号となる中間周波数信号(IF信号)を生成し、これをフィルタ1043に向けて出力する。フィルタ1043は、ベースバンド信号からノイズを除去し、ノイズが除去されたベースバンド信号をアナログ-デジタル変換器1044に入力する。
アナログ-デジタル変換器1044は、アナログ信号であるベースバンド信号をデジタルベースバンド信号に変換し、得られたデジタルベースバンド信号を受信制御部1102に入力する。上記で得たデジタルベースバンド信号は、式(6)に記載の同相成分(In-phase信号)I(t)に相当する。
受信制御部1102は、同相成分I(t)にヒルベルト変換を掛けて、直交成分(Quadrature信号)Q(t)を生成する。更に、受信制御部1102は、同相成分I(t)と直交成分Q(t)とから式(8)に従って、複素ベースバンド信号r(t)を生成する。生成された複素ベースバンド信号r(t)は、データ処理部1093へと受け渡される。なお、上述したとおり、ミキサ1042の代わりに、直交変調器が用いられて、直交成分Q(t)が生成されても良い。
データ処理部1093は、受け渡された複素ベースバンド信号r(t)に対して、図2~図8を用いて説明した処理、即ち、受信した電波1007の到来方向の推定処理を実行する。更に、データ処理部1106は、対象物1001のイメージング処理(画像生成)も実行する。その後、データ処理部1106は、処理の結果、即ち、推定した到来方向と生成した画像とを、出力部1094に出力する。出力部1094は、例えば、表示装置であり、画面上に、処理の結果を表示する。
図9で示した例では、送信部1091と受信部1092はそれぞれ一つずつ示されているが、本実施の形態1は、この例に限定されない。本実施の形態1では、物体検知装置1000は、送信部1091と受信部1092とを、それぞれ複数備えていても良い。また、データ処理部1093及び出力部1094は、送信部1091又は受信部1092に内蔵されていても良い。
[装置動作]
次に、本発明の実施の形態1における物体検知装置1000の動作について図7を用いて説明する。図7は、本発明の実施の形態1における物体検知装置100の動作を示すフロー図である。以下の説明においては、適宜図1~図8を参酌する。また、本実施の形態1では、物体検知装置を動作させることによって、物体検知方法が実施される。よって、本実施の形態1における物体検知方法の説明は、以下の物体検知装置1000の動作説明に代える。
次に、本発明の実施の形態1における物体検知装置1000の動作について図7を用いて説明する。図7は、本発明の実施の形態1における物体検知装置100の動作を示すフロー図である。以下の説明においては、適宜図1~図8を参酌する。また、本実施の形態1では、物体検知装置を動作させることによって、物体検知方法が実施される。よって、本実施の形態1における物体検知方法の説明は、以下の物体検知装置1000の動作説明に代える。
図10に示すように、最初に、送信部1091において、送信制御部1104は、現在のサンプリング時間tmを特定し、送信アンテナ1003が送出するRF信号の周波数(fmin+αtm)を算出する(ステップA1)。
次に、送信制御部1104は、周波数(fmin+αtm)のRF信号が送信アンテナ1003から送出されるように、発振器1103の制御信号を生成し、これを出力することによって、送信アンテナ1003から、周波数が(fmin+αtm)のRF信号を送出させる(ステップA2)。
具体的には、送信制御部1104は、発振器1103の出力周波数が(fmin+αtm)となるように制御信号を発振器1103に向けて送出し、発振器1103はキャリア周波数が(fmin+αtm)のRF信号を出力する。これにより、同RF信号は、電力増幅器1071で増幅され、送信アンテナ1003から送出される。
また、発振器1103が出力したRF信号は、カプラ1075を経由して、受信部1092内のミキサ1042に対しても送出される。
次に、受信部1092において、受信アンテナ1004が、対象物1001から反射された電波(RF信号)1007を受信する(ステップA3)。
次に、受信制御部1102は、受信されたRF信号から得られたベースバンド信号の同相成分I(t)から、複素ベースバンド信号r(t)を算出する(ステップA4)。
具体的には、ステップA4では、まず、受信アンテナ1004で受信されたRF信号1007は、低雑音増幅器1041で増幅された後、ミキサ1042に入力される。ミキサ1042は、低雑音増幅器1041で増幅された受信RF信号に、カプラ1075経由で送信部1091から出力されてきたRF信号をLO信号としてミキシングして、ベースバンド信号(同相成分I(t))を生成する。ベースバンド信号(同相成分I(t))は、フィルタ1043を経由して、アナログ-デジタル変換器1044に入力され、そこでデジタル信号に変換される。受信制御部1102は、このデジタル変換されたベースバンド信号(同相成分I(t))から、複素ベースバンド信号r(t)を算出する。
次に、データ処理部1093は、複素ベースバンド信号r(t)を用いて、受信した電波1007の到来方向を推定し、更に、推定結果を用いて、対象物1001のイメージング処理を実行する(ステップA5)。
また、本実施の形態1では、ステップA1~A5は繰り返し実行され、繰り返し行なわれた処理の結果は、出力部1094によって、画面上に、表示される。
[実施の形態1による効果]
以上のように、本実施の形態1によれば、従来のように、大量の受信アンテナ及び受信機を用意することなく、物体を精度良く検知できる。また、受信アンテナの数を増やす必要がないので、装置コスト、サイズ、及び重量の増大化が抑制される。
以上のように、本実施の形態1によれば、従来のように、大量の受信アンテナ及び受信機を用意することなく、物体を精度良く検知できる。また、受信アンテナの数を増やす必要がないので、装置コスト、サイズ、及び重量の増大化が抑制される。
加えて、本実施の形態1では、電波の送信及び受信の方式として、FM-CW方式が採用されている。このため、受信部1092に発振器を設ける必要がなく、この点でも装置コストが削減できる。更に、受信部1092で発振器が不要であるため、送信部1091内の発振器1103と受信部1092内の発振器で同期を取る必要が無く、結果として送信部1091と受信1092の間の同期エラーとそれに起因する検知精度の劣化も発生しない。
なお、実施の形態1における物体検知装置1000は、後述する実施の形態2及び実施の形態3において利用される。実施の形態1で行なわれる処理は、実施の形態2における対象物1001の位置(特に1次元の方向)を推定する処理、実施の形態3における対象物1001の配置状況及び形状を2次元画像で表示する処理に用いられる。これらの処理もまたデータ処理部1093において実施される。
(実施の形態2)
続いて、本発明の実施の形態2における物体検知装置及び物体検知方法について、図11~図13を参照しながら説明する。
続いて、本発明の実施の形態2における物体検知装置及び物体検知方法について、図11~図13を参照しながら説明する。
本実施の形態2は、実施の形態1で示した物体検知装置1000を用いて、対象物の位置、特に一次元の方向を推定する例を示している。従って、本実施の形態2においても、物体検知装置は、図1及び図9に示した、送信部1091、受信部1092、及びデータ処理部1093を備えている。但し、本実施の形態2は、受信部1092の個数の点で、実施の形態1と異なっている。以下、具体的に説明する。
図11は、本発明の実施の形態2における物体検知装置の構成及び動作原理を示す図である。まず、本実施の形態2においては、物体検知装置は、1つの送信部に対して、N個の受信部を備えている。従って、図11に示すように、本実施の形態2では、物体検知装置は、一本の送信アンテナ1003と、N本の受信アンテナ10041、・・・、1004n、・・・、1004Nとを備えている。なお、以下において、特定の受信アンテナを示さない場合では、「受信アンテナ1004」と表記することとする。
また、図11に示すように、本実施の形態2では、各受信アンテナは、送信アンテナを基準にした一方向に沿って設置されている。具体的には、送信アンテナ1003と各受信アンテナ1004とはx軸上(z=0)に設置されている。送信アンテナ1003の位置は(x,z)座標で(d0,0)とする。また、N本の受信アンテナ1004の位置をそれぞれ(dx1,0),(dx2,0),・・・,(dxN,0)とする。
なお、物体検知装置は、受信アンテナの数Nが最小の1であっても動作可能である。但し、ここでは、理論に一般性を持たせるため、N本の受信アンテナの場合を扱う。また、対象物1001はz=z0の軸上において、D個の位置(x1,z0),(x2,z0),・・・,(xD,z0)に設置されるものとする。また、説明を簡単にするため、送信アンテナ1003、受信アンテナ1004、及び対象物1001の位置は、上記の位置に固定されているものとする。
そして、このような構成において、データ処理部は、各受信アンテナ1004で受信されたベースバンド信号の測定値から、電波の到来方向を推定する。また、データ処理部は、推定した電波の到来方向に基づいて、電波の強度分布を特定し、特定した強度分布に基づいて、対象物1001の一方向における位置を検知する。
また、データ処理部は、サンプリング時間毎のベースバンド信号の測定値から時間仮想アレイを構築し、前記時間仮想アレイの相関行列を算出する。より詳細には、データ処理部は、サンプリング時間が異なるベースバンド信号の測定値から時間仮想アレイのサブアレイを構築し、前記サブアレイ毎の相関行列を算出し、前記サブアレイ毎の相関行列の平均値を算出する。そして、データ処理部は、前記相関行列の平均値に基づいて、対象物1001の位置を反映する評価関数を求め、求めた評価関数から対象物1001の画像を生成する。以下に、本実施の形態2における物体検知装置で行なわれる処理について具体的に説明する。
まず、本実施の形態2においても、本実施の形態1と同じく、送信アンテナ1003からはFM-CW信号が送出される。
受信アンテナ1004は、対象物1001からの反射波1007を受信する。
ここで、d番目(d=1,2,・・・,D)の対象物1001dで反射され、n番目の受信アンテナ1004nで受信した反射波1007から得られる複素ベースバンド信号をsxn(xd,tm)とする。添字の「xn」は、x軸方向に配置されたn番目の受信アンテナ1004nで受信された信号であることを意味している。また、ここではサンプリング時間tm(m=1,2,・・・,M)における複素ベースバンド信号sxn(xd,tm)が取得すべきデータとなる。
各受信アンテナ1004nで実際に受信される信号は、それぞれ全ての対象物1001d(d=1,2,・・・,D)からの反射波1007の合成であり、個別対象からの反射波1007の複素振幅sxn(xd,tm)は未知数である。受信アンテナ1004nで実際に測定される信号の複素振幅をsxn’(tm)とすると、sxn’(tm)とsxn(xd,tm)との間に以下の関係がある。
なお、上記式(18)におけるsxn’(tm)は、実施の形態1で説明した式(8)の複素ベースバンド信号r(t)に相当するものである。
次に、各対象物1001d(d=1,2,・・・,D)から反射され、n番目の受信アンテナ1004nで受信された反射波1007の複素振幅sxn(xd,tm)について、詳細解析する。送信アンテナ1003と対象物1001dまでの距離L0(xd)と、n番目の受信アンテナ1004nと対象物1001dまでの距離Lxn(xd)は、以下の式(19)と(20)で与えられる。
送信アンテナ1003から送出されるRF信号1010の複素振幅s0と、n番目の受信アンテナ1004nで受信された反射波1007から得られる複素振幅sxn(xd,tm)との間には、以下の関係がある。
式(21)において、σ(xd)は対象物1001dの反射率を表す未知数である。式(21)右辺内の指数項は、対象物1001d経由で送信アンテナ1003から受信アンテナ1004nに至るまでの経路で生じる電波の位相シフトを表している。式(21)を式(18)に代入する事で、以下の式(22)が得られる。
続いて、データ処理部で行なわれる処理(解析)について説明するが、その前に、いくつかの信号を以下に定義する。式(11)左辺の信号sxn’(tm)(n=1,2,・・・,N,m=1,2,・・・,M)を用いて、測定信号ベクトルsxを以下の式(23)によって定義する。
添字[]Tはベクトル又は行列の転置を表す。次に、式(11)右辺内に含まれる指数項を用いて、方向行列Aを、以下の式(24)によって定義する。
式(24)において、行列AのサイズはMN×D、行列AnのサイズはM×D、ベクトルan(xd)のサイズはM×1となる。なお、本明細書では行列のサイズを縦×横の要素数で表記する。また、式(11)右辺内の変数s0とσ(xd)を用いて、所望信号ベクトルsを以下の式(25)によって定義する。
また、本実施の形態2では、受信アンテナ1004による測定で所望信号ベクトルsのxd依存性(即ち、σ(xd))を反映した評価関数を決定する事が目的となる。所望信号ベクトルsのxd依存性から、対象物1001の分布及び形状が検知される。上記の式(22)の関係は、測定信号ベクトルsx、方向行列A、所望信号ベクトルsを用いて、以下の式(26)のように表現できる。
なお、式(22)から式(26)への展開にあたり、式(26)の右辺にノイズ(乱数)を要素とするMN×1次のベクトルn(t)が新たに付加されている。このノイズ(乱数)n(t)の不可は、データ処理部において人為的に行なわれる。また、一つのサンプリング時間tmに対し、n(t)を定義する時間tの点数(スナップショット数)は1よりも大きい。
後述するようにMUSIC法の適用条件として行列Aはフルランクである事が要求される。ノイズベクトルn(t)を付加する事は、行列A内の列ベクトルないし行ベクトルの従属性を実効的に破壊し、行列Aをフルランクに近付ける効果がある。
本実施の形態2においては、式(23)で定義された測定信号ベクトルsx(t)が、受信アンテナ1004で受信される。そして、データ処理部は、受信された測定信号ベクトルsxを用いて、以下の式(27)に示す相関行列Rxを計算する。
式(27)内におけるE[]は、ノイズ(乱数)ベクトルn(t)を定義する時間tの点数(スナップショット数)に渡る平均を表す。
式(27)に示す相関行列Rxの定義に、上記の式(26)を代入する事で、以下の式(28)によって相関行列Rxと方向行列Aとの関係が導かれる。
式(28)において、PNはノイズ電力、IはMN×MN次の単位行列である。添字Hは複素共役転置を表す。相関行列Rx、行列A、行列SのサイズはそれぞれMN×MN次、MN×D次、D×D次となる。
ところで、非特許文献1に記載されているように、式(26)と式(28)とが成立する系に対しMUSIC法を適用する事で、所望信号ベクトルsの強度のx依存性(即ち、σ(x))を反映した評価関数PMU(x)を計算できる事が知られている。
但し、MUSIC法の適用条件として、式(28)内の行列Aと行列Sとがフルランクである事が要求される。フルランクとは、行列の階数が行列のサイズ(行数又は列数のうちの少ない方)と一致する事であり、行列内の全ての行ベクトル及び列ベクトルが全て線形独立である事と定義される。
方向行列Aは、各列ベクトルが異なる位置xdの関数であるので、各列ベクトルは独立でありフルランクとなる。行列Sの要素を見ると、σ(xi)=σ(xj)(i≠j)の場合、行列Sの第i行と第j行の行ベクトルが同じ値となり線形従属となるため、階数が一つ下がりフルランクでは無くなる。式(17)は連立方程式と見なせるが、行列Sの階数が減る事は、独立な方程式の数が減る事と等価であり、所望の未知数σ(xd)(d=1,2,・・・,D)の情報を得る事が困難になる。
以下では、サブアレイの概念を用いて行列Sをフルランクに戻す手法を示す。本実施の形態1で示したとおり、本実施の形態2においても、一つの周波数を一つのアンテナに見立てて仮想アレイが構築される。
本実施の形態2においては、図12に示すように、サンプリング時間を変えて測定した全てのデータを全体アレイ、サンプリング時間毎のデータをグループに分けてまとめたものをサブアレイと見なす。図12は、本発明の実施の形態2における物体検知装置で用いられるサブアレイの概念を説明する図である。
また、図12に示すように、全体アレイはM0個の周波数の測定データで構成され、サブアレイはM個(M0>M)の周波数の測定データで構成されている。サブアレイの数をQとすると、Q=M0-M+1の関係がある。サブアレイq(q=1,2,・・・,Q)の測定信号ベクトルsxqは、以下の式(29)によって定義される。
この時、式(29)のサブアレイqの測定信号ベクトルsxqには、式(24)の方向行列Aと式(14)の所望信号ベクトルsとの間において、以下の式(30)で与えられる関係がある。
ここで、サンプリング時間t1,t2,・・・,tMは等間隔であり、その間隔(サンプリング周期)をΔtとしている。すなわち、tm=m・Δt,(m=1,2,・・・,M)とする。サブアレイqの相関行列Rxqは、以下の式(31)のように計算される。
式(31)において、相関行列Rxq、行列A’、行列S’のサイズは、それぞれNM×NM次、NM×ND次、ND×ND次となる。次に全てのサブアレイq(q=1,2,・・・,Q)の相関行列の平均Rx’を計算する。全サブアレイ平均の相関行列Rx’と方向行列Aの関係は以下の式(32)のように計算される。
式(32)内の相関行列Rx’とは、式(17)の相関行列と同じくA’S”A’Hの形を持つ。そこで、行列A’とS”がフルランクであれば、相関行列Rx’にMUSIC法を適用して所望信号ベクトルsの強度のx依存性(すなわちσ(x))を反映した評価関数PMU(x)を計算できる。
方向行列A1,A2,・・・,ANは、それぞれ独立かつフルランクであるので、式(31)で与えられる行列A’もまたフルランクである。
次に行列S”について考察する。式(17)において、全ての対象物の反射率が同じ状況、すなわちσを定数としてσ=σ(x1)=σ(x2)=・・・=σ(xD)となっている状況を考える。この時、行列Sの階数は1となり、MUSIC法を適用する上では最も厳しい状況となる。このような状況においても、条件を満たせば式(21)の行列S’’がフルランクになる事を示す。σ=σ(x1)=σ(x2)=・・・=σ(xD)の場合に、式(32)の行列S’を計算した結果は、以下の式(33)となる。
行列Ciにおいて、biu=biv(u≠v)であれば、行列Cの第u行と第v行の行ベクトルが同じ値となり線形従属となるため、階数が一つ下がりフルランクでは無くなる。一方、式(30)で見られるように、bidは距離L0(xd)とLx(xd)の関数であり、位置xdが異なればこれらの距離は異なる値を取るので、biu=biv(u≠v)が満たされる事はなく、Ciはフルランクとなる。
Ciの行列サイズはD×Qであるので、Ciの階数はDとQの小さい方となる。したがってQ≧DであればCiのランクはDとなり、S”ijの階数もDとなりフルランクの条件が満たされる。また、各S’’ijは独立であるので、S”はフルランクとなる。
式(28)の行列Sは、位置xdが異なっても反射率σ(xd)は同じ値を取り得るという条件から、フルランクにならない場合があった。一方で行列S”は、位置xdが変化すれば距離L0(xd)とLx(xd)も必ず変化するという性質から、フルランクになる事が保証されている。
Q<Dの状況においてS”の階数はQになり、サブアレイの数Qを一つ増やす毎にS”の階数も一つ増える。この事は、各サブアレイは互いに独立な信号集合であり、サブアレイの数Qを一つ増やす事で独立な信号集合が一つ増えるので、行列S”の階数も一つ増える、と解釈できる。
なお、Q=M0-M+1の関係とMUSIC法のもう一つの適用条件MN≧D+1も含めて考えると、必要となる周波数の個数M0の条件は、以下の式(34)で与えられる。すなわち、必要となる周波数の個数M0は、検知すべき位置の数Dに比例して増大する。
非特許文献1では、一般的なアレイアンテナの相関行列に対しMUSIC法を適用する事で、到来方向推定と行っている。本実施の形態2では、式(21)で計算した全サブアレイ平均の相関行列Rx’に対し、(形式的に一般的なアレイアンテナに適用するのと同じ方式で)MUSIC法を適用する事で、所望信号ベクトルsの強度のx依存性(即ち、σ(x))を反映した評価関数PMU(x)が計算される。この時、評価関数PMU(x)は、以下の式(35)で与えられる。
ここで、a(x)は式(34)で定義された方向行列Aの列ベクトルである。また、ENは以下の式(36)で与えられる。
ここでベクトルek(k=D+1,D+2,・・・,MN)は、相関行列Rx’の固有ベクトルの内、その固有値がノイズ電力に等しいものである。MUSIC法によれば、式(35)の評価関数PMU(x)は、対象物1001d(d=1,2,・・・,D)の位置xdにおいてピークを与える。
したがって、評価関数PMU(x)がピーク値を与える位置xから、対象物1001d(d=1,2,・・・,D)の位置xdを割り出す事ができる。MUSIC法を適用する場合、(MN-D)個のノイズ空間の固有ベクトル{eD+1,eD+2,・・・,eMN}が利用されるが、それが最低1個必要であるので、MN-D≧1、すなわちMN≧D+1を満たす必要がある。
上述の例では、MUSIC法を用いて対象物1001d(d=1,2,・・・,D)の位置xdが検知されている。但し、本実施の形態2では、相関行列Rx’に対し、(形式的に一般的なアレイアンテナに適用するのと同じ方式で非特許文献1に記載の)ビームフォーマ法、Capon法、線形予測法を適用する事で、所望信号ベクトルs(t)の強度のx依存性(即ち、σ(x))を反映した評価関数を計算する事もできる。本実施の形態2におけるビームフォーマ法に基づく評価関数PBF(x)は、以下の式(37)で与えられる。
また、本実施の形態2におけるCapon法に基づく評価関数PCP(x)は、以下の式(38)で与えられる。
また、本実施の形態2における線形予測法に基づく評価関数PLP(x)は、以下の式(39)で与えられる。
上記の評価関数PBF(x),PCP(x),PLP(x)も、MUSIC法により得られる評価関数PMU(x)と同じく物対象物1001d(d=1,2,・・・,D)の位置xdにおいてピーク値を取る。従って、評価関数がピーク値を与える位置xから、対象物1001d(d=1,2,・・・,D)の位置xdを割り出す事ができる。
上述の本実施の形態2において開示した処理、即ち、反射波の測定結果から評価関数を算出し、その評価関数から対象物の位置を割り出す処理は、図9に示したデータ処理部1093によって実行される。また、本実施の形態2における評価関数を算出して評価関数のピークを探索する過程は、実施の形態1における、移相器1031と加算器1032とによる制御を行って、受信信号強度が最大になるビーム方向を探索する過程に対応している。
また、本実施の形態2では、送信部と受信部とを結ぶ方向の座標(即ち、x軸)の位置情報xd(即ち、一次元方向の位置)のみを検知する事ができる。何故なら、送信部と受信部とを備えた物体検知装置は、x軸を軸とした回転対称性があるため、対象物1001のx軸以外の座標値が異なっていても区別ができないからである。X軸以外の座標の位置情報も検知する方法については、実施の形態3において後述する。
続いて、図13を用いて、本実施の形態2における物体検知装置の動作について説明する。図13は、本発明の実施の形態2における物体検知装置の動作を示すフロー図である。また、本実施の形態2においても、物体検知装置を動作させることによって、物体検知方法が実施される。よって、本実施の形態2における物体検知方法の説明は、以下の物体検知装置1000の動作説明に代える。
図13に示すように、最初に、物体検知装置において、送信部が対象物に向け周波数を変化させながら電波を照射する(ステップB1)。
次に、複数の受信部それぞれは、対象物からの各周波数の反射波を、対応する受信アンテナによって受信する(ステップB2)。各受信アンテナは、送信部から見て一つの方向に配置されている。
次に、データ処理部は、q番目からq+M番目までのサンプリング時間の受信信号を用いて相関行列Rxq(q=1,2,・・・,Q,Q=M0-M+1)を計算する(ステップB3)。
次に、データ処理部は、計算したQ個の相関行列Rxq(q=1,2,・・・,Q)を平均した相関行列Rx’を計算し(ステップB4)、更に、相関行列Rx’から対象物の位置を反映する評価関数を計算する(ステップB5)。
その後、データ処理部は、評価関数のピークから対象物の位置を算出する(ステップB6)。算出結果は、出力部に出力される。
以上のように、本実施の形態2によれば、大量の受信アンテナを用意することなく、対象物の一次元の方向を推定することができる。また、本実施の形態2においても、実施の形態1で述べた効果を得ることができる。
(実施の形態3)
続いて、本発明の実施の形態3における物体検知装置及び物体検知方法について、図14~図18を参照しながら説明する。
続いて、本発明の実施の形態3における物体検知装置及び物体検知方法について、図14~図18を参照しながら説明する。
本実施の形態3は、実施の形態1で示した物体検知装置1000による仮想アレイの概念に基づいて、対象物の配置及び形状を識別するための二次元画像を生成する例を示している。従って、本実施の形態3においても、物体検知装置は、図1及び図9に示した、送信部1091、受信部1092、及びデータ処理部1093を備えている。但し、本実施の形態3は、受信部1092の個数の点で、実施の形態1と異なっている。以下、具体的に説明する。
図14は、本発明の実施の形態3における物体検知装置の構成及び動作原理を示す図である。また、図14には、各アンテナと対象物との位置関係が示されている。まず、本実施の形態3における物体検知装置においては、受信アンテナ1004は、送信部の送信アンテナ1003を基準にしたN(N=2,3,・・・)方向に沿って設置されている。また、データ処理部は、複数の受信部それぞれが生成したベースバンド信号の積を算出し、算出した積に基づいて、N方向を座標軸とするN次元の座標空間における、対象物1001の位置を検知する。
具体的には、送信アンテナ1003が座標の原点の位置に設置され、受信部の受信アンテナ1004(x)と受信アンテナ1004(y)とがそれぞれx軸上とy軸上とに設置されている。この場合、N=2である。
本実施の形態3では、送信アンテナ1003と受信アンテナ1004(x)を結ぶ方向と、送信アンテナ1003と受信アンテナ1004(y)を結ぶ方向とが、互いに異なる方向である事(平行でない事)が、2次元の画像を得る上で望ましい態様である。なお、送信アンテナ1003と受信アンテナ1004(x)を結ぶ方向と、送信アンテナ1003と受信アンテナ1004(y)を結ぶ方向とが直交している必要は必ずしも無い。
送信アンテナ1003からRF信号(電波)1010が、焦平面1002上に存在する対象物1001に向けて照射される。RF信号1010が対象物1001に照射された後、対象物1001からの反射波1007(x)と反射波1007(y)とが、それぞれ受信アンテナ1004(x)と受信アンテナ1004(y)とにおいて受信される。本実施の形態3においても、実施の形態1及び2と同様に、送信アンテナ1003が出力するRF信号1010のキャリア周波数は、時間の経過と共に連続的に変化する。
また、図14に示す実施の形態3は、図26で示したMills cross法における2つのアレイアンテナ201を、それぞれ周波数仮想アレイに置換する事を意図している。具体的には、図27で示したMills cross法における2つのアレイアンテナ201は、送信アンテナ1003と受信アンテナ1004(x)との組で構成される仮想アレイと、送信アンテナ1003と受信アンテナ1004(y)との組で構成される仮想アレイとに置換される。従って、本実施の形態3において、2次元画像の生成に必要なアンテナの数は最小の場合3本で良いという事になる。
次に、図15及び図16を用いて、本実施の形態3における物体検知装置による2次元画像の生成処理の詳細について説明する。図15及び図16は、本発明の実施の形態3における2次元周波数仮想アレイの相関行列の計算方法を説明する説明図である。また、図15及び図16では、2次元画像生成の動作解析のための計算モデルが示されている。
図15及び図16に示すように、本実施の形態3の計算モデルでは、x軸に1本の送信アンテナ1003(x0)とN本の受信アンテナ1004(x1)、・・・、1004(xN)が設置されている。更に、本実施の形態3の計算モデルでは、y軸にも、1本の送信アンテナ1003(y0)とN本の受信アンテナ1004(y1)、・・・、1004(yN)が設置されている。
xyz軸座標で、x軸上の送信アンテナ1003(x0)の位置を(dx0,0,0),n番目の受信アンテナ1004(xn)の位置を(dxn,0,0)とする。また、y軸上の送信アンテナ1003(y0)の位置を(0,dy0,0),n番目の受信アンテナ1004(yn)の位置を(0,dyn,0)とする。
また、対象物1001はz=z0の平面上でD個の位置(x1,y1,z0),(x2,y2,z0),・・・,(xD,yD,z0)に設置されるものとする。説明を簡単にするため、物体検知装置(送信アンテナ1003及び受信アンテナ1004)と対象物1001との位置関係は上記の位置関係に固定されているものとする。
また理論計算上では、図15で示すようにx軸上の送信アンテナ1003(x0)が送信している時はx軸上の受信アンテナ1004(x1)、・・・、1004(xN)のみが受信を行い、y軸上の送信アンテナ1003(y0)が送信している時はy軸上の受信アンテナ1004(y1)、・・・、1004(yN)のみが受信を行うものとする。
また、図15及び図16の例では、x軸とy軸とで別々に送信アンテナ1003(x0)と送信アンテナ1003(y0)とが配置されているが、これは理論の説明に一般性を持たせるためである。実用上は、送信アンテナ1003(x0)と送信アンテナ1003(y0)とは、一本の送信アンテナで構成されていても良く、その場合は、この一本の送信アンテナが送信している時に、x軸上にある受信アンテナとy軸上にある受信アンテナとで同時に受信を行えばよい。
また、本実施の形態3においても、実施の形態1及び2と同じく、送信アンテナ1003(x0)及び送信アンテナ1003(y0)は、M個のキャリア周波数αt1,αt2,・・・,αtMのRF信号1010を送信する。RF信号1010の変調は、本実施の形態3においても、上述したFM-CW方式によって行なわれる。
対象1001d(d=1,2,・・・,D)から反射され、x軸上のn番目の受信アンテナ1004(xn)で受信されたRF信号1007のサンプリング時間tmにおける複素振幅をsxn(xd,yd,tm)とする。また、x軸上のn番目の受信アンテナ1004(xn)で実際に測定される受信信号(各対象からの反射波の合成)の複素振幅をsx(tm)とする。sxn(tm)とsxn(xd,yd,tm)との間には、以下の式(40)に示す関係がある。
また、y軸上のn番目の受信アンテナ1004(yn)についても同様の信号syn(tm)およびsyn(xd,yd,tm)を定義すると、この場合も、以下の式(41)に示すように、式(29)と同様の関係が成立する。
また、対象物1001dとx軸上の送信アンテナ1003(x0)との距離Lxo(xd,yd)は、以下の式(42)で与えられる。また、対象物1001dとx軸上のn番目の受信アンテナ1004(x0)との距離Lxn(xd,yd)は、以下の式(43)で与えられる。
y軸上の送信アンテナ1003(y0)およびn番目の受信アンテナ1004(yn)に関しても同様に対象1001dとの距離をそれぞれLyo(xd,yd)とLyn(xd,yd)とすると、それらは以下の式(44)と(45)で与えられる。
送信アンテナ1003(x0)から送出されるRF信号の複素振幅s0と、x軸上のn番目の受信アンテナ1004(xn)で受信されたRF信号から得られる複素振幅sx(xd,yd,tm)との間には、以下の式(46)に示す関係がある。
式(46)において、σ(xd,yd)は、対象物1001d(d=1,2,・・・,D)の反射率を表す未知数である。また、y軸上の受信アンテナ1004(yn)についても、以下の式(47)に示すように、同様の関係が成立する。
式(47)を式(40)に代入することで、以下の式(48)が得られ、式(48)を式(41)に代入する事で、以下の式(49)が得られる。
次に、x軸上のn番目の受信アンテナ1004(xn)(n=1,2,・・・,N)における測定信号sxn(tm)を用いて、以下の式(50)に示すように、測定信号ベクトルsxを定義する。
y軸方向の受信アンテナ1004(yn)(n=1,2,・・・,N)における測定信号についても、以下の式(51)に示すように同様に定義する。
次に、Mills cross法の手法に従い、上記の式(50)のx軸方向測定ベクトルsxの要素と上記の式(51)のy軸方向測定ベクトルsyの要素との全ての組み合わせについて積を算出すると、以下の式(52)に示す直積ベクトルsxyが生成される。なお、ここでいう「積」は、上述した「ベースバンド信号の積」に該当する。
式(52)において、nとvはそれぞれx方向とy方向に配置されたアンテナ番号、mとwはそれぞれx方向とy方向に配置されたアンテナで受信した信号の周波数番号を表す添字である。次に、以下の式(53)によって方向行列Aを定義する。
式(53)において、方向行列Aのサイズは(MN)2×D、行列AnvのサイズはM2×D、ベクトルanv(xd,yd)のサイズはM2×1となる。行列Anvは、n番目のx方向アンテナ1004(xn)とv番目のy方向アンテナ1004(yv)が関与する方向行列である。系全体の方向行列Aは、全てのアンテナ番号の組(n,v)の方向行列Anvをまとめたものとなる。
ここで上述した1次元到来方向推定の場合と同じく、複素振幅s0と反射率σ(xd,yd)を用いて、以下の式(54)によって所望信号ベクトルsを定義する。
式(48)及び(49)から、式(52)の測定信号ベクトルsxy(t)と、式(53)の方向行列Aと、式(54)の所望信号ベクトルsとの間には、以下の式(55)に示す関係式が得られる。式(55)では、ノイズ(乱数)を要素とするベクトルn(t)を付加している。
次に、測定で得た式(52)の測定信号ベクトルsxyを用いて、相関行列Rxyを計算する。式(55)の関係から、相関行列Rxyと方向行列Aの関係は、以下の式(56)で与えられる。
式(56)において、PNはノイズ項n(t)の平均電力、Iは(MN)2×(MN)2次の単位行列である。相関行列Rxy、行列A、行列Sのサイズはそれぞれ(MN)2×(MN)2次、(MN)2×D次、D×D次となる。
式(55)と式(56)とは、実施の形態2で議論した1次元到来方向推定における式(26)と式(28)と同型である。よって、1次元到来方向推定と同じ手順で相関行列Rxyに対しMUSIC法を適用する事で、σ(xd,yd)を反映した評価関数PMU(x,y)を計算できる。
但し、1次元到来方向推定の場合と同じく、MUSIC法の適用条件として式(56)内の行列Aと行列Sとがフルランクである事が要求される。そして上述の説明と同じく、方向行列Aはフルランクであるが、行列Sはσ(xi)=σ(xj)(i≠j)となる場合においてフルランクではない。そのため、サブアレイ法によって行列Sがフルランクになるように処理を行う必要がある。
2次元画像生成の場合においても、本実施の形態2で説明した1次元到来方向推定におけるサブアレイ法と同じ手順で、M個の周波数で一つのサブアレイが構築され、合計Q個のサブアレイが構築される。全体のサンプリング時間の個数をM0とすると、Q=M0-M+1の関係がある。q番目のサブアレイ信号は、以下の式(57)によって定義される。信号ベクトルsxyの成分sxy(nv)(mw)のサンプリング時間を表す添字mとwとを同時に+(q-1)個シフトしたものが、q番目のサブアレイ信号となる。
式(57)のサブアレイ信号sxy
qと、式(42)の方向行列との間には、以下式(58)に示す関係式が成立する。
サブアレイqの相関行列Rx
qは、以下の式(59)のように計算される。
式(59)において、相関行列Rxy
q、行列A’、行列S’のサイズは、それぞれ(NM)2×(NM)2次、(NM)2×N2D次、N2D×N2D次となる。次に全てのサブアレイq(q=1,2,・・・,Q)の相関行列の平均Rxy’を計算する。全サブアレイ平均の相関行列Rxy’と方向行列A’との関係は以下の式(60)のように計算される。
上述の実施の形態2で示した1次元到来方向推定の場合と同様にして、以下の事が示される。
(1)行列A’とS”がフルランクであれば、相関行列Rxy’にMUSIC法を適用してσ(xd,yd)を反映した評価関数PMU(x,y)を計算できる。
(2)行列A’については、方向行列A11,A12,・・・,A1N,・・・,AN1,・・・,ANNはそれぞれ独立かつフルランクであるので、式(59)で与えられるA’もまたフルランクである。
(3)行列S’’は、Q≧Dであればフルランクとなる。1次元到来方向推定におけるMUSIC法の適用条件MN≧D+1は、2次元画像生成では(MN)2≧D+1となる。これとサブアレイにおける条件Q=M0-M+1とQ≧Dを考慮すると、必要となるサンプリング時間の個数(周波数の数)M0の条件は、以下の式(61)で与えられる。即ち、必要となるサンプリング時間の個数M0は、検知すべき位置の数Dに概ね比例して増大する。
(1)行列A’とS”がフルランクであれば、相関行列Rxy’にMUSIC法を適用してσ(xd,yd)を反映した評価関数PMU(x,y)を計算できる。
(2)行列A’については、方向行列A11,A12,・・・,A1N,・・・,AN1,・・・,ANNはそれぞれ独立かつフルランクであるので、式(59)で与えられるA’もまたフルランクである。
(3)行列S’’は、Q≧Dであればフルランクとなる。1次元到来方向推定におけるMUSIC法の適用条件MN≧D+1は、2次元画像生成では(MN)2≧D+1となる。これとサブアレイにおける条件Q=M0-M+1とQ≧Dを考慮すると、必要となるサンプリング時間の個数(周波数の数)M0の条件は、以下の式(61)で与えられる。即ち、必要となるサンプリング時間の個数M0は、検知すべき位置の数Dに概ね比例して増大する。
次に、式(60)で計算した全サブアレイ平均の相関行列Rxy’にMUSIC法を適用する事で、σ(xd,yd)を反映した評価関数PMU(x,y)を計算する。その結果、以下の式(62)に示す評価関数が得られる。
ここで、a(x,y)は式(42)で定義された方向行列Aの列ベクトルである。また、ENは以下の式(63)で与えられる。
ここでベクトルek(k=D+1,D+2,・・・,(MN)2)は、相関行列Rsxy’の固有ベクトルの内、その固有値がノイズ電力に等しいものである。
評価関数PMU(x,y)は対象物1001dの位置(xd,yd)(d=1,2,・・・,D)においてピークを与える。従って、評価関数PMU(x,y)から対象物1001dの位置情報(xd,yd)(d=1,2,・・・,D)を検知し、そこから対象物1001の分布ないし形状を検知する事ができる。
上記ではMUSIC法を用いて対象物1001d(d=1,2,・・・,D)の位置を検知したが、相関行列Rsxy’に対し、(形式的に一般的なアレイアンテナに適用するのと同じ方式で非特許文献1に記載の)ビームフォーマ法、Capon法、線形予測法を適用する事で、各方式の評価関数を計算する事もできる。
上記の考察に従い、本実施の形態3におけるビームフォーマ法に基づく評価関数PBF(x,y)は、以下の式(64)で与えられる。
また、本実施の形態3におけるCapon法に基づく評価関数PCP(x,y)は、以下の式(65)で与えられる。
また、本実施の形態3における線形予測法に基づく評価関数PLP(x,y)は、以下の式(66)で与えられる。
上記の評価関数PBF(x,y),PCP(x,y),PLP(x,y)も、MUSIC法により得られる評価関数PMU(x,y)と同じく物対象物1001d(d=1,2,・・・,D)の位置(xd,yd)においてピーク値を取る。従って、評価関数がピーク値を与える位置(x,y)から、対象物1001d(d=1,2,・・・,D)の位置xdを割り出す事ができる。
本実施の形態3において開示された処理、即ち、反射波の測定結果から評価関数を算出し、その評価関数から対象物の位置を割り出す処理も、実施の形態2と同様に、図9に示したデータ処理部1093によって実行される。また、本実施の形態3における評価関数を算出して評価関数のピークを探索する過程は、実施の形態1における、移相器1031と加算器1032とによる制御を行って、受信信号強度が最大になるビーム方向を探索する過程に対応している。
続いて、図17を用いて、本実施の形態3における物体検知装置の動作について説明する。図17は、本発明の実施の形態3における物体検知装置の動作を示すフロー図である。また、本実施の形態3においても、物体検知装置を動作させることによって、物体検知方法が実施される。よって、本実施の形態3における物体検知方法の説明は、以下の物体検知装置1000の動作説明に代える。
図17に示すように、最初に、物体検知装置において、送信部が対象物に向けRFキャリア周波数を変化させながら電波を照射する(ステップC1)。
次に、複数の受信部それぞれは、対象物からの反射波を、対応する受信アンテナによって受信する(ステップC2)。各受信アンテナは、送信部から見て2つの方向に配置されている。
次に、データ処理部は、q番目からq+M番目までのサンプリング時間の受信信号を用いて相関行列Rxy
q(q=1,2,・・・,Q,Q=M0-M+1)を計算する(ステップC3)。
次に、データ処理部は、計算したQ個の相関行列Rxy
q(q=1,2,・・・,Q)を平均した相関行列Rxy’を計算し(ステップC4)、更に、相関行列Rxy’から対象物1001の位置を反映する評価関数を計算する(ステップC5)。
その後、データ処理部は、評価関数のピークから対象物の位置を算出し、更に、対象物の配置及び形状を、二次元画像として出力部に出力する(ステップC6)。
ここで、図18を用いて、本実施の形態3における物体検知装置により得られる2次元画像の例について説明する。図18は、本発明の実施の形態3における物体検知装置から出力された画像の一例を示す図である。
図18の例では、対象物1001は、(x,y,z)座標表示で、(-20cm,-20cm,100cm)、(0cm,0cm,100cm)、(20cm,20cm,100cm)の3箇所に配置されているとする。
また、送信アンテナ1003は(-100cm,-100cm,0cm)の位置に配置されている。受信アンテナ1004は(0cm,-100cm,0cm)の位置と、(-100cm,0cm,0cm)の位置とに配置されているとする。
また、送信アンテナ1003は、76GHzから81GHzの間(帯域幅BW=5GHz)でキャリア周波数を変化させるFM-CW変調を掛けたRF信号1010を対象物1001に向けて照射しているとする。1チャープ周期(Tchirp)内におけるサンプリング時間の個数(全ての周波数の数)M0は21、サブアレイの数Qは10、一つのサブアレイあたりの個数(周波数の数)Mは12であるとする。サンプリング周期Δtとサンプリング時間の周波数変化率αは、αΔt=250MHzとなるように設定している。このような条件下においては、図18に示すように、実際に3箇所に配置されている対象物1001が検知される。そして、各対象物1001は、2次元画像上で表示される。
また、図14~図18の例は、受信アンテナが2方向(N=2)に沿って設置されている場合について説明しているが、本実施の形態3は、受信アンテナが3方向(N=3)以上に沿って設置されている場合にも、適用できる。特に、受信アンテナが直交する3方向に沿って設置されている場合は、対象物の3次元空間での位置を特定できる。
(実施の形態4)
続いて、本発明の実施の形態4における物体検知装置及び物体検知方法について、図19~図21を参照しながら説明する。
続いて、本発明の実施の形態4における物体検知装置及び物体検知方法について、図19~図21を参照しながら説明する。
図19は、本発明の実施の形態4における物体検知装置の概略構成を示す図である。図18に示すように、本実施の形態4では、物体検知装置1200は、複数の物体検知ユニット1202p(p=1,2,・・・,P)を備えている。また、各物体検知ユニット1202p(p=1,2,・・・,P)は、送信部1091pと受信部1092pとの組を備えている。ここで、Pは物体検知ユニット1202の個数を示している。
そして、物体検知ユニット1202pにおいて、送信部1091pは対象物1201p1、1201p2、・・・、1201pQに電波を照射し、受信部1092pは対象物1201p1、1201p2、・・・、1201pQからの反射波を受信する。これにより、対象物1201p1、1201p2、・・・、1201pQの状態が検知される。Qは対象物1201の個数である。
また、物体検知装置1200は、対象物1201p1、1201p2、・・・、1201pQが人である場合は、人(1201p1、1201p2、・・・、1201pQ)が着用している衣服を透過した電波によって、衣服下にある物品の存在を検知することができる。
更に、物体検知装置1200は、対象物1201p1、1201p2、・・・、1201pQが物体(特に誘電体)である場合は、物体(1201p1、1201p2、・・・、1201pQ)を透過した電波によって、物体(1201p1、1201p2、・・・、1201pQ)の内部構造を検知することができる。
物体検知装置1200は、物体検知ユニット1202pにより、対象物が流れ作業の対象である場合に、対象物1201p1、1201p2、・・・、1201pQの状態を順番に検知することもできる。
また、図19の例では、一つの対象物1201の検知又は検査に対して、一つの物体検知ユニット1202pが割り当てられている。但し、本実施の形態4は、これに限定されず、一つの対象物1201の検知又は検査に対して、複数の物体検知ユニット1202pが割り当てられていても良い。また、本実施の形態4では、複数の対象物1201の検知又は検査に対し、一つの物体検知ユニット1202pが割り当られていても良い。
このように、本実施の形態では、物体検知ユニット1202は小型かつ低コストで実現されるので、物体検知ユニット1202の個数Pを容易に増やす事ができる。従って、図19の例で示した物体検知装置1200においては、対象物1201p1、1201p2、・・・、1201pQ(p=1,2,・・・,P)の検査速度を、物体検知ユニット1202の個数Pに比例して引き上げる事ができる。
ところで、図19に示す物体検知装置1200においては、物体検知ユニット1202p(p=1,2,・・・,P)それぞれ同士での干渉によって、誤動作が生じる可能性がある。即ち、送信部1091pから受信部1092r(p≠r)への電波の回り込みが、誤動作を生じさせる干渉の要因になる。この問題を回避するための構成及び動作について、図20を用いて説明する。
図20は、本発明の実施の形態4における物体検知装置の構成を具体的に示すブロック図である。図20に示すように、本実施の形態4における物体検知装置1200は、複数の物体検知ユニット1200pに加えて、データ制御部1203を備えている。データ制御部1203は、物体検知ユニット1202p(p=1,2,・・・,P)を構成する送信部1091pと受信部1092pとを制御する。具体的には、データ制御部1203は、各物体検知ユニット1202pに対して、実施の形態1~3におけるデータ制御部1093と同様の処理を行なう。
また、データ制御部1203は、物体検知ユニット毎に、使用される電波の周波数が異なるように、物体検知ユニットそれぞれを動作させる。具体的には、データ制御部1203は、物体検知ユニット1202pのRF周波数fpと物体検知装置1202rのRF周波数frとを異なる値とする制御を行なっている(p,r=1,2,・・・,P,かつp≠r)。そして、このような制御により、互いに異なる物体検知ユニット1202pと物体検知ユニット1202r(p≠r)とは、異なるRF周波数で動作する。このため、物体検知ユニット1202pと物体検知ユニット1202r(p≠r)との間での干渉の発生が抑制される。
ここで、図21を用いて、各物体検知ユニット1202p(p=1,2,・・・,P)におけるRF周波数の制御について説明する。図21は、本発明の実施の形態4において行なわれる周波数制御の一例を示す図である。
図3に示したように、本実施の形態4においても、実施の形態1~3と同様に、各物体検知ユニット1202p(p=1,2,・・・,P)では、サンプリング時間t1、t2、・・・、tMにおいて、キャリア周波数(RF周波数)fはfminからfmaxまで連続的に変化する。つまり、各物体検知ユニット1202pにおいて、RF周波数の時間変化がチャープ状に制御される。
そして、図21に示すように、本実施の形態4では、データ制御部1203は、各物体検知ユニット1202pのRF周波数の時間変化が互いにずれるように制御を行なう。この結果、異なる物体検知ユニット1202pと物体検知ユニット1202r(p≠r)とが、同じRF周波数で動作する事は無い。
このように、本実施の形態4では、各物体検知ユニット1202pに対して、上述したデータ処理部として機能しながら、各各物体検知ユニット1202pにおけるRF周波数が異なるように制御を行なっている。具体的には、データ制御部1203は、図10に示したステップA1~A5を実行し、その際、ステップA2において、物体検知ユニット毎に、異なるαtmを設定する。
(実施の形態による効果)
以下において、本実施の形態における効果を要約する。一般的なアレイアンテナ方式と本実施の形態1~4とを比較した場合、アレイアンテナ方式は多数のアンテナを必要とする。その一方で、本実施の形態1~4では、実際のアンテナの数を増やす代わりに周波数の数を増やす事で仮想的なアンテナを増やす事ができる。その結果、本実施の形態では少なくとも1本の送信アンテナと1方向あたり1本の受信アンテナで、一般的なアレイアンテナ方式と同等の機能を実装でき、実際のアンテナ本数を一般的なアレイアンテナ方式と比べて大幅に削減できる。
以下において、本実施の形態における効果を要約する。一般的なアレイアンテナ方式と本実施の形態1~4とを比較した場合、アレイアンテナ方式は多数のアンテナを必要とする。その一方で、本実施の形態1~4では、実際のアンテナの数を増やす代わりに周波数の数を増やす事で仮想的なアンテナを増やす事ができる。その結果、本実施の形態では少なくとも1本の送信アンテナと1方向あたり1本の受信アンテナで、一般的なアレイアンテナ方式と同等の機能を実装でき、実際のアンテナ本数を一般的なアレイアンテナ方式と比べて大幅に削減できる。
合成開口レーダー方式と本実施の形態とを比較した場合、合成開口レーダー方式は受信機を機械的に動かす必要があり、これが物体の検知及び検査のための時間が長くなるという問題があった。一方、本実施の形態では、受信機の位置ではなく受信周波数を電子的に走査すればよいので、合成開口レーダー方式に比べて物体の検知及び検査のための時間を短縮できる。
すなわち、本実施の形態における物体検知装置及び物体検知方法においては、一般的なアレイアンテナ方式よりも必要なアンテナおよびそれに付随する受信機の数を削減する事ができるので、装置のコスト、サイズ、重量を削減できるという効果を奏する。また、本実施の形態における物体検知装置及び物体検知方法においては、一般的な合成開口レーダー方式と異なり、装置を機械的に動かす必要がないため、物体検知及び検査の時間を短縮できるという効果も奏する。
本実施の形態では、サンプリング時間毎に異なるRF周波数の電波を検知対象物に照射し、対象物で反射された電波、又は対象物から放射される電波を検知する事で、検知対象物の画像を生成することができる。従って、本実施の形態によれば、従来よりも必要なアンテナおよび受信部の数を減らし、かつ移動させる必要も無く、高速な走査による画像生成を実現することができる。
以上、実施の形態を参照して本願発明を説明したが、本願発明は上記実施の形態に限定されるものではない。また、上述の各特許文献等に開示されている内容は、本願発明に引用をもって繰り込むことも可能とする。本願発明の全開示(特許請求の範囲を含む)の枠内において、さらにその基本的技術思想に基づいて、実施の形態の変更・調整が可能である。また、本願発明の特許請求の範囲の枠内において種々の開示要素の多様な組み合わせあるいは選択も可能である。すなわち、本願発明は、特許請求の範囲を含む全開示、技術的思想にしたがって、当業者であればなし得ることが可能な各種変形、修正を含むことは勿論である。
以上のように本発明によれば、電波を用いた物体の検知において、精度を向上させつつ、装置コスト、サイズ、及び重量の増大化を抑制することができる。本発明は、衣服の下に隠されている物品又は鞄の中の物品等を画像化して検査する場合に有用である。
1000 物体検知装置
1001、1201 対象物(検知対象となる物体)
1002 焦平面
1003 送信アンテナ
1004 受信アンテナ
1007、1010 電波(RF信号)
1041 低雑音増幅器
1042 ミキサ
1043 フィルタ
1044 アナログ-デジタル変換器
1075 カプラ
1091 送信部
1092 受信部
1093 データ受信部
1094 出力部
1103 発振器
1102 受信制御部
1104 送信制御部
1202 物体検知ユニット
1200 物体検知装置(実施の形態4)
1203 データ制御部
1001、1201 対象物(検知対象となる物体)
1002 焦平面
1003 送信アンテナ
1004 受信アンテナ
1007、1010 電波(RF信号)
1041 低雑音増幅器
1042 ミキサ
1043 フィルタ
1044 アナログ-デジタル変換器
1075 カプラ
1091 送信部
1092 受信部
1093 データ受信部
1094 出力部
1103 発振器
1102 受信制御部
1104 送信制御部
1202 物体検知ユニット
1200 物体検知装置(実施の形態4)
1203 データ制御部
Claims (10)
- 電波によって物体を検知するための物体検知装置であって、
時間の経過と共に周波数が連続的に変化する電波を、送信信号として放射する、送信部と、
前記送信信号を取得し、前記物体からの前記電波を受信信号として受信し、更に、受信した前記受信信号に、取得した前記送信信号をミキシングしてベースバンド信号を生成する、受信部と、
サンプリング時間毎の前記ベースバンド信号の測定値から、前記電波の到来方向を推定し、推定した前記電波の到来方向に基づいて、前記電波の強度分布を特定し、特定した前記強度分布に基づいて、前記物体を検知する、データ処理部と、
を備えている、ことを特徴とする物体検知装置。 - 前記受信部が、前記送信部に接続されたミキサと、フィルタとを備え、前記ミキサによって、前記受信信号に、取得した前記送信信号をミキシングし、前記フィルタによって、ミキシングで得られた信号の所望周波数以外の成分を除去することによって、前記ベースバンド信号を生成する、
請求項1に記載の物体検知装置。 - 前記送信部が、送信アンテナを備え、
前記受信部が、前記送信アンテナを基準にした一方向に沿って設置された受信アンテナを備え、
前記データ処理部が、前記強度分布に基づいて、前記物体の前記一方向における位置を検知する、
請求項1または2に記載の物体検知装置。 - 前記送信部が、送信アンテナを備え、
前記受信部が、複数備えられ、複数の前記受信部それぞれは受信アンテナを備え、
前記受信アンテナは、前記送信アンテナを基準したN方向に沿って設置され、
前記データ処理部は、複数の前記受信部それぞれが生成した前記ベースバンド信号の積を算出し、算出した積に基づいて、前記N方向を座標軸とするN次元の座標空間における、前記物体の位置を検知する、
請求項1または2に記載の物体検知装置。 - 前記データ処理部が、サンプリング時間毎の前記ベースバンド信号の測定値から相関行列を算出し、更に前記相関行列から前記物体の位置を反映する評価関数を求め、求めた前記評価関数から前記物体の画像を生成する、
請求項1~4のいずれかに記載の物体検知装置。 - 前記データ処理部が、前記サンプリング時間の範囲が異なる前記ベースバンド信号の測定値から、各々の前記サンプリング時間の範囲に対応する前記相関行列を算出し、更に各々の前記サンプリング時間の範囲に対応する前記相関行列の平均値を算出し、更に前記相関行列の平均値に基づいて、前記物体の位置を反映する評価関数を求め、求めた前記評価関数から前記物体の画像を生成する、
請求項5に記載の物体検知装置。 - 前記データ処理部が、サンプリング時間毎の前記ベースバンド信号の測定値にノイズを付加し、前記ベースバンド信号の測定値に前記ノイズを付加した信号から、前記相関行列を算出する。
請求項5~6に記載の物体検知装置。 - 前記送信部において、予め設定された可視領域と分解能との比、又は予め設定された画素数に応じて、送信される信号の時間長、又はサンプリング時間の数が選択される、
請求項1~7のいずれかに記載の物体検知装置。 - 前記送信部及び前記受信部が、それぞれ複数備えられ、
少なくとも一つの前記送信部と少なくとも一つの前記受信部とが組となって物体検知ユニットを構成し、
前記データ処理部は、前記物体検知ユニット毎に、前記電波の周波数が異なるように、前記物体検知ユニットそれぞれを動作させる、
請求項1~8のいずれかに記載の物体検知装置。 - 電波によって物体を検知するための方法であって、
(a)送信機によって、時間の経過と共に周波数が連続的に変化する電波を、送信信号として放射する、ステップと、
(b)受信機によって、前記送信信号を取得し、前記物体からの前記電波を受信信号として受信し、更に、受信した前記受信信号に、取得した前記送信信号を加算して、ベースバンド信号を生成する、ステップと、
(c)データ処理装置によって、サンプリング時間毎の前記ベースバンド信号の測定値から、前記電波の到来方向を推定し、推定した前記電波の到来方向に基づいて、前記電波の強度分布を特定し、特定した前記強度分布に基づいて、前記物体を検知する、ステップと、
を有する、ことを特徴とする物体検知方法。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018531727A JP6911861B2 (ja) | 2016-08-05 | 2016-08-05 | 物体検知装置および物体検知方法 |
PCT/JP2016/073205 WO2018025421A1 (ja) | 2016-08-05 | 2016-08-05 | 物体検知装置および物体検知方法 |
US16/323,031 US20190187266A1 (en) | 2016-08-05 | 2016-08-05 | Object detection apparatus and object detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2016/073205 WO2018025421A1 (ja) | 2016-08-05 | 2016-08-05 | 物体検知装置および物体検知方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2018025421A1 true WO2018025421A1 (ja) | 2018-02-08 |
Family
ID=61073706
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2016/073205 WO2018025421A1 (ja) | 2016-08-05 | 2016-08-05 | 物体検知装置および物体検知方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20190187266A1 (ja) |
JP (1) | JP6911861B2 (ja) |
WO (1) | WO2018025421A1 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019180767A1 (ja) * | 2018-03-19 | 2019-09-26 | 日本電気株式会社 | 物体検知装置、物体検知方法、及びコンピュータ読み取り可能な記録媒体 |
JP2020122720A (ja) * | 2019-01-31 | 2020-08-13 | 古河電気工業株式会社 | 位置推定装置および位置推定方法 |
DE102020115998A1 (de) | 2020-06-17 | 2021-12-23 | Tk Elevator Innovation And Operations Gmbh | Aufzugsanlage |
JP7166496B1 (ja) * | 2022-02-24 | 2022-11-07 | 三菱電機株式会社 | レーダシステム |
JP7464126B2 (ja) | 2020-07-16 | 2024-04-09 | 日本電気株式会社 | 物体検知装置及び物体検知方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20180016820A (ko) * | 2016-08-08 | 2018-02-20 | 엘지전자 주식회사 | 재실 감지 장치 및 그 제어 방법 |
DE112020002763T5 (de) * | 2019-06-25 | 2022-02-24 | Murata Manufacturing Co., Ltd. | Empfänger, radarvorrichtung einschliesslich eines empfängers, fahrzeug einschliesslich eines empfängers und kommunikationssystem einschliesslich eines empfängers |
US11129042B2 (en) * | 2019-09-16 | 2021-09-21 | Nec Corporation Of America | Path crossing detection based on wireless signals travel time |
EP4067935A4 (en) * | 2019-11-26 | 2023-12-13 | Kyocera Corporation | ELECTRONIC DEVICE, ELECTRONIC DEVICE CONTROL METHOD AND PROGRAM |
CN111144267B (zh) * | 2019-12-23 | 2023-09-26 | 上海联影智能医疗科技有限公司 | 设备运行状态检测方法、装置、存储介质及计算机设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04313091A (ja) * | 1991-04-10 | 1992-11-05 | Honda Motor Co Ltd | Fmレーダ装置 |
JP2003185726A (ja) * | 2001-12-21 | 2003-07-03 | Toshiba Corp | 受信装置及び検出装置 |
JP2007285912A (ja) * | 2006-04-18 | 2007-11-01 | Toyota Central Res & Dev Lab Inc | 物体検知装置 |
JP2010071958A (ja) * | 2008-09-22 | 2010-04-02 | Denso Corp | レーダ装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013096887A (ja) * | 2011-11-02 | 2013-05-20 | Panasonic Corp | レーダ装置 |
US20130169471A1 (en) * | 2011-12-28 | 2013-07-04 | Hrl Laboratories, Llc | Coded aperture beam analysis method and apparatus |
JP6028388B2 (ja) * | 2012-05-11 | 2016-11-16 | 富士通株式会社 | 到来方向推定装置、及び到来方向推定方法 |
US9541638B2 (en) * | 2014-11-11 | 2017-01-10 | Nxp B.V. | MIMO radar system |
US10359511B2 (en) * | 2014-12-29 | 2019-07-23 | Sony Corporation | Surveillance apparatus having a radar sensor |
-
2016
- 2016-08-05 WO PCT/JP2016/073205 patent/WO2018025421A1/ja active Application Filing
- 2016-08-05 JP JP2018531727A patent/JP6911861B2/ja active Active
- 2016-08-05 US US16/323,031 patent/US20190187266A1/en not_active Abandoned
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04313091A (ja) * | 1991-04-10 | 1992-11-05 | Honda Motor Co Ltd | Fmレーダ装置 |
JP2003185726A (ja) * | 2001-12-21 | 2003-07-03 | Toshiba Corp | 受信装置及び検出装置 |
JP2007285912A (ja) * | 2006-04-18 | 2007-11-01 | Toyota Central Res & Dev Lab Inc | 物体検知装置 |
JP2010071958A (ja) * | 2008-09-22 | 2010-04-02 | Denso Corp | レーダ装置 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019180767A1 (ja) * | 2018-03-19 | 2019-09-26 | 日本電気株式会社 | 物体検知装置、物体検知方法、及びコンピュータ読み取り可能な記録媒体 |
JPWO2019180767A1 (ja) * | 2018-03-19 | 2021-02-18 | 日本電気株式会社 | 物体検知装置、物体検知方法、及びプログラム |
US11906650B2 (en) | 2018-03-19 | 2024-02-20 | Nec Corporation | Object detection apparatus, object detection method, and computer-readable recording medium |
JP2020122720A (ja) * | 2019-01-31 | 2020-08-13 | 古河電気工業株式会社 | 位置推定装置および位置推定方法 |
JP7251993B2 (ja) | 2019-01-31 | 2023-04-04 | 古河電気工業株式会社 | 位置推定装置および位置推定方法 |
DE102020115998A1 (de) | 2020-06-17 | 2021-12-23 | Tk Elevator Innovation And Operations Gmbh | Aufzugsanlage |
JP7464126B2 (ja) | 2020-07-16 | 2024-04-09 | 日本電気株式会社 | 物体検知装置及び物体検知方法 |
JP7166496B1 (ja) * | 2022-02-24 | 2022-11-07 | 三菱電機株式会社 | レーダシステム |
Also Published As
Publication number | Publication date |
---|---|
US20190187266A1 (en) | 2019-06-20 |
JPWO2018025421A1 (ja) | 2019-05-23 |
JP6911861B2 (ja) | 2021-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2018025421A1 (ja) | 物体検知装置および物体検知方法 | |
US11670866B2 (en) | Radar apparatus using transmitting array antenna and receiving array antenna | |
JP7001069B2 (ja) | 推論用知識生成装置、推論用知識生成方法、及びプログラム | |
US9121943B2 (en) | Beam forming device and method | |
JP5617334B2 (ja) | レーダ装置及び目標探知方法 | |
WO2017159521A1 (ja) | 物体検知装置および物体検知方法 | |
US20130002488A1 (en) | Wideband beam forming device; wideband beam steering device and corresponding methods | |
JP6838658B2 (ja) | 物体検知装置、物体検知方法、及びプログラム | |
US20220003834A1 (en) | Radar device and transmitting/receiving array antenna | |
RU2593595C1 (ru) | Способ измерения угловых координат в нелинейном радиолокаторе | |
JP6849100B2 (ja) | 物体検知装置、物体検知方法及びプログラム | |
US20220187424A1 (en) | Methods and systems for processing radar signals | |
JP6939981B2 (ja) | 物体検知装置、及び物体検知方法 | |
JP6415392B2 (ja) | 信号処理装置 | |
Leino et al. | Millimeter-wave imaging method based on frequency-diverse subarrays | |
JP2020056592A (ja) | レーダ装置 | |
US20240219517A1 (en) | Radar system | |
EP4343365A1 (en) | Radar system and inspection method | |
JP2023043026A (ja) | レーダ装置、送受信方法及びレーダシステム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 16911680 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2018531727 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 16911680 Country of ref document: EP Kind code of ref document: A1 |