WO2021250866A1 - レーダ装置およびレーダ画像生成方法 - Google Patents

レーダ装置およびレーダ画像生成方法 Download PDF

Info

Publication number
WO2021250866A1
WO2021250866A1 PCT/JP2020/023091 JP2020023091W WO2021250866A1 WO 2021250866 A1 WO2021250866 A1 WO 2021250866A1 JP 2020023091 W JP2020023091 W JP 2020023091W WO 2021250866 A1 WO2021250866 A1 WO 2021250866A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
antennas
transmission
transmitting
radar
Prior art date
Application number
PCT/JP2020/023091
Other languages
English (en)
French (fr)
Inventor
尭之 北村
啓 諏訪
Original Assignee
三菱電機株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to JP2020555255A priority Critical patent/JP6921339B1/ja
Priority to PCT/JP2020/023091 priority patent/WO2021250866A1/ja
Priority to DE112020007038.7T priority patent/DE112020007038T5/de
Publication of WO2021250866A1 publication Critical patent/WO2021250866A1/ja
Priority to US17/966,041 priority patent/US20230039504A1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/08Systems for measuring distance only
    • G01S13/32Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
    • G01S13/325Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of coded signals, e.g. P.S.K. signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/42Simultaneous measurement of distance and other co-ordinates
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/03Details of HF subsystems specially adapted therefor, e.g. common to transmitter and receiver
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves

Definitions

  • This disclosure relates to a radar device and a radar image generation method.
  • a body check device using a radar device is a device that images a subject by transmitting, for example, a radio wave in a millimeter wave band from the radar device and receiving a reflected wave reflected by the subject. Since the radio waves in the millimeter wave band pass through human clothes and the amount of radiation to humans is small, the body check device using the radio waves in the millimeter wave band can distinguish between metal and non-metal contraband.
  • Non-Patent Document 1 millimeter waves are transmitted in time divisions from a plurality of transmitting antennas while a subject is moving at a constant speed, and a two-dimensional receiving array receives the reflected waves of the millimeter waves. Describes a radar device that performs three-dimensional millimeter-wave imaging.
  • the radar device described in Non-Patent Document 1 adopts a time division multiplexing-multiple input multiple output (TDM-MIMO) method in which the transmitting antenna is switched by time division in order to widen the virtual aperture by using a plurality of transmitting antennas. ing.
  • TDM-MIMO radar device if the subject moves while switching the transmitting antenna, the received signal of the reflected wave of the radio wave transmitted from a different transmitting antenna cannot be coherently added. As a result, there is a problem that a three-dimensional imaging image (three-dimensional radar image) is blurred and a high-resolution image cannot be generated.
  • the present disclosure is to solve the above-mentioned problems, and an object of the present invention is to obtain a radar device and a radar image generation method capable of increasing the resolution of a three-dimensional radar image.
  • the radar device includes a transmitting unit that transmits a transmission signal to space using a plurality of transmitting antennas arranged linearly, and a plurality of receiving antennas arranged linearly in the same direction as the plurality of transmitting antennas.
  • a receiver that receives a reflected signal from which a transmission signal is reflected in space and acquires digital data of the received signal, and a transmitter that simultaneously transmits transmission signals from multiple transmitting antennas and multiple receivers.
  • a control unit that causes the transmitter and receiver to repeatedly execute a series of processes that receive reflected signals from the receiving antenna and acquire digital data, and a two-dimensional array of digital data that is sequentially acquired by the series of processes that are repeatedly executed.
  • it includes a signal processing unit that generates a three-dimensional radar image of a detection object that has moved in a direction crossing the antenna arrangement direction of the plurality of transmitting antennas and the plurality of receiving antennas.
  • a transmitting unit transmits a transmission signal to space using a plurality of transmitting antennas arranged linearly, and a plurality of receiving antennas arranged linearly in the same direction as the plurality of transmitting antennas are used.
  • the receiving unit receives the reflected signal from which the transmitted signal is reflected in space, transmits the transmitted signal from multiple transmitting antennas at the same time, receives the reflected signal from multiple receiving antennas, and acquires digital data.
  • the transmitter and receiver are made to execute repeatedly. By using the digital data sequentially acquired in a series of repeatedly executed processes as two-dimensional array data, the three-dimensional detection target moves in the direction crossing the antenna array direction of the multiple transmitting antennas and the plurality of receiving antennas. Generate a radar image.
  • the radar device can improve the resolution of the three-dimensional radar image.
  • FIG. It is a block diagram which shows the structure of the radar apparatus which concerns on Embodiment 1.
  • FIG. It is a waveform diagram which shows the modulation pattern of the transmission signal transmitted from four transmission antennas provided in a transmitter. It is a schematic diagram which shows the positional relationship between a one-dimensional array antenna and a target. It is a yz plan view which shows the image
  • FIG. 5 is an image diagram showing a cross-section of the k y '-k z' plane parallel to two three-dimensional sphere.
  • Is a schematic diagram showing the frequency region in k y '-k z' plane about the wave number k x.
  • In k y '-k z' plane about obtain the wave number k x by two transmitting antennas is an image diagram showing a frequency domain.
  • It is a flowchart which shows the radar image generation method which concerns on Embodiment 1.
  • It is a schematic diagram which shows an example of the positional relationship between a one-dimensional array antenna and a target.
  • 3D spectrum S (k x, k y, k) is a spectrum diagram showing a 2-dimensional spectrum which is projected in a two-dimensional plane of k x -k.
  • FIG. 12A is a cross-sectional view showing a cross section of an image intensity distribution calculated based on a received signal derived from the transmitted signal of the transmitting antenna Tx1 along the yz plane
  • FIG. 12B is a cross-sectional view showing the image intensity of FIG. 12A.
  • FIG. 12C is a cross-sectional view showing a cross section of the distribution cut along the xy plane
  • FIG. 12C is a cross-sectional view showing a cross section of the image intensity distribution of FIG. 12A cut on the xy plane
  • FIG. 12D is a target position. It is a distribution map showing the image intensity distribution of FIG. 12A on a straight line along the y-axis.
  • FIG. 13A is a cross-sectional view showing a cross section of an image intensity distribution calculated based on the received signals derived from the transmitted signals of the transmitting antennas Tx1 to Tx4 along the yz plane
  • FIG. 13B is a cross-sectional view of FIG. 13A.
  • 13C is a cross-sectional view showing a cross section of the image intensity distribution cut at the xy plane
  • FIG. 13C is a cross-sectional view showing a cross section of the image intensity distribution of FIG. 13A cut at the xy plane
  • FIG. 13D is a cross-sectional view.
  • FIG. 14A is a block diagram showing a hardware configuration that realizes the function of the radar signal processor included in the radar device according to the first embodiment
  • FIG. 14B is a block diagram showing the radar signal processing included in the radar device according to the first embodiment.
  • It is a block diagram which shows the hardware configuration which executes the software which realizes the function of a device.
  • FIG. 18A is a characteristic diagram showing the characteristics of the signal related to the point diffusion function of the point target when the radar image of the point target is reproduced using the one-dimensional array antenna in which all the transmitting antennas are arranged in the center of the substrate.
  • FIG. 18B is a characteristic diagram showing the characteristics of the signal related to the point diffusion function of the point target at which the radar image is reproduced by using the one-dimensional array antenna in which at least one transmitting antenna is arranged at both ends of the substrate.
  • FMCW frequency-modulated continuous wave
  • SFW step frequency continuous wave
  • FIG. 1 is a block diagram showing a configuration of a radar device 1 according to a first embodiment, and shows an FMCW type radar device.
  • the radar device 1 uses a one-dimensional array antenna to transmit a transmission signal as a radio wave in space, and the transmitted radio wave is detected by receiving a reflected wave that is reflected by a detection object existing in the space and returned. Generates a 3D radar image of the object.
  • the radar device 1 includes a radar signal processor 2, a transmission signal generator 3, a transmitter 4, a receiver 5, a beat signal generator 6, and an A / D converter 7.
  • the transmitter 4 includes M transmitting antennas Tx1 to TxM, and the receiver 5 includes N receiving antennas Rx1 to RxN.
  • the radar signal processor 2 controls the transmission signal generator 3 to transmit a transmission signal, and generates a three-dimensional radar image by processing the received signal of the reflected signal reflected by the transmission signal.
  • the transmission signal generator 3 is controlled by the radar signal processor 2 to generate a transmission signal, and outputs the transmission signal to the transmitter 4.
  • the transmitter 4 transmits the transmission signal output from the transmission signal generator 3 to the space by using M transmission antennas Tx1 to TxM.
  • the transmission signal generator 3 and the transmitter 4 function as a transmitter that transmits a transmission signal to space.
  • the receiver 5 uses N receiving antennas Rx1 to RxN to receive the reflected signal that is reflected in space and returned to the radar device 1.
  • the beat signal generator 6 generates a beat signal, which is an intermediate frequency signal (IF signal), using the received signal received by the receiver 5 and the transmitted signal transmitted by the transmitter 4, and filters the generated beat signal.
  • the processed beat signal is output to the A / D converter 7.
  • the A / D converter 7 is controlled by the radar signal processor 2 to perform analog / digital conversion of the beat signal output from the beat signal generator 6. For example, the A / D converter 7 generates digital data for each receiving antenna by sampling the beat signals corresponding to each of the receiving antennas Rx1 to RxN at a preset sampling frequency and sampling point.
  • the receiver 5, the beat signal generator 6, and the A / D converter 7 function as a receiver for receiving the reflected signal in which the transmitted signal is reflected in space.
  • the M transmitting antennas Tx1 to TxM and the N receiving antennas Rx1 to RxN constitute an array antenna (hereinafter referred to as a one-dimensional array antenna).
  • the M transmitting antennas Tx1 to TxM are arranged linearly, and the N receiving antennas Rx1 to RxN are linearly arranged in the same direction as the transmitting antenna.
  • the radar signal processor 2 includes a control unit 21, a data storage unit 22, and a signal processing unit 23.
  • the control unit 21 simultaneously transmits transmission signals from M transmission antennas Tx1 to TxM, and the reception antennas Rx1 to RxN receive the reflected signals reflected by the transmission signals, and a series of processes for acquiring digital data of the reception signals. Is repeatedly executed by the transmitting unit and the receiving unit.
  • the control unit 21 outputs the VCO 31, the phase control unit 32, the power amplifier 33, and the A / D conversion by outputting the VCO control signal, the phase control signal, the transmission control signal, the A / D control signal, and the signal processing control signal. Control the vessel 7.
  • the data storage unit 22 stores digital data of the voltage value of each beat signal A / D converted by the A / D converter 7.
  • the data storage unit 22 may be a component provided in a storage device provided separately from the radar signal processor 2.
  • the radar signal processor 2 includes a control unit 21 and a signal processing unit 23, and the signal processing unit 23 reads digital data of the voltage value of each beat signal from the data storage unit 22 provided in the storage device.
  • the signal processing unit 23 uses the digital data sequentially acquired for the detection target object that has moved in the direction crossing the antenna arrangement direction of the one-dimensional array antenna according to the signal processing control signal as the two-dimensional array data, whereby the detection target object is used. Generate a 3D radar image.
  • the transmission signal generator 3 includes a voltage controlled oscillator (VCO) 31, a phase control unit 32, and a power amplifier 33.
  • VCO 31 generates a VCO control signal output from the control unit 21, that is, a transmission signal modulated so that the frequency changes with time according to the control voltage of the VCO 31. Further, the VCO 31 outputs the generated transmission signal to the phase control unit 32 and the beat signal generator 6.
  • the phase control unit 32 sets an arbitrary initial phase for each transmission channel with respect to the transmission signal output from the VCO 31 according to the phase control signal output from the control unit 21, and the transmission channel in which the initial phase is set. Each transmission signal is output to the power amplifier 33.
  • the power amplifier 33 amplifies the transmission signal output from the phase control unit 32 for each transmission channel to a preset intensity according to the transmission control signal output from the control unit 21, and transmits the amplified transmission signal to the transmitter 4. Output.
  • the transmitter 4 includes M transmitting antennas Tx1 to TxM, and uses the transmitting antennas Tx1 to TxM to transmit the amplified transmission signal input from the power amplifier 33 to space as an electromagnetic wave.
  • the receiver 5 includes N receiving antennas Rx1 to RxN, and uses the receiving antennas Rx1 to RxN to receive the electromagnetic waves transmitted by the transmitter 4 reflected on the space side and returned. Each receiving antenna outputs the received electromagnetic wave as a receiving signal to the beat signal generator 6.
  • the beat signal generator 6 includes a low noise amplifier (LNA) 61, a distribution circuit 62, a mixer 63, and a filtering circuit 64.
  • the LNA 61 amplifies the reception signal for each reception channel received by the reception antennas Rx1 to RxN.
  • the distribution circuit 62 distributes the transmission signal output from the VCO 31 to the mixer 63 for each reception channel.
  • the mixer 63 generates a beat signal for each reception channel by using the reception signal amplified for each reception channel and the transmission signal distributed for each reception channel by the distribution circuit 62, and outputs the beat signal to the filtering circuit 64.
  • the filtering circuit 64 includes a bandpass filter (BPF) and an amplifier.
  • BPF bandpass filter
  • the filtering circuit 64 filters the beat signal output from the mixer 63 for each reception channel, and outputs the filtered beat signal for each reception channel to the A / D converter 7.
  • Each beat signal output to the A / D converter 7 suppresses low frequency components and high frequency components unnecessary for the detection process by the radar device 1, and is amplified to a preset intensity by the amplifier.
  • the digital data of the beat signal for each reception channel A / D converted by the A / D converter 7 is stored in the data storage unit 22.
  • the control unit 21 When the digital data of the received signal corresponding to the transmission signal is stored in the data storage unit 22, the control unit 21 outputs the signal processing control signal to the signal processing unit 23.
  • the signal processing unit 23 generates a three-dimensional image of a detection target (hereinafter referred to as a target) using the digital data stored in the data storage unit 22 according to the signal processing control signal from the control unit 21. ..
  • FIG. 2 is a waveform diagram showing a modulation pattern of transmission signals transmitted from the four transmission antennas Tx1 to Tx4 included in the transmitter 4.
  • the horizontal axis is time and the vertical axis is frequency.
  • the modulation pattern of the transmission signal includes, for example, FMCW modulation using time-continuous up-chirp modulation.
  • the FMCW-modulated transmission signal is appropriately described as an FMCW signal.
  • the FMCW signals (chirp signals) transmitted by the four transmitting antennas Tx1 to Tx4 are simultaneously transmitted at regular intervals.
  • the amount of phase rotation between the chirp signals of the initial phase of the transmitted chirp signal is set to be different for each transmitting antenna by the phase control unit 32.
  • the phase rotation amount between the chirp signals in the transmitting antenna Tx1 is 0 (rad)
  • the phase rotation amount between the chirp signals in the transmitting antenna Tx2 is ⁇ / 2 (rad).
  • the phase rotation amount between the chirp signals in the transmitting antenna Tx3 is ⁇ (rad)
  • the phase rotation amount between the chirp signals in the transmitting antenna Tx4 is ⁇ / 2 (rad).
  • FIG. 3 is a schematic diagram showing the positional relationship between the one-dimensional array antenna 8 and the target 9.
  • the substrate 8a is, for example, a rectangular substrate as shown in FIG. 3, and the longitudinal direction is arranged along the y direction in the xyz coordinate system of FIG.
  • the transmitting antennas Tx1 to Tx4 are each fed via the substrate 8a, and the receiving antennas Rx1 to Rx16 output the received signal to the receiver 5 via the substrate 8a.
  • the transmitting antennas Tx1 to Tx4 constitute a transmitting linear array 41 linearly arranged along the y direction on one surface of the substrate 8a.
  • the position of the connection point where the transmitting antenna Tx2 is connected to the substrate 8a is the coordinates ( ⁇ x t2 , y 0 + ⁇ y t2 , z 0 ).
  • the position of the connection point where the transmitting antenna Tx3 is connected to the substrate 8a is the coordinates ( ⁇ x t3 , y 0 + ⁇ y t3 , z 0 ).
  • the position of the connection point where the transmitting antenna Tx4 is connected to the substrate 8a is the coordinates ( ⁇ x t4 , y 0 + ⁇ y t4 , z 0 ).
  • the 16 receiving antennas Rx1 to Rx16 constitute a receiving linear array 51 linearly arranged along the y direction on one of the above surfaces of the substrate 8a.
  • the position of the center MP of the reception linear array 51 is the coordinates (0, y 0 , z 0 ).
  • the four transmitting antennas Tx1 to Tx4 irradiate electromagnetic waves until a predetermined fixed irradiation time elapses. Further, it is assumed that the target 9 moves linearly at a constant velocity in the ⁇ x direction and passes through the coordinate origin O at the intermediate time of the irradiation time.
  • the transmission signal s Tx1 (t) transmitted from the transmission antenna Tx1 can be expressed by the following equation (1).
  • f 0 is the sweep start frequency
  • is the modulation slope
  • ⁇ 0 is the initial phase
  • T is the modulation time, and the relationship of 0 ⁇ t ⁇ T is established.
  • the received signal s Rx (t) is represented by the following equation (2).
  • c is the speed of light.
  • the mixer 63 mixes the detected reception signal s Rx (t) and the transmission signal s Tx1 (t) to generate a beat signal s IF (t) represented by the following equation (3).
  • LPF [ ⁇ ] represents the low-pass filter processing in the filtering circuit 64
  • j is an imaginary unit.
  • s Tx1 + 90 deg (t) is a signal obtained by advancing the phase of the transmission signal s Tx1 (t) by 90 (deg).
  • IQ detection is assumed, it is possible to generate a similar s IF (t) by Hilbert transform even in non-IQ detection.
  • a pseudo-time signal s quasi ( ⁇ ) represented by the following equation (4), such as a received signal of a pulse radar, is obtained.
  • the phase of the beat signal s IF (t) is rotated by multiplying it by exp ⁇ j2 ⁇ (2R ⁇ / c) (T / 2) ⁇ .
  • the sampling frequency of the received signal s Rx (t) is F s
  • the range of the pseudo time ⁇ is ⁇ F s / 2 ⁇ ⁇ ⁇ F s / 2.
  • the pseudo-time signal s quasi ( ⁇ ) has an extra phase term exp ⁇ -j2 ⁇ (2R 2 ⁇ / c 2 ) ⁇ for the distance R, so the following equation (5) By canceling this phase term by the phase correction process using), a phase-corrected pseudo-time signal s quasi-comp ( ⁇ ) is generated.
  • the spectrum S (f) converted into the frequency domain is calculated as shown in the following equation (6). ..
  • the range of the frequency f is ⁇ N / 2F s ⁇ f ⁇ N / 2F s .
  • the target is a point target, but in reality, a spatially expansive target is the target of imaging. Therefore, the target is represented by a function g (x, y, z) representing the reflection intensity at each coordinate (x, y, z) in the three-dimensional space shown in FIG. It is also necessary to correctly reflect each coordinate in the geometry shown in FIG. 3 for the round-trip distance 2R of the electromagnetic wave.
  • the spectrum S (x', y', k) of the reflected signal in which the transmitted signal transmitted by the transmitting antenna Tx1 is reflected by the target g (x, y, z) is represented by the following equation (7).
  • x' is the x coordinate of the center position of the target when the transmission signal is transmitted
  • y' corresponds to the value of y coordinate ⁇ y 0 of the receiving antenna that received the reflected signal of the transmission signal.
  • l (x', y') is represented by the following formula (8)
  • l T (x') is represented by the following formula (9)
  • l R (x', y') is represented by the following formula (x', y'). It is represented by the formula (10).
  • k y ' is defined by the following equation (12)
  • k z' is defined by the following equation (13)
  • k z is defined by the following formula (14).
  • ky and inc are defined by the following equation (15)
  • k z and inc are defined by the following equation (16).
  • k yz is defined by the following equation (17).
  • cos ⁇ inc is defined by the following equation (18)
  • sin ⁇ inc is defined by the following equation (19).
  • ⁇ in the above equation (11) is defined by the following equation (20).
  • R 1 is defined by the following equation (21)
  • R 2 is defined by the following equation (22).
  • Frequency signal S (k x, k y, k) with respect to, performing phase correction of the bulk compression function and bistatic arrangement correction term described above. Thereafter, the coordinate system by Stolt interpolation, by performing (k x, k y ', k z') 3 -D inverse FFT after resampling at equal intervals grid in the coordinate, the target g (x, y, z) is reproduced.
  • FIG. 4 is a yz plan view showing an image target area.
  • (K x, k y, k z) and (k x, k y ', k z') relational expression is approximate coordinate center (y cnt in the y-z plane when the approximate plane wave transmission path, z cnt ) Depends on. In this case, as shown in FIG.
  • the imaged area is divided into two-dimensional cells that can be approximated by a plane wave in the yz plane, and image reproduction is performed for each xyz three-dimensional area corresponding to each cell. Need to do (mosaic polar format).
  • ⁇ ino corresponds to an angle seen from the + y direction to the center of each cell.
  • the above equations (20) to (22) are bistatic placement correction terms, and also in these equations, the value of the center (y ct , z ct ) of each cell shown in FIG. 4 is (y, z). It is used as a representative value of.
  • This correction term is a term for correcting the extra phase rotation caused by the position offset in the x-axis direction between the transmission linear array 41 and the reception linear array 51.
  • FIG. 5 shows a three-dimensional sphere S1 centered on the origin O with an upper limit value k max of the wave number k determined by the transmission frequency band in the three-dimensional wave number space as a radius, and a radius centered on the origin O with a lower limit value kmin as a radius.
  • W is the bandwidth of the wave number k.
  • Figure 6 is an image diagram showing the wavenumber k through the point C corresponding to the x, k y '-k z' plane parallel to two three-dimensional sphere of the cross section of FIG.
  • Section S1a passes the point corresponding to the wave number k x in the support range C (k x / 2,0,0), a cross section of the k y '-k z' plane parallel to the three-dimensional sphere S1.
  • Section S2a passes the point C (k x / 2,0,0), a cross section of the k y '-k z' plane parallel to the three-dimensional sphere S2.
  • the radius of the circle D1 is the distance kyz, min between the point C and the intersection A, and the radius of the circle D2 is the distance kyz, max between the point C and the intersection B.
  • k x / 2 axially in the region E which fanned the frequency at k y '-k z' plane about wavenumber k x Area support.
  • Figure 7 is an image diagram showing a frequency region in the k y '-k z' plane about the wave number k x.
  • the support of the transmitted signal transmitted by each transmitting antenna may have a different shift amount on the Doppler frequency axis.
  • different initial phase changes may be set for each chirp for the transmission signals transmitted by different transmission antennas.
  • the target g (x, y, z) is restored using the reflected signal S (x', y', k).
  • Table in the above formulas (11) to the equation (22) by changing the part relating to the wave number k x the wave number k x ', target g (x, y, z) is the following formula (24) Will be done.
  • k x ' is defined by the following equation (25)
  • k y' is defined by the above formula (12)
  • Kz is defined by the above equation (14).
  • ky and inc are defined by the above equation (15)
  • k z and inc are defined by the above equation (16).
  • k yz is defined by the following equation (26).
  • cos ⁇ inc is defined by the above equation (18)
  • sin ⁇ inc is defined by the above equation (19).
  • ⁇ in the above formula (24) is represented by the following formula (27).
  • R 1 is defined by the above equation (21), and R 2 is defined by the above equation (22).
  • Stolt [ ⁇ ] function is a function indicating the Stolt interpolation in the (k x, k y, k z) ⁇ (k x ', k y', k z ') There is.
  • the equations (24) to (27) above mean that the three-dimensional support shifts by the initial phase change amount a (rad / m) on the Doppler frequency axis.
  • Figure 8 is an image diagram showing a frequency region in the k y '-k z' plane about the wave number k x obtained by the transmission antenna Tx1 and the transmission antenna Tx2.
  • the transmitting antenna Tx1 and the transmitting antenna Tx2 did not change the initial phase between the chirps, and the transmitting antenna Tx2 changed the initial phase between the chirps by the initial phase change amount a (rad / m).
  • the frequency support is the fan-shaped regions E1 and E2 shown in FIG. 8 in the k x / 2-k z'two-dimensional plane.
  • the initial phase set for the signal transmitted by the transmission antenna Tx2 is rotated by ⁇ (rad) between chirps, the speed of movement along the x direction of the target is v (m / s), and the chirp interval is set.
  • T (s) the initial phase change amount a (rad / m) is represented by the following equation (29).
  • the initial phase of is changed as shown in FIG. As a result, 4-transmission DDM-MIMO is realized.
  • FIG. 9 is a flowchart showing the radar image generation method according to the first embodiment, showing the generation of a target three-dimensional radar image using the 4-transmission DDM-MIMO processing in the geometry of FIG.
  • the control unit 21 simultaneously transmits a transmission signal from the four transmission antennas Tx1 to Tx4, the reception antenna Rx receives the reflected signal reflected by the transmission signal, and the digital data of the reception signal is received.
  • the transmitter and receiver are controlled so as to repeat a series of processes to acquire.
  • the signal processing unit 23 regards the digital data sequentially acquired for the target 9 that has moved across the antenna array direction of the one-dimensional array antenna 8 as the two-dimensional array data, and performs the inverse composite aperture processing. That is, a series of processes from step ST1 to step ST9 is an inverse synthetic aperture process, and this inverse synthetic aperture process generates a three-dimensional radar image of the target.
  • the control unit 21 simultaneously transmits a transmission signal having a different initial phase change amount a for each transmission antenna from the transmission antennas Tx1 to Tx4, thereby multiplexing the transmission signal on the Doppler frequency axis.
  • the signal processing unit 23 generates a target three-dimensional radar image for each transmitting antenna based on the component of the reflected signal corresponding to the transmitting signal for each transmitting antenna, and adds the three-dimensional radar image for each transmitting antenna. , Generate the final 3D radar image of the target.
  • the signal processing unit 23 performs a three-dimensional FFT on the time signal received by each receiving antenna for the reflected signal in which the transmission signal is reflected at each movement position of the target moving along the x direction.
  • Data S (x', y', k) is generated (step ST1).
  • the signal processing unit 23 performs S (k x , ky, k) by performing a two-dimensional Fourier transform on the three-dimensional data S (x', y', k) in the x'direction and the y'direction.
  • Generate step ST2).
  • Signal four transmission antennas Tx1 ⁇ Tx4 is transmitted, because they are separated on k x axis, for the signal the transmission antennas Tx1 ⁇ Tx4 transmits, processes independently of one another from step ST3 ⁇ step ST8 following It is done.
  • the signal processing unit 23 among the respective y-z space of a plurality of two-dimensional cells, the selected y-z space as a processing target, S '(k x, k y, k) with respect to, the transmission path by performing bulk compression, S '' (k x, k y, k
  • y cnt, z cnt) the three-dimensional inverse Fourier transform on (step ST7).
  • y ct , z ct ) in the yz space of the two-dimensional cell selected as the processing target is generated.
  • y ct , z ct ) is appropriately described as a reproduced image g (x, y, z
  • the signal processing unit 23 combines a plurality of reproduced images g (x, y, z
  • step ST3 to step ST8 the processing for each transmitting antenna (step ST3 to step ST8) can be performed independently of each other. Further, the processing for each yz small cell space of the two-dimensional cell (steps ST4 to ST7) can also be processed independently of each other. These points can be processed at high speed by parallel processing using a multi-core CPU or GPU.
  • FIG. 10 is a schematic diagram showing an example of the positional relationship between the one-dimensional array antenna 8 and the target 9.
  • the image reproduction result of the point target by the transmission antenna Tx1 and the image reproduction result by the 4-transmission DDM-MIMO using the transmission antennas Tx1 to Tx4 are compared.
  • the movement of the target 9 is the same as that in FIG. 3, and the distance from the antenna surface of the substrate 8a to the movement locus of the target 9 is 2 (m).
  • the four transmitting antennas Tx1 to Tx4 form a linear array in which they are arranged at equal intervals, and the distance between the transmitting antennas is 6 (mm).
  • the 16 receiving antennas Rx1, Rx2, ..., Rx16 also form a linear array in which they are arranged at equal intervals, and the intervals between the receiving antennas are 8 (mm).
  • the one-dimensional array antenna 8 has a MIMO configuration of 4 transmissions and 16 receptions, and is regarded as a virtual array in which the element antennas are spaced by 2 (mm) and the total aperture is about 12 (cm) when viewed from a sufficient distance.
  • the modulation method is linear FMCW modulation, and the frequency band is 77 to 81 (GHz). Further, the horizontal directivity width of the transmission / reception antenna is assumed to be ⁇ 13 (deg), and signals from a wider angle range are not received.
  • the initial phase change sequence for each chirp set for each transmitting antenna is the same as in FIG.
  • Figure 11 is a three-dimensional spectrum S (k x, k y, k) is a spectrum diagram showing a 2-dimensional spectrum which is projected in a two-dimensional plane of k x -k, lower part of FIG. 11, upper The figure shows a spectrum portion to which the reference numeral F is added.
  • signals transmitted from the transmission antennas Tx1 ⁇ Tx4 is separated for each transmission antenna in the k x axis, the radar device 1, the fourth transmission DDM-MIMO is realized. Further, as shown in the lower part of FIG. 11, a two-dimensional spectrum corresponding to each transmitting antenna is generated.
  • FIG. 12A is a cross-sectional view showing a cross section of an image intensity distribution calculated based on a received signal derived from a transmitted signal of the transmitting antenna Tx1 cut along an yz plane.
  • the image intensity distribution shown in FIG. 12A is a reproduction result of a three-dimensional radar image.
  • FIG. 12B is a diagram showing a cross section of the image intensity distribution of FIG. 12A cut along the xz plane.
  • FIG. 12C is a cross-sectional view showing a cross section of the image intensity distribution of FIG. 12A cut along the xy plane.
  • FIG. 12D is a distribution map showing the image intensity distribution of FIG. 12A on a straight line passing through the target position and along the y-axis.
  • FIG. 13A is a cross-sectional view showing a cross section of an image intensity distribution calculated based on the received signals derived from the transmitted signals of the transmitting antennas Tx1 to Tx4, cut along the yz plane.
  • FIG. 13B is a diagram showing a cross section of the image intensity distribution of FIG. 13A cut along the xz plane.
  • FIG. 13C is a cross-sectional view showing a cross section of the image intensity distribution of FIG. 13A cut along the xy plane.
  • FIG. 13D is a distribution map showing the image intensity distribution of FIG. 13A on a straight line passing through the target position and along the y-axis.
  • the integrated sidelobe ratio (ISLR) has also improved by about 4 dB.
  • the resolution of the three-dimensional radar image generated by the radar device 1 does not change.
  • the resolution of the three-dimensional radar image is improved by arranging the transmitting antenna and the receiving antenna so that the spacing between the element antennas (MIMO virtual array spacing) is increased and the total aperture is widened. It is possible.
  • FIG. 14A is a block diagram showing a hardware configuration that realizes the functions of the radar signal processor 2.
  • FIG. 14B is a block diagram showing a hardware configuration for executing software that realizes the functions of the radar signal processor 2.
  • the input / output interface 100 inputs A / D converted digital data by the A / D converter 7, relays the control signal from the control unit 21, and outputs the control signal to the transmission unit and the reception unit.
  • the storage device 101 is a storage device that stores data for which the input / output interface 100 has received an input, and is a data storage unit 22.
  • the radar signal processor 2 includes a processing circuit that executes the processes from step ST1 to step ST9 shown in FIG.
  • the processing circuit may be dedicated hardware, or may be a CPU that executes a program stored in the memory.
  • the processing circuit 102 may be, for example, a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, or an ASIC (Application Specific Integrated Circuit). ), FPGA (Field-Programmable Gate Array), or a combination of these.
  • the functions of the control unit 21 and the signal processing unit 23 included in the radar signal processor 2 may be realized by separate processing circuits, or these functions may be collectively realized by one processing circuit.
  • the processing circuit is the processor 103 shown in FIG. 14B
  • the functions of the control unit 21 and the signal processing unit 23 included in the radar signal processor 2 are realized by software, firmware, or a combination of software and firmware.
  • the software or firmware is described as a program and stored in the memory 104.
  • the processor 103 realizes the functions of the control unit 21 and the signal processing unit 23 by reading and executing the program stored in the memory 104.
  • the radar signal processor 2 includes a memory 104 that stores a program in which the processes of steps ST1 to ST9 shown in FIG. 9 are executed as a result when the radar signal processor 2 is executed by the processor 103. These programs cause the computer to execute the procedure or method of the control unit 21 and the signal processing unit 23.
  • the memory 104 may be a computer-readable storage medium in which a program for causing the computer to function as the control unit 21 and the signal processing unit 23 is stored.
  • the memory 104 may be, for example, a RAM (Random Access Memory), a ROM (Read Only Memory), a flash memory, an EPROM (Erasable Programmable Read Only Memory), an EEPROM (Electrically-volatile), an EEPROM (Electricularly-volatile), or the like. This includes disks, flexible disks, optical disks, compact disks, mini disks, DVDs, and the like.
  • a part of the functions of the control unit 21 and the signal processing unit 23 may be realized by dedicated hardware, and the remaining part may be realized by software or firmware.
  • the function of the control unit 21 is realized by the processing circuit 102, which is dedicated hardware, and the function of the signal processing unit 23 is realized by the processor 103 reading and executing the program stored in the memory 104. ..
  • the processing circuit can realize the above-mentioned functions by hardware, software, firmware or a combination thereof.
  • the round-trip distance of the transmission signal is l (x', y'), and is expressed by the above equation (8).
  • the outward distance of the transmission signal is l T (x'), which is expressed by the above equation (9), and the return distance of the reflected signal of the transmission signal is l R (x', y'), which is the above equation (10). ).
  • the signal processing unit 23, S (x ', y' , k) a, x 'and y' frequency signals is calculated by two-dimensional Fourier transform with respect to S (k x, k y, k) and g (x, The relationship of the following equation (31) holds with y, z).
  • the final integration part on the right side of the above equation (31) can be calculated by the stationary phase method (MOSP) and is represented by the following equation (32).
  • a plane wave approximation is introduced to remove this square root for -jk yz ⁇ (y-y 0- ⁇ y t ) 2 + (z-z 0 ) 2 ⁇ 1/2 on the right side of the above equation (48).
  • the square root portion can be approximated by a plane wave as shown by the following equation (49).
  • ky and inc are defined by the above equation (15)
  • k z and inc are defined by the above equation (16)
  • cos ⁇ inc is defined by the above equation (18)
  • sin ⁇ inc is defined by the above equation (19).
  • the radar device 1 transmits a transmission signal to the space using the transmission antennas Tx1 to TxM arranged linearly, and is arranged linearly in the same direction as the transmission antenna.
  • the receiving antennas Rx1 to RxN are used to receive the reflected signal reflected in space, the transmitting antennas Tx1 to TxM simultaneously transmit the transmitting signal, and the receiving antennas Rx1 to RxN receive the reflected signal to transmit digital data. Repeatedly execute a series of acquisition processes.
  • the target 9 3 has moved in a direction crossing the antenna arrangement direction of the transmitting antennas Tx1 to TxM and the receiving antennas Rx1 to RxN. Generate a dimensional radar image.
  • the signal transmission intervals of the transmission antennas Tx1 to TxM are shortened, so that the influence of the movement of the target 9 within the time of the signal transmission interval is reduced.
  • the radar device 1 of the DDM-MIMO system has one transmitting antenna.
  • the signal transmission interval of can be shortened to T / N (s). Therefore, the received signals derived from the transmitted signals from different transmitting antennas can be coherently added, and the SN ratio of the 3D radar image becomes N times, so that the radar device 1 enhances the resolution of the 3D radar image. be able to.
  • FIG. 15 is a block diagram showing the configuration of the radar device 1A according to the second embodiment, and shows the FMCW type radar device.
  • the same components as those in FIG. 1 are designated by the same reference numerals, and the description thereof will be omitted.
  • the radar device 1A transmits a transmission signal as a radio wave to the space using a one-dimensional array antenna, and receives the reflected wave returned after the transmitted radio wave is reflected by a target existing in the space. By doing so, a three-dimensional radar image of the target is generated.
  • the radar device 1A is used in a millimeter-wave security system for body-checking a detection target that continuously arrives in real time without hindering its movement.
  • the radar device 1A will be described as a millimeter-wave security system that reproduces three-dimensional radar images of moving detection objects one after another.
  • the radar device 1A includes a radar signal processor 2A, a transmission signal generator 3, a transmitter 4, a receiver 5, a beat signal generator 6, and an A / D converter 7.
  • the radar signal processor 2A includes a control unit 21, a data storage unit 22, and a signal processing unit 23A.
  • the control unit 21 sets the pulse repetition interval of the transmission signals simultaneously transmitted from the transmission antennas Tx1 to TxM to be constant, and sets the initial phase change amount between the transmission signals of the pulse repetition intervals to a value different for each transmission antenna.
  • the signal processing unit 23A sequentially Fourier transforms the digital data of the reflected signal sequentially received by the receiving antennas Rx1 to RxN at pulse repetition intervals into a signal in the Doppler frequency region, thereby performing a three-dimensional radar image of the object to be detected. To generate.
  • FIG. 16 is a flowchart showing a radar image generation method according to the second embodiment.
  • the processing from step ST4a to step ST10a is the same as the processing from step ST3 to step ST9 shown in FIG.
  • the signal processing unit 23A generates received signals s (y', k
  • the radar device 1A can increase the calculation speed required for generating the three-dimensional radar image. For example, when reproducing a 3D radar image that is continuous in time, the chirp number of the latest transmission signal is set to h, and the chirp numbers (h-N + 1) to h are received signals for N chapters. Reproduce the 3D radar image.
  • the signal processing unit 23A performs a one-dimensional Fourier transform on the received signal s (y', k
  • the signal processing unit 23A receives signals for the latest N chirps from the chirp numbers (h—N + 1) to h. Perform a three-dimensional Fourier transform. Thus, three-dimensional frequency signal S (k x, k y, k
  • the signal processing unit 23A is represented by the following formula (54), the following equation (55) and the following equation (56) all k y, by applying against k, executes the successive FFT (step ST3a).
  • each of the S transmit antennas Tx1 ⁇ Tx4 (k x, k y, k) is generated.
  • the following processing is the same as in step ST3 and subsequent steps of FIG.
  • the signal processing unit 23A sequentially transfers the digital data of the reflected signal sequentially received by the receiving antenna Rx at the pulse repetition interval to the signal in the Doppler frequency domain.
  • a three-dimensional radar image of the target 9 is generated by Fourier transform.
  • the radar device 1A can generate a three-dimensional radar image of the target 9 at high speed.
  • FIG. 17 is a schematic diagram showing the positional relationship between the three one-dimensional array antennas 8 (1) to 8 (3) included in the radar device according to the third embodiment and the target 9.
  • the virtual opening of MIMO can be realized by arranging a plurality of substrates 8a shown in FIG. 3 closely side by side in the y direction.
  • the boards 8a (1) to 8a (3) are arranged in the y direction.
  • large gaps in which virtual elements do not exist occur at a plurality of locations.
  • all the received signals received by the one-dimensional array antennas 8 (1) to (3) provided on each of the three boards 8a (1) to 8a (3) are coherently added. As a result, a large ripple occurs in the reproduced 3D radar image.
  • the first end portion (end portion in the + y direction) of each substrate 8a and the second end portion facing the first end portion are opposed to each other.
  • At least one transmitting antenna is arranged at the end of (the end in the ⁇ y direction).
  • a MIMO virtual array composed of one-dimensional array antennas 8 (1) to 8 (3) has three substrates 8a in which at least one transmitting antenna is arranged at a first end and a second end. , Arranged in the y direction.
  • FIG. 18A shows a signal related to the point diffusion function of the target 9 when the three-dimensional radar image of the target 9 is reproduced by using the one-dimensional array antenna 8 in which all the transmitting antennas Tx1 to TxM are arranged in the center of the substrate 8a.
  • the characteristic H1 is the target 9 when the three-dimensional radar image of the target 9 is reproduced by the radar device 1 having one one-dimensional array antenna 8, that is, one substrate 8a shown in FIG.
  • the characteristic H2 is a signal related to the point diffusion function of the target 9 when the three-dimensional radar image of the target 9 is reproduced by the radar device having the one-dimensional array antennas 8 (1) to 8 (3) shown in FIG. It is a characteristic of.
  • FIG. 18B is a characteristic diagram showing the characteristics of the signal related to the point diffusion function of the target 9 that reproduces the three-dimensional radar image by using the one-dimensional array antenna 8 in which at least one transmitting antenna Tx is arranged at both ends of the substrate 8a.
  • the characteristic I1 is a characteristic of a signal related to the point diffusion function of the target 9 when the three-dimensional radar image of the target 9 is reproduced by the radar device 1 having one substrate 8a shown in FIG. .
  • the characteristic I2 is a signal related to the point diffusion function of the target 9 when the three-dimensional radar image of the target 9 is reproduced by the radar device having the one-dimensional array antennas 8 (1) to 8 (3) shown in FIG. It is a characteristic of.
  • the one-dimensional array antenna 8 (1) is arranged as shown in FIG. 18B.
  • the ripple of the main lobe is suppressed and the coherent peak is obtained.
  • the radar device includes three one-dimensional array antennas 8 (1) to 8 (3), and the one-dimensional array antennas 8 (1) to 8 (3) are linear. It is arranged in a shape.
  • the one-dimensional array antennas 8 (1) to 8 (3) at least one of the transmitting antennas Tx has the transmitting antennas at the first end and the second end of the substrates 8a (1) to 8a (3). Is located in.
  • the radar device shown in the first to third embodiments can be used in a system for detecting a person on an escalator or a moving walkway. For example, by placing one one-dimensional array antenna 8 on each side of the escalator or moving walkway, the radar device can determine whether a person on the escalator or moving walkway has a metal or non-metal forbidden item. Can be monitored from both the left and right sides.
  • the radar devices shown in the first to third embodiments can be used not only for people but also for baggage inspection.
  • a baggage inspection device at a station, an airport, or the like luggage is placed on a belt conveyor, and the belt conveyor operates in that state.
  • the belt conveyor operates in that state.
  • the radar device shown in the first to third embodiments by arranging the one-dimensional array antenna on the side of an outdoor road, a car, a bicycle or a person passing through the road can be made of metal or non-metal. It is possible to monitor whether or not you have prohibited items.
  • the radar device according to the present disclosure can be used, for example, as an imaging radar device.

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

レーダ装置(1)は、直線状に配置された送信アンテナ(Tx1~TxM)を用いて空間に送信信号を送信し、送信アンテナ(Tx1~TxM)と同一方向に直線状に配置された受信アンテナ(Rx1~RxN)を用いて空間で送信信号が反射された反射信号を受信し、送信アンテナ(Tx1~TxM)から同時に送信信号を送信し、受信アンテナ(Rx1~RxN)によって反射信号を受信してデジタルデータを取得する一連の処理を繰り返し実行させる制御部(21)と、繰り返し実行された一連の処理で順次取得されたデジタルデータを、2次元アレイデータとして用いることにより、送信アンテナ(Tx1~TxM)および受信アンテナ(Rx1~RxN)のアンテナ配列方向を横切る方向に移動したターゲット(9)の3次元レーダ画像を生成する信号処理部(23)とを備える。

Description

レーダ装置およびレーダ画像生成方法
 本開示は、レーダ装置およびレーダ画像生成方法に関する。
 ボディチェック装置に用いられるレーダ装置が知られている。ボディチェック装置は、空港、鉄道駅または港といった出入りする人のセキュリティチェックが必要な場所に配置される。レーダ装置を用いたボディチェック装置は、レーダ装置から、例えば、ミリ波帯の電波を送信し、被験者で電波が反射された反射波を受信することにより、被験者を撮像する装置である。ミリ波帯の電波は、人の服を透過し、人に対する放射量が少ないので、ミリ波帯の電波を用いたボディチェック装置は、金属および非金属の禁制品の識別が可能である。
 ミリ波帯の電波を用いる従来のボディチェック装置には、被験者が静止する必要があるものと、被験者が静止する必要がないものとの2種類がある。ただし、ボディチェックの効率化の観点から、被験者が停止する必要があるものよりも、停止する必要がないボディチェック装置の要望が高まっている。
 例えば、非特許文献1には、被験者が一定の速度で移動している間に、複数の送信アンテナから時分割でミリ波を送信し、2次元受信アレイがミリ波の反射波を受信することにより、3次元ミリ波イメージングを行うレーダ装置が記載されている。
Borja Gonzalez-Valdes et al., "Millimeter wave imaging architecture for cnthe-move whole body imaging," IEEE Trans. Ant. Propag., vol.64, no 5, May 2016.
 非特許文献1に記載されるレーダ装置は、複数の送信アンテナを用いて仮想開口を広げるために、時分割で送信アンテナを切り替える時分割多重-多重入力多重出力(TDM-MIMO)方式を採用している。TDM-MIMO方式のレーダ装置では、送信アンテナを切り替える間に被験者が移動すると、異なる送信アンテナから送信された電波の反射波の受信信号をコヒーレントに加算することができない。これにより、3次元イメージング画像(3次元レーダ画像)にぼけが発生して高分解能な画像を生成できないという課題があった。
 本開示は、上記課題を解決するものであり、3次元レーダ画像の分解能を高めることができるレーダ装置およびレーダ画像生成方法を得ることを目的とする。
 本開示に係るレーダ装置は、直線状に配置された複数の送信アンテナを用いて空間に送信信号を送信する送信部と、複数の送信アンテナと同一方向に直線状に配置された複数の受信アンテナを用いて空間で送信信号が反射された反射信号を受信し、受信した信号のデジタルデータを取得する受信部と、送信部が複数の送信アンテナから同時に送信信号を送信し、受信部が複数の受信アンテナによって反射信号を受信してデジタルデータを取得する一連の処理を、送信部および受信部に繰り返し実行させる制御部と、繰り返し実行された一連の処理で順次取得されたデジタルデータを2次元アレイデータとして用いることにより、複数の送信アンテナおよび複数の受信アンテナのアンテナ配列方向を横切る方向に移動した検知対象物の3次元レーダ画像を生成する信号処理部とを備える。
 本開示によれば、直線状に配置された複数の送信アンテナを用いて送信部が空間に送信信号を送信し、複数の送信アンテナと同一方向に直線状に配置された複数の受信アンテナを用いて受信部が空間で送信信号が反射された反射信号を受信し、複数の送信アンテナから同時に送信信号を送信し、複数の受信アンテナによって反射信号を受信してデジタルデータを取得する一連の処理を、送信部および受信部に繰り返し実行させる。繰り返し実行された一連の処理で順次取得されたデジタルデータを、2次元アレイデータとして用いることにより、複数の送信アンテナおよび複数の受信アンテナのアンテナ配列方向を横切る方向に移動した検知対象物の3次元レーダ画像を生成する。これにより、送信アンテナの信号送信間隔が短くなり、信号送信間隔の時間内に検知対象物が移動することによる影響が低減されて、異なる送信アンテナからの送信信号に由来する受信信号をコヒーレントに加算できるので、本開示に係るレーダ装置は、3次元レーダ画像の分解能を高めることができる。
実施の形態1に係るレーダ装置の構成を示すブロック図である。 送信器が備える4つの送信アンテナから送信された送信信号の変調パターンを示す波形図である。 1次元アレイアンテナとターゲットとの位置関係を示す概要図である。 画像化対象エリアを示すy-z平面図である。 3次元波数空間における、送信周波数帯域によって決定された波数kの上限値kmaxを半径とし原点Oを中心とした3次元球と、下限値kminを半径とし原点Oを中心とした3次元球を示すイメージ図である。 図5における波数kに対応した点Cを通り、k’-k’平面と平行な2つの3次元球の断面を示すイメージ図である。 波数kに関するk’-k’平面における周波数領域を示すイメージ図である。 2個の送信アンテナによって得られる波数kに関するk’-k’平面における周波数領域を示すイメージ図である。 実施の形態1に係るレーダ画像生成方法を示すフローチャートである。 1次元アレイアンテナとターゲットとの位置関係の一例を示す概要図である。 3次元スペクトラムS(k,k,k)がk-kの2次元平面に投影された2次元スペクトラムを示すスペクトラム図である。 図12Aは、送信アンテナTx1の送信信号に由来する受信信号に基づいて算出された画像強度分布を、y-z面で切った断面を示す断面図であり、図12Bは、図12Aの画像強度分布を、x-z面で切った断面を示す図であり、図12Cは、図12Aの画像強度分布を、x-y面で切った断面を示す断面図であり、図12Dは、ターゲット位置を通り、y軸に沿った直線上の図12Aの画像強度分布を示す分布図である。 図13Aは、送信アンテナTx1~Tx4の送信信号に由来する受信信号に基づいて算出された画像強度分布を、y-z面で切った断面を示す断面図であり、図13Bは、図13Aの画像強度分布を、x-z面で切った断面を示す図であり、図13Cは、図13Aの画像強度分布を、x-y面で切った断面を示す断面図であり、図13Dは、ターゲット位置を通り、y軸に沿った直線上の図13Aの画像強度分布を示す分布図である。 図14Aは、実施の形態1に係るレーダ装置が備えるレーダ信号処理器の機能を実現するハードウェア構成を示すブロック図であり、図14Bは、実施の形態1に係るレーダ装置が備えるレーダ信号処理器の機能を実現するソフトウェアを実行するハードウェア構成を示すブロック図である。 実施の形態2に係るレーダ装置の構成を示すブロック図である。 実施の形態2に係るレーダ画像生成方法を示すフローチャートである。 実施の形態3に係るレーダ装置が備える3つの1次元アレイアンテナとターゲットとの位置関係を示す概要図である。 図18Aは、基板の中央に全ての送信アンテナが配置された1次元アレイアンテナを用いて点ターゲットのレーダ画像を再生した場合における、点ターゲットの点拡散関数に係る信号の特性を示す特性図であり、図18Bは、基板の両端に少なくとも一つの送信アンテナが配置された1次元アレイアンテナを用いてレーダ画像を再生した点ターゲットの点拡散関数に係る信号の特性を示す特性図である。
 以下の説明においては、周波数変調連続波(FMCW)レーダ装置が使用される場合を例示する。ただし、FMCWレーダ装置の替わりに、パルスレーダ装置が使用されてもよいし、ステップ周波数連続波(SFCW)レーダ装置が使用されてもよい。
実施の形態1.
 図1は、実施の形態1に係るレーダ装置1の構成を示すブロック図であって、FMCW方式のレーダ装置を示している。レーダ装置1は、1次元アレイアンテナを用いて、空間に送信信号を電波として送信し、送信した電波が空間に存在する検知対象物で反射されて戻ってきた反射波を受信することによって、検知対象物の3次元レーダ画像を生成する。レーダ装置1は、レーダ信号処理器2、送信信号生成器3、送信器4、受信器5、ビート信号生成器6およびA/D変換器7を備える。送信器4は、M個の送信アンテナTx1~TxMを備え、受信器5は、N個の受信アンテナRx1~RxNを備えている。
 レーダ信号処理器2は、送信信号生成器3を制御して送信信号を送信させ、送信信号が反射された反射信号の受信信号を信号処理することにより3次元レーダ画像を生成する。送信信号生成器3は、レーダ信号処理器2によって制御されて送信信号を生成し、送信信号を送信器4へ出力する。送信器4は、M個の送信アンテナTx1~TxMを用いて、送信信号生成器3から出力された送信信号を空間に送信する。送信信号生成器3および送信器4は、空間に送信信号を送信する送信部として機能する。
 受信器5は、N個の受信アンテナRx1~RxNを用いて、空間で送信信号が反射されてレーダ装置1側に戻ってきた反射信号を受信する。ビート信号生成器6は、受信器5が受信した受信信号と送信器4が送信した送信信号とを用いて、中間周波数信号(IF信号)であるビート信号を生成し、生成したビート信号をフィルタリング処理したビート信号をA/D変換器7に出力する。
 A/D変換器7は、レーダ信号処理器2によって制御されて、ビート信号生成器6から出力されたビート信号をアナログ/デジタル変換する。例えば、A/D変換器7は、受信アンテナRx1~RxNのそれぞれに対応するビート信号から、予め設定されたサンプリング周波数およびサンプリング点数でサンプリングを行うことにより、受信アンテナごとのデジタルデータを生成する。受信器5、ビート信号生成器6およびA/D変換器7は、空間において送信信号が反射された反射信号を受信する受信部として機能する。
 M個の送信アンテナTx1~TxMおよびN個の受信アンテナRx1~RxNは、アレイアンテナ(以下、1次元アレイアンテナと呼称する)を構成する。1次元アレイアンテナにおいて、M個の送信アンテナTx1~TxMは、直線状に配置され、N個の受信アンテナRx1~RxNは、送信アンテナと同一方向に直線状に配列されている。
 レーダ信号処理器2は、制御部21、データ記憶部22および信号処理部23を備えている。制御部21は、M個の送信アンテナTx1~TxMから同時に送信信号を送信し、送信信号が反射された反射信号を受信アンテナRx1~RxNが受信し、受信信号のデジタルデータを取得する一連の処理を、送信部および受信部に繰り返し実行させる。例えば、制御部21は、VCO制御信号、位相制御信号、送信制御信号、A/D制御信号および信号処理制御信号を出力することにより、VCO31、位相制御部32、パワーアンプ33およびA/D変換器7を制御する。
 データ記憶部22は、A/D変換器7によってA/D変換された各ビート信号の電圧値のデジタルデータを記憶する。なお、データ記憶部22は、レーダ信号処理器2とは別に設けられた記憶装置が備える構成要素であってもよい。この場合、レーダ信号処理器2は制御部21および信号処理部23を備え、信号処理部23が、上記記憶装置に設けられたデータ記憶部22から各ビート信号の電圧値のデジタルデータを読み出す。
 信号処理部23は、信号処理制御信号に従って、1次元アレイアンテナのアンテナ配列方向を横切る方向に移動した検知対象物について順次取得されたデジタルデータを2次元アレイデータとして用いることにより、検知対象物の3次元レーダ画像を生成する。
 送信信号生成器3は、電圧制御発振器(VCO)31、位相制御部32およびパワーアンプ33を備える。VCO31は、制御部21から出力されたVCO制御信号、すなわちVCO31の制御電圧に応じて、時間的に周波数が変化するように変調された送信信号を生成する。また、VCO31は、生成した送信信号を位相制御部32およびビート信号生成器6に出力する。
 位相制御部32は、制御部21から出力された位相制御信号に応じて、VCO31から出力された送信信号に対して送信チャネルごとに任意の初期位相を設定し、初期位相が設定された送信チャネルごとの送信信号をパワーアンプ33に出力する。パワーアンプ33は、制御部21から出力された送信制御信号に従って、位相制御部32から送信チャネルごとに出力された送信信号を予め設定された強度に増幅し、増幅した送信信号を送信器4に出力する。
 送信器4は、M個の送信アンテナTx1~TxMを備え、送信アンテナTx1~TxMを用いて、パワーアンプ33から入力した増幅後の送信信号を、電磁波として空間へ送信する。受信器5は、N個の受信アンテナRx1~RxNを備え、受信アンテナRx1~RxNを用いて、送信器4によって送信された電磁波が空間側で反射されて戻ってきた電磁波を受信する。各受信アンテナは、受信した電磁波を受信信号としてビート信号生成器6に出力する。
 ビート信号生成器6は、低雑音増幅器(LNA)61、分配回路62、ミキサ63およびフィルタリング回路64を備える。LNA61は、受信アンテナRx1~RxNにより受信された受信チャネルごとの受信信号を増幅する。分配回路62は、VCO31から出力された送信信号を、受信チャネルごとのミキサ63に分配する。ミキサ63は、受信チャネルごとに増幅された受信信号と分配回路62によって受信チャネルごとに分配された送信信号とを用いて、受信チャネルごとのビート信号を生成してフィルタリング回路64に出力する。
 フィルタリング回路64は、帯域通過フィルタ(BPF)およびアンプを備える。フィルタリング回路64は、ミキサ63から受信チャネルごとに出力されたビート信号をフィルタリングし、フィルタリングされた受信チャネルごとのビート信号をA/D変換器7に出力する。A/D変換器7に出力された各ビート信号は、レーダ装置1による検知処理に不要な低い周波数成分と高い周波数成分が抑圧され、アンプによって予め設定された強度に増幅される。
 A/D変換器7は、フィルタリング回路64から出力された各ビート信号の電圧値を、制御部21から出力されたA/D制御信号に従ってA/D変換し、A/D変換後の信号をレーダ信号処理器2に出力する。A/D変換器7によってA/D変換された受信チャネルごとのビート信号のデジタルデータは、データ記憶部22に記憶される。
 制御部21は、送信信号に対応する受信信号のデジタルデータがデータ記憶部22に記憶されると、信号処理制御信号を信号処理部23に出力する。信号処理部23は、制御部21からの信号処理制御信号に従って、データ記憶部22に記憶されたデジタルデータを用いて、検知対象物(以下、ターゲットと呼称する)の3次元イメージ画像を生成する。
 図2は、送信器4が備える4つの送信アンテナTx1~Tx4から送信された送信信号の変調パターンを示す波形図である。図2において、横軸は時間であり、縦軸は周波数である。送信アンテナTx1~Tx4は、M=4、すなわち、送信アンテナが4個である。送信信号の変調パターンには、例えば、時間連続的なアップチャープ変調を用いたFMCW変調がある。以下、FMCW変調された送信信号は、FMCW信号と適宜記載される。
 4個の送信アンテナTx1~Tx4が送信するFMCW信号(チャープ信号)は、一定の周期ごとに同時に送信される。送信されるチャープ信号の初期位相のチャープ信号間での位相回転量は、位相制御部32によって送信アンテナごとに異なるように設定される。
 図2において、送信アンテナTx1におけるチャープ信号間の位相回転量は0(rad)であり、送信アンテナTx2におけるチャープ信号間の位相回転量はπ/2(rad)である。さらに、送信アンテナTx3におけるチャープ信号間の位相回転量はπ(rad)であり、送信アンテナTx4におけるチャープ信号間の位相回転量は-π/2(rad)である。送信アンテナごとにチャープ信号間の位相回転量を変更することにより、レーダ装置1において、ドップラ分割多重-多重入力多重出力(DDM-MIMO)方式が実現される。
 図3は、1次元アレイアンテナ8とターゲット9との位置関係を示す概要図である。図3において、1次元アレイアンテナ8は、例えば、基板8a、4(M=4)個の送信アンテナTx1~Tx4および16(N=16)個の受信アンテナRx1~Rx16を備えている。基板8aは、例えば、図3に示すような長方形の基板であり、長手方向が、図3のxyz座標系におけるy方向に沿って配置されている。送信アンテナTx1~Tx4は、基板8aを介してそれぞれ給電され、受信アンテナRx1~Rx16は、基板8aを介して受信信号を受信器5に出力する。
 送信アンテナTx1~Tx4は、基板8aの一方の面においてy方向に沿って、直線状に配列された送信リニアアレイ41を構成している。送信アンテナTxi(i=1~4)が基板8aに接続された接続点の位置は、座標(Δxti,y+Δyti,z)である。すなわち、送信アンテナTx1が基板8aに接続された接続点の位置は、座標(Δxt1,y+Δyt1,z)である。送信アンテナTx2が基板8aに接続された接続点の位置は、座標(Δxt2,y+Δyt2,z)である。送信アンテナTx3が基板8aに接続された接続点の位置は、座標(Δxt3,y+Δyt3,z)である。送信アンテナTx4が基板8aに接続された接続点の位置は、座標(Δxt4,y+Δyt4,z)である。
 16個の受信アンテナRx1~Rx16は、基板8aの上記一方の面においてy方向に沿って直線状に配列された受信リニアアレイ51を構成している。受信リニアアレイ51の中心MPの位置は、座標(0,y,z)である。4個の送信アンテナTx1~Tx4は、予め設定された一定の照射時間が経過するまで電磁波を照射する。また、ターゲット9は、-x方向に等速直線運動し、上記照射時間の中間時刻で座標原点Oを通過するものと仮定する。
 送信アンテナTx1から送信される送信信号sTx1(t)は、下記式(1)で表すことができる。下記式(1)において、fは掃引開始周波数であり、μは変調傾きであり、φは初期位相である。Tは変調時間であり、0≦t≦Tの関係が成り立つ。

Figure JPOXMLDOC01-appb-I000001
 送信アンテナTx1によって空間に照射され、空間に存在する点ターゲットで反射されて受信アンテナに受信されるまでの電磁波の往復距離が2Rである場合に、点ターゲットで反射されて戻ってきた反射波の受信信号sRx(t)は、下記式(2)で表される。cは光速である。

Figure JPOXMLDOC01-appb-I000002
 ミキサ63は、検波した受信信号sRx(t)と送信信号sTx1(t)をミキシングすることにより、下記式(3)で表されるビート信号sIF(t)を生成する。ただし、LPF[・]は、フィルタリング回路64におけるローパスフィルタ処理を表しており、jは虚数単位である。sTx1+90deg(t)は、送信信号sTx1(t)の位相を、90(deg)進めた信号である。IQ検波を仮定したが、非IQ検波においてもヒルベルト変換により同様のsIF(t)を生成することが可能である。

Figure JPOXMLDOC01-appb-I000003
 ビート信号sIF(t)に対して距離FFT処理が行われると、パルスレーダの受信信号のような、下記式(4)に示す疑似時間信号squasi(τ)が得られる。下記式(4)において、ビート信号sIF(t)は、exp{j2π(2Rμ/c)(T/2)}が乗算されて位相が回転されている。

Figure JPOXMLDOC01-appb-I000004
 sinc関数は、sinc(x)=sin(πx)/πxで定義される。また、受信信号sRx(t)のサンプリング周波数がFである場合に、疑似時間τの範囲は、-F/2≦τ≦F/2である。疑似時間信号squasi(τ)には、上記式(4)に示すように、距離Rに関する余分な位相項exp{-j2π(2Rμ/c)}が存在するので、下記式(5)を用いた位相補正処理によって、この位相項をキャンセルすることにより、位相補正済みの疑似時間信号squasi-comp(τ)が生成される。

Figure JPOXMLDOC01-appb-I000005
 位相補正済みの疑似時間信号squasi-comp(τ)に対して距離FFT処理が行われると、下記式(6)に示すように、周波数領域に変換されたスペクトルS(f)が算出される。受信信号sRx(t)のサンプリング数がNである場合、周波数fの範囲は、-N/2F≦f≦N/2Fである。また、波数kは、k=2π(f+μf)/cで表しことができる。

Figure JPOXMLDOC01-appb-I000006
 これまでの説明では、ターゲットが点ターゲットであると仮定したが、実際には、空間的に広がりのあるターゲットがイメージングの対象である。そこで、ターゲットは、図3に示す3次元空間の各座標(x,y,z)における反射強度を表す関数g(x,y,z)で表される。電磁波の往復距離2Rについても、図3に示すジオメトリにおける各座標を正しく反映させる必要がある。
 送信アンテナTx1によって送信された送信信号がターゲットg(x,y,z)で反射された反射信号のスペクトルS(x’,y’,k)は、下記式(7)で表される。なお、x’は、送信信号が送信されたときのターゲットの中心位置のx座標であり、y’は、送信信号の反射信号を受信した受信アンテナのy座標-yの値に対応している。ただし、l(x’,y’)は、下記式(8)で表され、l(x’)は、下記式(9)で表され、l(x’,y’)は、下記式(10)で表される。

Figure JPOXMLDOC01-appb-I000007
 スペクトルS(x’,y’,k)がx’およびy’に関して2次元フーリエ変換されて得られた周波数信号がS(k,k,k)である場合、ターゲットg(x,y,z)は、下記式(11)を用いて推定される。導出に関しては後に別途説明を行う。

Figure JPOXMLDOC01-appb-I000008
 上記式(11)において、k’は、下記式(12)で定義され、k’は、下記式(13)で定義され、kは、下記式(14)で定義される。ky,incは、下記式(15)で定義され、kz,incは、下記式(16)で定義される。kyzは、下記式(17)で定義される。cosθincは、下記式(18)で定義され、sinθincは、下記式(19)で定義される。

Figure JPOXMLDOC01-appb-I000009
 上記式(11)におけるΦは、下記式(20)で定義される。ただし、Rは下記式(21)で定義され、Rは下記式(22)で定義される。

Figure JPOXMLDOC01-appb-I000010
 上記式(11)の右辺におけるexp{-j(k(Δx/2)+k+k)}は、受信経路のバルク圧縮関数である。上記式(11)の右辺におけるexp{-j(ky,inc(y+Δy)+kz,inc)}は、平面波近似した送信経路のバルク圧縮関数である。上記式(11)の右辺におけるexp{jΦ(k,k,k|ycnt,zcnt)}は、x軸上でのバイスタティック配置補正項である。Stolt[・]関数は、(k,k,k)→(k,k’,k’)へのStolt補間処理を表している。
 周波数信号S(k,k,k)に対して、前述したバルク圧縮関数およびバイスタティック配置補正項の位相補正を行う。この後、Stolt補間によって座標系を、(k,k’,k’)座標での等間隔グリッドにリサンプリングしてから3次元逆FFTを行うことにより、ターゲットg(x,y,z)が再生される。
 (k,k,k)と(k,k’,k’)との関係式は、上記式(12)から上記式(19)に示す通りである。図4は、画像化対象エリアを示すy-z平面図である。(k,k,k)と(k,k’,k’)の関係式は、送信経路を平面波近似したときのy-z平面における近似座標中心(ycnt,zcnt)に依存している。この場合、図4に示すように、y-z平面で画像化対象エリアを平面波近似が可能な細かなセルに2次元的に区切り、各セルに対応するxyzの3次元エリアごとに画像再生を行う必要がある(モザイクポーラーフォーマット)。上記式(15)、上記式(16)、上記式(18)および上記式(19)における、θinoは、+y方向から各セルの中心までを見込んだ角度に対応している。
 上記式(20)から上記式(22)は、バイスタティック配置補正項であり、これらの式においても、図4に示す各セルの中心(ycnt,zcnt)の値を(y,z)の代表値として用いている。この補正項は、送信リニアアレイ41と受信リニアアレイ51とのx軸方向の位置オフセットに起因する余分な位相回転を補正するための項である。
 また、上記式(12)から上記式(19)までによって表される(k,k,k)と(k,k’,k’)との関係式についての幾何学的なイメージを説明する。図5は、3次元波数空間における、送信周波数帯域によって決定された波数kの上限値kmaxを半径とし原点Oを中心とした3次元球S1と、下限値kminを半径とし原点Oを中心とした3次元球S2とを示すイメージ図である。Wは、波数kの帯域幅である。
 k/2軸において、信号のサポートは、特定範囲に制限される。例えば、図3に示すx方向(アジマス方向)のアンテナ素子指向性の3dB半値幅を±θ3dBとし、この範囲でのみ、信号が得られると仮定した場合、波数kのサポート範囲は、k∈[kx,min,kx,max]=[-kmaxsinθ3dB,kmaxsinθ3dB]によって制限される。図5には、このサポート範囲がk/2軸と平行な矢印Rで描かれている。
 図6は、図5における波数kに対応した点Cを通り、k’-k’平面と平行な2つの3次元球の断面を示すイメージ図である。断面S1aは、サポート範囲にある波数kに対応した点C(k/2,0,0)を通り、k’-k’平面と平行な3次元球S1の断面である。断面S2aは、点C(k/2,0,0)を通り、k’-k’平面と平行な3次元球S2の断面である。
 図6に示すように、+k’方向に対する角度がθincで、点Cから、断面S1aおよびS2aの各外周側に向けて延びた直線Lと断面S2aの外周円との交点をAとし、直線Lと断面S1aの外周円との交点をBとする。これにより、交点Aを中心として点Cを通る円D1と、交点Bを中心として点Cを通る円D2が定義される。
 円D1の半径は、点Cと交点Aとの距離kyz,minであり、円D2の半径は、点Cからと交点Bとの距離kyz,maxである。円D1の外周側にありかつ円D2の内周側にある領域のうち、k/2軸方向に扇状に広がった領域Eは、波数kに関するk’-k’平面での周波数領域サポートである。
 図7は、波数kに関するk’-k’平面における周波数領域を示すイメージ図である。扇状の領域Eは、図7に示すように、k/2軸方向において、k/2∈[kx,min/2,kx,max/2]=[-(kmax/2)sinθ3dB,(kmax/2)sinθ3dB]の範囲で制限されている。これは、k/2軸方向、すなわち、ドップラ周波数軸の空いた周波数領域を利用することによって、送信信号を多重化し、MIMOを実現できる可能性があることを意味する。
 ドップラ周波数軸で信号を多重化するためには、各送信アンテナが送信する送信信号のサポートを、ドップラ周波数軸上で異なるシフト量を持つようにすればよい。それを実現するためには、異なる送信アンテナが送信する送信信号に対して、チャープごとに異なる初期位相変化量を設定すればよい。
 ターゲットが単位距離を移動する間の初期位相変化量をa(rad/m)とした場合、空間的に広がりのあるターゲットg(x,y,z)において送信信号が反射して得られる反射信号S(x’,y’,k)は、下記式(23)で表される。上記式(7)で表されるS(x’,y’,k)との違いは、グリーン関数に対してe-jax’が乗算されている部分である。なお、下記式(23)において、l(x’,y’)は、上記式(7)で表され、l(x’)は、上記式(9)で表され、l(x’,y’)は、上記式(10)で表される。

Figure JPOXMLDOC01-appb-I000011
 ターゲットg(x,y,z)は、反射信号S(x’,y’,k)を用いて復元される。例えば、上記式(11)から上記式(22)までにおいて、波数kに関する部分を波数k’に変更することで、ターゲットg(x,y,z)は、下記式(24)で表される。ただし、下記式(24)において、k’は、下記式(25)で定義され、k’は、上記式(12)で定義され、k’は、上記式(13)で定義され、kは、上記式(14)で定義される。ky,incは、上記式(15)で定義され、kz,incは、上記式(16)で定義される。kyzは、下記式(26)で定義される。cosθincは、上記式(18)で定義され、sinθincは、上記式(19)で定義される。

Figure JPOXMLDOC01-appb-I000012
 上記式(24)におけるΦは、下記式(27)で表される。Rは、上記式(21)で定義され、Rは、上記式(22)で定義される。

Figure JPOXMLDOC01-appb-I000013
 上記式(24)においては、Stolt[・]関数が、(k,k,k)→(k’,k’,k’)へのStolt補間処理を示す関数になっている。上記式(24)から上記式(27)の式は、ドップラ周波数軸上で3次元サポートが初期位相変化量a(rad/m)だけシフトすることを意味する。
 図8は、送信アンテナTx1と送信アンテナTx2によって得られる波数kに関するk’-k’平面における周波数領域を示すイメージ図である。送信アンテナTx1および送信アンテナTx2のうち、送信アンテナTx1は、チャープ間の初期位相を変化させず、送信アンテナTx2は、チャープ間で初期位相を初期位相変化量a(rad/m)で変化させた場合、周波数サポートは、k/2-k’の2次元平面において、図8に示す扇状の領域E1とE2である。このように、複数の送信アンテナから同時送信する際に、送信アンテナごとに異なる初期位相変化量を与えた場合、各送信アンテナが送信した信号は、ドップラ周波数軸上で分離され、DDM-MIMOが実現可能になる。
 図8に示すように、送信アンテナTx1が送信した信号と、送信アンテナTx2が送信した信号とを、ドップラ周波数軸上で綺麗に分離するためには、下記の関係式(28)を満たす必要がある。

Figure JPOXMLDOC01-appb-I000014
 例えば、送信アンテナTx2が送信した信号に設定された初期位相を、チャープ間でφ(rad)ずつ回転させ、ターゲットのx方向に沿った移動の速度をv(m/s)とし、チャープ間隔をT(s)とした場合に、初期位相変化量a(rad/m)は、下記式(29)で表される。

Figure JPOXMLDOC01-appb-I000015
 レーダ信号処理器2においてDDM-MIMOを実現するためには、チャープ間の位相回転量φ[rad]に関して、以下の条件を示す関係式(30)を満たす必要がある。

Figure JPOXMLDOC01-appb-I000016
 図3に示したジオメトリにおいて、4個の送信アンテナTx1~Tx4のそれぞれから送信された信号の周波数サポートが互いに重ならない条件を満たしつつ、送信アンテナTx1~Tx4のそれぞれから送信された信号のチャープ系列の初期位相を、図2に示したように変化させる。これにより、4送信DDM-MIMOが実現される。
 次に、レーダ信号処理器2による4送信DDM-MIMO処理の詳細は、以下の通りである。図9は、実施の形態1に係るレーダ画像生成方法を示すフローチャートであって、図2のジオメトリにおける4送信DDM-MIMO処理を利用した、ターゲットの3次元レーダ画像の生成を示している。
 図9の処理の前提として、制御部21が、4個の送信アンテナTx1~Tx4から同時に送信信号を送信し、送信信号が反射された反射信号を受信アンテナRxが受信し、受信信号のデジタルデータを取得する一連の処理を繰り返すように送信部および受信部を制御する。信号処理部23が、1次元アレイアンテナ8のアンテナ配列方向を横切って移動したターゲット9について順次取得されたデジタルデータを、2次元アレイデータとみなして逆合成開口処理を行う。すなわち、ステップST1からステップST9の一連の処理が逆合成開口処理であり、この逆合成開口処理によって、ターゲットの3次元レーダ画像が生成される。
 制御部21が、送信アンテナTx1~Tx4から初期位相変化量aが送信アンテナごとに異なる送信信号を同時に送信することにより、送信信号を、ドップラ周波数軸で多重化する。信号処理部23が、送信アンテナごとの送信信号に対応する反射信号の成分に基づいて、ターゲットの3次元レーダ画像を送信アンテナごとに生成し、送信アンテナごとの3次元レーダ画像を加算することにより、ターゲットの最終的な3次元レーダ画像を生成する。
 例えば、信号処理部23は、x方向に沿って移動するターゲットの各移動位置において送信信号が反射された反射信号を各受信アンテナが受信した時間信号に対して距離FFTを行うことにより、3次元データS(x’,y’,k)を生成する(ステップST1)。信号処理部23は、3次元データS(x’,y’,k)に対してx’方向およびy’方向に2次元フーリエ変換を行うことにより、S(k,k,k)を生成する(ステップST2)。4個の送信アンテナTx1~Tx4が送信する信号は、k軸上で分離されているため、各送信アンテナTx1~Tx4が送信する信号について、以下のステップST3~ステップST8までの処理が互いに独立して行われる。
 信号処理部23は、送信アンテナTx1~Tx4が送信した各信号のS(k,k,k)に対して、受信経路のバルク圧縮処理を行うことにより、送信アンテナごとにS’(k,k,k)を生成する(ステップST3)。これ以降、ステップST4からステップST7までの処理は、送信信号の平面波近似が必要である。このため、これらの処理は、y-z空間を区分した複数の2次元セルにおける各2次元空間で互いに独立して行われる。
 信号処理部23は、複数の2次元セルのそれぞれのy-z空間のうち、処理対象として選択したy-z空間において、S’(k,k,k)に対して、送信経路のバルク圧縮を行うことにより、S’’(k,k,k|ycnt,zcnt)を生成する(ステップST4)。続いて、信号処理部23は、S’’(k,k,k|ycnt,zcnt)に対して、バイスタティック配置補正を行うことで、S’’’(k,k,k|ycnt,zcnt)を生成する(ステップST5)。
 次に、信号処理部23は、S’’’(k,k,k|ycnt,zcnt)に対して、Stolt補間を行うことにより、S’’’’(k’,k’,k’|ycnt,zcnt)を生成する(ステップST6)。信号処理部23は、S’’’’(k’,k’,k’|ycnt,zcnt)に対して3次元逆フーリエ変換を行う(ステップST7)。これにより、処理対象として選択された2次元セルのy-z空間における3次元レーダ画像g(x,y,z|ycnt,zcnt)が生成される。以下、3次元レーダ画像g(x,y,z|ycnt,zcnt)は、再生画像g(x,y,z|ycnt,zcnt)と適宜記載される。
 信号処理部23は、y-z空間を区分した2次元セルごとに算出された複数の再生画像g(x,y,z|ycnt,zcnt)を結合する(ステップST8)。これにより、送信アンテナごとの送信信号に対応した再生画像g(x,y,z|Tx)が生成される。次に、信号処理部23は、送信アンテナTx1~Tx4から送信された全ての送信信号に関する送信アンテナごとの再生画像g(x,y,z|Tx)をコヒーレントに加算し、最終的な再生画像g(x,y,z)を出力する(ステップST9)。
 図9において、送信アンテナごとの処理(ステップST3からステップST8)は、互いに独立して行うことが可能である。また、2次元セルのy-z小セル空間ごとの処理(ステップST4からステップST7)についても互いに独立して処理することが可能である。これらの箇所は、マルチコアCPUまたはGPUを用いた並列処理によって高速な処理が可能である。
 図10は、1次元アレイアンテナ8とターゲット9との位置関係の一例を示す概要図である。図10に示すジオメトリにおいて、送信アンテナTx1による点ターゲットの画像再生結果と、送信アンテナTx1~Tx4を用いた4送信DDM-MIMOによる画像再生結果との比較を行う。ターゲット9の動きは、図3と同じであり、基板8aのアンテナ面からターゲット9の移動軌跡までの距離は2(m)である。
 4個の送信アンテナTx1~Tx4は、それぞれが等間隔に並んだリニアアレイを構成し、送信アンテナの間隔は、6(mm)である。16個の受信アンテナRx1,Rx2,・・・,Rx16においても、それぞれが等間隔に並んだリニアアレイを構成し、受信アンテナの間隔は、8(mm)である。1次元アレイアンテナ8は、4送信16受信のMIMO構成であり、十分遠方からみたときに、素子アンテナが2(mm)間隔で、総開口が約12(cm)の仮想アレイとみなされる。
 変調方式は、リニアFMCW変調であるものとし、周波数帯域は77~81(GHz)である。また、送受信アンテナの水平指向性幅は±13(deg)であるものとし、それよりも広角の範囲からの信号は受信されない。各送信アンテナに設定されたチャープごとの初期位相の変化シーケンスは、図2と同じである。
 図11は、3次元スペクトラムS(k,k,k)がk-kの2次元平面に投影された2次元スペクトラムを示すスペクトラム図であり、図11の下段の図は、上段の図で符号Fを付与したスペクトラム部分を示している。図11における2次元スペクトラムは、図9のステップST2において、信号処理部23によって算出された3次元スペクトラムS(k,k,k)を、k-kの2次元平面に全投影した2次元スペクトラムである。図11の上段に示すように、送信アンテナTx1~Tx4から送信された信号は、k軸において送信アンテナごとに分離されており、レーダ装置1では、4送信DDM-MIMOが実現される。また、図11の下段に示すように、個々の送信アンテナに対応する2次元スペクトラムが生成されている。
 図12Aは、送信アンテナTx1の送信信号に由来する受信信号に基づいて算出された画像強度分布を、y-z面で切った断面を示す断面図である。図12Aに示す画像強度分布は、3次元レーダ画像の再生結果である。図12Bは、図12Aの画像強度分布を、x-z面で切った断面を示す図である。図12Cは、図12Aの画像強度分布を、x-y面で切った断面を示す断面図である。図12Dは、ターゲット位置を通り、y軸に沿った直線上の図12Aの画像強度分布を示す分布図である。
 図13Aは、送信アンテナTx1~Tx4の送信信号に由来する受信信号に基づいて算出された画像強度分布を、y-z面で切った断面を示す断面図である。図13Bは、図13Aの画像強度分布を、x-z面で切った断面を示す図である。図13Cは、図13Aの画像強度分布を、x-y面で切った断面を示す断面図である。図13Dは、ターゲット位置を通り、y軸に沿った直線上の図13Aの画像強度分布を示す分布図である。
 図12A、図12B、図12Cおよび図12Dと図13A、図13B、図13Cおよび図13Dとから明らかなように、送信アンテナTx1~Tx4のそれぞれから送信された全ての信号に由来する受信信号に基づいて送信アンテナごとに算出された画像強度分布をコヒーレントに加算することにより、送信アンテナTx1の送信信号に由来する受信信号に基づいて算出された画像強度分布に比べてピーク電力が12dB(=20log104)程度改善している。これは、SN比が送信アンテナ数倍(=4倍)改善されたことを意味する。また、統合サイドローブ比(ISLR)も約4dB改善している。
 図10に示したオドメトリにおいて、MIMOであることから、総開口は変わらない。このため、レーダ装置1によって生成される3次元レーダ画像の分解能は変化しない。
 1次元アレイアンテナ8において、素子アンテナの間隔(MIMO仮想アレイ間隔)を大きくしつつ、かつ総開口が広がるように、送信アンテナおよび受信アンテナを配置することにより、3次元レーダ画像の分解能を向上させることが可能である。
 レーダ信号処理器2の機能を実現するハードウェア構成は、以下の通りである。
 図14Aは、レーダ信号処理器2の機能を実現するハードウェア構成を示すブロック図である。図14Bは、レーダ信号処理器2の機能を実現するソフトウェアを実行するハードウェア構成を示すブロック図である。図14Aおよび図14Bにおいて、入出力インタフェース100は、A/D変換器7によってA/D変換されたデジタルデータを入力し、制御部21からの制御信号を中継して送信部および受信部へ出力するインタフェースである。記憶装置101は、入出力インタフェース100が入力を受け付けたデータを記憶する記憶装置であり、データ記憶部22である。
 レーダ信号処理器2が備える制御部21および信号処理部23の機能は処理回路により実現される。すなわち、レーダ信号処理器2は、図9に示したステップST1からステップST9までの処理を実行する処理回路を備える。処理回路は、専用のハードウェアであってもよいが、メモリに記憶されたプログラムを実行するCPUであってもよい。
 処理回路が図14Aに示す専用のハードウェアの処理回路102である場合、処理回路102は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ASIC(Application Specific Integrated Circuit)、FPGA(Field-Programmable Gate Array)、または、これらを組み合わせたものが該当する。レーダ信号処理器2が備える制御部21および信号処理部23の機能は、別々の処理回路で実現されてもよく、これらの機能がまとめて1つの処理回路で実現されてもよい。
 処理回路が図14Bに示すプロセッサ103である場合、レーダ信号処理器2が備える制御部21および信号処理部23の機能は、ソフトウェア、ファームウェアまたはソフトウェアとファームウェアとの組み合わせにより実現される。なお、ソフトウェアまたはファームウェアは、プログラムとして記述されてメモリ104に記憶される。
 プロセッサ103は、メモリ104に記憶されたプログラムを読み出して実行することにより、制御部21および信号処理部23の機能を実現する。例えば、レーダ信号処理器2は、プロセッサ103によって実行されるときに、図9に示したステップST1からステップST9の処理が結果的に実行されるプログラムを記憶するメモリ104を備える。これらのプログラムは、制御部21と信号処理部23の手順または方法をコンピュータに実行させる。メモリ104は、コンピュータを、制御部21および信号処理部23として機能させるためのプログラムが記憶されたコンピュータ可読記憶媒体であってもよい。
 メモリ104は、例えば、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ、EPROM(Erasable Programmable Read Only Memory)、EEPROM(Electrically-EPROM)などの不揮発性または揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、DVDなどが該当する。
 制御部21および信号処理部23の機能の一部が専用のハードウェアで実現され、残りの一部がソフトウェアまたはファームウェアによって実現されてもよい。例えば、制御部21は、専用のハードウェアである処理回路102によって機能が実現され、信号処理部23は、プロセッサ103がメモリ104に記憶されたプログラムを読み出して実行することにより機能が実現される。このように、処理回路は、ハードウェア、ソフトウェア、ファームウェアまたはこれらの組み合わせによって、上記機能を実現することができる。
 上記式(11)から上記式(22)の導出は、以下の通りである。
 図3に示したジオメトリにおいて、空間的に広がりのあるターゲットg(x,y,z)が存在する場合、送信アンテナTx1が送信した送信信号の反射信号に関するスペクトルS(x’,y’,k)は、上記式(7)で表される。ただし、x’は、送信アンテナTx1によって送信信号が送信されたときのターゲット9の中心のx座標であり、y’は、この送信信号の反射信号を受信した受信アンテナRxのy座標-yの値に対応しており、送信信号の往復距離がl(x’,y’)であり、上記式(8)で表される。送信信号の往路距離がl(x’)であり、上記式(9)で表され、当該送信信号の反射信号の復路距離がl(x’,y’)であり、上記式(10)で表される。
 信号処理部23が、S(x’,y’,k)を、x’およびy’に関して2次元フーリエ変換して算出される周波数信号S(k,k,k)とg(x,y,z)との間は、下記式(31)の関係が成り立つ。

Figure JPOXMLDOC01-appb-I000017
 上記式(31)の右辺の最後の積分部分は、停留位相法(MOSP)によって算出することができ、下記式(32)で表される。

Figure JPOXMLDOC01-appb-I000018
 上記式(32)を、上記式(31)の右辺の最後の積分部分に代入することによって、S(k,k,k)とg(x,y,z)との間には下記式(33)の関係が成り立つ。

Figure JPOXMLDOC01-appb-I000019
 上記式(33)の右辺の第1項の最後の積分部分、第3項および第4項は、MOSPを用いて算出される。これらの指数部を-jφ(x’)とすると、φ(x’)は、下記式(34)のように、φ(x’)とφ(x’)に分解される。

Figure JPOXMLDOC01-appb-I000020
 φ(x’)を2次の項までテーラー近似した場合、dφ(x’)/dx’=0となるx’をx’とすると、x’は、下記式(35)で表される。ただし、kyzは、上記式(17)で表され、Rは、下記式(36)で表される。

Figure JPOXMLDOC01-appb-I000021
 上記式(35)により、φ(x’)は、下記式(37)で近似される。

Figure JPOXMLDOC01-appb-I000022
 同様に、φ(x’)を2次の項までテーラー近似した場合、dφ(x’)/dx’=0となるx’をx’とすると、x’は、下記式(38)で表される。ただし、kは、下記式(39)で表され、Rは、下記式(40)で表される。

Figure JPOXMLDOC01-appb-I000023
 上記式(38)により、φ(x’)は、下記式(41)で近似される。

Figure JPOXMLDOC01-appb-I000024
 従って、φ(x’)は、下記式(42)によって近似される。

Figure JPOXMLDOC01-appb-I000025
 上記式(42)の右辺の第3項および第4項を1次の項までテーラー近似した場合に、これらの項におけるx’に関する微分が0となるx’をx’とすると、x’は、下記式(43)で表される。ただし、Aは、下記式(44)で表され、Bは、下記式(45)で表される。

Figure JPOXMLDOC01-appb-I000026
 上記式(42)の右辺の第3項および第4項は、下記式(46)で表される。

Figure JPOXMLDOC01-appb-I000027
 従って、φ(x’)は、最終的に、下記式(47)で近似される。

Figure JPOXMLDOC01-appb-I000028
 φ(x’)が上記式(47)で表されることにより、S(k,k,k)とg(x,y,z)との間には、下記式(48)の関係が成り立つ。

Figure JPOXMLDOC01-appb-I000029
 上記式(48)の右辺における-jkyz{(y-y-Δy+(z-z1/2について、この平方根を取り除くために平面波近似が導入される。例えば、ターゲットg(x,y,z)のy-z平面での座標中心を(ycnt,zcnt)とした場合、上記平方根部分は、下記式(49)で示すように平面波近似できる。ただし、ky,incは、上記式(15)で定義され、kz,incは、上記式(16)で定義される。cosθincは、上記式(18)で定義され、sinθincは、上記式(19)で定義される。

Figure JPOXMLDOC01-appb-I000030
 上記式(49)を、上記式(48)の右辺における平方根部分である-jkyz{(y-y-Δy+(z-z1/2に代入して整理することにより、S(k,k,k)とg(x,y,z)との間には、最終的に、下記の近似関係式(50)が成り立つ。ただし、k’は、上記式(12)で定義され、k’は、上記式(13)で定義される。

Figure JPOXMLDOC01-appb-I000031
 従って、g(x,y,z)は、S(k,k,k)から、上記式(11)から上記式(22)までによって推定することが可能である。
 以上のように、実施の形態1に係るレーダ装置1は、直線状に配置された送信アンテナTx1~TxMを用いて空間に送信信号を送信し、送信アンテナと同一方向に直線状に配置された受信アンテナRx1~RxNを用いて空間で送信信号が反射された反射信号を受信し、送信アンテナTx1~TxMから同時に送信信号を送信し、受信アンテナRx1~RxNによって反射信号を受信してデジタルデータを取得する一連の処理を繰り返し実行させる。繰り返し実行された一連の処理で順次取得されたデジタルデータを、2次元アレイデータとして用いることにより、送信アンテナTx1~TxMおよび受信アンテナRx1~RxNのアンテナ配列方向を横切る方向に移動したターゲット9の3次元レーダ画像を生成する。これにより、送信アンテナTx1~TxMの信号送信間隔が短くなるので、信号送信間隔の時間内にターゲット9が移動することによる影響が低減される。例えば、送信アンテナ数をNとし、TDM-MIMO方式の従来のレーダ装置における一つの送信アンテナの信号送信間隔をT(s)とした場合、DDM-MIMO方式のレーダ装置1では、一つの送信アンテナの信号送信間隔を、T/N(s)に短くすることができる。従って、異なる送信アンテナからの送信信号に由来する受信信号をコヒーレントに加算することができ、3次元レーダ画像のSN比がN倍になるので、レーダ装置1は、3次元レーダ画像の分解能を高めることができる。
実施の形態2.
 図15は、実施の形態2に係るレーダ装置1Aの構成を示すブロック図であり、FMCW方式のレーダ装置を示している。図15において、図1と同一の構成要素には、同一の符号が付され、その説明が省略される。レーダ装置1Aは、レーダ装置1と同様に、1次元アレイアンテナを用いて、空間に送信信号を電波として送信し、送信した電波が空間に存在するターゲットで反射されて戻ってきた反射波を受信することにより、ターゲットの3次元レーダ画像を生成する。例えば、レーダ装置1Aは、リアルタイムに連続的にやってくる検知対象物を、その移動を阻害することなく、ボディチェックするミリ波セキュリティシステムに利用される。以下、レーダ装置1Aが、移動する検知対象物の3次元レーダ画像の再生を次々に行うミリ波セキュリティシステムであるものとして説明する。
 レーダ装置1Aは、レーダ信号処理器2A、送信信号生成器3、送信器4、受信器5、ビート信号生成器6およびA/D変換器7を備える。レーダ信号処理器2Aは、制御部21、データ記憶部22および信号処理部23Aを備えている。制御部21は、送信アンテナTx1~TxMから同時に送信する送信信号のパルス繰り返し間隔を一定とし、かつパルス繰り返し間隔の送信信号間の初期位相変化量を送信アンテナごとに異なる値とする。信号処理部23Aは、受信アンテナRx1~RxNによってパルス繰り返し間隔で逐次的に受信された反射信号のデジタルデータを、ドップラ周波数領域の信号へ逐次フーリエ変換することにより、検知対象物の3次元レーダ画像を生成する。
 実施の形態2に係るレーダ画像生成方法は、以下の通りである。
 図16は、実施の形態2に係るレーダ画像生成方法を示すフローチャートである。図16において、ステップST4aからステップST10aまでの処理は、図9に示したステップST3からステップST9までの処理と同様である。信号処理部23Aは、1列ずつ受信信号s(y’,k|h)を生成する(ステップST1a)。
 例えば、1次元フーリエ変換における逐次離散フーリエ変換は、離散信号x(n)(n=0,1,2,・・・)に対する、N点の短時間離散フーリエ変換が、下記式(51)で表される。ただし、ωは、下記式(52)で表される。

Figure JPOXMLDOC01-appb-I000032
 上記式(51)を変形することにより、逐次離散フーリエ変換の下記式(53)を導くことができる。

Figure JPOXMLDOC01-appb-I000033
 上記式(51)および上記式(53)を用いることにより、レーダ装置1Aは、3次元レーダ画像の生成に要する演算速度が高速化される。例えば、時間的に連続した3次元レーダ画像を再生する場合、直近の送信信号のチャープ番号をh番とし、チャープ番号(h-N+1)番からh番までのNチャープ分の受信信号を用いて、3次元レーダ画像を再生する。信号処理部23Aは、チャープ番号hの受信信号s(y’,k|h)を、y’に関して1次元フーリエ変換し、2次元周波数信号s(y’,k|h)を生成する(ステップST2a)。信号処理部23Aは、2次元周波数信号s(y’,k|h)の生成と並行して、チャープ番号(h-N+1)番からh番までの直近のNチャープ分の受信信号に対して3次元フーリエ変換を行う。これにより、3次元周波数信号S(k,k,k|h)が生成される。
 上記式(53)で表される逐次離散フーリエ変換を行うことにより、s(y,k|h)とS(k,k,k|h)との間には、下記式(54)、下記式(55)および下記式(56)の関係が成り立つ。ただし、下記式(54)、下記式(55)および下記式(56)において、Sij/h(k)=S(k,k=i,k=j|h)、sij/h=s(k=i,k=j|h)と略記している。信号処理部23Aは、下記式(54)、下記式(55)および下記式(56)を全てのk,kに対して適用することで、逐次FFTを実行する(ステップST3a)。これにより、送信アンテナTx1~Tx4のそれぞれのS(k,k,k)が生成される。以下の処理は、図9のステップST3以降と同様である。

Figure JPOXMLDOC01-appb-I000034
 以上のように、実施の形態2に係るレーダ装置1Aにおいて、信号処理部23Aは、受信アンテナRxによってパルス繰り返し間隔で逐次的に受信された反射信号のデジタルデータを、ドップラ周波数領域の信号へ逐次フーリエ変換することにより、ターゲット9の3次元レーダ画像を生成する。これにより、レーダ装置1Aは、ターゲット9の3次元レーダ画像を高速に生成することができる。
実施の形態3.
 図17は、実施の形態3に係るレーダ装置が備える3つの1次元アレイアンテナ8(1)~8(3)とターゲット9との位置関係を示す概要図である。検知対象物の3次元レーダ画像におけるy方向(垂直方向)の分解能を向上させるためには、MIMOの仮想開口長をy方向に長くする必要がある。MIMOの仮想開口は、図3に示した基板8aを、複数枚、y方向に密に並べて配置することにより実現できる。
 基板8a(1)~8a(3)の中央部分に、1つの1次元アレイアンテナにおける全ての送信アンテナが集中して配置されると、基板8a(1)~8a(3)をy方向に並べたMIMOの仮想アレイにおいて、仮想素子が存在しない大きなギャップが、複数箇所発生してしまう。このギャップが存在すると、例えば、3枚の基板8a(1)~8a(3)のそれぞれに設けられた1次元アレイアンテナ8(1)~(3)が受信した全ての受信信号をコヒーレントに加算することにより、再生された3次元レーダ画像に大きなリップルが生じてしまう。
 これに対して、実施の形態3に係るレーダ装置では、図17に示すように、各基板8aにおける第1の端部(+y方向の端部)と、第1の端部に対向する第2の端部(-y方向の端部)に少なくとも一つの送信アンテナが配置されている。1次元アレイアンテナ8(1)~8(3)によって構成されるMIMOの仮想アレイは、第1の端部と第2の端部に少なくとも一つの送信アンテナが配置された3枚の基板8aを、y方向に並べたものである。この構造を有することにより、レーダ装置1Aは、例えば、仮想素子が存在しない大きなギャップがMIMO仮想アレイ内に発生することを防止できる。
 図18Aは、基板8aの中央に全ての送信アンテナTx1~TxMが配置された1次元アレイアンテナ8を用いてターゲット9の3次元レーダ画像を再生した場合における、ターゲット9の点拡散関数に係る信号の特性を示す特性図である。図18Aにおいて、特性H1は、1次元アレイアンテナ8が一つ、すなわち、図3に示した基板8aが1枚であるレーダ装置1によってターゲット9の3次元レーダ画像を再生した場合における、ターゲット9の点拡散関数に係る信号の特性である。特性H2は、図17に示した1次元アレイアンテナ8(1)~8(3)を有したレーダ装置によってターゲット9の3次元レーダ画像を再生した場合における、ターゲット9の点拡散関数に係る信号の特性である。
 基板8a(1)~8a(3)の各中央部分に全ての送信アンテナTx1~TxMを配置した場合、図18Aに示すように、基板8a(1)~8a(3)に設けられた1次元アレイアンテナ8(1)~8(3)によって受信された全ての信号をコヒーレントに加算しても綺麗な点像が再現されず、メインローブに大きなリップルが生じてしまう。
 図18Bは、基板8aの両端に少なくとも一つの送信アンテナTxが配置された1次元アレイアンテナ8を用いて3次元レーダ画像を再生したターゲット9の点拡散関数に係る信号の特性を示す特性図である。図18Aにおいて、特性I1は、図3に示した基板8aが1枚であるレーダ装置1によってターゲット9の3次元レーダ画像を再生した場合における、ターゲット9の点拡散関数に係る信号の特性である。特性I2は、図17に示した1次元アレイアンテナ8(1)~8(3)を有したレーダ装置によってターゲット9の3次元レーダ画像を再生した場合における、ターゲット9の点拡散関数に係る信号の特性である。
 基板8a(1)~8a(3)のそれぞれの第1の端部と第2の端部に少なくとも一つの送信アンテナTxを配置すると、図18Bに示すように、1次元アレイアンテナ8(1)~8(3)によって受信された全ての信号をコヒーレントに加算することにより、メインローブのリップルが抑制され、コヒーレントが大きなピークが得られる。
 以上のように、実施の形態3に係るレーダ装置は、3個の1次元アレイアンテナ8(1)~8(3)を備え、1次元アレイアンテナ8(1)~8(3)は、直線状に配列される。1次元アレイアンテナ8(1)~8(3)は、送信アンテナTxのうちの少なくとも一つの送信アンテナが、基板8a(1)~8a(3)における第1の端部および第2の端部に配置されている。この構成を有することにより、メインローブのリップルを抑制しつつ、コヒーレントが大きなピークを得ることができる。
 実施の形態1から実施の形態3に示したレーダ装置は、エスカレータまたは動く歩道に乗っている人を検知対象とするシステムに利用可能である。例えば、エスカレータまたは動く歩道の両サイドに1つずつ、1次元アレイアンテナ8を配置することで、レーダ装置は、エスカレータまたは動く歩道にいる人が金属または非金属の禁制品を所持しているか否かを左右両側からモニタリングすることができる。
 また、実施の形態1から実施の形態3に示したレーダ装置は、人だけでなく手荷物検査に利用可能である。例えば、駅、空港などにおける手荷物検査装置では、荷物が、ベルトコンベアに置かれ、その状態でベルトコンベアが動作する。ベルトコンベアの両サイドに1次元アレイアンテナを配置することにより、ベルトコンベア上で移動している荷物のうちに金属または非金属の禁制品が含まれているか否かを、左右両側からモニタリングすることができる。
 さらに、実施の形態1から実施の形態3に示したレーダ装置は、1次元アレイアンテナを屋外の道路脇に配置することにより、その道路を通過した車、自転車または人が、金属または非金属の禁制品を所持しているか否かをモニタリングすることができる。
 なお、各実施の形態の組み合わせまたは実施の形態のそれぞれの任意の構成要素の変形もしくは実施の形態のそれぞれにおいて任意の構成要素の省略が可能である。
 本開示に係るレーダ装置は、例えば、イメージングレーダ装置に利用可能である。
 1,1A レーダ装置、2,2A レーダ信号処理器、3 送信信号生成器、4 送信器、5 受信器、6 ビート信号生成器、7 A/D変換器、8,8(1)~8(3) 1次元アレイアンテナ、8a,8a(1)~8a(3) 基板、9 ターゲット、21 制御部、22 データ記憶部、23,23A 信号処理部、31 VCO、32 位相制御部、33 パワーアンプ、41 送信リニアアレイ、51 受信リニアアレイ、61 LNA、62 分配回路、63 ミキサ、64 フィルタリング回路、100 入出力インタフェース、101 記憶装置、102 処理回路、103 プロセッサ、104 メモリ、Rx,Rx1~RxN,Rx1~Rx16 受信アンテナ、Tx,Tx1~TxM,Tx1~Tx4 送信アンテナ。

Claims (9)

  1.  直線状に配置された複数の送信アンテナを用いて空間に送信信号を送信する送信部と、
     複数の前記送信アンテナと同一方向に直線状に配置された複数の受信アンテナを用いて前記空間で前記送信信号が反射された反射信号を受信し、受信した信号のデジタルデータを取得する受信部と、
     前記送信部が複数の前記送信アンテナから同時に前記送信信号を送信し、前記受信部が複数の前記受信アンテナによって前記反射信号を受信して前記デジタルデータを取得する一連の処理を、前記送信部および前記受信部に繰り返し実行させる制御部と、
     繰り返し実行された前記一連の処理で順次取得された前記デジタルデータを2次元アレイデータとして用いることにより、複数の前記送信アンテナおよび複数の前記受信アンテナのアンテナ配列方向を横切る方向に移動した検知対象物の3次元レーダ画像を生成する信号処理部と、
     を備えたことを特徴とするレーダ装置。
  2.  前記送信部は、複数の前記送信アンテナから初期位相変化量が前記送信アンテナごとに異なる前記送信信号を同時に送信して前記送信信号をドップラ周波数軸で多重化し、
     前記信号処理部は、前記送信信号に対応する前記反射信号の成分に基づいて、前記検知対象物の前記3次元レーダ画像を前記送信アンテナごとに生成し、前記送信アンテナごとに生成した前記3次元レーダ画像を加算することにより、前記検知対象物の最終的な前記3次元レーダ画像を生成すること
     を特徴とする請求項1に記載のレーダ装置。
  3.  前記制御部は、複数の前記送信アンテナから同時に送信する前記送信信号のパルス繰り返し間隔を一定とし、かつ前記パルス繰り返し間隔の前記送信信号間の初期位相変化量を前記送信アンテナごとに異なる値とすること
     を特徴とする請求項2に記載のレーダ装置。
  4.  複数の前記受信アンテナは、等間隔に配列されていること
     を特徴とする請求項3に記載のレーダ装置。
  5.  前記信号処理部は、前記送信アンテナごとの前記検知対象物の前記3次元レーダ画像の生成を並列処理することにより、前記検知対象物の最終的な前記3次元レーダ画像を生成すること
     を特徴とする請求項2から請求項4のいずれか1項に記載のレーダ装置。
  6.  前記信号処理部は、前記受信アンテナによって前記パルス繰り返し間隔で逐次的に受信された反射信号の前記デジタルデータを、ドップラ周波数領域の信号へ逐次フーリエ変換することにより、前記検知対象物の3次元レーダ画像を生成すること
     を特徴とする請求項3に記載のレーダ装置。
  7.  複数の前記送信アンテナのうちの少なくとも一つの前記送信アンテナは、基板における第1の端部および前記第1の端部に対向する第2の端部に配置されていること
     を特徴とする請求項1記載のレーダ装置。
  8.  複数の前記送信アンテナおよび複数の前記受信アンテナが前記基板に配置された複数のアレイアンテナを備え、
     複数の前記アレイアンテナは、直線状に配列されていること
     を特徴とする請求項7に記載のレーダ装置。
  9.  直線状に配置された複数の送信アンテナを用いて空間に送信信号を送信する送信部と、
     複数の前記送信アンテナと同一方向に直線状に配置された複数の受信アンテナを用いて前記空間で前記送信信号が反射された反射信号を受信し、受信した信号のデジタルデータを取得する受信部とを備えたレーダ装置のレーダ画像生成方法であって、
     制御部が、前記送信部が複数の前記送信アンテナから同時に前記送信信号を送信し、前記受信部が複数の前記受信アンテナによって前記反射信号を受信して前記デジタルデータを取得する一連の処理を、前記送信部および前記受信部に繰り返し実行させるステップと、
     信号処理部が、繰り返し実行された前記一連の処理で順次取得された前記デジタルデータを2次元アレイデータとして用いることにより、複数の前記送信アンテナおよび複数の前記受信アンテナのアンテナ配列方向を横切る方向に移動した検知対象物の3次元レーダ画像を生成するステップと、
     を備え、
     前記送信部は、複数の前記送信アンテナから初期位相変化量が前記送信アンテナごとに異なる前記送信信号を同時に送信して前記送信信号をドップラ周波数軸で多重化し、
     前記信号処理部は、前記送信信号に対応する前記反射信号の成分に基づいて、前記検知対象物の前記3次元レーダ画像を前記送信アンテナごとに生成し、前記送信アンテナごとに生成した前記3次元レーダ画像を加算することにより、前記検知対象物の最終的な前記3次元レーダ画像を生成すること
     を特徴とするレーダ画像生成方法。
PCT/JP2020/023091 2020-06-11 2020-06-11 レーダ装置およびレーダ画像生成方法 WO2021250866A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2020555255A JP6921339B1 (ja) 2020-06-11 2020-06-11 レーダ装置およびレーダ画像生成方法
PCT/JP2020/023091 WO2021250866A1 (ja) 2020-06-11 2020-06-11 レーダ装置およびレーダ画像生成方法
DE112020007038.7T DE112020007038T5 (de) 2020-06-11 2020-06-11 Radareinrichtung und radarbild-generierungsverfahren
US17/966,041 US20230039504A1 (en) 2020-06-11 2022-10-14 Radar device and radar image generation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2020/023091 WO2021250866A1 (ja) 2020-06-11 2020-06-11 レーダ装置およびレーダ画像生成方法

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US17/966,041 Continuation US20230039504A1 (en) 2020-06-11 2022-10-14 Radar device and radar image generation method

Publications (1)

Publication Number Publication Date
WO2021250866A1 true WO2021250866A1 (ja) 2021-12-16

Family

ID=77269521

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/023091 WO2021250866A1 (ja) 2020-06-11 2020-06-11 レーダ装置およびレーダ画像生成方法

Country Status (4)

Country Link
US (1) US20230039504A1 (ja)
JP (1) JP6921339B1 (ja)
DE (1) DE112020007038T5 (ja)
WO (1) WO2021250866A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7205003B1 (ja) * 2022-03-17 2023-01-16 三菱電機株式会社 検知装置、センサ、検知システム、および、検知方法

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023074275A1 (ja) * 2021-10-29 2023-05-04 パナソニックIpマネジメント株式会社 レーダ装置
WO2023152884A1 (ja) * 2022-02-10 2023-08-17 三菱電機株式会社 画像生成装置及びシステム
WO2023162188A1 (ja) * 2022-02-28 2023-08-31 三菱電機株式会社 信号処理装置、信号処理方法及びレーダ装置
JP2023141038A (ja) * 2022-03-23 2023-10-05 パナソニックIpマネジメント株式会社 レーダ装置
WO2024100884A1 (ja) * 2022-11-11 2024-05-16 三菱電機株式会社 物体透視装置および物体透視方法
CN116594014B (zh) * 2023-04-28 2024-02-09 南京航空航天大学 一种稀疏Mosaic成像模式的实现方法及实现装置
CN116541799B (zh) * 2023-05-06 2023-11-07 航科院中宇(北京)新技术发展有限公司 航道结构数据精度提升方法、系统、电子设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106546981A (zh) * 2016-10-24 2017-03-29 复旦大学 运动人体安检成像系统和方法
JP2017223575A (ja) * 2016-06-16 2017-12-21 日本信号株式会社 物体検出装置
WO2019024170A1 (zh) * 2017-07-31 2019-02-07 深圳市无牙太赫兹科技有限公司 毫米波实时成像安全检测系统及安全检测方法
JP2019052952A (ja) * 2017-09-15 2019-04-04 株式会社デンソー レーダ装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9733340B2 (en) * 2014-11-21 2017-08-15 Texas Instruments Incorporated Techniques for high arrival angle resolution using multiple nano-radars

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017223575A (ja) * 2016-06-16 2017-12-21 日本信号株式会社 物体検出装置
CN106546981A (zh) * 2016-10-24 2017-03-29 复旦大学 运动人体安检成像系统和方法
WO2019024170A1 (zh) * 2017-07-31 2019-02-07 深圳市无牙太赫兹科技有限公司 毫米波实时成像安全检测系统及安全检测方法
JP2019052952A (ja) * 2017-09-15 2019-04-04 株式会社デンソー レーダ装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
ZHU, ZHANYU ET AL.: "Demonstration of 3-D Security Imaging at 24 GHz With a 1-D Sparse MIMO Array", IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, vol. 17, no. 12, December 2020 (2020-12-01), pages 1 - 5, XP011822534, DOI: 10.1109/LGRS.2019.2962618 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7205003B1 (ja) * 2022-03-17 2023-01-16 三菱電機株式会社 検知装置、センサ、検知システム、および、検知方法
WO2023175812A1 (ja) * 2022-03-17 2023-09-21 三菱電機株式会社 検知装置、センサ、検知システム、および、検知方法

Also Published As

Publication number Publication date
JPWO2021250866A1 (ja) 2021-12-16
DE112020007038T5 (de) 2023-01-26
JP6921339B1 (ja) 2021-08-18
US20230039504A1 (en) 2023-02-09

Similar Documents

Publication Publication Date Title
WO2021250866A1 (ja) レーダ装置およびレーダ画像生成方法
CN110622026B (zh) 用于环境检测的方法和装置
Yanik et al. Development and demonstration of MIMO-SAR mmWave imaging testbeds
US20210184367A1 (en) Radar apparatus
US11656325B2 (en) Methods and apparatus to realize scalable antenna arrays with large aperture
US10795009B2 (en) Digital beamforming for radar sensing using wireless communication chipset
JP4665590B2 (ja) 干渉型レーダ
CN111175750A (zh) 一种合成孔径雷达的成像方法及装置、设备、存储介质
Bleh et al. A 100 GHz FMCW MIMO radar system for 3D image reconstruction
CA3056070A1 (en) Systems and methods for performing synthetic aperture radar imaging using superchirps
Jeon et al. W-band FMCW MIMO radar system for high-resolution multimode imaging with time-and frequency-division multiplexing
Garry et al. Investigations toward multistatic passive radar imaging
Kitamura et al. Doppler division multiplexed multiple-input–multiple-output imaging using cascaded millimeter-wave radars
Gottinger et al. Coherent signal processing for loosely coupled bistatic radar
KR102053881B1 (ko) 지상배열안테나 시스템 및 지상배열안테나 시스템에서 상공의 탐지 영역의 이미지를 획득하는 방법
Herschel et al. MIMO imaging for next generation passenger security systems
Rankin et al. Radar imaging: Conventional and MIMO
Wojaczek et al. Range compression strategies for passive radar on airborne platforms
Soumekh Bistatic synthetic aperture radar imaging using wide-bandwidth continuous-wave sources
WO2023162188A1 (ja) 信号処理装置、信号処理方法及びレーダ装置
Zhang et al. A novel high-resolution and wide-swath SAR imaging mode
Garry Imaging methods for passive radar
Matsuda et al. DDMA-MIMO Observations With the MU Radar: Validation by Measuring a Beam Broadening Effect
Rhoads et al. A scalable W-band imaging radar development platform
Nan et al. Piecewise constant Doppler algorithm: Performance analysis, further simplification, and motion compensation

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2020555255

Country of ref document: JP

Kind code of ref document: A

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20940117

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 20940117

Country of ref document: EP

Kind code of ref document: A1