KR100941697B1 - Apparatus and method for estimating location of single source in uniform circular array - Google Patents
Apparatus and method for estimating location of single source in uniform circular array Download PDFInfo
- Publication number
- KR100941697B1 KR100941697B1 KR1020090062620A KR20090062620A KR100941697B1 KR 100941697 B1 KR100941697 B1 KR 100941697B1 KR 1020090062620 A KR1020090062620 A KR 1020090062620A KR 20090062620 A KR20090062620 A KR 20090062620A KR 100941697 B1 KR100941697 B1 KR 100941697B1
- Authority
- KR
- South Korea
- Prior art keywords
- sound source
- source signal
- incident
- distance
- equation
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
- G01S5/20—Position of source determined by a plurality of spaced direction-finders
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S1/00—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
- G01S1/72—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith using ultrasonic, sonic or infrasonic waves
- G01S1/74—Details
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Description
본 발명은 위치 추정 장치 및 방법에 관한 것으로, 특히 균일 원 배열(Uniform Circular Array) 센서의 중앙 대칭성(Centrosymmetry)과 음원의 비순환성(Noncircularity)을 이용하여 음원 신호의 위치를 추정하는 균일 원 배열 센서에서 단일 음원 신호의 위치 추정 장치 및 방법에 관한 것이다.BACKGROUND OF THE
일반적으로 이동통신 서비스의 욕구가 점차 다양해지면서 그 수요도 크게 늘어남에 따라 기존의 기술만으로는 증대되는 고성능, 고용량의 필요를 충족시키기는 점차 어려워지고 있다. 그 외에도 여러 패킷 및 영상 신호 전송을 위한 고성능 데이터 및 영상 서비스 시스템에 대한 필요도 크게 부각되고 있는 것이 현 실정이어서, 향후 이동통신시스템은 기존의 셀룰러 및 개인휴대통신과 비교할 때 고품질에 훨씬 높은 용량을 필요로 하는 멀티미디어 통신서비스가 될 것이며, 음질조차도 유선 통화음질 정도나 그 이상의 고품질 음성 서비스를 요구받게 된다. 이러한 간섭신호 및 채널 특성에 의한 성능 열화 현상에 대한 해결책으로 상용화 개발가치가 가장 높은 유망 핵심기술로 평가되어지고 있는 것이 스마트 안테나 기술이다. 기존의 두 개 디버시티 안테나를 사용하여 다중경로 신호를 결합하는 경우와는 달리 스마트 안테나 기술은 안테나 배열과 첨단 고성능 디지털 신호처리 기술을 이용하여 RF 신호환경의 변화에 따른 적응적 안테나 빔 패턴 제어에 의해 송, 수신 성능 및 용량의 극대화를 가능케 하는 첨단 신호처리 및 안테나 기술이다. 즉, 전 방향으로 방사 빔을 형성하는 대신 해당 가입자에게만 지향성의 빔을 방사함으로써 섹터에서 활동하고 있는 전 가입자에게 신호간섭 효과를 최소화함으로써 통신품질과 시스템 채널용량을 그 만큼 높일 수 있는 기술이다In general, as the needs of mobile communication services are diversified, the demands are greatly increased, making it difficult to meet the needs of high performance and high capacity, which are only increased by existing technologies. In addition, the need for high-performance data and video service systems for the transmission of multiple packets and video signals has been highlighted. As a result, future mobile communication systems will have higher quality and much higher capacities than conventional cellular and personal cellular communications. It will be the multimedia communication service that is needed, and even the sound quality will be required for the high quality voice service of the wire call quality or more. As a solution to the performance degradation caused by the interference signal and the channel characteristics, the smart antenna technology has been evaluated as a promising core technology with the highest commercial development value. Unlike the case of combining multipath signals using two conventional diversity antennas, smart antenna technology uses an antenna array and advanced high-performance digital signal processing technology to control the adaptive antenna beam pattern according to the change of RF signal environment. Advanced signal processing and antenna technology enables maximum transmission and reception performance and capacity. In other words, instead of forming a radiation beam in all directions, by radiating a directional beam only to the subscriber, it is possible to increase the communication quality and system channel capacity by minimizing the signal interference effect to all subscribers operating in the sector.
이러한 스마트 안테나는 안테나 배열 크기에 따라 기지국 및/또는 단말기에 위치할 수 있다. 스마트 안테나가 기지국에 위치할 경우, 상향 링크에서 적응적으로 원하는 방향으로부터 수신하고 하향 링크에서는 적응적으로 원하는 방향으로 송신하도록 하여 원하는 사용자에게 안테나 이득과 다양성 이득을 증가시키는 동시에 상향 링크에서 다른 방향으로부터 수신되는 간섭 신호나 하향 링크에서 다른 방향으로 송신되는 간섭 신호를 줄일 수 있다. 이렇게 스마트 안테나를 사용하면, 더 많은 사용자를 수용할 수 있으며 그에 따라 시스템 용량을 늘릴 수 있다. 이러한 스마트 안테나의 구성 중에는 안테나 배열(Antenna Array)이 있다.Such a smart antenna may be located in the base station and / or the terminal according to the antenna array size. When the smart antenna is located in the base station, it is adaptively received in the uplink on the uplink and adaptively transmitted in the downlink on the downlink, thereby increasing antenna gain and diversity gain to the desired user, while simultaneously increasing the antenna gain and diversity gain on the uplink. The interference signal received or the interference signal transmitted in the other direction in the downlink can be reduced. This smart antenna can accommodate more users and increase system capacity accordingly. Among the smart antennas, there is an antenna array.
안테나 배열은 다수의 안테나로 구성되며, 원하는 안테나 빔 패턴을 생성하는데, 안테나 수가 증가할수록 빔 패턴이 좁아져서 성능이 증가하지만 시스템의 복잡도를 고려하여 보통 4~12개가 사용된다. 안테나 빔 패턴에 따라 균일 선 배열(Uniform Linear Array) 센서 및 균일 원 배열(Uniform Circular Array) 센서 등 이 있다. The antenna array is composed of a plurality of antennas, and generates the desired antenna beam pattern. As the number of antennas increases, the beam pattern is narrowed and the performance increases, but 4 to 12 are usually used in consideration of the system complexity. According to the antenna beam pattern, there are a uniform linear array sensor and a uniform circular array sensor.
종래 균일 원 배열 센서에 신호가 입사되는 경우, 음원의 위치를 추정하는 과정을 설명하면 다음과 같다.When a signal is incident on the conventional uniform circle array sensor, a process of estimating the position of the sound source will be described.
M+1개의 센서로 이루어진 반지름 R인 균일 원 배열센서는 원주상에 M개의 센서가 일정한 간격으로 형성되어 있으며, 원 중앙에도 1개의 센서가 구비되어 있다. 이러한 균일 원 배열 센서에 음원으로부터 신호가 입사되면, m번째 센서에서의 n번째 시간에 수신되는 신호는 아래 <수학식 1>로 표현된다.In the uniform circle array sensor having a radius R composed of M + 1 sensors, M sensors are formed at regular intervals on the circumference, and one sensor is provided at the center of the circle. When a signal is input from the sound source to the uniform circle array sensor, the signal received at the nth time in the mth sensor is expressed by
상기 <수학식 1>에서 s(n)은 BPSK(Binary Phase Shift Keying), ASK(Amplitude Shift Keying), AM(Amplitude Modulation) 변조, PAM(Pulse Amplitude Modulation) 변조 신호등과 같이 비순환성 신호이다. 이러한 s(n)은 일반적으로 상관관계(Correlation)와 공액 상관관계(Conjugate Correlation)가 각각 과 으로 표현된다. 또한 입사하는 음원신호의 초기 위상 은 미지(unknown)의 임의의 값으로 가정한다. 는 특성을 가지는 가우시안(Gaussian) 백색잡음으로 가정한다. 여기서, 는 방위각을 나타내며, 는 고각을 나타내며, r은 거리를 나타내며, 는 기대치 연산자(Expectation operator)을 나타내며, 는 음원 신호의 파워를 나타내며, 는 잡음신호의 파워를 나타내며, 는 크로네커 델타 함수(Kronecker delta function)를 나타내며, 는 입사되는 신호의 파장을나타내며, R은 균일 원 배열 센서의 반지름을 나타내며, M은 센서의 개수를 나타내며, m은 m번째 센서를 나타내며, 은 음원 신호의 2차원 입사각을 나타낸다.In
그리고, 상기 <수학식 1>을 벡터 형태로 표현하면 아래 <수학식 2>로 표현된다.When
상기 <수학식 2>에서 T는 벡터의 전치연산(Transpose)이다. 그리고, x(n)은 공분산 행렬(Covariance Matrix)과 고유치분해(Eigendecomposition)는 아래 <수학식 3>과 같이 표현된다.In
상기 <수학식 3>에서 는 R의 고유치(Eigenvalue)이며, 는 R의 고유벡터(Eigenvector)이며, 는 신호 부공간(Signal Subspace)이며, 는 잡음 부공간(Noise Subspace)이며, H는 행렬의 공액전치연산(Conjugate Transpose)이다.3차원 MUSIC(Multiple Signal Classification)은 신호부공간과 잡음 부공간의 직교성(Othorgonality)을 이용하여, 아래 <수학식 4>와 같은 음원 스펙트럼을 방위각 θ, 고각 φ, 거리 r의 3개의 변수를 동시에 최대화 시키는 의 탐색을 통해 음원의 위치를 추정한다.In <
상기 <수학식 4>에서 의 m번째 항은 아래 <수학식 5>와 같이 표현된다.In Equation 4 above The m-th term of is expressed as
상기 <수학식 5>에서 이다.In
그런데, 종래 3차원 MUSIC등과 같이 탐색 기반의 음원 위치 추정 기법은 최적의 성능을 제공하지만, 탐색격자 설정기준의 모호하고, 입사된 음원의 방위각, 고각, 거리를 함께 구하기 때문에, 3차원 탐색에 따른 고연산량이 요구되는 문제점이 있다.However, the search-based sound source position estimation method, like the conventional three-dimensional MUSIC, provides the best performance, but since the search lattice criteria are ambiguous, and the azimuth, elevation, and distance of the incident sound source are obtained together, There is a problem that a high amount of computation is required.
본 발명은 상술한 종래 문제점을 해결하기 위해서 착안된 것으로서, 상기한 바를 달성하기 위한 본 발명은 균일 원 배열 장치의 중앙 대칭성과 음원의 비순환성을 결합하여, 3차원 탐색이 필요 없는 저연산량의 간단한 처리 구조를 가지는 균일 원 배열 센서에서 음원 위치 추정 장치 및 방법을 제공한다.The present invention has been conceived to solve the above-mentioned conventional problems, and the present invention for achieving the above-described simple combination of the central symmetry of the uniform circle arrangement device and the acyclic nature of the sound source, so that a simple calculation of a low calculation amount without the need for three-dimensional search Provided are a sound source position estimation apparatus and method in a uniform circle array sensor having a processing structure.
상술한 바를 달성하기 위한 본 발명은 균일 원 배열 센서에서 음원 신호의 위치 추정 방법에 있어서, 음원 신호가 입사되면, 상기 원 배열에서 대칭관계에 있는 센서간의 상관 관계를 추정하여 상기 입사된 음원 신호의 입사각을 계산하는 제 1 과정과, 상기 대칭관계에 있는 센서간의 공액 상관관계를 추정하여 상기 입사된 음원 신호의 입사각 및 거리를 계산하는 제2 과정과, 상기 제1 과정에서 계산된 입사각 및 제2 과정에서 계산된 입사각 및 거리를 이용하여, 상기 입사된 음원 신호의 거리를 계산하는 과정과, 상기 계산된 거리를 이용하여 상기 입사된 음원 신호의 위치를 추정하는 과정을 포함한다.According to an aspect of the present invention, there is provided a method of estimating a position of a sound source signal in a uniform circle array sensor. When a sound source signal is incident, the correlation between the sensors having a symmetric relationship in the circle array is estimated. A first process of calculating an angle of incidence, a second process of estimating the angle of incidence and distance of the incident sound source signal by estimating a conjugate relationship between the symmetrical sensors, and the angle of incidence and the second computed in the first process Comprising a step of calculating the distance of the incident sound source signal using the incident angle and the distance calculated in the process, and estimating the position of the incident sound source signal using the calculated distance.
또한, 상술한 바를 달성하기 위한 본 발명은 원 배열(Uniform Circular Array) 센서에서 음원 신호의 위치 추정 장치에 있어서, 음원 신호가 입사되면, 상기 균일 원 배열에서 대칭관계에 있는 센서간의 상관 관계를 추정하여 상기 입사된 음원 신호의 제1 입사각을 계산하고, 상기 대칭관계에 있는 센서간의 공액 상관관계를 추정하여 상기 입사된 음원 신호의 제2 입사각 및 거리를 계산하고, 상기 계산된 제 1 및 제 2입사각과 상기 거리를 이용하여 상기 입사된 음원 신호의 거리를 계산하는 것을 특징으로 하는 음원 신호의 위치 추정 장치를 포함한다. In addition, the present invention for achieving the above-mentioned in the apparatus for estimating the position of the sound source signal in a uniform circular array sensor, when the sound source signal is incident, the correlation between the sensors in the symmetrical relationship in the uniform circle array is estimated The first incident angle of the incident sound source signal is calculated, the conjugate angle between the sensors in the symmetric relationship is estimated, the second incident angle and distance of the incident sound source signal are calculated, and the calculated first and second And an apparatus for estimating the distance of the incident sound source signal using the incident angle and the distance.
상술한 바와 같이, 본 발명에 따른 균일 원 배열 장치는 입사된 음원 신호의 2차원 입사각(방위각과 고각)을 먼저 구하고, 입사된 음원 신호의 거리를 분리하여 계산함으로써 3차원 MUSIC에 비해 구조가 간단하며, 적은 연산량으로 음원의 위치 추정할 수 있는 효과가 있다.As described above, the uniform circle array device according to the present invention obtains a two-dimensional incidence angle (azimuth and elevation) of the incident sound source signal first, and calculates by separating the distance of the incident sound source signal, which is simpler than the three-dimensional MUSIC. In addition, there is an effect that can estimate the position of the sound source with a small amount of calculation.
이하, 첨부된 도면을 참조하여 본 발명의 바람직한 실시예에 대한 동작 원리 를 상세히 설명한다. 하기에서 본 발명을 설명함에 있어 관련된 공지 기능 또는 구성에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우에는 그 상세한 설명을 생략할 것이다. 그리고 후술되는 용어들은 본 발명에서의 기능을 고려하여 정의된 용어들로서 이는 사용자, 사용자의 의도 또는 관례 등에 따라 달라질 수 있다. 그러므로 그 정의는 본 명세서 전반에 걸친 내용을 토대로 내려져야 할 것이다.Hereinafter, with reference to the accompanying drawings will be described in detail the operating principle of the preferred embodiment of the present invention. In the following description of the present invention, detailed descriptions of well-known functions or configurations will be omitted if it is determined that the detailed description of the present invention may unnecessarily obscure the subject matter of the present invention. In addition, terms to be described below are terms defined in consideration of functions in the present invention, and may be changed according to a user, a user's intention or custom. Therefore, the definition should be made based on the contents throughout the specification.
도 1은 본 발명의 실시예에 따른 균일 원 배열과 음원의 기하학적 구조를 나타낸 예시도이다.1 is an exemplary view showing a uniform circle arrangement and the geometry of the sound source according to an embodiment of the present invention.
도시된 바와 같이, M+1개의 센서로 이루어진 반지름 R인 균일 원 배열센서는 x, y, z축의 3차원 기하학적인 구조이며, 반지름이 R인 원주상에 M개의 센서가 일정한 간격으로 형성되어 있으며, 원 중앙에도 1개의 센서가 구비되어 있다. 이러한, M+1개의 센서 중 각각은 원점을 기준으로 서로 대칭적으로 형성되어 있다. 이러한, 반지름이 R인 균일 원 배열 센서에 음원으로부터 신호가 입사되면, 각각의 센서들의 상관 관계와 공액 상관 관계는 아래 <수학식 6>과 <수학식 7>을 통해 구한다.As shown, the uniform circle array sensor of radius R consisting of M + 1 sensor is a three-dimensional geometric structure of the x, y, z axis, M sensors are formed at regular intervals on the circumference of radius R, One sensor is also provided in the center of the circle. Each of the M + 1 sensors is formed symmetrically with respect to the origin. When a signal is input from the sound source to the uniform circular array sensor having a radius of R, the correlations and conjugate correlations of the respective sensors are obtained through Equations 6 and 7 below.
상기 <수학식 6>에서 *는 공액 연산자(Conjugate)이다.In Equation 6, * is a conjugate operator (Conjugate).
상기 상관 관계라 함은 대칭관계에 있는 센서에 입사되는 음원 신호의 닮은 비를 의미하는 것이다. 예를 들어, M=12인 경우 도 1에서 2번째 센서의 대칭 센서는 8번째 센서이며, 이들 2번째와 8번째 대칭관계라 함은 2번째 센서에서 수신되는 신호와 8번째 센서에서 수신되는 신호가 얼마나 닮아 있는지를 나타내는 닮은 비를 말한다. 그리고, 상기 <수학식 6> 및 <수학식 7>을 이용하여 입사된 신호의 입사각과 음원의 거리를 구할 수 있으며, 이하, <수학식 6>을 이용하여 2차원 입사각을 추정하는 과정을 설명하면 다음과 같다.The correlation means a similar ratio of sound source signals incident on a symmetrical sensor. For example, in the case of M = 12, the symmetric sensor of the second sensor in FIG. 1 is the eighth sensor, and the second and eighth symmetry relationships are the signals received from the second sensor and the signals received from the eighth sensor. Resembles a rain that represents how much resemblance is. Then, the angle of incidence of the incident signal and the distance of the sound source can be obtained using Equations 6 and 7, and a process of estimating a two-dimensional angle of incidence using Equation 6 will now be described. Is as follows.
균일 원배열센서의 중앙 대칭성 과 음원의 비순환성 을 이용하면, 상기 <수학식 6>과 <수학식 7>은 아래 <수학식 8>과 <수학식 9>로 정리된다.Central Symmetry of Uniform Circular Array Sensors Cyclicity of Music and Sound Using Equation 6 and Equation 7 are summarized as Equation 8 and Equation 9 below.
상기 <수학식 8>에서 의 위상은 음원의 2차원 입사각 의 함수이고, 의 위상은 2차원 입사각 과 거리 r의 함수인 점을 이용하여, 음원의 3차원 위치를 2차원 입사각과 거리로 분리하여 각각을 순차적으로 추정한다. 그리고, 상기 <수학식 8> 및 <수학식 9>에서 는 음원 신호의 파워를 나타내며, 는 입사되는 신호의 파장을 나타내며, R은 균일 원 배열 센서의 반지름을 나타내며, 은 음원 신호의 2차원 입사각을 나타내며, M은 센서의 개수를 나타내며, m은 m번째 센서를 나타낸다.In Equation 8 above The phase of the two-dimensional incident angle of the sound source Is a function of The phase of the two-dimensional incident angle Using the points as a function of and the distance r, the three-dimensional position of the sound source is separated into two-dimensional incidence angle and distance, and each is sequentially estimated. And, in Equation 8 and Equation 9 Indicates the power of the sound source signal, Is the wavelength of the incident signal, R is the radius of the uniformly circular array sensor, Denotes a two-dimensional incident angle of the sound source signal, M denotes the number of sensors, and m denotes the m-th sensor.
이하, 본 발명의 실시예에 따른 2차원 입사각을 추정하는 기법과 음원의 거리를 추정하는 기법을 설명하면 다음과 같다.Hereinafter, a technique for estimating a two-dimensional incident angle and a technique for estimating a distance of a sound source according to an embodiment of the present invention will be described.
1. 본 발명의 실시예에 따른 2차원 입사각 추정 기법1. Two-dimensional incident angle estimation technique according to an embodiment of the present invention
<수학식 8>에서 의 위상을 아래 <수학식 10>으로 정의한다.In Equation 8 The phase of is defined by
이때, 의 위상의 모호성을 제거하기 위해, 즉 을 만 족시키기 위해 균일 원 배열 센서의 반지름은 으로 가정한다. 상기 <수학식 10>에서 은 N개의 신호 샘플로부터 의 계산을 통해 의 추정치 를 구할 수 있으며, 이를 벡터 형태로 표현하면 아래 <수학식 11>과 같다.At this time, To eliminate the ambiguity of the phase of, i.e. In order to satisfy the Assume In
상기 <수학식 11>에서 이고, 이고, 이다.In Equation 11 above ego, ego, to be.
상기 <수학식 11>을 이용하여 최소자승법을 이용하여 b의 해를 구하면 아래 <수학식 12>로 표현된다.The solution of b using the least square method using Equation 11 is expressed by Equation 12 below.
상기 <수학식 12>를 통해서 과 의 관계를 이용하여 음원의 2차원 입사각을 추정하면, 아래 <수학식 13>과 같다.Through Equation 12 above and Estimate the two-dimensional angle of incidence of the sound source using the relationship of Equation 13 below.
상기 <수학식 14>에서 의 4분면은 과 의 부호로 결정한다.In Equation 14 above Quadrant of and Determined by the sign of.
상술한 바와 같이, 본 발명에 따른 균일 원 배열 센서의 2차원 입사각은 상술한 <수학식 6> 및 <수학식 10> ~ <수학식 13>을 통해 구할 수 있다.As described above, the two-dimensional incidence angle of the uniform circular array sensor according to the present invention can be obtained through the above-described Equation 6 and
2. 본 발명의 실시예에 따른 거리 추정 기법2. Distance Estimation Technique According to an Embodiment of the Present Invention
본 발명의 실시예에 따른 균일 원 배열 센서에서 입사된 음원의 거리는 상기 <수학식 7>을 이용한다. The distance of the sound source incident from the uniform circle array sensor according to the embodiment of the present invention uses Equation (7).
음원의 거리 추정은 2차원 입사각 추정에서와 같이 최소 자승법을 통한 접근법을 이용할 수 있지만, 상기 <수학식 7>에서 의 위상 에 포함되어 있는 음원 신호의 임의의 초기 위상 에 의해 거리 추정에 필요한 위상 계산 시, 일관되지 않는 오차가 발생되면, 정확한 음원의 거리 추정을 위해서는 먼저 음원 신원의 초기 위상 을 추정하여 이를 보상처리 하는 전처리 과정을 거친 후, 음원의 거리추정을 수행한다.The distance estimation of the sound source may use the least square approach as in the 2D incidence angle estimation, but in Equation 7 Phase of Any initial phase of the sound source signal contained in the Required for distance estimation by If inconsistencies occur in the phase calculation, the initial phase of the sound source's identity must first be used for accurate distance estimation of the sound source. After the preprocessing process to estimate and compensate for this, distance estimation of the sound source is performed.
즉, 음원의 초기 위상의 추정은 균일 원 배열 센서 중심에 위치한 기준 센서 을 이용하여 다음과 같이 <수학식 14>를 통해 수행된다.That is, the estimation of the initial phase of the sound source is a reference sensor located at the center of the uniform circle array sensor. Using Equation 14 is performed as follows.
따라서, 음원 신호의 초기 위상이 보상된 공액 상관관계 와 이의 위상 은 다음과 같이 <수학식 15> 및 <수학식 16>를 통해서 추정할 수 있다.Thus, the conjugated correlation where the initial phase of the sound source signal is compensated And its phase Can be estimated through
상기 <수학식 16>을 벡터 형태로 표현하면 아래 <수학식 17>과 같이 표현할 수 있다.If Equation 16 is expressed in a vector form, it may be expressed as Equation 17 below.
상기 <수학식 17>에서 이고, 이고, 이다.In Equation 17 above ego, ego, to be.
상기 <수학식 17>을 이용하여 최소자승법의 해를 구하면 아래 <수학식 18>과 같다.Using the equation (17), the solution of the least squares method is shown in Equation (18) below.
최종적으로 입사하는 음원의 거리는 아래 <수학식 19>를 통해 구할 수 있다.The distance of the last incident sound source can be obtained through Equation 19 below.
도 2는 본 발명의 실시예에 따른 환형배열센서에서 단일 음원의 3차원 위치 추정 방법을 나타낸 순서도이다.2 is a flowchart illustrating a three-dimensional position estimation method of a single sound source in an annular array sensor according to an embodiment of the present invention.
도시된 바와 같이, 본 발명의 실시 예에 따른 균일 원 배열 센서에서 단일 음원의 3차원 위치 추정 방법은 균일 원 배열 센서에 음원 신호가 입사되면, 대칭관계에 있는 센서간의 상관관계를 추정한다(S212). 즉, 반지름이 R인 균일 원 배열 센서의 둘레에 구비된 임의의 센서와 대칭관계에 있는 센서의 상관관계를 추정한다. 상기 추정된 상관관계를 통해서, 입사된 음원 신호의 2차원 입사각을 파악할 수 있다. 그리고, 대칭 관계에 있는 센서간의 공액 상관관계를 추정한다(S214). 상기 과정(S212)에서 상관 관계가 추정된 두 센서의 공액 상관관계를 추정함으로써, 상기 입사된 음원 신호의 2차원 입사각과 거리를 추정한다. 상기 추정된 상관관계 값에서 상기 입사된 음원 신호의 2차원 입사각을 추정한다(S216). 그리고, 추정된 공액 상관관계 값에서 입사된 음원 신호의 2차원 입사각 및 거리를 추정한다(S218). 상기 218과정(S218)에서 추정된 값들을 연산하여 입사된 음원 신호의 거리를 계산한다(S220). 그리고, 상기 216과정(S216)에서 추정된 입사각과 상기 220과정(S220)에서 계산된 거리를 이용하여 입사된 음원 신호의 위치를 추정한다.As shown, in the method of estimating a three-dimensional position of a single sound source in the uniform circle array sensor according to an exemplary embodiment of the present invention, when a sound source signal is incident on the uniform circle array sensor, the correlation between the sensors having a symmetric relationship is estimated (S212). ). That is, the correlation of the sensor which is symmetrical with the arbitrary sensor provided in the periphery of the uniform circular array sensor whose radius is R is estimated. Through the estimated correlation, the 2D incident angle of the incident sound source signal may be detected. The conjugate correlation between the sensors in the symmetrical relationship is estimated (S214). By estimating the conjugate correlation of two sensors whose correlation is estimated in step S212, the two-dimensional incidence angle and distance of the incident sound source signal are estimated. A two-dimensional incident angle of the incident sound source signal is estimated from the estimated correlation value (S216). Then, the 2D incidence angle and distance of the sound source signal incident on the estimated conjugate correlation value are estimated (S218). The distances of the incident sound source signals are calculated by calculating the values estimated in step 218 (S220). Then, the position of the incident sound source signal is estimated using the incident angle estimated in step 216 (S216) and the distance calculated in step 220 (S220).
그런데, 상기 과정(S218)에서 입사된 음원 신호의 2차원 입사각 및 거리 외에 추정된 값에는 상기 입사된 음원 신호의 초기 위상값이 포함될 수 있다. 입사된 음원 신호의 초기 위상값이 포함되면, 균일 원 배열 센서의 중심에 위치한 기준 센서 를 이용하여 음원 신호의 초기 위상이 보상된 공액 상관관계와 이의 위 상을 추정한다.However, the estimated value other than the two-dimensional incident angle and distance of the incident sound source signal in step S218 may include the initial phase value of the incident sound source signal. If the initial phase value of the incident sound source signal is included, the reference sensor located at the center of the uniform circle array sensor We estimate the conjugated correlation and its phase by which the initial phase of the sound source signal is compensated using.
도 3a, 3b, 3c 및 3d는 본 발명의 실시예에 따른 균일 원 배열 센서에서 단일 음원의 3차원 위치 추정 방법과 3차원 MUSIC의 신호대 잡음비에 대한 RMSE(Root Mean Squared Error)를 비교한 예시도이다.3A, 3B, 3C, and 3D are exemplary diagrams comparing a root mean squared error (RMSE) with respect to a signal-to-noise ratio of a three-dimensional MUSIC and a three-dimensional position estimation method of a single sound source in a uniform circle array sensor according to an exemplary embodiment of the present invention. to be.
도시된 바와 같이, 본 발명의 실시예에 따른 균일 원 배열 센서에서 단일 음워의 3차원 위치 추정 방법의 성능을 알아보기 위해 벤치마크인 3차원 MUSIC의 음원 위치 추정 결과를 비교하였다. 균일 원 배열 센서는 M=12이고, 반지름이 이웃한 두 센서간의 간격이 음원 신호의 반파장에 해당하는 로 설정하였다. 음원은 에 위치하며, 음원 신호는 임의의 초기 위상이 인 BPSK 신호로 설정하였다. N=200인 신호 샘플을 이용하여 500번의 Monte Carlo 시뮬레이션을 수행하였다. 이때 탐색이 필요한 3차원 MUSIC의 경우 격자 크기를 로 설정하였다. 도 3a에서 도 3c의 본 발명에 따른 음원의 3차원 위치 추정 방법은 3차원 MUSIC에 비해 고연산량의 격자 탐색을 수행하지 않고도 넓은 신호대 잡음비 구간에서 3차원 MUSIC과 유사한 성능을 확인 할 수 있다. 그리고, 도 3d는 거리 추정에 필요한 음원 신호의 초기 위상 추정이 성공적으로 이루어지고 있음을 알 수 있다. 초기 위상 추정치의 RMSE는 로 정의하였다. 이때, 와 는 각각 k번째 Monte Carlo 시뮬레이션에서의 의 초기위상과 위상 추정치이다.As shown, in order to determine the performance of the three-dimensional position estimation method of a single tone in the uniform circle array sensor according to an embodiment of the present invention, the sound source position estimation results of the three-dimensional MUSIC benchmark was compared. The uniform circle array sensor has M = 12, and the distance between two adjacent neighboring sensors corresponds to the half wavelength of the sound source signal. Set to. Sound source , The sound source signal has any initial phase Set to the BPSK signal. 500 Monte Carlo simulations were performed using a signal sample with N = 200. In this case, for 3D MUSIC that needs to be searched, Set to. In FIG. 3A, the 3D position estimation method of the sound source according to the present invention of FIG. 3C can confirm performance similar to that of the 3D MUSIC in a wide signal-to-noise ratio section without performing a higher computational lattice search than the 3D MUSIC. 3D shows that the initial phase estimation of the sound source signal necessary for the distance estimation has been successfully performed. The RMSE of the initial phase estimate is Defined as. At this time, Wow Are each kth Monte Carlo simulations Is the initial phase and phase estimate of.
따라서, 본 발명에 따른 음원의 3차원 위치 추정 방법은 상관관계와 공액상관관계 수행에 복소곱셈 (M+1)N과 복소덧셈(M+1)(N-1)의 연산량이 필요하며, 최소자승법 수행에 O(M)의 연산량이 필요하다. 반면에, 3차원 MUSIC은 공분산행렬 계산에 복소곱셈 과 복소덧셈 의 연산량이 필요하며, 고유치분해 계산에 의 연산량이 필요하다. 이와 더불어 3 차원 MUSIC은 전체 연산량 중에서 대부분을 차지하는 격자탐색에 따른 연산량을 추가적으로 필요로 한다. 따라서 제안한 기법은 3 차원 MUSIC에 비해 요구하는 연산량이 적고 구조가 간단한 장점과 더불어 3 차원 MUSIC과 유사한 음원의 위치추정 성능을 보유하고 있음을 확인할 수 있다.Therefore, the three-dimensional position estimation method of the sound source according to the present invention requires a computation amount of complex multiplication (M + 1) N and complex addition (M + 1) (N-1) to perform correlation and conjugate correlation. O (M) calculation amount is required to perform the square method. On the other hand, three-dimensional MUSIC is complex multiplication for covariance matrix calculation And complex addition Of computations required for eigenvalue decomposition The amount of computation is required. In addition, the 3D MUSIC needs an additional amount of calculation according to the grid search, which takes up most of the total amount of computation. Therefore, the proposed technique has the advantages of less computation and simple structure than 3D MUSIC, and it has similar location estimation performance with 3D MUSIC.
도 1은 본 발명의 실시예에 따른 균일 원 배열과 음원의 기하학적 구조를 나타낸 예시도.1 is an exemplary view showing a uniform circle arrangement and the geometry of the sound source according to an embodiment of the present invention.
도 2는 본 발명의 실시예에 따른 환형배열센서에서 단일 음원의 3차원 위치 추정 방법을 나타낸 순서도.Figure 2 is a flow chart showing a three-dimensional position estimation method of a single sound source in the annular array sensor according to an embodiment of the present invention.
도 3a, 3b, 3c 및 3d는 본 발명의 실시예에 따른 균일 원 배열 센서에서 단일 음원의 3차원 위치 추정 방법과 3차원 MUSIC의 신호대 잡음비에 대한 RMSE를 비교한 예시도.3A, 3B, 3C, and 3D are exemplary diagrams comparing the RMSE of the signal-to-noise ratio of the three-dimensional MUSIC and the three-dimensional position estimation method of a single sound source in a uniform circle array sensor according to an embodiment of the present invention.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020090062620A KR100941697B1 (en) | 2009-07-09 | 2009-07-09 | Apparatus and method for estimating location of single source in uniform circular array |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020090062620A KR100941697B1 (en) | 2009-07-09 | 2009-07-09 | Apparatus and method for estimating location of single source in uniform circular array |
Publications (1)
Publication Number | Publication Date |
---|---|
KR100941697B1 true KR100941697B1 (en) | 2010-02-12 |
Family
ID=42083263
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020090062620A KR100941697B1 (en) | 2009-07-09 | 2009-07-09 | Apparatus and method for estimating location of single source in uniform circular array |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR100941697B1 (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR0175348B1 (en) * | 1995-12-22 | 1999-04-01 | 양승택 | Receive Direction Detection Method |
US20060007043A1 (en) | 2003-06-25 | 2006-01-12 | Jingmin Xin | Method and device for tracking the directions-of-arrival of radio waves |
JP2006234445A (en) | 2005-02-22 | 2006-09-07 | Advanced Telecommunication Research Institute International | Arrival direction estimation apparatus |
-
2009
- 2009-07-09 KR KR1020090062620A patent/KR100941697B1/en active IP Right Grant
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR0175348B1 (en) * | 1995-12-22 | 1999-04-01 | 양승택 | Receive Direction Detection Method |
US20060007043A1 (en) | 2003-06-25 | 2006-01-12 | Jingmin Xin | Method and device for tracking the directions-of-arrival of radio waves |
JP2006234445A (en) | 2005-02-22 | 2006-09-07 | Advanced Telecommunication Research Institute International | Arrival direction estimation apparatus |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6347234B1 (en) | Practical space-time radio method for CDMA communication capacity enhancement | |
US10382230B2 (en) | System and method for channel estimation in mmWave communications exploiting joint AoD-AoA angular spread | |
US20200112360A1 (en) | Method for fast beam sweeping and device discovery in 5g millimeter wave and upper centimeter-wave systems | |
Huang et al. | Frequency-domain AoA estimation and beamforming with wideband hybrid arrays | |
CN104935366B (en) | Beam search method in millimetre-wave attenuator | |
KR100604822B1 (en) | Combined beamforming-diversity wireless fading channel de-modulator using sub-array grouped adaptive array antennas, portable telecommunication receiving system comprising it and method thereof | |
KR100663525B1 (en) | Interference power measurement apparatus and method required space-time beam forming | |
KR101268691B1 (en) | Apparatus and method for receiving data by beamforming in a smart antenna system | |
JP2007159130A (en) | Uplink receiving method and device in distributed antenna mobile communication system | |
EP3387805A1 (en) | Method and computer implemented method for decoding symbol transmitted over mmwave channel, receiver for receiving and decoding symbols transmitted over mmwave channel | |
CN112929962B (en) | Positioning method, positioning device, computer equipment and storage medium | |
Khan et al. | Antenna beam-forming for a 60 Ghz transceiver system | |
Roger et al. | Fast channel estimation in the transformed spatial domain for analog millimeter wave systems | |
CN109411890B (en) | Single-station high-precision three-dimensional positioning method based on multiple antennas | |
Vesa et al. | Direction-of-Arrival estimation for uniform sensor arrays | |
Thazeen et al. | An efficient reconfigurable optimal source detection and beam allocation algorithm for signal subspace factorization | |
KR100941697B1 (en) | Apparatus and method for estimating location of single source in uniform circular array | |
Alam et al. | Design and Performance Analysis of Smart Antenna System for DECT Radio Base Station in Wireless Local Loop. | |
JP2004061468A (en) | Method for estimating arrival direction of multiplex wave using spatial characteristic and receiving beam formation device using the same | |
KHUMANE et al. | Smart Antenna: DOA Estimation and Beamforming using MUSIC & LMS Algorithms | |
El Zooghby et al. | A neural‐network‐based linearly constrained minimum variance beamformer | |
KR20020041558A (en) | Apparatus and Method for forming beam using direction of arrival estimation in a mobile communication system | |
JP2004061467A (en) | Method for estimating arrival direction of period steady signal in multi-path propagation environment and receiving beam formation device using the same | |
Elmurtada et al. | Adaptive smart antennas in 3G networks and beyond | |
KR20130013182A (en) | Smart antenna system for rejection of coherent and incoherent interferences |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A201 | Request for examination | ||
A302 | Request for accelerated examination | ||
E902 | Notification of reason for refusal | ||
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant | ||
FPAY | Annual fee payment |
Payment date: 20130204 Year of fee payment: 4 |
|
FPAY | Annual fee payment |
Payment date: 20140128 Year of fee payment: 5 |
|
FPAY | Annual fee payment |
Payment date: 20150130 Year of fee payment: 6 |
|
FPAY | Annual fee payment |
Payment date: 20160128 Year of fee payment: 7 |
|
FPAY | Annual fee payment |
Payment date: 20170201 Year of fee payment: 8 |
|
FPAY | Annual fee payment |
Payment date: 20180126 Year of fee payment: 9 |
|
FPAY | Annual fee payment |
Payment date: 20200203 Year of fee payment: 11 |