電力系統(tǒng)計(jì)算機(jī)輔助分析試驗(yàn)報(bào)告.doc

_附錄： 學(xué)生實(shí)驗(yàn)報(bào)告表頭格式電力工程 學(xué)院學(xué)生實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)課程名稱： 電力系統(tǒng)計(jì)算機(jī)輔助分析 開課實(shí)驗(yàn)室： 計(jì)算中心408 2014 年 12 月 11 日年級(jí)、專業(yè)、班學(xué)號(hào)姓名成績實(shí)驗(yàn)項(xiàng)目 名 稱電力系統(tǒng)計(jì)算機(jī)輔助分析指 導(dǎo) 教 師教 師 評 語 教師簽名：年 月 日實(shí)驗(yàn)一 MATLAB 軟件的基本操作及程序的基本結(jié)構(gòu)一、程序流程、源碼及仿真結(jié)果啟動(dòng)MATLAB軟件，識(shí)別出常用的四個(gè)窗口命令窗口（Command Window）、歷史命令窗口（Command History）、工作空間瀏覽器（Workspace Browser）。熟悉Matlab工作環(huán)境。2、上機(jī)求解下面的電路問題，腳本文件已經(jīng)給出。已知某個(gè)電路的端電壓的表達(dá)式為：區(qū)間時(shí)，試?yán)L制電壓的波形。提示：在MATLAB的M文件編輯器中鍵入以下命令語句，并保存為ex.m：Clear,clc,closet=0:0.5:10;u=6*exp(-2*t); plot(t,u);title(RC電路的電壓響應(yīng)曲線);xlabel(時(shí)間/s);ylabel(電壓/v)grid%5、 對于如下電路圖，已知V，自己編寫程序，畫出和的波形（四個(gè)周期）。提示：寫出的向量表達(dá)式用于的計(jì)算（寫出串聯(lián)段和并聯(lián)段的阻抗，利用分壓公式），寫出和的幅值和相角就能寫出波形表達(dá)式，設(shè)定波形的長度（四個(gè)周波）和采樣頻率（兩點(diǎn)間隔）后作圖。注意畫圖時(shí)使用標(biāo)題，橫縱坐標(biāo)及網(wǎng)格線等命令。3.%clear,clc,closemUsm=200; % 的幅值aUsm=(pi/180)*(+90); % 的相角Usm=mUsm*exp(j*aUsm); % 用于計(jì)算w=100*pi;r1=22;l=0.05;r2=33;c=0.0001;z1=r1+j*w*l;z2=r2/(j*w*c*r2+1);Ucm=Usm*z2/(z1+z2);mUcm=abs(Ucm) %求幅值aUcm=angle(Ucm) %求相角 T=4*2*pi/w; %用w表示 ，T表示四個(gè)周波N=100; %畫圖的點(diǎn)數(shù)dt=T/N; %兩點(diǎn)間隔，即采樣頻率t=0:dt:T;for k=1:101 %畫波形us(k)=mUsm*cos(w*t(k)+aUsm);uc(k)=mUcm*cos(w*t(k)+aUcm);endplot(t,us,t,uc) %注意區(qū)分兩條曲線的線型和顏色grid ontitle()xlabel(,fontsize,16) % fontsize,16表示橫坐標(biāo)的字體大小為16ylabel()legend(us,uc)%一、 討論1、 MATLAB軟件是什么？有什么特點(diǎn)？答： Matlab是一個(gè)高級(jí)的矩陣/陣列語言，它包含控制語句、函數(shù)、數(shù)據(jù)結(jié)構(gòu)、輸入和輸出和面向?qū)ο缶幊烫攸c(diǎn)。用戶可以在命令窗口中將輸入語句與執(zhí)行命令同步，也可以先編寫好一個(gè)較大的復(fù)雜的應(yīng)用程（M文件）后再一起運(yùn)行。 特點(diǎn)：1. 編程效率高，2. 用戶使用方便，3. 擴(kuò)充能力強(qiáng)，交互性好，4. 移植性和開放性很好，5. 語句簡單，內(nèi)涵豐富，6. 高效方便的矩陣和數(shù)組運(yùn)算，7. 方便的繪圖功能。3. 如何獲得命令或函數(shù)的幫助文檔？ 答：MATLAB 的各個(gè)函數(shù)，不管是內(nèi)建函數(shù)、M文件函數(shù)、還是MEX文件函數(shù)等，一般它們都有M文件的使用幫助和函數(shù)功能說明，各個(gè)工具箱通常情況下也具有一個(gè)與工具箱名相間的M文件用來說明工具箱的構(gòu)成內(nèi)容等，在MATLAB命令窗口中，可以通過指令來獲取這些純文本的幫助信息。 通常能夠起到幫助作用，獲取幫助信息的指令有help、lookfor、which、doc、get、type等。實(shí)驗(yàn)二電力系統(tǒng)計(jì)算中常用的數(shù)值算法及電力網(wǎng)絡(luò)的數(shù)學(xué)模型1、 用因子表法求解電路方程中的節(jié)點(diǎn)電壓； A= 1/10+1/15 -1/10 -1/15 -1 -1/10 1/10+1/50+1/30 -1/50 0 -1/15 -1/50 1/15+1/50+1/30 0 1 0 0 0 ;B=0 0 0 10;n,m=size(A);for i=1:n A(i,i)=1/A(i,i); for j=i+1:n A(i,j)=A(i,j)*A(i,i); end for k=i+1:n for j=i+1:n A(k,j)=A(k,j)-A(k,i)*A(i,j); end endenddisp(矩陣A的因子表為：);disp(A)for i=1:n B(i)=B(i)*A(i,i); for j=i+1:n B(j)=B(j)-A(j,i)*B(i); endendfor i=n-1:-1:1 for j=i+1:-1:2 B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1); endenddisp(在因子表的基礎(chǔ)上求解線性方程組的解為：x=);disp(B)2、將圖2.1中的電壓源分別改為2V、4V、6V、8V時(shí)，用因子表進(jìn)行求解，體會(huì)該算法帶來的益處； 電壓源改為2V電壓源改為4V電壓源改為6V電壓源改為8V3、 上機(jī)驗(yàn)證課本例17.617.7的計(jì)算結(jié)果；例17.6例17.74、 設(shè)計(jì)系統(tǒng)變更時(shí)修改導(dǎo)納矩陣的函數(shù)（參考課本P.323），并調(diào)試通過，給出導(dǎo)納矩陣修正函數(shù)的調(diào)用規(guī)則及其主要的實(shí)現(xiàn)代碼。% n=input(請輸入節(jié)點(diǎn)數(shù):n=);% nl=input(請輸入支路數(shù):nl=);% B=input(請輸入由支路參數(shù)形成的矩陣:B=);% X=input(請輸入由節(jié)點(diǎn)號(hào)及其對地阻抗形成的矩陣:X=);n=4;nl=4;B=1 2 0.08+0.4i 0 1 0;2 3 0.1+0.4i 0 1 0;3 4 0.3i 0 1 0;1 3 0.12+0.5i 0 1 0;X=1 0;2 0;3 0;4 0Y=zeros(n);for i=1:n if X(i,2)=0; p=X(i,1); Y(p,p)=1./X(i,2); end endfor i=1:nl if B(i,6)=0 p=B(i,1); q=B(i,2); else p=B(i,2); q=B(i,1); end Y(p,q)=Y(p,q)-1./(B(i,3)*B(i,5); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B(i,3)*B(i,5)2)+B(i,4)./2; Y(p,p)=Y(p,p)+1./B(i,3)+B(i,4)./2;end disp(導(dǎo)納矩陣Y=:); disp(Y) B1=input(請輸入追加支路參數(shù)形成的矩陣:B1=);% B1=3 4 0.3i 0 1.2 0;% B1=1 2 0.3i 0 1 0;% B1=1 4 0.3i 0 1 0;% B1=3 5 0.3i 0 1 0;z=B1(3);% if B1(6)=0 % p=B1(1); q=B1(2);% else % p=B1(2); q=B1(1);% endp=B1(1); q=B1(2);if B1(6)=1 B1(5)=1/(B1(5);end if B1(2)n Y(p,p)=Y(p,p)+1/z; Y(p,q)=-1/z; Y(q,p)=-1/z; Y(q,q)=(1/z)+B1(4);endif Y(p,q)=0 Y(p,p)=Y(p,p)+1/z; Y(q,q)=Y(q,q)+1/z; Y(p,q)=-1/z; Y(q,p)=-1/z;endif Y(p,q)=0&B1(2)=n for i=1:nl if B(i,1)=B1(1)&B(i,2)=B1(2) if B(i,5)=B1(5) Y(p,p)=Y(p,p)+1/z; Y(q,q)=Y(q,q)+1/z; Y(p,q)=Y(p,q)-1/z; Y(q,p)=Y(q,p)-1/z; else K1=B1(5);K=B(i,5); Y(p,p)=Y(p,p)+(K12-K2)*(1/z);Y(q,q)=Y(q,q); Y(p,q)=Y(p,q)-(K1-K)*(1/z); Y(q,p)=Y(q,p)-(K1-K)*(1/z); end end endend if B1(2)n n=n+1; Y1=zeros(n); for i=1:n for j=1:n Y1(i,j)=Y(i,j); end endenddisp(新的導(dǎo)納矩陣Y1=); disp(Y1);仿真結(jié)果：X = 1 0 2 0 3 0 4 0導(dǎo)納矩陣Y=: 0.9346 - 4.2949i -0.4808 + 2.4038i -0.4539 + 1.8911i 0 -0.4808 + 2.4038i 1.0690 - 4.7568i -0.5882 + 2.3529i 0 -0.4539 + 1.8911i -0.5882 + 2.3529i 1.0421 - 7.5773i 0 + 3.3333i 0 0 0 + 3.3333i 0 - 3.3333i請輸入追加支路參數(shù)形成的矩陣:B1=3 4 0.3i 0 1.2 0;新的導(dǎo)納矩陣Y1= 0.9346 - 4.2949i -0.4808 + 2.4038i -0.4539 + 1.8911i 0 -0.4808 + 2.4038i 1.0690 - 4.7568i -0.5882 + 2.3529i 0 -0.4539 + 1.8911i -0.5882 + 2.3529i 1.0421 - 9.0440i 0 + 4.0000i 0 0 0 + 4.0000i 0 - 3.3333i思考題：1、高斯消去法與因子表法有何異同？答： 因子表法是以高斯消去法為基礎(chǔ)求解線性方程組的一種方法，當(dāng)線性方程組的系數(shù)矩陣不變，而常數(shù)項(xiàng)發(fā)生變化時(shí)，采用因子表法可避免系數(shù)矩陣的重復(fù)運(yùn)算。因子表法首先要建立因子表，然后用因子表對常數(shù)項(xiàng)進(jìn)行消去、回代，從而求出線性方程組的解。2、節(jié)點(diǎn)導(dǎo)納陣與節(jié)點(diǎn)阻抗陣之間有何關(guān)系，二者之間如何轉(zhuǎn)換？電力系統(tǒng)的網(wǎng)絡(luò)方程式一般都用節(jié)點(diǎn)方程式表示，節(jié)點(diǎn)方程有導(dǎo)納型和阻抗型兩種。節(jié)點(diǎn)導(dǎo)納矩陣具有對稱、稀疏、可以根據(jù)網(wǎng)絡(luò)接線圖直接寫出、易于修改的特點(diǎn)；節(jié)點(diǎn)阻抗矩陣是節(jié)點(diǎn)導(dǎo)納矩陣的逆陣，它具有對稱、滿陣的特點(diǎn)，它不能由網(wǎng)絡(luò)接線圖直接寫出，通常由支路追加法形成。 實(shí)驗(yàn)三短路電流計(jì)算合并課本第18章例程。即要求實(shí)現(xiàn)計(jì)算指定節(jié)點(diǎn)各種故障條件下，各支路的故障電流分布情況。其中對于實(shí)用算法，可只提供計(jì)算電抗。提示：盡量使用先前實(shí)驗(yàn)中已封裝好的函數(shù)。2.用新程序驗(yàn)證第18章的各例題。（注意短路電流曲線繪制時(shí)，至少包含57個(gè)點(diǎn)）W=3;switch W case 1 clcclear allNF=1for i1=1:NFclear;n1=3n2=3n0=3nl1=5nl2=5nl0=5Lf=1 %單相接地短路為1 橫向故障為0f=3If=0zf=0zg=0B1=0 1 0.15i 0 1 0;0 2 0.075i 0 1 0;1 2 0.1i 0 1 0;1 3 0.1i 0 1 0;2 3 0.1i 0 1 0B2=0 1 0.15i 0 1 0;0 2 0.075i 0 1 0;1 2 0.1i 0 1 0;1 3 0.1i 0 1 0;2 3 0.1i 0 1 0B0=0 1 0.05i 0 1 0;0 2 0.025i 0 1 0;1 2 0.2i 0 1 0;1 3 0.2i 0 1 0;2 3 0.2i 0 1 0A1=1 1 1;1 -1./2-i*sqrt(3)/2 -1./2+i*sqrt(3)/2;1 -1./2+i*sqrt(3)/2 -1./2-i*sqrt(3)/2Vcs=1;1;1Z1=zeros(n1);Z2=zeros(n2);Z0=zeros(n0);Y1=zeros(n1);Y2=zeros(n2);Y0=zeros(n0);V1=zeros(n1,1);V2=zeros(n2,1);V0=zeros(n0,1);I1=zeros(nl1,1);I2=zeros(nl2,1);I0=zeros(nl0,1);for m1=1:3m=0;if m1=1nl=nl1;B=B1;Z=Z1;n=n1;elseif m1=2nl=nl2;B=B2;Z=Z2;n=n2;elseif m1=3nl=nl0;B=B0;Z=Z0;n=n0;endm=0;Z=zeros(n);V=zeros(n);I=zeros(nl);for k1=1:nlp=B(k1,1);q=B(k1,2);if B(k1,6)=0k=1./B(k1,5);else k=B(k1,5);endif p=0if qm %追加接地樹支Z(q,q)=B(k1,3);m=m+1;else %追加接地連支for i=1:mZ(i,m+1)=-Z(i,q);Z(m+1,i)=-Z(q,i);endZ(m+1,m+1)=Z(q,q)+B(k1,3);for i=1:mfor j=1:mZ(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1);endZ(i,m+1)=0endfor i=1:m+1Z(m+1,i)=0endendelse if qm %追加不接地樹支for i=1:mZ(i,q)=Z(i,p)*k;Z(q,i)=Z(p,i)*k;endZ(q,q)=k2*Z(p,p)+k2*B(k1,3);m=m+1;elsefor i=1:m %追加不接地連支Z(i,m+1)=k*Z(i,p)-Z(i,q);Z(m+1,i)=k*Z(p,i)-Z(q,i);endZ(m+1,m+1)=k2*Z(p,p)+Z(q,q)-2*k*Z(p,q)+k2*B(k1,3);for i=1:m for j=1:mZ(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1);endZ(i,m+1)=0;endfor i=1:m+1Z(m+1,i)=0;end end end end if m1=1Z1=Z;elseif m1=2Z2=Z; elseif m1=3Z0=Z;endend %以上程序求出正,負(fù),零序?qū)?yīng)的阻抗陣if If=0 %求出橫向故障的阻抗陣Z1(f,f)=Z1(f,f);Z2(f,f)=Z2(f,f);Z0(f,f)=Z0(f,f);else %求出縱向故障的阻抗陣Z1(f,f)=Z1(f,f)+Z1(1,1)-2*Z1(f,1);Z2(f,f)=Z2(f,f)+Z2(1,1)-2*Z2(f,1);Z0(f,f)=Z0(f,f)+Z0(1,1)-2*Z0(f,1);end %-故障點(diǎn)自阻抗-if Lf=1 %求出正負(fù)零序短路節(jié)點(diǎn)的電流標(biāo)幺值Z4=Z2(f,f)+Z0(f,f)+3*zf;k2=1;k0=1;I1(f)=Vcs(f)./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);elseif Lf=2 %兩相短路Z4=Z2(f,f)+2*zf;k2=-1;k0=0;I1(f)=Vcs(f)./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);elseif Lf=3 %兩相接地短路Z4=zf+(Z2(f,f)+zf)*(Z0(f,f)+zf+3*zg)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg);k2=-(Z0(f,f)+zf+3*zg)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg);k0=-(Z2(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg);I1(f)=Vcs(f)./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);elseif Lf=4 %單相斷線Z4=zf+(Z2(f,f)+zf)*(Z0(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf);k2=-(Z0(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf);k0=-(Z2(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf);I1(f)=Vcs(f)./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);elseif Lf=5 %兩相斷線Z4=Z2(f,f)+Z0(f,f)+3*zf;k2=1;k0=1;I1(f)=Vcs(f)./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);enddisp(故障點(diǎn)處的正序電流標(biāo)幺值為:);disp(a1);disp(故障點(diǎn)處的負(fù)序電流標(biāo)幺值為:);disp(a2);disp(故障點(diǎn)處的零序電流標(biāo)幺值為:);disp(a0);for i=1:n1 %求出各節(jié)點(diǎn)的電壓標(biāo)幺值 V1(i)=Vcs(i)-I1(f)*Z1(i,f); %正序endfor i=1:n2V2(i)=-I2(f)*Z2(i,f); %負(fù)序endfor i=1:n0V0(i)=-I0(f)*Z0(i,f); %零序enddisp(各節(jié)點(diǎn)的正序電壓標(biāo)幺值(節(jié)點(diǎn)號(hào)從小到大排):);disp(V1);disp(各節(jié)點(diǎn)的負(fù)序電壓標(biāo)幺值(節(jié)點(diǎn)號(hào)從小到大排):);disp(V2);disp(各節(jié)點(diǎn)的零序電壓標(biāo)幺值(節(jié)點(diǎn)號(hào)從小到大排):);disp(V0);for i=1:3 %求各支路的各序電流標(biāo)幺值if i=1B=B1;V=V1;nl=nl1;elseif i=2B=B2;V=V2;nl=nl2;elseif i=3B=B0;V=V0;nl=nl0;endfor j=1:nlif B(j,6)=0k=B(j,5);else k=1./B(j,5);endp=B(j,1);q=B(j,2);if p=0e=0;b=B(j,3);I(j)=(e-V(q)./k)./b; %接地支路電流elseI(j)=(V(p)-V(q)./k)./B(j,3); %不接地支路電流endendif i=1I1=I;elseif i=2I2=I;elseif i=3I0=I;endenddisp(各支路的正序電流標(biāo)幺值為(順序同您輸入B時(shí)一樣):);for i=1:nldisp(I1(i);enddisp(各支路的負(fù)序電流標(biāo)幺值為(順序同您輸入B時(shí)一樣):);for i=1:nldisp(I2(i);enddisp(各支路的零序電流標(biāo)幺值為(順序同您輸入B時(shí)一樣):);for i=1:nldisp(I0(i);endfor i=1:3 %求故障點(diǎn)A,B,C三相的總電流,各個(gè)節(jié)點(diǎn)A,B,C三相各相的總電壓,各支路A,B,C三相各相的總電流if i=1 %故障點(diǎn)電流s0=a0;s1=a1;s2=a2;elseif i=2 %各序支路電壓s0=V0;s1=V1;s2=V2;else %各序支路電流s0=I0;s1=I1;s2=I2;endfor j=1:3 %不對稱系數(shù)CA=s0*A1(j,1)+s1*A1(j,2)+s2*A1(j,3);if j=1A=CA;elseif j=2B=CA;elseif j=3C=CA;endendif i=1Iadl=A;Ibdl=B;Icdl=C;elseif i=2Vajy=A;Vbjy=B;Vcjy=C;else Iazl=A;Ibzl=B;Iczl=C;endenddisp(故障點(diǎn)處A相電流Iadl為:);disp(Iadl);disp(故障點(diǎn)處B相電流Ibdl為:);disp(Ibdl);disp(故障點(diǎn)處C相電流Icdl為:);disp(Icdl);disp(各節(jié)點(diǎn)處A相電壓Vajy分別為(節(jié)點(diǎn)號(hào)從小到大排):);disp(Vajy);disp(各節(jié)點(diǎn)處B相電壓Vbjy分別為(節(jié)點(diǎn)號(hào)從小到大排):)disp(Vbjy);disp(各節(jié)點(diǎn)處C相電壓Vcjy分別為(節(jié)點(diǎn)號(hào)從小到大排):);disp(Vcjy);disp(各支路A相電流Iazl分別為(順序同你輸入B時(shí)一樣):);for i=1:nldisp(Iazl(i);enddisp(各支路B相電流Ibzl分別為(順序同你輸入B時(shí)一樣):);for i=1:nldisp(Ibzl(i);enddisp(各支路C相電流Iczl分別為(順序同你輸入B時(shí)一樣):);for i=1:nldisp(Iczl(i);endendcase 2 NF=1;n=4;nl=5;B=0 1 0.2i 0 1 0;0 2 4i 0 1 0;1 3 0.51i 0 1 0;2 3 0.59i 0 1 0;3 4 1.43i 0 1 0 ;D=4,0;V0=1;1;1;1;m=0;Z=zeros(n);V=zeros(n);I=zeros(nl);for k1=1:nl p=B(k1,1);q=B(k1,2); if B(k1,6)=0 k=1./B(k1,5); else k=B(k1,5); end if p=0 if qm Z(q,q)=B(k1,3);m=m+1; else for i=1:m, Z(i,m+1)=-Z(i,q);Z(m+1,i)=-Z(q,i); end Z(m+1,m+1)=Z(q,q)+B(k1,3); for i=1:m for j=1:m Z(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1); end Z(i,m+1)=0; end for i=1:m+1 Z(m+1,i)=0; end end else if qm for i=1:m Z(i,q)=Z(i,p)*k;Z(q,i)=Z(p,i)*k; end Z(q,q)=k2*Z(p,p)+k2*B(k1,3); m=m+1; else for i=1:m Z(i,m+1)=k*Z(i,p)-Z(i,q);Z(m+1,i)=k*Z(p,i)-Z(q,i); end Z(m+1,m+1)=k2*Z(p,p)+Z(q,q)-2*k*Z(p,q)+k2*B(k1,3); for i=1:m for j=1:m Z(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1); end Z(i,m+1)=0; end for i=1:m+1 Z(m+1,i)=0; end end endendfor k=1:NF I(D(k,1),D(k,1)=V0(D(k,1),1)./(Z(D(k,1),D(k,1)+D(k,2);ft=num2str(D(k,1);ts1=(點(diǎn)短路時(shí));ts2=(電流的標(biāo)么值If=); dn=strcat(ft,ts1,ts2);disp(dn);disp(I(D(k,1),D(k,1); for i=1:n V(i,i)=V0(i,1)-I(D(k,1),D(k,1)*Z(i,D(k,1); end for i=1:nl if B(i,6)=0 K=B(i,5); else k=1./B(i,5); end p=B(i,1);q=B(i,2); if p=0 e=0;b=B(i,3); I(i,i)=(e-V(q,q)./k)./b; else I(i,i)=(e-V(p,p)-V(q,q)./k)./B(i,3); end end disp(各節(jié)點(diǎn)的電壓標(biāo)么值U為：); for i=1:n disp(V(i,i); end disp(各支路短路電流的標(biāo)么值I為：); for i=1:nl disp(I(i,i); endend case 3 clcNF=1;n=3;nl=5;Sb=100;D=3 0;vb=115 115 115;V0=1;1;1;B=0 1 0.15i 0 1 0 100 1 1;0 2 0.075i 0 1 0 200 1 1;1 2 0.1i 0 1 0 0 0 2;1 3 0.1i 0 1 0 0 0 2;2 3 0.1i 0 1 0 0 0 2;N=5;m=0;Z=zeros(n);sq3=sqrt(3);s=zeros(N,1);for k1=1:nl p=B(k1,1);q=B(k1,2); if B(k1,6)=0 k=1./B(k1,5); else k=B(k1,5); end if p=0 if qm Z(q,q)=B(k1,3);m=m+1; else for i1=1:m, Z(i1,m+1)=-Z(i1,q);Z(m+1,i1)=-Z(q,i1); end Z(m+1,m+1)=Z(q,q)+B(k1,3); for i1=1:m for j=1:m Z(i1,j)=Z(i1,j)-Z(i1,m+1)*Z(m+1,j)./Z(m+1,m+1); end Z(i1,m+1)=0; end for i1=1:m+1 Z(m+1,i1)=0; end end else if qm for i1=1:m Z(i1,q)=Z(i1,p)*k;Z(q,i1)=Z(p,i1)*k; end Z(q,q)=k2*Z(p,p)+k2*B(k1,3); m=m+1; else for i1=1:m Z(i1,m+1)=k*Z(i1,p)-Z(i1,q);Z(m+1,i1)=k*Z(p,i1)-Z(q,i1); end Z(m+1,m+1)=k2*Z(p,p)+Z(q,q)-2*k*Z(p,q)+k*B(k1,3); for i1=1:m for j=1:m Z(i1,j)=Z(i1,j)-Z(i1,m+1)*Z(m+1,j)./Z(m+1,m+1); end Z(i1,m+1)=0; end for i1=1:m+1 Z(m+1,i1)=0; end end endendVb=zeros(1,n);V=zeros(1,n);Vd=zeros(1,n);Ib=zeros(1,nl);for r=1:NF Idb=V0(D(r,1),1)./(Z(D(r,1),D(r,1)+D(r,2); ts1=(以下是短路點(diǎn)); ft=NUm2str(D(r,1); ts2=(各時(shí)刻的結(jié)果); dn=strcat(ts1,ft,ts2); disp(dn); for k=1:n Vb(k)=V0(k,1)-Idb*Z(k,D(r,1); end for i1=1:nl p=B(i1,1);q=B(i1,2); if p=0&B(i1,8)=0 if B(i1,6)=0 k=B(i1,5); V(i1)=Vb(p)-Vb(q)./k else k=1./B(i1,5); V(i1)=k*Vb(p)-Vb(q); end else V(i1)=1-Vb(q); end Ib(i1)=V(i1)./B(i1,3); end disp(0s時(shí)短路點(diǎn)的電流為); disp(Idb); disp(0s時(shí)各節(jié)點(diǎn)的電壓標(biāo)么值為); disp(Vb); disp(0s時(shí)各支路的電流標(biāo)么值為); disp(Ib) for a=1:N Is=zeros(1,nl);Isy=zeros(1,nl);Idzs=0; s(a)=input(請輸入短路后幾秒（若要終止本程序請輸入-1）s=); if s(a)=-1 disp(本程序已被終止，謝謝使用！); break end for k=1:nl if B(k,1)=0&B(k,8)=1 Sn(k)=B(k,7); cf=Idb./Ib(k);Z1(k)=cf*Z(D(r,1),D(r,1);X(k)=abs(imag(Z1(k); Xjs1(k)=X(k)*Sn(k)./Sb; if s4|Xjs1(k)3.45|B(k,9)=0 Is(k)=1./Xjsl(k); if B(k,9)=0 Isy(k)=Is(k)*Sb./(sq3*vb(D(r,1); tsl=(無窮大電源); ft=num2str(k); ts2=(對短路點(diǎn)的計(jì)算電抗Xjsl，對應(yīng)的電流標(biāo)么值，有名值分別為); dn=strcat(tsl.ft,ts2);disp(dn); disp(Xjsl(k);disp(Is(k);disp(Isy(k); else Isy(k)=Is(k)*Sn(k)./(sq3*vb(D(r,1); tsl=(等值電源); ft=num2str(k); ts2=(對短路點(diǎn)的計(jì)算電抗Xjsl，對應(yīng)的電流標(biāo)么值，有名值分別為); dn=strcat(tsl,ft,ts2); disp(dn); disp(Xjsl(k);disp(Is(k);disp(Isy(k); end else Xjs=Xjs1(k); tsl=(等值電源); ft=num2str(k); ts2=(對短路點(diǎn)的計(jì)算電抗Xjsl為); dn=strcat(tsl,ft,ts2);disp(dn); disp(Xjs) Is(k)=input(查與計(jì)算電抗Xjs對應(yīng)的電流標(biāo)么值I=); Isy(k)=Is(k)*Sn(k)./(sq3*vb(D(r,1); disp(對應(yīng)的電流標(biāo)么值、有名值（單位：KA）分別為);disp(Is(k);disp(Isy(k); end Idzs=Idzs+Isy(k); end end Id(a)=Idzs; sj=num2str(s(a);ts=(秒對應(yīng)的短路點(diǎn)短路總電流為（單位：kA）);jg=strcat(sj,ts); disp(jg) disp(Id(a); endplot(s,Id);xlabel(時(shí)間（秒）);ylabel(短路電流大小的標(biāo)么值);title(短路電流曲線);endotherwisebreakend思考題1.計(jì)算短路電流的目的是什么？答：選擇電氣設(shè)備,計(jì)算各種保護(hù)的定值,制定運(yùn)行方案等 2. 在計(jì)算短路電流時(shí)，是否用到了節(jié)點(diǎn)阻抗陣中的所有元素？據(jù)此，可以如何簡化短路計(jì)算程序？ 答：不需要，可以利用高斯消元法消去部分中間節(jié)點(diǎn)和聯(lián)絡(luò)節(jié)點(diǎn)，化簡成等效網(wǎng)絡(luò)