電力系統(tǒng)分析大作業(yè)-matlab三機九節(jié)點潮流計算報告-文庫吧資料
2025-01-24 23:31本頁面
【正文】 2(k,3)*(real(Y(i,k))*sin(B2(i,4)B2(k,4))imag(Y(i,k))*cos(B2(i,4)B2(k,4)))。 else Jacobian(n11+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)))。 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)))。 else Jacobian(i,n11+j)=0B2(i,3)*(real(Y(i,j))*cos(B2(i,4)B2(j,4))+imag(Y(i,j))*sin(B2(i,4)B2(j,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)))。 else Jacobian(i,j)=0B2(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)))。 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)))。 Mismatch_power Jacobian=zeros(l+m*2)?！?失配功率:】39。)。 end disp(39。 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(i,1)=real(B2(i,1))real(B2(i,2))Pj。while(max(Mismatch_power)Accuracy) for i=1:(n11) Pj=0。 end Mismatch_power(k,1)=imag(B2((kn1+1),1))imag(B2((kn1+1),2))Qj。 end for k=n1:(l+m*2) Qj=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))))。 endendl %PV節(jié)點個數(shù)Mismatch_power=zeros(l+m*2,1)。 endendm %PQ節(jié)點個數(shù)l=0。Y %顯示導(dǎo)納矩陣m=0?！?導(dǎo)納矩陣:】39。 Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3))。 Y(q,p)=Y(p,q)。 q=B1(i,2)。 Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)。 Y(q,p)=Y(p,q)。 q=B1(i,2)。B2 %顯示節(jié)點參數(shù)矩陣% 以上為從excel中導(dǎo)入初值的程序Y=zeros(n1)。【 節(jié)點參數(shù)矩陣:】39。)。enddisp(39。 B2(j,4)=Data_B2(j,6)。 B2(j,2)=T2(j,2)。 B1(j,6)=Data_B1(j,7)。 B1(j,4)=T1(j
研究報告相關(guān)推薦
鄂ICP備17016276號-1