昆明理工大學(xué)計算機(jī)輔助分析課程設(shè)計及matlab自動實現(xiàn)潮流調(diào)節(jié)與計算講解

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

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

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

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

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

6、進(jìn)行綜合評定。原始數(shù)據(jù)與思考題目標(biāo)電網(wǎng)單線圖系統(tǒng)參數(shù)表1.線路參數(shù)表線路編號線路型號線路長度（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.變壓器參數(shù)表線路編號變壓器型號變壓器變比（kV）短路電壓白分?jǐn)?shù)（）2-7SSPL-220000242±3X 2.5%/20103-9SSPL-120000242±3X 2.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ī)母線名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.負(fù)荷數(shù)據(jù)表節(jié)點號有功負(fù)荷（MW）無功負(fù)荷（MVA）51255069030810035思考題1、牛頓拉夫遜法與PQ分解法有哪些聯(lián)系？有哪些區(qū)別？二者的計算性能如何？2、為什么在用計算機(jī)對某網(wǎng)絡(luò)初次進(jìn)行潮流計算時往往是要調(diào)潮流，而并非任何情 況下只一次送入初始值算出結(jié)果就行呢？要考慮什么條件？各變量是如何劃分 的？哪些可調(diào)？哪些不可調(diào)？3、發(fā)電機(jī)節(jié)點的注入無功為負(fù)

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

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

11、%參數(shù)計算的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; % 線路各段電抗(均為標(biāo)幺值，取 Ub=220KV,SB=100MVA , 下同)disp(X1);for i=1:6R(i)=100/220A2*L(i)*r(i);%加分號避免出現(xiàn)過程項，下同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(線路正序、負(fù)序阻抗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('變壓器正、負(fù)序電抗值Xt： ');disp(XT);disp('變壓器零序電抗值 Xt0： ');disp(XT0);x=0.0608,0.1198,0.1813;%發(fā)電機(jī)暫態(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ā)電機(jī)暫態(tài)電抗disp('發(fā)電機(jī)暫態(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%負(fù)荷節(jié)點阻抗disp('負(fù)荷節(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;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('支數(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,

16、0,0,2;0,1+0.35i,1,0,0,2;0,0,1,0,0,2;disp('支數(shù)參數(shù)形成的矩陣：B2='); disp(B2);%輸出支數(shù)參數(shù)形成的矩陣B2%潮流自動調(diào)節(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: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ā)電機(jī)2出口電壓：');disp(B2(2,4);dispC發(fā)電機(jī)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('滿足各點電壓要求的數(shù)據(jù)組共有：);di

29、sp(mm);LOL=min(ZDSS);for dota=1:mmif ZDSS(dota)=LOL mm=dota;end end)；disp('當(dāng)加入線路損耗時，得到其中線路損耗最小的一組組號及關(guān)鍵數(shù)據(jù)為: disp(mm);disp。發(fā)電機(jī)2出口電壓：');disp(ZB224(mm);disp。發(fā)電機(jī)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('用此組數(shù)據(jù)算的的潮流結(jié)果為：');disp('節(jié)點數(shù)：');disp(n);disp('節(jié)點支路數(shù)：');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é)點數(shù)：');disp(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: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(迭代次數(shù));disp(ICTI);dispC沒有達(dá)到精度要求的個數(shù)');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é)點的實際電壓標(biāo)幺值 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('迭代次數(shù)'),ylabel('電壓'),title('電壓迭代次數(shù)曲線');最后是課程設(shè)計報告反刃理）大老課程設(shè)計課程名稱：電力系統(tǒng)分析設(shè)計題目：基于Matlab計算程序的電力系統(tǒng)運(yùn)行分析學(xué)院：電力工程學(xué)院專業(yè)：電氣工程自動化年 級：2009學(xué)生姓名易軍（你可以稱

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

44、系統(tǒng)中可能發(fā)生的各種故障中，危害最大且發(fā)生概率較高的首推短路故障。產(chǎn)生短 路故障的主要原因是電力設(shè)備絕緣損壞。短路故障分為三相短路、兩相短路、單相接地短路及兩 相接地短路。其中三相短路時三相電流仍然對稱，其余三類短路統(tǒng)成為不對稱短路。短路故障大 多數(shù)發(fā)生在架空輸電線路。電力系統(tǒng)設(shè)計與運(yùn)行時，要采取適當(dāng)?shù)拇胧┙档投搪饭收系陌l(fā)生概率。短路計算可以為設(shè)備的選擇提供原始數(shù)據(jù)。目標(biāo)電網(wǎng)接線圖、電網(wǎng)模型的建立設(shè)計中，選取 SB=100MVA, U b=220KV ,將所有支路的參數(shù)都折算到220KV電壓等級側(cè),本文利用matlab計算參數(shù)其源程序為：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;%線路各段電抗(均為標(biāo)幺值，取 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、;線路正序、負(fù)序阻抗 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('變壓器正、負(fù)序電抗值 Xt： ');disp(XT);disp(

47、'變壓器零序電抗值Xt0： ');disp(XT0);Ub=220KV,SB=100MV A)%加分號避免出現(xiàn)過程項，下同%線路各段電阻%輸出R%線路各容納之半%輸出B%發(fā)電機(jī)暫態(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ā)電機(jī)暫態(tài)電抗disp('發(fā)電機(jī)暫態(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%負(fù)荷節(jié)點阻抗disp('負(fù)荷節(jié)點阻抗 Zl='); disp(ZL);0.0101 + 0.0860i0.0168 + 0.0909i0.0160 + 0.0810i0.0310 + 0.1347i0.0085 + 0.0727i計算結(jié)果為:線路電抗X=0.08600.09090.08100.13470.07270.1008線路電阻R=0.01010.01680.01600.03100.00850.0118線路正序、負(fù)序阻抗 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變壓器正、負(fù)序電抗值Xt：0 + 0.0555i0.0550i0.0605i變壓器零序電抗值Xt

50、0 :0 + 0.0555i0.0550i0.0605i發(fā)電機(jī)暫態(tài)電抗 Xdg=0 + 0.0270i負(fù)荷節(jié)點阻抗Zl=0.0447i0.1152i0.8909 + 0.3118i0.6897 + 0.2759i1.0000 + 0.3333i支數(shù)參數(shù)形成的矩陣：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支數(shù)參數(shù)形成的矩陣：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)等值電路圖的繪制根據(jù)以上計算結(jié)果，得到系統(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)參數(shù)的設(shè)置設(shè)計中要求所有結(jié)點電壓不得低于1.0p.u.,也不得高于1.05p.u.,若電壓不符合該條件，可采取下面的方法進(jìn)行調(diào)壓：(1) 改變發(fā)電機(jī)的機(jī)端電壓(2) 改變變壓器的變比(即改變分接頭)(3) 改變發(fā)電機(jī)的出力(4) 在電壓不符合要求的結(jié)點處增加無功補(bǔ)