数值分析复习要点201312
数值分析复习提纲(修改完)
第一章 绪论【考点1】绝对误差概念。
近似数的绝对误差(误差):()a =x a E -,如果()δa E ≤则称δ为a 的绝对误差限(误差限)。
【考点2】相对误差限的概念。
近似数a 的相对误差:()()/x a x =a E r -,实际运算()()/a a x a E r -=,a r /δδ=。
【考点3】有效数字定义。
设*x 的近似值a 可表示为n m a a .a a= 21010⨯±,m 为整数,其中1a 是1到9中的一个整数,n a a 2为0到9中的任意整数,若使()n m a||=|x a |E -*⨯≤-1021成立,则a 称近似*x 有位有效数字。
例:设256010002560,00256702.×=.a .=x -*=,则4-10×21=0.00005a -x ≤*。
因为,2-m=所以2n=,a 有2位有效数字。
若257.01000257.02⨯==-a ,则5102100000500000030-≤×=..=x-a ,因为2-=m ,所以3=n ,a 有3位有效数字。
例:设000018.x=,则00008.a=具有五位有效数字。
41021000010-≤×.=x-a ,因为1=m ,所以5=n ,即a 具有五位有效数字。
例:若3587.64=x *是x 的具有六位有效数字的近似值,求x 的绝对误差限。
410×0.358764=x *,即4=m ,6=n ,0.005=1021x -x 6-4⨯≤*【考点4】四舍五入后得到的近似数,从第一位非零数开始直到末位,有几位就称该近似数有几位有效数字。
【考点5】有效数字与相对误差的关系。
设x 的近似数为n m a a .a ×a= 21010±,)(a 01≠如果a 具有n 位有效数字,则的相对误差限为()111021--≤n r ×a δ,反之,若a 的相对误差限为()()1110121--+≤n r ×a δ,则a 至少具有n 位有效数字。
数值分析主要知识点
第三章
非线性方程的数值解法
二分法的思想以及其中对分次数的计算;
不动点迭代法、迭代格式的收敛性判定方法、
误差估计式;
Newton迭代法及其收敛性; 割线法迭代格式;
迭代加速方法。
第四章
线性方程组的直接解法
Gauss消去法与列主元素Gauss消去法; 三角分解(LU)法; 平方根方法(Cholesky分解); 向量与矩阵范数; 条件数与病态方程组求解。
第五章
曲线拟合与最小二乘问题
拟合与插值的异同点、矛盾方程组的最小二乘解; 满秩分解、法方程组、可化为线性拟合的非线性拟合;
(极小)最小二乘解的存在唯一性、广义逆与极小
最小二乘解;
GS与MGS正交化与最小二乘解;
Householder正交化与最小二乘解。
第六章代法与Gauss-Seidel迭代法及其收敛性;
SOR迭代法及其收敛的必要条件、最佳松弛因子; 解非线性方程组的Newton迭代法与拟Newton思想。
第七章
最优化方法与共轭梯度法
与方程组等价的变分问题、线性寻查(线搜索)法;
最速下降法; 解线性方程组的共轭梯度法。
写、不得打印、不得复印,纸上签有姓名和学号;
可以携带计算器(考试期间不允许互借)。
《数值分析》复习主要知识点 第一章
绪论 基本概念:误差的分类(截断误差、舍入误差)、 绝对误差和相对误差、有效数字;
数值稳定性; 误差分析的原则:1)尽量避免相近的数相减,2)
尽量避免绝对值小的数做除数,3)防止大数吃小数, 4)先化简再计算,5)选用数值稳定的算法;
浮点数系统特征(四个整数表征)。
第八章
数值微分与数值积分
数值分析期末知识点总结
数值分析期末知识点总结一、引言数值分析是一门研究如何使用计算机提高数学模型数值计算精度和效率的学科。
它是计算数学的一个重要分支,涉及到数值计算、数值逼近和误差分析等一系列内容。
在数值分析课程中,我们将学习到数值解微分方程、线性代数问题的求解、插值与拟合、积分等一系列内容。
本文将对数值分析期末知识点进行总结,以便帮助大家复习。
二、常见数值计算方法1. 插值与拟合插值与拟合是数值分析中重要的内容,它们用于在给定数据点集上构造一个函数,以便在其他点上进行求值。
插值是通过一些已知数据点来求得一个函数,使得这个函数能够通过这些点,而拟合则是通过已知数据点来求得一个函数,使得这个函数在这些点附近能够比较好地拟合数据。
常见的插值方法包括线性插值、拉格朗日插值、牛顿插值等;而拟合方法包括最小二乘法拟合、多项式拟合等。
2. 数值解微分方程数值解微分方程是数值分析的一个重要内容,它讨论如何使用计算机对微分方程进行数值求解。
微分方程是自然界中描述变化的数学方程,它们在物理学、化学、生物学等领域都有着重要的应用。
数值解微分方程的方法包括欧拉法、中点法、四阶龙格-库塔法等。
3. 数值线性代数数值线性代数是数值分析领域的另一个重要内容,它讨论如何使用数值方法解决线性代数问题。
原始的线性代数问题可能非常大或者非常复杂,因此我们常常需要使用计算机进行数值计算。
数值线性代数的方法包括高斯消元法、LU分解、Jacobi迭代法、Gauss-Seidel 迭代法等。
4. 数值积分数值积分是数值分析的一个重要内容,它讨论如何使用数值方法对积分进行数值求解。
在实际问题中,有很多积分问题是无法解析求解的,因此我们需要使用数值方法进行近似求解。
数值积分的方法包括复合辛普森法、复合梯形法、龙贝格积分法等。
三、数值分析的误差分析在数值计算过程中,我们会遇到误差的问题。
这些误差可能来自于测量、舍入、截断等各种原因。
因此,误差分析是数值分析中一个非常重要的内容。
(整理)《数值分析》期末复习纲要.
《数值分析》期末复习纲要 第一章 数值计算中的误差分析主要内容(一)误差分析 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. 知道产生误差的主要来源。
2. 了解绝对误差和绝对误差限、相对误差和相对误差限和有效数字等概念以及它们之间的关系。
3. 知道四则运算中的误差传播公式。
4. 避免误差危害的若干原则第二章插值法(一) 考核知识点插值函数,插值多项式,被插值函数,节点;拉格朗日插值多项式:插值基函数;均差及其性质,牛顿插值多项式;分段线性插值、线性插值基函数。
(二)复习要求1. 了解插值函数,插值节点等概念。
2. 熟练掌握拉格朗日插值多项式的公式,知道拉格朗日插值多项式余项。
3. 掌握牛顿插值多项式的公式,了解均差概念和性质,掌握均差表的计算,知道牛顿插值多项式的余项。
4. 掌握分段线性插值的方法和线性插值基函数的构造。
第三章函数逼近(一) 考核知识点函数逼近的基本概念,内积,范数,勒让德与切比雪夫正交多项式,最佳一次一致逼近,最佳平方逼近,曲线拟合的最小二乘法(二)复习要求1. 熟练掌握内积,范数等基本概念。
2. 熟练掌握勒让德与切比雪夫正交多项式的性质。
3. 掌握用多项式做最佳平方逼近的方法。
4. 最小二乘法及其计算方法。
第四章数值积分与数值微分(一) 考核知识点数值求积公式,求积节点,求积系数,代数精度;插值型求积公式,牛顿―科特斯求积公式,牛顿―科特斯系数及其性质,(复合)梯形求积公式,(复合)Simpson求积公式;高斯型求积公式,高斯点,(二点、三点)高斯―勒让德求积公式;(二) 复习要求1. 熟练掌握数值积分和代数精度等基本概念。
2. 熟练掌握牛顿−科特斯求积公式和科特斯系数的性质。
熟练掌握并推导(复合)梯形求积公式和(复合)Simpson求积公式。
3. 知道高斯求积公式和高斯点概念。
会用高斯−勒让德求积公式求定积分的近似值。
(完整版)数值分析考试复习总结汇总,推荐文档
10
100
误差估计:
f
max | f (x) fh (x) |
(x ih) (x (i 1)h) . 2! ixx(i1)h
□
第三章
最佳一致逼近:(了解) 最佳平方逼近
主要分两种情形:
1. 连续意义下
在空间 L2[a,b]中讨论
2. 离散意义下
在 n 维欧氏空间 Rn 中讨论,只要求提供 f 的样本值
n (x)
(x
xi
)
n
(xi
)
ji
n
n
其中: n (x) (x x j ), n xi (xi x j ) .
j0
j0
ji
例 1 n=1 时,线性插值公式
P1 ( x)
y0
(x x1) (x0 x1)
y1
(x x0 ) (x1 x0 )
,
例 2 n=2 时,抛物插值公式
P2 (x)
可得: L3 (x) x 2 (x 1 2)
方法二. 令
L3 (x) x(x 1 2) ( Ax B)
由
L3
(1)
3 2
,
L3 (1)
1, 2
定 A,B
(称之为待定系数法)
□
15.设 f (x) x2 ,求 f (x) 在区间[0,1] 上的分段线性插值函数 fh (x) ,并估计误差, 取等距节点,且 h 1/10 .
(2)
2x ( x 1 x
x 1 x) .
(3) 1 cos x sin 2 x sin x .
□
x
x(1 cos x) 1 cos x
第二章
拉格朗日插值公式(即公式(1))
数值分析复习要点
1.设矩阵A
2
1
,
用Schmidt正交化方法,
1 2
对A作正交分解A QR.
2.设矩阵A
2 0
1 3
7 10
,
用Householder变换法,
0 4 5
对A作正交分解A QR.
3.已知一组线性无关的向量
u1 (1,1,1)T , u2 (1, 0, 1)T , u1 (0,1,1)T , 由此向量组,按Schmidt正交化方法,求一组对应的
Gauss变换阵
1
Lj
1
l j1, j 1
ln, j
1
对x x1,..., x j ,..., xn T 0, x j 0 构造Gauss变换阵G,使Gx x1,..., x j ,0,...,0 T
奇异值与奇异值矩阵
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) min ( AT A)
y0
1 yn n 5 yn1
n 1, 2,...
计算yn,试分析算法的稳定性
习题:p15 10
数值计算中应注意的问题
(1) 防止相近的两数相减 (2) 防止大数吃小数 (3) 防止接近零的数做除数 (4) 注意计算步骤的简化,减小运算次数
数值分析期末复习要点总结省公开课获奖课件市赛课比赛一等奖课件
15
Lagrange插值
Lagrange插值基函数
设 lk(x) 是 n 次多项式,在插值节点 x0 , x1 , … , xn 上满足
1, j k lk ( x j ) 0, j k
则称 lk(x) 为节点 x0 , x1 , … , xn 上旳拉格朗日插值基函数
16
线性与抛物线插值
两种特殊情形
x0 ƒ(x0)
x1 ƒ(x1) ƒ[x0, x1]
x2 ƒ(x2) ƒ[x1, x2] ƒ[x0, x1, x2]
x3 ƒ(x3) ƒ[x2, x3] ƒ[x1, x2, x3] ƒ[x0, x1, x2, x3]
…
xn ƒ(xn) ƒ[xn-1, ƒ[xn-2, xn-1, ƒ[xn-3, xn-2, xn-1, … ƒ[x0, x1,2…7 ,
ln 0.54 旳精确值为:-0.616186···
可见,抛物线插值旳精度比线性插值要高
Lagrange插值多项式简朴以便,只要取定节点就可写 出基函数,进而得到插值多项式,易于计算机实现。
19
Lagrange插值
lk(x) 旳体现式 由构造法可得
lk (
x)
( x x0 ) ( xk x0 )
Rn(x)
n1
Nn( x) a0 a1( x x0 ) a2( x x0 )( x x1 ) an ( x xi )
i 1
其中 a0 f ( x0 ), ai f [x0 ,, xi ], i 1,2,, n
Nn(x) 是 n 次多项式
Rn( x) f [x, x0 , ... , xn]( x x0 )...( x xn1)( x xn )
若
e(x*) x x*
数值分析2012总结
x ( k +1) = Mx ( k ) + g
(k=0,1,2,⋯ )
产生向量序列{ x ( k ) },当k充分大时,以 x ( k )作为 方程组Ax = b的近似解,这就是求解线性方程组 的单步定常线性迭代法。 M 称为迭代矩阵。
收敛问题:
(Jacobi) 迭代法 雅可比 雅可比(Jacobi) (Jacobi)迭代法 高斯 —塞德尔迭代法 高斯—
� b 有扰动,A 无扰动
|| δ x || || δ b || ≤ (|| A |||| A−1 ||) || x || || b ||
� A 有扰动,b 有扰动
|| δ x ||≤|| ( I + A−1δ A)−1 |||| A−1 |||| δ A |||| x || || δ x || ≤ ( I + A−1δ A) −1 A−1 || x || || A−1 |||| A || || A−1 |||| δ A || δA ≤ 1− || A−1δ A ||
∞
A− 1 A
−1
∞
2
2
λmax ( AT A) = λmin ( AT A)
结论: 当条件数很大时 ,方程组 Ax = b 是病态问题 ;当条件数较小时 , 当条件数很大时, 是病态问题; 当条件数较小时, 方程组 Ax = b 是良态问题。
� || Ax || p � || A || p = max = max || Ax || p � � � � x≠0 || x || = 1 p || x || p
特别有: || A ||∞ = max ∑ | aij | (行和范数) 1≤ i ≤ n
j =1
n
n
|| A ||1 = max ∑| aij |
数值分析期末复习总结(优选.)
线性插值多项式(一次插值多项式)
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、四舍五入的数如何求有效数字?已知绝对 误差的数如何求有效数字?
2、常用向量、矩阵范数的计算 3、一般范数的定义,如何证明一个函数为范数? 4、数值计算中注意的五个原则
二、线性方程组的解法
1、Gauss消去法(顺序、列主元) 2、矩阵的条件数 3、常用的三种迭代法及其收敛条件
Gauss型求积公式. 第一步:找高斯点
构造[a,b]上带权 (x) 的n次正交多项式,并求其
零点 x1, x2, , xn 作为高斯点。
由于首项系数并不影响正交性,不妨把首项系数均定为1.
1)待定系数法
设 g0 (x) 1, g1(x) x a, g2 (x) x2 bx c,
由正交性确定待定系数a,b,c,…... 2)利用递推公式
在 g0 (x), g1(x) 已知的情况下,才能用递推公式.
第二步:确定求积系数 Ak
4、 Gauss型求积公式的构造
第一步:找高斯点
第二步: 确定求积系数: Ak
1)解线性方程组
分别取 f (x) 1, x, x2, , xn1 使积分公式精确成立:
y f (t, y), t0 t T
y(t0
)
y0
(1)基于数值积分方法
(7.1) (7.2)
y(t)dt f (t, y(t))dt
解决两个问题:
1、积分上下限的选择 2、右端被积函数的近似
y(t)dt f (t, y(t))dt
解决两个问题:
H5(xi ) f (xi )
P158
H5(x) p3(x) q1(x)4(x)
数值分析期末复习要点总结
数值分析期末复习要点总结数值分析是一门研究用数值方法来解决数学问题和科学工程问题的学科。
它包括数值计算、数值逼近、数值求解以及数值模拟等内容。
本文将从数值计算的基础知识、数值逼近方法、数值求解方法以及数值模拟方法等方面进行复习要点总结。
一、数值计算的基础知识1. 计算误差:绝对误差、相对误差、有效数字、舍入误差等等。
2. 机器精度:机器数、舍入误差、截断误差等等。
3. 数值稳定性:条件数、病态问题等等。
4. 误差分析:前向误差分析、后向误差分析等等。
二、数值逼近方法1. 插值方法:拉格朗日插值、Newton插值、Hermite插值等等。
2. 曲线拟合:最小二乘法、Chebyshev逼近等等。
3. 数值微分:前向差分、后向差分、中心差分等等。
4. 数值积分:梯形法则、Simpson法则等等。
三、数值求解方法1. 非线性方程求解:二分法、牛顿迭代法、弦截法等等。
2. 线性方程组求解:直接法(Gauss消元法、LU分解法)和迭代法(Jacobi法、Gauss-Seidel法)。
3. 特征值和特征向量:幂法、反幂法、QR分解法等等。
4. 非线性最优化问题:牛顿法、拟牛顿法、梯度下降法等等。
四、数值模拟方法1. 常微分方程数值解法:Euler法、改进Euler法、Runge-Kutta法等等。
2. 偏微分方程数值解法:差分法、有限元法、有限差分法等等。
3. 数值优化方法:线性规划、非线性规划、整数规划等等。
五、数值计算软件1. MATLAB基础:向量、矩阵、符号计算等等。
2. MATLAB数值计算工具箱:插值与拟合工具箱、符号计算工具箱等等。
3. 其他数值计算软件:Python、R、Octave等等。
总结数值分析是一门重要的数学学科,它为解决实际问题提供了有效的数值方法。
在数值计算的基础知识中,我们需要了解计算误差、机器精度和数值稳定性等概念,同时也需要掌握误差分析的方法。
数值逼近方法包括插值、曲线拟合、数值微分和数值积分等内容,其中插值和拟合是常见的逼近方法。
《数值分析复习提纲》word版
数值分析第一部分线性方程组的数值解法一、基本要求1、掌握每一种解法的基本思想,适用范围,收敛条件,计算公式以及误差估计.2、在应用中不同解法的异同、优劣,加深对算法的理解,最好能上机计算.二、主要概念及结果主要概念定义1.1 对于方程通过某种方法建立了迭代法(2.1.1)如果对于任何使得极限成立,则称该迭代法是收敛的.定义1.2 如果,对于,都有成立,则称A是严格对角占优的.主要算法与定理高斯(Gauss)消去法假设A的所有顺序主子式都不等于零,原来的方程组为计算步骤为1) 把上面的第一个方程除以,在分别乘上后与第k 个方程相加(),得到于是我们从第2到第n 个方程中消去了.2) 把上面的第二个方程除以,再分别乘上后与第k 个方程相加()得到于是我们从第3到第n 个方程中消去了.3) 继续这个过程直到我们得到4) 由上面的最后一个方程很容易得到,然后按相反次序回代逐一计算出方程的解.高斯(Gauss)列主元消去法 假设A 的所有顺序主子式都不等于零,原来的方程组为(1) 消元过程.对,进行以下运算: 1) 选主元.找行号,使得; 2) 交换中的ki k ,两行;3) 消元:对于; 对.(2) 回代过程.按下述公式;回代求解即可得到方程组的解.定理1.1 对于,如果A 的所有顺序主子式都不为零,则存在唯一的上三角矩阵U 和对角元素为1的下三角矩阵L,使得Doolittle 分解 根据定理1.1,对于,如果A 的所有顺序主子式都不为零,则存在唯一的上三角矩阵U 和对角元素为1的下三角矩阵L,使得.可以直接计算分解式中的诸元素.为此,我们假设⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=-11111,21323121n n n n l l l l l l L,⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=-----nn n n n n n n n n u u u u u u u u u u U ,11,121,22211,11211用U 的第k 列()乘L ,然后与A 的相应列比较,可以逐列(逐行)计算出L(U)的元素.定理1.2 设A 是一个对称正定矩阵,则存在唯一的下三角阵L ,其对角元素都是正的,使得定理1.3 设A 是一个对称正定矩阵,则存在一个单位下三角阵L和对角矩阵D,使得定理1.4 迭代法对于任意收敛的充分必要条件是,其中是迭代矩阵的谱半径.如果及假设A的对角元素,令A=D-L-U,其中D是A 的对角部分构成的矩阵.L和U分别是A的严格下(上)三角矩阵,则有以下几个具体算法:雅可比迭代法高斯-赛德尔迭代法关于这两个算法的收敛性有如下定理:定理1.5 如果方程组Ax=b的系数矩阵是严格对角占优的,则雅可比迭代法和高斯-赛德尔迭代法都收敛.定理1.6 如果方程组Ax=b的系数矩阵是对称正定的,则高斯-赛德尔迭代法收敛.第二部分非线性方程的数值解法一、基本要求掌握每种方法的基本思想、迭代公式、收敛条件以及与其他方法的差异.二、主要概念及结果主要概念定义2.1 对于方程,通过某种方法建立了迭代法(2.1)如果存在使得极限,则称该迭代法是收敛的.主要算法与定理定理2.1 设有方程,如迭代函数在有根区间[a,b]上满足:(1)当时,;(2)在[a,b]上可导,且有,则有:(1)方程在[a,b]上有唯一的根*x;(2)对任意初值,迭代公式产生的数列收敛于方程的唯一根*x,即;(3)误差估计定理2.2 设*x是方程的根,在*x的某个邻域内连续,且有,则必存在*x的一个邻域,对于任意选取的初值,迭代公式产生的数列收敛于方程的根*x.二分法假设的隔根区间为,取,计算.如果,则取,否则取.继续这个过程直到取见足够的小,就可以把最后区间的中点作为方程的近似根.此法称为二分法.牛顿法计算公式定理 2.3 如果,且在*x的某个邻域内连续,则牛顿法是局部收敛的.弦截法计算公式第三部分插值法一、基本要求1、在算法上要求熟练掌握拉格朗日插值法,等距节点插值法,牛顿插值法.2、要求能按所给条件,选用适当的近似公式求出近似函数或计算出函数的近似值,并会估计其误差.二、主要概念及结果主要概念定义3.1 设在区间上有定义,且在上的个不同的点的函数值为,若存在一个代数多项式(3.1)其中为实数,使得成立,则称为函数的插值多项式,点称为插值节点.主要算法与定理定理3.1 在个互异节点上满足插值条件的次数不高于的插值多项式存在且唯一.拉格朗日插值多项式的一般形式 其中为插值基函数, 插值余项为其中是区间中的某一个值,且和x 有关,所以牛顿插值多项式及余项)())(](,,,[))(](,,[)](,[)()(11010102100100----++--+-+=n n n x x x x x x x x x f x x x x x x x f x x x x f x f x N余项牛顿前插公式牛顿后插公式第四部分数值积分与数值微分一、基本要求掌握梯形求积公式、辛普森求积公式以及复化的梯形公式、复化的辛普森公式和龙贝格公式的构造方法.二、主要概念及结果主要概念定义4.1 若求积公式对于任意不高于次的代数多项式都准确成立,而对于次多项式却不能准确成立,则称该求积公式具有次代数精度.定义 4.2 将个节点的具有次代数精度的插值型求积公式称为高斯型求积公式,节点称为高斯点,称为高斯系数.主要算法与定理插值型求积公式其中牛顿-柯特斯公式其中梯形公式辛普森公式柯特斯公式其中复化梯形公式复化辛普森公式复化柯特斯公式其中龙贝格求积公式定理4.1 节点为高斯点的充分必要条件是以这些点为零点的多项式与任意次数不大于的多项式在上正交,即.第五部分常微分方程的数值解法一、基本要求掌握欧拉公式、经典的龙格-库塔公式二、主要概念及结果主要算法和定理显式欧拉方法隐式欧拉方法梯形公式预报-校正方法预估校正龙格-库塔方法二阶龙格-库塔公式经典的四阶龙格-库塔公式。
数值分析复习总结
数值分析复习总结数值分析课本重点知识点第一章P4定义一P5定义二P6定理1P7例题3P10条件数(1)绝对误差(限)和相对误差(限)公式(2)有效数字(3)条件数及其公式第二章P26定理2(以及余项推导过程)P36两个典型的埃尔米特插值(1)拉格朗日插值多项式(包括其直线公式和抛物线公式)(2)插值余项推导及误差分析(估计)(3)两个典型的埃尔米特插值(4)三次样条插值的概念第三章P63例题3(1)最佳平方逼近公式的计算(2)T3(x)的表达式第四章P106复合梯形公式P107复合辛普森求积公式P108例题3(1)复合公式及其余项(2)判断一个代数的精确度第五章P162定义3向量的范数P165定理17P169定义8(1)左中右矩形公式(2)LU分解(3)谱半径和条件数(4)向量的范数第六章P192定理9第1条P192例题8第七章P215不动点和不动点迭代法P218定理3P228弦截法P229定理6第九章P280欧拉法与后退欧拉法P283改进欧拉公式数值分析课后点题答案第一章数值分析误差第二章插值法第三章函数逼近所以无解19。
观测物体的直线运动,得出以下数据:时间t(s) 0 0.9 1.9 3.0 3.9 5.0 距离s(m)10305080110求运动方程。
解:被观测物体的运动距离与运动时间大体为线性函数关系,从而选择线性方程 s a bt =+ 令{}1,span t Φ=22012201016,53.63,(,)14.7,(,)280,(,)1078,s s =====则法方程组为614.728014.753.631078a b = ??? ?从而解得7.85504822.25376a b =-??=? 故物体运动方程为22.253767.855048S t =-20。
已知实验数据如下:i x 19 25 31 38 44 j y19.032.349.073.397.8用最小二乘法求形如2s a bx =+的经验公式,并计算均方误差。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
caik sail , saik cail ,
bij aij ,
1i n 1i n j k,l
设 A (aij ) , B J (k,l) A (bij ) ,便有:
bkj cakj salj , blj sakj calj , bij aij ,
1 j n 1 j n i, j k,l
变换法, 对A作正交分解A QR.
数值分析复习要点
一. 基本概念
二. Gauss变换与矩阵的三角分解 三. Householder变换 四. Givens变换 五.矩阵的正交分解 六.正交相似化简 七.解线性方程组Ax=b的直接法和迭代法 八. 构造正交多项式
九. 连续函数的最佳平方逼近 十. 离散数据的最小二乘曲线拟合 十一. 函数插值 十二. 数值积分 十三. 数值微分 十四.非线性方程的数值解法 十五. 常微分方程的数值解法 十六.数值计算的基本思想 十七.读程序 ,写程序
五.矩阵的正交分解
(1)Schmidt正交化法
(2) Householder变换法
习题 :
(3) Givens变换法
1 2
1.设矩阵A
2
1
,
用Schmidt正交化方法,
1 2
对A作正交分解A QR.
2 1 7
2.设矩阵A
0
3
10
,
用Householder变换法或Givens
0 4 5
Lj , 其中li, j
xi xj
,
i j 1, j 2,...n
LU分解
习题 : 1.设x (2,1, 1, 3)T ,求一Gauss变换阵L,
使Lx (2, 0, 0, 0)T .
2.已知矩阵
1 2 3
A
2
6
2
,
3 1 5
对矩阵A作三角分解,即A LU .
三ห้องสมุดไป่ตู้ Householder变换
Householder变换阵 H I 2wwT ,其中|| w ||2 1
定理 : 设n维向量x, y, x y, 但 || x ||2 || y ||2 , u x y, 则存在Householder变换阵 H I 2wwT , w u ,
|| u ||2 使Hx y.
习题
1.已知向量x (2, 0, 2,1)T , 试构造Householder阵H
求x, y之间的距离( x, y).
返回
二. Gauss变换与矩阵的三角分解
Gauss变换阵
1
1
Lj
l j1, j 1
ln, j
1
对x
T
x1,..., x j ,..., xn 0,
xj 0
构造Gauss变换阵G,使Gx
T
x1,..., x j , 0,..., 0
解:G
使Hx ke3,其中e3 0, 0,1, 0T , k R.
2.已知向量x (1, 2,1, 2)T , 试构造Householder阵H
使Hx (1, 2, 0, 0)T .
3.已知向量x (1, 2, 2)T , y (0,3, 4)T , 试构造一个 Householder阵H使Hx ky,k R。给出k值和 变换阵H .
1
(
b
|
f
( x)
g(x) |p
dx) p
a
p 1, 2, ,
p2627
习题 :
3 2
1.已知矩阵A
2
3
1 0
(1) A的谱半径( A),
1 0 , 试计算 3
(2) A的谱条件数cond ( A)2, (3) A的奇异值
2.已知向量x (1, 4, 3, 0)T , y (3, 6,1, 2)T ,
距离概念
向量空间的距离
n
1
(x, y) || x y ||p ( | xi yi |p ) p
i 1
矩阵空间的距离
( A, B) || A B ||p p 1, 2, , F
p 1, 2, , p25
连续函数空间的距离
( f ( x), g( x)) || f ( x) g( x) ||p
xk2 xl2 r
yk cxk sxl r xk2 xl2 ,
yl 0
y J(k,l)x
=
x1, ..., xk-1, r, xk+1, ..., xl-1, 0, xl+1, ..., xn
T
.
A (aij ) , B AJ T (k, l) (bij ) ,则有:
bbiikl
一. 基本概念
绝对误差,相对误差,有效数字,数值稳定性等. 1、设x和y的相对误差为0.001,则x*y的相对误差约为
____________.
绝对误差,相对误差,有效数字,数值稳定性等.
3、递推公式
y0
2
,
yn 10 yn1 1, n 1, 2, ...
如果取y0 2 1.41作计算,则计算到y10时, 误差为_______;这个计算公式数值稳定不稳定?
条件数 cond( A)p || A1 ||p|| A ||p ,
p F,1, 2,
谱条件数
cond( A)2 || A1 ||2|| A ||2
max ( AT A) min ( AT A)
1 n
A对称时,
cond( A)2
||
A1
||2||
A
||2
| |
( A) |max ( A) |min
4、计算 y = ln x。若 x 20,则取 x 的几位有效数 字可保证 y 的相对误差 < 0.1% ?
数值计算中应注意的问题
(1) 防止相近的两数相减 (2) 防止大数吃小数 (3) 防止接近零的数做除数 (4) 注意计算步骤的简化,减小运算次数
向量范数
矩阵范数 ①矩阵A的1 范数(又称为A的列范数)
四. Givens变换
定义4-1 n阶方阵
1
1
cos
1
J(k,l, )
sin
称为Givens变换阵。
sin
1
cos
1
k
l
1 (l k)
构造Givens变换阵J(k, l),使y J(k, l)x的第l个 分量为零.
c cos xk xk , s sin xl xl
xk2 xl2 r
n
A
1
max 1 jn
i 1
aij
②矩阵A的 范数(又称为A的行范数)
n
A
max
1i n
j 1
aij
③矩阵A的2 范数(又称为A的谱范数)
A 2
max ( AT A) 1
其中max ( AT A)表示( AT A)的最大特征值,
1为A的最大奇异值.
谱半径 ( A) | ( A) |max