SAMPLED-DATA H FILTERING FOR ROBUST KINEMATICS ESTIMATION: APPLICATIONS TO BIOMECHANICS-BASED CARDIAC IMAGE ANALYSIS
|
|
- Clifton Murphy
- 5 years ago
- Views:
Transcription
1 SAMPLED-DATA H FILTERING FOR ROBUST KINEMATICS ESTIMATION: APPLICATIONS TO BIOMECHANICS-BASED CARDIAC IMAGE ANALYSIS Shan Tong 1, Albert Sinusas 2, and Pengcheng Shi 1,3 1 Biomedical Research Laboratory, Department of Electrical and Electronic Engineering Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong 2 Section of Cardiology, Yale University School of Medicine, New Haven, CT 06520, U.S.A. 3 College of Biomedical Engineering, Southern Medical University, Guangzhou, China ABSTRACT A sampled-data H filtering strategy is proposed for the estimation of cardiac kinematic functions from periodic medical image sequences. Given the biomechanics-based myocardial dynamics, stochastic multi-frame filtering frameworks are constructed to deal with the parameter uncertainty of the biomechanical constraining model and the noisy nature of the imaging data in a coordinated fashion. As robustness is of paramount importance in cardiac wall motion estimation, especially for clinical applications, this mini-max H strategy is particulary powerful for real-world problems where the types and levels of model uncertainties and data disturbances are not available a priori. For the hybrid cardiac analysis system with continuous dynamics and discrete measurements, the state estimates are predicted according to the continuous-time state equation between observation time points, and then updated with the new measurements obtained at discrete time instants, yielding physically more meaningful and more accurate estimation results for the continuously evolving cardiac dynamics. The strategy is validated through synthetic data experiments to illustrate its advantages and on canine MR phase contrast images to show its clinical relevance. Index Terms H filtering Robust kinematics estimation, sampled-data are dealt with in such a coordinated effort, that patient-specific kinematics/material estimates can be obtained in some optimal senses [4, 5]. In [4, 6], Gaussian assumption with known statistics is made on system and data uncertainties, so that Kalman filter can be applied to obtain the minimum-mean-squared-error estimate for cardiac kinematics. In practical situations, however, the Gaussian assumption is typically unrealistic and the noise statistics is usually not available a priori. In order to relax such restrictions, H -based strategy has been proposed for the motion recovery problem [5]. However, as cardiac dynamics is a physiological process evolving continuously in time, the state estimation should be performed on the hybrid cardiac analysis system with continuous-time dynamics and discrete-time measurements. In [4, 5], the system dynamics is first converted to discrete time to obtain the predicted state estimates. Such conversion requires the assumptions of piecewise constant system input and time-invariant system, which are unrealistic in most real cases, especially for large sampling interval T. A transition matrix also needs to be calculated for the discretizing conversion, which is computationally very expensive for such systems with large degrees of freedom as in cardiac image analysis. 1. INTRODUCTION Quantitative and noninvasive estimation of cardiac kinematic functions has significant physiological and clinical implications, and accordingly there have been abundant efforts devoted to cardiac motion and deformation recovery from medical image sequences [1]. Given a set of imaging-derived, sparse, noisy measurements on cardiac kinematics, additional constraining models of mathematical or biomechanical nature are needed to constrain the ill-posed problem for a unique solution [1, 2, 3]. However, as all these constraining models are required to be known exactly as prior information, which is almost impossible in practical situations and especially for pathological cases, system uncertainties should be incorporated into the constraining models to obtain the flexibility of dealing with subject-dependent data sets. On the other hand, the imaging/imaging-derived measurements are corrupted by noises of various sources, so data uncertainties also need to be properly considered in the estimation process. Filtering strategies from control and estimation literature have been applied to couple the constraining model with the noisy observations, in which the system uncertainty of a priori constraining models and the data uncertainty of a posteriori noisy measurements In this paper, we present a sampled-data H filtering strategy for kinematics estimation of the hybrid cardiac analysis system. Instead of converting system dynamics to discrete time, state estimates are predicted according to the original continuous-time state equation and then updated with new measurements at discrete observation time points, which is physically more meaningful for the continuously evolving cardiac dynamics. As approximation errors in the discretization are avoided, estimates of higher accuracy are obtained, and such advantages of sampled-data filtering have been validated in its H 2 filtering counterpart [6]. Meanwhile, sampled-data H filtering differs from H 2-based estimation as follows. 1) No a priori knowledge of noise statistics is required. 2) The mini-max estimation criterion is to minimize the worst possible effects of the disturbances (process and measurement noises) on the state estimation errors, which will ensure that if the disturbances are small (in energy), the estimation errors will be as small as possible (in energy), regardless the noise types. These two aspects make H filtering more appropriate for such real-world problems as cardiac image analysis, where the noise statistics is unknown and robustness of the estimation is highly important. Experiments have been conducted on synthetic data to validate the advantages of this strategy, and also on canine MR images to show its clinical relevance.
2 Method 20dB 30dB 20dB 30dB (Gaussian) (Gaussian) (Poisson) (Poisson) KF 2.91± ± ± ±0.61 H 2.92± ± ± ±0.51 (a) Original configuration and material composition of the object frames #4, #8, #12, (b) Deformed geometry at selected #16 Fig. 1. Generation of the 16-frame synthetic data sequence. 2. METHODOLOGY 2.1. Biomechanics-Based Myocardial Dynamics For computational simplicity, our current implementation adopts the linear isotropic continuum material for the myocardium in order to illustrate the basic ideas of the sampled-data H filtering strategy. More complicated material models and deformation formulation can also be incorporated into this filtering framework [7]. Under spatio-temporal biomechanical constrains, the governing dynamic equation for the myocardium is derived by applying the principle of minimum potential energy [4]: MÜ + C U + KU = R (1) where M, C and K are the mass, damping and stiffness matrices respectively, R is the load vector, and U is the displacement vector. M is a known function of the material density and is temporally and spatially constant. K is a function of the material constitutive law, and is related to the material-specific Young s modulus and Poisson s ratio. C is frequency dependent, and we assume small proportional Rayleigh damping with C = αm + βk in our implementation State Space Representation for Hybrid Cardiac Analysis System In order to apply filtering strategies to our estimation problem, which are typically seen in estimation and control literature [8, 9], the above dynamic equation (1) is transformed into a state-space representation of a continuous-time system: ẋ(t) = A(t)x(t) + B(t)w(t) + ṽ(t) (2) where the state vector x(t), the input vector w(t), the system matrix A(t), and the input gain B(t) are as follows: A(t) =» x(t) = [U(t), U(t)], w(t) = [0, R(t)],, B(t) = 0 I M 1 K M 1 C» M 1, and ṽ(t) is the process noise describing the disturbances/uncertainties in cardiac dynamics. As the material model parameters and the geometry of heart will change over time, this state-space equation represents a time-varying system. Although cardiac dynamics has a continuously evolving nature, in medical image analysis and other real life systems, the system states can be observed only at discrete points in time, yielding available measurements only at discrete observation time instants. As a result, the observations provided by the imaging/imaging-derived data are expressed in a discrete-time measurement equation: y(k) = Dx(k) + e(k) (3) Table 1. Differences between the ground truth and the KF/H estimated nodal positions, (mean error ± standard deviation) Noise Type (20dB) Average Positional Error ( 10 2 ) Gaussian 2.92±0.63 Poisson 2.93±0.52 Uniform 2.91±0.55 Rayleigh 2.88±0.58 Exponential 2.89±0.62 Table 2. Comparison of average nodal positional errors from sampled-data H filter results under various types of noise. where D is the measurement matrix, e(k) is the discrete-time measurement noise accounting for the imaging data uncertainties, and k is used to denote the observation time instant t = kt. In summary, the realistic system model for cardiac image analysis is a hybrid system with continuous-time dynamics and discretetime measurements, which is represented by equations (2) and (3) Sampled-Data H Filtering for Hybrid Cardiac Analysis System In order to achieve robustness in kinematics estimation, sampleddata H filtering is applied to the hybrid system described by equations (2) and (3), which aims to estimate the states of a continuoustime system given sampled measurements of the output. The performance of the sampled-data H filter is measured by the size of the estimation error relative to the sizes of the process noise, the measurement noise, and the uncertainty in the initial state, which is defined as follows: J = x(t) ˆx(t) 2 ṽ(t) 2 + e(k) 2 + (x o ˆx o) R 0(x 0 ˆx o) Unlike the discrete H filter in [5], the performance measure of the sampled-data H filter is defined directly in terms of the continuous-time system state x(t) and disturbance ṽ(t), and intersample behavior of the system is thus taken into account. The denominator of J can be regarded as mixed L 2/l 2 norm on the uncertain disturbances affecting the system. The weighting matrix R 0 is a measure of the relative importance of the uncertainty due to the unknown initial state x o. A larger R 0 reflects smaller uncertainty in x o relative to the uncertainty in ṽ(t) and e(k). Given a prescribed noise attenuation level γ > 0, the sampleddata H filter will search ˆx(t) such that the optimal estimate of x(t) among all possible ˆx(t) should satisfy (4) sup J γ 2 (5) where the supremum is taken over all possible disturbances and initial states. The sampled-data H filter can be interpreted as a minimax problem where the estimator strategy plays against the exogenous disturbances ṽ(t), e(k) and the uncertainty in the initial state x o. The problem formulation in equation (5) guarantees the bounded estimation error over all possible disturbances of finite energy, regardless of the noise statistics. As a result, the filter achieves greater
3 Fig. 2. Comparison of displacement magnitude maps. From top to bottom: ground truth, estimation from H filter, estimation from KF. From left to right, frames #4, #8, #12 and the color scale. (Poisson noise, SNR = 20dB) robustness to disturbance variations and is well suited to such realworld problems as in cardiac analysis, where the types and levels of system disturbances and data uncertainties are not available a priori. The sampled-data H filtering algorithm for the hybrid cardiac analysis system of equations (2) and (3) is given as follows [9]: ˆx(t) = A(t)ˆx(t) + B(t)w(t) (6) ˆx(kT ) = ˆx(kT ) + P (kt )D [y(k) Dˆx(kT )] (7) where ˆx(kT ) = lim ε 0 ˆx(kT ε), and P (kt ) is the stabilizing solution to the following Riccati equation with jumps: P (t) = A(t)P (t) + P (t)a(t) + P (t)2 + I (8) γ 2 P (kt ) = P (kt )[I + D DP (kt )] 1 (9) with the initial condition P (0 ) = R 1 0. The filter given above is a linear system with finite jumps at discrete instants of time, which also has an intuitively appealing structure. In between the sampling instants, the state estimate evolves according to the continuous-time system dynamics, and the predicted state ˆx(kT ) is obtained by solving equation (6) on the time interval [(k 1)T, kt ], with the previous state estimate ˆx((k 1)T ) as the initial condition of the differential equation. Then at the observation time t = kt, the new measurement y(k) is used to update the estimate with the filter gain being P (kt )D. The filtering Riccati equation (8) is a Riccati differential equation quite similar to the continuous filtering Riccati differential equation in [10] with the measurement matrix equal to zero, as there is no measurement between observation time instants. P (kt ) is also obtained by solving the differential equation (8) with P ((k 1)T ) as the initial condition. It is necessary to point out that solving equation (8) is quite nontrivial and numerical integration is usually required. We adopt the Möbius schemes proposed in [11] for its ability of dealing with numerical instability, and the details are omitted here Discussions In [4, 5], the continuous-time state equation (2) in the hybrid system is first converted to discrete time and discrete Kalman or H filters are applied for kinematics estimation. This approach is reviewed briefly here for comparison purpose. After discretization, the state equation becomes: x(k + 1) = F (k + 1, k)x(k) + G(k)w(k) + v(k) (10) where the transition matrix F (k+1, k) = e AT, G(k) = A 1 (e AT I)B, and v(k) is the discrete-time process noise. Fig. 3. Comparison of strain maps at frame #8. From top to bottom: ground truth, estimation from H, estimation from KF and the color scales. From left to right: x-strain, y-strain, xy-strain. (Poisson noise, SNR = 20dB) It is necessary to point out that there are several implicit approximations/assumptions in the above converting process. First, the assumption that the input is piecewise constant (i.e. constant during each sampling interval) is needed for the conversion, which is not satisfied for large sampling interval T. Second, evaluation of the transition matrix F (k + 1, k) is nontrivial as it in general has no explicit form. If and only if the commutativity property and the time-invariant system assumption are satisfied, F (t + T, t) = e AT. From the above analysis, we can see that the accuracy of the estimate from discrete Kalman and H filters largely depends on the extent to which the discretized state equation (10) approximates the true continuous-time dynamics of equation (2). For large sampling interval T, the assumptions of piecewise constant input and timeinvariant system become unrealistic, and thus equation (10) is not a faithful representation of the original continuous-time state equation. In contrast, the sampled-data H filtering strategy evolves the state estimate and the filter gain according to the original continuoustime system dynamics through equations (6) and (8), and thus has several advantages over the system discretization approach. As cardiac dynamics has a continuously evolving nature, predicting in continuous time is closer to the underlying physical/physiological process, and thus physically more meaningful estimation results could be obtained. Without discretizing the state equation, the assumptions of piecewise constant input and time-invariant system are not required. In this way, the approximation errors introduced in the conversion process and propagated with the inaccurate transition matrix no longer exist, and more accurate estimates could be obtained. Such advantages of sampled-data H filtering over the discretization approach have been validated on its H 2 filtering counterpart [6], and we will focus on demonstrating the robustness of this sampled-data H strategy in next section due to space limit. And from a computational point of view, although the matrix exponential is only an approximation of the transition matrix in the system discretization approach, its calculation is nontrivial and computationally very expensive. In our strategy, such calculation is naturally avoided. Meanwhile, the sampled-data H filtering scheme can also deal with the nonuniform sampling case, offering more flexibility in the estimation process.
4 Fig. 4. Top row: Canine MRI phase contrast data (from left to right): MR intensity, x-velocity, y-velocity. Bottom row (from left to right): displacement constraints on the endo- and epi-cardial boundaries, phase contrast velocity field, and TTC-stained post mortem myocardium with infarcted tissue highlighted. 3. EXPERIMENTS AND VALIDATION 3.1. Validation with Synthetic Data To validate the accuracy and robustness of this sampled-data H filtering strategy, synthetic experiments have been performed on an object undergoing deformation in the vertical direction with the bottom being fixed. As shown in Fig.1(a), the object consists of two different materials with Young s Modulus E inside = 105, E outside = 75, and Poisson s ratio ν = sampling frames of the cyclic motion are acquired (Fig.1(b)), with displacements for all nodal points as the ground truth data set. Then only the displacements at the boundary nodes are selected and added with different types and levels of noise to generate the partial and noisy measurements. Both sampled-data H filter and its H 2 counterpart continuousdiscrete Kalman filter (KF) are implemented for the kinematics estimation. The material parameters are fixed as E = 75, ν = 0.47 in the prior constraining model. For data frames #4, #8, #12, Fig.2 shows the displacement magnitude maps of the ground truth, and the estimation results of H filter and KF from the noisy data (Poisson noise, SNR=20dB). The strain maps at frame #8 are presented in Fig.3. Point-by-point positional errors under different levels of Gaussian and Poisson noises are computed as shown in Table.1 for quantitative assessment and comparison of the H and KF results. Overall, the KF results for Poisson-corrupted data are not satisfactory, which indicates that if the Gaussian assumptions on the noise statistics are violated, small noise errors may lead to large estimation errors for KF. On the other hand, H filter produces accurate and very similar results for two sets of data contaminated by different levels of Gaussian and Poisson noises, showing its desired robustness for real-world problems. Table.2 further demonstrates the robustness of sampled-data H filter, where essentially the same results are obtained for data with various types of additive noise Canine Image Application Fig.4 demonstrates the MR phase contrast image data set. Myocardial boundaries and frame-to-frame boundary displacements are extracted using a unified active region model strategy [12], and the displacements are used as our system input data along with the midwall phase contrast velocities. The infarcted tissue is highlighted in the triphenyl tetrazolium chloride (TTC) stained post mortem myocardium (Fig.4), which provides the clinical gold standard for the assessment of the image analysis results. Fig. 5. Estimated displacement magnitude, radial, circumferential, and RC shear strain maps for frame #9 (with respect to frame #1). The estimated radial (R), circumferential (C), and R-C shear strain maps are shown in Fig.5, with the material parameters in the model set to Young s modulus E = 75000P ascal, Poisson s ratio ν = Abnormality at the infarct region can be identified from the strain maps, and the most obvious difference is observed in the RC shear strain map, where the lower-right quarter of the myocardium has much larger strain than other normal tissues. These patterns are in good agreement with the highlighted histological results of TTC stained postmortem myocardium in Fig.4, demonstrating the clinical relevance of the filtering strategy. Acknowledgment Thanks to IBM PhD fellowship for supporting Shan Tong. This work is supported in part by Hong Kong Research Grants Council CERG- HKUST6151/03E and the 973 Program of China 2003CB REFERENCES [1] A.J. Frangi, W.J. Niessen, and M.A. Viergever, Threedimensional modeling for functional analysis of cardiac images: A review, IEEE Transactions on Medical Imaging, vol. 20(1), pp. 2 25, [2] J. Park, D.N. Metaxas, and L. Axel, Analysis of left ventricular wall motion based on volumetric deformable models and mri-spamm, Medical Image Analysis, vol. 1(1), pp , [3] P. Shi, A.J. Sinusas, R.T. Constable, and J.S. Duncan, Volumetric deformation analysis using mechanics-based data fusion: Applications in cardiac motion recovery, International Journal of Computer Vision, vol. 35(1), pp , [4] P. Shi and H. Liu, Stochastic finite element framework for simultaneous estimation of cardiac kinematic functions and material parameters, Medical Image Analysis, vol. 7, pp , [5] E.W.B. Lo, H. Liu, and P. Shi, H filtering and physical modeling for robust kinematics estimation, in IEEE International Conference on Image Processing, 2003, p [6] S. Tong, A. Sinusas, and P. Shi, Continuous-discrete filtering for cardiac kinematics estimation under spatio-temporal biomechanical constrains, to appear in International Conference on Pattern Recognition, 2006.
5 [7] C.L.K. Wong and P. Shi, Finite deformation guided nonlinear filtering for multiframe cardiac motion analysis, in MICCAI, 2004, pp [8] Y. Bar-Sharlom, X.R. Li, and T. Kirubarajan, Estimation with applications to tracking and navigation, Wiley, [9] W. Sun, K. M. Nagpal, and P. Khargonekar, H control and filtering for sampled-data systems, IEEE Transactions on Automatic Control, vol. 38, pp , [10] K. M. Nagpal and P. Khargonekar, Filtering and smoothing in an H setting, IEEE Transactions on Automatic Control, vol. 36, pp , [11] J. Schiff and S. Shnider, A natural approach to the numerical integration of riccati differential equations, SIAM Journal on Numerical Analysis, vol. 36(5), pp , [12] L.N. Wong, H. Liu, A. Sinusas, and P. Shi, Spatio-temporal active region model for simultaneous segmentation and motion estimation of the whole heart, in IEEE Workshop on Variational, Geometric and Level Set Methods in Computer Vision, 2003, pp
Disturbance Attenuation Properties for Discrete-Time Uncertain Switched Linear Systems
Disturbance Attenuation Properties for Discrete-Time Uncertain Switched Linear Systems Hai Lin Department of Electrical Engineering University of Notre Dame Notre Dame, IN 46556, USA Panos J. Antsaklis
More information5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,
More informationStructural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake
Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake Seung Keun Park and Hae Sung Lee ABSTRACT This paper presents a system identification (SI) scheme
More informationOPTIMAL FUSION OF SENSOR DATA FOR DISCRETE KALMAN FILTERING Z. G. FENG, K. L. TEO, N. U. AHMED, Y. ZHAO, AND W. Y. YAN
Dynamic Systems and Applications 16 (2007) 393-406 OPTIMAL FUSION OF SENSOR DATA FOR DISCRETE KALMAN FILTERING Z. G. FENG, K. L. TEO, N. U. AHMED, Y. ZHAO, AND W. Y. YAN College of Mathematics and Computer
More informationOptimal control and estimation
Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationwith Application to Autonomous Vehicles
Nonlinear with Application to Autonomous Vehicles (Ph.D. Candidate) C. Silvestre (Supervisor) P. Oliveira (Co-supervisor) Institute for s and Robotics Instituto Superior Técnico Portugal January 2010 Presentation
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationKalman Filter and Parameter Identification. Florian Herzog
Kalman Filter and Parameter Identification Florian Herzog 2013 Continuous-time Kalman Filter In this chapter, we shall use stochastic processes with independent increments w 1 (.) and w 2 (.) at the input
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationPrediction, filtering and smoothing using LSCR: State estimation algorithms with guaranteed confidence sets
2 5th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC) Orlando, FL, USA, December 2-5, 2 Prediction, filtering and smoothing using LSCR: State estimation algorithms with
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationNonlinear Observer Design for Dynamic Positioning
Author s Name, Company Title of the Paper DYNAMIC POSITIONING CONFERENCE November 15-16, 2005 Control Systems I J.G. Snijders, J.W. van der Woude Delft University of Technology (The Netherlands) J. Westhuis
More informationMathematical Model of Blood Flow in Carotid Bifurcation
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Mathematical Model of Blood Flow in Carotid Bifurcation E. Muraca *,1, V. Gramigna 1, and G. Fragomeni 1 1 Department of Experimental Medicine
More informationDiscriminative Direction for Kernel Classifiers
Discriminative Direction for Kernel Classifiers Polina Golland Artificial Intelligence Lab Massachusetts Institute of Technology Cambridge, MA 02139 polina@ai.mit.edu Abstract In many scientific and engineering
More informationA New Approach to Tune the Vold-Kalman Estimator for Order Tracking
A New Approach to Tune the Vold-Kalman Estimator for Order Tracking Amadou Assoumane, Julien Roussel, Edgard Sekko and Cécile Capdessus Abstract In the purpose to diagnose rotating machines using vibration
More informationEvent-Triggered Decentralized Dynamic Output Feedback Control for LTI Systems
Event-Triggered Decentralized Dynamic Output Feedback Control for LTI Systems Pavankumar Tallapragada Nikhil Chopra Department of Mechanical Engineering, University of Maryland, College Park, 2742 MD,
More informationPerformance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis
Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan
More informationESTIMATOR STABILITY ANALYSIS IN SLAM. Teresa Vidal-Calleja, Juan Andrade-Cetto, Alberto Sanfeliu
ESTIMATOR STABILITY ANALYSIS IN SLAM Teresa Vidal-Calleja, Juan Andrade-Cetto, Alberto Sanfeliu Institut de Robtica i Informtica Industrial, UPC-CSIC Llorens Artigas 4-6, Barcelona, 88 Spain {tvidal, cetto,
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationNONLINEAR SAMPLED-DATA OBSERVER DESIGN VIA APPROXIMATE DISCRETE-TIME MODELS AND EMULATION
NONLINEAR SAMPLED-DAA OBSERVER DESIGN VIA APPROXIMAE DISCREE-IME MODELS AND EMULAION Murat Arcak Dragan Nešić Department of Electrical, Computer, and Systems Engineering Rensselaer Polytechnic Institute
More information1330. Comparative study of model updating methods using frequency response function data
1330. Comparative study of model updating methods using frequency response function data Dong Jiang 1, Peng Zhang 2, Qingguo Fei 3, Shaoqing Wu 4 Jiangsu Key Laboratory of Engineering Mechanics, Nanjing,
More informationAN EXTENSION OF GENERALIZED BILINEAR TRANSFORMATION FOR DIGITAL REDESIGN. Received October 2010; revised March 2011
International Journal of Innovative Computing, Information and Control ICIC International c 2012 ISSN 1349-4198 Volume 8, Number 6, June 2012 pp. 4071 4081 AN EXTENSION OF GENERALIZED BILINEAR TRANSFORMATION
More informationObservers for Bilinear State-Space Models by Interaction Matrices
Observers for Bilinear State-Space Models by Interaction Matrices Minh Q. Phan, Francesco Vicario, Richard W. Longman, and Raimondo Betti Abstract This paper formulates a bilinear observer for a bilinear
More informationENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material
ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a
More informationEvery real system has uncertainties, which include system parametric uncertainties, unmodeled dynamics
Sensitivity Analysis of Disturbance Accommodating Control with Kalman Filter Estimation Jemin George and John L. Crassidis University at Buffalo, State University of New York, Amherst, NY, 14-44 The design
More informationEDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d
Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi
More informationParameter Estimation in a Moving Horizon Perspective
Parameter Estimation in a Moving Horizon Perspective State and Parameter Estimation in Dynamical Systems Reglerteknik, ISY, Linköpings Universitet State and Parameter Estimation in Dynamical Systems OUTLINE
More informationLinear Discrete-time State Space Realization of a Modified Quadruple Tank System with State Estimation using Kalman Filter
Journal of Physics: Conference Series PAPER OPEN ACCESS Linear Discrete-time State Space Realization of a Modified Quadruple Tank System with State Estimation using Kalman Filter To cite this article:
More informationRiccati difference equations to non linear extended Kalman filter constraints
International Journal of Scientific & Engineering Research Volume 3, Issue 12, December-2012 1 Riccati difference equations to non linear extended Kalman filter constraints Abstract Elizabeth.S 1 & Jothilakshmi.R
More informationLocation Prediction of Moving Target
Location of Moving Target Department of Radiation Oncology Stanford University AAPM 2009 Outline of Topics 1 Outline of Topics 1 2 Outline of Topics 1 2 3 Estimation of Stochastic Process Stochastic Regression:
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationIncorporation of Time Delayed Measurements in a. Discrete-time Kalman Filter. Thomas Dall Larsen, Nils A. Andersen & Ole Ravn
Incorporation of Time Delayed Measurements in a Discrete-time Kalman Filter Thomas Dall Larsen, Nils A. Andersen & Ole Ravn Department of Automation, Technical University of Denmark Building 326, DK-2800
More informationSTRAIN RATE IMAGING BY DOPPLER ULTRASOUND
STRAI RATE IMAGIG BY DOPPER UTRASOUD 43.8.Qf Heimdal 1, Andreas; D hooge 2, Jan 1 Dept. of Informatics, University of Oslo; P.O. Box 18 Blindern; -316 Oslo; ORWAY; Tel:+47 22 85 27 71; Fax: +47 22 85 24
More informationHere represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.
19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that
More informationExam in Automatic Control II Reglerteknik II 5hp (1RT495)
Exam in Automatic Control II Reglerteknik II 5hp (1RT495) Date: August 4, 018 Venue: Bergsbrunnagatan 15 sal Responsible teacher: Hans Rosth. Aiding material: Calculator, mathematical handbooks, textbooks
More informationPerformance Analysis of an Adaptive Algorithm for DOA Estimation
Performance Analysis of an Adaptive Algorithm for DOA Estimation Assimakis K. Leros and Vassilios C. Moussas Abstract This paper presents an adaptive approach to the problem of estimating the direction
More informationDiscrete-Time H Gaussian Filter
Proceedings of the 17th World Congress The International Federation of Automatic Control Discrete-Time H Gaussian Filter Ali Tahmasebi and Xiang Chen Department of Electrical and Computer Engineering,
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 informationBiomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationAccelerated MRI Image Reconstruction
IMAGING DATA EVALUATION AND ANALYTICS LAB (IDEAL) CS5540: Computational Techniques for Analyzing Clinical Data Lecture 15: Accelerated MRI Image Reconstruction Ashish Raj, PhD Image Data Evaluation and
More informationStochastic structural dynamic analysis with random damping parameters
Stochastic structural dynamic analysis with random damping parameters K. Sepahvand 1, F. Saati Khosroshahi, C. A. Geweth and S. Marburg Chair of Vibroacoustics of Vehicles and Machines Department of Mechanical
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 informationSimulation of mixing of heterogeneous HE components
Chapter Simulation of mixing of heterogeneous HE components The majority on high explosives (HEs) used are blend ones. Properties of components differ that produces interaction on the grain scale (mesoprocesses).
More informationI. INTRODUCTION /01/$10.00 c 2001 IEEE IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS VOL. 37, NO. 1 JANUARY
Correspondence Comparison of Two Measurement Fusion Methods for Kalman-Filter-Based Multisensor Data Fusion Currently there exist two commonly used measurement fusion methods for Kalman-filter-based multisensor
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More informationReconstructing conductivities with boundary corrected D-bar method
Reconstructing conductivities with boundary corrected D-bar method Janne Tamminen June 24, 2011 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data,
More informationSparsity in system identification and data-driven control
1 / 40 Sparsity in system identification and data-driven control Ivan Markovsky This signal is not sparse in the "time domain" 2 / 40 But it is sparse in the "frequency domain" (it is weighted sum of six
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 informationTopology Optimization of Compliant Mechanism with Geometrical Advantage
610 Topology Optimization of Compliant Mechanism with Geometrical Advantage Seungjae MIN and Younggi KIM A compliant mechanism is a mechanism that produces its motion by the flexibility of some or all
More informationTowards a Mathematical Theory of Super-resolution
Towards a Mathematical Theory of Super-resolution Carlos Fernandez-Granda www.stanford.edu/~cfgranda/ Information Theory Forum, Information Systems Laboratory, Stanford 10/18/2013 Acknowledgements This
More informationECG Noise Filtering Using Online Model-Based Bayesian Filtering Techniques
ECG Noise Filtering Using Online Model-Based Bayesian Filtering Techniques by Aron Su A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master
More informationRobust Model Predictive Control for Autonomous Vehicle/Self-Driving Cars
Robust Model Predictive Control for Autonomous Vehicle/Self-Driving Cars Che Kun Law, Darshit Dalal, Stephen Shearrow A robust Model Predictive Control (MPC) approach for controlling front steering of
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationMyocardial Deformation from Local Frequency Estimation in Tagging MRI
Myocardial Deformation from Local Frequency Estimation in Tagging MRI L.C.M. Bruurmijn 1, H.B. Kause 2, O.G. Filatova 2, R. Duits 1,3, A. Fuster 1,3, L.M.J. Florack 1,3, and H.C. van Assen 2 1 Department
More informationANSYS Mechanical Basic Structural Nonlinearities
Lecture 4 Rate Independent Plasticity ANSYS Mechanical Basic Structural Nonlinearities 1 Chapter Overview The following will be covered in this Chapter: A. Background Elasticity/Plasticity B. Yield Criteria
More informationLarge deflection analysis of planar solids based on the Finite Particle Method
yuying@uiuc.edu 10 th US National Congress on Computational Mechanics Large deflection analysis of planar solids based on the Finite Particle Method 1, 2 Presenter: Ying Yu Advisors: Prof. Glaucio H. Paulino
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationFixed-Order Robust H Filter Design for Markovian Jump Systems With Uncertain Switching Probabilities
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 54, NO. 4, APRIL 2006 1421 Fixed-Order Robust H Filter Design for Markovian Jump Systems With Uncertain Switching Probabilities Junlin Xiong and James Lam,
More informationDiscrete-time Consensus Filters on Directed Switching Graphs
214 11th IEEE International Conference on Control & Automation (ICCA) June 18-2, 214. Taichung, Taiwan Discrete-time Consensus Filters on Directed Switching Graphs Shuai Li and Yi Guo Abstract We consider
More informationFinite element simulation of residual stresses in laser heating
IAS-2008-66-546ST Finite element simulation of residual stresses in laser heating G. H. Farrahi 1, M. Sistaninia 2, H. Moeinoddini 3 1,2-School of Mechanical Engineering, Sharif University of Technology,
More informationState observers for invariant dynamics on a Lie group
State observers for invariant dynamics on a Lie group C. Lageman, R. Mahony, J. Trumpf 1 Introduction This paper concerns the design of full state observers for state space systems where the state is evolving
More informationStress-Strain Analysis of Abdominal Aortic Wall: A Case of 3D Geometry Simulation
Energy Research Journal 1 (2): 165-170, 2010 ISSN 1949-0151 2010 Science Publications Stress-Strain Analysis of Abdominal Aortic Wall: A Case of 3D Geometry Simulation P. Khamdaengyodtai, P. Sakulchangsatjatai
More informationConditions for Suboptimal Filter Stability in SLAM
Conditions for Suboptimal Filter Stability in SLAM Teresa Vidal-Calleja, Juan Andrade-Cetto and Alberto Sanfeliu Institut de Robòtica i Informàtica Industrial, UPC-CSIC Llorens Artigas -, Barcelona, Spain
More informationA Hierarchy of Suboptimal Policies for the Multi-period, Multi-echelon, Robust Inventory Problem
A Hierarchy of Suboptimal Policies for the Multi-period, Multi-echelon, Robust Inventory Problem Dimitris J. Bertsimas Dan A. Iancu Pablo A. Parrilo Sloan School of Management and Operations Research Center,
More informationUnravel Faults on Seismic Migration Images Using Structure-Oriented, Fault-Preserving and Nonlinear Anisotropic Diffusion Filtering
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2019 SGP-TR-214 Unravel Faults on Seismic Migration Images Using Structure-Oriented,
More informationSuper-resolution via Convex Programming
Super-resolution via Convex Programming Carlos Fernandez-Granda (Joint work with Emmanuel Candès) Structure and Randomness in System Identication and Learning, IPAM 1/17/2013 1/17/2013 1 / 44 Index 1 Motivation
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationASIGNIFICANT research effort has been devoted to the. Optimal State Estimation for Stochastic Systems: An Information Theoretic Approach
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 42, NO 6, JUNE 1997 771 Optimal State Estimation for Stochastic Systems: An Information Theoretic Approach Xiangbo Feng, Kenneth A Loparo, Senior Member, IEEE,
More informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationConstitutive model of brain tissue suitable for finite element analysis of surgical procedures
Journal of Biomechanics 32 (1999 531 537 Technical Note Constitutive model of brain tissue suitable for finite element analysis of surgical procedures Karol Miller* Department of Mechanical and Materials
More informationState-Estimation Techniques for a Simple 3DOF Structure
State-Estimation Techniques for a Simple 3DOF Structure Alex Mead Zeshi Zheng I. Abstract Structural health monitoring (SHM) is a relatively new field of study in structural engineering, with the goal
More informationNonlinear Analysis of Reinforced Concrete Bridges under Earthquakes
6 th International Conference on Advances in Experimental Structural Engineering 11 th International Workshop on Advanced Smart Materials and Smart Structures Technology August 1-2, 2015, University of
More informationFIR Filters for Stationary State Space Signal Models
Proceedings of the 17th World Congress The International Federation of Automatic Control FIR Filters for Stationary State Space Signal Models Jung Hun Park Wook Hyun Kwon School of Electrical Engineering
More information4 Finite Element Method for Trusses
4 Finite Element Method for Trusses To solve the system of linear equations that arises in IPM, it is necessary to assemble the geometric matrix B a. For the sake of simplicity, the applied force vector
More informationLinear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013
Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Abstract As in optimal control theory, linear quadratic (LQ) differential games (DG) can be solved, even in high dimension,
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 information6.4 Kalman Filter Equations
6.4 Kalman Filter Equations 6.4.1 Recap: Auxiliary variables Recall the definition of the auxiliary random variables x p k) and x m k): Init: x m 0) := x0) S1: x p k) := Ak 1)x m k 1) +uk 1) +vk 1) S2:
More informationScalable robust hypothesis tests using graphical models
Scalable robust hypothesis tests using graphical models Umamahesh Srinivas ipal Group Meeting October 22, 2010 Binary hypothesis testing problem Random vector x = (x 1,...,x n ) R n generated from either
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 informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationOPTIMAL CONTROL AND ESTIMATION
OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION
More informationRobotics 2 Target Tracking. Giorgio Grisetti, Cyrill Stachniss, Kai Arras, Wolfram Burgard
Robotics 2 Target Tracking Giorgio Grisetti, Cyrill Stachniss, Kai Arras, Wolfram Burgard Linear Dynamical System (LDS) Stochastic process governed by is the state vector is the input vector is the process
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationFAST AND ACCURATE DIRECTION-OF-ARRIVAL ESTIMATION FOR A SINGLE SOURCE
Progress In Electromagnetics Research C, Vol. 6, 13 20, 2009 FAST AND ACCURATE DIRECTION-OF-ARRIVAL ESTIMATION FOR A SINGLE SOURCE Y. Wu School of Computer Science and Engineering Wuhan Institute of Technology
More informationEstimation for Nonlinear Dynamical Systems over Packet-Dropping Networks
Estimation for Nonlinear Dynamical Systems over Packet-Dropping Networks Zhipu Jin Chih-Kai Ko and Richard M Murray Abstract Two approaches, the extended Kalman filter (EKF) and moving horizon estimation
More informationMulti-Robotic Systems
CHAPTER 9 Multi-Robotic Systems The topic of multi-robotic systems is quite popular now. It is believed that such systems can have the following benefits: Improved performance ( winning by numbers ) Distributed
More informationNUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS
IGC 009, Guntur, INDIA NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS Mohammed Younus Ahmed Graduate Student, Earthquake Engineering Research Center, IIIT Hyderabad, Gachibowli, Hyderabad 3, India.
More information(q 1)t. Control theory lends itself well such unification, as the structure and behavior of discrete control
My general research area is the study of differential and difference equations. Currently I am working in an emerging field in dynamical systems. I would describe my work as a cross between the theoretical
More informationROBUST CONSTRAINED ESTIMATION VIA UNSCENTED TRANSFORMATION. Pramod Vachhani a, Shankar Narasimhan b and Raghunathan Rengaswamy a 1
ROUST CONSTRINED ESTIMTION VI UNSCENTED TRNSFORMTION Pramod Vachhani a, Shankar Narasimhan b and Raghunathan Rengaswamy a a Department of Chemical Engineering, Clarkson University, Potsdam, NY -3699, US.
More informationEVALUATING SYMMETRIC INFORMATION GAP BETWEEN DYNAMICAL SYSTEMS USING PARTICLE FILTER
EVALUATING SYMMETRIC INFORMATION GAP BETWEEN DYNAMICAL SYSTEMS USING PARTICLE FILTER Zhen Zhen 1, Jun Young Lee 2, and Abdus Saboor 3 1 Mingde College, Guizhou University, China zhenz2000@21cn.com 2 Department
More informationStructures in Seismic Zones. J. Georey Chase 2. This paper presents the ndings of a study devoted to a comparison of the eectiveness
Comparison of LQR and H1 Algorithms for Vibration Control of Structures in Seismic Zones Abstract H. Allison Smith 1 (Assoc. Member) J. Georey Chase 2 This paper presents the ndings of a study devoted
More informationH State-Feedback Controller Design for Discrete-Time Fuzzy Systems Using Fuzzy Weighting-Dependent Lyapunov Functions
IEEE TRANSACTIONS ON FUZZY SYSTEMS, VOL 11, NO 2, APRIL 2003 271 H State-Feedback Controller Design for Discrete-Time Fuzzy Systems Using Fuzzy Weighting-Dependent Lyapunov Functions Doo Jin Choi and PooGyeon
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationConstrained State Estimation Using the Unscented Kalman Filter
16th Mediterranean Conference on Control and Automation Congress Centre, Ajaccio, France June 25-27, 28 Constrained State Estimation Using the Unscented Kalman Filter Rambabu Kandepu, Lars Imsland and
More informationRecursive Generalized Eigendecomposition for Independent Component Analysis
Recursive Generalized Eigendecomposition for Independent Component Analysis Umut Ozertem 1, Deniz Erdogmus 1,, ian Lan 1 CSEE Department, OGI, Oregon Health & Science University, Portland, OR, USA. {ozertemu,deniz}@csee.ogi.edu
More informationCorrespondence. Sampled-Data Filtering with Error Covariance Assignment
666 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 49, NO. 3, MARCH 21 Correspondence Sampled-Data Filtering with Error Covariance Assignment Zidong Wang, Biao Huang, and Peijun Huo Abstract In this correspondence,
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 information