计算方法考试试卷及答案
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
《计算方法》试卷(A 卷)
一、填空题(每空3分,共27分)
1、若15.3=x 是π的的近似值,则误差限是 0.05 ,有 2 位有效数字。
2、方程013
=--x x 在区间]2,1[根的牛顿迭代格式为
1
312131-)()(2
3
2
31
-+=---='-=+k k k k k k k k k k x x x x x x x f x f x x 。 3、对252)(23-+-=x x x f ,差商 =]3,3,3,3[432f -2 ,=]3,3,3,3,3[5
432f 0 。 4、数值积分中的梯形公式为
)]()([2
)(b f a f a
b dx x f b
a
+-≈
⎰
,Simpson 公式为 )]()2
(4)([6)(b f b
a f a f a
b dx x f b
a
+++-≈
⎰
。 5、求解微分方程初值问题⎩⎨
⎧==∈=5
.01
)0(]1,0['
h y x xy y 用欧拉公式计算得到=1y 1 ,用改进
的欧拉公式计算得到=1y 1.125 。
二、已知方程14
-=x x 在区间]2,0[内有根 (1)用二分法求该方程的根,要求误差不超过0.5。
(2)写出求解方程的一种收敛的简单迭代格式,并说明收敛原因。 解:(1)由题意,令分。3....,.........013)2(,01)0(,1)(4
<-=>=+-=f f x x x f 列表如下:
所以取1满足误差不超过0.5。...........................................7 分 (2) 原方程等价变形为4
1+=
x x ,迭代函数41)(+=x x ϕ,……………………….2分
则43
)1(41)(+=
'x x ϕ且在区间]2,0[上14
1
)1(41)(043
<<
+=
' <=≤+=≤=<ϕϕϕx x ,.7分 符合简单收敛的全局收敛条件, 所以收敛的简单迭代格式可构造为:3 15+=+k k x x .............................................8 分 三、利用x x f sin )(=在点2, 6, 0π π的函数值:(1)建立其拉格朗日插值多项式,并进行 误差分析;(2)构造差商表,建立牛顿插值多项式。 解:(1)由题意知,.1,2 ;21,6 ;0,0221100=== = ==y x y x y x π π ..........................2 分 所以过这三个点的拉格朗日插值函数为: )02 )(62()6)(0(1)26)(06() 2)(0(21)20)(60()2)(6(0) ()()()(2211002----⋅+----⋅+----⋅ =++=ππππ ππππππππx x x x x x x l y x l y x l y x L ........................ 5 分 x x x x x x π πππ ππ 27 3 )6 (6)2 (182 1 22 2222+ - =- + - - =.....................................7分 误差分析: )2 )(6)(0()!12()()()12(2π πξ---+=+x x x f x R .........................8 分 由题知:1|-cos ||)(| ,sin )() 3(≤=∴=ξξf x x f |)2 )(6(|61|)2)(6)(0()!12()(||)(|)12(2π πππξ--≤---+=+x x x x x x f x R ...........................10 分 (2)建立差商表: ....... ..... 4 分 牛顿插值多项式为 分 6)....)(](,,[)](,[)()(1021001002x x x x x x x f x x x x f x f x N --+-+= 分8 (27) 3)6)(0(3 )0(3 022 2x x x x x ππ πππ +-=--- -+ = 专业班级: 姓名: 学号: ……………密………………………………封………………………………线………………………… 四、确定以下求积公式中的系数,使其代数精度尽可能高,并确定代数精度。 )1()0()1()(1 1 Cf Bf Af dx x f ++-=⎰ - 解:令2 ,,1)(x x x f =,代入上面式子并令其等号两端相等得到:......2分 ⎪⎪⎩⎪⎪⎨⎧+==+-==++==⎰⎰⎰---1 1 21 11 132 021C A dx x C A xdx C B A dx .......6 分 解方程组可得:.31 ,34 ,31=== C B A ..............8 分 所以此求积公式为:)1(3 1 )0(34)1(31)(11f f f dx x f ++-=⎰-,...............10 分 将3 )(x x f =,代入上式可验证等式成立,但是将4 )(x x f =代入上式等式不成立, 所以此求积公式的代数精度为3. ......................12 分 五、求解矛盾方程组⎪⎪ ⎪⎪⎪⎭ ⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛3251011110 2121x x 的最小二乘解。 解,由题知,此矛盾方程组的系数矩阵分右端向量2............................. 3251b ,0112 1101⎪⎪⎪⎪ ⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=A 求得分4................... 6333 0112 1101011211 01⎪⎪⎭⎫ ⎝ ⎛=⎪⎪⎪⎪⎪ ⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛=A A T 分6................................ 96325101121101⎪⎪⎭ ⎫ ⎝⎛=⎪⎪⎪⎪⎪ ⎭ ⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛=b A T 求解分可得最小二乘解为分,10......................1,1 8..........2121===⎪⎪⎭ ⎫ ⎝⎛x x b A x x A A T T 六、(1)用高斯消去法求解方程组⎪⎩⎪ ⎨⎧=++-=+--=+-11 2123454321 321321x x x x x x x x x ; (2)对于方程组⎪⎩⎪ ⎨⎧=++=-+=--16 8374725 45321 321321x x x x x x x x x 用雅克比迭代格式及其对应的赛德尔迭代格式求解 是否收敛?并说明理由;写出两种迭代格式,取)0,0,0() 0(=x ,计算)1(x 。 解:(1) 高斯消去过程为: 分4...........550082104111191308210411111112123454111231312325-⎪⎪⎪⎭ ⎫ ⎝⎛----−→−⎪⎪⎪⎭⎫ ⎝⎛----−−→−⎪⎪⎪⎭⎫ ⎝⎛------r r r r r r 回代可得解为:分 6.........................3 ,6 ,1123==-=x x x (2) 由题此方程组的系数矩阵不是严格对角占优的,下用迭代矩阵的范数判断收敛性。 可知雅各比迭代矩阵为: ⎪⎪⎪⎪ ⎪⎪⎭ ⎫ ⎝⎛ -- =081 83740 72 -51540 J B ....................4 分,计算得140 37 1 <= J B .........................6 分 所以可以判断出解此方程的雅各比迭代格式及其相应的赛德尔迭代格式收敛。...............8 分 雅各比迭代格式为: 赛德尔迭代格式为: 分10.........)163(81)742(71)54(5 1) (2)(1)1(3 )(3)(1)1(2) (3)(2)1(1 ⎪⎪⎪ ⎩⎪⎪ ⎪⎨⎧+--=++-=++=+++k k k k k k k k k x x x x x x x x x 分14..........)163(81)742(71)54(5 1) 1(2)1(1)1(3 ) (3)1(1)1(2 ) (3)(2)1(1 ⎪⎪⎪ ⎩ ⎪⎪ ⎪⎨⎧+--=++-=++=++++++k k k k k k k k k x x x x x x x x x 分可得取)(12....)2,1,1(,)0,0,0()1(0T T x x == 分可得取) (16....)28 43,75,1(,)0,0,0()1(0T T x x == ………封………………………………线………………