资源描述:
《数值分析2012总结》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、数值分析2012年复习总结任课教师王希云第二章数值分析基本概念教学内容:1.误差与有效数字误差、误差限、相对误差、相对误差限和有效数字的定义及相互关系;误差的来源和误差的基本特性;误差的计算(估计)的基本方法。2.算法的适定性问题数值分析中的病态和不稳定性问题;病态问题和不稳定算法的实例分析。3.数值计算的几个注意问题数值计算的基本概念�误差概念和分析误差的定义:设x是精确值,p是近似值,则定义两者之差是绝对误差:∆=x−pa由于精确值一般是未知的,因而Δ不能求出来,但可以根据测量误差或计算情况估计它的上限
2、-xp
3、<εε称为绝对误差限。相对误
4、差定义为绝对误差与精确值之比∆a∆=rx∆a∆=<η称为相对误差限rx�误差的来源:舍入误差将无限位字长的精确数处理成有限位字长近似数的处理方法称为舍入方法。带来舍人误差。截断误差用数值法求解数学模型时,往往用简单代替复杂,或者用有限过程代替无限过程所引起的误差。�有效数字对于a=a0a1…am.am+1…am+n(a0≠0)的近似数,若
5、Δ
6、≤0.5x10-n,则称a为具有m+n+1位有效数字的有效数,其中每一位数字都叫做a的有效数字。有效数和可靠数的最末位数字称为可疑数字有效数位的多少直接影响到近似值的绝对误差与相对误差的大小。推论1对于给
7、出的有效数,其绝对误差限不大于其最末数字的半个单位。mx=±0.aa⋯a×1012n∗1mn−∆=xx−≤×102推论2对于给出的一个有效数,其相对误差限可估计如下:mx=±0.aa⋯a×1012n5−n∆()x≤×10ra1例:计算y=lnx。若x≈20,则取x的几位有效数字可保证y的相对误差<0.1%?�数值计算的算法问题“良态”问题和“病态”问题在适定的情况下,若对于原始数据很小的变化δX,对应的参数误差δy也很小,则称该数学问题是良态问题;若δy很大,则称为病态问题。病态问题中解对于数据的变化率都很大,因此数据微小变化必将导致参数模型精
8、确解的很大变化。数学问题的性态完全取决于该数学问题本身的属性,在采用数值方法求解之前就存在,与数值方法无关。稳定算法和不稳定算法如果用数值方法计算时,误差在计算过程中不扩散的算法称为稳定算法。否则称为不稳定算法。�数值计算应注意的问题避免相近二数相减;避免小分母;避免大数吃小数;选用稳定的算法。绝对误差的运算:ε(x±x)=ε(x)+ε(x)1212ε(xx)≈xε(x)+xε(x)121221x1x1ε(x2)+x2ε(x1)ε()=2x2x2ε(f(x))≈fx′()(x)ε第三章线性方程组求解的数值方法教学内容:1.高斯消元法消元法的实现
9、过程;主元问题。2.矩阵分解矩阵LU分解的一般计算公式;利用LU分解的线性方程组求解方法;Cholesky分解;Matlab的Cholesky分解函数。3.向量范数与矩阵范数向量范数及其性质;矩阵函数及其性质;常用范数形式。4.线性方程组的迭代法求解迭代求解的思路;Jacobi迭代法;高斯_赛德尔迭代法;迭代法的收敛性。5.方程组的病态问题与误差分析线性方程组解的误差分析;条件数和方程组的病态性。消元法:问题:消去法是按照系数矩阵的主对角线上的元素(主元)进行消元。从而可能出现:(1)某个主元为零,导致消元过程无法进行。(2)当某个主元的绝对值
10、很小时,计算结果误差很大。定理:若A的所有顺序主子式均不为0,则高斯消元无需换行即可进行到底,得到唯一解。全主元消去法每一步选绝对值最大的元素为主元素。
11、a
12、=max
13、a
14、≠0;ijkkijk≤ij,≤n列主元消去法省去换列的步骤,每次仅选一列中最大的元。
15、a
16、max
17、=a
18、0≠ikk,ikkin≤≤矩阵三角分解法由Gauss消去法加上列主元(或全主元)有LU分解:A=LU⎛1⎞⎛u11u12⋯u1n⎞⎜⎟⎜⎟⎜l211⎟⎜u22⋯u2n⎟L=⎜l31l321⎟U=⎜⎟⎜⎟⎜⎟⎜⋯⋯⎟⎜⋯⋯⋯⎟⎜lll1⎟⎜u⎟⎝n1n2nn(−1)⎠⎝nn⎠
19、计算公式:u1j=a1j,j1,=⋯,nli1=aui111,i=2,⋯,n对k=2,3,⋯,n计算k-1ukj=akj−∑lukmmjj=k,⋯,nm1=k-1lik=(aik−∑luimmk)/ukki=+k1,⋯,nm1=算法:(1)解:Yy1=b1i-1yi=bi−∑lijyj,i=2,3,⋯,nj=1(2)解x:xn=ynunn⎛i+1⎞xi=⎜yi−∑uxijj⎟uii,⎝jn=⎠i=−n1,⋯,1Cholesky分解:定理:设矩阵A对称正定,则存在非奇异下三角阵L使得TA=LL。若限定L对角元为正,则分解唯一。Matlab中的Ch
20、olesky分解函数:chol()向量和矩阵的范数为了研究线性方程组近似解的误差估计和迭代法的收敛性,引进向量(矩阵)的范数的概念。向量范数定义:Rn