볼테라 적분방정식

1. 개요

볼테라 적분방정식 Volterra integral equation
제2종	$u(t)=f(t)+\int_0^t K(t,s)\,u(s)\,ds$
제1종	$\int_0^t K(t,s)\,u(s)\,ds=f(t)$
특징	적분 상한이 $t$ — 과거만 본다
연산자	준멱영, 스펙트럼 반지름 0
이산화	하삼각 연립방정식 → 시간 행진
이름	비토 볼테라 (1896)

적분 상한 하나 바꿨을 뿐인데 문제의 성격이 통째로 바뀐다.

볼테라 적분방정식(Volterra integral equation)은 미지함수가 적분 기호 안에 들어 있고, 그 적분의 상한이 독립변수 $t$ 자신인 적분방정식이다. 제2종은

$$u(t)=f(t)+\int_0^t K(t,s)\,u(s)\,ds,$$

제1종은 좌변의 $u(t)$ 가 없는 $\int_0^t K(t,s)u(s)\,ds=f(t)$ 다. 적분 구간이 고정된 $[a,b]$ 인 것은 프레드홀름 적분방정식이라 부른다.

상한이 $t$ 라는 건 현재 값이 과거 값의 가중합으로만 정해진다는 뜻이다. 그래서 볼테라 방정식은 “기억이 있는 상미분방정식”처럼 행동한다. 실제로 $u'=g(t,u)$, $u(0)=u_0$ 를 적분하면 $u(t)=u_0+\int_0^t g(s,u(s))\,ds$ 라는 (비선형) 볼테라 방정식이 되고, 피카르 반복이 바로 이 꼴에 대한 고정점 반복이다. 이탈리아 수학자 비토 볼테라가 1896년 일련의 논문에서 이 방정식의 역을 다뤘다.

이 문서는 왜 제2종 볼테라 방정식이 항상 잘 풀리는지, 제1종은 왜 까다로운지, 그리고 시간 행진형 수치해법과 그 $O(N^2)$ 비용 문제를 다룬다.

2. 제2종은 왜 항상 풀리나 — 스펙트럼 반지름 0

연산자 $(\mathcal K u)(t)=\int_0^t K(t,s)u(s)\,ds$ 의 거듭제곱을 보자. $\lvert K\rvert\le M$ 이면 귀납법으로

$$\big\lvert(\mathcal K^n u)(t)\big\rvert\le \frac{(Mt)^n}{n!}\,\lVert u\rVert_\infty$$

가 나온다. 적분 구간이 매번 $[0,t]$ 로 잘리면서 단체(simplex)의 부피 $t^n/n!$ 가 생기기 때문이다. $n!$ 은 어떤 지수보다 빨리 크므로 $\lVert\mathcal K^n\rVert^{1/n}\to0$, 즉 스펙트럼 반지름이 0 이다(이런 연산자를 준멱영이라 한다). 따라서 노이만 급수

$$u=\sum_{n=0}^{\infty}\mathcal K^n f$$

가 커널의 크기와 구간 길이에 상관없이 항상 수렴한다. 프레드홀름 방정식이 $\lambda$ 값에 따라 해가 없거나 유일하지 않을 수 있는 것과 대조적이다 — 볼테라에는 공명이 없다. 스펙트럼 반지름이 0 이라는 사실은 이산화하면 더 노골적으로 보인다. 행렬이 순하삼각(대각 아래만 0 아님)이 되어 멱영 행렬이 되고, 대각까지 포함하면 하삼각 행렬이라 전진 대입 한 번으로 풀린다.

단, “항상 수렴”이 “빨리 수렴”은 아니다. $(Mt)^n/n!$ 은 $n\approx Mt$ 까지 커졌다가 줄어들므로, 긴 구간에서 급수를 그대로 더하는 건 비효율적이고 상쇄 오차도 생긴다. 실전에서는 급수 대신 시간 행진을 한다.

3. 수치해법 — 시간 행진

격자 $t_n=nh$ 에서 적분을 사다리꼴 규칙으로 바꾸면

$$u_n=f_n+h\Big[\tfrac12K_{n,0}u_0+\sum_{j=1}^{n-1}K_{n,j}u_j+\tfrac12K_{n,n}u_n\Big]$$

이고, $u_n$ 이 양변에 있지만 선형이면 바로 나눠서 구한다. $u_0=f_0$ 에서 출발해 한 칸씩 전진한다. 검산 예제 $u(t)=1+\int_0^t u(s)\,ds$ (해 $e^t$) 에서 $t=1$ 의 오차는

$h$	오차	비율
0.1	$2.27\times10^{-3}$	—
0.05	$5.67\times10^{-4}$	4.0
0.025	$1.42\times10^{-4}$	4.0

로 2차 수렴이 깔끔하게 보인다(이 예제에서 사다리꼴 볼테라 해법은 상미분방정식의 사다리꼴 공식 $u_n=\big((1+h/2)/(1-h/2)\big)^n$ 과 정확히 같아진다).

일반론과 실무 포인트:

• 비용이 $O(N^2)$. 매 스텝 과거 전체를 더해야 하므로 $N$ 스텝이면 $N^2/2$ 번 커널 평가다. 메모리도 이력 전체 $O(N)$. 상미분방정식의 $O(N)$ 과 비교하면 기억의 대가다.
• 합성곱 커널 $K(t-s)$ 면 이력 합이 이산 합성곱이라 FFT 로 블록 단위 가속이 가능하다. 하이러·루비히·슐리히테(1985)의 방법이 $O(N\log^2 N)$ 이다.
• 고차 방법은 수치적분 규칙을 그대로 가져온 직교법(심프슨, 그레고리)과 배치법(collocation)이 있다. 브루너의 교과서가 배치법 수렴 이론을 정리한다. 단, 고차 직교법 일부는 시작값 문제와 불안정성이 있어 강성 방정식 쪽 안정성 분석을 그대로 해야 한다.
• 약특이 커널 $K\sim(t-s)^{-\alpha}$ 에서는 해가 $t=0$ 근처에서 매끄럽지 않아($u\sim t^{1-\alpha}$ 같은 항) 균등 격자의 수렴 차수가 떨어진다. 이때는 $t=0$ 쪽으로 몰린 등급 격자나 곱 적분(product integration)을 쓴다.
• 루비히의 합성곱 직교법(1988)은 커널의 라플라스 변환만 알면 되는 방식이라, 커널 자체가 닫힌 꼴로 없을 때도 쓴다. 경계요소법의 시간영역 판본과 분수계 도함수 이산화의 표준 도구다.

4. 제1종 — 미분이 숨어 있다

제1종 $\int_0^t K(t,s)u(s)\,ds=f(t)$ 를 $t$ 로 미분하면

$$K(t,t)\,u(t)+\int_0^t \partial_t K(t,s)\,u(s)\,ds=f'(t)$$

가 되어, $K(t,t)\ne0$ 이면 제2종으로 돌아온다. 대신 $f$ 를 한 번 미분한 셈이라 데이터 잡음이 증폭된다. 이것이 제1종 볼테라 방정식이 “약하게” 불량조건인 이유다 — 프레드홀름 제1종(매끄러운 커널)처럼 지수적으로 나쁘지는 않지만, $f$ 를 측정값으로 받는다면 티호노프 정규화나 미래 시점 정보를 쓰는 벡의 방법 같은 장치가 필요하다. 벽면 온도 이력에서 표면 열유속을 되찾는 역 열전도 문제가 대표적인 역문제다.

고전 예제는 아벨 방정식

$$\int_0^t\frac{u(s)}{\sqrt{t-s}}\,ds=f(t)\quad\Longrightarrow\quad u(t)=\frac1\pi\frac{d}{dt}\int_0^t\frac{f(s)}{\sqrt{t-s}}\,ds$$

다. 1823년 아벨이 등시강하(tautochrone) 곡선을 구하며 풀었고, 반차 적분·반차 미분의 원형이다.¹ 축대칭 물체의 투영 데이터에서 반경 분포를 되찾는 아벨 역변환(플라스마 진단, 연소 분광)도 같은 식이다.

5. 어디서 나오나

• 유전 효과·점탄성. 볼츠만 중첩 원리 $\sigma(t)=\int_0^t G(t-s)\,\dot\varepsilon(s)\,ds$ 가 볼테라 적분이다. 응력 이력에서 변형률을 되찾는 크리프 문제는 제1종을 푸는 일이다(점탄성).
• 인구·역학. 로트카의 갱신 방정식 $B(t)=g(t)+\int_0^t\beta(a)B(t-a)\,da$, 커맥-매켄드릭의 1927년 원 논문의 감염 연령 모형이 모두 제2종 볼테라다. 흔히 보는 SIR 상미분방정식은 감염 기간을 지수분포로 가정한 특수한 경우다.²
• 분수계 미분방정식. 카푸토 도함수는 약특이 커널을 가진 볼테라 연산자라, 분수계 미적분 문제는 결국 볼테라 방정식 수치해법이다.
• 일반화 랑주뱅 방정식. 기억 커널이 붙은 마찰항 $\int_0^t\gamma(t-s)v(s)\,ds$ 는 거친 입자 모형(coarse-graining)에서 나온다. 기억이 지수 감쇠면 보조 변수 하나로 마르코프화되어 오른슈타인-울렌벡 과정과 연결된다.
• 전자기·음향의 시간영역 적분방정식. 지연 퍼텐셜을 쓰는 시간영역 경계적분은 이력 합이 볼테라 구조다. 긴 시간에 발생하는 늦은 시간 불안정이 이 분야의 오래된 숙제다.³

6. 관련 문서

• 노이만 급수 · 스펙트럼 반지름 · 바나흐 공간
• 프레드홀름 적분방정식 · 경계요소법 · 그린 함수
• 상미분방정식 · 룽게-쿠타법 · 수치적분 · 강성 방정식
• 역문제 · 티호노프 정규화 · 점탄성
• 분수계 미적분 · 오른슈타인-울렌벡 과정

7. Footnotes

1. 어느 높이에서 놓든 바닥까지 같은 시간에 도착하는 곡선이 사이클로이드라는 건 하위헌스가 이미 알았다. 아벨은 그걸 “도착 시간 함수가 주어지면 곡선을 찾아라”라는 역문제로 일반화했고, 그 과정에서 반차 미분을 꺼냈다. 역문제의 원조가 진자시계에서 나온 셈. ↩

2. 실제 감염 기간은 지수분포가 아니라 꽤 뾰족하다. 그래서 SIR 상미분방정식은 전파 속도는 맞춰도 유행 곡선 모양이 어긋나곤 하는데, 원 논문은 처음부터 적분방정식이었다. 후대가 단순화하고 원전을 잊는 건 학계 국룰. ↩

3. 시간영역 적분방정식 코드를 오래 돌리면 아무 이유 없이 해가 지수적으로 불어나는 현상이 있다. 원인은 이산화된 볼테라 연산자의 미세한 불안정 모드이고, 수십 년째 “이번엔 진짜 고쳤다”는 논문이 나온다. ↩