情報科学 › 偏微分方程式の離散化 › 波動方程式と Schrödinger 方程式の時間発展

leapfrog 法の離散エネルギー

定理 4.1(「波動方程式と Schrödinger 方程式の時間発展」§4) ・ ✓公理まで証明済み

AA を実対称な正定値行列、θ=Δt2λd/4\theta = \Delta t^2\lambda_d/4 とする。実ベクトルの列 unu^n が、与えられた rnr^n について

un+1−2un+un−1=−Δt2Aun+Δt2rn(n≥1)u^{n+1} - 2u^n + u^{n-1} = -\Delta t^2Au^n + \Delta t^2r^n\qquad(n \ge 1)

を満たすとし、δn=(un+1−un)/Δt\delta^n = (u^{n+1} - u^n)/\Delta t、uˉn=(un+1+un)/2\bar u^n = (u^{n+1} + u^n)/2、

En=12∥δn∥2+12⟨un+1,Aun⟩E^n = \frac12\|\delta^n\|^2 + \frac12\langle u^{n+1}, Au^n\rangle

とおく。

  1. En−En−1=Δt2⟨rn,δn+δn−1⟩E^n - E^{n-1} = \frac{\Delta t}{2}\langle r^n, \delta^n + \delta^{n-1}\rangle。とくに rn=0r^n = 0 なら En=E0E^n = E^0(すべての nn)である。
  2. En=12∥δn∥2−Δt28⟨Aδn,δn⟩+12⟨Auˉn,uˉn⟩E^n = \frac12\|\delta^n\|^2 - \frac{\Delta t^2}{8}\langle A\delta^n, \delta^n\rangle + \frac12\langle A\bar u^n, \bar u^n\rangle で、 1−θ2∥δn∥2+12⟨Auˉn,uˉn⟩≤En≤12∥δn∥2+12⟨Auˉn,uˉn⟩.\frac{1-\theta}{2}\|\delta^n\|^2 + \frac12\langle A\bar u^n, \bar u^n\rangle \le E^n \le \frac12\|\delta^n\|^2 + \frac12\langle A\bar u^n, \bar u^n\rangle.
  3. θ<1\theta \lt 1 のとき En≤E0+Δt2(1−θ)∑m=1n∥rm∥\sqrt{E^n} \le \sqrt{E^0} + \dfrac{\Delta t}{\sqrt{2(1-\theta)}}\displaystyle\sum_{m=1}^{n}\|r^m\|。rn=0r^n = 0 なら、すべての nn で ∥δn∥2≤2E01−θ,∥un∥≤2E0λ1+Δt22E01−θ.\|\delta^n\|^2 \le \frac{2E^0}{1-\theta},\qquad \|u^n\| \le \sqrt{\frac{2E^0}{\lambda_1}} + \frac{\Delta t}{2}\sqrt{\frac{2E^0}{1-\theta}}.

本文で読む(証明と前後の説明)→

閉包の状態

✓ 公理まで証明済み

証明が使う定理を 前提 としてたどると、Rikai の中の証明だけで定義・公理(物理では基本仮定)までつながっています。

印の意味: ✓公理まで証明済み ・ △証明あり(閉包未完)・ ※主張のみ。テーマごとの集計は 閉包の状況。

前提をたどる 証明が使う定理・定義

証明は Rikai の他の定理を参照していません(定義と論理だけで示されています)。

これを使っている定理 証明で参照している 1 件

関連するノック 1 問