Acessibilidade / Reportar erro

MONITORING THE STATOR CURRENT IN INDUCTION MACHINES FOR POSSIBLE FAULT DETECTION: A FUZZY/BAYESIAN APPROACH FOR THE PROBLEM OF TIME SERIES MULTIPLE CHANGE POINT DETECTION

ABSTRACT

This paper addresses the problem of fault detection in stator winding of induction machine by a multiple change points detection approach in time series. To handle this problem a new fuzzy/Bayesian approach is proposed which differs from previous approaches since it does not require prior information as: the number of change points or the characterization of the data probabilistic distribution. The approach has been applied in the monitoring the current of the stator winding induction machine. The good results obtained by proposed methodology illustrate its efficiency.

Keywords:
Fuzzy/Bayesian; multiple change points detection; fault monitoring

1 INTRODUCTION

Fault detection and analysis is a very important strategy that is commonly employed in the industry with the purpose of allowing a cost-effective maintenance policy, keeping productivity standards and ensuring safety. The fault analysis gives support for the design of corrective actions, system redundancies, and safety policies in order to mitigate the effects of a fault3333 ISERMANN R & BALLE P. 1997. Trends in the application of model-based fault detection and diagnosis of technical processes. Control Engineering Practice, 5(5): 707-719.. A fault diagnosis procedure is typically divided into three tasks:

  1. fault detection, indicating the occurrence of some fault in a monitored system;

  2. fault isolation, establishing the type and/or location of the fault; and

  3. fault identification, determining the magnitude of the fault.

After a fault has been detected and diagnosed, in some applications it is required that the fault be self-corrected, usually through controller reconfiguration. This is usually referred to as fault accommodation.

The literature presents several classes of strategies to deal with fault detection and isolation (FDI)1414 CHEN J & PATTON RJ. 1999. Robust model-based fault diagnosis for dynamic systems. Dordrecht: Kluwer Academic Publishers.. These strategies can be, in general, divided in approaches based on quantitative models6868 VENKATASUBRAMANIAN V, RENGASWAMY R, YIN K & KAVURI SN. 2003. A review of process fault detection and diagnosis - part I: Quantitative model-based methods. Computers and Chemical Engineering, 27: 293-311. and on qualitative models6666 VENKATASUBRAMANIAN V, RENGASWAMY R & KAVURI SN. 2003. A review of process fault detection and diagnosis - part II: Qualitative models and search strategies. Computers and Chemical Engineering, 27: 313-326.), (6767 VENKATASUBRAMANIAN V, RENGASWAMY R, KAVURI SN & YIN K. 2003. A review of process fault detection and diagnosis - part III: Process history based methods. Computers and Chemical Engineering, 27: 327-346..

Most of the quantitative model-based approaches are based on the knowledge of mathematical models of the plant. Many survey papers with different emphasis on various model-based approaches have been published over the past years. The main approaches in this context are based on (unknown input) observers1414 CHEN J & PATTON RJ. 1999. Robust model-based fault diagnosis for dynamic systems. Dordrecht: Kluwer Academic Publishers.), (1515 CHEN W & SAIF M. 2007. Observer-based strategies for actuator fault detection, isolation and estimation for certain class of uncertain nonlinear systems. IET Control Theoty and Applications, 1(6): 1672-1680.), (1616 CAMINHAS WM & TAKAHASHI RHC. 2001. Dynamic system failure detection and diagnosis employing sliding mode observers and fuzzy neural networks. In Proceedings of the Joint 9th IFSA and 20th NAFIPS, pages 304-309, Vancouver.), (5252 PUIG V, STANCU A, ESCOBET T, NEJJARI F, QUEVEDO J & PATTON RJ. 2006. Passive robust fault detection using interval observers: Application to the DAMADICS benchmark problem. Control Engineering Practice, 14(6): 621-633.), (6262 TAKAHASHI RHC & PERES PLD. 1999. Unknown input observers for uncertain systems: A unifying approach. European Journal of Control, 5(2-4): 261-275.), (6363 TAKAHASHI RHC, PALHARES RM & PERES PLD. 1999. Discrete-time singular observers: H2/H∞ optimality and unknown inputs. International Journal of Control, 72(6): 481-492., parity relations1414 CHEN J & PATTON RJ. 1999. Robust model-based fault diagnosis for dynamic systems. Dordrecht: Kluwer Academic Publishers.), (5050 PLOIX S & ADROT O. 2006. Parity relations for linear uncertain dynamic systems. Automatica, 42(9): 1553-1562. and Kalman or robust filters11 ALAVI SMM, IZADI-ZAMANABADI R & HAYES MJ. 2008. Robust fault detection and isolation technique for single-input/single-output closed-loop control systems that exhibit actuator and sensor faults. IET Control Theoty and Applications, 2(11): 951-965.), (1414 CHEN J & PATTON RJ. 1999. Robust model-based fault diagnosis for dynamic systems. Dordrecht: Kluwer Academic Publishers.), (2828 GAO H, CHEN T & WANG L. 2008. Robust fault detection with missing measurements. International Journal of Control, 81(5): 804-819.), (2929 GAO Z, WANG H & CHAI T. 2007. A robust fault detection filtering for stochastic distribution systems via descriptor estimator and parametric gain design. IET Control Theoty and Applications, 1(5): 1286-1293.), (7272 ZOLGHADRI A. 1996. An algorithm for real-time failure detection in Kalman filters. IEEE Transactions on Automatic Control, 41(10): 1537-1539.. The requirement of a mathematical model of the plant can lead to several difficulties in the implementation of these approaches, for instance due to factors such as system complexity, high dimensionality, nonlinearities and parametric uncertainties.

On the other hand, most of the qualitative model approaches are based on some pattern analysis of the historical process data. The main related approaches are: signed directed graph55 BOUKHOBZA T, HAMELIN F & CANITROT S. 2008. A graph-theoretic approach to fault detection and isolation for structured bilinear systems. International Journal of Control, 81(4): 661-678.), (1313 CHENG H, NIKUS M & JAMSA-JOUNELA S. 2008. Fault diagnosis of the paper machine short circulation process using novel dynamic causal digraph reasoning. Journal of Process Control, 18(7-8): 676-691.), (4747 MAURYA MR, RENGASWAMY R & VENKATASUBRAMANIAN V. 2006. A signed directed graph-based systematic framework for steady-state malfunction diagnosis inside control loops. Chemical Engineering Science, 61(6): 1790-1810., fault tree2323 DUTUIT Y & RAUZY A. 2005. Approximate estimation of system reliability via fault trees. Reliability Engineering & System Safety, 87(2): 163-172., fuzzy systems1818 DETROJA KP, GUDI RD & PATWARDHAN SC. 2006. A possibilistic clustering approach to novel fault detection and isolation. Journal of Process Control, 16(10): 1055-1073.), (2525 EL-SHAL SM & MORRIS AS. 2000. A fuzzy expert system for fault detection in statistical process control of industrial processes. IEEE Transactions on Systems, Man and Cybernetics, Part C, 30(2): 281-289.), (5454 RAGOT J & MAQUIN D. 2006. Fault measurement detection in an urban water supply network. Journal of Process Control, 16(9): 887-902., qualitative trend analysis2020 D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC, LOSCHI RH, BACCARINI LMR & CAMINHAS WM. 2011. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian change point detection approach. Applied Soft Computing, 11(1): 179-192.), (2222 D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC & LOSCHI RH. 2011. Fuzzy/Bayesian change point detection approach to incipient fault detection. IET Control Theory & Applications, 5(4): 539-551.), (2626 FONG KF, LOH AP & TAN WW. 2008. A frequency domain approach for fault detection. International Journal of Control, 81(2): 264-276.), (4848 MAURYA MR, RENGASWAMY R & VENKATASUBRAMANIAN V. 2007. Fault diagnosis using dynamic trend analysis: A review and recent developments. Engineering Applications of Artificial Intelligence, 20(2): 133-146., mutual information6969 VERRON S, TIPLICA T & KOBI A. 2008. Fault detection and identification with a new feature selection based on mutual information. Journal of Process Control, 18(5): 479-490., neural networks1212 CALADO JMF, KORBICZ J, PATAN K, PATTON RJ & DA COSTA JMGS. 2001. Soft computing approaches to fault diagnosis for dynamic systems. European Journal of Control, 7(2-3): 248-286.), (1717 D'ANGELO MFSV & COSTA PP. 2001. Detection of shorted turns in the field winding of turbogenerators using the neural network mlp. In Proceedings of the IEEE International Conference on Systems, Man, and Cybernetics, pages 1930-1935, Tucson. (neural networks also can be used as observer), artificial immune systems4343 LAURENTYS CA, PALHARES RM & CAMINHAS WM. 2010. Design of an artificial immune system based on danger model for fault detection. Expert Systems with Applications, 37(7): 5145-5152.), (4444 LAURENTYS CA, PALHARES RM & CAMINHAS WM. 2011. A novel artificial immune system for fault behavior detection. Expert Systems with Applications, 38(11): 6957-6966.), (5757 SILVA GC, PALHARES RM & CAMINHAS WM. 2012. Immune inspired fault detection and diagnosis: A fuzzy-based approach of the negative selection algorithm and participatory clustering. Expert Systems with Applications, 39(16): 12474-12486., Bayesian networks6565 VERRON S, LI J & TIPLICA T. 2010. Fault detection and isolation of faults in a multivariate process with Bayesian network. Journal of Process Control, 20(8): 902-911.), (7070 XU BG. 2012. Intelligent fault inference for rotating flexible rotors using Bayesian belief network. Expert Systems with Applications, 39(1): 816-822. and the combination of techniques1919 D'ANGELO MFSV, PALHARES RM, COSME LB, AGUIAR LA, FONSECA FS & CAMINHAS WM. 2014. Fault detection in dynamic systems by a fuzzy/bayesian network formulation. Applied Soft Computing, 21: 647-653.), (3838 LAURENTYS CA, BOMFIM CHM, MENEZES BR & CAMINHAS WM. 2011. Design of a pipeline leakage detection using expert system: A novel approach. Applied Soft Computing, 11(1): 1057-1066..

In this paper, a new quantitative approach for fault monitoring is presented. This new approach is based on a fuzzy/Bayesian representation and was extended to the detection of multiple change points, not just one or two change points as in2020 D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC, LOSCHI RH, BACCARINI LMR & CAMINHAS WM. 2011. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian change point detection approach. Applied Soft Computing, 11(1): 179-192. and4646 MOREIRA FS, D'ANGELO MFSV, PALHARES RM & CAMINHAS WM. 2010. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian two change points detection approach. In 9th IEEE/IAS International Conference on Industry Applications., respectively. To illustrate the efficiency of the proposed methodology, the problem of monitoring the stator current of induction machine, for possible fault detection, has been solved.

Induction motors are the most important electric machinery for different industrial applications. Faults in the stator windings of three-phase induction motor represent a significant part of the failures that arise during the motor lifetime. When these motors are fed through an inverter, the situation tends to become even worse due to the voltage stresses imposed by the fast switching of the inverter1111 CRUZ SMA & CARDOSO JM. 2004. Diagnosis of stator interturn short circuits in dtc induction motor drives. IEEE Transactions on Industry Applications, 40(5): 1349-1360.. From a number of surveys, it can be realized that, for the induction motors, stator winding failures account for approximately 30% of all failures22 ARKAN M, PEROVIC DK & UNSWORTH P. 2001. Online stator fault diagnosis in induction motors. IEE Proceedings Electronic Power Application, 148(6): 537-547, November.), (3434 IRAVANI MR & KARIMI-GHARTEMANI M. 2003. Online estimation of steady state and instantaneous symmetrical components. IEE Proceedings Generation, Transmission and Distribution, 150(5): 616-622, September..

The stator winding of induction machine is subject to stress induced by a variety of factors, which include thermal overload, mechanical vibrations and voltage spikes. Deterioration of winding insulation usually begins as an inter-turn short circuit in one of the stator coils. The increased heating due to this short circuit will eventually cause turn to turn and turn to ground faults which finally lead the stator to break down6060 TALLAM RM, LEE SB, STONE G, KLIMAN GB, YOO J, HABETLER TJ & HARLEY RG. 2003. A survey of methods for detection of stator faults in induction machines. In Proceedings of SDEMPED, Diagnostics for Electric Machines, Power Eletronics and Drives, pages 35-46..

Although there is no experimental data that indicate the time delay between inter-turn and ground insulation failure, it is believed that the transition between the two states is not instantaneous. Therefore, early detection of inter turn short circuit during motor operation can be of great significance as it would eliminate subsequent damage to adjacent coils and the stator core, reducing repairing cost and motor outage time66 BOQIANG X, HEMING L & SUN LILING. 2003. Apparent impedance angle based detection of stator winding interturn short circuit fault in induction motors. In Proceedings of the Industry Application Conference, pages 1118-1125.), (5959 THOMSON WT & FENGER M. 2001. Current signature analysis to detect induction motor faults. IEEE Industry Applications Magazine, 7: 26-34, July/August..

However, early stages of deterioration are difficult to detect. In general, most of the previousreferences present approaches for dealing with abrupt faults in the stator winding, which are easier to be detected than incipient faults. In spite of these difficulties, a great deal of progress has been made on induction machine stator-winding incipient fault detection. Methods that use voltage and current measurements offer several advantages over test procedures that require machine to be taken off line or techniques that require special sensors to be mounted on the motor5858 SOTTILE J, TRUTT FC & KOHLER JL. 2000. Experimental investigation of on-line methods for incipient fault detection in induction motors. In Proceedings of the Industry Application Conference, pages 2682-2687.. Other methods, in the context of fault related to the stator-winding, can be found in99 BACCARINI LMR, MENEZES BR, GUIMARÃES HN & CAMINHAS WM. 2006. A method for early detection of stator winding faults. In Proceedings of VII International Conference on Industrial Applications, pages 1-6, Recife/Brazil.), (1010 BALLAL MS, SURYAWANSHI HM & MISHRA MK. 2008. Detection of incipient faults in induction motors using FIS, ANN and ANFIS techniques. Journal of Power Electronics, 8(2): 181-191.), (2424 DA SILVA AM, POVINELLI RJ & DEMERDASH NAO. 2008. Induction machine broken bar and stator short-circuit fault diagnostics based on three-phase stator current envelopes. IEEE Transactions on Industrial Electronics, 55(3): 1310-1318.), (5353 RODRIGUEZ PVJ & ARKKIO A. 2008. Detection of stator winding fault in induction motor using fuzzy logic. Applied Soft Computing, 8: 1112-1120.), (6161 TALLAM RM, LEE SB, STONE GC, KLIMAN GB, YOO J, HABETLER TG & HARLEY RG. 2007. A survey of methods for detection of stator-related faults in induction machines. IEEE Transactions on Industry Applications, 43(4): 920-933.. There exist other type of approaches to deal with different faults in induction machine, unlike the one considered in this paper, as, for example, dynamic eccentricity, unbalanced rotors, bearing defects and broken rotor bars (see5656 RAJAGOPALAN S, RESTREPO JA, ALLER JM, HABETLER TG & HARLEY RG. 2008. Nonstationary motor fault detection using recent quadratic time-frequency representations. IEEE Transactions on Industry Applications, 44(3): 735-744. for details and further references).

This paper is organized as follows. Section 2 shows the new fuzzy/Bayesian approach. Section 3 presents and analyzes the induction machine modeling considering the case of turn-to-turn short circuit in stator winding and shows the results for monitoring stator current of induction machine. Finally, section 4 presents the concluding remarks.

2 FUZZY/BAYESIAN APPROACH FOR MULTIPLE CHANGE POINTS DETECTION

Traditional methods of building models for time series are based on statistical techniques, aiming to select models that satisfactorily explain its dynamic behavior. However, some questions can be raised:

  • There is regime change in the series?

  • Can a single model to portray this dynamic for the whole data set?

If there are significant changes in the time series, it seems natural to find turning points before the entire modeling process. There are several papers approaching the problem of detecting points of change, as in financial series4949 OH KJ, ROH TH & MOON MS. 2005. Developing time-based clustering neural networks to use change-point detection: Application to financial time series. Asia-Pacific Journal of Operational Research, 22(1): 51-70., ecological series77 BECKAGE B, JOSEPH L, BELISLE P, WOLFSON DB & PLATT WJ. 2007. Bayesian change-point analyses in ecology. New Phytologist, 174(2): 456-467., crime rate4141 LOSCHI RH, GONÇALVES FB & CRUZ FRB. 2005. Avaliação de uma medida de evidência de um ponto de mudança e sua utilização na identificação de mudanças na taxa de criminalidade em belo horizonte. Pesquisa Operacional, 25(3): 459-463., fault detection2020 D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC, LOSCHI RH, BACCARINI LMR & CAMINHAS WM. 2011. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian change point detection approach. Applied Soft Computing, 11(1): 179-192.), (2222 D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC & LOSCHI RH. 2011. Fuzzy/Bayesian change point detection approach to incipient fault detection. IET Control Theory & Applications, 5(4): 539-551., etc. The main techniques for change points detection presented in the literature are based on statistical tests and Bayesian analysis. The most common statistical test in change points detection problem is the CUSUM and in the context of Bayesian analysis are widely used the MCMC methods. The CUSUM test proposed by3131 HINKEY DV. 1971. Inference about the change point from cumulative sum test. Biometria, 26: 279-284. is widely used in detecting change points. Some applications, modifications and extensions of this method can be seen in3232 HADJILIADIS O & MOUSTAKIDES V. 2006. Optimal and asymptotically optimal cusum rules for change point detection in the brownian motion model with multiple alternatives. Theory of Probability and its Applications, 50(1): 75-85.), (4242 LEE S, NISHIYAMA Y & YOSHIDA N. 2006. Test for parameter change in diffusion processes by cusum statistics based on one-step estimators. Annals of the Institute of Statistical Mathematics, 58(2): 211-222. and4545 LEE S, PARK S, MAEKAWA K & KAWAI K. 2006. Test for parameter change in arima models. Communications in Statistics: Simulation and Computation, 35(2): 429-439.. The context of Bayesian analysis has as a reference44 BARRY D & HARTIGAN JA. 1993. A Bayesian analysis for change point problems. Journal of the American Statistical Association, 88(421): 309-319., which uses a partition product model (MPP) for identifying points of change in the average data with a normal distribution. In3939 LOSCHI RH & CRUZ FRB. 2005. Bayesian identification of multiple change points in poisson data. Advances in Complex Systems, 8: 465-482., the posteriori probability of a time to be a turning point is used as a measure of evidence that the behavior of a data sequence has changed at some point. In4040 LOSCHI RH & CRUZ FRB. 2005. Extension to the product partition model: computing the probability of a change. Computational Statistics and Data Analysis, 48(2): 255-268., the effectiveness of the measure is evaluated in identifying changes, proposed by3939 LOSCHI RH & CRUZ FRB. 2005. Bayesian identification of multiple change points in poisson data. Advances in Complex Systems, 8: 465-482., in the rate of the Poisson distribution in sequentially observed data, and a comparison being done with the one proposed by3030 HARTIGAN JA. 1990. Partition models. Communication in Statistics-Theory and Methods, 19(8): 2745-2756.. However, all these researches require some a priori knowledge of the statistical behavior of the time series, for example, what type of distribution represents better its dynamic behavior. In this work, an approach that is independent of such a priori knowledge about the time series will be used, considering the empirical demonstration presented in2020 D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC, LOSCHI RH, BACCARINI LMR & CAMINHAS WM. 2011. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian change point detection approach. Applied Soft Computing, 11(1): 179-192. that series, after being transformed using fuzzy operations, can be adequately approximated by series with beta distribution. Thus, the parametrization of the beta distribution is used to replace a priori knowledge about the time series. Moreover, the method was extended to the detection of multiple change points, not just one or two change points as in2020 D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC, LOSCHI RH, BACCARINI LMR & CAMINHAS WM. 2011. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian change point detection approach. Applied Soft Computing, 11(1): 179-192. and4646 MOREIRA FS, D'ANGELO MFSV, PALHARES RM & CAMINHAS WM. 2010. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian two change points detection approach. In 9th IEEE/IAS International Conference on Industry Applications., respectively.

2.1 Time Series Transformation by Fuzzy Sets

The fuzzy set theory, proposed by7171 ZADEH LA. 1965. Fuzzy sets. Information and Control, 8(3): 338-353., has received much attention recently, not only in the context of theoretical developments, but also in applications. One of its main applications is in clustering methods, whereas classical methods of grouping the data into k separate categories, and in many cases some elements may not belong to a specific category, belonging to two or more categories simultaneously. Using fuzzy clustering methods is a good way to solve thisproblem because, unlike the classical approach, an element can belong to more than one category simultaneously.

Below it's proposed an alternative quantization of a time series, through fuzzy clustering, for use in the Metropolis-Hastings algorithm.

Definition (Fuzzy Clustering): Let y(t) be a time series, and consider a positive integer k. Define the set C = {Ci |min{y(t)} ≤ Ci ≤ max{y(t)}, i = 1, 2, ..., k} such that it solves the minimization problem:

(1)

so C = {Ci , i = 1, 2, ..., k} is the set of centers of the time series y(t). In (1),

(2)

is the fuzzy membership degree of y(t) for each center Ci .

Notice that, given a set C of cluster centers, it is an easy task to measure the distance of each point in the time series y(t) to each center Ci . The problem of finding the centers can be solved, for instance, via K-means3636 KAUFMAN L & ROUSSEEUW PJ. 1990. Finding groups in data: An introduction to cluster analysis. John Wiley & Sons., C-means33 BEZDEK JC. 1981. Pattern recognition with fuzzy objective function algorithms. Plenum Press., or Kohonen network3535 KOHONEN T. 2001. Self-organizing maps. Springer Series in Information Sciences. Springer., and in all these approaches, prior knowledge of the number of groups is needed, a fact that justified the use of an adaptive algorithm, where this information is not required. The optimization of the Average Silhouette Width, initially proposed by5555 ROUSSEOUW PJ. 1987. Silhouettes: a Graphical Aid to the interpretation and Validation of Cluster Analysis. Computacion and Applied Mathematics, 20: 53-65., was used in this paper.

For illustration purposes, in this paper, the time series (3) is used:

(3)

being p 1 the first operation point, p 2 the second operation point, pk the k - th operation point, ε(t) a noise with distribution π(.) in the instant t, and mi the change points.

Figure 1 illustrates the time series (3) with p 1 = 1, p 2 = 2, p 3 = 3, p 4 = 4, p 5 = 5, m 1 = 30, m 2 = 60, m 3 = 90, m 4 = 120 and ε(t)~U(0, 1) for 150 observations of the time series, Figure 2 shows the centers of time series and the membership degrees μi (t) is illustrated by Figure 3. Note that the first and last membership degrees have a single change, which will be identified using the approach proposed by2020 D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC, LOSCHI RH, BACCARINI LMR & CAMINHAS WM. 2011. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian change point detection approach. Applied Soft Computing, 11(1): 179-192., and the intermediate membership degrees have two changes to be identified by the approach proposed by4646 MOREIRA FS, D'ANGELO MFSV, PALHARES RM & CAMINHAS WM. 2010. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian two change points detection approach. In 9th IEEE/IAS International Conference on Industry Applications..

Figure 1
Time series with p1 = 1, p2 = 2, p3 = 3, p4 = 4, p5 = 5, m1 = 30, m2 = 60, m3 = 90, m4 = 120 and ε(t)~U(0, 1) for 150 observations.

Figure 2
Centers of time series of Figure 1.

The proposed fuzzy clustering to transform a given time series into a new one is described below:

  1. Input the time series y(t);

  2. Find the set of k centers, C = {Ci |min{y(t)} ≤ Ci ≤ max{y(t)}, i = 1, 2, ..., k}, that minimizes the Euclidean distance as (1), considering, for example, the time series (3) (Fig. 1).

  3. Compute the fuzzy membership degree given in (2), for each sample of the time series, y(t), with respect to each center Ci (Fig. 3).

Figure 3
Membership functions for each center of time series y(t).

2.2 Formulation of the Metropolis-Hastings algorithm

The goal of the Metropolis-Hastings algorithm2727 GAMERMAN D. 1997. Markov chain Monte Carlo: stochastic simulation for Bayesian inference. Chapman & Hall. is to construct a Markov chain that has a certain equilibrium distribution π. Define a Markov chain as follows. If Xi - 1 = xi - 1, then generate a candidate value X * from a distribution with density fX *| X (y) = 1(xi - 1, x *). The function q(.) is known as the core transition of the Markov chain. The candidate value X * is accepted or rejected with a probability of acceptance:

(4)

If the candidate is accepted, do xi = Y, otherwise do xi = Xi - 1. Thus, if the candidate is rejected, the Markov chain repeats the sequence. It is possible to show that, under generalconditions, the sequence X 0, X 1, X 2, ... is a Markov chain with equilibrium distribution π. In practical terms, the Metropolis- Hastings algorithm can be specified by the following steps:

  • 1. Choose an initial value x 0, the number of simulations, R, and make the simulations counter r = 1;

  • 2. Generate a candidate value y ~ q(xi , .);

  • 3. Calculate the probability of acceptance as in (4) and generate u ~ U(0, 1);

  • 4. Calculate the new value of the current state:

  • 5. If r<R, return to step 2. Otherwise, stop.

Note that, as discussed in2121 D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC & LOSCHI RH. 2011. A fuzzy/bayesian approach for the time series change point detection problem. Pesquisa Operacional, 31: 217-234., the quantization technique generates a time series transformed to the first and last membership degree with the following probability distribution:

μi (t) ~ Beta(a, b), for t = 1, ..., m

μi (t) ~ Beta(c, d), for t = m + 1, ..., n

The parameters to be estimated by Metropolis-Hastings algorithm are a, b, c, d and the change point m. In this type of algorithm, usually, uninformative priors are used, for example:

The Gamma distribution with parameters shape and scale equal to 0.1 was chosen because it is uninformative, with the goal of sweeping the entire parameter space. Considering the intermediate membership degree, we have the following probability distribution:

μi (t) ~ Beta(a, b), for t = 1, ..., mi - 1

μi (t) ~ Beta(c, d), for t = mi - 1 + 1, ..., mi

μi (t) ~ Beta(e, f), to mi = t + 1, ..., n

The parameters to be estimated by Metropolis- Hastings algorithm are a, b, c, d, e, f and the second point of change, mi , whereas the first switch point, mi - 1 is identified in the previous step. In this type of algorithm uninformative priors are commonly used, for example:

The transition kernels of the Markov chain for the model with a single point of change and for the model with two change points are presented in appendices A APPENDIX A This appendix is intended to build the probabilities of acceptance for the posterior distribution of the parameters a, b, c, d and m described in the Metropolis-Hastings algorithm. The reference distributions used in the work are the prior distributions, i.e., the Gamma(0.1, 0.1), that have been chosen to be uninformative. 1. For the parameter a: 2. For the parameter b: 3. For the parameter c: 4. For the parameter d: 5. For the parameter m: and B APPENDIX B This appendix is intended to show the probabilities of acceptance for the posterior distributions of the parameters a, b, c, d, e, f, m 1 and m 2 described in the Metropolis-Hastings algorithm for the model with two changes and remembering that the point m 1 will always be obtained by the previous series, it is not necessary then that will be treasured. The reference distributions used here are his own prior distributions, which in this work distributions Gamma(0.1,0.1) were chosen to be uninformative. 1. For the parameter a: 2. For the parameter b: 3. For the parameter c: 4. For the parameter d: 5. For the parameter e: 6. For the parameter f: 7. For the parameter m 2: , respectively.

2.3 Example of multiple change point detection in time series presented in Figure 1

The proposed fuzzy/Bayesian multiple change points detection in time series is described below:

  1. Input the time series y(t) (Fig. 1);

  2. Find the set of k centers, C = {Ci |min{y(t)} ≤ Ci ≤ max{y(t)}, i = 1, 2, ..., k}, given by Figure 2.

  3. Compute the fuzzy membership degree given in (2), for each sample of the time series, y(t), with respect to each center Ci , illustrated by Figure 3.

  4. Compute Metropolis-Hastings Algorithm in each membership function for change point detection as illustrated in Figure 4 (for first and last membership function is executed the detection of the one change point, and in the others membership functions is executed the detection of the two change points). The final analysis is performed as: the change point, mi , is obtained by checking where the maximum of mi histogram occurs.

Figure 4
Results of proposed methodology.

3 CASE STUDY: MONITORING INDUCTION MACHINE WITH TURN-TO-TURN SHORT CIRCUIT IN STATOR WINDING

Many studies have shown that a large proportion of induction machine faults are related to the stator-winding1010 BALLAL MS, SURYAWANSHI HM & MISHRA MK. 2008. Detection of incipient faults in induction motors using FIS, ANN and ANFIS techniques. Journal of Power Electronics, 8(2): 181-191.), (2424 DA SILVA AM, POVINELLI RJ & DEMERDASH NAO. 2008. Induction machine broken bar and stator short-circuit fault diagnostics based on three-phase stator current envelopes. IEEE Transactions on Industrial Electronics, 55(3): 1310-1318.), (5353 RODRIGUEZ PVJ & ARKKIO A. 2008. Detection of stator winding fault in induction motor using fuzzy logic. Applied Soft Computing, 8: 1112-1120.), (6161 TALLAM RM, LEE SB, STONE GC, KLIMAN GB, YOO J, HABETLER TG & HARLEY RG. 2007. A survey of methods for detection of stator-related faults in induction machines. IEEE Transactions on Industry Applications, 43(4): 920-933.. The induction machine stator-winding is subject to stress due to many factors, which include thermal overload, mechanical vibration and peak voltage caused by a speed controller. The deterioration of insulation usually begins as a short circuit fault of the stator-winding. This section describes the model that is employed here for the simulation of inter-turn short circuits in the stator windings of induction machines.

This work employs a generic model for the machine88 BACCARINI LMR, MENEZES BR & CAMINHAS WM. 2010. Fault induction dynamic model, suitable for computer simulation: simulation results and experimental validation. Mechanical Systems and Signal Processing, 24: 300-311., valid for any dq (direct and quadrature) axis speed obtained by the Park's transformation3737 KRAUSE PC. 1986. Analysis of Electric Machinery. McGraw-Hill.. Representing the currents, voltages and electromagnetic flows by i, v and λ, the resistance, leakage and mutual inductance by r, Li and Lm , the phases a, b and c by indexes a, b and c, the windings of the stator and rotor by indexes s and r, the stator and rotor voltages equations become:

(5)

(6)

where

s ] = [υas 1 υas 2 υbs υcs ]T

r ] = [υar υbr υcr ]T

[is ] = [ias ias - if ibs ics ]T

[ir ] = [iar ibr icr ]T

s ] = [λas 1 λas 2 λbs λcs ]T

r ] = [λar λbr λcr ]T

In the above, the index as 2 represents the shorted turns and if is the current in the shorted turns. Figure 5 represents the schematic diagram of a motor with an inter-turn short circuit.

Figure 5
Representation of stator windings of a motor with inter-turn short circuit.

In the model proposed in reference88 BACCARINI LMR, MENEZES BR & CAMINHAS WM. 2010. Fault induction dynamic model, suitable for computer simulation: simulation results and experimental validation. Mechanical Systems and Signal Processing, 24: 300-311., the stator windings voltages are given by:

(7)

(8)

(9)

The rotor circuit equations are the same as for traditional symmetrical model.

The stator and the rotor electromagnetic flows of stator in dq axis, are given by:

(10)

(11)

(12)

(13)

(14)

The voltage and the induced electromagnetic flow in the short-circuit turns are given by:

(15)

(16)

The electromagnetic torque is given by:

(17)

The induction machine stator current simulation results for a fault beginning which 2% of turns in the phase a become in short circuit after the time 2s and increase of 2% every time of 2s to reach a level of 8% of short-circuit is shown in the Figure 6. The root mean square (rms) current values is illustrated in Figure 7. The proposed monitoring for future fault detection approachhere is based on finding such events of non-balanced changes in the current rms values.

Figure 6
Current of phase a.

Figure 7
rms current of phase a.

3.1 Results of Proposed Methodology

The simulation results of current monitoring in the stator winding of an induction machine with star connection are shown in Figure 8, considering become in short circuit after the time 2s and increase of 2% every time of 2s to reach a level of 8% of short-circuit, was shown in the Figure 6. These results have been obtained by simulation of the induction machine using the model described in this Section and the change point detection methodology presented in Section 2. The centers in the induction machine currents are found by optimization of the Average Silhouette Width5555 ROUSSEOUW PJ. 1987. Silhouettes: a Graphical Aid to the interpretation and Validation of Cluster Analysis. Computacion and Applied Mathematics, 20: 53-65..

Figure 8
rms current of phase a and the results of proposed methodology.

4 CONCLUSIONS

In this paper a novel fuzzy/Bayesian methodology for multiple change points detection in time series has been used to treat the on-line current monitoring of the induction machine stator-winding. The methodology is based on a two-step formulation: Firstly, a fuzzy clustering generates a transformed time series with beta distribution. In the second step, a Metropolis-Hastings algorithm is used to detect the probability of the occurrence of one change point or two change points in the transformed time series. This two-step formulation allows a systematic efficient way to solve a change point detection problem, which is employed for detecting incipient faults as change points that occur in some system signal. The proposed methodology has as advantages, compared to other techniques for the FDI problem, the fact that no mathematical models of the motor and no previous knowledge of signal statistical distributions are necessary, and also an enhanced resilience against false alarms, combined with a good sensitivity that allows the detection of rather small fault signals. This methodology has been successfully applied: Simulation results have been presented as evidences of the effectiveness of the proposed methodology, even in the case of faults that cause very low level disturbances. In the future, new methods of data clustering and other uninformative priors can be used to increase the efficiency of the proposed methodology.

ACKNOWLEDGEMENTS

This work has been supported in part by the Brazilian agencies CNPq and FAPEMIG.

REFERENCES

  • 1
    ALAVI SMM, IZADI-ZAMANABADI R & HAYES MJ. 2008. Robust fault detection and isolation technique for single-input/single-output closed-loop control systems that exhibit actuator and sensor faults. IET Control Theoty and Applications, 2(11): 951-965.
  • 2
    ARKAN M, PEROVIC DK & UNSWORTH P. 2001. Online stator fault diagnosis in induction motors. IEE Proceedings Electronic Power Application, 148(6): 537-547, November.
  • 3
    BEZDEK JC. 1981. Pattern recognition with fuzzy objective function algorithms. Plenum Press.
  • 4
    BARRY D & HARTIGAN JA. 1993. A Bayesian analysis for change point problems. Journal of the American Statistical Association, 88(421): 309-319.
  • 5
    BOUKHOBZA T, HAMELIN F & CANITROT S. 2008. A graph-theoretic approach to fault detection and isolation for structured bilinear systems. International Journal of Control, 81(4): 661-678.
  • 6
    BOQIANG X, HEMING L & SUN LILING. 2003. Apparent impedance angle based detection of stator winding interturn short circuit fault in induction motors. In Proceedings of the Industry Application Conference, pages 1118-1125.
  • 7
    BECKAGE B, JOSEPH L, BELISLE P, WOLFSON DB & PLATT WJ. 2007. Bayesian change-point analyses in ecology. New Phytologist, 174(2): 456-467.
  • 8
    BACCARINI LMR, MENEZES BR & CAMINHAS WM. 2010. Fault induction dynamic model, suitable for computer simulation: simulation results and experimental validation. Mechanical Systems and Signal Processing, 24: 300-311.
  • 9
    BACCARINI LMR, MENEZES BR, GUIMARÃES HN & CAMINHAS WM. 2006. A method for early detection of stator winding faults. In Proceedings of VII International Conference on Industrial Applications, pages 1-6, Recife/Brazil.
  • 10
    BALLAL MS, SURYAWANSHI HM & MISHRA MK. 2008. Detection of incipient faults in induction motors using FIS, ANN and ANFIS techniques. Journal of Power Electronics, 8(2): 181-191.
  • 11
    CRUZ SMA & CARDOSO JM. 2004. Diagnosis of stator interturn short circuits in dtc induction motor drives. IEEE Transactions on Industry Applications, 40(5): 1349-1360.
  • 12
    CALADO JMF, KORBICZ J, PATAN K, PATTON RJ & DA COSTA JMGS. 2001. Soft computing approaches to fault diagnosis for dynamic systems. European Journal of Control, 7(2-3): 248-286.
  • 13
    CHENG H, NIKUS M & JAMSA-JOUNELA S. 2008. Fault diagnosis of the paper machine short circulation process using novel dynamic causal digraph reasoning. Journal of Process Control, 18(7-8): 676-691.
  • 14
    CHEN J & PATTON RJ. 1999. Robust model-based fault diagnosis for dynamic systems. Dordrecht: Kluwer Academic Publishers.
  • 15
    CHEN W & SAIF M. 2007. Observer-based strategies for actuator fault detection, isolation and estimation for certain class of uncertain nonlinear systems. IET Control Theoty and Applications, 1(6): 1672-1680.
  • 16
    CAMINHAS WM & TAKAHASHI RHC. 2001. Dynamic system failure detection and diagnosis employing sliding mode observers and fuzzy neural networks. In Proceedings of the Joint 9th IFSA and 20th NAFIPS, pages 304-309, Vancouver.
  • 17
    D'ANGELO MFSV & COSTA PP. 2001. Detection of shorted turns in the field winding of turbogenerators using the neural network mlp. In Proceedings of the IEEE International Conference on Systems, Man, and Cybernetics, pages 1930-1935, Tucson.
  • 18
    DETROJA KP, GUDI RD & PATWARDHAN SC. 2006. A possibilistic clustering approach to novel fault detection and isolation. Journal of Process Control, 16(10): 1055-1073.
  • 19
    D'ANGELO MFSV, PALHARES RM, COSME LB, AGUIAR LA, FONSECA FS & CAMINHAS WM. 2014. Fault detection in dynamic systems by a fuzzy/bayesian network formulation. Applied Soft Computing, 21: 647-653.
  • 20
    D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC, LOSCHI RH, BACCARINI LMR & CAMINHAS WM. 2011. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian change point detection approach. Applied Soft Computing, 11(1): 179-192.
  • 21
    D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC & LOSCHI RH. 2011. A fuzzy/bayesian approach for the time series change point detection problem. Pesquisa Operacional, 31: 217-234.
  • 22
    D'ANGELO MFSV, PALHARES RM, TAKAHASHI RHC & LOSCHI RH. 2011. Fuzzy/Bayesian change point detection approach to incipient fault detection. IET Control Theory & Applications, 5(4): 539-551.
  • 23
    DUTUIT Y & RAUZY A. 2005. Approximate estimation of system reliability via fault trees. Reliability Engineering & System Safety, 87(2): 163-172.
  • 24
    DA SILVA AM, POVINELLI RJ & DEMERDASH NAO. 2008. Induction machine broken bar and stator short-circuit fault diagnostics based on three-phase stator current envelopes. IEEE Transactions on Industrial Electronics, 55(3): 1310-1318.
  • 25
    EL-SHAL SM & MORRIS AS. 2000. A fuzzy expert system for fault detection in statistical process control of industrial processes. IEEE Transactions on Systems, Man and Cybernetics, Part C, 30(2): 281-289.
  • 26
    FONG KF, LOH AP & TAN WW. 2008. A frequency domain approach for fault detection. International Journal of Control, 81(2): 264-276.
  • 27
    GAMERMAN D. 1997. Markov chain Monte Carlo: stochastic simulation for Bayesian inference. Chapman & Hall.
  • 28
    GAO H, CHEN T & WANG L. 2008. Robust fault detection with missing measurements. International Journal of Control, 81(5): 804-819.
  • 29
    GAO Z, WANG H & CHAI T. 2007. A robust fault detection filtering for stochastic distribution systems via descriptor estimator and parametric gain design. IET Control Theoty and Applications, 1(5): 1286-1293.
  • 30
    HARTIGAN JA. 1990. Partition models. Communication in Statistics-Theory and Methods, 19(8): 2745-2756.
  • 31
    HINKEY DV. 1971. Inference about the change point from cumulative sum test. Biometria, 26: 279-284.
  • 32
    HADJILIADIS O & MOUSTAKIDES V. 2006. Optimal and asymptotically optimal cusum rules for change point detection in the brownian motion model with multiple alternatives. Theory of Probability and its Applications, 50(1): 75-85.
  • 33
    ISERMANN R & BALLE P. 1997. Trends in the application of model-based fault detection and diagnosis of technical processes. Control Engineering Practice, 5(5): 707-719.
  • 34
    IRAVANI MR & KARIMI-GHARTEMANI M. 2003. Online estimation of steady state and instantaneous symmetrical components. IEE Proceedings Generation, Transmission and Distribution, 150(5): 616-622, September.
  • 35
    KOHONEN T. 2001. Self-organizing maps. Springer Series in Information Sciences. Springer.
  • 36
    KAUFMAN L & ROUSSEEUW PJ. 1990. Finding groups in data: An introduction to cluster analysis. John Wiley & Sons.
  • 37
    KRAUSE PC. 1986. Analysis of Electric Machinery. McGraw-Hill.
  • 38
    LAURENTYS CA, BOMFIM CHM, MENEZES BR & CAMINHAS WM. 2011. Design of a pipeline leakage detection using expert system: A novel approach. Applied Soft Computing, 11(1): 1057-1066.
  • 39
    LOSCHI RH & CRUZ FRB. 2005. Bayesian identification of multiple change points in poisson data. Advances in Complex Systems, 8: 465-482.
  • 40
    LOSCHI RH & CRUZ FRB. 2005. Extension to the product partition model: computing the probability of a change. Computational Statistics and Data Analysis, 48(2): 255-268.
  • 41
    LOSCHI RH, GONÇALVES FB & CRUZ FRB. 2005. Avaliação de uma medida de evidência de um ponto de mudança e sua utilização na identificação de mudanças na taxa de criminalidade em belo horizonte. Pesquisa Operacional, 25(3): 459-463.
  • 42
    LEE S, NISHIYAMA Y & YOSHIDA N. 2006. Test for parameter change in diffusion processes by cusum statistics based on one-step estimators. Annals of the Institute of Statistical Mathematics, 58(2): 211-222.
  • 43
    LAURENTYS CA, PALHARES RM & CAMINHAS WM. 2010. Design of an artificial immune system based on danger model for fault detection. Expert Systems with Applications, 37(7): 5145-5152.
  • 44
    LAURENTYS CA, PALHARES RM & CAMINHAS WM. 2011. A novel artificial immune system for fault behavior detection. Expert Systems with Applications, 38(11): 6957-6966.
  • 45
    LEE S, PARK S, MAEKAWA K & KAWAI K. 2006. Test for parameter change in arima models. Communications in Statistics: Simulation and Computation, 35(2): 429-439.
  • 46
    MOREIRA FS, D'ANGELO MFSV, PALHARES RM & CAMINHAS WM. 2010. Incipient fault detection in induction machine stator-winding using a fuzzy-Bayesian two change points detection approach. In 9th IEEE/IAS International Conference on Industry Applications.
  • 47
    MAURYA MR, RENGASWAMY R & VENKATASUBRAMANIAN V. 2006. A signed directed graph-based systematic framework for steady-state malfunction diagnosis inside control loops. Chemical Engineering Science, 61(6): 1790-1810.
  • 48
    MAURYA MR, RENGASWAMY R & VENKATASUBRAMANIAN V. 2007. Fault diagnosis using dynamic trend analysis: A review and recent developments. Engineering Applications of Artificial Intelligence, 20(2): 133-146.
  • 49
    OH KJ, ROH TH & MOON MS. 2005. Developing time-based clustering neural networks to use change-point detection: Application to financial time series. Asia-Pacific Journal of Operational Research, 22(1): 51-70.
  • 50
    PLOIX S & ADROT O. 2006. Parity relations for linear uncertain dynamic systems. Automatica, 42(9): 1553-1562.
  • 51
    PATAN K & PARISINI T. 2005. Identification of neural dynamic models for fault detection and isolation: the case of a real sugar evaporation process. Journal of Process Control, 15(1): 67-79.
  • 52
    PUIG V, STANCU A, ESCOBET T, NEJJARI F, QUEVEDO J & PATTON RJ. 2006. Passive robust fault detection using interval observers: Application to the DAMADICS benchmark problem. Control Engineering Practice, 14(6): 621-633.
  • 53
    RODRIGUEZ PVJ & ARKKIO A. 2008. Detection of stator winding fault in induction motor using fuzzy logic. Applied Soft Computing, 8: 1112-1120.
  • 54
    RAGOT J & MAQUIN D. 2006. Fault measurement detection in an urban water supply network. Journal of Process Control, 16(9): 887-902.
  • 55
    ROUSSEOUW PJ. 1987. Silhouettes: a Graphical Aid to the interpretation and Validation of Cluster Analysis. Computacion and Applied Mathematics, 20: 53-65.
  • 56
    RAJAGOPALAN S, RESTREPO JA, ALLER JM, HABETLER TG & HARLEY RG. 2008. Nonstationary motor fault detection using recent quadratic time-frequency representations. IEEE Transactions on Industry Applications, 44(3): 735-744.
  • 57
    SILVA GC, PALHARES RM & CAMINHAS WM. 2012. Immune inspired fault detection and diagnosis: A fuzzy-based approach of the negative selection algorithm and participatory clustering. Expert Systems with Applications, 39(16): 12474-12486.
  • 58
    SOTTILE J, TRUTT FC & KOHLER JL. 2000. Experimental investigation of on-line methods for incipient fault detection in induction motors. In Proceedings of the Industry Application Conference, pages 2682-2687.
  • 59
    THOMSON WT & FENGER M. 2001. Current signature analysis to detect induction motor faults. IEEE Industry Applications Magazine, 7: 26-34, July/August.
  • 60
    TALLAM RM, LEE SB, STONE G, KLIMAN GB, YOO J, HABETLER TJ & HARLEY RG. 2003. A survey of methods for detection of stator faults in induction machines. In Proceedings of SDEMPED, Diagnostics for Electric Machines, Power Eletronics and Drives, pages 35-46.
  • 61
    TALLAM RM, LEE SB, STONE GC, KLIMAN GB, YOO J, HABETLER TG & HARLEY RG. 2007. A survey of methods for detection of stator-related faults in induction machines. IEEE Transactions on Industry Applications, 43(4): 920-933.
  • 62
    TAKAHASHI RHC & PERES PLD. 1999. Unknown input observers for uncertain systems: A unifying approach. European Journal of Control, 5(2-4): 261-275.
  • 63
    TAKAHASHI RHC, PALHARES RM & PERES PLD. 1999. Discrete-time singular observers: H2/H optimality and unknown inputs. International Journal of Control, 72(6): 481-492.
  • 64
    UPPAL FJ, PATTON RJ & WITCZAK M. 2006. A neuro-fuzzy multiple-model observer approach to robust fault diagnosis based on the DAMADICS benchmark problem. Control Engineering Practice, 14(6): 699-717.
  • 65
    VERRON S, LI J & TIPLICA T. 2010. Fault detection and isolation of faults in a multivariate process with Bayesian network. Journal of Process Control, 20(8): 902-911.
  • 66
    VENKATASUBRAMANIAN V, RENGASWAMY R & KAVURI SN. 2003. A review of process fault detection and diagnosis - part II: Qualitative models and search strategies. Computers and Chemical Engineering, 27: 313-326.
  • 67
    VENKATASUBRAMANIAN V, RENGASWAMY R, KAVURI SN & YIN K. 2003. A review of process fault detection and diagnosis - part III: Process history based methods. Computers and Chemical Engineering, 27: 327-346.
  • 68
    VENKATASUBRAMANIAN V, RENGASWAMY R, YIN K & KAVURI SN. 2003. A review of process fault detection and diagnosis - part I: Quantitative model-based methods. Computers and Chemical Engineering, 27: 293-311.
  • 69
    VERRON S, TIPLICA T & KOBI A. 2008. Fault detection and identification with a new feature selection based on mutual information. Journal of Process Control, 18(5): 479-490.
  • 70
    XU BG. 2012. Intelligent fault inference for rotating flexible rotors using Bayesian belief network. Expert Systems with Applications, 39(1): 816-822.
  • 71
    ZADEH LA. 1965. Fuzzy sets. Information and Control, 8(3): 338-353.
  • 72
    ZOLGHADRI A. 1996. An algorithm for real-time failure detection in Kalman filters. IEEE Transactions on Automatic Control, 41(10): 1537-1539.

APPENDIX A

This appendix is intended to build the probabilities of acceptance for the posterior distribution of the parameters a, b, c, d and m described in the Metropolis-Hastings algorithm. The reference distributions used in the work are the prior distributions, i.e., the Gamma(0.1, 0.1), that have been chosen to be uninformative.

  • 1. For the parameter a:

  • 2. For the parameter b:

  • 3. For the parameter c:

  • 4. For the parameter d:

  • 5. For the parameter m:

APPENDIX B

This appendix is intended to show the probabilities of acceptance for the posterior distributions of the parameters a, b, c, d, e, f, m 1 and m 2 described in the Metropolis-Hastings algorithm for the model with two changes and remembering that the point m 1 will always be obtained by the previous series, it is not necessary then that will be treasured. The reference distributions used here are his own prior distributions, which in this work distributions Gamma(0.1,0.1) were chosen to be uninformative.

  • 1. For the parameter a:

  • 2. For the parameter b:

  • 3. For the parameter c:

  • 4. For the parameter d:

  • 5. For the parameter e:

  • 6. For the parameter f:

  • 7. For the parameter m 2:

Publication Dates

  • Publication in this collection
    May-Aug 2016

History

  • Received
    27 Jan 2016
  • Accepted
    18 May 2016
Sociedade Brasileira de Pesquisa Operacional Rua Mayrink Veiga, 32 - sala 601 - Centro, 20090-050 Rio de Janeiro RJ - Brasil, Tel.: +55 21 2263-0499, Fax: +55 21 2263-0501 - Rio de Janeiro - RJ - Brazil
E-mail: sobrapo@sobrapo.org.br