重庆大学研究生数值分析试题解析
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
(1)写出Jacobi法和SOR法的迭代格式(分量形式);
(2)讨论这两种迭代法的收敛性.
(3)取初值x(0)=(0,0,0)T,若用Jacobi迭代法计算时,
预估误差x*-x(10) (取三位有效数字).
解 (1)Jacobi法和SOR法的迭代格式分别为
x
( 1
k
1
)
x
解 只要取(x)=x3-a ,或(x)=1-x3/a. 5.设(x)=x3+x2-3,则差商[3,32,33,34]=__1_____.
6.设l0(x),l1(x),l2(x),l3(x)是以x0,x1,x2,x3为互异节点
的三次插值基函数,则
3
l j
(x)(xj
2)3=___(_x_-_2_)_3____.
72
2 1
,所以‖A‖1=5,‖A-1‖1=5/7.
2.设矩阵A=
1 a
a 1
a 0
,当a取______值时,A可以唯一分解
a 0 1
为GGT,其中G为下三角矩阵.
1aa
解
令
1 a
a1a20,a 1
a
1 0
012a20, 得:1 a 1
1
2
2
3.向量x=(x1,x2,x3)T,试问|x1|+|2x2|+|x3|是不是一种向 量范数__是____,而|x1|+|2x2+x3|是不是一种向量范数_不_是___.
4.求 3 a 的Newton迭代格式为_x _k _1 __xk _ _x _3 k 3 _x k _2a _或 _x _k_ 1_ _3 2 _x _k_ _a 3 _x_k 2 _.
于是 H3(x)=-(x-1)2(x-2)-3x(x-2)+2.5x(x-1)2 –0.5x(x-1)(x-2) =x3-2.5x2 +2.5x+2
由于,R(0)=R(1)=R(2)=R(1)=0, 故可设 R(x)=C(x)x(x-1)2(x-2)
构造函数(t)=(t)-H3(t)-C(x)t(t-1)2(t-2) 于是,存在x,使(4)(x)=0,即(4)(x)-4!C(x)=0
x2(k1)
x(k) 2
(1
5
x(k1) 1
x(k) 2
1 5
x(k) 3
2) 5
x3(k1)
x(k) 3
(1
3
x(k1) 1
1 6
x(k1) 2
x(k) 3
1) 2
(2)因为A是严格对角占优矩阵,但不是正定矩阵,故
Jacobi法收敛,SOR法当0<1时收敛.
1 0 .75
四、(13分)已知(0)=2,(1)=3,(2)=5,(1)=0.5, (1)试建立一个三次插值多项式H3(x),使满足插值条件:
H3(0)=2,H3(1)=3,H3(2)=5,H3(1)=0.5; (2)设y=(x)在[0,2]上四次连续可微,试确定插值余项
(3)由(1)可见B=3/4,且取x(0)=(0,0,0)T,经计算可
得x(1)=(1/4,-2/5,1/2)T,于是x(1)-x(0)=1/2,所以有
B k x * x (1)0 x (1 ) x (0 )
0 .715 0 0 .5 0 .113
1 B
R(x)=(x)-H3(x).
解 (1)由y0=2,y1=3,y2=5,y1=0.5,得 H3(x)=20(x)+31(x)+52(x)+0.51(x) 令0(x)=c(x-1)2(x-2),可得0(x)=-0.5(x-1)2(x-2), 令1(x)=x(x-2)(ax+b),可得1(x)=-x(x-2), 令2(x)=cx(x-1)2,可得2(x)=0.5x(x-1)2; 令1(x)=cx(x-1)(x-2),可得1(x)=-x(x-1)(x-2),
取初值x0=1.5, 计算得x1=1.41333, x2=1.40983,由于 |x2-x1|=0.0035<10-2 , 故可取根的近似值x2=1.40983.
(3)因为0<</2,所以()co /2 s1sin 0
故,此迭代法线性收敛(收敛阶为1).
三、(14分)设线性方程组 4x1 x2 2x3 1 x1 5x2 x3 2 2x1 x2 6x3 3
j0
7.设S(x)=
x3x2 2x3b2xcx1
0x1 是以0,1,2为节
1x2
点的三次样条函数,则b=__-__2____c=___3______. 解 由2=b+c+1,5=6+2b+c,8=12+2b,可得
二、(13分)设函数(x)=x2-sinx-1 (1)试证方程(x)=0有唯一正根; (2)构造一种收敛的迭代格式xk=(xk),k=0,1,2,…计
( 2
k
1
)
1 4
x (k) 2
1 2
x (k) 3
1 4
1 5
x (k) 1
1 5
x (k) 3
2 5
x
( 3
k
1
)
1 3
x (k) 1
1 6
x (k) 2
1 2
x1(k1)
x(k) 1
(x1(k)
1 4
x(k) 2
1 2
x(k) 3
1) 4
考试题解析
一、填空题(每空3分,共30分)
1.设矩阵A= 1
2
2 3
Байду номын сангаас
,则(A)=____7 ___,Cond(A)1=___2_5/_7__.
解 由于
AE1
2
2 2 4 7 0
3
得特征值: 1 2 3 i,2 2 3 i
又A-1= 1 3
算精度为=10-2的近似根; (3)此迭代法的收敛阶是多少?说明之. 解 (1)因为0<x1时,(x)<0,x2时,(x)>0,所以(x)
仅在(1,2)内有零点,而当1<x<2时,(x)>0,故(x)单调. 因此方程(x)=0有唯一正根,且在区间(1,2)内.
(2)构造迭代格式: x k 1 1 sx ikn k 0 ,1 ,2 ,... 由于|(x)|=| co x/s 21sixn |<1,故此迭代法收敛.
(2)讨论这两种迭代法的收敛性.
(3)取初值x(0)=(0,0,0)T,若用Jacobi迭代法计算时,
预估误差x*-x(10) (取三位有效数字).
解 (1)Jacobi法和SOR法的迭代格式分别为
x
( 1
k
1
)
x
解 只要取(x)=x3-a ,或(x)=1-x3/a. 5.设(x)=x3+x2-3,则差商[3,32,33,34]=__1_____.
6.设l0(x),l1(x),l2(x),l3(x)是以x0,x1,x2,x3为互异节点
的三次插值基函数,则
3
l j
(x)(xj
2)3=___(_x_-_2_)_3____.
72
2 1
,所以‖A‖1=5,‖A-1‖1=5/7.
2.设矩阵A=
1 a
a 1
a 0
,当a取______值时,A可以唯一分解
a 0 1
为GGT,其中G为下三角矩阵.
1aa
解
令
1 a
a1a20,a 1
a
1 0
012a20, 得:1 a 1
1
2
2
3.向量x=(x1,x2,x3)T,试问|x1|+|2x2|+|x3|是不是一种向 量范数__是____,而|x1|+|2x2+x3|是不是一种向量范数_不_是___.
4.求 3 a 的Newton迭代格式为_x _k _1 __xk _ _x _3 k 3 _x k _2a _或 _x _k_ 1_ _3 2 _x _k_ _a 3 _x_k 2 _.
于是 H3(x)=-(x-1)2(x-2)-3x(x-2)+2.5x(x-1)2 –0.5x(x-1)(x-2) =x3-2.5x2 +2.5x+2
由于,R(0)=R(1)=R(2)=R(1)=0, 故可设 R(x)=C(x)x(x-1)2(x-2)
构造函数(t)=(t)-H3(t)-C(x)t(t-1)2(t-2) 于是,存在x,使(4)(x)=0,即(4)(x)-4!C(x)=0
x2(k1)
x(k) 2
(1
5
x(k1) 1
x(k) 2
1 5
x(k) 3
2) 5
x3(k1)
x(k) 3
(1
3
x(k1) 1
1 6
x(k1) 2
x(k) 3
1) 2
(2)因为A是严格对角占优矩阵,但不是正定矩阵,故
Jacobi法收敛,SOR法当0<1时收敛.
1 0 .75
四、(13分)已知(0)=2,(1)=3,(2)=5,(1)=0.5, (1)试建立一个三次插值多项式H3(x),使满足插值条件:
H3(0)=2,H3(1)=3,H3(2)=5,H3(1)=0.5; (2)设y=(x)在[0,2]上四次连续可微,试确定插值余项
(3)由(1)可见B=3/4,且取x(0)=(0,0,0)T,经计算可
得x(1)=(1/4,-2/5,1/2)T,于是x(1)-x(0)=1/2,所以有
B k x * x (1)0 x (1 ) x (0 )
0 .715 0 0 .5 0 .113
1 B
R(x)=(x)-H3(x).
解 (1)由y0=2,y1=3,y2=5,y1=0.5,得 H3(x)=20(x)+31(x)+52(x)+0.51(x) 令0(x)=c(x-1)2(x-2),可得0(x)=-0.5(x-1)2(x-2), 令1(x)=x(x-2)(ax+b),可得1(x)=-x(x-2), 令2(x)=cx(x-1)2,可得2(x)=0.5x(x-1)2; 令1(x)=cx(x-1)(x-2),可得1(x)=-x(x-1)(x-2),
取初值x0=1.5, 计算得x1=1.41333, x2=1.40983,由于 |x2-x1|=0.0035<10-2 , 故可取根的近似值x2=1.40983.
(3)因为0<</2,所以()co /2 s1sin 0
故,此迭代法线性收敛(收敛阶为1).
三、(14分)设线性方程组 4x1 x2 2x3 1 x1 5x2 x3 2 2x1 x2 6x3 3
j0
7.设S(x)=
x3x2 2x3b2xcx1
0x1 是以0,1,2为节
1x2
点的三次样条函数,则b=__-__2____c=___3______. 解 由2=b+c+1,5=6+2b+c,8=12+2b,可得
二、(13分)设函数(x)=x2-sinx-1 (1)试证方程(x)=0有唯一正根; (2)构造一种收敛的迭代格式xk=(xk),k=0,1,2,…计
( 2
k
1
)
1 4
x (k) 2
1 2
x (k) 3
1 4
1 5
x (k) 1
1 5
x (k) 3
2 5
x
( 3
k
1
)
1 3
x (k) 1
1 6
x (k) 2
1 2
x1(k1)
x(k) 1
(x1(k)
1 4
x(k) 2
1 2
x(k) 3
1) 4
考试题解析
一、填空题(每空3分,共30分)
1.设矩阵A= 1
2
2 3
Байду номын сангаас
,则(A)=____7 ___,Cond(A)1=___2_5/_7__.
解 由于
AE1
2
2 2 4 7 0
3
得特征值: 1 2 3 i,2 2 3 i
又A-1= 1 3
算精度为=10-2的近似根; (3)此迭代法的收敛阶是多少?说明之. 解 (1)因为0<x1时,(x)<0,x2时,(x)>0,所以(x)
仅在(1,2)内有零点,而当1<x<2时,(x)>0,故(x)单调. 因此方程(x)=0有唯一正根,且在区间(1,2)内.
(2)构造迭代格式: x k 1 1 sx ikn k 0 ,1 ,2 ,... 由于|(x)|=| co x/s 21sixn |<1,故此迭代法收敛.