2011数值分析试题及答案

合集下载

2011级数值分析试卷

2011级数值分析试卷

菏泽学院数学系2011级 2013-2014学年第一学期数学与应用数学专业《数值分析和计算方法》期末试卷(A )(110分钟)题号 一 二 三 四 五 总分得分 阅卷人一.选择题(将正确选项前的代号写在题号前的括号内,每小题3分,共15分)( )1.若用最小刻度为0.5mm 的刻度尺测量物体,其误差限为( )A.0.25mmB.1.0mmC.0.5mmD.0mm ( )2.下列具有最高代数精度的求积公式是( )A.龙贝格求积公式B.复合辛普森求积公式C.牛顿-科特斯求积公式D.高斯求积公式( )3.已知2,1,0,,1)(==-=i i x x x f i i i 。

则函数)(x f 的插值多项式为( )A. 145412-+x x B.1-xC.-145412-+x x D.2+-x( )4.下列给出的是用不动点迭代法求032=-x 的根3*=x 的迭代函数,则相应的迭代方法局部收敛的是A.x x 3=)(ϕ B.3)(2-+=x x x ϕC.2321)(2-+=x x x ϕD.)3(21)(xx x +=ϕ( )5.线性方程组AX=b 能用高斯消元法求解的充要条件是( )A.A 为对称矩阵B.A.为实矩阵C.A 的各阶顺序主子式不为零D.0≠A得分 阅卷人二.填空题(请将正确答案填写在每小题的横线上,每空4分,共20分)1.计算积分⎰b adx x f )(的梯形公式为 。

2.设向量T n x )2,1,0( =,则=∞x 。

3.用牛顿法求方程0)(=x f 的根的公式为 。

4.已知n=3时的牛顿-科特斯系数83,83,81)3(2)3(1)3(0===C C C ,则=)3(3C 。

5.已知点,5,4,3,2,1,1=-=i i x i 则二阶差分=∆32x 。

三.判断题(对的在题前括号内划√,错的划×,每题2分,共10分)( )1.高斯求积公式的系数都是正的,故计算总是稳定的。

数值分析测试题答案

数值分析测试题答案
(3) 解:
构造另一种迭代公式: 1 xk + 1 = xk - (xk2 - 3),k = 0,1,2 4 1 1 ∴ψ ( x) = x - (x 2 - 3),ψ' ( x) = 1 - x, 4 2 ψ' ( x * ) =ψ'
( 3) = 1-
3 = 0.134 < 1,收敛。 2
5π π π ( ) sin f x = sin x , x ∈ [ , ] 5. 已知函数 18 ,给出误差估计。 6 3 ,在函数图像上取型值点做二次插值多项式,计算
∴ N (x) = -11 +14(x + 2) -12(x + 2)(x +1) +12x(x +1)(x + 2)
7. 插值与下面表列 i xi yi 的自然样条定义如下:
3 s x a x 1 b x 1 1,1 x 2 0 S x 3 2 3 s x c x 2 x 2 d x 2 1, 2 x 3 1 4
i 0 i 0
4
4
于是得到如下方程: 8a0 22a1 47 22a0 74a1 145.5 解得 a0 2.77, a1 1.13 于是所求的最小二乘法拟合曲线为
* y s1 x 2.77 1.13x
均方误差为
2
=0.9424
12/1/2011 12:23:43 PM
数值分析测试题参考答案
1. 写出有效数字的概念,并计算 2.718 近似 e 有几位有效数字。(e = 2.71828182…) 解:(1)设 x*是 x 的一个近似数,表示为

福州大学2010-2011年数值分析考题及答案1

福州大学2010-2011年数值分析考题及答案1
得分 评卷人
1、若向量 x (4, 2,3) ,则
T
x 2 =___ 29 _________
=____ 6 ____,A 的
2、
1 1 A , 则 A 的谱半径 -5 1
=____6____
3、 确定求积公式 尽量高,则 A0=_

1
1
f ( x)dx A0 f (1) A1 f (0) A2 f '(1) 中的待定参数,使其代数精度
0 2 0 5、设 B 2 1 2 ,试用平面旋转矩阵对矩阵 A 进行 QR 分解,其中 Q 为正交 0 2 1
矩阵,R 为上三角阵(8 分)
4
记A1 A, 先将A的第一列变得与e1平行 cos = 0 2 0,sin = 1 04 04 0 1 0 0 1 0 0 P A 2 P A1 1 12 12 0 0 0 1
3、
h 用二步法 yn1 yn [ f ( xn , yn ) f ( xn1 , yn1 )] 求解一阶常微分方程初值问题 2
y f ( x, y ) 问:如何选择参数 , 的值,才使该方法的阶数尽可能地高?写出 y ( x0 ) y0
此时的局部截断误差主项,并说明该方法是几阶的。 证明:局部截断误差为:
( x x )l ( x) 等于
i 0 i i
4
( a ) 1 (c) 2 (d) 4
(a)
0
(b)
3、设 f ( x) 3x5 4 x 4 x 2 1 和节点 xk k / 2, k 0,1 则差商 f [ x0 , x1 x5 ] (a) 4 (b) 2 (c) 3 (d) 1 ( ( c ) c )

数值分析试卷及答案

数值分析试卷及答案

数值分析试卷及答案数值分析试卷一、选择题(共10题,每题2分,共计20分)1. 数值分析的研究内容主要包括以下哪几个方面?A. 数值计算方法B. 数值误差C. 数值软件D. 数学分析答:A、B、C2. 下列哪种方法不属于数值积分的基本方法?A. 插值法B. 微积分基本公式C. 数值微积分D. 数值积分公式答:A3. 数值积分的目的是求解什么?A. 函数的导数B. 函数的原函数C. 函数的极值D. 函数的积分答:D4. 数值微分的目的是求解什么?A. 函数的导数B. 函数的原函数C. 函数的极值D. 函数的积分答:A5. 数值微分的基本方法有哪几种?A. 前向差分B. 后向差分C. 中心差分D. 插值法答:A、B、C6. 用数值方法求解方程的基本方法有哪几种?A. 迭代法B. 曲线拟合法C. 插值法D. 数值积分法答:A、B、C7. 用迭代法求方程的根时,当迭代结果满足何条件时可停止迭代?A. 当迭代结果开始发散B. 当迭代结果接近真实解C. 当迭代次数超过一定阈值D. 当迭代结果在一定范围内波动答:B8. 下列哪种插值方法能够确保经过所有给定数据点?A. 拉格朗日插值B. 牛顿插值C. 三次样条插值D. 二次插值答:A、B、C9. 数值解线性方程组的基本方法有哪几种?A. 直接法B. 迭代法C. 插值法D. 拟合法答:A、B10. 下列哪种方程求解方法适用于非线性方程?A. 直接法B. 迭代法C. 插值法D. 曲线拟合法答:B二、填空题(共5题,每题4分,共计20分)1. 数值积分的基本公式是_________。

答:牛顿-科特斯公式2. 数值微分的基本公式是_________。

答:中心差分公式3. 数值积分的误差分为_________误差和_________误差。

答:截断、舍入4. 用插值法求解函数值时,通常采用_________插值。

答:拉格朗日5. 数值解线性方程组的常用迭代法有_________方法和_________方法。

2011年秋季工学硕士研究生学位课程(数值分析)真题试卷B

2011年秋季工学硕士研究生学位课程(数值分析)真题试卷B

2011年秋季工学硕士研究生学位课程(数值分析)真题试卷B(总分:28.00,做题时间:90分钟)一、填空题(总题数:6,分数:12.00)1.填空题请完成下列各题,在各题的空处填入恰当的答案。

(分数:2.00)__________________________________________________________________________________________ 解析:2.设|x|>>1______(分数:2.00)__________________________________________________________________________________________正确答案:()解析:3.求积分∫ a b f(x)dx的两点Gauss公式为______(分数:2.00)__________________________________________________________________________________________正确答案:()解析:4.设∞ =______,‖A‖ 2 =______.(分数:2.00)__________________________________________________________________________________________正确答案:()解析:5.给定f(x)=x 4,以0为三重节点,2为二重节点的f(x)的Hermite插值多项式为______.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:x 4)解析:6.己知差分格式r≤______时,该差分格式在L ∞范数下是稳定的.(分数:2.00)__________________________________________________________________________________________正确答案:()解析:二、计算题(总题数:2,分数:4.00)7.给定方程lnx-x 2+4=0,分析该方程存在几个根,并用迭代法求此方程的最大根,精确至3位有效数字.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:令f(x)=lnx-x 2 +4,则f"(x)= -2x,当x= 时,f"(x)=0. 注意到f(0.01)=-0.6053<0,f(1)=3>0,f(3)=-3.9014<0,而当时,f"(x)>0,当时,f"(x)<0,所以方程f(x)=0有两个实根,分别在(0.01,1)和(1,3)内.方程的最大根必在(1,3)内,用Newton迭代格式取x 0 =2,计算得x 1 =2.1980,x 2 =2.1)解析:8.用列主元Gauss(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:求得x 1 =3,x 2 =1,x 3 =5.)解析:三、综合题(总题数:6,分数:12.00)9.设α,β表示求解方程组.Ax=b的Jacobi迭代法与Gauss-Seidel迭代法收敛的充分必要条件.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:Jacobi迭代格式的迭代矩阵特征方程为展开得500λ3—15αβλ=0或者λ(500λ2—15αβ)=0,解得λ=0或λ2 = 则Jacobi格式收敛的充要条件为|αβ|<Gauss-Seidel格式迭代矩阵的特征方程为展开得500λ3—15αβλ2 =0或者λ2(500λ-15αβ)=0,解得λ=0或λ则Gauss-Seidel格式收敛的充)解析:10.设x 0,x 1,x 2为互异节点,a,b,m为已知实数.试确定x 0,x 1,x 2的关系,使满足如下三个条件p(x 0 )=a, p"(x 1 )=m,p(x 2 )=b的二次多项式p(x)存在且唯一,并求出这个插值多项式p(x).(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:由条件p(x 0 )=a,p(x 2 )=b确定一次多项式p 1 (x),有所以p(x)-P 1(x)=A(x—x 0 )(x—x 2 ),p"(x)=p" 1 (x)+A(x—x 0 +x—x 2 ),p"(x 1+A(2x 1 -x 0 -x 2) 解析:11.求y=|x|在[-1,1]上形如c 0 +c 1 x 2的最佳平方逼近多项式.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:取φ0 (x)=1,φ1 (x)=x 2,则(φ0,φ0)=∫ -11 =2,(φ0,φ1)=∫ -11 x 2)1 x 2,(φ1,φ1)=∫ -1解析:12.已知函数f(x)∈C 3 [0,3],试确定参数A,B,C,使下面的求积公式数精度尽可能高,并给出此时求积公式的截断误差表达式.(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:当f(x)=1时左=∫ 03 1dx=3,右=A+B+C,当f(x)=x时左=∫ 03 xdx= ,右=B+2C 当f(x)=x 2时左=∫ 03 x 2 dx=9,右=B+4C.要使公式具有尽可能高的代数精度,则而当f(x)=x 3时,左=∫ 03 x 3)解析:13.给定常微分方程初值问题取正整数n,并记h=a/n,x i =a+ih,0≤i≤n.证明:用梯形公式求解该初值问题所得的数值解为且当h→0时,y n收敛于y(a).(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:梯形公式应用于方程有y i+1=y i+ (-y i—y i+1),即有所以i=1,2,….当h→0时,n→∞我们有而由方程知解析解y=e -x则y(a)=e -a,所以)解析:14.Ω={0<x<3,0<y<3).试用五点差分格式求u(1,1),u(1,2),u(2,1),u(2,2)的近似值.(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:五点差分格式为根据要求,可取h= ,将(1,1),(2,1),(1,2),(2,2)处的差分格式列成方程组有或者解得u 11=15.8750,u 21=22.6250,u 12=15.8750,u 22 =22.6250.)解析:。

中国石油大学《数值分析》2011年考试试题A卷及答案

中国石油大学《数值分析》2011年考试试题A卷及答案

f (4)(x)
1 2880
1 n
4
6
1 2
104
,
仅要 n 4 1 101 2.54 ,取 n 3 即对将[1,2] 作 6 等分,则有 240
(8 分)
2
1 ln xdx
1 [0 4(ln 7 ln 3 ln 11) 2(ln 4 ln 5) ln 2] 0.38628716327880 .
0.000040074
( 4 分)
七、(10 分)(1)牛顿迭代格式
x(k 1)
x(k)
f f
(x(k ) ) '(x(k) )
x(k)
x(k) 1 (2
(x(k) )2 )(x(k) )1
1
(1 (2
)(
x( )(
)k ) 2 x(k ) )1
(2)
x(k 1)
lim
k
x(k)
1 1
fgdx
,取( x) ax bx3 , f ( x) sin x ,则法方程为
(0 ,0 )
(1
,
0
)
(0 ,1) (1 , 1 )
a b
( (
f f
,0 ,1
) )
( 4 分)
其中 0,0
1
x xdx
1
2, 3
0 ,1
(1 )(x(k) )2
lim
k
1
(2
)(x(k ) )1
c0
2
c 1
(5 分) (5 分)
1
x(k) 2
x(k) 3
1
x(k) 1
x(k) 3
/2
x3( k
1)

2011年哈工大(数值分析)试题

2011年哈工大(数值分析)试题

y0 , y1 , (h为步长) 。
(1)确定方法中的局部误差主项,并指出方法的阶数; (2)讨论该方法的收敛性和绝对稳定性。 ( 在 线 性 多 步 法 的 局 部 截 断
Cr
p 1 p p 1 (i ) r 1 bi , r 2,3, ) i a r i r ! i 0 i 1
1 2
a 2 1 x1 1 3、已知方程组 2 a 2 x2 2 , 1 2 a x3 1
(1)写出求解此方程组的 Jacobi 迭代格式; (2)用已知结论说明,当 a 4 时,该迭代格式收敛;
0 2
T
(1)求 A , x0 , x1 使求积公式具有尽可能高的代数精度,并指出此求积公式的 代数精度是多少?(2)并用此公式计算积分 x 4 dx 。 (计算结果保留四位小数)
0 2
ρ 。 6、试用共轭梯度法(cg 法)求解线性方程组。 (初始值取 x (0) (0, 0, 0)T )
a x 2 x3 , 0 x 1 ,具有连续二 2、 (1)求 a 及不超过二次多项式 P( x) 使 S ( x) P( x) , 1 x 2
阶导数且满足 P (2) 0 ; (2) 当 f ( x) 用满足条件 f (1) P(1), f (2) P(2), f ' (1) P ' (1) 的插值多项式近似时, 求 f ( x)dx 。
2011 年哈工大《数值分析》考试
1、 设 f ( x) x3 5 。 (1)应用 Newton 迭代法于方程 f ( x) 0 ,导出 3 5 的迭
2
代公式; 并讨论迭代公式的收敛速度; (2) 尝试把导出的迭代公式加以改进, 提高迭代公式的收敛速度, 并用改进后的迭代公式计算 3 5(取初值 x0 1.0 , 计算三步,结果保留四位小数) 。

西安石油大学研究生数值分析10 11年试题

西安石油大学研究生数值分析10 11年试题
2010/2011学年第I学期数值分析考试题(卷)
一、填空题(每题2分,共20分) 1.近似数 x =0.231关于真值x=0.229有

位有效数字。 。
n
2.求方程 f ( x) 0 的根时,对应的牛顿切线法迭代公式为 3.设 l i ( x) (i=0,1,2,…,n)是n次拉格朗日插值基函数,则
4 0 x1 5 2 3 1 1 x 2 9 2 2 0 x 3 3
四、(12分)写出解线性方程组
4 x1 2 x3 4 x1 4 x 2 2 x3 1 的高斯—赛德尔迭代法的迭代格式,并判断其收敛性。 3 x 5 x x 2 2 3 1
l ( x) =
i 1 i

4.求解微分方程初值问题
y ' f ( x, y ) 时,设x节点步长为h,则欧拉预估— y ( x0 ) y 0
迭代法和
校正方法的局部截断误差为 。 5.若线性方程组AX=b的系数矩阵A为严格对角占优矩阵,则 迭代法收敛。 6.差商与向前差分满足关系: 差商与向后差分满足关系: 7.用数值方法求积分 。 。
五、(12分)已知一组观察数据为 i 0 1 2 2 3 3 4
xi
1
yi
0
-5
-6
3
试用此组数据构造3次牛顿插值多项式 N 3 ( x) ,并计算 N 3 (1.5) 的值。 六、(12分)试确定经验公式 y ae 中的参数a和b(a为正数),使该函数曲线与下列数
bx
据按最小二乘原则相拟合(至少保留ห้องสมุดไป่ตู้位小数)。 1 2
xi
3 20

2011年秋研究生数值分析试题A卷答案

2011年秋研究生数值分析试题A卷答案

2011年秋研究生数值分析期末考试试题答案一、单选题(4*5=20分)1、B;2、D ;3、D ;4、B ;5、C 。

二、填空题(4*5=20)1、2;2、()()1k k k k f x x x f x +=-',平方收敛;3、8,8;4、9; 5、a <。

三、(10分)解:构造3次Lagrange 插值多项式3001001201()()(,)()(,,)()()L x f x f x x x x f x x x x x x x =+-+--0123012(,,,)()()()f x x x x x x x x x x +--- 3’利用待定系数法,令430123()()()()()()H x L x A x x x x x x x x =+----, 5’同时, '''14131101213()()()()()()f x H x L x A x x x x x x ==+--- 7’解出A 即可。

8’ 考虑余项4()()()E x f x H x =-,如果5()[,],,0,1,2,3i f x C a b a x b i ∈≤≤=,那么,当a x b ≤≤时()()5240123()()()()()()()5!f E x f x H x x x x x x x x x ξ=-=----. 0 10’ 四、(10分)解:设所求多项式为23202)(x C x C C x P ++=,10=ϕ,x =1ϕ,22x =ϕ,11),(10++==⎰+k j dx e k j k j ϕϕ,1),(100-==⎰e dx e f x ϕ, 1),(101==⎰dx xe f xϕ,2),(1022-==⎰e dx e x f x ϕ 5’ 所以有⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡21151413141312131211210e e C C C ,求解得到 8’ ⎪⎩⎪⎨⎧===83917.085114.001299.1321C C C ,所求最佳平方逼近多项式为:2283917.085114.001299.1)(x x x P ++=。

2011数值分析试题及答案

2011数值分析试题及答案
122446
由于f(x)二si nx的4阶导数在[0,二]上的最大值为:M4=1,所以
5
误差为:|I-S2|::——44=0.006641
2880x24
6.求解初值问题」y=sin(x+2y),0兰x兰2的改进Euler方法是否收敛?为什
.y(0) = 1
么?
解:由于|sin(x 2y)-sin(x 2y)|二| 2cos(x 2 )(y-y) 2 | y-y |
5.设f(x) = 4x33x-5,求差商f[0,1], f[1,2,3,4]和f[1,2,3,4,5]。
f(D…f(0)
解:f[0,1]==2-(-5) = 7
1-0
f [1,2,3,4^4,f[1,2,3,4,5]=0
3.解线性方程组丿X1-2忑=2的Jacobi迭代法是否收敛,为什么?
+9x2=3
即,函数f(x, y)二sin(x•2y)连续,且关于变量y满足Lipschitz条件,所以,改 进Euler方法收敛。
所以,a=0, b=5/6,拟合曲线为:y=5/6x2
3.求满足条件f(0)=1,f(1)=2,f(2) =0,f(1)=0的三次插值多项式Ha(x)
的表达式。
解:设H3(x)二(^2)(ax2bx c),则有:
1213
所以,H3(x) (x-2)(x2x 1) (x-3x-2)。
22
11
4.确定求积公式Jf(x)dx痒三f(-1)+Af(0)+A2f(1)中的待定系数,使其代数精 度尽可能高,并问此公式是不是插值型求积公式.
解:令公式对f(x) = 1,x都精确成立,得:A,・A2= 3/2, A2= 1/2,
o
• • •

2011年秋季工学硕士研究生学位课程(数值分析)真题试卷A

2011年秋季工学硕士研究生学位课程(数值分析)真题试卷A

2011年秋季工学硕士研究生学位课程(数值分析)真题试卷A(总分:28.00,做题时间:90分钟)一、填空题(总题数:6,分数:12.00)1.填空题请完成下列各题,在各题的空处填入恰当的答案。

(分数:2.00)__________________________________________________________________________________________ 解析:2.已知x 1 =0.724,x 2 =1.25均为有效数,则|e r (x 1 x 2 )|≤______|e(x 1/x 2 )|≤_______.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:0.469×10 -2,0.272×10 -2)解析:3.设∞ =______,cond(A) 2 =______.(分数:2.00)__________________________________________________________________________________________正确答案:()解析:4.超定方程组x 1 =______.x 2 =_______.(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:-0.8333或-0.6667)解析:5.用Simpson______.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:1.4757)解析:6.0,1,2为节点的三次样条函数,则a=_____,b=_____.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:3,-3)解析:二、计算题(总题数:3,分数:6.00)7.给定方程e x,分析此方程有几个实根,并用迭代法求此方程的正根,精确至3位有效数字.(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:设当x=ln0.5时,f"(x)=0;当x∈(-∞,ln0.5)时,f"(x)<0;当x∈(ln0.5,+∞)时,f"(x)>0.再注意到f(-4)>0,f(-3)<0,f(1)<0,f(2)>0,则该方程存在两个实根,分别在[-4,-3]和[1,2]内.构造迭代格式x *∈[1,2],取x 0 =1.5,计算得x 1 =1.0651,x 2 =0.9116,x 3 =0.8953,x )解析:8.用列主元Guass(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:求得x 1 =-2,x 2 =1,x 3 =-1.)解析:9.给定求解线性方程组Ax=b的迭代格式Bx (k+1) +ωCx k =b,其中ω的值使上述迭代格式收敛.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:方法1:由Bx (k+1) +ωCx (k) =b得x (k+1) =-ωB -1 Cx (k) +B -1 b 上述格式收敛的充要条件为ρ(-ωB -1 C)<1.迭代矩阵-ωB -1 C的特征方程为|λI+ωB -1 C|=0,可变形为|B -1||λB+ωC|=0,即展开得16λ2—8λω)解析:三、综合题(总题数:5,分数:10.00)10.作一个3次多项式H(x),使得H(a)=b 3,H(b)=a 3,H"(a)=6b,H"(b)=6a.(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:方法1:根据H"(a)=6b,H"(6)=6a可知(x-a)=6b—6(x-a),两边积分得H(x)=6b(x—a)-3(x-a) 2 +c,H(x)=3b(x-a) 2 -(x—a) 3 +c(x—a)+d.由H(a)=b 3得d=b 3,再由H(b)=a 3有c=-3b 2,所以H(x)=-(x-a) 3 +3b(x-a) )解析:11.求函数y(x)=x 4在区间[0,1]上的一次最佳一致逼近多项式p(x).(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:设p(x)=a+bx.由f"(x)=4x 3,f"(x)=12x 2知,当x∈(0,1)时,f"(x)恒大于零.则f(x)-p(x)在[0,1]上有三个交错偏差点:0,x 1,1,且满足即求解得所以)解析:12.已知函数f(x)∈C 4[a,b],I(f)=∫ a b f(x)dx 1)写出以a,b为二重节点所建立的f(x)的3次Hermite 插值多琐式H(x)及插值余项; 2)根据f(x)≈H(x)建立一个求解I(f)的数值求积公式I H (x),并分析该公式的截断误差和代数精度.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:1)由条件H(a)=f(a),H"(a)=f"(a),H(b)=f(b),H"(b)=f"(b),作差商表:所以2)根据题意,有I(f)≈∫ a b H(x)dx,下面求代数精度.由插值余项知,当f(x)=1,x,x 2,x 3时,插值余项为零,I H (f)精确求积;当f(x)=x 4时此时b 5系数为I H) 解析:13.给定常微分方程初值问题n,并记h=(b—a)/n,x i=a+ih,0≤i≤n.试确定参数A,B,C,使求解公式y i+1 =Ay i +(1-A)y i-1 +h[Bf(x i+1,y i+1 )+Cf(x i,y i )]的局部截断误差R i+1的阶数达到最高,指出所达剑的最高阶数并给出局部截断误差表达式.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:局部截断误差为R i+1=y(x i+1)-Ay(x i)-(1-A)y(x i-1)-h[By"(x i+1)+Cy"(x i)]=y(xi )+hy"(x 1 )+ y"(x i )+ y""(x i y (4) (xi )+O(h5)解析:14.给定如下抛物方程初边值问题:取步长用古典隐格式计算u(x,t)(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:求解该问题的古典隐格式为记则差分格式可写为(1+2r)u i k-r(u i+1k +ui-1k )=uik-1 +τ(3—3xi ),用方程组表示为k=1.2.因为所以,当k=1时,方程为或解得u 11 =0.7870,u 2<) 解析:。

[考研类试卷]2011年工程硕士研究生学位课程(数值分析)真题试卷B.doc

[考研类试卷]2011年工程硕士研究生学位课程(数值分析)真题试卷B.doc

[考研类试卷]2011年工程硕士研究生学位课程(数值分析)真题试卷B
1 设x=1.231,y=0.5122是由四舍五入法得到的近似值,试计算函数e xy的绝对误差限和相对误差限.
2 给定方程x3+2x-1=0,判别该方程有几个实根,并用迭代法求出方程所有实根,精确到4位有效数字.
3 用列主元Gauss 消去法求下面线性方程组的解:
4 给定线性方程组写出求解上述方程组的Gauss-Seidel 迭代格式,并分析收敛性.
5 已知f(x)=xe x,求一个3次多项式H(x),使之满足H(0)=f(0),H(1)=f(1),
H'(0)=f'(0),H"(1)=f"(1).
6 求a,b ,使得积分取最小值.
7 试用Simpson 公式计算积分的近似值,精确到4位有效数字.
8 给定常微分方程初值问题取正整数n,记h=(b—a)/n,
x i=a+ih,i=0,1,2,…,n;y i≈y(x i),1≤i≤n,y0=η.求常数A,B,使数值求解公
式y i+1=y i十h[A,(x i+1,y i+1)+f(x i,y i)+Bf(x i-1,y i-1)],1≤i≤n-1的阶数尽可能高,并
求出公式的阶数和局部截断误差表达式.
答案见麦多课文库。

华南农业大学2011-2012第2学期数值分析试题及答案

华南农业大学2011-2012第2学期数值分析试题及答案

华南农业大学期末考试试卷(A 卷)2011-2012学年第 2 学期 考试科目: 数值分析 考试类型:(闭卷)考试 考试时间: 120 分钟 学号 姓名 年级专业一、 填空题(本大题共5小题,每小题6分,共30分)1、用四舍五入得到的近似数1.55,有___位有效数字,其相对误差限是________。

2、用二分法求方程1x xe =在区间[0,1]内的根,进行一步后根所在区间为____,进行二步后根所在区间为____。

3、采用牛顿迭代法求正实数a 的开平方,迭代公式为________。

4、设有矩阵131211122A ⎡⎤⎢⎥=-⎢⎥⎢⎥-⎣⎦,则A ∞=____,F A =____。

5、非线性方程迭代求解的敛散性与初始值的选取____;(选填:有关或无关)线性方程组迭代求解的敛散性与初始值的选取____。

(选填:有关或无关)试用4[4,4]P 和4[4,4,4]P 计算'(4)f 和''(4)f 的近似值。

(本题共10分)三、给定方程30(0)x e x x --=>(1) 分析方程存在几个解,并找出解的范围;(2) 将方程改为3x x e =-,写出相应的迭代公式,并说明能不能用该公式迭代求原方程的解;(3) 如果不能,试将方程改写为能用迭代法求解的形式,并说明理由。

(本题共15分)四、分别讨论用雅可比迭代法和高斯-赛德尔迭代法解方程AX b =的收敛性,其中122111221A -⎛⎫ ⎪=-- ⎪ ⎪--⎝⎭。

(本题共16分)五、已知函数3y x =的函数表如下:(1) 求3y x =的3次拉格朗日插值多项式;(2) 求3y x =的3次牛顿插值多项式。

(本题共14分)六、采用龙贝格法计算140I x dx =⎰的值。

(本题共15分)华南农业大学期末考试试卷参考答案(A 卷)2010-2011学年第 2 学期 考试科目: 数值分析 一、填空题(本大题共5小题,每小题6分,共30分)1、用四舍五入得到的近似数1.55,有3位有效数字,其相对误差限是0.323%。

2011级硕士研究生《数值分析》试卷(A)

2011级硕士研究生《数值分析》试卷(A)

合肥工业大学2011级硕士研究生《数值分析》试卷(A)班级 姓名 学号 成绩一、判断题 (下列各题,你认为正确的,请在题后的括号内打“√ ”,错误的打“×”,每题2分,共10分) 1. 设函数f 具有5阶导数,则(5)[0,1,2,3,4,5]()f f ξ=,其中ξ介于0,1,2,3,4,5之间,[0,1,2,3,4,5]f 是()f x 关于节点0,1,2,3,4,5的5阶差商。

( )2. 若方阵A 是严格对角占优的,则可用Gauss 消去法直接求解方程组=Ax b ,无须选主元素。

( )3. 若()()0f a f b <,则方程()0f x =在区间(,)a b 内至少有一个根。

( )4. 若函数()f x 是多项式,则它的Lagrange 插值多项式()()p x f x ≡. ( )5. 解常微分方程初值问题的四阶Runge-Kutta 方法的局部截断误差是5()O h ,其中h 是步长。

( )二、填空题 (每空2分,共10分)1. 近似数*3.200x =关于准确值 3.200678x =有 位有效数字。

2. 设2435A =⎡⎤⎢⎥⎣⎦,则1Cond()A = . 3. 设函数(2.6)13.4673,(2.7)14.8797,(2.8)16.4446f f f ===, 用三点数值微分公式计算(2.7)f '= 14.8865 .4. 设函数sin 2()x f x =, 2()p x 是()f x 的以1,2,3为节点的二次Lagrange 插值多项式,则余项2()()f x p x -= .5. 二元函数(,)f x y 在区域D 上关于y 满足Lipschitz 条件是:.三 (本题满分12分) 对下列方程组1231231235212,4220,23103x x x x x x x x x ++=-⎧⎪-++=⎨⎪-+=⎩ 建立Jacobi 迭代格式(4分)和Gauss –Seidel 迭代格式(4分),写出Jacobi 迭代格式的迭代矩阵,并用迭代矩阵的范数判断所建立的Jacobi 迭代格式是否收敛(4分)。

数值分析、计算方法试题库及答案

数值分析、计算方法试题库及答案
4、设y=sinx,当取x0=1.74,x1=1.76,x2=1.78建立拉格朗日插值公式计算x=1.75的函数值时,函数值y0,y1,y2应取几位小数?
5、已知单调连续函数y=f(x)的如下数据:
xi
-0.11
0.00
1.50
1.80
f(xi)
-1.23
-0.10
1.17
1.58
若用插值法计算,x约为多少时f(x)=1。(计算时小数点后保留5位)。
《 计算方法》课程试卷A参考答案及评分标准 第2页 共3页
七(8分)、证: 的元素为 ,
因此 为对称矩阵。
记 ,则
对任意n-1维非零向量 ,作 ,记 ,则 ,
而 ,从而 为正定矩阵。
《 计算方法》课程试卷A参考答案及评分标准 第3页 共3页
课程编号:12000044北京理工大学2010-2011学年第一学期
2、用牛顿——埃尔米特插值法求满足下列表中插值条件的四次插值多项式P4(x),并写出其截断误差的表达式(设f(x)在插值区间上具有直到五阶连续导数)。
xi
0
1
2
f(xi)
1
-1
3
f’(xi)
1
5
3、对下面的线性方程组变化为等价的线性方程组,使之应用雅克比迭代法和高斯——赛德尔迭代法均收敛,写出变化后的线性方程组及雅克比迭代法和高斯——赛德尔迭代法的迭代公式,并简单说明收敛的理由。
四、填空题(2 0×2′)
15.设x=0.231是精确值x*=0.229的近似值,则x有2位有效数字。
16.设 ,‖A‖∞=___5____,‖X‖∞=__3_____,
‖AX‖∞≤_15___。
17.非线性方程f(x)=0的迭代函数x=(x)在有解区间满足|’(x)| <1,则使用该迭代函数的迭代解法一定是局部收敛的。

数值分析(2011)试题A卷 参考答案

数值分析(2011)试题A卷  参考答案

装订线年 级 学 号 姓 名 专 业一、填空题(本题40分, 每空4分)1.设),,1,0()(n j x l j =为节点n x x x ,,,10 的n 次基函数,则=)(i j x l 1,0,1,,0i j i j n i j=⎧=⎨≠⎩ 。

2.已知函数1)(2++=x x x f ,则三阶差商]4,3,2,1[f = 0 。

3.当n=3时,牛顿-柯特斯系数83,81)3(2)3(1)3(0===C C C ,则=)3(3C 1/8 。

4.用迭代法解线性方程组Ax=b 时,迭代格式 ,2,1,0,)()1(=+=+k f Bx x k k 收敛的充分必要条件是 ()1B ρ< 。

5.设矩阵⎥⎦⎤⎢⎣⎡=1221A ,则A 的条件数2)(A Cond = 3 。

6.正方形的边长约为100cm ,则正方形的边长误差限不超过 0.005 cm才能使其面积误差不超过12cm 。

(结果保留小数)7.要使求积公式)()0(41)(111x f A f dx x f +≈⎰具有2次代数精确度,则 =1x23 , =1A 34。

8. 用杜利特尔(Doolittle )分解法分解LUA =,⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=1359 45- 279 126 0 945- 0 45 1827- 9 189A 其中,则=L 10002100121023113⎛⎫⎪ ⎪ ⎪-⎪ ⎪- ⎪⎝⎭=U 918927091890281540009-⎛⎫⎪-⎪ ⎪-⎪⎝⎭。

二、计算题(10分)已知由数据(0,0),(0.5,y ),(1,3)和(2,2)构造出的三次插值多项式)(3x P 的3x 的系数是6,试确定数据y 。

2011级数值分析 试题 A 卷 2011 ~ 2012学年,第 1 学期一 二 三 四 五 六 七 八 九 十 总分年 级2011级研究生 份 数 拟题人 王吉波 审核人装 订线年级 学 号 姓 名 专 业三、计算题(15分)试导出计算)0(1>a a的Newton迭代格式,使公式中(对n x )既无开方,又无除法运算,并讨论其收敛性。

2011年工程硕士研究生学位课程(数值分析)真题试卷B

2011年工程硕士研究生学位课程(数值分析)真题试卷B

2011年工程硕士研究生学位课程(数值分析)真题试卷B(总分:16.00,做题时间:90分钟)一、计算题(总题数:3,分数:6.00)1.设x=1.231,y=0.5122是由四舍五入法得到的近似值,试计算函数e xy的绝对误差限和相对误差限.(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:根据题意,可知|e(x)|≤ ×10 -3,|e(y)|≤ ×10 -4,|e(e xy)|≈|ye xy e(x)+xe xy e(y)|≤e xy (y|e(x)|+x|e(y)|)≤0.5967×10 -3,)解析:2.给定方程x 3 +2x-1=0,判别该方程有几个实根,并用迭代法求出方程所有实根,精确到4位有效数字.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:改写方程为x 3 =-2x+1,作函数y=x 3和y=-2x+1的图像(见下图),由图像知方程有一个实根x *∈ 构造Newton迭代格式:x k+1=x k k=0,1,2,…,取初值x 0=0.25,计算得x 1 =0.47143,x 2 =0.45357,x 3 =0.45340,x 4 =0)解析:3.用列主元Gauss(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:求得x 1 =1,x 2 =1,x 3 =-1,x 4 =-1.)解析:二、综合题(总题数:5,分数:10.00)4.Gauss-Seidel迭代格式,并分析收敛性.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:Gauss-Seidel迭代格式为Gauss—Seidel迭代矩阵的特征方程为展开得4λ3—4λ2—2λ+8λ2 +2λ2—2λ2 =0,即λ(2λ2 +2λ-1)=0,求得λ1 =0,因为所以Gauss-Seidel迭代发散.)解析:5.已知f(x)=xe x,求一个3次多项式H(x),使之满足H(0)=f(0),H(1)=f(1),H"(0)=f"(0),H"(1)=f"(1).(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:作2次插值多项式p(x),满足p(0)=f(0) p(1)=f(1),P"(0)=f"(0),则p(x)=f(0)+f[0,0]x+y[0,0,1]x 2.列表求差商:可得p(x)=x+(e-1)x 2.由插值条件易知H(x)=p(x)+Ax 2 (x-1),其中A为待定系数.由条件H"(1)=f"(1)得2(e-1)+4A=3e,求得A= 所以H(X)=x+(e-1)x 2x 2 ()解析:6.求a,b(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:取φ0 (x)=1,φ1 (x)=x 2,则(φ0,φ0)=∫ -11 1dx=2,(φ0,φ1)=∫-11 x 2,(φ1,φ0)=∫ -11 x 2)解析:7.试用Simpson4位有效数字.(分数:2.00)__________________________________________________________________________________________正确答案:(正确答案:由复化Simpson公式得S 1(f)= [f(1)+4f(1.5)+f(2)]= (e -1+4e -1.52+e -22)=0.13463,S 2(f)= [f(1)+4f(1.25)+2f(1.5)+4f(1.75)+f(2)]=0.13521,因为|S 2 (f)-S 1 (f)|=0.38675×10)解析:8.给定常微分方程初值问题取正整数n,记h=(b—a)/n,x i=a+ih,i=0,1,2,…,n;y i≈y(xi ),1≤i≤n,y 0 =η.求常数A,B,使数值求解公式y i+1 =y i十h[A,(x i+1 ,y i+1f(x i,y i )+Bf(x i-1 ,y i-1 )],1≤i≤n-1的阶数尽可能高,并求出公式的阶数和局部截断误差表达式.(分数:2.00)__________________________________________________________________________________________ 正确答案:(正确答案:求解公式的局部截断误差为 R i+1 =y(x i+1 )-y(x i )-Ahf(x i+1,y(x i+1 ))-hy(x i,y(x i ))-Bhf(x i-1,y(x i-1 ))=y(x i+1 )-y(x i )-Ahy"(x i+1hy"(x )解析:。

数值分析A卷(2011年秋)

数值分析A卷(2011年秋)

三、 (10 分)设 f ( x) 在 [ x0 , x3 ] 上有 5 阶连续导数,且 x0 x1 x2 x3 , (1)试作一个次数不高于 4 次的多项式 H 4 ( x) ,满足条件
H 4 ( x j ) f ( x j ) , j 0,1, 2,3, ' ' H 4 ( x1 ) f ( x1 ) ;
(D)、它不是插值型求积公式。
y f x, y 1、求解常微分方程初值问题 的改进的欧拉法是 y x0 y0
阶方法。 ,其 收 敛 阶
2 、 解 非 线性方程 f x 0 的单根 的牛顿法 格 式为 为 。
4 1 2 2 , x 1 , 则 Ax 3 2 1 1 3 3、设矩阵 A 1 2
(2)写出 E( x) f ( x) H 4 ( x) 的表达式。
四、 (10 分)求 f ( x) e x 在 [0,1] 上的二次最佳平方逼近,权为 1。
3
五、 (10 分)用 n 2, 3 的高斯-勒让德公式计算积分 e x sin xdx
1
六、 (10 分)已知
0 1 2 1 , b , 2 5 8 1 请用 Doolittle 三角分解法求解线性方程组 Ax b 。 1 3 7 7 1 3 9 9
,对于其产生的数列 k 0 , 1, 2 , xk ,下列说法正确的是
(A)、若数列 xk 收敛,则迭代函数 x 唯一; (B)、若对于 x a, b , x 1 ,则 xk 收敛; (C)、若对于 x a, b , x 1 ,则 xk 收敛; (D)、若对于 x a, b , x L 1 ,则 xk 收敛。 3、对矩阵 A 采用幂法迭代,如果该方法收敛,则其收敛速度取决于( (A)、模最大特征值和模次最大特征值的比值; (B)、模最大特征值和模次最大特征值的模的比值; (C)、模次最大特征值和模最大特征值的比值; (D)、模次最大特征值和模最大特征值的模的比值。 ) 。

2011年下学期数值分析考试试卷答案(A)

2011年下学期数值分析考试试卷答案(A)

2011年下学期数值分析考试试卷答案(A)D222223221()()(1)(2)(1)21(45)2P x H x Ax x x x x x x x x =+-=-+-=-+余项为 R(x)=(5)22()(1)(2)5!f x x x ξ-- ……………………………12分解法2:构造带重节点的Newton 差商表 0 0 0 0 0 1 1 1 1 1 1 1 0 -1 2211/2 ………………………8分2222221()00(0)1(0)1(0)(1)(0)(1)21(45)2N x x x x x x x x x x =+-+----+--=-+…………………12分三、 (12分) 求()xf x e -= 在区间[-1,1]上的最佳平方逼近2次多项式. (用勒让德正交多项式2121{(),(),()}{1,,(31)}2P x P x P x x x =-) 解:用勒让德多项式20121{(),(),()}{1,,(31)}2P x P x P x x x =-,2(,)21iiP P i =+ …………………………………………………………………………………..3分计算:11101(,)( 2.3504)x f P e dx e e ---==-≈⎰,1111(,)20.7358x f P xe dx e---==-≈-⎰121211(,)(31)70.143132x f P x e dx e e ---=-=-≈⎰…………………………………………………………………………………..8分111101010011(,)(,)2* 1.1752,*3 1.1036(,)2(,)2/3 f P f P e e e a a e P P P P ----==≈==-=-≈-12222(,)7*0.3578(,)2/5f P e e a P P --==≈故最优平方逼近函数为:11112225351()3(31)22211.1752 1.10360.3758(31)20.5367 1.10360.9963e e e e p x e x x x x x x -----=-+⋅-≈-+⋅-=-+。

2010-2011学年第一学期末数值分析考试试题A参考答案与评分标准

2010-2011学年第一学期末数值分析考试试题A参考答案与评分标准

2010-2011学年第一学期末数值分析考试试题A参考答案与评分标准D第2页,共5页第3页,共5页第4页,共5页()()()1121212,,n n n n n n h y y k k k f t y k f t h y hk +⎧=++⎪⎪=⎨⎪=++⎪⎩,它是 2 阶方法。

二、解答下列各题(每小题12分,共24分) 1. 用LU分解法求解线性方程组1234230156715112133721182x x x x -⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪=⎪ ⎪ ⎪- ⎪ ⎪ ⎪--⎝⎭⎝⎭⎝⎭;解:() 2 3 0 1 -53 -2 1 2 41 2 1 -2 41 2 -1 2 1 -11 -2 2 -3015671511|2 1 133 7201182LU A b ⎛⎫-⎛⎫ ⎪⎪ ⎪-⎪ ⎪=−−−→ ⎪- ⎪ ⎪ ⎪--⎝⎭ ⎪⎝⎭分解1 0 0 023 0 13 1 0 0 0 -2 1 21 2 1 0 0 0 1 -21 2 -1 1 0 0 0 154,,41L U y -⎛⎫⎛⎫⎛ ⎪ ⎪⎪ ⎪ === ⎪ ⎪ ⎪ ⎪ -⎝⎭⎝⎭⎝ 1-2 2-1,x ⎫⎛⎫⎪ ⎪⎪ ⎪=⎪ ⎪⎪ ⎪⎭⎝⎭.. LU 矩阵(或对应元素每算对两个给1分) 2. x19 25 31 38 44 y 19.0 32.3 49.0 73.3 97.8 求形如2y a bx =+的拟合函数及误差平方和,并估计变量y 在21x =处的值。

第5页,共5页解:令a cb ⎛⎫= ⎪⎝⎭,1361162519611144411936A ⎛⎫ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭, ........................... 3分对应的正规方程组TTA Ac A y =为55327271.453277277699369321.5a b ⎛⎫⎛⎫⎛⎫= ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭....................... 7分 解此方程得0.972578660.05003512a b ⎛⎫⎛⎫= ⎪ ⎪⎝⎭⎝⎭...................................... 9分 即最终的拟合函数 (2)0.972578660.05003512y x =+ 10分误差平方和为 ................................ ()()4220.1502kkk y x y σ==-=∑ 11分2(21)0.972578660.050035122123.0380683y =+⨯= ....... 12分三(本题12分)(1)Gauss 求积公式(()11585()d 0.600.6999f x x f f f -≈-++⎰的代数精度是多少?(2)利用上述公式计算定积分111d 3x x-+⎰,并与真值进行比较;0.60.7745967≈)(3)将区间[1,1]-等分成两个小区间,对每个小区间分别利用(1)中的Gauss 求积公式进行计算,并将最后得到的111d 3x x-+⎰的近似值与(2)进行比较。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
解:由 x( k 1) Mx( k ) g 和 x* Mx* g 可得:
x( k 1) x* M ( x( k ) x* ) , k 0,1,2,...
递推的: x( k ) x* M k ( x(0) x* ) 设 y 是矩阵 M 属于特征值 的特征向量,取 x(0) y x* ,则有:
3 5x f ( xk ) xk a a , xk 1 k 2 , k 0,1,2,... xk xk 2 6 f ( xk ) 6 xk 6 xk
一、解答下列各题: (每题 5 分,共 30 分) 1.设近似值 x 具有 5 位有效数字,则 x 的相对误差限为多少? 解:记 x 0.a1a2 ...10 ,则 x 的相对误差为:
1 1
1 0 0 0 0.3 0.2 0 0.3 0.2 0 0. 4 G ( D L) U 0 1 0 0 0 0.4 0 1 0 1 0 0 0 0 0. 3 0. 2
h2 h3 y ( xn ) ( y ( xn ) O(h 4 ) 2 6
y( xn1 ) y ( xn ) hy ( xn )
yn f n h
h 2 f n f n h3 ( fn ) y ( xn ) O(h 4 ) 2 x y 6
2
… … … … ○ … … … … 密 … … … … ○ … … … … 封 … … … … ○ … … … 线 … … … … … … … …
1 三、 (9 分)说明方程 2 x sin x 2 0 在区间 [ , ] 内有唯一根,并建立一个 2 2 1 收敛的迭代格式,使对任意初值 x0 [ , ] 都收敛,说明收敛理由和收敛阶。 2 2 1 1 解:记 f ( x) 2 x sin x 2 ,则 f ( x) C[ , ] ,且 f ( ) 0, f ( ) 0 ,而且, 2 2 2 2 1 ] 内有唯一根。 f ( x) 2 cos x 0 ,所以,方程 2 x sin x 2 0 在区间 [ , 2 2 1 建立迭代格式: xk 1 sin xk 1, k 0,1,2,... 2 1 1 由于,迭代函数 (x) sin x 1 在区间 [ , ] 上满足条件: 2 2 2 1 3 1 1 1 ( x) , | ( x) || cos x | 1 2 2 2 2 2 1 所以,此迭代格式对任意初值 x0 [ , ] 都收敛。 2 2 1 又由于, ( ) cos 0 ,所以,此迭代格式 1 阶收敛。 2
五、 (4 分)设矩阵 M 是 n 阶方阵, M 有一个绝对值小于 1 的特征值 ,且方程 组 x Mx g 有 唯 一 解 x * , 证 明 : 存 在 初 始 向 量 x ( 0 ) 使 迭 代 格 式 :
x ( k 1) Mx ( k ) g , k 0,1,2,...产生的序列 {x ( k ) } 收敛到 x * .
么? 解:由于 | sin(x 2 y) sin(x 2 y ) | | 2cos( x 2 )( y y ) | 2 | y y | 即,函数 f ( x, y) sin(x 2 y) 连续,且关于变量 y 满足 Lipschitz 条件,所以,改 进 Euler 方法收敛。
* m
0.5 105 x* x 0.5 10m5 0.5 104 m 0.1 x* 0.a1 a 2 ... 10
xk 1
即,相对误差限为: 0.5 10 .
4 2 0 2.问 a , b 满足什么条件时, 矩阵 A 2 5 a 有分解式 A GGT , 并求 a b 2 时 0 b 5
x1 0.3 x 2 0.2 x3 1 1.用 Gauss-Saidel 迭代法求解方程组 x 2 0.4 x3 2 ,如果取初值 x x 1 3 1
1 1 所以, H 3 ( x) ( x 2)( x 2 2 x 1) ( x 3 3x 2) 。 2 2 1 1 4. 确定求积公式 f ( x)dx f (1) A1 f (0) A2 f (1) 中的待定系数, 使其代数精 1 2
(k ) (k ) x1( k 1) 0.3x 2 0.2 x3 1 ( k 1) (k ) 0.4 x3 2 由于 Gauss-Saidel 迭代格式为: x 2 ,所以, x ( k 1) x ( k 1) 1 1 3源自解: I S 2

G
10
1 G
x
(1)
x
0.510 1 1 2 8 0.00390625 0.5 256 2
5M 4
2880 2 4
=0.006641
2.给定离散数据 xi yi -1 2 0 -1 1 1 2 3
y sin(x 2 y) , 0 x 2 6.求解初值问题 的改进 Euler 方法是否收敛?为什 y(0) 1
0
0
二、解答下列各题: (每题 8 分,共 48 分)
2c 1 , a b c 2 , a c 0 ,解得:
a c 1 / 2, b 1,
… … … … ○ … … … … 密 … … … … ○ … … … … 封 … … … … ○ … … … 线 … … … … … … … …
f [1,2,3,4] 4 , f [1,2,3,4,5] 0
6. 设 p2 ( x) 是 区 间 [0, 1] 上 权 函 数 为 x 的 二 次 正 交 多 项 式 , 计 算 积 分
x 2 x2 2 3.解线性方程组 1 的 Jacobi 迭代法是否收敛,为什么? 2 x1 9 x2 3
试求形如 y a bx2 的拟合曲线。 解:由于 0 ( x) 1,1 ( x) x 2 ,所以 0 (1,1,1,1)T ,1 (1,0,1,4)T ,
4a 6b 5 , f (2,1,1,3)T ,所以,正则方程组为: 6a 18b 15






… … … … ○ … … … … 密 … … … … ○ … … … … 封 … … … … ○ … … … 线 … … … … … … … …
东 北 大 学 研 究 生 院 考 试 试 卷 2011 —2012 学年第 数值分析 一 学期
总分





课程名称:
(共 3 页)
4.对方程 f ( x) ( x 3 a) 2 0 建立 Newton 迭代格式,并说明此迭代格式是否收 敛?若收敛,收敛阶是多少? 解:Newton 迭代格式为:
f( 1 ) f( 0 ) 2 (5) 7 1 0
所以,当 2 5 a b 2 5 时有分解式 A GGT , a b 2 时有:
4 2 0 2 0 0 2 1 0 A 2 5 2 1 2 0 0 2 1 0 2 5 0 1 2 0 0 2
x( k ) x* M k y k y ,于是有: x ( k ) x*
k
k
y
所以, lim x ( k ) x* 0 ,即序列 {x ( k ) } 收敛到 x * .
y f ( x, y) , x [a, b] 四、 (9 分) 已知求解常微分方程初值问题: 的差分公式: y ( a)
1
1 1 f (1) f (0) f (1) 代数精度最高. 2 2
又由于公式对 f ( x) x 2 不能精确成立,所以,代数精度为 1,不是插值型求 积公式。 5.利用复化 Simpson 公式 S 2 计算定积分 I sinxdx 的近似值,并估计误差。
0
所以, G

0.5 ,
h h y n 1 y n hf ( x n , y n f ( x n , y n )) 2 2 y 0
求此差分公式的阶。 解:由于
y n1
h 2 f n f n h3 2 f n 2 f 2 f yn f n h ( fn ) ( 2 2 f n 2 f n2 ) O(h 4 ) 2 x y 8 x xy y

12
[sin 0 sin 2 sin

2
4 sin

4
4 sin
3 ] (1 2 2 ) 2.00456 4 6
x (1) (1,2,2)T , x (1) x ( 0)
x10 x
*

2 ,于是
( 0)
由于 f ( x) sin x 的 4 阶导数在 [0, ] 上的最大值为: M 4 1,所以 误差为: | I S 2 |
4
由于迭代函数为: ( x)
5x a 2 ,方程根为: 3 a ,所以, 6 6x
( )
1 5 a 1 3 1 ,且 ( ) 0 2 6 3 2
所以,此迭代格式收敛,收敛阶是 1.
的分解式(其中 G 是对角线元素大于零的下三角形矩阵).
x0 (0,0,0)T ,试估计迭代 10 步的误差 x10 x *
解:由于 Gauss-Saidel 迭代矩阵为:
1

.
度尽可能高,并问此公式是不是插值型求积公式. 解:令公式对 f ( x) 1, x 都精确成立,得: A1 A2 3 / 2, A2 1 / 2 , 所以, A1 1, A2 1 / 2 时,公式 f ( x)dx
相关文档
最新文档