(完整word版)基于matlab--psat软件的电力系统潮流计算课程设计
(完整word版)电力系统潮流计算课程设计
一、问题重述课程设计要求1、在读懂程序的基础上画出潮流计算根本流程图2、经过输入数据，进行潮流计算输出结果3、对不同样的负荷变化，解析潮流分布，写出解析说明。
4、对不同样的负荷变化，进行潮流的调治控制，并说明调治控制的方法，并列表表示调治控制的参数变化。
5、打印利用 DDRTS 进行潮流解析绘制的系统图，以及潮流分布图。
课程设计题目1、系统图：两个发电厂分别经过变压器和输电线路与四个变电所相连。
变电所 1变电所 2变电所 3变电所 435kV 母线10kV 母线35kV 母线10kV 母线一次侧电压 220kV一次侧电压 220kV线路长为 60km线路长为 80km线路长为 100km线路长为 80km线路长为 80km线路长为 100km母线 1母线 2。
母线 32*QFQ-50 -22*QFS-50-22*TQN-100 -22*TQN-100 -2电厂一电厂二2、发电厂资料：母线 1 和 2 为发电厂高压母线，发电厂一总装机容量为〔400MW〕，母线 3 为机压母线，机压母线上装机容量为〔100MW〕，最大负荷和最小负荷分别为 50MW和 30MW；发电厂二总装机容量为〔 200MW〕。
3、变电所资料：〔1〕变电所 1、2、3、4 低压母线的电压等级分别为： 10KV 35KV 10KV35KV 〔2〕变电所的负荷分别为：50MW 40MW 50MW60MW〔3〕每个变电所的功率因数均为cosφ=0.85 ；〔4〕变电所 2 和变电所 4 分别配有两台容量为 75MVA的变压器，短路耗费414KW，短路电压〔 %〕=16.7 ；变电所 1 和变电所 3 分别配有两台容量为63MVA 的变压器，短路耗费为 245KW，短路电压〔 %〕=10.5 ；4、输电线路资料：发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为，单位长度的电抗为，单位长度的电纳为 2.78 * 10 -6 S 。
基于MATLAB的电力系统潮流计算
%输出计算结果
disp('节点电压为：');
通过这个程序，我们可以得到电力系统的节点电压向量。同样地，我们也可 以用节点电流法或迭代算法来求解潮流计算问题。
对于不同的算法，它们的优缺点也不尽相同。节点电压法具有计算量小、收 敛速度快等优点，但需要已知各节点的电压初始值。节点电流法相对于节点电压 法而言，收敛速度较慢，但不需要知道电压初始值。迭代算法具有普适性，可以 处理各种复杂的
基于MATLAB的电力系统潮流计算
目录
01 引言
03 Matlab工具
02 背景 04 潮流计算方法
05 结果分析
07 参考内容
目录
06 结论
引言
电力系统潮流计算是电力工程领域中非常重要的分析工具之一，用于研究电 力系统中电压、电流、功率等参数的分布和分配情况。准确地进行电力系统潮流 计算能够为电力系统的设计和运行提供重要的参考依据。本次演示将介绍使用 Matlab进行电力系统
2、利用Matlab的仿真功能，设置计算迭代的步长和算法类型等参数。
3、调用电力系统潮流计算函数， 开始计算并得到潮流结果。
4、对潮流结果进行分析和优化，为电力系统的设计和运行提供参考。
潮流计算方法
电力系统潮流计算的方法主要包括以下几个步骤：
1、网络拓扑分析：根据电力系统的结构，分析其网络拓扑关系，确定电力 系统的运行状态。
电力系统，但需要设定合适的迭代步长和初始值。
在未来研究中，我们可以进一步探索混合潮流计算方法，将不同的算法进行 组合，以获得更好的计算性能。此外，随着智能电网技术的发展，我们可以考虑 将潮流计算与优化、控制相结合，实现电力系统的智能化运行。
综上所述，基于MATLAB的电力系统潮流计算在电力工程领域具有广泛的应用 前景。通过深入研究和不断创新，我们可以进一步提高潮流计算的精度和效率， 为电力系统的稳定和经济运行提供更好的支持。
用matlab电力系统潮流计算(课程设计)
4设计图及源程序4.1MATLAB仿真相关的原始数据输入格式如下：1、B1是支路参数矩阵，第一列和第二列是节点编号。
节点编号由小到大编写。
2、对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号，将变压器的串联阻抗置于低压侧处理，第三列为支路的串列阻抗参数，第四列为支路的对地导纳参数，第五烈为含变压器支路的变压器的变比，第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，“0”为不含有变压器。
3、B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中“1”为平衡节点，“2”为PQ节点，“3”为PV 节点参数。
4、X为节点号和对地参数矩阵。
其中第一列为节点编号，第二列为节点对地参数。
4.2潮流计算源程序close allclear alln=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=input('请输入平衡节点号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入支路参数:B1=');B2=input('请输入节点参数:B2=');X=input('节点号和对地参数:X=');Y=zeros(n);Times=1; %置迭代次数为初始值%创建节点导纳矩阵for i=1:n1if B1(i,6)==0 %不含变压器的支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/B1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);else %含有变压器的支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5));Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3);YOrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1); %将OrgS、DetaS初始化%创建OrgS，用于存储初始功率参数h=0;j=0;for i=1:n %对PQ节点的处理if i~=isb&B2(i,6)==2h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j) )*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*im ag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendendfor i=1:n %对PV节点的处理，注意这时不可再将h初始化为0if i~=isb&B2(i,6)==3h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j) )*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*im ag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendendOrgS%创建PVU 用于存储PV节点的初始电压PVU=zeros(n-h-1,1);t=0;for i=1:nif B2(i,6)==3t=t+1;PVU(t,1)=B2(i,3);endend%创建DetaS，用于存储有功功率、无功功率和电压幅值的不平衡量h=0;for i=1:n %对PQ节点的处理if i~=isb&B2(i,6)==2h=h+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);endendt=0;for i=1:n %对PV节点的处理，注意这时不可再将h初始化为0if i~=isb&B2(i,6)==3h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;endendDetaS%创建I，用于存储节点电流参数i=zeros(n-1,1);h=0;for i=1:nif i~=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));endendI%创建Jacbi(雅可比矩阵)Jacbi=zeros(2*n-2);h=0;k=0;for i=1:n %对PQ节点的处理if B2(i,6)==2h=h+1;for j=1:nif j~=isbk=k+1;if i==j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendendk=0;for i=1:n %对PV节点的处理if B2(i,6)==3h=h+1;for j=1:nif j~=isbk=k+1;if i==j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=2*imag(B2(i,3));Jacbi(2*h,2*k)=2*real(B2(i,3));else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendendJacbi%求解修正方程，获取节点电压的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU%修正节点电压j=0;for i=1:n %对PQ节点处理if B2(i,6)==2j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfor i=1:n %对PV节点的处理if B2(i,6)==3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendB2%开始循环**********************************************************************while abs(max(DetaU))>prOrgS=zeros(2*n-2,1); %!!!初始功率参数在迭代过程中是不累加的，所以在这里必须将其初始化为零矩阵h=0;j=0;for i=1:nif i~=isb&B2(i,6)==2h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j) )*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*im ag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendendfor i=1:nif i~=isb&B2(i,6)==3h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j) )*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*im ag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendendOrgS%创建DetaSh=0;for i=1:nif i~=isb&B2(i,6)==2h=h+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);endendt=0;for i=1:nif i~=isb&B2(i,6)==3h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;endendDetaS%创建Ii=zeros(n-1,1);h=0;for i=1:nif i~=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));endendI%创建JacbiJacbi=zeros(2*n-2);h=0;k=0;for i=1:nif B2(i,6)==2h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1)k=0;endendendendendk=0;for i=1:nif B2(i,6)==3h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=2*imag(B2(i,3));Jacbi(2*h,2*k)=2*real(B2(i,3));elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1)k=0;endendendendendJacbiDetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU%修正节点电压j=0;for i=1:nif B2(i,6)==2j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfor i=1:nif B2(i,6)==3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendB2Times=Times+1; %迭代次数加1endTimes。
基于MATLAB的电力系统潮流计算_毕业设计论文
基于MATLAB的电⼒系统潮流计算_毕业设计论⽂基于MATLAB的电⼒系统潮流计算摘要潮流计算是电⼒系统最基本、最常⽤的计算。
根据系统给定的运⾏条件、⽹络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值及相⾓），各元件中流过的功率、整个系统的功率损耗等。
潮流计算是实现电⼒系统安全经济发供电的必要⼿段和重要⼯作环节。
因此潮流计算在电⼒系统的规划设计、⽣产运⾏、调度管理及科学研究中都有着⼴泛的应⽤。
本次设计的主要⽬的就是⾯向⼀般的电⼒⽹络，形成节点导纳矩阵，确定合适的算法，编写通⽤的计算程序，得到计算结果。
设计中主要介绍了⽜顿拉夫逊和PQ分解两种算法，PQ分解法虽然在结构上⽐⽜顿法更加简化，但是针对⼀般⽹络现代计算机在存储空间及计算速度上已经⼗分强⼤，鉴于对⽜顿法的熟悉与其算法的直观性，本次设计在编程时采⽤了⽜顿拉夫逊法的直⾓坐标形式。
解⽅程的过程利⽤Matlab的强⼤计算功能，编写M语⾔，合理设置变量，实现通⽤计算功能。
关键词: 电⼒系统，潮流计算，⽜顿—拉夫逊法，Matlab。
AbstractPower system load flow calculation is the most basic and commonly used calculations. Given according to the system operating conditions, the network connection and device parameters can be determined by power flow calculation of the bus voltage (magnitude and phase angle), the power flowing through the components, overall system power consumption and so on. Flow calculation is to achieve economic development of power system supply the necessary means and important part of the work. Therefore flow calculation in power system planning and design, production and operation, scheduling management, and scientific research have a wide range of applications.The main purpose of this design is for the general electricity network, the formation of the node admittance matrix, determine the appropriate method, the preparation of general-purpose computer program to get results. Introduces the design and the PQ decomposition Newton Raphson two algorithms, PQ decomposition although the structure is more streamlined than the Newton method, but for the general network of modern computer storage space and computing speed has been very strong, in view of the Newton Familiar with its intuitive algorithm, this design in programming using Newton Raphson polar form. The process of solving equations using matlab powerful computing capabilities, the preparation of M language, a reasonable set variables, to achieve general-purpose computing functions.Keywords: power system, power flow calculation, Newton - Raphson method, Matlab.⽬录摘要 (I)Abstract ................................................................................................................................................ II ⽬录.................................................................................................................................................... I II 1 引⾔ .. (1)1.1 潮流计算⽬的 (1)1.2 潮流计算意义 (1)1.3 潮流计算发展史 (1)1.4基于MATLAB 的电⼒系统潮流计算发展前景 (2)2简单电⼒系统潮流计算 (4)2.1简单辐射⽹络的潮流计算 (4)2.1.1简单⽀路的潮流分布和电压降落 (4)2.1.2 辐射型⽹络的⼿⼯潮流计算⽅法 (6)2.2 简单环⽹的潮流计算 (7)2.2.1两端电压相等 (7)2.2.2两端电压不相等 (8)3 复杂电⼒系统潮流计算的计算机算法 (10)3.1电⼒⽹络⽅程及等值电路 (10)3.2节点导纳矩阵形成及修改 (11)3.3节点的分类 (14)3.3.1 PQ节点 (14)3.3.2 PV节点 (14)3.3.3 平衡节点 (14)3.4潮流计算的约束条件 (15)3.5⽜顿-拉夫逊法（直⾓坐标） (15)3.5.1⽜顿-拉夫逊法的推导过程 (15)3.5.2潮流计算时的修正⽅程（直⾓坐标） (17) 3.5.3雅可⽐矩阵的特点: (19)3.5.4⽜顿-拉夫逊法计算步骤 (19)3.6 P-Q分解法潮流计算 (20)3.6.1 P-Q分解法潮流计算概述 (20)3.6.2 P-Q分解法的潮流计算步骤 (20)3.6.3 P-Q分解法的特点 (21)4 Matlab概述 (22)4.1Matlab简介 (22)4.2 Matlab中的变量 (22)4.3 Matlab编程 (23)4.3.1矩阵的输⼊ (23)4.3.2矩阵的运算 (24)4.3.3 MatLab的控制流 (24)5 ⽜顿法潮流计算程序设计及实例 (26)5.1⼿算 (26)5.2计算机算法的数据输⼊ (29)5.3潮流计算程序 (30)5.3 计算结果分析 (36)结论 (37)参考⽂献 (38)附录A 程序流程图 (39)附录B Matlab仿真 (40)致谢 (1)1 引⾔1.1 潮流计算⽬的电⼒系统潮流计算是研究电⼒系统稳态运⾏情况的⼀种基本电⽓计算。
用matlab电力系统潮流计算课程设计
湖北民族学院信息工程学院课程设计报告书题目: 电力系统潮流计算专业：电气工程及其自动化班级：学号：学生姓名：指导教师：2011年 6 月 10 日2011年6月12日信息工程学院课程设计成绩评定表摘要电力系统稳态分析包括潮流计算和静态安全分析。
本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。
本位就是运用潮流计算具体分析，并有MATLAB仿真。
关键词：电力系统潮流计算 MATLAB仿真AbstractElectric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation.Key words: Power system; Flow calculation; MATLAB simulation目录1 任务提出与方案论证 (2)2 总体设计 (3)2.1潮流计算等值电路 (3)2.2建立电力系统模型 (3)2.3模型的调试与运行 (3)3 详细设计 (4)3.1 计算前提 (4)3.2手工计算 (7)4设计图及源程序 (11)4.1MATLAB仿真 (11)4.2潮流计算源程序 (11)5 总结 (19)参考文献 (20)1 任务提出与方案论证潮流计算是在给定电力系统网络结构、参数和决定系统运行状态的边界条件的情况下确定系统稳态运行状态的一种基本方法,是电力系统规划和运营中不可缺少的一个重要组成部分。
基于MATLAB的电力系统潮流计算课程设计
基于MATLAB的电力系统潮流计算课程设计课程设计报告学生姓名:学号：学院: 电气工程学院班级:题目: 基于MATLAB的电力系统潮流计算指导教师：职称:职称：2014年 1 月 10 日目录：一、电力系统潮流计算概述 (5)二、潮流计算方法概述 (5)2．1节点分类 (5)2.2牛顿—拉夫逊法概要 (6)2.2.1牛顿—拉夫逊法迭代原理 (6)2．2.1 牛顿法的框图及求解过程 (7)三、课程设计任务 (9)3.1题目原始数据 (10)3.2课程设计要求 (11)四、初步分析 (11)4.1节点设置及分类 (11)4.2 参数计算 (11)4.3等值电路绘制 (13)五、潮流计算 (14)5.1 给定负荷下的潮流计算 (14)5.1.1 B1/B2矩阵的形成 (14)5.1.2 潮流调整及分析 (14)5.1.3 仿真比较 (16)5.2 变电所负荷变化时的潮流计算 (18)5.2.1 4个变电所负荷同时以2%的比例增大 (18)5.2.2 4个变电所负荷成比例2%下降 (20)5.2.3 1、4号负荷下降，2、3号负荷上升 (23)5.2.4 仿真比较 (25)5.3断线潮流计算 (27)5.5.1 断开1、3节点间的一条支路 (27)5.3.2 断开1、5支路的一条线 (29)5.3.3 断开1、7支路的一条线 (31)5.3.4 断开2、9支路的一条线 (33)5.3.5 断开7、9支路的一条线 (35)5.3.6 仿真比较 (37)六、潮流分析总结 (39)七、心得体会 (41)参考文献： (42)附录： (42)潮流计算课程设计摘要电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压。
各元件中流过的功率，系统的功率损耗等等。
在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量的分析比较供电方案或运行方式的合理性、可靠性和经济性。
电力系统潮流计算的MATLAB辅助程序设计,潮流计算程序(精编文档).doc
【最新整理，下载后即可编辑】电力系统潮流计算的MATLAB辅助程序设计潮流计算，通常指负荷潮流，是电力系统分析和设计的主要组成部分，对系统规划、安全运行、经济调度和电力公司的功率交换非常重要。
此外，潮流计算还是其它电力系统分析的基础，比如暂态稳定，突发事件处理等。
现代电力系统潮流计算的方法主要:高斯法、牛顿法、快速解耦法和MATLAB的M语言编写的MATPOWER4.1，这里主要介绍高斯法、牛顿法和快速解耦法。
高斯法的程序是lfgauss，其与lfybus、busout和lineflow程序联合使用求解潮流功率。
lfybus、busout和lineflow程序也可与牛顿法的lfnewton程序和快速解耦法的decouple程序联合使用。
(读者可以到MATPOWER主页下载MATPOWER4.1,然后将其解压到MATLAB目录下，即可使用该软件进行潮流计算)一、高斯-赛德尔法潮流计算使用的程序：高斯-赛德法的具体使用方法读者可参考后面的实例，这里仅介绍各程序的编写格式:lfgauss：该程序是用高斯法对实际电力系统进行潮流计算，需要用到busdata和linedata两个文件。
程序设计为输入负荷和发电机的有功MW和无功Mvar，以及节点电压标幺值和相角的角度值。
根据所选复功率为基准值将负荷和发电机的功率转换为标幺值。
对于PV节点，如发电机节点，要提供一个无功功率限定值。
当给定电压过高或过低时，无功功率可能超出功率限定值。
在几次迭代之后(高斯-塞德尔迭代为10次)，需要检查一次发电机节点的无功出力，如果接近限定值，电压幅值进行上下5%的调整，使得无功保持在限定值内。
lfybus：这个程序需要输入线路参数、变压器参数以及变压器分接头参数。
并将这些参数放在名为linedata的文件中。
这个程序将阻抗转换为导纳，并得到节点导纳矩阵。
busout：该程序以表格形式输出结果，节点输出包括电压幅值和相角，发电机和负荷的有功和无功功率，以及并联电容器或电抗器的有功和无功功率。
《matlab潮流计算》word版
%线路数据Line=[1-节点编号2-线路首端节点号3-线路末端节点号4-支路电阻5-支路电抗6-支路电纳(注意：此处取的是B/2)]Line=[1 1 2 0.04 0.25 0.252 13 0.1 0.35 03 2 3 0.08 0.30 0.25 ];%变压器数据transform=[1-支路编号2-支路首节点编号3-支路末节点编号4-支路电阻(p.u.) 5-支路电抗(p.u.) 6-变压器变比(p.u.)]transform=[ 1 2 4 0 0.015 1.052 3 5 0 0.03 1.05];%数据预处理Nbus=5; % 节点数nline=size(Line,1); % 线路个数ntrans=size(transform,1); % 变压器个数slack=5; % 平衡节点号Npq=4; % PQ节点的个数Ss=0;%计算节点导纳矩阵Y=zeros(Nbus);if nline>=1%判断是否存在线路for k=1:nline %以下处理线路t1=Line(k,2); t2=Line(k,3); b2=Line(k,6);%分别取出线路的首端节点编号t1、末端节点编号t2和对地电纳b2Yl=1/(Line(k,4)+j*Line(k,5));%计算线路的支路电导YlY(t1,t1)=Y(t1,t1)+Yl+j*b2;%修正第k条线路首端节点的自导纳Y(t1,t2)=Y(t1,t2)-Yl;%修正第k条线路首端节点与末端节点之间的互导纳Y(t2,t1)=Y(t2,t1)-Yl;%修正第k条线路末端节点与首端节点之间的互导纳Y(t2,t2)=Y(t2,t2)+Yl+j*b2;%修正第k条线路末端节点的自导纳endendif ntrans>=1%判断是否存在变压器for k=1:ntrans %以下处理变压器t1=transform(k,2);t2=transform(k,3);t3=transform(k,6);%分别取出变压器的首端节点编号t1、末端节点编号t2和变比t3Yt=1/(transform(k,4)+j*transform(k,5));Yt1=Yt/t3;Yt2=Yt*(1-t3)/(t3*t3);Yt3=Yt*(t3-1)/t3;Y(t1,t1)=Y(t1,t1)+Yt1+Yt2;Y(t1,t2)=Y(t1,t2)-Yt1;Y(t2,t1)=Y(t2,t1)-Yt1;Y(t2,t2)=Y(t2,t2)+Yt1+Yt3;endendG=real(Y);B=imag(Y); %区分节点导纳矩阵的实部和虚部GB%赋初值delt(1)=0;delt(2)=0;delt(3)=0;delt(4)=0;u(1)=1;u(2)=1;u(3)=1;u(4)=1;p(1)=-0.30;q(1)=-0.18;p(2)=-0.55;q(2)=-0.13;p(3)=0;q(3)=0;p(4)=0.5;q(4)=1.10;p(5)=0.8;q(5)=0.50;k=0;precision=1;Npq=4; %Npq分别是网络中的PQ节点数%[Unbalance]=-[Jacobi][Correction]while precision>0.00001 %设定误差上限，判断是否继续迭代u(5)=1.06; delt(5)=0; %设定平衡节点电压相角与幅值k; u; delt;for m=1:Npqfor n=1:Nbuspt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n))); %由节点电压求得的PQ节点注入有功功率qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n))); %由节点电压求得的PQ节点注入无功功率endUnbalance(2*m-1)=p(m)-sum(pt); %计算PQ节点有功功率不平衡量Unbalance(2*m )=q(m)-sum(qt); %计算PQ节点无功功率不平衡量end %[Unbalance]是节点不平衡量矩阵for m=1:Npqfor n=1:Nbush0(n)= u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));l0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));endH(m,m)=sum(h0)-u(m)^2*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));N(m,m)=sum(n0)+u(m)^2*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)))-2*u(m)^2*G(m,m);J(m,m)=sum(j0)+u(m)^2*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)));L(m,m)=sum(l0)+u(m)^2*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)))+2*u(m)^2*B(m,m);Jacobi(2*m-1,2*m-1)=H(m,m);Jacobi(2*m-1,2*m )=N(m,m);Jacobi(2*m ,2*m-1)=J(m,m);Jacobi(2*m ,2*m )=L(m,m);end %计算m=n情况下的Jacobi矩阵中的子矩阵元素for m=1:Npqfor n=1:Npqif m==nelseH(m,n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));J(m,n)= u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));N(m,n)=-J(m,n);L(m,n)= H(m,n);Jacobi(2*m-1,2*n-1)=H(m,n);Jacobi(2*m-1,2*n )=N(m,n);Jacobi(2*m ,2*n-1)=J(m,n);Jacobi(2*m ,2*n )=L(m,n);endendend %计算m≠n情况下的Jacobi矩阵中的子矩阵元素Correction=-Jacobi\(Unbalance'); %计算电压相角和幅值的修正量precision = max(abs(Correction)); %取误差最大值for m=1:Npqdelt(m)=delt(m)+Correction(2*m-1); %修正PQ节点电压相角u(m)= u(m)+Correction(2*m ); %修正PQ节点电压幅值endk=k+1; %迭代轮数+1endk, u, delt, Jacobi, precisionfor m=1:NbusU(m)=u(m)*(cos(delt(m))+j*sin(delt(m))); %采用直角坐标系表示电压I(m)=Y(Nbus,m)*U(m); %计算注入平衡节点的电流endSslack=U(Nbus)*sum(conj(I)) %计算注入平衡节点的功率for m=1:Nbusfor n=1:NbusS(m,n)=U(m)*(conj(U(m))-conj(U(n)))*conj(-Y(m,n)); %计算线路功率Ploss(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));endendS %显示线路功率矩阵，S(m,n)表示假设功率由节点m流向节点n%若数值为+则说明实际功率流向与假设方向相同，若数值为-则说明实际功率流向与假设方向相反。
Matlab电力系统潮流仿真计算课程设计报告
课程设计报告学生姓名:学号：学院:班级:题目:指导教师：职称: 指导教师：职称:2019年月日目录1 绪论 (1)1.1 潮流计算在电力系统分析的作用及意义 (1)1.1.1 潮流计算的作用 (1)1.1.2潮流计算的意义 (1)1.2 潮流计算的主要方法 (2)1.2.2牛顿-拉夫逊法 (2)1.2.3 P-Q分解法 (3)1.3 课程设计主要工作 (6)2 计算原理与流程 (7)2.1 潮流计算数学模型 (7)2.1.1变量及节点分类 (7)2.1.2各元件数学模型 (7)2.1.3潮流方程 (8)2.2 计算方法的数学原理 (9)2.2.1基本原理 (9)2.2.2修正方程 (10)2.2.3牛顿—拉夫逊法的求解过程及结构框图 (14)2.3 潮流计算的约束条件 (15)3 仿真工具说明 (17)3.1 MATLAB说明 (17)3.2 MATPOWER说明 (17)3.2.1 电力常规潮流运算 (17)3.2.2 数据文件格式 (17)3.2.3 潮流计算 (19)4 仿真工具应用说明 (20)4.1 IEEE标准数据格式 (20)4.2 潮流计算结果 (21)4.2.1 牛顿-拉夫逊法 (22)4.2.2 P-Q分解法 (23)4.2.3 方法比较 (23)4.3潮流情况 (23)4.3.1电压情况 (24)4.3.2线路情况 (25)4.4电压调整 (27)4.1.1改变发电机端电压调压 (27)4.4.2改变变压器的变比调压 (28)4.5 提高节点电压水平 (29)5 总结 (30)参考文献 (31)1 绪论1.1 潮流计算在电力系统分析的作用及意义1.1.1 潮流计算的作用根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线的电压、各支路流过的功率、整个系统的功率损耗，用以研究系统规划和运行中提出的各种问题。
对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。
matalb潮流计算课程设计
matalb潮流计算课程设计一、教学目标本课程旨在通过MATLAB潮流计算的学习，使学生掌握电力系统潮流计算的基本原理和方法，能够熟练使用MATLAB软件进行潮流计算，并分析计算结果。
具体目标如下：1.掌握电力系统潮流计算的基本概念和原理。
2.理解电力系统潮流计算的数学模型和算法。
3.熟悉MATLAB软件的基本操作和编程方法。
4.能够使用MATLAB软件进行电力系统潮流计算。
5.能够分析潮流计算结果，判断系统运行状态。
6.能够根据实际情况，调整计算参数，优化系统运行。
情感态度价值观目标：1.培养学生的科学思维和解决问题的能力。
2.培养学生的团队合作和创新精神。
3.培养学生的责任感和使命感，使其能够为电力系统的安全稳定运行做出贡献。
二、教学内容本课程的教学内容主要包括以下几个部分：1.电力系统潮流计算的基本概念和原理。
2.电力系统潮流计算的数学模型和算法。
3.MATLAB软件的基本操作和编程方法。
4.使用MATLAB软件进行电力系统潮流计算的实践操作。
5.潮流计算结果的分析方法和技巧。
6.电力系统潮流计算的基本概念和原理（2课时）。
7.电力系统潮流计算的数学模型和算法（2课时）。
8.MATLAB软件的基本操作和编程方法（3课时）。
9.使用MATLAB软件进行电力系统潮流计算的实践操作（4课时）。
10.潮流计算结果的分析方法和技巧（2课时）。
三、教学方法为了提高学生的学习兴趣和主动性，本课程将采用多种教学方法相结合的方式进行教学，包括：1.讲授法：用于讲解电力系统潮流计算的基本概念、原理和算法。
2.讨论法：用于讨论MATLAB软件的使用方法和编程技巧。
3.案例分析法：用于分析实际电力系统的潮流计算案例。
4.实验法：用于学生亲自动手进行MATLAB软件的操作和潮流计算实践。
四、教学资源为了支持教学内容和教学方法的实施，丰富学生的学习体验，我们将准备以下教学资源：1.教材：《电力系统潮流计算》相关教材。
(完整word版)基于MATLAB进行潮流计算
基于MATLAB进行潮流计算学生：王仕龙2011148213指导老师：李咸善摘要：电力系统潮流计算方法有两类，即手算潮流和计算机潮流计算。
手算潮流主要借助于形成简化的等值电路来实现，这种方法尤其适用于规模不大的辐射型电力潮流计算。
计算机潮流计算的实现有两种途径：其一是编程实现网络方程的迭代求解；其二是借助与电力系统分析仿真软件，搭建系统模型来完成潮流计算。
MATLAB具有强大的矩阵运算功能，同时其具有电力系统仿真平台也为直观地实现潮流计算提供了更便捷的手段[1]。
本文是基于MATLAB软件，采用极坐标形式牛顿─拉夫逊法进行潮流计算，为其他形式的潮流计算有借鉴的作用。
关键词：电力系统；计算机潮流计算；MATLAB ；牛顿─拉夫逊法Abstract:The power flow calculation method has two kinds,which are the hand calculation of tidal current and computer power flow calculation.Hand calculation tidal current is mainly realized by means of the formation of simplified equivalent circuit.This method is especially suitable for small scale radiation power flow calculation.There are two ways to realize the computer power flow calculation.The first one is through the programming iteration for solving network equation,the second one is with the help of analysis of power system simulation software to build the system model complete the power flow calculation.The software of MATLAB has strong matrix function,.At the same time,It’s power system simulation platform provides a more convenient means to realize power flow calculation intuitively[1].This paper is based on the software of MATLAB to calculate the power flow calculation by adopting the form of Newton-Raphson method of power flow calculation of polar coordinates.And it can be the role of reference of other forms of power flow calculation.Key words: power system computer; power flow calculation;MATLAB;Newton-Raphson1.计算原理电力系统潮流是指系统中所有运行参数的总体，包括各个母线电压的大小和相位，各个发电机和负荷的功率及电流，以及各个变压器和线路等元件所通过的功率,电流和其中的损耗。
基于matlab的潮流计算
课程设计任务书第学期学院：信息商务学院专业：电气工程及其自动化学生姓名：学号课程设计题目：基于MATLAB的潮流计算起迄日期：月日～月日课程设计地点：电气工程系工程实验室指导教师：系主任：课程设计任务书课程设计任务书1 引言潮流计算是电力系统最基本最常用的计算。
根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值和相角），各支路流过的功率，整个系统的功率损耗。
潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。
因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。
潮流计算在数学上是多元非线性方程组的求解问题，牛顿—拉夫逊Newton-Raphson法是数学上解非线性方程组的有效方法，有较好的收敛性。
运用电子计算机计算一般要完成以下几个步骤：建立数学模型，确定解算方法，制订计算流程，编制计算程序。
1.1 设计任务与要求利用MATLAB/SIMULINK仿真软件，搭建多电压等级电力网络潮流计算模型。
掌握各种潮流计算方法，掌握MATLAB/SPS工具箱使用方法及其在潮流计算中的应用。
MATLAB作为先进的仿真软件，在电力系统潮流计算方面，尤其是矩阵运算上具有突出的优秀性。
本设计要求利用MATLAB编写程序，输入数据并完成潮流分布计算。
1.2实用价值与理论意义在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。
要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。
为了保证电力系统的稳定运行，要进行潮流调节。
随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。
电力系统潮流计算是电力系统分析计算中最基本的内容，也是电力系统运行及设计中必不可少的工具。
根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线电压的幅值及相角、各元件中流过的功率、整个系统的功率损耗等。
matlab实验电力系统潮流计算
实验一 电力系统潮流计算一、一元非线性方程求解例1-1 试求非线性方程 )(x f =0 的解。
解：（1）取一个合理的初值)0(x 作为方程)(x f =0的解，如果正好0)()0(=x f ，则方程的解*x =)0(x 。
否则做下一步。
（2）取)0()0(x x ∆+为第一次修正值。
)0(x ∆充分小，将)()0()0(x x f ∆+在)0(x 附近展开成泰勒级数，并且将的高次项略去，取其线性部分，得到)()0()0(x x f ∆+≈)0()0()0()()(x x f x f ∆'+=0 （1-1）上式表明，在)0(x 处把非线性方程0)(=x f 线性化，变成求)0(x 附近修正量)0(x ∆的线性方程，这个方程也称为修正方程式。
从而可求得)()()0()0()0(x f x f x'-=∆ （1-2） 所以，可以确定第一次修正值)0()0()1(x x x ∆+=。
若0)()1(=x f ，则)1(*x x =。
（3）若0)()1(≠x f ，则用步骤（2）阐述的方法由)1(x 确定出第二次修正值)2(x。
如此迭代下去，在第)1(+k 次迭代时，)1(+k x应为)()()1(k k k xxx∆+=+=)()()()()(k k k x f x f x'- （1-3） 式中k 为迭代次数。
如果ε<+)()1(k x f （ε是预设的一个小的正数，如510-=ε），则方程的解)1(*+=k xx ，迭代停止。
例1-2 应用牛顿—拉夫逊法求解非线性方程0122)(23=-+-=x x x x f解：设初始近似解0.2)0(=x，首先根据（1-1）计算)()0(x f10)()0(-=x f然后计算)()0(x f '5)()0(='x f根据（1-2）式计算)0(x ∆2510)()()0()0()0(=--='-=∆x f x f x再根据（1-3）式计算)1(x ∆422)0()0()1(=+=∆+=x x x重复以上计算直到5)1(10)(-+<k x f ，得到的计算过程量和结果见表1-1。
基于MATLAB和PSASP的电力系统潮流分析与计算
基于MATLAB和PSASP的电力系统潮流分析与计算
本报告将介绍一个基于MATLAB和PSASP的电力系统潮流分析和计算的过程。
电力系统潮流分析主要用于分析电力系统的物理特性。
本文将描述如何使用MATLAB和PSASP软件进行潮流计算和分析，并对计算结果进行解释。
MATLAB是一种高级编程语言，它通常用于编写数学程序和模拟。
它可以用来分析电力系统的潮流情况。
MATLAB可以帮助用户构建模型，并分析其潮流状态。
MATLAB可以根据电力系统中每个组件的电压、电流和功率数据，来计算整个系统的潮流情况。
PSASP是一个用于电力系统潮流分析和计算的强大工具。
它可以用来分析电力系统的潮流、电压和功率因素。
PSASP可以用于计算电力系统的潮流，并分析电力系统的安全性和优化性能。
本文将使用MATLAB和PSASP软件分析和计算电力系统的潮流。
使用MATLAB编写的程序，可以轻松构建电力系统的模型，并计算每个组件的电压、电流和功率数据。
通过将这些数据导入PSASP，可以对电力系统进行潮流分析。
由于PSASP可以分析潮流、电压和功率因素，因此它可以用来诊断电力系统的故障，以提高系统的安全性和运行效率。
因此，通过使用MATLAB和PSASP，可以对电力系统进行有效的潮流分析和计算。
通过这种分析，可以帮助用户更好地了
解电力系统的特性，更有效地优化系统的运行，并保持电力系统的安全性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
东北电力大学课程设计改革试用任务书:
电力系统潮流计算课程设计任务书
设计名称：电力系统潮流计算课程设计
设计性质：理论计算，计算机仿真与验证
计划学时：两周
一、设计目的
1.培养学生独立分析问题、解决问题的能力；
2.培养学生的工程意识，灵活运用所学知识分析工程问题的能力
3.编制程序或利用电力系统分析计算软件进行电力系统潮流分析。
二、原始资料
1、系统图：IEEE14节点。
2、原始资料：见IEEE14节点标准数据库
三、课程设计基本内容：
1.采用PSAT仿真工具中的潮流计算软件计算系统潮流；
1)熟悉PSAT仿真工具的功能；
2)掌握IEEE标准数据格式内容；
3)将IEEE标准数据转化为PSAT计算数据；
2.分别采用NR法和PQ分解法计算潮流，观察NR法计算潮流中雅可比矩阵的变化情况，
分析两种方法计算潮流的优缺点；
3.分析系统潮流情况，包括电压幅值、相角，线路过载情况以及全网有功损耗情况。
4.选择以下内容之一进行分析：
1)找出系统中有功损耗最大的一条线路，给出减小该线路损耗的措施，比较各种措施
的特点，并仿真验证；
2)找出系统中电压最低的节点，给出调压措施，比较各种措施的特点，并仿真验证；
3)找出系统中流过有功功率最大的一条线路，给出减小该线路有功功率的措施，比较
各种措施的特点，并仿真验证；
5.任选以下内容之一作为深入研究：（不做要求）
1)找出系统中有功功率损耗最大的一条线路，改变发电机有功出力，分析对该线路有
功功率损耗灵敏度最大的发电机有功功率，并进行有效调整，减小该线路的损耗；
2)找出系统中有功功率损耗最大的一条线路，进行无功功率补偿，分析对该线路有功
功率损耗灵敏度最大的负荷无功功率，并进行有效调整，减小该线路的损耗；
3)找出系统中电压最低的节点，分析对该节点电压幅值灵敏度最大的发电机端电压，
并有效调整发电机端电压，提高该节点电压水平；
四、课程设计成品基本要求：
1.绘制系统潮流图，潮流图应包括：
1)系统网络参数
2)节点电压幅值及相角
3)线路和变压器的首末端有功功率和无功功率
2.撰写设计报告，报告内容应包括以下几点：
1)本次设计的目的和设计的任务；
2)电力系统潮流计算的计算机方法原理，分析NR法和PQ分解法计算潮流的特点；
3)对潮流计算结果进行分析，评价该潮流断面的运行方式安全性和经济性；
4)找出系统中运行的薄弱环节，如电压较低点或负载较大线路，给出调整措施；
5)分析各种调整措施的特点并比较它们之间的差异；
6)结论部分以及设计心得；
五、考核形式
1.纪律考核：学生组织出勤情况和工作态度等；
2.书面考核：设计成品的完成质量、撰写水平等；
3.答辩考核：参照设计成品，对计算机方法进行电力系统潮流计算的相关问题等进行答辩；
4.采用五级评分制：优、良、中、及格、不及格五个等级。
原版任务书:
电力系统潮流计算课程设计任务书
一 .题目原始资料
1、系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。
2、发电厂资料：
母线1和2为发电厂高压母线，发电厂一总装机容量为（ 400MW ），母线3为机压母线，机压母线上装机容量为（ 100MW ），最大负荷和最小负荷分别为50MW 和30MW ；发电厂二总装机容量为（ 200MW ）。
3、变电所资料：
(1) 变电所1、2、3、4低压母线的电压等级分别为：10KV 35KV 10KV 35KV
(2) 变电所的负荷分别为：
50MW 70MW 40MW 70MW
(3) 每个变电所的功率因数均为cos φ=0.85；
(4) 变电所2和变电所4分别配有两台容量为75MVA 的变压器，短路损耗414KW ，短路电
压（%）=16.7；变电所1和变电所3分别配有两台容量为63MVA 的变压器，短路损耗为245KW ，短路电压（%）=10.5；
4、输电线路资料：
发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为 17.0，单位长
变电所1 变电所2 变电所3 变电所4 35kV 母线 10kV 母线
35kV 母线 10kV 母线 一次侧电压220kV 一次侧电压220kV 母线1 母线3 母线2 线路长为100km 线路长为80km 线路长为60km
线路长为80km
线路长为100km 2*QFQ -50-2 2*QFS -50-2 2*TQN -100-2
2*TQN -100-2 。
线路长为80km 。
电厂一 电厂二
度的电抗为 0.402，单位长度的电纳为S -610*2.78。
六、 课程设计基本内容：
4. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。
5. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮流计
算，并对计算结果进行分析。
6. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。
1) 4个变电所的负荷同时以2%的比例增大；
2) 4个变电所的负荷同时以2%的比例下降
3) 1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以2%
的比例上升；
7. 在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要求，进行电压的
调整。
（变电所低压母线电压10KV 要求调整范围在9.5-10.5之间；电压35KV 要求调整范围在35-36之间）
8. 轮流断开环网一回线，分析潮流的分布。
9. 利用DDRTS 软件，进行绘制系统图进行上述各种情况潮流的分析，并进行结果的比较。
10. 最终形成课程设计成品说明书。
三、课程设计成品基本要求：
1. 在读懂程序的基础上画出潮流计算基本流程图
2. 通过输入数据，进行潮流计算输出结果
3. 对不同的负荷变化，分析潮流分布，写出分析说明。
4. 对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并列表表示调节控
制的参数变化。
5. 打印利用DDRTS 进行潮流分析绘制的系统图，以及潮流分布图。
以下是基于原图，将线路由轻载改为重载后的课程设计：。