数值分析(计算方法)期末试卷及参考答案
(完整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 极限存在。
数值分析、计算方法试题库及答案
19.区间[a,b]上的三次样条插值函数S(x)在[a,b]上具有直到2阶的连续导数。
20.当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商公式的前插公式，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的后插公式；如果要估计结果的舍入误差，应该选用插值公式中的拉格朗日插值公式。
2、用牛顿——埃尔米特插值法求满足下列表中插值条件的四次插值多项式P4(x)，并写出其截断误差的表达式(设f(x)在插值区间上具有直到五阶连续导数)。
xi
0
1
2
f(xi)
1
-1
3
f’(xi)
1
5
3、对下面的线性方程组变化为等价的线性方程组，使之应用雅克比迭代法和高斯——赛德尔迭代法均收敛，写出变化后的线性方程组及雅克比迭代法和高斯——赛德尔迭代法的迭代公式，并简单说明收敛的理由。
16、解线性方程组的的平方根直接解法适用于任何线性方程组AX＝b。(×)
17、迭代解法的舍入误差估计要从第一步迭代计算的舍入误差开始估计,直到最后一步迭代计算的舍入误差。(×)
18、数值计算中的总误差如果只考虑截断误差和舍入误差，则误差的最佳分配原则是截断误差＝舍入误差。()
10、插值计算中避免外插是为了减少舍入误差。(×)
寂涯网络xx~xx学年第1学期 《计算方法》课程试卷A 第2页 共4页
四、设方程 在 [ ] 内有实根 ，试写出迭代公式
使 ，并说明迭代公式的收敛性。（10分）
五、设有线性方程组 ，其中
（1）求 分解;(2)求方程组的解(3)判断矩阵 的正定性（14分）
数值计算期末试题及答案
数值计算期末试题及答案1. 题目：求方程 f(x) = x^3 - 3x^2 + 2x - 1 = 0 在区间 [0, 2] 上的根。
解答：为了求解方程 f(x) = 0 在给定区间上的根，可以使用二分法或者牛顿法等数值方法。
这里我将采用二分法进行求解。
首先，观察方程在区间 [0, 2] 上的图像，可以发现 f(0) = -1，f(2) = 1，即方程在区间 [0, 2] 上存在根。
接下来，我们可以通过二分法逼近此根的位置。
二分法的基本思路是不断将给定区间一分为二，并判断根的位置在前半部分还是后半部分，然后继续在包含根的那一半区间内进行二分，直到达到所需的精确度为止。
具体的二分法迭代过程如下：1. 初始化区间左边界 a = 0，右边界 b = 2，以及精确度 eps。
2. 当 (b - a) > eps 时，执行以下步骤：a. 计算区间中点 c = (a + b) / 2。
b. 如果 f(c) 等于 0 或者在所需的精确度 eps 内，返回 c。
c. 否则，根据 f(c) 和 f(a) 的符号判断根的位置：- 如果 f(c) 和 f(a) 的符号相同，说明根在区间 [c, b] 中，更新 a = c。
- 否则，根在区间 [a, c] 中，更新 b = c。
3. 返回最终得到的近似根 c。
根据上述算法，我们可以得到方程 f(x) = 0 在区间 [0, 2] 上的近似根为c ≈ 1.521。
2. 题目：使用梯形法则计算定积分∫[0, π] sin(x) dx。
解答：定积分的数值计算可以通过数值积分方法来实现。
其中，梯形法则是一种常用的数值积分方法。
梯形法则的基本思路是将定积分区间划分成多个小梯形，然后计算各个小梯形的面积之和作为近似解。
具体的步骤如下：1. 初始化定积分区间的左边界 a = 0，右边界b = π，以及划分的小梯形数量 n。
2. 计算每个小梯形的宽度 h = (b - a) / n。
数值分析期末复习题答案
数值分析期末复习题答案一、选择题1. 以下哪个算法是用于求解线性方程组的直接方法？A. 牛顿法B. 高斯消元法C. 共轭梯度法D. 辛普森积分法答案：B2. 插值法中，拉格朗日插值法和牛顿插值法的主要区别是什么？A. 插值点的选取不同B. 插值多项式的构造方式不同C. 计算复杂度不同D. 适用的函数类型不同答案：B3. 在数值积分中，梯形法则和辛普森法则的主要区别是什么？A. 精度不同B. 适用的积分区间不同C. 计算方法不同D. 稳定性不同答案：A二、简答题1. 解释什么是数值稳定性，并举例说明。
答案：数值稳定性指的是数值方法在计算过程中对于舍入误差的敏感程度。
例如，在求解线性方程组时，如果系数矩阵的条件数很大，则该方程组的数值解对舍入误差非常敏感，即数值稳定性差。
2. 说明数值微分与数值积分的区别。
答案：数值微分是估计函数在某一点的导数，而数值积分是估计函数在某个区间上的积分。
数值微分通常用于求解函数的局部变化率，而数值积分用于求解函数在一定区间内的累积效果。
三、计算题1. 给定一组数据点：(1, 2), (2, 3), (3, 5), (4, 6)，请使用拉格朗日插值法构造一个三次插值多项式。
答案：首先写出拉格朗日插值基函数，然后根据数据点构造插值多项式。
具体计算过程略。
2. 给定函数 f(x) = x^2，使用牛顿-科特斯公式中的辛普森积分法在区间 [0, 1] 上估计积分值。
答案：首先确定区间划分，然后应用辛普森积分公式进行计算。
具体计算过程略。
四、论述题1. 论述数值分析中误差的来源及其控制方法。
答案：误差主要来源于舍入误差和截断误差。
舍入误差是由于计算机在进行浮点数运算时的精度限制造成的，而截断误差是由于数值方法的近似性质导致的。
控制误差的方法包括使用高精度的数据类型、选择合适的数值方法、增加计算步骤等。
五、综合应用题1. 给定一个线性方程组 Ax = b，其中 A 是一个 3x3 的矩阵，b 是一个列向量。
数值分析计算方法试题集及答案
数值分析复习试题第一章 绪论 一. 填空题 1.*x为精确值x 的近似值；()**x f y =为一元函数()x f y =1的近似值；()**,*y x f y =为二元函数()y x f y ,2=的近似值，请写出下面的公式：**e x x =-：***r x xe x -=()()()*'1**y f x x εε≈⋅ ()()()()'***1**r r x f x y x f x εε≈⋅()()()()()**,**,*2**f x y f x y y x y x yεεε∂∂≈⋅+⋅∂∂()()()()()****,***,**222r f x y e x f x y e y y x y y y ε∂∂≈⋅+⋅∂∂ 2、 计算方法实际计算时，对数据只能取有限位表示，这时所产生的误差叫 舍入误差 。
3、 分别用2.718281，2.718282作数e 的近似值，则其有效数字分别有6 位和7 位；又取 1.73≈-211.73 10 2≤⨯。
4、 设121.216, 3.654x x ==均具有3位有效数字，则12x x 的相对误差限为0.0055 。
5、 设121.216, 3.654x x ==均具有3位有效数字，则12x x +的误差限为0.01 。
6、 已知近似值 2.4560A x =是由真值T x 经四舍五入得到,则相对误差限为0.0000204 .7、递推公式,⎧⎪⎨⎪⎩0n n-1y =y =10y -1,n =1,2,如果取0 1.41y ≈作计算,则计算到10y 时,误差为8110 2⨯;这个计算公式数值稳定不稳定 不稳定 . 8、精确值 14159265.3*=π，则近似值141.3*1=π和1415.3*2=π分别有 3位和 4 位有效数字。
9、若*2.71828x e x =≈=,则x 有 6 位有效数字，其绝对误差限为1/2*10-5。
10、 设x*的相对误差为2％，求(x*)n的相对误差0.02n11、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字；12、计算方法主要研究( 截断 )误差和( 舍入 )误差； 13、为了使计算 ()()2334610111y x x x =++---- 的乘除法次数尽量地少，应将该表达式改写为11,))64(3(10-=-++=x t t t t y ，为了减少舍入误差，应将表达式19992001-改写为199920012+。
数值分析期末考试题带答案
湖北民族学院2012年秋季期末试卷A或BA卷课程数值分析使用班级0210403、4、5、6 制卷份数86 考生姓名命题人刘波课程负责人单位审核人答题纸数班级题号一二三四五六七八九十合计学号评分分数阅卷人注意：所有答案必须填写在答题纸上！ 一、填空题（4分⨯10=40分）1、向量T x )3,2,1(-=的范数1x = ，∞x = ，2x 。
2、已知,3)2(,1)1(==f f 那么)(x f y =以2,1=x 为节点的拉格朗日线性差值多项式为 。
3、设矩阵A 是对称正定矩阵，则用 迭代法接线性方程组,b AX =其迭代解数列一定收敛。
4、辛普森公式： 。
5、牛顿-柯特斯求积分公式的系数和=∑=nk n k C 0)( 。
6、,1)(2+=x x f 则=]3,2,1[f ,=]4,3,2,1[f 。
7、积分公式)42(32)21(31)41(32)(10f f f dx x f +-≈⎰具有 次代数精度。
二、计算题（10分⨯3=30分） 1、求01162=+-x x 的小正根。
2、给定形如)0()1()0()('01010f B f A f A dx x f ++≈⎰的求积公式，试确定系数,,,010B A A 使公式具有尽可能高的代数精确度。
3、求⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=242422221A 的特征值及普半径。
三、证明题（20分⨯1=20分） 1、用直接三角分解法解⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛201814513252321321x x x四、讨论题（10分⨯1=10分）1、用4点（n=3）的高斯——勒让德求积公式计算xdx x I cos 22⎰=π答案：一：1: 6，3，14 解： ∞x=||max 1i ni x ≤≤；1x =∑=n i i x 1||；2x =2112)(∑=ni ix ；向量的p 范数：p x =pni p ix 11)||(∑=2： 2x-1 3、高斯-赛德尔4、)]()2(4)([6)(b f b a f a f a b dx x f b a +++-≈⎰5、16、1,07、3二：1：解：6381+=x ，*2206.094.78638x x ==-≈-=,*2x 只有一位有效数字，若改用0627.094.15163816382≈≈+=-=x ，具有三位有效数字。
数值分析试题（卷）与答案解析
数值分析试题（卷）与答案解析数值分析试题一、填空题（2 0×2′）1.322A1, X23设 x=0.231 是精确值 x*=0.229 的近似值，则 x 有 2 位有效数字。
2. 若 f(x)=x7－ x3＋ 1 ，则 f[20 ,21,2 2,23 ,24,25,26,2 7]= 1 ，f[2 0,2 1,22,23 ,24,25 ,26 ,27,28 ]= 0 。
3.设，‖A‖∞＝ ___5 ____，‖X‖∞＝__ 3_____，‖AX ‖∞≤ 15_ __ 。
4.非线性方程 f x)=0的迭代函数 x x 在有解区间满足’x)| <1，则使用该迭(=( ) | (代函数的迭代解法一定是局部收敛的。
5.区间 [a,b]上的三次样条插值函数 S(x)在[a,b]上具有直到 2 阶的连续导数。
6.当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商公式的前插公式，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的后插公式；如果要估计结果的舍入误差，应该选用插值公式中的拉格朗日插值公式。
7.拉格朗日插值公式中f(xi)的系数 ai(x)的特点是：a i ( x)1；所以i当系数 ai(x)满足ai(x)>1 ，计算时不会放大f(xi)的误差。
8.要使 20 的近似值的相对误差小于0.1% ，至少要取 4 位有效数字。
9.对任意初始向量(0) 及任意向量g ，线性方程组的迭代公式x(k+1)=Bx(k)+ (=0,1, ?)X g k收敛于方程组的精确解x*的充分必要条件是(B)<1 。
10.由下列数据所确定的插值多项式的次数最高是 5 。
x 0 0.5 1 1.5 2 2.5y=f(x) -2 -1.75 -1 0.25 2 4.2511.牛顿下山法的下山条件为|f(xn+1)|<|f(xn)| 。
12. 线性方程组的松弛迭代法是通过逐渐减少残差ri ( i=0,1, ? ,n)来实现的，其中的残差r ＝ (b -a x-a x-? -a x )/aii， (i=0,1,,n)。
数值分析计算方法试题集及答案
数值分析复习试题第一章绪论一.填空题1.为精确值的近似值；为一元函数的近似值；*xx ()**x f y =()x f y =1为二元函数的近似值，请写出下面的公式：：()**,*y x f y =()y x f y ,2=**e x x =-***r x xe x -=()()()*'1**y f x x εε≈⋅()()()()'***1**r r x f x y x f x εε≈⋅()()()()()**,**,*2**f x y f x y y x y x yεεε∂∂≈⋅+⋅∂∂()()()()()****,***,**222r f x y e x f x y e y y x y y y ε∂∂≈⋅+⋅∂∂2、计算方法实际计算时，对数据只能取有限位表示，这时所产生的误差叫舍入误差。
3、分别用2.718281，2.718282作数e 的近似值，则其有效数字分别有 6 位和7（三位有效数字）。
1.73≈-211.73 10 2-≤⨯4、设均具有3位有效数字，则的相对误差限为 0.0055 。
121.216, 3.654x x ==12x x 5、设均具有3位有效数字，则的误差限为 0.01 。
121.216, 3.654x x ==12x x +6、已知近似值是由真值经四舍五入得到,则相对误差限为0.0000204 .2.4560A x =T x 7、递推公式如果取作计算,则计算到时,误差为,⎧⎪⎨⎪⎩0n n-1y =y =10y -1,n =1,2,0 1.41y =≈10y ;这个计算公式数值稳定不稳定 不稳定 .8110 2⨯8、精确值，则近似值和分别有 3 位和14159265.3*=π141.3*1=π1415.3*2=π4 位有效数字。
9、若,则x 有 6 位有效数字，其绝对误差限为1/2*10-5 。
*2.71828x e x =≈=10、 设x*的相对误差为2％，求(x*)n 的相对误差0.02n11、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字；12、计算方法主要研究( 截断 )误差和( 舍入 )误差；13、为了使计算 的乘除法次数尽量地少，应将该表达式()()2334610111y x x x =++----改写为11,))64(3(10-=-++=x t t t t y ，为了减少舍入误差，应将表达式19992001-改写为199920012+。
(完整)数值分析学期期末考试试题与答案(A),推荐文档
期末考试试卷（A 卷）2007学年第二学期 考试科目： 数值分析 考试时间：120 分钟学号 姓名 年级专业一、判断题（每小题2分，共10分）1. 用计算机求1000100011n n=∑时，应按照n 从小到大的顺序相加。
（ ）2. 为了减少误差,进行计算。
（ ）3. 用数值微分公式中求导数值时，步长越小计算就越精确。
（ ）4. 采用龙格－库塔法求解常微分方程的初值问题时，公式阶数越高，数值解越精确。
（ ）5. 用迭代法解线性方程组时，迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有关，与常数项无关。
（ ）二、填空题（每空2分，共36分）1. 已知数a 的有效数为0.01，则它的绝对误差限为________，相对误差限为_________.2. 设1010021,5,1301A x -⎡⎤⎡⎤⎢⎥⎢⎥=-=-⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦则1A =_____，2x =______，Ax ∞=_____.3. 已知53()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---= .4. 为使求积公式11231()((0)f x dx A f A f A f -≈++⎰的代数精度尽量高，应使1A = ，2A = ，3A = ，此时公式具有 次的代数精度。
5. n 阶方阵A 的谱半径()A ρ与它的任意一种范数A 的关系是 .6. 用迭代法解线性方程组AX B =时，使迭代公式(1)()(0,1,2,)k k XMX N k +=+=K 产生的向量序列{}()k X收敛的充分必要条件是 .7. 使用消元法解线性方程组AX B =时，系数矩阵A 可以分解为下三角矩阵L 和上三角矩阵U 的乘积，即.A LU = 若采用高斯消元法解AX B =，其中4221A -⎡⎤=⎢⎥⎣⎦，则L =_______________，U =______________；若使用克劳特消元法解AX B =，则11u =____；若使用平方根方法解AX B =，则11l 与11u 的大小关系为_____（选填：>，<，=，不一定）。
数值分析计算方法试题集及答案
数值分析复习试题第一章 绪论 一. 填空题 1.*x 为精确值x 的近似值;()**x f y =为一元函数()x f y =1的近似值；()**,*y x f y =为二元函数()y x f y ,2=的近似值，请写出下面的公式：**e x x =-:***r x xe x -=()()()*'1**y f x x εε≈⋅ ()()()()'***1**r r x f x y x f x εε≈⋅()()()()()**,**,*2**f x y f x y y x y x yεεε∂∂≈⋅+⋅∂∂()()()()()****,***,**222r f x y e x f x y e y y x y y y ε∂∂≈⋅+⋅∂∂ 2、 计算方法实际计算时，对数据只能取有限位表示，这时所产生的误差叫 舍入误差.3、 分别用2.718281，2.718282作数e 的近似值，则其有效数字分别有 6位和 7 1.73≈（三位有效数字）-211.73 10 2≤⨯。
4、 设121.216, 3.654x x ==均具有3位有效数字,则12x x 的相对误差限为 0.0055 .5、 设121.216, 3.654x x ==均具有3位有效数字，则12x x +的误差限为 0。
01 。
6、 已知近似值 2.4560A x=是由真值T x 经四舍五入得到,则相对误差限为 0。
0000204 。
7、 递推公式,⎧⎪⎨⎪⎩0n n-1y =y =10y -1,n =1,2,如果取01.41y ≈作计算，则计算到10y 时，误差为8110 2⨯;这个计算公式数值稳定不稳定 不稳定 。
8、 精确值 14159265.3*=π，则近似值141.3*1=π和1415.3*2=π分别有 3 位和 4 位有效数字。
9、 若*2.71828x e x =≈=,则x 有 6 位有效数字，其绝对误差限为1/2＊10—5。
数值分析试卷及答案
数值分析试卷及答案**注意：以下是一份数值分析试卷及答案，试卷和答案分别按照题目和解答的格式排版，以确保整洁美观，语句通顺。
**---数值分析试卷一、选择题（每题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：从第一行开始，选取第一个非零元素作为主元，然后通过行变换将其它行的该列元素消去。
数值分析计算方法期末考试(一) +答案
4. 计算积分 1 x dx , 取 4 位有效数字, 用梯形公式计算求得的值为_______,用辛普森 0.5 公式计算求得的值为_____________.
5. 设 f (x) 可微, 求方程 x f (x) 根的牛顿迭代格式是______________________.
=
2
x
3 k
+1 ,
(k
3xk2 − 3
= 0,
1,
2,) ,
计算得
x0 = 2, x1 = 1.8889, x2 = 1.8795, 由 | x2 − x* |≈ 1.148 ×10−4 < 5 ×10−5 ,
所以 x ≈ 1.8795
(2) 用弦截法:
xk +1
= xk
−
f
(
xk
f )
(xk −f
1 x2 式,并用它分别计算 x 1 处的值.
3 2. 求三个不同的节点 x0 , x1, x2 和常数 c ,使求积公式
1
1 f (x) dx c[ f (x0 ) f (x1) f (x2 ) ]
具有尽可能高的代数精确度.
3. 用最小二乘法求下列数据的线性拟合函数 y ax b
xi -2
Lagrange 插值多项式为
f (−1) = 1 , f (0) = 1, f (1) = 1 . 由 此 确 定 的 二 次
2
2
L2 (x)
=
x(x −1) (−1)(−1 −1)
f
(−1) +
(x + 1)(x −1) (0 + 1)(0 −1)
数值分析报告期末考试复习题及其问题详解
数值分析期末考试复习题及其答案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 也是复数，上述结论依旧成立。
(完整)数值分析学期期末考试试题与答案(A),推荐文档
期末考试试卷( A 卷)2007 学年第二学期 考试科目： 数值分析 考试时间： 120 分钟学号 姓名 年级专业100011. 用计算机求11000时，应按照 n 从小到大的顺序相加。
n1n2. 为了减少误差 ,应将表达式 2001 1999 改写为 2进行计算。
( )2001 19993. 用数值微分公式中求导数值时，步长越小计算就越精确。
( )4. 采用龙格－库塔法求解常微分方程的初值问题时， 公式阶数越高，数值解越精确。
( )5. 用迭代法解线性方程组时， 迭代能否收敛与初始向量的选择、 系数矩阵及其演变方式有关，与常数项无关。
( ) 二、填空每空 2 分，共 36 分)1. 已知数 a 的有效数为 0.01 ，则它的绝对误差限为 _______ ，相对误差限为 _1 0 1 02. 设 A0 2 1 ,x 5 ,则 A 1____________________________ _， x 2 ______ ，Ax1 3 0 13. 已知 f (x) 2x 54x 35x,则 f[ 1,1,0] , f[ 3, 2, 1,1,2,3] .14. 为使求积公式 f (x)dx A 1f ( 3) A 2f (0) A 3f ( 3)的代数精度尽量高，应使13 3A 1 ， A 2 ， A 3，此时公式具有 次的代数精度。
5. n 阶方阵 A 的谱半径 ( A)与它的任意一种范数 A 的关系是 .6. 用迭代法解线性方程组 AX B 时，使迭代公式 X (k 1)MX (k)N (k 0,1,2,K )产 生的向量序列X (k)收敛的充分必要条件是 .7. 使用消元法解线性方程组AX B时，系数矩阵A可以分解为下三角矩阵L 和上三角矩阵U 的乘积，即A LU. 若采用高斯消元法解AX B，其中A 4 2，则21L ___________ ，U ____________ ；若使用克劳特消元法解AX B ，则u11 _______ ；若使用平方根方法解AX B，则l11与u11的大小关系为（选填：＞，＜，=，不一定）。
研究生数值分析期末考试试卷参考答案
研究生数值分析期末考试试卷参考答案太原科技大学硕士研究生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分。
数值分析计算方法试题集及答案
第一章绪论 一. 填空题 x* -x4、 设x j —1.216, x^ - 3.654均具有3位有效数字，则 X 1X 2的相对误差限为 0.0055 ________ 。
5、 设为=1.216,x 2 =3.654均具有3位有效数字，则 为 他 的误差限为 0.01 _______________ 。
6、 已知近似值X A = 2.4560是由真值x T 经四舍五入得到，则相对误差限为0.0000204 .7、 递推公式二血 ,如果取y 0=逅".41作计算,则计算到 血时,误差为y n = 10y n-i -1, n = 1,2,11 8-10 ;这个计算公式数值稳定不稳定 不稳定 .28、 精确值二* =3.14159265…，则近似值=3.141和二2* =3.1415分别有_J3 _____________ 位和4 _ 位有效数字。
9、 若x 二e 、2.71828二x ,则x 有_6_位有效数字，其绝对误差限为 1/2*10 -。
10、 设x*的相对误差为2%,求（x*） n 的相对误差0.02n数值分析复习试题1. x 为精确值x 的近似值；y * = f x *为一元函数y1二 f X 的近似值；二f x*, y*为二元函数y2f x, y 的近似值，请写出下面的公式:*e 二 x* -x :f x*x*；r yi* :x* f f x*x*；r x*y2*f x*, y*:x（X*）+ 硏（x*,y*）:y舍入误差 _______ 。
6 _____ 位和 7位; 又取 乔“.73 （三位有效数字），则1.73＜丄汉10-22计算方法实际计算时，对数据只能取有限位表示，这时所产生的误差叫 分别用2.718281，2.718282作数e 的近似值，则其有效数字分别有ry2*：11、近似值x'0.231关于真值x =0.229有(2 )位有效数字;2..2001 • ”1999 。
重庆大学《数值分析》期末考试真题及答案
重庆大学《数值分析》期末考试真题及答案(总8页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--一.填空题: 1． 若求积公式对任意不超过 m 次的多项式精确成立,而对 m+1 次多项式不成立,则称此公式的代数精度为m 次. 2． 高斯消元法求解线性方程组的的过程中若主元素为零会发生 计算中断 ;.主元素的绝对值太小会发生 误差增大 .3． 当A 具有对角线优势且 不可约 时，线性方程组Ax=b 用简单迭代法和塞德尔迭代法均收敛.4． 求解常微分方程初值问题的欧拉方法是 1 阶格式; 标准龙格库塔法是 4 阶格式.5． 一个n 阶牛顿-柯特斯公式至少有 n 次代数精度,当n 偶数时,此公式可以有n+1 次代数精度.6． 相近数 相减会扩大相对误差,有效数字越多,相对误差 越大 .二计算题: 1. 线性方程组:⎪⎩⎪⎨⎧-=++-=+-=++5.1526235.333321321321x x x x x x x x x 1) 对系数阵作LU 分解,写出L 阵和U 阵；⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎭⎫⎝⎛--=79/123/5413314/33/113/11U L 2) 求出此方程组的解.)5.0,1,2('-=x2. 线性方程组:⎪⎩⎪⎨⎧=++-=++=++332212325223321321321x x x x x x x x x 1)对系数阵作LU 分解,写出L 阵和U 阵；⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫⎝⎛=573235223152321321//////U L2)求出此方程组的解.),,('-=133x3) 此方程组能否用用简单迭代法和高斯塞德尔迭代法求解.0732223222305322303>=>=>,,A 对称正定,用高斯-塞德尔迭代法收敛;..,.,//////)(,6667033331027163432323232323232131=-==+-=-⎪⎪⎪⎭⎫ ⎝⎛-=+-=-λλλλλJ J B I U L D B 用简单迭代法不收敛3. 设f (x )= x 4, 以-1,0,1,2为插值节点,1) 试写出f (x )的三次拉格朗日插值多项式P 3(x )及其插值余项R 3(x );6)2)(1())()(())()(()(3020103210---=------=x x x x x x x x x x x x x x x x l 2)2)(1)(1())()(())()(()(3121013201--+=------=x x x x x x x x x x x x x x x x l 2)2)(1())()(())()(()(3212023102-+-=------=x x x x x x x x x x x x x x x x l 6)1)(1())()(())()(()(2313032103-+=------=x x x x x x x x x x x x x x x x l )(8)()()(3203x l x l x l x P ++=())2)(1)(1()2)(1()1(!4)()4(43--+=--+=x x x x x x x x x x R 2) 求出f 的近似值,并估计误差.0625.55.1)5.1(4==f-0.93755.05.05.25.1)2)(1)(1()5.1(3=-⨯⨯⨯=--+=x x x x R 6)9375.0(0625.5)5.1(3=--=P或:0.3125610.9375 0625.0)5.1(8)5.1()5.1()5.1(3203⨯++=++=l l l P =6 -0.937560625.5)5.1()5.1()5.1(33=-=-=P f R4 设x x f ln )(=, 以1,2,3为插值节点,1) 试写出f (x )的二次拉格朗日插值多项式P 2(x )及其插值余项R 2(x );2322010210))(())(())(()(--=----=x x x x x x x x x x x l ))(())(())(()(312101201---=----=x x x x x x x x x x x l2211202102))(())(())(()(--=----=x x x x x x x x x x x l98080124711438009861693102212...)(.)(.)(-+-=+=x xx l x l x P 23112312333ln ()()()()()()()!R x x x x x x x ξξ'''=---=---2) 求出)(ln e p e 2≈的近似值,与精确值1比较,并用误差公式估计误差限.0135010135122.,ln ,.)(===R e e p231123123331171830718302817011593ln ()()()()()()()!..(.).R e e e e e e e ξξ'''=---=---≤⨯⨯⨯-=5 有积分公式()()2)0(2)(33f c f b f a dx x f ⨯+⨯+-⨯=⎰-，c b a ,,是待定参数，试确定c b a ,,，使得上述公式有尽可能高的代数精度,并确定代数精度为多少.⎰⎰⎰---==+==+-==++==332333318)(40)(2612,1,0,)(dx x b a xdx b a dx c b a k x x f k)]()()([)(/,/33023343234933f f f dx x f c b a ++-====∴⎰- 至少有2次代数精度.[][]10872072435486,024024430,)(33433343=++≠==++-===⎰⎰--dx x dx x x x x f此公式代数精度为3. 6 有积分公式)]2(3)0(2)2(3[43)(33f f f dx x f ++-=⎰- 1) 试确定代数精度为多少;2) 用它计算⎰-33dx e x,精确到2位小数,与3333---=⎰e e dx e x 作比较.[][][][][]10872072435486,02402443012012431860643032343614,3,2,1,0,)(3343333323333=++≠==++-==++==++-==++====⎰⎰⎰⎰⎰-----dx x dx x dx x xdx dx k x x f k代数精度为3.04.2043.18]323[43333320332=-==++≈⎰⎰----e e dx e e e e dx e x x7. 某企业产值与供电负荷增长情况如下表：1) 试用一次多项式拟合出经验公式bx a y +=；⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛415521186062625..b a 解之: 0006101811.,.ab =-=0061018110..-=x y8. 测试某型号水泵得到扬程(米)和出水量(立米/小时)的对照表如下:1)试用一次多项式拟合出经验公式x ba y +=；bX a y x X +==,/1⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛12365491404514515....b a 解之: 59953043864.,.ba ==-38644953059./.-=x y2) 计算拟合值填入上表的空格，看是否与实际值基本吻合； 3) 某用户使用此型号水泵时扬程为米,试估计此时出水量67183864462953059.../.=-=y9 方程01=-+-x xe x有一个实根：1)用区间对分法搜索确定根所在的区间 (a,b )，使 b-a ≤;,1) 用某种迭代法求出此正根,精确到5位有效数字65905.0*≈x10 方程x e x-= 1) 证明它在(0,1)区间有且只有一个实根； 2) 证明 ,,,101==-+k e x k x k ，在(0,1)区间内收敛；3) 用牛顿迭代法求出此根,精确到5位有效数字1),.)(,)(,)(063201100>=-==-=-f f e x x f x(0,1)区间有一个实根；)(,)(x f e x f x 011>>+='-是严格增函数，只有一个实根。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
代入求解得：
5.解，设A可以三解分解，即
由矩阵的乘法及矩阵相等可得：
，
令
求解三角方程组： ，得：
求解三角方程组： ，得：
所以，原方程组的解为：
三．证明
证明：分别将 ， ， 在 处用Taylor公式展开得：
将以上三式代入线性二步法中，得：
又方程的真解的Taylor展式为：
所以，，
2.解：设 则
牛顿迭代公式为：
将 代入上式，得 ， ，
所以，方程的近似根为：
3.解：设 时，左 ，右 ，左＝右得：
时，左 ，右 ，左＝右得：
时，左 ，右 ，左＝右得：
时，左 ，右 ，左＝右得：
联立上述四个方程，解得：
时，左 ，右 ，左 右
所以，该求积公式的代数精度是3
4.解：Euler公式是：
标准答案
一．填空
1.舍入误差2. 729，1，0 3. 5
4. 5. 6，3， ,9
二．计算
1.构造重节点的差商表：
n
x
y
一阶
二阶
0
1
2
1
1
2
0
2
2
3
1
1
所以，要求的Hermite插值为：
2.
证明：由题意可知
由插值条件知：
所以，可设： （＃）
构造函数：
易知： 时， ，且 至少有一个根 ，即
对（＃）式求三阶导，并代入得：