数值分析期末复习(整理版)
期末数值分析重点总结
期末数值分析重点总结第一部分:数值逼近(Approximation)数值逼近是数值分析的基础,主要研究如何利用有限的计算资源得到逼近数学问题的有效算法。
数值逼近的主要内容包括多项式逼近、插值和最小二乘等。
1. 多项式逼近多项式逼近是指用一个多项式函数来逼近给定函数的值。
通过选择合适的多项式次数和插值点,可以使得多项式逼近误差最小化。
其中最常用的方法是最小二乘法,它可以通过最小化残差来得到最佳的多项式逼近。
多项式逼近在信号处理、图像处理和计算机图形学等领域中有广泛的应用。
2. 插值插值是指通过已知数据点的函数值来估计在其他点的函数值。
常用的插值方法有拉格朗日插值和牛顿插值。
拉格朗日插值通过构造一个满足插值条件的多项式来逼近给定函数。
牛顿插值则利用差商的概念来构造插值多项式。
插值方法在数值微分和数值积分中有广泛的应用。
3. 最小二乘最小二乘是一种在一组离散数据点上拟合曲线的方法。
通过最小化数据点与拟合曲线之间的欧几里得距离,可以得到最佳拟合曲线。
最小二乘法可以用于曲线拟合、参数估计和数据关联等问题。
第二部分:数值解方程(Numerical Solution of Equations)数值解方程是数值分析的重要内容之一,研究如何通过数值计算来求解非线性方程组和线性方程组。
数值解方程的主要方法有迭代法、常微分方程数值解和偏微分方程数值解等。
1. 迭代法迭代法是求解非线性方程组的常用方法之一。
通过不断迭代逼近方程的根,可以得到方程组的数值解。
常用的迭代法有牛顿迭代法和弦截法。
迭代法在计算机辅助设计、优化和数据分析等领域中有广泛的应用。
2. 常微分方程数值解常微分方程数值解研究如何通过数值计算来求解常微分方程。
常微分方程数值解的主要方法有Euler方法、Runge-Kutta方法和线性多步法等。
常微分方程数值解在物理学、工程学和生物学等领域中有广泛的应用。
3. 偏微分方程数值解偏微分方程数值解研究如何通过数值方法来求解偏微分方程。
(完整版),数值分析笔记期末复习汇总,推荐文档
x
*n )
e(x *1)
f
(x *1,
x *2 ,, xn
x *n
)
e(x *n )
n i 1
f
(x *1, x *2 ,, x *n ) xi
e(x *i )
9、加减乘除运算的误差估计
加法
绝
对 误
e(x1 x2 ) e(x1) e(x2 )
差
绝
对
误 (x1 x2 ) (x1) (x2 )
x1
b
sign(b) 2a
b2 4ac 109
x1
x2
c a
x2
c a x1
109 109
1
求和时从小到大相加,可使和的误差减小。若干数相加,采用绝对值较小者先加的算法,
结果的相对误差限较小
y 54321100 0.4100 0.3100 0.4100 54322
(三) 注意简化计算步骤,减少运算次数,避免误差积累(秦九韶)
则称 r (x*) 为近似值 x*的相
对误差限。 (2)性质:
当|| er (x*) | 较小时,可用下
是有量纲的。 (2)绝对误差限是正的,有无穷
常取
er
( x*)
e( x*) x*
式计算
绝对误差是误差的绝对值? 多个【则比 * 大的任意正数均
(错)
是绝对误差
限】
r
( x*)
(x*) | x |
取
x2* =3.14
作为 π 的近似值,则 | e2
| 0.00159
1 102 :三个有效数字 2
取
x3* =3.1416 作为 π 的近似值,则 | e3
| 0.00000734
数值分析期末复习
有效数字
x=p =3.14159265…
按照四舍五入的原则可以得到x的前几位.
取3位x3*=3.14, e3*≤0.002, 取5位x5*=3.1416, e5*≤0.000008,
一般地,误差小于末位数字的半个单位
|p-3.14|≤0.5×10-2, |p -3.1416|≤0.5×10-4
28
引言
如果要求近似函数取给定的离散数据,则称之 为插值函数.实用上,我们常取结构相对比较简 单的代数多项式作为插值函数,这就是所谓的 代数插值. 本章先讨论代数插值,然后在此基础上进一步 研究所谓的样条插值.
29
引例
三角函数表的构造 如何用初等方法给出sin x在一系列点处的函 数值? 已知sin x在x=0,p/6,p/4,p/3,p/2等处的函数值. 将这些点处的函数值在图形上标出.
38
线性插值
问题 求作一次式 L1(x),满足条件 L1(xk)=yk, L1(xk+1)=yk+1, 从几何图形上看,y= L1(x)表示过两点(xk,yk), (xk+1,yk+1)的直线,因此可表为如下对称形式: L1(x)=yklk(x)+yk+1lk+1(x) 其中
En In In (1 nIn1) (1 nIn1)
n( In1 In1 )
nEn1
因此有En=(-1)nn!E0. 这说明第n步的误差是初始误差的n!倍. 下面的表给出了实际的误差.
12
例1的误差分析
n
In
En
n In
En
0 0.6321 2.056 10–5 5 0.1480 –2.467 10–3
数值分析笔记期末复习
第一章引论1、数值分析研究对象:数值分析是计算数学的一个主要部分,计算数学是数学科学的一个分支,它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现。
2、数值分析特点:①面向计算机,要根据计算机特点设计切实可行的有效算法②有可靠的理论分析,能任意逼近并达到精度要求,对近似计算要保证收敛性和数值稳定性③要有好的计算复杂性,时间复杂性好是指节省时间,空间复杂性好是指节省存贮量,这也是建立算法要研究的问题。
④要有数值试验,即任何一个算法除了从理论上要满足上述三点外,还要通过数值试验证明是行之有效的。
3、数值分析实质:是以数学问题为研究对象,不像纯数学那样只研究数学本身的理论,而是把理论与计算紧密结合,着重研究数学问题的数值方法及理论。
4、用计算机解决科学计算问题通常经历以下过程实际问题--数学模型(应用数学)--数值计算方法--程序设计--上机计算结果(计算数学)5、误差来源及分类1.模型误差——从实际问题中抽象出数学模型2.观测误差——通过测量得到模型中参数的值(通常根据测量工具的精度,可以知道这类误差的上限值。
)要用数值计算方法求它的近似解,由此产生的误差称为(截断误差)或(方法误差)原始数据的输入及浮点数运算过程中都有可能产生误差,这样产生的误差称为舍入误差6、五个关于误差的概念5.有效数字(1)定义:若近似值x*的绝对误差限是某一位的半个单位,该位到x*的第一位非零数字一共有n 位,则称近似值x*有n 位有效数字,或说x*精确到该位。
注意:近似值后面的零不能随便省去!2≤⨯1102≤⨯10.00000734102≤⨯(3)性质:(1)有效数字越多,则绝对误差越小 (2)有效数字越多,则相对误差越小有效数字的位数可刻画近似数的精确度! 6、一元函数的误差估计问题:设y =f (x ),x 的近似值为x *,则y 的近似值 y *的误差如何计算?(*)(*)(*)(*)e y dy f x dx f x e x ''≈=≈ (*)(*)(*)e y f x e x '≈ *(*)(*)(*)(*)r r x e y f x e x f x '≈故相应的误差限计算如下(*)(*)(*)y f x x εε'≈ *(*)(*)(*)(*)r r x y f x x f x εε'≈ 7、二元函数的误差估计问题:设y=f(x1, x2), x1, x2的近似值为x1*, x2* ,则y 的误差如何计算?**121212(*)*(,)(,)(*,*)e y y y f x x f x x df x x =-=-≈12121212(*,*)(*,*)(*)(*)f x x f x x e x e x x x ∂∂=+∂∂(*)(*)*(*)(*)(*)(*)(*)r r dy f x e x x e y f x e x y f x f x ''≈≈=1212121212121212(,)(,)(*,*)(*,*)(*)()()(*)(*)f x x f x x f x x f x x e y e x e x e x e x x x x x ∂∂∂∂=+≤⋅+⋅∂∂∂∂故绝对误差限为12121212(*,*)(*,*)(*)(*)(*)f x x f x x y x x x x εεε∂∂=+∂∂8、多元函数的误差估计121211121(*,*,,*)(*,*,,*)(*)*(*)(*)(*,*,,*)(*)n n n nnn i i i f x x x f x x x e y y y e x e x x x f x x x e x x =∂∂=-=++∂∂∂=∂∑9、加减乘除运算的误差估计(1)定义:初始数据的误差或计算中的舍入误差在计算过程中的传播,因算法不同而异。
(完整版)数值分析整理版试题及答案,推荐文档
9
1
xdx T4
h[ 2
f
1
3
2 k 1
f
xk
f
9]
2[ 1 2 3 5 7 9] 2
17.2277
(2)用 n 4 的复合辛普森公式
由于 h 2 , f x
x
,
xk
1
2k k
1, 2,3,
x
k
1
2
2k k
0,1, 2,3,所以,有
2
3
9
1
xdx S4
h[ 6
f
1
若 span1, x,则0 (x) 1 ,1(x) x ,这样,有
2
1
0 ,0 1dx 1
0
1,1
1 0
x2dx
1 3
0
,1
1,0
1
0
xdx
1 2
1
f ,0 exdx 1.7183
0
1
f ,1 xexdx 1
0
所以,法方程为
1
1
1
2 1
a0
a1
1.7183 1
1 0
1
23
2 1
a0
a1
6 1
12
3
再回代解该方程,得到
a1
4
,
a0
11 6
故,所求最佳平方逼近多项式为
S1*
(
x)
11 6
4x
例 3、 设 f (x) ex , x [0,1] ,试求 f (x) 在[0, 1]上关于 (x) 1 , span1, x的最
佳平方逼近多项式。 解:
1
4
x1
1 5
数值分析期末知识点总结
数值分析期末知识点总结一、引言数值分析是一门研究如何使用计算机提高数学模型数值计算精度和效率的学科。
它是计算数学的一个重要分支,涉及到数值计算、数值逼近和误差分析等一系列内容。
在数值分析课程中,我们将学习到数值解微分方程、线性代数问题的求解、插值与拟合、积分等一系列内容。
本文将对数值分析期末知识点进行总结,以便帮助大家复习。
二、常见数值计算方法1. 插值与拟合插值与拟合是数值分析中重要的内容,它们用于在给定数据点集上构造一个函数,以便在其他点上进行求值。
插值是通过一些已知数据点来求得一个函数,使得这个函数能够通过这些点,而拟合则是通过已知数据点来求得一个函数,使得这个函数在这些点附近能够比较好地拟合数据。
常见的插值方法包括线性插值、拉格朗日插值、牛顿插值等;而拟合方法包括最小二乘法拟合、多项式拟合等。
2. 数值解微分方程数值解微分方程是数值分析的一个重要内容,它讨论如何使用计算机对微分方程进行数值求解。
微分方程是自然界中描述变化的数学方程,它们在物理学、化学、生物学等领域都有着重要的应用。
数值解微分方程的方法包括欧拉法、中点法、四阶龙格-库塔法等。
3. 数值线性代数数值线性代数是数值分析领域的另一个重要内容,它讨论如何使用数值方法解决线性代数问题。
原始的线性代数问题可能非常大或者非常复杂,因此我们常常需要使用计算机进行数值计算。
数值线性代数的方法包括高斯消元法、LU分解、Jacobi迭代法、Gauss-Seidel 迭代法等。
4. 数值积分数值积分是数值分析的一个重要内容,它讨论如何使用数值方法对积分进行数值求解。
在实际问题中,有很多积分问题是无法解析求解的,因此我们需要使用数值方法进行近似求解。
数值积分的方法包括复合辛普森法、复合梯形法、龙贝格积分法等。
三、数值分析的误差分析在数值计算过程中,我们会遇到误差的问题。
这些误差可能来自于测量、舍入、截断等各种原因。
因此,误差分析是数值分析中一个非常重要的内容。
《数值分析》期末复习题(1)
《数值分析》期末复习题一、单项选择题1. 数值x *的近似值x =0.32502×10-1,若x 有5位有效数字,则≤-*x x ( ).(A)21×10-3 (B) 21×10-4 (C) 21×10-5 (D) 21×10-6 2. 设矩阵A =10212104135⎡⎤⎢⎥⎢⎥⎣⎦,那么以A 为系数矩阵的线性方程组A X =b 的雅可比迭代矩阵为( )(A)00.20.10.200.40.20.60--⎡⎤⎢⎥--⎢⎥--⎣⎦(B)10.20.10.210.40.20.61⎡⎤⎢⎥⎢⎥⎣⎦(C) 00.20.10.200.40.20.60⎡⎤⎢⎥⎢⎥⎣⎦ (D)021204130⎡⎤⎢⎥⎢⎥⎣⎦3. 已知(1)1,(2)4,(3)9f f f ===,用拉格朗日2次插值,则(2.5)f =( )(A) 6.15 (B) 6.25 (C) 6.20 (D) 6.10 4. 抛物形求积公式的代数精度是( )A. 1,B. 2 ,C. 3,D. 45. 改进欧拉格式的 局部截断误差是( ). (),A O h 2. (),B O h 3. (),C O h 4. ().D O h二、填空题1、以722作为π的近似值,它有( )位有效数字; 2、经过)1,2( ),2,1( ),1,0(C B A 三个节点的插值多项式为( ); 3、用高斯-赛德尔迭代法解方程组⎩⎨⎧-=+-=+,10,232121x bx bx x 其中b 为实数,则方法收敛的充分条件是b 满足条件( );4、取步长为1.0=h ,用欧拉法计算初值问题22',(0)0,y x y y ⎧=+⎨=⎩的解函数)(x y ,它在3.0=x 的近似值为( );5、已知方程0sin 1=--x x 在)1,0(有一个根,使用二分法求误差不大于41021-⨯的近似解至少需要经过( )次迭代。
(整理)《数值分析》期末复习纲要.
《数值分析》期末复习纲要 第一章 数值计算中的误差分析主要内容(一)误差分析 1、误差的基本概念:(1)绝对误差:设x 是精确值, *x 是其近似值,则称()E x x x*=-是近似值*x 的绝对误差,简称误差。
特点:可正可负,带量纲。
(2)相对误差:称()r x x E x x *-=是近似值*x 的相对误差,若精确值x 未知,则定义()r x x E x x **-=。
注: 由四舍五入得到的近似值,误差不超过最末位的半个单位(准确到最末位)。
2、有效数字的概念:P6;3、算法的数值稳定性:数值稳定的算法:初始数据所带有的误差在计算的过程中能得到有效控制,不至于因误差的过度增长影响计算结果的精度。
数值不稳定的算法:初始数据所带有的误差在计算的过程中得不到有效控制,以至于因误差的过度增长而使计算结果的精度大大降低。
P11:例子(二)算法设计的基本准则P11-15 应用实例:课堂练习,作业基本要求1、掌握误差、有效数字等基本概念2、熟记算法设计准则,并能依据算法设计准则构造或选择计算公式。
(参见课堂练习、作业)第二章 线性代数方程组的数值解法直接法:不计初始数据的误差和计算过程中的舍入误差,经过有限步四则运算求得方程组的精确解。
迭代法:先给出方程组解的某一初始值,然后按照一定的迭代法则(公式)进行迭代,经过有限次迭代,求得满足精度要求的方程组的近似解。
主要内容(一)直接法的基本模式:高斯顺序消去法基本思想:按照各方程的自然排列顺序(不交换方程),通过按列消去各未知元,将方程组化为同解的三角形方程组来求解求解过程:⎩⎨⎧回代过程消元过程应用实例:课堂例题;练习 (二)高斯列主元消去法基本思想:按列消元,但每次按列消元之前,先选取参与消元的 方程首列系数,选取绝对值最大者,通过交换方程,使之成为主元,再进行消元。
(每一步消元之前先按列选取主元) 应用实例:课堂例题,作业(三)迭代法基本原理:(1)将原方程组b Ax =改写成如下等价形式:f Bx x += (2)构造相应的迭代公式:f Bx x m m +=-)1()((3)任取一初始向量)0(x代入上述迭代公式,经迭代得到向量序列{}Tm n m m m x x x x ),,,()()(2)(1)( =,如果该向量序列{})(m x 收敛于某一向量Tn x x x x ),,,(21****= ,即),,2,1(lim )(n i x x i m i m ==*∞→Tn x x x x ),,,(21****= 即为原方程组的解。
数值分析期末复习题
数值分析期末复习题⼀、填空题1.设真值x=983350,则其近似值y=98000的有效数字的位数,绝对误差为,相对误差为。
2.x=0.1062,y=0.947,计算x+y 其有效数字的位数为。
3.对f(x)=x 3+x+1,差商f[0,1,2,3]= ;f[0,1,2,3,4]= 。
4.设f(x)可微,求⽅程x=f(x)根的⽜顿迭代法格式是。
5.设⽅程x=?(x)有根x *,且设?(x)在含x *的区间(a,b)内可导,设x 0∈(a,b)则迭代格式x k+1=?(x k )收敛的充要条件为。
6.求解线性⽅程组Ax=b 的迭代格式x (k+1)=Jx (k)+f 收敛的充要条件为。
7.=011001001001....A ,||A||∝= ,cond(A)∝= 。
8.n 次Legendre 多项式的最⾼次项系数为。
9.中矩形公式:)()2()(a b b a f dx x f b a -+=?的代数精度为。
10.求积公式:)1(21)0()(10f f dx x f '+≈?的代数精度为。
11.在区间[1,2]上满⾜插值条件??==3)2(1)1(P P 的⼀次多项式P(x)= 。
12.设∑==n k k k n x f A f I 0)()(是函数f(x)在区间[a,b]上的插值型型求积公式,则 ∑=n k k A0= 。
13.梯形公式和改进的Euler 公式都是阶精度的。
⼆、计算题1.利⽤矩阵的⾼斯消元法,解⽅程组=++=++=++2053182521432321321321x x x x x xx x x2.设有函数值表试求各阶差商,并写出Newton 插值多项式。
3.求解超定⽅程组= ?43231211121x x的最⼩⼆乘解。
4.给定下列函数值表:求3次⾃然样条插值函数5.给定x x f =)(在x=100, 121, 144 三点处的值,试以这三点建⽴f(x)的⼆次(抛物)插值公式,利⽤插值公式求115的近似值并估计误差。
数值分析总复习
A
4
5
4,
X
x2
,
8 4 22
x3
解: l11 a11 16 4,
l21 a21 l11 4 4 1,
l31 a31 l11 2,
4
b
3
.
10
l11 a11 , l21 a21 l11 , l22 a22 l221 ,
l31 a31 l11 , l32 a32 l21l31 l22 , l33 a33 l321 l322 . 19 第20页/共36页
l11 a11 , l21 a21 l11 , l22 a22 l221 ,
l31 a31 l11 , l32 a32 l21l31 l22 , l33 a33 l321 l322 . 18 第19页/共36页
一. 用平方根法求线性方程组AX=b, 其中
16 4 8
x1
26
第27页/共36页
六. 确定求解初值问题
y' f ( x, y), a x b,
y(a)
y0 .
的二步隐式Adams方法
yn1
yn
h 12
(5
fn1
fn
fn1 )
中的参数, 使该方法成为三阶方法, 并写出其局部截断误差主项.
可用数值积分方法或Taylor展开方法
8,
Rn1
1 24
h4
解 (1) 由已知, 当 f (x)分别为1, x, x2时, 求积公式等号成立. 即
11x3dx 1
0 1dx 14
11 2
((1x13
1)x23
)
2
故该公式具有3次代数精确度.
1 xdx 1
0
数值分析总复习
样条插值;整体连续光滑,且不需知导数值。
插值问题提法:已知
x y f(x)
x0 y
x1 y
xn y
0
1
n
求一个三次分段函数 S(x) 使
1,
S(
xi
)
y i
x x 2, 在 [ , ] 上是三次多项式
i
i 1
C 3, S(x) 2 ( a,b )
i 0, 1, , n
计算三次样条算法
由边界条件 i , i , , i 0 ,1,, n
插值基函数方法
插值问题解的一般形式 :
n (x) a0 a1 x an xn
(1 )
实质上是在求多项式的 自然基底 Bn Span{1, x , ,xn}
张成的线性空间中的一 个点 —一个多项式 (1) ,由(2 18)
式知,解存在唯一 ,只要解方程组求出线 性组合系数 {ai}
就可以了 , 但计算量太大 .
定理2.5(余项) .
(2 - 35)
设H (x)是过 x0 , x1 的 Hermite 插值多项式 , C f f(x) 3 , ( 4 )(x)在 (a,b) 内存在, (a,b)是
(a,b)
含点 x0 , x1 的任一区间, 则对任意给定的
x (a,b) 总存在一点ξ (x)使
R(x)
f(x) H(x)
f
( 4 )(ξ
4!
)
(x
x0
)2(x
x1
)2
分段三次 Hermite 插值多项式及余项
∑ y h m H n
H (x) [ (x)
( x)]
i0
ii
ii
定理2.7(余项) :
数值分析期末复习总结(优选.)
线性插值多项式(一次插值多项式)
n=2
L2 ( x) =
y0
(x ( x0
− −
x1 )( x − x2 ) x1 )( x0 − x2 )
+
y1
(x ( x1
− −
x0 )( x − x2 ) x0 )( x1 − x2 )
+
y2
(x ( x2
− −
x0 )( x − x1 ) x0 )( x2 − x1 )
f ( x=) f ( x0 ) + ( x − x0 ) f [x, x0]
1
f [ x, x0 ] = f [ x0 , x1] + ( x − x1 ) f [ x, x0 , x1]
2
……
f [ x, x0 , ... , xn−1] = f [ x0 , ... , xn ] + ( x − xn ) f [ x, x0 , ... , xn ] n−1
19
Newton 插值
为什么 Newton 插值
Lagrange 插值简单易用,但若要增加一个节点时,全部基函
数 lk(x) 都需重新计算,不太方便。
解决办法
设计一个可以逐次生成插值多项式的算法,即 n 次插值多项式 可以通过 n-1 次插值多项式生成 —— Newton 插值法
20
新的基函数
设插值节点为 x0 , … , xn ,考虑插值基函数组 ϕ0(x) = 1 ϕ1( x)= x − x0 ϕ2( x) = ( x − x0 )( x − x1 )
18
插值余项
几点说明
余项公式只有当 f(x) 的高阶导数存在时才能使用
ξx 与 x 有关,通常无法确定, 实际使用中通常是估计其上界
数值分析期末复习题
一、填空题1.设真值x=983350,则其近似值y=98000的有效数字的位数 ,绝对误差为 , 相对误差为 。
2.x=0.1062,y=0.947,计算x+y 其有效数字的位数为 。
3.对f(x)=x 3+x+1,差商f[0,1,2,3]= ;f[0,1,2,3,4]= 。
4.设f(x)可微,求方程x=f(x)根的牛顿迭代法格式是 。
5.设方程x=ϕ(x)有根x *,且设ϕ(x)在含x *的区间(a,b)内可导,设x 0∈(a,b)则迭代格式x k+1=ϕ(x k )收敛的充要条件为 。
6.求解线性方程组Ax=b 的迭代格式x (k+1)=Jx (k)+f 收敛的充要条件为 。
7.⎪⎪⎭⎫ ⎝⎛=011001001001....A ,||A||∝= ,cond(A)∝= 。
8.n 次Legendre 多项式的最高次项系数为 。
9.中矩形公式:)()2()(a b b a f dx x f b a -+=⎰的代数精度为 。
10.求积公式:)1(21)0()(10f f dx x f '+≈⎰的代数精度为 。
11.在区间[1,2]上满足插值条件⎩⎨⎧==3)2(1)1(P P 的一次多项式P(x)= 。
12.设∑==n k k k n x f A f I 0)()(是函数f(x)在区间[a,b]上的插值型型求积公式,则 ∑=n k k A0= 。
13.梯形公式和改进的Euler 公式都是 阶精度的。
二、计算题1.利用矩阵的高斯消元法,解方程组⎪⎩⎪⎨⎧=++=++=++2053182521432321321321x x x x x xx x x2.设有函数值表试求各阶差商,并写出Newton 插值多项式。
3.求解超定方程组⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛43231211121x x的最小二乘解。
4.给定下列函数值表:求3次自然样条插值函数5.给定x x f =)(在x=100, 121, 144 三点处的值,试以这三点建立f(x)的二次(抛物)插值公式,利用插值公式求115的近似值并估计误差。
数值分析期末复习
xn 0.1 0.2 0.3 0.4 0.5
yn 1.0959 1.1ቤተ መጻሕፍቲ ባይዱ41 1.2662 1.3434 1.4164
在计算时,迭代终止的时间可以用上式判别
例. 判别下列方程组用J法和G-S法求解是否收敛
1 1 2 2 1 2 - 2 x1 1 1 x2 = 1 1 x3 1
ttttaluyuxx???????????????????????????求解242649615269186151840??????????????用lu直接三角分解法求解方程组axb其中a9232247????????????b11242621123121363321119ay???????????????????????????????????????????????????????????????????234212312122332147953105231ttylybyyyuxyx?求解即得求解得
所以Gauss-Seidel迭代法发散
说明G-S法发散时而J法却收敛 因此,不能说G-S法比J法更好
例: 已知x=1,4,9的平方根为1,2,3,利用牛顿基本差商
公式求 7 的近似值。 解: x i
xi
1
2
f [ xi , xi +1 ]
f [ xi , xi +1 , xi + 2 ]
1 4
9
3
第三章 函数逼近的基本概念
曲线拟合的最小二乘法 已知一组实验数据,求它的拟合曲线。 线性化 建立法方程组 求解未知变量 给出拟合曲线函数
第四章 数值积分与数值微分
数值积分的基本思想 代数精度 牛顿-科特斯公式 等距离节点的求积公式 n=1 梯形公式 代数精度为1 n=2 Simpson公式 代数精度为3 n=4 Cotes公式 代数精度为5 偶阶求积公式的代数精度
数值分析期末复习知识点
第一章(有效数字位数)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*邻近局部收敛。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Chapter 1 误差误差限计算、有效数字分析•绝对课差址t洵准确俏”*为工的-个近似偵「称T —工対近似偵.T '的絶村谋差,简厳供邛*可简记为E.|g(T)|=| T —*|兰£(/)数值貞门称为T的11绐对误差限或误差限*l『*、F(x ) x —x E© ) = —=——为近似值/的担zt溟誉可简{己址•有效数字若才作加的近tilt其鲍对误差的绝对值不超过某一位数字的半个单恆,而该位数字到F的第—位非零数字共有斤位關称用F近恤时具有血有效做字'简称丫有畀位有效数字.Chapter 2插值法差值条件(唯一性)1、拉格朗日差值a) 插值基函数b) 差值余项2.2拉格朗曰抽值2.2.1基函数考虑最简单、晟舉本的骼值问起+ 求押次插值家项式『低)…肋,便加滿足播值条伸可知,除斗点外.其余都星”.巧的零点■械可诛< (A) ^.4(X 一%[…(-V址 d 為"* <A -A;)X)=A(X - J- (A- - \_, )(.Y -J)其中M为常數.由&工戶1町得』=-------------------- -----------------(閔円)心7冷K%-咖卜-a -斗)和対讼>:T^V为准确血"为玄的一个近似伉称relativeerror称之为拉厳朗LI垒曲绘都是M次帝项武.. 2.1.2拉榕朗n插佢雾项式利用拉辭朗H皋啦数/态人构造次数不趙过"的雾项式£(巧二必机朗+^( v) + •…I J;/,(.v) = £昭(曰可知其搆足7韩为拉格阴Id插说饕砂式.再由插菽牟嘶的唯亠杵“ 鲁 D I特别地*造时又叫钱件擂僮其几何童又为过两点的直级-当*匸2时又叫拋物<线)掩值•具几何鳶义为过三点的拋物线.滾丘阖淘若取人1).伸伏=札1*…飒由插痕参项式的唯一性有£址工)# =x\ k= 0」厂』特别当k-OfiL就得到£佃-1□则铉格朗U的丄抚抽值雾项式为V)= j^(j(X> + I'Jj (x> + j/2(.v) * MQO=(2)弓…仗扣讪—协-町H^)xll(A + l)(r-JX^ 4}+3x —(x H)(x-LXx-3) 8 15■裁1M T-3X V-4)+^X HX A-1M A4)+ l(.v+lX.v-lXr-3)+ 3)a 1已知$ =五,耳=4眄=S.用皴件插值f即一次插惟藝坝如历的近似值.解片=2・曲=3•菇函数付别为:t-9 1 x-4 I4(J)=——=—(x-9j, Zjx)=——= -{x -4)砂14-9 5尸门9-4 5播債孝项式为V)-片fj.i) +」'占(巧-2x^(.v 夕”:(* 4)---(.V 4 J -4)(- (X + fr))所以乔金厶⑺二空R点5使2求过啟-1,-毎川』人(乱-创*(4」)的抛物线播值(即三次插値务项式).蔦-U 斗=-t t A|二L x2=3»A3- 4以为苗点加墓函.数分别为:厶何」匸迪住1±J (.r +lXA -3}(x-4)1(1 ► 1)(1-3)(1- 4J 12心)」:十汽-1年¥二Uw心一ncz (34-1X3-1X3-4) K=⑴】心-叭7= *十叫讣7】(4 + IX4-1X4-3) 152.23極値肇项M tt'r滾^Ji n(x)=f(x)兀糾也称为"次1川甘"叱插伯赛境式的余坝。
以下为脸幡朗日冬顼進理"闻 $ 设川沖=-/tn点七二2».V| = 2r5.-v; -'l.求f(x)定理2设/何在区问S期上禅在卄1阶导暫fflttWiffirjJi aj.fi il 11/⑶師诉組怖并怙汁谡熱x t G風切(QM…+* ■)为沪1F互异if点*则对住何M 片=/Ul- WAji - = /U) -Q.25.诈眩』]*有插恢窑项式为Z = Rv) - L(.v) = 歐“的3「%(2 Z5M1 41 (2.5 邛"4)+ 2弘僅际询n(4-2M*-tS>具中啦M町=11(击-和«石2上)H.与k有黄)J-4=U r<ISr1(L42&T41.15(11.2? KIWIU?坏.ln .…b . flt25-]rtM11J5-IU+" ' * {1J-ID)(12-11)1.4M26MHllOjahnx 的二阶羽淌1卿沁血%(!】站片屿 |(Jl.*S-!«)1lU5 LIXLL£S-l2)K^Mri5实脈上」讥1』5吃剔)36 &|Rt (1L 25) M. 000056.V J (A > (I.4KI75+1.116q i -(L4ri)-t (X2R (HXx - n.4fl)(x-0.55)植/<0L S^6J * v (fts»6) ariiioioXf\^t x t ^x s ]-<iA970吋再竝前四点的三汝牛帧抽位痒琐式.V-, | A-) = .V ; (JU) I U J 177O(.r-Q J 4UXJV甜)(*一口-各力枚七| - H.6JWH5/1.切*'-、乜1可得心(A 〉的:施断谋莖\Jt s (.x )| t |IMB44(A' - 0-0.55)(A -IL 础R H - OLflO>|X, 10.59ti *| = 03 = 10」2、牛顿插值构造差商表TSy ⑶斗⑶"3圧因为 =U. = max | r (x)冃厂⑵ | £ x唄网f*故 | 也)恰善 |(x-2X-V-2.5X^-4) |<;^|(X -2H X -25X A -4)|I R ⑶ H f(3)-L(3)l^ 1 11(3-2X3-23X3 fl IQ S= 0.03125性质」若爪)倒砧]上存在刊阶导数且节点5曲”屮 斗EM 』]捌至少存在•点胆S 』】满足下式严忆)八辛朋讨,片]片号」例1 /■何=—6円加一地求/[1J …”9]及/1W ,…,10].例」蛤屯函數表 用二氏厨fl OlulL 药的诉他血并估计谣*.例]比尸/UAH/vR’ 1品盒九乱工用三试插佰率项 式求/U ②玻/U 曲的圻似値.斛相应鸽函散値及逐分表如下:求川用悴頓插公成*且由弗須A;(1.2}=171N2U+ l.TdMJ x(L4+U密iU利〕4-1)求皿成湘斗顿牯插公式山由人S=JTL5f*得匸4U 用逼库再⑴}|求fl由邮?]一—* £10962 .-2UU«5SJb7.!XJJU5M + —, -( (14)x( <14 i 1) [•22356 A +—(-0.4)x(-ft.4+1X^+2) -lS76aQS713*r3、埃尔米特插值构造三次埃尔米特插值多项式如下2.5坝尔来特(Hermite)插ffl2.S.1三狀埃尔米特插值寒項式设J弓懐)是区间他“]上的实函数,.畑山屣血列上相昴吋点,且比VgJ T/C V)也J L的曲數值HJ ―齢諭数値分别为FWfSJ 利舰严广(对(HM K求三次峯项式J/J-x),使其满足=I H 心J = y tJ t OHMX町称为三次埃尔;3^特插Ifi爭項式.定砰3満足条祥式丹£儿)=片.丹;(叫)=叫0 = 0.巧的三次娓尔米特插值多项式存在fl喉一”构倉三抉埃尔米特插值多项式如下;兀©"儿啊(时+ y^i('>+%A(v)*叫舛W由flf.Qr"-%(.、•[ ) = 0可将卩;写成% (算)= [a+b(x-切]"-屯Frtl <X(J(A'o) = lr '得口■=飞凤--VJ再由ff;(-V,) = 0T得* = ----- -- •所以(#。
一科)K O(.V)=|1+2^丑]•甘厂斗屿-Jfj同理I将也< >A\)盘心2[中0](二^曲一%曲一&同样-o,可令舛(絆=创大一科]X 片一 A \)-再由 /^(^) — 1* c —'C 點(> —'升)A (小=(兀-叭X "旳)・ 冷-V J角(划=1工—巧X 竺五尸H 厂*0曲一眄用盘)-U A?X1)血才[中4](二几*一叼 Jj -A t闻町={.、叼)(巧一观即^(x)=(l >l^(xM(jr)A<Jf) = (-K-X,rf(*)A (x> = U l -A l K I (x)* !>,/, ( v)为 U ( •”*•)・(v,,y t )插直点的T 吗mnftp—次基憐数.可得滿足条件的三钦城尔耒特播備多项式为H 侶= J^(.V )+ 帕何 4 i^A (X )+ 叫恥)5當強5中駛*4、分段低次插值15.2蟆差佔计定评4颐舒在血含局・卫的区间内存在四附 导議,时肴余项 ^)=f (x )妆(Q 二吕严运臼-阳愉-就(紀(以)口与侑为 设叫=黑二严(蝌则当涎鵬曲财 余頂有如下怙计式g 蛙限》&⑴生色孑丨貞A JS4 科2己知血戶怦及艮一■阶导数的Off.下表弼埃尔 来特捕值公式 013511的近似值•并佔讣叹載斷误差* 解瓦的■市(1-71巧匕-144)壮盯+ ———,(Jt-12!}(.¥-14Jf +——-_ Li - J44W.v-lil¥ I : 23 x A J 24 23 ' " f— 121、 im 」• r I I -V —144 '+12 M 1 + 2 -------- I 121-144J1 ( .r-UlWJ.Y) =llx 14-2A in 144u 12 rm 1/21 1/24 + -X IIr- L44121 144y-12l 144-121.t- 1-M Y 144-121 AIM' 144J x-121 V )44-121J疋“绚| =而石忌衣了 乞 1519- a : o.oooo n384 1113 >11[叫(用 A n X X/ * 科一叫I- ' 构造函SW-v)-1avrti 1 SX'SI0上的数表、应如何选取步伽才能便利用数表进行分段插值时逞差不超过心5X1解f H(x)二一丄,M 二tUA\/*(x)| - 1,耳. 1牡切1 1欲使I ^加_ P(X)|玄y醉|/*3|吟兮曲旷得A <2x10^即进行分段线性插值时,应取A0X1U%课垦不超过0.5X11H。