진폭의 탄생과 임계 이완

1. 회전하고 있다는 것과 진동을 유지한다는 것

복소수 z=x+iy를 평면 위의 화살표로 생각하자. 화살표의 길이는 진폭이고 방향은 위상이다. 화살표가 회전하면서 짧아진다면 진동은 보이지만 결국 사라진다. 반대로 길이가 일정한 양수에 가까워지면서 계속 회전하면 지속적인 진동을 얻는다. 유한한 관측 창에서 몇 번 원을 그렸다는 사실만으로 두 상황을 구별할 수 없다.

이 단원은 진폭의 성장·포화와 위상의 회전을 분리한다. 특히 “임계에서 느려진다”가 회전 자체의 정지를 뜻하는지, 진폭 교란의 회복이 느려진다는 뜻인지 관측량을 명시한다. 뒤의 위상 동기화 단원에 필요한 진폭 안정화 조건도 여기서 준비한다.

2. 모형과 단위부터 고정하기

사용할 결정론적 모형은 다음과 같다. 진폭은 무차원, 시간은 s이며 μ,ω,g,c의 단위는 s⁻¹이다. g는 양수로 제한한다.

$$\dot z=(\mu+i\omega)z-(g+ic)|z|^2z,\qquad g>0.$$

μ는 작은 진폭의 성장률, g는 성장을 억제하는 계수, ω는 원점 근처의 각주파수다. c는 진폭에 따라 회전 속도가 달라지는 정도다. c=0이어도 진폭의 성장과 포화는 남는다. 이 다항식 자체를 실험 모형으로 사용한다. 실제 물리계가 이 식으로 근사되는 영역과 오차는 별도의 정당화가 필요하다.

$$\dot x=(\mu-g\rho)x-(\omega-c\rho)y,\qquad \dot y=(\omega-c\rho)x+(\mu-g\rho)y.$$

3. 극좌표로 읽되 원점을 잊지 않기

r>0에서 z=r exp(iθ)를 대입해 실수부와 허수부를 나누면 다음을 얻는다.

$$\dot r=\mu r-gr^3,\qquad \dot\theta=\omega-cr^2.$$

반지름 방정식은 θ를 몰라도 풀 수 있다. 그러나 위상 방정식에는 반지름이 들어간다. 이것이 진폭과 위상을 완전히 독립적인 것으로 취급할 수 없는 이유다. r=0에서는 방향이 없으므로 θ와 실제 위상 속도를 미제공으로 남긴다. 복소수 방정식에 ω가 있어도 영벡터가 관측 가능한 방향을 갖게 되는 것은 아니다.

4. 제곱 진폭은 Riccati 방정식이다

ρ=r²로 두면 제곱근 없이 진폭의 동역학을 계산할 수 있다. ρ는 확률이나 에너지가 아니라 여기서 선택한 무차원 관측량이다.

$$\dot\rho=2\mu\rho-2g\rho^2,\qquad u=\rho^{-1}\Longrightarrow\dot u=-2\mu u+2g.$$

두 번째 식의 변수는 u다. 역수 변환은 ρ>0에서만 사용한다. 비선형 방정식이 u의 선형 비동차 방정식으로 바뀌므로 닫힌 해와 독립 행렬 지수 검산을 모두 만들 수 있다. 원점의 해는 이 변환에 억지로 넣지 않고 따로 처리한다.

5. 임계값을 통과하는 정확한 해

초기 제곱 진폭을 ρ₀라 하면 다음 해를 얻는다.

$$\rho(t)= \begin{cases} \dfrac{\rho_0e^{2\mu t}}{1+(g\rho_0/\mu)(e^{2\mu t}-1)},&\mu\ne0,\\[6pt] \dfrac{\rho_0}{1+2g\rho_0t},&\mu=0. \end{cases}$$

이 식은 닫힌 모형의 해이며 시간 간격을 작게 잡아 얻은 근삿값이 아니다. 구현에서는 양의 큰 지수를 직접 만들지 않도록 μ의 부호에 따라 식을 다시 쓴다. μ=0은 정확한 정수·유리수 조건으로 분기한다. μ=10⁻⁹와 μ=0이 짧은 시간 동안 비슷하게 보일 수 있어도 같은 모형으로 분류하지 않는다.

6. 위상은 진폭의 누적 이력을 기억한다

A(t)를 제곱 진폭의 시간 적분으로 정의하자. D는 앞 절 해의 분모이며 μ=0에서는 1+2gρ₀t다.

$$A(t)=\int_0^t\rho(s)\,ds=\frac{\ln D(t)}{2g},\qquad \theta(t)-\theta(0)=\omega t-cA(t).$$

현재 반지름이 같아도 앞서 다른 진폭 이력을 거쳤다면 누적 위상은 다를 수 있다. 이것은 c≠0일 때의 기억 효과다. 계산하는 위상 변화는 감지 않은 누적 각도다. 이를 주기 안으로 접은 그림 좌표와 혼동하지 않는다. x,y는 초기 벡터를 이 각도만큼 회전시키고 r/r₀만큼 늘려 재구성하며, 매 단계 임의로 원 위에 투영하지 않는다.

7. 임계 아래에서는 무엇이 얼마나 빨리 줄어드나

μ<0이면 원점이 안정하다. 충분히 작은 진폭에서 r은 exp(μt), ρ는 exp(2μt)에 비례해 줄어든다. 따라서 같은 궤적이라도 무엇을 관측하는지에 따라 선형 이완시간이 다르다.

관측량·상태	선형 성장률	양의 이완시간
μ<0의 작은 진폭 r	μ	1/|μ|
μ<0의 작은 제곱 진폭 ρ	2μ	1/(2|μ|)
μ>0 원 주변의 반지름 교란	−2μ	1/(2μ)

이 표는 국소 선형화의 시간척도다. 유한 진폭이 정확히 반으로 줄어드는 시간을 보편적으로 이 값과 같다고 하지 않는다. μ>0의 원점은 불안정하므로 그 점에 “양의 회복시간”을 붙이지 않는다.

8. 임계에서는 지수 대신 대수적 감쇠

μ=0이고 ρ₀>0이면 다음의 서로 다른 반감시간을 얻는다.

$$r(t)=\frac{r_0}{\sqrt{1+2g\rho_0t}},\qquad t_{\rho/2}=\frac1{2g\rho_0},\qquad t_{r/2}=\frac3{2g\rho_0}.$$

선형 고윳값의 실수부가0이어도 비선형 항이 진폭을 줄인다. 원점은 약하게 끌어당기며 장시간 진폭은 t⁻¹ᐟ²로 줄어든다. 선형 이완시간이 발산한다는 말이 모든 유한 교란의 반감시간이 무한이라는 뜻은 아니다. 교란 크기와 μ를 어떤 순서로0에 보내는지도 구별해야 한다.

9. 임계 위의 원과 그 안정성

μ>0이면 양의 일정 진폭을 갖는 불변 원이 존재한다. 반지름 방향과 원을 따라가는 방향은 다른 안정성을 갖는다.

$$\rho_* =\frac\mu g,\qquad \Omega_* =\omega-c\frac\mu g,\qquad \lambda_{\rm radial}=-2\mu.$$

Ω*≠0일 때 이 원은 주기 궤도이며 주기는 2π/|Ω*|다. 전체 시간 원점이 자유로우므로 위상 방향은 중립적이다. 반지름이 안정하다는 사실이 모든 초기 위상이 같은 한 점으로 수렴한다는 뜻은 아니다. Ω*=0이면 원 위의 모든 점이 평형점이며 회전 주기를 정의하지 않는다. 음의 Ω*는 회전 방향의 반전이지 음의 주기가 아니다.

10. 실제 Hopf 조건과 중요한 예외

원점에서의 선형화는 다음과 같다.

$$J_0=\begin{pmatrix}\mu&-\omega\\\omega&\mu\end{pmatrix},\qquad \lambda_\pm=\mu\pm i\omega.$$

통상적인 Hopf 분기를 말하려면 μ=0에서 ω≠0이어야 한다. 이 계수 경로에서는 실수부가 μ와 함께 횡단하고 g>0의 포화가 초임계 형태를 준다. ω=0이면 선형 고윳값 두 개가0이므로 같은 이름으로 분류하지 않는다. Hopf의 정의·비퇴화 조건은 Kuznetsov의 원문 설명을 참고한다. 이 장의 해·예제·관측량 비교는 주어진 모형에서 직접 유도했다.

또 하나의 예외는 정확한 초기값 z₀=0이다. μ>0이어도 잡음 없는 미분방정식의 해는 계속0이다. 불안정성은 작은 교란의 성장을 뜻한다. 모형에 없는 잡음을 넣어 진동을 시작한 것처럼 표시하지 않는다.

11. 임계 이완과 위상 축약의 적용 조건

μ가0에 가까워지면 반지름 회복이 느려진다. 양의 μ에서도 회복률2μ가 외부 결합·구동의 변화율보다 충분히 빠르지 않으면 r을 즉시 일정하다고 둘 근거가 약해진다. 뒤의 위상 모형은 이 분리를 가정한 별도의 모형으로 읽어야 한다. 진폭이 사라지는 경계에서 위상만의 표현이 자동으로 유효하다고 주장하지 않는다.

이 단원에서의 임계 이완은 유한 차원 정상형의 동역학이다. 다수 진동자의 열역학적 극한이나 잡음에 의한 집단 전이와 같은 말로 합치지 않는다. 회전 주파수가 유한한 채로 진폭 회복시간만 길어질 수도 있다.

12. 손으로 확인하는 임계 예제

μ=0,g=1,c=1,ω=2,z₀=1을 선택하면 ρ=1/(1+2t), 누적 위상은 2t−ln(1+2t)/2다. 진폭은 줄지만 회전 속도는 2−ρ로 올라간다. 나선의 모양만 보고 진폭 회복이 빨라졌다고 해석하지 말자.

t [s]	ρ	r	위상 변화 [rad]
0	1	1	0
1/2	1/2	1/√2	1−ln2/2
1	1/3	1/√3	2−ln3/2
3/2	1/4	1/2	3−ln2

초기 진폭을 바꾼 실험의 민감도도 따로 구할 수 있다.

$$\frac{\partial\rho(t)}{\partial\rho_0} =\left(\frac{\rho(t)}{\rho_0}\right)^2e^{-2\mu t},\qquad\rho_0>0.$$

이 값은 선택한 한 초기 상태 근처의 민감도다. 모형 계수 전체의 식별 가능성이나 학습 알고리즘의 안정성을 한 번에 판정하는 값이 아니다.

13. 계산하고 설명하는 여섯 문제

1. 임계 예제에서 진폭과 제곱 진폭의 반감시간을 구하라. 풀이: ρ가1/2이 되는 시각은1/2s, r이1/2이 되는 시각은3/2s다. r의 절반은 ρ의1/4에 해당한다.

2. μ=−1,g=1,z₀=1인 경우의 ρ(t)와 작은 진폭 이완시간은? 풀이: ρ=1/(2e²ᵗ−1)이며 장시간 r은 상수×e⁻ᵗ다. r의 국소 선형 시간은1s이고 ρ의 시간은1/2s다.

3. μ=1,g=1,c=2,ω=2일 때 끌어당기는 원은 주기 궤도인가? 풀이: ρ*=1,Ω*=0이다. 원 위 점들은 평형점이며 회전 주기를 제공하지 않는다. 반지름 교란은 여전히2/s의 국소 속도로 줄어든다.

4. 같은 조건에서 c=3이면? 풀이: Ω*=−1rad/s이므로 반대 방향으로 회전한다. 주기는2πs이며 음수가 아니다. 반지름 안정성은 c 때문에 바뀌지 않는다.

5. μ>0인데 정확히 z₀=0에서 시작하면 왜 원으로 나가지 않는가? 풀이: 우변이0이므로 z(t)=0이 해다. 미분방정식의 유일성과 불안정성은 양립한다. 작은 양의 진폭을 별도로 주면 성장하지만, 그것은 다른 초기값이다.

6. 한 불변 원 궤적만으로 μ,g,ω,c를 모두 식별할 수 있는가? 풀이: 아니다. 관측은 μ/g와 ω−cμ/g의 조합만 제한한다. 다른 초기 진폭의 회복과 위상 속도 변화를 측정하는 설계가 추가 정보를 준다. 짧은 관측에서 오차가 작다는 사실만으로 임계값의 정확한 부호도 보장하지 못한다.

14. 수치 근거와 다음 연구로의 연결

생산 계산은 정확 유리수의 부호·동일성, 방향성256bit 지수·로그·제곱근·삼각함수 구간을 사용한다. 큰 누적 각도도 접거나 상태 벡터를 재정규화하지 않고 계산한다. x²+y²−ρ와 로그 진폭 수지를 별도 구간으로 확인하며, 대수적으로 성립하는 항등식과 반올림 잔차를 구별한다. 시간 기록 수는 물리적 해를 바꾸지 않는다.

독립 참조는 제품 모듈을 import하지 않고 역제곱 진폭의2×2행렬 지수를220자리에서 계산한다. 유한 사례 검산이 모든 입력에 대한 기계 증명은 아니다. 아직 읽기 화면·Worker·시각화 등록을 하지 않았으며 두 진동자와 다수 진동자의 가족도 남아 있다.

AI 연구에서는 학습된 잠재 동역학이 회전하는지, 진폭을 유지하는지, 교란 뒤 회복하는지 관측을 분리할 수 있다. 여기의 작은 모형을 임의의 신경망·실측계에 그대로 적용하지 않는다. 다음에는 반지름이 안정화된 위상 모형의 상대 위상과, 집단 질서 관측량의 정보 손실을 다룬다.