From 708115e7f929ddb0b8efc7bd4b5f47d5057e0d19 Mon Sep 17 00:00:00 2001 From: "facat@lab.com" Date: Mon, 3 Nov 2014 17:14:18 +0800 Subject: [PATCH] =?UTF-8?q?=E5=88=A0=E6=8E=89=E5=A4=9A=E4=BD=99=E6=B3=A8?= =?UTF-8?q?=E9=87=8A?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: facat@lab.com --- OPF.m | 22 ---------------------- 1 file changed, 22 deletions(-) diff --git a/OPF.m b/OPF.m index a42e3f1..367410a 100644 --- a/OPF.m +++ b/OPF.m @@ -75,27 +75,6 @@ while(k<=kmax && maxD> EPS) iterPhaseASpotLoadQ=phaseASpotLoadQ+imag(SA); iterPhaseBSpotLoadQ=phaseBSpotLoadQ+imag(SB); iterPhaseCSpotLoadQ=phaseCSpotLoadQ+imag(SC); -% %全部转换为负荷电流 -% CurpA=conj((iterPhaseASpotLoadP+1j*iterPhaseASpotLoadQ)./VoltpA); -% CurpB=conj((iterPhaseBSpotLoadP+1j*iterPhaseBSpotLoadQ)./VoltpB); -% CurpC=conj((iterPhaseCSpotLoadP+1j*iterPhaseCSpotLoadQ)./VoltpC); -% %转换为序电流 -% f012=Tp2f*conj([CurpA';CurpB';CurpC']); -% %把三序电流分离出来 -% If0=conj(f012(1,:)'); -% If1=conj(f012(2,:)'); -% If2=conj(f012(3,:)'); -% %试着算一下正序电流 -% fsY11*V1; -% %形成负荷序电流的测量值 -% mIf0=If0; -% mIf1=If1; -% mIf1(3)=-mIf1(2); -% mIf2=If2; -% %计算 -% fsY11=fsY11+sparse(Balance,Balance,ones(length(Balance),1),busNum,bus -% Num);%这里要置0，置1，否则是奇异的 - %%做最小二乘法 [dP, dQ, YdotSinVolt, YdotCosVolt, diag_Volt_YdotSin, diag_Volt_YdotCos]=Unbalance(Balance,busNum, ... PQi,PG,QG,QGi,iterPD,iterQD,Vmf1,Vaf1,fsY1amp,fsY1ang,r,c,Vf2,If2,Vf0,If0);%不平衡量 maxD=max(abs([dP;dQ;])); @@ -217,7 +196,6 @@ while(abs(Gap)>Precision) if KK>kmax break; end -% plotGap(KK+1)=Gap; Init_u=Gap/2/RestraintCount*CenterA; AngleIJMat=0; %% 开始计算OPF