东北大学10数值分析A(研)答案

  • 格式:doc
  • 大小:205.01 KB
  • 文档页数:3

1

总分 1--3 4--6 7--9 10--12 13-15

班 级

学 号

姓 名

…………○…………密…………○…………封…………○………线…………………… 东 北 大 学 研 究 生 院 考 试 试 卷

2010 — 2011 学年第 一 学期

课程名称: 数值分析(A)

1.(5分)设近似值23.25x近似*x的相对误差限为0.0003,问x至少具有几位有效数字。

解 由于绝对误差限为:25.230.0003=0.007569<0.510-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/3B, 2分

又由于,)4/3,3/2,3/1(1Tx所以4/3)0()1(xx 4分

所以,1689.04311011)0()1(10*)10(xxBBxx 8分

5.(9分)说明方程033xx在区间[1, 2]内有唯一根,并建立一个收敛的迭代格式,使对任意初值]2,1[0x都收敛,说明收敛理由。

解 由013)(,03)2(,03)1(2xxfff,知有唯一根。 3分

又由于253)(41333xx,13/1)3(313/2x 7分

所以,,...2,1,0,331kxxkk对任意初值]2,1[0x都收敛。 9分

6.(6分)设,...2,1,0,231kcbxaxxxkkkk是求方程根1的迭代法,试确定参数cba,,使迭代法的收敛阶尽可能高,并指出阶是多少?

解 由1得cba11 2分

令026)1(,023)1(aba得0,3,3cba, 5分

此时,迭代法3阶收敛。 6分 2

…………○…………密…………○…………封…………○………线…………………… 7.(6分)设324)(23xxxf,求差商]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,12cacbac 4分

解得:3,2/1bca 6分

所以,1213421)16)(2(21)(2323xxxxxxxH 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(10f 3分

正则方程组为:15186564baba 4分

解得:6/5,0ba,

拟合曲线为: 265xy 7分 10.(5分)设求积公式)2(,)()(0nxfAdxxfkbankk是插值型求积公式,求nkkkxA02.

解 由于插值型求积公式代数精度至少是n, 2分

所以,)(3133220abdxxxAbaknkk 5分

11.(6分)对积分10)(dxxf建立两点Gauss公式。

解 由于 21),(),()(,1)(000010xPPPxPxxPxP,

61),(),(),(),()(2111210002022xxPPPxPPPPxPxxP 2分

Gauss点为:63211x,63212x 4分

积分公式为:)]6321()6321([21)(10ffdxxf 6分

12.(9分)利用复化Simpson公式2S计算定积分10cosxdxI的近似值,并估计误差。

解 ]43cos441cos421cos21cos0[cos1212SI 3分

=0.841489382 5分

1cosmax)(max)4(4xxfM 6分

000021701.0228801|)(|42fRSI 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时,不一定有0A, 例如0010A, 2分

所以,)(A不满足范数的非负性,不是范数。 5分