数值分析第六章数值积分
- 格式:ppt
- 大小:2.32 MB
- 文档页数:37
数值分析原理习题答案数值分析原理习题答案【篇一:数值分析习题】学号班级习题主要考察点:有效数字的计算、计算方法的比较选择、误差和误差限的计算。
1 若误差限为0.5?10,那么近似数0.003400有几位有效数字?(有效数字的计算)2 ??3.14159?具有4位有效数字的近似值是多少?(有效数字的计算)3 已知a?1.2031,b?0.978是经过四舍五入后得到的近似值,问a?b,a?b有几位有效数字?(有效数字的计算)4 设x?0,x的相对误差为?,求lnx的误差和相对误差?(误差的计算)**5测得某圆柱体高度h的值为h?20cm,底面半径r的值为r?5cm,已知5|h?h*|?0.2cm,|r?r*|?0.1cm,求圆柱体体积v??rh的绝对误差限与相对误差限。
(误差限的计算)6 设x的相对误差为a%,求y?xn的相对误差。
(函数误差的计算)7计算球的体积,为了使体积的相对误差限为1%,问度量半径r时允许的相对误差限为多大?(函数误差的计算)128 设in?e1nxx?edx,求证: 0(1)in?1?nin?1(n?0,1,2?)(2)利用(1)中的公式正向递推计算时误差逐步增大;反向递推计算时误差逐步减小。
(计算方法的比较选择)第二章插值法姓名学号班级习题主要考察点:拉格朗日插值法的构造,均差的计算,牛顿插值和埃尔米特插值构造,插值余项的计算和应用。
1 已知f(?1)?2,f(1)?1,f(2)?1,求f(x)的拉氏插值多项式。
(拉格朗日插值)2 已知y?x,x0?4,x1?9,用线性插值求7的近似值。
(拉格朗日线性插值)3 若xj(j?0,1,...n)为互异节点,且有lj(x)?试证明(x?x0)(x?x1)?(x?xj?1)(x?xj?1)?(x?xn)(xj?x0)(xj?x1)?(xj?xj?1)(xj?xj?1)?(xj?xn)xlj?0nkjj(拉格朗日插值基函数的性质) (x)?xk(k?0,1,...n)。
第1篇一、实验目的本次实验旨在通过数值分析的方法,研究几种常见的数值积分方法,包括梯形法、辛普森法、复化梯形法和龙贝格法,并比较它们在计算精度和效率上的差异。
通过实验,加深对数值积分理论和方法的理解,提高编程能力和实际问题解决能力。
二、实验内容1. 梯形法梯形法是一种基本的数值积分方法,通过将积分区间分割成若干个梯形,计算梯形面积之和来近似积分值。
实验中,我们选取了几个不同的函数,对积分区间进行划分,计算积分近似值,并与实际积分值进行比较。
2. 辛普森法辛普森法是另一种常见的数值积分方法,它通过将积分区间分割成若干个等距的区间,在每个区间上使用二次多项式进行插值,然后计算多项式与x轴围成的面积之和来近似积分值。
实验中,我们对比了辛普森法和梯形法的计算结果,分析了它们的精度差异。
3. 复化梯形法复化梯形法是对梯形法的一种改进,通过将积分区间分割成多个小区间,在每个小区间上使用梯形法进行积分,然后计算所有小区间积分值的和来近似积分值。
实验中,我们对比了复化梯形法和辛普森法的计算结果,分析了它们的精度和效率。
4. 龙贝格法龙贝格法是一种通过外推加速提高计算精度的数值积分方法。
它通过比较使用不同点数(n和2n)的积分结果,得到更高精度的积分结果。
实验中,我们使用龙贝格法对几个函数进行积分,并与其他方法进行了比较。
三、实验步骤1. 编写程序实现梯形法、辛普森法、复化梯形法和龙贝格法。
2. 选取几个不同的函数,对积分区间进行划分。
3. 使用不同方法计算积分近似值,并与实际积分值进行比较。
4. 分析不同方法的精度和效率。
四、实验结果与分析1. 梯形法梯形法在计算精度上相对较低,但当积分区间划分足够细时,其计算结果可以接近实际积分值。
2. 辛普森法辛普森法在计算精度上优于梯形法,但当积分区间划分较细时,计算量较大。
3. 复化梯形法复化梯形法在计算精度上与辛普森法相当,但计算量较小。
4. 龙贝格法龙贝格法在计算精度上优于复化梯形法,且计算量相对较小。
实验09 数值微积分与方程数值求解(第6章 MATLAB 数值计算)一、实验目的二、实验内容1. 求函数在指定点的数值导数232()123,1,2,3026x x x f x x xx x==2. 用数值方法求定积分(1) 210I π=⎰的近似值。
程序及运行结果:《数学软件》课内实验王平(2) 2221I dx x π=+⎰程序及运行结果:3. 分别用3种不同的数值方法解线性方程组6525494133422139211x y z u x y z u x y z u x y u +-+=-⎧⎪-+-=⎪⎨++-=⎪⎪-+=⎩ 程序及运行结果:4. 求非齐次线性方程组的通解1234123412342736352249472x x x x x x x x x x x x +++=⎧⎪+++=⎨⎪+++=⎩5. 求代数方程的数值解(1) 3x +sin x -e x =0在x 0=1.5附近的根。
程序及运行结果(提示:要用教材中的函数程序line_solution ):(2) 在给定的初值x 0=1,y 0=1,z 0=1下,求方程组的数值解。
23sin ln 70321050y x y z x z x y z ⎧++-=⎪+-+=⎨⎪++-=⎩6. 求函数在指定区间的极值(1) 3cos log ()xx x x xf x e ++=在(0,1)内的最小值。
(2) 33212112122(,)2410f x x x x x x x x =+-+在[0,0]附近的最小值点和最小值。
7. 求微分方程的数值解,并绘制解的曲线2250(0)0'(0)0xd y dyy dx dx y y ⎧-+=⎪⎪⎪=⎨⎪=⎪⎪⎩程序及运行结果(注意:参数中不能取0,用足够小的正数代替):令y 2=y,y 1=y ',将二阶方程转化为一阶方程组:'112'211251(0)0,(0)0y y y x x y y y y ⎧=-⎪⎪=⎨⎪==⎪⎩8. 求微分方程组的数值解,并绘制解的曲线123213312123'''0.51(0)0,(0)1,(0)1y y y y y y y y y y y y =⎧⎪=-⎪⎨=-⎪⎪===⎩程序及运行结果:三、实验提示四、教程:第6章 MATLAB 数值计算(2/2)6.2 数值微积分 p155 6.2.1 数值微分1. 数值差分与差商对任意函数f(x),假设h>0。
第一章 非线性方程和方程组的数值解法 1)二分法的基本原理,误差:~12k b ax α+--<2)迭代法收敛阶:1lim0i pi ic εε+→∞=≠,若1p =则要求01c <<3)单点迭代收敛定理:定理一:若当[],x a b ∈时,[](),x a b ϕ∈且'()1x l ϕ≤<,[],x a b ∀∈,则迭代格式收敛于唯一的根;定理二:设()x ϕ满足:①[],x a b ∈时,[](),x a b ϕ∈, ②[]121212,,, ()(),01x x a b x x l x x l ϕϕ∀∈-≤-<<有 则对任意初值[]0,x a b ∈迭代收敛,且:110111i i iii x x x ll x x x lαα+-≤---≤--定理三:设()x ϕ在α的邻域内具有连续的一阶导数,且'()1ϕα<,则迭代格式具有局部收敛性;定理四:假设()x ϕ在根α的邻域内充分可导,则迭代格式1()i i x x ϕ+=是P 阶收敛的 ()()()0,1,,1,()0j P j P ϕαϕα==-≠L (Taylor 展开证明)4)Newton 迭代法:1'()()i i i i f x x x f x +=-,平方收敛 5)Newton 迭代法收敛定理:设()f x 在有根区间[],a b 上有二阶导数,且满足: ①:()()0f a f b <; ②:[]'()0,,f x x a b ≠∈;③:[]'',,f x a b ∈不变号④:初值[]0,x a b ∈使得''()()0f x f x <;则Newton 迭代法收敛于根α。
6)多点迭代法:1111111()()()()()()()()()i i i i i i i i i i i i i i i f x f x f x x x x x f x f x f x f x f x f x x x -+-----=-=+----收敛阶:P =7)Newton 迭代法求重根(收敛仍为线性收敛),对Newton 法进行修改 ①:已知根的重数r ,1'()()i i i i f x x x rf x +=-(平方收敛) ②:未知根的重数:1''()(),()()()i i i i u x f x x x u x u x f x +=-=,α为()f x 的重根,则α为()u x 的单根。
数值分析智慧树知到课后章节答案2023年下湖南师范大学第一章测试1.在数值计算中因四舍五入产生的误差称为()A:观测误差 B:方法误差 C:舍入误差 D:模型误差答案:舍入误差2.当今科学活动的三大方法为()。
A:科学计算 B:实验C:数学建模 D:理论答案:科学计算;实验;理论3.计算过程中如果不注意误差分析,可能引起计算严重失真。
A:错 B:对答案:对4.算法设计时应注意算法的稳定性分析。
A:对 B:错答案:对5.在进行数值计算时,每一步计算所产生的误差都是可以准确追踪的。
A:错 B:对答案:错第二章测试1.A: B: C: D:答案:2.某函数过(0,1),(1,2)两点,则其关于这两点的一阶差商为A:3 B:0 C:2 D:1 答案:13.A: B: C: D:答案:4.下列说法不正确的是A:高次多项式插值不具有病态性质 B:分段线性插值逼近效果依赖于小区间的长度 C:分段线性插值的导数一般不连续D:分段线性插值的几何图形就是将插值点用折线段依次连接起来答案:分段线性插值的几何图形就是将插值点用折线段依次连接起来5.下列关于分段线性插值函数的说法,正确的是A:对于光滑性不好的函数优先用分段线性插值 B:对于光滑性较好的函数优先用分段线性插值 C:一次函数的分段线性插值函数是该一次函数本身 D:二次函数的分段线性插值函数是该二次函数本身答案:对于光滑性不好的函数优先用分段线性插值;一次函数的分段线性插值函数是该一次函数本身6.A: B: C:D:答案:;;7.同一个函数基于同一组插值节点的牛顿插值函数和拉格朗日插值函数等价。
A:错 B:对答案:对第三章测试1.A: B:C:D:答案:2.以下哪项是最佳平方逼近函数的平方误差A: B: C:D:答案:3.当区间为[-1,1],Legendre多项式族带权 ( ) 正交。
A: B: C: D:答案: 4.n次Chebyshev多项式在 (-1,1) 内互异实根的个数为A:n+1 B:n-1 C:nD:n+2 答案:n5.用正交函数族做最小二乘法有什么优点A:每当逼近次数增加1时,系数需要重新计算 B:得到的法方程非病态C:不用解线性方程组,系数可简单算出 D:每当逼近次数增加1时,之前得到的系数不需要重新计算答案:得到的法方程非病态;不用解线性方程组,系数可简单算出;每当逼近次数增加1时,之前得到的系数不需要重新计算6.用正交多项式作基求最佳平方逼近多项式,当n较大时,系数矩阵高度病态,舍入误差很大。
学科分类号110.3420本科毕业论文题目几种常用数值积分方法的比较姓名潘晓祥学号1006020540200院(系)数学与计算机科学学院专业数学与应用数学年级2010 级指导教师雍进军职称讲师二〇一四年五月贵州师范学院本科毕业论文(设计)诚信声明本人郑重声明:所呈交的本科毕业论文(设计),是本人在指导老师的指导下,独立进行研究工作所取得的成果,成果不存在知识产权争议,除文中已经注明引用的内容外,本论文不含任何其他个人或集体已经发表或撰写过的作品成果。
对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。
本人完全意识到本声明的法律结果由本人承担。
本科毕业论文作者签名:年月日贵州师范学院本科毕业论文(设计)任务书毕业设计题目几种常用数值积分方法的比较作者姓名潘晓祥学号1006020540200 年级2010级所属学院数学与计算机科专业数学与应用数学班级四班指导教师签名雍进军讲师职称讲师开题日期2013年7月10日主要目标1.了解什么数值积分基本思想和一些常用的数值积分方法;2.对各种数值积分方法的误差以及代数精度进行分析;3.对各积分方法进行比较总结出优缺点。
主要要求通过对几种常用的数值积分方法进行了的分析,并用这几种方法对被积函数是普通函数做了数值积分,并在计算机上进行实验。
数值积分是计算方法或数值分析理论中非常重要的内容,数值积分方法也是解决实际计算问题的重要方法,对几种常用数值积分方法的分析很必要。
主要内容本文通过对复化求积公式, Newton—Cotes求积公式, Romberg求积公式,高斯型求积公式进行分析讨论并在计算机上积分实验,从代数精度,求积公式误差等角度对这些方法进行分析比较,并总结出每种求积分法的优缺点以及实用性。
贵州师范学院本科毕业论文(设计)开题报告书论文题目几种常用数值积分方法的比较作者姓名潘晓祥学号1006020540200 年级2010级数学与计算机所属学院专业数学与应用数学班级数本(4)班科学学院指导教师姓名雍进军职称讲师预计字数5000.00字题目性质应用研究日期2013年7月05 日选题的原由:研究意义:数值积分是数学上的重要课题之一,是数值分析中的重要内容之一,也是数学的研究重点.并在实际问题及应用中有着广泛的应用.常用于科学与工程的计算中,如涉及到积分方程,工程计算,计算机图形学,金融数学等应用科学领域都有着相当重要的应用,所以研究数值积分问题有很重要的意义.数值积分是研究如何求出一个积分的数值.这一课题的起源可追溯到古代,其中一个突出的例子是希腊人用内接与外接正多边形推算出圆面积的方法.也正是此法使阿基米德得以求出π值得上界与下界,若干世纪以来,尤其是十六世纪后,已提出了多种数值积分方法,其中有矩形求积法,内插求积法,牛顿科特斯公式,复化求积公式,龙贝格求积公式,高斯型求积公式.但各种方法都有特点,在不同的情况下试用程度不同,我们将着重从求积公式的代数精度和余项等角度对这些方法进行分析比较. 研究动态:这些年来,有关数值积分的研究已经成为一个很活跃的研究领域,历史上,阿基米德,牛顿,欧拉,高斯,切比雪夫等人都对此有过贡献.研究出各种各样的数值求积公式,但一个好的数值求积公式应该满足:计算简单,误差小,代数精度高.我们将对矩形求积法,内插求积法,牛顿科特斯公式,化求积公式,贝格求积公式,斯型求积公式进行比较.对数值求积公式能有进一步的了解和学习.主要内容:1 数值积分方法的基本思想2 几类常用数值积分方法的基本分析2.1 Newton—Cotes求积公式2.2 复化求积公式2.3 Romberg求积公式2.4 高斯型求积公式3 几类数值积分方法的简单比较评述4利用MATLAB编程应用对几类求积算法的分析比较研究方法:本论文主要通过对相关文献和书籍的参考,合自己的见解,复化求积公式,Newton—Cotes求积公式,Romberg求积公式,高斯型求积公式进行讨论并进行上机实验,从代数精度,求积公式误差等角度对这些方法进行分析比较.完成期限和采取的主要措施:本论文计划用6个月的时间完成,阶段的任务如下:(1)7月份查阅相关书籍和文献;(2)8月份完成开题报告并交老师批阅;(3)9月份完成论文初稿并交老师批阅;(4)10月份完成论文二搞并交老师批阅;(5)11月份完成论文三搞;(6)12月份定稿.主要措施:考相关书籍和文献,合自己的见解,老师的指导下和同学的帮助下完成主要参考文献及资料名称:[1] 关治. 陆金甫. 数学分析基础(第二版)[M]. 北京:等教育出版社.2010.7[2] 胡祖炽. 林源渠. 数值分析[M] 北京:等教育出版社.1986.3[3] 薛毅. 数学分析与实验[M] 北京:业大学出版社2005.3[4] 徐士良. 数值分析与算法[M]. 北京:械工业出版社2007.1[5] 王开荣. 杨大地. 应用数值分析[M] 北京:等教育出版社2010.7[6] 杨一都. 数值计算方法[M]. 北京:等教育出版社 . 2008.4[7] 韩明. 王家宝. 李林. 数学实验(MATLAB)版[M]. 上海:济大学出版社2012.1[8] 圣宝建. 关于数值积分若干问题的研究[J]. 南京信息工程大学. 2009.05.01. : 42[9] 刘绪军. 几种求积公式计算精确度的比较[J]. 南京职业技术学院. 2009.[10] 史万明.吴裕树.孙新.数值分析[M]. 北京理工大学出版社.2010.4.开题报告会纪要时间2013年8月26日地点宁静楼229教师办公室与会人员姓名职务(职称)姓名职务(职称)姓名职务(职称)雍进军导师(讲师)邓喜才副教授李晟副教授龙林林组长指导教师意见:签名:年月日会议记录摘要:指导小组针对课题《二次函数性质的应用》提问了以下问题以及报告人的回答:雍老师问:选择此题目的目的?潘晓祥答:随着计算机和计算方法的飞速发展,几乎所有学科都走向定量化和精确化,计算数学中的数值计算方法则是解决“计算”问题的桥梁和工具。
重点考察内容第一章:基本概念第二章:Gauss消去法,Lu分解法第三章:题型:具体题+证明,误差分析三个主要迭代法,条件误差估计,范数的小证明第四章:掌握三种插值方法:拉格朗日,牛顿,厄尔米特,误差简单证明,构造复合函数第五章:最小二乘法计算第六章:梯形公式,辛普森(抛物线)公式,高斯公式三个重要公式,误差分析高斯求积公式的构造第七章:几种常用的迭代格式构造,收敛性证明第九章:基本概念(收敛阶,收敛条件,收敛区域等)简单欧拉法第一章误差1. 科学计算中的误差来源有4个,分别是 _________ , ________ , ________ , ________ 。
2. 用Taylor 展开近似计算函数f (x ) :、f (x 0) f'(x 0)(x-x 0),这里产生是什么误差?3. 0.7499作3的近似值,是位有效数字,65.380是舍入得到的近似值,有 几4位有效数字,相对误差限为 _______ . 0.0032581是四舍五入得到的近似值,有 ________ 位有效数字.4. 改变下列表达式,使计算结果比较精确:(1) —|x|=1( 2) +J 1-丄,|x|=11 +2x 1 +x Y x Y x1「cosx(3), x=0,|x| 1. (4) sin : -sin :, 一—■x5. 采用下列各式计算(、、2-1)6时,哪个计算效果最好?并说明理由。
1 1(1) 6 ( 2) 99-70,2( 3) (3-2、月)6( 4) 3(V2+1)6(3 + 2问36. 已知近似数x *有4位有效数字,求其相对误差限。
上机实验题:kx匸 Xx1、 利用Taylor 展开公式计算 e,编一段小程序,上机用单精度计算 e 的函数k£k !值.分别取x =1, 5, 10, 20, -1,-5,-10,-15,-20,观察所得结果是否合理,如不合 理请分析原因并给出解决方法.1 n2、 已知定积分I n— dx,n =0,1,2,…,20,有如下的递推关系 ‘° x +6可建立两种等价的计算公式11(1) I n 61 nd ,取 I 。