《Matlab及C语言在潮流计算运用》由会员分享，可在线阅读，更多相关《Matlab及C语言在潮流计算运用（18页珍藏版）》请在金锄头文库上搜索。
1、Matlab及C语言在潮流计算中的运用（三峡大学电气信息学院20071096班）摘要：潮流计算是研究电力系统稳态运行情况的一种基本电气计算，常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。关键词：电力系统分析 潮流计算 牛顿-拉夫逊法 C语言 Matlab一，潮流计算算法原理：牛顿拉夫逊法的基本原理 牛顿-拉夫逊法是一种求解非线性方程的数值解法，由于便于编写程序用计算机求解，应用较广。下面以一元非线性代数方程的求解为例，来说明牛顿-拉夫逊法的基本思想。 设欲
2、求解的非线性代数方程为 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
3、)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*为止。二，节点电压用极坐标牛顿拉夫逊法潮流计算 节点的功率方程写成 其中式中，是两节点电压的相角差。把节点功率表示为节点电压的幅值和相角的函数。在有n个节点的系统中，假定第1m号节点为节点，第m+1n-1号节点为节点，第n号节点为平衡节点。Un和n是给定的，PV节点的电压幅值Um+1Un-1也是给定的。因此，只剩下n-1个节
4、点的电压相角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是mm阶方阵，其元素为 。在这里把节点不平衡功率对节点电压幅值的偏导数都乘以该节点电压，相应地把节点电压的修正量都除以该节点的电压幅值，这样，雅可比矩阵的表达式就具有比较整齐的形式。三，matlab及c
5、语言运算过程的具体实例例：网络接线如图，各支路阻抗和各节点功率均以标幺值标于图中，其中节点2连接的实际是发额定功率的发电厂，设节点1 的电压保持为1.06，用牛顿-拉夫逊法计算系统中的潮流分布。0.02+j0.060.45+j0.150.08+j0.240.01+j0.030.4+j0.050.08+j0.240.04+j0.120.06+j0.180.06+j0.180.2+j0.2GG52134-0.02+j0.16四，程序运行结果：Matlab程序运行结果如下：Y = 3.7500 -11.2500i -2.5000 + 7.5000i 0 -1.2500 + 3.7500i 0 -2.
6、5000 + 7.5000i 10.8333 -32.5000i -1.6667 + 5.0000i -1.6667 + 5.0000i -5.0000 +15.0000i 0 -1.6667 + 5.0000i 12.9167 -38.7500i -10.0000 +30.0000i -1.2500 + 3.7500i -1.2500 + 3.7500i -1.6667 + 5.0000i -10.0000 +30.0000i 12.9167 -38.7500i 0 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 6.2500 -18.7500icount =
7、 4U = 1.0555 1.0572 1.0319 1.0331 1.0600Angle = 0.0417 0.0055 -0.0261 -0.0236 0Error_Angle_U = 1.0e-008 * -0.5980 -0.0366 0.1334 0.1379 -0.8207 -0.4688 -0.5584 -0.6216S = 0.6000 - 0.1000i 0.2000 + 0.2000i -0.4500 - 0.1500i -0.4000 - 0.0500i 0.0677 + 0.1531iLine_S = 0 0.2905 - 0.1333i 0 0.2991 - 0.01
8、63i 0 -0.2956 + 0.1229i 0 0.2185 + 0.0764i 0.2028 + 0.0746i 0.0764 - 0.0760i 0 -0.2109 - 0.0812i 0 -0.0918 - 0.0132i -0.1389 - 0.0789i -0.2931 - 0.0031i -0.1960 - 0.0786i 0.0919 + 0.0130i 0 0 0 -0.0771 + 0.0758i 0.1448 + 0.0773i 0 0 Line_I = 0 0.2697 + 0.1377i 0 0.2824 + 0.0273i 0 -0.2790 - 0.1178i
9、0 0.2071 - 0.0711i 0.1922 - 0.0695i 0.0719 + 0.0723i 0 -0.2022 + 0.0840i 0 -0.0886 + 0.0151i -0.1326 + 0.0799i -0.2836 + 0.0097i -0.1879 + 0.0806i 0.0886 - 0.0146i 0 0 0 -0.0727 - 0.0715i 0.1366 - 0.0729i 0 0 C语言运行结果如下：节点导纳矩阵：3.750000+j-11.250000 -2.500000+j7.500000 -0.000000+j-0.000000 -1.250000+j3
10、.750000 -0.000000+j-0.000000-2.500000+j7.500000 10.833333+j-32.500000 -1.666667+j5.000000 -1.666667+j5.000000 -5.000000+j15.000000-0.000000+j-0.000000 -1.666667+j5.000000 12.916667+j-38.750000 -10.000000+j30.000000 -1.250000+j3.750000-1.250000+j3.750000 -1.666667+j5.000000 -10.000000+j30.000000 12.9
11、16667+j-38.750000 -0.000000+j-0.000000-0.000000+j-0.000000 -5.000000+j15.000000 -1.250000+j3.750000 -0.000000+j-0.000000 6.250000+j-18.750000迭代次数：count=4节点相角分别为:=0.041710 0.005480 -0.026114 -0.023594 0.000000节点电压分别为:U=1.055525 1.057178 1.031931 1.033072 1.060000节点功率为:S=0.600000+j-0.100000 0.200000+j
12、0.200000 -0.450000+j-0.150000 -0.400000+j-0.050000 0.067700+j0.153099请按任意键继续. . .经过比较matlab和c语言运行的结果是一样的！最终的潮流分布如下图：0.0764 - 0.0760i0.45+j0.15-0.1389 - 0.0789i-0.0918 - 0.0132i0.4+j0.050.2905 - 0.1333i0.2028 + 0.0746i0.2+j0.2GG125340.2956 - 0.1229i-0.06+0.1i0.2931 + 0.0031i0.0677 + 0.1531i0.2185 + 0.0764i0.2991 - 0.0163i0.2109 + 0.0812i-0.0918 - 0.0132i-0.092 - 0.013i-0.0918 - 0.0132i0.0771 - 0.0758i-0.1448 - 0.0773i0.1960 + 0.0786i五，小结在计算机潮流计算中，使用牛顿-拉夫逊法进行
