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

平衡节点电压固定 

Signed-off-by: facat <dmy@dmy-PC.(none)>
This commit is contained in:
facat 2012-05-22 11:41:52 +08:00
parent cb34aec845
commit 703802da3f
4 changed files with 16 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Modification.m
View File

@ -17,6 +17,8 @@ t=deltX(size(PVi,1)+size(PGi,1)+1:ContrlCount)';
t(Busnum+Balance)=0;
%Volt=Volt+AlphaP*t(2:2:2*Busnum);暂时改一下20111227
%UAngel=UAngel+AlphaP*t(1:2:2*Busnum);暂时改一下20111227
balVolt=Volt(Balance);
Volt=Volt+AlphaP*t(1:Busnum);
Volt(Balance)=balVolt;
UAngel=UAngel+AlphaP*t(Busnum+1:2*Busnum);
end

6
OPF.m
View File

@ -1,12 +1,12 @@
tic
clear
%[kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB,Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,LineLimti,LineLimtj,LinePLimt,PG,QG,PD,QD,CenterA,LineCount,PGi,PVQU,PVQL]=pf('5sj.txt');
[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]=pf('ieee3001.dat');
[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]=pf('ieee10471.dat');
GB=full(GB);
%PVi电压节点序号
%PVu电压节点电压标幺值
Volt
UAngel*180/3.1415926
Volt;
UAngel*180/3.1415926;
%sprintf('%f\n',Volt);
%sprintf('%f\n',Angel);
%% 初值

4
OPF_Init.m
View File

@ -2,8 +2,8 @@ function [Volt,UAngel,Init_Z,Init_W,Init_L,Init_U,Init_Y,PG,QG,RestraintCount]=O
RestraintCount=size(PVi,1)+size(PGi,1)+Busnum; %约束条件数
t_Bal_volt=Volt(Balance);
Volt=ones(1,Busnum);
%Volt(Balance)=t_Bal_volt;
Volt(Balance)=1;
Volt(Balance)=t_Bal_volt;
%Volt(Balance)=1;
UAngel=zeros(1,Busnum);
Init_Z=ones(1,RestraintCount);
Init_W=-1*ones(1,RestraintCount);

14
SolveIt.m
View File

@ -8,15 +8,19 @@ aa=[
deltH',zeros(size(Init_Y,2));
];
yy=[LxComa;-Ly];
% t=size(PVi,1)+size(PGi,1);
% aa(t+2*Balance-1,:)=0;
% aa(:,t+2*Balance-1)=0;
% aa(t+2*Balance-1,t+2*Balance-1)=1;
%暂时改一下
%% 平衡节点电压不变
t=size(PVi,1)+size(PGi,1);
aa(t+Balance,:)=0;
aa(:,t+Balance)=0;
aa(t+Balance,t+Balance)=1;
deltG(t+Balance,:)=0;
%%
t=size(PVi,1)+size(PGi,1)+Busnum;
aa(t+Balance,:)=0;
aa(:,t+Balance)=0;
aa(t+Balance,t+Balance)=1;
deltG(t+Balance,:)=0;
%%
dxdy=aa\yy;
dX=dxdy(1:ContrlCount);
dY=dxdy(ContrlCount+1:ContrlCount+2*Busnum);