diff --git a/OPF.m b/OPF.m index 8728029..22241d9 100644 --- a/OPF.m +++ b/OPF.m @@ -11,7 +11,7 @@ nodeMaxDVolt=zeros(33,1); nodeMaxDVAngle=zeros(33,1); maxDVoltPos=0; maxDAnglePos=0; -for badDataNode=1:33 aaa +for badDataNode=1:33 loopN=1; maxDVolt=0; maxDVAngle=0; @@ -61,9 +61,9 @@ for badDataNode=1:33 aaa maxDVAngle_=max(abs((UAngel(2:33)-rUAngel(2:33)))); if maxDVAngle_>maxDVAngle maxDVAngle=maxDVAngle_; - maxDVoltPos_=find(maxDVolt==abs((UAngel(2:33)-rUAngel(2:33)))); - if isempty(maxDVoltPos_)==0 - maxDVoltPos=maxDVoltPos_; + maxDVAnglePos_=find(maxDVAngle==abs((UAngel(2:33)-rUAngel(2:33)))); + if isempty(maxDVAnglePos_)==0 + maxDVAnglePos=maxDVAnglePos_; end end % nodeMaxDVolt_t=abs((rVolt-Volt))';