# Solar irradiance forecasting for Chulalongkorn University location using time series models

Save this PDF as:

Size: px
Start display at page:

## Transcription

1 Senior Project Proposal Year 2016 Solar irradiance forecasting for Chulalongkorn University location using time series models Vichaya Layanun ID Advisor: Assist. Prof. Jitkomut Songsiri Department of Electrical Engineering Faculty of Engineering, Chulalongkorn University

2 Why solar forecasting is important? Because, the unreliability of solar power generating has made the entire electrical power generation difficult for power management. Solar forecasting is a widely-common approach to deal with the problem. An improved accuracy in the forecast can provide a better management of electrical power production. 2

3 Forecasting Models Persistence Forecasts Time Series Models Numerical Weather Prediction (NWP) Artificial Intelligence (AI) Stationary Models Non-stationary Models AR MA ARIMA ARIMAX ARMA ARMAX There are many forecasting models to predict the future solar irradiance. We focus on time series models. 3

4 Objective To study the relevant variables of solar irradiance forecasting. To apply ARIMA models, a Seasonal ARIMA models and ARIMAX models to forecast solar irradiance. To validate results of forecasting performance among the models using RMSE, MAE and a sample autocovarince function of the residual as model validation criterion. To solve the practical issues on data pre-processing. 4

5 Essential Elements Relevant variable Forecasting Horizon Forecasting Performance Evaluation Measures 5

6 Relevant Variables Global Horizontal Irradiance (GHI) is the considered variable which is the geometric sum of Direct Normal Irradiance (DNI) and Diffuse Horizontal Irradiance (DHI). This study predicted GHI and we will use the notation of I(t) throughout this project. 6

7 Forecasting Horizon To decide how far to do prediction, forecasting horizon is the key. Using time series model for solar forecasting is suitable and widely used for very-short-term forecasting. Different forecasting horizons lead to different suitable forecasting method. Terms of forecasting horizon. Forecasting Horizon Time Vert-short Short Medium Long From a few second to 6 hours Up to hours ahead Up to week ahead Up to month or year 7

8 Forecasting Performance Evaluation Measures Common evaluation measures are Root Mean Squared Error (RMSE) and Mean Absolute Error (MAE) which are used to validate a forecasting method. RMSE and MAE can be defined as RMSE = N 1 N (I t I (t)) 2 t=1 (1) N MAE = 1 N I t I (t) t=1 (2) where N is the length of the time horizon. Desired value of RMSE and MAE is minimized. 8

9 Forecasting Performance Evaluation Measures The residual analysis can also be used for validation. The residuals is defined as e t = I t I (t) (3) If I (t) can capture the dynamic in the true model, then e t should behave like white noise. Because we do not know the true model, a sample autocovariance function (ACV) of the residual is determined. ACV of residual is defined as R e τ = 1 N N e t e t τ. t= τ (4) Because ACV of white noise is zero after τ > 0, a desired model is the model of which R e τ is small for τ > 0. 9

10 Time Series Models ARIMAX description Estimation of integrated term Seasonal ARIMA model 10

11 ARIMAX description An autoregressive integrated moving average model with an exogenous input (ARIMAX) is employed to predict the future solar irradiance. The model ARIMAX(p,d,q) is defined by A q 1 1 q 1 d y t = B q 1 u t + C q 1 v(t) (5) where A q 1 = I (a 1 q 1 + a 2 q a p q p ) B q 1 C q 1 = B 1 q 1 + B 2 q B m q p = I + c 1 q 1 + c 2 q c q q p q 1 is lag operator, y t is the solar irradiance (GHI), u t are exogenous inputs, v t is white noise with zero mean and σ 2 variance. In this project, the exogenous inputs consist of the local temperature, relative humidity, wind speed, and air pressure. (1 q 1 ) d is integrated term where d is degree of differencing. 11

12 Estimation of integrated term A q 1 y t = C q 1 v(t) (6) A q 1 1 q 1 d y t = C q 1 v(t) (12) After differencing for certain order, if the resulting ACF behaves similarly to that of a stationary process in the table, then the actual process is ARIMA. In other words, if ACF decrease apparently slowly, the transfer function has pole lying near or on unit circle. Hence, the model needs more differencing. 12

13 Seasonal ARIMA model Decomposition of Solar Irradiance in January 2014 in Bangkok. 13

14 Seasonal ARIMA model GHI (W/m 2 ) Time GHI after removing seasonal trend in January 2014 in Bangkok. We imply that GHI can be described as ARMA models containing s season: A q 1 y t = s t + α + C q 1 v(t) (7) 14

15 Seasonal ARIMA model In this study, we used a seasonal ARIMA models to remove a seasonal trend. The seasonal term and constant are removed by using this transformation. This method is called a Seasonal ARIMA (p, d, q)(p, D, Q) T models which can be defined as A q T A q 1 (1 q T ) D 1 q 1 d y t = C q T C q 1 v(t) (8) where A q 1 = I (a 1 q T + a 2 q 2T + + a p q pp ) C q 1 = I + c 1q T + c 2q 2T + + c qq qq T is a seasonal period and D is integrated seasonal order. 15

16 Preliminary Results Model Selection Model Validation 16

17 Model selection The goal of the experiment is to apply the Akaike information criterion (AIC) and Bayesian information criterion (BIC) to find the optimum order of ARIMA models and Seasonal ARIMA models. AIC and BIC are measures of the performance of the candidate models. The best performance model gives the lowest AIC or BIC. 17

18 AIC scores of the candidate models 18

19 Model validation After we select the order of the candidate models, we would compare their performance on the unseen data sets. We used RMSE and MAE to validate the performance of the candidate models. The validation data is historical GHI data in year

20 Comparison among the models 20

21 Comparison among the models GHI (W/m 2 ) Time 21

22 Comparison among the models 22

23 Project overview Scope of work Expected outcomes 23

24 Scope of work We focus on solar irradiance in the area of Chulalongkorn University location by using ARIMAX and seasonal ARIMAX models. The exogenous inputs consist of the local temperature, relative humidity, wind speed and air pressure will be included in the models. We will conduct experiments to verify our approach by using data obtained from Thai Meteorological Department (TMD). We will conduct an experiment by using data obtained from department of electrical engineering if they are available. We will study the consequence of seasons of the year in the location. 24

25 Expected outcomes Comparison results of forecasting performance among Persistence forecast, ARIMA models, a Seasonal ARIMA models and ARIMAX models using Root Mean Squared Error (RMSE), Mean Absolute Error and a sample autocovariance function (ACV) of the residual as model validation criterion. Schemes for solving practical issues on data pre-processing which are 1) missing data and 2) timely asynchronous data. 25

26 Q&A 26

27 Backup 27

28 Relevant Variables Global Horizontal Irradiance (GHI) is the considered variable which is the geometric sum of Direct Normal Irradiance (DNI) and Diffuse Horizontal Irradiance (DHI). GHI = DHI + DNI cos(θ) (9) where θ is the solar zenith angle. The unit of GHI, DNI and DHI are W/m 2. This study predicted GHI and we will use the notation of I(t) throughout this project. 28

29 Persistence Forecasts Persistence forecast as a baseline prediction method is used for a comparison to more advanced methods. There are many equations of persistence forecast in accordance with each study. I (t + h) = I(t) (10) 1) Persistence model supposes that solar irradiance at time t + h can be predicted by it value at time t. 29

30 Persistence Forecasts 2) Clearness persistence forecast supposes that clearness index is constant, i.e. K t + h = K t, where clearness index is defined as : K t = I(t) I EX (t). (11) I EX (t) = I 0 cosθ(t) where I 0 is an average value of the extraterrestrial irradiance and is suggested to be 1360 and θ(t) is solar zenith angle at time t. Thus, I (t + h) = K t I EX (t + h) (12) 30

31 Persistence Forecasts 2) Clear sky persistence forecast supposes that clear sky index is constant, i.e. k t + h = k t, where clearness index is defined as : k t = I(t) I clr (t). (13) I clr (t) is the value from a clear sky irradiance model at time t which can be widely estimated by Ineichen model. Thus, I (t + h) = k t I clr (t + h) (14) 31

32 Seasonal ARIMA model An estimate of the daily seasonal effect s t can be obtained by averaging s (t) which is obtained by subtracting m t : s t = N 1 s (t + kk) N + 1 k=0 where T is the seasonal period, s t = y t m (t) m (t) = 1 2 y t T 2 + y t T y t + T y t + T 2 T 32

33 Seasonal ARIMA model In this study, we used a seasonal ARIMA models to remove a seasonal trend. Giving the following equation as an example: We subtract y(t) by y(t T). y t = s t + α + v(t) (15) where s(t) = s(t kk), k is an integer and α is a constant. y t y t T y t y t T = s t + α + v t s t T α - v t T = v t v t T 33

34 Model Selection The AIC is defined as The BIC is defined as AIC = 2L + 2d (16) where L is the loglikelyhood function and d is the number of effective parameters. where N is the number of sample. BIC = 2L + dlogn (17) 34

35 Autocorrelation Function An autocorrelation function is defined by ACF = R(τ) R(0) (18) where the sample autocovariance function (ACVF) is defined by R(τ) = 1 N N y t y(t τ) t=τ (19) where the expectation of y is zero and the variance of y which are the same at all time t. 35

36 BIC scores of the candidate models 36

37 The price of a solar panel have obviously dropped from \$ in 1975 to \$0.447 in 2016 and global solar panel installations have risen from 2 MW in 1975 to 64,892 MW in Source:

38 Total installed renewable power plant capacity from AEDP in Thailand. 38

### Solar irradiance forecasting for Chulalongkorn University location using time series models

Senior Project Proposal 2102499 Year 2016 Solar irradiance forecasting for Chulalongkorn University location using time series models Vichaya Layanun ID 5630550721 Advisor: Assist. Prof. Jitkomut Songsiri

### Lecture 2: Univariate Time Series

Lecture 2: Univariate Time Series Analysis: Conditional and Unconditional Densities, Stationarity, ARMA Processes Prof. Massimo Guidolin 20192 Financial Econometrics Spring/Winter 2017 Overview Motivation:

### Probabilistic forecasting of solar radiation

Probabilistic forecasting of solar radiation Dr Adrian Grantham School of Information Technology and Mathematical Sciences School of Engineering 7 September 2017 Acknowledgements Funding: Collaborators:

### TIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA

CHAPTER 6 TIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA 6.1. Introduction A time series is a sequence of observations ordered in time. A basic assumption in the time series analysis

### The Center for Renewable Resource Integration at UC San Diego

The Center for Renewable Resource Integration at UC San Diego Carlos F. M. Coimbra ccoimbra@ucsd.edu; solarwind.ucsd.edu Jan Kleissl and Byron Washom UCSD Center of Excellence in Renewable Resources and

### An Adaptive Multi-Modeling Approach to Solar Nowcasting

An Adaptive Multi-Modeling Approach to Solar Nowcasting Antonio Sanfilippo, Luis Martin-Pomares, Nassma Mohandes, Daniel Perez-Astudillo, Dunia A. Bachour ICEM 2015 Overview Introduction Background, problem

### Solar Nowcasting with Cluster-based Detrending

Solar Nowcasting with Cluster-based Detrending Antonio Sanfilippo, Luis Pomares, Daniel Perez-Astudillo, Nassma Mohandes, Dunia Bachour ICEM 2017 Oral Presentation 26-29June 2017, Bari, Italy Overview

### Forecasting using R. Rob J Hyndman. 2.4 Non-seasonal ARIMA models. Forecasting using R 1

Forecasting using R Rob J Hyndman 2.4 Non-seasonal ARIMA models Forecasting using R 1 Outline 1 Autoregressive models 2 Moving average models 3 Non-seasonal ARIMA models 4 Partial autocorrelations 5 Estimation

### Prashant Pant 1, Achal Garg 2 1,2 Engineer, Keppel Offshore and Marine Engineering India Pvt. Ltd, Mumbai. IJRASET 2013: All Rights are Reserved 356

Forecasting Of Short Term Wind Power Using ARIMA Method Prashant Pant 1, Achal Garg 2 1,2 Engineer, Keppel Offshore and Marine Engineering India Pvt. Ltd, Mumbai Abstract- Wind power, i.e., electrical

### Time Series I Time Domain Methods

Astrostatistics Summer School Penn State University University Park, PA 16802 May 21, 2007 Overview Filtering and the Likelihood Function Time series is the study of data consisting of a sequence of DEPENDENT

### Univariate ARIMA Models

Univariate ARIMA Models ARIMA Model Building Steps: Identification: Using graphs, statistics, ACFs and PACFs, transformations, etc. to achieve stationary and tentatively identify patterns and model components.

### Firstly, the dataset is cleaned and the years and months are separated to provide better distinction (sample below).

Project: Forecasting Sales Step 1: Plan Your Analysis Answer the following questions to help you plan out your analysis: 1. Does the dataset meet the criteria of a time series dataset? Make sure to explore

### Time Series Analysis -- An Introduction -- AMS 586

Time Series Analysis -- An Introduction -- AMS 586 1 Objectives of time series analysis Data description Data interpretation Modeling Control Prediction & Forecasting 2 Time-Series Data Numerical data

### FORECASTING SUGARCANE PRODUCTION IN INDIA WITH ARIMA MODEL

FORECASTING SUGARCANE PRODUCTION IN INDIA WITH ARIMA MODEL B. N. MANDAL Abstract: Yearly sugarcane production data for the period of - to - of India were analyzed by time-series methods. Autocorrelation

### 2. An Introduction to Moving Average Models and ARMA Models

. An Introduction to Moving Average Models and ARMA Models.1 White Noise. The MA(1) model.3 The MA(q) model..4 Estimation and forecasting of MA models..5 ARMA(p,q) models. The Moving Average (MA) models

### Some Time-Series Models

Some Time-Series Models Outline 1. Stochastic processes and their properties 2. Stationary processes 3. Some properties of the autocorrelation function 4. Some useful models Purely random processes, random

### NANYANG TECHNOLOGICAL UNIVERSITY SEMESTER II EXAMINATION MAS451/MTH451 Time Series Analysis TIME ALLOWED: 2 HOURS

NANYANG TECHNOLOGICAL UNIVERSITY SEMESTER II EXAMINATION 2012-2013 MAS451/MTH451 Time Series Analysis May 2013 TIME ALLOWED: 2 HOURS INSTRUCTIONS TO CANDIDATES 1. This examination paper contains FOUR (4)

### A Comparison of the Forecast Performance of. Double Seasonal ARIMA and Double Seasonal. ARFIMA Models of Electricity Load Demand

Applied Mathematical Sciences, Vol. 6, 0, no. 35, 6705-67 A Comparison of the Forecast Performance of Double Seasonal ARIMA and Double Seasonal ARFIMA Models of Electricity Load Demand Siti Normah Hassan

### Statistics 349(02) Review Questions

Statistics 349(0) Review Questions I. Suppose that for N = 80 observations on the time series { : t T} the following statistics were calculated: _ x = 10.54 C(0) = 4.99 In addition the sample autocorrelation

### INTRODUCTION TO TIME SERIES ANALYSIS. The Simple Moving Average Model

INTRODUCTION TO TIME SERIES ANALYSIS The Simple Moving Average Model The Simple Moving Average Model The simple moving average (MA) model: More formally: where t is mean zero white noise (WN). Three parameters:

### FE570 Financial Markets and Trading. Stevens Institute of Technology

FE570 Financial Markets and Trading Lecture 5. Linear Time Series Analysis and Its Applications (Ref. Joel Hasbrouck - Empirical Market Microstructure ) Steve Yang Stevens Institute of Technology 9/25/2012

### Short-Term Load Forecasting Using ARIMA Model For Karnataka State Electrical Load

International Journal of Engineering Research and Development e-issn: 2278-67X, p-issn: 2278-8X, www.ijerd.com Volume 13, Issue 7 (July 217), PP.75-79 Short-Term Load Forecasting Using ARIMA Model For

### Empirical Approach to Modelling and Forecasting Inflation in Ghana

Current Research Journal of Economic Theory 4(3): 83-87, 2012 ISSN: 2042-485X Maxwell Scientific Organization, 2012 Submitted: April 13, 2012 Accepted: May 06, 2012 Published: June 30, 2012 Empirical Approach

### MCMC analysis of classical time series algorithms.

MCMC analysis of classical time series algorithms. mbalawata@yahoo.com Lappeenranta University of Technology Lappeenranta, 19.03.2009 Outline Introduction 1 Introduction 2 3 Series generation Box-Jenkins

### Analysis. Components of a Time Series

Module 8: Time Series Analysis 8.2 Components of a Time Series, Detection of Change Points and Trends, Time Series Models Components of a Time Series There can be several things happening simultaneously

### Transformations for variance stabilization

FORECASTING USING R Transformations for variance stabilization Rob Hyndman Author, forecast Variance stabilization If the data show increasing variation as the level of the series increases, then a transformation

### { } Stochastic processes. Models for time series. Specification of a process. Specification of a process. , X t3. ,...X tn }

Stochastic processes Time series are an example of a stochastic or random process Models for time series A stochastic process is 'a statistical phenomenon that evolves in time according to probabilistic

### SOLAR POWER FORECASTING BASED ON NUMERICAL WEATHER PREDICTION, SATELLITE DATA, AND POWER MEASUREMENTS

BASED ON NUMERICAL WEATHER PREDICTION, SATELLITE DATA, AND POWER MEASUREMENTS Detlev Heinemann, Elke Lorenz Energy Meteorology Group, Institute of Physics, Oldenburg University Workshop on Forecasting,

### Box-Jenkins ARIMA Advanced Time Series

Box-Jenkins ARIMA Advanced Time Series www.realoptionsvaluation.com ROV Technical Papers Series: Volume 25 Theory In This Issue 1. Learn about Risk Simulator s ARIMA and Auto ARIMA modules. 2. Find out

### Lab: Box-Jenkins Methodology - US Wholesale Price Indicator

Lab: Box-Jenkins Methodology - US Wholesale Price Indicator In this lab we explore the Box-Jenkins methodology by applying it to a time-series data set comprising quarterly observations of the US Wholesale

### Modelling using ARMA processes

Modelling using ARMA processes Step 1. ARMA model identification; Step 2. ARMA parameter estimation Step 3. ARMA model selection ; Step 4. ARMA model checking; Step 5. forecasting from ARMA models. 33

### Importance of Numerical Weather Prediction in Variable Renewable Energy Forecast

Importance of Numerical Weather Prediction in Variable Renewable Energy Forecast Dr. Abhijit Basu (Integrated Research & Action for Development) Arideep Halder (Thinkthrough Consulting Pvt. Ltd.) September

### Forecasting Area, Production and Yield of Cotton in India using ARIMA Model

Forecasting Area, Production and Yield of Cotton in India using ARIMA Model M. K. Debnath 1, Kartic Bera 2 *, P. Mishra 1 1 Department of Agricultural Statistics, Bidhan Chanda Krishi Vishwavidyalaya,

### Lecture 3: Autoregressive Moving Average (ARMA) Models and their Practical Applications

Lecture 3: Autoregressive Moving Average (ARMA) Models and their Practical Applications Prof. Massimo Guidolin 20192 Financial Econometrics Winter/Spring 2018 Overview Moving average processes Autoregressive

### Automatic Forecasting

Automatic Forecasting Summary The Automatic Forecasting procedure is designed to forecast future values of time series data. A time series consists of a set of sequential numeric data taken at equally

### Classic Time Series Analysis

Classic Time Series Analysis Concepts and Definitions Let Y be a random number with PDF f Y t ~f,t Define t =E[Y t ] m(t) is known as the trend Define the autocovariance t, s =COV [Y t,y s ] =E[ Y t t

Advanced Econometrics Marco Sunder Nov 04 2010 Marco Sunder Advanced Econometrics 1/ 25 Contents 1 2 3 Marco Sunder Advanced Econometrics 2/ 25 Music Marco Sunder Advanced Econometrics 3/ 25 Music Marco

### Minitab Project Report - Assignment 6

.. Sunspot data Minitab Project Report - Assignment Time Series Plot of y Time Series Plot of X y X 7 9 7 9 The data have a wavy pattern. However, they do not show any seasonality. There seem to be an

### Scenario 5: Internet Usage Solution. θ j

Scenario : Internet Usage Solution Some more information would be interesting about the study in order to know if we can generalize possible findings. For example: Does each data point consist of the total

### Suan Sunandha Rajabhat University

Forecasting Exchange Rate between Thai Baht and the US Dollar Using Time Series Analysis Kunya Bowornchockchai Suan Sunandha Rajabhat University INTRODUCTION The objective of this research is to forecast

### Ch 6. Model Specification. Time Series Analysis

We start to build ARIMA(p,d,q) models. The subjects include: 1 how to determine p, d, q for a given series (Chapter 6); 2 how to estimate the parameters (φ s and θ s) of a specific ARIMA(p,d,q) model (Chapter

### Estimation and application of best ARIMA model for forecasting the uranium price.

Estimation and application of best ARIMA model for forecasting the uranium price. Medeu Amangeldi May 13, 2018 Capstone Project Superviser: Dongming Wei Second reader: Zhenisbek Assylbekov Abstract This

### The ARIMA Procedure: The ARIMA Procedure

Page 1 of 120 Overview: ARIMA Procedure Getting Started: ARIMA Procedure The Three Stages of ARIMA Modeling Identification Stage Estimation and Diagnostic Checking Stage Forecasting Stage Using ARIMA Procedure

### STAT 436 / Lecture 16: Key

STAT 436 / 536 - Lecture 16: Key Modeling Non-Stationary Time Series Many time series models are non-stationary. Recall a time series is stationary if the mean and variance are constant in time and the

### EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 011 MODULE 3 : Stochastic processes and time series Time allowed: Three Hours Candidates should answer FIVE questions. All questions carry

### A SARIMAX coupled modelling applied to individual load curves intraday forecasting

A SARIMAX coupled modelling applied to individual load curves intraday forecasting Frédéric Proïa Workshop EDF Institut Henri Poincaré - Paris 05 avril 2012 INRIA Bordeaux Sud-Ouest Institut de Mathématiques

### Analysis of Violent Crime in Los Angeles County

Analysis of Violent Crime in Los Angeles County Xiaohong Huang UID: 004693375 March 20, 2017 Abstract Violent crime can have a negative impact to the victims and the neighborhoods. It can affect people

### Empirical Market Microstructure Analysis (EMMA)

Empirical Market Microstructure Analysis (EMMA) Lecture 3: Statistical Building Blocks and Econometric Basics Prof. Dr. Michael Stein michael.stein@vwl.uni-freiburg.de Albert-Ludwigs-University of Freiburg

### Modelling Wind Farm Data and the Short Term Prediction of Wind Speeds

Modelling Wind Farm Data and the Short Term Prediction of Wind Speeds An Investigation into Wind Speed Data Sets Erin Mitchell Lancaster University 6th April 2011 Outline 1 Data Considerations Overview

### Topic 4 Unit Roots. Gerald P. Dwyer. February Clemson University

Topic 4 Unit Roots Gerald P. Dwyer Clemson University February 2016 Outline 1 Unit Roots Introduction Trend and Difference Stationary Autocorrelations of Series That Have Deterministic or Stochastic Trends

### Comparing the Univariate Modeling Techniques, Box-Jenkins and Artificial Neural Network (ANN) for Measuring of Climate Index

Applied Mathematical Sciences, Vol. 8, 2014, no. 32, 1557-1568 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.4150 Comparing the Univariate Modeling Techniques, Box-Jenkins and Artificial

### 5 Autoregressive-Moving-Average Modeling

5 Autoregressive-Moving-Average Modeling 5. Purpose. Autoregressive-moving-average (ARMA models are mathematical models of the persistence, or autocorrelation, in a time series. ARMA models are widely

### at least 50 and preferably 100 observations should be available to build a proper model

III Box-Jenkins Methods 1. Pros and Cons of ARIMA Forecasting a) need for data at least 50 and preferably 100 observations should be available to build a proper model used most frequently for hourly or

### A SOLAR AND WIND INTEGRATED FORECAST TOOL (SWIFT) DESIGNED FOR THE MANAGEMENT OF RENEWABLE ENERGY VARIABILITY ON HAWAIIAN GRID SYSTEMS

ALBANY BARCELONA BANGALORE ICEM 2015 June 26, 2015 Boulder, CO A SOLAR AND WIND INTEGRATED FORECAST TOOL (SWIFT) DESIGNED FOR THE MANAGEMENT OF RENEWABLE ENERGY VARIABILITY ON HAWAIIAN GRID SYSTEMS JOHN

### Spectral Analysis. Al Nosedal University of Toronto. Winter Al Nosedal University of Toronto Spectral Analysis Winter / 71

Spectral Analysis Al Nosedal University of Toronto Winter 2016 Al Nosedal University of Toronto Spectral Analysis Winter 2016 1 / 71 essentially, all models are wrong, but some are useful George E. P.

### Basics: Definitions and Notation. Stationarity. A More Formal Definition

Basics: Definitions and Notation A Univariate is a sequence of measurements of the same variable collected over (usually regular intervals of) time. Usual assumption in many time series techniques is that

Loughborough University Institutional Repository Adaptive diurnal prediction of ambient dry-bulb temperature and solar radiation This item was submitted to Loughborough University's Institutional Repository

### Quantitative Finance I

Quantitative Finance I Linear AR and MA Models (Lecture 4) Winter Semester 01/013 by Lukas Vacha * If viewed in.pdf format - for full functionality use Mathematica 7 (or higher) notebook (.nb) version

### arxiv: v1 [stat.me] 5 Nov 2008

arxiv:0811.0659v1 [stat.me] 5 Nov 2008 Estimation of missing data by using the filtering process in a time series modeling Ahmad Mahir R. and Al-khazaleh A. M. H. School of Mathematical Sciences Faculty

### Chapter 2: Unit Roots

Chapter 2: Unit Roots 1 Contents: Lehrstuhl für Department Empirische of Wirtschaftsforschung Empirical Research and undeconometrics II. Unit Roots... 3 II.1 Integration Level... 3 II.2 Nonstationarity

### A time series is called strictly stationary if the joint distribution of every collection (Y t

5 Time series A time series is a set of observations recorded over time. You can think for example at the GDP of a country over the years (or quarters) or the hourly measurements of temperature over a

### Lesson 13: Box-Jenkins Modeling Strategy for building ARMA models

Lesson 13: Box-Jenkins Modeling Strategy for building ARMA models Facoltà di Economia Università dell Aquila umberto.triacca@gmail.com Introduction In this lesson we present a method to construct an ARMA(p,

### Seasonal Autoregressive Integrated Moving Average Model for Precipitation Time Series

Journal of Mathematics and Statistics 8 (4): 500-505, 2012 ISSN 1549-3644 2012 doi:10.3844/jmssp.2012.500.505 Published Online 8 (4) 2012 (http://www.thescipub.com/jmss.toc) Seasonal Autoregressive Integrated

### Optimization of the forecasting of wind energy production Focus on the day ahead forecast

Optimization of the forecasting of wind energy production Focus on the day ahead forecast Market Processes, EnBW Transportnetze AG Mirjam Eppinger, Dietmar Graeber, Andreas Semmig Quantitative Methods,

### Explanatory Information Analysis for Day-Ahead Price Forecasting in the Iberian Electricity Market

Energies 2015, 8, 10464-10486; doi:10.3390/en80910464 Article OPEN ACCESS energies ISSN 1996-1073 www.mdpi.com/journal/energies Explanatory Information Analysis for Day-Ahead Price Forecasting in the Iberian

### Lesson 2: Analysis of time series

Lesson 2: Analysis of time series Time series Main aims of time series analysis choosing right model statistical testing forecast driving and optimalisation Problems in analysis of time series time problems

### FORECASTING THE INVENTORY LEVEL OF MAGNETIC CARDS IN TOLLING SYSTEM

FORECASTING THE INVENTORY LEVEL OF MAGNETIC CARDS IN TOLLING SYSTEM Bratislav Lazić a, Nebojša Bojović b, Gordana Radivojević b*, Gorana Šormaz a a University of Belgrade, Mihajlo Pupin Institute, Serbia

### Chapter 3: Regression Methods for Trends

Chapter 3: Regression Methods for Trends Time series exhibiting trends over time have a mean function that is some simple function (not necessarily constant) of time. The example random walk graph from

### Introduction to Signal Processing

to Signal Processing Davide Bacciu Dipartimento di Informatica Università di Pisa bacciu@di.unipi.it Intelligent Systems for Pattern Recognition Signals = Time series Definitions Motivations A sequence

### MODELING INFLATION RATES IN NIGERIA: BOX-JENKINS APPROACH. I. U. Moffat and A. E. David Department of Mathematics & Statistics, University of Uyo, Uyo

Vol.4, No.2, pp.2-27, April 216 MODELING INFLATION RATES IN NIGERIA: BOX-JENKINS APPROACH I. U. Moffat and A. E. David Department of Mathematics & Statistics, University of Uyo, Uyo ABSTRACT: This study

### Purdue University Meteorological Tool (PUMET)

Purdue University Meteorological Tool (PUMET) Date: 10/25/2017 Purdue University Meteorological Tool (PUMET) allows users to download and visualize a variety of global meteorological databases, such as

### Short- term solar forecas/ng with sta/s/cal models and combina/on of models

Short- term solar forecas/ng with sta/s/cal models and combina/on of models Philippe Lauret, Mathieu David PIMENT Laboratory, University of La Réunion L. Mazorra Aguiar University of Las Palmas de Gran

### Current best practice of uncertainty forecast for wind energy

Current best practice of uncertainty forecast for wind energy Dr. Matthias Lange Stochastic Methods for Management and Valuation of Energy Storage in the Future German Energy System 17 March 2016 Overview

### Stat 565. (S)Arima & Forecasting. Charlotte Wickham. stat565.cwick.co.nz. Feb

Stat 565 (S)Arima & Forecasting Feb 2 2016 Charlotte Wickham stat565.cwick.co.nz Today A note from HW #3 Pick up with ARIMA processes Introduction to forecasting HW #3 The sample autocorrelation coefficients

### Homework 4. 1 Data analysis problems

Homework 4 1 Data analysis problems This week we will be analyzing a number of data sets. We are going to build ARIMA models using the steps outlined in class. It is also a good idea to read section 3.8

### Econometrics I: Univariate Time Series Econometrics (1)

Econometrics I: Dipartimento di Economia Politica e Metodi Quantitativi University of Pavia Overview of the Lecture 1 st EViews Session VI: Some Theoretical Premises 2 Overview of the Lecture 1 st EViews

### Design of Time Series Model for Road Accident Fatal Death in Tamilnadu

Volume 109 No. 8 2016, 225-232 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu Design of Time Series Model for Road Accident Fatal Death in Tamilnadu

### The log transformation produces a time series whose variance can be treated as constant over time.

TAT 520 Homework 6 Fall 2017 Note: Problem 5 is mandatory for graduate students and extra credit for undergraduates. 1) The quarterly earnings per share for 1960-1980 are in the object in the TA package.

### 4. MA(2) +drift: y t = µ + ɛ t + θ 1 ɛ t 1 + θ 2 ɛ t 2. Mean: where θ(l) = 1 + θ 1 L + θ 2 L 2. Therefore,

61 4. MA(2) +drift: y t = µ + ɛ t + θ 1 ɛ t 1 + θ 2 ɛ t 2 Mean: y t = µ + θ(l)ɛ t, where θ(l) = 1 + θ 1 L + θ 2 L 2. Therefore, E(y t ) = µ + θ(l)e(ɛ t ) = µ 62 Example: MA(q) Model: y t = ɛ t + θ 1 ɛ

### Improving the accuracy of solar irradiance forecasts based on Numerical Weather Prediction

Improving the accuracy of solar irradiance forecasts based on Numerical Weather Prediction Bibek Joshi, Alistair Bruce Sproul, Jessie Kai Copper, Merlinde Kay Why solar power forecasting? Electricity grid

### A HYBRID MODEL OF SARIMA AND ANFIS FOR MACAU AIR POLLUTION INDEX FORECASTING. Eason, Lei Kin Seng (M-A ) Supervisor: Dr.

A HYBRID MODEL OF SARIMA AND ANFIS FOR MACAU AIR POLLUTION INDEX FORECASTING THESIS DISSERATION By Eason, Lei Kin Seng (M-A7-6560-7) Supervisor: Dr. Wan Feng In Fulfillment of Requirements for the Degree

### SEASONAL MODELING OF HOURLY SOLAR IRRADIATION SERIES

SEASONAL MODELING OF HOURLY SOLAR IRRADIATION SERIES M. PAULESCU 1, N. POP 2, N. STEFU 1, E. PAULESCU 1, R. BOATA 3, D. CALINOIU 2 1 West University of Timisoara, Physics Department, V. Parvan 4, 300223

### Prof. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis

Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation

### LOAD FORECASTING APPLICATIONS for THE ENERGY SECTOR

LOAD FORECASTING APPLICATIONS for THE ENERGY SECTOR Boris Bizjak, Univerza v Mariboru, FERI 26.2.2016 1 A) Short-term load forecasting at industrial plant Ravne: Load forecasting using linear regression,

### 5 Transfer function modelling

MSc Further Time Series Analysis 5 Transfer function modelling 5.1 The model Consider the construction of a model for a time series (Y t ) whose values are influenced by the earlier values of a series

### Exercises - Time series analysis

Descriptive analysis of a time series (1) Estimate the trend of the series of gasoline consumption in Spain using a straight line in the period from 1945 to 1995 and generate forecasts for 24 months. Compare

### AR, MA and ARMA models

AR, MA and AR by Hedibert Lopes P Based on Tsay s Analysis of Financial Time Series (3rd edition) P 1 Stationarity 2 3 4 5 6 7 P 8 9 10 11 Outline P Linear Time Series Analysis and Its Applications For

### Unit root problem, solution of difference equations Simple deterministic model, question of unit root

Unit root problem, solution of difference equations Simple deterministic model, question of unit root (1 φ 1 L)X t = µ, Solution X t φ 1 X t 1 = µ X t = A + Bz t with unknown z and unknown A (clearly X

### The autocorrelation and autocovariance functions - helpful tools in the modelling problem

The autocorrelation and autocovariance functions - helpful tools in the modelling problem J. Nowicka-Zagrajek A. Wy lomańska Institute of Mathematics and Computer Science Wroc law University of Technology,

### FORECASTING OF COTTON PRODUCTION IN INDIA USING ARIMA MODEL

FORECASTING OF COTTON PRODUCTION IN INDIA USING ARIMA MODEL S.Poyyamozhi 1, Dr. A. Kachi Mohideen 2. 1 Assistant Professor and Head, Department of Statistics, Government Arts College (Autonomous), Kumbakonam

### Ross Bettinger, Analytical Consultant, Seattle, WA

ABSTRACT USING PROC ARIMA TO MODEL TRENDS IN US HOME PRICES Ross Bettinger, Analytical Consultant, Seattle, WA We demonstrate the use of the Box-Jenkins time series modeling methodology to analyze US home

### AR(p) + I(d) + MA(q) = ARIMA(p, d, q)

AR(p) + I(d) + MA(q) = ARIMA(p, d, q) Outline 1 4.1: Nonstationarity in the Mean 2 ARIMA Arthur Berg AR(p) + I(d)+ MA(q) = ARIMA(p, d, q) 2/ 19 Deterministic Trend Models Polynomial Trend Consider the

### Frequency Forecasting using Time Series ARIMA model

Frequency Forecasting using Time Series ARIMA model Manish Kumar Tikariha DGM(O) NSPCL Bhilai Abstract In view of stringent regulatory stance and recent tariff guidelines, Deviation Settlement mechanism

TEST #1 STA 4853 March 4, 215 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. There are 31 questions. Circle

### CARLOS F. M. COIMBRA (PI) HUGO T. C. PEDRO (CO-PI)

HIGH-FIDELITY SOLAR POWER FORECASTING SYSTEMS FOR THE 392 MW IVANPAH SOLAR PLANT (CSP) AND THE 250 MW CALIFORNIA VALLEY SOLAR RANCH (PV) PROJECT CEC EPC-14-008 CARLOS F. M. COIMBRA (PI) HUGO T. C. PEDRO

### TMA4285 December 2015 Time series models, solution.

Norwegian University of Science and Technology Department of Mathematical Sciences Page of 5 TMA4285 December 205 Time series models, solution. Problem a) (i) The slow decay of the ACF of z t suggest that

### Final Examination 7/6/2011

The Islamic University of Gaza Faculty of Commerce Department of Economics & Applied Statistics Time Series Analysis - Dr. Samir Safi Spring Semester 211 Final Examination 7/6/211 Name: ID: INSTRUCTIONS: