행렬 기반 이산 분수 미분 방정식 접근법 II: 부분 분수 미분 방정식

초록

본 논문은 Podlubny의 행렬 접근법을 확장하여, 시간·공간에 대한 임의 실수 차수의 분수 미분과 지연 항을 포함하는 편미분 방정식을 손쉽게 이산화하는 방법을 제시한다. MATLAB 구현 코드를 포함한 다섯 가지 예제로 방법의 일반성 및 효율성을 검증한다.

상세 분석

본 연구는 기존 Podlubny(2000)의 행렬 기반 분수 미분 연산자를 다변수 상황에 적용하기 위한 이론적·수치적 확장을 수행한다. 핵심 아이디어는 1차원 구간을 N개의 격자로 나누고, 각 격점에서의 함수값을 열벡터 u 로 표현한 뒤, 분수 차수 α 에 대해 Grünwald‑Letnikov(또는 Caputo) 정의를 행렬 Aα 로 전개하는 것이다. 이때 Aα 는 Toeplitz 구조를 가지며, 고정된 격자 간격 h 에 대해 계수 c_k = (−1)^k · (α choose k) 를 이용해 구성한다. 다변수 경우, 시간과 공간 각각에 대해 독립적인 행렬 Aα_t, Aβ_x 를 만든 뒤, 크로네커 곱(Kronecker product) ⊗ 을 사용해 전체 연산자를 구축한다. 이렇게 하면 2‑D 혹은 3‑D 분수 편미분 연산자를 하나의 큰 희소 행렬로 표현할 수 있어, 기존 유한 차분법에서 발생하는 복잡한 가중치 계산을 일괄적으로 행렬 연산으로 대체한다.

시간‑공간 분리형 PDE ∂^α_t u = κ ∂^β_x u + f 에 대해, 시간 행렬 Aα_t 와 공간 행렬 Aβ_x 를 각각 전진·후진 차분 스킴에 맞게 구성하고, 초기·경계 조건을 벡터 u₀, g 으로 삽입한다. 전체 시스템은 (I − Δt^α Aα_t) u^{n+1} = u^{n} + Δt^α κ Aβ_x u^{n} + Δt^α f^{n} 형태가 되며, 이는 선형(또는 비선형) 시스템을 직접 풀 수 있게 만든다.

지연 항이 포함된 경우, 예를 들어 ∂^α_t u(t) = κ ∂^β_x u(t − τ) + f(t) 에서는 지연 τ 를 격자 간격 Δt 의 정수배 m 으로 가정하고, 이전 m 스텝의 해를 저장하는 버퍼 벡터 U_delay 를 도입한다. 행렬식은 변하지 않으며, 오른쪽 항에 U_delay 를 삽입함으로써 동일한 연산 흐름을 유지한다.

수치 실험에서는 (1) 정수‑정수, (2) 분수‑정수, (3) 정수‑분수, (4) 분수‑분수 조합의 1‑D 및 2‑D 확산 방정식을 풀었으며, 각 경우에 대해 정확해와 비교한 오차 수렴률이 이론적 O(h^{min(α,β)}) 에 부합함을 확인했다. 또한, 지연 사례에서는 시간 지연 τ 가 0.5·T 인 경우에도 안정적인 수렴을 보이며, 기존 방법에 비해 구현 코드 라인이 30 % 정도 감소하고, 행렬 희소성 덕분에 메모리 사용량이 크게 절감되는 장점을 입증했다.

MATLAB 구현은 genFracMatrix(α,N,Δt) 함수로 행렬을 자동 생성하고, solveFracPDE 함수로 전체 시스템을 풀도록 설계되었으며, 샘플 코드와 함께 공개되어 재현성을 높였다. 이 접근법은 복잡한 경계 조건, 비선형 반응항, 다중 지연 등으로 확장 가능함을 논의하고, 향후 고성능 병렬 구현 및 비정규 격자 적용 가능성도 제시한다.