미분방정식과 시간 이산화

변화 법칙·초기 조건·수치해를 나누기

미분방정식은 순간 변화율을, 초기 조건은 궤적을 정합니다. 수치해법은 그 궤적을 이산 시간에서 근사합니다. 세 대상은 다릅니다. Neural ODE나 미분 가능한 시뮬레이터에서도 연속 법칙을 학습한다고 해법 오차가 사라지는 것은 아닙니다.

$x'(t)=-kx(t)$, $k>0$, $x(0)=x_0$를 씁니다. 시간 단위를 정하면 $k$는 그 역수 단위입니다. 상태는 무차원일 수 있지만 $k\Delta t$는 반드시 무차원입니다.

연속 문제 풀기

미분하면 방정식, 0을 대입하면 초기 조건을 확인합니다. 우변이 상태에 대해 전역 Lipschitz이므로 이 초기값 문제의 해는 유일합니다. 양의 초기 상태는 양수를 유지하며 0으로 감쇠합니다.

이것은 일차 완화 모형이며 마찰 없는 Newton 진동자가 아닙니다. 이차함수의 gradient flow도 수학적 대응 아래 같은 식을 쓰지만 물리 해석에는 단위·소산·계의 정의가 추가로 필요합니다.

전진 Euler 만들기

현재 미분으로 짧은 시간의 변화를 예측합니다.

$x_0=1,k=1,\Delta t=0.25$이면 첫 수치값은 0.75이고 정확한 값은 약 0.778801입니다. 20회 후 시간은 5이며 정확한 값은 약 0.006738, Euler는 $0.75^{20}\approx0.003171$입니다.

기울기를 이전 상태에서 구합니다. 새 미지의 상태로 바꾸면 후진 Euler라는 다른 방법이며 안정성도 달라집니다.

안정 조건 유도

한 단계 배율은 $r=1-k\Delta t$입니다. $|r|<1$, 즉 $0<k\Delta t<2$일 때 수치값이 0으로 갑니다. $0<k\Delta t\le1$이면 음이 아닌 상태를 유지합니다. 1과 2 사이에서는 부호를 바꾸며 작아져 정확한 양의 해와 다른 거동을 보입니다.

$k\Delta t=2$에서는 배율 $-1$로 크기가 줄지 않습니다. 그보다 크면 수치 크기가 증가합니다. 불안정 예제는 $k=3,\Delta t=1$이어서 매번 $-2$를 곱합니다. 참 감쇠 법칙이 아니라 수치 불안정입니다.

정확도는 같은 시간에서 비교하기

매끄러운 해의 전진 Euler는 적절한 안정성 아래 한 단계 국소 절단 오차가 $\Delta t^2$, 고정 유한 구간의 전역 오차가 $\Delta t$ 차수입니다. 차수는 세분화 과정의 성질이지 임의의 두 실험 비교에 대한 보장은 아닙니다.

실험은 20회 고정입니다. 간격을 줄이면 최종 시간도 달라집니다. 다른 시각에서 오차가 작아졌다는 사실만으로 수렴 차수를 검증하지 마세요. 연구 검산에서는 $T$를 고정하고 $h,h/2,h/4$의 보폭에 맞춰 단계 수를 늘려 같은 $T$에서 비교합니다.

그림과 실패 원인 읽기

파랑은 정확한 해, 주황은 같은 시점의 Euler 값입니다. 안정적으로 보이게 값을 잘라내지 않습니다. 수치해가 커지면 자동 축 때문에 참해가 평평해 보일 수 있으므로 최종 오차와 무차원 배율을 함께 읽으세요.

연속 법칙이 실제 계를 설명하는지, 수치해법이 그 법칙을 잘 푸는지, 추정한 매개변수가 관측에 맞는지는 다른 질문입니다. 거친 해법의 오차를 잘못된 연속 매개변수로 보상하여 적합도가 좋아질 수도 있습니다.

AI와 연구 연결

학습 모형 안의 해법도 계산 목적함수의 일부가 됩니다. 간격·허용 오차·보간·기울기 계산이 결과에 영향을 줄 수 있습니다. 한 해상도의 손실만 보지 말고 해법 세분화에 따른 매개변수와 예측도 비교하세요. 이산 해법의 기울기인지 연속 민감도를 따로 유도한 것인지 기록합니다.

이 단원은 선형 ODE 하나를 다룹니다. 강직계·PDE·경계 조건·적응 해법에는 별도 분석이 필요합니다. 작은 예제 통과가 범용 과학 시뮬레이터의 검증은 아닙니다.

연습과 근거

1. $k=2,h=0.25$의 배율은? 답: 0.5입니다.
2. $k=1,h=1.5$에서 양의 초기값이 수치적으로 부호를 바꾸는 이유는?
3. 최종 시간을 고정한 세분화 검사를 설계하고 일차 수렴의 오차 비율을 예측하세요.
4. 거친 해법이 추정한 물리 상수에 편향을 만들 수 있는 이유는?

배경은 MIT Differential Equations를 참고하세요. 위 정확해와 안정성은 지정한 시험 방정식에서 직접 계산합니다.