2020年昆明理工大学计算机辅助分析课程设计及matlab自动实现潮流调节与计算

昆明理工大学计算机辅助分析课程设计及

matlab自动实现潮流调节与计算

资料仅供参考

重要声明：

1、本文只是提供一定参考，有错勿怪。

2.、本文采用的潮流调节及计算只针对本实验，其中有很多固化

的数据与参数，以及特殊的公式、请根据要求的不同进行一定

的调整（潮流计算的程序适用任何情况、不同的只是系统参数）。

3、本文潮流调节为自动调节生成，调节变量有发电机2、3电

压，三台变压器分接头，没有进行无功补偿，请根据具体情况

加入或删减变量。

4、自动调节潮流程序中，控制变量为各点电压首次迭代值、最

终各点电压。最终选择标准为网损最低的一组数据。

5、本实验报告为草版，有错误，但懒得改。

6、鉴于个人所有文档免费得原则，本文为免费文档，请尽量下

载。

7、最后！都免费了，要求不要太高。

资料仅供参考

票求及原始数据

电力系统潮流计算是电力系统分析中的一种最基本的计算,是对复杂电力系统正

常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行

状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是

否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行

和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都

是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮

流等对潮流计算的模型和方法有直接影响。

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

障。产生短路故障的主要原因是电力设备绝缘损坏。短路故障分为三相短路、两相

短路、单相接地短路及两相接地短路。其中三相短路时三相电流依然对称，其余三

类短路统成为不对称短路。短路故障大多数发生在架空输电线路。电力系统设计与

运行时，要采取适当的措施降低短路故障的发生概率。短路计算能够为设备的选择

提供原始数据。

一、课程设计的目的与任务：

该课程设计为电气工程及其自动化专业的主要专业课设计之一，是学习电力系

统分析基础和电力系统计算机辅助分析课程后的一个重要的实践性教学环节。其目

的在于巩固和加深对电力系统潮流和短路电流计算基本原理的理解，学习和掌握应

用计算机进行电力系统设计和计算的方法，培养学生独立分析和解决问题的能力。

二、课程设计的基本要求：

掌握电力系统等值模型和参数计算，以及潮流和短路计算的基本原理，学会应

用计算机计算系统潮流分布和短路电流的方法。

三、课程设计选题原则：

该课程设计是根据电力系统分析基础和电力系统计算机辅助分析课程内容，结

合实际工程和科研的电力系统网络进行系统的潮流和短路电流计算。

四、课程设计内容：

资料仅供参考

1、课程设计地点：机房和教室

2、设计内容

（1）电力系统潮流计算（2）电力系统短路电流计算

具体内容包括：

1）系统等值模型和参数计算。

2）利用现有的潮流和短路电流计算程序进行计算，或编制电力系统潮流和短路计算

程序。

3）计算结果分析。

4）编写说明书。

3、时间安排

1、系统等值模型和参数计算（2天）

2、读懂现有的算法程序或视需要编制算法程序并进行计算（4天）

3、计算结果整理和分析（1天）

4、编写说明书（2天）

5、成绩考核（1天）

五、课程设计主要参考资料：

1、李光琦.电力系统暂态分析.北京：水利电力出版社.第二版

2、陈玲.电力系统稳态分析.北京：水利电力出版社.第二版

3、孟祥萍.基于MATLAB的电力系统设计与计算.长春：吉林人民出版社

4、孟祥萍.电力系统分析.北京：高等教育出版社

5、电力系统潮流和短路计算课程设计指导书.自编

六、课程设计成绩考核：

根据学生的平时表现，说明书，以及设计完成后的口试成绩等进行综合评定。

课程结束时，要求学生有课程设计说明书。课程设计成绩以平时成绩（占50%）、

设计说明书（占50%）进行综合评定。

资料仅供参考

原始数据与思考题

目标电网单线图

bus-2必-7bus-9,

bus-8bus-3o

bus-5bus-6

bus-4

bus-1

系统参数

表1.线路参数表

线线号线路长线路电线路线

路度（km）阻访半

编{Q领{S/km}

号/km）{Q

/km}

4-5LGJ-240/301040.0471.78X

ICY

资料仅供参考

4-6LGJ-120/701100.0741.47X

io-6

5-7LGJ-120/25980.0791.60X

0.4io-6

6-9LGJ-95/551630.0921.80X

10-6

7-8LGJ-240/30880.0471.78X

io-6

8-9LGJ-240/301220.0471.78X

io-6

说明：线路零序电抗为正序电抗3倍。

表2.变压器参数表

线牌号变压器型号变压器变比短路电压

(kV)百分数(％)

2-7SSPL-20242+3X10

2.5%/20

3-9SSPL-10242±3X6

2.5%/15

1-4SSPL-240000242±3X11

2.5%/17.5

说明：变压器零序电抗与正序电抗相等，且均为△

资料仅供参考

/Yo接法。

表3.发电机参数表1

发电机额定功率额定电压额定功率

{MW}{kV}因数

120016.50.85

2180180.85

310013.80.85

表4.发电机参数表2

发母兀{。Z{Qx,{QTq。%

电线))(S)))(S)(S)

机名

110.1460.0608.960.0960.096047.2

08998

220.8950.1196.000.8640.1960.5312.8

885950

331.3120.1818.591.2570.2500.606.02

53800

表5.负荷数据表

有功负荷无功负荷

节点号

(MW)(MVA)

资料仅供参考

512550

69030

810035

思考题

1、牛顿拉夫逊法与PQ分解法有哪些联系？有哪些区别？二者的计算性能如何？

2、为什么在用计算机对某网络初次进行潮流计算时往往是要调潮流，而并非任何情

况下只一次送入初始值算出结果就行呢？要考虑什么条件？各变量是如何划分

的？哪些可调？哪些不可调？

3、发电机节点的注入无功为负值说明了什么？

4、线路有功潮流最有可能的流向是？线路无功潮流最有可能的流向是？

5、负荷功率因数对系统潮流有什么影响？

6、潮流控制的主要手段有哪些？

7、如何降低系统网损？

8、如何校验主设备的短路容量？

9、利用节点阻抗矩阵计算短路电流和利用导纳矩阵计算短路电流的算法有何异

同？

10、减小短路电流的措施有哪些？

设计主要内容：

1、完成调试潮流计算与短路计算程序调试；

计算程序以《电力系统计算机辅助分析》所用教材中的例程为基础，要求掌握

算法原理，建议进行适当改写以利于本次设计应用。

2、根据设计资料准备计算数据；

整理原始数据，使其符合计算程序的输入条件。要求比较采用变压器II型模型

前后两种情况的差异。

3、为指定电网安排合理的稳态运行方式；

资料仅供参考

合理设置目标电网中的控制变量，以1号机为平衡机，将潮流分布调整为可行

的运行方式（要求所有节点的电压不得低于1.0P.U.,也不能高于1.05p.u.）绘制

潮流分布图。分析该运行方式的特点，提出进一步降低网损的改进的措施，并利用

计算结果简单证明之。具体内容参考设计资料所附思考题。

4、列表分析节点7发生各种短路故障时的短路电流，以及运行方式对故障电流的影

响；

5、完成设计说明书。要求结构清晰，即从电网模型建立，到初态潮流分布计算，到

问题分析，再到调整后的方案计算。故障电流计算也类似。不必逐一回答所有思考

题。

参数计算、潮流调节、与潮流计算

matlab源程序

%参数计算的

资料仅供参考

clear;

L=[104,110,98,163,88,122];

r=[0.047,0.074,0.079,0.092,0.047,0.047];

b=[1.78e-006,1.47e-006,1.60e-006,1.80e-006,1.78e-006,1.78e-006];

dispC线路电抗X=');

Xl=100/220A2*0.4*L;%线路各段电抗(均为标幺值，取

Ub=220KV,SB=100MVA,下同)

disp(Xl);

fori=l:6

R(i)=100/220A2*L(i)*r(i);%加分号避免出现过程项，下同

end%线路各段电阻

dispC线路电阻R=');disp(R);%输出R

fori=l:6

Z(i)=R(i)+j*Xl(i);

Z0(i)=R(i)+3*j*Xl(i);

end

dispC线路正序、负序阻抗Z1=Z2=');disp(Z);

fori=l:6

BHi)=220A2/100*b⑴*L⑴；

B(i)=0+j*Bl(i);

end%线路各容纳之半

dispC线路容纳之半B/2=');disp(B);%输出B

dispC线路零序阻抗Z0:');

disp(ZO);

U0=[ll,10,6];

SN=[240,220,120];

UN=[242,242,242];

fori=l:3;

XT(i)=j*U0(i)*242A2/SN(i)/220A2;

XT0(i)=XT(i);

end;

dispC变压器正、负序电抗值Xt：，);disp(XT);

dispC变压器零序电抗值XtO：');disp(XT0);

x=[0.0608,0.H98,0.1813];%发电机暂态电抗有名值

UGN=[16.5,18,13.8];

fori=l:3

K(i)=242/UGN(i);

end

fori=l:3

Xdg(i)=j*100/220A2*x(i)*K(i)*K(i);

end%发电机暂态电抗

dispC发电机暂态电抗Xdg=)disp(Xdg);%输出Xdg

PL=[125,90,100];

QL=[50,30,35];

资料仅供参考

fori=l:3

ZL(i)=100/(PL(i)A2+QL(i)A2)*(PL(i)+j*QL(i));

end%负荷节点阻抗

dispC负荷节点阻抗Zl=');disp(ZL);%输出ZL

%支数参数形成的矩阵Bl、支数参数形成的矩阵B2、节点号及其对地阻抗形

成的矩阵X的形成

Bl=[l,4,XT(l),0,l,0;2,7,XT(2),0,l,0;3,9,XT(3),0,l,0；4,5,R(l)+Xl(l)*j,B(l),l,0；4,6

,取2)+

Xl(2)*j3(2),l,0；5,7,R(3)+Xl(3)*j,B(3),l,0；6,9,R(4)+Xl(4)*j,B(4),l,0;7,8,R(5)+Xl(5)*j

,B(5),l,0；8,9,R(6)+Xl(6)*j,B(6),l,0;];

dispC支数参数形成的矩阵：Bl=);disp(Bl);%输出支数参数形成的矩阵

B1

P=[200,180,100];

W=(P/100)+((P/0.85)*sqrt(l-0.85*0.85)/100)*j;

B2=[W(1),0,105,1,0,1；W⑵,0,1,1,03W(3),0,l,l,0,3;0,0,l,0,0,2；0,

1.25+0.51,l,0,0,2；0,0.9+0.3i,l,0,0,2；0,0,l,0,0,2；0,l+0.35i,l,0,0,2；0,0,l,0,0,2；l；

dispC支数参数形成的矩阵：B2=);disp(B2);%输出支数参数形成的矩阵

B2

%潮流自动调节的与自动计算的

n=9;

nl=9;

isb=l;

pr=0.0001;

mm=0;

X=[l,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0;9,0];

'dispC只考虑各点电压水平，得到的满足要求的数据如下：,);

foraa=0:5

B2(2,4)=l+0.01*aa;

forab=0:5

B2(3,4)=l+0.01*ab;

forba=0:5

Bl(l,5)=0.95+0.025米ba;

forbb=0:5

Bl(2,5)=0.95+0.025*bb;

forbc=0:5

Bl(3,5)=0.95+0.025*bc;

Y=zeros(n);e=zeros(l,n);f=zeros(l,n);V=zeros(l,n);O=zeros(l,n);Sl=zeros(nl);

资料仅供参考

fori=l:n

ifX(i,2)~=0;

P=X(i,l);

Y(p,p)=l./X(i,2);

end

end

fori=l:nl

ifBl(i,6)==0

P=Bl(i,l)；q=Bl(i,2);

elsep=Bl(i,2);q=Bl(i,l)；

end

Y(p,q)=Y(p,q)-l./(Bl(i,3)*Bl(i,5));

Y(q,p)=Y(p,q);

Y(q,q)=Y(q,q)+l./(Bl(i,3)*Bl(i,5)A2)+Bl(i,4)./2；

Y(p,p)=Y(p,p)+l./Bl(i,3)+Bl(i,4)./2;

end

G=real(Y);B=imag(Y);

fori=l:n

e(i)=real(B2(i,3));

f(i)=imag(B2(i,3))；

V(i)=B2(i,4);

end

fori=l:n

S(i)=B2(i,l)-B2(i,2);

B(i,i)=B(i,i)+B2(i,5);

end

P=real(S);Q=imag(S);

ICTl=0;IT2=l;N0=2*n;N=N0+l;a=0;

whileIT2~=0

IT2=0;a=a+l;

fori=l:n

ifi~=isb

C(i)=0;

D(i)=0;

forjl=l:n

C(i)=C(i)+G(ijl)*eai)-B(iJl)*f(jl)；

D(i)=D(i)+G(Ml)*f(jl)+B(ijl)*e(jl)；

end

Pl=C(i)*e(i)+f(i)*D(i);

Ql=f(i)*C(i)-D(i)*e(i);

V2=e(i)A2+f(i)A2;

ifB2(i,6)~=3

DP=P(i)-Pl;

DQ=Q(i)-Ql;

资料仅供参考

forjl=l:n

ifjl~=isb&jl~=i

Xl=-G(ijl)*e(i)-B(ijl)*f(i);

X2=B(iJl)*e(i)-G(iJl)*f(i);

X3=X2;

X4=-X1;

p=2*i-l;q=2*jl-l;J(p,q)=X3;J(p,N)=DQ;m=p+l;

J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X4;J(m,q)=X2;

elseifjl==i&jl~=isb

Xl=«i)-G(i,i)*e(i)-B(i,i)*f(i);

X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);

X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i)；

X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);

p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;

J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X4;J(m,q)=X2;

end

end

else

DP=P(i)-Pl;

DV=V(i)A2-V2;

forjl=l:n

ifjl~=isb&jl~=i

Xl=-G(ijl)*e(i)-B(ijl)*f(i);

X2=B(ijl)*e(i)-G(iJl)*f(i);

X5=0;

X6=0;

p=2*i-l;q=2*jl-l;J(p,q)=X5;J(p,N)=DV;m=p+l;

J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X6;J(m,q)=X2;

elseifjl==i&jl~=isb

Xl=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);

X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i)；

X5=-2*e(i);

X6=-2*f(i);

p=2*i-l;q=2*jl-l;J(p,q)=X5;J(p,N)=DV;m=p+l;

J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X6;J(m,q)=X2;

end

end

end

end

end

fork=3:N0

kl=k+l;Nl=N;

fork2=kl:Nl

J(k,k2)=J(k,k2)./J(k,k);

资料仅供参考

end

J(k,k)=l;

ifk~=3

k4=k-l;

fork3=3:k4

fork2=kl:Nl

J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

ifk==NO,break;end

fork3=kl:N0

fork2=kl:Nl

J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

else

fork3=kl:N0

fork2=kl:Nl

J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

end

end

fork=3:2:N0-l

L=(k+l)./2;

e(L)=e(L)-J(k,N);

kl=k+l;

f(L)=f(L)-J(kl,N);

end

fork=3:N0

DET=abs(J(k,N));

ifDET>=pr

IT2=IT2+1;

end

end

ICT2(a)=IT2;

ICT1=ICT1+1;

fork=l:n

dy(k)=sqrt(e(k)A2+f(k)A2);

end

fori=l:n

资料仅供参考

Dy(ICTl,i)=dy(i);

end

forvb=l:n

ifDy(l,vb)<=1.05&Dy(l,vb)>=1.00;

zy(vb)=l;

elsezy(vb)=0;

end

end

end%用高斯消去法解"w=・J*V”

fork=l:n

V(k)=sqrt(e(k)A2+f(k)A2);

O(k)=atan(f(k)./e(k))*180Jpi;

end

E=e+f*j;

forva=l:n

ifV(va)<=1.05&V(va)>=1.000;

yj(va)=l;

elseyj(va)=0;

end

end

ifyj(l)*yj(2)*yj(3)*yj(4)*yj(5)*yj(6)*yj(7)*yj(8)*yj(9)==l

ifzy(l)*zy(2)*zy(3)*zy(4)*zy(5)*zy(6)*zy(7)*zy(8)*zy(9)==l;

dispC发电机2出口电压：();

disp(B2(2,4));

dispC发电机3出口电压：,);

disp(B2(3,4));

dispC变压器1接头电压：,);

disp(Bl(l,5));

dispC变压器2接头电压：，)；

disp(B1(2,5));

dispC变压器3接头电压：，)；

disp(B1(3,5));

dispC各节点的电压大小V为(节点号从小到大排列)：上

disp(V);

forp=l:n

C(p)=0;

forq=l:n

C(p)=C(p)+conj(Y(p,q))*conj(E(q));

end

S(p)=E(p)*C(p);

end

fori=l:nl

ifBl(i,6)==0

p=Bl(i,l);q=Bl(i,2);

资料仅供参考

elsep=Bl(i,2);q=Bl(i,l)；

end

Si(p,q)=E(p)*(conj(E(p))*conj(Bl(i,4)./2)+(conj(E(p)*Bl(i,5))-conj(E(q)))*conj(l

・/(Bl(i,3)*Bl(i,5))));

end

fori=l:nl

ifBX(i,6)==0

p=Bl(M)；q=Bl(i,2);

elsep=Bl(i,2);q=Bl(iJ);

end

Sj(q,p)=E(q)*(conj(E(q))*conj(Bl(i,4)./2)+(conj(E(q)./Bl(i,5))-conj(E(p)))*conj(l

./(Bl(i,3)*Bl(i,5))));

end

ZDS=0;

fori=l:nl

ifBl(i,6)==0

p=Bl(i4)；q=Bl(i,2);

elsep=Bl(i,2);q=Bl(i,l)；

end

DS(i)=Si(p,q)+Sj(q,p);

ZDS=DS(i)+ZDS;

end

mm=mm+l;

ZDSS(mm)=sqrt(real(ZDS)A2+imag(ZDS)A2);

fori=l:n

ZB224(mm)=B2(2,4)；

ZB234(mm)=B2(3,4)；

ZB115(mm)=Bl(l,5);

ZB125(mm)=Bl(2,5);

ZB135(mm)=Bl(3,5);

end

end

end

end

end

end

end

end

'dispC满足各点电压要求的数据组共有：，)；

disp(mm);

资料仅供参考

LOL=min(ZDSS);

fordota=l:inni

ifZDSS(dota)==LOL

mm=dota;

end

end

dispC当加入线路损耗时，得到其中线路损耗最小的一组组号及关键数据为：

disp(mm);

dispC发电机2出口电压：,);

disp(ZB224(mm));

dispC发电机3出口电压：，)；

disp(ZB234(mm));

dispC变压器1接头电压：,);

disp(ZB115(mm));

dispC变压器2接头电压：');

disp(ZB125(mm));

dispC变压器3接头电压：，)；

disp(ZB115(mm));

B2(2,4)=ZB224(mm);

B2(3,4)=ZB234(mm);

Bl(l,5)=ZB115(mm);

Bl(2,5)=ZB125(mm);

Bl(3,5)=ZB135(mm);

dispC用此组数据算的的潮流结果为：，)；

dispC节点数：');

disp(n);

dispC节点支路数：,)；

disp(nl);

dispC平衡节点：,);

disp(isb);

dispC精度：,)；

disp(pr);

disp('Bl=');

disp(Bl);

disp('B2=');

disp(B2);

disp('X=');

disp(X);

Y=zeros(n);e=zeros(l,n);f=zeros(l,n);V=zeros(l,n);0=zeros(l,n);Sl=zeros(nl);

dispC节点数：，)；

disp(n);

dispC节点支路数：)

disp(nl);

dispC平衡节点：,);

资料仅供参考

disp(isb);

dispC精度：');

disp(pr);

mm=O;

X=[l,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0;9,0];

fori=l:n

ifX(i,2)~=0;

P=X(i,l)；

Y(p,p)=l./X(i,2);

end

end

fori=l:nl

ifBl(i,6)==0

P=Bl(i,l)；q=Bl(i,2);

elsep=Bl(i,2);q=Bl(i,l);

end

Y(p,q)=Y(p,q)-l./(Bl(i,3)*Bl(i,5));

Y(q,p)=Y(p,q);

Y(q,q)=Y(q,q)+l./(Bl(i,3)*Bl(i,5)^2)+Bl(i,4)./2;

Y(p,p)=Y(p,p)+l./Bl(i,3)+Bl(i,4)./2;

end

G=real(Y);B=imag(Y);

fori=l:n

e(i)=real(B2(i,3));

f(i)=imag(B2(i,3))；

V(i)=B2(i,4);

end

fori=l:n

S(i)=B2(i,l)-B2(i,2);

B(i,i)=B(i,i)+B2(i,5);

end

P=real(S);Q=imag(S);

ICTl=0;IT2=l;N0=2*n;N=N0+l;a=0;

whileIT2~=0

IT2=0;a=a+l;

fori=l:n

ifi~=isb

C(i)=0;

D(i)=0;

forjl=l:n

C(i)=C(i)+G(ijl)*e(jl)-B(iJl)*f(jl);

D⑴=5i)+G(i,jl)*f(jl)+B(iJl)*e(jl);

end

Pl=C(i)*e(i)+f(i)*D(i);

资料仅供参考

Ql=f(i)*C(i)-D(i)*e(i);

V2=e(i)A2+f(i)A2;

ifB2(i,6)~=3

DP=P(i)-Pl;

DQ=Q(i)-Ql;

forjl=l:n

ifjl~=isb&jl~=i

Xl=-G(ijl)*e(i)-B(ijl)*f(i);

X2=B(ijl)*e(i)-G(ijl)*f(i);

X3=X2;

X4=-X1;

p=2*i-l;q=2*jl-l;J(p,q)=X3;J(p,N)=DQ;m=p+l;

J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X4;J(m,q)=X2;

elseifjl==i&jl~=isb

Xl=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);

X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);

X3=D(i)+B(i,i)*e(i)-G(i4)*f(i)；

X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);

p=2*i-l;q=2*jl-l;J(p,q)=X3;J(p,N)=DQ;m=p+l;

J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X4;J(m,q)=X2;

end

end

else

DP=P(i)-Pl;

DV=V(i)A2-V2;

forjl=l:n

ifjl~=isb&jl~=i

Xl=-G(ijl)*e(i)-B(ijl)*f(i);

X2=B(ijl)*e(i)-G(ijl)*f(i);

X5=0;

X6=0;

p=2*i-l;q=2*jl-l;J(p,q)=X5;J(p,N)=DV;m=p+l;

J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X6;J(m,q)=X2;

elseifjl==i&jl~=isb

Xl=-C(i)-G(i,i)*e(i)-B(i,i)*f(i)；

X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i)；

X5=-2*e(i);

X6=-2*f(i);

p=2*i-l;q=2*jl-l;J(p,q)=X5;J(p,N)=DV;m=p+l;

J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X6;J(m,q)=X2;

end

end

end

end

资料仅供参考

end

fork=3:N0

kl=k+l;Nl=N;

fork2=kl:Nl

J(k,k2)=J(k,k2)./J(k,k);

end

J(k,k)=l；

ifk~=3

k4=k-l;

fork3=3:k4

fork2=kl:Nl

J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

ifk==NO,break;end

fork3=kl:N0

fork2=kl:Nl

J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

else

fork3=kl:N0

fork2=kl:Nl

J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);

end

J(k3,k)=0;

end

end

end

fork=3:2:N0-l

L=(k+l)./2;

e(L)=e(L)-J(k,N);

kl=k+l;

f(L)=f(L)-J(kl,N);

end

fork=3:N0

DET=abs(J(k,N));

ifDET>=pr

IT2=IT2+1;

end

end

ICT2(a)=IT2;

资料仅供参考

ICT1=ICT1+1;

fork=l:n

dy(k)=sqrt(e(k)A2+f(k)A2);

end

fori=l:n

Dy(ICTl,i)=dy(i);

end

end%用高斯消去法解"w=-J*V"

dispC迭代次数，)；

disp(ICTl);

dispC没有达到精度要求的个数，)；

disp(ICT2);

fork=l:n

V(k)=sqrt(e(k)A2+f(k)A2);

O(k)=atan(f(k)./e(k))*180Jpi;

end

E=e+f*j;

dispC各节点的实际电压标幺值E为(节点号从小到大排列)：，)；

disp(E);

dispC各节点的电压大小V为(节点号从小到大排列)：,);

disp(V);

dispC各节点的电压角。为(节点号从小到大排列)：，)；

disp(O);

forp=l:n

C(p)=0;

forq=l:n

C(p)=C(p)+conj(Y(p,q))*conj(E(q));

end

S(p)=E(p)*C(p);

end

dispC各节点的功率S为(节点号从小到大排列)：，)；

disp(S);

dispC各条支路的首段功率Si为(顺序同您输入B1时一样)：力

fori=l:nl

ifBl(i,6)==0

p=Bl(i,l)；q=Bl(i,2);

elsep=Bl(i,2);q=Bl(i,l);

end

Si(p,q)=E(p)*(conj(E(p))*conj(Bl(i,4)./2)+(conj(E(p)*Bl(i,5))-conj(E(q)))*conj(l

J(Bl(i,3)*Bl(i,5))));

disp(Si(p,q));

end

dispC各条支路的末段功率Sj为(顺序同您输入Bl时一样)：，)；

fori=l:nl

资料仅供参考

ifBX(i,6)==0

p=Bl(M);q=Bl(i,2);

elsep=B1(i,2);q=B1(i,l);

end

Sj(q,p)=E(q)*(conj(E(q))*conj(Bl(i,4)y2)+(conj(E(q)./Bl(i,5))-conj(E(p)))*conj(l

J(Bl(i,3)*Bl(i,5))));

disp(Sj(q,p));

end

dispc各条支路的功率损耗DS为(顺序同您输入B1时一样)：工

fori=l:nl

ifBl(i,6)==0

p=Bl(i,l);q=Bl(i,2);

elsep=Bl(i,2);q=Bl(i,l);

end

DS(i)=Si(p,q)+Sj(q,p);

disp(DS(i));

end

fori=l:ICTl

Cs(i)=i;

end

dispC以下是各节点每次迭代

后的电压值(如图所示),);

plot(Cs,Dy),xlabel(，迭代次数lylabelC电压上title。电压迭代次数曲线，)；

最后是课程设计报告

二刃理N大孝

课程设计

资料仅供参考

课程名称：________电力系统分析__________________

设计题目：基于Matlab计算程序的电力系统运行分析

学院：___________电力工程学院________________

专业：___________电气工程自动化______________

年级：____________________________________

学生姓名易军（你能够称我易大神，可惜dota、1。1都是菜

鸟顶级水平）

指导教师：_________宋琪、谢静、单节杉___________

日期：_____________zl-13_____________________

教务处制

前S

电力系统潮流计算是电力系统分析中的一种

最基本的计算，是对复杂电力系统正常和故障条件

下稳态运行状态的计算。潮流计算的目标是求取电

资料仅供参考

力系统在给定运行状态的计算。即节点电压和功率

分布，用以检查系统各元件是否过负荷,各点电压是

否满足要求，功率的分布和分配是否合理以及功率

损耗等。对现有电力系统的运行和扩建，对新的电

力系统进行规划设计以及对电力系统进行静态和

暂态稳定分析都是以潮流计算为基础。潮流计算结

果可用如电力系统稳态研究，安全估计或最优潮流

等对潮流计算的模型和方法有直接影响。

在电力系统中可能发生的各种故障中，危害

最大且发生概率较高的首推短路故障。产生短路故

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

三相短路、两相短路、单相接地短路及两相接地短

路。其中三相短路时三相电流依然对称，其余三类

短路统成为不对称短路。短路故障大多数发生在架

空输电线路。电力系统设计与运行时，要采取适当

的措施降低短路故障的发生概率。短路计算能够为

设备的选择提供原始数据。

目标电网接线图

资料仅供参考

一、电网模型的建立

设计中，选取S=100MVA,S=220KV,将所有支

路的参数都折算到220KV电压等级侧，

本文利用matlab计算参数

其源程序为：

clc;

clear;

L=[104,H0,98,163,88,122];

r=[0.047,0.074,0.079,0.092,0.047,0.047];

b=[1.78e-006,1.47e-006,1.60e-006,1.80e-006,1.78e-0

资料仅供参考

06,1.78e-006];

dispC线路电抗X=9；

Xl=100/220八2*0.4*L;%线路各段电抗(均为

标幺值，Ub=220KV,SB=100MVA)

disp(Xl);

fori=l:6

R(i)=100/220A2*L(i)*r(i);

%加分号避免出现过程项，下同

end

%线路各段电阻

dispC线路电阻RR;

disp(R);

%输出R

fori=l:6

Z(i)=R(i)+j*Xl(i);

Z0(i)=R(i)+3*j*Xl(i);

end

dispC线路正序、负序阻抗Z1=Z2=);

disp(Z);

fori=l:6

Bl(i)=220A2/100*b(i)*L(i);

资料仅供参考

B(i)=O+j«Bl(i);

end

%线路各容纳之半

dispC线路容纳之半B/2+);

disp(B);%输出B

dispC线路零序阻抗Z0:);

disp(ZO);

U0=[ll,10,6];

SN=[240,220,120];

UN=[242,242,242];

fori=l:3;

XT⑴二产UO⑴*242八2/SN⑴/220八2;

XT0(i)=XT(i);

end;

dispC变压器正、负序电抗值Xt：，);disp(XT);

dispC变压器零序电抗值XtO：，);disp(XT0);

x=[0.0608,0.1198,0.1813];

%发电机暂态电抗有名值

UGN=[16.5,18,13.8]；

fori=l:3

K(i)=242/UGN(i);

end

资料仅供参考

fori=l:3

Xdg(i)=j*100/220A2«x(i)*K(i)*K(i);

end

%发电机暂态电抗

dispC发电机暂态电抗Xdg=*);disp(Xdg);

PL=[125,90,100];

QL=[50,30,35];

fori=l:3

ZL(i)=100/(PL(i)A2+QL(i)A2)*(PL(i)+j*QL(i));

end

%负荷节点阻抗

dispC负荷节点阻抗Zl=*);disp(ZL);

计算结果为：

线路电抗X=

0.08600.09090.08100.1347

0.07270.1008

线路电阻R二

0.01010.01680.01600.0310

0.00850.0118

资料仅供参考

线路正序、负序阻抗Z1=Z2=

Columns1through5

0.0101+0.0860i0.0168+0.0909i0.0160+

0.0810i0.0310+0.1347i0.0085+0.0727i

Column6

0.0118+0.1008i

线路容纳之半B/2=

Columns1through5

0+0.0896i0+0.0783i0+

0.0759i0+0.1420i0+0.0758i

Column6

0+0.1051i

线路零序阻抗zo：

Columns1through5

0.0101+0.2579i0.0168+0.2727i0.0160+

0.2430i0.0310+0.4041i0.0085+0.2182i

Column6

0.0118+0.3025i

变压器正、负序电抗值Xt：

0+0.0555i0+0.0550i0+

0.0605i

变压器零序电抗值xto：

资料仅供参考

0+0.0555i0+0.0550i0+

0.0605i

发电机暂态电抗Xdg=

0+0.0270i0+0.0447i0+

0.1152i

负荷节点阻抗Zl=

0.6897+0.2759i1.0000+0.3333i0.8909+

0.3118i

支数参数形成的矩阵：Bl=

Columns1through5

1.00004.0000

0+0.0555i01.0000

2.00007.0000

0+0.0550i01.0000

3.00009.0000

0+0.0605i01.0000

4.00005.0000

0.0101+0.0860i0+0.0896i1.0000

4.00006.0000

0.0168+0.0909i0+0.0783i1.0000

5.00007.0000

资料仅供参考

0.0160+0.0810i0+0.0759i1.0000

6.00009.0000

0.0310+0.1347i0+0.1420i1.0000

7.00008.0000

0.0085+0.0727i0+0.0758i1.0000

8.00009.0000

0.0118+0.1008i0+0.1051i1.0000

Column6

0

0

0

0

0

0

0

0

0

支数参数形成的矩阵：B2=

Columns1through5

2.0000+1.2395i01.0500

1.00000

1.8000+1.1155i01.0000

资料仅供参考

1.00000

1.0000+0.6197i01.0000

1.00000

001.0000

00

01.2500+0.5000i1.0000

00

00.9000+0.3000i1.0000

00

001.0000

00

01.0000+0.3500i1.0000

00

001.0000

00

Column6

1.0000

3.0000

3.0000

2.0000

2.0000

资料仅供参考

2.0000

2.0000

2.0000

2.0000

.•系统等值电路图的绘制

根据以上计算结果，得到系统等值电路图如下:

0.8909+J0.3118

0.0310+J0.07270.0118+J0.1008

J0.0447J0.0550

J0.0605J0.1152

E3

_J0.0758J0,0758J0.1051J0,1051

ii-r