电力系统计算机辅助潮流计算实验报告.docVIP

（ 二 〇 一 一 年 十 二 月 （2011-2012学年第一学期） 1、实验目的： 了解计算机潮流分析的基本原理、主要步骤；掌握节点导纳矩阵形成和修改的方法，掌握数据处理的基本方法；熟悉Matlab运行环境，了解Matlab基本编程语句和语法；运用潮流分析程序对给定网络的运行方式做潮流分析，并初步分析计算结果 2、实验要求： 通过预习，对计算机潮流分析基本理论有深入了解；为程序准备必要的、准确的原始数据；熟悉Matlab运行环境，输入潮流程序，上机独立完成程序的调试，给出潮流分析的结果并按要求绘制潮流分布图 3、实验内容： 输入网络参数，包括节点号、节点导纳矩阵、节点功率等；输入潮流程序、调试并输出结果，绘制潮流分布图 4、实验步骤： 1、熟悉原始资料：根据计算要求，整理数据，包括：计算网络中线路、变压器的参数、形成节点导纳矩阵；表示各节点的注入功率。（以上数据均采用有名值计算） 2、读通潮流程序：完成程序的解释和说明，必要时附加对应的公式和程序语言的说明 3、上机调试：熟悉Matlab的运行环境，准确输入原始数据、节点编号、节点注入功率等信息 4、整理计算结果：根据计算结果作电网潮流分布图 原始网络： 5、实验数据及处理： 一、实验程序： clear G(1,1)=3.75;B(1,1)=-11.25;G(1,2)=-2.5;B(1,2)=7.5;G(1,3)=0;B(1,3)=0; G(1,4)=-1.25;B(1,4)=3.75;G(1,5)=0;B(1,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,1)=0;B(3,1)=0;G(3,2)=-1.667;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)=-1.25;B(3,5)=3.75; G(4,1)=-1.25;B(4,1)=3.75;G(4,2)=-1.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,1)=0;B(5,1)=0;G(5,2)=-5;B(5,2)=15;G(5,3)=-1.25;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(1)=1.0;u(2)=1.0;u(3)=1.0;u(4)=1.0; ps(1)=-0.6;qs(1)=-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; %设迭代初值 k=1;precision=1 %设迭代次数和精度 N1=4; %PQ节点数 while precision0.0000001 %判断是否满足精度要求 delt(5)=0;u(5)=1.06; %给定平衡节点编号 for m=1:N1 for n=1:N1+1 pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n))); qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n))); end pi(m)=sum(pt);qi(m)=sum(qt); %计算PQ节点的注入功率 dp(m)=ps(m)-pi(m); dq(m)=qs(m)-qi(m); %计算P