第六章数值研究模型§非线性方程求根§估计水塔的水流量

1、课题第六章数值分析模型 6.2非线性方程求根- 6.6估计水塔的水流量教案内容1、求非线性方程根的方法：不动点迭代法，牛顿法，弦截法与抛物线法,2、数值分析模型：孩子成长和学生考试成绩问题，估计水塔的水流量教案目标1、掌握不动点迭代法，牛顿法，弦截法与抛物线法2、了解数值分析模型的建立、求解。教案重点求非线性方程根的方法教案难点数值分析模型的建立、求解。双语教案 内容、安 排Numerical a nalysis model数值分析模型教案手段、措施以板演为主，多媒体与课堂讨论为辅作业、后 记讨论题：P163: T3、T4教案过程及教案设计备注 6.2非线性方程求根对教一、问题的提出本节主要讨

2、论单变量非线性方程巨06-24)案 容 欲内 及 达求根冋题，这 在科学与工程文里1。扛十算中有大量方程求根问题，其中一类特殊的问题是多项式方程=6-目 的、讲授25)方法其中系数I II为实数。加 说 明)以1、函数方程的解通常称为方程的根或函数的零点。特别地，如果函数可因式分解为1日1 乂 1则称也 是函数的旦 重零占或方程的 M重根。丿 7?八 7 r J -/Jy 八、-八八 1-2、对于充分可微的函数旦，回是函数的回重零点的充分必要条件是二、二分法1、原理：设函数在上J上连续，且1 一 1。不妨设据连续函数存在定理知道，在E上一定有实根。1 ，根2、计算步骤1）计算区间中点 a对于预

3、先给定的小量耳，若-1 ，那么|就是所求的根若，此时根据a继续下去，得到I的符号确定新的分割区间。，显然有，故方程的根的近似值为三、不动点迭代法i、不动点迭代法 将方程6-24 ）改写成等价的形式.II6-26）若要求 满足 I ，则 I;反之亦然，称为函数 E 的一个不动点。求的零点就等价于求E 的不动点，选择一个初始近似值71，将它代入6-26）右端，即可求得 I可以如此反复迭代计算丨， ，6-27 ）称为迭代函数。如果对任何:，由6-27 ）得到的序列冋有极限 二则称迭代方程 6-27 ）收敛，且 ! 为一的不动点，故称 6-27 ）为不动点迭代法。2、不动点的存在性与迭代法的收敛性1

4、）不动点的存在唯一性。定理1设 I 满足以下两个条件：1 ）对任意的 有 2 ）存在正常，使对任意 I 都有6-28 ）则一在一上存在唯一的不动点。在 E 的不动点存在唯一的情况下，可得到迭代法6-27 ）收敛的一个充分条件。定理2设 满足定理1中的两个条件，则对任意:，由6-27 ）得到的迭代序列叵收敛到的不动点，并有误差估计6-29)对定理1和定理2中的条件2），在使用时如果 1 且对任意 一 有_16-30 ）则由中值定理可知对 I 有它表明定理中的条件 2 ）可用6-30 ）代替。3、局部收敛性与收敛阶上面给出了迭代序列丽在区间上的收敛性，通常称为全局收敛性。有时不易检验定理的条件，实

5、际应用时通常只在不动点的邻近考察其收敛性，即局部收敛性。定义1设n有不动点，如果存在.：1的某个邻域:r ，对任意，迭代6-26 ）产生的序列收敛且收敛到，则称迭代法6-26 ）局部收敛。定理3设 E1有不动点L 在H的某个邻域连续，且Wf，则迭代法6-26）局部收敛。定义2设迭代过程：-：收敛于方程的根I，如果迭代误差I匕U满足下列渐进关系式则称该迭代过程是阶收敛的，特别地，匚|定理4对于迭代过程时称线性收敛，|_|时称超线性收敛,在所求根的邻近连续，并且I ，如果则该迭代过程在点邻近是阶收敛的。上述定理告诉我们，迭代过程的收敛速度依赖于迭代函数 E 的选取，如果当时 ，则该迭代过程只可能是

6、线性收敛。 6.3牛顿法及其收敛性1、牛顿法设已知方程厂I有近似根E 假定 ），将函数 -I在点展开，有于是方程可近似地表示为则可得近似计算公式为6-32)6-33)这就是牛顿Newton ）法。关于牛顿法6-33 ）的收敛性，可直接由定理 4得到，牛顿法在根的邻近是平方收敛的。又因|，可得6-34)2、简化牛顿法与牛顿下山法（1）简化牛顿法，也称平行弦法，其迭代公式为)6-35 )I ，I6-35迭代函数 1。若IT，即取 1 * 。在根呂附近成立，则迭代法在6-35 )中取 IX1局部收敛。，则称为简化牛顿法，这类方法计算量省，但只有线性收敛，其几何意义是用平行弦与轴交点作为的近似。2）牛

7、顿下山法。为了防止迭代发散，我们对迭代过程再附加一项要求，即具有单调性:6-36)满足这项要求的算法称为下山法。我们将牛顿法和下山法结合起来使用，即在下山法保证函数值稳定下降的前提下，用牛 顿法加快收敛速度。为此，我们将牛顿法的计算结果与前一步的近似值适当加权平均作为新的改进值6-37)其中称为下山因子,6-37 )即为6-38)称为牛顿下山法。选择下山因子时从亠 开始，逐次将减半进行试算，直到能使下降条件6-36 ）成立为止。 6.4弦截法与抛物线法用牛顿法求方程根，每步除计算外还要算二一.，而计算二一.往往较困难,为此下面介绍两种常用方法。1、弦截法设.=丨是 的近似根，我们利用I构造一次

8、插值多项式I_|并用 |的根作为 的新的近似根I。由于6-39 )因此有6-40）2、抛物线法 密勒Mu Iler）法）插值多项式0有两个零点J6-41）式中0为了从6-41 ）中定出一个值丨，我们需要讨论根式前正负号的取舍问题。在 三个近似根中，自然假定更接近所求的根，这时，为了保证精度，我们选式6-41 ）中较接近，：1的值作为新的近似根丨。为此，只要取根式前的符号与可的符号相同。 6.5孩子成长和学生考试成绩问题1）孩子成长问题一个男孩在11岁长到21岁过程中，身高的变化如表5-1所示，试找一个最佳的函数曲线来表示这个男孩的成长过程。见表5-1ti从11岁起年龄）0 0.81.42.02

9、.43.24.0增长高度hicm）0 0.742.255.258.2515.0021.38ti从11岁起年龄）4.8 5.46.07.08.0 10.0增长高度hicm）26.2528.8830.6032.25 33352）学生考试成绩问题:对8个学生测量其智商Iq和课后复习某门课时间t及该门课考试成绩g,得表5-2，试研究 该门课考试成绩与智商和课后复习时间之间的关系。见表5-2。智商Iq)105110120116122130 114102复习时间t)1012613168 2015考试成绩g)7579688591799876以上两个问题，均可以抽象为：通过实验等方法观测到反映某个函数.F 的

10、数据要求利用这些数据构造出f 的近似表达式 F。上一节介绍的插值法就是寻求近似函数的方法之一。但由于实验观测数据不可避免地带有误差，甚至是较大的误差， 所以使用插值法是不合适的，它会保留数据的误差。因此，不必要求近似函数L-1满足I，而只要求偏差按某种标准最小，以反映所给数据的总体趋势，消除局部波动的影响，这就是曲线拟合问题。 6.6估计水塔的水流量某居民区的民用自来水是由一个圆柱形的水塔提供，一般可以通过测量其水位来估计水 的流量。但问题是，当水塔水位下降到设定的最低水位时水泵自动启动向水塔供水，到设定 的最高水位时停止供水，在水泵自动加水期间无法测量水塔的水位和水泵的供水量。通常水 泵每天

11、供水两次，每次约两小时。水塔是一个高12.2M，直径17.4M的正圆柱，当水塔的水位降至最低水位，约8.2M时，水泵自动启动供水当水塔的水位升高到一个最高水位，约 10.8M时，水泵停止工作。表 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.

12、92510.95412.032水位8.5258.3888.220/10.82010.500时刻t12.95413.87514.98215.90316.82617.93119.037水位10.2109.9369.6539.4099.1808.9218.662时刻t19.95920.83922.01522.95823.88024.98625.908水位8.4338.220/10.59110.35410.180模型假设及符号说明1 ）模型假设假设水塔中流出的水流量只受社区的日常生活需要的影响，水的消耗每天大致差不多。 由Torricelli定律知，从水塔流出的最大流速正比于水位高度的平方根，题目中给

13、出水塔 的最高和最低水位分别为 10.82M和8.22M，所以对于这两种高度，最大水流速度的比约 为I，这表明我们可以假设水塔中水位对水流速度影响忽略不计。水泵工作时单位时间的供水量大致为常数，这个常数大于单位时间内从水塔中流出的水 流的最大流速，这是因为居民区内一直需要用水，不允许水塔中的水用光。水塔中水流量是时间的连续光滑函数，与水泵工作与否无关。这是因为虽然就个别用户 而言可能用水量有较大的变化，但由于个人的用水量与整个居民区用水量相比是非常小 的，从统计意义上来讲，不太可能同时整个社区的用水量增长或减少。水泵工作与否完全取决于水塔内水位的高度，且每次加水的工作时间为2小时，根据表6-1

14、中的数据可知，水泵第一次供水时间段为8.967,10.954，第二次供水时间段为20.839, 23.880。?符号说明：测量的时刻；：水位的高度；：水塔中水的体积；I亠：水塔中水流速度，即流量； ）:第i时段用水量，即没有供水时间段用水量；:第1和第2供水时间段用水量之和； ：一天中总用水量。2 ）问题分析问题要求任意时刻的用水率，即求单位时间流出的水的体积，一般称为水流速度或流量。由于水塔是一个圆柱体，体积口 |可以很容易地通过水位高度h计算出来，这样在水泵不工作的时间段，水流速度就可以从体积对时间的导数计算出来，由于没有水的体积 关于时间的函数表达式，而只能利用问题中给定的原始数据表6-

15、1，通过体积公式计算出离散的在测量时刻的体积 V,因此可以考虑用差商代替微商，也就是用离散代替连续的思想。为提高计算精度，采用二阶差商，即=。由于所有数据被水泵两次供水分割成三组数据对每组数据的中间数据采用中心差商，前 后两个数据不能采用中心差商，改用向前差商、向后差商或用中点公式进行差商。中点公式:以上分析了水泵不工作的时段，用水率的计算。对于水泵供水时段的用水率，计算难度较 大，我们只好用供水时间段前后的用水率进行插值或拟合而得到。有了任何时刻的用水率， 可以采用数值积分计算一天的总用水量。3）模型建立及求解通过以上对问题的分析，现在的问题已转化为根据某一天已测量的时刻水塔中水的流 速，产

16、生在整个区间24小时）上的函数或函数值，一般来说插值和拟合是两种最常用的方 法。分两步进行：计算水流速度并画出散点图计算水塔中水的体积 4）模型的建立通过对不同插值方法的比较，结合假设，考虑到流速应该是时间的连续光滑函数，所 以采用三次样条插值模型计算用水率函数。计算用水率首先用三次样条插值计算用水率函数回,MATLAB程序执行后，可得样条插值下的流速图计算一天的总用水量用三次样条插值模型得到的函数回 在时间区间0, 24上做数值积分可得一天的总用水量。程序执行后，可得一天的总用水量1266.1立方M。见图6-25）模型稳定性分析用不同时刻作为起始点，使用三次样条插值模型得到的用水率函数，在 长度为24小时区间上进行数值积分，计算程序为：执行文件后，所得结果列入表6-2。单位：立方M）从表6-2易见，选择不冋的起始点，所得结果比较稳定。=618.25第三时间段，即时的实际用水量为：237.787 X （10.591-10.18=97.73 其次用三次样条插值模型得到的函数回 分别在水泵未工作的三段时间内作数值积分，执行程序后，可得水泵未工作的三段时间内，用模型计算出的用水量见表6-3 ）。分三段的实际用水量与模型用水量的比较列成表6-3。并且从表6-3