2009南华大学硕研课程数值分析试题答案及评分标准
- 格式:doc
- 大小:220.50 KB
- 文档页数:5
1 南华大学 2009 级硕士研究生课程考试试题答案
及评分标准
考试科目: 数值分析 所属学院 考试时间
考生姓名: 考生学号 任课教师 考试成绩
一. (15分)已知函数1|)(|xf且)(xf的离散值如下表:
i 0 1 2
ix 0.32 0.34 0.36
)(ixf 0.314567 0.333487 0.352274
试用拉格朗日二次插值公式求)3367.0(f(保留小数点后6位)并估计截断误差。
(提示:拉格朗日n次插值公式:nkkknyxlxL0)()(,nkjjjkjkxxxxxl0)(
误差公式:njjnnnxxnfxLxfxR0)1()()!1()()()()()
评分标准:列出计算式5分,计算结果5分;误差估计5分;计算有误酎情扣1-2分;
用拉格朗日插值公式计算:
352274.0)34.036.0)(32.036.0()34.03367.0)(32.03367.0(333487.0)36.034.0)(32.034.0()36.03367.0)(32.03367.0(314567.0)36.032.0)(34.032.0()36.03367.0)(34.03367.0()3367.0(f(5分)
352274.0)02.0)(04.0()0033.0)(0167.0(333487.0)02.0)(02.0()0233.0)(0167.0(314567.0)04.0)(02.0()0233.0()0033.0(
330374.0330374364.00242672751.0324407818.00302338208.0(10分)
用牛顿插值公式计算:
00101013367.0)3367.0()3367.0(xxxyyyPf
330365.032.03367.032.034.0314567.0333487.0314567.0(5分) 2 100201011212123367.03367.0/)3367.0()3367.0()3367.0(xxxxxxyyxxyyPPf
330374.0)0033.0(0167.004.0946.093935.0330365.0(10分)
7210210102.20233.00033.00167.061||!3|)(|||xxxxxxfxR(15分)
二.(15分)求4x在[-1,1]上次数不超过3的最佳一致逼近多项式。
(提示:第一类切比雪夫多项式188)(244xxxT)。
评分标准:给出关系式5分;移项代入比雪夫多项式5分;计算5分;计算有误酎情扣1-2分
解:设最佳一致逼近多项式为)(*xP,则有
414*421)(TxPx,(10分)
811888181)(224444*xxxxTxxP(15分)
三.(15分)已知实验数据如下:
Xi 19 25 31 38 44
Yi 19.0 32.3 49.0 73.3 97.8
求形如cbxaxy2的经验公式,并估计均方误差。
(拟合公式提示:
nijjxaxS0)()(,nkdakjnjjk,,1,0,,0,21022miiiyxS
),,1,0,(,,0njkxxximijikijk,nkxxfxdmiikiik,,1,0,0)
评分标准:给出法方程5分;求出系数5分;给出均方误差5分;计算有误酎情扣1-2分;
令
144441383813131125251191922222A,cbax,8.973.730.493.3219Y
解 YAx,YAAxATT, 3 5 157 5327 157 5327 192331 5327 192331 7277699AAKT,0.00270.09783.6932005+1.0eYAfT(5分)
0.68820.01930.0497\fKx,0.6882 c0.0193, b0.0497, a
0.6882 0.0193x 0.0497x 2y(5分)
,0.0209 0.0636,- 0.0801, 0.0460,- 0.0086, YAUd,
ddsT,0.1145s(5分)
四.(15分)用四点高斯-勒让德公式求定积分1021xdxI(要求小数点后4位)。
(提示:四点公式系数:6521452.0,3478548.0,3399810.0,8611363.0kkAx)
评分标准:给出变换式5分;正确带入系数5分;正确计算5分;计算有误差酎情扣1-2分;
1,0ba令21222tbatabx,得5221211210ttdtxdxI(5分)
用四点高斯-勒让德公式求得
26521452.05)3399810.0(23399810.0153399810.023399810.0123478548.05)8611363.0(28611363.0158611363.028611363.012222I2}6521452.04.435615.795513478548.04.019317.46381{(5分)
26521452.00.22544730.17254623478548.00.24880060.1339795
20.65214520.39799360.34785480.382780
0.7854020.392701520.25954960.1331519(5分)
(精确值为33974480.785398164) 4 五.(15分)对方程8951251032241321xxx构造一种高斯-塞德尔迭代公式,
(1) 说明其收敛的理由;
(2) 以Tx000)0(为初值,迭代到2)()1(10kkxx。
(提示:高斯-塞德尔迭代:
bAx,ULDA,fBxxbLDfULDBkk)()1(11,,)
评分标准:给出变换后的对角占优式5分;判断收敛5分;迭代计算5分,每步1分;计算有误酎情扣1-2分;
910325241812521,238125910325241rrrr,1032241125A,958b,
1045D,321L,212U,
D=diag(diag(A)),L=-tril(A,-1),U=-triu(A,1),B=inv(D-L)*U, f=inv(D-L)*b.
125.005.0055.01.002.04.00)(1ULDB,
075.1650.1600.1)(1bLDf(5分)
λ1=0,λ2=-0.1125 + 0.1654i, λ3=-0.1125 - 0.1654i,
谱半径10.6512B,所以迭代
fBxxkk)()1(收敛, (5分)。取Tx0,0,0)0(
fBxxkk)()1(,(k=0,1,2,3,4,5)得下表
k Xk ||Xk-Xk-1||∞ ||Ax-b||2
0 (0 0 0) 13.0384
1 (1.6 1.65 1.075) 2.5373 4.8747
2 (0.7250 0.8937 1.0231) 1.1577 1.5678
3 (1.0379 0.9979 0.9918) 0.3312 0.1877
4 (1.0025 1.0047 1.0009) 0.0372 0.0292
5 (0.9979 0.9990 1.0001) 0.0073 0.0123
6 (1.0004 1.0000 0.9999) 0.0026 0.0019
7 (1.000000282,1.000032093, .000009572) 3.7424e-004 1.6531e-004
5 实际取k=5即可. (5分)注精确解为 Tx111
六.(25分)对于初值问题
2)1(2138yxyy
(1) 证明用改进的欧拉格式(梯形公式)是绝对稳定的;
(2) 取步长2.0h,用(1)的格式计算其数值解(小数点后保留5位)。
(提示:梯形公式:),(),(2111nnnnnnyxfyxfhyy)
评分标准:给出递推式10分,判断绝对稳定性5分;计算10分,每步2分;计算有误差酎情扣2-4分;
1138382nnnnyyhyy,解得2/3182/312/311hhyhhynn(10分)
对任意的0h,12/312/31rhh,故有nnnr1,误差是缩小的,所以格式绝对稳定。(5分)取2.0h,则,13161371nnyy ,由2)0(0yy计算得:
30769.22.11yy,(2分)
47337.24.12yy,(2分)
56258.26.13yy,(2分)
61062.28.14yy,(2分)
63649.20.25yy,(2分)