基于MATLAB的五节点牛拉潮流计算程序.docx

1、基于MATLAB的五节点牛拉潮流计算程序Smax请输入节点数:n=5请输入支路数:n1=5请输入平衡节点号:isb=5请输入误差精度:pr=0.00000001请输入支路参数:B1=1 2 0.04+0.25i 0.5i 1 0;1 3 0.1+0.35i 0 1 0;2 3 0.08+0.30i 0.5i 1 0;4 2 0.015i 0 1.05 1;5 3 0.03i 0 1.05 1请输入节点参数:B2=0 -0.8-0.4i 1.05 0 0 2;0 -2-1i 1.1 0 0 2;0 -3.7-3.5i 1.1 0 0 2;0 5+0i 1.05 1.05 0 3;0 0 1.05

2、 1.05 0 1节点号和对地参数:X=1 0;2 0;3 0;4 0;5 0Y = 1.3787 - 6.2917i -0.6240 + 3.9002i -0.7547 + 2.6415i 0 0 -0.6240 + 3.9002i 1.4539 -66.9808i -0.8299 + 3.1120i 0 +63.4921i 0 -0.7547 + 2.6415i -0.8299 + 3.1120i 1.5846 -35.7379i 0 0 +31.7460i 0 0 +63.4921i 0 0 -66.6667i 0 0 0 0 +31.7460i 0 0 -33.3333iOrgS =

3、-0.0724 -0.6191 0.0343 -0.5568 0.0415 -0.2404 0 0.1667PVU = 1.0500DetaS = -0.7276 0.2191 -2.0343 -0.4432 -3.7415 -3.2596 5.0000 0I = -0.0689 + 0.5896i 0.0312 + 0.5062i 0.0377 + 0.2185i 0 - 0.1587iJacbi = 7.1958 1.3787 -4.0952 -0.6552 -2.7736 -0.7925 0 0 -1.5166 6.0167 0.6552 -4.0952 0.7925 -2.7736 0

4、 0 -4.2902 -0.6864 74.1851 1.6305 -3.4232 -0.9129 -69.8413 0 0.6864 -4.2902 -1.5681 73.1727 0.9129 -3.4232 0 -69.8413 -2.9057 -0.8302 -3.4232 -0.9129 39.5302 1.7808 0 0 0.8302 -2.9057 0.9129 -3.4232 -1.7053 39.0931 0 0 0 0 -66.6667 0 0 0 69.8413 0 0 0 0 0 0 0 0 2.1000DetaU = 0.1182 -0.0229 0.4280 -0

5、.0034 -0.0450 -0.0998 0.4801 0B2 = 0 -0.8000 - 0.4000i 1.0271 + 0.1182i 0 0 2.0000 0 -2.0000 - 1.0000i 1.0966 + 0.4280i 0 0 2.0000 0 -3.7000 - 3.5000i 1.0002 - 0.0450i 0 0 2.0000 0 5.0000 1.0500 + 0.4801i 1.0500 0 3.0000 0 0 1.0500 1.0500 0 1.0000 OrgS = -0.8067 -0.4946 -1.8146 -1.8742 -3.4957 -3.04

6、99 4.8950 2.7173DetaS = 0.0067 0.0946 -0.1854 0.8742 -0.2043 -0.4501 0.1050 -0.2305I = -0.8298 + 0.3860i -2.0150 + 0.9228i -3.3513 + 3.2001i 4.8345 - 0.3775iJacbi = 7.0114 -0.1571 -4.0798 -0.1801 -2.8024 -0.4631 0 0 -1.5024 6.2393 0.1801 -4.0798 0.4631 -2.8024 0 0 -4.5438 0.9848 74.9932 -29.0852 -3.

7、7677 0.4218 -69.6225 27.1714 -0.9848 -4.5438 25.0552 73.1475 -0.4218 -3.7677 -27.1714 -69.6225 -2.6080 -0.8736 -3.0752 -0.9699 38.8722 -0.1595 0 0 0.8736 -2.6080 0.9699 -3.0752 -6.5431 32.4720 0 0 0 0 -66.6667 30.4818 0 0 69.6225 -27.1714 0 0 0 0 0 0 0.9602 2.1000DetaU = -0.0014 -0.0665 0.0125 -0.10

8、35 -0.0086 -0.0311 0.0135 -0.1159B2 = 0 -0.8000 - 0.4000i 0.9606 + 0.1168i 0 0 2.0000 0 -2.0000 - 1.0000i 0.9931 + 0.4405i 0 0 2.0000 0 -3.7000 - 3.5000i 0.9691 - 0.0535i 0 0 2.0000 0 5.0000 0.9341 + 0.4936i 1.0500 0 3.0000 0 0 1.0500 1.0500 0 1.0000 OrgS = -0.7973 -0.4046 -2.0020 -1.0797 -3.6971 -3

9、.4795 5.0032 1.7083DetaS = -0.0027 0.0046 0.0020 0.0797 -0.0029 -0.0205 -0.0032 -0.0136I = -0.8684 + 0.3156i -2.0875 + 0.1613i -3.6057 + 3.7897i 4.9425 + 0.7831iJacbi = 6.5206 -0.2786 -3.8195 -0.1440 -2.6257 -0.4165 0 0 -1.4581 5.8893 0.1440 -3.8195 0.4165 -2.6257 0 0 -4.1481 1.0981 67.3200 -30.1461

10、 -3.4561 0.5466 -63.0536 27.9657 -1.0981 -4.1481 25.9710 66.9973 -0.5466 -3.4561 -27.9657 -63.0536 -2.5195 -0.8728 -2.9714 -0.9708 38.3381 -0.1571 0 0 0.8728 -2.5195 0.9708 -2.9714 -7.0544 30.7588 0 0 0 0 -59.3052 31.3412 0 0 63.0536 -27.9657 0 0 0 0 0 0 0.9872 1.8681DetaU = -0.0001 -0.0044 0.0014 -

11、0.0067 -0.0003 -0.0018 0.0011 -0.0079B2 = 0 -0.8000 - 0.4000i 0.9563 + 0.1166i 0 0 2.0000 0 -2.0000 - 1.0000i 0.9864 + 0.4418i 0 0 2.0000 0 -3.7000 - 3.5000i 0.9673 - 0.0538i 0 0 2.0000 0 5.0000 0.9262 + 0.4947i 1.0500 0 3.0000 0 0 1.0500 1.0500 0 1.0000 OrgS = -0.8000 -0.4000 -2.0002 -1.0005 -3.700

12、0 -3.4999 5.0002 1.6182DetaS = 1.0e-003 * -0.0277 0.0166 0.2361 0.4592 -0.0048 -0.0536 -0.2356 -0.0630I = -0.8746 + 0.3116i -2.0673 + 0.0882i -3.6126 + 3.8190i 4.9264 + 0.8842iJacbi = 6.4889 -0.2900 -3.8023 -0.1418 -2.6140 -0.4136 0 0 -1.4592 5.8656 0.1418 -3.8023 0.4136 -2.6140 0 0 -4.1229 1.1077 6

13、6.8024 -30.2283 -3.4365 0.5564 -62.6304 28.0537 -1.1077 -4.1229 26.0937 66.6260 -0.5564 -3.4365 -28.0537 -62.6304 -2.5146 -0.8721 -2.9657 -0.9701 38.3040 -0.1578 0 0 0.8721 -2.5146 0.9701 -2.9657 -7.0675 30.6659 0 0 0 0 -58.8059 31.4096 0 0 62.6304 -28.0537 0 0 0 0 0 0 0.9894 1.8524DetaU = 1.0e-004

14、* -0.0370 -0.1960 0.0504 -0.2868 -0.0118 -0.0645 0.0010 -0.3406B2 = 0 -0.8000 - 0.4000i 0.9562 + 0.1166i 0 0 2.0000 0 -2.0000 - 1.0000i 0.9864 + 0.4419i 0 0 2.0000 0 -3.7000 - 3.5000i 0.9673 - 0.0538i 0 0 2.0000 0 5.0000 0.9262 + 0.4947i 1.0500 0 3.0000 0 0 1.0500 1.0500 0 1.0000 OrgS = -0.8000 -0.4

15、000 -2.0000 -1.0000 -3.7000 -3.5000 5.0000 1.6176DetaS = 1.0e-007 * -0.0091 0.0009 0.1099 0.0776 -0.0002 -0.0058 -0.1073 -0.0116I = -0.8746 + 0.3116i -2.0670 + 0.0879i -3.6127 + 3.8191i 4.9261 + 0.8847iJacbi = 6.4888 -0.2900 -3.8023 -0.1418 -2.6139 -0.4136 0 0 -1.4592 5.8655 0.1418 -3.8023 0.4136 -2

16、.6139 0 0 -4.1228 1.1077 66.8002 -30.2284 -3.4364 0.5565 -62.6286 28.0540 -1.1077 -4.1228 26.0944 66.6244 -0.5565 -3.4364 -28.0540 -62.6286 -2.5146 -0.8721 -2.9657 -0.9701 38.3038 -0.1578 0 0 0.8721 -2.5146 0.9701 -2.9657 -7.0675 30.6656 0 0 0 0 -58.8037 31.4096 0 0 62.6286 -28.0540 0 0 0 0 0 0 0.98

17、94 1.8523DetaU = 1.0e-009 * -0.1293 -0.4166 0.1004 -0.5398 -0.0250 -0.1106 -0.0701 -0.5889B2 = 0 -0.8000 - 0.4000i 0.9562 + 0.1166i 0 0 2.0000 0 -2.0000 - 1.0000i 0.9864 + 0.4419i 0 0 2.0000 0 -3.7000 - 3.5000i 0.9673 - 0.0538i 0 0 2.0000 0 5.0000 0.9262 + 0.4947i 1.0500 0 3.0000 0 0 1.0500 1.0500 0

18、 1.0000 Times = 5请输入节点数:n=5请输入支路数:n1=5请输入平衡节点号:isb=5请输入误差精度:pr=0.00000001请输入支路参数:B1=1 2 0.04+0.25i 0.5i 1 0;1 3 0.1+0.35i 0 1 0;2 3 0.08+0.30i 0.5i 1 0;4 2 0.015i 0 1.05 1;5 3 0.03i 0 1.05 1请输入节点参数:B2=0 -0.8-0.4i 1.05 0 0 2;0 -2-1i 1.1 0 0 2;0 -3.7+0.9i 1.1 0 0 2;0 5+0i 1.05 1.05 0 3;0 0 1.05 1.05 0

19、 1节点号和对地参数:X=1 0;2 0;3 0;4 0;5 0Y = 1.3787 - 6.2917i -0.6240 + 3.9002i -0.7547 + 2.6415i 0 0 -0.6240 + 3.9002i 1.4539 -66.9808i -0.8299 + 3.1120i 0 +63.4921i 0 -0.7547 + 2.6415i -0.8299 + 3.1120i 1.5846 -35.7379i 0 0 +31.7460i 0 0 +63.4921i 0 0 -66.6667i 0 0 0 0 +31.7460i 0 0 -33.3333iOrgS = -0.0724

20、 -0.6191 0.0343 -0.5568 0.0415 -0.2404 0 0.1667PVU = 1.0500DetaS = -0.7276 0.2191 -2.0343 -0.4432 -3.7415 1.1404 5.0000 0I = -0.0689 + 0.5896i 0.0312 + 0.5062i 0.0377 + 0.2185i 0 - 0.1587iJacbi = 7.1958 1.3787 -4.0952 -0.6552 -2.7736 -0.7925 0 0 -1.5166 6.0167 0.6552 -4.0952 0.7925 -2.7736 0 0 -4.29

21、02 -0.6864 74.1851 1.6305 -3.4232 -0.9129 -69.8413 0 0.6864 -4.2902 -1.5681 73.1727 0.9129 -3.4232 0 -69.8413 -2.9057 -0.8302 -3.4232 -0.9129 39.5302 1.7808 0 0 0.8302 -2.9057 0.9129 -3.4232 -1.7053 39.0931 0 0 0 0 -66.6667 0 0 0 69.8413 0 0 0 0 0 0 0 0 2.1000DetaU = 0.1359 0.0390 0.4559 0.0061 -0.0

22、450 0.0171 0.5067 0B2 = 0 -0.8000 - 0.4000i 1.0890 + 0.1359i 0 0 2.0000 0 -2.0000 - 1.0000i 1.1061 + 0.4559i 0 0 2.0000 0 -3.7000 + 0.9000i 1.1171 - 0.0450i 0 0 2.0000 0 5.0000 1.0500 + 0.5067i 1.0500 0 3.0000 0 0 1.0500 1.0500 0 1.0000 OrgS = -0.8624 -0.4949 -1.8242 -1.9413 -3.7112 1.0705 5.1957 2.

23、2129DetaS = 0.0624 0.0949 -0.1758 0.9413 0.0112 -0.1705 -0.1957 -0.2568I = -0.8357 + 0.3501i -2.0281 + 0.9192i -3.3553 - 0.8231i 4.8385 + 0.2275iJacbi = 7.3891 -0.1895 -4.3320 -0.1494 -2.9792 -0.4628 0 0 -1.4818 6.6888 0.1494 -4.3320 0.4628 -2.9792 0 0 -4.5984 1.0877 75.6684 -30.9537 -3.8205 0.5007

24、-70.2275 28.9433 -1.0877 -4.5984 26.8975 73.8300 -0.5007 -3.8205 -28.9433 -70.2275 -2.9169 -0.9619 -3.4391 -1.0671 39.0286 0.0228 0 0 0.9619 -2.9169 1.0671 -3.4391 -6.7334 40.6748 0 0 0 0 -66.6667 32.1732 0 0 70.2275 -28.9433 0 0 0 0 0 0 1.0135 2.1000DetaU = -0.0192 -0.0564 -0.0341 -0.0934 -0.0081 -

25、0.0161 -0.0357 -0.1051B2 = 0 -0.8000 - 0.4000i 1.0326 + 0.1168i 0 0 2.0000 0 -2.0000 - 1.0000i 1.0127 + 0.4217i 0 0 2.0000 0 -3.7000 + 0.9000i 1.1010 - 0.0531i 0 0 2.0000 0 5.0000 0.9449 + 0.4711i 1.0500 0 3.0000 0 0 1.0500 1.0500 0 1.0000 OrgS = -0.8005 -0.4034 -1.9749 -1.0664 -3.7007 0.9030 4.9840

26、 0.9512DetaS = 0.0005 0.0034 -0.0251 0.0664 0.0007 -0.0030 0.0160 -0.0123I = -0.8091 + 0.2992i -2.0357 + 0.2053i -3.3930 - 0.6566i 4.6265 + 1.2997iJacbi = 6.9567 -0.1202 -4.1000 -0.1889 -2.8157 -0.4708 0 0 -1.4980 6.3584 0.1889 -4.1000 0.4708 -2.8157 0 0 -4.2127 1.0129 68.6473 -28.8118 -3.5014 0.472

27、1 -64.2960 26.7771 -1.0129 -4.2127 24.7404 68.2368 -0.4721 -3.5014 -26.7771 -64.2960 -2.8682 -0.9711 -3.3822 -1.0789 38.6059 0.2486 0 0 0.9711 -2.8682 1.0789 -3.3822 -7.0345 39.9191 0 0 0 0 -59.9964 29.9082 0 0 64.2960 -26.7771 0 0 0 0 0 0 0.9421 1.8899DetaU = -0.0004 -0.0033 -0.0003 -0.0053 -0.0003

28、 -0.0008 -0.0002 -0.0064B2 = 0 -0.8000 - 0.4000i 1.0292 + 0.1163i 0 0 2.0000 0 -2.0000 - 1.0000i 1.0074 + 0.4215i 0 0 2.0000 0 -3.7000 + 0.9000i 1.1002 - 0.0533i 0 0 2.0000 0 5.0000 0.9385 + 0.4708i 1.0500 0 3.0000 0 0 1.0500 1.0500 0 1.0000 OrgS = -0.8000 -0.4000 -1.9999 -1.0004 -3.7000 0.9000 5.00

29、00 0.8751DetaS = 1.0e-003 * -0.0067 0.0064 -0.0577 0.3536 -0.0006 -0.0030 0.0407 -0.0411I = -0.8108 + 0.2970i -2.0432 + 0.1382i -3.3947 - 0.6534i 4.6299 + 1.3904iJacbi = 6.9331 -0.1237 -4.0868 -0.1886 -2.8066 -0.4695 0 0 -1.4980 6.3391 0.1886 -4.0868 0.4695 -2.8066 0 0 -4.1919 1.0152 68.2247 -28.808

30、7 -3.4847 0.4756 -63.9593 26.7597 -1.0152 -4.1919 24.7224 67.9483 -0.4756 -3.4847 -26.7597 -63.9593 -2.8659 -0.9712 -3.3796 -1.0790 38.5804 0.2550 0 0 0.9712 -2.8659 1.0790 -3.3796 -7.0445 39.8873 0 0 0 0 -59.5895 29.8949 0 0 63.9593 -26.7597 0 0 0 0 0 0 0.9417 1.8771DetaU = 1.0e-004 * -0.0201 -0.1093 0.0017 -0.1626 -0.0086 -0.0234 -0.006