From 7409723aae502864fa67d123c08b40ed53da8665 Mon Sep 17 00:00:00 2001
From: "dugg@lab-desk" <dugg@lab-desk>
Date: Wed, 29 Oct 2014 21:26:39 +0800
Subject: [PATCH] =?UTF-8?q?=E6=B7=BB=E5=8A=A0=E6=8A=8A=E7=BB=93=E6=9E=9C?=
 =?UTF-8?q?=E8=BD=AC=E6=8D=A2=E4=BC=9A=E4=B8=89=E7=9B=B8=E5=80=BC=20Signed?=
 =?UTF-8?q?-off-by:=20dugg@lab-desk=20<dugg@lab-desk>?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 run.m | 13 ++++++++++++-
 1 file changed, 12 insertions(+), 1 deletion(-)

diff --git a/run.m b/run.m
index bbe94f5..ae61756 100644
--- a/run.m
+++ b/run.m
@@ -298,8 +298,8 @@ sigmaI0=normrnd(0,sigma,length(Loadi),1);
 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);
 % fprintf('目标值 %f\n',full(f));
 %状态量
@@ -319,3 +319,14 @@ wI0i=abs( imag(mIf0(Loadi).*sigmaI0) );
 % SEQD=sparse(zeros(busNum,1));
 
 %检查目标函数
+%% 把三序合成三相
+% 三相电压
+SEVf0=V0r+1j*V0i;
+SEVf1=V1r+1j*V1i;
+SEVf2=V2r+1j*V2i;
+SEVoltpABC=Tf2p*conj([ SEVf0'; SEVf1'; SEVf2']);
+% 三序电流
+SEIf0=I0r+1j*I0r;
+SEIf1=I1r+1j*I1r;
+SEIf2=I2r+1j*I2r;
+SEIpABC=Tf2p*conj([SEIf0';SEIf1';SEIf2']);
\ No newline at end of file