电力系统分析潮实验报告.docVIP

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
电力系统分析潮实验报告

南昌大学实验报告 学生姓名: 学 号: 专业班级: 实验类型:□ 验证 □ 综合 ■ 设计 □ 创新 实验日期: 实验成绩: 电力系统潮流计算实验 一、实验目的: 本实验通过对电力系统潮流计算的计算机程序的编制与调试,获得对复杂电力系统进行潮流计算的计算机程序,使系统潮流计算能够由计算机自行完成,即根据已知的电力网的数学模型(节点导纳矩阵)及各节点参数,由计算程序运行完成该电力系统的潮流计算。通过实验教学加深学生对复杂电力系统潮流计算计算方法的理解,学会运用电力系统的数学模型,掌握潮流计算的过程及其特点,熟悉各种常用应用软件,熟悉硬件设备的使用方法,加强编制调试计算机程序的能力,提高工程计算的能力,学习如何将理论知识和实际工程问题结合起来。 二、实验内容: 编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电力网的数学模型(节点导纳矩阵)及各节点参数,完成该电力系统的潮流计算,要求计算出节点电压、功率等参数。 在各种潮流计算的算法中选择一种,按照计算方法编制程序。 将事先编制好的电力系统潮流计算的计算程序原代码由自备移动存储设备导入计算机。 在相应的编程环境下对程序进行组织调试。 应用计算例题验证程序的计算效果。 三、实验程序: function [e,f,p,q]=flow_out(g,b,kind,e,f) %计算潮流后efpq的终值 s=flow(g,b,kind,e,f); k=0; while max(abs(s))10^-5 J=J_out(g,b,kind,e,f); J_ni=inv(J); dv=J_ni*s; l=length(dv)/2; for i=1:l e(i)=e(i)-dv(2*i-1); f(i)=f(i)-dv(2*i); end s=flow(g,b,kind,e,f); end l=length(e); for i=1:l s1=0; s2=0; for j=1:l s1=s1+g(i,j)*e(j)-b(i,j)*f(j); s2=s2+g(i,j)*f(j)+b(i,j)*e(j); end p(i)=e(i)*s1+f(i)*s2; q(i)=f(i)*s1-e(i)*s2; end function s=flow(g,b,kind,e,f) %计算当前ef与规定的pqv的差值 l=length(e); s=zeros(2*l-2,1); for i=1:(l-1) s1=0; s2=0; for j=1:l s1=s1+g(i,j)*e(j)-b(i,j)*f(j); s2=s2+g(i,j)*f(j)+b(i,j)*e(j); end s(2*i-1)=kind(2,i)-e(i)*s1-f(i)*s2; if kind(1,i)==1 s(2*i)=kind(3,i)-f(i)*s1+e(i)*s2; else s(2*i)=kind(3,i)^2-f(i)^2-e(i)^2; end end function J=J_out(g,b,kind,e,f) %计算节点的雅克比矩阵 l=length(e); J=zeros(2*l-2,2*l-2); for i=1:(l-1); if kind(1,i)==1 s=PQ_out(g,b,e,f,i); for j=1:(2*l-2) J(2*i-1,j)=s(1,j); J(2*i,j)=s(2,j); end else s=PV_out(g,b,e,f,i); for j=1:(2*l-2) J(2*i-1,j)=s(1,j); J(2*i,j)=s(2,j); end end end function pq=PQ_out(g,b,e,f,i) %计算pq节点的雅克比矩阵 l=length(e); pq=zeros(2,2*l-2); for j=1:(l-1) if j==i s=0; for k=1:l s=s-(g

您可能关注的文档

文档评论(0)

mk808606 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档