武汉大学数值分析期末考试题目和答案

合集下载

数值分析试题与答案

一. 填空题（本大题共4小题，每小题4分，共16分）1.设有节点012,,x x x ，其对应的函数()y f x =的值分别为012,,y y y ，则二次拉格朗日插值基函数0()l x 为。

2.设()2f x x =，则()f x 关于节点0120,1,3x x x ===的二阶向前差分为。

3.设110111011A -⎡⎤⎢⎥=--⎢⎥⎢⎥-⎣⎦，233x ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，则1A ＝，1x = 。

4. 1n +个节点的高斯求积公式的代数精确度为。

二．简答题（本大题共3小题，每小题8分，共24分）1. 哪种线性方程组可用平方根法求解？为什么说平方根法计算稳定？2. 什么是不动点迭代法？()x ϕ满足什么条件才能保证不动点存在和不动点迭代序列收敛于()x ϕ的不动点？3. 设n 阶矩阵A 具有n 个特征值且满足123n λλλλ>≥≥≥，请简单说明求解矩阵A 的主特征值和特征向量的算法及流程。

三．求一个次数不高于3的多项式()3P x ，满足下列插值条件：i x 1 2 3 i y 2 4 12 i y '3并估计误差。

（10分）四．试用1,2,4n =的牛顿-科特斯求积公式计算定积分1011I dx x=+⎰。

（10分）五．用Newton 法求()cos 0f x x x =-=的近似解。

（10分）六．试用Doolittle 分解法求解方程组：12325610413191963630x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦⎣⎦ （10分）七．请写出雅可比迭代法求解线性方程组123123123202324812231530x x x x x x x x x ++=⎧⎪++=⎨⎪-+=⎩ 的迭代格式，并判断其是否收敛？（10分）八．就初值问题0(0)y yy y λ'=⎧⎨=⎩考察欧拉显式格式的收敛性。

（10分）《数值分析》（A ）卷标准答案（2009－2010－1）一．填空题（每小题3分，共12分） 1. ()1200102()()()()x x x x l x x x x x --=--； 2.7；3. 3，8；4. 2n+1。

数值分析期末试题及答案

数值分析期末试题及答案一、选择题（每题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)。

数值分析期末复习题答案

数值分析期末复习题答案一、选择题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 是一个列向量。

数值分析试题及答案解析

数值分析试题及答案解析数值分析试题一、填空题（2 0×2′）1.-=?-=32,1223X A 设x =0.231是精确值x *=0.229的近似值，则x 有 2位有效数字。

2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ，f [20,21,22,23,24,25,26,27,28]= 0 。

3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____，‖AX ‖∞≤_15_ __。

4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足|?’(x )| <1 ，则使用该迭代函数的迭代解法一定是局部收敛的。

5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。

6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商公式的前插公式，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的后插公式；如果要估计结果的舍入误差，应该选用插值公式中的拉格朗日插值公式。

7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=ni i x a 0)( 1 ；所以当系数a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。

8. 要使20的近似值的相对误差小于0.1%，至少要取 4 位有效数字。

9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收敛于方程组的精确解x *的充分必要条件是ρ(B)<1 。

10. 由下列数据所确定的插值多项式的次数最高是 5 。

11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。

12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的，其中的残差r i ＝ (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ，(i =0,1,…,n )。

数值分析复习题参考答案

x1 )
h
2
x 0 x x1 6
4
所以， R x
h 10
2
8
解得， h 0 . 000383
4. 习题(第二章) 7
5. 习题(第二章) 9
6. 习题(第二章) 11
7. 习题(第二章) 13
8. 习题(第二章) 14
9. 习题(第二章) 20
10. 习题(第四章) 1
2
, k 0 ,1, 2 2 3 2a 3x
3
此时， ( x )
2x a 3x
, '( x) 2a

所以， ' ( 3 a )
2 3
3(
3
a)
3
0 1, 所以该迭代公式收敛。
21. 习题(第七章) 13
本题没有给出精度要求，但x3与x2之间的差为已经很小了，足以满足精度。
[ f ( x n , y n ) f ( x n 1 , y n 1 )]
( 3 ) 基于 Taylor 展开法：
y ( x n 1 ) y ( x n h ) y ( x n ) y ' ( x n ) h
h
2
2
y ''( xn )
取 y ( x n 1 ) y ( x n ) y ' ( x n ) h ，即 y n 1 y n hf ( x n , y n )
k 个点的值
求解隐式：先用欧拉公求解多步法：单步法开
式求得一个初步的近似表头，然后预报
修正校正修正。
( 其实只要给出公式会用
就行！！ )

数值分析习题（含答案）

数值分析习题（含答案）第一章绪论姓名学号班级习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。

1 若误差限为5105.0-?,那么近似数0.003400有几位有效数字？（有效数字的计算）解：2*103400.0-?=x ，325*10211021---?=?≤-x x 故具有3位有效数字。

2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算）解：10314159.0?= π，欲使其近似值*π具有4位有效数字，必需41*1021-?≤-ππ，3*310211021--?+≤≤?-πππ，即14209.314109.3*≤≤π即取（3.14109 , 3.14209）之间的任意数，都具有4位有效数字。

3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ?有几位有效数字？（有效数字的计算）解：3*1021-?≤-aa ，2*1021-?≤-b b ，而1811.2=+b a ，1766.1=?b a 2123****102110211021)()(---?≤?+?≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。

2123*****10210065.01022031.1102978.0)()(---?≤=?+?≤-+-≤-b b a a a b ba ab 故b a ?至少具有2位有效数字。

4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算）解：已知δ=-**xx x ，则误差为δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=，已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v2π=的绝对误差限与相对误差限。

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

数值分析期末考试题及答案一、选择题（每题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 \) 处插值，并求出插值多项式。

武汉大学数值分析期末考试(05-11年)

x1 + 2 x 2 − 2 x3 = 1, x1 + x 2 + x3 = −1, 2 x + 2 x + x = 0. 2 3 1 2a a 0 二、（8 分）若矩阵 A = 0 a 0 ，说明对任意实数 a ≠ 0 ，方程组 AX = b 都是 0 0 a
dy = f (t , y ) 的单步法：八、（14 分）对于下面求解常微分方程初值问题 dt y (t 0 ) = y 0
λ
y n +1 = y n + hk 2 k1 = f (t n , y n ) 1 1 k 2 = f (t n + h, y n + hk1 ) 2 2 （1）验证它是二阶方法；（2）确定此单步法的绝对稳定区域。
1 1 y n +1 = y n + h( 2 k1 + 2 k 2 ) k1 = f ( x n , y n ) k = f ( x + h, y + hk ) n n 1 2 （1）验证它是二阶方法；（2）确定此单步法的绝对稳定域。
武汉大学
2006~2007 学年第一学期硕士研究生期末考试试题（A 卷）
非病态的（范数用 ⋅ ∞ ）四、（15 分）已知 y = f ( x) 的数据如下： xi f ( xi ) f ′( xi )
0
0
1 2 1
2
6
求 f ( x) 的 Hermite 插值多项式 H 3 ( x) ，并给出截断误差 R ( x) = f ( x) − H 3 ( x) 。五、（10 分）在某个低温过程中，函数 y 依赖于温度 x（℃）的试验数据为 xi yi 1 0．8 2 1．5 3 1．8 4 2．0

武汉大学数值分析期末考试题目和答案

（2 分）
（2 分）（1 分）
注：1、教师命题时题目之间不留空白； 2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考生须在试题图上作解答，请另附该试题图。3、请在试卷类型、考试方式后打勾注明。（第 5 页）
八．证明题（本大题共 2 小题，每小题 7 分，共 14 分） 1. 证：该问题的精确解为 y( x) y0e
六．试用 Doolittle 分解法求解方程组：
5 6 x1 1 0 2 4 1 3 1 9 x 1 9（ 10 分） 2 6 3 6 x3 3 0 20 x1 2 x2 3x3 24 七．请写出雅可比迭代法求解线性方程组 x1 8 x2 x3 12 的迭代格式，并 2 x 3x 15 x 30 2 3 1
故有 B 1.25 1 ，因而雅可比迭代法不收敛。（2）对于方程组，Gauss-Seidel 迭代法迭代矩阵为
0 0.5 0.5 B 0 0.5 0.5 0 0.5 0
其特征值为 1 0, 2 3 0.5 故有 B 0.5 1 ，因而雅可比迭代法收敛。
判断其是否收敛？（10 分）
y y 八．就初值问题考察欧拉显式格式的收敛性。（10 分） y (0) y0
注：1、教师命题时题目之间不留空白； 2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考生须在试题图上作解答，请另附该试题图。3、请在试卷类型、考试方式后打勾注明。（第 2 页）
步 6：若 k<N,置 k:=k+1, μ:=mk，转 3；否则输出计算失败信息，停止三．解：（1）利用插值法加待定系数法：设 p2 x 满足 p2 1 2, p2 2 4, p2 3 12, 则 p2 x 3x 7 x 6, （3 分）

武汉大学研究生课程数值分析期末考试

�� 5.2 Lagrange 插值多项式：��n (��) = ∑�� =0 �� (��)�� = ∑��=0 [{∏��=0 �� −�� } �� ]; ��+1 则��n (��) = ∑�� =0 (��−�� )��′ ��
武汉大学研究生课程《数值分析》半开卷考试资料
姓名：学号：第 1 章绪论 1.1 误差的基本概念绝对误差：∆ �� = �� ∗ − �� ; ∆f(x) = ��(��) = �� ′ (��)�� ; 绝对误差：∆�� = 有效数字：�� = (0. ��1 ��2 … �� × 10−�� ) × 10�� ,则有 n 位有效数字。 1 1 误差限：|∆ ��| = |�� ∗ − ��| ≤ × 10��−�� ; |∆�� | ≤ × 10−(��−1)
|�� ∗ −��+1 | ��→∞ |�� ∗ −��|��
= C (对于收敛的迭代格式，当|��′(��)| = 0，则是线性收敛)
��(�� )
��
若碰到求收敛阶：迭代公式是个方程，准确解带进去是个方程，两方程相减，然后适当变形利用微分中值定理。 3.3 Newton 法 (二阶收敛)：��+1 = �� − ��′(�� ) ; 假设�� ∗ 是 f(x) = 0 的单根, f(x)在�� ∗ 的邻域内具有连续的二阶导数且 f ′(�� ∗ ) ≠ 0 , 则牛顿公式具有局部收敛性；若 f′′(�� ∗ ) ≠ 0 且��0 ≠ �� ∗ , 则序列{�� }是平方收敛。第 4 章矩阵特征值特征向量 (略) Householder 变换（H=I-2wwT）、 Givens 变换、幂法第 5 章插值与逼近 5.1 插值多项式的唯一性 Pn (x) = ��0 + ��1 �� + ��2 �� 2 + ⋯ + �� ; |��| = ∏(�� − �� ) ≠ 0

数值分析题库

一．单项选择题（每小题2分，共10分）1．在下列四个数中，有一个数具有4位有效数字，且其绝对误差限为 51021-⨯，则该数是（） A 0.001523 B 0.15230 C 0.01523 D 1.52300 2．设方阵A 可逆，且其n 个特征值满足：n λλλ>≥> (21)，则1-A 的主特征值是（）A11λ B nλ1 C1λ或n λ D 11λ或nλ13．设有迭代公式→→+→+=fxB x k k )()1(。

若||B|| > 1，则该迭代公式（）A 必收敛B 必发散C 可能收敛也可能发散4．常微分方程的数值方法，求出的结果是（）A 解函数B 近似解函数C 解函数值D 近似解函数值 5．反幂法中构造向量序列时，要用到解线性方程组的（） A 追赶法 B LU 分解法C 雅可比迭代法D 高斯—塞德尔迭代法二．填空题（每小题4分，共20分）1．设有方程组⎪⎩⎪⎨⎧=+-=+-=+02132432132132x x x x x x x x ，则可构造高斯—塞德尔迭代公式为⎪⎩⎪⎨⎧2．设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=111112101A ，则=∞A3．设1)0(,2'2=+=y y x y ，则相应的显尤拉公式为=+1n y4．设1)(+=ax x f ,2)(x x g =。

若要使)(x f 与)(x g 在[0，1]上正交，则a =5．设T x )1,2,2(--=→，若有平面旋转阵P ，使P →x 的第3个分量为0，则P =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡ 三．计算题（每小题10分，共50分）1．求27的近似值。

若要求相对误差小于0.1%，问近似值应取几位有效数字？2．设42)(x x x f -=,若在[-1，0]上构造其二次最佳均方逼近多项式,请写出相应的法方程。

3．设有方程组⎪⎩⎪⎨⎧=++=++=-+1221122321321321x x x x x x x x x ，考察用雅可比迭代解此方程组的收敛性。

数值分析试卷及答案

数值分析试卷及答案**注意：以下是一份数值分析试卷及答案，试卷和答案分别按照题目和解答的格式排版，以确保整洁美观，语句通顺。

**---数值分析试卷一、选择题（每题2分，共20分）1. 数值分析是研究如何用计算机处理数值计算问题的一门学科。

以下哪个选项不是数值分析的应用领域？A. 金融风险评估B. 天气预测C. 数据挖掘D. 图像处理2. 在数值计算中，稳定性是指算法对于输入数据的微小扰动具有较好的性质。

以下哪个算法是稳定的？A. 高斯消元法B. 牛顿迭代法C. 不动点迭代法D. 雅可比迭代法二、填空题（每题3分，共30分）1. 下面关于插值多项式的说法中，不正确的是：一般情况下，插值多项式的次数等于插值点的个数减1。

2. 线性方程组中，如果系数矩阵A是奇异的，则该方程组可能无解或有无穷多解。

......三、解答题（共50分）1. 请给出用割线法求解非线性方程 f(x) = 0 的迭代格式，并选择合适的初始值进行计算。

解：割线法的迭代公式为：x_(k+1) = x_k - f(x_k) * (x_k - x_(k-1)) / (f(x_k) - f(x_(k-1)))选择初始值 x0 = 1，x1 = 2 进行计算：迭代1次得到：x2 = x1 - f(x1) * (x1 - x0) / (f(x1) - f(x0))迭代2次得到：x3 = x2 - f(x2) * (x2 - x1) / (f(x2) - f(x1))继续迭代直至满足精度要求。

2. 对于一个给定的线性方程组，高斯消元法可以用来求解其解空间中的向量。

请简要描述高斯消元法的基本思想并给出求解步骤。

高斯消元法的基本思想是通过一系列的行变换将线性方程组化为上三角形式，然后再通过回代求解方程组的未知数。

求解步骤如下：步骤1：将方程组表示为增广矩阵形式，即将系数矩阵和常数向量连接在一起。

步骤2：从第一行开始，选取第一个非零元素作为主元，然后通过行变换将其它行的该列元素消去。

武汉大学06-10年(缺08-09)研究生数值分析考试试卷

武汉大学2006~2007学年第一学期硕士研究生期末考试试题（A 卷）科目名称：数值分析学生所在院：学号：姓名：注意：所有的答题内容必须答在答题纸上，凡答在试题或草稿纸上的一律无效。

一、（12分）设方程组b Ax =为⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛37111221x x （1）用Doolittle 分解法求解方程组；（2）求矩阵A 的条件数∞)(A Cond二、（12分）设A 为n 阶对称正定矩阵，A 的n 个特征值为n λλλ≤≤≤ 21，为求解方程组b Ax =，建立迭代格式 )()()()1(k k k Ax b x x -+=+ω，求出常数ω的取值范围，使迭代格式收敛。

三、（12分）已知数据试用二次多项式c bx ax x p ++=2)(拟合这些数据。

四、（14分）已知 )(x f y = 的数据如下：（1）求)(x f 的Hermite 插值多项式)(3x H ；（2）为求⎰31)(dx x f 的值，采用算法：R dx x H dx x f +=⎰⎰31331)()(试导出截断误差R五、（12分）确定常数 a ，b 的值，使积分dx e b ax b a I x 210)(),(⎰-+=取得最小值。

六、（12）确定常数i A ，使求积公式)2()1()0()(32120f A f A f A dx x f ++≈⎰的代数精度尽可能高，并问是否是Gauss 型公式。

七、（12分）设)(x ϕ导数连续，迭代格式)(1k k x x ϕ=+一阶局部收敛到点*x 。

对于常数λ，构造新的迭代格式：)(1111k k k x x x ϕλλλ+++=+问如何选取λ，使新迭代格式有更高的收敛阶，并问是几阶收敛。

八、（14分）对于下面求解常微分方程初值问题 ⎪⎩⎪⎨⎧==00)(),(y t y y t f dt dy的单步法：⎪⎪⎩⎪⎪⎨⎧++==+=+)21,21(),(12121hk y h t f k y t f k hk y y n n n n n n （1）验证它是二阶方法；（2）确定此单步法的绝对稳定区域。

数值分析报告期末考试复习题及其问题详解

数值分析期末考试复习题及其答案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.5 (k,k)= dx x ⎰-114=0 (m,y)= dx x ⎰-11=1(n,y)=dx x⎰-112=0 (k,y)= dx x ⎰-113=0.5得方程组：⎪⎩⎪⎨⎧==+=5.05.005.011201a a a a 3分解之得c a a c a 2,1,210-=== （c 为任意实数，且不为零）即二次最佳平方逼近多项式222)(cx x c x P -+= 1分平方误差：32),(22222222=-=-=∑=i i i y a fp f ϕδ 2分10. 已知如下数据：用复合梯形公式，复合Simpson 公式计算⎰+=10214dx x π的近似值（保留小数点后三位）（8分）解：用复合梯形公式：)}1()]87()43()85()21()83()41()81([2)0({1618f f f f f f f f f T ++++++++==3.139 4分用复合Simpson 公式： )}1()]43()21()41([2)]87()85()83()81([4)0({2414f f f f f f f f f S ++++++++==3.142 4分11. 计算积分⎰=20sin πxdx I ，若用复合Simpson 公式要使误差不超过51021-⨯，问区间]2,0[π要分为多少等分?若改用复合梯形公式达到同样精确度，区间]2,0[π应分为多少等分? （10分）解： ①由Simpson 公式余项及x x f x x f sin )(,sin )()4(==得544)4(2041021)1()4(360)(max )4(1802)(-≤≤⨯≤=≤n x f n f R x n πππππ 2分即08.5,6654≥≥n n ，取n=6 2分即区间]2,0[π分为12等分可使误差不超过51021-⨯ 1分②对梯形公式同样1)(''max 20≤≤≤x f x π，由余项公式得51021)2(122)(-⨯≤≤n f R n ππ2分即255,2.254=≥n n 取 2分即区间]2,0[π分为510等分可使误差不超过51021-⨯ 1分12. 用改进Euler 格式求解初值问题：⎩⎨⎧==++1)1(0sin 2'y x y y y 要求取步长h 为0.1，计算y（1.1）的近似值（保留小数点后三位）[提示：sin1=0.84,sin1.1=0.89] （6分）解：改进Euler 格式为：⎪⎩⎪⎨⎧++=+=+-++-+)],(),([2),(1111n n n n n n n n n n y x f y x f hy y y x hf y y 2分于是有⎪⎩⎪⎨⎧+++-=+-=+-++-+-+)sin sin (05.0)sin (1.012112121n n n n n n n n n n n n n x y y x y y y y x y y y y （n=0,1,2……） 2分由y(1)=0y =1，计算得⎪⎩⎪⎨⎧=≈=+-=-838.0)1.1(816.0)1sin 11(1.01121y y y 2分即y(1.1)的近似值为0.83813. ][],[],,[lim ],[),,(],,[)(0'000000'x f x x f x x f x x f b a x b a C x f x x ==∈∈→证明：定义：设（4分）证明：]['],[],[],[lim ][][lim]['00000000000x f x x f x x f x x f x x x f x f x f x x x x ===--=→→故可证出 4分14. 证明：设nn RA ⨯∈，⋅为任意矩阵范数，则A A ≤)(ρ （6分）证明：设λ为A 的按模最大特征值，x 为相对应的特征向量，则有Ax=λx 1分且λρ=)(A ，若λ是实数，则x 也是实数，得Ax x =λ 1分而x x ⋅=λλ x A x ,⋅≤⋅⋅≤λ故x A Ax 2分由于A x 0x ≤≠λ得到，两边除以 1分故A A ≤)(ρ 1分当λ是复数时，一般来说x 也是复数，上述结论依旧成立。

数值分析试题与答案

一、单项选择题（每小题3分，共15分）1. 和分别作为π的近似数具有（）和（）位有效数字.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 -=-单项选择题答案二、填空题（每小题3分，共15分）1. 设TX )4,3,2(-=, 则=1||||X ，2||||X = .2. 一阶均差()01,f x x =3. 已知3n =时，科茨系数()()()33301213,88C C C ===，那么()33C =4. 因为方程()420x f x x =-+=在区间[]1,2上满足，所以()0f x =在区间内有根。

5. 取步长0.1h =，用欧拉法解初值问题()211yy yx y ⎧'=+⎪⎨⎪=⎩的计算公式 .填空题答案0,1,2三、计算题（每题15分，共60分）1. 已知函数211y x =+的一组数据：求分段线性插值函数，并计算()1.5f 的近似值.计算题1.答案1. 解[]0,1x ∈，()1010.510.50110x x L x x --=⨯+⨯=---[]1,2x ∈，()210.50.20.30.81221x x L x x --=⨯+⨯=-+--所以分段线性插值函数为()[][]10.50,10.80.31,2x x L x x x ⎧-∈⎪=⎨-∈⎪⎩()1.50.80.3 1.50.35L =-⨯=2. 已知线性方程组1231231231027.21028.35 4.2x x x x x x x x x --=⎧⎪-+-=⎨⎪--+=⎩（1）写出雅可比迭代公式、高斯－塞德尔迭代公式；（2）对于初始值()()00,0,0X =，应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算()1X（保留小数点后五位数字）.计算题2.答案1.解原方程组同解变形为 1232133120.10.20.720.10.20.830.20.20.84x x x x x x x x x =++⎧⎪=-+⎨⎪=++⎩ 雅可比迭代公式为 ()()()()()()()()()1123121313120.10.20.720.10.20.830.20.20.84m m m m m m m m m x x x x x x x x x +++⎧=++⎪⎪=-+⎨⎪=++⎪⎩(0,1...)m = 高斯－塞德尔迭代法公式()()()()()()()()()1123112131113120.10.20.720.10.20.830.20.20.84m m m m m m m m m x x x x x x x x x ++++++⎧=++⎪⎪=-+⎨⎪=++⎪⎩ (0,1...)m = 用雅可比迭代公式得()()10.72000,0.83000,0.84000X = 用高斯－塞德尔迭代公式得()()10.72000,0.90200,1.16440X =3. 用牛顿法求方程3310x x --=在[]1,2之间的近似根（1）请指出为什么初值应取2（2）请用牛顿法求出近似根，精确到.计算题3.答案方程的根 1.87939x *≈4. 写出梯形公式和辛卜生公式，并用来分别计算积分1011dx x +⎰.计算题4.答案4 解梯形公式四、证明题（本题10分）确定下列求积公式中的待定系数，并证明确定后的求积公式具有3次代数精确度()()()()1010hhf x dx A f h A f A f h --=-++⎰证明题答案证明：求积公式中含有三个待定系数，即101,,A A A -，将()21,,f x x x =分别代入求积公式，一、填空（共20分，每题2分）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 。

数值分析试题及答案

数值分析试题及答案一、选择题（每题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. 求解线性方程组时，迭代法的基本思想是从一个初始近似解出发，通过不断________来逼近精确解。

数值分析期末试题及答案

数值分析期末考试试题一、单项选择题（每小题3分，共15分）1. 3.142和3.141分别作为π的近似数具有（）和（）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．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.A2.D3.D4.C5.B得分评卷人二、填空题（每小题3分，共15分）1. 设TX )4,3,2(-=, 则=1||||X ，2||||X = .2. 一阶均差()01,f x x =3. 已知3n =时，科茨系数()()()33301213,88C C C ===，那么()33C = 4. 因为方程()420x f x x =-+=在区间[]1,2上满足，所以()0f x =在区间内有根。

5. 取步长0.1h =，用欧拉法解初值问题()211yy yx y ⎧'=+⎪⎨⎪=⎩的计算公式 .填空题答案1. 9和292.()()0101f x f x x x --3. 18 4. ()()120f f < 5. ()1200.11.1,0,1,210.11k k y y k k y +⎧⎛⎫⎪ ⎪=+⎪ ⎪=+⎨⎝⎭⎪=⎪⎩得分评卷人三、计算题（每题15分，共60分）1. 已知函数211y x =+的一组数据：求分段线性插值函数，并计算()1.5f 的近似值.计算题1.答案()101x L x -=-()212x L x -=⨯-所以分段线性插值函数为()10.50.80.3x x L x x x ⎧-∈⎪=⎨-⎪⎩()1.50.8L =2. 已知线性方程组1231231231027.21028.35 4.2x x x x x x x x x --=⎧⎪-+-=⎨⎪--+=⎩（1）写出雅可比迭代公式、高斯－塞德尔迭代公式；（2）对于初始值()()00,0,0X =，应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算()1X（保留小数点后五位数字）.计算题2.答案1.解原方程组同解变形为 1232133120.10.20.720.10.20.830.20.20.84x x x x x x x x x =++⎧⎪=-+⎨⎪=++⎩雅可比迭代公式为()()()()()()()()()1123121313120.10.20.720.10.20.830.20.20.84m m m m m m m m m x x x x x x x x x +++⎧=++⎪⎪=-+⎨⎪=++⎪⎩(0,1...)m = 高斯－塞德尔迭代法公式()()()()()()()()()1123112131113120.10.20.720.10.20.830.20.20.84m m m m m m m m m x x x x x x x x x ++++++⎧=++⎪⎪=-+⎨⎪=++⎪⎩ (0,1...)m =用雅可比迭代公式得()()10.72000,0.83000,0.84000X =用高斯－塞德尔迭代公式得()()10.72000,0.90200,1.16440X =3. 用牛顿法求方程3310x x --=在[]1,2之间的近似根（1）请指出为什么初值应取2？（2）请用牛顿法求出近似根，精确到0.0001.计算题3.答案4. 写出梯形公式和辛卜生公式，并用来分别计算积分1011dx x +⎰.计算题4.答案确定下列求积公式中的待定系数，并证明确定后的求积公式具有3次代数精确度()()()()1010hhf x dx A f h A f A f h --=-++⎰证明题答案一、填空（共20分，每题2分）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 。

武汉大学数值分析期末考试题目和答案.pdf

注：1、教师命题时题目之间不留空白； 2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考n λ≥，
注：1、教师命题时题目之间不留空白；2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考
注：1、教师命题时题目之间不留空白；2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考
注：1、教师命题时题目之间不留空白；2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考
注：1、教师命题时题目之间不留空白；2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考
注：1、教师命题时题目之间不留空白；2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考。