硕士研究生数值分析试卷
研究生《数值分析》试卷
研究生“数值分析”试题一, 填空(20分)1,n +1个互异节点插值型数值求积公式的代数精度为________次,最高为________次。
2,SOR 方法收敛的必要条件:松弛因子ω满足条件_________。
3,对于插值型求积公式∑⎰=-≈nk k k x f A dx x f 011)()(,其节点),,1,0(n k x k =是高斯点的充分必要条件是_________。
4,设)(ij a A =为n ×n 矩阵,则1A =________,∞A =________。
5,设解方程组b Ax =的迭代法为d Bx x k k +=+1,则迭代收敛的充分必要条件是________。
6,判断下面的函数是否为三次样条函数(填是或否)(1)211001)1(0)(233≤≤<≤<≤⎪⎩⎪⎨⎧-+=x x x x x x x f - (2)⎩⎨⎧≤≤<≤-++++=100112212)(33x x x x x x x f二,(10分)在22-≤≤-x 上给出x e x f -=)(等距节点函数运用二次插值求x e -的近似值,要使误差不超过610-,问使用函数表的步长应取多大?三,(10分)四,(10分)设)(x f 在[]30,x x 上有三阶连续导数,且3210x x x x <<<,试作一个次数不高于四次的多项式)(x p ,满足条件)()(j j x f x p ==j 0,1,2,3)(')('11x f x p = 推导它的余项)()()(x p x f x E -=的表达式五,(10分)试用Romberg (龙贝格)方法,计算积分⎰311dx x,并精确到小数点后4位。
六,(10分)利用数值积分的Simpson (辛甫生)公式,导出公式)''4'(31111-+-++++=n n n n n y y y h y y 并指出次方法的阶七,(10分)设0)(=x f 的单根α,)(x F x =是0)(=x f 的等价方程,则:)(x F 可表为)()()(x f x m x x F -=证明: 当1)]('[)(-≠ααf m 时,)(x F 是一阶的。
研究生数值分析试题
一、选择题(四个选项中仅有一项符合题目要求,每小题 3 分,共计 15 分)
1、一般用高斯消元法解线性代数方程组要采用的技术是(
)
(1)调换方程位置; (2)选主元; (3)直接求解; (4)化简方程组。
⎛ 2 2 3⎞ ⎛ 1 0 0 ⎞⎛2 2 3⎞
2、设矩阵
A
为初值迭代一步。
四、(12 分)应用牛顿法于方程
f (x) =
xn
−a
Байду номын сангаас
=
0和
f (x) =1−
a xn
= 0 ,分别导出求 n
a
的
迭代公式,并求极限 lim n a − xk+1 。 k→∞ ( n a − xk )2
五 、 ( 12 ) 方 程 x3 − 6 x − 8 = 0 在 x = 3 附 近 有 根 , 把 方 程 写 成 三 种 不 同 的 等 价 形 式
零, A = LU 为 Doolitte 分解,则上三角矩阵 U 的上半带宽为
。
5、设对称正定矩阵
A
=
(aij
)∈
Rn×n , a11
≠
0
,经过一次
Gauss
消元得到形如
A
=
⎛ ⎜ ⎝
a11 0
∗⎞
A1
⎟ ⎠
的
矩阵,则 A1 是
矩阵。
三、(12 分)试用高斯列主元素法求解线性方程组
⎡ 1 3 −2 −4 ⎤ ⎡ x1 ⎤ ⎡3 ⎤
3、设矩阵 A ∈ Rn×n , Q ∈ Rn×n ,且 QT Q = E ,则下列关系式不成立的是(
)
(1) A = AQ ;(2) QA = A ;(3) Qx = x ,其中 x ∈ Rn ;
太原理工大学11级硕士研究生数值分析期末考试题
11级(12/07/03)
一、基础题(40分)
(一)、单项选择(2×5=10分)
1、求解常微分方程的预估—校正法的局部截断误差为( )。
2、过
3
4
5、
(二)
1、是一日插值基函数在节点上的取值是______________。
2、设分段多项式,
,
是以0,2,3为节点的三次样条
函数。
则a =____________,b =____________, c =____________。
3、设,则关于节点,,的二阶向前差分为_________。
4、5个节点的牛顿—科特斯求积公式的代数精度为________,5个节点的求积公式最高代数精度为________。
5、设,则a的取值范围为________A可分解为A = LL T,且当L满足________,分解是唯一的。
6、设是切比雪夫正交多项式系,则的正交区间为________,它的权
7迭8。
1
2
(1)
(2)
3、用二步法求解一阶常微分方程初值问题
,
,问:如何选择、的值,才能使该方法的阶数尽可能高?写出此时的局部截断误差主项。
三、计算题(15×2=30分)
1、(1)设,,是区间[-1,1]上权函数为的最高项系数为1的正交多项式组,其中,,求。
(2)利用,,求函数在[-1,1]上的二次最佳平方逼近多项式。
2、已知求解方程组Ax = b的分量迭代格式:,
,,,;,,,
(1)试求出矩阵格式及迭代矩阵。
(2)证明当A为严格对角占优矩阵,时,该迭代格式收敛。
数值分析(研究生)试卷
华中科技大学研究生课程考试试卷课程名称:_______________________ 课程类别考核形式数值分析学生类别______________考试日期______________学生所在院系_______________ □公共课□专业课√□开卷□√闭卷2009.5.6学号__________________姓名__________________任课教师___________________ 一、填空题(每空2分,共20分) 1. 为避免有效数字的损失,应将,1,ln )1ln(>>-+x x x 改写为_____________。
2. 设其三阶差商,200720082009)(3++=x x x f =]3,2,1,0[f _____________,四阶差商____________。
=]4,3,2,1,0[f 3. 设是上带权b x x x +-=22)(?]1,0[1)(=x ρ的正交多项式,则=b ___________。
4. 对于常微分方程数值解,若某算法的局部截断误差为,则称该算法有_____________阶精度;显式欧拉法有____________阶精度。
)O(h1p+5. 设是的二重根。
*x 0)(=x f )(x f ′′在邻近连续,则用迭代公式________________*x 求此根的近似值所产生的序列至少具有二阶收敛性。
6. ,当a 满足条件___________时,A 可作LU 分解,当a 满足条件__________时,必有分解式,这种分解唯一吗? _____________ ??????+=1221a A TL L A ?=二、(10分)函数在上有三阶连续导数,作一个不高于二次的多项式满足)(x f ],[10x x )(x P .)()(),()(),()(110000x f x P x f x p x f x P =′=′=证明其唯一性,并写出它的余项的表达式。
Ch1-Ch4(2009-2011级硕士研究生《数值分析》试卷)
2009级一、判断题 (每题2分)3. 若n 阶方阵A 是严格对角占优的,则解方程组A =x b 的Jacobi 迭代法收敛。
( √ )4. 设是方程的根,则求的Newton 迭代法至少是平方收敛的。
( ) *x 0)(=x f *x二、填空题 (每空2分)1. 近似数关于准确值* 3.120x = 3.12065x =有 位有效数字,相对误差是 . 4. 设2543A −⎡⎤=⎢⎥−⎣⎦,则1A = ,A ∞= ,1Cond()A = .五(本题满分10分) 对于下列方程组1231231234222633245x x x x x x x x x ,,,−+=⎧⎪++=⎨⎪++=⎩ 建立Gauss–Seidel 迭代公式,写出相应的迭代矩阵,并用迭代矩阵的范数判断所建立的Gauss–Seidel 迭代公式是否收敛。
七(本题满分10分) 已知方程在10x xe −=00.5x =附近有一个实根.*x (1) 取初值00.5x =,用Newton 迭代法求(只迭代两次)。
*x (2) 取初值010.5,0.6x x ==,用弦截法求(只迭代两次)。
*x2010级一、填空题 (每空2分,共20分)1. 近似数关于准确值*2.315x = 2.31565x =有 位有效数字,相对误差是 .4. 设2345A −⎡⎤=⎢⎥−⎣⎦,则1A = ,Cond()A ∞= .5. 设是方程的3重实根,则求的改进的Newton 迭代公式为 *x 0)(=x f *x .二 (本题满分8分) 对下列方程组1231231232633245,422x x x x x x x x x ,++=⎧⎪++=⎨⎪−+=⎩ 建立收敛的Jacobi 迭代公式和收敛的Gauss–Seidel 迭代公式,并说明理由。
五(本题满分10分) 已知方程在3210x x −−=0 1.5x =附近有一个实根.*x (1) 取初值0 1.5x =,用Newton 迭代法求(只迭代两次)。
硕士研究生数值分析试卷
数值分析(研究生,2008-12-15)1.(10分)求函数⎩⎨⎧≤≤++<≤-+=10,101,1sin )(2x x x x x x f 在区间[-1,1]上的最佳平方逼近式x e a x a a x 210)(++=φ。
2.(15分)利用乘幂法计算下列矩阵的主特征值和相应的特征向量⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----110141012,初始向量为T x ]0,0,1[0=(要求结果有三位有效数字)。
同时计算该矩阵的1-条件数和谱条件数。
3.(15分)已知函数x x f sin )(=在36.0,34.0,32.0210===x x x 处的值分别为352274.0,333487.0,314567.0210===y y y 。
用Lagrange 插值多项式对3167.0=x 的函数值进行近似计算,并估计近似计算的误差界。
4.(15分)用Newton 迭代法求方程0ln 2=+x x 在区间(0,2π)内的解,选择你认为合适的初始点,计算方程的根,使得近似解具有四位有效数字。
请从理论上估计达到所需精度所需的迭代次数。
5.(15分)用Gauss-Seidel 迭代法解方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---542834*********x x x 取初始近似向量0[0,0,0]Tx =,估计达到4位有效数字需要的迭代次数,并实际计算之。
就该具体问题分析计算过程中总的乘除法计算量。
6. (10分)应用拟牛顿法解非线性方程组⎪⎩⎪⎨⎧=-+=-+.12,2322112221x x x x x x 取T x ]1,0[)0(= ,终止容限210-=ε。
7.(10分) 求解矛盾方程组⎪⎪⎩⎪⎪⎨⎧=++=++=++=++232328.12221321321321321x x x x x x x x x x x x8. (10分)用复合Simpson 公式计算积分⎰=21sin )(xdx f I 讨论在误差要求不超过410-的条件下的步长。
060708研究生数值分析试卷(A).doc
武汉大学2006〜2007学年第一学期硕士研究生期末考试试题(A 卷)科H 名称:数值分析 学生所在院: 学号: 姓名:注意:所有的答题内容必须答在答题纸上,凡答在试题或草稿纸上的一律无效。
一、(12分)设方程组Ax = 0为■1、 (1\J 1>(1)用Doolittle 分解法求解方程组;(2) 求矩阵A 的条件数Cwd(A)g 二、(12分)设A 为n 阶对称正定矩阵,A的n 个特征值为山 < 心< .•. V 九,为 求解方程组Ax = b,建立迭代格式求出常数s 的取 值范围,使迭代格式收敛。
三、(12分)已知数据试用二次多项式p ⑴=ax 1 2+hx + c 拟合这些数据。
四、(14分)已知y = /(x)的数据如下:取得最小值。
六、 (12)确定常数片,使求积公式1求f (x)的Hermite 插值多项式W 3(x);2 为求\\f{x)dx 的值,采用算法:•⑴必:=「久3)击+ R 试导出截断误差R五、(12分)确定常数。
,b 的值,使积分r I.2I(a,b) = J 0(czx + /?-/) dxc 2^f{x)dx a A/(0) + A2/(l) + A3/(2)的代数精度尽可能高,并问是否是Gauss型公式。
七、(12分)设伊⑴导数连续,迭代格式x M =(p{x k)—阶局部收敛到点x*。
对于常数人,构造新的迭代格式:A 1 ,、队=一从+ 一心)1 +2 1 + 人问如何选取人,使新迭代格式有更高的收敛阶,并问是儿阶收敛。
八、(14分)对于下面求解常微分方程初值问题」方= 的单步法:Mo) = JoA)'〃+】=儿 + hk2< k、=(1)验证它是二阶方法;(2)确定此单步法的绝对稳定区域。
武汉大学2007~2008学年第一学期硕士研究生期末考试试题科目名称:数值分析学生所在院:学号:姓名:注意:所有的答题内容必须答在答题纸上,凡答在试题或草稿纸上的一律无效。
【研究生】数值分析试题
文档来源为:从网络收集整理.word 版本可编辑.欢迎下载支持. 1文档来源为:从网络收集整理.word 版本可编辑.华南理工大学研究生课程考试 《数值分析》试卷A (2015年1月9日) 1. 考前请将密封线内各项信息填写清楚; 所有答案请按要求填写在本试卷上; 课程代码:S0003004; 4. 考试形式:闭卷; 5. 考生类别:硕士研究生; 本试卷共八大题,满分100分,考试时间为150分钟。
一.(12分)解答下列问题 1.欲计算下式: ()13(1)2(1)(2)7(1)(2)(3)6(1)(2)(3)(4),P x x x x x x x x x x x =+-+------+---- 2.设有递推公式 0161,1,2,n n y y y n -⎧=⎪⎨=-=⎪⎩ *001.732y y = 作实际计算,问计算到10y 时误差为初始误差*00y y -的多少 这个计算过程数值稳定吗 ? . (14分)解答下列问题 1. 若2()63f x x =+,则[1,2,3]f 和[1,2,34]f ,的值分别是多少? 2. 1012 . (10分) 设f 在互易节点i x 上的值()()0,1,....i i f f x i n ==。
试证明:f 在节点i x 上n 次最小二乘拟合多项式()n p x 与f 在节点i x 上的n 次Lagrange 插值多项式()n L x 一致,()()=n n p x L x 。
. (12分) 按代数精度的定义,构造下列形式的求积公式(即确定参数,A B ,α): Gauss 型求积公式。
. (14分) 已知线性代数方程组Ax=b 为: ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-----n n n n n n n n b b b b x x x x d u u u v d v d v d 121121121112211000000文档来源为:从网络收集整理.word 版本可编辑.欢迎下载支持.2文档来源为:从网络收集整理.word 版本可编辑. (1) 用顺序高斯消去法求解方程组Ax=b ;(2) 先由(1)的消元过程直接写出A 的LU 分解,再利用该LU 分解求解方程组Ax=b 。
研究生考试数值分析试题
研究⽣考试数值分析试题研究⽣2002级数值分析⼀(12分)、对于积分=+1,2,1,0,999n dx x x n。
(1)试推导递推公式 ,2,1,19991=+-=-n nI I n n ;(2)分析上述算法的数值稳定性;(3)若上⾯算法不稳定,请选择合适的算法,并分析其稳定性。
⼆(12分)、解⽅程组= 00001.8800001.626221x x 和?=00002.8800001.626221x x ,就所观察到的现象进⾏分析。
三(12分)、设⽅程组=--=+-=+-7989783212121x x x x x x x ;(1)适当调整⽅程的排列顺序,使得⽤Gauss-Seidel 迭代法求解时收敛?说明收敛原因。
(2)取初始向量()()Tx 0,0,00=,⽤Gauss-Seidel 迭代求近似解()2x,并求其()()k k x x-+1误差。
四(12分)、(1)已知函数()4xe xf =,在[0,1]内三点0,1/2,1的函数值,求其⼆次插值的余项;(2)三个节点如何安排能使其余项达最⼩,此时⼈余项为多少?五(12分)、对于⽅程()02ln =+-x x ,若求[-1.9,-1]内的根,分别选取迭代⽅程()2ln +=x x 和2-=x e x ,它们的收敛性如何?再写出⽜顿迭代公式。
六(10分)、设()?=>+-='100,5y x x y y ,解析解xe x y -+-=25262515,分别取45.0,4.0,2.0,1.0=h ,利⽤Euler ⽅法计算得y(10)的近似值分别为1.96,1.96,5.2851,142.8863,对此现象进⾏分析。
七(10分)、设()x e x f =,分别取步长0001.0,01.0,5.0=h ,⽤中⼼差商公式计算()0f '的近似值并求出误差,对结果作分析⽐较。
⼋(10分)、求不超过2次的多项式()x P 2,使其满⾜条件:()21=f ,()32=f ,()12='f ,并写出其误差估计。
2010-2011学年北京理工大学硕士研究生数值分析期末试卷
T
1 2 Байду номын сангаас 2
1 2 3 3
1 2 ,b=(4,3,2,1) 3 4
1. 用平方根法解线性方程组。 2. 对上述方程组构造收敛的迭代格式,说明其收敛原因,取初始值 X(0)=(0,0,0,0)T 用所给的迭代格式计算迭代序列的前两项(用分数表示) 。 二(15 分)已知sin (0.32)=0.314567,sin (0.34)=0.333487 均具有 6 位有效 数字。 1. 请用线性插值求sin (0.33)的近似值。 2. 证明在区间[0.32,0.34]上用线性插值求sin x的近似值时至少有 4 位有效数字。 三(20 分)长半轴为 2,短半轴为 1 的椭圆的周长 s 为 s=8
h
π 2 0
1−
3 cos x 4
2
dx,
用数值积分的方法求其近似值(要求计算结果具有四位有效数字) 。 2 四(15 分)用迭代法求 x +10x-18=0 在[1,2]内的根,取初值为 1.5 1. 构造一个收敛的迭代格式,并证明此格式的收敛性。 2. 先用上述迭代格式计算 2 步,然后采用 Aitken 加速算法再计算一步是否能得 到更精确的近似值?计算过程中小数点后保留 4 位。 五(10 分)求函数 ex 在区间[0,1]上的一次最佳平方逼近多项式。 y ′ + y = 0, xϵ[0,1] 六(20 分)对初值问题 y 0 =1 1. 求此微分方程的精确解。 2. 证明:用格式yn+1 = yn + 2 (−yn − yn+1 )所求得的近似解在步长 h0 时收敛 于精确解。 3. 写出上述格式的 Matlab 程序源代码, 要求: 输出近似解曲线图和误差曲线图。
硕士课程—数值分析题集(附答案).docx
2009-2010数值分析第一章绪论 (1)第二章函数插值 (2)第三章函数逼近 (5)第四章数值积分与数值微分 (10)第五章解线性方程组的直接解法 (12)第六章解线性方程组的迭代解法 (16)第七章非线性方程求根 (19)第九章常微分方程初值问题的数值解法 (21)第一章绪论1.1要使胸的相对误差不超过0.1%,应取几位有效数字?解:面的首位数字%=4。
设/有n位有效数字,由定理知相对误差限k(.r*)|<—xlO1^ =-xl0^1 r 1 2x4 84-xio1-" <0.1%, 8解得〃Z3.097,即需取四位有效数字.1.2 序列{/}满足关系式y,,=10y,_]-l(n = l,2,...),若y0=V2«1.41,计算到M。
,误差有多大?这个算法稳定吗?解:y0 = V2,y* =1.41,|y0 -y*| <^-xl0-2=5 ,于是|/i 一川=|1。
》0 —IT。
〉;+1| = 1。
|光 - 司 < 1。
5卜2-》;| = |10》1一1一10》;+1| = 10卜1一酣〈10逆, 一般地|儿一司<103 因此计算到Mo其误差限为1010^,可见这个计算过程是不稳定的。
1. 3计算球的体积,要使相对误差限为1%,问测量半径R时允许的相对误差限是多少?解:5,、九兀K ~-7tK R_R* R2+R*R + R*2R_R* 37?2R_R*。
,“ ,(v)= _2 ---------- 2 «■«.____________ = _____ 3 = 1% ' 4 f RR- R R 2 R-7lR 3》=一' ,即测量半径R 时允许的相对误差限是一、。
R 300300第二章函数插值2.1、利用如下函数值表构造差商表,并写出牛顿插值多项式。
进而得牛顿多项式为 地⑴=f (.%) + /■氏次』吼⑴+ /[.r (p x 1,.r 2]<»2(.r) + /[.r (p x 1,.r 2,.r 3]<»3(.r)1 1 33A^3 (x) = 3 + — (x -1) + — (x -1)(尤)-2(x- l)(x )x2. 2、已知f(-2) = 2, f(-1) = 1, f (0) = 2, f (0.5) = 3试选用合适的插值节点利用Lagrange 二次插值多项式计算f (-o.5)的近似值,使之精度 尽可能高。
硕士研究生《数值分析》试卷2013(A)
硕士研究生《数值分析》试卷2013(A)一、判断题 (下列各题,你认为正确的,请在题后的括号内打“√ ”,错误的打“×”,每题2分,共10分) 1. 近似数*3.200x =关于准确值 3.200678x =有4位有效数字。
( ) 2. 设(0,1,2,3)i x i =是互异的点,()(0,1,2,3)i l x i =是Lagrange 插值基函数,则3224()4i ii x l x x==∑. ( )3. 设73()32f x x x =-+,则差商1234567[2,2,2,2,2,2,2]1f =。
( ) 4. 设A 是n 阶非奇异方阵,则解方程组A =x b 的迭代法收敛的充要条件是A 的谱半径()1A ρ<。
( )5. 解常微分方程初值问题的四阶Runge-Kutta 方法的整体截断误差是4()O h ,其中h 是步长。
( )二、填空题 (每空2分,共16分) 1. 设T(2,1,3,4)=-x ,2543A -⎡⎤=⎢⎥-⎣⎦. 则 1||||x = , Cond()A ∞= .2. 设20()d I f x x =⎰,若用梯形求积公式计算I ,结果是4;用Simpson 求积公式计算I ,结果是2. 则(1)f = .3. 设S 是函数f 在区间[0,3]上满足第一类边界条件的的三次样条:()()22,01,()111,13,2x x S x x a x b x ⎧≤≤⎪=⎨-+-+≤≤⎪⎩ 则a = ,b = ,(3)f '= .4. 设函数(0.8) 1.2,(0.9) 1.4,(1) 1.0,(1.1)0.2,(1.2)0.5f f f f f =-=-=-==, 步长0.2h =,则用三点数值微分公式计算(1)f '的近似值为 .5. 设函数()f x 是最高次项系数为1-的3次多项式,2()p x 是()f x 在节点1,0,1-上的Lagrange 插值多项式, 则余项2()()f x p x -= .三(本题满分8分)的近似值*x 的相对误差限是0.01%,求*x 至少应具有几位有效数字?四(本题满分10分) 对下列方程组分别建立收敛的Jacobi 和Gauss-Seidel 迭代格式,并说明理由。
河海大学硕士数值分析试卷08~09
1. 填空(10分,每空2分)1) 为了减少运算次数,应将表达式113141817162345---++x x x x x 改写为 .为了减少舍入误差的影响,应将表达式20002001-改写为 .2) 用二分法求方程0152)(3=--=x x x f 在区间[1,3]内的根,进行一步后根所在区间为,进行二步后根所在区间为 .3) 在高斯顺序消去法中,)1(-k kk a 称为第k 步主元.为使消去法得以顺利进行,必须 .2. 选择题(10分,每题2分)(1)设有求方程1=xxe 根的迭代公式kx k ex -+=1,取初值5.00=x ,则迭代公式A. 发散B. 敛散不定C. 收敛D. 不确定 (2)辛普森(Simpson)公式)]()2(4)([6)(b f ba f a f ab dx x f ba+++-≈⎰可由 A. 分段线性插值导出 B. 抛物插值导出 C. 线性插值导出 D. 分段抛物型插值导出(3)矩阵A 满足什么条件时,LR A =且分解唯一(L 是单位下三角阵,R 是上三角阵)A. 无限制B. 对称C. 可逆D. 严格对角占优 (4)为什么要把解常微分方程的欧拉(Euler)方法发展为改进的欧拉方法?A. 提高精度B. 便于计算C. 提高精度和便于计算D. 稳定性需要 (5) 当0)(=x f 有m 重根时,牛顿(Newton)迭代公式中的迭代函数应为 A. )()()(x f x f x x '-=ϕ B. )()()(x f x f m x x '-=ϕ C. )()()(x f m x f x x '-=ϕ D. )()()(x f x f mx x '-=ϕ4. (10分)选取求积公式中的常数a ,使)]()0([)]()0([2)(0h f f a h f f hdx x f h'-'++≈⎰的代数精度尽量高。
试指出最高代数精度,并估计误差。
060708研究生数值分析试卷(A).doc
武汉大学2006〜2007学年第一学期硕士研究生期末考试试题(A 卷)科H 名称:数值分析 学生所在院: 学号: 姓名:注意:所有的答题内容必须答在答题纸上,凡答在试题或草稿纸上的一律无效。
一、(12分)设方程组Ax = 0为■1、 (1\J 1>(1)用Doolittle 分解法求解方程组;(2) 求矩阵A 的条件数Cwd(A)g 二、(12分)设A 为n 阶对称正定矩阵,A的n 个特征值为山 < 心< .•. V 九,为 求解方程组Ax = b,建立迭代格式求出常数s 的取 值范围,使迭代格式收敛。
三、(12分)已知数据试用二次多项式p ⑴=ax 1 2+hx + c 拟合这些数据。
四、(14分)已知y = /(x)的数据如下:取得最小值。
六、 (12)确定常数片,使求积公式1求f (x)的Hermite 插值多项式W 3(x);2 为求\\f{x)dx 的值,采用算法:•⑴必:=「久3)击+ R 试导出截断误差R五、(12分)确定常数。
,b 的值,使积分r I.2I(a,b) = J 0(czx + /?-/) dxc 2^f{x)dx a A/(0) + A2/(l) + A3/(2)的代数精度尽可能高,并问是否是Gauss型公式。
七、(12分)设伊⑴导数连续,迭代格式x M =(p{x k)—阶局部收敛到点x*。
对于常数人,构造新的迭代格式:A 1 ,、队=一从+ 一心)1 +2 1 + 人问如何选取人,使新迭代格式有更高的收敛阶,并问是儿阶收敛。
八、(14分)对于下面求解常微分方程初值问题」方= 的单步法:Mo) = JoA)'〃+】=儿 + hk2< k、=(1)验证它是二阶方法;(2)确定此单步法的绝对稳定区域。
武汉大学2007~2008学年第一学期硕士研究生期末考试试题科目名称:数值分析学生所在院:学号:姓名:注意:所有的答题内容必须答在答题纸上,凡答在试题或草稿纸上的一律无效。
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分)。
数值分析试题及答案
数值分析试题及答案一、单项选择题(每题3分,共30分)1. 线性代数中,矩阵A的逆矩阵记作()。
A. A^TB. A^-1C. A^+D. A*答案:B2. 插值法中,拉格朗日插值多项式的基函数是()。
A. 多项式B. 指数函数C. 正弦函数D. 余弦函数答案:A3. 在数值积分中,梯形规则的误差是()阶的。
A. O(h^2)B. O(h^3)C. O(h)D. O(1/h)答案:A4. 求解线性方程组时,高斯消元法的基本操作不包括()。
A. 行交换B. 行乘以非零常数C. 行加行D. 行除以非零常数答案:D5. 非线性方程f(x)=0的根的迭代法中,收敛的必要条件是()。
A. f'(x)≠0B. f'(x)=0C. |f'(x)|<1D. |f'(x)|>1答案:C6. 利用牛顿法求解非线性方程的根时,需要计算()。
A. 函数值B. 函数值和导数值C. 函数值和二阶导数值D. 函数值、一阶导数值和二阶导数值答案:B7. 矩阵的特征值和特征向量是()问题中的重要概念。
A. 线性方程组B. 特征值问题C. 线性规划D. 非线性方程组答案:B8. 在数值分析中,条件数是衡量矩阵()的量。
A. 稳定性B. 可逆性C. 正交性D. 稀疏性答案:A9. 利用龙格现象说明,高阶插值多项式在区间端点附近可能产生()。
A. 振荡B. 收敛C. 稳定D. 单调答案:A10. 雅可比迭代法和高斯-塞德尔迭代法都是求解线性方程组的()方法。
A. 直接B. 迭代C. 精确D. 近似答案:B二、填空题(每题4分,共20分)11. 线性代数中,矩阵A的行列式记作________。
答案:det(A) 或 |A|12. 插值法中,牛顿插值多项式的基函数是________。
答案:差商13. 在数值积分中,辛普森规则的误差是________阶的。
答案:O(h^4)14. 求解线性方程组时,迭代法的基本思想是从一个初始近似解出发,通过不断________来逼近精确解。
研究生数值分析期末考试试卷参考答案
研究生数值分析期末考试试卷参考答案太原科技大学硕士研究生2012/2013学年第1学期《数值分析》课程试卷参考答案一、填空题(每小题3分,共30分)1、x x ++11;2、2;3、20;4、6;5、kk k k k x x x x x cos 11sin 1----=+ ( ,1,0=k ); 6、12121)(2++=x x x f ;7、311+=+k k x x ( ,1,0=k );8、12-n ;9、2; 10、+++++++--100052552452552052552525524;二、(本题满分10分)解:Gauss-Seidel 迭代方法的分量形式为+--=+--=++-=++++++3221522)1(2)1(1)1(3)(3)1(1)1(2)(3)(2)1(1k k k k k k k k k x x x x x x x x x -----5分取初始向量T x )0,0,0()0(=时,则第一次迭代可得===315)1(3)1(2)1(1x x x ,--------------7分答案有错误第二次迭代可得=-==7119)2(3)2(2)2(1x x x ,-----------9分所以T x )7,11,9()2(-=.---------------10分三、(本题满分10分)解:构造正交多项式:取)()()()(,)(,1)(01112010x x x x x x x ?β?α?α??--=-==,1)()(402040200=∑∑===i i i i i x x x ??α,1)()(402140211=∑∑===i i i i i x x x ??α,2)()(402040211=∑∑===i i i i x x ??β;所以点集{}1,0,1,2,3-上的正交多项式为12)(,1)(,1)(2210--=-==x x x x x x .-------------------------5分则矩阵???????? ?-----=221111*********A , ??=14000100005A A T ,????? ??=3915y A T ;法方程=????? ??????? ??391514000100005210c c c ----------------8分解得===1431093210c c c ;--------9分所以要求的二次多项式为35667033143)12(143)1(109322++=--+-+=x x x x x y .-----------10分四、(本题满分10分)解:取基函数210)(,1)(x x x ==??,则1),(1000=?=dx ??,31),(10201=?=dx x ??, 51),(10411=?=dx x ?? ππ?2sin ),(100=?=xdx f , 3102141sin ),(πππ?-=?=xdx x f----------------------------------6分法方程-=???? ???????? ??34125131311πππb a -----------------8分解得-=+=33454151543ππππb a .---------------9分所以最佳平方逼近多项式233)45415(1543)(x x ππππ?-++=.---------10分五、(本题满分10分)解:在区间[]1,+n n x x 上对微分方程),(y x f dxdy =进行积分得 ??=++11),(n n n n x x x x dx y x f dx dxdy 即=-+n n y y 1?+1),(n n xx dx y x f -------2分对上式等号右边的积分采用梯形公式进行求解,即+1),(n n x x dx y x f []n n f f h +=+12-------5分所以原微分方程初值问题的数值求解公式为11()2n n n n h y y f f ++=++.-------6分上述数值求解公式的截断误差为 ))](,())(,([2)()(1111n n n n n n n x y x f x y x f h x y x y R +--=++++---8分而又由泰勒公式得)()()()(2'1h O x hy x y x y n n n ++=+;)())(,())(,(11h O x y x f x y x f n n n n +=++;所以))](,()())(,([2)()()()(2'1n n n n n n n n x y x f h O x y x f h x y h O x hy x y R ++--++=+ )()())(,()(22'h O h O x y x hf x hy n n n =+-= 故该方法是一阶的方法.-----------------10分六、(本题满分20分)解:(1)构造的差商表如下:x )(x f 一阶差商二阶差商三阶差商 1 22 4 23 5 1 21- 4 8 3 121 -----------------------------15分(2)取2、3、4作为插值点,----------------------------------------------------17分构造的二次牛顿插值多项式为84)3)(2()2(4)(22+-=--+-+=x x x x x x P -----19分所以25.6)5.3()5.3(2=≈P f .------------------------------20分七、(本题满分10分)解:由泰勒公式可得)2)(()2()('b a x f b a f x f +-++=ξ,),(b a ∈ξ. 把上式代入积分公式?b a dx x f )(可得dx b a x f b a f dx x f b a b a+-++=?)2)(()2()('ξ ?+-++-=b a dx b a x f b a f a b )2)(()2()('ξ 故求积公式的截断误差表达式为?+-b a dx b a x f )2)(('ξ,),(b a ∈ξ.-----------5分当1)(=x f 时,求积公式左边=右边=a b -.当x x f =)(时,求积公式左边=右边=222a b -. 当2)(x x f =时,求积公式左边=333a b -,右边=()()92a b a b +-,左边≠右边. -----8分所以求积公式具有一次代数精度.-------------------------- -----10分。
研究生数值分析试卷
2005~2006学年第一学期硕士研究生期末考试试题(A 卷)科目名称:数值分析 学生所在院: 学号: 姓名: 注意:所有的答题内容必须答在答题纸上,凡答在试题或草稿纸上的一律无效。
一、(15分)设求方程 0cos 2312=+-x x 根的迭代法k k x x cos 3241+=+(1) 证明对R x ∈∀0,均有*lim x x k k =∞→,其中*x 为方程的根.(2) 此迭代法收敛阶是多少? 证明你的结论.二、(12分)讨论分别用Jacobi 迭代法和Gauss-Seidel 迭代法求解下列方程组的收敛性。
⎪⎩⎪⎨⎧=++-=++=-+.022,1,122321321321x x x x x x x x x三、(8分)若矩阵⎪⎪⎪⎭⎫ ⎝⎛=a a a a A 000002,说明对任意实数0≠a ,方程组b AX =都是非病态的。
(范数用∞⋅)四、(15求)(x f 的Hermite 插值多项式)(3x H ,并给出截断误差)()()(3x H x f x R -=。
五、(10分)在某个低温过程中,函数 y 依赖于温度x (℃)的试验数据为已知经验公式的形式为 2bx ax y += ,试用最小二乘法求出 a ,b 。
六、(12分)确定常数 a ,b 的值,使积分[]dx x b ax b a I 2112),(⎰--+=取得最小值。
七、(14分)已知Legendre(勒让德)正交多项式)(x L n 有递推关系式:⎪⎪⎩⎪⎪⎨⎧=+-++===-+),2,1()(1)(112)()(,1)(1110 n x L n n x xL n n x L x x L x L n n n 试确定两点的高斯—勒让德(G —L )求积公式⎰-+≈112211)()()(x f A x f A dx x f的求积系数和节点,并用此公式近似计算积分⎰=211dx e I x八、(14分)对于下面求解常微分方程初值问题 ⎪⎩⎪⎨⎧==00)(),(y x y y x f dx dy的单步法:⎪⎪⎩⎪⎪⎨⎧++==++=+),(),()2121(121211hk y h x f k y x f k k k h y y n n n n n n (1) 验证它是二阶方法; (2) 确定此单步法的绝对稳定域。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
数值分析(研究生,2008-12-15)
( 分)求函数⎩⎨⎧≤≤++<≤-+=1
0,101,1sin )(2x x x x x x f 在区间☯ , 上的最佳平方逼近式
x e a x a a x 210)(++=φ。
.( 分)利用乘幂法计算下列矩阵的主特征值和相应的特征向量
⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡----110141012,初始向量为T x ]0,0,1[0=(要求结果有三位有效数字)。
同时计算该矩阵的 条件数和谱条件数。
( 分)已知函数x x f sin )(=在36.0,34.0,32.0210===x x x 处的值分别为352274.0,333487.0,314567.0210===y y y 。
用☹♋♑❒♋⏹♑♏插值多项式对3167.0=x 的函数值进行近似计算,并估计近似计算的误差界。
( 分)用☠♏♦♦☐⏹迭代法求方程0ln 2=+x x 在区间( ,2
π)内的解,选择你认为合适的初始点,计算方程的根,使得近似解具有四位有效数字。
请从理论上估计达到所需精度所需的迭代次数。
( 分)用☝♋◆♦♦♏♓♎♏●迭代法解方程组
⎥⎥⎥⎦⎤⎢
⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---542834*********x x x 取初始近似向量0[0,0,0]T
x =,估计达到 位有效数字需要的迭代次数,并实际计算之。
就该具体问题分析计算过程中总的乘除法计算量。
( 分)应用拟牛顿法解非线性方程组
⎪⎩⎪⎨⎧=-+=-+.
12,2322112221x x x x x x 取T x ]1,0[)0(= ,终止容限210
-=ε。
( 分) 求解矛盾方程组
⎪⎪⎩⎪⎪⎨⎧=++=++=++=++2
32328
.12221
321321321321x x x x x x x x x x x x
( 分)用复合 ♓❍☐♦☐⏹公式计算积分
⎰=2
1sin )(xdx f I 讨论在误差要求不超过410-的条件下的步长。
数值分析(研究生, )
( 分)求函数x x f =)(在区间☯, 上关于权函数3
1)(x x w =的最佳逼近多项式.)(2cx bx a x ++=φ
.( 分)利用乘幂法计算下列矩阵的主特征值和相应的特征向量
⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡---110121012,初始向量为T x ]0,0,1[0=(要求结果有三位有效数字)。
同时计算该矩阵的∞ 条件数和谱条件数。
( 分)假定给出函数x x f sin )(=的等距点函数表)11(≤≤-x 。
若用二次抛物线插值求x sin 的近似值,要使截断误差不超过310-,问使用函数表的步长h 应取多少?
( 分)用☠♏♦♦☐⏹迭代法求方程032=-x e x 在区间( , )内的解,选择你认为合适的初始点,计算方程的根,使得近似解具有四位有效数字。
请从理论上估计达到所需精度所需的迭代次数。
( 分)用 迭代法解方程组
⎥⎥⎥⎦⎤⎢
⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---042830261016321x x x 取初始近似向量0[0,0,0]T
x =,估计达到精度达到310-需要的迭代次数,并实际计算之。
就该具体问题分析计算过程中总的乘除法计算量。
( 分)应用拟牛顿法解非线性方程组 ⎪⎩
⎪⎨⎧+=+=++=+32,1222y z x e e x y xyz z xy
取T ]1,1,1[作为初始值,终止容限210-=ε。
( 分) 求数据( , ),( , ),( , ),( , )的最小二乘拟合.)(210x a a x +=φ
( 分)用复合 ♓❍☐♦☐⏹公式计算积分 ⎰=1
0cos )(xdx e f I x 讨论在误差要求不超过210-的条件下的步长,并比较实际计算结果与精确结果。