[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

KR102001394B1 - Method of estimating DOA of received signals based on logarithmic-domain antenna array interpolation, and apparatus for the same - Google Patents

Method of estimating DOA of received signals based on logarithmic-domain antenna array interpolation, and apparatus for the same Download PDF

Info

Publication number
KR102001394B1
KR102001394B1 KR1020170088392A KR20170088392A KR102001394B1 KR 102001394 B1 KR102001394 B1 KR 102001394B1 KR 1020170088392 A KR1020170088392 A KR 1020170088392A KR 20170088392 A KR20170088392 A KR 20170088392A KR 102001394 B1 KR102001394 B1 KR 102001394B1
Authority
KR
South Korea
Prior art keywords
antenna array
array elements
interpolated
received signals
interpolation
Prior art date
Application number
KR1020170088392A
Other languages
Korean (ko)
Other versions
KR20190007221A (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 KR1020170088392A priority Critical patent/KR102001394B1/en
Publication of KR20190007221A publication Critical patent/KR20190007221A/en
Application granted granted Critical
Publication of KR102001394B1 publication Critical patent/KR102001394B1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14Systems for determining direction or deviation from predetermined direction
    • G01S3/46Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
    • G01S3/48Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems the waves arriving at the antennas being continuous or intermittent and the phase difference of signals derived therefrom being measured
    • 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
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/04Details
    • G01S3/12Means for determining sense of direction, e.g. by combining signals from directional antenna or goniometer search coil with those from non-directional antenna
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0257Control of position or course in two dimensions specially adapted to land vehicles using a radar

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법과 이를 위한 장치가 개시된다. 원래의 안테나 어레이 요소들의 제1 스티어링 행렬과, 원래의 안테나 어레이 요소들에 기초하여 보간한 보간 안테나 어레이 요소들의 제2 스티어링 행렬 각각에 대해 로그를 취한다. 로그 영역에서 제1 스티어링 행렬을 제2 스티어링 행렬로 변환하는 변환행렬을 선형 최소 제곱법과 같은 최적화 기법을 이용하여 구한다. 원래의 안테나 어레이 요소들의 수신신호들의 위상에 대하여 상기 변환행렬을 적용하여 보간 안테나 어레이 요소들에 대응되는 위상을 갖는 보간 수신신호들로 변환한다. 소정의 도래각(DOA) 추정 알고리즘에 기초하여 보간 안테나 어레이 요소들의 보간 수신신호들의 DOA를 추정한다. 보간 안테나 어레이 요소들의 보간 수신신호들 간의 전력 레벨의 편차를 보상하여 보간 수신신호들의 전력 레벨을 실질적으로 균등하게 보정함으로써 DOA 추정의 정확도를 더 높일 수 있다. 이러한 DOA 추정 방법은 프로그램으로 구현되어 디지털 신호처리기와 같은 컴퓨팅 장치를 사용하여 실시할 수 있다.A method for estimating the arrival angle of a received signal based on log-area antenna array interpolation and an apparatus therefor are disclosed. A first steering matrix of the original antenna array elements and a second steering matrix of interpolated antenna array elements based on the original antenna array elements. A transformation matrix for transforming the first steering matrix into the second steering matrix in the logarithmic region is obtained by using an optimization technique such as a linear least squares method. And applies the transformation matrix to the phases of the received signals of the original antenna array elements to convert them into interpolated reception signals having phases corresponding to the interpolation antenna array elements. Estimates the DOA of the interpolated received signals of the interpolation antenna array elements based on a predetermined DOA estimation algorithm. The accuracy of the DOA estimation can be further improved by compensating for the power level deviation between the interpolated received signals of the interpolation antenna array elements to substantially evenly compensate the power level of the interpolated received signals. Such a DOA estimation method can be implemented as a program and can be implemented using a computing device such as a digital signal processor.

Description

로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법과 이를 위한 장치 {Method of estimating DOA of received signals based on logarithmic-domain antenna array interpolation, and apparatus for the same}Field of the Invention The present invention relates to a method of estimating an arrival angle of a received signal based on log-area antenna array interpolation,

본 발명은 안테나 기술 분야에 관한 것으로서, 보다 상세하게는 균일 선형 어레이(uniform linear array: ULA) 안테나 어레이를 기반으로 하여 수신신호들의 도래각(direction of arrival: DOA)을 효과적이고 정확하게 추정하는 방법과 이를 위한 장치에 관한 것이다.BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an antenna technology, and more particularly, to a method of effectively and accurately estimating a direction of arrival (DOA) of received signals based on a uniform linear array (ULA) And a device for this.

자율 주행 차량 시스템에서 가장 중요한 문제 중 하나는 사람의 안전이다. 예기치 않은 자동차 사고를 예방하기 위해 카메라, 무선 탐지 및 거리 측정 (radio detection and ranging: RADAR) 및 광 탐지 및 거리 측정 (light detection and ranging: LiDAR)과 같은 많은 센서가 자율 차량에 사용되고 있다. 다양한 표적 감지 센서들 중에서 악천후에 사용될 수 있는 레이더가 널리 연구되고 있다.One of the most important problems in autonomous vehicle systems is human safety. Many sensors such as cameras, radio detection and ranging (RADAR) and light detection and ranging (LiDAR) are used in autonomous vehicles to prevent unexpected car accidents. Among various target detection sensors, radar that can be used in bad weather is widely studied.

예를 들어, 전방에 2대의 차량이 레이더 센서로부터 같은 거리에 서로 인접해 있을 때, 레이더는 전방의 차량이 1대가 아닌 2대가 존재한다고 인식할 수 있어야 한다. 표적들의 속도 및 거리(range)가 비슷한 경우 여러 표적을 분리하기 위해 표적의 DOA을 정확하게 추정하는 것이 중요하다. 따라서 다중 신호 분류 (multiple signal classification: MUSIC) 알고리즘과 회전 불변 기법을 통한 신호 파라미터의 추정(estimation of signal parameters via rotational invariance technique: ESPRIT) 알고리즘과 같은 고해상도 DOA 추정 알고리즘이 지난 수십 년 동안 연구되어 왔다. 또한 최근 몇 년 동안에는 MUSIC 알고리즘보다 신호 대 잡음비(SNR)의 영향을 덜 받는 바틀렛(Bartlett) 알고리즘이 관심을 받고 많이 사용되고 있다.For example, when two cars in front are adjacent to each other at the same distance from the radar sensor, the radar must be able to recognize that there are two cars ahead of them. If the speeds and ranges of the targets are similar, it is important to accurately estimate the DOA of the target in order to separate multiple targets. Thus, high-resolution DOA estimation algorithms such as the multiple signal classification (MUSIC) algorithm and the estimation of signal parameters via rotational invariance technique (ESPRIT) algorithms have been studied for decades. Also, in recent years, the Bartlett algorithm, which is less affected by the signal-to-noise ratio (SNR) than the MUSIC algorithm, has received much attention and is being widely used.

도착 방향 (DOA) 추정 알고리즘의 정확성을 향상시키기 위해, 레이더의 안테나 어레이 보간 기술도 지속적으로 연구되어 오고 있다. 안테나 어레이 요소들을 원래 위치에서 원하는 위치로 이동하려면 변환행렬(transformation matrix)을 적용해야 한다. 일반적으로, 적절한 변환행렬을 찾기 위해 선형 최소 제곱법 (linear least squares: LLS)이 널리 사용된다. 그러나 LLS 방법으로 얻은 변환행렬은 안테나 어레이 요소들을 보간하기 위한 최적의 솔루션이 아니다. 이 변환행렬을 적용할 때, 원래의 안테나 어레이 요소들의 선형 조합에 의해 보간된 안테나 어레이 요소들이 생성된다. 따라서 보간된 안테나 어레이 요소들의 진폭은 원래의 안테나 어레이 요소들의 진폭과 다를 수 있다. 안테나 어레이 요소들 사이에 진폭 차이가 존재하면 DOA 추정 알고리즘의 성능이 저하된다. 또한, 진폭 및 위상의 차이를 동시에 최소화 하는 과정에서 LLS법의 해답이 얻어지므로, 보간된 안테나 어레이 요소들의 위상 정보가 정확하게 체계적으로 나타내어지지 못한다. 이런 점들 때문에 DOA 추정 성능이 좋지 않게 된다.In order to improve the accuracy of the DOA estimation algorithm, the antenna array interpolation technique of radar has been continuously studied. To move the antenna array elements from their original position to the desired location, a transformation matrix must be applied. In general, linear least squares (LLS) are widely used to find appropriate transform matrices. However, the transformation matrix obtained by the LLS method is not an optimal solution for interpolating antenna array elements. When applying this transformation matrix, the interpolated antenna array elements are generated by a linear combination of the original antenna array elements. Thus, the amplitude of the interpolated antenna array elements may be different from the amplitude of the original antenna array elements. If there is an amplitude difference between the antenna array elements, the performance of the DOA estimation algorithm degrades. In addition, since the solution of the LLS method is obtained in the process of minimizing the amplitude and phase difference at the same time, the phase information of the interpolated antenna array elements can not be accurately represented systematically. These points lead to poor DOA estimation performance.

향상된 안테나 어레이 보간 방법에 관한 연구는 거의 없다. 균일한 원형 안테나 어레이에서 보간된 안테나 어레이 요소들을 생성하기 위해 테일러(Taylor) 시리즈 근사법을 사용하는 방법이 알려져 있다(아래 비특허문헌 참조). 그 방법은 더 나은 DOA 추정 성능을 달성할 수 있지만, 시리즈의 차수가 안테나 어레이 요소들의 최대 수보다 작은 것으로 제한되는 단점을 갖는다. 그렇기 때문에, 그 방법은 작은 개수의 안테나 어레이 요소들(예를 들어, 4 개의 안테나 어레이 요소들)을 사용하는 자동차 레이더 시스템에 대한 근사화 성능은 보장하지 못한다. There are few studies on the improved antenna array interpolation method. A method using a Taylor series approximation to generate interpolated antenna array elements in a uniform circular antenna array is known (see non-patent literature below). The method can achieve better DOA estimation performance, but has the disadvantage that the order is limited to less than the maximum number of antenna array elements. As such, the method does not guarantee approximate performance for a car radar system using a small number of antenna array elements (e.g., four antenna array elements).

M.-Y. Cao, L. Huang, W.-X. Xie, and H. C. So, "Interpolation array technique for direction finding via Taylor series fitting," IEEE China Summit and International Conference on Signal and Information Processing (ChinaSIP), pp. 736-740, July 2015.M.-Y. Cao, L. Huang, W.-X. Xie, and H. C. So, "Interpolation Array Technique for Direction Finding via Taylor series fitting," IEEE China Summit and International Conference on Signal and Information Processing (ChinaSIP), pp. 736-740, July 2015.

안테나 시스템을 통하여 전방에 위치한 타겟을 추정할 때, 두 타겟이 매우 인접해 있는 경우, 기존의 안테나 어레이를 통한 타겟 각도 추정 알고리즘으로는 그 두 타겟들을 구분하는 것이 힘들 수 있다.  When estimating a target located ahead through the antenna system, if the two targets are very close together, it may be difficult to distinguish between the two targets by a target angle estimation algorithm over existing antenna arrays.

따라서 본 발명의 일 목적은 원래의 안테나 어레이 요소들을 로그 영역에서 원하는 위치로 보간하여 그 보간된 안테나 어레이 요소들에 대한 수신신호들의 위상을 정확히 생성하여 DOA 추정의 성능을 향상시키기 위한 로그-영역 안테나 어레이 보간 기반 수신신호 DOA 추정 방법과 이를 위한 장치를 제공하는 것이다.It is therefore an object of the present invention to provide a log-area antenna for interpolating original antenna array elements from a logarithmic domain to a desired location to accurately generate the phases of the received signals for the interpolated antenna array elements, And to provide a method for estimating received signal DOA based on array interpolation and an apparatus therefor.

본 발명의 다른 목적은 로그 영역에서 보간된 안테나 어레이 요소들의 수신신호들 간의 전력 편차를 최소화하는 보정을 함으로써 수신신호의 DOA 추정의 정확성을 높일 수 있는 로그-영역 안테나 어레이 보간 기반 수신신호 DOA 추정 방법과 이를 위한 장치를 제공하는 것이다.Another object of the present invention is to provide a method and apparatus for estimating received signal DOA based on a log-area antenna array interpolation method capable of enhancing the accuracy of DOA estimation of a received signal by correcting to minimize a power deviation between received signals of antenna array elements interpolated in a logarithmic region And an apparatus therefor.

또한, 본 발명의 또 다른 목적은 위와 같은 DOA 추정 방법을 컴퓨터에서 실행하기 위한 컴퓨터 프로그램이 기록된 컴퓨터 판독 가능 기록 매체를 제공하는 것이다.It is still another object of the present invention to provide a computer readable recording medium on which a computer program for executing the DOA estimation method as described above is recorded.

본 발명의 목적은 상술한 것들에 한정되는 것이 아니며, 본 발명의 사상 및 영역으로부터 벗어나지 않는 범위에서 다양하게 확장될 수 있을 것이다.The object of the present invention is not limited to the above-described embodiments, but may be variously modified without departing from the spirit and scope of the present invention.

타겟의 각도 정보는 대부분 수신신호의 위상에 담겨 있다. 따라서 수신신호의 위상 정보에 초점을 맞추어서 수신신호의 각도를 추정할 수 있다. 이런 관점에서, 본 발명의 상기 목적을 달성하기 위한 실시예들에 따른 로그-영역 안테나 어레이 보간에 기반한 수신신호의 DOA 추정 방법은, 균일 선형 어레이(ULA) 형태로 배열된 안테나 어레이 요소들에 있어서, 원래의 안테나 어레이 요소들의 제1 스티어링 행렬과, 상기 원래의 안테나 어레이 요소들에 기초하여 보간한 보간 안테나 어레이 요소들의 제2 스티어링 행렬 각각에 대해 로그를 취하는 단계를 포함한다. 상기 DOA 추정 방법은 로그를 취한 상기 제1 스티어링 행렬을 로그를 취한 상기 제2 스티어링 행렬로 변환하는 변환행렬을 최적화 기법을 이용하여 생성하는 단계를 포함한다. 또한, 상기 DOA 추정 방법은 상기 원래의 안테나 어레이 요소들의 수신신호들의 위상에 대하여 상기 변환행렬을 적용하여 상기 보간 안테나 어레이 요소들에 대응되는 위상을 갖는 보간 수신신호들로 변환하는 단계를 포함한다. 나아가, 상기 DOA 추정 방법은 소정의 DOA 추정 알고리즘에 기초하여 상기 보간 안테나 어레이 요소들의 상기 보간 수신신호들의 DOA를 추정하는 단계를 포함한다. The target angle information is mostly contained in the phase of the received signal. Therefore, the angle of the received signal can be estimated by focusing on the phase information of the received signal. In this regard, a DOA estimation method of a received signal based on log-area antenna array interpolation according to embodiments of the present invention to achieve the above object, in an antenna array element arranged in a uniform linear array (ULA) And taking a log for each of a first steering matrix of the original antenna array elements and a second steering matrix of interpolated antenna array elements based on the original antenna array elements. The DOA estimation method includes a step of generating a transformation matrix for converting the first steering matrix obtained by taking the log into the second steering matrix obtained by using the optimization technique. In addition, the DOA estimation method may include applying the transformation matrix to the phases of the received signals of the original antenna array elements to convert them into interpolated reception signals having phases corresponding to the interpolation antenna array elements. Further, the DOA estimation method includes estimating a DOA of the interpolated received signals of the interpolation antenna array elements based on a predetermined DOA estimation algorithm.

예시적인 실시예들에 있어서, 상기 최적화 기법은 상기 수신신호들 간의 관계에 기초하여 상기 원래의 안테나 어레이 요소들과 상기 보간 안테나 어레이 요소들 간의 위상차를 최소화 하는 선형 최소 제곱법(Linear Least Square: LLS)일 수 있다.In exemplary embodiments, the optimization scheme may include a linear least squares (LLS) method that minimizes the phase difference between the original antenna array elements and the interpolation antenna array elements based on the relationship between the received signals. ).

예시적인 실시예들에 있어서, 상기 변환행렬은 상기 수신신호들의 진폭은 그대로 유지하면서 상기 수신신호들의 위상만을 로그 영역에서 조합하여 상기 보간 수신신호의 위상으로 변환할 수 있다.In the exemplary embodiments, the transformation matrix may combine only the phases of the received signals in the log domain to convert the phase of the received signals to the phase of the received signal while maintaining the amplitudes of the received signals intact.

예시적인 실시예들에 있어서, 상기 DOA 추정 방법은, 상기 안테나 어레이 요소들에 관해 설정된 관측 시야에 맞는 상기 원래의 안테나 어레이 요소들의 상기 제1 스티어링 행렬과, 상기 원래의 안테나 어레이 요소들을 선형 조합하여 원하는 위치로 보간한 상기 보간 안테나 어레이 요소들의 상기 제2 스티어링 행렬을 구하는 단계를 더 포함할 수 있다.In exemplary embodiments, the DOA estimation method further comprises: linearly combining the original antenna array elements with the first steering matrix of the original antenna array elements matching the observation field set with respect to the antenna array elements And obtaining the second steering matrix of the interpolation antenna array elements interpolated to the desired position.

예시적인 실시예들에 있어서, 상기 DOA 추정 방법은 상기 보간 수신신호들의 도래각(DOA) 추정의 정확도를 높이기 위해, 상기 보간 안테나 어레이 요소들의 상기 보간 수신신호들 간의 전력 레벨의 편차를 보상하여 상기 보간 수신신호들의 전력 레벨을 실질적으로 균등하게 보정하는 단계를 더 포함할 수 있다.In exemplary embodiments, the DOA estimation method comprises compensating for a deviation in the power level between the interpolated received signals of the interpolation antenna array elements to increase the accuracy of the DOA estimation of the interpolated received signals, And correcting the power level of the interpolation received signals substantially equally.

예시적인 실시예들에 있어서, 상기 보간 수신신호들의 위상은 상기 전력 레벨을 보정하기 전후로 동일하게 유지될 수 있다.In exemplary embodiments, the phase of the interpolated received signals may remain the same before and after correcting the power level.

예시적인 실시예들에 있어서, 상기 로그를 취하는 단계는, 상기 원래의 안테나 어레이 요소들의 관측 시야를 설정하는 단계와, 설정된 관측 시야에 맞는 상기 원래의 안테나 어레이 요소들의 상기 제1 스티어링 행렬과 상기 보간 안테나 어레이 요소들의 상기 제2 스티어링 행렬을 생성하는 단계와, 상기 제1 및 제2 스티어링 행렬의 각 성분에 대하여 로그를 취하는 단계를 포함할 수 있다.In the exemplary embodiments, the step of taking the log further comprises the steps of: setting an observation field of view of the original antenna array elements; interpolating the first steering matrix of the original antenna array elements Generating the second steering matrix of antenna array elements, and taking a log for each component of the first and second steering matrices.

본 발명의 실시예들에 따르면, 앞에서 언급한 수신신호의 DOA 추정 방법을 컴퓨터에서 실행하기 위한 컴퓨터 프로그램이 기록된 컴퓨터 가독형 기록 매체가 제공된다.According to embodiments of the present invention, a computer readable recording medium on which a computer program for executing the above-mentioned DOA estimation method of a received signal on a computer is recorded.

한편, 본 발명의 상기 목적들을 달성하기 위한 실시예들에 따른 로그-영역 안테나 어레이 보간 기반 수신신호의 도래각 추정 장치는 수신 안테나부, 수신부, 그리고 신호처리부를 포함한다. 상기 수신 안테나부는 일렬로 등간격으로 배열되고 전방에서 입사되는 무선신호를 수신하는 복수 개의 원래의 안테나 어레이 요소들을 포함한다. 상기 수신부는 상기 수신 안테나부의 각 안테나를 통해 수신된 복수 개의 수신신호로부터 소정의 신호를 추출하여 복수 개의 디지털 수신신호로 변환한다. 상기 신호처리부는 원래의 안테나 어레이 요소들의 제1 스티어링 행렬과, 상기 원래의 안테나 어레이 요소들에 기초하여 보간한 보간 안테나 어레이 요소들의 제2 스티어링 행렬 각각에 대해 로그를 취하는 기능과, 로그를 취한 상기 제1 스티어링 행렬을 로그를 취한 상기 제2 스티어링 행렬로 변환하는 변환행렬을 최적화 기법을 이용하여 생성하는 기능과, 상기 수신부가 제공하는 상기 복수 개의 디지털 수신신호들의 위상에 대하여 상기 변환행렬을 적용하여 상기 보간 안테나 어레이 요소들에 대응되는 위상을 갖는 보간 수신신호들로 변환하는 기능과, 소정의 도래각(DOA) 추정 알고리즘에 기초하여 상기 보간 안테나 어레이 요소들의 상기 보간 수신신호들의 도래각(DOA)을 추정하는 기능을 포함한다. According to another aspect of the present invention, there is provided an apparatus for estimating the arrival angle of a received signal based on log-area antenna array interpolation, the apparatus including a receiving antenna unit, a receiving unit, and a signal processing unit. The receive antenna portion includes a plurality of original antenna array elements arranged at equal intervals in a row and receiving a radio signal incident from the front side. The receiving unit extracts a predetermined signal from a plurality of reception signals received through the antennas of the reception antenna unit and converts the predetermined signals into a plurality of digital reception signals. Wherein the signal processing unit has a function of taking a log for each of a first steering matrix of original antenna array elements and a second steering matrix of interpolated antenna array elements based on the original antenna array elements, A function of generating a transformation matrix for transforming the first steering matrix into the logarithm of the second steering matrix using an optimization technique and a function of applying the transformation matrix to the phases of the plurality of digital received signals provided by the receiver (DOA) of the interpolated received signals of the interpolation antenna array elements based on a predetermined arrival angle (DOA) estimation algorithm, and a function to convert the interpolated received signals into an interpolated received signal having a phase corresponding to the interpolated antenna array elements. As shown in FIG.

예시적인 실시예들에 있어서, 상기 최적화 기법은 상기 수신신호들 간의 관계에 기초하여 상기 원래의 안테나 어레이 요소들과 상기 보간 안테나 어레이 요소들 간의 위상차를 최소화 하는 선형 최소 제곱법(Linear Least Square: LLS)일 수 있다. In exemplary embodiments, the optimization scheme may include a linear least squares (LLS) method that minimizes the phase difference between the original antenna array elements and the interpolation antenna array elements based on the relationship between the received signals. ).

예시적인 실시예들에 있어서, 상기 변환행렬은 상기 수신신호들의 진폭은 그대로 유지하면서 상기 수신신호들의 위상만을 로그 영역에서 조합하여 상기 보간 수신신호의 위상으로 변환할 수 있다.In the exemplary embodiments, the transformation matrix may combine only the phases of the received signals in the log domain to convert the phase of the received signals to the phase of the received signal while maintaining the amplitudes of the received signals intact.

예시적인 실시예들에 있어서, 상기 신호처리부는 상기 안테나 어레이 요소들에 관해 설정된 관측 시야에 맞는 상기 원래의 안테나 어레이 요소들의 상기 제1 스티어링 행렬과, 상기 원래의 안테나 어레이 요소들을 선형 조합하여 원하는 위치로 보간한 상기 보간 안테나 어레이 요소들의 상기 제2 스티어링 행렬을 구하는 기능을 더 포함할 수 있다.In exemplary embodiments, the signal processing unit may be configured to linearly combine the original antenna array elements with the first steering matrix of the original antenna array elements matching the viewing field set with respect to the antenna array elements, To obtain the second steering matrix of the interpolation antenna array elements interpolated by the interpolation antenna array elements.

예시적인 실시예들에 있어서, 상기 신호처리부는, 상기 보간 수신신호들의 도래각(DOA) 추정의 정확도를 높이기 위해, 상기 보간 안테나 어레이 요소들의 상기 보간 수신신호들 간의 전력 레벨의 편차를 보상하여 상기 보간 수신신호들의 전력 레벨을 실질적으로 균등하게 보정하는 기능을 더 포함할 수 있다.In exemplary embodiments, the signal processing unit may compensate for a deviation in the power level between the interpolated received signals of the interpolation antenna array elements to increase the accuracy of the DOA estimation of the interpolated received signals, And a function of substantially evenly correcting the power level of the interpolation received signals.

예시적인 실시예들에 있어서, 상기 보간 수신신호들의 위상은 상기 전력 레벨을 보정하기 전후로 동일하게 유지될 수 있다.In exemplary embodiments, the phase of the interpolated received signals may remain the same before and after correcting the power level.

본 발명은 DOA 추정의 정확도를 높이기 위해 어레이 보간에 사용되는 새로운 로그 영역 변환행렬을 제공한다. 본 발명에 따른 변환행렬은 원래의 안테나 어레이 요소와 보간된 안테나 어레이 요소들 간의 위상 차이를 최소화하여, 종래의 방법에 비해 보간 오차를 줄여 훨씬 더 정교한 위상으로 보간할 수 있다. 보간된 안테나 어레이 요소들에 대해 정확한 위상을 생성할 수 있다. 이러한 변환행렬을 이용하여 수신신호의 DOA를 추정하면, 위상 왜곡을 최소화 하여 DOA 추정 성능을 크게 개선할 수 있다. 따라서 본 발명에 따른 방법이 적용된 차량용 레이더 센서에 적용되면, 전방 표적 차량으로부터 반사되어 수신되는 수신신호에 대해 보다 향상된 각도 분해능을 가질 수 있다.The present invention provides a new log domain transformation matrix used in array interpolation to increase the accuracy of DOA estimation. The transformation matrix according to the present invention minimizes the phase difference between the original antenna array element and the interpolated antenna array elements, thereby interpolating into a much more sophisticated phase by reducing the interpolation error compared to the conventional method. It is possible to generate an accurate phase for the interpolated antenna array elements. Estimating the DOA of the received signal using these transform matrices minimizes the phase distortion and can greatly improve the DOA estimation performance. Therefore, when applied to a radar sensor for a vehicle to which the method according to the present invention is applied, it is possible to have an improved angular resolution for a received signal reflected and received from a front target vehicle.

또한, 본 발명의 실시예에 따른 안테나 어레이 변환은 보간된 안테나 어레이 요소들의 진폭에 영향을 미치지 않으며 변형 후에도 보존된다. In addition, the antenna array conversion according to embodiments of the present invention does not affect the amplitude of the interpolated antenna array elements and is preserved after deformation.

이러한 변환행렬은 사전에 계산되어 저장되므로 (오프라인 상태에서) 실시간으로 계산할 필요가 없다. Since these transformation matrices are calculated and stored in advance, they do not need to be computed in real time (off-line).

또한, 본 발명에 의해 로그 영역에서의 안테나 어레이 요소들의 보간을 위한 변환행렬을 사용할 때, 보간된 수신신호의 전력은 모든 안테나 어레이 요소들에 걸쳐 균일하지 않다. 이 경우, 안테나 어레이 보간 효과 및 DOA 추정의 성능은 완전하지 않을 수 있다. 보간된 안테나 어레이 요소들의 수신신호들 간의 전력 편차가 존재하는 경우, 그 수신신호들 간의 전력 편차를 줄여 보간된 수신신호의 전력을 비슷한 수준으로 조정하면 DOA 추정 알고리즘의 정확도를 더 높일 수 있다. Further, when using the transformation matrix for interpolation of antenna array elements in the log domain by the present invention, the power of the interpolated received signal is not uniform across all antenna array elements. In this case, the antenna array interpolation effect and DOA estimation performance may not be complete. If there is a power deviation between the received signals of the interpolated antenna array elements, adjusting the power of the interpolated received signal to a similar level by reducing the power deviation between the received signals can further improve the accuracy of the DOA estimation algorithm.

시뮬레이션 결과로부터 본 발명에 따른 로그 영역에서의 안테나 어레이 보간을 이용한 DOA 추정 방법이 기존의 안테나 어레이 보간법을 이용한 DOA 추정보다 우수한 DOA 추정 성능을 보이는 것을 확인할 수 있었다. 또한 차량용 레이더로 얻은 실제 측정 데이터를 통해 향상된 각도 분해능 및 추정 성능을 보여주는 것을 확인할 수 있었다. From the simulation results, it can be confirmed that the DOA estimation method using the antenna array interpolation in the log domain according to the present invention has better DOA estimation performance than the DOA estimation using the conventional antenna array interpolation method. In addition, the actual measurement data obtained from the vehicle radar shows the improved angular resolution and estimated performance.

도 1은 N개의(단, N은 2이상의 자연수) 안테나가 ULA의 형태로 배열된 ULA 안테나의 무선신호 수신을 도시한다.
도 2는 본 발명의 예시적인 실시예에 따른 도래각(DOA)의 추정 방법을 나타내는 흐름도이다.
도 3은 도 2의 흐름도의 단계 S100의 구체적인 실행 절차를 나타낸다.
도 4는 안테나 어레이 보간 방법을 설명하기 위한 도면이다.
도 5는 도 2의 흐름도의 단계 S300의 구체적인 실행 절차를 나타낸다.
도 6은 본 발명에 따른 방법을 실시하기 위한 시스템의 구성을 예시한다.
도 7은 도 6의 시스템의 DSP에서 수신신호로부터 도래각을 추정하는 처리를 수행하는 프로그램의 구성을 나타내는 블록도이다.
도 8은 본 발명의 알고리즘을 차량용 레이더 장치에 적용하여 전방의 2대의 차량을 인식하는 상황을 나타낸다.
도 9는 4개의 안테나 어레이 요소를 사용하여 본 발명에 따른 보간을 하여 DOA를 추정하는 방법의 성능에 관한 시뮬레이션을 하는 것을 설명하기 위한 도면이다.
도 10은 두 변환행렬

Figure 112017066647630-pat00001
Figure 112017066647630-pat00002
에 관해 FOV의 구간 크기에 따른 보간오차를 나타낸다.
도 11은 [θ1, θ2] = [-3.5˚, 2.5˚] 의 두 방향에 위치한 두 개의 인접한 표적에 관해 정규화한 바틀렛 의사스펙트럼을 예시적으로 나타낸 도면이다.
도 12의 (a)는 SNR에 따른 해상도 확률 P r 을 나타내며, (b)는 SNR에 따른 평균 제곱근 오차(Root Mean Square Error: RMSE)를 나타낸다.
도 13의 (a)는 시간 샘플 개수에 따른 해상도 확률을 나타내며, (b)는 시간 샘플 개수에 따른 RMSE를 나타낸다.
도 14는 [θ1, θ2, θ3] = [-8˚, 1.5˚, 7.5˚]의 세 방향에 위치한 세 개의 인접한 표적들에 관해 정규화한 바틀렛 의사스펙트럼을 예시적으로 나타낸다.
도 15의 (a)는 SNR에 따른 해상도 확률을 나타내고, (b)는 SNR에 따른 RMSE를 나타낸다(N=3인 경우).
도 16의 (a)는 SNR에 따른 해상도 확률을 나타내고, (b)는 SNR에 따른 RMSE를 나타낸다(N=5 인 경우).1 shows radio signal reception of a ULA antenna in which N (where N is a natural number of 2 or more) antennas are arranged in the form of a ULA.
2 is a flow chart illustrating a method of estimating an angle of incidence (DOA) in accordance with an exemplary embodiment of the present invention.
FIG. 3 shows a concrete execution procedure of step S100 of the flowchart of FIG.
4 is a diagram for explaining an antenna array interpolation method.
FIG. 5 shows a concrete execution procedure of step S300 of the flowchart of FIG.
Figure 6 illustrates the configuration of a system for implementing the method according to the present invention.
7 is a block diagram showing a configuration of a program for performing processing for estimating an incoming angle from a received signal in the DSP of the system of Fig.
8 shows a situation in which two vehicles ahead are recognized by applying the algorithm of the present invention to a radar system for a vehicle.
9 is a diagram for explaining simulation of performance of a method of estimating a DOA by performing interpolation according to the present invention using four antenna array elements.
FIG.
Figure 112017066647630-pat00001
Wow
Figure 112017066647630-pat00002
The interpolation error according to the segment size of the FOV.
Figure 11 is an illustration of a normalized Bartlet pseudo-spectrum for two adjacent targets located in two directions [θ 1 , θ 2 ] = [-3.5 °, 2.5 °].
12 (a) shows the resolution probability P r according to the SNR, and FIG. 12 (b) shows the root mean square error (RMSE) according to the SNR.
13 (a) shows the resolution probability according to the number of time samples, and (b) shows the RMSE according to the number of time samples.
14 illustrates an exemplary normalized Bartlet pseudo-spectrum for three adjacent targets located in three directions [θ 1 , θ 2 , θ 3 ] = [-8 °, 1.5 °, 7.5 °].
15 (a) shows the resolution probability according to the SNR, and (b) shows the RMSE according to the SNR (when N = 3).
16 (a) shows the resolution probability according to the SNR, and (b) shows the RMSE according to the SNR (when N = 5).

본문에 개시되어 있는 본 발명의 실시예들에 대해서, 특정한 구조적 내지 기능적 설명들은 단지 본 발명의 실시예를 설명하기 위한 목적으로 예시된 것으로, 본 발명의 실시예들은 다양한 형태로 실시될 수 있으며 본문에 설명된 실시예들에 한정되는 것으로 해석되어서는 아니 된다.For the embodiments of the invention disclosed herein, specific structural and functional descriptions are set forth for the purpose of describing an embodiment of the invention only, and it is to be understood that the embodiments of the invention may be practiced in various forms, The present invention should not be construed as limited to the embodiments described in Figs.

본 발명은 다양한 변경을 가할 수 있고 여러 가지 형태를 가질 수 있는 바, 특정 실시예들을 도면에 예시하고 본문에 상세하게 설명하고자 한다. 그러나 이는 본 발명을 특정한 개시 형태에 대해 한정하려는 것이 아니며, 본 발명의 사상 및 기술 범위에 포함되는 모든 변경, 균등물 내지 대체물을 포함하는 것으로 이해되어야 한다. 각 도면을 설명하면서 유사한 참조부호를 구성요소에 대해 사용하였다.The present invention is capable of various modifications and various forms, and specific embodiments are illustrated in the drawings and described in detail in the text. It is to be understood, however, that the invention is not intended to be limited to the particular forms disclosed, but on the contrary, is intended to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention. Similar reference numerals have been used for the components in describing each drawing.

제1, 제2 등의 용어는 다양한 구성요소들을 설명하는 데 사용될 수 있지만, 상기 구성요소들은 상기 용어들에 의해 한정되어서는 안 된다. 상기 용어들은 하나의 구성요소를 다른 구성요소로부터 구별하는 목적으로만 사용된다. The terms first, second, etc. may be used to describe various components, but the components should not be limited by the terms. The terms are used only for the purpose of distinguishing one component from another.

본 발명에서 사용한 용어는 단지 특정한 실시예를 설명하기 위해 사용된 것으로, 본 발명을 한정하려는 의도가 아니다. 단수의 표현은 문맥상 명백하게 다르게 뜻하지 않는 한, 복수의 표현을 포함한다. 본 출원에서 "포함하다" 또는 "가지다" 등의 용어는 명세서 상에 기재된 특징, 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것이 존재함을 지정하려는 것이지, 하나 또는 그 이상의 다른 특징들이나 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것들의 존재 또는 부가 가능성을 미리 배제하지 않는 것으로 이해되어야 한다.The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The singular expressions include plural expressions unless the context clearly dictates otherwise. It is to be understood that the terms "comprises" or "having" in this application are intended to specify the presence of stated features, integers, steps, operations, elements, parts, or combinations thereof, But do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, parts, or combinations thereof.

다르게 정의되지 않는 한, 기술적이거나 과학적인 용어를 포함해서 여기서 사용되는 모든 용어들은 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자에 의해 일반적으로 이해되는 것과 동일한 의미를 가지고 있다. 일반적으로 사용되는 사전에 정의되어 있는 것과 같은 용어들은 관련 기술의 문맥 상 가지는 의미와 일치하는 의미를 가지는 것으로 해석되어야 하며, 본 명세서에서 명백하게 정의하지 않는 한, 이상적이거나 과도하게 형식적인 의미로 해석되지 않는다. Unless defined otherwise, all terms used herein, including technical or scientific terms, have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Terms such as those defined in commonly used dictionaries are to be interpreted as having a meaning consistent with the contextual meaning of the related art and are to be interpreted as ideal or overly formal in the sense of the art unless explicitly defined herein Do not.

본문에 개시되어 있는 본 발명의 실시예들에 대해서, 특정한 구조적 내지 기능적 설명들은 단지 본 발명의 실시예를 설명하기 위한 목적으로 예시된 것으로, 본 발명의 실시예들은 다양한 형태로 실시될 수 있으며 본문에 설명된 실시예들에 한정되는 것으로 해석되어서는 안 된다. 즉, 본 발명은 다양한 변경을 가할 수 있고 여러 가지 형태를 가질 수 있는 바, 특정 실시예들을 도면에 예시하고 본문에 상세하게 설명하고자 한다. 그러나 이는 본 발명을 특정한 개시 형태에 대해 한정하려는 것이 아니며, 본 발명의 사상 및 기술 범위에 포함되는 모든 변경, 균등물 내지 대체물을 포함하는 것으로 이해되어야 한다. For the embodiments of the invention disclosed herein, specific structural and functional descriptions are set forth for the purpose of describing an embodiment of the invention only, and it is to be understood that the embodiments of the invention may be practiced in various forms, And should not be construed as limited to the embodiments described in the foregoing description. That is, the present invention is capable of various modifications and various forms, and specific embodiments are illustrated in the drawings and described in detail in the following description. It is to be understood, however, that the invention is not intended to be limited to the particular forms disclosed, but on the contrary, is intended to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention.

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

이하에서는 첨부한 도면을 참조하면서 본 발명의 실시예들에 관해 좀 더 상세히 설명하기로 한다.Hereinafter, embodiments of the present invention will be described in more detail with reference to the accompanying drawings.

도 1은 N개의 안테나가 일렬로 등간격 d 로 이격 배치된 ULA 안테나에서 전방의 표적으로부터 반사되어 되돌아오는 무선신호를 수신하는 것을 나타낸다. FIG. 1 shows that N antennas receive a radio signal reflected from a target ahead in a ULA antenna spaced apart by an equal distance d in a row.

도 1을 참조하면, ULA 형태의 안테나에서, θ 1 , θ 2 , ..., θ K 의 K개의 방향에서 나오는 무선 신호가 N개의 안테나가 일렬로 등간격 d 로 이격 배치된 ULA 안테나 요소들에 입사하는 경우를 고려한다. 안테나 어레이 요소들의 위치는 [d 1 , d 2 , ..., d N ]이며 안테나 어레이의 조준 방향으로부터 θ k (k = 1, 2, ..., K)가 정의된다. 원거리 협대역 신호원을 가정하면, 시간 t에서 실제 안테나 어레이에 수신된 신호 벡터 X (t)는 식 (1)과 같이 표현될 수 있다. In 1, in the form of a ULA antenna, θ 1, θ 2, ..., θ K the radio signals coming from the direction K of N antennas arranged in a ULA antenna elements spaced at equal intervals d in a line of And the incident light is considered. The positions of the antenna array elements are [ d 1 , d 2 , ..., d N ] and θ k (k = 1, 2, ..., K) is defined from the direction of sight of the antenna array. Assuming a distant narrowband signal source, the signal vector X (t) received at the actual antenna array at time t can be expressed as Equation (1).

Figure 112017066647630-pat00003
(1)
Figure 112017066647630-pat00003
(One)

여기서,

Figure 112017066647630-pat00004
는 벡터 전치연산자이고, A = [a(θ 1 ), a(θ 2 ), ..., a(θ K )]는 아래 식 (2)로 주어지는 스티어링 벡터(steering vector)이다.here,
Figure 112017066647630-pat00004
Is a vector permutation operator and A = [a ( ? 1 ), a ( ? 2 ), ..., a ( ? K )] is a steering vector given by the following equation (2).

Figure 112017066647630-pat00005
(2)
Figure 112017066647630-pat00005
(2)

여기서, K는 표적의 수, λ는 안테나의 수신신호의 파장, d i 는 첫 번째 안테나에서 i 번째 안테나까지의 거리를 나타낸다. 본 발명의 실시예에서는 ULA가 사용되기 때문에

Figure 112017066647630-pat00006
이고, 여기서, d 는 안테나 간격이다. Here, K is the number of the target, λ is the wavelength, d i of the received signal of the antenna is the distance from the first antenna to the i-th antenna. In the embodiment of the present invention, since ULA is used
Figure 112017066647630-pat00006
, Where d is the antenna spacing.

Figure 112017066647630-pat00007
는 시간 t에서의 Lx1 입사 신호 벡터를 나타낸다. N(t)는 Nx1 제로 평균 백색 가우시안 잡음 벡터(zero mean white Gaussian noise vector)를 나타낸다.
Figure 112017066647630-pat00008
는 ULA 안테나 어레이의 연장선과 직각을 이루는 직선과 각 안테나로 수신되는 입사 신호 벡터가 이루는 각을 나타낸다.
Figure 112017066647630-pat00007
Represents the Lx1 incident signal vector at time t. N (t) represents the Nx1 zero mean white Gaussian noise vector.
Figure 112017066647630-pat00008
Represents the angle formed by a straight line perpendicular to the extension line of the ULA antenna array and an incident signal vector received by each antenna.

도 2는 본 발명의 예시적인 실시예에 따른 로그 도메인 상에서의 보간에 기초한 수신신호의 도래각 추정 방법의 실행 절차를 나타낸다. 2 shows an execution procedure of a method of estimating the arrival angle of a received signal based on interpolation on the log domain according to an exemplary embodiment of the present invention.

우선, 로그 영역에서 변환행렬을 산출한다(S100 단계). 도 3은 예시적인 실시예에 따른 로그 영역에서의 변환행렬 산출 절차(S100 단계)를 구체적으로 나타낸다. 도 4는 자동차 레이더에 적용된 일반적인 안테나 어레이 보간 방법을 설명하기 위한 도면이다.First, a transformation matrix is calculated in the log area (S100). FIG. 3 specifically shows a transformation matrix calculation procedure (step S100) in the log domain according to the exemplary embodiment. 4 is a diagram for explaining a general antenna array interpolation method applied to an automobile radar.

먼저, 자동차 레이더의 관측 시야

Figure 112017066647630-pat00009
(field of view: FOV)를 설정한다(S110 단계). 자동차 레이더의 관측시야
Figure 112017066647630-pat00010
는 식 (3)으로 표현될 수 있다.First, the observation field of a car radar
Figure 112017066647630-pat00009
a field of view (FOV) is set (step S110). Observation field of a car radar
Figure 112017066647630-pat00010
Can be expressed by equation (3).

Figure 112017066647630-pat00011
(3)
Figure 112017066647630-pat00011
(3)

여기서, θL 및 θR은 도 4에 도시된 것처럼 FOV의 왼쪽 및 오른쪽 경계를 나타내는 각도이다. P는 각도 스템의 개수(number of angle step)를 나타낸다.Here,? L and? R are angles indicating the left and right boundaries of the FOV as shown in Fig. P represents the number of angle steps.

LLS 기법에서 유도된 변환행렬을 사용하는 기존의 안테나 어레이 보간 기법에 따르면, 원래의 안테나 어레이 요소들을 보간된 안테나 어레이 요소들로 변환하는 적절한 행렬을 찾는다. 이 범위에서 원래의 안테나 어레이 요소들의 스티어링 행렬

Figure 112017066647630-pat00012
은 식 (4)와 같이 나타낼 수 있다(S120 단계). 아래 식에서,
Figure 112017066647630-pat00013
은 각도 스텝의 크기(angle step size)이다. According to the conventional antenna array interpolation method using the transformation matrix derived from the LLS technique, an appropriate matrix for converting the original antenna array elements into the interpolated antenna array elements is found. In this range, the steering matrix of the original antenna array elements
Figure 112017066647630-pat00012
Can be expressed as Equation (4) (Step S120). In the equation below,
Figure 112017066647630-pat00013
Is the angle step size.

Figure 112017066647630-pat00014
(4)
Figure 112017066647630-pat00014
(4)

그런 다음 [g 1 , g 2 , ..., g M ] 위치의 안테나 어레이 요소들을 보간하면 보간된 안테나 어레이 요소들의 스티어링 행렬

Figure 112017066647630-pat00015
은 식 (5)와 같이 정해진다(S120 단계). Then, interpolating the antenna array elements at positions [ g 1 , g 2 , ..., g M ] results in a steering matrix of the interpolated antenna array elements
Figure 112017066647630-pat00015
Is determined as shown in equation (5) (step S120).

Figure 112017066647630-pat00016
(5)
Figure 112017066647630-pat00016
(5)

여기서,

Figure 112017066647630-pat00017
, (p = 1, 2, ..., P)이다. here,
Figure 112017066647630-pat00017
, ( p = 1, 2, ..., P).

원래의 스티어링 행렬

Figure 112017066647630-pat00018
을 보간된 스티어링 행렬
Figure 112017066647630-pat00019
로 변환하는 변환행렬
Figure 112017066647630-pat00020
가 있다고 가정하면, 그 변환행렬
Figure 112017066647630-pat00021
는 식 (6)과 같이 나타낼 수 있다.The original steering matrix
Figure 112017066647630-pat00018
The interpolated steering matrix < RTI ID =
Figure 112017066647630-pat00019
Conversion matrix
Figure 112017066647630-pat00020
, The transformation matrix
Figure 112017066647630-pat00021
Can be expressed as Equation (6).

Figure 112017066647630-pat00022
(6)
Figure 112017066647630-pat00022
(6)

적절한 변환행렬

Figure 112017066647630-pat00023
를 구하기 위해, 최소 자승법이 식 (7)과 같이 사용된다. 아래 식에서,
Figure 112017066647630-pat00024
는 Frobenius 행렬 표준을 나타낸다. Appropriate conversion matrix
Figure 112017066647630-pat00023
The least squares method is used as Equation (7). In the equation below,
Figure 112017066647630-pat00024
Represents the Frobenius matrix standard.

Figure 112017066647630-pat00025
(7)
Figure 112017066647630-pat00025
(7)

그런 다음, LLS 방법에 기초하여, 기존의 안테나 어레이 보간의 변환행렬

Figure 112017066647630-pat00026
는 식 (8)과 같이 구할 수 있다. 아래 식에서,
Figure 112017066647630-pat00027
Figure 112017066647630-pat00028
의 은닉 행렬(Hermitian matrix)을 나타낸다. Then, based on the LLS method, the transformation matrix of the existing antenna array interpolation
Figure 112017066647630-pat00026
Can be obtained as shown in equation (8). In the equation below,
Figure 112017066647630-pat00027
The
Figure 112017066647630-pat00028
(Hermitian matrix).

Figure 112017066647630-pat00029
(8)
Figure 112017066647630-pat00029
(8)

마지막으로 식 (6)과 식 (8)에서 보간된 스티어링 행렬

Figure 112017066647630-pat00030
의 추정은 식 (9)와 같이 주어진다.Finally, the steering matrix interpolated in Eqs. (6) and (8)
Figure 112017066647630-pat00030
Is given by Eq. (9).

Figure 112017066647630-pat00031
(9)
Figure 112017066647630-pat00031
(9)

또한, 행렬

Figure 112017066647630-pat00032
의 (m, p) 번째 구성요소는 식 (10)과 같이 나타낼 수 있다. 아래 식에서,
Figure 112017066647630-pat00033
Figure 112017066647630-pat00034
는 기존의 변환행렬
Figure 112017066647630-pat00035
와 원래의 스티어링 행렬
Figure 112017066647630-pat00036
의 (m, n) 번째와 (n, p) 번째 요소를 나타낸다.Also,
Figure 112017066647630-pat00032
The (m, p) th component of (10) can be expressed by Equation (10). In the equation below,
Figure 112017066647630-pat00033
and
Figure 112017066647630-pat00034
The conventional transformation matrix
Figure 112017066647630-pat00035
And the original steering matrix
Figure 112017066647630-pat00036
(M, n) -th element and (n, p)

Figure 112017066647630-pat00037
(10)
Figure 112017066647630-pat00037
(10)

이 변환행렬을 사용하여, 기존의 어레이 보간법으로 보간된 안테나 어레이 요소들의 수신신호

Figure 112017066647630-pat00038
들이 또한 생성될 수 있고, 이들은 식 (11)과 같이 주어진다. 기존에는 이러한 보간된 수신신호를 이용하여 DOA 추정을 수행하였다.Using this transformation matrix, the received signal of the antenna array elements interpolated by the existing array interpolation
Figure 112017066647630-pat00038
Can also be generated, and these are given by Equation (11). In the past, DOA estimation was performed using the interpolated received signal.

Figure 112017066647630-pat00039
(11)
Figure 112017066647630-pat00039
(11)

그런데, 안테나 어레이 보간을 위해 기존의 변환행렬

Figure 112017066647630-pat00040
을 사용할 때 몇 가지 중요한 문제가 발생한다. LLS 방법으로부터 얻어진 변환행렬에 기초하여, 보간된 안테나 어레이 요소들은 원래의 안테나 어레이 요소들의 선형 조합에 의해 생성된다. 이 경우 보간된 안테나 어레이 요소들의 진폭은 원래의 안테나 어레이 요소들의 진폭과 같지 않을 수 있다. 다시 말해, 원래의 안테나 어레이 요소들과 보간된 안테나 어레이 요소들이 수신한 신호의 진폭에 있어서 항상 일치되는 것은 아닐 수 있다. 각 안테나 어레이 요소들의 수신신호의 진폭이 안테나 어레이 전체에 걸쳐 균일하지 않은 경우, DOA 추정의 성능은 저하된다. 또한, 기존의 LLS 방법의 솔루션을 통해 보간된 안테나 어레이 요소들의 위상은 정확하게 생성되지 않을 수 있다. However, in order to interpolate the antenna array,
Figure 112017066647630-pat00040
There are some important problems when using. Based on the transformation matrix obtained from the LLS method, the interpolated antenna array elements are generated by a linear combination of the original antenna array elements. In this case, the amplitude of the interpolated antenna array elements may not be equal to the amplitude of the original antenna array elements. In other words, the original antenna array elements and the interpolated antenna array elements may not always coincide in amplitude with the received signal. If the amplitudes of the received signals of the respective antenna array elements are not uniform across the antenna array, the performance of the DOA estimation is degraded. In addition, the phase of the interpolated antenna array elements through the solution of the existing LLS method may not be precisely generated.

DOA 추정을 정확하게 하기 위해서는 보간된 안테나 어레이 요소들의 위상 정보가 중요하다. 따라서 안테나 어레이 요소들의 진폭을 동일하게 유지하면서, 원래의 안테나 어레이 요소들과 보간된 안테나 어레이 요소들 간의 위상차를 최소화 할 수 있는 보다 효과적인 안테나 어레이 보간 방법이 필요하다.The phase information of the interpolated antenna array elements is important for accurate DOA estimation. Therefore, there is a need for a more efficient antenna array interpolation method that can minimize the phase difference between the original antenna array elements and the interpolated antenna array elements while maintaining the same amplitude of the antenna array elements.

이를 위해, 본 발명의 예시적인 일 실시예에 따르면, 식 (6)의 양측에 있는 스티어링 행렬 즉, 원래의 안테나 어레이 요소들의 스티어링 행렬

Figure 112017066647630-pat00041
과 보간된 안테나 어레이 요소들의 스티어링 행렬
Figure 112017066647630-pat00042
의 원소들에 대해 식 (12)와 같은 로그를 취한다(S130 단계).To this end, according to an exemplary embodiment of the present invention, the steering matrix on both sides of Equation (6), i.e. the steering matrix of the original antenna array elements
Figure 112017066647630-pat00041
And the steering matrix of the interpolated antenna array elements
Figure 112017066647630-pat00042
(12) is obtained for the elements of < / RTI >

Figure 112017066647630-pat00043
(12)
Figure 112017066647630-pat00043
(12)

여기서 LOG(ㆍ)는 행렬의 각 원소에 대해 로그를 취하는 연산자이고,

Figure 112017066647630-pat00044
는 행렬
Figure 112017066647630-pat00045
의 (m, p) 번째 원소를 나타낸다.
Figure 112017066647630-pat00046
Figure 112017066647630-pat00047
행렬의 모든 원소들은 순수한 허수 값을 가진다. Here, LOG () is an operator that takes a log for each element of the matrix,
Figure 112017066647630-pat00044
The matrix
Figure 112017066647630-pat00045
(M, p) < th >
Figure 112017066647630-pat00046
Wow
Figure 112017066647630-pat00047
All elements of the matrix have pure imaginary values.

그런 다음, 로그를 취한 제1 스티어링 행렬

Figure 112017066647630-pat00048
를 로그를 취한 제2 스티어링 행렬
Figure 112017066647630-pat00049
로 변환하는 적절한 변환행렬
Figure 112017066647630-pat00050
를 찾아서 식 (13)으로 표현한다.Then, the first steering matrix obtained by taking the log
Figure 112017066647630-pat00048
The second steering matrix
Figure 112017066647630-pat00049
A suitable transformation matrix to transform
Figure 112017066647630-pat00050
And express it by Eq. (13).

Figure 112017066647630-pat00051
(13)
Figure 112017066647630-pat00051
(13)

원래의 도메인에서와 같이, 적절한 변환행렬

Figure 112017066647630-pat00052
는 최적화 기법을 사용하여 구할 수 있다. 그 최적화 기법의 대표적인 예로는 LLS 방법을 들 수 있다. 그 변환행렬
Figure 112017066647630-pat00053
는 식 (14)와 같이 주어진다(S140 단계).As in the original domain, the appropriate transformation matrix
Figure 112017066647630-pat00052
Can be obtained using optimization techniques. A typical example of the optimization technique is the LLS method. The transformation matrix
Figure 112017066647630-pat00053
Is given as Equation (14) (Step S140).

Figure 112017066647630-pat00054
......(14)
Figure 112017066647630-pat00054
(14)

이 변환행렬

Figure 112017066647630-pat00055
는 원래의 안테나 어레이 요소들의 위상들을 보간된 안테나 어레이 요소들의 위상들로 변환한다. 그런데 변환행렬
Figure 112017066647630-pat00056
은 로그 도메인에서 정의되기 때문에 식 (9)와 같이 원래의 안테나 어레이 요소들에 직접 적용될 수 없다. 달리 말하면, 이 변환행렬
Figure 112017066647630-pat00057
은 선형 연산자로 표현될 수는 없다. 그 대신, 원래의 안테나 어레이 요소들을 써서 식 (15)와 같이 나타낼 수 있다.This transformation matrix
Figure 112017066647630-pat00055
Converts the phases of the original antenna array elements into the phases of the interpolated antenna array elements. However,
Figure 112017066647630-pat00056
Can not be applied directly to the original antenna array elements as in Equation (9) because they are defined in the log domain. In other words,
Figure 112017066647630-pat00057
Can not be expressed as a linear operator. Instead, it can be expressed as Eq. (15) using the original antenna array elements.

Figure 112017066647630-pat00058
......(15)
Figure 112017066647630-pat00058
(15)

여기서

Figure 112017066647630-pat00059
는 새롭게 보간된 스티어링 매트릭스
Figure 112017066647630-pat00060
의 (m, p) 번째 요소이며,
Figure 112017066647630-pat00061
는 변환행렬
Figure 112017066647630-pat00062
의 (m, n) 번째 요소를 나타낸다.here
Figure 112017066647630-pat00059
Lt; RTI ID = 0.0 >
Figure 112017066647630-pat00060
(M, p) th element of < RTI ID = 0.0 >
Figure 112017066647630-pat00061
The transformation matrix
Figure 112017066647630-pat00062
(M, n) < th >

다시 도 2를 참조하면, 로그 영역에서 변환행렬을 구해놓은 상태에서, 수신 안테나 어레이 요소들이 레이더 신호를 시간 영역에서 수신할 수 있다(S200 단계). Referring again to FIG. 2, in a state where the transformation matrix is obtained in the logarithmic region, the reception antenna array elements can receive the radar signal in the time domain (step S200).

그리고 그 시간 도메인의 수신신호를 미리 산출해둔 변환행렬

Figure 112017066647630-pat00063
을 적용하여 로그 도메인에서 보간된 안테나의 어레이 요소들의 수신신호로 변환한다(S300 단계). 도 5는 이 단계를 좀 더 구체적으로 나타낸다.Then, a conversion matrix in which the reception signal of the time domain is calculated in advance
Figure 112017066647630-pat00063
(Step S300), and converts the received signal into a reception signal of the array elements of the interpolated antenna in the log domain. Figure 5 shows this step in more detail.

식 (8)의 기존의 변환행렬

Figure 112017066647630-pat00064
은 보간된 안테나 어레이 요소들을 원래의 안테나 어레이 요소들의 선형 조합으로 체계적으로 나타낸다(formulate). 그런데, 이 새로운 변환행렬
Figure 112017066647630-pat00065
은 원래의 안테나 어레이 요소들의 위상들을 조합하여 보간된 안테나 어레이 요소들의 위상 정보만을 생성한다(S310 단계). 다시 말하면, 본 발명에 따른 변환에 기초하여, 상기 새로운 변환행렬
Figure 112017066647630-pat00066
은 보간된 안테나 어레이 요소들 내의 원래의 안테나 어레이 요소들의 진폭들을 보존한다. 왜냐하면,
Figure 112017066647630-pat00067
의 크기가 식 (16)으로 나타내지기 때문이다.The existing transformation matrix of equation (8)
Figure 112017066647630-pat00064
Formulate the interpolated antenna array elements into a linear combination of the original antenna array elements. However, this new transformation matrix
Figure 112017066647630-pat00065
Generates only the phase information of the interpolated antenna array elements by combining the phases of the original antenna array elements (step S310). In other words, based on the transform according to the invention, the new transform matrix
Figure 112017066647630-pat00066
Preserves the amplitudes of the original antenna array elements in the interpolated antenna array elements. because,
Figure 112017066647630-pat00067
(16). ≪ / RTI >

Figure 112017066647630-pat00068
......(16)
Figure 112017066647630-pat00068
(16)

따라서 본 발명의 바람직한 실시예에 따른 어레이 변환은 보간된 안테나 어레이 요소들의 위상에만 영향을 주며, 보간된 안테나 어레이 요소들에 대해 더 정확한 위상을 생성한다. 기존의 변환행렬

Figure 112017066647630-pat00069
는 원래의 안테나 어레이 요소들의 진폭을 보존하지는 않는다. 왜냐하면
Figure 112017066647630-pat00070
가 항상 1(unity)이 아니기 때문이다. 따라서 새로운 변환행렬
Figure 112017066647630-pat00071
로부터 도출된 보간의 정확도는 기존의 변환행렬
Figure 112017066647630-pat00072
에로부터 도출된 보간의 정확도보다 높다고 할 수 있다.Thus, the array transformation according to the preferred embodiment of the present invention affects only the phase of the interpolated antenna array elements and produces a more accurate phase for the interpolated antenna array elements. Conventional transformation matrix
Figure 112017066647630-pat00069
Does not preserve the amplitude of the original antenna array elements. because
Figure 112017066647630-pat00070
Is always 1 (unity). Therefore,
Figure 112017066647630-pat00071
The accuracy of the interpolation derived from the conventional transformation matrix
Figure 112017066647630-pat00072
Is higher than the accuracy of the interpolation derived from Eq.

기존의 변환행렬

Figure 112017066647630-pat00073
을 통한 보간된 안테나 어레이 요소들의 수신신호가 식 (11)로 표현되는 것과 마찬가지로, 본 발명에 따른 새로운 변환행렬
Figure 112017066647630-pat00074
로 보간된 안테나 어레이 요소들의 수신신호
Figure 112017066647630-pat00075
는 식 (17)과 같이 표현될 수 있다(S320 단계).Conventional transformation matrix
Figure 112017066647630-pat00073
The received signal of the interpolated antenna array elements via the new transform matrix < RTI ID = 0.0 >
Figure 112017066647630-pat00074
≪ / RTI > of the antenna array elements < RTI ID =
Figure 112017066647630-pat00075
Can be expressed as equation (17) (step S320).

Figure 112017066647630-pat00076
......(17)
Figure 112017066647630-pat00076
(17)

이 보간된 수신신호 벡터

Figure 112017066647630-pat00077
를 사용하여 소정의 DOA 추정 알고리즘에 기초하여 DOA 추정을 수행할 수 있다(S400 단계). DOA 추정 알고리즘의 대표적인 예로는 바틀렛(Bartlett) 방법을 들 수 있지만, 그 밖의 다른 알려진 DOA 추정 알고리즘을 사용할 수 있음은 물론이다. This interpolated received signal vector
Figure 112017066647630-pat00077
The DOA estimation can be performed based on a predetermined DOA estimation algorithm (step S400). A typical example of the DOA estimation algorithm is a Bartlett method, but it is needless to say that other known DOA estimation algorithms can be used.

본 발명에 따라 로그 도메인 상에서의 보간 수신신호 벡터

Figure 112017066647630-pat00078
를 사용하여 DOA를 추정하는 것은 기존의 방법 즉, 시간 도메인 상에서의 보간을 통해 얻은 수신신호 벡터
Figure 112017066647630-pat00079
를 사용한 DOA 추정에 비해 훨씬 향상된 성능을 보인다. In accordance with the present invention, the interpolated received signal vector
Figure 112017066647630-pat00078
Is used to estimate the DOA by using the conventional method, that is, the received signal vector obtained through the interpolation in the time domain
Figure 112017066647630-pat00079
Which is much better than DOA estimation.

본 발명의 예시적인 실시예에 따르면, DOA 추정을 더 정확하게 하기 위해 수신신호의 전력도 고려할 수 있다. 보간된 수신신호 벡터

Figure 112017066647630-pat00080
Figure 112017066647630-pat00081
를 사용할 때, 그 보간된 수신신호들 간에는 전력의 편차가 존재할 수 있다. 즉, 아래 식 (18)은 항상 성립하지는 않을 수 있다.According to an exemplary embodiment of the present invention, the power of the received signal may also be considered to make the DOA estimation more accurate. The interpolated received signal vector
Figure 112017066647630-pat00080
Wow
Figure 112017066647630-pat00081
There may be a power variation between the interpolated received signals. That is, the following equation (18) may not always be established.

Figure 112017066647630-pat00082
......(18)
Figure 112017066647630-pat00082
(18)

왜냐하면, 아래 식 (19)가 성립하기 때문이다. 아래 식에서

Figure 112017066647630-pat00083
는 복소수의 복소 공액을 나타낸다. This is because the following equation (19) holds. In the equation below
Figure 112017066647630-pat00083
Represents a complex conjugate of a complex number.

Figure 112017066647630-pat00084
......(19)
Figure 112017066647630-pat00084
(19)

이러한 전력 불균일(power imbalance)은 DOA 추정에서 성능 저하를 일으킬 수 있다. 본 발명의 예시적인 실시예에 따른 방법에 의하면, 보간된 안테나 어레이 요소들

Figure 112017066647630-pat00085
의 진폭은 모든 안테나 어레이 요소들에 대해 동일하다. 그러나 이는 보간된 수신신호들의 전력과 직접 관련이 없으며, 보간된 수신신호들 간에는 전력 차가 존재한다. This power imbalance can cause performance degradation in DOA estimation. According to a method according to an exemplary embodiment of the present invention, interpolated antenna array elements
Figure 112017066647630-pat00085
Is the same for all antenna array elements. However, this is not directly related to the power of the interpolated received signals, and there is a power difference between the interpolated received signals.

그러므로 본 발명은 이러한 점을 고려하여, 보간된 수신신호들 간의 전력 차를 완화하기 위하여 단계 S320에서 새롭게 생성된 수신신호들의 진폭을 보정할 수 있다(S330 단계). 즉, 위에서 언급한 본 발명에 따른 위상 보간 방법의 효과를 유지하면서도, 보간된 안테나 어레이 요소들 각각의 수신신호들을 서로 거의 동일한 전력 레벨을 갖도록 체계적으로 나타내기(formulate) 위한 효과적인 보상 방법을 사용할 수 있다. 즉, 전력 레벨이 보상된 수신신호

Figure 112017066647630-pat00086
는 식 (20)과 같이 표현할 수 있다.Therefore, in order to mitigate the power difference between the interpolated received signals, the present invention may correct the amplitudes of the newly generated received signals in step S320 (step S330). That is, an effective compensation method can be used to systematically form the received signals of each of the interpolated antenna array elements to have substantially the same power level while maintaining the effect of the phase interpolation method according to the present invention described above have. That is, when the power level of the received signal
Figure 112017066647630-pat00086
Can be expressed as Equation (20).

Figure 112017066647630-pat00087
......(20)
Figure 112017066647630-pat00087
(20)

여기서,

Figure 112017066647630-pat00088
은 아래 식 (21)과 같다.here,
Figure 112017066647630-pat00088
Is given by the following equation (21).

Figure 112017066647630-pat00089
......(21)
Figure 112017066647630-pat00089
(21)

Figure 112017066647630-pat00090
은 집합
Figure 112017066647630-pat00091
의 집합원의 개수(cardinality) 즉, 보간된 안테나 어레이 요소들 각각을 생성하는 데 관여된 원래의 안테나 어레이 요소들의 개수를 나타낸다. 전력 레벨의 보상을 하지 않은 수신신호
Figure 112017066647630-pat00092
와 전력 레벨이 보상된 수신신호
Figure 112017066647630-pat00093
를 비교할 때,
Figure 112017066647630-pat00094
의 보간된 위상은
Figure 112017066647630-pat00095
의 보간된 위상과 같다. 그러므로 변환행렬
Figure 112017066647630-pat00096
로부터의 위상 보간 효과가 그대로 유지될 수 있다. 또한, 전력레벨 보상된 보간 수신신호를 사용할 때, 아래 식 (22)이 항상 성립한다. 왜냐하면
Figure 112017066647630-pat00097
의 크기의 제곱이 식 (23)과 같기 때문이다.
Figure 112017066647630-pat00090
Is a set
Figure 112017066647630-pat00091
Represents the cardinality of the original antenna array elements involved in generating each of the interpolated antenna array elements. Receive signal without power level compensation
Figure 112017066647630-pat00092
And the received signal whose power level is compensated
Figure 112017066647630-pat00093
When compared,
Figure 112017066647630-pat00094
The interpolated phase of
Figure 112017066647630-pat00095
≪ / RTI > Therefore,
Figure 112017066647630-pat00096
The phase interpolation effect can be maintained as it is. Further, when using the power level compensated interpolation received signal, the following equation (22) always holds. because
Figure 112017066647630-pat00097
Is the same as the formula (23).

Figure 112017066647630-pat00098
......(22)
Figure 112017066647630-pat00098
(22)

Figure 112017066647630-pat00099
......(23)
Figure 112017066647630-pat00099
(23)

즉, 전력레벨 보상된 보간 수신신호들의 전력

Figure 112017066647630-pat00100
은 안테나 어레이 요소들 간에 거의 균등하다. 따라서 DOA 추정을 위해 전력레벨 보상된 보간 수신신호 벡터
Figure 112017066647630-pat00101
가 사용되면, 시간 도메인 상에서의 보간 수신신호 벡터
Figure 112017066647630-pat00102
및 로그 도메인 상에서의 보간 수신신호 벡터
Figure 112017066647630-pat00103
를 사용할 때에 비해 더욱 향상된 성능을 달성할 수 있다.That is, the power of the power level compensated interpolation receiving signals
Figure 112017066647630-pat00100
Is approximately equal between the antenna array elements. Therefore, the power level-compensated interpolation received signal vector
Figure 112017066647630-pat00101
Is used, the interpolated received signal vector < RTI ID = 0.0 >
Figure 112017066647630-pat00102
And an interpolated received signal vector on the log domain
Figure 112017066647630-pat00103
The performance can be further improved as compared with the case of using.

한편, 위에서 설명한 본 발명의 실시예들에 따른 DOA 추정방법은 프로그램으로 구현하여 컴퓨터 가독형 불휘발성 메모리 장치 등과 같은 기록 매체에 저장해 두었다가, 필요 시 신호처리 프로세서 등과 같은 컴퓨팅 장치를 통해 그 프로그램을 실행하여 DOA를 산출하는 방식으로 실시할 수 있다. Meanwhile, the DOA estimation method according to embodiments of the present invention described above may be implemented as a program and stored in a recording medium such as a computer-readable nonvolatile memory device, and the program may be executed through a computing device such as a signal processing processor And calculating the DOA.

도 6은 본 발명에 따른 수신신호 도래각 추정 방법을 실시하기 위한 예시적인 실시예에 따른 ULA 안테나 기반 레이더 시스템(10)의 구성을 나타낸다. FIG. 6 shows a configuration of a ULA antenna-based radar system 10 according to an exemplary embodiment for implementing a method for estimating a received signal arrival angle according to the present invention.

도 6을 참조하면, 이 레이다 시스템(10)은 ULA 수신 안테나(60)와 ULA 송신 안테나(70)에 연결된 무선통신 모듈(RF module)(20), 패스밴드부(30), 디지털신호처리부(Digital Signal Processor: DSP, 40), 사용자 인터페이스(User Interface, 50)를 포함한 구성일 수 있다. 6, the radar system 10 includes a radio module 20, a passband unit 30, a digital signal processor (Digital Signal Processor) 30 connected to a ULA receive antenna 60 and a ULA transmit antenna 70, A signal processor (DSP) 40, and a user interface (50).

ULA 수신 안테나(60)와 ULA 송신 안테나(70)는 각각 다수의 안테나를 포함할 수 있다. 특히, ULA 수신 안테나(60)는 다수의 안테나가 등간격으로 일렬로 배치된 ULA 형태의 안테나 어레이를 가질 수 있다. 송신 안테나(70)를 통해 송출된 무선주파수 신호가 전방의 표적에서 반사되어 되돌아오는 반사 신호를 ULA 수신 안테나(60)가 수신할 수 있다.The ULA receive antenna 60 and the ULA transmit antenna 70 may each include a plurality of antennas. In particular, the ULA receiving antenna 60 may have a ULA type antenna array in which a plurality of antennas are arranged in a line at equal intervals. The ULA receiving antenna 60 can receive a reflected signal that is reflected from a target ahead of the radio frequency signal transmitted through the transmitting antenna 70 and returned.

무선통신모듈(20)은 ULA 수신 안테나(60)에 연결되어 ULA 수신 안테나(60)가 잡은 미약한 신호를 증폭시키는 저잡음 증폭기(Low-Noise Amplifier: LNA)(22)를 포함한다. 또한, 디지털 신호 처리부(40)가 제공하는 디지털 출력신호에 기초하여 소정의 아날로그 파형을 갖는 신호를 생성하는 파형발생기(Waveform Generator, 24), 생성된 신호를 무선 송출하기 위해 무선주파수(RF) 신호로 변환하는 발진기(26), 그리고 발진기(26)에서 출력되는 신호를 송출에 필요한 출력으로 증폭하여 송신 안테나(70)에 제공하는 전력증폭기(Power Amplifier: PA, 28)를 포함한다. 발진기(26)는 예컨대 전압제어발진기(Voltage Control Oscillator: VCO)로 구성될 수 있다.The wireless communication module 20 includes a Low-Noise Amplifier (LNA) 22 connected to the ULA receiving antenna 60 and amplifying a weak signal captured by the ULA receiving antenna 60. In addition, a waveform generator (Waveform Generator) 24 for generating a signal having a predetermined analog waveform based on the digital output signal provided by the digital signal processing unit 40, a radio frequency (RF) signal And a power amplifier (PA) 28 for amplifying a signal output from the oscillator 26 to an output necessary for transmission and providing the amplified signal to a transmission antenna 70. [ The oscillator 26 may be composed of, for example, a voltage controlled oscillator (VCO).

패스밴드부(30)는 저잡음 증폭기(22)에서 출력되는 수신신호와 발진기(26)에서 출력되는 발진신호를 혼합하는 혼합기(Mixer, 32)와, 이 혼합기(32)의 출력신호에 포함된 상기 수신신호와 상기 발진신호의 합 주파수를 제거하여 중간주파수(IF)의 신호출력을 추출하기 위한 저역통과필터(Low-Pass Filter: LPF, 34), 그리고 LPF(34)를 통해 얻어지는 중간주파수 신호를 디지털 신호로 변환하기 위한 아날로그-디지털 변환기(Analog-to-Digital Converter: ADC, 36)를 포함한다.The passband unit 30 includes a mixer 32 for mixing a received signal output from the low noise amplifier 22 and an oscillation signal output from the oscillator 26, Pass filter (LPF) 34 for extracting the signal output of the intermediate frequency (IF) by removing the sum frequency of the signal and the oscillation signal, and an intermediate frequency signal obtained through the LPF 34, And an analog-to-digital converter (ADC) 36 for converting the analog signal into a digital signal.

DSP(40)는 패스밴드부(30)에서 제공되는 디지털 중간주파수 신호를 위에서 설명한 방법에 따라 처리하여 안테나 어레이 수신신호의 DOA를 추정한다. 또한, 표적으로 송출할 신호를 생성하여 RF 모듈(20)에 제공한다.The DSP 40 processes the digital intermediate frequency signal provided by the passband unit 30 according to the method described above to estimate the DOA of the received signal of the antenna array. Further, a signal to be transmitted as a target is generated and provided to the RF module 20.

사용자 인터페이스(UI, 50)는 DSP(40)의 처리 결과를 표시하거나 또는 사용자의 지시를 DSP(40)에 전달한다. The user interface (UI) 50 displays the processing result of the DSP 40 or delivers the user's instruction to the DSP 40. [

예시적인 실시예에 따르면, 본 발명에 따른 수신신호의 DOA 추정방법은 프로그램으로 구현될 수 있다. 그 프로그램은 DSP(40)에 내장되어 실행될 수 있다. 도 7은 이 프로그램의 구성을 나타내는 블록도이다. According to an exemplary embodiment, a DOA estimation method of a received signal according to the present invention can be implemented by a program. The program can be embedded in the DSP 40 and executed. 7 is a block diagram showing the configuration of this program.

예시적인 실시예에 따르면, 이 프로그램은 위상 정보 생성부(310), 수신신호 생성부(320), 진폭 보정부(330), 그리고 DOA 추정부(400)를 포함할 수 있다.According to an exemplary embodiment, the program may include a phase information generator 310, a received signal generator 320, an amplitude corrector 330, and a DOA estimator 400.

위상 정보 생성부(310), 수신신호 생성부(320), 그리고 진폭 보정부(330)는 수신 안테나 소자(60)로 수신하여 디지털 신호로 변환된 시간 도메인의 수신신호를 미리 산출해둔 변환행렬

Figure 112017066647630-pat00104
을 적용하여 변환하는 단계 S300의 처리를 수행한다. 구체적으로, 위상 정보 생성부(310)는 패스밴드부(30)의 ADC(36)로부터 디지털화된 안테나 어레이 수신신호를 제공받아서, 미리 구해둔 변환행렬
Figure 112017066647630-pat00105
을 적용하여 위상 정보를 생성하는 단계 S310의 처리를 수행한다. 수신신호 생성부(320)는 위상 정보 생성부(310)가 구한 위상 정보를 갖는 수신신호를 생성하는 단계 S320의 처리를 수행한다. 그리고 진폭 보정부(330)는 수신신호 생성부(320)가 생성한 수신신호의 진폭을 보정하는 단계 S330의 처리를 수행한다.The phase information generation unit 310, the reception signal generation unit 320 and the amplitude correction unit 330 receive the reception signal from the reception antenna element 60,
Figure 112017066647630-pat00104
The process of the step S300 of applying the conversion is carried out. Specifically, the phase information generation unit 310 receives the digitized antenna array reception signal from the ADC 36 of the passband unit 30,
Figure 112017066647630-pat00105
And performs the process of step S310 of generating phase information. The reception signal generation unit 320 performs the process of step S320 of generating a reception signal having the phase information obtained by the phase information generation unit 310. [ The amplitude corrector 330 performs the process of step S330 for correcting the amplitude of the reception signal generated by the reception signal generator 320. [

DOA 추정부(400)는 진폭 보정부(330)를 통해 얻어지는 로그 도메인에서 보간된 안테나 어레이 요소들의 수신신호를 이용하여 DOA 추정 알고리즘을 수행한다. 이를 통해, 수신 안테나 소자(60)의 안테나 어레이 요소들에 수신되는 신호들의 도래각을 추정하는 단계 S400의 처리를 수행한다. The DOA estimation unit 400 performs a DOA estimation algorithm using the received signals of the antenna array elements interpolated in the log domain obtained through the amplitude correction unit 330. [ Thereby, the processing of step S400 of estimating the angle of arrival of the signals received at the antenna array elements of the receiving antenna element 60 is performed.

한편, 발명자들은 본 발명의 예시적인 실시예에 따른 보간법을 사용하여 얻어지는 DOA 추정의 정확도를 검증하기 위해, 몇 가지 시뮬레이션을 해보았다. 한 가지 시뮬레이션을 위해, 정확히 동일한 개구들(apertures)을 유지하면서, 원래의 안테나 어레이 요소들을 최소 중복(redundancy) 선형 안테나 어레이로 변환하였다. 일반적으로 최소 중복 선형 안테나 어레이는 안테나 어레이의 중복 공간 수를 최소화하여 주어진 수의 안테나 어레이 요소들에 대한 최대 해상도를 보여준다. 더욱이, 비균일 선형 어레이는 동일한 개구들을 갖는 균일한 선형 어레이보다 더 나은 DOA 추정 성능을 보인다고 알려져 있다. 따라서 시뮬레이션에서 원래의 안테나 선형 어레이로 변환함으로써, DOA 추정의 성능 향상을 분석하였다.Meanwhile, the inventors have conducted some simulations to verify the accuracy of the DOA estimation obtained using the interpolation method according to the exemplary embodiment of the present invention. For one simulation, the original antenna array elements were converted into a minimal redundant linear antenna array while maintaining exactly the same apertures. Generally, the least redundant linear antenna array minimizes the number of redundant spaces in the antenna array and thus shows the maximum resolution for a given number of antenna array elements. Moreover, non-uniform linear arrays are known to exhibit better DOA estimation performance than uniform linear arrays with identical apertures. Therefore, the performance improvement of the DOA estimation is analyzed by converting to the original antenna linear array in the simulation.

본 발명에 따른 로그-영역 어레이 보간에 기반한 레이더 수신신호의 DOA 추정 방법은 표적 검출용 레이더 장치에 응용하여 DOA를 고해상도로 추정할 수 있다. 도 8은 본 발명에 따른 레이더 수신신호의 DOA 추정 알고리즘을 차량(100)용 레이더 장치(10)에 적용하여 전방의 2대의 차량(200a, 200b)을 인식하는 상황을 나타낸다. The DOA estimation method of the radar received signal based on the log-area array interpolation according to the present invention can be applied to the radar device for target detection to estimate the DOA at high resolution. FIG. 8 shows a situation in which two vehicles 200a and 200b ahead are recognized by applying the DOA estimation algorithm of the radar received signal according to the present invention to the radar device 10 for the vehicle 100. FIG.

시뮬레이션에서는 도 9에 도시된 것처럼 자동차 장거리 레이더 (LRR)에서 널리 사용되는 4개의 안테나 어레이 요소(X1, X2, X3, X4) (N = 4)를 사용하였다. 원래의 안테나 어레이 요소들의 위치는 [d1, d2, d3, d4] = [0, 2λ, 4λ, 6λ]이고, 4개의 안테나 어레이 요소들의 최소 중복 선형 어레이 위치는 [0, 1, 4, 6]로 하였다. 그래서 도 9에 예시된 것처럼 안테나 어레이 변환행렬을 사용하여 [g1, g2, g3, g4] = [0, 1λ, 4λ, 6λ] 위치의 안테나 어레이 요소들을 보간하였다. 여기서, 두 대의 표적 차량(200a, 200b)은 [θ1, θ2] = [-3.5˚, 2.5˚] 방향에 있다고 가정한다. DOA 추정 알고리즘으로 바틀렛 알고리즘을 사용하였다. 또한, 안테나 어레이 요소들에서의 신호 대 잡음비(SNR)는 10dB로 설정되며, 수천 개의 시간 샘플을 사용하여 바틀렛 알고리즘에 사용된 상관 행렬을 구성하였다. 또한, FOV는 Θ = {θp | LRR의 FOV와 등가인 θp = -10˚ + (p-1) x 0.1˚, p = 1, 2, ..., 201}이다. 안테나 어레이 요소들 수와 FOV가 주어지면 기존 기술에 따른 변환행렬

Figure 112017066647630-pat00106
와 본 발명에 따른 변환행렬
Figure 112017066647630-pat00107
가 한 번 계산되어 저장되기 때문에, 다른
Figure 112017066647630-pat00108
Figure 112017066647630-pat00109
를 찾을 필요 없이 저장된 값을 반복적으로 사용할 수 있다.The four antenna array elements in the simulation as shown in Figure 9, which are widely used in automobile long-range radar (LRR) (X 1, X 2, X 3, X 4) (N = 4) were used. And the minimum redundant linear array position of the four antenna array elements is [0, 1, 2, 4, 6], and the positions of the original antenna array elements are [d 1 , d 2 , d 3 , d 4 ] = [ 4, 6]. Thus, the antenna array elements at positions [g 1 , g 2 , g 3 , g 4 ] = [0, 1λ, 4λ, 6λ] are interpolated using the antenna array transformation matrix as illustrated in FIG. Here, it is assumed that the two target vehicles 200a and 200b are in the [θ 1 , θ 2 ] = [-3.5 °, 2.5 °] direction. The Bartlett algorithm was used as the DOA estimation algorithm. In addition, the signal-to-noise ratio (SNR) in the antenna array elements was set to 10 dB and the correlation matrix used in the Bartlett algorithm was constructed using thousands of time samples. Also, the FOV can be expressed as? = {? P | Equivalent to FOV of LRR, θ p = -10 ° + (p-1) x 0.1 deg, p = 1, 2, ..., 201}. Given the number of antenna array elements and the FOV,
Figure 112017066647630-pat00106
And the transformation matrix according to the present invention
Figure 112017066647630-pat00107
Is calculated and stored once,
Figure 112017066647630-pat00108
Wow
Figure 112017066647630-pat00109
You can use the stored value repeatedly without having to look up the value.

먼저, 이러한 시뮬레이션 조건을 사용하여 식 (24)와 같이 두 가지 유형의 보간 오차를 계산한다.First, using these simulation conditions, two types of interpolation errors are calculated as in equation (24).

Figure 112017066647630-pat00110
(24)
Figure 112017066647630-pat00110
(24)

상기 식(24)에 기초하여 계산된 오차 값이 작을수록, 더욱 정확하게 안테나 어레이의 보간이 수행될 수 있다. 두 변환행렬

Figure 112017066647630-pat00111
Figure 112017066647630-pat00112
에 대해서, FOV의 구간 크기 변화에 따른 보간 오차를 계산할 수 있다. 그 결과는 도 10에서 확인할 수 있듯이 본 발명에 따른 변환행렬
Figure 112017066647630-pat00113
를 적용하였을 때의 보간 오차와 위상 보간 오차는 거의 0에 가깝다. 이에 비해, 기존 기술에 따른 변환행렬
Figure 112017066647630-pat00114
를 적용한 경우에는 보간 오차와 위상 보간 오차는 FOV가 커짐에 따라 큰 폭으로 증가한다.The smaller the error value calculated based on the above equation (24), the more accurately the interpolation of the antenna array can be performed. Two transformation matrices
Figure 112017066647630-pat00111
Wow
Figure 112017066647630-pat00112
The interpolation error according to the variation of the section size of the FOV can be calculated. As a result, as shown in FIG. 10, the transformation matrix
Figure 112017066647630-pat00113
The interpolation error and the phase interpolation error are close to zero. On the other hand,
Figure 112017066647630-pat00114
The interpolation error and the phase interpolation error increase greatly as the FOV increases.

또한 LRR의 FOV (즉, FOV의 길이는 20˚)에 대해 오차들은 다음과 같이 주어진다: For the FOV of the LRR (i.e., the length of the FOV is 20 degrees), the errors are given by:

Figure 112017066647630-pat00115
Figure 112017066647630-pat00116
.
Figure 112017066647630-pat00115
Wow
Figure 112017066647630-pat00116
.

그러므로 위의 두 가지 유형의 보간 오차의 결과로부터 판단해볼 때, 본 발명의 방법에 따른 어레이 변환행렬

Figure 112017066647630-pat00117
Figure 112017066647630-pat00118
에 비해
Figure 112017066647630-pat00119
와 더 유사하다. 다시 말하면, 본 발명에 따른 보간 방법을 채용할 때, 보간된 안테나 어레이 요소들이 더 정확하게 생성될 수 있다.Therefore, judging from the results of the above two types of interpolation errors, the array transformation matrix according to the method of the present invention
Figure 112017066647630-pat00117
this
Figure 112017066647630-pat00118
Compared to
Figure 112017066647630-pat00119
. In other words, when employing the interpolation method according to the present invention, the interpolated antenna array elements can be more accurately generated.

이러한 변환행렬을 사용하여 수신신호를 체계화하고 DOA 추정을 수행할 수 있다. 도 11은 [θ1, θ2] = [-3.5˚, 2.5˚]의 두 방향에 위치한 두 개의 인접한 표적들에 관해 정규화한 바틀렛 의사스펙트럼을 예시적으로 나타낸 도면이다. Such a transformation matrix can be used to organize the received signal and perform DOA estimation. Figure 11 is an illustration of a normalized Bartlet pseudo-spectrum for two adjacent targets located in two directions [θ 1 , θ 2 ] = [-3.5 °, 2.5 °].

도 11의 바틀렛 의사스펙트럼에 따르면, 바틀렛 방법은 원래의 수신신호에서 두 표적을 분간할 수 없으며 추정 DOA는 -0.1˚이다(파란색 그래프 참조). 일반적으로 2λ 간격으로 배치된 4개의 안테나 어레이 요소들을 사용하면, 반 전력 빔폭(half-power beamwidth)은 6.5˚가 된다. 그러므로 주어진 DOA를 구별할 수 없는 것이 어쩌면 합리적인 결과이다. 기존 기술에 따른 변환행렬

Figure 112017066647630-pat00120
을 적용하여 얻은 보간된 수신신호로도 두 개의 서로 다른 DOA가 추정되지는 않으며, 추정된 DOA는 1.2˚(정확한 값은 아님)이다(주황색 그래프 참조).According to the Bartlett pseudo-spectrum of FIG. 11, the Bartlett method can not distinguish between two targets in the original received signal and the estimated DOA is -0.1 (see blue graph). In general, using four antenna array elements arranged at 2 [lambda] spacings, the half-power beamwidth is 6.5 [deg.]. Therefore, it is perhaps a reasonable outcome to be unable to distinguish a given DOA. Conversion matrix according to existing technology
Figure 112017066647630-pat00120
, The estimated DOA is 1.2 ° (not the exact value) (see the orange graph).

이에 비해, 본 발명에 따른 변환행렬

Figure 112017066647630-pat00121
을 적용하여 얻은 보간된 수신신호들의 경우, Bartlett 방법은 향상된 각도 분해능을 보여주며, [-2.8˚, 2.0˚]와 같은 두 가지 다른 DOA를 얻게 해주는 것을 알 수 있다(노란색 그래프 참조). 또한, 전력레벨 보상된 보간 수신신호 벡터
Figure 112017066647630-pat00122
를 사용하면, 가장 좋은 DOA 추정 결과를 얻을 수 있으며, 추정된 DOA 값은 실제 DOA 값에 거의 근접한 [-3.1˚, 2.2˚]이다(보라색 그래프 참조). In contrast, the transformation matrix
Figure 112017066647630-pat00121
, The Bartlett method shows an improved angular resolution and gives two different DOAs such as [-2.8 °, 2.0 °] (see yellow graph). Also, the power level compensated interpolated received signal vector
Figure 112017066647630-pat00122
, The best DOA estimation result is obtained, and the estimated DOA value is [-3.1˚, 2.2˚] which is close to the actual DOA value (see purple graph).

통계적 성능 평가를 위해 기존의 Bartlett 방법과 본 발명에 따른 로그 영역 상의 어레이 보간을 적용한 Bartlett 방법의 해상도 확률 P r 을 계산하는 것도 의미가 있다. 또한 다른 시뮬레이션 조건을 그대로 유지하면서 안테나 어레이 SNR 값만을 변화시키면서 해상도 확률 P r 과 RMSE를 산출할 수 있다. 이 계산의 결과가 도 12에 표시되어 있다. 도 12는 SNR에 따른 해상도 확률 P r 을 나타내며, (b)는 SNR에 따른 평균 제곱근 오차(RMSE)를 나타낸다. 이 해상도 확률 P r 은 다음과 같이 정의될 수 있다.For the statistical performance evaluation, it is also meaningful to calculate the resolution probability P r of the Bartlett method and the Bartlett method applying the array interpolation on the log domain according to the present invention. Also, the resolution probability P r and RMSE can be calculated while changing only the antenna array SNR value while maintaining other simulation conditions. The results of this calculation are shown in Fig. 12 shows the resolution probability P r according to the SNR, and (b) shows the mean square root error (RMSE) according to the SNR. This resolution probability P r can be defined as follows.

Figure 112017066647630-pat00123
(25)
Figure 112017066647630-pat00123
(25)

여기서, N r 은 수신된 신호들로부터 두 개의 별개의 DOA들이 추출되는 횟수를 나타내고, N t 는 시뮬레이션들의 수를 나타낸다. 동일한 조건에서 이 시뮬레이션을 1000번 수행하므로, N t 는 1000이다. 아래 식 (26)과 같이 정의되는 평균 제곱근 오차(Root Mean Square Error: RMSE)를 계산한다.Where N r represents the number of times two separate DOAs are extracted from the received signals and N t represents the number of simulations. We perform this simulation 1000 times under the same conditions, so N t is 1000. A root mean square error (RMSE) defined by the following equation (26) is calculated.

Figure 112017066647630-pat00124
(26)
Figure 112017066647630-pat00124
(26)

여기서,

Figure 112017066647630-pat00125
는 q번째 (q = 1, 2, ..., Nt) 시뮬레이션에서 θ k (k = 1, 2)의 추정된 값이다. 결과는 표 1과 같다. here,
Figure 112017066647630-pat00125
The q-th (q = 1, 2, ... , Nt) from the simulation θ k (k = 1, 2). The results are shown in Table 1.

[-3.5˚, 2.5˚]에 위치한 두 인접 표적에 대한 해상도 확률 P r 및 RMSEThe resolution probabilities P r and RMSE for two adjacent targets located at [-3.5˚, 2.5˚] DOA 추정 DOA estimation P r (%) P r (%) RMSE (˚)RMSE (°) (a) 기존의 바틀렛 방법 수행 시(a) When carrying out the conventional Bartlett method 0.0 0.0 4.284.28 (b)

Figure 112017066647630-pat00126
Figure 112017066647630-pat00127
로 바틀렛 방법 수행 시(b)
Figure 112017066647630-pat00126
Wow
Figure 112017066647630-pat00127
When performing the Bartlett method 74.974.9 2.672.67 (c)
Figure 112017066647630-pat00128
Figure 112017066647630-pat00129
로 바틀렛 방법 수행 시
(c)
Figure 112017066647630-pat00128
Wow
Figure 112017066647630-pat00129
When performing the Bartlett method
99.499.4 0.640.64
(d)
Figure 112017066647630-pat00130
Figure 112017066647630-pat00131
로 바틀렛 방법 수행 시
(d)
Figure 112017066647630-pat00130
Wow
Figure 112017066647630-pat00131
When performing the Bartlett method
99.999.9 0.450.45

해상도 확률 P r 과 RMSE를 고려할 때, 본 발명에 따라 로그 영역 상의 어레이 보간을 적용한 Bartlett 방법(표 1의 케이스 (c)와 (d))은 기존의 단순한 바틀렛 방법(표 1의 케이스 (a) ) 또는 기존의 변환행렬

Figure 112017066647630-pat00132
를 적용하는 바틀렛 방법(표 1의 케이스 (b))에 비해 더 우수한 DOA 추정 성능을 보여줌을 알 수 있다. 도 12는 본 발명에 따른 방법은 안테나 어레이의 SNR 값이 달라지더라도 우수한 추정 결과를 내놓을 수 있음을 보여준다.Considering the resolution probability P r and RMSE, the Bartlett method (Case (c) and (d) in Table 1) applying the array interpolation on the logarithmic region according to the present invention is similar to the conventional simple Bartlett method (Case (a) ) Or an existing transformation matrix
Figure 112017066647630-pat00132
(Case (b) in Table 1), which is a better performance than the Bartlet method. FIG. 12 shows that the method according to the present invention can produce an excellent estimation result even if the SNR value of the antenna array changes.

또한, 상관행렬 구축에 사용된 샘플 수를 변경하면서 보간 방법들 간의 성능 비교를 수행한 결과가 도 13에 도시되어 있다. 도 13의 (a)는 시간 샘플 개수에 따른 해상도 확률을 나타내며, (b)는 시간 샘플 개수에 따른 RMSE를 나타낸다. 적은 수의 시간 샘플이 사용 되더라도 본 발명에 따른 안테나 어레이 변환방법은 향상된 예측 성능을 보여줌을 알 수 있다. Also, the results of performing performance comparisons between interpolation methods while changing the number of samples used in the construction of the correlation matrix are shown in FIG. 13 (a) shows the resolution probability according to the number of time samples, and (b) shows the RMSE according to the number of time samples. It can be seen that the antenna array conversion method according to the present invention shows improved prediction performance even when a small number of time samples are used.

나아가, 상기 레이더의 FOV 내에 3개의 표적 차량이 존재하는 경우에 대한 시뮬레이션을 수행해보았다. 이 시뮬레이션은 표적 정보를 제외하고는 도 10과 동일한 시뮬레이션 조건을 유지하면서 수행하였다. 도 14는 [θ1, θ2, θ3] = [-8˚, 1.5˚, 7.5˚]의 세 방향에 위치한 세 개의 인접한 표적들에 관해 정규화한 바틀렛 의사스펙트럼을 예시적으로 나타낸다. 도 14에 도시된 바와 같이 기존의 바틀렛 방법 (도 14의 파란색 그래프 참조) 및 기존의 변환행렬

Figure 112017066647630-pat00133
를 적용하는 바틀렛 방법(도 14의 주황색 그래프 참조)은 단지 2개의 DOA [-8.4˚, 3.2˚]와 [-8.6˚, 2.7˚]만을 각각 추정한다. 따라서 이들 두 방법은 [θ2, θ3] = [1.5˚, 7.5˚]에 놓인 표적을 분간하지 못한다. Further, simulation was performed for the case where three target vehicles exist in the FOV of the radar. This simulation was performed while maintaining the same simulation conditions as in Fig. 10 except for the target information. 14 illustrates an exemplary normalized Bartlet pseudo-spectrum for three adjacent targets located in three directions [θ 1 , θ 2 , θ 3 ] = [-8 °, 1.5 °, 7.5 °]. As shown in FIG. 14, the conventional Bartlet method (see the blue graph in FIG. 14) and the existing transformation matrix
Figure 112017066647630-pat00133
(See the orange graph in FIG. 14) estimates only two DOAs [-8.4 DEG, 3.2 DEG] and [-8.6 DEG, 2.7 DEG], respectively. Therefore, these two methods do not distinguish the target placed at [θ 2 , θ 3 ] = [1.5 °, 7.5 °].

그러나 본 발명에 따른 변환행렬

Figure 112017066647630-pat00134
을 적용하는 경우, 세 가지 다른 DOA를 식별할 수 있다(도 14의 노란색 그래프 참조). 또한, 전력 보정된 보간된 수신신호 벡터
Figure 112017066647630-pat00135
를 이용하는 경우, 추정된 DOA는 가장 정확한 추정값 인 [-8.6˚, 1.8˚, 8.8˚]로 계산된다(도 14의 보라색 그래프 참조). However,
Figure 112017066647630-pat00134
, Three different DOAs can be identified (see the yellow graph in Fig. 14). Further, the power-corrected interpolated received signal vector
Figure 112017066647630-pat00135
, The estimated DOA is calculated to be the most accurate estimation values [-8.6 DEG, 1.8 DEG, 8.8 DEG] (see the purple graph in Fig. 14).

또한 본 발명에 따른 DOA 추정 방법의 성능을 고해상도 DOA 추정 알고리즘으로 알려진 다중 신호 분류 (MUSIC) 알고리즘과 비교한다. MUSIC 알고리즘을 적용하기 위해서는 AIC(Akaike information criterion) 또는 MDL(minimum description length)을 이용하여 표적의 수를 미리 추정해야한다. 표적의 수가 잘 추정된다면

Figure 112017066647630-pat00136
, 가장 정확한 성능을 보여준다.The performance of the DOA estimation method according to the present invention is also compared with a multi-signal classification (MUSIC) algorithm known as a high-resolution DOA estimation algorithm. To apply the MUSIC algorithm, the number of targets must be estimated in advance using Akaike information criterion (AIC) or minimum description length (MDL). If the number of targets is well estimated
Figure 112017066647630-pat00136
, The most accurate performance.

그러나

Figure 112017066647630-pat00137
= 1 또는
Figure 112017066647630-pat00138
= 2와 같이 정확하게 추정되지 않으면 추정 성능이 현저히 저하되어 사용할 수 없다. 또한, MUSIC 알고리즘이 잡음 고유 벡터에 의해 확장된 행렬의 고유값 분해 및 행렬을 수행하기 때문에, 기존의 빔 형성 알고리즘 즉, 바틀렛 알고리즘과 비교하여 추가적인 계산 복잡도 O (N3 + N2 x (2N - 2L - 1))가 발생한다. 게다가 바틀렛 알고리즘은 MUSIC 알고리즘보다 잡음 분산에 더 강한 특성을 갖는다. 따라서 자동차 레이더의 경우, 바틀렛 알고리즘은 시끄러운 도로 환경에서 타겟의 DOA를 안정적으로 추정하는 데 더 적합 할 수 있다.But
Figure 112017066647630-pat00137
= 1 or
Figure 112017066647630-pat00138
= 2, the estimation performance is significantly degraded and can not be used. In addition, since the MUSIC algorithm performs the eigenvalue decomposition and matrix of the matrix extended by the noise eigenvectors, the additional computational complexity O (N 3 + N 2 x (2N - 2L - 1)). In addition, the Bartlett algorithm has a stronger noise variance than the MUSIC algorithm. Thus, in the case of automotive radar, the Bartlett algorithm may be more suitable for stably estimating the DOA of the target in a noisy road environment.

또한, 시뮬레이션은 4 개의 안테나 어레이 요소들뿐만 아니라 3개 또는 5개와 같은 다른 개수의 안테나 어레이 요소들에 대해서도 수행하였다. 어레이 요소들의 수가 3 일 때, 어레이 요소들의 원래 위치는 [d1, d2, d3] = [0, 1.5λ, 3λ]로 주어진다. 이 어레이는 최소 중복 어레이로 변환되며, 보간된 어레이 요소들은 [g1, g2, g3] = [0, 1λ, 3λ]에 위치한다. 또한 표적의 위치는 [θ1, θ2] = [-4˚, 6.5˚]이며, FOV의 범위는 -15˚ ~ 15˚라고 가정한다. 주어진 어레이의 반 전력 빔폭이 12이므로, 어레이는 매우 낮은 각도 분해능을 가지며, 주어진 DOA는 기존의 바틀렛 알고리즘과 구별되지 않았다. In addition, the simulation was performed for four antenna array elements as well as for a different number of antenna array elements, such as three or five. When the number of array elements is 3, the original position of the array elements is given by [d 1 , d 2 , d 3 ] = [0, 1.5λ, 3λ]. This array is converted to the least redundant array, and the interpolated array elements are located at [g 1 , g 2 , g 3 ] = [0, 1λ, 3λ]. The position of the target is assumed to be [θ 1 , θ 2 ] = [-4 °, 6.5 °], and the FOV range is -15 ° to 15 °. Because the half-power beam width of a given array is 12, the array has very low angular resolution, and the given DOA is indistinguishable from the conventional Bartlett algorithm.

또한 5개의 안테나 어레이 요소들에 대해, 원래의 안테나 어레이 요소들의 위치는 [d1, d2, d3, d4, d5] = [0, 2.25λ, 4.5λ, 6.75λ, 9λ]로 주어지고, 위치 [g1, g2, g3, g4, g5] = [0, 1λ, 4λ, 7λ, 9λ]로 변환된다. 이 경우, FOV는 LRR의 그것과 동일하게 주어지고, 표적들은 [θ1, θ2] = [-1˚, 3˚]에 위치된다. 이 DOA들은 주어진 어레이의 반 전력 빔폭이 4.5˚이므로, 기존의 바틀렛 알고리즘을 사용하여서는 구분할 수 없다. 3개와 5개의 안테나 어레이 요소들을 갖는 두 가지 경우, 어레이 SNR을 0 dB에서 10 dB로 증가시켰을 때의 해상도 확률과 RMSE가 도 15와 도 16에 각각 표시되어 있다. 도면에서 알 수 있듯이, 본 발명에 따른 방법은 안테나 어레이 요소의 개수가 3개인 경우나 5개인 경우에도 더 나은 성능을 보여준다.Also for the five antenna array elements, the positions of the original antenna array elements are [d 1 , d 2 , d 3 , d 4 , d 5 ] = [0, 2.25λ, 4.5λ, 6.75λ, 9λ] And are converted into positions [g 1 , g 2 , g 3 , g 4 , g 5 ] = [0, 1λ, 4λ, 7λ, 9λ]. In this case, the FOV is given the same as that of the LRR, and the targets are located at [θ 1 , θ 2 ] = [-1 °, 3 °]. These DOAs can not be distinguished using the conventional Bartlett algorithm because the half-power beam width of a given array is 4.5˚. In both cases with three and five antenna array elements, the resolution probability and RMSE when the array SNR is increased from 0 dB to 10 dB are shown in FIGS. 15 and 16, respectively. As can be seen from the figure, the method according to the present invention shows better performance when the number of antenna array elements is three or five.

위와 같은 시뮬레이션 결과를 검증하기 위해, 자동차 시험장에서 자동차 LRR을 사용하여 실제 측정을 수행해본 결과가 표 2에 정리되어 있다. In order to verify the above simulation results, the actual measurement results using the automobile LRR in the automobile test center are summarized in Table 2.

[-3.5˚, 2.5˚]에 위치한 두 인접 표적에 대한 해상도 확률 P r 및 RMSEThe resolution probabilities P r and RMSE for two adjacent targets located at [-3.5˚, 2.5˚] DOA 추정 DOA estimation P r (%) P r (%) RMSE (˚)RMSE (°) (a) 기존의 바틀렛 방법 수행 시(a) When carrying out the conventional Bartlett method 0.0 0.0 5.725.72 (b)

Figure 112017066647630-pat00139
Figure 112017066647630-pat00140
로 바틀렛 방법 수행 시(b)
Figure 112017066647630-pat00139
Wow
Figure 112017066647630-pat00140
When performing the Bartlett method 0.670.67 6.206.20 (c)
Figure 112017066647630-pat00141
Figure 112017066647630-pat00142
로 바틀렛 방법 수행 시
(c)
Figure 112017066647630-pat00141
Wow
Figure 112017066647630-pat00142
When performing the Bartlett method
46.746.7 4.654.65
(d)
Figure 112017066647630-pat00143
Figure 112017066647630-pat00144
로 바틀렛 방법 수행 시
(d)
Figure 112017066647630-pat00143
Wow
Figure 112017066647630-pat00144
When performing the Bartlett method
68.068.0 3.903.90

시뮬레이션 결과와 유사하게, 본 발명에 따른 변환행렬

Figure 112017066647630-pat00145
를 이용한 DOA 추정 방법(표 2의 (c) 참조)은 기존의 바틀렛 방법이나 기존의 변환행렬
Figure 112017066647630-pat00146
를 이용한 바틀렛 방법(표 2의 (a)와 (b) 참조)에 비해 더 우수한 각도 분해능 및 추정 정확도를 나타낸다. 또한, 수신 전력 조정 신호는 다른 신호보다 높은 분해능 및 추정 성능을 나타낸다(표 2의 (d) 참조). 결과적으로 본 발명이 제안하는 로그 도메인 상에서의 보간 방법을 이용한 DOA 추정 방법은 실제 실험 데이터를 통해서도 성능이 검증되었다.Similar to the simulation results, the transformation matrix according to the present invention
Figure 112017066647630-pat00145
(See Table 2 (c)) can be obtained by a conventional Bartlett method or a conventional transformation matrix
Figure 112017066647630-pat00146
(See (a) and (b) of Table 2). In addition, the received power adjustment signal shows higher resolution and estimation performance than other signals (see Table 2 (d)). As a result, the DOA estimation method using the interpolation method on the log domain proposed by the present invention has been verified through actual experimental data.

위와 같이 시뮬레이션과 측정 결과로부터, 본 발명에 따른 로그 영역에서의 변환행렬을 이용하는 수신신호의 DOA 추정 방법이 LLS의 방법을 이용하여 도출된 종래의 변환행렬을 사용하는 DOA 추정보다 훨씬 우수한 각도 분해능 및 추정 정확도를 나타냄을 확인할 수 있다.From the above simulation and measurement results, it can be seen that the DOA estimation method of the received signal using the transformation matrix in the log domain according to the present invention is superior to the DOA estimation using the conventional transformation matrix derived using the LLS method It can be confirmed that the estimation accuracy is shown.

본 발명은 안테나 어레이 장치에서 수신신호의 DOA를 측정하는 데 이용될 수 있다. 예를 들어 차량용 레이더 장치에 적용하여 전방의 표적 차량을 정확하게 인식하는 데 효과적이다. The present invention can be used to measure DOA of a received signal in an antenna array device. For example, it is effective to accurately recognize a target vehicle ahead by applying it to a radar apparatus for a vehicle.

10: ULA 안테나 기반 레이더 시스템 20: 송신부
24: 파형 발생기(Waveform Generator) 26: 발진기
28: 전력 증폭기(PA) 30: 수신부
32: 저잡음 증폭기(LNA) 34: 혼합기
36: 저역통과필터(LPF) 38: AD변환기
40: 신호 처리부 50: 사용자 인터페이스 (User Interface)
60: ULA 수신 안테나 70: ULA 송신 안테나
10: ULA antenna-based radar system 20:
24: Waveform generator 26: Oscillator
28: Power amplifier (PA) 30: Receiver
32: Low Noise Amplifier (LNA) 34: Mixer
36: Low-pass filter (LPF) 38: AD converter
40: a signal processor 50: a user interface
60: ULA receiving antenna 70: ULA transmitting antenna

Claims (14)

균일 선형 어레이(ULA) 형태로 배열된 안테나 어레이 요소들에 있어서,
원래의 안테나 어레이 요소들의 제1 스티어링 행렬과, 상기 원래의 안테나 어레이 요소들에 기초하여 보간한 보간 안테나 어레이 요소들의 제2 스티어링 행렬 각각에 대해 로그를 취하는 단계;
로그를 취한 상기 제1 스티어링 행렬을 로그를 취한 상기 제2 스티어링 행렬로 변환하는 변환행렬을 최적화 기법을 이용하여 생성하는 단계;
상기 원래의 안테나 어레이 요소들의 수신신호들의 위상에 대하여 상기 변환행렬을 적용하여 상기 보간 안테나 어레이 요소들에 대응되는 위상을 갖는 보간 수신신호들로 변환하는 단계; 및
소정의 도래각(DOA) 추정 알고리즘에 기초하여 상기 보간 안테나 어레이 요소들의 상기 보간 수신신호들의 도래각(DOA)을 추정하는 단계를 포함하는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법.
For antenna array elements arranged in the form of a uniform linear array (ULA)
Taking a log for a first steering matrix of the original antenna array elements and a second steering matrix of interpolated antenna array elements based on the original antenna array elements;
Generating a transformation matrix for transforming the first steering matrix obtained by taking the log into the second steering matrix obtained by using an optimization technique;
Applying the transform matrix to the phase of the received signals of the original antenna array elements to transform them into interpolated received signals having a phase corresponding to the interpolated antenna array elements; And
Estimating the DOA of the interpolated received signals of the interpolated antenna array elements based on a predetermined DOA estimation algorithm. Arrival angle estimation method.
제1항에 있어서, 상기 최적화 기법은 상기 수신신호들 간의 관계에 기초하여 상기 원래의 안테나 어레이 요소들과 상기 보간 안테나 어레이 요소들 간의 위상차를 최소화 하는 선형 최소 제곱법(Linear Least Square: LLS)인 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법.2. The method of claim 1, wherein the optimization technique is a linear least squares (LLS) method that minimizes the phase difference between the original antenna array elements and the interpolation antenna array elements based on the relationship between the received signals. And estimating the arrival angle of the received signal based on the log-area antenna array interpolation. 제1항에 있어서, 상기 변환행렬은 상기 수신신호들의 진폭은 그대로 유지하면서 상기 수신신호들의 위상만을 로그 영역에서 조합하여 상기 보간 수신신호의 위상으로 변환하는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법.2. The method of claim 1, wherein the transform matrix combines only the phases of the received signals in the log domain while maintaining the amplitudes of the received signals intact, and converts the combined signals into a phase of the interpolated received signal. A method for estimating an arrival angle of a received signal based on a received signal. 제1항에 있어서, 상기 안테나 어레이 요소들에 관해 설정된 관측 시야에 맞는 상기 원래의 안테나 어레이 요소들의 상기 제1 스티어링 행렬과, 상기 원래의 안테나 어레이 요소들을 선형 조합하여 원하는 위치로 보간한 상기 보간 안테나 어레이 요소들의 상기 제2 스티어링 행렬을 구하는 단계를 더 포함하는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법.2. The method of claim 1, further comprising: combining the first steering matrix of the original antenna array elements in accordance with an observation field set with respect to the antenna array elements, the first steering matrix of the original antenna array elements, Further comprising: obtaining the second steering matrix of array elements. ≪ Desc / Clms Page number 21 > 제1항에 있어서, 상기 보간 수신신호들의 도래각(DOA) 추정의 정확도를 높이기 위해, 상기 도래각(DOA)을 추정하는 단계 이전에, 상기 변환하는 단계에서 얻어지는 상기 보간 수신신호들 간의 전력 레벨의 편차를 보상하여 상기 보간 수신신호들의 전력 레벨을 실질적으로 균등하게 보정하는 단계를 더 포함하는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법.2. The method of claim 1, further comprising: before the step of estimating the DOA, to increase the accuracy of the DOA estimation of the interpolated received signals, Compensating the deviation of the received signal to substantially equalize the power level of the interpolated received signals. ≪ Desc / Clms Page number 21 > 제5항에 있어서, 상기 보간 수신신호들의 위상은 상기 전력 레벨을 보정하기 전후로 동일하게 유지되는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법.6. The method of claim 5, wherein the phase of the interpolated received signals remains the same before and after correcting the power level. 제1항에 있어서, 상기 로그를 취하는 단계는, 상기 원래의 안테나 어레이 요소들의 관측 시야를 설정하는 단계; 설정된 관측 시야에 맞는 상기 원래의 안테나 어레이 요소들의 상기 제1 스티어링 행렬과 상기 보간 안테나 어레이 요소들의 상기 제2 스티어링 행렬을 생성하는 단계; 상기 제1 및 제2 스티어링 행렬의 각 성분에 대하여 로그를 취하는 단계를 포함하는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 방법.2. The method of claim 1, wherein taking the log further comprises: setting an observation field of view of the original antenna array elements; Generating the first steering matrix of the original antenna array elements and the second steering matrix of the interpolation antenna array elements corresponding to the set observation field; And taking a log for each component of the first and second steering matrices. ≪ Desc / Clms Page number 21 > 제1항 내지 제7항 중 어느 하나의 항에 따른 수신신호의 도래각 추정 방법을 컴퓨터에서 실행하기 위한 컴퓨터 프로그램이 기록된 컴퓨터 가독형(computer-readable) 기록 매체.A computer-readable recording medium on which a computer program for executing a method of estimating an arrival angle of a received signal according to any one of claims 1 to 7 is recorded in a computer. 일렬로 등간격으로 배열되고 전방에서 입사되는 무선신호를 수신하는 복수 개의 원래의 안테나 어레이 요소들을 포함하는 수신 안테나부;
상기 수신 안테나부의 각 안테나를 통해 수신된 복수 개의 수신신호로부터 소정의 신호를 추출하여 복수 개의 디지털 수신신호로 변환하는 수신부; 및
원래의 안테나 어레이 요소들의 제1 스티어링 행렬과, 상기 원래의 안테나 어레이 요소들에 기초하여 보간한 보간 안테나 어레이 요소들의 제2 스티어링 행렬 각각에 대해 로그를 취하는 기능과, 로그를 취한 상기 제1 스티어링 행렬을 로그를 취한 상기 제2 스티어링 행렬로 변환하는 변환행렬을 최적화 기법을 이용하여 생성하는 기능과, 상기 수신부가 제공하는 상기 복수 개의 디지털 수신신호들의 위상에 대하여 상기 변환행렬을 적용하여 상기 보간 안테나 어레이 요소들에 대응되는 위상을 갖는 보간 수신신호들로 변환하는 기능과, 소정의 도래각(DOA) 추정 알고리즘에 기초하여 상기 보간 안테나 어레이 요소들의 상기 보간 수신신호들의 도래각(DOA)를 추정하는 기능을 포함하는 신호처리부를 포함하는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 장치.
A receiving antenna comprising a plurality of original antenna array elements arranged at regular intervals in a row and receiving a radio signal incident from the front;
A receiving unit for extracting a predetermined signal from a plurality of received signals received through each antenna of the receiving antenna unit and converting the predetermined signal into a plurality of digital receiving signals; And
A function of taking a log for each of a first steering matrix of original antenna array elements and a second steering matrix of interpolated antenna array elements interpolated based on said original antenna array elements, To the second steering matrix obtained by taking the log by using an optimization technique; and a function of generating a transformation matrix by using the transformation matrix for the phase of the plurality of digital received signals provided by the receiver, (DOA) estimation algorithm based on a predetermined angle of arrival (DOA) estimation algorithm, a function of estimating an angle of arrival (DOA) of the interpolated received signals of the interpolation antenna array elements And a signal processing unit including a signal processing unit Each estimator arrival of the received signal based on the interpolation.
제9항에 있어서, 상기 최적화 기법은 상기 수신신호들 간의 관계에 기초하여 상기 원래의 안테나 어레이 요소들과 상기 보간 안테나 어레이 요소들 간의 위상차를 최소화 하는 선형 최소 제곱법(Linear Least Square: LLS)인 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 장치.10. The method of claim 9, wherein the optimization technique is a linear least squares (LLS) method that minimizes the phase difference between the original antenna array elements and the interpolation antenna array elements based on the relationship between the received signals. Wherein the received-signal arrival angle estimator is based on log-area antenna array interpolation. 제9항에 있어서, 상기 변환행렬은 상기 수신신호들의 진폭은 그대로 유지하면서 상기 수신신호들의 위상만을 로그 영역에서 조합하여 상기 보간 수신신호의 위상으로 변환하는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 장치.10. The method of claim 9, wherein the transform matrix combines only the phases of the received signals in the log domain while maintaining the amplitudes of the received signals intact, and converts the combined signals into a phase of the interpolated received signal. Based on the received signal. 제9항에 있어서, 상기 신호처리부는 상기 안테나 어레이 요소들에 관해 설정된 관측 시야에 맞는 상기 원래의 안테나 어레이 요소들의 상기 제1 스티어링 행렬과, 상기 원래의 안테나 어레이 요소들을 선형 조합하여 원하는 위치로 보간한 상기 보간 안테나 어레이 요소들의 상기 제2 스티어링 행렬을 구하는 기능을 더 포함하는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 장치.10. The apparatus of claim 9, wherein the signal processing unit is configured to linearly combine the original steering antenna array elements and the first steering matrix of the original antenna array elements matching an observation field of view with respect to the antenna array elements, Further comprising a function of obtaining the second steering matrix of one of the interpolation antenna array elements based on the logarithmic antenna array interpolation. 제9항에 있어서, 상기 신호처리부는, 상기 보간 수신신호들의 도래각(DOA) 추정의 정확도를 높이기 위해, 상기 보간 안테나 어레이 요소들의 상기 보간 수신신호들 간의 전력 레벨의 편차를 보상하여 상기 보간 수신신호들의 전력 레벨을 실질적으로 균등하게 보정하는 기능을 더 포함하는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 장치.10. The apparatus of claim 9, wherein the signal processor compensates for a deviation in the power level between the interpolated received signals of the interpolation antenna array elements to increase the accuracy of the DOA estimation of the interpolated received signals, Further comprising a function of substantially evenly correcting the power level of the received signals based on the log-area antenna array interpolation. 제13항에 있어서, 상기 보간 수신신호들의 위상은 상기 전력 레벨을 보정하기 전후로 동일하게 유지되는 것을 특징으로 하는 로그-영역 안테나 어레이 보간에 기반한 수신신호의 도래각 추정 장치.14. The apparatus of claim 13, wherein the phase of the interpolated received signals remains the same before and after correcting the power level.
KR1020170088392A 2017-07-12 2017-07-12 Method of estimating DOA of received signals based on logarithmic-domain antenna array interpolation, and apparatus for the same KR102001394B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020170088392A KR102001394B1 (en) 2017-07-12 2017-07-12 Method of estimating DOA of received signals based on logarithmic-domain antenna array interpolation, and apparatus for the same

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020170088392A KR102001394B1 (en) 2017-07-12 2017-07-12 Method of estimating DOA of received signals based on logarithmic-domain antenna array interpolation, and apparatus for the same

Publications (2)

Publication Number Publication Date
KR20190007221A KR20190007221A (en) 2019-01-22
KR102001394B1 true KR102001394B1 (en) 2019-07-18

Family

ID=65287894

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020170088392A KR102001394B1 (en) 2017-07-12 2017-07-12 Method of estimating DOA of received signals based on logarithmic-domain antenna array interpolation, and apparatus for the same

Country Status (1)

Country Link
KR (1) KR102001394B1 (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4071507A4 (en) * 2019-12-26 2023-03-08 Huawei Technologies Co., Ltd. Angle of arrival measurement method and device
KR102247740B1 (en) * 2020-01-10 2021-05-04 홍익대학교 산학협력단 Apparatus for estimating direction-of-arrival based on artificial neural network and method thereof
KR102491770B1 (en) * 2021-06-03 2023-01-26 한국항공대학교산학협력단 System for estimating direction of arrival based on esprit and me-esprit and method thereof
CN113589224A (en) * 2021-08-03 2021-11-02 宜宾电子科技大学研究院 DOA estimation method based on enhanced nested array
CN114217266B (en) * 2021-11-30 2024-06-11 清华大学 Processing method for reducing channel error based on antenna data time domain fitting

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004061468A (en) * 2002-07-31 2004-02-26 Fujitsu Ltd Method for estimating arrival direction of multiplex wave using spatial characteristic and receiving beam formation device using the same
JP2004061467A (en) * 2002-07-31 2004-02-26 Fujitsu Ltd Method for estimating arrival direction of period steady signal in multi-path propagation environment and receiving beam formation device using the same
JP2006258529A (en) * 2005-03-16 2006-09-28 Fujitsu Ten Ltd Device and method for estimating radio wave incoming direction
JP2011053056A (en) * 2009-09-01 2011-03-17 Fujitsu Ltd Apparatus for estimating direction of arrival and method thereof
KR101274554B1 (en) * 2012-03-02 2013-06-13 국방과학연구소 Method for estimating direction of arrival and array antenna system using the same
KR101375425B1 (en) * 2013-02-28 2014-03-17 충북대학교 산학협력단 Method for estimating direction of arrival in antenna array

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2988463B2 (en) * 1998-03-24 1999-12-13 日本電気株式会社 Direction finding device and measurement result processing device therefor
KR101137793B1 (en) * 2010-08-09 2012-04-18 국방과학연구소 System, apparatus and method for detecting direction
KR101414574B1 (en) * 2012-11-05 2014-07-03 한남대학교 산학협력단 A Spatially Smoothed Auto-convolution based DOA Estimating Method for Correlated Signals in the Presence of Correlated Antenna Noises
KR101524793B1 (en) * 2013-02-21 2015-06-03 부산대학교 산학협력단 Apparatus and method for estimating direction of arrival using array antenna

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004061468A (en) * 2002-07-31 2004-02-26 Fujitsu Ltd Method for estimating arrival direction of multiplex wave using spatial characteristic and receiving beam formation device using the same
JP2004061467A (en) * 2002-07-31 2004-02-26 Fujitsu Ltd Method for estimating arrival direction of period steady signal in multi-path propagation environment and receiving beam formation device using the same
JP2006258529A (en) * 2005-03-16 2006-09-28 Fujitsu Ten Ltd Device and method for estimating radio wave incoming direction
JP2011053056A (en) * 2009-09-01 2011-03-17 Fujitsu Ltd Apparatus for estimating direction of arrival and method thereof
KR101274554B1 (en) * 2012-03-02 2013-06-13 국방과학연구소 Method for estimating direction of arrival and array antenna system using the same
KR101375425B1 (en) * 2013-02-28 2014-03-17 충북대학교 산학협력단 Method for estimating direction of arrival in antenna array

Also Published As

Publication number Publication date
KR20190007221A (en) 2019-01-22

Similar Documents

Publication Publication Date Title
KR102001394B1 (en) Method of estimating DOA of received signals based on logarithmic-domain antenna array interpolation, and apparatus for the same
AU2022259835B2 (en) Direction of arrival estimation
KR102449214B1 (en) Method for estimating direction of arrival and apparatus for estimating direction of arrival using the same
KR102323655B1 (en) Method and apparatus for estimating direction of arrival using generation of virtual received signals
JP5911331B2 (en) Monopulse angle measuring device and monopulse angle measuring method
RU2382379C2 (en) Method for multi-signal location of radio-frequency radiation sources at one frequency
CN105137409B (en) The sane space-time adaptive processing method of echo signal mutually constrained based on width
CN111948598B (en) Method and device for detecting space domain interference signal
CN111983556B (en) Device and method for estimating angle of arrival
CN110426670B (en) Super-resolution DOA estimation method for external radiation source radar based on TLS-CS
CN110824414A (en) Device and method for estimating angle of arrival
CN103926555A (en) Method for testing amplitude and phase response of antenna array receiver through non-circular signals
JP3918573B2 (en) Radar equipment
Park et al. The fast correlative interferometer direction finder using I/Q demodulator
CN113985345A (en) Ultra-wideband correction amplitude-phase fusion direction finding method and system
CN111458677A (en) Double-channel single-pulse amplitude-comparison angle measurement method and device
Lee et al. Two‐stage DOA estimation method for low SNR signals in automotive radars
Kikuchi et al. Autocalibration algorithm for robust Capon beamforming
CN103248412B (en) A kind of method of satellite navigation Multibeam synthesis
Biswas et al. New high resolution direction of arrival estimation using Compressive Sensing
CN107621623B (en) Signal direction detection method and beam former using same
JP4294634B2 (en) Calibration table creation method and arrival direction estimation method
JP4119719B2 (en) Mobile station direction estimation method and apparatus
JP4660562B2 (en) Mobile station direction estimation method and apparatus
Ping et al. Study of 2D DOA estimation for uniform circular array in wireless location system

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right