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