diff --git a/@Opti/obj.m b/@Opti/obj.m
index 1698e78..57abde4 100644
--- a/@Opti/obj.m
+++ b/@Opti/obj.m
@@ -26,6 +26,8 @@ Y=this.Y;
 VMatrix=sparse(r,c,SEVAngel(r)-SEVAngel(c)-Angel,busNum,busNum);
 PD_=PG-diag(SEVolt)*(Y.*cos(VMatrix))*SEVolt;
 QD_=QG-diag(SEVolt)*(Y.*spfun(@sin,VMatrix))*SEVolt;
+% wPD=ones(length(wPD),1);
+% wQD=ones(length(wPD),1);
 t4=(wPD(PDi).*(PD_(PDi)-PD0(PDi))).^2;
 t5=(wQD(QDi).*(QD_(QDi)-QD0(QDi))).^2;
 % loadCurrent=LoadCurrent( SEVolt,SEVAngel,PD0(PDi),QD0(QDi),PDi,QDi );
@@ -34,7 +36,7 @@ t5=(wQD(QDi).*(QD_(QDi)-QD0(QDi))).^2;
 % t7=(wLoadCurrent.*(loadCurrent-mLoadCurrent)).^2;
 % t7=wLoadCurrent.*((loadCurrent-mLoadCurrent)./mLoadCurrent).^2;
 % output_args=sum(t4)+sum(t5)+sum(t6);
-output_args=sum(t6)+sum(t4)+sum(t5);
+output_args=sum(t6)+sum(t5)+sum(t4);
 output_args=full(output_args);
 
 
diff --git a/OPF_Init.m b/OPF_Init.m
index 9a97f3d..b75e6bc 100644
--- a/OPF_Init.m
+++ b/OPF_Init.m
@@ -37,7 +37,7 @@ wQG=1*ones(size(PVi,1),1);
 %randPDind=randInt(1:10);
 randPDind=0;
 % load('I');
-sigma=0.03;
+sigma=0.01;
 % sigma=sigma*I/200;
 % wPD=1/(1)*ones(length(PD),1);
 % wPD(11)=wPD(11)*1;
diff --git a/Run_YALMIP2.m b/Run_YALMIP2.m
index d5b4310..ebc1270 100644
--- a/Run_YALMIP2.m
+++ b/Run_YALMIP2.m
@@ -9,7 +9,7 @@ for I=1:1
     close all;
     [kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB, ...
         Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,PG,QG,PD,QD,CenterA,PGi,PVQU,PVQL, ...
-        Liner,Linex,Lineb,Transforr,Transforx,Transfork0]=pf('E:/算例/东际911_2751267_2012-09-05/newFIle20-使用.txt');
+        Liner,Linex,Lineb,Transforr,Transforx,Transfork0]=pf('E:\算例\feeder33\feeder33.txt');
 %     'E:\算例\feeder33\feeder33.txt'
     %pf('C:\bpa\长虹世纪线_556844_2013-10-28\长虹世纪线_556844_2013-10-28_iPso_newFil
     %e.txt');
@@ -68,7 +68,7 @@ for I=1:1
 %     QD0(11)=(1.+0.01*I)*rQD(11);
     % % load('mLoadCurrent');
     % load('I');
-    sigma=0.03;
+    sigma=0.1;
     %
     [ mVolt,PD0,QD0,mLoadCurrent ] = GetMeasure(sigma,rVolt,rPD,rQD,rLoadCurrent );
     % sigma=sigma*(I-1)/50;
@@ -90,12 +90,13 @@ for I=1:1
     opts.maxtime=30000;
     opts.maxfeval=85000;
     opts.maxnodes=85000;
-    opts.tolrfun=1e-4;
-    opts.tolafun=1e-4;
+    opts.tolrfun=1e-3;
+    opts.tolafun=1e-3;
     opts.warnings='all';
     opts.display='off';
-    x0=[ones(length(Volt),1); ...
-        zeros(length(Volt),1)];
+%     x0=[0.95*ones(length(Volt),1); ...
+%         zeros(length(Volt),1)];
+    x0=[rVolt;rVAngel];
     % x0=[PD(PDi);QD(QDi);xVolt';xUAngel'];
     [~,seOpti]=seOpti.equ(x0);
     cl=seOpti.Getcl();