KR101564192B1 - 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법 및 그 장치. - Google Patents
슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법 및 그 장치. Download PDFInfo
- Publication number
- KR101564192B1 KR101564192B1 KR1020150009759A KR20150009759A KR101564192B1 KR 101564192 B1 KR101564192 B1 KR 101564192B1 KR 1020150009759 A KR1020150009759 A KR 1020150009759A KR 20150009759 A KR20150009759 A KR 20150009759A KR 101564192 B1 KR101564192 B1 KR 101564192B1
- Authority
- KR
- South Korea
- Prior art keywords
- obtaining
- matrix
- linear array
- arrival angle
- vector
- Prior art date
Links
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01Q—ANTENNAS, i.e. RADIO AERIALS
- H01Q3/00—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system
- H01Q3/26—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01Q—ANTENNAS, i.e. RADIO AERIALS
- H01Q21/00—Antenna arrays or systems
- H01Q21/06—Arrays of individually energised antenna units similarly polarised and spaced apart
- H01Q21/08—Arrays of individually energised antenna units similarly polarised and spaced apart the units being spaced along or adjacent to a rectilinear path
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
본 발명은 입사하는 신호의 도래각을 추정하는 방법 및 그 장치에 관한 것으로서, 출력신호 열벡터를구하는 제1 단계와, 최초의 추정행렬을 구하는 제2 단계와, 향상된 추정행렬을 구하는 제3 단계와, 2N-1 차원 열벡터 를 정의하는 제4 단계와, 2N-K차원 부분행렬을 추출하는 제5 단계와, K개의 정규직교 벡터를 구하는 제6 단계와, 목적함수를 구한 다음 그 결과를 이용하여 도래각 추정치를 구하는 제7 단계로 구성되는 도래각 추정 방법과, 상기 방법으로부터 도래각을 추정하는 장치로 구성되어, 기존의 저복잡도 고성능 도래각 추정기법에 대비하여 크게 다르지 않은 복잡도를 유지하면서도 훨씬 향상된 도래각 추정 성능을 가능하게 하는 효과가 있다.
Description
본 발명은 입사하는 신호의 도래각을 추정하는 방법 및 그 장치에 관한 것으로서, 보다 상세하게는 등간격 선형 어레이(Uniform Linear Array) 안테나 시스템에 입사하는 신호원들의 도래각(direction-of-arrival, DOA)을 추정하는 방법 및 그 장치에 관한 것이다.
일반적으로 안테나에 입사하는 신호의 도래각을 추정하는 방법은 빔형성 알고리즘을 적용하였다. 그러나, 이러한 방법은 도래각의 추정 성능이 낮아 두 개 이상의 신호가 인접한 방향에서 도래하는 경우에는 정확한 추정이 어렵다. R. O. Schmdit(선행기술문헌 중 비특허문헌 1 참조)는 이러한 분해성능의 한계점을 근본적으로 극복할 수 있는 MUSIC (multiple signal classification) 기법을 제안하였다. 그러나, MUSIC 기법을 이용할 경우, 높은 분해능의 장점에도 불구하고 복잡도가 매우 큰 고유치 분해(Eigenvalue Decomposition)와 같은 과정을 필요로 한다. 이러한 복잡도 문제를 극복하고자, N. Xi와 L. Liping은 2014년 고유치 분해의 과정을 필요로 하지 않으면서도 비교적 높은 성능을 보여 주는 새로운 도래각 추정기법을 제안하였다.(선행기술문헌 중 비특허문헌 2 참조) 그러나, Xi와 Liping의 기법은 빔형성 알고리즘의 성능보다는 훨씬 뛰어나지만, MUSIC 기법의 성능에는 미치지 못한다는 한계가 있었다.
IEEE Trans. Antennas Propagation, Multiple Emitter Location and Signal Parameter Estimation : Schmidt, R. O., 1986, Vol.34, ISSN : 0018-926X, 276-280쪽
IEEE Signal Process. Lett. A Computationally Efficient Subspace Algorithm for 2-D DOA Estimation with L-shpaed Array : Xi, N. and Liping, L., 2014, Vol.21, ISSN : 1070-9908, 971-974쪽
종래의 저복잡도 고성능 추정기법 (비특허문헌 2)에 비하여, 복잡도는 크게 다르지 않으면서도 훨씬 향상된 성능을 지니는 도래각 추정 방법과 장치를 제안하는 것을 목적으로 한다.
목적을 달성하기 위해 출력신호 열벡터를 구하는 제1 단계와, 최초의 추정행렬을 구하는 제2 단계와, 향상된 추정행렬을 구하는 제3 단계와, 2N-1 차원 열벡터를 정의하는 제4 단계와, 2N-K차원 부분행렬을 추출하는 제5 단계와, K개의 정규직교 벡터를 구하는 제6 단계와, 목적함수를 구한 다음 그 결과를 이용하여 도래각 추정치를 구하는 제7 단계로 구성되는 도래각 추정 방법이다.
본 발명의 다른 특징으로는 상기 방법으로부터 도래각을 추정하는 장치이다.
상기한 바와 같이, 기존의 저복잡도 고성능 도래각 추정기법에 대비하여 크게 다르지 않은 복잡도를 유지하면서도 훨씬 향상된 도래각 추정 성능을 가능하게 하는 효과가 있다.
도 1은 기존 발명인 도래각 추정방법 및 이를 이용한 배열 안테나 시스템의 대표도.
도 2는 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 개념도.
도 3은 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 순서도.
도 4는 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 구성요소인 제4 단계의 세부설명도.
도 5는 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 구성요소인 제5 단계의 세부설명도.
도 6은 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 제1 시뮬레이션 결과표.
도 7은 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 제2 시뮬레이션 결과표.
도 2는 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 개념도.
도 3은 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 순서도.
도 4는 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 구성요소인 제4 단계의 세부설명도.
도 5는 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 구성요소인 제5 단계의 세부설명도.
도 6은 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 제1 시뮬레이션 결과표.
도 7은 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 제2 시뮬레이션 결과표.
본 발명은 선형 어레이 안테나 시스템에 입사하는 신호원들의 도래각을 추정하는 방법에 있어서, 출력신호 열벡터 x(t)(20)를 구하는 제1 단계(S10);와, 최초의 추정행렬 (30)을 구하는 제2 단계(S20);와, 향상된 추정행렬 (40)을 구하는 제3 단계(S30);와, 상기 제3 단계(S30)에서 구한 추정치를 이용하여 2N-1 차원 열벡터 r (50)를 정의하는 제4 단계(S40);와, 상기 2N-1 차원 열벡터 r (50)의 2N-K차원 부분행렬 (60)을 추출하는 제5 단계(S50);와, 상기 제5 단계(S50)에서 추출한 2N-K차원 부분행렬 (60)로부터 K개의 정규직교 벡터(70)를 구하는 제6 단계(S60);와, 상기 제6 단계(S60)에서 구한 K개의 정규직교 벡터(70)를 이용하여 목적함수 (80)의 값을 구한 다음 그 결과를 이용하여 도래각 추정치 (90)를 구하는 제7 단계(S70);로 구성되는 것을 특징으로 하는 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법이다.
이하 첨부된 도면을 참고로 상세하게 설명한다.
먼저, 등간격 선형 어레이 안테나의 센서의 개수를 N, 입사하는 신호의 개수를 K라 정의하겠다.
도 2는 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 개념도로서, N개의 안테나 센서(200)가 일정한 간격 d(210)를 두고 설치되어 있다. 이때, 상기 안테나 센서(200)로 입사신호(220)가 들어오는데 상기 입사신호(220)가 입사하는 각도인 입사각(230)을 도 2와 같이 정의한다.
도 3은 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 순서도로서, 먼저 선형 안테나 어레이의 N개 센서들에서 관측되는 출력신호(10) 를 묶어 N차원 출력신호 열벡터 x(t)(20)를 구하는 제1 단계(S10)를 실행한다.
이후 상기 출력신호 열벡터 x (t)(20)의 자기상관행렬의 신호부분, 즉,
상기 <수학식 1>에서 는 행렬 A의 수반행렬을 나타내며, Q는 잡음의 자기상관행렬 (즉, N개의 센서의 잡음들을 열벡터로 구성한 벡터를 w (t) 라 두면 로 정의되는 행렬)을 의미한다.
보다 구체적으로, 자기상관행렬 R(100)의 n번째 행, m번째 열의 원소 R nm 의 최초 추정치 (즉, 최초의 추정행렬 (30)의 n번째 행, m번째 열의 원소)는 시간에서 측정하여 얻은 T개의 안테나 출력신호 벡터 를 이용하여 아래의 <수학식 2>와 같이 구한다.
한편, 상기의 자기상관행렬 R(100)은 에르미트 행렬(Hermitian)이면서 동시에 토에플리츠(Toeplitz) 행렬이 된다. 따라서, 이러한 성질을 이용하여 도 4에 정의된 것과 같이 상기의 자기상관행렬 R(100)의 보다 향상된 추정행렬 (40)을 구하는 제3 단계(S30)를 다음 단계로 실행한다. 구체적으로 설명하면, 상기 제3 단계(S30)에서는 상기 향상된 추정행렬 (40)을 구성하는 원소 을 아래의 <수학식 3>과 같이 구한다.
이후, 상기 제5 단계(S50)에서 추출한 2N-K차원 부분행렬 (60)에 Gram-Schmidt 직교정규화과정(Orthonormalization Process)을 적용하여, 크기가 1이고 서로 직교하는 K개의 정규직교 벡터(70)를 구하는 제6 단계(S60)를 수행한다.
이후, 를 에서 까지 일정간격 또는 임의의 간격으로 변화시켜 가면서 목적함수 (80)에 피크값을 주는 K개의 도래각 추정치 (90)를 구하는 제7 단계(S70)가 수행된다. 상기 도래각 추정치 (90)는 로 표현된다.
도 4는 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 구성요소인 제4 단계(S40)의 세부설명도로서, 추정치로부터 얻어진 향상된 추정행렬 (40)의 첫번째 행과 첫번째 열의 총 2(N-1)개의 원소를 추출하는 방법을 좀 더 자세하게 알 수 있다. 첫번째 행이 부터 까지 N개 이고, 첫번째 열이 부터 까지 N개 이지만 이 반복되기 때문에 총 개수는 2(N-1)개가 된다.
도 5는 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 구성요소인 제5 단계(S50)의 세부설명도로서, 2N-1 차원 열벡터 r (50)의 2N-K차원 부분행렬 (60)을 추출하는 방법을 더 직관적으로 나타내 주고 있다.
도 6, 7은 본 발명인 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법의 제 1, 2 시뮬레이션 결과표로서, 도 6, 7의 x축은 신호대 잡음비(signal to noise ratio)의 dB을 뜻하고, y축은 평균 제곱근 오차(Root Mean Squared Error)를 뜻한다. 제 1, 2 시뮬레이션 결과표를 살펴보면 본 발명의 알고리즘(B)가 상기 선행기술(비특허문헌 2)의 알고리즘(C)보다 훨씬 낮은 신호대 잡음비(signal to noise ratio)에서도 사용이 가능하며 오차가 훨씬 낮음을 알 수 있다.
본 발명은 특정의 실시 예 및 적용 예와 관련하여 도시 및 설명하였지만, 첨부된 특허청구범위에 의해 나타난 발명의 사상 및 영역으로부터 벗어나지 않는 한도 내에서 다양한 개조 및 변화가 가능 하다는 것을 당 업계에서 통상의 지식을 가진 자라면 누구나 쉽게 알 수 있을 것이다.
10. 출력신호 20. 출력신호 열벡터 x (t)
30. 최초의 추정행렬 40. 향상된 추정행렬
50. 2N-1 차원 열벡터 r 60. 2N-K차원 부분행렬
70. K개의 정규직교 벡터
80. 목적함수 90. 도래각 추정치
100. 자기상관행렬 R 200.안테나 센서
210.간격 d 220.입사신호
230.입사각
B. 본 발명의 알고리즘 C. 선행기술의 알고리즘
N. 안테나의 센서의 개수 K. 입사하는 신호의 개수
S10.제1 단계 S20.제2 단계
S30.제3 단계 S40.제4 단계
S50.제5 단계 S60.제6 단계
S70.제7 단계
30. 최초의 추정행렬 40. 향상된 추정행렬
50. 2N-1 차원 열벡터 r 60. 2N-K차원 부분행렬
70. K개의 정규직교 벡터
80. 목적함수 90. 도래각 추정치
100. 자기상관행렬 R 200.안테나 센서
210.간격 d 220.입사신호
230.입사각
B. 본 발명의 알고리즘 C. 선행기술의 알고리즘
N. 안테나의 센서의 개수 K. 입사하는 신호의 개수
S10.제1 단계 S20.제2 단계
S30.제3 단계 S40.제4 단계
S50.제5 단계 S60.제6 단계
S70.제7 단계
Claims (6)
- 선형 어레이 안테나 시스템에 입사하는 신호원들의 도래각을 도래각 추정장치를 통하여 추정하는 방법에 있어서,
등간격 선형 어레이를 구성하는 각 센서들로 얻어지는 신호를 하나의 열벡터로 묶은 출력신호 열벡터 x(t)(20)를구하는 제1 단계(S10);와,
상기 출력신호 열벡터 x(t)(20)의 자기상관행렬의 신호부분, 즉,
의 식을 이용하여 최초의 추정행렬 (30)을 구하는 제 2 단계(S20);와,
의 식을 이용하여 향상된 추정행렬 (40)을 구하는 제3 단 계(S30);와,
상기 제3 단계(S30)에서 구한 추정치를 이용하여 2N-1 차원 열벡터 r (50)를 식을 이용하여 정의하는 제4 단계(S40);와,
상기 2 N-1 차원 열벡터 r (50)의 2N-K차원 부분행렬 (60)을 식을 통하여 추출하는 제5 단계(S50);와,
상기 제5 단계(S50)에서 추출한 2N-K차원 부분행렬 (60)로부터 K개의 정규직 교 벡터 (70)를 구하는 제6 단계(S60);와,
는 행렬 , 는 입사신호의 파장일 경우, 상기 제6 단계(S60)에서 구한 K개의 정규직교 벡터 (70)와 식을 이용하여 목적함수 (80)의 값을 구한 다음 그 결과를 이용하여 도래 각 추정치 (90)를 구하는 제7 단계(S70);로 구성되는 것을 특징으로 하는 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법. - 삭제
- 삭제
- 삭제
- 삭제
- 선형 어레이 안테나 시스템에서 입사하는 신호원들의 도래각을 추정하는 장치에 있어서,
상기 제1 항에 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법으로 도래각을 추정하는 것을 특징으로 하는 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 장치.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020150009759A KR101564192B1 (ko) | 2015-01-21 | 2015-01-21 | 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법 및 그 장치. |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020150009759A KR101564192B1 (ko) | 2015-01-21 | 2015-01-21 | 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법 및 그 장치. |
Publications (1)
Publication Number | Publication Date |
---|---|
KR101564192B1 true KR101564192B1 (ko) | 2015-10-28 |
Family
ID=54429106
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020150009759A KR101564192B1 (ko) | 2015-01-21 | 2015-01-21 | 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법 및 그 장치. |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR101564192B1 (ko) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327292A (zh) * | 2020-10-11 | 2021-02-05 | 南京理工大学 | 一种二维稀疏阵列doa估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002243826A (ja) | 2000-12-12 | 2002-08-28 | Matsushita Electric Ind Co Ltd | 電波到来方向推定装置、電波到来方向推定方法及び指向性可変送受信装置 |
JP2011053056A (ja) | 2009-09-01 | 2011-03-17 | Fujitsu Ltd | 到来方向推定装置及び方法 |
JP2013205399A (ja) | 2012-03-29 | 2013-10-07 | Panasonic Corp | レーダ装置 |
KR101375425B1 (ko) | 2013-02-28 | 2014-03-17 | 충북대학교 산학협력단 | 배열 안테나에서 전파의 도래각 추정 방법 |
-
2015
- 2015-01-21 KR KR1020150009759A patent/KR101564192B1/ko not_active IP Right Cessation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002243826A (ja) | 2000-12-12 | 2002-08-28 | Matsushita Electric Ind Co Ltd | 電波到来方向推定装置、電波到来方向推定方法及び指向性可変送受信装置 |
JP2011053056A (ja) | 2009-09-01 | 2011-03-17 | Fujitsu Ltd | 到来方向推定装置及び方法 |
JP2013205399A (ja) | 2012-03-29 | 2013-10-07 | Panasonic Corp | レーダ装置 |
KR101375425B1 (ko) | 2013-02-28 | 2014-03-17 | 충북대학교 산학협력단 | 배열 안테나에서 전파의 도래각 추정 방법 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327292A (zh) * | 2020-10-11 | 2021-02-05 | 南京理工大学 | 一种二维稀疏阵列doa估计方法 |
CN112327292B (zh) * | 2020-10-11 | 2024-01-09 | 南京理工大学 | 一种二维稀疏阵列doa估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4722132B2 (ja) | 到来波数推定方法、到来波数推定装置及び無線装置 | |
Bazzi et al. | Online angle of arrival estimation in the presence of mutual coupling | |
JP6395677B2 (ja) | 到来方向推定装置 | |
CN107255796B (zh) | 一种非均匀噪声下窄带近场信号源定位方法 | |
KR101958337B1 (ko) | 신호의 도래각을 추정하는 방법 및 장치 | |
CN102385048A (zh) | 基于均匀线阵的混合信号方向估计方法 | |
JP6106571B2 (ja) | 音源位置推定装置、方法及びプログラム | |
KR101559270B1 (ko) | 자기상관행렬을 이용한 슬라이딩벡터 기반의 저계산량 등간격 선형 어레이 도래각 추정 방법 및 그 장치. | |
CN110941980B (zh) | 一种密集环境中基于压缩感知的多径时延估计方法及装置 | |
KR101557259B1 (ko) | 슬라이딩벡터 기반의 등간격 선형 어레이 도래각의 향상된 추정 방법 및 그 장치 | |
KR101564192B1 (ko) | 슬라이딩벡터 기반의 등간격 선형 어레이 도래각 추정 방법 및 그 장치. | |
JP2014173940A (ja) | 分散開口レーダ装置及びレーダ信号処理方法 | |
Albagory et al. | MUSIC 2D-DOA estimation using split vertical linear and circular arrays | |
KR101680742B1 (ko) | 저중복 저계산량 건너뛰기 벡터 슬라이딩 등간격 선형 어레이 도래각 추정 방법 및 그 장치. | |
KR101683827B1 (ko) | 샘플링 기법을 이용한 슬라이딩벡터 기반의 저계산량 등간격 선형 어레이 도래각 추정방법 및 그 장치. | |
TW201342364A (zh) | 於旋積混合盲信號源中將多個訊號源進行分離之方法 | |
KR101181317B1 (ko) | 순방향/역방향 벡터를 이용한 신호의 도래방향 탐지 방법 | |
CN105303009A (zh) | 基于压缩感知与正则mfocuss的超分辨谱估计方法 | |
Das et al. | Performance analysis of TLS-Esprit and QR TLS-Esprit algorithm for Direction of Arrival estimation | |
Vijayamohanan et al. | Convolutional neural networks for radio source detection | |
Cordill et al. | Mutual coupling calibration using the Reiterative Superresolution (RISR) algorithm | |
Du et al. | Three stochastic measurement schemes for direction-of-arrival estimation using compressed sensing method | |
CN110888106A (zh) | 一种角度与频率联合估计的增广doa矩阵方法 | |
Satheesh et al. | Compressive sensing for array signal processing | |
Harley et al. | Matched field processing localization with random sensor topologies |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AMND | Amendment | ||
E601 | Decision to refuse application | ||
AMND | Amendment | ||
X701 | Decision to grant (after re-examination) | ||
GRNT | Written decision to grant | ||
LAPS | Lapse due to unpaid annual fee |