第三章数值分析
(完整版)数值分析答案第三章习题
解:(1)建立假设 H0 : 0.973
n=100, x 2.62 ,s=0.06
(2)在 H0 成立的前提下,构造统计量
u x 0 ~ N (0,1)
s/ n
(3)给定 0.05 ,查得 u 1.96 ,使
2
p u u
2
(4)由样本计算,
x 62.24,s*2 404.77,问这天保险丝融化时间
分散度与通常有无明显差异( 1%)?假定融
化时间是正态母体。
解:(1)建立假设H0: 2 02 400
(2)在H
成立前提下,构造统计量
0
2
(n
1)s*2
02
~
2 (n 1)
(3)给定显著水平 0.01,查得
20.00(5 24) 45.559, 20.99(5 24) 9.886,使
:
2
2 0
1.62
s* 2.296
(2)在 H0 成立的前提下,构造统计量
2 (n 1)s* ~ 2 (n 1) 0
(3)给定 0.05 ,查得
2 (5) 14.49
2
1 2 (5) 1.237 2
(4)由样本计算,
2
6 2.296 2 1.62
12.355
2 (6) 14.49
确,收集到一组使用新安眠药的睡眠时间为
26.7,22.0,24.1,21.0,27.2,25.0,
23.4
试问:从这组数据能否说明新安眠药的睡眠时
间已达到新的疗效(假定睡眠时间服从正态分
布,取 0.05 )?
解:
x1 ~ (20.8,1.62 ) x2 ~ (, 22 )
数值分析 第三章 基于MATLAB的科学计算—线性方程组1概述
科学计算—理论、方法及其基于MATLAB的程序实现与分析 三、 解线性方程组(线性矩阵方程)解线性方程组是科学计算中最常见的问题。
所说的“最常见”有两方面的含义:1) 问题的本身是求解线性方程组;2) 许多问题的求解需要或归结为线性方程组的求解。
关于线性方程组B A x B Ax 1-=⇒=(1)其求解方法有两类:1) 直接法:高斯消去法(Gaussian Elimination ); 2) 间接法:各种迭代法(Iteration )。
1、高斯消去法1) 引例考虑如下(梯形)线性方程组:()⎪⎩⎪⎨⎧==+==+-=⇒⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--⇔⎪⎩⎪⎨⎧==-=+-5.0141315.3221122004301211214322332321321332321x x x x x x x x x x x x x x x 高斯消去法的求解思路:把一般的线性方程组(1)化成(上或下)梯形的形式。
2)高斯消去法——示例考虑如下线性方程组:⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---⇔⎪⎩⎪⎨⎧=++-=-+-=+-306015129101.2001.221113060129501.2001.221321321321321x x x x x x x x x x x x1) 第一个方程的两端乘12加到第二个方程的两端,第一个方程的两端乘-1加到第三个方程的两端,得⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--3060031110001.0001.00111321x x x2) 第二个方程的两端乘001.010-加到第三个方程的两端,得 ⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--60600311010001.0001.00111321x x x3) 从上述方程组的第三个方程依此求解,得()⎪⎩⎪⎨⎧==+-==+-=600300001.03100024011332321x x x x x x 3)高斯消去法的不足及其改进——高斯(全、列)主元素消去法在上例中,由于建模、计算等原因,系数2.001而产生0.0005的误差,实际求解的方程组为⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---306015129101.20005.22111321x x x ⎪⎩⎪⎨⎧===⇒70.4509.30142.2565321x x x 注:数值稳定的算法高斯列主元素消去法就是在消元的每一步选取(列)主元素—一列中绝对值最大的元取做主元素,高斯列主元素消去法是数值稳定的方法。
数值分析第三章
称为1 − 范数 , 称为 2 − 范数 .
(
b 2 ∫a f ( x )dx
),
1 2
三、内积与内积空间
R n中向量x及y定义内积 : ( x, y ) = x1 y1 + L + x n y n .
定义3 上的线性空间, 定义3 设X是数域 K ( R或C)上的线性空间,对 ∀u, v ∈ X, 中一个数与之对应, 并满足条件: 有K中一个数与之对应,记 为( u, v ),并满足条件: (1) ( u,v ) = (v , u), ∀u,v ∈ X ; (2) (αu,v ) = α ( u,v ), α ∈ R; (3) ( u + v , w ) = ( u,w ) + (v,w ), ∀u,v,w ∈ X ; (4) ( u, u) ≥ 0, 当且仅当 u = 0时, , u) = 0. (u 则称( u, v )为X上的u与v的内积. 定义了内积的线性空间 称 的共轭, 为内积空间. (v , u)为( u,v )的共轭,当 K = R时 (v , u) = ( u,v ).
2)
j =1
∑ α ju j = 0 ⇔ ( ∑ α ju j , ∑ α ju j ) = 0
j =1 n j =1
n
n
n
⇔ ( ∑ α j u j , uk ) = 0, k = 1,L, n.
j =1
∴ G非奇异 ⇒ u1 , u2 ,L, un线性无关 (反证法 );反之亦然 .
在内积空间X上可以由内积导出一种范数, 即对u ∈ X , 记 || u ||= (u , u ), Cauchy − Schwarz不等式得出. (1.10) 易证它满足范数定义的正定性和齐次性, 而三角不等式由
数值分析第三章线性方程组解法
数值分析第三章线性方程组解法在数值分析中,线性方程组解法是一个重要的主题。
线性方程组是由一组线性方程组成的方程组,其中未知数的次数只为一次。
线性方程组的解法包括直接解法和迭代解法两种方法。
一、直接解法1.1矩阵消元法矩阵消元法是求解线性方程组的一种常用方法。
这种方法将方程组转化为上三角矩阵,然后通过回代求解得到方程组的解。
1.2LU分解法LU分解法是将系数矩阵A分解为一个下三角矩阵L和一个上三角矩阵U的乘积,然后通过解两个三角方程组求解线性方程组。
这种方法可以减少计算量,提高计算效率。
1.3 Cholesky分解法Cholesky分解法是对称正定矩阵进行分解的一种方法。
它将系数矩阵A分解为一个下三角矩阵L和它的转置的乘积,然后通过解两个三角方程组求解线性方程组。
Cholesky分解法适用于对称正定矩阵的求解,具有较高的精度和稳定性。
二、迭代解法2.1 Jacobi迭代法Jacobi迭代法是一种迭代求解线性方程组的方法。
它通过分解系数矩阵A为一个对角矩阵D和一个余项矩阵R,然后通过迭代更新未知数的值,直至达到一定精度要求为止。
Jacobi迭代法简单易懂,容易实现,但收敛速度较慢。
2.2 Gauss-Seidel迭代法Gauss-Seidel迭代法是一种改进的Jacobi迭代法。
它通过使用新计算出的未知数值代替旧的未知数值,达到加快收敛速度的目的。
Gauss-Seidel迭代法是一种逐步逼近法,每次更新的未知数值都会被用于下一次的计算,因此收敛速度较快。
2.3SOR迭代法SOR迭代法是一种相对于Jacobi和Gauss-Seidel迭代法更加快速的方法。
它引入了一个松弛因子,可以根据迭代的结果动态地调整未知数的值。
SOR迭代法在理论上可以收敛到线性方程组的解,而且收敛速度相对较快。
三、总结线性方程组解法是数值分析中的一个重要内容。
直接解法包括矩阵消元法、LU分解法和Cholesky分解法,可以得到线性方程组的精确解。
数值分析第三章线性方程组迭代法
数值分析第三章线性方程组迭代法线性方程组是数值分析中的重要问题之一,涉及求解线性方程组的迭代法也是该领域的研究重点之一、本文将对线性方程组迭代法进行深入探讨。
线性方程组的一般形式为AX=b,其中A是一个n×n的系数矩阵,x和b是n维向量。
许多实际问题,如电路分析、结构力学、物理模拟等,都可以归结为求解线性方程组的问题。
然而,当n很大时,直接求解线性方程组的方法计算量很大,效率低下。
因此,我们需要寻找一种更高效的方法来求解线性方程组。
线性方程组迭代法是一种基于迭代思想的求解线性方程组的方法。
其基本思想是通过构造一个序列{xn},使得序列中的每一项都逼近解向量x。
通过不断迭代,可以最终得到解向量x的一个近似解。
常用的线性方程组迭代法有雅可比迭代法、高斯-赛德尔迭代法和逐次超松弛迭代法等。
雅可比迭代法是其中的一种较为简单的迭代法。
其基本思想是通过分解系数矩阵A,将线性方程组AX=b转化为x=Tx+c的形式,其中T是一个与A有关的矩阵,c是一个常向量。
然后,通过不断迭代,生成序列xn,并使序列中的每一项都逼近解向量x。
高斯-赛德尔迭代法是雅可比迭代法的改进方法。
其核心思想是利用当前迭代步骤中已经求得的近似解向量的信息。
具体而言,每次迭代时,将前一次迭代得到的近似解向量中已经计算过的分量纳入计算,以加速收敛速度。
相比于雅可比迭代法,高斯-赛德尔迭代法的收敛速度更快。
逐次超松弛迭代法是高斯-赛德尔迭代法的改进方法。
其核心思想在于通过引入一个松弛因子ω,将高斯-赛德尔迭代法中的每次迭代变为x[k+1]=x[k]+ω(d[k+1]-x[k])的形式,其中d[k+1]是每次迭代计算得到的近似解向量的一个更新。
逐次超松弛迭代法可以根据问题的特点调整松弛因子的值,以获得更好的收敛性。
除了以上提到的三种迭代法,还有一些其他的线性方程组迭代法,如SOR迭代法、共轭梯度法等。
这些方法都具有不同的特点和适用范围,可以根据问题的具体情况选择合适的迭代法。
数值分析第三章插值法
( n 1) 又 ( n1) ( t ) Rn ( t ) K ( x )[( t x0 )( t x1 ) ( t xn )]( n1)
19
插值误差举例
f ( 3) ( ) 抛物线插值:R2 ( x ) ( x x0 )( x x1 )( x x2 ) 3!
x0=0.4, x1=0.5, x2=0.6, (0.4, 0.6)
f ( 3) ( ) 2 3 31.25
31.25 R2 (0.54) (0.54 0.4)(0.54 0.5)(0.54 0.6) 3! 0.00175 R1 (0.54) 0.048
试估计线性插值和抛物线插值计算 ln 0.54 的误差
解 线性插值
f ( 2) ( ) R1 ( x ) ( x x0 )( x x1 ) 2
f ( 2) ( ) 2 4
x0=0.5, x1=0.6, (0.5, 0.6)
R1 (0.54) 2(0.54 0.5)(0.54 0.6) 0.0048
ln 0.54 的精确值为:-0.616186···
可见,抛物线插值的精度比线性插值要高 Lagrange插值多项式简单方便,只要取定节点就可写 出基函数,进而得到插值多项式,易于计算机实现。
11
Lagrange插值
lk(x) 的表达式
由构造法可得
( x x0 ) ( x xk 1 )( x xk 1 ) ( x xn ) lk ( x ) ( xk x0 ) ( xk xk 1 )( xk xk 1 ) ( xk xn )
数值分析--第三章--迭代法
数值分析--第三章--迭代法迭代⼀般⽅程:本⽂实例⽅程组:⼀.jacobi迭代法从第i个⽅程组解出xi。
线性⽅程组Ax=b,先给定⼀组x的初始值,如[0,0,0],第⼀次迭代,⽤x2=0,x3=0带⼊第⼀个式⼦得到x1的第⼀次迭代结果,⽤x1=0,x3=0,带⼊第⼆个式⼦得到x2的第⼀次迭代结果,⽤x1=0,x2=0带⼊第三个式⼦得到x3的第⼀次迭代结果。
得到第⼀次的x后,重复第⼀次的运算。
转化成⼀般的形式:(其中L是A的下三⾓部分,D是A的对⾓元素部分,U 是上三⾓部分)得到迭代公式:其中的矩阵B和向量f如何求得呢?其实,矩阵B的计算也很简单,就是每⾏的元素/该⾏上的对⾓元素⼆.Gauss-Seidel迭代法【收敛速度更快】这个可以和jacobi法对⽐进⾏理解,我们以第⼆次迭代为例(这⾥的第⼀次迭代结果都⽤⼀样的,懒得去换)从上表对⽐结果可以看出,Jacobi⽅法的第⼆次迭代的时候,都是从第⼀次迭代结果中,获取输⼊值。
上⼀次迭代结果[2.5,3.0,3.0],将这个结果带⼊上⾯式⼦1,得到x1=2.88,;将[2.5,3.0,3.0]替换成[2.88,3.0,3.0]带⼊第⼆个式⼦的运算,这⾥得到x2=1.95,所以把[2.88,3.0,3.0]替换成[2.88,1.95,3.0]输⼊第三个式⼦计算X3=1.0.这就完成了这⼀次的迭代,得到迭代结果[2.88,1.95,1.0],基于这个结果,开始下⼀次迭代。
特点:jacobi迭代法,需要存储,上⼀次的迭代结果,也要存储这⼀次的迭代结果,所以需要两组存储单元。
⽽Gauss-Seidel迭代法,每⼀次迭代得到的每⼀个式⼦得到的值,替换上⼀次迭代结果中的值即可。
所以只需要⼀组存储单元。
转化成⼀般式:注意:第⼆个式⼦中的是k+1次迭代的第⼀个式⼦的值,不是第k次迭代得值。
计算过程同jacobi迭代法的类似三.逐次超松弛法SOR法上⾯仅仅通过实例说明,Jacobi和Seidel迭代的运算过程。
数值分析(颜庆津)第三章 学习小结
第三章 矩阵特征值与特征向量的计算--------学习小结一、 本章学习体会通过本章的学习,我们学到了四种矩阵特征值和特征向量的计算方法,分别是幂法、反幂法、Jacobi 方法和QR 方法。
四种方法各有其特点和适用范围。
幂法主要用于计算矩阵按模最大的特征值及其相应的特征向量;反幂法主要用于计算矩阵按模最小的特征值及其相应的特征向量;Jacobi 方法用于求实对称矩阵的全部特征值和特征向量的方法;QR 方法则适用于计算一般实矩阵的全部特征值,尤其适用于计算中小型实矩阵的全部特征值。
归结起来,这四种方法亦有其共同点,那就是都是用了迭代的方法来求矩阵的特征值和特征向量。
此外,用MA TLAB 自带的解法求解特征值和特征向量也非常快速,而且不用编辑函数建立m 文件。
其自带函数Eig 功能强大,即便得到结果是虚数也可以算出,并且结果自动正交化。
二、 本章知识梳理本章对于矩阵的特征值和特征向量的算法提出了新的思路,如幂法和反幂法、Jacobi 、QR 方法等。
本章的小结主要从方法的思想,以及一些定理展开。
2.1各种方法的运用范围1、幂法:主要用于计算矩阵按模最大的特征值和其相应的特征向量;2、反幂法:主要计算矩阵按模最小的特征值以及其相应的特征向量;3、Jacobi 方法:用于求实对称矩阵的全部特征值和特征向量的方法;4、QR 方法:适用于计算一般实矩阵的全部特征值,尤其适用于计算中小型实矩阵的全部特征值。
2.2各种方法的基本思想以及迭代公式 1.幂法幂法的基本思想:设n×n 实矩阵A 具有n 个线性无关的特征向量n x x x ,....,,21,其相应的特征值,,...,21n λλλ满足不等式n λλλλ≥≥>....321,其中iix i Ax λ=)...,3,2,1(n i =。
任取一n 维非零向量u 0,从u 0出发,按照如下的递推公式...)2,1(1===k Au u k k因n 维向量组n x x x ,....,,21线性无关,故对向量u 0必存在唯一的不全为0的数组a 1,a 2,...,a n ,使得n n x a x a x a u +++= (22110)])(...)([ (1)212211122211122110221n n n k n kn n k k n k n k k k k k k x a x a x a x a x a x a x A a x A a x A a u A u A Au u λλλλλλλλ+++=+++=+++=====--设a 1≠0,由上式可以看出,当k 充分大时有111x a u kkλ≈得迭代公式:9u A u k k = (1)从实际中来看,为了避免迭代向量u k 的模过大,(当11>λ)或过小(当11<λ),通常对u kj 进行归一化,使其范数等于1。
数值分析第三章小结
第三章矩阵特征值与特征向量的计算--------学习小结一、本章学习体会本章我们学习了矩阵特征值与特征向量的计算方法即幂法、反幂法、Jacobi方法和QR方法。
下边介绍一下四种方法各自的特点和适用范围。
幂法:主要用于计算矩阵按模最大的特征值及其相应的特征向量;反幂法:主要用于计算矩阵按模最小的特征值及其相应的特征向量;Jacobi法:用于求实对称矩阵的全部特征值和特征向量的方法;QR法:则适用于计算一般实矩阵的全部特征值,尤其适用于计算中小型实矩阵的全部特征值。
归结起来,这四种方法有一个共同的特点,即都是用了迭代的方法来求矩阵的特征值和特征向量。
还有利用用MATLAB自带的解法求解特征值和特征向量,其自带函数Eig即得到结果是虚数也可以算出,并且结果自动正交化。
二、本章知识梳理在工程技术中,计算矩阵的特征值和特征向量主要使用数值解法。
本章将阐述幂法、反幂法、Jacobi 方法、和QR 方法,并且只限于讨论实矩阵的情况。
3.1 幂法和反幂法(1)幂法幂法主要用于计算矩阵的按模为最大的特征值和相应的特征向量,其思想是迭代。
设n ⨯n 实矩阵A 具有n 个线性无关的特征向量,,...,,321n x x x x 其相应的特征值n λλλ...21,,满足如下不等式 n λλλλ≥≥≥> (321)其中i i i x Ax λ= )。
(n i ,...2,1=现在要求出1λ和相应的特征向量。
任取一n 维非零向量0u ,从0u 出发,按照如下的递推公式 1-=k k Au u ),,(...21=k 因n 维向量组n x x x ,...,21线性无关,故对于向量0u ,必存在唯一的不全为零的数组n ααα,...,21,使得n n x x x u ααα...22110++=n k n k k k k k k x A x A x A u A u A Au u ααα+++=====--......22110221=⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝⎛+=+++n kn n k kn k n n k k x x x x x x 12122111222111......λλαλλααλλαλαλα 设01≠α。
数值分析第3章
20
定义了内积的线性空间称为内积空间. 定义中(1)的右端 (u称, v为) 的(u共,轭v), 当K为实数域R时 (u, v) .(v, u) 如果 (u, v,) 则 称0 与 正交u ,这v 是向量相互垂 直概念的推广.
b a
f
2
(
x)dx
2
33
若 0 ,1,,n是 C[a, b]中的线性无关函数族,记 span{0 ,1,,n}, 它的格拉姆矩阵为
G G(0 ,1,,n )
(0 ,0 ) (0 ,1) (0 ,n )
(1
,
0
)
(1 , 1 )
(1
,
n
)
(n ,0 )
(n ,1 )
(
n
,
n
)
(1.17)
Hn span{1, x,, xn},
且 (a0 , a1,, an ) 是 p(x) 的坐标向量,H n 是 n 1维的.
8
对连续函数 f (x) C[a,b],它不能用有限个线性无关的 函数表示,故 C[a,b]是无限维的,但它的任一元素 f (x) 均可用有限维的 p(x) Hn逼近,使误差
与数的乘法构成实数域上的线性空间, 记作 R n,称为 n维
向量空间.
4
对次数不超过 n( n为正整数)的实系数多项式全体,
按通常多项式与多项式加法及数与多项式乘法也构成数域
R上一个线性空间,用
H
表示,称为多项式空间.
n
所有定义在 [a,b] 上的连续函数集合,按函数加法和
数值分析第三章小结
矩阵的特征值和特征向量的计算线性代数中对于x Ax λ=,解该方程的特征值λ和特征向量x 的方法主要是使用数值解法,本章学习另外的方法用MATLAB 来编程解某个实矩阵的特征值和特征向量. 一、幂法和反幂法 1、乘幂法幂法主要用于计算矩阵的按模为最大的特征值和对应的特征向量。
(1)思想为: ,2,1,1==-k u A u k kn n X X X u ααα+++= 22110])([2111101∑=-+===ni i ki i kkk k X X u A u A u λλααλ当k 取得足够大时,特征值向量得计算公式为:111X u kk αλ≈特征值为:111111111)()()()(λαλαλ=≈++ik i k ik i k X X u u迭代格式为之一⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=====∈-------- ,2,1111111110k u y y A u uy u u R u k Tk k k k k k k k Tk k n βηη任取初始向量迭代格式之二⎪⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎪⎨⎧=======------≤≤- ,2,1)sgn(),,,(max ),,,(任取初始向量)()1()()(2)(11)1(11)1(1)1()0()0(2)0(10k h h h h h y A u h u y h h h h h u k r k r k Tk n k k k k k r k k k j nj k r Tn β两种迭代格式相比较,格式一编程容易,迭代一次所需时间也短,迭代格式二迭代时间长,但它在计算过程中舍入误差的影响较格式一小。
幂法的缺点是如果矩阵A 的特征根有重根时不能用。
2、反幂法目的同乘幂法,用于计算矩阵的按模为最大的特征值和对应的特征向量。
反幂法的迭代格式为⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=====∈-------- ,2,1任取初始向量111111110k u y y A u uy u u R u k T k k k k k k k k Tk k n βηη3、带原点位移的反幂法迭代格式为 ,2,1)max(11=⎪⎪⎩⎪⎪⎨⎧===--k m y u y m u A y k kkk k k k三、Jacobi 方法和QR 方法Jacobi 方法主要用于求实对称矩阵的全部特征值和特征向量的一种方法,所以个人觉得雅克比法更为现实更为有用。
数值分析第三章Ch3
则 ρ(x) 称为 [a, b] 上的一个权函数。
数值逼近
数值分析
. . . .... .... .... . .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..
15/88
. .. . .. .. ..
则 ρ(x) 称为 [a, b] 上的一个权函数。
10/88
. .. . .. .. ..
设 S 是一个内积空间,u1, . . . , un ∈ S,
矩阵
G
=
((uu11,,...
u1) u2)
(u2, u1)
(u2, u2) ...
··· ···
(un,
(un, ...
uu12))
(u1, un) (u2, un) · · · (un, un) 称为格拉姆 (Gram) 矩阵。
间
[−π, π]族。 . . . .... .... .... .
.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..
数值逼近 数值分析
19/88
. .. . .. .. ..
称多项式序列{φn}∞ n=0 是带权 ρ(x) 正 交的,如果每个 φn 是首项系数 an ̸= 0 的 n 次多项式,且 {φn}∞ n=0 是带权 ρ(x) 的正交函数族。
. .. . .. .. ..
. .四、最佳逼近
函数逼近主要讨论给定 f ∈ C[a, b],求最
佳逼近多项式。若 P ∗(x) ∈ Hn,使误差 ∥f (x) − P ∗(x)∥ = min ∥f (x) − P (x)∥,则
P ∈Hn
称 P ∗(x) 是 f (x) 在 [a, b] 上的n 次最佳逼
数值分析(第三章)实验报告
L1 ( x)
L2 ( x)
( x 0)( x 0.6) (100*x*(x - 3/5))/27 (0.9 0)(0.9 0.6)
30)*(x 9/10))/27 + -
P2 ( x) L0 ( x) cos 0 L1 ( x) cos 0.6 L2 ( x) cos 0.9 =((50*x
可以预测 1930,1965,2010 年的人口分别是 169649,1.9177e+005,171351
EXERCISE SET 3.2 4、
P131
a) 根据 Algorithm 3.2,利用课本作者网站上的关于本书的 MATLAB 程序 ALG032.M 运行该程序,在命令行窗口出现如下: Warning: Could not find an exact (case-sensitive) match for 'ALG032'. E:\ 个人 \ 工作 \ 高教数值分析 \ 第八版英文程序 \Matlab-Programs\matlab\m1\ALG032.M is a case-insensitive match and willbe used instead. You can improve the performance of your code by using exact name matches and we therefore recommend that you update your usage accordingly. Alternatively, you can disable this warning using warning('off','MATLAB:dispatcher:InexactCaseMatch'). This warning will become an error in future releases. Newtons form of the interpolation polynomial Choice of input method: 1. Input entry by entry from keyboard 2. Input data from a text file 3. Generate data using a function F Choose 1, 2, or 3 please 1 Input n 4 Input X(0) and F(X(0)) on separate lines 0.0 -6.00000 Input X(1) and F(X(1)) on separate lines 0.1 -5.89483 Input X(2) and F(X(2)) on separate lines 0.3 -5.65014 Input X(3) and F(X(3)) on separate lines 0.6 -5.17788 Input X(4) and F(X(4)) on separate lines 1.0 -4.28172 Select output destination 1. Screen 2. Text file Enter 1 or 2 1 NEWTONS INTERPOLATION POLYNOMIAL Input data follows: X(0) = 0.00000000 F(X(0)) =
《数值分析》课件-第三章
数值分析
A为正交矩阵的充要条件是下列条件之一成立:
1 A1 AT ;
2 AAT AT A E;
3 A的列向量是两两正交的单位向量; 4 A的行向量是两两正交的单位向量.
正交矩阵的性质
(1) A是正交阵,则A是非奇异的,det( A) 1;
(2) A是正交阵,则AT A1也是正交阵,且同阶 正交阵的乘积仍是正交阵;
( , ) 由Schwarz不 等 式 可 以 证 明 内 积 范数 公 理 中 的 三 角 不 等 式.
证明:三角不等式
证:在内积空间V中, , V , 有 2 ( , ) ( , ) 2( , ) ( , ) 2 2 2 ( )2 所以
数值分析
数值分析
定理2-6 : (Cauchy Schwarz不等式)
设 , 是内积空间V中任意两个向量,则有 ( , )2 ( , )( , )
等号只有当且仅当和 是线性相关时才成立.
证明 : 任取实数k, 考虑内积
( k , k ) ( , ) 2k( , ) k 2( , ) 0 利用一元二次方程根的判别式, 有4( , )2 4( , )( , ) 0 所以有( , )2 ( , )( , ) 当 k (k R,非 零),显 然 定 理 中 等 号 成 立;反 之, 如 果 等 号 成 立,则 , 必 线 性 相 关.因 为 若 , 线 性 无 关,则k R, 非 零, 都 有 k 0.从 而( k , k ) 0 所 以 等 号 不 成 立,矛 盾.
数值分析
第三章
数值分析
数值分析
第一节 内积空间与内积空间中的正交系
一、内积和内积空间的基本概念
定义2-12:设V是实数域R上的线性空间,如果α,β V
数值分析第三章线性方程组的迭代法课件
§ 3.3.2 Gauss—Seidel 迭代法的矩阵表示
将A分裂成A =D+L+U,则Ax b 等价于
(D+L+U )x = b
于是,则高斯—塞德尔迭代过程
Dx(k1) Lx(k1) Ux(k) b
因为 D 0 ,所以 D L D 0
故
(D L)x(k1) Ux(k) b
x(k1) (D L)1Ux(k) (D L)1b
e(k) x(k) x* Gx(k1) d (Gx* d) G(x(k1) x* ) Ge(k1)
于是 e(k) Ge(k1) G 2e(k2) Gk e(0)
由于 e (0)可以是任意向量,故 e(k) 收敛于0当且仅
故 (D L)x(k1) (1)D U x(k) b
显然对任何一个ω值,(D+ωL)非奇异, (因为假设 aii 0,i 1,2,, n )于是超松弛迭代公式为
x(k1) (D L)1 (1)D U x(k) (D L)1b
令 L (D L)1 (1)D U
f (D L)1b
则超松弛迭代 公式可写成
称为雅可比迭代公式, B称为雅可比迭代矩阵
雅可比迭代矩阵表示法,主要是用来讨论其收敛 性,实际计算中,要用雅可比迭代法公式的分量 形式。即
x1(k 1)
1 a11
(a12 x2(k )
a13 x3(k )
a1n xn(k )
b1 )
x2(k 1)
1 a 22
(a21 x1(k )
a23 x3(k )
§ 3.4.2超松弛迭代法的矩阵表示 设线性方程组 Ax=b 的系数矩阵A非奇异,且主对角
元素 aii 0(i 1,2,, n) , 则将A分裂成
第三章 Biot固结方程的数值分析
3.10
Q u t H hn1 F2 n F2 t 1 H t hn
当把权函数w取为常数时,
3.11
,这样的时间差分格式即相当于 w 时, 中心差分格式,常称为Crank-Nicholson格式。当 2 / 3 ,是目前常用的一种差分格式。经验表明,这时即使两 个时段之间 t 发生突变,误差也不明显增大。 1相当于全隐
式中
w d / t wdr
0
r
r
W为积分时所用的权函数。 如果在 时刻体系满足式(3.2),式(3.6)右端括号中的值 为零,并把劲度矩阵[K]换成切线矩阵 Kt ,同时去掉初切应变 的等价力F0 ,即令
可得
T
Kt u K u F0 Kt u Q h F1
3.1
本章限于讨论荷载缓慢变化条件下惯性项可以忽略的情况,这 时离散化后的Biot固结方程可以简化为 3.2 K u Qh F1 F0
H h F S h Q u 2
T
3.3
上式为时间的一次微分方程,常用两点递进格式求解。设 t n 和 tn 1 为时域上两点,相应的变量为un 和un1 及hn 和hn1 ,则 h / t ,并采用下 u /t 和h 在 t tn1 tn 时刻内可令 u 列线性内插公式,
3.8
Q u S t H h F2 n F2 t H t hn 3.9 对饱和土 S 0 的情况,许多人习惯用hn1 作为未知量,即把
式(3.8)和式(3.9)改写为
T
Kt u Qhn1 F1 Qhn
''' 3 2 y ' ih yi 1 yi 1 2 yi h 3!
03 材料科学研究中常用的数值分析方法
解决这类问题通常有两种途径:(1)对方 程和边界条件进行简化,从而得到问题在简 化条件下的解答;(2)采用数值解法。 第一种方法只在少数情况下有效,因为过多 的简化会引起较大的误差,甚至得到错误的 结论。 目前,常用的数值解法大致可以分为两类: 有限差分法和有限元法。
数值模拟通常由前处理、数值计算、后处理三 部分组成 前处理 实体造型、物性赋值、定义单元类型、网格 划分 数值计算 施加载荷、设定时间步、确定计算控制条件、 求解计算 后处理 显示和分析计算结果、分析计算误差
1.差分方程的建立
合理选择网格布局及步长 将离散后各相邻离散点之间的距离,或者离散 化单元的长度称为步长。
y
(i,j+1)
dy
(i-1,j) (i,j)
(i+1,j)
(i,j-1)
dx x
将微分方程转化为差分方程
向前差分
T T (i 1, j ) T (i, j ) x x
2
1 1 T (i, j ) T (i, j ) T (i, j 1) 2T (i, j ) T (i, j 1) T 2 2 2 y y y y 2
2
差分格式的物理意义
y
dT dx T(x+dx)-T(x) dx
2u 2u 2 0,0 x 0.5,0 y 0.5 2 x y u (0, y ) u ( x,0) 0 u ( x,0.5) 200x u (0.5, y ) 200y
3.3 有限单元法
有限元法(FEMA)也称为有限单元法或有限元素 法,基本思想是将求解区域离散为一组有限个且按 一定方式相互连接在一起的单元的组合体。它是随 着电子计算机的发展而迅速发展起来的一种现代计 算方法。 把物理结构分割成不同大小、不同类型的区域,这 些区域称为单元。 根据不同分析科学,推导出每一个单元的作用力方 程,组集成整个结构的系统方程,最后求解该系统 方程,就是有限元法。
第三章习题解答 _数值分析
第三章习题解答1、 利用Gauss 变换阵,求矩阵A 2113113112⎡⎤⎢⎥⎢⎥=⎢⎥-⎢⎥-⎣⎦的LU 分解(要求写出分解过程)。
解:因为,k k i k i k k k m a a ()(),,,/=可得,L 111/2111⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥⎣⎦,A 1215/2113112⎡⎤⎢⎥⎢⎥=⎢⎥-⎢⎥-⎣⎦L 2112/511⎡⎤⎢⎥⎢⎥=⎢⎥-⎢⎥⎣⎦,A 2215/2113/5112⎡⎤⎢⎥⎢⎥=⎢⎥-⎢⎥-⎣⎦ L 31115/131⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥-⎣⎦,U 215/2113/5121/13⎡⎤⎢⎥⎢⎥=⎢⎥-⎢⎥-⎣⎦ L L L L 12311/212/515/131⎡⎤⎢⎥-⎢⎥==⎢⎥-⎢⎥-⎣⎦9、112221,m mm B C A B C A A A C A B -⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦已知为块三对角阵非奇异,=, i m m mm mm B i m A LU B C L I U A B C R L I U R C U A B R L I 11112222223111,2,,),--=⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦其中均为方程(设有分块分解式i i i i i i i i i R A i m L B U L C L B A U i m U L C i m 111111111(2,3,...,)23(2,3,...,)4(2,3,...,).---=====-===试证明:(),;(),;();() 证:1122111111111111,,,,,,(2,3,...,)(2,3,...,)(2,3,...,)i i i i i i i i i i i i i i i i i iB L A RC L U B R U L A R C L U L B U L C R A i m L B A U i m U L C i m ----====+==∴=====-===11、已知A 1234⎡⎤=⎢⎥⎣⎦,试计算cond A cond A cond A 12(),(),()∞。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
若集合V是线性空间, 若集合 是线性空间,ϕ 0 ( x), ϕ1 ( x),L , ϕ n ( x) 是线性空间 线性无关, 线性无关,则 ϕ ( x) 可以表示为
ϕ ( x) = c0ϕ 0 ( x) + c1ϕ1 ( x) + L + cnϕ n ( x)
若为线性空间V的一组基, 若为线性空间 的一组基,则 的一组基 是一个n+1维线性空间 维线性空间 是一个
1 2
1 2
证明:对内积空间上的任意元素 、 和任意 证明:对内积空间上的任意元素f、g和任意 实数t, 实数 ,有
0 ≤ ( f + tg , f + tg ) = ( f , f ) + 2 t ( f , g ) + t 2 ( g , g )
固定f和 ,右端是关于t的一元二次多项式 的一元二次多项式, 固定 和g,右端是关于 的一元二次多项式 且该多 项式函数值非负, 项式函数值非负,利用二次多项式的判别式有
1≤ i ≤ n
= max{ vi }
n
赋范线性空间 ( R , • ∞ )
无穷范数) 连续函数空间 (无穷范数 无穷范数
定义于区间[a,b]上连续函数的集合 上连续函数的集合C[a,b]是一线性空 定义于区间 上连续函数的集合 是一线性空 间。定义
f
∞
= max f ( x) , ∀f ∈ C[a, b]
定义于区间[a,b]上连续函数的集合 上连续函数的集合C[a,b]是一线性空间。 是一线性空间。 定义于区间 上连续函数的集合 是一线性空间 1 定义 b f 2 = [ ∫ f 2 ( x)dx] 2 ≥ 0∀f ∈ C[a, b]
a
• 2是线性空间 是线性空间C[a,b]上的一种范数。C[a,b]关 上的一种范数。 上的一种范数 关
f ( x ) g ( x ) ρ ( x ) dx
Remark: :
区间端点可以是无穷大,此时为广义积分。 区间端点可以是无穷大,此时为广义积分。 常简记 (•, •) ρ 为 (•, •) 没有确切指出权函数时,约定ρ(x)=1。 没有确切指出权函数时,约定 。 在理论证明和公式推导过程中, 在理论证明和公式推导过程中 如没有明确 权函数具体形式, 则表示对任意权函数均成立。 权函数具体形式 则表示对任意权函数均成立。
第三章 函数逼近
一、问题的提出
逼 f ( x) 的话, 近 的话 运用插值函数有时就要失败。 pn ( x) ,运用插值函数有时就要失败。 级数的前n+1项部分和,其 项部分和, 设 为f(x) Taylor级数的前 级数的前 项部分和 截断误差为 f ( n +1) (ξ ) Rn ( x) = f ( x) − pn ( x) = ( x − x0 ) n +1ξ ∈ ( x0, x)
rr r r n ( R , • ∞ )上d(u, v)= u − v
2
i =1
∞
= max ui − vi
1≤i ≤ n
∞
(C[a, b], • ∞ )上d(f,g)= f − g
(C[a, b], • 2 )上d(f,g)= f − g
2
= max f ( x) − g ( x)
a ≤ x ≤b
( f , g) =
∫
b
a
f ( x ) g ( x ) dx
满足内积的三个条件 线性空间 C[ a, b] 关于上式所规定的内积是一 内积空间。 内积空间。
2、性质 、
满足Cauchy不等式 ①内积空间上任意两元素f和g满足 内积空间上任意两元素 和 满足 不等式
f , g ) ( f , f ) ( g , g ) ( ≤
(n + 1)!
如果要求解在[a,b]区间的每一点都“很好地” 区间的每一点都“很好地” 如果要求解在 区间的每一点都
的点且误差也较小, 为了照顾到远离 的点且误差也较小,常将阶 取得很大。 数n取得很大。这样做费事、多占存贮单元。因此往 取得很大 这样做费事、多占存贮单元。 往要求在给定精度下 求形式简单的计算公式, 在给定精度下, 往要求在给定精度下,求形式简单的计算公式,使 f ( x) 这就是函数逼近要解决的问题。 其均匀地逼近 。这就是函数逼近要解决的问题。
1 2 1 2 1 2
1 2 1 2
1 2 1 2
1 2
2、距离 、
对于赋范线性空间 (V , • ) 上的任意两个元 素f 和g,它们之间的距离为 d ( f , g ) = f − g1 , n rr r r ( R n , • 2 )上d(u, v)= u − v = [∑ (ui − vi ) 2 ] 2
∆ = 4( f , g ) − 4( g , g )( f , f ) ≤ 0
2
≤ 得到Cauchy不等式: f , g ) ( f , 得到 不等式: ( 不等式
f ) (g, g)
1 2
1 2
内积空间上的任意两元素f和 满足三角不等式 ② 内积空间上的任意两元素 和g满足三角不等式 (Schwarz不等式 : 不等式): 不等式
四、函数逼近问题
为被逼近函数。 设 f ( x ) ∈ C [ a , b ] 为被逼近函数。 Φ为赋范线性空间 (C [ a , b ], • ) 的一个子集合 为赋范线性空间 范数 • 可以是 • ∞ 或 • 2 等。 问题: 称问题:求 ϕ * ( x) ∈ Φ 使得
x0
二、函数逼近问题
已知复杂函数 f ( x) ,或仅知道函数 f ( x) 在某些 采样点处的函数值,在某函数集合V中 采样点处的函数值,在某函数集合 中,寻找 f ( x) 最好” ϕ 的“最好”近似* ( x ) ∈ V 函数逼近问题: 函数逼近问题:对集合中给定的函数 f ( x) , 要求在另一类较简单的便于计算的函数集合V 要求在另一类较简单的便于计算的函数集合V中, * ϕ * ( x) 与 f ( x)之差在某 求函数 ϕ ( x ) ∈ V ⊂ V ,使得 种度量意义下最小。 种度量意义下最小。
( f , f ) ≥ 0 当且仅当 f = 0 时有 ( f , f ) = 0
是线性空间V上的一种内积。 上的一种内积 则称实值函数 (•, •) 是线性空间 上的一种内积。 并称线性空间V关于实值函数 内积空间。 并称线性空间 关于实值函数 (•, •) 是内积空间。 对于线性空间 C[ a, b] , 如下定义的实值函数
b 1 2 a
= ( ∫ ( f ( x) − g ( x)) 2 dx)
Remark
d(f,0)= f − 0 = f
二、内积空间
1、定义 、
为一线性空间, 设V为一线性空间,若定义实值函数 ( •, • ) ,对 为一线性空间 任意 f , g , h ∈ V 满足 对称性) ①(对称性 ( f , g ) = ( g , f ) 对称性 线性性) ②(线性性 (r1 f + r2 h, g ) = r1 ( f , g ) + r2 (h, g ), ∀r1 , r2 ∈ R 线性性 非负性) ③(非负性 非负性
∫
b a
ρ ( x ) dx > 0
Байду номын сангаас
∫
b a
x k ρ ( x ) dx ( k = 0,1, 2, L ) 存在
为区间[a,b]的一个权函数。 的一个权函数 则称 ρ ( x ) 为区间 的一个权函数。
2、带权的内积 C[a,b]带权的内积: 带权的内积: 带权的内积
( f , g )ρ =
∫
b a
( f + g, f + g) ≤ ( f , f ) (g, g)
证明: 利用Cauchy不等式有 证明 利用 不等式有
无法显示图像。计算机可能没有足够的内存以打开该图像,也可能是该图像已损坏。请重新启动计算机,然后重新打开该文件。如果仍然显示红色 “x”,则可能需要删除该图像,然后重新将其插入。
a ≤ x ≤b
∞
=0⇔ f =0
∞
rf
∞
= max rf ( x) = r max f ( x) = r f
a ≤ x ≤b
f +g
a ≤ x ≤b
∞
= max f ( x) + g ( x)
a ≤ x ≤b
≤ max f ( x) + max g ( x) = f
a ≤ x ≤b
∞
+ g
∞
范数) 连续函数空间 (Euclid范数 范数
V = span{ϕ 0 , ϕ1 ,L , ϕ n }
背景:在某一函数集合中找最好的近似。 背景:在某一函数集合中找最好的近似。 最好的近似 赋范空间、内积空间、 赋范空间、内积空间、正交多项式 最佳平方逼近 曲线最小二乘拟合 最佳一致逼近(工科研究生不要求) 最佳一致逼近(工科研究生不要求)
§1 预备知识与函数逼近问题
1 2
1 2
1 2
≤ ( f , f ) + 2( f , f ) ( g , g ) + ( g , g )
= [( f , f ) + ( g , g ) ]
1 2 1 2 2
1 2
1 2
两边开平方, 三角不等式得证。 两边开平方 三角不等式得证。
三、权函数
1、定义
定义在[a,b]上的实值函数 ρ ( x ) ,如果满足 上的实值函数 定义在 ① ρ ( x) ≥ 0, ∀x ∈ [a, b] ② ③
一、赋范线性空间
1、定义
为定义于线性空间V上的实值函数 并满足: 上的实值函数, 设 • 为定义于线性空间 上的实值函数,并满足: 非负性) 当且仅当g=0时有 g = 0 ①(非负性 g ≥ 0 , ∀ g ∈ V 当且仅当 非负性 时有 ②(齐次性 rg = r g , ∀r ∈ R, g ∈ V 齐次性) 齐次性 三角不等式) ③(三角不等式 f + g ≤ f + g , ∀f , g ∈ V 则称 • 是 三角不等式 线性空间V上的范数。并称线性空间V为赋范线性空间, 上的范数 线性空间 上的范数。并称线性空间 为赋范线性空间, 记为 (V , • ) Remark: 子空间 V1 ⊂ V , V上的范数 • 也是 V1 上的范 上的范数 数。