[go: up one dir, main page]

KR102641300B1 - Target detection method and system based on phase information of MIMO radar signal - Google Patents

Target detection method and system based on phase information of MIMO radar signal Download PDF

Info

Publication number
KR102641300B1
KR102641300B1 KR1020210171594A KR20210171594A KR102641300B1 KR 102641300 B1 KR102641300 B1 KR 102641300B1 KR 1020210171594 A KR1020210171594 A KR 1020210171594A KR 20210171594 A KR20210171594 A KR 20210171594A KR 102641300 B1 KR102641300 B1 KR 102641300B1
Authority
KR
South Korea
Prior art keywords
signal
phase
sequence
target
antenna
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
KR1020210171594A
Other languages
Korean (ko)
Other versions
KR20230083525A (en
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 KR1020210171594A priority Critical patent/KR102641300B1/en
Publication of KR20230083525A publication Critical patent/KR20230083525A/en
Application granted granted Critical
Publication of KR102641300B1 publication Critical patent/KR102641300B1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/52Discriminating between fixed and moving objects or between objects moving at different speeds
    • G01S13/536Discriminating between fixed and moving objects or between objects moving at different speeds using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves
    • 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/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/583Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
    • 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/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/292Extracting wanted echo-signals
    • G01S7/2923Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
    • G01S7/2927Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods by deriving and controlling a threshold value
    • 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/35Details of non-pulse systems
    • G01S7/352Receivers
    • G01S7/356Receivers involving particularities of FFT processing

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

MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법 및 시스템이 개시된다. 복수의 송신 안테나를 통해 송출된 레이더 송신신호가 표적에 의해 반사되어 복수의 수신 안테나를 통해 수신될 때, 각각의 수신신호들에 대한 소정의 신호처리를 통해 각 수신신호의 거리-속도 검출 결과에서의 피크를 찾아내어 해당 피크의 위상정보를 추출한다. 복수의 수신 안테나 전체의 수신신호들로부터 추출된 거리-속도 피크들의 위상 값들을 등차수열을 이루도록 배열하여 그 등차수열로 배열된 위상 값 시퀀스에 기반하여 상기 수신신호들을 각 송신 안테나별로 구분한다. 등차수열로 배열된 위상 값을 이용하여 표적의 올바른 속도를 추정하고 각도 추정 또한 동시에 수행할 수 있다.A target detection method and system based on phase information of MIMO radar signals are disclosed. When a radar transmission signal transmitted through a plurality of transmitting antennas is reflected by a target and received through a plurality of receiving antennas, the distance-speed detection result of each received signal is obtained through predetermined signal processing for each received signal. Find the peak and extract the phase information of the peak. The phase values of the distance-velocity peaks extracted from the received signals of all the plurality of receiving antennas are arranged to form an arithmetic sequence, and the received signals are classified for each transmitting antenna based on the phase value sequence arranged in the arithmetic sequence. Using phase values arranged in an arithmetic sequence, the correct speed of the target can be estimated and angle estimation can be performed simultaneously.

Description

MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법 및 시스템 {Target detection method and system based on phase information of MIMO radar signal}Target detection method and system based on phase information of MIMO radar signal}

본 발명은 레이더 시스템을 이용한 표적 탐지 기술에 관한 것으로, 보다 상세하게는 다중 입력 다중 출력(Multi Input Multi Output: MIMO) 안테나 기반 주파수 변조 연속파(frequency-modulated continuous wave: FMCW) 레이더 시스템에서 각 송신 안테나로부터 수신된 신호들을 송신 안테나별로 구분하여 표적을 정확하게 탐지할 수 있도록 하는 기술에 관한 것이다. The present invention relates to target detection technology using a radar system, and more specifically, to each transmitting antenna in a frequency-modulated continuous wave (FMCW) radar system based on a multiple input multiple output (MIMO) antenna. This relates to a technology that allows accurate detection of targets by classifying signals received from each transmitting antenna.

FMCW 레이더는 거리, 속도, 각도와 같은 주변 대상에 대한 정보를 얻기 위한 자동차 센서로 널리 사용된다. 또한, MIMO 안테나 시스템은 제한된 하드웨어 크기로 높은 각도 분해능을 얻기 위해 활용되고 있다. MIMO 기반 FMCW 레이더 시스템에 있어서, 올바르게 동작하기 위해서는 다중 송신(Tx) 안테나에서 전송된 신호들을 수신(Rx) 안테나 측에서 구별할 수 있어야 한다. 이를 위해 다양한 신호 다중화 방식이 사용되고 있다. 예컨대, 시분할 다중화(Time Division Multiplexing: TDM), 주파수 분할 다중화(Frequency Division Multiplexing: FDM), 코드 분할 다중화(Code Division Multiplexing: CDM)와 같은 다양한 신호 다중화 방법이 전송된 신호의 직교성을 유지하는 데 사용할 수 있다. 그 중에서도 각 송신 안테나가 시간 슬롯을 나누어 다른 시간에 송신하는 TDM 방식이 널리 활용되고 있다. FMCW radar is widely used as an automotive sensor to obtain information about surrounding objects such as distance, speed, and angle. Additionally, MIMO antenna systems are being utilized to obtain high angular resolution with limited hardware size. In a MIMO-based FMCW radar system, in order to operate correctly, the receive (Rx) antenna must be able to distinguish between signals transmitted from multiple transmit (Tx) antennas. For this purpose, various signal multiplexing methods are used. For example, various signal multiplexing methods such as Time Division Multiplexing (TDM), Frequency Division Multiplexing (FDM), and Code Division Multiplexing (CDM) can be used to maintain orthogonality of transmitted signals. You can. Among them, the TDM method, in which each transmission antenna divides time slots and transmits at different times, is widely used.

TDM 방식은 구현이 용이한 장점이 있다. 그러나 TDM 방식의 가장 큰 단점은 처프 반복 시간이 길어져 송신 안테나의 개수에 비례하여 검출 가능한 최대 속도 즉, 최대 탐지 속도가 감소한다는 점이다. 또한, 각 송신 안테나는 교대로 신호를 전송하기 때문에 검출하고자 하는 표적이 움직이면 그 표적의 움직임에 의한 도플러 효과로 인해 각도 추정 성능이 저하될 수 있는 단점도 있다.The TDM method has the advantage of being easy to implement. However, the biggest disadvantage of the TDM method is that the chirp repetition time becomes longer, which reduces the maximum detectable speed, that is, the maximum detection speed, in proportion to the number of transmit antennas. Additionally, because each transmitting antenna transmits signals alternately, there is a disadvantage that if the target to be detected moves, angle estimation performance may deteriorate due to the Doppler effect caused by the movement of the target.

이러한 이슈를 해결하기 위해 TDM-MIMO 레이더 시스템의 도플러-각도 모호성(Doppler-angle ambiguity)을 해결하려고 몇 가지 제안이 있었다. 아래 비특허문헌 1에서는 통상적인 처프 시퀀스가 각 처프마다 다른 시작 주파수를 할당하도록 수정되었다. 유사하게, 아래 비특허문헌 2에서는 여러 종류의 처프 반복 주파수가 할당되었고 원래의 최대 감지 가능한 속도는 Chinese remainder theorem을 사용하여 복구되었다. 그러나 이러한 방법과 관련된 주요 단점은 기존 하드웨어를 수정해야 한다는 것이다.To solve this issue, several proposals have been made to resolve the Doppler-angle ambiguity of TDM-MIMO radar systems. In Non-Patent Document 1 below, the typical chirp sequence has been modified to assign a different start frequency to each chirp. Similarly, in Non-Patent Document 2 below, several types of chirp repetition frequencies were assigned and the original maximum detectable rate was recovered using the Chinese remainder theorem. However, the main disadvantage associated with these methods is that existing hardware must be modified.

또 다른 접근 방식은 각 송신 안테나에 특정 코드를 할당하고 신호를 동시에 전송하는 것이다. 특히 이진 위상천이변조 (Binary Phase-Shift Keying: BPSK) 방식은 전송 신호의 위상을 0 또는 π로 변조하는 방식으로 널리 사용된다. 이 방법은 송신된 신호의 직교성을 이용하여 각 송신 안테나의 신호를 분리할 수 있다. TDM-MIMO 시스템과 달리 BPSK-MIMO 시스템은 각 송신 안테나로부터의 신호를 분리하기 위해 추가적인 디코딩 과정이 필요하다. 그러나 송신 안테나의 완전한 전송 기능을 사용하여 더 높은 신호 대 잡음비(SNR)를 달성할 수 있다. Another approach is to assign a specific code to each transmit antenna and transmit the signals simultaneously. In particular, Binary Phase-Shift Keying (BPSK) is a widely used method of modulating the phase of a transmission signal to 0 or π. This method can separate the signals of each transmit antenna by using the orthogonality of the transmitted signal. Unlike the TDM-MIMO system, the BPSK-MIMO system requires an additional decoding process to separate the signals from each transmit antenna. However, higher signal-to-noise ratios (SNR) can be achieved by using the full transmit capability of the transmit antenna.

그런데, BPSK-MIMO 시스템 역시 전송 안테나의 수에 비례하여 속도 추정에 사용되는 블록 시간이 길어지기 때문에 도플러 각도 모호성(Doppler-angle ambiguity)을 겪는다. 즉, BPSK-MIMO 레이더 시스템의 거리-속도 탐지 결과에서 각 표적은 송신 안테나 수만큼 중복되어 나타난다. 달리 말하면, N개의 송신 안테나로 레이더 신호를 송출한 경우, 하나의 표적은 수신 신호의 피크 검출에서 하나의 실제 표적에 해당하는 피크와 (N-1)개의 고스트 표적에 해당하는 피크가 감지된다. BPSK 방식을 사용하는 경우 그 고스트 표적에 해당하는 피크 때문에 최대 탐지속도가 절반으로 줄어드는 문제가 발생한다. However, the BPSK-MIMO system also experiences Doppler-angle ambiguity because the block time used for speed estimation becomes longer in proportion to the number of transmission antennas. In other words, in the distance-speed detection results of the BPSK-MIMO radar system, each target appears overlapping as many as the number of transmitting antennas. In other words, when a radar signal is transmitted with N transmitting antennas, a peak corresponding to one actual target and a peak corresponding to (N-1) ghost targets are detected in the peak detection of the received signal. When using the BPSK method, a problem occurs in which the maximum detection speed is reduced by half due to the peak corresponding to the ghost target.

아래 비특허문헌 3에서는 다른 명확한 속도를 가진 다중 프레임이 사용되었고 원래의 최대 감지 가능한 속도는 Chinese remainder theorem 또는 클러스터링 알고리즘을 사용하여 복구되었다. 그러나 이 방법은 비특허문헌 1과 2에서 사용된 방법과 유사하게 기존 하드웨어를 수정해야 한다. 비특허문헌 3에서는 또한 추정된 속도 값을 사용하여 수신 신호의 위상을 보상하여 원래의 최대 감지 가능한 속도를 복원하려고 시도했다. 그러나 이런 식으로 위상을 직접 보정하는 방법은 추정된 속도 값에 작은 오차가 있는 경우 큰 오차를 초래할 수 있기 때문에 매우 불안정하다.In Non-Patent Document 3 below, multiple frames with different definite rates were used and the original maximum detectable rate was recovered using the Chinese remainder theorem or clustering algorithm. However, this method requires modifying existing hardware similar to the method used in Non-Patent Documents 1 and 2. Non-patent Document 3 also attempted to restore the original maximum detectable speed by compensating the phase of the received signal using the estimated speed value. However, this method of directly correcting the phase is very unstable because small errors in the estimated velocity value can lead to large errors.

또한, 미국 특허출원공개번호 2018-0011170(A1) (발명의 명칭: Methods and Apparatus for Velocity Detection in MIMO Radar Including Velocity Ambiguity Resolution)은 표적의 움직임으로 인해 발생한 도플러 위상을 추정한 뒤, 이 값만큼 각 송신 안테나의 위상 값을 보정해준 다음, 신호 처리 과정을 통해 표적의 속도 및 각도 정보를 추정하는 기술을 개시한다. 이 기술은 위상 값을 직접 보정해주는 경우 보정해주는 위상 값에 약간의 오차가 있는 경우에도 큰 오류가 발생할 수 있다.In addition, U.S. Patent Application Publication No. 2018-0011170(A1) (title of the invention: Methods and Apparatus for Velocity Detection in MIMO Radar Including Velocity Ambiguity Resolution) estimates the Doppler phase caused by the movement of the target and then adjusts the angle by this value. After correcting the phase value of the transmitting antenna, a technology for estimating the target's speed and angle information through signal processing is initiated. When this technology directly corrects the phase value, a large error may occur even if there is a slight error in the corrected phase value.

1. 미국 특허출원공개번호 2018-0011170 (A1)1. U.S. Patent Application Publication No. 2018-0011170 (A1)

1. K. Thurn, D. Shmakov, G. Li, S. Max, M. Meinecke, and M. Vossiek, "Concept and implementation of a PLL-controlled interlaced chirp sequence radar for optimized range-Doppler measurements," IEEE Trans. Microw. Theory Tech., vol. 64, no. 10, pp. 3280-3289, Oct. 2016.1. K. Thurn, D. Shmakov, G. Li, S. Max, M. Meinecke, and M. Vossiek, “Concept and implementation of a PLL-controlled interlaced chirp sequence radar for optimized range-Doppler measurements,” IEEE Trans . Microw. Theory Tech., vol. 64, no. 10, pp. 3280-3289, Oct. 2016. 2. Y. Li, C. Xu, X. Yan, and Q. Liu, "An improved algorithm for Doppler ambiguity resolution using multiple pulse repetition frequencies," in Int. Conf. Wirel. Commun. Signal Process., Nanjing, China, 2017, pp. 1-5.2. Y. Li, C. Xu, X. Yan, and Q. Liu, “An improved algorithm for Doppler ambiguity resolution using multiple pulse repetition frequencies,” in Int. Conf. Wirel. Commun. Signal Process., Nanjing, China, 2017, pp. 1-5. 3. H. A. Gonzalez, C. Liu, B. Vogginger, and C. G. Mayr, "Doppler ambiguity resolution for binary-phase-modulated MIMO FMCW radars," in Int.3. H. A. Gonzalez, C. Liu, B. Vogginger, and C. G. Mayr, “Doppler ambiguity resolution for binary-phase-modulated MIMO FMCW radars,” in Int. Radar Conf., Toulon, France, 2019, pp. 1-6.Radar Conf., Toulon, France, 2019, pp. 1-6.

이러한 문제를 해결하기 위해 본 발명은 MIMO 레이더 시스템에서 수신신호의 위상 정보를 이용하여 고스트 표적 신호를 배제하고 각 송신 안테나별로 실제 표적의 신호를 구별해내어 최대 탐지 속도의 감소 없이 도플러-각도 모호성을 해결하여 원래의 탐지 속도로 표적의 속도를 추정할 수 있음과 동시에, 표적의 움직임 여부에 관계없이 표적의 각도를 정확하고 효율적으로 추정할 수 있는 MIMO 레이더 기반 표적 탐지 방법 및 시스템을 제공하기 위한 것이다.To solve this problem, the present invention uses the phase information of the received signal in a MIMO radar system to exclude ghost target signals and distinguish the actual target signal for each transmit antenna, eliminating Doppler-angle ambiguity without reducing the maximum detection speed. The purpose is to provide a MIMO radar-based target detection method and system that can estimate the target's speed at the original detection speed and at the same time accurately and efficiently estimate the target's angle regardless of whether the target is moving. .

본 발명이 해결하고자 하는 과제는 상술한 과제들에 한정되는 것이 아니며, 본 발명의 사상 및 영역으로부터 벗어나지 않는 범위에서 다양하게 확장될 수 있을 것이다.The problems to be solved by the present invention are not limited to the above-described problems, and may be expanded in various ways without departing from the spirit and scope of the invention.

상기 본 발명의 목적을 실현하기 위한 실시예들에 따른 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법은, MIMO 레이더 시스템을 이용한 표적 탐지 방법에 있어서, 복수의 송신 안테나를 통해 송출된 레이더 송신신호가 표적에 의해 반사되어 되돌아오는 신호를 복수의 수신 안테나를 통해 동시에 수신하는 단계; 상기 복수의 수신 안테나 각각의 수신신호들에 대한 소정의 신호처리를 통해 각 수신신호의 거리-속도 검출 결과에서의 피크를 찾아내어 해당 피크의 위상정보를 추출하는 단계; 및 상기 복수의 수신 안테나 전체의 수신신호들로부터 추출된 거리-속도 피크들의 위상 값들을 등차수열을 이루도록 배열하여 그 등차수열로 배열된 위상 값 시퀀스에 기반하여 상기 수신신호들을 각 송신 안테나별로 구분하는 단계를 포함한다. A target detection method based on phase information of a MIMO radar signal according to embodiments for realizing the purpose of the present invention is a target detection method using a MIMO radar system, where a radar transmission signal transmitted through a plurality of transmission antennas is a target. Simultaneously receiving signals reflected and returned through a plurality of receiving antennas; finding a peak in the distance-velocity detection result of each received signal through predetermined signal processing on the received signals from each of the plurality of receiving antennas and extracting phase information of the corresponding peak; and arranging the phase values of the distance-velocity peaks extracted from the received signals of all the plurality of receiving antennas to form an arithmetic sequence and classifying the received signals for each transmitting antenna based on the phase value sequence arranged in the arithmetic sequence. Includes steps.

예시적인 실시예에 따르면, 상기 위상정보를 추출하는 단계는 상기 수신신호의 신호처리를 통해 얻어지는 cos 채널(I 채널)의 중간 주파수(IF) 대역 신호와 sin 채널(Q 채널)의 중간 주파수(IF) 대역 신호에 대하여 2차원 고속푸리에변환(2D FFT)을 적용하여 각 표적에 관한 거리-속도 검출 결과를 얻는 단계; 및 상기 거리-속도 검출 결과에 대하여 일정 오경보율 (constant false alarm rate: CFAR) 알고리즘을 적용하여 각 표적에 해당하는 거리-속도 피크들과 그 피크들의 위상값을 추출하는 단계를 포함할 수 있다.According to an exemplary embodiment, the step of extracting the phase information includes the intermediate frequency (IF) band signal of the cos channel (I channel) and the intermediate frequency (IF) of the sin channel (Q channel) obtained through signal processing of the received signal. ) Obtaining distance-velocity detection results for each target by applying a two-dimensional fast Fourier transform (2D FFT) to the band signal; and applying a constant false alarm rate (CFAR) algorithm to the distance-speed detection result to extract distance-velocity peaks corresponding to each target and phase values of the peaks.

예시적인 실시예에 따르면, 상기 레이더 송신신호는 PSK 변조방식으로 변조된 신호일 수 있다.According to an exemplary embodiment, the radar transmission signal may be a signal modulated by PSK modulation.

예시적인 실시예에 따르면, 상기 복수의 송신 안테나와 상기 복수의 수신 안테나는 ULA 형태로 배열된 안테나 어레이일 수 있다.According to an exemplary embodiment, the plurality of transmitting antennas and the plurality of receiving antennas may be an antenna array arranged in a ULA form.

예시적인 실시예에 따르면, 상기 각 송신 안테나별로 구분하는 단계에서, 상기 송신 안테나가 두 개일 때, 등차수열을 이루는 거리-속도 피크의 위상 시퀀스 중 선행 서브 위상 시퀀스(former phase sequence)는 제1 송신 안테나로부터 송신된 신호들에 대응하고, 후행 서브 위상 시퀀스(latter phase sequence)는 제2 송신 안테나로부터 송신된 신호들에 대응하는 것으로 결정될 수 있다. 여기서 각 서브 위상 시퀀스는 수신 안테나의 개수와 동일한 수의 위상들의 시퀀스이다.According to an exemplary embodiment, in the step of classifying each transmission antenna, when there are two transmission antennas, the preceding sub-phase sequence among the phase sequences of the distance-velocity peaks forming the arithmetic sequence is the first transmission Corresponding to signals transmitted from the antenna, a latter phase sequence may be determined to correspond to signals transmitted from the second transmit antenna. Here, each sub-phase sequence is a sequence of phases equal to the number of receiving antennas.

예시적인 실시예에 따르면, 상기 각 송신 안테나별로 구분하는 단계에서, 상기 송신 안테나가 N(N은 3이상의 자연수)개 있는 경우, 등차수열을 이루는 거리-속도 피크들의 위상 시퀀스 중 k 번째(여기서, k는 1부터 N까지의 자연수) 서브 위상 시퀀스는 제k 송신 안테나로부터 송신된 신호들에 대응하는 것으로 결정될 수 있다. 여기서 각 서브 위상 시퀀스는 수신 안테나의 개수와 동일한 수의 위상들의 시퀀스이다. According to an exemplary embodiment, in the step of classifying each transmission antenna, if there are N transmission antennas (N is a natural number of 3 or more), the kth phase sequence of distance-velocity peaks forming an arithmetic sequence (here, k is a natural number from 1 to N) The sub-phase sequence may be determined to correspond to signals transmitted from the kth transmission antenna. Here, each sub-phase sequence is a sequence of phases equal to the number of receiving antennas.

예시적인 실시예에 따르면, 상기 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법은, 상기 구분하는 단계 후, 각 송신 안테나별로 구분된 수신 신호를 이용하여 각 표적의 속도와 각도에 관한 추정정보를 산출하는 단계를 더 포함할 수 있다.According to an exemplary embodiment, the target detection method based on phase information of the MIMO radar signal calculates estimated information about the speed and angle of each target using the received signal divided for each transmit antenna after the classification step. Additional steps may be included.

예시적인 실시예에 따르면, 각 표적의 각도는 각 표적별 위상 시퀀스에 고속 푸리에변환을 적용하여 추정될 수 있다.According to an exemplary embodiment, the angle of each target may be estimated by applying fast Fourier transform to the phase sequence for each target.

예시적인 실시예에 따르면, 상기 레이더 송신신호는 시간에 따라 주파수가 증가하는 처프 신호일 수 있다.According to an exemplary embodiment, the radar transmission signal may be a chirp signal whose frequency increases with time.

한편, 상기 본 발명의 목적을 실현하기 위한 실시예들에 따른 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 시스템은, 복수의 송신 안테나, 복수의 수신 안테나, 송신부, 수신부, 디지털 신호 처리부, 거리/속도 추정부, 그리고 각도 추정부를 포함할 수 있다. 상기 송신부는 레이더 송신신호를 생성하여 상기 복수의 송신 안테나를 통해 무선 송신하도록 구성된다. 상기 수신부는 전방의 표적에 의해 반사되어 상기 복수의 수신 안테나를 통해 수신된 수신신호와 상기 송신신호 및 상기 송신신호의 90도 위상천이된 신호를 이용하여 cos 채널(I 채널)의 중간 주파수(IF) 대역 신호와 sin 채널(Q 채널)의 중간 주파수(IF) 대역 신호를 생성하도록 구성된다. 상기 디지털 신호 처리부는 상기 cos 채널(I 채널)의 중간 주파수 대역 신호와 sin 채널(Q 채널)의 중간 주파수 대역 신호에 대하여 필터링 및 고속 푸리에변환(FFT) 처리를 수행하도록 구성된다. 상기 거리/속도 추정부는 상기 디지털 신호 처리부에서 처리된 신호를 이용하여 거리-속도 피크들을 검출하고 각 피크의 위상 값을 추출하는 처리, 추출된 위상 값들을 등차수열을 이루도록 배열하고, 등차수열로 배열된 위상 값 시퀀스에 기반하여 수신신호들과 각 송신 안테나들 간의 대응관계를 결정하여 수신신호를 송신안테나별로 구분하는 처리, 그리고 각 송신안테나별로 구분된 수신신호를 이용하여 표적의 거리와 속도를 추정하는 처리를 수행하도록 구성된다. Meanwhile, a target detection system based on phase information of a MIMO radar signal according to embodiments for realizing the purpose of the present invention includes a plurality of transmitting antennas, a plurality of receiving antennas, a transmitter, a receiver, a digital signal processor, and a distance/speed detector. It may include a government, and an angle estimator. The transmitter is configured to generate a radar transmission signal and transmit it wirelessly through the plurality of transmission antennas. The receiver uses the received signal reflected by the target in front and received through the plurality of receiving antennas, the transmitted signal, and a 90-degree phase-shifted signal of the transmitted signal to an intermediate frequency (IF) of the cos channel (I channel). ) band signal and an intermediate frequency (IF) band signal of the sin channel (Q channel). The digital signal processing unit is configured to perform filtering and fast Fourier transform (FFT) processing on the intermediate frequency band signal of the cos channel (I channel) and the intermediate frequency band signal of the sin channel (Q channel). The distance/speed estimator detects distance-speed peaks using the signal processed by the digital signal processor and extracts the phase value of each peak, arranges the extracted phase values to form an arithmetic sequence, and arranges them in an arithmetic sequence. Based on the phase value sequence, the correspondence between the received signals and each transmitting antenna is determined, the received signal is classified by transmitting antenna, and the target distance and speed are estimated using the received signal classified by each transmitting antenna. It is configured to perform processing.

예시적인 실시예에 따르면, 거리/속도 추정부는 상기 디지털 신호 처리부에서 고속 푸리에변환(FFT) 처리된 신호에 대하여 일정 오경보율 (constant false alarm rate: CFAR)알고리즘을 적용하여 각 표적에 해당하는 거리-속도 피크들과 그 피크들의 위상값을 추출하도록 구성될 수 있다.According to an exemplary embodiment, the distance/speed estimator applies a constant false alarm rate (CFAR) algorithm to the signal processed by the fast Fourier transform (FFT) in the digital signal processor to determine the distance- It may be configured to extract velocity peaks and phase values of the peaks.

예시적인 실시예에 따르면, 상기 거리/속도 추정부는, 상기 송신 안테나가 두 개일 때, 등차수열을 이루는 거리-속도 피크의 위상 시퀀스 중 선행 서브 위상 시퀀스(former phase sequence)는 제1 송신 안테나로부터 송신된 신호들에 대응하고, 후행 서브 위상 시퀀스(latter phase sequence)는 제2 송신 안테나로부터 송신된 신호들에 대응하는 것으로 결정하도록 구성될 수 있다. 여기서 각 서브 위상 시퀀스는 수신 안테나의 개수와 동일한 수의 위상들의 시퀀스이다. According to an exemplary embodiment, when there are two transmission antennas, the distance/speed estimator is configured to transmit a preceding sub-phase sequence among the phase sequences of the distance-velocity peaks forming an arithmetic sequence from the first transmission antenna. Corresponding to the signals transmitted from the second transmitting antenna, a latter phase sequence may be configured to determine that corresponding to the signals transmitted from the second transmit antenna. Here, each sub-phase sequence is a sequence of phases equal to the number of receiving antennas.

예시적인 실시예에 따르면, 상기 거리/속도 추정부는, 상기 송신 안테나가 N(N은 3이상의 자연수)개 있는 경우, 등차수열을 이루는 거리-속도 피크들의 위상 시퀀스 중 k 번째(여기서, k는 1부터 N까지의 자연수) 서브 위상 시퀀스는 제k 송신 안테나로부터 송신된 신호들에 대응하는 것으로 결정하도록 구성될 수 있다. 여기서 각 서브 위상 시퀀스는 수신 안테나의 개수와 동일한 수의 위상들의 시퀀스이다.According to an exemplary embodiment, the distance/speed estimator, when there are N transmitting antennas (N is a natural number of 3 or more), the kth phase sequence of distance-speed peaks forming an arithmetic sequence (where k is 1) to N) may be configured to determine that the sub-phase sequence corresponds to signals transmitted from the kth transmit antenna. Here, each sub-phase sequence is a sequence of phases equal to the number of receiving antennas.

예시적인 실시예에 따르면, 상기 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 시스템은 각 송신안테나별로 구분된 수신신호를 이용하여 소정의 각도 추정 알고리즘을 적용하여 각 표적의 각도정보를 산출하도록 구성된 각도 추정부를 더 포함할 수 있다.According to an exemplary embodiment, the target detection system based on phase information of the MIMO radar signal includes an angle estimation unit configured to calculate the angle information of each target by applying a predetermined angle estimation algorithm using received signals divided for each transmission antenna. More may be included.

예시적인 실시예에 따르면, 상기 각도 추정부는 각 표적별 위상 시퀀스에 고속 푸리에변환을 적용하여 각 표적의 각도를 추정하도록 구성될 수 있다.According to an exemplary embodiment, the angle estimator may be configured to estimate the angle of each target by applying fast Fourier transform to the phase sequence for each target.

예시적인 실시예에 따르면, 상기 복수의 송신 안테나와 상기 복수의 수신 안테나는 ULA 형태로 배열된 안테나 어레이일 수 있다. According to an exemplary embodiment, the plurality of transmitting antennas and the plurality of receiving antennas may be an antenna array arranged in a ULA form.

예시적인 실시예에 따르면, 상기 레이더 송신신호는 시간에 따라 주파수가 증가하는 처프 신호일 수 있다.According to an exemplary embodiment, the radar transmission signal may be a chirp signal whose frequency increases with time.

BPSK-MIMO 레이더 시스템의 거리-속도 탐지 결과에서 각 표적은 송신 안테나 수만큼 중복되어 나타난다. 즉, 하나의 실제 표적과 여러 개의 고스트 표적이 감지된다. BPSK-MIMO 레이더 시스템의 수신 신호에서 실제 표적과 고스트 표적을 구별하기 위해 실제 표적과 고스트 표적에 해당하는 피크 신호의 위상 값이 등차수열의 형태를 가져야 한다는 아이디어를 사용한다. 따라서 수신신호에서 이러한 피크의 위상 값을 등차수열로 제대로 배열하면 각 표적이 어느 송신 안테나로부터 생성되었는지를 결정할 수 있다. 즉, 수신 신호에서 실제 표적에 해당하는 피크와 고스트 표적에 해당하는 피크를 구별한 후, 배열된 위상 값을 이용하여 표적의 올바른 속도를 추정하고 각도 추정 또한 동시에 수행할 수 있다.In the range-speed detection results of the BPSK-MIMO radar system, each target appears overlapping as many as the number of transmitting antennas. That is, one real target and several ghost targets are detected. In order to distinguish between real targets and ghost targets in the received signals of the BPSK-MIMO radar system, the idea is used that the phase values of the peak signals corresponding to the real targets and ghost targets must have the form of an arithmetic sequence. Therefore, by properly arranging the phase values of these peaks in the received signal into an arithmetic sequence, it is possible to determine which transmitting antenna each target was generated from. That is, after distinguishing between the peak corresponding to the actual target and the peak corresponding to the ghost target in the received signal, the correct speed of the target can be estimated using the arranged phase values, and angle estimation can also be performed at the same time.

이와 같은 본 발명의 예시적인 실시예들에 따르면, 추가적인 하드웨어의 설치 없이 기존 레이더 시스템의 디지털 신호 프로세서에서 구현될 수 있다. 신호 처리 과정을 통해 표적의 속도 및 각도 추정 성능을 개선할 수 있다. 즉, 기존의 디지털 신호 처리기에 간단히 구현할 수 있다. 거리-속도 추정을 위한 2차원(2D) 고속 푸리에 변환(Fast Fourier Transform: FFT)은 FMCW 레이더 시스템에서 수행해야 하는 필수 신호 처리이며, 위상 값을 정렬하는 알고리즘만 추가로 구현하면 되기 때문이다. According to exemplary embodiments of the present invention, it can be implemented in a digital signal processor of an existing radar system without installing additional hardware. The target's speed and angle estimation performance can be improved through signal processing. In other words, it can be simply implemented in an existing digital signal processor. Two-dimensional (2D) Fast Fourier Transform (FFT) for distance-velocity estimation is an essential signal processing that must be performed in an FMCW radar system, and only an algorithm for sorting phase values needs to be additionally implemented.

본 발명은 복수 개의 송신 안테나를 사용했을 때 발생하는 최대 탐지 속도 저하 문제를 해결할 수 있다. The present invention can solve the problem of lower maximum detection speed that occurs when using multiple transmit antennas.

본 발명은 또한, 각 수신 안테나에 동시에 수신된 신호들을 사용하기 때문에 표적의 움직임으로 인해 발생하는 각도 추정 성능 저하 문제를 해결할 수 있다. 즉, 본 발명은 신호를 교대로 수신하는 TDM 방식과 달리 동시에 수신되는 신호를 사용하므로, 서로 다른 시간 슬롯 사이에서 대상의 이동으로 인한 각도 추정의 부정확성 문제를 해결할 수 있다. The present invention can also solve the problem of angle estimation performance degradation caused by target movement because it uses signals received simultaneously at each receiving antenna. That is, the present invention uses signals received simultaneously, unlike the TDM method in which signals are received alternately, and thus can solve the problem of angle estimation inaccuracy due to movement of the object between different time slots.

도 1은 본 발명의 예시적인 실시예에 따른 BPSK-MIMO FMCW 레이더 시스템의 구성을 나타내는 블록도이다.
도 2는 FMCW 레이더 시스템의 송신 신호 파형을 예시한다.
도 3은 ULA(Uniform Linear Array) 안테나 시스템을 이용한 표적에 대한 각도 추정 방법을 도식적으로 나타낸다.
도 4는 MIMO ULA 안테나 시스템을 사용한 각도 추정을 예시한다.
도 5는 2개의 송신 안테나를 사용하는 BPSK-MIMO 레이더 시스템의 송신 신호를 보여준다.
도 6은 2개의 송신 안테나를 구비한 BPSK-MIMO 레이더 시스템을 이용한 시뮬레이션에서 거리-속도 추정 결과를 보여준다.
도 7은 2개의 송신 안테나 대신 4개의 송신 안테나를 구비한 BPSK-MIMO 레이더 시스템을 이용한 시뮬레이션의 거리-속도 추정 결과를 보여준다.
도 8은 정지 표적(A)과 비정지 표적(A)의 경우에 대한 각도 추정 결과를 보여준다.
도 9는 본 발명의 예시적인 실시예에 따른 수신신호의 송신안테나별 구분 방법을 적용하는 전체적인 절차를 나타낸다.
도 10은 본 발명의 예시적인 실시예에 따라 프로세서를 이용한 신호처리를 통해 수신신호들을 송신안테나별로 구분하여 표적의 속도 및 각도를 추정하는 과정을 도시한다.
도 11은 본 발명의 예시적인 실시예에 따라 2개의 송신 안테나와 4개의 수신 안테나를 사용하여 시뮬레이션 한 거리-속도 추정 결과를 보여준다: (A)는 3차원도이고, (B)는 위에서 내려본 2차원도이다.
도 12는 위상 시퀀스 [Pa Pb] 및 [Pd Pc]에 해당하는 8개의 가상 안테나의 래핑되지 않은(unwrapped) 위상 값을 보여준다.
도 13은 등차수열로 배열된 위상 시퀀스에 FFT를 적용하여 제1 및 제2 표적의 각도를 추정한 결과를 보여준다.
도 14는 다양한 SNR에 대한 8개의 가상 안테나의 래핑되지 않은 위상 값을 보여준다.
도 15는 3개의 송신 안테나를 사용하는 3-PSK MIMO 레이더 시스템의 전송 신호를 보여준다.
도 16은 본 발명의 예시적인 실시예에 따라 3개의 송신 안테나와 4개의 수신 안테나를 갖는 3-PSK MIMO 레이더 시스템을 사용하여 시뮬레이션 한 거리-속도 추정 결과를 보여준다: (A)는 3차원도이고, (B)는 위에서 내려본 2차원도이다.
1 is a block diagram showing the configuration of a BPSK-MIMO FMCW radar system according to an exemplary embodiment of the present invention.
Figure 2 illustrates the transmitted signal waveform of the FMCW radar system.
Figure 3 schematically shows an angle estimation method for a target using a Uniform Linear Array (ULA) antenna system.
Figure 4 illustrates angle estimation using a MIMO ULA antenna system.
Figure 5 shows the transmitted signal of a BPSK-MIMO radar system using two transmit antennas.
Figure 6 shows the distance-speed estimation results from a simulation using a BPSK-MIMO radar system equipped with two transmit antennas.
Figure 7 shows the distance-speed estimation results of a simulation using a BPSK-MIMO radar system with four transmit antennas instead of two transmit antennas.
Figure 8 shows angle estimation results for the case of a stationary target (A) and a non-stationary target (A).
Figure 9 shows the overall procedure of applying the method for classifying received signals by transmission antenna according to an exemplary embodiment of the present invention.
Figure 10 shows a process for estimating the speed and angle of a target by dividing received signals by transmission antenna through signal processing using a processor according to an exemplary embodiment of the present invention.
Figure 11 shows simulated distance-velocity estimation results using two transmit antennas and four receive antennas according to an exemplary embodiment of the present invention: (A) is a three-dimensional view, and (B) is a top-down view. It is a two-dimensional diagram.
Figure 12 shows the unwrapped phase values of eight virtual antennas corresponding to the phase sequences [P a P b ] and [P d P c ].
Figure 13 shows the results of estimating the angles of the first and second targets by applying FFT to the phase sequence arranged in an arithmetic sequence.
Figure 14 shows the unwrapped phase values of eight virtual antennas for various SNRs.
Figure 15 shows the transmission signal of a 3-PSK MIMO radar system using three transmit antennas.
Figure 16 shows range-velocity estimation results simulated using a 3-PSK MIMO radar system with three transmit antennas and four receive antennas according to an exemplary embodiment of the present invention: (A) is a three-dimensional diagram; , (B) is a two-dimensional diagram viewed from above.

이하, 첨부한 도면들을 참조하여, 본 발명의 바람직한 실시예를 보다 상세하게 설명하고자 한다. 도면상의 동일한 구성요소에 대해서는 동일한 참조부호를 사용하고 동일한 구성요소에 대해서 중복된 설명은 생략한다.Hereinafter, preferred embodiments of the present invention will be described in more detail with reference to the attached drawings. The same reference numerals are used for the same components in the drawings, and duplicate descriptions for the same components are omitted.

본 발명에서 사용한 용어는 단지 특정한 실시예를 설명하기 위해 사용된 것으로, 본 발명을 한정하려는 의도가 아니다. 단수의 표현은 문맥상 명백하게 다르게 뜻하지 않는 한, 복수의 표현을 포함한다. 본 출원에서, "포함하다" 또는 "가지다" 등의 용어는 명세서 상에 기재된 특징, 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것이 존재함을 지정하려는 것이지, 하나 또는 그 이상의 다른 특징들이나 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것들의 존재 또는 부가 가능성을 미리 배제하지 않는 것으로 이해되어야 한다. 또한, 제1, 제2 등의 용어는 다양한 구성요소들을 설명하는 데 사용될 수 있지만, 상기 구성요소들은 상기 용어들에 의해 한정되어서는 안 된다. 상기 용어들은 하나의 구성요소를 다른 구성요소로부터 구별하는 목적으로만 사용된다. The terms used in the present invention are only used to describe specific embodiments and are not intended to limit the present invention. Singular expressions include plural expressions unless the context clearly dictates otherwise. In this application, terms such as “comprise” or “have” are intended to designate the presence of features, numbers, steps, operations, components, parts, or combinations thereof described in the specification, but are not intended to indicate the presence of one or more other features. It should be understood that this does not exclude in advance the possibility of the existence or addition of elements, numbers, steps, operations, components, parts, or combinations thereof. Additionally, terms such as first, second, etc. may be used to describe various components, but the components should not be limited by the terms. The above terms are used only for the purpose of distinguishing one component from another.

도 1은 본 발명의 예시적인 실시예에 따른 BPSK-MIMO FMCW 레이더 시스템의 구성을 나타내는 블록도이다.1 is a block diagram showing the configuration of a BPSK-MIMO FMCW radar system according to an exemplary embodiment of the present invention.

도 1을 참조하면, BPSK-MIMO FMCW 레이더 시스템(10)은 송신부(20)와 이에 연결된 송신 안테나(30), 수신부(40)와 이에 연결된 수신 안테나(50), 그리고 디지털 신호 처리부(60), 거리/속도 추정부(70), 각도 추정부(80)를 포함할 수 있다.Referring to FIG. 1, the BPSK-MIMO FMCW radar system 10 includes a transmitter 20 and a transmit antenna 30 connected thereto, a receiver 40 and a receive antenna 50 connected thereto, and a digital signal processor 60, It may include a distance/speed estimator 70 and an angle estimator 80.

송신부(20)는 신호 발생기(22)와 발진기(24)를 포함할 수 있다. 신호 발생기(22)와 발진기(24)는 시간에 따라 주파수가 증가하는 처프 신호를 생성하여 송신 안테나(30)를 통해 무선 송신하도록 구성될 수 있다. 구체적으로 파형 발생기(22)는 제공되는 송신정보(디지털)에 기초하여 원하는 주기와 모양의 아날로그 신호를 생성할 수 있다. 예를 들어 신호발생기(22)는 삼각파형의 변조된 신호(삼각파)를 발진기(24)에 제공할 수 있다. 발진기(24)는 파형발생기(24)가 제공하는 신호를 무선 송출하기 위해 시간에 따라 주파수가 증가하는 처프 신호로 변환할 수 있다. 그 처프 신호는 무선송출에 필요한 출력으로 증폭되어 송신 안테나(30)를 통해 송신될 수 있다. 발진기(24)는 예컨대 전압제어발진기(Voltage Control Oscillator: VCO)로 구성될 수 있다. The transmitter 20 may include a signal generator 22 and an oscillator 24. The signal generator 22 and the oscillator 24 may be configured to generate a chirp signal whose frequency increases with time and transmit it wirelessly through the transmission antenna 30. Specifically, the waveform generator 22 can generate an analog signal with a desired period and shape based on the provided transmission information (digital). For example, the signal generator 22 may provide a modulated signal (triangular wave) in a triangular waveform to the oscillator 24. The oscillator 24 can convert the signal provided by the waveform generator 24 into a chirp signal whose frequency increases with time for wireless transmission. The chirp signal can be amplified to the output required for wireless transmission and transmitted through the transmitting antenna 30. The oscillator 24 may be configured as, for example, a voltage control oscillator (VCO).

이렇게 송신된 신호가 주변의 표적에 의해 반사되어 수신 안테나(50)를 통해 수신될 수 있다. 수신부(40)의 구성요소들은 수신 안테나(50)에 연결된 cos 채널(I 채널)과 sin 채널(Q 채널)을 형성할 수 있다. cos 채널(I 채널)은 제1 믹서(41), 제1 저역통과필터(LPF, 44), 제1 아날로그-디지털 변환기(ADC, 46)를 포함할 수 있고, sin 채널(Q 채널)도 마찬가지로 제2 믹서(43), 제2 저역통과필터(45), 제2 아날로그-디지털 변환기(ADC, 47)를 포함할 수 있다. cos 채널(I 채널)에서, 제1 믹서(41)는 수신 안테나(50)와 발진기(24)에 연결되어 수신 신호를 송신에 사용한 처프 신호와 곱해준다. 제1 믹서(41)에서 두 개의 정현파 신호를 곱하면 고주파 신호와 저주파 신호가 생성되는데, 제1 LPF(44)를 통해 고주파 신호를 제거해줄 수 있다. 제1 ADC(46)는 제1 LPF(44)의 출력신호의 샘플링을 통해 아날로그 신호를 디지털 신호로 변환할 수 있다. sin 채널(Q 채널)에서도 cos 채널(I 채널)에서와 동일한 신호처리가 이루어질 수 있다. 다만, 제2 믹서(43)에는, 송신에 사용한 처프 신호(즉, 발진기(24) 출력신호)가 그대로 제공되는 것이 아니라, 위상 쉬프터(42)를 통해 90° 만큼 위상 쉬프트한 처프 신호가 제공된다는 점이 cos 채널(I 채널)에서의 신호처리와 다르다. 이처럼, 제1 믹서(41)는 수신신호에 송신용 처프 신호를 그대로 곱해주고, 또 다른 제2 믹서(43)는 수신신호에 처프 신호를 90도만큼 위상 쉬프트 한 뒤 곱해줌으로써, cos 채널(I 채널) 신호와 sin 채널(Q 채널) 신호를 형성할 수 있다.The signal transmitted in this way may be reflected by a surrounding target and received through the receiving antenna 50. Components of the receiving unit 40 may form a cos channel (I channel) and a sin channel (Q channel) connected to the receiving antenna 50. The cos channel (I channel) may include a first mixer (41), a first low-pass filter (LPF, 44), and a first analog-to-digital converter (ADC, 46), and the sin channel (Q channel) may also include It may include a second mixer 43, a second low-pass filter 45, and a second analog-to-digital converter (ADC, 47). In the cos channel (I channel), the first mixer 41 is connected to the receiving antenna 50 and the oscillator 24 and multiplies the received signal with the chirp signal used for transmission. When two sinusoidal signals are multiplied in the first mixer 41, a high-frequency signal and a low-frequency signal are generated, and the high-frequency signal can be removed through the first LPF 44. The first ADC 46 can convert an analog signal into a digital signal through sampling of the output signal of the first LPF 44. The same signal processing as in the cos channel (I channel) can be performed in the sin channel (Q channel). However, the chirp signal used for transmission (i.e., the output signal of the oscillator 24) is not provided as is to the second mixer 43, but a chirp signal phase-shifted by 90° through the phase shifter 42 is provided. This is different from signal processing in the cos channel (I channel). In this way, the first mixer 41 multiplies the received signal by the chirp signal for transmission, and the second mixer 43 phases the received signal by shifting the chirp signal by 90 degrees and then multiplies it, thereby producing a cos channel (I channel) signal and sin channel (Q channel) signal can be formed.

디지털 신호 처리부(60)는 수신부(40)의 제1 및 제2 ADC(46, 47)를 통과한 신호에 여러 디지털 신호처리(윈도우잉, 필터링, FFT 등)를 수행할 수 있다. 특히, 디지털화 된 cos 채널(I 채널) 신호와 sin 채널(Q 채널) 신호에 대하여 2D FFT를 수행할 수 있다.The digital signal processing unit 60 may perform various digital signal processing (windowing, filtering, FFT, etc.) on the signal that has passed through the first and second ADCs 46 and 47 of the receiving unit 40. In particular, 2D FFT can be performed on digitized cos channel (I channel) signals and sin channel (Q channel) signals.

디지털 신호 처리부(60)에서 처리된 디지털 신호는 거리/속도 추정부(70)에 제공되어 표적의 거리와 속도를 추정할 수 있고, 또한 각도 추정부(80)에 제공되어 표적의 각도를 추정할 수 있다. 특히, 거리/속도 추정부(70)는 2D FFT 처리된 신호에 대하여 CFAR 알고리즘을 적용하여 거리-속도 도메인에서의 피크들을 검출하고 각 피크의 위상 값을 추출하는 처리, 추출된 위상 값들이 등차수열을 이루도록 배열하고, 그렇게 등차수열로 배열된 위상 값 시퀀스에 기반하여 수신신호들과 각 송신 안테나들 간의 대응관계를 결정하여 수신신호를 송신안테나별로 구분하는 처리, 그리고 각 송신안테나별로 구분된 수신신호를 이용하여 표적의 거리와 속도를 추정하는 처리를 수행할 수 있다. 각도 추정부(80) 역시 각 송신안테나별로 구분된 수신신호를 이용하여 각도 추정 알고리즘을 적용하여 각 표적의 각도정보를 산출할 수 있다.The digital signal processed in the digital signal processing unit 60 is provided to the distance/speed estimator 70 to estimate the distance and speed of the target, and is also provided to the angle estimator 80 to estimate the angle of the target. You can. In particular, the distance/speed estimation unit 70 detects peaks in the distance-speed domain by applying the CFAR algorithm to the 2D FFT processed signal and extracts the phase value of each peak, and the extracted phase values are converted into an arithmetic sequence. Arranged to form an arithmetic sequence, determine the correspondence relationship between the received signals and each transmitting antenna based on the phase value sequence arranged in an arithmetic sequence, process to classify the received signals for each transmitting antenna, and receive signals classified for each transmitting antenna. You can use to perform processing to estimate the distance and speed of the target. The angle estimation unit 80 can also calculate the angle information of each target by applying an angle estimation algorithm using received signals classified for each transmission antenna.

디지털 신호 처리부(60), 거리/속도 추정부(70), 그리고 각도 추정부(80)는 처리장치(90)와 같은 하드웨어와 그 처리장치(90)에서 실행될 수 있는 프로그램(이 프로그램은 후술하는 본 발명의 실시예에 따른 방법의 알고리즘을 컴퓨터 프로그램으로 구현한 것임)으로 구현될 수 있다. 처리장치(90)는 예를 들어 프로세서, 콘트롤러, ALU(arithmetic logic unit), 디지털 신호 프로세서(digital signal processor), 마이크로컴퓨터, FPA(field programmable array), PLU(programmable logic unit), 마이크로프로세서, 또는 명령(instruction)을 실행하고 응답할 수 있는 다른 어떠한 장치와 같이, 하나 이상의 범용 컴퓨터 또는 특수 목적 컴퓨터를 이용하여 구현될 수 있다. 처리 장치(90)는 운영 체제(OS) 및 상기 운영 체제 상에서 수행되는 하나 이상의 소프트웨어 애플리케이션을 수행할 수 있다. 또한, 처리 장치는 소프트웨어의 실행에 응답하여, 데이터를 접근, 저장, 조작, 처리 및 생성할 수도 있다.The digital signal processing unit 60, the distance/speed estimating unit 70, and the angle estimating unit 80 include hardware such as the processing unit 90 and a program that can be executed on the processing unit 90 (this program will be described later). The algorithm of the method according to the embodiment of the present invention is implemented as a computer program. The processing unit 90 may include, for example, a processor, a controller, an arithmetic logic unit (ALU), a digital signal processor, a microcomputer, a field programmable array (FPA), a programmable logic unit (PLU), a microprocessor, or Like any other device capable of executing and responding to instructions, it may be implemented using one or more general-purpose or special-purpose computers. Processing device 90 may execute an operating system (OS) and one or more software applications running on the operating system. Additionally, a processing device may access, store, manipulate, process, and generate data in response to the execution of software.

다음으로, BPSK-MIMO FMCW 레이더 시스템(10)의 기본적인 동작 원리를 좀 더 구체적으로 설명한다. FMCW 레이더는 각 주파수가 시간에 따라 선형적으로 증가하는 주기적 처프 신호를 전송한다. 단일 처프에 대해 전송된 신호는 다음과 같이 표현될 수 있다.Next, the basic operating principle of the BPSK-MIMO FMCW radar system 10 will be described in more detail. FMCW radar transmits periodic chirp signals where each frequency increases linearly with time. The transmitted signal for a single chirp can be expressed as:

......(1) ......(One)

여기서 fc, B 및 Tc는 각각 처프 신호의 캐리어 주파수, 대역폭 및 지속 시간을 각각 나타낸다. where f c , B and T c represent the carrier frequency, bandwidth and duration of the chirp signal, respectively.

도 2는 FMCW 레이더 시스템의 송신 신호 파형을 나타낸다. 도 2에 도시된 바와 같이, 송신부(20)에 의한 송신은 L개의 연속적인 처프에 대해 반복된다. 그런 다음, 표적에 의해 반사되어 수신부(40)에 수신된 신호는 전송된 신호와 곱해져서 반송파 주파수 성분을 제거한다. 그 결과로 생성된 중간 주파수(IF) 대역 신호는 다음과 같이 2D 지수 함수로 표현될 수 있다.Figure 2 shows the transmission signal waveform of the FMCW radar system. As shown in Figure 2, transmission by the transmitter 20 is repeated for L consecutive chirps. Then, the signal reflected by the target and received by the receiver 40 is multiplied by the transmitted signal to remove the carrier frequency component. The resulting intermediate frequency (IF) band signal can be expressed as a 2D exponential function as follows.

......(2) ......(2)

여기서, n과 N은 각 처프 내 샘플링 인덱스와 샘플링 개수를 나타내고, l과 L은 각 처프 내 처프 인덱스와 처프 개수를 나타낸다. 또한, R과 v는 표적의 거리와 속도를 나타내고, Ts는 샘플링 주기, c는 빛의 속도, φ는 고정 위상 성분을 나타낸다. 이 IF 대역 신호는 n-축과 l-축을 따라 각각 2BR/Tcc와 2v/λ의 주파수를 갖는다. 따라서 식 (2)에 2D 푸리에 변환을 적용하면 표적의 거리와 속도를 추정할 수 있다.Here, n and N represent the sampling index and sampling number within each chirp, and l and L represent the chirp index and number of chirps within each chirp. Additionally, R and v represent the distance and speed of the target, T s represents the sampling period, c represents the speed of light, and ϕ represents the fixed phase component. This IF band signal has frequencies of 2BR/T c c and 2v/λ along the n- and l-axes, respectively. Therefore, by applying the 2D Fourier transform to equation (2), the distance and speed of the target can be estimated.

또한, 각도 추정은 균일 선형 어레이(Uniform Linear Array: ULA) 안테나 시스템을 사용하여 달성할 수 있다. 도 3은 ULA 안테나 시스템(송신 안테나 1개, 수신 안테나 8개)을 이용하여 표적의 각도 추정을 하는 것을 나타낸다. 도 3에 예시된 ULA 안테나 시스템(100)에서, 다수의 수신 안테나들(Rx1 - Rx8)의 안테나 간격이 d이고, 하나의 송신 안테나(Tx1)에서 송신되어 표적에 의해 반사되어 수신 안테나(Rx1 - Rx8)를 통해 수신되는 수신신호의 파장과 입사각이 각각 λ과 θ인 경우, 인접 수신 안테나 간의 위상차(ω)는 ω= 2π(d sinθ)/λ과 같이 나타낼 수 있다. 이 위상 항 ω은 식 (2)에 통합되어 다음과 같이 표현될 수 있는 3차원 IF 대역 신호를 만들 수 있다. Additionally, angle estimation can be achieved using a Uniform Linear Array (ULA) antenna system. Figure 3 shows the angle estimation of the target using the ULA antenna system (1 transmitting antenna, 8 receiving antennas). In the ULA antenna system 100 illustrated in FIG. 3, the antenna spacing of the plurality of receiving antennas (Rx1 - Rx8) is d, and is transmitted from one transmitting antenna (Tx1) and reflected by the target to receive the receiving antennas (Rx1 - When the wavelength and incidence angle of the received signal received through Rx8) are λ and θ, respectively, the phase difference (ω) between adjacent receiving antennas can be expressed as ω= 2π(d sinθ)/λ. This phase term ω can be integrated into equation (2) to create a three-dimensional IF band signal that can be expressed as:

......(3) ......(3)

여기서, p는 안테나 인덱스이고, P는 수신 안테나(Rx1 - Rx8)의 수이다. 이 신호는 일반적으로 레이더 데이터 큐브(radar data cube)로 알려져 있다. 레이더 데이터 큐브의 각 차원에 푸리에 변환을 적용하여 표적의 거리, 속도 및 각도를 추정할 수 있다.Here, p is the antenna index, and P is the number of receiving antennas (Rx1 - Rx8). This signal is commonly known as a radar data cube. By applying the Fourier transform to each dimension of the radar data cube, the target's range, speed, and angle can be estimated.

도 4는 MIMO ULA 안테나 시스템을 사용한 각도 추정을 예시한다. MIMO 레이더 시스템에서는 레이더의 탐지 성능을 향상시키기 위해 ULA 형태로 배열된 여러 개의 송신안테나 및 수신 안테나가 사용될 수 있다. 예를 들어, 도 4는 2개의 송신 안테나(Tx1, Tx2)와 4개의 수신 안테나(Rx1 - Rx4)가 있는 MIMO 레이더 시스템의 MIMO ULA 안테나 구성을 보여준다. 제2 송신 안테나(Tx2)에서 전송된 신호는 제1 송신 안테나(Tx1)에서 전송된 신호에 비해 4d sinθ의 거리를 더 이동하므로, 제2 송신 안테나(Tx2)에서 전송된 수신 신호의 위상은 제1 송신 안테나(Tx1)에서 전송된 수신신호의 위상보다 2π(4d sinθ)/λ만큼 더 높다. 결과적으로 도 4에 도시된 안테나 구성의 각도 추정 성능은 도 3에 도시된 안테나 구성의 그것과 유사하다. MIMO 안테나 방식의 주요 장점은 전체 안테나 소자의 수를 줄일 수 있다는 점이다.Figure 4 illustrates angle estimation using a MIMO ULA antenna system. In a MIMO radar system, multiple transmitting and receiving antennas arranged in a ULA form can be used to improve radar detection performance. For example, Figure 4 shows the MIMO ULA antenna configuration of a MIMO radar system with two transmit antennas (Tx1, Tx2) and four receive antennas (Rx1 - Rx4). Since the signal transmitted from the second transmitting antenna (Tx2) moves a distance of 4d sinθ compared to the signal transmitted from the first transmitting antenna (Tx1), the phase of the received signal transmitted from the second transmitting antenna (Tx2) is 1 The phase of the received signal transmitted from the transmitting antenna (Tx1) is higher by 2π(4d sinθ)/λ. As a result, the angle estimation performance of the antenna configuration shown in FIG. 4 is similar to that of the antenna configuration shown in FIG. 3. The main advantage of the MIMO antenna method is that the overall number of antenna elements can be reduced.

MIMO 안테나 시스템을 최대한 활용하려면 서로 다른 송신 안테나에서 전송되는 신호를 구별해야 한다. 예를 들어, 도 4의 경우 제1 송신 안테나(Tx1)에서 어떤 신호 성분이 전송되었고, 제2 송신 안테나(Tx2)에서는 어떤 신호 성분이 전송되는지를 수신신호에서 구별할 필요가 있다. 이러한 문제를 해결하기 위해 TDM, FDM, CDM 등 다양한 다중화 방식이 존재한다. 이러한 다양한 기술 중 BPSK는 미리 할당된 코드를 이용하여 송신 신호의 위상을 변조하는 방식이다. 이 코드는 전송된 신호들이 서로 직교하도록 특별히 설계된다.To get the most out of a MIMO antenna system, signals transmitted from different transmit antennas must be distinguished. For example, in the case of FIG. 4, it is necessary to distinguish in the received signal which signal component is transmitted from the first transmission antenna (Tx1) and which signal component is transmitted from the second transmission antenna (Tx2). To solve this problem, various multiplexing methods such as TDM, FDM, and CDM exist. Among these various technologies, BPSK is a method of modulating the phase of a transmission signal using a pre-assigned code. This code is specifically designed so that the transmitted signals are orthogonal to each other.

도 5는 2개의 송신 안테나를 사용하는 BPSK-MIMO 레이더 시스템의 송신 신호를 보여준다. (A)와 (B)의 신호는 제1 송신 안테나와 제2 송신 안테나를 통해 각각 송신될 수 있다. 처프(chirp)에 코드 1이 할당되면, 전송 신호의 위상이 0으로 변조되는데, 이는 아래 식 (4)로 표시되는 신호가 전송됨을 의미한다. Figure 5 shows the transmitted signal of a BPSK-MIMO radar system using two transmit antennas. Signals (A) and (B) may be transmitted through the first and second transmission antennas, respectively. When code 1 is assigned to the chirp, the phase of the transmission signal is modulated to 0, which means that the signal represented by equation (4) below is transmitted.

......(4) ......(4)

반면에 처프에 0의 코드가 할당되면, 전송신호의 위상은 π만큼 변조되고 아래 식 (5)의 신호가 전송된다. On the other hand, if a code of 0 is assigned to the chirp, the phase of the transmission signal is modulated by π and the signal in equation (5) below is transmitted.

......(5) ......(5)

도 5와 같이 제1 송신 안테나(Tx1)가 1111의 코드로 송신하고, 제2 송신 안테나(Tx2)가 1010의 코드로 송신할 때, 홀수 인덱스의 수신 신호는 rodd(t) = r1(t) + r2(t)로 표현될 수 있고, 짝수 인덱스의 수신 신호는 reven(t) = r1(t) - r2(t)로 표현될 수 있다. 여기서, r1(t)는 제1 송신 안테나(Tx1)로부터 수신된 신호 성분을 나타내고, 와 r2(t)는 제2 송신 안테나(Tx2)로부터 수신된 신호 성분을 나타낸다. 따라서 각 송신 안테나로부터의 수신 신호는 다음과 같은 디코딩 과정을 통해 분해될 수 있다.As shown in FIG. 5, when the first transmission antenna (Tx1) transmits with a code of 1111 and the second transmission antenna (Tx2) transmits with a code of 1010, the received signal with an odd index is r odd (t) = r 1 ( It can be expressed as t) + r 2 (t), and the received signal of an even index can be expressed as r even (t) = r 1 (t) - r 2 (t). Here, r 1 (t) represents the signal component received from the first transmission antenna (Tx1), and r 2 (t) represents the signal component received from the second transmission antenna (Tx2). Therefore, the received signal from each transmit antenna can be decomposed through the following decoding process.

......(6) ......(6)

......(7) ......(7)

한편, BPSK-MIMO 레이더 시스템을 사용하여 표적을 탐지하는 것은 주로 표적의 거리(range), 속도(velocity), 각도(angle)를 추정하는 것을 포함한다. BPSK-MIMO 레이더 시스템에서 거리 추정은 일반적인 FMCW 레이더 시스템과 유사하게 수행된다. 그러나 기존의 속도 및 각도 추정 방법을 사용하는 경우 MIMO 시스템의 도플러-각도 모호성으로 인해 문제가 발생한다. Meanwhile, detecting a target using a BPSK-MIMO radar system mainly involves estimating the target's range, velocity, and angle. In the BPSK-MIMO radar system, distance estimation is performed similarly to a typical FMCW radar system. However, when using traditional velocity and angle estimation methods, problems arise due to the Doppler-angle ambiguity of MIMO systems.

FMCW 레이더 시스템에서 표적의 속도 추정은 여러 처프 신호 간의 위상차를 사용하여 수행된다. 예를 들어, 기존의 FMCW 레이더 시스템에서는 식 (2)와 식 (3)에서 관찰된 바와 같이 인접한 처프들 간에 위상이 2π(2v/λ)Tc 만큼 차이가 나고, 속도는 다중 처프 신호들 간의 위상 증분을 이용하여 추정된다. 그러나 BPSK-MIMO 레이더 시스템에서 수신된 신호는 다중 송신 안테나로부터의 다중 신호 성분으로 구성되며, 위상 증가분은 단순히 2π(2v/λ)Tc 로 표현될 수 없다. ........(8)In FMCW radar systems, target velocity estimation is performed using the phase difference between several chirp signals. For example, in the existing FMCW radar system, as observed in Equations (2) and (3), the phase between adjacent chirps differs by 2π(2v/λ)T c , and the speed is the difference between multiple chirp signals. It is estimated using phase increments. However, the signal received in the BPSK-MIMO radar system consists of multiple signal components from multiple transmit antennas, and the phase increment cannot be simply expressed as 2π(2v/λ)T c . ........(8)

식 (8)은 도 4와 같은 안테나 구성을 설정하고, 도 5와 같이 송신 신호를 변조한 경우, 제1 수신 안테나(Rx1)의 수신 신호를 나타낸다. 식 (8)에서, 첫 번째 항은 2π(2v/λ)Tc의 위상 증분을 갖는 일반적인 수신 신호를 나타낸다. 그러나 두 번째 항은 부호 반전을 포함하므로 위상이 2π(2v/λ)Tc+π 만큼씩 증가한다. 그 결과 IF대역 신호에 2D 푸리에 변환을 적용하면 거리-속도 검출 결과에서 2 개의 피크가 발생한다.Equation (8) represents the received signal of the first receiving antenna (Rx1) when the antenna configuration as shown in FIG. 4 is set and the transmission signal is modulated as shown in FIG. 5. In equation (8), the first term represents a typical received signal with a phase increment of 2π(2v/λ)T c . However, the second term includes a sign inversion, so the phase increases by 2π(2v/λ)T c +π. As a result, when 2D Fourier transform is applied to the IF band signal, two peaks occur in the distance-speed detection result.

레이더의 시야에서 3개의 표적이 움직인다고 가정하고, 앞서 언급한 안테나 구성을 사용하여 시뮬레이션을 수행했다. 도 6은 2개의 송신 안테나를 구비한 BPSK-MIMO 레이더 시스템을 이용한 시뮬레이션에서 거리-속도 추정 결과를 보여준다. 이 시뮬레이션을 통해 얻은 거리-속도 추정 결과를 보여준다. 시뮬레이션 파라미터는 기존의 자동차 레이더 파라미터를 사용하는 표 1과 같이 설정하였다. 그 결과 ㅁλ/(4Tc) = ㅁ19.48 m/s의 최대 탐지가능 속도(maximum detectable velocity)가 나타났다. 또한 SNR은 0dB로 설정하였고, 표적들의 거리와 속도는 (30m, 12m/s), (50m, 4m/s), (60m, -3m/s)로 각각 설정하였다. Simulations were performed using the previously mentioned antenna configuration, assuming three targets moving in the radar's field of view. Figure 6 shows the distance-speed estimation results from a simulation using a BPSK-MIMO radar system equipped with two transmit antennas. It shows the distance-speed estimation results obtained through this simulation. The simulation parameters were set as shown in Table 1 using existing automotive radar parameters. As a result, the maximum detectable velocity of ㅁλ/(4T c ) = ㅁ19.48 m/s was obtained. Additionally, the SNR was set to 0dB, and the distance and speed of the targets were set to (30m, 12m/s), (50m, 4m/s), and (60m, -3m/s), respectively.

도 6에서 볼 수 있듯이 각 표적에 대해 거리-속도 검출 결과에서 두 개의 피크가 발생하여 도플러 모호성이 발생함을 알 수 있다. 이는 하나의 피크는 제1 송신 안테나(Tx1)의 신호로 인해 생성되고, 나머지 하나의 피크는 제2 송신 안테나(Tx2)의 신호로 인해 생성되기 때문이다. 예를 들어 30m 떨어진 지점에 위치한 표적의 경우, 제1 송신 안테나(Tx1)로부터의 송신신호로 인해 12m/s에서 피크가 발생하였고, 다른 하나의 피크는 제2 송신 안테나(Tx2)의 신호로 인해 12-19.48 = -7.48 m/s에서 발생했다. 나머지 다른 표적에 대해서도 유사한 결과가 관찰되었으며, 표적의 속도는 유일하게(uniquely) 결정되지 않았다. 그 결과 최대 탐지가능 속도가, 19.48m/s에서 9.74m/s로, 절반으로 감소했다.As can be seen in Figure 6, two peaks occur in the distance-velocity detection results for each target, resulting in Doppler ambiguity. This is because one peak is generated by the signal of the first transmission antenna (Tx1), and the other peak is generated by the signal of the second transmission antenna (Tx2). For example, in the case of a target located 30m away, a peak occurred at 12m/s due to the transmission signal from the first transmission antenna (Tx1), and another peak occurred due to the signal from the second transmission antenna (Tx2). Occurred at 12-19.48 = -7.48 m/s. Similar results were observed for the remaining targets, whose speeds were not uniquely determined. As a result, the maximum detectable speed was reduced by half, from 19.48 m/s to 9.74 m/s.

앞의 결과는 2개 이상의 송신 안테나로 구성된 BPSK-MIMO 레이더 시스템으로 확장될 수 있다. 도 7은 2개의 송신 안테나 대신 4개의 송신 안테나를 구비한 BPSK-MIMO 레이더 시스템을 이용한 시뮬레이션의 거리-속도 추정 결과를 보여준다. 송신 안테나의 수를 제외하고 레이더 시스템 파라미터와 타겟 정보는 도 6과 동일하게 설정하였다. 또한 4개의 송신 안테나에서 전송되는 신호의 위상은 Hadamard 코딩 방법을 사용하여 코딩하였다. 이 방법은 4개의 송신 안테나가 (1, 1, 1, 1), (1, 0, 1, 0), (1, 1, 0, 0), (1, 0, 0, 1)의 코드로 전송한다. 그리고 각 송신안테나로부터의 수신 신호들은 코드들 사이의 직교성을 이용하여 분해될 수 있다. The previous results can be extended to BPSK-MIMO radar systems consisting of two or more transmit antennas. Figure 7 shows the distance-speed estimation results of a simulation using a BPSK-MIMO radar system with four transmit antennas instead of two transmit antennas. Except for the number of transmitting antennas, the radar system parameters and target information were set the same as in FIG. 6. Additionally, the phases of signals transmitted from the four transmit antennas were coded using the Hadamard coding method. This method uses four transmit antennas with codes of (1, 1, 1, 1), (1, 0, 1, 0), (1, 1, 0, 0), and (1, 0, 0, 1). send. And the received signals from each transmission antenna can be decomposed using orthogonality between codes.

도 7에 도시된 것과 같이 각각의 표적에 대해 거리-속도 검출 결과에서 4개의 피크가 발생하였고, 최대 탐지가능 속도는 4분의 1로 감소하였다. 일반적으로 타임 슬롯을 나누어 다중화를 하는 경우, 송신 안테나의 개수만큼 최대 탐지가능 속도가 감소한다. 이는 속도 추정 성능을 상당히 저하시킨다. 이는 특히 표적에 대하여 매우 빠른 탐지 속도가 요구되는 자동차 애플리케이션에 바람직하지 않을 수 있다. 따라서 도플러 모호성을 해결하고 표적에 대한 최대 탐지가능 속도를 유지하기 위해서는 피크와 송신 안테나 간의 대응 관계를 결정해야 한다.As shown in Figure 7, four peaks occurred in the distance-velocity detection results for each target, and the maximum detectable speed was reduced to one fourth. In general, when multiplexing is performed by dividing time slots, the maximum detectable speed decreases by the number of transmission antennas. This significantly degrades speed estimation performance. This may be undesirable, especially in automotive applications where very fast detection speeds for targets are required. Therefore, in order to resolve Doppler ambiguity and maintain the maximum detectable speed for the target, the correspondence between the peak and the transmit antenna must be determined.

다음으로, BPSK-MIMO 레이다 시스템을 이용한 기존의 각도 추정 방법은 다음과 같다. MIMO 레이더 시스템에서는 앞서 설명한 것처럼 표적의 각도를 추정할 때 각 송신 안테나의 신호를 구분해야 한다. 이를 위해 BPSK-MIMO 레이더 시스템에서 수신된 신호에 식 (6)과 (7)의 디코딩 과정을 적용할 수 있다. 그러나 표적이 비정지(non-stationary) 표적일 때 각 송신 안테나의 신호는 속도 유도 위상 항(Velocity-induced phase term)으로 인해 제대로 구분될 수 없다. 예를 들어, 디코딩 프로세스가 식 (8)에 표시된 수신 신호에 적용될 때 결과 신호는 다음과 같다.Next, the existing angle estimation method using the BPSK-MIMO radar system is as follows. In the MIMO radar system, as described above, the signal from each transmitting antenna must be distinguished when estimating the target angle. To this end, the decoding process of equations (6) and (7) can be applied to the signal received from the BPSK-MIMO radar system. However, when the target is a non-stationary target, the signals from each transmitting antenna cannot be properly distinguished due to the velocity-induced phase term. For example, when the decoding process is applied to the received signal shown in equation (8), the resulting signal is:

......(9) ......(9)

......(10) ......(10)

속도 유도 위상 항 2π(2v/λ)Tc 로 인해, 불필요한 요소가 상쇄되지 못하고 각 송신 안테나로부터의 신호가 제대로 구분될 수 없다. Due to the velocity-induced phase term 2π(2v/λ)T c , unnecessary elements cannot be canceled out and signals from each transmit antenna cannot be properly distinguished.

도 8은 정지 표적(A)과 비정지 표적(A)의 경우에 대한 각도 추정 결과를 보여준다. 도 8을 참조하면, 각도는 -20°로 동일하게 설정하고, 속도는 각각 0m/s 및 15m/s로 설정하였다. 도시된 것과 같이 표적이 정지해 있을 때 표적의 각도가 적절하게 추정되었다. 각도 추정에는 FFT 방법을 사용할 수 있다. FFT 방법을 사용하는 경우, 추정된 각도에 약간의 오차가 있을 수 있다. MUSIC 또는 ESPRIT와 같은 보다 복잡한 알고리즘을 사용하여 추정 정확도를 향상시킬 수 있지만, 빠른 계산을 위해서는 FFT 방법이 유유용하다. 하지만, 표적이 움직이는 경우, 표적의 각도 추정에 큰 오차가 생긴다.Figure 8 shows angle estimation results for the case of a stationary target (A) and a non-stationary target (A). Referring to Figure 8, the angle was set to -20°, and the speed was set to 0 m/s and 15 m/s, respectively. As shown, the angle of the target was appropriately estimated when the target was stationary. The FFT method can be used for angle estimation. When using the FFT method, there may be some error in the estimated angle. Although more complex algorithms such as MUSIC or ESPRIT can be used to improve estimation accuracy, the FFT method is useful for fast calculations. However, when the target moves, a large error occurs in estimating the angle of the target.

이 문제를 해결하기 위해, 수신 신호의 위상을 보상하여 시간 영역 또는 주파수 영역에서 속도 유도 위상 항을 제거할 수 있다. 시간 영역과 주파수 영역에서 목표물의 속도를 먼저 추정한 다음, 추정한 속도 값을 사용하여 수신 신호의 위상을 -(2v/λ)Tc 만큼 보상했다. 그러나 이러한 위상 보상 방법은 추정된 속도 값이 실제 값과 약간 다를 경우에도 큰 오차가 발생할 수 있다. 예를 들어, 속도 분해능이 vres인 경우, 속도 bin은 vres의 간격으로 나뉘며 추정된 속도는 노이즈에 관계없이 실제 값에서vres/2 만큼 벗어날 수 있다. 속도 분해능은 FMCW 레이더 시스템에서 λ/(2LTc)로 정의될 수 있다. 따라서 보상된 위상 값은 만큼 오프셋될 수 있다. 이 값은 와 비슷한 차수이므로, 위상 보상을 적용하더라도 표적의 각도를 확실하게 추정할 수 없다.To solve this problem, the phase of the received signal can be compensated to remove the velocity-induced phase term in the time domain or frequency domain. The target's speed was first estimated in the time domain and frequency domain, and then the phase of the received signal was compensated by -(2v/λ)T c using the estimated speed value. However, this phase compensation method can cause large errors even if the estimated speed value is slightly different from the actual value. For example, if the velocity resolution is v res , the velocity bins are divided into intervals of v res and the estimated velocity may deviate from the actual value by v res /2, regardless of noise. Velocity resolution can be defined as λ/(2LT c ) in FMCW radar systems. Therefore, the compensated phase value is It can be offset by This value is Since it is of a similar order to , the target angle cannot be estimated with certainty even if phase compensation is applied.

본 발명의 예시적인 실시예는 표적의 속도와 각도를 추정할 때 전술한 문제들 특히, 도플러-각도 모호성을 해결하기 위해, 수신된 레이더 신호를 각 송신 안테나별로 효율적으로 구별하는 방법을 제공한다. 이하에서는 그 방법을 구체적으로 설명한다. Exemplary embodiments of the present invention provide a method for efficiently distinguishing received radar signals for each transmit antenna to solve the above-mentioned problems, especially Doppler-angle ambiguity, when estimating the speed and angle of a target. Below, the method is described in detail.

각 송신 안테나의 신호를 구별하기 위해, 발명자들은 수신된 신호의 거리-속도 검출 결과에서 피크의 위상이 도 3과 같이 등차수열을 형성한다는 사실에 주목한다. 예컨대 도 4에 도시된 것처럼, 두 개의 송신 안테나(Tx1, Tx2)와 4개의 수신 안테나(Rx1 - Rx4)가 ULA 형태로 구비된 MIMO FMCW 레이더 시스템에서, 제1 송신 안테나(Tx1)에서 전송되고 제1 수신 안테나(Rx1)에서 수신된 신호의 위상(Pa1) 값이 φ일 때, 이를 포함한 나머지 다른 수신 신호들의 위상(Pa1 - Pa4, Pb1 - Pb4)은 다음과 같이 표현할 수 있다. In order to distinguish the signal of each transmitting antenna, the inventors pay attention to the fact that the phase of the peak in the distance-velocity detection result of the received signal forms an arithmetic sequence as shown in FIG. 3. For example, as shown in FIG. 4, in a MIMO FMCW radar system in which two transmit antennas (Tx1, Tx2) and four receive antennas (Rx1 - Rx4) are provided in the form of a ULA, the first transmit antenna (Tx1) transmits and the second transmits. 1 When the phase (P a1 ) value of the signal received at the receiving antenna (Rx1) is ϕ, the phases (P a1 - P a4 , P b1 - P b4 ) of the remaining received signals including this can be expressed as follows. .

Tx1 → Rx1 : Pa1 = φTx1 → Rx1 : P a1 = ϕ

Tx1 → Rx2 : Pa2 = φ+ωTx1 → Rx2 : P a2 = ϕ+ω

Tx1 → Rx3 : Pa3 = φ+2ωTx1 → Rx3 : P a3 = ϕ+2ω

Tx1 → Rx4 : Pa4 = φ+3ωTx1 → Rx4 : P a4 = ϕ+3ω

Tx2 → Rx1 : Pb1 = φ+4ωTx2 → Rx1 : P b1 = ϕ+4ω

Tx2 → Rx2 : Pb2 = φ+5ωTx2 → Rx2 : P b2 = ϕ+5ω

Tx2 → Rx3 : Pb3 = φ+6ωTx2 → Rx3 : P b3 = ϕ+6ω

Tx2 → Rx4 : Pb4 = φ+7ω ......(11)Tx2 → Rx4 : P b4 = ϕ+7ω ......(11)

이 위상값들은 등차수열을 이룬다. 예컨대, 위의 예에서 두 개의 송신 안테나(Tx1, Tx2)로부터 전송되어 수신된 수신신호의 거리-속도 평면에서의 피크들의 두 위상 시퀀스가 등차수열을 형성할 때, 그 위상 시퀀스 중 선행 서브 위상 시퀀스(former phase sequence)는 제1 송신 안테나(Tx1)로부터 송신된 신호들에 대응하고, 후행 서브 위상 시퀀스(latter phase sequence)는 제2 송신 안테나(Tx2)로부터 송신된 신호들에 대응한다. 여기서, 각 서브 위상 시퀀스는 수신 안테나의 개수와 동일한 수의 위상들의 시퀀스이다. 예컨대, 수신신호 위상 시퀀스 [Pa1 Pa2 Pa3 Pa4 Pb1 Pb2 Pb3 Pb4]가 등차수열을 형성할 때, 선행 위상 시퀀스 [Pa1 Pa2 Pa3 Pa4]는 제1 송신 안테나(Tx1)로부터 송신된 신호들에 대응하고, 후행 위상 시퀀스 [Pb1 Pb2 Pb3 Pb4]는 제2 송신 안테나(Tx2)로부터 송신된 신호들에 대응한다. 만약 위상 시퀀스 [Pb1 Pb2 Pb3 Pb4 Pa1 Pa2 Pa3 Pa4]가 등차수열을 이룬다면, 선행 위상 시퀀스(제1 서브 위상 시퀀스) [Pb1 Pb2 Pb3 Pb4]가 제1 송신 안테나(Tx1)로부터 송신된 신호들에 대응하고, 후행 위상 시퀀스(제2 서브 위상 시퀀스) [Pa1 Pa2 Pa3 Pa4]가 제2 송신 안테나(Tx2)로부터 송신된 신호들에 대응한다. 이러한 등차수열로 배열된 수신신호들의 위상정보는 수신신호들과 송신 안테나들 간의 대응관계를 결정하는 데 활용될 수 있다.These phase values form an arithmetic sequence. For example, in the above example, when the two phase sequences of peaks in the distance-velocity plane of the received signal transmitted and received from the two transmit antennas (Tx1, Tx2) form an arithmetic sequence, the preceding sub-phase sequence among the phase sequences The (former phase sequence) corresponds to signals transmitted from the first transmission antenna (Tx1), and the subsequent sub-phase sequence (latter phase sequence) corresponds to signals transmitted from the second transmission antenna (Tx2). Here, each sub-phase sequence is a sequence of phases equal to the number of receiving antennas. For example, when the received signal phase sequence [P a1 P a2 P a3 P a4 P b1 P b2 P b3 P b4 ] forms an arithmetic sequence, the preceding phase sequence [P a1 P a2 P a3 P a4 ] is the first transmitting antenna. Corresponds to signals transmitted from (Tx1), and the trailing phase sequence [P b1 P b2 P b3 P b4 ] corresponds to signals transmitted from the second transmit antenna (Tx2). If the phase sequence [P b1 P b2 P b3 P b4 P a1 P a2 P a3 P a4 ] forms an arithmetic sequence, then the preceding phase sequence (first sub-phase sequence) [P b1 P b2 P b3 P b4 ] is the first. 1 corresponds to the signals transmitted from the transmission antenna (Tx1), and the trailing phase sequence (second sub-phase sequence) [P a1 P a2 P a3 P a4 ] corresponds to the signals transmitted from the second transmission antenna (Tx2) do. The phase information of the received signals arranged in this arithmetic sequence can be used to determine the correspondence between the received signals and the transmitting antennas.

도 9는 예시적인 실시예에 따른 수신신호의 송신안테나별 구분 방법을 적용하는 전체적인 절차를 나타낸다. 이 방법은 도 1에 도시된 MIMO FMCW 레이더 시스템(10)을 이용하여 수행될 수 있다. Figure 9 shows the overall procedure of applying a method for classifying received signals by transmission antenna according to an exemplary embodiment. This method can be performed using the MIMO FMCW radar system 10 shown in FIG. 1.

도 9를 참조하면, 먼저 MIMO FMCW 레이더 시스템에서 송신부(20)는 송신할 레이더 신호를 생성하여 송신 안테나 어레이(30)를 송신할 수 있다(S100). 송신되는 레이더 신호는 예컨대 PSK 방식으로 변조된 신호일 수 있으며, 대표적인 예로 BPSK 변조방식으로 변조된 신호일 수 있다. Referring to FIG. 9, first, in the MIMO FMCW radar system, the transmitter 20 can generate a radar signal to be transmitted and transmit it to the transmit antenna array 30 (S100). The transmitted radar signal may be, for example, a signal modulated using the PSK method, and as a representative example, it may be a signal modulated using the BPSK modulation method.

송신된 레이더 신호는 시야각 범위 내에 존재하는 표적들에 의해 반사되어 MIMO FMCW 레이더 시스템(10)의 수신 안테나 어레이(50)를 통해 수신될 수 있다. 수신된 신호들은 소정의 신호처리를 거친 다음, 각 수신신호의 위상 정보를 추출하는 처리가 수행될 수 있다(S200). The transmitted radar signal may be reflected by targets existing within the viewing angle range and received through the receiving antenna array 50 of the MIMO FMCW radar system 10. After the received signals undergo predetermined signal processing, processing to extract phase information of each received signal may be performed (S200).

단계 S200에서, MIMO FMCW 레이더 시스템(10)의 수신부(40)는 수신 안테나 어레이(50)로부터 각 수신 신호를 전달받고, 송신부(20)로부터는 전송신호를 제공받는다. cos 채널(I 채널) 신호를 얻기 위해, 제1 믹서(41)에서 수신신호와 전송신호를 곱하고, 그 곱해진 신호에 포함된 고주파 신호를 제1 LPF(44)로 제거한 다음, 제1 ADC(46)를 이용하여 디지털 신호로 변환한다. sin 채널(Q 채널) 신호를 얻기 위해, 송신부(20)로부터 제공되는 전송신호를 90도만큼 위상 천이시킨 후 제2 믹서(43)에서 수신신호와 곱한다. 곱해진 신호에 포함된 고주파 신호를 제2 LPF(45)로 제거한 다음, 제2 ADC(47)를 이용하여 디지털 신호로 변환한다. 이와 같은 신호처리를 통해 수신부(40)에서 cos 채널(I 채널) 신호와 sin 채널(Q 채널) 신호가 얻어져 디지털 신호처리부(60)에 제공될 수 있다. 디지털 신호처리부(60)에서는, 제1 및 제2 ADC(46, 47)를 통과한 신호에 대하여 여러 가지 디지털 신호처리(윈도우잉, 필터링, FFT) 등)를 수행할 수 있다. 이와 같이 신호처리된 수신신호는 거리/속도 추정부(70)에 제공되어, 각 수신신호의 위상 정보가 추출될 수 있다. In step S200, the receiving unit 40 of the MIMO FMCW radar system 10 receives each received signal from the receiving antenna array 50 and receives a transmitted signal from the transmitting unit 20. To obtain a cos channel (I channel) signal, the received signal and the transmitted signal are multiplied in the first mixer 41, the high frequency signal included in the multiplied signal is removed by the first LPF 44, and then the first ADC ( 46) is used to convert it into a digital signal. To obtain a sin channel (Q channel) signal, the transmission signal provided from the transmitter 20 is phase-shifted by 90 degrees and then multiplied by the received signal in the second mixer 43. The high frequency signal included in the multiplied signal is removed by the second LPF (45) and then converted into a digital signal using the second ADC (47). Through this signal processing, a cos channel (I channel) signal and a sin channel (Q channel) signal can be obtained from the receiving unit 40 and provided to the digital signal processing unit 60. The digital signal processing unit 60 can perform various digital signal processing (windowing, filtering, FFT, etc.) on the signals that have passed through the first and second ADCs 46 and 47. The received signal processed in this way is provided to the distance/speed estimator 70, and the phase information of each received signal can be extracted.

도 10은 본 발명의 예시적인 실시예에 따라 처리장치(90)를 이용한 신호처리를 통해 수신신호들을 송신안테나 별로 구분하여 표적의 속도 및 각도를 추정하는 과정을 도시한다.Figure 10 shows a process of estimating the speed and angle of the target by dividing received signals by transmission antenna through signal processing using the processing device 90 according to an exemplary embodiment of the present invention.

도 10을 참조하면, 수신 안테나 어레이(50)를 통해 수신되어 수신부(40)에서 소정의 신호처리를 거친 레이더 수신 신호들의 위상 정보를 추출하기 위해, 먼저 각 수신 안테나에서 수신되는 IF 대역 신호에 2D FFT를 적용할 수 있다(S210). 이 과정을 통해 예컨대 BPSK 변조의 특성으로 인해 각 표적이 복제되는 (즉, 각 표적의 실제 표적 외에 고스트 표적이 검출되는) 거리-속도 검출 결과를 얻을 수 있다. Referring to FIG. 10, in order to extract phase information of radar reception signals received through the reception antenna array 50 and subjected to predetermined signal processing in the reception unit 40, a 2D signal is first applied to the IF band signal received from each reception antenna. FFT can be applied (S210). Through this process, for example, due to the characteristics of BPSK modulation, it is possible to obtain distance-velocity detection results in which each target is duplicated (i.e., a ghost target is detected in addition to the actual target of each target).

그런 다음 예컨대, 일정 오경보율 (constant false alarm rate: CFAR) 알고리즘을 적용하여 각 표적에 해당하는 거리-속도 피크값과 그 피크값의 위상 정보를 추출할 수 있다. 예를 들어, CFAR 알고리즘을 이용한 거리-속도 피크값들의 위상 추출 결과, 제1 수신 안테나(Rx1), 제2 수신 안테나(Rx2), 제3 수신 안테나(Rx3), 제4 수신 안테나(Rx4)를 통해 각각 거리-속도 피크값들의 위상 정보는 [Pa1, Pb1], [Pa2, Pb2], [Pa3, Pb3], [Pa4, Pb4]일 수 있다. 추출된 피크값의 위상 정보는 저장수단에 기록될 수 있다(S220).Then, for example, the constant false alarm rate (CFAR) algorithm can be applied to extract the distance-velocity peak value corresponding to each target and the phase information of the peak value. For example, as a result of phase extraction of distance-velocity peak values using the CFAR algorithm, the first receiving antenna (Rx1), the second receiving antenna (Rx2), the third receiving antenna (Rx3), and the fourth receiving antenna (Rx4) are The phase information of each distance-velocity peak value may be [P a1 , P b1 ], [P a2 , P b2 ], [P a3 , P b3 ], and [P a4 , P b4 ]. The phase information of the extracted peak value may be recorded in a storage means (S220).

그런 다음, 추출된 거리-속도 피크들의 위상 값들을 등차수열을 이루도록 배열할 수 있다(S230). 식 (10)의 8개의 위상 값 즉, 추출된 거리-속도 피크들의 위상 값들(Pa1 - Pa4, Pb1 - Pb4)이 등차수열을 이루도록 배열한다 (S230). 위 예에서, 제1 내지 제4 수신 안테나(Rx1 - Rx4)를 통해 수신된 신호들에서 추출된 위상 값들의 조합은 [Pa1 Pa2 Pa3 Pa4 Pb1 Pb2 Pb3 Pb4]과 [Pb1 Pb2 Pb3 Pb4 Pa1 Pa2 Pa3 Pa4]이 존재할 수 있다. 이 중에서 한 가지 조합만이 등차수열을 이룰 수 있다.Then, the phase values of the extracted distance-velocity peaks can be arranged to form an arithmetic sequence (S230). The eight phase values of Equation (10), that is, the phase values (P a1 - P a4 , P b1 - P b4 ) of the extracted distance-velocity peaks, are arranged to form an arithmetic sequence (S230). In the above example, the combination of phase values extracted from the signals received through the first to fourth receiving antennas (Rx1 - Rx4) is [P a1 P a2 P a3 P a4 P b1 P b2 P b3 P b4 ] and [ P b1 P b2 P b3 P b4 P a1 P a2 P a3 P a4 ] may exist. Among these, only one combination can form an arithmetic sequence.

등차수열을 이루도록 배열된 위상값 시퀀스가 얻어지면, 그 등차수열로 배열된 위상값 시퀀스에 기초하여 수신 안테나 어레이(50)를 통해 수신된 수신신호들을 각 송신 안테나 별로 구분할 수 있다(S300). 즉, 어떤 거리-속도 피크가 어떤 송신 안테나에 해당하는지 결정할 수 있다. 특히, 첫 번째 서브 위상 시퀀스는 첫 번째 송신 안테나로부터의 신호에 해당하고, 두 번째 서브 위상 시퀀스는 두 번째 송신 안테나로부터의 신호에 해당하며, 세 번째 이하의 서브 위상 시퀀스들도 이런 식으로 세 번째 이하의 송신안테나로부터의 신호에 속해야 한다. 여기서, 각 서브 위상 시퀀스는 수신 안테나의 개수와 동일한 위상들의 시퀀스이다. 즉, 송신 안테나가 N(N은 2이상의 자연수)개 있는 경우, k 번째(여기서, k는 1부터 N까지의 자연수) 위상 시퀀스는 제k 송신 안테나(Txk)로부터의 신호에 해당한다. 이런 관계를 이용하면, 실제 표적에 해당하는 피크와 고스트 표적에 해당하는 피크를 결정할 수 있고, 결국 수신신호들을 각 송신 안테나 별로 구분할 수 있다. 예컨대, 위상 시퀀스 [Pa1 Pa2 Pa3 Pa4 Pb1 Pb2 Pb3 Pb4]가 등차수열을 이루는 경우, 선행하는 4개의 위상 [Pa1 Pa2 Pa3 Pa4] 즉, 제1 서브 위상 시퀀스는 제1 송신 안테나(Tx1)로부터 전송되어 제1 내지 제4 수신 안테나(Rx1 - Rx4)를 통해 수신된 4개의 수신신호의 위상을 나타내고, 후행하는 4개의 위상 [Pb1 Pb2 Pb3 Pb4] 즉, 제2 서브 위상 시퀀스는 제2 송신 안테나(Tx2)에서 전송되어 제1 내지 제4 수신 안테나(Rx1 - Rx4)를 통해 수신된 4개의 수신신호의 위상을 나타낸다. 여기서, 수신 신호의 파장과 입사각이 각각 λ과 θ인 경우 인접 안테나 간의 위상차(ω)는 이다.Once the phase value sequence arranged to form an arithmetic sequence is obtained, the received signals received through the receive antenna array 50 can be distinguished for each transmit antenna based on the phase value sequence arranged as an arithmetic sequence (S300). In other words, it is possible to determine which distance-velocity peak corresponds to which transmit antenna. In particular, the first sub-phase sequence corresponds to the signal from the first transmit antenna, the second sub-phase sequence corresponds to the signal from the second transmit antenna, and so on for the third and subsequent sub-phase sequences. It must belong to signals from the following transmitting antennas. Here, each sub-phase sequence is a sequence of phases equal to the number of receiving antennas. That is, when there are N transmitting antennas (N is a natural number of 2 or more), the kth (where k is a natural number from 1 to N) phase sequence corresponds to a signal from the kth transmitting antenna (Txk). Using this relationship, the peak corresponding to the actual target and the peak corresponding to the ghost target can be determined, and ultimately the received signals can be distinguished for each transmitting antenna. For example, if the phase sequence [P a1 P a2 P a3 P a4 P b1 P b2 P b3 P b4 ] forms an arithmetic sequence, the four preceding phases [P a1 P a2 P a3 P a4 ], that is, the first sub-phase The sequence represents the phases of four received signals transmitted from the first transmitting antenna (Tx1) and received through the first to fourth receiving antennas (Rx1 - Rx4), and the following four phases [P b1 P b2 P b3 P b4 ] That is, the second sub-phase sequence represents the phases of four received signals transmitted from the second transmitting antenna (Tx2) and received through the first to fourth receiving antennas (Rx1 - Rx4). Here, when the wavelength and incident angle of the received signal are λ and θ, respectively, the phase difference (ω) between adjacent antennas is am.

위와 같은 원리에 기초하여, 수신 안테나 어레이(50)를 통해 수신된 신호들을 각 송신 안테나별로 신호별로 구분한 다음, 각 송신 안테나별로 구분된 수신 신호에 대하여 속도 및 각도 추정 알고리즘을 적용하여 각 표적의 속도와 각도에 관한 추정정보를 산출할 수 있다(S400). 각 수신 안테나에 수신된 신호가 어느 송신 안테나로부터 송신된 신호인지를 알 수 있으므로, 속도 및 각도 추정에 있어서 이른바 도플러각 모호성이 해소되어 정확한 추정이 이루어질 수 있다. Based on the above principle, the signals received through the receiving antenna array 50 are classified into signals for each transmitting antenna, and then a speed and angle estimation algorithm is applied to the received signals classified for each transmitting antenna to identify each target. Estimated information about speed and angle can be calculated (S400). Since it is possible to know which transmitting antenna the signal received at each receiving antenna is from, the so-called Doppler angle ambiguity in speed and angle estimation is resolved and accurate estimation can be made.

도 11은 본 발명의 예시적인 실시예에 따라 2개의 송신 안테나와 4개의 수신 안테나를 사용하여 시뮬레이션 한 거리-속도 추정 결과를 보여준다: (A)는 3차원도이고, (B)는 위에서 내려본 2차원도이다. Figure 11 shows simulated distance-velocity estimation results using two transmit antennas and four receive antennas according to an exemplary embodiment of the present invention: (A) is a three-dimensional view, and (B) is a top-down view. It is a two-dimensional diagram.

도 11을 참조하면, 시뮬레이션에는 2개의 송신 안테나(Tx1, Tx2)와 4개의 수신 안테나(Rx1 - Rx4)를 사용하였고, 시뮬레이션 파라미터는 표 1에 따라 설정되었다. 동일한 Range bin에 여러 표적이 있을 때 제안하는 방법을 적용할 수 있음을 입증하기 위해 2개의 표적의 거리, 속도, 각도를 (30m, -15m/s, 20°)와 (30m, 8m/s, -12°)로 각각 설정하였다. CFAR 알고리즘을 사용하여 얻어진 거리-속도 피크의 위상이 도면에 표시되어 있다. 수신 안테나(Rx1 - Rx4)의 수신신호에 있어서, 도시된 것과 같이 각 표적에 대해 2개의 피크가 발생하여 총 4개의 거리-속도 피크가 발생했다. 도플러 모호성으로 인해 최대 감지 가능한 속도의 절반만큼 이동된 또 다른 피크가 발생하기 때문에 (거리, 속도) = (30m, -15m/s) 및 (30m, 4.48m/s)가 제1 표적에 의한 것이고, (거리, 속도) = (30m, -11.48m/s) 및 (30m, 8m/s)는 제2 표적으로 인한 것이다. 따라서 두 표적의 도플러 모호성을 독립적으로 해결할 수 있다.Referring to FIG. 11, two transmitting antennas (Tx1, Tx2) and four receiving antennas (Rx1 - Rx4) were used in the simulation, and the simulation parameters were set according to Table 1. To prove that the proposed method can be applied when there are multiple targets in the same range bin, the distance, speed, and angle of the two targets are (30m, -15m/s, 20°) and (30m, 8m/s, -12°), respectively. The phase of the distance-velocity peak obtained using the CFAR algorithm is shown in the figure. In the received signals of the receiving antennas (Rx1 - Rx4), two peaks occurred for each target as shown, resulting in a total of four distance-velocity peaks. Since Doppler ambiguity results in another peak shifted by half the maximum detectable velocity, (distance, velocity) = (30m, -15m/s) and (30m, 4.48m/s) are due to the first target. , (distance, speed) = (30m, -11.48m/s) and (30m, 8m/s) are due to the secondary target. Therefore, the Doppler ambiguity of the two targets can be resolved independently.

먼저, 제1 표적의 실제 속도가 -15m/s 또는 4.48m/s인지 조사한다. 피크 위치 (30m, -15m/s)에 해당하는 위상 시퀀스는 Pa = [40.23° 101.43° 166.01° -134.15°]이고, 피크 위치(30m, 4.48m/s)에 해당하는 위상 시퀀스는 Pb = [-72.99° -10.48° 51.01° 110.53°]이었다. 이 두 위상 시퀀스를 배열함으로써 우리는 캐스케이드 시퀀스 (cascaded sequence) [Pa Pb]=[40.23° 101.43° 166.01° -134.15° -72.99° -10.48° 51.01° 110.53°]가 실질적인 등차수열을 형성한다는 것을 알 수 있다. 따라서 Pa는 제1 송신 안테나(Tx1)의 신호에 해당하고 Pb는 제2 송신 안테나(Tx2)의 신호에 해당한다고 결정할 수 있다. 또한 제1 송신 안테나(Tx1)의 신호는 일반 FMCW 신호이고, 제2 송신 안테나(Tx2)의 신호는 부호 반전 FMCW 신호이다. 따라서 -15m/s는 표적의 실제 속도를 나타내고, 4.48m/s는 제2 송신 안테나(Tx2)의 부호 반전으로 인한 고스트 표적에 해당한다는 결론을 내릴 수 있다. 이와 같이 각 송신 안테나의 신호를 구별함으로써 도플러 모호성이 해결되었고 최대 감지 속도를 절반으로 줄이지 않고 표적의 속도가 고유하게 결정되었다.First, check whether the actual speed of the first target is -15m/s or 4.48m/s. The phase sequence corresponding to the peak position (30m, -15m/s) is P a = [40.23° 101.43° 166.01° -134.15°], and the phase sequence corresponding to the peak position (30m, 4.48m/s) is P b = [-72.99° -10.48° 51.01° 110.53°]. By arranging these two topological sequences, we can see that the cascaded sequence [P a P b ]=[40.23° 101.43° 166.01° -134.15° -72.99° -10.48° 51.01° 110.53°] forms a substantial arithmetic sequence. You can see that Therefore, it can be determined that P a corresponds to the signal of the first transmission antenna (Tx1) and P b corresponds to the signal of the second transmission antenna (Tx2). Additionally, the signal of the first transmission antenna (Tx1) is a normal FMCW signal, and the signal of the second transmission antenna (Tx2) is a sign-inverted FMCW signal. Therefore, it can be concluded that -15m/s represents the actual speed of the target, and 4.48m/s corresponds to a ghost target due to sign inversion of the second transmitting antenna (Tx2). By distinguishing the signal from each transmit antenna in this way, Doppler ambiguity was resolved and the target's velocity was uniquely determined without halving the maximum detection rate.

다음으로 제2 표적의 실제 속도가 -11.48m/s인지 또는 8m/s인지 확인하기 위해 동일한 과정을 반복한다. 피크 위치(30m, -11.48m/s)에 해당하는 위상 시퀀스는 Pc = [3.26° -34.72° -75.04° -111.81°]이고, 피크 위치(30m, 4.48 m/s)에 해당하는 위상 시퀀스는 Pd = [152.18° 112.91° 75.75° 40.42°]였다. 캐스캐이드 시퀀스 [Pd Pc]가 등차수열을 형성하기 때문에, 위상 시퀀스 Pd 가 제1 송신 안테나 (Tx1)의 신호에 해당하고, 표적의 실제 속도가 8m/s라고 결정할 수 있다. Next, the same process is repeated to determine whether the actual speed of the second target is -11.48 m/s or 8 m/s. The phase sequence corresponding to the peak position (30 m, -11.48 m/s) is P c = [3.26° -34.72° -75.04° -111.81°], and the phase sequence corresponding to the peak position (30 m, 4.48 m/s) was P d = [152.18° 112.91° 75.75° 40.42°]. Since the cascade sequence [P d P c ] forms an arithmetic sequence, it can be determined that the phase sequence P d corresponds to the signal of the first transmitting antenna (Tx1), and the actual speed of the target is 8 m/s.

도 12는 위상 시퀀스 [Pa Pb] 및 [Pd Pc]에 해당하는 8개의 가상 안테나의 래핑되지 않은(unwrapped) 위상 값을 보여준다. 연속된 위상 값의 차이를 180° 미만으로 만들어 불연속성을 방지하기 위해 위상 값을 래핑하지 않았다(unwrap). 올바르게 캐스캐이드된 위상 시퀀스가 등차수열을 형성한다는 것은 도면으로부터 분명하다.Figure 12 shows the unwrapped phase values of eight virtual antennas corresponding to the phase sequences [P a P b ] and [P d P c ]. To prevent discontinuity by making the difference between consecutive phase values less than 180°, the phase values were not wrapped (unwrap). It is clear from the figure that a correctly cascaded phase sequence forms an arithmetic sequence.

배열된 위상 시퀀스를 사용하여 각도 추정을 동시에 수행할 수도 있다. 각 표적의 각도는 각 표적별 위상 시퀀스에 고속 푸리에변환을 적용하여 추정될 수 있다. 각 송신 안테나로부터의 신호들이 분리되어 있기 때문에 움직이는 표적이나 위상 보상과 같은 문제가 없다. 예를 들어, (30m, -15m/s)에 해당하는 제1 표적의 각도는 위상 시퀀스 [Pa Pb]에 FFT를 적용하여 추정할 수 있다. 마찬가지로 FFT를 위상 시퀀스 [Pc Pd]에 적용하여 (30m, 8m/s)에 해당하는 제2 표적의 각도를 추정할 수 있다. Angle estimation can also be performed simultaneously using an arrayed phase sequence. The angle of each target can be estimated by applying fast Fourier transform to the phase sequence for each target. Because the signals from each transmit antenna are separated, there are no problems such as moving targets or phase compensation. For example, the angle of the first target corresponding to (30m, -15m/s) can be estimated by applying FFT to the phase sequence [P a P b ]. Similarly, FFT can be applied to the phase sequence [P c P d ] to estimate the angle of the second target corresponding to (30 m, 8 m/s).

도 13은 등차수열로 배열된 위상 시퀀스에 FFT를 적용하여 제1 및 제2 표적의 각도를 추정한 결과를 보여준다. Figure 13 shows the results of estimating the angles of the first and second targets by applying FFT to the phase sequence arranged in an arithmetic sequence.

도 13을 참조하면, 제1 및 제2 표적의 각도는 20°와 -12°의 각도로 추정되었다. 이 추정된 각도는 제1 및 제2 표적의 실제 각도와 정확하게 일치함을 알 수 있다. Referring to FIG. 13, the angles of the first and second targets were estimated to be 20° and -12°. It can be seen that this estimated angle exactly matches the actual angles of the first and second targets.

또한, SNR을 달리하여 제안한 방법의 성능을 분석하였다. 도 14는 다양한 SNR에 대한 8개의 가상 안테나의 래핑되지 않은 위상 값을 보여준다. SNR이 0dB일 때, 비래핑(unwrapped) 위상은 직선을 형성하여 위상 시퀀스가 등차수열을 형성했음을 나타낸다. SNR이 실제 시나리오에서의 SNR보다 훨씬 낮은 -20dB인 경우에도 비래핑 위상 값은 유사한 결과를 보여주었다. 반면에, SNR이 -40dB이고 표적의 각도가 정확하게 추정되지 않는 경우 큰 오차가 발생하였다.Additionally, the performance of the proposed method was analyzed by varying the SNR. Figure 14 shows the unwrapped phase values of eight virtual antennas for various SNRs. When the SNR is 0dB, the unwrapped phase forms a straight line, indicating that the phase sequence forms an arithmetic progression. Even when the SNR was -20 dB, which was much lower than the SNR in real scenarios, the non-wrapping phase values showed similar results. On the other hand, when the SNR was -40dB and the target angle was not accurately estimated, a large error occurred.

본 발명은 송신 안테나의 개수가 3개 이상인 경우에도 확대 적용될 수 있다. 이를 위해 송신 안테나의 개수에 따라 전송되는 신호의 위상을 변조하는 일반적인 PSK(Phase-Shift Keying) 변조 방식을 사용할 수 있다. 이하에서는 편의상 3개의 송신 안테나 시스템을 예로 하여 설명하지만, 결과는 임의의 갯수의 송신 안테나로 확장될 수 있다.The present invention can be expanded and applied even when the number of transmission antennas is three or more. For this purpose, a general PSK (Phase-Shift Keying) modulation method can be used, which modulates the phase of the transmitted signal according to the number of transmission antennas. Below, for convenience, a three transmit antenna system is used as an example, but the results can be extended to any number of transmit antennas.

도 15는 3개의 송신 안테나를 사용하는 3-PSK MIMO 레이더 시스템의 전송 신호를 보여준다.Figure 15 shows the transmission signal of a 3-PSK MIMO radar system using three transmit antennas.

도 15를 참조하면, 처프에 1, 2, 3의 코드가 할당되면 전송 신호의 위상은 0, 2π/3 및 4π/3으로 각각 변조된다. 각기. 다시 말해, 아래의 신호들이 사용된다.Referring to FIG. 15, when codes of 1, 2, and 3 are assigned to the chirp, the phases of the transmission signal are modulated to 0, 2π/3, and 4π/3, respectively. each. In other words, the signals below are used.

......(12) ......(12)

......(13) ......(13)

......(14) ......(14)

도 15에 도시된 바와 같이, 제1 송신 안테나(Tx1)는 111111의 코드로 송신하고, 제2 송신 안테나(Tx2)는 123123의 코드로 송신하고, 제3 송신 안테나(Tx3)는 132132의 코드로 송신한다. 이 경우, 수신된 신호는 아래 식 (15)와 같이 표현될 수 있다. As shown in FIG. 15, the first transmission antenna (Tx1) transmits with a code of 111111, the second transmission antenna (Tx2) transmits with a code of 123123, and the third transmission antenna (Tx3) transmits with a code of 132132. Send. In this case, the received signal can be expressed as equation (15) below.

......(15)......(15)

첫 번째 항은 2π(2v/λ)Tc만큼 위상이 증가하는 제1 송신 안테나(Tx1)에서 전송되는 신호에 해당한다. 또한, 두 번째 항과 세 번째 항은 제2 송신 안테나(Tx2)와 제3 송신 안테나(Tx3)에서 각각 전송되는 신호에 해당하며, 위상이 2π(2v/λ) Tc+2π/3와 2π(2v/λ)Tc+4π/3만큼 각각 증가한다. 그 결과, IF 대역 신호에 2D FFT를 적용하면 거리-속도 검출 결과에서 3개의 거리-속도 피크가 발생한다. 하나의 피크는 표적의 실제 속도에 해당하고, 나머지 두 피크는 속도가 최대 탐지가능 속도의 1/3과 2/3만큼 이동한 고스트 표적에 해당한다.The first term corresponds to a signal transmitted from the first transmitting antenna (Tx1) whose phase increases by 2π(2v/λ)T c . In addition, the second and third terms correspond to signals transmitted from the second transmission antenna (Tx2) and the third transmission antenna (Tx3), respectively, and the phases are 2π(2v/λ) T c +2π/3 and 2π (2v/λ)T c +4π/3 each increases. As a result, when 2D FFT is applied to the IF band signal, three distance-velocity peaks occur in the distance-velocity detection result. One peak corresponds to the actual speed of the target, and the remaining two peaks correspond to ghost targets whose speeds have moved by 1/3 and 2/3 of the maximum detectable speed.

도 16은 본 발명의 예시적인 실시예에 따라 3개의 송신 안테나와 4개의 수신 안테나를 갖는 3-PSK MIMO 레이더 시스템을 사용하여 시뮬레이션 한 거리-속도 추정 결과를 보여준다: (A)는 3차원도이고, (B)는 위에서 내려본 2차원도이다. Figure 16 shows range-velocity estimation results simulated using a 3-PSK MIMO radar system with three transmit antennas and four receive antennas according to an exemplary embodiment of the present invention: (A) is a three-dimensional diagram; , (B) is a two-dimensional diagram viewed from above.

도 16을 참조하면, 도플러-각도 모호성을 해결하는 방법을 보여준다. 3개의 송신 안테나와 4개의 수신 안테나가 있다고 가정하고, 전송된 신호는 도 15와 같이 변조되었다. 또한 표적의 거리, 속도, 각도는 (50m, 14m/s, 15°)로 설정되었다. 도 16에 도시된 것처럼, 3-PSK 변조의 특성으로 인해 거리-속도 평면에서 3개의 피크가 발생하였다. 3개의 피크 각각이 어느 송신 안테나에 해당하는지 결정하기 위해, 획득한 위상 값에 대해 등차수열 체크를 수행했다. 피크 위치(50m, -11.97m/s)에 해당하는 위상 시퀀스는 Pa = [-44.54° 2.21° 47.86° 94.98°]였으며, 피크 위치(50m, 1.01m/s)에 해당하는 위상 시퀀스는 Pb = [140.71° -172.53° -125.80° -79.10°]이었고, 피크 위치(50m, 14m/s)에 해당하는 위상 시퀀스는 Pc = [128.34° 174.43° -138.15° -91.51°]이었다. 이 세 가지 위상 시퀀스 Pa, Pb, Pc 를 배열함으로써 캐스캐이드 위상 시퀀스 [Pa Pb Pc]가 등차수열을 형성한다는 것을 알 수 있다. 따라서 Pa, Pb, Pc는 제1 송신 안테나(Tx1), 제2 송신 안테나(Tx2) 및 제3 송신 안테나(Tx3)로부터 송신된 신호에 대응한다고 결론을 내릴 수 있다. 제1 송신 안테나 (Tx1)의 신호는 위상이 2π(2v/λ)Tc만큼 증가하는 일반적인 FMCW 신호이므로 표적의 실제 속도는 14m/s로 추정할 수 있다. Referring to Figure 16, a method for resolving Doppler-angle ambiguity is shown. Assuming there are three transmitting antennas and four receiving antennas, the transmitted signal was modulated as shown in FIG. 15. Additionally, the target's distance, speed, and angle were set to (50m, 14m/s, 15°). As shown in Figure 16, three peaks occurred in the distance-velocity plane due to the characteristics of 3-PSK modulation. To determine which transmitting antenna each of the three peaks corresponded to, an arithmetic sequence check was performed on the obtained phase values. The phase sequence corresponding to the peak location (50 m, -11.97 m/s) was P a = [-44.54° 2.21° 47.86° 94.98°], and the phase sequence corresponding to the peak location (50 m, 1.01 m/s) was P a b = [140.71° -172.53° -125.80° -79.10°], and the phase sequence corresponding to the peak position (50 m, 14 m/s) was P c = [128.34° 174.43° -138.15° -91.51°]. By arranging these three phase sequences P a , P b , and P c , we can see that the cascade phase sequence [P a P b P c ] forms an arithmetic sequence. Therefore, it can be concluded that P a , P b , and P c correspond to signals transmitted from the first transmission antenna (Tx1), the second transmission antenna (Tx2), and the third transmission antenna (Tx3). Since the signal of the first transmitting antenna (Tx1) is a general FMCW signal whose phase increases by 2π(2v/λ)T c , the actual speed of the target can be estimated to be 14m/s.

또한, 위상 시퀀스 [Pc Pa Pb]에 FFT를 적용하여 도 13과 유사한 방식으로 표적의 각도 추정을 수행할 수 있다. 따라서 본 발명의 예시적인 실시예에 따른 방법은 3-PSK 변조 방식을 사용하는 3 송신 안테나 시스템에도 효과적으로 적용되어 도플러-각도 모호성을 해결할 수 있다. 본 발명의 방법은 송신 안테나의 개수가 4개 이상인 경우에도 확대 적용될 수 있다. 즉, 일반화하면, 본 발명의 방법은 N-PSK MIMO 레이더 시스템(단, N은 2 이상의 자연수)에 적용되어 도플러-각도 모호성을 해결하면서 표적들에 대한 거리, 속도, 각도를 정확하고 효과적으로 추정할 수 있다.Additionally, angle estimation of the target can be performed in a manner similar to FIG. 13 by applying FFT to the phase sequence [P c P a P b ]. Therefore, the method according to an exemplary embodiment of the present invention can be effectively applied to a three-transmission antenna system using the 3-PSK modulation method to resolve Doppler-angle ambiguity. The method of the present invention can be expanded and applied even when the number of transmission antennas is four or more. In other words, in general terms, the method of the present invention can be applied to an N-PSK MIMO radar system (where N is a natural number of 2 or more) to accurately and effectively estimate the distance, speed, and angle to targets while resolving Doppler-angle ambiguity. You can.

위와 같은 시뮬레이션 결과는 FMCW 레이더 센서를 이용한 실제 실험을 통해서도 확인할 수 있었다. 실험에 사용된 FMCW 레이더 센서는 2개의 송신 안테나와 4개의 수신 안테나로 MIMO 모드를 지원하며, 전송되는 신호에 BPSK 변조를 적용할 수 있다. 정지해 있는 표적이나 움직이는 표적 모두에 대하여 측정된 거리-속도 평면에서의 피크들의 캐스케이드 위상 시퀀스가 등차수열을 형성하였고, 최대 탐지가능 속도의 저하 없이 그 등차수열을 이용하여 표적의 속도와 각도를 추정한 값이 실제의 속도와 각도와 정확하게 일치함을 확인할 수 있었다. 따라서 본 발명에 따라 송신 신호의 위상을 변조하고 수신 신호의 위상을 배열하여 각 송신 안테나의 신호를 구별하는 방법은 도플러 각도 모호성을 해결하면서 표적의 명확한 속도와 각도를 추정할 수 있고, 결과적으로 MIMO 레이더 시스템의 성능을 향상시킬 수 있다. The above simulation results could also be confirmed through actual experiments using the FMCW radar sensor. The FMCW radar sensor used in the experiment supports MIMO mode with two transmitting antennas and four receiving antennas, and BPSK modulation can be applied to the transmitted signal. A cascade phase sequence of peaks in the distance-velocity plane measured for both stationary and moving targets formed an arithmetic sequence, and the arithmetic sequence was used to estimate the target's speed and angle without reducing the maximum detectable speed. It was confirmed that the value exactly matched the actual speed and angle. Therefore, the method of distinguishing the signal of each transmitting antenna by modulating the phase of the transmitted signal and arranging the phase of the received signal according to the present invention can estimate a clear speed and angle of the target while resolving Doppler angle ambiguity, and as a result, MIMO The performance of the radar system can be improved.

이상에서 설명한 바와 같이, BPSK-MIMO FMCW 레이더 시스템에서 각 표적에 해당하는 거리와 속도 값은 송신 안테나 개수만큼 복제되어 고스트 표적으로 탐지되어 도플러-각도 모호성이 생기는데, 본 발명이 제시하는 방법은 이를 해결할 수 있다. 즉, 각 수신 안테나의 거리-속도 검출 결과에서 동일한 표적에 해당하는 신호의 위상 값은 등차수열의 형태를 가져야 한다는 새로운 인식에 기반하여, 각 표적에서 파생된 위상 값을 올바르게 정렬하여 도플러 모호성을 해결할 수 있다.다중 송신 안테나를 사용하는 경우에도 본 발명에 따른 방법을 사용하여 최대 검출 가능한 속도를 유지할 수 있고, 또한 표적의 움직임으로 인한 각도 추정의 부정확성을 해결할 수 있다. 이와 같은 본 발명에 따른 방법은 시뮬레이션 데이터와 측정 데이터를 모두 사용하여 우수한 성능을 나타냄을 확인할 수 있었다. As explained above, in the BPSK-MIMO FMCW radar system, the distance and speed values corresponding to each target are duplicated as many as the number of transmitting antennas and are detected as ghost targets, resulting in Doppler-angle ambiguity. The method proposed by the present invention solves this problem. You can. In other words, based on the new recognition that the phase values of signals corresponding to the same target in the distance-velocity detection results of each receiving antenna must have the form of an arithmetic sequence, Doppler ambiguity can be resolved by correctly aligning the phase values derived from each target. Even when using multiple transmit antennas, the maximum detectable speed can be maintained using the method according to the present invention, and inaccuracy in angle estimation due to target movement can also be resolved. It was confirmed that the method according to the present invention shows excellent performance using both simulation data and measurement data.

위에서 설명된 본 발명의 실시예에 따른 방법은 처리장치(90)를 통하여 수행될 수 있는 프로그램 명령 형태로 구현되어 컴퓨터 판독 가능 매체에 기록될 수 있다. 컴퓨터 판독 가능 기록 매체의 예에는 하드 디스크, 플로피 디스크 및 자기 테이프와 같은 자기 매체(magnetic media), CD-ROM, DVD와 같은 광기록 매체(optical media), 플롭티컬 디스크(floptical disk)와 같은 자기-광 매체(magneto-optical media), 및 롬(ROM), 램(RAM), 플래시 메모리 등과 같은 프로그램 명령을 저장하고 수행하도록 특별히 구성된 하드웨어 장치가 포함된다. The method according to the embodiment of the present invention described above may be implemented in the form of program instructions that can be executed through the processing device 90 and recorded on a computer-readable medium. Examples of computer-readable recording media include magnetic media such as hard disks, floppy disks, and magnetic tapes, optical media such as CD-ROMs and DVDs, and magnetic media such as floptical disks. -Includes optical media (magneto-optical media) and hardware devices specifically configured to store and execute program instructions, such as ROM, RAM, flash memory, etc.

본 발명에 따른 방법과 시스템은 MIMO 안테나를 이용하는 레이더 시스템에 사용될 수 있다. 특히, MIMO 기반 FMCW 레이더 시스템에서 표적의 속도 및 각도 추정 성능을 높이는 데 활용될 수 있다.The method and system according to the present invention can be used in radar systems using MIMO antennas. In particular, it can be used to improve target speed and angle estimation performance in MIMO-based FMCW radar systems.

이상과 같이 실시예들이 비록 한정된 도면에 의해 설명되었으나, 해당 기술 분야의 숙련된 당업자는 하기의 특허 청구 범위에 기재된 본 발명의 사상 및 영역으로부터 벗어나지 않는 범위 내에서 본 발명을 다양하게 수정 및 변경시킬 수 있음을 이해할 수 있을 것이다. 예를 들어, 설명된 기술들이 설명된 방법과 다른 순서로 수행되거나, 및/또는 설명된 시스템, 구조, 장치, 회로 등의 구성요소들이 설명된 방법과 다른 형태로 결합 또는 조합되거나, 다른 구성요소 또는 균등물에 의하여 대치되거나 치환되더라도 적절한 결과가 달성될 수 있다. 그러므로 다른 구현들, 다른 실시예들 및 특허청구범위와 균등한 것들도 후술하는 특허청구범위의 범위에 속한다.Although the embodiments have been described with limited drawings as described above, those skilled in the art can make various modifications and changes to the present invention without departing from the spirit and scope of the present invention as set forth in the following patent claims. You will understand that it is possible. For example, the described techniques are performed in a different order than the described method, and/or components of the described system, structure, device, circuit, etc. are combined or combined in a different form than the described method, or other components are used. Alternatively, appropriate results may be achieved even if substituted or substituted by an equivalent. Therefore, other implementations, other embodiments, and equivalents to the claims also fall within the scope of the claims described below.

10: BPSK-MIMO FMCW 레이더 시스템
20: 송신부 30: 송신 안테나
40: 수신부 50: 수신 안테나
60: 디지털 신호 처루비 70: 거리/속도 추정부
80: 각도 추정부 90: 처리 장치
10: BPSK-MIMO FMCW radar system
20: Transmitting unit 30: Transmitting antenna
40: Receiving unit 50: Receiving antenna
60: Digital signal processing unit 70: Distance/speed estimation unit
80: angle estimation unit 90: processing device

Claims (19)

MIMO 레이더 시스템을 이용한 표적 탐지 방법에 있어서,
복수의 송신 안테나를 통해 송출된 레이더 송신신호가 표적에 의해 반사되어 되돌아오는 신호를 복수의 수신 안테나를 통해 동시에 수신하는 단계;
상기 복수의 수신 안테나 각각의 수신신호들에 대한 소정의 신호처리를 통해 각 수신신호의 거리-속도 검출 결과에서의 피크를 찾아내어 해당 피크의 위상정보를 추출하는 단계; 및
상기 복수의 수신 안테나 전체의 수신신호들로부터 추출된 거리-속도 피크들의 위상 값들을 등차수열을 이루도록 배열하여 그 등차수열로 배열된 위상 값 시퀀스에 기반하여 상기 수신신호들을 각 송신 안테나별로 구분하는 단계를 포함하는 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법.
In a target detection method using a MIMO radar system,
Simultaneously receiving a radar transmission signal transmitted through a plurality of transmitting antennas reflected by a target and returning a signal through a plurality of receiving antennas;
finding a peak in the distance-velocity detection result of each received signal through predetermined signal processing on the received signals from each of the plurality of receiving antennas and extracting phase information of the corresponding peak; and
Arranging the phase values of the distance-velocity peaks extracted from the received signals of all the plurality of receiving antennas to form an arithmetic sequence and classifying the received signals for each transmitting antenna based on the phase value sequence arranged in the arithmetic sequence. A target detection method based on phase information of a MIMO radar signal, comprising:
제1항에 있어서, 상기 위상정보를 추출하는 단계는 상기 수신신호의 신호처리를 통해 얻어지는 cos 채널(I 채널)의 중간 주파수(IF) 대역 신호와 sin 채널(Q 채널)의 중간 주파수(IF) 대역 신호에 대하여 2차원 고속푸리에변환(2D FFT)을 적용하여 각 표적에 관한 거리-속도 검출 결과를 얻는 단계; 및 상기 거리-속도 검출 결과에 대하여 일정 오경보율 (constant false alarm rate: CFAR) 알고리즘을 적용하여 각 표적에 해당하는 거리-속도 피크들과 그 피크들의 위상값을 추출하는 단계를 포함하는 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법. The method of claim 1, wherein the step of extracting the phase information includes the intermediate frequency (IF) band signal of the cos channel (I channel) and the intermediate frequency (IF) of the sin channel (Q channel) obtained through signal processing of the received signal. Obtaining distance-velocity detection results for each target by applying a two-dimensional fast Fourier transform (2D FFT) to the band signal; And applying a constant false alarm rate (CFAR) algorithm to the distance-speed detection result to extract distance-velocity peaks corresponding to each target and phase values of the peaks. A target detection method based on phase information of MIMO radar signals. 제1항에 있어서, 상기 레이더 송신신호는 PSK 변조방식으로 변조된 신호인 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법. The method of claim 1, wherein the radar transmission signal is a signal modulated using PSK modulation. 제1항에 있어서, 상기 복수의 송신 안테나와 상기 복수의 수신 안테나는 ULA 형태로 배열된 안테나 어레이인 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법. The method of claim 1, wherein the plurality of transmitting antennas and the plurality of receiving antennas are antenna arrays arranged in a ULA form. 제1항에 있어서, 상기 각 송신 안테나별로 구분하는 단계에서, 상기 송신 안테나가 두 개일 때, 등차수열을 이루는 거리-속도 피크의 위상 시퀀스 중 선행 서브 위상 시퀀스(former phase sequence)는 제1 송신 안테나로부터 송신된 신호들에 대응하고, 후행 서브 위상 시퀀스(latter phase sequence)는 제2 송신 안테나로부터 송신된 신호들에 대응하는 것으로 결정되며, 여기서 각 서브 위상 시퀀스는 수신 안테나의 개수와 동일한 수의 위상들의 시퀀스인 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법. The method of claim 1, wherein in the step of classifying each transmission antenna, when there are two transmission antennas, the preceding sub-phase sequence among the phase sequences of the distance-velocity peaks forming an arithmetic sequence is the first transmission antenna. corresponds to signals transmitted from the second transmit antenna, and a latter phase sequence is determined to correspond to signals transmitted from the second transmit antenna, where each sub-phase sequence has a number of phases equal to the number of receive antennas. A target detection method based on phase information of a MIMO radar signal, characterized in that it is a sequence of. 제1항에 있어서, 상기 각 송신 안테나별로 구분하는 단계에서, 상기 송신 안테나가 N(단, N은 3이상의 자연수)개 있는 경우, 등차수열을 이루는 거리-속도 피크들의 위상 시퀀스 중 k 번째(여기서, k는 1부터 N까지의 자연수) 서브 위상 시퀀스는 제k 송신 안테나로부터 송신된 신호들에 대응하는 것으로 결정되고, 여기서 각 서브 위상 시퀀스는 수신 안테나의 개수와 동일한 수의 위상들의 시퀀스인 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법. The method of claim 1, wherein in the step of classifying each transmission antenna, if there are N transmission antennas (where N is a natural number of 3 or more), the kth phase sequence of distance-velocity peaks forming an arithmetic sequence (where , k is a natural number from 1 to N) The sub-phase sequence is determined to correspond to signals transmitted from the kth transmitting antenna, where each sub-phase sequence is a sequence of phases with the same number of phases as the number of receiving antennas. A target detection method based on phase information of MIMO radar signals. 제1항에 있어서, 상기 구분하는 단계 후, 각 송신 안테나별로 구분된 수신 신호를 이용하여 각 표적의 속도와 각도에 관한 추정정보를 산출하는 단계를 더 포함하는 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법. The phase of the MIMO radar signal according to claim 1, further comprising calculating estimated information about the speed and angle of each target using the received signal divided for each transmit antenna after the step of classifying. Information-based target detection method. 제7항에 있어서, 각 표적의 각도는 각 표적별 위상 시퀀스에 고속 푸리에변환을 적용하여 추정되는 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법. The method of claim 7, wherein the angle of each target is estimated by applying fast Fourier transform to the phase sequence for each target. 제1항에 있어서, 상기 레이더 송신신호는 시간에 따라 주파수가 증가하는 처프 신호인 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법. The method of claim 1, wherein the radar transmission signal is a chirp signal whose frequency increases with time. 복수의 송신 안테나;
복수의 수신 안테나;
레이더 송신신호를 생성하여 상기 복수의 송신 안테나를 통해 무선 송신하도록 구성되는 송신부;
전방의 표적에 의해 반사되어 상기 복수의 수신 안테나를 통해 수신된 수신신호와 상기 송신신호 및 상기 송신신호의 90도 위상천이된 신호를 이용하여 cos 채널(I 채널)의 중간 주파수(IF) 대역 신호와 sin 채널(Q 채널)의 중간 주파수(IF) 대역 신호를 생성하도록 구성된 수신부;
상기 cos 채널(I 채널)의 중간 주파수 대역 신호와 sin 채널(Q 채널)의 중간 주파수 대역 신호에 대하여 필터링 및 고속 푸리에변환(FFT) 처리를 수행하도록 구성된 디지털 신호 처리부; 및
상기 디지털 신호 처리부에서 처리된 신호를 이용하여 거리-속도 피크들을 검출하고 각 피크의 위상 값을 추출하는 처리, 추출된 위상 값들을 등차수열을 이루도록 배열하고, 등차수열로 배열된 위상 값 시퀀스에 기반하여 수신신호들과 각 송신 안테나들 간의 대응관계를 결정하여 수신신호를 송신안테나별로 구분하는 처리, 그리고 각 송신안테나별로 구분된 수신신호를 이용하여 표적의 거리와 속도를 추정하는 처리를 수행하도록 구성된 거리/속도 추정부를 구비하는 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 시스템.
a plurality of transmitting antennas;
a plurality of receiving antennas;
a transmitter configured to generate a radar transmission signal and transmit it wirelessly through the plurality of transmission antennas;
An intermediate frequency (IF) band signal of the cos channel (I channel) using the received signal reflected by the target in front and received through the plurality of receiving antennas, the transmitted signal, and a 90-degree phase-shifted signal of the transmitted signal. A receiver configured to generate intermediate frequency (IF) band signals of the and sin channels (Q channels);
a digital signal processor configured to perform filtering and fast Fourier transform (FFT) processing on the intermediate frequency band signal of the cos channel (I channel) and the intermediate frequency band signal of the sin channel (Q channel); and
Process of detecting distance-velocity peaks and extracting the phase value of each peak using the signal processed by the digital signal processor, arranging the extracted phase values to form an arithmetic sequence, and based on the sequence of phase values arranged as an arithmetic sequence. It is configured to perform processing to determine the correspondence between received signals and each transmitting antenna, classify the received signal by transmitting antenna, and estimate the distance and speed of the target using the received signal classified for each transmitting antenna. A target detection system based on phase information of a MIMO radar signal, comprising a distance/speed estimator.
제10항에 있어서, 거리/속도 추정부는 상기 디지털 신호 처리부에서 고속 푸리에변환(FFT) 처리된 신호에 대하여 일정 오경보율 (constant false alarm rate: CFAR)알고리즘을 적용하여 각 표적에 해당하는 거리-속도 피크들과 그 피크들의 위상값을 추출하도록 구성된 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 시스템.The method of claim 10, wherein the distance/speed estimator applies a constant false alarm rate (CFAR) algorithm to the signal processed by the fast Fourier transform (FFT) in the digital signal processor to determine the distance-speed corresponding to each target. A target detection system based on phase information of a MIMO radar signal, characterized in that it is configured to extract peaks and phase values of the peaks. 제10항에 있어서, 상기 거리/속도 추정부는, 상기 송신 안테나가 두 개일 때, 등차수열을 이루는 거리-속도 피크의 위상 시퀀스 중 선행 서브 위상 시퀀스(former phase sequence)는 제1 송신 안테나로부터 송신된 신호들에 대응하고, 후행 서브 위상 시퀀스(latter phase sequence)는 제2 송신 안테나로부터 송신된 신호들에 대응하는 것으로 결정하도록 구성되며, 여기서 각 서브 위상 시퀀스는 수신 안테나의 개수와 동일한 수의 위상들의 시퀀스인 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 시스템.The method of claim 10, wherein, when there are two transmitting antennas, the preceding sub-phase sequence (former phase sequence) of the phase sequences of the distance-velocity peaks forming an arithmetic sequence is transmitted from the first transmitting antenna. corresponding to the signals, and is configured to determine that a later sub-phase sequence corresponds to signals transmitted from the second transmit antenna, where each sub-phase sequence has a number of phases equal to the number of receive antennas. A target detection system based on phase information of MIMO radar signals, characterized in that the sequence. 제10항에 있어서, 상기 거리/속도 추정부는, 상기 송신 안테나가 N(N은 3이상의 자연수)개 있는 경우, 등차수열을 이루는 거리-속도 피크들의 위상 시퀀스 중 k 번째(여기서, k는 1부터 N까지의 자연수) 서브 위상 시퀀스는 제k 송신 안테나로부터 송신된 신호들에 대응하는 것으로 결정하도록 구성되고, 여기서 각 서브 위상 시퀀스는 수신 안테나의 개수와 동일한 수의 위상들의 시퀀스인 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 시스템.The method of claim 10, wherein, when there are N transmitting antennas (N is a natural number of 3 or more), the kth phase sequence of distance-velocity peaks forming an arithmetic sequence (where k starts from 1) up to N natural numbers) is configured to determine the sub-phase sequence corresponding to signals transmitted from the kth transmit antenna, where each sub-phase sequence is a sequence of phases with a number equal to the number of receive antennas. Target detection system based on phase information of radar signals. 제10항에 있어서, 각 송신안테나별로 구분된 수신신호를 이용하여 소정의 각도 추정 알고리즘을 적용하여 각 표적의 각도정보를 산출하도록 구성된 각도 추정부를 더 구비하는 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 시스템.11. The phase information of the MIMO radar signal according to claim 10, further comprising an angle estimation unit configured to calculate the angle information of each target by applying a predetermined angle estimation algorithm using received signals separated for each transmission antenna. based target detection system. 제14항에 있어서, 상기 각도 추정부는 각 표적별 위상 시퀀스에 고속 푸리에변환을 적용하여 각 표적의 각도를 추정하도록 구성되는 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 시스템.The target detection system according to claim 14, wherein the angle estimation unit is configured to estimate the angle of each target by applying a fast Fourier transform to the phase sequence of each target. 제10항에 있어서, 상기 복수의 송신 안테나와 상기 복수의 수신 안테나는 ULA 형태로 배열된 안테나 어레이인 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 시스템.The target detection system based on phase information of a MIMO radar signal according to claim 10, wherein the plurality of transmitting antennas and the plurality of receiving antennas are antenna arrays arranged in a ULA form. 제10항에 있어서, 상기 레이더 송신신호는 시간에 따라 주파수가 증가하는 처프 신호인 것을 특징으로 하는 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 시스템.The target detection system based on phase information of a MIMO radar signal according to claim 10, wherein the radar transmission signal is a chirp signal whose frequency increases with time. 제1항 내지 제9항 중 어느 한 항에 기재된 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법을 수행하기 위하여 컴퓨터 판독 가능한 기록 매체에 저장된 컴퓨터 실행가능 프로그램.A computer-executable program stored in a computer-readable recording medium to perform the target detection method based on phase information of a MIMO radar signal according to any one of claims 1 to 9. 제1항 내지 제9항 중 어느 한 항에 기재된 MIMO 레이더 신호의 위상 정보 기반 표적 탐지 방법을 수행하기 위한 컴퓨터 프로그램이 기록된 컴퓨터 판독 가능한 기록 매체.A computer-readable recording medium on which a computer program for performing the target detection method based on phase information of a MIMO radar signal according to any one of claims 1 to 9 is recorded.
KR1020210171594A 2021-12-03 2021-12-03 Target detection method and system based on phase information of MIMO radar signal Active KR102641300B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020210171594A KR102641300B1 (en) 2021-12-03 2021-12-03 Target detection method and system based on phase information of MIMO radar signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020210171594A KR102641300B1 (en) 2021-12-03 2021-12-03 Target detection method and system based on phase information of MIMO radar signal

Publications (2)

Publication Number Publication Date
KR20230083525A KR20230083525A (en) 2023-06-12
KR102641300B1 true KR102641300B1 (en) 2024-02-27

Family

ID=86770146

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020210171594A Active KR102641300B1 (en) 2021-12-03 2021-12-03 Target detection method and system based on phase information of MIMO radar signal

Country Status (1)

Country Link
KR (1) KR102641300B1 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102587986B1 (en) * 2023-07-04 2023-10-11 포항공과대학교 산학협력단 Multipath signal detection apparatus based on MPSK-MIMO FMCW radar and method thereof
KR102587987B1 (en) * 2023-07-21 2023-10-11 포항공과대학교 산학협력단 Target detection apparatus using angle estimation in MPSK-MIMO FMCW radar method thereof
CN117641644B (en) * 2023-12-28 2024-09-20 重庆南天智能设施有限公司 Intelligent on-demand illumination control system based on radar sensing

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003130944A (en) 2001-10-22 2003-05-08 Mitsubishi Electric Corp On-vehicle radar device
US20210096234A1 (en) 2019-09-30 2021-04-01 Qualcomm Incorporated Radar coordination for multi-radar coexistence

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DK201570293A1 (en) 2015-05-19 2016-12-12 Nel Hydrogen As Diaphragm compressor with an oblong shaped chamber
KR102143948B1 (en) * 2018-07-23 2020-08-12 재단법인대구경북과학기술원 Apparatus and method for processing radar signal for a multi channel antenna
KR102182293B1 (en) * 2019-03-11 2020-11-24 현대모비스 주식회사 Apparatus and mehtod for estimating direction of arrival in mimo system
KR20210101957A (en) * 2020-02-11 2021-08-19 삼성전자주식회사 Method and apparatus for detecting velocity of object in radar system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003130944A (en) 2001-10-22 2003-05-08 Mitsubishi Electric Corp On-vehicle radar device
US20210096234A1 (en) 2019-09-30 2021-04-01 Qualcomm Incorporated Radar coordination for multi-radar coexistence

Also Published As

Publication number Publication date
KR20230083525A (en) 2023-06-12

Similar Documents

Publication Publication Date Title
KR102641300B1 (en) Target detection method and system based on phase information of MIMO radar signal
US12019142B2 (en) Frequency modulated continuous wave radar system
CN110824444B (en) MIMO radar decoding for resolving velocity ambiguities
US11906618B2 (en) MIMO radar apparatuses and MIMO radar methods
JP6300959B2 (en) Radar system and method for determining range, relative velocity, and orientation of an object using continuous wave and chirp signals
US9971028B2 (en) Method and apparatus for detecting target using radar
CN110535481A (en) Use the combined radar and communication system of shared signal waveform
CN106291559B (en) Method for operating a radar system having at least two transmitting antennas and at least one receiving antenna
CN108885254B (en) Object detection device
US10921434B2 (en) Radar system
US7548193B2 (en) Radar device
US7961147B1 (en) Long baseline phase interferometer ambiguity resolution using frequency differences
US20180143297A1 (en) Decentralised radar system
US20130229300A1 (en) On-board radar apparatus, object detection method, and object detection program
US9568601B1 (en) Successive-MFCW modulation for ultra-fast narrowband radar
EP3270180B1 (en) Signal processing apparatus for generating a range-doppler map
CN114002654A (en) A method for resolving velocity ambiguity in millimeter-wave radar
KR101125276B1 (en) Method and Apparatus for Detecting Range and Velocity of Target by Using Radar and Computer Readable Recording Medium for Recording Program Therefor
KR102169873B1 (en) Mehtod for estimating number of signals for vehicle radar using relative pomparison of eigenvalues
Schroeder et al. X-band FMCW radar system with variable chirp duration
JP7345099B2 (en) Radar device and radar method
US20240019565A1 (en) Motion compensation for fast target detection in automotive radar
KR20200053222A (en) Radar apparatus and radar signal processing method for precise measurement of distance, angular velocity
KR20180101002A (en) Angle Estimation Method and Apparatus for Automotive Radars
EP4390454A1 (en) Slow phase-time division multiple access (sp-tdma) modulation scheme in a mimo radar

Legal Events

Date Code Title Description
PA0109 Patent application

Patent event code: PA01091R01D

Comment text: Patent Application

Patent event date: 20211203

PA0201 Request for examination
PG1501 Laying open of application
E701 Decision to grant or registration of patent right
PE0701 Decision of registration

Patent event code: PE07011S01D

Comment text: Decision to Grant Registration

Patent event date: 20240130

GRNT Written decision to grant
PR0701 Registration of establishment

Comment text: Registration of Establishment

Patent event date: 20240222

Patent event code: PR07011E01D

PR1002 Payment of registration fee

Payment date: 20240223

End annual number: 3

Start annual number: 1

PG1601 Publication of registration