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

合集下载

基于MATLAB的电力系统复杂潮流分析课程设计

基于MATLAB的电力系统复杂潮流分析课程设计

内蒙古科技大学本科生课程设计说明书题目:基于MATLAB的电力系统复杂潮流分析学生姓名:XXX学号:专业:电气工程及其自动化班级:电气2011—1班指导教师:XXX目录摘要--------------------------------------------------------------------------------1Abstract------------------------------------------------------------------------- 2内蒙古科技大学课程设计任务书-------------------------------------- 31 概述--------------------------------------------------------------------------51.1待设计电气设备系统图-----------------------------------------51.2电力系统概述------------------------------------------------------51.3 潮流计算的概述--------------------------------------------------61.4电力系统稳态分析潮流计算总结----------------------71.5 MATLAB的相关知识与学习----------------------------------82程序设计-------------------------------------------------------------------142.1设计思路----------------------------------------------------------142.2设计流程思路----------------------------------------------------152.3设计程序----------------------------------------------------------163 设计结果----------------------------------------------------------------204 总结-----------------------------------------------------------------------38参考文献---------------------------------------------------------------------39基于MATLAB的电力系统复杂潮流分析摘要电力系统潮流计算是电网分析的基础应用,是对复杂电力系统正常和故障条件下稳态运行状态的计算。

matalab电力系统课程设计

matalab电力系统课程设计

matalab电力系统课程设计一、课程目标知识目标:1. 学生能够理解电力系统的基础理论,掌握Matlab在电力系统中的应用。

2. 学生能够运用Matlab软件进行电力系统的基本分析和计算,如潮流计算、短路计算等。

3. 学生能够掌握电力系统中常用的Matlab函数和工具箱。

技能目标:1. 学生能够熟练操作Matlab软件,进行电力系统的模拟和仿真。

2. 学生能够利用Matlab解决电力系统中的实际问题,如优化、稳定分析等。

3. 学生能够通过Matlab编程实现电力系统的自动化分析,提高计算效率。

情感态度价值观目标:1. 学生对电力系统和Matlab产生浓厚的兴趣,培养主动探索和学习的积极性。

2. 学生在团队协作中培养良好的沟通能力和合作精神,增强解决问题的自信心。

3. 学生能够认识到电力系统对社会发展的重要性,增强对电力行业的责任感。

课程性质分析:本课程为电力系统专业课程,结合Matlab软件进行教学,旨在培养学生的实际操作能力和解决问题的能力。

学生特点分析:学生已具备一定的电力系统基础知识,对Matlab软件有一定了解,但实际操作能力较弱。

教学要求:1. 理论与实践相结合,注重培养学生的动手能力。

2. 采用案例教学,引导学生运用Matlab解决实际问题。

3. 强化团队合作,提高学生的沟通能力和协作精神。

二、教学内容1. 电力系统基础理论回顾:包括电力系统的基本概念、元件和模型,为后续Matlab应用打下基础。

- 教材章节:第一章 电力系统基本概念与元件- 内容:电力系统的组成、基本运行原理,主要电气设备的工作原理。

2. Matlab在电力系统中的应用介绍:介绍Matlab软件在电力系统分析中的应用,如潮流计算、短路计算等。

- 教材章节:第二章 Matlab在电力系统中的应用- 内容:Matlab软件的安装与配置,电力系统工具箱的介绍。

3. 潮流计算及稳定性分析:利用Matlab进行电力系统的潮流计算和稳定性分析。

用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的电力系统潮流计算课程设计课程设计报告学生姓名:学号:学院: 电气工程学院班级:题目: 基于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的电力系统潮流计算

《电力系统分析》课程设计基于Matlab的电力系统潮流计算专业电气工程及其自动化班级 1402班姓名郭娟学号 14039208同组成员孙帆郭娟陈朝辉苏新波吕应发指导教师楚冰清完成时间 2016 .12 .30目录1概述 (1)1.1设计目的 (1)1.2设计题目 (1)1.3设计内容 (1)2电力系统潮流计算概述 (2)2.1电力系统简介 (2)2.2潮流计算简介 (2)2.3潮流计算意义及其发展 (2)3 Matlab概述 (3)3.1 Matlab简介 (3)3.2 Matlab应用 (4)4 潮流计算过程 (4)4.1系统图及参数 (5)4.1.1 系统图 (5)4.1.2 各元件参数 (5)4.2电网潮流计算思路 (5)4.3潮流计算过程 (5)4.3.1各元件参数计算 (5)4.3.2绘制等效电路图 (7)4.3.3功率分布计算 (7)5 Matlab程序及运行结果 (10)5.1Matlab程序 (10)5.2运行结果 (20)6设计心得与体会 (23)参考文献 (25)1 概述1.1 设计目的在如今的社会,电力已经成为人们必不可少的需求,而建立结构合理的大型电力系统不仅便于电能生产与消费的集中管理、统一调度和分配,减少总装机容量节省动力设施投资,且有利于地区能源资源的合理开发利用,更大限度地满足地区国民经济日益增长的用电需要。

电力系统建设往往是国家及地区国民经济发展规划的重要组成部分。

电力系统的出现,使高效、无污染、使用方便、易于调控的电能得到广泛应用,推动了社会生产各个领域的变化,开创了电力时代,发生了第二次技术革命。

电力系统的规模和技术水准已成为一个国家经济发展水平的标志之一。

电力系统稳态分析包括潮流计算(或潮流分析)和静态安全分析。

潮流计算针对电力革统各正常运行方式,而静态安全分析则要研究各种运行方式下个别系统元件退出运行后系统的状况。

其目的是校验系统是否能安全运行,即是否有过负荷的元件或电压过低的母线等。

matlab电力系统分析课程设计--基于Matlab计算程序的电力系统运行分析{修}

matlab电力系统分析课程设计--基于Matlab计算程序的电力系统运行分析{修}

课程设计课程名称:电力系统分析设计题目:基于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倍。

基于MATLAB的电力系统潮流计算_毕业设计论文

基于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电力系统潮流计算课程设计

用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 任务提出与方案论证潮流计算是在给定电力系统网络结构、参数和决定系统运行状态的边界条件的情况下确定系统稳态运行状态的一种基本方法,是电力系统规划和运营中不可缺少的一个重要组成部分。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 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 母线 电厂一 电厂二
度的电抗为 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 进行潮流分析绘制的系统图,以及潮流分布图。

以下是基于原图,将线路由轻载改为重载后的课程设计:。

相关文档
最新文档