-
Notifications
You must be signed in to change notification settings - Fork 0
/
AMvsYears.m
69 lines (63 loc) · 2.71 KB
/
AMvsYears.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
% |**********************************************************************;
% * Project : Examinations of Assumptions of
% * Environmental Load Models
% *
% * Program name : AMvsYears.m
% *
% * Author : Maulin Amin
% *
% * Date created : November 22, 2018
% *
% * Purpose : This program plots the values of Annual Maximum
% Wind Speed vs Year.
%
% * Note :
% *
% * Revision History :
% *
% * Date Author Ref Revision (Date in YYYYMMDD format)
% * YYYY/MM/DD M Amin 1 Description of Revision
% *
% |**********************************************************************;
%==================================
% Basic Initialization instructions
%==================================
clear; % Clear Memory
clc; %Clear Command Window
%===================================
% Inputs
%===================================
% Location of Trenton CSV files obtained from Environment Canada
folder='C:\Users\Maulin Amin\OneDrive - University of Waterloo\Waterloo\Winter 2018\Environment Canada\Wind&SnowData\CSV\Trenton';
%===================================
% Outputs
%===================================
data = CollectMaxDataPointsForAStation(folder);
figure(1);
plot(data.Date_Time,data.SpdOfMaxGust_km_h_);
ylabel({'Annual Maximum Wind Speed'}); xlabel('Year');
% set(gca,'Ylim',[120 180],'Xlim',[2018 2027]) % Adjust Y limits of "current axes"
set(gca,'FontName','Times');
set(gcf,'Units','inches') % Set figure size units of "current figure"
set(gcf,'Color','white');
set(gcf,'Position',[0,0,3,2.5]) % Set figure width (6 in.) and height (4 in.)
print -deps2c 602.eps % Save as PDF
movefile('602.eps','C:\Users\Maulin Amin\OneDrive - University of Waterloo\Waterloo\Winter 2018\Environment Canada\Thesis\Latex\plots');
years = [1955:2018]';
data.years = years;
figure(2);
p = polyfit(data.years,data.SpdOfMaxGust_km_h_,1);
f = polyval(p,data.years);
plot(data.years,data.SpdOfMaxGust_km_h_,'.',data.years,f,'-')
grid on;
dim = [0.2 0.5 0.3 0.3];
mdl = fitlm(data.years,data.SpdOfMaxGust_km_h_);
ylabel('Annual Maximum Wind Speed'); xlabel('Years');
legend('Wind Speed','Linear Fit','Location','northeast')
set(gca,'Xlim',[1955 2018]) % Adjust Y limits of "current axes"
set(gca,'FontName','Times');
set(gcf,'Units','inches') % Set figure size units of "current figure"
set(gcf,'Color','white');
set(gcf,'Position',[0,0,3,2.5]) % Set figure width (6 in.) and height (4 in.)
print -deps2c 603.eps % Save as PDF
movefile('603.eps','C:\Users\Maulin Amin\OneDrive - University of Waterloo\Waterloo\Winter 2018\Environment Canada\Thesis\Latex\plots');