DOCSLIB.ORG

A New Nonlinear Filtering Technique for Source Localization Frederik Beutler, Uwe D

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
A New Nonlinear Filtering Technique for Source Localization Frederik Beutler, Uwe D A New Nonlinear Filtering Technique for Source Localization Frederik Beutler, Uwe D. Hanebeck Intelligent Sensor–Actuator–Systems Institute of Computer Design and Fault Tolerance Universitat¨ Karslruhe Kaiserstrasse 12, 76128 Karlsruhe, Germany [email protected], [email protected] Abstract A method for source direction estimation based on hemi- A new model-based approach for estimating the pa- sphere sampling is presented in [2], where a nonlinear trans- rameters of an arbitrary transformation between two formation is used to convert the correlation vectors of the discrete-time sequences will be introduced. One se- GCC for all TDoA measurements to a common coordinate quence is interpreted as part of a nonlinear measure- system. ment equation, the other sequence is typically measured In [3] a fast Bayesian acoustic localization method is pre- sequentially. Based on every measured value, the prob- sented, a unifying framework for localization can be found ability density function of the parameters is updated in [1]. These approaches are similar to the nonlinear filtering using a Bayesian approach. For the evolution of the technique presented here, but they use a window for source system over time, a system equation is included. The estimation from received samples. Furthermore, a motion new approach provides a high update rate for the de- model is not included. sired parameters up to the sampling rate with high ac- This paper introduces a new model-based approach for esti- curacy. It will be demonstrated for source localization mating the state, e.g. the position of a source, directly from of a speaker, where the parameters describe the position time sequences. Hence, a model for the evolution of the of the source. desired state over time can be used. The novelty is the new interpretation of a known reference sequence as part of a mea- Keywords surement equation of a nonlinear system. This new approach Source Localization, Nonlinear Filter Technique, Bayesian can be used as a framework for a multitude of procedures. Approach An important application is the use of this new approach in sensor-actuator-networks, because the received sequence INTRODUCTION does not have to be stored on the sensor node, since every Estimating the position of a target object is a typical problem received measured value will be immediately processed. arising in acoustic localization systems. A sound source is In comparison to the conventional method there are sev- placed in the environment and emits signals that are picked eral advantages. In the conventional method the estimate of up by a sensor node equipped with microphones. If the signal the time delay from the time sequences is the result from is additionally transmitted through a different medium to the processing overlapping or consecutive time windows. This sensor node, e.g. wireless radio, the position of the source disadvantage disappears in this new approach, since every can be calculated. The conventional method for this problem measured value will be used immediately. Thus the compu- consists of two steps. In the first step, the Times of Arrival tational complexity is lower and the processing delay for state (ToA) are estimated for example by means of generalized estimation is shorter compared to the conventional method, cross correlation (GCC) [6]. These estimates are then used because for every measured value the approach provides a in the second step, which converts Times of Arrival to object result for the desired state. In addition, the structure and positions based on the known microphone arrangement. This the density of measurement noise and process noise can be is in general a nonlinear optimization procedure [5]. exploited. The search complexity compared to the conven- When the signal emitted by the sound source is not known tional beamforming algorithm is reduced. In addition, prior a priori and no timing reference is given, the position of the knowledge can be exploited. source has to be calculated on the basis of Time Differences The structure of the paper is as follows. The problem of esti- of Arrival (TDoA). In [8] a Sequential Monte Carlo (SMC) mating the parameters of a transformation between two time method for source tracking based on Auxiliary Particle Filter sequences is formulated in the next section. For efficient (APF) is presented for that purpose. The method makes use estimation based on sequential measurements, a nonlinear of a model for speaker motion. Other frameworks for source filtering technique is presented. This approach is then spe- tracking based on particle filtering algorithms can be found cialized for the case of source localization. The performance in [7, 9, 10]. of the proposed new approach is evaluated in comparison to v(k) where f e(x(k)) is the estimated density from the filter step T x(k) h (x(k)) g(s(k),h(x(k)),v(k)) and f (x(k +1)) is the transition density, which depends on h(). g(). y(k) the system equation. s(k) Calculation of Point Estimate The point estimate of the state can be calculated according Figure 1. Diagram of the measurement equation. to E(x(k)) = x(k)f e(x(k))dx(k) , the conventional method. Conclusions and some hints on R future investigations are given in the last section. where E(x(k)) is the expected value of the estimated density. PROBLEM FORMULATION SPECIAL CASE: SOUND SOURCE LOCALIZATION We consider a nonlinear discrete-time system with a system We consider the far field of an acoustic point source located T equation, which describes the evolution of the system state at x =[xyz] , e.g. a speaker. The source emits a known over time as signal s(k), which is detected by several microphones. The T microphones are placed at known points M i =[xi yi zi] , x(k +1)=a(x(k),w(k)) , i =1...N in the environment. For microphone i the where a(.) is the nonlinear system equation, x(k) is the sys- received signal yi(k) according to (1) is given as tem state at time tk and w(k) is the process noise. The main yi(k)=g(s(k),hi(x(k),Mi),vi(k)) . problem of state estimation based on sequences is the recon- struction of the state x(k) from different measured values The speaker signal is delayed and attenuated depending y(k). In Fig. 1 a block diagram for the discrete-time mea- on the distance from the microphones. Assuming additive surement equation is shown. The parameters of an arbitrary measurement noise gives transformation h(.) represent the system state x(k), which yi(k)=hi(k, x(k),Mi) s(k)+vi(k) , is mapped to an intermediate state. The nonlinear function g(.) maps the sequence s(.) depending on this intermediate where stands for convolution. We assume that there is no state on the sequence y(.). The sequence y(.) is typically reflection. The time delay and attenuation of the signal can measured sequentially and can be described as be described by means of a dirac-function y(k)=g(s(k),h(x(k)),v(k)) , (1) hi(k)=Di(x(k)) · δ(k − τi(x(k))) , where v(k) is the measurement noise. where the attenuation factor and the time delay are given as 1 Di(x(k)) = NEW MODEL-BASED APPROACH FOR M i − x(k) DIRECT NONLINEAR STATE ESTIMATION and The nonlinear state estimation technique for estimating the M i − x(k) τi(x(k)) = , desired transformation parameters will be presented in this c section. where c is the acoustic velocity. The received signal at Measurement Update, Filter Step microphone i depending on the position of the microphone In the measurement update, or the filter step, the estimated M i and of the source x is given by e density f (x(k)) will be calculated according to Bayes’ law yi(k)=Di(x(k)) · s(k − τi(x(k))) + vi(k) e J p f (x(k)) = c(k) · f (x(k)) · f (x(k)) , or J where f (x(k)) is the joint density, which depends on the 1 M i − x(k) yi(k)= ·s k − +vi(k) . likelihoods corresponding to the measurement equation. M i − x(k) c f p(x(k)) is the predicted density resulting from the previous (2) estimation step and c(k) is a normalization constant. Measurement Equation Prediction Step For the derivation of the joint density, according to the mea- In the prediction step the new predicted density is given by surement update we use (2) and can thus define the likelihood f L as p T · e f (x(k + 1)) = f (x(k + 1)) f (x(k))dx(k) , L v − · − R f (x(k)) = f (yi(k) Di(x(k)) s(k τi(x(k)))) , New algorithm respectively. The speaker signal is delayed depending on 0.4 Conventional method the distance between microphone and speaker. The signals 0.2 are generated with zero-mean white Gaussian noise with the RMSE in meters 0 0 50 100 150 200 variance 1. In addition, the signals are disturbed by adding Samples zero-mean white Gaussian noise of appropriate variance of 2. The estimated probability density function for the new al- Figure 2. Comparison of the new algorithm and the gorithm is calculated on a two-dimensional grid comprising conventional method: RMSE versus number of received 50 × 50 grid points and ranging from 0 mto1 m. For every samples. received sample the position is calculated. In the conven- tional method, the Time of Arrival is estimated by means of where s(.) is piecewise constant. cross correlation, where the block length for the received se- We assume the density f v to be Gaussian, resulting in quence is set to 100.
Load more

Recommended publications
  • Linearisation of Analogue to Digital and Digital to Analogue Converters Can Be Achieved by Threshold Tracking
    LINEARISATION OF ANALOGUE TO DIGITAL AND DIGITAL TO ANALOGUE CONVERTERS ALAN CHRISTOPHER DENT, B.Sc, AMIEE, MIEEE A thesis submitted for the degree of Doctor of Philosophy, to the Faculty of Science, at the University of Edinburgh, May 1990. 7. DECLARATION OF ORIGINALITY I hereby declare that this thesis and the work reported herein was composed and originated entirely by myself in the Department of Electrical Engineering, at the University of Edinburgh, between October 1986 and May 1990. Iffl ABSTRACT Monolithic high resolution Analogue to Digital and Digital to Analogue Converters (ADC's and DAC's), cannot currently be manufactured with as much accuracy as is desirable, due to the limitations of the various fabrication technologies. The tolerance errors introduced in this way cause the converter transfer functions to be nonlinear. This nonlinearity can be quantified in static tests measuring integral nonlinearity (INL) and differential nonlinearity (DNL). In the dynamic testing of the converters, the transfer function nonlinearity is manifested as harmonic distortion, and intermodulation products. In general, regardless of the conversion technique being used, the effects of the transfer function nonlinearities get worse as converter resolution and speed increase. The result is that these nonlinearities can cause severe problems for the system designer who needs to accurately convert an input signal. A review is made of the performance of modern converters, and the existing methods of eliminating the nonlinearity, and of some of the schemes which have been proposed more recently. A new method is presented whereby code density testing techniques are exploited so that a sufficiently detailed characterisation of the converter can be made.
    [Show full text]
  • Active Non-Linear Filter Design Utilizing Parameter Plane Analysis Techniques
    1 ACTIVE NONLINEAR FILTER DESIGN UTILIZING PARAMETER PLANE ANALYSIS TECHNIQUES by Wi 1 ?am August Tate It a? .1 L4 ULUfl United States Naval Postgraduate School ACTIVE NONLINEAR FILTER DESIGN UTILIZING PARAMETER PLANE ANALYSIS TECHNIQUES by William August Tate December 1970 ThU document ka6 bzzn appnjovtd ^oh. pubtic. K.Z.- tza&c and 6 ate; Ltb dUtAtbuutlon -U untunitzd. 13781 ILU I I . Active Nonlinear Filter Design Utilizing Parameter Plane Analysis Techniques by William August Tate Lieutenant, United States Navy B.S., Leland Stanford Junior University, 1963 Submitted in partial fulfillment of the requirement for the degrees of ELECTRICAL ENGINEER and MASTER OF SCIENCE IN ELECTRICAL ENGINEERING from the NAVAL POSTGRADUATE SCHOOL December 1970 . SCHOSf NAVAL*POSTGSi&UATE MOMTEREY, CALIF. 93940 ABSTRACT A technique for designing nonlinear active networks having a specified frequency response is presented. The technique, which can be applied to monolithic/hybrid integrated-circuit devices, utilizes the parameter plane method to obtain a graphical solution for the frequency response specification in terms of a frequency-dependent resistance. A limitation to this design technique is that the nonlinear parameter must be a slowly-varying quantity - the quasi-frozen assumption. This limitation manifests itself in the inability of the nonlinear networks to filter a complex signal according to the frequency response specification. Within the constraints of this limitation, excellent correlation is attained between the specification and measurements of the frequency response of physical networks Nonlinear devices, with emphasis on voltage-controlled nonlinear FET resistances and nonlinear device frequency controllers, are realized physically and investigated. Three active linear filters, modified with the nonlinear parameter frequency are constructed, and comparisons made to the required specification.
    [Show full text]
  • Chapter 3 FILTERS
    Chapter 3 FILTERS Most images are a®ected to some extent by noise, that is unexplained variation in data: disturbances in image intensity which are either uninterpretable or not of interest. Image analysis is often simpli¯ed if this noise can be ¯ltered out. In an analogous way ¯lters are used in chemistry to free liquids from suspended impurities by passing them through a layer of sand or charcoal. Engineers working in signal processing have extended the meaning of the term ¯lter to include operations which accentuate features of interest in data. Employing this broader de¯nition, image ¯lters may be used to emphasise edges | that is, boundaries between objects or parts of objects in images. Filters provide an aid to visual interpretation of images, and can also be used as a precursor to further digital processing, such as segmentation (chapter 4). Most of the methods considered in chapter 2 operated on each pixel separately. Filters change a pixel's value taking into account the values of neighbouring pixels too. They may either be applied directly to recorded images, such as those in chapter 1, or after transformation of pixel values as discussed in chapter 2. To take a simple example, Figs 3.1(b){(d) show the results of applying three ¯lters to the cashmere ¯bres image, which has been redisplayed in Fig 3.1(a). ² Fig 3.1(b) is a display of the output from a 5 £ 5 moving average ¯lter. Each pixel has been replaced by the average of pixel values in a 5 £ 5 square, or window centred on that pixel.
    [Show full text]
  • The Unscented Kalman Filter for Nonlinear Estimation
    The Unscented Kalman Filter for Nonlinear Estimation Eric A. Wan and Rudolph van der Merwe Oregon Graduate Institute of Science & Technology 20000 NW Walker Rd, Beaverton, Oregon 97006 [email protected], [email protected] Abstract 1. Introduction The Extended Kalman Filter (EKF) has become a standard The EKF has been applied extensively to the field of non- technique used in a number of nonlinear estimation and ma- linear estimation. General application areas may be divided chine learning applications. These include estimating the into state-estimation and machine learning. We further di- state of a nonlinear dynamic system, estimating parame- vide machine learning into parameter estimation and dual ters for nonlinear system identification (e.g., learning the estimation. The framework for these areas are briefly re- weights of a neural network), and dual estimation (e.g., the viewed next. Expectation Maximization (EM) algorithm) where both states and parameters are estimated simultaneously. State-estimation This paper points out the flaws in using the EKF, and The basic framework for the EKF involves estimation of the introduces an improvement, the Unscented Kalman Filter state of a discrete-time nonlinear dynamic system, (UKF), proposed by Julier and Uhlman [5]. A central and (1) vital operation performed in the Kalman Filter is the prop- (2) agation of a Gaussian random variable (GRV) through the system dynamics. In the EKF, the state distribution is ap- where represent the unobserved state of the system and proximated by a GRV, which is then propagated analyti- is the only observed signal. The process noise drives cally through the first-order linearization of the nonlinear the dynamic system, and the observation noise is given by system.
    [Show full text]
  • Improved Audio Coding Using a Psychoacoustic Model Based on a Cochlear Filter Bank Frank Baumgarte
    IEEE TRANSACTIONS ON SPEECH AND AUDIO PROCESSING, VOL. 10, NO. 7, OCTOBER 2002 495 Improved Audio Coding Using a Psychoacoustic Model Based on a Cochlear Filter Bank Frank Baumgarte Abstract—Perceptual audio coders use an estimated masked a band and the range of spectral components that can interact threshold for the determination of the maximum permissible within a band, e.g., two sinusoids creating a beating effect. This just-inaudible noise level introduced by quantization. This es- interaction plays a crucial role in the perception of whether a timate is derived from a psychoacoustic model mimicking the properties of masking. Most psychoacoustic models for coding sound is noise-like which in turn corresponds to a significantly applications use a uniform (equal bandwidth) spectral decomposi- more efficient masking compared with a tone-like signal [2]. tion as a first step to approximate the frequency selectivity of the The noise or tone-like character is basically determined by human auditory system. However, the equal filter properties of the the amount of envelope fluctuations at the cochlear filter uniform subbands do not match the nonuniform characteristics of outputs which widely depend on the interaction of the spectral cochlear filters and reduce the precision of psychoacoustic mod- eling. Even so, uniform filter banks are applied because they are components in the pass-band of the filter. computationally efficient. This paper presents a psychoacoustic Many existing psychoacoustic models, e.g., [1], [3], and [4], model based on an efficient nonuniform cochlear filter bank employ an FFT-based transform to derive a spectral decom- and a simple masked threshold estimation.
    [Show full text]
  • Median Filter Performance Based on Different Window Sizes for Salt and Pepper Noise Removal in Gray and RGB Images
    International Journal of Signal Processing, Image Processing and Pattern Recognition Vol.8, No.10 (2015), pp.343-352 http://dx.doi.org/10.14257/ijsip.2015.8.10.34 Median Filter Performance Based on Different Window Sizes for Salt and Pepper Noise Removal in Gray and RGB Images Elmustafa S.Ali Ahmed1, Rasha E. A.Elatif2 and Zahra T.Alser3 Department of Electrical and Electronics Engineering, Red Sea University, Port Sudan, Sudan [email protected], [email protected], [email protected] Abstract Noised image is a major problem of digital image systems when images transferred between most electronic communications devices. Noise occurs in digital images due to transmit the image through the internet or maybe due to error generated by noisy sensors. Other types of errors are related to the communication system itself, since image needed to be transferred from analogue to digital and vice versa, also to be transmitted in most of the communication systems. Impulse noise added to image signal due to process of converting image signal or error from communication channel. The noise added to the original image by changes the intensity of some pixels while other remain unchanged. Salt-and-pepper noise is one of the impulse noises, to remove it a simplest way used by windowing the noisy image with a conventional median filter. Median filters are the most popular nonlinear filters extensively applied to eliminate salt- and-pepper noise. This paper evaluates the performance of median filter based on the effective median per window by using different window sizes and cascaded median filters. The performance of the proposed idea has been evaluated in MATLAB simulations on a gray and RGB images.
    [Show full text]
  • Using Nonlinear Kalman Filtering to Estimate Signals
    Using Nonlinear Kalman Filtering to Estimate Signals Dan Simon Revised September 10, 2013 It appears that no particular approximate [nonlinear] filter is consistently better than any other, though ... any nonlinear filter is better than a strictly linear one.1 The Kalman filter is a tool that can estimate the variables of a wide range of processes. In mathematical terms we would say that a Kalman filter estimates the states of a linear system. There are two reasons that you might want to know the states of a system: • First, you might need to estimate states in order to control the system. For example, the electrical engineer needs to estimate the winding currents of a motor in order to control its position. The aerospace engineer needs to estimate the velocity of a satellite in order to control its orbit. The biomedical engineer needs to estimate blood sugar levels in order to regulate insulin injection rates. • Second, you might need to estimate system states because they are interesting in their own right. For example, the electrical engineer needs to estimate power system parameters in order to predict failure probabilities. The aerospace engineer needs to estimate satellite position in order to intelligently schedule future satellite activities. The biomedical engineer needs to estimate blood protein levels in order to evaluate the health of a patient. The standard Kalman filter is an effective tool for estimation, but it is limited to linear systems. Most real-world systems are nonlinear, in which case Kalman filters do not directly apply. In the real world, nonlinear filters are used more often than linear filters, because in the real world, systems are nonlinear.
    [Show full text]
  • Nonlinear Matched Filter
    A MUTUAL INFORMATION EXTENSION TO THE MATCHED FILTER Deniz Erdogmus1, Rati Agrawal2, Jose C. Principe2 1 CSEE Department, Oregon Graduate Institute, OHSU, Portland, OR 97006, USA 2 CNEL, ECE Department, University of Florida, Gainesville, FL 32611, USA Abstract. Matched filters are the optimal linear filters demonstrated the superior performance of information for signal detection under linear channel and white noise theoretic measures over second order statistics in signal conditions. Their optimality is guaranteed in the additive processing [4]. In this general framework, we assume an white Gaussian noise channel due to the sufficiency of arbitrary instantaneous channel and additive noise with an second order statistics. In this paper, we introduce a nonlinear arbitrary distribution. Specifically in the case of linear filter for signal detection based on the Cauchy-Schwartz channels, we will consider Cauchy noise, which is a member quadratic mutual information (CS-QMI) criterion. This filter of the family of symmetric α-stable distributions, and is an is still implementing correlation but now in a high impulsive noise source. These types of noise distributions are dimensional transformed space defined by the kernel utilized known to plague signal detection [5, 6]. in estimating the CS-QMI. Simulations show that the Specifically, the proposed nonlinear signal detection nonlinear filter significantly outperforms the traditional filter is based on the Cauchy-Schwartz Quadratic Mutual matched filter in nonlinear channels, as expected. In the Information (CS-QMI) measure that has been proposed by linear channel case, the proposed filter outperforms the Principe et al. [7]. This definition of mutual information is matched filter when the received signal is corrupted by preferred because of the existence of a simple nonparametric impulsive noise such as Cauchy-distributed noise, but estimator for the CS-QMI, which in turn forms the basis of performs at the same level in Gaussian noise.
    [Show full text]
  • ECG Filtering
    12/1/13 ECG Filtering Willem Einthoven’s EKG machine, 1903 ECG Filtering n Three common noise sources q Baseline wander q Power line interference q Muscle noise n When filtering any biomedical signal care should be taken not to alter the desired information in any way n A major concern is how the QRS complex influences the output of the filter; to the filter they often pose a large unwanted impulse n Possible distortion caused by the filter should be carefully quantified 1 12/1/13 Baseline Wander Baseline Wander n Baseline wander, or extragenoeous low- frequency high-bandwidth components, can be caused by: q Perspiration (effects electrode impedance) q Respiration q Body movements n Can cause problems to analysis, especially when exmining the low-frequency ST-T segment n Two main approaches used are linear filtering and polynomial fitting 2 12/1/13 BW – Linear, time-invariant filtering n Basically make a highpass filter to cut of the lower- frequency components (the baseline wander) n The cut-off frequency should be selected so as to ECG signal information remains undistorted while as much as possible of the baseline wander is removed; hence the lowest-frequency component of the ECG should be saught. n This is generally thought to be definded by the slowest heart rate. The heart rate can drop to 40 bpm, implying the lowest frequency to be 0.67 Hz. Again as it is not percise, a sufficiently lower cutoff frequency of about 0.5 Hz should be used. n A filter with linear phase is desirable in order to avoid phase distortion that can alter various temporal realtionships in the cardiac cycle n Linear phase response can be obtained with finite impulse response, but the order needed will easily grow very high (approximately 2000) q Figure shows filters of 400 (dashdot) and 2000 (dashed) and a 5th order forward-backward filter (solid) n The complexity can be reduced by for example forward- backward IIR filtering.
    [Show full text]
  • Neural Networks for Signal Processing
    Downloaded from orbit.dtu.dk on: Oct 05, 2021 A neural architecture for nonlinear adaptive filtering of time series Hoffmann, Nils; Larsen, Jan Published in: Proceedings of the IEEE Workshop Neural Networks for Signal Processing Link to article, DOI: 10.1109/NNSP.1991.239488 Publication date: 1991 Document Version Publisher's PDF, also known as Version of record Link back to DTU Orbit Citation (APA): Hoffmann, N., & Larsen, J. (1991). A neural architecture for nonlinear adaptive filtering of time series. In Proceedings of the IEEE Workshop Neural Networks for Signal Processing (pp. 533-542). IEEE. https://doi.org/10.1109/NNSP.1991.239488 General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Users may download and print one copy of any publication from the public portal for the purpose of private study or research. You may not further distribute the material or use it for any profit-making activity or commercial gain You may freely distribute the URL identifying the publication in the public portal If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. A NEURAL ARCHITECTURE FOR NONLINEAR ADAPTIVE FILTERING OF TIME SERIES Nils Hoffmann and Jan Larsen The Computational Neural Network Center Electronics Institute, Building 349 Technical University of Denmark DK-2800 Lyngby, Denmark INTRODUCTION The need for nonlinear adaptive filtering may arise in different types of fil- tering tasks such as prediction, system identification and inverse modeling [17].
    [Show full text]
  • Adaptive Method for 1-D Signal Processing Based on Nonlinear Filter Bank and Z-Parameter
    ADAPTIVE METHOD FOR 1-D SIGNAL PROCESSING BASED ON NONLINEAR FILTER BANK AND Z-PARAMETER V. V. Lukin 1, A. A. Zelensky 1, N. O. Tulyakova 1, V. P. Melnik 2 1 State Aerospace University (Kharkov Aviation Institute), Kharkov, Ukraine E-mail: [email protected] 2 Digital Media Institute, Tampere University of Technology, Tampere, Finland E-mail: [email protected] ABSTRACT nonlinear filter and vise versa [6]. One way out is to apply an adaptive nonlinear smoother. A basic idea put behind the locally adaptive aproach is to process the noisy signal fragment by An approach to synthesis of adaptive 1-D filters based on means of nonlinear filter suited in the best (optimal) manner for nonlinear filter bank and the use of Z-parameter is put forward. situation at hand. In other words, the local signal and noise The nonlinearity of elementary filters ensures the predetermined properties and the priority of requirements of fragmentary data robustness of adaptive procedure with respect to impulsive noise filtering should be taken into consideration. By minimizing the and outliers. In turn, a local adaptation principle enables to output local errors one obtains the minimization of the total error minimize the total output error being a sum of the residual in MSE or MAE sense. fluctuation component and dynamic errors. The use of Z- There already exist many locally-adaptive (data dependent) parameter as a local activity indicator permits to “recognize ” the nonlinear 1-D filters. As examples, let us mention ones proposed signal and noise properties for given fragment quite surely and, by R.
    [Show full text]
  • Kalman Filters for Nonlinear Systems and Heavy-Tailed Noise
    Linköping studies in science and technology. Licentiate Thesis. No. 1613 Kalman Filters for Nonlinear Systems and Heavy-Tailed Noise Michael Roth LERTEKN REG IK AU L TO RO MATIC CONT LINKÖPING Division of Automatic Control Department of Electrical Engineering Linköping University, SE-581 83 Linköping, Sweden http://www.control.isy.liu.se [email protected] Linköping 2013 This is a Swedish Licentiate’s Thesis. Swedish postgraduate education leads to a Doctor’s degree and/or a Licentiate’s degree. A Doctor’s Degree comprises 240 ECTS credits (4 years of full-time studies). A Licentiate’s degree comprises 120 ECTS credits, of which at least 60 ECTS credits constitute a Licentiate’s thesis. Linköping studies in science and technology. Licentiate Thesis. No. 1613 Kalman Filters for Nonlinear Systems and Heavy-Tailed Noise Michael Roth [email protected] www.control.isy.liu.se Department of Electrical Engineering Linköping University SE-581 83 Linköping Sweden ISBN 978-91-7519-535-3 ISSN 0280-7971 LIU-TEK-LIC-2013:47 Copyright © 2013 Michael Roth Printed by LiU-Tryck, Linköping, Sweden 2013 Dedicated to all vegetarians. Abstract This thesis is on filtering in state space models. First, we examine approximate Kalman filters for nonlinear systems, where the optimal Bayesian filtering recur- sions cannot be solved exactly. These algorithms rely on the computation of cer- tain expected values. Second, the problem of filtering in linear systems that are subject to heavy-tailed process and measurement noise is addressed. Expected values of nonlinearly transformed random vectors are an essential in- gredient in any Kalman filter for nonlinear systems, because of the required joint mean vector and joint covariance of the predicted state and measurement.
    [Show full text]