마티외 방정식

1. 개요[편집]

마티외 방정식 Mathieu equation
표준형 (DLMF)	y″ + (a − 2q cos 2t) y = 0
계수 주기	π
유래	É. Mathieu (1868), 타원형 막의 진동
일반화	힐 방정식 (임의 주기계수)
안정성 도구	플로케 이론 · |tr M| > 2
불안정 혀	a ≈ n², n = 1, 2, 3, …
주기해	cen(t, q) · sen(t, q)

마티외 방정식(Mathieu equation)은 복원 강성이 시간에 대해 주기적으로 흔들리는 선형 진동자를 기술하는 2계 상미분방정식으로, 표준형은

$$\frac{d^2y}{dt^2} + \bigl(a - 2q\cos 2t\bigr)\,y = 0$$

이다. 에밀 마티외가 1868년 타원형 막의 진동을 분리변수로 풀다가 각 방향 성분으로 얻은 것이 시초이지만, 오늘날 이 식이 유명한 이유는 막이 아니라 파라메트릭 여기(parametric excitation) 때문이다. 외력이 방정식의 우변에 더해지는 보통의 강제진동과 달리, 여기서는 계수 자체가 흔들린다. 우변이 0이므로 $y\equiv0$은 언제나 해다 — 문제는 그 해가 안정한가이고, 답이 “아니오”인 파라미터 영역이 $(a,q)$ 평면에서 혀 모양으로 뻗어 나온다.

이 문서는 표준형의 관례, 혀가 왜 $a\approx n^2$에서 나오는지, 그리고 이온 트랩·유체 표면·구조물에서 이 식이 어떻게 실물로 나타나는지를 다룬다. 주기계수 선형계의 일반론 — 모노드로미 행렬, 플로케 승수, 랴푸노프-플로케 변환 — 은 플로케 이론 문서가 소스이고, 여기서는 그 도구를 마티외에 실제로 대는 쪽에 집중한다.

2. 표준형의 관례 — 2q 인가 q 인가[편집]

이 방정식은 분야마다 다른 옷을 입고 나타나서, 논문 두 편을 나란히 놓으면 계수가 안 맞는 일이 흔하다. 이 문서는 NIST DLMF §28.2.1의 관례를 따른다.

$$y'' + (a - 2q\cos 2t)\,y = 0$$

앞에 붙은 $2$와 코사인 안의 $2t$가 이 관례의 특징이다. 굳이 이렇게 쓰는 이유가 있다. $\cos 2t$의 주기가 $\pi$라 플로케 해석의 기본 주기가 $\pi$가 되고, 그 대가로 불안정 혀의 꼭짓점이 $a = n^2$ ($n$은 정수)이라는 깔끔한 자리에 앉는다. $2q$의 2는 삼각함수 곱의 반각에서 나오는 $1/2$을 상쇄해 급수 계수를 정수로 만들려는 장치다.

공학 문헌에서는 흔히 다음 형태를 쓴다.

$$\ddot x + (\delta + \varepsilon\cos\tau)\,x = 0$$

플로케 이론 문서가 쓰는 것이 이쪽이다. 두 관례는 $\tau = 2t$ 한 줄로 연결된다.

대응	DLMF	공학 관례
시간	$t$	$\tau = 2t$
계수 주기	$\pi$	$2\pi$
상수항	$a$	$\delta = a/4$
변조 진폭	$2q$	$\varepsilon = -q/2$
혀 꼭짓점	$a = n^2$	$\delta = n^2/4$

$\varepsilon$의 부호는 물리적으로 무의미하다. $\tau$를 반주기만큼 옮기면 코사인 부호가 뒤집히므로, 안정선도는 $q$에 대해 좌우대칭이다. 그래서 보통 $q \ge 0$만 그린다. 계수가 $\cos 2t$가 아니라 임의의 $\pi$-주기 함수로 일반화된 것이 힐 방정식이고, 아래 논의의 대부분은 그쪽에도 그대로 넘어간다.

3. 그네와 거꾸로 선 진자 — 파라메트릭 여기의 얼굴[편집]

가장 익숙한 예가 서서 구르는 그네다. 앉아서 발로 굴리는 것은 외력이지만, 서서 무릎을 굽혔다 펴는 것은 무게중심의 위치, 즉 진자의 실효 길이를 바꾸는 행위다. 길이가 바뀌면 고유진동수 $\sqrt{g/l}$이 바뀌므로 계수가 흔들리는 셈이고, 이때 한 번의 그네 왕복 동안 무릎을 두 번 굽혀야 진폭이 커진다. 구동 진동수가 고유진동수의 두 배 — 그것이 $n=1$ 주공진이다.

물리적으로 더 놀라운 것이 **카피차 진자**다. 지지점을 연직으로 빠르게 흔들면 원래 불안정한 거꾸로 선 평형점이 안정해진다. 이것을 마티외의 언어로 옮겨 보면 안정선도가 왜 그렇게 생겼는지가 정량적으로 확인된다. 길이 $l$, 지지점 변위 $A\cos\omega t$, 위쪽 연직에서 잰 각 $y$에 대한 선형화 방정식은

$$\ddot y - \frac{g}{l}y + \frac{A\omega^2}{l}\cos(\omega t)\,y = 0$$

이고, $t \to \omega t/2$로 시간을 재척도하면 표준형의

$$a = -\frac{4g}{l\omega^2}, \qquad q = \frac{2A}{l}$$

이 된다(변환은 $q$에 마이너스를 붙여 주지만 앞서 말한 대칭성 때문에 부호는 버려도 된다). $a$가 음수라는 것이 “원래 불안정”을 뜻한다. 그런데 $a<0$ 영역이라고 전부 불안정한 것은 아니다. 뒤에 나올 $a_0(q) \approx -q^2/2$ 곡선 위쪽은 안정이며, 그 조건을 대입하면

$$-\frac{4g}{l\omega^2} > -\frac{1}{2}\left(\frac{2A}{l}\right)^2 \quad\Longleftrightarrow\quad (A\omega)^2 > 2gl$$

정확히 이중 평균법이 유효 퍼텐셜로 유도하는 카피차 조건이 나온다. 한쪽은 빠른 진동을 평균해서, 다른 쪽은 마티외 특성곡선의 2차 항에서 — 같은 부등식이 두 경로로 떨어진다.¹

4. 플로케로 읽기 — 2×2 행렬 하나[편집]

계수가 $\pi$-주기이므로 플로케 이론이 그대로 적용된다. 초기조건 $(y,y') = (1,0)$과 $(0,1)$을 각각 한 주기 $[0,\pi]$만 적분해 결과를 열로 세우면 모노드로미 행렬 $M$이 나오고, 장기 거동은 $M^k$가 정한다.

$\operatorname{tr}A = 0$이므로 리우빌 공식에 의해 $\det M = 1$이다. 따라서 특성방정식이

$$\rho^2 - (\operatorname{tr}M)\,\rho + 1 = 0$$

이고 두 플로케 승수는 $\rho, 1/\rho$로 짝을 이룬다. 판정이 여기서 한 줄로 끝난다.

$$\lvert\operatorname{tr}M\rvert < 2 \;\Rightarrow\; \text{안정(유계)}, \qquad \lvert\operatorname{tr}M\rvert > 2 \;\Rightarrow\; \text{불안정(지수 성장)}$$

경계 $\lvert\operatorname{tr}M\rvert = 2$가 정확히 특성곡선이다. 관례상 $\operatorname{tr}M = 2\cos(\pi\nu)$로 두고 $\nu$를 특성지수(characteristic exponent)라 부르는데, $\nu$가 실수면 안정, 복소면 불안정이다. $\nu$가 정수인 자리가 $\lvert\operatorname{tr}M\rvert = 2$, 즉 혀의 경계이며, $\nu$가 짝수면 승수 $+1$($\pi$-주기 해), 홀수면 승수 $-1$($2\pi$-주기 해)이다.

보존계라서 점근 안정은 없다. 승수의 곱이 1이니 둘 다 단위원 안에 들어갈 수 없다. 감쇠가 없는 마티외의 “안정”은 언제나 유계 진동일 뿐 감쇠 수렴이 아니다. 안정선도를 처음 보는 사람이 흰 영역을 “잦아든다”로 오독하는 것이 흔한 실수다.

5. 인스-스트럿 선도 — 혀는 왜 a ≈ n² 에서 나오는가[편집]

$(q, a)$ 평면을 훑으며 위 판정을 칠한 것이 인스-스트럿 안정선도(Ince-Strutt diagram)다.² $q=0$ 축에서 뾰족한 쐐기(혀)들이 $a = 1, 4, 9, 16, \ldots$에 꼭짓점을 두고 위로 뻗는다.

이유는 공진 조건이다. $q=0$이면 해가 $\cos(\sqrt a\,t)$라 고유진동수가 $\omega_0 = \sqrt a$이고, 계수는 진동수 $\omega_d = 2$로 흔들린다. 파라메트릭 공진은 구동 진동수가 고유진동수의 정수분의 2배일 때 일어난다.

$$\omega_d = \frac{2\omega_0}{n} \quad\Longleftrightarrow\quad 2 = \frac{2\sqrt a}{n} \quad\Longleftrightarrow\quad a = n^2$$

$n=1$이 구동 = 고유의 2배인 주공진이다.

혀의 굵기가 급격히 다른 이유도 같은 계산의 연장이다. 해를 $y\sim\cos(nt)$로 두고 섭동항 $-2q\cos2t\cdot y$를 곱해 보면

$$\cos 2t\,\cos nt = \tfrac12\bigl[\cos(n+2)t + \cos(n-2)t\bigr]$$

인데, 되먹임이 원래 진동수 $n$으로 곧장 돌아오려면 $\lvert n-2\rvert = n$, 즉 $n=1$이어야 한다. 주공진만 $q$의 1차에서 자기 자신을 되먹인다는 뜻이다. $n=2$는 두 번, $n=3$은 세 번 튕겨야 돌아오므로 성장률이 각각 $q^2$, $q^3$ 차수로 떨어진다. 실제 특성곡선의 섭동 전개가 이를 확인해 준다.

혀	경계 곡선	폭
$n=0$ (아래쪽)	$a_0 = -q^2/2 + 7q^4/128 - \cdots$	—
$n=1$	$a_1 = 1 + q - q^2/8 - \cdots$, $b_1 = 1 - q - q^2/8 + \cdots$	$2q$
$n=2$	$a_2 = 4 + 5q^2/12 - \cdots$, $b_2 = 4 - q^2/12 + \cdots$	$q^2/2$
$n=3$	$a_3 = 9 + q^2/16 + \cdots$, $b_3 = 9 + q^2/16 - \cdots$	$O(q^3)$

$q=0.1$이면 첫 혀의 폭이 0.2, 둘째 혀는 0.005다. 40배 차이다. 실험에서 주공진만 보이는 이유가 이것이고, 고차 혀는 존재는 하되 잡음과 감쇠에 먼저 먹힌다.

6. 감쇠 — 혀가 공중에 뜬다[편집]

실물에는 감쇠가 있다. 점성 감쇠를 넣은

$$y'' + 2\gamma\,y' + (a - 2q\cos 2t)\,y = 0$$

는 $y = e^{-\gamma t}u$로 놓으면 정확히 감쇠 없는 마티외로 되돌아간다.

$$u'' + \bigl(a - \gamma^2 - 2q\cos 2t\bigr)u = 0$$

근사가 아니라 항등변환이다. 결론이 두 줄로 나온다. 첫째, 감쇠계의 혀는 무감쇠 혀를 $a$축으로 $\gamma^2$만큼 평행이동한 것이다. 둘째, 불안정하려면 무감쇠 성장률이 $\gamma$를 이겨야 한다. 무감쇠 주공진의 성장률이 대략 $q/2$ 정도이므로 문턱은

$$q \gtrsim 2\gamma$$

가 되고, 혀 꼭짓점이 $q$축에서 위로 들려 올라간다. 파라메트릭 여기는 구동이 문턱을 넘어야만 켜지는 현상이라는 실무 감각이 여기서 나온다. 진폭이 0에서 자라나는 불안정이라 선형 주파수응답 함수에는 아예 나타나지 않는다는 점도 같이 기억해 둘 만하다 — 주파수응답만 보고 설계하면 이 모드를 통째로 놓친다.

7. 혀의 경계에서 실제로 일어나는 일[편집]

안정선도의 경계선은 단순한 색깔 경계가 아니라 분기선이다. 승수가 단위원을 뚫는 방식이 그 자리에서 무슨 분기가 일어나는지를 정한다(플로케 이론의 분기 표가 그대로 적용된다).

• 홀수 $n$ 혀($a\approx1, 9, \ldots$): 경계 위의 주기해가 $2\pi$-주기, 즉 계수 주기의 두 배다. 승수가 $-1$을 통과하므로 주기배가 분기다. 파라메트릭 공진에서 응답이 구동의 절반 진동수로 나오는 부조화 응답이 바로 이것의 물리적 얼굴이다.
• 짝수 $n$ 혀($a\approx4, 16, \ldots$): 경계 위의 주기해가 $\pi$-주기라 승수가 $+1$을 통과한다. 조화(harmonic) 응답이고 분기 종류도 다르다.

선형 방정식만 보면 불안정 영역에서 진폭이 무한히 자라지만, 실물은 그렇지 않다. 비선형 항 — 예컨대 진자의 $\sin y$를 $y - y^3/6$까지 편 마티외-더핑 꼴 —

$$y'' + (a - 2q\cos 2t)\,y + \beta y^3 = 0$$

이 붙으면 진폭이 커질수록 유효 고유진동수가 이동해 계가 스스로 공진에서 빠져나온다. 그 결과 불안정한 자명해에서 유한 진폭의 주기해가 갈라져 나와 자리를 잡는다. $\beta>0$이면 초임계 — 진폭이 문턱 너머에서 부드럽게 자라고, $\beta<0$이면 아임계 — 진폭이 튀며 이력(hysteresis)이 생긴다. 그네가 아무리 굴러도 진폭이 유한한 이유, 그리고 파라메트릭 롤이 어느 진폭에서 멈추는 이유가 여기 있다. 선형 마티외는 언제 켜지는지만 말해 주고 얼마나 커지는지는 비선형이 정한다.³

8. 마티외 함수 ce 와 se[편집]

특성곡선 위에서만 주기해가 존재한다. 이 주기해에 붙은 이름이 마티외 함수다.

• $\mathrm{ce}_n(t,q)$ — cosine-elliptic. 짝함수, $a = a_n(q)$ 위에서 존재.
• $\mathrm{se}_n(t,q)$ — sine-elliptic. 홀함수, $a = b_n(q)$ 위에서 존재.

$n$이 짝수면 주기 $\pi$, 홀수면 $2\pi$다. $q\to0$에서 $\mathrm{ce}_n \to \cos nt$, $\mathrm{se}_n \to \sin nt$로 삼각함수에 붙는다. 즉 마티외 함수는 주기적으로 흔들리는 세계의 사인·코사인이고, 특성값 $a_n(q), b_n(q)$가 그 세계의 고유값이다.

원래의 마티외 문제로 돌아가면 의미가 더 분명하다. 타원형 막의 헬름홀츠 방정식을 타원 좌표에서 분리하면 각 방향은 위 방정식(주기 마티외 방정식), 반지름 방향은 $\cos 2t \to \cosh 2t$로 바뀐 변형 마티외 방정식이 된다. 원형 막에서 베셀 함수와 삼각함수가 하는 역할을 타원형 막에서 이 두 짝이 대신한다. 막의 고유진동수가 이산적인 이유가 곧 $a_n(q)$가 이산적인 이유다.

9. 수치적으로 특성값 구하기[편집]

두 갈래가 있다.

푸리에 전개 → 삼중대각 고유값 문제. $\mathrm{ce}_{2m}$처럼 $\pi$-주기 짝함수 해를 $y = \sum_{k\ge0} A_{2k}\cos 2kt$로 넣고 계수를 맞추면 세 항 점화식이 나온다.

$$a A_0 - q A_2 = 0, \qquad (a-4)A_2 - q(2A_0 + A_4) = 0$$ $$(a - 4k^2)A_{2k} - q(A_{2k-2} + A_{2k+2}) = 0 \qquad (k \ge 2)$$

$k=0,1$ 행에만 2가 붙는 것은 $\cos 0 = 1$이라 항이 겹쳐 세어지기 때문이다. $A_0 \to \sqrt2 A_0$로 재척도하면 대칭 삼중대각 행렬이 되고, 적당한 차수에서 잘라 QR 알고리즘에 넣으면 $a_{2m}(q)$가 한꺼번에 나온다. $4k^2$가 $a$와 $q$를 확실히 압도하는 영역부터 계수 $A_{2k}$가 초지수적으로 죽으므로, 절단 차수를 그 지점보다 넉넉히 잡고 값이 더는 변하지 않는지 확인하면 된다. $q$가 커질수록 필요한 차수도 함께 커진다. 홀수 차수와 $\mathrm{se}$ 계열도 사인 전개로 같은 구조가 나온다.

연분수. 같은 점화식을 비 $A_{2k+2}/A_{2k}$에 대해 풀면 $a$를 미지수로 하는 연분수 방정식이 되고, 이를 뉴턴법으로 푼다. 고전적인 블랜치의 방법으로, 특성값 하나만 고정밀로 원할 때 여전히 유용하다.

선도만 필요하면 둘 다 필요 없다. $(q,a)$ 격자를 훑으며 $[0,\pi]$를 룽게-쿠타법으로 두 번 적분해 $\operatorname{tr}M$을 재는 것이 가장 짧다. 특성값의 정확한 자릿수는 못 얻지만 안정·불안정 이진 판정에는 충분하고, 감쇠·비선형 항이 붙어도 코드가 그대로 굴러간다.

10. 어디서 실물로 나타나는가[편집]

• 파울 트랩과 질량분석. 정전기장만으로는 하전입자를 3차원으로 가둘 수 없다(라플라스 방정식이 안장점만 허용한다). 사중극 전극에 직류 $U$와 각진동수 $\Omega$의 RF $V$를 함께 걸면 운동방정식이 정확히 마티외가 되고, $\xi = \Omega t/2$로 두었을 때 $a_x = 8eU/(m r_0^2\Omega^2)$, $q_x = 4eV/(m r_0^2\Omega^2)$다. 두 값 모두 질량에 반비례하므로, $U/V$를 고정한 채 전압을 키우면 서로 다른 $m/z$의 이온이 $(a,q)$ 평면의 같은 직선 위에 질량 순으로 늘어선다. 그 직선을 첫 안정영역의 꼭짓점($a \approx 0.237$, $q \approx 0.706$)에 아슬아슬하게 걸치면 좁은 질량 창만 살아남고 나머지는 발산해 전극에 부딪힌다. 사중극 질량분석기가 바로 이 원리이며, 직선의 기울기가 곧 질량 분해능 손잡이다. RF만 걸면($a=0$) 조건은 $q < 0.908$ 하나로 줄어든다. $q\ll1$인 깊은 안정영역에서의 유사퍼텐셜 근사는 이중 평균법이 다룬다.
• 패러데이 파. 액체를 담은 그릇을 연직으로 진동시키면 어느 진폭부터 표면에 정재파 무늬가 돋는다. 패러데이가 1831년에 보고했고, 벤저민과 어셀이 1954년에 각 표면 모드의 진폭이 감쇠 마티외 방정식을 따른다는 것을 보였다. 응답이 구동 진동수의 절반으로 나오는 부조화(subharmonic) 응답이 주공진의 서명이다. 어느 모드가 먼저 터지는지는 표면장력과 중력이 정하는 분산관계와 혀의 위치가 겹치는 지점이 결정한다.
• 동적 좌굴과 구조물. 축력이 $P_0 + P_t\cos\theta t$로 맥동하는 기둥은 모드별 운동방정식이 마티외-힐 꼴이 되고, $\theta \approx 2\omega_1$에서 횡방향 진동이 폭발한다. 정적 좌굴 하중보다 훨씬 작은 축력으로도 무너질 수 있다는 것이 이 현상의 무서운 점이며, 볼로틴의 1964년 저작이 이 분야의 고전이다. 대형 선박의 파라메트릭 횡동요도 같은 기제인데, 종파를 만나 수선면적이 주기적으로 변하며 복원 강성이 흔들리는 경우다(내항성 참고).
• 주기 매질과 양자계. 계수의 주기가 시간이 아니라 공간이면 같은 식이 블로흐 정리의 무대가 된다. 광격자에 갇힌 냉원자의 운동방정식이 문자 그대로 마티외이고, 그 특성값 $a_n(q)$가 밴드 가장자리, 불안정 혀가 밴드 갭이다. 시간축의 불안정 혀와 공간축의 금지대가 같은 그림이라는 사실은 매번 봐도 신기하다.

11. 관련 문서[편집]

• 플로케 이론 · 힐 방정식 · 블로흐 정리 · 주기경계조건
• 이중 평균법 · 카피차 진자 · 분기 이론 · 극한 순환
• 파울 트랩 · 패러데이 파 · 표면장력 · 슬로싱
• 좌굴 · 내항성 · 모드 해석 · 감쇠 · 정규형 이론
• 고유값 문제 · QR 알고리즘 · 연분수 · 룽게-쿠타법 · 수치 연속법
• 랴푸노프 지수 · 동역학계

12. Footnotes[편집]

1. 두 유도가 겹치는 범위도 정직하게 맞는다. 카피차 근사는 $A\ll l$, $\omega\gg\sqrt{g/l}$을 요구하는데 이는 각각 $q\ll1$, $\lvert a\rvert\ll1$이다. 즉 평균화는 원점 근처에서만 유효하고, $q$가 0.5쯤 되면 $a_0(q)$의 4차 항부터 어긋나기 시작한다. 근사법이 자기 유효범위를 정확히 자백하는 드문 사례. ↩

2. 이름은 특성값 표를 계산한 E. L. 인스(1927)와 선도를 널리 퍼뜨린 M. J. O. 스트럿(1928)에서 왔다. 컴퓨터가 없던 시절 이 표를 손으로 만들었다는 사실이 지금 보면 더 인상적이다. 요즘은 노트북에서 열 줄짜리 코드가 몇 초 만에 같은 그림을 그려 주고, 그래서 이 선도를 처음 보는 학생들은 그것이 한때 사람의 일생이 걸린 계산이었다는 것을 잘 모른다. ↩

3. 그래서 “선형 해석으로 불안정이 나왔으니 이 설계는 끝났다”는 결론은 절반만 맞다. 문턱을 넘느냐(선형)와 넘은 뒤 진폭이 견딜 만한가(비선형)는 별개 질문이고, 후자를 답하려면 결국 시간 적분을 돌려야 한다. 반대로 “시뮬 돌려 보니 안 터지던데요”도 위험하다 — 아임계라면 초기 교란 크기에 따라 답이 갈린다. ↩