Schrödinger 方程的有限差分方法(笔记)

Schrödinger 方程的有限差分方法(笔记)

复值函数预备知识

$$\begin{array}{rl}& u=u_1+u_{2}i\\ & \bar{u}=u_1-u_{2}i\\ & u\cdot \bar{u}=\bar{u}\cdot u=u_1^2+u_2^2=|u{|}^2\\ & |u{|}_1^2={\int }_L{\left|u_x\right|}^{2}dx\\ & \Vert u{\Vert }^2=(u,u)={\int }_{L}u\cdot \bar{u}dx={\int }_L|u{|}^{2}dx.\end{array}$$

经典不等式

$H\ddot{o}lder$ 不等式(离散形式)
$$\sum _{i=1}^n{a}_i{b}_i⩽{\left(\sum _{i=1}^n{a}_i^p\right)}^{\frac{1}{p}}{\left(\sum _{i=1}^n{b}_i^q\right)}^{\frac{1}{q}}$$
$H\ddot{o}lder$ 不等式(连续形式)
$${\int }_{\Omega }f(x)\cdot g(x)dx⩽{\left({\int }_{\Omega }|f(x){|}^{p}dx\right)}^{\frac{1}{p}}\cdot {\left({\int }_{\Omega }|g(x){|}^{q}dx\right)}^{\frac{1}{q}}$$
柯西-斯瓦茨不等式
$$|⟨x,y⟩|⩽\Vert x\Vert \cdot \Vert y\Vert$$
其中 $⟨x,y⟩$ 为 $x$ 与 $y$ 的内积, $\mid \mid \cdot \mid \mid$ 为内积导出范数.

引理

通过泰勒展开式建立差分格式的过程中, 不可避免的需要讨论余项 $R_j^k$ 的阶. 有时为了完成收敛性的证明, 还需要讨论余项的一阶差商 ${\delta }_t{R}_j^{k+\frac{1}{2}}$ 的阶.

例如: 在点 $\left(x_j,t_{k+\frac{1}{2}}\right)$ 处考虑控制方程, 有
$$\begin{array}{c}\mathrm{i}{u}_t\left(x_j,t_{k+\frac{1}{2}}\right)+u_{xx}\left(x_j,t_{k+\frac{1}{2}}\right)+q{\left|u\left(x_j,t_{k+\frac{1}{2}}\right)\right|}^{2}u\left(x_j,t_{k+\frac{1}{2}}\right)=0,\\ 1⩽j⩽m-1,0⩽k⩽n-1.\end{array}$$
应用 Taylor 展开式及微分公式, 有
$$\begin{array}{r}\mathrm{i}{\delta }_t{U}_j^{k+\frac{1}{2}}+{\delta }_x^2{U}_j^{k+\frac{1}{2}}+\frac{q}{2}\left({\left|U_j^k\right|}^2+{\left|U_j^{k+1}\right|}^2\right)U_j^{k+\frac{1}{2}}=R_j^k,\\ 1⩽j⩽m-1,0⩽k⩽n-1,\end{array}$$
且存在正常数 $c_1$ 使得
$$\begin{array}{rl}& \left|R_j^k\right|⩽c_1\left({\tau }^2+h^2\right),1⩽j⩽m-1,0⩽k⩽n-1,\\ & \left|{\delta }_t{R}_j^{k+\frac{1}{2}}\right|⩽c_1\left({\tau }^2+h^2\right),1⩽j⩽m-1,0⩽k⩽n-2.\end{array}$$

要得到余项一阶差商的阶, 需要用到带积分余项的 Taylor 展开式.

---

出于简明及通用性的考虑, 此处仅讨论空间方向 $n$ 阶泰勒展开式的余项.

证明

首先给出在 $x_0$ 处的积分型余项、 拉格朗日型余项和柯西型余项, 它们分别是
$$\begin{array}{rl}{R}_n(x)& =\frac{1}{n!}{\int }_{x_0}^x{f}^{(n+1)}(s)(x-s{)}^{n}ds\\ R_n(x)& =\frac{1}{(n+1)!}{f}^{(n+1)}(\xi )(x-x_0{)}^{n+1},\xi \text{在}{x}_0\text{与}x\text{之间}\\ R_n(x)& =\frac{1}{n!}{f}^{(n+1)}(\theta x)(1-\theta {)}^n(x-x_0{)}^{n+1},0⩽\theta ⩽1\end{array}$$
于是, 将 $f\left(x_j+h\right)$ 在点 $x_j$ 处展开同时$f\left(x_j+2h\right)$ 在点 $x_j+h$ 处展开的积分型余项分别为
$$\begin{array}{rl}{R}_j& =\frac{1}{n!}{\int }_{x_j}^{x_j+h}{f}^{(n+1)}\left(s^{\prime }\right){\left(x_j+h-s^{\prime }\right)}^{n}d{s}^{\prime }.\text{令}{s}^{\prime }=sh+x_j.\\ R_{j+1}& =\frac{1}{n!}{\int }_{x_j+h}^{x_j+2h}{f}^{(n+1)}\left(s^{\prime }\right){\left(x_j+2h-s^{\prime }\right)}^{n}d{s}^{\prime }.\text{令}{s}^{\prime }=sh+x_j+h.\end{array}$$
可得
$$\begin{array}{rl}{R}_j& =\frac{1}{n!}{\int }_0^1{f}^{(n+1)}\left(sh+x_j\right)(1-s{)}^n{h}^{n+1}ds.\\ R_{j+1}& =\frac{1}{n!}{\int }_0^1{f}^{(n+1)}\left(sh+x_j+h\right)(1-s{)}^n{h}^{n+1}ds.\end{array}$$
于是有
$$R_{i+1}-R_j=\frac{1}{n!}{\int }_0^1\left[f^{(n+1)}\left(sh+x_j+h\right)-f^{(n+1)}\left(sh+x_j\right)\right](1-s{)}^n{h}^{n+1}ds.\text{\hspace{1em}(1)}$$
根据泰勒展式, 将 $f^{(n+1)}\left(sh+x_j+h\right)$ 在点 $sh+x_j$ 处展开得
$$f^{(n+1)}\left(sh+x_j+h\right)=f^{(n+1)}\left(sh+x_j\right)+h{f}^{(n+2)}\left(sh+x_j\right)+o\left(h^2\right)$$
移项得
$$f^{(n+1)}\left(sh+x_j+h\right)-f^{(n+1)}\left(sh+x_j\right)=h{f}^{(n+2)}\left(sh+x_j\right)+o\left(h^2\right)\text{\hspace{1em}(2)}$$
将 $(2)$ 带入 $(1)$ 得到
$$\begin{array}{rl}{R}_{i+1}-R_j& =\frac{1}{n!}{\int }_0^1\left[h{f}^{(n+2)}\left(sh+x_j\right)+o\left(h^2\right)\right](1-s{)}^n{h}^{n+1}ds\\ & =\frac{1}{n!}{\int }_0^{1}h{f}^{(n+2)}\left(sh+x_j\right)(1-s{)}^n{h}^{n+1}ds+\frac{1}{n!}{\int }_0^{1}o\left(h^2\right)(1-s{)}^n{h}^{n+1}ds\\ & =o\left(h^{n+2}\right)+o\left(h^{n+3}\right)\\ & =o\left(h^{n+2}\right)\end{array}.$$
于是
$${\delta }_x{R}_{j+\frac{1}{2}}=\frac{R_{j+1}-R_j}{h}=\frac{o\left(h^{n+2}\right)}{h}=o\left(h^{n+1}\right)$$
综上所述, 余项与其一阶差商同阶.

证明完毕.