Lampiran 5

Data Bus Sistem Kelistrikan PT. PLN Wilayah SULSELTRABAR

function busdt=busdatas(num)

%
|Bus |
QLi |
Qmin |
busdata35 = [ 1
1.916
-90.09
2
1.861
0
3
1.439
0
4
1.936
-0.72
5
-0.079
0
6
0
-61.78
7
0.692
0
8
5.385
0
9
0.00
0
10
18.76
-94.24
11
19.271
0
12
6.49
0
13
21.9
0
14
2.639
0
15
1.7
0
16
6.7
0
17
6.8
0
18
7.06
0
19
6.4
0
20
5.296
-139.43
21
5.86
0
22
7.57
0
23
5.7
0

Type | Vsp
Qmax |
1
1.015
90.09;
2
1.033
0;
3
1.037
0;
3
0.999
0.72;
2
0.992
0;
3
0.992
61.78;
2
0.986
0;
2
0.965
0;
2
0.981
0;
3
0.967
94.24;
2
0.986
0;
3
0.951
0;
2
0.987
0;
2
0.973
0;
2
0.969
0;
2
0.992
0;
2
1.013
0;
2
1.000
0;
2
1.013
0;
3
1.009
139.43;
2
0.990
0;
2
0.96
0;
2
0.959
0;

theta

PGi

62.1

QGi

PLi

-4.8

4.316

5.425

4.765

0.6

5.552

-4.497

51.8

16.5

14.9

1.196

15.787

18.76

138

81.2

44.12

24.002

42.751

10.639

10.86

30.8

3.59

19.97

18.6

6.8

14.773

189.1

46

18.637

4.5

15.8

24
9.6

9.6

2.8

0.014

25.91

15.57
0.0
0.8

12.8

0.992

0.12

0.987

19.6

0.972

-0.013

0.989

32.866

0.994
14.3;
2
0.95
0;
3
0.994
0;
2
0.966
0;
2
0.994
28.22;
3
1.00
0;
3
1.00
0];

-22.667

12

0.0

40.951

4.8

7.16

2
2
2
0;

28

0.0

0;
27

-6.142

0;
26

0.0

0.96
0;

25

-1.32

29
-14.3
30
0
31
0
32
0
33
-28.22
34
0
35
0

2
0;
3

switch num
case 35
busdt = busdata35;
end

10.34

0.0
0.158

Lampiran 6
Data Saluran Sistem Kelistrikan SULSELTRABAR
function linedt = linedatas(num)

%
| From | To
%
| Bus | Bus
linedata35 =[1
2
1
4
2
3
2
5
4
5
5
7
5
8
6
5
7
8
8
9
8
10
8
25
1.001
9
10
10
11
10
12
10
13
1.001
10
30
10
34
11
32
1.001
12
10
13
14
15
13
16
17
16
5
17
18
18
17
19
16
19
21
20
19
21
23
21
22
22
23
23
24
24
15
25
26
25
27
25
28
26
25
27
28
28
29
30
28
30
29
31
30

|
|

R
|
X
|
pu
| pu
|
0.02627
0.09440
0.03076
0.11023
0.02630
0.09451
0.03663
0.13159
0.01388
0.04974
0.02314
0.08290
0.04732
0.16958
0.00393
0.01413
0.02419
0.08667
0.01090
0.03919
0.02382
0.08535
0.0
0.13164

B/2 | X'mer |
pu
| TAP (a) |
0.00743
1
0.01012
1
0.00744
1
0.01819
1
0.00670
1
0.00670
1
0.02283
1
0.00111
1
0.01167
1
0.00493
1
0.01149
1
0.0

0.01683
0.00363
0.00236
0.00192

0.06049
0.01300
0.00848
0.01318

0.00761
0.00175
0.00067
0.00249

1
1
1

0.0
0.0
0.0

0.20794
0.55350
0.41587

0.0
0.0
0.0

0.99
0.99

0.00236
0.00354
0.00485
0.03137
0.01002
0.01959
0.01959
0.02821
0.02289
0.01053
0.07195
0.04064
0.03131
0.02431
0.03333
0.01638
0.08188
0.18159
0.01638
0.13631
0.03420
0.05828
0.02408
0.06069

0.00848
0.02128
0.03324
0.18876
0.03599
0.07039
0.07039
0.10138
0.08153
0.06335
0.25851
0.14603
0.11249
0.08733
0.11974
0.03006
0.15032
0.33335
0.03006
0.25024
0.06278
0.10699
0.04421
0.11141

0.00067
0.00271
0.00627
0.02406
0.00283
0.00554
0.00554
0.00964
0.00804
0.00807
0.02035
0.01149
0.00885
0.00687
0.00942
0.00942
0.00045
0.00100
0.00942
0.00075
0.00019
0.00032
0.00013
0.00034

1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1

32

33

0.02023

0.03714

0.00011

31
35

25
30

0.1885
0.0

0.3292
0.55350

0.0
0.0

1.001
0.99];
switch num
case 35
linedt = linedata35;
end

Lampiran 7
Data Pengukuran
% Traditional Measurement Data..
% Vi - 1, Pi - 2, Qi - 3, Pij - 4, Qij - 5, Iij - 6;
% Institute of Electrical and Electronic Engineering (IEEE)
function zdt = zdatas(num)
%
zdata35

|Msnt |Type | Value | From | To | Rii |

= [ %---- Voltage Magnitude ------------%
1
1
1.015
1
0
1e-4;
2
1
1.033
2
0
1e-4;
2
1
0.999
4
0
1e-4;
3
1
0.986
7
0
1e-4;
4
1
0.967
10
0
1e-4;
5
1
0.986
13
0
1e-4;
6
1
0.976
27
0
1e-4;
7
1
0.994
29
0
1e-4;
8
1
0.964
33
0
1e-4;
%-----------------------------------%
%---- Real Power Injection ---------%
3
2
0.00
5
0
1e-4;
4
2
0.00
6
0
1e-4;
5
2
0.00
8
0
1e-4;
6
2
0.00
10
0
1e-4;
7
2
0.00
11
0
1e-4;
8
2
0.00
13
0
1e-4;
9
2
0.00
14
0
1e-4;
10
2
0.00
15
0
1e-4;
11
2
0.00
16
0
1e-4;
12
2
0.00
18
0
1e-4;
13
2
0.00
20
0
1e-4;
14
2
0.00
21
0
1e-4;
15
2
0.00
24
0
1e-4;
16
2
0.00
25
0
1e-4;
17
2
-0.0350 26
0
1e-4;
18
2
0.00
28
0
1e-4;
19
2
-0.0240 29
0
1e-4;
%------------------------------------%
%---- Reative Power Injection -------%
20
3
0.00
4
0
1e-4;
21
3
0.1538
5
0
1e-4;
22
3
0.00
6
0
1e-4;
23
3
0.2096
8
0
1e-4;
24
3
-0.0200 10
0
1e-4;
25
3
0.2066 11
0
1e-4;
26
3
0.1483 13
0
1e-4;
27
3
-0.0160 14
0
1e-4;
28
3
-0.0250 15
0
1e-4;
29
3
-0.0180 16
0
1e-4;
30
3
-0.0090 18
0
1e-4;
31
3
-0.0070 20
0
1e-4;
32
3
-0.1120 21
0
1e-4;
33
3
-0.0670 24
0
1e-4;

34
3
0.4214 25
0
1e-4;
35
3
-0.0230 26
0
1e-4;
36
3
0.00
28
0
1e-4;
37
3
0.2857 29
0
1e-4;
36
3
0.1457 30
0
1e-4;
37
3
0.3057 32
0
1e-4;
%------------------------------------%
%------ Real Power Flow ------------- %
38
4
0.0767
2
3
64e-6;
39
4
0.3392
2
5
64e-6;
42
4
0.5115
4
5
64e-6;
43
4
0.2868
5
7
64e-6;
44
4
0.4148
5
8
64e-6;
45
4
0.8233
7
8
64e-6;
46
4
0.1660
6
5
64e-6;
47
4
0.6611
8
9
64e-6;
48
4
0.3113
8
10
64e-6;
48
4
0.0953
8
25
64e-6;
48
4
0.4473
9
10
64e-6;
49
4
0.1733 10
11
64e-6;
50
4
0.1198 10
12
64e-6;
51
4
-0.0552 10
13
64e-6;
51
4
0.5167 10
30
64e-6;
51
4
-0.2933 10
34
64e-6;
50
4
0.3689 11
32
64e-6;
50
4
0.0397 13
14
64e-6;
52
4
0.0707 15
13
64e-6;
53
4
0.1123 16
17
64e-6;
54
4
0.0795 16
5
64e-6;
54
4
0.0308 17
18
64e-6;
55
4
0.2395 19
16
64e-6;
55
4
0.3140 19
21
64e-6;
56
4
0.4272 20
19
64e-6;
58
4
0.4868 21
23
64e-6;
59
4
0.5436 21
22
64e-6;
60
4
0.405
22
23
64e-6;
60
4
0.1969 23
24
64e-6;
61
4
0.2139 24
15
64e-6;
62
4
0.1151 25
26
64e-6;
63
4
0.0771 25
27
64e-6;
63
4
0.0212 25
28
64e-6;
63
4
0.0464 27
28
64e-6;
64
4
-0.5778 28
29
64e-6;
65
4
0.4810 30
28
64e-6;
65
4
1.3246 30
29
64e-6;
66
4
0.3035 31
30
64e-6;
67
4
0.3678 32
33
64e-6;
69
4
0.002
35
30
64e-6;
%------------------------------------%
%------ Real Power Flow ------------- %
70
5
-0.0037
2
3
64e-6;
71
5
-0.0291
2
5
64e-6;
74
5
-0.057
4
5
64e-6;
75
5
0.0844
5
7
64e-6;
76
5
0.0019
5
8
64e-6;
77
5
0.0680
7
8
64e-6;

78
5
0.0304
6
5
64e-6;
79
5
0.0265
8
9
64e-6;
80
5
0.0097
8
10
64e-6;
80
5
-0.0655
8
25
64e-6;
80
5
0.0053
9
10
64e-6;
81
5
0.0328 10
11
64e-6;
82
5
0.0461 10
12
64e-6;
83
5
0.0387 10
13
64e-6;
83
5
-1.7446 10
30
64e-6;
83
5
-0.3467 10
34
64e-6;
82
5
0.0816 11
32
64e-6;
83
5
0.0182 13
14
64e-6;
84
5
0.0031 15
13
64e-6;
85
5
0.0057 16
17
64e-6;
86
5
-0.0133 16
5
64e-6;
86
5
0.0064 17
18
64e-6;
87
5
0.0261 19
16
64e-6;
87
5
0.0494 19
21
64e-6;
88
5
0.0916 20
19
64e-6;
89
5
0.0359 21
23
64e-6;
90
5
0.0586 21
22
64e-6;
91
5
0.0016 22
23
64e-6;
92
5
0.0021 23
24
64e-6;
93
5
-0.0116 24
15
64e-6;
94
5
0.0895 25
26
64e-6;
95
5
1.0896 25
27
64e-6;
95
5
0.0523 25
28
64e-6;
95
5
0.1157 27
28
64e-6;
96
5
-0.1125 28
29
64e-6;
97
5
-0.0193 30
28
64e-6;
97
5
-0.3225 30
29
64e-6;
98
5
-0.2632 31
30
64e-6;
99
5
0.0496 32
33
64e-6;
99
5
-0.2625 31
25
64e-6;
101
5
-0.004
35
30
64e-6];
%--------------------------------------%
switch num
case 35
zdt = zdata35;
end
end

Lampiran 8
% Line Data for B-Bus (Shunt Admittance)Formation.
% Institute of Electrical and Electronic Engineering (IEEE)
function bbus = bbusppg(num)
linedata = linedatas(num);
fb = linedata(:,1);
tb = linedata(:,2);
b = linedata(:,5);
nbus = max(max(fb),max(tb));
nbranch = length(fb);
bbus = zeros(nbus,nbus);

% Returns B-bus..

% no. of buses...
% no. of branches...

for k=1:nbranch
bbus(fb(k),tb(k)) = b(k);
bbus(tb(k),fb(k)) = bbus(fb(k),tb(k));
end

Lampiran 9
% Program to form Admittance And Impedance Bus Formation....
% with Transformer Tap setting..
% Praviraj P G, MTech I Year, EE Dept., IIT Roorkee, India, Email
:[email protected]
% Institute of Electrical and Electronic Engineering (IEEE)
function ybus = ybusppg(num)

% Returns ybus

linedata = linedatas(num); % Calling "linedata6.m" for Line Data...

fb = linedata(:,1);
% From bus number...
tb = linedata(:,2);
% To bus number...
r = linedata(:,3);
% Resistance, R...
x = linedata(:,4);
% Reactance, X...
b = linedata(:,5);
% Ground Admittance, B/2...
a = linedata(:,6);
% Tap setting value..
z = r + i*x;
% Z matrix...
y = 1./z;
% To get inverse of each element...
b = i*b;
% Make B imaginary...
nbus = max(max(fb),max(tb));
nbranch = length(fb);
ybus = zeros(nbus,nbus);

% no. of buses...
% no. of branches...
% Initialise YBus...

% Formation of the Off Diagonal Elements...

for k=1:nbranch
ybus(fb(k),tb(k)) = ybus(fb(k),tb(k))-y(k)/a(k);
ybus(tb(k),fb(k)) = ybus(fb(k),tb(k));
end
% Formation of Diagonal
for m =1:nbus
for n =1:nbranch
if fb(n) == m
ybus(m,m) =
elseif tb(n) ==
ybus(m,m) =
end
end
end
%ybus;
%zbus = inv(ybus);

Elements....

ybus(m,m) + y(n)/(a(n)^2) + b(n);

m
ybus(m,m) + y(n) + b(n);

% Bus Admittance Matrix

% Bus Impedance Matrix

Lampiran 10
% Polar to Rectangular Conversion
% Praviraj P G, MTech I Year, EE Dept., IIT Roorkee, India, Email
:[email protected]
% [RECT] = RECT2POL(RHO, THETA)
% RECT - Complex matrix or number, RECT = A + jB, A = Real, B =
Imaginary
% RHO - Magnitude
% THETA - Angle in radians
% Institute of Electrical and Electronic Engineering (IEEE)
function rect = pol2rect(rho,theta)
rect = RHO.*cos(theta) + j*RHO.*sin(theta);

Lampiran 11
% Polar to Rectangular Conversion
% Praviraj P G, MTech I Year, EE Dept., IIT Roorkee, India, Email
:[email protected]
% [RHO THETA] = RECT2POL(X)
% X - Complex matrix or number, X = A + jB
% RHO - Magnitude
% THETA - Angle in radians
% Institute of Electrical and Electronic Engineering (IEEE)
function [rho theta] = rect2pol(x)
rho = sqrt(real(x).^2 + imag(x).^2);
theta = atan(imag(x)./real(x));

Lampiran 12 1
Program State Estimasi dengan menggunakan Weighted Least Square
% Power System State Estimation using Weighted Least Square Method..
% Praviraj P G, MTech I Year, EE Dept., IIT Roorkee, India, Email
:[email protected]
% Institute of Electrical and Electronic Engineering (IEEE)
num = 35;
ybus = ybusppg(num); % Get YBus..
zdata = zdatas(num); % Get Measurement data..
bpq = bbusppg(num); % Get B data..
nbus = max(max(zdata(:,4)),max(zdata(:,5))); % Get number of buses..
type = zdata(:,2); % Type of measurement, Vi - 1, Pi - 2, Qi - 3,
Pij - 4, Qij - 5, Iij - 6..
z = zdata(:,3); % Measuement values..
fbus = zdata(:,4); % From bus..
tbus = zdata(:,5); % To bus..
Ri = diag(zdata(:,6)); % Measurement Error..
V = ones(nbus,1); % Initialize the bus voltages..
del = zeros(nbus,1); % Initialize the bus angles..
E = [del(2:end); V];
% State Vector..
G = real(ybus);
B = imag(ybus);
vi = find(type == 1); % Index of voltage magnitude measurements..
ppi = find(type == 2); % Index of real power injection
measurements..
qi = find(type == 3); % Index of reactive power injection
measurements..
pf = find(type == 4); % Index of real powerflow measurements..
qf = find(type == 5); % Index of reactive powerflow measurements..
nvi = length(vi); % Number of Voltage measurements..
npi = length(ppi); % Number of Real Power Injection measurements..
nqi = length(qi); % Number of Reactive Power Injection
measurements..
npf = length(pf); % Number of Real Power Flow measurements..
nqf = length(qf); % Number of Reactive Power Flow measurements..
iter = 1;
tol = 5;
while(tol > 1e-3)
%Measurement Function, h
h1 = V(fbus(vi),1);
h2 = zeros(npi,1);
h3 = zeros(nqi,1);
h4 = zeros(npf,1);
h5 = zeros(nqf,1);
for i = 1:npi
m = fbus(ppi(i));
for k = 1:nbus

h2(i) = h2(i) + V(m)*V(k)*(G(m,k)*cos(del(m)-del(k)) +

B(m,k)*sin(del(m)-del(k)));
end
end
for i = 1:nqi
m = fbus(qi(i));
for k = 1:nbus
h3(i) = h3(i) + V(m)*V(k)*(G(m,k)*sin(del(m)-del(k)) B(m,k)*cos(del(m)-del(k)));
end
end
for i = 1:npf
m = fbus(pf(i));
n = tbus(pf(i));
h4(i) = -V(m)^2*G(m,n) - V(m)*V(n)*(-G(m,n)*cos(del(m)del(n)) - B(m,n)*sin(del(m)-del(n)));
end
for i = 1:nqf
m = fbus(qf(i));
n = tbus(qf(i));
h5(i) = -V(m)^2*(-B(m,n)+bpq(m,n)) - V(m)*V(n)*(G(m,n)*sin(del(m)-del(n)) + B(m,n)*cos(del(m)-del(n)));
end
h = [h1; h2; h3; h4; h5];
% Residue..
r = z - h;
% Jacobian..
% H11 - Derivative of V with respect to angles.. All Zeros
H11 = zeros(nvi,nbus-1);
% H12 - Derivative of V with respect to V..
H12 = zeros(nvi,nbus);
for k = 1:nvi
for n = 1:nbus
if n == k
H12(k,n) = 1;
end
end
end
% H21 - Derivative of Real Power Injections with Angles..
H21 = zeros(npi,nbus-1);
for i = 1:npi
m = fbus(ppi(i));
for k = 1:(nbus-1)
if k+1 == m
for n = 1:nbus
H21(i,k) = H21(i,k) + V(m)* V(n)*(G(m,n)*sin(del(m)-del(n)) + B(m,n)*cos(del(m)-del(n)));

Lampiran 12 - 2
end
H21(i,k) = H21(i,k) - V(m)^2*B(m,m);
else
H21(i,k) = V(m)* V(k+1)*(G(m,k+1)*sin(del(m)del(k+1)) - B(m,k+1)*cos(del(m)-del(k+1)));
end
end
end
% H22 - Derivative of Real Power Injections with V..
H22 = zeros(npi,nbus);
for i = 1:npi
m = fbus(ppi(i));
for k = 1:(nbus)
if k == m
for n = 1:nbus
H22(i,k) = H22(i,k) + V(n)*(G(m,n)*cos(del(m)del(n)) + B(m,n)*sin(del(m)-del(n)));
end
H22(i,k) = H22(i,k) + V(m)*G(m,m);
else
H22(i,k) = V(m)*(G(m,k)*cos(del(m)-del(k)) +
B(m,k)*sin(del(m)-del(k)));
end
end
end
% H31 - Derivative of Reactive Power Injections with Angles..
H31 = zeros(nqi,nbus-1);
for i = 1:nqi
m = fbus(qi(i));
for k = 1:(nbus-1)
if k+1 == m
for n = 1:nbus
H31(i,k) = H31(i,k) + V(m)*
V(n)*(G(m,n)*cos(del(m)-del(n)) + B(m,n)*sin(del(m)-del(n)));
end
H31(i,k) = H31(i,k) - V(m)^2*G(m,m);
else
H31(i,k) = V(m)* V(k+1)*(-G(m,k+1)*cos(del(m)del(k+1)) - B(m,k+1)*sin(del(m)-del(k+1)));
end
end
end
% H32 - Derivative of Reactive Power Injections with V..
H32 = zeros(nqi,nbus);
for i = 1:nqi
m = fbus(qi(i));
for k = 1:(nbus)
if k == m
for n = 1:nbus
H32(i,k) = H32(i,k) + V(n)*(G(m,n)*sin(del(m)del(n)) - B(m,n)*cos(del(m)-del(n)));

Lampiran 12 3
end
H32(i,k) = H32(i,k) - V(m)*B(m,m);
else
H32(i,k) = V(m)*(G(m,k)*sin(del(m)-del(k)) B(m,k)*cos(del(m)-del(k)));
end
end
end
% H41 - Derivative of Real Power Flows with Angles..
H41 = zeros(npf,nbus-1);
for i = 1:npf
m = fbus(pf(i));
n = tbus(pf(i));
for k = 1:(nbus-1)
if k+1 == m
H41(i,k) = V(m)* V(n)*(-G(m,n)*sin(del(m)-del(n)) +
B(m,n)*cos(del(m)-del(n)));
else if k+1 == n
H41(i,k) = -V(m)* V(n)*(-G(m,n)*sin(del(m)-del(n)) +
B(m,n)*cos(del(m)-del(n)));
else
H41(i,k) = 0;
end
end
end
end
% H42 - Derivative of Real Power Flows with V..
H42 = zeros(npf,nbus);
for i = 1:npf
m = fbus(pf(i));
n = tbus(pf(i));
for k = 1:nbus
if k == m
H42(i,k) = -V(n)*(-G(m,n)*cos(del(m)-del(n)) B(m,n)*sin(del(m)-del(n))) - 2*G(m,n)*V(m);
else if k == n
H42(i,k) = -V(m)*(-G(m,n)*cos(del(m)-del(n)) B(m,n)*sin(del(m)-del(n)));
else
H42(i,k) = 0;
end
end
end
end
% H51 - Derivative of Reactive Power Flows with Angles..
H51 = zeros(nqf,nbus-1);
for i = 1:nqf
m = fbus(qf(i));
n = tbus(qf(i));
for k = 1:(nbus-1)
if k+1 == m

Lampiran 12 4
H51(i,k) = -V(m)* V(n)*(-G(m,n)*cos(del(m)-del(n)) B(m,n)*sin(del(m)-del(n)));
else if k+1 == n
H51(i,k) = V(m)* V(n)*(-G(m,n)*cos(del(m)-del(n)) B(m,n)*sin(del(m)-del(n)));
else
H51(i,k) = 0;
end
end
end
end
% H52 - Derivative of Reactive Power Flows with V..
H52 = zeros(nqf,nbus);
for i = 1:nqf
m = fbus(qf(i));
n = tbus(qf(i));
for k = 1:nbus
if k == m
H52(i,k) = -V(n)*(-G(m,n)*sin(del(m)-del(n)) +
B(m,n)*cos(del(m)-del(n))) - 2*V(m)*(-B(m,n)+ bpq(m,n));
else if k == n
H52(i,k) = -V(m)*(-G(m,n)*sin(del(m)-del(n)) +
B(m,n)*cos(del(m)-del(n)));
else
H52(i,k) = 0;
end
end
end
end
% Measurement Jacobian, H..
H = [H11 H12; H21 H22; H31 H32; H41 H42; H51 H52];
% Gain Matrix, Gm..
Gm = H'*inv(Ri)*H;
%Objective Function..
J = sum(inv(Ri)*r.^2);
% State Vector..
dE = inv(Gm)*(H'*inv(Ri)*r);
E = E + dE;
del(2:end) = E(1:nbus-1);
V = E(nbus:end);
iter = iter + 1;
tol = max(abs(dE));
end
CvE = diag(inv(H'*inv(Ri)*H)); % Covariance matrix..
Del = 180/pi*del;
E2 = [V Del]; % Bus Voltages and angles..

Lampiran 12 - 5
disp('-------- State Estimation ------------------');
disp('--------------------------');
disp('| Bus |
V
| Angle | ');
disp('| No |
pu
| Degree | ');
disp('--------------------------');
for m = 1:n
fprintf('%4g', m); fprintf(' %8.4f', V(m)); fprintf('
Del(m)); fprintf('\n');
end
disp('---------------------------------------------');

%8.4f',