基于MATLAB实现的牛顿拉夫逊法的 IEEE33巴士辐射状配电系统负荷潮流程序知道总线数据和加载数据的任何系统+使用说明文档

basemva = 100; accuracy = 0.001; maxiter = 1000; % IEEE 33-BUS TEST SYSTEM % Bus Bus Voltage Angle ----Load---- -------Generator----- Injected % No code Mag. Degree MW Mvar MW Mvar Qmin Qmax Mvar busdata=[1 1 1.0 0.0 0.0 0.0 0.0 0.0 0 0 0 2 0 1.0 0.0 0.10000 0.06000 0.0 0.0 0 0 0 3 0 1.0 0.0 0.09000 0.04000 0.0 0.0 0 0 0 4 0 1.0 0.0 0.12000 0.08000 0.0 0.0 0 0 0 5 0 1.0 0.0 0.06000 0.03000 0.0 0.0 0 0 0 6 0 1.0 0.0 0.06000 0.02000 0.0 0.0 0 0 0 7 0 1.0 0.0 0.20000 0.10000 0.0 0.0 0 0 0 8 0 1.0 0.0 0.20000 0.10000 0.0 0.0 0 0 0 9 0 1.0 0.0 0.06000 0.02000 0.0 0.0 0 0 0 10 0 1.0 0.0 0.06000 0.02000 0.0 0.0 0 0 0 11 0 1.0 0.0 0.04500 0.03000 0.0 0.0 0 0 0 12 0 1.0 0.0 0.06000 0.03500 0.0 0.0 0 0 0 13 0 1.0 0.0 0.06000 0.03500 0.0 0.0 0 0 0 14 0 1.0 0.0 0.12000 0.08000 0.0 0.0 0 0 0 15 0 1.0 0.0 0.06000 0.01000 0.0 0.0 0 0 0 16 0 1.0 0.0 0.06000 0.02000 0.0 0.0 0 0 0 17 0 1.0 0.0 0.06000 0.02000 0.0 0.0 0 0 0 18 0 1.0 0.0 0.09000 0.04000 0.0 0.0 0 0 0 19 0 1.0 0.0 0.09000 0.04000 0.0 0.0 0 0 0 20 0 1.0 0.0 0.09000 0.04000 0.0 0.0 0 0 0 21 0 1.0 0.0 0.09000 0.04000 0.0 0.0 0 0 0 22 0 1.0 0.0 0.09000 0.04000 0.0 0.0 0 0 0 23 0 1.0 0.0 0.09000 0.05000 0.0 0.0 0 0 0 24 0 1.0 0.0 0.42000 0.20000 0.0 0.0 0 0 0 25 0 1.0 0.0 0.42000 0.20000 0.0 0.0 0 0 0 26 0 1.0 0.0 0.06000 0.02500 0.0 0.0 0 0 0 27 0 1.0 0.0 0.06000 0.02500 0.0 0.0 0 0 0 28 0 1.0 0.0 0.06000 0.02000 0.0 0.0 0 0 0 29 0 1.0 0.0 0.12000 0.07000 0.0 0.0 0 0 0 30 0 1.0 0.0 0.20000 0.60000 0.0 0.0 0 0 0 31 0 1.0 0.0 0.15000 0.07000 0.0 0.0 0 0 0 32 0 1.0 0.0 0.21000 0.10000 0.0 0.0 0 0 0 33 0 1.0 0.0 0.06000 0.04000 0.0 0.0 0 0 0]; % Line code % Bus bus R X 1/2 B = 1 for lines % nl nr p.u. p.u. p.u. > 1 or < 1 tr. tap at bus nl linedata=[1 2 0.057500 0.029300 0 1 2 3 0.307600 0.156600 0 1 3 4 0.228400 0.116300 0 1 4 5 0.237800 0.121100 0 1 5 6 0.511000 0.441100 0 1 6 7 0.116800 0.386100 0 1 7 8 1.068000 0.771000 0 1 8 9 0.642600 0.461700 0 1 9 10 0.651400 0.461700 0 1 10 11 0.122700 0.040600 0 1 11 12 0.233600 0.077200 0 1 12 13 0.915900 0.720600 0 1 13 14 0.337900 0.444800 0 1 14 15 0.368700 0.328200 0 1 15 16 0.465600 0.340000 0 1 16 17 0.804200 1.073800 0 1 17 18 0.456700 0.358100 0 1 2 19 0.102300 0.097600 0 1 19 20 0.938500 0.845700 0 1 20 21 0.255500 0.298500 0 1 21 22 0.442300 0.584800 0 1 3 23 0.281500 0.192400 0 1 23 24 0.560300 0.442400 0 1 24 25 0.559000 0.437400 0 1 6 26 0.126700 0.064500 0 1 26 27 0.177300 0.090300 0 1 27 28 0.660700 0.582600 0 1 28 29 0.501800 0.437100 0 1 29 30 0.316600 0.161300 0 1 30 31 0.608000 0.600800 0 1 31 32 0.193700 0.225800 0 1 32 33 0.212800 0.330800 0 1]; % lfybus % form the bus admittance matrix % lfnewton % Load flow solution by Newton-Raphson method % busout % Prints the power flow solution on the screen % lineflow % computes and displays the line flow and losses