数值计算方法期末复习答案终结版
一、 名词解释
1．误差：设*x 为准确值x 的一个近似值，称**()e x x x =-为近似值*x 的绝对误差，简称误差。
2．有效数字：有效数字是近似值的一种表示方法，它既能表示近似值的大小，又能表示其精确程度。如果近似值*x 的误差限是1
102
n -?，则称*x 准确到小数点后n 位，
并从第一个不是零的数字到这一位的所有数字均称为有效数字。
3. 算法：是指解题方案的准确而完整的描述，是一系列解决问题的清晰指令，算法代表着用系统的方法描述解决问题的策略机制。计算一个数学问题，需要预先设计好由已知数据计算问题结果的运算顺序，这就是算法。
4. 向量范数：设对任意向量n x R ∈ ，按一定的规则有一实数与之对应，记为||||x ，若||||x
满足
（1）||||0x ≥ ，且||||0x =
当且仅当0x = ；
（2）对任意实数α，都有||||||x αα= ||||x
； （3）对任意,n x y R ∈
，都有||||||||||||x y x y +≤+ 则称||||x
为向量x 的范数。
5. 插值法：给出函数()f x 的一些样点值，选定一个便于计算的函数形式，如多项式、分段
线性函数及三角多项式等，要求它通过已知样点，由此确定函数()x ?作为()f x 的近似的方法。
6相对误差：设*x 为准确值x 的一个近似值，称绝对误差与准确值之比为近似值*x 的相对误
差，记为*
()r e x ，即**
()()r e x e x x
=
7. 矩阵范数：对任意n 阶方阵A ，按一定的规则有一实数与之对应，记为||||A 。若||||A 满足
（1）||||0A ≥，且||||0A =当且仅当0A =； （2）对任意实数α，都有||||||A αα=||||A ；
（3）对任意两个n 阶方阵A,B,都有||||||||||||A B A B +≤+； （4）||||||||AB A =||||B 称||||A 为矩阵A 的范数。
8． 算子范数：设A 为n 阶方阵，||||?是n R 中的向量范数，则0
||||
||||||||max x Ax A x ≠= 是一种矩
阵范数，称其为由向量范数||||?诱导出的矩阵范数，也称算子范数。
9. 矩阵范数与向量范数的相容性：对任意n 维向量x
，都有
||||||||Ax A ≤ ||||x
这一性质称为矩阵范数与向量范数的相容性。
10. 1-范数，∞-范数和2-范数：
（1）1-范数 11
||||||
n
i i x x ==∑
（2）∞-范数 1||||m a x {||}
i i n
x x ∞
≤≤=
（3）2-范数
2||||x =
二、简答题
1．高斯消元法的思想是：先逐次消去变量，将方程组化成同解的上三角形方程组，此过程称为消元过程。然后按方程相反顺序求解上三角形方程组，得到原方程组的解，此过程称为回代过程。
2. 迭代法的基本思想是：构造一串收敛到解的序列，即建立一种从已有近似解计算新的近似解得规则，由不同的计算规则得到不同的迭代法。
3. 雅可比（Jacobi ）迭代法的计算过程（算法）：
（1）输入()ij A a =，1(,,)n b b b = ，维数n ，(0)(0)(0)(0)12(,,,)n x x x x =
，ε，最大容许迭代次数N 。 （2）置1k =
（3）对1,2,,i n = (0)
1
()/n
i i i j j i i
j j i x b a x a
=≠=-∑ （4）若(0)x x ε-<，输出x 停机；否则转5。
（5）k N <，置(0)1,(1,2,,)i i k k x x i n +??= ，转3，否则，输出失败信息，停机。 4. 插值多项式的误差估计：（P102）
由(1)(1)101()()()()()()()(1)!(1)!
n n n n n f f R x x x x x x x x n n ξξω+++==---++
当(0,1,,)i x x i n == 时，上式自然成立，因此，上式对[,]a b 上的任意点都成立，这就叫插值多项式的误差估计。
5. 反幂法的基本思想：设A 为阶非奇异矩阵，λ，u
为A 的特征值和相应的特征向量，
则1A - 的特征值是A 的特征值的倒数，而相应的特征向量不变，即
11A u u λ
-=
因此，若对矩阵1A -用幂法，，即可计算出1A -的按模最大的特征值，其倒数恰为A 的按模最小的特征值。
6. 雅可比（Jacobi ）迭代法是：选取初始向量(0)x
代入迭代公式
(1)()k k i x Bx g +=+
(0,1,2,k =
产生向量序列(){}k x
，由上述计算过程所给出的迭代法。 7. 数值计算中应注意的问题是： （1）避免两个相近的数相减 （2）避免大数“吃”小数的现象
（3）避免除数的绝对值远小于被除数的绝对值 （4）要简化计算，减少运算次数，提高效率 （5）选用数值稳定性好的算法
8. 高斯消去法的计算量：由消去法步骤知，在进行第k 次消元时，需作除法n k -次，乘法
()n k -(1)n k -+次，故消元过程中乘除运算总量为
乘法次数1
21()(1)(1)3n k n n k n k n -=--+=-∑ 除法次数1
1
()(1)2n k n
n k n -=-=-∑
在回代过程中，计算k x 需要(1)n k -+次乘除法，整个回代过程需要乘除运算的总量为
1
(1)(1)2
n
k n
n k n =-+=
+∑，所以，高斯消去法的乘除总运算量为 322(1)(1)(1)32233
n n n n n N n n n n =-+-++=+-
9. 迭代法的收敛条件：对任意初始向量(0)x 和右端项g
，由迭代格式
(1)()k k x Mx g +=+
(0,1,2,k =
产生的向量序列(){}k x
收敛的充要条件是()1M ρ<。
10. 迭代法的误差估计：设有迭代格式(1)()k k x Mx g +=+，若||||1M <，(){}k x
收敛于*x ，则有误差估计式()
*
(1)(0)||||||||||||1||||
K
k M x x x x M -≤--。
二、 计算题
1.假定运算中数据都精确到两位小数，试求*
1.21 3.659.81x =?-的绝对误差限和相对误差
限，计算结果有几位有效数字？
解：由式121212
12
121212()()()()()()r r r e x x e x e x x x e x x e x e x x x x x ±=±?
?
?±=±?±±?
和1221121212()()()()()()r r r e x x x e x x e x e x x e x e x ≈+??≈+?得 *() 3.65(1.21) 1.21(3.65)(9.81)e x e e e =?+?-
因为式中数据都精确到两位小数，即其误差限均为21
102
-?，故有
*|()| 3.65|(1.21)| 1.21|(3.65)||(9.81)|e x e e e ≤?+?+
**
*|()|0.0293
|()|0.0054|| 5.3935
r e x e x x =≤=
所以，*x 的绝对误差限为0.0293，相对误差限为0.0054，计算结果有两位有效数字。
2.求矩阵223477245A ??
=????-??
的三角分解。
解：由式111
111(1,2,,)(2,,,,,)()/(1,2,,1,1,,)
j j i ij ij ik kj k j ij ij ik kj jj
k u a j n u a l u i n j i n l a l u u j n i j n -=-=??==?
?
=-==?
?
?
=-=-=+??
∑∑
，12122u a ==，13133u a ==
2121114/22l a u ==
=，3131112
/12
l a u -===- 222221127223u a l u =-=-?=，232321137231u a l u =-=-?=
3232311222()/[4(1)2]/32l a l u u =-=--?= 333331133223()5[(1)321]6u a l u l u =-+=--?+?=
所以
100223210031121006A ????
????=????
????-????
21
(3.65 1.211)100.0293
2
-≤++??=
3．用幂法（2k =）求矩阵210021012A -??
=-????-??
的按模最大的特征值和相应的特征向量。取
(0)(0,0,0)T x =. （P 77） 解：(0)(0)(0,0,1)T y x ==
(1)(0)(0,1,2)T x Ay ==-, 2α=
(1)
(1)
(0,0.5,1)T x y
α
=
=-
(2)(1)(0.5,2,2.5)T x Ay ==-， 2.5
α=
4. 已知函数ln y x =，x 的值是10,11,12,13,14对应的ln y x =的值分别是 2.3026，2.3979, 2.4849, 2.5649, 2.6391。用Lagrange 线性插值求ln11.5的近似值。 解：取两个节点011x =，112x =，插值基函数为
1
001()(12)x x l x x x x -=
=--- 0110
()11x x l x x x x -==-- 由式01
10
1
0110
()x x x x x y y x x x x ?--=+--得 1() 2.3979(12) 2.4849(11)L x x x =--+-
将x=11.5代入，即得
1ln11.5(11.5) 2.39790.5 2.48490.5 2.4414L ≈=?+?=
按式(1)1()
()()(1)!
n n n f R x x n ξω++=+ (,)
a b ξ∈得 "1(ln )()(11)(12)2!
x R x x x ξ
=
--
因为"2
1
(ln )x x =-
，ξ在11和12之间，故 "2
2
1
1
|(ln )|0.008264511x ξξ=
≤
= 于是
311
|(11.5)|0.00826450.50.5 1.03306102
R -≤???=?
5. 用Jacobi 迭代法（1k =）求解线性方程组12312312
31027210283542
x x x x x x x x x --=??
-+-=??--+=? .
解：由Jacobi 迭代法得计算公式(1)
()11n
k k i
i
ij
j
j ii
ii
j i
b x
a x
a a +=≠=-
+
∑得 (1)()()
123(1)()()
213(1)()()3
120.10.27.20.10.28.30.20.28.4
k k k k k k k k k x x x x x x x x x +++?=++?=++??=++? 取(0)(0,0,0)T x =，代入上式得
(1)17.2x = (1)28.3x = (1)
38.4x = (2)
10.18.30.28.47.29.71x =?+?+= (2)20.17.20.28.48.310.70x =?+?+= (2)30.27.20.28.38.411.50x =?+?+=
6. 设有方程组Ax b =，其中1
11221
112211122
A ?
?????
=?
????????
?
，讨论用Jacobi 迭代法求解的收敛性。 解：因为A 为对称矩阵，且其各阶主子式皆大于零，故A 为对称正定矩阵，A 不是弱对角
占优阵，故不能判别Jacobi 迭代的收敛性。易算出Jacobi 迭代法的迭代矩阵为
11
102
21
10221102
2
B I D A -??--??????=-=-
-??????--????
其特征方程
31
12
21
113||2244112
2
I B λλλ
λλλ???????
?-==+-??????????
21
()(1)02
λλ=-+=
有根121
2
λλ==
，31λ=-，因而()1B ρ=。由向量序列(){}k x 收敛的充要条件是()1B ρ<，故Jacobi 迭代法不收敛。
7．用反幂法（1k =）求矩阵210021012A -??
=-????-??
接近2.93的特征值，并求相应的特征向量，取
(0)(0,0,0)T x =.
解：对 2.93A I -作三角分解得
0.93
102.9300.931010.93A I --????-=--??
??--??
1000.93100
1000.9311101000.930.930.93????????--?
???=--????????-+???????
?
8. 已知函数ln y x =，x 的值是10,11,12,13,14对应的ln y x =的值分别是 2.3026，2.3979, 2.4849, 2.5649, 2.6391。用Lagrange 抛物线插值求ln11.5的近似值。 解：取011x =，112x =，213x =，插值多项式为
2(12)(13)(11)(13)(11)(12)
() 2.3979
2.4849 2.5649(1112)(1113)(1211)(1213)(1311)(1312)
x x x x x x L x ------=++------
1.19895(12)(13)
2.4849(11)(13) 1.28245(11)(12)x x x x x x =-----+--
所以
2ln11.5(11.5)L ≈
1.19895(0.5)( 1.5)
2.48490.5( 1.5) 1.282450.5(0.5) 2.442275=?-?--??-+??-=
因为"'3
2
(ln )x x =
，于是 "'2
31113
2max |(ln )|0.15031011
x x -≤≤≤
=? 因此用抛物线插值法计算的误差为
"'2|(ln )||(11.5)||(11.511)(11.512)(11.513)|3!
x R ξ=
---
251
0.1503100.50.5 1.59.3938106
--≤?????=? 查表可得ln11.5 2.442347= 三、 证明题
1. 若x 的近似值
x *=1210.10(0)m n a a a a ±?≠…有n 位有效数字，则11
1102n a -+?为其相对误差限。反之，若x *
的相对误差限r ε满足111
102(1)
n r a ε-+≤
?+，则x *至少具有n 位有效数字。
证明：由式*1
||102
m n x x --≤?得
**1
|()|||102
m n e x x x -=-≤?
从而有
*
*1*121
1
10()12|()|||100.102m n
n r m n e x e x x a a a a --+?=≤≤?? 所以
11
1
102n a -+?是*x 的相对误差限。 若1
11102(1)
n r a ε-+≤?+，由式***()|()|||r r e x e x x ε=≤得
***12|()||()|0.10m r n r e x x e x a a a ε=≤?
111111
(1)1010102(1)2
m n m n a a --+-≤+??
?=?+
由式*1
||102
m n x x --≤?，*x 至少有n 位有效数字。
2. 设01,,,n x x x …为1n +个互异节点，(),(0,1,)i l x i =…,n 为这组点上的Lagrange 插值基函数，试证明0()1n
i i l x =≡∑。
证明：上式的左端为插值基函数的线性组合，其组合系数均为1。显然，函数()1f x ≡在这
n+1个节点处取值均为1，即()1i i y f x == (0,1,,i n =
，由式0()()n
n i i i L x y l x ==∑知，它的n 次Lagrange 插值多项式为
0()()n
n i i L x l x ==∑
对任意x ，插值余项为
(1)1()
()()()()0(1)!
n n n n f R x f x L x x n ξω++=-=≡+
所以 0()()()1
n
n i i L x l x f x ==≡=∑
3 设A 为任意n 阶方阵，?为任意由向量范数诱导出的矩阵范数，则()A A ρ≤
证明：对A 的任一特征值i λ及相应的特征向量i u
，都有
||i λ||||||||||||||||i i i i u u Au A λ==≤ ||||i u
因为i u
为非零向量，于是有 ||||||
i A λ≤ 由i λ的任意性即得 ()||||A A
ρ≤
4. 设A 为n 阶方阵，则lim 0k k A →∞
=的充分必要条件为()1A ρ<。
证明：必要性。若
lim 0k k A →∞
=
由相关定义得 l i m ||||k k A
→∞
= 而 0()[()]||K K K A A A ρρ≤=≤
于是由极限存在准则，有 l i m [()]k k A ρ
→∞
= 所以()1A ρ<。
充分性。若()1A ρ<，取1()
02
A ρε-=
>，由||||()A A ρε≤+,存在一种矩阵范数?，使得
1()
||||()12
A A A ρρε+≤+=
< 而||||||||k k A A ≤，于是 l i m ||||l i m |||
k k k k A
A →∞
→∞
== 所以 l i m
0k k A →∞
=
五、应用题
1.平面桁架是由刚性元件通过结点互相联结而组成的力学结构，它通常出现在桥梁结构和其他需要力学支撑的结构中。如图是一个简单的静力桁架结构，其中刚性元件（5m =）通过
结点,,,A B C D 相连。求各个结点的合力方程，并求出当,3
6
π
π
αβ=
=
外部负荷
12250,1500g N g N ==时，求各个节点内力。
解：设五个刚性元件的内力为125{,,,}f f f ，它们都处理为压力，如果解是负的，表明该力是张力。桁架的左边由固定结点A 支撑，右边由滑轮D 支撑，678,,f f f 是外部支撑力，12,g g 是外部负荷。由于在静力平衡时，每个结点处的水平方向合力与垂直方向的合力为零，那么有
结点A 12617cos 0sin 0f f f f f αα+-=??+=? 结点B 141134cos cos 0
sin sin 0f f g f f f αβαβ-++=??---=?
结点C 253200f f f g -+=??-=? 结点D 454
8cos 0
sin 0f f f f ββ--=??+=?
设f
表示未知力向量，上述方程组可用矩阵表示为
12cos 100
01
000sin 00000100cos 00cos 0000sin 01sin 0000001001000000100000000cos 10000000sin 000
10g f g αααβαβββ-????
?????
???????--?
???---????=????-????????
????--?????????
??
?
若取,3
6
π
π
αβ==
，外部负荷12250,1500g N g N ==。采用列主元素法，得各结点的内力如
下：
(1174,837,1500,966.5,837,250,1017,483.3)T f =--
数值计算方法试题及答案
【 数值计算方法试题一 一、 填空题（每空1分，共17分） 1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数，需对分（ ）次。 2、迭代格式)2(2 1-+=+k k k x x x α局部收敛的充分条件是α取值在（ ）。 3、已知?????≤≤+-+-+-≤≤=31)1()1()1(211 0)(2 33x c x b x a x x x x S 是三次样条函数， 则 a =( )， b =（ ）， c =（ ）。 4、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数，则 ∑== n k k x l 0)(( )， ∑== n k k j k x l x 0 )(( )，当2≥n 时 = ++∑=)()3(20 4x l x x k k n k k ( )。 ； 5、设1326)(2 47+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 和=?07 f 。 6、5个节点的牛顿-柯特斯求积公式的代数精度为 ，5个节点的求积公式最高代数精度为 。 7、{}∞ =0)(k k x ?是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x ?，则?= 1 4)(dx x x ? 。 8、给定方程组?? ?=+-=-2211 21b x ax b ax x ，a 为实数，当a 满足 ，且20<<ω时，SOR 迭代法收敛。 9、解初值问题 00 (,)()y f x y y x y '=?? =?的改进欧拉法 ??? ??++=+=++++)],(),([2),(] 0[111] 0[1n n n n n n n n n n y x f y x f h y y y x hf y y 是 阶方法。
北师大网络教育 数值分析 期末试卷含答案
注：1、教师命题时题目之间不留空白； 2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考北师大网络教育——数值分析——期末考试卷与答案 一．填空题（本大题共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 =的牛顿-科特斯求积公式计算定积分1 01 1I dx x =+? 。（10分） 五．用Newton 法求()cos 0f x x x =-=的近似解。（10分） 六．试用Doolittle 分解法求解方程组：
注：1、教师命题时题目之间不留空白； 2、考生不得在试题纸上答题，教师只批阅答题册正面部分，若考 12325610413191963630 x x x -?????? ??????-=?????? ??????----?????? （10分） 七．请写出雅可比迭代法求解线性方程组1231231 23202324 812231530 x x x x x x x x x ++=?? ++=??-+=? 的迭代格式，并 判断其是否收敛？（10分） 八．就初值问题0(0)y y y y λ'=??=?考察欧拉显式格式的收敛性。（10分）
数值计算方法I上机实验考试题
数值计算方法I 上机实验考试题（两题任选一题） １．小型火箭初始质量为900千克，其中包括600千克燃料。火箭竖直向上发射时燃料以15千克/秒的速率燃烧掉，由此产生30000牛顿的恒定推力．当燃料用尽时引擎关闭。设火箭上升的整个过程中，空气阻力与速度平方成正比，比例系数为0.4（千克/米）．重力加速度取9.8米/秒2. A. 建立火箭升空过程的数学模型（微分方程）； B. 求引擎关闭瞬间火箭的高度、速度、加速度，及火箭到达最高点的时间和高度． ２．小型火箭初始质量为1200千克，其中包括900千克燃料。火箭竖直向上发射时燃料以15千克/秒的速率燃烧掉，由此产生40000牛顿的恒定推力．当燃料用尽时引擎关闭。设火箭上升的整个过程中，空气阻力与速度平方成正比，比例系数记作k ，火箭升空过程的数学模型为 0)0(,0,01222==≤≤-+?? ? ??-==t dt dx x t t mg T dt dx k dt x d m 其中)(t x 为火箭在时刻t 的高度，m =1200-15t 为火箭在时刻t 的质量，T （=30000牛顿）为推力，g (=9.8米/秒2)为重力加速度, t 1 (=900/15=60秒)为引擎关闭时刻． 今测得一组数据如下（t ~时间（秒），x ~高度（米），v ~速度（米/秒））： 现有两种估计比例系数k 的方法： 1．用每一个数据(t,x,v )计算一个k 的估计值（共11个），再用它们来估计k 。 2．用这组数据拟合一个k ． 请你分别用这两种方法给出k 的估计值，对方法进行评价，并且回答，能否认为空气阻力系数k=0.5（说明理由）．
现代数值计算方法习题答
现代数值计算方法习题答案 习 题 一 1、解：根据绝对误差限不超过末位数的半个单位，相对误差限为绝对误差限除以 有效数字本身，有效数字的位数根据有效数字的定义来求.因此 49×10 -2 ：E = 0.005； r E = 0.0102； 2位有效数字. 0.0490 ：E = 0.00005；r E = 0.00102； 3位有效数字. 490.00 ：E = 0.005； r E = 0.0000102；5位有效数字. 2、解： 7 22 = 3.1428 …… ， π = 3.1415 …… ， 取它们的相同部分3.14，故有3位有效数字. E = 3.1428 - 3.1415 = 0.0013 ；r E = 14 .3E = 14 .30013.0 = 0.00041. 3、解：101的近似值的首位非0数字1α = 1,因此有 |)(*x E r |) 1(10 1 21--??=n < = 2 1× 10 -4 , 解之得n > = 5，所以 n = 5 . 4、证：) ()(1)()(1)(* 1 1* * 1 1 * * x x x n x E x n x E n n n -= ≈ -- )(11)()(1) ()(* * * * * 1 1 ** * * x E n x x x n x x x x n x x E x E r n n n n n r = -= -≈ = - 5、解：(1)因为=20 4.4721…… ， 又=)(*x E |*x x -| = |47.420-| = 0.0021 < 0.01, 所以 =*x 4.47. (2)20的近似值的首位非0数字1α = 4,因此有 |)(*x E r |) 1(10 4 21--??= n < = 0.01 , 解之得n > = 3 .所以，=*x 4.47. 6、解：设正方形的边长为x ,则其面积为2x y =，由题设知x 的近似值为*x = 10 cm . 记*y 为y 的近似值，则
数值计算方法试题及答案
数值计算方法试题一 一、填空题（每空1分，共17分） 1、如果用二分法求方程在区间内的根精确到三位小数，需对分（）次。 2、迭代格式局部收敛的充分条件是取值在（）。 3、已知是三次样条函数，则 =( )，=（），=（）。 4、是以整数点为节点的Lagrange插值基函数，则 ( )，( )，当时( )。 5、设和节点则 和。 6、5个节点的牛顿-柯特斯求积公式的代数精度为，5个节点的求积公式最高代数精度为。 7、是区间上权函数的最高项系数为1的正交多项式族，其中，则。 8、给定方程组，为实数，当满足，且时，SOR迭代法收敛。 9、解初值问题的改进欧拉法是 阶方法。 10、设，当（）时，必有分解式，其中为下三角阵，当其对角线元素满足（）条件时，这种分解是唯一的。 二、二、选择题（每题2分） 1、解方程组的简单迭代格式收敛的充要条件是（）。（1）, (2) , (3) , (4) 2、在牛顿-柯特斯求积公式：中，当系数是负值时，公式的稳定性不能保证，所以实际应用中，当（）时的牛顿-柯特斯求积公式不使用。 （1），（2），（3），（4）， （1）二次；（2）三次；（3）四次；（4）五次 4、若用二阶中点公式求解初值问题，试问为保证该公式绝对稳定，步长的取值范围为（）。 (1), (2), (3), (4)
三、1、 2、（15 (1)(1) 试用余项估计其误差。 （2）用的复化梯形公式（或复化 Simpson公式）计算出该积分的近似值。 四、1、（15分）方程在附近有根，把方程写成三种不同的等价形式（1）对应迭代格式；(2)对应迭代格式；（3）对应迭代格式。判断迭代格式在的收敛性，选一种收敛格式计算附近的根，精确到小数点后第三位。选一种迭代格式建立Steffensen迭代法，并进行计算与前一种结果比较，说明是否有加速效果。 2、（8分）已知方程组，其中 ， （1）（1）列出Jacobi迭代法和Gauss-Seidel迭代法的分量形式。 （2）（2）求出Jacobi迭代矩阵的谱半径，写出SOR 迭代法。 五、1、（15分）取步长，求解初值问题用改进的欧拉法求的值；用经典的四阶龙格—库塔法求的值。 2、（8分）求一次数不高于4次的多项式使它满足 ,,,, 六、（下列2题任选一题，4分） 1、1、数值积分公式形如 （1）（1）试确定参数使公式代数精度尽量高；（2）设，推导余项公式，并估计误差。 2、2、用二步法 求解常微分方程的初值问题时，如何选择参数使方法阶数尽可能高，并求局部截断误差主项，此时该方法是几阶的。 数值计算方法试题二 一、判断题：（共16分，每小题２分） １、若是阶非奇异阵，则必存在单位下三角阵和上三角阵，使唯一成立。（）
数值分析学期期末考试试题与答案(A)
期末考试试卷（A 卷） 2007学年第二学期 考试科目： 数值分析 考试时间：120 分钟 学号 姓名 年级专业 一、判断题（每小题2分，共10分） 1. 用计算机求 1000 1000 1 1 n n =∑时，应按照n 从小到大的顺序相加。 （ ） 2. 为了减少误差,进行计算。 （ ） 3. 用数值微分公式中求导数值时，步长越小计算就越精确。 （ ） 4. 采用龙格－库塔法求解常微分方程的初值问题时，公式阶数越高，数值解越精确。（ ） 5. 用迭代法解线性方程组时，迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有 关，与常数项无关。 （ ） 二、填空题（每空2分，共36分） 1. 已知数a 的有效数为0.01，则它的绝对误差限为________，相对误差限为_________. 2. 设1010021,5,1301A x -????????=-=-????????-???? 则1A =_____，2x =______，Ax ∞ =_____. 3. 已知5 3 ()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---= . 4. 为使求积公式 1 1231 ()()(0)33 f x dx A f A f A f -≈- ++? 的代数精度尽量高，应使1A = ，2A = ，3A = ，此时公式具有 次的代数精度。 5. n 阶方阵A 的谱半径()A ρ与它的任意一种范数A 的关系是 . 6. 用迭代法解线性方程组AX B =时，使迭代公式(1) ()(0,1,2,)k k X MX N k +=+=产 生的向量序列{ }() k X 收敛的充分必要条件是 . 7. 使用消元法解线性方程组AX B =时，系数矩阵A 可以分解为下三角矩阵L 和上三角矩
《数值计算方法》上机实验报告
《数值计算方法》上机实验报告华北电力大学 实验名称数值il?算方法》上机实验课程名称数值计算方法专业班级：电力实08学生姓名：李超然学号:200801001008 成绩: 指导教师:郝育黔老师实验日期：2010年04月华北电力大学实验报告数值计算方法上机实验报吿一. 各算法的算法原理及计算机程序框图1、牛顿法求解非线性方程 *对于非线性方程，若已知根的一个近似值，将在处展开成一阶 xxfx ()0, fx ()xkk 泰勒公式 "f 0 / 2 八八,fxfxfxxxxx 0 0 0 0 0 kkkk2! 忽略高次项，有 ，fxfxfxxx 0 ()()()，，, kkk 右端是直线方程，用这个直线方程来近似非线性方程。将非线性方程的 **根代入，即fx ()0, X ，* fxfxxx 0 0 0 0, ,, kkk fx 0 fx 0 0,
解出 fX 0 *k XX,, k' fx 0 k 水将右端取为，则是比更接近于的近似值，即xxxxk, Ik, Ik fx ()k 八XX, Ikk* fx()k 这就是牛顿迭代公式。 ,2,计算机程序框图:，见, ,3,输入变量、输出变量说明: X输入变量:迭代初值，迭代精度，迭代最大次数，\0 输出变量:当前迭代次数，当前迭代值xkl ,4,具体算例及求解结果: 2/16 华北电力大学实验报吿 开始 读入 l>k /fx()0?,0 fx 0 Oxx，，01* fx ()0 XX,,，?10 kk, ,1，kN, ?xx, 10 输出迭代输出X输出奇异标志1失败标志
,3,输入变量、输出变量说明: 结束 例：导出计算的牛顿迭代公式，并il ?算。（课本P39例2-16） 115cc （0）, 求解结果: 10. 750000 10.723837 10. 723805 10. 723805 2、列主元素消去法求解线性方程组，1,算法原理: 高斯消去法是利用现行方程组初等变换中的一种变换，即用一个不为零的数乘 -个 方程后加只另一个方程，使方程组变成同解的上三角方程组，然后再自下而上 对上三角 3/16 华北电力大学实验报告方程组求解。 列选主元是当高斯消元到第步时，从列的以下（包括）的各元素中选出绝 aakkkkkk 对值最大的，然后通过行交换将其交换到的位置上。交换系数矩阵中的 两行（包括常ekk 数项），只相当于两个方程的位置交换了，因此，列选主元不影响求解的结 ,2,计算机程序框图:，见下页, 输入变量:系数矩阵元素，常向量元素baiji 输出变量:解向量元素bbb,,12n
现代数值分析
研硕16《化工数值方法及Matlab应用》试题 班级姓名成绩 1．（15分）数值计算方法的主要研究对象有哪些？其常用基本算法主要包括哪三个方面？举例说明Matlab在解决化工数值计算问题方面有什么样实用价值？答：（1）数值计算方法的主要研究对象为非线性方程求根，插值法、曲线拟合、数值积分、常微分方程（组）、初值问题求解、线性和非线性方程组求解。（2）基本算法包括①离散化方法：用差商代替导数、差分代替微分等，将连续的数学问题转化为离散问题。②逼近方法：用简单函数的值近似代替求解困难或形式未知的复杂函数的值。③迭代法：用一个固定公式反复计算，对较为粗糙的根的近似值进行加工直到满足精度要求的方法。 （3）Matlab在解决化工数值计算问题的实用价值有：数值计算和符号计算功能；图形功能；MATLAB语言；功能性和学科性工具箱。 2．（10分）数值计算中的“曲线拟合”，一般有哪些方法？请至少指出四种，并简述各自的基本特点。 答：（1）拉格朗日插值：，优点在于不要求数据点事等间隔的，缺点是数据点不易过多，当数据比较多时，差值函数有偏离原函数的风险； （2）牛顿插值法：它不仅克服了“增加一个节点时整个计算工作必须重新开始”的缺点，而且可以节省乘、除法运算次数。同时，在牛顿插值多项式中用到的差分与差商等概念，又与数值计算的其他方面有着密切的关系。
（3）牛顿迭代法：牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根，此时线性收敛，但是可通过一些方法变成超线性收敛。 （4）区间二分法：优点：算法简单，容易理解，且总是收敛的。缺点：收敛速度太慢，浪费时间，二分法不能求复根跟偶数重根。 （5）最小二乘法：通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。 3. （15分）在298K 下，化学反应 2OF 2=O 2+2F 2 的平衡常数为0.410 atm ，如在298K 下将OF 2 通入容器，当t=0 时为1 atm ，问最后总压是多少？取计算精度为10-3 。 解：首先写出求解问题的数学方程式。 假设气体是理想气体，由反应的化学计量式可知， 22222F O OF += 设氧的分压为p ，平衡时有p 21- p p 2。 平衡时，有()410.02142 3=-p p 整理得 0410.064.1640.1423=-+-p p p 函数关系为 ()0410.064.1640.1423=-+-=p p p p f
数值分析期末考试复习题及其答案.doc
数值分析期末考试复习题及其答案 1. 已知325413.0,325413* 2* 1==X X 都有6位有效数字，求绝对误差限。（4分） 解： 由已知可知,n=6 5.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==A 3. 设3 2 )()(a x x f -= （6分） ① 写出f(x)=0解的Newton 迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （k=0,1……）产生的序列{}k x 收敛于2 解： ①Newton 迭代格式为： x a x x x a x 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)()(')(2 2 32 1 += +=---=-=+? 3分
②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-=a a x a x ?? 3分 4. 给定线性方程组Ax=b ，其中：? ??=1 3A ??? 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)(数值分析习题与答案
第一章绪论 习题一 1.设x>0,x*的相对误差为δ，求f(x)=ln x的误差限。解：求lnx的误差极限就是求f(x)=lnx的误差限，由公式(1. 2.4)有 已知x*的相对误差满足，而 ，故 即 2.下列各数都是经过四舍五入得到的近似值，试指出它们有几位有效数字，并给出其误差限与相对误差限。 解：直接根据定义和式(1.2.2)(1.2.3)则得 有5位有效数字，其误差限，相对误差限 有2位有效数字， 有5位有效数字， 3.下列公式如何才比较准确？ （1） （2）
解：要使计算较准确，主要是避免两相近数相减，故应变换所给公式。 （1） （2） 4.近似数x*=0.0310,是 3 位有数数字。 5.计算取，利用：式计算误差最小。 四个选项： 第二、三章插值与函数逼近 习题二、三 1. 给定的数值表 用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解：仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计（5.8）。线性插值时，用0.5及0.6两点，用Newton插值 误差限，因
，故 二次插值时，用0.5，0.6，0.7三点，作二次Newton插值 误差限 ，故 2. 在-4≤x≤4上给出的等距节点函数表，若用二次插值法求的近似值，要使误差不超过，函数表的步长h 应取多少? 解：用误差估计式（5.8）， 令 因 得 3. 若，求和.
解：由均差与导数关系 于是 4. 若互异，求 的值，这里p≤n+1. 解：，由均差对称性 可知当有 而当P＝n＋1时 于是得 5. 求证. 解：解：只要按差分定义直接展开得 6. 已知的函数表
数值计算方法》试题集及答案
《计算方法》期中复习试题 一、填空题： 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：2.367，0.25 2、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2 x 的系数为 ，拉 格朗日插值多项式为 。 答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )； 答案 )(1)(1n n n n n x f x f x x x '--- =+ 5、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 6、计算方法主要研究( 截断 )误差和( 舍入 )误差； 7、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为 ( 1 2+-n a b )； 8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f )，代数精度 为( 5 )； 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该表达 式改写为 11 ,))64(3(10-= -++=x t t t t y ，为了减少舍入误差，应将表达式1999 2001-
MATLAB与数值分析实验报告一
MATLAB与数值分析实验报告 报告人：秦旸照 学号： 2015020901033 时间： 2016.4.8 电子科技大学电子工程学院
一、实验目的 实验一：MATLAB软件平台与程序设计实验 二、实验原理 1.熟练掌握矩阵的生成、加、减、乘、除、转置、行列式、逆、范数等运算操作。（用.m文件和Matlab函数编写一个对给定矩阵进行运算操作的程序） 2. 熟练掌握算术符号操作和基本运算操作，包括矩阵合并、向量合并、符号转换、展开符号表达式、符号因式分解、符号表达式的化简、代数方程的符号解析解、特征多项式、函数的反函数、函数计算器、微积分、常微分方程的符号解、符号函数的画图等。（用.m文件编写进行符号因式分解和函数求反的程序） 3. 掌握Matlab函数的编写规范。 4.掌握Matlab常用的绘图处理操作，包括：基本平面图、图形注释命令、三维曲线和面的填充、三维等高线等。（用.m文件编写在一个图形窗口上绘制正弦和余弦函数的图形，并给出充分的图形注释） 5. 熟练操作MATLAB软件平台，能利用M文件完成MATLAB的程序设计。 三、实验方案 1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以 x，y为坐标显示图像 x(n+1) = a*x(n)-b*(y(n)-x(n)^2); y(n+1) = b*x(n)+a*(y(n)-x(n)^2) 2. 编程实现奥运5环图，允许用户输入环的直径。 3. 实现对输入任意长度向量元素的冒泡排序的升序排列。 不允许使用sort函数。 四、实验结果 1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以 x，y为坐标显示图像
现代数值分析复习题
复习题(一) 一、填空题： 1、求方程0.5x2 101x 1 0的根，要求结果至少具有6位有效数字。已知 V10203 101.0099，贝卩两个根为x1 _____________________________ ， X2 ________________________________ .(要有计算过程和结果) 4 1 0 A A 1 4 1 2、0 1 4，则A的LU分解为。 1 2 A 3、 3 5，贝卩(A) ____________ ，A __________ . 4、已知f(1)「Q f(2)「2，f(3) =3，则用抛物线(辛卜生)公式计算求 3 得1 f(x)dx -------------------- ，用三点式求得f (1) ________________ . 5、f(1) 1，f(2) 2，f(3) 1，则过这三点的二次插值多项式中x2的系数 为_____ ，拉格朗日插值多项式为 _________________________ . 二、单项选择题： 1、Jacobi迭代法解方程组Ax b的必要条件是( ). A. A的各阶顺序主子式不为零 B. (A) 1 C a ii 0,i 1,2, ,n D|| A 1 2、设f(x) 3x99 5x 7，均差f[1,2,22, ,299]=(). D. 3
4、三点的高斯求积公式的代数精度为 （ ). A.3 B. -3 C. 5 D.0 2 2 3 A 0 5 1 3、设 0 0 7 ，则 (A )为( ). A. 2 B. 5 C. 7
分别用拉格朗日插值法和牛顿插值法求 f （x ）的三次插值多项式P 3（x ），并 求f （2）的近似值（保留四位小数）. 4、 取步长h 0.2,用预估-校正法解常微分方程初值问题 y 2x 3y y （0） 1 （0 x 1） 5、 已知 A. 2 B.5 C. 3 D. 4 5、幕法的收敛速度与特征值的分布 A.有关 B.不一定 C. 无关 三、计算题: 1、用高斯-塞德尔方法解方程组 4X ! 2X 2 X 3 11 X 1 4X 2 2X 3 18 2X ! X 2 5X 3 22 （°） /c c c\T ，取 x （°，°，°），迭 四次（要求按五位有效数字计算 ). 1 2、求A 、B 使求积公式 1 f （X ）dX A[f( 1) f (1)] 1 B [f (2)f （2）］ 的代数精 度尽量高，并求其代数精度；利用此公式求 I 21dx 1 x （保留四位小 数）。 3、已知
数值分析作业答案
数值分析作业答案 插值法 1、当x=1，-1，2时，f(x)=0，-3，4，求f(x)的二次插值多项式。 （1）用单项式基底。 （2）用Lagrange插值基底。 （3）用Newton基底。 证明三种方法得到的多项式是相同的。 解：（1）用单项式基底 设多项式为: ， 所以： 所以f(x)的二次插值多项式为： （2）用Lagrange插值基底 Lagrange插值多项式为： 所以f(x)的二次插值多项式为： (3) 用Newton基底: 均差表如下： xk f(xk) 一阶均差二阶均差 1 0 -1 -3 3/2 2 4 7/ 3 5/6 Newton插值多项式为： 所以f(x)的二次插值多项式为： 由以上计算可知，三种方法得到的多项式是相同的。 6、在上给出的等距节点函数表，若用二次插值求ex的近似值，要使截断误差不超过10-6，问使用函数表的步长h应取多少？ 解：以xi-1,xi,xi+1为插值节点多项式的截断误差，则有 式中 令得 插值点个数
是奇数，故实际可采用的函数值表步长 8、，求及。 解：由均差的性质可知，均差与导数有如下关系： 所以有： 15、证明两点三次Hermite插值余项是 并由此求出分段三次Hermite插值的误差限。 证明：利用[xk，xk+1]上两点三次Hermite插值条件 知有二重零点xk和k+1。设 确定函数k(x): 当或xk+1时k(x)取任何有限值均可； 当时，，构造关于变量t的函数 显然有 在[xk,x][x,xk+1]上对g(x)使用Rolle定理，存在及使得 在，，上对使用Rolle定理，存在，和使得 再依次对和使用Rolle定理，知至少存在使得 而，将代入，得到 推导过程表明依赖于及x 综合以上过程有： 确定误差限： 记为f(x)在[a,b]上基于等距节点的分段三次Hermite插值函数。在区间[xk，xk+1]上有 而最值 进而得误差估计： 16、求一个次数不高于4次的多项式，使它满足，，。
数值计算方法期末考试题
一、单项选择题（每小题3分，共15分） 1. 3.142和3.141分别作为的近似数具有（ ）和（ ）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和4 2. 已知求积公式 ，则＝（ ） A ． B ． C ． D ． 3. 通过点 的拉格朗日插值基函数满足（ ） A ． ＝0， B ． ＝0， C ．＝1， D ． ＝1， 4. 设求方程 的根的牛顿法收敛，则它具有（ ）敛速。 A ．超线性 B ．平方 C ．线性 D ．三次 5. 用列主元消元法解线性方程组 作第一次消元后得到的第3个方程（ ）. A ． B ． C ． D ． π()()2 1 121 1()(2)636f x dx f Af f ≈ ++? A 1613122 3()()0011,,,x y x y ()()01,l x l x ()00l x ()110l x =() 00l x ()111 l x =() 00l x ()111 l x =() 00l x ()111 l x =()0 f x =12312312 20 223332 x x x x x x x x ++=?? ++=??--=?232 x x -+=232 1.5 3.5 x x -+=2323 x x -+=
单项选择题答案 1.A 2.D 3.D 4.C 5.B 二、填空题（每小题3分，共15分） 1. 设, 则 ， . 2. 一阶均差 3. 已知时，科茨系数 ，那么 4. 因为方程 在区间 上满 足 ，所以 在区间内有根。 5. 取步长，用欧拉法解初值问题 的计算公 式 . 填空题答案 230.5 1.5 x x -=-T X )4,3,2(-==1||||X 2||||X =()01,f x x = 3n =()()() 33301213,88C C C === () 3 3C =()420 x f x x =-+=[]1,2()0 f x =0.1h =()211y y y x y ?'=+?? ?=?
数值计算方法实验报告(例)讲解
实验报告 一、实验目的 二、实验内容 三、实验环境 四.实验方法 五、实验过程 1实验步骤 2 关键代码及其解释 3 调试过程 六、实验总结 1．遇到的问题及解决过程 2．产生的错误及原因分析 3．体会和收获。 七、程序源代码： 八、教师评语
实验报告 一.试验目的：练习用数值方法求解给定的非线性方程。 二.实验内容：求解人口方程: )1(5 .43e 1004.156-+ =λλλ e 要求误差小于410-。 三.实验环境：PC 计算机，FORTRAN 、C 、C ++、VB 任选一种。 四.实验方法：牛顿法 牛顿法简述：牛顿法是一种特殊的迭代法，其迭代公式为： ,2,1,0,) () (1='- =+k x f x f x x k k k k ， 当数列{}k x 收敛时，其极限值x 即为方程的解。 定理：给定方程],[,0)(b a x x f ∈= 1）设0)()(''x f x f ； 则牛顿法产生的序列{}k x 收敛于0)(=x f 在],[b a 内的唯一解x 。 五.实验过程: 1．编程: 用C 语言编出牛顿法的源程序。 2. 开机, 打开C 语言编译程序，键入所编程序源代码. 3. 调试程序, 修改错误至能正确运行. 六．实验总结: （1）牛顿法收敛速度快，但初值不容易确定，往往由于初值取得
不当而使迭代不收敛或收敛慢，但若能保证)()(1+>K K x f x f （称为下山条件），则有可能收敛。把新的近似值看作初值的话会比原来的取得好，有可能落入局部收敛的邻域。 （2）牛顿法要求)(x f '在x 附近不为零。亦即x 只能是单根, 不能求重根。可用重根加速收敛法求重根。 （3）牛顿法的每一步迭代中，都要计算一次导数值，若计算)(x f '比计算函数的近似值要麻烦的多。为了避免求导数，可用差商近似代替微商 1 1) ()()(----='K K K K K x x x f x f x f 此时牛顿迭代法改为 )() ()() (111--+--- =K K K K K K K x x x f x f x f x x . (4) 由于人口方程来源于实际问题, λ代表人口增长率, 其真实 值不会太大, 初值不应取得过大.否则会得到该方程的另外一个解 七、程序源代码： #include 1.已知 ln(2.0)=0.6931;ln(2.2)=0.7885,ln(2.3)=0 .8329,试用线性插值和抛物插值计算.ln2.1的值并估计误差 2.已知x=0,2,3,5对应的函数值分别为y=1,3,2,5.试求三次多项式的插值 3. 分别求满足习题1和习题2 中插值条件的Newton插值 (1) (2) 3()1(2)(2)(3) 310 N x x x x x x x =+--+--4. 给出函数f(x)的数表如下,求四次Newton 插值多项式,并由此计算f(0.596)的值 解： 5.已知函数y=sinx的数表如下,分别用前插和后插公式计算sin0.57891的值 6.求最小二乘拟合一次、二次和三次多项式，拟合如下数据并画出数据点以及拟合函数的图形。 (a) (b) 7.试分别确定用复化梯形、辛浦生和中矩形 求积公式计算积分2 14dx x +?所需的步长h ，使得精度达到5 10 -。 8.求A 、B 使求积公式 ?-+-++-≈1 1)] 21()21([)]1()1([)(f f B f f A dx x f 的 代数精度尽量高,并求其代数精度；利用 此公式求? =2 1 1dx x I (保留四位小数)。 9.已知 分别用拉格朗日插值法和牛顿插值法求 ) (x f 的三次插值多项式)(3 x P ，并求)2(f 的近 似值（保留四位小数）。 10.已知 求)(x f 的二次拟合曲线)(2 x p ，并求)0(f 的近似值。 11.已知x sin 区间[0.4，0.8]的函数表 《数值计算方法》复习试题 一、填空题： 1、?? ??? ?????----=410141014A ，则A 的LU 分解为 A ??? ?????????=????????????。 答案： ?? ????????--??????????--=1556141501 4115401411A 2、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：2.367，0.25 3、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2 x 的系数为 ， 拉格朗日插值多项式为 。 答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 5、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )； 答案 )(1)(1n n n n n x f x f x x x '--- =+ 6、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 7、计算方法主要研究( 截断 )误差和( 舍入 )误差； 8、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为 ( 1 2+-n a b )； 9、求解一阶常微分方程初值问题y '= f (x ,y )，y (x 0)=y 0的改进的欧拉公式为 ( )] ,(),([2111+++++=n n n n n n y x f y x f h y y )； 吉林大学研究生公共数学课程 教学大纲 课程编号： 课程名称：现代数值计算方法 课程英文名称：Modern numerical method 学时/学分：64/3（硕士）/32/2（博士） 课程类别：研究生公共课程 课程性质：必修课 适用专业：理、工、经、管等专业 开课学期：第Ⅰ或第Ⅱ学期 考核方式：考试（闭卷） 执笔人：李永海 制定日期：2011年5月 吉林大学研究生公共数学课程教学大纲 课程编号： 课程名称：现代数值计算方法 课程英文名称：Modern numerical method 学时/学分：64/3（硕士）/32/2（博士） 课程类别：研究生教育课程 课程性质：必修课 适用专业：理、工、经、管等专业 开课学期：第Ⅰ或第Ⅱ学期 考核方式：考试（闭卷） 一、本课程的性质、目的和任务 本课程属于非数学类研究生数学公共基础课程之一，数值计算方法作为一种基本的数学工具，在数学学科与其他科学技术领域诸如力学、电磁学、化学、生物、系统工程等学科都有广泛应用。电子计算机及计算技术的发展也为数值计算方法的应用开辟了更广阔的前景。因此，学习和掌握现代数值计算方法，对于将来从事工程技术工作的工科研究生来说是必不可少的。通过该门课程的学习，期望学生能深刻地理解现代数值计算方法的基本知识和数学思想，掌握有关的计算方法及技巧，提高学生的数学素质，提高科研能力，掌握现代数值计算方法在物理、电子、化学、生物、工程等领域的许多应用。 二、本课程教学基本要求 1. 线性代数方程组直接法 理解线性代数方程组直接法求解算法原理，了解算法收敛性结果；理解算法应用条件；掌握用软件实现一般线性代数方程组直接法的求解步骤。 2. 线性代数方程组迭代法 理解线性代数方程组迭代法求解算法原理，了解算法收敛性结果；理解算法应用条件；掌握用软件实现一般线性代数方程组迭代法的求解步骤。 3. 矩阵特征值与特征向量计算 理解乘幂法和反幂法算法原理，了解实对称矩阵的Jacobi方法；理解算法应用条件；掌握用软件实现一般矩阵特征值与特征向量计算。 4. 非线性方程（组）求根 理解二分法和牛顿法原理，了解解非线性方程组的牛顿法和拟牛顿法；理解算法应用条件；掌握用软件实现非线性方程（组）求根计算。 5. 函数插值 理解一般函数插值公式原理，了解三次样条插值；理解算法应用条件；掌握用软件实现函数插值计算。 6. 数值积分 数值计算方法考试试题 一、选择题（每小题4分，共20分） 1. 误差根据来源可以分为四类，分别是（ A ） A. 模型误差、观测误差、方法误差、舍入误差； B. 模型误差、测量误差、方法误差、截断误差； C. 模型误差、实验误差、方法误差、截断误差； D. 模型误差、建模误差、截断误差、舍入误差。 2. 若132)(3 56++-=x x x x f ，则其六阶差商 =]3,,3,3,3[6210 f （ C ） A. 0； B. 1； C. 2； D. 3 。 3. 数值求积公式中的Simpson 公式的代数精度为 （ D ） A. 0； B. 1； C. 2； D. 3 。 4. 若线性方程组Ax = b 的系数矩阵A 为严格对角占优矩阵，则解方程组的Jacobi 迭代法和Gauss-Seidel 迭代法 （ B ） A. 都发散； B. 都收敛 C. Jacobi 迭代法收敛，Gauss-Seidel 迭代法发散； D. Jacobi 迭代法发散，Gauss-Seidel 迭代法收敛。 5. 对于试验方程y y λ='，Euler 方法的绝对稳定区间为（ C ） A. 02≤≤-h ； B. 0785.2≤≤-h ； C. 02≤≤-h λ； D. 0785.2≤≤-h λ ； 二、填空题（每空3分，共18分） 1. 已知 ? ??? ??--='-=4321,)2,1(A x ，则 =2 x 5，= 1Ax 16 ，=2A 22115+ 2. 已知 3)9(,2)4(==f f ，则 f (x )的线性插值多项式为)6(2.0)(1+=x x L ,且用线性插值可得f (7)= 2.6 。 3. 要使 20的近似值的相对误差界小于0.1%，应至少取 4 位有效数字。 三、利用下面数据表， 1. 用复化梯形公式计算积分 dx x f I )(6 .28 .1? =的近似值； 解：1.用复化梯形公式计算 取 2.048 .16.2,4=-= =h n 1分 分 分分7058337 .55))6.2()2.08.1(2)8.1((22.04)) ()(2)((231 1 1 4=+++=++=∑∑=-=f k f f b f x f a f h T k n k k 10.46675 8.03014 6.04241 4.42569 3.12014 f (x ) 2.6 2.4 2.2 2.0 1.8 x吉林大学 研究生 数值计算方法期末考试 样卷
数值计算方法试题集及答案
现代数值计算方法
《数值计算方法》试题及答案