diff --git a/FormLwstate1.m b/FormLwstate1.m
index 6e11624..d552ff4 100644
--- a/FormLwstate1.m
+++ b/FormLwstate1.m
@@ -13,6 +13,14 @@ upper(nIr)=0.8*Ir(nIr);
 upper(pIi+length(Ir))=1.2*Ii(pIi);
 upper(nIi+length(Ir))=0.8*Ii(nIi);
 
+%太小的数都要放宽一些
+tooSmall=find(abs(Ir)<0.05);
+% upper(tooSmall)=0.99*abs(Ir(tooSmall));
+upper(tooSmall)=0.2;
+tooSmall=find(abs(Ii)<0.05);
+% upper(tooSmall+length(Ir))=0.99*abs(Ii(tooSmall));
+upper(tooSmall+length(Ir))=0.2;
+
 % 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..2adaea9 100644
--- a/FormLzstate1.m
+++ b/FormLzstate1.m
@@ -12,6 +12,14 @@ lower(nIr)=1.2*Ir(nIr);
 lower(pIi+length(Ir))=0.8*Ii(pIi);
 lower(nIi+length(Ir))=1.2*Ii(nIi);
 
+%太小的数都要放宽一些
+tooSmall=find(abs(Ir)<0.05);
+% lower(tooSmall)=-0.99*abs(Ir(tooSmall));
+lower(tooSmall)=-0.2;
+tooSmall=find(abs(Ii)<0.05);
+% lower(tooSmall+length(Ir))=-.99*abs(Ii(tooSmall));
+lower(tooSmall+length(Ir))=-0.2;
+
 % lower=-ones(length(Ir)*2,1);
 % lower([4,5,6])=[-1;-1;-1];
 Lz=Mat_G-Init_L-lower;
diff --git a/feeder13/data1.txt b/feeder13/data1.txt
index a873457..1c37398 100644
--- a/feeder13/data1.txt
+++ b/feeder13/data1.txt
@@ -6,9 +6,9 @@
 0
 632 0 0 0
 0
-630 80 70 60 50 40 30
-631 40 10 11 50 30 30
-632 120 120 120 110 90 90
+630 80 70 80 60 80 70
+631 40 10 40 10 42 10
+632 120 120 120 120 120 120
 0
 
 
diff --git a/run.m b/run.m
index 9c845c9..05910a9 100644
--- a/run.m
+++ b/run.m
@@ -238,8 +238,8 @@ wI2i=abs(imag(I2measurement)).*sigma;
 % f=sum(([real(I2measurement);imag(I2measurement)]-[I2r;I2i]).^2)+sum((real(rV2)-V2r).^2)+sum((imag(rV2)-V2i).^2);
 % fprintf('目标值 %f\n',full(f));
 %% 算零序
-BalI0r=real(-sum(mIf0));
-BalI0i=imag(-sum(mIf0));
+BalI0r=real(-sum(rIf0));
+BalI0i=imag(-sum(rIf0));
 %电压
 %制作量测量
 % mfsY00=fsY00;
diff --git a/theoryTest.m b/theoryTest.m
new file mode 100644
index 0000000..154f306
--- /dev/null
+++ b/theoryTest.m
@@ -0,0 +1,30 @@
+clc
+clear
+a=exp(1j*2*pi/3);
+Tp2f=1/3*[1 1 1;
+    1 a a^2;
+    1 a^2 a];
+sigma=0.03;
+rLoadP=[10;16;12];
+rLoadQ=[9;9;9];
+eLoadP=normrnd(0,sigma,length(rLoadP),1);
+mLoadP=rLoadP.*(1+eLoadP);
+eLoadQ=normrnd(0,sigma,length(rLoadQ),1);
+mLoadQ=rLoadQ.*(1+eLoadQ);
+rVoltP=[1;1*exp(1j*-120/180*pi);1*exp(1j*+120/180*pi)];
+mVoltP=rVoltP.*(1+normrnd(0,sigma,length(rVoltP),1));
+rI=conj((rLoadP+1j*rLoadQ)./rVoltP);
+mI=conj((mLoadP+1j*mLoadQ)./mVoltP);
+rFI=Tp2f*rI
+mFI=Tp2f*mI
+%计算不平衡度
+LmP=(sum(mLoadP.^4))/(sum(mLoadP.^2)*sum(mLoadP.^2));
+LmQ=(sum(mLoadQ.^4))/(sum(mLoadQ.^2)*sum(mLoadQ.^2));
+ubmP=sqrt( (1-sqrt(3-6*LmP))/(1+sqrt(3-6*LmP) ) )*100
+ubmQ=sqrt( (1-sqrt(3-6*LmQ))/(1+sqrt(3-6*LmQ) ) )*100
+LrP=(sum(rLoadP.^4))/(sum(rLoadP.^2)*sum(rLoadP.^2));
+LrQ=(sum(rLoadQ.^4))/(sum(rLoadQ.^2)*sum(rLoadQ.^2));
+ubrP=sqrt( (1-sqrt(3-6*LrP))/(1+sqrt(3-6*LrP) ) )*100
+ubrQ=sqrt( (1-sqrt(3-6*LrQ))/(1+sqrt(3-6*LrQ) ) )*100
+ramI=conj((rLoadP+1j*rLoadQ)./mVoltP)+conj((rLoadP.*eLoadP+1j*rLoadQ.*eLoadQ)./mVoltP);
+ramFI=Tp2f*conj((rLoadP+1j*rLoadQ)./mVoltP)+Tp2f*conj((rLoadP.*eLoadP+1j*rLoadQ.*eLoadQ)./mVoltP);
\ No newline at end of file