- 1、本文档共9页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
基于MATLAB进行潮流计算
本文介绍了基于MATLAB软件的潮流计算方法。电力系
统潮流计算方法分为手算潮流和计算机潮流计算两类。手算潮
流主要适用于规模较小的辐射型电力潮流计算,而计算机潮流
计算有两种途径:编程实现网络方程的迭代求解和借助电力系
统分析仿真软件搭建系统模型完成潮流计算。MATLAB具有
强大的矩阵运算功能和电力系统仿真平台,可以为实现潮流计
算提供更便捷的手段。本文采用极坐标形式牛顿─拉夫逊法进
行潮流计算,为其他形式的潮流计算提供借鉴。
Abstract:Thepowerflownmethodcanbedividedintotwo
categories:handnoftidalcurrentandcomputerpowerflowXXX
simplifiedequivalentcircuits。makingitXXX:programming
XXXns。orusingpowersystemXXXsystemmodelforpower
flown。MATLABarehasstrongmatrixnsanditspowersystem
XXX-Raphsonmethodofpowerflowninpolarcoordinateswith
MATLABare。andcanserveasareferenceforotherformsof
powerflown.
1.电力系统中的牛顿法潮流计算是一种常用的电力系统分
析方法。该方法基于节点电压的相等条件和潮流方程的等式条
件,通过迭代求解电压和相位的不平衡量,最终得到各节点的
电压、相位和功率等参数。
2.牛顿法潮流计算的步骤包括输入系统原始数据、形成节
点导纳矩阵、给定各节点电压初值、计算功率偏差向量、判断
收敛条件、计算雅克比矩阵、解修正方程、计算节点电压和相
位的修正值、迭代计算直至满足收敛条件、计算各节点功率等
参数并输出计算结果。
3.以一个简单电力系统为例,演示了牛顿法潮流计算的具
体步骤。该系统中包括4个节点,其中1、2为PQ节点,3为
PV节点,4为平衡节点。已知各节点的有功、无功功率和电
压等参数,通过迭代计算得到各节点的电压、相位和功率等参
数,并输出计算结果。
本次使用MATLAB软件进行了潮流计算,初步掌握了该
软件,并且坚信只要有信念,就没有学不会的东西。在大三时
研究结构力学时,也曾使用一个新的软件来分析结构的弯矩、
转矩和剪力,这次的研究也是如此。任何软件都不难,难的是
不愿意动手去实践,不愿意向他人请教。
在编程时,我按照一本书上的思路进行编程。开始时遇到
了很多不懂的问题,于是向老师和懂得MATLAB编程的同学
请教,最终通过自己的努力得到了与书上提供的结果一致的结
果。
在编程时,一定要特别细心,注意符号的使用,循环的次
数以及最终得到的结果。虽然编程过程痛苦,但结果出来的那
一刻却是无比喜悦的。
k=0;
n=1;
whilen。0.
delt(4)=0;
v(4)=1.05.%平衡节点的电压和相位
form=1:4
Pt(m)=0;
forn=1:4
Pt(m)=Pt(m)+v(m)*v(n)*(G(m,n)*cos(delt(m)-delt(n))
+B(m,n)*sin(delt(m)-delt(n)));
end
pp(m)=P(m)-Pt(m)。%节点有功功率不平衡量
end
form=1:4
Qt(m)=0;
forn=1:4
Qt(m)=Qt(m)+v(m)*v(n)*(G(m,n)*sin(delt(m)-delt(n))
-B(m,n)*cos(delt(m)-delt(n)));
end
qq(m)=Q(m)-Qt(m)。%节点无功功率不平衡量
end
form=1:3
h0(m)
文档评论(0)