数值计算和最优化(lecture1)计算方法第一章.pptVIP

要求: 1. 课前预习和课后复习,并及时完成课后作业.（掌握构造方法的原理和思想） 2. 课余将书本上的程序进行上机计算, 掌握并熟练运用Matlab进行数值计算。（掌握运用数值计算方法在计算机上解决各类工程问题） 课程成绩： 平时成绩（含考勤、作业及平时表现）30% 期末考试成绩 70% “备”则“倍” 有准备、有规划的人生更精彩！ 四、误差的传播 1、数据误差的传播 由多元函数的Taylor展开公式可得， 的绝对误差为： 相对误差为： 称为 f 的条件数，其绝对值的大小可反映函数值对数据的敏感程度 利用上面的误差估计公式，可以得到两个数的和、差、积、商的误差估计 2、舍入误差的传播 因舍入导致的相对误差限仅与计算机的字长有关，通常 称相对误差限 为计算机的相对精度。 即 在计算机中，数需首先转化为机器数，比如浮点数，在 运算器中参与运算后仍需将运算结果转化成浮点数的形 式进行存储。 由上面的讨论可以看出，为了求得满意的计算解，在选 用计算公式和设计算法时，都应注意如下普遍原则： (1) 防止大数吃小数 主要由计算机的位数引起 选用算法应遵循的原则 计算机中数的计算特点： 加法先对阶，后运算，再舍入。 乘法先运算，再舍入。 不在计算机数系中的数做四舍五入处理。 作一个有效数字为4位的连加运算 而如果将小数放在前面计算 在作连加时，为防止大数吃小数，应从小到大进行相加， 如此，精度将得到适当改善。当然也可采取别的方法。 例 (2) 作减法时应避免两个相近数相减 两个相近的数相减，会使有效数字的位数严重损失！ 例1.2.10 用四位浮点数计算 解 只有一位有效数字,有效数字大量损失，造成相对 误差扩大。 结果仍然有四位有效数字。 这说明了算法设计的重要性。在算法设计中，若可能出现两个相近数相减，则改变计算公式，如使用三角变换、有理化等等。 (3) 避免小数作除数和大数作乘数 小数作除数或大数作乘数会产生溢出错误，因而产生大的误差。 在算法设计时，要避免这类情况在计算公式中出现。此时可以 根据一些具体情况, 把某些算式改写成另一种等价的形式，如 分母有理化等。 根据误差传播的估计式 §3.算法的稳定性 如前所述，由于各种误差的存在，计算机往往只能近似地求解实际问题，因而计算时会冒风险。 一、问题的性态 如把方程组的系数 舍入成两位有效数字 它的精确解为x1 = -6.222... x2= 38.25… x3= -33.65... 例 求解线性方程组 其精确解为 x1=x2=x3=1. 若对方程组的系数和中间结果均取3位10进制有效数字，然后用Gauss消元法求解，得到计算解为： 显然，该计算解的精度较差。 同样用Gauss消元法求解方程组： 也取3位10进制有效数字，得到计算解为： 容易验证，它是方程组的精确解。 上述例子表明，数值问题计算解的精度，与数值问题本 身的性态有关。 定义1.3.1 在数值问题中，如果输出数据对输入数据的 扰动（如误差）很敏感，即若输入数据（如原始数据） 有较小的变化，会引起输出数据（如计算解）的较大变 化，称这类数值问题为病态问题或坏条件问题。非病态 问题又称为良态问题。问题输出变量的相对误差与输入 变量的相对误差的商称为问题的条件数 二、算法的稳定性与设计原则 例1.3.3 计算定积分 解 一个程序往往要进行大量的四则运算才能得出结果，每一步的运算均可能会产生舍入误差。在大量计算中,舍入误差是积累还是能控制,这与算法有关。 误差放大 5千倍! 并假设计算过程中不产生新的舍入误差 。 误差会放大 由公式 可推出： * 主讲教师： 白敏茹 副教授 Email: minru-bai@163.com Tel $1.1 引 言 第一章 三大科学方法： 理论研究、科学实验、科学计算 科学计算核心技术：数值计算技术和计算机技术 计算数学、图像处理、统计分析、计算机仿真。 科学计算包含： 一个科学计算过程主要包括如下几个环节： (1) 数学建模：将工程问题数学化 工程中的数学模型一般可分为三类： (2) 算法设计：将数学问题数值化(指标：速度和精度) 连续型（确定型） 离散型（统计型） 不确定型（随机型） 本书重点讨论 例1.1.4 求解线性方程组 求解二次方程 是数值问题 求解微分方程 不是数值问题 将其变成数值问题，即将其“离散化” “离散化”是将非数值问题的数学模型化为数值问题 的主要方法，这也是计算方法的任务之一 (3) 程序设计：将数值问题机器