数值分析习题.doc

习题11 填空题(1) 为便于算法在计算机上实现，必须将一个数学问题分解为 的 运算；(2) 在数值计算中为避免损失有效数字，尽量避免两个 数作减法运算；为避免误差的扩大，也尽量避免分母的绝对值 分子的绝对值；(3) 误差有四大来源,数值分析主要处理其中的 和 ；(4) 有效数字越多,相对误差越 ；2. 用例1.4的算法计算,迭代3次,计算结果保留4位有效数字.3. 推导开平方运算的误差限公式,并说明什么情况下结果误差不大于自变量误差.4. 以下各数都是对准确值进行四舍五入得到的近似数,指出它们的有效数位、误差限和相对误差限.5. 证明1.2.3之定理1.1.6. 若钢珠的的直径d的相对误差为1.0%,则它的体积V的相对误差将为多少。（假定钢珠为标准的球形）7. 若跑道长的测量有0.1%的误差,对400m成绩为60s的运动员的成绩将会带来多大的误差和相对误差.8. 为使的近似数相对误差小于0.05%,试问该保留几位有效数字.9. 一个园柱体的工件,直径d为10.250.25mm,高h为40.001.00mm,则它的体积V的近似值、误差和相对误差为多少10 证明对一元函数运算有并求出时的值，从而说明在时是病态问题11. 定义多元函数运算求出的表达式，并说明全为正数时，计算是稳定的，有正有负时，误差难以控制12. 下列各式应如何改进,使计算更准确：习题21. 填空题(1) Gauss消元法求解线性方程组的的过程中若主元素为零会发生 ;. 主元素的绝对值太小会发生 ;(2) Gauss消元法求解线性方程组的计算工作量以乘除法次数计大约为 . 平方根法求解对称正定线性方程组的计算工作量以乘除法次数计大约为 ;(3) 直接LU分解法解线性方程组时的计算量以乘除法计为 , 追赶法解对角占优的三对角方程组时的计算量以乘除法计为 ;(4) , , ;(5) , ;(6) , ;2用Gauss消元法求解下列方程组 , 3用列主元消元法解下列方程组 4. 用GaussJordan消元法求：5用直接分解方法求1题中两个矩阵的分解，并求解此二方程组6用平方根法解方程组 7 用追赶法解三对角方程组 8证明：(1)单位下三角阵的逆仍是单位下三角阵(2)两个单位下三角阵的乘积仍是单位下三角阵9由，(见(2.18)式)，证明:10证明向量范数有下列等价性质：11求下列矩阵的 12求 13证明：(1)若是正交矩阵，即, 则;(2)若是对称正定阵， 是的最大特征值， 是最小特征值，则.习题31. 填空题:(1) 当A具有严格对角线优势或具有对角优势且 时，线性方程组Ax=b用Jacobi迭代法和GaussSeidel迭代法均收敛;(2) 当线性方程组的系数矩阵A对称正定时, 迭代法收敛.(3) 线性方程组迭代法收敛的充分必要条件是迭代矩阵的 小于1; SOR法收敛的必要条件是 ;(4) 用迭代法求解线性方程组,若q = r (B), q 时不收敛, q接近 时收敛较快, q接近 时收敛较慢;(5) ; ; ; 2用Jacobi迭代法和GaussSeidel迭代法求解方程组(1) ； (2) 各分量第三位稳定即可停止3用SOR法解方程组，取，与取 (即Gauss-Seidel法)作比较4下面是一些方程组的系数阵，试判断它们对Jacobi迭代法，Gauss-Seidel迭代法的收敛性(1)； (2)；(3)； (4)；(5) ； (6) 5方程组证明用Jacobi迭代法收敛的充要条件是：6设 (1)若正定，的取值范围；(2)若Jacobi迭代法收敛，的取值范围习题41. 填空题：(1) 幂法主要用于求一般矩阵的 特征值，Jacobi旋转法用于求对称矩阵的 特征值；(2) 古典的Jacobi法是选择 的一对 元素将其消为零;(3) QR方法用于求 矩阵的全部特征值，反幂法加上原点平移用于一个近似特征值的 和求出对应的 2用幂法求矩阵， 按模最大的特征值和对应的特征向量，精确到小数三位3已知： 取t =15，作原点平移的幂法，求按模最大特征值4 用反幂法加原点平移求最接近12的特征值与相应的特征向量，迭代三次5若的特征值为是一实数，证明：是的特征值，且特征向量不变6已知求平面反射阵使，即使的1，3两个分量化零7 试用Jacobi旋转法求作一次旋转，消去最大的非对角元，写出旋转矩阵，求出角和结果8设 已知是的特征值，相应的特征向量为，证明也是的特征值，相应的特征向量为9 证明定理4.510 证明(421)中的和相似习题51填空题(1) 用二分法求方程在0,1内的根，迭代一次后,根的存在区间为 ，迭代两次后根的存在区间为 ；(2) 设可微，则求方程根的Newton迭代格式为 ；(3) ，若要使迭代格式局部收敛到，则C取值范围为 ；(4) 用迭代格式求解方程的根，要使迭代序列是二阶收敛，则= ；(5) 迭代格式收敛于根= ，此迭代格式是 阶收敛的2证明Newton迭代格式(5.10)满足3. 方程的根全正实根，试用逐次扫描法(h=1)，找出它的全部实根的存在区间，并用二分法求出最大实根，精确到0.014用二分法求下列方程的根，精度(1) (2) 5用迭代法求的正根，简略判断以下三种迭代格式：(1) ； (2) ； (3) 在附近的收敛情况，并选择收敛的方法求此根精度6. 方程 (1) 证明它在(0,1)区间有且只有一个实根；(2) 证明，在(0,1)区间内收敛；(3) 用Newton迭代法求出此根,精确到5位有效数字7对方程，分别用(1) Newton法；(2) 割线法求其根精度8用迭代法求下列方程的最小正根(1) ； (2) ； (3) 9设有方程 (1) 以，找出根的全部存在区间；(2) 验证在区间0,1上Newton法的区间收敛定理条件不成立；(3) 验证取, 用Newton法不收敛；(4) 用Newton下山法，取求出根的近似值，精度10分别用Jacobi法，GaussSeidel法求解非线性方程组 在(1.5,0.7)附近的根,精确到11分别用Newton法，简化Newton法求解非线性方程组 在(0,1)附近的根,精确到习题61填空题(1) 设，则 ，= ，= ； (2) 设是以节点0,1,2,n的Lagrange插值基函数，则 ； (3) 设 ， ，的二次Newton插值多项式为 2已知函数的数据如下-0.6-0.4-0.200.20.40.60.6976760.8521140.96078910.9607890.8521140.697676试用二次，三次插值计算=0.35，=0.55的近似函数值，使其精度尽量地高3利用在及处的值，求的近似值，并估计误差4利用数据00.20.40.60.81.000.199560.396460.588130.772100.94608计算积分, 当=0.45时的的取值5试用Newton插值求经过点(-3,-1),(0,2),(3,-2),(6,10)的三次插值多项式6求满足及的次数不超过2次的插值多项式，并给出其误差表达式7设是互异节点，是Lagrange插值基函数()，证明(1) ； (2) ()；(3) () 8设有如下数据0123436111827试计算此表中函数的差分表，并分别利用Newton向前，向后插值公式求出它的插值多项式9试构造一个三次Hermite插值多项式使其满足10已知函数的数据表0.00.20.40.60.81.00001.221401.491821.822122.22554分别用Newton向前插值公式和向后插值公式求=0.05，=0.42，=0.75的近似值11对函数进行分段线性插值，要求误差不超过，问步长h应如何选取12设有数据0.250.300.390.450.530.50000.54770.62450.67080.7280用三转角插值法求满足下述条件的三次样条插值函数(1) ，(2) ， 13. 证明定理6.6.习题81填空题(1) 个点的插值型数值积分公式的代数精度至少是 ，最高不超过 (2) 梯形公式有 次代数精度，Simpson公式有 次代数精度(3) 求积公式中的参数 时，才能保证该求积公式的代数精度达到最高，最高代数精度为 2确定下列求积公式的求积系数和求积节点，使其代数精度尽量高，并指出其最高代数精度(1) (2) (3) (4) (5) 3分别利用复化梯形公式，复化Simpson公式，复化Cotes公式计算下列积分(1) (=8)(2) (=10)(3) (=10) (4) (=6)(5) (=8)4用Romberg公式计算积分(1) (精度要求)(2) (精度要求)5分别取节点数为2，3，4利用GaussLegendre求积公式计算积分(1) ， (2) ， (3) 6利用Gauss型求积公式，分别取节点数2，3，4计算积分(1) ， (2) 7用节点数为4的GaussLaguerre求积公式和GaussHermite求积公式计算积分的近似值，并与准确值作比较8分别用两点公式与三点公式求在=1.0，=1.2的导数值，并估计误差，其中的数据由下表给出9已知的数据如下取=0.1，=0.2，分别用二点、三点公式计算=2.7处的一阶和二阶导数值习题91填空题(1) 解初值问题的Euler法是 阶方法，梯形方法是 阶方法，标准RK方法是 阶方法(2) 解初值问题时，为保证计算的稳定性，若用经典的四阶RK方法，步长 采用Euler方法，步长h的取值范围为 ，若采用Euler梯形方法，步长h的取值范围为 若采用Adams外推法，步长h的范围为 ，若采用Adams内插法，步长h的取值范围为 (3) 求解初值问题Euler方法的局部截断误差为 Euler梯形方法的局部截断误差为 ， Adams外推法的局部截断误差为 Adams内插法的局部截断误差为 .2对初值问题 试用Euler法取步长=0.1和=0.2计算其近似解，并与准确解进行比较3利用Euler预测校正法和四阶经典RK方法，取