Matlab求解非线性方程组及在潮流计算中的运用.PDFVIP

Matlab 求解非线性方程组及在潮流计算中的运用 引言 潮流计算是研究电力系统稳态运行情况的一种基本电气计算，常规潮流计算 的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上 的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是 电力系统稳定计算和故障分析的基础。 一 潮流计算算法原理：牛顿—拉夫逊法的基本原理 牛顿-拉夫逊法是一种求解非线性方程的数值解法，由于便于编写程序用计算 机求解，应用较广。下面以一元非线性代数方程的求解为例，来说明牛顿-拉夫 逊法的基本思想。 设欲求解的非线性代数方程为 f(x)=o 设方程的真实解为x* ，则必有f(x*)=0 。用牛顿-拉夫逊法求方程真实解x* 的步骤如下： 首先选取余割合适的初始估值x°作为方程f(x)=0 的解，若恰巧有f(x°)=0 ， 则方程的真实解即为x*= x°若f(x°)_0 ，则做下一步。 取x¹=x°+Äx°为第一次的修正估值，则 f(x¹)=f(x°+Äx°) 其中Äx°为初始估值的增量，即Äx°=x¹-x° 。设函数f(x)具有任意阶导数， 即可将上式在x° 的邻域展开为泰勒级数，即： f(x¹)=f(x°+Äx°)=f(x°)+f(x°)Äx°+[f(x°)(Äx°)2]/2+… 若所取的|Äx°|足够小，则含(Äx°)²的项及其余的一切高阶项均可略去，并使 其等于零，即： f(x¹)_f(x°)+f(x°)Äx°=0 故得 Äx°=-f(x°)/f(x°) 从而 x¹= x°-f(x°)/f(x°) 可见，只要f(x°)_0 ，即可根据上式求出第一次的修正估值x¹ ，若恰巧有 f(x¹)=0 ，则方程的真实解即为x*=x¹ 。若f(x¹)_0 ，则用上述方法由x¹再确定 第二次的修正估值x² 。如此反复叠代下去，直到求得真实解x*为止。 二 节点电压用极坐标牛顿—拉夫逊法潮流计算 节点的功率方程写成 其中 ij  i j 式中, ij  i j是两节点电压的相角差。 把节点功率表示为节点电压的幅值和相角的函数。在有n 个节点的系统 中，假定第1～m 号节点为节点，第m+1～n-1 号节点为节点，第n 号节点 为平衡节点。Un 和θn 是给定的，PV 节点的电压幅值Um+1～Un-1 也是 给定的。因此，只剩下n-1 个节点的电压相角θ1，θ2，⋯，θn-1 和m 个节点的电压幅值U1，U2，⋯，Um 是未知量。 实际上，对于每一个PQ 节点或每一个PV 节点都可以列写一个有功功 率不平衡量方程式，一共包含了n-1+m 个方程式，正好同未知数的数目相 同，而比直角坐标形式的方程少了n-1-m 个。 对于方程式可以写出修正方程式如下：  H 是(n-1) ×(n-1)阶方阵，其元素为 ；N 是(n-1) ×m 阶矩阵，其元素为 ； K 是m ×(n-1)阶矩阵，其元素为 ；L 是m ×m 阶方阵，其元素为 。 在这里把节点不平衡功率对节点电压幅值的偏导数都乘以该节点电压， 相应地把节点电压的修正量都除以该节点的电压幅值，这样，雅可比矩阵 的表达式就具有比较整齐的形式。 三 matlab运算过程的具体实例 例：网络接线如图，各支路阻抗和各节点功率均以标幺值标于图中，其 中节点2 连接的实际是发额定功率的发电厂，设节点1 的电压保持为1.06， 用牛顿-拉夫逊法计算系统中的潮流分布。  四 程序运行结果： 最终的潮流分布如下图： 五 小结 在计算机潮流计算中，使用牛顿-拉夫逊法进行潮流计算，非常快捷。在程序 编写过程中，对公式的由来应该非常清楚 ，不然在编写过程中易出错，比如for 循环的次数，一些细节问题都应该反复琢磨。在程序调试过程