KR100335873B1 - 적응배열안테나에서조타행렬및신호행렬검출방법 - Google Patents

적응배열안테나에서조타행렬및신호행렬검출방법 Download PDF

Info

Publication number
KR100335873B1
KR100335873B1 KR1020000002798A KR20000002798A KR100335873B1 KR 100335873 B1 KR100335873 B1 KR 100335873B1 KR 1020000002798 A KR1020000002798 A KR 1020000002798A KR 20000002798 A KR20000002798 A KR 20000002798A KR 100335873 B1 KR100335873 B1 KR 100335873B1
Authority
KR
South Korea
Prior art keywords
matrix
signal
steering
steering matrix
array antenna
Prior art date
Application number
KR1020000002798A
Other languages
English (en)
Other versions
KR20010075887A (ko
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 KR1020000002798A priority Critical patent/KR100335873B1/ko
Publication of KR20010075887A publication Critical patent/KR20010075887A/ko
Application granted granted Critical
Publication of KR100335873B1 publication Critical patent/KR100335873B1/ko

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B30PRESSES
    • B30BPRESSES IN GENERAL
    • B30B15/00Details of, or accessories for, presses; Auxiliary measures in connection with pressing
    • B30B15/28Arrangements for preventing distortion of, or damage to, presses or parts thereof
    • B30B15/285Arrangements for preventing distortion of, or damage to, presses or parts thereof preventing a full press stroke if there is an obstruction in the working area
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V8/00Prospecting or detecting by optical means
    • G01V8/10Detecting, e.g. by using light barriers
    • G01V8/12Detecting, e.g. by using light barriers using one transmitter and one receiver

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

본 발명은 다중 사용자 간섭 환경하에서 송신 신호와 조타 행렬을 결정하는 적응 배열 안테나에서의 알고리즘에 관한 것으로서, 본 발명에서는 변조된 신호의 특성 뿐만 아니라 배열 안테나의 배열(geometry) 성질도 이용하므로써 각 사용자의 조타 벡터간의 직교성을 증대 시키고, 따라서, 빔 형성 능력과 방향 추정 능력에 있어서 향상된 성능을 보인다.

Description

적응 배열 안테나에서 조타 행렬 및 신호 행렬 검출 방법{METHOD FOR DETECTING STEERING MATRIX AND SIGNAL MATRIX IN ADAPTIVE ARRAY ANTENNA}
본 발명은 적응 배열 안테나에 관한 것으로서, 더욱 상세하게는 다중 사용자 간섭 환경하에서의 조타(steer) 행렬 및 신호 행렬을 검출하는 방법에 관한 것이다.
셀룰라 방식의 이동 통신 환경에서 여러 사용자들이 같이 사용하는 경우 간섭 신호 환경하에서 원하는 신호를 찾아서 추정해내는 것이 중요한 문제가 되어 왔다. 간섭 신호 환경하에서 효과적으로 원하는 신호를 찾아내는 일은 셀당 사용자 증가를 가져온다. 최근에 배열 안테나에 적용해 블라인드(blind)방식으로 신호를 추정하는 방법에 관하여 많은 연구가 진행되어 왔는데 블라인드 신호 추정 방식은 신호의 몇가지 특성을 이용하는 방법이다. 이 중에는 수신 신호의 자동 상관성 메트릭스(autocorrelation matrix)를 이용하는 전력 최대화(power maximization) 방법, 디지털 변조된 신호가 원 고정성(cyclostationarity)을 가진다는 사실을 이용하는 LS-SCORE(least square spectral self coherence restoral), 변조된 디지털 통신 신호의 특성을 이용하는 결정 관리(decision-directed) 방식, FSK, PSK, QPSK 등의 변조된 신호가 동일한 엔벨로프(constant envelope)를 가진다는 사실을 이용하는 CMA(constant modulus algorithm) 등이 있는데 그 중에서도 CMA는 디지틀 변조된 신호가 동일한 엔벨로프를 가진다는 사실을 이용하는 방법으로 동일 채널내에 간섭이 있는 상황하에서도 안정적으로 동작 한다.
ILSP(Iterative least square projection) 방식의 알고리즘은 최대 추정(maximum likelihood(ML)) 기준값을 두 단계로 나누어서 최적화하므로 국부적인 최선의 방법(sub-optimum)이라 할 수 있다. 또한 ILSP방식의 알고리즘은 반복하여 최소 제곱(least square) 문제를 해결하는 동시에 추정 방식을 이용한다. ILSP-CMA는 CMA를 실제로 적용하기 위해 나온 방법중 하나로 안테나로 입사하는 신호의 방향을 추정하는 DF(direction finding)뿐만 아니라 BF(beamforming)도 같이 수행하는 장점을 가지고 있다. 하지만, CMA는 신호의 CM(constant modulus) 특성만 이용하기 때문에 조타 벡터(steer vector)들 사이에 종속적(dependent)인 성질이 발생하여 반복 횟수가 증가하면서 잡음 향상(noise enhancement) 효과가 나타나고, 사용자 구분을 모호하게 하는 단점이 있다.
본 발명은 이러한 문제를 해결하기 위한 것으로서 본 발명의 목적은 조타 벡터에 대해서도 배열 안테나의 배열(geometric) 특성에 맞추어 제한을 주므로써 좀더 정확한 DF/BF를 할 수 있는 방법을 제공하는데 있다.
이러한 목적을 달성하기 위하여 본 발명은, X=[x(0),...,x(n),...,x(N-1)]는 M×N 크기의 배열 안테나로 들어오는 신호이고, A=[a(θ1),...,a(θK),...,a(θD)]는 M×D 조타 행렬이고, V=[v(0),...,v(n),...,v(N-1)]는 M×N 크기의 잡음 행렬, S=[s(0),...,s(n),...,s(N-1)]는 D×N 크기의 신호 행렬이며, X=AS+V인 적응 배열 안테나에서 조타 행렬 (A) 및 신호 행렬(S)을 검출하는 방법으로서, ⅰ) 신호 행렬 (S0,) 및 n=0으로 설정하는 단계; ⅱ) a) n을 1 증가시키고을 수행하여 조타 행렬(An)를 검출하는 단계; b) 조타 행렬( An)을 배열 안테나의 배열 특성에 맞게 맵핑하여 추정하는 단계(선형 배열 안테나의 경우 단위원에 매핑하여 추정); c) 추정된 조타 행렬(An)의 각 요소들을 각열의 첫 번째 요소로 정규화 하는 단계; d) 정규화된 조타 행렬(An)을 이용하여을 수행하므로써 신호 행렬(Sn)을 검출하는 단계; ⅲ) 조타 행렬(An과)과 조타 행렬(An-1)이 충분히 가까워질때까지 ⅱ)의 과정을 반복하는 단계를 구비한다.
본 발명은 또한 X=[x(0),...,x(n),...,x(N-1)]는 M×N 크기의 배열 안테나로 들어오는 신호이고, A=[a(θ1),...,a(θK),...,a(θD)]는 M×D 조타 행렬이고, V=[v(0),...,v(n),...,v(N-1)]는 M×N 크기의 잡음 행렬, S=[s(0),...,s(n),...,s(N-1)]는 D×N 크기의 신호 행렬 이며, X=AS+V인 적응 배열 안테나에서 상기 조타 행렬 (A) 및 신호 행렬(S)을 검출하는 방법으로서, ⅰ) 신호 행렬(S0,) 및 n=0으로 설정하는 단계; ⅱ) a) n을 1 증가시키고을 수행하여 조타 행렬(An)를 검출하는 단계; b) 조타 행렬( An)을 배열 안테나의 배열 특성에 맞게 맵핑하여 추정하는 단계(선형 배열 안테나의 경우 단위원에 맵핑하여 추정); c) 추정된 조타 행렬( An)의 각 요소들을 각열의 첫 번째 요소로 정규화 하는 단계; d) 정규화된 조타 행렬( An)을 이용하여을 수행하므로써 신호 행렬(Sn)을 검출하는 단계; e) 검출된 신호 행렬 (Sn)을 가장 가까운 단위원상의 점으로 추정하는 단계; ⅲ) 조타 행렬( An)과 조타 행렬( An-1)이 충분히 가까워질때까지 상기 ⅱ)의 과정을 반복하는 단계를 구비한다.
도 1은 배열 안테나를 이용하는 시스템의 개략 블록도,
도 2는 종래 ILSP-CMA 방식을 설명하는 흐름도,
도 3은 본 발명에 따른 ILSP-SVM 방식을 설명하는 흐름도,
도 4는 본 발명에 따른 ILSP-SVMCMA 방식을 설명하는 흐름도,
도 5는 기준 조타 벡터와 추정된 조타 벡터와의 차이를 비교 도시한 도면,
도 6은 안테나의 수에 따른 SOI의 측정치를 도시한 도면,
도 7은 도 2, 3, 4의 방식에 따른 사용자의 빔 모양을 도시한 도면,
<도면의 주요부분에 대한 부호의 설명>
A1-An 안테나 1 : 처리기
이하, 본 발명을 상세히 설명한다.
먼저, 설명의 편이를 위하여 일반적인 배열 안테나 시스템에 의해 수신된 신호는 다음과 같이 모델링할 수 있다.
도 1은 배열 안테나를 이용하는 시스템의 개략 블록도로서, 도면부호(A1-AM)는 안테나를 의미하며, 도면부호(1)은 배열 안테나의 신호를 처리하기 위한 처리기를 의미한다.
도 1과 같이 M개의 동일한 안테나가 등 간격(d)으로 선형 배열되어 있고, 신호원이 충분히 먼 거리(far-field)에 있다고 가정하면 일차원 구조만 고려하면 되고, 도 1에서 Dis는 신호가 겪는 각 안테나 간의 거리 지연으로 dsinθ 이고, 안테나(A1-AM)간의 시간 지연(τ)은 수학식 1과 같으며, 입사 신호의 파장이이고, 기준 안테나로부터 m 번째 안테나의 거리를 dm,안테나간의 거리를라 하면 dm,은 수학식 2로서 구할 수 있고, N개의 배열중 m번째 안테나에서 수신한 신호는 수학식 3과 같이 주어진다.
여기서, c는 빛 속도를 의미한다.
여기서,는 k번째 전송원(source)에서 전송된 신호 이고,는 k번째 전송원의 입사각이며,는 m번째 안테나에 첨가는 되는 백색 잡음 신호이다. 이때, 안테나 어레이에 의해 D명의 사용자의 신호가 수신되고 신호 스냅샷을 블록 단위로 처리하는 경우 수신 신호는 수학식 4와 같이 표현할 수 있다.
X=AS+V
여기서, X=[x(0),...,x(n),...,x(N-1)]는 M×N 크기의 배열 안테나로 들어오는 신호이고, A=[a(θ1),...,a(θK),...,a(θD)]는 M×D 조타 행렬이고, V=[v(0),...,v(n),...,v(N-1)]는 M×N 크기의 잡음 행렬, S=[s(0),...,s(n),...,s(N-1)]는 D×N 크기의 신호 행렬이다. 여기서 x(n),a(θk),v(n),s(n)은 수학식 5, 6, 7, 8과 같이 주어진다.
x(n)=[x1(n),x2(n),...,xM(n)]T
v(n)=[V1(n),V2(n),...,VM(n)]T
s(n)=[S1(n),S2(n),...,SD(n)]T
여기서, x(n)은 n번째 시간에 배열 안테나로 수신한 신호 벡터이고, a(θk)는 D개의 신호중 k번째 신호의 입사 방향을 나타내는 조타 벡터이고, v(n)은 n번째 시간에 수신되는 신호에 더해지는 백색 가우시안 잡음 벡터이며, s(n)은 n번째 시간에 송신되는 여러 사용자들의 벡터 신호이다.
상술한 모델링에서 입사 신호(X)로부터 잡음 성분을 제외한 A와 S를 추정하는 ML(maximum likelihood)방식의 추정 방법을 이하에서 설명한다.
ML 추정 방식은 신호 구조 모델(signal structure model)을 이용하는 서브 스패이스(subspace)방식의 알고리즘에 비해 신호 사이에 상관성이 클 때 비교적 향상된 성능을 보여준다. 신호의 수는 알고 있다고 가정하면, 배열 안테나로 수신한 신호의 로그 가능성(log likelihood) 함수는 수학식 9와 같으며, 수학식10을 만족하는 A와 S를 찾으면 된다.
여기서,은 노이즈 전력이며,은 제곱 프로베니어스 식(squared Frobenius norm)이다. 여기서 비용 함수를 수학식 11과 같이 두고 비용 함수를 ILSP-CMA에 의해 최소화 한다.
ILSP-CMA는 수학식 11을 최소화하는데 효과적인 알고리즘으로, 이를 도 2를 참조하여 상세히 설명한다. 도시된 바와 같이 n을 0로 설정하고, 이때의 An(n이 0이므로 결국 A0)에 대한 임의의 초기 추정치를 설정한다.(S1)
이후, n을 1 증가시키고, Sn을 수학식 12로서 산출한다.(S2, S3)
여기서, H는 Hermitian 변환을 의미한다.
즉, 수학식 12는 이전 An-1값에 의하여 Sn을 산출함을 의미한다. 그리고, 검출된 Sn을 자기 자신의 크기로 나누어 주므로써 단위원상의 가장 가까운 점으로 추정한다.(S4)
그리고, 단계(S4)에서 추정된 Sn을 이용하여 최소 제곱 방법인 수학식 13을 행하므로 최초의 추정치 A0보다 나은 추정치(An)를 얻을 수 있다(S5).
이후에, An의 각 요소들을 각 열의 첫 번째 요소로 정규화하고(S6), 이 An의 값이 이전 An-1의 값과의 차가 소정 임계값 이하가 될 때까지 상술한 단계(S2-S7)들을 연속하여 수행한다(S7). 즉, An의 값이 이전 An-1의 값과의 차가 소정 임계값 이하이면, 구해진 An및 Sn으로 수신 신호와 방향을 결정하는 것이다(S8).
즉, 종래의 방법에서는 A에 대한 초기 추정치로 A'를 먼저 임의로 선택해서 f(A',S)를 최소화 하는 S를 구한다. 여기서 구해진 S를 자기 자신의 크기로 나누어 주어 단위원상의 신호로 바꾼후, S를이용해서 다시 비용 함수 f(A,S)를 최소화 하는 A를 구하고 다시 A를 대입해서 문제를 해결하는 식을 반복하므로써 A와 S를 향상시키며, 이 과정은 A가 어떤 값으로 수렴할 때까지 수행한다. 이를 간략하게 정리하면 다음과 같다.
1.A0,n=0 설정
2. n=n+1
[Sn] 을 가장 가까운 단위원상의 점으로 추정한다.
An의 각 요소들을 각열의 첫 번째 요소로 정규화 한다.
3. An과 An-1이 충분히 가까워질때까지 2의 과정을 반복 한다
그러나, ILSP-CMA방법은 상술한 바와 같이 신호의 CM 특성만 이용하기 때문에 조타 벡터(steer vector)들 사이에 종속적(dependent)인 성질이 생기게 되어 반복 횟수가 증가하면서 잡음 향상(noise enhancement) 효과가 나타나고, 사용자 구분을 모호하게 하는 단점이 있다.
본 발명자는 이러한 문제를 해결하기 위하여 ILSP-SVM(iterative least square projection based steering vector mapping), ILSP-SVMCMA(iterative least square projection based steering vector mapping and constant modulus approach)방식을 발명하였으며, 본 발명은 ILSP-CMA에 비해 향상된 성능을 보여주며 사용자 사이에 주파수를 공유해서 쓰는 CDMA 시스템에 적용하기에 적합하다.
먼저, ILSP-SVM를 설명하면 다음과 같다.
이 방식에서도 종래의 방식과 마찬가지로 수학식 11을 최소화한다. 그러나 여기서는 배열 안테나의 배열(geometric)특성을 이용하여 반복하여 구하는 조타 벡터를 배열 안테나의 배열 특성에 맞게 맵핑(mapping)한다. 예를 들어 선형 배열 안테나의 경우 입력 신호의 조타 벡터들은 단위 원(unit circle)상에 위치하므로 조타 벡터를 가장 가까운 단위 원상에 맵핑하면 신호가 하나 이상 일 때 각 신호를 구분할 수 있는 능력이 더 커지게 된다. 그래서, ILSP-SVM은 A가 단위원상에 표현된다는 특성을 이용한다. CMA방식은 들어오는 신호의 변조 특성만 이용하므로 두 개의 입사 신호는 대개 같은 변조 방식을 사용하기 때문에 두 조타 벡터는 서로 종속적일 가능성이 커지게 된다. 제안하는 방식에서는 S 에 대한 초기 추정치로 S'를 먼저 임의로 선택해서 f(A, S')를 최소화 하는 A를 구한다. 여기서 구해진 A를 배열 안테나의 배열 특성에 맞게 맵핑 한 후 A를 이용해서 다시 비용 함수 f(A,S)를 최소화 하는 S를 구하고 다시 S를 대입해서 문제를 해결하는 식을 반복하여 A와 S를 향상시켜 나간다. 이 과정을 A가 어떤 값으로 수렴할 때까지 수행한다. 이를 도 3의 흐름도를 참조하여 설명하면 다음과 같다.
도시된 바와 같이 n을 0로 설정하고, 이때의 Sn(n이 0이므로 결국 S0)에 대한 임의의 초기 추정치(S')를 설정한다.(S11)
이후, n을 1 증가시키고, Sn을 수학식 14로서 산출한다.(S12, S13)
수학식 14에 의하여 검출된 An을 배열 안테나의 배열 특성에 맞게 맵핑한다. 즉, 각 요소들을 자기 자신의 크기로 나누어 주므로서 단위원상의 가장 가까운 점으로 추정한다.(S14)
이후에, An의 각 요소들을 각 열의 첫 번째 요소로 정규화하고(S15), 수학식 15를 행하므로써 Sn을 검출한다. (S16)
그리고, 단계(S15)에 검출한 An의 값이 이전 An-1의 값과의 차가 소정 임계값 이하가 될 때까지 상술한 단계(S12-S17)들을 연속하여 수행한다(S17). 즉, An의 값이 이전 An-1의 값과의 차가 소정 임계값 이하이면, 구해진 An및 Sn으로 수신 신호와 방향을 결정하는 것이다(S18).
상술한 흐름도를 간략하게 정리하면 다음과 같다.
1. S0, n=0로 설정
2. n=n+1
[An] 을 배열 안테나의 배열 특성에 맞게 맵핑한다. 선형 배열 안테나의 경우는 조타 행렬의 각 요소들을 자기 자신의 크기로 나누어 가장 가까운 단위원상의 점으로 추정한다.
An의 각 요소들을 각열의 첫 번째 요소로 정규화 한다.
3. An과 An-1이 충분히 가까워질때까지 2의 과정을 반복 한다
다음으로 본 발명자가 발명한 다른 방법 즉, ILSP-SVMCMA 에 대하여 설명하면 다음과 같다.
ILSP-SVMCMA에서는 변조된 신호가 가지는 진폭 특성인 CM과 배열 안테나의배열 특성에 맞게 입력 신호의 조타 벡터를 맵핑하는 특성을 동시에 이용한다. 따라서 이 방식에서는 신호의 진폭 특성과 배열 안테나의 배열 특성을 동시에 이용하므로 종래의 방법에 비해 수렴 속도도 빠를뿐 아니라 더 향상된 성능을 지니게 된다. 이 방법에 대한 흐름도가 도 4에 도시되어 있으며, 도시된 바와 같이 도 3의 흐름도와 동일함을 알 수 있다. 다만 단계(S6)의 수행에 의하여 검출된 Sn에 대하여 자기 자신의 크기로 나누어 주어 단위원상의 신호로 바꾸어줌(S7)을 알 수 있다. 이 과정을 간략하게 설명하면 다음과 같다.
1. S0, n=0으로 설정
2. n=n+1
[An] 을 배열 안테나의 배열 특성에 맞게 맵핑한다. 선형 배열
안테나의 경우는 조타 행렬의 각 요소들을 자기 자신의 크기로
나누어 가장 가까운 단위원상의 점으로 추정한다.
An의 각 요소들을 각열의 첫 번째 요소로 정규화 한다.
[Sn]ij을 가장 가까운 단위원상의 점으로 추정한다.
3. An과 An-1이 충분히 가까워질때까지 2의 과정을 반복한다
이 방법의 변형으로 S 와 A를 구하는 순서를 바꾸어도 유사한 특성을 보인다.
본 발명자는 두 사용자의 신호들과 잡음이 섞여 입사하는 신호에서 각 사용자의 신호를 추출해 내는 정도를 측정하기 위해 조타 행렬(steering matrix) A 로부터 요구되는 기준 조타 벡터와 추정된 조타 벡터와의 차이를 측정해보고 신호 복원 추정치 S로부터 SOI(signal over interference ratio)를 측정해 보았다. 도 5 a는 두 사용자가 20도, 50도의 각도로 10(M=10)개의 안테나로 구성된 배열 안테나로 입사하고 기지국으로부터 같은 거리에 있으며, 사용자 사이의 파워가 동일한 경우(P1=P2)에 세가지 알고리즘에 대해 사용자 1의 조타 벡터가 얼마나 정확하게 추정되는지를 보여준다. 특정 방향으로 조타 벡터를 먼저 구성한 것을 기준으로 하고 입사된 신호로 조타 벡터를 추정해 내어 기준 조타 벡터와 추정된 조타 벡터사이의 차이값을 구해본 것이다. 도 5 b는 두 사용자가 기지국으로부터 거리가 다르다는 가정하에 사용자 1의 조타 벡터의 정확도를 추정한 것이다. 그래서, 사용자1의 파워가 사용자2의 파워의 2배(P1=2P2)라는 것을 제외하고, 다른 조건은 도 5 a 경우와 동일하다. 조타 벡터의 정확도는 사용자 신호의 방향을 추정하는데 그대로 반영되므로 두 사용자 신호를 동시에 수신해서 각 사용자를 구분해내는 정도의 척도로 활용할 수 있다. 사용자와 기지국간의 거리가 같은 경우에 추정된 조타 벡터와 기준 조타 벡터의 차이가 두 사용자 (사용자1 (S1)과 사용자2 (S2)) 각각에 대해 구해본 결과 각 사용자 사이의 동일한 신호 파워의 영향으로 성능이 비슷하다. 제안한 두 알고리즘의 신호대 잡음비가 20dB일 때, 조타 벡터의 오류 정도를 보면 매우 작아서 조타 벡터 추정이 상당히 정확하게 이루어진 것을 볼 수 있다.
도 5의 b에서와 같이 기지국에서 각 사용자 사이의 거리가 다른 경우에는 기지국에서 거리가 가까운 사용자의 조타 벡터의 추정이 좀 더 정확하다. 도 5의 조타 벡터의 정확도 비교에서도 알 수 있는 바와 같이 ILSP-SVM이나 ILSP-SVMCMA방식은 조타 벡터에 제한을 주므로 배열 안테나 구성에 대한 정보도 이용할 수 있게 되어 신호의 특성만을 이용하는 ILSP-CMA방식에 비해 사용자들 사이의 조타 벡터가 각각의 신호들 사이에 같아지게 될 가능성이 작아지므로 좀 더 정확한 입사신호의 방향 추정이 가능한 것이다.
사용자 신호를 추출해내는 정도를 측정하는 또 다른 방법인 SOI에 관해 살펴본다. 추정된 조타 벡터, 추정된 신호, 입사 신호를 이용해서 각 사용자의 신호세기와 간섭 신호 세기를 측정해 보았다. 도 6에 20도, 50도의 입사각을 가지는 두 신호가 기지국으로 입사하는 경우 도 6 a는 기지국에서 각 사용자 사이의 거리가 같을 때 사용자1 의 안테나 수에 따른 SOI 측정치를 신호대 잡음비가 20dB일 때 알고리즘에 따라 나타내었고, 도 6 b는 기지국에서 각 사용자 사이의 거리가 달라서 사용자1의 파워가 사용자2의 파워의 2배 라는 것을 제외하곤 다른 조건은 도 6 a의 경우와 동일 하고, 안테나 수에 따른 SOI 측정치가 신호대 잡음비가 20dB일 때 알고리즘에 따라 나타내었다. 안테나의 수가 어느 정도 감소해도 간섭 신호를 제거할 수 있는 정도는 안테나 수의 감소에 따라 계속해서 급격히 저하되지 않으며 제안된 알고리즘의 간섭 신호 제거 능력도 기존의 ILSP-CMA와 비교해 볼 때 상당히 좋으므로 하나의 기지국으로 여러 사용자들이 억세스(access)하는 실제 BS시스템에 적합하다.
도 7은 사용자1이 20도로 입사하고 사용자2가 50도로 배열 안테나로 입사하며 기지국에서 각 사용자까지의 거리는 같다고 가정하고 신호대 잡음비가 20dB일 때, 알고리즘에 따라 각 사용자에 형성되는 조타 벡터를 이용해서 10개의 안테나로 이루어진 배열 안테나의 빔을 사용자1에 대해 그려본 것이다. 도 7의 ILSP-CMA에 의해 형성된 빔 모양을 보면 사용자1의 조타 벡터가 사용자2의 신호인 50도 부분 에서도 ILSP-SVM이나, ILSP-SVMCMA와 달리 상당한 크기의 이득을 가지게 되므로 사용자 사이의 구분 능력이 상당히 떨어짐을 알 수 있다. 사용자2의 조타 벡터도 사용자1의 조타 벡터와 비슷한 경향을 가지고 있다.
이와 같이 본 발명에서는 조타 벡터에 제한을 주는 방식을 채택하므로 ILSP-CMA가 갖는 단점을 극복하면서 신호의 CMA 특성 또한 이용할 수 있게 되어 향상된 성능을 갖는다는 효과가 있다.

Claims (7)

  1. X=[x(0),...,x(n),...,x(N-1)]는 M×N 크기의 배열 안테나로 들어오는 신호이고, A=[a(θ1),...,a(θK),...,a(θD)]는 M×D 조타 행렬이고, V=[v(0),...,v(n),...,v(N-1)]는 M×N 크기의 잡음 행렬, S=[s(0),...,s(n),...,s(N-1)]는 D×N 크기의 신호 행렬이며, 상기 X=AS+V인 적응 배열 안테나에서 상기 조타 행렬 (A) 및 신호 행렬(S)을 검출하는 방법으로서,
    ⅰ) 상기 신호 행렬(S0,) 및 n=0으로 설정하는 단계;
    ⅱ) a) 상기 n을 1 증가시키고 하기식을 수행하여 조타 행렬(An)를 검출하는 단계;
    b) 상기 조타 행렬( An)을 배열 안테나의 배열 특성에 맞게 맵핑하여 추정하는 단계;
    c) 상기 추정된 조타 행렬(An)의 각 요소들을 각열의 첫 번째 요소로 정규화하는 단계;
    d) 상기 정규화된 조타 행렬(An)을 이용하여 하기식을 수행하므로써 신호 행렬(Sn)을 검출하는 단계;
    ⅲ) 상기 조타 행렬(An)과 조타 행렬(An-1)이 충분히 가까워질때까지 상기 ⅱ)의 과정을 반복하는 단계를 구비하는 적응 배열 안테나에서 조타 행렬 및 신호 행렬 검출 방법.
  2. X=[x(0),...,x(n),...,x(N-1)]는 M×N 크기의 배열 안테나로 들어오는 신호이고, A=[a(θ1),...,a(θK),...,a(θD)]는 M×D 조타 행렬이고, V=[v(0),...,v(n),...,v(N-1)]는 M×N 크기의 잡음 행렬, S=[s(0),...,s(n),...,s(N-1)]는 D×N 크기의 신호 행렬이며, 상기 X=AS+V인 적응 배열 안테나에서 상기 조타 행렬 (A) 및 신호 행렬(S)을 검출하는 방법으로서,
    ⅰ) 상기 신호 행렬(S0,) 및 n=0으로 설정하는 단계;
    ⅱ) a) 상기 n을 1 증가시키고 하기식을 수행하여 조타 행렬(An)를 검출하는 단계;
    b) 상기 추정된 조타 행렬(An)의 각 요소들을 각열의 첫 번째 요소로 정규화하는 단계;
    c)상기 정규화된 조타 행렬(An)을 배열 안테나의 배열 특성에 맞게 맵핑하여 추정하는 단계;
    d) 상기 정규화된 조타 행렬(An)을 이용하여 하기식을 수행하므로써 신호 행렬(Sn)을 검출하는 단계;
    ⅲ) 상기 조타 행렬(An)과 조타 행렬(An-1)이 충분히 가까워질때까지 상기 ⅱ)의 과정을 반복하는 단계를 구비하는 적응 배열 안테나에서 조타 행렬 및 신호 행렬 검출 방법.
  3. X=[x(0),...,x(n),...,x(N-1)]는 M×N 크기의 배열 안테나로 들어오는 신호이고, A=[a(θ1),...,a(θK),...,a(θD)]는 M×D 조타 행렬이고, V=[v(0),...,v(n),...,v(N-1)]는 M×N 크기의 잡음 행렬, S=[s(0),...,s(n),...,s(N-1)]는 D×N 크기의 신호 행렬이며, 상기 X=AS+V인 적응 배열 안테나에서 상기 조타 행렬 (A) 및 신호 행렬(S)을 검출하는 방법으로서,
    ⅰ) 상기 신호 행렬(S0,) 및 n=0으로 설정하는 단계;
    ⅱ) a) 상기 n을 1 증가시키고 하기식을 수행하여 조타 행렬(An)를 검출하는 단계;
    b)상기 조타 행렬(An)을 배열 안테나의 배열 특성에 맞게 맵핑하여 추정하는 단계;
    c) 상기 추정된 조타 행렬(An)의 각 요소들을 각열의 첫 번째 요소로 정규화하는 단계;
    d) 상기 정규화된 조타 행렬(An)을 이용하여 하기식을 수행하므로써 신호 행렬(Sn)을 검출하는 단계;
    e) 상기 검출된 신호 행렬 (Sn)을 가장 가까운 단위원상의 점으로 추정하는 단계;
    ⅲ) 상기 조타 행렬(An)과 조타 행렬(An-1)이 충분히 가까워질때까지 상기 ⅱ)의 과정을 반복하는 단계를 구비하는 적응 배열 안테나에서 조타 행렬 및 신호 행렬 검출 방법.
  4. X=[x(0),...,x(n),...,x(N-1)]는 M×N 크기의 배열 안테나로 들어오는 신호이고, A=[a(θ1),...,a(θK),...,a(θD)]는 M×D 조타 행렬이고, V=[v(0),...,v(n),...,v(N-1)]는 M×N 크기의 잡음 행렬, S=[s(0),...,s(n),...,s(N-1)]는 D×N 크기의 신호 행렬이며, 상기 X=AS+V인 적응 배열 안테나에서 상기 조타 행렬 (A) 및 신호 행렬(S)을 검출하는 방법으로서,
    ⅰ) 상기 신호 행렬(S0,) 및 n=0으로 설정하는 단계;
    ⅱ) a) 상기 n을 1 증가시키고 하기식을 수행하여 조타 행렬(An)를 검출하는 단계;
    b) 상기 조타 행렬(An)의 각 요소들을 각열의 첫 번째 요소로 정규화하는 단계;
    c)상기 정규화된 조타 행렬(An)을 배열 안테나의 배열 특성에 맞게 맵핑하여 추정하는 단계;
    d) 상기 추정된 조타 행렬(An)을 이용하여 하기식을 수행하므로써 신호 행렬(Sn)을 검출하는 단계;
    e) 상기 검출된 신호 행렬 (Sn)을 가장 가까운 단위원상의 점으로 추정하는 단계;
    ⅲ) 상기 조타 행렬(An)과 조타 행렬(An-1)이 충분히 가까워질때까지 상기 ⅱ)의 과정을 반복하는 단계를 구비하는 적응 배열 안테나에서 조타 행렬 및 신호 행렬 검출 방법.
  5. X=[x(0),...,x(n),...,x(N-1)]는 M×N 크기의 배열 안테나로 들어오는 신호이고, A=[a(θ1),...,a(θK),...,a(θD)]는 M×D 조타 행렬이고, V=[v(0),...,v(n),...,v(N-1)]는 M×N 크기의 잡음 행렬, S=[s(0),...,s(n),...,s(N-1)]는 D×N 크기의 신호 행렬이며, 상기 X=AS+V인 적응배열 안테나에서 상기 조타 행렬 (A) 및 신호 행렬(S)을 검출하는 방법으로서,
    ⅰ) 상기 조타 행렬(A0) 및 n=0으로 설정하는 단계;
    ⅱ) a) 상기 n을 1 증가시키고 하기 식을 수행하여 신호 행렬(Sn)을 검출하는 단계;
    b) 상기 검출된 신호 행렬(Sn)을 가장 가까운 단위원상의 점으로 추정하는 단계;
    c)상기 추정된 신호 행렬(Sn)을 이용하여 하기 식을 수행하므로써 조타 행렬(An)을 검출하는 단계;
    d) 상기 검출된 조타 행렬(An)의 각 요소들을 각열의 첫 번째 요소로 정규화하는 단계;
    e)상기 정규화된 조타 행렬( An)을 배열 안테나의 배열 특성에 맞게 맵핑하여 추정하는 단계;
    ⅲ) 상기 조타 행렬(An)과 조타 행렬(An-1)이 충분히 가까워질때까지 상기 ⅱ)의 과정을 반복하는 단계를 구비하는 적응 배열 안테나에서 조타 행렬 및 신호 행렬 검출 방법.
  6. X=[x(0),...,x(n),...,x(N-1)]는 M×N 크기의 배열 안테나로 들어오는 신호이고, A=[a(θ1),...,a(θK),...,a(θD)]는 M×D 조타 행렬이고, V=[v(0),...,v(n),...,v(N-1)]는 M×N 크기의 잡음 행렬, S=[s(0),...,s(n),...,s(N-1)]는 D×N 크기의 신호 행렬이며, 상기 X=AS+V인 적응 배열 안테나에서 상기 조타 행렬 (A) 및 신호 행렬(S)을 검출하는 방법으로서,
    ⅰ) 상기 조타 행렬(A0) 및 n=0으로 설정하는 단계;
    ⅱ) a) 상기 n을 1 증가시키고 하기 식을 수행하여 신호 행렬(Sn)을 검출하는 단계;
    b) 상기 검출된 신호 행렬(Sn)을 가장 가까운 단위원상의 점으로 추정하는 단계;
    c)상기 추정된 신호 행렬(Sn)을 이용하여 하기 식을 수행하므로써 조타 행렬(An)을 검출하는 단계;
    d)상기 검출된 조타 행렬(An)을 배열 안테나의 배열 특성에 맞게 맵핑하여 추정하는 단계;
    e) 상기 추정된 조타 행렬(An)의 각 요소들을 각열의 첫 번째 요소로 정규화하는 단계;
    ⅲ) 상기 조타 행렬(An)과 조타 행렬(An-1)이 충분히 가까워질때까지 상기 ⅱ)의 과정을 반복하는 단계를 구비하는 적응 배열 안테나에서 조타 행렬 및 신호 행렬 검출 방법.
  7. 제 1 항 내지 6항중 어느 한항에 있어서,
    상기 조타 행렬(An)을 배열 안테나의 배열 특성에 맞게 맵핑하여 추정하기 위하여 선형 배열 안테나에서 상기 조타 행렬의 각 요소들을 자기 자신의 크기로 나누어 가장 가까운 단위원상의 점으로 추정하는 적응 배열 안테나에서 조타 행렬 및 신호 행렬 검출 방법.
KR1020000002798A 2000-01-21 2000-01-21 적응배열안테나에서조타행렬및신호행렬검출방법 KR100335873B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020000002798A KR100335873B1 (ko) 2000-01-21 2000-01-21 적응배열안테나에서조타행렬및신호행렬검출방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020000002798A KR100335873B1 (ko) 2000-01-21 2000-01-21 적응배열안테나에서조타행렬및신호행렬검출방법

Publications (2)

Publication Number Publication Date
KR20010075887A KR20010075887A (ko) 2001-08-11
KR100335873B1 true KR100335873B1 (ko) 2002-05-08

Family

ID=19640241

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020000002798A KR100335873B1 (ko) 2000-01-21 2000-01-21 적응배열안테나에서조타행렬및신호행렬검출방법

Country Status (1)

Country Link
KR (1) KR100335873B1 (ko)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100390072B1 (ko) * 2001-02-17 2003-07-04 (주)텔레시스테크놀로지 배열 스마트 안테나 시스템에서 적응 배열 스마트안테나의 신호 처리 방법
CN110808760A (zh) * 2019-10-25 2020-02-18 重庆邮电大学 基于ilsp-cma的同步ds-cdma信号伪码序列和信息序列联合盲估计

Also Published As

Publication number Publication date
KR20010075887A (ko) 2001-08-11

Similar Documents

Publication Publication Date Title
US7539273B2 (en) Method for separating interfering signals and computing arrival angles
US20050276361A1 (en) Interference power measurement apparatus and method for space-time beam forming
US7068221B2 (en) Technique for direction-of-arrival estimation without eigendecomposition and its application to beamforming at base station
US7477190B2 (en) Smart antenna beamforming device in communication system and method thereof
Choi et al. Design of an adaptive antenna array for tracking the source of maximum power and its application to CDMA mobile communications
US20060007043A1 (en) Method and device for tracking the directions-of-arrival of radio waves
US20050271016A1 (en) Beam forming apparatus and method for an array antenna system
US20050259006A1 (en) Beam forming apparatus and method using interference power estimation in an array antenna system
US20020057660A1 (en) Method and apparatus for searcher beamforming in CDMA base station system using array antenna
US7391370B2 (en) Method and apparatus for adaptive direction estimation of radio waves
WO2007007390A1 (ja) 到来波数推定方法、到来波数推定装置及び無線装置
CA2426650A1 (en) Method and apparatus for producing a channel estimate with reduced rank
US20020027526A1 (en) Method and device for using array antenna to estimate location of source in near field
CN107317781B (zh) 一种水下通信中多普勒扩展抑制方法
US6411257B1 (en) Antenna angle-of-arrival estimation using uplink weight vectors
US20050259621A1 (en) Space-time multipath searcher
KR100335873B1 (ko) 적응배열안테나에서조타행렬및신호행렬검출방법
EP1443692B1 (en) Method and apparatus for determining an operating condition in a communications system using antenna arrays
US6904107B1 (en) Space-weighted communication path estimation
US6600447B1 (en) Apparatus and method for determining signal direction from an estimated signal medium response for a ray component of a radio signal
US11764857B2 (en) Systems and methods for multiple signal reception using receiver diversity
CN110299982B (zh) 基于有限带宽的Wi-Fi设备的多径信道侧写方法
US6807517B2 (en) Method to estimate the parameters of a propagation channel
Xin et al. Minimum MSE-based detection of cyclostationary signals in array processing
CN117452442B (zh) 一种基于空频自适应处理的高精度卫星导航抗干扰方法

Legal Events

Date Code Title Description
A201 Request for examination
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment
FPAY Annual fee payment
FPAY Annual fee payment

Payment date: 20160330

Year of fee payment: 15

FPAY Annual fee payment

Payment date: 20170330

Year of fee payment: 16

FPAY Annual fee payment

Payment date: 20180329

Year of fee payment: 17

FPAY Annual fee payment

Payment date: 20190327

Year of fee payment: 18