불안정 스프링의 초기 상태·민감도·역복원

위치 사진 한 장으로 미래를 알 수 있을까

언덕 꼭대기 근처의 구슬을 같은 위치에서 촬영했다고 생각해 봅시다. 하나는 꼭대기를 향해 움직이고 다른 하나는 바깥쪽으로 움직일 수 있습니다. 힘의 법칙은 현재 위치에서의 가속도를 정하지만, 사진에 없는 초기 속도를 대신 정해 주지 않습니다. 뉴턴의 이계 방정식을 시간에 따른 상태 전이로 읽으려면 위치와 속도를 함께 지정해야 합니다.

이 단원에서는 세 질문을 나눕니다. 완전한 상태를 알면 해가 하나인가, 가까운 상태들은 얼마나 벌어지는가, 나중 상태를 기록하면 과거를 얼마나 정확하게 복원할 수 있는가입니다. 해가 유일하면서 초기 오차가 크게 증폭될 수 있고, 역함수가 존재하면서 관측 반올림에 민감할 수도 있습니다.

불안정 스프링의 법칙과 단위

질량 m>0, 불안정 강성 k≥0에 대해 다음 법칙을 사용합니다.

$$m\ddot x=kx, \qquad z=(x, v), \qquad \dot z=\begin{pmatrix}v\\(k/m)x\end{pmatrix}.$$

힘은 원점에서 멀어지는 방향입니다. 보통 스프링의 복원력 −kx와 부호가 다릅니다. 위치는 미터, 속도는 미터/초, 질량은 킬로그램, k는 뉴턴/미터 단위입니다. 퍼텐셜의 원점을 꼭대기에 두면 U=−kx²/2입니다.

실제 언덕에서는 먼 위치까지 이차 퍼텐셜이 유지되지 않을 수 있습니다. 여기서는 표시 구간 전체에 이상화한 법칙을 적용합니다. 큰 변위의 지수 성장을 실제 장치의 무한한 에너지 공급으로 해석하지 않습니다.

완전한 상태와 투영된 교차

현재 벡터장은 선형이고 매끄럽습니다. m과 k를 고정하고 x(0), v(0)을 함께 주면 하나의 전역 해가 정해집니다. 위치만 주면 여러 초기 속도가 남습니다. 반대로 가속도는 이미 (k/m)x(0)으로 결정되므로 독립적인 추가 초기값이 아닙니다.

위치–시간 그래프 두 개가 같은 시각에 만났더라도 속도가 다르면 완전한 상태는 다릅니다. 위상 공간에서는 다른 점입니다. 이 매끄러운 자율계에서 같은 시각의 완전한 상태가 일치한 두 해가 이후 서로 다른 미래로 갈라질 수는 없습니다. 그래프의 투영과 상태 공간의 교차를 구별해야 합니다.

성장 모드와 감쇠 모드를 분해하기

k>0이면 λ=√(k/m)으로 둡니다. λ의 단위는 초의 역수입니다.

$$x(t)=c_+e^{\lambda t}+c_-e^{-\lambda t}, \qquad v(t)=\lambda c_+e^{\lambda t}-\lambda c_-e^{-\lambda t},$$ $$c_+=\frac{x_0+v_0/\lambda}{2}, \qquad c_-=\frac{x_0-v_0/\lambda}{2}.$$

두 계수의 단위는 길이입니다. v0=λx0는 성장 방향, v0=−λx0는 감쇠 방향입니다. 일반적인 초기 상태에는 두 성분이 함께 있습니다. 감쇠하는 방향 하나가 있다고 평형점이 안정한 것은 아닙니다. 임의로 작은 주변에도 성장 성분이 있는 초기 상태가 존재합니다.

수치 단위에서 m=k=1로 두면 (x0, v0)=(1, −1)은 x=e^(−t), (1, 1)은 x=e^t입니다. 두 초 뒤 위치는 약 0.1353m와 7.3891m입니다. 같은 위치에서 출발하지만 속도까지 지정하면 어느 미래인지 모호하지 않습니다.

에너지는 상태 전체를 대체하지 못한다

보존되는 역학적 에너지는 다음과 같습니다.

$$E=\frac{m v^2}{2}-\frac{k x^2}{2}.$$

양의 운동에너지와 음의 퍼텐셜이 상쇄될 수 있습니다. 앞의 두 초기 상태는 모두 E=0입니다. 따라서 위치와 에너지를 함께 알아도 속도의 부호를 알 수 없습니다. 에너지 준위는 상태 공간의 여러 점으로 이루어진 집합입니다.

에너지 그래프가 0에 머문다는 사실은 물체가 정지했거나 초기 오차에 둔감하다는 뜻이 아닙니다. 크게 성장한 상태에서는 큰 양수와 큰 음수를 부동소수점으로 빼므로 작은 에너지 잔차가 남을 수 있습니다. 초기 에너지는 제어 코드로부터 정확한 유리수로 계산하고, 이후 에너지와 절대·상대 규모 잔차는 수치 결과로 구분합니다.

단위를 맞춘 상태 오차와 상계

위치 오차 제곱과 속도 오차 제곱을 곧바로 더하면 단위가 다릅니다. 기준 시간 T=1초를 선언하고 y=(x, Tv)를 사용하면 두 성분이 모두 길이입니다.

$$\delta y(t)=M_T(t)\delta y(0), \qquad M_T(t)=\begin{pmatrix} \cosh(\lambda t)&\sinh(\lambda t)/(\lambda T)\\ T\lambda\sinh(\lambda t)&\cosh(\lambda t) \end{pmatrix}.$$

선형계이므로 두 해의 차이도 같은 행렬로 전파됩니다. 최대 특이값은 다음 유클리드 노름 상계를 줍니다.

$$\|\delta y(t)\|_2\leq\sigma_{\max}(M_T(t))\|\delta y(0)\|_2.$$

오른쪽은 모든 방향에 적용하는 상계입니다. 실제 증폭률이 모든 오차 방향에서 이 값과 같은 것은 아닙니다. λT=1이면 성장 방향의 길이는 e^(λt)배, 감쇠 방향의 길이는 e^(−λt)배입니다. T를 바꾸면 오차를 재는 기준과 유한 시간 특이값이 달라지지만 물리적 궤적은 바뀌지 않습니다.

면적 보존과 불안정성은 양립한다

전이 행렬의 해석적 행렬식은 1입니다. 한 방향이 늘어나는 만큼 다른 방향이 줄어들어 작은 위상 공간 면적이 보존됩니다. 두 특이값은 서로 역수이지만 모든 길이나 좌표 불확실성이 보존되는 것은 아닙니다.

이 예제는 한 방향의 지수적 분리를 보이지만, 해가 명시적인 비유계 선형계입니다. 지수 성장 하나만으로 혼돈을 입증하지 못합니다. 비선형 혼합이나 유계 영역의 복잡한 재귀 운동을 이 방정식에 덧붙여 해석하지 않습니다. 여기서 배우는 대상은 불안정 방향과 조건수입니다.

강성이 0인 자유 입자 극한

k=0에서는 λ로 나누는 모드 식을 사용하면 안 됩니다. 올바른 극한은 다음과 같습니다.

$$x(t)=x_0+v_0t, \qquad v(t)=v_0, \qquad M_T(t)=\begin{pmatrix}1&t/T\\0&1\end{pmatrix}.$$

속도 오차는 시간이 지날수록 위치 오차를 만듭니다. λ=0이라고 해서 상태 오차의 보편적 상계가 e^(λt)=1이 되지는 않습니다. 이 전단 행렬의 최대 특이값은 일반적으로 1보다 큽니다. λ=0에서 성장·감쇠 모드 계수를 억지로 유한한 수로 표시하지 않습니다.

역함수와 관측 정밀도는 별개다

해석적 역전이는 M_T(−t)입니다. 정확한 나중 상태를 넣으면 정확한 초기 상태가 복원됩니다. 여기서는 위치와 속도를 각각 표시 단위의 0.001 간격으로 반올림한 가상 관측도 별도로 비교합니다. 이는 관측 연산이지 원래 운동이 정보를 잃었다는 가정이 아닙니다.

관측 오차가 η이면 초기 복원 오차는 M_T(−t)η입니다. 이 행렬의 최대 특이값으로 오차 상계를 얻습니다. 현재 행렬식이 1인 이차원계에서는 최대 역특이값과 최대 순방향 특이값이 같습니다. T=1초일 때 반올림 칸 전체를 고려한 최악 상계는 σmax√2×0.0005m입니다.

감쇠 해의 5초 상태는 약 (0.00673795, −0.00673795)입니다. 이를 (0.007, −0.007)로 반올림하면 역전이 결과가 약 (1.03889, −1.03889)가 됩니다. 반올림 없는 수치 역전이의 잔차도 따로 확인합니다. 해석적 가역성, 부동소수점 오차, 관측 오차는 서로 다른 주장입니다.

AI 모델에 적용하기

위치만 입력받는 동역학 모델은 숨은 속도에 대한 조건부 평균을 예측할 수 있습니다. 그 평균이 하나의 올바른 뉴턴 궤적이라는 보장은 없습니다. 에너지를 정확히 보존하더라도 잘못된 속도 가지를 선택할 수 있습니다. 상태 전이 검증과 에너지 잔차 검증은 다른 실패를 찾습니다.

연구에서는 상태 변수, 관측 정밀도, 시간 구간, 단위와 노름을 함께 보고해야 합니다. 한 숫자의 증폭률만 제시하지 말고 방향 의존성도 확인하세요. 계산식 자체의 잔차뿐 아니라 별도 행렬 지수와 ODE 적분으로 결과를 대조해야 합니다.

해설이 있는 연습 문제

1. 한 위치에서 두 미래. m=k=1, x0=1에서 v0=−1과 v0=1을 비교하세요. 해설: 모드 계수는 각각 (0, 1)과 (1, 0)입니다. 2초 위치는 e^(−2)와 e²이고 속도도 각각 −e^(−2)와 e²입니다. 위치만 같으며 완전한 초기 상태는 다릅니다.

2. 에너지로 가지를 고를 수 있을까. 두 상태의 초기 에너지를 구하세요. 해설: 모두 E=(1−1)/2=0입니다. 위치와 에너지로 속도의 부호를 결정할 수 없습니다. 영 에너지 준위에는 성장·감쇠 방향이 모두 들어 있습니다.

3. 오차의 방향. m=k=T=1에서 초기 오차 (0.1, 0.1)과 (0.1, −0.1)의 3초 뒤 길이를 비교하세요. 해설: 초기 길이는 같지만 첫째는 e³배, 둘째는 e^(−3)배입니다. 두 경우 모두 최대 특이값의 상계를 만족하고 성장 방향만 그 상계에 도달합니다.

4. 자유 입자의 전단. m=2, k=0에서 A=(1, 0), B=(1, 1), t=2로 두세요. 해설: A의 위치는 1, B의 위치는 3이고 속도 차이는 1입니다. 상태 차이 길이는 √5입니다. 초기 길이 1에 최대 특이값 1+√2를 곱한 상계보다 작습니다.

5. 반올림 뒤 역복원. 감쇠 해의 5초 상태를 (0.007, −0.007)로 반올림하세요. 해설: 같은 감쇠 방향이므로 복원 진폭은 0.007e^5≈1.03889입니다. 두 성분을 합친 초기 오차는 약 0.0550m이며 역전이 오차 상계와 양립합니다. 해의 비유일성을 뜻하지 않습니다.

6. 검증 기준 선택. 학습 모델이 E=0과 x0=1을 맞추지만 감쇠 자료에 성장 궤적을 출력합니다. 해설: 에너지 검사는 두 가지를 모두 허용합니다. v0=−1을 포함한 완전한 초기 상태와 뉴턴 상태 전이를 검사해야 합니다. 초기 상태 불일치와 수치 잔차를 구분해 보고하세요.

참고 자료와 다음 범위

위스콘신 대학교 해의 유일성 강의노트는 Lipschitz 조건과 해의 분리 상계를 설명합니다. Sutherland의 ODE 교재는 완전한 상태 공간에서의 유일성을 다룹니다. 이 단원의 불안정 스프링과 반올림 실험은 직접 유도한 예제입니다. 비선형 혼돈, 확률 동역학, 제약 상태 복원과 모델 오차는 이후 심화에서 다룹니다.