[논문 리뷰] On Wavelet-Galerkin methods for Semilinear Parabolic Equations with Additive Noise

Mihály Kovács, Stig Larsson외 1명|arXiv (Cornell University)|2012. 08. 02.

Image and Signal Denoising Methods참고 문헌 13인용 수 6

한 줄 요약

이 논문은 추가적인 소음을 가진 비선형 확률 파라볼릭 방정식을 해결하기 위해 파형-갈레르킨 방법의 새로운 접근을 제안한다. 문제는 선형 확률 방정식(비적응형 파형-갈레르킨으로 해결)과 비선형 랜덤 진화 방정식(적응형 파형-갈레르킨으로 해결)으로 분해된다. 주요 기여는 소음과 초기 자료의 정규성 조건 하에 최적 수렴 차수 $ \tau^{\beta/2} $ 를 보여주는 평균 제곱 오차 추정치를 제공하는 것이다. $ 0 \leq \beta < 1 $ 이다.

ABSTRACT

We consider the semilinear stochastic heat equation perturbed by additive noise. After time-discretization by Euler's method the equation is split into a linear stochastic equation and a non-linear random evolution equation. The linear stochastic equation is discretized in space by a non-adaptive wavelet-Galerkin method. This equation is solved first and its solution is substituted into the nonlinear random evolution equation, which is solved by an adaptive wavelet method. We provide mean square estimates for the overall error.

연구 동기 및 목표

• 추가 소음에 의해 영향을 받는 비선형 확률 파라볼릭 방정식을 위한 수치적 방법을 개발한다. 특히 소음이 거친 경우나 마르코프가 아닌 경우를 고려한다.
• 선형 확률 부분과 비선형 결정론적 부분을 분리하여 효율적으로 이러한 방정식을 해결하는 데 도전한다.
• 선형 확률 성분에 대해 비적응형 파형-갈레르킨을, 비선형 진화에 대해 적응형 파형-갈레르킨을 적용하여 서로 다른 정규성 성질을 활용한다.
• 전체 수치적 방법에 대해 엄밀한 평균 제곱 오차 추정치를 수립하며, 시간 및 공간 이산화 오차를 포함한다.
• 소음과 초기 자료의 정규성에 따라 수렴 속도를 도출하며, 공분산 연산자의 Hilbert-Schmidt 노름 조건에서 파arameter $ \beta $ 를 통해 정량화한다.

제안 방법

• Rothe 방법을 적용한다: 먼저 시간에 대해 역 오일러 스킴을 사용해 확률 편미분방정식을 이산화하고, 해를 $ u^n = v^n + w^n $ 로 분리한다. 여기서 $ w^n $ 은 선형 확률 문제를, $ v^n $ 은 비선형 문제를 각각 해결한다.
• 공간에 대해 선형 확률 방정식 $ w^n $ 을 비적응형 파형-갈레르킨 방법으로 이산화한다. 다중해상도 공간 $ S_J $ 의 순서 $ m $ 을 기반으로 하며, 안정성과 수렴성을 보장한다.
• 비선형 방정식 $ v^n $ 은 적응형 파형-갈레르킨 방법으로 이산화한다. 이는 해의 특징을 효율적으로 포착하기 위해 공간 기저를 동적으로 정밀화한다.
• SPDE의 온화한 해 표현식을 사용한다. 이를 통해 확률 적분 $ w(t) = \int_0^t e^{-(t-s)A} dW(s) $ 과 $ v $ 를 위한 랜덤 진화 방정식으로 분리하여 순차적 해법이 가능하도록 한다.
• 시간에 대해 반세미군자 $ e^{-tA} $ 를 근사하기 위해 해석 연산자 $ r^n(\tau A) $ 를 적용하고, 파형-갈레르킨 사영의 안정성 및 근사 성질을 활용한다.
• 이토 등식과 스펙트럼 함수 해석학을 활용해 확률 적분의 오차를 제한한다. Hilbert-Schmidt 노름 조건 $ \|A^{(\beta-1)/2}Q^{1/2}\|_{\mathrm{HS}} < \infty $ 을 활용한다.

실험 결과

연구 질문

• RQ1 파형-갈레르킨 방법이 추가 소음을 가진 비선형 확률 파라볼릭 방정식에 대해 최적 수렴 속도를 달성할 수 있는가?
• RQ2 선형 확률 성분과 비선형 결정론적 성분을 분리함으로써 전체 오차와 수치적 안정성에 어떤 영향을 미치는가?
• RQ3 소음과 초기 자료의 정규성 조건 하에 시간 간격 $ \tau $ 에 따른 파형-갈레르킨 방법의 최적 수렴 속도는 무엇인가?
• RQ4 이러한 유형의 S-PDE에 대해 비적응형과 적응형 파형-갈레르킨 접근 방식은 오차 제어 및 계산 효율성 측면에서 어떻게 비교되는가?
• RQ5 공분산 연산자 $ Q $ 와 초기 자료에 대해 어떤 조건이면 수치적 방법이 평균 제곱 수렴성을 유지하는가?

주요 결과

• 소음과 초기 자료의 정규성 조건 하에 $ \|A^{(\beta-1)/2}Q^{1/2}\|_{\mathrm{HS}} < \infty $ 이면 전체 수치적 해의 평균 제곱 오차는 $ \mathbb{E}\|e^n\|^2 \leq C\tau^\beta \mathrm{e}^{Ct_n} $ 를 만족하며, $ 0 \leq \beta < 1 $ 이다.
• 평균 제곱 노름에서 수렴 속도는 $ \tau^{\beta/2} $ 이며, 이는 주어진 소음과 초기 자료의 정규성 조건 하에 최적이 된다.
• 조건 $ \tau L_f \leq 1/2 $ 하에 비선형 성분의 오차 $ \bar{v}^n - v^n $ 는 파형 근사에서의 확률 적분 오차 $ w_J^n - w^n $ 의 상수 배수로 제한된다.
• 시간 이산화된 확률 적분의 오차는 $ \mathbb{E}\|e_{21}\|^2 \leq C\tau^\beta \|A^{(\beta-1)/2}Q^{1/2}\|_{\mathrm{HS}}^2 $ 로 제한되며, 이는 소음의 정규성에 따른 수렴 영향을 보여준다.
• 반세미군의 시간 이산화 오차는 비율 근사 $ r^n(\tau A) $ 의 안정성에 의해 제어되며, 이는 이산 그론월의 부등식을 통해 유도된 경계를 포함한다.
• 분석 결과 비선형 성분에 대해 적응형 파형-갈레르킨 방법이 효과적임을 확인하였으며, 오차 전파가 선형적으로 비선형 성분의 오차에 의존하여 안정성과 수렴성을 보장한다.