电力系统分析潮流计算最终完整版

合集下载

电力系统潮流计算完整c语言程序(含网损计算的最终版)

for(j=0;j<Bus_Num;j++)
{
ia[i]=ia[i]+gY_G[n][j]*ge[j]-gY_B[n][j]*gf[j];
ib[i]=ib[i]+gY_G[n][j]*gf[j]+gY_B[n][j]*ge[j];
}
}
for(i=0,n=1;i<Bus_Num-1;i++,n++)
{
gDelta_PQ[2*i]=gDelta_P[i];
gDelta_PQ[2*i+1]=gDelta_Q[i];
}
if((fp=fopen("C:\\Documents and Settings\\Zorro\\桌面\\1\\data\\unbalance.txt","w"))==NULL)
if(gBus[n].Type==1)
gDelta_Q[i]=gDelta_Q[i]-gf[n]*(gY_G[n][j]*ge[j]-gY_B[n][j]*gf[j])+ge[n]*(gY_G[n][j]*gf[j]+gY_B[n][j]*ge[j]);
}
}
for(i=0;i<Bus_Num-1;i++)
{
gY_G[i][j]=0.0;
gY_B[i][j]=0.0;
}
for(l=0;l<Line_Num;l++)
{
i=gLine[l].No_I-1;
j=gLine[l].No_J-1;
r=gLine[l].R;
x=gLine[l].X;

电力系统潮流计算

3.2.1 节点电压方程与节点导纳矩阵和阻抗矩阵
将节点电压法应用于电力系统潮流计算，变量为节点电压与节
点注入电流。通常以大地作为电压幅值的参考（U0 = 0），以
系统中某一指定母线的电压角度作为电压相角的参考，以支路
导纳作为电力网的参数进行计算。节点注入电流规定为流向网
络为正，流出为负。
Pmax P
表征年有功负荷曲线特点的两个指标
0
年最大负荷利用小时数 Tmax
t Tmax 8760
根据年负荷曲线，可求得全年所需电能：
8760
A 0
Pdt MWh
定义年最大负荷（最大值 Pmax）利用小时： Tmax
A Pmax
h
Tmax 越大，负荷曲线越平坦
负荷曲线为一水平线时， Tmax 达到最大值8760 (h)
2
1 ZT1
2
Zl
T2
34
3
ZT2 4
YT3
Yl /2
YT2
已知末端功率和电压，计算网上潮流分布。
1 ZT1 2 Zl
3 ZT2 4
已知始端功率和电压，计算网上潮流分布。
Y20
Y30
已知末端功率和始端电压，计算网上的潮流。
不管哪种情况，先作等值电路
3.1.3 辐射形网络的分析计算
1）已知末端功率、电压利用前面的方法，从末端逐级往上推算，直至求得各要求的量。
Pm(t)
损耗称年电能损耗，是电网运行经
济性的指标。
Pmi
1）年电能损耗的准确计算方法
已知各负荷的年有功和无功负荷曲线时，理论上可准确计算年电能损耗。
8760小时分为 n 段，第 i 时段时间为 Dti (h)，全网功率损耗为DPi (MW)，则全网年电能损耗为

电力系统分析潮流计算最终完整版

电力系统分析潮流计算最终完整版电力系统潮流计算是电力系统运行的基础，它对电力系统的稳定运行和安全运行具有重要意义。

本文将介绍电力系统潮流计算的主要内容和步骤，并阐述其在电力系统运行中的应用。

电力系统潮流计算是指对电力系统中各节点的电压和功率进行计算和分析的过程。

它主要用于确定电力系统中各个节点的电压和相应的功率，以评估电力系统的稳定性和安全性。

潮流计算的结果可以用于电力系统的规划、调度和运行等各个环节。

潮流计算的主要步骤主要包括：建立电力系统潮流模型、制定潮流计算方程、选择潮流计算方法和求解潮流计算方程。

建立电力系统潮流模型是潮流计算的第一步，它主要包括确定电力系统的拓扑结构、电气参数和发电机和负荷模型等。

通过建立电力系统的拓扑结构和电气参数，可以确定电力系统中各个节点之间的连接关系和传输条件。

发电机和负荷模型则用于描述电力系统中的发电机和负荷之间的相互作用。

制定潮流计算方程是潮流计算的第二步，它主要是根据电力系统的拓扑结构和电气参数，建立潮流计算的数学模型。

潮流计算方程主要包括功率方程、节点电压方程和变压器方程等。

功率方程用于描述发电机和负荷之间的功率平衡关系，节点电压方程用于描述电力系统中各个节点的电压平衡关系，变压器方程用于描述变压器的运行状况。

选择潮流计算方法是潮流计算的第三步，它主要是选择合适的方法来求解潮流计算方程。

常见的方法包括直接迭代法、高斯-赛德尔迭代法、牛顿-拉夫逊迭代法和快速迭代法等。

不同的方法在精度和收敛速度上有所差异，根据实际情况选择合适的方法。

求解潮流计算方程是潮流计算的最后一步，它主要是通过迭代计算，求解潮流计算方程得到电力系统各个节点的电压和功率值。

在求解过程中，需要根据实际情况设置迭代的初始值和收敛条件，以保证计算结果的准确性和稳定性。

电力系统潮流计算在电力系统运行中具有广泛的应用。

它可以用于电力系统规划，通过计算电力系统中各个节点的电压和功率，评估电力系统的输电能力和供电质量，为电力系统的扩容和优化提供指导。

(完整)电力系统潮流计算方法分析

电力系统潮流分析—基于牛拉法和保留非线性的随机潮流姓名：**＊学号:＊**1 潮流算法简介1.1 常规潮流计算常规的潮流计算是在确定的状态下.即：通过已知运行条件（比如节点功率或网络结构等)得到系统的运行状态（比如所有节点的电压值与相角、所有支路上的功率分布和损耗等）。

常规潮流算法中的一种普遍采用的方法是牛顿-拉夫逊法.当初始值和方程的精确解足够接近时，该方法可以在很短时间内收敛.下面简要介绍该方法。

1.1。

1牛顿拉夫逊方法原理对于非线性代数方程组式(1-1),在待求量x 初次的估计值(0)x 附近，用泰勒级数（忽略二阶和以上的高阶项）表示它，可获得如式(1-2)的线性化变换后的方程组，该方程组被称为修正方程组。

'()f x 是()f x 对于x 的一阶偏导数矩阵,这个矩阵便是重要的雅可比矩阵J 。

12(,,,)01,2,,i n f x x x i n ==（1-1）(0)'(0)(0)()()0f x f x x +∆=（1—2)由修正方程式可求出经过第一次迭代之后的修正量(0)x ∆，并用修正量(0)x ∆与估计值(0)x 之和，表示修正后的估计值(1)x ，表示如下（1—4）.(0)'(0)1(0)[()]()x f x f x -∆=-（1—3)(1)(0)(0)x x x =+∆（1-4)重复上述步骤.第k 次的迭代公式为： '()()()()()k k k f x x f x ∆=-（1—5）(1)()()k k k x x x +=+∆（1-6）当采用直角坐标系解决潮流方程，此时待解电压和导纳如下式:i i i ij ij ijV e jf Y G jB =+=+ （1-7）假设系统的网络中一共设有n 个节点，平衡节点的电压是已知的，平衡节点表示如下.n n n V e jf =+（1-8）除了平衡节点以外的所有2(1)n -个节点是需要求解的量。

电力系统分析潮流计算最终完整版

电力系统分析潮流计算实验报告姓名：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软件有了初步的了解，对它强大的矩阵运算能力有了更深的体会，同时掌握了设置断点和断点调试的一般方法，结合课本上的程序流程图和参考资料上的例子单步跟踪调试，再一次的熟悉了牛顿拉夫逊法潮流计算的一般方法和步骤，对计算机计算潮流计算有了更进一步的认识，在学习潮流计算时，虽然依次学习了节点导纳矩阵，功率方程、雅可比矩阵，但不能将它们联系起来，更不知道其中的原委，通过程序的编写，知道了其中的联系，也知道了每个方程、矩阵在计算中的作用。

电力系统分析潮流计算

电力系统分析潮流计算电力系统分析是对电力系统运行状态进行研究、分析和评估的一项重要工作。

其中，潮流计算是电力系统分析的一种重要方法，用于计算电力系统中各节点的电压、功率和电流等参数。

本文将详细介绍电力系统潮流计算的原理、方法和应用。

一、电力系统潮流计算的原理电力系统潮流计算是基于潮流方程的求解，潮流方程是描述电力系统各节点电压和相角之间的关系的一组非线性方程。

潮流方程的基本原理是基于电力系统的等效导纳矩阵和节点电压相位差的关系，通过潮流计算可以得到电力系统各节点的电压和功率等参数。

电力系统潮流方程的一般形式如下：\begin{align*}P_i &= \sum_{j=1}^{n}(V_iV_j(G_{ij}\cos(\theta_i-\theta_j)+B_{ij}\sin(\theta_i-\theta_j))) \\Q_i &= \sum_{j=1}^{n}(V_iV_j(G_{ij}\sin(\theta_i-\theta_j)-B_{ij}\cos(\theta_i-\theta_j)))\end{align*}其中，$n$为节点数，$P_i$和$Q_i$表示第i个节点的有功功率和无功功率。

$V_i$和$\theta_i$表示第i个节点的电压和相角。

$G_{ij}$和$B_{ij}$表示节点i和节点j之间的等效导纳。

二、电力系统潮流计算的方法电力系统潮流计算的方法主要包括直接法、迭代法和牛顿-拉夫逊法等。

1.直接法：直接法是一种适用于小规模电力系统的潮流计算方法，它通过直接求解潮流方程来计算电力系统的潮流。

直接法的计算速度快，但对系统规模有一定的限制。

2.迭代法：迭代法是一种常用的潮流计算方法，通常使用高尔顿法或牛顿法。

迭代法通过迭代求解潮流方程来计算电力系统的潮流。

迭代法相对于直接法来说，可以适用于大规模电力系统，但计算时间较长。

3.牛顿-拉夫逊法：牛顿-拉夫逊法是一种高效的潮流计算方法，它通过求解潮流方程的雅可比矩阵来进行迭代计算，可以有效地提高计算速度。

电力系统的潮流计算

V1

有效值：
V1、V2间的相位角
2019/4/26
V2＝ (V1－V1)2 (V1)2
arctg
V1
6
V1－V1
电力系统分析第十一章电力系统的潮流计算
注意：
P '' R Q '' X V2 V2 P'' X Q '' R V2 V2
≠
V1
P ''2 Q ''2 S ＝ 2 ( R jX ) S T T V 2

S ＝ S S T S 0
∝与负荷2
2019/4/26
与负荷无关，∝V2
13
电力系统分析第十一章电力系统的潮流计算
三、实际计算 1. 已知末端功率与电压，求另一端功率和电压
S ＝ P jQ S 2 2 2 LD
2019/4/26
电力系统分析第十一章电力系统的潮流计算
15
第十一章电力系统潮流计算
定义根据给定的运行条件求取给定运行条件下的节点电压和功率分布意义电力系统分析计算中最基本的一种：规划、扩建、运行方式安排
2019/4/26
电力系统分析第十一章电力系统的潮流计算
16
所需知识
V V V A G 1 2
当输电线路不长，首末两端的相角差不大时，近似有：
V1
δ O
B
AG≈AD
A
V2
D
G

V V 1 2 百分数表示： V% 100 V N
2019/4/26
电力系统分析第十一章电力系统的潮流计算

电力系统分析计算公式

电力系统分析计算公式1.电力系统潮流计算电力系统潮流计算是一种用于确定电力系统各个节点电压和功率的方法。

常用的电力系统潮流计算公式包括：- 节点功率方程：P = V * I * cos(theta) + V * U * sin(theta) - 节点电流方程：I = V * I * sin(theta) - V * U * cos(theta)其中，P为节点有功功率，V为节点电压，I为节点电流，theta为节点相角，U为无功功率系数。

2.短路电流计算短路电流计算是用于评估电力系统短路故障时电流的大小和方向的方法。

常用的短路电流计算公式包括：- 对称短路电流公式：Isc = V / Zs其中，Isc为短路电流，V为电压，Zs为短路阻抗。

3.电力系统电压稳定性计算电力系统电压稳定性计算是为了评估电力系统节点电压的稳定性。

常用的电力系统电压稳定性计算公式包括：-V/Q稳定器灵敏度公式：dV/dQ=-Ry*dQ/dP+Xy*(dQ/dQ+dV/dV)其中，V为节点电压，Q为节点无功功率，P为节点有功功率，Ry为负荷灵敏度，Xy为发电机灵敏度。

4.功率系统频率计算功率系统频率计算是为了评估电力系统频率的稳定性。

常用的功率系统频率计算公式为：- 系统频率变化率公式：df/dt = (P - Pd) / (2 * H)其中，df/dt为频率变化率，P为实际功率，Pd为负荷功率，H为系统等效惯量。

5.电力系统稳定裕度计算电力系统稳定裕度计算是为了评估电力系统在各种故障情况下的稳定性。

常用的电力系统稳定裕度计算公式包括：- 稳定裕度指标公式：S ω = (δmax - δmin) / δfc其中，Sω为稳定裕度指标，δmax为最大转子转角，δm in为最小转子转角，δfc为临界转子转角。

以上是一些常用的电力系统分析计算公式，这些公式是电力系统工程师进行电力系统设计和运行评估的重要依据。

电力系统分析计算的结果可以帮助工程师评估电力系统的稳定性，指导运维工作，并制定相应的措施以确保电力系统的安全、可靠和高效运行。

电力系统分析第7章电力系统的潮流计算

Sa 1 Sb 2
( Z12 Z b 2 )S1 Z b 2 S2 Za 2 Z12 Z b 2 Za1 S1 ( Za1 Z12 )S2 Za 2 Z12 Z b 2
* * * * * * * * *
*
*
*

(U a U b )U N Z a 2 Z12 Z b 2
4.运算功率发电厂的等值电源功率减去发电厂输出母线上所有相连线路的充电功率的一半称为发电厂的运算电源功率，简称运算功率。
电力系统分析 7.3.2开式网络的潮流计算方法
第7章电力系统的潮流计算
计算步骤：
1、计算电力网各元件参数，作电力网等值电路。 2、计算变电所的运算负荷和发电厂的运算功率，并将它们接在相应的节点上，从而组成了只包括运算负荷和运算功率及网络参数的等值网络。 3、如果已知电源电压和末端负荷，由末端向首端逐段计算功率损耗，这种情况由于各点电压未知，可用电网额定电压代替实际电压，求取电力网的功率分布。求得电源功率后，再运用已知电源电压和求得的首端功率向末端逐段求电压降落，计算出各点电压。此过程不必重新计算功率损耗，在110kv的高压电网中也可忽略电压降落的横分量。 4如果已知末端电压和负荷，从末端开始逐段交替计算电压降落和功率损耗。向电源端推算功率分布和各节点电压。如果有变压器，还应进行电压归算。
P2 X U 2 ( U 2 U 2 )
电力系统分析结论
第7章电力系统的潮流计算
电压降落的纵分量取决于所输送的无功功率的大小；电压降落的横分量主要取决于所输送的有功功率的大小。纵分量主要影响电压的大小，横分量主要影响电压的相角。（2）电压损耗：电力线路首末端或电力网任意两节点间电压的代数差。电压损耗近似等于电压的纵分量大小

电力系统分析电力系统潮流计算

对于电力系统中的每个节点，要确定其运行状态，需要有四个变量：有功注入无功注入电压模值电压相角 n个节点总共有4n个运行变量要确定。再观察式 (1-6)或式(1-7)，总共包括n个复数方程式，如果将实部与虚部分开，则形成2n个实数方程式，由此仅可以解得2n个未知运行变量。为此在计算潮流以前，必须将另外2n个变量作为已知量而预先给以指定。也即对每个节点，要给定其两个变量的值作为已知条件，而另两个变量作为待求量。
一些实际用于生产的潮流程序往往在上述基本潮流的框架内再加入模拟实际系统运行控制特点的自动调整计算功能，如潮流控制，分接头调整等，这部分内容将在本章第八节中予以介绍。
60年代中期，结合电力系统经济调度工作的开展，针对经典的经济调度方法的不足，开辟了一个新的研究领域，称之为最优潮流问题。这种以非线性规划作为计算模型的潮流问题能够统筹兼顾电力系统的经济性、安全性和电能质量，因而受到很大的重视，发展很快，其应用领域正在不断扩大。我们将在本章第九节中以较大的篇幅讨论这种潮流问题。
缺点：
本算法的主要缺点是收敛速度很慢。病态条件系统，计算往往会发生收敛困难
节点间相角差很大的重负荷系统；包含有负电抗支路(如某些三绕组变压器或线路串联电容等)的系统；具有较长的辐射形线路的系统；长线路与短线路接在同一节点上，而且长短线路的长度比值又很大的系统。此外，平衡节点所在位置的不同选择，也会影响到收敛性能。目前高斯一塞德尔法已很少使用
结合电力系统的特点，对这样的线性网络进行分析，普遍采用的是节点法，节点电压与节点电流之间的关系为：
其展开式分别为：
但是在工程实际中，已知的节点注入量往往不是节点电流而是节点功率，为此必须应用联系节点电流和节点功率的关系式

电力系统分析潮流计算

电力系统分析潮流计算报告目录一．配电网概述 (3)1.1 配电网的分类 (3)1.2 配电网运行的特点及要求 (4)1.3 配电网潮流计算的意义 (5)二．计算原理及计算流程 (5)2.1 前推回代法计算原理 (5)2.2 前推回代法计算流程 (8)2.3主程序清单： (11)2.4 输入文件清单： (13)2.5计算结果清单： (15)三．前推回代法计算流程图 (17)参考文献 (18)一．配电网概述1.1 配电网的分类在电力网中重要起分配电能作用的网络就称为配电网；配电网按电压等级来分类，可分为高压配电网（35—110KV），中压配电网（6—10KV，苏州有20KV的），低压配电网（220/380V）；在负载率较大的特大型城市，220KV电网也有配电功能。

按供电区的功能来分类，可分为城市配电网，农村配电网和工厂配电网等。

在城市电网系统中,主网是指110KV及其以上电压等级的电网，主要起连接区域高压（220KV及以上）电网的作用。

配电网是指35KV及其以下电压等级的电网，作用是给城市里各个配电站和各类用电负荷供给电源。

从投资角度看，我国与国外先进国家的发电、输电、配电投资比率差异很大，国外基本上是电网投资大于电厂投资，输电投资小于配电投资。

我国刚从重发电轻供电状态中转变过来，而在供电投资中，输电投资大于配电投资。

从我国城网改造之后，将逐渐从输电投资转入配电建设为主。

本文是基于前推回代法的配电网潮流分析计算的研究，研究是是以根节点为10kV的电压等级的配电网。

1.2 配电网运行的特点及要求配电系统相对于输电系统来说，由于电压等级低、供电范围小，但与用户直接相连，是供电部门对用户服务的窗口，因而决定了配电网运行有如下特点和基本要求：（1）10kV中压配电网在运行中，负荷节点数多，一般无表计实时记录负荷，无法应用现在传统潮流程序进行配电网的计算分析，要求建立新的数学模型和计算方法。

（2）随着铁道电气化和用户电子设备的大量使用，配电网运行中有大量的谐波源、三相电压不平衡、电压闪变等污染，要求准确测量与计算配电网中的谐波分布，从而采取有效措施抑制配电网运行中的谐波危害。

电力系统分析P-Q分解法潮流计算终稿

《电力系统分析》课程设计报告题目：复杂电力系统潮流分析程序设计所在学院电气工程学院专业班级12级电气工程及其自动化一班学生姓名陈剑秋学生学号************指导教师房大中提交日期 2014年 12 月 16日摘要电力系统潮流计算是电力系统分析课程基本计算的核心部分之一。

它既有本身的独立意义，又是电力系统规划设计，运行分析和理论研究的基础。

电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算，是电力系统计算分析中的一种最基本的计算。

潮流计算是电力系统的各种计算的基础，同时它又是研究电力系统的一项重要分析功能，是进行故障计算，继电保护鉴定，安全分析的工具。

电力系统潮流计算是计算系统动态稳定和静态稳定的基础。

在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用电力系统潮流计算来定量的比较供电方案或运行方式的合理性、可靠性和经济性本课程设计采用P-Q分解法对三机九节点系统进行潮流计算。

计算结果的得出是通过MATLAB软件编程实现。

关键词：电力系统潮流计算，P-Q分解法，三机九节点系统AbstractPower flow calculation is the core part of the power system analysis course basic computing. It has both the independent meaning itself, and is the power system planning and design, the basic operation analysis and theoretical study.Power flow calculation is the calculation of the steady state operation of normal and fault conditions of complex power system, is one of the most basic calculation of power system calculation and analysis of. Power flow calculation is the basis of all kinds of power system calculation, and it is also an important function of power system analysis study, is the fault calculation, relay protection and identification, safety analysis tool. Power flow calculation is the basis of calculation system of dynamic and static stabilization. In the research of operation mode of power system planning and the existing power system, the rationality of power supply, reliability and economy of the program or run mode comparison requires the use of power flow calculation to quantitativeThe curriculum design using P-Q decomposition method for power flow calculation of the three machine nine bus system. It is concluded that the calculation result is achieved by MATLAB software programming.Keywords:power system power flow calculation, the P-Q decomposition method, three machine nine node system目录摘要 (I)Abstract (II)一、绪论 (1)1.1本课题的目的和意义 (1)1.2国内外发展现状 (1)1.2.1高斯-赛德尔迭代法 (2)1.2.2牛顿-拉夫逊法和P-Q分解法 (2)1.2.3基于MATLAB 的电力系统潮流计算发展前景 (3)二、设计目的 (4)三、设计要求和设计指标 (4)四、设计内容 (4)4.1 选题内容 (4)4.2 基础资料 (5)4.3算法原理 (6)4.3.1 节点导纳矩阵形成的计算机方法 (6)4.3.2 电力网络的潮流方程 (9)4.3.3 电力网络极坐标形式的潮流方程 (10)4.3.4 极坐标牛顿潮流算法的修正方程雅和可比矩阵 (10)4.3.5 PQ分解法潮流计算 (11)五、程序设计 (13)5.1 主函数 (13)5.2 子函数 (16)六、输入与输出结果 (21)6.1 输入数据 (21)6.2 输出结果 (22)6.3 当bus4-6的发生故障被切断后，系统的运行情况 (25)6.4 两机五节点网络 (30)七、结果分析 (34)7.1 三机九节点系统正常运行结果分析 (34)7.2 三机九节点系统故障切除后结果分析 (35)结论 (35)参考文献 (36)一、绪论1.1本课题的目的和意义电力系统潮流计算是对复杂电力系统正常和故障条件下稳态运行状态的计算。

(最新整理)(完整版)电力系统交直流电力系统潮流计算

2021/7/26
2
一、直流系统的标幺制
潮流计算时，交流系统通常采用标幺制，因此直流系统也应采用标幺制。因此，需要将换流器的基本方程化为标幺制下的形式以与交流系统相连接。本书选取直流系统的基准功率和基准电压与交流系统相等，即
S dcB S B
V dcB V B
式中： SdcB、SB 分别是直流系统和交流系统的基准功率；VdcB、VB分别为直流系统和交流系统的基准电压。
⑵ 直流系统运行在控制方式二、即整流侧定最小触发角、逆变侧定
电流控制。即
min
Id Ids Id
2021/7/26
13
在Kr和Ki已知的条件下，由于触发角已知，故由整流侧向逆变侧作直流电量计算。计算顺序如下：
首先计算整流侧电量
然后计算逆变侧直流电量
2021/7/26
14
Pd、 r Pd、 i Qd、 r Qd作i 为输出，将用于交流潮流的下一次迭代中。
SVC的基本元件为晶闸管控制的电抗器 TCR和晶闸管投切的电容器TSC。为了降低 SVC的造价，大多数SVC通过降压变压器接入系统。由于阀的控制作用，SVC将产生谐波电流，为降低谐波污染，SVC中还要有滤波器。
2021/7/26
17
1、TCR的工作原理及数学模型
TCR支路由电抗器与两个背靠背连接的晶闸管相串联构成。通过控制晶闸管的触发延迟角，可以控制每个周波内电感L接入系统的时间长短，从而改变TCR的等值电抗。
2021/7/26
9
交直流电力系统的潮流问题可按照牛顿法求解传统潮流的计算流程求解。
三、交直流潮流的顺序解法
顺序解法的基本思想是：迭代计算过程中，将交流系统潮流方程和直流系统潮流方程分别单独进行求解。在求解交流系统方程时，将直流系统换流站处理成接在相应交流节点上的一个等效 P.Q 负荷。而在求解直流系统方程时，将交流系统模拟成加在换流站交流母线上的一个恒定电压。

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

电力系统分析潮流计算实验报告姓名：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;endendendendendfor 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);end% 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;endendendendend% 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:nC(p)=0;for 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、3710i 0、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、0611i1、0000 + 0、0000i 3、0000 + 0、0000i -0、2427 - 0、1970i 0、2556 + 0、1952i 0、0129 - 0、0018i1、0000 + 0、0000i 4、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软件有了初步得了解，对它强大得矩阵运算能力有了更深得体会，同时掌握了设置断点与断点调试得一般方法，结合课本上得程序流程图与参考资料上得例子单步跟踪调试，再一次得熟悉了牛顿拉夫逊法潮流计算得一般方法与步骤，对计算机计算潮流计算有了更进一步得认识，在学习潮流计算时，虽然依次学习了节点导纳矩阵，功率方程、雅可比矩阵，但不能将它们联系起来，更不知道其中得原委，通过程序得编写，知道了其中得联系，也知道了每个方程、矩阵在计算中得作用。