牛拉潮流程序-資料下載頁
2025-01-09 22:00本頁面
【正文】 9。,num2str(num),39。 次時的Delta總的失配量為：39。]) % DEL %繼續(xù)判斷最大值 MAX =max(abs(DEL))。 Theta_first=Theta_first39。end %%%%%%%%%%%%%%%%%%求平衡節(jié)點的有功功率和無功功率%%%%%%%%%%%%%%%%%%%%%%%%%%%Ps0=0。i=n。for t=1:n theij=Theta_first(i)Theta_first(t)。 Ps0=Ps0+U_first(t)*(G(i,t)*cos(theij)+B(i,t)*sin(theij))。endPs0=U_first(i)*Ps0。Z(i,7)=Ps0。Qs0=0。i=n。for t=1:n theij=Theta_first(i)Theta_first(t)。 Qs0=Qs0+U_first(t)*(G(i,t)*sin(theij)B(i,t)*cos(theij))。endQs0=U_first(i)*Qs0。Z(i,8)=Qs0。%%%%%%%%%%%%%%%%%%%%%計算PV節(jié)點的無功功率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Qv0=zeros(1,n)。for i=nm:n1 for t=1:n theij=Theta_first(i)Theta_first(t)。 Qv0(i)=Qv0(i)+U_first(t)*(G(i,t)*sin(theij)B(i,t)*cos(theij))。 endQv0(i)=U_first(i)*Qv0(i)。endfor i=nm:n1 Z(i,8)=Qv0(i)。endZj=sqrt(1)。for i=1:n Angle(i)=Theta_first(i)*180/pi。 P_in(i)=Z(i,7)Z(i,5)+j*(Z(i,8)Z(i,6))。endfor i=1:n Flow(i,1)=Z(i,1)。 Flow(i,2)=U_first(i)。 Flow(i,3)=Angle(i)。 Flow(i,4)=P_in(i)。endFlow(:,1)=Z2(:,1)。disp(39。節(jié)點電壓和注入功率39。)disp(39。 節(jié)點號 節(jié)點電壓幅值 節(jié)點電壓相位 節(jié)點注入功率39。)Flowdisp(39。=======================完*************成======================39。)%Bus 信息 %1編號 2類型 3U *4 5 P(L) 6Q(L) 7P(G) 8Q(G) 9Uo 10Qmax 11Qmin 12(串并聯(lián)電容電抗器)NODE=[1 0 0 0 0 1 0 0 0 0。 2 0 0 0 0 1 0 0 0 0。 3 0 0 0 0 0 0 1 0 0 0 0。 4 2 0 0 0 0 0 0 0 0。 5 3 0 0 0 0 0 0 0 0 0。]%% Branch 信息 %i j R X B K(1)Branch=[1 2 。 1 3 。 2 3 。 4 2 。 5 3 。] %第 9 頁 共 9 頁
公司管理相關(guān)推薦
鄂ICP備17016276號-1