FORECASTING ENERGY INTENSITY WITH FOURIER RESIDUAL MODIFIED GREY MODEL: AN EMPIRICAL STUDY IN TAIWAN

Similar documents
Intermittent demand forecasting by using Neural Network with simulated data

As metioed earlier, directly forecastig o idividual product demads usually result i a far-off forecast that ot oly impairs the quality of subsequet ma

Grey Correlation Analysis of China's Electricity Imports and Its Influence Factors Hongjing Zhang 1, a, Feng Wang 1, b and Zhenkun Tian 1, c

A Study of the Collaborative Degree between Economic Growth and Carbon Reduction Targets under International Comparative Perspective

A VERHULST MODEL ON TIME SERIES ERROR CORRECTED FOR PORT THROUGHPUT FORECASTING

Presenting A Framework To Study Linkages Among Tqm Practices, Supply Chain Management Practices, And Performance Using Dematel Technique

Study on Coal Consumption Curve Fitting of the Thermal Power Based on Genetic Algorithm

Forecasting foreign tourist arrivals in India using time series models

A New Multivariate Markov Chain Model with Applications to Sales Demand Forecasting

Interval Intuitionistic Trapezoidal Fuzzy Prioritized Aggregating Operators and their Application to Multiple Attribute Decision Making

Electricity consumption forecasting method based on MPSO-BP neural network model Youshan Zhang 1, 2,a, Liangdong Guo2, b,qi Li 3, c and Junhui Li2, d

Research on real time compensation of thermal errors of CNC lathe based on linear regression theory Qiu Yongliang

New Exponential Strengthening Buffer Operators and Numerical Simulation

A Study on Travel Time Value in Railway Transportation Corridors

A New Solution Method for the Finite-Horizon Discrete-Time EOQ Problem

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Signal Processing in Mechatronics

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

Properties and Hypothesis Testing

Orthogonal Gaussian Filters for Signal Processing

Molecular Mechanisms of Gas Diffusion in CO 2 Hydrates

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS

Power Weighted Quantile Regression and Its Application

Excellent Performances of The Third-level Disturbed Chaos in The Cryptography Algorithm and The Spread Spectrum Communication

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

Research Article A Unified Weight Formula for Calculating the Sample Variance from Weighted Successive Differences

Modeling and Estimation of a Bivariate Pareto Distribution using the Principle of Maximum Entropy

Signal Processing. Lecture 02: Discrete Time Signals and Systems. Ahmet Taha Koru, Ph. D. Yildiz Technical University.

Analysis of the No-Load Characteristic of the Moving Coil Linear Compressor

EE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course

Prediction of returns in WEEE reverse logistics based on spatial correlation

Stopping oscillations of a simple harmonic oscillator using an impulse force

Research Article Single-Machine Group Scheduling Problems with Deterioration to Minimize the Sum of Completion Times

The Space Redundant Robotic Manipulator Chaotic Motion Dynamics Control Algorithm

The Mathematical Model and the Simulation Modelling Algoritm of the Multitiered Mechanical System

Direction of Arrival Estimation Method in Underdetermined Condition Zhang Youzhi a, Li Weibo b, Wang Hanli c

Approximating the ruin probability of finite-time surplus process with Adaptive Moving Total Exponential Least Square

Time-Domain Representations of LTI Systems

PAijpam.eu ON TENSOR PRODUCT DECOMPOSITION

Four-dimensional Vector Matrix Determinant and Inverse

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan

Generalization of Two Types of Improper Integrals

Uniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations

A proposed discrete distribution for the statistical modeling of

DG Installation in Distribution System for Minimum Loss

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

DISTRIBUTION LAW Okunev I.V.

Polynomial Functions and Their Graphs

11 Correlation and Regression

Control chart for number of customers in the system of M [X] / M / 1 Queueing system

Newton Homotopy Solution for Nonlinear Equations Using Maple14. Abstract

Mathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution

Weighted Correlation Coefficient with a Trigonometric Function Entropy of Intuitionistic Fuzzy Set in Decision Making

New Particle Swarm Neural Networks Model Based Long Term Electrical Load Forecasting in Slovakia

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

Approximate Confidence Interval for the Reciprocal of a Normal Mean with a Known Coefficient of Variation

THE SEISMIC RESPONSE CHARACTERISTICS OF A NEW STRUCTURAL CONFIGURATION

Research Article Research on Application of Regression Least Squares Support Vector Machine on Performance Prediction of Hydraulic Excavator

AN INTERIM REPORT ON SOFT SYSTEMS EVALUATION

Short Term Load Forecasting Using Artificial Neural Network And Imperialist Competitive Algorithm

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY

ADVANCED SOFTWARE ENGINEERING

Double Stage Shrinkage Estimator of Two Parameters. Generalized Exponential Distribution

Research Article Some E-J Generalized Hausdorff Matrices Not of Type M

Research on the Optimization Model of Energy Efficiency of Industrial Park

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

ECONOMIC OPERATION OF POWER SYSTEMS

Stability Evaluation of Surfactant-Grafted Polyacrylamide Used in Oilfields by Evolution of Ammonia

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Research on the Rework Strategies of RMS Based on Process Reliability

A Cobb - Douglas Function Based Index. for Human Development in Egypt

ELEC1200: A System View of Communications: from Signals to Packets Lecture 3

A study of Forecasting Methods for the Cross-Straits Trade

Senvion SE Franz-Lenz-Str. 1, Osnabrück, Germany

Solution of Differential Equation from the Transform Technique


Bayesian and E- Bayesian Method of Estimation of Parameter of Rayleigh Distribution- A Bayesian Approach under Linex Loss Function

ENGI 4421 Confidence Intervals (Two Samples) Page 12-01

A statistical method to determine sample size to estimate characteristic value of soil parameters

Chapter 12 Correlation

A Block Cipher Using Linear Congruences

NEW IDENTIFICATION AND CONTROL METHODS OF SINE-FUNCTION JULIA SETS

Evapotranspiration Estimation Using Support Vector Machines and Hargreaves-Samani Equation for St. Johns, FL, USA

1 Introduction to reducing variance in Monte Carlo simulations

Estimating the Change Point of Bivariate Binomial Processes Experiencing Step Changes in Their Mean

Stability Analysis of the Euler Discretization for SIR Epidemic Model

Statistics 511 Additional Materials

OBJECTIVES. Chapter 1 INTRODUCTION TO INSTRUMENTATION FUNCTION AND ADVANTAGES INTRODUCTION. At the end of this chapter, students should be able to:

CHAPTER 4 BIVARIATE DISTRIBUTION EXTENSION

Symmetric Division Deg Energy of a Graph

Noah Williams Economics 312. University of Wisconsin Spring Midterm Examination Solutions 1 FOR GRADUATE STUDENTS ONLY

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Chapter 9: Numerical Differentiation

Chimica Inorganica 3

Progress In Electromagnetics Research, PIER 51, , 2005

Grey forecast model for torpedo development cost based on neural network

Some New Iterative Methods for Solving Nonlinear Equations

Recurrence Relations

Product Mix Problem with Radom Return and Preference of Production Quantity. Osaka University Japan

Summary: CORRELATION & LINEAR REGRESSION. GC. Students are advised to refer to lecture notes for the GC operations to obtain scatter diagram.

Transcription:

ORECASTING ENERGY INTENSITY WITH OURIER RESIDUAL MODIIED GREY MODEL: AN EMPIRICAL STUDY IN TAIWAN Thah-Lam Nguye, Yig-ag Huag Natioal Kaohsiug Uiversity of Applied Scieces, Kaohsiug 80778, Taiwa ABSTRACT: Eergy itesity is defied as the eergy cosumptio for producig every uit of real GDP i a certai time frame. Studies i forecastig the eergy itesity have ot well positioed due to the difficulty i collectig relevat data o the determiats affectig the eergy cosumptio ad GDP. Therefore, i this study, it is proposed to use the Grey forecastig model GM(,) to predict the eergy cosumptio ad real GDP before the itesity is forecasted. To ehace the accuracy level of the forecastig models, their residuals are the modified with ourier series. I the case of Taiwa, the modified models resulted i very low values of mea of absolute percetage error (MAPE) of 0.33% ad 0.58%, respectively to the eergy cosumptio ad real GDP. Hece, the modified model is strogly suggested to forecast the eergy itesity i Taiwa from 0-05. Keywords: GM(,), GM(,), Grey forecastig, ourier series, Eergy itesity. I. INTRODUCTION Eergy is the core of most ecoomic, evirometal ad developmetal issues aroud the globe. It has bee well proved that there is a close relatioship betwee the eergy cosumptio ad ecoomic developmet. As per the defiitio offered by the Departmet of Ecoomic ad Social Affairs of the Uited Natios Secretariat, eergy itesity is defied as the eergy cosumptio for producig every uit of real GDP i a certai time frame which meas that the lower the eergy itesity of a ecoomy is, the better the ecoomy performs. Eergy itesity idicates the total eergy used to support a wide rage of productio ad cosumptio (ecoomic ad social) activities []. Therefore, it is usually cosidered as oe of the measures of sustaiable developmet. A coutry with highly ecoomical productivity, pleasat weather, geographically well-allocated work places, fuel efficiet vehicles, mass trasportatio, etc., will have a far lower eergy itesity ad vice versa. May researches have bee coducted ad it has bee foud out that eergy cosumptio ad GDP are positively correlated though the correlatio coefficiets may be differet from coutry to coutry []. Chages i the ecoomy structure may result i ug less additioal eergy; however, total eergy cosumptio is still icreag [3]. I reducig the emissios through reducig the eergy cosumptio, it was poited out that the developed coutries ted to be more affected by such policy rather tha developig oes [4]. May differet researchers [5-5] have focused o the aalysis of the relatioship betwee the eergy cosumptio ad GDP. Nevertheless, the umber of studies i forecastig the eergy itesity is actually limited due to the fact that collectig relevat data o the determiats affectig the eergy cosumptio ad GDP rus ito a lot of difficulties. Therefore, i this study, it is proposed to use the covetioal Grey forecastig model GM(,), which has bee widely used i differet areas due to its ability to deal with the problems of ucertaity with few data poits ad/or partial kow, partial ukow iformatio, to predict the eergy cosumptio ad real GDP before the itesity is forecasted. To ehace the accuracy level of the forecastig model, its residuals are the modified with ourier series. A empirical study i Taiwa is ivestigated as a example for this improved model.

II. LITERATURE REVIEWS. Grey Model Grey theory offers a ew approach to deal maily with the problems of ucertaity with few data poits ad/or poor iformatio which is said to be partial kow, partial ukow [6]. The core of Grey theory is the Grey dyamic model which is usually called Grey model (GM). The Grey model is used to execute the short-term forecastig operatio with o strict hypothesis for the distributio of the origial data series [7]. The geeral GM model has the form of GM(d,v), d is the rak of differetial equatio ad v is the umber of variables appeared i the equatio. The basic model of Grey model is GM(,), a first-order differetial model with oe iput variable. The procedure to obtai GM(,) is as the followig: Step : Suppose a origial series with (0) etries is x : (0) (0) (0) (0) x x,, x ( k), x ( ) x (0) ( k ) is the value at time k k,. Step : rom the origial series x (0), a ew series x ca be geerated by oe time accumulated geeratig operatio (- AGO), which is x x,, x ( k),, x ( ) k (0) j x ( k) x ( j) Step 3: A first-order differetial equatio with oe variable is expressed as: ax b 3 dx dt a is called a developig coefficiet ad b is called a grey iput coefficiet. These two coefficiets ca be determied by the least square method as show below:, T T T a b B B B Y 4 x x () / x () x (3) / B x ( ) x ( ) / (0) (0) (0) Y x (), x (3),, x ( ) T Therefore, the solutio of equatio (3) is expressed as: b at b x t x e a 5 a Equatio (5) is also kow as time respose fuctio of the equatio (3). rom equatio (5), the time respose fuctio of the GM(,) is give by: (0) b a k ˆ b x k x e a a 6 k, Based o the operatio of oe time iverse accumulated geeratig operatio (0) (-IAGO), the predicted series ˆx ca be obtaied as the followig: (0) (0) (0) (0) xˆ xˆ,, xˆ k,, xˆ 7 (0) xˆ xˆ (0) xˆ k xˆ k xˆ k k,. ourier Residual Modificatio I order to improve the accuracy of forecastig models, the ourier series has bee widely ad successfully applied i modifyig the residuals i Grey forecastig model GM(,) which reduces the values of RMSE, MAE, MAPE, etc., [8-]. The overall procedure to obtai the modified model is as the followigs: Let x is the orgial series of etries ad ˆx is the predicted series obtaied from GM(,). Based o the predicted series ˆx, a residual series amed is defied as:, 3, 4,, k,, 8 ( k) x( k) xˆ ( k) k, Expressed i ourier series, k is rewritte as: k a0 ai D bi D 9 i D ik / ( ) k,

/ called the rewritte as: miimum deploymet frequecy of ourier PC. 0 series [] ad oly take iteger umber [8-0]. Ad therefore, the residual series is 3 3 3 P 3 C a, a, b, a, b,, a, b 0 The parameters a0, a, b, a, b,, a, b are obtaied by ug the ordiary least squares method (OLS) which results i the equatio of: T T T C P P P Oce the parameters are calculated, the modified residual series ˆ is the achieved based o the followig expressio: ˆ k a0 ai D bi D i rom the predicted series ˆx ad ˆ, the ourier modified series x of series ˆx is determied by: x x, x,, x k,, x x xˆ x k xˆ k ˆ k k, To evaluate the model accuracy, there are four importat idexes to be cosidered, such as: The mea absolute percetage error (MAPE) [9,, 3]: x( k) v( k) MAPE k, x( k) k v(k) is the forecasted value of kth etry from the model ( v( k) xˆ ( k) i T GM(,) or v( k) x( k) i GM(,)). The post-error ratio C [4, 5]: S C S S x ( k ) x x x ( k ) S k k k ( k) ( k) x( k) v( k) ( k) k The smaller the C value is, the higher accuracy the model has. The small error probability P [4, 5]: k P p 0.6745 S The higher the P value is, the higher accuracy the model has. The forecastig accuracy [5]: MAPE The above four idexes are used to classify the grades of forecastig accuracy as i Table. Table. our grades of forecastig accuracy Grade level MAPE C P I (Excellet) < 0.0 < 0.35 > 0.95 > 0.95 II (Good) < 0.05 < 0.50 > 0.80 > 0.90 III (Qualified) < 0.0 < 0.65 > 0.70 > 0.85 IV (Uqualified) 0.0 0.65 0.70 0.85

III. EMPIRICAL RESULTS The data of eergy cosumptio from 999 0 i Taiwa are obtaied from the Bureau of Eergy of Miistry of Ecoomic Affairs of Taiwa [6]; as the data of the Taiwa GDP from the same period are collected from Iteratioal Moetary ud [7]. Oly data from 999-00 are used to build relevat GM(,) ad GM(,) models. The data i 0 is used to compare with the forecasted value from the selected model to further affirm its forecastig power. 3.. orecastig model for the eergy cosumptio Based o the algorithm expressed i sectio., the fudametal Grey forecastig model for the eergy cosumptio amed GM(,) E is foud as the followig: 0.03073( k) xˆ ( k) 3377369.76 e 385397.6 The residual series attaied from GM(,) E is the modified with ourier series, which results i the modified model GM(,) E as per the algorithm stated i sectio.. The evaluatio idexes of GM(,) E ad GM(,) E are summarized as i Table. Table clearly showed that betwee GM(,) E ad GM(,) E, GM(,) E is selected because it has a lower value of MAPE ad a better forecastig power. So, GM(,) E is used to forecast the eergy cosumptio i 0. The forecasted value is the compared with the actual cosumptio i order to further affirm its forecastig power as show i Table 3. The MAPE value of 5.55% idicates that GM(,) E ca be appropriately used to forecast the cosumptio i 0 05. The forecasted values i this period are show i Table 4. 3.. orecastig model for the GDP Similarly, the fudametal Grey forecastig model for the GDP amed GM(,) G is foud as the followig: 0.0309( ) ˆ ( ) 8943.0 k x k e 8667.89 GM(,) G is accordigly obtaied based o sectio.. It is also selected because it outperforms GM(,) G i term of low MAPE value as show i Table. Its forecasted value of GDP i 0 show i Table 3 has a MAPE value of.83% idicatig that it ca be used to forecast the GDP i 0 05. Its relevat forecasted values are also show i Table 4. Table. Summary of evaluatio idexes of model accuracy Idex orecastig MAPE S S C P Model power GM(,) E 0.035 478.86 503.84 0.34.00 0.9648 Good GM(,) E 0.0033 478.86 475.76 0.03.00 0.9967 Excellet GM(,) G 0.044 4. 5.94 0.38.00 0.9586 Good GM(,) G 0.0058 4..3 0.05.00 0.994 Excellet Table 3. orecasted eergy cosumptio ad GDP i 0 Model Uit Actual value orecasted value MAPE GM(,) E 0 3 KLOE 3,83.50 3948.40 0.0555 GM(,) G 0 9 USD 430.58 4.7 0.083 Table 4. orecasted eergy cosumptio ad GDP from 0-05 Item Uit 0 03 04 05 Eergy 0 3 KLOE 4859.70 4759.90 5744.00 57530.00 GDP 0 9 USD 447.45 454.54 465.6 479.98 Eergy itesity KLOE/0 6 USD 33.9 34.7 38.30 38.0

Table 4 shows that there is a small decrease i the eergy itesity idex of Taiwa i the comig years. This could be explaied as a outcome of the past ad curret ivestmet i the ew productio techology as well as the moder facilities i the trasportatio, services ad residetial sectors. Besides, the movig of its maufacturig factories to other coutries icludig Chia, Vietam, Idoesia, Malaysia, Laos, Thailad, etc., as well as the ehacig of its service idustries make Taiwa ot oly cosume less eergy but also produce higher GDP, which sigificatly cotribute to the decrease of the eergy itesity idex of Taiwa. IV. CONCLUSION The accuracy level of the traditioal Grey forecastig model GM(,) ca be well improved if the model is modified with ourier series. I the case of eergy itesity of Taiwa, with the ourier modified Grey forecastig model GM(,), it was foud out that the eergy itesity of Taiwa becomes lower ad lower represetig a better & stable developmet of the coutry. This result plays as a excellet motivatio for the authorities to assert that they are o the right way to develop Taiwa i geeral ad its ecoomy i particular. Other coutries could refer to this as a good example to focus o research & developmet as well as ivest ad apply advaced techology i most of their activities. [] Departmet of Ecoomic ad Social Affairs of the Uited Natios Secretariat, Idicators of sustaiable developmet: Guidelies ad Methodologies, 3 rd Editio, pp.87, Uited Natios, New York, 007. [] R. Haesso, Eergy ad GDP growth, Iteratioal Joural of Eergy Sector Maagemet, Vol. 3, Issue, 009, pp.57 70. [3] Europea Evirometal Agecy, Total primary eergy itesity (CSI 08/ENER 07), EEA, 006. [4] J. Chotaawat, L. C. Hut ad R. Pierse, Casuality betwee eergy cosumptio ad GDP: Evidece from 30 OECD ad 78 No-OECD coutries, Surrey Eergy Ecoomics Cetre (SEEC), Uiversity of Surrey, Uited Kigdom, 006, pp.4. [5] D. I. Ster, A multivariate Coitegratio Aalysis of the role of the eergy i US, Eergy Ecoomics, Vol., 000, pp.67 83. [6]. Climet ad A. Pardo, Decouplig factors o the eergy-output likage: The Spaish case, Eergy Policy, Vol. 35, 007, pp.5 58. REERENCES [7] K. Ghali ad M. El-Sakka, Eergy use ad output growth i Caada: A multivariate Coitegratio aalysis, Eergy Ecoomics, Vol. 6, 004, pp.5 38. [8] U. Soytas ad R. Sari, The relatioship betwee eergy ad productio: Evidece from Turkish maufacturig idustry, Eergy Ecoomics, Vol. 9, pp.5 65. [9] B. S. Warr ad R. U. Ayres, Evidece of causality betwee the quatity ad quality of eergy cosumptio ad ecoomic growth, Eergy, Vol. 35, 00, pp.688 693. [0] C. Lee ad C. Chag, Structural breaks, eergy cosumptio, ad ecoomic growth revisited: Evidece from Taiwa, Eergy Ecoomics, Vol. 7, 005, pp.857 87. [] J. E. Paye, US disaggregate fossil fuel cosumptio ad real GDP: A empirical ote, Eergy Sources, Part B: Ecoomics, Plaig ad Policy, Vol. 6, 0, pp.63 68. [] J. Asafu-Adjaye, The relatioship betwee eergy cosumptio, eergy prices ad ecoomic growth: Time

series evidece from Asia developig coutries, Eergy Ecoomics, Vol., 000, pp.65 65. [3] K. atai, L. Oxley ad. G. Scrimgeour, Modellig the causal relatioship betwee eergy cosumptio ad GDP i New Zealad, Australia, Idia, Idoesia, Philipies, ad Thailad, Mathematics ad Computers i Simulatio, Vol. 64, 004, pp.43 445. [4]. Climet ad A. Pardo, Decouplig factors o the eergy-output likage: The Spaish case, Eergy Policy, Vol. 35, 007, pp.5 58. [5] J. L. Hu ad C. H. Li, Disaggregated eergy cosumptio ad GDP i Taiwa: A threshold co-itegratio aalysis, Eergy Ecoomics, Vol. 30, 008, pp.34 358. [6] S.. Liu, Advaces i Grey system theory ad its applicatios, Proceedigs of 007 IEEE Iteratioal Coferece o Grey Systems ad Itelliget Services, Chia, 007, pp. 6. [7] S. J. Wag, W. L. Wag, C. T. Huag ad S. C. Che, Improvig ivetory effectiveess i RID-eabled global supply chai with Grey forecastig model, Joural of Strategic Iformatio Systems, Vol. 0, 0, pp.307 3. [8] M. Askari ad A. etaat, Log-term load forecastig i power system: Grey system predictio-based models, Joural of Applied Scieces, Vol., 0, pp.3034 3038. [9] Z. Guo, X. Sog ad J. Ye, A Verhulst model o time series error corrected for port throughput forecastig, Joural of the Easter Asia Society for Trasportatio Studies, Vol. 6, 005, pp.88 89. [0] L.C. Hsu, Applyig the grey predictio model to the global itegrated circuit idustry, Techological orecastig ad Social Chage, Vol. 70, 003, pp.563 574. [] Y. L. Huag ad Y. H. Lee, Accurately forecastig model for the Stochastic Volatility data i tourism demad, Moder Ecoomy, Vol., 0, pp.83 89. [] M.L. Ka, Y. B. Lee ad W. C. Che, Apply grey predictio i the umber of Tourist, The 4 th iteratioal coferece o Geetic ad Evolutioary Computig, 00, pp.48 484. [3] C. C. Hsu ad C. Y. Che, Applicatios of improved grey predictio model for power demad forecastig, Eergy Coversio ad Maagemet, Vol. 44, 003, pp.4 49. [4] J. Hua ad Q. Liag, Applicatio research of the Grey forecast i the Logistics demad forecast, irst iteratioal workshop o Educatio Techology ad Computer Sciece, 009, pp. 36 363. [5] H. Ma ad Z. Zhag, Grey predictio with Markov-Chai for Crude oil productio ad cosumptio i Chia, Advaces i Itelliget ad Soft Computig, Vol. 56, 009, pp.55 56. [6] Bureau of Eergy of Miistry of Ecoomic Affairs of Taiwa, Eergy Statistics. Retrieved o April 5, 0 from www.moeaboe.gov.tw/eglish/ Statistics/EStatistics.aspx [7] Iteratioal Moetary ud, Coutry iformatio. Retrieved o April 0, 0 from www.tradigecoomics.com/ taiwa/gdp Correspodig author: Thah-Lam Nguye Graduate Istitute of Mechaical ad Precisio Egieerig, Natioal Kaohsiug Uiversity of Applied Scieces 45, Chie Kug Rd., Kaohsiug 80778, Taiwa, R.O.C. Email: gree4rest.v@gmail.com