完整word版基于matlab--psat软件的电力系统潮流计算课程设计

课程设计课程名称：电力系统分析设计题目：基于Matlab计算程序的电力系统运行分析学院：电力工程学院专业：电气工程自动化年级： 2011教务处制目录任务书 (1)前言 (2)第一章电网模型的建立 (3)第二章潮流计算 (8)第三章故障电流计算 (18)第四章思考题 (31)总结体会 (34)参考文献 (35)前言电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。

潮流计算的目标是求取电力系统在给定运行状态的计算。

即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。

对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。

潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。

在电力系统中可能发生的各种故障中，危害最大且发生概率较高的首推短路故障。

产生短路故障的主要原因是电力设备绝缘损坏。

短路故障分为三相短路、两相短路、单相接地短路及两相接地短路。

其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路。

短路故障大多数发生在架空输电线路。

电力系统设计与运行时，要采取适当的措施降低短路故障的发生概率。

短路计算可以为设备的选择提供原始数据。

第一章电网模型的建立一、目标电网接线图二、目标电网数据系统参数表1. 线路参数表线路编号线路型号线路长度（km）线路电阻{Ω/km} 线路正序电抗{Ω/km}线路容纳之半{S/km}4-5 LGJ-240/30 113 0.0470.4 1.78×4-6 LGJ-120/70 120 0.0741.47×610-5-7 LGJ-120/25 165 0.0791.60×610-6-9 LGJ-95/55 166 0.0921.80×610-7-8 LGJ-240/30 92 0.0471.78×610-8-9 LGJ-240/30 122 0.0471.78×610-说明：线路零序电抗为正序电抗3倍。

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的电力系统潮流计算指导教师：职称:职称：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的电⼒系统潮流计算_毕业设计论⽂基于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 潮流计算⽬的电⼒系统潮流计算是研究电⼒系统稳态运⾏情况的⼀种基本电⽓计算。

湖北民族学院信息工程学院课程设计报告书题目: 电力系统潮流计算专业：电气工程及其自动化班级：学号：学生姓名：指导教师：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的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：该程序以表格形式输出结果，节点输出包括电压幅值和相角，发电机和负荷的有功和无功功率，以及并联电容器或电抗器的有功和无功功率。

课程设计报告学生姓名:学号：学院:班级:题目:指导教师：职称: 指导教师：职称: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 潮流计算的作用根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线的电压、各支路流过的功率、整个系统的功率损耗，用以研究系统规划和运行中提出的各种问题。

对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。

东北电力大学课程设计改革试用任务书:电力系统潮流计算课程设计任务书设计名称：电力系统潮流计算课程设计设计性质：理论计算，计算机仿真与验证计划学时：两周一、设计目的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.采用五级评分制：优、良、中、及格、不及格五个等级。

基于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的潮流计算起迄日期：月日～月日课程设计地点：电气工程系工程实验室指导教师：系主任：课程设计任务书课程设计任务书1 引言潮流计算是电力系统最基本最常用的计算。

根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值和相角），各支路流过的功率，整个系统的功率损耗。

潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。

因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。

潮流计算在数学上是多元非线性方程组的求解问题，牛顿—拉夫逊Newton-Raphson法是数学上解非线性方程组的有效方法，有较好的收敛性。

运用电子计算机计算一般要完成以下几个步骤：建立数学模型，确定解算方法，制订计算流程，编制计算程序。

1.1 设计任务与要求利用MATLAB/SIMULINK仿真软件，搭建多电压等级电力网络潮流计算模型。

掌握各种潮流计算方法，掌握MATLAB/SPS工具箱使用方法及其在潮流计算中的应用。

MATLAB作为先进的仿真软件，在电力系统潮流计算方面，尤其是矩阵运算上具有突出的优秀性。

本设计要求利用MATLAB编写程序，输入数据并完成潮流分布计算。

1.2实用价值与理论意义在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。

要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。

为了保证电力系统的稳定运行，要进行潮流调节。

随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。

电力系统潮流计算是电力系统分析计算中最基本的内容，也是电力系统运行及设计中必不可少的工具。

根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线电压的幅值及相角、各元件中流过的功率、整个系统的功率损耗等。