diff --git a/FormLwstate1.m b/FormLwstate1.m
index 6e11624..bc33469 100644
--- a/FormLwstate1.m
+++ b/FormLwstate1.m
@@ -8,10 +8,10 @@ pIi=find(Ii>0);
 nIi=find(Ii<0);
 % lower=-0.2*sparse(ones(length(Loadi)*2,1));
 upper=ones(length(Loadi)*2,1);
-upper(pIr)=1.2*Ir(pIr);
-upper(nIr)=0.8*Ir(nIr);
-upper(pIi+length(Ir))=1.2*Ii(pIi);
-upper(nIi+length(Ir))=0.8*Ii(nIi);
+upper(pIr)=1.1*Ir(pIr);
+upper(nIr)=0.9*Ir(nIr);
+upper(pIi+length(Ir))=1.1*Ii(pIi);
+upper(nIi+length(Ir))=0.9*Ii(nIi);
 
 % upper([4,5,6])=[1;1;1];
 % upper([4,5,6])=0.2*Ii(nIi);
diff --git a/FormLzstate1.m b/FormLzstate1.m
index dce009d..103dd3e 100644
--- a/FormLzstate1.m
+++ b/FormLzstate1.m
@@ -7,10 +7,10 @@ pIi=find(Ii>0);
 nIi=find(Ii<0);
 % lower=-0.2*sparse(ones(length(Loadi)*2,1));
 lower=ones(length(Loadi)*2,1);
-lower(pIr)=0.8*Ir(pIr);
-lower(nIr)=1.2*Ir(nIr);
-lower(pIi+length(Ir))=0.8*Ii(pIi);
-lower(nIi+length(Ir))=1.2*Ii(nIi);
+lower(pIr)=0.9*Ir(pIr);
+lower(nIr)=1.1*Ir(nIr);
+lower(pIi+length(Ir))=0.9*Ii(pIi);
+lower(nIi+length(Ir))=1.1*Ii(nIi);
 
 % lower=-ones(length(Ir)*2,1);
 % lower([4,5,6])=[-1;-1;-1];
diff --git a/run.m b/run.m
index 6f0c2e3..356f77b 100644
--- a/run.m
+++ b/run.m
@@ -148,7 +148,7 @@ ub=checkSSatisfied(Balance,phaseABCY,VoltpABC, ...
     phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ );
 fprintf('最大不平衡量为%f\n\n',full(max(abs(ub))))
 %% 开始进入状态估计
-clear PD QD PG QG;
+% clear PD QD PG QG;
 %准备量测量
 iterPhaseASpotLoadP=phaseASpotLoadP;
 iterPhaseBSpotLoadP=phaseBSpotLoadP;
@@ -157,6 +157,9 @@ iterPhaseASpotLoadQ=phaseASpotLoadQ;
 iterPhaseBSpotLoadQ=phaseBSpotLoadQ;
 iterPhaseCSpotLoadQ=phaseCSpotLoadQ;
 %全部转换为负荷电流
+VoltpA=sparse(ones(busNum,1));
+VoltpB=sparse(ones(busNum,1)).*exp(1j*-120/180*pi);
+VoltpC=sparse(ones(busNum,1)).*exp(1j*+120/180*pi);
 CurpA=conj((iterPhaseASpotLoadP+1j*iterPhaseASpotLoadQ)./VoltpA);
 CurpB=conj((iterPhaseBSpotLoadP+1j*iterPhaseBSpotLoadQ)./VoltpB);
 CurpC=conj((iterPhaseCSpotLoadP+1j*iterPhaseCSpotLoadQ)./VoltpC);
@@ -249,9 +252,12 @@ 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']);
+SEIf0=I0r+1j*I0i;
+SEIf1=I1r+1j*I1i;
+SEIf2=I2r+1j*I2i;
+SEIpABC=full(Tf2p*conj([SEIf0';SEIf1(1:3)';SEIf2']));
 %看一下差多少
-(full(abs(VoltpABC))-abs(SEVoltpABC))./abs(SEVoltpABC)*100
\ No newline at end of file
+(full(abs(VoltpABC))-abs(SEVoltpABC))./abs(SEVoltpABC)*100
+
+% 三相负荷
+SEVoltpABC(:,setxor(1:size(SEVoltpABC,2),Balance)).*conj(-SEIpABC)
\ No newline at end of file