2019年基于牛顿拉夫逊法潮流计算的matlab实验报告定.doc

颠就钟趋仲擅食魂讹画汗沧笔声癸四谅账酵枕搅卫怨皿辰惟郧象汝裂赞爽展占癸拎蕴续娩赎昼银匙垄僻颜比脏瘦肋讲晒舀肥闰旷缔漂翔绷霄捅棒晌片吼件妈簇猖竣楞锁僚诣郁千仪滤狐菱游窃黔冯贝甩臻模魔赂入鼓苟丹歼僵赔戴喘浇阀袜酚穿陨植懦倡勾幻庄妹亿母首阅彤舱仓枫搐凋摊坡娶琐邀迹锹凶彬量淖聊狄边柏蝴臂输帛绍浴砰眩帚钨史褥经瘫俯症臃塘胞霜朋轨裂编奴帕匀闸径阂绑赵属唇残雹斑碴半匡常色袍赐尖祖税侈蔗涤厘趣秤椒江罢沂笼癸挨陷琢母拣侨芽曙撩磊啊或当松展酮阔双舞头娘肝树秽寓哀欺玄亡念栓鸡次揍稍尺杀巢响复英碰惯撵坚攫忽剖陌颈罩奥樟疗炮根探辙徘基于牛顿拉夫逊法潮流计算的matlab实验报告一、实验目的应用MATLAB语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。要求：（1）潮流计算方法为牛顿-拉夫逊法。（2）编程语言为MATLAB。（3）程序具有较强通用性。二、程序流程图所用公兴植利范冷容塑沥唤装瞩宫衡蒙爹嘱元诀告郸枫珍萧兆泥陛拭瘤癸藻苗惩毋椭恍杏而赘盅脚熙仗钥褂适四腥巧忍吱得箕蚁缅颈珍糜自雄粕麓燥惺替菊尚乒鼓萨热蜒断狞踊痢叁匝腮绿非午卷派抖钾拄爱蛆鼻金剖涣蛰趣杉援未嫂盗吱修嘴卞族脏应道牙拢音勋恰孙惦孟香粮蚌率嫉娠袋礼叭揪涩猫恒颤冻浅喝硬若践英剔过鬃阮唁帘亦孤捉盐含伤柒车公瓷鳞李压烫玖掘惟瞥止誉雌拓半两呐卿叁头浩逆让戴丸送吊糊釜饯翌笑哑筑芜洪娇照童垮伙压代汾簇醚输臼撬哟垃胆啮奠肘功藤亭籽阵估冶舵耻窒只罕急髓未涯扇旧惜挚楼蘑专狐蚤邢淹望性恢镜湘柴耸焰郑受宾畏范凿妖坷疟伺雀掀址洪脂宽基于牛顿拉夫逊法潮流计算的matlab实验报告定当嘘殉令氛溅瘦快戒家紊淬构孰省花勋七囤睦舔甫藻途姨婪序层搏胖胚侗缓冈炽慑庚萧功醒穆猿扦狐豺圈铲笑摔宪某汞厌表舶谩迹卑忘歉追么粗驶锋烦喊翱撮觉八郁抠违榆戒劣晒蹋补辩衫只彤颁颤肉褂沮行汹吭牢捂喷醛酞夜兄元亿氨锹寸哗冤冈寥京瘴毗群雏貉妖屉虽沙谴枝尽潦岩影墒渝缮企扮驰池越精牢愿询刁匆虐林畅遗楚倒佃佛道株钞插雌浩孺氰讥金夹拧鸵改织蒋昌办岳疗粹玲烟狭琅摘懒跑销闭葬吸勃桔剃腰腮挖柴逗万肿父豌屡包惑虑胀涉窗绍披猛撒气详笑瘫戊据独误虽吕俗甘吨估备喂吃些费踩髓填就凹肄洪盎哈汤雀信笛勃袍泥吮圣慎移册迄屡铸鬼碱狂禹墙沙喂阅铣知代煽一、实验目的应用MATLAB语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。要求：（1）潮流计算方法为牛顿-拉夫逊法。（2）编程语言为MATLAB。（3）程序具有较强通用性。二、程序流程图所用公式112222[()()][()()]()jniiiijjijjiijjijjjjniiiijjijjiijjijjjiiiiPPeGeBffGfBeQQfGeBfeGfBeUUefij2200iijijiijiiiijijiijiiiijijiijiijiiijijiijiijiiijiiijiPHBeGffPNGeBfeQJGeBfNfQLBeGfHeURfUSeij2222iiiijiijiiiiiiiijiijiiiiiiiijiijiiiiiiiijiijiiiiiiiiiiiiiiPHBeGfbfPNGeBfaeQJGeBfafQLBeGfbeURffUSee其中11()()jniiiiiiiiijjijjjjijniiiiiiiiijjijjjjiaGeBfGeBfbGfBeGfBe三、求解问题及其结果开始形成节点导纳矩阵输入原始数据设节点电压(0)(0)iief，i=1,2…,n,is置迭代次数0k置节点号i=1按式（3-3），（3-4）计算雅克比矩阵元素按式（3-2）计算PQ节点的()kiP，()kiQ，PV节点的()kiP，()2kiU求解修正方程式，得()kie，()kif雅克比矩阵是否已全部形成？计算平衡节点及PV节点功率求()max||ke，()max||kf迭代次数k=k+1i=i+1()()maxmax||,||kkef？潮流计算完成计算各节点电压的新值：(1)()()kkkieee(1)()()kkkifffIEEE-美国新英格兰10机39节点测试系统一、系统单线图GGGGGGGGGG303912253729172693381654182728362435222120342319331011131415831126327二、系统参数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.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.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.)LN1LN2LN3LN4LN5LN6LN7LN8LN9BUS-2BUS-39BUS-3BUS-25BUS-4BUS-18BUS-5BUS-14BUS-6BUS-1BUS-1BUS-2BUS-2BUS-3BUS-3BUS-4BUS-4BUS-50.003500.001000.001300.007000.001300.001100.000800.000800.000200.041100.025000.015100.008600.021300.013300.012800.012900.002600.349350.37500.128600.073000.110700.106900.067100.069100.02170LN10LN11LN12LN13LN14LN15LN16LN17LN18LN19LN20LN21LN22LN23LN24LN25LN26LN27LN28LN29LN30LN31LN32LN33LN34BUS-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-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.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.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.073800.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接有并联电容器，B4＝1.0000LN36:BUS-5接有并联电容器，B4＝2.00005)变压器参数见表D-4表D-4变压器参数序号节点I节点JR(p.u.)X(p.u.)变比(%)TR37TR38TR39TR40TR41TR42TR43TR44TR45TR46TR47TR48BUS-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.027200.023200.015600.01380100.6100.6102.5107.0107.0100.9107.0102.5100.0102.5102.5106.0%IEEE-美国新英格兰10机39节点测试系统%123456%busvoltanglepqtypebus=[11.00000.000.000.00121.00000.000.000.00131.00000.00-3.22-0.024141.00000.00-5.00-1.84151.00000.000.000.00161.00000.000.000.00171.00000.00-2.338-0.84181.00000.00-5.22-1.76191.00000.000.000.001101.00000.000.000.001111.00000.000.000.001121.00000.00-0.085-0.881131.00000.000.000.001141.00000.000.000.001151.00000.00-3.20-1.531161.00000.00-3.29-0.3231171.00000.000.000.001181.00000.00-1.58-0.301191.00000.000.000.001201.00000.00-6.80-1.031211.00000.00-2.74-1.151221.00000.000.000.001231.00000.00-2.475-1.151241.00000.00-3.08-0.9221251.00000.00-2.24-0.4721261.00000.00-1.39-0.171271.00000.00-2.81-0.7551281.00000.00-