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

1.不要电容了 

2.减小线路长度，得到一个比较小的阻抗。

3.用这个数据和参数。

Signed-off-by: dmy@lab <dmy@lab.lab>
This commit is contained in:
dmy@lab 2015-04-05 22:04:02 +08:00
parent 920fcc1cd5
commit 8c4500f3db
4 changed files with 19 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
FormLw.m
View File

@ -26,8 +26,8 @@ function Lw=FormLw(Mat_G,Init_U,Loadi,rPD3P,rQD3P,mVoltABCV)
% QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF; % QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF;
% PDU=1.8*ones(length(Loadi),1); % PDU=1.8*ones(length(Loadi),1);
% QDU=1.8*ones(length(Loadi),1); % QDU=1.8*ones(length(Loadi),1);
PDU=1.25*rPD3P.*(1+unifrnd(-0.2,0.2,length(rPD3P),1)); PDU=1.3*rPD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1));
QDU=1.25*rQD3P.*(1+unifrnd(-0.2,0.2,length(rPD3P),1)); QDU=1.3*rQD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1));
% t1=([PDU',QDU',1*ones(length(Volt),1)'])'; % t1=([PDU',QDU',1*ones(length(Volt),1)'])';
t1=([PDU',QDU'])'; t1=([PDU',QDU'])';
t2=Mat_G+Init_U'-t1; t2=Mat_G+Init_U'-t1;

4
FormLz.m
View File

@ -23,8 +23,8 @@ function Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P,mVoltABCV)
% QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF; % QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF;
% PDL=0*ones(length(Loadi),1); % PDL=0*ones(length(Loadi),1);
% QDL=0*ones(length(Loadi),1); % QDL=0*ones(length(Loadi),1);
PDL=0.75*rPD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1)); PDL=0.70*rPD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1));
QDL=0.75*rQD3P.*(1+unifrnd(-0.15,0.15,length(rQD3P),1)); QDL=0.70*rQD3P.*(1+unifrnd(-0.15,0.15,length(rQD3P),1));
% t1=([PDL',QDL',0.88*ones(length(Volt),1)'])'; % t1=([PDL',QDL',0.88*ones(length(Volt),1)'])';
t1=([PDL',QDL'])'; t1=([PDL',QDL'])';
t2=Mat_G-Init_L'-t1; t2=Mat_G-Init_L'-t1;

6
OPF.m
View File

@ -225,10 +225,10 @@ wQD=1./(abs(mQD3P*sigma).^2);
%% 没有量测量 %% 没有量测量
% noLoadi=[1,5,6,10,11]; % noLoadi=[1,5,6,10,11];
% noLoadi=[1,8,9,10,12]; noLoadi=[1,8,10,12];
% noLoadi=Loadi; % noLoadi=Loadi;
% noLoadi=[1,11]; % noLoadi=[1,11];
noLoadi=[]; % noLoadi=[];
noPQi3P=zeros(length(noLoadi)*3,1); noPQi3P=zeros(length(noLoadi)*3,1);
noPQi3P(1:3:end)=(noLoadi-1)*3+1; noPQi3P(1:3:end)=(noLoadi-1)*3+1;
noPQi3P(2:3:end)=(noLoadi-1)*3+2; noPQi3P(2:3:end)=(noLoadi-1)*3+2;
@ -255,7 +255,7 @@ Precision=1e-5;
CenterA=0.1; CenterA=0.1;
%% 加误差 %% 加误差
Volt=1*Vp3m; Volt=1*Vp3m;
Volt=rVoltABCV; % Volt=mVoltABCV;
UAngel=Vp3a; UAngel=Vp3a;
while(abs(Gap)>Precision) while(abs(Gap)>Precision)
if KK>kmax if KK>kmax

27
feeder13/data1.txt
View File

@ -1,21 +1,18 @@
650 4.16 650 4.16
0 0
600 632 645 500 600 632 645 50
602 632 633 500 602 632 633 50
601 633 634 600 601 633 634 60
600 645 646 300 600 645 646 30
601 650 632 2000 601 650 632 200
602 684 652 800 602 684 652 80
601 632 671 2000 601 632 671 200
600 671 684 300 600 671 684 30
600 671 692 300 600 671 692 30
602 684 611 300 602 684 611 30
601 692 675 500 601 692 675 50
0 0
611 100 100 100 692 0 0 0
652 200 200 200
684 200 200 200
692 200 200 200
0 0
634 120 90 120 90 160 110 634 120 90 120 90 160 110
645 40 30 50 40 60 50 645 40 30 50 40 60 50