System Diagnostics using Kalman Filter Estimation Error

Similar documents
A Simplified Derivation of Scalar Kalman Filter using Bayesian Probability Theory

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015

Statistical Pattern Recognition

Clustering. CM226: Machine Learning for Bioinformatics. Fall Sriram Sankararaman Acknowledgments: Fei Sha, Ameet Talwalkar.

Classification with linear models

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

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

Research on Dependable level in Network Computing System Yongxia Li 1, a, Guangxia Xu 2,b and Shuangyan Liu 3,c

Machine Learning.

Outline. Linear regression. Regularization functions. Polynomial curve fitting. Stochastic gradient descent for regression. MLE for regression

ANALYSIS OF EXPERIMENTAL ERRORS

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

Title: Damage Identification of Structures Based on Pattern Classification Using Limited Number of Sensors

Lecture 11 and 12: Basic estimation theory

First Year Quantitative Comp Exam Spring, Part I - 203A. f X (x) = 0 otherwise

R. van Zyl 1, A.J. van der Merwe 2. Quintiles International, University of the Free State

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1

Statistical Fundamentals and Control Charts

Probability and MLE.

State Space Representation

Expectation-Maximization Algorithm.

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5

Empirical Process Theory and Oracle Inequalities

Signal Processing in Mechatronics

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

Random Signals and Noise Winter Semester 2017 Problem Set 12 Wiener Filter Continuation

MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND.

HMM-Based Semantic Learning for a Mobile Robot

Tracking Bearing Degradation using Gaussian Wavelets


Research Article Health Monitoring for a Structure Using Its Nonstationary Vibration

Chapter 11 Output Analysis for a Single Model. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation

Chapter 12 EM algorithms The Expectation-Maximization (EM) algorithm is a maximum likelihood method for models that have hidden variables eg. Gaussian

A STUDY OF VIBRATION MEASURING AND FATIGUE ANALYSIS FOR CANTILEVER BEAMS

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

MECHANICAL INTEGRITY DESIGN FOR FLARE SYSTEM ON FLOW INDUCED VIBRATION

Lectures on Stochastic System Analysis and Bayesian Updating

Damped Vibration of a Non-prismatic Beam with a Rotational Spring

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Chapter 10 Advanced Topics in Random Processes

Regression and generalization

Outline. CSCI-567: Machine Learning (Spring 2019) Outline. Prof. Victor Adamchik. Mar. 26, 2019

3/8/2016. Contents in latter part PATTERN RECOGNITION AND MACHINE LEARNING. Dynamical Systems. Dynamical Systems. Linear Dynamical Systems

Change Point Estimation of Location Parameter in Multistage Processes

ME NUMERICAL METHODS Fall 2007

CDS 101: Lecture 8.2 Tools for PID & Loop Shaping

1. Introduction (Received 11 February 2013; accepted 3 June 2013)

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

APPENDIX A EARLY MODELS OF OXIDE CMP

Probability and Information Theory for Language Modeling. Statistical Linguistics. Statistical Linguistics: Adult Monolingual Speaker

The Bayesian Learning Framework. Back to Maximum Likelihood. Naïve Bayes. Simple Example: Coin Tosses. Given a generative model

PhD Preliminary Oral Exam CHARACTERIZATION AND PREDICTION OF CFD SIMULATION UNCERTAINITIES. by Serhat Hosder

UNBALANCED MACHINE FAULT DETECTION USING INSTANTANEOUS FREQUENCY

Chapter 10: Power Series

Direct Position Determination Algorithm Based on Multi-array in the Presence of Gain-phase Errors

SOME RESULTS IN MODEL BASED TRANSVERSE CRACK IDENTIFICATION IN ROTOR SYSTEMS

Exam. Notes: A single A4 sheet of paper (double sided; hand-written or computer typed)

SCORE. Exam 2. MA 114 Exam 2 Fall 2017

A Note on Effi cient Conditional Simulation of Gaussian Distributions. April 2010

Chapter 6 Sampling Distributions

MATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4

Dynamic Response of Second Order Mechanical Systems with Viscous Dissipation forces

Natural frequencies, mode shapes and modal damping values are the most important parameters to

Economics 241B Relation to Method of Moments and Maximum Likelihood OLSE as a Maximum Likelihood Estimator

CSIE/GINM, NTU 2009/11/30 1

Uncertainty. Variables. assigns to each sentence numerical degree of belief between 0 and 1. uncertainty

Sequential Monte Carlo Methods - A Review. Arnaud Doucet. Engineering Department, Cambridge University, UK

ESTIMATING UNCERTAINTIES in ANTENNA MEASUREMENTS

IE 230 Probability & Statistics in Engineering I. Closed book and notes. No calculators. 120 minutes.

Numerical Methods in Fourier Series Applications

CS284A: Representations and Algorithms in Molecular Biology

Maximum likelihood estimation from record-breaking data for the generalized Pareto distribution

1.010 Uncertainty in Engineering Fall 2008

11 Hidden Markov Models

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

Research Article Spherical Simplex-Radial Cubature Quadrature Kalman Filter

CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE

CSE 527, Additional notes on MLE & EM

Lab(8) controller design using root locus

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Pixel Recurrent Neural Networks

CHAPTER 2. Mean This is the usual arithmetic mean or average and is equal to the sum of the measurements divided by number of measurements.

TERRAIN REFERENCED UAV NAVIGATION WITH LIDAR A COMPARISON OF SEQUENTAIL PROCESSING AND BATCH PROCESSING ALGORITHMS

Error & Uncertainty. Error. More on errors. Uncertainty. Page # The error is the difference between a TRUE value, x, and a MEASURED value, x i :

THE KALMAN FILTER RAUL ROJAS

Aircraft Gas Turbine Engines' Technical Condition Identification System

GUIDELINES ON REPRESENTATIVE SAMPLING

Random Variables, Sampling and Estimation

ON-LINE ESTIMATION OF THE TRANSPIRATION IN GREENHOUSE HORTICULTURE

TABLES AND FORMULAS FOR MOORE Basic Practice of Statistics

Accuracy assessment methods and challenges

Calculus II Homework: The Comparison Tests Page 1. a n. 1 n 2 + n + 1. n= n. n=1

Problem 1. Problem Engineering Dynamics Problem Set 9--Solution. Find the equation of motion for the system shown with respect to:

Statistics 300: Elementary Statistics

6.3 Testing Series With Positive Terms

Fitting an ARIMA Process to Data

Measurement uncertainty of the sound absorption

Controlo Em Espaço de Estados. First Test

Semi-supervised Inference for Explained Variance in High-dimensional Linear Regression and Its Applications

Transcription:

System Diagostics usig Kalma Filter Estimatio Error Prof. Seugchul Lee, Seugtae Park, Heechag Kim, Hyusuk Huh ICMR2015 (11/25/2015)

Machie Health Diagostics Desig Desig DNA + + Blue Prit phagocytes lymphocytes Huma Body Egie - Data-drive - Model-based - Expert system Huma Immue System Egieerig Immue System 2

Geeral Methods for Diagostics Geeral framework of data-drive diagostics Time Tred Frequecy Tred Bode plot Orbit Aalysis 3

Model-based Diagostics Possible oly for simple systems Aalytical Computer simulatio But, expesive 4

Assumptios Model-based Diagostics system ca be approximated as state space represetatio If system dyamics are chaged (due to fault) x 1 Ax y Cx Real-time diagostics via estimatig matrix A From x ad y From A x 1 Ax y Cx x x y : state : observatio (measuremet) A: system matrix C : measuremet matrix : system oise : obervatio oise Ax 1 y Cx 5

Kalma Filter for Model-based PHM Widely used i guidace, avigatio ad cotrol of vehicles, particularly aircraft ad spacecraft, trajectory optimizatio Optimally estimate ukow states x usig a series of measuremets y observed over time, cotaiig statistical oise ad ucertaities x늿 x K y C x? 1 1 Kalma gai: T K A C ( C C V ) T t t1 t t1 1 6

Machie Learig vs. Kalma Filter Machie Learig (data-drive): Statistical approach o data Data-drive method Supervised Kalma Filter (model-based) Based o physical model + oise distributio Optimal estimator for White oise Kow covariace Kow model 7

positio Kalma Gai as Features? Questio: Moitorig the Kalma gai for system diagostics? Tested with a sprig ad mass system x(t) k 0 m x k m x x Ac x 1 0 x x Covariace matrix is statioary Kalma Gai is ot a fuctio of observatios Caot be used for system diagostics T K A C ( C C V ) T t t1 t t1 8 1

Estimatio Error: Keep moitorig e Kalma Filter Estimatio Error e y y늿 y Cx 1 Ca be used for a Health Idex (i.e., Coditio Idex) to determie the degree o coditio moitorig 9

Doe? or Ay Problems? Must kow system matrix A to apply the Kalma filter System modelig or System idetificatio problem Difficult to idetify matrix A for complicated systems Equally difficult (PHM System Idetificatio) Pseudo Matrix A p 10

Pseudo Matrix A p Positio ad velocity Real physical system positio 1 dt positio velocity 0 1 velocity k1 k k 0 km A 1 0 x 1 Ax y Cx xk 1 xk vk t k v v k1 k 1 txk 0 1 v k x k Ap v k x 1 AP x y Cx A p 1 t 0 1 11

Simulatio Validatio At time 5 sec Itetioally chaged m ad k Estimatio error x(t) k m Ap 12

System Idetificatio System + + z 1 x c v Aˆ X X ( X X ) T T 1 1 A y  K - z Maximum Likelihood z 1 xˆ 1 + c yˆ 1 ( 1) z A P xˆ 1 z 1 Kalma filter with Pseudo A p System idetificatio with MLE 13

System Dyamics: Pole Locatios System idetificatio via MLE x 1 Ax y Cx x 1 AP x y Cx x ˆ 1 Ax y Cx Compariso of pole locatios 14

Rotatig Machiery: Misaligmet Vibratio measuremet Iduce misaligmet durig operatio Kalma Filter Estimatio error Data-drive ML classificatio (SVM) Cotiuous No traiig step required 15

Bearig Fault Detectio Kalma filter with o a prior kowledge of system dyamics of matrix A Fault bearig itroduced 16

Coclusio Proposed a model-based progostics method usig Kalma filter Does ot require high fidelity models for complex systems Health idex Estimatio error ca be used for a health idex of system deviatio from kow state Future work Observability Diagosability (?) 17

Questios 18

19