基于牛顿拉夫逊法潮流计算的matlab实验报告(共24页).doc

《基于牛顿拉夫逊法潮流计算的matlab实验报告(共24页).doc》由会员分享，可在线阅读，更多相关《基于牛顿拉夫逊法潮流计算的matlab实验报告(共24页).doc（24页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、精选优质文档-倾情为你奉上基于牛顿拉夫逊法潮流计算的matlab实验报告一、 实验目的和要求1.学习掌握matlab的基本用法2.应用MATLAB语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。要求：（1）潮流计算方法为牛顿-拉夫逊法。（2）编程语言为MATLAB。（3）程序具有较强通用性。二、程序流程图1.程序流程图开始形成节点导纳矩阵输入原始数据，节点重新编号设节点电压初值，i=1,2,n,is置迭代次数P=0置节点号i=1计算雅克比矩阵元素按公式计算节点的，节点的，求解修正方程式，得，雅克比矩阵是否已全部形成？节点编号逆回，算平衡节点及PV节点功率求，迭代次数P=P+1i=i+1？潮

2、流计算完成计算各节点电压的新值： 三、 求解问题及其结果1. 求解问题：IEEE-美国新英格兰10机39节点测试系统1）系统单线图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- 23B

3、US- 24BUS- 25BUS- 26BUS- 27BUS- 28BUS- 29BUS- 39274.0247.5308.0224.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.001

4、000.00175.90103.3596.881.047501.012301.049301.027801.026501.03000平衡节点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-

5、18BUS-5BUS-14BUS-6BUS-8BUS-7BUS-11BUS-8BUS-9BUS-39BUS-11BUS-13BUS-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-17B

6、US-17BUS-21BUS-22BUS-23BUS-25BUS-26BUS-26BUS-26BUS-280.003500.001000.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.

7、025000.015100.008600.021300.013300.012800.012900.002600.011200.009200.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.

8、056500.069450.039000.190200.600000.036450.036450.086150.183000.085500.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.)变比(%)TR37TR38TR39TR40TR41TR42TR43TR4

9、4TR45TR46TR47TR48BUS-11BUS-13BUS-30BUS-31BUS-32BUS-34BUS-33BUS-35BUS-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

10、.027200.023200.015600.01380100.6100.6102.5107.0107.0100.9107.0102.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

11、.00 0.00 0.00 1 7 1.0000 0.00 -2.338 -0.84 1 8 1.0000 0.00 -5.22 -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

12、.0000 0.00 0.00 0.00 1 18 1.0000 0.00 -1.58 -0.30 1 19 1.0000 0.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

13、.81 -0.755 1 28 1.0000 0.00 -2.06 -0.276 1 29 1.0000 0.00 -2.835 -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.

14、0265 0.00 8.30 0.00 2 39 1.0300 0.00 -1.04 0.00 2; % 1 2 3 4 5 6 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.1069

15、0 0; 5 40.00080 0.01280 0 0.06710 0; 14 40.00080 0.01290 0 0.06910 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 1

16、00.00040 0.00430 0 0.03645 0; 13 100.00040 0.00430 0 0.03645 0; 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

17、; 18 170.00070 0.00820 0 0.06595 0; 27 170.00130 0.01730 0 0.16080 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.5

18、1450 0; 29 280.00140 0.01510 0 0.12450 0; 4 0 0 0 0 1.0000 0; 5 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.01

19、800 0 0 100.90000/100 ; 33 190.00070 0.01420 0 0 107.00000/100 ; 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节点 节点

20、电压 节点相角（角度） 节点注入功率 1 1. -8. 0. + j 0. 2 1. -6. 0. + j 0. 3 1. -9. -3. + j -0. 4 1. -10. -5. + j -1. 5 1. -8. 0. + j 0. 6 1. -8. 0. + j 0. 7 1. -10. -2. + j -0. 8 1. -10. -5. + j -1. 9 1. -10. 0. + j 0.10 1. -6. 0. + j 0.11 1. -6. 0. + j 0.12 1. -6. -0. + j -0.13 1. -6. 0. + j 0.14 1. -8. 0. + j 0.15

21、 1. -8. -3. + j -1.16 1. -7. -3. + j -0.17 1. -8. 0. + j 0.18 1. -8. -1. + j -0.19 1. -2. 0. + j 0.20 0. -3. -6. + j -1.21 1. -4. -2. + j -1.22 1. -0. 0. + j 0.23 1. -0. -2. + j -1.24 1. -6. -3. + j -0.25 1. -4. -2. + j -0.26 1. -6. -1. + j -0.27 1. -8. -2. + j -0.28 1. -2. -2. + j -0.29 1. 0. -2. +

22、 j -0.30 1. 1. 6. + j 1.31 0. 2. 6. + j 1.32 1. 7. 5. + j 0.33 1. -3. 2. + j 1.34 1. 1. 5. + j 1.35 1. 4. 6. + j 2.36 1. 1. 5. + j -0.37 1. 7. 8. + j 0.38 1. -10. -1. + j -2.39 1. 0. 5. + j 1.线路计算结果：n节点I 节点J 线路功率S(I,J) 线路功率S(J,I) 线路损耗dS(I,J) 2 1 1. + j -0. -1. + j -0. 0. + j -0.39 1 6. + j -2. -6. +

23、 j 2. 0. + j 0. 3 2 -3. + j -0. 3. + j 0. 0. + j -0.25 2 2. + j -1. -2. + j 0. 0. + j -0. 4 3 -0. + j -0. 0. + j 0. 0. + j -0.18 3 0. + j -0. -0. + j 0. 0. + j -0. 5 4 1. + j 0. -1. + j -0. 0. + j -0.14 4 2. + j -0. -2. + j 0. 0. + j -0. 6 5 4. + j -0. -4. + j 0. 0. + j 0. 8 5 -3. + j -1. 3. + j 0.

24、0. + j -0. 7 6 -4. + j -0. 4. + j 1. 0. + j 0.11 6 3. + j -0. -3. + j 0. 0. + j -0. 8 7 -1. + j -0. 1. + j 0. 0. + j -0. 9 8 0. + j 0. -0. + j -0. 0. + j -0.39 9 7. + j -1. -7. + j 2. 0. + j 0.11 10 -3. + j -0. 3. + j 0. 0. + j -0.13 10 -3. + j -0. 3. + j 0. 0. + j -0.14 13 -2. + j -0. 2. + j 0. 0.

25、+ j -0.15 14 -0. + j -0. 0. + j 0. 0. + j -0.16 15 2. + j 0. -2. + j -0. 0. + j -0.17 16 -2. + j 0. 2. + j -1. 0. + j -0.19 16 4. + j 0. -4. + j -0. 0. + j 0.21 16 3. + j -0. -3. + j 0. 0. + j -0.24 16 0. + j -0. -0. + j 0. 0. + j -0.18 17 -1. + j 0. 1. + j -0. 0. + j -0.27 17 -0. + j 0. 0. + j -0.

26、0. + j -0.22 21 6. + j 1. -6. + j -0. 0. + j 0.23 22 -0. + j -1. 0. + j 0. 0. + j -0.24 23 -3. + j -0. 3. + j 0. 0. + j 0.26 25 -0. + j -0. 0. + j -0. 0. + j -0.27 26 -2. + j -0. 2. + j 0. 0. + j -0.28 26 1. + j -0. -1. + j -0. 0. + j -0.29 26 1. + j -0. -1. + j -0. 0. + j -0.29 28 3. + j -0. -3. +

27、j 0. 0. + j -0. 4 0 0. + j -1. 0. + j 0. 0. + j -1. 5 0 0. + j -2. 0. + j 0. 0. + j -2.11 12 0. + j 0. -0. + j -0. 0. + j 0.13 12 0. + j 0. -0. + j -0. 0. + j 0.30 2 7. + j -0. -7. + j 1. 0. + j 1.31 6 7. + j 1. -7. + j 0. 0. + j 1.32 10 12. + j 5. -12. + j -2. 0. + j 3.34 20 5. + j 1. -5. + j -1. 0

28、. + j 0.33 19 -1. + j 5. 1. + j -4. 0. + j 0.35 22 6. + j 2. -6. + j -1. 0. + j 0.36 23 1. + j -0. -1. + j 0. 0. + j 0.37 25 9. + j 0. -9. + j 1. 0. + j 2.38 29 -12. + j 2. 12. + j 0. 0. + j 2.20 19 -1. + j 0. 1. + j -0. 0. + j 0. 结果分析：此程序的运行结果和试验程序给出的结果是一致的。说明程序无误，但在精确度上有微小差异，这主要是和导纳矩阵的精确度以及显示精度有关。

29、心得：本程序分模块进行，先是排序，再是求导纳阵，然后求雅阁比，再进行迭代运算，程序本身很简洁明了，运行的时候只需要在matlab里输入main就行了，然后打开BUS和line所在的.m文件，结果就会自动存在result文件中了，通过编写牛顿拉夫逊法matlab潮流计算程序复习了潮流计算的知识，也实现了计算机算法附录：实验源程序：Main函数：cleardfile,pathname=uigetfile(*.m,Select Data File);if pathname = 0 error( you must select a valid data file)else lfile =length(

30、dfile); % strip off .m eval(dfile(1:lfile-2);end nb,mb=size(bus);nl,ml=size(line); % 计算bus和line矩阵的行数和列数 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-1

31、0; % 设定误差精度for t = 1:100 % 开始迭代计算，设定最大迭代次数为100，以便不收敛情况下及时跳出dP,dQ = dPQ(Y,bus,nPQ,nPV); % 计算功率偏差dP和dQ的子程序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; %