diff --git a/OPF.m b/OPF.m
index 8319786..e7321e4 100644
--- a/OPF.m
+++ b/OPF.m
@@ -17,7 +17,7 @@ for badDataNode=1:1
         close
         [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,Branchi,Branchg,Branchb,Transfork0]= ...
             pf('E:\算例\feeder33\feeder33.txt');
-        sigma=0.001;
+        sigma=0.01;
         RealPD=PD;
         RealQD=QD;
         rVolt=Volt;
@@ -39,7 +39,7 @@ for badDataNode=1:1
 %         mVolt(badDataNode)=rVolt(badDataNode)*(1-sigma*6);
         %% 画Case A的图
         % figure('Color',[1 1 1]);
-        [Busnum,Loadi,Volt,PD,QD,rVolt,UAngel,RealPD,RealQD,rUAngel,Vbi,PDbi,QDbi,plotGapA,isConverge]=subOPF([1:33],PD0,QD0,mVolt,sigma);%全部有
+        [Busnum,Loadi,Volt,PD,QD,rVolt,UAngel,RealPD,RealQD,rUAngel,Vbi,PDbi,QDbi,plotGapA,isConverge]=subOPF([22],PD0,QD0,mVolt,sigma);%全部有
         %         badDataResult(I,badDataNode)=sum(Vbi);
         %         badDataLocation(1:33,I)=Vbi;
         %         badDataLocation(34,I)=sum(abs((rVolt-Volt)./rVolt./length(rVolt)))+sum(abs( (UAngel(2:33)-rUAngel(2:33))./rUAngel(2:33)./length(rUAngel(2:33))));
diff --git a/openfile2.m b/openfile2.m
index 53952b6..10ac88e 100644
--- a/openfile2.m
+++ b/openfile2.m
@@ -56,6 +56,12 @@ QD=QD/Base;
 PD=sparse(PD);
 QD=sparse(QD);
 %QD=PD*sqrt(1-.85^2)/.85;
+
+%%
+%增加大负荷
+PD(22)=PD(22)*65;
+%%
+
 PG=sparse(PG);
 QG=sparse(QG);
 %% pv节点功率参数矩阵
diff --git a/subOPF.m b/subOPF.m
index a6fee55..04591da 100644
--- a/subOPF.m
+++ b/subOPF.m
@@ -46,7 +46,7 @@ Gap=(Init_L*Init_Z'-Init_U*Init_W');
 KK=0;
 plotGap=zeros(1,60);
 ContrlCount=size(Loadi,1)*2+Busnum*2+Busnum+length(Loadi)*2;
-kmax=150;
+kmax=350;
 Precision=Precision/1;
 %% 加误差
 %找DG
@@ -88,6 +88,8 @@ uQD(DGi)=abs(mQD(DGi)*3*sigma*0.3);
 lVolt(DGi)=abs(mVolt(DGi)'*3*sigma*0.3);
 uVolt(DGi)=abs(mVolt(DGi)'*3*sigma*0.3);
 %% 不测量的数据
+% mPD(22)=mPD(22)/65;
+
 lPD(noMeasurei)=0.15*mPD(noMeasurei);%15%
 uPD(noMeasurei)=0.15*mPD(noMeasurei);
 lQD(noMeasurei)=0.15*mQD(noMeasurei);
@@ -95,10 +97,10 @@ uQD(noMeasurei)=0.15*mQD(noMeasurei);
 
 
 %电压伪量测错误
-mVolt(noMeasurei)=0.95;
+mVolt(noMeasurei)=.95;
 
-lVolt(noMeasurei)=0.7*mVolt(noMeasurei);%电压在0.93~1.07
-uVolt(noMeasurei)=0.7*mVolt(noMeasurei);
+lVolt(noMeasurei)=0.07*mVolt(noMeasurei);%电压在0.93~1.07
+uVolt(noMeasurei)=0.07*mVolt(noMeasurei);
 
 
 
@@ -138,7 +140,7 @@ while(abs(Gap)>Precision*10)
     %% 形成方程矩阵
     Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1);
     Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1);
-    bigM=1;
+    bigM=0.4;
     Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD,mQD,uPD,lPD,uQD,lQD,uVolt,lVolt);
     Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi);
     Ly=Mat_H;
@@ -191,8 +193,8 @@ while(abs(Gap)>Precision*10)
     %% 形成方程矩阵
     Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1);
     Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1);
-    bigM=1;
-    eps=Gap*0.1;
+    bigM=0.4;
+    eps=Gap*0.001;
     Mat_G=FormG(Volt,PD,QD,Loadi,bigM,mVolt,rVolt,sigma,Vbi,PDbi,QDbi,mPD,mQD,uPD,lPD,uQD,lQD,uVolt,lVolt);
     Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi);
     Ly=Mat_H;