[논문 리뷰] Stability of densities for perturbed Diffusions and Markov Chains

Valentin Konakov, Anna Kozhina외 1명|arXiv (Cornell University)|2015. 06. 29.

Advanced Mathematical Modeling in Engineering참고 문헌 31인용 수 3

한 줄 요약

이 논문은 비퇴화 확산 과정과 마르코프 체인의 전이 밀도에 대해 계수의 변화에 대한 정량적 안정성 경계를 수립한다. 파라메트릭스 방법과 헬더 연속성 가정을 사용하여 계수 차이의 함수로 명시적인 오차 추정치를 도출하며, 이는 비선형성에 따라 드리프트 및 확산 계수의 정규성에 따라 수렴 속도가 결정된다. 특히 헬더 또는 다항형 감쇠 조건 하에서 성립한다.

ABSTRACT

We are interested in studying the sensitivity of diffusion processes or their approximations by Markov Chains with respect to a perturbation of the coefficients.

연구 동기 및 목표

• 비퇴화 확산 및 마르코프 체인의 전이 밀도가 드리프트 및 확산 계수의 변화에 얼마나 민감한지 정량화하는 것.
• 계수에 대한 최소한의 정규성 가정 하에서 진짜 밀도와 편향된 밀도 간의 차이에 대해 명시적인 오차 경계를 제공하는 것.
• 특히 옐러 방법 약오차 분석의 맥락에서 확산 방정식의 이산 시간 마르코프 체인 근사에 안정성 결과를 확장하는 것.
• 모델 부정확성, 매개변수 추정, 수치적 방법 수렴성과 같은 적용 분야에서 계수 변화의 영향을 분석하는 것.
• 헬더 연속성 및 다항 감쇠 조건 하에서 밀도 차이의 수렴 속도를 수립하고, 이는 편향 크기와 정규성 매개변수에 명시적인 의존성을 가진다.

제안 방법

• 확산 및 마르코프 체인의 전이 밀도를 가우시안 커널과 반복 적분을 포함하는 급수 전개를 통해 표현하기 위해 파라메트릭스 방법을 사용한다.
• 맥케인–싱어 접근법을 파라메트릭스 전개에 적용하여 시간에 대한 스무스함 부족에 강건하며 이산 시간 마르코프 체인에 적응 가능하도록 한다.
• 특히 두 번째 및 세 번째 도함수에 중점을 두어 파라메트릭스 급수 항의 변화를 분석함으로써 밀도 차이에 대한 경계를 도출한다.
• 특히 헬더 연속성 가정 하에서 드리프트 및 확산 계수의 차이를 제어하기 위해 헬더 부등식과 가우시안 유형 경계를 적용한다.
• 다항 감쇠 성질을 가진 밀도(예: $ q_m(t,x) riangleq c^d (t)^{-d/2} Q_{M-(d+5+ u)}(x / t^{1/2}) $)에 대한 컨볼루션 분석을 도입하여 컨볼루션 적분에서 대각선 및 비대각선 영역을 다룬다.
• 귀납법과 감마 함수 항등식을 사용하여 파라메트릭스 전개에서 반복된 컨볼루션의 성장을 통제하고, 이로 인해 $ au^{rac{r heta}{2}} $ 형태의 명시적 수렴 속도를 도출한다.

실험 결과

연구 질문

• RQ1 비퇴화 확산 과정의 전이 밀도는 드리프트 및 확산 계수에 작은 변화가 가해질 경우 어떻게 변화하는가?
• RQ2 계수에 변화가 가해졌을 때, 마르코프 체인 근사의 밀도가 진짜 확산 밀도로 수렴하는 속도는 어떠한가?
• RQ3 계수에 대한 헬더 연속성과 다항 감쇠 조건은 밀도의 안정성에 어떻게 영향을 미치는가?
• RQ4 편향된 계수를 가진 옐러 방법의 약오차에 대해 파라메트릭스 방법을 사용하여 명시적인 오차 경계를 유도할 수 있는가?
• RQ5 고차 도함수는 밀도 차이에 어떤 역할을 하는가? 그리고 파라메트릭스 전개에서 농도 손실은 어떻게 통제되는가?

주요 결과

• 계수의 헬더 연속성 하에서, 진짜 밀도와 편향된 밀도 간의 차이는 $ \frac{\tilde{\rho}_{\text{sup}}}{(t_j - t_k)^{1 - \frac{\nu}{2}}} $로 경계지며, 여기서 $ \tilde{\rho}_{\text{sup}} $는 계수 차이의 상한 노름을 나타낸다.
• 마르코프 체인 근사의 경우, 밀도 차이는 $ |\tilde{p}^h \boxtimes H^{h,(r)} - \tilde{p}_\nu^h \boxtimes H_\nu^{h,(r)}| \triangleq \text{error} \triangleq (r+1)\rho_{\nu,\nu,\text{sup}} \frac{\big\bracevert (1 \bigvee T^{(1-\nu)/2})c_1 \big\bracevert^{r+1} [\text{B}(\nu/2)]^r}{\text{B}(1 + r\nu/2)} \frac{c^d}{(t_j - t_i)^{d/2}} Q_{M-(d+5+\nu)}\big( \frac{y-x}{(t_j - t_i)^{1/2}/c} \big) (t_j - t_i)^{r\nu/2} $로 표현되며, 이는 $ (t_j - t_i)^{r\nu/2} $ 에 대해 명시적인 수렴을 보여준다.
• 다항 감쇠 성질을 가진 밀도의 컨볼루션 분석 결과, $ I_{t_k}^1(t_i,t_j,x,y) \triangleq \big| \big( q_m(t_k - t_i, \bullet - x) * q_m(t_j - t_k, \bullet - y) \big)(y) \big| \triangleq \text{convolution} \triangleq \bar{c} q_m(t_j - t_i, y - x) $로, $ i < k < j $ 에 대해 일관되게 성립하며, 이는 $ m > d $ 조건 하에서 성립한다.
• 파라메트릭스 전개에서 농도 손실을 통제하기 위해 영역을 분할하고 $ |w|^3 $ 및 $ |w|^4 $ 항을 통합함으로써, 적분 가능성 확보를 위해 $ m > d+5+\nu $ 가 필요하다는 것을 밝혀냈다.
• 밀도 차이의 수렴 속도는 $ O((t_j - t_i)^{r\nu/2}) $이며, 이는 가정된 헬더 정규성 $ \nu $ 하에서 최적이며, 이 경계는 공간과 시간에 대해 일관되다.
• 결과는 $ \rho_{\nu,\nu,\text{sup}} $의 관점에서 계수 변화에 대해 강건하며, 이는 $ (b, \nu) $와 $ (b_\nu, \nu_\nu) $ 간의 차이의 상한 노름을 측정하며, 이 경계는 이 변화에 대해 선형적으로 스케일링된다.