facat
/
pes2014-wronglowervoltagebound
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-03-11 20:53:51 +08:00
parent 1f9f3ff999
commit 66fd308278
2 changed files with 20 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
OPF.m
View File

@ -9,7 +9,9 @@ badDataResult=zeros(10,33);
badDataLocation=zeros(34,10);
nodeMaxDVolt=zeros(33,1);
nodeMaxDVAngle=zeros(33,1);
for badDataNode=1:1
maxDVoltPos=0;
maxDAnglePos=0;
for badDataNode=1:33
loopN=1;
maxDVolt=0;
maxDVAngle=0;
@ -51,20 +53,28 @@ for badDataNode=1:1
maxDVolt_=max(abs((rVolt-Volt)));
if maxDVolt_>maxDVolt
maxDVolt=maxDVolt_;
maxDVoltPos_=find(maxDVolt==abs((rVolt-Volt)));
if isempty(maxDVoltPos_)==0
maxDVoltPos=maxDVoltPos_;
end
end
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_;
end
end
nodeMaxDVolt_t=abs((rVolt-Volt))';
nodeMaxDVAngle_t=abs((UAngel-rUAngel))';
nodeMaxDVolt_t([1:18,23:33])=0;
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);
% nodeMaxDVolt_t=abs((rVolt-Volt))';
% nodeMaxDVAngle_t=abs((UAngel-rUAngel))';
% nodeMaxDVolt_t([1:18,23:33])=0;
% 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;
nodeMaxDVolt(badDataNode)=maxDVolt;
nodeMaxDVAngle(badDataNode)=maxDVAngle;
break;
end

2
subOPF.m
View File

@ -237,7 +237,7 @@ dQ(Balance)=0;
mdP=max(dP);
mdQ=max(dQ);
%判断是否收敛
if min(abs(Vbi-.1))>5e-005
if min(abs(Vbi-.1))>0.1+5e-005
isConverge=0;
end
if KK>=kmax