Matlab으로 구현한 마티외 함수 툴박스

본 논문은 타원형 원통 좌표계에서 파동 방정식을 풀기 위해 필수적인 마티외 함수의 계산을 목표로, Matlab 기반의 전용 툴박스를 설계·구현하였다. 서론에서는 마티외 함수가 1868년 에밀 마티외에 의해 도입된 이래, 타원형 구조물(예: 타원형 파이프, 전자기 스캐터링 등)에서 해석적 해법을 제공하지만, 기존의 수치 라이브러리는 부족하거나 표기법·정규화가 일관되지 않아 활용이 제한돼 왔음을 지적한다. 이를 해소하기 위해 저자는 Stratton과 Stamnes의 표기법을 기본으로 하면서, 자체적인 간소화와 통일된 정규화 체계를 제시한다. II 장에서는 타원형 원통 좌표계(ξ, η, z)와 그 스케일 팩터 h_ξ, h_η를 정의하고, 파동 방정식(∇² + k²)U = 0을 변수분리하여 Z(z), S(v), R(u) 형태로 분리한다. 여기서 S(v)와 R(u)는 각각 각도·반경 마티외 방정식(6,7)으로, q = k²τ f²/4, a는 분리 상수이다. III 장에서는 각도 마티외 함수의 네 가지 주기해(짝‑짝, 짝‑홀, 홀‑짝, 홀‑홀)를 Sₚₘ(v,q,n) 형태로 정의하고, 전개계수 Aₖ를 재귀식(9‑12)으로 도출한다. 이 재귀식을 행렬 형태(13‑16)로 변환함으로써 고유값 문제를 명시하고, Matlab의 eig 함수를 이용해 특성값 a와 전개계수 벡터를 동시에 구한다. 특히 “1 짝‑짝” 카테고리에서 비대칭 행렬이 발생함을 언급하고, 필요 시 대칭화 변환을 적용할 수 있음을 제시한다. 정규화와 직교성은 Sₚₘ(0,q,n)=1 및 dSₚₘ/dv|_{v=π/2}=0 조건을 사용해 정의한다(17‑20). 이를 통해 각 카테고리별 정규화 상수 Nₚₘ을 구하고, 다른 문헌과의 비교 시 발생할 수 있는 스케일 차이를 명확히 한다. 또한 서로 다른 매질 영역에서 q가 달라질 때의 상관계수 Cₚₘ(q,q′,n)를 전개계수의 내적 형태(21‑22)로 정의하여, 경계 조건 적용 시 필요한 변환 계수를 쉽게 계산하도록 한다. IV 장에서는 반경 마티외 함수를 각도 함수와 동일한 전개계수를 사용해 Bessel 함수 Jₙ, Yₙ의 곱 형태로 전개한다. 이는 급수 수렴성을 크게 개선한다. 구체적으로 첫 번째 종류(Jₚₘ)와 두 번째 종류(Yₚₘ)를 (30‑32)식으로 제시하고, 결합계수 gₚₘ을 (24‑29)식으로 정의한다. u=0에서 Sₚₘ(0)=1, Jₚₘ(0)=1/(√2π gₚₘ) 관계를 이용해 gₚₘ을 전개계수와 특성값으로 직접 계산한다. 세 번째·네 번째 종류는 복소수 조합 H₁=J+iY, H₂=J−iY 로 정의되어, 파동 전파와 스캐터링 문제에서 복소수 파동함수로 바로 활용 가능하다. V 장에서는 Matlab 구현 세부 사항을 설명한다. K_F 파라미터를 통해 네 카테고리 중 하나를 선택하고, 전개계수 개수를 25로 고정한다. eig_spm 함수는 입력 q와 K_F를 받아 특성값 벡터 v_a, 전개계수 행렬 mc, 실제 차수 t를 반환한다. 반환된 전개계수는 정규화 조건(18)을 만족하도록 사후 처리된다. 이후 사용자는 제공된 함수들을 호출해 Sₚₘ(v,q,n), Rₚₘ(u,q,n), 그 미분, 정규화 상수, 상관계수 등을 손쉽게 얻을 수 있다. 논문 말미에는 다양한 q 값에 대한 수치표와 검증 결과가 제시되어, 구현의 정확성을 확인한다. 전체적으로 이 논문은 마티외 함수의 이론적 정의, 고유값 기반 전개계수 계산, 정규화·상관·미분 공식, Bessel 기반 급수 전개, 그리고 Matlab 코드 구현까지 전 과정을 일관된 프레임워크로 정리하였다. 이는 타원형 구조물의 전자기·음향·양자역학적 해석에 필요한 마티외 함수를 손쉽게 활용하고, 기존 문헌 간의 표기·정규화 차이로 인한 혼란을 최소화하려는 연구자들에게 매우 유용한 도구가 될 것이다.