From 4c2f16c8c7e8b07e8aeea50d2783a523e81eeac3 Mon Sep 17 00:00:00 2001
From: "dugg@lab-desk" <dugg@lab-desk>
Date: Wed, 29 Oct 2014 21:56:38 +0800
Subject: [PATCH] =?UTF-8?q?=E5=85=88=E4=B8=8D=E8=80=83=E8=99=91=E7=94=B5?=
 =?UTF-8?q?=E5=AE=B9?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
---
 run.m | 20 ++++++++++++--------
 1 file changed, 12 insertions(+), 8 deletions(-)

diff --git a/run.m b/run.m
index ae61756..d0f5942 100644
--- a/run.m
+++ b/run.m
@@ -65,6 +65,10 @@ while(k<=kmax && maxD> EPS)
     SA=VoltpA.*conj(VoltpA.*sparse(cap.capNode,1,1j*cap.capB(:,1),busNum,1));
     SB=VoltpB.*conj(VoltpA.*sparse(cap.capNode,1,1j*cap.capB(:,2),busNum,1));
     SC=VoltpC.*conj(VoltpA.*sparse(cap.capNode,1,1j*cap.capB(:,3),busNum,1));
+    %先不要电容
+    SA=0;
+    SB=0;
+    SC=0;
     iterPD=PD+real(SA+SB+SC)/3;
     iterQD=QD+imag(SA+SB+SC)/3;
     iterPhaseASpotLoadP=phaseASpotLoadP+real(SA);
@@ -231,7 +235,6 @@ mIf2=-If2;
 % fsY11(Balance,:)=0;
 % fsY11=fsY11+sparse(Balance,Balance,ones(length(Balance),1),busNum,busNum);
 % mIf1(3)=1;
-
 %% 先算正序的
 %平衡节点电流
 fprintf('正序\n');
@@ -252,10 +255,8 @@ sigmaI1=normrnd(0,sigma,length(Loadi),1);
 I1measurement=mIf1(Loadi).*(1+sigmaI1);%测量值是等效发电机电流
 wI1r=abs( real(mIf1(Loadi).*sigmaI1) );
 wI1i=abs( imag(mIf1(Loadi).*sigmaI1) );
-
 [ V1r,V1i,I1r,I1i ]=IPMLoop(V1measurement,wV1r,wV1i,I1measurement,wI1r,wI1i,BalI1r,BalI1i,busNum,Loadi,fsY11,Balance,1 );
-% f=sum(([real(Imeasurement);imag(Imeasurement)]-[-I1r;-I1i]).^2)+sum((real(Vmeasurement)-V1r(Loadi)).^2)+sum((imag(Vmeasurement)-V1i(Loadi)).^2);
-% f=sum(([real(-mIf1(Loadi));imag(-mIf1(Loadi))]-[-I1r;-I1i]).^2)+sum((real(rV1(Loadi))-V1r(Loadi)).^2)+sum((imag(rV1(Loadi))-V1i(Loadi)).^2);
+f=sum(([real(I1measurement);imag(I1measurement)]-[I1r;I1i]).^2)+sum((real(rV1)-V1r).^2)+sum((imag(rV1)-V1i).^2);
 % fprintf('目标值 %f\n',full(f));
 %% 算负序的
 fprintf('负序\n');
@@ -277,7 +278,7 @@ I2measurement=mIf2(Loadi).*(1+sigmaI2);%
 wI2r=abs( real(mIf2(Loadi).*sigmaI2) );
 wI2i=abs( imag(mIf2(Loadi).*sigmaI2) );
 [ V2r,V2i,I2r,I2i ]=IPMLoop(V2measurement,wV2r,wV2i,I2measurement,wI2r,wI2i,BalI2r,BalI2i,busNum,Loadi,fsY22,Balance,0 );
-% f=sum(([real(Imeasurement);imag(Imeasurement)]-[-I1r;-I1i]).^2)+sum((real(Vmeasurement)-V2r(Loadi)).^2)+sum((imag(Vmeasurement)-V2i(Loadi)).^2);
+f=sum(([real(I2measurement);imag(I2measurement)]-[I2r;I2i]).^2)+sum((real(rV2)-V2r).^2)+sum((imag(rV2)-V2i).^2);
 % fprintf('目标值 %f\n',full(f));
 %% 算零序
 fprintf('零序\n');
@@ -299,9 +300,9 @@ I0measurement=mIf0(Loadi).*(1+sigmaI0);%
 wI0r=abs( real(mIf0(Loadi).*sigmaI0) );
 wI0i=abs( imag(mIf0(Loadi).*sigmaI0) );
 [ V0r,V0i,I0r,I0i ]=IPMLoop(V0measurement,wV0r,wV0i,I0measurement,wI0r,wI0i,BalI0r,BalI0i,busNum,Loadi,fsY00,Balance,0 );
-%% 
-% f=sum(([real(Imeasurement);imag(Imeasurement)]-[-I1r;-I1i]).^2)+sum((real(Vmeasurement)-V0r(Loadi)).^2)+sum((imag(Vmeasurement)-V0i(Loadi)).^2);
+f=sum(([real(I0measurement);imag(I0measurement)]-[I0r;I0i]).^2)+sum((real(rV0)-V0r).^2)+sum((imag(rV0)-V0i).^2);
 % fprintf('目标值 %f\n',full(f));
+
 %状态量
 % SEVoltpA=sparse(ones(busNum,1));
 % SEVoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi);
@@ -325,8 +326,11 @@ SEVf0=V0r+1j*V0i;
 SEVf1=V1r+1j*V1i;
 SEVf2=V2r+1j*V2i;
 SEVoltpABC=Tf2p*conj([ SEVf0'; SEVf1'; SEVf2']);
+% SEVoltpABC2=Tf2p*conj([ rV0'; rV1'; rV2']);
 % 三序电流
 SEIf0=I0r+1j*I0r;
 SEIf1=I1r+1j*I1r;
 SEIf2=I2r+1j*I2r;
-SEIpABC=Tf2p*conj([SEIf0';SEIf1';SEIf2']);
\ No newline at end of file
+SEIpABC=Tf2p*conj([SEIf0';SEIf1';SEIf2']);
+%看一下差多少
+(full(abs(VoltpABC))-abs(SEVoltpABC))./abs(SEVoltpABC)*100
\ No newline at end of file