電力系統(tǒng)分析大作業(yè)-matlab三機(jī)九節(jié)點(diǎn)潮流計(jì)算報(bào)告(存儲(chǔ)版)
【正文】 c=U_It(n,1)*(conj(U_It(n,1))*conj(B1(i,4))*+(conj(U_It(n,1))conj(U_It(m,1)))/conj(B1(i,3)))。 %顯示線路損耗無功分量 else %帶變壓器支路(同以上內(nèi)容) m=B1(i,1)。 Branch_It(i,5)=sqrt(real(b)^2+imag(b)^2)。 Branch_It(i,10)=imag(e)。,Node_U_It,1,39。 % xlswrite(39。E339。四、程序運(yùn)行結(jié)果五、手算結(jié)果（第一次迭代）六、總結(jié) 通過采用計(jì)算機(jī)和手算進(jìn)行潮流計(jì)算，我對(duì)潮流計(jì)算的計(jì)算過程和MATLAB軟件的使用有了更深層次的了解。)。39。C339。 % xlswrite(39。 e=a+c。 Branch_It(i,4)=imag(a)。 Branch_It(i,9)=real(e)。 %顯示Qij b=U_It(m,1)*B1(i,4)*+(U_It(m,1)U_It(n,1))/B1(i,3)。 if B1(i,6)==0。Node_U_It=abs(U_It) %顯示迭代收斂后各節(jié)點(diǎn)的電壓幅值disp(39。times %顯示迭代次數(shù)U_It=zeros(n1,1)。 end disp(39。 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)))。 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)))?！?失配功率:】39。 end Mismatch_power(i,1)=real(B2(i,1))real(B2(i,2))Pj。 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))))?！?導(dǎo)納矩陣:】39。 Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)?！?節(jié)點(diǎn)參數(shù)矩陣:】39。 B2(j,2)=T2(j,2)。B2=zeros(n1,5)。%節(jié)點(diǎn)參數(shù)T1=zeros(n,2)。)。Accuracy=xlsread(39。)?！?節(jié)點(diǎn)數(shù):】39。 顯然。取ε= 。然后轉(zhuǎn)入下一次迭代。39。,39。B4:B439。,39。format shortfor j=1:n T1(j,1)=Data_B1(j,3)+Data_B1(j,4)*1i。 B1(j,3)=T1(j,1)。 B2(j,5)=Data_B2(j,7)。for i=1:n if B1(i,6)==0 %不含變壓器的支路 p=B1(i,1)。 Y(p,q)=Y(p,q)1/(B1(i,3)*B1(i,5))。for i=1:n1 if B2(i,5)==2 m=m+1。 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。%雅克比矩陣7*7 %————————————————————————————————————
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1