電力系統(tǒng)分析大作業(yè)-matlab三機九節(jié)點潮流計算報告-wenkub
2023-02-02 23:31:05 本頁面
【正文】 39。 Y(p,q)=Y(p,q)1/(B1(i,3)*B1(i,5))。 Y(p,p)=Y(p,p)+1/B1(i,3)+*B1(i,4)。for i=1:n if B1(i,6)==0 %不含變壓器的支路 p=B1(i,1)。B1 %顯示支路參數(shù)矩陣disp(39。 B2(j,5)=Data_B2(j,7)。endfor j=1:n1 B2(j,1)=T2(j,1)。 B1(j,3)=T1(j,1)。endB1=zeros(n,6)。format shortfor j=1:n T1(j,1)=Data_B1(j,3)+Data_B1(j,4)*1i。%支路參數(shù)Data_B2=node。,39。E4:K1239。B4:B439。)。,39。【 支路數(shù):】39。39。disp(39。然后轉(zhuǎn)入下一次迭代。 第五步，計算失配功率 =0，=，=，=0，=1，=0，=， =； =，=，=，=，=，=。取ε= 。二、計算步驟 第一步，為了方便編程，修改節(jié)點的序號，將平衡節(jié)點放在最后。 顯然。經(jīng)三次迭代后?！?節(jié)點數(shù):】39。,39。)。B3:B339。Accuracy=xlsread(39。)%精度[branch]=xlsread(39。)。M4:S1239。%節(jié)點參數(shù)T1=zeros(n,2)。 T1(j,2)=Data_B1(j,5)*1i。B2=zeros(n1,5)。 B1(j,4)=T1(j,2)。 B2(j,2)=T2(j,2)。enddisp(39?！?節(jié)點參數(shù)矩陣:】39。 q=B1(i,2)。 Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)。 Y(q,p)=Y(p,q)?！?導(dǎo)納矩陣:】39。 endendm %PQ節(jié)點個數(shù)l=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(k,1)=imag(B2((kn1+1),1))imag(B2((kn1+1),2))Qj。 end Mismatch_power(i,1)=real(B2(i,1))real(B2(i,2))Pj。 end disp(39?！?失配功率:】39。 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)))。 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)))。 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)B2(k,4)))。 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(i,4)B2(k,4)))。 %初始化電壓角度變化量 S=inv(Jacobian)*(0Mismatch_power)。 end disp(39。endtimes=times1。times %顯示迭代次數(shù)U_It=zeros(n1,1)。endNode_S_It=U_It.*(conj(Y)*conj(U_It))。Node_U_It=abs(
研究報告相關(guān)推薦
鄂ICP備17016276號-1