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

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

2、系統(tǒng)潮流計算是電力系統(tǒng)分析中的一種最基本的計算，是對復雜電力系統(tǒng)正 常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。潮流計算的目標是求取電力系統(tǒng)在給定運行狀 態(tài)的計算。即節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷.各點電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對現有電力系統(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)成為不對稱短路。短路故障大多數發(fā)生在架空輸電線路。電力系統(tǒng)設計與運行時， 要采取適當的措施降低短路故障的發(fā)生概率。短路計算可以為設備的選擇提供原始數 據。一、課程設計的目的與任務：該課程設計為電氣工程及其自動化專業(yè)的主要專業(yè)課設計之一， 是學習電力系統(tǒng) 分析基礎和電力系統(tǒng)計算機輔助分析課程后的一個重要的實踐性教學環(huán)節(jié)。 其目的在 于鞏固和加深對電力系統(tǒng)潮流和短路電流計算基本原理的理解， 學習和掌握應用計算 機進行電力系統(tǒng)設計和計算的方法，培養(yǎng)學生獨立分析和解決問題的能力。二、課程設計的基本

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

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

6、進行綜合評定。原始數據與思考題目標電網單線圖系統(tǒng)參數表1.線路參數表線路編號線路型號線路長度（km）線路電阻 Q/km線路正序電抗Q/kn）線路容納之半s/km4-5LGJ-240/301040.0470.4/ c -61.78 X 104-6LGJ-120/701100.074/ C -61.47 X 105-7LGJ-120/25980.0791.60X 10"6-9LGJ-95/551630.0921.80X 10“7-8LGJ-240/30880.047/ C -61.78 X 108-9LGJ-240/301220.047,-61.78 X 10說明：線路零序電抗為正序電抗

7、 3倍。表2.變壓器參數表線路編號變壓器型號變壓器變比（kV）短路電壓白分數（）2-7SSPL-220000242±3X 2.5%/20103-9SSPL-120000242±3X 2.5%/1561-4SSPL-240000242±3X2.5%/17.511說明：變壓器零序電抗與正序電抗相等，且均為 A/Y。接法。表3.發(fā)電機參數表1發(fā)電機額定功率MW額定電壓kV額定功率因數120016.50.852180180.85310013.80.85表4.發(fā)電機參數表2發(fā)電機母線名XdQ-, XdQTd0（S）XqQ一'XqQ一'Tq0 (S)Tj (S

8、)110.14600.06088.960.09690.0969047.28220.89580.11986.000.86450.19690.53512.80331.31250.18138.591.25780.25000.6006.02表5.負荷數據表節(jié)點號有功負荷（MW）無功負荷（MVA）51255069030810035思考題1、牛頓拉夫遜法與PQ分解法有哪些聯系？有哪些區(qū)別？二者的計算性能如何？2、為什么在用計算機對某網絡初次進行潮流計算時往往是要調潮流，而并非任何情 況下只一次送入初始值算出結果就行呢？要考慮什么條件？各變量是如何劃分 的？哪些可調？哪些不可調？3、發(fā)電機節(jié)點的注入無功為負

9、值說明了什么？4、線路有功潮流最有可能的流向是？線路無功潮流最有可能的流向是？5、負荷功率因數對系統(tǒng)潮流有什么影響？6、潮流控制的主要手段有哪些？7、如何降低系統(tǒng)網損？8、如何校驗主設備的短路容量？9、利用節(jié)點阻抗矩陣計算短路電流和利用導納矩陣計算短路電流的算法有何異同？10、減小短路電流的措施有哪些?設計主要內容：1、完成調試潮流計算與短路計算程序調試;計算程序以電力系統(tǒng)計算機輔助分析所用教材中的例程為基礎，要求掌握算 法原理，建議進行適當改寫以利于本次設計應用。2、根據設計資料準備計算數據；整理原始數據，使其符合計算程序的輸入條件。要求比較采用變壓器口型模型前 后兩種情況的差異。3、為指定

10、電網安排合理的穩(wěn)態(tài)運行方式；合理設置目標電網中的控制變量， 以1號機為平衡機，將潮流分布調整為可行的 運行方式（要求所有節(jié)點的電壓不得低于 1.0 p.u.,也不能高于1.05p.u.）繪制潮 流分布圖。分析該運行方式的特點，提出進一步降低網損的改進的措施， 并利用計算 結果簡單證明之。具體內容參考設計資料所附思考題。4、列表分析節(jié)點7發(fā)生各種短路故障時的短路電流，以及運行方式對故障電流的影 響；5、完成設計說明書。要求結構清晰，即從電網模型建立，到初態(tài)潮流分布計算，到 問題分析，再到調整后的方案計算。故障電流計算也類似。不必逐一回答所有思考題。參數計算、潮流調節(jié)、與潮流計算matlab源程序

11、%參數計算的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/220A2*0.4*L; % 線路各段電抗(均為標幺值，取 Ub=220KV,SB=100MVA , 下同)disp(X1);for i=1:6R(i)=100/220A2*L(i)*r(i);%加分號避免出現過程項，下同end%線路各段電阻disp('線路電阻

12、 R='); disp(R); %輸出 Rfor i=1:6Z(i尸R(i)+j*X1(i);Z0(i)=R(i)+3*j*X1(i);enddisp(線路正序、負序阻抗Z1=Z2='); disp(Z);for i=1:6Bl(i)=220A2/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;fo門=1:3;XT(i尸j*U0

13、(i)*242A2/SN(i)/220A2;XT0(i)=XT(i);end;disp('變壓器正、負序電抗值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:3K(i)=242/UGN(i);endfor i=1:3Xdg(i)=j*100/220A2*x(i)*K(i)*K(i);end%發(fā)電機暫態(tài)電抗disp('發(fā)電機暫態(tài)電抗Xdg='); disp(Xdg);%輸出XdgPL

14、=125,90,100;QL=50,30,35;for i=1:3ZL(i)=100/(PL(i)A2+QL(i)A2)*(PL(i)+j*QL(i);end%負荷節(jié)點阻抗disp('負荷節(jié)點阻抗Zl='); disp(ZL);%輸出ZL%支數參數形成的矩陣B1、支數參數形成的矩陣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;57R(3)+X1(3)*j,B(3),1,0;

15、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(6)+X1(6)*j,B(6),1,0;disp('支數參數形成的矩陣：B1='); disp(B1);%輸出支數參數形成的矩陣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,

16、0,0,2;0,1+0.35i,1,0,0,2;0,0,1,0,0,2;disp('支數參數形成的矩陣：B2='); disp(B2);%輸出支數參數形成的矩陣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('只考慮各點電壓水平，得到的滿足要求的數據如下：)；for aa=0:5B2(2,4)=1+0.01*aa;for ab=0:5B2(3,4)=1+0.01*ab;for ba=0:5B1(1,5)=0.95+

17、0.025*ba;for bb=0:5B1(2,5)=0.95+0.025*bb;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:nif X(i,2) =0;p=X(i,1);Y(p,p)=1./X(i,2);endendfor i=1:nlif B1(i,6)=0p=B1(i,1); q=B1(i,2);else p=B1(i,2); q=B1(i,1);endY(p,q尸Y(p,q)-1./(B1(i,3)

18、*B1(i,5);Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)A2)+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:ne(i尸real(B2(i,3);f(i尸imag(B2(i,3);V(i尸B2(i,4);endfor i=1:nS(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

19、 =0IT2=0;a=a+1;for i=1:nif i=isbC(i)=0；D(i)=0；for j1=1:nC(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);endP1=C(i)*e(i)+f(i)*D(i);Q1=f(i)*C(i)-D(i)*e(i);V2=e(i)A2+f(i)A2;if B2(i,6)=3DP=P(i)-P1;DQ=Q(i)-Q1;for j1=1:nif j1=isb & j1=iX1=-G(i,j1)*e(i)-B(i,j1)*f(i);X2=B(i,j1)

20、*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=isbX1=-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*j1-1;J(p,q)

21、=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)-P1;DV=V(i)A2-V2;for j1=1:nif j1=isb&j1=iX1=-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;J(m,q)=X2; elseif j1=i&

22、amp;j1=isbX1=-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;endendendendendfor k=3:N0k1=k+1;N1=N;for k2=k1:N1J(k,k2)=J(k,k2)./J(k,k);endJ(k,k)=1;if k=3k4=k-1;for k3=3:k4for

23、k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endif k=N0,break;endfor k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endelsefor k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendendfor k=3:2:N0-1L=(k+1).e(L尸e(L)-J(k,N);k1=k+1;f(L尸f(L)-J(k1,N);endfo

24、r k=3:N0DET=abs(J(k,N);if DET>=prIT2=IT2+1;endendICT2(a)=IT2;ICT1=ICT1+1;for k=1:ndy(k尸sqrt(e(k)A2+f(k)A2);endfor i=1:nDy(ICT1,i尸dy(i);endfor vb=1:nif Dy(1,vb)<=1.05&Dy(1,vb)>=1.00;zy(vb)=1;else zy(vb)=0;endendend %用高斯消去法解"w=-J*V"for k=1:nV(k尸sqrt(e(k)A2+f(k)A2);O(k尸atan(f(k).

25、/e(k)*180./pi;endE=e+f*j;for va=1:nif 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*yj(8)*yj(9)=1if zy(1)*zy(2)*zy(3)*zy(4)*zy(5)*zy(6)*zy(7)*zy(8)*zy(9)=1;dispC發(fā)電機2出口電壓：');disp(B2(2,4);dispC發(fā)電機3出口電壓：');disp(B2(3,4);disp('變壓器1接頭

26、電壓：');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:nC(p)=0；for q=1:nC(p)=C(p)+conj(Y(p,q)*conj(E(q);endS(p尸E(p)*C(p);endfor i=1:nlif B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endSi(p,q)=E(

27、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);endfor i=1:nlif B1(i,6)=0p=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);endZDS=0;for i=1:nlif B1(i,6)=0p=B1(i,1);q=B

28、1(i,2);else p=B1(i,2);q=B1(i,1); endDS(i尸Si(p,q)+Sj(q,p);ZDS=DS(i)+ZDS;endmm=mm+1;ZDSS(mm尸sqrt(real(ZDS)A2+imag(ZDS)A2);disp('*)；for i=1:nZB224(mm尸B2(2,4);ZB234(mm尸B2(3,4);ZB115(mm尸B1(1,5);ZB125(mm尸B1(2,5);ZB135(mm尸B1(3,5);endendendend end endendenddisp('*');disp('滿足各點電壓要求的數據組共有：);di

29、sp(mm);LOL=min(ZDSS);for dota=1:mmif ZDSS(dota)=LOL mm=dota;end end)；disp('當加入線路損耗時，得到其中線路損耗最小的一組組號及關鍵數據為: disp(mm);disp。發(fā)電機2出口電壓：');disp(ZB224(mm);disp。發(fā)電機3出口電壓：');disp(ZB234(mm);disp('變壓器1接頭電壓：');disp(ZB115(mm);disp('變壓器2接頭電壓：');disp(ZB125(mm);disp('變壓器3接頭電壓：');

30、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('用此組數據算的的潮流結果為：');disp('節(jié)點數：');disp(n);disp('節(jié)點支路數：');disp(nl);disp('平衡節(jié)點：');disp(isb);disp(精度：)；disp(pr);disp('B1=');disp(B1);disp('B2=');dis

31、p(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é)點數：');disp(n);disp('節(jié)點支路數：');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:nif X(i,2)=0;p=X(i,1);Y(p

32、,p)=1./X(i,2);endendfor i=1:nlif B1(i,6)=0p=B1(i,1); q=B1(i,2);else p=B1(i,2); q=B1(i,1); endY(p,q尸Y(p,q)-1./(B1(i,3)*B1(i,5);Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)A2)+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:ne(i尸real(B2(i,3);f(i尸imag(B2(i,3);V(i尸B2(i,4)

33、;endfor i=1:nS(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 -=0IT2=0;a=a+1;for i=1:n if i=isbC(i)=0;D(i)=0;for j1=1:nC(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);endP1=C(i)*e(i)+f(i)*D(i);Q1=f(i)*C(i)-D(i)*e(i);

34、V2=e(i)A2+f(i)A2;if B2(i,6)=3DP=P(i)-P1;DQ=Q(i)-Q1;for j1=1:nif j1=isb & j1=iX1=-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=isbX1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);

35、 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*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)A2-V2;for j1=1:nif j1=isb&j1=iX1=-G(i,j1)*e(i)-B(i,j1)*f(i);X2=B(i,j1)*e(i)

36、-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;J(m,q)=X2; elseif j1=i&j1=isbX1=-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)=

37、X6;J(m,q)=X2; end endendendendfor k=3:N0k1=k+1;N1=N;for k2=k1:N1J(k,k2)=J(k,k2)./J(k,k);endJ(k,k)=1;if k=3k4=k-1;for k3=3:k4for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endif k=N0,break;endfor k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endelsefor k3=k1:N0for k2=

38、k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendendfor k=3:2:N0-1L=(k+1).e(L尸e(L)-J(k,N);k1=k+1;f(L尸f(L)-J(k1,N);endfor k=3:N0DET=abs(J(k,N);if DET>=prIT2=IT2+1;endendICT2(a)=IT2;ICT1=ICT1+1;for k=1:ndy(k尸sqrt(e(k)A2+f(k)A2);endfor i=1:nDy(ICT1,i尸dy(i);endend %用高斯消去法解"w=-J*V"d

39、isp(迭代次數);disp(ICTI);dispC沒有達到精度要求的個數');disp(ICT2);for k=1:nV(k尸sqrt(e(k)A2+f(k)A2);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é)點的電壓角。為(節(jié)點號從小到大排列)：)；disp(O);for p=1:nC(p)=0；for q=1:nC(p)=C(p)+conj

40、(Y(p,q)*conj(E(q);endS(p尸E(p)*C(p);enddisp('各節(jié)點白功率S為(節(jié)點號從小到大排列)：');disp(S);disp('各條支路的首段功率 Si為(順序同您輸入B1時一樣)：)；for i=1:nlif B1(i,6)=0p=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);disp(Si(p,q)

41、;enddisp('各條支路的末段功率 Sj為(順序同您輸入B1時一樣)：)；for i=1:nlif B1(i,6)=0p=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);enddisp('各條支路的功率損耗DS為(順序同，輸入B1時一樣)：)；for i=1:nlif B1(i,6)=0p=B1(i,1);q=B1

42、(i,2);else p=B1(i,2);q=B1(i,1);endDS(i尸Si(p,q)+Sj(q,p);disp(DS(i);endfor i=1:ICT1Cs(i)=i;enddisp('以下是各節(jié)點每次迭代后的電壓值(如圖所示);plot(Cs,Dy),xlabel('迭代次數'),ylabel('電壓'),title('電壓迭代次數曲線');最后是課程設計報告反刃理）大老課程設計課程名稱：電力系統(tǒng)分析設計題目：基于Matlab計算程序的電力系統(tǒng)運行分析學院：電力工程學院專業(yè)：電氣工程自動化年 級：2009學生姓名易軍（你可以稱

43、我易大神，可惜 dota、lol都是菜鳥頂級水平）指導教師：宋琪、謝靜、單節(jié)杉日 期： 2013-1-13教務處制電力系統(tǒng)潮流計算是電力系統(tǒng)分析中的一種最基本的計算，是對復雜電力系統(tǒng)正常和故障 條件下穩(wěn)態(tài)運行狀態(tài)的計算。潮流計算的目標是求取電力系統(tǒng)在給定運行狀態(tài)的計算。即節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷.各點電壓是否滿足要求， 功率的分布和分配是否 合理以及功率損耗等。對現有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設計以及對電力 系統(tǒng)進行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計算為基礎。潮流計算結果可用如電力系統(tǒng)穩(wěn)態(tài)研究, 安全估計或最優(yōu)潮流等對潮流計算的模型和方法有直接影響。在電力

44、系統(tǒng)中可能發(fā)生的各種故障中，危害最大且發(fā)生概率較高的首推短路故障。產生短 路故障的主要原因是電力設備絕緣損壞。短路故障分為三相短路、兩相短路、單相接地短路及兩 相接地短路。其中三相短路時三相電流仍然對稱，其余三類短路統(tǒng)成為不對稱短路。短路故障大 多數發(fā)生在架空輸電線路。電力系統(tǒng)設計與運行時，要采取適當的措施降低短路故障的發(fā)生概率。短路計算可以為設備的選擇提供原始數據。目標電網接線圖、電網模型的建立設計中，選取 SB=100MVA, U b=220KV ,將所有支路的參數都折算到220KV電壓等級側,本文利用matlab計算參數其源程序為：clc;clear;L=104,110,98,163,8

45、8,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/220a2*0.4*L;%線路各段電抗(均為標幺值，取 disp(X1);for i=1:6R(i)=100/220A2*L(i)*r(i);enddisp('線路電阻R=');disp(R)；for i=1:6Z(i尸R(i)+j*X1(i);Z0(i)=R(i)+3*j*X1(i);enddisp('

46、;線路正序、負序阻抗 Z1=Z2=');disp(Z)；for i=1:6Bl(i)=220A2/100*b(i)*L(i);B(i)=0+j*Bl(i);enddisp('線路容納之半 B/2='); disp(B);disp('線路零序阻抗Z0:');disp(Z0);U0=11,10,6;SN=240,220,120;UN=242,242,242;fo門=1:3;XT(i)=j*U0(i)*242A2/SN(i)/220A2;XT0(i尸XT(i);end;disp('變壓器正、負序電抗值 Xt： ');disp(XT);disp(

47、'變壓器零序電抗值Xt0： ');disp(XT0);Ub=220KV,SB=100MV A)%加分號避免出現過程項，下同%線路各段電阻%輸出R%線路各容納之半%輸出B%發(fā)電機暫態(tài)電抗有名值x=0.0608,0.1198,0.1813;UGN=16.5,18,13.8;for i=1:3K(i)=242/UGN(i);endfor i=1:3Xdg(i)=j*100/220A2*x(i)*K(i)*K(i);end%發(fā)電機暫態(tài)電抗disp('發(fā)電機暫態(tài)電抗Xdg='); disp(Xdg);PL=125,90,100;QL=50,30,35;for i=1:3Z

48、L(i)=100/(PL(i)A2+QL(iF2)*(PL(i)+j*QL(i);end%負荷節(jié)點阻抗disp('負荷節(jié)點阻抗 Zl='); disp(ZL);0.0101 + 0.0860i0.0168 + 0.0909i0.0160 + 0.0810i0.0310 + 0.1347i0.0085 + 0.0727i計算結果為:線路電抗X=0.08600.09090.08100.13470.07270.1008線路電阻R=0.01010.01680.01600.03100.00850.0118線路正序、負序阻抗 Z1=Z2=Columns 1 through 5Column

49、60.0118 + 0.1008i線路容納之半B/2=Columns 1 through 50 + 0.0896i0 + 0.0783i0 + 0.0759i0 + 0.1420i0 + 0.0758iColumn 60 + 0.1051i線路零序阻抗Z0：Columns 1 through 50.0101 + 0.2579i0.0168 + 0.2727i0.0160 + 0.2430i0.0310 + 0.4041i0.0085 + 0.2182iColumn 60.0118 + 0.3025i變壓器正、負序電抗值Xt：0 + 0.0555i0.0550i0.0605i變壓器零序電抗值Xt

50、0 :0 + 0.0555i0.0550i0.0605i發(fā)電機暫態(tài)電抗 Xdg=0 + 0.0270i負荷節(jié)點阻抗Zl=0.0447i0.1152i0.8909 + 0.3118i0.6897 + 0.2759i1.0000 + 0.3333i支數參數形成的矩陣：B1 =Columns 1 through 51.00004.00000+0.0555i01.00002.00007.00000+0.0550i01.00003.00009.00000+0.0605i01.00004.00005.00000.0101+0.0860i0+0.0896i1.00004.00006.00000.0168+0

51、.0909i0+0.0783i1.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.0000Column 6000000000支數參數形成的矩陣：B2=Columns 1 through 52.0000 + 1.2395i01.05001.000001.8000+1.1155i01.00001.000001.0000+0

52、.6197i01.00001.00000001.00000001.2500 +0.5000i1.00000000.9000 +0.3000i1.000000001.00000001.0000 +0.3500i1.000000001.000000Column 61.00003.00003.00002.00002.00002.00002.00002.00002.0000二.系統(tǒng)等值電路圖的繪制根據以上計算結果，得到系統(tǒng)等值電路圖如下:J0XJ4470.0310+JDJ07E7OOI10+JD.1OQ8JD.U5gn-iH 一? JOLIET®C.116EHjaC®10r JD.

53、I051 加05】工Ja17&S F"加739HH"OU033104JIM34JD.14E0HH-0.6S97+J0.2759IIJ0L1I頜OAim-i-joins&oJOJ06%TIII J0,D719 001gg"Hong 嶺 旭通HIJQJ0555J0J0270E1二、潮流計算1 .系統(tǒng)參數的設置設計中要求所有結點電壓不得低于1.0p.u.,也不得高于1.05p.u.,若電壓不符合該條件，可采取下面的方法進行調壓：(1) 改變發(fā)電機的機端電壓(2) 改變變壓器的變比(即改變分接頭)(3) 改變發(fā)電機的出力(4) 在電壓不符合要求的結點處增加無功補