资源描述：
《电力系统计算机辅助潮流计算实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、学号：200911202007电力系统计算机辅助潮流计算实验报告学生姓名：学院:电力学院系别:电力系专业:电气工程及其自动化班级:授课教师:二o年十二月（2011-2012学年第一学期）1、实验目的：了解计算机潮流分析的基本原理、主要步骤；掌握节点导纳矩阵形成和修改的方法,掌握数据处理的基本方法;熟悉Matlab运行环境,了解Matlab基本编程语句和语法;运用潮流分析程序对给定网络的运行方式做潮流分析，并初步分析计算结杲2、实验要求：通过预习，对计算机潮流分析基本理论有深入了解；为程序准备必要的、准确的原始数据；熟悉Matlab运行环
2、境，输入潮流程序，上机独立完成程序的调试，给出潮流分析的结果并按要求绘制潮流分布图3、实验内容：输入网络参数，包括节点号、节点导纳矩阵、节点功率等；输入潮流程序、调试并输出结果，绘制潮流分布图4、实验步骤：1、熟悉原始资料：根据计算要求，整理数据，包括：计算网络屮线路、变压器的参数、形成节点导纳矩阵；表示各节点的注入功率。（以上数据均采用有名值计算）2、读通潮流程序：完成程序的解释和说明，必要时附加对应的公式和程序语言的说明3、上机调试：熟悉Matlab的运行环境，准确输入原始数据、节点编号、节点注入功率等信息4、整理计算结果：根据计算
3、结果作电网潮流分布图原始网络:5.实验数据及处理:一、实验程序：clearG(l,1)=3.75;B(1,1)=-11.25;G(1,2)二一2.5;B(1,2)=7.5;G(1,3)=0;B(l,3)=0;G(l,4)=-l.25;B(1,4)=3.75;G(1,5)=0;B(l,5)=0;G(2,1)二一2.5;B(2,1)二7.5;G(2,2)=10.834;B(2,2)=-32.5;G(2,3)二一1.667;B(2,3)=5;G(2,4)二-1・667;B(2,4)=5;G(2,5)二-5;B(2,5)=15;G(3,l)=0
4、;B(3,l)=0;G(3,2)7667;B(3,2)=5;G(3,3)=12.917;B(3,3)=-38.75;G(3,4)=-10;B(3,4)=30;G(3,5)=-l.25;B(3,5)=3.75;G(4,1)=-l.25;B(4,1)=3.75;G(4,2)=-l.667;B(4,2)二5;G(4,3)=-10;B(4,3)=30;G(4,4)=12.917;B(4,4)=-38.750;G(4,5)=0;B(4,5)=0;G(5,l)=0;B(5,l)=0;G(5,2)二一5;B(5,2)=15;G(5,3)=-l.25;
5、B(5,3)=3.75;G(5,4)=0;B(5,4)=0;G(5,5)=6.25;B(5,5)=-18.75;Y二G+j*B%形成节点导纳矩阵delt(1)=0;delt(2)=0;delt(3)=0;delt(4)=0;u(l)=l.0;u(2)=l.0;u(3)=l.0;u(4)=l.0;ps(1)=-0.6;qs(l)=-0.10;ps(2)=0.2;qs(2)=0.2;ps(3)二一0.45;qs(3)二一0.15;ps(4)二一0・4;qs(4)二一0.05;%设迭代初值%设迭代次数和精度%PQ节点数Nl=4;delt(5)
6、二0;u⑸二1.06;%判断是否满足精度要求%给定平衡节点编号form=l:N1forn=l:Nl+lpt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-deIt(n)));qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));pi(m)=sum(pt);qi(m)=sum(qt);dp(m)二ps(m)-pi(m);%计算PQ节点的注入功率dq(m)=qs(m)一qi
7、(m);%计算PQ节点的功率不平衡量endform=l:N1forn=l:Nlifm==nH(m,m)=-qi(m)-u(m厂2*B(m,m);N(m,m)=pi(m)+u(m厂2*G(m,m);J(m,m)=pi(m)-u(m厂2*G(m,m);L(m,m)=qi(m)-u(m厂2*B(m,m);JJ(2*m~l,2*m-l)=H(m,m);JJ(2*mT,2*m)=N(m,m);JJ(2*m,2*m-l)=J(m,m);JJ(2*m,2*m)=L(m,m);elseH(m,n)=u(m)*u(n)*(G(m,n)*sin(delt(
8、m)-delt(n))-B(m,n)*cos(deIt(m)-delt(n)));J(m,n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt
