2025년 12월 14일 일요일

등방성 안테나의 불가능성(Impossibility of Isotropic Antenna)

[경고] 아래 글을 읽지 않고 "등방성 안테나의 불가능성"을 보면 바보로 느껴질 수 있습니다.


[그림 1] 한 방향으로 빗질한 털공(출처: wikipedia.org)

전자기학은 물리학의 일부이고, 물리학은 수학이란 도구로 기술된다. 그리하여 자연계를 설명하는 근원은 수학보다는 물리학일 수 있다. 이런 평범한 관념을 무너뜨리는 재미있는 개념 중에 등방성 안테나(isotropic antenna)가 있다. 등방성 안테나는 모든 방향으로 전자파 전력을 골고루 복사시켜서 방향도(directivity) 혹은 안테나 이득(antenna gain)을 1 혹은 0 dBi로 만든 안테나이다. 그러나 우리가 안테나 설계를 아무리 잘 하더라도 등방성 안테나를 만들 수 없다[1]. 등방성 안테나의 불가능성(impossibility of isotropic antenna)은 물리학이 아니고 벡터 함수(vector function) 혹은 벡터장(vector field)에 대한 수학으로 설명해야 한다. 그래서 물리학을 도와주는 수학이 아니고, 물리학의 한계를 엄격히 제한하는 우리 사고의 틀로서의 수학을 등방성 안테나에서 만날 수 있다. 물론 물리학과 수학의 관계에 대해서는 다양한 관점이 존재할 수 있다. 하지만 등방성 안테나의 불가능성은 전자기학이 속한 물리학보다 벡터(vector)를 다루는 수학이 더 근원적이며 물리학에 앞서 수학이 실존한다는 증거일 것이다.

[그림 2] 스칼라 함수의 권선수(출처: wikipedia.org)

등방성 안테나가 불가능한 이유는 전자기파가 편파(polarization) 혹은 벡터 특성을 가지기 때문이다. 더 원초적 이유는 벡터 함수의 권선수(winding number) 성질 때문이다[2], [3]. 권선수는 스칼라(scalar) 혹은 벡터 함수가 닫힌 곡선(closed curve)의 반시계 방향을 따라 돈 회수라서 항상 정수로 나온다. [그림 2]처럼 스칼라 함수의 권선수는 우리가 선택한 곡선의 회전수로 헤아리면 된다. 닫힌 곡선 $c$에 대한 권선수 $W(c)$는 다음과 같이 정의한다.

                          (1)

여기서 $\phi_c(s)$는 극형식(polar form)으로 결정한 $c$의 연속적 편각(argument)[예를 들어, $2 \pi$를 0으로 바꾸지 않고 계속 증가시킨다.], $s$는 $c$에 대한 호의 길이(arc length), $c$는 $s$ = $0$에서 시작하고 $s$ = $1$에서 마치며, $c$는 닫힌 곡선이라서 $s$ = $0$과 $1$은 3차원 공간 $\mathbb{R}^3$에서 같은 위치를 가리킨다. 특정점 $\bar p$에 대한 권선수는 $W(c; \bar p)$로 표현한다. 이때 $c$의 편각 $\phi_c$는 위치 벡터(position vector) $\bar R(s)$로 헤아린다.

                          (2a)

여기서 $\bar R$ = $\bar r(s) - \bar p$; $\bar r(s)$는 $c$상의 $s$에 대응하는 3차원 점이다. 만약 $\bar p$가 곡선의 바깥에 있으면, 편각은 한 바퀴를 온전히 채우지 못해서 항상 $W(c; \bar p)$ = $0$이다. 점 $\bar p$가 내부점이면, 닫힌 곡선의 모양에 따라 권선수는 $\pm 1, \pm 2$ 등이 될 수 있다.

[그림 3] 벡터 함수의 지표: 초록색 점은 중심점(center or center point), 빨간색 점은 안장점(saddle or saddle point)(출처: wikipedia.org)

[곡선의 지표 이론(index theory on curves)]

스칼라 함수와 다르게 벡터 함수의 권선수는 정의하기가 까다롭다. 식 (1)을 벡터 함수에도 적용하려면, $\phi_c(s)$가 구성되는 원리를 더 확장해야 한다. 고정된 기준 단위 벡터(reference unit vector) $\hat e_r$에 대해 $\phi_c(s)$를 벡터 내적으로 재정의한다.

                          (2b)

여기서 $\hat v(s)$는 벡터 함수의 $s$에 딸려 존재하는 벡터 $\bar v(s)$의 단위 벡터[그림 3에서 화살표]이다. 벡터 함수의 권선수 $W(c; \bar p)$와 [그림 3]에 보인 이산적인 고립된 특이점(isolated singularity)이 만날 때는 주로 지표(index) 혹은 회전 지표(rotation index) $I_c(\bar p)$로 이름 붙인다. 왜냐하면 권선수는 단순한 양이지만 곡선이나 표면의 전체적인 성질을 나타내는 중요한 숫자라서, 밋밋한 권선수보다 형태를 분류하는 지표란 어감이 더 좋기 때문이다. 여기서 지표는 권선수의 일종이므로 정수이며, [그림 3]에 나온 고립된 특이점은 닫힌 곡선 $c$를 아무리 줄여도 지표가 0이 될 수 없는 곡선의 내부점이다.

                          (3)

여기서 $\bar p$는 고립된 특이점 혹은 일반점이다. 연속 벡터 함수(continuous vector function)의 지표는 복소 함수(complex function)와 비슷한 성질을 가지고 있다.

[연속 벡터 함수의 지표(index of continuous vector function) 성질] [4]
단순 닫힌 곡선(simple closed curve) $c$상에서 연속 벡터 함수가 비소멸(nonvanishing)혹은 언제나 0이 아닐 경우는 다음 성질이 성립한다.
(a) 일반점: $I_c(\bar p)$ = $0$
(b) 안장점(saddle or saddle point): $I_c(\bar p)$ = $-n$, $n$은 양수
(c) 안장점이 아닌 영(zero): $I_c(\bar p)$ = $+n$, $n$은 양수
(d) 닫힌 곡선의 변형: $I_c(\bar p)$ = $I_{c'}(\bar p)$
(e) 닫힌 곡선의 분해 혹은 병합: $I_{c_1 + c_2}(\bar p)$ = $I_{c_1} (\bar p)+ I_{c_2}(\bar p)$
(f) 여러 고정점(fixed point)$I_c(\bar p_1, \bar p_2)$ = $I_c(\bar p_1)$ + $I_c(\bar p_2)$, $\bar p_1, \bar p_2$는 고정점
(g) 모든 고정점을 포함하는 궤적(trajectory): $I_c(\bar p)$ = $\sum_{k=1}^m I_c (\bar p_k)$ = $+1$

[명제 (a)의 증명]

일반점 $\bar p$를 내부에 두고 일주하는 닫힌 경로 $c$에서는 편각 $\phi_c$가 한 바퀴 이상을 완전히 돌지 못하여서 권선수는 항상 0이 된다. 여기서 경로상의 벡터 함수는 비소멸이라서 편각이 정확히 정의된다.
혹은 명제 (d)를 쓰면 더 쉽게 증명된다. 곡선 $c$ 안에는 고정점이 없어서 더 작은 곡선 $c'$으로 줄여도 지표가 유지된다. 이를 한없이 반복한 경우, 줄어든 곡선은 거의 $\bar p$만 회전한다. 이 점은 상수처럼 보여서 편각이 바뀌지 않으므로 지표는 항상 0으로 나온다.

[명제 (b)의 증명]
[그림 3]에 나온 빨간색 안장점을 포함하는 닫힌 곡선을 새로 정해서 벡터의 회전 각도를 모두 더하면 항상 $-2 \pi$가 된다. 일반적으로 이러한 회전은 여러 바퀴가 될 수 있으므로 지표는 $-n$으로 표현된다. 여기서 $n$은 0이 아닌 양수이다. 벡터 함수의 안장점(saddle or saddle point)은 방향에 따라 인력(attraction)과 척력(repulsion)이 동시에 존재하는 불안정점(unstable point)이다. [그림 3]의 빨간색 안장점 주변을 보면, 45˚와 -135˚ 방향에서는 벡터가 안장점으로 모이며[인력] 135˚와 -45˚에서는 안장점 바깥으로 벡터가 나간다[척력]. 더 구체적으로 안장점 근방의 벡터 분포를 인력이 작용하는 방향은 ∪로, 척력 방향으로는 ∩로 간주하면, 안장점 근방은 말 안장과 비슷한 모양이며 안장점에서 조금만 벗어나도 굴러떨어진다. 안장점의 대표적 예는 쌍곡선(hyperbola)스러운 $f(x, y)$ = $x^2 - y^2$이다. 이 스칼라 함수에 구배(gradient)를 취하면 벡터 함수는 $\bar \nabla f$ = $2x \hat x - 2y \hat y$로 나온다. 그러면 안장점 $(0, 0)$에서 힘은 0이 되지만 안장점을 벗어나면 방향에 따라 인력과 척력이 동시에 발생한다. 이러한 안장점은 벡터 함수의 영점(zero)이며 고정점(fixed point)이다. 벡터 함수의 영점은 벡터 크기가 0이 되는 점이다. 비슷하게 벡터 함수의 고정점은 벡터가 사라져서(vanishing) 어느 방향으로도 움직일 수 없는 점이다. 그래서 움직임이 없는 고정점은 크기가 0인 영점과 실질적으로 같다. 반대로 고정점이 아닌 일반점은 보통 비고정점(non-fixed point)으로 부른다.

     
[그림 4] 양전하 ($+$)와 음전하 ($-$)의 전기력선(출처: wikipedia.org)

[명제 (c)의 증명]
명제 (b)의 증명처럼 초록색 중심점을 내부에 담도록 닫힌 곡선을 만들어서 권선수를 재면 항상 0보다 큰 $+n$이 나온다. 여기서 중심점(center or center point)은 주변 벡터가 [그림 3]처럼 회전하는 점이다. 안장점이 아닌 영점은 [그림 3]의 초록색 중심점이나 [그림 4]와 같을 수 있다. 어느 경우든지 기준 단위 벡터 $\hat e_r$을 정해놓고 벡터가 돌아가는 각도를 추적한다. 그러면 신기하게도 벡터가 어느 방향으로 회전하든지 혹은 벡터의 원천(source)이나 수채(sink)가 되든지, 모든 영점의 지표는 언제나 $+1$이다. 이런 영점에 [그림 4]의 음전하와 같은 인력이 작용하면 안장점이 아닌 영점은 안정점(stable point)이다. 반면에 양전하의 전기력선처럼 척력이 있는 경우는 불안정점(unstable point)이 된다. 안장점이 아닌 원점의 예시는 (circle) 방정식인 $f(x, y)$ = $x^2 + y^2$으로 만들 수 있다. 이 함수에 구배를 적용한 벡터 함수는 $\bar \nabla f$ = $2x \hat x + 2y \hat y$이다. 이때 $\bar \nabla f$의 영점은 $(0, 0)$이며 벡터의 원천 역할을 한다.

[명제 (d)의 증명]
고정점을 내부에 두면서 변형한 닫힌 곡선 $c'$의 지표를 $I_{c'}(\bar p)$로 둔다. 그러면 우리가 고려하는 연속 벡터 함수의 지표는 끊어지지 않고 $I_{c}(\bar p)$에서 $I_{c'}(\bar p)$로 바뀐다. 하지만 지표는 권선수이므로 정수만 될 수 있다. 연속인 정수 함수는 상수 함수만 가능하므로 곡선을 수정해도 지표는 그대로이다.

[그림 5] 변형을 통한 닫힌 곡선의 분해와 병합

[명제 (e)의 증명]
[그림 5]에 그린 단순 닫힌 곡선 $c$를 또 다른 닫힌 곡선 $c_1, c_2$로 분해할 때는 복소 함수론의 경로 적분(contour integral)에서 유추한 방식을 쓴다. 먼저 명제 (d)를 써서 $c$의 왼쪽과 오른쪽 경로를 조금씩 내부로 넣는다. 이 과정을 지속해서 내부에 들어간 경로가 서로 만나도록 함으로써 $c_1, c_2$로 닫힌 곡선이 나누어지도록 한다. 이때 바뀐  경로에도 고정점이 없어야 한다. 따라서 고정점을 피해 닫힌 곡선 $c$를 [그림 5]처럼 변형해도 지표는 바뀌지 않으며, 닫힌 곡선 $c_1, c_2$의 지표는 $I_{c_1}(\bar p)$와 $I_{c_2}(\bar p)$의 합으로 계산할 수 있다.

[그림 6] 극점(pole)을 피하는 경로 적분의 변형

[명제 (f)의 증명]
명제 (e)의 증명처럼 복소 함수론에서 실마리를 찾아서 [그림 6]처럼 고정점을 둘러싸는 경로를 가지도록 닫힌 곡선을 변형한다. 그러면 명제 (e)로 인해 각 고정점을 둘러싸는 닫힌 경로의 지표 합으로 전체 지표 $I_c(\bar p_1, \bar p_2)$가 결정된다.

[그림 7] 구면에 형성된 벡터 함수의 궤적: 영점 혹은 고정점이 하나(출처: wikipedia.org)

[명제 (g)의 증명]
연속 벡터 함수(continuous vector function)궤적(trajectory)은  [그림 7]과 같이 특정 위치의 벡터가 향하는 방향을 따라 연속적으로 연결한 곡선이다. 이런 궤적을 이해한 상태에서 [그림 3]을 보면 궤적의 특징이 명확히 보인다. 바로 등고선처럼 연결한 선이 [그림 3]에 나오는 벡터 함수의 궤적이다. 이 궤적을 따라가는 단순 닫힌 곡선의 벡터 편각을 모두 추적하면 [그림 3]이든 [그림 7]이든 항상 $2 \pi$이기 때문에, 모든 궤적의 지표는 언제나 $+1$로 나온다. 또한 명제 (f)의 결과에 의해 궤적의 지표는 내부에 있는 여러 고정점으로 분리되어 계산될 수 있다. 예시적으로 [그림 3]을 다시 본다. 초록색 중심점 $\bar p_1, \bar p_2$의 지표는 각각 $I_c(\bar p_1)$ = $I_c(\bar p_2)$ = $1$이고, 빨간색 안장점 $\bar p_3$가 만드는 지표는 $-1$이다. 이를 모두 더하면 $1+1-1$ = $1$이 되므로, 모든 궤적의 지표는 항상 $1$이라는 명제 (g)의 정확성을 재확인할 수 있다.
______________________________

이제 닫힌 곡선의 지표 이론(index theory)이 우리 사고의 범주에 있기 때문에 더 심오한 정리인 털공 정리(hairy ball theorem, HBT)로 나갈 수 있다. [그림 1]의 소개처럼 털공의 털을 모든 영역에서 가지런하게 빗을 수는 없다는 일상적인 관찰을 수학으로 증명한 정리라서 편하게 털공 정리로 부른다.

[털공 정리(hairy ball theorem, HBT)] [2], [3]
구면(spherical surface)에 접한 벡터 함수(tangent vector function)는 반드시 한 곳 이상에서 영점이 생긴다. 
 
[증명]
[그림 7]과 같은 임의의 벡터 함수를 고려한다. 이 벡터 함수는 구의 표면에 접해있는 접선 벡터(tangent vector)의 모임이다. 먼저 이 벡터 함수의 영점은 없다고 가정한다. 이를 활용해서 [그림 7]에 나온 임의의 벡터 하나를 뽑는다. 언제나 영점이 없기 때문에 이 벡터는 방향이 딱 하나로 결정된다. 또한 벡터 함수는 연속이라서 끊기지 않고 이 방향을 계속 따라갈 수 있다. 계속 진행하면 원래 출발한 벡터와 만나서, 이 벡터의 경로는 궤적을 이룬다. 그러면 지표 성질의 명제 (g)에 따라 우리가 편성한 궤적의 지표는 1로 나온다. 지표가 1인 상태에서 지표 성질의 명제 (d)를 이용함으로써 궤적의 둘레 길이를 계속 줄일 수 있다. 왜냐하면 영점이 없기 때문에 닫힌 곡선의 변형은 자유롭기 때문이다. 궤적이 지속적으로 축소되면 결국 한 점에 이를 것이다. 지표가 1이며 한 점을 따라 도는 궤적은 [그림 3]에 보인 중심점 주변부이기 때문에 영점이 반드시 생긴다. 결국 영점이 없다는 전제가 잘못이라서, 구면에 만든 접선 벡터의 연속된 집합은 한 곳 이상에서 영점이 꼭 생겨야 한다.
______________________________

털공 정리 관점으로 [그림 1, 7]을 관찰한다. 구 좌표계(spherical coordinate system)처럼 [그림 1]은 극축(pole axis)에서 영점이 2개 생긴다. 하지만 영점이 반드시 2개일 필요는 없다. [그림 7]과 같이 [그림 1]의 두 영점이 한 점으로 모일 수도 있다. 이 경우 영점은 달랑 1개이다. 그러나 표면에 지표가 1인 영점이 2개 혹은 지표 1인 영점이 모여서 하나가 생기더라도, 연속 벡터 함수가 생성하는 구면의 지표는 항상 $1+1$ = $2$로 정해진다.

[그림 8] 다이폴 안테나(dipole antenna)로 만든 전방향 복사 패턴(omnidirectional radiation pattern)(출처: wikipedia.org)

이제는 모든 방향으로 전자파를 복사하는 전방향 안테나의 불가능성을 증명할 차례이다. [그림 8]처럼 다이폴 안테나(dipole antenna)를 아무리 작게 만들어도 방위각 $\phi$방향으로만 골고루 복사가 되는 전방향 안테나(omnidirectional antenna)만 설계 가능하다. 여분의 각도인 극고도각 $\theta$ 방향에서는 전자파가 나가지 않는 영점 혹은 공점(null)이 존재한다. 털공 정리는 우리가 온갖 노력을 다해도 등방성 안테나를 만들 수 없다고 말한다.

[등방성 안테나의 불가능성(impossibility of isotropic antenna)] [1]
모든 각도로 골고루 복사되는 등방성 안테나는 구현할 수 없다. 
 
[증명: 털공 정리]
안테나를 둘러싸는 구면을 하나 가정한다. 털공 정리에 의해 이 구면상에는 접선 전기장이나 자기장이 0인 지점이 하나 이상 생긴다. 우리는 가상 구면의 반지름을 마음대로 계속 늘릴 수 있기 때문에, 안테나에서 한없이 멀어지더라도 접선 전자기장이 0인 각도는 꼭 생긴다.

[증명: 헬름홀츠의 정리]
안테나 외부에 복사장이 존재하는 구면을 하나 잡는다. 구면 내부의 전자기장은 표면 등가의 원리(surface equivalence principle)를 준용해서 0으로 만든다. 이 구면에 벡터의 영점이 없다는 말은 [그림 3, 4]와 같은 접선 벡터 모양이 생기지 않는다는 뜻이다. 발산(divergence)회전(curl) 연산자 관점에서 [그림 3, 4]의 영점은 각각 회전 및 발산 원천에 해당하므로, 복사 전기장 $\bar E$에 대해 발산과 회전 원천이 없는 $\bar \nabla \times \bar E$ = $0$, $\bar \nabla \cdot \bar E$ = $0$을 도출한다. 그러면 헬름홀츠의 정리(Helmholtz' theorem)에 의해 이 전기장은 상수 벡터 $\bar E$ = $\bar E_0$만 가능하다. 상수 전기장 $\bar E_0$은 전자기파를 만들 수 없기 때문에, 구면 어딘가에 벡터의 발산이나 회전 원천이 반드시 있어야 한다. 이런 원천은 당연히 그 지점에서 [그림 3, 4]처럼 생긴 영점을 만든다.
______________________________

등방성 안테나는 현실에서 구현할 수 없기 때문에 우리가 연구할 때는 최대한 등방성 안테나에 가깝게 만드는 준등방성 안테나(quasi-isotropic antenna)를 지향한다[5], [6].

[참고문헌]
[1] H. F. Mathis, "A short proof that an isotropic antenna is impossible," Proc. IRE, vol. 39, no. 8, p. 970, Aug. 1951.
[2] P. McGrath, "An extremely short proof of the hairy ball theorem," Amer. Math. Monthly, vol. 123, no. 5, pp. 502–503, May 2016.
[3] E. Curtin, "Another short proof of the hairy ball theorem," Amer. Math. Monthly, vol. 125, no. 5, pp. 462–463, May 2018.
[5] J.-H. Kim and S. Nam, "A compact quasi-isotropic antenna based on folded split-ring resonators," IEEE Antennas Wirel. Propag. Lett., vol. 16, pp. 294–297, 2017.
[6] S. I. Hussain Shah, S. M. Radha, P. Park, and I.-J. Yoon, "Recent advancements in quasi-isotropic antennas: A review," IEEE Access, vol. 9, pp. 146296–146317, 2021.

2025년 4월 11일 금요일

배열 안테나의 합성(Synthesis of Array Antenna)

[경고] 아래 글을 읽지 않고 "배열 안테나의 합성"을 보면 바보로 느껴질 수 있습니다.

[확인] 본 페이지는 exp(jωt) 시간 약속을 사용하고 있습니다.


[그림 1] 원소 안테나의 위상을 조정해서 자유롭게 빔 조향(beam steering)(출처: wikipedia.org)

배열 안테나(array antenna)는 원소 안테나(element antenna)의 여기(excitation) 진폭과 위상을 조절해서 주빔(main beam) 방향, 안테나 이득(antenna gain), 부엽 수준(sidelobe level, SLL), 격자엽(grating lobe) 등의 안테나 규격을 만족시킨다[1]. 원소 안테나의 진폭을 마음대로 바꾸기는 어려워서, 주로 원소 안테나의 위상을 위상 천이기(phase shifter)로 변경한다. 왜냐하면 증폭기로 진폭을 바꾸려면 가격이 비싸고 정밀한 제어가 어렵고, 가변 감쇠기(variable attenuator)는 RF 전력을 감소시켜서 선호되지 않기 때문이다. 그리하여 배열 안테나의 기본형으로 여기 진폭을 1로 둔 균일 선형 배열(uniform linear array, ULA)을 선택한다. ULA는 겉보기에 간단해 보이지만, 경우의 수 관점에서 ULA 설계는 매우 까다롭다. 예시로서 통상적인 위상 천이기의 위상 분해능(phase resolution)은 대략 6비트(6-bits) 혹은 5.6˚ 정도임을 고려한다. 그러면 50개 원소를 가진 ULA의 개별 위상을 바꾸는 경우의 수는 $(2^6)^{50}$ $\approx$ $10^{90}$이 나온다. 우주에 존재하는 모든 원자(atom) 개수를 약 $10^{80}$으로 예상하므로, 원소 개수가 $N$ = $50$인 ULA 설계는 우주 속의 모든 원자보다 많은 경우의 수를 가진다. 따라서 배열 안테나 설계는 단순한 최적화로는 달성할 수 없고 지능적 전략과 이에 걸맞는 최고 수준의 최적화 기법이 필요하다.
배열 안테나의 복사 패턴을 합성하는 출발점은 푸리에 변환(Fourier transform)이다. 안테나 가까이에서 측정한 근역장(near-field) $\bar E_\text{nf}(x', y')$을 복사 패턴이 정의되는 원역장(far-field) $\bar E_\text{ff}(\bar r)$로 바꾸는 근역장-원역장 변환(near-field-to-far-field transformation, NF-FFT, NFFFT, 엔에프에프에프티 혹은 니프트)은 푸리에 역변환(inverse Fourier transform)으로 공식화된다.

                  (1)

배열 안테나의 설계에는 식 (1)을 간략화한 형태를 활용한다[12].

                          (2)

여기서 $\zeta$ = $k \cos \theta$, $k$ = $\omega \sqrt{\mu \epsilon}$; $\text{SF}(\theta)$는 연속 원천의 배열 인자(array factor)공간 인자(space factor), $I(z')$ $z$축을 따라 흐르는 전류, $\mathfrak{F}^{-1}[\cdot]$는 푸리에 역변환이다. 배열 안테나의 원소가 이산적으로 나누어진 경우에는 푸리에 변환 대신 이산 푸리에 변환(discrete Fourier transform, DFT)을 도입한다.

                          (3a)

여기서 $\zeta_m$ = $m \Delta \zeta$ = $k \cos \theta_m$, $z_n'$ = $n \Delta z'$, $\Delta \zeta$ = $2 \pi \mathbin{/} (N \Delta z')$; $\Delta z'$는 원소 안테나 사이의 고정된 간격이다. 우리가 얻기 원하는 배열 인자 $\text{AF}(\theta_m)$이 정해지면, DFT의 원리를 이용한 식 (3a)의 둘째식으로 $I(z_n')$을 결정한다. 위치 $z'$와 파수 $\zeta$의 원점의 급수 중앙에 두고 싶을 때는 대칭 이산 푸리에 변환(symmetric DFT, SDFT)을 이용한다.

                          (3b)

여기서 $N_e$ = $\lfloor N/2 \rfloor$, $N_o$ = $\lfloor (N+1)/2 \rfloor$, $\zeta_m$ = $[m+0.5 \Delta_e(N)] \Delta \zeta$ = $k \cos \theta_m$, $z_n'$ = $[n+0.5 \Delta_e(N)] \Delta z'$; $\lfloor x \rfloor$은 바닥 함수(floor function), $\Delta_e(N)$은 짝수 판별식(discriminant of even number)이다. 따라서 원소 안테나의 개수 $N$이 짝수이든 홀수이든 식 (3b)의 둘째식으로 개별 원소 안테나의 여기 계수(excitation coefficient) $a_n$ = $I(z_n')$을 도출할 수 있다.


   1. 쉘쿠노프 다항식 방법(Schelkunoff polynomial method)    [2]

배열 인자(array factor, AF)의 영점(zero) 혹은 무효화(nulling)를 활용하는 방법으로 쉘쿠노프 다항식 방법(Schelkunoff polynomial method)이 있다. 이 방법은 특정 방항에서 오는 신호 간섭(signal interference)을 완전히 없애기 위한 영점 조향 혹은 널 조향(null steering)에 자주 쓰인다. 쉘쿠노프Sergei Alexander Schelkunoff(1897–1992)는 1936년쉘쿠노프 39세, 일제 식민지 시절에 쉘쿠노프 등가 원리(Schelkunoff equivalence principle)를 제안한 수학자 출신 전파 공학자이다.
쉘쿠노프 다항식 방법은 대수학의 기본 정리(fundamental theorem of algebra)에 기반을 두고 있다. 우리가 $M$개 각도에서 복사 패턴의 영점을 만들고 싶다면, 이 위치에서 근을 가지는 다항식 $P(z)$를 배열 인자 $\text{AF}(\theta)$로 활용한다.

                          (1.1)

여기서 $z$ = $e^{j \psi}$, $\psi$ = $kd \cos \theta$, $k$ = $2 \pi \mathbin{/} \lambda$; 원소 안테나 개수는 $N$ = $M+1$, $d$는 원소 안테나의 간격, $z_n$ = $e^{j kd \cos \theta_n}$는 $\theta$ = $\theta_n$에서 복사 패턴이 0이 되는 영점, $a_n$은 $n$번 안테나의 여기 계수(excitation coefficient), $a_M$은 $z^M$의 여기 계수이다. 식 (1.1)에 나온 $P(z)$를 전개한 후, 여기 계수 $a_n$을 각 $z^n$ 항의 계수로 정한다. 쉘쿠노프 다항식 방법은 영점을 설계하는 좋은 접근법이지만 원소 안테나의 진폭까지 조정해야 하는 약점이 있다. 최근에는 위상 배열 안테나(phased array antenna, PAA)처럼 위상만(phase-only) 조정해서 영점을 결정하는 방법이 최적화를 도구로 다양하게 연구되고 있다[3]–[5].
식 (1.1)처럼 배열 인자는 전기장 기준의 복사 패턴인 전기장 패턴(field pattern)을 만들 때 사용된다. 간혹 편리한 최적화를 위해 전기장 패턴인 식 (1.1)을 제곱해서 전력 패턴(power pattern)으로 처리하기도 한다[13].

                          (1.2)

간편성을 위해 식 (1.2)에서 $G_0$ = $0$ dB로 자주 선택한다. 전력 패턴을 결정하는 $|z-z_m|^2$은 더 간단해질 수 있다[13].

                          (1.3)

여기서 $z_m$ = $e^{\alpha_m + j \beta_m}$; $\alpha_m, \beta_m$은 실수이다. 영점 크기를 조절하는 $\alpha_m$은 보통 0이지만, 부엽 수준(sidelobe level)을 제어하기 위해 0에서 벗어날 수 있다. 예를 들면, $\beta_m$ = $\psi$에서 식 (1.3)은 $|z-z_m|^2$ = $(1-e^{\alpha_m})^2$이다. 만약 $\alpha_m$ $\approx$ $0$인 경우, 식 (1.3)의 크기는 $\alpha_m^2$에 근사적으로 비례하여 영점 채우기(null-filling)가 가능해진다.


   2. 이항 배열(binomial array)    [6], [7]

이항 정리(binomial theroem)를 이용해서 배열 안테나의 복사 패턴이 거듭제곱으로 표현되게 함으로써 아예 부엽을 만들지 않는 배열은 이항 배열(binomial array)로 불린다. 각 원소 안테나가 가져야 되는 가중치는 이항 계수(binomial coefficient)와 연결된다.

                          (2.1)

여기서 $\psi$ = $kd \cos \theta$, $a_0$ = $a_M$; 원소 안테나 개수는 $N$, 이웃한 원소 안테나간의 위상차 $\phi_0$은 0인 측면 발사 배열(broadside array)을 가정, 다항식의 차수는 $M$ = $N-1$, 이항 계수 ${}_MC_n$의 성질로부터 $a_n$ = $a_{M-n}$이다. 이항 배열의 배열 인자는 복사 패턴의 영점이 $z_i$ = $-1$로 동일한 식 (1.1)로 간주할 수 있다. 이와 같이 이항 배열은 배열 인자가 간단한 거듭제곱으로 표현되는 장점이 있지만, 각 여기 계수 $a_n$은 이항 계수 ${}_MC_n$을 따르기 때문에 중심과 양끝 원소 안테나의 진폭 차이는 매우 크다.
식 (2.1)에 제시한 배열 인자의 크기는 다음과 같다.

                          (2.2)

식 (2.2)를 활용해서 균일 선형 배열(uniform linear array, ULA)에서 도입한 정규화 배열 인자(normalized array factor)와 비슷한 이항 배열의 정규화 배열 인자를 정의한다.

                          (2.3)

여기서 $\psi$ = $0$에서 $\text{NAF}(\theta)$ = $1$이다. 이항 배열은 균일 선형 배열에 속하기 때문에, 원소 안테나간의 간격 $d$가 $\lambda$를 넘어가면 격자엽(grating lobe)이 반드시 생긴다. 이런 격자엽을 없앨 때는 꼭 $d < \lambda$ 조건을 지켜야 한다. 더불어 부엽(sidelobe) 자체가 없는 이항 배열은 $d \le \lambda/2$인 요건이 필요하다. 왜냐하면 $|\psi/2|$ $\le$ $\pi d / \lambda$ $\le$ $\pi/2$가 되어야 $\theta$에 따라 주빔이 작아지다가 커지지 않기 때문이다.
식 (2.3)을 써서 이항 배열의 반전력 빔폭(half power beamwidth, HPBW)을 구한다.

                          (2.4a)

                          (2.4b)

여기서 $\theta$ = $\pi/2 \pm \theta_\text{3dB}$, $\text{HPBW}$ = $\theta_\text{bw}$ = $2 \theta_\text{3dB}$이다. 식 (2.4b)는 정확하지만 실제로 쓰기는 복잡해서, $M \gg 1$인 조건으로 식 (2.4b)를 근사화한다.

                          (2.4c)

여기서 $y$ = $\cos x$ $\approx$ $1-x^2 / 2$, $x$ = $\cos^{-1} y$ $\approx$ $\sqrt{2 (1-y)}$, $e^x$ $\approx$ $1-x$이다. 이항 배열의 방향도(directivity) $D$를 구하기 위해, 이 배열의 전체 복사 전력(total radiated power, TRP) $P_\text{rad}$를 계산한다.

                          (2.5a)

만약 $\pi d / \lambda \le \pi/2$ 혹은 $d \le \lambda/2$라면 식 (2.5a)는 불완전 베타 함수(incomplete beta function) $B(x; a, b)$로 간략화된다.

                          (2.5b)

불완전 베타 함수의 성질 $B(1; a, b)$ = $B(a, b)$를 쓰기 위해 $d$ = $\lambda/2$를 선택해본다.

                          (2.5c)

여기서 $B(a, b)$는 베타 함수(beta function), $(\cdot)!!$는 이중 계승(double factorial)이다. 그러면 식 (2.5)를 방향도의 정의에 대입해서 $D$를 결정한다.

                          (2.6a)

                          (2.6b: $d$ = $\lambda/2$)

이항 배열의 차수 $M$이 매우 커지면 스털링의 공식(Stirling's formula)을 써서 식 (2.6b)를 더 간단한 모양으로 바꿀 수 있다.

                          (2.6c)

원소 안테나의 간격 $d$가 반파장이 넘어갈 때는 반파장 단위로 식 (2.5a)를 끊고 적분해서 방향도 $D$를 구한다.


   3. 돌프–체비셰프 배열(Dolph–Chebyshev array)    [8]–[10]

복사 패턴의 빔폭(beamwidth)과 부엽(sidelobe)을 동시에 설계할 수 있는 기교로 돌프–체비셰프 배열(Dolph–Chebyshev array)이 널리 알려져 있다. 돌프–체비셰프 배열은 삼각 함수(trigonometric function)쌍곡선 함수(hyperbolic function)를 오가는 체비셰프 다항식(Chebyshev polynomial) $T_n(x)$의 특성을 이용한다. 체비셰프 다항식 $T_n(x)$는 $x$ = $[-1, 1]$에서는 삼각 함수처럼 변해서 함수값은 $-1$~$1$ 사이에 있으며, $x$가 $[-1, 1]$ 범위를 벗어나면 쌍곡선 함수(hyperbolic function)가 되어서 $|T_n(x)|$는 지수 함수처럼 빠르게 커진다.

                          (3.1)

여기서 $x$ = $\cos \theta$이다.

[그림 3.1] 배열 원소의 배치 구조 (a) 홀수 개수 (b) 짝수 개수(출처: [11])

돌프–체비셰프 배열을 구성할 때 필요한 코사인 함수를 만들기 위해, 배열 원소 개수를 짝수와 홀수로 구분해서 배열 인자를 정의한다.

                          (3.2a)

                          (3.2b)

여기서 $\psi$ = $kd \cos \theta$ = $2u$; $a_0'$는 0번 원소 안테나의 여기 계수(excitation coefficient), 나머지 $a_n$은 $\pm n$번 원소 안테나용 여기 계수, 각 원소간의 위상차는 $\theta_0$ = $0$이므로 돌프–체비셰프 배열은 측면 발사 배열(broadside array)이다. 지금까지 과정은 아무나 생각할 수 있고, 이 다음이 중요하다. 변수 $z$ = $z_0 \cos u$ = $z_0 x$로 두고 배열 인자를 체비셰프 다항식으로 대체한다.

                          (3.3a)

                          (3.3b)

여기서 $T_{K}(z_0)$ = $R_0$, $K$ = $N_\text{tot}-1$는 $2N$ 혹은 $2N-1$, $N_\text{tot}$는 원소 안테나의 총개수; $b_m$은 체비셰프 다항식 $T_K(x)$에서 $x^m$의 계수, $R_0$은 부엽 수준(sidelobe level, SLL)의 역수, 부엽 수준은 전기장 기준으로 $\text{SLL}$ = $1/R_0$, $\text{SLL (dB)}$ = $-20 \log_{10} R_0$이다. 식 (3.3)에 나온 $x^{2m}$과 $x^{2m-1}$은 다시 체비셰프 다항식으로 공식화할 수 있다.

                  (3.4)

그러면 여기 계수 $a_n$은 $R_0$, $b_m$, 식 (3.4)의 계수로 완벽하게 표현된다.

[그림 3.2] 돌프–체비셰프 배열의 복사 패턴

[그림 3.2]는 식 (3.3)으로 선택한 배열 인자가 만드는 복사 패턴을 보여준다. 측면 발사 배열의 주빔 방향(main beam direction)인 $\theta$ = $90^\circ$에서 $u$ = $0$이 되어서  $x$ = $1$, $z$ = $z_0$가 나온다. 이 경우에  $T_{K}(z_0)$ = $R_0$으로 설정하기 때문에 다른 각도보다 값이 훨씬 큰 빔이 나오고 마땅히 주빔이 된다. 또한 주빔을 제외하고 큰 복사 패턴을 만드는 부엽(sidelobe)의 크기는 체비셰프 다항식의 특성으로 인해 1이다. 그래서 돌프–체비셰프 배열의 부엽 수준은 $\text{SLL}$ = $1/R_0$으로 계산된다. 부엽 수준의 역수 $R_0$를 알 때, 필요한 $z$ = $z_0$는 역쌍곡 함수(inverse hyperbolic function)로 유도한다.

                  (3.5a)

                          (3.5b)

여기서 $R_0 \ge 1$, $K$ = $N_\text{tot}-1$이다. 다만 $1/R_0$가 되기 위해서는 [그림 3.2]처럼 모든 부엽의 크기는 항상 1이어야 한다. 그래서 원소 간격 $d$가 가질 수 있는 최대값 $d_\max$가 제한된다.

                  (3.6)

식 (3.5)와 비슷한 과정을 거쳐서 HPBW, $\theta_\text{bw}$ = $2 \theta_\text{3dB}$도 명시적으로 공식화한다.

                  (3.7a)

                          (3.7b)

식 (3.7b)는 정확하지만 계산하기가 번거로워서 $R_0 \gg 1$ 및 $K \gg 1$로 가정해서 $\theta_\text{bw}$를 근사화한다[10].

                  (3.7c)

여기서 $\cosh^{-1} x$ $\approx$ $\log (2x)$, $\cosh x$ $\approx$ $1 + x^2/2$, $\cos^{-1} x$ $\approx$ $\sqrt{2(1-x)}$, $\sin^{-1} x$ $\approx$ $x$; $L$ = $Kd$는 배열 안테나의 전체 길이이다. 측면 발사 배열의 방향도 $D_b$와 반전력 빔폭 $\theta_\text{bw}^b$의 관계식을 활용해서 돌프–체비셰프 배열의 방향도 $D$도 근사화한다.

                          (3.8a)

                          (3.8b)

여기서 $K \gg 1$, $s_\text{hp}$ = 1.39155737825$\cdots$이다. 식 (3.8b)를 정확한 방향도에 더 근접시키기 위해, $K$를 조금씩 바꾼다. 그러면 다음 근사가 식 (3.8b)보다 좀더 나은 결과를 도출한다.

                          (3.8c)

여기서 $N_\text{tot}$ = $K+1$은 원소 안테나의 총개수이다. 이상의 산출물을 종합해서 돌프–체비셰프 배열의 설계법을 확립한다.

[설계법]
  • 목표로 하는 $D$와 $\text{SLL}$을 선택
  • $R_0$ = $1/ \text{SLL}$ 계산
  • 식 (3.8c)로 $N_\text{tot} d$를 얻음
  • 식 (3.6)을 유지하면서 가능한 한 $N_\text{tot}$ = $K+1$은 작게 $d$는 크게 채택
  • 식 (3.5b)로 $z_0$가 정해짐
  • 식 (3.3)을 전개해서 $a_n$을 확정
______________________________

체비셰프 다항식이 있어서 돌프–체비셰프 배열은 설계가 까다롭지만, 근사 공식에 기반한 설계법으로 누구나 여기 계수 $a_n$을 거뜬히 획득할 수 있다.


   4. 우드워드–로슨 방법(Woodward–Lawson method, WLM)    [12]

배열 안테나의 초기 설계를 할 때 주로 사용하는 기법으로 우드워드–로슨 방법(Woodward–Lawson method, WLM)이 있다. 우드워드–로슨 방법은 DFT 설계법인 식 (3b)를 살짝 변형한다. 원래 DFT는 공간 영역(spatial domain)과 스펙트럼 영역(spectral domain)의 표본수가 같아야 한다. 하지만 배열 안테나 설계에서는 이러한 제한이 불편하므로, 원소 안테나의 개수 $N$과 배열 인자의 표본수 $M$을 다르게 해도 DFT 관계가 성립하도록 관계식을 살짝 바꾼다. 공식화를 편하게 하기 위해 $N,M$을 홀수라 가정한다. $N,M$이 짝수일 때는 식 (3b)를 그대로 써도 된다. 일단 배열 인자의 표본점을 $b_m$ = $\text{AF}(\theta_m)$로 둔다.

                          (4.1)

여기서 $N$ = $2N_h + 1$, $\zeta_m$ = $m \Delta \zeta$ = $k \cos \theta_m$, $z_n$ = $n d$, $\Delta \zeta$ = $2 \pi \mathbin{/} L$, $L$ = $Nd$; $a_n$은 $n$번 원소 안테나의 여기 계수, $L$은 배열 안테나의 전체 길이이다. 식 (3b)처럼 $b_m$으로 $a_n$을 얻을 수 있도록 DFT 기반 급수를 만든다. 다만 표본수를 $N$ 대신 $M$으로 대체한다.

                          (4.2)

여기서 $M$ = $2M_h+1$이다. 식 (4.2)가 복사 패턴의 표본 $b_m$을 정확히 형성하는지 알기 위해 식 (4.2)를 식 (4.1)에 대입한다.

                          (4.3)

여기서 $\text{NAF}(m-m')$은 기수 급수(cardinal series)의 성질을 가진 정규화 배열 인자(normalized array factor)이며, $m-m'$ = $0$에서만 $1$이고 나머지는 $0$이다. 결과적으로 식 (4.2)는 식 (4.1)을 잘 생성하므로 우드워드–로슨 방법은 타당하다. 추가적으로 $m$은 $\theta_m$과 연결되어야 하므로, 모든 관측각(observation angle)을 보는 $M_h$의 한계가 유도된다.

                          (4.4)

여기서 $\cos \theta_m > 1$을 만드는 $m$은 버린다. 그러므로 복사 패턴의 표본수 $M$을 증가시키려면, 배열 안테나의 전체 길이 $L$을 파장 $\lambda$보다 크게 늘려야 한다. 다른 관점으로 안테나의 반전력 빔폭(half power beamwidth, HPBW)은 $\Theta$ $\approx$ $\lambda / L$이다. HPBW로 식 (4.4)를 다시 쓰면 $M_h \cdot \Theta$ $\ge$ $1$인 관계도 얻는다. 식 (4.4)는 나이퀴스트 기준(Nyquist criterion)인 $f_s \ge 2 f$와 등가이다.

                          (4.5)

여기서 $v$ = $f \lambda$ = $f_s L_s$; $f_s$와 $f$는 각각 표본화(sampling) 및 신호 주파수, $L_s$는 배열 안테나의 표본화 간격이다.


[참고문헌]
[1] E. K. Miller, "Synthesizing linear-array patterns via matrix computation of element currents,"  IEEE Antennas Propag. Mag., vol. 55, no. 5, pp. 85–96, Oct. 2013.
[2] S. A. Schelkunoff, "A mathematical theory of linear arrays," Bell Syst. Tech. J., vol. 22, no. 1, pp. 80–107, Jan. 1943.
[3] R. L. Haupt, "Phase-only adaptive nulling with a genetic algorithm," IEEE Trans. Antennas Propag., vol. 45, no. 6, pp. 1009–1015, Jun. 1997.
[4] V. Zuniga, A. T. Erdogan, and T. Arslan, "Adaptive radiation pattern optimization for antenna arrays by phase perturbations using particle swarm optimization," NASA/ESA Conference on Adaptive Hardware and Systems, pp. 209–214, 2010.
[5] A. F. Morabito, A. Massa, P. Rocca, and T. Isernia, "An effective approach to the synthesis of phase-only reconfigurable linear arrays," IEEE Trans. Antennas Propag., vol. 60, no. 8, pp. 3622–3631, Aug. 2012.
[6] J. S. Stone, "Directive antenna array," US Patent No. 1643323, Sep. 1927.
[7] P. Apostolov, "An addition to binomial array antenna theory," Prog. Electromagn. Res. Lett., vol. 113, pp. 113–117, 2023.
[9] R. S. Elliott, "Beamwidth and directivity of large scanning arrays," Microw. J., pp. 53–60, Dec. 1963.
[10] C. J. Drane, "Useful approximations for the directivity and beamwidth of large scanning Dolph-Chebyshev arrays," Proc. IEEE, vol. 56, no. 11, pp. 1779–1787, Nov. 1968.
[11] M. Khalaj-Amirhosseini, "Synthesis of low sidelobe level antenna arrays through only main lobe assumption," Sci. Rep., vol. 11, Nov. 2021, art. no. 22797.
[12] C. A. Balanis, Antenna Theory: Analysis and Design, 3rd ed., New Jersey, USA: Wiley-Interscience, 2005.
[13] H.J. Orchard, R.S. Elliott, and G.J. Stern, "Optimising the synthesis of shaped beam antenna patterns," IEE Proc. Pt. H, vol. 132, no. 1, pp. 63 68, Feb. 1985.