昆明理工大學計算機輔助分析課程設計及matlab自動實現(xiàn)潮流調節(jié)與計算

1、重要聲明： 1、本文只是提供一定參考，有錯勿怪。2.、本文采用的潮流調節(jié)及計算只針對本實驗，其中有很多固化的數(shù)據(jù)與參數(shù)，以及特殊的公式、請根據(jù)要求的不同進行一定的調整（潮流計算的程序適用任何情況、不同的只是系統(tǒng)參數(shù)）。3、本文潮流調節(jié)為自動調節(jié)生成，調節(jié)變量有發(fā)電機2、3電壓，三臺變壓器分接頭，沒有進行無功補償，請根據(jù)具體情況加入或刪減變量。4、自動調節(jié)潮流程序中，控制變量為各點電壓首次迭代值、最終各點電壓。最終選擇標準為網損最低的一組數(shù)據(jù)。5、本實驗報告為草版，有錯誤，但懶得改。6、鑒于個人所有文檔免費得原則，本文為免費文檔，請盡量下載。7、最后！都免費了，要求不要太高。要求及原始數(shù)據(jù)電力系

2、統(tǒng)潮流計算是電力系統(tǒng)分析中的一種最基本的計算，是對復雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。潮流計算的目標是求取電力系統(tǒng)在給定運行狀態(tài)的計算。即節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷.各點電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對現(xiàn)有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設計以及對電力系統(tǒng)進行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計算為基礎。潮流計算結果可用如電力系統(tǒng)穩(wěn)態(tài)研究，安全估計或最優(yōu)潮流等對潮流計算的模型和方法有直接影響。在電力系統(tǒng)中可能發(fā)生的各種故障中，危害最大且發(fā)生概率較高的首推短路故障。產生短路故障的主要原因是電力設備絕緣損壞。短路故障分為三相短路、

3、兩相短路、單相接地短路及兩相接地短路。其中三相短路時三相電流仍然對稱，其余三類短路統(tǒng)成為不對稱短路。短路故障大多數(shù)發(fā)生在架空輸電線路。電力系統(tǒng)設計與運行時，要采取適當?shù)拇胧┙档投搪饭收系陌l(fā)生概率。短路計算可以為設備的選擇提供原始數(shù)據(jù)。一、課程設計的目的與任務:該課程設計為電氣工程及其自動化專業(yè)的主要專業(yè)課設計之一，是學習電力系統(tǒng)分析基礎和電力系統(tǒng)計算機輔助分析課程后的一個重要的實踐性教學環(huán)節(jié)。其目的在于鞏固和加深對電力系統(tǒng)潮流和短路電流計算基本原理的理解，學習和掌握應用計算機進行電力系統(tǒng)設計和計算的方法，培養(yǎng)學生獨立分析和解決問題的能力。二、課程設計的基本要求:掌握電力系統(tǒng)等值模型和參數(shù)計算，

4、以及潮流和短路計算的基本原理，學會應用計算機計算系統(tǒng)潮流分布和短路電流的方法。三、課程設計選題原則:該課程設計是根據(jù)電力系統(tǒng)分析基礎和電力系統(tǒng)計算機輔助分析課程內容，結合實際工程和科研的電力系統(tǒng)網絡進行系統(tǒng)的潮流和短路電流計算。四、課程設計內容:1、課程設計地點：機房和教室2、設計內容（1）電力系統(tǒng)潮流計算（2）電力系統(tǒng)短路電流計算具體內容包括：1）系統(tǒng)等值模型和參數(shù)計算。2）利用現(xiàn)有的潮流和短路電流計算程序進行計算，或編制電力系統(tǒng)潮流和短路計算程序。3）計算結果分析。4）編寫說明書。3、時間安排1、系統(tǒng)等值模型和參數(shù)計算（2天）2、讀懂現(xiàn)有的算法程序或視需要編制算法程序并進行計算（4天）3、

5、計算結果整理和分析（1天）4、編寫說明書（2天）5、成績考核（1天）五、課程設計主要參考資料：、李光琦.電力系統(tǒng)暫態(tài)分析.北京：水利電力出版社.第二版、陳珩.電力系統(tǒng)穩(wěn)態(tài)分析.北京：水利電力出版社.第二版、孟祥萍.基于的電力系統(tǒng)設計與計算.長春：吉林人民出版社、孟祥萍.電力系統(tǒng)分析.北京：高等教育出版社、電力系統(tǒng)潮流和短路計算課程設計指導書.自編六、課程設計成績考核：根據(jù)學生的平時表現(xiàn)，說明書，以及設計完成后的口試成績等進行綜合評定。課程結束時，要求學生有課程設計說明書。課程設計成績以平時成績（占50%）、設計說明書（占50%）進行綜合評定。原始數(shù)據(jù)與思考題目標電網單線圖系統(tǒng)參數(shù)表1. 線路參

6、數(shù)表線路編號線路型號線路長度（km）線路電阻/km線路正序電抗/km線路容納之半S/km4-5LGJ-240/301040.0470.41.78×4-6LGJ-120/701100.0741.47×5-7LGJ-120/25980.0791.60×6-9LGJ-95/551630.0921.80×7-8LGJ-240/30880.0471.78×8-9LGJ-240/301220.0471.78×說明：線路零序電抗為正序電抗3倍。表2. 變壓器參數(shù)表線路編號變壓器型號變壓器變比（kV）短路電壓百分數(shù)(%)2-7SSPL-2200002

7、42±3×2.5%/20103-9SSPL-120000242±3×2.5%/1561-4SSPL-240000242±3×2.5%/17.511說明：變壓器零序電抗與正序電抗相等，且均為/Y0接法。表3. 發(fā)電機參數(shù)表1發(fā)電機額定功率MW額定電壓kV額定功率因數(shù)120016.50.852180180.85310013.80.85 表4. 發(fā)電機參數(shù)表2發(fā)電機母線名（S）（S）（S）110.14600.06088.960.09690.0969047.28220.89580.11986.000.86450.19690.53512.803

8、31.31250.18138.591.25780.25000.6006.02表5. 負荷數(shù)據(jù)表節(jié)點號有功負荷（MW）無功負荷（MVA）51255069030810035思考題1、 牛頓拉夫遜法與PQ分解法有哪些聯(lián)系？有哪些區(qū)別？二者的計算性能如何？2、 為什么在用計算機對某網絡初次進行潮流計算時往往是要調潮流，而并非任何情況下只一次送入初始值算出結果就行呢？要考慮什么條件？各變量是如何劃分的？哪些可調？哪些不可調？3、 發(fā)電機節(jié)點的注入無功為負值說明了什么？4、 線路有功潮流最有可能的流向是？線路無功潮流最有可能的流向是？5、 負荷功率因數(shù)對系統(tǒng)潮流有什么影響？6、 潮流控制的主要手段有哪些？

9、7、 如何降低系統(tǒng)網損？8、 如何校驗主設備的短路容量？9、 利用節(jié)點阻抗矩陣計算短路電流和利用導納矩陣計算短路電流的算法有何異同？10、減小短路電流的措施有哪些？設計主要內容：1、完成調試潮流計算與短路計算程序調試；計算程序以電力系統(tǒng)計算機輔助分析所用教材中的例程為基礎，要求掌握算法原理，建議進行適當改寫以利于本次設計應用。2、根據(jù)設計資料準備計算數(shù)據(jù)；整理原始數(shù)據(jù)，使其符合計算程序的輸入條件。要求比較采用變壓器型模型前后兩種情況的差異。3、為指定電網安排合理的穩(wěn)態(tài)運行方式；合理設置目標電網中的控制變量，以1號機為平衡機，將潮流分布調整為可行的運行方式（要求所有節(jié)點的電壓不得低于1.0 p.

10、u.，也不能高于1.05p.u.）繪制潮流分布圖。分析該運行方式的特點，提出進一步降低網損的改進的措施，并利用計算結果簡單證明之。具體內容參考設計資料所附思考題。4、列表分析節(jié)點7發(fā)生各種短路故障時的短路電流，以及運行方式對故障電流的影響；5、完成設計說明書。要求結構清晰，即從電網模型建立，到初態(tài)潮流分布計算，到問題分析，再到調整后的方案計算。故障電流計算也類似。不必逐一回答所有思考題。參數(shù)計算、潮流調節(jié)、與潮流計算matlab源程序%參數(shù)計算的clc; clear;L=104,110,98,163,88,122;r=0.047,0.074,0.079,0.092,0.047,0.047;b=

11、1.78e-006,1.47e-006,1.60e-006,1.80e-006,1.78e-006,1.78e-006;disp('線路電抗X=');X1=100/2202*0.4*L; %線路各段電抗（均為標幺值，取Ub=220KV,SB=100MVA，下同）disp(X1);for i=1:6 R(i)=100/2202*L(i)*r(i); %加分號避免出現(xiàn)過程項，下同end %線路各段電阻disp('線路電阻R='); disp(R); %輸出R for i=1:6Z(i)=R(i)+j*X1(i);Z0(i)=R(i)+3*j*X1(i); enddi

12、sp('線路正序、負序阻抗Z1=Z2='); disp(Z); for i=1:6 Bl(i)=2202/100*b(i)*L(i); B(i)=0+j*Bl(i); end %線路各容納之半disp('線路容納之半B/2='); disp(B); %輸出Bdisp('線路零序阻抗Z0:');disp(Z0);U0=11,10,6;SN=240,220,120;UN=242,242,242;for i=1:3; XT(i)=j*U0(i)*2422/SN(i)/2202; XT0(i)=XT(i);end;disp('變壓器正、負序電抗值

13、Xt：');disp(XT);disp('變壓器零序電抗值Xt0：');disp(XT0);x=0.0608,0.1198,0.1813; %發(fā)電機暫態(tài)電抗有名值UGN=16.5,18,13.8;for i=1:3 K(i)=242/UGN(i);endfor i=1:3 Xdg(i)=j*100/2202*x(i)*K(i)*K(i);end %發(fā)電機暫態(tài)電抗disp('發(fā)電機暫態(tài)電抗Xdg='); disp(Xdg); %輸出XdgPL=125,90,100;QL=50,30,35;for i=1:3 ZL(i)=100/(PL(i)2+QL(i)2

14、)*(PL(i)+j*QL(i);end %負荷節(jié)點阻抗disp('負荷節(jié)點阻抗Zl='); disp(ZL); %輸出ZL%支數(shù)參數(shù)形成的矩陣B1、支數(shù)參數(shù)形成的矩陣B2、節(jié)點號及其對地阻抗形成的矩陣X的形成B1=1,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(2)+ X1(2)*j,B(2),1,0;5,7,R(3)+X1(3)*j,B(3),1,0;6,9,R(4)+X1(4)*j,B(4),1,0;7,8,R(5)+X1(5)*j,B(5),1,0;8,9,R

15、(6)+X1(6)*j,B(6),1,0;disp('支數(shù)參數(shù)形成的矩陣：B1='); disp(B1); %輸出支數(shù)參數(shù)形成的矩陣B1P=200,180,100;W=(P/100)+(P/0.85)*sqrt(1-0.85*0.85)/100)*j;B2=W(1),0,1.05,1,0,1; W(2),0,1,1,0,3; W(3),0,1,1,0,3;0,0,1,0,0,2;0, 1.25+0.5i,1,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;disp('支數(shù)參數(shù)形成的矩

16、陣：B2='); disp(B2); %輸出支數(shù)參數(shù)形成的矩陣B2%潮流自動調節(jié)的與自動計算的n=9;nl=9;isb=1;pr=0.0001;mm=0; X=1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0; 9,0;disp('*');disp('只考慮各點電壓水平，得到的滿足要求的數(shù)據(jù)如下：');for aa=0:5 B2(2,4)=1+0.01*aa;for ab=0:5B2(3,4)=1+0.01*ab;for ba=0:5B1(1,5)=0.95+0.025*ba;for bb=0:5B1(2,5)=0.95+0.025*bb;

17、for bc=0:5B1(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);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 B1(i,6)=0 p=B1(i,1); q=B1(i,2); else p=B1(i,2); q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5); Y(q,p)=Y(p,q); Y(q,q)

18、=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2; Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;endG=real(Y);B=imag(Y); for i=1:n e(i)=real(B2(i,3); f(i)=imag(B2(i,3); V(i)=B2(i,4); endfor i=1:n S(i)=B2(i,1)-B2(i,2); B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;while IT2 =0 IT2=0;a=a+1;

19、 for i=1:n if i=isb C(i)=0; D(i)=0; for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1); D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1); end P1=C(i)*e(i)+f(i)*D(i); Q1=f(i)*C(i)-D(i)*e(i); V2=e(i)2+f(i)2; if B2(i,6)=3 DP=P(i)-P1; DQ=Q(i)-Q1; for j1=1:n if j1=isb & j1=i X1=-G(i,j1)*e(i)-B(i,j1)*f(i); X2=B(i,

20、j1)*e(i)-G(i,j1)*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; elseif j1=i & j1=isb X1=-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-1;q=2

21、*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; end end else DP=P(i)-P1; DV=V(i)2-V2; for j1=1:n if j1=isb&j1=i X1=-G(i,j1)*e(i)-B(i,j1)*f(i); X2=B(i,j1)*e(i)-G(i,j1)*f(i); X5=0; X6=0; 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

22、;J(m,q)=X2; elseif j1=i&j1=isb X1=-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; end end end endendfor k=3:N0 k1=k+1;N1=N; for k2=k1:N1 J(k,k2)=J(k,k2)./J(k,k); en

23、d J(k,k)=1; if k=3 k4=k-1; for k3=3:k4 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; end if k=N0,break;end for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; endelse for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; end endendfor k=

24、3:2:N0-1 L=(k+1)./2; e(L)=e(L)-J(k,N); k1=k+1; f(L)=f(L)-J(k1,N);endfor k=3:N0 DET=abs(J(k,N); if DET>=pr IT2=IT2+1; endendICT2(a)=IT2;ICT1=ICT1+1;for k=1:n dy(k)=sqrt(e(k)2+f(k)2);end for i=1:n Dy(ICT1,i)=dy(i); end for vb=1:n if Dy(1,vb)<=1.05&Dy(1,vb)>=1.00; zy(vb)=1; else zy(vb)=0;

25、endendend %用高斯消去法解"w=-J*V"for k=1:n V(k)=sqrt(e(k)2+f(k)2); O(k)=atan(f(k)./e(k)*180./pi;endE=e+f*j;for va=1:n if V(va)<=1.05&V(va)>=1.000; yj(va)=1; else yj(va)=0; endendif yj(1)*yj(2)*yj(3)*yj(4)*yj(5)*yj(6)*yj(7)*yj(8)*yj(9)=1 if zy(1)*zy(2)*zy(3)*zy(4)*zy(5)*zy(6)*zy(7)*zy(8)

26、*zy(9)=1;disp('發(fā)電機2出口電壓：');disp(B2(2,4);disp('發(fā)電機3出口電壓：');disp(B2(3,4);disp('變壓器1接頭電壓：');disp(B1(1,5);disp('變壓器2接頭電壓：');disp(B1(2,5);disp('變壓器3接頭電壓：');disp(B1(3,5);disp('各節(jié)點的電壓大小V為(節(jié)點號從小到大排列)：');disp(V);for p=1:n C(p)=0; for q=1:n C(p)=C(p)+conj(Y(p,q)

27、*conj(E(q); end S(p)=E(p)*C(p);endfor i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); endSi(p,q)=E(p)*(conj(E(p)*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5)-conj(E(q)*conj(1./(B1(i,3)*B1(i,5);end for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); endSj(q,p)=E(q)*(conj(

28、E(q)*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5)-conj(E(p)*conj(1./(B1(i,3)*B1(i,5); endZDS=0;for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end DS(i)=Si(p,q)+Sj(q,p); ZDS=DS(i)+ZDS;endmm=mm+1;ZDSS(mm)=sqrt(real(ZDS)2+imag(ZDS)2);disp('*');for i=1:nZB224(mm)=B2(2,4);ZB234(mm

29、)=B2(3,4);ZB115(mm)=B1(1,5);ZB125(mm)=B1(2,5);ZB135(mm)=B1(3,5);endendendendendendendenddisp('*');disp('滿足各點電壓要求的數(shù)據(jù)組共有：');disp(mm);LOL=min(ZDSS);for dota=1:mm if ZDSS(dota)=LOL mm=dota; endenddisp('當加入線路損耗時，得到其中線路損耗最小的一組組號及關鍵數(shù)據(jù)為：');disp(mm);disp('發(fā)電機2出口電壓：');disp(ZB22

30、4(mm);disp('發(fā)電機3出口電壓：');disp(ZB234(mm);disp('變壓器1接頭電壓：');disp(ZB115(mm);disp('變壓器2接頭電壓：');disp(ZB125(mm);disp('變壓器3接頭電壓：');disp(ZB115(mm);B2(2,4)=ZB224(mm);B2(3,4)=ZB234(mm);B1(1,5)=ZB115(mm);B1(2,5)=ZB125(mm);B1(3,5)=ZB135(mm);disp('用此組數(shù)據(jù)算的的潮流結果為：');disp('

31、;節(jié)點數(shù)：');disp(n);disp('節(jié)點支路數(shù)：');disp(nl);disp('平衡節(jié)點：');disp(isb);disp('精度：');disp(pr);disp('B1='); disp(B1); 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);S1=zeros(nl);disp('節(jié)點數(shù)：');d

32、isp(n);disp('節(jié)點支路數(shù)：');disp(nl);disp('平衡節(jié)點：');disp(isb);disp('精度：');disp(pr);mm=0; X=1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0; 9,0;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 B1(i,6)=0 p=B1(i,1); q=B1(i,2); else p=B1(i,2); q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B

33、1(i,3)*B1(i,5); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2; Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;endG=real(Y);B=imag(Y); for i=1:n e(i)=real(B2(i,3); f(i)=imag(B2(i,3); V(i)=B2(i,4); endfor i=1:n S(i)=B2(i,1)-B2(i,2); B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICT1=0;IT2=1;N0=2*n

34、;N=N0+1;a=0;while IT2 =0 IT2=0;a=a+1; for i=1:n if i=isb C(i)=0; D(i)=0; for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1); D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1); end P1=C(i)*e(i)+f(i)*D(i); Q1=f(i)*C(i)-D(i)*e(i); V2=e(i)2+f(i)2; if B2(i,6)=3 DP=P(i)-P1; DQ=Q(i)-Q1; for j1=1:n if j1=isb & j1=i X

35、1=-G(i,j1)*e(i)-B(i,j1)*f(i); X2=B(i,j1)*e(i)-G(i,j1)*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; elseif j1=i & j1=isb X1=-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)

36、+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)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2; end end else DP=P(i)-P1; DV=V(i)2-V2; for j1=1:n if j1=isb&j1=i X1=-G(i,j1)*e(i)-B(i,j1)*f(i); X2=B(i,j1)*e(i)-G(i,j1)*f(i); X5=0; X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1

37、; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2; elseif j1=i&j1=isb X1=-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; end end end endendfor k=3:N0 k1=k+1;N1=N; for

38、 k2=k1:N1 J(k,k2)=J(k,k2)./J(k,k); end J(k,k)=1; if k=3 k4=k-1; for k3=3:k4 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; end if k=N0,break;end for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; endelse for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k

39、,k2); end J(k3,k)=0; end endendfor k=3:2:N0-1 L=(k+1)./2; e(L)=e(L)-J(k,N); k1=k+1; f(L)=f(L)-J(k1,N);endfor k=3:N0 DET=abs(J(k,N); if DET>=pr IT2=IT2+1; endendICT2(a)=IT2;ICT1=ICT1+1;for k=1:n dy(k)=sqrt(e(k)2+f(k)2);end for i=1:n Dy(ICT1,i)=dy(i); endend %用高斯消去法解"w=-J*V"disp('迭代次數(shù)

40、');disp(ICT1);disp('沒有達到精度要求的個數(shù)');disp(ICT2);for k=1:n V(k)=sqrt(e(k)2+f(k)2); O(k)=atan(f(k)./e(k)*180./pi;endE=e+f*j;disp('各節(jié)點的實際電壓標幺值E為(節(jié)點號從小到大排列)：');disp(E);disp('各節(jié)點的電壓大小V為(節(jié)點號從小到大排列)：');disp(V);disp('各節(jié)點的電壓角O為(節(jié)點號從小到大排列)：');disp(O);for p=1:n C(p)=0; for q=1:

41、n C(p)=C(p)+conj(Y(p,q)*conj(E(q); end S(p)=E(p)*C(p);enddisp('各節(jié)點的功率S為(節(jié)點號從小到大排列)：');disp(S);disp('各條支路的首段功率 Si為(順序同您輸入B1時一樣)：');for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); endSi(p,q)=E(p)*(conj(E(p)*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5)-conj(E(q)*conj(1./(B

42、1(i,3)*B1(i,5); disp(Si(p,q);enddisp('各條支路的末段功率 Sj為(順序同您輸入B1時一樣)：'); for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); endSj(q,p)=E(q)*(conj(E(q)*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5)-conj(E(p)*conj(1./(B1(i,3)*B1(i,5); disp(Sj(q,p);end disp('各條支路的功率損耗DS為(順序同您輸入B1時一樣)：');for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end DS(i)=Si(p,q)+Sj(q,p); disp(DS(i);endfor i=1:ICT1 Cs(i)=i;enddisp('以下是各節(jié)點每次迭代后的電壓值（如圖所示）');plot(Cs,Dy),xlabel('迭代次數(shù)'),ylabel('電壓'),title('電壓迭代次數(shù)曲線