电力系统潮流计算的MATLAB辅助程序设计潮流计算程序Word文件下载.docx

1、节点信息输入格式为单行输入，输入的数据形成一个矩阵，叫做矩阵。第一列为节点号；第二列为节点类型；第三列和第四列分别为节点电压幅值（标幺值）和相角（单位为度）；第五列和第六列分别为负荷的有功功率和无功功率；第七列到十列分别为发电机的有功功率、无功功率、最小无功出力和最大无功出力；最后一列为并联电容器注入无功功率。第二列的编码用0、1、2来区分节点、平衡节点和节点：0表示节点，输入正的有功功率（）和无功功率（），并且要设定节点电压初始估计值，一般幅值和相角分别设为1和0，若已经给定初始值，则用其给定值来代替1和0。1表示平衡节点，且已知该节点的电压幅值和相角。2表示节点，要设定该节点的节点电压幅值

2、和发电机的有功功率（），并设定发电机的无功最小出力和最大出力（）。2.线路数据文件 线路数据用节点对的方法来确定，数据包含在称为的矩阵中。第一列和第二列为节点号码，第三列到第五列为线路电阻、电抗及该线路电纳值的一半，以标幺值表示。最后一列为变压器分接头设定值，对线路来说，需要输入1。线路输入为无输入顺序，对变压器来说，左侧的节点号设为分接头端。3是线路数据输入变量，包括四项，前两项是节点编号，后两项是线路电阻和电抗，均以标幺值表示，函数返回节点导纳矩阵。三、潮流计算的程序清单1. 程序清单 % 0; 0; =0; = （:,1）; ; 1（k,1）;（n）（k,2）; （n）（k,3）; （n

3、）（k, 4）;（n）（k,5）; （n）（k,6）; （n）（k,7）; （n） = （k,8）;（n）（k, 9）; （n）（k, 10）;（n）（k, 11）; （n） = & = 1; n = 1; = 0*0; L = 1; （L） n & （L） 1）, （L）; = + （）*V（k）; （L）1）, （L）; = （V（n）*（）*V（n） + ） ; = （）; （n） = P（n） - （）; （n） = Q（n） - （）; （n） 1 S（n） ; P（n） = （）; Q（n） = （）; （n） =0; （n）=0; （n） = V（n）; （n） 2 （n） 0

4、= Q（n）* + （n） - （n）; （n） = .005 &= 10= 0.045 （n）, （n） = （n） + 0.005; （n） = （n）+.005; （n） = （n） - 0.005; （n）（n）+.005; （n） = （S（n）（V（n） - ）/ （）; V（n） = V（n） + *（n）（n）; = （n）; = （n）2 - 2）; （n） = + j*; V（n） = V（n） + *（n） （n）; （ （）, （） ）; &: ）, ）, （ .nn n , , ）; n = 1 （n） = （V（n）; （n） = （V（n）*180; S（n）（n

5、）*Q（n）; （n） = P（n）* + （n）; （n） = Q（n）* + （n） - （n）; 1; （k）（n）;（n） = （n）- j*（n）*（n）/（*（n）2）;（:,3）; （:,4） L 2程序清单（-1）; i = （-1）; = （:,1）;,2）; R = （:,3）;X = （:,4）; = j*（:,5）; a = （:, 6）;（: = （）, （）;Z = R + j*X; （,1）; %支路导纳 a（n） = 0 a（n） = 1;（）; % 将初始化为0 %非对角元素的数值 （k）（k）（k）（k）（k）（k）; （k）（k）（k）（k）;% 对角元素

6、的数值 （k） （） = （）（k）/（a（k）2） + （k）; （） = （）（k） （k）;3. 程序清单% a % .（）（ = n, ） . = nn = . . ; %5g, n）, （ %7.3f, （n）, %8.3f, （n）, （ %9.3f %9.3f %8.3fn, （n）, （, ）, %9.3fnn4程序清单n nn & 0%6g, P（n）*）%9.3f, Q（n）*）, （%9.3fn, （S（n）*） （L） k = （L）; = （V（n） - a（L）*V（k）*y（L）（L）2 + （L）（L）2*V（n）; = （V（k） - V（n）（L）*y（L）

7、 + （L）*V（k）; = V（n）*（）*; = V（k）*（）*; = + ; = （V（n） - V（k）（L）*y（L） + （L）*V（n）; = （V（k） - a（L）*V（n）*y（L）（L）2 + （L）（L）2*V（k）; （L） | （L）%12g, k）, （）, （, （）, （）, （L） & a（L） 1, a（L） , （ = 2;5程序清单% （k） 1, = 1; （k） 2 = 1;（k） = ; t = （）;2*2*;%雅可比矩阵 A J A（）=0; %初始化雅可比矩阵, （n）;（n）（n）;J11=0; J22=0; J33=0; J44=0;

8、 （）1 （） n | （） n （） n , l = （）; J1111+ （n）*（l）*（）*（t（）- （n） + （l）; J3333+ （n）*（l）*（）*（t（）- （n） + （l）; （n）1 J2222+ （l）*（）*（t（）- （n） + （l）; J4444+ （l）*（）*（t（）- （n） + （l）; （n） 1 & （l） 1 = （l）（l）; = l （l）; % J1的非对角元素 A（, ） （n）*（l）*（）*（t（）- （n） + （l）; （l） 0 % J2的非对角元素 A（, ） （n）*（）*（t（）- （n） + （l） （n） 0 %

9、 J3的非对角元素 A（, ） （n）*（l）*（）*（t（）- （n）（l）; （n） 0 & （l） 0 % J4的非对角元素 = （n）2*（）*（t（）33; = （n）2*（）*（t（）11; （n） 1 P（n）; Q（n） = ; % P （n） 2 Q（n）;= 7 2 （n） = （n） + 0.01; （n） = （n） - 0.01 A（） = J11; % J1对角元素 （） = P（n）; A（） = 2*（n）*（）*（t（）22; % J2对角元素 A（）= J33; % J3对角元素 A（） 2*（n）*（）*（t（）44; % J4对角元素 （） = Q（n）

10、; （n）; （n）（n）; （n） = （n）（）; （n）（n）（）; （）;V = .*（）*.*（）;180*; S（n）= P（n）*Q（n）;6程序清单 （n） 1, = 1;（n） = ; （） 0 | （） 2 B1（）（）; （） 0 B2（）（）;B1（B1）; B2 = （B2）;% 开始迭代= % 检验不平衡功率 （） n, l = （）;= 20 % 1 & 6 = 10 % （n）, % . （n） % 0.05 （n）, % （n） = （n） - 0.005 % . （） = （P（n）（n）; （） = （Q（n）（n）;1200; （）（）; A; ; J1

11、1 J33 B1 B1 B2 B2 四、30节点电力系统计算实例潮流计算时，必须将前面的六个程序保存在目录下格式为的文件，然后在的命令窗口输入如下命令： = 1.8;% 30节点电力系统% 母线母线 电压 相角 负载 发电机 注入功率% 编号节点 幅值 角度 有功 无功 有功 无功 无功最小值 无功最大值 无功1 1 1.06 0.0 0.0 0.0 0.0 0.0 0 0 0 2 2 1.043 0.0 21.70 12.7 40.0 0.0 -40 50 0 3 0 1.0 0.0 2.4 1.2 0.0 0.0 0 0 0 4 0 1.06 0.0 7.6 1.6 0.0 0.0 0 0

12、 0 5 2 1.01 0.0 94.2 19.0 0.0 0.0 -40 40 0 6 0 1.0 0.0 0.0 0.0 0.0 0.0 0 0 0 7 0 1.0 0.0 22.8 10.9 0.0 0.0 0 0 0 8 2 1.01 0.0 30.0 30.0 0.0 0.0 -30 40 0 9 0 1.0 0.0 0.0 0.0 0.0 0.0 0 0 0 10 0 1.0 0.0 5.8 2.0 0.0 0.0 -6 24 19 11 2 1.082 0.0 0.0 0.0 0.0 0.0 0 0 0 12 0 1.0 0 11.2 7.5 0 0 0 0 0 13 2 1.0

13、71 0 0 0.0 0 0 -6 24 0 14 0 1 0 6.2 1.6 0 0 0 0 0 15 0 1 0 8.2 2.5 0 0 0 0 0 16 0 1 0 3.5 1.8 0 0 0 0 0 17 0 1 0 9.0 5.8 0 0 0 0 0 18 0 1 0 3.2 0.9 0 0 0 0 0 19 0 1 0 9.5 3.4 0 0 0 0 0 20 0 1 0 2.2 0.7 0 0 0 0 0 21 0 1 0 17.5 11.2 0 0 0 0 0 22 0 1 0 0 0.0 0 0 0 0 0 23 0 1 0 3.2 1.6 0 0 0 0 0 24 0 1

14、0 8.7 6.7 0 0 0 0 4.3 25 0 1 0 0 0.0 0 0 0 0 0 26 0 1 0 3.5 2.3 0 0 0 0 0 27 0 1 0 0 0.0 0 0 0 0 0 28 0 1 0 0 0.0 0 0 0 0 0 29 0 1 0 2.4 0.9 0 0 0 0 0 30 0 1 0 10.6 1.9 0 0 0 0 0;% 线路数据% R X 1/2 B 1 1 2 0.0192 0.0575 0.02640 1 1 3 0.0452 0.1852 0.02040 1 2 4 0.0570 0.1737 0.01840 1 3 4 0.0132 0.0379 0.00420 1 2 5 0.0472 0.1983 0.02090 1 2 6 0.0581 0.1763 0.01870 1 4 6 0.0119 0.0414 0.00450 1 5 7 0.0460 0.1160 0.01020 1 6 7 0.0267 0.0820 0.00850 1 6 8 0.0120 0.0420 0.00450 1 6 9 0.0 0.2080 0.0 0.978 6 10 0