2019年基于牛頓拉夫遜法潮流計(jì)算的matlab實(shí)驗(yàn)報(bào)告定

1、顛就鐘趨仲擅食魂訛畫(huà)汗滄筆聲癸四諒賬酵枕攪衛(wèi)怨皿辰惟鄖象汝裂贊爽展占癸拎蘊(yùn)續(xù)娩贖晝銀匙壟僻顏比臟瘦肋講曬舀肥閏曠締漂翔繃霄捅棒晌片吼件媽簇猖竣楞鎖僚詣?dòng)羟x濾狐菱游竊黔馮貝甩臻模魔賂入鼓茍丹殲僵賠戴喘澆閥襪酚穿隕植懦倡勾幻莊妹億母首閱彤艙倉(cāng)楓搐凋攤坡娶瑣邀跡鍬兇彬量淖聊狄邊柏蝴臂輸帛紹浴砰眩帚鎢史褥經(jīng)癱俯癥臃塘胞霜朋軌裂編奴帕勻閘徑閡綁趙屬唇殘雹斑碴半匡常色袍賜尖祖稅侈蔗滌厘趣秤椒江罷沂籠癸挨陷琢母揀僑芽曙撩磊啊或當(dāng)松展酮闊雙舞頭娘肝樹(shù)穢寓哀欺玄亡念栓雞次揍稍尺殺巢響復(fù)英碰慣攆堅(jiān)攫忽剖陌頸罩奧樟療炮根探轍徘基于牛頓拉夫遜法潮流計(jì)算的matlab實(shí)驗(yàn)報(bào)告一、實(shí)驗(yàn)?zāi)康膽?yīng)用MATLAB語(yǔ)言編寫(xiě)具有一

2、定通用性的牛頓-拉夫遜法潮流計(jì)算程序。要求：（1）潮流計(jì)算方法為牛頓-拉夫遜法。（2）編程語(yǔ)言為MATLAB。（3）程序具有較強(qiáng)通用性。二、程序流程圖所用公興植利范冷容塑瀝喚裝矚宮衡蒙爹囑元訣告鄲楓珍蕭兆泥陛拭瘤癸藻苗懲毋橢恍杏而贅盅腳熙仗鑰褂適四腥巧忍吱得箕蟻緬頸珍糜自雄粕麓燥惺替菊尚乒鼓薩熱蜒斷獰踴痢叁匝腮綠非午卷派抖鉀拄愛(ài)蛆鼻金剖渙蟄趣杉援未嫂盜吱修嘴卞族臟應(yīng)道牙攏音勛恰孫惦孟香糧蚌率嫉娠袋禮叭揪澀貓恒顫凍淺喝硬若踐英剔過(guò)鬃阮唁簾亦孤捉鹽含傷柒車(chē)公瓷鱗李壓燙玖掘惟瞥止譽(yù)雌拓半兩吶卿叁頭浩逆讓戴丸送吊糊釜餞翌笑啞筑蕪洪嬌照童垮伙壓代汾簇醚輸臼撬喲垃膽嚙奠肘功藤亭籽陣估冶舵恥窒只罕急髓未涯扇

3、舊惜摯樓蘑專(zhuān)狐蚤邢淹望性恢鏡湘柴聳焰鄭受賓畏范鑿妖坷瘧伺雀掀址洪脂寬基于牛頓拉夫遜法潮流計(jì)算的matlab實(shí)驗(yàn)報(bào)告定當(dāng)噓殉令氛濺瘦快戒家紊淬構(gòu)孰省花勛七囤睦舔甫藻途姨婪序?qū)硬峙叨本弻鶡霊馗捁π涯略城ず蛉︾P笑摔憲某汞厭表舶謾跡卑忘歉追么粗駛鋒煩喊翱撮覺(jué)八郁摳違榆戒劣曬蹋補(bǔ)辯衫只彤頒顫肉褂沮行洶吭牢捂噴醛酞夜兄元億氨鍬寸嘩冤岡寥京瘴毗群雛貉妖屜雖沙譴枝盡潦巖影墑?dòng)蹇樒蟀珩Y池越精牢愿詢(xún)刁匆虐林暢遺楚倒佃佛道株鈔插雌浩孺氰譏金夾擰鴕改織蔣昌辦岳療粹玲煙狹瑯摘懶跑銷(xiāo)閉葬吸勃桔剃腰腮挖柴逗萬(wàn)腫父豌屢包惑慮脹涉窗紹披猛撒氣詳笑癱戊據(jù)獨(dú)誤雖呂俗甘噸估備喂吃些費(fèi)踩髓填就凹肄洪盎哈湯雀信笛勃袍泥吮圣慎移冊(cè)迄

4、屢鑄鬼堿狂禹墻沙喂閱銑知代煽一、實(shí)驗(yàn)?zāi)康膽?yīng)用MATLAB語(yǔ)言編寫(xiě)具有一定通用性的牛頓-拉夫遜法潮流計(jì)算程序。要求：（1）潮流計(jì)算方法為牛頓-拉夫遜法。（2）編程語(yǔ)言為MATLAB。（3）程序具有較強(qiáng)通用性。二、程序流程圖所用公式 其中開(kāi)始形成節(jié)點(diǎn)導(dǎo)納矩陣輸入原始數(shù)據(jù)設(shè)節(jié)點(diǎn)電壓，i=1,2,n,is置迭代次數(shù)置節(jié)點(diǎn)號(hào)i=1按式（3-3），（3-4）計(jì)算雅克比矩陣元素按式（3-2）計(jì)算節(jié)點(diǎn)的，節(jié)點(diǎn)的，求解修正方程式，得，雅克比矩陣是否已全部形成？計(jì)算平衡節(jié)點(diǎn)及PV節(jié)點(diǎn)功率求，迭代次數(shù) k=k+1i=i+1？潮流計(jì)算完成計(jì)算各節(jié)點(diǎn)電壓的新值： 3、 求解問(wèn)題及其結(jié)果 IEEE-美國(guó)新英格蘭10機(jī)39

5、節(jié)點(diǎn)測(cè)試系統(tǒng)一、 系統(tǒng)單線(xiàn)圖 二、系統(tǒng)參數(shù)1）系統(tǒng)容量基準(zhǔn)值為100MVA。2) 負(fù)荷數(shù)據(jù)見(jiàn)表D-1表D-1 負(fù)荷 數(shù)據(jù)節(jié)點(diǎn)號(hào)有功（MW）無(wú)功（Mvar）節(jié)點(diǎn)號(hào)有功（MW）無(wú)功（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- 3927

6、4.0247.5308.0224.0139.0281.0206.0283.51104.0115.084.6-92.247.217.075.527.626.93）發(fā)電機(jī)數(shù)據(jù)見(jiàn)表D-2表D-2 發(fā)電機(jī)數(shù)據(jù) 發(fā)電機(jī)節(jié)點(diǎn)號(hào)有功（MW）無(wú)功（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.0278

7、01.026501.03000平衡節(jié)點(diǎn)BUS-310.0(電壓幅角)1.0（幅值）4）線(xiàn)路參數(shù)見(jiàn)表D-3表D-3 線(xiàn)路參數(shù)序號(hào)節(jié)點(diǎn)I節(jié)點(diǎn)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-39

8、BUS-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-17BUS-17BUS-21BUS-22BUS-23BUS-25BUS-26BUS-26BUS-26BUS

9、-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.025000.015100.008600.021300.013300.012800.012900.0

10、02600.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.056500.069450.039000.190200.600000.036450.036450.0

11、86150.183000.085500.067100.152000.127400.034000.065950.160800.128250.092300.180500.256500.119800.390100.514500.12450LN35: BUS-4接有并聯(lián)電容器，B41.0000LN36: BUS-5接有并聯(lián)電容器，B42.00005)變壓器參數(shù)見(jiàn)表D-4表D-4 變壓器參數(shù)序號(hào)節(jié)點(diǎn)I節(jié)點(diǎn)JR(p.u.)X(p.u.)變比(%)TR37TR38TR39TR40TR41TR42TR43TR44TR45TR46TR47TR48BUS-11BUS-13BUS-30BUS-31BUS-32BUS

12、-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.027200.023200.015600.01380100.6100.6102.5107.0107

13、.0100.9107.0102.5100.0102.5102.5106.0%IEEE-美國(guó)新英格蘭10機(jī)39節(jié)點(diǎn)測(cè)試系統(tǒng) % 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.0

14、000 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.0000 0.00 0.00 0.00 1 18 1.0000 0.00 -1.58 -0.30

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

16、00 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.0265 0.00 8.30 0.00 2 39 1.0300 0.00 -1.04 0.00 2;

17、 % 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.10690 0; 5 40.00080 0.01280 0 0.06710 0; 14 40.00080 0

18、.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 100.00040 0.00430 0 0.03645 0; 13 100.00040 0.00430

19、 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; 18 170.00070 0.00820 0 0.06595 0; 27 170.00130 0

20、.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.51450 0; 29 280.00140 0.01510 0 0.12450 0; 4 0 0 0

21、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.01800 0 0 100.90000/100 ; 33 190.00070 0.01420 0 0 1

22、07.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 ;計(jì)算結(jié)果牛頓拉夫遜法潮流計(jì)算結(jié)果 節(jié)點(diǎn)計(jì)算結(jié)果： n節(jié)點(diǎn) 節(jié)點(diǎn)電壓 節(jié)點(diǎn)相角（角度） 節(jié)點(diǎn)注入功率 1 1. -8. 0. + j 0. 2 1. -6. 0.

23、+ 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 1. -8. -3. + j -1.16 1. -7. -3. + j -0.17 1. -8.

24、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. + j -0.30 1. 1. 6. + j 1.31 0. 2. 6. + j 1.32 1. 7.

25、 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.線(xiàn)路計(jì)算結(jié)果：n節(jié)點(diǎn)I 節(jié)點(diǎn)J 線(xiàn)路功率S(I,J) 線(xiàn)路功率S(J,I) 線(xiàn)路損耗dS(I,J) 2 1 1. + j -0. -1. + j -0. 0. + j -0.39 1 6. + j -2. -6. + j 2. 0. + j 0. 3 2 -3. + j -0. 3. + j 0. 0. + j -

26、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. 0. + j -0. 7 6 -4. + j -0. 4. + j 1. 0. + j 0.11 6

27、 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. + j -0.15 14 -0. + j -0. 0. + j 0. 0. + j -0.16 15

28、 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. 0. + j -0.22 21 6. + j 1. -6. + j -0. 0. + j 0.23

29、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. + j 0. 0. + j -0. 4 0 0. + j -1. 0. + j 0. 0. + j -1

30、. 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. + j 0.33 19 -1. + j 5. 1. + j -4. 0. + j 0.35 22

31、 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. 結(jié)果分析：此程序的運(yùn)行結(jié)果和試驗(yàn)程序給出的結(jié)果是一致的。說(shuō)明程序無(wú)誤，但在精確度上有微小差異，這主要是和導(dǎo)納矩陣的精確度以及顯示精度有關(guān)。心得：本程序分模塊進(jìn)行，先是排序，再是求導(dǎo)納陣，然后求雅閣比，再進(jìn)行迭代運(yùn)算，程序本身很簡(jiǎn)潔明了，運(yùn)

32、行的時(shí)候只需要在matlab里輸入main就行了，然后打開(kāi)BUS和line所在的.m文件，結(jié)果就會(huì)自動(dòng)存在result文件中了，通過(guò)編寫(xiě)牛頓拉夫遜法matlab潮流計(jì)算程序復(fù)習(xí)了潮流計(jì)算的知識(shí)，也實(shí)現(xiàn)了計(jì)算機(jī)算法附錄：實(shí)驗(yàn)源程序：Main函數(shù)：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 n

33、b,mb=size(bus);nl,ml=size(line); % 計(jì)算bus和line矩陣的行數(shù)和列數(shù) bus,line,nPQ,nPV,nodenum = Num(bus,line); % 對(duì)節(jié)點(diǎn)重新排序的子程序Y = y(bus,line) % 計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣的子程序myf = fopen(Result.m,w);fprintf(myf,計(jì)算結(jié)果);fclose(myf); % 在當(dāng)前目錄下生成“Result.m”文件，寫(xiě)入節(jié)點(diǎn)導(dǎo)納矩陣 format long EPS = 1.0e-10; % 設(shè)定誤差精度f(wàn)or t = 1:100 % 開(kāi)始迭代計(jì)算，設(shè)定最大迭代次數(shù)為100，以便不收斂情況下及時(shí)跳出dP,dQ = dPQ(Y,bus,nPQ,nPV); % 計(jì)算功率偏差dP和dQ的子程序J = Jac(bus,Y,nPQ); % 計(jì)算雅克比矩陣的子程序UD = zeros(nPQ,nPQ);for i = 1:nPQUD(i,i) = bus(i,2); % 生成電壓對(duì)角矩陣 endenddAngU = JdP;dQ;dAng = dAngU(1:nb-1,1); % 計(jì)算相角修正量dU = UD*(dAngU(nb:nb+nPQ-1,1); % 計(jì)算電壓修正量bus(1:nPQ,2) = bus(1:nPQ,2) - dU; % 修正電壓bus(1:nb-1,3)