브론스키 행렬식

1. 개요[편집]

브론스키 행렬식 Wronskian
정의	W(f₁,…,fn) = det[ fj(i−1) ]
이름	Józef Hoene-Wroński (1812), 명명은 Thomas Muir (1882)
핵심 항등식	아벨 항등식 W′ = −p(x) W
따름정리	한 점에서 0이면 구간 전체에서 0
쓰이는 곳	매개변수 변환법 · 그린 함수 · 스투름 진동 정리
수치적으로는	직접 계산 금지 — 지수적으로 폭발/소멸

브론스키 행렬식(Wronskian)은 $n$ 개의 함수 $f_1,\dots,f_n$ 을 세로로 세우고 도함수를 계속 쌓아 만든 $n\times n$ 행렬식이다.

$$W(f_1,\dots,f_n)(x) \;=\; \begin{vmatrix} f_1 & f_2 & \cdots & f_n \\ f_1' & f_2' & \cdots & f_n' \\ \vdots & \vdots & & \vdots \\ f_1^{(n-1)} & f_2^{(n-1)} & \cdots & f_n^{(n-1)} \end{vmatrix}$$

즉 $(i,j)$ 성분이 $f_j^{(i-1)}$ 이다. 만들어진 목적은 명확하다 — 함수들이 선형독립인지 판정하기 위해서다. 함수 $n$ 개의 선형결합이 0이라는 조건을 $n-1$ 번 미분해 얻은 연립방정식의 계수행렬이 위 행렬이고, 그 행렬식이 0이 아니면 계수가 전부 0일 수밖에 없다.

문제는 이 판정이 한쪽 방향으로만 성립한다는 것이다. 교과서가 ”$W \ne 0$ 이면 독립”까지만 조심스럽게 쓰고 역을 얼버무리는 데는 이유가 있으며, 그 이유를 정확히 아는 것이 이 문서의 절반이다. 나머지 절반은 아벨 항등식이 브론스키 행렬식을 직접 계산하지 않고도 알 수 있게 만든다는 사실, 그리고 그럼에도 수치적으로는 이 값을 절대 직접 굴리면 안 된다는 실무적 결론이다.

2. 선형독립 판정 — 정확히 어디까지 참인가[편집]

항상 참인 방향. $f_1,\dots,f_n$ 이 구간 $I$ 에서 선형종속이면 $W \equiv 0$ 이다. 자명하다 — 계수 $c_j$ 로 $\sum c_j f_j \equiv 0$ 이면 미분해도 성립하므로 행렬의 열들이 같은 계수로 종속이다. 대우를 취하면 어떤 한 점에서라도 $W \ne 0$ 이면 독립이고, 이것이 실무에서 쓰는 유일한 방향이다.

거짓인 방향. $W \equiv 0$ 이라고 종속인 것은 아니다. 표준 반례는 $I = (-1,1)$ 위의

$$f_1(x) = x^2, \qquad f_2(x) = x|x| = \begin{cases} x^2 & x \ge 0\\ -x^2 & x < 0\end{cases}$$

둘 다 $C^1$ 이고, $x \ge 0$ 에서는 $f_2 = f_1$, $x<0$ 에서는 $f_2 = -f_1$ 이므로 각 반구간에서 종속이라 $W \equiv 0$ 이다. 그런데 전체 구간에서는 독립이다 — $c_1 f_1 + c_2 f_2 \equiv 0$ 을 $x=1$ 과 $x=-1$ 에 대입하면 $c_1+c_2=0$ 과 $c_1-c_2=0$ 이 나와 $c_1=c_2=0$ 이다. 종속성이 “구간마다 다른 계수로” 성립하기 때문에 국소적으로만 보는 브론스키 행렬식이 이를 놓친다. 페아노가 1889년에 이 함정을 지적했다.¹

역이 살아나는 두 경우가 있고, 실전에서는 대개 이 둘 중 하나다.

• 해석함수. $f_j$ 가 구간에서 해석적이면 $W \equiv 0 \iff$ 선형종속이다. 위 반례의 $x|x|$ 는 원점에서 두 번 미분이 안 되므로 해석적이지 않고, 해석성이 “국소 정보가 전역을 결정한다”를 보증하기 때문에 반례가 원천 차단된다.
• 같은 선형 미분방정식의 해. $f_j$ 가 계수가 연속인 $n$ 계 선형 동차 방정식 $$y^{(n)} + a_{n-1}(x)\,y^{(n-1)} + \cdots + a_0(x)\,y = 0$$ 의 해라면, $W \equiv 0 \iff$ 어떤 한 점에서 $W = 0 \iff$ 선형종속이다. 이유는 초기값 문제의 유일성이다 — $W(x_0)=0$ 이면 $x_0$ 에서 초기데이터가 종속인 해가 있고, 유일성에 의해 그 해는 전 구간에서 0이다.

두 번째 항목이 미분방정식 교과서가 브론스키 행렬식을 편하게 쓰는 근거다. 기본해계인지 확인하려면 아무 점이나 하나 잡아 $W$ 를 계산하면 된다. 임의의 함수 뭉치에 대해서는 못 하는 일을, 해 공간 안에서는 마음껏 한다.

3. 아벨 항등식 — W를 풀지 않고 안다[편집]

2계 방정식 $y'' + p(x)y' + q(x)y = 0$ 의 두 해에 대해 $W = y_1y_2' - y_2y_1'$ 를 미분하면 $y_1y_2''-y_2y_1''$ 만 남고, 방정식을 대입하면 $q$ 항이 상쇄되어

$$W'(x) = -p(x)\,W(x) \qquad\Longrightarrow\qquad W(x) = W(x_0)\,\exp\!\left(-\int_{x_0}^{x} p(t)\,dt\right)$$

를 얻는다. 이것이 아벨 항등식(Abel’s identity, 아벨-야코비-리우빌 공식)이다. $n$ 계로 올리면 $W' = -a_{n-1}W$ 이고, 1계 연립계 $\mathbf y' = A(t)\mathbf y$ 의 기본행렬 $\Phi$ 에 대해서는

$$\frac{d}{dt}\det\Phi(t) = \operatorname{tr}A(t)\,\det\Phi(t) \qquad\Longrightarrow\qquad \det\Phi(t) = \det\Phi(t_0)\exp\!\left(\int_{t_0}^{t}\operatorname{tr}A(s)\,ds\right)$$

가 된다(야코비 공식). $n$ 계 방정식을 동반형 1계계로 바꾸면 $\operatorname{tr}A = -a_{n-1}$ 이라 두 식이 정확히 일치한다. 플로케 이론에서 모노드로미 행렬이 항상 정칙임을 보증하는 것이 바로 이 공식이고, 해밀토니안계에서 $\operatorname{tr}A = 0$ 이라 $\det\Phi \equiv 1$ 이 되는 것이 리우빌 정리의 선형화 판이다.

여기서 나오는 따름정리가 실용적으로 가장 중요하다.

$W$ 는 어디서도 0이 아니거나, 항등적으로 0이거나 둘 중 하나다. 중간은 없다.

지수함수는 절대 0이 되지 않으므로 $W(x_0)=0$ 이면 전체가 0, 아니면 전체가 0이 아니다. “구간 어딘가에서만 0”인 브론스키 행렬식은 (같은 방정식의 해에 대해서는) 존재할 수 없다. 그래서 기본해계 판정에 계산하기 제일 편한 점 하나를 골라도 되는 것이다.

스투름-리우빌 형식에서는 더 예쁜 꼴이 된다. $(p(x)y')' + q(x)y = 0$ 을 표준형으로 바꾸면 1계 계수가 $p'/p$ 이므로 아벨 항등식이 $W' = -(p'/p)W$, 즉

$$\bigl(p(x)\,W(x)\bigr)' = 0 \qquad\Longrightarrow\qquad p(x)\,W(x) = \text{상수}$$

가 된다. 아래에서 볼 그린 함수의 분모에 등장하는 그 상수다.

4. 매개변수 변환법과 그린 함수[편집]

브론스키 행렬식이 계산에 실제로 쓰이는 첫 번째 자리. 비동차 방정식 $y'' + p y' + q y = g$ 의 특수해를 구할 때, 기본해 $y_1,y_2$ 에 대해 $y_p = u_1y_1 + u_2y_2$ 로 놓고 $u_1'y_1+u_2'y_2=0$ 이라는 편의 조건을 붙이면

$$\begin{bmatrix} y_1 & y_2 \\ y_1' & y_2' \end{bmatrix} \begin{bmatrix} u_1' \\ u_2' \end{bmatrix} = \begin{bmatrix} 0 \\ g \end{bmatrix}$$

라는 $2\times2$ 선형계가 나오고, 크라메르 법칙으로 즉시 풀린다.

$$u_1' = -\frac{y_2\,g}{W},\qquad u_2' = \frac{y_1\,g}{W} \qquad\Longrightarrow\qquad y_p(x) = -y_1(x)\!\int\!\frac{y_2 g}{W}\,dx \;+\; y_2(x)\!\int\!\frac{y_1 g}{W}\,dx$$

매개변수 변화법의 그 외우기 싫은 부호 배치는 크라메르 법칙의 부호 그 자체이고, 분모의 $W$ 가 0이 아니라는 조건이 곧 $y_1,y_2$ 가 기본해계라는 조건이다. 공식을 외우는 대신 $2\times2$ 계를 세우는 습관을 들이면 부호를 틀릴 일이 없다.

같은 구조가 경계값 문제의 그린 함수로 이어진다. 스투름-리우빌 연산자 $L y = (py')'+qy$ 에 대해 좌·우 경계조건을 각각 만족하는 해 $y_1, y_2$ 를 잡으면

$$G(x,\xi) = \frac{y_1(x_<)\,y_2(x_>)}{p(\xi)\,W(\xi)}, \qquad x_< = \min(x,\xi),\ x_> = \max(x,\xi)$$

이고, 앞 절에서 본 대로 분모 $pW$ 가 상수라 $\xi$ 에 무관하다. 그린 함수 유도에서 ”$pW$ 는 상수이므로”라는 문장이 아무 설명 없이 지나가는 것을 본 적이 있다면, 그 근거가 아벨 항등식이다. $y_1$ 과 $y_2$ 가 종속이면 $W=0$ 이 되어 그린 함수가 존재하지 않는데, 이는 곧 동차 문제에 자명하지 않은 해가 있다는 뜻이고 원래 경계값 문제가 해를 갖지 않거나 유일하지 않다는 뜻이다. 프레드홀름 택일이 분모 하나로 드러난다.

5. 스투름 진동 이론 — 영점의 배치[편집]

브론스키 행렬식이 정성적 정보를 주는 대표 사례. $y'' + q(x)y = 0$ 의 선형독립인 두 해 $y_1, y_2$ 를 잡으면 $W$ 는 부호를 바꾸지 않는 상수(아벨 항등식, $p \equiv 0$)다. 이제 $y_1$ 의 연속한 두 영점 $a < b$ 를 잡으면

$$W(a) = y_1(a)y_2'(a) - y_2(a)y_1'(a) = -y_2(a)\,y_1'(a), \qquad W(b) = -y_2(b)\,y_1'(b)$$

인데, $a$ 와 $b$ 가 연속한 영점이므로 $y_1'(a)$ 와 $y_1'(b)$ 는 부호가 반대다. $W(a)=W(b)\ne0$ 이려면 $y_2(a)$ 와 $y_2(b)$ 도 부호가 반대여야 하고, 따라서 중간값 정리로 $(a,b)$ 안에 $y_2$ 의 영점이 있다. 이것이 스투름 분리 정리 — 독립인 두 해의 영점은 서로 엇갈려 배치된다. 한 줄짜리 행렬식이 진동의 구조를 결정한다.

여기서 스투름 비교 정리로, 다시 스투름-리우빌 이론의 고유값 문제로 이어진다. 고유함수 $\phi_k$ 의 영점 개수가 정확히 $k$ 개라는 진동 정리, 고유값이 아래로 유계이고 무한대로 발산한다는 정리, 그리고 그 고유함수계가 완비 직교계를 이룬다는 정리가 전부 이 계보 위에 있다. 물리적으로는 “고차 모드일수록 마디가 많다”는 익숙한 사실이고, 수치적으로는 고유값 문제 풀이에서 부호 변화 개수를 세어 몇 번째 고유값인지 판정하는 슈투름 수열 기법으로 나타난다.

6. 수치적으로는 계산하지 않는다[편집]

지금까지가 이론이고, 여기서부터가 실무의 결론이다. 브론스키 행렬식을 부동소수점으로 직접 적분해서 쓰는 코드는 거의 전부 틀린다.

원인은 아벨 항등식이 그대로 알려 준다.

$$W(x) = W(x_0)\exp\!\left(\int_{x_0}^{x}\operatorname{tr}A\right)$$

지수 안의 적분이 조금만 커도 $W$ 는 지수적으로 폭발하거나 소멸한다. 강성이 큰 문제에서 $\operatorname{tr}A$ 가 $-10^3$ 급이면 몇 단위 구간만 가도 $W$ 가 배정밀도 최소 정규수 아래로 내려가 언더플로하고, 부호가 반대면 오버플로해서 inf 가 된다. 값이 살아 있어도 문제다 — $W$ 는 0에 “거의” 가까울 수 없다는 이론적 사실이 있지만, 부동소수점에서는 $10^{-300}$ 과 0을 구별할 방법이 사실상 없다.

더 근본적인 병은 해들이 서로 평행해진다는 것이다. 기본행렬 $\Phi$ 의 열들을 각각 적분하면 모든 열이 지배적 성장률 $\lambda_1$ 방향으로 정렬되어, $\det\Phi$ 는 (이론적으로는) 0이 아닌데 계산된 행렬은 수치적으로 랭크 1이 된다. 조건수가 $\exp\bigl((\lambda_1-\lambda_n)t\bigr)$ 처럼 자라므로, 유효숫자 16자리는 지수가 $\ln 10^{16} \approx 37$ 을 넘는 순간 전부 소진된다. 랴푸노프 지수 계산에서 접선 벡터들이 최대 지수 방향으로 정렬되어 버리는 그 현상과 완전히 같은 병이다.

그래서 실무는 브론스키 행렬식을 우회한다.

• 연속 직교화(연속 QR). 기본행렬을 그대로 적분하는 대신 $\Phi = QR$ 로 두고 $Q$ 와 $R$ 을 각각 적분하거나, 일정 간격마다 그람-슈미트·QR 분해로 재직교화한다. 성장은 $R$ 의 대각에 로그로 흡수되고 $Q$ 는 항상 직교라 조건수가 1이다. 이 방법이 경계값 문제의 슈팅법에서는 고두노프-콘테 직교화라는 이름으로, 카오스 쪽에서는 벤티틴 알고리즘이라는 이름으로 각각 독립 발견됐다.
• 리카티 변환. 해공간을 좌표 대신 “부분공간”으로 추적한다. $\Phi = \begin{bmatrix}Y_1 \\ Y_2\end{bmatrix}$ 를 $S = Y_2Y_1^{-1}$ 로 바꾸면 $S$ 가 행렬 리카티 방정식 $S' = A_{21} + A_{22}S - SA_{11} - SA_{12}S$ 를 만족하고, 크기 정보가 소거되어 지수적 성장이 사라진다. 대가는 $Y_1$ 이 특이해지는 지점에서 $S$ 가 발산한다는 것이라, 그때마다 좌표를 갈아타는 재시작 전략이 필요하다. 최적제어의 리카티 방정식과 같은 물건이 미분방정식 수치해에도 등장하는 것이다.
• 복합행렬법. 해 부분공간의 $k\times k$ 소행렬식(플뤼커 좌표)들을 미지수로 삼아 직접 적분한다. 크기를 정규화하기 쉽고 부분공간 정보만 남으므로 안정적이며, 유체 안정성의 오어-조머펠트 방정식처럼 성장률 차이가 극단적인 문제의 고전적 해법이다. 자세한 구성은 복합행렬 쪽 이야기다.

선형독립을 수치적으로 확인하고 싶을 때도 브론스키 행렬식은 쓰지 않는다. 표본점에서 함수값 행렬을 만들어 QR 분해나 특이값을 보는 쪽이 훨씬 낫다. 그람 행렬 $G_{ij} = \int f_i f_j$ 의 최소 고유값을 보는 방법도 있는데, 이쪽은 “얼마나 독립인가”라는 정량적 답까지 준다. 브론스키 행렬식의 답은 “0인가 아닌가”라는 이분법인데, 부동소수점 세계에서 0 판정은 원래 성립하지 않는 질문이다.²

그럼에도 아벨 항등식은 코드에 하나 넣어 둘 가치가 있다. 적분한 기본행렬의 행렬식이 $\exp(\int\operatorname{tr}A)$ 와 몇 자리나 어긋나는지 재 보면 적분기의 정확도를 공짜로 감시할 수 있다. 심플렉틱계라면 $\det\Phi \equiv 1$ 이므로 더 간단하다. 세상에 이렇게 값싼 단위 테스트가 흔치 않다.³

7. 관련 문서[편집]

• 행렬식 · 크라메르 법칙 · 매개변수 변화법
• 스투름-리우빌 이론 · 그린 함수 · 고유값 문제
• 플로케 이론 · 리우빌 정리 · 랴푸노프 지수
• 리카티 방정식 · 복합행렬 · 오어-조머펠트 방정식
• QR 분해 · 그람-슈미트 · 조건수
• 룽게-쿠타법 · 직교 배치법 · 편미분방정식

8. Footnotes[편집]

1. 이름의 주인 유제프 호에네-브론스키는 폴란드 출신으로, 수학보다 자기만의 “절대 철학” 체계를 세우는 데 더 열심이었고 급수 전개에 관한 주장 상당수가 틀린 것으로 판명됐다. 정작 이 행렬식에 그의 이름을 붙인 사람은 70년 뒤의 토머스 뮤어(1882)다. 브론스키 본인은 자기 이름이 붙은 이 물건으로 선형독립을 판정할 생각이 없었다는 점에서, 스티글러의 명명 법칙이 두 겹으로 적용된 사례라 할 만하다. ↩

2. 그래서 기호계산에서는 사정이 정반대다. 컴퓨터 대수 시스템이나 유한체 위의 정확 산술에서는 ”$W$ 가 0인가”가 판정 가능한 질문이고, 미분 갈루아 이론이나 홀로노믹 함수 계산에서 브론스키 행렬식이 실제 알고리즘의 일부로 돌아간다. 이 문서가 “쓰지 말라”고 하는 것은 어디까지나 부동소수점 이야기다. ↩

3. 반대로 이 검사를 통과했다고 적분이 정확한 것도 아니다. 행렬식 하나는 스칼라라서 $n^2$ 개 성분의 오차를 다 잡아내지 못하고, 특히 오차가 부피를 보존하는 방향으로만 쌓이면 통과한다. 심플렉틱 적분기가 에너지는 잘 지키면서 위상은 서서히 어긋나는 것과 같은 종류의 눈속임이다. 값싼 테스트는 값싼 만큼만 알려 준다. ↩