水力学中常用的基本计算方法.docVIP

水力学中常用的基本计算方法 水力学中经常会遇到一些高次方程，微分方程的求解问题。多年来，求解复杂高次方程的基本方法便是试算法，或查图表法，对于简单的微分方程尚可以用积分求解，而边界条件较为复杂的微分方程的求解就存在着较大的困难，但随着计算数学的发展及计算机的广泛使用，一门新的水力学分支《计算水力学》应运而生，但用计算机解决水力学问题，还需要了解一些一般的计算方法。在水力学课程中常用的有以下几种，现分述于后。 一、高次方程式的求解方法： （一）二分法 1、二分法的基本内容：在区间[X1，X2]上有一单调连续函数F（x）=0，则可绘出F（x）～Ｘ关系曲线。如果在两端点处函数值异号即F（x１）·F（x2）＜0，（见图（一）），则方程F（x）=0，在区间[X1，X2]之间有实根存在，其根的范围大致如下：取 1°若F（x2）·F（x3）＞0， 则解ξ∈[X1，X3] 2°若F（x2）·F（x3）＜0， 则解ξ∈[X3，X2] 3°若F（x2）·F（x3）=0， 则解ξ=X3 对情况1°，可以令x2=x3，重复计算。 对情况2°，可以令x1=x3，重复计算。 当规定误差ε之后，只要|x1-x2|≤ε，则x1(或x2)就是方程F(x)=0的根。 显然，二分法的理论依据就是高等数学中的连续函数介值定理。 它的优点是思路清晰，计算简单，其收敛速度与公比为的等比级数相同；它的局限性在于只能求实根，而不能求重根。 2、二分法的程序框图(以求解明渠均匀流正常水深为例) 最后必须说明，二分法要求x2值必须足够大，要保证F1·F2＜0，否则计算得不到正确结果。为了避免x2值不够大，产生计算错误，在程序中加入了判别条件F1·F2＞0。也可以给定xJ及步长△x，让计算机选择x2(x2=x1+△x)。 (二)牛顿法， 1、牛顿法的基本内容：设有连续函数F(x)=0，则可以绘出F(x)～x关系曲线，选取初值xo，过点（xo·F(xo)）作一切线，其斜率为辅F＇（xo），切线与x轴的交点是x1， 则有： 再过(x1，F(x1)作切线，如此类推得到牛顿法的一个迭代序列：xn+l=xn-F(xn)/F＇(xn)，令xn=xn+1，重复计算，直至满足给定的精度要求，即|xn+1-xn|≤，从而得到方程F(x)=0的根。 牛顿法具有平方收敛速度，比较快，但计算工作量大，每次运算除计算函数值外，还要计算微商值。对于牛顿法来讲，只要F(x)在零点附近存在连续的二阶微商，ξ是F(x)的一重零点，且初值xo充分接近于ξ，那么牛顿迭代就一定收敛。 2、牛顿迭代法的程序图框可以参看临界水深计算程序。 (三)迭代法 1、迭代法的基本内容：设有一方程F（x）=0，可以写成等价形式x=φ(x)，再做迭代式xn+1=φ(xn),只要给定初始值xo，用式子 进行反复运算，就可以得到了个序列{xn}。 如果序列{xn}收敛，并假设ξ= xn+1=φ（xn）=φ(ξ)也就是说ξ是方程F(x)=0的根。一般可以根据问题的性质人为地规定精度ε，认为当|xn+1-xn|≤ε时，xn+1或xn就是方程的根。 从上面可以看出，迭代法的实质就是把方程F(x)=0的根问题转化为求迭代序列{xn}的极限问题，迭代过程的几何意义就是求 两曲线的交点问题，交点处的横坐标正好就是方程的解。但是，用计算机进行迭代计算，最关键的问题就是要构造出正确的迭代计算公式；如果构造的迭代计算公式不适当，或者初值给的不合适，则都可能导致计算不收敛或者发生计算错误。例如，泄水建筑物下游收缩水深的计算，其基本数学公式是： 若进行迭代计算选用如下公式： （hc1为初值） 则只要初值hc1∈[0，Eo]，则计算结果均收敛于正确解。如若选用迭代公式，。只要初值hc1＞0，则计算均收敛，但收敛于下游一个较大的淹没水深，并非正确解，发生计算错误。在非棱柱体明渠水面线计算中，同样存在这一问题。可见迭代公式选择的重要性。 迭代公式的收敛性可以用下述定理描述：把方程F(x)=0改写成x=φ(x)，如果φ(x)满足lipschitz条件：即对任意的x1和x2，都有|φ（x1）-φ（x2）|≤L|x1-x2|，其中L为一与x1和x2均无关的正常数(简称李氏常数)；若L＜1，则迭代收敛，收敛的准确解为ξ，其收敛速度为：|xn-ξ|≤ 高次方程的解往往不是一个，因此构造出迭代公式，并判定其收敛之后，尚要根据问题的背景及实际情况对收敛解进行分析讨论，确认为是真解方可，否则，尚需修改初值或重新构造迭代公式。 2、迭代法的程序框图可以参看水跃水力计算程序。 本程序集中迭代法应用较多，主要原因是迭代法思路简单明了，收敛较快，程序也短。 二、插值