플로케 이론

계산물리 수치해석 마지막 수정: 2026-08-17 04:44:52

1. 개요[편집]

플로케 이론 Floquet Theory
정립	Gaston Floquet (1883)
대상	$\dot{\mathbf x} = A(t)\mathbf x$, $A(t+T)=A(t)$
핵심 대상	모노드로미 행렬 $M = \Phi(T)$
판정	모든 플로케 승수가 단위원 안 ⟹ 점근 안정
자율계 특례	승수 하나는 항상 정확히 1
공간판	블로흐 정리

계수가 주기적으로 흔들리는 선형계는, 한 주기만 제대로 보면 나머지는 거듭제곱이다.

플로케 이론(Floquet theory)은 계수가 시간에 대해 주기적인 선형 상미분방정식

\dot{\mathbf x} = A(t)\,\mathbf x,\qquad A(t+T) = A(t)

의 해 구조와 안정성을 다루는 이론이다. 1883년 가스통 플로케가 정리했고, 요지는 딱 두 줄이다. (1) 이런 계의 해는 “주기적인 것 × 지수적인 것”으로 항상 쪼개진다. (2) 따라서 안정성은 한 주기 동안의 흐름을 나타내는 행렬 하나, 즉 모노드로미 행렬의 고유값만 보면 끝난다.

이게 왜 필요하냐면, 주기계수는 상수계수와 근본적으로 다르게 굴기 때문이다. 상수계수 $\dot{\mathbf x} = A\mathbf x$ 라면 $A$ 의 고유값 실수부만 보면 되지만, $A(t)$ 가 흔들리면 매 순간의 고유값이 전부 좌반평면에 있어도 해가 발산할 수 있다. 마르쿠스와 야마베가 1960년에 내놓은 2×2 예가 유명한데, $A(t)$ 의 고유값이 시간과 무관하게 $-\tfrac14 \pm \tfrac{\sqrt7}{4}i$ 로 고정돼 있는데도 해가 $e^{t/2}$ 로 자란다. “얼린 고유값”(frozen-time eigenvalue)으로 안정성을 판정하면 안 되는 이유이고, 플로케 이론이 존재하는 이유이기도 하다.¹

응용은 두 갈래다. 하나는 계수 자체가 주기적인 물리계 — 매개 여기(parametric excitation), 회전 기계, 주기 구동 양자계. 다른 하나는 극한 순환의 안정성 판정인데, 주기궤도 둘레의 변분방정식이 자동으로 주기계수 선형계가 되기 때문이다. 후자 덕분에 수치 연속법 패키지가 전부 플로케를 내장하고 있다.

2. 기본행렬과 모노드로미[편집]

$\Phi(t)$ 를 $\dot\Phi = A(t)\Phi$ , $\Phi(0) = I$ 를 만족하는 기본행렬(fundamental matrix)이라 하자. 모든 해가 $\mathbf x(t) = \Phi(t)\mathbf x(0)$ 로 적히므로 $\Phi$ 를 아는 것이 곧 계를 아는 것이다.

여기서 핵심 관찰. $A$ 가 $T$ -주기이므로 $\Phi(t+T)$ 도 같은 방정식의 해행렬이다. 두 기본해집합은 상수행렬 하나로 연결되므로

\Phi(t+T) = \Phi(t)\,M,\qquad M \equiv \Phi(T)

이고, 이 $M$ 을 모노드로미 행렬(monodromy matrix)이라 부른다. 위 관계를 반복하면

\Phi(t + kT) = \Phi(t)\,M^{k}

즉 한 주기의 행동만 알면 $k$ 주기 뒤는 $M$ 의 거듭제곱이다. 장기 거동이 $M^k$ 의 거동으로 환원되고, 그건 $M$ 의 스펙트럼 문제다.

$M$ 은 항상 정칙이다. 아벨-야코비-리우빌 공식

\det \Phi(t) = \exp\!\left(\int_0^{t}\operatorname{tr}A(s)\,ds\right) \quad\Longrightarrow\quad \det M = \exp\!\left(\int_0^{T}\operatorname{tr}A(s)\,ds\right) \neq 0

이 그것을 보증한다. 이 식은 나중에 승수들의 곱을 손으로 확인하는 데도 요긴하게 쓰인다.

3. 플로케 정리 — 주기계를 상수계수로 편다[편집]

플로케 정리. $\Phi(t)$ 는 다음과 같이 분해된다.

\Phi(t) = P(t)\,e^{Bt},\qquad P(t+T) = P(t),\quad P(0) = I,\quad e^{BT} = M

여기서 $B = \tfrac{1}{T}\log M$ 은 행렬 로그이고, $P$ 는 $T$ -주기 정칙 행렬함수다. 증명은 $B$ 를 위와 같이 정의하고 $P(t) \equiv \Phi(t)e^{-Bt}$ 로 두면 $P(t+T) = \Phi(t)Me^{-BT}e^{-Bt} = P(t)$ 로 바로 확인된다.

여기서 실용적으로 결정적인 따름정리가 나온다. 변수변환 $\mathbf x = P(t)\mathbf y$ 를 넣으면

\dot{\mathbf y} = B\,\mathbf y

주기계수 방정식이 상수계수 방정식으로 바뀐다. 이것을 랴푸노프-플로케 변환이라 부른다. 즉 주기계수 선형계는 “본질적으로 새로운 것”이 아니라 주기적으로 흔들리는 좌표계에서 본 상수계수계일 뿐이다. 어렵게 느껴지던 대상이 좌표 선택의 문제로 내려앉는다.

단, 실수에서는 잔가시가 하나 붙는다. 실행렬 $M$ 이 음의 실고유값을 가지면 실수 로그가 존재하지 않을 수 있어서 $B$ 가 복소가 된다. 이 경우 주기를 $2T$ 로 늘리면 ( $M^2$ 은 항상 실로그를 가진다) 실수 판 플로케 형이 얻어진다. 주기배가 현상이 이 대목과 무관하지 않다는 것이 이론적으로 재미있는 지점이다.

4. 승수, 지수, 그리고 안정성[편집]

$M$ 의 고유값 $\mu_1,\dots,\mu_n$ 을 플로케 승수(Floquet multiplier), $\mu_i = e^{\lambda_i T}$ 를 만족하는 $\lambda_i$ 를 플로케 지수(Floquet exponent)라 한다. 지수는 $2\pi i/T$ 의 정수배만큼 모호하다 — 로그의 다가성이 그대로 넘어온 것이라 어쩔 수 없다. 그래서 실무에서는 대개 승수를 그대로 쓴다.

판정은 단위원 기준이다.

승수의 위치	결론
전부 단위원 안	점근 안정 (지수적으로 0으로)
하나라도 단위원 밖	불안정
전부 단위원 안 또는 위, 원 위의 것은 반단순	안정 (유계, 점근 아님)
단위원 위에 조르당 블록이 있음	불안정 (다항식 성장)

지수 언어로 옮기면 $\operatorname{Re}\lambda_i < 0$ 이 점근 안정이다. 상수계수 문제에서 고유값 실수부를 보던 것이 그대로 옮겨온 형태이고, 그래서 플로케 지수를 “유효 고유값”이라 부르기도 한다. 최대 실수부는 랴푸노프 지수의 주기계 판본에 해당한다.

또 하나 유용한 항등식. 승수들의 곱은 앞서의 리우빌 공식으로

\prod_{i=1}^{n}\mu_i = \det M = \exp\!\left(\int_0^{T}\operatorname{tr}A(s)\,ds\right)

이다. 해밀턴계처럼 $\operatorname{tr}A \equiv 0$ 이면 곱이 정확히 1이라, 승수가 $\mu$ 와 $1/\mu$ 로 짝지어 나타난다. 보존계에서는 점근 안정이 원리적으로 불가능하다는 사실이 이 한 줄에서 바로 읽힌다.

5. 마티외 방정식과 안정성 선도[편집]

플로케 이론의 교과서 예제이자 물리적 본체가 **마티외 방정식**이다.

\ddot x + \bigl(\delta + \varepsilon\cos t\bigr)\,x = 0

강성이 시간에 따라 주기적으로 흔들리는 진동자다. 마티외가 1868년 타원형 막의 진동을 다루며 얻었고, 힐이 1877년 달의 근지점 운동을 계산하며 더 일반적인 힐 방정식(임의 주기함수 계수)으로 확장했다. 이 계열의 이론적 뼈대가 그대로 플로케다.

$(\delta,\varepsilon)$ 평면에서 안정·불안정 영역을 칠한 것이 안정성 선도(stability chart)다. 그림의 특징이 뚜렷하다. $\varepsilon = 0$ 축의

\delta = \frac{n^2}{4},\qquad n = 1, 2, 3, \dots

에서 쐐기 모양 불안정 영역이 혀처럼 뻗어 나온다. $n=1$ 인 $\delta = 1/4$ 이 주공진(principal parametric resonance)이며, 물리적 의미가 명확하다. 고유진동수가 $\sqrt\delta = 1/2$ 인데 계수는 진동수 1로 흔들리니, 구동 진동수가 고유진동수의 정확히 두 배일 때 가장 위험하다는 것이다. 그네를 탈 때 무릎을 한 주기에 두 번 굽히는 것이 정확히 이 공진이다.

혀의 굵기도 규칙적이다. 첫 혀의 경계는 섭동 전개로

\delta = \frac14 \pm \frac{\varepsilon}{2} + O(\varepsilon^2)

이라 폭이 $O(\varepsilon)$ 이고, $n$ 번째 혀의 폭은 $O(\varepsilon^{n})$ 으로 급격히 가늘어진다. 고차 공진은 존재하지만 실질적으로 잡기 어렵다는 뜻이고, 실험에서 주공진만 잘 관측되는 이유다.

용어 하나만 정리하고 가자. 이 쐐기들을 아널드 혀(Arnold tongue)라 부르는 문헌이 많은데, 엄밀히 아널드 혀는 원사상의 위상 잠금 영역이나 나이마크-사커 분기의 공진 쐐기를 가리킨다. 마티외 선도의 것은 보통 불안정 혀(instability tongue) 또는 공진 혀라 부른다. 둘 다 매개변수 평면에서 정수비 공진을 꼭짓점으로 뻗는 쐐기라 이름이 섞여 쓰이지만, 논문에 쓸 때는 구분하는 편이 안전하다.

감쇠 $c\dot x$ 를 넣으면 혀가 $\varepsilon$ 축에서 위로 들려 올라간다 — 즉 구동이 어느 문턱을 넘어야만 불안정해진다. 실제 계에 감쇠가 항상 있으므로 이 문턱이 매개 여기 설계의 실무 기준이 된다.

6. 극한주기해에 붙일 때 — 승수 하나는 항상 1[편집]

비선형 자율계 $\dot{\mathbf x} = \mathbf f(\mathbf x)$ 가 주기 $T$ 의 주기해 $\mathbf x^{*}(t)$ 를 가진다 하자. 그 둘레의 작은 교란 $\boldsymbol\xi$ 는 변분방정식

\dot{\boldsymbol\xi} = D\mathbf f\bigl(\mathbf x^{*}(t)\bigr)\,\boldsymbol\xi

를 따르고, $\mathbf x^{*}$ 가 $T$ -주기이므로 계수행렬 $D\mathbf f(\mathbf x^{*}(t))$ 도 $T$ -주기다. 비선형 문제의 안정성 판정이 정확히 플로케의 무대로 들어온다.

여기서 자율계 특유의 사실이 하나 나온다. 원 방정식 $\dot{\mathbf x} = \mathbf f(\mathbf x)$ 를 시간에 대해 미분하면

\frac{d}{dt}\dot{\mathbf x}^{*} = D\mathbf f(\mathbf x^{*})\,\dot{\mathbf x}^{*}

즉 속도벡터 $\dot{\mathbf x}^{*}(t)$ 자신이 변분방정식의 $T$ -주기 해다. 주기해이므로 대응하는 승수가 정확히 1이다. 기하학적으로는 궤도를 따라 앞뒤로 밀어도 같은 궤도라는 사실 — 위상 방향으로는 복원력도 발산력도 없다는 뜻이다.

따라서 자율계 극한주기해의 안정성은 자명한 승수 1을 제외한 나머지 $n-1$ 개로 판정한다. 전부 단위원 안이면 궤도적 점근 안정. 이 $n-1$ 개는 푸앵카레 단면 위에서 정의한 복귀사상의 야코비안 고유값과 정확히 일치한다.

승수가 단위원을 어떻게 뚫고 나가느냐가 곧 다음 분기 이론의 종류를 정한다.

탈출 방식	분기
$\mu \to +1$	순환의 접힘(안장-마디), LPC
$\mu \to -1$	주기배가, PD
켤레쌍이 단위원 통과	나이마크-사커(토러스), NS

평면계( $n=2$ )에서는 자명한 승수 1 하나를 빼면 남는 게 하나뿐이라 판정이 명시적이다.

\mu_2 = \exp\!\left(\int_0^{T}\nabla\!\cdot\mathbf f\bigl(\mathbf x^{*}(t)\bigr)\,dt\right)

한 주기 동안의 발산 적분이 음수면 안정. 벤딕손-뒬락 판정법과 같은 양이 여기서 다시 나오는 것이 우연이 아니다.

주기 구동계(비자율)에서는 자명한 승수가 없다. 강제항의 위상이 시간을 고정해 버리므로 $n$ 개 승수 전부가 의미를 갖고, 스트로보스코픽 사상의 고유값이 그대로 승수다. 자율계 코드를 비자율계에 그대로 돌리며 “승수 1이 안 나온다”고 당황하는 것이 흔한 통과의례다.

7. 수치적으로 모노드로미를 구하는 법[편집]

정석은 원 방정식과 변분방정식을 함께 적분하는 것이다.

\dot{\mathbf x} = \mathbf f(\mathbf x),\qquad \dot\Phi = D\mathbf f(\mathbf x(t))\,\Phi,\qquad \Phi(0) = I

미지수가 $n + n^2$ 개로 늘고, 한 주기 뒤 $M = \Phi(T)$ 를 얻어 고유값 문제를 푼다. 자코비안 행렬은 자동 미분이나 해석적 미분으로 주는 것이 좋다 — 유한차분으로 때우면 승수의 유효숫자가 뭉텅 날아간다.

여기서 밟는 함정이 셋이다.

자명한 승수를 정확도 계기로 쓴다. 자율계라면 계산된 승수 중 하나가 1에서 얼마나 벗어났는지가 곧 그 계산 전체의 오차 지표다. $1 \pm 10^{-3}$ 이 나왔다면 나머지 승수도 그 정도로만 믿어야 한다. 적분 허용오차와 주기 $T$ 의 정밀도가 이 값에 그대로 반영된다.
강한 수축이 섞이면 작은 승수는 계산할 수 없다. 승수가 $1$ 과 $10^{-16}$ 처럼 벌어지면 $M$ 의 조건수가 그 비율만큼 커지고, 작은 승수는 반올림 잡음에 완전히 묻힌다. 다행히 실무에서 그 값이 정확히 얼마인지는 대개 중요하지 않다 — “아주 안정”이면 충분하다. 중요한 것은 그 잡음이 큰 승수를 오염시키지 않게 하는 것이다.
행렬 곱을 그대로 쌓지 마라. 스텝별 전이행렬을 곱해 $M = M_K\cdots M_1$ 을 명시적으로 만들면 위 문제가 최악으로 나타난다. 처방은 주기 슈어 분해(periodic Schur decomposition, Bojanczyk–Golub–Van Dooren 1992)로, 곱을 형성하지 않은 채 인자들에 QR 분해류의 직교변환을 번갈아 적용해 곱의 고유값을 직접 뽑는다. MatCont나 AUTO-07p 계열은 여기서 한 발 더 나가, 직교 배치법으로 이산화된 경계값 문제의 야코비 구조를 그대로 소거해 승수를 얻는다. 궤도를 다시 적분할 필요조차 없다는 것이 이 방식의 장점.

주기 구동 선형계라면 더 쉽다. 단위벡터 $n$ 개를 초기조건으로 한 주기씩 적분한 결과를 열로 쌓으면 그게 곧 $M$ 이다. 마티외 안정성 선도를 그리는 코드가 실제로 이렇게 짧다 — $(\delta,\varepsilon)$ 격자를 훑으며 2×2 행렬 하나씩 만들고 $\lvert\operatorname{tr}M\rvert > 2$ 인지 보면 끝이다.²

8. 어디서 쓰이나[편집]

매개 공진. 아래위로 흔들리는 지지점 위의 역진자가 안정화되는 카피차 진자, 이중 평균법에서 다루는 폴 이온 트랩의 $a$ - $q$ 안정 선도, 유체 표면의 패러데이 파. 전부 마티외/힐 방정식의 혀 구조가 실물로 나타난 것이다.
회전 기계. 헬리콥터 로터나 풍력 터빈 블레이드는 회전좌표계 때문에 고정좌표계에서 본 운동방정식의 계수가 주기적이 된다. 등방 로터라면 다중블레이드 좌표변환으로 계수를 상수화할 수 있지만, 전진 비행처럼 방위각 의존성이 남으면 플로케 전이행렬법이 사실상 유일한 안정성 판정 수단이다. 플러터·모드 해석의 상수계수 상식이 그대로는 통하지 않는 영역.
절삭 채터. 밀링은 날이 주기적으로 들어왔다 나가므로 계수가 주기적이고, 게다가 이전 회전의 흔적 때문에 시간 지연까지 붙는다. 반이산화(semi-discretization)로 유한차원 모노드로미 연산자를 만들어 승수를 보는 것이 안정성 로브 선도를 그리는 표준 방법이다.
주기궤도의 안정성. 궤도 역학의 헤일로 궤도, 화학 진동 반응의 극한 순환, 신경 발화의 주기해. 구조적 안정성이 요구하는 “모든 주기궤도가 쌍곡”이라는 조건도 승수 언어로 적혀 있다.
주기 구동 양자계. 시간 주기 해밀토니안의 준에너지가 플로케 지수의 양자판이고, 여기서 플로케 공학·플로케 위상물질이 나왔다. 공간 주기성에 같은 논증을 적용한 것이 블로흐 정리이며, 주기경계조건을 쓰는 밴드 계산은 전부 이 구조 위에 있다.³

정리하면 플로케 이론은 “주기적으로 흔들리는 것”을 다루는 최소 공통 언어다. 흔들리는 대상이 강성이든 회전 위치든 해밀토니안이든, 한 주기의 전이행렬 하나로 환원한 뒤 그 스펙트럼을 보는 절차는 바뀌지 않는다.

9. 관련 문서[편집]

이 반례는 제어공학 강의에서도 단골이다. “시변계는 매 순간 극점을 보면 된다”는 직관이 얼마나 위험한지를 보여주는 최단 경로이기 때문. 참고로 마르쿠스-야마베는 원래 이것보다 강한 예상(자코비안 고유값이 항상 좌반평면이면 대역 점근 안정)을 내놓았는데, 그 예상 자체는 3차원 이상에서 반례가 나와 무너졌다. ↩
2차원 마티외에서 $\lvert\operatorname{tr}M\rvert > 2$ 가 불안정 조건인 이유는 리우빌 공식 덕분이다. $\operatorname{tr}A = 0$ 이므로 $\det M = 1$ 이고, 그러면 두 승수가 $\mu, 1/\mu$ 이라 특성방정식이 $\mu^2 - (\operatorname{tr}M)\mu + 1 = 0$ 이 된다. 판별식 부호가 곧 안정성이다. 선도 그리는 코드가 열 줄이면 끝나는 이유. ↩
이름 순서가 헷갈리는데 시간 쪽(플로케 1883)이 공간 쪽(블로흐 1928)보다 먼저다. 그래서 광결정·포노닉 분야에서는 아예 “블로흐-플로케”라고 두 이름을 붙여 쓴다. 물리학과에서 블로흐로 배운 사람이 기계공학과 회전익 세미나에서 같은 정리를 플로케로 다시 만나 당황하는 일이 종종 있다. ↩