求函数近似解的一种计算方法
§2.4.2 求函数的零点近似解的一种计算方法——二分法
教学目标:
知识与技能:通过具体实例理解二分法的概念及其适用条件,了解二分法是求方程近似解的常用方法,从中体会函数与方程之间的联系及其在实际问题中的应用.过程与方法:能借助计算器用二分法求方程的近似解,并了解这一数学思想,为学习算法做准备.
情感、态度、价值观:体会数学逼近过程,感受精确与近似的相对统一.
教学重点、难点:
重点通过用二分法求方程的近似解,体会函数的零点与方程根之间的联系,初步形成用函数观点处理问题的意识.
难点恰当地使用信息技术工具,利用二分法求给定精确度的方程的近似解.
教学程序与环节设计:
由二分查找及高次多项式方程的求问题引入.
初步应用二分法解
.二分法为什么可以逼近零点的再分析;
.追寻阿贝尔和伽罗瓦.
数值计算方法试题及答案
【 数值计算方法试题一 一、 填空题(每空1分,共17分) 1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数,需对分( )次。 2、迭代格式)2(2 1-+=+k k k x x x α局部收敛的充分条件是α取值在( )。 3、已知?????≤≤+-+-+-≤≤=31)1()1()1(211 0)(2 33x c x b x a x x x x S 是三次样条函数, 则 a =( ), b =( ), c =( )。 4、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数,则 ∑== n k k x l 0)(( ), ∑== n k k j k x l x 0 )(( ),当2≥n 时 = ++∑=)()3(20 4x l x x k k n k k ( )。 ; 5、设1326)(2 47+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 和=?07 f 。 6、5个节点的牛顿-柯特斯求积公式的代数精度为 ,5个节点的求积公式最高代数精度为 。 7、{}∞ =0)(k k x ?是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族,其中1)(0=x ?,则?= 1 4)(dx x x ? 。 8、给定方程组?? ?=+-=-2211 21b x ax b ax x ,a 为实数,当a 满足 ,且20<<ω时,SOR 迭代法收敛。 9、解初值问题 00 (,)()y f x y y x y '=?? =?的改进欧拉法 ??? ??++=+=++++)],(),([2),(] 0[111] 0[1n n n n n n n n n n y x f y x f h y y y x hf y y 是 阶方法。
框架结构内力与位移计算
《高层建筑结构与抗震》辅导材料四 框架结构内力与位移计算 学习目标 1、熟悉框架结构在竖向荷载和水平荷载作用下的弯矩图形、剪力图形和轴力图形; 2、熟悉框架结构内力与位移计算的简化假定及计算简图的确定; 3、掌握竖向荷载作用下框架内力的计算方法——分层法; 4、掌握水平荷载作用下框架内力的计算方法——反弯点法和D值法,掌握框架结构的侧移计算方法。 学习重点 1、竖向荷载作用下框架结构的内力计算; 2、水平荷载作用下框架结构的内力及侧移计算。 框架在结构力学中称为刚架,刚架的内力和位移计算方法很多,可分为精确算法和近似算法。精确法是采用较少的计算假定,较为接近实际情况地考虑建筑结构的内力、位移和外荷载的关系,一般需建立大型的代数方程组,并用电子计算机求解;近似算法对建筑结构引入较多的假定,进行简化计算。由于近似计算简单、易于掌握,又能反映刚架受力和变形的基本特点,因此近似的计算方法仍为工程师们所常用。 本章内容主要介绍框架结构在荷载作用下内力与位移的近似计算方法。其中分层法用于框架结构在竖向荷载作用下的内力计算,反弯点法和D值法用于框架结构在水平荷载作用下的内力计算。既然是近似计算,就需要熟悉框架结构的计算简图和各种计算方法的简化假定。 一、框架结构计算简图的确定 一般情况下,框架结构是一个空间受力体系,可以按照第四章所述的平面结构假定的简化原则,忽略结构纵向和横向之间的空间联系,忽略各构件的抗扭作用,将框架结构简化为沿横方向和纵方向的平面框架,承受竖向荷载和水平荷载,进行内力和位移计算。 结构设计时一般取中间有代表性的一榀横向框架进行分析,若作用于纵向框架上的荷载各不相同,则必要时应分别进行计算。 框架结构的节点一般总是三向受力的,但当按平面框架进行结构分析时,则节点也相应地简化。在常见的现浇钢筋混凝土结构中,梁和柱内的纵向受力钢筋都将穿过节点或锚入节点区,这时节点应简化为刚
计算方法的课后答案
《计算方法》习题答案 第一章 数值计算中的误差 1.什么是计算方法?(狭义解释) 答:计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算,以便在计算机上编程上机,求出问题的数值解,并对算法的收敛性、稳定性和误差进行分析、计算。 2.一个实际问题利用计算机解决所采取的五个步骤是什么? 答:一个实际问题当利用计算机来解决时,应采取以下五个步骤: 实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果 4.利用秦九韶算法计算多项式4)(5 3 -+-=x x x x P 在3-=x 处的值,并编程获得解。 解:400)(2 3 4 5 -+?+-?+=x x x x x x P ,从而 所以,多项式4)(5 3 -+-=x x x x P 在3-=x 处的值223)3(-=-P 。 5.叙述误差的种类及来源。 答:误差的种类及来源有如下四个方面: (1)模型误差:数学模型是对实际问题进行抽象,忽略一些次要因素简化得到的,它是原始问题的近似,即使数学模型能求出准确解,也与实际问题的真解不同,我们把数学模型与实际问题之间存在的误差称为模型误差。 (2)观测误差:在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的,由于仪器的精密性,实验手段的局限性,周围环境的变化以及人们的工作态度和能力等因素,而使数据必然带有误差,这种误差称为观测误差。 (3)截断误差:理论上的精确值往往要求用无限次的运算才能得到,而实际运算时只能用有限次运算的结果来近似,这样引起的误差称为截断误差(或方法误差)。 (4)舍入误差:在数值计算过程中还会用到一些无穷小数,而计算机受机器字长的限制,它所能表示的数据只能是一定的有限数位,需要把数据按四舍五入成一定位数的近似的有理数来代替。这样引起的误差称为舍入误差。 6.掌握绝对误差(限)和相对误差(限)的定义公式。 答:设* x 是某个量的精确值,x 是其近似值,则称差x x e -=* 为近似值x 的绝对误差(简称误差)。若存在一个正数ε使ε≤-=x x e * ,称这个数ε为近似值x 的绝对误差限(简称误差限或精度)。 把绝对误差e 与精确值* x 之比* **x x x x e e r -==称为近似值x 的相对误差,称
框剪结构剪力墙合理数量的确定
RC 框—剪结构剪力墙合理数量的确定 陈烈火 杨建明 (厦门中福元建筑设计研究院,厦门,361009) 提 要:剪力墙是框一剪结构的主要的抗侧力构件.本文根据框架与剪力墙协调作用、层间位移的限值推导出框一剪结构剪力墙数量的简化确定方法,以供工程设计时参考。 关键词:RC 框-剪结构、剪力墙数量、协调工作变形、地震力、风荷载。 一、前言 RC 框剪结构是高层最常见的结构体系之一,整体变形呈剪弯型;就水平剪力而言剪力墙承担了几乎80%的水平剪力;就顶点位移而言,框剪结构刚度特征值λ=1.5~2.0时,水平荷载作用下的结构顶点位移约为对应结构中的剪力墙单独承受水平荷载作用下的位移的 40~55%[1] 。而相同受力条件的框剪结构即使剪力墙刚度增加1倍其顶点位移比值及层间位 移比值的减少也仅13~19%;而剪力墙刚度增加1倍,地震力增大20%左右[2] 。剪力墙多了则不经济,少了则影响结构安全,因此剪力墙数量的确定是框剪结构设计的重要环节。 以往的剪力墙数量的估算大多数以剪力墙的顶点位移或剪力墙所承担的首层剪力来估算框剪结构的剪力墙数量;作为剪弯型变形的框剪结构的最大层间变形往往出现在(1/3~2/3)H (H 为结构总高度)范围的楼层。因此本文引入协调变形工作原理,取最接近最大层间变形处(0.5~0.6H )的楼层的框架抗侧刚度作为框架平均抗侧刚度估算综合框架抗推刚度参与整体协调变形计算,再由变形限值条件反算剪力墙的数量,同时引入风荷载和地震作用的首层剪力比对估算得到的剪力墙数量进行调整;这样使得估算值与计算机计算值更为接近,以达到更为经济安全的目的。下面我们逐步介绍计算过程。 二、估算框架综合抗推刚度C f : 取0.5~0.6H (H 为结构总高度)处的楼层框架抗侧刚度为D ;在计算调整系数α时假定层高为h ,梁高/12()b h l l =为框架梁柱网跨度、梁宽/2.5b b b h =,A γ为楼层面容重,G E 为总重力荷载代表值,G E ≥H/2 为0.5H 以上重力荷载代表值,柱截面面积Ac ,且柱子为方柱, 柱截面边长为a ,柱轴压比限值为μc =0.85,;楼层建筑建筑面积为Az 。以下是公式的推导: 梁线刚度:33//12 1.607b b b b b b b K E I l E b h l E l -6 ?10=== 由轴压比控制的柱截面尺寸: a =柱线刚度:443//12/34.68c c c c c K E I h E a h l E h 22A γH === 系数K : 333442 1.60734.68 1.11522b b b c b c c c K K E l E h E E h K K l l -6-4A A ??10??10K γH γH ====∑ 中柱抗侧刚度调整系数α: 3 3 1.11522 1.115b c b c K E E h K l E E h -422-4A ?10α?γH ?10==++ 由于强柱弱梁时,K 值一般在0.3~0.5之间,α=0.13~0.2 ;设梁柱刚度比为i ,则取
数值计算方法试题及答案
数值计算方法试题一 一、填空题(每空1分,共17分) 1、如果用二分法求方程在区间内的根精确到三位小数,需对分()次。 2、迭代格式局部收敛的充分条件是取值在()。 3、已知是三次样条函数,则 =( ),=(),=()。 4、是以整数点为节点的Lagrange插值基函数,则 ( ),( ),当时( )。 5、设和节点则 和。 6、5个节点的牛顿-柯特斯求积公式的代数精度为,5个节点的求积公式最高代数精度为。 7、是区间上权函数的最高项系数为1的正交多项式族,其中,则。 8、给定方程组,为实数,当满足,且时,SOR迭代法收敛。 9、解初值问题的改进欧拉法是 阶方法。 10、设,当()时,必有分解式,其中为下三角阵,当其对角线元素满足()条件时,这种分解是唯一的。 二、二、选择题(每题2分) 1、解方程组的简单迭代格式收敛的充要条件是()。(1), (2) , (3) , (4) 2、在牛顿-柯特斯求积公式:中,当系数是负值时,公式的稳定性不能保证,所以实际应用中,当()时的牛顿-柯特斯求积公式不使用。 (1),(2),(3),(4), (1)二次;(2)三次;(3)四次;(4)五次 4、若用二阶中点公式求解初值问题,试问为保证该公式绝对稳定,步长的取值范围为()。 (1), (2), (3), (4)
三、1、 2、(15 (1)(1) 试用余项估计其误差。 (2)用的复化梯形公式(或复化 Simpson公式)计算出该积分的近似值。 四、1、(15分)方程在附近有根,把方程写成三种不同的等价形式(1)对应迭代格式;(2)对应迭代格式;(3)对应迭代格式。判断迭代格式在的收敛性,选一种收敛格式计算附近的根,精确到小数点后第三位。选一种迭代格式建立Steffensen迭代法,并进行计算与前一种结果比较,说明是否有加速效果。 2、(8分)已知方程组,其中 , (1)(1)列出Jacobi迭代法和Gauss-Seidel迭代法的分量形式。 (2)(2)求出Jacobi迭代矩阵的谱半径,写出SOR 迭代法。 五、1、(15分)取步长,求解初值问题用改进的欧拉法求的值;用经典的四阶龙格—库塔法求的值。 2、(8分)求一次数不高于4次的多项式使它满足 ,,,, 六、(下列2题任选一题,4分) 1、1、数值积分公式形如 (1)(1)试确定参数使公式代数精度尽量高;(2)设,推导余项公式,并估计误差。 2、2、用二步法 求解常微分方程的初值问题时,如何选择参数使方法阶数尽可能高,并求局部截断误差主项,此时该方法是几阶的。 数值计算方法试题二 一、判断题:(共16分,每小题2分) 1、若是阶非奇异阵,则必存在单位下三角阵和上三角阵,使唯一成立。()
建筑结构抗震试卷答案
一、填空题:(每空1 分,共20 分) 1、高层建筑结构平面不规则分为、、几种 类型。 2、高层建筑按不同情况分别采用、和的地 震作用计算方法。 3、框架结构水平荷载作用近似手算方法包括、 。 4、高层建筑框架结构柱反弯点高度应考虑、 、的影响。 5、高层建筑剪力墙分为、、 、四种类型。 6、高层建筑框剪结构中的梁分为、 、三种类型。 7、高层建筑框剪结构剪力墙布置应遵循、分散、对 称的原则。 1、扭转不规则、凹凸不规则、楼板局部不连续。 2、底部剪力法、振型分解反应谱法、弹性时程分析法。 3、反弯点法、D 值法。 4、梁柱线刚度比及层数、层次,上、下横梁线刚度比,层高变化。 5、整体剪力墙、整体小开口墙、联肢剪力墙、壁式框架。 6、普通框架梁、剪力墙间连梁、墙技一框架柱连梁。 7、均匀、周边 二、单项选择题:(每题2 分,共10 分) 1、高层建筑采用()限制来保证结构满足舒适度要求。 A、层间位移 B、顶点最大位移 C、最大位移与层高比值 D、顶点最大加速度 2、地震时必须维持正常使用和救灾需要的高层建筑物为()。 A、一类建筑 B、甲类建筑 C、二类建筑 D、乙类建筑
3、计算框架结构梁截面惯性矩I时考虑楼板影响,对现浇楼盖,中框架取I=()。 A、2I0 B、1.5I0 C、1.2I0 D、I0 4、联肢剪力墙计算宜选用()分析方法。 A、材料力学分析法 B、连续化方法 C、壁式框架分析法 D、有限元法 5、框架-剪力墙结构侧移曲线为()。 A、弯曲型 B、剪切型 C、弯剪型 D、复合型 1、D 2、D 3、A 4、B 5、C 三、多项选择题:(将正确的答案的编号填入括弧中,完全选对才得分,否则不得分,每小题4分,共20分) 1、抗震设防结构布置原则() A、合理设置沉降缝 B、合理选择结构体系 C、足够的变形能力 D、增大自重 E、增加基础埋深 2、框架梁最不利内力组合有() A、端区-Mmax,+Mmax,Vmax B、端区Mmax及对应N,V C、中间+Mmax D、中间Mmax及对应N,V E、端区Nmax及对应M,V 3、双肢和多肢剪力墙内力和位移计算中假定() A、连梁反弯点在跨中 B、各墙肢刚度接近 C、考虑D 值修正 D、墙肢应考虑轴向变形影响 E、考虑反弯点修正 4、高层建筑剪力墙类型判别以()为主要特征。 A、等效刚度 B、墙肢宽度 C、整体性能 D、反弯点出现情况 E、墙肢高宽比 5、高层建筑基础埋置深度主要考虑()。 A、稳定性 B、施工便利性 C、抗震性 D、沉降量 E、增加自重
数值分析习题与答案
第一章绪论 习题一 1.设x>0,x*的相对误差为δ,求f(x)=ln x的误差限。解:求lnx的误差极限就是求f(x)=lnx的误差限,由公式(1. 2.4)有 已知x*的相对误差满足,而 ,故 即 2.下列各数都是经过四舍五入得到的近似值,试指出它们有几位有效数字,并给出其误差限与相对误差限。 解:直接根据定义和式(1.2.2)(1.2.3)则得 有5位有效数字,其误差限,相对误差限 有2位有效数字, 有5位有效数字, 3.下列公式如何才比较准确? (1) (2)
解:要使计算较准确,主要是避免两相近数相减,故应变换所给公式。 (1) (2) 4.近似数x*=0.0310,是 3 位有数数字。 5.计算取,利用:式计算误差最小。 四个选项: 第二、三章插值与函数逼近 习题二、三 1. 给定的数值表 用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解:仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计(5.8)。线性插值时,用0.5及0.6两点,用Newton插值 误差限,因
,故 二次插值时,用0.5,0.6,0.7三点,作二次Newton插值 误差限 ,故 2. 在-4≤x≤4上给出的等距节点函数表,若用二次插值法求的近似值,要使误差不超过,函数表的步长h 应取多少? 解:用误差估计式(5.8), 令 因 得 3. 若,求和.
解:由均差与导数关系 于是 4. 若互异,求 的值,这里p≤n+1. 解:,由均差对称性 可知当有 而当P=n+1时 于是得 5. 求证. 解:解:只要按差分定义直接展开得 6. 已知的函数表
数值计算方法答案
数值计算方法习题一(2) 习题二(6) 习题三(15) 习题四(29) 习题五(37) 习题六(62) 习题七(70) 2009.9,9
习题一 1.设x >0相对误差为2%4x 的相对误差。 解:由自变量的误差对函数值引起误差的公式: (())(())'()()()() f x x f x f x x f x f x δδ?= ≈得 (1)()f x = 11 ()()*2%1% 22x x δδδ≈ ===; (2)4 ()f x x =时 44 4 ()()'()4()4*2%8%x x x x x x δδδ≈ === 2.设下面各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出他们各有几位有效数字。 (1)12.1x =;(2)12.10x =;(3)12.100x =。 解:由教材9P 关于1212.m n x a a a bb b =±型数的有效数字的结论,易得上面三个数的有效 数字位数分别为:3,4,5 3.用十进制四位浮点数计算 (1)31.97+2.456+0.1352; (2)31.97+(2.456+0.1352) 哪个较精确? 解:(1)31.97+2.456+0.1352 ≈2 1 ((0.3197100.245610)0.1352)fl fl ?+?+ =2 (0.3443100.1352)fl ?+ =0.3457210? (2)31.97+(2.456+0.1352) 2 1 (0.319710(0.245610))fl fl ≈?+? = 21 (0.3197100.259110)fl ?+? =0.34562 10? 易见31.97+2.456+0.1352=0.3456122 10?,故(2)的计算结果较精确。 4.计算正方形面积时,若要求面积的允许相对误差为1%,测量边长所允许的相对误差限为多少?
高层建筑结构设计苏原第5章习题
第五章 5.1 平面结构和楼板在自身平面内具有无限刚性这两个基本假定是什么意义,在 框架、剪力墙、框架-剪力墙结构的近似计算中为什么要用这两个假定? 答:(1)假定一,一片框架或一片剪力墙可以抵抗在本身平面内的侧向力,而在平面外的刚度很小,可以忽略。因而整个结构可以划分成若干个平面结构共同抵抗与平面结构平行的侧向荷载,垂直于该方向的结构不参加受力。 假定二,楼板在其自身平面内刚度无限大,楼板平面外刚度很小,可以忽略。 因而在侧向力作用下,楼板可作刚体平移或转动,各个平面抗侧力结构之间通过楼板互相联系并协同工作。 上述两个基本假定的意义在于:近似方法将结构分成独立的平面结构单元,内力分析解决两个问题,第一,水平荷载在各片抗侧力结构之间的分配。荷载分配与抗侧力单元的刚度有关,要计算抗侧力单元的刚度,然后按刚度分配水平力,刚度愈大,分配的荷载也愈多。第二,计算每片平面结构在所分到的水平荷载作用下的内力和位移。 (2)在框架、剪力墙、框架-剪力墙结构的近似计算中要用这两个假定,这三大 结构体系的抗侧力构件均为平面构件,可以简化为平面结构,同时是为了简化计算,在不考虑扭转效应下,对计算的精度不会产生大的影响。 5.2分别画出一片三跨4层框架在垂直荷载(各层各跨满布均布荷载)和水平荷 载作用下的弯距图形、剪力图形和轴力图形。 5.3 刚度系数D和d的物理意义是什么?有什么区别?为什么?应用的条件是什么?应用时有哪些不同? 答:(1)D的物理意义:当柱端有转角时使柱端产生单位水平位移所需施加
的水平推力。d的物理意义:当柱端固定时使柱端产生单位水平位移所需施加的水平推力。 (2)抗侧刚度D值小于d值,即梁刚度较小时,柱的抗侧刚度减小了。因为 当梁的刚度较小时,对柱的约束作用减小,从而使柱的抗侧刚度减小。 (3)当梁比柱的抗弯刚度大很多时,刚度修正系数α值接近1,可近似认为α =1,此时第i层柱的侧移刚度为d值,在剪力分配公式中可用d值代替D i i 时可采用反 值,即反弯点法。工程中用梁柱线刚度比判断,当/35 b c 弯点法,反之,则采用D值法。 5.4 影响水平荷载下柱反弯点位置的主要因素是什么?框架顶层、底层和中部各 层反弯点位置有什么变化?反弯点高度比大于1的物理意义是什么? 答:(1)影响水平荷载下柱反弯点位置的主要因素:结构的总层数及该层所在位置;梁柱线刚度比;荷载形式;上层梁与下层梁刚度比;上下层层高比。 (2)在框架顶层反弯点位置在顶层柱中点以上;底层反弯点位置在2h/3高度 处(h是底层柱的高度);中部各层反弯点位置在各柱中点。 (3)反弯点高度比大于1的物理意义是表示柱下端的约束弯矩远大于柱上端 的约束弯矩,使得反弯点超过了柱的上端,使该柱中没有反弯点。 5.5梁柱杆件的弯曲变形和柱轴向变形对框架侧移有什么影响?框架为什么具有剪切型侧移曲线? 答:(1)框架总位移由杆件弯曲变形产生的侧移和柱轴向变形产生的侧移两部分叠加而成。由杆件弯曲变形引起的“剪切型侧移”,可由D值计算,为框架侧移的主要部分;由柱轴向变形产生的“弯曲型侧移”,可由连续化方法作近似估算。后者产生的侧移变形很小,多层框架可以忽略,当结构高度增大
数值计算方法》试题集及答案
《计算方法》期中复习试题 一、填空题: 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ,则用辛普生(辛卜生)公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案:2.367,0.25 2、1)3(,2)2(,1)1(==-=f f f ,则过这三点的二次插值多项式中2 x 的系数为 ,拉 格朗日插值多项式为 。 答案:-1, )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字; 4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( ); 答案 )(1)(1n n n n n x f x f x x x '--- =+ 5、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 ); 6、计算方法主要研究( 截断 )误差和( 舍入 )误差; 7、用二分法求非线性方程 f (x )=0在区间(a ,b )内的根时,二分n 次后的误差限为 ( 1 2+-n a b ); 8、已知f (1)=2,f (2)=3,f (4)=5.9,则二次Newton 插值多项式中x 2系数为( 0.15 ); 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f ),代数精度 为( 5 ); 12、 为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少,应将该表达 式改写为 11 ,))64(3(10-= -++=x t t t t y ,为了减少舍入误差,应将表达式1999 2001-
数值分析作业答案
数值分析作业答案 插值法 1、当x=1,-1,2时,f(x)=0,-3,4,求f(x)的二次插值多项式。 (1)用单项式基底。 (2)用Lagrange插值基底。 (3)用Newton基底。 证明三种方法得到的多项式是相同的。 解:(1)用单项式基底 设多项式为: , 所以: 所以f(x)的二次插值多项式为: (2)用Lagrange插值基底 Lagrange插值多项式为: 所以f(x)的二次插值多项式为: (3) 用Newton基底: 均差表如下: xk f(xk) 一阶均差二阶均差 1 0 -1 -3 3/2 2 4 7/ 3 5/6 Newton插值多项式为: 所以f(x)的二次插值多项式为: 由以上计算可知,三种方法得到的多项式是相同的。 6、在上给出的等距节点函数表,若用二次插值求ex的近似值,要使截断误差不超过10-6,问使用函数表的步长h应取多少? 解:以xi-1,xi,xi+1为插值节点多项式的截断误差,则有 式中 令得 插值点个数
是奇数,故实际可采用的函数值表步长 8、,求及。 解:由均差的性质可知,均差与导数有如下关系: 所以有: 15、证明两点三次Hermite插值余项是 并由此求出分段三次Hermite插值的误差限。 证明:利用[xk,xk+1]上两点三次Hermite插值条件 知有二重零点xk和k+1。设 确定函数k(x): 当或xk+1时k(x)取任何有限值均可; 当时,,构造关于变量t的函数 显然有 在[xk,x][x,xk+1]上对g(x)使用Rolle定理,存在及使得 在,,上对使用Rolle定理,存在,和使得 再依次对和使用Rolle定理,知至少存在使得 而,将代入,得到 推导过程表明依赖于及x 综合以上过程有: 确定误差限: 记为f(x)在[a,b]上基于等距节点的分段三次Hermite插值函数。在区间[xk,xk+1]上有 而最值 进而得误差估计: 16、求一个次数不高于4次的多项式,使它满足,,。
计算方法模拟试题及答案
计算方法模拟试题 一、 单项选择题(每小题3分,共15分) 1.近似值210450.0?的误差限为( )。 A . 0.5 B. 0.05 C . 0.005 D. 0.0005. 2. 求积公式)2(3 1 )1(34)0(31)(2 0f f f dx x f ++≈ ?的代数精确度为( )。 A. 1 B. 2 C. 3 D. 4 3. 若实方阵A 满足( )时,则存在唯一单位下三角阵L 和上三角阵R ,使LR A =。 A. 0det ≠A B. 某个0 det ≠k A C. )1,1(0det -=≠n k A k D. ),,1(0det n k A k =≠ 4.已知?? ?? ? ?????=531221112A ,则=∞A ( )。 A. 4 B. 5 C. 6 D 9 5.当实方阵A 满足)2(,221>>-=i i λλλλ,则乘幂法计算公式1e =( )。 A. 1+k x B. k k x x 11λ++ C. k x D. k k x x 11λ-+ 二、填空题(每小题3分,共15分) 1. 14159.3=π,具有4位有效数字的近似值为 。 2. 已知近似值21,x x ,则=-?)(21x x 。 3.已知1)(2-=x x f ,则差商=]3,2,1[f 。 4.雅可比法是求实对称阵 的一种变换方法。