,')次差值:'];个人收集整理勿做商业用途%----------------求取各个节点地功率及功率偏差及PV节点地电压偏差--------------------个人收集整理勿做商业用途fori=1:n%n个节点2n行(每节点两个方程P和Q或U)p=2*i-1;m=p+1;C(i)=0;D(i)=0;forj1=1:n%第i行共n列(n个节点间互导纳及节点电压相乘即电流)C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)个人收集整理勿做商业用途D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bij*ej)个人收集整理勿做商业用途end%求i节点有功和无功功率P',Q'地计算值P1=C(i)*e(i)+f(i)*D(i);%节点功率P计算eiΣ(Gij*ej-Bij*fj)+fiΣ(Gij*fj+Bij*ej)个人收集整理勿做商业用途Q1=C(i)*f(i)-e(i)*D(i);%节点功率Q计算fiΣ(Gij*ej-Bij*fj)-eiΣ(Gij*fj+Bij*ej)个人收集整理勿做商业用途V2=e(i)^2+f(i)^2; %电压模平方%===求取功率差及PV节点电压模平方差=========ifi~=isb %非平衡节点(PQ或PV节点)ifB2(i,5)~=3 %非PV节点(只能是PQ节点)J(m,N1)=P(i)-P1; %PQ节点有功功率差J(m,N1)扩展列△PJ(p,N1)=Q(i)-Q1;%PQ节点无功功率差J(p,N1)扩展列△Qelse%PV节点==================J(m,N1)=P(i)-P1; %PV节点有功功率差J(m,N1)扩展列△PJ(p,N1)=V(i)^2-V2; %PV节点电压模平方差J(p,N1)扩展列△Uend
全文预览
