计算物理1102.pdfVIP

Computational Physics Bo Li Department of Physics Chongqing Universtiy 内容 • 绪论 • 第一章 计算误差 • 第二章 常用数据处理方法 • 第三章 有限差分法 • 第四章 蒙特卡洛方法 绪论 • 什么是计算物理？ – 物理可分为理论物理、实验物理和计算物理 • 计算物理的特点 – 利用强有力的计算机工具 – 验证理论物理给出的结论 – 模拟现阶段不能达到的实验条件 • 计算物理的研究手段 – 数值计算、数值模拟 – 确定性方法与随机方法结合 第一章 计算误差 • 误差的种类 – 算法近似误差、物理模型误差、观测误差 （不讨论） – 舍入误差 round-off errors • e ≈2.7182818284590452353602874713527 • π ≈3.1415926535897932384626433832795 – 截断误差 truncation errors • Taylor expansions f 0 f 0 2 f n 0 n f x =f 0 x  x ...  x ... 1 ! 2 ! n ! f n1 y  n f x =P n x  y 0y x  n 1 ! 第一章 计算误差 • 误差的衡量 – 绝对误差与相对误差 • 假设 x 为精确值， x* 为近似值，则 • 绝对误差 e, e = x* - x ( 简称误差， e 可正可负 ) • 绝对误差限 ε, x = x* ±ε ( ε 为正 ) • 为了比较绝对误差限，引入相对误差限 εr , εr=ε/x*≈ε/x • 设 x=2.18 是由精确值经过四舍五入得到的，那么它的绝对误 差限和相对误差限是多少？ – 有效数字 • 若 x 绝对误差限是它某一数位的半个单位，并且 x 从左起第一 个非零数字到该位共 n 位数字，则称 x 的有效数字有 n 位 • 科学计算法与有效数字 x =±0. a a a ...a ∗10m 1 2 3 k 第一章 计算误差 • 数值运算中的几个典型例子 – 选用数值稳定的算法 1 • 计算积分 20 x −1 I =∫x e dx 0 • 可以利用递推公式 A =1−n A