计算方法与实习(第五版)期末复习资料.docVIP

第 PAGE 8 页 共 NUMPAGES 16 页 《计算机在材料科学中的应用》习题课 第一章 误差等概念 误差来源：模型误差、观测误差、截断误差、舍入误差 绝对误差（限）：e=x*-x，|e|＝|x*-x|≤ε 相对误差（限）：er＝(x*-x)/x，|er|＝|x*-x|/|x|≤εr 有效数字：|e|≤ 防止误差的危害：避免两相近数相减，多数作乘数或小数作除数，大数“吃”小数 第二章 方程求根 根的存在及隔离 二分法：误差是 迭代法： 加速法： 牛顿迭代法： 第三章 方程组求解 消去法：高斯消去法，列主元消去法，高斯-约当法，消元因子 消元公式 回代公式 矩阵直接分解：紧凑格式 追赶法 迭代法：收敛条件 ①雅可比法迭代格式：②高斯-赛德尔法迭代格式： 第四章 插值法 插值多项式 拉格朗日插值： 插值基函数 差商： 牛顿插值公式 f(x)=f(x0)+f[x0,x1](x-x0)+f[x0,x1,x2](x-x0)(x-x1)+… +f[x0,x1,…,xn](x-x0)(x-x1)…(x-xn-1) 差分（等间距节点） 牛顿前插公式 样条插值：三次样条插值，要求光滑、连续 第五章 曲线拟合 最小二乘原理 列表计算累加和如下 j … … 1 2 … n 从正规方程组中解得拟合多项式的各个系数值ai (i=0,1,…,m)。 第六章 数值积分、微分 积分的有限过程  插值型求积公式用插值多项式代替被积函数， 从而在有两个求积节点时得到梯形公式 有三个等距求积节点时得到Simpson公式 柯特斯公式（等距节点情况）：①柯特斯系数 ②柯特斯求积公式（有五个等距求积节点时） 复化求积①复化梯形公式 ②复化Simpson公式③复化柯特斯公式 步长自适应 龙贝格求积公式 数值微分①二点公式 ②三点公式  第七章 常微分方程的数值解法 边值问题和初值问题①边值问题 ②初值问题  求解满足上述两式的近似值yi ，即在a￡x0￡x1￡…￡xn￡b上的y(xi)的近 似值yi (i=0,1,2,…,n)。通常取等距节点，即h=xi+1-xi，有xi=x0+ih (i=0,1,2,…,n)。初值问题的数值解法特点：按节点顺序依次推进，由已知的y0 , y1 , …, yi ，求出yi+1，这可以通过递推公式得到。 单步法①欧拉折线法 ②改进欧拉法 ③四阶龙格-库塔法 多步法：阿当姆斯内插、外推预测校正公式 一阶微分方程组求解 高阶微分方程求解 例1.求解超越方程 在区间〔0，1〕内的根，要求有二位有效数字。解：设f(x)= ,该函数是连续的，且f(0)=1，f(1)=e-1-10, 因f(0)×f(1) 0， f’(x)=-e-x-π/2cos（πx/2）＝-（e-x+π/2cos（πx/2））0 〔0，1〕, 单调递减，所以方程在〔0，1〕内只有1个根。 用二分法求解，先求二分次数 要求满足 1/2k+1 (1-0) 10-2/2 则解得k后整k=6，用二分法求解过程如下： k ak bk xk f(xk) 0 0 1 0.5 -0.100576 1 0 0.5 0.25 0.396117 2 0.25 0.5 0.375 0.131719 3 0.375 0.5 0.4375 0.011255 4 0.4375 0.5 0.46875 -0.045775 5 0.4375 0.46875 0.4453125 -0.003208 6 0.4375 0.4453125 0最后求得根为0.44。 例2.用牛顿迭代法求出f(x)= x41+x3+1＝0在x0＝-1附近的实根，要求迭代至 | xk+1 -xk|0.001时停止。 解：迭代是否收敛可用f(x0)f’’(x0)0判断， f(x0)＝f(-1)＝（-1）41+（-1）3+1＝-10 f’(x)=41x40+3x2, f’’(x)=41*40x39+3*2x, f’’ (x0) =41*40*（-1）39+3*2*（-1）0, 所以f(x0)f’’(x0)0，迭代是收敛的。牛顿迭代格式为 ，取x0＝-1,计算如下 k xk 0 -1 1 -0.9773 2 -0.9605 3 -0.9534 4 -0.9525 5 -0.95