电力系统短路故障分析的MATLAB辅助程序设计短路计算程序Word格式文档下载.docx

电力系统短路故障分析的MATLAB辅助程序设计短路计算程序Word格式文档下载.docx

《电力系统短路故障分析的MATLAB辅助程序设计短路计算程序Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《电力系统短路故障分析的MATLAB辅助程序设计短路计算程序Word格式文档下载.docx（41页珍藏版）》请在冰豆网上搜索。

程序sｙmfauｌt（zdata,Zbｕｓ,Ｖ）用来计算三相对称故障，程序要求输入ｚｄata和Ｚｂｕs两个矩阵,第三个参数V是可选的。

如果V不存在，程序将默认故障前所有的节点电压标幺值为１.0，如果变量Ｖ存在,那么V包括节点编号和复数电压值。

电压向量V也可以由潮流计算程序自动生成。

当syｍfａult程序运行时,用户要输入故障节点编号和故障阻抗,运行可得到总的故障电流，节点电压幅值以及故障情况下的线路电流。

在三相短路计算中，zｂｕild和ｓｙmfaｕｌt程序，ｚbｕildpi和ｓymｆaｕｌt程序都可以进行计算，下面是三相短路计算使用的程序代码:

（1）Ｚbｕｉｌｄ.ｍ程序代码：

function[Zｂｕs]　=　zbuｉld（linｅdaｔa）

nl=　linedata（:

１）;

　nr＝　ｌinｅdaｔa（：

２）;

R=lｉnｅｄatａ（:

，３）;

X　=lｉneｄatａ（:

4）;

nbr=ｌｅnｇth（liｎedａta（：

1））;

　nbuｓ＝max（maｘ（nｌ）,ｍaｘ（ｎr））;

forｋ＝1:

ｎbr

iｆ　R（k）　==inf|　X（k）==inｆ

R（k）　=999９99９99;

Ｘ（ｋ）　=99９９999９9;

％无穷

　　else，end

end

ZB＝R　+　j*Ｘ;

Ｚｂuｓ=zeros（ｎbus,nbuｓ）；

tｒee=0;

%从参考总线0上添加一个分支

　ｆorI=1:

ｎtｒee（I）＝１;

　　ifnl（I）＝=0|nr（I）＝＝　0

　　　if　nl（I）=＝0　　　n　=　ｎr（I）；

　　ｅlseifnｒ（Ｉ）　==0n　=ｎｌ（I）;

　　ｅｎｄ

　　　　ifaｂs（Zbuｓ（ｎ，　ｎ））==　0　　Ｚｂus（n，ｎ）　=ZB（I）;

tｒeｅ＝ｔｒeｅ＋1;

％％ｎｅw

　　eｌseZｂus（n,n）　=Zbｕｓ（n,n）＊ZＢ（I）/（Zbｕs（n,ｎ）　+　ZB（I））;

　　end

　　　　ntrｅｅ（Ｉ）=２；

　　　　else，ｅｎd

enｄ

%添加一个新总线分支到现有总线上

whilｅ　ｔree<

nbuｓ

　ｆorn=1:

ｎbus

　　　　nａｄｄ=　１;

　　ｉfabs（Zbus（n，n））==0

　　for　I=1:

　　　　　ｉfnadd==１；

　　　　ｉfnｌ（I）==n｜ｎｒ（I）　＝=n

　　　　　　ifｎｌ（I）=＝n　ｋ　＝nｒ（I）;

　　　elｓｅifｎr（I）==n　k=　nl（I）;

　　　　　end

　　　　　　ifabs（Zbus（k，k））　～=　0

　　　　　　　fｏｒm　=1:

nbus

　　　　　iｆ　ｍ～=n

　　　　　　　　Zbｕｓ（m,n）=Zbｕs（m，k）;

　　　　　　　　Zbus（ｎ，m）=Zbus（m,ｋ）；

　　　　　　　　elｓe,end

　　　　eｎd

　　　　　　　Zbｕs（n,n）=　Ｚbus（k,k）+ＺＢ（I）;

　ｔrｅe=ｔｒｅe+1;

　　　　　　　　　　　　　　nａｄｄ=　２;

ｎtｒee（I）=２；

　　　　　else,　end

　　　　　else,end

　　　　ｅlｓe，ｅnd

　ｅｎd

　elsｅ，ｅnd

eｎd

%增加两个原有总线间的支路阻抗

fｏr　n=１:

ｎｂus

　forI＝１：

ｎbｒ

ifntｒee（I）　＝=１

　　　　ｉfｎｌ（I）==　n|nr（I）==n

　　　ifｎl（Ｉ）　==n　　ｋ=　nr（I）;

　　　elsｅifｎｒ（I）　==　nk=ｎl（Ｉ）;

　　　end

　　　　　DＭ=Zbus（n,n）＋Zbｕｓ（k,k）　＋　ZB（I）-２*Zbｕs（n,k）；

　　　　ｆor　jj=１:

nｂus

　　　　　　　AP　=　Zbuｓ（jj,ｎ）-Zbｕs（jj，k）;

　　　　foｒkk=1：

nbｕs

　　　　　　ＡT　=　Zｂus（ｎ，kk）　-Zbus（k,ｋk）;

　　　　　ＤＥLZ（jj，ｋk）=AＰ＊AT/DＭ；

　　　　　　　　eｎｄ

　　　end

　　　　　　Ｚbuｓ＝Ｚｂｕｓ-DELＺ;

　　　　　ntrｅe（I）=2;

　　　eｌse,ｅnd

　　　　eｌｓe，end

　　enｄ

　　ｅnd

（2）Zbuilｄpi.m程序代码:

%Thisprｏgｒamformｓthｅcoｍpleｘbｕsｉmpｅｄanｃe　matｒixｂythe　mｅtｈoｄ

%ｏfbuｉldingalgorｉｔhm.　　Buszero　istaken　aｓreｆerence.

%　Thispｒｏgramｉscoｍｐatible　wｉtｈ　powerfｌow　ｄaｔａ.

fｕncｔion[Zbus,liｎeｄaｔa]=　zbuildpi（lｉnedａta,gendatａ,ｙloａd）

nｇ=length（geｎｄatａ（:

，1））;

nlg=genｄatａ（:

nrg=gendata（：

2）;

zｇ＝ｇｅndａta（:

2）+　ｊ*geｎdａｔａ（：

3）;

nｌ＝linedata（:

1）；

　nr=ｌineｄata（：

R＝linｅｄata（:

，3）；

Ｘ＝　liｎedatａ（:

４）；

nｂｒ＝ｌenｇth（ｌinedaｔａ（:

nbus＝　ｍax（mａx（nl），maｘ（nｒ））;

nc=leｎｇｔh（liｎｅdａｔａ（1,:

））；

for　ｋ=1:

nbr

　ifR（k）＝＝　ｉnf|　X（k）＝=inf

　Ｒ（k）＝9９９9999９;

Ｘ（k）＝　9９999999;

　else,　end

ｅｎd

if　nc＞4

ＢC=ｌinedatａ（：

5）；

　　for　ｎ＝１：

　　yc（n）　=0;

　　ｎｌc（n）=　0;

nrc（ｎ）=n；

　fork　＝１:

　　　ifnｌ（ｋ）==n|nr（k）　==　n

　　　　　yc（ｎ）=yｃ（n）+ｊ*BC（k）;

　　　　　　　ｅlse，ｅnd

　end

eｎd

elsｅｉf　nc=＝4　yc=　zeros（1,nbr）;

　end

ｎｌc＝nlc＇;

ｎrc=nrc'

;

ｙｃ=yｃ.'

；

ＺB=R+j*X;

iｆ　ｅxist（'

yloａd'

）=＝　1

　ylｏａd＝ｙｌｏaｄ．＇；

ｙc　=yc＋　yload;

elｓe,ｅnｄ

m　=0；

ｆorｎ=　1:

ｉｆabs（yc（n））　~=0

m=m＋1;

　nlcc（ｍ）=ｎlc（ｎ）;

　ｎrｃc（m）＝　nrｃ（n）；

　zc（m）＝1/yc（n）；

　else,eｎd

nlcｃ=ｎlcc'

nrcc＝ｎrcc＇；

zc＝zc.＇;

nｌ=[ｎlg;

nlcc；

nl];

nr=　[nrｇ；

nrcc；

nr];

　ＺB=[zg;

　zc;

　ZB];

lｉneｄatａ=［ｎｌｎr　real（ＺB）imag（ZＢ）]；

nbｒ=ｌengｔh（nｌ）；

Zbus＝zｅrｏｓ（nｂuｓ,nbus）;

tｒee=０;

%%从参考总线0上添加一个分支

fｏr　I＝1:

nbｒ

　ｎtrｅe（I）=　１;

ifnl（I）　==０|nｒ（Ｉ）　==0

　　ifnl（I）＝=０　　　n=nr（I）；

　　elseifnｒ（Ｉ）　=＝0ｎ=nl（I）;

　　　end

　　ifaｂs（Zｂｕs（n,ｎ））==0Ｚｂｕs（n,n）=ZＢ（I）;

tｒee=ｔｒee+1;

　　　　elsｅZbuｓ（n,n）=　Zbｕｓ（n,ｎ）*ＺB（I）/（Zbus（n，n）　+ZB（I））；

　　　　　　eｎd

　ｎｔｒee（I）=2；

　　　eｌｓｅ,enｄ

ｅnd

%%　添加一个新总线分支到现有总线上

whiｌｅｔrｅｅ<

　nbus

forn=1:

　　nadｄ＝　１；