数值分析期末复习知识点
数值分析 知识点总结

数值分析知识点总结一、数值分析的基本概念1. 数值分析的对象数值分析的对象是现实生活中的数字数据和信息。
这些数据和信息可以来自各个领域,包括自然科学、社会科学、技术工程等。
例如,物理实验中测得的实验数据、经济管理中的统计信息、天气观测中的气象数据等,都是数值分析的对象。
2. 数值分析的目的数值分析的主要目的是通过对数值数据和信息的定量分析,发现其中的规律,提取有用的信息,做出科学的预测和决策。
例如,通过对某种药物的临床试验数据进行数值分析,可以得出这种药物的疗效和毒性情况,为临床医生的治疗决策提供依据。
3. 数值分析的方法数值分析采用数学和计算机科学的方法对数值数据和信息进行处理和分析。
它涉及的具体方法包括数值计算、插值与逼近、数值微分和积分、常微分方程数值解、数值线性代数等。
二、数值分析的基本内容1. 数值计算数值计算是数值分析的基本方法之一,它包括离散化、数值稳定性、误差分析等内容。
离散化是将连续问题转化为离散问题,这是数值计算的基本工作方式。
数值稳定性研究的是数值方法对误差的敏感程度,是评价数值方法好坏的重要指标。
误差分析则研究数值计算中产生的误差的成因和大小。
2. 插值与逼近插值与逼近是数值分析的重要内容之一,它研究如何通过已知的数值数据估计未知函数的值。
插值是通过已知的离散数据点构造一个连续函数,使得这个函数通过这些数据点;逼近则是通过已知的离散数据点构造一个近似函数,使得这个函数与原函数的差尽量小。
3. 数值微分和积分数值微分和积分是数值分析的又一重要内容,它研究如何通过已知的函数值计算函数的导数和定积分值。
数值微分是通过函数值计算函数的导数值;数值积分则是通过函数值计算函数的定积分值。
这两项工作在科学计算中有着广泛的应用。
4. 常微分方程数值解常微分方程数值解也是数值分析的重要内容之一,它研究如何通过数值方法计算常微分方程的近似解。
常微分方程是自然界和技术工程中经常出现的数学模型,因此其数值解的研究有着广泛的应用价值。
数值分析知识点总结

数值分析知识点总结数值分析是计算数值解的方法和理论,它研究的是如何利用计算机对数学问题进行数值计算和数值逼近。
数值分析包括了数值方法的设计、分析和实现,以及误差分析和计算复杂性分析等方面。
下面是数值分析的一些重要知识点的总结。
1.数值算法:数值算法是解决数学问题的计算方法,它由一系列具体的计算步骤组成。
常见的数值算法有插值、数值积分、数值微分、常微分方程数值解法等。
2.数值稳定性:数值稳定性是指数值算法在计算过程中对误差的敏感程度。
一个数值算法如果对输入数据的微小扰动具有较大的响应,就称为不稳定算法;反之,如果对输入数据的微小扰动具有较小的响应,就称为稳定算法。
3.四舍五入误差:在浮点数计算中,由于计算机表示的限制,涉及舍入运算的计算可能会引入误差。
四舍五入误差是指在进行舍入运算时,取最近的浮点数近似值所引入的误差。
4.条件数:条件数是用来衡量数值问题的不稳定性的一个指标。
它描述了输入数据的微小扰动在计算结果中的放大程度。
条件数的大小决定了数值算法的数值稳定性,通常越大表示问题越不稳定。
5.插值:插值是基于已知数据点,构造插值函数来近似未知数据点的方法。
常用的插值方法有线性插值、多项式插值和样条插值等。
6. 数值积分:数值积分是用数值方法进行积分计算的一种方法。
常见的数值积分方法有梯形法则、Simpson法则和Gauss-Legendre积分法等。
7.数值微分:数值微分是通过数值方法来计算函数的导数的一种方法。
常用的数值微分方法有中心差分法和前向差分法等。
8. 常微分方程数值解法:常微分方程数值解法用于求解常微分方程的近似解。
常用的常微分方程数值解法有Euler法、Runge-Kutta法和Adams法等。
9.误差分析:误差分析是对数值算法计算结果误差的研究。
可以通过理论分析或实验方法来估计误差,并找到减小误差的方法。
10.计算复杂性分析:计算复杂性分析是对数值算法运行时间和计算资源的需求进行评估的方法。
数值计算期末复习指南(整理版)

数值计算期末复习指南(整理版)
本文档旨在为数值计算的期末复提供指导。
以下是一些重要的复要点,帮助您进行有针对性的复。
1. 数值计算的基础知识
- 理解计算机中的数值表示方法和数值精度问题。
- 掌握计算机中数值的舍入误差和截断误差概念。
- 理解机器精度和有效数字的概念,并能计算有效数字。
2. 数值计算中的误差分析
- 熟悉误差分析的基本概念和方法。
- 掌握绝对误差和相对误差的计算方法。
- 理解截断误差和舍入误差在数值计算中的作用和影响。
3. 插值与逼近
- 理解插值和逼近的基本概念。
- 掌握常见的插值方法,如拉格朗日插值、牛顿插值等。
- 了解最小二乘逼近的原理和方法。
4. 数值积分
- 掌握数值积分的基本方法,如梯形公式、辛普森公式等。
- 理解数值积分的误差分析方法。
- 了解自适应积分和复化积分规则。
5. 数值微分方程的求解
- 熟悉常见数值求解常微分方程的方法,如欧拉法、龙格-库塔法等。
- 了解常微分方程初值问题和边值问题的数值求解方法。
- 掌握求解偏微分方程的基本方法。
请注意,本文档仅提供了数值计算复的基本要点,建议您结合教材和课堂笔记进行综合复。
祝您期末复顺利!。
数值分析考试复习总结

第一章1误差相对误差和绝对误差得概念 例题:当用数值计算方法求解一个实际的物理运动过程时 ,一般要经历哪几个阶段?在哪些阶 段将有哪些误差产生?答:实际问题-数学模型-数值方法-计算结果在这个过程中存在一下几种误差:建立数学模型过程中产生:模型误差 参数误差 选用数值方法产生:截断误差 计算过程产生:舍入误差传播误差6•设a =0.937关于精确数x 有3位有效数字,估计a 的相对误差.对于f(x^ .j_x ,估计f(a)对于f(x)的误差和相对误差I l /£、I I 匚 . a-x I .(2^10 . _ _3 | E( f)冃心 _x —G —a |= ------------ _,=] < ------------------ =10、1—x +H — a| 2 沃 0.25| E r(f)|E10,1 -a =4 10;.□2有效数字基本原则:1两个很接近的数字不做减法:2: 不用很小得数做分母(不用很大的数做分子)例题:4 •改变下列表达式使计算结果比较精确:a 的相对误差:由于1 _3x —a|E(x)|< x — <-10 .E r (X )=—2 XE r (x) < 12 7 2 1 2 10 =— 10 .18(Th1)解 f(a)对于f(x)的误差和相对误差第二章拉格朗日插值公式(即公式(1))插值基函数(因子)可简洁表示为n其中:n (x) - JI. 1n(X - X j ),nX i =「 (X i - X j )j /j工j料例1 n=1时,线性插值公式R(x) = y ° x( )+y 1 (x-X 0)X ----------------- ?(xo-xj ' (X 1 -X o )例2 n=2时,抛物插值公式 牛顿(Newton )插值公式由差商的引入,知(1) 过点x 0, x 1的一次插值多项式 为 其中(2) 过点x 0, x 1, x 2的二次插值多项式为其中重点是分段插值: 例题:(1)-1 01/2 1-3 -1/2 0 1 (2)-1 0 1/2 1-3/21/2解⑵:方法一.由Lagrange 插值公式(1) (2) 1 1 - x 1 2x 1 x(3)解⑴⑶1 - COS Xx对 x 0,| x 卜:::1.2x1(1x)(1 2x).⑵2 x(\ X 1 X 、X - 1 x)21 -cosx sin xsin x------------------ = ------------------------------ a s ---------------------x x(1 cosx) 1 cosx1 x可得:L3(x)=x2(x -1 2)方法二•令3 1由L3(-1)=-3,L3(1)=—,定A, B (称之为待定系数法)□2 215.设f(x) =x2,求f(x)在区间[0,1]上的分段线性插值函数f h(x),并估计误差,取等距节点,且h =1/10.解f(x) =x2,人=ih ,i =0,1, ,10,110设X i乞X乞X i 1 ,贝U:误差估计:第三章最佳一致逼近:(了解) 最佳平方逼近 主要分两种情形:1. 连续意义下在空间L 2[a,b ]中讨论2. 离散意义下在n 维欧氏空间R n 中讨论,只要求提供f 的样本值1. 最佳逼近多项式的法方程组设 L 2[a,b ]的 n 1 维子空间 P =span {1,x,x 2, x n },其中1,x,x 2…,x n 是L 2[a,b ]的线性无关多项式系. n对-f • L 2[a,b ],设其最佳逼近多项式''可表示为:''二a i x i i=0由(f - *, )=o, - P n即 n (x i ,x j )a * =(f ,x i ),0(1) n(*2)j=o其中称(*2)式为最佳逼近多项式的法方程组(或正规方程组) . 由{x [二的线性无关性,可证明G 正定,即 上述法方程组的解存在且唯一.11、求f(x)二cos 二X , X- [0,1]的一次和二次最佳平方逼近多项式. 解: 设 P ;(X) =a 0 a 1x , P ;(x)二 b 0 b,x b 2x 2 分别为f (x)的一次、二次最佳平方逼近多项式。
数值分析考试复习总结汇总

第一章1 误差相对误差和绝对误差得概念 例题:当用数值计算方法求解一个实际的物理运动过程时, 一般要经历哪几个阶段? 在哪些阶段将有哪些误差产生? 答: 实际问题-数学模型-数值方法-计算结果 在这个过程中存在一下几种误差:建立数学模型过程中产生:模型误差 参数误差选用数值方法产生:截断误差 计算过程产生:舍入误差 传播误差6.设937.0=a 关于精确数x 有3位有效数字,估计a 的相对误差. 对于x x f -=1)(,估计)(a f 对于)(x f 的误差和相对误差.解 a 的相对误差:由于 31021|)(|-⋅≤-≤a x x E . x ax x E r -=)(, 221018110921)(--⋅=⨯≤x E r . (1Th ))(a f 对于)(x f 的误差和相对误差.|11||)(|a x f E ---==()25.021011321⨯⋅≤-+---ax x a =310-33104110|)(|--⨯=-≤a f E r . □2有效数字基本原则:1 两个很接近的数字不做减法:2: 不用很小得数做分母(不用很大的数做分子)例题:4.改变下列表达式使计算结果比较精确:(1) ;1||,11211<<+--+x xxx 对(2) ;1,11>>--+x xx xx 对(3)1||,0,cos 1<<≠-x x xx对.解 (1) )21()1(22x x x ++. (2) )11(2x x x x x-++.(3) xxx x x x x cos 1sin )cos 1(sin cos 12+≈+=-. □ 第二章拉格朗日插值公式(即公式(1))∑==ni i i n x l y x p 0)()(插值基函数(因子)可简洁表示为)()()()()()(0i n i n nij j j i j i x x x x x x x x x l ωω'-=--=∏≠= 其中: ()∏∏≠==-='-=nij j j i i nnj jn x x x xx x 00)(,)()(ωω. 例1 n=1时,线性插值公式 )()()()()(010110101x x x x y x x x x y x P --⨯+--⨯=, 例2 n=2时,抛物插值公式))(())(())(())(())(())(()(1202102210120120102102x x x x x x x x y x x x x x x x x y x x x x x x x x y x P ----⨯+----⨯+----⨯= 牛顿(Newton )插值公式由差商的引入,知(1) 过点10,x x 的一次插值多项式为)()()(0101x x c x f x p -+=其中],[)()(1001011x x f x x x f x f c =--=⇒ )](,[)()(01001x x x x f x f x p -+=(2) 过点210,,x x x 的二次插值多项式为))(()()(10212x x x x c x p x p --+=其中],,[)()()()(21002010112122x x x f x x x x x f x f x x x f x f c =------=⇒ ))(](,,[)()(1021012x x x x x x x f x p x p --+=))(](,,[)](,[)(102100100x x x x x x x f x x x x f x f --+-+=重点是分段插值: 例题:1. 利用Lagrange 插值公式求下列各离散函数的插值多项式(结果要简化):(1)(2)解(2):方法一. 由 Lagrange 插值公式)()()()()(332211003x l f x l f x l f x l f x L ⋅+⋅+⋅+⋅= 可得: )21()(23-=x x x L 方法二. 令)()21()(3B Ax x x x L +-=由 23)1(3-=-L , 21)1(3=L , 定A ,B (称之为待定系数法) □15.设2)(x x f =,求)(x f 在区间]1,0[上的分段线性插值函数)(x f h ,并估计误差,取等距节点,且10/1=h .解 2)(x x f =, ih x i = , 10,,1,0 =i , 101=h设 1+≤≤i i x x x ,则: ii ii i i i i h x x x x x f x x x x x f x f --+--⋅=++++1111)()()(h ihx h i h h i x h i -++-+-⋅=22))1(()1()( 100)1(10)12(+-+=i i x i 误差估计: ))1(()(!2|)()(|max)1(h i x ih x f x f x f hi x ix h +--''≤-+≤≤. □第三章最佳一致逼近:(了解) 最佳平方逼近 主要分两种情形:1. 连续意义下在空间],[2b a L 中讨论2. 离散意义下在n 维欧氏空间n R 中讨论,只要求提供f 的样本值1. 最佳逼近多项式的法方程组设],[2b a L 的1+n 维子空间 n P =span },,,1{2n x x x , 其中 n x x x ,,,12 是],[2b a L 的线性无关多项式系.对],[2b a L f ∈∀,设其最佳逼近多项式*φ可表示为: ∑==ni i i x a 0**φ由 n P f ∈∀=-φφφ ,0),(*⇒ ∑===-ni j i i n j x x a f 0*)1(0,0),(即 ∑===nj ij j i n i x f a x x 0*)1(0),,(),((*2) 其中⎰⎰⎰⋅==⋅=+b ab abai iji jijidx x x f x f dx x dx x x x x)(),( ,),(称(*2)式为最佳逼近多项式的法方程组(或正规方程组). 由n i i x 0}{=的线性无关性,可证明G 正定,即 上述法方程组的解存在且唯一 .11、 求x x f πcos )(= ,]1,0[∈x 的一次和二次最佳平方逼近多项式. 解: 设 x a a x P 10*1)(+= , 2210*2)(x b x b b x P ++= 分别为)(x f 的一次、二次最佳平方逼近多项式。
数值分析考试复习总结

第一章1 误差相对误差和绝对误差得概念 例题:当用数值计算方法求解一个实际的物理运动过程时, 一般要经历哪几个阶段? 在哪些阶段将有哪些误差产生?答: 实际问题-数学模型-数值方法-计算结果 在这个过程中存在一下几种误差:建立数学模型过程中产生:模型误差 参数误差选用数值方法产生:截断误差 计算过程产生:舍入误差 传播误差6.设937.0=a 关于精确数x 有3位有效数字,估计a 的相对误差. 对于x x f -=1)(,估计)(a f 对于)(x f 的误差和相对误差.解 a 的相对误差:由于31021|)(|-⋅≤-≤a x x E . x ax x E r -=)(,221018110921)(--⋅=⨯≤x E r . (1Th ))(a f 对于)(x f 的误差和相对误差.|11||)(|a x f E ---==()25.021011321⨯⋅≤-+---ax x a =310-33104110|)(|--⨯=-≤a f E r . □2有效数字基本原则:1 两个很接近的数字不做减法:2: 不用很小得数做分母(不用很大的数做分子) 例题:4.改变下列表达式使计算结果比较精确:(1) ;1||,11211<<+--+x xxx 对(2);1,11>>--+x xx xx 对(3)1||,0,cos 1<<≠-x x xx对.解 (1) )21()1(22x x x ++. (2) )11(2x x x x x-++.(3) xxx x x x x cos 1sin )cos 1(sin cos 12+≈+=-. □第二章拉格朗日插值公式(即公式(1))插值基函数(因子)可简洁表示为其中: ()∏∏≠==-='-=nij j j i i nnj jn x x x xx x 0)(,)()(ωω. 例1 n=1时,线性插值公式 )()()()()(010110101x x x x y x x x x y x P --⨯+--⨯=, 例2 n=2时,抛物插值公式 牛顿(Newton )插值公式由差商的引入,知(1) 过点10,x x 的一次插值多项式为其中(2) 过点210,,x x x 的二次插值多项式为其中重点是分段插值:例题:1. 利用Lagrange 插值公式求下列各离散函数的插值多项式(结果要简化):(1) (2) 解(2):方法一. 由 Lagrange 插值公式 可得: )21()(23-=x x x L 方法二. 令由 23)1(3-=-L , 21)1(3=L , 定A ,B (称之为待定系数法) □15.设2)(x x f =,求)(x f 在区间]1,0[上的分段线性插值函数)(x f h ,并估计误差,取等距节点,且10/1=h .解 2)(x x f =, ih x i = , 10,,1,0Λ=i , 101=h设 1+≤≤i i x x x ,则: 误差估计: ))1(()(!2|)()(|max)1(h i x ih x f x f x f hi x ix h +--''≤-+≤≤. □第三章最佳一致逼近:(了解) 最佳平方逼近 主要分两种情形:1. 连续意义下在空间],[2b a L 中讨论2. 离散意义下在n 维欧氏空间n R 中讨论,只要求提供f 的样本值1. 最佳逼近多项式的法方程组设],[2b a L 的1+n 维子空间 n P =span },,,1{2n x x x Λ, 其中 n x x x ,,,12Λ是],[2b a L 的线性无关多项式系.对],[2b a L f ∈∀,设其最佳逼近多项式*φ可表示为: ∑==ni i i x a 0**φ由 n P f ∈∀=-φφφ ,0),(*即 ∑===nj ij j i n i x f a x x 0*)1(0),,(),((*2) 其中称(*2)式为最佳逼近多项式的法方程组(或正规方程组). 由n i i x 0}{=的线性无关性,可证明G 正定,即 上述法方程组的解存在且唯一 .11、 求x x f πcos )(= ,]1,0[∈x 的一次和二次最佳平方逼近多项式. 解: 设 x a a x P 10*1)(+= , 2210*2)(x b x b b x P ++= 分别为)(x f 的一次、二次最佳平方逼近多项式。
数值分析与试验期末复习资料
1. MA TLAB 计算中,在命令窗口运行语句f=polyval([2 3 1],2),返回结果f= . 2.设1)(2+=x x f ,则=]4,2[f .设函数13)(47+++=x x x x f ,则7阶差商]3,,3,3[710 f = . 8阶差商]3,,3,3[810 f = .若3219()5767f x x x =++,则()f x 的一阶差商[0,1]f = ,32阶差商0132[3,3,,3]f = .设(1)4,(2)6,(3)7f f f ===,则()f x 的二阶差商[1,2,3]f = . 3.在MATLAB 软件中,用于绘制平面数据散点图的函数为 .4.在求解方程组b AX =时,迭代格式f BX X)()1(+=+k k 对于任意初始向量)0(X 及任意f 都收敛的充要条件是 . 5.在函数22)(x x f =上任取四个互异点,通过这四个点的lagrange 插值多项式为 .在函数1)(2++=x x x f 上任取三个互异的点,通过这三个点的lagrange 插值多项式为 . 6.设4[1,1]()f x C -∈,已知节点1,21,0,13210===-=x x x x ,其相应的函数值为31(),0,0,22f x =-,则()f x 的三次Lagrange 插值多项式3()p x = .当2,1,1-=x 时, 4,0,3)(-=x f 则)(x f 的二次Lagrange 插值多项式为 .设(3)()[0,2]∈f x C ,已知节点0120,1,2===x x x ,其相应的函数值为()2,1,2f x =--,则()f x 的二次Lagrange 插值多项式的插值基函数1()l x = ,插值余项2()R x = . 7. MATLAB的值,可在命令窗口命令提示符后输入 .8.在MATLAB 操作中,把变量x ,y 定义为符号变量的语句为 . 9.设(0,1,,)j x j n =为互异节点,则lagrange 插值基函数满足0()nj j l x ==∑ .10.在MATLAB 软件中,进行MATLAB 操作的最主要的窗口称为 .11. 当1x >>改写为 . 12.用MATLAB 对一组数据进行多项式拟合的函数为 .13.用牛顿迭代法求2()1150=-=f x x 的正根时,迭代公式为 .14.MATLAB 中,用命令polyval 计算多项式13)(23++=x x x f 在100,,2,1,0 =x 时的值, 可在命令窗口中输入 . 用命令polyval 计算多项式124)(33+++=x x x x f 在2,1,1-=x 时的值,可在命令窗口中输入 . 用命令feval 要计算函数f1.m 在0x 处的值,在命令窗口中应输入 . 15.梯形求积公式的代数精度是 ,辛普森求积公式的代数精度是 . 16.求积公式2141()(0)(1)(2)333f x dx f f f ≈++⎰的代数精度是 . 17. 误差的来源大体可分为观测误差、 、 、 等四类.18.用二分法求3()251f x x x =--=0在[1,3]内的实根时,进行一步后根所在的区间为 ,进行二步后根所在的区间为 .19. 若x 的相对误差为3%,则nx 的相对误差为 .1.求积公式)1()1()(11f f dx x f +-≈⎰-在]1,1[-上具有( )次代数精确度.A. 1B. 2C. 3D. 4 2.下面对可进行LU 分解的矩阵A 的描述不正确的是( ).A. 分解所得的L 为单位下三角矩阵B. 分解所得的U 为上三角矩阵C. A 的顺序主子式可以等于零 D. 这种分解是唯一的3.通过点),(k k y x 、),(11++k k y x 的拉格朗日插值基函数)(x l k 、)(1x l k +满足( ).A. 0)(,0)(11==++k k k k x l x lB. 1)(,1)(11==++k k k k x l x lC. 0)(,1)(11==++k k k k x l x lD. 1)(,0)(11==++k k k k x l x l4.应用牛顿迭代法于方程03=-a x ,导出的求立方根3a 的迭代公式为( ).A. a x a x x x k k k k ---=+2313 B. ax ax x x k k k k --+=+2313 C. 2313k k k k x a x x x --=+ D. 2313k k k k x ax x x -+=+ 5.MA TLAB 命令窗口中,运行语句A=[1 2;3 4];A(1,2)^A(2,1),所得结果为( )A .6 B. 16 C. 8 D. 96.辛普森求积公式的代数精度为( ).A. 1B. 2C. 3D. 47.用雅可比迭代法求解方程组⎪⎩⎪⎨⎧=++=++=-+5223122321321321x x x x x x x x x ,则迭代矩阵)(B =.A. ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----022101220 B. ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---022110220 C. ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----522311122 D. ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-1221112218.通过点),(k k y x 、),(11++k k y x 的拉格朗日插值基函数)(x l k 、)(1x l k +满足( ).A. 0)(,0)(11==++k k k k x l x lB. 1)(,1)(11==++k k k k x l x lC. 0)(,1)(11==++k k k k x l x lD. 1)(,0)(11==++k k k k x l x l 9.关于用MA TLAB 库函数对方阵A 的操作下面叙述不正确的是( ).A. 运行diag (A ),可得一列向量B. 运行diag (diag (A )),可得一对角阵C. 运行triu (A )可得一上三角矩阵 D. 运行triu (A )可得一下三角矩阵 10.下面表达式是MATLAB 软件中合法变量名的是( )A. 3a_xB. ab_34C. a%3eD. bn+x 三、简答题1.请给出MATLAB 中M 函数文件的格式. 2.请给出数值积分中代数精度的概念. 3.请给出算法稳定的概念.4.请给出MATLAB 软件中分号、圆括号、方括号的功能.5.请给出数值分析中截断误差、舍入误差的概念. 6.请给出Matlab 软件中合法的变量名的命名规则.7.请给出数值计算中避免误差危害的至少四条原则. 8.请给出用迭代法解非线性方程时收敛阶的概念.四、解答题1. 线性方程组b AX =的系数矩阵A =11(0)1a a a a a a a ⎡⎤⎢⎥>⎢⎥⎢⎥⎣⎦,求能使解方程组的Jacobic 迭代法收敛的a 的范围.2. 用梯形公式,辛普森公式计算21sin ⎰xdx x,21ln(1+⎰dx ,10⎰x e dx .(小数点后保留四位).(已知:9975.0)5.1sin(,9093.02sin ,8415.01sin === ) (已知:ln2 0.6931,ln(1 0.8814,ln(10.7996===)(已知:122.7183, 1.6487==e e )3. 已知⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=2222222222222222A 求1A 、2A 、∞A .4. 叙述压缩映象不动点定理,并证明不动点的存在性.5. 正方形的边长大约为100cm ,应怎样测量才能使得其面积误差不超过1cm 2?6.用列主元消去法解线性方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--421641284247321x x x 7.用简单迭代法理论分析:对于任意的]4,0[0∈x ,由迭代格式 2,1,0,21=+=+k x x k k 得到的序列∞=0}{k k x 均收敛于同一个数*x .8.计算球体积要使相对误差限为1%,问度量半径R 所允许的相对误差限是多少? 9. 计算函数[]1(),0,12=-∈f x x x 关于[0,1]C 的1,f f ∞与2f .10. 已知数据表用复化梯形公式计算2.42.0()f x dx ⎰.11.(8分)已给数据表(取步长0.1h =) 用复化辛普森公式计算⎰4.10.1)(dx x f12. 设有解线性方程组b AX =的迭代格式g BX X k K +=+)()1(,其中A I B -=,如果A B ,的特征值全为正数,分析该迭代格式是否收敛.13. 设线性方程组b AX =的系数矩阵A =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111a a a a a a , 证明:当121<<-a 时,高斯-塞德尔迭代法收敛.14.取初始向量(0)(0,0,0)=T X ,用Jacobic 迭代法求解如下方程组,给出(1)(2)(3),,.X X X1231231232213225x x x x x x x x x +-=⎧⎪++=⎨⎪++=⎩ 15.(7分)计算61)f =≈1.4,若利用等式f =计算,试分析所得近似值的相对误差是多少? 16.(8分)设线性方程组b AX =的系数矩阵A =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-a a a 232131, 试求能使Jacobic 方法收敛的a 的范围.17.(8分)求方程0123=--x x 在5.10=x 附近的一个根,分析下列迭代公式是否收敛:在区间]6.1,3.1[上考察(1)2111kK x x +=+ (2)211)1(1-=+k K x x (已知4648.06.023=)18.(8分)若A 是n 级实矩阵,证明:F FA A An≤≤2119.(10分)设(0,1,,)j x j n =为互异节点,),...,1,0)((n j x l j =为拉格朗日插值基函数.证明: (1) 1)(0≡∑=nj j x l(2)⎩⎨⎧===∑=nk k x l kjnj j ,...,2,1001)0(0四、程序设计题(10分)1. 在某次阻尼振荡试验中测得如下表所列的9组数据点,已知阻尼振荡对应的函数模型为w x e )kx cos(a )x (f =,利用已知数据,求拟合函数的待定参数w k a ,,.(请写出用MATLAB 软件编程求解该题的代码)2. 用最小二乘法确定经验公式x be a y +=中的参数b a ,,使该曲线拟合下面的数据。
数值分析期末复习要点总结86页PPT
ቤተ መጻሕፍቲ ባይዱ 66、节制使快乐增加并使享受加强。 ——德 谟克利 特 67、今天应做的事没有做,明天再早也 是耽误 了。——裴斯 泰洛齐 68、决定一个人的一生,以及整个命运 的,只 是一瞬 之间。 ——歌 德 69、懒人无法享受休息之乐。——拉布 克 70、浪费时间是一桩大罪过。——卢梭
数值分析期末复习要点总结
6、法律的基础有两个,而且只有两个……公平和实用。——伯克 7、有两种和平的暴力,那就是法律和礼节。——歌德
8、法律就是秩序,有好的法律才有好的秩序。——亚里士多德 9、上帝把法律和公平凑合在一起,可是人类却把它拆开。——查·科尔顿 10、一切法律都是无用的,因为好人用不着它们,而坏人又不会因为它们而变得规矩起来。——德谟耶克斯
新数值分析复习要点
在c[a,b]中构造正交多项式 :
由1,
x,
x
2
,...,
x
n
,...)构造首1的正交多项式
n
(
x)
n0
0 (x) 1
k
1
(x)
(x
k
1 )
k
(x)
k
k
1 (
x)
k 0,1,2,...
k 1
(xk (x),k (x)) (k (x),k (x))
k 0,1,2,...
0
0, k
(k (x),k (x)) ( k 1 ( x), k 1 ( x))
2 j 1 2
1 1
f (x)Pj (x)dx,
j 0,1,..., n
n
最佳平方逼近元s(x) c j Pj (x) j0
在任意区间[a, b]上, 利用标准的正交多项式
例 : 利用Legendra正交多项式
x [a,b] t [1,1]
x batba,
2
2
t 1 (2x a b) ba
奇异值与奇异值矩阵
i
i ( AT A) 0,
i 1,..., r,
r
0
0 0
条件数 cond(A) p || A1 ||p|| A ||p , p F,1,2,
谱条件数 cond ( A)2 || A1 ||2|| A ||2
max ( AT A) 1 min ( AT A) n
A对称时,cond ( A)2 || A1 ||2|| A ||2
返回
四.矩阵的正交分解
Schmidt正交化法 P65-例2-34
Householder变换法 P67-例2-35
厦门大学 数值分析 期末考试复习整理
拉格朗日插值余项(余项定理):(1)()()()()()(1)!n nn n i i f R x f x L x x x n ξ+==-=-+∏n 次牛顿(Newton)插值公式为[][])())((,,)(,)()(110100100----++-+=n n n x x x x x x x x x f x x x x f x f x N由插值多项式的唯一性可知 N n (x ) ≡ L n (x ),故其余项也相同。
定理:Newton 插值多项式的余项为 R n (x )= f[x 0 ,x 1,… x n , x] ωn+1(x )其中ωn+1 (x )=(x - x 0)(x - x 1 )(x - x 2 )…(x - x n )注:一般当 x 靠近 x 0 时用前插,靠近 x n 时用后插,故两种公式亦称为表初公式和表末公式。
Newton 向前差分插值公式020000()()(1)(1)(1)1!2!!n n nN x N x th t t t t t t n f f f f n =+---+=+∆+∆++∆Newton 向后差分插值公式22()()(1)(1)2!(1)((1)(1)!n n n n n n n n n N x N x th t t f t f f t t t n f n =--=-∇+-∇+--++-∇法方程?,),(),(),( ),(),(),(),(),(),(),(),(),(1010101110101000⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡∴n n n n n n n n f f f a a a ϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕ曲线拟合(1) 直线拟合已知数据点:()m i y x i i ,,2,1,, =,设拟合直线为:x a a x y 10)(+=,则正规方程为: ⎪⎪⎩⎪⎪⎨⎧=+=+∑∑∑∑∑=====mi i i mi m i i i mi i m i i y x x a x ay x am a 1110211110(2) 多项式拟合 对于给定的一组数据(),,1,2,,i i x y i m=,寻求次数不超过n (n<<m ) 的多项式,2012nn y a a x a x a x=++++正规方程组0121011201n i n i i n i i n ii i n n n n i i n i i i a m a x a x y a x a x a x x y a x a x a x x y ++⎧+++=⎪+++=⎪⎨⎪⎪+++=⎩∑∑∑∑∑∑∑∑∑∑∑数值积分——插值型判断是否是插值型求积公式Newton-Cotes 公式)()()(0)(x C j ba nj n j f a b dxx f ⎰∑=-≈• 柯特斯系数解线性方程组的直接法 1)列主消元法 2)三角分解法迭代法矩阵的谱半径就是指矩阵的特征值中绝对值最大的那个,谱半径是矩阵的函数,但非矩阵范数.对任一矩阵范数有如下关系:ρ(A)≤║A║第七章非线性方程与方程组的数值解法1二分法2 迭代法不动点迭代法及其收敛性构造函数,时刻保持)(1k k x g x =+,不能单独考虑)(k x g 的导数<12.牛顿迭代法第9章 常微分方程初值问题数值解法 向前欧拉(Euler )公式),(1i i i i y x hf y y +=+9.2.2 梯形公式 []),(),(2111+++++=i i i i i i y x f y x f h y y。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第一章(有效数字位数)
1、经四舍五入取近似值,其绝对误差限不超过末尾数字的半个单位。
2、设X*为准确值,X为近似值,称e=X*-X为近似值X的绝对误差,简称误差(显然e可正可负,准确值X*未知,因此e的准确值无法求出)
3、|e|=|X-X*|≤ŋ,则称ŋ为近似值X的绝对误差限,简称误差限。
4、e r=e/X*称为相对误差,由于准确值X*总是未知的,所以也把e r*=e/X称为近似值X的相对误差
5、|e r*|=|e/X|≤ŋ*,则称ŋ*为近似值X的相对误差限
6、设X是X*的近似值,如果|X*-X|≤1/2×10-k,则称用X近似值表示X*时准确到小数点后第k位,并称从小数点后第k位起,直到最左边的非零数字之间的所有数字为有效数字,称有效数字的位数为有效数位。
7、设X是X*的近似值,X=±10m×0.a1a2…,其中a i(i=2,3…)是0到9之间的自然数,a1≠0,m为整数,如果|X*-X|≤1/2×10m-n,那么称近似值有n位有效数字。
8、四舍五入所得到的数均为有效数字,但并不是说非四舍五入所得到的数不能为有效数字。
第二章、非线性方程求根(不动点迭代、牛顿法、弦截法、快速弦截法、局部收敛、全局收敛、收敛阶)
1、不动点迭代法(迭代法)(单根区间求解方法):将非线性方程f(x)=0化为一个同解方程x=ø(x),若要求f(x*)=0,则x*=ø(x*),称x*为f(x)的零点,为ø(x)的一个不动点。
2、定理:设迭代函数ø(x)在【a,b】上连续,且满足(1)当x∈【a,b】时,a≤ø(x)≤b,(2)存在一正数L,满足0<L<1,且∀x∈【a,b】,有|ø/(x)|≤L<1。
则1、方程x=ø(x)在【a,b】内有唯一解x*。
2、对于任意初值x0∈【a,b】,迭代法
x k+1=ø(x k)均收敛x*
3、设ø(x)有不动点x*,如果存在x*的一个邻域 S:|X*-X|< ŋ,对任意初值x0∈S,迭代过程x k+1=ø(x k)均收敛,则称迭代过程在根x*邻近局部收敛。
4、设x*为ø(x)的不动点,ø/(x)在x*的邻域S连续,|ø/(x*)|<1,则迭代法局部收敛。
5、若存在实数p≥1和c≠0(c为常数)满足lim k→∞|e k+1|/|e k|p=c,则称迭代法p阶收敛,p=1为线性收敛,p=2为平方阶收敛。
6、定理:如果迭代法迭代函数ø(x)在根x*附近满足:(1)ø(x)存在p阶导数且连续,(2)ø/(x*)=ø//(x*)=…=ø(p-1)(x*)=0,而ø(p)(x*)≠0,则迭代法x k+1=ø(x k)收敛阶是p
7、牛顿迭代法(将非线性方程线性化,利用泰勒公式),迭代公式:x k+1=x k-f(x k)/ f/(x k)
8、牛顿法在单根附近至少平方阶收敛(f(x)=(x-x*)g(x)且g(x*)≠0)当x*是方程m重根时(m>1)则牛顿迭代法是线性收敛的。
(f(x)=(x-x*)m g(x)且g(x*))9、弦截法:牛顿迭代公式中f/(x k)用f(x k)-f(x0)/x k-x0代替,公式为x k+1=x k- (x k-x0)/ (f(x k)-f(x0))*f(x k)
10、快速弦截法:牛顿迭代公式中f/(x k)用f(x k)-f(x k-1)/x k-x k-1代替
第三章、线性方程组的迭代法{范数(向量范数、矩阵范数)、谱半径、严格对角占优、迭代法收敛的充分条件、迭代法收敛的充要条件、雅可比迭代法、高斯赛德尔迭代法}
1、常用向量范数
(1)||x||2=(|x1|2+|x2|2+…+|x n|2)1/2
(2)||x||1=|x1|+|x2|+…+|x n|
(3) ||x||∞=max|x i|(1≤i≤n)
(4)||x||p=(|x1|p+|x2|p+…+|x n|p)1/p
2、矩阵范数
(1)设A=(a ij)∈R n*n,则(1)||A||1=max1≤j≤n∑i=1~n|a ij|,A的每列绝对值之和的最大值,称为A的列范数
(2)||A||∞,A的每行绝对值之和的最大值,称为A的行范数
(3)||A||2=√λmax(A T A), λmax(A T A)为A T A的特征值的绝对值的最大值
3、设A∈R n*n的特征值为λ1,λ2,λ3,…λn,称ρ(A)=max{|λ1|,|λ2|,…|λn|}为矩阵A的谱半径
4、严格对角占优:主对角线元素绝对值>同行其他元素绝对值之和
5、迭代法收敛的充分条件:迭代矩阵A的某种范数<1
6、迭代法收敛的充要条件:(1)ρ(A)<1、(2)矩阵符合严格对角占优
7、雅可比迭代法(迭代矩阵,迭代公式)
8、高斯赛德尔迭代法(迭代矩阵,迭代公式)
第五章、拉格朗日插值与插值余项定理、牛顿插值与差商、分段线性插值(插值区间的选取)
拉格朗日插值
1、线性插值:p1(x)=l0(x)y0+l1(x)y1
记l0=(x-x1)/x0-x1 l1=(x-x0)/x1-x0
2、抛物插值公式
P2(x)= l0(x)y0+l1(x)y1+l2(x)y2
记l0(x)=(x-x1)(x-x2)/ (x0-x1)(x0-x2) l1(x)=(x-x0)(x-x2)/ (x1-x0)(x1-x2) l2=(x-x0)(x-x1)/(x2-x0)(x2-x1)
3、p n(x)=Σi=0·n y i l i(x)
插值余项定理 R n (x )=f (x )-P n (x )={f (n+1)(ξ)/(n+1)!}*∏j=0·n (x-x j ) ξ∈(a ,b ) 牛顿插值
1、 线性插值
P 1(x )=p 0(x )+c 1(x-x 0)
记c 1={f (x 1)-f (x 0)}/(x 1-x 0)
2、 牛顿插值公式
P n (x )=f (x 0)+f[x 0,x 1](x-x 0)+…+f[x 0,x 1,…x n ](x-x 0)…(x -x n-1)
差商
1、一阶差商:f[x 0,x 1]=f(x 1)-f(x 0)/x 1-x 0
2、二阶差商:f[x 0,x 1,x 2]=f[x 1,x 2]-f[x 0,x 1]/x 2-x 0
3、n 阶差商:f[x 0,x 1,…,x n ]=f[x 1,x 2,…x n ]-f[x 0,x 1,…x n-1]/x n -x 0
差商性质
f[x 0,x 1,…x n ]=f (n)(ξ)/n!
分段线性插值(插值区间的选取)
第六章(代数精度、求积公式{左矩形、右矩形、中矩形、梯形、辛普生、复化梯形、复化辛普生、变步长梯形}、求积余项)
1、代数精度:m 次准确成立,m+1次不准确,则称此求积公式具有m 阶代数精度
2、求积公式
左矩形:∫f (x )dx b
a
=f (a )(b-a ) 右矩形:∫f (x )dx b a
=f (b )(b-a ) 中矩形:∫f (x )dx b a
=f (a+b/2)(b-a ) 梯形求积公式:∫f (x )dx b a
=(b-a )*[f(a)+f(b)]/2 (1次代数精度) 辛普生求积公式:∫f (x )dx b a
=(b-a )[f(a)+4f(a+b/2)+f(b)]/6 复化梯形公式:∫f (x )dx b a
=h/2(f (a )+2Σk=1·(n-1)f (x k )+f (b )) h=b-a/n 复化辛普生公式: ∫f (x )dx xn x0
=h/6(f (a )+4Σk=0·n-1f (x k+1/2)+2Σk=1·n-1f (x k )+f (b ))
变步长梯形求积公式h=b-a/n。