3.2 五点差分格式 3.2.1 五点差分格式的建立 (1) 建立差分格式 将区间 \([a,b]\) 做 \(m\) 等分，记 \[h_1 = \frac{b-a}{m}, \quad x_i = a + ih_1, \quad i=0,\cdots,m \] 将区间 \([c ...

4.5 Crank-Nicolson 格式 本节对于定解问题 \((3.1.1) \sim (3.1.3)\) 建立一个具有 \(O(\tau^2 + h^2)\) 精度的无条件稳定的差分格式。 注意，对各个符号取上标 \(k+\frac{1}{2}\) 和取下标 \(k+\frac ...

实验目的 用Matlab实现欧拉法、后退欧拉法、梯形方法和改进欧拉公式 实验要求 1. 给出欧拉法、后退欧拉法、梯形方法和改进欧拉公式算法 2. 用Matlab实现欧拉法、后退欧拉法、梯形方法 ...

MATLAB常微分方程的数值解法 作者：凯鲁嘎吉 - 博客园http://www.cnblogs.com/kailugaji/ 一、实验目的 科学技术中常常要求解常微分方程的定解问题，所谓数值解法就是求未知函数在一系列离散点处的近似值。 二、实验原理 三、实验程序 1. 尤拉 ...

1.2 Euler 方法及其改进方法 1.2.1 Euler 方法 用 \(f(x_n, y_n)\) 代替式 \((1.2)\) 中的 \(\varphi_n\)，得到差分方程初值问题： \[\left\{ \begin{align*} & y_{n+1} = y_{n ...

10 常微分方程初值问题的数值解法 10.1 引言 包含自变量、未知函数以及未知函数导数或微分的方程称为微分方程。在微分方程中，如果自变量的个数只有一个，就称为常微分方程；如果自变量个数两个及以上，就称为偏微分方程。微分方程中出现的未知函数最高阶导数的阶称为微分方程的阶。如果未知函数\(y ...

1.5 相容性、收敛性与稳定性 1.5.1 相容性与收敛性 定义相容性。（非数学性质严格） 定义 1.5.1 相容性 当步长 \(h \to 0\) 时，差分方程是否无限逼近微分方程。 定义收敛性。（非数学性质严格） 定义 1.5.2 收敛性 ...

常微分方程的差分方法重点回顾： 差分方法是一类重要的数值解法，这类方法是要寻求一系列离散结点上的近似解h为步长。一般来说，假定h为定数。 能求解的常微分方程的条件。A.两个方程B.满足李普希兹条件C.f(x,y)适当光滑。这样可以保证解存在且唯一。 数值解法的第一步是设法消除其导数 ...