포트 해밀턴 시스템의 반균일 안정성 정량화와 디오판틴 근사

초록

본 논문은 1차원 포트‑해밀턴 시스템에서 발생하는 C₀‑반연속 반군 semigroup의 반균일 안정성을 정량적으로 규명한다. 해석자는 resolvent 연산자 ‖R(i t,−A)‖의 성장함수 M(η)를 시스템 행렬의 역행렬 ‖Tₜ⁻¹‖와 동등하게 추정함으로써, 에너지 감쇠율을 정확히 제어한다. 또한 α∈ℝ⁺ 파라미터에 따라 강안정성은 α가 무리수인지에 달려 있음을 보이고, α의 디오판틴 근사 정도에 따라 t^{-β} (β<½) 형태의 임의의 다항 감쇠율을 구현하는 보편적인 예시를 제시한다.

상세 분석

논문은 먼저 포트‑해밀턴 연산자 A를 정의하고, A가 강연속 수축 semigroup T(t)=e^{-tA}를 생성함을 확인한다. 핵심은 iℝ이 A의 resolvent 집합에 포함될 때, 즉 강안정성(스펙트럼이 허수축에 없을 때)이라면, 실수 t에 대해 ‖R(i t,−A)‖을 두 개의 행렬량 ‖Φ_t‖와 ‖T_t⁻¹‖로 정확히 상·하한을 잡을 수 있다는 정리(2.3)를 증명한다. 여기서 Φ_t는 시스템의 기본 행렬이며, T_t = W Φ_t(b)I_d는 경계 행렬 W와 Φ_t의 곱으로 정의된다. Lemma 2.5와 Lemma 2.6을 통해 ‖R(i t,−A)‖ ≤ C(‖T_t⁻¹‖+1)와 ‖T_t⁻¹‖ ≤ C_t(‖R(i t,−A)‖+1) 를 얻어, 두 함수 M(η)=sup_{|t|≤η}‖R(i t,−A)‖와 m(η)=sup_{|t|≤η}‖T_t⁻¹‖가 서로 동등함을 보인다. 이 결과는 기존의 지수 안정성 분석을 일반화하여, “양의 증가” 성질을 갖는 함수 γ에 대해 γ가 양의 증가이면 M도 양의 증가함을 즉시 얻는다(Remark 2.4).

다음 단계에서는 구체적인 예시를 제시한다. H는 상수 행렬 diag(1,α)⁻¹, P₁=I₂, P₀=0, W은 특정 고정 행렬로 잡는다. 이 경우 Φ_t = e^{i t H^{-1}}이며, T_t = f_W e^{i t H^{-1}}. 행렬식 det T_t = 1 + ½(e^{i t}+e^{i α t}) 로부터, α가 유리수이면 허수축에 고유값이 존재해 강안정성이 깨지고, α가 무리수이면 스펙트럼이 허수축을 피하므로 강안정성을 갖는다(정리 3.1). 그러나 어떤 α에 대해서도 지수적 감쇠는 일어나지 않는다(정리 3.1, 명제 3.1).

핵심은 α의 디오판틴 근사 성질이다. α가 “badly approximable” (즉, |α−p/q| ≥ c/q²)이면 ‖T_t⁻¹‖는 t^{½} 정도까지만 성장하고, 따라서 M(η)≈η^{½}가 된다. 이때 반균일 안정성에 대한 일반적인 정리(부록 A)로부터 ‖T(t)A^{-1}‖ = O(t^{-½}) 를 얻는다. 반대로 α를 연속분수 전개에서 부분수(p_n/q_n)로 선택하면, q_n가 원하는 속도로 성장하도록 조정할 수 있다. 이렇게 하면 ‖T_t⁻¹‖는 η^{β} (0<β<½) 수준으로 제한되며, 결과적으로 semigroup은 t^{-β} 속도의 다항 감쇠를 보인다. 따라서 “t^{-½}보다 느린 임의의 감쇠율”을 구현하는 보편적인 예시가 존재함을 증명한다(정리 3.2, 3.3).

이러한 결과는 포트‑해밀턴 시스템이 복잡한 PDE(예: 파동, 유체역학)의 경계 제어 모델에 적용될 때, 에너지 감쇠를 정량적으로 예측할 수 있는 강력한 도구가 된다. 특히, 시스템 매개변수(여기서는 α)의 수론적 특성이 물리적 감쇠율에 직접적인 영향을 미친다는 점은 기존 PDE 안정성 이론과는 다른 새로운 관점을 제공한다.