電力系統(tǒng)分析實(shí)驗(yàn)報(bào)告

南昌大學(xué)實(shí)驗(yàn)報(bào)告學(xué)生 號(hào) 實(shí)驗(yàn)類型：□驗(yàn)證□綜合■設(shè)計(jì)□創(chuàng) 電力系統(tǒng)故障分析計(jì)23每組計(jì)算機(jī)1臺(tái)、軟件1套將事先編制好的形成電力網(wǎng)數(shù)學(xué)模型的計(jì)算程序原代碼由自備移動(dòng)設(shè)1、自編程序如下%%%%%%%主程序clearformatshort;switchcase'兩相短路接地'case單相斷開(kāi)case單相短路case兩相短路

case'兩相斷開(kāi)'Vc=exp(2i*pi/3)*Vf1+exp(-%%%%%%%%短路故障計(jì)算函數(shù)function[Eeq,Zff]=xql(h)ifh==1disp('請(qǐng)輸入正序網(wǎng)絡(luò)參數(shù)B1=input('請(qǐng)輸入支路參數(shù)矩陣elseifh==2fori=1:nforifforifB1(i,k)~=0for

ififY(i,j)=-

iflength(find(B4~=0))~=0forifelseifB4(i1)==1form=1:n2forifm==pifm~=p

%%%%%%%%斷線故障計(jì)算函數(shù)function[Eeq,Zff]=xqt(h)ifh==1disp('請(qǐng)輸入正序網(wǎng)絡(luò)參數(shù)elseifh==2fori=1:nforifforifB1(i,k)~=0for

ififY(i,j)=-

iflength(find(B2~=0))~=0forifelseifB4(i1)==1form=1:n2forifm==pifm~=p

2、應(yīng)用計(jì)算例題例驗(yàn)證程序的計(jì)算效運(yùn)行輸入?yún)?shù)如請(qǐng)輸入支路參數(shù)矩陣:B1=[00.510.51000.590請(qǐng)輸入電勢(shì)源參數(shù)矩陣:B2=[1.080.8;00;1.2請(qǐng)輸入節(jié)點(diǎn)參數(shù)矩陣:B3=[0.21.17;00;4請(qǐng)輸入節(jié)點(diǎn)電勢(shì)源數(shù):B4=[20運(yùn)行結(jié)果如下If1If2=-If0=-Vf1Vf2Vf0Ia=-1.1102e-Ib=-1.1814-1.7460iIc=-1.1814+1.7460iVa=0.9237Vb=5.5511e-17-1.1102e-16i運(yùn)行輸入?yún)?shù)如請(qǐng)輸入故障點(diǎn):f=[23]請(qǐng)輸入支路參數(shù)矩陣:B1=[00.5100.5100000000.59請(qǐng)輸入電勢(shì)源參數(shù)矩陣:B2=[1.080.8;00;00;1.2請(qǐng)輸入節(jié)點(diǎn)參數(shù)矩陣:B3=[0.21.17;00;00;4請(qǐng)輸入節(jié)點(diǎn)電勢(shì)源數(shù):B4=[200請(qǐng)輸入故障點(diǎn):f=[1請(qǐng)輸入支路參數(shù)矩陣:B1=[00;0請(qǐng)輸入節(jié)點(diǎn)電勢(shì)源數(shù):B4=[11]運(yùn)行結(jié)果如下If1If2=-If0=-Vf1Vf2Vf0=0.0297Ia=6.9389e-18Ib=-0.0282-0.0362iIc=-0.0282+0.0362iVa=0.0892Vb=1.3878e-17Vc=1.3878e-（說(shuō)明：If1、If2、If0分別是特殊相的正序、負(fù)序、零序電流分量，Vf1、Vf2、Vf0分別是特殊相的正序、負(fù)序、零序電壓分量，Ia、Ib、Ic分別是a、b、c三相的短路電流，Va、Vb、Vc分別是a、b、c三相的短路電壓。需要注意的是，由于的精度問(wèn)題，運(yùn)行結(jié)果中的Ia、Vb、Vc的值為無(wú)窮小，其值可If1If2=-If0=-Vf1Vf2Vf0IaIb=-1.18-1.74iIc=-1.18+1.74iVa=0.923VbVcIf1If2=-If0=-Vf1Vf2Vf0IaIb=-0.028-Ic=-0.028+VbVc理解上講述的各種不對(duì)稱故障的各序網(wǎng)絡(luò)圖