基于牛顿拉夫逊法潮流计算的matlab实验报告

《基于牛顿拉夫逊法潮流计算的matlab实验报告》由会员分享，可在线阅读，更多相关《基于牛顿拉夫逊法潮流计算的matlab实验报告（26页珍藏版）》请在装配图网上搜索。

1、一、实验目的应用MATLAB语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。要求：（1）潮流计算方法为牛顿-拉夫逊法。（2）编程语言为MATLAB。（3）程序具有较强通用性。二、程序流程图所用公式 其中开始形成节点导纳矩阵输入原始数据设节点电压，i=1,2,n,is置迭代次数置节点号i=1按式（3-3），（3-4）计算雅克比矩阵元素按式（3-2）计算节点的，节点的，求解修正方程式，得，雅克比矩阵是否已全部形成？计算平衡节点及PV节点功率求，迭代次数 k=k+1i=i+1？潮流计算完成计算各节点电压的新值： 3、 求解问题及其结果 IEEE-美国新英格兰10机39节点测试系统一、 系统单线图

2、 二、系统参数1）系统容量基准值为100MVA。2) 负荷数据见表D-1表D-1 负荷 数据节点号有功（MW）无功（Mvar）节点号有功（MW）无功（Mvar）BUS- 3BUS- 4BUS- 7BUS- 8BUS- 12BUS- 15BUS- 16BUS- 18BUS- 20322.0500.0233.8522.08.5320.0329.0158.0680.02.4184.084.0176.088.0153.032.330.0103.0BUS- 21BUS- 23BUS- 24BUS- 25BUS- 26BUS- 27BUS- 28BUS- 29BUS- 39274.0247.5308.02

3、24.0139.0281.0206.0283.51104.0115.084.6-92.247.217.075.527.626.93）发电机数据见表D-2表D-2 发电机数据 发电机节点号有功（MW）无功（Mvar）电压（p.u.）PVPQPQPVPVPQPVPVPVBUS-30BUS-32BUS-33BUS-34BUS-35BUS-36BUS-37BUS-38BUS-39250.00650.00632.00508.00650.00560.00540.00830.001000.00175.90103.3596.881.047501.012301.049301.027801.026501.0300

4、0平衡节点BUS-310.0(电压幅角)1.0（幅值）4）线路参数见表D-3表D-3 线路参数序号节点I节点JR(p.u.)X(p.u.)B/2(p.u.)LN1LN2LN3LN4LN5LN6LN7LN8LN9LN10LN11LN12LN13LN14LN15LN16LN17LN18LN19LN20LN21LN22LN23LN24LN25LN26LN27LN28LN29LN30LN31LN32LN33LN34BUS-2BUS-39BUS-3BUS-25BUS-4BUS-18BUS-5BUS-14BUS-6BUS-8BUS-7BUS-11BUS-8BUS-9BUS-39BUS-11BUS-13BU

5、S-14BUS-15BUS-16BUS-17BUS-19BUS-21BUS-24BUS-18BUS-27BUS-22BUS-23BUS-24BUS-26BUS-27BUS-28BUS-29BUS-29BUS-1BUS-1BUS-2BUS-2BUS-3BUS-3BUS-4BUS-4BUS-5BUS-5BUS-6BUS-6BUS-7BUS-8BUS-9BUS-10BUS-10BUS-13BUS-14BUS-15BUS-16BUS-16BUS-16BUS-16BUS-17BUS-17BUS-21BUS-22BUS-23BUS-25BUS-26BUS-26BUS-26BUS-280.003500.00

6、1000.001300.007000.001300.001100.000800.000800.000200.000800.000600.000700.000400.002300.001000.000400.000400.000900.001800.000900.000700.001600.000800.000300.000700.001300.000800.000600.002200.003200.001400.004300.005700.001400.041100.025000.015100.008600.021300.013300.012800.012900.002600.011200.0

7、09200.008200.004600.036300.025000.004300.004300.010100.021700.009400.008900.019500.013500.005900.008200.017300.014000.009600.035000.032300.014700.047400.062500.015100.349350.37500.128600.073000.110700.106900.067100.069100.021700.073800.056500.069450.039000.190200.600000.036450.036450.086150.183000.0

8、85500.067100.152000.127400.034000.065950.160800.128250.092300.180500.256500.119800.390100.514500.12450LN35: BUS-4接有并联电容器，B41.0000LN36: BUS-5接有并联电容器，B42.00005)变压器参数见表D-4表D-4 变压器参数序号节点I节点JR(p.u.)X(p.u.)变比(%)TR37TR38TR39TR40TR41TR42TR43TR44TR45TR46TR47TR48BUS-11BUS-13BUS-30BUS-31BUS-32BUS-34BUS-33BUS-3

9、5BUS-36BUS-37BUS-38BUS-20BUS-12BUS-12BUS-2BUS-6BUS-10BUS-20BUS-19BUS-22BUS-23BUS-25BUS-29BUS-190.001600.001600.000000.000000.000000.000900.000700.000000.000500.000600.000800.000700.043500.043500.018100.025000.020000.018000.014200.014300.027200.023200.015600.01380100.6100.6102.5107.0107.0100.9107.010

10、2.5100.0102.5102.5106.0%IEEE-美国新英格兰10机39节点测试系统 % 1 2 3 4 5 6 % bus volt angle p q typebus= 1 1.0000 0.00 0.00 0.00 1 2 1.0000 0.00 0.00 0.00 1 3 1.0000 0.00 -3.22 -0.024 1 4 1.0000 0.00 -5.00 -1.84 1 5 1.0000 0.00 0.00 0.00 1 6 1.0000 0.00 0.00 0.00 1 7 1.0000 0.00 -2.338 -0.84 1 8 1.0000 0.00 -5.22

11、 -1.76 1 9 1.0000 0.00 0.00 0.00 1 10 1.0000 0.00 0.00 0.00 1 11 1.0000 0.00 0.00 0.00 1 12 1.0000 0.00 -0.085 -0.88 1 13 1.0000 0.00 0.00 0.00 1 14 1.0000 0.00 0.00 0.00 1 15 1.0000 0.00 -3.20 -1.53 1 16 1.0000 0.00 -3.29 -0.323 1 17 1.0000 0.00 0.00 0.00 1 18 1.0000 0.00 -1.58 -0.30 1 19 1.0000 0.

12、00 0.00 0.00 1 20 1.0000 0.00 -6.80 -1.03 1 21 1.0000 0.00 -2.74 -1.15 1 22 1.0000 0.00 0.00 0.00 1 23 1.0000 0.00 -2.475 -1.15 1 24 1.0000 0.00 -3.08 -0.922 1 25 1.0000 0.00 -2.24 -0.472 1 26 1.0000 0.00 -1.39 -0.17 1 27 1.0000 0.00 -2.81 -0.755 1 28 1.0000 0.00 -2.06 -0.276 1 29 1.0000 0.00 -2.835

13、 -0.269 1 30 1.0475 0.00 2.50 0.00 2 31 1.0000 0.00 0.00 0.00 3 32 1.0000 0.00 6.50 1.759 1 33 1.0000 0.00 6.32 1.0335 1 34 1.0123 0.00 5.08 0.00 2 35 1.0493 0.00 6.50 0.00 2 36 1.0000 0.00 5.60 0.9688 1 37 1.0278 0.00 5.40 0.00 2 38 1.0265 0.00 8.30 0.00 2 39 1.0300 0.00 -1.04 0.00 2; % 1 2 3 4 5 6

14、 7% line: from bus to bus R, X, G, B/2 Kline= 2 10.00350 0.04110 0 0.34935 0; 39 10.00100 0.02500 0 0.37500 0; 3 20.00130 0.01510 0 0.12860 0; 25 20.00700 0.00860 0 0.07300 0; 4 30.00130 0.02130 0 0.11070 0; 18 30.00110 0.01330 0 0.10690 0; 5 40.00080 0.01280 0 0.06710 0; 14 40.00080 0.01290 0 0.069

15、10 0; 6 50.00020 0.00260 0 0.02170 0; 8 50.00080 0.01120 0 0.07380 0; 7 60.00060 0.00920 0 0.05650 0; 11 60.00070 0.00820 0 0.06945 0; 8 70.00040 0.00460 0 0.03900 0; 9 80.00230 0.03630 0 0.19020 0; 39 90.00100 0.02500 0 0.60000 0; 11 100.00040 0.00430 0 0.03645 0; 13 100.00040 0.00430 0 0.03645 0;

16、14 130.00090 0.01010 0 0.08615 0; 15 140.00180 0.02170 0 0.18300 0; 16 150.00090 0.00940 0 0.08550 0; 17 160.00070 0.00890 0 0.06710 0; 19 160.00160 0.01950 0 0.15200 0; 21 160.00080 0.01350 0 0.12740 0; 24 160.00030 0.00590 0 0.03400 0; 18 170.00070 0.00820 0 0.06595 0; 27 170.00130 0.01730 0 0.160

17、80 0; 22 210.00080 0.01400 0 0.12825 0; 23 220.00060 0.00960 0 0.09230 0; 24 230.00220 0.03500 0 0.18050 0; 26 250.00320 0.03230 0 0.25650 0; 27 260.00140 0.01470 0 0.11980 0; 28 260.00430 0.04740 0 0.39010 0; 29 260.00570 0.06250 0 0.51450 0; 29 280.00140 0.01510 0 0.12450 0; 4 0 0 0 0 1.0000 0; 5

18、0 0 0 0 2.0000 0; 11 120.00160 0.04350 0 0 100.60000/100; 13 120.00160 0.04350 0 0 100.60000/100; 30 20.00000 0.01810 0 0 102.50000/100 ; 31 60.00000 0.02500 0 0 107.00000/100 ; 32 100.00000 0.02000 0 0 107.00000/100 ; 34 200.00090 0.01800 0 0 100.90000/100 ; 33 190.00070 0.01420 0 0 107.00000/100 ;

19、 35 220.00000 0.01430 0 0 102.50000/100 ; 36 230.00050 0.02720 0 0 100.00000/100 ; 37 250.00060 0.02320 0 0 102.50000/100 ; 38 290.00080 0.01560 0 0 102.50000/100 ; 20 190.00070 0.01380 0 0 106.00000/100 ;计算结果牛顿拉夫逊法潮流计算结果 节点计算结果： n节点 节点电压 节点相角（角度） 节点注入功率 1 1.049185 -8.874991 0.000000 + j 0.000000 2

20、1.053167 -6.367180 0.000000 + j 0.000000 3 1.041493 -9.207297 -3.220000 + j -0.024000 4 1.036574 -10.042585 -5.000000 + j -1.840000 5 1.044652 -8.959237 0.000000 + j 0.000000 6 1.043883 -8.293104 0.000000 + j 0.000000 7 1.032645 -10.342431 -2.338000 + j -0.840000 8 1.031177 -10.811816 -5.220000 + j

21、-1.760000 9 1.042715 -10.595648 0.000000 + j 0.00000010 1.046426 -6.010476 0.000000 + j 0.00000011 1.044322 -6.792462 0.000000 + j 0.00000012 1.030736 -6.795388 -0.085000 + j -0.88000013 1.042351 -6.675491 0.000000 + j 0.00000014 1.036310 -8.232337 0.000000 + j 0.00000015 1.018517 -8.519794 -3.20000

22、0 + j -1.53000016 1.025492 -7.051856 -3.290000 + j -0.32300017 1.032750 -8.077118 0.000000 + j 0.00000018 1.034779 -8.936485 -1.580000 + j -0.30000019 1.044862 -2.382169 0.000000 + j 0.00000020 0.988148 -3.811032 -6.800000 + j -1.03000021 1.024926 -4.596980 -2.740000 + j -1.15000022 1.042650 -0.0705

23、12 0.000000 + j 0.00000023 1.032952 -0.245457 -2.475000 + j -1.15000024 1.021125 -6.906503 -3.080000 + j -0.92200025 1.060163 -4.952002 -2.240000 + j -0.47200026 1.052697 -6.205207 -1.390000 + j -0.17000027 1.037683 -8.217337 -2.810000 + j -0.75500028 1.050444 -2.695196 -2.060000 + j -0.27600029 1.0

24、50163 0.063077 -2.835000 + j -0.26900030 1.004392 1.594781 6.500000 + j 1.75900031 0.991632 2.892572 6.320000 + j 1.03350032 1.050539 7.797786 5.600000 + j 0.96880033 1.047500 -3.957598 2.500000 + j 1.21117434 1.012300 1.385774 5.080000 + j 1.82635935 1.049300 4.925324 6.500000 + j 2.63756636 1.0278

25、00 1.819476 5.400000 + j -0.10822437 1.026500 7.125579 8.300000 + j 0.21422538 1.030000 -10.390696 -1.040000 + j -2.29163939 1.000000 0.000000 5.628660 + j 1.384403线路计算结果：n节点I 节点J 线路功率S(I,J) 线路功率S(J,I) 线路损耗dS(I,J) 2 1 1.178698 + j -0.360055 -1.174311 + j -0.360481 0.004386 + j -0.72053639 1 6.405845

26、 + j -2.096152 -6.361848 + j 2.408287 0.043997 + j 0.312135 3 2 -3.633961 + j -0.542613 3.649983 + j 0.446577 0.016021 + j -0.09603625 2 2.370242 + j -1.109311 -2.328681 + j 0.997356 0.041562 + j -0.111955 4 3 -0.750370 + j -0.307172 0.751094 + j 0.080014 0.000724 + j -0.22715918 3 0.337560 + j -0.6

27、63855 -0.337133 + j 0.438599 0.000427 + j -0.225256 5 4 1.635254 + j 0.499000 -1.633054 + j -0.609119 0.002200 + j -0.11011914 4 2.621711 + j -0.216428 -2.616576 + j 0.150777 0.005135 + j -0.065651 6 5 4.826035 + j -0.675350 -4.821682 + j 0.684607 0.004353 + j 0.009257 8 5 -3.178130 + j -1.041836 3.

28、186428 + j 0.998989 0.008297 + j -0.042847 7 6 -4.249274 + j -0.969559 4.259899 + j 1.010657 0.010625 + j 0.04109811 6 3.465003 + j -0.270003 -3.457273 + j 0.209136 0.007730 + j -0.060866 8 7 -1.909893 + j -0.196732 1.911274 + j 0.129559 0.001381 + j -0.067173 9 8 0.132235 + j 0.116464 -0.131977 + j

29、 -0.521432 0.000258 + j -0.40496839 9 7.617154 + j -1.902126 -7.557438 + j 2.142687 0.059717 + j 0.24056111 10 -3.483660 + j -0.203064 3.488121 + j 0.171352 0.004461 + j -0.03171213 10 -3.008372 + j -0.730489 3.011879 + j 0.688680 0.003508 + j -0.04180914 13 -2.934129 + j -0.411463 2.941429 + j 0.30

30、7264 0.007300 + j -0.10419915 14 -0.311115 + j -0.998556 0.312417 + j 0.627891 0.001303 + j -0.37066516 15 2.896296 + j 0.430232 -2.888885 + j -0.531444 0.007411 + j -0.10121217 16 -2.048841 + j 0.950740 2.052282 + j -1.049122 0.003441 + j -0.09838319 16 4.542969 + j 0.681545 -4.511670 + j -0.625873

31、 0.031300 + j 0.05567221 16 3.324778 + j -0.302389 -3.316338 + j 0.177006 0.008440 + j -0.12538324 16 0.410793 + j -0.811601 -0.410571 + j 0.744757 0.000222 + j -0.06684418 17 -1.917560 + j 0.363855 1.920087 + j -0.475208 0.002527 + j -0.11135327 17 -0.128621 + j 0.132648 0.128754 + j -0.475531 0.00

32、0133 + j -0.34288422 21 6.093176 + j 1.070437 -6.064778 + j -0.847611 0.028398 + j 0.22282623 22 -0.406149 + j -1.116063 0.406824 + j 0.928040 0.000675 + j -0.18802424 23 -3.490793 + j -0.110399 3.516516 + j 0.138837 0.025723 + j 0.02843826 25 -0.771398 + j -0.442881 0.773189 + j -0.111580 0.001791

33、+ j -0.55446027 26 -2.681379 + j -0.887648 2.691475 + j 0.731900 0.010096 + j -0.15574828 26 1.416063 + j -0.565082 -1.408178 + j -0.210747 0.007885 + j -0.77583029 26 1.921038 + j -0.679443 -1.901899 + j -0.248272 0.019138 + j -0.92771529 28 3.491624 + j -0.395924 -3.476063 + j 0.289082 0.015561 +

34、j -0.106842 4 0 0.000000 + j -1.074485 0.000000 + j 0.000000 0.000000 + j -1.074485 5 0 0.000000 + j -2.182596 0.000000 + j 0.000000 0.000000 + j -2.18259611 12 0.018656 + j 0.473066 -0.018327 + j -0.464126 0.000329 + j 0.00894013 12 0.066943 + j 0.423225 -0.066673 + j -0.415874 0.000270 + j 0.00735

35、130 2 7.897633 + j -0.731582 -7.897633 + j 1.860277 0.000000 + j 1.12869531 6 7.506817 + j 1.371343 -7.506817 + j 0.109153 0.000000 + j 1.48049632 10 12.260592 + j 5.296517 -12.260592 + j -2.064007 0.000000 + j 3.23250934 20 5.080000 + j 1.826359 -5.054406 + j -1.314473 0.025594 + j 0.51188633 19 -1

36、.716763 + j 5.348910 1.736896 + j -4.940504 0.020133 + j 0.40840535 22 6.500000 + j 2.637566 -6.500000 + j -1.998477 0.000000 + j 0.63908936 23 1.402814 + j -0.195113 -1.401865 + j 0.246763 0.000949 + j 0.05165037 25 9.586236 + j 0.419689 -9.533808 + j 1.607517 0.052428 + j 2.02720638 29 -12.165903

37、+ j 2.106593 12.280860 + j 0.135062 0.114957 + j 2.24165520 19 -1.745594 + j 0.284473 1.747837 + j -0.240265 0.002242 + j 0.044208 结果分析：此程序的运行结果和试验程序给出的结果是一致的。说明程序无误，但在精确度上有微小差异，这主要是和导纳矩阵的精确度以及显示精度有关。心得：本程序分模块进行，先是排序，再是求导纳阵，然后求雅阁比，再进行迭代运算，程序本身很简洁明了，运行的时候只需要在matlab里输入main就行了，然后打开BUS和line所在的.m文件，结果就会自

38、动存在result文件中了，通过编写牛顿拉夫逊法matlab潮流计算程序复习了潮流计算的知识，也实现了计算机算法附录：实验源程序：Main函数：cleardfile,pathname=uigetfile(*.m,Select Data File);if pathname = 0 error( you must select a valid data file)else lfile =length(dfile); % strip off .m eval(dfile(1:lfile-2);end nb,mb=size(bus);nl,ml=size(line); % 计算bus和line矩阵的行数

39、和列数 bus,line,nPQ,nPV,nodenum = Num(bus,line); % 对节点重新排序的子程序Y = y(bus,line) % 计算节点导纳矩阵的子程序myf = fopen(Result.m,w);fprintf(myf,计算结果);fclose(myf); % 在当前目录下生成“Result.m”文件，写入节点导纳矩阵 format long EPS = 1.0e-10; % 设定误差精度for t = 1:100 % 开始迭代计算，设定最大迭代次数为100，以便不收敛情况下及时跳出dP,dQ = dPQ(Y,bus,nPQ,nPV); % 计算功率偏差dP和dQ

40、的子程序J = Jac(bus,Y,nPQ); % 计算雅克比矩阵的子程序UD = zeros(nPQ,nPQ);for i = 1:nPQUD(i,i) = bus(i,2); % 生成电压对角矩阵 endenddAngU = JdP;dQ;dAng = dAngU(1:nb-1,1); % 计算相角修正量dU = UD*(dAngU(nb:nb+nPQ-1,1); % 计算电压修正量bus(1:nPQ,2) = bus(1:nPQ,2) - dU; % 修正电压bus(1:nb-1,3) = bus(1:nb-1,3) - dAng; % 修正相角 if (max(abs(dU)EPS)&

41、(max(abs(dAng)EPS)break end % 判断是否满足精度误差，如满足则跳出，否则返回继续迭代endbus = PQ(bus,Y,nPQ,nPV); % 计算每个节点的有功和无功注入的子程序bus,line = ReNum(bus,line,nodenum); % 对节点恢复编号的子程序YtYm = YtYm(line); % 计算线路的等效Yt和Ym的子程序，以计算线路潮流 bus_res = bus_res(bus); % 计算节点数据结果的子程序S_res = S_res(bus,line,YtYm); % 计算线路潮流的子程序 myf = fopen(Result.m

42、,a);fprintf(myf,牛顿拉夫逊法潮流计算结果 节点计算结果：n节点 节点电压 节点相角（角度） 节点注入功率n);for i = 1:nbfprintf(myf,%2.0f ,bus_res(i,1);fprintf(myf,%10.6f ,bus_res(i,2); fprintf(myf,%10.6f ,bus_res(i,3); fprintf(myf,%10.6f + j %10.6fn,real(bus_res(i,4),imag(bus_res(i,4);endfprintf(myf,n线路计算结果：n节点I 节点J 线路功率S(I,J) 线路功率S(J,I) 线路损耗

43、dS(I,J)n);for i = 1:nlfprintf(myf,%2.0f ,S_res(i,1);fprintf(myf,%2.0f ,S_res(i,2);fprintf(myf,%10.6f + j %10.6f ,real(S_res(i,3),imag(S_res(i,3);fprintf(myf,%10.6f + j %10.6f ,real(S_res(i,4),imag(S_res(i,4);fprintf(myf,%10.6f + j%10.6fn,real(S_res(i,5),imag(S_res(i,5);endfclose(myf); % 迭代结束后继续在“Res

44、ult.m”写入节点计算结果和线路计算结果 程序结束Num.m 作用为对节点重排序，并修改相应的线路数据function bus,line,nPQ,nPV,nodenum = Num(bus,line)nb,mb=size(bus);nl,ml=size(line);nSW = 0; % number of swing bus counternPV = 0; % number of PV bus counternPQ = 0; % number of PQ bus counterfor i = 1:nb, % nb为总节点数type= bus(i,6);if type = 3,nSW = nS

45、W + 1; % increment swing bus counterSW(nSW,:)=bus(i,:);elseif type = 2,nPV = nPV +1; % increment PV bus counterPV(nPV,:)=bus(i,:);elsenPQ = nPQ + 1; % increment PQ bus counterPQ(nPQ,:)=bus(i,:);endend;bus=PQ;PV;SW;newbus=1:nb;nodenum=newbus bus(:,1);bus(:,1)=newbus;for i=1:nlfor j=1:2for k=1:nbif line(i,j)=nodenum(k,2)line(i,j)=nodenum(k,1);breakendendendendy.m 作用为计算节点导纳矩阵function Y = y(bus,line)nb,mb=size(bus);nl,ml=size(line);Y=zeros(nb,nb);for k=1:nlI=line(k,1);