電力系統(tǒng)分析大作業(yè)-matlab三機(jī)九節(jié)點(diǎn)潮流計(jì)算報(bào)告-文庫(kù)吧
2025-01-03 23:31 本頁(yè)面
【正文】 a_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,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(39?！?支路參數(shù)矩陣:】39。)。B1 %顯示支路參數(shù)矩陣disp(39?！?節(jié)點(diǎn)參數(shù)矩陣:】39。)。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/B1(i,3)+*B1(i,4)。 Y(q,q)=Y(q,q)+1/B1(i,3)+*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(39?！?導(dǎo)納矩陣:】39。)。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。for 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:n11 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((kn1+1),3)*B2(j,3)*(real(Y((kn1+1),j))*sin(B2((kn1+1),4)B2(j,4))imag(Y((kn1+1),j))*cos(B2((kn1+1),4)B2(j,4)))。 end Mismatch_power(k,1)=imag(B2((kn1+1),1))imag(B2((kn1+1),2))Qj。 end% Mismatch_power %計(jì)算失配功率times=0。while(max(Mismatch_power)Accuracy) for i=1:(n11) 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)
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1