复习试题
一、填空题：
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 )；
10、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )；
11、 两点式高斯型求积公式?1
0d )(x x f ≈(?++-≈1
)]
321
3()3213([21d )(f f x x f )，代数精
度为( 5 )；
12、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为(A 的各阶顺序主子式均
不为零)。
13、 为了使计算
32)1(6
)1(41310--
-+-+
=x x x y 的乘除法次数尽量地少，应将该表
达式改写为
11
,))64(3(10-=
-++=x t t t t y ，为了减少舍入误差，应将表达式
19992001-
14、 用二分法求方程01)(3
=-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间
为 0.5，1 ,进行两步后根的所在区间为 0.5，0.75 。 15、 计算积分?1
5
.0d x
x ,取4位有效数字。用梯形公式计算求得的近似值为 0.4268 ，
用辛卜生公式计算求得的近似值为 0.4309 ，梯形公式的代数精度为 1 ，辛卜生公式的代数精度为 3 。
16、 求解方程组??
?=+=+042.01532121x x x x 的高斯—塞德尔迭代格式为
?????-=-=+++20/3/)51()1(1)1(2)(2)1(1k k k k x x x x ，该迭
代格式的迭代矩阵的谱半径)(M ρ= 121
。
17、 设46)2(,16)1(,0)0(===f f f ,则=)(1x l )2()(1--=x x x l ，)(x f 的二次牛顿
插值多项式为 )1(716)(2-+=x x x x N 。
18、 求积公式
?∑=≈b
a k n
k k x f A x x f )(d )(0
的代数精度以( 高斯型 )求积公式为最高，具
有( 12+n )次代数精度。
19、 已知f (1)=1,f (3)=5,f (5)=-3,用辛普生求积公式求?5
1
d )(x
x f ≈( 12 )。
20、 设f (1)=1， f (2)=2，f (3)=0，用三点式求≈')1(f ( 2.5 )。
21、如果用二分法求方程043
=-+x x 在区间]2,1[内的根精确到三位小数，需对分（ 10 ）次。
22、已知?????≤≤+-+-+-≤≤=31)1()1()1(2110)(2
33x c x b x a x x x x S 是三次样条函数，则
a =( 3 )，
b =（ 3 ），
c =（ 1 ）。
23、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数，则
∑==
n
k k
x l
)(( 1 )，∑==
n
k k j
k x l
x 0
)((
j
x )，当2≥n 时=
++∑=)()3(20
4
x l x x
k k n
k k ( 32
4++x x )。
24、解初值问题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 是
2 阶方法。
25、区间[]b a ,上的三次样条插值函数)(x S 在[]b a ,上具有直到_____2_____阶的连续导数。 26、改变函数f x x x ()=
+-1 (x >>1)的形式，使计算结果较精确
()x x x f ++=
11
。
27、若用二分法求方程()0=x f 在区间[1,2]内的根，要求精确到第3位小数，则需要对分 10
次。
28、设
()???≤≤+++≤≤=21,10,22
3
3x c bx ax x x x x S 是3次样条函数，则 a= 3 , b= -3 , c= 1 。 29、若用复化梯形公式计算
?
10
dx
e x ，要求误差不超过6
10-，利用余项公式估计，至少用 477
个求积节点。
30、写出求解方程组
??
?=+-=+2
4.016.12121x x x x 的Gauss-Seidel 迭代公式
()()
()() ,1,0,4.026.111112211=???+=-=+++k x x x x k k k k ，迭代矩阵为
????
??--64.006.10，此迭代法是否收敛 收敛 。
31、设
A =?? ?
??
5443,则=∞A 9 。
32、设矩阵
482257136A ????=??????的A LU =，则U = 4820161002U ??
????
=??
??-??
?? 。
33、若4
321()f x x x =++，则差商2481632[,,,,]f = 3 。
34、数值积分公式1
12
18019()[()()()]f x dx f f f -'≈-++?的代数精度为 2 。
35、
线性方程组121015112103x ????
????????=?????
???????的最小二乘解为
11??
?
?? 。
36、设矩阵
32120
4135A ??
??=??????分解为A LU =，则U = 321410033
21002??
????
??-???????
? 。 二、单项选择题：
1、 Jacobi 迭代法解方程组b x =A 的必要条件是（ C ）。 A ．A 的各阶顺序主子式不为零 B ． 1)( 2、设 ?? ????????--=700150322A ，则)(A ρ为( C )． A ． 2 B ． 5 C ． 7 D ． 3 3、三点的高斯求积公式的代数精度为( B )。 A ． 2 B ．5 C ． 3 D ． 4 4、求解线性方程组A x =b 的LU 分解法中，A 须满足的条件是( B )。 A ． 对称阵 B ． 正定矩阵 C ． 任意阵 D ． 各阶顺序主子式均不为零 5、舍入误差是( A )产生的误差。 A. 只取有限位数 B ．模型准确值与用数值方法求得的准确值 C ． 观察与测量 D ．数学模型准确值与实际值 6、3.141580是π的有( B )位有效数字的近似值。 A ． 6 B ． 5 C ． 4 D ． 7 7、用 1+x 近似表示e x 所产生的误差是( C )误差。 A ． 模型 B ． 观测 C ． 截断 D ． 舍入 8、解线性方程组的主元素消去法中选择主元的目的是( A )。 A ．控制舍入误差 B ． 减小方法误差 C ．防止计算时溢出 D ． 简化计算 9、用1+3x 近似表示3 1x +所产生的误差是( D )误差。 A ． 舍入 B ． 观测 C ． 模型 D ． 截断 10、-324．7500是舍入得到的近似值，它有( C )位有效数字。 A ． 5 B ． 6 C ． 7 D ． 8 11、设f (-1)=1,f (0)=3,f (2)=4,则抛物插值多项式中x 2的系数为( A )。 A ． –0．5 B ． 0．5 C ． 2 D ． -2 12、三点的高斯型求积公式的代数精度为( C )。 A ． 3 B ． 4 C ． 5 D ． 2 13、( D )的3位有效数字是0.236×102。 (A) 0.0023549×103 (B) 2354.82×10－2 (C) 235.418 (D) 235.54×10－1 14、用简单迭代法求方程f(x)=0的实根，把方程f(x)=0表示成x=?(x)，则f(x)=0的根是 ( B )。 (A) y=?(x)与x 轴交点的横坐标 (B) y=x 与y=?(x)交点的横坐标 (C) y=x 与x 轴的交点的横坐标 (D) y=x 与y=?(x)的交点 15、用列主元消去法解线性方程组??? ??-=+--=-+-=+-1 340921433 21321321x x x x x x x x x ，第1次消元，选择主元为 ( A ) 。 (A) －4 (B) 3 (C) 4 (D)－9 16、拉格朗日插值多项式的余项是( B ),牛顿插值多项式的余项是( C ) 。 (A) f(x,x0,x1,x2,…,xn)(x －x1)(x －x2)…(x －xn －1)(x －xn)， (B) )!1() ()()()()1(+= -=+n f x P x f x R n n n ξ (C) f(x,x0,x1,x2,…,xn)(x －x0)(x －x1)(x －x2)…(x －xn －1)(x －xn)， (D) ) ()!1() ()()()(1)1(x n f x P x f x R n n n n +++= -=ωξ 17、等距二点求导公式f '(x1) ≈( A )。 1011 0101 0010 101)()() D ()()() C ()()() B ()()() A (x x x f x f x x x f x f x x x f x f x x x f x f +--+---- 18、用牛顿切线法解方程f(x)=0，选初始值x0满足( A ),则它的解数列{xn}n=0,1,2,… 一定收敛到方程f(x)=0的根。 0)()()D (0 )()()C (0 )()()B (0 )()()A (0000<'<''>'>''x f x f x f x f x f x f x f x f 19、为求方程x3―x2―1=0在区间[1.3,1.6]内的一个根，把方程改写成下列形式，并建 立相应的迭代公式，迭代公式不收敛的是(A )。 (A) 1 1:,1 1 12-=-=+k k x x x x 迭代公式 (B) 21211:,11k k x x x x +=+ =+迭代公式 (C) 3/12123)1(:,1k k x x x x +=+=+迭代公式 (D) 11:,12 2 1 2 3+++==-+k k k k x x x x x x 迭代公式 20、求解初值问题?? ?=='0 0y x y y x f y )(),(欧拉法的局部截断误差是();改进欧拉法的局部截断误差 是();四阶龙格－库塔法的局部截断误差是( A ) (A)O(h2) (B)O(h3) (C)O(h4) (D)O(h5) 21、解方程组b Ax =的简单迭代格式g Bx x k k +=+) ()1(收敛的充要条件是（ ）。 （1）1)(A ρ, (4) 1)(>B ρ 22、在牛顿-柯特斯求积公式： ?∑=-≈b a n i i n i x f C a b dx x f 0 )() ()()(中，当系数) (n i C 是负值时，公式的 稳定性不能保证，所以实际应用中，当（ ）时的牛顿-柯特斯求积公式不使用。 （1）8≥n ， （2）7≥n ， （3）10≥n ， （4）6≥n ， （1）二次； （2）三次； （3）四次； （4）五次 24、若用二阶中点公式 )),(2,2(1n n n n n n y x f h y h x hf y y ++ +=+求解初值问题1)0(,2=-='y y y ， 试问为保证该公式绝对稳定，步长h 的取值范围为（ ）。 (1)10≤ 25 1732.≈ 计算4 1)x =，下列方法中哪种最好？（ ） (A)28-； (B)24(-； (C ) ； (D) 。 26、已知 33 0221224()()()x x S x x a x b x ?≤≤=?-+-+≤≤?是三次样条函数，则,a b 的值为( ) (A )6，6； (B)6，8； (C)8，6； (D)8，8。 (A)； (B)4； (C) ； (D ) 2。 28、形如112233()()()() b a f x dx A f x A f x A f x ≈++? 的高斯（Gauss ）型求积公式的代数精度为 （ ） (A)9； (B)7； (C ) 5； (D) 3。 29Newton 迭代格式为( ) (A) 132k k k x x x += +；(B )1322k k k x x x +=+；(C) 122k k k x x x +=+；(D) 133k k k x x x +=+。 30、用二分法求方程32 4100x x +-=在区间12[,]内的实根，要求误差限为31102ε-=?，则对分 次数至少为( ) (A )10； (B)12； (C)8； (D)9。 31、经典的四阶龙格—库塔公式的局部截断误差为 ( ) (A)4 ()O h ； (B)2 ()O h ； (C ) 5 ()O h ； (D) 3 ()O h 。 32、设()i l x 是以019(,,,)k x k k == 为节点的Lagrange 插值基函数，则9 ()i k kl k == ∑( ) (A)x ； （B ）k ； （C ）i ； （D ）1。 33、5个节点的牛顿-柯特斯求积公式，至少具有( )次代数精度 (A )5； (B)4； (C)6； (D)3。 34、已知 33 0221224()()()x x S x x a x b x ?≤≤=?-+-+≤≤?是三次样条函数，则,a b 的值为( ) (A )6，6； (B)6，8； (C)8，6； (D)8，8。 35、已知方程3 250x x --=在2x =附近有根，下列迭代格式中在02x =不收敛的是( ) (A)1k x += (B)1k x += (C )315k k k x x x +=--； (D) 3 1225 32k k k x x x ++=-。 (A ) 4； (B)2； (C)1； (D)3。 37、5个节点的Gauss 型求积公式的最高代数精度为( ) (A)8； (B )9； (C)10； (D)11。 四、计算题： 1、用高斯-塞德尔方法解方程组 ??? ??=++=++=++225218241124321321321x x x x x x x x x ，取T )0,0,0()0(=x ，迭代四次(要求按五位有效数字计算)。 答案：迭代格式 ??? ??? ???--=--=--=++++++)222(51)218(41)211(41) 1(2)1(1)1(3 ) (3)1(1)1(2 ) (3)(2)1(1k k k k k k k k k x x x x x x x x x 2、求A 、B 使求积公式 ?-+-++-≈1 1)]21 ()21([)]1()1([)(f f B f f A dx x f 的代数精度尽量高,并求其代数精度；利用此公式求 ? =2 1 1dx x I (保留四位小数)。 答案：2 ,,1)(x x x f =是精确成立，即 ??? ??=+=+32212222B A B A 得98,91==B A 求积公式为 )]21 ()21([98)]1()1([91)(1 1f f f f dx x f +-++-=?- 当3)(x x f =时，公式显然精确成立；当4 )(x x f =时，左=52，右=31。所以代 数精度为3。 69286.0140 97 ] 3 21132/11[98]311311[9131111322 1 ≈= +++-++++-≈+=??--=dt t dx x x t 3、已知 分别用拉格朗日插值法和牛顿插值法求)(x f 的三次插值多项式)(3x P ，并求)2(f 的近似值（保留四位小数）。 答案： )53)(43)(13() 5)(4)(1(6 )51)(41)(31()5)(4)(3(2 )(3------+------=x x x x x x x L )45)(35)(15() 4)(3)(1(4 )54)(34)(14()5)(3)(1(5 ------+------+x x x x x x 差商表为 )4)(3)(1(41 )3)(1()1(22)()(33---+ ----+==x x x x x x x N x P 5.5)2()2(3=≈P f 4、取步长2.0=h ，用预估-校正法解常微分方程初值问题 ?? ?=+='1)0(32y y x y )10(≤≤x 答案：解： ? ????+++?+=+?+=++++)]32()32[(1.0) 32(2.0)0(111)0(1n n n n n n n n n n y x y x y y y x y y 即 04.078.152.01++=+n n n y x y 5、已知 求)(x f 的二次拟合曲线)(2x p ，并求)0(f '的近似值。 答案：解： 正规方程组为 ??? ? ?=+==+41 34103101510520120a a a a a 1411,103,710210=== a a a 221411103710)(x x x p ++= x x p 711 103)(2+=' 103 )0()0(2 ='≈'p f 6、已知x sin 区间[0.4，0.8]的函数表 如用二次插值求63891.0sin 的近似值，如何选择节点才能使误差最小？并求该近似 值。 答案：解： 应选三个节点，使误差 |)(|!3|)(|33 2x M x R ω≤ 尽量小，即应使|)(|3x ω尽量小，最靠近插值点的三个节点满足上述要求。即取节点 }7.0,6.0,5.0{最好，实际计算结果 596274.063891.0sin ≈， 且 4 1055032.0)7.063891.0)(6.0963891.0)(5.063891.0(! 31 596274 .063891.0sin -?≤----≤ - 7、构造求解方程0210=-+x e x 的根的迭代格式 ,2,1,0),(1==+n x x n n ?，讨论其收敛 性，并将根求出来，4 110||-+<-n n x x 。 答案：解：令 010)1(, 02)0(, 210e )(>+=<-=-+=e f f x x f x . 且010e )(>+='x x f )(∞+-∞∈?， 对x ，故0)(=x f 在(0,1)内有唯一实根.将方程0)(=x f 变形为 )e 2(101 x x -= 则当)1,0(∈x 时 )e 2(101 )(x x -= ?， 1 10 e 10e |)(|<≤-='x x ? 故迭代格式 )e 2(101 1n x n x -= + 收敛。取5.00=x ，计算结果列表如下： 且满足 6671095000000.0||-<≤-x x .所以008525090.0*≈x . 8﹑利用矩阵的LU 分解法解方程组 ??? ??=++=++=++20 53182521432321321321x x x x x x x x x 。 答案：解： ?? ????????--??????????-==244132 11531 21LU A 令b y =L 得T )72,10,14(--=y ，y x =U 得T )3,2,1(=x . 9﹑对方程组 ??? ??=-+=--=++8 41025410151023321321321x x x x x x x x x （1） 试建立一种收敛的Seidel 迭代公式，说明理由； （2） 取初值T )0,0,0()0(=x ，利用（1）中建立的迭代公式求解，要求 3)()1(10||||-∞+<-k k x x 。 解：调整方程组的位置，使系数矩阵严格对角占优 ??? ??=++=-+=--15 1023841025410321321321x x x x x x x x x 故对应的高斯—塞德尔迭代法收敛.迭代格式为 ??? ??? ???+--=++-=++=++++++)1523(101)842(101)54(101)1(2)1(1)1(3) (3)1(1)1(2 )(3)(2)1(1k k k k k k k k k x x x x x x x x x 取T )0,0,0() 0(=x ,经7步迭代可得： T )010000.1,326950999.0,459991999.0()7(*=≈x x . 10、已知下列实验数据 试按最小二乘原理求一次多项式拟合以上数据。 解：当0 ，则 e )(≤''x f ，且x x d e 1 0?有一位整数. 要求近似值有5位有效数字，只须误差 4) (11021 )(-?≤ f R n . 由 )(12)()( 2 3 ) (1ξf n a b f R n ''-≤，只要 4 22) (1102112e 12e ) e (-?≤≤≤n n R x n ξ 即可，解得 ???=?≥ 30877.67106e 2n 所以 68=n ，因此至少需将 [0,1] 68等份。 11、用列主元素消元法求解方程组 ??????????--=???????????????? ????--11124112345111321x x x 。 解： ?? ? ?? ?????----???→???????? ????----111124111123451111212345411121r r ????? ? ?? ????????-----???→?????????? ??? ?????- -----???→?-5852510579515130123 4 5 57951513058525101234 5 5 2 51 321312r r r r r r ??? ?? ? ?? ??????? ?----???→?+135 1350579515130 123 4 5131 23r r 回代得 3,6,1123==-=x x x 。 12、取节点1,5.0,0210===x x x ,求函数x x f -=e )(在区间[0,1]上的二次插值多项式 )(2x P ,并估计误差。 解： )15.0)(05.0() 1)(0()10)(5.00()1)(5.0()(5.002----? +----? =--x x e x x e x P )5.0(2)1(4)1)(5.0(2) 5.01)(01() 5.0)(0(15.01-+----=----? +---x x e x x e x x x x e 又 1 |)(|max ,)(,)(] 1,0[3='''=-='''=∈--x f M e x f e x f x x x 故截断误差 |)1)(5.0(|!31 |)(||)(|22--≤ -=-x x x x P e x R x 。 13、用欧拉方法求 ?-=x t t x y 0 d e )(2 在点0.2,5.1,0.1,5.0=x 处的近似值。 解： ?-=x t t x y 0 d e )(2 等价于 ?????=='-0)0(e 2 y y x (0>x ) 记2 e ),(x y x f -=,取5.0=h ,0.2,5.1,0.1,5.0,043210=====x x x x x . 则由欧拉公式 ?? ?=+=+0),(01y y x hf y y n n n n , 3,2,1,0=n 可得 88940 .0)0.1(, 5.0)5.0(21≈==≈y y y y , 12604.1)0.2(,07334.1)5.1(43≈==≈y y y y 14、给定方程01e )1()(=--=x x x f 1) 分析该方程存在几个根； 2) 用迭代法求出这些根，精确到5位有效数字； 3) 说明所用的迭代格式是收敛的。 解：1）将方程 01e )1(=--x x （1） 改写为 x x -=-e 1 （2） 作函数1)(1-=x x f ，x x f -=e )(2的图形（略）知（2）有唯一根)2,1(*∈x 。 2) 将方程（2）改写为 x x -+=e 1 构造迭代格式 ?? ?=+=-+5.1e 101x x k x k ),2,1,0( =k 计算结果列表如下： 3) x x -+=e 1)(?，x x --='e )(? 当]2,1[∈x 时，]2,1[)]1(),2([)(?∈???x ，且 1e |)(|1<≤'-x ? 所以迭代格式 ),2,1,0()(1 ==+k x x k k ?对任意]2,1[0∈x 均收敛。 15、用牛顿(切线)法求3的近似值。取x 0=1.7, 计算三次，保留五位小数。 解：3是03)(2 =-=x x f 的正根，x x f 2)(='，牛顿迭代公式为 n n n n x x x x 232 1-- =+， 即 ) ,2,1,0(2321 =+=+n x x x n n n 取x 0=1.7, 列表如下： 16、已知f (-1)=2，f (1)=3，f (2)=-4，求拉格朗日插值多项式)(2x L 及f (1，5)的近似值，取五位小数。 解： )12)(12() 1)(1(4)21)(11()2)(1(3)21)(11()2)(1(2)(2-+-+? --+-+?+------? =x x x x x x x L )1)(1(34 )2)(1(23)2)(1(32-+--+---= x x x x x x 04167 .0241 )5.1()5.1(2≈=≈L f 17、n =3,用复合梯形公式求x x d e 10?的近似值（取四位小数）,并求误差估计。 解： 7342.1]e )e e (2e [3201d e 1210 310 ≈+++?-= ≈?T x x x x x f x f e )(,e )(=''=，10≤≤x 时，e |)(|≤''x f 05.0025.0108e 312e |e |||2 3≤==?≤ -= T R x 至少有两位有效数字。 18、用Gauss-Seidel 迭代法求解线性方程组 ????? ? ?--411131103????? ??321x x x =? ???? ??--815， 取x (0)=(0,0,0)T ,列表计算三次，保留三位小数。 解：Gauss-Seidel 迭代格式为： ??? ??? ???-+-=----=+-=++++++)8(41)1(31)5(31)1(2)1(1)1(3) (3)1(1)1(2 )(3)1(1k k k k k k k k x x x x x x x x 系数矩阵????? ? ????--411131103严格对角占优，故Gauss-Seidel 迭代收敛. 取x (0)=(0,0,0)T ，列表计算如下: 19、用预估—校正法求解?? ?=+='1)0(y y x y (0≤x ≤1)，h =0。2，取两位小数。 解：预估—校正公式为 ??? ?? ???? ++==++=+),(),()(2 1121211 k y h x hf k y x hf k k k y y n n n n n n ,2,1,0=n 其中y x y x f +=),(，10=y ，h =0.2，4,3,2,1,0=n ，代入上式得： 20、（82 bx a y +=解： },1{x span =Φ ??? ???=22 22 38312519 1111T A []3.730.493.320.19=T y 解方程组 y A AC A T T = 其中 ??????=3529603339133914A A T ??????=7.1799806.173y A T 解得： ??????=0501025.09255577.0C 所以 9255577.0=a ， 0501025.0=b 21、（15分）用8=n 的复化梯形公式（或复化 Simpson 公式）计算dx e x ? -1 时，试用余项估计其误 差。用8=n 的复化梯形公式（或复化 Simpson 公式）计算出该积分的近似值。 解： 001302.07681 81121)(12][022==??≤''-- =e f h a b f R T η ] )()(2)([2)8(7 1∑=++=k k b f x f a f h T ]36787947.0)41686207.047236655.05352614.060653066.07788008.08824969.0(21[16 1 ++++++?+= 6329434.0= 22、（15分）方程013 =--x x 在5.1=x 附近有根，把方程写成三种不同的等价形式（1）3 1 += x x 对应迭代格式311+=+n n x x ；(2)x x 11+=对应迭代格式n n x x 111+=+；（3）13-=x x 对应 迭代格式13 1-=+n n x x 。判断迭代格式在5.10=x 的收敛性，选一种收敛格式计算5.1=x 附近的根， 精确到小数点后第三位。 解：（1）32 1(31 )(-+='）x x ?，118.05.1<='）（?，故收敛； （2）x x x 1 121 )(2+ - ='?，117.05.1<='）（?，故收敛； （3）23)(x x ='?，15.135.12>?='）（?，故发散。 选择（1）：5.10=x ，3572.11=x ，3309.12=x ，3259.13=x ，3249.14=x ， 32476.15=x ，32472.16=x 23、（8分）已知方程组f AX =，其中 ??????????--=4114334A ，?? ??? ?????-=243024f （1） 列出Jacobi 迭代法和Gauss-Seidel 迭代法的分量形式。 （2） 求出Jacobi 迭代矩阵的谱半径。 解：Jacobi 迭代法：?? ???????=+-=+-=-=+++ ,3,2,1,0)24(41)330(41)324(41) (2)1(3)(3)(1)1(2) (2)1(1k x x x x x x x k k k k k k k Gauss-Seidel 迭代法：?? ???????=+-=+-=-=+++++ ,3,2,1,0)24(41)330(41)324(41) 1(2)1(3)(3)1(1)1(2) (2)1(1k x x x x x x x k k k k k k k ????? ? ????? ?--=+-=-0430430 430430)(1 U L D B J ， 790569 .0)4 10 (85)(==或J B ρ 24、1、（15分）取步长1.0=h ，求解初值问题?????=+-=1 )0(1y y dx dy 用改进的欧拉法求)1.0(y 的值；用经 典的四阶龙格—库塔法求)1.0(y 的值。 解：改进的欧拉法：??? ??+=++=+=+=++++095.0905.0)],(),([21.09.0),() 0(111) 0(1n n n n n n n n n n n n y y x f y x f h y y y y x hf y y 所以1)1.0(1==y y ； 经典的四阶龙格—库塔法： ?? ??? ?? ??? ? ++=++=++==++++=+),()2,2()2,2(),(]22[6 342312143211 hk y h x f k k h y h x f k k h y h x f k y x f k k k k k h y y n n n n n n n n n n 04321====k k k k ，所以1)1.0(1==y y 。 25、数值积分公式形如 ?'+'++=≈1 )1()0()1()0()()(f D f C Bf Af x S dx x xf 试确定参数D C B A ,,,使公式代数精度尽 量高；（2）设]1,0[)(4 C x f ∈，推导余项公式?-=1 ) ()()(x S dx x xf x R ，并估计误差。 解：将3 2 ,,,1)(x x x x f =分布代入公式得： 201,301,207,203-==== D B B A 构造Hermite 插值多项式)(3x H 满足?? ? ='='=1,0)()()()(33i x f x H x f x H i i i i 其中1,010==x x 则有：?=1 03)()(x S dx x xH ， 22)4(3)1(!4)()()(-=-x x f x H x f ξ dx x x f dx x S x f x x R 21 03 )4(1 0)1(!4)(])()([)(-=-=??ξ 1440)(60!4)()1(!4)()4()4(1023)4(ηηηf f dx x x f = ?=-=? 26、用二步法 )],()1(),([111101---+-+++=n n n n n n n y x f y x f h y y y θθαα 求解常微分方程的初值问题? ? ?=='00)() ,(y x y y x f y 时，如何选择参数θαα,,10使方法阶数尽可能高，并求局 部截断误差主项，此时该方法是几阶的 解： ] )(!3)(!2)()()(1()([) )(! 3)(!2)()(()()(!3)(!2)()()() 4(323 2103 211, +-'''+''-'-+'-+'''-''+'---+'''+''+'+=-=++n n n n n n n n n n n n n n n n h n x y h x y h x y h x y x y h x y h x y h x y h x y x y x y h x y h x y h x y y x y R θθαα ) ()()21661()()1221() ()11()()1(41312110h O x y h x y h x y h x y n n n n +'''--++''-+-+'+-+--=θαθαααα 所以?????? ?=-+-==--012210011110θαααα ???? ???===?230110θαα 主项：) (1253 n x y h ''' 该方法是二阶的。 27、（10分）已知数值积分公式为： )]()0([)]()0([2)(''20 h f f h h f f h dx x f h -++≈ ? λ，试确定积分公式中的参数λ，使其代数精 确度尽量高，并指出其代数精确度的次数。 解：1)(=x f 显然精确成立； x x f =)(时， ] 11[]0[22220 -++==? h h h h xdx h λ； 2)(x x f =时，12122]20[]0[2332 230 2 = ?-=-++==?λλλh h h h h h h dx x h ； 3)(x x f =时，] 30[121 ]0[24223403 h h h h h dx x h -++==?； 4)(x x f =时，6]40[121]0[2553 2450 4 h h h h h h dx x h = -++≠=?； 所以，其代数精确度为3。 28、（8分）已知求)0(>a a 的迭代公式为： 2,1,00)(2101=>+= +k x x a x x k k k 证明：对一切a x k k ≥=,,2,1 ，且序列{}k x 是单调递减的， 从而迭代过程收敛。 证明： 2,1,0221)(211==???≥+= +k a x a x x a x x k k k k k 《数值计算方法》复习试题 一、填空题： 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 。 答案：， 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 )； 10、已知f (1)＝2，f (2)＝3，f (4)＝，则二次Newton 插值多项式中x 2系数为( )； 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f )，代数精 度为( 5 )； 12、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为(A 的各阶顺序主子式均 不为零)。 13、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该表 达式改写为 11 ,))64(3(10-= -++=x t t t t y ，为了减少舍入误差，应将表达式 19992001-改写为 199920012 + 。 14、 用二分法求方程01)(3 =-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间 为 ，1 ,进行两步后根的所在区间为 ， 。 15、 、 16、 计算积分?1 5 .0d x x ,取4位有效数字。用梯形公式计算求得的近似值为 ，用辛卜 生公式计算求得的近似值为 ，梯形公式的代数精度为 1 ，辛卜生公式的代数精度为 3 。 17、 求解方程组?? ?=+=+042.01532121x x x x 的高斯—塞德尔迭代格式为 ?????-=-=+++20/3/)51()1(1)1(2)(2)1(1 k k k k x x x x ，该迭 代格式的迭代矩阵的谱半径)(M ρ= 121 。 18、 设46)2(,16)1(,0)0(===f f f ,则=)(1x l )2()(1--=x x x l ，)(x f 的二次牛顿 插值多项式为 )1(716)(2-+=x x x x N 。 19、 求积公式 ?∑=≈b a k n k k x f A x x f )(d )(0 的代数精度以( 高斯型 )求积公式为最高，具 有( 12+n )次代数精度。 目录 第一章非线性方程求根 (3) 1.1迭代法 (3) 1.2牛顿法 (4) 1.3弦截法 (5) 1.4二分法 (6) 第二章插值 (7) 2.1线性插值 (7) 2.2二次插值 (8) 2.3拉格朗日插值 (9) 2.4分段线性插值 (10) 2.5分段二次插值 (11) 第三章数值积分 (13) 3.1复化矩形积分法 (13) 3.2复化梯形积分法 (14) 3.3辛普森积分法 (15) 3.4变步长梯形积分法 (16) 第四章线性方程组数值法 (17) 4.1约当消去法 (17) 4.2高斯消去法 (18) 4.3三角分解法 (20) 4.4雅可比迭代法 (21) 4.5高斯—赛德尔迭代法 (23) 第五章常积分方程数值法 (25) 5.1显示欧拉公式法 (25) 5.2欧拉公式预测校正法 (26) 5.3改进欧拉公式法 (27) 5.4四阶龙格—库塔法 (28) 数值计算方法 第一章非线性方程求根 1.1迭代法 程序代码： Private Sub Command1_Click() x0 = Val(InputBox("请输入初始值x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = (Exp(2 * x0) - x0) / 5 If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例：求f(x)=e2x-6x=0在x=0.5附近的根（ep=10-10） 1.2牛顿法 程序代码： Private Sub Command1_Click() b = Val(InputBox("请输入被开方数x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = x0 - (x0 ^ 2 - b) / (2 * b) If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例：求56的值。（ep=10-10） 《计算方法》期中复习试题 一、填空题： 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- 六年级数学简便计算专项练习题（附答案＋计算方法汇总） 小学阶段（高年级）的简便运算，在一定程度上突破了算式原来的运算顺序，根据运算定律、性质重组运算顺序。如果学生没真正理解运算定律、性质，他只能照葫芦画瓢。在实际解题的过程当中，学生的思路不清晰，常出现这样或那样的错误。因此，培养学生思维的灵活性就显得尤为重要。 下面，为大家整理了8种简便运算的方法，希望同学们在理解的基础上灵活运用，不提倡死记硬背哟！ 1.提取公因式 这个方法实际上是运用了乘法分配律，将相同因数提取出来，考试中往往剩下的项相加减，会出现一个整数。 注意相同因数的提取。 例如： 0.92×1.41＋0.92×8.59 =0.92×（1.41+8.59） 2.借来借去法 看到名字，就知道这个方法的含义。用此方法时，需要注意观察，发现规律。还要注意还哦,有借有还，再借不难。 考试中，看到有类似998、999或者1.98等接近一个非常好计算的整数的时候，往往使用借来借去法。 例如： 9999+999+99+9 =9999+1+999+1+99+1+9+1－4 3.拆分法 顾名思义，拆分法就是为了方便计算把一个数拆成几个数。这需要掌握一些“好朋友”，如：2和5，4和5，2和2.5，4和2.5，8和1.25等。分拆还要注意不要改变数的大小哦。 例如： 3.2×12.5×25 =8×0.4×12.5×25 =8×12.5×0.4×25 4.加法结合律 注意对加法结合律 (a＋b)＋c=a＋(b＋c) 的运用，通过改变加数的位置来获得更简便的运算。 例如： 5.76＋13.67＋4.24＋ 6.33 =（5.76＋4.24）＋(13.67＋6.33) 5.拆分法和乘法分配律结合 这种方法要灵活掌握拆分法和乘法分配律，在考卷上看到99、101、9.8等接近一个整数的时候，要首先考虑拆分。 例如： 34×9.9 = 34×(10－0.1) 案例再现：57×101=？ 6.利用基准数 在一系列数种找出一个比较折中的数字来代表这一系列的数字，当然要记得这个数字的选取不能偏离这一系列数字太远。 例如： 2072+2052+2062+2042+2083 参考答案 第一章 1 *1x =1.7； * 2x =1.73； *3x =1.732 。 2． 3． （1） ≤++)(* 3*2*1x x x e r 0.00050； （注意：应该用相对误差的定义去求） （2） ≤)(*3*2*1x x x e r 0.50517； （3） ≤)/(*4*2x x e r 0.50002。 4．设6有n 位有效数字，由6≈2.4494……，知6的第一位有效数字1a =2。 令3)1()1(1* 102 1 102211021)(-----?≤??=?= n n r a x ε 可求得满足上述不等式的最小正整数n =4，即至少取四位有效数字，故满足精度要求可取6≈2.449。 5. 答：（1）*x (0>x )的相对误差约是* x 的相对误差的1/2倍； （2）n x )(* 的相对误差约是* x 的相对误差的n 倍。 6． 根据******************** sin 21)(cos 21sin 21)(sin 21sin 21)(sin 21)(c b a c e c b a c b a b e c a c b a a e c b S e r ++≤ =* *****) ()()(tgc c e b b e a a e ++ 注意当20* π < 7．设20= y ，41.1*0 =y ，δ=?≤--2* 00102 1y y 由 δ1* 001*111010--≤-=-y y y y ， δ2*111*221010--≤-=-y y y y M δ10*991*10101010--≤-=-y y y y 即当0y 有初始误差δ时，10y 的绝对误差的绝对值将减小10 10-倍。而110 10 <<-δ，故计算过程稳定。 8. 变形后的表达式为： （1）)1ln(2--x x =)1ln(2-+-x x （2）arctgx x arctg -+)1(=) 1(11 ++x x arctg （3） 1ln )1ln()1(ln 1 --++=? +N N N N dx x N N =ΛΛ+-+- +3 2413121)1ln(N N N N 1ln )11ln()1(-++ +=N N N N =1)1ln()1 1ln(-+++N N N （4）x x sin cos 1-=x x cos 1sin +=2x tg 习题一 1. 什么叫数值方法？数值方法的基本思想及其优劣的评价标准如何？ 数值方法是利用计算机求解数学问题近似解的方法 x max x i , x ( x 1 , x 2 , x n ) T R n 及 A n R n n . 2. 试证明 max a ij , A ( a ij ) 1 i n 1 i n 1 j 证明： （ 1）令 x r max x i 1 i n n p 1/ p n x i p 1/ p n x r p 1/ p 1/ p x lim( x i lim x r [ ( ] lim x r [ lim x r ) ) ( ) ] x r n p i 1 p i 1 x r p i 1 x r p 即 x x r n p 1/ p n p 1/ p 又 lim( lim( x r x i ) x r ) p i 1 p i 1 即 x x r x x r ⑵ 设 x (x 1,... x n ) 0 ，不妨设 A 0 ， n n n n 令 max a ij Ax max a ij x j max a ij x j max x i max a ij x 1 i n j 1 1 i n j 1 1 i n j 1 1 i n 1 i n j 1 即对任意非零 x R n ，有 Ax x 下面证明存在向量 x 0 0 ，使得 Ax 0 ， x 0 n ( x 1,... x n )T 。其中 x j 设 j a i 0 j ，取向量 x 0 sign(a i 0 j )( j 1,2,..., n) 。 1 n n 显然 x 0 1 且 Ax 0 任意分量为 a i 0 j x j a i 0 j ， i 1 i 1 n n 故有 Ax 0 max a ij x j a i 0 j 即证。 i i 1 j 1 3. 古代数学家祖冲之曾以 355 作为圆周率的近似值，问此近似值具有多少位有效数字？ 113 解： x 325 &0.314159292 101 133 x x 355 0.266 10 6 0.5 101 7 该近似值具有 7 为有效数字。 2020年奥鹏吉大网络教育《计算方法》大作业解答 (说明：前面是题目，后面几页是答案完整解答部分，注意的顺序。) 一、解线性方程 用矩阵的LU分解算法求解线性方程组 用矩阵的Doolittle分解算法求解线性方程组 用矩阵的Doolittle分解算法求解线性方程组 用高斯消去法求解线性方程组 用高斯消去法求解线性方程组 用主元素消元法求解线性方程组 用高斯消去法求解线性方程组 利用Doolittle分解法解方程组Ax=b，即解方程组 1、用矩阵的LU分解算法求解线性方程组 X1+2X2+3X3 = 0 2X1+2X2+8X3 = -4 -3X1-10X2-2X3 = -11 2、用矩阵的Doolittle分解算法求解线性方程组 X1+2X2+3X3 = 1 2X1– X2+9X3 = 0 -3X1+ 4X2+9X3 = 1 3、用矩阵的Doolittle分解算法求解线性方程组 2X1+X2+X3 = 4 6X1+4X2+5X3 =15 4X1+3X2+6X3 = 13 4、用高斯消去法求解线性方程组 2X 1- X 2+3X 3 = 2 4X 1+2X 2+5X 3 = 4 -3X 1+4X 2-3X 3 = -3 5、用无回代过程消元法求解线性方程组 2X 1- X 2+3X 3 = 2 4X 1+2X 2+5X 3 = 4 -3X 1+4X 2-3X 3 = -3 6、用主元素消元法求解线性方程组 2X 1- X 2+3X 3 = 2 4X 1+2X 2+5X 3 = 4 -3X 1+4X 2-3X 3 = -3 7、用高斯消去法求解线性方程组 123123123234 4272266 x x x x x x x x x -+=++=-++= 8、利用Doolittle 分解法解方程组Ax=b ，即解方程组 12341231521917334319174262113x x x x -? ????? ???? ??-??????=? ? ????--?????? --???? ?? 《数值计算方法》复习试题 一、填空题: 1、????? ?????----=410141014A ，则A 的ＬU 分解为 A ??? ?????????=? ?????????? ?。 答案： ?? ????????--??????????--=1556141501 4115401411A 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 )位有效数字； ５、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( ）； 答案 )(1)(1n n n n n x f x f x x x '--- =+ ６、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 )，=]4,3,2,1,0[f ( ０ )； 7、计算方法主要研究（ 截断 ）误差和( 舍入 )误差; 8、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时,二分n 次后的误差限为 （ 1 2+-n a b ); １０、已知f （１）＝2,f (2)=３,ｆ(４)＝5.９，则二次Ne ｗton 插值多项式中x 2系数为 ( ０.15 )； 11、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为（A 的各阶顺序主子式均 不为零)。 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该 计算方法上机报告 姓名： 学号： 班级： 目录 题目一------------------------------------------------------------------------------------------ - 4 - 1.1题目内容 ---------------------------------------------------------------------------- - 4 - 1.2算法思想 ---------------------------------------------------------------------------- - 4 - 1.3Matlab源程序----------------------------------------------------------------------- - 5 - 1.4计算结果及总结 ------------------------------------------------------------------- - 5 - 题目二------------------------------------------------------------------------------------------ - 7 - 2.1题目内容 ---------------------------------------------------------------------------- - 7 - 2.2算法思想 ---------------------------------------------------------------------------- - 7 - 2.3 Matlab源程序---------------------------------------------------------------------- - 8 - 2.4计算结果及总结 ------------------------------------------------------------------- - 9 - 题目三----------------------------------------------------------------------------------------- - 11 - 3.1题目内容 --------------------------------------------------------------------------- - 11 - 3.2算法思想 --------------------------------------------------------------------------- - 11 - 3.3Matlab源程序---------------------------------------------------------------------- - 13 - 3.4计算结果及总结 ------------------------------------------------------------------ - 14 - 题目四----------------------------------------------------------------------------------------- - 15 - 4.1题目内容 --------------------------------------------------------------------------- - 15 - 4.2算法思想 --------------------------------------------------------------------------- - 15 - 4.3Matlab源程序---------------------------------------------------------------------- - 15 - 4.4计算结果及总结 ------------------------------------------------------------------ - 16 - 题目五----------------------------------------------------------------------------------------- - 18 - 计算方法实验报告 班级： 学号： 姓名： 成绩： 1 舍入误差及稳定性 一、实验目的 （1）通过上机编程，复习巩固以前所学程序设计语言及上机操作指令； （2）通过上机计算，了解舍入误差所引起的数值不稳定性 二、实验内容 1、用两种不同的顺序计算10000 21n n -=∑，分析其误差的变化 2、已知连分数() 1 01223//(.../)n n a f b b a b a a b =+ +++，利用下面的算法计算f ： 1 1 ,i n n i i i a d b d b d ++==+ (1,2,...,0 i n n =-- 0f d = 写一程序，读入011,,,...,,,...,,n n n b b b a a 计算并打印f 3、给出一个有效的算法和一个无效的算法计算积分 1 041 n n x y dx x =+? (0,1,...,1 n = 4、设2 2 11N N j S j == -∑ ，已知其精确值为1311221N N ?? -- ?+?? （1）编制按从大到小的顺序计算N S 的程序 （2）编制按从小到大的顺序计算N S 的程序 （3）按两种顺序分别计算10001000030000,,,S S S 并指出有效位数 三、实验步骤、程序设计、实验结果及分析 1、用两种不同的顺序计算10000 2 1n n -=∑，分析其误差的变化 （1）实验步骤： 分别从1~10000和从10000~1两种顺序进行计算，应包含的头文件有stdio.h 和math.h （2）程序设计： a.顺序计算 #include 《现代设计方法》课程习题集 西南科技大学成人、网络教育学院 版权所有 习题 【说明】：本课程《现代设计方法》（编号为09021）共有单选题,计算题,简答题, 填空题等多种试题类型，其中，本习题集中有[ 填空题,单选题]等试题类型未进入。 一、计算题 1. 用黄金分割法求解以下问题（缩小区间三次）。 342)(m in 2+-=x x x f ，给定初始区间[][]3,0,=b a ，取1.0=ε。 2. 用黄金分割法求解以下问题（缩小区间三次） 32)(m in 2+=x x f ，给定[][],1,2a b =-，取1.0=ε 3. 用黄金分割法求解以下问题（缩小区间三次） 432+=x )x (f min ，给定[][]40,b ,a =，取10.=ε。 4. 用黄金分割法求解以下问题（缩小区间三次）。 12)(m in 3+-=x x x f ，给定初始区间[][]3,0,=b a ，取5.0=ε 5. 用黄金分割法求解以下问题（缩小区间三次）。 107)(m in 2+-=x x x f ，给定初始区间[][]3,0,=b a ，取1.0=ε 6. 用梯度法求解无约束优化问题： 168)(m in 22221+-+=x x x X f ，取初始点[]T X 1,1)0(= ，计算精度1.0=ε。 7. 用梯度法求解96)(m in 12221+-+=x x x X f ，[]T X 1,1)0(= ，1.0=ε。 8. 用梯度法求解44)(m in 22221+-+=x x x X f ，[]T X 1,1)0(=，1.0=ε 。 9. 用梯度法求解无约束优化问题：1364)(m in 222 121+-+-=x x x x X f ，取初始点 : 第一次作业 1.下列各数都是经过四舍五入得到的近似数，指出他们有几位有效数字，并写出绝对误差限。 9800107480.566.385031.01021.1*65*5*4*3*2*1=?=====x x x x x x 解： 1* 11011021.01021.1?==x ，有5位有效数字，绝对误差限为4-5-1105.0105.0?=?； 1-* 2 1031.0031.0?==x ，有2位有效数字，绝对误差限为3-2-1-105.0105.0?=?； 3* 3103856.06.385?==x ；有4位有效数字，绝对误差限为-14-3105.0105.0?=?； 2* 41056480.0480.56?==x ；有5位有效数字，绝对误差限为3-5-2105.0105.0?=?； ; 65* 5 107.0107?=?=x ；有1位有效数字，绝对误差限为51-6105.0105.0?=?； 4* 6 109800.09800?==x ；有4位有效数字，绝对误差限为5.0105.04-4=?。 2.要使20的近似值的相对误差限小于%1.0，要取几位有效数字 解：由于110447213595.047213595.420??=?=，设要取n 位有效数字，则根据 定理，有()()%1.01081 1021111 计算方法模拟试题 一、 单项选择题（每小题3分，共15分） 1．近似值210450.0?的误差限为（ ）。 A ． 0.5 B. 0.05 C ． 0.005 D. 0.0005. 2. 求积公式)2(3 1 )1(34)0(31)(2 0f f f dx x f ++≈ ?的代数精确度为（ ）。 A. 1 B. 2 C. 3 D. 4 3. 若实方阵A 满足（ ）时，则存在唯一单位下三角阵L 和上三角阵R ，使LR A =。 A. 0det ≠A B. 某个0 det ≠k A C. )1,1(0det -=≠n k A k D. ),,1(0det n k A k =≠ 4.已知?? ?? ? ?????=531221112A ，则=∞A （ ）。 A. 4 B. 5 C. 6 D 9 5.当实方阵A 满足)2(,221>>-=i i λλλλ，则乘幂法计算公式1e =（ ）。 A. 1+k x B. k k x x 11λ++ C. k x D. k k x x 11λ-+ 二、填空题（每小题3分，共15分） 1. 14159.3=π，具有4位有效数字的近似值为 。 2. 已知近似值21,x x ，则=-?)(21x x 。 3．已知1)(2-=x x f ，则差商=]3,2,1[f 。 4．雅可比法是求实对称阵 的一种变换方法。《数值计算方法》试题集及答案
数值计算方法大作业
数值计算方法》试题集及答案
六年级数学简便计算专项练习题(附答案+计算方法汇总)
西工大计算方法作业答案
(完整word版)计算方法习题集及答案.doc
2020年奥鹏吉大网络教育《计算方法》大作业解答
《数值计算方法》试题集及答案
西安交通大学计算方法B大作业
计算方法上机实习题大作业(实验报告).
现代设计方法复习题集含答案
计算方法作业参考答案(不断更新)
计算方法模拟试题及答案