유한 집단의 질서와 동조

1. 같은 방향과 같은 속도는 다른 질문이다

사람들이 원형 트랙에서 걷는 모습을 생각하자. 모두 한곳에 모여 있어도 걷는 속도가 다르면 곧 흩어진다. 반대로 트랙 양쪽에 나뉘어 같은 속도로 걸으면 두 무리는 계속 떨어져 있다. 첫 장면의 순간적인 정렬과 둘째 장면의 공통 속도는 서로 다른 관측이다.

이번 단원에서는 개별 진동자의 진폭을 일정하게 두고 위상만 추적한다. 앞 단원의 Hopf 진폭 모형에서 진폭 회복이 충분히 빠르고 결합이 적절히 약할 때 위상 축약을 생각할 수 있지만, 여기서 사용하는 사인 결합이 모든 진동계에서 자동으로 유도되는 것은 아니다. 결합식과 적용 조건을 먼저 지정한다.

2. 전체 결합을 N으로 나누는 이유

진동자 i의 위상을 θᵢ, 고유 각주파수를 ωᵢ라 하자. 각주파수는 결합이 없을 때의 속도다. 결합을 켠 뒤 실제 속도는 일반적으로 ωᵢ와 다르다. 이번에는 자기 결합 항이0인 완전 연결 집단을 쓴다.

$$\dot\theta_i=\omega_i+\frac{K}{N}\sum_{j=1}^{N}\sin(\theta_j-\theta_i),\qquad K\geq0.$$

위상은 무차원 각도이고 시간 단위는 s다. ωᵢ와 K의 단위는 s⁻¹이다. N으로 나누면 구성원이 늘어도 한 사람이 받는 총 결합의 규모가 K 정도로 유지된다. 이를 빼면 N을 바꾸는 실험이 결합 강도까지 함께 바꾸게 된다. 실습 계약은 N=1…32, |ωᵢ|≤8, K≤16, 관측 시간0…4s다. 이는 계산 범위이며 일반 이론의 제한이 아니다.

3. 복소수 평균이 보여주는 그림

각 위상을 단위원 위 화살표로 바꾸고 평균한다. 평균 화살표의 길이는 정렬 정도, 방향은 평균 위상이다. 긴 화살표들이 서로 상쇄될 수 있다는 점이 핵심이다.

$$Z_m=\frac1N\sum_{j=1}^N e^{im\theta_j},\qquad R_m=|Z_m|,\qquad \Psi_m=\arg Z_m\quad(Z_m\ne0).$$

m=1에서는 원래 화살표를, m=2에서는 각도를 두 배 한 화살표를 평균한다. 따라서 반대 방향의 두 집단은 첫 평균에서는 상쇄되고 둘째 평균에서는 겹친다. 정확한 Z₁=0에는 평균 방향이 없다. 화면이 마지막 방향을 남기거나 작은 부동소수점 잔차의 방향을 유의미한 각도로 보여주면 오해를 만든다. 일반 수치의 작은 R₁을 정확한0이라고 판정하지도 않는다.

4. R₁ 하나로는 집단을 식별할 수 없다

다음 네 진동자의 고유 주파수는 모두1이고 K=3이다. 표의 위상은 초기값이며, 이후에는 모든 위상에 t가 더해진다. 표의0과1은 균형 구조에서 유도한 정확한 항등식이다.

초기 위상	R₁	R₂	다른 관측	모든 실제 속도
0, 0, 0, 0	1	1	R₄=1	1
0, 0, π, π	0	1	R₄=1	1
0, π/2, π, 3π/2	0	0	R₄=1	1

둘째와 셋째는 R₁만 보면 같지만 집단 구조는 다르다. 셋째에서는 R₂까지 보아도 규칙적인 사분점 배열을 놓친다. 원 위 점들과 높은 조화 변수를 함께 보아야 한다. 이 표가 안정성도 같다는 뜻은 아니다. 정확한 대칭 해가 존재하는 것과 작은 교란 뒤 그 해로 돌아오는 것은 별개의 문제다.

5. 순간 속도와 시간 창 평균

고유 주파수, 순간 속도, 유한 관측 창의 평균 속도를 구분한다. 위상을 감아 표시하면 한 바퀴를 돈 정보가 사라지므로, 평균 속도에는 연속적으로 누적한 위상을 사용한다.

$$\overline\Omega_i(T)=\frac{\theta_i(T)-\theta_i(0)}{T}\ (T>0),\qquad S_{\rm inst}(t)=\max_i\dot\theta_i(t)-\min_i\dot\theta_i(t).$$

T=0의 평균 속도는 정의되지 않으므로 빈값이다. 짧은 창에서 속도가 비슷해 보여도 장시간 동조를 증명하지 않는다. 예를 들어 초기 위상이 모두0이면 R₁=1이지만 고유 주파수가 다를 때 초기 결합 항은 모두0이라 즉시 서로 다른 속도로 움직인다. 반대로 표의 균형 집단은 R₁=0이어도 정확한 해에서 속도가 같다. 실제 실험에는 측정 오차와 관측 창도 명시해야 한다.

6. 전체 평균 속도가 보존되는 이유

한 쌍의 사인 항은 서로 반대다. 모든 방정식을 더하면 내부 결합이 상쇄된다. 이는 개별 속도가 고유 주파수와 같다는 주장이 아니다.

$$\frac1N\sum_i\dot\theta_i=\frac1N\sum_i\omega_i,\qquad \frac1N\sum_i\theta_i(t)=\frac1N\sum_i\theta_i(0)+t\frac1N\sum_i\omega_i.$$

구현은 각 쌍의 힘을 한 번 계산해 두 진동자에 더하고 뺀다. 그럼에도 부동소수점 합계에는 잔차가 남을 수 있다. 정확한 평균 주파수 항등식과 실제 계산의 평균 잔차를 따로 반환한다. 잔차를 숨기기 위해 매 단계 평균 위상을 다시 맞추지 않는다. 모든 ωᵢ에 같은 상수를 더하면 공통 회전만 바뀌고 상대 위상과 Rₘ은 유지되어야 한다.

7. 정확한 대칭과 계산의 대칭 오차

동일 주파수의 균형 m개 집단은 각 집단을 같은 수의 점으로 채운 정다각형이다. 복소수 합이0이므로 각 점이 받는 총 결합도0이다. 모든 위상은 같은 속도로 회전하고 상대 배열이 유지된다. m=1의 완전 정렬도 같은 방식으로 유지된다.

$$\theta_j(t)=\theta_j(0)+\omega_0t,\qquad \sum_{q=0}^{m-1}e^{2\pi i q/m}=0\ (m>1).$$

이 항등식은 입력 위상을 π의 정확한 유리수 배로 지정하고 정수 관계를 검사해 적용한다. 가까워 보이는 위상을 대칭 상태로 반올림하지 않는다. 작은 교란은 다른 초기 조건이다. K=0이면 주파수가 달라도 각 진동자의 자유 회전을 직접 계산할 수 있다.

별도로 수정하지 않은 RK4 궤적도 보관한다. ω=(0,0), θ=(0,π), K=16, T=4의 정확한 해는 정지한 반대 위상이다. 하지만 기계의 sin(π) 잔차가 불안정 방향으로 증폭되면 수치 궤적은 크게 벗어날 수 있다. 이 비교는 물리적 잡음을 넣은 실험이 아니다. 실제 잡음을 논하려면 확률 모형과 크기를 따로 지정해야 한다.

8. 세 격자의 차이는 오차의 증명이 아니다

일반 초기 조건은 고전적 RK4로 n, 2n, 4n 간격을 각각 적분하고 가장 촘촘한 결과를 주 궤적으로 사용한다. 같은 물리 시각의 감지 않은 위상을 비교한다.

$$d_1=\max_{i,k}|\theta_i^{(n)}(t_k)-\theta_i^{(2n)}(t_k)|,\quad d_2=\max_{i,k}|\theta_i^{(2n)}(t_k)-\theta_i^{(4n)}(t_k)|,\quad E_{\rm formal}=d_2/15.$$

충분히 작은 간격에서 선도 오차가 h⁴에 비례하면 두 차이의 비가 약16이 되고 가장 촘촘한 오차를 d₂/15로 추정할 수 있다. 그러나 해당 근사 영역에 들어갔는지는 별도 문제다. 독립 DOP853 참조와 비교한 의도적인 거친 격자 예제에서는 실제 위상 오차가1rad보다 크고 형식적 추정의10배도 넘는다. 추정 수치가 나온다는 사실만으로 결과를 신뢰하지 않는다.

실습 기록 수 S는 계산 간격 수와 다르다. 표시 기록 i는 기본 격자의 floor(in/S)번째 시각을 선택한다. 나누어떨어지지 않으면 표시 간격은 일정하지 않을 수 있으므로 실제 시간을 표시한다. 기록 수를 바꿔도 끝점과 전체 격자 진단은 같아야 한다. 독립 참조 역시 고정밀 허용오차를 사용한 부동소수점 수치해이며 전 입력의 엄밀한 오차 구간은 아니다.

9. 같은 주파수일 때 내려가는 함수

모든 ωᵢ가 같을 때 그 공통 회전을 제거한 좌표 φᵢ를 쓰면 결합은 다음 함수의 기울기 흐름이다.

$$V=-\frac{K}{2N}\sum_{i,j}\cos(\phi_j-\phi_i)=-\frac{KN}{2}R_1^2,\qquad \dot\phi_i=-\frac{\partial V}{\partial\phi_i},\qquad \dot V=-\sum_i\dot\phi_i^2\leq0.$$

미분해 사인 항의 부호와 N 계수를 확인해보자. V가 감소한다는 사실은 모든 초기 조건이 완전 정렬로 간다는 뜻이 아니다. 앞의 정확한 균형 배열도 결합 속도가0인 정지 상대 배열이다. 주파수가 서로 다르면 위의 마지막 등식을 그대로 쓸 수 없다. 이 위상 모형의 V는 동역학적 퍼텐셜이며 지정하지 않은 물리적 열·에너지·엔트로피로 해석하지 않는다.

10. 두 진동자 모형과의 연결

N=2로 줄여 앞 단원과 계수를 비교한다. 이때 각 간선의 계수는 K/2다.

$$\delta=\theta_1-\theta_2,\quad \Delta=\omega_1-\omega_2,\qquad \dot\delta=\Delta-K\sin\delta,\qquad R_1^2=\frac{1+\cos\delta}{2}.$$

K>|Δ|의 안정·불안정 잠김, K=|Δ|의 임계 병목, K<|Δ|의 미끄러짐을 이전의 해석해와 대조한다. R₁은 상대 위상의 코사인만 알려주므로 δ의 부호나 몇 바퀴 돌았는지는 복원하지 못한다. 유한 집단 구현이 두 진동자의 K/2를 놓치면 잠김 임계와 회복률이 모두 달라지므로 이 비교는 중요한 정규화 검사다.

11. Lorentzian 무한 집단은 별도의 비교 모형이다

이제 유한 배열과 다른 가정을 도입한다. 무한 집단의 고유 주파수 밀도를 Lorentzian으로 지정하고, 중심 Ωc와 폭 γ>0을 준다. 이 분포의 중심은 산술 평균이 아니다. Cauchy 분포의 평균은 정의되지 않는다. 앞에서 계산한 유한 배열의 경험적 평균과 혼동하지 않는다.

$$g(\omega)=\frac{\gamma}{\pi[(\omega-\Omega_c)^2+\gamma^2]},\qquad \dot R=\left(\frac K2-\gamma\right)R-\frac K2R^3.$$

이 축약식에는 무한 집단과 특별한 불변 다양체 및 해석성 조건이 필요하다. Ott–Antonsen 원 논문의 II–III절은 위상 밀도의 Fourier 계수가 거듭제곱으로 묶이는 형태와 복소 주파수 평면에서의 조건을 명시한다. 임의의 작은 유한 집단이 이 식을 정확히 따른다는 뜻은 아니다. 실습의 γ는 따로 지정한 비교 모형의 값이며 유한 주파수 배열에 맞춘 추정치가 아니다.

12. 임계값과 회복시간을 직접 유도하기

위의 1차원 방정식에서 μ=K/2−γ라 쓰고 양의 고정점을 구한다. 비영 가지가 존재하는 쪽에서 그 주변을 선형화한다.

$$K_c=2\gamma,\quad R_*^2=1-\frac{2\gamma}{K}\ (K>2\gamma),\quad f'(R_*)=-(K-2\gamma),\quad \tau_+=\frac1{K-2\gamma}.$$

γ=1인 비교표다. τ는 명시한 가지에서의 작은 진폭 교란에 대한 지수 이완시간이다.

K	μ	양의 고정점 R*²	해당 지수 이완시간
1	−1/2	없음	원점에서 2
2	0	없음	임계에서 정의하지 않음
3	1/2	1/3	양의 가지에서 1

원점 아래쪽의 감쇠율은 γ−K/2이고 양의 가지의 감쇠율은 K−2γ다. 두 값을 혼용하면2배 오류가 난다. K=2γ에서는 선형률이0이고 Ṙ=−γR³이므로 양의 초기값에서 대수적으로 감소한다. 지수 이완시간 칸을0이나 임의의 큰 수로 채우지 않는다. K와2γ의 비교는 정확한 유리수로 하며, 두 숫자가 기계 표시에서 같아 보여도 부호를 바꾸지 않는다.

13. AI에서 무엇을 관측하고 추론하는가

많은 개체의 상태를 하나의 통계량으로 줄이는 일은 학습된 표현과 닮아 있다. R₁만 입력받는 모델은 표의 두 균형 배열을 구분할 수 없다. 이는 데이터가 부족한 것이 아니라 관측 사상이 정보를 버린 경우다. R₂와 R₄, 점 배열, 속도와 시간 창을 추가해 어떤 가설이 구별되는지 확인하자.

Lorentzian 식의 집단 질서 R과 Hopf 정상형의 개별 진동자 진폭 r은 같은 관측이 아니다. 두 방정식에 삼차 포화 항이 있어도 적용 대상을 바꾸어 쓸 수 없다. 중심 회전이0인 질서 매개변수의 분기를 통상적인 Hopf 분기라고 자동으로 부르지도 않는다. 그래프·관성·지연·잡음·일반 결합과 유한 크기 효과, 다른 주파수 분포, 안정성과 극한 순서의 증명은 연구 후속 범위다. 이 작은 모형의 검사 통과로 그 범위를 완료했다고 세지 않는다.

14. 풀이와 함께 확인하기

1. 반대 집단의 숨은 질서. 위상0,0,π,π의 R₁과 R₂를 구하라. 풀이: 첫 화살표 합은1+1−1−1=0이고 각도를 두 배 하면 모두1이다. 따라서 R₁=0,R₂=1이다. 평균 첫 위상은 정의되지 않는다.

2. 정렬과 속도. θ=(0,0), ω=(1,3)일 때 처음의 R₁과 실제 속도는? 풀이: R₁=1이지만 모든 사인 항이0이어서 속도는1과3이다. 순간 정렬에서 주파수 동조를 결론낼 수 없다.

3. 평균 보존과 쌍 환원. ω=(3,0), N=2, K=5의 평균 속도와 상대 위상식을 구하라. 풀이: 평균은3/2이고 상대식은 δ̇=3−5sinδ다. 두 간선에 각각5를 넣으면 상대식의 계수가10이 되어 다른 모형이 된다.

4. 격자 진단. 세 격자의 차이가 d₁=0.016,d₂=0.001rad이면 무엇을 말할 수 있는가? 풀이: 비는16이고 형식적 추정은1/15000rad다. 4차 오차 가정에 부합하는 징후이지만 보증은 아니다. 독립 참조나 더 작은 간격으로 확인한다. 표시 기록만 늘리는 것은 이 확인을 대신하지 못한다.

5. 분기 양쪽의 시간 척도. γ=1,K=3의 R와 회복시간을 구하고 K=1과 비교하라. 풀이: R=1/√3, 양의 가지 τ=1이다. K=1에서는 원점의 감쇠율1/2이므로 τ=2다. 같은 공식을 양쪽에 복사하면 틀린다.

6. 유한 배열과 연속 비교의 차이. γ=1,K=3이면 모든 유한 집단이 R₁=1/√3로 가는가? 풀이: 아니다. 동일 주파수의 정확한 사분점 배열은 R₁=0,R₄=1을 유지한다. Lorentzian 무한 집단의 양의 가지는 별도 가정 아래의 결과다. 부동소수점 궤적이 그 배열에서 벗어나더라도 그것만으로 물리적 잡음이나 유한 집단 임계값을 측정한 것은 아니다.