Aerospace Science and Technology
|
|
- Edmund Cooper
- 5 years ago
- Views:
Transcription
1 Aerospace Science and Technology 16 (2012) Contents lists available at ScienceDirect Aerospace Science and Technology Aircraft engine health management via stochastic modelling of flight data interrelations D. Dimogianopoulos, J. Hios, S. Fassois Stochastic Mechanical Systems & Automation (SMSA) Laboratory, Department of Mechanical Engineering & Aeronautics, University of Patras, GR Patras, Greece article info abstract Article history: Received 20 November 2009 Received in revised form 1 February 2011 Accepted 2 March 2011 Availableonline5March2011 Keywords: Health management Fault detection Fault isolation Stochastic modelling Nonlinear models Engine systems Decision maing Statistical inference A novel engine health management (EHM) scheme is introduced. It uses flight-level, instead of thermodynamic, data to cost-effectively augment the onboard EHM redundancy. For a nominal healthy aircraft, fault-sensitive interrelations among flight data are globally modelled inside a flight regime via Constant-Coefficient Pooled Nonlinear AutoRegressive with exogenous (CCP-NARX) excitation representations. Single or sequential engine faults perturb these interrelations. Statistically evaluating the perturbation-induced effects draws reliable conclusions on the engine s health. Validation and comparisons with Kalman filter-based alternatives are made throughout the regime under various operational conditions Elsevier Masson SAS. All rights reserved. 1. Introduction In modern aircraft, and even more so in future pilot-less versions, high reliability and safety should be cost-effectively obtained. The multiplication of existing critical hardware (the hardware redundancy principle), coupled to voting schemes for health management (HM) purposes, has its limits due to added weight and cost [26]. Significant improvements on reliability/safety may, hence, only result from additional analytical redundancy based on the smart use of available system-level signals and hardware [18]. Introducing engine HM (EHM) analytical redundancy by modelling part of the healthy engine dynamics using physics-based principles, and then detecting fault-induced trends in model parameters (or specific functions of them) is commonly used in either on-board [6] or on-the-service-bay [34] versions. Other schemes use the gas path analysis (GPA) approach to identify highly accurate models on typical health-state-related data along the gas path (afterfan total temperature and pressure, high pressure compressor temperature and pressure, and so on): In [25] linear stochastic models (AutoRegressive with exogenous excitation * Corresponding author. Tel.: ; fax: Since March 2010, the author is Assistant Professor at the Technological Educational Institute of Piraeus, 250 Thivon Av., Egaleo Athens, Greece. Tel: addresses: dimogian@mech.upatras.gr (D. Dimogianopoulos), hiosj@mech.upatras.gr (J. Hios), fassois@mech.upatras.gr (S. Fassois). ARX, Output Error OE or Prediction Error PE) are augmented with model-error models, designed following the H principle and aiming at minimizing modelling uncertainties around an envelope location. In [27] the model augmentation uses a neural networ-based part. In both schemes fault-induced trends between the model output and that of the actual engine provides HM results. A majority of approaches rely, however, on models of the healthy engine dynamics identified via GPA measurements and coupled to fuzzy logic rules [19], or Kalman filters (KF) in standard or adaptive form [24] for estimating fault-induced changes on health parameters. An interesting comprehensive study [30] compares schemes using KFs in linear (LKF), extended (EKF) and unscented (UKF) versions. Both EKF- and UKF-based schemes show superior performance at the price of (one and two, respectively, orders of magnitude) higher online computational effort than the LKF. This is due to necessary onboard transformations (which potentially involve badly-conditioned matrices). Other EHM schemes rely on multiple engine models or observers to describe the engine dynamics in various health states, and compare the output of each model to that of the engine. Then, the correct model and the associated health state are deduced based on the model output being close to that currently measured from the engine. Multiple health-state-related models are built using either KFs [21,23], Taagi Sugeno fuzzy rules [12], or neural networs [11,28]. Multiple health-state-related observers (designed /$ see front matter 2011 Elsevier Masson SAS. All rights reserved. doi: /j.ast
2 D. Dimogianopoulos et al. / Aerospace Science and Technology 16 (2012) Fig. 1. Theenginehealthmanagementschemeanddetailoftheaircraftsimulatorwiththeenginesubsystem. for enhanced robustness) are reported in [10,29], or, in nonlinear versions in [35]. Finally, non-model based schemes use thresholds on engine health state parameters [32], pattern recognition techniques [3], or even fuzzy inference [2] on groups of engine parameters to detect fault-induced trends. These wors are comprehensively reviewed in [33]. The proposed onboard EHM scheme innovates by operating on flight-level, rather than engine thermodynamic, data. Since available data are used, the expensive/complex development of new physical sensors or hardware is avoided. The resulting (unrelated to GPA) EHM redundancy is intended to complement the existing GPA-related one, in order to enhance the available decision-maing capabilities. The approach extends the generic ideas of relevant wor in [13], and relies on stochastic nonlinear modelling of the interrelations among flight data such as acceleration, thrust and so on. Hence, the use of (often unavailable) physics-based or empirical models, as in [6], is no longer needed. These interrelations are fault-sensitive, that is, valid exclusively for a nominal healthy aircraft, and are modelled via Constant-Coefficient Pooled Nonlinear AutoRegressive with exogenous (CCP-NARX) excitation representations. Potential single or sequentially occurring engine faults perturb these nominal (CCP-NARX modelled) interrelations. Then, custom-built statistical hypothesis tests evaluate the perturbationinduced effects in order to detect abnormal engine operation, and (when possible) isolate causes of engine malfunctions. Hence, additional ability to judge the GPA-based internal engine diagnosis is provided. The paper is organized as follows: Section 2 presents the aircraft (a nonlinear simulator) and the considered engine faults. Section 3 describes the CCP-NARX modelling procedure and the design of the statistical tests. Section 4 shows some modelling and HM results, along with brief comparisons with LKF-based schemes. Finally, Section 5 presents some concluding remars. 2. The aircraft and the faults A simulator of an autonomous twin-engined small passenger aircraft is used in this study. It is a 6 degree-of-freedom nonlinear aircraft model, with the inputs (stic, wheel, pedal, throttle) provided from a specifically designed autopilot bloc according to a predetermined flight trajectory (Fig. 1). The input and the resulting attitude signals (including the angle-of-attac α, the sideslip angle β, surface moments, accelerations and angular rates) are available. The wind and turbulence effects (conforming to MIL- F-8785C) are considered as system disturbances. The considered faults affect the left engine, and are simulated by acting on selected signals fed into the engine subsystem, as shown in Fig. 1. Table 1 The faults considered. Type Description Magnitude F A Part of throttle input (% oftotal) = 25 fed to engine F B Unstable thrust & flameout = 550 F C Recurrent engine thrust reduction (% ofnominal) = 30 F D F A fault shortly followed (in t s) 1,2,t 1 1 = 25, 2 = 550, by F B fault t = 10, 20 s 2 The HM scheme monitors the engine operational level, rather than any internal components (fans, compressors and so on) as in GPAbased schemes. Thus, a relatively simple engine model is considered, since no GPA signals (other than the total thrust) have to be used. The simulated faults involve: a) Partial reduction of the throttle input signal fed into the left engine control module, with the measurements of altitude and Mach signals (designated as Alt[t] and Mach[t], respectively) being unaffected. The engine operation is substantially altered, as if the lin from throttle to engine suffered a severe, but not fatal, failure (see Fig. 2(a) from instant 150 s onwards). This fault is referred to as F A with A indicating the fault type and its considered magnitude (% of the throttle input normally entering the engine control module, see Table 1). b) Unstable left engine thrust output simulating inconsistent compressor airflow before flameout (see Fig. 2(b)), that is, engine shutdown. This situation may be due to severe successive internal incidents (such as blade loss) causing fast deterioration. This fault is referred to as F B with B indicating the fault type, and its magnitude, corresponding to the colored noise s standard deviation σ = added to the thrust output to simulate the fault. c) Recurrent left engine thrust reduction (see Fig. 2(c)). The left engine thrust repeatedly drops to lower values for short time intervals, possibly due to unstable compressor airflow. This situation may be due to either physical damage, or control software issues. The fault is referred to as F C,withC indicating the fault type and the % periodic reduction of the current left engine thrust nominal value. d) Multiple sequentially occurring faults affecting the left engine. The occurrence of F A faults is followed some seconds later by F B faults affecting the already faulty engine. Clearly, this sequential occurrence is not equivalent to superposing the effects of single F A and single F B faults, meaning that a separate F D class is necessary. These sequential faults are 1, 2,t referred to as F D, with D indicating the fault type and 1, 2,t
3 72 D. Dimogianopoulos et al. / Aerospace Science and Technology 16 (2012) Fig. 2. Total engine thrust response for faults occurring at 150 s: (a) Healthy versus F A 25 ;(b)healthyversusf B 550 ;(c)healthyversusf C 30 ;(d)healthyversusf D 25,550,10 affected aircraft. 1, 2 the magnitudes of the initial F A and the subsequent F B 1 2 faults, respectively. Finally, t indicates the time interval in s between the F A 1 and F B 2 fault occurrence (see a typical case in Fig. 2(d)). Both F A and F B are considered as modelled faults, meaning that their possible occurrence has been accounted for during the HM scheme s design. The same is valid for F D faults, since they 1, 2,t consist of single modelled but sequentially occurring ones. On the contrary, F C faults are considered as unmodelled, since many possible incidents cause a recurrent partial loss of engine thrust. As shown in Section 4.2, an abnormal engine operation engine may always be detected, but the originating fault may be isolated only if it belongs to the set of modelled faults. Remar 1. Shortly after F B or F C fault occurrence at 150 s, the affected engine s thrust in Figs. 2(b) and (c), respectively, is largely similar to that from a healthy engine. Changes are noted only after the flameout (160 s onwards) in the F B case. This is due to the choice of F B and F C magnitude, intended to produce minor fault-induced effects on the total thrust, in order to chec the EHM scheme s sensitivity to fault occurrence. 3. The engine health management scheme 3.1. Overview of the scheme s design and operation The attitude of a standard passenger aircraft is designed to be linear under the (auto) pilot s action, and may, hence, be modelled via linear representations. Nevertheless, abrupt maneuvering, significant turbulence/wind and other disturbances may induce nonlinearities in relationships between specific flight variables (see Section 4.1). Such factors are also responsible for extra noise in measured signals. Finally, aircraft operation in broad envelope areas may also induce modelling uncertainty, meaning that models with varying parameters are required to accurately represent the dynamics at each area s point. Instead of using multiple, often adaptive, local models (associated to high onboard computational burden) for globally describing the aircraft dynamics as in [21], the proposed scheme relies on a single stochastic pooled nonlinear model of the considered dynamics (see Fig. 1). Specifically, a Multi-Input Single-Output (MISO) CCP-NARX representation (Section 3.2) describes the relationship among four signals obtained from healthy aircraft; three taen as inputs (lateral acceleration A y, side-slip angle β and rudder moment) and one as output (total engine thrust). Thrust from one engine may also be used, for compatibility with vectoring strategies. Note that, instead of examining separate signals, the joint dependence between all inputs and the thrust is considered. The reason is that, for instance, changes in roll rate may be caused by either F A faults or manoeuvering, faults in control surfaces and so on. Such changes must be related to the engine behavior, in order to obtain reliable EHM results. This particular input and output signal selection is made for two reasons: first, because these input signals are random functions of time due to turbulence and other factors and have, thus, varying amplitude and rich frequency content even under reduced pilot activity. This provides persistent excitation, which is the prerequisite for effective CCP-NARX identification [20]. Second, because several fault-induced engine thrust variations may be associated to specific combined changes of the acceleration A y, β angle and rudder moment. Naturally, other signal choices are possible, provided that the previous conditions are met. The current MISO CCP-NARX model addresses all previous modelling issues, because: The model s stochastic nonlinear form explicitly accounts both for noise-in-measurement issues and for nonlinearities in the dynamics, without requiring nown disturbance bounds as when linear [9] or nonlinear deterministic models are used [35,36]. Via the pooling technique (Section 3.2) relevant recorded past data may contribute to the CCP-NARX identification (see Fig. 3). Thus, modelling uncertainties from operation in broad envelope areas are globally reduced, instead of being minimized around given points via H or other robust approaches [25]. The identified CCP-NARX model is then used during the operational (diagnostic) phase (Section 3.3), which taes place onboard an aircraft operating with engines in unnown health state. Using data collected during any such flight and the identified CCP-NARX
4 D. Dimogianopoulos et al. / Aerospace Science and Technology 16 (2012) Fig. 3. Schematic of the pooled representation s identification problem, showing the flight envelope, the excitation/response signals from different flights and the CCP-NARX structure. model, the engine s current condition is assessed as follows: For healthy engines, the obtained residual sequence e[t] (that is, the difference between the real thrust output and its CCP-NARX resulting prediction, see Fig. 1) features certain statistical properties. For faulty engines, these statistical properties are no longer valid. Then, statistical decision-maing hypothesis tests may chec and evaluate the fault-related information in e[t] for EHM purposes. Remar 2. The proposed EHM scheme is designed to operate along (and interact with) similar HM schemes, monitoring critical systems such as control surfaces (see [15]). Hence, the analytical decision-maing ability is enhanced and further complements that based on GPA or other approaches for cost-effectively promoting safety. Finally, the interaction between HM schemes enables crosschecing of HM decisions, so that faults affecting a given circuit may not influence HM decisions in other monitored systems Baseline modelling phase (on-the-ground tuning phase) The CCP-NARX model representing the interrelations among the considered input/output signals admits the form: L y j [t]= θ i p i, j [t]+e j [t] { i=0 E e j [t] e s [t τ ] } = γ e [ j, s]δ[τ ] e j [t] ( NID 0,σe 2 ( j)) (1) with t designating normalized discrete time, y j [t] and e j [t] the representation s output and prediction error (or residual, assumed to be a zero-mean uncorrelated sequence with variance σe 2( j)) signals of the j-th flight, respectively. The input signals u m, j [t] (m = 1, 2,...) are not explicitly shown in (1), but are involved in the regressors p i, j [t] (i = 0,...,L and j = 1,...,M). These are nonlinear terms such as products between lagged values of the output y j [t] and the m-th input u m, j [t], orevenpowersofthese signal values (with p 0, j [t] = 1, by definition). The nonlinearity degree of the representation, that is, the sum of powers of the signal values involved in each regressor is less than or equal to nl. The coefficient corresponding to the i-th regressor is designated as θ i and is constant for all flights. The maximum lags of signals y j [t], u m, j [t] in (1) (model orders) are n y and n um, respectively. The statistical expectation is designated as E{ }, whereas NID(.,.) stands for Normally Independently Distributed (with the indicated mean and variance), δ[τ ] the Kronecer delta (δ[τ ]=1 when τ = 0 and δ[τ ]=0 when τ 0), and γ e [ j, s] the cross covariance (see [8, pp ]). Note that, from (1), the residuals from two different flights j and s may potentially be correlated. Rewriting (1) in a linear regression form yields 1 : y j [t]=φ T j [t] θ + e j[t] (2) with φ j [t] =[p 0, j [t]...p L, j [t]] T the regressor and θ =[θ 0...θ L ] T the vector of coefficients, respectively. If M flights (each of N recorded data samples) are available for the modelling phase, the following matrix equation is obtained: y = Φ θ + e (3) where y = [y 1... y M ] T R [NM 1], Φ = [Φ 1...Φ M ] T R [NM (L+1)], e =[e 1...e M ] T R [NM 1]. The vectors y j, e j R [N 1] include the N data samples from the j-th flight, while Φ j =[φ T j [1]...φT j [N]]T R N (L+1). Finally, note that when a single flight is considered, the previous CCP-NARX representation becomes a simple NARX one. During the baseline phase, the objective is to identify the CCP- NARX(n y, n um ) model, that is: a) To choose the regressors p i, j of nonlinearity degree up to nl, which most accurately describe the interrelations between the considered signals (structure selection), and b) to compute the associated vector of CCP-NARX coefficients. The structure selection involves determination of nl for the mapping between the I/O signals, based on a ey property of linear systems: Their response to a sinusoidal input (of arbitrary amplitude but given frequency) is also a sinusoidal signal of the same 1 Lower case/capital bold face symbols designate column vector/matrix quantities, respectively.
5 74 D. Dimogianopoulos et al. / Aerospace Science and Technology 16 (2012) frequency (and shifted phase). The ratio of the output over the input amplitudes is constant and only related to the frequency used. On the contrary, for nonlinear systems this ratio may change for a given input frequency according to the input amplitude. Hence, determining nl involves exciting the nonlinear system with input signals of different amplitudes at a given frequency, and creating a plot of the output amplitude versus that of the input: The value of nl equals the degree of the polynomial, which approximates sufficiently well the plot curve. The procedure is repeated for as many frequencies as possible inside the frequency range of interest, for finding the most representative nl with respect to the considered frequencies. The determination of the current system s nonlinearity degree is presented in Section 4.1. For a given nl value, the regressors are selected using either the Forward Orthogonal Search (FOS) algorithm [7], or bootstrapping techniques [22], or even genetic algorithms. Past experience with CCP-NARX identification [13] shows that the best results are achieved by combining the FOS algorithm with the bootstrapping technique. In the next paragraphs, this procedure is briefly outlined, with more details (along with presentations of the FOS algorithm and the bootstrapping technique) given in [13]. i) From the original data set of M flights, each data sequence (of N samples) is sectioned to a number of N s -sample long overlapping segments. ii) New pooled data sets, each one containing M N s samples, are formed and for each such set the FOS algorithm performs structure selection. The regressors identified for each such set represent the dynamics of M flights during the short time segment of N s samples. Due to the overlap, the selected regressors account effectively for changing dynamics between two successive segments. iii) The union of all previously selected regressor forms the Extended Regressor Set (ERS), which is highly representative of the dynamics of all M flights throughout the considered time period. Via the bootstrapping technique [22] the ERS is considerably reduced without noticeably deteriorating the modelling accuracy. iv) Finally, the number of regressors obtained in step iii) may be reduced again via term clustering techniques [1]. Each cluster contains regressors formed by the same combination of signals (but of different lagged values). The clusters are successively removed from the CCP-NARX model (one at each iteration) and the accuracy of the remaining CCP-NARX model part is assessed by computing the Residual e Sum of Squares over the Signal y Sum of Squares (RSS/SSS) ratio, with signals e and y found in (3). The omitted clusters causing significant RSS/SSS increases are retained, thus yielding a compact final set of regressors. The corresponding CCP-NARX model is validated via the tests in [7]: The successful CCP-NARX model must offer the best compromise between low RSS/SSS values, almost uncorrelated residuals e j [t], and limited cross-correlation among the inputs u m, j [t] and the residuals e j [t], for each flight j. For the selected structure, the estimation of CCP-NARX coefficients is made via an Ordinary Least Squares (OLS) or Weighted Least Squares (WLS) criterion [16]. Remar 3. The model orders n y and n um in (1), are determined prior to the global CCP-NARX structure selection, via a simple trialand-error procedure: For many individual flights, local NARX(n y, n um ) models are identified (as in the CCP-NARX case) with the maximum initial n y and n um values supported by the given computers. Then, local NARX(n y, n um ) models are identified and validated for reducing n y and n um values, until obtaining a low order but accurate model. Those model orders are then considered for the global CCP-NARX model. Remar 4. The identified CCP-NARX model is valid for a healthy aircraft of given weight and center of gravity position. However, due to fuel consumption the accuracy of the initial CCP-NARX modelmaybealteredbytheendoftheflight.thisproblemmay be trivially solved by creating in the baseline phase CCP-NARX models valid for weight ranges (for instance, full tan, two-thirds, and so on). The tas is trivial since the original CCP-NARX regressors are ept, and only the associated coefficients are instantly re-estimated with flight data from an aircraft operating with the designated weight. Switching among the CCP-NARX models allows for consistently obtaining accurate system representations at the price of the HM scheme being off during the switch. Similar reestimations could also deal with (healthy) engine aging issues. Note that, a definitive solution to this problem involves allowing the (currently constant) CCP-NARX coefficients to be functions of the aircraft weight. Then, the CCP-NARX model becomes a Timedependent Functional Pooled-NARX (TFP-NARX) one, since the aircraft weight is a function of time. Recently, TFP-NARX representations have been used for HM of control surface actuators, with promising results [14]. Remar 5. The computational effort required for a CCP-NARXbased HM scheme is analyzed in [15]: Compared with schemes based on multiple models [21], EKF or UKF[30] any considerable effort is related to the on-the-ground baseline phase (the CCP-NARX identification), rather than the onboard diagnostic one. Crucially, insignificant onboard effort is needed for operating the statistical tests Operational (diagnostic) phase (onboard) For a flight conducted under unnown aircraft engine state (healthy/faulty), the residuals e j [t] (Fig. 1) are collected using the CCP-NARX model identified during the baseline modelling phase. For simplicity, the index j will be dropped hereafter. If the aircraft operates with healthy engines, the residual variance σe 2 is bounded within specific limits. However, this changes when certain engine faults occur: The residual variance increases substantially. Based on this property, a statistical hypothesis test may be designed for the online detection of faults affecting the engine(s). For this purpose, an l-sample long vector [e[t (l 1)]...e[t]], referred to as sliding window, is considered at each time instant t, and the hypothesis testing problem is designed as follows [17]: H 0 : σu 2 [t] σ 2 0 [t] healthy engine H 1 : σu 2 [t] > σ 2 0 [t] faulty engine (4) with H 0 and H 1 designating the null and alternative hypothesis, respectively. The residual variances for the e[t] sequences resulting from the current (unnown) engine state and the healthy state are designated as σu 2[t] and σ 2 0 [t], respectively. The current (unnown health state) variance estimate 2 ˆσ u 2 [t] iscomputedfromthe l-sample long sliding window. Since the current flight is assumed to start with healthy aircraft engines, an estimate of the healthy state variance ˆσ 2 0 is computed using l 0 of the initial samples of the e[t] residuals. The statistic F[t]= ˆσ u 2[t] follows F distribution, and ˆσ 2 0 the following test is designed at ris level α (that is false alarm probability equal to α): F[t] F u H 0 accepted, no fault at time t Else H 1 accepted, fault at time t (5) 2 Symbols with hat designate estimators/estimates.
6 D. Dimogianopoulos et al. / Aerospace Science and Technology 16 (2012) Fig. 4. Statistical hypothesis tests used in the HM scheme: (a) F-test for checing the variance of the residuals e[t] (fault detection); (b) chi-square test for checing the correlation of the residuals e[t] (fault isolation), and (c) normal test for the PACF coefficients (fault isolation). The shaded areas correspond to accepting the null hypothesis H 0 at ris level α. with F u = F (1 α) (l 1,l 0 1) indicating the distribution s (1 α) critical point (see Fig. 4(a)). Following fault detection, isolation is also performed by means of specifically designed hypothesis tests. These statistical tests utilize the lac of correlation property of the residuals e[t], when they result from an aircraft with either healthy engines or engines affected by faults unrelated to those designated for isolation. An l-sample long sliding window [e[t (l 1)]...e[t]] T is formed at each time instant t, and a hypothesis testing problem is designed [31]: H 0 : ρ i = 0 i = 1,...,r (engine not affected by considered fault) H 1 : ρ i 0 for some i (engine affected by considered fault) (6) The terms ρ i designate the autocorrelation coefficients of the sequence e[t] [8, p. 55]. Under the null hypothesis the Q statistic below follows χ 2 (chi-square) distribution with (r s + 1) degrees of freedom: r Q = l ˆρ 2 i χ 2 (r s + 1) (7) i=s with l designating the window length and s indicating the lower sum limit, usually chosen equal to 1 as in [31, pp ]. The hypothesis test at the ris level α, is then formulated as: Q < χ 2 (r s + 1) (1 α) H 0 is accepted Else H 1 is accepted (8) with χ 2 (r s + 1) (1 α) being the value of a statistical quantity following χ 2 distribution with (r s + 1) degrees of freedom and (1 α) critical point (see Fig. 4(b)). The test (8) (referred to as Test A) utilizes(r s + 1) coefficients of the autocorrelation estimate ˆρ i (insteadofthetrueunnownρ i ). An alternative test for fault isolation, based upon the Partial Auto-Correlation Function (PACF) Φ ii of the residual sequence, is also designed. The PACF is the correlation between e[t] and e[t i] after their mutual linear dependency on the intervening variables e[t 1], e[t 2],...,e[t i + 1] has been removed (see [8, pp ] for more details). To set up the test, the variable δφ ii = Φ ii Φ 0 ii is considered, with Φ ii being the partial autocorrelation coefficient at lag i for a data sequence (that is, an l sample long sliding window of e[t]) andφ 0 the empirical mean partial autocorrelation value at lag i when the system is affected by the fault ii designated for isolation. This value is obtained from a number of flights (other than those used for testing) conducted with an aircraft with engines affected by the considered fault. The following composite hypothesis testing problem for the true (but unnown) partial autocorrelation coefficient is then designed: H 0 : δφ ii = 0 (the considered fault occurred) H 1 : δφ ii 0 (another fault occured) (9) In practice, the PACF is obtained by fitting successive AutoRegressive (AR) representations of orders 1, 2,... by OLS and eeping the last coefficient of each regression [8, pp ]. Then ˆΦ ii is equal to the last coefficient of the identified i-th order AR model. Moreover, its variance is given by the (i, i)-th element of the covariance matrix of the estimated AR coefficients. It may be shown that the AR coefficient estimates follow normal distribution [31, pp ]. Hence, treating the computed variance of ˆΦ ii as a fixed quantity, the statistic Z = δσ δ ˆΦ ii 2 follows N (0, 1) distribution under the H 0 hypothesis. Then, a test (referred to as Test B) characterized by ris level α is formed: Z α Z Z 2 1 α H 0 is accepted 2 Else H 1 is accepted (10) with δσ 2 designating the variance of δ ˆΦ ii (which is equal to that of ˆΦ ii ) and Z α the standard normal distribution s α critical point, as in Fig. 4(c).
7 76 D. Dimogianopoulos et al. / Aerospace Science and Technology 16 (2012) Engine health management results A CCP-NARX model is identified upon data from 84 flights of 300 s each, conducted in the clean flight regime under low turbulence with healthy engines. As seen in Fig. 5, these flights evolve in almost all areas of the clean flight regime. The aim is to obtain data representative of the aircraft dynamics in the majority of clean flight operating points. Thus, the identified CCP-NARX model is globally valid throughout this regime. Then, the EHM scheme is extensively tested via 1008 flights conducted in the clean flight regime as follows: For each one of the six health states 84 flights (of 300 s each) are conducted under low and again under moderate turbulence. All single faults (or the first of the sequential faults) occur at around the same time (t = 150 s), that is, well after the beginning but also well before the end of the flight. In all cases the acceleration A y, β angle, rudder moment and thrust signals are recorded at 25 Hz maing for flights of N = 7501 samples. All flights are conducted under slowly changing wind, with commanded heading slowly varying around zero degrees, and with commanded roll tracing a reference signal r[t] = sin(2π 0.01 t), corresponding to slowly changing cruising conditions. Finally, in order to simulate as realistic conditions as possible, the signal measurements are corrupted with noise with its variance being equal to ,5 10 5, and 1 for the A y, β, rudder and thrust signal measurements, respectively (values typically used in relevant studies [9]) Baseline (on-the-ground) phase modelling results Fig. 5. The aircraft clean flight envelope with the trajectories of the 84 flights used in CCP-NARX identification. When identifying a CCP-NARX model, the first tas to be undertaen is related to determining the nonlinearity degree nl. For this purpose, frequency analysis is performed on the set of 84 modelling flights, and those flights with A y, β angle, rudder moment and thrust signals featuring the most intense amplitude variations and rich frequency content are chosen. Then, following Section 3.2, plotsoftheoutput(thrust) amplitude versus that of each of the three inputs (A y, β and rudder) at each of the considered frequencies are formed: at a given frequency, the square root of the output s power spectral density G y (estimated via the Welch method) is plotted for growing values of each input s G um, m = 1, 2, 3, with G y and G um being the amplitudes of the output and the m-th input, respectively. Fig. 6 presents such a plot for the dominant thrust frequency of Hz, identified by the pre- Fig. 6. Estimation of the nonlinearity degree nl via excitation response plots at f = Hz: (a) thrust versus A y,(b)thrust versus β, and(c)thrust versus rudder moment ( 2nd order, 3rd order and 4th order polynomial fits).
8 D. Dimogianopoulos et al. / Aerospace Science and Technology 16 (2012) Table 2 Regressors of the CCP-NARX(7, 4) model with 33 terms. 1 y[t 1] y[t 4] y[t 6] y[t 2] y[t 6] y[t 6] y[t 3] u 3 [t 4] u 3 [t 4] y[t 4] u 2 [t 4] y[t 7] y[t 5] u 1 [t 4] y[t 1] y[t 1] y[t 6] y[t 3] y[t 3] y[t 3] y[t 7] y[t 3] y[t 4] y[t 4] u 3 [t 2] y[t 3] y[t 4] y[t 5] u 3 [t 3] y[t 3] y[t 4] y[t 7] y[t 1] y[t 1] y[t 4] y[t 4] y[t 4] y[t 1] y[t 2] y[t 4] y[t 4] y[t 6] y[t 1] y[t 3] y[t 4] y[t 5] y[t 5] y[t 1] y[t 4] y[t 4] y[t 5] y[t 6] y[t 1] y[t 5] y[t 5] y[t 5] y[t 5] y[t 2] y[t 5] y[t 5] y[t 5] y[t 6] y[t 4] y[t 5] y[t 6] y[t 7] y[t 7] vious frequency analysis. Since G y Gum corresponds to the modulus of the empirical transfer function estimate H y,um at a given frequency (see [31, p. 57]), had the system been linear, H y,um would have been constant for growing input amplitude levels. This is not thecaseinfig. 6 for all three input/output relationships: in fact, the curves fitted correspond equally well to third- or fourth-order polynomials. Hence, according to Section 3.2, nonlinear CCP-NARX regressors with value nl up to three may be considered. For nl 3, the model orders n y and n um (m = 1, 2, 3) are selected via the trial-and-error procedure of Section 3.2. Starting from n y = 20 and n um = 10 (the maximum supported by the given computers), tests with individual flights lead to local NARX(7, 4) models offering the best compromise between low order structure and accurate description of the system dynamics. Nevertheless, for a CCP-NARX(7, 4) with nl 3 the number of possible regressors p i, j [t] in (1), as resulting from all possible combinations of lagged input and output signals, is almost In order to reduce this clearly excessive value without significantly compromising the CCP-NARX model s ability of representing the system dynamics, the four-step structure selection procedure of Section 3.2 is used. Finally a CCP-NARX(7, 4) model with 33 regressors (shown in Table 2) is identified and successfully validated using the correlation tests in [7] On-board operational (diagnostic) phase results and discussion The on-the-ground identified CCP-NARX(7, 4) model with 33 regressors is implemented on-board and used for online obtaining the residual sequence e[t], as specified in Section 3.1. Compromised engine operation is detected by using the statistical hypothesis test (5) on the e[t] sequence. Subsequently, the faults compromising the engine operation are isolated by using either Test A in (8) (for F B faults) or Test B in (10) (for F A faults), also applied to the e[t] sequence. For the case of F D faults, each of 1, 2,t the sequentially occurring single faults are separately detected and isolated via the previous tests. A sliding window of l = 200 samples and a fixed window of l 0 = 200 samples are defined, while the ris level is set to α = The choice of the window length has to optimize the trade-off between minimal delay of detecting malfunctions and the scheme s low sensitivity to false alarms. The latter is achieved via the rather long sliding window currently selected, as minimizing false alarms is a priority. Choosing a long sliding window is possible thans to the CCP-NARX model s capability of accounting for fast emerging nonlinearities in the modelled dynamics. These may occur when the aircraft operates close to the envelope boundaries and/or under significant wind/turbulence levels. With respect Table 3 Fault detection and isolation results for low[moderate] turbulence flights. Aircraft condition Correct detections/ flights Correct detection rates (%) Correct isolations/ flights Correct isolation rates (%) Healthy 82[81]/ [96.4] n/a[n/a] n/a[n/a] F A 25 83[79]/ [94] 83/83[78/79] 100[98.7] F B [80]/ [95.2] 82/82[80/80] 100[100] F C 30 82[81]/ [96.4] n/a[n/a] n/a[n/a] F D 25,550,10 82[80]/ [95.2] 82/82[79/80] 100[98.8] F D 25,550,20 82[79]/ [94] 82/82[78/79] 100[98.7] Regime: clean flight Heading around 0 deg Turbulence: low[moderate] Flights: 504[504] Ris level α: 10 6 Commanded roll r[t] N/a: not applicable. to false alarms due to such nonlinearities, the CCP-NARX based scheme is designed to perform better than linear model based schemes, especially when long sliding windows are used. Obviously, the downside of a long sliding window is the longer fault detection or isolation time. Robustness of the HM scheme is also affected by the choice of the signal part used to form the fixed window during the initial seconds of the flight. It has been noted that the variance of e[t] changes more over the few seconds following the scheme s initialization, than it does later on. In the present case, this initial phase is limited to almost 400 samples (or 16 s). For enhanced robustness, the scheme monitors all successive 200-sample long windows during these 16 s and selects that with the largest variance for computing ˆσ 2 0. Thus, the scheme is operational 16 s after being switched on Detection phase Detection results obtained via the test (5) are presented in Table 3. Fig. 7 shows a typical case with engine malfunction-causing faults detected when the test quantity exceeds the statistical limits corresponding to the α = 10 6 ris level. The HM scheme performs quite well, since the majority of faults are detected and very few false alarms ( 3.5%) are noted. As expected, the scheme s performance for low turbulence flights is better than that obtained for moderate turbulence flights. Obviously, the turbulence increase is a source of extra noise in all signals, which mass part of the fault-induced effects on the left engine dynamics. But even then, the detection performance is not severely degraded, as shown by the only slightly worse false alarm rates in Table 3 resulting from the comparison of low versus moderate turbulence flights: In general, the CCP-NARX(7, 4) based scheme is quite accurate and robust with respect to induced noise with successful detection rates between 94% and 98.8% Isolation phase Isolation results of F A 25 faults obtained via Test B in (10) are presented in Table 3. Note that computing the Z statistic in (10) means that a typical (empirical) value for Φ 0 11 (see Section 3.3) must be nown. The latter is computed from 10 representative flights conducted with a F A 25-affected aircraft. Fig. 8(b) presents a typical case of F A 25 being isolated when the Z statistic enters the isolation zone (delimited by the horizontal dashed lines designating statistical limits at the α = 10 6 ris level). Generally, F A faults show as abrupt increases of the e[t] values followed a few samples later by equally abrupt reductions (spies). This behavior is caused by the autopilot compensating immediately after F A 25 fault occurrence. The PACF coefficient ˆΦ 11 abruptly increases when the sliding window covers the region of e[t] sequence where these spies are present. Finally, since the ˆΦ 11 values stay away from the isolation
9 78 D. Dimogianopoulos et al. / Aerospace Science and Technology 16 (2012) Fig. 7. Fault detection for typical flights under low turbulence and commanded roll tracing of the reference signal r[t]: (a) Healthy; (b) F A 25 affected; (c) F B 550 affected; (d) F C 30 affected; (e) F D 25,550,10 affected aircraft (the horizontal lines designate statistical limits at the α = 10 6 ris level). Fig. 8. Fault isolation with Test A (Q statistic) and Test B (Z statistic) for typical flights under low turbulence and commanded roll tracing of the reference signal r[t]: (a) (b) F A 25 fault isolated via Test B; (c) (d) F B 550 fault isolated via Test A; (e) (f) F C 30 fault not isolated; (g) (h) F D 25,550,10 fault isolated via Test A and Test B (the horizontal lines designate statistical limits at the α = 10 6 ris level). zone when the engine is healthy or affected by F B 550 or F C 30 faults (see Figs. 8(d), (f)), the F A 25 faults cannot be mistaen for any other fault types. Isolation results of F B 550 faults obtained via Test A in (8) are also presented in Table 3. TheQ statistic, involving all lags from the 6th to the 9th, is used in Test A. Fig. 8(c)presentsatypicalcase of F B 550 being isolated when the Q statistic crosses the horizontal dashed line (with the isolation area being over the line designating statistical limit at the α = 10 6 ris level). Generally, F B faults show up as abrupt increases of the e[t] sequence s variance followed by its abrupt reduction samples later, when the left engine shuts down. During these samples, the values of the Q statistic cross the limits (see Fig. 8(c)), and this is exclusively due to the effects of F B fault occurrence: Indeed, for healthy, F A 25 or F C 30 affected engines the Q values are always out of the isolation area (see Figs. 8(a), (e)). Thus, F B 550 faults cannot be mistaen for other fault types. Finally, isolation results for the sequential faults F D 25,550,10 and F D 25,550,20 obtained via Tests A and B are also presented in 3, with typical cases in Figs. 8(g), (h). Note that the isolation of such multiple faults is not a trivial tas: The effects of F B 550 faults on already F A 25 -affected engines (a situation corresponding to F D 25,550,10 or F D 25,550,20 fault occurrence) are quite different than the effects of F B 550 faults on healthy engines. The reason is that the autopilot compensation follows immediately after F A 25 fault occurrence, and, thus, the system dynamics are altered. The difficulty of isolating sequential faults is reflected in the isolation rates for F D 25,550,t
10 D. Dimogianopoulos et al. / Aerospace Science and Technology 16 (2012) faults under demanding circumstances induced by moderate turbulence levels: These are lower (though only slightly) than the rates for pure F B 550 faults. Nonetheless, the CCP-NARX(7, 4) based scheme still remains quite robust leading to good isolation of sequential faults under all circumstances, with rates between 98.7 and 100%. This concludes on the EHM scheme s robustness with respect to external disturbances Comparison with Kalman filter based approaches Due to operating on cruise flight instead of GPA-derived data, the proposed scheme complements rather than replaces the traditional EHM schemes. These, furthermore, focus on detecting and isolating faults affecting specific internal engine parts, whereas the current scheme detects engine malfunctions and isolates possible causes on macroscopic level. A direct comparison between such schemes would, thus, be meaningless. However, meaningful comparisons may be done between schemes sharing the operational concept based on cruise flight data, but featuring different modelling and decision-maing approaches. Using the guidelines in [21] (except for the part of employing GPA data), such schemes may operate with Kalman filters (instead of CCP-NARX modelling) and thresholding of residuals (instead of statistical hypothesis tests). In the absence of an explicit (physics-based or other) model representing the fault-sensitive interrelation among acceleration A y, side-slip angle β, rudder moment and thrust, the standard linear version (LKF) is chosen. This is due to the fact that a linear state-space model may easily be identified via OLS upon the same flight data used for CCP-NARX identification. On the contrary, the use of Extended or Unscented Kalman filters would require a linearized approximation at each operating point of the true nonlinear model of the fault-sensitive interrelation, which, as stated, is not nown. For the same reason, using particle filters based, for instance, on the Sequential Importance Resampling (SIR) for nonlinear state and output estimation via Bayesian techniques (see [5] or [4]) are hardly applicable. Although quite effective, these methods require among others nowledge of the nonlinear (as seen in Section 4.1) state and output generating functions along with the probability density function of the noise affecting the state [4], all of which are currently unavailable. Besides, these methods rely on Monte Carlo recursive techniques to simulate N s vector samples of system state and to execute O(N s ) operations for resampling: All these operations should be carried out online at each time instant [4], and since N s is large for good accuracy, the online computational effort could be significant. In contrast, any considerable computational effort for identifying the CCP-NARX models is invested off-line, and the online effort is insignificant, since mainly associated to running the statistical hypothesis tests for obtaining HM results. Following the approach in [21], the data from the 84 modelling flights are used for identifying via OLS a linear state-space model admitting the same inputs (lateral acceleration A y, side-slip angle β and rudder moment) and covering the same envelope area as the CCP-NARX model. The state variables involve four health parameters, which are measurable sensor outputs and related to potentially isolable faults: The throttle position (related to F A faults), the thrust value (related to F B faults) and the Alt and Mach measurements. Note, that only the first and second measurements are related to modelled faults, which are those considered in the CCP- NARX based case. Equivalently, the unmodelled F C faults are not isolable. The HM scheme utilizes bans of Kalman filters, with each one predicting all state variables using all but one measured sensor output. As in [21], the filter not using the output of a faulty sensor gives the less distorted state prediction, with respect to all other bans. The sensor whose output is not used by that filter is, then, concluded as faulty. Obviously, unlie the CCP-NARX based scheme, sequential faults cannot be trivially isolated, unless the standard KF based scheme is suitably modified. In Fig. 9 the residuals between the measured and the predicted throttle/thrust value, are shown for two different aircraft missions. The first evolves at 7500 ft Mach 0.5 under low turbulence (Figs. 9(a) (d)), and the second at ft Mach 0.5 under moderate turbulence (Figs. 9(e) (h)). Figs. 9(a), (b) and (e), (f) show residuals from the first KF ban which performs state prediction not using throttle sensor measurements, whereas Figs. 9(c), (d) and (g), (h) show residuals from the second KF ban which performs state prediction not using thrust measurements. During the few critical time samples after (abrupt) F A 25 fault occurrence at t = 150 s, the second ban residuals exhibit abnormally large spies with respect to their previous time history. These admit maximum values in Figs. 9(d) and (h), with spies 9 and 5 times larger, respectively, than each signal s past maximum pea. Thus, a threshold on the spies relative magnitude as previously defined, may provide F A 25 fault isolation results under both turbulence levels. However, more turbulence leads to lower relative magnitude of F A 25-caused spies (see also Figs. 9(c) and (g)), and slightly compromised fault isolation. Following F B 550 fault occurrence, the first ban residuals become more noisy than those from the second ban. The noisy evolution is very obvious under low turbulence, as shown in Fig. 9(a). Accordingly, under such conditions F B 550 fault isolation may be achieved when an appropriately selected frequency threshold is exceeded. However, the moderate turbulence induces severe noise even before the F B 550 fault occurrence at t = 150 s, as shown in Figs. 9(e) and (f). Hence, isolating the F B 550 fault via the previous frequency thresholds may prove ineffective for moderate turbulence flights. Even if a frequency threshold specific to the turbulence level could be determined, the problem of distinguishing between noise-induced uncertainty and F B 550 fault occurrence is not trivial. This behavior shows that the LKF is affected by turbulenceinduced noise, or modelling uncertainties, which partly mas the fault effects and undermine the HM scheme s capabilities. In comparison, the CCP-NARX model is more accurate in representing the dynamics throughout the flight regime under noise or modelling uncertainty. Thus, all fault-induced information may be retrieved by the (computationally simple) onboard statistical hypothesis tests. Again, if explicit models of the fault-sensitive interrelation were available, the EKF and UKF could outperform LKF based schemes. Even then, however, the UKF operation requires the square root of the (often ill-conditioned) covariance matrix to be timely and reliably computed onboard at each time instant. Under the current circumstances, the CCP-NARX based scheme proves effective for single and sequentially occurring faults, even in the absence of explicit system model. 5. Conclusions An HM scheme dealing with single and sequentially occurring faults in aircraft engines has been introduced. The scheme uses available cruise flight data, instead of data from gas path analysis, to create additional redundancy, unrelated to that already implemented onboard. Hence, without new HM-dedicated sensors, increased onboard HM decision-maing capabilities for enhanced safety are cost-effectively achieved. The scheme relies on a CCP- NARX(7, 4) representation with 33 regressors to model the MISO relationship among the lateral acceleration A y, side-slip angle β, rudder moment and thrust signals. The CCP-NARX model is identified on-the-ground via relevant past available data. Thus, unlie schemes using extended/unscented Kalman (or even particle) filters, the scheme may use any suitable fault-sensitive relationship for HM purposes even if a (physics-based or other) model
Onboard Engine FDI in Autonomous Aircraft Using Stochastic Nonlinear Modelling of Flight Signal Dependencies
Onboard Engine FDI in Autonomous Aircraft Using Stochastic Nonlinear Modelling of Flight Signal Dependencies Dimitrios G. Dimogianopoulos, John D. Hios and Spilios D. Fassois Stochastic Mechanical Systems
More informationNONLINEAR INTEGRAL MINIMUM VARIANCE-LIKE CONTROL WITH APPLICATION TO AN AIRCRAFT SYSTEM
NONLINEAR INTEGRAL MINIMUM VARIANCE-LIKE CONTROL WITH APPLICATION TO AN AIRCRAFT SYSTEM D.G. Dimogianopoulos, J.D. Hios and S.D. Fassois DEPARTMENT OF MECHANICAL & AERONAUTICAL ENGINEERING GR-26500 PATRAS,
More informationVibration Based Health Monitoring for a Thin Aluminum Plate: Experimental Assessment of Several Statistical Time Series Methods
Vibration Based Health Monitoring for a Thin Aluminum Plate: Experimental Assessment of Several Statistical Time Series Methods Fotis P. Kopsaftopoulos and Spilios D. Fassois Stochastic Mechanical Systems
More informationStationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study
Stationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study Andriana S. GEORGANTOPOULOU & Spilios D. FASSOIS Stochastic Mechanical Systems & Automation
More informationParametric Output Error Based Identification and Fault Detection in Structures Under Earthquake Excitation
Parametric Output Error Based Identification and Fault Detection in Structures Under Earthquake Excitation J.S. Sakellariou and S.D. Fassois Department of Mechanical & Aeronautical Engr. GR 265 Patras,
More informationthe Functional Model Based Method
Multi-Site Damage Localization via the Functional Model Based Method Christos S. Sakaris, John S. Sakellariou and Spilios D. Fassois Stochastic Mechanical Systems & Automation (SMSA) Laboratory Department
More informationIn-Flight Engine Diagnostics and Prognostics Using A Stochastic-Neuro-Fuzzy Inference System
In-Flight Engine Diagnostics and Prognostics Using A Stochastic-Neuro-Fuzzy Inference System Dan M. Ghiocel & Joshua Altmann STI Technologies, Rochester, New York, USA Keywords: reliability, stochastic
More informationT.-C.J. Aravanis, J.S. Sakellariou and S.D. Fassois
Vibration based fault detection under variable non-measurable, operating conditions via a stochastic Functional Model method and application to railway vehicle suspensions T.-C.J. Aravanis, J.S. Sakellariou
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi Lecture 13-29 May, 2013 Courtesy of Prof. S. Fassois & Dr. F. Kopsaftopoulos, SMSA Group, University of Patras Statistical methods for SHM courtesy of Prof. S. Fassois & Dr. F. Kopsaftopoulos,
More informationVibration-Response-Based Damage Detection For Wind Turbine Blades Under Varying Environmental Conditions
Vibration-Response-Based Damage Detection For Wind Turbine Blades Under Varying Environmental Conditions Ana Gómez González Spilios D. Fassois Stochastic Mechanical Systems & Automation (SMSA) Laboratory
More informationA GPC-based Auto-Throttle with Robust Fuel-efficient Operation in 4D Flights
A GPC-based Auto-Throttle with Robust Fuel-efficient Operation in 4D Flights Dimitrios Dimogianopoulos a,, Vangelis Petratos b, Fotis Kopsaftopoulos b, Spilios Fassois b a Technological Education Institute
More informationIntegral Minimum Variance-Like Control for Pooled Nonlinear Representations with Application to an Aircraft System
International Journal of Control Vol. 00, No. 00, DD Month 200x, 1 18 This is a preprint of an article whose final and definitive form has been published in the International Journal of Control 2007 copyright
More informationMini-Course 07 Kalman Particle Filters. Henrique Massard da Fonseca Cesar Cunha Pacheco Wellington Bettencurte Julio Dutra
Mini-Course 07 Kalman Particle Filters Henrique Massard da Fonseca Cesar Cunha Pacheco Wellington Bettencurte Julio Dutra Agenda State Estimation Problems & Kalman Filter Henrique Massard Steady State
More informationAdaptive Unscented Kalman Filter with Multiple Fading Factors for Pico Satellite Attitude Estimation
Adaptive Unscented Kalman Filter with Multiple Fading Factors for Pico Satellite Attitude Estimation Halil Ersin Söken and Chingiz Hajiyev Aeronautics and Astronautics Faculty Istanbul Technical University
More informationA WAVELET BASED FLIGHT DATA PREPROCESSING METHOD FOR FLIGHT CHARACTERISTICS ESTIMATION AND FAULT DETECTION
A WAVELET BASED FLIGHT DATA PREPROCESSING METHOD FOR FLIGHT CHARACTERISTICS ESTIMATION AND FAULT DETECTION Masaru Naruoka Japan Aerospace Exploration Agency Keywords: Flight analyses, Multiresolution Analysis,
More informationStationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study
6th International Symposium on NDT in Aerospace, 12-14th November 2014, Madrid, Spain - www.ndt.net/app.aerondt2014 More Info at Open Access Database www.ndt.net/?id=16938 Stationary or Non-Stationary
More informationFailure Prognostics with Missing Data Using Extended Kalman Filter
Failure Prognostics with Missing Data Using Extended Kalman Filter Wlamir Olivares Loesch Vianna 1, and Takashi Yoneyama 2 1 EMBRAER S.A., São José dos Campos, São Paulo, 12227 901, Brazil wlamir.vianna@embraer.com.br
More informationIdentification of Stochastic Systems Under Multiple Operating Conditions: The Vector Dependent FP ARX Parametrization
Identification of Stochastic Systems Under Multiple Operating Conditions: The Vector Dependent FP ARX Parametrization Fotis P Kopsaftopoulos and Spilios D Fassois Abstract The problem of identifying stochastic
More informationFAULT DETECTION AND FAULT TOLERANT APPROACHES WITH AIRCRAFT APPLICATION. Andrés Marcos
FAULT DETECTION AND FAULT TOLERANT APPROACHES WITH AIRCRAFT APPLICATION 2003 Louisiana Workshop on System Safety Andrés Marcos Dept. Aerospace Engineering and Mechanics, University of Minnesota 28 Feb,
More informationAFAULT diagnosis procedure is typically divided into three
576 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 47, NO. 4, APRIL 2002 A Robust Detection and Isolation Scheme for Abrupt and Incipient Faults in Nonlinear Systems Xiaodong Zhang, Marios M. Polycarpou,
More informationA Robust Approach to Estimating Production Functions: Replication of the ACF procedure
A Robust Approach to Estimating Production Functions: Replication of the ACF procedure Kyoo il Kim Michigan State University Yao Luo University of Toronto Yingjun Su IESR, Jinan University August 2018
More informationA. Poulimenos, M. Spiridonakos, and S. Fassois
PARAMETRIC TIME-DOMAIN METHODS FOR NON-STATIONARY RANDOM VIBRATION IDENTIFICATION AND ANALYSIS: AN OVERVIEW AND COMPARISON A. Poulimenos, M. Spiridonakos, and S. Fassois DEPARTMENT OF MECHANICAL & AERONAUTICAL
More informationIntroduction to System Identification and Adaptive Control
Introduction to System Identification and Adaptive Control A. Khaki Sedigh Control Systems Group Faculty of Electrical and Computer Engineering K. N. Toosi University of Technology May 2009 Introduction
More information3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE
3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE 3.0 INTRODUCTION The purpose of this chapter is to introduce estimators shortly. More elaborated courses on System Identification, which are given
More information2 Introduction of Discrete-Time Systems
2 Introduction of Discrete-Time Systems This chapter concerns an important subclass of discrete-time systems, which are the linear and time-invariant systems excited by Gaussian distributed stochastic
More informationVARIANCE COMPUTATION OF MODAL PARAMETER ES- TIMATES FROM UPC SUBSPACE IDENTIFICATION
VARIANCE COMPUTATION OF MODAL PARAMETER ES- TIMATES FROM UPC SUBSPACE IDENTIFICATION Michael Döhler 1, Palle Andersen 2, Laurent Mevel 1 1 Inria/IFSTTAR, I4S, Rennes, France, {michaeldoehler, laurentmevel}@inriafr
More informationCOMBINED ADAPTIVE CONTROLLER FOR UAV GUIDANCE
COMBINED ADAPTIVE CONTROLLER FOR UAV GUIDANCE B.R. Andrievsky, A.L. Fradkov Institute for Problems of Mechanical Engineering of Russian Academy of Sciences 61, Bolshoy av., V.O., 199178 Saint Petersburg,
More informationInteracting multiple-models, state augmented Particle Filtering for Fault Diagnostics
Interacting multiple-models, state augmented Particle Filtering for Fault Diagnostics Michele Compare 1,2, Piero Baraldi 1, Pietro Turati 1, Enrico Zio 1,2,3,* 1 Energy Department, Politecnico di Milano,
More informationProf. 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
More informationComparision of Probabilistic Navigation methods for a Swimming Robot
Comparision of Probabilistic Navigation methods for a Swimming Robot Anwar Ahmad Quraishi Semester Project, Autumn 2013 Supervisor: Yannic Morel BioRobotics Laboratory Headed by Prof. Aue Jan Ijspeert
More informationOUTPUT-ONLY STATISTICAL TIME SERIES METHODS FOR STRUCTURAL HEALTH MONITORING: A COMPARATIVE STUDY
7th European Workshop on Structural Health Monitoring July 8-11, 2014. La Cité, Nantes, France More Info at Open Access Database www.ndt.net/?id=17198 OUTPUT-ONLY STATISTICAL TIME SERIES METHODS FOR STRUCTURAL
More informationModel-based Stochastic Fault Detection and Diagnosis for Lithium-ion Batteries
Article Model-based Stochastic Fault Detection and Diagnosis for Lithium-ion Batteries Jeongeun Son, Yuncheng Du* Department of Chemical & Biomolecular Engineering, Clarkson University, Potsdam NY 13676,
More informationA nonlinear filtering tool for analysis of hot-loop test campaings
A nonlinear filtering tool for analysis of hot-loop test campaings Enso Ikonen* Jenő Kovács*, ** * Systems Engineering Laboratory, Department of Process and Environmental Engineering, University of Oulu,
More informationRAO-BLACKWELLISED PARTICLE FILTERS: EXAMPLES OF APPLICATIONS
RAO-BLACKWELLISED PARTICLE FILTERS: EXAMPLES OF APPLICATIONS Frédéric Mustière e-mail: mustiere@site.uottawa.ca Miodrag Bolić e-mail: mbolic@site.uottawa.ca Martin Bouchard e-mail: bouchard@site.uottawa.ca
More informationSet-valued Observer-based Active Fault-tolerant Model Predictive Control
Set-valued Observer-based Active Fault-tolerant Model Predictive Control Feng Xu 1,2, Vicenç Puig 1, Carlos Ocampo-Martinez 1 and Xueqian Wang 2 1 Institut de Robòtica i Informàtica Industrial (CSIC-UPC),Technical
More informationA NONLINEARITY MEASURE FOR ESTIMATION SYSTEMS
AAS 6-135 A NONLINEARITY MEASURE FOR ESTIMATION SYSTEMS Andrew J. Sinclair,JohnE.Hurtado, and John L. Junkins The concept of nonlinearity measures for dynamical systems is extended to estimation systems,
More informationEconomics 536 Lecture 7. Introduction to Specification Testing in Dynamic Econometric Models
University of Illinois Fall 2016 Department of Economics Roger Koenker Economics 536 Lecture 7 Introduction to Specification Testing in Dynamic Econometric Models In this lecture I want to briefly describe
More informationTAKEHOME FINAL EXAM e iω e 2iω e iω e 2iω
ECO 513 Spring 2015 TAKEHOME FINAL EXAM (1) Suppose the univariate stochastic process y is ARMA(2,2) of the following form: y t = 1.6974y t 1.9604y t 2 + ε t 1.6628ε t 1 +.9216ε t 2, (1) where ε is i.i.d.
More informationFault Diagnosis and Estimation of Dynamical Systems with Application to Gas Turbines
Fault Diagnosis and Estimation of Dynamical Systems with Application to Gas Turbines Esmaeil Naderi A Thesis in The Department of Electrical and Computer Engineering Presented in Partial Fulfillment of
More informationBox-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
More informationIdentifying the Monetary Policy Shock Christiano et al. (1999)
Identifying the Monetary Policy Shock Christiano et al. (1999) The question we are asking is: What are the consequences of a monetary policy shock a shock which is purely related to monetary conditions
More informationAvailable online at ScienceDirect. Procedia Engineering 119 (2015 ) 13 18
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 119 (2015 ) 13 18 13th Computer Control for Water Industry Conference, CCWI 2015 Real-time burst detection in water distribution
More informationImpeller Fault Detection for a Centrifugal Pump Using Principal Component Analysis of Time Domain Vibration Features
Impeller Fault Detection for a Centrifugal Pump Using Principal Component Analysis of Time Domain Vibration Features Berli Kamiel 1,2, Gareth Forbes 2, Rodney Entwistle 2, Ilyas Mazhar 2 and Ian Howard
More informationBootstrapping Heteroskedasticity Consistent Covariance Matrix Estimator
Bootstrapping Heteroskedasticity Consistent Covariance Matrix Estimator by Emmanuel Flachaire Eurequa, University Paris I Panthéon-Sorbonne December 2001 Abstract Recent results of Cribari-Neto and Zarkos
More informationVEHICLE WHEEL-GROUND CONTACT ANGLE ESTIMATION: WITH APPLICATION TO MOBILE ROBOT TRACTION CONTROL
1/10 IAGNEMMA AND DUBOWSKY VEHICLE WHEEL-GROUND CONTACT ANGLE ESTIMATION: WITH APPLICATION TO MOBILE ROBOT TRACTION CONTROL K. IAGNEMMA S. DUBOWSKY Massachusetts Institute of Technology, Cambridge, MA
More informationMultivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8]
1 Multivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8] Insights: Price movements in one market can spread easily and instantly to another market [economic globalization and internet
More informationBasics of Uncertainty Analysis
Basics of Uncertainty Analysis Chapter Six Basics of Uncertainty Analysis 6.1 Introduction As shown in Fig. 6.1, analysis models are used to predict the performances or behaviors of a product under design.
More informationRobotics 2 Target Tracking. Kai Arras, Cyrill Stachniss, Maren Bennewitz, Wolfram Burgard
Robotics 2 Target Tracking Kai Arras, Cyrill Stachniss, Maren Bennewitz, Wolfram Burgard Slides by Kai Arras, Gian Diego Tipaldi, v.1.1, Jan 2012 Chapter Contents Target Tracking Overview Applications
More informationSequential Monte Carlo methods for filtering of unobservable components of multidimensional diffusion Markov processes
Sequential Monte Carlo methods for filtering of unobservable components of multidimensional diffusion Markov processes Ellida M. Khazen * 13395 Coppermine Rd. Apartment 410 Herndon VA 20171 USA Abstract
More informationHW VS SW SENSOR REDUNDANCY: FAULT DETECTION AND ISOLATION OBSERVER BASED APPROACHES FOR INERTIAL MEASUREMENT UNITS
HW VS SW SENSOR REDUNDANCY: FAULT DETECTION AND ISOLATION OBSERVER BASED APPROACHES FOR INERTIAL MEASUREMENT UNITS Abstract Immacolata Notaro*, Marco Ariola**, Egidio D Amato*, Massimiliano Mattei* *Dipartimento
More informationGAS TURBINES FAULT IDENTIFICATION BY FUSING VIBRATION TRENDING AND GAS PATH ANALYSIS
GAS TURBINES FAULT IDENTIFICATION BY FUSING VIBRATION TRENDING AND A. Kyriazis Research Assistant A. Tsalavoutas Research Associate K. Mathioudakis Professor Laboratory of Thermal Turbomachines National
More informationKeywords: Multimode process monitoring, Joint probability, Weighted probabilistic PCA, Coefficient of variation.
2016 International Conference on rtificial Intelligence: Techniques and pplications (IT 2016) ISBN: 978-1-60595-389-2 Joint Probability Density and Weighted Probabilistic PC Based on Coefficient of Variation
More informationUsing the Kalman Filter to Estimate the State of a Maneuvering Aircraft
1 Using the Kalman Filter to Estimate the State of a Maneuvering Aircraft K. Meier and A. Desai Abstract Using sensors that only measure the bearing angle and range of an aircraft, a Kalman filter is implemented
More informationNonlinear State Estimation Methods Overview and Application to PET Polymerization
epartment of Biochemical and Chemical Engineering Process ynamics Group () onlinear State Estimation Methods Overview and Polymerization Paul Appelhaus Ralf Gesthuisen Stefan Krämer Sebastian Engell The
More informationEEG- Signal Processing
Fatemeh Hadaeghi EEG- Signal Processing Lecture Notes for BSP, Chapter 5 Master Program Data Engineering 1 5 Introduction The complex patterns of neural activity, both in presence and absence of external
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More informationA PARAMETRIC INVESTIGATION OF THE DIAGNOSTIC ABILITY OF PROBABILISTIC NEURAL NETWORKS ON TURBOFAN ENGINES
A PARAMETRIC INVESTIGATION OF THE DIAGNOSTIC ABILITY OF PROBABILISTIC NEURAL NETWORKS ON TURBOFAN ENGINES C. Romessis Research Assistant A. Stamatis Research associate K. Mathioudakis Assistant Professor
More informationChapter 12: An introduction to Time Series Analysis. Chapter 12: An introduction to Time Series Analysis
Chapter 12: An introduction to Time Series Analysis Introduction In this chapter, we will discuss forecasting with single-series (univariate) Box-Jenkins models. The common name of the models is Auto-Regressive
More informationEvolution Strategies Based Particle Filters for Fault Detection
Evolution Strategies Based Particle Filters for Fault Detection Katsuji Uosaki, Member, IEEE, and Toshiharu Hatanaka, Member, IEEE Abstract Recent massive increase of the computational power has allowed
More informationA BAYESIAN APPROACH FOR PREDICTING BUILDING COOLING AND HEATING CONSUMPTION
A BAYESIAN APPROACH FOR PREDICTING BUILDING COOLING AND HEATING CONSUMPTION Bin Yan, and Ali M. Malkawi School of Design, University of Pennsylvania, Philadelphia PA 19104, United States ABSTRACT This
More informationA NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1
A NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1 Jinglin Zhou Hong Wang, Donghua Zhou Department of Automation, Tsinghua University, Beijing 100084, P. R. China Control Systems Centre,
More informationA COMPARISON OF TWO METHODS FOR STOCHASTIC FAULT DETECTION: THE PARITY SPACE APPROACH AND PRINCIPAL COMPONENTS ANALYSIS
A COMPARISON OF TWO METHODS FOR STOCHASTIC FAULT DETECTION: THE PARITY SPACE APPROACH AND PRINCIPAL COMPONENTS ANALYSIS Anna Hagenblad, Fredrik Gustafsson, Inger Klein Department of Electrical Engineering,
More information1 Motivation for Instrumental Variable (IV) Regression
ECON 370: IV & 2SLS 1 Instrumental Variables Estimation and Two Stage Least Squares Econometric Methods, ECON 370 Let s get back to the thiking in terms of cross sectional (or pooled cross sectional) data
More informationDiscussion of Bootstrap prediction intervals for linear, nonlinear, and nonparametric autoregressions, by Li Pan and Dimitris Politis
Discussion of Bootstrap prediction intervals for linear, nonlinear, and nonparametric autoregressions, by Li Pan and Dimitris Politis Sílvia Gonçalves and Benoit Perron Département de sciences économiques,
More informationAPPLICATION OF A KERNEL METHOD IN MODELING FRICTION DYNAMICS
APPLICATION OF A KERNEL METHOD IN MODELING FRICTION DYNAMICS Yufeng Wan, Chian X. Wong, Tony J. Dodd, Robert F. Harrison Department of Automatic Control and Systems Engineering, The University of Sheffield,
More informationFAULT - TOLERANT PROCEDURES FOR AIR DATA ELABORATION
25 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES FAULT - TOLERANT PROCEDURES FOR AIR DATA ELABORATION Alberto Calia, Eugenio Denti, Roberto Galatolo, Francesco Schettini University of Pisa Department
More informationA SELF-TUNING KALMAN FILTER FOR AUTONOMOUS SPACECRAFT NAVIGATION
A SELF-TUNING KALMAN FILTER FOR AUTONOMOUS SPACECRAFT NAVIGATION Son H. Truong National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC) Greenbelt, Maryland, USA 2771 E-mail:
More informationQuadratic Extended Filtering in Nonlinear Systems with Uncertain Observations
Applied Mathematical Sciences, Vol. 8, 2014, no. 4, 157-172 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ams.2014.311636 Quadratic Extended Filtering in Nonlinear Systems with Uncertain Observations
More informationState Estimation by IMM Filter in the Presence of Structural Uncertainty 1
Recent Advances in Signal Processing and Communications Edited by Nios Mastorais World Scientific and Engineering Society (WSES) Press Greece 999 pp.8-88. State Estimation by IMM Filter in the Presence
More informationLecture 9. Time series prediction
Lecture 9 Time series prediction Prediction is about function fitting To predict we need to model There are a bewildering number of models for data we look at some of the major approaches in this lecture
More informationModule 3. Descriptive Time Series Statistics and Introduction to Time Series Models
Module 3 Descriptive Time Series Statistics and Introduction to Time Series Models Class notes for Statistics 451: Applied Time Series Iowa State University Copyright 2015 W Q Meeker November 11, 2015
More informationState estimation of linear dynamic system with unknown input and uncertain observation using dynamic programming
Control and Cybernetics vol. 35 (2006) No. 4 State estimation of linear dynamic system with unknown input and uncertain observation using dynamic programming by Dariusz Janczak and Yuri Grishin Department
More informationHowever, reliability analysis is not limited to calculation of the probability of failure.
Probabilistic Analysis probabilistic analysis methods, including the first and second-order reliability methods, Monte Carlo simulation, Importance sampling, Latin Hypercube sampling, and stochastic expansions
More informationChapter 2 Review of Linear and Nonlinear Controller Designs
Chapter 2 Review of Linear and Nonlinear Controller Designs This Chapter reviews several flight controller designs for unmanned rotorcraft. 1 Flight control systems have been proposed and tested on a wide
More informationA Multiple-Observer Scheme for Fault Detection, Isolation and Recovery of Satellite Thrusters
Proceedings of the EuroGNC 3, nd CEAS Specialist Conference on Guidance, Navigation & Control, Delft University of Technology, Delft, The Netherlands, April -, 3 FrBT.3 A Multiple-Observer Scheme for Fault
More informationAdvanced Adaptive Control for Unintended System Behavior
Advanced Adaptive Control for Unintended System Behavior Dr. Chengyu Cao Mechanical Engineering University of Connecticut ccao@engr.uconn.edu jtang@engr.uconn.edu Outline Part I: Challenges: Unintended
More informationMMSE-Based Filtering for Linear and Nonlinear Systems in the Presence of Non-Gaussian System and Measurement Noise
MMSE-Based Filtering for Linear and Nonlinear Systems in the Presence of Non-Gaussian System and Measurement Noise I. Bilik 1 and J. Tabrikian 2 1 Dept. of Electrical and Computer Engineering, University
More informationThe Unscented Particle Filter
The Unscented Particle Filter Rudolph van der Merwe (OGI) Nando de Freitas (UC Bereley) Arnaud Doucet (Cambridge University) Eric Wan (OGI) Outline Optimal Estimation & Filtering Optimal Recursive Bayesian
More informationFAULT DETECTION for SPACECRAFT ATTITUDE CONTROL SYSTEM. M. Amin Vahid D. Mechanical Engineering Department Concordia University December 19 th, 2010
FAULT DETECTION for SPACECRAFT ATTITUDE CONTROL SYSTEM M. Amin Vahid D. Mechanical Engineering Department Concordia University December 19 th, 2010 Attitude control : the exercise of control over the orientation
More informationLecture 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:
More information7. Integrated Processes
7. Integrated Processes Up to now: Analysis of stationary processes (stationary ARMA(p, q) processes) Problem: Many economic time series exhibit non-stationary patterns over time 226 Example: We consider
More informationNew Statistical Model for the Enhancement of Noisy Speech
New Statistical Model for the Enhancement of Noisy Speech Electrical Engineering Department Technion - Israel Institute of Technology February 22, 27 Outline Problem Formulation and Motivation 1 Problem
More informationEfficient Monitoring for Planetary Rovers
International Symposium on Artificial Intelligence and Robotics in Space (isairas), May, 2003 Efficient Monitoring for Planetary Rovers Vandi Verma vandi@ri.cmu.edu Geoff Gordon ggordon@cs.cmu.edu Carnegie
More informationEconometrics of Panel Data
Econometrics of Panel Data Jakub Mućk Meeting # 9 Jakub Mućk Econometrics of Panel Data Meeting # 9 1 / 22 Outline 1 Time series analysis Stationarity Unit Root Tests for Nonstationarity 2 Panel Unit Root
More informationA CUSUM approach for online change-point detection on curve sequences
ESANN 22 proceedings, European Symposium on Artificial Neural Networks, Computational Intelligence and Machine Learning. Bruges Belgium, 25-27 April 22, i6doc.com publ., ISBN 978-2-8749-49-. Available
More informationModel-Based Diagnosis of Chaotic Vibration Signals
Model-Based Diagnosis of Chaotic Vibration Signals Ihab Wattar ABB Automation 29801 Euclid Ave., MS. 2F8 Wickliffe, OH 44092 and Department of Electrical and Computer Engineering Cleveland State University,
More informationA Non-Parametric Approach of Heteroskedasticity Robust Estimation of Vector-Autoregressive (VAR) Models
Journal of Finance and Investment Analysis, vol.1, no.1, 2012, 55-67 ISSN: 2241-0988 (print version), 2241-0996 (online) International Scientific Press, 2012 A Non-Parametric Approach of Heteroskedasticity
More informationComplexity Metrics. ICRAT Tutorial on Airborne self separation in air transportation Budapest, Hungary June 1, 2010.
Complexity Metrics ICRAT Tutorial on Airborne self separation in air transportation Budapest, Hungary June 1, 2010 Outline Introduction and motivation The notion of air traffic complexity Relevant characteristics
More informationTesting an Autoregressive Structure in Binary Time Series Models
ömmföäflsäafaäsflassflassflas ffffffffffffffffffffffffffffffffffff Discussion Papers Testing an Autoregressive Structure in Binary Time Series Models Henri Nyberg University of Helsinki and HECER Discussion
More informationUsing all observations when forecasting under structural breaks
Using all observations when forecasting under structural breaks Stanislav Anatolyev New Economic School Victor Kitov Moscow State University December 2007 Abstract We extend the idea of the trade-off window
More informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
More informationSimultaneous state and input estimation with partial information on the inputs
Loughborough University Institutional Repository Simultaneous state and input estimation with partial information on the inputs This item was submitted to Loughborough University's Institutional Repository
More informationTrajectory planning and feedforward design for electromechanical motion systems version 2
2 Trajectory planning and feedforward design for electromechanical motion systems version 2 Report nr. DCT 2003-8 Paul Lambrechts Email: P.F.Lambrechts@tue.nl April, 2003 Abstract This report considers
More informationNoise Robust Isolated Words Recognition Problem Solving Based on Simultaneous Perturbation Stochastic Approximation Algorithm
EngOpt 2008 - International Conference on Engineering Optimization Rio de Janeiro, Brazil, 0-05 June 2008. Noise Robust Isolated Words Recognition Problem Solving Based on Simultaneous Perturbation Stochastic
More informationMech 6091 Flight Control System Course Project. Team Member: Bai, Jing Cui, Yi Wang, Xiaoli
Mech 6091 Flight Control System Course Project Team Member: Bai, Jing Cui, Yi Wang, Xiaoli Outline 1. Linearization of Nonlinear F-16 Model 2. Longitudinal SAS and Autopilot Design 3. Lateral SAS and Autopilot
More information7. Integrated Processes
7. Integrated Processes Up to now: Analysis of stationary processes (stationary ARMA(p, q) processes) Problem: Many economic time series exhibit non-stationary patterns over time 226 Example: We consider
More informationPrediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother
Prediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother Simo Särkkä, Aki Vehtari and Jouko Lampinen Helsinki University of Technology Department of Electrical and Communications
More informationEUSIPCO
EUSIPCO 3 569736677 FULLY ISTRIBUTE SIGNAL ETECTION: APPLICATION TO COGNITIVE RAIO Franc Iutzeler Philippe Ciblat Telecom ParisTech, 46 rue Barrault 753 Paris, France email: firstnamelastname@telecom-paristechfr
More informationAR, 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
More informationMultivariate Time Series
Multivariate Time Series Notation: I do not use boldface (or anything else) to distinguish vectors from scalars. Tsay (and many other writers) do. I denote a multivariate stochastic process in the form
More information