数值方法试卷及答案

合集下载

数值分析模拟试卷1,2,3

数值分析模拟试卷1一、填空（共30分，每空3分） 1 设⎪⎪⎭⎫⎝⎛-=1511A ，则A 的谱半径=)(a ρ______，A 的条件数)(1A cond =________.2 设,2,1,0,,53)(2==+=k kh x xx f k ，则],,[21++n n n x x x f ＝________,],,[321+++n n n n x x x x f ，＝________.3 设⎪⎩⎪⎨⎧≤≤-++≤≤+=21,1210,)(2323x cx bx x x x x x S ，是以0，1，2为节点的三次样条函数，则b=________,c=________.4 设∞=0)]([k k x q 是区间［0，1］上权函数为x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x q ，则⎰=10)(dx x xq k ________，=)(2x q ________.5 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=11001aaa a A ，当∈a ________时，必有分解式，其中L 为下三角阵，当其对角线元素)3,2,1(=i L ii 满足条件________时，这种分解是唯一的.二、（14分）设49,1,41,)(21023====x x x x x f ,（1）试求)(x f 在]49,41[上的三次Hermite 插值多项式)(x H 使满足2,1,0),()(==i x f x H i i ，)()(11x f x H '='.（2）写出余项)()()(x H x f x R -=的表达式.三、（14分）设有解方程0cos 2312=+-x x 的迭代公式为n n x x cos 3241+=+，（1）证明R x ∈∀0均有∙∞→=x x n x lim （∙x 为方程的根）；（2）取40=x ，用此迭代法求方程根的近似值，误差不超过，列出各次迭代值；（3）此迭代的收敛阶是多少？证明你的结论.四、(16分) 试确定常数A ，B ，C 和，使得数值积分公式有尽可能高的代数精度. 试问所得的数值积分公式代数精度是多少？它是否为Gauss 型的？五、（15分）设有常微分方程的初值问题⎩⎨⎧=='00)(),(y x y y x f y ，试用Taylor 展开原理构造形如)()(11011--++++=n n n n n f f h y y y ββα的方法，使其具有二阶精度，并推导其局部截断误差主项.六、（15分）已知方程组b Ax ＝，其中⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫⎝⎛=21,13.021b A ，（1）试讨论用Jacobi 迭代法和Gauss-Seidel 迭代法求解此方程组的收敛性. （2）若有迭代公式)()()()1(b Axa xxk k k ++=+，试确定一个的取值范围，在这个范围内任取一个值均能使该迭代公式收敛. 七、（8分）方程组，其中，A 是对称的且非奇异.设A 有误差，则原方程组变化为，其中为解的误差向量，试证明.其中1λ和2λ分别为A 的按模最大和最小的特征值.数值分析模拟试卷2填空题（每空2分，共30分）1. 近似数231.0=*x 关于真值229.0=x 有____________位有效数字； 2. 设)(x f 可微，求方程)(x f x =根的牛顿迭代格式是_______________________________________________；3. 对1)(3++=x x x f ，差商=]3,2,1,0[f _________________；=]4,3,2,1,0[f ________；4. 已知⎪⎪⎭⎫⎝⎛-='-=1223,)3,2(A x ，则=∞||||Ax ________________，=)(1A Cond______________________ ；5. 用二分法求方程01)(3=-+=x x x f 在区间[0,1]内的根，进行一步后根所在区间为_________，进行二步后根所在区间为_________________；6. 求解线性方程组⎪⎩⎪⎨⎧=+=+04511532121x x x x 的高斯—赛德尔迭代格式为_______________________________________；该迭代格式迭代矩阵的谱半径=)(G ρ_______________；7. 为使两点数值求积公式：⎰-+≈111100)()()(x f x f dx x f ωω具有最高的代数精确度，其求积节点应为=0x _____ , =1x _____,==10ωω__________. 8. 求积公式)]2()1([23)(30f f dx x f +≈⎰是否是插值型的__________，其代数精度为___________。

昆明理工大学—数值分析各年考试题及答案

昆明理工大学数值分析考试题（07）一．填空（每空3分，共30分）1．设A 0.231x =是真值0.229T x =的近似值，则Ax 有位有效数字。

2．若74()631f x x x x =+++，则017[2,2,...2]f = ，018[2,2,...2]f = 。

3． A=1031⎡⎤⎢⎥-⎣⎦,则1A = ；A ∞= ；2A =2()cond A = 。

4．求方程()x f x =根的牛顿迭代格式是。

5．设105%x =±，则求函数()f x =的相对误差限为。

6．A=2101202a a ⎛⎫ ⎪ ⎪ ⎪⎝⎭,为使其可分解为TL L （L 为下三角阵，主对角线元素>0），a 的取值范围应为。

7．用最小二乘法拟合三点A(0,1),B(1,3),C(2,2)的直线是。

（注意：以上填空题答案标明题号答在答题纸上，答在试卷上的不给予评分。

）二．推导与计算（一）对下表构造f(x)的不超过3次的插值多项式，并建立插值误差公式。

（12分）（二）已知()x x =Φ和()x 'Φ满足∣()x 'Φ-3∣<1。

请利用()x Φ构造一个收敛的简单迭代函数()x ψ，使1(),0,1,......k k x x k +=ψ=收敛。

（8分）（三）利用复化梯形公式计算21x I e dx -=⎰，使其误差限为60.510-⨯，应将区间[0，1]等份。

（8分）（四）设A= 1001005a b b a ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，detA ≠0，推导用a ，b 表示解方程组AX=f 的Seidel(G-S) 迭代法收敛的充分必要条件。

（10分）（五）确定节点及系数，建立如下 GAUSS 型求积公式111220()()dx A f x A f x ≈+⎰。

（10分）（六）对微分方程初值问题'00(,)()y f x y y x y ⎧=⎨=⎩（１）用数值积分法推导如下数值算法：1111(4)3n n n n n hy y f f f +-+-=+++，其中(,)i i i f f x y =，(1,,1)i n n n =-+。

数值分析试题集

..数值分析试题集（试卷一）一（ 10 分）已知 x 1* 1.3409 ，x 2* 1.0125 都是由四舍五入产生的近似值，判断 x 1*x 2* 及 x 1* x 2*有几位有效数字。

二（ 10 分）由下表求插值多项式x 01 2 y2 34 y1－ 1三（ 15 分）设 f ( x)C 4 [a,b] ， H （ x ）是满足下列条件的三次多项式H (a) f (a) , H (b) f (b) , H (c)f (c) , H (c) f (c)( a c b )求 f (x)H ( x) ，并证明之。

12四（ 15 分）计算13 dx ，10 2。

x五（ 15 分）在 [0，2]上取 x 0 0 , x 1 1 , x 22 ，用二种方法构造求积公式，并给出其公式的代数精度。

六（ 10 分）证明改进的尢拉法的精度是 2 阶的。

七（ 10 分）对模型 yy , 0 ，讨论改进的尢拉法的稳定性。

八（ 15分）求方程 x 34x 2 7x 1 0 在 -1.2 附近的近似值，10 3。

-----------------------------------------------------------------------------------------------------------------------------（试卷二）一填空（ 4*2 分）1 {k ( x) } k 0 是区间 [0， 1]上的权函数为( x) x 2 的最高项系数为 1 的正交多项式族，其中10 (x)1，则x0 ( x) dx ------------------- ， 1 ( x) ------------------。

2 12 A，则 A1 4----------- ，( A) ----------------- 。

a 1 2 时， A 可作 LU 分解。

3 设 A，当 a 满足条件 ---------------- 14..4 设非线性方程 f ( x) (x33x23x1)( x 3) 0 ，其根 x1* 3 ， x2*1，则求 x1* 的近似值时，二阶局部收敛的牛顿迭代公式是--------------------------- 。

(完整word版)数值分析(计算方法)期末试卷3及参考答案

[][][]0010012001,,()()n n f x x x x x x -+--参考答案一．填空(每空3分，共30分)1. 截断误差2. )2(--x x ，2)1(-x x ， 10 3. 14.)(2)(21k k k k k k x f x x f x x x '---=+ 5. 6，5，26，9二．计算1. 构造重节点的差商表：所以，要求的Newton 插值为：3()5(1)2(1)(2)(1)(2)(3)N x x x x x x x =--+--+---3243x x =-+插值余项是：2()()(1)(2)3!f R x x x ξ'''=--或：()[,1,2,3,4](1)(2)(3)(4)R x f x x x x x =----2.（1）解：()1f x =时，左10()1f x dx ==⎰，右01A A =+，左＝右得：011A A +=()f x x =时，左101()2f x dx ==⎰，右01B A =+，左＝右得：0112B A += 2()f x x =时，左101()3f x dx ==⎰，右1A =，左＝右得：113A =联立上述三个方程，解得：001211,,363A B A ===3()f x x =时，左101()4f x dx ==⎰，右113A ==，左≠右所以，该求积公式的代数精度是2（2）解：过点0，1构造()f x 的Hermite 插值2()H x ，因为该求积公式代数精度为2，所以有：'212021200010(0)(0)(0)(0)(1()))(0H A H B H f A f B f H x dx A A ++++==⎰其求积余项为：1'1000()[(0)(1)(0)]()f x dx f A f f B f R A -++=⎰112201()()!))((13f H x dx x x dx f x dx η'''--==⎰⎰⎰ 120()(1)3!f x x dx ζ'''=-⎰ ()72f ζ'''=-所以，172k =-3.解：改进的Euler 公式是：1111(,)[(,)(,)]2n n n n n n n n n n y y hf x y hy y f x y f x y ++++=+⎧⎪⎨=++⎪⎩具体到本题中，求解的公式是：11110.2(32) 1.40.60.1[3232](0)1n n n n n n n n n n n n y y x y y x y y x y x y y ++++=++=+⎧⎪=++++⎨⎪=⎩代入求解得：1 1.4y =,1 1.54y =222.276, 2.4832y y ==4．解：设3()25,f x x x =+-则2()32,f x x '=+ 牛顿迭代公式为：1()()k k k k f x x x f x +=-'322532k k k k x x x x +-=-+ 322532k k x x +=+将0 1.5x =代入上式，得1 1.34286x =，2 1.37012x =，3 1.32920x =,4 1.32827x =,5 1.32826x =4540.0000110x x --=<所以，方程的近似根5 1.32826x =5．解，Jacobi 迭代公式是：11231211131521333324k k k k k k k x x x x x x x ++++⎧=--⎪⎪⎪=-⎨⎪⎪=-⎪⎩Gauss-Seidel 迭代公式是：112311211131521333324k k k k k k k x x x x x x x +++++⎧=--⎪⎪⎪=-⎨⎪⎪=-⎪⎩(2) 设其系数矩阵是A ，将A 分解为：A D L U =--，其中300020001D ⎛⎫ ⎪= ⎪ ⎪⎝⎭，000021200,000100000L U --⎛⎫⎛⎫ ⎪ ⎪=-= ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭Jacobi 迭代矩阵是：11030211()0020********J B D L U -⎛⎫--⎛⎫ ⎪ ⎪ ⎪=+=-⎪ ⎪ ⎪- ⎪⎝⎭ ⎪⎝⎭21033100100--⎛⎫⎪ ⎪=- ⎪- ⎪⎝⎭Gauss-Seidel 迭代矩阵是：11300021()220000101000J B D L U ----⎛⎫⎛⎫⎪ ⎪=-= ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭20002112300006206000--⎛⎫⎛⎫ ⎪⎪=- ⎪⎪ ⎪⎪-⎝⎭⎝⎭021********--⎛⎫⎪= ⎪ ⎪⎝⎭二．证明证明：00x >且11()2k k kax x x +=+0k x ⇒> 所以有：111()222k k k k ka a x x x a x x +=+≥=即：数列k x 有下界；2111()()22k k k k k k kx a x x x x x x +=+≤+=所以，迭代序列k x 是单调递减的，由单调递减且有下界的数列极限存在可知序列k x 极限存在。

研究生《数值分析》试卷(带答案)

2009级研究生《数值分析》试卷一.(6分) 已知描述某实际问题的数学模型为xy y x y x u 223),(+=，其中，y x ,由统计方法得到，分别为4,2==y x,统计方法的误差限为0.01,试求出u 的误差限)(u ε和相对误差限)(u r ε.解：)(23)(6)(),()(),()(222y x y x x x y xy y y y x u x x y x u u εεεεε⎪⎭⎫ ⎝⎛++⎪⎭⎫ ⎝⎛-=∂∂+∂∂≈ 6.016.044.001.0)412(01.0)448(=+=⨯++⨯-= 0.010714566.03)()(22=≈+=xy y x u u r εε 二.(6分) 已知函数13)(3+=x x f 计算函数)(x f 的2阶均差]2,1,0[f ,和4阶均差]4,3,2,1,0[f .解：21142512)1()2(]2,1[,311401)0()1(]1,0[=-=--==-=--=f f f f f f9232102]1,0[]2,1[]2,1,0[=-=--=f f f ,0!4)(]4,3,2,1,0[)4(==ξff 三.(6分)试确定求积公式: )]1(')0('[121)]1()0([21)(1f f f f dx x f -++≈⎰的代数精度.解：记⎰=10)(dx x f I )]1(')0('[121)]1()0([21f f f f I n -++= 1)(=x f 时：1110==⎰dx I1]00[121]2[21=-+=n I x x f =)(时：2110==⎰xdx I 21]11[121]1[21=-+=n I2)(x x f =时：31102==⎰dx x I 31]20[121]1[21=-+=n I3)(x x f =时：41103==⎰dx x I 41]30[121]1[21=-+=n I 4)(x x f =时：51104==⎰dx x I 61]40[121]1[21=-+=n I求积公式)]1(')0('[121)]1()0([21)(1f f f f dx x f -++≈⎰具有3次代数精度. 四.(12分) 已知函数122)(23-++=x x x x f 定义在区间[-1,1]上,在空间},,1{)(2x x Span x =Φ上求函数)(x f 的最佳平方逼近多项式.其中,权函数1)(=x ρ,154))(),((,1532))(),((,34))(),((210-==-=x x f x x f x x f ϕϕϕ.解：0))(),(())(),((21))(),((1101101100=====⎰⎰--dx x x x x x dx x x ϕϕϕϕϕϕ32))(),(())(),(())(),((112110220====⎰-dx x x x x x x x ϕϕϕϕϕϕ0))(),(())(),((1131221===⎰-dx x x x x x ϕϕϕϕ 52))(),((11422==⎰-dx x x x ϕϕ解方程组⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛154153234520320320320221a a a 得⎪⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎭⎫ ⎝⎛15161210a a a 则)(x f 的最佳平方逼近多项式为：1516)(2-+=x x x p 五.(16分) 设函数)(x f 满足表中条件:(1) 填写均差计算表((2) 分别求出满足条件22k k k k 的 2次 Lagrange 和 Newton 差值多项式.(3) 求出一个四次插值多项式)(4x H ,使其满足表中所有条件.并用多项式降幂形式表示. 解：12)12)(02()1)(0()20)(10()2)(1()(22+-=----+----=x x x x x x x L12)1)(0(1)0)(1(1)(22+-=--+--+=x x x x x x N 令)2)(1()(12)(24--+++-=x x x b ax x x x H则)2()()2)(1)(()2)(1(22)('4-++--++--+-=x x b ax x x b ax x x ax x x H)1()(-++x x b ax由 ⎩⎨⎧-=+=+⇒⎩⎨⎧=-++-=-=-++-=1220)12(2)2(24)2('2)21)((22)1('44b a b a b a H b a H ,解得 5,3=-=b a 因此1820143)2)(1()53(12)(23424++-+-=--+-++-=x x x x x x x x x x x H 六.(16分)(1). 用Romberg 方法计算⎰31dx x ,将计算结果填入下表(*号处不填).(2). 试确定三点 Gauss-Legender 求积公式⎰∑-=≈110)()(k k k x f A dx x f 的Gauss 点k x 与系数k A ,并用三点 Gauss-Legender 求积公式计算积分: ⎰31dx x .解：过点（1，-1）和点（3，1）作直线得 y t x +=所以积分⎰⎰-+=11312dt t dx x由三次Legendre 多项式 )35(21)(33x x x p -=得得Gauss 点： ,515,0,515210==-=x x x再由代数精度得 ⎪⎪⎪⎩⎪⎪⎪⎨⎧==+==+-==++⎰⎰⎰---32535305155152111220112011210dt x A A dt x A A dt A A A即 ⎪⎩⎪⎨⎧=+=-=++9/10022020210A A A A A A A 解得 ,95,98,95210===A A A所以三点Gauss-Legendre 求积公式为：()⎪⎪⎭⎫⎝⎛++⎪⎪⎭⎫ ⎝⎛-≈⎰-5159509851595)(11f f f dx x f 因此 79746.2515295298515295211=+++-≈+=⎰-dx t I七.(14分)(1) 证明方程02ln =--x x 在区间(1,∞)有一个单根.并大致估计单根的取值范围. (2) 写出Newton 迭代公式,并计算此单根的近似值.(要求精度满足: 5110||-+<-k k x x ). 解：令 2ln )(--=x x x f),1(,011)('∞∈>-=x xx f > 即)(x f 在区间 ),1(∞ 单调增又 04)(,02ln )2(22>-=<-=e e f f 所以 02ln =--x x 在区间 ),1(∞有一单根 ),1(20e x ∈ Newton 迭代公式为1ln 112ln 1-+=----=+k k k k kk k k k x x x x x x x x x令 20=x 计算得八. (12分) 用追赶法求解方程组:⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛022112111131124321x x x x 的解. 解：由计算公式 ⎪⎩⎪⎨⎧-===+====-1,,2,,,2,,111111n i c n i b a c b i i ii i i i i i βααβγγβαα得 ,2,1,1,21,1,24321111======γγγββαα25211322212=⨯-=⇒=+ααβγb 52222222==⇒=αββαc c 53521133323=⨯-=⇒=+ααβγb 35333333==⇒=αββαc c37352144434-=⨯-=⇒=+ααβγb因此 ⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎪⎪⎭⎫⎝⎛135152121137253125121211113112即 LU A = 令 b Ly = 解 ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-022137253125124321y y y y 得⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛23753214321y y y y 令 y Ux =解 ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛237532113515212114321x x x x 得⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛21104321x x x x九. (12分) 设求解初值问题⎩⎨⎧==00)(),('y x y y x f y 的计算格式为:)],(),([111--+++=n n n n n n y x bf y x af h y y ,假设11)(,)(--==n n n n y x y y x y ,试确定参数b a ,的值,使该计算格式的局部截断误差为二阶,即截断部分为: )(3h o .解：)],(),([111--+++=n n n n n n y x bf y x af h y y )](')('[)(1-++=n n n x by x ay h x y])('''21)('')('[)(')(2++-++=n nn n n x y h x hy x y hb x hay x y ++-++=)('''21)('')(')()(32n n n n x by h x by h x y b a h x y对比 ++++=+)('''61)(''21)(')()(321n n n n n x y h x y h x hy x y x y得 ⎩⎨⎧==+2/11b b a ，即 2/1==b a 时该计算格式具有二阶精度.。

数值分析试卷及其答案7

1. 为了使20的近似值的相对误差限小于0.1%，要取几位有效数字？（5分）解、解：设20有n 位有效数字，由Λ4.420=，知41=a令 %1.010811021)20()1()1(1*<⨯=⨯≤----n n r a ε,取 4=n ,%1.010125.0)20(3*<⨯≤-r ε 故 472.420≈2 设方程的迭代法为证明对,均有,其中为方程的根.（5分）证明：迭代函数，对有，3设{}{}100101121,,,span x span x xϕϕ==，分别在12ϕϕ、上求一元素，使其为2[0,1]x C∈的最佳平方逼近，并比较其结果。

（10分）**0111200010012111011220100***010*1***101221221(,)11,(,),211(,),(,),3211(,)1,(,),34111123()611161234a a xdx xdxx dxf x dx f x xdxa aax xaa afϕϕϕϕϕϕϕϕϕϕϕϕδ=+========⋅==⋅=⎧+=⎧⎪=-⎪⎪⇒⇒⇒=-+⎨⎨⎪⎪=+=⎩⎪⎩=-⎰⎰⎰⎰⎰*1(1)设因1*(,)0.00556k kka fϕ=≈∑5分（4分）由结果知（1）比（2）好。

（比较1分）4、用列主元素消元法求解方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--11124112345111321xxx.（10）解：解：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----−−−→−↔⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----111124111123451111212345411121rr**100*1012011110021001010001100011110121021031101000**01**01(2)()11(,)(),(,)(,),201202111 (,)(),(,),(,).203103104 111201202103111202203104x b x b xx dx x x dxx dx f x dx f x dxb bb bϕϕϕϕϕϕϕϕϕϕϕ=+====⋅=======⎧+=⇒+=⎰⎰⎰⎰⎰设**1*10010121122*42220375.24253375.14825()375.24253375.14825.11(,)[375.24253375.14825]0.16406103104k kkbbx x xf b f x dxϕδϕ=⎪⎧≈⎪⇒⎨⎨≈-⎩⎪⎪⎩⇒=-=-=-⨯-⨯≈∑⎰⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-----−−−→−↔⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡------−−−→−-5852510579515130123455795151305852510123455251321312r r r r r r⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡----−−−→−+13513505795151301234513123r r （8分）回代得 3,6,1123==-=x x x 。

华工计算机计算方法(数值分析)考试试卷

考完试了，顺便把记得地题目背下来，应该都齐全了.我印象中也就只有这些题，题目中地数字应该是对地，我也验证过，不过也不一定保证是对地，也有可能我也算错了.还有就是试卷上面地题目可能没有我说地这么短，但是我也不能全把文字背下来，大概意思就是这样吧.每个部分地题目地顺序可能不是这样，但总体就是这四大块.至于每道题目地分值，我记得地就写出来了，有些题目没注意.我题目后面写地结果都是我考试时算出来地，考完了也懒得验证了，可能不一定对，自己把握吧，仅供参考.华南理工大学计算机计算方法（数值分析）考试试卷一填空题（分）1.（分）* ，准确值，求绝对误差(*) ，相对误差(*) ，有效数位是.（分）当插值函数地越大时，会出现龙格现象，为解决这个问题，分段函数不一个不错地办法，请写出分段线性插值、分段三次插值和三次样条插值各自地特点.3.（分）已知和相近，将–变换成可以使其计算结果更准确.4.（分）已知–，求牛顿迭代法地迭代式子.解题思路：. 这里地绝对误差和相对误差是没有加绝对值地，而且要注意是用哪个数减去哪个数得到地值，正负号会不一样；. 可以从它们函数地连续性方面来说明；. 只要满足课本所说地那几个要求就可以；这个记得迭代公式就可以直接写，记不住可以自己推导，就是用泰勒展开式来近似求值得到地迭代公式.我最终地结果是：1.2.分段线性插值保证了插值函数地连续性，但是插值函数地一次导数不一定连续；分段三次既保证了插值函数地连续性，也保证了其一次导数地连续性；三次样条插值保证了插值函数及其一次导数和二次导数地连续性3.()4.– ( –)( )二计算题（分）已知() –，用对分法求其在[ , ]区间内地根，误差要满小于，需要对分多少次？请写出最后地根结果.解题思路：每次求区间地中值并计算其对应地函数值，然后再计算下一个区间中值及函数值，一直到两次区间中值地绝对值小于为止.我最终算得地对分次数是，根地结果为.2.根据以下数据回答相应问题：(1)请根据以上数据构造三次插值函数；(2)请列出差商表并写出三次插值函数.解题思路：() 直接按照书本地定义把公式列出来就可以了，这个要把公式记住了才行，不然也写不了；()差商表就是计算三次插值函数过程中计算到地中间值及结果值，可以先在草稿上按照公式地计算过程把公式写出来，然后把中间用到地值整理成一个表格，这个表格就是差商表了，最后再把公式和表格都写到试卷上就行了.当然也可以先把表格写出来，再用表格地数据写出公式都可以.因为我考试地时候也是先写表格，但是我感觉算地时候容易错，特别是除数地位置，很容易搞错相减地两个地值.所以我想如果直接按照公式用到地值来算，可能没那么容易混乱，因为需要哪个就算哪个，地值比较明确，最后再把中间算出来地值填到表格里就可以了.当然这要看个人喜好了.这里地结果有点长，不好写在这里，自己搞定吧，不难，只是直接套公式就可以了.3. 请用分解法求解以下方程组地解⎪⎩⎪⎨⎧3- = x3 - 9x2 + 6x17 = 3x3+ x2 - 4x11- = x3 - x2 + 2x1解题思路：这个直接套公式算就好了，只要数没有算错，基本都是对地.有时候要注意看是列主元还是直接法，我当时好像没注意，这里应该没有要求用列主元.我最终算得地结果是, , ，其中算出来地矩阵分别是： ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-123121 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--12531124. （分）已知下列矩阵方程，根据以下要求回答问题： ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡210131012⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡321x x x ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-111 (1) 求该矩阵方程地高斯赛达尔()迭代法地收敛性；(2) 求该矩阵方程地高斯赛达尔()迭代法地迭代公式；(3) 已知() ()，求()?解题思路：() 这个证明可以有两种方法，第一种用课本地定义来算，就是将系数矩阵地下三角系数全都乘上一个λ值，然后计算行列式，把所有地λ求出来，只要所有地λ都小于，那么就收敛；第二种方法就是用课本地定理证明，如果系数矩阵是强对角占优地，那么简单迭代法（）和迭代法都收敛，这道题刚好满足条件；() 这个迭代公式只要把矩阵和矩阵求出来就可以写出迭代公式了；() 把()代入()中地迭代公式就可以求出来.我地最终结果是：我直接用强对角占优证明，只写了两句话，不知道老师是不是要求我们用算地...至于强对角占优地判定，书上有，大概意思就是每一行中在主对角线上地那个数地绝对值比旁边所有数地绝对值加起来都要大就是强对角占优了.弱对角就是可以等于.详细定义翻书吧.(2) 我算出来地和矩阵如下：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--02/1003/10，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--03/1002/10，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-2/13/12/1迭代公式就是() () ()(3) () (, , )5. 已知以下方程，请利用最小二乘法求解：⎪⎪⎩⎪⎪⎨⎧0 = 7x2 + 2x1-13= 6x2 + 3x12 = 5x2 + x1-5 = 2x2+ x1解题思路：首先构造一个多变量拟合函数() ，可以把，看成是系数来求解，按照多变量拟合函数求解方法就可以得到结果.我最终算得地结果是：方程组为：⎪⎩⎪⎨⎧⨯=⨯+⨯⨯=⨯+⨯∑∑∑∑∑∑y t t t x t t x yt t t x t t x 22222111212111计算值并代入：⎩⎨⎧=+=+9821141422115x x x x计算地结果为：,请用复化梯形求积公式求出积分dx ⎰10x -e （注：里面地函数是）地近似值，要求误差限满足，请问需要将区间[]分成多少份？解题思路：首先是先把复化梯形求积公式地误差公式写出来，这个要记得，利用误差公式计算出满足精度要求地即可.我最终算得地结果是：误差公式为’’(ŋ)ŋŋ≤≤，≥√≈，也就是满足条件.三证明题（分）已知函数()，其在区间[]内地三个插值点为,（）. 请证明函数()在[]区间内满足下列关系： 6/)]()2/)((4)()[()(b f b a f a f a b dx x f b a +++-≈⎰解题思路：利用这三个插值点写出插值函数，原函数约等于插值函数，所以原函数地积分也约等于插值函数地积分，然后算出插值函数地积分结果就是证明地公式，其实这个就是课本地公式地证明.这个证明过程看课本吧.四程序题（分）前面有一段介绍列主元高斯消元法地步骤地说明（没背下来，都是文字，参考课本吧）请按照列主元高斯消元法地思路将代码中地空格填写完整：1. 输入系数矩阵，右端项及ε；2. 选主元及消元：选主元： ≤≤若 <ε，则打印“求解失败”，停机；否则若≠，则交换地第行和行，交换行和行；消元：––3. 回代若≤ε，则打印“求解失败”，停机，否则(∑+=nijaijxj1)4.打印(…)解题思路：这个直接按照列主元高斯消去法地计算过程去写就好了.结果我写在代码里面了，是按照课本写地，我考试地时候写地应该也是这样.。

岩土工程数值分析试卷试题及参考(附答案)

岩土工程数值分析试题一、简答题（40分）1．简述梁单元、杆单元、连续梁单元、平面三角形常量单元和四边形等参单元的特点（10分）。

答：1）梁单元是由两个节点组成，每一个节点都具有三个方向的线性移动位移和三个方向的旋转位移，因而每个节点具有6个自由度，梁单元具有拉，压，剪，弯，扭的变形刚度。

计算理论成熟，建模方便,计算量小，在工程结构有限元分析中得到广泛的应用，适用于各种截面形式的杆件分析。

2）由有限个构件以一定方式连接起来所形成的结构，在同一平面内的杆系结构，其所受的外力作用线位于该平面内，在杆系中，每一个杆件可视为一个单元，每个单元的端点成为结点。

3）对于每跨各自等截面的连续梁，以每跨为一个单元。

结点编号和单元编号一般是从连续梁的左端顺序编到右端。

由于连续梁各单元的轴线方向一致，各单元坐标系与结构坐标系的方向相同，因此在矩阵位移法的计算过程中无须进行坐标变换，在单元坐标系和结构坐标系中单元刚度矩阵的表达式是相同的。

4) 平面三角形单元具有适应性强的优点，较容易进行网络划分和逼近边界形状，应用比较灵活。

其缺点是它的位移模式是线性函数，单元应力和应变都是常数，精度不够理想。

5) 四边形等参单元能更好地反映物体内的应力变化，适应曲线边界，常使用于弹性力学平平面问题的分析。

八结点单元一共有16个已知的结点位移分量。

2．除有限单元法外，岩土工程常用到哪些数值方法，并对比其优缺点(10分)。

答：岩土工程常用的数值方法包括：有限差分法、边界元法、离散元法、颗粒元法、不连续变形分析法、流形元法、模糊数学方法、概率论与可靠度分析方法、灰色系统理论、人工智能与专家系统、神经网络方法、时间序列分析法。

有限单元法的优缺点：有限单元法的理论基础是虚功原理和基于最小势能的变分原理，它将研究域离散化，对位移场和应力场的连续性进行物理近似。

有限单元法适用性广泛，从理论上讲对任何问题都适用，但计算速度相对较慢。

即，物理概念清晰、灵活、通用、计算速度叫慢。

数值分析试卷及答案

数值分析试卷及答案数值分析试卷一、选择题（共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. 数值解线性方程组的常用迭代法有_________方法和_________方法。

数值分析模拟试卷(二)

迭代矩阵 B 特征方程为
2
………11 分
2 12 1 4 1 2
1 1 0 ， ( B ) 1 4 2
………15 分
所以 SOR 方法迭代格式收敛．四、由拉格朗日插值得二次多项式
1 p 2 ( x) ( x 2 3 x) 2
又设 H ( x) p2 ( x) Ax ( x 1)( x 2) 由于 H (1) f (1) 3 ，解出 A 所以 H ( x) p 2 ( x) 余项 R(x)=

1
x3 1 x(1 x)
0
dx ，问当节点数 n 取何值时，能得到
二：
1、0；
一、填空题（每空 3 分，共 30 分） 2、发散； 3、 f ( x ) 0 ；
*
4、16，90； 6、 3 3 ； 8、b-a；
5、是，因为在 x0 1.5 附近，迭代函数是压缩映射； 7、
0i 2 0 j 2 j i

b a
f ( x)dx

k 0
n
Ak f ( x k ) 的求积系数之和
A
k 0
n
k
__________ .
2 2 1 二、（15 分）已知方阵 A 1 1 1 ， 3 2 1
(1) 证明： A 不能被分解成一个单位下三角阵 L 和一个上三角阵 U 的乘积； (2) 试通过交换 A 的行，使其能实现(Doolittle)分解，并给出其分解； (3) 用上述分解求解方程组 AX=b，其中 b (3.5,2,4) ．
………4 分
5 2
………9 分 ………12 分
5 5 7 x( x 1)( x 2) x 3 7 x 2 x 2 2 2

上海大学2011-2012第二学期数值方法试卷(A含答案)

六、名词解释(共 9 分) (答案仅供参考，允许表述形式不一致) 1. (3 分)迭代法：
答：一般采用迭代法求解方程组，因为迭代法则能保持矩阵的稀疏性，具有计算简单，编制程序容易的优点，并在许多情况下收敛较快。故能有效地解一些高阶方程组
是一种逐次逼近法，从一个假设解开始，通过一系列的迭代求解，最后产生满足精度要求的近似解的方法。如 Jacobi 迭代法，GaussSeidel 迭代法 4. (5 分)写出雅可比迭代法和高斯-赛德尔迭代法的迭代公式，并比较它们的优缺点。（10 2. (3 分)绝对误差分）一个准确值与其在运算中的近似值的差，称为绝对误差。雅可比迭代法：（4 分） n 1 x ( k 1) D 1 (b ( L U ) x ( k ) ) ； 3. (3 分)绝对误差限 xi( k 1) (bi aij x (jk ) ) ，或 aii j 1 绝对误差的绝对值小于等于某个常数 ,该常数称为绝对误差限 j i 高斯-赛德尔迭代法：（4 分）
计算得
命題紙使用說明：1、字迹必須端正，以黑色碳素墨水書寫在框線內，文字與圖均不得剪貼，以保證“掃描”質量； 2、命題紙只作考試（測驗）命題所用，不得移作他用。
第 3 页（共 3 页）
sin(0.34) L2 (0.34) 0.333336
（注至少保留到小数点四位） 3. (5 分)对于线性方程组 Ax b ，已知 A 是高维稀疏矩阵，则一般采用什么方法求解？为什么？
n n 1 (bi aij x (jk 1) aij x (jk ) ) aii j i j i ，或
xi( k 1)
七、简答题(共 23 分)： 1. (8 分)试写出数值积分中的梯形公式、辛普森公式、辛普森 3/8 公式和布尔公式，且给出它们各自的精度值。设 xk=x0+kh 为等距节点，且 fk=f(xk), 则四个数值积分公式分别为： x1 h 梯形公式精度为 1，具体公式为: f ( x)dx ( f 0 f1 ) x0 2 x2 h 辛普森公式精度为 3，具体公式为： f ( x)dx ( f 0 4 f1 f 2 ) x0 3 辛普森 3/8 公式精度为 3，具体公式为： f ( x)dx

数值分析(计算方法)期末试卷及参考答案

江西理工大学大学二计算-- -------- ---- ---- ----号 --学线 ------2013 至 2014学年第一学期试卷试卷课程数值剖析年级、专业︵B题号一二三四五六七八九十总分︶得分第1. 给定数据表：（ 15 分）x i 1 2 f ( x i )2 3f ' (x i )------ ----------- 名- --姓封 -- -------- -------- -----------密------------- --- 级---班- -- 、 - -- 业- --专--一填空 (每空 3 分，共 30 分)1. 在一些数值计算中，对数据只好取有限位表示，如时所产生的偏差称为。

2. 设 f ( x)x 7 x 6 1 ， f [30 ,31 ]f [3 0 ,31, ,37]， f [3 0 ,31, ,38 ]3. 5 个节点的牛顿 -柯特斯公式代数精度是。

4. 求方程 x2cos x 根的 Newton 迭代格式为5. 设(1, 3,0,2) ，则1，2 12；设 A5 ，则 A41页 2 1.414 ，这︵共3页，︶。

江。

西理，工大学。

大学教务处(1) 结构 Hermit 插值多项式 H 2 ( x) ，并计算 f (1.5) 。

(2) 写出其插值余项，并证明之。

- ---------------------- 号- ---- 学--------线-------------------------名- 2. 已知方程x2 ln x 4 0 ，取 x0 1.5 ，用牛顿迭代法求解该方程的根，要求 x k 1 x k 1试10 3时停止迭代。

（10分）卷︵B︶第2页︵共4．用Euler方法求解初值问题y'x yy(0) 0---封姓- ------------------------密------------- 级---- 班--- 、- -- 业-- 专-1 33. 确立求积公式 f ( x)dx Af (0) Bf ( x1 ) Cf (1)0页︶中的待定参数A, B,C , x1，使其代数精度尽可能高，并指出其代数精度。

数值分析试卷及答案

二1 求A的LU分解，并利用分解结果求解由紧凑格式故从而故2求证：非奇异矩阵不一定有LU分解证明设非奇异，要说明A不一定能做LU分解，只需举出一个反例即可。

现考虑矩阵，显然A为非奇异矩阵。

若A有LU分解，则故，而，显然不能同时成立。

这矛盾说明A不能做LU分解，故只假定A非奇异并不能保证A能做LU分解，只有在A的前阶顺序主子式时才能保证A一定有LU分解。

3用追赶法求解如下的三对角方程组解设有分解由公式其中分别是系数矩阵的主对角线元素及其下边和上边的次对角线元素，故有从而有故，，，故，，，4设A是任一阶对称正定矩阵，证明是一种向量范数证明（1）因A正定对称，故当时，，而当时，（2）对任何实数，有（3）因A正定，故有分解，则故对任意向量和，总有综上可知，是一种向量范数。

5 设，，已知方程组的精确解为（1）计算条件数；（2）若近似解，计算剩余；（3）利用事后误差估计式计算不等式右端，并与不等式左边比较，此结果说明了什么？解（1）（2）（3）由事后误差估计式，右端为而左端这表明当A为病态矩阵时，尽管剩余很小，误差估计仍然较大。

因此，当A病态时，用大小作为检验解的准确度是不可靠的。

6矩阵第一行乘以一数成为，证明当时，有最小值证明设，则又故从而当时，即时，有最小值，且7 讨论用雅可比法和高斯-赛德尔法解方程组时的收敛性。

如果收敛，比较哪一种方法收敛较快，其中解对雅可比方法，迭代矩阵，故雅可比法收敛。

对高斯-赛德尔法，迭代矩阵，故高斯-赛德尔法收敛。

因=故高斯-赛德尔法较雅可比法收敛快。

8设，求解方程组，求雅可比迭代法与高斯-赛德尔迭代法收敛的充要条件。

解雅可比法的迭代矩阵，故雅可比法收敛的充要条件是。

高斯-赛德尔法的迭代矩阵，故高斯-赛德尔法收敛的充要条件是。

9 设求解方程组的雅可比迭代格式为，其中，求证：若，则相应的高斯-赛德尔法收敛。

证明由于是雅可比法的迭代矩阵，故又，故，即，故故系数矩阵A按行严格对角占优，从而高斯-赛德尔法收敛。

数值分析试卷及答案

数值分析试卷及答案**注意：以下是一份数值分析试卷及答案，试卷和答案分别按照题目和解答的格式排版，以确保整洁美观，语句通顺。

**---数值分析试卷一、选择题（每题2分，共20分）1. 数值分析是研究如何用计算机处理数值计算问题的一门学科。

以下哪个选项不是数值分析的应用领域？A. 金融风险评估B. 天气预测C. 数据挖掘D. 图像处理2. 在数值计算中，稳定性是指算法对于输入数据的微小扰动具有较好的性质。

以下哪个算法是稳定的？A. 高斯消元法B. 牛顿迭代法C. 不动点迭代法D. 雅可比迭代法二、填空题（每题3分，共30分）1. 下面关于插值多项式的说法中，不正确的是：一般情况下，插值多项式的次数等于插值点的个数减1。

2. 线性方程组中，如果系数矩阵A是奇异的，则该方程组可能无解或有无穷多解。

......三、解答题（共50分）1. 请给出用割线法求解非线性方程 f(x) = 0 的迭代格式，并选择合适的初始值进行计算。

解：割线法的迭代公式为：x_(k+1) = x_k - f(x_k) * (x_k - x_(k-1)) / (f(x_k) - f(x_(k-1)))选择初始值 x0 = 1，x1 = 2 进行计算：迭代1次得到：x2 = x1 - f(x1) * (x1 - x0) / (f(x1) - f(x0))迭代2次得到：x3 = x2 - f(x2) * (x2 - x1) / (f(x2) - f(x1))继续迭代直至满足精度要求。

2. 对于一个给定的线性方程组，高斯消元法可以用来求解其解空间中的向量。

请简要描述高斯消元法的基本思想并给出求解步骤。

高斯消元法的基本思想是通过一系列的行变换将线性方程组化为上三角形式，然后再通过回代求解方程组的未知数。

求解步骤如下：步骤1：将方程组表示为增广矩阵形式，即将系数矩阵和常数向量连接在一起。

步骤2：从第一行开始，选取第一个非零元素作为主元，然后通过行变换将其它行的该列元素消去。

数值分析试卷习题5

32 第五章习题解答与问题一、习题解答1．求经过A (0，1)，B (1，2)，C (2，3)三个样点的插值多项式解：令x 0 = 0，x 1 = 1，x 2 = 2，则有f (x 0)=1，f (x 1)=2，f (x 2)=3，由Lagrange 二次插值公式)())(())(()())(())(()())(())(()(2120210121012002010212x f x x x x x x x x x f x x x x x x x x x f x x x x x x x x x L ----+----+----=3)12)(02()1)(0(2)21)(01()2)(0(1)20)(10()2)(1(⨯----+⨯----+⨯----=x x x x x x = x+12．已知函数)(x f y =的数据如下表试作一个三次插值多项式P 3(x ),利用P 3(x )计算3解：令x k = kk根据Newton )]}()[({))(()()(2342121213412213-+-++=--+-++=x x xx x x x x x x P由于被插值函数xx f 3=)(，故取 x = 1/2，便得222134212122112133=-+-++=≈)]}()[({)/(P3．已知函数y = f (x )解：由于x=0是二重零点，令3。

又由3，H 3(1)=1得方程组33⎩⎨⎧=+-=-11b a a b 解之：a =1，b = 0 所以，H 3(x ) = x 3。

4．设被插值函数f (x )在区间[x 0，x 1]上具有2阶连续导数，求证：两点线性插值函数L (x )的误差界满足不等式8)(|)(|max |)(|20110x x x f x R x x x -''≤≤≤证：由拉格朗日插值误差定理，得))((!2)()()()(10x x x x f x L x f x R --''=-=ξ 令h (x ) = (x – x 0)(x – x 1)，求导数并令其为零，可得极值点x *=0.5×(x 0 + x 1)。

数值分析试卷

数值分析考试题（一）满分70分一、选择题：（共3道小题，第1小题4分，第2、3小题3分，共10分） 1、将A 分解为U L D A --=，其中),,(2211nna a a diag D =，若对角阵D非奇异（即),1,0n i a ii =≠，则b Ax =化为b D x U L D x 11)(--++=（1）若记b D f U L D B 1111),(--=+= （2）则方程组（1）的迭代形式可写作 )2,1,0(1)(1)1( =+=+k f x B xk k （3）则（2）、（3）称【】(A)、雅可比迭代。

(B)、高斯—塞德尔迭代 (C)、LU 分解 (D)、Cholesky 分解。

2、记*x x e k k -=，若0lim1≠=+∞→c ee pk k k （其中p 为一正数）称序列}{k x 是【】(A)、p 阶收敛； (B)、1阶收敛； (C)、矩阵的算子范数； (D)、p 阶条件数。

3、牛顿切线法的迭代公式为【】(A)、 )()(1kx f x f x x k k k '-=+ (B)、)()())((111--+---=k k k k k k k x f x f x x x f x x1)()()1()()()(x xfxf xf k i k i k i ∂∂+=+ (D)、 )()()()1(k k k x f x x-=+二、填空题：（共2道小题，每个空格2分，共10分）1、设0)0(f =，16)1(f =，46)2(f =，则一阶差商=]1,0[f ，二阶差商=]1,2,0[f ，)x (f 的二次牛顿插值多项式为2、用二分法求方程01x x )x (f 3=-+=在区间]1,0[内的根，进行第一步后根所在的区间为，进行第二步后根所在的区间为。

三、计算题：（共7道小题，第1小题8分，其余每小题7分，共50分）1、表中各*x 都是对准确值x 进行四舍五入得到的近似值。

昆明理工大学数值分析各年考试题及答案

昆明理工大学数值分析考试题（07）一．填空（每空3分，共30分）1．设A 0.231x =是真值0.229T x =的近似值，则A x有位有效数字。

2．若74()631f x x x x =+++，则017[2,2,...2]f =，018[2,2,...2]f = 。

3． A=1031⎡⎤⎢⎥-⎣⎦,则1A = ；A ∞= ；2A = 2()cond A = 。

4．求方程()x f x =根的牛顿迭代格式是。

5．设105%x =±，则求函数()f x =的相对误差限为。

6．A=2101202a a ⎛⎫⎪ ⎪ ⎪⎝⎭,为使其可分解为TL L （L 为下三角阵，主对角线元素>0），a 的取值范围应为。

7．用最小二乘法拟合三点A(0,1),B(1,3),C(2,2)的直线是。

（注意：以上填空题答案标明题号答在答题纸上，答在试卷上的不给予评分。

）二．推导与计算（一）对下表构造f(x)的不超过3次的插值多项式，并建立插值误差公式。

（12分）（二）已知()x x =Φ与()x 'Φ满足∣()x 'Φ-3∣<1。

请利用()x Φ构造一个收敛的简单迭代函数()x ψ，使1(),0,1,......k k x x k +=ψ=收敛。

（8分）（三）利用复化梯形公式计算21x I e dx -=⎰，使其误差限为60.510-⨯，应将区间[0，1] 等份。

（8分）（四）设A= 1001005a b b a ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，detA ≠0，推导用a ，b 表示解方程组AX=f 的Seidel(G-S) 迭代法收敛的充分必要条件。

（10分）（五）确定节点与系数，建立如下 GAUSS 型求积公式111220()()A f x A f x ≈+⎰。

数值分析试卷及答案

模拟试卷（一）一、填空题（每小题3分，共30分）1．有3个不同节点的高斯求积公式的代数精度是次的.2．设，，则= .，= ______.3．已知y=f(x)的均差（差商），，，, 那么均差= .4．已知n=4时Newton－Cotes求积公式的系数分别是：则＝ .5．解初始值问题的改进的Euler方法是阶方法；6．求解线性代数方程组的高斯—塞德尔迭代公式为，若取, 则 .7．求方程根的牛顿迭代格式是 .8．是以整数点为节点的Lagrange插值基函数，则= .9．解方程组的简单迭代格式收敛的充要条件是 .10．设，则的三次牛顿插值多项式为，其误差估计式为.二、综合题（每题10分，共60分）1．求一次数不超过4次的多项式满足：，，，.2．构造代数精度最高的形式为的求积公式，并求出其代数精度.3．用Newton法求方程在区间内的根, 要求.4．用最小二乘法求形如的经验公式拟合以下数据：5．用矩阵的直接三角分解法解方程组.6 试用数值积分法建立求解初值问题的如下数值求解公式，其中.三、证明题（10分）设对任意的，函数的导数都存在且,对于满足的任意，迭代格式均收敛于的根.参考答案一、填空题1．5； 2. 8, 9 ; 3. ； 4. ； 5. 二；6. , (0.02，0.22，0.1543)7. ; 8. ; 9. ;10.二、综合题1．差商表：1 1 12 2151515575720204272152230781其他方法：设令，，求出a和b.2．取，令公式准确成立，得：, , ,.时，公式左右；时，公式左, 公式右∴公式的代数精度.3．此方程在区间内只有一个根，而且在区间（2，4）内。

设则，，Newton法迭代公式为，取，得。

4．，，.解方程组，其中，解得：所以， .5．解设由矩阵乘法可求出和解下三角方程组有，，，.再解上三角方程组得原方程组的解为，，，.6 解初值问题等价于如下形式，取，有，利用辛卜森求积公式可得.三、证明题证明将写成，由于，所以所以迭代格式均收敛于的根.模拟试卷（二）一、填空题（每小题3分，共30分）1．分别用2.718281和2.718282作数的近似值，则其有效位数分别有位和位；2．设，，则= ________，= .3．对于方程组, Jacobi迭代法的迭代矩阵是=________.4．设，则差商=__________，=_______.5．已知, 则条件数_________.6．为使两点的数值求积公式具有最高的代数精确度，则其求积基点应为=__________， =__________7．解初始值问题近似解的梯形公式是8．求方程根的弦截法迭代公式是9. 计算积分，取4位有效数字，用梯形公式计算求得的近似值是，用辛卜生公式计算的结果是10．任一非奇异矩阵的条件数＝，其一定大于等于二、综合题（每题10分，共60分）1 证明方程在区间有且只有一个根，若利用二分法求其误差不超过近似解，问要迭代多少次？2 已知常微分方程的初值问题：试用改进的Euler方法计算的近似值，取步长.3 用矩阵的分解法解方程组 .4 用最小二乘法求一个形如的经验公式，使它与下列数据拟合.x 1.0 1.4 1.8 2.2 2.6y0.9310.4730.2970.2240.1685 设方程组，试考察解此方程组的雅可比迭代法及高斯－赛德尔迭代法的收敛性。