-
Notifications
You must be signed in to change notification settings - Fork 0
/
resultados.m
379 lines (322 loc) · 18.3 KB
/
resultados.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
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
function varargout = resultados(varargin)
% RESULTADOS M-file for resultados.fig
% RESULTADOS, by itself, creates a new RESULTADOS or raises the existing
% singleton*.
%
% H = RESULTADOS returns the handle to a new RESULTADOS or the handle to
% the existing singleton*.
%
% RESULTADOS('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in RESULTADOS.M with the given input arguments.
%
% RESULTADOS('Property','Value',...) creates a new RESULTADOS or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before resultados_OpeningFunction gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to resultados_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help resultados
% Last Modified by GUIDE v2.5 24-Mar-2010 22:02:28
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @resultados_OpeningFcn, ...
'gui_OutputFcn', @resultados_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
% --- Executes just before resultados is made visible.
function resultados_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to resultados (see VARARGIN)
% Choose default command line output for resultados
global matresultados matfilteredP matfilteredQ matP barra npodrioP npodrioQ datosporhoraP matdesvP mes
handles.output = hObject;
% Muestra la cantidad de datos de P y Q con errores de cálculo y los
% valores importados y finales por hora
set(handles.text3,'String',int2str(npodrioP));
set(handles.text4,'String',int2str(npodrioQ));
% Muestra los datos por hora importados y finales de P (inicialmente)
set(handles.uitable3,'ColumnWidth',{70 60});
set(handles.uitable3,'ColumnName',{'Importados', 'Finales'});
set(handles.uitable3,'Data',datosporhoraP);
% Cambia el nombre de la ventana
set(handles.figure1,'Name',strcat(['Resultados para la Subestación: ',barra]));
% Llena la tabla con los promedios de P y Q para las 25 horas
set(handles.uitable2,'ColumnWidth',{45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45});
set(handles.uitable2,'Data',matresultados);
% Gráfica la curva de carga promedio insertando el valor de la hora 25
% (pico) entre las horas 19 y 29 o 18 y 19 en dependencia de ser verano o
% invierno.
axes(handles.axes2);
horas = [1:1:25];
activa = [];
reactiva = [];
if ismember(mes,[3 4 5 6 7 8 9 10])
% Es Verano, el pico va entre las horas 19 y 20
activa = cat(2,matfilteredP(1:19),matfilteredP(25));
activa = cat(2,activa,matfilteredP(20:24));
reactiva = cat(2,matfilteredQ(1:19),matfilteredQ(25));
reactiva = cat(2,reactiva,matfilteredQ(20:24));
plot(horas, activa,'r.-',horas,reactiva,'b.-');xlabel('Hora');legend('P(MW)','Location','NorthWest','Q(Mvar)','Location','NorthWest');grid;
set(gca,'XTick',[0:1:25]);
set(gca,'XTickLabel',{'0';'';'';'';'';'5';'';'';'';'';'10';'';'';'';'';'15';'';'';'';'';'';'20';'';'';'';'24'});
%zoom on
else
% Es Invierno, el pico va entre las horas 18 y 19
activa = cat(2,matfilteredP(1:18),matfilteredP(25));
activa = cat(2,activa,matfilteredP(19:24));
reactiva = cat(2,matfilteredQ(1:18),matfilteredQ(25));
reactiva = cat(2,reactiva,matfilteredQ(19:24));
plot(horas, activa,'r.-',horas,reactiva,'b.-');xlabel('Hora');legend('P(MW)','Location','NorthWest','Q(Mvar)','Location','NorthWest');grid;
set(gca,'XTick',[0:1:25]);
set(gca,'XTickLabel',{'0';'';'';'';'';'5';'';'';'';'';'10';'';'';'';'';'15';'';'';'';'';'';'20';'';'';'';'24'});
%zoom on
end
% Gráfica la dispersión de P para la hora 1 (por defecto al crearse la
% ventana) y muestra la desviación estandar para el evento anterior.
axes(handles.axes3);
[valores horas]=size(matP);
clear horas;
axes(handles.axes3);
xrange=[1:1:valores];
fill([xrange(1) xrange(end) xrange(end) xrange(1)],[matfilteredP(1,1)-matdesvP(1) matfilteredP(1,1)-matdesvP(1) matfilteredP(1,1)+matdesvP(1) matfilteredP(1,1)+matdesvP(1)],[0.9 0.9 1]);hold on;
scatter(xrange,matP(:,1),'filled','bs');xlabel('Muestra No.');ylabel('MW','fontsize',7);grid;set(gca,'Layer','top');hold on;
set(handles.axes3,'xlim',[1 valores]);
line([xrange(1) xrange(end)],[matfilteredP(1,1) matfilteredP(1,1)],'LineWidth',2,'LineStyle','-','Color',[1 0 0]);hold off;
clear xrange
set(handles.text5,'String',num2str(matdesvP(1)));
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes resultados wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% --- Outputs from this function are returned to the command line.
function varargout = resultados_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --- Executes on selection change in popupmenu1.
function popupmenu1_Callback(hObject, eventdata, handles)
% hObject handle to popupmenu1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global matP matQ matdesvP matdesvQ matfilteredP matfilteredQ
horamenu = get(hObject,'Value');
potgrupo = get(get(handles.uipanel8,'SelectedObject'),'Tag');
% Actualiza el gráfico de la dispersión de datos
switch potgrupo
case 'radiobutton3'
[valores horas]=size(matP);
clear horas;
set(handles.text5,'String',num2str(matdesvP(horamenu)));
axes(handles.axes3);
xrange=[1:1:valores];
fill([xrange(1) xrange(end) xrange(end) xrange(1)],[matfilteredP(1,horamenu)-matdesvP(horamenu) matfilteredP(1,horamenu)-matdesvP(horamenu) matfilteredP(1,horamenu)+matdesvP(horamenu) matfilteredP(1,horamenu)+matdesvP(horamenu)],[0.9 0.9 1]);hold on;
scatter(xrange,matP(:,horamenu),'filled','bs');xlabel('Muestra No.');ylabel('MW','fontsize',7);grid;set(gca,'Layer','top');hold on;
set(handles.axes3,'xlim',[1 valores]);
line([xrange(1) xrange(end)],[matfilteredP(1,horamenu) matfilteredP(1,horamenu)],'LineWidth',2,'LineStyle','-','Color',[1 0 0]);hold off;
case 'radiobutton5'
[valores horas]=size(matQ);
clear horas;
set(handles.text5,'String',num2str(matdesvQ(horamenu)));
axes(handles.axes3);
xrange=[1:1:valores];
fill([xrange(1) xrange(end) xrange(end) xrange(1)],[matfilteredQ(1,horamenu)-matdesvQ(horamenu) matfilteredQ(1,horamenu)-matdesvQ(horamenu) matfilteredQ(1,horamenu)+matdesvQ(horamenu) matfilteredQ(1,horamenu)+matdesvQ(horamenu)],[0.9 0.9 1]);hold on;
scatter(xrange,matQ(:,horamenu),'filled','bs');xlabel('Muestra No.');ylabel('Mvar','fontsize',7);grid;set(gca,'Layer','top');hold on;
set(handles.axes3,'xlim',[1 valores]);
line([xrange(1) xrange(end)],[matfilteredQ(1,horamenu) matfilteredQ(1,horamenu)],'LineWidth',2,'LineStyle','-','Color',[1 0 0]);hold off;
end
clear xrange
% --- Executes during object creation, after setting all properties.
function popupmenu1_CreateFcn(hObject, eventdata, handles)
% hObject handle to popupmenu1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: popupmenu controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
% --- Executes when selected object is changed in uipanel8.
function uipanel8_SelectionChangeFcn(hObject, eventdata, handles)
% hObject handle to the selected object in uipanel8
% eventdata structure with the following fields (see UIBUTTONGROUP)
% EventName: string 'SelectionChanged' (read only)
% OldValue: handle of the previously selected object or empty if none was selected
% NewValue: handle of the currently selected object
% handles structure with handles and user data (see GUIDATA)
global matP matQ matdesvP matdesvQ matfilteredP matfilteredQ
switch get(hObject,'Tag') % Get Tag of selected object
case 'radiobutton3'
[valores horas]=size(matP);
clear horas;
horamenu=get(handles.popupmenu1,'Value');
set(handles.text5,'String',num2str(matdesvP(horamenu)));
axes(handles.axes3);
xrange=[1:1:valores];
fill([xrange(1) xrange(end) xrange(end) xrange(1)],[matfilteredP(1,horamenu)-matdesvP(horamenu) matfilteredP(1,horamenu)-matdesvP(horamenu) matfilteredP(1,horamenu)+matdesvP(horamenu) matfilteredP(1,horamenu)+matdesvP(horamenu)],[0.9 0.9 1]);hold on;
scatter(xrange,matP(:,horamenu),'filled','bs');xlabel('Muestra No.');ylabel('MW','fontsize',7);grid;set(gca,'Layer','top');hold on;
set(handles.axes3,'xlim',[1 valores]);
line([xrange(1) xrange(end)],[matfilteredP(1,horamenu) matfilteredP(1,horamenu)],'LineWidth',2,'LineStyle','-','Color',[1 0 0]);hold off;
case 'radiobutton5'
[valores horas]=size(matQ);
clear horas;
horamenu=get(handles.popupmenu1,'Value');
set(handles.text5,'String',num2str(matdesvQ(horamenu)));
axes(handles.axes3);
xrange=[1:1:valores];
fill([xrange(1) xrange(end) xrange(end) xrange(1)],[matfilteredQ(1,horamenu)-matdesvQ(horamenu) matfilteredQ(1,horamenu)-matdesvQ(horamenu) matfilteredQ(1,horamenu)+matdesvQ(horamenu) matfilteredQ(1,horamenu)+matdesvQ(horamenu)],[0.9 0.9 1]);hold on;
scatter(xrange,matQ(:,horamenu),'filled','bs');xlabel('Muestra No.');ylabel('Mvar','fontsize',7);grid;set(gca,'Layer','top');hold on;
set(handles.axes3,'xlim',[1 valores]);
line([xrange(1) xrange(end)],[matfilteredQ(1,horamenu) matfilteredQ(1,horamenu)],'LineWidth',2,'LineStyle','-','Color',[1 0 0]);hold off;
end
% --- Executes when selected object is changed in uipanel10.
function uipanel10_SelectionChangeFcn(hObject, eventdata, handles)
% hObject handle to the selected object in uipanel10
% eventdata structure with the following fields (see UIBUTTONGROUP)
% EventName: string 'SelectionChanged' (read only)
% OldValue: handle of the previously selected object or empty if none was selected
% NewValue: handle of the currently selected object
% handles structure with handles and user data (see GUIDATA)
global datosporhoraP datosporhoraQ
switch get(hObject,'Tag') % Get Tag of selected object
case 'radiobutton6'
set(handles.uitable3,'Data',datosporhoraP);
case 'radiobutton7'
set(handles.uitable3,'Data',datosporhoraQ);
end
% --- Executes when selected object is changed in uipanel12.
function uipanel12_SelectionChangeFcn(hObject, eventdata, handles)
% hObject handle to the selected object in uipanel12
% eventdata structure with the following fields (see UIBUTTONGROUP)
% EventName: string 'SelectionChanged' (read only)
% OldValue: handle of the previously selected object or empty if none was selected
% NewValue: handle of the currently selected object
% handles structure with handles and user data (see GUIDATA)
global matfilteredP matfilteredQ mes
horas=[1:1:25];
switch get(hObject,'Tag') % Get Tag of selected object
case 'radiobutton8'
axes(handles.axes2);
set(handles.uipanel2,'Title','Curva de Carga Promedio:');
horas=[1:1:25];
activa = [];
reactiva = [];
if ismember(mes,[3 4 5 6 7 8 9 10])
% Es Verano, el pico va entre las horas 19 y 20
activa = cat(2,matfilteredP(1:19),matfilteredP(25));
activa = cat(2,activa,matfilteredP(20:24));
reactiva = cat(2,matfilteredQ(1:19),matfilteredQ(25));
reactiva = cat(2,reactiva,matfilteredQ(20:24));
plot(horas, activa,'r.-',horas,reactiva,'b.-');xlabel('Hora');legend('P(MW)','Location','NorthWest','Q(Mvar)','Location','NorthWest');grid;
set(gca,'XTick',[0:1:25]);
set(gca,'XTickLabel',{'0';'';'';'';'';'5';'';'';'';'';'10';'';'';'';'';'15';'';'';'';'';'';'20';'';'';'';'24'});
else
% Es Invierno, el pico va entre las horas 18 y 19
activa = cat(2,matfilteredP(1:18),matfilteredP(25));
activa = cat(2,activa,matfilteredP(19:24));
reactiva = cat(2,matfilteredQ(1:18),matfilteredQ(25));
reactiva = cat(2,reactiva,matfilteredQ(19:24));
plot(horas, activa,'r.-',horas,reactiva,'b.-');xlabel('Hora');legend('P(MW)','Location','NorthWest','Q(Mvar)','Location','NorthWest');grid;
set(gca,'XTick',[0:1:25]);
set(gca,'XTickLabel',{'0';'';'';'';'';'5';'';'';'';'';'10';'';'';'';'';'15';'';'';'';'';'';'20';'';'';'';'24'});
end
case 'radiobutton9'
axes(handles.axes2);
set(handles.uipanel2,'Title','Curva de Carga Promedio:');
horas=[1:1:25];
activa = [];
reactiva = [];
if ismember(mes,[3 4 5 6 7 8 9 10])
% Es Verano, el pico va entre las horas 19 y 20
activa = cat(2,matfilteredP(1:19),matfilteredP(25));
activa = cat(2,activa,matfilteredP(20:24));
reactiva = cat(2,matfilteredQ(1:19),matfilteredQ(25));
reactiva = cat(2,reactiva,matfilteredQ(20:24));
stairs(horas, activa,'r.-');hold on;stairs(horas,reactiva,'b.-');hold off;xlabel('Hora');legend('P(MW)','Location','NorthWest','Q(Mvar)','Location','NorthWest');grid;
set(gca,'XTick',[0:1:25]);
set(gca,'XTickLabel',{'0';'';'';'';'';'5';'';'';'';'';'10';'';'';'';'';'15';'';'';'';'';'';'20';'';'';'';'24'});
else
% Es Invierno, el pico va entre las horas 18 y 19
activa = cat(2,matfilteredP(1:18),matfilteredP(25));
activa = cat(2,activa,matfilteredP(19:24));
reactiva = cat(2,matfilteredQ(1:18),matfilteredQ(25));
reactiva = cat(2,reactiva,matfilteredQ(19:24));
stairs(horas, activa,'r.-');hold on;stairs(horas,reactiva,'b.-');hold off;xlabel('Hora');legend('P(MW)','Location','NorthWest','Q(Mvar)','Location','NorthWest');grid;
set(gca,'XTick',[0:1:25]);
set(gca,'XTickLabel',{'0';'';'';'';'';'5';'';'';'';'';'10';'';'';'';'';'15';'';'';'';'';'';'20';'';'';'';'24'});
end
%stairs(horas, matfilteredP,'r-');hold on;stairs(horas,matfilteredQ,'b-');hold off;xlabel('Hora');legend('P(MW)','Location','NorthWest','Q(Mvar)','Location','NorthWest');grid;
case 'radiobutton10'
% Dibuja la curva del FP
axes(handles.axes2);
set(handles.uipanel2,'Title','Curva del Factor de Potencia Promedio:');
matFP = matfilteredP./sqrt(matfilteredP.^2 + matfilteredQ.^2);
horas = [1:1:25];
fp = [];
if ismember(mes,[3 4 5 6 7 8 9 10])
% Es Verano, el pico va entre las horas 19 y 20
fp = cat(2,matFP(1:19),matFP(25));
fp = cat(2,fp,matFP(20:24));
plot(horas, fp,'k.-');xlabel('Hora');legend('FP','Location','NorthWest');grid;
set(gca,'XTick',[0:1:25]);
set(gca,'XTickLabel',{'0';'';'';'';'';'5';'';'';'';'';'10';'';'';'';'';'15';'';'';'';'';'';'20';'';'';'';'24'});
else
% Es Invierno, el pico va entre las horas 18 y 19
fp = cat(2,matFP(1:18),matFP(25));
fp = cat(2,fp,matFP(19:24));
plot(horas, fp,'k.-');xlabel('Hora');legend('FP','Location','NorthWest');grid;
set(gca,'XTick',[0:1:25]);
set(gca,'XTickLabel',{'0';'';'';'';'';'5';'';'';'';'';'10';'';'';'';'';'15';'';'';'';'';'';'20';'';'';'';'24'});
end
end
% --------------------------------------------------------------------
function histfreq_Callback(hObject, eventdata, handles)
% hObject handle to histfreq (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global orgmatP orgmatQ
potgrupo = get(get(handles.uipanel8,'SelectedObject'),'Tag');
switch potgrupo
case 'radiobutton3'
figure('Name','Histograma de Frecuencias','NumberTitle','off');hist(orgmatP(:,get(handles.popupmenu1,'Value')));
case 'radiobutton5'
figure('Name','Histograma de Frecuencias','NumberTitle','off');hist(orgmatQ(:,get(handles.popupmenu1,'Value')));
end
% --------------------------------------------------------------------
function histdensity_Callback(hObject, eventdata, handles)
% hObject handle to histdensity (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global orgmatP orgmatQ
potgrupo = get(get(handles.uipanel8,'SelectedObject'),'Tag');
switch potgrupo
case 'radiobutton3'
dfittool(orgmatP(:,get(handles.popupmenu1,'Value')));
case 'radiobutton5'
dfittool(orgmatQ(:,get(handles.popupmenu1,'Value')));
end