电力系统潮流分析报告

合集下载

电力系统分析实验报告

电力系统分析实验报告

电力系统分析实验报告电力系统分析实验报告引言:电力系统是现代社会不可或缺的基础设施,它为我们的生活提供了稳定的电力供应。

为了确保电力系统的可靠性和安全性,对电力系统进行分析是非常重要的。

本实验旨在通过对电力系统的分析,探讨电力系统的性能和效能,以及可能存在的问题和改进措施。

一、电力系统的基本原理电力系统由发电厂、输电网和配电网组成。

发电厂负责将化学能、机械能等转化为电能,输电网将发电厂产生的电能输送到各个地区,配电网将电能供应给终端用户。

电力系统的基本原理是通过电压和电流的传输,实现电能的转换和分配。

二、电力系统的分析方法1. 潮流计算潮流计算是电力系统分析中最基本的方法之一。

通过潮流计算,可以确定电力系统中各节点的电压和电流分布情况,从而评估系统的稳定性和负载能力。

潮流计算需要考虑各个节点的功率平衡和电压平衡,以及各个元件的参数和状态。

2. 短路分析短路分析是评估电力系统安全性的重要手段。

通过短路分析,可以确定电力系统中各个节点和支路的短路电流,从而评估设备的额定容量和保护措施的有效性。

短路分析需要考虑系统的拓扑结构、设备参数和保护装置的动作特性。

3. 阻抗分析阻抗分析是评估电力系统稳定性和负载能力的重要方法。

通过阻抗分析,可以确定电力系统中各个节点和支路的阻抗,从而评估系统的电压稳定性和电力传输能力。

阻抗分析需要考虑系统的拓扑结构、设备参数和负载特性。

三、实验结果与讨论在本实验中,我们选取了一个具体的电力系统进行分析。

通过潮流计算,我们确定了系统中各个节点的电压和电流分布情况。

通过短路分析,我们评估了系统的安全性,并确定了保护装置的动作特性。

通过阻抗分析,我们评估了系统的稳定性和负载能力。

实验结果显示,系统中存在一些节点电压偏低的问题,可能会影响设备的正常运行。

为了解决这个问题,我们建议采取增加变压器容量、调整负载分配和优化配电网结构等措施。

此外,我们还发现系统中某些支路的短路电流超过了设备的额定容量,可能导致设备的损坏和安全事故。

电力系统 运行方式及潮流分析实验结题报告(实验一)讲解

电力系统 运行方式及潮流分析实验结题报告(实验一)讲解

Beijing Jiaotong University电力系统运行方式及潮流分析实验结题报告姓名:TYP班级:电气0906学号:********指导老师:***完成日期:2012.3.20一、实验目的1、掌握电力系统主接线电路的建立方法;2、掌握辐射形网络的潮流计算方法;3、比较计算机潮流计算与手算潮流的差异;4、掌握不同运行方式下潮流分布的特点。

二、实验内容1、辐射形网络的潮流计算;2、不同运行方式下潮流分布的比较分析。

三、实验步骤1、熟悉仿真环境及主接线系统的搭建打开仿真软件,根据教程熟练软件里各项工具的使用,并最终搭建起辐射形网络主接线系统。

系统中各个模块的参数设定如下:(1)升压变压器B1根据变比=18/110,Uk%=14.3%,Pk=230KW,P0=150KW,I0/In=1%等数据设定参数,此处应将绕组2额定电压设为1.1倍的2测网络额定电压。

具体参数设定如下图:(2)线路L1,L2根据长度:80km,电阻:0.21Ω/km,电抗:0.416Ω/km,电纳:2.74×10-6S/km等数据设定参数,其中需将线路额定电压设为110kV,以保证仿真成功。

同时要通过电纳的数值换算出正序电容的数值,具体公式为b=2πf N C,算出正序电容为8.721690881nF/km。

具体参数设定如下图:(3)降压变压器B2,B3根据Un=15MVA,变比=110/11 KV,Uk%=10.5%,Pk=128KW,P0=40.5KW,I0/In=3.5%等数据设定参数。

其中2测网络额定电压须设为10kV,绕组2额定电压设为11kV,因为国家电网规定的线路电压标准中只有10kV,而没有11kV。

具体参数设定如下图:(4)负载F1,F2根据F1:20+j15MVA; F2:30+j12MVA等数据设定两个负载的参数。

具体参数设定如下图:(5)同步发电机G1根据G1:300+j180MVA(平衡节点)等数据设定参数。

简单电力系统的潮流分析报告

简单电力系统的潮流分析报告

P22 Q22
U
2 2
(R
jX )
(2)并联导纳上的损耗
QB1

1 2
BU12
QB 2

1 2
BU
2 2
(3)电力线路中的功率计算
S~1
U1 B
j 2
S~1 R
jQB1 SZ
jXS~2
jQB2
S~2 U 2 jB 2
电力线路阻抗支路末端流出的功率为
S~2 S~2 ( jQB2) P2 j(Q2 QB2 )
3.4 开式网的潮流分析
一、开式网络的潮流分析 电力系统在运行时,在电源电势激励之下,电 流或功率从电源通过系统各元件流入负荷,分 布于电力网各处,称为潮流分布。
潮流分析的步骤:
1、作出网络的等值电路 2、作出简化等值电路 3、逐段推算功率和电压分布
1 、已知 首端电压、首端功率, 求 末端电压、末端功率。
cos 0.8
100kW
3.5 闭式网络的潮流分析
闭式网络:网络中任何一个负荷点都可以 从两个及以上方向供电的电力网,如两端 供电网络、环网等。
一、闭式网络的潮流分析
闭式网络的潮流分析分为两步: 1、初步潮流分布计算——忽略各段上的功率
损耗求近似功率分布,找出无功功率分点。 2、最终潮流分布——用之前得到的近似的功
jB 2
3 、已知 首端电压、末端功率, 求 末端电压、首端功率。
从已知功率端开始,先假设末端及供电支 路各点的电压为额定电压
S~1
U1 jB 2
S~1 R
jQB1
~ SZ
U
jXS~2
jQB2
S~2
U 2

最新电力系统分析潮流计算实验报告

最新电力系统分析潮流计算实验报告

电力系统分析潮流计算实验报告姓名:XXXXXX 学号:XXXXXXXXXX 班级:XXXXXXXX一、实验目的掌握潮流计算计算机算法的方法,熟悉MATLAB的程序调试方法。

二、实验准备根据课程内容,熟悉MATLAB软件的使用方法,自行学习MATLAB程序的基础语法,并根据所学知识编写潮流计算牛顿拉夫逊法(或PQ分解法) 的计算程序,用相应的算例在MATLAB上进行计算、调试和验证。

三、实验要求每人一组,在实验课时内,调试和修改运行程序,用算例计算输出潮流结果。

四、程序流程五、实验程序%本程序的功能是用牛拉法进行潮流计算%原理介绍详见鞠平著《电气工程》%默认数据为鞠平著《电气工程》例8.4所示数据%B1是支路参数矩阵%第一列和第二列是节点编号。

节点编号由小到大编写%对于含有变压器的支路,第一列为低压侧节点编号,第二列为高压侧节点编号%第三列为支路的串列阻抗参数,含变压器支路此值为变压器短路电抗%第四列为支路的对地导纳参数,含变压器支路此值不代入计算%第五烈为含变压器支路的变压器的变比,变压器非标准电压比%第六列为变压器是否是否含有变压器的参数,其中“1”为含有变压器,“0”为不含有变压器%B2为节点参数矩阵%第一列为节点注入发电功率参数%第二列为节点负荷功率参数%第三列为节点电压参数%第四列%第五列%第六列为节点类型参数,“1”为平衡节点,“2”为PQ节点,“3”为PV节点参数%X为节点号和对地参数矩阵%第一列为节点编号%第二列为节点对地参数%默认算例% n=4;% n1=4;% isb=4;% pr=0.00001;% B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0];% B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1];% X=[1 0;2 0.05i;3 0;4 0];clear;clc;num=input('是否采用默认数据?(1-默认数据;2-手动输入)');if num==1n=4;n1=4;isb=4;pr=0.00001;B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0];B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1];X=[1 0;2 0.05i;3 0;4 0];elsen=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=input('请输入平衡节点号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入支路参数:B1=');B2=input('请输入节点参数:B2=');X=input('节点号和对地参数:X=');endTimes=1; %迭代次数%创建节点导纳矩阵Y=zeros(n);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)-B1(i,5)/B1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+B1(i,5)/B1(i,3)+(1-B1(i,5))/B1(i,3);Y(q,q)=Y(q,q)+B1(i,5)/B1(i,3)+(B1(i,5)*(B1(i,5)-1))/B1(i,3);endendfor i=1:n1Y(i,i)=Y(i,i)+X(i,2); %计及补偿电容电纳enddisp('导纳矩阵为:');disp(Y); %显示导纳矩阵%初始化OrgS、DetaSOrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1);%创建OrgS,用于存储初始功率参数h=0;j=0;for i=1:n %对PQ节点的处理if i~=isb&B2(i,6)==2 %不是平衡点&是PQ点h=h+1;for j=1:n%公式8-74%Pi=ei*(Gij*ej-Bij*fj)+fi*(Gij*fj+Bij*ej)%Qi=fi*(Gij*ej-Bij*fj)-ei*(Gij*fj+Bij*ej)OrgS(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))*imag(B2(j ,3))+imag(Y(i,j))*real(B2(j,3)));endendendfor i=1:n %对PV节点的处理,注意这时不可再将h初始化为0if i~=isb&B2(i,6)==3 %不是平衡点&是PV点h=h+1;for j=1:n%公式8-75-a%Pi=ei*(Gij*ej-Bij*fj)+fi*(Gij*fj+Bij*ej)%Qi=fi*(Gij*ej-Bij*fj)-ei*(Gij*fj+Bij*ej)OrgS(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 ))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendend%创建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); %delPiDetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1); %delQiendendt=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,1))-OrgS(2*h-1,1); %delPiDetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2; %delUi endend% DetaS%创建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));%conj求共轭endend%创建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;endendendendenddisp('初始雅可比矩阵为:');disp(Jacbi);%求解修正方程,获取节点电压的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS; %inv矩阵求逆% 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);endend% B2%开始循环**********************************************************************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 ))*imag(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 ))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendend% OrgS%创建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,1)=real(B2(i,1))-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;endend% DetaS%创建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));endend% I%创建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;endendendend% JacbiDetaU=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);endend% B2Times=Times+1; %迭代次数加1enddisp('迭代次数为:');disp(Times);disp('收敛时电压修正量为::');disp(DetaU);for k=1:nE(k)=B2(k,3);e(k)=real(E(k));f(k)=imag(E(k));V(k)=sqrt(e(k)^2+f(k)^2);sida(k)=atan(f(k)./e(k))*180./pi;end%=============== 计算各输出量=========================== disp('各节点的实际电压标幺值E为(节点号从小到大排列):'); disp(E); %显示各节点的实际电压标幺值E用复数表示disp('-----------------------------------------------------')disp('各节点的电压大小V为(节点号从小到大排列):');disp(V); %显示各节点的电压大小V的模值disp('-----------------------------------------------------');disp('各节点的电压相角sida为(节点号从小到大排列):');disp(sida); %显示各节点的电压相for p=1:nfor q=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q)); %计算各节点的注入电流的共轭值endS(p)=E(p)*C(p); %计算各节点的功率S = 电压X 注入电流的共轭值enddisp('各节点的功率S为(节点号从小到大排列):');disp(S); %显示各节点的注入功率Sline=zeros(n1,5);disp('-----------------------------------------------------');disp('各条支路的首端功率Si为(顺序同您输入B1时一致):');for i=1:n1p=B1(i,1);q=B1(i,2);Sline(i,1)=B1(i,1);Sline(i,2)=B1(i,2);if B1(i,6)==0Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));Siz(i)=Si(p,q);elseSi(p,q)=E(p)*(conj(E(p))*((1-B1(i,5))/B1(i,3))+(conj(E(p))-conj(E(q)))*(B1(i,5)/B1(i,3)));Siz(i)=Si(p,q);endSSi(p,q)=Si(p,q);Sline(i,3)=Siz(i);ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];disp(ZF);enddisp('-----------------------------------------------------');disp('各条支路的末端功率Sj为(顺序同您输入B1时一致):');for i=1:n1p=B1(i,1);q=B1(i,2);if B1(i,6)==0Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));Sjy(i)=Sj(q,p);elseSj(q,p)=E(q)*(conj(E(q))*((B1(i,5)*(B1(i,5)-1))/B1(i,3))+(conj(E(q))-conj(E(p)))*(B1(i,5)/B1(i,3)));Sjy(i)=Sj(q,p);endSSj(q,p)=Sj(q,p);Sline(i,4)=Sjy(i);ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];disp(ZF);enddisp('-----------------------------------------------------');disp('各条支路的功率损耗DS为(顺序同您输入B1时一致):');for i=1:n1p=B1(i,1);q=B1(i,2);DS(i)=Si(p,q)+Sj(q,p);DDS(i)=DS(i);Sline(i,5)=DS(i);ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];disp(ZF);enddisp('-----------------------------------------------------');disp('各支路首端编号末端编号首端功率末端功率线路损耗');disp(Sline);六、运行结果及其分析是否采用默认数据?(1-默认数据;2-手动输入)1导纳矩阵为:2.9056 -11.5015i 0.0000 + 5.3173i -1.6606 +3.1617i -1.2450 + 2.3710i0.0000 + 5.3173i 0.0000 - 4.6633i 0.0000 + 0.0000i 0.0000 + 0.0000i-1.6606 + 3.1617i 0.0000 + 0.0000i 2.4904 - 4.7039i -0.8298 + 1.5809i-1.2450 + 2.3710i 0.0000 + 0.0000i -0.8298 + 1.5809i 2.0749 - 3.9089i初始雅可比矩阵为:11.1267 2.7603 -5.3173 0 -3.1617 -1.6606-3.0509 11.8762 0 -5.3173 1.6606 -3.1617-5.3173 0 5.3173 0 0 00 -5.3173 0 4.0092 0 0-3.3198 -1.7436 0 0 4.8217 2.69800 0 0 0 0 2.1000迭代次数为:4收敛时电压修正量为::1.0e-05 *0.0349-0.2445-0.0101-0.5713-0.0931-0.0073各节点的实际电压标幺值E为(节点号从小到大排列):0.9673 - 0.0655i 1.0252 - 0.1666i 1.0495 - 0.0337i 1.0500 + 0.0000i -----------------------------------------------------各节点的电压大小V为(节点号从小到大排列):0.9695 1.0387 1.0500 1.0500-----------------------------------------------------各节点的电压相角sida为(节点号从小到大排列):-3.8734 -9.2315 -1.8419 0各节点的功率S为(节点号从小到大排列):-0.0000 + 0.0000i -0.5000 - 0.3000i 0.2000 + 0.1969i 0.3277 + 0.0443i -----------------------------------------------------各条支路的首端功率Si为(顺序同您输入B1时一致):S(1,2)=-0.5-0.30713iS(1,3)=-0.24266-0.197iS(1,4)=-0.25734-0.11013iS(3,4)=-0.055551+0.0017528i-----------------------------------------------------各条支路的末端功率Sj为(顺序同您输入B1时一致):S(2,1)=0.5+0.24606iS(3,1)=0.25555+0.1952iS(4,1)=0.2712+0.1014iS(4,3)=0.056496-0.057061i-----------------------------------------------------各条支路的功率损耗DS为(顺序同您输入B1时一致):DS(1,2)=0-0.06107iDS(1,3)=0.012892-0.0018014iDS(1,4)=0.013863-0.0087295iDS(3,4)=0.00094545-0.055308i-----------------------------------------------------各支路首端编号末端编号首端功率末端功率线路损耗1.0000 + 0.0000i2.0000 + 0.0000i -0.5000 - 0.3071i 0.5000 + 0.2461i 0.0000 - 0.0611i 1.0000 + 0.0000i3.0000 + 0.0000i -0.2427 - 0.1970i 0.2556 + 0.1952i 0.0129 - 0.0018i 1.0000 + 0.0000i4.0000 + 0.0000i -0.2573 - 0.1101i 0.2712 + 0.1014i 0.0139 - 0.0087i3.0000 + 0.0000i4.0000 + 0.0000i -0.0556 + 0.0018i 0.0565 - 0.0571i 0.0009 - 0.0553i七、实验体会及感悟通过这次实验,首先让我对matlab软件有了初步的了解,对它强大的矩阵运算能力有了更深的体会,同时掌握了设置断点和断点调试的一般方法,结合课本上的程序流程图和参考资料上的例子单步跟踪调试,再一次的熟悉了牛顿拉夫逊法潮流计算的一般方法和步骤,对计算机计算潮流计算有了更进一步的认识,在学习潮流计算时,虽然依次学习了节点导纳矩阵,功率方程、雅可比矩阵,但不能将它们联系起来,更不知道其中的原委,通过程序的编写,知道了其中的联系,也知道了每个方程、矩阵在计算中的作用。

电力系统分析实验报告

电力系统分析实验报告

一、实验目的1. 了解电力系统的基本组成和运行原理;2. 掌握电力系统潮流计算的方法和步骤;3. 熟悉电力系统故障计算的方法和步骤;4. 培养分析电力系统问题的能力。

二、实验原理1. 电力系统潮流计算:通过求解电力系统中的潮流方程,得到系统中各节点的电压、电流、功率等参数,从而分析电力系统的运行状态。

2. 电力系统故障计算:通过求解电力系统中的故障方程,得到故障点附近的电压、电流、功率等参数,从而分析电力系统故障的影响。

三、实验仪器与设备1. 电力系统分析软件:如PSCAD/EMTDC、MATLAB等;2. 电力系统仿真设备:如电力系统仿真机、计算机等;3. 电力系统相关教材和资料。

四、实验步骤1. 建立电力系统模型:根据实验要求,利用电力系统分析软件建立电力系统模型,包括发电机、变压器、线路、负荷等元件。

2. 潮流计算:(1)设置初始条件:根据实验要求,设置电力系统运行状态,如电压、功率等;(2)求解潮流方程:利用电力系统分析软件求解潮流方程,得到系统中各节点的电压、电流、功率等参数;(3)分析潮流计算结果:根据计算结果,分析电力系统的运行状态,如电压分布、潮流分布等。

3. 故障计算:(1)设置故障条件:根据实验要求,设置电力系统故障,如短路、断路等;(2)求解故障方程:利用电力系统分析软件求解故障方程,得到故障点附近的电压、电流、功率等参数;(3)分析故障计算结果:根据计算结果,分析电力系统故障的影响,如电压波动、潮流变化等。

五、实验结果与分析1. 潮流计算结果分析:(1)电压分布:根据潮流计算结果,分析系统中各节点的电压分布情况,判断电压是否满足运行要求;(2)潮流分布:根据潮流计算结果,分析系统中各线路的潮流分布情况,判断潮流是否合理。

2. 故障计算结果分析:(1)故障点电压:根据故障计算结果,分析故障点附近的电压变化情况,判断电压是否满足运行要求;(2)故障点电流:根据故障计算结果,分析故障点附近的电流变化情况,判断电流是否过大;(3)故障点功率:根据故障计算结果,分析故障点附近的功率变化情况,判断功率是否过大。

电力系统潮流分析报告

电力系统潮流分析报告
(6)计算各节点电压的新值
(7)运用各节点电压的新值自第三步开始进入下一次迭代。
(8)计算平衡节点功率。
第三节 计算机算法程序
第三章
第一节 PQ分解法潮流计算时的修正方程式
PQ分解法潮流计算派生于以极坐标表示时的牛顿拉夫逊法,其修正方程式是牛顿拉夫逊修正方程式的简化。为说明这一简化,将牛孙拉弗逊极坐标表示时的雅可比矩阵重新排列,由于高压输电线路的电阻远远小于电抗,以致各节点电压相位角的改变主要影响各元件中的有功功率潮流从而影响各节点的注入有功功率;各节点电压大小的改变主要影响各元件中无功功率潮流从而影响各节点的注入无功功率,可将修正方程简化为
(3)网络中有n—m个PV节点,编号m+1,m+2,…n。
至此,就可以建立修正方程,如下所示
(2-5)
式中雅可比矩阵各个元素分别为
(2—6)
第二节 牛顿拉夫逊法潮流计算的基本步骤
(1)形成节点导纳矩阵
(2)设置各节点电压的初始值
(3)代入初始值,求修正方程式中的不平衡量
(4)求取雅可比矩阵中各元素
(5)解修正方程式,求各节点的修正量。
(7)解修正方程式(4—3b),求各节点相位角的修正量
(8)求取各节点电压大小的新值
(9)运用各节点电压的新值自第三部开始下一次迭代
(10)计算平衡节点功率
第三节 PQ分解法的MATLAB实现(见附件)
第四章
与牛顿拉弗逊相比,PQ分解法有如下特点:
(1)以一个(n-1)阶和一个(m-1)阶系数矩阵 、 代替原有的(n+m—2)阶矩阵J,提高了计算速度,降低了对内存的要求。
(2)以迭代过程中保持不变的系数矩阵 、 替代不对称的系数矩阵J,提高了计算速度。

电力系统运行方式及潮流分析实验报告

电力系统运行方式及潮流分析实验报告

电力系统第一次实验报告——电力系统运行方式及潮流分析实验实验1 电力系统运行方式及潮流分析实验一、实验目的1、掌握电力系统主接线电路的建立方法2、掌握辐射形网络的潮流计算方法;3、比较计算机潮流计算与手算潮流的差异;4、掌握不同运行方式下潮流分布的特点。

二、实验内容1、辐射形网络的潮流计算;2、不同运行方式下潮流分布的比较分析三、实验方法和步骤1.辐射形网络主接线系统的建立输入参数(系统图如下):G1:300+j180MV A(平衡节点)变压器B1:Sn=360MV A,变比=18/121,Uk%=14.3%,Pk=230KW,P0=150KW,I0/In=1%;变压器B2、B3:Sn=15MV A,变比=110/11 KV,Uk%=10.5%,Pk=128KW,P0=40.5KW,I0/In=3.5%;负荷F1:20+j15MV A;负荷F2:28+j10MV A;线路L1、L2:长度:80km,电阻:0.21Ω/km,电抗:0.416Ω/km,电纳:2.74×10-6S/km。

辐射形网络主接线图(1)在DDRTS中绘出辐射形网络主接线图如下所示:(2)设置各项设备参数:G1:300+j180MV A(平衡节点)变压器B1:Sn=360MV A,变比=18/121,Uk%=14.3%,Pk=230KW,P0=150KW,I0/In=1%;变压器B2、B3:Sn=15MV A,变比=110/11 KV,Uk%=10.5%,Pk=128KW,P0=40.5KW,I0/In=3.5%;负荷F1:20+j15MV A;负荷F2:28+j10MV A;线路L1、L2:长度:80km,电阻:0.21Ω/km,电抗:0.416Ω/km,电纳:2.74×10-6S/km。

2.辐射形网络的潮流计算(1)调节发电机输出电压,使母线A的电压为115KV,运行DDRTS进行系统潮流计算,在监控图页上观察计算结果(2)手算潮流:(3)计算比较误差分析通过比较可以看出,手算结果与计算机仿真结果相差不大。

电力系统运行方式及潮流分析实验报告记录

电力系统运行方式及潮流分析实验报告记录

电力系统运行方式及潮流分析实验报告记录实验目的:1.了解电力系统的运行方式2.掌握潮流分析的基本原理及方法3.学会利用软件进行电力系统潮流计算实验仪器和设备:1.个人计算机2.电力系统分析软件:PSASP实验内容:1.电力系统潮流分析的基本原理及方法学习根据电力系统的基本控制方程式,通过潮流计算方法,求出电力系统各节点的电压、功率及电流等参数,以及各支路上的潮流大小和方向等相关参数。

2.利用软件进行电力系统潮流计算利用电力系统分析软件PSASP,建立某电力系统的模型,输入各节点的电源电压、负荷等基本信息,进行潮流计算,得出各节点的电压、功率及电流等参数,以及各支路上的潮流大小和方向等相关参数。

3.电力系统运行方式的了解电力系统是由发电厂、变电站、输电线路、配电站及用户组成的一种大型电力供应体系。

其运行方式一般分为功率控制和电压控制两种方式,其中功率控制是指在保持电压稳定的情况下,调整电网内发电和负荷之间的功率平衡,即通过调节发电量、负荷等手段,使得电网内的功率平衡得以维持;而电压控制则是指在保持电网功率平衡的情况下,调节电网内的电压稳定状态,即通过调节发电机电压、变压器变比、补偿电抗器等手段,使得电网内各节点的电压处于稳定状态。

实验结果分析:通过潮流分析计算,得出某电力系统各节点的电压、功率及电流等参数,和各支路上的潮流大小和方向等相关参数。

同时,通过学习电力系统运行方式,了解到电力系统在运行过程中,既要保持功率平衡,又要保持电压稳定,以确保电力的供应可靠性和稳定性。

实验结论:通过本次实验,我掌握了潮流分析的基本原理及方法,学会了利用电力系统分析软件进行电力系统潮流计算,并深入了解了电力系统的运行方式和特点,从而提高了对电力系统的理解和认识。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

潮流计算的意义(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。

(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。

(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。

(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。

总结为在电力系统运行方式和规划方案的研究中,都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。

同时,为了实时监控电力系统的运行状态,也需要进行大量而快速的潮流计算。

因此,潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。

在系统规划设计和安排系统的运行方式时,采用离线潮流计算;在电力系统运行状态的实时监控中,则采用在线潮流计算。

编辑本段潮流计算的发展史利用电子计算机进行潮流计算从20世纪50年代中期就已经开始。

此后,潮流计算曾采用了各种不同的方法,这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的。

对潮流计算的要求可以归纳为下面几点:(1)算法的可靠性或收敛性(2)计算速度和存占用量(3)计算的方便性和灵活性电力系统潮流计算属于稳态分析畴,不涉及系统元件的动态特性和过渡过程。

因此其数学模型不包含微分方程,是一组高阶非线性方程。

非线性代数方程组的解法离不开迭代,因此,潮流计算方法首先要求它是能可靠的收敛,并给出正确答案。

随着电力系统规模的不断扩大,潮流问题的方程式阶数越来越高,目前已达到几千阶甚至上万阶,对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的。

这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法。

在用数字计算机求解电力系统潮流问题的开始阶段,人们普遍采用以节点导纳矩阵为基础的高斯-赛德尔迭代法(一下简称导纳法)。

这个方法的原理比较简单,要求的数字计算机的存量也比较小,适应当时的电子数字计算机制作水平和电力系统理论水平,于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法(以下简称阻抗法)。

20世纪60年代初,数字计算机已经发展到第二代,计算机的存和计算速度发生了很大的飞跃,从而为阻抗法的采用创造了条件。

阻抗矩阵是满矩阵,阻抗法要求计算机储存表征系统接线和参数的阻抗矩阵。

这就需要较大的存量。

而且阻抗法每迭代一次都要求顺次取阻抗矩阵中的每一个元素进行计算,因此,每次迭代的计算量很大。

阻抗法改善了电力系统潮流计算问题的收敛性,解决了导纳法无法解决的一些系统的潮流计算,在当时获得了广泛的应用,曾为我国电力系统设计、运行和研究作出了很大的贡献。

但是,阻抗法的主要缺点就是占用计算机的存很大,每次迭代的计算量很大。

当系统不断扩大时,这些缺点就更加突出。

为了克服阻抗法在存和速度方面的缺点,后来发展了以阻抗矩阵为基础的分块阻抗法。

这个方法把一个大系统分割为几个小的地区系统,在计算机只需存储各个地区系统的阻抗矩阵及它们之间的联络线的阻抗,这样不仅大幅度的节省了存容量,同时也提高了节省速度。

克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法(以下简称牛顿法)。

牛顿法是数学中求解非线性方程式的典型方法,有较好的收敛性。

解决电力系统潮流计算问题是以导纳矩阵为基础的,因此,只要在迭代过程中尽可能保持方程式系数矩阵的稀疏性,就可以大大提高牛顿潮流程序的计算效率。

自从20世纪60年代中期采用了最佳顺序消去法以后,牛顿法在收敛性、存要求、计算速度方面都超过了阻抗法,成为直到目前仍被广泛采用的方法。

在牛顿法的基础上,根据电力系统的特点,抓住主要矛盾,对纯数学的牛顿法进行了改造,得到了P-Q分解法。

P-Q 分解法在计算速度方面有显著的提高,迅速得到了推广。

牛顿法的特点是将非线性方程线性化。

20世纪70年代后期,有人提出采用更精确的模型,即将泰勒级数的高阶项也包括进来,希望以此提高算法的性能,这便产生了保留非线性的潮流算法。

另外,为了解决病态潮流计算,出现了将潮流计算表示为一个无约束非线性规划问题的模型,即非线性规划潮流算法。

近20多年来,潮流算法的研究仍然非常活跃,但是大多数研究都是围绕改进牛顿法和P-Q分解法进行的。

此外,随着人工智能理论的发展,遗传算法、人工神经网络、模糊算法也逐渐被引入潮流计算。

但是,到目前为止这些新的模型和算法还不能取代牛顿法和P-Q分解法的地位。

由于电力系统规模的不断扩大,对计算速度的要求不断提高,计算机的并行计算技术也将在潮流计算中得到广泛的应用,成为重要的研究领域。

编辑本段潮流计算的发展趋势通过几十年的发展,潮流算法日趋成熟。

近几年,对潮流算法的研究仍然是如何改善传统的潮流算法,即高斯-塞德尔法、牛顿法和快速解耦法。

牛顿法,由于其在求解非线性潮流方程时采用的是逐次线性化的方法,为了进一步提高算法的收敛性和计算速度,人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来,于是产生了二阶潮流算法。

后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点,提出了采用直角坐标的保留非线性快速潮流算法。

对于保留非线性算法典型论文有:1.文献[保留非线性的电力系统概率潮流计算]提出了它在电力系统概率潮流计算中的应用。

该文献提出了一种新的概率潮流计算方法,它保留了潮流方程的非线性,又利用了P-Q解耦方法,因而数学模型精度较高,且保留了P-Q解耦的优点,有利于大电网的随机潮流计算,用提出的方法对一个典型的系统进行了计算,其数值用MonteCarlo随机模拟作了验证,得到了满意的结果。

2.文献[基于系统分割的保留非线性的快速P-Q解耦潮流计算法]分析研究了保留非线性的P-Q解耦快速潮流计算法。

该文献提出了一种新的状态估计算法,既保留了量测方程非线性又利用了快速P-Q分解方法,因此数学模型精度高且保留了快速P-Q分解的优点,提高了状态估计的计算精度和速度.采用系统分割方法将大系统分割为多个小系统,分别对每个小系统进行状态估计,然后对各小系统的状态估计结果进行协调,得到整个系统具有同一参考节点的状态估计结果,这样可大大提高状态估计的计算速度,有利于进行大电网的状态估计.在18节点系统上进行的数字仿真实验验证了该方法的有效性。

岩本伸一等提出了一种保留非线性的快速潮流计算法,但用的是直角坐标系,因而没法利用P-Q解耦。

为了更有利于大电网的潮流计算,将此原理推广用于P-Q 解耦。

这样,既利用了保留非线性的快速算法,在迭代中使用常数雅可比矩阵,又保留了P-Q解耦的优点。

对于一些病态系统,应用非线性潮流计算方法往往会造成计算过程的振荡或者不收敛,从数学上讲,非线性的潮流计算方程组本来就是无解的。

这样,人们提出来了将潮流方程构造成一个函数,求此函数的最小值问题,称之为非线性规划潮流的计算方法。

优点是原理上保证了计算过程永远不会发散。

如果将数学规划原理和牛顿潮流算法有机结合一起就是最优乘子法。

另外,为了优化系统的运行,从所有以上的可行潮流解中挑选出满足一定指标要求的一个最佳方案就是最优潮流问题。

最优潮流是一种同时考虑经济性和安全性的电力网络分析优化问题。

OPF 在电力系统的安全运行、经济调度、可靠性分析、能量管理以及电力定价等方面得到了广泛的应用。

最优潮流方面的典型论文有:1.文献[电力系统最优潮流新算法的研究]以NCP 方法为基础,提出了一种新的求解最优潮流算法——投影渐近半光滑牛顿型算法。

该文献以NCP方法为基础,提出了一种新的求解OPF算法——投影渐近半光滑牛顿型算法。

针对电力系统的特点,本文的研究工作如下: 1.建立了与OPF问题的KKT系统等价的带界约束的半光滑方程系统。

与已有的NCP 方法相比,新的模型由于无需考虑界约束对应的对偶变量(乘子变量),降低了问题的维数,从而适用于解大规模的电力系统问题。

2.基于建立的新模型,本文提出了一类新的Newton型算法,该算法一方面保持界约束的相容性,另一方面有较好的全局与局部超线性收敛性,同时,算法结构简单,易于实现。

3.考虑到电力系统固有的弱耦合特性,受传统解耦最优潮流方法的启示,在所提出的新Newton型方法的基础上,本文又设计了一类分解方法。

新方法基于解耦——校正的策略实现算法,不仅充分利用了系统的弱耦合特性,同时保证分解算法在理论上的收敛性。

4.根据所提出的两种算法,用标准的IEEE电力测试系统进行数值实验,并与已有的其他方法进行比较。

结果显示新算法具有良好的收敛性和计算效果,在电力系统的规划与运行方面将有广阔的应用前景。

2.文献[基于可信域点法的最优潮流问题研究]介绍了OPF点法具有收敛性强、多项式时间复杂性等优点,是极具潜力的优秀算法之一。

电力系统不断发展,使得OPF算法跻身于极其困难、非凸的大规模非线性规划行列。

可信域和线性搜索方法是保证最优化算法全局收敛性能的两类技术,将点法和可信域、线性搜索方法有机结合,构造新的优化算法,是数学规划领域的研究热点。

此方面的典型文献有:1.文献[电力市场环境下基于最优潮流的输电容量充裕度研究]首先以最优潮流为工具,选取系统中的关键线路作为系统输电容量充裕度的研究对象,从电网运行的安全性、可靠性的角度系统地研究了输电线路稳定限额对输电容量充裕度的影响,指出稳定限额因子与影子价格的乘积可直接反应出稳定限额水平的经济价值,同时也可以较好的指示出系统运行相对安全、经济的稳定限额水平区间。

2.文献[电力市场环境下基于最优潮流的节点实时电价和购电份额研究]为了为配电公司最优购电模型提供价格参考依据,以发电成本最小为目标函数,考虑电力需求价格弹性的影响,建立了实时电价模型。

模型利用预测校正原对偶点法求解,以IEEE30节点系统为算例验证了模型的可行性。

3.文献[电力系统动态最优潮流的模型与算法研究]指出电力系统动态最优潮流是对调度周期的系统状态进行统一优化的有效工具,对保证电力系统安全经济运行具有重要的理论意义和现实意义。

文献结合点法和免疫遗传算法,对经典动态最优潮流问题和动态无功优化问题的算法进行了深入的研究,提出了新的算法;并建立了含电压稳定约束、含无功型离散变量,以及含机组启停变量的动态最优潮流模型,将新算法推广应用于各种新模型,拓展了动态最优潮流的研究领域。

对于一些特殊性质的潮流计算问题有直流潮流计算方法、随机潮流计算方法和三相潮流计算方法。

直流潮流计算方法,文献[基于改进布罗伊登法的交直流潮流计算]主要介绍在分析求解非线性方程组的布罗伊登法和一种改进的布罗伊登法的基础上,针对交直流混联系统,运用改进的布罗伊登法,提出了一种潮流计算的统一迭代法,设计了算法的具体实现步骤,并以一个IEEE9节点修改系统进行仿真计算,结果表明本文采用的改进布罗伊登法交直流潮流计算方法有效可行。

相关文档
最新文档