電力系統(tǒng)分析大作業(yè)-matlab三機(jī)九節(jié)點(diǎn)潮流計(jì)算報(bào)告_第1頁
電力系統(tǒng)分析大作業(yè)-matlab三機(jī)九節(jié)點(diǎn)潮流計(jì)算報(bào)告_第2頁
電力系統(tǒng)分析大作業(yè)-matlab三機(jī)九節(jié)點(diǎn)潮流計(jì)算報(bào)告_第3頁
電力系統(tǒng)分析大作業(yè)-matlab三機(jī)九節(jié)點(diǎn)潮流計(jì)算報(bào)告_第4頁
電力系統(tǒng)分析大作業(yè)-matlab三機(jī)九節(jié)點(diǎn)潮流計(jì)算報(bào)告_第5頁
已閱讀5頁，還剩12頁未讀， 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、.電力系統(tǒng)分析大作業(yè) 一、設(shè)計(jì)題目本次設(shè)計(jì)題目選自課本第五章例5-8,美國(guó)西部聯(lián)合電網(wǎng)WSCC系統(tǒng)的簡(jiǎn)化三機(jī)九節(jié)點(diǎn)系統(tǒng)，例題中已經(jīng)給出了潮流結(jié)果，計(jì)算結(jié)果可以與之對(duì)照。取=0.00001 。二、計(jì)算步驟 第一步，為了方便編程，修改節(jié)點(diǎn)的序號(hào)，將平衡節(jié)點(diǎn)放在最后。如下圖： 92132745683 第二步，這樣得出的系統(tǒng)參數(shù)如下表所示： 第三步，形成節(jié)點(diǎn)導(dǎo)納矩陣。 第四步，設(shè)定初值： ； ，。 第五步，計(jì)算失配功率 =0，=-1.25，=-0.9，=0，=-1，=0，=1.63， =0.85； =0.8614，=-0.2590，=-0.0420，=0.6275，=-0.1710，=0.7101。

2、顯然，。 第六步，形成雅克比矩陣（階數(shù)為14×14）第七步，解修正方程，得到：-0.0371，-0.0668，-0.0628，0.0732，0.0191，0.0422，0.1726，0.0908；0.0334，0.0084，0.0223，0.0372，0.0266，0.0400。從而-0.0371，-0.0668，-0.0628，0.0732，0.0191，0.0422，0.1726，0.0908； 1.0334，1.0084，1.0223，1.0372，1.0266，1.0400。然后轉(zhuǎn)入下一次迭代。經(jīng)三次迭代后。迭代過程中失配功率的變化情況如下表：k0123=P10-0.0106

3、0.00010.0000000845=P2-1.250.03790.00050.0000000896=P3-0.90.04390.00050.0000000778=P40-0.0421-0.0012-0.0000003421=P5-10.0610.00090.0000001845=P60-0.0269-0.0007-0.0000001631=P71.63-0.0579-0.0004-0.0000000278=P80.85-0.0336-0.0002-0.0000000103=Q10.8614-0.0501-0.0004-0.0000000561=Q2-0.259-0.0714-0.0012-0

4、.0000002774=Q3-0.042-0.0424-0.0006-0.0000001236=Q40.6275-0.1875-0.0021-0.0000003279=Q5-0.171-0.0241-0.0004-0.0000000805=Q60.7101-0.0828-0.0007-0.0000000799max1.630.0610.00090.0000001845迭代過程中節(jié)點(diǎn)電壓變化情況如下表：kU1U2U3U4U5U6011111111.0334 1.0084 1.0223 1.0372 1.0266 1.0400 21.0259 0.9958 1.0128 1.0259 1.0160

5、 1.0324 31.0258 0.9956 1.0127 1.0258 1.0159 1.0324 迭代收斂后各節(jié)點(diǎn)的電壓和功率：kUPQ11.0258 -2.2168 0.0000 0.0000 20.9956 -3.9888 -1.2500 -0.5000 31.0127 -3.6874 -0.9000 -0.3000 41.0258 3.7197 0.0000 0.0000 51.0159 0.7275 -1.0000 -0.3500 61.0324 1.9667 0.0000 0.0000 71.0250 9.2800 1.6300 0.0665 81.0250 4.6648 0.8

6、500 -0.1086 91.0400 0.0000 0.7164 0.2705 同課本上給出的潮流相比較，結(jié)果完全一致，證明計(jì)算過程與程序編寫正確。 最后得出迭代收斂后各支路的功率和功率損耗：ijPijQijIijPjiQjiIjiPLQL120.4094 0.2289 0.4572 -0.4068 -0.3869 0.5639 0.0026 -0.1579 130.3070 0.0103 0.2995 -0.3054 -0.1654 0.3430 0.0017 -0.1551 24-0.8432 -0.1131 0.8545 0.8662 -0.0838 0.8484 0.0230 -0.

7、1969 36-0.5946 -0.1346 0.6020 0.6082 -0.1807 0.6146 0.0135 -0.3153 450.7638 -0.0080 0.7447 -0.7590 -0.1070 0.7546 0.0048 -0.1150 56-0.2410 -0.2430 0.3368 0.2418 0.0312 0.2362 0.0009 -0.2118 910.7164 0.2705 0.7363 -0.7164 -0.2392 0.7363 0.0000 0.0312 741.6300 0.0665 1.5916 -1.6300 0.0918 1.5916 0.000

8、0 0.1583 860.8500 -0.1086 0.8360 -0.8500 0.1496 0.8360 0.0000 0.0410 三、源程序及注釋由于計(jì)算流程比較簡(jiǎn)單，所以編寫程序過程中沒有采用模塊化的形式，直接按順序一步步進(jìn)行。disp('【 節(jié)點(diǎn)數(shù):】');n1=xlsread('input.xls','A3:A3')%節(jié)點(diǎn)數(shù)disp('【 支路數(shù):】');n=xlsread('input.xls','B3:B3')%支路數(shù)disp('【 精度:】');Accuracy=

9、xlsread('input.xls','B4:B4')%精度branch=xlsread('input.xls','E4:K12');node=xlsread('input.xls','M4:S12');Data_B1=branch;%支路參數(shù)Data_B2=node;%節(jié)點(diǎn)參數(shù)T1=zeros(n,2);T2=zeros(n1,3);i=sqrt(-1);format shortfor j=1:n T1(j,1)=Data_B1(j,3)+Data_B1(j,4)*1i; T1(j,2)=Dat

10、a_B1(j,5)*1i;endfor j=1:n1 T2(j,1)=Data_B2(j,1)+Data_B2(j,2)*1i; T2(j,2)=Data_B2(j,3)+Data_B2(j,4)*1i;endB1=zeros(n,6);B2=zeros(n1,5);for j=1:n B1(j,1)=Data_B1(j,1); B1(j,2)=Data_B1(j,2); B1(j,3)=T1(j,1); B1(j,4)=T1(j,2); B1(j,5)=Data_B1(j,6); B1(j,6)=Data_B1(j,7);endfor j=1:n1 B2(j,1)=T2(j,1); B2(j

11、,2)=T2(j,2); B2(j,3)=Data_B2(j,5); B2(j,4)=Data_B2(j,6); B2(j,5)=Data_B2(j,7);enddisp('【 支路參數(shù)矩陣:】');B1 %顯示支路參數(shù)矩陣disp('【 節(jié)點(diǎn)參數(shù)矩陣:】');B2 %顯示節(jié)點(diǎn)參數(shù)矩陣% 以上為從excel中導(dǎo)入初值的程序Y=zeros(n1);for i=1:n if B1(i,6)=0 %不含變壓器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/B1(i,3); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1

12、/B1(i,3)+0.5*B1(i,4); Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4); else %含有變壓器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3); Y(q,q)=Y(q,q)+1/(B1(i,5)2*B1(i,3); endenddisp('【 導(dǎo)納矩陣:】');Y %顯示導(dǎo)納矩陣m=0;for i=1:n1 if B2(i,5)=2 m=m+1; endendm %PQ節(jié)點(diǎn)個(gè)數(shù)l=0;f

13、or i=1:n1 if B2(i,5)=1 l=l+1; endendl %PV節(jié)點(diǎn)個(gè)數(shù)Mismatch_power=zeros(l+m*2,1);for i=1:n1-1 Pj=0; for j=1:n1 Pj=Pj+(B2(i,3)*B2(j,3)*(real(Y(i,j)*cos(B2(i,4)-B2(j,4)+imag(Y(i,j)*sin(B2(i,4)-B2(j,4); end Mismatch_power(i,1)=real(B2(i,1)-real(B2(i,2)-Pj; end for k=n1:(l+m*2) Qj=0; for j=1:n1 Qj=Qj+B2(k-n1+

14、1),3)*B2(j,3)*(real(Y(k-n1+1),j)*sin(B2(k-n1+1),4)-B2(j,4)-imag(Y(k-n1+1),j)*cos(B2(k-n1+1),4)-B2(j,4); end Mismatch_power(k,1)=imag(B2(k-n1+1),1)-imag(B2(k-n1+1),2)-Qj; end% Mismatch_power %計(jì)算失配功率times=0;while(max(Mismatch_power)>Accuracy) for i=1:(n1-1) Pj=0; for j=1:n1 Pj=Pj+B2(i,3)*B2(j,3)*(r

15、eal(Y(i,j)*cos(B2(i,4)-B2(j,4)+imag(Y(i,j)*sin(B2(i,4)-B2(j,4); end Mismatch_power(i,1)=real(B2(i,1)-real(B2(i,2)-Pj; end for k=n1:(l+m*2) Qj=0; for j=1:n1 Qj=Qj+B2(k-n1+1),3)*B2(j,3)*(real(Y(k-n1+1),j)*sin(B2(k-n1+1),4)-B2(j,4)-imag(Y(k-n1+1),j)*cos(B2(k-n1+1),4)-B2(j,4); end Mismatch_power(k,1)=im

16、ag(B2(k-n1+1),1)-imag(B2(k-n1+1),2)-Qj; end disp('【 當(dāng)前迭代次數(shù):】'); times disp('【 失配功率:】'); Mismatch_power Jacobian=zeros(l+m*2);%雅克比矩陣7*7 %H for i=1:(n1-1) for j=1:(n1-1) if i=j P_H=0; for k=1:n1 P_H=P_H+B2(i,3)*B2(k,3)*(real(Y(i,k)*sin(B2(i,4)-B2(k,4)-imag(Y(i,k)*cos(B2(i,4)-B2(k,4); e

17、nd Jacobian(i,i)=P_H-B2(i,3)*B2(i,3)*(0-imag(Y(i,i); else Jacobian(i,j)=0-B2(i,3)*B2(j,3)*(real(Y(i,j)*sin(B2(i,4)-B2(j,4)-imag(Y(i,j)*cos(B2(i,4)-B2(j,4); end end end %N for i=1:(n1-1) for j=1:m if i=j P_N=0; for k=1:n1 P_N=P_N+B2(k,3)*(real(Y(i,k)*cos(B2(i,4)-B2(k,4)+imag(Y(i,k)*sin(B2(i,4)-B2(k,4

18、); end Jacobian(i,n1-1+i)=0-B2(i,3)*real(Y(i,i)-P_N; else Jacobian(i,n1-1+j)=0-B2(i,3)*(real(Y(i,j)*cos(B2(i,4)-B2(j,4)+imag(Y(i,j)*sin(B2(i,4)-B2(j,4); end end end %K for i=1:m for j=1:(n1-1) if i=j P_K=0; for k=1:n1 P_K=P_K+B2(i,3)*B2(k,3)*(real(Y(i,k)*cos(B2(i,4)-B2(k,4)+imag(Y(i,k)*sin(B2(i,4)-B

19、2(k,4); end Jacobian(n1-1+i,i)=0+B2(i,3)*B2(i,3)*real(Y(i,i)-P_K; else Jacobian(n1-1+i,j)=B2(i,3)*B2(j,3)*(real(Y(i,j)*cos(B2(i,4)-B2(j,4)+imag(Y(i,j)*sin(B2(i,4)-B2(j,4); end end end %L for i=1:m for j=1:m if i=j P_L=0; for k=1:n1 P_L=P_L+B2(k,3)*(real(Y(i,k)*sin(B2(i,4)-B2(k,4)-imag(Y(i,k)*cos(B2(

20、i,4)-B2(k,4); end Jacobian(n1-1+i,n1-1+i)=0-P_L+B2(i,3)*imag(Y(i,i); else Jacobian(n1-1+i,n1-1+j)=0-B2(i,3)*(real(Y(i,j)*sin(B2(i,4)-B2(j,4)-imag(Y(i,j)*cos(B2(i,4)-B2(j,4); end end end S=zeros(l+m*2,1); %初始化電壓角度變化量 S=inv(Jacobian)*(0-Mismatch_power); %求解修正方程 S=(Jacobian)(0-Mismatch_power); %求解修正方程

21、for i=1:(n1-1) %角度初值加變化量 B2(i,4)=B2(i,4)+S(i,1); end for i=1:m %電壓初值加變化量 B2(i,3)=B2(i,3)+S(n1-1+i,1); end disp('【 雅克比矩陣:】'); Jacobian %顯示雅克比矩陣% S=inv(Jacobian) times=times+1;endtimes=times-1;disp('【 共計(jì)迭代次數(shù):】');times %顯示迭代次數(shù)U_It=zeros(n1,1); %初始化電壓向量for i=1:n1 U_It(i,1)=B2(i,3)*cos(B2

22、(i,4)+B2(i,3)*sin(B2(i,4)*1j;endangle_It=zeros(n1,1); %將電壓角度的弧度值轉(zhuǎn)為角度值for i=1:n1 angle_It(i,1)=B2(i,4)*180/pi;endNode_S_It=U_It.*(conj(Y)*conj(U_It); %求解節(jié)點(diǎn)功率disp('【 迭代收斂后各節(jié)點(diǎn)的電壓幅值:】');Node_U_It=abs(U_It) %顯示迭代收斂后各節(jié)點(diǎn)的電壓幅值disp('【 迭代收斂后各節(jié)點(diǎn)的電壓角度:】');angle_It %顯示迭代收斂后各節(jié)點(diǎn)的電壓角度disp('【 迭代收

23、斂后各節(jié)點(diǎn)的功率:】');Node_S_It %顯示迭代收斂后各節(jié)點(diǎn)的功率Branch_It=zeros(n,10);for i=1:n; if B1(i,6)=0; %不帶變壓器支路 m=B1(i,1); %得到支路號(hào) n=B1(i,2); Branch_It(i,1)=m; %顯示支路號(hào) Branch_It(i,2)=n; a=U_It(m,1)*(conj(U_It(m,1)*conj(B1(i,4)*0.5+(conj(U_It(m,1)-conj(U_It(n,1)/conj(B1(i,3); Branch_It(i,3)=real(a); %顯示Pij Branch_It(

24、i,4)=imag(a); %顯示Qij b=U_It(m,1)*B1(i,4)*0.5+(U_It(m,1)-U_It(n,1)/B1(i,3); Branch_It(i,5)=sqrt(real(b)2+imag(b)2); %顯示Iij c=U_It(n,1)*(conj(U_It(n,1)*conj(B1(i,4)*0.5+(conj(U_It(n,1)-conj(U_It(m,1)/conj(B1(i,3); Branch_It(i,6)=real(c); %顯示Pji Branch_It(i,7)=imag(c); %顯示Qji d=U_It(n,1)*B1(i,4)*0.5+(U

25、_It(n,1)-U_It(m,1)/B1(i,3); Branch_It(i,8)=sqrt(real(d)2+imag(d)2); %顯示Iji e=a+c; Branch_It(i,9)=real(e); %顯示線路損耗有功分量 Branch_It(i,10)=imag(e); %顯示線路損耗無功分量 else %帶變壓器支路(同以上內(nèi)容) m=B1(i,1); n=B1(i,2); Branch_It(i,1)=m; Branch_It(i,2)=n; a=U_It(m,1)*(conj(U_It(m,1)/conj(B1(i,3)-conj(U_It(n,1)*conj(1/(B1(

26、i,5)*B1(i,3); Branch_It(i,3)=real(a); Branch_It(i,4)=imag(a); b=U_It(m,1)*(B1(i,5)-1)/B1(i,3)/B1(i,5)+(U_It(m,1)-U_It(n,1)/(B1(i,5)*B1(i,3); Branch_It(i,5)=sqrt(real(b)2+imag(b)2); c=U_It(n,1)*(conj(U_It(n,1)/(conj(B1(i,5)*B1(i,5)*B1(i,3)-conj(U_It(m,1)*conj(1/(B1(i,5)*B1(i,3); Branch_It(i,6)=real(c); Branch_It(i,7)=imag(c); d=U_It(n,1)*(1-B1(i,5)/B1(i

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間，僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請(qǐng)與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

最新文檔

評(píng)論

0/150

提交評(píng)論