-
Notifications
You must be signed in to change notification settings - Fork 0
/
TC1_draft.m
123 lines (94 loc) · 2.86 KB
/
TC1_draft.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
% DISCIPLINA: TIP7044 - ESTIMAÇÃO E IDENTIFICAÇÃO DE SISTEMAS (2022.2 - T01)
% PROFESSOR: GUILHERME DE ALENCAR BARRETO
% ALUNO(A): ANDREZA COSTA NASCIMENTO
% MATRÍCULA: 544957
% 1º TRABALHO COMPUTACIONAL
%% Parte 1 - Método Alternativo
clc;clear;
serie_vazoes = load('furnas.dat');
m_serie = mean(serie_vazoes);
dp_serie = std(serie_vazoes);
modulo_serie = reshape((serie_vazoes./m_serie)', 1,[]);
serie_vazoes = reshape((serie_vazoes)', 1,[]);
fig = figure; clf
h = plot(serie_vazoes); hold on
grid on
% define figure properties
opts.Colors = get(groot,'defaultAxesColorOrder');
opts.saveFolder = 'img/';
opts.width = 20;
opts.height = 10;
opts.fontType = 'Times';
opts.fontSize = 24.6;
% add axis labes and legend
axis tight
title('Série de vazões Furnas de Jan/1930 até Dez/1990 com Sazonalidade')
xlabel('Mês')
ylabel('Vazão média do reservatório (m³)')
%legend('1960','1946','1970','1972','1975')
% scaling
fig.Units = 'centimeters';
fig.Position(3) = opts.width;
fig.Position(4) = opts.height;
% set text properties
set(fig.Children, ...
'FontName', 'Times', ...
'FontSize', 23);
% remove unnecessary white space
set(gca,'LooseInset',max(get(gca,'TightInset'), 0.02))
hold off
fig = figure; clf
h = plot(modulo_serie); hold on
grid on
% define figure properties
opts.Colors = get(groot,'defaultAxesColorOrder');
opts.saveFolder = 'img/';
opts.width = 20;
opts.height = 10;
opts.fontType = 'Times';
opts.fontSize = 24.6;
% add axis labes and legend
axis tight
title('Série de vazões Furnas de Jan/1930 até Dez/1990 sem Sazonalidade')
xlabel('Mês')
ylabel('Vazão média do reservatório (m³)')
%legend('1960','1946','1970','1972','1975')
% scaling
fig.Units = 'centimeters';
fig.Position(3) = opts.width;
fig.Position(4) = opts.height;
% set text properties
set(fig.Children, ...
'FontName', 'Times', ...
'FontSize', 23);
% remove unnecessary white space
set(gca,'LooseInset',max(get(gca,'TightInset'), 0.02))
hold off
%HISTOGRAMA
fig = figure; clf
h = histogram(modulo_serie); hold on
grid on
% define figure properties
opts.Colors = get(groot,'defaultAxesColorOrder');
opts.saveFolder = 'img/';
opts.width = 20;
opts.height = 10;
opts.fontType = 'Times';
opts.fontSize = 24.6;
% add axis labes and legend
axis tight
title('Histograma da Série de Vazões Furnas sem Sazonalidade')
%xlabel('Vazão média do reservatório (m³)')
ylabel('Frequência (Meses)')
%legend('1960','1946','1970','1972','1975')
% scaling
fig.Units = 'centimeters';
fig.Position(3) = opts.width;
fig.Position(4) = opts.height;
% set text properties
set(fig.Children, ...
'FontName', 'Times', ...
'FontSize', 23);
% remove unnecessary white space
set(gca,'LooseInset',max(get(gca,'TightInset'), 0.02))
hold off