energy.asv

可进行电力系统多节点系统的优化潮流计算

function [enrgy baseMVA, V, gen, branch]=E(Qg,QL,TAP)

[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
    VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;

[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, ...
    MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] = idx_gen;

[F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, ...
    RATE_C, TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST] = idx_brch;
%%%%%%%%%%%newton%%%%%%%%%%%%%%%%%%%%
if nargin>10
end
for i=1:NNU
     for j=1:NNU
    f1(i,energy)=f1(i,enrgy)-1/2*BB(i,j)*abs(Vs(i))*abs(Vs(j))*cos(deltas(i)-deltas(j))...
 +1/2*BB(i,j)*abs(Vu(i))*abs(Vu(j))*cos(deltau(i)-deltau(j))...
-BG(i,j)*abs(Vs(i))*abs(Vs(j))*cos(deltas(i)-deltas(j))*(deltau(i)-deltas(i))...
 -((Vu(i)-Vs(i))/Vs(i))*BG(i,j)*Vs(i)*Vs(j)*sin(deltas(i)-deltas(j))...
 +1/2*BC(i0)*(Vs(i0)^2-Vu(i0)^2);
     else
 f1(i)=f1(i)-1/2*BB(i,j)*abs(Vs(i))*abs(Vs(j))*cos(deltas(i)-deltas(j))...
 +1/2*BB(i,j)*abs(Vu(i))*abs(Vu(j))*cos(deltau(i)-deltau(j))...
-BG(i,j)*abs(Vs(i))*abs(Vs(j))*cos(deltas(i)-deltas(j))*(deltau(i)-deltas(i))...
-((Vu(i)-Vs(i))/Vs(i))*BG(i,j)*Vs(i)*Vs(j)*sin(deltas(i)-deltas(j));
 %%%%vq是全局能量函数%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %%%%%%%%%%%%%%%%%%%%%全局能量函数%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 end
 end
 end
 v(i)=-PL(i)/SB*(deltau(i)-deltas(i))+f1(i)-QL(i)/SB*log2(Vu(i)/Vs(i));%%%%v(i)是单个节点的能量函数
 vq=vq+v(i);
end
Qg=
Ybus;
S
[V, converged, i] = newtonpf(Ybus, Sbus, V0, ref, pv, pq, mpopt)

Qg= 

QL=

TAP=
%%%%%%%%%%%%%%%%%%%change Qg%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
newton;

%%%%%%%%%%%%%%%%%%%%%%change QL%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
newton;
newtonpf;
%%%%%%%%%%%%change transformer tap%%%%%%%%%%%%%%%%%
energy