微分方程数值方法习题2.docVIP

微分方程数值方法 常微分方程初值问题习题一 对初值问题 分别写出Euler 法和改进的Euler 法的近似解 的表达式，并求它们与真解 的差 . 取步长 ，分别用Euler 法和改进的Euler 法求下列初值问题的解，并与真解相比较. （1） 真解 ； （2） 真解 ； （3） 真解 . 用Euler 法计算 在 ，0.2 的近似值. 取步长 ，用四阶 Runge-Kutta 法解 并与真解 相比较. 设 关于 满足Lipschitz 条件，证明N 级 Runge-Kutta 法中的增量函数 关于 也满足 Lipschitz 条件. 对初值问题 写出四阶 Taylor 级数法和四阶 Runge-Kutta 法的计算公式，它们是否相同. 证明改进的 Euler 法的绝对稳定区间是（-2，0）. 证明：当 满足 时，四阶 Runge-Kutta 法绝对稳定. 用 Tayor 展开确定下面多步法中的系数，使其阶尽可能高，并求局部截断误差的主项. （1）； （2）； （3）. 10. 对初值问题 确定求解公式 中的系数 ，， 与局部截断误差主项. 求公式 的阶和局部截断误差. 取步长 ，用四阶Adams 方法的的预测-校正公式（72）解初值问题 并与真解 相比较. 讨论下列公式的相容性、稳定性、和绝对稳定性. （1）； （2）. 14. 求 ，使线性多步法 是相容的和稳定的. 证明三阶Adams 内插公式 的绝对稳定区间是（-6,0）. 证明中点公式（90）是A-稳定的. 求下列方程组的刚性比. 若用四阶 Adams 内插公式求解时，最 大步长应小于多少？ （1） （2） 18. 把下列高阶方程化为一届方程组，并写出它们的Euler 公式和四阶 Runge-Kutta 公式. （1） （2） 椭圆型方程的差分法习题三 用五点差分格式解下列椭圆型方程边值问题. （1） 取 （真解）； （2） 取（真解）； （3） 取 ； （4） 取 （真解）； （5） 取 （真解）. 2. 证明公式（34）. 3. 证明：对矩形网格，用积分守恒形式（32）导出的差分方程与五点差分格式（8）相同. 4. 证明三角形网格的差分方程（33）满足条件（36）. 5. 证明：若 ，则差分方程 的解满足 ， 其中由式（35）定义. 6. 对椭圆型方程 ， 设，,，证明： 五点差分方程 的截断误差为 ； 差分方程的系数满足条件（36）. 抛物型方程的差分法习题四 1.（上机题）对下列定解问题 取分别用古典显格式、古典隐格式和Crank-Nicolson格式计算 时的近似解，并与精确解 比较. （上机题）对定解问题 取 用双层加权平均格式，分别取计算时的近似值，并与精确解比较. 对于扩散方程 试求 Du Fort-Frankel 格式 的截断误差. 试求加权三层差分格式 的截断误差，并证明当时，截断误差的阶最高. 试求双向加权对称格式 的截断误差. 试构造变系数抛物型方程 的一个二阶精度差分格式. 证明习题3（2）、（3）中给出的两种差分格式均是绝对稳定的. 证明求解习题3中的扩散方程的Saul’ev 格式(1957) 是绝对稳定的（）. 对于二维扩散方程 试证明向后差分格式 和 Crank-Nicolson 格式 均绝对稳定（）. 8.（上机题）分别用显隐交替的格式1和格式2（跳点格式）计算习题1中的定解问题（取习题1中相同的步长），并与习题1的计算结果进行比较.