数值计算方法课件--第一章--绪论.ppt

《数值计算方法课件--第一章--绪论.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、玉溪师范学院理学院物理系数值计算方法NumericalComputationalMethod提问：数值计算方法是做什么用的？研究对象：数值问题——有限个输入数据（问题的自变量、原始数据）与有限个输出数据（待求解数据）之间函数关系的一个明确无歧义的描述。如一阶微分方程初值问题求函数解析表达式数学问题求函数在某些点的近似函数值数值问题程序设计上机计算设计高效、可靠的数值方法数值问题求解近似结果输出重点讨论数值问题的来源：实际问题建立数学模型数值问题数值方法的设计原则收敛性：方法的可行性稳定性：初始数据等产生的误差对结果的影响便于编程实现：逻辑复杂度要小计算量要小：时间复杂度要小，运行时间要

2、短存贮量要尽量小：空间复杂度要小可靠性分析计算复杂性误差估计：运算结果不能产生太大的偏差且能够控制误差物理学中：Origin+fortran七八十年代物理学家Origin+C或C++现代物理学家数学中：MatlabMathematicsFortranCC++物理问题中：§1误差/*Error*/一、误差的来源与分类/*Source&Classification*/1、从实际问题中抽象出数学模型——模型误差/*ModelingError*/2、通过观测得到模型中某些参数（或物理量）的值——观测误差/*MeasurementError*/3、数学模型与数值算法之间的误差——方法误差(截断误

3、差/*TruncationError*/)4、由于机器字长有限，原始数据和计算过程会产生新的误差——舍入误差/*RoundoffError*/二、误差分析的基本概念/*BasicConcepts*/设为真值（精确值），为的一个近似值称为近似值的绝对误差，简称误差。注：误差可正可负，常常是无限位的绝对误差限/*accuracy*/——绝对值的上界如：绝对误差还不能完全表示近似值的好坏（绝对误差/*absoluteerror*/）近似值的误差与准确值的比值：称为近似值的相对误差，记作注：实际计算时，相对误差通常取因为（相对误差/*relativeerror*/）相对误差也可正

4、可负,常常是无限位的（有效数字/*SignificantDigits*/）相对误差限——相对误差的绝对值的上界/*relativeaccuracy*/如：3位6位若近似值与准确值的误差绝对值不超过某一位的半个单位，该位到的第一位非零数字共有位，则称有位有效数字有效数字(另外一种定义形式)注：0.2300有4位有效数字，而00023只有2位有效数字。12300如果写成0.123105，则表示最多只有3位有效数字。数字末尾的0不可随意省去！用科学计数法，记其中，若（即的截取按四舍五入规则），则称为有位有效数字，精确到。例1：问：有几位有效数字？请证明你的结论。证明：有位有效数字，精确

5、到小数点后第位。43规格化形式注：若的每一位都是有效数字，则称是有效数特别，经“四舍五入”得到的数均为有效数将的近似值表示为，若是有效数字，则相对误差不超过；反之，若已知相对误差，且有，则必为有效数字。例2：为使的相对误差小于0.001%,至少应取几位有效数字？解：假设*取到n位有效数字，则其相对误差上限为要保证其相对误差小于0.001%，只要保证其上限满足解不等式可得n>5.69，即n=6，应取*=3.14159。例3计算下列多项式的值为已知数据分析：输入数据为，输出数据为，若直接由算出，再乘相应的系数并相加，则要做次乘法和次加法，占用个存储单元。秦九韶方法，也称为Horne

6、r算法用递推公式表示为只用次乘法和次加法，并占用个存储单元三、数值算法及稳定性/*NumericalAlgorithmandStability*/11/e解法：将作Taylor展开后再积分S4R4/*Remainder*/取则称为截断误差/*TruncationError*/

7、舍入误差/*RoundoffError*/

8、=0.746824……由截去部分/*excludedterms*/引起由留下部分/*includedterms*/引起例4近似计算一个算法如果输入数据有扰动（即误差），而计算过程中舍入误差不增长,则称此算法是数值稳定的，否则此算法就称为不稳定的。（数值稳定性/*Nume

9、ricalStability*/）对数学问题本身如果输入数据有微小扰动，引起输出数据（即问题真解）的很大扰动，这就是病态问题。（病态问题/*ill-posedproblem*/）它是数学问题本身性质所决定的，与算法无关，也就是说对病态问题，用任何算法（或方法）直接计算都将产生不稳定性。此公式精确成立记为则初始误差Whathappened?!例5计算公式一：考察第n步的误差我们有责任改变。造成这种情况的是不稳定的算法/*unstablealgo