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

波動方程式の leapfrog 法の収束

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

u∈C4([0,1]×[0,T])u \in C^4([0,1]\times[0,T]) が utt=c2uxxu_{tt} = c^2u_{xx}、u(0,t)=u(1,t)=0u(0,t) = u(1,t) = 0、u(x,0)=g(x)u(x,0) = g(x)、ut(x,0)=g1(x)u_t(x,0) = g_1(x) を満たすとする。h=1/Mh = 1/M、ν=cΔt/h≤ν0<1\nu = c\Delta t/h \le \nu_0 \lt 1 とし、Uj0=g(xj)U^0_j = g(x_j)、

Uj1=g(xj)+Δt g1(xj)+c2Δt22⋅g(xj+1)−2g(xj)+g(xj−1)h2U^1_j = g(x_j) + \Delta t\,g_1(x_j) + \frac{c^2\Delta t^2}{2}\cdot\frac{g(x_{j+1}) - 2g(x_j) + g(x_{j-1})}{h^2}

を出発値とする leapfrog・中心差分の解を UnU^n とする。Bt=max⁡∣∂t4u∣B_t = \max|\partial_t^4u|、Bx=max⁡∣∂x4u∣B_x = \max|\partial_x^4u|、B3=max⁡∣∂t3u∣B_3 = \max|\partial_t^3u| とおくと、tn≤Tt_n \le T である nn について

∥u(⋅,tn)−Un∥h≤tn1−ν02(c2Δt h224Bx+Δt26B3)+tn21−ν02(Δt212Bt+c2h212Bx)\|u(\cdot, t_n) - U^n\|_h \le \frac{t_n}{\sqrt{1-\nu_0^2}}\Bigl(\frac{c^2\Delta t\,h^2}{24}B_x + \frac{\Delta t^2}{6}B_3\Bigr) + \frac{t_n^2}{1-\nu_0^2}\Bigl(\frac{\Delta t^2}{12}B_t + \frac{c^2h^2}{12}B_x\Bigr)

が成り立つ(左辺の u(⋅,tn)u(\cdot,t_n) は内部の格子点での値のベクトル)。右辺は O(Δt2+h2)O(\Delta t^2 + h^2) である。

また、格子 hZh\mathbb{Z} の上で cΔt=hc\Delta t = h(ν=1\nu = 1)のとき、u∈C2(R×[0,∞))u \in C^2(\mathbb{R}\times[0,\infty)) が utt=c2uxxu_{tt} = c^2u_{xx} を満たし、Uj0=u(xj,0)U^0_j = u(x_j, 0)、Uj1=u(xj,Δt)U^1_j = u(x_j, \Delta t) ならば、すべての n≥0n \ge 0 と jj で Ujn=u(xj,tn)U^n_j = u(x_j, t_n) である。

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

閉包の状態

✓ 公理まで証明済み

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

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

前提をたどる 証明が使う定理・定義(直接 4 件)

▸ を開くと、その項目の証明が使う前提が出ます。深い先は各項目のページでたどれます。

これを使っている定理

証明で明示的に参照している定理は、まだありません。