第六章 数值分析模型§62非线性方程求根-§66估计水塔的水流量.doc

数 学 建 模 电 子 教 案第 14 次课课 题第六章数值分析模型 6.2非线性方程求根-6.6估计水塔的水流量教学内容1、求非线性方程根的方法：不动点迭代法，牛顿法，弦截法与抛物线法，2、数值分析模型：孩子成长和学生考试成绩问题，估计水塔的水流量教学目标1、掌握不动点迭代法，牛顿法，弦截法与抛物线法2、了解数值分析模型的建立、求解。教学重点求非线性方程根的方法教学难点数值分析模型的建立、求解。双语教学内容、安排 Numerical analysis model 数值分析模型 教学手段、措施以板演为主，多媒体与课堂讨论为辅作业、后记讨论题：P163: T3、T4教学过程及教学设计备注6.2非线性方程求根一、问题的提出本节主要讨论单变量非线性方程 （6-24） 求根问题，这里。在科学与工程计算中有大量方程求根问题，其中一类特殊的问题是多项式方程 （6-25）其中系数为实数。1、函数方程的解通常称为方程的根或函数的零点。特别地，如果函数可因式分解为，且，则称是函数的重零点或方程的重根。2、对于充分可微的函数，是函数的重零点的充分必要条件是二、二分法1、原理：设函数在上连续，且。不妨设，根据连续函数存在定理知道，在上一定有实根。2、计算步骤（1）计算区间中点 对于预先给定的小量，若，那么就是所求的根 若，此时根据的符号确定新的分割区间。继续下去，得到，显然有，故方程的根的近似值为三、不动点迭代法1、不动点迭代法将方程（6-24）改写成等价的形式 （6-26）若要求满足，则；反之亦然，称为函数的一个不动点。求的零点就等价于求的不动点，选择一个初始近似值，将它代入（6-26）右端，即可求得。 可以如此反复迭代计算， （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），在使用时如果且对任意有 （6-30）则由中值定理可知对有。它表明定理中的条件（2）可用（6-30）代替。3、局部收敛性与收敛阶上面给出了迭代序列在区间上的收敛性，通常称为全局收敛性。有时不易检验定理的条件，实际应用时通常只在不动点的邻近考察其收敛性，即局部收敛性。定义1 设有不动点，如果存在的某个邻域，对任意，迭代（6-26）产生的序列收敛且收敛到，则称迭代法（6-26）局部收敛。定理3 设有不动点，在的某个邻域连续，且，则迭代法（6-26）局部收敛。定义2 设迭代过程 收敛于方程的根，如果迭代误差满足下列渐进关系式。则称该迭代过程是阶收敛的，特别地，时称线性收敛，时称超线性收敛，时称平方收敛。定理4 对于迭代过程，如果在所求根的邻近连续，并且 ，则该迭代过程在点 邻近是阶收敛的。上述定理告诉我们，迭代过程的收敛速度依赖于迭代函数的选取，如果当 时，则该迭代过程只可能是线性收敛。6.3 牛顿法及其收敛性1、牛顿法设已知方程有近似根（假定 ），将函数在点展开，有于是方程可近似地表示为 （6-32）则可得近似计算公式为 （6-33）这就是牛顿（Newton）法。关于牛顿法（6-33）的收敛性，可直接由定理4得到，牛顿法在根的邻近是平方收敛的。又因，可得。 （6-34） 2、简化牛顿法与牛顿下山法(1) 简化牛顿法，也称平行弦法，其迭代公式为， （6-35）迭代函数。若，即取。在根附近成立，则迭代法（6-35）局部收敛。在（6-35）中取 ，则称为简化牛顿法，这类方法计算量省，但只有线性收敛，其几何意义是用平行弦与轴交点作为的近似。（2）牛顿下山法。为了防止迭代发散，我们对迭代过程再附加一项要求，即具有单调性： （6-36） 满足这项要求的算法称为下山法。我们将牛顿法和下山法结合起来使用，即在下山法保证函数值稳定下降的前提下，用牛顿法加快收敛速度。为此，我们将牛顿法的计算结果 与前一步的近似值适当加权平均作为新的改进值 （6-37） 其中 称为下山因子，（6-37）即为 （6-38）称为牛顿下山法。选择下山因子时从开始，逐次将减半进行试算，直到能使下降条件（6-36）成立为止。6.4 弦截法与抛物线法用牛顿法求方程根，每步除计算外还要算，而计算 往往较困难，为此下面介绍两种常用方法。1、弦截法设是的近似根，我们利用构造一次插值多项式，并用的根作为的新的近似根。由于 （6-39） 因此有 （6-40）2、抛物线法（密勒（Mller）法）插值多项式有两个零点 （6-41）式中 为了从（6-41）中定出一个值，我们需要讨论根式前正负号的取舍问题。在三个近似根中，自然假定更接近所求的根，这时，为了保证精度，我们选式（6-41）中较接近的值作为新的近似根。为此，只要取根式前的符号与的符号相同。6.5 孩子成长和学生考试成绩问题（1）孩子成长问题一个男孩在11岁长到21岁过程中，身高的变化如表5-1所示，试找一个最佳的函数曲线来表示这个男孩的成长过程。见表5-1 ti（从11岁起年龄）0 0.8 1.4 2.0 2.4 3.2 4.0 增长高度hi（cm） 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 增长高度hi（cm） 26.25 28.88 30.60 32.25 33 35 （2）学生考试成绩问题：对8个学生测量其智商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 以上两个问题，均可以抽象为:通过实验等方法观测到反映某个函数的数据。要求利用这些数据构造出的近似表达式。上一节介绍的插值法就是寻求近似函数的方法之一。但由于实验观测数据不可避免地带有误差，甚至是较大的误差，所以使用插值法是不合适的，它会保留数据的误差。因此，不必要求近似函数满足，而只要求偏差按某种标准最小，以反映所给数据的总体趋势，消除局部波动的影响，这就是曲线拟合问题。6.6 估计水塔的水流量某居民区的民用自来水是由一个圆柱形的水塔提供，一般可以通过测量其水位来估计水的流量。但问题是，当 水塔水位下降到设定的最低水位时水泵自动启动向水塔供水，到设定的最高水位时停止供水，在水泵自动加水期间无法测量水塔的水位和水泵的供水量。通常水泵每天供水两次，每次约两小时。水塔是一个高12.2米，直径17.4米的正圆柱，当水塔的水位降至最低水位，约8.2米时，水泵自动启动供水当水塔的水位升高到一个最高水位，约10.8米时，水泵停止工作。表5-4是某一天的水位测量记录数据，测量了28个时刻，但是由于其中有4个时刻遇到水泵正在向水塔供水，而无水位记录，表5-4中用符号表示。试估计任何时刻（包括水泵正供水时）的用水率，及一天的总用水量。见表6-1模型假设及符号说明（1）模型假设 假设水塔中流出的水流量只受社区的日常生活需要的影响，水的消耗每天大致差不多。 由Torricelli定律知，从水塔流出的最大流速正比于水位高度的平方根，题目中给出水塔的最高和最低水位分别为10.82米和8.22米，所以对于这两种高度，最大水流速度的比约为 ，这表明我们可以假设水塔中水位对水流速度影响忽略不计。 水泵工作时单位时间的供水量大致为常数，这个常数大于单位时间内从水塔中流出的水流的最大流速，这是因为居民区内一直需要用水，不允许水塔中的水用光。 水塔中水流量是时间的连续光滑函数，与水泵工作与否无关。这是因为虽然就个别用户而言可能用水量有较大的变化，但由于个人的用水量与整个居民区用水量相比是非常小的，从统计意义上来讲，不太可能同时整个社区的用水量增长或减少。 水泵工作与否完全取决于水塔内水位的高度，且每次加水的工作时间为2小时，根据表6-1中的数据可知，水泵第一次供水时间段为 8.967,10.954，第二次供水时间段为20.839, 23.880。 符号说明：测量的时刻；：水位的高度；：水塔中水的体积；：水塔中水流速度，即流量；：第i时段用水量，即没有供水时间段用水量； ：第1和第2供水时间段用水量之和；：一天中总用水量。（2）问题分析 问题要求任意时刻的用水率，即求单位时间流出的水的体积，一般称为水流速度或流量。由于水塔是一个圆柱体，体积可以很容易地通过水位高度h计算出来，这样在水泵不工作的时间段，水流速度就可以从体积对时间的导数计算出来，由于没有水的体积关于时间的函数表达式，而只能利用问题中给定的原始数据表6-1，通过体积公式计算出离散的在测量时刻的体积V，因此可以考虑用差商代替微商，也就是用离散代替连续的思想。为提高计算精度，采用二阶差商，即。由于所有数据被水泵两次供水分割成三组数据对每组数据的中间数据采用中心差商，前后两个数据不能采用中心差商，改用向前差商、向后差商或用中点公式进行差商。 中心差商公式： 向前差商公式： 向后差商公式： 中点公式：以上分析了水泵不工作的时段，用水率的计算。对于水泵供水时段的用水率，计算难度较大，我们只好用供水时间段前后的用水率进行插值或拟合而得到。有了任何时刻的用水率，可以采用数值积分计算一天的总用水量。（3）模型建立及求解通过以上对问题的分析，现在的问题已转化为根据某一天已测量的时刻水塔中水的流速，产生在整个区间（24小时）上的函数或函数值，一般来说插值和拟合是两种最常用的方法。分两步进行： 计算水流速度并画出散点图 计算水塔中水的体积（4）模型的建立通过对不同插值方法的比较，结合假设，考虑到流速应该是时间的连续光滑函数，所以采用三次样条插值模型计算用水率函数。计算用水率首先用三次样条插值计算用水率函数，MATLAB程序执行后，可得样条插值下的流速图计算一天的总用水量用三次样条插值模型得到的函数在时间区间0，24上做数值积分可得一天的总用水量。程序执行后，可得一天的总用水量1266.1立方米。见图6-2（5）模型稳定性分析用不同时刻作为起始点，使用三次样条插值模型得到的用水率函数，在长度为24小时区间上进行数值积分，计算程序为：执行文件后，所得结果列入表6-2。（单位：立方米）从表6-2易见，选择不同的起始点，所得结果比较稳定。（6）模型检验水泵未工作时的检验首先计算水泵未工作的三段时间内的实际用水量第一时间段，即时的实际用水量为：第二时间段，即 时的实际用水量为：237.787(10.82-8.22)=618.25第三时间段，即时的实际用水量为：237.787(10.591-10.18)=97.73其次用三次样条插值模型得到的函数分别在水泵未工作的三段时间内作数值积分，执行程序后，可得水泵未工作的三段时间内，用模型计算出的用水量（见表6-3）。分三段的实际用水量与模型用水量的比较列成表6-3。并且从表6-3可见，用三次样条插值模型计算的用水量与实际