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