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

1.去掉了不不要的Loadi 

2.加了123节点算例。
3.可以收敛。可以算123节点。

Signed-off-by: facat@lab.com <facat@lab.com>
This commit is contained in:
facat@lab.com 2014-11-03 21:35:41 +08:00
parent 51c1aaeaff
commit 5943513ea7
5 changed files with 299 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
FormG.m
View File

@ -1,7 +1,7 @@
function Mat_G=FormG(PD,QD,Loadi) function Mat_G=FormG(PD,QD,Loadi)
Mat_G=[ Mat_G=[
sparse(PD(Loadi)); sparse(PD);
sparse(QD(Loadi)); sparse(QD);
]; ];
end end

17
OPF.m
View File

@ -1,10 +1,19 @@
tic tic
clc clc
clear clear
lineZ=readLineZ('feeder13\lineParameter.txt'); lineZ=readLineZ('feeder123\lineParameter.txt');
[ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ... [ fsY0, fsY1, fsY2,phaseASpotLoadP,phaseBSpotLoadP,phaseCSpotLoadP ...
phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY ... phaseASpotLoadQ,phaseBSpotLoadQ,phaseCSpotLoadQ,setIJ,nodeNum,Balance,phaseABCY ...
cap]=dataRead(lineZ,'feeder13\data1.txt'); cap]=dataRead(lineZ,'feeder123\data.txt');
phaseASpotLoadP(phaseASpotLoadP==0)=0.002;
phaseBSpotLoadP(phaseBSpotLoadP==0)=0.002;
phaseCSpotLoadP(phaseCSpotLoadP==0)=0.002;
phaseASpotLoadQ(phaseASpotLoadQ==0)=0.002;
phaseBSpotLoadQ(phaseBSpotLoadQ==0)=0.002;
phaseCSpotLoadQ(phaseCSpotLoadQ==0)=0.002;
%% ³±Á÷¼ÆËãbegin %% ³±Á÷¼ÆËãbegin
a=exp(1j*2*pi/3); a=exp(1j*2*pi/3);
Tp2f=1/3*[1 1 1; Tp2f=1/3*[1 1 1;
@ -207,7 +216,7 @@ mPD3P=rPD3P.*PD3PSigma;
QD3PSigma=(1+normrnd(0,sigma,length(rQD3P),1)); QD3PSigma=(1+normrnd(0,sigma,length(rQD3P),1));
mQD3P=rPD3P.*QD3PSigma; mQD3P=rPD3P.*QD3PSigma;
%Á¿²â·½²î %Á¿²â·½²î
wVolt=1./(abs(mVoltABCV*sigma).^2)*0; wVolt=1./(abs(mVoltABCV*sigma).^2)*1;
wPD=1./(abs(mPD3P*sigma).^2)*1; wPD=1./(abs(mPD3P*sigma).^2)*1;
wQD=1./(abs(mQD3P*sigma).^2)*1; wQD=1./(abs(mQD3P*sigma).^2)*1;
% wVolt=1./abs(VoltSigma).^2; % wVolt=1./abs(VoltSigma).^2;
@ -223,7 +232,7 @@ Init_Y=sparse(1,2*Busnum*3);%
Gap=(Init_L*Init_Z'-Init_U*Init_W'); Gap=(Init_L*Init_Z'-Init_U*Init_W');
KK=0; KK=0;
ContrlCount=size(Loadi,1)*2+Busnum*6; ContrlCount=size(Loadi,1)*2+Busnum*6;
kmax=30; kmax=300;
%% %%
Precision=1e-5; Precision=1e-5;
CenterA=0.1; CenterA=0.1;

204
feeder123/data.txt Normal file
View File

@ -0,0 +1,204 @@
1 4.16
0
6 1 2 175
6 1 3 250
1 1 7 300
6 3 4 200
6 3 5 325
6 5 6 250
1 7 8 200
6 8 12 225
6 8 9 225
1 8 13 300
6 9 14 425
6 13 34 150
2 13 18 825
6 14 11 250
6 14 10 250
6 15 16 375
6 15 17 350
6 18 19 250
2 18 21 300
6 19 20 325
6 21 22 525
2 21 23 250
6 23 24 550
2 23 25 275
6 25 26 350
2 25 28 200
6 26 27 275
6 26 31 225
6 27 33 500
2 28 29 300
2 29 30 350
6 31 32 300
6 34 15 100
6 35 36 650
1 35 40 250
6 36 37 300
6 36 38 250
6 38 39 325
6 40 41 325
1 40 42 250
6 42 43 500
1 42 44 200
6 44 45 200
1 44 47 250
6 45 46 300
4 47 48 150
4 47 49 250
4 49 50 250
4 50 51 250
1 52 53 200
1 53 54 125
1 54 55 275
3 54 57 350
1 55 56 275
6 57 58 250
3 57 60 750
6 58 59 250
5 60 61 550
6 60 62 250
6 62 63 175
6 63 64 350
6 64 65 425
6 65 66 325
6 67 68 200
3 67 72 275
3 67 97 250
6 68 69 275
6 69 70 325
6 70 71 275
6 72 73 275
3 72 76 200
6 73 74 350
6 74 75 400
6 76 77 400
3 76 86 700
6 77 78 100
6 78 79 225
6 78 80 475
6 80 81 475
6 81 82 250
6 81 84 675
6 82 83 250
6 84 85 475
6 86 87 450
6 87 88 175
6 87 89 275
6 89 90 225
6 89 91 225
6 91 92 300
6 91 93 225
6 93 94 275
6 93 95 300
6 95 96 200
3 97 98 275
3 98 99 550
3 99 100 300
6 101 102 225
3 101 105 275
6 102 103 325
6 103 104 700
6 105 106 225
3 105 108 325
6 106 107 575
6 108 109 450
6 109 110 300
6 110 111 575
6 110 112 125
6 112 113 525
6 113 114 325
6 13 52 100
6 18 35 100
6 60 67 100
6 97 101 100
0
0
1 40 20 0 0 0 0
2 0 0 20 10 0 0
4 0 0 0 0 40 20
5 0 0 0 0 20 10
6 0 0 0 0 40 20
7 20 10 0 0 0 0
9 40 20 0 0 0 0
10 20 10 0 0 0 0
11 40 20 0 0 0 0
12 0 0 20 10 0 0
16 0 0 0 0 40 20
17 0 0 0 0 20 10
19 40 20 0 0 0 0
20 40 20 0 0 0 0
22 0 0 40 20 0 0
24 0 0 0 0 40 20
28 40 20 0 0 0 0
29 40 20 0 0 0 0
30 0 0 0 0 40 20
31 0 0 0 0 20 10
32 0 0 0 0 20 10
33 40 20 0 0 0 0
34 0 0 0 0 40 20
35 40 20 0 0 0 0
37 40 20 0 0 0 0
38 0 0 20 10 0 0
39 0 0 20 10 0 0
41 0 0 0 0 20 10
42 20 10 0 0 0 0
43 0 0 40 20 0 0
45 20 10 0 0 0 0
46 20 10 0 0 0 0
47 35 25 35 25 35 25
48 70 50 70 50 70 50
49 35 25 70 50 35 20
50 0 0 0 0 40 20
51 20 10 0 0 0 0
52 40 20 0 0 0 0
53 40 20 0 0 0 0
55 20 10 0 0 0 0
56 0 0 20 10 0 0
58 0 0 20 10 0 0
59 0 0 20 10 0 0
60 20 10 0 0 0 0
62 0 0 0 0 40 20
63 40 20 0 0 0 0
64 0 0 75 35 0 0
65 35 25 35 25 70 50
66 0 0 0 0 75 35
68 20 10 0 0 0 0
69 40 20 0 0 0 0
70 20 10 0 0 0 0
71 40 20 0 0 0 0
73 0 0 0 0 40 20
74 0 0 0 0 40 20
75 0 0 0 0 40 20
76 105 80 70 50 70 50
77 0 0 40 20 0 0
79 40 20 0 0 0 0
80 0 0 40 20 0 0
82 40 20 0 0 0 0
83 0 0 0 0 20 10
84 0 0 0 0 20 10
85 0 0 0 0 40 20
86 0 0 20 10 0 0
87 0 0 40 20 0 0
88 40 20 0 0 0 0
90 0 0 40 20 0 0
92 0 0 0 0 40 20
94 40 20 0 0 0 0
95 0 0 20 10 0 0
96 0 0 20 10 0 0
98 40 20 0 0 0 0
99 0 0 40 20 0 0
100 0 0 0 0 40 20
102 0 0 0 0 20 10
103 0 0 0 0 40 20
104 0 0 0 0 40 20
106 0 0 40 20 0 0
107 0 0 40 20 0 0
109 40 20 0 0 0 0
111 20 10 0 0 0 0
112 20 10 0 0 0 0
113 40 20 0 0 0 0
114 20 10 0 0 0 0
0
0

78
feeder123/lineParameter.txt Normal file
View File

@ -0,0 +1,78 @@
1
0.4576 0.1560 0.1535
0.1560 0.4666 0.1580
0.1535 0.1580 0.4615
0
1.0780 0.5017 0.3849
0.5017 1.0482 0.4236
0.3849 0.4236 1.0651
0
5.6765 -1.8319 -0.6982
-1.8319 5.9809 -1.1645
-0.6982 -1.1645 5.3971
0
2
0.4666 0.1580 0.1560
0.1580 0.4615 0.1535
0.1560 0.1535 0.4576
0
1.0482 0.4236 0.5017
0.4236 1.0651 0.3849
0.5017 0.3849 1.0780
0
5.9809 -1.1645 -1.8319
-1.1645 5.3971 -0.6982
-1.8319 -0.6982 5.6765
0
3
0.4615 0.1535 0.1580
0.1535 0.4576 0.1560
0.1580 0.1560 0.4666
0
1.0651 0.3849 0.4236
0.3849 1.0780 0.5017
0.4236 0.5017 1.0482
0
5.3971 -0.6982 -1.1645
-0.6982 5.6765 -1.8319
-1.1645 -1.8319 5.9809
0
4
0.4615 0.1580 0.1535
0.1580 0.4666 0.1560
0.1535 0.1560 0.4576
0
1.0651 0.4236 0.3849
0.4236 1.0482 0.5017
0.3849 0.5017 1.0780
0
5.3971 -1.1645 -0.6982
-1.1645 5.9809 -1.8319
-0.6982 -1.8319 5.6765
0
5
0.4666 0.1560 0.1580
0.1560 0.4576 0.1535
0.1580 0.1535 0.4615
0
1.0482 0.5017 0.4236
0.5017 1.0780 0.3849
0.4236 0.3849 1.0651
0
5.9809 -1.8319 -1.1645
-1.8319 5.6765 -0.6982
-1.1645 -0.6982 5.3971
0
6
0.4576 0.1535 0.1560
0.1535 0.4615 0.1580
0.1560 0.1580 0.4666
0
1.0780 0.3849 0.5017
0.3849 1.0651 0.4236
0.5017 0.4236 1.0482
0
5.6765 -0.6982 -1.8319
-0.6982 5.3971 -1.1645
-1.8319 -1.1645 5.9809
0

4
func_deltF.m
View File

@ -1,6 +1,6 @@
function deltF=func_deltF(wVolt,wPD,wQD,PD0,PD,QD,QD0,Volt,mVolt,Busnum,Loadi) function deltF=func_deltF(wVolt,wPD,wQD,PD0,PD,QD,QD0,Volt,mVolt,Busnum,Loadi)
t3=2*wPD.*(PD(Loadi)-PD0(Loadi)); t3=2*wPD.*(PD-PD0);
t4=2*wQD.*(QD(Loadi)-QD0(Loadi)); t4=2*wQD.*(QD-QD0);
t5=2*wVolt.*(Volt-mVolt); t5=2*wVolt.*(Volt-mVolt);
deltF=[sparse(t3); deltF=[sparse(t3);
sparse(t4); sparse(t4);