《数值分析-第一章-学习小结》由会员分享,可在线阅读,更多相关《数值分析-第一章-学习小结(9页珍藏版)》请在金锄头文库上搜索。
1、数值分析第1章 绪论 -学习小结一、 本章学习体会 通过本章的学习,让我初窥数学的又一个新领域。数值分析这门课,与我之前所学联系紧密,区别却也很大。在本章中,我学到的是对数据误差计算,对误差的分析,以及关于向量和矩阵的范数的相关内容。误差的计算方法很多,对于不同的数据需要使用不同的方法,或直接计算,或用泰勒公式。而对于二元函数的误差计算亦有其独自的方法。无论是什么方法,其目的都是为了能够通过误差的计算,发现有效数字、计算方法等对误差的影响。而对误差的分析,则是通过对大量数据进行分析,从而选择出相对适合的算法,尽可能减少误差。如果能够找到一个好的算法,不仅能够减少计算误差,同时也可以减少计算次数
2、,提高计算效率。对于向量和矩阵的范数,我是第一次接触,而且其概念略微抽象。因此学起来较为吃力,仅仅知道它是向量与矩阵“大小”的度量。故对这部分内容的困惑也相对较多。本章的困惑主要有两方面。一方面是如何能够寻找一个可靠而高效的算法。虽然知道算法选择的原则,但对于很多未接触的问题,真正寻找一个好的算法还是很困难。另一方面困惑来源于范数,不明白范数的意义和用途究竟算什么。希望通过以后的学习能够渐渐解开自己的疑惑。二、 本章知识梳理绪论数值分析的研究对象误差知识与算法知识向量范数与矩阵范数方法的构造2.1 数值分析的研究对象研究对象求解过程的理论分析数值分析是计算数学的一个重要分支,研究各种数学问题的
3、数值解法,包括方法的构造和求解过程的理论分析。它致力于研究如何用数值计算的方法求解各种基本数学问题以及在求解过程中出现的收敛性,数值稳定性和误差估计等内容。2.2误差知识与算法知识 2.2.1误差来源误差按来源分为模型误差、观测误差、截断误差、舍入误差与传播误差五种。其中模型误差与观测误差属于建模过程中产生的误差,而截断误差、舍入误差与传播误差属于研究数值方法过程中产生的误差。2.2.2绝对误差、相对误差与有效数字1.(1)绝对误差e指的是精确值与近似值的差值。绝对误差:e=x-a绝对误差限:e,x=a(2)相对误差是指绝对误差在原数中所占的比例。相对误差:er=x-ax或er=x-aa相对误
4、差限:era=r结论:凡是经过四舍五入而得到的近似值,其绝对误差不超过该近似值末位的半个单位。(3)有效数字的定义有效数字的第一种定义:设a是x的近似值,如果a的误差绝对值不超过x的第k位小数的半个单位,即|x-a|1210-k则称近似值a准确到小数点后第k位。从小数点后的第k位数字直到最左边非零数字之间的所有数字都叫有效数字。有效数字第二种定义:设数x的近似值x*=0.x1x2xn10m其中m是整数,x1x2xn是0,1,2,9中的任意数,但x10,若x-x*1210m-k (1kn) 则x*具有k位有效数字。通过学习总结出下面几个结论:(1)若a是经过四舍五入而得到的近似值,则从它的末位数
5、字到第一位非零数字都是有效数字。(2)将任何数乘以10p(p=0,1,2,)等于移动该数的小数点,并不影响其有效数字。(3)有效数字相同的两个近似值的绝对误差不一定相同。(4)准确值被认为具有无穷多位有效数字。从有效数字的定义可以知道,由准确值经过四舍五入得到的近似值,从它的末位数字到第一位非零数字都是有效数字。2.(1)相对误差与有效数字的关系:若近似数x*=0.x1x2xn10m具有n位有效数字,则其相对误差|erx*|12x1101-n。若近似数x*=0.x1x2xn10m的相对误差|erx*|5x1+110-n 则该近似数至少具有n位有效数字。结论:有效数字位数越多,相对误差越小。(2
6、)绝对误差与有效数字的关系:若a=a1a2a3ak10m 其中m是整数,aii=1,2,k 是0到9中的一个数字,a10.如果a作为数x的近似值,且a具有n位有效数字,则|ea|1210m-n若a=a1a2a3ak10m 其中m是整数,aii=1,2,k 是0到9中的一个数字,a10.如果a作为数x的近似值,如果|ea|1210m-n(nk)则a 具有n位有效数字。结论:有效数字位数越多,绝对误差越小。2.2.3误差估计的基本方法1.(1)对于一元函数: efafaea(f(a)|fa|(a)(2)二元函数:efa,bf(a,b)xea+f(a,b)yeb(3)n元函数:ef(x1*,x2*,
7、x3*,xn*)inf(x1*,x2*,x3*,xn*)xie(xi*)f(x1*,x2*,x3*,xn*)inf(x1*,x2*,x3*,xn*)xi(xi*)设u=f(a)存在足够高阶的导数,a是自变量x的近似值,则u=f(a)是u=f(x)的近似值。如果fa=fafk-1a=0,fk(a)0,且比值fk+1afk(a)不是很大,则efafk(a)k!(e(a)(k)fafk(a)k!(a)(k)2.算数运算误差:(ab)=(a)+(b)(ab)a(b)+b(a)aba(b)+b(a)b2 ,b0r(a+b)=(a)+(b)a+br(a-b)=(a)+(b)a-br(ab)r(b)+r(a
8、)rabr(b)+r(a)2.2.4算法及计算复杂性在数值计算中,要注意遵循一些原则,以保证数值稳定性。(1)能控制舍入误差的传播。(2)合理安排量级相差悬殊数间的运算次序,防止大数将小数吃掉。(3)避免两个相近的数相减。(4)避免接近零的数做除数,防止溢出。(5)简化计算步骤,尽量减少运算次数。2.3向量范数与矩阵范数2.3.1 向量范数1.向量范数满足三个条件:(1) 正定性(2) 齐次性(3) 成立三角不等式2.对于Rn中的任一向量x=x1,x2,xnT则有1- 范数(列范数) x1=inxi2- 范数(欧氏范数)x2=inxi2P-范数xp=inxip1p limpxp=x-范数x=m
9、ax1inxi3.在空间Rn中可以引进各种向量范数,且它们都满足下述向量定理:设,是Rn上的任意两种向量范数,则存在与向量x无关的数m和M(0mM),使下列关系成立。mxxMx也就是说,向量x的某一范数可以任意小(大)时,该向量的其它任意一种范数也会任意小(大)。2.3.2矩阵范数1.定义在Rnn上的实值函数称为矩阵范数,如果对于Rnn中任意的矩阵A和B,阵范数满足下列条件:(1) 非负性(2) 齐次性(3) 成立三角不等式(4) 相容性2.当一个问题中需要向量范数和矩阵范数时,向量范数和矩阵范数应该是相容的。对于给定的向量范数和矩阵范数,如果对于任一个xRn,ARnn,满足AxAx,则所给的
10、向量范数和矩阵范数是相容的 。设在Rn中给定了一种向量范数,对任意矩阵ARnn,令A=maxx=1Ax,由此定义的矩阵范数与给定的向量范数相容,将这种范数称为从属于所给定的向量范数的矩阵范数。3.设A=aijRnn,则:矩阵A的列范数 A1=max1ininaij矩阵A的谱范数 A2=maxATA矩阵的行范数A=max1jnjnaij弗罗贝尼乌斯范数AF=i,j=1naij24.设矩阵ARnn的某种范数A1,则IA为非奇异矩阵,并且当这种范数为算子范数时,还有IA-111-A 成立。三、 本章思考题问题:向量和矩阵有多种范数,如1范数、2范数、范数。而作为向量和矩阵“大小”的度量,为什么要用这么多种范数来度量,而不是专门指定一种范数?个人理解:1. 对于不同向量和矩阵,从运算等方面考虑,某一种或几种范数在计算上较为简单方便;2. 对于不同领域,某一种或者几种范数,其应用价值和使用价值更高。四、 本章测验题.设试求 知识点:关于1-范数、-范数、佛罗比尼乌斯范数的概念及其计算。解:max(2+|-5|+3,3+|-6|+5,4+|-7|+9)=20max(2+3+4,|-5|+|-6|+|-7|,3+5+9)=189