diff --git a/OPF.m b/OPF.m
index b1f4374..71e521b 100644
--- a/OPF.m
+++ b/OPF.m
@@ -225,7 +225,8 @@ wQD=1./(abs(mQD3P*sigma).^2);
 
 %% 没有量测量
 % noLoadi=[1,5,6,10,11];
-noLoadi=Loadi;
+noLoadi=[1,8,9,10,12];
+% noLoadi=Loadi;
 % noLoadi=[1,11];
 % noLoadi=[];
 noPQi3P=zeros(length(noLoadi)*3,1);
@@ -234,9 +235,9 @@ noPQi3P(2:3:end)=(noLoadi-1)*3+2;
 noPQi3P(3:3:end)=(noLoadi-1)*3+3;
 % noLoadi=noPQi3P;
 wVolt(setdiff(1:length(wVolt),Loadi))=0;%只有负荷处才有电压量测。
-wVolt( noLoadi)=0;
-wPD(ismember( Loadi,noLoadi))=0;
-wQD(ismember( Loadi,noLoadi))=0;
+wVolt( noPQi3P)=0;
+wPD(ismember( Loadi,noPQi3P))=0;
+wQD(ismember( Loadi,noPQi3P))=0;
 %%
 % RestraintCount=size(Loadi,1)*2+length(rVoltABCV); %约束条件数
 RestraintCount=size(Loadi,1)*2; %约束条件数