数值分析复习题要答案

合集下载

数值分析复习题及答案

数值分析复习题及答案

数值分析复习题一、选择题1. 3.142和3.141分别作为π的近似数具有〔 〕和〔 〕位有效数字.A .4和3B .3和2C .3和4D .4和42. 求积公式()()211211()(2)636f x dx f Af f ≈++⎰,那么A =〔 〕A . 16B .13C .12D .233. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足〔 〕A .()00l x =0,()110l x = B .()00l x =0,()111l x = C .()00l x =1,()111l x = D . ()00l x =1,()111l x =4. 设求方程()0f x =的根的牛顿法收敛,那么它具有〔 〕敛速。

A .超线性B .平方C .线性D .三次5. 用列主元消元法解线性方程组1231231220223332x x x x x x x x ++=⎧⎪++=⎨⎪--=⎩ 作第一次消元后得到的第3个方程〔 〕.A .232x x -+= B .232 1.5 3.5x x -+= C .2323x x -+= D .230.5 1.5x x -=-二、填空1. 设2.3149541...x *=,取5位有效数字,那么所得的近似值x=.2.设一阶差商()()()21122114,321f x f x f x x x x --===---,()()()322332615,422f x f x f x x x x --===-- 那么二阶差商()123,,______f x x x =3. 设(2,3,1)TX =--, 那么2||||X = ,=∞||||X 。

4.求方程 21.250x x --= 的近似根,用迭代公式x =01x =, 那么 1______x =。

5.解初始值问题 00'(,)()y f x y y x y =⎧⎨=⎩近似解的梯形公式是 1______k y +≈。

数值分析期末考试和答案

数值分析期末考试和答案

数值分析期末考试和答案一、单项选择题(每题2分,共20分)1. 在数值分析中,下列哪个方法用于求解线性方程组?A. 插值法B. 迭代法C. 直接法D. 拟合法答案:C2. 以下哪个数值方法是用于求解非线性方程的?A. 高斯消元法B. 牛顿迭代法C. 线性插值法D. 拉格朗日插值法答案:B3. 在数值积分中,梯形法则的误差与下列哪个因素无关?A. 被积函数的二阶导数B. 积分区间的长度C. 积分区间的划分数量D. 被积函数的一阶导数答案:D4. 以下哪个数值方法是用于求解常微分方程的?A. 欧拉方法B. 牛顿迭代法C. 拉格朗日插值法D. 高斯消元法答案:A5. 在数值分析中,下列哪个方法用于求解特征值问题?A. 高斯消元法B. 幂迭代法C. 牛顿迭代法D. 梯形法则答案:B6. 以下哪个数值方法是用于求解线性最小二乘问题的?A. 高斯消元法B. 梯形法则C. 正交分解法D. 牛顿迭代法答案:C7. 在数值分析中,下列哪个方法用于求解非线性方程组?A. 高斯消元法B. 牛顿迭代法C. 线性插值法D. 欧拉方法答案:B8. 在数值分析中,下列哪个方法用于求解偏微分方程?A. 有限差分法B. 牛顿迭代法C. 线性插值法D. 梯形法则答案:A9. 在数值分析中,下列哪个方法用于求解优化问题?A. 高斯消元法B. 梯形法则C. 牛顿迭代法D. 单纯形法答案:D10. 在数值分析中,下列哪个方法用于求解插值问题?A. 高斯消元法B. 梯形法则C. 牛顿迭代法D. 拉格朗日插值法答案:D二、填空题(每题2分,共20分)1. 在数值分析中,求解线性方程组的直接法包括______消元法和______消元法。

答案:高斯;LU2. 牛顿迭代法的收敛速度是______阶的。

答案:二3. 梯形法则的误差与被积函数的______阶导数有关。

答案:二4. 欧拉方法是一种求解______阶常微分方程的数值方法。

答案:一5. 幂迭代法是求解______特征值问题的数值方法。

数值分析复习试题及参考答案

数值分析复习试题及参考答案

1、设115.80,1025.621≈≈x x 均具有5位有效数字,试估计由这些数据计算21x x 的绝对误差限。

解:由有效数字的定义得,()104121-⨯=x ε,()103221-⨯=x ε()07057.00005.0115.80005.01025.621=⨯+⨯≈x x ε2、设430.56,1021.12≈≈x x均具有5位有效数字,试估计由这些数据计算21x x +的绝对误差限。

解:由有效数字的定义得,()104121-⨯=x ε,()103221-⨯=x ε0055.0)()()(2121=+=+x x x x εεε3、简答题 (1)已知12622)(256+-+-=x xxxx f ,求]1,0[f 及]6,5,4,3,2,1,0[f 。

解:由f(0)=1,f(1)=5得 []()()41011,0=-=f f f因为最高阶差商只出现在最高次,所以[]26,5,4,3,2,1,0=f(2)求积公式[])1()0(121)]1()0([21)(1f f f f dx x f '-'++≈⎰的代数精度为多少? 解:令()xx f =,则()21211021==⎰xdx x f ,右边=21,左边=右边同理令()2xx f =,()3xx f =均准确成立,()4xx f =时,左边≠右边所以,上式具有3阶精度4、求满足下表条件的Hermit 插值多项式。

x0 1)(x f -1 0 )(x f '-210解:使用重节点差商表法x y 一阶二阶 三阶 0 -1 0 -1 -2 1 0 1 3 1 010 9 6()()1236163212322---=-++--=x x x x xx x x H5、已知函数)(x f y =的数据如下:x1 2 4 -5 )(x f3 4 1 0(1)求3次Lagrange 插值多项式; (2)求3次Newton 插值多项式; (3)写出插值余项。

数值分析期末试题及答案

数值分析期末试题及答案

数值分析期末试题及答案试题一:1. 简答题(共10分)a) 什么是数值分析?它的主要应用领域是什么?b) 请简要解释迭代法和直接法在数值计算中的区别。

2. 填空题(共10分)a) 欧拉方法是一种______型的数值解法。

b) 二分法是一种______法则。

c) 梯形法则是一种______型的数值积分方法。

3. 计算题(共80分)将以下函数进行数值求解:a) 通过使用二分法求解方程 f(x) = x^3 - 4x - 9 = 0 的近似解。

b) 利用欧拉方法求解微分方程 dy/dx = x^2 + 2x + 1, y(0) = 1 在 x = 1 处的解。

c) 使用梯形法则计算积分∫[0, π/4] sin(x) dx 的近似值。

试题二:1. 简答题(共10分)a) 请解释什么是舍入误差,并描述它在数值计算中的影响。

b) 请解释牛顿插值多项式的概念及其应用。

2. 填空题(共10分)a) 数值稳定性通过______号检查。

b) 龙格-库塔法是一种______计算方法。

c) 零点的迭代法在本质上是将方程______转化为______方程。

3. 计算题(共80分)使用牛顿插值多项式进行以下计算:a) 已知插值节点 (-2, 1), (-1, 1), (0, 2), (1, 4),求在 x = 0.5 处的插值多项式值。

b) 已知插值节点 (0, 1), (1, 2), (3, 7),求插值多项式,并计算在 x = 2 处的值。

c) 使用 4 阶龙格-库塔法求解微分方程 dy/dx = x^2 + 1, y(0) = 1。

答案:试题一:1. a) 数值分析是研究使用数值方法解决数学问题的一门学科。

它的主要应用领域包括数值微积分、数值代数、插值和逼近、求解非线性方程、数值积分和数值解微分方程等。

b) 迭代法和直接法是数值计算中常用的两种方法。

迭代法通过反复迭代逼近解,直到满足所需精度为止;而直接法则通过一系列代数运算直接得到解。

数值分析期末试题及答案

数值分析期末试题及答案

数值分析期末试题及答案一、选择题(每题5分,共20分)1. 在数值分析中,下列哪个算法不是用于求解线性方程组的?A. 高斯消元法B. 牛顿法C. 雅可比法D. 追赶法答案:B2. 插值法中,拉格朗日插值法属于:A. 多项式插值B. 样条插值C. 线性插值D. 非线性插值答案:A3. 以下哪个选项不是数值分析中的误差来源?A. 截断误差B. 舍入误差C. 计算误差D. 测量误差答案:C4. 在数值积分中,梯形法则的误差项是:A. O(h^2)B. O(h^3)C. O(h)D. O(1)答案:A二、填空题(每题5分,共20分)1. 牛顿插值法中,插值多项式的一般形式为:______。

答案:f(x) = a_0 + a_1(x-x_0) + a_2(x-x_0)(x-x_1) + ...2. 牛顿迭代法求解方程的根时,迭代公式为:x_{n+1} = x_n -f(x_n) / __________。

答案:f'(x_n)3. 在数值分析中,______ 用于衡量函数在区间上的近似积分值与真实积分值之间的差异。

答案:误差4. 线性方程组的解法中,______ 法是利用矩阵的LU分解来求解。

答案:克兰特三、解答题(每题10分,共60分)1. 给定函数f(x) = e^(-x),使用拉格朗日插值法,求x = 0.5时的插值值。

解答:首先选取插值节点x_0 = 0, x_1 = 0.5, x_2 = 1,对应的函数值分别为f(0) = 1, f(0.5) = e^(-0.5), f(1) = e^(-1)。

拉格朗日插值多项式为:L(x) = f(0) * (x-0.5)(x-1) / (0-0.5)(0-1) + f(0.5) * (x-0)(x-1) / (0.5-0)(0.5-1) + f(1) * (x-0)(x-0.5) / (1-0)(1-0.5)将x = 0.5代入得:L(0.5) = 1 * (0.5-0.5)(0.5-1) / (0-0.5)(0-1) + e^(-0.5) * (0.5-0)(0.5-1) / (0.5-0)(0.5-1) + e^(-1) * (0.5-0)(0.5-0.5) / (1-0)(1-0.5)计算得L(0.5) = e^(-0.5)。

(完整版)数值分析整理版试题及答案,推荐文档

(完整版)数值分析整理版试题及答案,推荐文档

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. 以下哪个算法是用于求解线性方程组的直接方法?A. 牛顿法B. 高斯消元法C. 共轭梯度法D. 辛普森积分法答案:B2. 插值法中,拉格朗日插值法和牛顿插值法的主要区别是什么?A. 插值点的选取不同B. 插值多项式的构造方式不同C. 计算复杂度不同D. 适用的函数类型不同答案:B3. 在数值积分中,梯形法则和辛普森法则的主要区别是什么?A. 精度不同B. 适用的积分区间不同C. 计算方法不同D. 稳定性不同答案:A二、简答题1. 解释什么是数值稳定性,并举例说明。

答案:数值稳定性指的是数值方法在计算过程中对于舍入误差的敏感程度。

例如,在求解线性方程组时,如果系数矩阵的条件数很大,则该方程组的数值解对舍入误差非常敏感,即数值稳定性差。

2. 说明数值微分与数值积分的区别。

答案:数值微分是估计函数在某一点的导数,而数值积分是估计函数在某个区间上的积分。

数值微分通常用于求解函数的局部变化率,而数值积分用于求解函数在一定区间内的累积效果。

三、计算题1. 给定一组数据点:(1, 2), (2, 3), (3, 5), (4, 6),请使用拉格朗日插值法构造一个三次插值多项式。

答案:首先写出拉格朗日插值基函数,然后根据数据点构造插值多项式。

具体计算过程略。

2. 给定函数 f(x) = x^2,使用牛顿-科特斯公式中的辛普森积分法在区间 [0, 1] 上估计积分值。

答案:首先确定区间划分,然后应用辛普森积分公式进行计算。

具体计算过程略。

四、论述题1. 论述数值分析中误差的来源及其控制方法。

答案:误差主要来源于舍入误差和截断误差。

舍入误差是由于计算机在进行浮点数运算时的精度限制造成的,而截断误差是由于数值方法的近似性质导致的。

控制误差的方法包括使用高精度的数据类型、选择合适的数值方法、增加计算步骤等。

五、综合应用题1. 给定一个线性方程组 Ax = b,其中 A 是一个 3x3 的矩阵,b 是一个列向量。

数值分析整理版试题及答案

数值分析整理版试题及答案

例1、 已知函数表求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式. 解:(1)插值基函数分别为()()()()()()()()()()1200102121()1211126x x x x x x l x x x x x x x ----===--------()()()()()()()()()()021*******()1211122x x x x x x l x x x x x x x --+-===-+---+-()()()()()()()()()()0122021111()1121213x x x x x x l x x x x x x x --+-===-+--+-故所求二次拉格朗日插值多项式为()()()()()()()()()()()2202()11131201241162314121123537623k k k L x y l x x x x x x x x x x x x x ==⎡⎤=-⨯--+⨯-+-+⨯+-⎢⎥⎣⎦=---++-=+-∑(2)一阶均差、二阶均差分别为[]()()[]()()[][][]010*********011201202303,11204,41234,,52,,126f x f x f x x x x f x f x f x x x x f x x f x x f x x x x x ---===-----===----===---故所求Newton 二次插值多项式为()()[]()[]()()()()()20010012012,,,35311126537623P x f x f x x x x f x x x x x x x x x x x x =+-+--=-++++-=+-例2、 设2()32f x xx =++,[0,1]x ∈,试求()f x 在[0, 1]上关于()1x ρ=,{}span 1,x Φ=的最佳平方逼近多项式.解:若{}span 1,x Φ=,则0()1x ϕ=,1()x x ϕ=,且()1x ρ=,这样,有()()()()()()()()1120011011201100012101,11,,3123,,,,32269,324dx x dx xdx f x x dx f x x x dx ϕϕϕϕϕϕϕϕϕϕ========++==++=⎰⎰⎰⎰⎰ 所以,法方程为01123126119234a a ⎡⎤⎡⎤⎢⎥⎢⎥⎡⎤=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦⎣⎦,经过消元得01231162110123a a ⎡⎤⎡⎤⎢⎥⎢⎥⎡⎤=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦⎣⎦再回代解该方程,得到14a =,0116a =故,所求最佳平方逼近多项式为*111()46S x x =+ 例3、 设()xf x e =,[0,1]x ∈,试求()f x 在[0, 1]上关于()1x ρ=,{}span 1,x Φ=的最佳平方逼近多项式. 解:若{}span 1,x Φ=,则0()1x ϕ=,1()x x ϕ=,这样,有()()()()()()100012110101100100110,111,31,,2, 1.7183,1x x dx x dx xdx f e dx f xe dx ϕϕϕϕϕϕϕϕϕϕ===========⎰⎰⎰⎰⎰所以,法方程为0111 1.7183211123a a ⎡⎤⎢⎥⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎢⎥⎣⎦解法方程,得到00.8732a =,1 1.6902a =, 故,所求最佳平方逼近多项式为*1()0.8732 1.6902S x x =+例4、 用4n =的复合梯形和复合辛普森公式计算积分1⎰。

数值分析期末考试题及答案

数值分析期末考试题及答案

数值分析期末考试题及答案一、选择题(每题2分,共20分)1. 在数值分析中,下列哪个算法用于求解线性方程组?A. 牛顿法B. 高斯消元法C. 插值法D. 傅里叶变换答案:B2. 以下哪个选项不是数值分析中的误差类型?A. 舍入误差B. 截断误差C. 测量误差D. 累积误差答案:C3. 多项式插值中,拉格朗日插值法的特点是:A. 插值点必须等距分布B. 插值多项式的次数与插值点的个数相同C. 插值多项式是唯一的D. 插值多项式在插值点处的值都为1答案:B4. 在数值分析中,下列哪个方法用于求解非线性方程?A. 辛普森法则B. 牛顿迭代法C. 欧拉法D. 龙格-库塔法答案:B5. 以下哪个是数值稳定性的指标?A. 收敛性B. 收敛速度C. 条件数D. 误差传播答案:C二、简答题(每题10分,共20分)1. 简述高斯消元法求解线性方程组的基本原理。

答案:高斯消元法是一种直接解法,通过行变换将增广矩阵转换为上三角形式,然后通过回代求解线性方程组。

它包括三个基本操作:行交换、行乘以非零常数、行相加。

2. 解释什么是数值稳定性,并举例说明。

答案:数值稳定性是指数值解对输入数据小的扰动不敏感的性质。

例如,某些数值方法在计算过程中可能会放大舍入误差,导致结果不可靠,这样的方法就被认为是数值不稳定的。

三、计算题(每题15分,共30分)1. 给定线性方程组:\[\begin{align*}x + 2y - z &= 4 \\3x - y + 2z &= 1 \\-x + y + z &= 2\end{align*}\]使用高斯消元法求解该方程组,并给出解。

答案:首先将增广矩阵转换为上三角形式,然后回代求解,得到\( x = 1, y = 2, z = 1 \)。

2. 给定函数 \( f(x) = x^2 - 3x + 2 \),使用拉格朗日插值法在\( x = 0, 1, 2 \) 处插值,并求出插值多项式。

数值分析考试题和答案

数值分析考试题和答案

数值分析考试题和答案一、单项选择题(每题2分,共20分)1. 在数值分析中,插值法的主要目的是()。

A. 求解线性方程组B. 求解非线性方程C. 构造一个多项式来近似一个函数D. 求解微分方程答案:C2. 线性方程组的高斯消元法中,主元为零时,应采取的措施是()。

A. 停止计算B. 回代求解C. 转置矩阵D. 行交换答案:D3. 以下哪种方法不是数值积分方法()。

A. 梯形规则B. 辛普森规则C. 牛顿法D. 复合梯形规则答案:C4. 以下哪种方法用于求解非线性方程的根()。

A. 欧几里得算法B. 牛顿迭代法C. 高斯消元法D. 线性插值法答案:B5. 在数值分析中,最小二乘法主要用于()。

A. 求解线性方程组B. 求解非线性方程C. 曲线拟合D. 微分方程数值解答案:C6. 以下哪种方法不是数值微分方法()。

A. 前向差分B. 后向差分C. 中心差分D. 欧拉方法答案:D7. 以下哪种方法用于求解常微分方程的初值问题()。

A. 欧拉方法B. 龙格-库塔方法C. 牛顿迭代法D. 高斯消元法答案:B8. 在数值分析中,矩阵的特征值问题可以通过()方法求解。

A. 高斯消元法B. 幂迭代法C. 牛顿迭代法D. 梯形规则答案:B9. 以下哪种方法不是数值稳定性分析中的方法()。

A. 绝对稳定性B. 相对稳定性C. 条件数D. 牛顿法答案:D10. 在数值分析中,条件数用于衡量()。

A. 算法的效率B. 算法的稳定性C. 算法的准确性D. 算法的复杂度答案:B二、填空题(每题2分,共20分)1. 在数值分析中,插值多项式的次数最高为______,其中n是插值点的个数。

答案:n-12. 线性方程组的高斯消元法中,如果某行的主元为零,则需要进行______。

答案:行交换3. 梯形规则的误差与被积函数的______阶导数有关。

答案:二4. 牛顿迭代法中,每次迭代需要计算______。

答案:函数值和导数值5. 最小二乘法中,残差平方和最小化时,对应的系数向量是______。

数值分析试题及答案

数值分析试题及答案

数值分析试题及答案一、选择题(每题2分,共20分)1. 以下哪个算法是数值分析中用于求解线性方程组的直接方法?A. 牛顿法B. 高斯消元法C. 梯度下降法D. 蒙特卡洛方法答案:B2. 插值法中,拉格朗日插值法和牛顿插值法的共同点是:A. 都是多项式插值B. 都使用差商C. 都只适用于等距节点D. 都需要预先知道所有数据点答案:A3. 在数值积分中,辛普森(Simpson)公式比梯形公式的误差:A. 更大B. 更小C. 相同D. 无法比较答案:B4. 以下哪个是数值稳定性分析中常用的方法?A. 条件数B. 收敛性C. 收敛速度D. 误差分析答案:A5. 在求解常微分方程的数值解时,欧拉方法属于:A. 单步法B. 多步法C. 隐式方法D. 显式方法答案:A6. 以下哪个是数值分析中求解非线性方程的迭代方法?A. 高斯-约当消元法B. 牛顿-拉弗森方法C. 雅可比迭代法D. 高斯-赛德尔迭代法答案:B7. 线性插值公式中,如果给定两个点\( (x_0, y_0) \)和\( (x_1, y_1) \),插值多项式是:A. \( y = y_0 + \frac{y_1 - y_0}{x_1 - x_0}(x - x_0) \)B. \( y = y_0 + \frac{y_1 - y_0}{x_0 - x_1}(x - x_0) \)C. \( y = y_0 + \frac{x - x_0}{x_1 - x_0}(y_1 - y_0) \)D. \( y = y_1 + \frac{x_1 - x}{x_1 - x_0}(y_0 - y_1) \)答案:C8. 以下哪个是数值分析中用于求解特征值问题的算法?A. 幂法B. 共轭梯度法C. 牛顿法D. 欧拉法答案:A9. 在数值微分中,使用有限差分法来近似导数时,中心差分法的误差:A. 与步长成正比B. 与步长的平方成正比C. 与步长的立方成正比D. 与步长的四次方成正比答案:B10. 以下哪个是数值分析中用于求解线性最小二乘问题的算法?A. 梯度下降法B. 牛顿法C. 奇异值分解法D. 共轭梯度法答案:C二、简答题(每题10分,共30分)1. 简述数值分析中病态问题的特点及其对算法的影响。

数值分析练习题附答案

数值分析练习题附答案

目录一、绪论------------------------------------------------------------------------------------- 2-2二、线性方程组直接解法列主元高斯LU LDL T GG T-------------------- 3-6二、线性方程组迭代法----------------------------------------------------------------- 7-10 三、四、非线性方程组数值解法二分法不动点迭代---------------------- 11-13五、非线性方程组数值解法牛顿迭代下山弦截法----------------- 14-15六、插值线性插值抛物线插值------------------------------------------------ 16-18七、插值Hermite插值分段线性插值-----------------------------------------19-22八、拟合------------------------------------------------------------------------------------ 23-24九、数值积分----------------------------------------------------------------------------- 25-29十、常微分方程数值解法梯形欧拉改进----------------------------------- 30-32 十一、常微分方程数值解法龙格库塔------------------------------------------ 33-35绪论1-1 下列各数都是经过四舍五入得到的近似值 ,试分别指出它们的绝对误差限,相对误差限和有效数字的位数.X 1 =5.420, X 2 =0.5420, X 3 =0.00542, X 4 =6000, X 5 =0.6×105注:将近似值改写为标准形式X 1 =(5*10-1+4*10-2+2*10-3+0*10-4)*101 即n=4,m=1 绝对误差限|△X 1|=|X *1-X 1|≤ 12×10m-n =12×10-3 相对误差限|△r X 1|= |X∗1−X1||X∗1|≤|X∗1−X1||X1|= 12×10-3/5.4201-2 为了使101/2 的相对误差小于0.01%, 试问应取几位有效数字?1-3 求方程x 2 -56x+1=0的两个根, 使它们至少具有4位有效数字( √783≈27.982)注:原方程可改写为(x-28)2=783线性方程组解法(直接法)2-1用列主元Gauss消元法解方程组解:回代得解:X1=0 X2=-1 X3=12-2对矩阵A进行LU分解,并求解方程组Ax=b,其中解:(注:详细分解请看课本P25)A=(211132122)→(211(1/2)5/23/2(1/2)3/23/2)→(2111/25/23/21/2(3/5)3/5)即A=L×U=(11/211/23/51)×(2115/23/23/5)先用前代法解L y=P b 其中P为单位阵(原因是A矩阵未进行行变换)即L y=P b 等价为(11/211/23/51)(y1y2y3)=(111)(465)解得 y 1=4 y 2=4 y 3=35再用回代解Ux =y ,得到结果x即Ux =y 等价为(2115/23/23/5)(x 1x 2x 3)=(y 1y 2y 3)=(443/5) 解得 x 1=1 x 2=1 x 3=1即方程组Ax=b 的解为x =(111)2-3 对矩阵A 进行LDL T 分解和GG T 分解,求解方程组Ax=b,其中A=(164845−48−422) , b =(123)解:(注:课本 P 26 P 27 根平方法)设L=(l i j ),D=diag(d i ),对k=1,2,…,n,其中d k =a kk -∑l kj 2k−1j=1d jl ik =(a ik −∑l ij l kj k−1j=1d j )/ d k 即d 1=a 11-∑l 1j 20j=1d j =16-0=16因为 l 21=(a 21−∑l 2j l 1j 0j=1d j )/ d 1=a 21/ d 1=416=14 所以d 2=a 22-∑l 2j 21j=1d j =5-(14)2d 1=4同理可得d 3=9 即得 D=(1649)同理l 11=(a 11−∑l ij l 1j 0j=1d j )/ d 1=1616=1=l 22=l 33 l 21=(a 21−∑l 2j l 1j 0j=1d j )/ d 1=416=14 l 31=(a 31−∑l 3j l 1j 0j=1d j )/ d 1=816=12 l 32=(a 32−∑l 3j l 2j 1j=1d j )/ d 2=−4−12×14×164=−64=-32即L=(114112−321) L T=(114121−321) 即LDL T分解为A=(114112−321)(1649)(114121−321)解解:A=(164845−48−422)→(41212−32−33)故得GG T分解:A=(4122−33)(4122−33) LDL T分解为A=(114112−321)(1649)(114121−321) 由(114112−321)(y 1y 2y 3)=(123) ,得(y 1y 2y 3)=(0.250.8751.7083)再由(4122−33)(x 1x 2x 3)=(0.250.8751.7083) ,得(x 1x 2x 3)=(−0.54511.29160.5694)2-4 用追赶法求解方程组:解:(4−1−14−1−14−1−14−1−14)→(4−14−1154−415−15615−1556−120956−56209−1780209)由(4−1154−15615−120956−1780209)(y1y2y3y4y5)=(100200),得(y1y2y3y4y5)=(256.66671.785700.4784753.718)再由(1−141−4151−15561−562091)(x1x2x3x4x5)=(256.66671.785700.4784753.718),得(x1x2x3x4x5)=(27.0518.20525.769314.87253.718)线性方程组解法(迭代法)2-1 设线性方程组{4x 1−x 2+2x 3=1−x 1−5x 2+x 3=22x 1+x 2+6x 3=3(1) 写出Jacobi 法和SOR 法的迭代格式(分量形式) (2) 讨论这两种迭代法的收敛性(3) 取初值x (0)=(0,0,0)T ,若用Jacobi 迭代法计算时,预估误差 ||x*-x (10)||∞ (取三位有效数字)解:(1)Jacobi 法和SOR 法的迭代格式分别为Jacobi 法迭代格式SOR(2)因为A 是严格对角占优矩阵,但不是正定矩阵,故Jacobi 法收敛,SOR 法当0<ω≤1时收敛.⎪⎪⎪⎩⎪⎪⎪⎨⎧+--=-+-=+-=+++216131525151412141)(2)(1)1(3)(3)(1)1(2)(3)(2)1(1k k k k k k k k k x x x x x x xx x ⎪⎪⎪⎩⎪⎪⎪⎨⎧-++-=+-+-=+-+-+=++++++)216131()525151()412141()(3)1(2)1(1)(3)1(3)(3)(2)1(1)(2)1(2)(3)(2)(1)(1)1(1k k k k k k k k k k k k k k k x x x x x x x x x x x x x x x ωωω(3)由(1)可见||B ||∞=3/4,且取x (0)=(0,0,0)T ,经计算可得x (1)=(1/4,-2/5,1/2)T ,于是||x (1)-x (0)||∞=1/2,所以有2-2 设方程组为{5x 1+2x 2+x 3=−12−x 1+4x 2+2x 3=202x 1−3x 2+10x 3=3试写出其Jacobi 分量迭代格式以及相应的迭代矩阵,并求解。

数值分析考试试题及答案

数值分析考试试题及答案

数值分析考试试题及答案一、单项选择题(每题2分,共10分)1. 以下哪个算法是用于求解线性方程组的?A. 欧几里得算法B. 高斯消元法C. 牛顿迭代法D. 傅里叶变换答案:B2. 插值法中,拉格朗日插值和牛顿插值的主要区别是什么?A. 计算复杂度不同B. 插值点的选取不同C. 插值多项式的表达形式不同D. 适用的函数类型不同答案:C3. 在数值分析中,以下哪个概念用于描述数值解的误差?A. 收敛性B. 稳定性C. 精度D. 条件数答案:C4. 以下哪个数值方法是用于求解非线性方程的?A. 雅可比迭代B. 高斯-塞德尔迭代C. 牛顿法D. 欧拉法答案:C5. 以下哪个数值积分方法是基于梯形规则的?A. 辛普森规则B. 梯形规则C. 龙格-库塔方法D. 蒙特卡洛方法答案:B二、填空题(每题2分,共10分)1. 在数值分析中,条件数是衡量一个数学问题_______的量度。

答案:解对输入数据变化的敏感性2. 线性方程组的迭代法中,如果迭代矩阵的谱半径小于1,则该迭代法是_______的。

答案:收敛3. 函数f(x)在区间[a, b]上的积分可以通过_______方法来近似计算。

答案:数值积分4. 矩阵A的特征值和特征向量可以通过_______方法求解。

答案:幂迭代5. 在求解微分方程的数值方法中,_______方法是一种常用的显式方法。

答案:欧拉三、计算题(每题10分,共20分)1. 给定线性方程组\(\begin{cases} x+y=5 \\ 2x-y=1\end{cases}\),使用高斯消元法求解x和y的值。

答案:首先写出增广矩阵:\[\begin{bmatrix}1 & 1 & | & 5 \\2 & -1 & | & 1\end{bmatrix}\]然后进行行变换,将第二行减去第一行的两倍:\[\begin{bmatrix}1 & 1 & | & 5 \\0 & -3 & | & -9\end{bmatrix}\]解得y=3,代入第一个方程得到x=2。

数值分析考试及答案

数值分析考试及答案

数值分析考试及答案一、单项选择题(每题3分,共15分)1. 以下哪个方法用于求解线性方程组的直接法?A. 高斯消元法B. 雅可比迭代法C. 牛顿迭代法D. 梯度下降法答案:A2. 在数值分析中,插值与逼近的主要区别是什么?A. 插值要求函数值完全匹配,逼近则不需要B. 逼近要求函数值完全匹配,插值则不需要C. 插值和逼近都要求函数值完全匹配D. 插值和逼近都不需要函数值完全匹配答案:A3. 以下哪个选项是数值积分的常用方法?A. 梯形法则B. 牛顿迭代法C. 欧拉方法D. 龙格-库塔方法答案:A4. 以下哪个方法用于求解非线性方程?A. 高斯消元法B. 牛顿迭代法C. 梯形法则D. 欧拉方法答案:B5. 在数值分析中,以下哪个方法用于求解线性最小二乘问题?A. 高斯消元法B. 梯度下降法C. 正交投影法D. 牛顿迭代法答案:C二、填空题(每题4分,共20分)6. 线性方程组的高斯消元法中,如果遇到主元为零,则需要进行________。

答案:行交换7. 在数值分析中,插值多项式的拉格朗日形式中,每个基函数的形式为________。

答案:\( L_i(x) = \prod_{\substack{0 \le j \le n \\ j \ne i}} \frac{x - x_j}{x_i - x_j} \)8. 梯形法则的误差与被积函数的________阶导数有关。

答案:二阶9. 牛顿迭代法求解非线性方程时,每次迭代的公式为________。

答案:\( x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)} \)10. 线性最小二乘问题的正交投影法中,最小二乘解是________的正交投影。

答案:观测向量三、简答题(每题10分,共30分)11. 简述数值稳定性的概念及其重要性。

答案:数值稳定性是指数值算法在计算过程中对输入数据的微小变化不会产生过大的输出变化。

数值稳定性的重要性在于,它保证了算法的可靠性和准确性,尤其是在处理大规模数据或者要求高精度的科学计算中,数值稳定性是算法设计和选择的重要考虑因素。

数值分析试题及答案

数值分析试题及答案

数值分析试题及答案一、选择题(每题5分,共20分)1. 以下哪个算法是用于求解非线性方程的?A. 高斯消元法B. 牛顿迭代法C. 雅可比方法D. 追赶法答案:B2. 线性代数中,矩阵的秩是指:A. 矩阵的行数B. 矩阵的列数C. 矩阵中线性无关的行(或列)的最大数目D. 矩阵对角线上元素的乘积答案:C3. 在数值分析中,插值法主要用于:A. 求解微分方程B. 求解积分方程C. 通过已知数据点估计未知数据点的值D. 优化问题答案:C4. 以下哪个方法不是数值积分的方法?A. 梯形规则B. 辛普森规则C. 高斯消元法D. 龙贝格积分答案:C二、填空题(每题5分,共20分)1. 线性方程组的解的存在性和唯一性可以通过矩阵的________来确定。

答案:秩2. 牛顿迭代法中,迭代公式为 x_{n+1} = x_n - f(x_n)/f'(x_n),其中 f'(x_n) 表示函数 f(x) 在 x_n 处的________。

答案:导数3. 在数值积分中,如果积分区间被划分得越细,那么数值积分的误差通常会________。

答案:减小4. 用泰勒级数展开函数 f(x) 时,如果 x0 是展开点,那么 f(x) 可以表示为 f(x) = f(x0) + f'(x0)(x-x0) + ________。

答案:更高阶项三、简答题(每题10分,共20分)1. 简述数值稳定性的概念及其在数值分析中的重要性。

答案:数值稳定性是指数值算法在计算过程中对初始数据微小扰动的敏感程度。

一个数值稳定的算法能够保证在有限精度的计算下,算法的结果不会因为数值误差而产生较大的偏差。

在数值分析中,数值稳定性是算法设计的重要考虑因素,它直接影响到算法的可靠性和准确性。

2. 描述一下在数值分析中,如何使用龙贝格积分法来提高数值积分的精度。

答案:龙贝格积分法是一种基于梯形规则的数值积分方法,它通过递归地增加梯形规则的子区间数量来提高积分的精度。

数值分析习题(含标准答案)

数值分析习题(含标准答案)

数值分析习题(含标准答案)
一、选择题(每题5分,共20分)
1. 下列哪个选项不属于数值分析的研究范畴?
A. 数值微分
B. 数值积分
C. 数值逼近
D. 数据库管理
答案:D
2. 在数值分析中,求解线性方程组常用的方法有?
A. 高斯消元法
B. 迭代法
C. 拉格朗日乘数法
D. 上述所有方法
答案:D
3. 下列哪种方法适用于求解非线性方程组?
A. 牛顿法
B. 梯度下降法
C. 高斯消元法
D. 上述所有方法
答案:D
4. 在数值积分中,下列哪种方法具有最高的精度?
A. 梯形法则
B. 辛普森法则
C. 高斯求积法
D. 上述所有方法
答案:C
二、填空题(每题5分,共20分)
1. 数值分析的主要目的是通过有限步骤的运算,对数学问题进行近似求解。

2. 在数值微分中,常用的差分公式有前向差分、后向差分和中心差分。

3. 数值逼近的主要方法包括插值法和逼近法。

4. 在数值积分中,常用的方法有梯形法则、辛普森法则和高斯求积法。

三、解答题(每题10分,共30分)
1. 已知函数 f(x) = e^x,求其在 x = 0.5 处的导数。

答案:f'(0.5) ≈ 1.6487
2. 求解线性方程组 2x + 3y = 5,4x y = 1。

答案:x ≈ 0.625,y ≈ 1.25
3. 已知函数 f(x) = x^3 3x^2 + 4,求其在区间 [0, 2] 上的积分。

答案:f(x) 在区间 [0, 2] 上的积分≈ 3.6667。

数值分析试题及答案

数值分析试题及答案

数值分析试题及答案一、选择题(每题3分,共30分)1. 下列关于数值分析的说法,错误的是()。

A. 数值分析是研究数值方法的科学B. 数值分析是研究数值方法的数学理论C. 数值分析是研究数值方法的误差分析D. 数值分析是研究数值方法的数学理论、误差分析及数值方法的实现答案:B2. 在数值分析中,插值法主要用于()。

A. 求解微分方程B. 求解积分方程C. 求解线性方程组D. 通过已知数据点构造一个多项式答案:D3. 线性方程组的解法中,高斯消元法属于()。

A. 直接方法B. 迭代方法C. 矩阵分解方法D. 特征值方法答案:A4. 牛顿法(Newton's method)是一种()。

A. 插值方法B. 拟合方法C. 迭代方法D. 优化方法答案:C5. 在数值分析中,下列哪种方法用于求解非线性方程的根?A. 高斯消元法B. 牛顿法C. 雅可比方法D. 斯托尔-温格尔方法答案:B6. 下列关于误差的说法,正确的是()。

A. 绝对误差总是大于相对误差B. 相对误差总是小于绝对误差C. 误差是不可避免的D. 误差总是可以消除的答案:C7. 在数值分析中,下列哪个概念与数值稳定性无关?A. 条件数B. 截断误差C. 舍入误差D. 插值多项式的阶数答案:D8. 用泰勒级数展开函数f(x)=e^x,下列哪一项是正确的?A. f(x) = 1 + x + x^2/2! + x^3/3! + ...B. f(x) = 1 - x + x^2/2! - x^3/3! + ...C. f(x) = x + x^2/2 + x^3/6 + ...D. f(x) = x - x^2/2 + x^3/6 - ...答案:A9. 插值多项式的次数最多为()。

A. n-1B. nC. n+1D. 2n答案:B10. 下列关于数值积分的说法,错误的是()。

A. 梯形法则是一种数值积分方法B. 辛普森法则是一种数值积分方法C. 龙格法则是数值积分方法中的一种D. 数值积分方法总是精确的答案:D二、填空题(每题3分,共15分)1. 在数值分析中,条件数是衡量问题的______。

数值分析试题及答案

数值分析试题及答案

数值分析试题及答案一、单项选择题(每题3分,共30分)1. 线性代数中,矩阵A的逆矩阵记作()。

A. A^TB. A^-1C. A^+D. A*答案:B2. 插值法中,拉格朗日插值多项式的基函数是()。

A. 多项式B. 指数函数C. 正弦函数D. 余弦函数答案:A3. 在数值积分中,梯形规则的误差是()阶的。

A. O(h^2)B. O(h^3)C. O(h)D. O(1/h)答案:A4. 求解线性方程组时,高斯消元法的基本操作不包括()。

A. 行交换B. 行乘以非零常数C. 行加行D. 行除以非零常数答案:D5. 非线性方程f(x)=0的根的迭代法中,收敛的必要条件是()。

A. f'(x)≠0B. f'(x)=0C. |f'(x)|<1D. |f'(x)|>1答案:C6. 利用牛顿法求解非线性方程的根时,需要计算()。

A. 函数值B. 函数值和导数值C. 函数值和二阶导数值D. 函数值、一阶导数值和二阶导数值答案:B7. 矩阵的特征值和特征向量是()问题中的重要概念。

A. 线性方程组B. 特征值问题C. 线性规划D. 非线性方程组答案:B8. 在数值分析中,条件数是衡量矩阵()的量。

A. 稳定性B. 可逆性C. 正交性D. 稀疏性答案:A9. 利用龙格现象说明,高阶插值多项式在区间端点附近可能产生()。

A. 振荡B. 收敛C. 稳定D. 单调答案:A10. 雅可比迭代法和高斯-塞德尔迭代法都是求解线性方程组的()方法。

A. 直接B. 迭代C. 精确D. 近似答案:B二、填空题(每题4分,共20分)11. 线性代数中,矩阵A的行列式记作________。

答案:det(A) 或 |A|12. 插值法中,牛顿插值多项式的基函数是________。

答案:差商13. 在数值积分中,辛普森规则的误差是________阶的。

答案:O(h^4)14. 求解线性方程组时,迭代法的基本思想是从一个初始近似解出发,通过不断________来逼近精确解。

数值分析试题及答案汇总

数值分析试题及答案汇总

数值分析试题及答案汇总一、单项选择题(每题5分,共20分)1. 在数值分析中,下列哪个方法用于求解线性方程组?A. 牛顿法B. 插值法C. 迭代法D. 泰勒展开法答案:C2. 以下哪个选项是数值分析中用于求解非线性方程的迭代方法?A. 高斯消元法B. 牛顿法C. 多项式插值D. 辛普森积分法答案:B3. 以下哪个选项是数值分析中用于数值积分的方法?A. 牛顿法B. 辛普森积分法C. 牛顿-拉弗森迭代D. 拉格朗日插值答案:B4. 在数值分析中,下列哪个方法用于求解常微分方程的初值问题?A. 欧拉法B. 牛顿法C. 辛普森积分法D. 高斯消元法答案:A二、填空题(每题5分,共20分)1. 插值法中,拉格朗日插值法的插值多项式的阶数是______。

答案:n2. 泰勒展开法中,如果将函数展开到第三阶,那么得到的多项式是______阶多项式。

答案:三3. 在数值分析中,牛顿法求解非线性方程的迭代公式为______。

答案:x_{n+1} = x_n - f(x_n) / f'(x_n)4. 辛普森积分法是将积分区间分为______等分进行近似计算。

答案:偶数三、简答题(每题10分,共30分)1. 请简述数值分析中插值法的基本原理。

答案:插值法的基本原理是根据一组已知的数据点,构造一个多项式函数,使得该函数在给定的数据点上与数据值相等,以此来估计未知数据点的值。

2. 解释数值分析中误差的概念,并说明它们是如何影响数值计算结果的。

答案:数值分析中的误差是指由于计算方法或计算工具的限制,导致计算结果与真实值之间的差异。

误差可以分为舍入误差和截断误差。

舍入误差是由于计算机表示数值的限制而产生的,而截断误差是由于计算方法的近似性质而产生的。

这些误差会影响数值计算结果的准确性和稳定性。

3. 请说明在数值分析中,为什么需要使用迭代法求解线性方程组。

答案:在数值分析中,迭代法用于求解线性方程组是因为对于大规模的方程组,直接方法(如高斯消元法)的计算成本很高,而迭代法可以在较少的计算步骤内得到近似解,并且对于稀疏矩阵特别有效。

数值分析期末考试复习题及其答案

数值分析期末考试复习题及其答案

数值分析期末考试复习题及其答案1. 已知325413.0,325413*2*1==X X 都有6位有效数字,求绝对误差限.(4分)解:由已知可知,n=65.01021,0,6,10325413.0016*1=⨯==-=⨯=ε绝对误差限n k k X 2分 620*21021,6,0,10325413.0-⨯=-=-=⨯=ε绝对误差限n k k X 2分2. 已知⎢⎢⎢⎣⎡=001A 220- ⎥⎥⎥⎦⎤440求21,,A A A ∞ (6分)解:{},88,4,1max 1==A 1分 {},66,6,1max ==∞A 1分 ()A A A T max 2λ= 1分⎢⎢⎢⎣⎡=001A A T 420 ⎥⎥⎥⎦⎤-420⎢⎢⎢⎣⎡001 220- ⎥⎥⎥⎦⎤440=⎢⎢⎢⎣⎡001 080 ⎥⎥⎥⎦⎤3200 2分 {}3232,8,1max )(max ==A A T λ 1分 24322==A3. 设32)()(a x x f -= (6分) ① 写出f(x )=0解的Newton 迭代格式② 当a 为何值时,)(1k k x x ϕ=+ (k=0,1……)产生的序列{}k x 收敛于2解:①Newton 迭代格式为:xa x x x ax a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(22321+=+=---=-=+ϕ 3分②时迭代收敛即当222,11210)2(',665)('2<<-<-=-=a a x a x ϕϕ 3分4. 给定线性方程组Ax=b ,其中:⎢⎣⎡=13A ⎥⎦⎤22,⎥⎦⎤⎢⎣⎡-=13b 用迭代公式)()()()1(k k k Ax b x x -+=+α(k=0,1……)求解Ax=b ,问取什么实数α,可使迭代收敛 (8分)解:所给迭代公式的迭代矩阵为⎥⎦⎤--⎢⎣⎡--=-=ααααα21231A I B 2分其特征方程为0)21(2)31(=----=-αλαααλλB I 2分即,解得αλαλ41,121-=-= 2分 要使其满足题意,须使1)(<B ρ,当且仅当5.00<<α 2分5. 设方程Ax=b,其中⎢⎢⎢⎣⎡=211A 212 ⎥⎥⎥⎦⎤-112,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=765b 试讨论解此方程的Jacobi 迭代法的收敛性,并建立Gauss —Seidel 迭代格式 (9分)解:U D L A ++=⎢⎢⎢⎣⎡--=+-=-210)(1U L D B J 202-- ⎥⎥⎥⎦⎤-012 3分0,03213=====-λλλλλJ B I 2分即10)(<=J B ρ,由此可知Jacobi 迭代收敛 1分 Gauss-Seidel 迭代格式:⎪⎩⎪⎨⎧--=--=+-=++++++)1(2)1(1)1(3)(3)1(1)1(2)(3)(2)1(12276225k k k k k k k k k x x x x x x x x x (k=0,1,2,3……) 3分6. 用Doolittle 分解计算下列3个线性代数方程组:i i b Ax =(i=1,2,3)其中⎢⎢⎢⎣⎡=222A 331 ⎥⎥⎥⎦⎤421,23121,,974x b x b b ==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡= (12分)解:①11b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=9741x A=⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211=LU 3分 由Ly=b1,即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡974 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡234 1分 由Ux1=y ,即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡234 得x1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 2分 ②22b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 由Ly=b2=x1,即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001 1分 由Ux2=y,即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001 得x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 2分③33b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0由Ly=b3=x2,即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-05.05.0 1分 由Ux3=y ,即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-05.05.0 得x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-025.0375.0 2分7. 已知函数y=f (x)有关数据如下:要求一次数不超过3的H 插值多项式,使'11'33)(,)(y x H y x H i i == (6分)解:作重点的差分表,如下:3分21021101011001003))(](,,,[))(](,,[)](,[][)(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 H --+--+-+= =-1+(x+1)-x (x+1)+2x.x(x+1)=232x x + 3分8. 有如下函数表:试计算此列表函数的差分表,并利用Newton 前插公式给出它的插值多项式 (7分)解:由已知条件可作差分表,3分i ih x x i =+=0 (i=0,1,2,3)为等距插值节点,则Newton 向前插值公式为: 033210022100003!3))()((!2))((!1)()(f h x x x x x x f h x x x x f h x x f x N ∆---+∆--+∆-+==4+5x+x (x-1)=442++x x 4分9. 求f (x )=x 在[-1,1]上的二次最佳平方逼近多项式)(2x P ,并求出平方误差 (8分)解:令22102)(x a x a a x P ++= 2分取m=1, n=x , k=2x ,计算得: (m ,m)=dx ⎰-111=0 (m,n )=dx x ⎰-11=1 (m,k)=dx x ⎰-112=0(n,k )=dx x ⎰-113=0。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第一章1、ln2=0.69314718…,精确到 10-3的近似值是多少?解 精确到 10-3=0.001,即绝对误差限是 ε=0.05%,故至少要保留小数点后三位才可以。

ln2≈0.693。

2、设115.80,1025.621≈≈x x 均具有5位有效数字,试估计由这些数据计算21x x ,21x x +的绝对误差限解:记126.1025, 80.115x x ==ﻩ则有1123241110, | 102|||2x x x x --≤⨯-≤⨯-所以 121212121212211122||||||||||||x x x x x x x x x x x x x x x x x x -=-+-+≤--341180.11610 6.101025220.007057-==⨯⨯+≤⨯⨯ﻩ1212112243|()|||11|10100.0005522|x x x x x x x x --≤≤⨯+⨯=+-+-+-3、一个园柱体的工件,直径d 为10.25±0.25mm,高h为40.00±1.00mm ,则它的体积V 的近似值、误差和相对误差为多少。

解:()()22222222431421025400000033006422102540000251025100243644433006243624360073873833006,.....;()()()......,..().()..%.r d hV d h V mm d h V dh d d h V mm V V V πππππεεεεε=≈=⨯⨯===+=⨯⨯⨯+⨯==±====第二章:1、分别利用下面四个点的Lag ra ng e插值多项式和Newton 插值多项式N 3(x ),计算L 3(0.5)及N3(-0.5)解:(1)先求332211003)()()()()(y x l y x l y x l y x l x L +++=(1分)=----+---+=------=)12)(02)(12()1)(0)(1())()(())()(()(3020103210x x x x x x x x x x x x x x x x l x x x )1)(1(61-+-, (2分)=----+---+=------=)11)(01)(21()1)(0)(2())()(())()(()(3121013201x x x x x x x x x x x x x x x x l x x x )1)(2(21-+ ﻩ(2分)=-++-++=------=)10)(10)(20()1)(1)(2())()(())()(()(3212023102x x x x x x x x x x x x x x x x l )1)(1)(2(21-++-x x x ﻩ(2分)=-++-++=------=)01)(11)(21()0)(1)(2())()(())()(()(2313032103x x x x x x x x x x x x x x x x l x x x )1)(2(61++ ﻩ(2分)x x x x x x x x x x L )1)(2(31)1)(2(21)1)(1(61)(3+++-++-+=x x x 212323-+=ﻩ(1分)所以 41)5.0(3=L ﻩ ﻩﻩ(1分)(2)再求Ne wton 插值多项式 列均差表如下:)(123221)(23100)(211)(12],,,[],,[],[222232103210分分分分x x x x x x x x f x x x f x x f y x k j i j i -----所以x x x x x x x N )1)(2()1)(2(23)2(21)(3+++++-++-=x x x 212323-+= (2分) 21)5.0(3=-N ﻩﻩ (1分)2、求过下面四个点的L agrang e插值多项式L 3(x )和Newto n插值多项式N 3(x)。

)解:(1)L3(x )=l o(x )yo +l 1(x )y1+l 2(x)y 2+l 3(x )y 3(1分))())(())(()())(()1)(()(1110110n i i i i i i i n i i i x x x x x x x x x x x x x x x x x x x x l ---------+=+-+-得出)1)(1(61)(-+-=x x x x l o(2分))1)(2(21)(1-+=x x x x l(2分))1)(1)(2(21)(2-++-=x x x x x l (2分))1)(2(61)(3++=x x x x l (2分)∴)1)(2(61)1)(1)(2(21)1)(2(21)1)(1(31)(3++--++--++-+=x x x x x x x x x x x x x L(1分)(2)))()(())(()()(21031020103x x x x x x a x x x x a x x a a x N ---+--+-+=(1分)2)(00-==x f a (2分) 3)()(10101=--=x x x f x f a ﻩ (2分)23)()()()(20212110102-=------=x x x x x f x f x x x f x f a (2分),613=a ﻩ (2分)∴x x x x x x x N )1)(2(61)1)(2(23)2(32)(3+++++-++-=(1分)第三章 1、令1x 1,e )x (f x≤≤-=,且设x a a )x (p 1+=,求1a,a 使得)x (p 为)x (f 在[-1,1]上的最佳平方逼近多项式。

2.已知数据对(7,3.1),(8,4.9),(9,5.3),(10,5.8),(11,6.1), (12,6.4),(13,5.9)。

试用二次多项式拟合这组数据。

解:y=-0.145x2+3.324x -12.794第四章:1.数据如下表用中心差分公式,分别取h = 0.01、0.02计算)02.1(f '.解:中心差分公式为 hh x f h x f x f 2)()()(--+≈' (2分)1)取h =0.01时, 302.012.318.302.0)01.1()03.1()02.1(=-=-≈'f f f ﻩ (4分)2)取h =0.02时, 5.304.010.324.304.0)00.1()04.1()02.1(=-=-≈'ff f (4分)2.(10分)根据如下函数表 用中心差分公式,分别取h =0.3,0.1计算)3.1(f '解:中心差分公式hh x f h x f x f 2)()()(--+='ﻩﻩ (2分)取h =0.3时,7233.16.0)3.03.1()3.03.1()(=--+≈'f f x f ﻩ (4分)取h =0.1时,7000.12.0)1.03.1()1.03.1()(=--+≈'f f x f ﻩ ﻩ(4分)3.分别用复合梯形公式T6和复合辛普森公式S 3计算定积分⎰+6.00d 11x x的值. 解:)2)()0((2116∑-=++-=n i i n y x f f n ab T (2分) )])5.0()4.0()3.0()2.0()1.0([2)6.0()0((6206.0f f f f f f f ++++++⨯-=470510739.0=ﻩﻩﻩﻩﻩﻩ(3分))]}]4.0()2.0([2)]5.0()3.0()1.0([4)6.0()0({63f f f f f f f nab S ++++++-=470006382.0= ﻩﻩﻩﻩﻩﻩﻩ(3分)f (0)=1,f (0.1)=0.9090,f (0.2)=.08333,f (0.3)=0.7692,f (0.4)=0.7142,f (0.5)=0.6667,f(0.6)=0.625ﻩﻩﻩ(7分)4、利用复合Simpso n公式S4计算积分⎰+102d 11x x (取小数点后4位)。

解:)]24()2()0([611212∑∑==-+++⨯-=n i ni i i n y y n f f n ab S (2分) 00000.4)0(=f ,93846.381=⎪⎭⎫ ⎝⎛f ,76470.382=⎪⎭⎫⎝⎛f ,50685.383=⎪⎭⎫ ⎝⎛f ,20000.384=⎪⎭⎫ ⎝⎛f ,87640.285=⎪⎭⎫ ⎝⎛f ,56000.286=⎪⎭⎫⎝⎛f ,26000.287=⎪⎭⎫ ⎝⎛f ,00000.2)1(=f ﻩﻩ(9分) ⎭⎬⎫⎩⎨⎧⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛++⨯-=8)684822878583814)1()0(46014f f f f f f f f f S1416.3=ﻩﻩ(4分)第五章:1、利用列主元消去法求解线性方程组⎪⎩⎪⎨⎧=+-=-=++-6557710462332121321x x x x x x x x (计算过程保留到小数点后四位).解:216515707104623r r ↔⎪⎪⎪⎭⎫ ⎝⎛---(1分)⎪⎪⎪⎭⎫⎝⎛---6515462370710(2分) ⎪⎪⎪⎭⎫⎝⎛--+-5.255.201.661.0070710103132112r r r r (2分)⎪⎪⎪⎭⎫ ⎝⎛-↔+2.62.6005.255.2070710235.21.032r r r r (2分)回代解得 13=x ,12-=x ,03=x ﻩﻩ (1分)2、用矩阵的L U分解法解方程组⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛---432214122012321x x x解:设⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫⎝⎛==332221131211323121000101001u u u u u u l ll LU A ﻩﻩﻩ(1分) ⎪⎪⎪⎭⎫⎝⎛-⎪⎪⎪⎭⎫ ⎝⎛-=100110012112011001LU ﻩﻩ(4分)LUX =b其中设UX =y,则Ly =b ⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛-432112*********y y y(2分)∴y =(2,-1,1)T U X=y ⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛-112100110012321x x x ﻩﻩ(2分)∴x =(0,-2,1)T ﻩﻩ(1分)5. 用追赶法解三对角方程组Ax=b ,其中解:用解对三角方程组的追赶法公式计算得ﻫ6. 用平方根法解方程组解:用分解直接算得ﻫﻫ由及求得第六章:1、用G auss-Sei del 迭代法求解方程组⎪⎩⎪⎨⎧=+-=++=++301532128243220321321321x x x x x x x x x ,取初值T)0()0,0,0(=x ,写出G aus s-Se ide l迭代格式,求出)1(x ,)2(x,计算∞-)2()1(x x,并根据原方程组的系数矩阵说明该迭代格式是否收敛.2、对方程组⎪⎩⎪⎨⎧=+--=-+-=--1052151023210321321321x x x x x x x x x (1)写出其Jaco bi 迭代格式,并据迭代矩阵的范数,说明该迭代格式收敛。

相关文档
最新文档