다차원 라플라스 행렬 다항식 역연산 방법

초록

본 논문은 대칭 양정(positive definite) 라플라스 행렬의 다항식 형태에 대한 정확한 역행렬을 고유값·고유벡터 전개를 이용해 구하는 절차를 제시한다. 1차원 라플라스 행렬에 대한 기존 결과와 일치함을 보이며, 고차원 및 다항식 연산에 대한 일반화된 공식도 도출한다.

상세 분석

논문은 먼저 대칭 양정 행렬 A가 갖는 정규 직교 고유벡터 집합 {vₖ}와 실수 고유값 {λₖ}를 이용해 선형 시스템 Ax = b의 해를 고유벡터 전개 형태로 표현한다. 구체적으로 b를 고유벡터 기저에 대한 선형 결합 b = Σβₖ vₖ 로 쓰고, βₖ = (b, vₖ)/(vₖ, vₖ) 로 구한다. 이후 Ax = b 를 만족하는 x는 x = Σ (βₖ/λₖ) vₖ 로 얻어지며, 이는 고유값이 0이 아닌 경우에만 정의된다.

다음 단계에서는 A⁻¹ 자체를 고유값·고유벡터를 이용해 행렬 형태로 전개한다. A⁻¹의 k번째 열 gₖ는 Agₖ = eₖ 를 만족하므로, 위의 해식에 대입해 gₖ = Σ (vⱼₖ / λⱼ) vⱼ 로 얻는다. 여기서 vⱼₖ는 j번째 고유벡터의 k번째 성분이다. 따라서 A⁻¹의 (i,k) 원소는 (A⁻¹)ᵢₖ = Σ (vⱼₖ vⱼᵢ) / λⱼ 로 표현된다. 이 식은 A가 대칭 양정이므로 A⁻¹ 역시 대칭 양정임을 보장한다.

핵심적인 일반화는 행렬 다항식 P(A) = Σ a_m A^m 에 대해 고유값이 P(λₖ) 로 변환된다는 정리이다. 증명은 (1) A^m vₖ = λₖ^m vₖ 를 귀납적으로 보이고, (2) 선형 결합 형태의 다항식에 대해 P(A)vₖ = P(λₖ)vₖ 를 직접 계산함으로써 이루어진다. 따라서 P(A) 역시 원래 A와 동일한 고유벡터를 공유하고, 고유값은 다항식 P에 의해 변환된다.

이 결과를 이용하면 P(A)⁻¹ 의 원소는 (P(A)⁻¹)ᵢₖ = Σ (vⱼₖ vⱼᵢ) / P(λⱼ) 로 바로 구할 수 있다. 즉, 원래 행렬 A의 고유값·고유벡터만 알면, 어떠한 차수의 다항식 형태라도 역행렬을 명시적으로 구성할 수 있다.

논문은 1차원 라플라스 행렬에 대한 구체적인 사례를 제시한다. 연속 라플라스 연산자의 고유함수 ϕⱼ(x)=√2 sin(jπx) 를 격자점 xⱼ=jh (h=1/(n+1))에 투사하면 이산 라플라스 행렬 A의 고유벡터가 된다. 고유값은 λⱼ = 2h⁻²(1−cos(jπh)) = 4h⁻² sin²(jπh/2) 로 얻어진다. 이를 이용해 A⁻¹ 의 원소를 두 가지 형태, 즉 고유벡터 전개식과 기존에 알려진 삼각식 형태(예: (A⁻¹)ᵢₖ = h²/(n+1)·min(i,k)·(n+1−max(i,k)))가 동일함을 증명한다.

고차원 확장에 대해서는 기본해가 로그 형태가 되며, 전통적인 근사법(기본해 기반)에서는 높은 차수 미분이 남아 오차가 발생한다. 반면 고유값 전개 방식은 이러한 오차를 완전히 없앨 수 있다. 또한 Neumann 경계조건이 적용된 경우 A가 특이행렬이 되므로 역행렬이 존재하지 않지만, 다항식 형태의 연산에서는 여전히 유효한 해를 얻을 수 있다.

전체적으로 논문은 고유값·고유벡터 전개를 이용한 행렬 다항식 역연산의 이론적 기반을 명확히 하고, 실제 계산에서의 효율성과 정확성을 동시에 확보할 수 있음을 보여준다.