数值分析试卷及答案
数值分析试题及答案
数值分析试题及答案一、选择题1. 下列哪个方法不适合用于求解非线性方程的根?A. 二分法B. 牛顿法C. 弦截法D. 正割法2. 当使用二分法求解非线性方程的根时,需要满足的条件是:A. 函数f(x)在区间[a, b]上连续B. 函数f(x)在区间[a, b]上单调递增C. 函数f(x)在区间[a, b]上存在根D. 函数f(x)在区间[a, b]上可导3. 数值积分是通过将定积分转化为求和的方法来近似计算积分值的过程。
下列哪个方法是常用的数值积分方法?A. 矩形法则B. 辛普森规则C. 梯形规则D. 高斯-勒让德法则4. 龙格-库塔法是常用于求解常微分方程的数值解法。
以下哪个选项是描述龙格-库塔法的特点?A. 该方法是一种多步法B. 该方法是一种多项式插值法C. 该方法是一种单步法D. 该方法是一种数值积分法5. 用有限差分法求解偏微分方程时,通常需要进行网格剖分。
以下哪个选项是常用的网格剖分方法?A. 多边形剖分法B. 三角剖分法C. 矩形剖分法D. 圆形剖分法二、解答题1. 将函数f(x) = e^x 在区间[0, 1]上用复化梯形规则进行数值积分,分为6个子区间,求得的近似积分值为多少?解:将区间[0, 1]等分为6个子区间,每个子区间的长度为h = (1-0)/6 = 1/6。
根据复化梯形规则的公式,近似积分值为:I ≈ (1/2) * h * [f(0) + 2f(1/6) + 2f(2/6) + 2f(3/6) + 2f(4/6) + 2f(5/6) +f(1)]≈ (1/2) * (1/6) * [e^0 + 2e^(1/6) + 2e^(2/6) + 2e^(3/6) + 2e^(4/6) +2e^(5/6) + e^1]2. 使用二分法求解方程 x^3 - 3x + 1 = 0 在区间[1, 2]上的根。
要求精确到小数点后三位。
解:首先需要判断方程在区间[1, 2]上是否存在根。
数值分析整理版试题及答案
例 1、已知函数表 求 f(1) 由题可知 (x - x 1 )(x - x 2 ) = (x -1)(x - 2) = 1 ( )( )(x- x 0 )(x - x 2 ) (x + 1)(x - 2) 1 (x 1 - x 0 )(x 1 - x 2 ) (1 + 1)(1 - 2) l 1(x ) = = = - (x + 1)(x - 2) (x + 1)(x -1) = 1 x -1 x + 1(x - x 0 )(x - x 1 ) (x 2 - x 0 )(x 2 - x 1 ) (2 + 1)(2 -1) 3 L 2 (x ) = ∑ y k l k (x ) (x -1)(x - 2)+ 0 ⨯ ⎡⎢- (x + 1)(x - 2)⎤⎥ + 4 ⨯ (x + 1)(x -1)= - (x -1)(x - 2)+ (x + 1)(x -1) 5 2 x - f [x 0 , x 1]= = f (x 0 )- f (x 1 ) -3 - 0 3f [x 1, x 2 ]== = 4 - 4 f [x 0 , x 1]- f [x 1, x 2 ] f [x 0 , x 1, x 2 ]= = 2解:插值基函数分别为l 0 (x ) =(x 0 - x 1 )(x 0 - x 2 ) (-1 -1)(-1 - 2) 62l 2 (x ) = = ( )( )故所求二次拉格朗日插值多项式为2k =0= -3⨯ 1 1 1 6 ⎣ 2 ⎦ 31 423 = 6 x + 3 72 3(2)一阶均差、二阶均差分别为=x 0 - x 1 -1 -1 2 f (x 1 )- f (x 2 ) 0 - 4x 1 - x 2 1 - 23x 0 - x 2 -1 - 2=5 61(x + 1)+ (x + 1)(x -1)5 2 x - (f ,ϕ 0 ) = ⎰ (x 2 + 3x + 2)dx =1 232 6( f ,ϕ 1 ) = ⎰ x (x 2 + 3x +2)dx = ⎢ 1 ⎢1 ⎢ ⎥ ⎢ ⎥ = ⎢ 9 ⎥ ,经过消元得 ⎢⎢0 ⎥⎢ ⎥ =⎢ 1 ⎥⎣ ⎦ ⎣ ⎦ ⎣故所求Newton 二次插值多项式为P 2 (x ) = f (x 0 )+ f [x 0 , x 1](x - x 0 )+ f [x 0 , x 1, x 2 ](x - x 0 )(x - x 1 )= -3 + 3 52 63 7= x + 62 3 例 2、 设 f (x ) = x2 + 3x + 2 , x ∈[0,1] ,试求 f (x ) 在[0,1]上关于 ρ (x ) = 1 ,Φ = span {1, x }的最佳平方逼近多项式。
数值分析试题与答案
一. 填空题(本大题共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 3 i y 2 4 12 i y '3并估计误差。
(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 λ'=⎧⎨=⎩考察欧拉显式格式的收敛性。
(10分)《数值分析》(A )卷标准答案(2009-2010-1)一. 填空题(每小题3分,共12分) 1. ()1200102()()()()x x x x l x x x x x --=--; 2.7;3. 3,8;4. 2n+1。
数值分析期末试题及答案
数值分析期末试题及答案试题一: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) 迭代法和直接法是数值计算中常用的两种方法。
迭代法通过反复迭代逼近解,直到满足所需精度为止;而直接法则通过一系列代数运算直接得到解。
数值分析题库答案(含详细解题步骤)
第 1 页/共 22 页1. 正方形的边长大约为100cm ,应怎样测量才干使面积误差不超过1cm 22. 已测得某场地长l 的值为110=*l m ,宽d 的值为80=*d m ,已知 2.0≤-*l l m,1.0≤-*d d m, 试求面积ld s =的绝对误差限与相对误差限.3.为使π的相对误差小于0.001%,至少应取几位有效数字?4.设x的相对误差界为δ,求n x的相对误差界.5.设有3个近似数a=2.31,b=1.93,c=2.24,它们都有3位有效数字,试计算p=a+bc的误差界和相对误差界,并问p的计算结果能有几位有效数字?第 3 页/共 22 页6. 已知333487.034.0sin ,314567.032.0sin ==,请用线性插值计算3367.0sin 的值,并预计截断误差.7. 已知sin0.32=0.314567, sin0.34=0.333487, sin0.36= 0.352274,用抛物插值计算sin0.3367的值, 并预计误差.8. 已知16243sin ,sin πππ===请用抛物插值求sin50的值,并预计误差9. . .6,8,7,4,1)(,5,4,3,2,1求四次牛顿插值多项式时设当==i i x f x第 5 页/共 22 页10. 已知4)2(,3)1(,0)1(=-=-=f f f , 求函数)(x f 过这3点的2次牛顿插 值多项式.11. 设x x f =)(,并已知483240.1)2.2(,449138.1)1.2(,414214.1)0.2(===f f f ,试用二次牛顿插值多项式计算(2.15)f 的近似值,并研究其误差12. 设],[)(b a x f 在上有四阶延续导数,试求满意条件)2,1,0()()(==i x f x P i i 及)()(11x f x P '='的插值多项式及其余项表达式.13. 给定3201219(),,1,,44f x x x x x ====试求()f x 在1944⎡⎤⎢⎥⎣⎦,上的三次埃尔米特插值多项式()P x ,使它满意11()()(0,1,2),()(),i i P x f x i P x f x ''===并写出余项第 7 页/共 22 页表达式.14. 设],1,0[,23)(2∈++=x x x x f 试求)(x f 在]1,0[上关于,,1{,1)(x span x =Φ=ρ}2x 的最佳平方逼近多项式15.已知实验数据如下:用最小二乘法求形如y=a+bx2的拟合曲线,并计算均方误差.16.已知数据表如下第 9 页/共 22 页x i 1 2 3 4 5 y iωi4 4.56 8 8.5 2 1 3 1 1试用最小二乘法求多项式曲线与此数据组拟合17. .1)(},1{span ,1]41[)(的最佳平方逼近多项式中的关于上的在在求==Φ=x x x x f ρ18. 决定求积公式⎰++≈10110)1()(32)0()(f A x f f A dx x f 中的待定参数110,,A x A , 使其代数精度尽量高,并指出所决定的求积公式的代数精度.19. 用复化辛普森公式计算积分⎰=10dx e I x , 问区间[0,1]应分多少等分才干使截断误差不超过?10215-⨯第 11 页/共 22 页20. 利用下表中给出的数据,分离用复化梯形公式和复化辛甫生公式计算定积分dx x I ln 21⎰=的近似值(要求结果保留到小数点后六位)21. 用复化梯形公式和复化辛甫生公式计算积分⎰=6.28.1)(dx x f I ,函数)(x f 在某些节点上的值如下图:(本题共14分)22. 决定公式⎰+≈101100)()()(x f A x f A dx x f x 的系数1010,,,x x A A ,使其具有最高代数精度23. 决定求积公式⎰++≈1110)1()(32)0()(f A x f f A dx x f 中的待定参数110,,A x A ,使其代数精度尽量高,并指出所决定的求积公式的代数精度第 13 页/共 22 页24.用LU 分解法求解以下方程组 (10分)123123142521831520x x x ⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭25.用LU 分解法求解以下方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛8892121514131615141321x x x26. 用LU 分解法求解以下方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫⎝⎛542631531321321x x x27. 设方程组b Ax =,其中⎪⎪⎪⎭⎫⎝⎛-=220122101A ,Tb ⎪⎭⎫ ⎝⎛-=32,31,21, 已知它有解Tx ⎪⎭⎫⎝⎛-=0,31,21,若右端有小扰动61021-∞⨯=bδ,试预计由此引起的解的相对误差.第 15 页/共 22 页28. 设方程组b Ax =,其中212 1.0001A -⎛⎫= ⎪-⎝⎭,11.0001b -⎛⎫= ⎪⎝⎭,当右端向量b 有误差00.0001δ⎛⎫= ⎪⎝⎭b 时,试预计由此引起的解的相对误差(用∞范数计算)29. 给定b Ax =,其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111a a a a a a A 证实:(1) 当121<<-a 时,A 对称正定,从而GS 法收敛. (2) 惟独当2121<<-a 时,J 法收敛.30. 对于线性方程组⎪⎩⎪⎨⎧-=+-=-+=+1242043 16343232121x x x x x x x ,列出求解此方程组的Jacobi 迭代格式,并判断是否收敛。
数值分析试卷及答案
数值分析试卷及答案数值分析试卷一、选择题(共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. 数值解线性方程组的常用迭代法有_________方法和_________方法。
数值分析期末考试题及答案
数值分析期末考试题及答案一、选择题(每题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 \) 处插值,并求出插值多项式。
数值分析练习题附答案
目录一、绪论------------------------------------------------------------------------------------- 2-2二、线性方程组直接解法列主元高斯LU LDL T GG T-------------------- 3-6二、线性方程组迭代法----------------------------------------------------------------- 7-10 三、四、非线性方程组数值解法二分法不动点迭代---------------------- 11-13五、非线性方程组数值解法牛顿迭代下山弦截法----------------- 14-15六、插值线性插值抛物线插值------------------------------------------------ 16-18七、插值Hermite插值分段线性插值-----------------------------------------19-22八、拟合------------------------------------------------------------------------------------ 23-24九、数值积分----------------------------------------------------------------------------- 25-29十、常微分方程数值解法梯形欧拉改进----------------------------------- 30-32 十一、常微分方程数值解法龙格库塔------------------------------------------ 33-35绪论1-1 下列各数都是经过四舍五入得到的近似值 ,试分别指出它们的绝对误差限,相对误差限和有效数字的位数.X 1 =5.420, X 2 =0.5420, X 3 =0.00542, X 4 =6000, X 5 =0.6×105注:将近似值改写为标准形式X 1 =(5*10-1+4*10-2+2*10-3+0*10-4)*101 即n=4,m=1 绝对误差限|△X 1|=|X *1-X 1|≤ 12×10m-n =12×10-3 相对误差限|△r X 1|= |X∗1−X1||X∗1|≤|X∗1−X1||X1|= 12×10-3/5.4201-2 为了使101/2 的相对误差小于0.01%, 试问应取几位有效数字?1-3 求方程x 2 -56x+1=0的两个根, 使它们至少具有4位有效数字( √783≈27.982)注:原方程可改写为(x-28)2=783线性方程组解法(直接法)2-1用列主元Gauss消元法解方程组解:回代得解:X1=0 X2=-1 X3=12-2对矩阵A进行LU分解,并求解方程组Ax=b,其中解:(注:详细分解请看课本P25)A=(211132122)→(211(1/2)5/23/2(1/2)3/23/2)→(2111/25/23/21/2(3/5)3/5)即A=L×U=(11/211/23/51)×(2115/23/23/5)先用前代法解L y=P b 其中P为单位阵(原因是A矩阵未进行行变换)即L y=P b 等价为(11/211/23/51)(y1y2y3)=(111)(465)解得 y 1=4 y 2=4 y 3=35再用回代解Ux =y ,得到结果x即Ux =y 等价为(2115/23/23/5)(x 1x 2x 3)=(y 1y 2y 3)=(443/5) 解得 x 1=1 x 2=1 x 3=1即方程组Ax=b 的解为x =(111)2-3 对矩阵A 进行LDL T 分解和GG T 分解,求解方程组Ax=b,其中A=(164845−48−422) , b =(123)解:(注:课本 P 26 P 27 根平方法)设L=(l i j ),D=diag(d i ),对k=1,2,…,n,其中d k =a kk -∑l kj 2k−1j=1d jl ik =(a ik −∑l ij l kj k−1j=1d j )/ d k 即d 1=a 11-∑l 1j 20j=1d j =16-0=16因为 l 21=(a 21−∑l 2j l 1j 0j=1d j )/ d 1=a 21/ d 1=416=14 所以d 2=a 22-∑l 2j 21j=1d j =5-(14)2d 1=4同理可得d 3=9 即得 D=(1649)同理l 11=(a 11−∑l ij l 1j 0j=1d j )/ d 1=1616=1=l 22=l 33 l 21=(a 21−∑l 2j l 1j 0j=1d j )/ d 1=416=14 l 31=(a 31−∑l 3j l 1j 0j=1d j )/ d 1=816=12 l 32=(a 32−∑l 3j l 2j 1j=1d j )/ d 2=−4−12×14×164=−64=-32即L=(114112−321) L T=(114121−321) 即LDL T分解为A=(114112−321)(1649)(114121−321)解解:A=(164845−48−422)→(41212−32−33)故得GG T分解:A=(4122−33)(4122−33) LDL T分解为A=(114112−321)(1649)(114121−321) 由(114112−321)(y 1y 2y 3)=(123) ,得(y 1y 2y 3)=(0.250.8751.7083)再由(4122−33)(x 1x 2x 3)=(0.250.8751.7083) ,得(x 1x 2x 3)=(−0.54511.29160.5694)2-4 用追赶法求解方程组:解:(4−1−14−1−14−1−14−1−14)→(4−14−1154−415−15615−1556−120956−56209−1780209)由(4−1154−15615−120956−1780209)(y1y2y3y4y5)=(100200),得(y1y2y3y4y5)=(256.66671.785700.4784753.718)再由(1−141−4151−15561−562091)(x1x2x3x4x5)=(256.66671.785700.4784753.718),得(x1x2x3x4x5)=(27.0518.20525.769314.87253.718)线性方程组解法(迭代法)2-1 设线性方程组{4x 1−x 2+2x 3=1−x 1−5x 2+x 3=22x 1+x 2+6x 3=3(1) 写出Jacobi 法和SOR 法的迭代格式(分量形式) (2) 讨论这两种迭代法的收敛性(3) 取初值x (0)=(0,0,0)T ,若用Jacobi 迭代法计算时,预估误差 ||x*-x (10)||∞ (取三位有效数字)解:(1)Jacobi 法和SOR 法的迭代格式分别为Jacobi 法迭代格式SOR(2)因为A 是严格对角占优矩阵,但不是正定矩阵,故Jacobi 法收敛,SOR 法当0<ω≤1时收敛.⎪⎪⎪⎩⎪⎪⎪⎨⎧+--=-+-=+-=+++216131525151412141)(2)(1)1(3)(3)(1)1(2)(3)(2)1(1k k k k k k k k k x x x x x x xx x ⎪⎪⎪⎩⎪⎪⎪⎨⎧-++-=+-+-=+-+-+=++++++)216131()525151()412141()(3)1(2)1(1)(3)1(3)(3)(2)1(1)(2)1(2)(3)(2)(1)(1)1(1k k k k k k k k k k k k k k k x x x x x x x x x x x x x x x ωωω(3)由(1)可见||B ||∞=3/4,且取x (0)=(0,0,0)T ,经计算可得x (1)=(1/4,-2/5,1/2)T ,于是||x (1)-x (0)||∞=1/2,所以有2-2 设方程组为{5x 1+2x 2+x 3=−12−x 1+4x 2+2x 3=202x 1−3x 2+10x 3=3试写出其Jacobi 分量迭代格式以及相应的迭代矩阵,并求解。
《数值分析》练习题及答案解析
《数值分析》练习题及答案解析一、单选题1. 以下误差公式不正确的是( D )A .()1212x x x x ∆-≈∆-∆B .()1212x x x x ∆+≈∆+∆C .()122112x x x x x x ∆≈∆+∆D .1122()x x x x ∆≈∆-∆ 2. 已知等距节点的插值型求积公式()()352kkk f x dx A f x =≈∑⎰,那么3kk A==∑( C )A .1 B. 2 C.3 D. 4 3.辛卜生公式的余项为( c )A .()()32880b a f η-''-B .()()312b a f η-''-C .()()()542880b a f η--D .()()()452880b a f η--4. 用紧凑格式对矩阵4222222312A -⎡⎤⎢⎥=-⎢⎥⎢⎥--⎣⎦进行的三角分解,则22r =( A ) A .1 B .12C .–1D .–25. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足( D ) A .()00l x =0,()110l x = B . ()00l x =0,()111l x = C .()00l x =1,()111l x = D . ()00l x =1,()111l x =6. 用二分法求方程()0f x =在区间[],a b 上的根,若给定误差限ε,则计算二分次数的公式是n ≥( D )A .ln()ln 1ln 2b a ε-++ B. ln()ln 1ln 2b a ε-+-C.ln()ln 1ln 2b a ε--+ D. ln()ln 1ln 2b a ε--- 7.若用列主元消去法求解下列线性方程组,其主元必定在系数矩阵主对角线上的方程组是( B )A .123123123104025261x x x x x x x x x -+=⎧⎪-+=⎨⎪-+=-⎩ B 。
数值分析试题及答案
数值分析试题及答案一、选择题(每题2分,共20分)1. 以下哪个算法是数值分析中用于求解线性方程组的直接方法?A. 牛顿法B. 高斯消元法C. 梯度下降法D. 蒙特卡洛方法答案:B2. 插值法中,拉格朗日插值法和牛顿插值法的共同点是:A. 都是多项式插值B. 都使用差商C. 都只适用于等距节点D. 都需要预先知道所有数据点答案:A3. 在数值积分中,辛普森(Simpson)公式比梯形公式的误差:A. 更大B. 更小C. 相同D. 无法比较答案:B4. 以下哪个是数值稳定性分析中常用的方法?A. 条件数B. 收敛性C. 收敛速度D. 误差分析答案:A5. 在求解常微分方程的数值解时,欧拉方法属于:A. 单步法B. 多步法C. 隐式方法D. 显式方法答案:A6. 以下哪个是数值分析中求解非线性方程的迭代方法?A. 高斯-约当消元法B. 牛顿-拉弗森方法C. 雅可比迭代法D. 高斯-赛德尔迭代法答案:B7. 线性插值公式中,如果给定两个点\( (x_0, y_0) \)和\( (x_1, y_1) \),插值多项式是:A. \( y = y_0 + \frac{y_1 - y_0}{x_1 - x_0}(x - x_0) \)B. \( y = y_0 + \frac{y_1 - y_0}{x_0 - x_1}(x - x_0) \)C. \( y = y_0 + \frac{x - x_0}{x_1 - x_0}(y_1 - y_0) \)D. \( y = y_1 + \frac{x_1 - x}{x_1 - x_0}(y_0 - y_1) \)答案:C8. 以下哪个是数值分析中用于求解特征值问题的算法?A. 幂法B. 共轭梯度法C. 牛顿法D. 欧拉法答案:A9. 在数值微分中,使用有限差分法来近似导数时,中心差分法的误差:A. 与步长成正比B. 与步长的平方成正比C. 与步长的立方成正比D. 与步长的四次方成正比答案:B10. 以下哪个是数值分析中用于求解线性最小二乘问题的算法?A. 梯度下降法B. 牛顿法C. 奇异值分解法D. 共轭梯度法答案:C二、简答题(每题10分,共30分)1. 简述数值分析中病态问题的特点及其对算法的影响。
数值分析试卷及答案
数值分析试卷及答案**注意:以下是一份数值分析试卷及答案,试卷和答案分别按照题目和解答的格式排版,以确保整洁美观,语句通顺。
**---数值分析试卷一、选择题(每题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:从第一行开始,选取第一个非零元素作为主元,然后通过行变换将其它行的该列元素消去。
数值分析试题及答案
1. 和分别作为π的近似数具有( )和( )位有效数字. A .4和3 B .3和2 C .3和4 D .4和42. 已知求积公式()()211211()(2)636f x dx f Af f ≈++⎰,则A =( )A . 16B .13C .12D .233. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足( )A .()00l x =0,()110l x = B .()00l x =0,()111l x =C .()00l x =1,()111l x = D .()00l x =1,()111l x =4. 设求方程()0f x =的根的牛顿法收敛,则它具有( )敛速。
A .超线性B .平方C .线性D .三次5. 用列主元消元法解线性方程组1231231220223332x x x x x x x x ++=⎧⎪++=⎨⎪--=⎩ 作第一次消元后得到的第3个方程( ).A .232x x -+= B .232 1.5 3.5x x -+=C .2323x x -+= D .230.5 1.5x x -=-单项选择题答案二、填空题(每小题3分,共15分)1. 设TX )4,3,2(-=, 则=1||||X ,2||||X = .2. 一阶均差()01,f x x =3. 已知3n =时,科茨系数()()()33301213,88C C C ===,那么()33C = 4. 因为方程()420x f x x =-+=在区间[]1,2上满足 ,所以()0f x =在区间内有根。
5. 取步长0.1h =,用欧拉法解初值问题()211yy yx y ⎧'=+⎪⎨⎪=⎩的计算公式 .填空题答案1. 9和292.()()0101f x f x x x --3. 18 4. ()()120f f < 5. ()1200.11.1,0,1,210.11k k y y k k y +⎧⎛⎫⎪ ⎪=+⎪ ⎪=+⎨⎝⎭⎪=⎪⎩得 分 评卷人三、计算题(每题15分,共60分)1. 已知函数211y x =+的一组数据:求分段线性插值函数,并计算()1.5f 的近似值.计算题1.答案()101x L x -=-()12x L x -=-()10.8L x ⎧-⎪=⎨⎪⎩()1.50.8L =2. 已知线性方程组1231231231027.21028.35 4.2x x x x x x x x x --=⎧⎪-+-=⎨⎪--+=⎩(1) 写出雅可比迭代公式、高斯-塞德尔迭代公式;(2) 对于初始值()()00,0,0X =,应用雅可比迭代公式、高斯-塞德尔迭代公式分别计算()1X(保留小数点后五位数字).计算题2.答案1.解 原方程组同解变形为 1232133120.10.20.720.10.20.830.20.20.84x x x x x x x x x =++⎧⎪=-+⎨⎪=++⎩雅可比迭代公式为()()()()()()()()()1123121313120.10.20.720.10.20.830.20.20.84m m m m m m m m m x x x x x x x x x +++⎧=++⎪⎪=-+⎨⎪=++⎪⎩(0,1...)m =高斯-塞德尔迭代法公式()()()()()()()()()1123112131113120.10.20.720.10.20.830.20.20.84m m m m m m m m m x x x x x x x x x ++++++⎧=++⎪⎪=-+⎨⎪=++⎪⎩ (0,1...)m =用雅可比迭代公式得()()10.72000,0.83000,0.84000X =用高斯-塞德尔迭代公式得()()10.72000,0.90200,1.16440X =3. 用牛顿法求方程3310x x --=在[]1,2之间的近似根 (1)请指出为什么初值应取2 (2)请用牛顿法求出近似根,精确到.计算题3.答案4. 写出梯形公式和辛卜生公式,并用来分别计算积分1011dx x +⎰.计算题4.答案四、证明题(本题10分)确定下列求积公式中的待定系数,并证明确定后的求积公式具有3次代数精确度()()()()1010hhf x dx A f h A f A f h --=-++⎰证明题答案证明:求积公式中含有三个待定系数,即101,,A A A -,将()21,,f x x x =分别代入求积公式,并令其左右相等,得一、 填空(共20分,每题2分)1. 设2.3149541...x *=,取5位有效数字,则所得的近似值x= .2.设一阶差商()()()21122114,321f x f x f x x x x --===---,()()()322332615,422f x f x f x x x x --===--则二阶差商()123,,______f x x x =3. 设(2,3,1)TX =--, 则2||||X = ,=∞||||X 。
数值分析试题答案
数值分析试题答案一、选择题1. 以下哪个数值方法用于求解线性方程组?A. 牛顿法B. 高斯消元法C. 欧几里得算法D. 拉格朗日插值法答案:B2. 在数值分析中,舍入误差通常是由什么引起的?A. 人为计算错误B. 计算机表示数字的限制C. 测量误差D. 数据输入错误答案:B3. 插值和拟合的区别在于:A. 插值通过所有数据点,而拟合不通过B. 拟合通过所有数据点,而插值不通过C. 插值是线性的,拟合是非线性的D. 插值是精确的,拟合是近似的答案:A4. 以下哪种方法最适合求解非线性方程?A. 雅可比迭代法B. 牛顿-拉弗森方法C. 托马斯算法D. 布雷尔-史密斯算法答案:B5. 在数值分析中,条件数用于衡量什么?A. 方程组解的存在性B. 方程组解的唯一性C. 方程组解的稳定性D. 方程组解的精确性答案:C二、填空题1. 在数值分析中,__________误差指的是由于计算机舍入而产生的误差,而__________误差指的是由于数据不精确或截断而产生的误差。
答案:截断;舍入2. 线性方程组的矩阵表示为__________,其中A是系数矩阵,x是变量向量,b是常数向量。
答案:Ax = b3. 牛顿法求解非线性方程时,需要计算函数的__________。
答案:导数4. 拉格朗日插值法通过构建一个多项式来近似数据点,该多项式的每一段都与数据点的__________相匹配。
答案:切线5. 为了减少数值分析中的误差,通常采用__________方法来提高计算的精度。
答案:增量三、简答题1. 请简述高斯消元法的基本思想及其在求解线性方程组中的应用。
高斯消元法的基本思想是通过行变换将系数矩阵转化为阶梯形矩阵,进而简化方程组的求解过程。
在求解线性方程组时,首先将增广矩阵进行行变换,使得主元下方的元素为零,然后通过回代过程逐步求解出未知数。
2. 描述牛顿-拉弗森方法求解非线性方程的迭代过程。
牛顿-拉弗森方法是一种迭代求解非线性方程的方法。
数值分析试卷及答案
模 拟 试 卷(一)一、填空题(每小题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 设方程组,试考察解此方程组的雅可比迭代法及高斯-赛德尔迭代法的收敛性。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
二
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按行严格对角占优,从而高斯-赛德尔法收敛。
10设A为对称正定矩阵,考虑迭代格式
求证:(1)对任意初始向量,收敛;
(2)收敛到的解。
证明(1)所给格式可化为
这里存在是因为,由A对称正定,,故也对称正定。
设迭代矩阵的特征值为,为相应的特征向量,则
与做内积,有
因正定,故,从而,格式收敛。
(2)设收敛到,则即,
即收敛到的解。
三
1 设且.求证:
证明以和为插值节点建立的不超过一次的插值多项式
应用插值余项公式有
2求一个次数不高于4次的多项式,使它满足
.
解法一(待定参数法)满足的Hermite插值多项式为
设,令得
于是
解法二(带重节点的Newton插值法)建立如下差商表
这样可以写出Newton插值公式
3设,在上取,按等距节点求分段线性插值函数,计算各节点间中点处与的值,并估计误差.
解步长,.在区间上的线性插值函数
分段线性插值函数定义如下
,
各区间中点的函数值及插值函数值如表所示
估计误差:在区间上
而
令得的驻点,于是
故有结论
,
右端与无关,于是有
,
四
1确定参数和,使得积分取得最小值,并计算该最小值.
解本题实质上是求,关于权函数的二次最佳平方逼近多项式.
选切比雪夫多项式为基函数进行计算:
于是得的二次最佳平方逼近多项式
进而有参数.
最小值就是平方误差:
2 对彗星1968Tentax 的移动在某个极坐标系下有如表所示的观察数据.
假设忽略来自行星的干扰,坐标应满足 其中
为参数,
为离心率,试用最小二乘法拟合
和
,并给出平方误差.
解 由于关于参数和是非线性的,变形为
,这样有下表的数据.
记,得拟合模型.
求解法方程组
得
进而有,拟合方程为
平方误差为
3 求函数在指定区间上关于的最佳平方逼近多项式.
解对做线性变换,即
利用勒让德正交多项式为基建立的一次最佳平方逼近多项式
的最佳平方逼近为
五
1确定中的待定参数,使其代数精确度尽量高,并指明求积公式所具有的代数精确度。
解令,代入公式两端并令其相等,得
解得
令,得
令,得故求积公式具有3次代数精确度。
2计算积分,若复化梯形公式,问区间应分多少等份才能使截断误差不超过?若改用复化辛普森公式,要达到同样精确度,区间应分多少等份?
解由于,故对复化梯形公式,要求
即。
取,即将区间分为213等份时,用复化梯形
公式计算,截断误差不超过。
用复化辛普森公式,要求
即。
取,即将区间等分为8等份时,复化辛普森公式可
达精度。
3确定求积公式
中的系数,使代数精确度尽量高,并给出的表达式。
公式中。
解这是一个带权的且带导数值的求积公式。
为了积分方便,设该求积公式对准确成立,得
化简得
解得
又因为
故求积公式
具有3次代数精确度。
下面估计求积公式的余项。
设在上三次插值多项式为,即满足。
因前述求积公式具有3次代数精确度,故它对于是准确成立的,且
因此有
注意到在上不变号,故余项
4已知。
(1)推导以这3个点作为求积节点在上的插值型求积公式;(2)指明求积公式所具有的代数精确度;
(3)用所求公式计算。
解(1)过这3个点的插值多项式
故
其中
故所求的插值型求积公式为
(2)上述求积公式是由二次插值函数积分而来,故至少具有2次代数精确度。
再将代入上述求积公式,有
故上述求积公式具有3次代数精确度。
(3)
由于该求积公式具有3次代数精确度,从而为的精确度。
5设。
求证:(1)
(2)
(提示:直接使用泰勒展开即可得证)
七
1 对于迭代函数,试讨论:
(1)当为何值时,产生的序列收敛于;
(2)取何值时收敛最快?
(3)分别取计算的不动点,要求
解(1),根据定理7.3,当,亦即时迭代收敛。
(2)由定理7.4知,当,即时迭代至少是二阶收敛的,收敛最快。
(3)分别取,并取,迭代计算结果如表7-4所示。
0 1.2 0 1.2
此时都达到。
事实上,
2(牛顿迭代法收敛性定理)设在上具有二阶连续导数,且满足条件
(1);
(2)在上;
(3)满足。
则由牛顿迭代法产生的序列单调收敛于在内的唯一实根,并且是平方收敛的。
证明因在上连续,由条件(1)知,方程在内有根。
又由条件(2)知在上恒正或恒负,所以在上严格单调,因而是
在内的唯一实根。
条件(1)(2)共有四种情形:
(1)
(2)
(3)
(4)
仅就(1)进行定理证明,其余三种情况证明方法类似。
由,可知,再由知单增且。
又由牛顿迭代法知
由台劳展开的
其中介于,之间。
利用得
由以及前面证明的有
一般地,设,则必有且
再由台劳
及,得
根据归纳法原理数列单调下降有下界,因此有极限。
设,对迭代式两端取的极限,并利用,的连续性知即。
由上述证明知,有关系式,即对于单根,牛顿迭代法是平方收敛的。
3给定函数,对于一切,存在且,证明对于范围内的任意定数,迭代过程均收敛于的根.证明由于,为单调增函数,故方程的根是唯一的(假定方程有根)。
迭代函数,。
由
及得,故。
由此可得
即
4设,试确定函数和,使求解且以为迭代函数的迭代法至少三阶收敛。
解要求三阶收敛到的根,根据定理7.4,应有
于是由
得
故取
即迭代至少三阶收敛。