天津理工大学数值计算第一章 戴敏.ppt

计算机科学与工程系 计算机科学与工程系 数值计算方法 第一章 数值计算引论 1.1 数值计算方法 数值计算方法是应用数学的一个分支，又称为数值分析或计算方法，它是研究用数字计算机求解各种数学问题的数值方法及其理论的一门学科，是程序设计和对数值结果进行分析的依据和基础。 1.1 数值计算方法 利用电子计算机解题的一般步骤： 提出实际问题 构造数学模型——将实际问题归纳为明确的数学问题 选择计算方法——对数学问题，选择运算简单、工作量节省、并能保证精确度要求的计算方法，确定计算步骤 计算过程的程序设计——用一定的程序设计语言将计算过程编成程序 将计算程序和原始数据输入，上机计算，最后计算机输出计算结果 1.1 数值计算方法 数值计算方法（数值算法）：把对数学问题的解法归结为有加、减、乘、除等基本运算，并对运算顺序有完整而准确的描述的算法 数值算法具有的特点 面向计算机 有可靠的理论分析 要有好的计算复杂性 1.1 数值计算方法 总之，对于给定的数学模型所提出的可行、有效的算法应该是符合计算机的要求，在理论上收敛、稳定，在实际计算中精度高、计算复杂性小、能通过实验验证的数值方法 1.2 算法 算法的概念 人工手算过程 计算机运算过程：先将所制定的解题方案告诉计算机，令计算机按照人们所规定的计算顺序去自动执行 程序设计：用机器所能接受的“语言”来描述解题步骤 算法：指解题方案的准确而完整的描述 算法的描述：日常语言，数学语言，算法语言，算法框图 1.2 算法 举例：求解二元一次方程组 首先判断 是否为0 1.2 算法 算法框图 1.2 算法 算法的优劣 计算量的大小是衡量算法优劣的一个重要标准 例如：行列式解法的克莱姆法则，解一个n阶方程组，要算n+1个n阶行列式的值，总共需要n!(n-1)(n+1)次乘法。如果求解一个20阶左右的方程组，需要工作很多年才能完成 尽量节约存储量，也是设计算法时需要考虑的一个因素 尽量利用原有的工作单元进行累算。算式采用“动态”形式，如将y=x2表示为 ,即将单元x的内容自乘，结果仍存放在单元x中 1.2 算法 设计算法所要考虑的另一个因素就是逻辑结构的复杂性问题 算法的递推性 算法的递推化：将一个复杂的计算过程归结为简单过程的多次重复 例如：给定x的n次多项式 1.2 算法 令 上式可以表示为递推的公式： 初值为： 1.2 算法 多项式的秦九韶计算方法 将多项式按降幂的次序排列为： 将多项式表示为嵌套形式： 设用 表示第k层的值，那么第k层的结果 显然等于第k-1层的结果 乘上x再加上系数 1.2 算法 初值为： 1.2 算法 1.2 算法 例：用秦九韶方法求多项式 在x=-0.2的值 1.2 算法 1.3 误差的来源 数可以分为两类： 精确数（准确数、真值）：精确的反映实际情况的数。如教室里面有42名学生 近似数（某准确数的近似值）：只能近似的反映实际情况。如测量得到桌子的长度为950mm，这个测量值是不能精确的反映桌子实际长度的近似值 误差：将一个数的准确值与其近似值之差称为误差 1.3 误差的来源 误差产生的原因 模型误差（描述误差） 观测误差（测量误差） 截断误差（方法误差） 无穷过程用有限过程近似引起的误差 例如：指数函数 可展开成下列幂级数形式 但在实际计算中，只能截取有限项求出 1.3 误差的来源 根据泰勒余项定理， 关于 的截断误差是 例：假定|x|≤1，用泰勒多项式计算 的值，设要求截断误差不超过0.005，问公式该取多少项？ 舍入误差（计算误差） 无论用计算机、计算器计算还是笔算，都只能用有限位小数来代替无穷小数或用位数较少的小数来代替位数较多的有限小数，产生舍入误差 例如：л＝3.1415926…， 1.4 近似数的误差表示法 绝对误差 设x*是准确值x的一个近似值，则 称为近似值x*的绝对误差，简称误差。 定义 如果 那么 叫做近似数x*的绝对误差限，用它反映近似数的精度。 1.4 近似数的误差表示法 通常我们用x =x*+ε*来表示近似数x*的精确度或准确值所在的范围，ε*应该取得尽可能小 例x＝4.3762816…，取近似数x*＝4.376，则x-x*=0.0002816…，有 |x-x*|=00003=0.3×10－3 同样 |x-x*|=000