강성·빠른 모드·느린 모드

천천히 보고 싶은데 빠르게 계산해야 하는 이유

일부 물리계에는 빠르게 사라지는 과도 현상과 오래 지속되는 느린 변화가 함께 있습니다. 관측자는 느린 변화만 보고 싶어도 명시적 적분기는 빠른 감쇠율 때문에 작은 step을 요구할 수 있습니다. 이것이 강성을 이해하는 출발점입니다. 단순히 그래프가 가파르다거나 행렬 원소가 크다는 설명만으로는 부족합니다.

두 독립 모드 s′=-s, q′=-kq를 생각합니다. 정확해는 s(t)=exp(-t), q(t)=(A+ε)exp(-kt)입니다. 초기 교란 ε는 오직 시작할 때 더하며 매 step 새 잡음을 넣지 않습니다. k가 크면 느린 관측 시간에 비해 빠른 시간척도 1/k가 매우 짧습니다.

빠른 모드의 증폭을 계산하기

실수 양의 k에 대해 Euler의 비증폭 조건은 0≤kh≤2입니다. kh=2의 multiplier -1은 감쇠하지 않고 진동하므로 실제 감쇠를 정확히 재현하는 것과 다릅니다. k=40, h=0.1이면 -3이 되어 작은 교란도 step마다 세 배 커집니다. Backward Euler의 multiplier는 1/5라 안정적이지만 빠른 정확해의 exp(-4)와는 크게 다릅니다.

0인 모드는 불안정성을 숨길 수 있다

A=ε=0이면 정확한 산술의 Euler는 모든 step에서 q=0을 유지합니다. 화면이 조용하다는 사실은 방법이 모든 교란에 안정하다는 증거가 아닙니다. ε를 아주 작게 바꾸면 같은 multiplier가 잠재된 민감도를 드러냅니다. 이 페이지는 초기 교란의 영향과 이론적 증폭을 함께 표시합니다.

파랑은 느린 수치 모드이고 주황은 빠른 수치 모드와 빠른 정확해 기준입니다. 큰 fast 값이 y축을 지배하면 표에서 slow 오차를 따로 확인하세요. 독립 좌표를 합한 하나의 물리량으로 오해하지 않는 것도 중요합니다.

손계산과 직접 풀기

k=40, T=1, N=10, A=1이면 explicit fast 값은 (-3)^{10}=59049입니다. 반면 정확값은 exp(-40)으로 거의 0입니다. Slow Euler는 0.9^{10}≈0.348678이며 정확한 exp(-1)과 비교할 수 있습니다. 큰 fast 오차와 상대적으로 작은 slow 오차는 서로 다른 현상입니다.

직접 풀기: k=80에서 Euler가 비증폭이려면 h는 얼마 이하여야 하나요? 풀이: h≤2/80=0.025입니다. 하지만 경계에서 감쇠가 없으므로 정확한 모드 감쇠를 원하면 더 작은 h가 필요합니다. 안정성 한계에 딱 맞추는 것과 오차 목표를 만족시키는 것은 다릅니다.

AI와 다중 시간척도

반응계, 확산 이산화, 연속 깊이 모델은 서로 다른 고유값 시간척도를 만들 수 있습니다. 이 대각 예제는 비정규 행렬의 과도 증폭이나 비선형 강성 전체를 설명하지 않습니다. 실제 연구에서는 Jacobian 구조, 비선형 풀이 비용, 관측 시간, 허용 오차와 계산 예산을 함께 기록해야 합니다. Implicit을 선택했다는 이름만으로 정확도와 효율이 보장되지는 않습니다.

참고: SciPy 적분기와 강성 방법 설명.