-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathdata33.m
More file actions
99 lines (88 loc) · 6.62 KB
/
data33.m
File metadata and controls
99 lines (88 loc) · 6.62 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
function [BD,LD,TL]=data33()
Vbase=12.66; %%---Base Voltage in kV---%%
Sbase=1; %%---Base Power in MVA---%%
Zbase=Vbase*Vbase/Sbase; %%---Base impedance in ohms---%%
%% ----Line Data (or Branch data) of the system----- %%
%%--- Line From To R X Line Charging Tap Angle----%%
%%---Number Bus Bus (in ohms) (in ohms) (in ohms) Ratio Shift----%%
LD = [ 1 1 2 0.0922 0.0470 0 1.0 0
2 2 3 0.4930 0.2511 0 1.0 0
3 3 4 0.3660 0.1864 0 1.0 0
4 4 5 0.3811 0.1941 0 1.0 0
5 5 6 0.8190 0.7070 0 1.0 0
6 6 7 0.1872 0.6188 0 1.0 0
7 7 8 0.7114 0.2351 0 1.0 0
8 8 9 1.0300 0.7400 0 1.0 0
9 9 10 1.0440 0.7400 0 1.0 0
10 10 11 0.1966 0.0650 0 1.0 0
11 11 12 0.3744 0.1238 0 1.0 0
12 12 13 1.4680 1.1550 0 1.0 0
13 13 14 0.5416 0.7129 0 1.0 0
14 14 15 0.5910 0.5260 0 1.0 0
15 15 16 0.7463 0.5450 0 1.0 0
16 16 17 1.2890 1.7210 0 1.0 0
17 17 18 0.7320 0.5740 0 1.0 0
18 2 19 0.1640 0.1565 0 1.0 0
19 19 20 1.5042 1.3554 0 1.0 0
20 20 21 0.4095 0.4784 0 1.0 0
21 21 22 0.7089 0.9373 0 1.0 0
22 3 23 0.4512 0.3083 0 1.0 0
23 23 24 0.8980 0.7091 0 1.0 0
24 24 25 0.8960 0.7011 0 1.0 0
25 6 26 0.2030 0.1034 0 1.0 0
26 26 27 0.2842 0.1447 0 1.0 0
27 27 28 1.0590 0.9337 0 1.0 0
28 28 29 0.8042 0.7006 0 1.0 0
29 29 30 0.5075 0.2585 0 1.0 0
30 30 31 0.9744 0.9630 0 1.0 0
31 31 32 0.3105 0.3619 0 1.0 0
32 32 33 0.3410 0.5302 0 1.0 0];
%% ----Bus Data of the system---- %%
%%--- Type 1 - PQ Bus ; Type 2 - PV Bus; Type 3 - Swing Bus
%%----Bus Bus V theta Pg Qg Pl Ql Qgmax Qgmin----%%
%%---Number Type (in pu) (in rad) (in kW)(in kVAR)(in kW)(in kVAR)(in pu) (in pu)----%%
BD = [ 1 3 1.0 0 0 0 0 0 0 0 2
2 1 1.0 0 0 0 100 60 0 0 1
3 1 1.0 0 0 0 90 40 0 0 2
4 1 1.0 0 0 0 120 80 0 0 3
5 1 1.0 0 0 0 60 30 0 0 1
6 1 1.0 0 0 0 60 20 0 0 2
7 1 1.0 0 0 0 200 100 0 0 3
8 1 1.0 0 0 0 200 100 0 0 3
9 1 1.0 0 0 0 60 20 0 0 2
10 1 1.0 0 0 0 60 20 0 0 3
11 1 1.0 0 0 0 45 30 0 0 3
12 1 1.0 0 0 0 60 35 0 0 1
13 1 1.0 0 0 0 60 35 0 0 3
14 1 1.0 0 0 0 120 80 0 0 1
15 1 1.0 0 0 0 60 10 0 0 3
16 1 1.0 0 0 0 60 20 0 0 2
17 1 1.0 0 0 0 60 20 0 0 3
18 1 1.0 0 0 0 90 40 0 0 2
19 1 1.0 0 0 0 90 40 0 0 1
20 1 1.0 0 0 0 90 40 0 0 3
21 1 1.0 0 0 0 90 40 0 0 2
22 1 1.0 0 0 0 90 40 0 0 1
23 1 1.0 0 0 0 90 50 0 0 3
24 1 1.0 0 0 0 420 200 0 0 3
25 1 1.0 0 0 0 420 200 0 0 3
26 1 1.0 0 0 0 60 25 0 0 3
27 1 1.0 0 0 0 60 25 0 0 2
28 1 1.0 0 0 0 60 20 0 0 3
29 1 1.0 0 0 0 120 70 0 0 3
30 1 1.0 0 0 0 200 600 0 0 3
31 1 1.0 0 0 0 150 70 0 0 1
32 1 1.0 0 0 0 210 100 0 0 1
33 1 1.0 0 0 0 60 40 0 0 3];
%% ----Tie Line Data of the system---- %%
%%--- Line From To R X
%%---Number Bus Bus (in ohms) (in ohms)
TL = [ 33 8 21 2.0000 2.0000 0 1.0 0
34 9 15 2.0000 2.0000 0 1.0 0
35 12 22 2.0000 2.0000 0 1.0 0
36 18 33 0.5000 0.5000 0 1.0 0
37 25 29 0.5000 0.5000 0 1.0 0];
Sl=[4.6 4.1 2.9 2.9 2.9 1.5 1.05 1.05 1.05 1.05 1.05 0.5 0.45 0.3 0.25 0.25 0.1 0.5 0.5 0.21 0.11 1.05 1.05 0.5 1.5 1.5 1.5 1.5 1.5 0.5 0.5 0.1 0.5 0.5 0.5 0.5 0.1]./Sbase; %% Line MVA Limit
BD(:,5:8)=BD(:,5:8)/(1000*Sbase); %%---Conversion of P and Q in pu quantity---%%
LD(:,4:5)=LD(:,4:5)/Zbase; %%----Conversion of R and X in pu quantity---%%
TL(:,4:5)=TL(:,4:5)/Zbase; %%----Conversion of R and X in pu quantity---%%