10-11数值分析期末试卷 (3)
数值分析试题集
..数值分析试题集（试卷一）一（ 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 极限存在。
数值分析期末试题及答案
数值分析期末试题及答案试题一：1. 简答题（共10分）a) 什么是数值分析？它的主要应用领域是什么？b) 请简要解释迭代法和直接法在数值计算中的区别。
2. 填空题（共10分）a) 欧拉方法是一种______型的数值解法。
b) 二分法是一种______法则。
c) 梯形法则是一种______型的数值积分方法。
3. 计算题（共80分）将以下函数进行数值求解：a) 通过使用二分法求解方程 f(x) = x^3 - 4x - 9 = 0 的近似解。
b) 利用欧拉方法求解微分方程 dy/dx = x^2 + 2x + 1, y(0) = 1 在 x = 1 处的解。
c) 使用梯形法则计算积分∫[0, π/4] sin(x) dx 的近似值。
试题二：1. 简答题（共10分）a) 请解释什么是舍入误差，并描述它在数值计算中的影响。
b) 请解释牛顿插值多项式的概念及其应用。
2. 填空题（共10分）a) 数值稳定性通过______号检查。
b) 龙格－库塔法是一种______计算方法。
c) 零点的迭代法在本质上是将方程______转化为______方程。
3. 计算题（共80分）使用牛顿插值多项式进行以下计算：a) 已知插值节点 (-2, 1), (-1, 1), (0, 2), (1, 4)，求在 x = 0.5 处的插值多项式值。
b) 已知插值节点 (0, 1), (1, 2), (3, 7)，求插值多项式，并计算在 x = 2 处的值。
c) 使用 4 阶龙格－库塔法求解微分方程 dy/dx = x^2 + 1, y(0) = 1。
答案：试题一：1. a) 数值分析是研究使用数值方法解决数学问题的一门学科。
它的主要应用领域包括数值微积分、数值代数、插值和逼近、求解非线性方程、数值积分和数值解微分方程等。
b) 迭代法和直接法是数值计算中常用的两种方法。
迭代法通过反复迭代逼近解，直到满足所需精度为止；而直接法则通过一系列代数运算直接得到解。
数值分析期末试题及答案
数值分析期末试题及答案一、选择题（每题5分，共20分）1. 在数值分析中，下列哪个算法不是用于求解线性方程组的？A. 高斯消元法B. 牛顿法C. 雅可比法D. 追赶法答案：B2. 插值法中，拉格朗日插值法属于：A. 多项式插值B. 样条插值C. 线性插值D. 非线性插值答案：A3. 以下哪个选项不是数值分析中的误差来源？A. 截断误差B. 舍入误差C. 计算误差D. 测量误差答案：C4. 在数值积分中，梯形法则的误差项是：A. O(h^2)B. O(h^3)C. O(h)D. O(1)答案：A二、填空题（每题5分，共20分）1. 牛顿插值法中，插值多项式的一般形式为：______。
答案：f(x) = a_0 + a_1(x-x_0) + a_2(x-x_0)(x-x_1) + ...2. 牛顿迭代法求解方程的根时，迭代公式为：x_{n+1} = x_n -f(x_n) / __________。
答案：f'(x_n)3. 在数值分析中，______ 用于衡量函数在区间上的近似积分值与真实积分值之间的差异。
答案：误差4. 线性方程组的解法中，______ 法是利用矩阵的LU分解来求解。
答案：克兰特三、解答题（每题10分，共60分）1. 给定函数f(x) = e^(-x)，使用拉格朗日插值法，求x = 0.5时的插值值。
解答：首先选取插值节点x_0 = 0, x_1 = 0.5, x_2 = 1，对应的函数值分别为f(0) = 1, f(0.5) = e^(-0.5), f(1) = e^(-1)。
拉格朗日插值多项式为：L(x) = f(0) * (x-0.5)(x-1) / (0-0.5)(0-1) + f(0.5) * (x-0)(x-1) / (0.5-0)(0.5-1) + f(1) * (x-0)(x-0.5) / (1-0)(1-0.5)将x = 0.5代入得：L(0.5) = 1 * (0.5-0.5)(0.5-1) / (0-0.5)(0-1) + e^(-0.5) * (0.5-0)(0.5-1) / (0.5-0)(0.5-1) + e^(-1) * (0.5-0)(0.5-0.5) / (1-0)(1-0.5)计算得L(0.5) = e^(-0.5)。
复旦大学《数值分析》2016-2017学年第二学期期末试题(B卷)
课程编号：A072121复旦大学2016-2017学年第二学期2016级《数值分析》期末试卷(B 卷)（本试卷共6页、八个大题，满分100分；答题前请检查是否有漏印、缺页和印刷不清楚的情况，如有此种情况，请及时向监考教师反映）一、求解下列各题（每小题6分）1.已知直线3221:+==-z m y x L 与平面02:=++-D z y x π平行，且L 到π的距离为6,求m 与D 的值.2.设)()(1y x y xy f xz ++=ϕ,其中ϕ,f 二阶可导,求y x z ∂∂∂2.3.计算第二类曲线积分dy y xdx y x I L ⎰+= 2 ，其中L 是曲线x y =上从点)1 , 1(A 到点)2 , 4(B 的弧段.4.设有级数)11ln(1)1(11n n n pn +-∑∞=-,指出p 在什么范围内取值时级数绝对收敛,在什么范围内取值时级数条件收敛,在什么范围内取值时级数发散(要说明理由).二、解下列各题（每小题7分）1.已知n是曲面1222=+-z y x 在点)1 , 2 , 2(处指向z 增大方向的单位法向量,z z xy u ln 2-=,求)1 ,2 , 2(nu∂∂.2.将函数231)(2++=x x x f 展开成)1(-x 的幂级数,并求收敛区间及)1()5(f 的值.3.计算三重积分⎰⎰⎰Ω=zdV x I 2,其中Ω是由柱面2x y =与平面1=y ,0=z ，2=z 所围成的立体.4．求二元函数y x y x x y x f z 293),(223+---==的极值点与极值.三、(8分)设1)(2+=x x f ，ππ≤≤-x ,将)(x f 展开成以π2为周期的傅里叶级数.四、（8分）设V 是由曲面222y x z --=与22y x z +=围成的立体,求V 的表面积.五、(8分)计算第二类曲面积分⎰⎰++=Sdxdy dzdx y dydz x I 33,其中S 是曲面22y x z +=)10(≤≤z 的下侧.六、（8分）求幂级数∑∞=+12)(n n x n n 的收敛域与和函数.七、(8分)已知在半平面0>x 内dy y x y x dx y x y x λλ))(())((2222++++-为二元函数),(y x f 的全微分.(1)求λ的值；(2)求)0 , 2()3 , 1(f f -的值.八、（8分）设}|),,{()(2222t z y x z y x t ≤++=Ω，其中0>t .已知)(x f 在) , 0[∞+内连续，又设⎰⎰⎰Ω++=)( 222)()(t dxdydz z y x f t F .（1）求证：)(t F 在) , 0(∞+内可导，并求)(t F '的表达式；（2）设0)0(≠f ，求证：级数∑∞=-'111(n n F n λ在0>λ时收敛，0≤λ时发散.(此页纸不够时可写到背面)2016级《数值分析》期末试卷(B 卷)参考答案与评分标准一．求解下列各题1.直线过点(1,0,-2)，方向向量}3,,2{m s = ，平面法向量}2,1,1{-=n--------------------2分062}3,,2{}2,1,1{=+-=⋅-⇒⊥m m s n8=⇒m ------------------------------------4分6411|401|=+++--=D d 3,9-=⇒D ------------------------------------------------------6分或过)2,0,1(-与π,L 垂直的直线方程为22111+=-=-z y x 与π交点：321,615,69-=-=-=D z D y D x 6)2321()65()169(222=+-+-+--=D D D d .9,3-=⇒D 2.)()()(12y x y xy f x yxy f x x z +'+'+-=∂∂ϕ------------------------------------------------3分)()()(2y x y y x xy f y yx z+''++'+''=∂∂∂ϕϕ------------------------------------------------6分3.dy y x dx y x L+⎰2=dxxx x x x ⎰⋅+41221(-------------------------------------------------3分10139]2152[4125=+=x x -------------------------------------6分4.11)11ln(1lim 1=++∞→p pn n n n∑∑∞=∞=++-11)11ln(1|)11ln(1)1(|n pn p nn n n n 与有相敛散性----------------------------------2分1)0>p ，绝对收敛；2)01≤<-p ，条件收敛；3)1-<p ，发散。
数值分析试卷及答案
数值分析试卷及答案数值分析试卷一、选择题（共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. 数值解线性方程组的常用迭代法有_________方法和_________方法。
数值分析试题(卷)和答案解析
【试题__2009___年~__2010___年第 一学期课程名称： 数值分析 专业年级： 2009级（研究生） 考生学号： 考生姓名： 试卷类型： A 卷 √ B 卷 □ 考试方式: 开卷 √ 闭卷 □………………………………………………………………………………………………………一. 填空题（本大题共4小题，每小题4分，共16分）1.设有节点012,,x x x ，其对应的函数()y f x =的值分别为012,,y y y ，则二次拉格朗日插值基函数0()l x 为 。
-2.设()2f x x =，则()f x 关于节点0120,1,3x x x ===的二阶向前差分为 。
3.设110111011A -⎡⎤⎢⎥=--⎢⎥⎢⎥-⎣⎦，233x ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，则1A ＝ ，1x = 。
4. 1n +个节点的高斯求积公式的代数精确度为 。
二．简答题（本大题共3小题，每小题8分，共24分）1. 哪种线性方程组可用平方根法求解为什么说平方根法计算稳定2. 什么是不动点迭代法()x ϕ满足什么条件才能保证不动点存在和不动点迭代序列收敛于()x ϕ的不动点3. 设n 阶矩阵A 具有n 个特征值且满足123n λλλλ>≥≥≥，请简单说明求解矩阵A 的主特征值和特征向量的算法及流程。
三．求一个次数不高于3的多项式()3P x ，满足下列插值条件：。
i x 1 2 3i y2 4 12 <3i y '并估计误差。
（10分）四．试用1,2,4n =的牛顿-科特斯求积公式计算定积分1011I dx x=+⎰。
（10分） 五．用Newton 法求()cos 0f x x x =-=的近似解。
（10分） 六．试用Doolittle 分解法求解方程组： ，12325610413191963630x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦⎣⎦（10分）七．请写出雅可比迭代法求解线性方程组123123123202324812231530x x x x x x x x x ++=⎧⎪++=⎨⎪-+=⎩ 的迭代格式，并判断其是否收敛（10分）八．就初值问题0(0)y yy y λ'=⎧⎨=⎩考察欧拉显式格式的收敛性。
《数值分析》期末复习题(1)
《数值分析》期末复习题一、单项选择题1. 数值x *的近似值x =0.32502×10-1，若x 有5位有效数字，则≤-*x x ( ).(A)21×10-3 (B) 21×10-4 (C) 21×10-5 (D) 21×10-6 2. 设矩阵A ＝10212104135⎡⎤⎢⎥⎢⎥⎣⎦，那么以A 为系数矩阵的线性方程组A X ＝b 的雅可比迭代矩阵为( )(A)00.20.10.200.40.20.60--⎡⎤⎢⎥--⎢⎥--⎣⎦(B)10.20.10.210.40.20.61⎡⎤⎢⎥⎢⎥⎣⎦(C) 00.20.10.200.40.20.60⎡⎤⎢⎥⎢⎥⎣⎦ (D)021204130⎡⎤⎢⎥⎢⎥⎣⎦3. 已知(1)1,(2)4,(3)9f f f ===，用拉格朗日2次插值，则(2.5)f =( )(A) 6.15 (B) 6.25 (C) 6.20 (D) 6.10 4. 抛物形求积公式的代数精度是( )A. 1，B. 2 ，C. 3，D. 45. 改进欧拉格式的 局部截断误差是( ). (),A O h 2. (),B O h 3. (),C O h 4. ().D O h二、填空题1、以722作为π的近似值，它有（ ）位有效数字； 2、经过)1,2( ),2,1( ),1,0(C B A 三个节点的插值多项式为（ ）； 3、用高斯-赛德尔迭代法解方程组⎩⎨⎧-=+-=+,10,232121x bx bx x 其中b 为实数，则方法收敛的充分条件是b 满足条件（ ）；4、取步长为1.0=h ，用欧拉法计算初值问题22',(0)0,y x y y ⎧=+⎨=⎩的解函数)(x y ，它在3.0=x 的近似值为（ ）；5、已知方程0sin 1=--x x 在)1,0(有一个根，使用二分法求误差不大于41021-⨯的近似解至少需要经过（ ）次迭代。
(完整)数值分析历年考题
数值分析A 试题2007.1第一部分:填空题10⨯51.设3112A ⎛⎫= ⎪⎝⎭,则A ∞=___________ 2()cond A =___________2.将4111A ⎛⎫= ⎪⎝⎭分解成TA LL =，则对角元为正的下三角阵L =___________，请用线性最小二乘拟合方法确定拟合函数()bx f x ae =中的参数:a = ___________ b =___________4.方程13cos 2044x x π--=在[0,1]上有 个根，若初值取00.95x =,迭代方法113cos 244k k x x π+=-的收敛阶是5.解方程2210x x -+=的Newton 迭代方法为___________,其收敛阶为___________6。
设()s x = 3232323,[0,1]31,[1,2]ax x x x x x bx x +-+∈--+∈为三次样条函数，则a = ___________ b =___________ 7。
要想求积公式：1121()(()f x dx A f f x -≈+⎰的代数精度尽可能高，参数1A = ___________ 2x =___________此时其代数精度为：___________8.用线性多步法2121(0.50.5)n n n n n y y h f f f ++++-=-+来求解初值问题00'(,),(),y f x y y x y ==其中(,)n n n f f x y =，该方法的局部截断误差为___________，设,0,f y μμ=〈其绝对稳定性空间是___________9。
用线性多步法2121()n n n n n y ay by h f f ++++-+=-来求解初值问题00'(,),(),y f x y y x y ==其中(,)n n n f f x y =,希望该方法的阶尽可能高，那么a = ___________ b =___________,此时该方法是几阶的：___________10。
数值分析考试卷及详细答案解答
数值分析考试卷及详细答案解答姓名班级学号⼀、选择题1.()2534F,,,-表⽰多少个机器数(C ).A 64B 129C 257D 2562. 以下误差公式不正确的是( D)A ．()()()1212x *x *x *x *εεε-≈+B ．()()()1212x *x *x *x *εεε+≈+C ．()()()122112x *x *x *x *x x *εεε?≈+ D ．()()()1212x */x *x *x *εεε≈-3. 设)61a =, 从算法设计原则上定性判断如下在数学上等价的表达式，哪⼀个在数值计算上将给出a 较好的近似值？(D )A6)12(1+ B 27099- C 3)223(- D3)223(1+4. ⼀个30阶线性⽅程组, 若⽤Crammer 法则来求解, 则有多少次乘法? ( A )A 31×29×30!B 30×30×30!C 31×30×31!D 31×29×29!5. ⽤⼀把有毫⽶的刻度的⽶尺来测量桌⼦的长度, 读出的长度1235mm, 桌⼦的精确长度记为（ D ）A 1235mmB 1235-0.5mmC 1235+0.5mmD 1235±0.5mm⼆、填空1.构造数值算法的基本思想是近似替代、离散化、递推化。
2.⼗进制123.3转换成⼆进制为1111011.01001。
3.⼆进制110010.1001转换成⼗进制为 50.5625 。
4. ⼆进制0101.转换成⼗进制为57。
5.已知近似数x*有两位有效数字，则其相对误差限 5% 。
6. ln2=0.69314718…，精确到310-的近似值是 0.693 。
7.31415926x .π==，则131416*x .=，23141*x .=的有效数位分别为5 和3 。
8.设200108030x*.,y*.==-是由精确值x y 和经四舍五⼊得到的近似值，则x*y*+的误差限 0.55×10-3 。
数值分析期末考试题及答案
数值分析期末考试题及答案一、选择题（每题2分，共20分）1. 在数值分析中，下列哪个算法用于求解线性方程组？A. 牛顿法B. 高斯消元法C. 插值法D. 傅里叶变换答案：B2. 以下哪个选项不是数值分析中的误差类型？A. 舍入误差B. 截断误差C. 测量误差D. 累积误差答案：C3. 多项式插值中，拉格朗日插值法的特点是：A. 插值点必须等距分布B. 插值多项式的次数与插值点的个数相同C. 插值多项式是唯一的D. 插值多项式在插值点处的值都为1答案：B4. 在数值分析中，下列哪个方法用于求解非线性方程？A. 辛普森法则B. 牛顿迭代法C. 欧拉法D. 龙格-库塔法答案：B5. 以下哪个是数值稳定性的指标？A. 收敛性B. 收敛速度C. 条件数D. 误差传播答案：C二、简答题（每题10分，共20分）1. 简述高斯消元法求解线性方程组的基本原理。
答案：高斯消元法是一种直接解法，通过行变换将增广矩阵转换为上三角形式，然后通过回代求解线性方程组。
它包括三个基本操作：行交换、行乘以非零常数、行相加。
2. 解释什么是数值稳定性，并举例说明。
答案：数值稳定性是指数值解对输入数据小的扰动不敏感的性质。
例如，某些数值方法在计算过程中可能会放大舍入误差，导致结果不可靠，这样的方法就被认为是数值不稳定的。
三、计算题（每题15分，共30分）1. 给定线性方程组：\[\begin{align*}x + 2y - z &= 4 \\3x - y + 2z &= 1 \\-x + y + z &= 2\end{align*}\]使用高斯消元法求解该方程组，并给出解。
答案：首先将增广矩阵转换为上三角形式，然后回代求解，得到\( x = 1, y = 2, z = 1 \)。
2. 给定函数 \( f(x) = x^2 - 3x + 2 \)，使用拉格朗日插值法在\( x = 0, 1, 2 \) 处插值，并求出插值多项式。
数值分析期末考试题之经典例题
题型一:有效数字1,确定113的首位数字x 1,要使113的近似值x *的相对误差不超过0.5×10-5，至少要保留几位有效数字.(2010-2011)1*1151211||10100.5102226n n r x n e x n ---=≤⨯=⨯≤⨯⨯≥=解答：设至少要保留位有效数字，则有解得， n 5.7取位有效数字.2，要使112的相对误差不超过0.5×10-4，至少要保留几位有效数字?(2009-2010) 3,已知21.787654为有效数，确定其绝对误差界与相对误差界.(2007-2008)*6*118711||102111||1010102224n r e e x ----=⨯=⨯=⨯=⨯⨯解答：4,已知30.49876为有效数，确定其绝对误差界.(2006-2007B)5,设有效数x=12.4567,确定x 的绝对误差界.(2004-2005)题型二：插值多项式1,已知f(x)的函数值：f(0)=-2, f(1)=1, f(2)=5, 用反插值法求f(x)=0在[0，2]内的近似根x *.(2010-2011)11111202012012010210122021()()()()()()()()()()()()()()()()()(2)(5)(2)(1)012(12)(15)(52)(51)2991422884y y y y y y y y y y y y y L y f y f y f y y y y y y y y y y y y y y y y y y ----------=⋅+⋅+⋅------+-+-=+⨯+⨯+-+-=+-解答：对y=f(x)的反函数x=f 进行二次插值2*229(0)42y x L ≈=故，2,已知f(x)的如下函数值及导数值：f(-1)=1, f(0)=2, f ’(0)=3, f(1)=7; (1),建立不超过3次的埃尔米特插值多项式H 3(x);(2),x ∈[-1,1], 确定用H 3(x)代替f(x)的误差界（已知|f (4)(x)|≤M 4，x ∈[-1,1]).(2010-2011)32001001201232233)),(0,1,2)()()[,]()[,,]()()1(1)2(1)(0)232()()(1)(0)(1)232()'(i i H x f x i N x f x f x x x x f x x x x x x x x x x x x H x N x k x x x x x k x x H ===+-+--=++++-=++=++--=+++-解答：(1),满足插值条件（（的二次插值多项式为：也可用拉格朗日插值法满足题设插值条件的插值多项式为：2323(4)23443)43(31)'(0)'(0)3()232()(2),(1)(0)(1),(1,1)4!1||=4!496x x k x H f H x x x f R x x x M M R ζζ=++-===+++--∈-≤⨯由得：k=0故：误差（x)=则误差界（x)3,已知f(x)的函数值：f(0)=2, f(1)=4, f(2)=9, 写出二次拉格朗日插值多项式及余项.(2009-2010) 4,已知f(x)的如下函数值及导数值：f(1)=1, f(2)=2, f ’(1)=3, f(3)=9; (1),建立不超过3次的埃尔米特插值多项式；(2)计算f(1.6)的近似值；若M 4=0.5，估计f(1.6)的误差界.(已知|f (4)(x)|≤M 4).(2009-2010)5,写出满足条件H(0)=1, H(1)=0, H ’(1)=1, H(2)=1的三次插值多项式，并给出误差估计式.(2008-2009B)6,已知一组数据，求函数f(x)=0的根.(2008-2009B)x i -1 0 2 3 f(x i )-7-1177,已知f(x)的如下函数值及导数值：f(0)=1, f(1)=3, f ’(1)=1, f(2)=9, (1),建立不超过3次的埃尔米特插值多项式，写出误差估计式；(2),计算f(1.8)的近似值：若M 4=1，估计f(1.8)的误差界.(已知|f (4)(x)|≤M 4).(2007-2008) 8,已知f(x)的如下函数值及导数值：f(1)=2, f(2)=4, f ’(2)=5, f(3)=8, (1),建立不超过3次的埃尔米特插值多项式；(2),计算f(2.5)的近似值：若M 4=0.5,估计f(2.5)的误差界.(已知|f (4)(x)|≤M 4).(2006-2007） 9,已知f(x)的如下函数值表x i 0.1 0.2 0.3 0.4 f(x i )1.122.652.811.68选取合适的插值节点，用二次插值多项式计算f(0.35)的近似值.(2005-2006) 10,已知f(x)=sinx 的如下函数值表x i 1.0 1.5 2.0 sinx i0.84150.99750.9093用插值多项式计算sin1.8, 并估计误差界.(2004-2005)11,用f(x)的关于互异节点集112{}{}n ni i i i x x -==和的插值多项式g(x)和h(x)构造出关于节点集1{}ni i x =的插值多项式.(2005-2006)（课后习题）-11111121111{}(),()(){}(),()()()()))()())]()n n i i i i n n n n n n n n n n n n n n q x q x g x x x x x x x x x g A x x g x ==------=----=-解答：法一：设关于节点集x 的插值多项式为则与有共同插值节点x ，则设：q(x)=g(x)+Aw w f(x (x )由q(x )=f(x 得，w w 故：q(x)=g(x)+[f(x (x )w 法二：设q(x)=g(x)+1-122311111()()(){}()()()()(),01()=()[()()]()[()()]()()()()()()[()()]=-n n i i n n n n n n n n n n x g x h x B g x h x B x x x x x x B x x x g x h x BAx x g x h x Bq x f x h x Ah x g x x x g x h x BA B -=---=---≠----===+--Aw 由于和有共同插值节点x ，则存在常数，使得则，w 故：q(x)=g(x)+由得得1111()[()()]()n n x x x x h x g x x x ----则：q(x)=g(x)+12,(1),已知f(x)的如下函数值：f(0)=1,f(1)=3,f(3)=5,写出二次拉格朗日插值多项式L 2(x); (2),若同时已知：f ’(1)=1,用待定系数法求埃尔米特插值多项式H 3(x)； (3),当(3)(4)1|()|2|()|4,[0,3]fx fx x ≤≤≤≤∈及3时，x 不取节点，[0,3]x ∈，求32()()||()()f x H x f x L x --的上界.(2011-2012)题型三：最佳平方逼近多项式及最小二乘法1,已知函数值表：x -2 -1 0 1 2 y121用二次多项式y=C 0+C 1X+C 2X 2按最小二乘法拟合改组数据,并求平方逼近误差.(2010-2011)（2005-2006）()000102030410111213142021222324012()()()()()11111()()()()()21012()()()()()4101401210,5010010010034T T T T x x x x x A x x x x x x x x x x y A AC A y c c c ϕϕϕϕϕϕϕϕϕϕϕϕϕϕϕ⎛⎫⎛⎫ ⎪ ⎪==-- ⎪ ⎪⎪ ⎪⎝⎭⎝⎭==⎛⎫⎛ ⎪ ⎪ ⎪⎝⎭⎝解答：法一：线性拟合的法方程组为：即()()01222*20000100011402583,0,3575833570581358||||=(y,y)-Y 01210402023531701(,)0,(,)(T c c y x C x xx δϕϕϕαϕαϕϕϕα⎫⎛⎫⎪ ⎪=⎪ ⎪⎪ ⎪⎭⎝⎭===-=-⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪=-=⎪ ⎪ ⎪ ⎪- ⎪ ⎪⎝⎭⎝⎭====解得:c 则平方逼近误差：法二：构造首项系数为的正交多项式：（x)=1（x)=x-111211021100002*22022220,)0(,)(,)2,()()2(,)(,)46583()()0(2)(,)514357(,)8||||=(y,y)-(,)35i i i i i i i i i x x y x x x x y ϕϕϕϕϕβααϕβϕϕϕϕϕϕϕϕϕδϕϕ======----==++-=-=∑∑（x)（x)=x 则,平方逼近误差：2,求21()1f x x=+在区间[0,1]上的一次最佳平方逼近多项式及平方逼近误差（去权函数ρ(x)=x).(2009-2010) 3,通过实验获得以下数据：x i 0 1 2 3 y i13610请用最小二乘法求形如y=a+bx 2的经验公式.(2008-2009)T T A AC A y =解析：4,利用正交多项式的性质构造首项系数为1的正交多项式1{()}i i g x ∞=,有下列公式：010111()1()()()()(),(1,2,...)k k k k k g x g x x g x x g x g x k ααβ+--==-=--=其中：111(,),(0,1,2...)(,)(,),(1,2...)(,)k k k k k k k k k k xg g k g g g g k g g αβ---====(1),求[0,1]上首项系数为1的正交多项式（权函数ρ(x)=1),g 0(x),g 1(x),g 2(x)(2),以上述正交多项式为基，求sinx 在区间[0,1]上的二次最佳平方逼近多项式，并求平方逼近误差.(2008-2009B)(2004-2005)010000110001201111211021102110000*010001(1),()1(,)11,()(,)221()(,)121(,)2()2(,)11,()()()()(,)126(,)(,)(2),()(,)(g x xdx xg g g x x x g g dx x x dx xg g g g x dx g g g x x g x g x x x g g g f g f x g g g g αααβαβϕ=====-=--===-===--=-+=+⎰⎰⎰⎰解答：21212211120020111222000222*220(,),)(,)11()sin ()sin sin 11621()()1126()()260.00746 1.09130.23546(,)||||(,)0.000623.(,)i i i i g f g g g g g x x xdx x xdx xdx x x x dx x dx x x dx x x f g f f f g g ϕ=+-+-=⋅+⋅-+⋅-+--+=-+--=-=⎰⎰⎰⎰⎰⎰∑平方逼近误差：5,以正交多项式为基，求函数21()1f x x=+在区间[0,1]上的二次最佳平方逼近多项式，并求平方逼近误差.(2007-2008)(权函数ρ(x)=x,(2011-2012))20120122201201()1,(),(),111()2,()1,()2242211112234211113454111112224561.0656,0.503x x x x x f In f f In C F In c c c In ϕϕϕπϕϕϕπ=====-=-=⎛⎫⎛⎫⎪ ⎪⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪⋅=- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭ ⎪ ⎪- ⎪ ⎪⎝⎭⎝⎭==-解答：法一：取解得，，，正规方程组为：H 即：解得：c c 2*222*00001000111110110002,0.07423() 1.06560.503020.07423=(f,f)-F 0.000029041()11(,)223,()1(,)332(,)8(,)1,(,)15(,)T n p x x x C g x xg g g x x x g g xg g g g g g g g δαααβ=-=--======-=-====c 故二次最佳平方逼近多项式：平方逼近误差：法二：构造首项系数为的正交多项式：221100*201220120011222*1882163()()()()()()15318510(,)(,)(,)()()()() 1.06560.503020.07423(,)(,)(,)=(f,f)-F 0.00002904T n g x x g x g x x x x x f g f g f g p x g x g x g x x x g g g g g g C αβδ=--=---=-+=++=--=则：平方逼近误差：6,通过实验获得以下数据：u i 0 1 9 16 v i11/21/31/4请用最小二乘法求形如011v c c u=+的经验公式，并求平方误差.(2006-2007)011:c c u v=+解答转化题型四：代数精确度1,确定参数α，使求积公式20()[(0)()]['(0)'()]2hhf x dx f f h h f f h α≈++-⎰的代数精确度尽可能高，并求其代数精确度.(2010-2011)23322442320()1,,()1(),=121()()(0)(03)2121()()0+)(04)212()[(0)()]['(0)'()]2h h h f x x f x f x x h f x x f x dx h h h h f x x f x dx h h h hf x dx f f h h f f h αα====++-=≠+-≈++-⎰⎰⎰解答：令显然成立令得又时：时：（故具有三次代数精确度.2,确定参数A 1，A 2，使求积公式12()()(0)()3hhhf x dx A f h A f f h -≈-++⎰的代数精确度尽可能高，并求其代数精确度.(2009-2010) 3,建立高斯型求积公式1211221()()()x f x dx A f x A f x -≈+⎰.(2009-2010)231212113112211224112211335112211212000010001,23025031,53()1(,)0,()(,)x A A x dx A x A x x dx A x A x x dx A x A x x dx x A A g x xg g g x x xg g ααα----+==+==+==+===-=-======-=⎰⎰⎰⎰解答：法一：已知求积公式有3次代数精确度，令f(x)=1,x,x 得解上述方程组得：x 法二：构造二次正交多项式11110110022110021211222112111221121(,)(,)30,(,)(,)53()()()()53()0,511,33133()[()()]355xg g g g g g g g g x x g x g x x g x x x x x x A x dx A x dx x x x x x f x dx f f βαβρ---=====--=-==-=---=⋅==⋅=--≈-+⎰⎰⎰令得高斯点： x 故高斯型求积公式为：方法三：设[-1,1]上权(x)2221221122122121122221122331122212121().223()0,+0,5352()0,0,053().52:3250()()(),(g x x ax b b x g x dx b a x xg x dx a g x x A A A x A x A x A x A x A x x x x x x x c x c x ϕϕ--=++===-⋅====-+=+=+=+==--=++⎰⎰=x ,首项系数为1的二次正交多项式为则有：即即所以剩下步骤同法二.法四显然222221122111122212211221112221222332211122211221112221122112)()0()()()()()()()2230,535()()()()()20,053(),5x A x A x A x c x c A x c x c A x A x c A x A x c A A c c A x x A x x A x A x c A x A x c A x A x c c x x ϕϕϕϕϕϕ==+=+++++=+++++=+==-+=+++++====-剩下步骤同法二.4,确定求积公式()()(0)()hhf x dx Af h Bf Cf h -≈-++⎰中的参数A,B,C ，使其代数精度尽量高，并指出其代数精确度.(2008-2009B) 5,确定求积公式1211123()()()()343234f x dx f f f ≈-+⎰的代数精确度.(2006-2007B) 6,确定下列求积公式中的参数，使求积公式的代数精确度尽可能高，并求出代数精确度10120113()()()()424f x dx A f A f A f ≈++⎰.(2005-2006)7,确定下列求积公式中的参数，使求积公式的代数精确度尽可能高，并求出代数精确度101()()(0)()hhf x dx A f h A f A f h --≈-++⎰.(2004-2005)8,已知h>0，建立高斯型求积公式：21122()()()hhx f x dx A f x A f x -≈+⎰.(2011-2012)题型五：求积公式的最少节点数1,设定积分32x e dx -⎰，问用复化辛普森(Simpson)求积公式进行计算，要求误差小于10-6，所需要的最少节点数为多少？(2010-2011)(4)2244(4)461(),()16301[]||()|101801801696017.0519.x xS f x e fx eb a h f h f h b ahη---==--=-≤⋅=<-=解答：复化辛普森公式截断误差：|R 解得：h<0.176,n>故应取个节点2,设定积分13x edx -⎰，问用复化梯形求积公式进行计算，要求误差小于10-6，所需要的最少节点数为多少？（2009-2010）(2)3322(2)261(),()9101[]||()|10121891622.8.x x T f x e f x e b a h f h f h b ahη---==--=-≤⋅=<-=解答：复化梯形公式截断误差：|R 解得：h<0.357,n>故应取4个节点3,给定积分2cos2xdx ⎰，问用复化梯形求积公式和复化辛普森(Simpson)求积公式进行计算，要求误差小于10-6，所需要的最少节点数各为多少？（注：2(2)4(4)[](),[](),[,]122880T S b a b a R f h f R f h f a b ηηη--=-=-∈）(2008-2009B) 4,给定积分14x edx -⎰，问用复化梯形求积公式和复化辛普森(Simpson)求积公式进行计算，要求误差小于10-6，所需要的最少节点数各为多少？（2007-2008） 5，给定积分21Inxdx ⎰，问用复化梯形求积公式和复化辛普森(Simpson)求积公式进行计算，要求误差小于10-6，所需要的最少节点数各为多少？ （已知：2(2)4(4)1212[](),[](),,(,)12180T S b a b a R f h f R f h f a b ηηηη--=-=-∈）（2006-2007） 6,用积分82122dx In x=⎰计算In2，要使所得近似值具有7位有效数字，问用复化辛普森求积公式至少需要取多少个节点？（2005-2006）4(4)8(4)52(4)-744(4)4-7[](),[2,8]18011122,(),()223|()|,[2,8]817[]102631[]||()|101801808802820.04472,S S S b a R f h f In dx f x f x x x xf x x R f b a h R f h f h h n hηηη-=-∈===≤∈≤⨯-=-≤⋅=≤⨯-≤≥=⎰解答：复化辛普森公式截断误差公式：则使所得的近似值具有位有效数字，即令：|134.2137故至少需要取个节点.7，用积分6213dx In x=⎰计算In3，要使所得近似值具有5位有效数字，问用复化梯形求积公式至少需要取多少个节点？（2004-2005） 8,对于定积分1()If x dx =⎰，当M 2=1/8,M 4=1/32,用11点的复化辛普森(Simpson)求积公式求I 的截断误差为R s [f],用n 个节点的复化梯形求积公式求I 的截断误差为R T [f],要使R T [f]≤R s [f]，n 至少是多少？（M 2=max|f ”(x)|,M 4=max|f (4)(x)|,[0,1]x ∈).(2011-2012)题型六：Doolittle 分解及方程组求解1,求矩阵212454635⎛⎫ ⎪ ⎪ ⎪-⎝⎭的Doolittle 分解.(2010-2011) 212100212454210030635321001LU ⎛⎫⎛⎫⎛⎫ ⎪ ⎪⎪== ⎪ ⎪⎪ ⎪ ⎪⎪---⎝⎭⎝⎭⎝⎭解答：A=2,求矩阵114103241⎛⎫ ⎪- ⎪ ⎪⎝⎭的Doolittle 分解.(2009-2010) 3,设线性方程组123410135114152410162116x x x x ---⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪⋅= ⎪ ⎪ ⎪- ⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭ (1),对方程组的系数矩阵A 作Doolittle 分解；(2),用所得的Doolittle 分解求该线性方程组的解.(2007-2008&2005-2006)1234123410001013101311000132114124100013224101119162116210001313191,,,)(5,0,11,)13,,,)(1,1,1,1).T TT T A LU LY b y y y UX Y x x x --⎛⎫⎛⎫--⎛⎫ ⎪⎪- ⎪ ⎪⎪-⎪=== ⎪⎪--- ⎪ ⎪⎪ ⎪ ⎪⎪-- ⎪⎪⎝⎭⎝⎭⎝⎭==---==--解答：由得：(y 由得：(x4,设线性方程组123411415101312410762118x x x x -⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪-- ⎪⎪ ⎪⋅=⎪ ⎪ ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭ (1),对方程组的系数矩阵A 作Doolittle 分解；(2),用所得的Doolittle 分解求该线性方程组的解.(2006-2007)5,设线性方程组：12312312323153478113x x x x x x x x x ++=+-=-++=-(1),对方程组的系数矩阵A 作Doolittle 分解；(2),利用上述分解结果求解该线性方程组.(2004-2005)6,用高斯顺序消去法求解线性方程组：13241234242532431737x x x x x x x x x x +=+=+++=+=.(2010-2011)432110205102051020*******101301013=124317022312002160103701037000242,2,1, 1.x x x ⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪→→⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭====解答：增广矩阵回代求解：x 7,用高斯顺序消去法求解线性方程组:1231231233472212320x x x x x x x x x -+=-+-=---=.(2009-2010)题型七：条件数及范数1,求线性方程组1212391078981510x x x x x --=+==的系数矩阵A 的条件数cond 1(A),并说明其含义.（2010-2011）1111191008900015910089010015()||||||||19193611A A cond A A A A b ----⎛⎫ ⎪= ⎪⎪⎝⎭⎛⎫ ⎪-- ⎪= ⎪⎪⎪⎝⎭==⨯=解答：系数矩阵条件数远大于，这说明当和有小扰动时会引起解的较大误差，即该方程组是病态的.2,设矩阵15000910089A ⎛⎫ ⎪=-- ⎪ ⎪⎝⎭，求cond ∞(A).（2009-2010） 3,设三阶对称矩阵A 的特征值分别为：-2,1,3，求||A||2及cond 2(A).(2007-2008)222max max max 111-122-12max max max 1222||||()()()3||||(())()=()=1()|||||||| 3.T T A A A A A A A A A A cond A A A λλλλλλ----========解答：（）则：4,若n 元线性方程组Ax=b 为病态的，可以得到关于系数矩阵A 的什么性质.(2006-2007)5,若111123124A ⎛⎫⎪= ⎪ ⎪⎝⎭,求cond 1(A).（2005-2006）求cond ∞(A).(2004-2005) 6,设1231032475A -⎛⎫ ⎪=-- ⎪ ⎪-⎝⎭,求1||||||||A A ∞与.(2007-2008)7,若1234A ⎛⎫= ⎪⎝⎭,求谱半径()A ρ.(2005-2006)5332ρ+解答：最大特征值：（A)=题型八：雅可比迭代与高斯-赛德尔迭代1，写出求解方程组1231231237321241021534818x x x x x x x x x -+=--=--=的雅可比迭代公式，并说明其收敛性.(2010-2011)(1)()()123(1)()()213(1)()()312(0)1(3212)71(4215)101(3418)87324102348.k k k k k k k k k J x x x x x x x x +++=-+=--++=--++-⎛⎫ ⎪=-- ⎪ ⎪--⎝⎭解答：雅可比迭代公式为：x 雅可比迭代法迭代矩阵：B 严格对角占优，故求解该方程组的雅可比迭代法关于任意初始向量x 收敛2，设有方程组：132********2112212x x x x x x x -=+=-++=，讨论用雅可比迭代法和高斯-赛德尔迭代法解此方程组的收敛性.(2010-2011)112330200030000202100002000121221000200020031()002110211||0,=0=-=-12J J J L D U B D L U E B B λλλλρ---⎛⎫⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪=++=++ ⎪ ⎪ ⎪ ⎪⎪ ⎪ ⎪ ⎪--⎝⎭⎝⎭⎝⎭⎝⎭⎛⎫ ⎪ ⎪⎪=-+=- ⎪ ⎪⎪- ⎪⎝⎭-=解答：A=雅可比迭代矩阵：得，（）<1,故用雅可比迭代法解答此方程组对任意(0)1123(0)20031-()00211001211||0,=012-S S S B D L U E B B λλλλρ-⎛⎫ ⎪ ⎪⎪=-+=- ⎪ ⎪⎪ ⎪⎝⎭-===初始向量x 都收敛.高斯赛德尔迭代矩阵：得，（）<1,故用高斯赛格尔迭代法解答此方程组对任意初始向量x 都收敛.3,写出求解方程组：123123123532124721535818x x x x x x x x x -+=--=--=的高斯-赛德尔迭代公式，并说明收敛性.(2009-2010)4,用雅可比迭代法求解以313132323A ⎛⎫⎪= ⎪ ⎪-⎝⎭为系数矩阵的线性方程组时，确定其收敛性.(2009-2010)5,设线性方程组123123123221162222x x x x x x x x x -+=-+-=--+=-，讨论分别用雅可比迭代法和高斯-赛德尔迭代法解此线性方程组的收敛性，若收敛，请给出迭代格式.(2008-2009B)6,设线性方程组：1231231232215202225x x x x x x x x x +-=-++=++=-(1),证明求解该方程组的雅可比迭代法关于任意初始向量收敛；相应的高斯-赛德尔迭代法不是关于任意初始向量收敛；(2),取(0)(0,0,0)T x =，用雅可比迭代法进行求解，要求(1)()5||||10k k xx +--<.(2007-2008)11231123022()101220||0,===0)1022()023002||0,0,2,)1-J J J S S S D L U E B D L U E λλλλρλλλλρ---⎛⎫ ⎪=-+=-- ⎪⎪--⎝⎭-=<-⎛⎫⎪=-+=- ⎪⎪⎝⎭-====>解答：(1):B B 解得：，（B B 解得：（B 所以用雅可比迭代法解此方程组对任意初始向量都收敛,而用高斯赛德尔迭代法解此方程组不是对任意初始向量都收敛.(2):(1)()()123(1)()()213(1)()()312(0)(1)(2)(3)(4)2215202225(0,0,0)(15,20,25)(105,60,35)(205,160,65)(205,160,65)k k k k k k k k k T T T TTx x x xx x x x x x x x x +++=-+-=--+=---==--=--=-=-雅可比迭代公式：x 当时，计算得：（精确解）.7,设线性方程组:123123123821027325431111x x x x x x x x x ++=--++=-+=-(1),写出求解该方程组的雅可比迭代法的迭代公式和高斯-赛德尔迭代法的迭代公式，并确定其收敛性； (2),取(0)(0,0,0)T x=，用高斯-赛德尔迭代法计算x(3).(2006-2007)8,设线性方程组Ax=b 的系数矩阵232131t A t t ⎛⎫⎪= ⎪ ⎪-⎝⎭，其中t<0，问t 取何值时雅可比迭代法关于任意初始向量都收敛.(2006-2007)12122223021()0310422||()0=0=-,=)12||<1,t<-2,or t>20, 2.J J J t t D L U t t t t E B t t ttt t λλλλλλρ-⎛⎫-- ⎪ ⎪ ⎪=-+=-- ⎪ ⎪ ⎪- ⎪⎝⎭-=-=<<<-解答：雅可比迭代矩阵B 得，，雅可比迭代法对于任意初始向量都收敛，则（B 即：得又故9,1),设线性方程组：121232343243430424x x x x x x x +=+-=-+=-写出求解该方程组的雅可比迭代法的迭代公式，并确定该迭代法的收敛性；2),设线性方程组：123123123104413410811481025x x x x x x x x x ++=++=++=写出求解该方程组的高斯-赛德尔迭代法的迭代公式，并确定该迭代法的收敛性.(2004-2005)10,给定方程组：1231231232251223x x x x x x x x x +-=++=++=(1),用三角分解法解此方程组；(2),写出解此方程组的雅可比迭代公式，说明收敛性；取初始向量x 0=(0,0,0)T,当21||||10k kx x -+-<时，求其解.(2011-2012)11,设()21253sin 3421sincos 4134tan 5k k k k k k k Ak k k kkk⎛⎫- ⎪+ ⎪ ⎪= ⎪+ ⎪ ⎪+ ⎪⎝⎭，求()lim k k A →∞.(2007-2008)()020lim 021205K k A →∞⎛⎫⎪= ⎪ ⎪⎝⎭解答：12,若()()11,lim 1sin sin k k k k k k AA k k k k →∞⎛⎫⎪+=⎪ ⎪⎪⎝⎭求.(2004-2005)()01lim 10K k A→∞⎛⎫= ⎪⎝⎭解答： 题型九：非线性迭代1,设计一个算法求125的值.(2008-2009B)101125(),0.2k k kx x x +=+>解答：牛顿迭代公式：x2,给出用牛顿法求6170的近似值的迭代公式，并确定初值的取值范围.(2010-2011)6661556'5"4"*600066601050517017001701170[5]66()170,()60,()300170()()0,.1170170170(5)17061170()(5)6k k k k k kx x x x x x x x f x x f x x f x x x f x f x x x x x x g x x x +=-=-=-=+=-=>=>>⋅><-=+-=+-解答：转化为方程的正根.由牛顿迭代法得迭代公式：当时，故此时收敛到当0<时，设66'6666611*60170,(0,170)1850()(5)0,(0,170),()(170)0,6:1700,170,(0,170),.0.x g x x g x g xx x x x x ∈=-<∈>=->>∈>故故回到前段.所以当迭代公式也收敛到综上：3,给出用牛顿法求5140近似值的迭代公式，并给出初值的取值范围.(2009-2010)解答：方法同上.4,设φ(x)=x+c(x 2-5),当c 为何值时，x k+1=φ(x k )，(k=0,1,2…)产生的序列{x k }收敛于5；又c 为何值时收敛最快？(2010-2011)2''**1**'*5),||<1,||<1110,=50;51.25k k cx x c ϕϕϕϕϕ+-=-<<-<<解答：(x)=x+c(x (x)=1+2cxx (x )收敛，则有(x )即1+2cx 又，则当(x )=0,即c=-时，收敛最快5,设2()(3)x x c x ϕ=+-，应如何选取常数c 才能使迭代1(),(0,1,2)k k x x k ϕ+==具有局部收敛性？C 取何值时，这个迭代收敛最快？取x 0=2,123c =-计算()x ϕ的不动点，要求当61||10k k x x -+-<时结束迭代.(2004-2005)****21*2'****'**1(),(3)3,()(3)()|133|12|1,11,3,-0,,0.333(2),()0+0,636(3),k k k x x x x c x x x x c x x cx cx x c or c x x ϕϕϕϕ++==+-=±=+-<+<-<<=±<<<<==±±解答：(1),令x 收敛于则故要局部收敛，即|又得根据收敛阶定理，当时，迭代至少二阶收敛，即12cx 得c=故c=时，迭代收敛最快.迭代公式为：2012346*431(3)2321.7113248651.7319268031.7320508041.732050808|10,: 1.732050808.k k x x x x x x x x x x -=--=====-<=又因为|故6,方程x 3-3x-1=0在x=2附近有一根，构造一个局部收敛的不动点迭代法，并说明收敛的理由.(2009-2010)3'3223132(1.5) 1.765174168,(2.5) 2.040827551[1.5,2.5]()[1.5,2.5]11()|||0.33,(13) 5.5xx x x x x ϕϕϕϕϕ+===∈∈=≤<+解答：(x)=取的邻域[1.5,2.5]当时，又因为|故迭代在[1.5,2.5]上整体收敛.7,已知方程42()440f x x x =-+=有一个两重根02x =,请以初值x 0=1.5,用m 重根的牛顿迭代法计算其近似值，要求51||10k k x x -+-<.(2008-2009B)(P204例7.7）8,(1),已知方程240xex +-=在0.6附近有一根x ，迭代法214,0,1,2kx k x ek +=-=是否局部收敛？如果不收敛，试构造一个局部收敛的不动点迭代法，并说明收敛的理由.(2),取x 0=0.6,用你所构造的不动点迭代法求解该方程，迭代至x 5. (3),给出牛顿法求120的近似值的迭代公式，并给出初值的取值范围.(2007-2008)2'2'**1'''1(1):()4,()2|()|1,(0),1(4)211(4),()22(4)1(0)2,(1)3()[0,1]21()||(1)|161(4)2x xk k k k x e x e x x x In x In x x x In In x x x In x ϕϕϕϕϕϕϕϕϕϕ++=-=->>=---=-==∈≤=<=-解答：故该迭代公式不是局部收敛的.构造：理由：取邻域[0,1](x)=故又|故迭代式在[0,1]上整体收敛11021324354101(2),(4),21(4)0.61188771521(4)0.61013645921(4)0.61039483321(4)0.61035672221(4)0.61036234421120(3),(),0.2k k k k kx In x x In x x In x x In x x In x x In x x x x x ++=-=-==-==-==-==-==+>.则9,给定方程x 2+x-2=0,[0,2]x ∈,采用迭代公式xk+1=x k +c(x k 2+x k -2),(k=0,1,2…)求其根，问当c 为何值时，迭代法收敛？又当c 为何值时，迭代法收敛最快?(2011-2012)*2'''1,()(2)()1(21)2(1)||1(21)|1,-0.31(1)=03x x x c x x x c x c c ϕϕϕϕ==++-=++=++<<<解答：当|即时，线性收敛当，即c=-时收敛最快.10,给定方程230x xe -=，[3,4]x ∈(1),构造一种线性收敛的不动点迭代公式求该方程的根（含迭代公式，初值取何值或何区间，迭代法收敛的原因）； (2),构造一种二次收敛的不动点迭代公式求该方程的根（含迭代公式，初值取何值或何区间，迭代法收敛的原因).(2011-2012)21111'12102'"0(1),()(3),3.29(3)()(4) 3.8712(),[3,4]23(3),(0,1,2,)[3,4].(2),()3,[3,4](3)0,(4)0()60,()60,[3,4]3k k x x x x In x x x x In x k x f x x e x f f f x x e f x e x x ϕϕϕϕϕ+==≤≤=≤≤∈==∈=-∈><=-<=-<∈=解答：故不动点迭代公式:x 对于任意初值收敛取初值时，牛顿213.6kkx kk k x k x ex x x e+-=--迭代法：收敛，且二次收敛11,方程x 3-x 2-1=0在x=1.5附近有根,建立一个收敛的迭代公式，并证明其收敛性.(2004-2005)122''33312111.51()1(1.3) 1.591715976,(1.6) 1.390625[1.3,1.6]()[1.3,1.6]222(),|()|||0.921.311k k k kx x x x x x x x x x x x x ϕϕϕϕϕϕ++=+==+==∈∈=-=-≤<=+解答：取的邻域[1.3,1.6]故当时，又故迭代公式：在[1.3,1.5]上整体收敛.12,(1),已知方程1020x e x +-=在0.09附近有一根x,迭代法1(210),(0,1,2)k k x In x k +=-=是否局部收敛？如果不收敛，请构造一个局部收敛的不动点迭代法，并说明收敛的理由；(2),取x 0=0.09,用局部收敛的迭代法计算x 5； (3),用牛顿法求3234的近似值，并给出初值的取值.(2006-2007)'''*1''5(1),()(210),()15|()|1,[0,1],|()|>1.11510111(),()51010(0)0.1,(0.12)0.087250323[0,0.12]()[0,0.12]()|kx k x xx In x x xx x x x e x e x e x x x ϕϕϕϕϕϕϕϕϕϕ+-=-=->∈=-=-=-==∈∈≤解答：显然故该迭代公式不是局部收敛的构造：因为取[0,0.12]邻域考察故当时，又|'0.12110.09010.09058257820.09051881530.0905241|(0.12)|||0.1131101151011(2),510110.09,0.090582578510110.090518815510110.09052579651011510k kx k x k e x e x e x x e x e x e x e ϕ++=-<<=-=-==-==-==-==-故迭代公式：在[0,0.12]上整体收敛.57960.09052503151200.090525031110.0905251155102117(3),()30.k k k x e x x x +==-==+>使用迭代公式：进行求解.初值:x13,设方程x 3-3x-1=0在x=2附近有根;1),证明该方程在区间[1.5,2.5]内有唯一根x *;2),确定迭代函数φ(x).当初始值x 0在何区间取值时，迭代公式x k+1=φ(x k ),(k=0,1,2…)收敛到x *,并说明理由. 3),写出求解该方程组的牛顿法迭代公式，当初始值x 0在何区间取值时，牛顿法迭代公式收敛到x,并说明理由.取x 0=1.8,用牛顿法迭代公式计算x,要求(1)()4||||10k k x x +--<.4),写出求解该方程的弦截法迭代公式，当初始值在何区间取值时，弦截法迭代公式收敛到x,并说明理由.(2005-2006)3'2'331223(1),()31,()33(1.5) 2.125,(2.5)7.125(1.5)(2.5)0,()0()0,[1.5,2.5][1.5,2.5].(2),3121(3),,3333()3k k k k k k k f x x x f x x f f f f f x f x x x x x x x x x f x x +=--=-=-=⋅<=>∈--+=-=--=-解答：证明：故在[1.5,2.5]内有根.又故方程在区间内有唯一根牛顿法迭代公式：'2"1,()33,()6x f x x f x x-=-=题型十:稳定算法1,对给定的x ，下列两式能否直接计算，说明理由；如果不能，请给出变换算式：（1）21x x +-，x 很大；（2）311x +-，|x|很小.(2010-2011)223331(1):1111=.1+1x x x x x x x +-=+++-+解答：不能直接计算，因为两个相近的数相减，会产生较大的误差：；2,为了提高计算精度，当正数x 很大时，计算1x x +-时应转化成什么形式.（2005-2006）3,给出计算积分1,(0,1,2,10)10nnx I dx n x ==+⎰的递推稳定算法和初值.(2010-2011) 1111111000-11110002010101101010101=101011111)11101010(1)11121[].2111)101)220(1)n n n n n n n n n n x x dx x dx x dx I I x x nn n x x x dx dx dx n x n n n n ----+-===-=-++-=<<=+++=+=+++⎰⎰⎰⎰⎰⎰解：I 该算法不稳定，变形得：I 因为（取初值I （（4,设计一种求1x n nI e x dx =⎰（n 为非负整数）稳定的递推算法，包括递推公式，初值的确定；当初值201221e I =⋅时，利用上述稳定的递推公式计算三个连续的积分值.(2011-2012)题型十一:部分证明题1,利用差分的性质证明：12+22+…n 2=n(n+1)(2n+1)/6222()12,g n n n =++证明：设函数对任意的建立差分表：g(n)(n+1)22n+3 2 g(n+1) (n+2)2 2n+5 2 g(n+2) (n+3)2 2n+7 g(n+3) (n+4)2 g(n+4)函数g(n)的三阶差分是与n 无关的非零常数，故g(n)是n 的三次多项式：3(1)1,(2)5,(3)14,(4)30111()()14521231(1)(2)(1)(2)(3)(1)(21)14521!2!3!6g g g g n n n g n N n n n n n n n n n n ====---⎛⎫⎛⎫⎛⎫==+++ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭------++=+⋅+⋅+⋅=按等距节点牛顿向前插值公式建立三次插值多项式，则2,证明：n+1个互异节点的插值型求积公式的代数精确度至少为n.(2010-2011)(1)0()(),.(1)!n nbi ai f x x dx n ζ+=-+∏⎰证明：截断误差R[f]=易证 3,若0{()}ni i l x =是关于互异节点0{}ni i x =的拉格朗日插值基函数组，函数0011()()()(),(1)n n f x x l x x l x x l x n =++≥，证明：f(x)≡x.（2009-2010）00110()()()()()()()()n n i i n n i f x L x f x l x x l x x l x x l x f x x=≈==+++≡∑证明：故：4,证明：0101'()[()()]"()2hf x f x f x f h ζ=--，其中h=x 1-x 0,01(,)x x ζ∈.(2009-2010)"'20000"'211001010'"010())()()()2!(),())()()()2!1()[()()]()2f f x x x x x f x x f x f x x x x x hf x f x f x f h ζζζ+-+-==+-+-=--证明：由泰勒公式得f(x)=f(x 令则f(x 整理得： 5,证明：关于互异节点0{}ni i x =的拉格朗日插值基函数0{()}ni i l x =满足恒等式012()()()()1n l x l x l x l x +++≡.(2008-2009B)(2006-2007B)(2004-2005)120(1)(1)1010()1,(),,1=L ()()()()()()()1,()0,()()0(1)!()()()()1n n n n i n i n n n n ni n i f x f x x x x x R x l x f x R x f f x fx R x W x n l x l x l x l x ζ=+++==+=+=≡==+=+++≡∑∑证明：令对在上进行拉格朗日插值，有因故故：6,证明求积公式()[()()]2bab af x dx f a f b -≈+⎰的截断误差：3"()[](),12f R f b a ηη=--∈其中：（a,b).(2007-2008) (1)001(2)(2)(2)33()()(1)!1,,()()()1"()()()()()()()2!2!2!612n nb i ai b b aa f x x dx n n x a xb f f f f x a x b dx x a x b dx a b b a ζζηηη+=-+===--=--=⋅-=--∏⎰⎰⎰证明：插值型求积公式截断误差R[f]=R[f]=7,设矩阵A 为可逆上三角阵，证明A -1仍为上三角阵，并导出求逆算法.(2006-2007B)8,设x k =a+kh(k=0,1,2;h>0),f(x)的三阶导数连续，证明：2(3)102021'()[()()](),(,)26h f x f x f x f x x h ζζ=-+-∈其中为中值.(2011-2012)001122120201201201021012202112020101222,),,),,)()()()()()()()()()()()()()()()()()()()()()()()()(22x y x y x y x x x x x x x x x x x x x f x f x f x x x x x x x x x x x x x x x x x x x x x x x x x f x f x f x h h h ------=++------------=-+证明：过（（（的拉格朗日插值多项式为：L 12'2102(3)201202(3)'''1210122'(3)10202)1()[()()]2()()()()()(),(,)3!()()()[()()()]3!1()[()()](),(,)26x x L x f x f x hf f x L x x x x x x x x x f f x L x x x x x x x h f x f x f x f x x h ηηηζζ==-+-=---∈-=---=-+-∈又故：。
中国石油大学《数值分析》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)
【免费下载】江西理工大学大学数值分析计算方法期末试卷3及参考答案
姓名
学号
------------------------------密-----------------------------封-------------------------- -线------------------------
课程
江西理工大学 大 学
2010 至 2011 学年第 一 学期试卷
2. 数值积分公式形如（15）
1
0
f
( x)dx
A0
f
(0)
(1) 试确定求积公式中的参数 A0 , A1, B0 ，使其代数精度尽可能高.
并求出其代数精度。
(2) 已知该求积公式余项 R[ f ] kf ''' ( ), (0,1), 试求出余项中的
参数 k 。
A1
f
(1)
B0
，
(
A)
)(xFra bibliotekx0
)
2
进行近似
， f [0,1]
。
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电，力根通保据过护生管高产线中工敷资艺设料高技试中术卷资，配料不置试仅技卷可术要以是求解指，决机对吊组电顶在气层进设配行备置继进不电行规保空范护载高与中带资负料荷试下卷高问总中题体资，配料而置试且时卷可，调保需控障要试各在验类最；管大对路限设习度备题内进到来行位确调。保整在机使管组其路高在敷中正设资常过料工程试况中卷下，安与要全过加，度强并工看且作护尽下关可都于能可管地以路缩正高小常中故工资障作料高；试中对卷资于连料继接试电管卷保口破护处坏进理范行高围整中，核资或对料者定试对值卷某，弯些审扁异核度常与固高校定中对盒资图位料纸置试，.卷保编工护写况层复进防杂行腐设自跨备动接与处地装理线置，弯高尤曲中其半资要径料避标试免高卷错等调误，试高要方中求案资技，料术编试5交写、卷底重电保。要气护管设设装线备备置4敷高、调动设中电试作技资气高，术料课中并3中试、件资且包卷管中料拒含试路调试绝线验敷试卷动槽方设技作、案技术，管以术来架及避等系免多统不项启必方动要式方高，案中为；资解对料决整试高套卷中启突语动然文过停电程机气中。课高因件中此中资，管料电壁试力薄卷高、电中接气资口设料不备试严进卷等行保问调护题试装，工置合作调理并试利且技用进术管行，线过要敷关求设运电技行力术高保。中护线资装缆料置敷试做设卷到原技准则术确：指灵在导活分。。线对对盒于于处调差，试动当过保不程护同中装电高置压中高回资中路料资交试料叉卷试时技卷，术调应问试采题技用，术金作是属为指隔调发板试电进人机行员一隔，变开需压处要器理在组；事在同前发一掌生线握内槽图部内纸故，资障强料时电、，回设需路备要须制进同造行时厂外切家部断出电习具源题高高电中中源资资，料料线试试缆卷卷敷试切设验除完报从毕告而，与采要相用进关高行技中检术资查资料和料试检，卷测并主处且要理了保。解护现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
北京理工大学2009级数值分析试题及答案
课程编号：12000044 北京理工大学2010-2011学年第一学期2009级计算机学院《数值分析》期末试卷A 卷班级 学号 姓名 成绩注意：① 答题方式为闭卷。
② 可以使用计算器。
请将填空题和选择题的答案直接填在试卷上，计算题答在答题纸上。
一、 填空题 （2 0×2′）1. 设x =0.231是精确值x *=0.229的近似值，则x 有 位有效数字。
2. 设⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡-=32,1223X A ，‖A ‖∞＝___ ____，‖X ‖∞＝__ _____，‖AX ‖∞≤____ ___ (注意：不计算‖AX ‖∞的值) 。
3. 非线性方程f (x )=0的迭代函数x =ϕ(x )在有解区间满足 ，则使用该迭代函数的迭代解法一定是局部收敛的。
4. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= ，f [20,21,22,23,24,25,26,27,28]= 。
5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 阶的连续导数。
6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商公式的 （填写前插公式、后插公式或中心差分公式），若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 （填写前插公式、后插公式或中心差分公式）；如果要估计结果的舍入误差，应该选用插值公式中的 。
7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=ni i x a 0)( ；所以当系数a i (x )满足 ，计算时不会放大f (x i )的误差。
8. 要使20的近似值的相对误差小于0.1%，至少要取 位有效数字。
9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收敛于方程组的精确解x *的充分必要条件是 。
数值分析报告期末考试复习题及其问题详解
数值分析期末考试复习题及其答案1. 已知325413.0,325413*2*1==X X 都有6位有效数字，求绝对误差限。
（4分）解：由已知可知,n=65.01021,0,6,10325413.0016*1=⨯==-=⨯=ε绝对误差限n k k X 2分 620*21021,6,0,10325413.0-⨯=-=-=⨯=ε绝对误差限n k k X 2分2. 已知⎢⎢⎢⎣⎡=001A 220- ⎥⎥⎥⎦⎤440求21,,A A A ∞ （6分） 解：{},88,4,1max 1==A 1分 {},66,6,1max ==∞A 1分 ()A A A T max 2λ= 1分⎢⎢⎢⎣⎡=001A A T 420 ⎥⎥⎥⎦⎤-420⎢⎢⎢⎣⎡001 220- ⎥⎥⎥⎦⎤440=⎢⎢⎢⎣⎡001 080 ⎥⎥⎥⎦⎤3200 2分 {}3232,8,1max )(max ==A A T λ 1分 24322==A3. 设32)()(a x x f -= （6分） ① 写出f(x)=0解的Newton 迭代格式② 当a 为何值时，)(1k k x x ϕ=+ （k=0,1……）产生的序列{}k x 收敛于2解：①Newton 迭代格式为：xa x x x ax a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(22321+=+=---=-=+ϕ 3分②时迭代收敛即当222,11210)2(',665)('2<<-<-=-=a a x a x ϕϕ 3分4. 给定线性方程组Ax=b ，其中：⎢⎣⎡=13A ⎥⎦⎤22，⎥⎦⎤⎢⎣⎡-=13b 用迭代公式)()()()1(k k k Ax b x x -+=+α（k=0,1……）求解Ax=b ，问取什么实数α，可使迭代收敛（8分）解：所给迭代公式的迭代矩阵为⎥⎦⎤--⎢⎣⎡--=-=ααααα21231A I B 2分其特征方程为 0)21(2)31(=----=-αλαααλλB I 2分即，解得αλαλ41,121-=-= 2分 要使其满足题意，须使1)(<B ρ，当且仅当5.00<<α 2分5. 设方程Ax=b ，其中⎢⎢⎢⎣⎡=211A 212 ⎥⎥⎥⎦⎤-112，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=765b 试讨论解此方程的Jacobi 迭代法的收敛性，并建立Gauss-Seidel 迭代格式 （9分）解：U D L A ++=⎢⎢⎢⎣⎡--=+-=-210)(1U L D B J 202-- ⎥⎥⎥⎦⎤-012 3分0,03213=====-λλλλλJ B I 2分即10)(<=J B ρ，由此可知Jacobi 迭代收敛 1分 Gauss-Seidel 迭代格式：⎪⎩⎪⎨⎧--=--=+-=++++++)1(2)1(1)1(3)(3)1(1)1(2)(3)(2)1(12276225k k k k k k k k k x x x x x x x x x （k=0,1,2,3……） 3分6. 用Doolittle 分解计算下列3个线性代数方程组：i i b Ax =（i=1,2,3）其中⎢⎢⎢⎣⎡=222A 331 ⎥⎥⎥⎦⎤421，23121,,974x b x b b ==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡= （12分）解：①11b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=9741x A=⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211=LU 3分 由Ly=b1，即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡974 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡234 1分 由Ux1=y ，即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡234 得x1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 2分 ②22b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 由Ly=b2=x1，即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001 1分 由Ux2=y ，即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001 得x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 2分 ③33b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0由Ly=b3=x2，即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-05.05.0 1分 由Ux3=y ，即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-05.05.0 得x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-025.0375.0 2分7. 已知函数y=f(x)有关数据如下：要求一次数不超过3的H 插值多项式，使'11'33)(,)(y x H y x H i i == （6分）解：作重点的差分表，如下：3分21021101011001003))(](,,,[))(](,,[)](,[][)(x x x x x x x x f x x x x x x x f x x x x f x f x H --+--+-+= =-1+(x+1)-x(x+1)+2x.x(x+1)=232x x + 3分8. 有如下函数表：试计算此列表函数的差分表，并利用Newton 前插公式给出它的插值多项式 （7分）解：由已知条件可作差分表，3分i ih x x i =+=0 （i=0,1,2,3）为等距插值节点，则Newton 向前插值公式为： 033210022100003!3))()((!2))((!1)()(f h x x x x x x f h x x x x f h x x f x N ∆---+∆--+∆-+==4+5x+x(x-1)=442++x x 4分9. 求f(x)=x 在[-1,1]上的二次最佳平方逼近多项式)(2x P ，并求出平方误差 （8分）解：令22102)(x a x a a x P ++= 2分取m=1, n=x, k=2x ,计算得： (m,m)=dx ⎰-111=0 (m,n)=dx x ⎰-11=1 (m,k)= dx x ⎰-112=0(n,k)= dx x ⎰-113=0.5 (k,k)= dx x ⎰-114=0 (m,y)= dx x ⎰-11=1(n,y)=dx x⎰-112=0 (k,y)= dx x ⎰-113=0.5得方程组：⎪⎩⎪⎨⎧==+=5.05.005.011201a a a a 3分解之得c a a c a 2,1,210-=== （c 为任意实数，且不为零）即二次最佳平方逼近多项式222)(cx x c x P -+= 1分 平方误差：32),(22222222=-=-=∑=i i i y a fp f ϕδ 2分10. 已知如下数据：用复合梯形公式，复合Simpson 公式计算⎰+=10214dx x π的近似值（保留小数点后三位） （8分）解：用复合梯形公式：)}1()]87()43()85()21()83()41()81([2)0({1618f f f f f f f f f T ++++++++==3.139 4分用复合Simpson 公式： )}1()]43()21()41([2)]87()85()83()81([4)0({2414f f f f f f f f f S ++++++++==3.142 4分11. 计算积分⎰=20sin πxdx I ，若用复合Simpson 公式要使误差不超过51021-⨯，问区间]2,0[π要分为多少等分?若改用复合梯形公式达到同样精确度，区间]2,0[π应分为多少等分? （10分）解： ①由Simpson 公式余项及x x f x x f sin )(,sin )()4(==得544)4(2041021)1()4(360)(max )4(1802)(-≤≤⨯≤=≤n x f n f R x n πππππ 2分即08.5,6654≥≥n n ，取n=6 2分即区间]2,0[π分为12等分可使误差不超过51021-⨯ 1分②对梯形公式同样1)(''max 20≤≤≤x f x π，由余项公式得51021)2(122)(-⨯≤≤n f R n ππ2分即255,2.254=≥n n 取 2分即区间]2,0[π分为510等分可使误差不超过51021-⨯ 1分12. 用改进Euler 格式求解初值问题：⎩⎨⎧==++1)1(0sin 2'y x y y y 要求取步长h 为0.1，计算y（1.1）的近似值 （保留小数点后三位）[提示：sin1=0.84,sin1.1=0.89] （6分）解：改进Euler 格式为：⎪⎩⎪⎨⎧++=+=+-++-+)],(),([2),(1111n n n n n n n n n n y x f y x f hy y y x hf y y 2分 于是有⎪⎩⎪⎨⎧+++-=+-=+-++-+-+)sin sin (05.0)sin (1.012112121n n n n n n n n n n n n n x y y x y y y y x y y y y （n=0,1,2……） 2分 由y(1)=0y =1，计算得⎪⎩⎪⎨⎧=≈=+-=-838.0)1.1(816.0)1sin 11(1.01121y y y 2分 即y(1.1)的近似值为0.83813. ][],[],,[lim ],[),,(],,[)(0'000000'x f x x f x x f x x f b a x b a C x f x x ==∈∈→证明：定义：设（4分）证明：]['],[],[],[lim ][][lim]['00000000000x f x x f x x f x x f x x x f x f x f x x x x ===--=→→故可证出 4分14. 证明：设nn RA ⨯∈，⋅为任意矩阵范数，则A A ≤)(ρ （6分）证明：设λ为A 的按模最大特征值，x 为相对应的特征向量，则有Ax=λx 1分 且λρ=)(A ，若λ是实数，则x 也是实数，得Ax x =λ 1分而x x ⋅=λλ x A x ,⋅≤⋅⋅≤λ故x A Ax 2分由于A x 0x ≤≠λ得到，两边除以 1分故A A ≤)(ρ 1分 当λ是复数时，一般来说x 也是复数，上述结论依旧成立。
数值分析期末复习题
一、填空题1.设真值x=983350,则其近似值y=98000的有效数字的位数 ，绝对误差为 ， 相对误差为 。
2.x=0.1062,y=0.947,计算x+y 其有效数字的位数为 。
3.对f(x)=x 3+x+1,差商f[0,1,2,3]= ;f[0,1,2,3,4]= 。
4.设f(x)可微，求方程x=f(x)根的牛顿迭代法格式是 。
5.设方程x=ϕ(x)有根x *,且设ϕ(x)在含x *的区间(a,b)内可导，设x 0∈(a,b)则迭代格式x k+1=ϕ(x k )收敛的充要条件为 。
6.求解线性方程组Ax=b 的迭代格式x (k+1)=Jx (k)+f 收敛的充要条件为 。
7.⎪⎪⎭⎫ ⎝⎛=011001001001....A ,||A||∝= ,cond(A)∝= 。
8.n 次Legendre 多项式的最高次项系数为 。
9.中矩形公式：)()2()(a b b a f dx x f b a -+=⎰的代数精度为 。
10.求积公式：)1(21)0()(10f f dx x f '+≈⎰的代数精度为 。
11.在区间[1，2]上满足插值条件⎩⎨⎧==3)2(1)1(P P 的一次多项式P(x)= 。
12.设∑==n k k k n x f A f I 0)()(是函数f(x)在区间[a,b]上的插值型型求积公式，则 ∑=n k k A0= 。
13．梯形公式和改进的Euler 公式都是 阶精度的。
二、计算题1.利用矩阵的高斯消元法，解方程组⎪⎩⎪⎨⎧=++=++=++2053182521432321321321x x x x x xx x x2.设有函数值表试求各阶差商，并写出Newton 插值多项式。
3.求解超定方程组⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛43231211121x x的最小二乘解。
4.给定下列函数值表：求3次自然样条插值函数5.给定x x f =)(在x=100, 121, 144 三点处的值，试以这三点建立f(x)的二次（抛物）插值公式，利用插值公式求115的近似值并估计误差。
数值分析试题及答案
数值分析试题及答案一、单项选择题（每题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. 求解线性方程组时，迭代法的基本思想是从一个初始近似解出发，通过不断________来逼近精确解。
数值分析整理版试题及答案
例1、已知函数表求f（x）的Lagrange 次插值多项式和Newton〔次插值多项式。
解：（1）由题可知H 3 33H z3 3插值基函数分别为/ （「曲）（兀_七）（x_l）（x_2）°「（mF）（心-勺厂（-1-1）（-1-2）=*（尤一1）（尤一2）A（x）=(x_.lo)(x_£)(X] -曲)("一勺)("1)(—2)(1 + 1)(1一2)，2（X）=(x_")(x_xj _(x+l)(xT)_l(勺-心)(七 -")(2 + 1)(2-1) 3故所求二次拉格朗日插值多项式为3劝=2＞山(“)k=0■ ■=-3x丄(兀-l)(x-2) + 0x - —(x + l)(x-2) +4x—(x + l)(x—1)6 2 34,=-1(X-1)(X-2)+^(A +1)(X-1)57 3 7=+ —X- —6 2 3（2）一阶均差、二阶均差分别为-3-0 3八g卜A._X1_!_! 2亦小号半么齢43_均差表为故所求Newton 二次插值多项式为35= -3 + -(X + 1) + -(X + 1)(A -1)5 2 3 7 =-x +—兀一一6 2 3例2、 设/(X ) =A -2+3X + 2, Ae[0,l],试求/(x)在[0, 1]上关于p(x) = 1 , 6=span{l,x}的最佳平方逼近多项式。
解：若 O = span {!.%},则 0o (x) = l, (p x (x) = x ,且 p(x) = \ T 这样,有1S) ®) ) = J Idx = 1’ o所以，法方程为 ■■ 1r■ ■231 1■ ■23 1 2一 ~6,经过消元得2~6 1 190丄1 _2 3..4.一 12J_3.再回代解该方程，得到5=4, n 0=H6 故，所求最佳平方逼近多项式为S ；(x) = - + 4x6例|3、 设 f(x) = e x, A e[0,l],试求/⑴在[0, 1]上关于/7(A) = 1 , 6=span{l,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 设方程组，试考察解此方程组的雅可比迭代法及高斯－赛德尔迭代法的收敛性。