牛顿迭代法在架空线路中的应用.doc

牛顿迭代法在架空线路中的应用 牛顿迭代法（Newtons method）又称为牛顿-拉夫逊方法（Newton-Raphson method），它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。多数方程不存在求根公式，因此求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根,此时线性收敛，但是可通过一些方法变成超线性收敛。另外该方法广泛用于计算机编程中。。而牛顿迭代法是用一种间接而特殊的方法来确定的。下面具体推到牛顿迭代公式。 假设是非线性方程为的一个近似根，把在处作泰勒展开： 若取前两项来近似代替（称为的线性化），则得近似的线性方程 设，令其解为，则得 （1） 这称为的牛顿迭代公式。 它对应的迭代方程为显然是的同解方程，故其迭代函数为 （） 在的根α的某个邻域内， 在α的邻域R内，对任意初值x0，应用由公式(1)来解方程的方法称为牛顿迭代法，它是解代数方程和超越方程的有效方法之一。 二、牛顿迭代法的几何意义 由(1)式知是点处的切线与x轴的交点的横坐标（如图1）。也就是说，新的近似值是用代替曲线的切线与x轴相交得到的。继续取点，再做切线与x轴相交，又可得，……。由图可见，只要初值取的充分靠近α，这个序列就会很快收敛于α。 由此可知，牛顿迭代法的集合意义是用一系列的切线与x轴的交点来逐步逼近曲线与x轴的交点。因此，牛顿迭代法也称切线法。 图1 牛顿迭代法原理 三、牛顿迭代法的步骤 ①准备。选定初始近似值及允许误差ε，计算，； ②迭代。按公式迭代一次，得到新的近似值，计算，； ③控制。如果满足，则终止迭代,以作为所求的根；否则转到步骤④； ④修改。如果迭代次数达到预定指定的次数N，或者则方法失败；否则以代替转步骤②继续迭代。 四、工程实例 1、工程概况 在架空输电线路的设计中不同气象条件下架空路线的弧垂、应力和线长占有十分重要的位置，是输电线路力学研究的主要内容。这是因为架空线的弧垂和应力直接影响着线路的正常安全运行，而架空线的微小变化和误差都会引起弧垂和应力相当大的改变。设计弧垂小，架空线的拉应力就大，震动现象加剧，安全系数小，同时杆塔载荷增大因而要求强度提高。设计弧垂过大，满足对地安全距离所需杆塔高度增加，线路投资增大，而且架空线的风摆、舞动和跳跃会造成线路停电事故，若加大塔头尺寸，必然会使投资再度提高。因此，设计合适的弧垂是十分重要的。 架空线的线长和弧垂是档距、高差和架空比载、应力的函数。当气象条件发生变化时，这些参数将会发生变化。气温的升降引起架空线的热胀冷缩，是线长、弧垂、应力发生相应变化。大风时的覆冰架空线上比载增加，应力增大由于弹性变形使架空线线长增加。不同气象条件下架空线的各参数之间存在着一定的关系。为了使用方便，常将各种气象条件下的架空线的应力和有关弧垂档距的变化用曲线表示出来，这种曲线称之为应力弧垂曲线，亦称力学特性曲线。 架空输电线路应力弧垂曲线为线路杆塔设计、计算、施工提供技术资料，针对传统的绘制方法存在效率低、质量差等缺点，通过对应力和弧垂的计算及过程分析，考虑了控制气象条件的主要因素，根据架空线状态方程的特点，采用牛顿迭代法，应用Matlab语言编制了架空线应力和弧垂的计算程序，依据该程序，输入线路的基本参数和迭代初值即可迅速而准确地计算出应力和弧垂值，结合实例绘制了应力弧垂曲线，结果表明该方法可以极大地提高应力和弧垂的计算精度，有效克服认为误差，简化线路设计工作，提高线路设计效率和质量。 2、应力弧垂曲线 2.1架空线的状态方程 当气温、载荷变化时，支撑与两悬挂点间的架空线应力可按状态方程求解。悬挂点等高时孤立档的状态方程为： （2） 令A=， B= 则式（2）为： （3） 式中：E为架空线的弹性系数，γ为导线半径，l为档距，t为气温。 式（3）为悬挂点等高时孤立档的状态方程。对于计算等高连续档得架空应力时，只要将孤立的档距用代表档距表示，状态方程就具有相同的形式；同样，对于悬挂点不等高时，只要引进一些相应系数，就可以表示与式（3）同样的结构形式。解出状态方程即可求出应力，进而求出弧垂。 2.2计算各气象条件下的应力和弧垂 求解不含一次项的一元三次方程，本文采用牛顿迭代法。令，其导数为： ，则牛顿迭代式为：