《2022年第六章数值研究模型§非线性方程求根§估计水塔的水流量》由会员分享,可在线阅读,更多相关《2022年第六章数值研究模型§非线性方程求根§估计水塔的水流量(7页珍藏版)》请在金锄头文库上搜索。
1、1 / 7 课 题第六章数值分析模型 6.2 非线性方程求根 - 6.6 估计水塔的水流量教案内容1、求非线性方程根的方法:不动点迭代法,牛顿法,弦截法与抛物线法,2、数值分析模型:孩子成长和学生考试成绩问题,估计水塔的水流量教案目标1、掌握不动点迭代法,牛顿法,弦截法与抛物线法2、了解数值分析模型的建立、求解。教案重点求非线性方程根的方法教案难点数值分析模型的建立、求解。双 语 教 案内 容 、 安排Numerical analysis model 数值分析模型教案手段、措施以板演为主,多媒体与课堂讨论为辅作 业 、 后记讨论题: P163: T3 、T4教 案 过 程 及 教 案 设 计备
2、注6.2 非线性方程求根一、问题的提出本节主要讨论单变量非线性方程 6-24)求根问题,这里。在科学与工程计算中有大量方程求根问题,其中一类特殊的问题是多项式方程 6-25)其中系数为实数。1、函数方程的解通常称为方程的根或函数的零点。特别地,如果函数可因式分解为,且,则称是函数的重零点或方程的重根。2、对于充分可微的函数,是函数的重零点的充分必要条件是二、二分法1、原理:设函数在上连续,且。不妨设,根据连续函数存在定理知道,在上一定有实根。2、计算步骤对教案 内容 及欲 达目的 、讲 授方 法加 以说明)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - -
3、-第 1 页,共 7 页2 / 7 1)计算区间中点 对于预先给定的小量,若,那么就是所求的根 若,此时根据的符号确定新的分割区间。继续下去,得到,显然有,故方程的根的近似值为三、不动点迭代法1、不动点迭代法将方程 6-24 )改写成等价的形式 6-26)若要求满足,则;反之亦然,称为函数的一个不动点。求的零点就等价于求的不动点,选择一个初始近似值,将它代入 6-26)右端,即可求得。可以如此反复迭代计算, 6-27)称为迭代函数。如果对任何,由 6-27 )得到的序列有极限。则称迭代方程6-27 )收敛,且为的不动点,故称6-27 )为不动点迭代法。2、不动点的存在性与迭代法的收敛性1)不动
4、点的存在唯一性。定理 1设满足以下两个条件:1)对任意的有。2)存在正常,使对任意都有 6-28 )则在上存在唯一的不动点。在的不动点存在唯一的情况下,可得到迭代法6-27 )收敛的一个充分条件。定理 2 设满足定理 1中的两个条件,则对任意,由 6-27 )得到的迭代序列收敛到的不动点,并有误差估计 6-29)对定理 1和定理 2中的条件 2),在使用时如果且对任意有6-30 )则由中值定理可知对有。它表明定理中的条件2)可用 6-30 )代替。3、局部收敛性与收敛阶上面给出了迭代序列在区间上的收敛性,通常称为全局收敛性。有时不易检验定理的条件,实际应用时通常只在不动点的邻近考察其收敛性,即
5、局部收敛性。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共 7 页3 / 7 定义 1 设有不动点,如果存在的某个邻域,对任意,迭代6-26 )产生的序列收敛且收敛到,则称迭代法 6-26 )局部收敛。定理 3 设有不动点,在的某个邻域连续,且,则迭代法 6-26 )局部收敛。定义 2 设迭代过程收敛于方程的根,如果迭代误差满足下列渐进关系式。则称该迭代过程是阶收敛的,特别地,时称线性收敛,时称超线性收敛,时称平方收敛。定理 4 对于迭代过程,如果在所求根的邻近连续,并且,则该迭代过程在点邻近是阶收敛的。上 述 定 理 告 诉 我 们
6、, 迭 代 过 程 的 收 敛 速 度 依 赖 于 迭 代 函 数的 选 取 , 如 果 当 时,则该迭代过程只可能是线性收敛。6.3 牛顿法及其收敛性1、牛顿法设已知方程有近似根假定),将函数在点展开,有于是方程可近似地表示为 6-32)则可得近似计算公式为 6-33)这就是牛顿 Newton )法。关于牛顿法 6-33 )的收敛性,可直接由定理4得到,牛顿法在根的邻近是平方收敛的。又因,可得。 6-34)2、简化牛顿法与牛顿下山法(1) 简化牛顿法,也称平行弦法,其迭代公式为, 6-35)迭代函数。若,即取。在根附近成立,则迭代法6-35 )局部收敛。在6-35 )中取,则称为简化牛顿法,
7、这类方法计算量省,但只有线性收敛,其几何意义是用平行弦与轴交点作为的近似。2)牛顿下山法。为了防止迭代发散,我们对迭代过程再附加一项要求,即具有单调性: 6-36)满足这项要求的算法称为下山法。我们将牛顿法和下山法结合起来使用,即在下山法保证函数值稳定下降的前提下,用牛顿法加快收敛速度。为此,我们将牛顿法的计算结果精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 7 页4 / 7 与前一步的近似值适当加权平均作为新的改进值 6-37)其中称为下山因子,6-37 )即为 6-38)称为牛顿下山法。选择下山因子时从开始,逐次将减半进行试算,直
8、到能使下降条件 6-36 )成立为止。6.4 弦截法与抛物线法用牛顿法求方程根,每步除计算外还要算,而计算往往较困难,为此下面介绍两种常用方法。1、弦截法设是的近似根,我们利用构造一次插值多项式,并用的根作为的新的近似根。由于 6-39)因此有 6-40)2、抛物线法 密勒 M ller )法)插值多项式有两个零点 6-41)式中为了从 6-41 )中定出一个值,我们需要讨论根式前正负号的取舍问题。在三个近似根中,自然假定更接近所求的根,这时,为了保证精度,我们选式 6-41 )中较接近的值作为新的近似根。为此,只要取根式前的符号与的符号相同。6.5 孩子成长和学生考试成绩问题1)孩子成长问题
9、一个男孩在11 岁长到 21 岁过程中,身高的变化如表5-1 所示,试找一个最佳的函数曲线来表示这个男孩的成长过程。见表5-1ti 从 11 岁起年龄)0 0.8 1.4 2.0 2.4 3.2 4.0 增长高度hicm)0 0.74 2.25 5.25 8.25 15.00 21.38 ti 从 11 岁起年龄)4.8 5.4 6.0 7.0 8.0 10.0 增长高度hicm)26.25 28.88 30.60 32.25 33 35 2)学生考试成绩问题:精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 7 页5 / 7 对 8个学
10、生测量其智商Iq 和课后复习某门课时间t 及该门课考试成绩g,得表 5-2,试研究该门课考试成绩与智商和课后复习时间之间的关系。见表 5-2 。智商Iq)105 110 120 116 122 130 114 102 复习时间t )10 12 6 13 16 8 20 15 考试成绩g)75 79 68 85 91 79 98 76 以上两个问题,均可以抽象为: 通过实验等方法观测到反映某个函数的数据。要求利用这些数据构造出的近似表达式。上一节介绍的插值法就是寻求近似函数的方法之一。但由于实验观测数据不可避免地带有误差,甚至是较大的误差,所以使用插值法是不合适的,它会保留数据的误差。因此,不必
11、要求近似函数满足,而只要求偏差按某种标准最小,以反映所给数据的总体趋势,消除局部波动的影响,这就是曲线拟合问题。6.6 估计水塔的水流量某居民区的民用自来水是由一个圆柱形的水塔提供,一般可以通过测量其水位来估计水的流量。但问题是,当水塔水位下降到设定的最低水位时水泵自动启动向水塔供水,到设定的最高水位时停止供水,在水泵自动加水期间无法测量水塔的水位和水泵的供水量。通常水泵每天供水两次,每次约两小时。水塔是一个高12.2M ,直径 17.4M 的正圆柱,当水塔的水位降至最低水位,约8.2M 时,水泵自动启动供水当水塔的水位升高到一个最高水位,约10.8M 时,水泵停止工作。表5-4是某一天的水位
12、测量记录数据,测量了28个时刻,但是由于其中有 4个时刻遇到水泵正在向水塔供水,而无水位记录,表5-4中用符号表示。试估计任何时刻 包括水泵正供水时)的用水率,及一天的总用水量。见表6-1 模型假设及符号说明时刻t0 0.921 1.843 2.949 3.871 4.978 5.900 水位9.677 9.479 9.308 9.125 8.982 8.814 8.686 时刻t7.006 7.928 8.967 9.9811 10.925 10.954 12.032 水位8.525 8.388 8.220 10.820 10.500 时刻t12.954 13.875 14.982 15.9
13、03 16.826 17.931 19.037 水位10.210 9.936 9.653 9.409 9.180 8.921 8.662 时刻t19.959 20.839 22.015 22.958 23.880 24.986 25.908 水位8.433 8.220 10.591 10.354 10.180 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 7 页6 / 7 1)模型假设假设水塔中流出的水流量只受社区的日常生活需要的影响,水的消耗每天大致差不多。由Torricelli 定律知,从水塔流出的最大流速正比于水位高度的平方根,
14、题目中给出水塔的最高和最低水位分别为10.82M 和8.22M ,所以对于这两种高度,最大水流速度的比约为,这表明我们可以假设水塔中水位对水流速度影响忽略不计。水泵工作时单位时间的供水量大致为常数,这个常数大于单位时间内从水塔中流出的水流的最大流速,这是因为居民区内一直需要用水,不允许水塔中的水用光。水塔中水流量是时间的连续光滑函数,与水泵工作与否无关。这是因为虽然就个别用户而言可能用水量有较大的变化,但由于个人的用水量与整个居民区用水量相比是非常小的,从统计意义上来讲,不太可能同时整个社区的用水量增长或减少。水泵工作与否完全取决于水塔内水位的高度,且每次加水的工作时间为2小时,根据表6-1中
15、的数据可知,水泵第一次供水时间段为8.967,10.954,第二次供水时间段为20.839, 23.880。?符号说明:测量的时刻;:水位的高度;:水塔中水的体积;:水塔中水流速度,即流量;:第 i时段用水量,即没有供水时间段用水量;:第 1和第 2供水时间段用水量之和;:一天中总用水量。2)问题分析问题要求任意时刻的用水率,即求单位时间流出的水的体积,一般称为水流速度或流量。由于水塔是一个圆柱体,体积可以很容易地通过水位高度h计算出来,这样在水泵不工作的时间段,水流速度就可以从体积对时间的导数计算出来,由于没有水的体积关于时间的函数表达式,而只能利用问题中给定的原始数据表6-1 ,通过体积公
16、式计算出离散的在测量时刻的体积V,因此可以考虑用差商代替微商,也就是用离散代替连续的思想。为提高计算精度,采用二阶差商,即。由于所有数据被水泵两次供水分割成三组数据对每组数据的中间数据采用中心差商,前后两个数据不能采用中心差商,改用向前差商、向后差商或用中点公式进行差商。中心差商公式:向前差商公式:向后差商公式:中点公式:以上分析了水泵不工作的时段,用水率的计算。对于水泵供水时段的用水率,计算难度较大,我们只好用供水时间段前后的用水率进行插值或拟合而得到。有了任何时刻的用水率,可以采用数值积分计算一天的总用水量。3)模型建立及求解通过以上对问题的分析,现在的问题已转化为根据某一天已测量的时刻水
17、塔中水的流速,产生在整个区间24小时)上的函数或函数值,一般来说插值和拟合是两种最常用的方法。分两步进行:计算水流速度并画出散点图计算水塔中水的体积4)模型的建立通过对不同插值方法的比较,结合假设,考虑到流速应该是时间的连续光滑函数,所以采用三次样条插值模型计算用水率函数。计算用水率首先用三次样条插值计算用精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页,共 7 页7 / 7 水率函数,MATLAB 程序执行后,可得样条插值下的流速图计算一天的总用水量用三次样条插值模型得到的函数在时间区间 0,24 上做数值积分可得一天的总用水量。程序执行后
18、,可得一天的总用水量1266.1 立方 M。见图 6-2 5)模型稳定性分析用不同时刻作为起始点,使用三次样条插值模型得到的用水率函数,在长度为 24小时区间上进行数值积分,计算程序为:执行文件后,所得结果列入表6-2。单位:立方 M)从表 6-2易见,选择不同的起始点,所得结果比较稳定。=618.25 第三时间段,即时的实际用水量为:237.787 (10.591-10.18=97.73其次用三次样条插值模型得到的函数分别在水泵未工作的三段时间内作数值积分,执行程序后,可得水泵未工作的三段时间内,用模型计算出的用水量见表 6-3 )。分三段的实际用水量与模型用水量的比较列成表6-3 。并且从
19、表 6-3可见,用三次样条插值模型计算的用水量与实际用水量吻合的很好。其次水泵工作时的检验水泵充水量=充水后的水量+充水期间的流出量- 充水前的水量。首先利用模型计算两次充水期间的流出量,执行后可得:第一次充水期间的用水量为117.95 ,第二次充水期间的用水量为179.06 。 进 一 步 可 计 算 第 一 次 水 泵 充 水 量 =237.787 10.82+117.95-237.7878.22=736.20 第二次水泵充水量=237.78710.591+179.06-237.7878.22=742.85 易见,水泵两次注水量相差不大。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 7 页