-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathgenPlotsForSolvedCase.m
253 lines (215 loc) · 6.15 KB
/
genPlotsForSolvedCase.m
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
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
function genPlotsForSolvedCase(casefile )
%UNTITLED Summary of this function goes here
% Detailed explanation goes here
%
%
%
%
v2struct(casefile);
if exist('figures')~=7
mkdir('figures');
end
cd('figures');
% plots
figx0=0;
figy0=1;
width=8;
height=5;
%
%
freqyMin=min(min(omegaVec))./(2*pi);
freqyMax=max(max(omegaVec))./(2*pi);
freqyOffSet=0;
figure1=figure('Units','inches',...
'Position',[figx0 figy0 width height],...
'PaperPositionMode','auto');
set(figure1, 'Name', 'genfreq');
plot(t,omegaVec./(2*pi),'lineWidth',2);
xlabel('Time (sec)', 'FontWeight','bold');
ylabel('$\frac{1}{2\pi}\boldmath{\omega}$ (Hz)');
axis([0 Tfinal freqyMin-freqyOffSet freqyMax+freqyOffSet]);
set(gca,'box','on');
set(gca,'fontSize',22);
set(0,'defaulttextinterpreter','latex')
grid on;
title(' Generator frequencies');
print -dpdf freq.pdf
print -depsc2 freq
%
%
%
%
%
%
angleyMin=min(min(deltaVec));
angleyMax=max(max(deltaVec));
angleyOffSet=0;
figure1=figure('Units','inches',...
'Position',[figx0 figy0 width height],...
'PaperPositionMode','auto');
set(figure1, 'Name', 'genangle');
plot(t,deltaVec-repmat(deltaS,1,length(t)),'lineWidth',2);
xlabel('Time (sec)', 'FontWeight','bold');
ylabel('$\boldmath{\delta}-\boldmath{\delta}^{s}$ (Rad)');
axis([0 Tfinal angleyMin-angleyOffSet angleyMax+angleyOffSet]);
set(gca,'box','on');
set(gca,'fontSize',22);
set(0,'defaulttextinterpreter','latex')
grid on;
title('Generator angles');
print -dpdf angles.pdf
print -depsc2 angles
eyMin=min(min(eVec));
eyMax=max(max(eVec));
eyOffSet=0;
figure1=figure('Units','inches',...
'Position',[figx0 figy0 width height],...
'PaperPositionMode','auto');
set(figure1, 'Name', 'gene');
plot(t,eVec-repmat(eS,1,length(t)),'lineWidth',2);
xlabel('Time (sec)', 'FontWeight','bold');
ylabel('$\boldmath{e}-\boldmath{e}^{s}$ (pu)');
axis([0 Tfinal eyMin-eyOffSet eyMax+eyOffSet]);
set(gca,'box','on');
set(gca,'fontSize',22);
set(0,'defaulttextinterpreter','latex')
grid on;
title('Generator EMF');
print -dpdf e.pdf
print -depsc2 e
myMin=min(min(mVec));
myMax=max(max(mVec));
myOffSet=0;
figure1=figure('Units','inches',...
'Position',[figx0 figy0 width height],...
'PaperPositionMode','auto');
set(figure1, 'Name', 'gene');
plot(t,mVec-repmat(mS,1,length(t)),'lineWidth',2);
xlabel('Time (sec)', 'FontWeight','bold');
ylabel('$\boldmath{m}-\boldmath{m}^{s}$ (pu)');
axis([0 Tfinal myMin-myOffSet myMax+myOffSet]);
set(gca,'box','on');
set(gca,'fontSize',22);
set(0,'defaulttextinterpreter','latex')
grid on;
title('Generator mechanical input');
print -dpdf m.pdf
print -depsc2 m
%
%
%
%
%
vyMin=min(min(vVec));
vyMax=max(max(vVec));
vyOffSet=0;
figure1=figure('Units','inches',...
'Position',[figx0 figy0 width height],...
'PaperPositionMode','auto');
set(figure1, 'Name', 'Voltage mags');
plot(t,vVec-repmat(vS,1,length(t)),'lineWidth',2);
xlabel('Time (sec)', 'FontWeight','bold');
ylabel('$\boldmath{v}-\boldmath{v}^{s}$ (pu)');
axis([0 Tfinal vyMin-vyOffSet vyMax+vyOffSet]);
set(gca,'box','on');
set(gca,'fontSize',22);
set(0,'defaulttextinterpreter','latex')
grid on;
title('Nodal voltage magnitudes');
print -dpdf VoltageMags.pdf
print -depsc2 VoltageMags
thetayMin=min(min(thetaVec));
thetayMax=max(max(thetaVec));
thetayOffSet=0;
figure1=figure('Units','inches',...
'Position',[figx0 figy0 width height],...
'PaperPositionMode','auto');
set(figure1, 'Name', 'voltage angles');
plot(t,thetaVec-repmat(thetaS,1,length(t)),'lineWidth',2);
xlabel('Time (sec)', 'FontWeight','bold');
ylabel('$\boldmath{\theta}-\boldmath{\theta}^{s}$ (Rad)');
axis([0 Tfinal thetayMin-thetayOffSet thetayMax+thetayOffSet]);
set(gca,'box','on');
set(gca,'fontSize',22);
set(0,'defaulttextinterpreter','latex')
grid on;
title('Nodal voltage angles');
print -dpdf VoltageAngles.pdf
print -depsc2 VoltageAngles
pgyMin=min(min(pgVec));
pgyMax=max(max(pgVec));
pgyOffSet=0;
figure1=figure('Units','inches',...
'Position',[figx0 figy0 width height],...
'PaperPositionMode','auto');
set(figure1, 'Name', 'voltage angles');
plot(t,pgVec-repmat(pgS,1,length(t)),'lineWidth',2);
xlabel('Time (sec)', 'FontWeight','bold');
ylabel('$\boldmath{p}_{g}-\boldmath{p}_{g}^{s}$ (pu)');
axis([0 Tfinal pgyMin-pgyOffSet pgyMax+pgyOffSet]);
set(gca,'box','on');
set(gca,'fontSize',22);
set(0,'defaulttextinterpreter','latex')
grid on;
title('Generators real electrical power output');
print -dpdf pg.pdf
print -depsc2 pg
qgyMin=min(min(qgVec));
qgyMax=max(max(qgVec));
qgyOffSet=0;
figure1=figure('Units','inches',...
'Position',[figx0 figy0 width height],...
'PaperPositionMode','auto');
set(figure1, 'Name', 'pg');
plot(t,qgVec-repmat(qgS,1,length(t)),'lineWidth',2);
xlabel('Time (sec)', 'FontWeight','bold');
ylabel('$\boldmath{q}_g- \boldmath{q}_{g}^{s}$ (pu)');
axis([0 Tfinal qgyMin-qgyOffSet qgyMax+qgyOffSet]);
set(gca,'box','on');
set(gca,'fontSize',22);
set(0,'defaulttextinterpreter','latex')
grid on;
title('Generators reactive power output');
print -dpdf qg.pdf
print -depsc2 qg
prefyMin=min(min(prefVec));
prefyMax=max(max(prefVec));
prefyOffSet=0;
figure1=figure('Units','inches',...
'Position',[figx0 figy0 width height],...
'PaperPositionMode','auto');
set(figure1, 'Name', 'pg');
plot(t,prefVec-repmat(prefS,1,length(t)),'lineWidth',2);
xlabel('Time (sec)', 'FontWeight','bold');
ylabel('$\boldmath{r}-\boldmath{r}^s$ (pu)');
axis([0 Tfinal prefyMin-prefyOffSet prefyMax+prefyOffSet]);
set(gca,'box','on');
set(gca,'fontSize',22);
set(0,'defaulttextinterpreter','latex')
grid on;
title('Generator governor signal');
print -dpdf pref.pdf
print -depsc2 pref
fyMin=min(min(fVec));
fyMax=max(max(fVec));
fyOffSet=0;
figure1=figure('Units','inches',...
'Position',[figx0 figy0 width height],...
'PaperPositionMode','auto');
set(figure1, 'Name', 'pg');
plot(t,fVec-repmat(fS,1,length(t)),'lineWidth',2);
xlabel('Time (sec)', 'FontWeight','bold');
ylabel('$\boldmath{f}-\boldmath{f}^{s}$ (pu)');
axis([0 Tfinal fyMin-fyOffSet fyMax+fyOffSet]);
set(gca,'box','on');
set(gca,'fontSize',22);
set(0,'defaulttextinterpreter','latex')
grid on;
title('Exciter field voltage');
print -dpdf f.pdf
print -depsc2 f
get(0, 'DefaultFigureVisible');
set(0, 'DefaultFigureVisible', 'on');
cd(CurrentDirectory);
end
end