《第六章数值研究模型§非线性方程求根§估计水塔的水流量 .docx》由会员分享，可在线阅读，更多相关《第六章数值研究模型§非线性方程求根§估计水塔的水流量 .docx（8页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精品名师归纳总结课 题第六章数值分析模型 6.2 非线性方程求根 - 6.6 估量水塔的水流量教案内容1、求非线性方程根的方法：不动点迭代法，牛顿法，弦截法与抛物线法， 2、数值分析模型：孩子成长和同学考试成果问题，估量水塔的水流量教案目标1、把握不动点迭代法，牛顿法，弦截法与抛物线法2、明白数值分析模型的建立、求解。教案重点求非线性方程根的方法教案难点数值分析模型的建立、求解。可编辑资料 - - - 欢迎下载精品名师归纳总结双 语 教 案内 容 、 支配教 案 手段、措施作 业 、 后记Numerical analysis model数值分析模型以板演为主，多媒体与课堂争论为辅争论题： P1
2、63: T3 、T4教 案 过 程 及 教 案 设 计备注可编辑资料 - - - 欢迎下载精品名师归纳总结一、问题的提出 6.2 非线性方程求根对教案 内可编辑资料 - - - 欢迎下载精品名师归纳总结本节主要争论单变量非线性方程6-24）求根问题，这里。在科学与工程运算中有大量方程求根问题，其中一类特殊的问题是多项式方程25）其中系数为实数。1、函数方程的解通常称为方程的根或函数的零点。特殊的，假如函数可因式分解容 及欲 达目6-的 、讲 授方 法加 以说可编辑资料 - - - 欢迎下载精品名师归纳总结为，且，就称是函数的重零点或方程的重根。明）2、对于充分可微的函数，是函数的重零点的充分必
3、要条件是二、二分法1 、原理：设函数在上连续，且。不妨设，依据连续函数存在定理知道，在上肯定有实根。2、运算步骤可编辑资料 - - - 欢迎下载精品名师归纳总结1）运算区间中点 对于预先给定的小量，如，那么就是所求的根 如，此时依据的符号确定新的分割区间。连续下去，得到，明显有，故方程的根的近似值为可编辑资料 - - - 欢迎下载精品名师归纳总结三、不动点迭代法1 、不动点迭代法将方程 6-24 ）改写成等价的形式26 ）6-可编辑资料 - - - 欢迎下载精品名师归纳总结如要求满意，就。反之亦然，称为函数的一个不动点。求的零点就等价于求的不动点，挑选一个初始近似值，将它代入 6-26 ）右端
4、，即可求得。可以如此反复迭代运算，6-27）称为迭代函数。假如对任何，由 6-27 ）得到的序列有极限。就称迭代方程 6-27 ）收敛，且为的不动点，故称6-27 ）为不动点迭代法。2、不动点的存在性与迭代法的收敛性1 ）不动点的存在唯独性。定理 1 设满意以下两个条件：1 ）对任意的有。2 ）存在正常，使对任意都有6-28 ） 就在上存在唯独的不动点。在的不动点存在唯独的情形下，可得到迭代法6-27 ）收敛的一个充分条件。定理 2设满意定理 1 中的两个条件，就对任意，由6-27 ）得到的迭代序列收敛到的不动点，并有误差估量6-29）对定理 1和定理 2 中的条件 2 ），在使用时假如且对任
5、意有6-30 ）就由中值定理可知对有。它说明定理中的条件 2 ）可用 6-30 ）代替。3 、局部收敛性与收敛阶上面给出了迭代序列在区间上的收敛性，通常称为全局收敛性。有时不易检验定理的条件，实际应用时通常只在不动点的邻近考察其收敛性，即局部收敛性。可编辑资料 - - - 欢迎下载精品名师归纳总结定义 1设有不动点，假如存在的某个邻域，对任意，迭代6-26 ）产生的序列收敛且收敛到，就称迭代法 6-26 ）局部收敛。定理 3 设有不动点，在的某个邻域连续，且，就迭代法 6-26 ）局部收敛。定义 2设迭代过程收敛于方程的根，假如迭代误差满意以下渐进关系式。就称该迭代过程是阶收敛的，特殊的，时称
6、线性收敛，时称超线性收敛， 时称平方收敛。定理 4 对于迭代过程，假如在所求根的邻近连续，并且，就该迭代过程在点邻近是阶收敛的。上 述 定 理 告 诉 我 们 ， 迭 代 过 程 的 收 敛 速 度 依 赖 于 迭 代 函 数的 选 取 ， 如 果 当 时，就该迭代过程只可能是线性收敛。 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）可编辑资料 - - - 欢迎下载精品名师归纳总结我们将牛顿法和下山法结合起来使用，即在下山法保证函数值稳固下降的前提下，用牛顿法加快收敛速度。为此，我们将牛顿法的运算结果可编辑资料 - - - 欢迎下载
8、精品名师归纳总结与前一步的近似值适当加权平均作为新的改进值6-37）其中称为下山因子，6-37 ）即为6-38）称为牛顿下山法。挑选下山因子时从下降条件 6-36 ）成立为止。开头，逐次将 6.4 弦截法与抛物线法减半进行试算，直到能使用牛顿法求方程根，每步除运算外仍要算，而运算往往较困难，为此下面介绍两种常用方法。1 、弦截法设是的近似根，我们利用构造一次插值多项式， 并用的根作为的新的近似根。由于6-39）因此有6-40）2 、抛物线法 密勒 M ller ）法） 插值多项式有两个零点6-41）式中为了从 6-41 ）中定出一个值，我们需要争论根式前正负号的取舍问题。在三个近似根中，自然假
9、定更接近所求的根，这时，为了保证精度，我们 选式6-41 ）中较接近的值作为新的近似根。为此，只要取根式前的符号与的符号相同。 6.5 孩子成长和同学考试成果问题1）孩子成长问题一个男孩在 11 岁长到 21 岁过程中，身高的变化如表5-1 所示，试找一个正确的函数曲线来表示这个男孩的成长过程。见表5-1ti 从 11 岁起年龄）00.81.42.02.43.24.0增长高度 hicm ）00.742.255.258.2515.0021.38ti 从 11 岁起年龄）4.8 5.46.07.08.010.0增长高度 hicm ）2 ）同学考试成果问题：26.2528.8830.6032.253
10、335可编辑资料 - - - 欢迎下载精品名师归纳总结对 8个同学测量其智商Iq 和课后复习某门课时间t 及该门课考试成果 g，得表 5-2 ，试争论该门课考试成果与智商和课后复习时间之间的关系。见表 5-2 。智商 105110120116122130114102Iq ）复习时间t ）10126131682015考试成果g）7579688591799876以上两个问题，均可以抽象为: 通过试验等方法观测到反映某个函数的数据。要求利用这些数据构造出的近似表达式。上一节介绍的插值法就是寻求近似函数的方法之一。但由于试验观测数据不行防止的带有误差，甚至是较大的误差，所以使用插值法是不合适的，它会保
11、留数据的误差。因此，不必要求近似函数满意，而只要求偏差按某种标准最小，以反映所给数据的总体趋势，排除局部波动的影响，这就是曲线拟合问题。 6.6 估量水塔的水流量某居民区的民用自来水是由一个圆柱形的水塔供应，一般可以通过测量其水位来估量水的流量。但问题是，当水塔水位下降到设定的最低水位时水泵自动启动向水塔供水，到设定的最高水位时停止供水，在水泵自动加水期间无法测量水塔的水位和水泵的供水量。通常水泵每天供水两次，每次约两小时。水塔是一个高12.2M ，直径 17.4M 的正圆柱，当水塔的水位降至最低水位，约 8.2M 时，水泵自动启动供水当水塔的水位上升到一个最高水位，约10.8M 时，水泵停止
12、工作。表5-4是某一天的水位测量记录数据，测量了28 个时刻，但是由于其中有 4个时刻遇到水泵正在向水塔供水，而无水位记录，表5-4 中用符号表示。试估量任何时刻 包括水泵正供水时）的用水率，及一天的总用水量。见表6-1时刻t00.9211.8432.9493.8714.9785.900水位9.6779.4799.3089.1258.9828.8148.686时刻t7.0067.9288.9679.981110.92510.95412.032水位8.5258.3888.22010.82010.500时刻t12.95413.87514.98215.90316.82617.93119.037水位1
13、0.2109.9369.6539.4099.1808.9218.662时刻t19.95920.83922.01522.95823.88024.98625.908水位8.4338.22010.59110.35410.180模型假设及符号说明可编辑资料 - - - 欢迎下载精品名师归纳总结1 ）模型假设假设水塔中流出的水流量只受社区的日常生活需要的影响，水的消耗每天大致差不多。由Torricelli 定律知，从水塔流出的最大流速正比于水位高度的平方根，题目中给出水塔的最高和最低水位分别为10.82M 和8.22M ，所以对于这两种高度，最大水流速度的比约为，这说明我们可以假设水塔中水位对水流速度影
14、响忽视不计。水泵工作时单位时间的供水量大致为常数，这个常数大于单位时间内从水塔中流出的水流的最大流速，这是由于居民区内始终需要用水，不答应水塔中的水用光。水塔中水流量是时间的连续光滑函数，与水泵工作与否无关。这是由于虽然就个别用户而言可能用水量有较大的变化，但由于个人的用水量与整个居民区用水量相比是特别小的，从统计意义上来讲，不太可能同时整个社区的用水量增长或削减。水泵工作与否完全取决于水塔内水位的高度，且每次加水的工作时间为 2 小时，依据表6-1 中的数据可知，水泵第一次供水时间段为 8.967,10.954 ，其次次供水时间段为20.839, 23.880 。. 符号说明 ：测量的时刻。
15、：水位的高度。：水塔中水的体积。：水塔中水流速度，即流量。：第i时段用水量，即没有供水时间段用水量。：第 1和第 2供水时间段用水量之和。：一天中总用水量。 2 ）问题分析问题要求任意时刻的用水率，即求单位时间流出的水的体积，一般称为水流速度或流量。由于水塔是一个圆柱体，体积可以很简单的通过水位高度h运算出来，这样在水泵不工作的时间段，水流速度就可以从体积对时间的导数运算出来，由于没有水的体积关于时间的函数表达式，而只能利用问题中给定的原始数据表6-1 ，通过体积公式运算出离散的在测量时刻的体积V，因此可以考虑用差商代替微商，也就是用离散代替连续的思想。为提高运算精度，采纳二阶差商，即。由于全
16、部数据被水泵两次供水分割成三组数据对每组数据的中间数据采纳中心差商，前后两个数据不能采纳中心差商，改用向前差商、向后差商或用中点公式进行差商。中心差商公式： 向前差商公式：向后差商公式：中点公式：以上分析了水泵不工作的时段，用水率的运算。对于水泵供水时段的用水率，运算难度较 大，我们只好用供水时间段前后的用水率进行插值或拟合而得到。有了任何时刻的用水率， 可以采纳数值积分运算一天的总用水量。3）模型建立及求解通过以上对问题的分析，现在的问题已转化为依据某一天已测量的时刻水塔中水的流 速，产生在整个区间 24小时）上的函数或函数值，一般来说插值和拟合是两种最常用的方法。分两步进行：运算水流速度并
17、画出散点图运算水塔中水的体积4 ）模型的建立通过对不同插值方法的比较，结合假设，考虑到流速应当是时间的连续光滑函数，所 以采纳三次样条插值模型运算用水率函数。运算用水率第一用三次样条插值运算用可编辑资料 - - - 欢迎下载精品名师归纳总结水率函数， MATLAB 程序执行后，可得样条插值下的流速图运算一天的总用水量用三次样条插值模型得到的函数在时间区间 0 ，24 上做数值积分可得一天的总用水量。程序执行后，可得一天的总用水量1266.1 立方 M。见图 6-25 ）模型稳固性分析用不同时刻作为起始点，使用三次样条插值模型得到的用水率函数，在长度为 24小时区间上进行数值积分，运算程序为：执
18、行文件后，所得结果列入表6-2 。 单位：立方 M）从表 6-2 易见，挑选不同的起始点，所得结果比较稳固。=618.25第三时间段，即时的实际用水量为： 237.787 10.591-10.18=97.73其次用三次样条插值模型得到的函数分别在水泵未工作的三段时间内作数值积分，执行程序后，可得水泵未工作的三段时间内，用模型运算出的用水量见表 6-3 ）。分三段的实际用水量与模型用水量的比较列成表6-3 。并且从表 6-3 可见，用三次样条插值模型运算的用水量与实际用水量吻合的很好。其次水泵工作时的检验水泵充水量 =充水后的水量 +充水期间的流出量- 充水前的水量。第一利用模型运算两次充水期间的流出量，执行后可得：第一次充水期间的用水量为117.95 ，其次次充水期间的用水量为179.06 。 进 一 步 可 计 算 第 一 次 水 泵 充 水 量 =237.787 10.82+117.95-237.7878.22=736.20第 二 次 水 泵 充 水 量 =237.787 10.591+179.06-237.787 8.22=742.85易见，水泵两次注水量相差不大。可编辑资料 - - - 欢迎下载