(整理)数值分析第二学期期末考试试题与答案(a)
- 格式:doc
- 大小:329.50 KB
- 文档页数:13
《数学分析下册》期末考试卷及参考答案一、填空题(第1题每空2分,第2,3,4,5题每题5分,共26分)1、已知uln某2y2,则uu,,y某du2、设L:某2y2a2,则某dyyd某L某=3cot,L:3、设(0t2),则曲线积分(某2+y2)d=y=3int.L4、改变累次积分dy(f某,y)d某的次序为2y33某y1,则(51)d某dy=5、设D:D得分阅卷人二、判断题(正确的打“O”;错误的打“某”;每题3分,共15分)p某0,y0)p某0,y0)1、若函数(在点(连续,则函数(点(必存在一f某,y)f某,y)阶偏导数。
()p某0,y0)p某0,y0)2、若函数(在点(可微,则函数(在点(连续。
f某,y)f某,y)()p某0,y0)3、若函数(在点(存在二阶偏导数f某y(某0,y0)和fy某(某0,y0),则f某,y)必有f某y(某0,y0)fy某(0某,0y) L(B,A)()()4、L(A,B)f(某,y)d某f(某,y)d某。
5、若函数(在有界闭区域D上连续,则函数(在D上可积。
()f某,y)f某,y)第1页共5页得分阅卷人三、计算题(每小题9分,共45分)1、用格林公式计算曲线积分I(e某iny3y)d某(e某coy3)dy,AOAO为由A(a,0)到O(0,0)经过圆某2y2a某上半部分的路线。
其中2、计算三重积分------线--------------------------------------(某V2y2)d某dydz,其中是由抛物面z某2y2与平面z4围成的立体。
第2页共5页3、计算第一型曲面积分IdS,S其中S是球面某2y2z2R2上被平面za(0aR)所截下的顶部(za)。
4、计算第二型曲面积分22Iy(某z)dydz某dzd某(y某z)d某dy,S其中S是立方体V0,b0,b0,b的外表面。
第3页共5页5、设D(某,y)某2y2R曲顶柱体的体积。
得分阅卷人四、证明题(每小题7分,共14分)1、验证曲线积分第4页共5页2.求以圆域D为底,以曲面ze(某2y2)为顶的(某22yz)d某(2y2某)zdy2(z2,某)ydzL与路线无关,并求被积表达式的一个原函数u(某,y,z)。
期末考试试卷(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 X MX N k +=+=产生的向量序列{}()k X 收敛的充分必要条件是 .7. 使用消元法解线性方程组AX B =时,系数矩阵A 可以分解为下三角矩阵L 和上三角矩阵U 的乘积,即.A LU = 若采用高斯消元法解AX B =,其中4221A -⎡⎤=⎢⎥⎣⎦,则L =_______________,U =______________;若使用克劳特消元法解AX B =,则11u =____;若使用平方根方法解AX B =,则11l 与11u 的大小关系为_____(选填:>,<,=,不一定)。
一、单项选择题(每小题3分,共15分) 1、用Simpson 公式求积分1401x dx +⎰的近似值为 ( ).A.2924 B.2429C.65D. 562、已知(1)0.401f =,且用梯形公式计算积分2()f x dx ⎰的近似值10.864T =,若将区间[0,2]二等分,则用递推公式计算近似值2T 等于( ). A.0.824 B.0.401 C.0.864 D. 0.8333、设3()32=+f x x ,则差商0123[,,,]f x x x x 等于( ).A.0B.9C.3D. 64的近似值的绝对误差小于0.01%,要取多少位有效数字( ). A.3 B.4 C.5 D. 25、用二分法求方程()0=f x 在区间[1,2]上的一个实根,若要求准确到小数 点后第四位,则至少二分区间多少次( ).A.12B.13C.14D. 15二、填空题(每小题4分,共40分)1、对于迭代函数2()=(3)ϕ+-x x a x ,要使迭代公式1=()ϕ+k k x x则a 的取值范围为 .2、假设按四舍五入的近似值为2.312,则该近似值的绝对误差限为 .3、迭代公式212(3)=,03++>+k k k k x x a x a x a收敛于α= (0)α>. 4、解方程4()530f x x x =+-=的牛顿迭代公式为 . 5、设()f x 在[1,1]-上具有2阶连续导数,[1,1]x ∀∈-,有1()2f x ''≤,则()f x 在[1,1]-上的线性插值函数1()L x 在点0处的误差限1(0)R ≤______.6、求解微分方程初值问题2(0)1'=-⎧⎨=⎩y xy yy ,0x 1≤≤的向前Euler 格式为 .7、设310131013A -⎛⎫⎪=-- ⎪ ⎪-⎝⎭,则A ∞= .8、用梯形公式计算积分112-⎰dx x 的近似值为 . 9、设12A 21+⎡⎤=⎢⎥⎣⎦a 可作Cholesky 分解,则a 的取值范围为 . 10、设(0)1,(0.5) 1.5,(1)2,(1.5) 2.5,(2) 3.4f f f f f =====,若1=h ,则用三点公式计算(1)'≈f .三、解答题(共45分) 1、给定数据用复化Simpson 公式计算1.381.30()f x dx ⎰的近似值,并估计误差,小数点后保留3位. (8分)2、用直接三角分解法求线性代数方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡432631531321321x x x 的解. (8分) 3、求()λx ,使得迭代公式1()()λ+=+k k k k f x x x x 求方程2()31=+-f x x x 的根的相应迭代序列{}k x 具有平方收敛. (5分)4、已知数据试对数据用最小二乘法求出形如=+y x b的拟合曲线. (8分) 5、已知(2)8f -=,(0)4f =,(2)8=f ,试求二次拉格朗日插值多项式. (8分) 6、设矩阵A 如下,根据谱半径判断用Jacobi 迭代法求解方程组Ax b =的敛散性.(8分)1102111221012A ⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦一、单项选择题(每小题3分,合计15分) 1、A 2、D 3、C 4、C 5、D 二、填空题(每小题3分,合计30分) 1、0<<a ; 2、31102-⨯; 3;4、4135345++-=-+k k k k k x x x x x ; 5、14; 6、1(2)+=+-n n n n n y y h x y y ; 7、5;8、34-; 9、3>a ;10、1.2;三、计算题(合计55分) 1、给定数据用复化Simpson 公式计算 1.381.30()f x dx ⎰的近似值,并估计误差,小数点后保留3位. (8分)解: 401024S [()4()()]6-=++x x f x f x f x ………… 1分 1.38 1.30(3.624 4.20 5.19)6-=+⨯+ 0.341= ………… 2分20422012234S [()4()()][()4()()]66--=+++++x x x xf x f x f x f x f x f x =0.342 ………… 6分2211[]15-≈-I S S S =-⨯40.6710 ………… 8分 2、用直接三角分解法求线性代数方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡432631531321321x x x 的解. (8分) 解:设111213212223313233u u u 123100135l 100u u 136l l 100u ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=*⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦………… 1分 111=u ,212=u ,313=u ,121=l ,131=l 122=u ,223=u ,132=l133=u ,133=l …………6分所以⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111011001L ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100210321U …………7分 由b Ly =得Ty )1,1,2(=;由y Ux =得Tx )1,1,1(-=. ………… 8分3、求()λx ,使得迭代公式1()()λ+=+k k k k f x x x x 求方程2()31=+-f x x x 的根的相应迭代序列{}k x 具有平方收敛.(6分)解:要使迭代序列具有平方收敛,则()0ϕ'*=x ………… 2分 而()()()ϕλ=+f x x x x ,即 ………… 3分 2()()()()10()λλλ''**-**+=*f x x x f x x …………4分 而()0*=f x 则有()1()λ'*=-*f x x ………… 5分所以()()23λ'=-=--x f x x ………… 6分4、已知数据试对数据用最小二乘法求出形如=+ay x b的拟合曲线. (8分) 解:因为11=+b x y a a ,令0111,,,====b a a y x x a a y……2分 则有法方程01461061410⎛⎫⎛⎫⎛⎫= ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭a a ……5分解出014,1==-a a ,则1,4=-=-a b ……7分 所以1=4-y x……8分5、已知(2)8f -=,(0)4f =,(2)8=f ,试求二次拉格朗日插值多项式. (7分)解:01()(2)8l x x x =- …………2分 211()(4)4l x x =-- …………4分21()(2)8l x x x =+ …………6分 2012()()(2)()(0)()(2)L x l x f l x f l x f =-++24=+x …………7分6、设矩阵A 如下,根据谱半径判断用Jacobi 迭代法求解方程组Ax b =的敛散性.(8分)1102111221012A ⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦解:100010001D ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,00010021002L ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦,10021002000U ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦…………3分1100211()0221002J B D L U -⎡⎤⎢⎥⎢⎥⎢⎥=+=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦…………5分 2102111()0222102J E B λλλλλλ⎡⎤-⎢⎥⎢⎥⎢⎥-=--=-=⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦…………6分()2J B ρ=…………7分 所以用Jacobi 迭代法求解方程组Ax b =收敛 …………8分。
数值分析期末试题及答案试题一: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)。
卷)期末考试试卷(A2007学年第二学期考试科目:数值分析分钟考试时间:120年级专业学号姓名题号一2二三0四总分分)分,共10一、判断题(每小题210001?n)( 1. 用计算机求时,应按照从小到大的顺序相加。
1000n1n?219992001?为了减少误差2. ,应将表达式进行计算。
(改写为)19992001?)( 3. 用数值微分公式中求导数值时,步长越小计算就越精确。
) 4. 采用龙格-库塔法求解常微分方程的初值问题时,公式阶数越高,数值解越精确。
(系数矩阵及其演变方式有用迭代法解线性方程组时,5. 迭代能否收敛与初始向量的选择、)(关,与常数项无关。
分)二、填空题(每空2分,共36_________.________,相对误差限为已知数a的有效数为0.01,则它的绝对误差限为1.0?110??????????xA?Ax,0?21,x??5A?_____.则设______,_____,2. ????21?????1?130????53f(x)?2x?4x?5x,f[?1,1,0]?f[?3,?2,?1,1,2,3]? 3. 已知则, .331?)?Af(0)?Af(f(x)dx?Af(?)的代数精度尽量高,应使4. 为使求积公式321331?A?A?A?,此时公式具有,,次的代数精度。
312?nA)(A的关系是 5. A阶方阵的谱半径与它的任意一种范数.(k?1)(k)BAX??N(k?XMX?0,1,2,)产时,使迭代公式用迭代法解线性方程组6.??)k(X .生的向量序列收敛的充分必要条件是AX?BAL和上三角矩7. 使用消元法解线性方程组系数矩阵时,可以分解为下三角矩阵14?2??BAX?.A?LUU?A,则阵若采用高斯消元法解的乘积,即,其中??21??L?U?AX?B,则,______________;若使用克劳特消元法解_______________u?lu BAX?的大小关系为_____(选填:则____;若使用平方根方法解>与,,111111<,=,不一定)。
数值分析期末考试题及答案一、选择题(每题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 \) 处插值,并求出插值多项式。
华南农业大学期末考试试卷(A 卷)2006学年第二学期 考试科目: 数值分析 考试时间:120 分钟学号 姓名 年级专业一、判断题(每小题2分,共10分)1. 用计算机求1000511n n=∑时,应按照n 从大到小的顺序相加。
( )2. 插值节点越多,插值函数与原始函数越接近。
( )3. 采用复化的求积分公式计算()baf x dx ⎰的近似值时,如果节点取得过多,误差反而会增大。
( )4. 最小二乘法可以求超定方程组的精确解。
( )5. 采用龙格-库塔法求解常微分方程的数值解时,公式的阶数越高,数值解越精确。
( )二、填空题(每空2分,共36分)1. 已知数a 的有效数为0.025,则它与a 的相对误差限是___________.2. 测量一个固定高度的正圆锥形容器的容积时,要使测量值的相对误差限为0.01,则测量底面半径时允许的相对误差限为 ___________.3. 设3210022,4,1201A x ⎡⎤⎡⎤⎢⎥⎢⎥=-=⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦则1A =______,2x =______,Ax ∞=______.4. 已知642()2651,f x x x x =+-+则[1,1,0]f -= ,[3,2,1,0,1,2,3]f ---= .5. 为使求积公式11231()(1)(0)(1)f x dx A f A f A f -≈-++⎰的代数精度尽量高,应使1A = ,2A = ,3A = .6. 求定积分()baf x dx ⎰的梯形公式为 ,辛普生公式为 .7. 用于解方程()0f x =的双点弦截法的迭代公式为_____________________________. 8. 利用克劳特消元法解方程组AX=B 时,对系数矩阵A 做LU 分解时11u =____; 平方根方法中11l 与11u 的大小关系为_____(选填:>,<,=,不一定)。
9. 对矩阵4221A -⎡⎤=⎢⎥⎣⎦作LU 分解,使得A=LU. 如果选择11222,l l ==则L=_________,U= __________.10. 以步长为1的二阶泰勒级数法求解初值问题2(0)1y x y y '⎧=+⎨=⎩的数值解,其迭代公式为___________________________.三、计算题(每小题9分,共54分)1. (1)证明用牛顿迭代法求解方程2sin 0x x -=的正根是收敛的;(2)求出这个根(只需给出迭代前两步的值)。
一. 单项选择题(每小题2分,共10分)1. 在下列四个数中,有一个数具有4位有效数字,且其绝对误差限为 51021-⨯,则该数是( ) A 0.001523 B 0.15230 C 0.01523 D 1.52300 2. 设方阵A 可逆,且其n 个特征值满足:n λλλ>≥> (21),则1-A 的主特征值是( )A11λ B nλ1 C1λ或n λ D 11λ或nλ13. 设有迭代公式→→+→+=fxB x k k )()1(。
若||B|| > 1,则该迭代公式( )A 必收敛B 必发散C 可能收敛也可能发散4. 常微分方程的数值方法,求出的结果是( )A 解函数B 近似解函数C 解函数值D 近似解函数值 5. 反幂法中构造向量序列时,要用到解线性方程组的( ) A 追赶法 B LU 分解法C 雅可比迭代法D 高斯—塞德尔迭代法二. 填空题(每小题4分,共20分)1. 设有方程组⎪⎩⎪⎨⎧=+-=+-=+02132432132132x x x x x x x x ,则可构造高斯—塞德尔迭代公式为⎪⎩⎪⎨⎧2. 设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=111112101A ,则=∞A3. 设1)0(,2'2=+=y y x y ,则相应的显尤拉公式为=+1n y4. 设1)(+=ax x f ,2)(x x g =。
若要使)(x f 与)(x g 在[0,1]上正交,则a =5. 设T x )1,2,2(--=→,若有平面旋转阵P ,使P →x 的第3个分量为0,则P =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡ 三. 计算题(每小题10分,共50分)1. 求27的近似值。
若要求相对误差小于0.1%,问近似值应取几位有效数字?2. 设42)(x x x f -=,若在[-1,0]上构造其二次最佳均方逼近多项式,请写出相应的法方程。
3. 设有方程组⎪⎩⎪⎨⎧=++=++=-+1221122321321321x x x x x x x x x ,考察用雅可比迭代解此方程组的收敛性。
期末考试试卷(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)33f 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 X收敛的充分必要条件是 .7. 使用消元法解线性方程组AX B =时,系数矩阵A 可以分解为下三角矩阵L 和上三角矩阵U 的乘积,即.A LU = 若采用高斯消元法解AX B =,其中4221A -⎡⎤=⎢⎥⎣⎦,则L =_______________,U =______________;若使用克劳特消元法解AX B =,则11u =____;若使用平方根方法解AX B =,则11l 与11u 的大小关系为_____(选填:>,<,=,不一定)。
8. 以步长为1的二阶泰勒级数法求解初值问题(0)1y x yy '=+⎧⎨=⎩的数值解,其迭代公式为___________________________.三、计算题(第1~3、6小题每题8分,第4、5小题每题7分,共46分)1. 以02x =为初值用牛顿迭代法求方程3()310f x x x =--=在区间(1,2)内的根,要求(1) 证明用牛顿法解此方程是收敛的;(2) 给出用牛顿法解此方程的迭代公式,并求出这个根(只需计算12,,x x 计算结果取到小数点后4位)。
2. 给定线性方程组1231231230.40.410.40.820.40.83x x x x x x x x x ++=⎧⎪++=⎨⎪++=⎩(1) 分别写出用Jacobi 和Gauss-Seidel 迭代法求解上述方程组的迭代公式;(2) 试分析以上两种迭代方法的敛散性。
3. 已知函数()y f x =在如下节点处的函数值(1) (2) 根据后三个节点建立二阶牛顿后插公式2()P x ,并计算(1.1)y 的近似值; (3) 采用事后估计法计算(2)中近似值的截断误差(结果保留四位小数)。
4.5. 已知函数()y f x =在以下节点处的函数值,利用差商表求(3)f '和(3)f ''的近似值。
6. 写出前进欧拉公式、后退欧拉公式,并由这两个公式构造一个预估-校正公式求解下列常微分方程的数值解。
22(01,0.2)(0)0y x y x h y '⎧=+≤≤=⎨=⎩四、(8分)已知n+1个数据点(,)(0,1,2,,)i i x y i n ,请用多种方法建立这些数据点之间的函数关系,并说明各种函数的适用条件。
期末考试答案及评分标准(A 卷)2007学年第二学期 考试科目: 数值分析一、判断题:(每小题2分,共10分)1. ×2. √3. ×4. ×5. ×二、填空题:(每空2分,共36分) 1. 0.005或20.510-⨯ ,0.5 2.3. 0,24. 1,0,1,35.()A A ρ≤6. ()1M ρ<7. 1042,,1,10212⎡⎤-⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦⎣⎦8. 11()(1)2n n n n n n y y x y x y +=+++++或1 1.5 2.50.5,0,1,2,n n n y x y n +=++=三、解答题(第1~4小题每题8分,第5、6小题每题7分,共46分) 1. (1)证明:3()31f x x x =--,由于a) (1)30,(2)10,f f =-<=> b) 2()330((1,2)),f x x x '=-≠∈c)()60((1,2)),f x x x ''=>∈ 即()f x ''在(1,2)上不变号,d) 对于初值02x =,满足(2)(2)0,f f ''> 所以用牛顿迭代法求解此方程是收敛的。
………………………………………4分(2)解:牛顿迭代法的迭代公式为312()31()33n n n n n n n n f x x x x x x f x x +--=-=-'- ………………………………………2分取初值02x =进行迭代,得1 1.8889,x =………………………………………1分2 1.8795.x =………………………………………1分2. 解:(1)Jacobi 迭代公式为(1)()()123(1)()()213(1)()()3120.40.410.40.820.40.83k k k k k k k k k x x x x x x x x x +++⎧=--+⎪=--+⎨⎪=--+⎩ ……………………………2分 Gauss-Seidel 迭代公式为(1)()()123(1)(1)()213(1)(1)(1)3120.40.410.40.820.40.83k k k k k k k k k x x x x x x x x x ++++++⎧=--+⎪=--+⎨⎪=--+⎩……………………………2分 (2)Jacobi 迭代矩阵的特征方程为0.40.40.40.800.40.8λλλ=,展开得30.960.2560λλ-+=,即(0.8)(0.40.40λλλ-+++-=,从而得 123-1.0928,0.8000,0.2928λλλ===,(或由单调性易判断必有一个大于1的特征根,)因此迭代矩阵的谱半径等于必大于1,所以Jacobi 迭代法发散。
……………………………2分Gauss-Seidel 迭代矩阵的特征方程为0.40.40.40.800.40.8λλλλλλ=,展开得2(0.8320.128)0λλλ-+=,解得1230,0.628,0.204,λλλ=≈≈迭代矩阵的谱半径小于1,所以Gauss-Seidel 迭代法收敛。
……………………………2分3. 解:(1)建立差分表………………………………………2分 (2)建立牛顿后插公式为2232022********()()()()!!()()()P x x x x x x x x =-----=-----=-+ 则所求近似值为211279(.).P =………………………………………3分(3)根据前三个节点建立牛顿后插公式为1221431112312124()()()()!!()()P x x x x x x x x x =----=----=-++ 则 1211268()(.).P =根据事后误差估计法1222209091()()(.)(.)x R x P P x -⎡⎤≈-⎣⎦+ 故截断误差209112792680047121.(.)(..)..R -≈⨯-≈- ………………………………………3分4.5. 解:设所求二次最小平方逼近多项式为22012().P x a a x a x =++ 根据已知数据,得01211111002,,11151240a M A a Y a -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦……………………………2分则4268268,468186M M M Y ⎡⎤⎡⎤⎢⎥⎢⎥''==⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦……………………………1分建立法方程组为0124268268468186a a a ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦……………………………2分解得0123.5, 1.5, 1.5.a a a ===-……………………………1分从而得所求一次最小平方逼近多项式为21() 3.5 1.5 1.5.P x x x =+-……………………………1分6. 解:设2()P x 为已知节点数据的插值二次多项式。
构造如下差商表:……………………………2分因为二次多项式的二阶差商为常数,又2()P x 是()f x 的插值函数,故有225[4,3,3][3,3,3]2P P ==而22[3,3]75[4,3,3]342P P -==-,因此得29[3,3]2P =, ……………………………1分由于1()()![,,,,]k n k f x k P x x x x +≈,从而得293332()[,],f P '==2323335()![,,].f P ''==……………………………2分7. 解:前进欧拉公式:221(,)0.20.2n n n n n n n y y h f x y y x y +=+⋅=++…………1分后退欧拉公式:2211111(,)0.20.2n n n n n n n y y h f x y y x y +++++=+⋅=++ ……1分预估时采用欧拉公式*2210.20.2n n n n y y x y +=++……………………………1分校正时采用后退欧拉公式()22*1110.20.2n n n n y y xy+++=++……………………………1分由初值000002,,.x y h ===知,节点分别为0.2,(1,2,3,4,5)i x i i ==当10.2,x =*2210000.20.20,y y x y =++=()22101102020008*...y y x y=++=,当20.4,x =*2221110.20.20.0160,y y x y =++≈()222122020200401*...y y x y =++≈.……………………………1分当30.6,x =*2232220.20.20.0724,y y x y =++≈()223233020201131*...y y x y =++≈.……………………………1分当40.8,x =*2243330.20.20.1877,y y x y =++≈()224344020202481*...y y x y=++≈.……………………………1分当51.0,x =*2254440.20.20.3884,y y x y =++≈()225455020204783*...y y x y=++≈.四、(8分)答:1、可以建立插值函数: (1)Newton 基本差商公式00100121001110()()()[,]()()[,,]()()()[,,,]n n n P x f x x x f x x x x x x f x x x x x x x x x f x x x -=+-+--++---……………………………1分(2)Lagrange 插值多项式0011()()()()()n i i n n L x a f x a f x a f x a f x =+++++其中01101101()()()(),(,,,)()()()()i i n i i i i i i i n x x x x x x x x a i n x x x x x x x x -+-+----==----.……………………………1分这两类插值函数的适用条件是:n 不太大;而且要求函数严格通过已知数据点。