《数值计算方法》PPT课件.ppt

《《数值计算方法》PPT课件.ppt》由会员分享，可在线阅读，更多相关《《数值计算方法》PPT课件.ppt（30页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、&教材教材 (Text Book)数值计算方法数值计算方法 郑慧娆等郑慧娆等 编著编著（武汉大学出版（武汉大学出版社）社）&参考书目参考书目(Reference)Numerical Analysis:Mathematics of Scientific Computing (Third Edition)数值分析数值分析 （英文版（英文版 第第3版版 ）David Kincaid&Ward Cheney（机械工业出版社机械工业出版社)Numerical Analysis (Seventh Edition)数值分析数值分析 （第七版（第七版 影印版）影印版）Richard L.Burden&J.Do

2、uglas Faires （高等教育出版社）高等教育出版社）学习方法学习方法考试方法考试方法1.闭卷考试占闭卷考试占70%.10%.3.上机实验占上机实验占15%，创新成绩，创新成绩5%.学习和了解科学计算的桥梁学习和了解科学计算的桥梁数值分析数值分析 能够做什么？IntroductionIntroduction 研究使用计算机求解各种数学问题的研究使用计算机求解各种数学问题的数值方法（近似方法），对求得的解的数值方法（近似方法），对求得的解的精度进行评估，以及如何在计算机上实精度进行评估，以及如何在计算机上实现求解等现求解等计算机解决实际问题的步骤计算机解决实际问题的步骤建立数学模型建立数学

3、模型选择数值方法选择数值方法编写程序编写程序上机计算上机计算举例1。求下列方程的根或零点：（第三章的内容：非线性方程的数值解法）Can you solveCan you solve举例2。怎么求解下列积分？（第七章的内容：数值积分）数值分析的特点1。近似：由此产生“误差”在计算数学和应用数学中一个有趣的问题：在计算数学和应用数学中一个有趣的问题：什么是什么是零零？原点附近原点附近在纯数学中，认为此矩阵为满秩矩阵但在计算数学中，它却是降秩矩阵？2。与计算机不能分离：上机实习（掌握一门语言：C语言，会用Matlab）1.2 误差误差 (Error)1 误差的背景介绍误差的背景介绍 (Introdu

4、ction)1.来源与分类来源与分类 (Source&Classification)u模型误差模型误差 (Modeling Error):从实际问题中抽象出数从实际问题中抽象出数学模型学模型 u观测误差观测误差 (Measurement Error):通过测量得到模型通过测量得到模型中参数的值中参数的值 u方法误差方法误差(截断误差截断误差 Truncation Error）:求近似解求近似解 u舍入误差舍入误差 (Roundoff Error):机器字长有限机器字长有限 误差与有效数字误差与有效数字(Error and Significant Digits)u 绝对误差绝对误差 (absol

5、ute error)其中其中 x*为精确值，为精确值，x为为x*的近似值。的近似值。例如：例如：工程上常记为工程上常记为的上限记为的上限记为 ,称为称为绝对误差限绝对误差限 (accuracy)，u 相对误差相对误差 (relative error)x 的的相对误差上限相对误差上限 定义为定义为u有效数字有效数字（significant digits）用科学计数法，记用科学计数法，记 (其中其中 ）若若 （即（即 的截取按四舍五入规的截取按四舍五入规则），则称则），则称 为有为有n 位有效数字，精确到位有效数字，精确到 。例：例：问：问：有几位有效数字？请证明你的结论。有几位有效数字？请证明你

6、的结论。证明：证明：有有4 位有效数字，精确到小数点后第位有效数字，精确到小数点后第 3 位。位。1.4 向量和矩阵范数向量和矩阵范数 向量范数向量范数 (vector norms)对任意对任意定义定义1：Rn空间的空间的向量范数向量范数|,对任意对任意 满足下列条件满足下列条件常用向量范数：常用向量范数：=niixx11|vv=niixx122|vv|max|1inixx =v主要主要性质性质性质性质1:1:-x=x-x=x性质性质2 2:x-yx-yx-yx-y性质性质3 3:向量范数向量范数xx是是R Rn n上向量上向量x x的连续函数的连续函数.范数等价范数等价:设设A A 和和B

7、B是是R R上任意两种范数，若存在上任意两种范数，若存在 常数常数 C C1 1、C C2 2 0 0 使得使得 ,则称则称 A A 和和B B 等价等价。定理定理1.4.1 Rn 上上一切一切范数都等价范数都等价。定义定义2:2:设设x xk k是是R Rn n上的向量序列，上的向量序列，令令 x xk k=(x=(xk1k1,x,xk2k2,，x xknkn),k=1,2 k=1,2，.,又设又设x x*（x x1 1*，x x2 2*，x xn n*)是是R Rn n上的向量上的向量.如果如果lim xlim xkiki=x=xi i对所有的对所有的i=1,2,i=1,2,，n n成立，

8、成立，那么那么,称向量称向量x x*是向量序列是向量序列x xk k的极限的极限 ，若一个向量序列有极限若一个向量序列有极限,称这个向量序列是称这个向量序列是收敛的收敛的.对任意一种向量范数对任意一种向量范数而言，向量而言，向量序列序列x xk k收敛于向量收敛于向量x x*的充分必要条件是的充分必要条件是 矩阵范数矩阵范数 (matrix norms)定义定义3:3:对任意对任意 ,称称|为为Rm n空间的空间的矩阵矩阵范数范数,指指|满足满足(1)-(3)(1)-(3)：对任意对任意(4)|AB|A|B|若还满足若还满足(4),(4),称为相容的矩阵范数称为相容的矩阵范数例5:设设A A(

9、a(aijij)M.)M.定义定义证明证明：这样定义的非负实数不是相容的矩阵范数这样定义的非负实数不是相容的矩阵范数.证明：设从而相容性相容性（1 1）矩阵范数与）矩阵范数与矩阵矩阵范数的相容范数的相容:ABABABAB（2 2）矩阵范数与）矩阵范数与向量向量范数范数设设AMAM，AA是矩阵范数是矩阵范数,xRxRn n，xx是是向量范数向量范数.如果满足不等式如果满足不等式:AxAxA Ax x则称矩阵范数则称矩阵范数AA与向量范数与向量范数xx相容相容.常用的算子范数常用的算子范数：由向量范数由向量范数|p 导出关于矩阵导出关于矩阵 A Rn n 的的 p p 范数范数:则则（行和范数行和

10、范数）（列和范数列和范数）（谱范数谱范数 (spectral norm)）利用利用Cauchy 不等式不等式 可证（例可证（例6）。）。可以证明可以证明,对方阵对方阵 和和 有有:，(向量向量|2 2的直接推广的直接推广)FrobeniusFrobenius范数范数:(operator norm),又称为从属的矩阵范数又称为从属的矩阵范数:算子算子范数范数定理定理1.4.6对任意算子范数对任意算子范数|有有:证明：证明：由算子范数的相容性，得到由算子范数的相容性，得到将任意一个特征根将任意一个特征根 所对应的特征向量所对应的特征向量 代入代入命题命题(P26，推论推论1）若若A A对称，则有对

11、称，则有:A对称对称证明：证明：若若 是是 A 的一个特征根，则的一个特征根，则 2 必是必是 A2 的特征根。的特征根。又：对称矩阵的特征根为实数，即又：对称矩阵的特征根为实数，即 2(A)为非负实数，为非负实数，故得证。故得证。对某个对某个 A 的特征根的特征根 成立成立所以所以2-范数亦称为范数亦称为谱范数谱范数。若矩阵若矩阵 A 对某个算子范数满足对某个算子范数满足|A|1，则必有则必有.可逆；可逆；.证明：证明：若不然，则若不然，则 有非零解，即存在非零向量有非零解，即存在非零向量 使得使得 1.5 线性方程组的性态（误差分析）线性方程组的性态（误差分析）(Error Analysi

12、s for Linear system of Equations)思考思考:求解求解 时时,A 和和 的误差对解的误差对解 有何影响有何影响？设设 A 精确，精确，有误差有误差 ，得到的解为，得到的解为 ，即，即绝对误差放大因子绝对误差放大因子又又相对误差放大因子相对误差放大因子 设设 精确，精确，A有误差有误差 ，得到的解为，得到的解为 ，即，即(只要只要 A充分小，使得充分小，使得 是关键是关键的误差放大因子，称为的误差放大因子，称为A的状态数的状态数(条件数条件数)，记为记为cond(A),注注:condcond(A A)与与 所取的范数有关所取的范数有关常用条件数有：常用条件数有：cond(A)2特别地，若特别地，若 A 对称，则对称，则cond(A)1=A1 1cond(A)=A 例：例：Hilbert 阵阵cond(H2)=27cond(H3)748cond(H6)=2.9 106注：注：现在用现在用MatlabMatlab数学软件可以很方便数学软件可以很方便求矩阵的状态数求矩阵的状态数!定义定义2:2:设线性方程组的系数矩阵是非奇异的,如果condcond(A A)越大越大,就称这个方程组越病态.反之,condcond(A A)越小越小,就称这个方程组越良态.