东北大学10数值分析A(研)答案
- 格式:doc
- 大小:205.01 KB
- 文档页数:3
1
总分 1--3 4--6 7--9 10--12 13-15
班 级
学 号
姓 名
…………○…………密…………○…………封…………○………线…………………… 东 北 大 学 研 究 生 院 考 试 试 卷
2010 — 2011 学年第 一 学期
课程名称: 数值分析(A)
1.(5分)设近似值23.25x近似*x的相对误差限为0.0003,问x至少具有几位有效数字。
解 由于绝对误差限为:25.230.0003=0.007569<0.510-1 2分
所以,x至少具有3位有效数字。 5分
2.(6分)写出矩阵241353242A的Crout分解式TMA.
解 由于
1521613122241353242A 3分
所以,1006101211521013002A 6分
3.(7分)解线性方程组1764333222321321321xxxxxxxxx的Gauss-Seidel迭代法是否收敛,为什么?
解 令764332210)67)(23(6700123021 3分
得17/6)(G, 5分
所以,Gauss-Seidel迭代法收敛。 7分 4.(8分)用Jacobi法解线性方程组3422313321321321xxxxxxxxx,取000)0(x,估计迭代10步的误差*)10(xx。
解 由于04/12/13/10133/13/10B,所以4/3B, 2分
又由于,)4/3,3/2,3/1(1Tx所以4/3)0()1(xx 4分
所以,1689.04311011)0()1(10*)10(xxBBxx 8分
5.(9分)说明方程033xx在区间[1, 2]内有唯一根,并建立一个收敛的迭代格式,使对任意初值]2,1[0x都收敛,说明收敛理由。
解 由013)(,03)2(,03)1(2xxfff,知有唯一根。 3分
又由于253)(41333xx,13/1)3(313/2x 7分
所以,,...2,1,0,331kxxkk对任意初值]2,1[0x都收敛。 9分
6.(6分)设,...2,1,0,231kcbxaxxxkkkk是求方程根1的迭代法,试确定参数cba,,使迭代法的收敛阶尽可能高,并指出阶是多少?
解 由1得cba11 2分
令026)1(,023)1(aba得0,3,3cba, 5分
此时,迭代法3阶收敛。 6分 2
…………○…………密…………○…………封…………○………线…………………… 7.(6分)设324)(23xxxf,求差商]5,4,3,2,1[],4,3,2,1[],1,0[fff。
解 f[0,1]=(3-(-3))/1=6 2分
f[1,2,3,4]=4 4分
f[1,2,3,4,5]=0 6分
8.(7分)求满足条件0)1(,0)2(,2)1(,1)0(ffff的三次插值多项式)(3xH的表达式。
解 令 ))(2()(23cbxaxxxH 2分
则有:0,2,12cacbac 4分
解得:3,2/1bca 6分
所以,1213421)16)(2(21)(2323xxxxxxxH 7分
9.(7分)给定离散数据
xi -1 0 1 2
yi 2 -1 1 3
试求形如2bxay的拟合曲线。
解 由于基函数为:210)(,1)(xxx 1分
于是:)3,1,1,2(),4,1,0,1(),1,1,1,1(10f 3分
正则方程组为:15186564baba 4分
解得:6/5,0ba,
拟合曲线为: 265xy 7分 10.(5分)设求积公式)2(,)()(0nxfAdxxfkbankk是插值型求积公式,求nkkkxA02.
解 由于插值型求积公式代数精度至少是n, 2分
所以,)(3133220abdxxxAbaknkk 5分
11.(6分)对积分10)(dxxf建立两点Gauss公式。
解 由于 21),(),()(,1)(000010xPPPxPxxPxP,
61),(),(),(),()(2111210002022xxPPPxPPPPxPxxP 2分
Gauss点为:63211x,63212x 4分
积分公式为:)]6321()6321([21)(10ffdxxf 6分
12.(9分)利用复化Simpson公式2S计算定积分10cosxdxI的近似值,并估计误差。
解 ]43cos441cos421cos21cos0[cos1212SI 3分
=0.841489382 5分
1cosmax)(max)4(4xxfM 6分
000021701.0228801|)(|42fRSI 9分 3
…………○…………密…………○…………封…………○………线…………………… 13.(5分)求解初值问题2)1(21yxyeyx的改进Euler方法是否收敛?
为什么?
解 由于xyeyxf),(关于y满足Lipschitz条件, 2分
所以,改进Euler法收敛。 5分
14.(9分)已知求解常微分方程初值问题:
],[,)(),(baxayyxfy
的差分公式:
0121211)2,2(),()3(4yhkyhxfkyxfkkkhyynnnnnn
求此差分公式的阶。
解 由于
)()2(2)(2322222222hOfyffyxfxfhfyfxfhfknnnnnnnnn 2分
)(221nnnnnnfyfxfhhfyy)()2(242222223hOfyffyxfxfhnnnnn 4分
)()(6)(2)()()(4321hOxyhxyhhxyxyxynnnnn 5分
)(22nnnnnfyfxfhhfy
)(])(2[6422222223hOfyfyfxffyffyxfxfhhnnnnnnnn 7分
于是,)()(311hOyxynn,此差分公式是2阶的。 9分
15.(5分)证明矩阵谱半径)(A不是矩阵范数。
证明 因为0)(A时,不一定有0A, 例如0010A, 2分
所以,)(A不满足范数的非负性,不是范数。 5分