数值分析习题(含答案)(共43页).doc
《数值分析习题(含答案)(共43页).doc》由会员分享,可在线阅读,更多相关《数值分析习题(含答案)(共43页).doc(43页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、精选优质文档-倾情为你奉上第一章 绪论姓名 学号 班级 习题主要考察点:有效数字的计算、计算方法的比较选择、误差和误差限的计算。1 若误差限为,那么近似数0.有几位有效数字?(有效数字的计算)解:,故具有3位有效数字。2 具有4位有效数字的近似值是多少?(有效数字的计算)解:,欲使其近似值具有4位有效数字,必需,即即取(3.14109 , 3.14209)之间的任意数,都具有4位有效数字。3 已知,是经过四舍五入后得到的近似值,问,有几位有效数字?(有效数字的计算)解:,而,故至少具有2位有效数字。故至少具有2位有效数字。4 设,的相对误差为,求的误差和相对误差?(误差的计算)解:已知,则误差
2、为 则相对误差为 5测得某圆柱体高度的值为,底面半径的值为,已知,求圆柱体体积的绝对误差限与相对误差限。(误差限的计算)解:绝对误差限为相对误差限为6 设的相对误差为,求的相对误差。(函数误差的计算)解:,7计算球的体积,为了使体积的相对误差限为,问度量半径时允许的相对误差限为多大?(函数误差的计算)解:球体积为 ,欲使,必须 。8 设,求证:(1)(2)利用(1)中的公式正向递推计算时误差逐步增大;反向递推计算时误差逐步减小。(计算方法的比较选择)解:如果初始误差为,若是向前递推,有可见,初始误差的绝对值被逐步地扩大了。如果是向后递推,其误差为可见,初始误差的绝对值被逐步减少了。第二章 插值
3、法姓名 学号 班级 习题主要考察点:拉格朗日插值法的构造,均差的计算,牛顿插值和埃尔米特插值构造,插值余项的计算和应用。1 已知,求的拉氏插值多项式。(拉格朗日插值)解法一(待定系数法):设,由插值条件,有解得:。故 。解法二(基函数法):由插值条件,有2 已知,用线性插值求的近似值。(拉格朗日线性插值)解:由插值节点与被插函数,可知,其线性插值函数为的近似值为。3 若为互异节点,且有试证明。(拉格朗日插值基函数的性质)解:考虑辅助函数,其中,。是次数不超过的多项式,在节点()处,有这表明,有n+1个互异实根。故,从而对于任意的均成立。4 已知,用抛物线插值计算的值并估计截断误差。(拉格朗日二
4、次插值)解:由插值条件,其抛物线插值函数为将代入,计算可得:。其余项为: 其中,故误差的上界为:。5 用余弦函数在,三个节点处的值,写出二次拉格朗日插值多项式, 并近似计算及其绝对误差与相对误差,且与误差余项估计值比较。(拉格朗日二次插值)解:由插值条件,二次拉格朗日插值多项式为绝对误差为:相对误差为:余项为:,其中,其余项的上界为:比较可知,实际计算所得的绝对误差较余项公式所估计出的值要小一些。6 已知函数值,求函数的四阶均差和二阶均差。(均差的计算)解:采用列表法来计算各阶均差,有xy一阶均差二阶均差三阶均差四阶均差0611043461814/34823661/362126529/311/
5、151/15从表中可查得:。xy一阶均差二阶均差48211072/3346186故。其实,根据均差的对称性,该值在第一个表中就可以查到。7 设求之值,其中,而节点互异。(均差的计算)解:由均差可以表示成为函数值的线性组合,有而 ,故。8 如下函数值表012419233建立不超过三次的牛顿插值多项式。(牛顿插值多项式的构造)解:先构造均差表xf(x)一阶均差二阶均差三阶均差0119822314343-10-8-11/4故 。9求一个次数小于等于三次多项式,满足如下插值条件:,。(插值多项式的构造)解法一(待定系数法):设,则,由插值条件,有解得:。故 解法二(带重节点的均差法):据插值条件,造差
6、商表xy一阶差商二阶差商三阶差商122422431312852故 10 构造一个三次多项式,使它满足条件(埃尔米特插值)。解:设,利用插值条件,有解得:。11 设。(1)试求在上的三次埃尔米特插值多项式,使得,以升幂形式给出。(2)写出余项的表达式。(埃尔米特插值及其余项的计算)。解:,设,解得:,。故 。,其中,。12 若,试证明:(插值余项的应用)解:以为插值条件,作线性插值多项式,有其余项为故 。13 设求使;又设 ,则估计余项的大小。(插值误差的估计)解:由插值条件,有解得:从而 其余项为第三章 函数逼近姓名 学号 班级 习题主要考察点:最小二乘法,最佳平方逼近,正交多项式的构造。1
7、设,求于上的线性最佳平方逼近多项式。(最佳平方逼近)解:,法方程组为解得:,线性最佳平方逼近多项式为:。2 令,且设,求使得为于 上的最佳平方逼近多项式。(最佳平方逼近)解:,法方程组为解得:,线性最佳平方逼近多项式为:。3证明:切比雪夫多项式序列在区间上带权正交。(正交多项式的证明)解:对于,有对于,有故,序列在-1,1上带权正交。4求矛盾方程组:的最小二乘解。(最小二乘法)解法一:求与,使得达到最小。于是,令即:,其最小二乘解为:。解法二:,记作,该矛盾方程组的最小二乘解,应满足以下方程组,即解之,得。5 已知一组试验数据22.53455.544.5688.59试用直线拟合这组数据. (计
8、算过程保留3位小数)。(最小二乘线性逼近)解:作矩阵,法方程为即解得:,。其直线拟合函数为。6 用最小二乘原理求一个形如的经验公式,使与下列数据相拟合.19253138 441932.34973.397.8(最小二乘二次逼近)解:等价于对数据表3616259611444 19361932.34973.397.8作线性拟合。其法方程组为:解得:,故经验公式为 。第四章 数值积分姓名 学号 班级 习题主要考察点:代数精度的计算,构造插值型求积公式(梯形,辛甫生公式),复化求积的计算,高斯公式的构造。1给定求积公式试确定使它的代数精度尽可能高。(代数精度的应用和计算)解:分别取,使上述数值积分公式准
9、确成立,有;解得:。故求积公式为。再取,左边=,右边=再取,左边=,右边=此求积公式的最高代数精度为3。2 求积公式,试确定系数,及,使该求积公式具有尽可能高的代数精确度,并给出代数精确度的次数。(代数精度的应用和计算)解:分别取,使求积公式准确成立,有解得:。求积公式为。再取,左边=右边故该求积公式的最高代数精度为2。3数值积分公式,是否为插值型求积公式,为什么?又该公式的代数精确度为多少?(插值型求积公式特征)解:令,故代数精度为1。由于求积节点个数为2,代数精度达到1次,故它是插值型的求积公式。4如果,证明用梯形公式计算积分所得到的结果比准确值大,并说明其几何意义。(梯形求积)解:梯形求
10、积公式是由过点,的线性插值函数在a,b上的定积分。注意到:在区间a,b上,而,有从而。其几何意义可作以下解释:在区间a,b上,故曲线下凹,直线位于曲线之上,因此,曲边梯形的面积小于梯形面积。5用的复化梯形公式计算积分,并估计误差。(复化梯形求积)解:,取求积节点为因,则误差大约为:。6设,则用复化辛甫生公式计算,若有常数使 ,则估计复化辛甫生公式的整体截断误差限。(复化辛甫生公式)解:7已知高斯求积公式 将区间0,1二等分,用复化高斯求积法求定积分的近似值。(高斯公式)解:对于作变量换,有对于作变量换,有8 试确定常数A,B,C和,使得数值积分公式有尽可能高的代数精度。试问所得的数值积分公式代
11、数精度是多少?它是否为高斯型的?(代数精度的应用和计算,高斯点的特征)解:分别取,使上述数值积分公式准确成立,有;整理得:解得:。数值求积公式为再取,左边=,右边=再取,左边=,右边=可见,该数值求积公式的最高代数精度为5。由于该公式中的节点个数为3,其代数精度达到了次,故它是高斯型的。9设是0,1区间上带权的最高次幂项系数为1的正交多项式系(1)求。(2)构造如下的高斯型求积公式。(高斯求积)解(1):采用施密特正交化方法,来构造带权且在0,1上正交的多项式序列取,设,且它与在0,1上带权正交,于是,故 。设,且它与、在0,1上带权正交,于是,解(2):的零点为:。设 分别取,使上述求积公式
12、准确成立,有,即解得:,。高斯型求积公式为第五章 非线性方程求根姓名 学号 班级 习题主要考察点:二分法、迭代法、牛顿法和弦截法求根,迭代法求根的收敛性和收敛速度的讨论。1用二分法求方程的正根,要求误差小于0.05。(二分法)解:,在0,2连续,故0,2为函数的有根区间。(1)计算,故有根区间为1,2。(2)计算,故有根区间为。(3)计算,故有根区间为。(4)计算,故有根区间为。(5)计算,故有根区间为。(6)计算,故有根区间为。(7)计算,故有根区间为。(8)若取中点作为取根的近似值,其误差小于取近似根,可满足精度要求。2说明方程 在区间1,2内有惟一根,并选用适当的迭代法求(精确至3位有效
13、数),并说明所用的迭代格式是收敛的。(迭代法)解: ,故函数单调增加,因此,该方程在(1,2)之间存在着惟一的实根。取迭代函数 显然,且故迭代 ()对任意初始值收敛。对于初值,其迭代值分别为,由于,故作为近似值,已精确到了3位有效数字。3设有解方程的迭代法 (1)证明均有(为方程的根)。(2)此迭代法的收敛阶是多少,证明你的结论。 (3) 取用此迭代法求方程根的近似值,误差不超过,列出各次迭代值。(和收敛性讨论)解(1):,(),故该迭代对任意初值均收敛于方程的根。解(2):由,故有。,故该迭代的收敛速度是1阶的。解(3):取,代入迭代式,可计算出以下结果:,由于,取可满足精度要求。4设,试证
14、明:由 ,得到的序列收敛于。(收敛性证明)证明:由知,方程有根。由,当时,有,即序列收敛于。5 设方程在0,1内的根为,若采用迭代公式,试证明:均有为方程的根);此迭代的收敛阶是多少,证明你的结论。(迭代法和收敛性讨论)解:迭代函数,当故迭代在区间上整体收敛。设,则,且故 故该迭代的收敛速度为1阶的。6 方程在附近有根,把方程写成3种不同的等价形式:(1) ,对应迭代格式:(2) ,对应迭代格式:(3) ,对应迭代格式:讨论这些迭代格式在时的收敛性。若迭代收敛,试估计其收敛速度,选一种收敛格式计算出附近的根到4位有效数字。(收敛速度的计算和比较)解:,故方程在上有根。,故方程在上有根。,故方程
15、在上有根。对于迭代式(1):, 而,故该迭代局部收敛,且收敛速度为1阶的。对于迭代式(2):在上,又,故该迭代在上整体收敛,且收敛速度为一阶的。对于迭代式(3):在1,2上的值域为,该迭代式不收敛。取迭代式,进行计算,其结果如下:,取为近似值具有4位有效数字。7设 (1) 写出解 的牛顿迭代格式;(2) 证明此迭代格式是线性收敛的。(牛顿迭代的构造与收敛速度)解:牛顿迭代式为 ,方程的根为,因,故迭代局部收敛。又因,故迭代收敛速度为1阶。8 设计一个计算的牛顿迭代法,且不用除法(其中)。(牛顿迭代法)解:考虑方程,而,该迭代局部收敛。9 用牛顿法求的近似值,取或11为初始值,计算过程保留4位小
16、数。(牛顿迭代的构造)解:考虑方程,取为初始值,计算其迭代值如下:,取为初始值,计算其迭代值如下:,10设是非线性方程的m重根,试证明:迭代法具有至少2阶的收敛速度。(收敛速度证明)解:设是非线性方程的m重根,则,且及,其牛顿迭代函数为牛顿迭代式故该迭代的收敛速度至少是2阶的。11设是非线性方程的m重根,证明:用牛顿迭代法求只是线性收敛。(收敛速度证明)解:设是非线性方程的m重根,则,且及,其牛顿迭代函数为牛顿迭代式故收敛速度为1阶的。12设,在附近有直到阶的连续导数,且,试证:迭代法在附近是阶收敛的。(收敛速度证明)解:将在点附近作泰勒展式,有,其中,在与之间。于是: ,其中,在与之间。由于
17、,故,从而。因此,迭代的收敛速度为p。第六章 常微分方程数值解姓名 学号 班级 习题主要考察点:欧拉方法的构造,单步法的收敛性和稳定性的讨论,线性多步法中亚当姆斯方法的构造和讨论。1 用改进的欧拉公式,求以下微分方程的数值解(取步长),并与精确解作比较。(改进的尤拉公式的应用)解:原方程可转化为 ,令,有解此一阶线性微分方程,可得 。利用以下公式求在节点处的数值解,其中,初值为。MATLAB程序如下:x(1)=0;%初值节点y(1)=1;%初值fprintf(x(%d)=%f,y(%d)=%f,yy(%d)=%fn,1,x(1),1,y(1),1,y(1);for i=1:5 yp=y(i)+
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数值 分析 习题 答案 43
限制150内