牛顿拉夫逊法在潮流计算中的应用.docxVIP

牛顿-拉夫逊法在潮流计算中的应用简介牛顿迭代法又称为牛顿-拉夫逊方法，它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。多数方程不存在求根公式，因此求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要。方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。由于便于编写程序用计算机求解，应用较广。下面以一元非线性代数方程的求解为例，来说明牛顿-拉夫逊法的基本思想。设欲求解的非线性代数方程为　　　　　 ?f(x)=o设方程的真实解为x*，则必有f(x*)=0。用牛顿-拉夫逊法求方程真实解x*的步骤如下：首先选取余割合适的初始估值x°作为方程f(x)=0的解，若恰巧有f(x°)=0，则方程的真实解即为x*= x°若f(x°)≠0，则做下一步。取x1=x°+Δx°为第一次的修正估值，则f(x1)=f(x°+Δx°)其中Δx°为初始估值的增量，即Δx°=x1-x°。设函数f(x)具有任意阶导数，即可将上式在x°的邻域展开为泰勒级数，即：　f(x1)=f(x°+Δx°)=f(x°)+f(x°)Δx°+[f(x°)(Δx°)2]/2+…若所取的|Δx°|足够小，则含(Δx°)2的项及其余的一切高阶项均可略去，并使其等于零，即：f(x1)≈f(x°)+f(x°)Δx°=0Δx°=-f(x°)/f(x°)x1= x°-f(x°)/f(x°)可见，只要f(x°)≠0，即可根据上式求出第一次的修正估值x1，若恰巧有f(x1)=0，则方程的真实解即为x*=x1。若f(x1)≠0，则用上述方法由x1再确定第二次的修正估值x2。如此反复叠代下去，直到求得真实解x*为止。设为第k次的估值， 为第(k+1)次的修正估值，则有：　 此式为一叠代表达式，称为牛顿-拉夫逊算式。函数f(x)在的一阶导数值在实际中只要足够小，即满足：或迭代即可结束。式中ε为预先指定的一个小正数,视需要而定。需要指出的是，用牛顿-拉夫逊法求解非线性代数方程时，巧妙地设定一个合适的初始值是十分重要的，若初始值设定的合适，不但可以保证迭代收敛（即向真实值逼近），而且可使迭代次数减少。相反，若初始值设定的不合适，则不但要增加迭代次数，而且还可能迭代发散（即远离真实值）或者循环（即在真实值附近往复变化）。为解决此问题，可在叠代前先粗略画出函数f(x)是曲线，在该曲线与x轴的焦点附近设定合适的初始估值会有裨益。?牛顿一拉夫逊法潮流计算是最基本的潮流算法。计算是在三元素牛顿一拉夫逊法潮流的基础上的进一步简化。这种算法的特点是：方法简单，占用内存少，迭代过程中消元计算元素少，且不影响其收敛特性及计算结果的准确性。一、牛顿—拉夫逊法修正方程式简介假设系统中共有竹个节点，PQ节点数为m，PV节点数为n一俄一1，最后一个为平衡节点。直角坐标牛顿一拉夫逊法修正方程式可写成其雅可比矩阵的各元素定义分别为相应地，极坐标牛顿一拉夫逊法的修正方程式为雅可比矩阵的各元素定义分别为二、三元素牛顿一拉夫逊法潮流简介通过对雅可比矩阵元素的计算，我们发现对角元的计算是至关重要的。另一方面，由于假设前提的使用，各元素的计算也存在着很大的近似性。因此，可对其进行各种简化以得到各种派生算法。根据此观点和电力系统的结构，可假设各节点只与2个相邻节点连接。且在取节点号时可假设节点1仅和节点2，3连接，节点2仅和节点3，4连接，以此类推。在这种情况下，其它节点之间的连接被忽略。下面按直角坐标和极坐标两种方式分别予以讨论。采用直角坐标时，按上述假设，其雅可比矩阵可简化为对式(3)进行高斯消去时，只需对偶数行消元，且只需消去一个元素．这样，消元计算量将大大减少。如果再忽略对角元素J和R，式(3)可简化为式(3——1)．此时无须进行消元，可直接进行迭代计算。若忽略所有的N，J，R元素，可得到其进一步的简化形式。对上述3式迭代计算节点电压新值时，适当选择加速因子可使计算速度大大加快。a-加速因子三、当采用极坐标时，其雅可比矩阵可简化为对（4）式进行高斯消去时，也只需对偶数行消元。忽略对角元素的J、R，可得其进一步的简化形式，如式（4-1）若忽略所有的N，．厂。R元素，式(4—1)还可简化为同样，计算节点电压新值时，适当选择加速因子可使计算速度大大加快四、结束语根据以上计算分析可见，三元素牛顿一拉夫逊法潮流计算及其简化形式方法简单，占用内存少，迭代过程中计算量小，且不影响收敛特性及计算结果的准确性．虽然迭代次数有所增加，但计算速度仍可接受．对简化法二而言，若作进一步简化并引入塞德尔迭代方式，计算速度还可有较大程度的提高。此外，适当选择加速