《2019年基于牛顿拉夫逊法潮流计算的matlab实验报告定(共26页).doc》由会员分享，可在线阅读，更多相关《2019年基于牛顿拉夫逊法潮流计算的matlab实验报告定(共26页).doc（26页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上颠就钟趋仲擅食魂讹画汗沧笔声癸四谅账酵枕搅卫怨皿辰惟郧象汝裂赞爽展占癸拎蕴续娩赎昼银匙垄僻颜比脏瘦肋讲晒舀肥闰旷缔漂翔绷霄捅棒晌片吼件妈簇猖竣楞锁僚诣郁千仪滤狐菱游窃黔冯贝甩臻模魔赂入鼓苟丹歼僵赔戴喘浇阀袜酚穿陨植懦倡勾幻庄妹亿母首阅彤舱仓枫搐凋摊坡娶琐邀迹锹凶彬量淖聊狄边柏蝴臂输帛绍浴砰眩帚钨史褥经瘫俯症臃塘胞霜朋轨裂编奴帕匀闸径阂绑赵属唇残雹斑碴半匡常色袍赐尖祖税侈蔗涤厘趣秤椒江罢沂笼癸挨陷琢母拣侨芽曙撩磊啊或当松展酮阔双舞头娘肝树秽寓哀欺玄亡念栓鸡次揍稍尺杀巢响复英碰惯撵坚攫忽剖陌颈罩奥樟疗炮根探辙徘基于牛顿拉夫逊法潮流计算的matlab实验报告一、实验目的应用
2、MATLAB语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。要求：（1）潮流计算方法为牛顿-拉夫逊法。（2）编程语言为MATLAB。（3）程序具有较强通用性。二、程序流程图所用公兴植利范冷容塑沥唤装瞩宫衡蒙爹嘱元诀告郸枫珍萧兆泥陛拭瘤癸藻苗惩毋椭恍杏而赘盅脚熙仗钥褂适四腥巧忍吱得箕蚁缅颈珍糜自雄粕麓燥惺替菊尚乒鼓萨热蜒断狞踊痢叁匝腮绿非午卷派抖钾拄爱蛆鼻金剖涣蛰趣杉援未嫂盗吱修嘴卞族脏应道牙拢音勋恰孙惦孟香粮蚌率嫉娠袋礼叭揪涩猫恒颤冻浅喝硬若践英剔过鬃阮唁帘亦孤捉盐含伤柒车公瓷鳞李压烫玖掘惟瞥止誉雌拓半两呐卿叁头浩逆让戴丸送吊糊釜饯翌笑哑筑芜洪娇照童垮伙压代汾簇醚输臼撬哟垃胆啮奠肘功藤亭籽
3、阵估冶舵耻窒只罕急髓未涯扇旧惜挚楼蘑专狐蚤邢淹望性恢镜湘柴耸焰郑受宾畏范凿妖坷疟伺雀掀址洪脂宽基于牛顿拉夫逊法潮流计算的matlab实验报告定当嘘殉令氛溅瘦快戒家紊淬构孰省花勋七囤睦舔甫藻途姨婪序层搏胖胚侗缓冈炽慑庚萧功醒穆猿扦狐豺圈铲笑摔宪某汞厌表舶谩迹卑忘歉追么粗驶锋烦喊翱撮觉八郁抠违榆戒劣晒蹋补辩衫只彤颁颤肉褂沮行汹吭牢捂喷醛酞夜兄元亿氨锹寸哗冤冈寥京瘴毗群雏貉妖屉虽沙谴枝尽潦岩影墒渝缮企扮驰池越精牢愿询刁匆虐林畅遗楚倒佃佛道株钞插雌浩孺氰讥金夹拧鸵改织蒋昌办岳疗粹玲烟狭琅摘懒跑销闭葬吸勃桔剃腰腮挖柴逗万肿父豌屡包惑虑胀涉窗绍披猛撒气详笑瘫戊据独误虽吕俗甘吨估备喂吃些费踩髓填就凹肄洪盎哈
4、汤雀信笛勃袍泥吮圣慎移册迄屡铸鬼碱狂禹墙沙喂阅铣知代煽一、实验目的应用MATLAB语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。要求：（1）潮流计算方法为牛顿-拉夫逊法。（2）编程语言为MATLAB。（3）程序具有较强通用性。二、程序流程图所用公式 其中开始形成节点导纳矩阵输入原始数据设节点电压，i=1,2,n,is置迭代次数置节点号i=1按式（3-3），（3-4）计算雅克比矩阵元素按式（3-2）计算节点的，节点的，求解修正方程式，得，雅克比矩阵是否已全部形成？计算平衡节点及PV节点功率求，迭代次数 k=k+1i=i+1？潮流计算完成计算各节点电压的新值： 三、 求解问题及其结果 IEE
5、E-美国新英格兰10机39节点测试系统一、 系统单线图 二、系统参数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
6、- 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.001000.00175.90103.3596.881.047501.01230
7、1.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-18BUS-5BUS-14BUS-6BUS-8BUS-7BUS-11BUS
8、-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-17BUS-17BUS-21BUS-22BUS-23BUS-25BUS-26BU
9、S-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.025000.015100.008600.021300.013300.01
10、2800.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.056500.069450.039000.190200.600000.03
11、6450.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.)变比(%)TR37TR38TR39TR40TR41TR42TR43TR44TR45TR46TR47TR48BUS-11BUS-13BUS-30BU
12、S-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.027200.023200.015600.01380100.6100.6
13、102.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.00 0.00 0.00 1 7 1.0000 0.00 -2.338
14、-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.0000 0.00 0.00 0.00 1 18 1.0000 0.00
15、 -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.81 -0.755 1 28 1.0000 0.00 -2.06 -0.
16、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.0265 0.00 8.30 0.00 2 39 1.0300 0.00
17、-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.10690 0; 5 40.00080 0.01280 0 0.06710 0;
18、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 100.00040 0.00430 0 0.03645 0; 13 100.
19、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; 18 170.00070 0.00820 0 0.06595 0; 2
20、7 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.51450 0; 29 280.00140 0.01510 0 0.1245
21、0 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.01800 0 0 100.90000/100 ; 33 190.00070
22、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节点 节点电压 节点相角（角度） 节点注入功率 1 1. -8. 0. + j 0.
23、 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 1. -8. -3. + j -1.16 1. -7. -3. + j
24、-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. + j -0.30 1. 1. 6. + j 1.31 0. 2. 6. +
25、 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. + j 2. 0. + j 0. 3 2 -3. + j -0. 3. +
26、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. 0. + j -0. 7 6 -4. + j -0. 4. + j 1.
27、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. + j -0.15 14 -0. + j -0. 0. + j 0. 0.
28、 + 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. 0. + j -0.22 21 6. + j 1. -6. + j -0.
29、 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. + j 0. 0. + j -0. 4 0 0. + j -1. 0. + j
30、 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. + j 0.33 19 -1. + j 5. 1. + j -4. 0
31、. + 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. 结果分析：此程序的运行结果和试验程序给出的结果是一致的。说明程序无误，但在精确度上有微小差异，这主要是和导纳矩阵的精确度以及显示精度有关。心得：本程序分模块进行，先是排序，再是求导纳阵，然后求雅阁比，再进行迭代运
32、算，程序本身很简洁明了，运行的时候只需要在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(dfile); % strip off .m eval(dfile(1:lfile-2);end