(整理)数值分析计算方法超级总结
大学数学易考知识点数值分析的基本方法和应用

大学数学易考知识点数值分析的基本方法和应用大学数学易考知识点:数值分析的基本方法和应用一、引言数值分析是现代数学在科学计算和工程实践中的应用研究领域,是研究数值计算方法和数值算法的理论与实践的学科。
在大学数学课程中,数值分析是一个重要的知识点,它涉及到数值计算的基本方法和应用。
本文将介绍数值分析的基本方法和应用,以帮助学生更好地理解和掌握这一易考的知识点。
二、数值分析的基本方法1. 插值和逼近插值与逼近方法是数值分析中常用的方法之一,它们用于通过已知数据点构造一个近似函数,以在给定范围内估计未知数据点的值。
常见的插值与逼近方法包括拉格朗日插值、牛顿插值、最小二乘逼近等。
2. 数值微积分数值微积分方法用于对函数进行数值积分和数值微分。
在实际计算中,往往难以通过解析方法求得函数的积分或导数,这时可以利用数值积分和数值微分方法来近似计算。
其中常见的数值积分方法包括梯形法则、辛普森法则等,数值微分方法包括中心差商法、向前差商法、向后差商法等。
3. 常微分方程的数值解法常微分方程数值解法用于求解无法通过解析方法得到解的常微分方程。
常见的常微分方程数值解法有欧拉法、改进欧拉法、龙格-库塔法等,它们根据不同的精度和稳定性要求,选择不同的数值解法来计算常微分方程的近似解。
4. 线性方程组的数值解法线性方程组数值解法是解决线性方程组问题的常见方法。
当线性方程组的规模较大时,无法通过直接求解的方法得到解,此时可以利用数值解法来近似求解。
常见的线性方程组数值解法包括高斯消元法、LU分解法、迭代法等。
三、数值分析的应用1. 插值与逼近的应用插值与逼近方法在科学计算和工程实践中有广泛的应用。
例如,在地理信息系统中,插值方法可以用于根据已知地理数据点生成等高线图;在图像处理中,逼近方法可以用于图像的平滑处理和边缘检测。
2. 数值积分的应用数值积分方法在物理学、经济学等领域的科学研究中有重要的应用。
例如,在物理学中,数值积分方法可以用于计算物体的质心、面积、弧长等物理量;在经济学中,数值积分方法可以用于计算经济指标、积分收益等。
数值分析各章重点公式整理

数值分析各章重点公式整理数值分析是计算数学的一个分支,主要涉及计算和分析数值近似解的方法。
本文将从数值分析的基本概念、插值与逼近、数值微分和数值积分、非线性方程数值解、线性方程组直接解法、线性方程组迭代解法和特征值问题等几个方面,对数值分析中的重点内容进行整理。
一、数值分析的基本概念数值分析是用数值方法解决实际问题的方法和技术。
其主要研究目标是通过一定的计算机运算来获取数学问题的近似解。
数值分析涉及到误差分析、收敛性分析、稳定性分析等概念,对于数值方法的正确性和可行性提供了理论依据。
二、插值与逼近插值是通过已知数据点构造一个函数,使得这个函数通过已知数据点。
常用的插值方法有拉格朗日插值和牛顿插值。
逼近是选择一种较为简单的函数来近似表示给定的复杂函数。
常用的逼近方法有最小二乘法和切比雪夫逼近。
三、数值微分和数值积分数值微分主要研究如何通过函数值的有限差分来估计导数值。
常用的数值微分方法有前向差分、后向差分和中心差分。
数值积分主要研究如何通过数值方法求出函数在一定区间上的定积分值。
常用的数值积分方法有梯形法则和 Simpson 法则。
四、非线性方程数值解非线性方程通常难以用解析方法求解,而数值方法则可以通过迭代来逼近方程的根。
常用的数值解法有二分法、牛顿法和割线法。
同时,对于多维非线性方程,也可以使用牛顿法的变形,牛顿下山法。
五、线性方程组直接解法线性方程组是数值分析中的一个重要问题。
直接解法主要有高斯消元法、LU 分解法和 Cholesky 分解法。
高斯消元法通过矩阵的初等行变换将线性方程组化为上三角方程组来求解。
LU 分解法将系数矩阵分解为一个下三角矩阵和一个上三角矩阵的乘积,然后通过回代求解。
Cholesky 分解法则适用于对称正定矩阵的解法。
六、线性方程组迭代解法线性方程组的迭代解法通过选取适当的初始解,通过迭代来逼近精确解。
常用的迭代解法有Jacobi迭代法、Gauss-Seidel迭代法和超松弛迭代法。
数值分析(计算方法)总结

第一章 绪论误差来源:模型误差、观测误差、截断误差(方法误差)、舍入误差ε(x )=|x −x ∗|是x ∗的绝对误差,e =x ∗−x 是x ∗的误差,ε(x )=|x −x ∗|≤ε,ε为x ∗的绝对误差限(或误差限) e r =ex =x ∗−x x为x ∗ 的相对误差,当|e r |较小时,令 e r =ex ∗=x ∗−x x ∗相对误差绝对值得上限称为相对误差限记为:εr 即:|e r |=|x ∗−x||x ∗|≤ε|x ∗|=εr绝对误差有量纲,而相对误差无量纲若近似值x ∗的绝对误差限为某一位上的半个单位,且该位直到x ∗的第一位非零数字共有n 位,则称近似值 x ∗有n 位有效数字,或说 x ∗精确到该位。
例:设x=π=3.1415926…那么x ∗=3,ε1(x )=0.1415926…≤0.5×100,则x ∗有效数字为1位,即个位上的3,或说 x ∗精确到个位。
科学计数法:记x ∗=±0.a 1a 2⋯a n ×10m (其中a 1≠0),若|x −x ∗|≤0.5×10m−n ,则x ∗有n 位有效数字,精确到10m−n 。
由有效数字求相对误差限:设近似值x ∗=±0.a 1a 2⋯a n ×10m (a 1≠0)有n 位有效数字,则其相对误差限为12a 1×101−n由相对误差限求有效数字:设近似值x ∗=±0.a 1a 2⋯a n ×10m (a 1≠0)的相对误差限为为12(a 1+1)×101−n 则它有n 位有效数字令x ∗、y ∗是x 、y 的近似值,且|x ∗−x|≤η(x )、|y ∗−y|≤η(y)1. x+y 近似值为x ∗+y ∗,且η(x +y )=η(x )+η(y )和的误差(限)等于误差(限)的和2. x-y 近似值为x ∗−y ∗,且η(x +y )=η(x )+η(y )3. xy 近似值为x ∗y ∗,η(xy )≈|x ∗|∗η(y )+|y ∗|∗η(x)4. η(xy )≈|x ∗|∗η(y )+|y ∗|∗η(x)|y ∗|21.避免两相近数相减2.避免用绝对值很小的数作除数 3.避免大数吃小数 4.尽量减少计算工作量 第二章 非线性方程求根1.逐步搜索法设f (a ) <0, f (b )> 0,有根区间为 (a , b ),从x 0=a 出发, 按某个预定步长(例如h =(b -a )/N )一步一步向右跨,每跨一步进行一次根的搜索,即判别f (x k )=f (a +kh )的符号,若f (x k )>0(而f (x k -1)<0),则有根区间缩小为[x k -1,x k ] (若f (x k )=0,x k 即为所求根), 然后从x k -1出发,把搜索步长再缩小,重复上面步骤,直到满足精度:|x k -x k -1|< 为止,此时取x *≈(x k +x k -1)/2作为近似根。
数值分析公式

2
xk
1 2
h k
0,1,, n
每个区间上用Simpson求积公式
b f x dx n a k 1
xk f x dx
xk 1
h 6
n k 1
f
xk1
4
f
x
k
1
2
f xk E~n f
Sn f
h 6
f
n
a 4
k 1
f
x k
1 2
2
n1 k 1
f
xk
f
b
E~n
f
b a h4 2880
b1 c1 a2 b2 c2
x1 d1
x2
d2
ai bi ci
xi
di
简记 Ax d.
an1 bn1 cn1 xn1 dn1
an bn xn dn
此系数矩阵的非零元素集中分布在主对角线及其相邻两次对角线
上,称为三对角矩阵。方程组称为三对角方程组。
sx a00 x a11x amm x 很好地拟合
如 数据大体规律为 y aebx , y x 等, 可通过一定 ax b
的变换的到线性模型
y aebx ln y ln a bx, y x 1 a 1 b
ax b
yx
正交函数
定义 : 设 f x, gx Ca,b, x 是 a,b上的权函数, 则
显式Euler公式:yn1 yn h f xn , yn
改 进E u ler公 式
yn1
yn
hf
2
xn ,
yn
f
xn 1 ,
yn1
梯形E
u
ler公式
现代数值计算方法公式总结

现代数值计算方法公式总结现代数值计算方法公式一、插值法1.拉格朗日(Lagrange)插值法a)两点一次:b)三点二次:2.牛顿(Newton)插值a)n次牛顿法多项式:其中一阶差二阶差商三阶差商四阶差商商b)向前差分:下减上c)向后差分:上减下3.三次埃米尔特(Hermite)插值二、拟合曲线(最小二乘)三、数值积分1.牛顿-柯特思(Newton-Cotes)公式梯形求积公式(2节点)复化梯形求积公式辛普生求积公式(3节点)复化辛普生求积公式2.高斯(Gauss)公式高斯-勒让德求积公式1.先用勒让德公式求解x i2.利用“高斯积分公式具有2n+1次代数精度”将x i带入求A i3.将xi、Ai带入公式求取积分、并计算误差。
普通积分化标准形式:积分区间[a,b]变换3.代数精度若求积公式对f(x)=1,x,x2,…x m时精确成立,而对f(x)=x m+1时不成立,则称此求积公式具有m次代数精确度四、解线性代数方程组的直接方法三角形分解法求解,先将A分解为,则原式变为,那么问题就变为了求解五、解线性代数方程的迭代法1.范数向量范数定义:设其中R为实数域、C为复数域,若某实值函数满足条件1)非负性,||x||=0当且仅当x=0成立2)其次行3)三角不等式称为域上的一个向量范数常见范数:矩阵范数定义:设其中R为实数域、C为复数域,若某实值函数满足条件1)非负性,||A||=0当且仅当A=0成立2)其次行3)三角不等式4)乘积性质称为域上的一个矩阵范数常见范数:行范数列范数为的最大按模特征值2.谱半径3.雅可比迭代向量:用第i个方程解出xi的方程,分量通式如下:矩阵:对于Ax=b,先将A拆分成对角线矩阵D减去下三角矩阵L,再减去上三角矩阵U。
其中4.高斯-塞德尔迭代向量:用第i个方程解出xi的方程,并将上式得到的带入下边的公式,分量通式如下:矩阵:对于Ax=b,先将A拆分成对角线矩阵D减去下三角矩阵L,再减去上三角矩阵U。
数值计算方法期末总结

数值计算方法期末总结导言数值计算是近年来发展迅速的一门学科,它研究如何利用数字近似计算数学方程和问题的解。
在科学计算、工程分析、金融建模等领域都有广泛应用。
本文将对数值计算方法进行总结,包括数值逼近、插值与外推、数值微积分、线性方程组解法、非线性方程解法、数值积分与数值微分以及随机数生成与蒙特卡洛方法。
通过总结这些方法的基本原理、优缺点和应用领域,可以帮助读者更好地理解和运用数值计算方法。
一、数值逼近数值逼近是指通过有限次数的计算,利用某一数列逐步逼近函数的值。
数值逼近可以分为插值和外推。
插值是在给定的有限个数据点之间找到一个函数,使得函数经过这些数据点。
而外推是利用已知数据点的决策逐渐增加,以获得更精确的近似值。
在实际应用中,数值逼近被广泛应用于数据处理和数据分析中,常用于构造曲线拟合、图像处理和信号处理中。
数值逼近的方法有拉格朗日插值、牛顿插值和埃尔米特插值等。
二、插值与外推插值与外推是数值计算中用于估计未知函数值的重要工具。
插值是在给定数据点之间构造一个模型函数,使得函数经过这些数据点。
外推是利用一些已知数据点的决策逐渐逼近未知函数的方向。
常用的插值与外推方法有多项式外推、样条插值、最小二乘法、有限差分法等。
它们可以用于函数逼近、数据拟合和数值求解等问题。
三、数值微积分数值微积分是一种利用数值方法来近似计算积分和求解微分方程的方法。
数值微积分广泛应用于工程计算、金融建模和科学研究等领域,是计算机辅助设计和分析的关键技术之一。
在数值微积分中,常用的方法有数值积分和数值微分。
数值积分主要用于求解曲线下面积和计算函数的平均值等问题,常用方法有复合梯形公式、复合辛普森公式、复合高斯公式等。
而数值微分主要用于近似计算函数的导数,常用方法有有限差分法、龙贝格公式和微分方程的数值解法等。
四、线性方程组解法线性方程组是科学计算中的重要问题之一,其求解方法的好坏直接影响到计算结果的精度和稳定性。
线性方程组的求解方法有直接法和迭代法两种。
数值分析实验报告总结

一、实验背景数值分析是研究数值计算方法及其理论的学科,是计算机科学、数学、物理学等领域的重要基础。
为了提高自身对数值分析理论和方法的理解,我们进行了数值分析实验,通过实验加深对理论知识的掌握,提高实际操作能力。
二、实验目的1. 理解数值分析的基本理论和方法;2. 掌握数值分析实验的基本步骤和技巧;3. 培养实验设计和数据分析能力;4. 提高编程和计算能力。
三、实验内容本次实验主要分为以下几个部分:1. 线性方程组求解实验:通过高斯消元法、LU分解法等求解线性方程组,并分析算法的稳定性和误差;2. 矩阵特征值问题计算实验:利用幂法、逆幂法等计算矩阵的特征值和特征向量,分析算法的收敛性和精度;3. 非线性方程求根实验:运用二分法、牛顿法、不动点迭代法等求解非线性方程的根,比较不同算法的优缺点;4. 函数插值实验:运用拉格朗日插值、牛顿插值等方法对给定的函数进行插值,分析插值误差;5. 常微分方程初值问题数值解法实验:运用欧拉法、改进的欧拉法、龙格-库塔法等求解常微分方程初值问题,比较不同算法的稳定性和精度。
四、实验过程1. 线性方程组求解实验:首先,编写程序实现高斯消元法、LU分解法等算法;然后,对给定的线性方程组进行求解,记录计算结果;最后,分析算法的稳定性和误差。
2. 矩阵特征值问题计算实验:编写程序实现幂法、逆幂法等算法;然后,对给定的矩阵进行特征值和特征向量的计算,记录计算结果;最后,分析算法的收敛性和精度。
3. 非线性方程求根实验:编写程序实现二分法、牛顿法、不动点迭代法等算法;然后,对给定的非线性方程进行求根,记录计算结果;最后,比较不同算法的优缺点。
4. 函数插值实验:编写程序实现拉格朗日插值、牛顿插值等方法;然后,对给定的函数进行插值,记录计算结果;最后,分析插值误差。
5. 常微分方程初值问题数值解法实验:编写程序实现欧拉法、改进的欧拉法、龙格-库塔法等算法;然后,对给定的常微分方程初值问题进行求解,记录计算结果;最后,比较不同算法的稳定性和精度。
数值分析学习公式总结

第一章1霍纳(Horner )方法: n a 1-n a 2-n a ……2a 1a 0a输入=c+ n b *c c b n *1- c b *3 c b *2 c b *1n b 1-n b 2-n b 2b 1b 0bAnswer P (x )=0b该方法用于解决多项式求值问题P (x )=n a n x +1-n a 1-n x +2-n a 2-n x +……+2a 2x +1a x +0a2 注:p ˆ为近似值绝对误差:|ˆ|pp E p -=相对误差:|||ˆ|p pp R p -=有效数字:210|||ˆ|1d p p pp R -<-= (d 为有效数字,为满足条件的最大整数) 3 Big Oh(精度的计算): O(h ⁿ)+O(h ⁿ)=O(h ⁿ);O(h m )+O(h n )=O(h r ) [r=min{p,q}]; O(h p )O(h q )=O(h s ) [s=q+p]; 第二章2.1 求解x=g(x)的迭代法 用迭代规则,可得到序列值{}。
设函数g 。
如果对于所有x ,映射y=g(x)的范围满足y , 则函数g 在内有一个不动点; 此外,设定义在内,且对于所有x ,存在正常数K<1,使得,则函数g 在内有唯一的不动点P 。
定理2.3 设有(i )g ,g ’,(ii )K 是一个正常数,(iii )。
如果对于所有如果对于所有x 在这种情况下,P 成为排斥不动点,而且迭代显示出局部发散性。
. 波尔查诺二分法(二分法定理)<收敛速度较慢>试值(位)法:<条件与二分法一样但改为寻求过点(a,f(a))和(b,f(b))的割线L 与x 轴的交点(c,0)>应注意越来越小,但可能不趋近于0,所以二分法的终止判别条件不适合于试值法.牛顿—拉夫森迭代函数:)(')()(1111-----==k k k k k p f p f p p g p 其中k=1,2,……证明:用泰勒多项式证明第三章线性方程组的解法 对于给定的解线性方程组Ax=b一Gauss Elimination (高斯消元法 )第一步Forward Elimination 第二步 BackSubstitution二LU Factorization第一步 A = LU 原方程变为LUx=y ;第二步 令Ux=y,则Ly = b 由下三角解出y ; 第三步 Ux=y,又上三角解出x ;三Iterative Methods (迭代法)2n n 22221211n n 1212111b x a x a x a b x a x a x a =+++=+++nn nn 22n 11n 2n n 22221211n n 1212111b x a x a x a b x a x a x a b x a x a x a =+++=+++=+++初始值四 Jacobi Method1.选择初始值2.迭代方程为五Gauss Seidel Method1.迭代方程为00201,,,n x x x 00201,,,n x x x nnk n nn k n k n n k n k nn k k kn n k k a x a x a x a bx a x a x a bx a x a x a b x )()()(1122111222121212111212111--++++++-=++-=++-=k k k kn n k k kn n k k a x a x a bx a x a x a bx )()(1112221121212111212111++++++++-=++-=2.选择初始值 判断是否能用Jacobi Method 或者GaussSeidel Method 的充分条件(绝对对角占优原则)第四章 插值与多项式逼近·第一节 泰勒级数和函数计算一些常用函数的泰勒级数展开:for all x for all x for all x -1 -1for00201,,,nx x x定理4.1(泰勒多项式逼近)设,而是固定值。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
工程硕士《数值分析》总复习题(2011年用)[由教材中的习题、例题和历届考试题选编而成,供教师讲解和学生复习用]一. 解答下列问题:1)下列所取近似值有多少位有效数字( 注意根据什么? ):a) 对 e = 2.718281828459045…,取*x = 2.71828b) 数学家祖冲之取 113355作为π的近似值.c) 经过四舍五入得出的近似值12345,-0.001, 90.55000, 它们的有效 数字位数分别为 位, 位, 位。
2) 简述下名词:a) 截断误差 (不超过60字) b) 舍入误差 (不超过60字)c) 算法数值稳定性 (不超过60字)3) 试推导( 按定义或利用近似公式 ): 计算3x 时的相对误差约等于x 的相对误差的3倍。
4) 计算球体积334r Vπ= 时,为使其相对误差不超过 0.3%,求半径r 的相对 误差的允许范围。
5) 计算下式3418)1(3)1(7)1(5)1(22345+-+---+---=x x x x x x P)(时,为了减少乘除法次数, 通常采用什么算法? 将算式加工成什么形式?6) 递推公式 ⎪⎩⎪⎨⎧=-==-,2,1,110210n y y y n n如果取*041.12y y =≈= ( 三位有效数字 ) 作近似计算, 问计算到10y 时误差为初始误差的多少倍? 这个计算过程数值稳定吗 ?二. 插值问题:1) 设函数)(x f 在五个互异节点 54321,,,,x x x x x 上对应的函数值为54321,,,,f f f f f ,根据定理,必存在唯一的次数 (A ) 的插值多项式)(x P ,满足插值条件 ( B ) . 对此,为了构造Lagrange 插值多项式 )(x L ,由5个节点作 ( C ) 个、次数均为 ( D ) 次的插值基函数)(x l i = _(E ) , 从而得Lagrange 插值多项式)(x L = (F ) ,而插值余项 )()()(x L x f x R -== (G ) 。
2 ) 试用三种方法求过三个离散点:A (0,1) 、B (1,2) 、C (2,3) 的插值多项式。
3) 求函数x e x f -=)( 在 [ 0 , 1 ]上的近似一次插值多项式。
4 ) 由函数值表: x : 1 2 3 x e - : 0.367879441 , 0.135335283 , 0.049787068求1.2-e的近似值.5) 利用插值方法推导 x i j i jx ni nij j =--∑∏=≠=][,0 三. 拟合问题:1) 对离散实验数据做最小二乘拟合的两个主要步骤是 ( A ) 和 ( B ) . 2) 对同一个量的多个近似值, 常取其算术平均作为该量的近似值, 这种做法的意义是什么? 3) 设有实验数据如下: x 1.36 1.73 1.95 2.28f14.094 16.844 18.475 20.963按最小二乘法求其拟合曲线。
4) 已知某试验过程中函数f 依赖于x 的试验数据如下:i x : 1 2 3 4i f : 0.8 1.5 1.8 2.0试按最小二乘法拟合出一个形如 2bx ax S += 的经验公式。
5 ) 设有实验数据如下: x 1 2 3 4f4 10 18 26按最小二乘法拟合出一个形如 2bx a S += 的经验公式 。
四. 数值求积:1) 写出数值求积公式的一般形式, 指出其特点, 并说明它对计算机的计算有什么意义?2) 简述数值求积公式的 ”代数精度” 的概念 3) 插值型求积公式()() nbk k ak f x dx A f x =≈∑⎰中,每个系数可用公式k A =( A ) 计算,它们之和∑=nk kA= ( B ) , 其代数精度 ( C ) .又Newton-Cotes 公式的一般形式为 ( D ) , 其主要特点是 ( E ) , 其 Cotes 系数之和∑=nk n kC)(= ( F ) , 其代数精度 ( G ) ;4) 考察数值求积公式⎰--++-≈11101)1()0()1()(f A f A f A dx x f ,直接指出: 它是什么类型的公式? 为使其精度尽可能高,101,,A A A -应取什么确值? 它是不是Gauss 型公式?5 ) 求dx xI⎰+=10311的近似值, 试写出使用11个等分点函数值的求积 公式( 要求只列出数值公式,不需要求出具体结果 )。
6 ) 利用复化Simpson 公式求积分 dx x I⎰=21的近似值(只需列出算式) 。
7) 利用现成函数表,分别用复化梯形公式n T 和复化Simpson 公式n S 计算积分ϕϕπd I ⎰-=62sin 4 ϕ ϕ2sin 4-0 2π 9981001.1362π 9924473.1 363π 9831825.1364π 9705386.1365π 9548386.1366π9364917.1五. 解线性代数方程组的直接法:1) Gauss 消去过程中引入选主元技巧的目的是下列中的哪一项或哪几项?A .提高计算速度;B .提高计算精度;C .简化计算公式;D .提高计算公式的数值稳定性;E .节省存储空间。
2) 采用“列主元Gauss 消去法” 解下列方程组:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡565331743532321x x x a) 用 ”列主元Gauss 消去过程” 将方程组约化成上三角方程组; b) 用 ”回代过程” 依次列式计算出方程组的解。
3) 设方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---6745150710623321x x x 现采用“列主元Gauss 消去法”求解,试回答: a ) 所用列主元Gauss 消去法包括哪两个过程?b ) 要用几步消元?c ) 每一步消元计算之前需做哪些工作(用简短、准确的文字叙述)?d ) 现经第1步消元结果, 上述方程组已被约化为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--251061321251017560710x x x 请你继续做消元计算, 直至约化成上三角方程组。
e )对所得上三角方程组依次列式计算出方程组的解。
六. 解线性代数方程组的迭代法: 1) 解线性代数方程组 fx B x+= 的基本型迭代公式,1,0,)()1(=+=+k f x B xk k其中B 称为什么? )0(x又称为什么? 如果迭代序列{})(k x 有极限*x(即迭代公式收敛),则极限*x 是什么? 2) 设解线性代数方程组b Ax =(其中n n R A ⨯∈非奇异,0≠b )的迭代公式为,1,0,)()()()1(=--=+k b Ax x xk k k λ则其迭代矩阵是什么? 此迭代公式对任意的初始向量)0(x 收敛的充分必要条件是什么? 又此迭代公式对任意的初始向量)0(x 收敛的一个充分条件是什么?3) 设线性方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡53411221x x , 试构造解此方程组的Jacobi 迭代公式和GS 迭代公式; 试问所作的两种迭代公式是否收敛,为什么? 试用初值 T x )0,0()0(= 计算GS 迭代公式的前三个值.4 ) 设方程组⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡--84195121x x 试构造解此方程组的收敛的Jacobi 迭代公式和收敛的Guass-Seidel 迭代公式, 并说明两者收敛的根据; 求出这两种迭代的迭代矩阵. 5) 设线性方程组3,,15.05.025.05.01,R b x a a A b Ax ∈⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----==请按便于计算的收敛充分条件, 求使J 法和GS 法均收敛的 a 的取值范围.七.一元方程求根: 1) 写出求方程013)(3=--=x x x f 在 [ 1,2 ]中的近似根的一个收敛的不动点迭代公式,并证明其收敛性。
2) 已知方程 )1(2ln>=-x x x 的有根区间 [ 3,4 ] .试写出求该方程在 [ 3 , 4 ] 中的根的一个不动点迭代公式; 证明所给出 的迭代公式是收敛的。
试设计其计算机算法. 3) 用Newton 迭代法求方程013)(3=--=x x x f 在20=x 附近的根,试写其Newton 迭代公式; 并说明其收敛情况。
4)的Newton 迭代公式,并说明其收敛情况。
八. 常微分方程初值问题:1) 常微分方程定解问题分为初值问题和 ( A ) 问题.初值问题是指由 (B ) 和(C ) 两部分联立起来构成的问题。
研究常微分方程初值问题时, 通常针对基本形式 (D ) 进行研究。
设函数)(x y 是某初值问题的解析解, 则该初值问题在n x 处的解为 ( E ) 而数值解(通常记)为 (F ) ,它们的关系是 ( G ) .若记)(1+n x y 是初值问题在点1+n x 处的解, 1+n y 是由某数值方法得出的1+n x 处的数值解,则该数值方法在1+n x 处的局部截断误差是指 (H ) .2) 设初值问题 ⎩⎨⎧=≤≤--='1)0(6.00,2y x y y x y试用Euler 方法取2.0=h ,求解上述初值问题的数值解。
3 ) 设初值问题 ⎩⎨⎧=≤≤-='2)1(21,38y x y y试用梯形方法求其解在两点 4.1,2.1=x 处的值)4.1(,)2.1(y y 的近似值。
4) 设初值问题 ⎩⎨⎧=<<++='1)0(10,122y x x y y试用改进的Euler 方法,并取1.0=h,设计一个求解上述初值问题数值解的求解方案 (或称计算机算法描述; 不必求出解的具体数值) 。
5 ) 设初值问题 ⎩⎨⎧=≤≤+='1)0(10),1/(3y x x y y试用4阶经典R-K 方法,并取1.0=h,设计一个求解上述初值问题数值解的求解方案 (或称计算机算法描述; 不必求出解的具体数值) 。
九、下列各小题任选其中已学过的小题作练习: 1) 设T x )3,2,0(=, 求,1x,2x,∞x;设 ⎥⎦⎤⎢⎣⎡=4321A ,求1A,∞A,2A, )(A ρ。
2 ) 用较简捷的方法分别求下列的插值多项式)(x H 和)(x p ,并写出其余项公式:a) 1)1(,0)0()0(,1)1(=='=-=-H H H H b) 2)2(,0)1()1(,1)0(=='==p p p p3 ) 用插值方法求在0=x处与x cos 相切 ,在2π=x 处与x cos 相交的二次多项式)(2x p ,并推导插值余项的估计式为|2|61)(cos 22π-≤-x x x p x 4 ) 试用最小二乘法原理求下列超定方程组的近似解:⎪⎪⎩⎪⎪⎨⎧=+=+=-=+7262353114221212121x x x x x x x x5 ) 要计算函数dt e x y xt ⎰-=02)( 在x = 0.2, 0.4, 0.6 三处的近似值,试用解初值问题的数值方法,设计其计算方案 (要求采用二阶精度的计算公式).6) 用追赶法解三对角方程组: ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡022112111131124321x x x x 7) 对方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡==21,13.021,b A b Ax 拟用迭代法,1,0,)()()()1(=-+=+k b Ax x x k k k α求解, 试确定 α 的取值范围,使得上述迭代公式收敛.8) 对迭代函数2()(5)x x x ϕλ=+-,试求使迭代公式,1,0),(1==+k x x k k ϕ,局部收敛于5=*x 的λ的取值范围。