計(jì)算機(jī)輔助分析課程設(shè)計(jì)及matlab自動(dòng)實(shí)現(xiàn)潮流調(diào)節(jié)與計(jì)算

重要聲明：1、本文只是提供一定參考，有錯(cuò)勿怪。2.、本文采用的潮流調(diào)節(jié)及計(jì)算只針對(duì)本實(shí)驗(yàn)，其中有很多固化的數(shù)據(jù)與參數(shù)，以及特殊的公式、請(qǐng)根據(jù)要求的不同進(jìn)行一定的調(diào)整（潮流計(jì)算的程序適用任何情況、不同的只是系統(tǒng)參數(shù)）。3、本文潮流調(diào)節(jié)為自動(dòng)調(diào)節(jié)生成，調(diào)節(jié)變量有發(fā)電機(jī)2、3電壓，三臺(tái)變壓器分接頭，沒有進(jìn)行無功補(bǔ)償，請(qǐng)根據(jù)具體情況加入或刪減變量。4、自動(dòng)調(diào)節(jié)潮流程序中，控制變量為各點(diǎn)電壓首次迭代值、最終各點(diǎn)電壓。最終選擇標(biāo)準(zhǔn)為網(wǎng)損最低的一組數(shù)據(jù)。5、本實(shí)驗(yàn)報(bào)告為草版，有錯(cuò)誤，但懶得改。6、鑒于個(gè)人所有文檔免費(fèi)得原則，本文為免費(fèi)文檔，請(qǐng)盡量下載。7、最后！都免費(fèi)了，要求不要太高。要求及原始數(shù)據(jù)電力系統(tǒng)潮流計(jì)算是電力系統(tǒng)分析中的一種最基本的計(jì)算，是對(duì)復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運(yùn)行狀態(tài)的計(jì)算。潮流計(jì)算的目標(biāo)是求取電力系統(tǒng)在給定運(yùn)行狀態(tài)的計(jì)算。即節(jié)點(diǎn)電壓和功率分布，用以檢查系統(tǒng)各元件是否過負(fù)荷.各點(diǎn)電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對(duì)現(xiàn)有電力系統(tǒng)的運(yùn)行和擴(kuò)建，對(duì)新的電力系統(tǒng)進(jìn)行規(guī)劃設(shè)計(jì)以及對(duì)電力系統(tǒng)進(jìn)行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計(jì)算為基礎(chǔ)。潮流計(jì)算結(jié)果可用如電力系統(tǒng)穩(wěn)態(tài)研究，安全估計(jì)或最優(yōu)潮流等對(duì)潮流計(jì)算的模型和方法有直接影響。在電力系統(tǒng)中可能發(fā)生的各種故障中，危害最大且發(fā)生概率較高的首推短路故障。產(chǎn)生短路故障的主要原因是電力設(shè)備絕緣損壞。短路故障分為三相短路、兩相短路、單相接地短路及兩相接地短路。其中三相短路時(shí)三相電流仍然對(duì)稱，其余三類短路統(tǒng)成為不對(duì)稱短路。短路故障大多數(shù)發(fā)生在架空輸電線路。電力系統(tǒng)設(shè)計(jì)與運(yùn)行時(shí),要采取適當(dāng)?shù)拇胧┙档投搪饭收系陌l(fā)生概率。短路計(jì)算可以為設(shè)備的選擇提供原始數(shù)據(jù)。一、課程設(shè)計(jì)的目的與任務(wù)：該課程設(shè)計(jì)為電氣工程及其自動(dòng)化專業(yè)的主要專業(yè)課設(shè)計(jì)之一,是學(xué)習(xí)電力系統(tǒng)分析基礎(chǔ)和電力系統(tǒng)計(jì)算機(jī)輔助分析課程后的一個(gè)重要的實(shí)踐性教學(xué)環(huán)節(jié)。其目的在于鞏固和加深對(duì)電力系統(tǒng)潮流和短路電流計(jì)算基本原理的理解，學(xué)習(xí)和掌握應(yīng)用計(jì)算機(jī)進(jìn)行電力系統(tǒng)設(shè)計(jì)和計(jì)算的方法，培養(yǎng)學(xué)生獨(dú)立分析和解決問題的能力。二、課程設(shè)計(jì)的基本要求：掌握電力系統(tǒng)等值模型和參數(shù)計(jì)算，以及潮流和短路計(jì)算的基本原理，學(xué)會(huì)應(yīng)用計(jì)算機(jī)計(jì)算系統(tǒng)潮流分布和短路電流的方法。三、課程設(shè)計(jì)選題原則：該課程設(shè)計(jì)是根據(jù)電力系統(tǒng)分析基礎(chǔ)和電力系統(tǒng)計(jì)算機(jī)輔助分析課程內(nèi)容，結(jié)合實(shí)際工程和科研的電力系統(tǒng)網(wǎng)絡(luò)進(jìn)行系統(tǒng)的潮流和短路電流計(jì)算。四、課程設(shè)計(jì)內(nèi)容：1、課程設(shè)計(jì)地點(diǎn)：機(jī)房和教室2,設(shè)計(jì)內(nèi)容（1）電力系統(tǒng)潮流計(jì)算（2）電力系統(tǒng)短路電流計(jì)算具體內(nèi)容包括：1）系統(tǒng)等值模型和參數(shù)計(jì)算。2）利用現(xiàn)有的潮流和短路電流計(jì)算程序進(jìn)行計(jì)算，或編制電力系統(tǒng)潮流和短路計(jì)算程序。3）計(jì)算結(jié)果分析。4）編寫說明書。3、時(shí)間安排1、系統(tǒng)等值模型和參數(shù)計(jì)算（2天）2、讀懂現(xiàn)有的算法程序或視需要編制算法程序并進(jìn)行計(jì)算（4天）3、計(jì)算結(jié)果整理和分析（1天）4、編寫說明書（2天）5、成績(jī)考核（1天）五、課程設(shè)計(jì)主要參考資料：1、李光琦.電力系統(tǒng)暫態(tài)分析.北京：水利電力出版社.第二版2、陳拓.電力系統(tǒng)穩(wěn)態(tài)分析.北京：水利電力出版社.第二版3、孟祥萍.基于MATLAB的電力系統(tǒng)設(shè)計(jì)與計(jì)算.長(zhǎng)春：吉林人民出版社4、孟祥萍.電力系統(tǒng)分析.北京：高等教育出版社5、電力系統(tǒng)潮流和短路計(jì)算課程設(shè)計(jì)指導(dǎo)書.自編六、課程設(shè)計(jì)成績(jī)考核：根據(jù)學(xué)生的平時(shí)表現(xiàn)，說明書，以及設(shè)計(jì)完成后的口試成績(jī)等進(jìn)行綜合評(píng)定。課程結(jié)束時(shí)，要求學(xué)生有課程設(shè)計(jì)說明書。課程設(shè)計(jì)成績(jī)以平時(shí)成績(jī)（占50%）、設(shè)計(jì)說明書（占50%）進(jìn)行綜合評(píng)定。

原始數(shù)據(jù)與思考題目標(biāo)電網(wǎng)單線圖系統(tǒng)參數(shù)表1.線路參數(shù)表繳編號(hào)(km)線犍阻{Q/km}型ITT刪九{fi/km}維循此半{SAm}4-5LGJ-240/301040.0470.41.78X10-64-6LGJ-120/701100.0741.47X10-65-7LGJ-120/25980.0791.60X10-66-9LGJ-95/551630.0921.80XI。/7-8LGJ-240/30880.0471.78XI。/8-9LGJ-240/301220.0471.78XI。/說明：線路零序電抗為正序由抗3倍。表2.變壓器參數(shù)表

線輟號(hào)變壓耀號(hào)變壓敏:匕（kV）短路電壓百分?jǐn)?shù)（％）2-7SSPL-220000242±3X2.5%/20103-9SSPL-120000242±3X2.5%/1561-4SSPL-240000242±3X2.5%/17.511說明：變壓器零序電抗與正序電抗相等，且均為A/Y。接法。表3.發(fā)電機(jī)參數(shù)表1發(fā)電機(jī)額定功率{MW}額定電壓{kV}額定功率因數(shù)120016.50.852180180.85310013.80.85表4.發(fā)電機(jī)參數(shù)表2發(fā)電機(jī)母線名X(/｛。｝x“｛。｝Go(S)X“｛。｝Tq。（S）小S)110.14600.06088.960.09690.0969047.28220.89580.11986.000.86450.19690.53512.80331.31250.18138.591.25780.25000.6006.02表5.負(fù)荷數(shù)據(jù)表節(jié)點(diǎn)號(hào)有功負(fù)荷（MW）無功負(fù)荷（MVA）51255069030810035思考題1、牛頓拉夫遜法與PQ分解法有哪些聯(lián)系？有哪些區(qū)別？二者的計(jì)算性能如何？2、為什么在用計(jì)算機(jī)對(duì)某網(wǎng)絡(luò)初次進(jìn)行潮流計(jì)算時(shí)往往是要調(diào)潮流，而并非任何情況下只一次送入初始值算出結(jié)果就行呢？要考慮什么條件？各變量是如何劃分的？哪些可調(diào)？哪些不可調(diào)？3、發(fā)電機(jī)節(jié)點(diǎn)的注入無功為負(fù)值說明了什么？4、線路有功潮流最有可能的流向是？線路無功潮流最有可能的流向是？5、負(fù)荷功率因數(shù)對(duì)系統(tǒng)潮流有什么影響？6、潮流控制的主要手段有哪些？7、如何降低系統(tǒng)網(wǎng)損？8、如何校驗(yàn)主設(shè)備的短路容量？9、利用節(jié)點(diǎn)阻抗矩陣計(jì)算短路電流和利用導(dǎo)納矩陣計(jì)算短路電流的算法有何異同？10、減小短路電流的措施有哪些?設(shè)計(jì)主要內(nèi)容：1、完成調(diào)試潮流計(jì)算與短路計(jì)算程序調(diào)試；計(jì)算程序以《電力系統(tǒng)計(jì)算機(jī)輔助分析》所用教材中的例程為基礎(chǔ)，要求掌握算法原理，建議進(jìn)行適當(dāng)改寫以利于本次設(shè)計(jì)應(yīng)用。2、根據(jù)設(shè)計(jì)資料準(zhǔn)備計(jì)算數(shù)據(jù)；整理原始數(shù)據(jù)，使其符合計(jì)算程序的輸入條件。要求比較采用變壓器n型模型前后兩種情況的差異。3、為指定電網(wǎng)安排合理的穩(wěn)態(tài)運(yùn)行方式；合理設(shè)置目標(biāo)電網(wǎng)中的控制變量，以1號(hào)機(jī)為平衡機(jī)，將潮流分布調(diào)整為可行的運(yùn)行方式（要求所有節(jié)點(diǎn)的電壓不得低于L0p.U.,也不能高于L05P.U.）繪制潮流分布圖。分析該運(yùn)行方式的特點(diǎn)，提出進(jìn)一步降低網(wǎng)損的改進(jìn)的措施，并利用計(jì)算結(jié)果簡(jiǎn)單證明之。具體內(nèi)容參考設(shè)計(jì)資料所附思考題。4、列表分析節(jié)點(diǎn)7發(fā)生各種短路故障時(shí)的短路電流，以及運(yùn)行方式對(duì)故障電流的影響；5、完成設(shè)計(jì)說明書。要求結(jié)構(gòu)清晰，即從電網(wǎng)模型建立，到初態(tài)潮流分布計(jì)算，到問題分析，再到調(diào)整后的方案計(jì)算。故障電流計(jì)算也類似。不必逐一回答所有思考題。參數(shù)計(jì)算、潮流調(diào)節(jié)、與潮流計(jì)算matlab源程序%參數(shù)計(jì)算的clc;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];disp(線路電抗X=');X1=100/220竺*0.4*L;%線路各段電抗(均為標(biāo)幺值,取Ub=220KV,SB=1OOMVA,下同)disp(Xl);fori=l:6R(i)=100/220A2*L(i)*r(i);%加分號(hào)避免出現(xiàn)過程項(xiàng)，下同end %線路各段電阻disp(線路電阻R=); disp(R);%輸出Rfori=l:6Z(i)=R(i)+j*Xl(i);Z0(i)=R(i)+3*j*Xl(i);enddisp('線路正序、負(fù)序阻抗Z1=Z2=);disp(Z);fori=l:6Bl(i)=220A2/100*b(i)*L(i);B(i)=O+j*Bl(i);end %線路各容納之半disp(線路容納之半B/2=');disp(B); %輸出Bdisp(線路零序阻抗Z0:');disp(ZO);U0=[l1,10,6];SN=[240,220,120];UN=f242,242,242];fori=l:3;XT(i)=j*U0(i)*242A2/SN(i)/220A2;XT0(i)=XT(i);end;disp,變壓器正、負(fù)序電抗值Xt：');disp(XT);dispC變壓器零序電抗值XtO：');disp(XT0);x=[0.0608,0.1198,0.1813]; %發(fā)電機(jī)暫態(tài)電抗有名值UGN=[16.5,18,13.8];fori=l:3K(i)=242/UGN(i);endfori=l:3Xdg(i)=j*100/220A2*x(i)*K(i)*K(i);end %發(fā)電機(jī)暫態(tài)電抗dispC發(fā)電機(jī)暫態(tài)電抗Xdg=');disp(Xdg);%輸出XdgPL=[125,90,100];QL=[5O,3O,35];fori=l:3ZL(i)=100/(PL(i)A2+QL(i)A2)*(PL(i)+j*QL(i));end %負(fù)荷節(jié)點(diǎn)阻抗dispC負(fù)荷節(jié)點(diǎn)阻抗ZIR; disp(ZL);%輸出ZL%支數(shù)參數(shù)形成的矩陣Bl、支數(shù)參數(shù)形成的矩陣B2、節(jié)點(diǎn)號(hào)及其對(duì)地阻抗形成的矩陣X的形成B1=r1,4,XT(1),0,1,0;2,7,XT(2),0,1,0;3,9,XT(3),0,1,0;4,5,R(1)+X1(1)*j,B(1),1,0;4,6,R⑵+Xl(2)*j,B(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支數(shù)參數(shù)形成的矩陣：Bl=');disp(Bl);%輸出支數(shù)參數(shù)形成的矩陣B1P=[200,180,100];W=(P/lOO)+((P/O.85)*sqrt(l-O.85*O.85)/lOO)*j;B2=[W(1),0,105,1,0,1； W(2),0,l,l,0,3; W(3),0,l,1,0,3；0,0,l,0,0,2;0,1.25+0.5i,l,0,0,2;0,0.9+0.3i,1,0,0,2;0,0,1,0,0,2;0,1+0.35i,1,0,0,2;0,0,1,0,0,2;];dispC支數(shù)參數(shù)形成的矩陣：B2=');disp(B2);%輸出支數(shù)參數(shù)形成的矩陣B2%潮流自動(dòng)調(diào)節(jié)的與自動(dòng)計(jì)算的n=9;nl=9;isb=l;pr=0.0001;mm=0;X=[1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0;9,0];dispC*******************************************************************、?disp(，只考慮各點(diǎn)電壓水平，得到的滿足要求的數(shù)據(jù)如下：)；foraa=0:5B2(2,4)=l+0.01*aa;forab=0:5B2(3,4)=l+0.01*ab;forba=0:5Bl(l,5)=0.95+0.025*ba;forbb=0:5Bl(2,5)=0.95+0.025*bb;forbc=0:5Bl(3,5)=0.95+0.025*bc;Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);S1=zeros(nl);fori=l:nifX(i,2)~=0;p=X(i,l);Y(p,p)=l./X(i,2);endendfori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endY(p,q)=Y(p,q)-l./(B1(i,3)*Bl(i,5));Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+l./(B1(i,3)*Bl(i,5)A2)+Bl(i,4),/2;Y(p,p)=Y(p,p)+1./B1(i,3)+Bl(i,4)./2;endG=real(Y);B=imag(Y);fori=l:ne(i)=real(B2(i,3))；f(i)=imag(B2(i,3))；V(i)=B2(i,4);endfori=l:nS(i)=B2(i,l)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICTl=0;IT2=l;N0=2*n;N=N0+l;a=0;whileIT2~=0IT2=0;a=a+l;fori=l:nifi~=isbC(i)=o；D(i)=O;forjl=l:nC(i)=C(i)+G(i,jl)*eGl)-B(i,jl)*f(jl);D(i)=D(i)+G(i,jl)*f(jl)+B(ijl)*e(j1);endPl=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)~=3DP=P(i)-Pl;DQ=Q(i)-Ql;forjl=l:nifjl~=isb&jl~=iXl=-G(i,jl)*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*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;elseifj1==i&j1~=isbXl=-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(i,i)*f(i);X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);p=2*i-l;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;endendelseDP=P(i)-Pl;DV=V(i)A2-V2;forjl=l:nifjl~=isb&jl~=iXl=-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*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+l;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;elseifjl==i&jl~=isbXl=-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*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=Xl;J(m,N)=DP;q=q+I;J(p,q)=X6;J(m,q)=X2;endendendendendfork=3:N0kl=k+l;Nl=N;fork2=kl:NlJ(k,k2)=J(k,k2)./J(k,k);endJ(k,k)=l;ifk~=3k4=k-l;fork3=3:k4fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endifk==NO,break;endfork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endelsefork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendendfork=3:2:N0-lL=(k+l)./2;e(L)=e(L)-J(k,N);kl=k+l;f(L)=f(L)-J(kl,N);endfork=3:N0DET=abs(J(k,N));ifDET>=prIT2=IT2+1;endendICT2(a)=IT2;ICT1=ICT1+1;fork=l:ndy(k)=sqrt(e(k)A2+f(k)A2);endfori=l:nDy(ICTl,i)=dy(i);endforvb=l:nifDy(l,vb)<=1.05&Dy(l,vb)>=1.00;zy(vb)=l;elsezy(vb)=0;endendend%用高斯消去法解"w=?J*V”fork=l:nV(k)=sqrt(e(k)A2+f(k)A2);O(k)=atan(f(k)./e(k))*180./pi;endE=e+f*j;forva=l:nifV(va)<=1.05&V(va)>=1.000;yj(va)=l;elseyj(va)=0;endendifyj(D*yj(2)*yj(3)*yj(4)*yj(5)*yj(6)*yj(7)*yj(8)*yj(9)==lifzy(1)*zy(2)*zy(3)*zy(4)*zy(5)*zy(6)*zy(7)*zy(8)*zy(9)==l;dispC發(fā)電機(jī)2出口電壓：)；disp(B2(2,4))；dispC發(fā)電機(jī)3出口電壓：上disp(B2(3,4))；dispC變壓器1接頭電壓：,);disp(Bl(l,5));dispC變壓器2接頭電壓：上disp(B1(2,5));dispC變壓器3接頭電壓：)disp(B1(3,5));dispC各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：)；disp(V);forp=l:nC(p)=0;forq=l:nC(p)=C(p)+conj(Y(p,q))*conj(E(q));endS(p)=E(p)*C(p);endfori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2)；q=Bl(i,l);endSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./(Bl(i,3)*Bl(i,5))));endfori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./(Bl(i,3)*Bl(i,5))));endZDS=0;fori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2)；q=Bl(i,l);endDS(i)=Si(p,q)+Sj(q,p);ZDS=DS(i)+ZDS;endmm=mm+l;ZDSS(mm)=sqrt(real(ZDS)A2+imag(ZDS)A2);dispC******************************************************************，＞fori=l:nZB224(mm)=B2(2,4)；ZB234(mm)=B2(3,4)；ZB115(mm)=Bl(l,5);ZB125(mm)=B1(2,5);ZB135(mm)=B1(3,5);endendendendendendendenddisp('******************************************************************)；'dispC滿足各點(diǎn)電壓要求的數(shù)據(jù)組共有：)；disp(mm);LOL=min(ZDSS);fordota=l:mmifZDSS(dota)==LOLmm=dota;endenddispC當(dāng)加入線路損耗時(shí)，得到其中線路損耗最小的一組組號(hào)及關(guān)鍵數(shù)據(jù)為：);disp(mm);disp('發(fā)電機(jī)2出口電壓：);disp(ZB224(mm));dispC發(fā)電機(jī)3出口電壓：');disp(ZB234(mm));dispC變壓器1接頭電壓：1);disp(ZB115(mm));dispf變壓器2接頭電壓：1);disp(ZB125(mm));dispC變壓器3接頭電壓：，)；disp(ZBl15(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用此組數(shù)據(jù)算的的潮流結(jié)果為：上disp('節(jié)點(diǎn)數(shù)：');disp(n);dispC節(jié)點(diǎn)支路數(shù)：);disp(nl);disp('平衡節(jié)點(diǎn)：上disp(isb);disp('精度：)；disp(pr);disp('Bl=');disp(Bl);disp('B2=');disp(B2);disp('X=');disp(X);Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);Sl=zeros(nl);dispC節(jié)點(diǎn)數(shù)：上disp(n);dispC節(jié)點(diǎn)支路數(shù)：,);disp(nl);disp('平衡節(jié)點(diǎn)：上disp(isb);disp('精度：上disp(pr);mm=O;X=[1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0;9,0];fori=l:nifX(i,2)~=0;p=X(i,l);Y(p,p)=l./X(i,2);endendfori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endY(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;endG=real(Y);B=imag(Y);fori=l:ne(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=B2(i,4)；endfori=l:nS(i)=B2(i,l)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICT1=O;IT2=1;N0=2*n;N=N0+l;a=0;whileIT2~=0IT2=0;a=a+l;fori=l:nifi~=isbC(i)=0;D(i)=0;forjl=l:nC(i)=C(i)+G(i,jl)*e(jl)-B(i,jl)*f(jl);D(i)=D(i)+G(i,jl)*f01)+B(i,jl)*e(jl);endPl=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)~=3DP=P(i)-Pl;DQ=Q(i)-Ql;forjl=l:nifjl~=isb&jl~=iXl=-G(ijl)*e(i)-B(ijl)*f(i)；X2=B(ijl)*e(i)-G(ijl)*f(i);X3=X2;X4=-X1;p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;elseifj1==i&j1~=isbXl=-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(i,i)*f(i);X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);p=2*i-l;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;endendelseDP=P(i)-Pl;DV=V(i)A2-V2;forjl=l:nifjl~=isb&jl~=iXl=-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+1;J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X6;J(m,q)=X2;elseifj1==i&j1~=isbXl=-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-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;endendendendendfork=3:N0kl=k+l;Nl=N;fork2=kl:NlJ(k,k2)=J(k,k2),/J(k,k);endJ(k,k)=l;ifk~=3k4=k-l;fork3=3:k4fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endifk==NO,break;endfork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endelsefork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendendfork=3:2:N0-lL=(k+1).⑵e(L)=e(L)-J(k,N);kl=k+l;f(L)=f(L)-J(kl,N);endfork=3:N0DET=abs(J(k,N));ifDET>=prIT2=IT2+1;endendICT2(a)=IT2;ICT1=ICT1+1;fork=l:ndy(k)=sqrt(e(k)A2+f(k)A2);endfori=l:nDy(ICTl,i)=dy(i);endend%用高斯消去法解"w=-J*V"disp('迭代次數(shù)上disp(ICTl);dispC沒有達(dá)到精度要求的個(gè)數(shù)工disp(ICT2);fork=l:nV(k)=sqrt(e(k)A2+f(k)A2);O(k)=atan(f(k)./e(k))*18O./pi;endE=e+f*j;dispC各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E為(節(jié)點(diǎn)號(hào)從小到大排列)：1disp(E);disp('各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：,)；disp(V);dispC各節(jié)點(diǎn)的電壓角。為(節(jié)點(diǎn)號(hào)從小到大排列)：)；disp(O);forp=l:nC(p)=0；forq=l:nC(p)=C(p)+conj(Y(p,q))*conj(E(q));endS(p)=E(p)*C(p);enddisp。各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：,);disp(S);dispC各條支路的首段功率Si為(順序同您輸入Bl時(shí)一樣)：)；fori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./(Bl(i,3)*Bl(i,5))));disp(Si(p,q));enddispC各條支路的末段功率Sj為(順序同您輸入Bl時(shí)一樣)：)；fori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./(Bl(i,3)*Bl(i,5))));disp(Sj(q,p));enddispC各條支路的功率損耗DS為(順序同您輸入B1時(shí)一樣)：,);fori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endDS(i)=Si(p,q)+Sj(q,p);disp(DS(i));endfori=l:ICTlCs(i)=i;enddispC以下是各節(jié)點(diǎn)每次迭代后的電壓值(如圖所示),);plot(Cs,Dy),xlabel。迭代次數(shù))ylabelC電壓)title。電壓迭代次數(shù)曲線,)；最后是課程設(shè)計(jì)報(bào)告二刃理N火彎課程設(shè)計(jì)課程名稱： 電力系統(tǒng)分析 設(shè)計(jì)題目：基于Matlab計(jì)算程序的電力系統(tǒng)運(yùn)行分析學(xué) 院： 電力工程學(xué)院 專業(yè)： 電氣工程自動(dòng)化 年 級(jí)： 2009 學(xué)生姓名易軍（你可以稱我易大神，可惜dota、lol都是菜鳥頂級(jí)水平）指導(dǎo)教師： 宋琪、謝靜、單節(jié)杉 日期： 2013T-13 教務(wù)處制電力系統(tǒng)潮流計(jì)算是電力系統(tǒng)分析中的一種最基本的計(jì)算，是對(duì)復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運(yùn)行狀態(tài)的計(jì)算。潮流計(jì)算的目標(biāo)是求取電力系統(tǒng)在給定運(yùn)行狀態(tài)的計(jì)算。即節(jié)點(diǎn)電壓和功率分布，用以檢查系統(tǒng)各元件是否過負(fù)荷.各點(diǎn)電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對(duì)現(xiàn)有用力系統(tǒng)的運(yùn)行和擴(kuò)建，對(duì)新的電力系統(tǒng)進(jìn)行規(guī)劃設(shè)計(jì)以及對(duì)用力系統(tǒng)進(jìn)行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計(jì)算為基礎(chǔ)。潮流計(jì)算結(jié)果可用如電力系統(tǒng)穩(wěn)態(tài)研究，安全估計(jì)或最優(yōu)潮流等對(duì)潮流計(jì)算的模型和方法有直接影響。在電力系統(tǒng)中可能發(fā)生的各種故障中，危害最大且發(fā)生概率較高的首推短路故障。產(chǎn)生短路故障的主要原因是電力設(shè)備絕緣損壞。短路故障分為三相短路、兩相短路、單相接地短路及兩相接地短路。其中三相短路時(shí)三相電流仍然對(duì)稱，其余三類短路統(tǒng)成為不對(duì)稱短路。短路故障大多數(shù)發(fā)生在架空輸電線路。電力系統(tǒng)設(shè)計(jì)與運(yùn)行時(shí)，要采取適當(dāng)?shù)拇胧┙档投搪饭收系陌l(fā)生概率。短路計(jì)算可以為設(shè)備的選擇提供原始數(shù)據(jù)。H標(biāo)電網(wǎng)接線圖

一、電網(wǎng)模型的建立設(shè)計(jì)中，選取Sb=100MVA,Ub=220KV,將所有支路的參數(shù)都折算到220KV電壓等級(jí)側(cè)本文利用matlab計(jì)算參數(shù)其源程序?yàn)椋篶lc;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=l00/220A2*0.4*L; %線路各段電抗(均為標(biāo)幺值，disp(Xl);fori=l:6R(i)=100/220A2*L(i)*r(i);enddispC線路電阻R=');disp(R);fori=l:6Z(i)=R(i)+j*Xl(i);Z0(i)=R(i)+3*j*Xl(i);enddispC線路正序、負(fù)序阻抗Z1=Z2=);disp(Z);fori=l:6Bl(i)=220A2/100*b(i)*L(i);B(i)=0+j*Bl(i);enddispC線路容納之半B/2=');disp(B);dispC線路零序阻抗Z0:?);disp(ZO);U0=[l1,10,6];SN=[240,220J20J;UN=[242,242,242];fori=l:3;XT(i)=j*U0(i)*242A2/SN(i)/220A2;XT0(i)=XT(i);end;dispC變壓器正、負(fù)序電抗值Xt：);disp(XT);dispC變壓器零序電抗值XtO：*);disp(XT0);取Ub=220KV,SB=100MVA)%加分號(hào)避免出現(xiàn)過程項(xiàng)，下同取Ub=220KV,SB=100MVA)%加分號(hào)避免出現(xiàn)過程項(xiàng)，下同%線路各段電阻%輸出R%線路各容納之半%輸出B%發(fā)電機(jī)暫態(tài)電抗有名值UGN=[16.5,18,13.8];fori=l:3K(i)=242/UGN(i);endfori=l:3Xdg(i)=j*l00/220A2*x(i)*K(i)*K(i);enddispC發(fā)電機(jī)暫態(tài)電抗Xdg=')；disp(Xdg);PL=[125,90,100];QL=[50,30,35];fori=l:3ZL(i)=i00/(PL(i)A2+QL(i)A2)*(PL(i)+j*QL(i));enddispC負(fù)荷節(jié)點(diǎn)阻抗Zl=');disp(ZL);%發(fā)電機(jī)暫態(tài)電抗%負(fù)荷節(jié)點(diǎn)阻抗計(jì)算結(jié)果為：線路電抗X=0.0860 0.0909 0.0810線路電阻!^=0.0101 0.0168 0.0160線路正序、負(fù)序阻抗Z1=Z2=Columns1through50.13470.03100.0101+0.0860iColumn60.0118+0.1008i線路容納之半B/2=Columns1through50+0.0896iColumn60+0.105H線路零序阻抗Z0:Columns1through50.0101+0.2579iColumn60.0118+0.3025i0.0168+0.0909i0+0.0783i0.0168+0.2727i變壓器正、負(fù)序電抗值Xt：0+0.0555i 0+0.0550i變壓器零序電抗值XtO：0+0.0555i 0+0.0550i發(fā)電機(jī)暫態(tài)電抗Xdg=0+0.0270i 0+0.0447i負(fù)荷節(jié)點(diǎn)阻抗Zl=0.6897+0.2759i 1.0000+O.3333i0.07270.00850.0160+0.08lOi0+0.0759i0.0160+0.2430i0+0.0605i0+0.0605i0+0.1152i0.8909+0.3118i0.10080.01180.0310+0.1347i0+0.1420i0.0310+0.404li0.0085+0.0727i0+0.0758i0.0085+0.2182i支數(shù)參數(shù)形成的矩陣：EH=Columns1through51.00004.0000 0+0.0555i 0 1.00002.00007.0000 0+0.0550i 0 1.00003.00009.0000 0+0.0605i 0 1.00004.00005.0000 0.0101+0.0860i 0+0.0896i 1.00004.00006.0000 0.0168+0.0909i 0+0.0783i 1.00005.00007.0000 0.0160+0.08lOi 0+0.0759i 1.00006.00009.0000 0.0310+0.1347i 0+0.1420i 1.00007.00008.0000 0.0085+0.0727i 0+0.0758i 1.00008.00009.0000 0.0118+0.1008i 0+0.1051i 1.0000Column6000000000支數(shù)參數(shù)形成的矩陣：B2=Columns1through52.0000+1.2395i0 1.0500 1.0000 01.8000+1.1155i0 1.0000 1.0000 01.0000+0.6197i0 1.0000 1.0000 000 1.0000 0 00 1.2500+0.5000i 1.0000 0 00 0.9000+0.3000i 1.0000 0 00 01.0000 0 00 1.0000+O.35OOi 1.0000 0 00 01.0000 0 0Column61.00003.00003.00002.00002.00002.00002.00002.00002.0000

系統(tǒng)等值電路圖的繪制根據(jù)以上計(jì)算結(jié)果，得到系統(tǒng)等值電路圖如下:0.8909+J0.31180.8909+J0.3118E2J0.0447jaosso0.0310+J0X)7270.0118+J0.1008J0.06051J04051J0.1051371工JO0758 J0Q75^1-1-1J0,1278E2J0.0447jaosso0.0310+J0X)7270.0118+J0.1008J0.06051J04051J0.1051371工JO0758 J0Q75^1-1-1J0,12780.0160*J0.0810.II-.J01708

0.03310+JO.134^J0.1152J0.0759TIJ0.0759TIJ0J420II-0.6897+J0.2759—iiJ0.1060a0101+J0.0860J0.08960.6897+J0.2759—iiJ0.1060a0101+J0.0860J0.0896HRII—JO07190.0168+J0.0909ja0783-HI—1.0000+J0^333J0.0555J0.0270ElJ0.0555J0.0270El二、潮流計(jì)算.系統(tǒng)參數(shù)的設(shè)置設(shè)計(jì)中要求所有結(jié)點(diǎn)電壓不得低于1.Op.U.,也不得高于1.05p.u.,若電壓不符合該條件，可采取下面的方法進(jìn)行調(diào)壓：改變發(fā)電機(jī)的機(jī)端電壓改變變壓器的變比（即改變分接頭）改變發(fā)電機(jī)的出力（4）在電壓不符合要求的結(jié)點(diǎn)處增加無功補(bǔ)償調(diào)壓方式應(yīng)屬于逆調(diào)壓。結(jié)點(diǎn)的分類：根據(jù)電力系統(tǒng)中各結(jié)點(diǎn)性質(zhì)的不同，將結(jié)點(diǎn)分為三類：PQ結(jié)點(diǎn)、PV結(jié)點(diǎn)和平衡結(jié)點(diǎn)，在潮流計(jì)算中，大部分結(jié)點(diǎn)屬于PQ結(jié)點(diǎn)，小部分結(jié)點(diǎn)屬于PV結(jié)點(diǎn)，一般只設(shè)一個(gè)平衡結(jié)點(diǎn)。對(duì)于平衡結(jié)點(diǎn)，給定其電壓的幅值和相位，整個(gè)系統(tǒng)的功率平衡由這一點(diǎn)承擔(dān)。本設(shè)計(jì)中，選1號(hào)節(jié)點(diǎn)為平衡節(jié)點(diǎn)；2、3號(hào)節(jié)點(diǎn)為P、U節(jié)點(diǎn)；4、5、6、7、8、9號(hào)結(jié)點(diǎn)為P、Q節(jié)點(diǎn)。采用PQ分解法進(jìn)行潮流的計(jì)算和分析。設(shè)計(jì)中，節(jié)點(diǎn)數(shù):n=9,支路數(shù):nl=9,平衡母線節(jié)點(diǎn)號(hào):isb=l,誤差精度:pr=O.OOOOl。由支路參數(shù)形成的矩陣：矩陣B1的每行是由卜列參數(shù)構(gòu)成的：1、某支路的首段號(hào)p；2、末端號(hào)q,且p<q；3、支路的阻抗（R+Xi）；4、支路的對(duì)地導(dǎo)納；5、支路的變比；6、折算到哪一側(cè)的標(biāo)志（支路首段p處于高壓側(cè)則輸入“1”，否則為"0”）；Bl=[l40.05546i010;270.055i010;390.0605i010;450.01+0.0859510.1792i10;460.0168+0.091i0.1566i10;570.016+0.081i0.1518i10;690.031+0.1347i0.284i10;780.0085+0.0727i0.1516i10;890.0119+0.1008i0.2102i10]矩陣B2的每行是由下列參數(shù)構(gòu)成的：1、節(jié)點(diǎn)所接發(fā)電機(jī)的功率SC；2、節(jié)點(diǎn)負(fù)荷功率S1;3、節(jié)點(diǎn)電壓初始值；4、PU節(jié)點(diǎn)電壓U的給定值：5、節(jié)點(diǎn)所接的無功補(bǔ)償設(shè)備容量；6、節(jié)點(diǎn)分類標(biāo)號(hào)（1—平衡節(jié)點(diǎn)；2—PQ節(jié)點(diǎn)；3—PU節(jié)點(diǎn)）本文采用的是自動(dòng)篩選法，利用matlab選出合理的數(shù)據(jù)組成，然后驗(yàn)證合理性。首先采用只調(diào)節(jié)發(fā)電機(jī)，結(jié)果沒有可選值。加上變壓器后通過限制初次迭代值和節(jié)點(diǎn)電壓值最終篩選出了結(jié)果。源程序如下clc;n=9;nl=9;isb=l;pr=O.0001;mm=0;Bl=[1.000004.00000+0.0555i01.00002.00007.00000+0.0550i01.000003.00009.00000+0.0605i01.000004.00005.00000.0101+0.0860i04-0.0896i1.000004.00006.00000.0168+0.0909i04-0.0783i1.000005.00007.00000.0160+0.0810i0+0.0759i1.000006.00009.00000.0310+0.1347i04-0.1420i1.000007.00008.00000.0085+0.0727i04-0.0758i1.000008.00009.00000.0118+0.1008i04-0.1051i1.00000]B2=[2.0000+1.2395i01.05001.000001.0000;1.8000+1.1155i01.00001.000003.0000;1.0000+0.6197i01.00001.000003.0000：001.0000002.0000;01.2500+0.5000i1.0000002.0000;00.9000+0.3000i1.0000002.0000;001.0000002.0000;01.0000+0.3500i1.0000002.0000;001.0000000000;];X=[l,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0;9,0];foraa=0:5B2(2,4)=l+0.01*aa;forab=0:5B2(3,4)=l+0.01*ab;forba=0:5Bl(l,5)=0.95+0.025*ba;forbb=O:5Bl(2,5)=0.95+0.025*bb;forbc=0:5Bl(3,5)=0.95+0.025*bc;Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);0二zeros(1,n);Sl=zeros(nl)生成初始矩陣Y為n*n的全0矩陣efV0均為Ixn的全0矩陣fori=l:nifX(i,2)~=0;P=X(i,1);Y(p,p)=l./X(i,2);endendfori=l:nlifBl(i,6)=0p=Bl(i,1);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,1);end(p,q)=Y(p,q)-l./(Bl(i,3)*Bl(i,5));(q,p)=Y(p,q);(q,q)=Y(q,q)+l./(Bl(i,3)*Bl(i,5廠2)+Bl(i,4)./2;(p,p)=Y(p,p)+l./Bl(i,3)+Bl(i,4)./2;endG=real(Y);B=imag(Y);fori=l:ne(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=B2(i,4)；%判斷X(i,2)不等于0%判斷是否有接地阻抗支路%不含變壓器的支路%含有變壓器的支路%計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣%電導(dǎo)%電納%電壓fori=l:nS(i)=B2(i,l)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICT1=O;IT2=1;N0=2*n;N=NO+1;a=0;whileIT2'=0IT2=0;a=a+l;fori=l:nifi"=isbCG)=0;D⑴=0;forjl=l:nC(i)=C(i)+G(i,jl)*e(jl)-B(i,jl)*f(jl);D(i)=D(i)+G(i,jl)*f(jl)+B(i,jl)*e(jl);endPl=C(i)*e(i)+f(i)*D(i);Ql=f(i)*C(i)-D(i)*e(i);V2二e(i廠2+f(i)C2;ifB2(i,6)DP二P⑴-Pl;DQ=Q(i)-Ql;forjl=l:nifjP=isb&jl"=iXl=-G(i,jl)*e(i)-B(i,jl)*f(i)；X2=B(i,jl)*e(i)-G(i,jl)*f(i);X3=X2;X4=-X1;p=2*i-1;q=2*jl-l;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;elseifjl==i&jl"=isbXl=-C(i)-G(i,i)*e(i)-B(i,i)*f⑴；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*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;endendelseDP=P(i)-Pl;DV=V(i)^2-V2;forjl=l:nifjT=isb&jr=iXl=-G(i,jl)*e(i)-B(i,jl)*f(i);X2=B(i,jl)*e(i)-G(i,jl)*f(i);X5=0;X6=0;P=2*i-1;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~=isbXl=-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-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X6;J(m,q)=X2;endendendendendfork=3:N0kl=k+l;Nl=N;fork2=kl:NlJ(k,k2)=J(k,k2)./J(k,k);endJ(k,k)=l;ifk、=3k4=k-l;fork3=3:k4fork2=kl:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endifk==N0,break;endfork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endelsefork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendendfork=3:2:N0-lL=(k+l)./2;e(L)=e(L)-J(k,N);kl=k+l;f(L)=f(L)-J(kl,N);endfork=3:N0DET=abs(J(k,N));ifDET>=prIT2=IT2+1;endendICT2(a)=IT2;ICT1=ICT1+1;fork=l:ndy(k)=sqrt(e(k)*2+f(k)2);endfori=l:nDy(ICTl,i)=dy(i);endforvb=l:nifDy(l,vb)<=1.05&Dy(l,vb)>=1.00;zy(vb)=l;elsezy(vb)=0;endendend%用高斯消去法解〃w=-J*V〃fork=l:nV(k)=sqrt(e(k)"2+f(k)2);0(k)=atan(f(k)./e(k))*180./pi;endE=e+f*j;forva=l:nifV(va)<=1.05&V(va)>=1.000;yj(va)=l;elseyj(va)=0;endendifyjd)*yj(2)*yj(3)*yj(4)*yj(5)*yj(6)*yj(7)*yj(8)*yj(9)==1ifzy(l)*zy(2)*zy(3)*zy(4)*zy(5)*zy(6)*zy(7)*zy(8)*zy(9)==1disp('發(fā)電機(jī)2出口電壓：')；disp(B2(2,4));disp('發(fā)電機(jī)3出口電壓：')；disp(B2(3,4));disp('變壓器1接頭電壓：')；disp(Bl(l,5));disp('變壓器2接頭電壓：')；disp(Bl(2,5));disp('變壓器3接頭電壓：')；disp(81(3,5));dispf各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：’)；disp(V);disp('#################################################################')；endendendendendendenddisp('滿足要求的數(shù)據(jù)組共有：’)；disp(mm)：dispC#################################################################')；運(yùn)行后共篩選出了41組數(shù)據(jù)，經(jīng)過改變電壓范圍最終確定了下面兩組數(shù)據(jù)。第一組數(shù)據(jù)：節(jié)點(diǎn)數(shù)： 9節(jié)點(diǎn)支路數(shù)：9平衡節(jié)點(diǎn)： 1精度：1.0000e-004Bl=Columns1through51.00004.00000+O.O555i01.02502.00007.00000+0.0550i01.02503.00009.00000+0.0605i00.97504.00005.00000.0101+0.0860i0+0.089611.00004.00006.00000.0168+0.0909i0+0.078311.00005.00007.00000.0160+0.0810i0+0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i0+0.0758i1.00008.00009.00000.0118+0.1008i0+0.1051i1.0000Column6000000000B2=Columns1through52.0000+1.2395i01.05001.000001.8000+1.1155i01.00001.030001.0000+0.6197i01.00001.03000001.00000001.2500+0.5000i1.00000000.9000+0.3000i1.000000001.00000001.0000+0.3500i1.000000001.000000Column61.00003.00003.00002.00002.000()2.00002.(XX)02.00002.0000迭代次數(shù)4沒有達(dá)到精度要求的個(gè)數(shù)16 16 1 0各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E為(節(jié)點(diǎn)號(hào)從小到大排列)：Columns1through51.0500 1.0176+0.1595i 1.0241+O.HOOi1.0402?0.0213i 1.0065-

Columns6through91.0062-0.0383i 1.0299+0.0617i 1.0072+0.0163i 1.0116+0.051Oi各節(jié)點(diǎn)的電壓大小V為（節(jié)點(diǎn)號(hào)從小到大排列）：Columns1through81.0500 1.0300 1.0300 1.0404 1.0069 1.0070 1.0317Column91.0129各節(jié)點(diǎn)的電壓角O為（節(jié)點(diǎn)號(hào)從小到大排列）：Columns1through80 8.9093 6.1298 -1.1722 -1.5602 -2.1784 3.4298Column92.8886各節(jié)點(diǎn)的功率S為（節(jié)點(diǎn)號(hào)從小到大排列）：Columns1through50.3928+0.6657i 1.8000+0.5248i 1.0000-0.1229i 0.0000+O.OOOOiColumns6through9-0.9000-0.3000i-0.0000+O.OOOOi-1.0000-0.3500i 0各條支路的首段功率Si為（順序同您輸入Bl時(shí)?樣）：0.3928+0.6657i1.8000+0.5248i1.0000-0.1229i0.1284+0.3424i1.00730.9292-1.2500-0.5000i0.2644+0.2933i1.00730.9292-1.2500-0.5000i-0.6384+0.0600i0.6569+0.2433i-0.3472-0.0628i各條支路的末段功率Sj為（順序同您輸入Bl時(shí)一樣）:-0.3928-0.6356i-1.8000-0.3426i-1.0000+0.1808i-0.1268-0.4228i-0.2616-0.3600i1.1431+0.0992i0.6514-0.1484i-0.6528-0.2872i0.3486-0.0324i各條支路的功率損耗DS為（順序同您輸入Bl時(shí)一樣）:0.0000+0.030li0+0.1822i-0.0000+0.0579i0.0016-0.0805i0.0028-0.0667i0.0199+0.0220i0.0130-0.0884i0.0041-0.0438i以下是各節(jié)點(diǎn)每次迭代后的電壓值（如圖所示）第二組數(shù)據(jù)：#################################################################節(jié)點(diǎn)數(shù)：9節(jié)點(diǎn)支路數(shù)：9平衡節(jié)點(diǎn)：1精度：Bl=Columns1through51.00004.00000+O.O555i01.02502.00007.00000+O.O55Oi01.02503.00009.00000+0.0605100.97504.00005.00000.0101+0.0860i0+0.0896i1.00004.00006.00000.0168+0.0909i0+0.0783i1.00005.00007.00000.0160+0.08lOi0+0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i0+0.0758i1.00008.00009.00000.0118+0.100810+0.105Hl.(XX)OColumn6000000000Columns1through52.0000+1.2395i 01.05001.000001.8000+1.1155i01.00001.03000l.(XX)0+0.6197i01.00001.04000001.00000001.2500+0.5000i1.00000000.9000+0.3000i1.00000000l.OO(X)000l.OOOO+O.35OOi1.000000001.000000Column61.00003.00003.000()2.00002.00002.00002.00002.00002.0000迭代次數(shù)沒有達(dá)到精度要求的個(gè)數(shù)16 16I0各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E為（節(jié)點(diǎn)號(hào)從小到大排列）：Columns1through51.0500 1.0177+0.1589i 1.0343+0.1086i 1.0418-0.0212i 1.0083-0.0275iColumns6through91.0101-0.0385i 1.0317+0.0613i 1.0112+0.0158i 1.0186+0.0499i各節(jié)點(diǎn)的電壓大小V為（節(jié)點(diǎn)號(hào)從小到大排列）：Columns1through81.0500 1.0300 1.0400 1.0420 1.0086 1.0108 1.0335 1.0113Column91.0198各節(jié)點(diǎn)的電壓角O為（節(jié)點(diǎn)號(hào)從小到大排列）：Columns1through80 8.8729 5.9939 -1.1681 -1.5623 -2.1817 3.4029 0.8959Column92.8057各節(jié)點(diǎn)的功率S為（節(jié)點(diǎn)號(hào)從小到大排列）：Columns1through50.3921+0.6364i 1.8000+0.492li1.0000-0.0747i-0.0000+O.OOOOi-1.2500-0.5000iColumns6through9-0.9000-0.3000i 0.0000+O.OOOOi-1.0000-O.35OOi0各條支路的首段功率Si為（順序同您輸入Bl時(shí)一樣）：0.3921+0.6364i1.0000-0.0747i0.1298+0.3405i0.2623+0.2678i-1.1218-0.0786i-0.6402+0.0365i0.6584+0.2118i-0.3456-0.0928i各條支路的末段功率Sj為（順序同您輸入Bl時(shí)一樣）:-0.3921-0.6083i-1.8000-0.3116i-1.0000+0.1310i-0.1282-0.4214i-0.2598-0.3365i1.1416+0.0998i0.6530-0.1273i-0.6544-0.2572i0.3470-0.0037i各條支路的功率損耗DS為（順序同您輸入Bl時(shí)一樣）:0.0000+0.028li-0.0000+0.1805i0+0.0562i0.0016-0.0809i

0.0198+0.0212i0.0128-0.0908i0.0040-0.0454i0.0014-0.0965i以下是各節(jié)點(diǎn)每次迭代后的電壓值（如圖所示）Figure1,|n|Figure1FileEditViewInsertToolsDesktopWindowHelp ，口安。昌?0爭(zhēng)@I變口國I?亙"電壓迭代次數(shù)曲線迭代次數(shù)DS1=O.0428-0.0825iDS2=0.0422-0.0813i以上兩組數(shù)據(jù)均滿足要求，從減小線路損耗的角度看，選擇了第二組數(shù)據(jù)。潮流分布圖:bus-8Ol391140.3295!bus-8Ol391140.3295!-0.3900*Jft3270對(duì)運(yùn)行結(jié)果的分析線路有功功率總是從電壓相位超前的結(jié)點(diǎn)流向滯后的結(jié)點(diǎn)，無功功率的總是從電壓幅值高的結(jié)點(diǎn)流向電壓幅值低的結(jié)點(diǎn)。功率數(shù)值前的符號(hào)若為正表示和箭頭的方向一致，若為負(fù)則表示和箭頭的方向相反。在潮流分布圖中，3號(hào)發(fā)電機(jī)發(fā)出的無功功率為負(fù)值，說明3號(hào)發(fā)電機(jī)從系統(tǒng)吸收無功。由于電網(wǎng)運(yùn)行時(shí)負(fù)荷等的變化，引起系統(tǒng)中各節(jié)點(diǎn)電壓的變化。由于發(fā)電機(jī)容量有限，而電壓的調(diào)整又必須服從對(duì)電壓質(zhì)量的要求，調(diào)整的幅度不能太大，另一方面從電網(wǎng)運(yùn)行經(jīng)濟(jì)性和安全等方面考慮，網(wǎng)絡(luò)中的潮流往往需要控制。在實(shí)際的網(wǎng)絡(luò)潮流控制中主要采用：串聯(lián)電容（作用以容抗抵償線路感抗）；串聯(lián)電抗（作用在于限流）；附加串聯(lián)加壓器（作用在于產(chǎn)生一環(huán)流或強(qiáng)制循環(huán)功率，使強(qiáng)制循環(huán)功率與自然分布功率的疊加可達(dá)到理想值）等手段控制潮流。三、故障電流計(jì)算1.三相短路電流的計(jì)算利用結(jié)點(diǎn)阻抗矩陣和導(dǎo)納矩陣都可以計(jì)算短路電流，其算法有所不同。利用結(jié)點(diǎn)阻抗陣時(shí),只要形成了阻抗陣，計(jì)算網(wǎng)絡(luò)中任意一點(diǎn)的對(duì)稱短路電流和網(wǎng)絡(luò)中電流、電壓的分布非常方便，計(jì)算工作量小，但是，形成阻抗陣的工作量大，網(wǎng)絡(luò)變化時(shí)的修改也比較麻煩，而且結(jié)點(diǎn)阻抗矩陣是滿陣，需要計(jì)算機(jī)存儲(chǔ)量較大。利用結(jié)點(diǎn)導(dǎo)納矩陣計(jì)算短路電流，實(shí)質(zhì)是先用它計(jì)算出與指定的短路點(diǎn)有關(guān)的結(jié)點(diǎn)阻抗矩陣的相應(yīng)列的元素，然后利用公式/f=―2一,q=(/,(0)-zj=Ui~Uj(z,為連接節(jié)點(diǎn)ijz.+z * jij 4和j的支路阻抗)進(jìn)行短路電流的計(jì)算。然而，導(dǎo)納陣是對(duì)稱、稀疏陣，極易形成，且網(wǎng)絡(luò)結(jié)構(gòu)變化時(shí)也易于修改。NF=1n=9nl=9B=Columns1through51.00004.00000+0.0555i01.02502.00007.00000+0.0550i01.02503.00009.00000+0.0605i00.97504.00005.00000.01014-0.0860i0+0.0896i1.00004.00006.00000.0168+0.0909i0+0.0783i1.00005.00007.00000.0160+0.0810i04-0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i04-0.0758i1.00008.00009.00000.0118+0.1008i0+0.1051i1.0000Column600000000vo=1.05001.03001.04001.04201.00861.01081.03351.01131.0198D=7 07點(diǎn)短路時(shí)電流的標(biāo)幺值If=0.5504-4.636li各節(jié)點(diǎn)的電壓標(biāo)幺值U為（節(jié)點(diǎn)號(hào)從小到大排）：1.05001.03001.04001.04200.3594-0.0296i1.01080+1.3878e-017i-0.0222+O.OOOOi-0.0137+0.0000i各支路短路電流的標(biāo)幺值I為（順序同您輸入B時(shí)一樣）:0-0.602li-0.0000-18.7273i-0.0000-17.4223i1.2588-7.7890i0.0613-0.3319i0.4920-4.3402i1.6624-7.2232i0.0352-0.3012i-0.0097+0.0832i2.單相短路電流的計(jì)算請(qǐng)輸入短路點(diǎn)的數(shù)目：NF=1請(qǐng)輸入短路類型（單相接地短路為1,兩相短路為2,兩相短路接地為3）:Lf=l正序、負(fù)序、零序節(jié)點(diǎn)數(shù)為：nl=9n2=9n0=9正序、負(fù)序、零序支路數(shù)為:nil=9nl2=9nlO=9短路節(jié)點(diǎn)號(hào):f=7縱橫故障標(biāo)志：If=0接地阻抗:zf=0接地阻抗:zg=0正序支路參數(shù):Bl=Columns1through51.00002.00004.00007.00000+0.0555i0+0.0550i001.02501.02503.00009.00000+0.0605i00.97504.00005.00000.0101+0.0860i0+0.0896i1.00004.00006.00000.0168+0.0909i0+0.0783i1.00005.00007.00000.0160+0.08lOi0+0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i0+0.0758i1.00008.00009.00000.0118+0.1008i0+0.1051i1.0000Column6000000000負(fù)序支路參數(shù)B2=Columns1through51.00004.00000+0.0555i01.02502.00007.00000+0.0550i01.02503.00009.00000+0.0605i00.97504.00005.00000.0101+0.0860i0+0.0896i1.00004.00006.00000.0168+0.0909i0+0.0783i1.00005.00007.00000.0160+0.08lOi0+0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i0+0.0758i1.0000

8.00009.00000.0118+0.1008i0+0.105H1.0000Column6000000000B0=Columns1through51.00004.0