计算方法-1引论.pptVIP

1 简介 参考书目 《计算方法与实习》第四版，袁慰平，东南大学出版社，2005 《计算方法》，李桂成，电子工业出版社，2005 《数值计算方法》，吕同富，清华大学出版社，2008 湍流的基本方程 数值解与解析解区别：解析解是精确解；数值解是近似解 * * * * 数值解与解析解区别：解析解是精确解；数值解是近似解 * 数值解与解析解区别：解析解是精确解；数值解是近似解 * * 以班上某位同学的体重为例，介绍误差及误差限 * * * F含有 1+2（ β-1） βt-1（U-L+1） 含义： 1表示机器数零； 2表示正负； 3（ β-1） βt-1表示ai的变化，原本有βt种变化，但因为a1是不为零的数，因此又要去除βt-1，于是有（βt -βt-1） 4 阶码的变化为U-L+1 * * * 数据误差的影响 对函数来说，自变量的误差必然引起函数值的误差，函数越复杂，这种数据误差的影响也随之复杂，通常用泰勒级数展开的方法来估计。 对二元函数 对函数 y=f(x1,x2),设x1，x2是近似值，x1*,x2*是准确值，函数准确值为：y*=f(x1*,x2*)，于是，函数y的绝对误差为： e(y) = y* - y = f(x1*,x2*) - f(x1,x2) 将f(x1*,x2*)在（x1,x2）处作泰勒展开，并取一阶泰勒多项式，得e(y)的近似表达式为： 函数的相对误差 几个常用的计算函数和、差、积、商的公式 e(x1+x2) ≈ e(x1) + e(x2) e(x1- x2) ≈ e(x1) - e(x2) e(x1 · x2) ≈ x2e(x1) + x1e(x2) e(x1 / x2) ≈ e(x1)/x2 - x1e(x2) /x22 er(x1+x2) ≈ er(x1) x1/(x1+x2)+ er(x2) x2/(x1+x2) er (x1- x2) ≈ er(x1) x1/(x1-x2)- er(x2) x2/(x1-x2) er (x1 · x2) ≈ er(x1) + er(x2) er (x1 / x2) ≈ er(x1) - er(x2) 例5 已测得某物体行程s*的近似值s=800m，所需时间t*的近似值t=35s。若已知|t* -t |≤0.5s ，试求平均速度v的绝对误差限和相对误差限。 解：公式 v=s/t的误差估计式为： 绝对误差 e(v) = e (s/t) ≈ e(s)/t – se(t)/t2 相对误差 er(v) = er (s/t) ≈ er (s) – er (t) | e(v) | ≈| e(s)/t – se(t)/t2| ≤| e(s)/t |+ | se(t)/s2|≤ 0.5/35+800*0.05/352 ≈0.0469≤ 0.05 | er (v) | ≈| e(s) – er (t)| ≤| er (s) |+ | er (t)|≤ 0.5/800+0.05/35 ≈0.00205 5 误差危害的防止 例7 建立积分 n=0,1,2,……20 的递推关系式，并研究它的误差传递 判据：若第n+1步的误差en+1与第n步的误差en满足 则称此计算公式是绝对稳定的，否则称为不是绝对稳定的 使用数值稳定的计算公式 尽量避免两相近数相减 由 可知，有两个相近的数相减时，可能会导致较大的误差出现。 尽量避免用绝对值很大的数做乘数 由 可知，若x1或x2很大时，可能导致|e(x1x2)|很大 防止大数吃掉小数 在数值运算时，若参加运算的数数量级相差很大，而计算机位数有限，又要作对阶处理，如果不注意运算次序，就可能导致小数被当作机器零处理，影响结果可靠性 注意简化步骤，减少运算次数 作业 P15 1.3，1.4 * 数值解与解析解区别：解析解是精确解；数值解是近似解 * 数值解与解析解区别：解析解是精确解；数值解是近似解 * 提问：方程如何求解（因式分解） 板书 * * * * 模型以班上某位同学的体重为例，介绍误差及误差限 * 以班上某位同学的体重为例，介绍误差及误差限 * 身高1.85m，误差等级1cm 地球到太阳1.5亿公里，误差等级公里 * * * 例1可以请同学上讲台或直接回答试做，然后对照答案 * 例1可以请同学上讲台或直接回答试做，然后对照答案 * * * * 可以请同学上讲台或直接回答试做，然后对照答案 * * * * * 计 算 方 法 主 讲： 孙 在 Email: sunzai@cjlu.edu.cn 2010 孙在，重庆人。 1991-1998