facat
/
pes2014
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

测试电压幅值下降错误的情况。 

Signed-off-by: dugg@lab-desk <dugg@lab-desk>
This commit is contained in:
dugg@lab-desk 2015-02-10 15:55:19 +08:00
parent 1c1e8f4b3c
commit 5a752ca4fb
2 changed files with 17 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
OPF.m
View File

@ -43,7 +43,7 @@ for badDataNode=1:1
% badDataResult(I,badDataNode)=sum(Vbi); % badDataResult(I,badDataNode)=sum(Vbi);
% badDataLocation(1:33,I)=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)))); % 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))));
break; % break;
if isConverge==0 if isConverge==0
continue; continue;
end end
@ -56,10 +56,15 @@ for badDataNode=1:1
if maxDVAngle_>maxDVAngle if maxDVAngle_>maxDVAngle
maxDVAngle=maxDVAngle_; maxDVAngle=maxDVAngle_;
end end
nodeMaxDVolt_t=abs((rVolt-Volt))';
if loopN>=1 nodeMaxDVAngle_t=abs((UAngel-rUAngel))';
nodeMaxDVolt(badDataNode)=maxDVolt; nodeMaxDVolt_t([1:18,23:33])=0;
nodeMaxDVAngle(badDataNode)=maxDVAngle; nodeMaxDVAngle_t([1:18,23:33])=0;
nodeMaxDVolt(nodeMaxDVolt<nodeMaxDVolt_t)=nodeMaxDVolt_t(nodeMaxDVolt<nodeMaxDVolt_t);
nodeMaxDVAngle(nodeMaxDVAngle<nodeMaxDVAngle_t)=nodeMaxDVAngle_t(nodeMaxDVAngle<nodeMaxDVAngle_t);
if loopN>=500
% nodeMaxDVolt(badDataNode)=maxDVolt;
% nodeMaxDVAngle(badDataNode)=maxDVAngle;
break; break;
end end

8
subOPF.m
View File

@ -49,7 +49,7 @@ Gap=(Init_L*Init_Z'-Init_U*Init_W');
KK=0; KK=0;
plotGap=zeros(1,60); plotGap=zeros(1,60);
ContrlCount=size(Loadi,1)*2+Busnum*2+Busnum+length(Loadi)*2; ContrlCount=size(Loadi,1)*2+Busnum*2+Busnum+length(Loadi)*2;
kmax=350; kmax=170;
Precision=Precision/1; Precision=Precision/1;
%% ¼ÓÎó²î %% ¼ÓÎó²î
%ÕÒDG %ÕÒDG
@ -99,6 +99,12 @@ lQD(noMeasurei)=0.15*mQD(noMeasurei);
uQD(noMeasurei)=0.15*mQD(noMeasurei); uQD(noMeasurei)=0.15*mQD(noMeasurei);
% lPD(22)=0.17*mPD(22);%15%
% uPD(22)=0.17*mPD(22);
% lQD(22)=0.17*mQD(22);
% uQD(22)=0.17*mQD(22);
%µçѹαÁ¿²â´íÎó %µçѹαÁ¿²â´íÎó
mVolt(noMeasurei)=.95; mVolt(noMeasurei)=.95;